The exact interaction energy of a many-electron system is determined by the electron pair density, which is not well-approximated in standard Kohn--Sham density functional models. Here we study the (complicated but well-defined) exact universal map from density to pair density. We survey how many common functionals, including the most basic version of the local density approximation (Dirac exchange with no correlation contribution), arise from particular approximations of this map. We develop an algorithm to compute the map numerically and apply it to one-parameter families
{αρ(αx)}α>0 of one-dimensional homogeneous and inhomogeneous single-particle densities. We observe that the pair density develops remarkable multiscale patterns which strongly depend on both the particle number and the “width”
α−1 of the single-particle density. The simulation results are confirmed by rigorous asymptotic results in the limiting regimes
α≫ and
α≪1. For one-dimensional homogeneous systems, we show that the whole spectrum of patterns is reproduced surprisingly well by a simple asymptotics-based ansatz which slowly smoothens out the “strictly correlated”
α=0 pair density while slowly turning on the
α=∞ “exchange” terms as
α increases. Our findings lend theoretical support to the celebrated semi-empirical idea [A. D. Becke, J. Chem. Phys., 98 (1993), pp. 5648--5652] to mix in a fractional amount of exchange, albeit not to assuming the mixing to be additive and taking the fraction to be a system-independent constant.
In pair density-functional theory, the only unknown piece of the energy is the kinetic energy
T as a functional of the pair density
P(x1,x2). Although
T[P] has a simpler structure than the Hohenberg-Kohn functional of conventional density-functional theory, computational requirements are still moderate. In the present work, a particularly convenient model system to represent many-electron pair densities is introduced. This ``boson pair model'' (BPM) approximately treats electron pairs as noninteracting bosons. The resulting explicit model for the kinetic energy
T2[P] is shown to be exact for two-electron systems and a lower bound to
T[P] for more than two electrons. The one- and two-particle density matrices obtained from the BPM yield upper bounds for the corresponding many-electron quantities. This suggests a partitioning
T[P]=T2[P]+Teff[P], where only the remainder ${T}_{\mathrm{eff}}\phantom{\rule{0.1em}{0ex}}[P]\ensuremath{\geqslant}0$ needs to be approximated. If the BPM is constrained to yield the exact ground-state pair density, a two-electron Schr\"odinger equation with an effective local two-particle potential results; the latter is identified as a sum of the bare Coulomb interaction and the functional derivative of
Teff[P]. This self-consistent scheme to minimize the energy with respect to
P is more efficient than previous procedures. Further information on the functional derivative of
Teff[P] is derived from a contracted Schr\"odinger equation. Since
Teff[P] is explicitly known in the two-electron and noninteracting (Hartree-Fock) limits, the present method provides an alternative to density-matrix functional theories, which can be exact in the same limits and are similar in computational cost.
Several explicit formulas for the kinetic energy of a many-electron system as a functional of the k-electron density are derived, with emphasis on the electron pair density. The emphasis is on general techniques for deriving approximate kinetic energy functionals and features generalized Weisacker bounds and methods using density-matrix reconstruction. Adapting results from statistical mechanics, a hierarchy of equations is derived that links electron pairs, triplets, quadruplets, etc.; this may be used to derive more accurate approximations. Several methods for defining the exact kinetic energy functional are presented, including the generalizations of the Levy and Lieb formulations of density-functional theory. Together with N-representability constraints on the k-density, this paper provides the basis for “generalized density functional theories” based on the electron pair density. There are also implications for conventional density-functional theory, notably regarding the development of more accurate density functionals for the kinetic energy.
A density-functional theory (DFT) of lattice fermion models is presented, which uses the single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ as basic variable. A simple, explicit approximation to the interaction-energy functional $W[\ensuremath{\gamma}]$ of the Hubbard model is derived from exact dimer results, scaling properties of $W[\ensuremath{\gamma}]$ and known limits. Systematic tests on the one-dimensional chain show a remarkable agreement with the Bethe-ansatz exact solution for all interaction regimes and band fillings. New results are obtained for the ground-state energy and charge-excitation gap in two dimensions. A successful description of strong electron correlations within DFT is achieved.
We provide a brief overview of approaches for calculating the density of states of quantum systems and random matrix Hamiltonians using the tools of free probability theory. For a given Hamiltonian of a quantum system or a generic random matrix Hamiltonian, which can be written as a sum of two non-commutating operators, one can obtain an expression for the density of states of the Hamiltonian from the known density of states of the two component operators by assuming that these operators are mutually free and by using the free additive convolution. In many examples of interacting quantum systems and random matrix models, this procedure is known to provide a reasonably accurate approximation to the exact numerical density of states. We review some of the examples that are known in the literature where this procedure works very well, and also discuss some of the limitations of this method in situations where the free probability approximation fails to provide a sufficiently accurate description of the exact density of states. Subsequently, we describe a perturbation scheme that can be developed from the subordination formulas for the Cauchy transform of the density of states and use it to obtain approximate analytical expressions for the density of states in various models, such as the Rosenzweig-Porter random matrix ensemble and the Anderson model with on-site disorder.
Density functional theory (DFT) is the de facto method for the electronic structure of weakly correlated systems. But for strongly correlated materials, common density functional approximations break down. Here, we derive a many-electron expansion (MEE) in DFT that accounts for successive one-, two-, three-, ... particle interactions within the system. To compute the correction terms, the density is first decomposed into a sum of localized, nodeless one-electron densities (${\ensuremath{\rho}}_{i}$). These one-electron densities are used to construct relevant two- (${\ensuremath{\rho}}_{i}+{\ensuremath{\rho}}_{j}$), three- (${\ensuremath{\rho}}_{i}+{\ensuremath{\rho}}_{j}+{\ensuremath{\rho}}_{k}$), ... electron densities. Numerically exact results for these few-particle densities can then be used to correct an approximate density functional via any of several many-body expansions. We show that the resulting hierarchy gives accurate results for several important model systems: the Hubbard and Peierls-Hubbard models in 1D and the pure Hubbard model in 2D. We further show that the method is numerically convergent for strongly correlated systems: applying successively higher order corrections leads to systematic improvement of the results. MEE thus provides a hierarchy of density functional approximations that applies to both weakly and strongly correlated systems.