Transport and mixing of r-process elements in neutron star binary merger blast waves

TL;DR

This study uses 3D hydrodynamic simulations to analyze how neutron star binary merger ejecta evolve and impact the interstellar medium, revealing similarities to supernova remnants and implications for galaxy feedback models.

Contribution

It provides the first detailed hydrodynamic modeling of NSBM ejecta evolution in galactic environments, connecting merger physics with galaxy-scale feedback processes.

Findings

NSBM remnants have complex early morphologies but similar late-time evolution to supernova remnants.

Shell formation occurs at about 10^3 solar masses of swept-up material, with significant r-process enrichment.

Results support using supernova feedback models for NSBM ejecta in galaxy simulations.

Abstract

The r-process nuclei are robustly synthesized in the material ejected during a neutron star binary merger (NSBM), as tidal torques transport angular momentum and energy through the outer Lagrange point in the form of a vast tidal tail. If NSBM are indeed solely responsible for the solar system r- process abundances, a galaxy like our own would require to host a few NSBM per million years, with each event ejecting, on average, about 5x10^{-2} M_sun of r-process material. Because the ejecta velocities in the tidal tail are significantly larger than in ordinary supernovae, NSBM deposit a comparable amount of energy into the interstellar medium (ISM). In contrast to extensive efforts studying spherical models for supernova remnant evolution, calculations quantifying the impact of NSBM ejecta in the ISM have been lacking. To better understand their evolution in a cosmological context, we perform a suite of three-dimensional hydrodynamic simulations with optically-thin radiative cooling of isolated NSBM ejecta expanding in environments with conditions adopted from Milky Way-like galaxy simulations. Although the remnant morphology is highly complex at early times, the subsequent radiative evolution that results from thermal instability in atomic gas is remarkably similar to that of a standard supernova blast wave. This implies that sub-resolution supernova feedback models can be used in galaxy-scale simulations that are unable to resolve the key evolutionary phases of NSBM blast waves. Among other quantities, we examine the radius, time, mass and kinetic energy content of the NSBM remnant at shell formation as well as the momentum injected to the ISM. We find that the shell formation epoch is attained when the swept-up mass is about 10^3 M_sun, at this point the mass fraction of r-process material is drastically enhanced up to two orders of magnitude in relation to a solar metallicity ISM.