Search NASA⌕ Search

SEARCH · Search NASA

Results for “Density function 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

DMFTwDFT: An open-source code combining Dynamical Mean Field Theory with various density functional theory packages

We report Dynamical Mean Field Theory (DMFT) is a successful method to compute the electronic structure of strongly correlated materials, especially when it is combined with density functional theory (DFT). Here, we present an open-source computational package (and a library) combining DMFT with various DFT codes interfaced through the Wannier90 package. The correlated subspace is expanded as a linear combination of Wannier functions introduced in the DMFT approach as local orbitals. In particular, we provide a library mode for computing the DMFT density matrix. This library can be linked and then internally called from any DFT package, assuming that a set of localized orbitals can be generated in the correlated subspace. The existence of this library allows developers of other DFT codes to interface with our package and achieve the charge-self-consistency within DFT+DMFT loops. To test and check our implementation, we computed the density of states and the band structure of well-known solid-state correlated materials, namely LaNiO, SrVO, and NiO. The obtained results are compared to those obtained from other DFT+DMFT implementations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗
DOE OSTI2022

Electrochemical CO2 Reduction over Metal-/Nitrogen-Doped Graphene Single-Atom Catalysts Modeled Using the Grand-Canonical Density Functional Theory

Renewably driven, electrochemical conversion of carbon dioxide into value-added products is expected to be a critical tool in global decarbonization. However, theoretical studies based on the computational hydrogen electrode largely ignore the nonlinear effects of the applied potential on the calculated results, leading to inaccurate predictions of catalytic behavior or mechanistic pathways. Here, we use grand canonical density functional theory (GC-DFT) to model electrochemical CO2 reduction (CO2R) over metal- and nitrogen-doped graphene catalysts (MNCs) and explicitly include the effects of the applied potential. We used GC-DFT to compute the CO2 to CO reaction intermediate energies at -0.3, -0.7, and -1.2 VSHE catalyzed by MNCs each doped with 1 of the 10 3d block metals coordinated by four pyridinic nitrogen atoms. Our results predict that Sc-, Ti-, Co-, Cu-, and Zn-N4Cs effectively catalyze CO2R at moderate to large reducing potentials (-0.7 to -1.2 VSHE). ZnN4C is a particularly promising electrocatalyst for CO2R to CO both at low and moderate applied potentials based on our thermodynamic analysis. Our findings also explain the observed pH independence of CO production over FeN4C and predict that the rate-determining step of CO2R over FeN4C is not *CO2- formation but rather *CO desorption. Additionally, the GC-DFT-computed density of states analysis illustrates how the electronic states of MNCs and adsorbates change non-uniformly with applied potential, resulting in a significantly increased *CO2- stability relative to other intermediates and demonstrating that the formation of the adsorbed *CO2- anion is critical to CO2R activation. This work demonstrates how GC-DFT paves the way for physically realistic and accurate theoretical simulations of reacting electrochemical systems.

CO2 reduction↗
DOE OSTI2021

Accurate Thermochemistry of Complex Lignin Structures via Density Functional Theory, Group Additivity, and Machine Learning

A molecular-level understanding of lignin structures and bond dissociation energies could facilitate depolymerization technologies. Still, this information is currently limited due to the lack of databases and the simplification of surrogate models. Here, substitution effects on seven common linkages in lignin polymers are systematically investigated. An automated reaction network generator is employed to create a database of structures. A new group additivity (GA) model based on principal component analysis (PCA) descriptors is introduced and trained on gas-phase density functional theory data of 4100 species at the M06-2X/6-311++G(d,p) level. Hydrogen bonds, local steric, and nonaromatic ring contributions are also incorporated. Lastly, we improve the accuracy of the group additivity model to reach the G4 theory by computing a data set of 770 species at this level and using a data fusion approach.

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

Bulk Diffusion of Cl through O Vacancies in α-Cr 2 O 3 : A Density Functional Theory Study

Many metals form passive oxides which are broken or weakened by anionic-induced degradation leading to material failure and anionic diffusion is an important step of this degradation process. The diffusion of anionic species through oxides involves a combination of diffusion along grain boundaries, cracks, and channels, as well as bulk diffusion via point defects which is the focus of this study. Using density functional theory, we study bulk diffusion of Cl through O vacancies in α -Cr 2 O 3 as a model system for passive metal oxides. Little is known about Cl diffusion and bonding characteristics, so we benchmark our work through comparison to numerous studies on O diffusion in α -Cr 2 O 3 to analyze similarities and differences between the O and Cl diffusion in the passivation and degradation process of α -Cr 2 O 3 respectively. Unlike O diffusion, the lowest diffusion barrier for Cl is cross-plane diffusion between two (0001) planes through a vacant cation site but the much shorter in-plane diffusion path within the same coordination polyhedron has 35% higher barrier. This work provides the basis for considering the contributions of Cl bulk-diffusion in the overall diffusion kinetics of Cl through α -Cr 2 O 3 .

Oware Sarfo, Kofi↗
DOE OSTI2020

Electronic Structures and Magnetism of Zr-, Th-, and U-based Metal-Organic Frameworks (MOFs) by Density Functional Theory

Metal-organic frameworks (MOFs) have recently gained wide interest as candidate materials for nuclear waste immobilization. While the fundamental thermodynamic properties, such as the substitution energies determine the favorability of radionuclide sequestration by utilization of a MOF matrix, the studies of MOF electronic structure reveal the role of d-, and/or f-electrons on changes in physical properties of actinide-containing materials. We use density functional theory (DFT) calculations to investigate the electronic structures of Zr-, Th-, and U-MOFs, including their electronic band structures and, where appropriate, their magnetic properties. We employ various DFT methods including DFT+U, collinear spin-polarization, spin-orbit coupling, and different flavors of exchange-correlation functionals to assess the robustness to the specific exchange-correlation functional. Unlike the Zr-, and Th-MOFs, the U-MOF is found to be sensitive to electron localization and spin; hence we explore the magnetic structure of the U-MOF in further detail.

Metal-Organic Framework, MOF, Actinides, DFT↗
DOE OSTI2024

Adsorption Properties of Au−Ni Surface Alloys with a Nonstoichiometric Moiré Structure: A Density Functional Theory Study

Due to the large lattice mismatch between gold and nickel, gold–nickel surface alloys can form unique nonstoichiometric overlayer structures characterized by a moiré pattern and subsurface defects. For this work, we performed density functional theory (DFT) calculations to study the adsorption of molecular oxygen, atomic hydrogen, and atomic carbon on a gold–nickel(111) surface alloy with 0.46 monolayer gold randomly distributed in the surface layer. We observed six distinct adsorption structures for molecular oxygen characterized by intramolecular stretching frequencies of <700, 729, 795, 857, 929, and 1004 cm –1 , which describe well the experimentally observed high-resolution electron energy-loss spectra. Surface atomic hydrogen adsorption is associated with adsorbate–surface modes in the ∼1000 cm –1 range, while subsurface hydrogen can have features as low as ∼400 cm –1 . We observed a unique adsorption structure for atomic carbon inside the surface dislocation loop defect, which explains the experimentally observed low carbon-surface mode at ∼340 cm –1 . Our study sheds light on the unique adsorption properties of the gold–nickel surface alloys and helps with rationalizing vibrational frequency experimental studies for this system.

adsorption↗
DOE OSTI2022

Predicting band gaps and band-edge positions of oxide perovskites using density functional theory and machine learning

Density functional theory (DFT) within the local or semilocal density approximations, i.e., the local density approximation (LDA) or generalized gradient approximation (GGA), has become a workhorse in the electronic structure theory of solids, being extremely fast and reliable for energetics and structural properties, yet remaining highly inaccurate for predicting band gaps of semiconductors and insulators. The accurate prediction of band gaps using first-principles methods is time consuming, requiring hybrid functionals, quasiparticle GW, or quantum Monte Carlo methods. Efficiently correcting DFT-LDA/GGA band gaps and unveiling the main chemical and structural factors involved in this correction is desirable for discovering novel materials in high-throughput calculations. In this direction, we, in this study, use DFT and machine learning techniques to correct band gaps and band-edge positions of a representative subset of ABO 3 perovskite oxides. Relying on the results of HSE06 hybrid functional calculations as target values of band gaps, we find a systematic band-gap correction of ~1.5 eV for this class of materials, where ~1 eV comes from downward shifting the valence band and ~0.5 eV from uplifting the conduction band. The main chemical and structural factors determining the band-gap correction are determined through a feature selection procedure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Conditional probability density functional theory

Here we present conditional probability (CP) density functional theory (DFT) as a formally exact theory. In essence, CP-DFT determines the ground-state energy of a system by finding the CP density from a series of independent Kohn-Sham (KS) DFT calculations. By directly calculating CP densities, we bypass the need for an approximate XC energy functional. In this paper we discuss and derive several key properties of the CP density and corresponding CP-KS potential. Illustrative examples are used throughout to help guide the reader through the various concepts and theory presented. We explore a suitable CP-DFT approximation and discuss exact conditions, limitations, and results for selected examples.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗
DOE OSTI2020

Density Functional Theory Study of the Adsorption and Dissociation of Copper(I) Acetamidinates on Ni(110): The Effect of the Substrate

In order to assess the role of a metal substrates on the thermal chemistry of adsorbed acetamidinate metalorganic compounds, we have studied the surface chemistry of copper(I)–N,N'-dimethylacetamidinate on Ni(110) using density functional theory and contrasted it with similar calculations we previously carried out on Cu(110). At low coverages, it was found that, in its most stable configuration, the molecular adsorption of copper(I)–N,N'-dimethylacetamidinate dimers occurs with the Cu atoms occupying surface hollow sites. In this work, the ligands reorient away from those metal centers, and the N atoms develop new direct bonds with surface Ni atoms. However, this configuration is not stable and decomposes by losing both ligands to the surface. In the final state, the two ligands bind via their N atoms to Ni sites one lattice space away from the sites where the Cu atoms remain, with their molecular planes perpendicular to that of the surface. This is in contrast with what happens in the case of adsorption on Cu(110), where the Cu atoms from the metalorganic complex still occupy hollow sites but where only one of the ligands breaks away and binds directly to the surface; the other remains on top of the two Cu ions. In terms of the energetics of adsorption and decomposition, the reactions on Ni(110) are much more exothermic than on Cu(110). Further analysis of the distribution of charge within the adsorbates shows a minor reduction of the Cu atoms of the dimer upon interaction with the surface; full reduction to metallic copper is complete only when both ligands have fully migrated to their new Ni surface sites.

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

DMC-ICE13 : Ambient and high pressure polymorphs of ice from diffusion Monte Carlo and density functional theory

Ice is one of the most important and interesting molecular crystals, exhibiting a rich and evolving phase diagram. Recent discoveries mean that there are now 20 distinct polymorphs; a structural diversity that arises from a delicate interplay of hydrogen bonding and van der Waals dispersion forces. This wealth of structures provides a stern test of electronic structure theories, with Density Functional Theory (DFT) often not able to accurately characterize the relative energies of the various ice polymorphs. Thanks to recent advances that enable the accurate and efficient treatment of molecular crystals with Diffusion Monte Carlo (DMC), we present here the DMC-ICE13 dataset; a dataset of lattice energies of 13 ice polymorphs. This dataset encompasses the full structural complexity found in the ambient and high-pressure molecular ice polymorphs, and when experimental reference energies are available, our DMC results deliver sub-chemical accuracy. Using this dataset, we then perform an extensive benchmark of a broad range of DFT functionals. Of the functionals considered, revPBE-D3 and RSCAN reproduce reference absolute lattice energies with the smallest error, while optB86b-vdW and SCAN+rVV10 have the best performance on the relative lattice energies. Our results suggest that a single functional achieving reliable performance for all phases is still missing, and that care is needed in the selection of the most appropriate functional for the desired application. The insights obtained here may also be relevant to liquid water and other hydrogen-bonded and dispersion-bonded molecular crystals.

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

Minimum-Energy Conical Intersections by Compressed Multistate Pair-Density Functional Theory

Compressed multistate pair-density functional theory (CMS-PDFT) is a multistate version of multiconfiguration pair-density functional theory that can capture the correct topology of coupled potential energy surfaces (PESs) around conical intersections. Here, in this work, we develop interstate coupling vectors (ISCs) for CMS-PDFT in the OpenMolcas and PySCF/mrh electronic structure packages. Yet, the main focus of this work is using ISCs to calculate minimum-energy conical intersections (MECIs) by CMS-PDFT. This is performed using the projected constrained optimization method in OpenMolcas, which uses ISCs to restrain the iterations to the conical intersection seam. We optimize the S 1 /S 0 MECIs for ethylene, butadiene, and benzene and show that CMS-PDFT gives smooth PESs in the vicinities of the MECIs. Furthermore, the CMS-PDFT MECIs are in good agreement with the MECI calculated by the more expensive XMS-CASPT2 method.

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

Substituted ALPO-5 Zeolites as Promising O 2 Sorption Pump Materials: A Density Functional Theory Study

Many chemical processes depend on having an environment that is low in oxygen partial pressure (P O2 < 100 Pa); sorption pumps are a promising route to establishing that environment by either oxygen pumping or oxygen separation from an inert gas. Near ambient sorption-based processes rely on either pressure or thermal swings, requiring no moving parts, no electricity, and neither very high nor very low temperatures. In this work, we use ab initio calculations to explore zeolites as a class of materials for sorption based oxygen pumping/separation. Our calculations indicate that while the neat ALPO-5 zeolite does not adsorb O 2 , O 2 does adsorb on zeolites selectively substituted with transition metals and metalloids and hence, can enable separation and pumping. ALPO-5 substituted with Si, Ge, Sn, Pd, Pt, Ti, V, Cr, Mn, Zr, Mo, Hf, W, Ce, and Pr provides adsorption energies ranging from -0.19 to -3.92 eV, (-) indicates exothermic process. Additionally, we provide a comprehensive understanding of what controls the adsorption energy: 1) the substi-tutions must be able to adopt an oxidation state that is more positive than the cation it replaces, and 2) the size of the pore into which the O 2 adsorbs to the wall. Additionally, we conduct a ther-modynamic analysis of a thermal swing cycle to approximate the optimal O 2 binding energy for 2 low energy O 2 pumping/separation. We find that the minimum energy cost likely occurs when the adsorption energy is in the range of 0.75 – 1.00 eV (72 – 97 kJ·mol -1 ), which corresponds to Ge, V, Pt, or Ce substituted ALPO-5.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Double and Charge-Transfer Excitations in Time-Dependent Density Functional Theory

Time-dependent density functional theory has emerged as a method of choice for calculations of spectra and response properties in physics, chemistry, and biology, with its system-size scaling enabling computations on systems much larger than otherwise possible. While increasingly complex and interesting systems have been successfully tackled with relatively simple functional approximations, there has also been increasing awareness that these functionals tend to fail for certain classes of approximations. Here I review the fundamental challenges the approximate functionals have in describing double excitations and charge-transfer excitations, which are two of the most common impediments for the theory to be applied in a black-box way. At the same time, I describe the progress made in recent decades in developing functional approximations that give useful predictions for these excitations.

Chemistry↗
DOE OSTI2023

The Predictive Power of Exact Constraints and Appropriate Norms in Density Functional Theory

Ground-state Kohn-Sham density functional theory provides, in principle, the exact ground-state energy and electronic spin densities of real interacting electrons in a static external potential. In practice, the exact density functional for the exchange-correlation (xc) energy must be approximated in a computationally efficient way. About 20 mathematical properties of the exact xc functional are known. In this work, we review and discuss these known constraints on the xc energy and hole. By analyzing a sequence of increasingly sophisticated density functional approximations (DFAs), we argue that (a) the satisfaction of more exact constraints and appropriate norms makes a functional more predictive over the immense space of many-electron systems and (b) fitting to bonded systems yields an interpolative DFA that may not extrapolate well to systems unlike those in the fitting set. We discuss both how the class of well-described systems has grown along with constraint satisfaction and the possibilities for future functional development.

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

Thermodynamic properties of metaschoepite predicted from density functional perturbation theory

Density functional perturbation theory (DFPT) calculations of the thermodynamic properties of metaschoepite, (UO 2 ) 8 0 2 (OH) 12 •10H 2 0, are reported. Using a recently revised crystal structure of metaschoepite, the predicted molar entropy and isobaric heat capacity are overall significantly smaller than previous calculations using an earlier orthorhombic crystal structure model. The present DFPT calculations also show large differences between the thermo-dynamic functions of metaschoepite and schoepite, which might reflect the change in phonon properties upon removal of two H 2 0 molecules per formula unit and alteration of the H-bonded interlayer water network from schoepite to metaschoepite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗