Nebular spectra of kilonovae with detailed recombination rates -- I. Light r-process composition

Abstract: To investigate spectra of kilonovae in the NLTE phase (t >=1week), we perform atomic calculations for dielectronic (DR) and radiative (RR) recombination rates for the light r-process elements: Se (Z = 34), Rb (Z = 37), Sr (Z = 38), Y (Z = 39), and Zr (Z = 40) using the HULLAC code. For the different elements, our results for the total rate coefficients for recombining from the ionization states of II to I, III to II, and IV to III vary between 10^{-13}-10^{-9} cm^3/s, 10^{-12}-10^{-10} cm^3/s, and 10^{-13}-10^{-10} cm^3/s, respectively, at a temperature of T = 10,000 K. We also provide fits to the ground state photoionization cross sections of the various ions, finding larger and more slowly declining values with energy in comparison to the hydrogenic approximation. Using this new atomic data, we study the impact on kilonova model spectra at phases of t = 10 days and t = 25 days using the spectral synthesis code SUMO. Compared to models using the previous treatment of recombination as a constant rate, the new models show significant changes in ionization and temperature, and correspondingly, in emergent spectra. With the new rates, we find that Zr (Z = 40) plays a yet more dominant role in kilonova spectra for light r-process compositions. Further, we show that previously predicted mid-infrared (e.g. [Se III] 4.55 \mum) and optical (e.g. [Rb I] 7802, 7949 A) lines disappear in the new model. Instead a strong [Se I] line is seen to be emerging at \lambda=5.03 \mum. These results demonstrate the importance of considering the detailed microphysics for modelling and interpreting the late-time kilonova spectra.