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:1012.4013 · Physical Review B

Fidelity study of the superconducting phase diagram in the two-dimensional single-band Hubbard model

Chunjing Jia, Brian Moritz, Cheng-Chien Chen, B. Sriram Shastry, Thomas Devereaux

Published Sep 9, 2011• 24 citations
arXivHTMLPDF

Abstract

Extensive numerical studies have demonstrated that the two-dimensional single-band Hubbard model contains much of the key physics in cuprate high-temperature superconductors. However, there is no definitive proof that the Hubbard model truly possesses a superconducting ground state or, if it does, of how it depends on model parameters. To answer these longstanding questions, we study an extension of the Hubbard model including an infinite-range $d$-wave pair field term, which precipitates a superconducting state in the $d$-wave channel. Using exact diagonalization on 16-site square clusters, we study the evolution of the ground state as a function of the strength of the pairing term. This is achieved by monitoring the fidelity metric of the ground state, as well as determining the ratio between the two largest eigenvalues of the $d$-wave pair/spin/charge-density matrices. The calculations show a $d$-wave superconducting ground state in doped clusters bracketed by a strong antiferromagnetic state at half filling controlled by the Coulomb repulsion $U$ and a weak short-range checkerboard charge ordered state at larger hole doping controlled by the next-nearest-neighbor hopping ${t}^{\ensuremath{'}}$. We also demonstrate that negative ${t}^{\ensuremath{'}}$ plays an important role in facilitating $d$-wave superconductivity.

DOI
10.1103/physrevb.84.125113

Related papers

Superconducting phases of the square-lattice extended Hubbard model

Wei-Chih Chen, Yao Wang, Cheng-Chien Chen

We study the square-lattice extended Hubbard model with on-site $U$ and nearest-neighbor $V$ interactions by exact diagonalization. We show that nonequilibrium quench dynamics can help determine the equilibrium phase transition boundaries, which agree with the calculations of the fidelity metric, dynamical structure factor, and correlation function. At half filling, the phase diagrams in the strong-coupling regime include spin density wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity at large positive $U$, charge density wave (extended ${s}^{*}$-wave superconductivity) at large positive (negative) $V$, and $s$-wave superconductivity at large negative $U$ with vanishing $V$. The energies of different particle sectors also help determine the phase separation region. With carrier doping, charge fluctuation results in a strong competition between different orders, making it more difficult to identify the leading instability on finite-size clusters. Nevertheless, the more exotic $p$-wave superconducting pairing is found to be enhanced when the system is heavily overdoped by $37.5%--50%$ holes, especially in the interaction parameter range relevant to the cuprate superconductors.

arXiv:2206.01119•Physical review. B./Physical review. B•Aug 18, 2023•23 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

Evolution of the superconductivity dome in the two-dimensional Hubbard model

K.-S. Chen, Zi Yang Meng, Shuxiang Yang +3 more

In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the $d$-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature ${T}_{c}$, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the $d$-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.

arXiv:1308.5946•Physical Review B•Dec 10, 2013•33 citations

Physics of cuprates with the two-band Hubbard model: The validity of the one-band Hubbard model

Alexandru Macridin, Mark Jarrell, Thomas Maier +1 more

We calculate the properties of the two-band Hubbard model using the dynamical cluster approximation. The phase diagram resembles the generic phase diagram of the cuprates, showing a strong asymmetry with respect to electron- and hole-doped regimes, in agreement with experiment. Asymmetric features are also seen in one-particle spectral functions and in the charge, spin, and $d$-wave pairing-susceptibility functions. We address the possible reduction of the two-band model to a low-energy single-band one, as it was suggested by Zhang and Rice. Comparing the two-band Hubbard model properties with the single-band Hubbard model ones, we have found similar low-energy physics provided that the next-nearest-neighbor hopping term ${t}^{\ensuremath{'}}$ has a significant value $({t}^{\ensuremath{'}}∕t\ensuremath{\approx}0.3)$. The parameter ${t}^{\ensuremath{'}}$ is the main culprit for the electron-hole asymmetry. However, a significant value of ${t}^{\ensuremath{'}}$ cannot be provided in a strict Zhang and Rice [Phys. Rev. B 37, R3759 (1988); 41, 7243 (1990)] picture where the extra holes added into the system bind to the existing Cu holes forming local singlets. We note that by considering approximate singlet states, such as plaquette states, reasonable values of ${t}^{\ensuremath{'}}$, which capture qualitatively the physics of the two-band model, can be obtained. We conclude that a single-band $t\text{\ensuremath{-}}{t}^{\ensuremath{'}}\text{\ensuremath{-}}U$ Hubbard model captures the basic physics of the cuprates concerning superconductivity, antiferromagnetism, pseudogap, and electron-hole asymmetry, but is not suitable for a quantitative analysis or to describe physical properties involving energy scales larger than about 0.5 eV.

arXiv:cond-mat/0411092•Physical Review B•Apr 29, 2005•95 citations

Global phase diagram of D-wave superconductivity in the square-lattice t-J model

Feng Chen, F. D. M. Haldane, D. N. Sheng

The Hubbard and closely related t-J models are exciting platforms for unconventional superconductivity (SC). Through state-of-the-art density matrix renormalization group calculations using the grand canonical ensemble, we address open issues regarding the ground-state phase diagram of the extended t-J model on a square lattice in the parameter regime relevant to cuprate superconductors. On large 8-leg cylinders, we demonstrate that the pure t-J model with only nearest-neighbor hoppings and superexchange interactions, for a wide range of doping ( δ = 0.1 to 0.2 ), hosts robust d-wave superconductivity possibly coexisting with weak unidirectional pair density wave. Furthermore, a small next nearest neighbor hopping t 2 suppresses pair and charge density waves, resulting in a uniform d-wave SC phase in both electron- and hole-doped cuprate model systems. Our work validates the t-J model as a proper minimum model for the emergence of superconductivity in cuprate superconductors.

arXiv:2311.15092•Proceedings of the National Academy of Sciences•Mar 21, 2025•27 citations

Stripes and spin-density waves in the doped two-dimensional Hubbard model: Ground state phase diagram

Hao Xu, Hao Shi, Ettore Vitali +2 more

We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near $1/8$ doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about $1/5$ doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.

arXiv:2112.02187•Physical Review Research•Mar 28, 2022•74 citations
Extensive numerical studies have demonstrated that the two-dimensional single-band Hubbard model contains much of the key physics in cuprate high-temperature superconductors. However, there is no definitive proof that the Hubbard model truly possesses a superconducting ground state or, if it does, of how it depends on model parameters. To answer these longstanding questions, we study an extension of the Hubbard model including an infinite-range ddd-wave pair field term, which precipitates a superconducting state in the ddd-wave channel. Using exact diagonalization on 16-site square clusters, we study the evolution of the ground state as a function of the strength of the pairing term. This is achieved by monitoring the fidelity metric of the ground state, as well as determining the ratio between the two largest eigenvalues of the ddd-wave pair/spin/charge-density matrices. The calculations show a ddd-wave superconducting ground state in doped clusters bracketed by a strong antiferromagnetic state at half filling controlled by the Coulomb repulsion UUU and a weak short-range checkerboard charge ordered state at larger hole doping controlled by the next-nearest-neighbor hopping t\ensuremath′{t}^{\ensuremath{'}}t\ensuremath′. We also demonstrate that negative t\ensuremath′{t}^{\ensuremath{'}}t\ensuremath′ plays an important role in facilitating ddd-wave superconductivity.
We study the square-lattice extended Hubbard model with on-site UUU and nearest-neighbor VVV interactions by exact diagonalization. We show that nonequilibrium quench dynamics can help determine the equilibrium phase transition boundaries, which agree with the calculations of the fidelity metric, dynamical structure factor, and correlation function. At half filling, the phase diagrams in the strong-coupling regime include spin density wave and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave superconductivity at large positive UUU, charge density wave (extended s∗{s}^{*}s∗-wave superconductivity) at large positive (negative) VVV, and sss-wave superconductivity at large negative UUU with vanishing VVV. The energies of different particle sectors also help determine the phase separation region. With carrier doping, charge fluctuation results in a strong competition between different orders, making it more difficult to identify the leading instability on finite-size clusters. Nevertheless, the more exotic ppp-wave superconducting pairing is found to be enhanced when the system is heavily overdoped by 37.537.5%--50%37.5 holes, especially in the interaction parameter range relevant to the cuprate superconductors.
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.
In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the ddd-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature Tc{T}_{c}Tc​, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the ddd-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.
We calculate the properties of the two-band Hubbard model using the dynamical cluster approximation. The phase diagram resembles the generic phase diagram of the cuprates, showing a strong asymmetry with respect to electron- and hole-doped regimes, in agreement with experiment. Asymmetric features are also seen in one-particle spectral functions and in the charge, spin, and ddd-wave pairing-susceptibility functions. We address the possible reduction of the two-band model to a low-energy single-band one, as it was suggested by Zhang and Rice. Comparing the two-band Hubbard model properties with the single-band Hubbard model ones, we have found similar low-energy physics provided that the next-nearest-neighbor hopping term ${t}^{\ensuremath{'}}$ has a significant value $({t}^{\ensuremath{'}}∕t\ensuremath{\approx}0.3)$. The parameter ${t}^{\ensuremath{'}}$ is the main culprit for the electron-hole asymmetry. However, a significant value of ${t}^{\ensuremath{'}}$ cannot be provided in a strict Zhang and Rice [Phys. Rev. B 37, R3759 (1988); 41, 7243 (1990)] picture where the extra holes added into the system bind to the existing Cu holes forming local singlets. We note that by considering approximate singlet states, such as plaquette states, reasonable values of ${t}^{\ensuremath{'}}$, which capture qualitatively the physics of the two-band model, can be obtained. We conclude that a single-band $t\text{\ensuremath{-}}{t}^{\ensuremath{'}}\text{\ensuremath{-}}U$ Hubbard model captures the basic physics of the cuprates concerning superconductivity, antiferromagnetism, pseudogap, and electron-hole asymmetry, but is not suitable for a quantitative analysis or to describe physical properties involving energy scales larger than about 0.5 eV.
The Hubbard and closely related t-J models are exciting platforms for unconventional superconductivity (SC). Through state-of-the-art density matrix renormalization group calculations using the grand canonical ensemble, we address open issues regarding the ground-state phase diagram of the extended t-J model on a square lattice in the parameter regime relevant to cuprate superconductors. On large 8-leg cylinders, we demonstrate that the pure t-J model with only nearest-neighbor hoppings and superexchange interactions, for a wide range of doping ( δ = 0.1 to 0.2 ), hosts robust d-wave superconductivity possibly coexisting with weak unidirectional pair density wave. Furthermore, a small next nearest neighbor hopping t 2 suppresses pair and charge density waves, resulting in a uniform d-wave SC phase in both electron- and hole-doped cuprate model systems. Our work validates the t-J model as a proper minimum model for the emergence of superconductivity in cuprate superconductors.
We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near 1/81/81/8 doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about 1/51/51/5 doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.