NLTE Spectra of Kilonovae

The electromagnetic transient following a binary neutron star merger is known as a kilonova (KN). Owing to rapid expansion velocities and small ejecta masses, KNe rapidly transition into the Non-Local Thermodynamic Equilibrium (NLTE) regime. In this study, we present synthetic NLTE spectra of KNe from 5 to 20 days after merger using the \texttt{SUMO} spectral synthesis code. We study three homogeneous composition, 1D multi-zone models with characteristic electron fractions of $Y_e \sim 0.35, 0.25$ and $0.15$. We find that emission features in the spectra tend to emerge in windows of reduced line blocking, as the ejecta are still only partially transparent even at 20 days. For the $Y_e \sim 0.35$ (lanthanide-free) ejecta, we find that the neutral and singly ionised species of Rb, Sr, Y and Zr dominate the spectra, all with good potential for identification. We directly test and confirm an impact of Sr on the 10000 angstrom spectral region in lanthanide-free ejecta, but also see that its signatures may be complex. We suggest the Rb I $\rm{5p^{1}}$- $\rm{5s^{1}}$ 7900 angstrom transition as a candidate for the $\lambda_0 \sim$ 7500--7900 angstrom P-Cygni feature in AT2017gfo. For the $Y_e \sim 0.25$ and $0.15$ compositions, lanthanides are dominant in the spectral formation, in particular Nd, Sm, and Dy. We identify key processes in KN spectral formation, notably that scattering and fluorescence play important roles even up to 20 days after merger, implying that the KN ejecta are not yet optically thin at this time.