Context. The astrophysical r-process site where about half of the elements, heavier than iron are produced, has been a puzzle for several decades. Here we discuss the role of one of the leading ideas ...– neutron star mergers (NSMs) – in the light of the first direct detection of such an event in both gravitational (GW) and electromagnetic (EM) waves. Aims. Our aim is to understand the implications of the first GW/EM observations of a NSM for cosmic nucleosynthesis. Methods. We analyse bolometric and NIR lightcurves of the first detected double NSM and compare them to nuclear reaction network-based macronova models. Results. The slope of the bolometric lightcurve is consistent with the radioactive decay of neutron star ejecta with Ye ≲ 0.3 (but not larger), which provides strong evidence for an r-process origin of the electromagnetic emission. This rules out in particular “nickel winds” as major source of the emission. We find that the NIR lightcurves can be well fitted either with or without lanthanide-rich ejecta. Our limits on the ejecta mass together with estimated rates directly confirm earlier purely theoretical or indirect observational conclusions that double neutron star mergers are indeed a major site of cosmic nucleosynthesis. If the ejecta mass was typical, NSMs can easily produce all of the estimated Galactic r-process matter, and – depending on the real rate – potentially even more. This could be a hint that the event ejected a particularly large amount of mass, maybe due to a substantial difference between the component masses. This would be compatible with the mass limits obtained from the GW-observation. Conclusions. The recent observations suggests that NSMs are responsible for a broad range of r-process nuclei and that they are at least a major, but likely the dominant r-process site in the Universe.
ABSTRACT
The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These ...observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy r-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound–bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion–opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and thus eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion–opacity formalism. The agreement between the line-binned and expansion–opacity results is explained as arising from the similarity in their opacities in the limit of low optical depth, where radiation transport is important in the ejecta. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semirelativistic approaches.
We follow the long-term evolution of the dynamic ejecta of neutron star mergers for up to 100 years and over a density range of roughly 40 orders of magnitude. We include the nuclear energy input ...from the freshly synthesized, radioactively decaying nuclei in our simulations and study its effects on the remnant dynamics. Although the nuclear heating substantially alters the long-term evolution, we find that running nuclear networks over purely hydrodynamic simulations (i.e. without heating) yields actually acceptable nucleosynthesis results. The main dynamic effect of the radioactive heating is to quickly smooth out inhomogeneities in the initial mass distribution, subsequently the evolution proceeds self-similarly and after 100 years the remnant still carries the memory of the initial binary mass ratio. We also explore the nucleosynthetic yields for two mass ejection channels. The dynamic ejecta very robustly produce 'strong' r-process elements with A > 130 with a pattern that is essentially independent of the details of the merging system. From a simple model we find that neutrino-driven winds yield 'weak' r-process contributions with 50 < A < 130 whose abundance patterns vary substantially between different merger cases. This is because their electron fraction, set by the ratio of neutrino luminosities, varies considerably from case to case. Such winds do not produce any 56Ni, but a range of radioactive isotopes that are long-lived enough to produce a second, radioactively powered electromagnetic transient in addition to the 'macronova' from the dynamic ejecta. While our wind model is very simple, it nevertheless demonstrates the potential of such neutrino-driven winds for electromagnetic transients and it motivates further, more detailed neutrino-hydrodynamic studies. The properties of the mentioned transients are discussed in more detail in a companion paper.
A long-standing paradigm in astrophysics is that collisions- or mergers-of two neutron stars form highly relativistic and collimated outflows (jets) that power Y-ray bursts of short (less than two ...seconds) duration. The observational support for this model, however, is only indirect. A hitherto outstanding prediction is that gravitational-wave events from such mergers should be associated with Y-ray bursts, and that a majority of these bursts should be seen off-axis, that is, they should point away from Earth. Here we report the discovery observations of the X-ray counterpart associated with the gravitational-wave event GW170817. Although the electromagnetic counterpart at optical and infrared frequencies is dominated by the radioactive glow (known as a 'kilonova') from freshly synthesized rapid neutron capture (r-process) material in the merger ejecta, observations at X-ray and, later, radio frequencies are consistent with a short Y-ray burst viewed off-axis. Our detection of X-ray emission at a location coincident with the kilonova transient provides the missing observational link between short Y-ray bursts and gravitational waves from neutron-star mergers, and gives independent confirmation of the collimated nature of the Y-ray-burst emission.
We study the optical and near-infrared luminosities and detectability of radioactively powered electromagnetic transients ('macronovae') occuring in the aftermath of binary neutron star and neutron ...star black hole mergers. We explore the transients that result from the dynamic ejecta and those from different types of wind outflows. Based on full nuclear network simulations we calculate the resulting light curves in different wavelength bands. We scrutinize the robustness of the results by comparing (a) two different nuclear reaction networks and (b) two macronova models. We explore in particular how sensitive the results are to the production of α-decaying trans-lead nuclei. We compare two frequently used mass models: the finite-range Droplet model (FRDM) and the nuclear mass model of Duflo and Zuker (DZ31). We find that the abundance of α-decaying trans-lead nuclei has a significant impact on the observability of the resulting macronovae. For example, the DZ31 model yields considerably larger abundances resulting in larger heating rates and thermalization efficiencies and therefore predicts substantially brighter macronova transients. We find that the dynamic ejecta from NSNS models can reach peak K-band magnitudes in excess of −15 while those from NSBH cases can reach beyond −16. Similar values can be reached by some of our wind models. Several of our models (both wind and dynamic ejecta) yield properties that are similar to the transient that was observed in the aftermath of the short GRB 130603B. We further explore the expected macronova detection frequencies for current and future instruments such as VISTA, ZTF and LSST.
We report the discovery and monitoring of the near-infrared counterpart (AT2017gfo) of a binary neutron-star merger event detected as a gravitational wave source by Advanced Laser Interferometer ...Gravitational-wave Observatory (LIGO)/Virgo (GW170817) and as a short gamma-ray burst by Fermi Gamma-ray Burst Monitor (GBM) and Integral SPI-ACS (GRB 170817A). The evolution of the transient light is consistent with predictions for the behavior of a "kilonova/macronova" powered by the radioactive decay of massive neutron-rich nuclides created via r-process nucleosynthesis in the neutron-star ejecta. In particular, evidence for this scenario is found from broad features seen in Hubble Space Telescope infrared spectroscopy, similar to those predicted for lanthanide-dominated ejecta, and the much slower evolution in the near-infrared K s -band compared to the optical. This indicates that the late-time light is dominated by high-opacity lanthanide-rich ejecta, suggesting nucleosynthesis to the third r-process peak (atomic masses A 195 ). This discovery confirms that neutron-star mergers produce kilo-/macronovae and that they are at least a major-if not the dominant-site of rapid neutron capture nucleosynthesis in the universe.
ABSTRACT We present a comprehensive nucleosynthesis study of the neutrino-driven wind in the aftermath of a binary neutron star merger. Our focus is the initial remnant phase when a massive central ...neutron star is present. Using tracers from a recent hydrodynamical simulation, we determine total masses and integrated abundances to characterize the composition of unbound matter. We find that the nucleosynthetic yields depend sensitively on both the life time of the massive neutron star and the polar angle. Matter in excess of up to 9 × 10−3 M becomes unbound until ∼200 ms. Due to electron fractions of Ye 0.2-0.4, mainly nuclei with mass numbers A < 130 are synthesized, complementing the yields from the earlier dynamic ejecta. Mixing scenarios with these two types of ejecta can explain the abundance pattern in r-process enriched metal-poor stars. Additionally, we calculate heating rates for the decay of the freshly produced radioactive isotopes. The resulting light curve peaks in the blue band after about 4 hr. Furthermore, high opacities due to heavy r-process nuclei in the dynamic ejecta lead to a second peak in the infrared after 3-4 days.
We have explored the Eu production in the Milky Way by means of a very detailed chemical evolution model. In particular, we have assumed that Eu is formed in merging neutron star (or neutron ...star-black hole) binaries as well as in Type II supernovae. We have tested the effects of several important parameters influencing the production of Eu during the merging of two neutron stars, such as (i) the time-scale of coalescence, (ii) the Eu yields and (iii) the range of initial masses for the progenitors of the neutron stars. The yields of Eu from Type II supernovae are very uncertain, more than those from coalescing neutron stars, so we have explored several possibilities. We have compared our model results with the observed rate of coalescence of neutron stars, the solar Eu abundance, the Eu/Fe versus Fe/H relation in the solar vicinity and the Eu/H gradient along the Galactic disc. Our main results can be summarized as follows: (i) neutron star mergers can be entirely responsible for the production of Eu in the Galaxy if the coalescence time-scale is no longer than 1 Myr for the bulk of binary systems, the Eu yield is around 3 × 10−7 M and the mass range of progenitors of neutron stars is 9-50 M; (ii) both Type II supernovae and merging neutron stars can produce the right amount of Eu if the neutron star mergers produce 2 × 10−7 M per system and Type II supernovae, with progenitors in the range 20-50 M, produce yields of Eu of the order of 10−8-10−9 M; (iii) either models with only neutron stars producing Eu or mixed ones can reproduce the observed Eu abundance gradient along the Galactic disc.
We present a detailed, three-dimensional hydrodynamic study of the neutrino-driven winds emerging from the remnant of a neutron star merger. Our simulations are performed with the Newtonian, Eulerian ...code fish, augmented by a detailed, spectral neutrino leakage scheme that accounts for neutrino absorption. Consistent with earlier two-dimensional studies, a strong baryonic wind is blown out along the original binary rotation axis within ≈100 ms. From this model, we compute a lower limit on the expelled mass of 3.5 × 10−3 M⊙, relevant for heavy element nucleosynthesis. Because of stronger neutrino irradiation, the polar regions show substantially larger electron fractions than those at lower latitudes. The polar ejecta produce interesting r-process contributions from A ≈ 80 to about 130, while the more neutron-rich, lower latitude parts produce elements up to the third r-process peak near A ≈ 195. We calculate the properties of electromagnetic transients powered by the radioactivity in the wind, in addition to the ‘macronova’ transient stemming from the dynamic ejecta. The polar regions produce ultraviolet/optical transients reaching luminosities up to 1041 erg s−1, which peak around 1 d in optical and 0.3 d in bolometric luminosity. The lower latitude regions, due to their contamination with high-opacity heavy elements, produce dimmer and more red signals, peaking after ∼2 d in optical and infrared.
A Broad Grid of 2D Kilonova Emission Models Wollaeger, R. T.; Fryer, C. L.; Chase, E. A. ...
Astrophysical journal/The Astrophysical journal,
09/2021, Letnik:
918, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Abstract
Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The ...characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of two-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1
M
⊙
and 0.05 to 0.3
c
, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (
Y
e
) ejecta with a robust
r
-process composition and either a spherical or lobed high-
Y
e
ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code
SuperNu
using a full suite of lanthanide and fourth-row element opacities. We examine the trends of these models with parameter variation, show how they can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.
PDF
