Augmented potential, energy densities, and virial relations in the weak- and strong-interaction limits of DFT

Large-Z atoms in the strong-interaction limit of DFT: Implications for gradient expansions and for the Lieb-Oxford bound

We numerically study the strong-interaction limit of the exchange-correlation functional for neutral atoms and Bohr atoms as the number of electrons increases. Using a compact representation, we analyze the second-order gradient expansion, comparing it with the one for exchange (weak interaction limit). The two gradient expansions, at strong and weak interaction, turn out to be very similar in magnitude but with opposite signs. We find that the point-charge plus continuum model is surprisingly accurate for the gradient expansion coefficient at strong coupling, while generalized gradient approximations, such as Perdew-Burke-Ernzerhof (PBE) and PBEsol, severely underestimate it. We then use our results to analyze the Lieb-Oxford bound from the point of view of slowly varying densities, clarifying some aspects on the bound at a fixed number of electrons.

Nonlocal Functionals Inspired by the Strongly Interacting Limit of DFT: Exact Constraints and Implementation

Capturing strong correlation effects remains a key challenge for the development of improved exchange-correlation (XC) functionals in density functional theory. The recently proposed multiple radii functional (MRF) [J. Phys. Chem. Lett. 2017, 8, 2799; J. Chem. Theory Comput. 2019, 15, 3580] was designed to capture strong correlation effects seamlessly, as its mathematical structure draws from that of the exact XC functional in the limit of infinite correlations. The MRF functional provides a framework for building approximations along the density-fixed adiabatic connection, delivers accurate XC energy densities in the standard DFT gauge (same as that of the exact exchange energy density), and is free of one-electron self-interaction errors. To facilitate the development of XC functionals based on the MRF, we examine the behavior of the MRF functional when applied to uniform and scaled densities and consider how it can be made exact for the uniform electron gas. These theoretical insights are then used to build improved forms for the fluctuation function, an object that defines XC energy densities within the MRF framework. We also show how the MRF fluctuation function for physical correlation can be easily readjusted to accurately capture the XC functional in the limit of infinite correlations, demonstrating the versatility of MRF for building approximations for different correlation regimes. We describe the implementation of MRF using densities expanded on Gaussian basis sets, which improves the efficiency of previous grid-based MRF implementations.

A Density Functional Theory for the Average Electron Energy

A formally exact density functional theory (DFT) determination of the average electron energy is presented. Our theory, which is based on a different accounting of energy functional terms, partially solves one well-known downside of conventional Kohn-Sham (KS) DFT: that electronic energies have but tenuous connections to physical quantities. Calculated average electron energies are close to experimental ionization potentials (IPs) in one-electron systems, demonstrating a surprisingly small effect of self-interaction and other exchange-correlation errors in established DFT methods. Remarkable agreement with ab initio quantum mechanical calculations of multielectron systems is demonstrated using several flavors of DFT, and we argue for the use of the average electron energy as a design criterion for density functional approximations.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Gradient Expansions for the Large-Coupling Strength Limit of the Møller-Plesset Adiabatic Connection

The adiabatic connection that has, as weak-interaction expansion, the Møller-Plesset perturbation series has been recently shown to have a large coupling-strength expansion, in terms of functionals of the Hartree-Fock density with a clear physical meaning. In this work, we accurately evaluate these density functionals and we extract second-order gradient coefficients from the data for neutral atoms, following ideas similar to the ones used in the literature for exchange, with some modifications. These new gradient expansions will be the key ingredient for performing interpolations that have already been shown to reduce dramatically MP2 errors for large noncovalent complexes. As a byproduct, our investigation of neutral atoms with large number of electrons N indicates that the second-order gradient expansion for exchange grows as N log(N) rather than as N, as often reported in the literature.

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.

MAP: An MP2 Accuracy Predictor for Weak Interactions from Adiabatic Connection Theory

Second-order Møller-Plesset perturbation theory (MP2) approximates the exact Hartree-Fock (HF) adiabatic connection (AC) curve by a straight line. Thus, by using the deviation of the exact curve from the linear behavior, we construct an indicator for the accuracy of MP2. We then use an interpolation along the HF AC to transform the exact form of our indicator into a highly practical MP2 accuracy predictor (MAP) that comes at a negligible additional computational cost. We show that this indicator is already applicable to systems that dissociate into fragments with a nondegenerate ground state, and we illustrate its usefulness by applying it to the S22 and S66 datasets.

Exchange-Correlation Energy Densities and Response Potentials: Connection between Two Definitions and Analytical Model for the Strong-Coupling Limit of a Stretched Bond

We analyze in depth two widely used definitions (from the theory of conditional probability amplitudes and from the adiabatic connection formalism) of the exchange-correlation energy density and of the response potential of Kohn–Sham density functional theory. We introduce a local form of the coupling-constant-dependent Hohenberg–Kohn functional, showing that the difference between the two definitions is due to a corresponding local first-order term in the coupling constant, which disappears globally (when integrated over all space), but not locally. We also design an analytic representation for the response potential in the strong-coupling limit of density functional theory for a model single stretched bond.

Kinetic Correlation Functionals from the Entropic Regularization of the Strictly Correlated Electrons Problem

In this work, we study the entropic regularization of the strictly correlated electrons formalism, discussing the implications for density functional theory and establishing a link with earlier works on quantum kinetic energy and classical entropy. We carry out a very preliminary investigation (using simplified models) on the use of the solution of the entropic regularized problem to build approximations for the kinetic correlation functional at large coupling strengths. We also analyze lower and upper bounds to the Hohenberg-Kohn functional using the entropic regularized strictly correlated electrons problem.

Local and global interpolations along the adiabatic connection of DFT: a study at different correlation regimes

Interpolating the exchange-correlation energy along the density-fixed adiabatic connection of density functional theory is a promising way to build approximations that are not biased toward the weakly correlated regime. These interpolations can be performed at the global (integrated over all spaces) or at the local level, using energy densities. Many features of the relevant energy densities as well as several different ways to construct these interpolations, including comparisons between global and local variants, are investigated here for the analytically solvable Hooke's atom series, which allows for an exploration of different correlation regimes. We also analyze different ways to define the correlation kinetic energy density, focusing on the peak in the kinetic correlation potential.

Response Potential in the Strong-Interaction Limit of Density Functional Theory: Analysis and Comparison with the Coupling-Constant Average

Using the formalism of the conditional amplitude, we study the response part of the exchange-correlation potential in the strong-coupling limit of density functional theory, analyzing its peculiar features and comparing it with the response potential averaged over the coupling constant for small atoms and for the hydrogen molecule. We also use a simple one-dimensional model of a stretched heteronuclear molecule to derive exact properties of the response potential in the strong-coupling limit. The simplicity of the model allows us to unveil relevant features also of the exact Kohn-Sham potential and its different components, namely the appearance of a second peak in the correlation kinetic potential on the side of the most electronegative atom.

Restoring Size Consistency of Approximate Functionals Constructed from the Adiabatic Connection

Approximate exchange-correlation functionals built by modeling in a nonlinear way the adiabatic connection (AC) integrand of density functional theory have many attractive features, being virtually parameter-free and satisfying different exact properties, but they also have a fundamental flaw: they violate the size-consistency condition, crucial to evaluate interaction energies of molecular systems. We show that size consistency in the AC-based functionals can be restored in a very simple way at no extra computational cost. Results on a large set of benchmark molecular interaction energies show that functionals based on the interaction strength interpolation approximations are significantly more accurate than second-order perturbation theory.

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.