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 ∼ 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.
.
The first neutron star (NS) merger observed by advanced LIGO and Virgo, GW170817, and its fireworks of electromagnetic counterparts across the entire electromagnetic spectrum marked the beginning ...of multi-messenger astronomy and astrophysics with gravitational waves. The ultraviolet, optical, and near-infrared emission was consistent with being powered by the radioactive decay of nuclei synthesized in the merger ejecta by the rapid neutron capture process (r-process). Starting from an outline of the inferred properties of this “kilonova” emission, I discuss possible astrophysical sites for r-process nucleosynthesis in NS mergers, arguing that the heaviest r-process elements synthesized in this event most likely originated in outflows from a post-merger accretion disk. I compare the inferred properties of r-process element production in GW170817 to current observational constraints on galactic heavy r-process nucleosynthesis and discuss challenges merger-only models face in explaining the r-process content of our galaxy. Based on the observational properties of GW170817 and recent theoretical progress on r-process nucleosynthesis in collapsars, I then show how GW170817 points to collapsars as the dominant source of r-process enrichment in the Milky Way. These rare core-collapse events arguably better satisfy existing constraints and overcome problems related to r-process enrichment in various environments that NS mergers face. Finally, I comment on the universality of the r-process and on how variations in light r-process elements can be obtained both in NS mergers and collapsars.
Merging binaries consisting of two neutron stars (NSs) or an NS and a stellar-mass black hole typically form a massive accretion torus around the remnant black hole or long-lived NS. Outflows from ...these neutrino-cooled accretion disks represent an important site for r-process nucleosynthesis and the generation of kilonovae. We present the first three-dimensional, general-relativistic magnetohydrodynamic (GRMHD) simulations including weak interactions and a realistic equation of state of such accretion disks over viscous timescales (380 ms). We witness the emergence of steady-state MHD turbulence, a magnetic dynamo with an ∼20 ms cycle, and the generation of a "hot" disk corona that launches powerful thermal outflows aided by the energy released as free nucleons recombine into -particles. We identify a self-regulation mechanism that keeps the midplane electron fraction low (Ye ∼ 0.1) over viscous timescales. This neutron-rich reservoir, in turn, feeds outflows that retain a sufficiently low value of Ye 0.2 to robustly synthesize third-peak r-process elements. The quasi-spherical outflows are projected to unbind 40% of the initial disk mass with typical asymptotic escape velocities of 0.1c and may thus represent the dominant mass ejection mechanism in NS-NS mergers. Including neutrino absorption, our findings agree with previous hydrodynamical -disk simulations that the entire range of r-process nuclei from the first to the third r-process peak can be synthesized in the outflows, in good agreement with observed solar system abundances. The asymptotic escape velocities and quantity of ejecta, when extrapolated to moderately higher disk masses, are consistent with those needed to explain the red kilonova emission following the NS merger GW170817.
The production of elements by rapid neutron capture (r-process) in neutron-star mergers is expected theoretically and is supported by multimessenger observations
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
delayed outflows from a remnant accretion disk formed around the newly born black hole
as the dominant source of heavy r-process material from that event
. 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
. 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
, 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.
The merger of binary neutron stars, or of a neutron star and a stellar-mass black hole, can result in the formation of a massive rotating torus around a spinning black hole. In addition to providing ...collimating media for γ-ray burst jets, unbound outflows from these disks are an important source of mass ejection and rapid neutron capture (r-process) nucleosynthesis. We present the first three-dimensional general-relativistic magnetohydrodynamic (GRMHD) simulations of neutrino-cooled accretion disks in neutron star mergers, including a realistic equation of state valid at low densities and temperatures, self-consistent evolution of the electron fraction, and neutrino cooling through an approximate leakage scheme. After initial magnetic field amplification by magnetic winding, we witness the vigorous onset of turbulence driven by the magnetorotational instability (MRI). The disk quickly reaches a balance between heating from MRI-driven turbulence and neutrino cooling, which regulates the midplane electron fraction to a low equilibrium value Y_{e}≈0.1. Over the 380-ms duration of the simulation, we find that a fraction ≈20% of the initial torus mass is unbound in powerful outflows with asymptotic velocities v≈0.1c and electron fractions Y_{e}≈0.1-0.25. Postprocessing the outflows through a nuclear reaction network shows the production of a robust second- and third-peak r process. Though broadly consistent with the results of previous axisymmetric hydrodynamical simulations, extrapolation of our results to late times suggests that the total ejecta mass from GRMHD disks is significantly higher. Our results provide strong evidence that postmerger disk outflows are an important site for the r process.
ABSTRACT Binary neutron star (BNS) mergers are the leading model to explain the phenomenology of short gamma-ray bursts (SGRBs). Recent observations of long-lasting X-ray afterglows of SGRBs ...challenge standard paradigms and indicate that in a large fraction of events a long-lived neutron star (NS) may be formed rather than a black hole. Understanding the mechanisms underlying these afterglows is necessary in order to address the open questions concerning the nature of SGRB central engines. However, recent theoretical progress has been hampered by the fact that the timescales of interest for the afterglow emission are inaccessible to numerical relativity simulations. Here we present a detailed model to bridge the gap between numerical simulations of the merger process and the relevant timescales for the afterglows, assuming that the merger results in a long-lived NS. This model is formulated in terms of a set of coupled differential equations that follow the evolution of the post-merger system and predict its electromagnetic (EM) emission in a self-consistent way, starting from initial data that can be extracted from BNS merger simulations. The model presented here also allows us to search for suitable EM counterparts for multimessenger astronomy, which is expected to become reality within the next few years thanks to ground-based GW detectors such as advanced LIGO and Virgo. This paper discusses the formulation and implementation of the model. In a companion paper, we employ this model to predict the EM emission from to after a BNS merger and discuss the implications in the context of SGRBs and multimessenger astronomy.
ABSTRACT Recent observations indicate that in a large fraction of binary neutron star (BNS) mergers a long-lived neutron star (NS) may be formed rather than a black hole. Unambiguous electromagnetic ...(EM) signatures of such a scenario would strongly impact our knowledge on how short gamma-ray bursts (SGRBs) and their afterglow radiation are generated. Furthermore, such EM signals would have profound implications for multimessenger astronomy with joint EM and gravitational-wave (GW) observations of BNS mergers, which will soon become reality thanks to the ground-based advanced LIGO/Virgo GW detector network. Here we explore such EM signatures based on the model presented in a companion paper, which provides a self-consistent evolution of the post-merger system and its EM emission up to ∼107 s. Light curves and spectra are computed for a wide range of post-merger physical properties. We present X-ray afterglow light curves corresponding to the "standard" and the "time-reversal" scenario for SGRBs (prompt emission associated with the merger or with the collapse of the long-lived NS). The light curve morphologies include single and two-plateau features with timescales and luminosities that are in good agreement with Swift observations. Furthermore, we compute the X-ray signal that should precede the SGRB in the time-reversal scenario, the detection of which would represent smoking-gun evidence for this scenario. Finally, we find a bright, highly isotropic EM transient peaking in the X-ray band at ∼102-104 s after the BNS merger with luminosities of LX ∼ 1046-1048 erg s−1. This signal represents a very promising EM counterpart to the GW emission from BNS mergers.
Abstract
Fast neutron-rich material ejected dynamically over ≲10 ms during the merger of a binary neutron star (BNS) can give rise to distinctive electromagnetic counterparts to the system’s ...gravitational-wave emission that serve as a “smoking gun” to distinguish between a BNS and an NS–black hole merger. We present novel ab initio modeling of the kilonova precursor and kilonova afterglow based on 3D general-relativistic magnetohydrodynamic simulations of BNS mergers with nuclear, tabulated, finite-temperature equations of state (EOSs), weak interactions, and approximate neutrino transport. We analyze dynamical mass ejection from 1.35–1.35
M
⊙
binaries, consistent with properties of the first observed BNS merger GW170817, using three nuclear EOSs that span the range of allowed compactness of 1.35
M
⊙
-neutron stars. Nuclear reaction network calculations yield a robust second-to-third-peak
r
-process. We find few ×10
−6
M
⊙
of fast (
v
> 0.6
c
) ejecta that give rise to broadband synchrotron emission on ∼years timescales, consistent with tentative evidence for excess X-ray/radio emission following GW170817. We find ≈2 × 10
−5
M
⊙
of free neutrons that power a kilonova precursor on ≲ hours timescale. A boost in early UV/optical brightness by a factor of a few due to previously neglected relativistic effects, with enhancements up to ≲10 hr post-merger, is promising for future detection with UV/optical telescopes like Swift or ULTRASAT. We find that a recently predicted opacity boost due to highly ionized lanthanides at ≳70,000 K is unlikely to affect the early kilonova based on the obtained ejecta structures. Azimuthal inhomogeneities in dynamical ejecta composition for soft EOSs found here (“lanthanide/actinide pockets”) may have observable consequences for both early kilonova and late-time nebular emission.
Short gamma-ray bursts (SGRBs) are among the most luminous explosions in the universe and their origin still remains uncertain. Observational evidence favors the association with binary neutron star ...or neutron star-black hole (NS-BH) binary mergers. Leading models relate SGRBs to a relativistic jet launched by the BH-torus system resulting from the merger. However, recent observations have revealed a large fraction of SGRB events accompanied by X-ray afterglows with durations ~10 super(2)-10 super(5) s, suggesting continuous energy injection from a long-lived central engine, which is incompatible with the short (<, ~1 s) accretion timescale of a BH-torus system. The formation of a supramassive NS, resisting the collapse on much longer spin-down timescales, can explain these afterglow durations, but leaves serious doubts on whether a relativistic jet can be launched at the merger. Here we present a novel scenario accommodating both aspects, where the SGRB is produced after the collapse of a supramassive NS. Early differential rotation and subsequent spin-down emission generate an optically thick environment around the NS consisting of a photon-pair nebula and an outer shell of baryon-loaded ejecta. While the jet easily drills through this environment, spin-down radiation diffuses outward on much longer timescales and accumulates a delay that allows the SGRB to be observed before (part of) the long-lasting X-ray signal. By analyzing diffusion timescales for a wide range of physical parameters, we find delays that can generally reach ~10 super(5) s, compatible with observations. The success of this fundamental test makes this "time-reversal" scenario an attractive alternative to current SGRB models.
Abstract
The merger of two neutron stars or a neutron star and a black hole typically results in the formation of a postmerger accretion disk. Outflows from disks may dominate the overall ejecta from ...mergers and be a major source of
r
-process nuclei in our universe. We explore the parameter space of such disks and their outflows and
r
-process yields by performing 3D general-relativistic magnetohydrodynamic simulations with weak interactions and approximate neutrino transport. We discuss the mapping between the initial binary parameters and the parameter space of the resulting disks, chiefly characterized by their initial accretion rate. We demonstrate the existence of an ignition threshold for weak interactions at around ∼10
−3
M
⊙
s
−1
for typical parameters by means of analytic calculations and numerical simulations. We find a degenerate, self-regulated, neutrino-cooled regime above the threshold and an advection-dominated regime below the threshold. Excess heating in the absence of neutrino cooling below the threshold leads to ≳60% of the initial disk mass being ejected in outflows, with typical velocities of ∼(0.1–0.2)
c
, compared to ≲40% at ∼(0.1–0.15)
c
above the threshold. While disks below the threshold show suppressed production of light
r
-process elements, disks above the threshold can produce the entire range of
r
-process elements, in good agreement with the observed solar system abundances. Disks below the ignition threshold may produce an overabundance of actinides seen in actinide-boost stars. As gravitational-wave detectors start to sample the neutron star merger parameter space, different disk realizations may be observable via their associated kilonova emission.