Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear astrophysics”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Self-consistent microscopic calculations for electron captures on nuclei in core-collapse supernovae

Calculations for electron capture rates on nuclei with atomic numbers between 𝑍 = 20 and 𝑍 = 52 are performed in a self-consistent finite-temperature covariant energy density functional theory within the relativistic quasiparticle random-phase approximation. Electron captures on these nuclei contribute most to reducing the electron fraction during the collapse phase of core-collapse supernovae. The rates include contributions from allowed (Gamow-Teller) and first-forbidden (FF) transitions, and it is shown that the latter become dominant at high stellar densities and temperatures. Temperature-dependent effects such as Pauli unblocking and transitions from thermally excited states are also included. The new rates are implemented in a spherically symmetric one-dimensional simulation of the core-collapse phase. The results indicate that the increase in electron capture rates, due to inclusion of FF transitions, leads to reductions of the electron fraction at nuclear saturation density, the peak neutrino luminosity, and enclosed mass at core bounce. The new rates reaffirm that the most relevant nuclei for the deleptonization situate around the 𝑁 = 50 and 82 shell closures, but, compared to previous simulations, nuclei are less proton rich. Here, the new rates developed in this work are available, and will be of benefit to improve the accuracy of multidimensional supernova simulations.

Electron & muon capture↗
DOE OSTI2025

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗
DOE OSTI2025

Symmetry-Energy Dependence of the Bulk Viscosity of Nuclear Matter

We clarify how the weak-interaction-driven bulk viscosity ζ and the bulk relaxation time τ Π of neutrino-transparent n p e matter depend on the nuclear symmetry energy. We show that, at saturation density, the equation-of-state dependence of these transport quantities is fully determined by the experimentally constrained nuclear symmetry energy S and its slope L . Variations of L can change the bulk viscosity by orders of magnitude, which can affect both the dissipative and the conservative tidal response of neutron stars. Furthermore, this suggests that both conservative and dissipative effects encoded in the gravitational-wave signatures of binary neutron star inspirals may help constrain nuclear symmetry-energy properties.

Neutron stars & pulsars↗
DOE OSTI2025

Bayesian quantification of observability and equation of state of twin stars

The possibility of discovering twin stars, two neutron stars (NSs) with the same mass but different radii, is usually studied in forward modelings by using a restricted number of NS matter equations of state (EOSs) encapsulating a first-order phase transition from hadronic to quark matter (QM). Informing our likelihood function with the NS radius data from GW170817 and using a metamodel with nine parameters capable of mimicking most NS EOSs available in the literature, we conduct a Bayesian quantification of the observability and underlying EOSs of twin stars. Of the accepted EOSs, between 12 and 18% yield twin stars, depending on the restrictions we place on the second branch. The possibility of twin stars remains robust even under recent observational constraints. Here, we show that many of these twin star scenarios are observable with currently available levels of accuracy in measuring NS radii. We also present the marginalized posterior probability density functions (PDFs) of every EOS parameter for each of four mass-radius correlation topologies. We find that the inferred EOS depends sensitively on not only whether twin stars are present, but also the category of twin stars, indicating that the observation of twin stars would provide a strong constraint on the underlying EOS. In particular, for two coexisting hybrid stars having QM cores at different densities, the PDF for QM speed of sound squared 𝑐$^2_{qm}$ has two peaks, one below and another above the conformal limit 𝑐$^2_{qm}$ = 1/3 predicted by perturbative QCD.

QCD in nuclear reactions↗
DOE OSTI2025

Measurements of Gamow-Teller transitions from 59 Co via the 59 Co ⁢(𝑡, 3 He +𝛾) charge-exchange reaction and its application to the stellar electron-capture rates

Electron-capture reactions on iron-group nuclei play a crucial role in the late stages of massive star evolution. Since stellar evolution simulations depend on accurate electron-capture rates—which are highly sensitive to the detailed Gamow-Teller (GT) strength distributions—reliable theoretical models are essential. However, experimental data on GT strength distributions are scarce. High-resolution measurements are therefore vital for benchmarking and improving these theoretical calculations. To provide high-resolution data on Gamow-Teller strength distributions of iron-group nuclei and to compare these results with theoretical calculations within this mass region. Differential cross sections for the 59 Co ⁢(𝑡, 3 He)⁢ 59 Fe charge-exchange reaction at 115 MeV/u were measured using the S800 spectrometer. Furthermore, to resolve individual levels that are not distinguishable in the S800 particle singles data, coincident 𝛾 rays from the 59 Fe residual nucleus were detected by using the Gamma-Ray Energy Tracking In-beam Nuclear Array 𝛾-ray tracking array. Here, the Gamow-Teller transition strength distribution from the ground state of 59 Co to 59 Fe was extracted up to an excitation energy of 10 MeV. Additionally, transition strengths for several low-lying states were determined from coincident 𝛾-ray measurements. Electron-capture rates calculated using the present data indicate that these low-lying states contribute significantly to the overall rates in relevant stellar environments. The experimental results show reasonable agreement with theoretical predictions based on both shell-model and projected shell-model calculations. High-resolution data on Gamow-Teller strength distributions—particularly for individual low-lying states—are essential for accurately determining electron-capture rates in iron-group nuclei. Coincident 𝛾-ray measurements provide a powerful tool for obtaining such detailed information. While the present work demonstrates that shell-model calculations successfully reproduce the experimental results, such comparisons are scarce and more experimental data are desirable.

59 ≤ A ≤ 89↗
DOE OSTI2025

Crust Composition and the Shallow Heat Source in KS 1731–260

The presence of a strong shallow heat source of unknown origin in accreting neutron star crusts has been inferred by analyzing x-ray observations of their cooling in quiescence. Here, we model the cooling of KS 1731–260 using realistic crust compositions and nuclear heating and cooling sources from detailed nuclear reaction network calculations. We find that the required strength of the shallow heat source in KS 1731–260 is reduced by more than a factor of 3 compared to previous analysis, a 5⁢𝜎 difference that alleviates the need for exotic solutions. Our analysis also suggests the existence of an impure nuclear pasta layer in the inner crust of KS 1731–260, though future observations will provide more stringent constraints. In addition, we obtain constraints on the dominant surface burning modes of KS 1731–260 over its history.

neutron star cooling↗
DOE OSTI2025

New Class of Three-Nucleon Forces and Their Implications

Here, we identify a new class of three-nucleon forces that arises in the low-energy effective theory of nuclear interactions including pions. We estimate their contribution to the energy of neutron and nuclear matter and find that it can be as important as the leading-order three-nucleon forces previously considered in the literature. The magnitude of this force is set by the strength of the coupling of pions to two nucleons and is presently not well constrained by experiments. The implications for nuclei, nuclear matter, and the equation of state of neutron matter are briefly discussed.

chiral perturbation theory↗
DOE OSTI2026

Enhanced antineutrino emission from 𝛽 decay in core-collapse supernovae with self-consistent weak decay rates

Nuclear weak-interaction rates are known to exert a prominent effect in the late stages of stellar collapse. Despite their importance, most studies to date on core-collapse supernovae (CCSNe) have focused primarily on the effects of electron captures, neglecting 𝛽 decay contributions. In this work, we present the first CCSNe simulation incorporating global 𝛽 decay rates from a microscopic theory. These are enabled by a large-scale evaluation of both electron capture and 𝛽 decay rates, obtained self-consistently utilizing the relativistic energy density functional theory and finite-temperature quasiparticle random-phase approximation. Including 𝛽 decay leads to a dramatic enhancement of the pre-bounce antineutrino signal as the antineutrino emissivity increases by more than two orders of magnitude and the luminosity by a factor of 50 relative to thermal emission alone, while the average antineutrino energy increases by over 1 MeV. Furthermore, it is expected that these new rates could help us constrain the model uncertainties related to weak-interaction processes, improving the prediction of antineutrino signal during the final stages of stellar death.

Beta decay↗
DOE OSTI2026

Is the trace anomaly at its minimum value at neutron star centers?

While the equation of state (EOS) 𝑃⁡(𝜖) of neutron star (NS) matter has been extensively studied, the EOS-parameter 𝜙 =𝑃/𝜖 or equivalently the dimensionless trace anomaly Δ =1/3 −𝜙, which quantifies the balance between pressure 𝑃 and energy density 𝜖, remains far less explored, especially in NS cores. Its bounds and density profile carry crucial information about the nature of superdense matter. Physically, the EOS-parameter 𝜙 represents the mean stiffness of matter accumulated from the stellar surface up to a given density. Based on the intrinsic structure of the Tolman-Oppenheimer-Volkoff equations, we show that 𝜙 decreases monotonically outward from the NS center, independent of any specific input NS EOS model. Furthermore, observational evidence of a peak in the SSS density profile near the center effectively rules out a valley and a subsequent peak in the radial profile of 𝜙 at similar densities, reinforcing its monotonic decrease. Finally, these model-independent relations impose strong constraints on the near-center behavior of the EOS-parameter 𝜙, particularly demonstrating that the mean stiffness (or equivalently Δ) reaches a local maximum (minimum) at the center.

nuclear astrophysics↗
DOE OSTI2020

Stringent upper limit on the direct 3α decay of the Hoyle state in 12 C

We investigate an implication of the most recent observation of a second J π = 2 + state in 12 C, which was measured using the 12 C(γ,α) 8 Be (g.s.) reaction. In addition to the dissociation of 12 C to an α-particle and 8 Be in its ground state, a small fraction of events (2%) were identified as direct decays and decays to excited states in 8 Be. This allowed a limit on the direct 3α partial decay width to be determined as Γ 3α < 32(4) keV. Since this 2 + state is predicted by all theoretical models to be a collective excitation of the Hoyle state, the 3α partial width of the Hoyle state is calculable from the ratio of 3α decay penetrabilities of the Hoyle and 2 + states. Here, this was calculated by using the semiclassical Wenzel-Kramers-Brillouin approach and we deduce a stringent upper limit for the direct decay branching ratio of the Hoyle state of Γ 3α /Γ < 5.7 × 10 –6 , over an order of magnitude lower than previously reported. This result places the direct measurement of this rare decay mode beyond current experimental capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

$β$-delayed-neutron studies of 135,136 Sb and 140 I performed with trapped ions

In this work, Beta-delayed-neutron $(βn)$ spectroscopy was performed using the Beta-decay Paul Trap and an array of radiation detectors. The $βn$ branching ratios and energy spectra for 135,136 Sb and 140 I were obtained by measuring the time of flight of recoil ions emerging from the trapped ion cloud. These nuclei are located at the edge of an isotopic region identified as having $βn$ branching ratios that impact the $r$-process abundance pattern around the A ≈ 130 peak. For 135,136 Sb and 140 I, $βn$ branching ratios of 14.6(13)%, 17.6(29)%, and 7.6(28)% were determined, respectively. The $βn$ energy spectra obtained for 135 Sb and 140 I are compared with results from direct neutron measurements, and the $βn$ energy spectrum for 136 Sb has been measured for the first time.

42 ENGINEERING↗
DOE OSTI2020

Low-energy 23 Al β-delayed proton decay and 22 Na destruction in novae

The radionuclide 22 Na is a target of γ-ray astronomy searches, predicted to be produced during thermonuclear runaways driving classical novae. The 22 Na (p, γ) 23 Mg reaction is the main destruction channel of 22 Na during a nova, hence, its rate is needed to accurately predict the 22 Na yield. However, experimental determinations of the resonance strengths have led to inconsistent results. In this Rapid Communication, we report a measurement of the branching ratios of the 23 Al β-delayed protons as a probe of the key 204-keV (center-of-mass) 22 Na (p, γ) 23 Mg resonance strength. We report a factor of 5 lower branching ratio compared to the most recent literature value. The variation in 22 Na yield due to nuclear data inconsistencies was assessed using a series of hydrodynamic nova outburst simulations and has increased to a factor of 3.8, corresponding to a factor of ≈ 2 uncertainty in the maximum detectability distance. Finally, this is the first reported scientific measurement using the Gaseous Detector with Germanium Tagging system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Primary fission fragment mass yields across the chart of nuclides

In this work, we have calculated a complete set of primary fission fragment mass yields, Y(A), for heavy nuclei across the chart of nuclides, including those of particular relevance to the rapid neutron capture process (r process) of nucleosynthesis. We assume that the nuclear shape dynamics are strongly damped, which allows for a description of the fission process via Brownian shape motion across nuclear potential-energy surfaces. The macroscopic energy of the potential was obtained with the Finite-Range Liquid-Drop Model (FRLDM), while the microscopic terms were extracted from the single-particle level spectra in the fissioning system by the Strutinsky procedure for the shell energies and the BCS treatment for the pairing energies. For each nucleus considered, the fission fragment mass yield, Y(A), is obtained from 50 000 to 500 000 random walks on the appropriate potential-energy surface. The full mass and charge yield, Y(Z,A), is then calculated by invoking the Wahl systematics. With this method, we have calculated a comprehensive set of fission-fragment yields from over 3800 nuclides bounded by 80 ≤ Z ≤ 130 and A ≤ 330; these yields are provided as an ASCII formatted database in the Supplemental Material. We compare our yields to known data and discuss general trends that emerge in low-energy fission yields across the chart of nuclides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Enhanced dynamics in fusion of neutron-rich oxygen nuclei at above-barrier energies

Above-barrier fusion cross sections for an isotopic chain of oxygen isotopes with A = 16 –19 incident on a C 12 target are presented. Experimental data are compared with both static and dynamical microscopic calculations. Here, these calculations are unable to explain the ≈ 37 % increase in the average above-barrier fusion cross section observed for O 19 as compared to β -stable oxygen isotopes. This result suggests that for neutron-rich nuclei existing time-dependent Hartree-Fock calculations underpredict the role of dynamics at near-barrier energies. High-quality measurement of above-barrier fusion for an isotopic chain of increasingly neutron-rich nuclei offers an effective means to probe this fusion dynamics.

6 ≤ A ≤ 19↗
DOE OSTI2020

Consistency of $^{16}O(n,α)$ cross sections

The evaluated $^{16}O(n,α)$ cross sections in the ENDF/B-VIII.0 nuclear data library remain uncertain because of systematic discrepancies in the measured data. In the energy region below the first excited state, R-matrix analyses rely heavily on the inverse reaction, and in particular, they rely on the measured $^{13}C(α,n)^{16}O$ cross sections reported by Bair et al. in 1973 and Harissopulos et al. in 2005. The Harissopulos cross section values are systematically lower than those previously reported by Bair. Here, drawing on the available experimental information, this paper briefly describes and demonstrates that the two sets of measured cross sections are consistent.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Fast neutrino cooling of nuclear pasta in neutron stars: Molecular dynamics simulations

The direct Urca process of rapid neutrino emission can occur in nonuniform nuclear pasta phases that are expected in the inner crusts of neutron stars. Here, the periodic potential for a nucleon in nuclear pasta allows momentum conservation to be satisfied for direct Urca reactions. We improve on earlier work by modeling a rich variety of pasta phases (gnocchi, waffle, lasagna, and anti-spaghetti) with large-scale molecular dynamics simulations. We find that the neutrino luminosity due to direct Urca reactions in nuclear pasta can be 3 to 4 orders of magnitude larger than that from the modified Urca process in the NS core. Thus neutrino radiation from pasta could dominate radiation from the core and this could significantly impact the cooling of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Fission and the r -process nucleosynthesis of translead nuclei in neutron star mergers

In this report we study the impact of fission on the production and destruction of translead nuclei during the r-process nucleosynthesis occurring in neutron-star mergers. Abundance patterns and rates of nuclear energy production are obtained for different ejecta conditions using three sets of stellar reaction rates, one of which is based on microscopic and consistent calculations of nuclear masses, fission barriers, and collective inertias. We show that the accumulation of fissioning material during the r process can strongly affect the free neutron abundance after the r-process freeze-out. This leads to a significant impact on the abundances of heavy nuclei that undergo α decay or spontaneous fission, affecting the radioactive energy production by the ejecta at timescales relevant for kilonova emission.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Quantifying uncertainties and correlations in the nuclear-matter equation of state

We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory (χ EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χ EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. Furthermore, we then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium χ EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗