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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Assessing cathode property prediction via exchange-correlation functionals with and without long-range dispersion corrections

In this work, we benchmark calculated interlayer spacings, average topotactic voltages, thermodynamic stabilities, and band gaps in layered lithium transition-metal oxides (TMOs) and their de-lithiated counterparts, which are used in lithium-ion batteries as positive electrode materials, against available experimental data. Specifically, we examine the accuracy of properties calculated within density functional theory (DFT) using eight different treatments of electron exchange-correlation: the strongly constrained and appropriately normed (SCAN) and Perdew–Burke–Ernzerhof (PBE) density functionals, Hubbard-U-corrected SCAN and PBE (i.e., SCAN+U and PBE+U), and SCAN(+U) and PBE(+U) with added long-range dispersion (D) interactions (i.e., DFT(+U)+D). van der Waals interactions are included respectively via the revised Vydrov-Van Voorhis (rVV10) for SCAN(+U) and the DFT-D3 for PBE(+U). We find that SCAN-based functionals predict larger voltages due to an underestimation of stability of the MO2 systems, while also predicting smaller interlayer spacings compared to their PBE-based counterparts. Furthermore, adding dispersion corrections to PBE has a greater effect on voltage predictions and interlayer spacings than with SCAN, indicating that DFT-SCAN – despite being a ground-state theory – fortuitously captures some short and medium-range dispersion interactions better than PBE. While SCAN-based and PBE-based functionals yield qualitatively similar band gap predictions, there is no significant quantitative improvement of SCAN-based functionals over the corresponding PBE-based versions. Finally, we expect SCAN-based functionals to yield more accurate property predictions than the respective PBE-based functionals for most TMOs, given SCAN's stronger theoretical underpinning and better predictions of systematic trends in interlayer spacings, intercalation voltages, and band gaps obtained in this work.

25 ENERGY STORAGE↗
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 OSTI2021

Magnetic interactions and spin excitations in van der Waals ferromagnet VI 3

Using a combination of density functional theory (DFT) and spin-wave theory methods, we investigate the magnetic interactions and spin excitations in semiconducting VI 3 . Exchange parameters of monolayer, bilayer, and bulk forms are evaluated by mapping the magnetic energies of various spin configurations, calculated using DFT + U, onto the Heisenberg model. The intralayer couplings remain largely unchanged in three forms of VI 3 , while the interlayer couplings show stronger dependence on the dimensionality of the materials. We calculate the spin-wave spectra within a linear spin-wave theory and discuss how various exchange parameters affect the magnon bands. Here, the magnon–magnon interaction is further incorporated, and the Curie temperature is estimated using a self-consistently renormalized spin-wave theory. To understand the roles of constituent atoms on magnetocrystalline anisotropy energy (MAE), we resolve MAE into sublattices and find that a strong negative V-I inter-sublattice contribution is responsible for the relatively small easy-axis MAE in VI 3 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Exploring metastable states in UO 2 using hybrid functionals and dynamical mean field theory

A detailed exploration of the f-atomic orbital occupancy space for UO2 is performed using a first principles approach based on density functional theory (DFT), employing a full hybrid functional within a systematic basis set. Specifically, the PBE0 functional is combined with an occupancy biasing scheme implemented in a wavelet-based algorithm which is adapted to large supercells. The results are compared with previous DFT + U calculations reported in the literature, while dynamical mean field theory is also performed to provide a further base for comparison. This work shows that the computational complexity of the energy landscape of a correlated f-electron oxide is much richer than has previously been demonstrated. The resulting calculations provide evidence of the existence of multiple previously unexplored metastable electronic states of UO 2 , including those with energies which are lower than previously reported ground states.

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 OSTI2022

Hubbard-corrected oxide formation enthalpies without adjustable parameters

A density functional theory (DFT) approach to computing transition metal oxide heat of formation without adjustable parameters is presented. Different degrees of d-electron localization in oxides are treated within the DFT+U approach with site-dependent, first-principles Hubbard U-parameters obtained from linear response theory, and delocalized states in the metallic phases are treated without Hubbard corrections. Comparison of relative stabilities of these differently treated phases is enabled by a local d-electron density matrix-dependent model, which was found by genetic programming against experimental reference formation enthalpies. This mathematically simple model does not explicitly depend on the Hubbard-corrected ionic species and is shown to reproduce the heats of formation of the Mott insulators Ca 2 RuO 4 and Y 2 Ru 2 O 7 within ~3% of experimental results, where the experimental training data did not contain Ru oxides. This newly developed method thus absolves from the need for element-specific corrections fitted to experiments in existing Hubbard-corrected approaches to the prediction of reaction energies of transition metal oxides and metals. The absence of fitting parameters opens up here the possibility to predict relative thermodynamic stabilities and reaction energies involving d-states of varying degree of localization at transition metal oxide interfaces and defects, where site-dependent U-parameters will be particularly important and devising a fitting scheme against experimental data with predictive power would be exceedingly difficult.

transition metal oxides↗
DOE OSTI2022

Theoretical study of the crystal and electronic properties of α - RuI 3

The material α-RuCl 3 , with a two-dimensional Ru honeycomb sublattice, has attracted considerable attention because it may be a realization of the Kitaev quantum spin liquid. Recently, a new honeycomb material, α-RuCl 3 , was prepared under moderately high pressure, and it is stable under ambient conditions. However, different from α-RuCl 3 , α-RuI 3 was reported to be a paramagnetic metal without long-range magnetic order down to 0.35 K. Here, the structural and electronic properties of the quasi-two-dimensional α-RuI 3 are theoretically studied. First, based on first-principles density functional theory calculations, the ABC stacking honeycomb-layer R$\bar{3}$ (No. 148) structure is found to be the most likely stacking order for α-RuI 3 along the c axis. Furthermore, both R$\bar{3}$ and P$\bar{3}$1c are dynamically stable because no imaginary frequency modes were obtained in the phononic dispersion spectrum without Hubbard U. Moreover, the different physical behavior of α-RuI 3 compared to α-RuCl 3 can be understood naturally. The strong hybridization between Ru 4d and I 5p orbitals decreases the “effective” atomic Hubbard repulsion, leading the electrons of RuI3 to be less localized than in RuCl 3 . As a consequence, the effective electronic correlation is reduced from Cl to I, leading to the metallic nature of α-RuI 3 . Based on the DFT+U (Ueff=2 eV) plus spin-orbital coupling, we obtained a spin-orbit Mott insulating behavior for α-RuCl 3 and, with the same procedure, a metallic behavior for α-RuI 3 , in good agreement with experimental results. Furthermore, when introducing large (unrealistic) U eff =6 eV, the spin-orbit Mott gap opens in α-RuI 3 as well, supporting the physical picture we are proposing. Our results provide guidance to experimentalists and theorists working on two-dimensional transition metal tri-iodide layered materials.

2-dimensional systems↗
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 OSTI2021

Ab Initio Study of Atomic and Electronic Structure of Promising Ba4XMn3O12 (X = Nb, Ce, Pr) Oxides for Solar Thermochemical Hydrogen Production

The two-step metal oxide water-splitting cycle is one of the most viable approach for Solar Thermochemical Hydrogen (STCH) production. Challenges exist in finding suitable oxides that can satisfy thermodynamics of the STCH redox cycle under viable range of temperatures and partial pressures. Recently, BXM family members, i.e., Ba4NbMn3O12 (BNM), Ba4CeMn3O12 (BCM), Ba4PrMn3O12 (BPM) are discovered to exhibit promising STCH performance. However, the magnetic degrees of freedom of their experimental crystal structures is not characterized, and only hypothetical models exist for their electronic structure. We performed Monte-Carlo sampling for the magnetic spin arrangement of Mn atoms to provide explicit atomic structure models for their ground state 12R polytype and determined the most stable spin configuration for each member after structure relaxation using density functional theory (DFT) calculations. We also elucidate the most stable charge configuration (Mn3+/Mn4+ arrangement) for BNM among the existing hypothesis in literature. We investigate these strongly correlated oxides using different levels of theory (DFT+U, Hybrid-DFT) to develop better understanding of their electronic structure. These calculations allowed us to develop a model based on molecular orbital interactions that help explain the observed changes in their experimental x-ray spectroscopy measurements. Thereby, help link differences in electronic structure of BXM family members to their varying water-splitting behavior. These structure models will further facilitate future in-depth defect studies to model STCH redox cycle of these oxides.

atomic and magnetic structure↗
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↗
DOE OSTI2021

Enhanced orbital anisotropy through the proximity to a SrTiO 3 layer in the perovskite iridate superlattices

We have used angle-dependent soft x-ray absorption spectroscopy (XAS) at the O K edge and first-principles calculations to investigate the electronic structures of iridate-based superlattices ( SrIr O 3 ) m / ( SrTi O 3 ) ( m = 1 , 2, 3, and ∞). We focus on the pre-edge Ir 5 d t 2 g – O 2 p orbital hybridization feature in the XAS spectra. By varying the measurement geometry relative to the incident photon polarization, we are able to extract the dichroic contrast and observe the systematic increase in the anisotropy of Ir 5 d orbitals as m decreases. First-principles calculations elucidate the orbital anisotropy coming mainly from the enhanced out-of-plane compression of Ir O 6 octahedra in the SrIr O 3 layers that are adjacent to the inserted SrTi O 3 layers. As m decreases, the increased volume fraction of these interfacial SrIr O 3 layers and their contact with the SrTi O 3 layers within the ( SrIr O 3 ) m / ( SrTi O 3 ) supercell lead to enhanced orbital anisotropy. Furthermore, the tilt and rotation of Ir O 6 octahedra are shown to be essential to understand the subtle orbital anisotropy in these superlattices, and constraining these degrees of freedom will give an incorrect trend. Overall, our results demonstrate that the structural constraint from the inserted SrTi O 3 layers, in addition to other electronic means such as polar interface and charge transfer, can serve as a knob to control the orbital degree of freedom in these iridate-based superlattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Dzyaloshinskii-Moriya interaction induced magnetoelectric coupling in a tetrahedral molecular spin-frustrated system

We have investigated magnetoelectric (ME) coupling in the single-molecule magnet Mn 4 Te 4 (PEt 3 ) 4 with tetrahedral spin frustration. Our density functional studies found that an electric dipole moment can emerge with various noncollinear spin orderings. The forms of spin-dependent dipole are determined and consistent with that in noncentrosymmetric magnets driven by the Dzyaloshinskii-Moriya interaction. Writing a parameterized spin Hamiltonian, after solving for eigenvalues and eigenstates, we quantified the ME coupling by calculating the thermal average of the electric and magnetic susceptibilities, which can be influenced by external magnetic and electric fields, respectively. The quadratic relations are expected to be observable in experiments.

36 MATERIALS SCIENCE↗