An official website of the United States government
Abstract The rapid neutron capture process (r-process) is one of the main mechanisms whereby elements heavier than iron are synthesized, and is entirely responsible for the natural production of the actinides. Kilonova emissions are modeled as being largely powered by the radioactive decay of species synthesized via ther-process. Given that ther-process occurs far from nuclear stability, unmeasured beta-decay rates play an essential role in setting the timescale for ther-process. In an effort to better understand the sensitivity of kilonova modeling to different theoretical global beta-decay descriptions, we incorporate these into nucleosynthesis calculations. We compare the results of these calculations and highlight differences in kilonova nuclear energy generation and light-curve predictions, as well as final abundances and their implications for nuclear cosmochronometry. We investigate scenarios where differences in beta-decay rates are responsible for increased nuclear heating on timescales of days that propagates into a significantly increased average bolometric luminosity between 1 and 10 days post-merger. We identify key nuclei, both measured and unmeasured, whose decay rates directly impact nuclear heating generation on timescales responsible for light-curve evolution. We also find that uncertainties in beta-decay rates significantly impact age estimates from cosmochronometry.
more »
« less
Element Formation in Radiation-hydrodynamics Simulations of Kilonovae
Abstract Understanding the details ofr-process nucleosynthesis in binary neutron star merger (BNSM) ejecta is key to interpreting kilonova observations and identifying the role of BNSMs in the origin of heavy elements. We present a self-consistent, two-dimensional, ray-by-ray radiation-hydrodynamic evolution of BNSM ejecta with an online nuclear network (NN) up to a timescale of days. For the first time, an initial numerical relativity ejecta profile composed of the dynamical component and spiral-wave and disk winds is evolved including detailedr-process reactions and nuclear heating effects. A simple model for the jet energy deposition is also included. Our simulation highlights that the common approach of relating in postprocessing the final nucleosynthesis yields to the initial thermodynamic profile of the ejecta can lead to inaccurate predictions. Moreover, we find that neglecting the details of the radiation-hydrodynamic evolution of the ejecta in nuclear calculations can introduce deviations of up to 1 order of magnitude in the final abundances of several elements, including very light and secondr-process peak elements. The presence of a jet affects element production only in the innermost part of the polar ejecta, and it does not alter the global nucleosynthesis results. Overall, our analysis shows that employing an online NN improves the reliability of nucleosynthesis and kilonova light-curve predictions.
more »
« less
- Award ID(s):
- 2011725
- PAR ID:
- 10546435
- Publisher / Repository:
- DOI PREFIX: 10.3847
- Date Published:
- Journal Name:
- The Astrophysical Journal Letters
- Volume:
- 974
- Issue:
- 1
- ISSN:
- 2041-8205
- Format(s):
- Medium: X Size: Article No. L5
- Size(s):
- Article No. L5
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We study the production of very light elements (Z< 20) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition-dependent nuclear equations of state, and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, ∼0.5–2 × 10−6M⊙of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10−6M⊙of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavyr-process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including nonlocal thermodynamic equilibrium effects. The latter could be crucial to observe helium lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long-lived remnant) the one required to explain early spectral features in the kilonova of GW170817.more » « less
-
ABSTRACT We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.more » « less
-
Abstract We present nucleosynthesis and light-curve predictions for a new site of the rapid neutron capture process (r-process) from magnetar giant flares (GFs). Motivated by observations indicating baryon ejecta from GFs, J. Cehula et al. proposed that mass ejection occurs after a shock is driven into the magnetar crust during the GF. We confirm using nuclear reaction network calculations that these ejecta synthesize moderate yields of third-peakr-process nuclei and more substantial yields of lighterr-nuclei, while leaving a sizable abundance of free neutrons in the outermost fastest expanding ejecta layers. The finalr-process mass fraction and distribution are sensitive to the relative efficiencies ofα-capture andn-capture freeze-outs. We use our nucleosynthesis output in a semianalytic model to predict the light curves of novae breves, the transients following GFs powered by radioactive decay. For a baryonic ejecta mass similar to that inferred of the 2004 Galactic GF from SGR 1806-20, we predict a peak UV/optical luminosity of ∼1039–1040erg s−1at ∼10–15 minutes, rendering such events potentially detectable to several Mpc following a gamma-ray trigger by wide-field transient monitors such as ULTRASAT/UVEX. The peak luminosity and timescale of the transient increase with the GF strength due to the larger ejecta mass. Although GFs likely contribute 1%–10% of the total Galacticr-process budget, their short delay-times relative to star formation make them an attractive source to enrich the earliest generations of stars.more » « less
-
Abstract This study investigates the origins of GW230529, delving into its formation from massive stars within isolated binary systems. Utilizing population-synthesis models, we present compelling evidence that the neutron star component forms second. However, the event’s low signal-to-noise ratio introduces complexities in identifying the underlying physical mechanisms driving its formation. Augmenting our analysis with insights from numerical relativity, we estimate the final black hole mass and spin to be approximately 5.3M⊙and 0.53, respectively. Furthermore, we employ the obtained posterior samples to calculate the ejecta mass and kilonova light curves resulting fromr-process nucleosynthesis. We find the ejecta mass to be within 0–0.06M⊙, contingent on the neutron star equation of state. The peak brightness of the kilonova light curves indicates that targeted follow-up observations with a Rubin-like observatory may have detected this emission.more » « less
