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

In Situ High-Temperature Raman Spectroscopy of UCl 3 : A Combined Experimental and Theoretical Study

Uranium trichloride (UCl 3 ) has received growing interest for its use in uranium-fueled molten salt reactors and in the pyrochemical processing of used fuel. In this paper, we report for the first time the experimentally determined Raman spectra of UCl 3 , at both ambient condition and in situ high temperatures up to 871 K. The frequencies of five of the Raman-active vibrational modes (v i ) of UCl 3 exhibit a negative temperature derivative ((∂ν i /∂T) P ) with increasing temperature. Here, this red-shift behavior is likely due to the elongation of U–Cl bonds. The average isobaric mode Grüneisen parameter (γ iP = 0.91 ± 0.02) of UCl 3 was determined through use of the coefficient of thermal expansion published in Vogel et al. (2021) and the (∂ νi /∂T) P values determined in this study. These results are in general agreement with those calculated here by density functional theory (DFT+U). Finally, a comparison of the ambient band positions of UCl 3 to those of isostructural lanthanide (La–Eu) and actinide chlorides (Am–Cf) has been made.

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

Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite

Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.

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

Correlation versus hybridization gap in CaMn$$_{2}$$Bi$$_{2}$$

Abstract We study the interplay between electronic correlations and hybridization in the low-energy electronic structure of CaMn $$_2$$ 2 Bi $$_2$$ 2 , a candidate hybridization-gap semiconductor. By employing a DFT+ U approach we find both the antiferromagnetic Néel order and band gap in good agreement with the corresponding experimental values. Under hydrostatic pressure, we find a crossover from hybridization gap to charge-transfer insulting physics due to the delicate balance of hybridization and correlations. Increasing the pressure above $$P_c=4$$ P c = 4 GPa we find a simultaneous pressure-induced volume collapse, plane-to-chain, insulator to metal transition. Finally, we have also analyzed the topology in the antiferromagnetic CaMn $$_2$$ 2 Bi $$_2$$ 2 for all pressures studied.

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 OSTI2023

Quantifying uncertainty in high-throughput density functional theory: A comparison of AFLOW, Materials Project, and OQMD

A central challenge in high-throughput density functional theory (HT-DFT) calculations is selecting a combination of input parameters and postprocessing techniques that can be used across all materials classes, while also managing accuracy-cost tradeoffs. To investigate the effects of these parameter choices, we consolidate three large HT-DFT databases: Automatic-FLOW (AFLOW), the Materials Project (MP), and the Open Quantum Materials Database (OQMD), and compare reported properties across each pair of databases for materials calculated using the same initial crystal structure. We find that HT-DFT formation energies and volumes are generally more reproducible than band gaps and total magnetizations; for instance, a notable fraction of records disagree on whether a material is metallic (up to 7%) or magnetic (up to 15%). The variance between calculated properties is as high as 0.105 eV/atom (median relative absolute difference, or MRAD, of 6%) for formation energy, 0.65 Å 3 / atom (MRAD of 4%) for volume, 0.21 eV (MRAD of 9%) for band gap, and 0.15 μB/ formula unit (MRAD of 8%) for total magnetization, comparable to the differences between DFT and experiment. Here, we trace some of the larger discrepancies to choices involving pseudopotentials, the DFT + U formalism, and elemental reference states, and argue that further standardization of HT-DFT would be beneficial to reproducibility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accuracy of DFT computed oxygen-vacancy formation energies and high-throughput search of solar thermochemical water-splitting compounds

The enthalpy change involved in metal oxide reduction is a key quantity in various processes related to energy conversion and storage, and is of particular interest for computational prediction. Often this prediction involves the simulation of a high temperature reduction process with a 0K methodology like density functional theory (DFT), and it is not infrequent for the high temperature and 0K stable crystal structures to differ. This introduces a conundrum with regards to the choice of crystal structure to utilize in the computation, with approaches in the literature varying and experimental validation remaining scarce. In this work we address both the crystal structure conundrum and the experimental validation, and then apply the insights we gain to guide a high-throughput search for new materials for solar thermochemical water-splitting applications. By computing the DFT+U oxygen vacancy formation energy (ΔE vf ) of a selection of ABO 3 compounds and comparing different crystal structures for each composition, we highlight the issues that arise when the structure utilized in the computation is dynamically unstable at 0K, namely the presence of an artificial lowering of ΔE vf , and the lack of convergence of ΔE vf with cell size. We solve these limitations by identifying and employing a suitable surrogate dynamically stable structure. We then validate the predictive power of our calculations against appositely generated experimental measurements of reduction enthalpy for a series of Hubbard U values, finding an accuracy ranging between 0.2-0.6 eV/O. In light of such conclusions, we revise and expand a previous a high-throughput DFT study on ABO 3 perovskite oxides. As a result, we provide a list of candidate STCH materials, highlight trends with redox-active cation and structural distortion, and identify Mn 4+ , Mn 3+ and Co 3+ as the most promising redox-active cations.

08 HYDROGEN↗
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↗
DOE OSTI2025

Selected magnetoelectrical phonon stiffening in Cr 2 ⁢O 3

A thorough understanding of lattice dynamics, particularly how they couple with spin through spin-phonon interactions, is crucial for pioneering new spin-caloritronic applications. Here, despite extensive studies on Cr 2 ⁢O 3 , the origin of nonlinear thermal expansions and the unusual stiffening of optical phonon modes is still elusive. Through inelastic neutron scattering and atomistic calculations, we found that these behaviors can be attributed to the renormalization of electron states owing to the magnetic transition, and we ruled out the effects from thermal expansion, phonon anharmonicity, magnetostriction, or electron-phonon interactions. Importantly, our quantitative modeling suggests that the common belief that dynamic spin-phonon interactions are the origin of anomalous phonon energy stiffening in magnetic compounds does not apply in Cr 2 ⁢O 3 .

DFT+U↗
DOE OSTI2025

High-resolution neutron diffraction determination of noncollinear antiferromagnetic order in the honeycomb magnetoelectric Fe 4 ⁢Nb 2⁢ O 9

Magnetoelectric systems offer potential for device applications exploiting coupled states between electric and magnetic properties. Among magnetoelectric materials, Fe 4 ⁢Nb 2 ⁢O 9 has attracted special attention because of its pronounced dielectric signal at high magnetic transition temperatures. However, the magnetic ground state, which is essential information for understanding its unusual magnetoelectricity, remains unclarified. Here, we report a noncollinear magnetic ground state of Fe 4 ⁢Nb 2 ⁢O 9 . To examine the magnetoelectric effect associated with sequential magnetic and structural transitions upon cooling, we conducted combined x-ray diffraction, magnetic susceptibility, magnetization, dielectric constant, and magnetodielectric experiments. Powder neutron diffraction experiments revealed a series of magnetic Bragg peaks and clear splitting of peaks via structural transition. Magnetic Rietveld refinements, combined with group theory analysis, determined a noncollinear antiferromagnetic structure including a significant 𝑐-axis moment component at 1.5 K. This article provides insights into the understanding of its magnetoelectric properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗