skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Tracer Particles for Core-collapse Supernova Nucleosynthesis: The Advantages of Moving Backward
Abstract After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multidimensional simulations. However, obtaining nucleosynthesis yields based on such multidimensional core-collapse supernova simulations is not straightforward. Frequently, tracer particles are employed. Tracer particles may be tracked in situ during the simulation, but often they are reconstructed in a post-processing step based on the information saved during the hydrodynamic simulation. Reconstruction can be done in a number of ways, and here we compare the approaches of backward and forward integration of the equations of motion to the results based on inline particle trajectories. We find that both methods agree reasonably well with the inline results for isotopes for which a large number of particles contribute. However, for rarer isotopes that are produced only by a small number of particle trajectories, deviations can be large. For our setup, we find that backward integration leads to better agreement with the inline particles by more accurately reproducing the conditions following freeze-out from nuclear statistical equilibrium, because the establishment of nuclear statistical equilibrium erases the need for detailed trajectories at earlier times. Based on our results, if inline tracers are unavailable, we recommend backward reconstruction to the point when nuclear statistical equilibrium was last applied, with an interval between simulation snapshots of at most 1 ms for nucleosynthesis post-processing.  more » « less
Award ID(s):
1913531
PAR ID:
10421044
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
950
Issue:
1
ISSN:
0004-637X
Format(s):
Medium: X Size: Article No. 34
Size(s):
Article No. 34
Sponsoring Org:
National Science Foundation
More Like this
  1. Context. The γ -process nucleosynthesis in core-collapse supernovae is generally accepted as a feasible process for the synthesis of neutron-deficient isotopes beyond iron. However, crucial discrepancies between theory and observations still exist: the average yields of γ -process nucleosynthesis from massive stars are still insufficient to reproduce the solar distribution in galactic chemical evolution calculations, and the yields of the Mo and Ru isotopes are a factor of ten lower than the yields of the other γ -process nuclei. Aims. We investigate the γ -process in five sets of core-collapse supernova models published in the literature with initial masses of 15, 20, and 25 M ⊙ at solar metallicity. Methods. We compared the γ -process overproduction factors from the different models. To highlight the possible effect of nuclear physics input, we also considered 23 ratios of two isotopes close to each other in mass relative to their solar values. Further, we investigated the contribution of C–O shell mergers in the supernova progenitors as an additional site of the γ -process. Results. Our analysis shows that a large scatter among the different models exists for both the γ -process integrated yields and the isotopic ratios. We find only ten ratios that agree with their solar values, all the others differ by at least a factor of three from the solar values in all the considered sets of models. The γ -process within C–O shell mergers mostly influences the isotopic ratios that involve intermediate and heavy proton-rich isotopes with A  > 100. Conclusions. We conclude that there are large discrepancies both among the different data sets and between the model predictions and the solar abundance distribution. More calculations are needed; particularly updating the nuclear network, because the majority of the models considered in this work do not use the latest reaction rates for the γ -process nucleosynthesis. Moreover, the role of C–O shell mergers requires further investigation. 
    more » « less
  2. 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
  3. We studied the detectability and reconstruction of gravitational waves from core-collapse supernova multidimensional models using simulated data from detectors predicted to operate in the late 2020s and early 2030s. We found that the detection range will improve by a factor of around two with respect to the second-generation gravitational-wave detectors, and the sky localization will significantly improve. We analyzed the reconstruction accuracy for the lower frequency and higher frequency portion of supernova signals with a 250 Hz cutoff. Since the waveform’s peak frequencies are usually at high frequencies, the gravitational-wave signals in this frequency band were reconstructed more accurately. 
    more » « less
  4. Nuclear clustering affects the nucleosynthesis occurring in a number of astrophysical environments. Highly-clusterized nuclear states typically occur near particle thresholds and therefore can produce dramatic impacts on the nuclear reaction rates. This is especially true for astrophysical explosions that are driven by the consumption of helium as fuel. Such burning can occur in X-ray bursts, supernovae type Ia, and core-collapse supernovae for instance. This article will focus on the explosive astrophysical events in which nuclear clustering is most important, will discuss the types of information and tools necessary to estimate the astrophysical reaction rates, and will discuss example experiments at Notre Dame and other facilities that have or will be performed to measure the critical nuclear data needed for such estimates. 
    more » « less
  5. Abstract The simulation of heavy element nucleosynthesis requires input from yet-to-be-measured nuclear properties. The uncertainty in the values of these off-stability nuclear properties propagates to uncertainties in the predictions of elemental and isotopic abundances. However, for any given astrophysical explosion, there are many different trajectories, i.e., temperature and density histories, experienced by outflowing material, and thus different nuclear properties can come into play. We consider combined nucleosynthesis results from 460,000 trajectories from a black hole accretion disk and find the spread in elemental predictions due solely to unknown nuclear properties to be a factor of a few. We analyze this relative spread in model predictions due to nuclear variations and conclude that the uncertainties can be attributed to a combination of properties in a given region of the abundance pattern. We calculate a cross-correlation between mass changes and abundance changes to show how variations among the properties of participating nuclei may be explored. Our results provide further impetus for measurements of multiple quantities on individual short-lived neutron-rich isotopes at modern experimental facilities. 
    more » « less