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

Ground-state properties of the Hubbard model by Lanczos diagonalizations

Federico Becca, Alberto Parola, Sandro Sorella

Published Jun 15, 2000• 27 citations
Publisher

Abstract

We investigate the Hubbard model on a two-dimensional cluster of 18 sites using Lanczos algorithm for different strengths of Coulomb repulsion in the low doping region. Energies and correlation functions are given for half-filling and two holes. In the strong coupling regime we compare our results with the $t\ensuremath{-}J$ model. We confirm the tendency of holes to repel, although a significant value of the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductive order parameter is measured. These data represent a benchmark for future quantum simulations.

DOI
10.1103/physrevb.61.r16287

Related papers

Ground-state properties of the two-dimensional Hubbard model

H. Q. Lin

Using all appropriate symmetries for a two-dimensional square system, we have been able to diagonalize exactly the Hubbard Hamiltonian on the 4\ifmmode\times\else\texttimes\fi{}4 and 18-site lattices. We calculate the ground-state energies, hole-hole interactions, and various correlation functions as functions of the Coulomb repulsion U and band filling, and compare the Hubbard-model results with those obtained from the t-J and strong-coupling models.

Physical review. B, Condensed matter•Oct 1, 1991•28 citations

Static and dynamical properties of doped Hubbard clusters

Elbio R. Dagotto, Adriana Moreo, F. Ortolani +2 more

We study the t-J and the Hubbard models at zero temperature using exact-diagonalization techniques on \ensuremath{\surd}10 \ifmmode\times\else\texttimes\fi{} \ensuremath{\surd}10 and 4\ifmmode\times\else\texttimes\fi{}4 sites clusters. Quantum Monte Carlo simulation results on larger lattices are also presented. All electronic fillings have been analyzed for the three models. We have measured equal-time correlation functions corresponding to various types of order (ranging from ``standard'' staggered spin order to more ``exotic'' possibilities like chiral order), as well as various dynamical properties of these models. Upper bounds for the critical hole doping (${\mathit{x}}_{\mathit{c}}$), where long-range antiferromagnetic order disappears, are presented. It was found that ${\mathit{x}}_{\mathit{c}}$ is very small in agreement with experiments for the high-${\mathit{T}}_{\mathit{c}}$ superconductors. For example, in the t-J model, ${\mathit{x}}_{\mathit{c}}$0.08 at J/t=0.4. However, short-distance spin correlations are important up to much higher dopings producing a sharp well-defined spin-wave-like peak in S(q=(\ensuremath{\pi},\ensuremath{\pi}),\ensuremath{\omega}). Regarding the possibility of phase separation in the Hubbard model, we have studied the behavior of the density of particles, 〈n〉, as a function of the chemical potential, using the Lanczos method on a 4\ifmmode\times\else\texttimes\fi{}4 Hubbard cluster, finding no indications of phase separation for any value of U/t. Then, we conclude that the t-J model at small J/t should not phase separate.In order to compare theoretical predictions with photoemission experiments, we evaluated the electronic density of states N(\ensuremath{\omega}), of the Hubbard and t-J models at several doping fractions. We found that upon doping the antiferromagnetic gap is filled for U\ensuremath{\sim}8t or smaller. The chemical potential moves across the insulating gap as one goes from electron to hole doping of the half-filled cluster, in agreement with x-ray absorption experiments but at variance with photoemission experiments. We have also calculated the optical conductivity, ${\mathrm{\ensuremath{\sigma}}}_{1}$(\ensuremath{\omega}), of the Hubbard and t-J models at all dopings on 4\ifmmode\times\else\texttimes\fi{}4 clusters. Results are compared with experiments and the weight of the Drude peak is presented as a function of couplings and dopings. Spectral weight found at small frequencies is associated with the mid-infrared band observed experimentally in ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$${\mathrm{Sr}}_{\mathit{x}}$${\mathrm{CuO}}_{4}$, and with the states filling the insulating gap in photoemission experiments. An overall good agreement with experiments in the normal state was found. Regarding the possibility of superconductivity n these models, we have studied s-, d-, and p-wave pairing correlations. Naively, the d-wave channel seems enhanced near half filling while the extended s-wave channel seems enhanced from half filling up to 40% doping.However, we found that the enhancement comes from short-distance effects and, thus, no numerical indications of superconductivity were found in these models. We emphasize the importance of analyzing the pairing correlations as a function of distance to distinguish between short- and long-distance effects in the susceptibilities. We also observed that spiral order is enhanced at small J/t and low doping. Uniform chiral rder is suppressed by dynamical holes while staggered chiral order may be enhanced, although with a small plaquette order parameter. We conclude with the observation that the simple one-band Hubbard model with intermediate values of U\ensuremath{\sim}8t may account for many of the ``anomalous'' properties of the normal state of the high-${\mathit{T}}_{\mathit{c}}$ superconductors.

Physical review. B, Condensed matter•May 1, 1992•201 citations

HOLE-HOLE CORRELATIONS IN THE 2D HUBBARD MODEL

Sandro Sorella, Alberto Parola, Erio Tosatti

The hole-hole correlation function in the 2D Hubbard model is investigated in the low doping regime both via Monte Carlo simulations for U=4 and Lanczos diagonalization on a 16 site cluster for very large U. The results are compared with available theoretical predictions for the strong coupling limit. There is a clear short range repulsion, but no signature of hole-hole attraction is found in either case. At U=4 correlations are still remarkably weak-coupling-like. For large U, they become rather similar to those of hard core bosons.

International Journal of Modern Physics B•Jan 1, 1991•3 citations

Fidelity study of the superconducting phase diagram in the two-dimensional single-band Hubbard model

Chunjing Jia, Brian Moritz, Cheng-Chien Chen +2 more

Extensive numerical studies have demonstrated that the two-dimensional single-band Hubbard model contains much of the key physics in cuprate high-temperature superconductors. However, there is no definitive proof that the Hubbard model truly possesses a superconducting ground state or, if it does, of how it depends on model parameters. To answer these longstanding questions, we study an extension of the Hubbard model including an infinite-range $d$-wave pair field term, which precipitates a superconducting state in the $d$-wave channel. Using exact diagonalization on 16-site square clusters, we study the evolution of the ground state as a function of the strength of the pairing term. This is achieved by monitoring the fidelity metric of the ground state, as well as determining the ratio between the two largest eigenvalues of the $d$-wave pair/spin/charge-density matrices. The calculations show a $d$-wave superconducting ground state in doped clusters bracketed by a strong antiferromagnetic state at half filling controlled by the Coulomb repulsion $U$ and a weak short-range checkerboard charge ordered state at larger hole doping controlled by the next-nearest-neighbor hopping ${t}^{\ensuremath{'}}$. We also demonstrate that negative ${t}^{\ensuremath{'}}$ plays an important role in facilitating $d$-wave superconductivity.

arXiv:1012.4013•Physical Review B•Sep 9, 2011•24 citations

Pairing in the two-dimensional Hubbard model: An exact diagonalization study

H. Q. Lin, J. E. Hirsch, Douglas J. Scalapino

We have studied the pair susceptibilities for all possible pair wave functions that fit on a two-dimensional (2D) eight-site Hubbard cluster by exact diagonalization of the Hamiltonian. Band fillings corresponding to four and six electrons were studied (two or four holes in the half-filled band) for a wide range of Hubbard interaction strengths and temperatures. Our results show that all pairing susceptibilities are suppressed by the Hubbard repulsion. We have also carried out perturbation-theory calculations which show that the leading-order ${U}^{2}$ contributions to the d-wave pair susceptibility suppresses d-wave pairing over a significant temperature range. These results are consistent with recent Monte Carlo results and provide further evidence suggesting that the 2D Hubbard model does not exhibit superconductivity.

Physical review. B, Condensed matter•May 1, 1988•80 citations

Superconductivity and antiferromagnetism in the two-dimensional Hubbard model: A variational study

David Eichenberger, D. Baeriswyl

A variational ground state of the repulsive Hubbard model on a square lattice is investigated numerically for an intermediate coupling strength $(U=8t)$ and for moderate sizes (from $6\ifmmode\times\else\texttimes\fi{}6$ to $10\ifmmode\times\else\texttimes\fi{}10$). Our ansatz is superior to other widely used variational wave functions. The results for order parameters and correlation functions provide insight into the antiferromagnetic state at half-filling as well as strong evidence for a superconducting phase away from half-filling.

arXiv:0708.2795•Physical Review B•Nov 30, 2007•71 citations
We investigate the Hubbard model on a two-dimensional cluster of 18 sites using Lanczos algorithm for different strengths of Coulomb repulsion in the low doping region. Energies and correlation functions are given for half-filling and two holes. In the strong coupling regime we compare our results with the t\ensuremath−Jt\ensuremath{-}Jt\ensuremath−J model. We confirm the tendency of holes to repel, although a significant value of the dx2\ensuremath−y2{d}_{{x}^{2}\ensuremath{-}{y}^{2}}dx2\ensuremath−y2​-wave superconductive order parameter is measured. These data represent a benchmark for future quantum simulations.
We study the t-J and the Hubbard models at zero temperature using exact-diagonalization techniques on \ensuremath{\surd}10 \ifmmode\times\else\texttimes\fi{} \ensuremath{\surd}10 and 4\ifmmode\times\else\texttimes\fi{}4 sites clusters. Quantum Monte Carlo simulation results on larger lattices are also presented. All electronic fillings have been analyzed for the three models. We have measured equal-time correlation functions corresponding to various types of order (ranging from ``standard'' staggered spin order to more ``exotic'' possibilities like chiral order), as well as various dynamical properties of these models. Upper bounds for the critical hole doping (xc{\mathit{x}}_{\mathit{c}}xc​), where long-range antiferromagnetic order disappears, are presented. It was found that xc{\mathit{x}}_{\mathit{c}}xc​ is very small in agreement with experiments for the high-Tc{\mathit{T}}_{\mathit{c}}Tc​ superconductors. For example, in the t-J model, ${\mathit{x}}_{\mathit{c}}$0.08 at J/t=0.4. However, short-distance spin correlations are important up to much higher dopings producing a sharp well-defined spin-wave-like peak in S(q=(\ensuremath{\pi},\ensuremath{\pi}),\ensuremath{\omega}). Regarding the possibility of phase separation in the Hubbard model, we have studied the behavior of the density of particles, 〈n〉, as a function of the chemical potential, using the Lanczos method on a 4\ifmmode\times\else\texttimes\fi{}4 Hubbard cluster, finding no indications of phase separation for any value of U/t. Then, we conclude that the t-J model at small J/t should not phase separate.In order to compare theoretical predictions with photoemission experiments, we evaluated the electronic density of states N(\ensuremath{\omega}), of the Hubbard and t-J models at several doping fractions. We found that upon doping the antiferromagnetic gap is filled for U\ensuremath{\sim}8t or smaller. The chemical potential moves across the insulating gap as one goes from electron to hole doping of the half-filled cluster, in agreement with x-ray absorption experiments but at variance with photoemission experiments. We have also calculated the optical conductivity, ${\mathrm{\ensuremath{\sigma}}}_{1}$(\ensuremath{\omega}), of the Hubbard and t-J models at all dopings on 4\ifmmode\times\else\texttimes\fi{}4 clusters. Results are compared with experiments and the weight of the Drude peak is presented as a function of couplings and dopings. Spectral weight found at small frequencies is associated with the mid-infrared band observed experimentally in ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$Srx{\mathrm{Sr}}_{\mathit{x}}Srx​CuO4{\mathrm{CuO}}_{4}CuO4​, and with the states filling the insulating gap in photoemission experiments. An overall good agreement with experiments in the normal state was found. Regarding the possibility of superconductivity n these models, we have studied s-, d-, and p-wave pairing correlations. Naively, the d-wave channel seems enhanced near half filling while the extended s-wave channel seems enhanced from half filling up to 40% doping.However, we found that the enhancement comes from short-distance effects and, thus, no numerical indications of superconductivity were found in these models. We emphasize the importance of analyzing the pairing correlations as a function of distance to distinguish between short- and long-distance effects in the susceptibilities. We also observed that spiral order is enhanced at small J/t and low doping. Uniform chiral rder is suppressed by dynamical holes while staggered chiral order may be enhanced, although with a small plaquette order parameter. We conclude with the observation that the simple one-band Hubbard model with intermediate values of U\ensuremath{\sim}8t may account for many of the ``anomalous'' properties of the normal state of the high-Tc{\mathit{T}}_{\mathit{c}}Tc​ superconductors.
The hole-hole correlation function in the 2D Hubbard model is investigated in the low doping regime both via Monte Carlo simulations for U=4 and Lanczos diagonalization on a 16 site cluster for very large U. The results are compared with available theoretical predictions for the strong coupling limit. There is a clear short range repulsion, but no signature of hole-hole attraction is found in either case. At U=4 correlations are still remarkably weak-coupling-like. For large U, they become rather similar to those of hard core bosons.
Extensive numerical studies have demonstrated that the two-dimensional single-band Hubbard model contains much of the key physics in cuprate high-temperature superconductors. However, there is no definitive proof that the Hubbard model truly possesses a superconducting ground state or, if it does, of how it depends on model parameters. To answer these longstanding questions, we study an extension of the Hubbard model including an infinite-range ddd-wave pair field term, which precipitates a superconducting state in the ddd-wave channel. Using exact diagonalization on 16-site square clusters, we study the evolution of the ground state as a function of the strength of the pairing term. This is achieved by monitoring the fidelity metric of the ground state, as well as determining the ratio between the two largest eigenvalues of the ddd-wave pair/spin/charge-density matrices. The calculations show a ddd-wave superconducting ground state in doped clusters bracketed by a strong antiferromagnetic state at half filling controlled by the Coulomb repulsion UUU and a weak short-range checkerboard charge ordered state at larger hole doping controlled by the next-nearest-neighbor hopping ${t}^{\ensuremath{'}}$. We also demonstrate that negative ${t}^{\ensuremath{'}}$ plays an important role in facilitating ddd-wave superconductivity.
We have studied the pair susceptibilities for all possible pair wave functions that fit on a two-dimensional (2D) eight-site Hubbard cluster by exact diagonalization of the Hamiltonian. Band fillings corresponding to four and six electrons were studied (two or four holes in the half-filled band) for a wide range of Hubbard interaction strengths and temperatures. Our results show that all pairing susceptibilities are suppressed by the Hubbard repulsion. We have also carried out perturbation-theory calculations which show that the leading-order U2{U}^{2}U2 contributions to the d-wave pair susceptibility suppresses d-wave pairing over a significant temperature range. These results are consistent with recent Monte Carlo results and provide further evidence suggesting that the 2D Hubbard model does not exhibit superconductivity.
A variational ground state of the repulsive Hubbard model on a square lattice is investigated numerically for an intermediate coupling strength (U=8t)(U=8t)(U=8t) and for moderate sizes (from $6\ifmmode\times\else\texttimes\fi{}6$ to $10\ifmmode\times\else\texttimes\fi{}10$). Our ansatz is superior to other widely used variational wave functions. The results for order parameters and correlation functions provide insight into the antiferromagnetic state at half-filling as well as strong evidence for a superconducting phase away from half-filling.