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

International Journal of Modern Physics C

NONEXISTENCE OF d-WAVE-SUPERCONDUCTIVITY IN THE QUANTUM MONTE CARLO SIMULATION OF THE HUBBARD MODEL

Hans‐Georg Matuttis, Nobuyasu Ito

Published Jun 1, 2005• 2 citations
Publisher

Abstract

For the existence of d-wave-superconductivity in the Hubbard model, previous quantum Monte Carlo results by other authors, which showed a power law increase of the d-wave susceptibility, seem to contradict a recently published theorem. We show those quantum Monte Carlo calculations were numerically contaminated, analyze the numerical problem and propose a numerically more stable computing scheme.

DOI
10.1142/s0129183105007571

Related papers

Quantum Monte Carlo and the Hubbard Model

Douglas J. Scalapino

Here we discuss quantum Monte Carlo and the Hubbard model. We note that while there is a fermion sign problem for the doped Hubbard model, Monte Carlo calculations have shown that the leading pairfield susceptibility has d-wave symmetry and that the pairing interaction is mediated by S = 1 particle-hole fluctuations.

Journal of Physics Conference Series•Aug 14, 2014•2 citations

Quantum Monte Carlo and the Hubbard model

Douglas J. Scalapino

Here we discuss quantum Monte Carlo and the Hubbard model. We note that while there is a fermion sign problem for the doped Hubbard model, Monte Carlo calculations have shown that the leading pairfield susceptibility has d-wave symmetry and that the pairing interaction is mediated by S = 1 particle-hole fluctuations.

•Jan 1, 1994•0 citations

Competing nodal d-wave superconductivity and antiferromagnetism: a Quantum Monte Carlo study

Xiao Yan Xu, Tarun Grover

Unconventional superconductors such as cuprates often host competing nodal superconductivity and antiferromagnetism. These systems are typically modeled as a repulsive Hubbard model whose unbiased simulation suffers from the fermion sign-problem at any non-zero doping. Here we will exploit the fact that neither of these phases, the nodal d-wave superconductor or the antiferromagnet, require any doping for their existence, and construct a sign-problem-free repulsive Hubbard model with an additional bosonic field which hosts both of these phases. Using Quantum Monte Carlo (QMC) simulations, supplemented with mean-field theory and continuum field-theory arguments, we find that it hosts three distinct phases: a nodal d-wave phase, an antiferromagnet, and an intervening phase which hosts coexisting antiferromagnetism and nodeless d-wave superconductivity. The transition from the coexisting phase to the antiferromagnet is described by the 2+1-D XY universality class, while the one from the coexisting phase to the nodal d-wave phase is described by the chiral Gross-Neveu-Heisenberg theory.

•Sep 16, 2020•3 citations

NUMERICAL STABILITY AND THE SIGN PROBLEM IN THE DETERMINANT QUANTUM MONTE CARLO METHOD

Eugene Loh, James E. Gubernatis, Richard T. Scalettar +3 more

A recent paper by Matuttis and Ito questions the numerical accuracy of a widely-used fermion Monte Carlo algorithm. They also claim that the increase in the d-wave pairfield susceptibility χd(T) of a doped 4×4 Hubbard model at low temperature, previously found using this algorithm, is an artifact due to numerical errors. Here, we provide tests which show that this algorithm is numerically accurate and show that the simulation of χd for a 2×2 lattice agrees with exact diagonalization results. We also provide more complete data for χd on a 4×4 lattice that is consistent with our previous results.

International Journal of Modern Physics C•Aug 1, 2005•31 citations

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

RENORMALISATION GROUP ANALYSIS OF THE TWO-DIMENSIONAL HUBBARD MODEL: SPIN CORRELATIONS AND D-WAVE SUPERCONDUCTIVITY

Walter Metzner

The repulsive two-dimensional Hubbard model is analysed within a functional renormalisation group (RG) approach. The flow equation for the effective interactions is evaluated on 1-loop level. The effective interactions diverge at a finite energy scale which is exponentially small for small bare interactions. To analyse the nature of the instabilities signalled by the diverging interactions, we compute the flow of the superconducting susceptibilities for various pairing symmetries and also charge and spin density susceptibilities. Depending on the choice of the model parameters (hopping amplitudes, interaction strength and band-filling) we find spin density wave instabilities or d-wave superconductivity as leading instability close to half-filling. This calculation establishes the existence of d-wave superconductivity in the Hubbard model beyond doubt.

International Journal of Modern Physics B•Nov 10, 2003•0 citations
For the existence of d-wave-superconductivity in the Hubbard model, previous quantum Monte Carlo results by other authors, which showed a power law increase of the d-wave susceptibility, seem to contradict a recently published theorem. We show those quantum Monte Carlo calculations were numerically contaminated, analyze the numerical problem and propose a numerically more stable computing scheme.
Unconventional superconductors such as cuprates often host competing nodal superconductivity and antiferromagnetism. These systems are typically modeled as a repulsive Hubbard model whose unbiased simulation suffers from the fermion sign-problem at any non-zero doping. Here we will exploit the fact that neither of these phases, the nodal d-wave superconductor or the antiferromagnet, require any doping for their existence, and construct a sign-problem-free repulsive Hubbard model with an additional bosonic field which hosts both of these phases. Using Quantum Monte Carlo (QMC) simulations, supplemented with mean-field theory and continuum field-theory arguments, we find that it hosts three distinct phases: a nodal d-wave phase, an antiferromagnet, and an intervening phase which hosts coexisting antiferromagnetism and nodeless d-wave superconductivity. The transition from the coexisting phase to the antiferromagnet is described by the 2+1-D XY universality class, while the one from the coexisting phase to the nodal d-wave phase is described by the chiral Gross-Neveu-Heisenberg theory.
A recent paper by Matuttis and Ito questions the numerical accuracy of a widely-used fermion Monte Carlo algorithm. They also claim that the increase in the d-wave pairfield susceptibility χd(T) of a doped 4×4 Hubbard model at low temperature, previously found using this algorithm, is an artifact due to numerical errors. Here, we provide tests which show that this algorithm is numerically accurate and show that the simulation of χd for a 2×2 lattice agrees with exact diagonalization results. We also provide more complete data for χd on a 4×4 lattice that is consistent with our previous results.
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.
The repulsive two-dimensional Hubbard model is analysed within a functional renormalisation group (RG) approach. The flow equation for the effective interactions is evaluated on 1-loop level. The effective interactions diverge at a finite energy scale which is exponentially small for small bare interactions. To analyse the nature of the instabilities signalled by the diverging interactions, we compute the flow of the superconducting susceptibilities for various pairing symmetries and also charge and spin density susceptibilities. Depending on the choice of the model parameters (hopping amplitudes, interaction strength and band-filling) we find spin density wave instabilities or d-wave superconductivity as leading instability close to half-filling. This calculation establishes the existence of d-wave superconductivity in the Hubbard model beyond doubt.