Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density functional theory”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Computation of Auger Electron Spectra in Organic Molecules with Multiconfiguration Pair-Density Functional Theory

Efficient and accurate computation of molecular Auger electron spectra for larger systems is limited by the rapid increase in the number of doubly ionized final states as the system size grows. Here, in this work, we benchmark the application of multiconfiguration pair-density functional theory with a restricted active space (RAS) reference wave function for computing the carbon K-edge decay spectra of 20 organic molecules. Decay rates are computed within the one-center approximation. We evaluate the performance of different basis sets and on-top functionals and find that multiconfiguration pair-density functional theory achieves accuracy comparable to RAS followed by second-order perturbation theory, but at significantly lower computational cost.

Fouda, Adam E. A. [Argonne National Laboratory (AN↗
DOE OSTI2020

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Relativistic density functional theory in the full potential linear muffin tin orbital method

In this write-up, we have outlined the general methodology for treating relativity in the FP-LMTO method. The Schrodinger equation is insufficient for describing the physics of heavy elements because it leaves out any relativistic effects. Therefore DFT based on the Schrodinger equation is not sufficient for studying heavy materials. It is, however, possible to incorporate all of the effects of relativity in DFT if starting with a reformulation of the Hohenberg-Kohn theorem that proceeds from the Dirac equation, rather than the Schrodinger equation. The corresponding theorem in relativistic density functional theory (RDFT) shows that the ground-state density of the system is a functional of not only the charge density as in Schrodinger-based DFT, but also of the current density. This leads to a more elaborate procedure in practice when solving the Dirac-Kohn-Sham equations, the relativistic analogs of the Kohn-Sham equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Excitation energy of fission fragments within nuclear time-dependent density functional theory

The number and properties of the neutrons and photons emitted in nuclear fission are directly related to the excitation energy of the fission fragments when they are formed at scission. Though not observable experimentally because of the extremely short timescales, the excitation energy of fission fragments can be predicted by microscopic theory based on time-dependent density functional theory (TDDFT). Initial results on the value of the total kinetic energy of fission reactions were very promising, but could not probe all possible fragmentations. Here, in this work, we perform large-scale TDDFT calculations in 240 Pu enabled by the development of a new TDDFT solver. We obtain TDDFT trajectories covering nearly all possible fragmentations. We find that the total kinetic energy is close to experimental values only for the most likely fission while it is severely underestimated at both small and large asymmetries. This conclusion seems rather independent of the parametrization of the energy functional, both in its particle-hole and particle-particle channels.

Physics - Nuclear physics and radiation physics↗
DOE OSTI2022

Zero-Field Splitting Calculations by Multiconfiguration Pair-Density Functional Theory

Zero-field splitting (ZFS) is a fundamental molecular property that is especially relevant for single-molecule magnets (SMMs), electron paramagnetic resonance spectra, and quantum computing. Developing a method that can accurately predict ZFS parameters can be very powerful for designing new SMMs. One of the challenges is to include external correlation in an inherently multiconfigurational open-shell species for the accurate prediction of magnetic properties. Previously available methods depend on expensive multireference perturbation theory calculations to include external correlation. In this paper, we present spin-orbit-inclusive multiconfiguration and multistate pair-density functional theory (MC-PDFT) calculations of ZFSs; these calculations have a cost comparable to complete-active-space self-consistent field (CASSCF) theory, but they include correlation external to the active space. We found that combining a multistate formulation of MC-PDFT, namely, compressed-state multistate pair-density functional theory, with orbitals optimized by weighted-state-averaged CASSCF, yields reasonably accurate ZFS results

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

GDB-9-Ex_TD-DFT-PBE0: Dataset containing Time Dependent Density Functional Theory (TDDFT) calculations for organic molecules of the GDB-9-Ex dataset.

This dataset contains data-intensive quantum chemical electronic structure calculations for 96,766 organic molecules of the GDB-9-Ex dataset. Calculations were performed using the Time Dependent Density Functional Theory (TDDFT) first principles method using the ORCA software. It provides UV-vis spectra calculations of molecules with a high level of accuracy. The optical spectra behavior was collected based on the optimized molecular geometries in the DFTB method with 3ob parameters. All calculations utilized the def2-TZVP basis sets with the auxiliary def2/J and def2-TZVP/C basis sets. The time-dependent density-functional theory (TDDFT) approach with the PBE0 exchange-correlation functional and ORCAs default integration grid was employed. For the excitation energy calculations, the lowest 50 excitation states were calculated.

AI dataset↗
DOE OSTI2023

Realizing the potentials of density functional theory (DFT) and of the materials genome initiative (MGI)

Abstract From 1964 and 1965 to present, the wide spread utilization of an incomplete density functional theory (DFT) has led to mixed results: The second theorem of the theory asserts that the energy functional reaches its minimum if the calculation employs the ground state charge density—without providing a mechanism for finding this density. Calculations purporting to employ DFT have mostly assumed that results obtained with a judiciously selected basis set, following self-consistent iterations, are those of the ground state. The state obtained with a single basis set is a stationary one, among an infinite number of such states, with no proven relation to the actual ground state of the material. Most failures or limitations of the incomplete DFT can be traced to this error. We present results from calculations using the completed DFT. They are in excellent agreement with experiment and portend the realization of the Materials Genome Initiative.

Bagayoko, Diola↗
DOE OSTI2022

Importance of imposing gauge invariance in time-dependent density functional theory calculations with meta-generalized gradient approximations

It has been known for more than a decade that the gauge variance of the kinetic energy density τ leads to additional terms in the magnetic orbital rotation Hessian used in linear-response time-dependent density functional theory (TDDFT), affecting excitation energies obtained with τ-dependent exchange–correlation functionals. While previous investigations found that a correction scheme based on the paramagnetic current density has a small effect on benchmark results, we report more pronounced effects here, in particular, for the popular M06-2X functional and for some other meta-generalized gradient approximations (mGGAs). In the first part of this communication, this is shown by a reassessment of a set of five Ni(II) complexes for which a previous benchmark study that did not impose gauge invariance has found surprisingly large errors for excitation energies obtained with M06-2X. These errors are more than halved by restoring gauge invariance. The variable importance of imposing gauge invariance for different mGGA-based functionals can be rationalized by the derivative of the mGGA exchange energy integrand with respect to τ. In the second part, a large set of valence excitations in small main-group molecules is analyzed. For M06-2X, several selected n → π* and π→π$^{*}_{⊥}$ excitations are heavily gauge-dependent with average changes of –0.17 and –0.28 eV, respectively, while π→π$^{*}_{∥}$ excitations are marginally affected (–0.04 eV). Similar patterns, but of the opposite signs, are found for SCAN0. Here, the results suggest that reevaluation of previous gauge variant TDDFT results based on M06-2X and other mGGA functionals is warranted.

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

Machine learning-enhanced design of lead-free halide perovskite materials using density functional theory

The investigation of emerging non-toxic perovskite materials has been undertaken to advance the fabrication of environmentally sustainable lead-free perovskite solar cells. This study introduces a machine learning methodology aimed at predicting innovative halide perovskite materials that hold promise for use in photovoltaic applications. The seven newly predicted materials are as follows: CsMnCl 4 , Rb 3 Mn 2 Cl 9 , Rb 4 MnCl 6 , Rb 3 MnCl 5 , RbMn 2 Cl 7 , RbMn 4 Cl 9 , and CsIn 2 Cl 7 . The predicted compounds are first screened using a machine learning approach, and their validity is subsequently verified through density functional theory calculations. CsMnCl 4 is notable among them, displaying a bandgap of 1.37 eV, falling within the Shockley-Queisser limit, making it suitable for photovoltaic applications. Here, through the integration of machine learning and density functional theory, this study presents a methodology that is more effective and thorough for the discovery and design of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unconventional Error Cancellation Explains the Success of Hartree–Fock Density Functional Theory for Barrier Heights

Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally efficient semilocal approximations for the exchange-correlation energy. The accuracy of a semilocal density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree–Fock electron densities, which has been known for ~30 years. Here, the conventional explanation is that the Hartree–Fock theory yields the more accurate density. This work presents a benchmark Kohn–Sham inversion of accurate coupled-cluster densities for the reaction H 2 + F → HHF → H + HF and finds a strong, understandable cancellation between positive (excessively overcorrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree–Fock density functional theory. This confirms earlier conclusions (Kaplan, A. D., et al. J. Chem. Theory Comput. 2023, 19, 532–543) based on 76 barrier heights and three less reliable, but less expensive, fully nonlocal density functional proxies for the exact density.

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

A parallel strategy for density functional theory computations on accelerated nodes

Using the Löwdin orthonormalization of tall-skinny matrices as a proxy-app for wavefunction-based Density Functional Theory solvers, we investigate a distributed memory parallel strategy focusing on Graphics Processing Unit (GPU)-accelerated nodes as available on some of the top ranked supercomputers at the present time. Here we present numerical results in the strong limit regime, as it is particularly relevant for First-Principles Molecular Dynamics. We also examine how matrix product-based iterative solvers provide a competitive alternative to dense eigensolvers on GPUs, allowing to push the strong scaling limit of these computations to a larger number of distributed tasks. Our strategy, which relies on replicated Gram matrices and efficient collective communications using the NCCL library, leads to a time-to-solution under 0.5 s for the Löwdin orthonormalization of a tall-skinny matrix of 3000 columns on Summit at Oak Ridge Leadership Facility (OLCF). Given the similarity in computational operations between one iteration of a DFT solver and this proxy-app, this shows the possibility of solving accurately the DFT equations well under a minute for 3000 electronic wave functions, and thus perform First-Principles molecular dynamics of physical systems much larger than traditionally solved on CPU systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Machine learned synthesizability predictions aided by density functional theory

Abstract A grand challenge of materials science is predicting synthesis pathways for novel compounds. Data-driven approaches have made significant progress in predicting a compound’s synthesizability; however, some recent attempts ignore phase stability information. Here, we combine thermodynamic stability calculated using density functional theory with composition-based features to train a machine learning model that predicts a material’s synthesizability. Our model predicts the synthesizability of ternary 1:1:1 compositions in the half-Heusler structure, achieving a cross-validated precision of 0.82 and recall of 0.82. Our model shows improvement in predicting non-half-Heuslers compared to a previous study’s model, and identifies 121 synthesizable candidates out of 4141 unreported ternary compositions. More notably, 39 stable compositions are predicted unsynthesizable while 62 unstable compositions are predicted synthesizable; these findings otherwise cannot be made using density functional theory stability alone. This study presents a new approach for accurately predicting synthesizability, and identifies new half-Heuslers for experimental synthesis.

Lee, Andrew (ORCID:0000000153014295)↗
DOE OSTI2025

Accuracy of Kohn–Sham density functional theory for warm- and hot-dense matter equation of state

We study the accuracy of Kohn–Sham density functional theory (DFT) for warm- and hot-dense matter (WDM and HDM). Specifically, considering a wide range of systems, we perform accurate ab initio molecular dynamics simulations with temperature-independent local/semilocal density functionals to determine the equations of state at compression ratios of 3x–7x and temperatures near 1 MK. We find very good agreement with path integral Monte Carlo benchmarks, while having significantly smaller error bars and smoother data, demonstrating the accuracy of DFT for the study of WDM and HDM at such conditions. In addition, using a Δ-machine learned force field scheme, we confirm that the DFT results are insensitive to the choice of exchange-correlation functional, whether local, semilocal, or nonlocal.

Suryanarayana, Phanish (ORCID:0000000151720049)↗
DOE OSTI2021

Localized Active Space Pair-Density Functional Theory

Accurate quantum chemical methods for the prediction of spin-state energy gaps for strongly correlated systems are computationally expensive and scale poorly with the size of the system. This makes calculations for many experimentally interesting molecules impractical even with abundant computational resources. Previous work has shown that the localized active space (LAS) self-consistent field (SCF) method can be an efficient way to obtain multiconfiguration SCF wave functions of comparable quality to the corresponding complete active space (CAS) ones. To obtain quantitative results, a post-SCF method is needed to estimate the complete correlation energy. One such method is multiconfiguration pair-density functional theory (PDFT), which calculates the energy based on the density and on-top pair density obtained from a multiconfiguration wave function. In this work, we introduce localized-active-space PDFT, which uses a LAS wave function for subsequent PDFT calculations. The method is tested by computing spin-state energies and gaps in conjugated organic molecules and a bimetallic compound and comparing to the corresponding CAS-PDFT values.

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

Mechanistic Insights into the Demethylation of Lignin‐Derived Structures Using Protic Ionic Liquids: A Density Functional Theory Study

Lignin valorization is restricted by the stability of its methoxy groups, creating a critical need for efficient demethylation strategies. Here, in this study, density functional theory (DFT) is employed to dissect the mechanistic pathways of demethylation in lignin model compounds, guaiacol and syringol, using protic ionic liquids (PILs) that act as both solvent and catalyst. Conductor‐like Screening Model for Real Solvents (COSMO‐RS) analysis identifies monoethanolammonium acetate ([MEOA][Ace]) as the most promising medium, attributed to its strong hydrogen bonding network and solvation ability. By integrating implicit and explicit solvation models, it is revealed that an acid‐catalyzed hydrolytic mechanism governs demethylation, with PILs stabilizing crucial transition states and intermediates. Complementary electronic structure evaluations, including highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO‐LUMO) gap analysis, charge distribution, and electrostatic potential mapping, demonstrate how PILs lower energetic barriers and enhance reactivity. To simulate realistic environments, this study is extended to lignin dimer complexes with varying water content, uncovering how water‐bridged solvation reverses demethylation preference from guaiacyl to syringyl units. This mechanistic shift aligns with experimental observations showing faster S‐unit reactivity in hydrated systems. Together, these findings provide atomic‐level insight into lignin demethylation dynamics and highlight how tuning acid concentration in PILs can accelerate kinetics, enabling a rational pathway toward next‐generation biomass conversion technologies.

COSMO-RS↗
DOE OSTI2021

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

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

Analysis and visualization of energy densities. I. Insights from real-time time-dependent density functional theory simulations

In this article, we report a scheme to analyze and visualize the energy density fluctuations during the real-time time-dependent density functional theory (RT-TDDFT) simulations. Using Ag 4 –N 2 complexes as examples, here we show that the grid-based Kohn–Sham energy density can be computed at each time step using a procedure from Nakai and coworkers. Then the instantaneous energy of each molecular fragment (such as Ag 4 and N 2 ) can be obtained by partitioning the Kohn–Sham energy densities using Becke or fragment-based Hirshfeld (FBH) scheme. A strong orientation-dependence is observed for the energy flow between the Ag 4 cluster and a nearby N 2 molecule in the RT-TDDFT simulations. Future applications of such an energy density analysis in electron dynamics simulations are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗