An interacting lattice model describing the subspace spanned by a set of strongly correlated bands is rigorously coupled to density functional theory to enable ab initio calculations of geometric and topological material properties. The strongly correlated subspace is identified from the occupation number band structure as opposed to a mean-field energy band structure. The self-consistent solution of the many-body model Hamiltonian and a generalized Kohn-Sham equation exactly incorporates momentum-dependent and crystal-symmetric correlations into electronic structure calculations in a way that does not rely on a separation of energy scales. Calculations for a multiorbital Hubbard model demonstrate that the theory accurately reproduces the many-body macroscopic polarization.
An interacting lattice model describing the subspace spanned by a set of strongly correlated bands is rigorously coupled to density functional theory to enable ab initio calculations of geometric and topological material properties. The strongly correlated subspace is identified from the occupation number band structure as opposed to a mean-field energy band structure. The self-consistent solution of the many-body model Hamiltonian and a generalized Kohn-Sham equation exactly incorporates momentum-dependent and crystal-symmetric correlations into electronic structure calculations in a way that does not rely on a separation of energy scales. Calculations for a multiorbital Hubbard model demonstrate that the theory accurately reproduces the many-body macroscopic polarization.
By introducing a set of auxiliary equations representing a many-body system, we have derived an extension of the Kohn-Sham scheme for the density functional theory. These equations consist of a Kohn-Sham-type equation determining single-particle orbitals and an eigen-value equation for an effective many-body problem. A variational method similar to the Kohn-Sham technique was utilized to derive effective interactions as well as effective potentials without artificial substitution of a Hubbard-type interaction and a mean-field correction in the energy functional. The second equation is described by an effective many-body Hamiltonian with both 2-body interactions and mean-field terms. Rigorous formulation of the extended Kohn-Sham equation is also given in accordance with the Hadjisavvas-Theophilou formulation. Our formulation can be interpreted as a way to define models of the strongly correlated electron systems, e.g. the Hubbard model
We propose a new, alternative method for ab-initio calculations of the electronic structure of solids, which has been specifically adapted to treat many-body effects in a more rigorous way than many existing ab-initio methods. We start from a standard band-structure calculation for an effective one-particle Hamiltonian approximately describing the material of interest. This yields a suitable set of one-particle basis functions, from which well localized Wannier functions can be constructed using a method proposed by Marzari and Vanderbilt. Within this (minimal) basis of localized Wannier functions the matrix elements of the non-interacting (one-particle) Hamiltonian as well as the Coulomb matrix elements can be calculated. The result is a many-body Hamiltonian in second quantization with parameters determined from first principles calculations for the material of interest. The Hamiltonian is in the form of a multi-band Hamiltonian in second quantization (a kind of extended, multi-band Hubbard model) such that all the standard many-body methods can be applied. We explicitly show how this approach can be solved in the simplest many-body approximation, the mean-field Hartree-Fock approximation (HFA), which takes into account exact exchange and corrects for self-interaction effects.
A density-functional theory for many-body lattice models is considered in which the single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ is the basic variable. Eigenvalue equations are derived for solving Levy's constrained search of the interaction energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is expressed as the sum of Hartree-Fock energy ${E}_{\mathrm{HF}}[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ and the correlation energy ${E}_{C}[{\ensuremath{\gamma}}_{\mathrm{ij}}].$ Exact results are obtained for ${E}_{C}({\ensuremath{\gamma}}_{12})$ of the Hubbard model on various periodic lattices, where ${\ensuremath{\gamma}}_{\mathrm{ij}}={\ensuremath{\gamma}}_{12}$ for all nearest neighbors i and j. The functional dependence of ${E}_{C}({\ensuremath{\gamma}}_{12})$ is analyzed by varying the number of sites Na, band filling Ne, and lattice structure. The infinite one-dimensional chain and one-, two-, or three-dimensional finite clusters with periodic boundary conditions are considered. The properties of ${E}_{C}({\ensuremath{\gamma}}_{12})$ are discussed in the limits of weak $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{0})$ and strong $({\ensuremath{\gamma}}_{12}\ensuremath{\simeq}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})$ electronic correlations, and in the crossover region $({\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}<~{\ensuremath{\gamma}}_{12}<~{\ensuremath{\gamma}}_{12}^{0}).$ Using an appropriate scaling we observe that ${\ensuremath{\varepsilon}}_{C}{(g}_{12}{)=E}_{C}{/E}_{\mathrm{HF}}$ has a pseudo-universal behavior as a function of ${g}_{12}=({\ensuremath{\gamma}}_{12}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})/({\ensuremath{\gamma}}_{12}^{0}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}).$ The fact that ${\ensuremath{\varepsilon}}_{C}{(g}_{12})$ depends weakly on Na,Ne, and lattice structure suggests that the correlation energy of extended systems could be obtained quite accurately from finite-cluster calculations. Finally, the behaviors of ${E}_{C}({\ensuremath{\gamma}}_{12})$ for repulsive (U>0) and attractive (U<0) interactions are contrasted.
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 (KS) 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 functional energy components, several limits of the exact functional (weak- and strong-coupling, symmetric and asymmetric), various adiabatic connection results, proofs of exact conditions for this model, and the origin of the Hubbard model from a minimal basis model for stretched H2.
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
We present and implement a self-consistent DΓA approach for multi-orbital models and ab initio materials calculations. It is applied to the one-band Hubbard model at various interaction strengths with and without doping, to the two-band Hubbard model with two largely different bandwidths, and to SrVO3. The self-energy feedback reduces critical temperatures compared to dynamical mean-field theory, even to zero temperature in two-dimensions. Compared to a one-shot, non-self-consistent calculation the non-local correlations are significantly reduced when they are strong. In case non-local correlations are weak to moderate as for SrVO3, one-shot calculations are sufficient.