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

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-princi...