Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0311470 · Physical Review B

Interaction-energy functional for lattice density functional theory: Applications to one-, two-, and three-dimensional Hubbard models

Román López-Sandoval, G. M. Pastor

Published Feb 13, 2004• 29 citations
arXivHTMLPDF

Abstract

The Hubbard model is investigated in the framework of lattice density functional theory (LDFT). The single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable of the many-body problem. A new approximation to the interaction-energy functional $W[\ensuremath{\gamma}]$ is proposed which is based on its scaling properties and which recovers exactly the limit of strong electron correlations at half-band filling. In this way, a more accurate description of W is obtained throughout the domain of representability of ${\ensuremath{\gamma}}_{\mathrm{ij}},$ including the crossover from weak to strong correlations. As examples of applications results are given for the ground-state energy, charge-excitation gap, and charge susceptibility of the Hubbard model in one-, two-, and three-dimensional lattices. The performance of the method is demonstrated by comparison with available exact solutions, with numerical calculations, and with LDFT using a simpler dimer ansatz for W. Goals and limitations of the different approximations are discussed.

DOI
10.1103/physrevb.69.085101

Related papers

Density-matrix functional study of the Hubbard model on one- and two-dimensional bipartite lattices

Matthieu Saubanère, G. M. Pastor

The inhomogeneous Hubbard model is investigated in one and two dimensions by using lattice density-functional theory (LDFT). The ground-state energy $E=T+W$ is regarded as a functional of the single-particle density matrix $\ensuremath{\gamma}$, where $T[\ensuremath{\gamma}]$ represents the kinetic and crystal-field energy, and $W[\ensuremath{\gamma}]$ the interaction energy. Besides the known functional $T[\ensuremath{\gamma}]$ we propose a simple scaling approximation to the interaction energy $W[\ensuremath{\gamma}]$, which is based on exact results for the Hubbard dimer and on a scaling hypothesis within the domain of representability of $\ensuremath{\gamma}$. As applications we consider the Hubbard model on one- and two-dimensional bipartite lattices. Several ground-state properties are determined including the kinetic and Coulomb energy, density distribution, nearest-neighbor bond order, and charge gap. Comparison with exact Lanczos diagonalizations shows that LDFT with the scaled dimer approximation to $W[\ensuremath{\gamma}]$ yields a quite accurate description of the interplay between correlations and charge redistributions, from weak to strong coupling regimes, and for all band fillings. Goals and limitations of the present approach are discussed.

Physical Review B•Jul 20, 2011•29 citations

Density-matrix functional theory of strongly correlated lattice fermions

Román López-Sandoval, G. M. Pastor

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.

arXiv:cond-mat/0207429•Physical review. B, Condensed matter•Oct 29, 2002•43 citations

Electronic properties of the dimerized one-dimensional Hubbard model using lattice density-functional theory

Román López-Sandoval, G. M. Pastor

The dimerized one-dimensional Hubbard model is studied in the framework of lattice density-functional theory (LDFT). The single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable. The corresponding interaction-energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is defined by Levy's constrained search. Exact numerical results are obtained for $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23}),$ where ${\ensuremath{\gamma}}_{12}={\ensuremath{\gamma}}_{i,i+1}$ for odd i and ${\ensuremath{\gamma}}_{23}={\ensuremath{\gamma}}_{i,i+1}$ for even i are the nearest-neighbor density-matrix elements along the chain. The domain of representability of ${\ensuremath{\gamma}}_{\mathrm{ij}}$ and the functional dependence of $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ are analyzed. A simple, explicit approximation to $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ is proposed, which is derived from scaling properties of W, exact dimer results, and known limits. Using this approximation, LDFT is applied to determine ground-state properties and charge-excitation gaps of finite and infinite dimerized chains as a function of the Coulomb-repulsion strength $U/t$ and of the alternation $\ensuremath{\delta}t$ of the hopping integrals ${t}_{\mathrm{ij}}$ ${(t}_{\mathrm{ij}}=t\ifmmode\pm\else\textpm\fi{}\ensuremath{\delta}t).$ The accuracy of the method is demonstrated by comparison with available exact solutions and accurate numerical calculations. Goals and limitations of the present approach are discussed particularly concerning its ability to describe the crossover from weak to strong electron correlations.

arXiv:cond-mat/0210658•Physical review. B, Condensed matter•Jan 24, 2003•24 citations

Development and applications of density matrix functional theory of generalized Hubbard Models

Matthieu Saubanère

In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.

•Jan 25, 2012•0 citations

Lattice density-functional theory of the attractive Hubbard model

Matthieu Saubanère, G. M. Pastor

The attractive Hubbard model is investigated in the framework of lattice density-functional theory (LDFT). The ground-state energy $E=T+W$ is regarded as a functional of the single-particle density matrix ${\ensuremath{\gamma}}_{ij}$ with respect to the lattice sites, where $T[\ensuremath{\gamma}]$ represents the kinetic and crystal-field energies and $W[\ensuremath{\gamma}]$ the interaction energy. Aside from the exactly known functional $T[\ensuremath{\gamma}]$, we propose a simple scaling approximation to $W[\ensuremath{\gamma}]$, which is based on exact analytic results for the attractive Hubbard dimer and on a scaling hypothesis within the domain of representability of $\ensuremath{\gamma}$. As applications, we consider one-, two-, and three-dimensional finite and extended bipartite lattices having homogeneous or alternating onsite energy levels. In addition, the Bethe lattice is investigated as a function of coordination number. Results are given for the kinetic, Coulomb, and total energies, as well as for the density distribution ${\ensuremath{\gamma}}_{ii}$, nearest-neighbor bond order ${\ensuremath{\gamma}}_{ij}$, and pairing energy $\ensuremath{\Delta}{E}_{p}$, as a function of the interaction strength $|U|/t$, onsite potential $\ensuremath{\varepsilon}/t$, and band filling $n={N}_{e}/{N}_{a}$. Remarkable even-odd and super-even oscillations of $\ensuremath{\Delta}{E}_{p}$ are observed in finite rings as a function of band filling. Comparison with exact Lanczos diagonalizations and density-matrix renormalization-group calculations shows that LDFT yields a very good quantitative description of the properties of the model in the complete parameter range, thus providing a significant improvement over the mean-field approaches. Goals and limitations of the method are discussed.

Physical Review B•Sep 17, 2014•12 citations

Scaling and transferability of the interaction-energy functional of the inhomogeneous Hubbard model

Matthieu Saubanère, G. M. Pastor

The inhomogeneous Hubbard model is investigated in the framework of lattice density-functional theory (LDFT) by considering the single-particle density matrix ${\ensuremath{\gamma}}_{ij}$ with respect to the lattice sites as the basic variable of the many-body problem. The domain of representability of ${\ensuremath{\gamma}}_{ij}$ is determined for charge-density wave states on finite bipartite lattices. Levy's constrained search of the interaction-energy functional $W[{\ensuremath{\gamma}}_{ij}]$ is numerically solved by applying the Lanczos method to an effective Hubbard-type model. The exact functional dependence of $W[{\ensuremath{\gamma}}_{ij}]$ is analyzed by varying systematically the charge transfer $\ensuremath{\Delta}n={\ensuremath{\gamma}}_{22}\ensuremath{-}{\ensuremath{\gamma}}_{11}$, the degree of electron delocalization ${g}_{12}$ between the sublattices, the number of sites ${N}_{a}$, and the band filling $n=({\ensuremath{\gamma}}_{11}+{\ensuremath{\gamma}}_{22})/2={N}_{e}/{N}_{a}$. For each $\ensuremath{\Delta}n$ the properties of $W$ 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, as well as in the crossover region $({\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}\ensuremath{\le}{\ensuremath{\gamma}}_{12}\ensuremath{\le}{\ensuremath{\gamma}}_{12}^{0})$. It is shown that $W$ follows quite closely a simple scaling behavior as a function of $\ensuremath{\Delta}n$ and ${g}_{12}=({\ensuremath{\gamma}}_{12}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})/({\ensuremath{\gamma}}_{12}^{0}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})$. The very good transferability of $W(\ensuremath{\Delta}n,{g}_{12})$ for different ${N}_{a}$, $n$ and lattice structure opens new possibilities of applying LDFT to inhomogeneous many-body models.

Physical Review B•Jun 1, 2009•17 citations
The Hubbard model is investigated in the framework of lattice density functional theory (LDFT). The single-particle density matrix \ensuremathγij{\ensuremath{\gamma}}_{\mathrm{ij}}\ensuremathγij​ with respect to the lattice sites is considered as the basic variable of the many-body problem. A new approximation to the interaction-energy functional W[\ensuremathγ]W[\ensuremath{\gamma}]W[\ensuremathγ] is proposed which is based on its scaling properties and which recovers exactly the limit of strong electron correlations at half-band filling. In this way, a more accurate description of W is obtained throughout the domain of representability of \ensuremathγij,{\ensuremath{\gamma}}_{\mathrm{ij}},\ensuremathγij​, including the crossover from weak to strong correlations. As examples of applications results are given for the ground-state energy, charge-excitation gap, and charge susceptibility of the Hubbard model in one-, two-, and three-dimensional lattices. The performance of the method is demonstrated by comparison with available exact solutions, with numerical calculations, and with LDFT using a simpler dimer ansatz for W. Goals and limitations of the different approximations are discussed.
The inhomogeneous Hubbard model is investigated in one and two dimensions by using lattice density-functional theory (LDFT). The ground-state energy E=T+WE=T+WE=T+W is regarded as a functional of the single-particle density matrix $\ensuremath{\gamma}$, where $T[\ensuremath{\gamma}]$ represents the kinetic and crystal-field energy, and $W[\ensuremath{\gamma}]$ the interaction energy. Besides the known functional $T[\ensuremath{\gamma}]$ we propose a simple scaling approximation to the interaction energy $W[\ensuremath{\gamma}]$, which is based on exact results for the Hubbard dimer and on a scaling hypothesis within the domain of representability of $\ensuremath{\gamma}$. As applications we consider the Hubbard model on one- and two-dimensional bipartite lattices. Several ground-state properties are determined including the kinetic and Coulomb energy, density distribution, nearest-neighbor bond order, and charge gap. Comparison with exact Lanczos diagonalizations shows that LDFT with the scaled dimer approximation to $W[\ensuremath{\gamma}]$ yields a quite accurate description of the interplay between correlations and charge redistributions, from weak to strong coupling regimes, and for all band fillings. Goals and limitations of the present approach are discussed.
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.
The dimerized one-dimensional Hubbard model is studied in the framework of lattice density-functional theory (LDFT). The single-particle density matrix ${\ensuremath{\gamma}}_{\mathrm{ij}}$ with respect to the lattice sites is considered as the basic variable. The corresponding interaction-energy functional $W[{\ensuremath{\gamma}}_{\mathrm{ij}}]$ is defined by Levy's constrained search. Exact numerical results are obtained for $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23}),$ where ${\ensuremath{\gamma}}_{12}={\ensuremath{\gamma}}_{i,i+1}$ for odd i and ${\ensuremath{\gamma}}_{23}={\ensuremath{\gamma}}_{i,i+1}$ for even i are the nearest-neighbor density-matrix elements along the chain. The domain of representability of ${\ensuremath{\gamma}}_{\mathrm{ij}}$ and the functional dependence of $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ are analyzed. A simple, explicit approximation to $W({\ensuremath{\gamma}}_{12},{\ensuremath{\gamma}}_{23})$ is proposed, which is derived from scaling properties of W, exact dimer results, and known limits. Using this approximation, LDFT is applied to determine ground-state properties and charge-excitation gaps of finite and infinite dimerized chains as a function of the Coulomb-repulsion strength U/tU/tU/t and of the alternation $\ensuremath{\delta}t$ of the hopping integrals tij{t}_{\mathrm{ij}}tij​ ${(t}_{\mathrm{ij}}=t\ifmmode\pm\else\textpm\fi{}\ensuremath{\delta}t).$ The accuracy of the method is demonstrated by comparison with available exact solutions and accurate numerical calculations. Goals and limitations of the present approach are discussed particularly concerning its ability to describe the crossover from weak to strong electron correlations.
In this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
The attractive Hubbard model is investigated in the framework of lattice density-functional theory (LDFT). The ground-state energy E=T+WE=T+WE=T+W is regarded as a functional of the single-particle density matrix ${\ensuremath{\gamma}}_{ij}$ with respect to the lattice sites, where $T[\ensuremath{\gamma}]$ represents the kinetic and crystal-field energies and $W[\ensuremath{\gamma}]$ the interaction energy. Aside from the exactly known functional $T[\ensuremath{\gamma}]$, we propose a simple scaling approximation to $W[\ensuremath{\gamma}]$, which is based on exact analytic results for the attractive Hubbard dimer and on a scaling hypothesis within the domain of representability of $\ensuremath{\gamma}$. As applications, we consider one-, two-, and three-dimensional finite and extended bipartite lattices having homogeneous or alternating onsite energy levels. In addition, the Bethe lattice is investigated as a function of coordination number. Results are given for the kinetic, Coulomb, and total energies, as well as for the density distribution ${\ensuremath{\gamma}}_{ii}$, nearest-neighbor bond order ${\ensuremath{\gamma}}_{ij}$, and pairing energy $\ensuremath{\Delta}{E}_{p}$, as a function of the interaction strength ∣U∣/t|U|/t∣U∣/t, onsite potential $\ensuremath{\varepsilon}/t$, and band filling n=Ne/Nan={N}_{e}/{N}_{a}n=Ne​/Na​. Remarkable even-odd and super-even oscillations of $\ensuremath{\Delta}{E}_{p}$ are observed in finite rings as a function of band filling. Comparison with exact Lanczos diagonalizations and density-matrix renormalization-group calculations shows that LDFT yields a very good quantitative description of the properties of the model in the complete parameter range, thus providing a significant improvement over the mean-field approaches. Goals and limitations of the method are discussed.
The inhomogeneous Hubbard model is investigated in the framework of lattice density-functional theory (LDFT) by considering the single-particle density matrix ${\ensuremath{\gamma}}_{ij}$ with respect to the lattice sites as the basic variable of the many-body problem. The domain of representability of ${\ensuremath{\gamma}}_{ij}$ is determined for charge-density wave states on finite bipartite lattices. Levy's constrained search of the interaction-energy functional $W[{\ensuremath{\gamma}}_{ij}]$ is numerically solved by applying the Lanczos method to an effective Hubbard-type model. The exact functional dependence of $W[{\ensuremath{\gamma}}_{ij}]$ is analyzed by varying systematically the charge transfer $\ensuremath{\Delta}n={\ensuremath{\gamma}}_{22}\ensuremath{-}{\ensuremath{\gamma}}_{11}$, the degree of electron delocalization g12{g}_{12}g12​ between the sublattices, the number of sites Na{N}_{a}Na​, and the band filling $n=({\ensuremath{\gamma}}_{11}+{\ensuremath{\gamma}}_{22})/2={N}_{e}/{N}_{a}$. For each $\ensuremath{\Delta}n$ the properties of WWW 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, as well as in the crossover region $({\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}}\ensuremath{\le}{\ensuremath{\gamma}}_{12}\ensuremath{\le}{\ensuremath{\gamma}}_{12}^{0})$. It is shown that WWW follows quite closely a simple scaling behavior as a function of $\ensuremath{\Delta}n$ and ${g}_{12}=({\ensuremath{\gamma}}_{12}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})/({\ensuremath{\gamma}}_{12}^{0}\ensuremath{-}{\ensuremath{\gamma}}_{12}^{\ensuremath{\infty}})$. The very good transferability of $W(\ensuremath{\Delta}n,{g}_{12})$ for different Na{N}_{a}Na​, nnn and lattice structure opens new possibilities of applying LDFT to inhomogeneous many-body models.