Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional approximations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

How Accurate Are Approximate Density Functionals for Noncovalent Interaction of Very Large Molecular Systems?

Noncovalent intermolecular interactions are very important in many research areas. Therefore, it is vital to understand the extent to which approximate density functionals give a proper description of noncovalent interactions. Previous research has demonstrated that some approximate density functionals can predict usefully accurate interaction energies for many noncovalent systems; however, most of that work is limited to small and moderate-sized molecules. Very recently though, accurate benchmarks have become available for some very large molecules. Here, the present work applies 21 approximate density functionals to compute the binding energies of seven large molecular systems that have a number of atoms ranging from 200 to 910. The results are judged by comparison to the recently published CIM-DLPNO-CCSD(T) results, which are assumed to provide a reliable benchmark. The five most accurate methods among those tested are found to be PW6B95-D4, PW6B95-D3(BJ), revM11, M06-L, and MN15.

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

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

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

Sensitivity of the electronic and magnetic structures of cuprate superconductors to density functional approximations

Abstract We discuss the crystal, electronic, and magnetic structures of La 2− x Sr x CuO 4 (LSCO) for x = 0.0 and x = 0.25 employing 13 density functional approximations, representing the local, semi-local, and hybrid exchange-correlation approximations within the Perdew–Schmidt hierarchy. The meta-generalized gradient approximation (meta-GGA) class of functionals is found to perform well in capturing the key properties of LSCO, a prototypical high-temperature cuprate superconductor. In contrast, the localspin-density approximation, GGA, and the hybrid density functional fail to capture the metal-insulator transition under doping.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mechanical properties of zeolite-templated carbons from approximate density functional theory calculations

Zeolite-templated carbon (ZTC) is a unique porous carbonaceous material whose structure is ordered at the nanometre scale, enabling a representative periodic description at the atomistic level. Utilizing an existing, well-defined reference model for ZTCs, a structural library of varying compositions was developed by refinement using density-functional tight-binding (DFTB) potentials parameterized for materials science applications. We first determined the quantum chemical-refined structures of models with CH, CHO, CHON, CHOB, and CHOBN compositions with various degrees of heteroatom substitution. These structural models comprise the characteristic morphological features of highly porous carbon materials, such as open-blade surfaces, edges, saddles, and closed-strut formations, spanning a range of curvatures and characteristic sizes. Second, we carried out alternating compression and expansion of the CHO model unit cell to determine the lowest energy structure as well as to obtain its bulk modulus in order to demonstrate a close connection between macroscopic observations and atomic-scale structures. Further, the agreement between experimental measurements and the computational model is remarkable and demonstrates the power of approximate density functional theory as a cost-effective computational tool with chemical accuracy for the investigation of structure/property relationships in real-world carbon-based solids.

03 NATURAL GAS↗
DOE OSTI2021

Opening band gaps of low-dimensional materials at the meta-GGA level of density functional approximations

The quasiparticle band structure can be properly described by Hedin's GW approximation (GW), at a high computational cost. For band gaps, semilocal density functionals up to the generalized gradient approximation (GGA) level cannot compete with the accuracy of hybrid-based approximations or GW. Meta-GGA density functionals with a strong dependence on the kinetic energy density ingredient can potentially give wider band gaps compared with GGAs. The recent TASK meta-GGA density functional from Aschebrock and Kümmel [ Phys. Rev. Research 1 , 033082 (2019) ], is constructed with an enhanced nonlocality in the generalized Kohn-Sham scheme and therefore harbors great opportunities for band gap prediction. Although this approximation was found to yield excellent band gaps of bulk solids, this accuracy cannot be straightforwardly transferred to low-dimensional materials. Additionally, the reduced screening of these materials results in larger band gaps compared with their bulk counterparts, as an additional barrier to overcome. In this paper we demonstrate how the alteration of this functional affects the band gaps of monolayers and nanoribbons and present accurate band gaps competing with the revised Heyd-Scuseria-Ernzerhof (HSE06) approximation. In order to achieve this goal, we have modified the TASK functional (a) by changing the tight upper bound for one- or two-electron systems ( h X 0 ) from 1.174 to 1.29 and (b) by changing the limit of the interpolation function f X ( α → ∞ ) of the TASK functional that interpolates the exchange enhancement factor F X ( s , α ) from α = 0 to 1. The resulting modified TASK (mTASK) was tested for various materials from three dimensions to two dimensions to one dimension (nanoribbons) and was compared with the results of the higher-level hybrid functional HSE06 or with the G 0 W 0 approximation within many-body perturbation theory. We find that mTASK systematically improves the band gaps and band structures of two-dimensional (2D) and 1D systems, without significantly affecting the accuracy of the original TASK for the bulk 3D materials, when compared with the Perdew-Burke-Ernzerhof (PBE) GGA and the strongly constrained and appropriately normed (SCAN) meta-GGA. We further demonstrate the applicability of mTASK by assessing the band structures of transition metal dichalcogenide nanoribbons with respect to various bending curvatures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

On the Sensitivity to Density-Functional Approximations for CO Binding Energies of Single-Atom Catalysts in Nitrogen-Doped Graphene

Density functional theory (DFT) methods are the working horse in screening new catalytic materials. They are widely used to predict trends in binding energies, which are then used to compare the activity of different materials. Here, the binding strength of CO is an important descriptor to the CO 2 reduction catalytic activity of the single transition metal atoms embedded on nitrogen-doped graphene (TM/NG). In this work, however, we show that CO binding strengths in different TM/NG has very different sensitivity to DFT methods. Specifically, Fe/NG CO binding energy changes dramatically with the percentage of exact exchange in the functional; Co/NG does less so, while Ni/NG nearly has no change. Such varying behaviors is a direct result of different local spin configurations, similar to the performance of DFT methods for metal porphyrin complexes. Therefore, caution should be exercised when using DFT binding energies for quantitative predictions in TM/NG single atom catalysis.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Self‐Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces

The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as “Jacob's ladder.” This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBEα framework with mathematical equation being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.

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

Machine learning for accuracy in density functional approximations

Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. In this paper, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.

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

Understanding the chemical bonding of ground and excited states of HfO and HfB with correlated wavefunction theory and density functional approximations

Knowledge of the chemical bonding of HfO and HfB ground and low-lying electronic states provides essential insights into a range of catalysts and materials that contain Hf–O or Hf–B moieties. Here, we carry out high-level multi-reference configuration interaction theory and coupled cluster quantum chemical calculations on these systems. We compute full potential energy curves, excitation energies, ionization energies, electronic configurations, and spectroscopic parameters with large quadruple-ζ and quintuple-ζ quality correlation consistent basis sets. We also investigate equilibrium chemical bonding patterns and effects of correlating core electrons on property predictions. Differences in the ground state electron configuration of HfB(X 4 Σ - ) and HfO(X 1 Σ + ) lead to a significantly stronger bond in HfO than HfB, as judged by both dissociation energies and equilibrium bond distances. We extend our analysis to the chemical bonding patterns of the isovalent HfX (X = O, S, Se, Te, and Po) series and observe similar trends. We also note a linear trend between the decreasing value of the dissociation energy (D e ) from HfO to HfPo and the singlet–triplet energy gap (ΔE S–T ) of the molecule. Finally, we compare these benchmark results to those obtained using density functional theory (DFT) with 23 exchange–correlation functionals spanning multiple rungs of “Jacob’s ladder.” When comparing DFT errors to coupled cluster reference values on dissociation energies, excitation energies, and ionization energies of HfB and HfO, we observe semi-local generalized gradient approximations to significantly outperform more complex and high-cost functionals.

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

Dynamic aspects of graphene deformation and fracture from approximate density functional theory

Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.

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

Erratum: “DFTB+, a software package for efficient approximate density functional theory based atomistic simulations” [J. Chem. Phys. 152, 124101 (2020)]

The implementation of the GPU support in DFTB+, as described in Sec. III C of the original publication,1 was developed based on a previous unpublished implementation by Jacek Jakowski. In order to acknowledge his work on this first implementation, the authors of the original publication wish to include J. Jakowski as co-author. The scientific content of the original publication is not affected.

42 ENGINEERING↗
DOE OSTI2022

First-principles wave-vector- and frequency-dependent exchange-correlation kernel for jellium at all densities

Here we propose a spatially and temporally nonlocal exchange correlation (XC) kernel for the spin-unpolarized fluid phase of ground-state jellium for use in time-dependent density functional and linear response calculations. The kernel is constructed to satisfy known properties of the exact XC kernel to accurately describe the correlation energies of bulk jellium and to satisfy frequency-moment sum rules at a wide range of bulk jellium densities, including those low densities that display strong correlation and symmetry breaking. These effects are easier to understand in the simple jellium model than in real systems. All exact constraints satisfied by the recent MCP07 kernel are maintained in the revised MCP07 (rMCP07) kernel, while others are added. The revision $f^{rMCP07}_{XC}$ (q, ω) differs from MCP07 only for nonzero frequencies ω. Only at densities much lower than those of real bulk metals is the frequency dependence of the kernel important for the correlation energy of jellium. As the wave vector q tends to zero, the kernel has a -4$πα(ω)/q^2$ divergence whose frequency-dependent ultranonlocality coefficient $α(ω)$ vanishes in jellium, and is predicted by rMCP07 to be extremely small for the real metals Al and Na.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-throughput predictions of metal–organic framework electronic properties: theoretical challenges, graph neural networks, and data exploration

Abstract With the goal of accelerating the design and discovery of metal–organic frameworks (MOFs) for electronic, optoelectronic, and energy storage applications, we present a dataset of predicted electronic structure properties for thousands of MOFs carried out using multiple density functional approximations. Compared to more accurate hybrid functionals, we find that the widely used PBE generalized gradient approximation (GGA) functional severely underpredicts MOF band gaps in a largely systematic manner for semi-conductors and insulators without magnetic character. However, an even larger and less predictable disparity in the band gap prediction is present for MOFs with open-shell 3 d transition metal cations. With regards to partial atomic charges, we find that different density functional approximations predict similar charges overall, although hybrid functionals tend to shift electron density away from the metal centers and onto the ligand environments compared to the GGA point of reference. Much more significant differences in partial atomic charges are observed when comparing different charge partitioning schemes. We conclude by using the dataset of computed MOF properties to train machine-learning models that can rapidly predict MOF band gaps for all four density functional approximations considered in this work, paving the way for future high-throughput screening studies. To encourage exploration and reuse of the theoretical calculations presented in this work, the curated data is made publicly available via an interactive and user-friendly web application on the Materials Project.

36 MATERIALS SCIENCE↗