Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory (DFT)”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Universal machine learning framework for defect predictions in zinc blende semiconductors

Our article introduces a universal predictive framework for point defect formation energies and charge transition levels in a wide chemical space of zinc blende semiconductors and possible impurity atoms selected from across the periodic table. This framework was developed by leveraging high-throughput quantum mechanical simulations benchmarked using some experimental data from the literature, as well as machine learning (ML)-based regressions techniques that map unique materials descriptors to computed defect properties and yield optimized and generalizable models. Furthermore, the power and utility of these models is revealed through quick predictions for thousands of new defects and screening of low-energy impurities, which may tune the equilibrium conductivity in the semiconductor. This work presents, to our knowledge, the largest density functional theory (DFT) dataset of defect properties in semiconductors and the largest DFT+ML-based screening of point defects in semiconductors to date.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Existence of La-site antisite defects in LaMO 3 (M = Mn, Fe, and Co) predicted with many-body diffusion quantum Monte Carlo

The properties of LaMO 3 (M: 3d transition metal) perovskite crystals are significantly dependent on point defects, whether introduced accidentally or intentionally. The most studied defects in La-based perovskites are the oxygen vacancies and doping impurities on the La and M sites. Here, we identify that intrinsic antisite defects, the replacement of La by the transition metal, M, can be formed under M-rich and O-poor growth conditions, based on results of an accurate many-body ab initio approach. Our fixed-node diffusion Monte Carlo (FNDMC) calculations of LaMO 3 (M = Mn, Fe, and Co) find that such antisite defects can have low formation energies and are magnetized. Complementary density functional theory (DFT)-based calculations show that Mn antisite defects in LaMnO 3 may cause the p-type electronic conductivity. These features could affect spintronics, redox catalysis, and other broad applications. Our bulk validation studies establish that FNDMC reproduces the antiferromagnetic state of LaMnO 3 , whereas DFT with PBE (Perdew–Burke–Ernzerhof), SCAN (strongly constrained and appropriately normed), and the LDA+U (local density approximation with Coulomb U) functionals all favor ferromagnetic states, at variance with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Enhancing high-fidelity neural network potentials through low-fidelity sampling

The efficacy of neural network potentials (NNPs) critically depends on the quality of the configurational datasets used for training. Prior research using empirical potentials has shown that well-selected liquid–solid transitional configurations of a metallic system can be translated to other metallic systems. This study demonstrates that such validated configurations can be relabeled using density functional theory (DFT) calculations, thereby enhancing the development of high-fidelity NNPs. Training strategies and sampling approaches are efficiently assessed using empirical potentials and subsequently relabeled via DFT in a highly parallelized fashion for high-fidelity NNP training. Our results reveal that relying solely on energy and force for NNP training is inadequate to prevent overfitting, highlighting the necessity of incorporating stress terms into the loss functions. To optimize training involving force and stress terms, we propose employing transfer learning to fine-tune the weights, ensuring that the potential surface is smooth for these quantities composed of energy derivatives. This approach markedly improves the accuracy of elastic constants derived from simulations in both empirical potential-based NNPs and relabeled DFT-based NNPs. Overall, this study offers significant insights into leveraging empirical potentials to expedite the development of reliable and robust NNPs at the DFT level.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High-throughput determination of Hubbard $U$ and Hund $J$ values for transition metal oxides via the linear response formalism

DFT+U provides a convenient, cost-effective correction for the self-interaction error (SIE) that arises when describing correlated electronic states using conventional approximate density functional theory (DFT). The success of a DFT+U(+J) calculation hinges on the accurate determination of its Hubbard U and Hund J parameters, and the linear response (LR) methodology has proven to be computationally effective and accurate for calculating these parameters. This study provides a high-throughput computational analysis of the U and J values for transition metal d-electron states in a representative set of over 1000 magnetic transition metal oxides (TMOs), providing a frame of reference for researchers who use DFT+U to study transition metal oxides. In order to perform this high-throughput study, an ATOMATE workflow is developed for calculating U and J values automatically on massively parallel supercomputing architectures. Here, to demonstrate an application of this workflow, the spin-canting magnetic structure and unit cell parameters of the multiferroic olivine LiNiPO4 are calculated using the computed Hubbard U and Hund J values for Ni-d and O-p states, and are compared with experiment. Both the Ni-d U and J corrections have a strong effect on the Ni-moment canting angle. Additionally, including a O-pU value results in a significantly improved agreement between the computed lattice parameters and experiment

36 MATERIALS SCIENCE↗
DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2023

Assessing the source of error in the Thomas–Fermi–von Weizsäcker density functional

We investigate the source of error in the Thomas–Fermi–von Weizsäcker (TFW) density functional relative to Kohn–Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn–Sham density, the corresponding energy deviates significantly from the Kohn–Sham value. Here, we show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn–Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn–Sham solution.

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

Lattice Structure and Dynamics of Sparse Molecular Crystals: OsO 4 and RuO 4

OsO 4 and RuO 4 are molecular oxides with unique tetrameric structures and rare +8 oxidation states. Accurately modeling their properties remains challenging for density functional theory (DFT) due to weak intertetramer interactions, which standard functionals fail to capture. Here, in this work, we show that the van der Waals (vdW)-corrected density functional (vdW-DF-optB86b) provides structural parameters that are much closer to experimental values than the standard generalized gradient approximation, with volume predictions that fall within the experimentally observed range. Phonon band structure analysis shows that the inclusion of vdW interactions stabilizes soft phonon modes, highlighting the importance of dispersion corrections for accurate predictions of lattice dynamics. Experimental measurements of the phonon density of states for OsO 4 , obtained via inelastic neutron scattering, demonstrate good agreement with our vdW-DF-optB86b calculations. These results validate OsO 4 and RuO 4 as valuable benchmarks for structural and vibrational calculations via vdW-corrected DFT methods and offer insights for studying the broader class of sparse molecular materials.

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

Transition Metal Taggants in UO 2 from First Principles

The incorporation of transition metals into nuclear fuel has gained attention both to improve fuel properties and as a possible nuclear forensics tool. Recent experimental studies by Ulrich et al. and Adorno Lopes et al. have investigated Ni and Fe as candidate transition metal dopants for potential nuclear forensics purposes and found that there is minimal alteration to key UO 2 fuel properties. In the present work, we performed density functional theory (DFT) investigations into possible defect structures for Ni and Fe incorporation. The dynamic stability of these defect structures was validated by calculating the phonon density of states. We found that transition metal incorporation likely occurs via substitution at U sites in the fluorite crystal structure with a nearby O vacancy for charge-balancing, which agrees with experimentally proposed structures. Additionally, this defect structure does not cause long-range alterations to the crystal parameters—an important consideration for use as a nuclear fuel taggant. Future needed work involves computational investigation of additional defect concentrations using larger supercells, additional transition metal charge states, and thermal effects. Such investigations can be carried forward into sintering models for a more complete understanding of the suitability of Ni and Fe as fuel taggants.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Building materials genome from ground‐state configuration to engineering advance

Individual phases are commonly considered as the building blocks of materials. However, the accurate theoretical prediction of properties of individual phases remains elusive. The top-down approach by decoding genomic building blocks of individual phases from experimental observations is nonunique. The density functional theory (DFT), as a state-of-the-art solution of quantum mechanics, prescribes the existence of a ground-state configuration at 0 K for a given system. It is self-evident that the ground-state configuration alone is insufficient to describe a phase at finite temperatures as symmetry-breaking non-ground-state configurations are excited statistically at temperatures above 0 K. Our multiscale entropy approach (recently terms as Zentropy theory) postulates that the entropy of a phase is composed of the sum of the entropy of each configuration weighted by its probability plus the configurational entropy among all configurations. Consequently, the partition function of each configuration in statistical mechanics needs to be evaluated by its free energy rather than total energy. The combination of the ground-state and symmetry-breaking non-ground-state configurations represents the building blocks of materials and can be used to quantitatively predict free energy of individual phases with the free energy of each configuration predicted from DFT as well as all properties derived from free energy of individual phases.

CALPHAD↗
DOE OSTI2020

High-pressure structural behavior and elastic properties of U3Si5: A combined synchrotron XRD and DFT study

We present an integrated experimental and theoretical study of the structural behavior of U 3 Si 5 at highpressure conditions using angle-dispersive synchrotron X-ray diffraction (XRD) in a diamond anvil cell (DAC) and density functional theory (DFT) calculations. On increasing pressure, the ambient hexagonal structure of U 3 Si 5 with space group P6/mmm remains stable up to 16.7 GPa, the maximum pressure tested with DAC. The bulk modulus and the a- and c-axial moduli of U 3 Si 5 were experimentally determined to be 126 ± 4 GPa, 173 ± 8 GPa and 79.7 ± 4.3 GPa, respectively. Thus an anisotropy in the axial compressibility of U 3 Si 5 is observed with its c-axis being more compressible than the a-axis. Our DFT calculation results are in general agreement with the experimental values, including reproducing the compressibility anisotropy. A comparison of the bulk modulus of U 3 Si 5 to those of other U-Si compounds reveals a general trend that the bulk modulus of U-Si decreases with increasing U/(U+Si) ratio.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tests of the DFT Ladder for the Fulminic Acid Challenge

Properties of the historically pivotal fulminic acid (HCNO) molecule have been computed with a panoply of 473 density functionals of all varieties, providing a snapshot of the performance of contemporary density functional theory (DFT) for a challenging chemical system. Exhaustive tabulations and statistical analyses have been carried out for geometric parameters, vibrational frequencies, barriers to linearity, and the HCN–O dissociation energy. As the DFT ladder is climbed, confusion rather than consensus ensues regarding the details of the distinctive, extremely flat H–C–N bending potential of fulminic acid and whether the equilibrium structure is linear or bent. While high-ranking DFT functionals produce the smallest errors for the HCN + O( 3 P) → HCNO reaction energy, lower rungs emerge as the best performers for many of the bond distances and harmonic vibrational frequencies. This research shows that the current DFT zoo of approximations does not constitute a transparent ladder of increasingly accurate methods that consistently converges on definitive predictions for various properties of HCNO. Additional analyses are performed on the side effects of popular dispersion corrections on the covalently bonded properties and thermochemistry of HCNO.

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

Fluorinated High-Voltage Electrolytes To Stabilize Nickel-Rich Lithium Batteries

As state-of-the-art (SOA) lithium-ion (Li-ion) batteries approach their specific energy limit (~250 Wh kg -1 ), layer-structured, nickel-rich (Ni-rich) lithium transition metal oxide-based cathode materials, e. g., LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811), have attracted great interest owing to their practical high specific capacities (>200 mAhg -1 ). Coupled with their high average discharge voltages (~4 V vs Li/Li + ), Ni-rich cathode-based lithium batteries possess a great potential to achieve much higher specific energies (>350 Wh kg -1 at the cell level) than the SOA Li-ion counterparts. In addition, Ni-rich oxides are low-cost battery cathode materials due to their low cobalt contents. However, Ni-rich cathode-based lithium batteries suffer quick capacity degradations upon cycling, particularly at high upper cutoff voltages (e.g., ≥= 4.5 V vs Li/Li + ), due to crystal structure changes of the active cathode materials and parasitic side reactions at the electrolyte/electrode interfaces. In this study, a fluorinated-solvent-based, high-voltage stable electrolyte (HVE), i.e., 1 M Li bis(trifluoromethanesulfonyl)imide (LiTFSI) in fluoroethylene carbonate (FEC), bis(2,2,2-trifluoroethyl) carbonate (FDEC), and methyl (2,2,2-trifluoroethyl) carbonate (FEMC) with Li difluoro(oxalate)borate (LiDFOB) additive, was formulated and evaluated in Li/NMC811 battery cells. Further, to the best of our knowledge, this class of electrolyte has not been investigated for Ni-rich cathode-based lithium batteries. Li/NMC811 cells with HVE exhibited a superior long-term cycle performance stability, maintaining ~80% capacity after ~500 cycles at a high cutoff voltage of 4.5 V (vs Li/Li + ) than a baseline carbonate-solvent-based electrolyte (BE). The superior cycle stability of the Li/NMC811 cells is attributed to the inherently high-voltage stability of HVE, supported by the physical and electrochemical analyses. This conclusion is supported by our density functional theory (DFT) modeling where HVE shows a less tendency of deprotonation/oxidation than BE, leading to the observed cycle stability. The findings in this study are important to help tackle the technical challenges facing Ni-rich cathode-based lithium batteries to realize their high energy density potentials with a long cycle life.

25 ENERGY STORAGE↗
DOE OSTI2020

Phonon Spectroscopy in Antimony and Tellurium Oxides

α-Sb 2 O 3 (senarmontite), β-Sb 2 O 3 (valentinite), and α-TeO 2 (paratellurite) are compounds with pronounced stereochemically active Sb and Te lone pairs. The vibrational and lattice properties of each have been previously studied but often lead to incomplete or unreliable results due to modes being inactive in infrared or Raman spectroscopy. Here in this paper, we present a study of the relationship between bonding and lattice dynamics of these compounds. Mössbauer spectroscopy is used to study the structure of Sb in α-Sb 2 O 3 and β-Sb 2 O 3 , whereas the vibrational modes of Sb and Te for each oxide are investigated using nuclear inelastic scattering, and further information on O vibrational modes is obtained using inelastic neutron scattering. Additionally, vibrational frequencies obtained by density functional theory (DFT) calculations are compared with experimental results in order to assess the validity of the utilized functional. Good agreement was found between DFT-calculated and experimental density of phonon states with a 7% scaling factor. The Sb–O–Sb wagging mode of α-Sb 2 O 3 whose frequency was not clear in most previous studies is experimentally observed for the first time at ~340 cm -1 . Softer lattice vibrational modes occur in orthorhombic β-Sb 2 O 3 compared to cubic α-Sb 2 O 3 , indicating that the antimony bonds are weakened upon transforming from the molecular α phase to the layer-chained β structure. The resulting vibrational entropy increase of 0.45 ± 0.1 kB/Sb 2 O 3 at 880 K accounts for about half of the α–β transition entropy. The comparison of experimental and theoretical approaches presented here provides a detailed picture of the lattice dynamics in these oxides beyond the zone center and shows that the accuracy of DFT is sufficient for future calculations of similar material structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2022

High-precision predictions of properties of chemically disordered crystals

Multiple scattering theory (MST) combined with density functional theory (DFT) allows to predict properties of chemically disordered materials from the first principles. However, such predictions often suffer from the systematic errors, which depend on crystal geometry. Each computed property of a particular crystal structure typically has a relatively small random error and a larger systematic error. Cancellation of systematic errors allows more accurate predictions. We propose a computational methodology based on the subtraction of the systematic errors in MST. To exemplify it, we apply it to the precipitated alloys. Considering precipitation strengthening in Ni superalloys, we compute the relative enthalpies of the competing Ni_3(Al_{1-x}Ti_x)_1 crystal structures with a chemical disorder on the Al+Ti sublattice. Such predicted composition-structure-property dependencies are useful for the guided design of the next-generation alloys with improved strength. Our predictions are validated by comparison with the results of other DFT methods (having a higher computational cost) and with experiment.

density functional theory↗
DOE OSTI2021

Thermodynamic properties of the Nd-Bi system via emf measurements, $\mathrm{DFT}$ calculations, machine learning, and $\mathrm{CALPHAD}$ modeling

Thermodynamic properties of the Nd-Bi system were investigated using a combination of experimental measurements, first-principles calculations based on density functional theory (DFT), data mining and machine learning (DM + ML) predictions, and calculation of phase diagrams (CALPHAD) modeling. The electromotive force (emf) of Nd-Bi alloys in molten LiCl-KCl-NdCl 3 at 773–973 K was measured via coulometric titration of Nd into Bi for the determination of thermochemical properties such as activity coefficients and solubilities of Nd in Bi. A new peritectic reaction of [liquid + NdBi 2 = Nd 3 Bi 7 ] at 774 K was confirmed using differential scanning calorimetry, structural (X-ray diffraction), and microstructural (scanning electron microscopy) analyses. The unknown crystal structure of NdBi2 was suggested to be a mixture of the anti-La 2 Sb configuration and the La 2 Te-type configuration based on ML predictions for over 26,000 data-mined AB 2 -type configurations together with DFT-based verifications. Using the newly acquired experimental data and DFT-based calculations, the thermodynamic description of the Nd-Bi system was remodeled, and a more complete Nd-Bi phase diagram was calculated, including the Nd 3 Bi 7 compound, invariant transition reactions, and liquidus temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanistic Study of 1,2-Dichloroethane Hydrodechlorination on Cu-Rich Pt–Cu Alloys: Combining Reaction Kinetics Experiments with DFT Calculations and Microkinetic Modeling

Cu-rich Pt–Cu bimetallic catalysts are among the most promising candidates for actively catalyzing the hydrodechlorination of 1,2-dichloroethane (1,2-DCA) toward ethylene production. Combining reaction kinetics experiments with density functional theory (DFT) calculations and mean-field microkinetic modeling, we present a systematic mechanistic study for 1,2-DCA hydrodechlorination on Cu-rich Pt–Cu alloy catalysts. Our DFT (PBE+(TS+SCS)) results suggest that increasing Cu content in the Pt–Cu alloy destabilizes C 2 -species adsorption while stabilizing the binding of atomic chlorine. The reaction energetics of all the elementary steps in the 1,2-DCA reaction network were calculated on a Pt 1 Cu 3 (111) model surface. The DFT results were then used to construct a microkinetic model, and the model-predicted reaction rates were compared with our reaction kinetics experimental results on a Cu-rich SiO 2 -supported Pt–Cu alloy catalyst through a parameter estimation procedure. Both the reaction kinetics experiments and the microkinetic model after parameter adjustments yielded 100% selectivity to ethylene. The microkinetic model pointed to a reaction pathway involving two sequential chlorine-removal steps on the Pt–Cu alloy catalyst, a mechanism distinct from the one previously identified on pure Pt/SiO 2 catalysts, which involved an initial hydrogen-removal step. Adjustments to the DFT-derived parameters indicate the possible formation of chlorine-induced Cu-enriched surface sites during 1,2-DCA hydrodechlorination conditions, sites that are more active than those encountered in the bulk Pt 1 Cu 3 (111) alloy surface. Furthermore, our study offers valuable initial insights on the 1,2-DCA hydrodechlorination reaction mechanism and the nature of the active sites on PtCu bimetallic catalysts.

1,2-Dichloroethane↗
DOE OSTI2022

Hydrogen trapping energetics at BCC iron-helium interfaces

Density functional theory (DFT) calculations have been performed to assess the trapping and segregation strength of hydrogen (H) to helium (He) interfaces in BCC (Fe) as a function of He density and surface orientation. The He density ranges from 0 to 2 He/V which corresponds to equilibrium bubbles that are expected to form in the structural components in a fusion reactor environment. The DFT calculations consist of a slab of BCC metal and an initial lattice of FCC He. Here, the report on the binding energy of H to these interfaces as well as the migration barriers into and away from the surface, which together provide information on the de-trapping energy for H at He bubbles. It was found that the binding energy of H to He-Fe interfaces decreases with increasing He density. Migration barriers are also modified due to the He density, which are sensitive to the surface orientation as well as the diffusion pathway. These results provide valuable first-principles energetics that are necessary for higher-order, mesoscale models, and provide insight into the extent to which He bubbles may trap tritium in fusion structural materials.

42 ENGINEERING↗