Generalized nonlocal kinetic energy density functionals based on the von Weizsäcker functional

We generalize the ideas behind the procedure for the construction of kinetic energy density functionals with a nonlocal term based on the structure of the von Weizs?cker functional, and present several types of nonlocal terms. In all cases, the functionals are constructed such that they reproduce the linear response function of the homogeneous electron gas. These functionals are designed by rewriting the von Weizs?cker functional with the help of a parameter β that determines the power of the electron density in the expression, a strategy we have previously used in the generalization of Thomas-Fermi nonlocal functionals. Benchmark calculations in localized systems have been performed with these functionals to test both their relative errors and the quality of their local behavior. We have obtained competitive results when compared to semilocal and previous nonlocal functionals, the generalized nonlocal von Weizs?cker functionals giving very good results for the total kinetic energies and improving the local behavior of the kinetic energy density. In addition, all the functionals discussed in this paper, when using an adequate reference density, can be evaluated as a single integral in momentum space, resulting in a quasilinear scaling for the computational cost.     

The AM05 density functional applied to solids

We show that the AM05 functional [Armiento and Mattsson, Phys. Rev. B 72, 085108 (2005)] has the same excellent performance for solids as the hybrid density functionals tested in Paier et al. [J. Chem. Phys. 124, 154709 (2006); 125, 249901 (2006)]. This confirms the original finding that AM05 performs exceptionally well for solids and surfaces. Hartree-Fock hybrid calculations are typically an order of magnitude slower than local or semilocal density functionals such as AM05, which is of a regular semilocal generalized gradient approximation form. The performance of AM05 is on average found to be superior to selecting the best of local density approximation and PBE for each solid. By comparing data from several different electronic-structure codes, we have determined that the numerical errors in this study are equal to or smaller than the corresponding experimental uncertainties.     

Assessment of the pairwise additive approximation and evaluation of many-body terms for water clusters

We have tested the ability of four commonly used density functionals (three of which are semilocal and one of which is nonlocal) to outperform accurate pairwise additive approximations in the prediction of binding energies for a series of water clusters ranging in size from dimer to pentamer. Comparison to results obtained with the Weizmann-1 (W1) level of wave function theory shows that while all density functionals are capable of outperforming the accurate pairwise data, the choice of basis set used is crucial to the performance of the method, and if a poor choice of basis set is made the errors obtained with density functional theory (DFT) can be larger than those obtained with the simple pairwise approximation. We have also compared the binding energies and many-body terms determined with DFT to those obtained with W1, and have found that all semilocal functionals have significant errors in the many-body components of the full interactions energy. Despite this limitation, however, we have found that, of the four functionals tested, PBE1W/MG3S is the most accurate for predicting the binding energies of the clusters.     

Many-electron self-interaction and spin polarization errors in local hybrid density functionals

Errors for systems with noninteger occupation have been connected to common failures of density functionals. Previously, global hybrids and pure density functionals have been investigated for systems with noninteger charge and noninteger spin state. Local hybrids have not been investigated for either of those systems to the best of our knowledge. This study intends to close this gap. We investigate systems with noninteger charge to assess the many-electron self-interaction error and systems with noninteger spin state to assess the spin polarization error of recently proposed local hybrids and their range-separated variants. We find that long-range correction is very important to correct for many-electron self-interaction error in cations, whereas most full-range local hybrids seem to be sufficient for anions, where long-range-corrected density functionals tend to overcorrect. On the other hand, while all hitherto proposed long-range-corrected density functionals show large spin polarization errors, the Perdew-Staroverov-Tao-Scuseria (PSTS) functional performs best of all local hybrids in this case and shows an outstanding behavior for the dependence of the energy on the spin polarization.     

Role and effective treatment of dispersive forces in materials: Polyethylene and graphite crystals as test cases

A semiempirical addition of dispersive forces to conventional density functionals (DFT-D) has been implemented into a pseudopotential plane-wave code. Test calculations on the benzene dimer reproduced the results obtained by using localized basis set, provided that the latter are corrected for the basis set superposition error. By applying the DFT-D/plane-wave approach a substantial agreement with experiments is found for the structure and energetics of polyethylene and graphite, two typical solids that are badly described by standard local and semilocal density functionals.     

Analysis of the Heyd-Scuseria-Ernzerhof density functional parameter space

The Heyd-Scuseria-Ernzerhof (HSE) density functionals are popular for their ability to improve upon the accuracy of standard semilocal functionals such as Perdew-Burke-Ernzerhof (PBE), particularly for semiconductor band gaps. They also have a reduced computational cost compared to hybrid functionals, which results from the restriction of Fock exchange calculations to small inter-electron separations. These functionals are defined by an overall fraction of Fock exchange and a length scale for exchange screening. We systematically examine this two-parameter space to assess the performance of hybrid screened exchange (sX) functionals and to determine a balance between improving accuracy and reducing the screening length, which can further reduce computational costs. Three parameter choices emerge as useful: "sX-PBE" is an approximation to the sX-LDA screened exchange density functionals based on the local density approximation (LDA); "HSE12" minimizes the overall error over all tests performed; and "HSE12s" is a range-minimized functional that matches the overall accuracy of the existing HSE06 parameterization but reduces the Fock exchange length scale by half. Analysis of the error trends over parameter space produces useful guidance for future improvement of density functionals.     

Hyper-generalized-gradient functionals constructed from the Lieb-Oxford bound: implementation via local hybrids and thermochemical assessment

In 2009 Odashima and Capelle (OC) showed a way to design a correlation-only density functional that satisfies a Lieb-Oxford bound on the correlation energy, without empirical parameters and even without additional theoretical parameters. However, they were only able to test a size-inconsistent version of it that employs total energies. Here, we show that their alternative size-consistent form that employs energy densities, when combined with exact or semilocal exchange, is a local hybrid (lh) functional. We test several variants of this nonempirical OC-lh functional on standard molecular test sets. Although no variant yields enthalpies of formation with the accuracy of the semilocal Tao-Perdew-Staroverov-Scuseria (TPSS) exchange-correlation, OC-lh correlation with exact exchange yields rather accurate energy barriers for chemical reactions. Our purpose here is not to advocate for a new density functional, but to explore a previously published idea. We also discuss the importance of near-self-consistency for fully nonlocal functionals.     

Electronic structure of copper phthalocyanine: a comparative density functional theory study

We present a systematic density functional theory study of the electronic structure of copper phthalocyanine (CuPc) using several different (semi)local and hybrid functionals and compare the results to experimental photoemission data. We show that semilocal functionals fail qualitatively for CuPc primarily because of underbinding of localized orbitals due to self-interaction errors. We discuss an appropriate choice of functional for studies of CuPc/metal interfaces and suggest the Heyd-Scuseria-Ernzerhof screened hybrid functional as a suitable compromise functional.     

Why are time-dependent density functional theory excitations in solids equal to band structure energy gaps for semilocal functionals, and how does nonlocal Hartree-Fock-type exchange introduce excitonic effects?

We examine the time-dependent density functional theory (TD-DFT) equations for calculating excitation energies in solids with Gaussian orbitals and analytically show that for semilocal functionals, their lowest eigenvalue collapses to the minimum band orbital energy difference. With the introduction of nonlocal Hartree-Fock-type exchange (as in hybrid functionals), this result is no longer valid, and the lowest TD-DFT eigenvalue reflects the appearance of excitonic effects. Previously reported "charge-transfer" problems with semilocal TD-DFT excitations in molecules can be deduced from our analysis by taking the limit to infinite lattice constant.     

Accurate Kohn–Sham ionization potentials from scaled‐opposite‐spin second‐order optimized effective potential methods

One important property of Kohn–Sham (KS) density functional theory is the exact equality of the energy of the highest occupied KS orbital (HOMO) with the negative ionization potential of the system. This exact feature is out of reach for standard density‐dependent semilocal functionals. Conversely, accurate results can be obtained using orbital‐dependent functionals in the optimized effective potential (OEP) approach. In this article, we investigate the performance, in this context, of some advanced OEP methods, with special emphasis on the recently proposed scaled‐opposite‐spin OEP functional. Moreover, we analyze the impact of the so‐called HOMO condition on the final quality of the HOMO energy. Results are compared to reference data obtained at the CCSD(T) level of theory. © 2016 Wiley Periodicals, Inc.     

Nonspherical model density matrices for Rung 3.5 density functionals

"Rung 3.5" exchange-correlation functionals for Kohn-Sham density functional theory depend linearly on the nonlocal one-particle density matrix of the noninteracting Kohn-Sham reference system. Rung 3.5 functionals also require a semilocal model for the one-particle density matrix. This work presents new model density matrices for Rung 3.5 functionals. The resulting functionals give reasonable predictions for total energies, molecular thermochemistry and kinetics, odd-electron bonds, and conjugated polymer bandgaps. Global-hybrid-like combinations of semilocal and Rung 3.5 exchange, and empirical density matrix models, also show promise.     

The performance of semilocal and hybrid density functionals in 3d transition-metal chemistry

We investigate the performance of contemporary semilocal and hybrid density functionals for bond energetics, structures, dipole moments, and harmonic frequencies of 3d transition-metal (TM) compounds by comparison with gas-phase experiments. Special attention is given to the nonempirical metageneralized gradient approximation (meta-GGA) of Tao, Perdew, Staroverov, and Scuseria (TPSS) [Phys. Rev. Lett. 91, 146401 (2003)], which has been implemented in TURBOMOLE for the present work. Trends and error patterns for classes of homologous compounds are analyzed, including dimers, monohydrides, mononitrides, monoxides, monofluorides, polyatomic oxides and halogenides, carbonyls, and complexes with organic pi ligands such as benzene and cyclopentadienyl. Weakly bound systems such as Ca(2), Mn(2), and Zn(2) are discussed. We propose a reference set of reaction energies for benchmark purposes. Our all-electron results with quadruple zeta valence basis sets validate semilocal density-functional theory as the workhorse of computational TM chemistry. Typical errors in bond energies are substantially larger than in (organic) main group chemistry, however. The Becke-Perdew'86 [Phys. Rev. A 38, 3098 (1988); Phys. Rev. B 33, 8822 (1986)] GGA and the TPSS meta-GGA have the best price/performance ratio, while the TPSS hybrid functional achieves a slightly lower mean absolute error in bond energies. The popular Becke three-parameter hybrid B3LYP underbinds significantly and tends to overestimate bond distances; we give a possible explanation for this. We further show that hybrid mixing does not reduce the width of the error distribution on our reference set. The error of a functional for the s-d transfer energy of a TM atom does not predict its error for TM bond energies and bond lengths. For semilocal functionals, self-interaction error in one- and three-electron bonds appears to be a major source of error in TM reaction energies. Nevertheless, TPSS predicts the correct ground-state symmetry in the vast majority of cases and rarely fails qualitatively. This further confirms TPSS as a general purpose functional that works throughout the periodic table. We also give workstation timing comparisons for the 645-atom protein crambin.     

Classical kinetic energy,quantum fluctuation terms and kinetic-energy functionals

We employ a recently formulated dequantization procedure to obtain an exact expression for the kinetic energy which is applicable to all kinetic-energy functionals. We express the kinetic energy of an N-electron system as the sum of an N-electron classical kinetic energy and an N-electron purely quantum kinetic energy arising from the quantum fluctuations that turn the classical momentum into the quantum momentum. This leads to an interesting analogy with Nelson’s stochastic approach to quantum mechanics, which we use to conceptually clarify the physical nature of part of the kinetic-energy functional in terms of statistical fluctuations and in direct correspondence with Fisher Information Theory. We show that the N-electron purely quantum kinetic energy can be written as the sum of the (one-electron) Weizsäcker term and an (N?1)-electron kinetic correlation term. We further show that the Weizsäcker term results from local fluctuations while the kinetic correlation term results from the nonlocal fluctuations. We then write the N-electron classical kinetic energy as the sum of the (one-electron) classical kinetic energy and another (N?1)-electron kinetic correlation term. For one-electron orbitals (where kinetic correlation is neglected) we obtain an exact (albeit impractical) expression for the noninteracting kinetic energy as the sum of the classical kinetic energy and the Weizsäcker term. The classical kinetic energy is seen to be explicitly dependent on the electron phase, and this has implications for the development of accurate orbital-free kinetic-energy functionals. Also, there is a direct connection between the classical kinetic energy and the angular momentum and, across a row of the periodic table, the classical kinetic energy component of the noninteracting kinetic energy generally increases as Z increases. Finally, we underline that, although our aim in this paper is conceptual rather than practical, our results are potentially useful for the construction of improved kinetic-energy functionals.     

libKEDF: An accelerated library of kinetic energy density functionals

Kinetic energy density functionals (KEDFs) approximate the kinetic energy of a system of electrons directly from its electron density. They are used in electronic structure methods that lack direct access to orbitals, for example, orbital‐free density functional theory (OFDFT) and certain embedding schemes. In this contribution, we introduce libKEDF, an accelerated library of modern KEDF implementations that emphasizes nonlocal KEDFs. We discuss implementation details and assess the performance of the KEDF implementations for large numbers of atoms. We show that using libKEDF, a single computing node or (GPU) accelerator can provide easy computational access to mesoscale chemical and materials science phenomena using OFDFT algorithms. © 2017 Wiley Periodicals, Inc.