The thermalization of $\gamma$-rays in radioactive expanding ejecta: A simple model and its application for Kilonovae and Ia SNe

Abstract

A semi-analytic approximation is derived for the time-dependent fraction $f_\gamma(t)$ of the energy deposited by radioactive decay $\gamma$-rays in a homologously expanding plasma of general structure. An analytic approximation is given for spherically symmetric plasma distributions. Applied to Kilonovae (KNe) associated with neutron stars mergers and Type Ia supernovae, our semi-analytic and analytic approximations reproduce, with a few percent and 10% accuracy, respectively, the energy deposition rates, $\dot{Q}_\text{dep}$, obtained in numeric Monte Carlo calculations. The time $t_\gamma$ beyond which $\gamma$-ray deposition is inefficient is determined by an effective frequency-independent $\gamma$-ray opacity $\kappa_{\gamma,\text{eff}}$, $t_\gamma = \sqrt{\kappa_{\gamma,\text{eff}}\langle\Sigma\rangle t^2}$, where $\langle\Sigma\rangle\propto t^{-2}$ is the average plasma column density. For $\beta$-decay dominated energy release, $\kappa_{\gamma,\text{eff}}$ is typically close to the effective Compton scattering opacity, $\kappa_{\gamma,\text{eff}} \approx 0.025~{\rm {cm}^{2}\,g^{-1}}$ with a weak dependence on composition. For KNe, $\kappa_{\gamma,\text{eff}}$ depends mainly on the initial electron fraction $Y_e$, $\kappa_{\gamma,\text{eff}} \approx 0.03(0.05)~{\rm {cm}^{2}\,g^{-1}}$ for $Y_e \gtrsim (\lesssim) 0.25$ (in contrast with earlier work that found $\kappa_{\gamma,\text{eff}}$ larger by 1-2 orders of magnitude for low $Y_e$), and is insensitive to the (large) nuclear physics uncertainties. Determining $t_\gamma$ from observations will therefore measure the ejecta $\langle\Sigma\rangle t^2$, providing a stringent test of models. For $\langle\Sigma\rangle t^2=2\times10^{11}~{\rm g\,{cm}^{-2}\,s^2}$, a typical value expected for KNe, $t_\gamma\approx1$ d.