3D Binary Neutron Star Merger Ejecta Evolution up to Seconds Timescale: Dynamics, Element Distribution, and Light Curves

Abstract

We present long-term, three-dimensional simulations of ejecta from four binary neutron star mergers to second-long time-scales. Numerical-relativity data serve as boundary conditions for a general-relativistic hydrodynamics evolution incorporating an equation of state valid outside nuclear statistical equilibrium and an effective nuclear-heating prescription based on reaction-network calculations. We investigate the ejecta's dynamical and geometrical properties, the impact of nuclear heating, the formation and spatial distribution of elements, and compute multi-angle kilonova light curves. % Nuclear heating significantly affects ejecta dynamics, delaying homologous expansion beyond second time-scales and reshapes the spatial distribution of heavy nuclei. This effect is largest for asymmetric binaries with long-lived remnants; extending the evolution from ${\sim}150$~ms to ${\sim}1$~s widens the angular polar region containing 90\% of the heavy-element mass (\eg, $Z=56$, $Z=79$) from $|θ|\lesssim15^{o}$ to $|θ|\lesssim30^{o}$. Our nucleosynthesis results confirm that the $^{56}$Ni$\rightarrow^{56}$Co$\rightarrow^{56}$Fe decay chain dominates the heating at $\sim$100 days, with cobalt decay producing gamma-ray lines at 846.77 and 1238.288~keV. % Comparing kilonova ray-by-ray light curves obtained from multi-angle 3D profiles and averaged 2D profiles, we find the latter approach broadly robust, although 2D light curves should be treated as upper limits. Increasing dimensionality generally lowers the bolometric luminosity, with binary asymmetry strengthening the viewing-angle dependence. For observers aligned with a lanthanide curtain's densest region, 3D emission can match its 2D counterpart in brightness. We conclude that increased dimensionality alone is unlikely to reconcile current theoretical models with AT2017gfo observations.