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:cond-mat/0101352 · Physical review. B, Condensed matter

Monte Carlo study of the superfluid weight in doped antiferromagnets

G. C. Psaltakis

Published Mar 23, 2001• 0 citations
arXivHTMLPDF

Abstract

The phase fluctuations of the condensate in doped antiferromagnets, described by a ${t\ensuremath{-}t}^{\ensuremath{'}}\ensuremath{-}J$ model and a suitable $1/N$ expansion, provide a mechanism for a Kosterlitz-Thouless (KT) type of transition to a superconducting state below ${T}_{c}.$ In this paper, we present a Monte Carlo study of the corresponding superfluid weight ${D}_{s}(T)$ in the classical (large-$N)$ limit, as a function of temperature and doping. Consistent with generic experimental trends, ${D}_{s}(T)$ exhibits a T-linear decrease at low temperatures, with the magnitude of the slope ${D}_{s}^{\ensuremath{'}}(0)$ increasing upon doping. Finite-size scaling in the underdoped regime predicts values for the dimensionless ratio ${A=k}_{B}{T}_{c}{/D}_{s}(0)$ of order unity, with $A=0.4435(5)$ in the half filled band limit, thus confirming ${D}_{s}(0)$ as the fundamental energy scale determining ${T}_{c}.$ Our Monte Carlo results for ${D}_{s}{(T)/D}_{s}(0)$ vs ${k}_{B}{T/D}_{s}(0),$ at 10% hole doping, are found to be in reasonable agreement with recent measurements on ${\mathrm{La}}_{2\ensuremath{-}x}{\mathrm{Sr}}_{x}{\mathrm{CuO}}_{4},$ with $x=0.10,$ throughout the temperature range below the theoretical KT transition temperature ${T}_{c}.$

DOI
10.1103/physrevb.63.144529

Related papers

Doping and temperature dependence of superfluid weight and spectral function: Universal scaling behavior of the pseudogap

Sung-Sik Lee, Sung‐Ho Suck Salk

Using the improved U(1) and SU(2) slave-boson approaches of the $t\ensuremath{-}J$ Hamiltonian [S.-S. Lee and Sung-Ho Suck Salk, Phys. Rev. B 64, 052501 (2001)] that we developed recently, we study the doping and temperature dependence of superfluid weight and spectral function and discuss our finding of the universal scaling behavior of the pseudogap. It is shown that at low hole doping concentrations $x$ and at low temperatures $T$, there exists a propensity of linear decrease in the superfluid weight ${n}_{s}∕{m}^{*}$ with temperature $T$, and a tendency of doping independence in the linearly decreasing slopes of $({n}_{s}∕{m}^{*})(x,T)$ with $T$ in qualitative agreement with the experimentally observed relation; that is, $({n}_{s}∕{m}^{*})(x,T)=({n}_{s}∕{m}^{*})(x,0)\ensuremath{-}\ensuremath{\alpha}T$, where $\ensuremath{\alpha}$ is a constant. It is also demonstrated that there exists a boomerang behavior; that is, both ${T}_{c}$ and ${n}_{s}∕{m}^{*}$ increase with hole doping concentration $x$ in the underdoped region, reach a saturation (maximum) at a hole doping above optimal doping, and decrease beyond the saturation point in the overdoped region, in agreement with $\ensuremath{\mu}\mathrm{SR}$ measurements. Based on our improved SU(2) slave-boson approach to the $t\ensuremath{-}J$ Hamiltonian, we further investigate the doping and temperature dependence of the spectral function and discuss our theoretical finding of a scaling behavior of the pseudogap. In addition, we discuss the cause of the hump and quasiparticle peak in the observed spectral functions of high ${T}_{c}$ cuprates. It is demonstrated that the sharpening of the observed quasiparticle peak below ${T}_{c}$ is attributed to the bose condensation of holon pair, in agreement with observations. From the predicted ratios of the pseudogap ${\ensuremath{\Delta}}_{0}$ to both the superconducting temperature ${T}_{c}$ and the pseudogap temperature ${T}^{*}$, i.e., $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}_{c}$ and $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}^{*}$, respectively, as a function of hole doping concentration $x$, we find that there exists a universal scaling behavior (sample independence) of these ratios with doping, by showing a nonlinearly decreasing behavior of the former; that is, $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}_{c}\ensuremath{\sim}{x}^{\ensuremath{-}\ensuremath{\alpha}}$ with $\ensuremath{\alpha}\ensuremath{\sim}2$ and a near doping independence of the latter; that is, $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}^{*}\ensuremath{\approx}4--6$.

arXiv:cond-mat/0410229•Physical Review B•Apr 26, 2005•7 citations

Superconductivity in the two-dimensional Hubbard model

J. Beenen, David M. Edwards

Quasiparticle bands of the two-dimensional Hubbard model are calculated using the Roth two-pole approximation to the one-particle Green's function. Excellent agreement is obtained with recent Monte Carlo calculations, including an anomalous volume of the Fermi surface near half-filling, which can possibly be explained in terms of a breakdown of Fermi liquid theory. The calculated bands are very flat around the (\ensuremath{\pi},0) points of the Brillouin zone in agreement with photoemission measurements of cuprate superconductors. With doping there is a shift in spectral weight from the upper band to the lower band. The Roth method is extended to deal with superconductivity within a four-pole approximation allowing electron-hole mixing. It is shown that triplet p-wave pairing never occurs. A self-consistent solution with singlet ${\mathit{d}}_{{\mathit{x}}^{2}\mathrm{\ensuremath{-}}{\mathit{y}}^{2}}$-wave pairing is found and optimal doping occurs when the van Hove singularity, corresponding to the flat band part, lies at the Fermi level. Nearest-neighbor antiferromagnetic correlations play an important role in flattening the bands near the Fermi level and in favoring superconductivity. However, the mechanism for superconductivity is a local one, in contrast to spin-fluctuation exchange models. For reasonable values of the hopping parameter the transition temperature ${\mathit{T}}_{\mathit{c}}$ is in the range 10--100 K. The optimum doping ${\mathrm{\ensuremath{\delta}}}_{\mathit{c}}$ lies between 0.14 and 0.25, depending on the ratio U/t. The gap equation has a BCS-like form and 2${\mathrm{\ensuremath{\Delta}}}_{\mathrm{max}}$/${\mathit{kT}}_{\mathit{c}}$\ensuremath{\simeq}4.

arXiv:cond-mat/9503054•Physical review. B, Condensed matter•Nov 1, 1995•111 citations

Superfluid weight vs. superconducting temperature based on a U(1) slave-boson approach to the t-J Hamiltonian

Sung-Sik Lee, Sung‐Ho Suck Salk

Based on an improved U(1) slave-boson approach to the t-J Hamiltonian, we invest igate a relationship between the superfluid weight $n_s/m^*$(the superconducting charge carrier density/the effective mass of the charge carrier) and the superc onducting temperature $T_c$. From the present study we find a linear increase of $n_s/m^*$ with $T_c$ with th e doping concentration in the underdoped region, a saturation around the optimal doping and a decrease in both $n_s/m^*$ and $T_c$ in the overdoped region. Such a trend of the `boomerang' shaped locus in $n_s/m^*$ vs. $T_c$ with increas ing doping concentration from the underdoped to the heavily overdoped region is predicted to be in complete agreement with muon-spin-relaxation measurements. The boomerang behavior is found to occur in correlation with reduction in the sp in singlet pairing(spinon pairing) order in the heavily overdoped region.

arXiv:cond-mat/0005120•May 6, 2000•0 citations

Static and dynamical properties of doped Hubbard clusters

Elbio R. Dagotto, Adriana Moreo, F. Ortolani +2 more

We study the t-J and the Hubbard models at zero temperature using exact-diagonalization techniques on \ensuremath{\surd}10 \ifmmode\times\else\texttimes\fi{} \ensuremath{\surd}10 and 4\ifmmode\times\else\texttimes\fi{}4 sites clusters. Quantum Monte Carlo simulation results on larger lattices are also presented. All electronic fillings have been analyzed for the three models. We have measured equal-time correlation functions corresponding to various types of order (ranging from ``standard'' staggered spin order to more ``exotic'' possibilities like chiral order), as well as various dynamical properties of these models. Upper bounds for the critical hole doping (${\mathit{x}}_{\mathit{c}}$), where long-range antiferromagnetic order disappears, are presented. It was found that ${\mathit{x}}_{\mathit{c}}$ is very small in agreement with experiments for the high-${\mathit{T}}_{\mathit{c}}$ superconductors. For example, in the t-J model, ${\mathit{x}}_{\mathit{c}}$0.08 at J/t=0.4. However, short-distance spin correlations are important up to much higher dopings producing a sharp well-defined spin-wave-like peak in S(q=(\ensuremath{\pi},\ensuremath{\pi}),\ensuremath{\omega}). Regarding the possibility of phase separation in the Hubbard model, we have studied the behavior of the density of particles, 〈n〉, as a function of the chemical potential, using the Lanczos method on a 4\ifmmode\times\else\texttimes\fi{}4 Hubbard cluster, finding no indications of phase separation for any value of U/t. Then, we conclude that the t-J model at small J/t should not phase separate.In order to compare theoretical predictions with photoemission experiments, we evaluated the electronic density of states N(\ensuremath{\omega}), of the Hubbard and t-J models at several doping fractions. We found that upon doping the antiferromagnetic gap is filled for U\ensuremath{\sim}8t or smaller. The chemical potential moves across the insulating gap as one goes from electron to hole doping of the half-filled cluster, in agreement with x-ray absorption experiments but at variance with photoemission experiments. We have also calculated the optical conductivity, ${\mathrm{\ensuremath{\sigma}}}_{1}$(\ensuremath{\omega}), of the Hubbard and t-J models at all dopings on 4\ifmmode\times\else\texttimes\fi{}4 clusters. Results are compared with experiments and the weight of the Drude peak is presented as a function of couplings and dopings. Spectral weight found at small frequencies is associated with the mid-infrared band observed experimentally in ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$${\mathrm{Sr}}_{\mathit{x}}$${\mathrm{CuO}}_{4}$, and with the states filling the insulating gap in photoemission experiments. An overall good agreement with experiments in the normal state was found. Regarding the possibility of superconductivity n these models, we have studied s-, d-, and p-wave pairing correlations. Naively, the d-wave channel seems enhanced near half filling while the extended s-wave channel seems enhanced from half filling up to 40% doping.However, we found that the enhancement comes from short-distance effects and, thus, no numerical indications of superconductivity were found in these models. We emphasize the importance of analyzing the pairing correlations as a function of distance to distinguish between short- and long-distance effects in the susceptibilities. We also observed that spiral order is enhanced at small J/t and low doping. Uniform chiral rder is suppressed by dynamical holes while staggered chiral order may be enhanced, although with a small plaquette order parameter. We conclude with the observation that the simple one-band Hubbard model with intermediate values of U\ensuremath{\sim}8t may account for many of the ``anomalous'' properties of the normal state of the high-${\mathit{T}}_{\mathit{c}}$ superconductors.

Physical review. B, Condensed matter•May 1, 1992•201 citations

Quantum Monte Carlo Simulations and Weak-Coupling Approximations for the Three-Band Hubbard Model

R. Putz, H. Endres, A. Muramatsu +1 more

We study the three-band Hubbard model for the cuprate superconductors using Quantum Monte Carlo simulations and the fluctuation exchange approximation including a self-consistent many-body renormalization of the single-particle propagator. The energy dispersion of the low-lying one-particle excitations is in very good agreement with spectroscopic measurements for a hole doping concentration δ = 0.25 of the CuO2 planes. At the same doping, we examine the spin susceptibility yielding a non-ordered phase with short-range incommensurate fluctuations near the antiferromagnetic wave vector. In the superconducting regime, the interaction vertex for the pairing correlation function in the extended s-wave channel shows a maximum for δ = 0.20, resembling the dependence of transition temperature Tc on doping.

Notes on numerical fluid mechanics•Jan 1, 1994•1 citation

Doping-driven antiferromagnetic insulator-superconductor transition: A quantum Monte Carlo study

Tianxing Ma, Da Wang, Congjun Wu

How superconductivity emerges in the vicinity of an antiferromagnetic insulating state is a long-standing issue of strong correlation physics. We study the transition from an antiferromagnetic insulator to a superconductor by hole doping based on a bilayer generalization of a Hubbard model. The projector quantum Monte Carlo simulations are employed, which are sign problem free both at and away from half filling. An anisotropic Ising antiferromagnetic Mott insulating phase occurs at half filling, which is weakened by hole doping. Below a critical doping value, antiferromagnetism coexists with the singlet superconductivity, which is a pairing across each rung with an extended $s$-wave symmetry. As doping further increases, the antiferromagnetic order vanishes, leaving only a superconducting phase. These results provide important information on how superconductivity appears upon doping the parent Mott-insulating state.

arXiv:1806.03652•Physical review. B./Physical review. B•Aug 12, 2022•10 citations
The phase fluctuations of the condensate in doped antiferromagnets, described by a t\ensuremath−t\ensuremath′\ensuremath−J{t\ensuremath{-}t}^{\ensuremath{'}}\ensuremath{-}Jt\ensuremath−t\ensuremath′\ensuremath−J model and a suitable 1/N1/N1/N expansion, provide a mechanism for a Kosterlitz-Thouless (KT) type of transition to a superconducting state below Tc.{T}_{c}.Tc​. In this paper, we present a Monte Carlo study of the corresponding superfluid weight Ds(T){D}_{s}(T)Ds​(T) in the classical (large-N)N)N) limit, as a function of temperature and doping. Consistent with generic experimental trends, Ds(T){D}_{s}(T)Ds​(T) exhibits a T-linear decrease at low temperatures, with the magnitude of the slope Ds\ensuremath′(0){D}_{s}^{\ensuremath{'}}(0)Ds\ensuremath′​(0) increasing upon doping. Finite-size scaling in the underdoped regime predicts values for the dimensionless ratio A=kBTc/Ds(0){A=k}_{B}{T}_{c}{/D}_{s}(0)A=kB​Tc​/Ds​(0) of order unity, with A=0.4435(5)A=0.4435(5)A=0.4435(5) in the half filled band limit, thus confirming Ds(0){D}_{s}(0)Ds​(0) as the fundamental energy scale determining Tc.{T}_{c}.Tc​. Our Monte Carlo results for Ds(T)/Ds(0){D}_{s}{(T)/D}_{s}(0)Ds​(T)/Ds​(0) vs kBT/Ds(0),{k}_{B}{T/D}_{s}(0),kB​T/Ds​(0), at 10% hole doping, are found to be in reasonable agreement with recent measurements on La2\ensuremath−xSrxCuO4,{\mathrm{La}}_{2\ensuremath{-}x}{\mathrm{Sr}}_{x}{\mathrm{CuO}}_{4},La2\ensuremath−x​Srx​CuO4​, with x=0.10,x=0.10,x=0.10, throughout the temperature range below the theoretical KT transition temperature Tc.{T}_{c}.Tc​.
Using the improved U(1) and SU(2) slave-boson approaches of the $t\ensuremath{-}J$ Hamiltonian [S.-S. Lee and Sung-Ho Suck Salk, Phys. Rev. B 64, 052501 (2001)] that we developed recently, we study the doping and temperature dependence of superfluid weight and spectral function and discuss our finding of the universal scaling behavior of the pseudogap. It is shown that at low hole doping concentrations xxx and at low temperatures TTT, there exists a propensity of linear decrease in the superfluid weight ns∕m∗{n}_{s}∕{m}^{*}ns​∕m∗ with temperature TTT, and a tendency of doping independence in the linearly decreasing slopes of (ns∕m∗)(x,T)({n}_{s}∕{m}^{*})(x,T)(ns​∕m∗)(x,T) with TTT in qualitative agreement with the experimentally observed relation; that is, $({n}_{s}∕{m}^{*})(x,T)=({n}_{s}∕{m}^{*})(x,0)\ensuremath{-}\ensuremath{\alpha}T$, where $\ensuremath{\alpha}$ is a constant. It is also demonstrated that there exists a boomerang behavior; that is, both Tc{T}_{c}Tc​ and ns∕m∗{n}_{s}∕{m}^{*}ns​∕m∗ increase with hole doping concentration xxx in the underdoped region, reach a saturation (maximum) at a hole doping above optimal doping, and decrease beyond the saturation point in the overdoped region, in agreement with $\ensuremath{\mu}\mathrm{SR}$ measurements. Based on our improved SU(2) slave-boson approach to the $t\ensuremath{-}J$ Hamiltonian, we further investigate the doping and temperature dependence of the spectral function and discuss our theoretical finding of a scaling behavior of the pseudogap. In addition, we discuss the cause of the hump and quasiparticle peak in the observed spectral functions of high Tc{T}_{c}Tc​ cuprates. It is demonstrated that the sharpening of the observed quasiparticle peak below Tc{T}_{c}Tc​ is attributed to the bose condensation of holon pair, in agreement with observations. From the predicted ratios of the pseudogap ${\ensuremath{\Delta}}_{0}$ to both the superconducting temperature Tc{T}_{c}Tc​ and the pseudogap temperature T∗{T}^{*}T∗, i.e., $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}_{c}$ and $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}^{*}$, respectively, as a function of hole doping concentration xxx, we find that there exists a universal scaling behavior (sample independence) of these ratios with doping, by showing a nonlinearly decreasing behavior of the former; that is, $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}_{c}\ensuremath{\sim}{x}^{\ensuremath{-}\ensuremath{\alpha}}$ with $\ensuremath{\alpha}\ensuremath{\sim}2$ and a near doping independence of the latter; that is, $2{\ensuremath{\Delta}}_{0}∕{k}_{B}{T}^{*}\ensuremath{\approx}4--6$.
Quasiparticle bands of the two-dimensional Hubbard model are calculated using the Roth two-pole approximation to the one-particle Green's function. Excellent agreement is obtained with recent Monte Carlo calculations, including an anomalous volume of the Fermi surface near half-filling, which can possibly be explained in terms of a breakdown of Fermi liquid theory. The calculated bands are very flat around the (\ensuremath{\pi},0) points of the Brillouin zone in agreement with photoemission measurements of cuprate superconductors. With doping there is a shift in spectral weight from the upper band to the lower band. The Roth method is extended to deal with superconductivity within a four-pole approximation allowing electron-hole mixing. It is shown that triplet p-wave pairing never occurs. A self-consistent solution with singlet ${\mathit{d}}_{{\mathit{x}}^{2}\mathrm{\ensuremath{-}}{\mathit{y}}^{2}}$-wave pairing is found and optimal doping occurs when the van Hove singularity, corresponding to the flat band part, lies at the Fermi level. Nearest-neighbor antiferromagnetic correlations play an important role in flattening the bands near the Fermi level and in favoring superconductivity. However, the mechanism for superconductivity is a local one, in contrast to spin-fluctuation exchange models. For reasonable values of the hopping parameter the transition temperature Tc{\mathit{T}}_{\mathit{c}}Tc​ is in the range 10--100 K. The optimum doping ${\mathrm{\ensuremath{\delta}}}_{\mathit{c}}$ lies between 0.14 and 0.25, depending on the ratio U/t. The gap equation has a BCS-like form and 2${\mathrm{\ensuremath{\Delta}}}_{\mathrm{max}}$/kTc{\mathit{kT}}_{\mathit{c}}kTc​\ensuremath{\simeq}4.
Based on an improved U(1) slave-boson approach to the t-J Hamiltonian, we invest igate a relationship between the superfluid weight ns/m∗n_s/m^*ns​/m∗(the superconducting charge carrier density/the effective mass of the charge carrier) and the superc onducting temperature TcT_cTc​. From the present study we find a linear increase of ns/m∗n_s/m^*ns​/m∗ with TcT_cTc​ with th e doping concentration in the underdoped region, a saturation around the optimal doping and a decrease in both ns/m∗n_s/m^*ns​/m∗ and TcT_cTc​ in the overdoped region. Such a trend of the `boomerang' shaped locus in ns/m∗n_s/m^*ns​/m∗ vs. TcT_cTc​ with increas ing doping concentration from the underdoped to the heavily overdoped region is predicted to be in complete agreement with muon-spin-relaxation measurements. The boomerang behavior is found to occur in correlation with reduction in the sp in singlet pairing(spinon pairing) order in the heavily overdoped region.
We study the t-J and the Hubbard models at zero temperature using exact-diagonalization techniques on \ensuremath{\surd}10 \ifmmode\times\else\texttimes\fi{} \ensuremath{\surd}10 and 4\ifmmode\times\else\texttimes\fi{}4 sites clusters. Quantum Monte Carlo simulation results on larger lattices are also presented. All electronic fillings have been analyzed for the three models. We have measured equal-time correlation functions corresponding to various types of order (ranging from ``standard'' staggered spin order to more ``exotic'' possibilities like chiral order), as well as various dynamical properties of these models. Upper bounds for the critical hole doping (xc{\mathit{x}}_{\mathit{c}}xc​), where long-range antiferromagnetic order disappears, are presented. It was found that xc{\mathit{x}}_{\mathit{c}}xc​ is very small in agreement with experiments for the high-Tc{\mathit{T}}_{\mathit{c}}Tc​ superconductors. For example, in the t-J model, ${\mathit{x}}_{\mathit{c}}$0.08 at J/t=0.4. However, short-distance spin correlations are important up to much higher dopings producing a sharp well-defined spin-wave-like peak in S(q=(\ensuremath{\pi},\ensuremath{\pi}),\ensuremath{\omega}). Regarding the possibility of phase separation in the Hubbard model, we have studied the behavior of the density of particles, 〈n〉, as a function of the chemical potential, using the Lanczos method on a 4\ifmmode\times\else\texttimes\fi{}4 Hubbard cluster, finding no indications of phase separation for any value of U/t. Then, we conclude that the t-J model at small J/t should not phase separate.In order to compare theoretical predictions with photoemission experiments, we evaluated the electronic density of states N(\ensuremath{\omega}), of the Hubbard and t-J models at several doping fractions. We found that upon doping the antiferromagnetic gap is filled for U\ensuremath{\sim}8t or smaller. The chemical potential moves across the insulating gap as one goes from electron to hole doping of the half-filled cluster, in agreement with x-ray absorption experiments but at variance with photoemission experiments. We have also calculated the optical conductivity, ${\mathrm{\ensuremath{\sigma}}}_{1}$(\ensuremath{\omega}), of the Hubbard and t-J models at all dopings on 4\ifmmode\times\else\texttimes\fi{}4 clusters. Results are compared with experiments and the weight of the Drude peak is presented as a function of couplings and dopings. Spectral weight found at small frequencies is associated with the mid-infrared band observed experimentally in ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$Srx{\mathrm{Sr}}_{\mathit{x}}Srx​CuO4{\mathrm{CuO}}_{4}CuO4​, and with the states filling the insulating gap in photoemission experiments. An overall good agreement with experiments in the normal state was found. Regarding the possibility of superconductivity n these models, we have studied s-, d-, and p-wave pairing correlations. Naively, the d-wave channel seems enhanced near half filling while the extended s-wave channel seems enhanced from half filling up to 40% doping.However, we found that the enhancement comes from short-distance effects and, thus, no numerical indications of superconductivity were found in these models. We emphasize the importance of analyzing the pairing correlations as a function of distance to distinguish between short- and long-distance effects in the susceptibilities. We also observed that spiral order is enhanced at small J/t and low doping. Uniform chiral rder is suppressed by dynamical holes while staggered chiral order may be enhanced, although with a small plaquette order parameter. We conclude with the observation that the simple one-band Hubbard model with intermediate values of U\ensuremath{\sim}8t may account for many of the ``anomalous'' properties of the normal state of the high-Tc{\mathit{T}}_{\mathit{c}}Tc​ superconductors.
We study the three-band Hubbard model for the cuprate superconductors using Quantum Monte Carlo simulations and the fluctuation exchange approximation including a self-consistent many-body renormalization of the single-particle propagator. The energy dispersion of the low-lying one-particle excitations is in very good agreement with spectroscopic measurements for a hole doping concentration δ = 0.25 of the CuO2 planes. At the same doping, we examine the spin susceptibility yielding a non-ordered phase with short-range incommensurate fluctuations near the antiferromagnetic wave vector. In the superconducting regime, the interaction vertex for the pairing correlation function in the extended s-wave channel shows a maximum for δ = 0.20, resembling the dependence of transition temperature Tc on doping.
How superconductivity emerges in the vicinity of an antiferromagnetic insulating state is a long-standing issue of strong correlation physics. We study the transition from an antiferromagnetic insulator to a superconductor by hole doping based on a bilayer generalization of a Hubbard model. The projector quantum Monte Carlo simulations are employed, which are sign problem free both at and away from half filling. An anisotropic Ising antiferromagnetic Mott insulating phase occurs at half filling, which is weakened by hole doping. Below a critical doping value, antiferromagnetism coexists with the singlet superconductivity, which is a pairing across each rung with an extended sss-wave symmetry. As doping further increases, the antiferromagnetic order vanishes, leaving only a superconducting phase. These results provide important information on how superconductivity appears upon doping the parent Mott-insulating state.