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

Variational Monte Carlo Indications ofd-Wave Superconductivity in the Two-Dimensional Hubbard Model

Takeshi Nakanishi, Yamaji K., Takashi Yanagisawa

Published Feb 15, 1997• 65 citations
Publisher

Abstract

The possibility of superconductivity in the ground state of the two-dimensional Hubbard model was investigated by means of the variational Monte Carlo method. With slight hole doping on the 6 ×6 and 10 ×10 lattices, a clear minimum was obtained on the curve of the total energy versus the amplitude of the d -wave gap function for large on-site Coulomb energies. Superconducting pair correlations were observed to be well developed. The-next-nearest-neighbor transfer was found to appreciably affect the minimum depth.

DOI
10.1143/jpsj.66.294

Related papers

Superconductivity and Inhomogeneous Charge-Ordered State in the Two-Dimensional Hubbard Model —Off-Diagonal Wave Function Monte Carlo Studies of Hubbard Model IV—

Takashi Yanagisawa

We investigate the ground-state properties of the two-dimensional Hubbard model, based on the off-diagonal wave function variational Monte Carlo method. We use an optimized wave function that is improved from an initial one-body wave function by multiplying by multiple correlation operators that are given by the form of exp($-S$)-type where $S$ is a suitable operator. We examine the inhomogeneous ground state at near 1/8 doping in the strongly correlated region where the on-site Coulomb interaction is larger than the bandwidth. The $d$-wave superconductivity with a spatially oscillating gap function can coexist with the charge ordering in the ground state without magnetic ordering, and superconducting condensation energy increases by this coexistence. We also show the $d$-wave pair correlation function as a function of lattice sites. The correlation function indicates that the long-range superconducting order indeed exists, and also that the strong electron correlation suppresses the pair correlation function.

arXiv:2506.12736•Journal of the Physical Society of Japan•Jun 20, 2025•0 citations

Superconductivity in the 2D Hubbard model: Electron doping is different

David Eichenberger, D. Baeriswyl

A variational Monte Carlo calculation is used for studying the ground state of the two-dimensional Hubbard model, including hopping between both nearest and next-nearest neighbor sites. Superconductivity with d-wave symmetry is found to be restricted to densities where the Fermi surface crosses the magnetic zone boundary. The condensate energy is much larger for hole doping than for electron doping. Superconductivity is kinetic energy driven for hole doping, but potential energy driven for electron doping. Our findings agree surprisingly well with experimental data for layered cuprates, both for electron- and hole-doped materials.

arXiv:0808.0433•Aug 4, 2008•1 citation

Enhanced pair correlation functions in the two-dimensional Hubbard model

Takashi Yanagisawa

In this study we have computed the pair correlation functions in the two-dimensional Hubbard model using a quantum Monte Carlo method. We employ a new diagonalization algorithm in the quantum Monte Carlo method which is free from the negative sign problem. We show that the d-wave pairing correlation function is indeed enhanced slightly for the positive on-site Coulomb interaction U when doping away from the half-filling. When the system size becomes large, the pair correlation function P d increases for U > 0 compared to the non-interacting case, while P d is suppressed for U > 0 when the system size is small. The enhancement ratio P d [ U ]/ P d [ U = 0] will give a criterion on the existence of superconductivity. The ratio P d [ U ]/ P d [ U = 0] increases almost linearly ∝ L when the system size L × L is increased. This increase is a good indication of the existence of a superconducting phase in the two-dimensional Hubbard model. There is, however, no enhancement of pair correlation functions in the half-filled case, which indicates the absence of superconductivity without hole doping.

arXiv:1304.1273•New Journal of Physics•Mar 11, 2013•26 citations

Variational ground states of the two-dimensional Hubbard model

D. Baeriswyl, David Eichenberger, M. Z. Menteshashvili

Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.

arXiv:0907.1593•New Journal of Physics•Jul 23, 2009•45 citations

Superconducting states of an extended Hubbard model

S. N. Coppersmith

The variational Monte Carlo method is used to study possible superconducting states of a two-dimensional extended Hubbard model on a copper-oxygen lattice. It is found that d-wave and extended s-wave superconducting variational wave functions can have lower energy than the corresponding normal state for different parameter regions, while simple s-wave superconductivity is never favored. Increasing the nearest-neighbor repulsion and decreasing the copper-oxygen level separation can promote extended s-wave superconductivity and inhibit d-wave superconductivity.

Physical review. B, Condensed matter•Aug 1, 1990•15 citations

Stripes and d-wave Superconductivity in the Two-dimensional Hubbard Model

Mitake Miyazaki, Yamaji K., Takashi Yanagisawa

In order to study the material dependence of the correlation between stripes and dx2-y2-wave superconductivity (SC) among cuprate superconductors, we have performed a variational Monte Carlo simulation on a two-dimensional Hubbard model with first-, second- and third-nearest-neighbor hopping terms (t, t’ and t”, respectively) by using the Gutzwiller-Jastrow-type wave function taking account of three-body correlations on a triangle, which is efficient for frustrated systems. The condensation energies of the vertical stripes and the homogeneous SC were computed as a function of t’/t with a fixed relation t”= –t’/2. We have found that, at the nearly optimal doping with moderately strong on-site Coulomb interaction, the state of maximum condensation energy switches from the vertical stripes to the homogeneous dSC as |t’/t| increases.

Physics Procedia•Jan 1, 2014•2 citations
The possibility of superconductivity in the ground state of the two-dimensional Hubbard model was investigated by means of the variational Monte Carlo method. With slight hole doping on the 6 ×6 and 10 ×10 lattices, a clear minimum was obtained on the curve of the total energy versus the amplitude of the d -wave gap function for large on-site Coulomb energies. Superconducting pair correlations were observed to be well developed. The-next-nearest-neighbor transfer was found to appreciably affect the minimum depth.
We investigate the ground-state properties of the two-dimensional Hubbard model, based on the off-diagonal wave function variational Monte Carlo method. We use an optimized wave function that is improved from an initial one-body wave function by multiplying by multiple correlation operators that are given by the form of exp(−S-S−S)-type where SSS is a suitable operator. We examine the inhomogeneous ground state at near 1/8 doping in the strongly correlated region where the on-site Coulomb interaction is larger than the bandwidth. The ddd-wave superconductivity with a spatially oscillating gap function can coexist with the charge ordering in the ground state without magnetic ordering, and superconducting condensation energy increases by this coexistence. We also show the ddd-wave pair correlation function as a function of lattice sites. The correlation function indicates that the long-range superconducting order indeed exists, and also that the strong electron correlation suppresses the pair correlation function.
A variational Monte Carlo calculation is used for studying the ground state of the two-dimensional Hubbard model, including hopping between both nearest and next-nearest neighbor sites. Superconductivity with d-wave symmetry is found to be restricted to densities where the Fermi surface crosses the magnetic zone boundary. The condensate energy is much larger for hole doping than for electron doping. Superconductivity is kinetic energy driven for hole doping, but potential energy driven for electron doping. Our findings agree surprisingly well with experimental data for layered cuprates, both for electron- and hole-doped materials.
In this study we have computed the pair correlation functions in the two-dimensional Hubbard model using a quantum Monte Carlo method. We employ a new diagonalization algorithm in the quantum Monte Carlo method which is free from the negative sign problem. We show that the d-wave pairing correlation function is indeed enhanced slightly for the positive on-site Coulomb interaction U when doping away from the half-filling. When the system size becomes large, the pair correlation function P d increases for U > 0 compared to the non-interacting case, while P d is suppressed for U > 0 when the system size is small. The enhancement ratio P d [ U ]/ P d [ U = 0] will give a criterion on the existence of superconductivity. The ratio P d [ U ]/ P d [ U = 0] increases almost linearly ∝ L when the system size L × L is increased. This increase is a good indication of the existence of a superconducting phase in the two-dimensional Hubbard model. There is, however, no enhancement of pair correlation functions in the half-filled case, which indicates the absence of superconductivity without hole doping.
Recent refinements of analytical and numerical methods have improved our understanding of the ground-state phase diagram of the two-dimensional (2D) Hubbard model. Here, we focus on variational approaches, but comparisons with both quantum cluster and Gaussian Monte Carlo methods are also made. Our own ansatz leads to an antiferromagnetic ground state at half filling with a slightly reduced staggered order parameter (as compared to simple mean-field theory). Away from half filling, we find d-wave superconductivity, but confined to densities where the Fermi surface passes through the antiferromagnetic zone boundary (if hopping between both nearest-neighbour and next-nearest-neighbour sites is considered). Our results agree surprisingly well with recent numerical studies using the quantum cluster method. An interesting trend is found by comparing gap parameters Δ (antiferromagnetic or superconducting) obtained with different variational wave functions. Δ varies by an order of magnitude and thus cannot be taken as a characteristic energy scale. In contrast, the order parameter is much less sensitive to the degree of sophistication of the variational schemes, at least at and near half filling.
The variational Monte Carlo method is used to study possible superconducting states of a two-dimensional extended Hubbard model on a copper-oxygen lattice. It is found that d-wave and extended s-wave superconducting variational wave functions can have lower energy than the corresponding normal state for different parameter regions, while simple s-wave superconductivity is never favored. Increasing the nearest-neighbor repulsion and decreasing the copper-oxygen level separation can promote extended s-wave superconductivity and inhibit d-wave superconductivity.
In order to study the material dependence of the correlation between stripes and dx2-y2-wave superconductivity (SC) among cuprate superconductors, we have performed a variational Monte Carlo simulation on a two-dimensional Hubbard model with first-, second- and third-nearest-neighbor hopping terms (t, t’ and t”, respectively) by using the Gutzwiller-Jastrow-type wave function taking account of three-body correlations on a triangle, which is efficient for frustrated systems. The condensation energies of the vertical stripes and the homogeneous SC were computed as a function of t’/t with a fixed relation t”= –t’/2. We have found that, at the nearly optimal doping with moderately strong on-site Coulomb interaction, the state of maximum condensation energy switches from the vertical stripes to the homogeneous dSC as |t’/t| increases.