No Arabic abstract
Modern density functional theory (DFT) calculations employ the Kohn-Sham (KS) system of non-interacting electrons as a reference, with all complications buried in the exchange-correlation energy (Exc). The adiabatic connection formula gives an exact expression for Exc. We consider DFT calculations that instead employ a reference of strictly-correlated electrons. We define a decorrelation energy that relates this reference to the real system, and derive the corresponding adiabatic connection formula. We illustrate this theory in three situations, namely the uniform electron gas, Hookes atom, and the stretched hydrogen molecule. The adiabatic connection for strictly-correlated electrons provides an alternative perspective for understanding density functional theory and constructing approximate functionals.
We investigate a number of formal properties of the adiabatic strictly-correlated electrons (SCE) functional, relevant for time-dependent potentials and for kernels in linear response time-dependent density functional theory. Among the former, we focus on the compliance to constraints of exact many-body theories, such as the generalised translational invariance and the zero-force theorem. Within the latter, we derive an analytical expression for the adiabatic SCE Hartree exchange-correlation kernel in one dimensional systems, and we compute it numerically for a variety of model densities. We analyse the non-local features of this kernel, particularly the ones that are relevant in tackling problems where kernels derived from local or semi-local functionals are known to fail.
In density functional theory (DFT), the exchange-correlation functional can be exactly expressed by the adiabatic connection integral. It has been noticed that as lambda goes to infinity, the lambda^(-1) term in the expansion of W(lambda) vanishes. We provide a simple but rigorous derivation to this exact condition in this work. We propose a simple parametric form for the integrand, satisfying this condition, and show that it is highly accurate for weakly-correlated two-electron systems.
The strong-interaction limit of the Hohenberg-Kohn functional defines a multimarginal optimal transport problem with Coulomb cost. From physical arguments, the solution of this limit is expected to yield strictly-correlated particle positions, related to each other by co-motion functions (or optimal maps), but the existence of such a deterministic solution in the general three-dimensional case is still an open question. A conjecture for the co-motion functions for radially symmetric densities was presented in Phys.~Rev.~A {bf 75}, 042511 (2007), and later used to build approximate exchange-correlation functionals for electrons confined in low-density quantum dots. Colombo and Stra [Math.~Models Methods Appl.~Sci., {bf 26} 1025 (2016)] have recently shown that these conjectured maps are not always optimal. Here we revisit the whole issue both from the formal and numerical point of view, finding that even if the conjectured maps are not always optimal, they still yield an interaction energy (cost) that is numerically very close to the true minimum. We also prove that the functional built from the conjectured maps has the expected functional derivative also when they are not optimal.
We present a new open-source program, DCore, that implements dynamical mean-field theory (DMFT). DCore features a user-friendly interface based on text and HDF5 files. It allows DMFT calculations of tight-binding models to be performed on predefined lattices as well as textit{ab initio} models constructed by external density functional theory codes through the Wannier90 package. Furthermore, DCore provides interfaces to many advanced quantum impurity solvers such as quantum Monte Carlo and exact diagonalization solvers. This paper details the structure and usage of DCore and shows some applications.
Starting from the (Hubbard) model of an atom, we demonstrate that the uniqueness of the mapping from the interacting to the noninteracting Greens function, $Gto G_0$, is strongly violated, by providing numerous explicit examples of different $G_0$ leading to the same physical $G$. We argue that there are indeed infinitely many such $G_0$, with numerous crossings with the physical solution. We show that this rich functional structure is directly related to the divergence of certain classes of (irreducible vertex) diagrams, with important consequences for traditional many-body physics based on diagrammatic expansions. Physically, we ascribe the onset of these highly non-perturbative manifestations to the progressive suppression of the charge susceptibility induced by the formation of local magnetic moments and/or RVB states in strongly correlated electron systems.