Late-time emission-line profiles from kilonova models

Abstract Numerical simulations suggest that neutron star mergers eject material with complex, non-spherical density and composition distributions. Here we use two-dimensional configurations of merger ejecta obtained from long-term hydrodynamic simulations to quantify the influence of such ejecta structure on the shapes of spectral lines in the optically thin limit. We consider three example elements of interest for kilonova modelling (selenium, tellurium and tungsten) and illustrate profile shapes for a sample of models and observer orientations. Many of our calculations yield complex profile shapes, including cases with multiple peaks and/or extended wings on scales large enough to be relevant to interpreting observations. For selenium and tellurium, our late-phase profile shapes are most sensitive to the structure of the low-velocity ejecta (≲ 0.1 c) launched after the merger from the relic black-hole torus system, while for heavier elements the contribution from the more rapidly expanding and more neutron-rich dynamical ejecta launched right after the merger is more significant and leads to broader line shapes. We also find that the dynamical influence of heating due to the decay of r-process elements can lead to considerably broader peaks than suggested by models that neglect this effect. Although idealised, our calculations demonstrate that line shapes are sensitive to the ejecta structure and could therefore constrain the polar observation angle or underlying properties of the merger that determine the spatial distributions of elements in the ejecta components, such as the binary mass ratio or even the equation of state of high-density matter.