Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nuclear densities”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Chiral Effective Field Theory and the High-Density Nuclear Equation of State

Born in the aftermath of core-collapse supernovae, neutron stars contain matter under extraordinary conditions of density and temperature that are difficult to reproduce in the laboratory. In recent years, neutron star observations have begun to yield novel insights into the nature of strongly interacting matter in the high-density regime where current theoretical models are challenged. At the same time, chiral effective field theory has developed into a powerful framework to study nuclear matter properties with quantified uncertainties in the moderate-density regime for modeling neutron stars. In this article, we review recent developments in chiral effective field theory and focus on many-body perturbation theory as a computationally efficient tool for calculating the properties of hot and dense nuclear matter. We also demonstrate how effective field theory enables statistically meaningful comparisons among nuclear theory predictions, nuclear experiments, and observational constraints on the nuclear equation of state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Finite-range pairing in nuclear density functional theory

Pairing correlations are ubiquitous in low-energy states of atomic nuclei. To incorporate them within nuclear density functional theory, often used for global computations of nuclear properties, pairing functionals that generate nucleonic pair densities and pairing fields are introduced. Many pairing functionals currently used can be traced back to zero-range nucleon-nucleon interactions. Unfortunately, such functionals are plagued by deficiencies that become apparent in large model spaces that contain unbound single-particle (continuum) states. In particular, the underlying computational schemes diverge as the single-particle space increases, and the results depend on how marginally occupied states are incorporated. These problems become more pronounced for pairing functionals that contain gradient-density dependence, such as in the Fayans functional. To remedy this, finite-range pairing functionals are introduced. In this study, this is done by folding the pair density with Gaussians. Here, we show that a folding radius of about 1 fm offers the best compromise between quality and stability, and substantially reduces the pathological behavior in different numerical applications.

Nuclear density functional theory↗
DOE OSTI2022

Building surrogate models of nuclear density functional theory with Gaussian processes and autoencoders

From the lightest Hydrogen isotopes up to the recently synthesized Oganesson (Z = 118), it is estimated that as many as about 8,000 atomic nuclei could exist in nature. Most of these nuclei are too short-lived to be occurring on Earth, but they play an essential role in astrophysical events such as supernova explosions or neutron star mergers that are presumed to be at the origin of most heavy elements in the Universe. Understanding the structure, reactions, and decays of nuclei across the entire chart of nuclides is an enormous challenge because of the experimental difficulties in measuring properties of interest in such fleeting objects and the theoretical and computational issues of simulating strongly-interacting quantum many-body systems. Nuclear density functional theory (DFT) is a fully microscopic theoretical framework which has the potential of providing such a quantitatively accurate description of nuclear properties for every nucleus in the chart of nuclides. Thanks to high-performance computing facilities, it has already been successfully applied to predict nuclear masses, global patterns of radioactive decay like β or γ decay, and several aspects of the nuclear fission process such as, e.g., spontaneous fission half-lives. Yet, predictive simulations of nuclear spectroscopy—the low-lying excited states and transitions between them—or of nuclear fission, or the quantification of theoretical uncertainties and their propagation to basic or applied nuclear science applications, would require several orders of magnitude more calculations than currently possible. However, most of this computational effort would be spent into generating a suitable basis of DFT wavefunctions. Such a task could potentially be considerably accelerated by borrowing tools from the field of machine learning and artificial intelligence. In this paper, we review different approaches to applying supervised and unsupervised learning techniques to nuclear DFT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Probing the high-density nuclear symmetry energy with the Ξ – / Ξ 0 ratio in heavy-ion collisions at s N N ≈ 3 GeV

Recent beam energy scan experiments at the BNL Relativistic Heavy Ion Collider by the STAR Collaboration found that hadronic interactions dominate the collective flow and the proton cumulant ratios are driven by baryon number conservation in a region of high baryon density in √ s NN = 3 GeV Au + Au reactions, indicating that the dense medium formed in such collisions is likely hadronic matter. Within an updated a relativistic transport model with momentum dependent isoscalar and isovector single-nucleon mean-field potentials corresponding to different symmetry energies at suprasaturation densities, the n/p, π – /π + , $K$ $^{0}_{s}$/K + , Σ – /Σ + , and Ξ – /Ξ 0 ratios are studied for central Au + Au collisions at √ s NN = 3 GeV, where the maximum central density reaches about (3.6–4.0) ρ 0 . The doubly strange Ξ – /Ξ 0 ratio is found to have the strongest sensitivity to the variation of high-density nuclear symmetry energy. Furthermore, the Ξ – /Ξ 0 ratio in relativistic heavy-ion reactions at √ s NN~3 GeV may help probe sensitively the poorly known symmetry energy of dense neutron-rich matter critically important for understanding various properties of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Imprints of High-Density Nuclear Symmetry Energy on Crustal Fraction of Neutron Star Moment of Inertia

The density dependence of nuclear symmetry energy E sym (ρ) remains the most uncertain aspect of the equation of state (EOS) of supradense neutron-rich nucleonic matter. Utilizing an isospin-dependent parameterization of the nuclear EOS, we investigate the implications of the observational crustal fraction of the neutron star (NS) moment of inertia ΔI/I for the E sym (ρ). We find that symmetry energy parameters significantly influence the ΔI/I, while the EOS of symmetric nuclear matter has a negligible effect. In particular, an increase in the slope L and skewness J sym of symmetry energy results in a larger ΔI/I, whereas an increase in the curvature K sym leads to a reduction in ΔI/I. Moreover, the ΔI/I is shown to have the potential for setting a lower limit of symmetry energy at densities exceeding 3 ρ 0 , particularly when L is constrained to values less than 60 MeV, thereby enhancing our understanding of supradense NS matter.

equation of state↗
DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗
DOE OSTI2023

Extracting model-independent nuclear level densities away from stability

The nuclear level density (NLD) is a fundamental measure of the complex structure of atomic nuclei at relatively high energies. Here, in this study, we present the first model-independent measurement of the absolute partial NLD for a short-lived nucleus. For this purpose we adapt the recently introduced “shape method” for β-decay experiments, providing the shape of the γ-ray strength function for exotic nuclei. In this work, we show that combining the shape method with the β-Oslo technique allows for the extraction of the NLD of the populated states without the need for theoretical input. This development opens the way for the extraction of experimental NLDs far from stability with major implications in astrophysical and other applications. We benchmark our approach using data for the stable 76 Ge nucleus, finding excellent agreement with previous experimental results. In addition, we present new experimental data and determine the absolute partial level density for the short-lived 88 Kr nucleus. Our results suggest a fivefold increase in the NLD for the case of 88 Kr, compared to the recommended values from semimicroscopic Hartree-Fock Bogoliubov calculations recommended by the RIPL3 nuclear data library. However, our results are in good agreement with other semimicroscopic level density models. We demonstrate the impact of our method on the 87 Kr(n, γ) neutron capture rate and show that our experimental uncertainties for NLDs fulfill the requirements needed for astrophysical calculations predicting r-process abundances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A new database website for nuclear level densities

We introduce a new open-access, web-based database (http://nld.ascsn.net), Current Archive of Nuclear Density of Levels (CANDL), that hosts experimental nuclear level density (NLD) datasets from a variety of techniques and energy ranges. Built using the Dash framework in Python, the database is designed to be interactive and user-friendly, allowing researchers to search, visualize, fit, and export NLD data with minimal effort. This resource includes data extracted from evaporation spectra, Oslo method variants, and other experimental techniques that cover excitation energies beyond the neutron resonance region. The database supports on-the-fly fitting with two widely-used phenomenological models—the Constant Temperature (CT) model and the Back-Shifted Fermi Gas (BSFG) model—selected for their simplicity and computational efficiency. Future versions aim to include additional datasets and model types, as well as easy-to-use interfaces to data science techniques. Here, this platform offers a vital tool for the nuclear physics, astrophysics, medicine, and reactor design communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Auxiliary function approach for determining symmetry energy at suprasaturation densities

Nuclear symmetry energy E sym (ρ) at density ρ is normally expanded or simply parameterized as a function of χ = (ρ – ρ 0 )/3ρ 0 in the form of Esym(ρ) ≈ S + Lχ + 2 –1 K sym χ 2 + 6 –1 J sym χ 3 + … using its magnitude S, slope L, curvature K sym and skewness J sym at the saturation density ρ0 of nuclear matter. Much progress has been made in recent years in constraining especially the S and L parameters using various terrestrial experiments and astrophysical observations. Furthermore, such kind of expansions/parameterizations do not converge at supra-saturation densities where χ is not small enough, hindering an accurate determination of high-density E sym (ρ) even if its characteristic parameters at ρ0 are all well determined by experiments/observations. By expanding the E sym (ρ) in terms of a properly chosen auxiliary function Π sym (χ,Θ sym ) with a parameter Θsym fixed accurately by an experimental E sym (ρ r ) value at a reference density ρr, we show that the shortcomings of the χ-expansion can be completely removed or significantly reduced in determining the high-density behavior of E sym (ρ). In particular, using two significantly different auxiliary functions, we show that the new approach effectively incorporates higher χ-order contributions and converges to the same E sym (ρ) much faster than the conventional χ-expansion at densities ≲3ρ 0 . Several quantitative demonstrations using Monte Carlo simulations are given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Optimized nuclear energy density functionals including long-range pion contributions

Nuclear energy density functionals successfully reproduce properties of nuclei across almost the entire nuclear chart. However, nearly all available functionals are phenomenological in nature and lack a rigorous connection to systematically improvable nuclear forces. This issue might be solved with an energy density functional obtained from first principles. As an intermediate step towards this goal we construct the GUDE family of functionals that is obtained from a hybrid scheme consisting of long-range pion-exchange contributions derived from chiral effective field theory at the Hartree-Fock level and a phenomenological Skyrme part. When including pion contributions beyond next-to-leading order in the chiral expansion, we find significant improvements over a reference Skyrme functional constructed following the same protocol. We analyze the importance of different pion contributions and identify which terms drive the observed improvements. Since pions are incorporated without adding further optimization parameters to the functionals, the improvements can be attributed to the functional form of these terms. Our work therefore suggests that the considered chiral contributions constitute useful ingredients for true energy density functionals. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Photon strength functions and nuclear level densities: invaluable input for nucleosynthesis

The pivotal role of nuclear physics in nucleosynthesis processes is being investigated, in particular the intricate influence of photon strength functions (PSFs) and nuclear level densities (NLDs) on shaping the outcomes of the i-, r- and p-processes. Exploring diverse NLD and PSF model combinations uncovers large uncertainties for (p, γ ), (n, γ ) and ( α , γ ) rates across many regions of the nuclear chart. These lead to potentially significant abundance variations of the nucleosynthesis processes and highlight the importance of accurate experimental nuclear data. Theoretical insights and advanced experimental techniques lay the ground work for profound understanding that can be gained of nucleosynthesis mechanisms and the origin of the elements. Recent results further underscore the effect of PSF and NLD data and its contribution to understanding abundance distributions and refining knowledge of the intricate nucleosynthesis processes. This article is part of the theme issue ‘The liminal position of Nuclear Physics: from hadrons to neutron stars’.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Bayesian inference of nuclear-matter density from proton scattering

Background: Proton elastic scattering at intermediate energy is widely employed as a tool for determining the matter radius of atomic nuclei. Here, the sensitivity of the approach relies on high-resolution measurements at small scattering angles and low-momentum transfer. Under these conditions, the Glauber multiple scattering theory accurately describes the proton-nucleus elastic cross section. Purpose: Investigate the sensitivity of the Glauber multiple scattering theory to uncertainties associated with input parameters such as the nuclear-matter density distribution and nucleon-nucleon data. Method: A joint Bayesian inference was performed using 12 angular distributions of elastic scattering at different energies on 58 Ni, 90 Zr, and 208 Pb targets. A Metropolis-Hastings algorithm was implemented to make an uncertainty quantification analysis for the input parameters used in the Glauber multiple scattering theory. Results: The experimental cross sections were fitted simultaneously using a joint Bayesian inference approach. Posterior probability density distributions of 42 input parameters were obtained from the analysis. A moderate correlation between the nuclear density parameters and the nucleon-nucleon cross sections was found. This correlation impacts the extraction of the nuclear-matter radius. Conclusions: The present analysis provided a consistent method for extracting the nuclear-matter density distribution of 58 Ni, 90 Zr, and 208 Pb from data across different incident energies. Due to the correlation of the nucleon-nucleon cross sections with the other input parameters, a constrained Bayesian inference using free nucleon-nucleon cross section data was performed. The nuclear-matter radii obtained from the analysis are in good agreement with multiple results reported in the literature.

190 ≤ A ≤ 219↗
DOE OSTI2025

HFBTHO-AD: Differentiation of a nuclear energy density functional code

The HFBTHO code implements a nuclear energy density functional solver to model the structure of atomic nuclei. HFBTHO has previously been used to calibrate energy functionals and perform sensitivity analysis by using derivative-free methods. To enable derivative-based optimization and uncertainty quantification approaches, we must compute the derivatives of HFBTHO outputs with respect to the parameters of the energy functional, which are a subset of all input parameters of the code. Here, we use the algorithmic/automatic differentiation (AD) tool Tapenade to differentiate HFBTHO. We compare the derivatives obtained using AD against finite-difference approximation and examine the performance of the derivative computation.

Algorithmic differentiation↗
DOE OSTI2024

Nuclear level density and γ -decay strength of Sr 93

This work presents the first experimentally determined nuclear level density and γ-ray strength function of the short-lived fission product 93 Sr, accomplished using the β-Oslo method. Direct measurement of the 92 Sr(n, γ) 93 Sr cross section is not currently possible, as the half-life of 2.66 hours is too short; instead, 93 Sr was formed through β decay of 93 Rb to excitation energies around the neutron separation energy. The γ-ray spectra were measured using a total absorption spectrometer at the National Superconducting Cyclotron Laboratory (NSCL) at Michigan State University (MSU). The statistical properties of the 93 Sr nucleus were experimentally determined, including the γ-ray strength function and nuclear level density. At low energies, the γ-ray strength function exhibits a constant γ-decay strength, rather than a slightly increasing strength with decreasing γ-ray energy as had been previously observed for several nuclei in this mid-mass region. Finally, these statistical properties were then implemented in the reaction code TALYS1.95 to calculate the 92 Sr(n, γ) 93 Sr cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Nuclear energy density functionals grounded in ab initio calculations

Here, we discuss the construction of a nuclear energy density functional (EDF) from ab initio computations and advocate the need for a methodical approach that is free from ad hoc assumptions. The equations of state (EoSs) of symmetric nuclear and pure neutron matter are computed using the chiral NNLO sat and the phenomenological AV4' + UIX c Hamiltonians as inputs to self-consistent Green's function (SCGF) and auxiliary field diffusion Monte Carlo (AFDMC) methods. We propose a convenient parametrization of the EoS as a function of the Fermi momentum and fit it on the SCGF and AFDMC calculations. We apply the ab initio based EDF to carry out an analysis of the binding energies and charge radii of different nuclei in the local density approximation. The NNLO sat -based EDF produces encouraging results, whereas the AV4' + UIX c -based one is farther from experiment. Possible explanations of these different behaviors are suggested, and the importance of gradient and spin-orbit terms is analyzed. Our paper paves the way for a practical and systematic way to merge ab initio nuclear theory and density functional theory, while shedding light on some critical aspects of this procedure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗