Strongly correlated electrons have been the subject of substantial theoretical interest for many years. Most work has focused on obtaining the energy in a low-cost fashion. Here, we show that even methods with good energies can yield significant "delocalization errors" that affect the orbitals and density, leading to large errors in predicting other important properties such as dipole moments. We illustrate this point by comparing existing state-of-art approaches with an accurate exchange correlation functional based on a generalised valence bond ansatz, in which orbitals and fractional occupations are treated as variational parameters via a common optimized effective potential (OEP). We show that the OEP exhibits step and peak features which, similar to the exact Kohn-Sham (KS) potential of DFT, are crucial to prevent charge delocalization. We further show that the step is missing in common approximations within reduced density matrix functional theory resulting in delocalization errors comparable to those found in DFT approximations. Finally, we explain the delocalisation error as coming from an artificial mixing of the ground state with a charge-transfer excited state which is avoided if occupation numbers exhibit discontinuities.
Strongly correlated electrons have been the subject of substantial theoretical interest for many years. Most work has focused on obtaining the energy in a low-cost fashion. Here, we show that even methods with good energies can yield significant "delocalization errors" that affect the orbitals and density, leading to large errors in predicting other important properties such as dipole moments. We illustrate this point by comparing existing state-of-art approaches with an accurate exchange correlation functional based on a generalised valence bond ansatz, in which orbitals and fractional occupations are treated as variational parameters via a common optimized effective potential (OEP). We show that the OEP exhibits step and peak features which, similar to the exact Kohn-Sham (KS) potential of DFT, are crucial to prevent charge delocalization. We further show that the step is missing in common approximations within reduced density matrix functional theory resulting in delocalization errors comparable to those found in DFT approximations. Finally, we explain the delocalisation error as coming from an artificial mixing of the ground state with a charge-transfer excited state which is avoided if occupation numbers exhibit discontinuities.
Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.
Recently, an approximate theoretical framework was introduced, called local reduced density matrix functional theory (local-RDMFT), where functionals of the one-body reduced density matrix (1-RDM) are minimized under the additional condition that the optimal orbitals satisfy a single electron Schrödinger equation with a local potential. In the present work, we focus on the character of these optimal orbitals. In particular, we compare orbitals obtained by local-RDMFT with those obtained with the full minimization (without the extra condition) by contrasting them against the exact NOs and orbitals from a density functional calculation using the local density approximation (LDA). We find that the orbitals from local-RMDFT are very close to LDA orbitals, contrary to those of the full minimization that resemble the exact NOs. Since local RDMFT preserves the good quality of the description of strong static correlation, this finding opens the way to a mixed density/density matrix scheme, where Kohn-Sham orbitals obtain fractional occupations from a minimization of the occupation numbers using 1-RDM functionals. This will allow for a description of strong correlation at a cost only minimally higher than a density functional calculation.
It is known that the separation of electrons into spinons and chargons, the spin-charge separation, plays a decisive role when describing one-dimensional (1D) strongly correlated systems [ Phys. Rev. B 2012 , 86 , 075132 ]. In this paper, within the density-functional theory (DFT) formalism, we extend the investigation by considering a model for the third electron fractionalization: the separation into spinons, chargons and orbitons, the last associated with the electronic orbital degree of freedom. Specifically, we deal with two exact constraints of exchange-correlation (XC) density-functionals: (i) the constancy of the highest occupied (HO) Kohn-Sham (KS) eigenvalues upon fractional electron numbers and (ii) their discontinuities at integers. By means of 1D discrete Hubbard chains and 1D H2 molecules in the continuum, we find that spin-charge separation yields almost constant HO KS eigenvalues, whereas the spin-orbital counterpart can be decisive when describing derivative discontinuities of XC potentials at strong correlations.
Abstract Density functional theory (DFT) is usually formulated in terms of the electron density as a function of position n(r). Here we discuss an alternative formulation of DFT in terms of the orbital occupation numbers {n α } associated with a local-orbital orthonormal basis set {ϕ α }. First, we discuss how the building blocks of DFT, namely the Hohenberg–Kohn theorems, the Levy–Lieb approach and the Kohn–Sham method, can be adapted for a description in terms of {n α }. In particular, the total energy is now a function of {n α }, E[{n α }], and a Kohn–Sham-like Hamiltonian is derived introducing the effects of the electron–electron interactions via effective potentials, { V α eff = ∂ E e e [ { n β } ] / ∂ n α } . In a second step we consider the Hartree and exchange energies and discuss how to describe them, in the spirit of a DFT approach, in terms of the orbital occupation numbers. In this contribution special attention is paid to the description of the (intra-atomic) correlation energy and corresponding correlation potentials {V corr,α }. For this purpose, a model system is analyzed in detail, whereby an atomic Hamiltonian interacts with the environment via a simplified model; the use of this model allows us to obtain the correlation energy and potentials (in terms of {n α }) for different cases corresponding to low, intermediate and high electron correlations.
Standard approximations for the exchange-correlation functional in Kohn-Sham density functional theory (KS-DFT) typically lead to unacceptably large errors when applied to strongly correlated electronic systems. Partition-DFT (PDFT) is a formally exact reformulation of KS-DFT in which the ground-state density and energy of a system are obtained through self-consistent calculations on isolated fragments, with a partition energy representing inter-fragment interactions. Here, we show how typical errors of the local density approximation (LDA) in KS-DFT can be largely suppressed through a simple approximation, the multi-fragment overlap approximation (MFOA), for the partition energy in PDFT. Our method is illustrated on simple models of one-dimensional strongly correlated linear hydrogen chains. The MFOA, when used in combination with the LDA for the fragments, improves LDA dissociation curves of hydrogen chains and produces results that are comparable to those of spin-unrestricted LDA, but without breaking the spin symmetry. MFOA also induces a correction to the LDA electron density that partially captures the correct density dimerization in strongly correlated hydrogen chains. Moreover, with an additional correction to the partition energy that is specific to the one-dimensional LDA, the approximation is shown to produce dissociation energies in quantitative agreement with calculations based on the density matrix renormalization group method.
This review explains the relationship between density functional theory and strongly correlated models using the simplest possible example, the two-site Hubbard model. The relationship to traditional quantum chemistry is included. Even in this elementary example, where the exact ground-state energy and site occupations can be found analytically, there is much to be explained in terms of the underlying logic and aims of Density Functional Theory. Although the usual solution is analytic, the density functional is given only implicitly. We overcome this difficulty using the Levy-Lieb construction to create a parametrization of the exact function with negligible errors. The symmetric case is most commonly studied, but we find a rich variation in behavior by including asymmetry, as strong correlation physics vies with charge-transfer effects. We explore the behavior of the gap and the many-body Green’s function, demonstrating the ‘failure ’ of the Kohn-Sham method to reproduce the fundamental gap. We perform benchmark calculations of the occupation and components of the KS potentials, the correlation kinetic energies, and the adiabatic connection. We test several approximate functionals (restricted and unrestricted Hartree-Fock and Bethe Ansatz Local Density Approximation) to show their successes and limitations. We also discuss and illustrate the concept of the derivative discontinuity. Useful appendices include analytic expressions for Density