Kinetic Monte Carlo simulations of excitation density dependent scintillation in Cs I and Cs I ( T l)

Nonlinear quenching of electron-hole pairs in the denser regions of ionization tracks created by [gamma]-ray and high-energy electrons is a likely cause of the light yield non-proportionality of many inorganic scintillators. Therefore, kinetic Monte Carlo (KMC) simulations were carried out to investigate the scintillation properties of pure and thallium-doped CsI as a function of electron-hole pair density. The availability of recent experimental data on the excitation density dependence of the light yield of CsI following ultraviolet excitation allowed for an improved parameterization of the interactions between self-trapped excitons (STE) in the KMC model via dipole-dipole Forster transfer. The KMC simulations reveal that nonlinear quenching occurs very rapidly (within a few picoseconds) in the early stages of the scintillation process. In addition, the simulations predict that the concentration of thallium activators can affect the extent of nonlinear quenching as it has a direct influence on the STE density through STE dissociation and electron scavenging. This improved model will enable more realistic simulations of the non-proportional [gamma]-ray and electron response of inorganic scintillators.