2018ApJ...868...65W

Radioactive energies from unstable nuclei made in the ejecta of neutron star mergers play principal roles in powering kilonovae. In previous studies, power-law-type heating rates (e.g., ∝t^–1.3) have frequently been used, which may be inadequate if the ejecta are dominated by nuclei other than the A ∼ 130 region. We consider, therefore, two reference abundance distributions that match the r-process residuals to the solar abundances for A >= 69 (light trans-iron plus r-process elements) and A >= 90 (r-process elements). Nucleosynthetic abundances are obtained by using free-expansion models with three parameters: expansion velocity, entropy, and electron fraction. Radioactive energies are calculated as an ensemble of weighted free-expansion models that reproduce the reference abundance patterns. The results are compared with the bolometric luminosity (> a few days since merger) of the kilonova associated with GW170817. We find that the former case (fitted for A >= 69) with an ejecta mass 0.06 M_☉ reproduces the light curve remarkably well, including its steepening at >=7 days, in which the mass of r-process elements is ≃0.01 M_☉. Two β-decay chains are identified: ⁶⁶Ni →⁶⁶Cu →⁶⁶Zn and ⁷²Zn →⁷²Ga →⁷²Ge with similar halflives of parent isotopes (≃2 days), which leads to an exponential-like evolution of heating rates during 1-15 days. The light curve at late times (>40 days) is consistent with additional contributions from the spontaneous fission of ²⁵⁴Cf and a few Fm isotopes. If this is the case, the GW170817 event is best explained by the production of both light trans-iron and r-process elements that originate from dynamical ejecta and subsequent disk outflows from the neutron star merger.