We present a systematic numerical relativity study of the dynamical ejecta, winds, and nucleosynthesis in neutron star (NS) merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOSs) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamic turbulence up to 100 ms postmerger. Spiral density waves propagating from the NS remnant to the disk trigger a wind with mass flux ∼0.1–0.5 M ⊙ s−1, which persists for the entire simulation as long as the remnant does not collapse to a black hole. This wind has average electron fraction ≳0.3 and average velocity ∼0.1–0.17 c and thus is a site for the production of weak r-process elements (mass number A < 195). Disks around long-lived remnants have masses ∼0.1–0.2 M ⊙, temperatures peaking at ≲10 MeV near the inner edge, and a characteristic double-peak distribution in entropy resulting from shocks propagating through the disk. The dynamical and spiral-wave ejecta computed in our targeted simulations are not compatible with those inferred from AT2017gfo using two-components kilonova models. Rather, they indicate that multicomponent kilonova models including disk winds are necessary to interpret AT2017gfo. The nucleosynthesis in the combined dynamical ejecta and spiral-wave wind in the long-lived mergers of comparable mass robustly accounts for all the r-process peaks, from mass number ∼75 to actinides in terms of solar abundances. Total abundances are weakly dependent on the EOS, while the mass ratio affects the production of first-peak elements.

We present new numerical relativity results of neutron star mergers with chirp mass 1.188M⊙ and mass ratios q = 1.67 and q = 1.8 using finite-temperature equations of state (EOS), approximate neutrino transport and a subgrid model for magnetohydrodynamics-induced turbulent viscosity. The EOS are compatible with nuclear and astrophysical constraints and include a new microphysical model derived from ab-initio calculations based on the Brueckner-Hartree-Fock approach. We report for the first time evidence for accretion-induced prompt collapse in high-mass-ratio mergers, in which the tidal disruption of the companion and its accretion onto the primary star determine prompt black hole formation. As a result of the tidal disruption, an accretion disc of neutron-rich and cold matter forms with baryon masses ∼0.15M⊙, and it is significantly heavier than the remnant discs in equal-masses prompt collapse mergers. Massive dynamical ejecta of order ∼0.01M⊙ also originate from the tidal disruption. They are neutron rich and expand from the orbital plane with a crescent-like geometry. Consequently, bright, red and temporally extended kilonova emission is predicted from these mergers. Our results show that prompt black hole mergers can power bright electromagnetic counterparts for high-mass-ratio binaries, and that the binary mass ratio can be in principle constrained from multimessenger observations.

Over the past three years, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. 2021 indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible with the emergence of a new component due to non-thermal emission from the fast tail of the dynamical ejecta of ab-initio binary neutron star (BNS) merger simulations. This provides a new avenue to constrain binary parameters. Specifically, we find that equal mass models with soft equation of state (EOS) and high mass ratio models with stiff EOS are disfavored as they typically predict afterglows that peak too early to explain the recent observations. Moderate stiffness and mass ratio models, instead, tend to be in good overall agreement with the data.

We present fitting formulae for the dynamical ejecta properties and remnant disk masses from the largest to date sample of numerical relativity simulations. The considered data include some of the latest simulations with microphysical nuclear equations of state (EOS) and neutrino transport as well as other results with polytropic EOS available in the literature. Our analysis indicates that the broad features of the dynamical ejecta and disk properties can be captured by fitting expressions, that depend on mass ratio and reduced tidal parameter. The comparative analysis of literature data shows that microphysics and neutrino absorption have a significant impact on the dynamical ejecta properties. Microphysical nuclear EOS lead to average velocities smaller than polytropic EOS, while including neutrino absorption results in larger average ejecta masses and electron fractions. Hence, microphysics and neutrino transport are necessary to obtain quantitative models of the ejecta in terms of the binary parameters.