Search NASA⌕ Search

SEARCH · Search NASA

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

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2025

Electron capture of superheavy nuclei with realistic lepton wave functions

The superheavy nuclei push the periodic table of the elements and the chart of the nuclides to their limits, providing a unique laboratory for studies of the electron-nucleus interactions. The most important weak decay mode in known superheavy nuclei is electron capture (EC). In the standard calculations of EC, the lepton wave functions are usually considered in the lowest-order approximation. In this work, we investigate the sensitivity of EC rates on the choice of the electron wave functions by (i) assuming the single-particle approximation for the electron wave functions, and (ii) carrying out Dirac-Hartree-Fock (DHF) calculations. The nuclear response is generated based on the state-of-the-art quasiparticle random phase approximation employing relativistic nuclear energy density functional theory. Here, we show that using the improved lepton wave functions reduces the EC rates up to 40% in the superheavy nucleus oganesson (𝑍=118). Interestingly, because of screening effects, the difference between the EC rates obtained with the DHF and single-particle calculations is fairly small.

Atomic orbital↗
DOE OSTI2025

Time-Dependent Density Functional Theory Description of 238 U⁡(n,f), 240,242 Pu⁢(n,f), and 237 Np(n,f) Reactions

In nuclei with an odd nucleon number the nonvanishing spin number density is the source of a pseudomagnetic field, which favors the splitting of the nucleon Cooper pairs. Such a pseudomagnetic field is generated always in the dynamics of any nucleus, but its effects on Cooper pairs are significantly enhanced in the dynamic evolution of nuclei with an odd number of nucleons. We present for the first time a microscopic study of the induced fission of the odd neutron compound nuclei 239 U, 241,243 Pu, and the odd proton, odd neutron compound nucleus 238 Np, performed within the time-dependent density functional theory extended to superfluid fermion systems, without any simplifying assumptions, with controlled numerical approximations, and for a very large number of initial conditions. Because of the presence of the unpaired odd nucleon(s), the time-reversal symmetry of the fission compound nucleus is spontaneously broken, an aspect routinely neglected in the most advanced microscopic approaches of the past. The emerging fission fragment properties are quite similar to the properties of fission fragments of neighboring even-even nuclei. The time from saddle-to-scission is often significantly longer in odd-odd or odd-mass nuclei than for even-even nuclei, since systems with unpaired nucleons are easier to excite and the potential energy surfaces of these nuclei have more structure, often resembling a very complicated obstacle course, rather than a more direct evolution of the nuclear shape from the top of the outer fission barrier to the scission configuration. The Pauli blocking approximation, often invoked in the literature, expected to inhibit the fission of nuclei with unpaired nucleons, is surprisingly strongly violated during the fission dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Charged Defects in UO 2 Bulk and Surface: A First-Principles Study

Uranium dioxide (UO 2 ) is the primary fuel used in nuclear reactors. Under the extreme heat and radiation inside a reactor, this material inevitably develops defects in its crystal structure. To investigate the nature and behavior of these defects, DFT+U calculations were employed to investigate charged point defects in both bulk UO 2 and its most stable surface, the (111) plane. The formation of defects and their impact on the electronic structure were systematically examined. The results reveal that these defects introduce localized electronic states, alter magnetic behavior, and modify the structural properties. In general, such defects act as deep traps capable of capturing and retaining charge carriers. The stability of these defects depends strongly on the chemical environment and the position of the Fermi level. Surface defect calculations reveal that oxygen vacancies form more readily at the surface than in the bulk over a wide range of electron chemical potential, with subsurface oxygen vacancies being more stable than those in the top layer. Overall, the findings demonstrate how charged defects influence magnetism, transport, and stability in UO 2 , providing insights that may guide improvements in the safety and efficiency of nuclear fuel.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Further steps toward the next generation of covariant energy density functionals

The present study aims at further development of covariant energy density functionals (CEDFs) towards more accurate description of binding energies across the nuclear chart. Infinite basis corrections to binding energies in the fermionic and bosonic sectors of the covariant density functional theory are taken into account in the fitting protocol within the covariant density functional theory. In addition, total electron binding energies are used in the conversion of atomic binding energies into nuclear ones. Their dependence on neutron excess is investigated across the nuclear chart within the atomic approach. Furthermore, these factors were disregarded in the previous generation of covariant energy density functionals, but their omission leads to substantial global calculation errors for physical quantities of interest. For example, these errors for binding energies are of the order of 0.8 MeV or higher for the three major classes of covariant energy density functionals.

Binding energy & masses↗
DOE OSTI2025

Constraints on magnetism and correlations in RuO 2 from lattice dynamics and Mössbauer spectroscopy

We provide experimental evidence for the absence of a magnetic moment in bulk RuO 2 , a candidate altermagnetic material, by using a combination of Mössbauer spectroscopy, nuclear forward scattering, inelastic X-ray and neutron scattering, and density functional theory calculations. Using complementary Mössbauer and nuclear forward scattering, we determine the Ru magnetic hyperfine splitting to be negligible. Inelastic X-ray and neutron scattering-derived lattice dynamics of RuO 2 are compared to density functional theory calculations of varying flavors. Comparisons among theory with experiments indicate that electronic correlations, rather than magnetic order, are key in describing the lattice dynamics.

Mössbauer spectroscopy↗
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 OSTI2024

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Abstract The improvement of Pt‐based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn : Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of Pt x Zn y ( x =1–3, y=0–3) active sites grafted on silanols of Silicalite‐1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt 1 Zn y sites, whose activity is similar for y=0–2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt 2 Zn 1 and Pt 3 Zn 1 sites are stable and exhibit activity similar to Pt 1 Zn 2 , but only Pt 1 Zn 2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Liu, Yilang↗
DOE OSTI2024

Propane Dehydrogenation on Pt x Zn y Active Sites in Silicalite‐1

Abstract The improvement of Pt‐based catalysts for propane dehydrogenation (PDH) has progressed by recent investigations that have identified Zn as a promising promoter for Pt subnanometer catalysts. It is desirable to gain insights into the structure, stability, and activity of such active sites and the factors that influence them, such as Zn : Pt ratio, Pt coordination and nuclearity. Here, we employ density functional theory and microkinetic simulations to investigate the stability of Pt x Zn y ( x =1–3, y=0–3) active sites grafted on silanols of Silicalite‐1 and the PDH activity of Pt. We find that the coordination of a Pt atom to a nest of grafted Zn(II) atoms increases the stability of the Pt 1 Zn y sites, whose activity is similar for y=0–2 and drops dramatically for y>2. We further demonstrate, via linear scaling relations and microkinetic simulations, that the turnover frequency obeys a volcano law as a function of propylene binding strength. The Pt 2 Zn 1 and Pt 3 Zn 1 sites are stable and exhibit activity similar to Pt 1 Zn 2 , but only Pt 1 Zn 2 manifests reaction kinetics consistent with experimental data, strongly suggesting the active site composition in the synthesized catalyst samples. The methodology presented here suggests a general strategy for deducing active site information such as composition through simple kinetic experiments.

Liu, Yilang↗
DOE OSTI2025

Speciation and diffusive dynamics in hydrated grain boundaries of complex oxide Gd2Ti2O7

Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.

B. Ghosh, Dipta↗
DOE OSTI2025

The Electron Thermal Conductivity of Pu and Zr Substituted $\mathcal{γ}$-U

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (𝑘) and fissile density. Limited experimental studies of the 𝑘 of U-Pu-Zr alloys in the range of 15 to 20 wt% Pu and 6 to 15 wt% Zr indicate that increasing the content of either Zr or Pu tends to lower 𝑘. However, which element has the greater effect on 𝑘, and the associated mechanisms, remains unclear. Here, in this study, the electron thermal conductivity (𝑘 𝑒 ) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the 𝑘 𝑒 of 𝛾-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate 𝑘 𝑒 ; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the more accurate method is that which employs the electron relaxation time of 𝛾-U, while if the elements are dissimilar, a mixed method that mixes several parameters associated with JNW_S⁢3033426825100132 from each element in the alloy is best. The introduction of all alloying elements decreases 𝑘 𝑒 ; however, in binary compounds, Pu and Zr have different effects. Pu flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases 𝑘 𝑒 more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel without sacrificing 𝑘. We also note that these 𝑘 𝑒 calculation methods can be applied to non-fuel alloys that require 𝑘 𝑒 predictions, such as cladding, heat exchanger, and structural materials.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Density functional theory calculations of the mixing enthalpy of ternary uranium carbide compounds

The high melting point of uranium-zirconium carbides (U,Zr)C makes them an ideal fuel for nuclear thermal propulsion (NTP) reactors. Gaps remain in the current understanding of the U-Zr-C system due to the difficulty of conducting thermodynamic experiments at NTP operation conditions. Density functional theory calculations using the Hubbard U model (DFT+U) were performed using orbital matrix occupation (OMC) to obtain the mixing enthalpy for UC and ZrC for (U,Zr)C ternary compounds. Similarly, DFT+U calculations were also carried out for the (U,Nb)C and (U,Ta)C systems. In conclusion, the DFT results are envisioned to be used in thermodynamic assessments of the uranium carbide systems based on the CALPHAD approach to supplement the lack of experimental data for the mixing thermodynamics.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Structure and Dynamics of Imidazolium in an Ionic Liquid-PEGDA Iongel via IR, 2D-IR, and NMR Spectroscopy

Fourier transform infrared (FTIR) absorption spectroscopy, ultrafast two-dimensional infrared (2D-IR) spectroscopy, nuclear magnetic resonance (NMR) spectroscopy, and density functional theory (DFT) calculations show that, in iongels composed of 1-ethyl-3-methylimidazolium bis(trifluoromethyl\-sulfonyl)imide ([C 2 C 1 Im][Tf 2 N]}) and cross-linked poly(ethylene glycol)diacrylate (cl-PEGDA), imidazolium cations solvate PEGDA chains in tight and loose complexes. H/D isotope substitution at the 2-position of the imidazolium ring simplifies the CH-stretching vibrational bands. Tight complexes are characterized by stronger, directional hydrogen bonding interactions between the 2-position of the imidazolium ring and the ethereal oxygens. The loose complexes resemble bulk ionic liquid (IL) non-specifically perturbed by the PEGDA chains.

ionic liquids, carbon dioxide, cross-linked PEGDA,↗
DOE OSTI2025

Tuning the Properties of Lignin‐Derived Deep Eutectic Solvents for Biomass Processing

Lignin-derived deep eutectic solvents (DESs) have been investigated as sustainable green media for biomass processing. However, the properties and processability of DESs have not been fully understood with the chemical structures of their constituents for biomass fractionation. Here, in this article, the properties of the phenolic DESs are discussed with different numbers of functional groups, such as –OCH 3 and –CHO in their hydrogen bond donor (HBD) structures. The formation of DES is significantly related to the hydrogen bond between its constituents, identified by nuclear magnetic resonance (NMR) analysis and density functional theory calculation (DFT). Lower viscosity and net basicity of DES are achieved with fewer –OCH 3 on HBD structures, resulting in enhanced processability and fractionation efficiency. The thermal stability of the DES is also influenced by the –OCH 3 and –CHO of HBD, as indicated by its onset temperature. The recyclability of the phenolic DES is confirmed by the fractionation performance of the recycled DES. Understanding the structural impacts of DES constituents on the properties and performance is crucial for designing solvents in biorefinery applications.

Ryu, Jiae [State University of New York (SUNY), Sy↗
DOE OSTI2026

Uncovering the linear boron environment in Na 3 BP 2 through solid-state 11 B NMR spectroscopy

Boron-based compounds exhibit a wide range of structural diversity, with potential applications spanning organic and inorganic chemistry. Herein, we focus on the characterization of the linear boron-phosphorus unit P═B═P in Na 3 BP 2 using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and density functional theory (DFT) calculations. High-resolution 11 B ssNMR spectra were recorded at two fields, and key parameters such as chemical shift anisotropy (CSA), quadrupolar coupling constants (C Q ), and electric field gradient (EFG) tensors were extracted. The 11 B NMR results revealed a distinct chemical environment for the two-coordinate boron atom, with a CSA span (Ω) of 280 ppm and a C Q of 3.0 MHz. These values were further validated through periodic plane-wave DFT calculations, which showed good agreement with experimental results. The obtained spectral parameters are compared to other linear boron units, such as the BO 2 motif, providing a broader context for understanding boron coordination in inorganic compounds. This work expands the body of NMR knowledge on boron-containing materials, particularly for linear boron motifs. The findings contribute to the growing field of boron chemistry and its potential applications in advanced materials.

Porter, Andrew P. [Ames Laboratory (AMES), Ames, I↗
DOE OSTI2024

Structural Studies of Alloyed and Nanoparticle Transition Metal Dichalcogenides by Selenium-77 Solid-State Nuclear Magnetic Resonance Spectroscopy

Layered transition metal dichalcogenides (TMDCs) such as MoS 2 , MoSe 2 , and WSe 2 are under intense investigation because they are atomically thin semiconductors with photophysical properties that can be tuned by changing their composition or morphology. Mechanochemical processing has been proposed as a method to obtain alloyed TMDCs in the series Mo 1–x W x S y Se 2–y (x = 0–1; y = 0–2). However, elucidating the chemical transformations occurring at the atomic scale following mechanochemical processing can be challenging because the products are often amorphous or microcrystalline. To address this challenge, we probe TMDC mixing and alloying by using a combination of powder X-ray diffraction, Raman spectroscopy, diffuse reflectance spectroscopy, 77 Se solid-state nuclear magnetic resonance (SSNMR) spectroscopy, and planewave density functional theory (DFT) calculations. The nature of the milling material and reaction atmosphere are shown to be essential factors in limiting the formation of undesired oxide byproducts. We demonstrate acquisition of 77 Se SSNMR spectra using different combinations of Carr-Purcell Meiboom-Gill acquisition (CPMG) pulse sequences, magic angle spinning (MAS), and MAS dynamic nuclear polarization. Further, the combination of SSNMR with the other characterization methods clearly demonstrates that high energy impact ball milling induces molecular level alloying of Mo, W and chalcogen atoms in the family Mo 1–x W x S y Se 2–y . Gauge including projector augmented wave DFT calculations yield accurate 77 Se chemical shift (CS) tensor components. 77 Se SSNMR spectroscopy was also applied to study the structure of WSe 2 nanocrystals intercalated with ethylenediamine. The intercalated WSe 2 nanocrystals exhibit a more positive isotropic 77 Se CS as compared to bulk WSe 2 , however, the 77 Se CS anisotropy is the same, confirming the WSe 2 layers have a similar structure as in their bulk counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Localizing tetrahedral aluminum in nitrate-bearing gibbsite to constrain defect-impurity coupling

The enhanced radiolytic stability of gibbsite (α-Al(OH) 3 ) containing trace nitrate (NO 3 − ) is a phenomenon in nuclear waste management, but its structural origins remain unresolved. Motivated by the detection of minority tetrahedral aluminum (T d ) defects in synthetic gibbsite, we hypothesized that these sites may participate in NO 3 − retention or mediate H 2 suppression. To evaluate this, we combined orthogonal techniques comprised of spatially selective solid-state 27 Al MAS NMR, comparative spectroscopy, and density functional theory (DFT) modeling. Paramagnetic editing and dynamic nuclear polarization (DNP) MAS NMR confirm that T d defects are confined to the particle interior. DFT calculations reveal no energetic stabilization of NO 3 − near T d sites. Comparative NMR analysis shows that T d is also present in chloride-bearing gibbsite, which exhibits high radiolytic hydrogen yields. These three independent disqualifications rule out T d as a structural contributor to nitrate-mediated suppression and narrow the scope of defect-driven explanations. The findings redirect mechanistic attention away from coordination defects and toward redox-active impurity pathways, providing a refined foundation for understanding radiation tolerance in Al(OH) 3 .

Graham, Trent R. [Pacific Northwest National Labor↗
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↗