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

Physical review. B, Condensed matter

Pairing correlations in two-dimensional large-size Cu-O clusters: Nonlocal world-line Monte Carlo algorithm

V. А. Kashurnikov

Published Mar 1, 1996• 7 citations
Publisher

Abstract

A Monte Carlo world-line algorithm for a Cu${\mathrm{O}}_{2}$ plane based on a breakup of the plane into five-site cells, Cu${\mathrm{O}}_{4}$ has been suggested. The algorithm has fast convergency and gives the possibility of investigating the superconducting state symmetry. The calculation of pairing correlation functions was made with an advent of additional imaginary-time slices to the Monte Carlo scheme. The characteristics of a two-dimensional Cu-O cluster with number of sites ${N}_{a}=768$ (16\ifmmode\times\else\texttimes\fi{}16 Cu${\mathrm{O}}_{2}$ cells) were calculated. It has been shown within the chosen interval of parameters of the Emery Hamiltonian (${U}_{d}=6t$, $\ensuremath{\varepsilon}=1\ensuremath{-}3t$, ${U}_{p}=V=0$), temperatures ($T>~0.125t$), and carrier concentrations ($0.7<~x<~1.5$) that (i) the long-range off-diagonal order associated with the superconducting state was not observed in the thermodynamic limit in any pairing channels and (ii) a tendency to divergence in ${s}^{*}$ channels and ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ channels, which is noted as the temperature decreases, is due to the antiferromagnetic ordering rather than the superconducting pairing. If the carrier concentration $x$ is equal to 1.0 (undoped dielectric state) the pairing correlations have maximum amplitude and demonstrate clear antiferromagnetic ordering in the copper sublattice, and the characteristic correlation length is close to the earlier-reported antiferromagnetic length.

DOI
10.1103/physrevb.53.5932

Related papers

Pairing correlation in CuO2 plane: new nonlocal world-line Monte Carlo algorithm

V. А. Kashurnikov

A new nonlocal Monte Carlo world-line algorithm is proposed for the 2D–CuO2 plane based on a breaking up of the plane into five-site CuO4 cells. The algorithm is rapidly convergent and makes it possible to study the symmetry of the superconducting state. It has been used to calculate for the first time the characteristics of the 2D Cu–O cluster with Na = 768(16 × 16CuO2 cells) sites. It is shown that within the chosen interval of Emery Hamiltonian parameters Ud = 3 – 8t, ε = 1 – 3t, Up = V = 0, temperatures T ⩾ 0.125t, and carrier concentration 0.7 ⩽ x ⩽ 1.5, (i) long-range off-diagonal order (indicating the superconducting state) is not observed in any of the pairing channels in the thermodynamic limit, (ii) a tendency towards the divergence of s*- and dx2−y2-pairing correlations, observed upon a decrease in temperature, is due to the antiferromagnetic ordering on the copper sublattice. For the carrier concentration x = 1.0 (insulating state), the pairing correlations have the maximum amplitude, the characteristic correlation length being close to the antiferromagnetic length.

Low Temperature Physics•May 1, 1996•0 citations

Long-range correlations in the two-dimensional spin-(1/2 antiferromagnetic Heisenberg model: A quantum Monte Carlo study

Efstratios Manousakis, Román Salvador

Using an extended version of Handscomb's Monte Carlo method we study the spin-(1/2 quantum antiferromagnetic Heisenberg model on square lattices of various sizes. The calculated correlation length associated with antiferromagnetic order grows very rapidly as the temperature is lowered suggesting an essential singularity. Our results are consistent with those recently observed by neutron scattering done on ${\mathrm{Lo}}_{2}$${\mathrm{CuO}}_{4}$.

Physical Review Letters•Feb 29, 1988•75 citations

Pairing in the two-dimensional Hubbard model: An exact diagonalization study

H. Q. Lin, J. E. Hirsch, Douglas J. Scalapino

We have studied the pair susceptibilities for all possible pair wave functions that fit on a two-dimensional (2D) eight-site Hubbard cluster by exact diagonalization of the Hamiltonian. Band fillings corresponding to four and six electrons were studied (two or four holes in the half-filled band) for a wide range of Hubbard interaction strengths and temperatures. Our results show that all pairing susceptibilities are suppressed by the Hubbard repulsion. We have also carried out perturbation-theory calculations which show that the leading-order ${U}^{2}$ contributions to the d-wave pair susceptibility suppresses d-wave pairing over a significant temperature range. These results are consistent with recent Monte Carlo results and provide further evidence suggesting that the 2D Hubbard model does not exhibit superconductivity.

Physical review. B, Condensed matter•May 1, 1988•80 citations

Quantum Monte Carlo study of hole binding and pairing correlations in the three-band Hubbard model

M. Guerrero, James E. Gubernatis, Shiwei Zhang

We simulated the three-band Hubbard model using the constrained-path Monte Carlo (CPMC) method in search for a possible superconducting ground state. The CPMC is a ground-state method that is free of the exponential scaling of computing time with system size. We calculated the binding energy of a pair of holes for systems up to $6\ifmmode\times\else\texttimes\fi{}4$ unit cells. We also studied the pairing correlation functions versus distance for both the $d$-wave and extended $s$-wave channels in systems up to $6\ifmmode\times\else\texttimes\fi{}6$. We found that holes bind for a wide range of parameters and that the binding increased as the system size is increased. However, the pairing correlation functions decay quickly with distance. For the extended $s$ channel, we found that as the Coulomb interaction ${U}_{d}$ on the Cu sites is increased, the long-range part of the correlation functions is suppressed and fluctuates around zero. For the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ channel, we found that the correlations decay rapidly with distance towards a small positive value. However, this value becomes smaller as the interaction ${U}_{d}$ or the system size is increased.

arXiv:cond-mat/9711125•Physical review. B, Condensed matter•May 15, 1998•55 citations

Detection of pairing correlation in the two-dimensional Hubbard model

Kazuhiko Kuroki, Hideo Aoki, Takashi Hotta +1 more

The quantum Monte Carlo method is used to reexamine superconductivity in the single-band Hubbard model in two dimensions. Instead of the conventional pairing, we consider a ``correlated pairing,'' 〈c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\uparrow}}}$c-${\mathrm{tilde}}_{{\mathrm{i}}^{|\mathrm{I}\mathrm{H}}\mathrm{\ensuremath{\downarrow}}}$〉 with c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$\ensuremath{\equiv}${\mathrm{c}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$(1-${\mathrm{n}}_{\mathrm{i}\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\sigma}}}$), which is inferred from the t-J model, the strong-coupling limit of the Hubbard model. The pairing in the d-wave channel is found to possess both a divergence like 1/T in the pairing susceptibility and a growth of the ground-state pairing correlation with sample size, indicating an off-diagonal long-range order near (but not exactly at) half-filling.

arXiv:cond-mat/9606143•Physical review. B, Condensed matter•Feb 1, 1997•10 citations

Finite-temperature superconducting correlations of the Hubbard model

Ehsan Khatami, Richard T. Scalettar, Rajiv R. P. Singh

We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square-lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the $d$-wave and extended $s$-wave pairing symmetries. The pairing structure factor shows a strong tendency to $d$-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended $s$-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the $d$-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the $d$-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer a unique window into pairing in the Hubbard model at strong coupling.

arXiv:1503.06213•Physical Review B•Jun 18, 2015•24 citations
A Monte Carlo world-line algorithm for a CuO2{\mathrm{O}}_{2}O2​ plane based on a breakup of the plane into five-site cells, CuO4{\mathrm{O}}_{4}O4​ has been suggested. The algorithm has fast convergency and gives the possibility of investigating the superconducting state symmetry. The calculation of pairing correlation functions was made with an advent of additional imaginary-time slices to the Monte Carlo scheme. The characteristics of a two-dimensional Cu-O cluster with number of sites Na=768{N}_{a}=768Na​=768 (16\ifmmode\times\else\texttimes\fi{}16 CuO2{\mathrm{O}}_{2}O2​ cells) were calculated. It has been shown within the chosen interval of parameters of the Emery Hamiltonian (Ud=6t{U}_{d}=6tUd​=6t, \ensuremathε=1\ensuremath−3t\ensuremath{\varepsilon}=1\ensuremath{-}3t\ensuremathε=1\ensuremath−3t, Up=V=0{U}_{p}=V=0Up​=V=0), temperatures (T> 0.125tT>~0.125tT> 0.125t), and carrier concentrations (0.7< x< 1.50.7<~x<~1.50.7< x< 1.5) that (i) the long-range off-diagonal order associated with the superconducting state was not observed in the thermodynamic limit in any pairing channels and (ii) a tendency to divergence in s∗{s}^{*}s∗ channels and dx2\ensuremath−y2{d}_{{x}^{2}\ensuremath{-}{y}^{2}}dx2\ensuremath−y2​ channels, which is noted as the temperature decreases, is due to the antiferromagnetic ordering rather than the superconducting pairing. If the carrier concentration xxx is equal to 1.0 (undoped dielectric state) the pairing correlations have maximum amplitude and demonstrate clear antiferromagnetic ordering in the copper sublattice, and the characteristic correlation length is close to the earlier-reported antiferromagnetic length.
A new nonlocal Monte Carlo world-line algorithm is proposed for the 2D–CuO2 plane based on a breaking up of the plane into five-site CuO4 cells. The algorithm is rapidly convergent and makes it possible to study the symmetry of the superconducting state. It has been used to calculate for the first time the characteristics of the 2D Cu–O cluster with Na = 768(16 × 16CuO2 cells) sites. It is shown that within the chosen interval of Emery Hamiltonian parameters Ud = 3 – 8t, ε = 1 – 3t, Up = V = 0, temperatures T ⩾ 0.125t, and carrier concentration 0.7 ⩽ x ⩽ 1.5, (i) long-range off-diagonal order (indicating the superconducting state) is not observed in any of the pairing channels in the thermodynamic limit, (ii) a tendency towards the divergence of s*- and dx2−y2-pairing correlations, observed upon a decrease in temperature, is due to the antiferromagnetic ordering on the copper sublattice. For the carrier concentration x = 1.0 (insulating state), the pairing correlations have the maximum amplitude, the characteristic correlation length being close to the antiferromagnetic length.
Using an extended version of Handscomb's Monte Carlo method we study the spin-(1/2 quantum antiferromagnetic Heisenberg model on square lattices of various sizes. The calculated correlation length associated with antiferromagnetic order grows very rapidly as the temperature is lowered suggesting an essential singularity. Our results are consistent with those recently observed by neutron scattering done on Lo2{\mathrm{Lo}}_{2}Lo2​CuO4{\mathrm{CuO}}_{4}CuO4​.
We have studied the pair susceptibilities for all possible pair wave functions that fit on a two-dimensional (2D) eight-site Hubbard cluster by exact diagonalization of the Hamiltonian. Band fillings corresponding to four and six electrons were studied (two or four holes in the half-filled band) for a wide range of Hubbard interaction strengths and temperatures. Our results show that all pairing susceptibilities are suppressed by the Hubbard repulsion. We have also carried out perturbation-theory calculations which show that the leading-order U2{U}^{2}U2 contributions to the d-wave pair susceptibility suppresses d-wave pairing over a significant temperature range. These results are consistent with recent Monte Carlo results and provide further evidence suggesting that the 2D Hubbard model does not exhibit superconductivity.
We simulated the three-band Hubbard model using the constrained-path Monte Carlo (CPMC) method in search for a possible superconducting ground state. The CPMC is a ground-state method that is free of the exponential scaling of computing time with system size. We calculated the binding energy of a pair of holes for systems up to $6\ifmmode\times\else\texttimes\fi{}4$ unit cells. We also studied the pairing correlation functions versus distance for both the ddd-wave and extended sss-wave channels in systems up to $6\ifmmode\times\else\texttimes\fi{}6$. We found that holes bind for a wide range of parameters and that the binding increased as the system size is increased. However, the pairing correlation functions decay quickly with distance. For the extended sss channel, we found that as the Coulomb interaction Ud{U}_{d}Ud​ on the Cu sites is increased, the long-range part of the correlation functions is suppressed and fluctuates around zero. For the ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ channel, we found that the correlations decay rapidly with distance towards a small positive value. However, this value becomes smaller as the interaction Ud{U}_{d}Ud​ or the system size is increased.
The quantum Monte Carlo method is used to reexamine superconductivity in the single-band Hubbard model in two dimensions. Instead of the conventional pairing, we consider a ``correlated pairing,'' 〈c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\uparrow}}}$c-${\mathrm{tilde}}_{{\mathrm{i}}^{|\mathrm{I}\mathrm{H}}\mathrm{\ensuremath{\downarrow}}}$〉 with c-${\mathrm{tilde}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$\ensuremath{\equiv}${\mathrm{c}}_{\mathrm{i}\mathrm{\ensuremath{\sigma}}}$(1-${\mathrm{n}}_{\mathrm{i}\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\sigma}}}$), which is inferred from the t-J model, the strong-coupling limit of the Hubbard model. The pairing in the d-wave channel is found to possess both a divergence like 1/T in the pairing susceptibility and a growth of the ground-state pairing correlation with sample size, indicating an off-diagonal long-range order near (but not exactly at) half-filling.
We utilize numerical linked-cluster expansions (NLCEs) and the determinantal quantum Monte Carlo algorithm to study pairing correlations in the square-lattice Hubbard model. To benchmark the NLCE, we first locate the finite-temperature phase transition of the attractive model to a superconducting state away from half filling. We then explore the superconducting properties of the repulsive model for the ddd-wave and extended sss-wave pairing symmetries. The pairing structure factor shows a strong tendency to ddd-wave pairing and peaks at an interaction strength comparable to the bandwidth. The extended sss-wave structure factor and correlation length are larger at higher temperatures but clearly saturate as temperature is lowered, whereas the ddd-wave counterparts, which start off lower at high temperatures, continue to rise near half filling. This rise is even more dramatic in the ddd-wave susceptibility. The convergence of NLCEs breaks down as the susceptibilities and correlation lengths become large, so we are unable to determine the onset of long-range order. However, our results extend the conclusion, previously restricted to only magnetic and charge correlations, that NLCEs offer a unique window into pairing in the Hubbard model at strong coupling.