Non-Thermal Ionization of Kilonova Ejecta: Observable Impacts

TL;DR

This paper introduces an approximate method to account for non-local thermodynamic equilibrium effects caused by high-energy electrons in kilonova ejecta, improving the accuracy of light curve modeling and spectral predictions.

Contribution

It presents a quasi-NLTE ionization calculation method that captures the impact of high-energy electrons, enhancing kilonova models and resolving previous modeling tensions.

Findings

Higher ionization levels reduce optical line blanketing.

Models require less ejecta mass for observed brightness.

Explains observed spectral features without fine-tuning.

Abstract

The characteristic rapid rise and decline at optical wavelengths of a kilonova is the product of the low ejecta mass ($\lesssim 0.05 M_\odot$) and high ejecta velocity ($\gtrsim 0.1$c). We show that, even at very early times ($\lesssim 2$ days), regions of ejecta fall below critical density and temperature thresholds at which non-local thermodynamic equilibrium (NLTE) effects become important. Here, we present an approximate method for calculating the ionization state of the ejecta that accounts for the NLTE impact of high-energy electrons produced in the beta decay of freshly synthesized $r$-process elements. We find that incorporating ionization from high-energy electrons produces an ``inverted" and ``blended" ionization structure, where the most highly ionized species are located in the fastest moving homologous ejecta and multiple ionization states coexist. In radiation transport calculations, the higher degree of ionization reduces line blanketing in optical bands, leading to improved agreement with the light curve properties of AT\,2017gfo such as the duration, decay rates, brightness, and colors. Our quasi-NLTE implementation helps to alleviate tensions in kilonova modeling: for high-velocity ($\sim 0.3c$) ejecta components our models require less mass for a given peak brightness in optical bands, by as much as a factor of 3; our models can explain the presence of observed features associated to Sr II, W III, Se III, and Te III under conditions where LTE models would predict only neutral species; and we naturally predict the coexistence of species like Sr II and Ce III without the need for fine-tuning of the ejecta properties.