An analytic description of electron thermalization in kilonovae ejecta

Abstract

Abstract A simple analytic description is provided of the rate of energy deposition by β-decay electrons in the homologously expanding radioactive plasma ejected in neutron star mergers, valid for a wide range of ejecta parameters- initial entropy, electron fraction {s0, Ye} and scaled (time-independent) density ρt3. The formulae are derived using detailed numerical calculations following the time-dependent composition and β-decay emission spectra (including the effect of delayed deposition). The deposition efficiency depends mainly on ρt3 and only weakly on {s0, Ye}. The time te at which the ratio between the rates of electron energy deposition and energy production drops to 1 − e−1, is given by $t_e=t_{0e}\Big (\frac{\rho t^3}{0.5(\rho t^3)_0}\Big )^a$, where $(\rho t^3)_0=\frac{0.05\, {\rm M}_{\odot }}{4\pi (0.2c)^3}$, t0e(s0, Ye) ≈ 17 days and 0.4 ≤ a(s0, Ye) ≤ 0.5. The fractional uncertainty in te due to nuclear physics uncertainties is $\approx 10~{{\%}}$. The result a ≤ 0.5 reflects the fact that the characteristic β-decay electron energies do not decrease with time (largely due to ‘inverted decay chains’ in which a slowly-decaying isotope decays to a rapidly-decaying isotope with higher end-point energy). We provide an analytic approximation for the time-dependent electron energy deposition rate, reproducing the numerical results to better than 50% (typically $<30~{{\%}}$, well within the energy production rate uncertainty due to nuclear physics uncertainties) over a 3-4 orders-of-magnitude deposition rate decrease with time. Our results may be easily incorporated in calculations of kilonovae light curves (with general density and composition structures), eliminating the need to numerically follow the time-dependent electron spectra. Identifying te, e.g. in the bolometric light curve, will constrain the ejecta density distribution.