Search NASA⌕ Search

SEARCH · Search NASA

Results for “nuclear energy density functional”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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 OSTI2022

Structure Functions and Parton Densities: a Session Summary

Studies of fragmentation and parton density functions are a core component of research in high energy particle and nuclear physics. These quantities are inherently interesting as a probe of the quantum nature of the strong force and are also essential ingredients to additional studies in high energy scattering experiments. These proceedings provide an overview of the state of the art in this area, as presented at the Deep-Inelastic Scattering Conference in the Spring of 2021.

Nachman, Benjamin↗
DOE OSTI2024

Using DAPPER to extract the photon strength function of 58Fe using the inverse Oslo and shape methods

The photon strength function of 58 Fe has been extracted using both the Oslo and Shape methods from particle-γ coincidence data measured using the Detector Array for Photons, Protons, and Exotic Residues, which probes nuclei utilizing (d,p) reactions in inverse kinematics. Four particle-γ coincidence matrices, each constructed with different treatments of the γ-ray energies, are explored in order to observe the impact on the resulting nuclear level density and photon strength. The final photon strength function reported is found to agree well with previous Oslo measurements of other iron isotopes. Systematic uncertainties are included, using different model parameters and their reported errors to perform the Oslo method normalization. The model-independent Shape method is explored and the functional form of the photon strength function obtained is in agreement with the Oslo method results. A low-energy enhancement is not reported for 58 Fe in this work given possible subtraction issues originating from strongly populated states.

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 OSTI2021

Optimization and supervised machine learning methods for fitting numerical physics models without derivatives

Here, we address the calibration of a computationally expensive nuclear physics model for which derivative information with respect to the fit parameters is not readily available. Of particular interest is the performance of optimization-based training algorithms when dozens, rather than millions or more, of training data are available and when the expense of the model places limitations on the number of concurrent model evaluations that can be performed. As a case study, we consider the Fayans energy density functional model, which has characteristics similar to many model fitting and calibration problems in nuclear physics. We analyze hyperparameter tuning considerations and variability associated with stochastic optimization algorithms and illustrate considerations for tuning in different computational settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Using DAPPER to extract the photon strength function of 58 Fe using the inverse Oslo and shape methods

The photon strength function of 58 Fe has been extracted using both the Oslo and Shape methods from particle–γ coincidence data measured using the Detector Array for Photons, Protons, and Exotic Residues, which probes nuclei utilizing (d,p) reactions in inverse kinematics. Four particle–γ coincidence matrices, each constructed with different treatments of the γ–ray energies, are explored in order to observe the impact on the resulting nuclear level density and photon strength. The final photon strength function reported is found to agree well with previous Oslo measurements of other iron isotopes. Systematic uncertainties are included, using different model parameters and their reported errors to perform the Oslo method normalization. The model-independent Shape method is explored and the functional form of the photon strength function obtained is in agreement with the Oslo method results. A low-energy enhancement is not reported for 58 Fe in this work given possible subtraction issues originating from strongly populated states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗
DOE OSTI2021

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

The TDHF code Sky3D version 1.2

The Sky3D code has been widely used to describe nuclear ground states, collective vibrational excitations, and heavy-ion collisions. The approach is based on Skyrme forces or related energy density functionals. The static and dynamic equations are solved on a three-dimensional grid, and pairing is been implemented in the BCS approximation. This updated version of the code aims to facilitate the calculation of nuclear strength functions in the regime of linear response theory, while retaining all existing functionality and use cases. The strength functions are benchmarked against available RPA codes, and the user has the freedom of choice when selecting the nature of external excitation (from monopole to hexadecapole and more). Some utility programs are also provided that calculate the strength function from the time-dependent output of the dynamic calculations of the Sky3D code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Material Effects on Electron-Capture Decay in Cryogenic Sensors

Several current searches for physics beyond the standard model are based on measuring the electron-capture (EC) decay of radionuclides implanted into cryogenic high-resolution sensors. The sensitivity of these experiments has already reached the level where systematic effects related to atomic state energy changes from the host material are a limiting factor. One example is a neutrino mass study based on the nuclear EC decay of 7 Be to 7 Li inside cryogenic Ta-based sensors. Here, to understand the material effects at the required level, we use density-functional theory to model the electronic structure of lithium atoms in different atomic environments of the polycrystalline Ta absorber film. The calculations reveal that the Li 1s binding energies can vary by more than 2 eV due to insertion at different lattice sites, at grain boundaries, in disordered Ta, and in the vicinity of various impurities. However, the total range of Li 1s shifts does not exceed 4 eV, even for extreme amorphous disorder. Furthermore, when investigating the effects on the Li 2s levels, we find broadening of more than 5 eV due to hybridization with the Ta band structure. Material effects are shown to contribute significantly to peak broadening in Ta-based sensors that are used to search for physics beyond the standard model in the EC decay of 7 Be, but they do not explain the full extent of observed broadening. Understanding these in-medium effects will be required for current- and future-generation experiments that observe low-energy radiation from the EC decay of implanted isotopes to evaluate potential limitations on the measurement sensitivity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fission with Exotic Nuclei (Abbreviated Report)

Nuclear fission is a key mechanism involved in the synthesis of heavy elements in the Cosmos and is the primary explanation for the stability of superheavy elements. Nevertheless, our knowledge of fission remains extremely fragmented. Most experiments have been conducted only on a tiny number of stable actinide nuclei and are often incomplete, leading to gaps in our basic understanding of the process. For many radioactive isotopes, basic fission data such as the charge or mass distribution of the fragments is unknown. These gaps cannot always be filled by simulation alone. Common fission models contain too many free parameters and lack predictive power. In contrast, the fundamental theory of fission under development at LLNL is much more predictive, but its current computational cost is too high to be used extensively for data evaluations. A unique window of opportunity to resolve these limitations has recently opened: the U.S. nuclear science community is ramping up major experimental programs at the Facility for Rare Isotope Beams (FRIB, the DOE flagship facility in low-energy nuclear science), and techniques from machine learning have shown great potential to simplify the use of a fundamental, quantum-mechanical theory of fission. This project has two components. On the experimental side, we acquired and deployed at the HIGS facility a new dual Frisch-Grid ionization chamber to measure correlated fragment-mass, kinetic energy, and angular distributions of fission fragments from induced fission. This new device was used to perform measurements of charge, mass and total kinetic energy of fission fragments in the photofission of 238 U and eight gamma-ray beam energies between 6.2 and 13 MeV, which allowed extracting high-precision independent yields for this reaction. The device was also used to perform measurements of the same quantities in the neutron-induced fission of 234 U with monoenergetic beams of energy between 5 and 8 MeV. In parallel, we collaborated with a team at Commissariat à l’énergie atomique et aux énergies alternatives (CEA) to perform a series of measurements of fission yields in inverse kinematics for the two isotopes of 236 U and 240 Pu. The experiment took place at the Grand Accélérateur National d’Ions Lourds in France in June 2023. The deployment of the VAMOS spectrometer with a new array called PISTA allowed determining the excitation energy of the fissioning system within 1 Mega-electronvolts. The second component of the project involved using deep neural networks to build fast and reliable emulators of our current fission models. In an invited paper published in Frontier in Physics, we showed that autoencoders could successfully compress nuclear wavefunctions in nuclear density functional theory. We achieved a dimensionality reduction of the order of two orders of magnitude while keeping the error in the total energy to less than 0.01%. In a second paper submitted to Physical Review Letters in June 2023 with our collaborators at CEA, we showed that variational autoencoders can learn the collective degrees of freedom driving the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Sensitivity of time-dependent density functional theory to initial conditions

Time-dependent density-functional theory is mathematically formulated through nonlinear coupled time-dependent three-dimensional partial differential equations, and it is natural to expect a strong sensitivity of its solutions to variations of the initial conditions, akin to the butterfly effect ubiquitous in classical dynamics. Since the Schrödinger equation for an interacting many-body system is, however, linear and mathematically the exact equations of the density-functional theory reproduce the corresponding one-body properties, it would follow that the Lyapunov exponents are also vanishing within a density-functional theory framework. Whether for realistic implementations of the time-dependent density-functional theory the question of the absence of the butterfly effect and whether the dynamics provided is indeed a predictable theory was never discussed. At the same time, since the time-dependent density-functional theory is a unique tool allowing us to study the nonequilibrium dynamics of strongly interacting many-fermion systems, the question of predictability of this theoretical framework is of paramount importance. Here our analysis, for a number of quantum superfluid many-body systems (unitary Fermi gas, nuclear fission, and heavy-ion collisions) with a classical equivalent number of degrees of freedom O(10 10 ) and larger, suggests that its maximum Lyapunov exponents are negligible for all practical purposes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗
DOE OSTI2020

Dirac oscillator: An alternative basis for nuclear structure calculations

The isotropic harmonic oscillator supplemented by a strong spin-orbit interaction has been the cornerstone of nuclear structure since its inception more than seven decades ago. In this paper we introduce—or rather re-introduce—the “Dirac oscillator,” a fully relativistic basis that has all the desired attributes of the ordinary harmonic oscillator while naturally incorporating a strong spin-orbit coupling. To assess the power and flexibility of the Dirac oscillator basis in the solution of nuclear structure problems within the framework of covariant density-functional theory. Here, self-consistent calculations of binding energies and ground-state densities for a selected set of doubly magic nuclei are performed using the Dirac oscillator basis and are then compared against results obtained with the often-used Runge-Kutta method. Results obtained using the Dirac oscillator basis reproduce with high accuracy those derived using the Runge-Kutta method and suggest a clear path for a generalization to systems with axial symmetry. Although the harmonic oscillator with spin-orbit corrections has been the staple of the nuclear shell model since the beginning, the Dirac oscillator is practically unknown among the nuclear physics community. In this paper we illustrate the power and flexibility of the Dirac oscillator and suggest extensions to the study of systems without spherical symmetry, as required in constrained calculations of nuclear excitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Description of the multinucleon transfer mechanism for Ca 48 + Pu 244 and Kr 86 + Pt 198 reactions in a quantal transport approach

Multinucleon transfer (MNT) reactions involving heavy projectile and target combinations stand as a promising method for synthesizing new neutron-rich exotic nuclei, which may not be possible using hot or cold fusion reactions or fragmentation. Exploring the mechanisms behind MNT reactions is essential and it requires a comprehensive theoretical framework that can explain the physical observables in these reactions. This work aims to show that the quantal diffusion approach based on the stochastic mean-field (SMF) theory is capable of explaining the reaction dynamics observed in MNT reactions. Primary product mass distributions in 48 Ca + 244 Pu reaction at E c.m. = 203.2 MeV and 86 Kr + 198 Pt reaction at E c.m. = 324.2 MeV are calculated and compared with the available experimental data. In this work, we utilize the time-dependent Hartree-Fock (TDHF) calculations to analyze the mean-field reaction dynamics computationally in the reactions 48 Ca + 244 Pu and 86 Kr + 198 Pt for a broad range of initial angular momenta. Quantal transport description based on the SMF approach is used to calculate quantal diffusion coefficients and mass variances in 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary products arising from quasifission reactions are described by joint probability distribution in the SMF approach and those arising from fusion-fission are estimated by using the statistical deexcitation code gemini + +. Mean values of charge and mass numbers, scattering angles of the primary reaction products, and the total kinetic energies after the collision are calculated within the TDHF framework for a broad range of initial angular momenta. Throughout all the collisions, drift toward the mass symmetry and large mass dispersion associated with this drift are observed. Here, the calculated primary fragment and mass distributions using the SMF approach successfully explain experimental observations for the 48 Ca + 244 Pu and 86 Kr + 198 Pt systems. The primary mass distributions, mean values of binary products, and mass dispersions are determined and results are compared with the available experimental data. The observed agreement between the experimental data and SMF results highlights the effectiveness of the quantal diffusion mechanism based on the SMF approach, which does not include any adjustable parameters other than standard parameters of Skyrme energy density functional.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Xenon mobility in γ-uranium and uranium–molybdenum alloys

Diffusion in bcc uranium and U–Mo alloys is of great interest because fission gas and other fission products impact the performance of nuclear fuels. We investigate the mobility of xenon and molybdenum in bcc uranium (⁠$γ$-U) and metallic U–Mo alloys by calculating the migration energies of xenon and molybdenum for various local compositions using density functional theory. We also calculate the solute–vacancy binding energies of different solutes to vacancies in bcc uranium. We find that the solute–vacancy binding energy in bcc uranium is significantly higher than it is in other bcc metals (e.g., Fe and W). We also find that the migration energy of molybdenum is substantially higher than the migration energy of xenon, indicating that xenon is much more mobile than molybdenum in bcc uranium. The presence of molybdenum in the nearest-neighbor shell around a xenon atom typically increases the migration energy of xenon, which indicates a reduction of xenon mobility in U–Mo alloys compared to pure bcc uranium.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Bridging nuclear physics across energy scales: from neutrinoless double-beta decay to high-energy heavy-ion collisions

This paper exemplifies how connecting methods at disparate energy scales can illuminate fundamental questions. By demonstrating that nuclear wave function properties governing rare decay processes also influence collective behavior at extreme temperatures and densities, the authors have opened a new pathway for constraining physics beyond the Standard Model. With multiple ton-scale 0νββ experiments under construction, any method reducing NME uncertainties will directly impact our ability to interpret discoveries or constrain neutrino properties. The general principle—that collective phenomena in high-energy collisions can illuminate subtle features of many-body correlations in the colliding nuclei—may find applications across nuclear and particle physics. Furthermore, this intersection of nuclear structure theory, heavy-ion physics, and fundamental symmetry tests represents fertile ground for future discoveries in modern physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗