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

Journal of the Physical Society of Japan

Superconducting Correlationof Two-Dimensional Hubbard Model near Half-Filling

Masatoshi Imada

Published Aug 15, 1991• 36 citations
Publisher

Abstract

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.

DOI
10.1143/jpsj.60.2740

Related papers

Detection of pairing correlation in the two-dimensional Hubbard model

Kazuhiko Kuroki, Hideo Aoki, Takashi Hotta +1 more

The quantum Monte Carlo method is used to reexamine superconductivity in the single-band Hubbard model in two dimensions. Instead of the conventional pairing, we consider a ``correlated pairing,'' 〈c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\uparrow}}}$c-${\mathrm{tilde}}_{{\mathrm{i}}^{|\mathrm{I}\mathrm{H}}\mathrm{\ensuremath{\downarrow}}}$〉 with c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$\ensuremath{\equiv}${\mathrm{c}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$(1-${\mathrm{n}}_{\mathrm{i}\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\sigma}}}$), which is inferred from the t-J model, the strong-coupling limit of the Hubbard model. The pairing in the d-wave channel is found to possess both a divergence like 1/T in the pairing susceptibility and a growth of the ground-state pairing correlation with sample size, indicating an off-diagonal long-range order near (but not exactly at) half-filling.

arXiv:cond-mat/9606143•Physical review. B, Condensed matter•Feb 1, 1997•10 citations

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Mingpu Qin, Hao Shi, Shiwei Zhang

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of $U/t$ from 2 to 8, for a range of densities including half-filling and $n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75$, and $0.875$. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.

arXiv:1605.09421•Physical review. B./Physical review. B•Aug 4, 2016•100 citations

Configuration-interaction approach to the two-dimensional Hubbard model near half-filling

Barry Friedman, Gregory Levine

The two-dimensional Hubbard model near half-filling is studied using a configuration-interaction method. Superconducting order-parameter and spin-spin correlation functions are calculated as functions of distance to probe for superconductivity and antiferromagnetism. Consistent with our previous results, these more extensive calculations point toward d-wave superconductivity in the ground state of the two-dimensional repulsive Hubbard model. The d-wave correlations tend to increase with interaction strength, size of the state space used in the calculation, and doping away from half-filling. The superconducting correlations on the Fermi surface are interpreted as d-wave pairing between spin density-wave quasiparticles.

Physical review. B, Condensed matter•Apr 15, 1997•4 citations

Numerical study of the two-dimensional Hubbard model

Steven R. White, Douglas J. Scalapino, Robert L. Sugar +3 more

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.

Physical review. B, Condensed matter•Jul 1, 1989•947 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

Correlations in the two-dimensional Hubbard model

Adriana Moreo, Douglas J. Scalapino

Using Monte Carlo simulations, we have studied the question of phase separation and pairing in the doped two-dimensional Hubbard model. For a range of U/t values greater than as well as less than the bandwidth, we find no evidence for phase separation when the system is doped away from half filling. In agreement with previous work, we find that the d-wave and extended s-wave pairing correlations are enhanced by the interaction vertex but that none of the eight different pairing correlations that we studied increases as a function of lattice size. We have also examined the dependence of the kinetic energy, which determines the total optical spectral weight, on U/t and filling 〈n〉.

Physical review. B, Condensed matter•Apr 1, 1991•54 citations
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.
The quantum Monte Carlo method is used to reexamine superconductivity in the single-band Hubbard model in two dimensions. Instead of the conventional pairing, we consider a ``correlated pairing,'' 〈c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\uparrow}}}$c-${\mathrm{tilde}}_{{\mathrm{i}}^{|\mathrm{I}\mathrm{H}}\mathrm{\ensuremath{\downarrow}}}$〉 with c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$\ensuremath{\equiv}${\mathrm{c}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$(1-${\mathrm{n}}_{\mathrm{i}\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\sigma}}}$), which is inferred from the t-J model, the strong-coupling limit of the Hubbard model. The pairing in the d-wave channel is found to possess both a divergence like 1/T in the pairing susceptibility and a growth of the ground-state pairing correlation with sample size, indicating an off-diagonal long-range order near (but not exactly at) half-filling.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/tU/tU/t from 2 to 8, for a range of densities including half-filling and n=0.3,0.5,0.6,0.75n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75n=0.3,0.5,0.6,0.75, and 0.8750.8750.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
The two-dimensional Hubbard model near half-filling is studied using a configuration-interaction method. Superconducting order-parameter and spin-spin correlation functions are calculated as functions of distance to probe for superconductivity and antiferromagnetism. Consistent with our previous results, these more extensive calculations point toward d-wave superconductivity in the ground state of the two-dimensional repulsive Hubbard model. The d-wave correlations tend to increase with interaction strength, size of the state space used in the calculation, and doping away from half-filling. The superconducting correlations on the Fermi surface are interpreted as d-wave pairing between spin density-wave quasiparticles.
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.
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.
Using Monte Carlo simulations, we have studied the question of phase separation and pairing in the doped two-dimensional Hubbard model. For a range of U/t values greater than as well as less than the bandwidth, we find no evidence for phase separation when the system is doped away from half filling. In agreement with previous work, we find that the d-wave and extended s-wave pairing correlations are enhanced by the interaction vertex but that none of the eight different pairing correlations that we studied increases as a function of lattice size. We have also examined the dependence of the kinetic energy, which determines the total optical spectral weight, on U/t and filling 〈n〉.