Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT Calculations”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
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 OSTI2020

Application of Self-Interaction Corrected Density Functional Theory to Early, Middle, and Late Transition States

Density functional theory (DFT)-based methods often significantly underpredict chemical reaction barriers compared with experiments because of the tendency of DFT to over-stabilize transition states with stretched bonds due to the impact of unphysical electron self-interaction. But many reactions have early or late transition states where the transition state geometry closely resembles the reactants or products, respectively. The role of self-interaction in those cases is not known. Here we compare the performance of DFT with and without self-interaction correction (SIC) for describing the hydrogenation of CO and CO 2 catalyzed by a Lewis acid-base pair incorporated onto an aromatic cluster, using CCSD(T) results for reference. The three elementary steps in these reactions consist of an early, a middle, and a late transition. Our results show that the Perdew-Zunger SIC (PZ-SIC), implemented in the Fermi-L Lowdin orbital SIC (FLO-SIC) approach, qualitatively improves the description of the forward and reverse reaction barriers relative to uncorrected DFT for the middle transition, but not the early or late transitions. By contrast, the local scaling SIC (LSIC) method, also implemented in the FLO-SIC framework, significantly improves the calculated barriers over DFT and PZ-SIC in all but one case. Lastly, the results also show how the FLO-SIC approach can provide insight into the bonding in aromatic systems.

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

Antimony stable isotope fractionation during adsorption onto birnessite: A molecular perspective from X-ray absorption spectroscopy and density functional theory

Sorption of antimony (Sb) onto birnessite significantly influences the fate of Sb in oceanic and terrestrial environments and fractionates Sb isotopes. Nevertheless, little is known about Sb isotopic fractionation during its adsorption on birnessite. Here, in this study, we show the value of Δ 123 Sb adsorbed-aqueous increases from −0.398 to −0.332 ‰ in 1 h and then decreases and stabilizes at −0.384 ‰ in 72 h. The enrichment of the light Sb isotope is predominantly due to the distortion of the octahedral symmetry. X-ray absorption spectroscopy results indicate Sb first forms a double-corner-sharing complex on birnessite and then transforms to a double-edge-sharing complex during adsorption. The optimized bond distances for double-corner-sharing (3.37 Å) and double-edge-sharing (2.90 Å) complexes calculated using density functional theory (DFT) fits well with the structure (3.41 and 3.00 Å) revealed by X-ray absorption spectroscopy, respectively. The fractionation derived from reduced partition function ratios calculated using DFT aligns well with the experimental results. Therefore, the variation in Sb isotopic fractionation during adsorption is attributed to the evolving structure of Sb complexes on birnessite. Our results demonstrate the isotopic fractionation of Sb during adsorption on birnessite and provide a molecular-scale understanding of Sb behavior, contributing to the correct reconstruction of the Sb isotope composition of ancient seawater using ferromanganese crusts and nodules, and efforts to trace Sb migration in epigenetic mining environments.

Adsorption↗
DOE OSTI2026

Accurate Electron-Phonon Interactions from Advanced Density Functional Theory

Electron-phonon coupling (EPC) is key to understanding charge transport, band renormalization, and superconductivity in energy materials, including correlated transition-metal oxides, ferroelectric perovskites, optoelectronic semiconductors, and phonon-mediated superconductors. Although first-principles density functional theory (DFT)-based EPC calculations are used widely, their predictive power is limited by the accuracy, transferability, and efficiency of the underlying exchange-correlation functionals. These limitations become exacerbated in complex 𝑑- and 𝑓-electron materials, where beyond-DFT approaches and additional corrections, such as the Hubbard 𝑈, are commonly invoked. Here, using the examples of Co⁢O and Ni⁢O, we show how the 𝑟 2 SCAN density functional correctly captures strong EPC effects in transition-metal oxides without requiring the introduction of Hubbard 𝑈 parameters. We also find that 𝑟 2 SCAN successfully describes the subtle interplay between ionic and covalent bonding, and strong EPC effects in the low-temperature rhombohedral phase of the prototypical ferroelectric Ba⁢Ti⁢O 3 , without requiring Hubbard 𝑈 and intersite 𝑉 corrections. We further demonstrate the ability of 𝑟 2 SCAN to accurately model the EPC of the main-group semiconductor Ga⁢As and the phonon-mediated superconducting properties of Mg⁢B 2 , with reliable electronic bands and phonons. Our study establishes 𝑟 2 SCAN-based EPC as a transferable, parameter-free, and computationally efficient framework for predictive material-specific modeling of EPC in energy materials, and opens a practical route toward high-throughput screening of superconductors, thermoelectrics, optoelectronic semiconductors, and oxide electrodes and catalysts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

High electron mobility of Al x Ga1− x N evaluated by unfolding the DFT band structure

We calculate the alloy-disorder-limited electron mobility of AlxGa1−xN from first principles. AlxGa1−xN is a technologically important ultra-wide-bandgap alloy with promise in light emitting diodes and high-power transistors. Alloying introduces statistical disorder, which causes electrons to scatter between different crystal-momentum states, leading to a reduction in mobility for intermediate alloy compositions. The corresponding lifetime, which appears as an energy broadening in the band structure, can be evaluated by unfolding the band structure from the supercell basis to the primitive-cell basis. We fit the first-principles band broadening with a model scattering potential and evaluate the low-field electron mobility using the semiclassical Boltzmann transport equation in the relaxation-time approximation. Our calculated mobility is in agreement with experimental values. We also find the lowest alloy-scattering electron mobility (total electron mobility) across the entire composition range to be 186 cm2/V s (136 cm2/V s), which is comparable to the highest electron mobility predicted in the competitor system, β-(AlxGa1−x)2O3. Our results elucidate the intrinsic limits imposed by alloy disorder on electron transport in AlxGa1−xN.

Pant, Nick (ORCID:0000000254506045)↗
DOE OSTI2020

Accelerating Finite-temperature Kohn-Sham Density Functional Theory\ with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accelerating Finite-Temperature Kohn-Sham Density Functional Theory with Deep Neural Networks

We present a numerical modeling workflow based on machine learning (ML) which reproduces the total energies produced by Kohn-Sham density functional theory (DFT) at finite electronic temperature to within chemical accuracy at negligible computational cost. Based on deep neural networks, our workflow yields the local density of states (LDOS) for a given atomic configuration. From the LDOS, spatially-resolved, energy-resolved, and integrated quantities can be calculated, including the DFT total free energy, which serves as the Born-Oppenheimer potential energy surface for the atoms. We demonstrate the efficacy of this approach for both solid and liquid metals and compare results between independent and unified machine-learning models for solid and liquid aluminum. Our machine-learning density functional theory framework opens up the path towards multiscale materials modeling for matter under ambient and extreme conditions at a computational scale and cost that is unattainable with current algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

CpFe(CO) 2 Radical Generated from Dinuclear [CpFe(CO) 2 ] 2 and Mononuclear (Cp)(CO) 2 Fe(H): Density Functional Theory Is Accurate for One, But Not Both

Density functional theory (DFT) methods remain the most practical approach to calculating properties and reaction mechanisms of transition metal complexes. While the accuracy of DFT methods has been evaluated for some properties of mononuclear organometallic complexes there has been a general lack of evaluation for dinuclear organometallic complexes, in particular bonding changes related to reaction mechanisms. Here, this work evaluated DFT and coupled cluster methods for the accuracy of calculating the CpFe(CO) 2 radical (Fp•) generated from dinuclear [CpFe(CO) 2 ] 2 (Fp 2 ) and mononuclear [(Cp)(CO) 2 Fe(H)] (Fp-H). This transition metal radical fragment was evaluated because dinuclear complexes built with it have recently shown a variety of unique reactions but has proven challenging to accurately calculate with DFT methods. Here we show that DFT methods provide a surprising wide range of fragmentation energies for Fp 2 and lower and mid rung DFT methods as well as DLPNO–CCSD(T) perform well for this dissociation energy. The highest rung double-hybrid methods have a large range in the Fp 2 dissociation energy, and the energy greatly depends on the amount of MP2 correlation energy included. For generating Fp• from Fp-H the lower and mid rung methods that worked well for Fp 2 showed significant error. Double-hybrid methods unfortunately are only accurate for the Fe–H bond if they are very inaccurate for the Fp 2 dissociation energy. While DLPNO–CCSD(T) is not perfect, and not close to chemically accurate for the Fe–H bond, it does provide reasonable accuracy for both Fp 2 and Fp-H dissociation energies.

density functional theory↗
DOE OSTI2025

Accurate point defect energy levels from non-empirical screened range-separated hybrid functionals: The case of native vacancies in ZnO

We use density functional theory (DFT) with non-empirically tuned screened range-separated hybrid (SRSH) functionals to calculate the electronic properties of native zinc and oxygen vacancy point defects in ZnO, and we predict their defect levels for thermal and optical transitions in excellent agreement with available experiments and prior calculations that use empirical hybrid functionals. Furthermore, the ability of this non-empirical first-principles framework to accurately predict quantities of relevance to both bulk- and defect-level spectroscopy enables high-accuracy DFT calculations with non-empirical hybrid functionals for defect physics, at a reduced computational cost.

Defects↗
DOE OSTI2022

Development of efficient solar cells using combination of QSPR and DFT approaches

The awarded research project aimed to design and determine the effective organic dye-sensitizers from the ensemble of chromophores for the dye-sensitized solar cell (DSSC) and fullerene-derivatives (FDs) as acceptor for polymer solar cell (PSC) followed by explore the electron transfer mechanism using the combined quantitative structure-property relationship (QSPR) analysis in conjunction with density functional theory (DFT) and time-dependent DFT (TDDFT)-based calculations. Based on our proposal we have employed different in silico approaches to design photo-efficient organic dye-sensitizers for DSSCs and FDs as acceptor for PSCs with exhaustive screening methodology.

14 SOLAR ENERGY↗
DOE OSTI2025

Relative phase stability of L1 2 and DO 22 /DO 23 structures in Al 3 Nb, Al 3 Zr and Al 3 V compounds

The relative stability of the different tri-aluminide (Al 3 M) phases in three binary systems (M = Zr, Nb and V) was assessed for their potential to form fine cubic L1 2 precipitates in additively manufactured alloys. Supersaturated thin films of Al-(8–30) at% M were sputtered and heat treated during in-situ x-ray diffraction (XRD) measurements to observe the temperature ranges of stability for each phase. As-sputtered films were then processed with laser tracks simulating additive manufacturing solidification conditions, and the formation of phases in the laser tracks was correlated with density functional theory (DFT) and nucleation rate calculations. We found that the metastable L1 2 structure is highly competitive with the stable DO 23 structure in the Al-Zr system, but much less stable than the DO 22 structure in the Al-Nb system, and both the DO 22 and Al 8 V 5 structure in the Al-V system. Furthermore, these experimental results were found to be in good agreement with the DFT and kinetic calculations, as we determined that the metastable L1 2 in Al-Zr only requires a small amount of undercooling to favor its nucleation over the stable DO 23 , suggesting additive manufacturing can be a viable pathway to develop Al-Zr alloys strengthened by a high volume fraction of L1 2 Al 3 Zr phase.

Perrin, Alice E. [Oak Ridge National Laboratory (O↗
DOE OSTI2025

Quick-and-Easy Validation of Protein–Ligand Binding Models Using Fragment-Based Semiempirical Quantum Chemistry

Electronic structure calculations in enzymes converge very slowly with respect to the size of the model region that is described using quantum mechanics (QM), requiring hundreds of atoms to obtain converged results and exhibiting substantial sensitivity (at least in smaller models) to which amino acids are included in the QM region. As such, there is considerable interest in developing automated procedures to construct a QM model region based on well-defined criteria. However, testing such procedures is burdensome due to the cost of large-scale electronic structure calculations. Here, we show that semiempirical methods can be used as alternatives to density functional theory (DFT) to assess convergence in sequences of models generated by various automated protocols. The cost of these convergence tests is reduced even further by means of a many-body expansion. We use this approach to examine convergence (with respect to model size) of protein–ligand binding energies. Fragment-based semiempirical calculations afford well-converged interaction energies in a tiny fraction of the cost required for DFT calculations. Two-body interactions between the ligand and single-residue amino acid fragments afford a low-cost way to construct a “QM-informed” enzyme model of reduced size, furnishing an automatable active-site model-building procedure. This provides a streamlined, user-friendly approach for constructing ligand binding-site models that needs neither a priori information nor manual adjustments. Extension to model-building for thermochemical calculations should be straightforward.

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

Examination of Solvent- Room Temperature Ionic Liquid Interactions using 2D NMR and Computational Chemistry

Rare earth metals (REMs) are the key functioning material found in microelectronics, advanced batteries, and high-efficiency electrical motors. Current methods for producing REMs require high temperatures (i.e. energy intensive) and produce an abundance of hazardous fluoride waste. Room temperature ionic liquids (RTILs) are an ideal alternative electrolyte because they allow for a low-temperature, energy-efficient electrochemical process. RTIL properties such as viscosity and conductivity determine how effective they can be in this process. Since viscosity is directly related to both intra and intermolecular forces, this study sets out to investigate two RTILs-Trihexyltetradecylphosphonium dicyanamide (Cyphos 105) and Trihexyltetradecylphosphonium bistrifimide (Cyphos 109) and their interactions. Results from 2-dimensional nuclear magnetic resonance (2D NMR) spectroscopy show that there is an interaction between the tetradecyl tail and the phosphorus atom of the Cyphos RTIL. In order to establish whether this is an intra or an inter-molecular interaction, DFT-based computational calculations were performed to determine the optimized geometry. Calculations were done both in the gas phase and with propylene carbonate as an augmenting solvent. Preliminary results show that there is minimal intramolecular interaction between the tetradecyl tail and the phosphorus atom of the RTILs, and that the interaction seen in the 2D NMR experiment is likely of an intermolecular nature instead.

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

Solvent manipulation of the pre-reduction metal–ligand complex and particle-ligand binding for controlled synthesis of Pd nanoparticles

Understanding how to control the nucleation and growth rates is crucial for designing nanoparticles with specific sizes and shapes. Herein, we show that the nucleation and growth rates are correlated with the thermodynamics of metal–ligand/solvent binding for the pre-reduction complex and the surface of the nanoparticle, respectively. To obtain these correlations, we measured the nucleation and growth rates by in situ small angle X-ray scattering during the synthesis of colloidal Pd nanoparticles in the presence of trioctylphosphine in solvents of varying coordinating ability. The results show that the nucleation rate decreased, while the growth rate increased in the following order, toluene, piperidine, 3,4-lutidine and pyridine, leading to a large increase in the final nanoparticle size (from 1.4 nm in toluene to 5.0 nm in pyridine). Using density functional theory (DFT), complemented by 31 P nuclear magnetic resonance and X-ray absorption spectroscopy, we calculated the reduction Gibbs free energies of the solvent-dependent dominant pre-reduction complex and the solvent-nanoparticle binding energy. The results indicate that lower nucleation rates originate from solvent coordination which stabilizes the pre-reduction complex and increases its reduction free energy. At the same time, DFT calculations suggest that the solvent coordination affects the effective capping of the surface where stronger binding solvents slow the nanoparticle growth by lowering the number of active sites (not already bound by trioctylphosphine). The findings represent a promising advancement towards understanding the microscopic connection between the metal–ligand thermodynamic interactions and the kinetics of nucleation and growth to control the size of colloidal metal nanoparticles.

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

Thermodynamic modeling of the Pd–Zn system with uncertainty quantification and its implication to tailor catalysts

Pd–Zn intermetallic catalysts show encouraging combinations of activity and selectivity on well-defined active site ensembles. Thermodynamic description of the Pd–Zn system, delineating phase boundaries and enumerating site occupancies within intermediate alloy phases, is essential to determining the ensembles of Pd–Zn atoms as a function of composition and temperature. Combining the present extensive first-principles calculations based on density functional theory (DFT) and available experimental data, the Pd–Zn system was remodeled using the CALculation of PHAse Diagrams (CALPHAD) approach. High throughput modeling tools with uncertainty quantification, i.e., ESPEI and PyCalphad, were incorporated in the phase analysis. Here, the site occupancies across the γ phase composition region were given special attention. A four-sublattice model was used for the γ phase owing to its four Wyckoff positions, i.e., the outer tetrahedral (OT) site 8c, the inner tetrahedral (IT) site 8c, the octahedral (OH) site 12e, and the cuboctahedral (CO) site 24g. The site fractions of Pd and Zn calculated from the present thermodynamic model show the occupancy preference of Pd in the OT and OH sublattices in agreement with experimental observations. The force constants obtained from DFT-based phonon calculations further supports the tendency of Pd occupying the OH sublattice compared with the IT and CO sublattice. The catalytic ensembles changing from Pd monomers (Pd 1 ) to trimers (Pd 3 ) on the surface of γ phase are attributed to the increase of Pd occupancy in the OH sublattice.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling dark- and light-induced crystal structures and single-crystal optical absorption spectra of ruthenium-based complexes that undergo SO 2 -linkage photoisomerization

A family of coordination complexes of the type [Ru(SO 2 )(NH 3 ) 4 X] m+ Y n – (m, n = 1 or 2) exhibit optical switching capabilities in their single-crystal states. This striking effect is caused by the light-induced formation of SO 2 -linkage photoisomers, which are metastable if kept at suitably cool temperatures. We modeled the dark- and light-induced states of these large crystalline complexes via plane-wave (PW)- and molecular-orbital (MO)-based density functional theory (DFT) and time-dependent DFT in order to calculate their structural and optical properties; the calculated results are compared with experimental data. We show that the PW-DFT-based periodic models replicate the structural properties of these complexes more effectively than the MO-DFT-based molecular-fragment models, observing only small deviations in key bond lengths relative to the experimentally derived crystal structures. The periodic models were also found to more effectively simulate trends seen in experimental optical absorption spectra, with optical absorbance and coverage of the visible region increasing with the formation of the photoinduced geometries. The contribution of the metastable photoisomeric species more heavily focuses on the lower-energy end of the spectra. Spectra generated from the molecular-fragment models are limited by the geometry of the fragment used and the number of excited-state roots considered in those calculations. In general, periodic models outperform the molecular-fragment models owing to their ability to better appreciate the periodic phenomena that are present in these crystalline materials as opposed to MO approaches, which are finite methods. We thus demonstrate that PW-DFT-based periodic models should be considered as a more than viable method for simulating the optical and electronic properties of these single-crystal optical switches.

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

Synthesis of Co‐Doped MoS 2 Monolayers with Enhanced Valley Splitting

Abstract Internal magnetic moments induced by magnetic dopants in MoS 2 monolayers are shown to serve as a new means to engineer valley Zeeman splitting (VZS). Specifically, successful synthesis of monolayer MoS 2 doped with the magnetic element Co is reported, and the magnitude of the valley splitting is engineered by manipulating the dopant concentration. Valley splittings of 3.9, 5.2, and 6.15 meV at 7 T in Co‐doped MoS 2 with Co concentrations of 0.8%, 1.7%, and 2.5%, respectively, are achieved as revealed by polarization‐resolved photoluminescence (PL) spectroscopy. Atomic‐resolution electron microscopy studies clearly identify the magnetic sites of Co substitution in the MoS 2 lattice, forming two distinct types of configurations, namely isolated single dopants and tridopant clusters. Density functional theory (DFT) and model calculations reveal that the observed enhanced VZS arises from an internal magnetic field induced by the tridopant clusters, which couples to the spin, atomic orbital, and valley magnetic moment of carriers from the conduction and valence bands. The present study demonstrates a new method to control the valley pseudospin via magnetic dopants in layered semiconducting materials, paving the way toward magneto‐optical and spintronic devices.

Zhou, Jiadong↗
DOE OSTI2022

Stable Anode–Free All–Solid–State Lithium Battery through Tuned Metal Wetting on the Copper Current Collector

A stable anode-free all-solid-state battery (AF-ASSB) with sulfide-based solid-electrolyte (SE) (argyrodite Li 6 PS 5 Cl) is achieved by tuning wetting of lithium metal on “empty” copper current-collector. Lithiophilic 1 µm Li 2 Te is synthesized by exposing the collector to tellurium vapor, followed by in situ Li activation during the first charge. The Li 2 Te significantly reduces the electrodeposition/electrodissolution overpotentials and improves Coulombic efficiency (CE). During continuous electrodeposition experiments using half-cells (1 mA cm –2 ), the accumulated thickness of electrodeposited Li on Li 2 Te–Cu is more than 70 µm, which is the thickness of the Li foil counter-electrode. Full AF-ASSB with NMC811 cathode delivers an initial CE of 83% at 0.2C, with a cycling CE above 99%. Additionally, cryogenic focused ion beam (Cryo-FIB) sectioning demonstrates uniform electrodeposited metal microstructure, with no signs of voids or dendrites at the collector-SE interface. Electrodissolution is uniform and complete, with Li 2 Te remaining structurally stable and adherent. By contrast, an unmodified Cu current-collector promotes inhomogeneous Li electrodeposition/electrodissolution, electrochemically inactive “dead metal,” dendrites that extend into SE, and thick non-uniform solid electrolyte interphase (SEI) interspersed with pores. Density functional theory (DFT) and mesoscale calculations provide complementary insight regarding nucleation-growth behavior. Unlike conventional liquid-electrolyte metal batteries, the role of current collector/support lithiophilicity has not been explored for emerging AF-ASSBs.

36 MATERIALS SCIENCE↗