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.

The AT2017gfo kilonova counterpart of the binary neutron star merger event GW170817 was characterized by an early-time bright peak in optical and UV bands. Such blue kilonova is commonly interpreted as a signature of weak r-process nucleosynthesis in a fast expanding wind whose origin is currently debated. Numerical-relativity simulations with microphysical equations of state, approximate neutrino transport, and turbulent viscosity reveal a new mechanism that can power the blue kilonova. Spiral density waves in the remnant generate a characteristic wind of mass ∼10 −2 M and velocity ∼0.2c. The ejected material has electron fraction mostly distributed above 0.25 being partially reprocessed by hydrodynamic shocks in the expanding arms. The combination of dynamical ejecta and spiral-wave wind can account for solar system abundances of r-process elements and early-time observed light curves.PACS numbers: 04.25.D-, 04.30.Db, 95.30.Sf, 95.30.Lz, 97.60.Jd arXiv:1907.04872v1 [astro-ph.HE]

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.

For the first ∼3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n < 0.01 cm−3. We report on observational evidence for an excess of X-ray emission at δt > 900 days after the merger. With L x ≈ 5 × 1038 erg s−1 at 1234 days, the recently detected X-ray emission represents a ≥3.2σ (Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance ≥3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented.