The shell model for the exchange-correlation hole in the strong-correlation limit

Natural determinant reference functional theory

The natural determinant reference (NDR) or principal natural determinant is the Slater determinant comprised of the N most strongly occupied natural orbitals of an N-electron state of interest. Unlike the Kohn-Sham (KS) determinant, which yields the exact ground-state density, the NDR only yields the best idempotent approximation to the interacting one-particle reduced density matrix, but it is well-defined in common atom-centered basis sets and is representation-invariant. We show that the under-determination problem of prior attempts to define a ground-state energy functional of the NDR is overcome in a grand-canonical ensemble framework at the zero-temperature limit. The resulting grand potential functional of the NDR ensemble affords the variational determination of the ground state energy, its NDR (ensemble), and select ionization potentials and electron affinities. The NDR functional theory can be viewed as an "exactification" of orbital optimization and empirical generalized KS methods. NDR functionals depending on the noninteracting Hamiltonian do not require troublesome KS-inversion or optimized effective potentials.

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.

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.

The factorization ansatz for non-local approximations to the exchange-correlation hole

Among the various types of approximations to the exchange-correlation energy (E_XC), the completely non-local approach is one of the lesser explored approximation schemes. It has not yet reached the predictive power of the widely used generalized gradient approximations, meta-generalized gradient approximations, hybrids, etc. In non-local functionals pursued here, the electron density at every point in space is employed to express the exchange-correlation energy per particle ϵ_XC(r) at a given position r. Here, we use the non-local, spherical-averaged density ρ(r,u)=∫dΩ_u4πρ(r+u) as a starting point to construct approximate exchange-correlation holes through the factorization ansatz ρ_XC(r, u) = f(r, u)ρ(r, u). We present upper and lower bounds to the exchange energy per particle ϵ_X(r) in terms of ρ(r, u). The factor f(r, u) is then designed to satisfy various conditions that represent important exchange and correlation effects. We assess the resulting approximations and find that the complex, oscillatory structure of ρ(r, u) makes the construction of a corresponding f(r, u) very challenging. This conclusion, identifying the main issue of the non-local approximation, is supported by a detailed analysis of the resulting exchange-correlation holes.

Noncovalent Interactions from Models for the Møller-Plesset Adiabatic Connection

Given the omnipresence of noncovalent interactions (NCIs), their accurate simulations are of crucial importance across various scientific disciplines. Here we construct accurate models for the description of NCIs by an interpolation along the Møller-Plesset adiabatic connection (MP AC). Our interpolation approximates the correlation energy, by recovering MP2 at small coupling strengths and the correct large-coupling strength expansion of the MP AC, recently shown to be a functional of the Hartree-Fock density. Our models are size consistent for fragments with nondegenerate ground states, have the same cost as double hybrids, and require no dispersion corrections to capture NCIs accurately. These interpolations greatly reduce large MP2 errors for typical π-stacking complexes (e.g., benzene-pyridine dimers) and for the L7 data set. They are also competitive with state-of-the-art dispersion enhanced functionals and can even significantly outperform them for a variety of data sets, such as CT7 and L7.

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.

Screened range-separated hybrid by balancing the compact and slowly varying density regimes: Satisfaction of local density linear response

Due to their quantitative accuracy and ability to solve several difficulties, screened range-separated hybrid exchange-correlation functionals are now a standard approach for ab initio simulation of condensed matter systems. However, the screened range-separated hybrid functionals proposed so far are biased either toward compact or slowly varying densities. In this paper, we propose a screened range-separated hybrid functional, named HSEint, which can well describe these density regimes, achieving good accuracy for both molecular and solid-state systems. The semilocal part of the proposed functional is based on the PBEint generalized gradient approximation [E. Fabiano et al., Phys. Rev. B 82, 113104 (2010)], constructed for hybrid interfaces. To improve the functional performance, we employ exact or nearly exact constraints in the construction of range-separated hybrid functional, such as recovering of the local density linear response and semiclassical atom linear response.

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.

Investigation of the Exchange-Correlation Potentials of Functionals Based on the Adiabatic Connection Interpolation

We have studied the correlation potentials produced by various adiabatic connection models (ACMs) for several atoms and molecules. The results have been compared to accurate reference potentials (coupled cluster and quantum Monte Carlo results) as well as to state-of-the-art ab initio DFT approaches. We have found that all the ACMs yield correlation potentials that exhibit a correct behavior, quite resembling scaled second-order Görling-Levy (GL2) potentials and including most of the physically meaningful features of the accurate reference data. The behavior and contribution of the strong-interaction limit potentials have also been investigated and discussed.

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.

Fermionic Statistics in the Strongly Correlated Limit of Density Functional Theory

Exact pieces of information on the adiabatic connection integrand, W_λ[ρ], which allows evaluation of the exchange-correlation energy of Kohn-Sham density functional theory, can be extracted from the leading terms in the strong coupling limit (λ → ∞, where λ is the strength of the electron-electron interaction). In this work, we first compare the theoretical prediction for the two leading terms in the strong coupling limit with data obtained via numerical implementation of the exact Levy functional in the simple case of two electrons confined in one dimension, confirming the asymptotic exactness of these two terms. We then carry out a first study on the incorporation of the Fermionic statistics at large coupling λ, both numerical and theoretical, confirming that spin effects enter at orders ∼e^-√λ.

Simple Fully Nonlocal Density Functionals for Electronic Repulsion Energy

From a simplified version of the mathematical structure of the strong coupling limit of the exact exchange-correlation functional, we construct an approximation for the electronic repulsion energy at physical coupling strength, which is fully nonlocal. This functional is self-interaction free and yields energy densities within the definition of the electrostatic potential of the exchange-correlation hole that are locally accurate and have the correct asymptotic behavior. The model is able to capture strong correlation effects that arise from chemical bond dissociation, without relying on error cancellation. These features, which are usually missed by standard density functional theory (DFT) functionals, are captured by the highly nonlocal structure, which goes beyond the "Jacob's ladder" framework for functional construction, by using integrals of the density as the key ingredient. Possible routes for obtaining the full exchange-correlation functional by recovering the missing kinetic component of the correlation energy are also implemented and discussed.