Effective homogeneity of the exchange–correlation and non-interacting kinetic energy functionals under density scaling

Effective homogeneity of Fermi-Amaldi-containing exchange-correlation functionals

Parr and Ghosh [Phys. Rev. A. 51 3564 (1995)] demonstrated that when near-exact electron densities and potentials are used, the exchange-correlation energies of first- and second-row atoms are well-described by a combination of the Fermi-Amaldi functional with a functional that is homogeneous of degree one under density scaling. Insight into this observation is provided by considering their work from the perspective of the effective homogeneity of the overall exchange-correlation functional. By considering a general form that combines the Fermi-Amaldi functional with a functional that is homogeneous of degree k, it is shown that for these atoms, the functional of Parr and Ghosh (k = 1) exhibits essentially optimal effective homogeneities on the electron-deficient side of the integer. Percentage errors in effective homogeneities are close to percentage errors in exchange-correlation energies.

Incorporation of the Fermi-Amaldi Term into Direct Energy Kohn-Sham Calculations

In direct energy Kohn-Sham (DEKS) theory, the density functional theory electronic energy equals the sum of occupied orbital energies, obtained from Kohn-Sham-like orbital equations involving a shifted Hartree exchange-correlation potential, which must be approximated. In the present study, the Fermi-Amaldi term is incorporated into approximate DEKS calculations, introducing the required -1/r contribution to the exchange-correlation component of the shifted potential in asymptotic regions. It also provides a mechanism for eliminating one-electron self-interaction error, and it introduces a nonzero exchange-correlation component of the shift in the potential that is of appropriate magnitude. The resulting electronic energies are very sensitive to the methodologies considered, whereas the highest occupied molecular orbital energies and exchange-correlation potentials are much less sensitive and are similar to those obtained from DEKS calculations using a conventional exchange-correlation functional.

Robust All-Electron Optimization in Orbital-Free Density-Functional Theory Using the Trust-Region Image Method

We present a Gaussian-basis implementation of orbital-free density-functional theory (OF-DFT) in which the trust-region image method (TRIM) is used for optimization. This second-order optimization scheme has been constructed to provide benchmark all-electron results with very tight convergence of the particle-number constraint, associated chemical potential, and electron density. It is demonstrated that, by preserving the saddle-point nature of the optimization and simultaneously optimizing the density and chemical potential, an order of magnitude reduction in the number of iterations required for convergence is obtained. The approach is compared and contrasted with a new implementation of the nested optimization scheme put forward by Chan, Cohen, and Handy. Our implementation allows for semilocal kinetic-energy (and exchange-correlation) functionals to be handled self-consistently in all-electron calculations. The all-electron Gaussian-basis setting for these calculations will enable direct comparison with a wide range of standard high-accuracy quantum-chemical methods as well as with Kohn-Sham density-functional theory. We expect that the present implementation will provide a useful tool for analyzing the performance of approximate kinetic-energy functionals in finite systems.

Simple DFT Scheme for Estimating Negative Electron Affinities

A simple density functional theory (DFT) scheme is proposed for estimating negative vertical electron affinities of neutral systems, based on a consideration of the integer discontinuity and density scaling homogeneity. The key feature is the derivation of two system-dependent exchange-correlation functionals, one appropriate for the electron deficient side of the integer and one appropriate for the electron abundant side. The electron affinity is evaluated as a linear combination of frontier orbital energies from self-consistent Kohn-Sham calculations on the neutral system using these functionals. For two assessments comprising a total of 43 molecules, the scheme provides electron affinities that are in good agreement with experimental values and which are an improvement over those from the DFT method of Tozer and De Proft [ J. Phys. Chem. A 2005 , 109 , 8923 ]. The scheme is trivial to implement in any Kohn-Sham program, and the computational cost is that of a series of generalized gradient approximation Kohn-Sham calculations. More generally, the study provides a prescription for performing low-cost, self-consistent Kohn-Sham calculations that yield frontier orbital energies that approximately satisfy the appropriate Koopmans conditions, without the need for exact exchange.

Exploiting Coordinate Scaling Relations To Accelerate Exact Exchange Calculations

Exact exchange is an important constituent in many state-of-the-art approximations to the exchange-correlation (xc) functional of Kohn-Sham DFT. However, its evaluation can be computationally intensive, which can be particularly prohibitive in DFT-based molecular dynamics (MD) simulations, often restricted to semilocal functionals. We derive a scheme based on the formal coordinate scaling properties of the exact xc functional that allows for a substantial reduction of the cost of the evaluation of both the exact exchange energy and potential. We show that within a plane-wave/pseudopotential framework, excellent accuracy is retained, while speed ups from up to ∼6 can be reached. The coordinate scaling thus accelerates hybrid-functional-based first-principles MD simulations by nearly one order of magnitude.

Approximating the Shifted Hartree-Exchange-Correlation Potential in Direct Energy Kohn-Sham Theory

Levy and Zahariev [Phys. Rev. Lett. 113 113002 (2014)] have proposed a new approach for performing density functional theory calculations, termed direct energy Kohn-Sham (DEKS) theory. In this approach, the electronic energy equals the sum of orbital energies, obtained from Kohn-Sham-like orbital equations involving a shifted Hartree-exchange-correlation potential, which must be approximated. In the present study, density scaling homogeneity considerations are used to facilitate DEKS calculations on a series of atoms and molecules, leading to three nonlocal approximations to the shifted potential. The first two rely on preliminary Kohn-Sham calculations using a standard generalized gradient approximation (GGA) exchange-correlation functional and the results illustrate the benefit of describing the dominant Hartree component of the shift exactly. A uniform electron gas analysis is used to eliminate the need for these preliminary Kohn-Sham calculations, leading to a potential with an unconventional form that yields encouraging results, providing strong motivation for further research in DEKS theory.

Charge-constrained auxiliary-density-matrix methods for the Hartree-Fock exchange contribution

Three new variants of the auxiliary-density-matrix method (ADMM) of Guidon, Hutter, and VandeVondele [J. Chem. Theory Comput. 6, 2348 (2010)] are presented with the common feature that they have a simplified constraint compared with the full orthonormality requirement of the earlier ADMM1 method. All ADMM variants are tested for accuracy and performance in all-electron B3LYP calculations with several commonly used basis sets. The effect of the choice of the exchange functional for the ADMM exchange-correction term is also investigated.

Comment on “Scaling properties of information-theoretic quantities in density functional reactivity theory” by C. Rong, T. Lu, P. W. Ayers, P. K. Chattaraj and S. Liu, Phys. Chem. Chem. Phys., 2015, 17, 4977–4988

The linear dependence between the density per particle σ and the electron density ρ facilitates the theoretical study of the N-scaling rules for quantum information functionals and their atomic partitions.

Simultaneous Determination of Structures, Vibrations, and Frontier Orbital Energies from a Self-Consistent Range-Separated Hybrid Functional

A self-consistent optimally tuned range-separated hybrid density functional (scOT-RSH) approach is developed. It can simultaneously predict accurate geometries, vibrational modes, and frontier orbital energies. This is achieved by optimizing the range-separation parameter, γ, to both satisfy the ionization energy theorem and minimize interatomic forces. We benchmark our approach against an established hybrid functional, B3LYP, using the G2 test set. scOT-RSH greatly improves the accuracy of occupied frontier orbital energies, with a mean absolute error (MAE) of only 0.2 eV relative to experimental ionization energies compared to 2.96 eV with B3LYP. Geometries do not change significantly compared to those obtained from B3LYP, with a bond length MAE of 0.012 Å compared to 0.008 Å for B3LYP, and a 6.5% MAE for zero-point energies, slightly larger than that of B3LYP (3.1%). scOT-RSH represents a new paradigm in which accurate geometries and ionization energies can be predicted simultaneously from a single functional approach.

Outer-valence Electron Spectra of Prototypical Aromatic Heterocycles from an Optimally Tuned Range-Separated Hybrid Functional

Density functional theory with optimally tuned range-separated hybrid (OT-RSH) functionals has been recently suggested [Refaely-Abramson et al. Phys. Rev. Lett.2012, 109, 226405] as a nonempirical approach to predict the outer-valence electronic structure of molecules with the same accuracy as many-body perturbation theory. Here, we provide a quantitative evaluation of the OT-RSH approach by examining its performance in predicting the outer-valence electron spectra of several prototypical gas-phase molecules, from aromatic rings (benzene, pyridine, and pyrimidine) to more complex organic systems (terpyrimidinethiol and copper phthalocyanine). For a range up to several electronvolts away from the frontier orbital energies, we find that the outer-valence electronic structure obtained from the OT-RSH method agrees very well (typically within ∼0.1–0.2 eV) with both experimental photoemission and theoretical many-body perturbation theory data in the GW approximation. In particular, we find that with new strategies for an optimal choice of the short-range fraction of Fock exchange, the OT-RSH approach offers a balanced description of localized and delocalized states. We discuss in detail the sole exception found—a high-symmetry orbital, particular to small aromatic rings, which is relatively deep inside the valence state manifold. Overall, the OT-RSH method is an accurate DFT-based method for outer-valence electronic structure prediction for such systems and is of essentially the same level of accuracy as contemporary GW approaches, at a reduced computational cost.

Revisiting the density scaling of the non-interacting kinetic energy

Two different density scaling approaches are compared and their prospects for use in functional development are reviewed.

Excitation Gaps of Finite-Sized Systems from Optimally Tuned Range-Separated Hybrid Functionals

Excitation gaps are of considerable significance in electronic structure theory. Two different gaps are of particular interest. The fundamental gap is defined by charged excitations, as the difference between the first ionization potential and the first electron affinity. The optical gap is defined by a neutral excitation, as the difference between the energies of the lowest dipole-allowed excited state and the ground state. Within many-body perturbation theory, the fundamental gap is the difference between the corresponding lowest quasi-hole and quasi-electron excitation energies, and the optical gap is addressed by including the interaction between a quasi-electron and a quasi-hole. A long-standing challenge has been the attainment of a similar description within density functional theory (DFT), with much debate on whether this is an achievable goal even in principle. Recently, we have constructed and applied a new approach to this problem. Anchored in the rigorous theoretical framework of the generalized Kohn-Sham equation, our method is based on a range-split hybrid functional that uses exact long-range exchange. Its main novel feature is that the range-splitting parameter is not a universal constant but rather is determined from first principles, per system, based on satisfaction of the ionization potential theorem. For finite-sized objects, this DFT approach mimics successfully, to the best of our knowledge for the first time, the quasi-particle picture of many-body theory. Specifically, it allows for the extraction of both the fundamental and the optical gap from one underlying functional, based on the HOMO-LUMO gap of a ground-state DFT calculation and the lowest excitation energy of a linear-response time-dependent DFT calculation, respectively. In particular, it produces the correct optical gap for the difficult case of charge-transfer and charge-transfer-like scenarios, where conventional functionals are known to fail. In this perspective, we overview the formal and practical challenges associated with gap calculations, explain our new approach and how it overcomes previous difficulties, and survey its application to a variety of systems.