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

A semi-analytic approximation is derived for the time-dependent fraction f_γ (t) of the energy deposited by radioactive decay γ-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 per cent and 10 per cent accuracy, respectively, the energy deposition rates, Q̇dep, obtained in numeric Monte Carlo calculations. The time t_γ beyond which γ-ray deposition is inefficient is determined by an effective frequency-independent γ-ray opacity κ γ ,eff, t_γ = √κ γ ,eff Σ t², where Σ ∝ t⁻² is the average plasma column density. For β-decay dominated energy release, κ γ ,eff is typically close to the effective Compton scattering opacity, κ γ ,eff ≈ 0.025 cm² g⁻¹ with a weak dependence on composition. For KNe, κ γ ,eff depends mainly on the initial electron fraction Yₑ, κ γ ,eff ≈ 0.03(0.05) cm² g⁻¹ for Yₑ () 0.25 (in contrast with earlier work that found κ γ ,eff larger by 1–2 orders of magnitude for low Yₑ), and is insensitive to the (large) nuclear physics uncertainties. Determining t_γ from observations will therefore measure the ejecta Σ t², providing a stringent test of models. For Σ t²=2× 10¹¹~ g\, cm⁻²\, s², a typical value expected for KNe, t_γ ≈ 1 d.