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/0205044 · Physical review. B, Condensed matter

Cluster perturbation theory for Hubbard models

David Sénéchal, Danny Pérez, Dany Plouffe

Published Aug 30, 2002• 238 citations
arXivHTMLPDF

Abstract

Cluster perturbation theory is a technique for calculating the spectral weight of Hubbard models of strongly correlated electrons, which combines exact diagonalizations on small clusters with strong-coupling perturbation theory at leading order. It is exact in both the strong- and weak-coupling limits and provides a good approximation to the spectral function at any wave vector. Following the paper by S\'en\'echal et al. [Phys. Rev. Lett. $84,$ 522 (2000)], we provide a more complete description and derivation of the method. We illustrate some of its capabilities, in particular regarding the effect of doping, the calculation of ground state energy and double occupancy, the disappearance of the Fermi surface in the $t\ensuremath{-}{t}^{\ensuremath{'}}$ Hubbard model, and so on. The method is applicable to any model with on-site repulsion only.

DOI
10.1103/physrevb.66.075129

Related papers

Spectral Weight of the Hubbard Model through Cluster Perturbation Theory

David Sénéchal, Danny Pérez, Michel Pioro-Ladrière

We calculate the spectral weight of the one- and two-dimensional Hubbard models by performing exact diagonalizations of finite clusters and treating intercluster hopping with perturbation theory. Even with relatively modest clusters (e.g., 12 sites), the spectra thus obtained give an accurate description of the exact results. Spin-charge separation (i.e., an extended spectral weight bounded by singularities dispersing with wave vector) is clearly recognized in the one-dimensional Hubbard model, and so is extended spectral weight in the two-dimensional Hubbard model.

arXiv:cond-mat/9908045•Physical Review Letters•Jan 17, 2000•402 citations

Perturbational treatment of correlation effects in the Hubbard model

G. Bulk, Rainer J. Jelitto

In the general frame of a modified many-body perturbation theory, the single-band Hubbard model is treated in the limit of weak electron-electron interaction. In this modification the decomposition of the Hamiltonian into an effective free system and an effective perturbation is rearranged iteratively in each order of the formalism, resulting in a successive reduction of the interaction part. In comparison with conventional perturbation techniques, this means an extended range of validity. Imposing thermodynamic self-consistency, we calculate the wave-vector-dependent electronic self-energy up to second order and analyze the influence of the Coulomb correlations. Several well-known physical properties of the Hubbard model are recovered in characteristic single-particle quantities. Therefore, the effect of the band occupation n and the electron-electron repulsion U on spectral densities, band structures, densities of states, magnetization, etc. is investigated in detail. Some of the most remarkable results are: first, a splitting of the free Bloch band into two quasiparticle subbands when the coupling strength reaches moderate values. These subbands are linked by a small damping-induced background in the density of states and are separated roughly by an energy of the order of U. In addition, the upper subband is shifted to higher energies in the spectrum proportionally to U while the lower one is nearly fixed. Second, we found a paramagnetic-ferromagnetic phase diagram that has a critical coupling strength ${U}_{c}$ and exhibits strong indications of a lower critical occupation number ${n}_{c}$. There is no spontaneous magnetization in the system for parameters smaller than these values. Compared with the Stoner model, one thus obtains striking improvements. In many cases our conclusions for the weak-coupling limit are qualitatively in good agreement with well-known results in the strong-correlation regime.

Physical review. B, Condensed matter•Jan 1, 1990•34 citations

Spectral function of electron-phonon models by cluster perturbation theory

Martin Hohenadler, Markus Aichhorn, Wolfgang von der Linden

Cluster perturbation theory in combination with the Lanczos method is used to compute the one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown that the method allows reliable calculations using relatively small clusters, and at the same time significantly reduces finite-size effects. Results are compared with exact data and the relation to existing work is discussed. We also use a strong-coupling perturbation theory---equivalent to the Hubbard I approximation---to calculate the spectral function of the quarter-filled Holstein model of spinless fermions, starting from the exact atomic-limit Green function. The results agree well with previous calculations within the many-body coherent potential approximation.

arXiv:cond-mat/0306740•Physical review. B, Condensed matter•Nov 19, 2003•56 citations

Applications of Cluster Perturbation Theory Using Quantum Monte Carlo Data

Fei Lin, Erik Schwartz Sørensen, Catherine Kallin +1 more

We study cluster perturbation theory [Phys. Rev. Lett. 84, 522 (2000)] when auxiliary field quantum Monte Carlo method is used for solving the cluster hamiltonian. As a case study, we calculate the spectral functions of the Hubbard model in one and two dimensions and compare our results for the spectral functions to results obtained using exact diagonalization to solve the cluster hamiltonian. The main advantage of using quantum Monte Carlo results as a starting point is that the initial cluster size can be taken to be considerably larger and hence potentially capture more of the relevant physics. The drawback is that quantum Monte Carlo methods yield results at imaginary times with stochastic errors.

arXiv:cond-mat/0512406•International Symposium on High Performance Computing Systems and Applications (HPCS)•Jan 1, 2006•0 citations

Antiferromagnetism in the Hubbard model using a cluster slave-spin method

Wei‐Cheng Lee, Ting-Kuo Lee

The cluster slave-spin method is introduced to systematically investigate the solutions of the Hubbard model including the symmetry-broken phases. In this method, the electron operator is factorized into a fermionic spinon describing the physical spin and a slave-spin describing the charge fluctuations. Following the $U(1)$ formalism derived by Yu and Si [Phys. Rev. B 86, 085104 (2012)], it is shown that the self-consistent equations to explore various symmetry-broken density wave states can be constructed in general with a cluster of multiple slave-spin sites. We employ this method to study the antiferromagnetic (AFM) state in the single band Hubbard model with the two- and four-site clusters of slave spins. While the Hubbard gap, the charge gap due to the doubly occupied states, scales with the Hubbard interaction $U$ as expected, the AFM gap $\mathrm{\ensuremath{\Delta}}$, the gap in the spinon dispersion in the AFM state, exhibits a crossover from the weak- to strong-coupling behaviors as $U$ increases. Our cluster slave-spin method reproduces not only the traditional mean-field behavior of $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}U$ in the weak-coupling limit, but also the behavior of $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}{t}^{2}/U$ predicted by the superexchange mechanism in the strong-coupling limit. In addition, the holon-doublon correlator as functions of $U$ and doping $x$ is also computed, which exhibits a strong tendency toward the holon-doublon binding in the strong coupling regime. We further show that the quasiparticle weight obtained by the cluster slave-spin method is in a good agreement with the generalized Gutzwiller approximation in both AFM and paramagnetic states, and the results can be improved beyond the generalized Gutzwiller approximation as the cluster is enlarged from a single site to four sites. Our results demonstrate that the cluster slave-spin method can be a powerful tool to systematically investigate the strongly correlated system.

arXiv:1707.03042•Physical review. B./Physical review. B•Sep 8, 2017•10 citations

Effect of hole doping on the electronic structure and the Fermi surface in the Hubbard model within norm-conserving cluster pertubation theory

Sergey V. Nikolaev, Sergei G. Ovchinnikov

The concentration dependences of the band structure, spectral weight, density of states, and Fermi surface in the paramagnetic state are studied in the Hubbard model within cluster pertubation theory with 2 × 2 clusters. Representation of the Hubbard X operators makes it possible to control conservation of the spectral weight in constructing cluster perturbation theory. The calculated value of the ground-state energy is in good agreement with the results obtained using nonperturbative methods such as the quantum Monte Carlo method, exact diagonalization of a 4 × 4 cluster, and the variational Monte Carlo method. It is shown that in the case of hole doping, the states in the band gap (in-gap states) lie near the top of the lower Hubbard band for large values of U and near the bottom of the upper band for small U . The concentration dependence of the Fermi surface strongly depends on hopping to second ( t ′) and third ( t ″) neighbors. For parameter values typical of HTSC cuprates, the existence of three concentration regions with different Fermi surfaces is demonstrated. It is shown that broadening of the spectral electron density with an energy resolution typical of contemporary ARPES leads to a pattern of arcs with a length depending on the concentration. Only an order-of-magnitude decrease in the linewidth makes it possible to obtain the true Fermi surface from the spectral density. The kinks associated with strong electron correlations are detected in the dispersion relation below the Fermi level.

Journal of Experimental and Theoretical Physics•Jan 1, 2012•12 citations
Cluster perturbation theory is a technique for calculating the spectral weight of Hubbard models of strongly correlated electrons, which combines exact diagonalizations on small clusters with strong-coupling perturbation theory at leading order. It is exact in both the strong- and weak-coupling limits and provides a good approximation to the spectral function at any wave vector. Following the paper by S'en'echal et al. [Phys. Rev. Lett. 84,84,84, 522 (2000)], we provide a more complete description and derivation of the method. We illustrate some of its capabilities, in particular regarding the effect of doping, the calculation of ground state energy and double occupancy, the disappearance of the Fermi surface in the t\ensuremath−t\ensuremath′t\ensuremath{-}{t}^{\ensuremath{'}}t\ensuremath−t\ensuremath′ Hubbard model, and so on. The method is applicable to any model with on-site repulsion only.
We calculate the spectral weight of the one- and two-dimensional Hubbard models by performing exact diagonalizations of finite clusters and treating intercluster hopping with perturbation theory. Even with relatively modest clusters (e.g., 12 sites), the spectra thus obtained give an accurate description of the exact results. Spin-charge separation (i.e., an extended spectral weight bounded by singularities dispersing with wave vector) is clearly recognized in the one-dimensional Hubbard model, and so is extended spectral weight in the two-dimensional Hubbard model.
In the general frame of a modified many-body perturbation theory, the single-band Hubbard model is treated in the limit of weak electron-electron interaction. In this modification the decomposition of the Hamiltonian into an effective free system and an effective perturbation is rearranged iteratively in each order of the formalism, resulting in a successive reduction of the interaction part. In comparison with conventional perturbation techniques, this means an extended range of validity. Imposing thermodynamic self-consistency, we calculate the wave-vector-dependent electronic self-energy up to second order and analyze the influence of the Coulomb correlations. Several well-known physical properties of the Hubbard model are recovered in characteristic single-particle quantities. Therefore, the effect of the band occupation n and the electron-electron repulsion U on spectral densities, band structures, densities of states, magnetization, etc. is investigated in detail. Some of the most remarkable results are: first, a splitting of the free Bloch band into two quasiparticle subbands when the coupling strength reaches moderate values. These subbands are linked by a small damping-induced background in the density of states and are separated roughly by an energy of the order of U. In addition, the upper subband is shifted to higher energies in the spectrum proportionally to U while the lower one is nearly fixed. Second, we found a paramagnetic-ferromagnetic phase diagram that has a critical coupling strength Uc{U}_{c}Uc​ and exhibits strong indications of a lower critical occupation number nc{n}_{c}nc​. There is no spontaneous magnetization in the system for parameters smaller than these values. Compared with the Stoner model, one thus obtains striking improvements. In many cases our conclusions for the weak-coupling limit are qualitatively in good agreement with well-known results in the strong-correlation regime.
Cluster perturbation theory in combination with the Lanczos method is used to compute the one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown that the method allows reliable calculations using relatively small clusters, and at the same time significantly reduces finite-size effects. Results are compared with exact data and the relation to existing work is discussed. We also use a strong-coupling perturbation theory---equivalent to the Hubbard I approximation---to calculate the spectral function of the quarter-filled Holstein model of spinless fermions, starting from the exact atomic-limit Green function. The results agree well with previous calculations within the many-body coherent potential approximation.
We study cluster perturbation theory [Phys. Rev. Lett. 84, 522 (2000)] when auxiliary field quantum Monte Carlo method is used for solving the cluster hamiltonian. As a case study, we calculate the spectral functions of the Hubbard model in one and two dimensions and compare our results for the spectral functions to results obtained using exact diagonalization to solve the cluster hamiltonian. The main advantage of using quantum Monte Carlo results as a starting point is that the initial cluster size can be taken to be considerably larger and hence potentially capture more of the relevant physics. The drawback is that quantum Monte Carlo methods yield results at imaginary times with stochastic errors.
The cluster slave-spin method is introduced to systematically investigate the solutions of the Hubbard model including the symmetry-broken phases. In this method, the electron operator is factorized into a fermionic spinon describing the physical spin and a slave-spin describing the charge fluctuations. Following the U(1)U(1)U(1) formalism derived by Yu and Si [Phys. Rev. B 86, 085104 (2012)], it is shown that the self-consistent equations to explore various symmetry-broken density wave states can be constructed in general with a cluster of multiple slave-spin sites. We employ this method to study the antiferromagnetic (AFM) state in the single band Hubbard model with the two- and four-site clusters of slave spins. While the Hubbard gap, the charge gap due to the doubly occupied states, scales with the Hubbard interaction UUU as expected, the AFM gap $\mathrm{\ensuremath{\Delta}}$, the gap in the spinon dispersion in the AFM state, exhibits a crossover from the weak- to strong-coupling behaviors as UUU increases. Our cluster slave-spin method reproduces not only the traditional mean-field behavior of $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}U$ in the weak-coupling limit, but also the behavior of $\mathrm{\ensuremath{\Delta}}\ensuremath{\sim}{t}^{2}/U$ predicted by the superexchange mechanism in the strong-coupling limit. In addition, the holon-doublon correlator as functions of UUU and doping xxx is also computed, which exhibits a strong tendency toward the holon-doublon binding in the strong coupling regime. We further show that the quasiparticle weight obtained by the cluster slave-spin method is in a good agreement with the generalized Gutzwiller approximation in both AFM and paramagnetic states, and the results can be improved beyond the generalized Gutzwiller approximation as the cluster is enlarged from a single site to four sites. Our results demonstrate that the cluster slave-spin method can be a powerful tool to systematically investigate the strongly correlated system.
The concentration dependences of the band structure, spectral weight, density of states, and Fermi surface in the paramagnetic state are studied in the Hubbard model within cluster pertubation theory with 2 × 2 clusters. Representation of the Hubbard X operators makes it possible to control conservation of the spectral weight in constructing cluster perturbation theory. The calculated value of the ground-state energy is in good agreement with the results obtained using nonperturbative methods such as the quantum Monte Carlo method, exact diagonalization of a 4 × 4 cluster, and the variational Monte Carlo method. It is shown that in the case of hole doping, the states in the band gap (in-gap states) lie near the top of the lower Hubbard band for large values of U and near the bottom of the upper band for small U . The concentration dependence of the Fermi surface strongly depends on hopping to second ( t ′) and third ( t ″) neighbors. For parameter values typical of HTSC cuprates, the existence of three concentration regions with different Fermi surfaces is demonstrated. It is shown that broadening of the spectral electron density with an energy resolution typical of contemporary ARPES leads to a pattern of arcs with a length depending on the concentration. Only an order-of-magnitude decrease in the linewidth makes it possible to obtain the true Fermi surface from the spectral density. The kinks associated with strong electron correlations are detected in the dispersion relation below the Fermi level.