Stefan Vuckovic

Hydrogen Molecule Dissociation Curve with Functionals Based on the Strictly Correlated Regime.

Using the dual Kantorovich formulation, we compute the strictly correlated electrons (SCE) functional (corresponding to the exact strong-interaction limit of density functional theory) for the hydrogen molecule along the dissociation curve. We use an exact relation between the Kantorovich potential and the optimal map to compute the comotion function, exploring corrections based on it. In particular, we analyze how the SCE functional transforms in an exact way the electron-electron distance into a one-body quantity, a feature that can be exploited to build new approximate functionals. We also show that the dual Kantorovich formulation provides in a natural way the constant in the Kohn-Sham potential recently introduced by Levy and Zahariev [Phys. Rev. Lett. 2014, 113, 113002] for finite systems.

Challenging the Lieb–Oxford bound in a systematic way

ABSTRACT The Lieb–Oxford bound, a nontrivial inequality for the indirect part of the many-body Coulomb repulsion in an electronic system, plays an important role in the construction of approximations in density functional theory. Using the wave function for strictly correlated electrons of a given density, we turn the search over wave functions appearing in the original bound into a more manageable search over electron densities. This allows us to challenge the bound in a systematic way. We find that a maximising density for the bound, if it exists, must have compact support. We also find that, at least for particle numbers N ≤ 60, a uniform density profile is not the most challenging for the bound. With our construction, we improve the bound for N = 2 electrons that was originally found by Lieb and Oxford, we give a new lower bound to the constant appearing in the Lieb–Oxford inequality valid for any N, and we provide an improved upper bound for the low-density uniform electron gas indirect energy.

Exchange–Correlation Functionals via Local Interpolation along the Adiabatic Connection

The construction of density-functional approximations is explored by modeling the adiabatic connection locally, using energy densities defined in terms of the electrostatic potential of the exchange-correlation hole. These local models are more amenable to the construction of size-consistent approximations than their global counterparts. In this work we use accurate input local ingredients to assess the accuracy of a range of local interpolation models against accurate exchange-correlation energy densities. The importance of the strictly correlated electrons (SCE) functional describing the strong coupling limit is emphasized, enabling the corresponding interpolated functionals to treat strong correlation effects. In addition to exploring the performance of such models numerically for the helium and beryllium isoelectronic series and the dissociation of the hydrogen molecule, an approximate analytic model is presented for the initial slope of the local adiabatic connection. Comparisons are made with approaches based on global models, and prospects for future approximations based on the local adiabatic connection are discussed.

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.

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

The augmented potential introduced by Levy and Zahariev [Phys. Rev. Lett. 113, 113002 (2014)] is shifted with respect to the standard exchange-correlation potential of the Kohn-Sham density functional theory by a density-dependent constant that makes the total energy become equal to the sum of the occupied orbital energies. In this work, we analyze several features of this approach, focusing on the limit of infinite coupling strength and studying the shift and the corresponding energy density at different correlation regimes. We present and discuss coordinate scaling properties of the augmented potential, study its connection to the response potential, and use the shift to analyze the classical jellium and uniform gas models. We also study other definitions of the energy densities in relation to the functional construction by local interpolations along the adiabatic connection. Our findings indicate that the energy density that is defined in terms of the electrostatic potential of the exchange-correlation hole is particularly well suited for this purpose.

Interpolated energy densities, correlation indicators and lower bounds from approximations to the strong coupling limit of DFT

We investigate the construction of approximated exchange-correlation functionals by interpolating locally along the adiabatic connection between the weak- and the strong-coupling regimes, focussing on the effect of using approximate functionals for the strong-coupling energy densities. The gauge problem is avoided by dealing with quantities that are all locally defined in the same way. Using exact ingredients at weak coupling we are able to isolate the error coming from the approximations at strong coupling only. We find that the nonlocal radius model, which retains some of the non-locality of the exact strong-coupling regime, yields very satisfactory results. We also use interpolation models and quantities from the weak- and strong-coupling regimes to define a correlation-type indicator and a lower bound to the exact exchange-correlation energy. Open problems, related to the nature of the local and global slope of the adiabatic connection at weak coupling, are also discussed.

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.

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.