Search NASA⌕ Search

SEARCH · Search NASA

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗
DOE OSTI2021

A deep learning interatomic potential developed for atomistic simulation of carbon materials

Interatomic potentials based on neural-network machine learning method have attracted considerable attention in recent years owing to their outstanding ability to balance the accuracy and efficiency in atomistic simulations. In this work, a neural-network potential (NNP) for carbon is generated to simulate the structural properties of various carbon structures. The potential is trained using a database consisting of crystalline and liquid structures obtained by the first-principles density functional theory (DFT) calculations. The developed potential accurately predicts the energies and forces in crystalline and liquid carbon structures, the energetic stability of defected graphene, and the structures of amorphous carbon as the function of density. As a result, the excellent accuracy and transferability of the NNP provide a promising tool for accurate atomistic simulations of various carbon materials with faster speed and much lower cost.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Learning local and semi-local density functionals from exact exchange-correlation potentials and energies

Finding accurate exchange-correlation (XC) functionals remains the defining challenge in density functional theory (DFT). Despite 40 years of active development, attaining general purpose chemical accuracy is still elusive with existing functionals. We present a data-driven pathway to learn the XC functional by using the exact density, XC energy, and XC potential. While the exact densities are obtained from accurate configuration interaction (CI), the exact XC energies and XC potentials are obtained via inverse DFT calculations on the CI densities. We demonstrate how simple neural network (NN)–based local density approximation (LDA) and generalized gradient approximation (GGA), trained on just five atoms and two molecules, provide remarkable improvement in total energies and densities. Particularly, the NN-based GGA functional attains similar accuracy as the higher rung SCAN meta-GGA on various thermochemistry datasets. These results underscore the promise of using the XC potential in modeling XC functionals and can pave the way for systematic learning of increasingly accurate XC functionals.

Science & Technology - Other Topics↗
DOE OSTI2022

First-Principles Calculations of the Structural, Electronic, Optical, and Mechanical Properties of 21 Pyrophosphate Crystals

Pyrophosphate crystals have a wide array of applications in industrial and biomedical fields. However, fundamental understanding of their electronic structure, optical, and mechanical properties is still scattered and incomplete. In the present research, we report a comprehensive theoretical investigation of 21 pyrophosphates A 2 M (H 2 P 2 O 7 ) 2 •2H 2 O with either triclinic or orthorhombic crystal structure. The molecule H 2 P 2 O 7 is the dominant molecular unit, whereas A = (K, Rb, NH 4 , Tl), M = (Zn, Cu, Mg, Ni, Co, Mn), and H 2 O stand for the cation elements, transition metals, and the water molecules, respectively. The electronic structure, interatomic bonding, partial charge distribution, optical properties, and mechanical properties are investigated by first-principles calculations based on density functional theory (DFT). Most of these 21 crystals are theoretically investigated for the first time. The calculated results show a complex interplay between A, M, H 2 P 2 O 7 , and H 2 O, resulting in either metallic, half-metallic, or semi-conducting characteristics. The novel concept of total bond order density (TBOD) is used as a single quantum mechanical metric to characterize the internal cohesion of these crystals to correlate with the calculated properties, especially the mechanical properties. This work provides a large database for pyrophosphate crystals and a road map for potential applications of a wider variety of phosphates.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Phase Stability in U-6Nb Alloy Doped with Ti from the First Principles Theory

First-principles calculations within the density-functional-theory (DFT) approach are conducted in order to explore and explain the effect of small amounts of titanium on phase stability in the U-6Nb alloy. During rapid quenching from high to room temperature, metastable phases α' (orthorhombic), α" (monoclinic), and γ 0 (tetragonal) can form, depending on Nb concentration. Important mechanical properties depend on the crystal structure and, therefore, an understanding of the effect of impurities on phase stability is essential. Insights on this issue are obtained from quantum-mechanical DFT calculations. The DFT framework does not rely on any material-specific assumptions and is therefore ideal for an unbiased investigation of the U-Nb system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Stochastic Vector Techniques in Ground-State Electronic Structure

Herein we review a suite of stochastic vector computational approaches for studying the electronic structure of extended condensed matter systems. These techniques help reduce algorithmic complexity, facilitate efficient parallelization, simplify computational tasks, accelerate calculations, and diminish memory requirements. While their scope is vast, we limit our study to ground-state and finite temperature density functional theory (DFT) and second-order many-body perturbation theory. More advanced topics, such as quasiparticle (charge) and optical (neutral) excitations and higher-order processes, are covered elsewhere. We start by explaining how to use stochastic vectors in computations, characterizing the associated statistical errors. Next, we show how to estimate the electron density in DFT and discuss effective techniques to reduce statistical errors. Finally, we review the use of stochastic vectors for calculating correlation energies within the second-order Møller-Plesset perturbation theory and its finite temperature variational form. Example calculation results are presented and used to demonstrate the efficacy of the methods.

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

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗
DOE OSTI2024

Theoretical Prediction of Thermal Expansion Anisotropy for Y 2 Si 2 O 7 Environmental Barrier Coatings Using a Deep Neural Network Potential and Comparison to Experiment

Environmental barrier coatings (EBCs) are an enabling technology for silicon carbide (SiC)-based ceramic matrix composites (CMCs) in extreme environments such as gas turbine engines. However, the development of new coating systems is hindered by the large design space and difficulty in predicting the properties for these materials. Density Functional Theory (DFT) has successfully been used to model and predict some thermodynamic and thermo-mechanical properties of high-temperature ceramics for EBCs, although these calculations are challenging due to their high computational costs. In this work, we use machine learning to train a deep neural network potential (DNP) for Y 2 Si 2 O 7 , which is then applied to calculate the thermodynamic and thermo-mechanical properties at near-DFT accuracy much faster and using less computational resources than DFT. We use this DNP to predict the phonon-based thermodynamic properties of Y 2 Si 2 O 7 with good agreement to DFT and experiments. We also utilize the DNP to calculate the anisotropic, lattice direction-dependent coefficients of thermal expansion (CTEs) for Y 2 Si 2 O 7 . Molecular dynamics trajectories using the DNP correctly demonstrate the accurate prediction of the anisotropy of the CTE in good agreement with the diffraction experiments. In the future, this DNP could be applied to accelerate additional property calculations for Y 2 Si 2 O 7 compared to DFT or experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of Sn addition on precipitation of Al–Zn–Mg alloy at early stage of natural aging

The effect of Sn addition on the natural aging behavior of Al–Zn–Mg alloy was investigated using atom probe tomography (APT) experiments and density functional theory (DFT) calculations. The results demonstrate that Sn retards the natural aging behavior by facilitating the preferential formation of Mg–Sn clusters, which suppress the formation of Guinier-Preston zones (I) during the earliest aging stage. Quantitative solute analysis conducted by APT has suggested increasing amount of Mg–Sn clusters in Sn-contained alloy due to their high binding energy, supported by DFT calculations on the Sn.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2025

Benchmarking DFT Accuracy in Predicting O 1s Binding Energies on Metals

X-ray photoelectron spectroscopy (XPS) is a powerful tool for probing the electronic structure and composition of materials, particularly metals and metal oxides of relevance to solar cells and catalysis. Density functional theory (DFT) is often used to support XPS peak assignments, but its reliability for predicting oxygen species is not well established. Here, we compile a large data set of experimental oxygen binding energies and evaluate corresponding DFT predictions. We find that as the binding energies of metal-bound atomic oxygen species increase, especially above ≈530 eV, there is a general decrease in the accuracy of DFTpredicted values. Thus, high-binding-energy atomic oxygen species, such as those proposed as active for selective Ag-catalyzed epoxidation, are less well represented. The chemical nature of the oxygen species also influences accuracy, with molecularly bound species more reliably captured across the entire range of energies. These findings illustrate the limitations of DFT for interpreting XPS spectra and provide a benchmark for improving computational methods.

Adsorption↗
DOE OSTI2021

Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies

Abstract A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Quasicubic model for metal halide perovskite nanocrystals

In this work, we present an analysis of quantum confinement of carriers and excitons, and exciton fine structure, in metal halide perovskite (MHP) nanocrystals (NCs). Starting with coupled-band k · P theory, we derive a nonparabolic effective mass model for the exciton energies in MHP NCs valid for the full size range from the strong to the weak confinement limits. We illustrate the application of the model to CsPbBr 3 NCs and compare the theory against published absorption data, finding excellent agreement. We then apply the theory of electron-hole exchange, including both short- and long-range exchange interactions, to develop a model for the exciton fine structure. We develop an analytical quasicubic model for the effect of tetragonal and orthorhombic lattice distortions on the exchange-related exciton fine structure in CsPbBr 3 , as well as some hybrid organic MHPs of recent interest, including formamidinium lead bromide (FAPbBr 3 ) and methylammonium lead iodide (MAPbI3). Testing the predictions of the quasicubic model using hybrid density functional theory (DFT) calculations, we find qualitative agreement in tetragonal MHPs but significant disagreement in the orthorhombic modifications. Moreover, the quasicubic model fails to correctly describe the exciton oscillator strength and with it the long-range exchange corrections in these systems. Introducing the effect of NC shape anisotropy and possible Rashba terms into the model, we illustrate the calculation of the exciton fine structure in CsPbBr 3 NCs based on the results of the DFT calculations and examine the effect of Rashba terms and shape anisotropy on the calculated fine structure.

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

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2026

Synergistic Doping and Stabilization of Magnetically Tunable LnTi3(Sb,Sn)4 (Ln: Ce–Gd) Kagome Metals

Here, we present our synthesis and characterization of the LnTi3(Sb,Sn)4 (Ln: Ce, Pr, Nd, Sm, Gd) family of cleavable kagome metals. While these materials are isostructural to the LnTi3Bi4 family, they only form as (Sb,Sn) solid-solutions with no corresponding LnTi3Sb4 or LnTi3Sn4 phases. We use a combination of first-principles density functional theory (DFT) and crystal orbital Hamilton population (COHP) calculations to show that (Sb,Sn) alloying has a stabilizing effect on the structure by adjusting the Fermi level, filling bonding states, depopulating antibonding states, and adjusting the density of states (DOS) toward local minima, an effect we call “synergistic doping.” Through a detailed characterization of the SmTi3(Sb,Sn)4 series, we further demonstrate that the tunable Fermi level has a profound effect on the physical properties. We observe multiple magnetic phases that stem from a competition between antiferromagnetic (AFM) and ferromagnetic-like (FM) ground states. Furthermore, the (Sb,Sn) ratio allows us to tune from direct AFM-FM competition toward a complex admixture with properties reminiscent of the recently discovered TbTi3Bi4. Ultimately, our work demonstrates how the concept of synergistic doping provides a means to stabilize new structures while developing systems with intrinsically tunable chemical, magnetic, and electronic properties.

Ortiz, Brenden [ORNL] (ORCID:0000000213337003)↗
DOE OSTI2024

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗
DOE OSTI2020

First Principle Investigation of Electronic, Transport, and Bulk Properties of Zinc-Blende Magnesium Sulfide

We have studied electronic, structural, and transport properties of zinc-blende magnesium sulfide (zb-MgS). We employed a local density approximation (LDA) potential and the linear combination of atomic orbitals (LCAO) method. Our computational method is able to reach the ground state of a material, as dictated by the second theorem of density functional theory (DFT). Consequently, our findings have the physical content of DFT and agree with available, corresponding experimental ones. The calculated band gap of zb-MgS, a direct gap equal to 4.43 eV, obtained at the experimental lattice constant of 5.620 Å, completely agrees with the experimental band gap of 4.45 ± 0.2 eV. We also report total (DOS) and partial (pDOS) densities of states, electron and hole effective masses, the equilibrium lattice constant, and the bulk modulus. The calculated pDOS also agree with the experiment for the description of the states at the top and the bottom of the valence and conduction bands, respectively.

42 ENGINEERING↗
DOE OSTI2020

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car–Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson’s equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. Furthermore, this is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500–1000 atoms (e.g., (H₂O)₂₅₆) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗