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:1203.1719 · Journal of the Physical Society of Japan

Variational Monte Carlo Study of Spin-Gapped Normal State and BCS–BEC Crossover in Two-Dimensional Attractive Hubbard Model

Shun Tamura, Hisatoshi Yokoyama

Published Jun 4, 2012• 9 citations
arXivHTMLPDF

Abstract

We study properties of normal, superconducting (SC) and CDW states for an attractive Hubbard model on the square lattice, using a variational Monte Carlo method. In trial wave functions, we introduce an interspinon binding factor, indispensable to induce a spin-gap transition in the normal state, in addition to the onsite attractive and intersite repulsive factors. It is found that, in the normal state, as the interaction strength $|U|/t$ increases, a first-order spin-gap transition arises at $|U_{\rm c}|\sim W$ ($W$: band width) from a Fermi liquid to a spin-gapped state, which is conductive through hopping of doublons. In the SC state, we confirm by analysis of various quantities that the mechanism of superconductivity undergoes a smooth crossover at around $|U_{\ma{co}}|\sim |U_{\rm c}|$ from a BCS type to a Bose-Einstein condensation (BEC) type, as $|U|/t$ increases. For $|U||U_{\ma{co}}|$, quantities such as the energy gain in the SC transition and superfluid stiffness, which is related to the cost of phase coherence, behave as $\sim t^2/|U|\propto T_{\rm c}$, as expected in a bosonic scheme. In this regime, the SC transition is induced by a gain in kinetic energy, in contrast with the BCS theory. We refer to the relevance to the pseudogap in cuprate superconductors.

DOI
10.1143/jpsj.81.064718

Related papers

Spin-gapped wave function as normal state in attractive Hubbard model

Shun Tamura, H. Yokoyama

With the pseudogap in cuprate superconductors in mind, we study an attractive Hubbard model on the square lattice, in which the mechanism of superconductivity undergoes a crossover at |Uco|/t from a BCS type to a Bose-Einstein condensation type, as the interaction strength |U|/t increases. Correspondingly, for |U|>|Uco|, a normatate, realized for T>Tc, is considered to exhibit gapped behavior. We try to construct such a state and study its properties using a variational Monte Carlo method, which treats the local correlation accurately and can connect the weak-and strong-correlation limits continuously. We introduce an intersite correlation factor between antiparallel spins in nearest-neighbor sites, in addition to the onsite (Gutzwiller) projection. This factor brings about a first-order transition from a metal to a spin-gapped phase at |U=|Uc|=W (bandwidth). For |U|>|Uc|, a gap opens in the spin degree of freedom, whereas the charge sector remains gapless. In this phase, conductivity is caused by a hopping of doublons.

Physics Procedia•Jan 1, 2012•0 citations

Crossover between BCS Superconductor and Doped Mott Insulator of d -Wave Pairing State in Two-Dimensional Hubbard Model

Hisatoshi Yokoyama, Masao Ogata, Yukio Tanaka +2 more

With high-Tc cuprates in mind, properties of correlated dx2-y2-wave superconducting (SC) and antiferromagnetic (AF) states are studied for the Hubbard (t-t'-U) model on square lattices, using a variational Monte Carlo method. We employ simple trial wave functions including only crucial parameters, such as a doublon-holon binding factor indispensable to describe correlated SC and normal states as doped Mott insulators. U/t, t'/t and δ(doping rate) dependence of relevant quantities are systematically calculated. As U/t increases, a sharp crossover of SC properties occurs at U_co/t \sim 10 from a conventional BCS type to a kinetic-energy-driven type for any t'/t. As δdecreases, U_co/t is smoothly connected to the Mott transition point at half filling. For U/t\lsim 5, steady superconductivity corresponding to the cuprates is not found, whereas the d-wave SC correlation function Pd^\infty rapidly increases for U/t\gsim 6 and becomes maximum at U=U_co. Comparing the δdependence of Pd^\infty with experimentally observed dome-shaped Tc and condensation energy, we find that the effective value of $U$ for the cuprates should be larger than the band width, for which the t-J model is valid. Analyzing the kinetic energy, we reveal that for U>U_co only doped holes (electrons) become charge carriers, which will make a small Fermi surface (hole pocket), but for UU_co is determined by two factors, i.e., the AF spin correlation, which creates singlet pairs (pseudogap), and the charge mobility dominated by Mott physics.

arXiv:1208.1102•Journal of the Physical Society of Japan•Jan 15, 2013•100 citations

Two-dimensional attractive Hubbard model and the BCS-BEC crossover

Rodrigo A. Fontenele, Natanael C. Costa, Raimundo R. dos Santos +1 more

Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak [Bardeen-Cooper-Schrieffer, (BCS)] and strong [Bose-Einstein condensation (BEC)] couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice from which several different quantities have been calculated and should be useful as a road map to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature ${T}_{c}$ in terms of the band filling $\ensuremath{\langle}n\ensuremath{\rangle}$ and interaction strength $U$ from which we pinpoint a somewhat wide region $|U|/t\ensuremath{\approx}5\ifmmode\pm\else\textpm\fi{}1$ ($t$ is the hopping amplitude) and $\ensuremath{\langle}n\ensuremath{\rangle}\ensuremath{\approx}0.79\ifmmode\pm\else\textpm\fi{}0.09$ leading to a maximum ${T}_{c}\ensuremath{\approx}0.16t$. Two additional temperature scales, namely, pairing ${T}_{p}$ and degeneracy ${T}_{d}$ have been highlighted: The former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), whereas the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.

Physical review. B./Physical review. B•May 2, 2022•34 citations

Variational Monte Carlo Studies of Attractive Hubbard Model. I: Normal-State Properties

H. Yokoyama

Normal states of the attractive Hubbard model, mainly in two dimensions, are studied in the context of a transition from a Fermi liquid to an insulating or gapped state. A series of variational Monte Carlo calculations with better statistics was carried out to estimate accurately the expectation values of several many-body wave functions. Although a relatively clear crossover is observed, even for the plain Gutzwiller wave function, the states in both regimes are metallic. Meanwhile, a substantial metal-insulator transition takes place at ∣U∣∼W (band width) for an improved wave function, into which an intersite correlation is introduced by taking account of virtual states in the second-order perturbation in the infinite ∣U∣ limit. The critical value compares favorably with recent results of the dynamical-mean-field approximation. In contrast, a conventional Jastrow-type wave function scarcely improves the normal state. In addition, the issue of the Brinkman-Rice metal-insulator transition is reconsidered with much larger systems.

arXiv:cond-mat/0205281•Progress of Theoretical Physics•Jul 1, 2002•45 citations

The 2D attractive Hubbard model and the BCS-BEC crossover

Rodrigo A. Fontenele, Natanael C. Costa, Raimundo R. dos Santos +1 more

Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak (Bardeen-Cooper-Schrieffer, BCS) and strong (Bose-Einstein condensation, BEC) couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice, from which several different quantities have been calculated and should be useful as a roadmap to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature, $T_c$, in terms of the band filling, $\ave{n}$, and interaction strength, $U$, from which we pinpoint a somewhat wide region $|U|/t \approx 5 \pm 1$ ($t$ is the hopping amplitude) and $\ave{n} \approx 0.79 \pm 0.09$ leading to a maximum $T_c \approx 0.16 t$. Two additional temperature scales, namely pairing, $T_p$, and degeneracy, $T_d$, have been highlighted: the former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), while the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites, for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.

arXiv:2201.02156•Jan 6, 2022•1 citation

Crossover of Superconducting Properties and Kinetic-Energy Gain in Two-Dimensional Hubbard Model

Hisatoshi Yokoyama, Yukio Tanaka, Masao Ogata +1 more

Superconductivity in the Hubbard model on a square lattice near half filling is studied using an optimization (or correlated) variational Monte Carlo method. Second-order processes of the strong-coupling expansion are considered in the wave functions beyond the Gutzwiller factor. Superconductivity of d_x^2-y^2-wave is widely stable, and exhibits a crossover around U=U_co\sim 12t from a BCS type to a new type. For U\gsim U_co (U\lsim U_co), the energy gain in the superconducting state is derived from the kinetic (potential) energy. Condensation energy is large and \propto exp(-t/J) [tiny] on the strong [weak] coupling side of U_co. Cuprates belong to the strong-coupling regime.

arXiv:cond-mat/0308264•Journal of the Physical Society of Japan•May 15, 2004•87 citations
We study properties of normal, superconducting (SC) and CDW states for an attractive Hubbard model on the square lattice, using a variational Monte Carlo method. In trial wave functions, we introduce an interspinon binding factor, indispensable to induce a spin-gap transition in the normal state, in addition to the onsite attractive and intersite repulsive factors. It is found that, in the normal state, as the interaction strength ∣U∣/t|U|/t∣U∣/t increases, a first-order spin-gap transition arises at ∣Uc∣∼W|U_{\rm c}|\sim W∣Uc​∣∼W (WWW: band width) from a Fermi liquid to a spin-gapped state, which is conductive through hopping of doublons. In the SC state, we confirm by analysis of various quantities that the mechanism of superconductivity undergoes a smooth crossover at around ∣U\maco∣∼∣Uc∣|U_{\ma{co}}|\sim |U_{\rm c}|∣U\maco​∣∼∣Uc​∣ from a BCS type to a Bose-Einstein condensation (BEC) type, as ∣U∣/t|U|/t∣U∣/t increases. For ∣U∣∣U\maco∣|U||U_{\ma{co}}|∣U∣∣U\maco​∣, quantities such as the energy gain in the SC transition and superfluid stiffness, which is related to the cost of phase coherence, behave as ∼t2/∣U∣∝Tc\sim t^2/|U|\propto T_{\rm c}∼t2/∣U∣∝Tc​, as expected in a bosonic scheme. In this regime, the SC transition is induced by a gain in kinetic energy, in contrast with the BCS theory. We refer to the relevance to the pseudogap in cuprate superconductors.
With the pseudogap in cuprate superconductors in mind, we study an attractive Hubbard model on the square lattice, in which the mechanism of superconductivity undergoes a crossover at |Uco|/t from a BCS type to a Bose-Einstein condensation type, as the interaction strength |U|/t increases. Correspondingly, for |U|>|Uco|, a normatate, realized for T>Tc, is considered to exhibit gapped behavior. We try to construct such a state and study its properties using a variational Monte Carlo method, which treats the local correlation accurately and can connect the weak-and strong-correlation limits continuously. We introduce an intersite correlation factor between antiparallel spins in nearest-neighbor sites, in addition to the onsite (Gutzwiller) projection. This factor brings about a first-order transition from a metal to a spin-gapped phase at |U=|Uc|=W (bandwidth). For |U|>|Uc|, a gap opens in the spin degree of freedom, whereas the charge sector remains gapless. In this phase, conductivity is caused by a hopping of doublons.
With high-Tc cuprates in mind, properties of correlated dx2-y2-wave superconducting (SC) and antiferromagnetic (AF) states are studied for the Hubbard (t-t'-U) model on square lattices, using a variational Monte Carlo method. We employ simple trial wave functions including only crucial parameters, such as a doublon-holon binding factor indispensable to describe correlated SC and normal states as doped Mott insulators. U/t, t'/t and δ(doping rate) dependence of relevant quantities are systematically calculated. As U/t increases, a sharp crossover of SC properties occurs at U_co/t \sim 10 from a conventional BCS type to a kinetic-energy-driven type for any t'/t. As δdecreases, U_co/t is smoothly connected to the Mott transition point at half filling. For U/t\lsim 5, steady superconductivity corresponding to the cuprates is not found, whereas the d-wave SC correlation function Pd^\infty rapidly increases for U/t\gsim 6 and becomes maximum at U=U_co. Comparing the δdependence of Pd^\infty with experimentally observed dome-shaped Tc and condensation energy, we find that the effective value of UUU for the cuprates should be larger than the band width, for which the t-J model is valid. Analyzing the kinetic energy, we reveal that for U>U_co only doped holes (electrons) become charge carriers, which will make a small Fermi surface (hole pocket), but for UU_co is determined by two factors, i.e., the AF spin correlation, which creates singlet pairs (pseudogap), and the charge mobility dominated by Mott physics.
Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak [Bardeen-Cooper-Schrieffer, (BCS)] and strong [Bose-Einstein condensation (BEC)] couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice from which several different quantities have been calculated and should be useful as a road map to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature Tc{T}_{c}Tc​ in terms of the band filling $\ensuremath{\langle}n\ensuremath{\rangle}$ and interaction strength UUU from which we pinpoint a somewhat wide region $|U|/t\ensuremath{\approx}5\ifmmode\pm\else\textpm\fi{}1$ (ttt is the hopping amplitude) and $\ensuremath{\langle}n\ensuremath{\rangle}\ensuremath{\approx}0.79\ifmmode\pm\else\textpm\fi{}0.09$ leading to a maximum ${T}_{c}\ensuremath{\approx}0.16t$. Two additional temperature scales, namely, pairing Tp{T}_{p}Tp​ and degeneracy Td{T}_{d}Td​ have been highlighted: The former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), whereas the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.
Normal states of the attractive Hubbard model, mainly in two dimensions, are studied in the context of a transition from a Fermi liquid to an insulating or gapped state. A series of variational Monte Carlo calculations with better statistics was carried out to estimate accurately the expectation values of several many-body wave functions. Although a relatively clear crossover is observed, even for the plain Gutzwiller wave function, the states in both regimes are metallic. Meanwhile, a substantial metal-insulator transition takes place at ∣U∣∼W (band width) for an improved wave function, into which an intersite correlation is introduced by taking account of virtual states in the second-order perturbation in the infinite ∣U∣ limit. The critical value compares favorably with recent results of the dynamical-mean-field approximation. In contrast, a conventional Jastrow-type wave function scarcely improves the normal state. In addition, the issue of the Brinkman-Rice metal-insulator transition is reconsidered with much larger systems.
Recent experiments with ultracold fermionic atoms in optical lattices have provided a tuneable and clean realization of the attractive Hubbard model (AHM). In view of this, several physical properties may be thoroughly studied across the crossover between weak (Bardeen-Cooper-Schrieffer, BCS) and strong (Bose-Einstein condensation, BEC) couplings. Here we report on extensive determinant Quantum Monte Carlo (DQMC) studies of the AHM on a square lattice, from which several different quantities have been calculated and should be useful as a roadmap to experiments. We have obtained a detailed phase diagram for the critical superconducting temperature, TcT_cTc​, in terms of the band filling, $\ave{n}$, and interaction strength, UUU, from which we pinpoint a somewhat wide region ∣U∣/t≈5±1|U|/t \approx 5 \pm 1∣U∣/t≈5±1 (ttt is the hopping amplitude) and $\ave{n} \approx 0.79 \pm 0.09$ leading to a maximum Tc≈0.16tT_c \approx 0.16 tTc​≈0.16t. Two additional temperature scales, namely pairing, TpT_pTp​, and degeneracy, TdT_dTd​, have been highlighted: the former sets the scale for pair formation (believed to be closely related to the scale for the gap of spin excitations in cuprates), while the latter sets the scale for dominant quantum effects. Our DQMC data for the distribution of doubly occupied sites, for the momentum distribution function, and for the quasiparticle weight show distinctive features on both sides of the BCS-BEC crossover, being also suggestive of an underlying crossover between Fermi- and non-Fermi liquid behaviors.
Superconductivity in the Hubbard model on a square lattice near half filling is studied using an optimization (or correlated) variational Monte Carlo method. Second-order processes of the strong-coupling expansion are considered in the wave functions beyond the Gutzwiller factor. Superconductivity of d_x^2-y^2-wave is widely stable, and exhibits a crossover around U=U_co\sim 12t from a BCS type to a new type. For U\gsim U_co (U\lsim U_co), the energy gain in the superconducting state is derived from the kinetic (potential) energy. Condensation energy is large and \propto exp(-t/J) [tiny] on the strong [weak] coupling side of U_co. Cuprates belong to the strong-coupling regime.