Jets from Neutron-Star Merger Remnants and Massive Blue Kilonovae

We perform three-dimensional general-relativistic magnetohydrodynamic simulations with weak interactions of binary neutron-star (BNS) mergers resulting in a long-lived remnant neutron star, with properties typical of galactic BNS and consistent with those inferred for the first observed BNS merger GW170817. We demonstrate self-consistently that within ≲30  ms postmerger magnetized (σ∼5-10) incipient jets emerge with asymptotic Lorentz factor Γ∼5-10, which successfully break out from the merger debris within ≲20  ms. A fast (v≲0.6c), magnetized (σ∼0.1) wind surrounds the jet core and generates a UV/blue kilonova precursor on timescales of hours, similar to the precursor signal due to free neutron decay in fast dynamical ejecta. Postmerger ejecta are quickly dominated by magnetohydrodynamically driven outflows from an accretion disk. We demonstrate that, within only 50 ms postmerger, ≳2×10^{-2}M_{⊙} of lanthanide-free, quasispherical ejecta with velocities ∼0.1-0.2c is launched, yielding a kilonova signal consistent with GW170817 on timescales of ≲5  d.

Self-Consistent Picture of the Mass Ejection from a One Second Long Binary Neutron Star Merger Leaving a Short-Lived Remnant in a General-Relativistic Neutrino-Radiation Magnetohydrodynamic Simulation

We perform a general-relativistic neutrino-radiation magnetohydrodynamic simulation of a one second-long binary neutron star merger on the Japanese supercomputer Fugaku using about 85 million CPU hours with 20 736 CPUs. We consider an asymmetric binary neutron star merger with masses of 1.2M_{⊙} and 1.5M_{⊙} and a "soft" equation of state SFHo. It results in a short-lived remnant with the lifetime of ≈0.017  s, and subsequent massive torus formation with the mass of ≈0.05M_{⊙} after the remnant collapses to a black hole. For the first time, we find that after the dynamical mass ejection, which drives the fast tail and mildly relativistic components, the postmerger mass ejection from the massive torus takes place due to the magnetorotational instability-driven turbulent viscosity in a single simulation and the two ejecta components are seen in the distributions of the electron fraction and velocity with distinct features.

Neutrino transport in general relativistic neutron star merger simulations

Numerical simulations of neutron star-neutron star and neutron star-black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino-matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Gigaelectronvolt emission from a compact binary merger

An energetic γ-ray burst (GRB), GRB 211211A, was observed on 11 December 20211,2. Despite its long duration, typically associated with bursts produced by the collapse of massive stars, the observation of an optical-infrared kilonova points to a compact binary merger origin3. Here we report observations of a significant (more than five sigma) transient-like emission in the high-energy γ-rays of GRB 211211A (more than 0.1 gigaelectronvolts) starting 103 seconds after the burst. After an initial phase with a roughly constant flux (about 5 × 10-10 erg per second per square centimetre) lasting about 2 × 104 seconds, the flux started decreasing and soon went undetected. Our detailed modelling of public and dedicated multi-wavelength observations demonstrates that gigaelectronvolt emission from GRB 211211A is in excess with respect to the flux predicted by the state-of-the-art afterglow model at such late time. We explore the possibility that the gigaelectronvolt excess is inverse Compton emission owing to the interaction of a late-time, low-power jet with an external source of photons, and find that kilonova emission can provide the seed photons. Our results open perspectives for observing binary neutron star mergers.

Thinking Outside the Box: Numerical Relativity with Particles

The observation of gravitational waves from compact objects has now become an active part of observational astronomy. For a sound interpretation, one needs to compare such observations against detailed Numerical Relativity simulations, which are essential tools to explore the dynamics and physics of compact binary mergers. To date, essentially all simulation codes that solve the full set of Einstein’s equations are performed in the framework of Eulerian hydrodynamics. The exception is our recently developed Numerical Relativity code SPHINCS_BSSN which solves the commonly used BSSN formulation of the Einstein equations on a structured mesh and the matter equations via Lagrangian particles. We show here, for the first time, SPHINCS_BSSN neutron star merger simulations with piecewise polytropic approximations to four nuclear matter equations of state. In this set of neutron star merger simulations, we focus on perfectly symmetric binary systems that are irrotational and have 1.3 M⊙ masses. We introduce some further methodological refinements (a new way of steering dissipation, an improved particle–mesh mapping), and we explore the impact of the exponent that enters in the calculation of the thermal pressure contribution. We find that it leaves a noticeable imprint on the gravitational wave amplitude (calculated via both quadrupole approximation and the Ψ4 formalism) and has a noticeable impact on the amount of dynamic ejecta. Consistent with earlier findings, we only find a few times 10−3M⊙ as dynamic ejecta in the studied equal mass binary systems, with softer equations of state (which are more prone to shock formation) ejecting larger amounts of matter. In all of the cases, we see a credible high-velocity (∼0.5…0.7c) ejecta component of ∼10−4M⊙ that is launched at contact from the interface between the two neutron stars. Such a high-velocity component has been suggested to produce an early, blue precursor to the main kilonova emission, and it could also potentially cause a kilonova afterglow.

Neutrino Fast Flavor Conversions in Neutron-Star Postmerger Accretion Disks

A compact accretion disk may be formed in the merger of two neutron stars or of a neutron star and a stellar-mass black hole. Outflows from such accretion disks have been identified as a major site of rapid neutron-capture (r-process) nucleosynthesis and as the source of "red" kilonova emissions following the first observed neutron-star merger GW170817. We present long-term general-relativistic radiation magnetohydrodynamic simulations of a typical postmerger accretion disk at initial accretion rates of M[over ˙]∼1  M_{⊙} s^{-1} over 400 ms postmerger. We include neutrino radiation transport that accounts for the effects of neutrino fast flavor conversions dynamically. We find ubiquitous flavor oscillations that result in a significantly more neutron-rich outflow, providing lanthanide and 3rd-peak r-process abundances similar to solar abundances. This provides strong evidence that postmerger accretion disks are a major production site of heavy r-process elements. A similar flavor effect may allow for increased lanthanide production in collapsars.

Electron capture in stars

Electron capture on nuclei plays an essential role in the dynamics of several astrophysical objects, including core-collapse and thermonuclear supernovae, the crust of accreting neutron stars in binary systems and the final core evolution of intermediate-mass stars. In these astrophysical objects, the capture occurs at finite temperatures and densities, at which the electrons form a degenerate relativistic electron gas. The capture rates can be derived from perturbation theory, where allowed nuclear transitions [Gamow-Teller (GT) transitions] dominate, except at the higher temperatures achieved in core-collapse supernovae, where forbidden transitions also contribute significantly to the capture rates. There has been decisive progress in recent years in measuring GT strength distributions using novel experimental techniques based on charge-exchange reactions. These measurements not only provide data for the GT distributions of ground states for many relevant nuclei, but also serve as valuable constraints for nuclear models which are needed to derive the capture rates for the many nuclei for which no data yet exist. In particular, models are needed to evaluate stellar capture rates at finite temperatures, where capture can also occur on nuclei in thermally excited states. There has also been significant progress in recent years in the modeling of stellar capture rates. This has been made possible by advances in nuclear many-body models as well as in computer soft- and hardware. Specifically, to derive reliable capture rates for core-collapse supernovae, a dedicated strategy has been developed based on a hierarchy of nuclear models specifically adapted to the abundant nuclei and astrophysical conditions present under various collapse conditions. In particular, for the challenging conditions where the electron chemical potential and the nuclearQvalues are of the same order, large-scale shell-model diagonalization calculations have proved to be an appropriate tool to derive stellar capture rates, often validated by experimental data. Such situations are relevant in the early stage of the core collapse of massive stars, for the nucleosynthesis of thermonuclear supernovae, and for the final evolution of the cores of intermediate-mass stars involving nuclei in the mass rangeA∼ 20-65. This manuscript reviews the experimental and theoretical progress recently achieved in deriving stellar electron capture rates. It also discusses the impact these improved rates have on our understanding of the various astrophysical objects.

Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model

Magnetohydrodynamic (MHD) turbulence in neutron star (NS) merger remnants can impact their evolution and multi-messenger signatures, complicating the interpretation of present and future observations. Due to the high Reynolds numbers and the large computational costs of numerical relativity simulations, resolving all the relevant scales of the turbulence will be impossible for the foreseeable future. Here, we adopt a method to include subgrid-scale turbulence in moderate resolution simulations by extending the large-eddy simulation (LES) method to general relativity (GR). We calibrate our subgrid turbulence model with results from very-high-resolution GRMHD simulations, and we use it to perform NS merger simulations and study the impact of turbulence. We find that turbulence has a quantitative, but not qualitative, impact on the evolution of NS merger remnants, on their gravitational wave signatures, and on the outflows generated in binary NS mergers. Our approach provides a viable path to quantify uncertainties due to turbulence in NS mergers.

Asymmetric mass ratios for bright double neutron-star mergers

The discovery of a radioactively powered kilonova associated with the binary neutron-star merger GW170817 remains the only confirmed electromagnetic counterpart to a gravitational-wave event^1,2. Observations of the late-time electromagnetic emission, however, do not agree with the expectations from standard neutron-star merger models. Although the large measured ejecta mass^3,4 could be explained by a progenitor system that is asymmetric in terms of the stellar component masses (that is, with a mass ratio q of 0.7 to 0.8)⁵, the known Galactic population of merging double neutron-star systems (that is, those that will coalesce within billions of years or less) has until now consisted only of nearly equal-mass (q > 0.9) binaries⁶. The pulsar PSR J1913+1102 is a double system in a five-hour, low-eccentricity (0.09) orbit, with an orbital separation of 1.8 solar radii⁷, and the two neutron stars are predicted to coalesce in [Formula: see text] million years owing to gravitational-wave emission. Here we report that the masses of the pulsar and the companion neutron star, as measured by a dedicated pulsar timing campaign, are 1.62 ± 0.03 and 1.27 ± 0.03 solar masses, respectively. With a measured mass ratio of q = 0.78 ± 0.03, this is the most asymmetric merging system reported so far. On the basis of this detection, our population synthesis analysis implies that such asymmetric binaries represent between 2 and 30 per cent (90 per cent confidence) of the total population of merging binaries. The coalescence of a member of this population offers a possible explanation for the anomalous properties of GW170817, including the observed kilonova emission from that event.

Kilonovae

The coalescence of double neutron star (NS-NS) and black hole (BH)-NS binaries are prime sources of gravitational waves (GW) for Advanced LIGO/Virgo and future ground-based detectors. Neutron-rich matter released from such events undergoes rapid neutron capture (r-process) nucleosynthesis as it decompresses into space, enriching our universe with rare heavy elements like gold and platinum. Radioactive decay of these unstable nuclei powers a rapidly evolving, approximately isotropic thermal transient known as a "kilonova", which probes the physical conditions during the merger and its aftermath. Here I review the history and physics of kilonovae, leading to the current paradigm of day-timescale emission at optical wavelengths from lanthanide-free components of the ejecta, followed by week-long emission with a spectral peak in the near-infrared (NIR). These theoretical predictions, as compiled in the original version of this review, were largely confirmed by the transient optical/NIR counterpart discovered to the first NS-NS merger, GW170817, discovered by LIGO/Virgo. Using a simple light curve model to illustrate the essential physical processes and their application to GW170817, I then introduce important variations about the standard picture which may be observable in future mergers. These include ∼ hour-long UV precursor emission, powered by the decay of free neutrons in the outermost ejecta layers or shock-heating of the ejecta by a delayed ultra-relativistic outflow; and enhancement of the luminosity from a long-lived central engine, such as an accreting BH or millisecond magnetar. Joint GW and kilonova observations of GW170817 and future events provide a new avenue to constrain the astrophysical origin of the r-process elements and the equation of state of dense nuclear matter.

Binary Neutron Star (BNS) Merger: What We Learned from Relativistic Ejecta of GW/GRB 170817A

Gravitational Waves (GW) from coalescence of a Binary Neutron Star (BNS) and its accompanying short Gamma-Ray Burst (GRB) GW/GRB 170817A confirmed the presumed origin of these puzzling transients and opened up the way for relating properties of short GRBs to those of their progenitor stars and their surroundings. Here we review an extensive analysis of the prompt gamma-ray and late afterglows of this event. We show that a fraction of polar ejecta from the merger had been accelerated to ultra-relativistic speeds. This structured jet had an initial Lorentz factor of about 260 in our direction, which was O ( 10 ∘ ) from the jet’s axis, and was a few orders of magnitude less dense than in typical short GRBs. At the time of arrival to circum-burst material the ultra-relativistic jet had a close to Gaussian profile and a Lorentz factor ≳ 130 in its core. It had retained in some extent its internal collimation and coherence, but had extended laterally to create mildly relativistic lobes—a cocoon. Its external shocks on the far from center inhomogeneous circum-burst material and low density of colliding shells generated slowly rising afterglows, which peaked more than 100 days after the prompt gamma-ray. The circum-burst material was somehow correlated with the merger. As non-relativistic outflows or tidally ejected material during BNS merger could not have been arrived to the location of the external shocks before the relativistic jet, circum-burst material might have contained recently ejected materials from resumption of internal activities, faulting and mass loss due to deformation and breaking of stars crusts by tidal forces during latest stages of their inspiral but well before their merger. By comparing these findings with the results of relativistic Magneto-Hydro-Dynamics (MHD) simulations and observed gravitational waves we conclude that progenitor neutron stars were most probably old, had close masses and highly reduced magnetic fields.

Collapsars as a major source of r-process elements

The production of elements by rapid neutron capture (r-process) in neutron-star mergers is expected theoretically and is supported by multimessenger observations^1-3 of gravitational-wave event GW170817: this production route is in principle sufficient to account for most of the r-process elements in the Universe⁴. Analysis of the kilonova that accompanied GW170817 identified^5,6 delayed outflows from a remnant accretion disk formed around the newly born black hole^7-10 as the dominant source of heavy r-process material from that event^9,11. Similar accretion disks are expected to form in collapsars (the supernova-triggering collapse of rapidly rotating massive stars), which have previously been speculated to produce r-process elements^12,13. Recent observations of stars rich in such elements in the dwarf galaxy Reticulum II¹⁴, as well as the Galactic chemical enrichment of europium relative to iron over longer timescales^15,16, are more consistent with rare supernovae acting at low stellar metallicities than with neutron-star mergers. Here we report simulations that show that collapsar accretion disks yield sufficient r-process elements to explain observed abundances in the Universe. Although these supernovae are rarer than neutron-star mergers, the larger amount of material ejected per event compensates for the lower rate of occurrence. We calculate that collapsars may supply more than 80 per cent of the r-process content of the Universe.