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

arXiv:2506.12736 · Journal of the Physical Society of Japan

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

Published Jun 20, 2025• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.7566/jpsj.94.074708

Related papers

Crossover from Weakly to Strongly Correlated Regions in the Two-dimensional Hubbard Model — Off-diagonal Wave Function Monte Carlo Studies of Hubbard Model II —

Takashi Yanagisawa

The ground state of the two-dimensional (2D) Hubbard model is investigated by adopting improved wave functions that take into account intersite electron correlation beyond the Gutzwiller ansatz. The ground-state energy is lowered considerably, giving the best estimate of the ground-state energy for the 2D Hubbard model. There is a crossover from weakly to strongly correlated regions as the on-site Coulomb interaction $U$ increases. The antiferromagnetic correlation induced by $U$ is reduced for hole doping when $U$ is large, being greater than the bandwidth, thus increasing the kinetic energy gain. The spin and charge fluctuations are induced in the strongly correlated region. These antiferromagnetic and kinetic charge fluctuations induce electron pairings, which would result in high-temperature superconductivity.

arXiv:1611.02391•Journal of the Physical Society of Japan•Oct 16, 2016•33 citations

Antiferromagnetism, Superconductivity and Phase Diagram in the Two-Dimensional Hubbard Model —Off-Diagonal Wave Function Monte Carlo Studies of Hubbard Model III—

Takashi Yanagisawa

We investigate the ground-state phase diagram of the two-dimensional Hubbard model based on the optimization variational Monte Carlo method. We use a wave function that is an off-diagonal type given as $ψ=\exp(-λK)P_Gψ_0$, where $ψ_0$ is a one-particle state, $P_G$ is the Gutzwiller operator, $K$ is the kinetic operator, and $λ$ is a variational parameter. The many-body effect plays an important role as an origin of spin correlation and superconductivity in correlated electron systems. We examine the competition between the antiferromagnetic state and superconducting state by varying the Coulomb repulsion $U$, the band parameter $t'$ and the electron density $n_e$. We show a phase diagram that includes superconducting and antiferromagnetic phases and that $t'=0$ is most favorable for superconductivity.

arXiv:1903.09368•Journal of the Physical Society of Japan•Apr 5, 2019•13 citations

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

Takeshi Nakanishi, Yamaji K., Takashi Yanagisawa

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.

Journal of the Physical Society of Japan•Feb 15, 1997•65 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

Correlated-Electron Systems and High-Temperature Superconductivity

Takashi Yanagisawa, Mitake Miyazaki, Yamaji K.

We present recent theoretical results on superconductivity in correlated-electron systems, especially in the two-dimensional Hubbard model and the three-band d-p model. The mechanism of superconductivity in high-temperature superconductors has been extensively studied on the basis of various electronic models and also electron-phonon models. In this study, we investigate the properties of superconductivity in correlated-electron systems by using numerical methods such as the variational Monte Carlo method and the quantum Monte Carlomethod. The Hubbard model is one of basic models for strongly correlated electron systems, and is regarded as the model of cuprate high temperature superconductors. The d-p model is more realistic model for cuprates. The superconducting condensation energy obtained by adopting the Gutzwiller ansatz is in reasonable agreement with the condensation energy estimated for YBa2Cu3O7. We show the phase diagram of the ground state using this method. We have further investigated the stability of striped and checkerboard states in the under-doped region. Holes doped in a half-filled square lattice lead to an incommensurate spin and charge density wave. The relationship of the hole density x and incommensurability δ, δ~x, is satisfied in the lower doping region, as indicated by the variationalMonte Carlocalculations for the two-dimensional Hubbard model. A checkerboard-like charge-density modulation with a roughly period has also been observed by scanning tunneling microscopy experiments in Bi2212 and Na-CCOC compounds. We have performed a variational Monte Carlo simulation on a two-dimensional t-t′-t″- U Hubbard model with a Bi-2212 type band structure and found that the period checkerboard spin modulation, that is characterized by multi Q vectors, is indeed stabilized. We have further performed an investigation by using a quantumMonte Carlomethod, which is a numerical method that can be used to simulate the behavior of correlated electron systems. We present a new algorithm of the quantum Monte Carlo diagonalization that is a method for the evaluation of expectation value without the negative sign problem. We compute pair correlation functions and show that pair correlation is indeed enhanced with hole doping.

arXiv:1307.2456•Journal of Modern Physics•Jan 1, 2013•23 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
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.
The ground state of the two-dimensional (2D) Hubbard model is investigated by adopting improved wave functions that take into account intersite electron correlation beyond the Gutzwiller ansatz. The ground-state energy is lowered considerably, giving the best estimate of the ground-state energy for the 2D Hubbard model. There is a crossover from weakly to strongly correlated regions as the on-site Coulomb interaction UUU increases. The antiferromagnetic correlation induced by UUU is reduced for hole doping when UUU is large, being greater than the bandwidth, thus increasing the kinetic energy gain. The spin and charge fluctuations are induced in the strongly correlated region. These antiferromagnetic and kinetic charge fluctuations induce electron pairings, which would result in high-temperature superconductivity.
We investigate the ground-state phase diagram of the two-dimensional Hubbard model based on the optimization variational Monte Carlo method. We use a wave function that is an off-diagonal type given as ψ=exp⁡(−λK)PGψ0ψ=\exp(-λK)P_Gψ_0ψ=exp(−λK)PG​ψ0​, where ψ0ψ_0ψ0​ is a one-particle state, PGP_GPG​ is the Gutzwiller operator, KKK is the kinetic operator, and λλλ is a variational parameter. The many-body effect plays an important role as an origin of spin correlation and superconductivity in correlated electron systems. We examine the competition between the antiferromagnetic state and superconducting state by varying the Coulomb repulsion UUU, the band parameter t′t't′ and the electron density nen_ene​. We show a phase diagram that includes superconducting and antiferromagnetic phases and that t′=0t'=0t′=0 is most favorable for superconductivity.
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.
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.
We present recent theoretical results on superconductivity in correlated-electron systems, especially in the two-dimensional Hubbard model and the three-band d-p model. The mechanism of superconductivity in high-temperature superconductors has been extensively studied on the basis of various electronic models and also electron-phonon models. In this study, we investigate the properties of superconductivity in correlated-electron systems by using numerical methods such as the variational Monte Carlo method and the quantum Monte Carlomethod. The Hubbard model is one of basic models for strongly correlated electron systems, and is regarded as the model of cuprate high temperature superconductors. The d-p model is more realistic model for cuprates. The superconducting condensation energy obtained by adopting the Gutzwiller ansatz is in reasonable agreement with the condensation energy estimated for YBa2Cu3O7. We show the phase diagram of the ground state using this method. We have further investigated the stability of striped and checkerboard states in the under-doped region. Holes doped in a half-filled square lattice lead to an incommensurate spin and charge density wave. The relationship of the hole density x and incommensurability δ, δ~x, is satisfied in the lower doping region, as indicated by the variationalMonte Carlocalculations for the two-dimensional Hubbard model. A checkerboard-like charge-density modulation with a roughly period has also been observed by scanning tunneling microscopy experiments in Bi2212 and Na-CCOC compounds. We have performed a variational Monte Carlo simulation on a two-dimensional t-t′-t″- U Hubbard model with a Bi-2212 type band structure and found that the period checkerboard spin modulation, that is characterized by multi Q vectors, is indeed stabilized. We have further performed an investigation by using a quantumMonte Carlomethod, which is a numerical method that can be used to simulate the behavior of correlated electron systems. We present a new algorithm of the quantum Monte Carlo diagonalization that is a method for the evaluation of expectation value without the negative sign problem. We compute pair correlation functions and show that pair correlation is indeed enhanced with hole doping.
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.