A high-energy electron in condensed matter deposits energy by creation of electron-hole pairs whose density generally increases as the electron slows, reaching the order of 10 20 eh/cm 3 near the end of its track. The subsequent interactions of the electrons and holes include nonlinear rate terms and transport as first hot and then thermalized carriers in the nanometer-scale radial dimension of the track. Charge separation and strong radial electric fields occur in a material such as CsI with contrasting diffusion rates of self-trapped holes and hot electrons. Eventual radiative recombination has a nonlinear relation to the primary electron energy because of these interactions. This socalled intrinsic nonproportionality of electron response limits achievable energy resolution of a given scintillation radiation detector material. We use a system of coupled transport and rate equations to describe a pure host (3 equations) and one dopant (4 more equations per dopant). Applying it first to the experimentally well-characterized system of CsI and CsI:Tl in this work, we use results of picosecond absorption spectroscopy, interband z-scan measurements of nonlinear rate constants, and other experiments and calculations to determine most of the more than 20 rate and transport coefficients required for modeling. The model is solved in a track environment approximated as cylindrical and is compared to the proportionality curve and total light yield of undoped CsI at temperatures of 295 K and 100 K, as well as thallium dopant in CsI:Tl at 295 K. With this degree of validation, the space-and time-distributions of carriers and excitons, both untrapped and trapped, are examined within the model to gain an understanding of the main competitions controlling the nonproportionality of response.