Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear density functional theory”

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 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

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 OSTI2020

Covariant Density Functional Theory in Nuclear Physics and Astrophysics

How does subatomic matter organize itself? Neutron stars are cosmic laboratories uniquely poised to answer this fundamental question that lies at the heart of nuclear science. Newly commissioned rare isotope facilities, telescopes operating across the entire electromagnetic spectrum, and ever more sensitive gravitational wave detectors will probe the properties of neutron-rich matter with unprecedented precision over an enormous range of densities. A coordinated effort between observation, experiment, and theoretical research is of paramount importance for realizing the full potential of these investments. Theoretical nuclear physics provides valuable insights into the properties of neutron-rich matter in regimes that are not presently accessible to experiment or observation. In particular, nuclear density functional theory is likely the only tractable framework that can bridge the entire nuclear landscape by connecting finite nuclei to neutron stars. This compelling connection is the main scope of the present review.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Fission of 240 Pu with Symmetry-Restored Density Functional Theory

Nuclear fission plays an important role in fundamental and applied science, from astrophysics to nuclear engineering, yet it remains a major challenge to nuclear theory. Theoretical methods used so far to compute fission observables rely on symmetry-breaking schemes where basic information on the number of particles, angular momentum, and parity of the fissioning nucleus is lost. In this letter, we analyze the impact of restoring broken symmetries in the benchmark case of 240 Pu.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Beyond the proton drip line: Bayesian analysis of proton-emitting nuclei

The limits of the nuclear landscape are determined by nuclear binding energies. Beyond the proton drip lines, where the separation energy becomes negative, there is not enough binding energy to prevent protons from escaping the nucleus. Predicting properties of unstable nuclear states in the vast territory of proton emitters poses an appreciable challenge for nuclear theory as it often involves far extrapolations. In addition, significant discrepancies between nuclear models in the proton-rich territory call for quantified predictions. With the help of Bayesian methodology, we mix a family of nuclear mass models corrected with statistical emulators trained on the experimental mass measurements, in the proton-rich region of the nuclear chart. Separation energies were computed within nuclear density functional theory using several Skyrme and Gogny energy density functionals. We also considered mass predictions based on two models used in astrophysical studies. Quantified predictions were obtained for each model using Bayesian Gaussian processes trained on separation-energy residuals and combined via Bayesian model averaging. Here, we obtained a good agreement between averaged predictions of statistically corrected models and experiment. In particular, we quantified model results for one- and two-proton separation energies and derived probabilities of proton emission. This information enabled us to produce a quantified landscape of proton-rich nuclei. The most promising candidates for two-proton decay studies have been identified. Finally, the methodology used in this work has broad applications to model-based extrapolations of various nuclear observables. It also provides a reliable uncertainty quantification of theoretical predictions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Quantified limits of the nuclear landscape

The chart of the nuclides is limited by particle drip lines beyond which nuclear stability to proton or neutron emission is lost. Predicting the range of particle-bound isotopes poses an appreciable challenge for nuclear theory as it involves extreme extrapolations of nuclear masses well beyond the regions where experimental information is available. Still, quantified extrapolations are crucial for a wide variety of applications, including the modeling of stellar nucleosynthesis. We use microscopic nuclear global mass models, current mass data, and Bayesian methodology to provide quantified predictions of proton and neutron separation energies as well as Bayesian probabilities of existence throughout the nuclear landscape all the way to the particle drip lines. Here, we apply nuclear density functional theory with several energy density functionals. We also consider two global mass models often used in astrophysical nucleosynthesis simulations. To account for uncertainties, Bayesian Gaussian processes are trained on the separation-energy residuals for each individual model, and the resulting predictions are combined via Bayesian model averaging. This framework allows to account for systematic and statistical uncertainties and propagate them to extrapolative predictions. We establish and characterize the drip-line regions where the probability that the nucleus is particle- bound decreases from 1 to 0. In these regions, we provide quantified predictions for one- and two-nucleon separation energies. According to our Bayesian model averaging analysis, 7759 nuclei with Z ≤ 119 have a probability of existence ≥ 0.5. The extrapolation results obtained in this study will be put through stringent tests when new experimental information on existence and masses of exotic nuclei becomes available. In this respect, the quantified landscape of nuclear existence obtained in this study should be viewed as a dynamical prediction that will be fine-tuned when new experimental information and improved global mass models become available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Variations in the charge radii of indium isotopes between 𝑁=52 and 82

Measurements of the 5⁢𝑝 2 𝑃 3/2 → 9⁢𝑠 2 𝑆 1/2 and 5⁢𝑝 2 𝑃 1/2 → 8⁢𝑠 2 𝑆 1/2 transitions in the indium atom, combined with new atomic physics calculations, were used to extract the changes in mean-square nuclear charge radii, 𝛿⁡⟨𝑟 2 ⟩, of the indium (𝑍 = 49) isotopes 101–111 In , 113–123 In , and 125–131 In . With a proton hole in the closed nuclear shell of 𝑍 = 50, indium provides a detailed study of the effect of unpaired nucleons adjacent to the proton-shell closure, allowing investigation into the charge radii for isotopes between the two major neutron-shell closures at 𝑁 = 50 and 𝑁 = 82. A study of the variations in charge radii between neighboring isotopes with neutron number (the ‘odd-even staggering') is presented and provides further insight and challenges for the theoretical description of the size of proton-hole nuclei. Two nuclear theories, density functional theory and the valence-space in-medium similarity renormalization group method, were employed to interpret the data. The new information obtained in this work provides valuable insights into the successes and shortcomings of the theoretical approaches employed.

ab initio calculations↗
DOE OSTI2020

Landscape of pear-shaped even-even nuclei

The phenomenon of reflection-asymmetric nuclear shapes is relevant to nuclear stability, nuclear spectroscopy, nuclear decays and fission, and the search for new physics beyond the standard model. Global surveys of ground-state octupole deformation, performed with a limited number of models, suggest that the number of pear-shaped isotopes is fairly limited across the nuclear landscape. Here, we carry out a global analysis of ground-state octupole deformations for particle-bound even-even nuclei with Z ≤ 110 and N ≤ 210 using nuclear density functional theory (DFT) with several nonrelativistic and covariant energy density functionals. In this way, we can identify the best candidates for reflection-asymmetric shapes. The calculations are performed in the frameworks of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory and relativistic Hartree-Bogoliubov theory using DFT solvers employing harmonic oscillator basis expansion. We consider five Skyrme and four covariant energy density functionals. We predict several regions of ground-state octupole deformation. In addition to the “traditional” regions of neutron-deficient actinide nuclei around 224 Ra and neutron-rich lanthanides around 146 Ba, we identified vast regions of reflection-asymmetric shapes in very neutron-rich nuclei around 200 Gd and 288 Pu, as well as in several nuclei around 112 Ba. Our analysis suggests several promising candidates with stable ground-state octupole deformation, primarily in the neutron-deficient actinide region, that can be reached experimentally. Detailed comparison between Skyrme and covariant models is performed. Octupole shapes predicted in this study are consistent with the current experimental information. This work can serve as the starting point of a systematic search for parity doublets in odd-mass and odd-odd nuclei, which will be of interest in the context of new physics searches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multimodal fission from self-consistent calculations

When multiple fission modes coexist in a given nucleus, distinct fragment yield distributions appear. Multimodal fission has been observed in a number of fissioning nuclei spanning the nuclear chart, and this phenomenon is expected to affect the nuclear abundances synthesized during the rapid neutron-capture process (𝑟-process). In this study, we generalize the previously proposed hybrid model for fission-fragment yield distributions to predict competing fission modes and estimate the resulting yield distributions. Here, our framework allows for a comprehensive large-scale calculation of fission-fragment yields suited for 𝑟-process nuclear network studies. Nuclear density functional theory is employed to obtain the potential energy and collective inertia tensor on a multidimensional collective space defined by mass multipole moments. Fission pathways and their relative probabilities are determined using the nudged elastic band method. Based on this information, mass and charge fission yields are predicted using the recently developed hybrid model. Fission properties of fermium isotopes are calculated in the axial quadrupole-octupole collective space for three energy density functionals (EDFs). Disagreement between the EDFs appears when multiple fission modes are present. Within our framework, the UNEDF⁢1 HFB EDF agrees best with experimental data. Calculations in the axial quadrupole-octupole-hexadecapole collective space improve the agreement with the experiment for SkM*. We also discuss the sensitivity of fission predictions on the choice of EDF for several superheavy nuclei. Fission-fragment yield predictions for nuclei with multiple fission modes are sensitive to the underlying EDF. For large-scale calculations in which a minimal number of collective coordinates is considered, UNEDF⁢1 HFB provides the best description of experimental data, though the sensitivity motivates robust quantification of the uncertainties of the theoretical model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Universal trend of charge radii of even-even Ca–Zn nuclei

Radii of nuclear charge distributions carry information about the strong and electromagnetic forces acting inside the atomic nucleus. Whereas the global behavior of nuclear charge radii is governed by the bulk properties of nuclear matter, their local trends are affected by quantum motion of proton and neutron nuclear constituents. The measured differential charge radii δ$\langle$$r^{2}_{c}$$\rangle$ between neutron numbers N = 28 and N = 40 exhibit a universal pattern as a function of n = N – 28 that is independent of the atomic number. In this work, we analyze this remarkable behavior in even-even nuclei from calcium to zinc using two state-of-the-art theories based on quantified nuclear interactions: the ab initio coupled cluster theory and nuclear density functional theory. Both theories reproduce the smooth rise of differential charge radii and their weak dependence on the atomic number. By considering a large set of isotopic chains, we show that this trend can be captured by just two parameters: the slope and curvature of δ$\langle$$r^{2}_{c}$$\rangle$(n). We demonstrate that these parameters show appreciable model dependence, and the statistical analysis indicates that they are not correlated with any single model property, i.e., they are impacted by both bulk nuclear properties as well as shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Microscopic calculation of fission product yields with particle-number projection

Fission fragments' charge and mass distribution is an important input to applications ranging from basic science to energy production or nuclear nonproliferation. In simulations of nucleosynthesis or calculations of superheavy elements, these quantities must be computed from models, as they are needed in nuclei where no experimental information is available. Until now, standard techniques to estimate these distributions were not capable of accounting for fine-structure effects, such as the odd-even staggering of the charge distributions. In this work, we combine a fully microscopic collective model of fission dynamics with a recent extension of the particle number projection formalism to provide the highest-fidelity prediction of the primary fission fragment distributions for the neutron-induced fission of 235 U and 239 Pu. Here, we show that particle-number projection is an essential ingredient to reproduce odd-even staggering in the charge yields and benchmark the performance of various empirical probability laws that could simulate its effect. This new approach also enables for the first time the realistic determination of two-dimensional isotopic yields within nuclear density functional theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Nuclear DFT analysis of electromagnetic moments in odd near doubly magic nuclei

We use the nuclear density functional theory to determine nuclear electric quadrupole and magnetic dipole moments in all one-particle and one-hole neighbours of eight doubly magic nuclei. We align angular momenta along the intrinsic axial-symmetry axis with broken time-reversal symmetry, which allows us to explore fully the self-consistent charge, spin, and current polarisation. Spectroscopic moments are determined for symmetry-restored wave functions and compared with available experimental data. We find that the obtained polarisations do not call for using quadrupole- or dipole-moment operators with effective charges or effective g-factors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Neural network emulation of spontaneous fission

Large-scale computations of fission properties are an important ingredient for nuclear reaction network calculations simulating rapid neutron-capture process (the 𝑟 process) nucleosynthesis. Due to the large number of fissioning nuclei potentially contributing to the 𝑟 process, a microscopic description of fission based on nuclear density functional theory (DFT) is computationally challenging. Here, we explore the use of neural networks (NNs) to construct DFT emulators capable of predicting potential energy surfaces and collective inertia tensors across the whole nuclear chart, starting from a minimal set of DFT calculations. We use constrained Hartree-Fock-Bogoliubov (HFB) calculations to predict the potential energy and collective inertia tensor in the axial quadrupole and octupole collective coordinates, for a set of nuclei in the 𝑟-process region. We then employ NNs to emulate the HFB energy and collective inertia tensor across the considered region of the nuclear chart. Least-action pathways characterizing spontaneous fission half-lives and fragment yields are then obtained by means of the nudged elastic band method. The potential energy predicted by NNs agrees with the DFT value to within a root-mean-square error of 500 keV, and the collective inertia components agree to within an order of magnitude. These results are largely independent of the NN architecture. The exit points on the outer turning line are found to be well emulated. For the spontaneous fission half-lives the NN emulation provides values that are found to agree with the DFT predictions within a factor of 10 3 across more than 70 orders of magnitude. Neural networks are able to emulate the potential energy and collective inertia well enough to reasonably predict physical observables. Future directions of study, such as the inclusion of additional collective degrees of freedom and active learning, will improve the predictive power of microscopic theory and further enable large-scale fission studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Quadrupole strength in isobaric triplets

The dependence of the 𝐸⁢2 matrix elements on isospin projection 𝑇 𝑧 is linked to the conservation of the isospin symmetry. To study this conjecture, we calculated the 𝐵⁡(𝐸⁢2 : 2 + →0 + ) rates for the even-even 𝑇=1 mirror nuclei with 42 ≤ 𝐴 ≤ 98 within nuclear density functional theory, employing the generalized Bohr Hamiltonian, and carrying out angular momentum projection. We demonstrated that collective effects are crucial for describing experimental data near the 𝑁=𝑍 line without invoking explicit beyond-Coulomb isospin symmetry-breaking corrections. We also determined the 𝐵⁡(𝐸⁢2↓) values for odd-odd 𝑇 𝑧 =0 nuclei 70 Br and 78 Y in doubly blocked configurations. We discussed the requirements for accurately describing isobaric analog states and emphasized how current theoretical results should be interpreted within the study of isospin symmetry across isospin triplets.

collective levels↗
DOE OSTI2026

Mass of 101 Sn and Bayesian extrapolations to the proton drip line

The favorable energy configurations of nuclei at magic numbers of 𝑁 neutrons and 𝑍 protons are fundamental for understanding the evolution of nuclear structure. The 𝑍 = 50 (tin) isotopic chain is a frontier for such studies, with particular interest at and around the doubly magic 100 Sn isotope, for which the mass is a topic of debate. Precise mass values for neutron-deficient isotopes provide necessary anchor points for mass models to test extrapolations near the proton drip line, where experimental studies remain out of reach. In this work, we report a Penning trap mass measurement of 101 Sn . The determined mass excess of −59889.89⁢(96) keV for 101 Sn represents a factor-of-300 improvement over the current precision and indicates that 101 Sn is less bound than previously thought. Mass predictions from a recently developed Bayesian model combination framework employing statistical machine learning and nuclear masses computed within seven global models based on nuclear density functional theory agree within 1⁢𝜎 with experimental masses from the 48 ≤ 𝑍 ≤ 52 isotopic chains. The framework's resilience to new mass data gave confidence in the extrapolation of tin masses down to 𝑁 = 46. Our calculations suggest that 96 Sn is a two-proton drip line nucleus and predict a mass excess of −58090⁢(800) keV for 100 Sn , showing a preference within 1⁢𝜎 for the mass of 100 Sn derived from the 𝛽-delayed 𝑄 value measured at GSI.

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↗
DOE OSTI2022

Evidence of Two-Source King Plot Nonlinearity in Spectroscopic Search for New Boson

Optical precision spectroscopy of isotope shifts can be used to test for new forces beyond the standard model, and to determine basic properties of atomic nuclei. We measure isotope shifts on the highly forbidden 2 S 1/2 → 2 F 7/2 octupole transition of trapped 168,170, 172, 174, 176 Yb ions. When combined with previous measurements in Yb + and very recent measurements in Yb, the data reveal a King plot nonlinearity of up to 240σ. Furthermore, the trends exhibited by experimental data are explained by nuclear density functional theory calculations with the Fayans functional. We also find, with 4.3σ confidence, that there is a second distinct source of nonlinearity, and discuss its possible origin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗