Skip to main content
U.S. Department of Energy
Office of Scientific and Technical Information

Monte Carlo simulation of electron thermalization in scintillator materials: Implications for scintillator nonproportionality

Journal Article · · Journal of Applied Physics
DOI:https://doi.org/10.1063/1.4998966· OSTI ID:1426365
 [1];  [2];  [3];  [4];  [1]
  1. Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
  2. Energy and Environment Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
  3. National Security Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA
  4. Department of Nuclear Engineering and Radiological Sciences, University of Michigan, Ann Arbor, Michigan 48109, USA

The lack of reliable quantitative estimates of the length and time scales associated with hot electron thermalization after a gamma-ray induced energy cascade obscures the interplay of various microscopic processes controlling scintillator performance and hampers the search for improved detector materials. We apply a detailed microscopic kinetic Monte Carlo model of the creation and subsequent thermalization of hot electrons produced by gamma irradiation of six important scintillating crystals to determine the spatial extent of the cloud of excitations produced by gamma rays and the time required for the cloud to thermalize with the host lattice. The main ingredients of the model are ensembles of microscopic track structures produced upon gamma excitation (including the energy distribution of the excited carriers), numerical estimates of electron-phonon scattering rates, and a calculated particle dispersion to relate the speed and energy of excited carriers. All these ingredients are based on first-principles density functional theory calculations of the electronic and phonon band structures of the materials. Details of the Monte Carlo model are presented along with results for thermalization time and distance distributions. These results are discussed in light of previous work. It is found that among the studied materials, calculated thermalization distances are positively correlated with measured nonproportionality. In the important class of halide scintillators, the particle dispersion is found to be more influential than the largest phonon energy in determining the thermalization distance.

Research Organization:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Sponsoring Organization:
USDOE
DOE Contract Number:
AC05-76RL01830
OSTI ID:
1426365
Report Number(s):
PNNL-SA-128275; DN2001000
Journal Information:
Journal of Applied Physics, Vol. 122, Issue 23; ISSN 0021-8979
Publisher:
American Institute of Physics (AIP)
Country of Publication:
United States
Language:
English

References (43)

Generalized Gradient Approximation Made Simple October 1996
Selection rules for second order infrared and raman processes. I. Caesium chloride structure and interpretation of the second order Raman spectra of CsBr and Csi January 1965
The light yield nonproportionality component of scintillator energy resolution June 1998
Calculation of energy relaxation rates of fast particles by phonons in crystals March 2015
Electron-phonon scattering rates in complex polar crystals September 2017
Lattice Dynamics and Ionic Deformation in Some Alkali Halides November 1968
Monte Carlo simulation of gamma-ray response of BaF 2 and CaF 2 November 2013
Gamma-ray interaction in Ge: A Monte Carlo simulation February 2007
Picosecond dynamics of hot carriers and phonons and scintillator non-proportionality September 2012
Sample-to-Sample Variation in Single Crystal YAP:Ce Non-Proportionality February 2014
First-principles computation of material properties: the ABINIT software project November 2002
Multiscale Approach to Estimation of Scintillation Characteristics February 2014
Strontium and barium iodide high light yield scintillators February 2008
Monte Carlo simulations of electron thermalization in alkali iodide and alkaline-earth fluoride scintillators July 2012
Lattice Dynamics of Alkali Halide Crystals. II. Experimental Studies of KBr and NaI August 1963
Study of nonproportionality in the light yield of inorganic scintillators July 2011
Separable dual-space Gaussian pseudopotentials July 1996
Toward a user's toolkit for modeling scintillator non-proportionality and light yield September 2014
Light yield nonproportionality of CsI(Tl), CsI(Na), and YAP June 1998
Recombination of Correlated Electron–Hole Pairs With Account of Hot Capture With Emission of Optical Phonons October 2012
Optimized norm-conserving Vanderbilt pseudopotentials August 2013
Crystallography Open Database – an open-access collection of crystal structures May 2009
Nonproportionality of Scintillator Detectors: Theory and Experiment. II December 2011
Computer simulation of electron thermalization in CsI and CsI(Tl) September 2011
A transport-based model of material trends in nonproportionality of scintillators June 2011
Role of hot electron transport in scintillators: A theoretical study September 2016
First-principles responses of solids to atomic displacements and homogeneous electric fields: Implementation of a conjugate-gradient algorithm April 1997
Monte Carlo method for simulating γ-ray interaction with materials: A case study on Si August 2007
ABINIT: First-principles approach to material and nanosystem properties December 2009
Scintillation and Optical Properties of ${\rm BaBrI}\!:\!{\rm Eu}^{2+}$ and ${\rm CsBa}_{2}{\rm I}_{5}\!:\!{\rm Eu}^{2+}$ December 2011
Density of states and group velocity of electrons in SiO 2 calculated from a full band structure November 2002
Dynamical matrices, Born effective charges, dielectric permittivity tensors, and interatomic force constants from density-functional perturbation theory April 1997
Host structure dependence of light yield and proportionality in scintillators in terms of hot and thermalized carrier transport July 2012
Coupled rate and transport equations modeling proportionality of light yield in high-energy electron tracks: CsI at 295 K and 100 K; CsI:Tl at 295 K September 2015
Electron-Hole Pairs Created by Photons and Intrinsic Properties in Detector Materials June 2008
Excited state electronic properties of sodium iodide and cesium iodide May 2013
An analytical model of nonproportional scintillator light yield in terms of recombination rates February 2009
Radiation response of inorganic scintillators: insights from Monte Carlo simulations September 2014
Efficient pseudopotentials for plane-wave calculations January 1991
Light yield nonproportionality of CsI(Tl), CsI(Na), and YAP January 1997
The light yield nonproportionality component of scintillator energy resolution January 1997
Picosecond dynamics of hot carriers and phonons and scintillator non-proportionality January 2012
Separable Dual Space Gaussian Pseudo-potentials January 1995

Cited By (3)