Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explosive nucleosynthesis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparing Early Evolution of Flames in X-Ray Bursts in Two and Three Dimensions

We explore the early evolution of flame ignition and spreading on the surface of a neutron star in three dimensions, in the context of X-ray bursts. We look at the nucleosynthesis and morphology of the burning front and compare to two-dimensional axisymmetric simulations to gauge how important a full three-dimensional treatment of the flame is for the early dynamics. Finally, we discuss the progress toward full-star resolved flame simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Potential Impact of Fast Flavor Oscillations on Neutrino-driven Winds and Their Nucleosynthesis

The wind driven by the intense neutrino emission from a protoneutron star (PNS) is an important site for producing nuclei heavier than the Fe group. Because of certain features in the neutrino angular distributions, the so-called fast flavor oscillations may occur very close to the PNS surface, effectively resetting the neutrino luminosities and energy spectra that drive the wind. Using the unoscillated neutrino emission characteristics from two core-collapse supernova simulations representative of relevant progenitors at the lower and higher mass end, we study the potential effects of fast flavor oscillations on neutrino-driven winds and their nucleosynthesis. We find that such oscillations can increase the total mass loss by factors up to ~1.5–1.7 and lead to significantly more proton-rich conditions. The latter effect can greatly enhance the production of 64 Zn and the so-called light p-nuclei 74 Se, 78 Kr, and 84 Sr. Implications for abundances in metal-poor stars, Galactic chemical evolution in general, and isotopic anomalies in meteorites are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

New 26 P( p, γ ) 27 S Thermonuclear Reaction Rate and Its Astrophysical Implications in the rp -process

Accurate nuclear reaction rates for 26 P(p, γ) 27 S are pivotal for a comprehensive understanding of the rp-process nucleosynthesis path in the region of proton-rich sulfur and phosphorus isotopes. However, large uncertainties still exist in the current rate of 26 P(p, γ) 27 S because of the lack of nuclear mass and energy level structure information for 27 S. We reevaluate this reaction rate using the experimentally constrained 27 S mass, together with the shell model predicted level structure. It is found that the 26 P(p, γ) 27 S reaction rate is dominated by a direct capture reaction mechanism despite the presence of three resonances at E = 1.104, 1.597, and 1.777 MeV above the proton threshold in 27 S. The new rate is overall smaller than the other previous rates from the Hauser–Feshbach statistical model by at least 1 order of magnitude in the temperature range of X-ray burst interest. In addition, we consistently update the photodisintegration rate using the new 27 S mass. The influence of new rates of forward and reverse reaction in the abundances of isotopes produced in the rp-process is explored by postprocessing nucleosynthesis calculations. The final abundance ratio of 27 S/ 26 P obtained using the new rates is only 10% of that from the old rate. The abundance flow calculations show that the reaction path 26 P(p, γ) 27 S(β + ,ν) 27 P is not as important as previously thought for producing 27 P. The adoption of the new reaction rates for 26 P(p, γ) 27 S only reduces the final production of aluminum by 7.1% and has no discernible impact on the yield of other elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Impact of Hypernova νp-process Nucleosynthesis on the Galactic Chemical Evolution of Mo and Ru

We calculate the Galactic Chemical Evolution of Mo and Ru by taking into account the contribution from νp-process nucleosynthesis. We estimate yields of p-nuclei such as 92,94 Mo and 96,98 Ru through the νp-process in various supernova progenitors based upon recent models. In particular, the νp-process in energetic hypernovae produces a large amount of p-nuclei compared to the yield in ordinary core-collapse SNe. Because of this, the abundances of 92,94 Mo and 96,98 Ru in the Galaxy are significantly enhanced at [Fe/H] = 0 by the νp-process. We find that the νp-process in hypernovae is the main contributor to the elemental abundance of 92 Mo at low metallicity [Fe/H] < −2. Our theoretical prediction of the elemental abundances in metal-poor stars becomes more consistent with observational data when the νp-process in hypernovae is taken into account.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Kilonova Emissions from Neutron Star Merger Remnants: Implications for the Nuclear Equation of State

Multimessenger observations of binary neutron star mergers can provide valuable information on the nuclear equation of state (EOS). Here, we investigate the extent to which electromagnetic observations of the associated kilonovae allow us to place constraints on the EOS. For this, we use state-of-the-art three-dimensional general-relativistic magnetohydrodynamics simulations and detailed nucleosynthesis modeling to connect properties of observed light curves to properties of the accretion disk, and hence, the EOS. Using our general approach, we use multimessenger observations of GW170817/AT2017gfo to study the impact of various sources of uncertainty on inferences of the EOS. We constrain the radius of a 1.4M ⊙ neutron star to lie within 10.30 ≤ R 1.4 ≤ 13.0 km and the maximum mass to be M TOV ≤ 3.06M ⊙ .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Gamma Rays from Kilonova: A Potential Probe of r-process Nucleosynthesis

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of r-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. Futhermore, we show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Elemental and Isotopic Yields from T Coronae Borealis: Predictions and Uncertainties

T Coronae Borealis (T CrB) is a symbiotic recurrent nova system expected to undergo its next outburst within the next 2 yr. Recent hydrodynamic simulations have predicted the nucleosynthetic yields for both carbon–oxygen (CO) and oxygen–neon (ONe) white-dwarf models, but without accounting for thermonuclear reaction rate uncertainties. We perform detailed Monte Carlo postprocessing nucleosynthesis calculations based on updated thermonuclear reaction rates and uncertainties from the 2025 evaluation. We quantify the resulting abundance uncertainties and identify the key nuclear reactions that dominate them. Our results show that both the CO and ONe nova models robustly produce characteristic CNO isotopes. More pronounced abundance differences emerge for elements with A ≥ 20. Sulfur is the most robust observational discriminator between the CO and ONe nova models, with a model-to-model difference of a factor of ≈30 and minimal sensitivity to reaction rate uncertainties. Neon, silicon, and phosphorus exhibit even larger abundance differences (factors of ≈150–250), providing strong diagnostic potential. While their predicted yields are subject to larger uncertainties, these remain smaller than the model-to-model differences, allowing these elements to serve as useful, though less precise, tracers of white-dwarf composition. Chlorine, argon, and potassium also differ between models, but the 1σ-abundance ranges for the CO and ONe models overlap, reducing their present usefulness as composition tracers. We find that only nine nuclear reactions dominate the abundance uncertainties of the most diagnostically important isotopes, and their influence is largely independent of the underlying white-dwarf composition. These results provide guidance for future experimental efforts and for interpreting ejecta compositions in the next eruption of T CrB.

Chemical Abundances↗
DOE OSTI2022

Production of Very Light Elements and Strontium in the Early Ejecta of Neutron Star Mergers

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 −6 M ⊙ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while ∼1.5–11 × 10 −6 M ⊙ of helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r -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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Multidimensional Parameter Study of Double Detonation Type Ia Supernovae Originating from Thin Helium Shell White Dwarfs

Despite the importance of Type Ia supernovae (SNe Ia) throughout astronomy, the precise progenitor systems and explosion mechanisms that drive SNe Ia are still unknown. An explosion scenario that has gained traction recently is double detonation, in which an accreted shell of He detonates and triggers a secondary detonation in the underlying white dwarf. Our research presents a number of high-resolution, multidimensional, full-star simulations of thin-He-shell, sub-Chandrasekhar-mass white dwarf progenitors that undergo a double detonation. This suite of thin-shell progenitors incorporates He shells that are thinner than those in previous multidimensional studies. We confirm the viability of the double detonation across a range of He-shell parameter space, as well as present bulk yields and ejecta profiles for each progenitor. The yields obtained are generally consistent with previous works and indicate the likelihood of producing observables that resemble SNe Ia. The dimensionality of our simulations allow us to examine features of the double detonation more closely, including the details of the off-center secondary ignition and asymmetric ejecta. We find considerable differences in the high-velocity extent of postdetonation products across different lines of sight. The data from this work will be used to generate predicted observables and may further support the viability of the double detonation scenario as an SN Ia channel, as well as show how the properties of the progenitor or viewing angle may influence trends in observable characteristics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Resolution Study of Thermonuclear Initiation in White Dwarf Tidal Disruption Events

Abstract We study the initiation of thermonuclear detonations in tidally disrupted white dwarf stars by intermediate-mass (10 3 M ⊙ ) black holes. The length scales required to resolve the initiation mechanism are not easily reached in 3D, so instead we have devised 2D proxy models, which, together with a logarithmic gridding strategy, can adequately capture detonation wave fronts as material undergoes simultaneous compression and stretching from tidal forces. We consider 0.15 and 0.6 M ⊙ white dwarf stars parameterized by tidal strengths in the range β = 4–23. High spatial resolution elucidates the manner and conditions leading to thermonuclear detonation, linking the initiation sequence to stellar composition and tidal strength. All of our models suffer sustained detonations triggered by a combination of adiabatic compression, mild thermonuclear preconditioning, and collisional heating, in degrees depending primarily on tidal strength. We find that many diagnostics, such as temperature, total released energy, and iron-group products, are fairly well converged (better than 10%) at resolutions below 10 km along the scale height of the orbital plane. The exceptions are intermediate-mass transients like calcium, which remain uncertain up to factors of 2, even at 1 km resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Comprehensive Analysis of the Neutrino Process in Core-collapsing Supernovae

We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Comparison of the Core-collapse Evolution of Two Nearly Equal-mass Progenitors

Abstract We compare the core-collapse evolution of a pair of 15.8 M ☉ stars with significantly different internal structures, a consequence of the bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79 M ☉ progenitors have core masses (masses interior to an entropy of 4 k B baryon −1 ) of 1.47 and 1.78 M ☉ and compactness parameters ξ 1.75 of 0.302 and 0.604, respectively. The core-collapse simulations are carried out in 2D to nearly 3 s postbounce and show substantial differences in the times of shock revival and explosion energies. The 15.78 M ☉ model begins exploding promptly at 120 ms postbounce when a strong density decrement at the Si–Si/O shell interface, not present in the 15.79 M ☉ progenitor, encounters the stalled shock. The 15.79 M ☉ model takes 100 ms longer to explode but ultimately produces a more powerful explosion. Both the larger mass accretion rate and the more massive core of the 15.79 M ☉ model during the first 0.8 s postbounce time result in larger ν e / ν ¯ e luminosities and RMS energies along with a flatter and higher-density heating region. The more-energetic explosion of the 15.79 M ☉ model resulted in the ejection of twice as much 56 Ni. Most of the ejecta in both models are moderately proton rich, though counterintuitively the highest electron fraction ( Y e = 0.61) ejecta in either model are in the less-energetic 15.78 M ☉ model, while the lowest electron fraction ( Y e = 0.45) ejecta in either model are in the 15.79 M ☉ model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Sensitivity of He Flames in X-Ray Bursts to Nuclear Physics

Through the use of axisymmetric 2D hydrodynamic simulations, we further investigate laterally propagating flames in X-ray bursts (XRBs). Our aim is to understand the sensitivity of a propagating helium flame to different nuclear physics. Using the Castro simulation code, we confirm the phenomenon of enhanced energy generation shortly after a flame is established by adding 12 C(p, γ) 13 N(α, p) 16 O to the network, in agreement with the past literature. This sudden outburst of energy leads to a short accelerating phase, causing a drastic alteration in the overall dynamics of the flame in XRBs. Furthermore, we investigate the influence of different plasma screening routines on the propagation of the XRB flame. We finally examine the performance of simplified spectral deferred correction, a novel approach to hydrodynamics and reaction coupling incorporated in Castro, as an alternative to operator splitting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, ${\mathcal O }(100)$ nuclei network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics time step. We discuss the application of this method to convection in massive stars leading up to core collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU-accelerated.

Explosive nucleosynthesis↗
DOE OSTI2024

Thermonuclear 28 P(p, γ ) 29 S reaction rate and astrophysical implication in ONe nova explosion

An accurate 28 P(p, γ) 29 S reaction rate is crucial to defining the nucleosynthesis products of explosive hydrogen burning in ONe novae. Using the recently released nuclear mass of 29 S, together with a shell model and a direct capture calculation, we reanalyzed the 28 P(p, γ) 29 S thermonuclear reaction rate and its astrophysical implication. We focus on improving the astrophysical rate for 28 P(p, γ) 29 S based on the newest nuclear mass data. Our goal is to explore the impact of the new rate and associated uncertainties on the nova nucleosynthesis. We evaluated this reaction rate via the sum of the isolated resonance contribution instead of the previously used Hauser-Feshbach statistical model. The corresponding rate uncertainty at different energies was derived using a Monte Carlo method. Nova nucleosynthesis is computed with the 1D hydrodynamic code SHIVA. The contribution from the capture on the first excited state at 105.64 keV in 28 P is taken into account for the first time. We find that the capture rate on the first excited state in 28 P is up to more than 12 times larger than the ground-state capture rate in the temperature region of 2.5 × 10 7 K to 4 × 10 8 K, resulting in the total 28 P(p, γ) 29 S reaction rate being enhanced by a factor of up to 1.4 at ~1 × 10 9 K. In addition, the rate uncertainty has been quantified for the first time. It is found that the new rate is smaller than the previous statistical model rates, but it still agrees with them within uncertainties for nova temperatures. The statistical model appears to be roughly valid for the rate estimation of this reaction in the nova nucleosynthesis scenario. Using the 1D hydrodynamic code SHIVA, we performed the nucleosynthesis calculations in a nova explosion to investigate the impact of the new rates of 28 P(p, γ) 29 S. Our calculations show that the nova abundance pattern is only marginally affected if we use our new rates with respect to the same simulations but statistical model rates. Finally, the isotopes whose abundance is most influenced by the present 28 P(p, γ) 29 S uncertainty are 28 Si, 33,34 S, 35,37 Cl, and 36 Ar, with relative abundance changes at the level of only 3% to 4%.

Astronomy & Astrophysics↗
DOE OSTI2022

First direct measurement of the 13 $\mathrm{N(α, ρ)}$ 16 $\mathrm{O}$ reaction relevant for core-collapse supernovae nucleosynthesis

We report that understanding the explosion mechanism of a core-collapse supernova (CCSN) is important to accurately model CCSN scenarios for different progenitor stars using model-observation comparisons. The uncertainties of various nuclear reaction rates relevant for CCSN scenarios strongly affect the accuracy of these stellar models. Out of these reactions, the 13 N(α, p) 16 O reaction has been found to affect various stages of a CCSN at varying temperatures. This work presents the first direct measurement of the 13 N(α, p) 16 O reaction performed using a 34.6 MeV beam of radioactive 13 N ions and the active-target detector MUSIC (MUlti-Sampling Ionization Chamber) at Argonne National Laboratory. The resulting total 13 N(α, p) 16 O reaction cross sections from this measurement in the center-of-mass energy range of 3.26 - 6.02 MeV are presented and compared with calculations using the Hauser-Feshbach formalism. Uncertainties in the reaction rate have been dramatically reduced at CCSN temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗