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

Weak decays in superheavy nuclei

Superheavy nuclei represent the extreme atoms and nuclides known at the limit of mass and charge. The observed superheavy nuclei are all proton-rich; they decay primarily by emitting 𝛼 particles and by fission with a possible small electron capture (EC) branch. Here, due to the huge atomic numbers and associated relativistic effects, EC decays of superheavy systems are expected to differ from what is known in lighter nuclei. In this letter, using the quantified relativistic nuclear density functional theory and the quasiparticle random-phase approximation with the interaction optimized to experimental EC/𝛽 ± -decay half-lives, and Gamow-Teller resonance energies, we study the EC/𝛽 ± -decays in 𝑍=101–118 nuclei. Both allowed (1 + ) and first-forbidden (0 − ,1 − and 2 − ) transitions are considered. We show that the first-forbidden 1 − transitions dominate the decay rates in almost all studied nuclei. For proton-rich nuclei, EC dominates over 𝛽 + decay. Based on calculations with two relativistic energy density functionals, we identify 45 candidate nuclei in which a competition between weak decays and 𝛼 decay and spontaneous fission is expected.

A ≥ 220↗
DOE OSTI2025

Tritium adsorption and absorption on (100) and (001) surfaces of pure and tin defective zirconium

Zirconium alloys such as zircaloy-4 are used as tritium (T) getter materials in tritium-producing burnable absorber rods (TPBARs) due to their ability to capture T, thereby forming metal hydrides. Developing an understanding of T adsorption onto zircaloy prior to diffusion into the subsurface is relevant for rational tritium getter and TPBAR design, to improve material properties for nuclear applications. Herein, density functional theory calculations revealed the preferred binding sites for T adsorption on Zr(001) and Zr(100). The energy barriers of T transfer, along the surface and from the surface to the subsurface were computed. The adsorption properties of Zr(001) were found to be superior to those of Zr(100). Surface tin impurities were found to strongly repel T. The presence of subsurface and surface tin resulted in higher absorption energy barriers for both the forward and reverse processes. Based on the calculated energy barriers, a surface to surface T diffusion coefficient of 9.53 × 10 -10 m 2 s -1 is expected for pristine Zr(001). Finally, a surface to subsurface T diffusion coefficient on the order of 10 -13 m 2 s -1 is predicted in pristine Zr, decreasing to 10 -19 m 2 s -1 for the transfer with a subsurface tin impurity.

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

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Xenon Trapping in Metal-Supported Silica Nanocages

Xenon (Xe) is a valuable and scarce noble gas used in various applications, including lighting, electronics, and anesthetics, among many others. It is also a volatile byproduct of the nuclear fission of uranium. In this work, a novel material architecture consisting of silicate nanocages in contact with a metal surface and an approach for trapping single Xe atoms in these cages is presented. The trapping is done at low Xe pressures and temperatures between 400 and 600 K, and the process is monitored in situ using synchrotron-based ambient pressure X-ray photoelectron spectroscopy. Release of the Xe from the cages occurs only when heating to temperatures above 750 K. A model that explains the experimental trapping kinetics is proposed and tested using Monte Carlo methods. Density functional theory calculations show activation energies for Xe exiting the cages consistent with experiments. This work can have significant implications in various fields, including Xe production, nuclear power, nuclear waste remediation, and nonproliferation of nuclear weapons. The results are also expected to apply to argon, krypton, and radon, opening an even more comprehensive range of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Nuclear Thermal Propulsion Material Trade for Additively Manufactured Regenerative Chambers

Nuclear thermal propulsion (NTP) is a game changing technology with advantages over chemical propulsion systems for deep space missions but comes with many engineering challenges. This paper focuses on the challenge of material selection for a regeneratively-cooled chamber for use in a nuclear environment with minimal design margins. Chemical compatibility with reactants and products following energy extraction (be it through combustion or neutron bombardment) must be considered. Similarly, mechanical and thermophysical properties (e.g., density, thermal conductivity, yield strength, etc.) as functions of temperature are critical. While hydrogen embrittlement environment (HEE) is a risk and consideration for traditional chemical combustion chambers, application in a nuclear thermal engine adds additional material compatibility and resilience requirements to long-term nuclear radiation exposure. Chamber materials in this trade included, but were not limited to, A-286, Haynes 230, Haynes 282, CRES 347, GRCop-42, GRCop-84, Inconel 625, JBK-75, and NASA HR-1. The material trade focuses on modern additive manufacturing techniques and hypothetical geometric constraints.

Nuclear Thermal Propulsion↗
NASA NTRSDate not supplied

Nuclear Thermal Propulsion Material Trade for Additively Manufactured Regenerative Chambers

Nuclear thermal propulsion (NTP) is a game changing technology with advantages over chemical propulsion systems for deep space missions but comes with many engineering challenges. This paper focuses on the challenge of material selection for a regeneratively-cooled chamber for use in a nuclear environment with minimal design margins. Chemical compatibility with reactants and products following energy extraction (be it through combustion or neutron bombardment) must be considered. Similarly, mechanical and thermophysical properties (e.g., density, thermal conductivity, yield strength, etc.) as functions of temperature are critical. While hydrogen embrittlement environment (HEE) is a risk and consideration for traditional chemical combustion chambers, application in a nuclear thermal engine adds additional material compatibility and resilience requirements to long-term nuclear radiation exposure. Chamber materials in this trade included, but were not limited to, A-286, Haynes 230, Haynes 282, CRES 347, GRCop-42, GRCop-84, Inconel 625, JBK-75, and NASA HR-1. The material trade focuses on modern additive manufacturing techniques and hypothetical geometric constraints.

Nuclear Thermal Propulsion↗
DOE OSTI2020

Solvent manipulation of the pre-reduction metal–ligand complex and particle-ligand binding for controlled synthesis of Pd nanoparticles

Understanding how to control the nucleation and growth rates is crucial for designing nanoparticles with specific sizes and shapes. Herein, we show that the nucleation and growth rates are correlated with the thermodynamics of metal–ligand/solvent binding for the pre-reduction complex and the surface of the nanoparticle, respectively. To obtain these correlations, we measured the nucleation and growth rates by in situ small angle X-ray scattering during the synthesis of colloidal Pd nanoparticles in the presence of trioctylphosphine in solvents of varying coordinating ability. The results show that the nucleation rate decreased, while the growth rate increased in the following order, toluene, piperidine, 3,4-lutidine and pyridine, leading to a large increase in the final nanoparticle size (from 1.4 nm in toluene to 5.0 nm in pyridine). Using density functional theory (DFT), complemented by 31 P nuclear magnetic resonance and X-ray absorption spectroscopy, we calculated the reduction Gibbs free energies of the solvent-dependent dominant pre-reduction complex and the solvent-nanoparticle binding energy. The results indicate that lower nucleation rates originate from solvent coordination which stabilizes the pre-reduction complex and increases its reduction free energy. At the same time, DFT calculations suggest that the solvent coordination affects the effective capping of the surface where stronger binding solvents slow the nanoparticle growth by lowering the number of active sites (not already bound by trioctylphosphine). The findings represent a promising advancement towards understanding the microscopic connection between the metal–ligand thermodynamic interactions and the kinetics of nucleation and growth to control the size of colloidal metal nanoparticles.

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

In-situ Synthesis of PuCl 3 and Corresponding Raman and Density Functional Theory Vibrational Modes

The experimental and calculated Raman spectrum of PuCl 3 has been reported for the first time. PuCl 3 is a primary species found in plutonium metal refinement, specifically in pyrochemical salt processes including multicycle direct oxide reduction, metal chlorination, and electrorefining. As such, Raman signatures of PuCl 3 could serve as potential forensic indicators of material process history. A novel technique for synthesizing PuCl 3 from the in situ chlorination of plutonium metal with HCl was developed to establish these signatures. Cerium metal surrogates were utilized to ensure optimization of the plutonium experiments and to minimize personnel exposure, and all experiments were carried out in a Raman reaction chamber designed for air-tight, high vacuum environments. In situ Raman spectroscopy was employed in conjunction with density functional theory (DFT) to investigate the vibrational modes of PuCl 3 . Associated mixed oxy and hydroxyl phases are also reported. The combined Raman and DFT results have eliminated inconsistencies in Raman mode assignments for the MX 3 family of metal chlorides having P6 3 /m symmetry, and IR modes derived from the DFT calculations are additionally presented. The data observed in this study are of potential interest to nuclear forensic analyses, nuclear sample aging, nuclear energy, plutonium processing, and stockpile stewardship.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Initializing BSQ with Open-Source ICCING

While it is well known that there is a significant amount of conserved charges in the initial state of nuclear collisions, the production of these due to gluon splitting has yet to be thoroughly investigated. The ICCING (Initial Conserved Charges in Nuclear Geometry) algorithm reconstructs these quark distributions, providing conserved strange, baryon, and electric charges, by sampling a given model for the g → qq¯ splitting function over the initial energy density, which is valid at top collider energies, even when µB = 0. The ICCING algorithm includes fluctuations in the gluon longitudinal momenta, a structure that supports the implementation of dynamical processes, and the c++ version is now open-source. A full analysis of parameter choices on the model has been done to quantify the effect these have on the underlying physics. We find there is a sustained difference across the different charges that indicates sensitivity to hot spot geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Role of isospin composition in low-energy nuclear fusion

We employ a microscopic approach that examines the impact of isospin dynamics on the process of low- energy nuclear fusion along an isotope chain and dependence on deformation. Our method utilizes the density constrained time-dependent Hartree-Fock theory (DC-TDHF), where isoscalar and isovector characteristics of the energy density functional (EDF) are examined in turn. This approach is applied to a series of fusion interactions of 176 Yb with increasingly neutron rich isotopes of calcium. By evaluating the contributions from the isoscalar and isovector components of the EDF, we look to quantify the influence of isospin composition on the conditions under which fusion is most likely to take place. Furthermore, our findings reveal that, in nonsymmetric systems, the isovector dynamics plays a significant role. Its typical effect is a reduction in the potential barrier, which turns into enhancement for neutron-rich systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Axially-deformed solution of the Skyrme-Hartree-Fock-Bogoliubov equations using the transformed harmonic oscillator basis (IV) HFBTHO (v4.0): A new version of the program

We describe the new version 4.0 of the code HFBTHO that solves the nuclear Hartree-Fock-Bogoliubov problem by using the deformed harmonic oscillator basis in cylindrical coordinates. In the new version, we have implemented the restoration of rotational, particle number, and reflection symmetry for even-even nuclei. The restoration of rotational symmetry does not require using bases closed under rotation. Furthermore, we added the SeaLL1 functional and improved the calculation of the Coulomb potential. Finally, we refactored the code to facilitate maintenance and future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Static and dynamic properties of atomic nuclei with high-resolution potentials

Here, we compute ground-state and dynamical properties of 4 He and 16 O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schrodinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the 4 He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of 16 O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Developments in Physics-Based Modeling and Machine Learning for Environmental Effects in Nuclear Materials [Slides]

This presentation is part of an annual review for the new Advanced Materials and Manufacturing Technologies program under the Office of Nuclear Energy within the Department of Energy. This presentation covers advances in physics-based modeling and machine learning for environmental effects in nuclear materials. It discusses efforts in combining density functional theory with machine learning to predict material properties, studying irradiation defect recombination, and improving MOOSE-based tools for crystal plasticity modeling and the Stochastic Tools Module.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Re-evaluation of the TSL for Yttrium Hydride

Yttrium hydride (YH x ) is of interest as a high-temperature moderator material in advanced nuclear reactor systems because of its superior ability to retain hydrogen at elevated temperatures. Thermal neutron scattering laws (TSL) for hydrogen bound in yttrium hydride (H-YH 2 ) and yttrium bound in yttrium hydride (Y-YH 2 ) were previously evaluated by Naval Nuclear Laboratory using the ab initio approach and released in ENDF/B-VIII.0. In that work, density functional theory, incorporating the generalized gradient approximation (GGA) for the exchange-correlation energy, was used to simulate the face-centered cubic structure of YH 2 and calculate the interatomic Hellmann-Feynman forces for a 2×2×2 supercell containing 96 atoms. Lattice dynamics calculations using PHONON were used to determine the phonon density of states. The calculated phonon density of states for H and Y in YH 2 were then used to prepare H-YH 2 and Y-YH 2 TSL evaluations, in the incoherent approximation, using the LEAPR module of NJOY. In addition, elastic scattering was assumed to be incoherent for both H and Y. While the incoherent elastic scattering approximation is appropriate for H-YH 2 , it introduces an undesirable approximation for Y-YH 2 . In this work, we re-evaluate the TSL for Y-YH 2 using FLASSH (Full Law Analysis Scattering System Hub). Y-YH 2 is evaluated using the FLASSH generalized coherent elastic scattering capability in order to capture the Bragg peaks associated with the YH 2 crystal structure which were neglected in the prior NJOY-based evaluation due to limitations in LEAPR. An experimental approach to validate the Y-YH 2 TSL using neutron transmission measurements is discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Strong Proton‐Phonon Coupling Drives Fast Ion Transport in Perovskites

Conduction of protons in solids is a cooperative process propelled by phonons, with molecular details obscured by the irregular movements in the thermal bath. It is shown that substitution with Y forms an imaginary phonon mode, instrumental for the function as proton conductor and effectively lowering the activation barrier for proton transport. To untangle the interplay in the exemplary proton conductor BaSn 0.9 Y 0.1 O 3 , its crystallographic structure is determined with high resolution neutron diffractometry and its phonon density of states with density functional theory calculations, experimentally validated by element specific nuclear resonant vibration spectroscopy. Based on phonon analysis, a quantitative transport model is present, which predicts the activation energy and performance by the ratio of ionic radii. Rather than individual vibrational modes, it is the oxygen sub-lattice which exerts its momentum on the protons. The extent of this momentum transfer is governed by the ratio of ionic radii. This model extends the transition state theory by the phonon-phonon interaction and complements the previously proposed idea that lattice dynamics is decisive for proton transport and specifies which properties of the material exactly define the vibration properties.

activation barriers↗
DOE OSTI2023

Evidence of Hexadecapole Deformation in Uranium-238 at the Relativistic Heavy Ion Collider

State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic 238 U + 238 U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding 238 U ions. Here, we show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like 238 U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of 238 U on high-energy collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗