Radioactivity and thermalization in the ejecta of compact object mergers and their impact on kilonova light curves

One promising electromagnetic signature of compact object mergers are kilonovae: approximately isotropic radioactively powered transients that peak days to weeks post-merger. Key uncertainties in kilonova modeling include the emission profiles of the radioactive decay products—non-thermal $\beta$-particles, $\alpha$-particles, fission fragments, and $\gamma$-rays—and the efficiency with which their kinetic energy is absorbed by the ejecta. The radioactive energy emitted, along with its thermalization efficiency, sets the luminosity budget and is therefore crucial for predicting kilonova light curves. We outline uncertainties in the radioactivity, describe the processes by which the decay products transfer energy to the ejecta, and calculate time-dependent thermalization efficiencies for each particle type. We determine the net thermalization efficiency and explore its dependence on r-process yields—in particular, the production of $\alpha$-decaying translead nuclei—and on ejecta mass, velocity, and magnetic fields. We incorporate our results into detailed radiation transport simulations, and calculate updated kilonova light curve predictions. Thermalization effects reduce kilonova luminosities by a factor of roughly 2 at peak, and by an order of magnitude at later times (15 days or more after explosion). We present analytic fits to time-dependent thermalization efficiencies, which can be used to improve light curve models. We revisit the putative kilonova that accompanied gamma-ray burst 130603B, and estimate the mass ejected in that event. We find later time kilonova light curves can be significantly impacted by $\alpha$-decay from translead isotopes; data at these times may therefore be diagnostic of ejecta abundances.