Global and local curvature in density functional theory
Abstract
Piecewise linearity of the energy with respect to fractional electron removal or addition is a requirement of an electronic structure method that necessitates the presence of a derivative discontinuity at integer electron occupation. Semilocal exchangecorrelation (xc) approximations within density functional theory (DFT) fail to reproduce this behavior, giving rise to deviations from linearity with a convex global curvature that is evidence of manyelectron, selfinteraction error and electron delocalization. Popular functional tuning strategies focus on reproducing piecewise linearity, especially to improve predictions of optical properties. In a divergent approach, Hubbard Uaugmented DFT (i.e., DFT+U) treats selfinteraction errors by reducing the local curvature of the energy with respect to electron removal or addition from one localized subshell to the surrounding system. Although it has been suggested that DFT+U should simultaneously alleviate global and local curvature in the atomic limit, no detailed study on real systems has been carried out to probe the validity of this statement. In this work, we show when DFT+U should minimize deviations from linearity and demonstrate that a “+U” correction will never worsen the deviation from linearity of the underlying xc approximation. However, we explain varying degrees of efficiency of the approach over 27 octahedral transition metal complexesmore »
 Authors:
 Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 (United States)
 (United States)
 Publication Date:
 OSTI Identifier:
 22678941
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Chemical Physics; Journal Volume: 145; Journal Issue: 5; Other Information: (c) 2016 Author(s); Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 37 INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL CHEMISTRY; AMMONIA; DENSITY FUNCTIONAL METHOD; ELECTRONIC STRUCTURE; ERRORS; EXPERIMENTAL DATA; OPTICAL PROPERTIES; TRANSITION ELEMENTS
Citation Formats
Zhao, Qing, Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, Ioannidis, Efthymios I., and Kulik, Heather J., Email: hjkulik@mit.edu. Global and local curvature in density functional theory. United States: N. p., 2016.
Web. doi:10.1063/1.4959882.
Zhao, Qing, Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, Ioannidis, Efthymios I., & Kulik, Heather J., Email: hjkulik@mit.edu. Global and local curvature in density functional theory. United States. doi:10.1063/1.4959882.
Zhao, Qing, Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, Ioannidis, Efthymios I., and Kulik, Heather J., Email: hjkulik@mit.edu. Sun .
"Global and local curvature in density functional theory". United States.
doi:10.1063/1.4959882.
@article{osti_22678941,
title = {Global and local curvature in density functional theory},
author = {Zhao, Qing and Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 and Ioannidis, Efthymios I. and Kulik, Heather J., Email: hjkulik@mit.edu},
abstractNote = {Piecewise linearity of the energy with respect to fractional electron removal or addition is a requirement of an electronic structure method that necessitates the presence of a derivative discontinuity at integer electron occupation. Semilocal exchangecorrelation (xc) approximations within density functional theory (DFT) fail to reproduce this behavior, giving rise to deviations from linearity with a convex global curvature that is evidence of manyelectron, selfinteraction error and electron delocalization. Popular functional tuning strategies focus on reproducing piecewise linearity, especially to improve predictions of optical properties. In a divergent approach, Hubbard Uaugmented DFT (i.e., DFT+U) treats selfinteraction errors by reducing the local curvature of the energy with respect to electron removal or addition from one localized subshell to the surrounding system. Although it has been suggested that DFT+U should simultaneously alleviate global and local curvature in the atomic limit, no detailed study on real systems has been carried out to probe the validity of this statement. In this work, we show when DFT+U should minimize deviations from linearity and demonstrate that a “+U” correction will never worsen the deviation from linearity of the underlying xc approximation. However, we explain varying degrees of efficiency of the approach over 27 octahedral transition metal complexes with respect to transition metal (Sc–Cu) and ligand strength (CO, NH{sub 3}, and H{sub 2}O) and investigate select pathological cases where the delocalization error is invisible to DFT+U within an atomic projection framework. Finally, we demonstrate that the global and local curvatures represent different quantities that show opposing behavior with increasing ligand field strength, and we identify where these two may still coincide.},
doi = {10.1063/1.4959882},
journal = {Journal of Chemical Physics},
number = 5,
volume = 145,
place = {United States},
year = {Sun Aug 07 00:00:00 EDT 2016},
month = {Sun Aug 07 00:00:00 EDT 2016}
}

Using vanderWaalscorrected density functional theory calculations, we explore the possibility of engineering the local structure and morphology of highsurfacearea graphenederived materials to improve the uptake of methane and carbon dioxide for gas storage and sensing. We test the sensitivity of the gas adsorption energy to the introduction of native point defects, curvature, and the application of strain. The binding energy at topological point defect sites is inversely correlated with the number of missing carbon atoms, causing Stone–Wales defects to show the largest enhancement with respect to pristine graphene (~20%). Improvements of similar magnitude are observed at concavely curved surfaces inmore »Cited by 24

Empirical ProtonNeutron Interactions and Nuclear Density Functional Theory: Global, Regional, and Local Comparisons
Calculations of nuclear masses, using nuclear density functional theory, are presented for eveneven nuclei spanning the nuclear chart. The resulting binding energy differences can be interpreted in terms of valence protonneutron interactions. These are compared globally, regionally, and locally with empirical values. Overall, excellent agreement is obtained. Discrepancies highlight neglected degrees of freedom and can point to improved density functionals. 
Empirical ProtonNeutron Interactions and Nuclear Density Functional Theory: Global, Regional, and Local Comparisons
Calculations of nuclear masses, using nuclear density functional theory, are presented for eveneven nuclei spanning the nuclear chart. The resulting binding energy differences can be interpreted in terms of valence protonneutron interactions. These are compared globally, regionally, and locally with empirical values. Overall, excellent agreement is obtained. Discrepancies highlight neglected degrees of freedom and can point to improved density functionals. 
Curvature and Frontier Orbital Energies in Density Functional Theory
Perdew et al. [Phys. Rev. Lett 49, 1691 (1982)] discovered and proved two different properties that exact KohnSham density functional theory (DFT) must obey: (i) The exact total energy versus particle number must be a series of linear segments between integer electron points; (ii) Across an integer number of electrons, the exchangecorrelation potential may ``jump’’ by a constant, known as the derivative discontinuity (DD). Here, we show analytically that in both the original and the generalized KohnSham formulation of DFT, the two are in fact two sides of the same coin. Absence of a derivative discontinuity necessitates deviation from piecewisemore » 
Localdensity approximation, hierarchy of equations, functional expansion, and adiabatic connection in currentdensityfunctional theory
Identities and approximations are proposed for the currentdensityfunctional theory. Based on the virial theorem and the scaling properties of the kinetic and exchange energy density functionals, when local and variableseparation assumptions are made, local formulas for these functionals are obtained. One also obtains hierarchies of equations, which combine to give series expansions of the functionals in terms of their functional derivatives with respect to both the density and the paramagnetic current density. In addition, a general formulation is investigated in the form of the constrained search of a coupled Hamiltonian. Subject to the validity of a Taylor series, both themore »