An official website of the United States government
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
An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks
Abstract Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.
more »
« less
- Award ID(s):
- 1852143
- PAR ID:
- 10400590
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 936
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 6
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract A promising astrophysical site to produce the lighter heavy elements of the first r -process peak ( Z = 38 − 47) is the moderately neutron-rich (0.4 < Y e < 0.5) neutrino-driven ejecta of explosive environments, such as core-collapse supernovae and neutron star mergers, where the weak r -process operates. This nucleosynthesis exhibits uncertainties from the absence of experimental data from ( α , xn ) reactions on neutron-rich nuclei, which are currently based on statistical model estimates. In this work, we report on a new study of the nuclear reaction impact using a Monte Carlo approach and improved ( α , xn ) rates based on the Atomki-V2 α optical model potential. We compare our results with observations from an up-to-date list of metal-poor stars with [Fe/H] < −1.5 to find conditions of the neutrino-driven wind where the lighter heavy elements can be synthesized. We identified a list of ( α , xn ) reaction rates that affect key elemental ratios in different astrophysical conditions. Our study aims to motivate more nuclear physics experiments on ( α , xn ) reactions using the current and new generation of radioactive beam facilities and also more observational studies of metal-poor stars.more » « less
-
We demonstrate that, for systems with spin–orbit coupling and an odd number of electrons, the standard fewest switches surface hopping algorithm does not conserve the total linear or angular momentum. This lack of conservation arises not so much from the hopping direction (which is easily adjusted) but more generally from propagating adiabatic dynamics along surfaces that are not time reversible. We show that one solution to this problem is to run along eigenvalues of phase-space electronic Hamiltonians H(R, P) (i.e., electronic Hamiltonians that depend on both nuclear position and momentum) with an electronic–nuclear coupling Γ · P [see Eq. (25)], and we delineate the conditions that must be satisfied by the operator Γ. The present results should be extremely useful as far as developing new semiclassical approaches that can treat systems where the nuclear, electronic orbital, and electronic spin degrees of freedom altogether are all coupled together, hopefully including systems displaying the chiral-induced spin selectivity effect.more » « less
-
We explore the Quantum Chromodynamics (QCD) phase diagram's complexities, including quark deconfinement transitions, liquid-gas phase changes, and critical points, using the chiral mean-field (CMF) model that is able to capture all these features. We introduce a vector meson renormalization within the CMF framework, enabling precise adjustments of meson masses and coupling strengths related to vector meson interactions. Performing a new fit to the deconfinement potential, we are able to replicate recent lattice QCD results, low energy nuclear physics properties, neutron star observational data, and key phase diagram features as per modern constraints. This approach enhances our understanding of vector mesons' roles in mediating nuclear interactions and their impact on the equation of state, contributing to a more comprehensive understanding of the QCD phase diagram and its implications for nuclear and astrophysical phenomena.more » « less
-
We present a robust and efficient method for simulating Lagrangian solid-fluid coupling based on a new operator splitting strategy. We use variational formulations to approximate fluid properties and solid-fluid interactions, and introduce a unified two-way coupling formulation for SPH fluids and FEM solids using interior point barrier-based frictional contact. We split the resulting optimization problem into a fluid phase and a solid-coupling phase using a novel time-splitting approach with augmentedcontact proxies, and propose efficient custom linear solvers. Our technique accounts for fluids interaction with nonlinear hyperelastic objects of different geometries and codimensions, while maintaining an algorithmically guaranteed non-penetrating criterion. Comprehensive benchmarks and experiments demonstrate the efficacy of our method.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3847/1538-4357/ac8478
