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

Numerical study of the two-dimensional Hubbard model

Steven R. White, Douglas J. Scalapino, Robert L. Sugar, Eugene Loh, James E. Gubernatis, Richard T. Scalettar

Published Jul 1, 1989• 947 citations
Publisher

Abstract

We report on a numerical study of the two-dimensional Hubbard model and describe two new algorithms for the simulation of many-electron systems. These algorithms allow one to carry out simulations within the grand canonical ensemble at significantly lower temperatures than had previously been obtained and to calculate ground-state properties with fixed numbers of electrons. We present results for the two-dimensional Hubbard model with half- and quarter-filled bands. Our results support the existence of long-range antiferromagnetic order in the ground state at half-filling and its absence at quarter-filling. Results for the magnetic susceptibility and the momentum occupation along with an upper bound to the spin-wave spectrum are given. We find evidence for an attractive effective d-wave pairing interaction near half-filling but have not found evidence for a phase transition to a superconducting state.

DOI
10.1103/physrevb.40.506

Related papers

Magnetism in the two-dimensional Hubbard model (abstract)

J. E. Hirsch, Hsueh‐Ju Lin

The Hubbard model has been widely used to describe magnetism due to itinerant electrons. However, except for the one-dimensional case, no exact solution of the model exists. We have studied the two-dimensional Hubbard model on a square lattice using a Monte Carlo simulation technique. Simulations were performed on lattices of up to 8×8 in spatial size. We will discuss results for thermodynamic properties and correlation functions as a function of temperature and interaction strength. For the half-filled band case, our results indicate that the system exhibits long-range antiferromagnetic order in the ground state for any nonzero interaction, although substantially reduced from the classical Néel state due to both charge and spin fluctuations. For the non-half-filled band case, our results suggest that no long-range magnetic order exists. We also don’t find any indication of short-ranged ferromagnetic correlations in the parameter range studied. We have also studied the predictions of various approximate treatments for this model, particularly the static approximation and the coherent potential approximation (CPA). By comparison with the Monte Carlo results, we discuss the validity of these schemes. Finally, we discuss the effect of changing the band structure on the magnetic properties of the system.

Journal of Applied Physics•Apr 15, 1985•0 citations

Antiferromagnetism in the Two-Dimensional Hubbard Model

J. E. Hirsch, Simon Y. Tang

Magnetic properties of the two-dimensional Hubbard-model are inferred from results of Monte Carlo simulations. Lattice sizes up to 8\ifmmode\times\else\texttimes\fi{}8 and temperatures down to $T=\frac{t}{20}$ ($t=\mathrm{hopping}$) were studied. The half-filled system is found to exhibit antiferromagnetic long-range order for all values of the Coulomb repulsion $U$. The low-temperature magnetic properties are found to be well described by spin-wave theory with renormalized local moment and spin-wave velocity. Numerical evidence presented suggests that when doped the system loses the long-range order immediately away from half filling.

Physical Review Letters•Jan 30, 1989•225 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

Two-dimensional Hubbard model: Numerical simulation study

J. E. Hirsch

We have studied the two-dimensional Hubbard model on a square lattice with nearest-neighbor hopping. We first discuss the properties of the model within the mean-field approximation: Because of the form of the band structure, some peculiar features are found. We then discuss the simulation algorithm used and compare simulation results with exact results for 6-site chains to test the reliability of the approach. We present results for thermodynamic properties and correlation functions for lattices up to 8\ifmmode\times\else\texttimes\fi{}8 in spatial size. The system is found to be an antiferromagnetic insulator for all values of the coupling constant at zero temperature in the half-filled-band case, but the long-range order is much smaller than predicted by mean-field theory. We perform a finite-size-scaling analysis to determine the character of the transition at zero coupling. For non-half-filled-band cases, our results suggest that the system is always paramagnetic, in contradiction with Hartree-Fock predictions. The system does not show tendency to ferromagnetism nor triplet superconductivity in the parameter range studied. We also discuss some properties of the attractive Hubbard model in the half-filled-band case.

Physical review. B, Condensed matter•Apr 1, 1985•1,082 citations

Superconducting Correlationof Two-Dimensional Hubbard Model near Half-Filling

Masatoshi Imada

Superconducting correlation of the Hubbard model on a square lattice is reinvestigated by quantum Monte Carlo method. Results from more than one order of magnitude larger Monte Carlo steps than before are presented to discuss the convergence of the data to the ground state. By examining the dependence on the choice of initial trial states as well as on the Trotter number in detail, the superconducting correlation in the ground state is reliably obtained up to 12×12 lattice away from the half-filling and up to 16×16 lattice at the half-filling, which substantially improve former results near the half-filling. Various channels of pairing in the ground state show no system size dependence and remain quite short ranged at any filling including the region close to the half-filling. Numerical data show further evidence that neither extended s - nor d -wave superconducting state is possible in the two-dimensional Hubbard model.

Journal of the Physical Society of Japan•Aug 15, 1991•36 citations

Monte Carlo Study of the Two-Dimensional Hubbard Model

J. E. Hirsch

A two-dimensional model of interacting electrons on a lattice, the Hubbard model, has been studied by means of a numerical simulation technique. Results for the half-filled-band sector are presented for lattices of up to 6\ifmmode\times\else\texttimes\fi{}6 spatial sites, as a function of temperature and electron-electron repulsion $U$. In the ground state, the system is found to exhibit antiferromagnetic long-range order for all values of $U$, although greatly reduced from the mean-field-theory predictions.

Physical Review Letters•Nov 14, 1983•77 citations
We report on a numerical study of the two-dimensional Hubbard model and describe two new algorithms for the simulation of many-electron systems. These algorithms allow one to carry out simulations within the grand canonical ensemble at significantly lower temperatures than had previously been obtained and to calculate ground-state properties with fixed numbers of electrons. We present results for the two-dimensional Hubbard model with half- and quarter-filled bands. Our results support the existence of long-range antiferromagnetic order in the ground state at half-filling and its absence at quarter-filling. Results for the magnetic susceptibility and the momentum occupation along with an upper bound to the spin-wave spectrum are given. We find evidence for an attractive effective d-wave pairing interaction near half-filling but have not found evidence for a phase transition to a superconducting state.
The Hubbard model has been widely used to describe magnetism due to itinerant electrons. However, except for the one-dimensional case, no exact solution of the model exists. We have studied the two-dimensional Hubbard model on a square lattice using a Monte Carlo simulation technique. Simulations were performed on lattices of up to 8×8 in spatial size. We will discuss results for thermodynamic properties and correlation functions as a function of temperature and interaction strength. For the half-filled band case, our results indicate that the system exhibits long-range antiferromagnetic order in the ground state for any nonzero interaction, although substantially reduced from the classical Néel state due to both charge and spin fluctuations. For the non-half-filled band case, our results suggest that no long-range magnetic order exists. We also don’t find any indication of short-ranged ferromagnetic correlations in the parameter range studied. We have also studied the predictions of various approximate treatments for this model, particularly the static approximation and the coherent potential approximation (CPA). By comparison with the Monte Carlo results, we discuss the validity of these schemes. Finally, we discuss the effect of changing the band structure on the magnetic properties of the system.
Magnetic properties of the two-dimensional Hubbard-model are inferred from results of Monte Carlo simulations. Lattice sizes up to 8\ifmmode\times\else\texttimes\fi{}8 and temperatures down to T=t20T=\frac{t}{20}T=20t​ (t=hoppingt=\mathrm{hopping}t=hopping) were studied. The half-filled system is found to exhibit antiferromagnetic long-range order for all values of the Coulomb repulsion UUU. The low-temperature magnetic properties are found to be well described by spin-wave theory with renormalized local moment and spin-wave velocity. Numerical evidence presented suggests that when doped the system loses the long-range order immediately away from half filling.
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.
We have studied the two-dimensional Hubbard model on a square lattice with nearest-neighbor hopping. We first discuss the properties of the model within the mean-field approximation: Because of the form of the band structure, some peculiar features are found. We then discuss the simulation algorithm used and compare simulation results with exact results for 6-site chains to test the reliability of the approach. We present results for thermodynamic properties and correlation functions for lattices up to 8\ifmmode\times\else\texttimes\fi{}8 in spatial size. The system is found to be an antiferromagnetic insulator for all values of the coupling constant at zero temperature in the half-filled-band case, but the long-range order is much smaller than predicted by mean-field theory. We perform a finite-size-scaling analysis to determine the character of the transition at zero coupling. For non-half-filled-band cases, our results suggest that the system is always paramagnetic, in contradiction with Hartree-Fock predictions. The system does not show tendency to ferromagnetism nor triplet superconductivity in the parameter range studied. We also discuss some properties of the attractive Hubbard model in the half-filled-band case.
Superconducting correlation of the Hubbard model on a square lattice is reinvestigated by quantum Monte Carlo method. Results from more than one order of magnitude larger Monte Carlo steps than before are presented to discuss the convergence of the data to the ground state. By examining the dependence on the choice of initial trial states as well as on the Trotter number in detail, the superconducting correlation in the ground state is reliably obtained up to 12×12 lattice away from the half-filling and up to 16×16 lattice at the half-filling, which substantially improve former results near the half-filling. Various channels of pairing in the ground state show no system size dependence and remain quite short ranged at any filling including the region close to the half-filling. Numerical data show further evidence that neither extended s - nor d -wave superconducting state is possible in the two-dimensional Hubbard model.
A two-dimensional model of interacting electrons on a lattice, the Hubbard model, has been studied by means of a numerical simulation technique. Results for the half-filled-band sector are presented for lattices of up to 6\ifmmode\times\else\texttimes\fi{}6 spatial sites, as a function of temperature and electron-electron repulsion UUU. In the ground state, the system is found to exhibit antiferromagnetic long-range order for all values of UUU, although greatly reduced from the mean-field-theory predictions.