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

Physical review. B, Condensed matter

Hubbard Hamiltonian: Highly correlated electrons from the perspective of the theory of disordered systems

G. Chiappe, Enrique Louis, E. V. Anda

Published Jan 1, 1995• 5 citations
Publisher

Abstract

A diagrammatic technique that treats the kinetic energy as a perturbation on the Coulomb interaction term is used to show the correspondence between the Hubbard Hamiltonian and a disordered alloy of singly and doubly occupied atoms. The method generalizes the alloy-analogy approach. Infinite and finite chains, and finite clusters of the square lattice, are considered. Finite systems are handled by means of a method that includes, in a self-consistent manner, only those configurations with the same number of particles. Total energies and momentum distributions are calculated. The results for chains up to 14 sites and 4\ifmmode\times\else\texttimes\fi{}4 clusters of the square lattice compare satisfactorily with exact results.

DOI
10.1103/physrevb.51.326

Related papers

The atomic method for the Hubbard dimer

Renan Lira Nunes, Marcos Sergio Figueira da Silva

The single-band Hubbard Hamiltonian [1] is the simplest model of interacting electrons in a lattice.It was developed by John Hubbard, who introduced the model to take into account the local electronic correlations in narrow energy bands.The model has a kinetic term that allows hopping of electrons between atomic sites and a term that considers the on-site interaction between electrons.This model is largely used in solid state physics to study magnetic properties of solids, insulator-metal transition (Mott transition) and high temperature superconductors.However, this is a problem that has exact analytical solution in a very few limiting cases and must be treated approximately or numerically.In this work, we propose a new methodology to analytically solve the Hubbard Hamiltonian, mapping it into a two-site model (the Hubbard dimer).To obtain the Green's function for the lattice, we employed the cumulant expansion technique, using as a 'seed' the exact two-site Green's function (the atomic method) [2].From the Green's function, we obtained the density of states [DOS] and the occupation numbers as functions of the external parameters of the model, and compared the results with the Hubbard I approximation [1].

Brazilian Workshop on Semiconductor Physics•Aug 19, 2017•1 citation

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

Anderson-Mott transition in a disordered Hubbard chain with correlated hopping

Francesca Battista, Alberto Camjayi, Liliana Arrachea

We study the ground-state phase diagram of the Anderson-Hubbard model with correlated hopping at half-filling in one dimension. The Hamiltonian has a local Coulomb repulsion $U$ and a disorder potential with local energies randomly distributed in the interval $(\ensuremath{-}W,+W)$ with equal probability, acting on the singly occupied sites. The hopping process which modifies the number of doubly occupied sites is forbidden. The hopping between nearest-neighbor singly occupied and empty sites or between singly occupied and doubly occupied sites has the same amplitude $t$. We identify three different phases as functions of the disorder amplitude $W$ and Coulomb interaction strength $U>0$. When $U<4t$ the system shows a metallic phase: (i) only when no disorder is present $W=0$ or an Anderson-localized phase, (ii) when disorder is introduced $W\ensuremath{\ne}0$. When $U>4t$ the Anderson-localized phase survives as long as disorder effects dominate on the interaction effects, otherwise a Mott-insulator phase (iii) arises. Phases (i) and (ii) are characterized by a finite density of doublons and a vanishing charge gap among the ground state and the excited states. Phase (iii) is characterized by the vanishing density of doublons and a finite gap for the charge excitations.

arXiv:1608.06866•Physical review. B./Physical review. B•Jul 12, 2017•4 citations

Theory of strongly correlated electron systems: Hubbard–Anderson models from an exact Hamiltonian, and perturbation theory near the atomic limit within a nonorthogonal basis set

I. Sandalov, Börje Johansson, Olle Eriksson

Abstract The theory of correlated electron systems is formulated in a form that allows the use as a reference point a density functional theory in the local density approximation (LDA DFT) for solids or molecules. The theory is constructed in two steps. As a first step, the total Hamiltonian is transformed into a correlated form. To elucidate the microscopic origin of the parameters of the periodic Hubbard–Anderson model (PHAM) the terms of the full Hamiltonian that have the operator structure of PHAM are separated. It is found that the matrix element of mixing interaction includes ion configuration‐ and number of particles‐dependent contributions from the Coulomb interaction. In a second step the diagram technique (DT) is developed by means of generalization of the Baym–Kadanoff method for correlated systems. The advantages of the method are that: (1) A nonorthogonal basis can be used, in particular the one generated by LDA DFT; (2) the equations for Green's functions (GFs) for the Fermi and Bose types of quasiparticles can be formulated in the form of a closed system of functional equations. The latter allows us to avoid the question of the nonunique decoupling procedure existing in previous versions of the DT and perform the expansion in terms of dressed GFs. Although the expressions for all interactions depend on the overlap matrix, it is shown that the theory is formally equivalent to one with orthogonal states with redefined interactions. When the PHAM is treated from the atomic‐limit side the vertexes are generated by kinematic interactions. The latter arise due to nontrivial commutation relations between X‐operators and come from the mixing, hopping, and overlap of states. The equations for GFs are derived within the nonorthogonal basis set in Hubbard‐I, one‐loop and random‐phase approximations with respect to kinematic interactions. The self‐consistent equation for “Hubbard Us” is derived. The technique developed is general, in particular its “bosonic” part can be used for description of spin systems with arbitrary anisotropy, systems with orbital ordering, or ordering of multipoles. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 113–143, 2003

arXiv:cond-mat/0011259•International Journal of Quantum Chemistry•Jan 1, 2003•28 citations

Density-matrix functional theory of the Hubbard model: An exact numerical study

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

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 ${N}_{a},$ band filling ${N}_{e},$ 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 ${N}_{a},$ ${N}_{e},$ 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.

arXiv:cond-mat/9910264•Physical review. B, Condensed matter•Jan 15, 2000•45 citations

Correlation between perturbation and variation methods in the study of strongly correlated electron systems

Robinson O. Okanigbuan, JO Idiodi

The ground–state wave function and energy are calculated for two electrons subject to a one-band Hubbard Hamiltonian on a one dimensional lattice containing N electronic sites, N = 2,3,4,5,6, and a 3 x 3 cluster of the square lattice, using perturbation and variational methods. The results from these two approximation methods are then compared with the result from exact calculational method. Journal of the Nigerian Association of Mathematical Physics Vol. 10 2006: pp. 583-594

Journal of the Nigerian Association of Mathematical Physics•May 22, 2008•0 citations
A diagrammatic technique that treats the kinetic energy as a perturbation on the Coulomb interaction term is used to show the correspondence between the Hubbard Hamiltonian and a disordered alloy of singly and doubly occupied atoms. The method generalizes the alloy-analogy approach. Infinite and finite chains, and finite clusters of the square lattice, are considered. Finite systems are handled by means of a method that includes, in a self-consistent manner, only those configurations with the same number of particles. Total energies and momentum distributions are calculated. The results for chains up to 14 sites and 4\ifmmode\times\else\texttimes\fi{}4 clusters of the square lattice compare satisfactorily with exact results.
The single-band Hubbard Hamiltonian [1] is the simplest model of interacting electrons in a lattice.It was developed by John Hubbard, who introduced the model to take into account the local electronic correlations in narrow energy bands.The model has a kinetic term that allows hopping of electrons between atomic sites and a term that considers the on-site interaction between electrons.This model is largely used in solid state physics to study magnetic properties of solids, insulator-metal transition (Mott transition) and high temperature superconductors.However, this is a problem that has exact analytical solution in a very few limiting cases and must be treated approximately or numerically.In this work, we propose a new methodology to analytically solve the Hubbard Hamiltonian, mapping it into a two-site model (the Hubbard dimer).To obtain the Green's function for the lattice, we employed the cumulant expansion technique, using as a 'seed' the exact two-site Green's function (the atomic method) [2].From the Green's function, we obtained the density of states [DOS] and the occupation numbers as functions of the external parameters of the model, and compared the results with the Hubbard I approximation [1].
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.
We study the ground-state phase diagram of the Anderson-Hubbard model with correlated hopping at half-filling in one dimension. The Hamiltonian has a local Coulomb repulsion UUU and a disorder potential with local energies randomly distributed in the interval $(\ensuremath{-}W,+W)$ with equal probability, acting on the singly occupied sites. The hopping process which modifies the number of doubly occupied sites is forbidden. The hopping between nearest-neighbor singly occupied and empty sites or between singly occupied and doubly occupied sites has the same amplitude ttt. We identify three different phases as functions of the disorder amplitude WWW and Coulomb interaction strength U>0U>0U>0. When U<4tU<4tU<4t the system shows a metallic phase: (i) only when no disorder is present W=0W=0W=0 or an Anderson-localized phase, (ii) when disorder is introduced $W\ensuremath{\ne}0$. When U>4tU>4tU>4t the Anderson-localized phase survives as long as disorder effects dominate on the interaction effects, otherwise a Mott-insulator phase (iii) arises. Phases (i) and (ii) are characterized by a finite density of doublons and a vanishing charge gap among the ground state and the excited states. Phase (iii) is characterized by the vanishing density of doublons and a finite gap for the charge excitations.
Abstract The theory of correlated electron systems is formulated in a form that allows the use as a reference point a density functional theory in the local density approximation (LDA DFT) for solids or molecules. The theory is constructed in two steps. As a first step, the total Hamiltonian is transformed into a correlated form. To elucidate the microscopic origin of the parameters of the periodic Hubbard–Anderson model (PHAM) the terms of the full Hamiltonian that have the operator structure of PHAM are separated. It is found that the matrix element of mixing interaction includes ion configuration‐ and number of particles‐dependent contributions from the Coulomb interaction. In a second step the diagram technique (DT) is developed by means of generalization of the Baym–Kadanoff method for correlated systems. The advantages of the method are that: (1) A nonorthogonal basis can be used, in particular the one generated by LDA DFT; (2) the equations for Green's functions (GFs) for the Fermi and Bose types of quasiparticles can be formulated in the form of a closed system of functional equations. The latter allows us to avoid the question of the nonunique decoupling procedure existing in previous versions of the DT and perform the expansion in terms of dressed GFs. Although the expressions for all interactions depend on the overlap matrix, it is shown that the theory is formally equivalent to one with orthogonal states with redefined interactions. When the PHAM is treated from the atomic‐limit side the vertexes are generated by kinematic interactions. The latter arise due to nontrivial commutation relations between X‐operators and come from the mixing, hopping, and overlap of states. The equations for GFs are derived within the nonorthogonal basis set in Hubbard‐I, one‐loop and random‐phase approximations with respect to kinematic interactions. The self‐consistent equation for “Hubbard Us” is derived. The technique developed is general, in particular its “bosonic” part can be used for description of spin systems with arbitrary anisotropy, systems with orbital ordering, or ordering of multipoles. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 94: 113–143, 2003
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,{N}_{a},Na​, band filling Ne,{N}_{e},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,{N}_{a},Na​, Ne,{N}_{e},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)(U>0)(U>0) and attractive (U<0)(U<0)(U<0) interactions are contrasted.