Search NASA⌕ Search

SEARCH · Search NASA

Results for “DFT+U”

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 OSTI2025

Xenon–metal pair formation in UO 2 investigated using DFT + U

A recent experimental study on a spent uranium dioxide (UO 2 ) fuel sample from Belgium Reactor 3 identified a unique pair structure formed by the noble metal phase (NMP) and fission gas [xenon (Xe)] precipitate. However, the fundamental mechanism behind this structure remains unclear. The present study aims to provide an understanding of the interaction between five different metal precipitates [molybdenum (Mo), ruthenium (Ru), palladium (Pd), technetium (Tc), and rhodium (Rh)] and the Xe fission gas atoms in UO 2 , by using density functional theory (DFT) in combination with the Hubbard U correction to compute the formation energies involved. All DFT + U calculations were performed with occupation matrix control to ensure antiferromagnetic ordering of UO 2 . The calculated formation and binding energies of the Xe and solid fission products in the NMP reveal that these metal precipitates form stable pair structures with Xe. Notably, the formation energy of Xe–metal pairs is lower than that of the isolated single defects in all instances, with Pd and Mo showing the most favorable binding energy, likely accounting for the observed pair structure formation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Correlation effects on coupled electronic and structural properties of doped rare-earth trihydrides

Rare-earth trihydride (𝑅⁢H 3 ) compounds exhibit intriguing coupled electronic and structural properties as a function of doping, hydrogen vacancies, and thermodynamic conditions. Theoretical studies of these materials typically rely on density functional theory (DFT), including the use of small supercells that may underestimate strong correlation effects and structural distortions which in turn may influence their metallicity. Here, we elucidate the roles of lattice distortions and correlation effects on the electronic properties of pristine and doped 𝑅⁢H 3 compounds by adopting DFT +U and quantum Monte Carlo (QMC) methods. Linear-response constrained DFT (LR-cDFT) methods find Hubbard U ≈ 2 eV for 𝑅 𝑑 orbitals and U ≈ 6 eV for H 𝑠 ⁡/N 𝑝 /O 𝑝 orbitals. The small U on Lu 𝑑 orbitals is consistent with QMC calculations on LuH 3 and LuH 2.875 ⁢N 0.125 . In pure face-centered-cubic (FCC) 𝑅⁢H 3 (𝑅 = Lu,Y) compounds, neither DFT nor DFT +U with the self-consistently determined U is enough to create a band gap, however a supercell with hydrogen distortions creates a small gap whose magnitude increases when performing DFT +U with self-consistently determined U values. Furthermore, correlation effects, in turn, have a moderate influence on the coupled structural and electronic properties of doped 𝑅⁢H 3 compounds and may be important when considering the competition between structural distortions and superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2024

Systematic Improvement of Quantum Monte Carlo Calculations in Transition Metal Oxides: sCI-Driven Wavefunction Optimization for Reliable Band Gap Prediction

Accurate determination of the electronic properties of correlated oxides remains a significant challenge for computational theory. Traditional Hubbard-corrected density functional theory (DFT+U) frequently encounters limitations in precisely capturing electron correlation, particularly in predicting band gaps. We introduce a systematic methodology to enhance the accuracy of diffusion Monte Carlo (DMC) simulations for both ground and excited states, focusing on LiCoO 2 as a case study. By employing a selected configuration interaction (sCI) approach, we demonstrate the capability to optimize wavefunctions beyond the constraints of single-reference DFT+U trial wavefunctions. Here, we show that the sCI framework enables accurate prediction of band gaps in LiCoO 2 , closely aligning with experimental values and substantially improving traditional computational methods. The study uncovers a nuanced mixed state of t 2g and e g orbitals at the band edges that is not captured by conventional single-reference methods, further elucidating the limitations of PBE+U in describing d-d excitations. Our findings advocate for the adoption of beyond-DFT methodologies, such as sCI, to capture the essential physics of excited-state wavefunctions in strongly correlated materials. The improved accuracy in band gap predictions and the ability to generate more reliable trial wavefunctions for DMC calculations underscore the potential of this approach for broader applications in the study of correlated oxides. This work not only provides a pathway for more accurate simulations of electronic structures in complex materials but also suggests a framework for future investigations of the excited states of other challenging systems.

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

Revealing short- and long-range Li-ion diffusion in Li 2 MnO 3 from finite-temperature dynamical mean field theory

Li 2 MnO 3 is a key component of Li-excess layered cathodes of the form (1 − x), LiMO 2 + x, Li 2 MnO 3 (M = Mn, Ni, Co, …), yet its role in setting Li-ion transport limitations remains under debate. Here, in this study, we combine DFT+U, finite-temperature DFT+DMFT with a continuous-time quantum Monte Carlo impurity solver, and nudged-elastic-band (NEB) calculations to study Li + migration in paramagnetic Li 2 MnO 3 in the presence of a single Li vacancy. Evaluating DMFT total energies along the DFT+U NEB geometries reveals that dynamical correlations strongly renormalize the lowest-barrier processes, reducing the activation energies to E a = 0.18 eV for the shortest-range hop and E a = 0.50 eV for the next-lowest (transport-controlling) step. The 0.18 eV barrier quantitatively reproduces the short-range activation energy from µ+SR, while the 0.50 eV barrier is consistent with the long-range transport scale extracted from ac-impedance measurements. This single-vacancy, paramagnetic DMFT description thus provides a unified interpretation of local and macroscopic probes without invoking clustered vacancy configurations or strong extrinsic disorder, consistent with nearly stoichiometric Li 2 MnO 3 powders. More broadly, our results highlight finite-temperature dynamical correlations as an essential ingredient for predicting ionic migration energetics in correlated oxide electrodes.

Lee, Alex Taekyung [University of Illinois, Chicag↗
DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

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

Structural, magnetic, optical, dielectric and electronic properties of R 2 NiIrO 6 (R = Pr and Nd): A comprehensive experimental and theoretical investigation

Double perovskites are highly promising materials capable of exhibiting a wide variety of phenomena. In this work, we perform a comprehensive experimental and theoretical study of polycrystalline R 2 NiIrO 6 (R = Pr and Nd) compounds. Both compounds were synthesised using the solid-state reaction method. Rietveld refinement confirmed a monoclinic structure with the P2 1 /n space group for both compounds. The scanning electron images showed the average grain sizes of 0.55 μm for R = Pr and 0.46 μm for R = Nd. Fourier transform infrared ra- diation spectra of the two compounds presented two intense bands at 470 cm -1 and 540 cm -1 . The optical measurements revealed that the band gaps of the compounds were in the visible absorption range. The field- cooled magnetisation - field hysteresis measurements indicated exchange bias properties in the synthesised compounds at low temperatures. Both temperature and frequency variation of dielectric constant and loss tangent measurements were conducted. The frequency-dependent ac conductivity measurements indicated that the conductivity increases with the increase of frequency as well as temperature. The Nd 2 NiIrO 6 compound showed lower ac conductivities compared to its isostructural Pr 2 NiIrO 6 compound. The atomic and electronic structures of Nd 2 NiIrO 6 and Pr 2 NiIrO 6 were explored using the spin-polarised calculations performed within the DFT+U method. Our results suggested that the inclusion of on-site correlations and repulsions for the d-states of atoms was necessary in order to obtain finite band gaps of Nd 2 NiIrO 6 and Pr 2 NiIrO 6 systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural Distortions and Uniaxial Negative Thermal Expansion in the Polar Dion–Jacobson Oxide RbNdTa 2 O 7

We provide deeper insight into the crystal structures, sequential structural phase transitions (I2cm → Cmce → I4/mcm → P4/mmm), thermal expansion, and electronic properties of the n = 2 Dion–Jacobson polar oxide RbNdTa 2 O 7 , through X-ray powder diffraction, neutron powder diffraction, Raman studies, and density functional theory calculations. We observed a uniaxial negative thermal expansion (NTE) across the first-order transition, I2cm → Cmce, where the unit cell contracts along the c-axis, which is driven by a contraction of the NdTa 2 O 6 layer. Here, this NTE occurs within the temperature range of the first-order phase transition and contrasts with the corkscrew mechanism typically observed in Ruddlesden–Popper phases. In RbNdTa 2 O 7 , the I2cm (hybrid improper ferroelectric) → Cmce (antipolar) transition involves crucial changes in the bond lengths of Nd and Ta polyhedra, coupled with polar to antipolar displacement of the Nd ions, leading to a net contraction in the NdTa 2 O 6 layer along the c-axis, while preserving the overall octahedral tilting magnitude. This transition highlights the intricate interplay between the Nd and Ta coordination and the associated TaO 6 distortions. Temperature-dependent Raman spectra analysis further confirms the first-order structural transition and associated NTE, providing evidence for increased bond stiffness across this transition. Additionally, using neutron powder diffraction, we have determined that the transition I4/mcm → P4/mmm occurs at approximately 1150 K. Finally, we have calculated from DFT + U, the partial density of states, the energy bandgaps, and effective masses of the charge carriers of the polar ground structure.

Chemical structure↗
DOE OSTI2025

Exploring the Coexistence of Spin States in [Fe(tpy-Ph) 2 ] 2+ Complexes on Au(111) Using DFT Calculations

In this work, we systematically study the electronic structure and stability of spin states of the [Fe-(tpy-ph) 2 ] 2+ molecule in both the gas phase and on a Au(111) substrate using density functional theory + U (DFT+ U ) calculations. We find that the stability of the Fe 2+ ion’s spin states predicted by the computations is significantly influenced by the Hubbard U parameter. In the gas phase, the low-spin (LS, S = 0) state is found to be energetically favorable for U (Fe) ≤ 3 eV, whereas the high-spin (HS, S = 2) state is stabilized for U (Fe) > 3 eV. Interaction with the Au(111) substrate is found to elevate the critical U for the spin-state transition to 3.5 eV. Additionally, we perform L-edge X-ray absorption spectroscopy (XAS) calculations for both HS and LS states. The calculated XAS suggests that the HS state more closely aligns with the experimental observations, indicating the potential coexistence of the HS state as the initial state during the X-ray excitation process. These findings enrich our understanding of spin-state dynamics in [Fe(tpy-Ph) 2 ] 2+ .

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

Excitons in Hematite Fe 2 O 3 : Short-Time Dynamics from TD-DFT and Non-Adiabatic Dynamics Theories

We present a first-principles study of the short-time dynamics of excitons in hematite Fe 2 O 3 . We used time-dependent density functional theory (TD-DFT) with an underlying DFT+U treatment of electron interactions to characterize the electronic structure of excitons and nonadiabatic molecular dynamics theory (NA-MD) to determine their recombination (electronic ground-state recovery) and relaxation dynamics. Decoherence-corrected trajectory surface hopping approaches in NA-MD simulations yielded recovery times of ~1.1 to 1.8 ns and “higher-lying” exciton relaxation times of ~60 to 70 fs, in accord with experimentally derived lifetimes. With hematite phonons in the range of ~100 to 700 cm –1 , higher-lying excitons relax within one or two oscillations of the phonons before getting trapped into an electron–hole pair Exc-3 structure on the first excited state potential energy surface. This structure resembles already a pair of polarons (electron polaron plus hole polaron) with associated lattice distortions three (3) basal planes away. On longer time scales, the electron–hole bipolaronic pair hops to structures Exc-5, then Exc-7, then Exc-9, ... with the electron polaron and hole polaron separated by 5, 7, 9, ... basal planes in a process of charge separation. The largest frequency phonon ~672 cm –1 for the Exc-3 exciton structure is associated with the electron polaron moiety of the exciton. This phonon is a good candidate for giving rise to the recently observed and reported postexcitation transient IR absorption peak.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles Assessment of ZnTe and CdSe as Prospective Tunnel Barriers at the InAs/Al Interface

Majorana zero modes are predicted to emerge in semiconductor/ superconductor interfaces, such as InAs/Al. Majorana modes could be utilized for fault tolerant topological qubits. However, their realization is hindered by materials challenges. The coupling between the superconductor and the semiconductor may be too strong for Majorana modes to emerge, due to effective doping of the semiconductor by the metallic contact. This could be mediated by adding a tunnel barrier of controlled thickness. We use density functional theory (DFT) with Hubbard U corrections, whose values are machine-learned via Bayesian optimization (BO), to assess ZnTe and CdSe as prospective tunnel barriers for the InAs/Al interface. The results of DFT +U(BO) for ZnTe are validated by comparison to angle resolved photoemission spectroscopy (ARPES). We then study bilayer interfaces of the three semiconductors with each other and with Al, as well as trilayer interfaces with a varying number of ZnTe or CdSe layers inserted between InAs and Al. We find that 16 atomic layers of either material completely insulate the InAs from metal induced gap states (MIGS). However, ZnTe and CdSe differ significantly in their band alignment, such that ZnTe forms an effective barrier for electrons, whereas CdSe forms a barrier for holes. Because of Fermi level pinning in the conduction band at the interface, only electron transport is relevant for InAs-based Majorana devices. Therefore, ZnTe is the better choice. Based on the results of our simulations, we suggest conducting experiments with ZnTe barriers in the thickness range of 6–18 atomic layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Importance of enforcing Hund’s rules in density functional theory calculations of rare earth magnetocrystalline anisotropy

Density functional theory (DFT) and its extensions, such as DFT+U and DFT+dynamical mean-field theory, are invaluable for studying magnetic properties in solids. However, rare-earth (R) materials remain challenging due to self-interaction errors and the lack of proper orbital polarization. We show how the orbital dependence of self-interaction error contradicts Hund’s rules and plagues magnetocrystalline anisotropy (MA) calculations, and how analyzing DFT states that respect Hund’s rules can mitigate this issue. We benchmark MA in RCo 5 , R 2 Fe 14 B, and RFe 12 , extending prior work on RMn 6 Sn 6 , achieving excellent agreement with experiments. Additionally, we illustrate a semi-analytical perturbation approach that treats crystal fields as a perturbation in the large spin-orbit coupling limit. Using Gd-4f crystal-field splitting, this method provides a microscopic understanding of MA and enables rapid screening of high-MA materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessment of minnealloy fabrication via three routes

Transformer cores exhibiting higher operating power and improved efficiency are of great interest to electrical utilities, industry, and the de-carbonization effort. Minnealloy, α″-Fe16(C,N)2, a martensite made only of iron, nitrogen, and carbon, has shown the largest saturation magnetization of any soft ferromagnet, 250 emu/g, and tunable magnetocrystalline anisotropy. Given this represents a significant increase in power transferred per cycle compared to legacy transformer core materials, we investigate three novel, industrially scalable routes for fabricating Minnealloy. The martensite phase content is investigated for each route. Vibrating sample magnetometry is used to investigate the change in saturation magnetization and coercivity with respect to the relative content of the desired phase and other iron, iron-nitride, and iron oxide phase impurities. The relationship between structure and magnetic properties of bulk α″-Fe16(C,N)2 is investigated using LDA, PBE, and PBEsol exchange-correlation functionals within the frameworks of Hubbard-corrected density functional theory (DFT+U).

36 MATERIALS SCIENCE↗
DOE OSTI2024

First-principles calculation of Hubbard U for Terbium metal under high pressure

Abstract Using density functional theory (DFT) and linear response approaches, we compute the on-site Hubbard interaction U of elemental Terbium (Tb) metal in the pressure range ∼ 0–65 GPa. The resulting first-principles U values with experimental crystal structures enable us to examine the magnetic properties of Tb using a DFT+U method. The lowest-energy magnetic states in our calculations for different high-pressure Tb phases—including hcp, α -Sm, and dhcp—are found to be compatible with the corresponding magnetic ordering vectors reported in experiments. The result shows that the inclusion of Hubbard U substantially improves the accuracy and efficiency in modeling correlated rare-earth materials. Our study also provides the necessary U information for other quantum many-body techniques to study Tb under extreme pressure conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamical correlations leading to site and orbital selective Mott insulator transition in hydrogen doped SmNiO 3

Electron doping induces metal -to -insulator transition (MIT) in SmNiO 3 as realized by experiments. While earlier density functional theory (DFT) studies with static correlations fell short of explaining the recent MIT observations at lower hydrogen concentrations, we present a comprehensive computational investigation employing an advanced approach. We combine DFT with dynamical mean field theory (DFT + DMFT) to efficiently analyze the insulating behavior of hydrogen -doped SmNiO 3 . In contrast to previous theoretical works, our calculations predict an insulator transition occurring at a reduced doping level of H:Ni = 0.5:1. Specifically, while the DFT + U method reveals a gap opening between p-to-d orbitals, the DMFT approach highlights a gap opening between d-to-d orbitals. Our findings uncover a selective Mott transition in site and orbital characteristics, with the Ni ions proximate to the doped hydrogen exhibiting Mott -like traits. Notably, DMFT calculations highlight a pronounced dependence on Hund's parameter J, implying the presence of Hundness in the Mott insulator. Finally, this study underscores the necessity of accounting for dynamical correlations to accurately describe the electronic structure of strongly correlated electron-doped rare-earth nickelates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Origin of metal-insulator transitions in the parent compounds of ruthenium-pnictide superconductors

Here we study the interplay of the structural phase transition, flat electronic band dispersion, and metal-to-insulator transition (MIT) in the parent compounds of the Ru-pnictide superconductors by using first-principles calculations. Our electron and phonon calculations reveal that Ru(P,As) undergo MIT accompanied by orthorhombic to monoclinic distortion at low temperature, but RuSb stays orthorhombic and metallic in agreement with the experimental findings. It is demonstrated that electronic correlation, as treated in DFT + U, DFT + Gutzwiller, and dynamical mean-field theory (DMFT), cannot induce MIT in the undistorted crystal structure. We find that, although small monoclinic distortion can remove the van Hove singularity at the Fermi level, it does not immediately gap out the Fermi surface and a large value of monoclinic distortion is necessary for a clear MIT, suggesting the possibility of an intermediate pseudogapped monoclinic metallic phase. Furthermore, we predict a light-induced two-step insulator-to-metal and structural transitions in the monoclinic phases of RuP and RuAs, which can be tested in future ultrafast pump-probe experiments as an alternative ideal playground to VO 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Supporting Data for "Constraints on magnetism and correlations in RuO2 from lattice dynamics and Mössbauer spectroscopy"

Contents of this DOI are data for the research paper "Constraints on magnetism and correlations in RuO2 from lattice dynamics and Mossbauer spectroscopy" by the authors: George Yumnam, et. al. The dataset contains data from inelastic x-ray scattering measurements of RuO2 single crystals performed at the Advanced Photon Source in Argonne National Lab. This dataset also contains data from inelastic neutron scattering experiments of RuO2 powder performed at the ARCS spectrometer at Spallation Neutron Source at ORNL. Supporting theoretical calculations based on density functional theory with r2SCAN and DFT+U methods is also included in this dataset. Please see the README.txt files for more information of the dataset and file contents. For comments and questions, please contact: George Yumnam (yumnamg@ornl.gov) and/or Raphael P Hermann (hermannrp@ornl.gov)

density functional theory↗