## Dominant p -wave pairing induced by nearest-neighbor attraction in the square-lattice extended Hubbard model

Zhangkai Cao, 1, ∗ Jianyu Li, 1, 2, ∗ Jiahao Su, 1, 2 Tao Ying, 3, † and Ho-Kin Tang 1, 2, ‡

1

School of Science, Harbin Institute of Technology, Shenzhen, 518055, China 2 Shenzhen Key Laboratory of Advanced Functional Carbon Materials Research and Comprehensive Application, Shenzhen 518055, China. 3 School of Physics, Harbin Institute of Technology, Harbin 150001, China

(Dated: January 20, 2025)

The two-dimensional (2D) Hubbard model is widely believed to capture the key ingredients of high-temperature superconductivity in cuprate materials. Here, we report a constrained path quantum Monte Carlo study of the square-lattice extended Hubbard model with on-site Coulomb repulsion U and nearest-neighbor (NN) electron attraction V . Upon doping δ = 0.125, we find that NN electron attraction V significantly promotes an exotic spin-triplet ( p -wave) pairing, with the p -wave pairing correlations growing stronger as V increases. In the meanwhile, the d x 2 -y 2 -wave ( d -wave) pairing in the intermediate coupling regime shows insignificant response to the increase of V . The spin density wave order, which is observed near half-filling in the particle-hole channel, is also almost unaffected by the NN electron attraction V , reinforcing the established relationship between d -wave superconductivity (SC) and spin order. As the doping becomes heavier (i.e., δ increases), the dominant region of p -wave pairing expands, further suppressing the presence of d -wave pairing. Our work suggests the p -wave pairing region can be induced and further broadened by the NN electron attraction V , offering a feasible mechanism to realize p -wave SC in realistic cuprate materials.

In almost four decades since the discovery of copperbased high-temperature ( T c ) superconductors [1], later identified with a d x 2 -y 2 -wave ( d -wave) pairing symmetry [2], reformed our understanding of strongly correlated electron systems, many of which defy Landau's Fermiliquid theory [3]. Lots of studies have pointed to the possibility of a phase transition [4] and symmetry breaking [5] especially in the underdoped region, resulting in many interweaving orders, including various forms of charge density wave (CDW) [6], spin density wave (SDW) [7], pair density wave (PDW) [8], electron nematic order [9]. The understanding of how these orders including superconducting, compete and cooperate with each other, has been the most significant problem in condensed matter physics for several decades remaining substantially unsolved [10, 11].

pairing order through various numerical simulation methods, reflecting the Hubbard model and its cousin t -J model exhibit extreme sensitivity to ground state configurations and low-lying excitations, as well as the competition between d -wave SC and other orders [20-25].

Although the ultimate theory for highT c superconductivity remains controversial, the strong electron repulsive interactions in Copper 3d orbitals are believed to play crucial roles [12]. The two-dimensional (2D) fermionic Hubbard model [13-16] successfully captures most of the important physics, like antiferromagnetism at half filling [17] and competing orders in doped cuprates [11, 18, 19]. However, the search for a d -wave superconductivity (SC) phase in the pure 2D Hubbard Model has not been quite so successful [20], which is still a lot of controversy. Recently, dome-like SC region is found in both the electronand hole-doped regimes of the 2D Hubbard model with next-nearest-neighbor (NNN) hopping t ′ [21]. Both positive and negative results have been observed for d -wave

∗ These authors contributed equally.

‡ denghaojian@hit.edu.cn

† taoying86@hit.edu.cn

The relation between superconductivity and other orders is also strongly influenced by the variations of the Hubbard model, such as the inclusion of the nearest-neighbor (NN) attractive interaction V . Moreover, recent angle-resolved photoemission spectroscopy (ARPES) experiments on one-dimensional (1D) cuprate chain Ba 2 -x Sr x CuO 3+ δ [26] have revealed an anomalously sizable attractive interaction between NN electrons, possibly mediated by phonons [26-28]. Such an effective attraction, which has been largely overlooked in previous studies, could play a crucial role in understanding highT c superconductivity. Although not so strong as the on-site Coulomb repulsion, this effective attractive interaction is comparable to the electron hopping integral ( V ∼ -t ), and thereby also should not be ignored in 2D systems [27, 29].

The effect of the NN electron attraction V in the 1D extended Hubbard chain has been examined recently [30], which can host dominant spin-triplet ( p -wave) pairing correlations and divergent superconductive susceptibility. This offers a feasible mechanism for realizing exotic p -wave SC in realistic cuprate chains. Beyond 1D, the density matrix renormalization group (DMRG) study of the extended Hubbard model on long four-leg cylinders on the square lattice finds that the NN electron attraction V can notably enhance the long-distance superconducting correlations while simultaneously suppressing the CDW correlations [29]. Numerically exact diagonalization study of the Hamiltonian is currently restricted to 4 × 4 systems [31], so whether the prediction of a p -

wave paired state could be stabilized in larger systems is doubtful.

Moreover, the proposed competition between spinsinglet ( s -wave and d -wave) pairing and spin-triplet ( p -wave) pairing warrants further investigation [32-35]. Although multiple possible triplet superconductors have been find, such as Sr RuO 2 4 [36, 37], UPt 3 [38], UTe 2 [39], WTe 2 [40], and Cu x Bi Se 2 3 [41], the pairing symmetry and microscopic superconducting mechanism of these materials are not yet completely understood. In the majority of unconventional superconductors, such as copper oxide [6, 10] and iron-based [42] highT c superconductors, it is generally believed that the electron Cooper pairs exist in a spin-singlet state with a total spin S = 0. Spin-triplet SC with S = 1 possesses internal properties, leading to rich physics such as further possible spinrotation symmetry breaking, collective modes of the order parameter and multiple phases of the condensate [3641]. Recently, the electrons with the same spin at the domain walls in the 1-UC Fe(T e , S )/STO were found to e form non-zero momentum Cooper pairs, thereby generating incommensurate PDW state [42]. This exotic triplet equal-spin pairing state can provide a new material platform to study the PDW state and its interplay with the topological electronic states and unconventional superconductivity. For p -wave superconductors, the disruption of time-reversal symmetry can give rise to topological superconductivity, which holds promising potential for applications in topological quantum computing [43, 44].

Considering the structural and quantum chemistry similarities among cuprates, the NN attraction effect should not be ignored in 2D systems. Therefore, investigating the 2D extended Hubbard model with NN attraction could be crucial for understanding the highT c pairing mechanism. Here, we using the constrained path quantum Monte Carlo (CPQMC) [45, 46] to systematically study the phase diagrams of the square-lattice extended Hubbard model with on-site Coulomb repulsion U and NN attraction V of doped cuprates, so called t -U -V model. CPQMC partially solves sign problem concerning only the ground state of the system at zero temperature. Firstly, we give a general phase diagram scan of the U -V space at the doping δ = 0.125, where the d -wave and exotic p -wave triplet pairing phases exist. We find that V can notably drive and enhance the exotic p -wave spintriplet pairing phase, which also leads to the competition between d -wave and p -wave pairing. Besides the pairing phase, a SDW order exists near the half-filling in the particle-hole channel. Moreover, as doping increases, the dominant region of p -wave pairing expands, further suppressing the presence of d -wave pairing. More technical details and additional results are included in the Supplemental Material [47].

We employ the CPQMC to investigate the groundstate properties of the hole-doped Hubbard model on the square lattice defined by the extended Hamiltonian

$$\text{is} \quad & \text{square lattice channel by the extended Hamiltonian} \\ H = & - \sum _ { i, j, \sigma } \left ( t _ { i, j } c _ { i, \sigma } ^ { \dagger } c _ { j, \sigma } + h. c. \right ) + U \sum _ { i } n _ { i, \uparrow } n _ { i, \downarrow } \\ p _ { i, j, \sigma, \sigma ^ { \prime } } ^ { \quad + V \sum _ { i, j, \sigma, \sigma ^ { \prime } } n _ { i, \sigma } n _ { j, \sigma ^ { \prime } } } \\ \text{e} _ { 2 } \quad & \text{,} \quad \text{$+$} \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdots$$

where c † i,σ ( c i,σ ) is electron creation (annihilation) operator with spin σ = ↑ ↓ , , and n i,σ = c † i,σ c i,σ is the electron number operator. We give an illustration of the t -U -V Hubbard model on 2D square lattice (see Fig. 1(a)), where the electron hopping amplitude t i,j = t if i and j are the NN sites, and t i,j = t ′ for NNN sites. Here, we assume t = 1 for simplicity. We use t ′ = -0 2 . for hole doping, according to band-structure calculations in cuprates [48, 49]. U &gt; 0 is the on-site Coulomb repulsion and V &lt; 0 is the NN electron attraction. We concentrate solely on the hole doped case, where the hole doping is given via δ = 1 -N /N e = 1 -n , where N e is the number of electrons in the system and N is the number of lattice sites.

We have defined the effective pair momentum distribution function

$$\text{$m$} _ { \substack { \text{$row$} _ { \substack { n$row$} _ { \substack { n$colon$row$} \\ \text{$the$} _ { \substack { n$colon$row$} \\ \text{$per$-} \\ \text{$i,j$} } } } } ( \mathbf k ) = ( 1 / N ) \sum _ { i, j } \exp [ i \mathbf k ( \mathbf r _ { i } - \mathbf r _ { j } ) ] C _ { \zeta - p a i r } ^ { \text{eff} } ( i, j ), \\ \text{$su$-} \\ \text{ad} -$$

FIG. 1. (Color online) Model and phase diagram. (a) Illustration the t -U -V Hubbard model with NN hopping t , NNN hopping t ′ , on-site repulsive U , and NN attractive V terms. (b) Schematic zero-temperature phase diagram at the underdoping case δ = 0.125, there are conventional d -wave and exotic p -wave triplet pairing phases. (c) The effective d -wave pair momentum distribution function N eff d x 2 -y 2 -pair ( k ) at U = 3 and V = -0 3, . and (d) The effective p -wave pair momentum distribution function N eff p -pair ( k ) at U = 2 and V = -0 9 on 12 . × 12 lattice.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-692573e7-1fea-4335-8ac4-fa05b965c41a.jpg)

where C eff ζ -pair is the effective real-space correlation of the different pairing, and ζ = s, d x 2 -y 2 , d xy , p . Further detailed definitions of correlation functions can be found in the Supplemental Material [47].

We focus on the electron density n =0.875 corresponding to the underdoping δ = 0.125, which is optimal for pairing to systematically explore its phase diagram. At zero temperature, we summarize our main findings in the phase diagram of the square-lattice extended Hubbard model in Fig. 1(b), where the conventional d -wave and exotic p -wave triplet pairing phases in the ground state is uncovered. At δ = 0.125, the d x 2 -y 2 -wave pairing phase dominates the system when there is no attraction V . As V increases, a slight enhancement of NN electron attraction V can drive the system to the p -wave triplet pairing phase, especially at relatively small U . But in the slightly larger U region, a larger attraction V is needed to induce the generation of p -wave pairing to suppress d -wave pairing. Notably, in specific parameter regimes, fermions exhibit the potential to form pairs with exotic p -wave pairing symmetries, further contributing to this intricate pairing interplay [30].This suggests the p -wave pairing order region is induced and further broadened by the NN electron attraction V .

Then we analyze the characteristic of the effective d -wave and p -wave pairing correlation function in momentum space. In Fig. 1(c), we show the effective d -wave pairing correlation function N eff d x 2 -y 2 -pair ( k ) at δ = 0.125 with U = 3 and V = -0 3, . which exhibits a singular nonzero condensation structure, and its maximum value occurs near nonzero momentum point Q = (2 π/ , 3 π ). This uneven distribution of d -wave pairing in momentum space is reminiscent of the PDW order, which is an inhomogeneous superconducting state whose Cooper pairs poss a finite momentum [8, 10, 11]. In other words, we demonstrate the presence of the d -wave PDW phase. Then, we focus on p -wave pairing state, which refers to the pairing between the same spin species. In Fig. 1(d), we show the effective p -wave pairing correlation function N eff p -pair ( k ) at δ = 0.125 with U = 2 and V = -0 9, which . condenses at point Q = (0 , 0) and locates in the dominant region of weak-coupling U and larger V , supporting the formation of Cooper-pair triplets with zero center-ofmass momentum and p -wave symmetry. We provide the effective d -wave and p -wave pairing correlation function in momentum space with the change of U V , and δ in the Supplemental Material [47]. In the region where d -wave pairing dominates, as discussed above, the d -wave pairing correlation function displays a singular nonzero condensation structure, indicating the presence of a d -wave PDWphase. In contrast, in the region where p -wave pairing dominates, the d -wave pairing gradually evolves into a Q = (0 , 0) condensation and has a lower strength than that of the p -wave pairing, with the latter consistently condenses at Q = (0 , 0). The p -wave pairing strength increases as V is raised, while it is suppressed by the increase of U . As doping increases, we observe that the d -wave pairing condensation point becomes more dispersed,

FIG. 2. (Color online) The competition and distinction of d -wave and p -wave pairing at δ = 0.125 on 12 × 12 lattice. (a)-(d) The effective pair momentum distribution function N eff d x 2 -y 2 -pair ( k ), N eff p -pair ( k ) and N eff d xy -pair ( k ) as a function of V , and fix (a) U = 1.0, (b) U = 2.0, (c) U = 3.0, (d) U = 4.0. (e) The effective correlation function of d -wave pairing mode in real space in the x and y direction with U = 3 and V = -0 3, which displays a different staggered behavior . with distance. (f) The effective correlation function of p -wave pairing mode in real space with U = 2 and V = -0 9, which . displays an exponential decay with distance.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-3d14901b-853d-49b9-9350-e0ae3bad3fa6.jpg)

with its intensity decreasing, while the p -wave pairing continues to condense at (0 , 0), maintaining a relatively stable intensity.

To show the competition and distinction of d -wave and p -wave pairing order, we present simulations along the vertical path in the phase diagram at δ = 0.125, namely, fixing U = 1 , 2 , 3 , 4 as shown in Fig. 2(a)-(d). According to the results, the d -wave pairing exhibits different behaviors depending on the strength of U . We found that in the absence of V , the strength of the d -wave pairing in the system increases with the increase of U . As V is enhanced, the strength of the d -wave pairing shows two types of changes: (i) in the weak coupling regime ( U = 1 , 2), the d -wave pairing slowly increases with the increase of V ; (ii) in the intermediate coupling regime ( U = 3 , 4), the d -wave pairing remains relatively unchanged with the variation of V . Although investigating stronger U may provide further insight, such simulations are challenging for the CPQMC method. Conse-

quently, our analysis primarily focuses on the intermediate U regime, where the Coulomb interaction is already large enough for correlated systems to prohibit two electrons from occupying the same site, and is believed to capture the behavior of cuprate superconductors.

The p -wave pairing is entirely absent when V = 0, indicating that the p -wave pairing is solely induced by the NNelectron attraction V . In the weakly coupling regime, the p -wave pairing strength increases sharply with the increase of V and quickly surpasses the strength of the d -wave pairing. In the intermediate coupling regime, the p -wave pairing is suppressed until V reaches a certain strength after which it rapidly increases. Overall, the region dominated by p -wave pairing gradually expands and the strength of p -wave pairing gets significantly enhanced as V increases. The competition between d -wave and p -wave pairing states at different V is primarily driven by the variation of the strength of the p -wave pairing. We also investigated the effective d xy -wave pairing correlation function and found that the strength of d xy -wave remains relatively weak at δ = 0.125, suggesting that the main competition is between the d -wave and p -wave pairing in the system.

Then we analyze the characteristic of d -wave and p -wave pairing correlation function in real space. In Fig. 2(e), we show the effective d -wave correlation function in real space, which displays different staggered behaviors in the x and y direction with distance. Such different staggered behaviors represent the uneven distribution of d -wave pairing in real space, which is consistent with the d -wave pairing condensation at nonzero momentum point Q = (2 π/ , 3 π ) in momentum space, indicating that the d -wave pairing order spontaneously breaks the C 4 rotational symmetry of the system. This may be related to the origin of electron nematic order, which corresponds to that the electronic structure preserves the translation symmetry but breaks the rotation symmetry of the CuO 2 plane. On the contrary, the real space p -wave pairing correlation function exhibits an exponential decay as a function of the distance and drops to 0 quickly, indicating that the p -wave pairing correlation is short-range, as shown in Fig. 2(f).

We also studied the spin density correlation in the particle-hole channel, which reflects the magnetic property of the system. We defined the spin structure factor as N S ( k ) = (1 /N ) ∑ i,j exp[ i k r ( i -r j )] 〈 S i · S j 〉 , and S i the spin operator at site i . At half filling or near half filling, the attractive Hubbard model is dominated by the ( π , π ) CDW. In contrast, a predominant positive U &gt; 0 is expected to support the ( π , π ) SDW. We analyze the characteristic of spin structure factor N S ( k ) in momentum space at δ = 0.125 with U = 4 and V = -1 0, . as shown in Fig. 3(a). The SDW order, characterized by the momentum ( π , π ), leds to a long-range order of staggered spin density in real space. To clearly show the impact of V on SDW, we present the simulations of N S ( k ) at δ = 0.125 as a function of V , namely, fixing U = 1 , 2 , 3 , 4 as shown in Fig. 3(b). The strength of spin

FIG. 3. (Color online) Spin density correlation. (a) The spin structure factor N S ( k ) at δ =0.125 with U = 4 and V = -1 0 . on 12 × 12 lattice. N S ( k ) is peaked at the antiferromagnetic ordering momentum, k max = ( π , π ). (b) N S ( k max ) at δ = 0.125 as a function of V , with U = 1.0, 2.0, 3.0, 4.0. Scaled spin structure factor N S ( k max ) L 2 β ν /N for different lattice sizes L = 8, 10, 12 and 14, with (c) V = 0 and (d) V = -1 0 at . δ ∼ 0.125. The crossing of the different data sets at U c /similarequal 1.7 and U c /similarequal 2.6 respectively, approximately gives the critical interaction U c identifying the phase transition. (Inset) Additionally, when scaling the onsite interaction U by ( U -U c ) L 1 ν /U c , the data sets for different system sizes collapse for critical exponents ν = 1.02(2) and β = 0.85(2).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-8c1d9481-9b2a-4c08-bd71-3b77fa831e4a.jpg)

correlation increases with the increase of U , which is consistent with the variation of d -wave pairing with U . As V is enhanced, the strength of spin correlation remains basically unchanged at small V , but slightly decreases at large V . Such behavior indicates that the spin correlation is not significantly affected by the change of V , which is again similar to the behavior of the d -wave pairing in the intermediate coupling regime, to some extent proves the validness of the antiferromagnetic spin fluctuations as the pairing glue of the d -wave superconductivity [50, 51].

As demonstrated in Fig. 3(c) and (d), we can approximately determined U c from the crossing point of the scaled spin structure factor N S ( k max ) L 2 β ν /N obtained for different lattice sizes L at δ ∼ 0.125 with V = 0 and V = -1 0, separately. . The critical value U = U c represents a quantum phase transition of the system from a metal phase to the SDW phase. The critical exponents ν and β in the scaling law are determined by collapsing the scaled N S ( k ) arising from different sizes (inset of Fig. 3(c) and (d)), consistent with previous research in squarelattice [52]. We found that at V = 0 and V = -1 0, the . critical interactions are U c /similarequal 1 7 and . U c /similarequal 2.4 respectively, indicating that, in the presence of an attractive V , the system requires a larger critical U to cause the phase transition.

FIG. 4. (Color online) Schematic zero-temperature phase diagram of different pairing modes under different doping: (a) δ = 0.014, (b) δ = 0.208, (c) δ = 0.319 on 12 × 12 lattice, d x 2 -y 2 -wave, d xy -wave and exotic p -wave triplet pairing phases are observed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-6e09615f-8a30-4711-bad5-42632e43a145.jpg)

Lastly, we explore the effect of doping on the pairing order. Three additional doping concentrations are considered, namely a lightly doping δ = 0.014 (near half filling), an overdoped case δ = 0.208, and a heavily overdoped case δ = 0.319. The zero-temperature phase diagrams are displayed in Fig. 4, with d -wave, d xy -wave and exotic p -wave triplet pairing phases observed. At near half filling ( δ = 0.014, Fig. 4(a)), underdoped ( δ = 0.125, Fig. 1b) and slightly overdoped ( δ = 0.208, Fig. 4(b)), the d -wave pairing phase dominates when there is no attraction V . However, in heavily overdoped case ( δ = 0.319, Fig. 4(c)), the system becomes dominated by d xy -wave pairing phase. This is consistent with Fig. S4, where we give the zero-temperature phase diagram of the dominant pairing channel in U -n space with V =0 [47].

Notably, there is no p -wave pairing phase in the absence of V for all the doping cases. However, when V presents, the p -wave triplet pairing order begins to competes with the d -wave pairing order (or the d xy -wave pairing for the heavily doping case δ = 0.319), and overwhelms with an appropriate increase of V . Moreover, as the doping increases, the dominant region of p -wave pairing expands, further suppressing the presence of d -wave (or the d xy -wave) pairing. Our results indicate that the p -wave pairing order is induced by the NN electron attraction V and further benefits when the doping is heavier. From the known results obtained in 1D cuprate chains,

- [1] J. G. Bednorz and K. A. M¨ller, Possible high u T c superconductivity in the Ba -La -Cu -O system, Z Physik B 64 , 189 (1986).
- [2] C. C. Tsuei and J. R. Kirtley, Pairing symmetry in cuprate superconductors, Reviews of modern physics 72 , 969 (2000).
- [3] L. D. Landau, The theory of a fermi liquid, Sov. Phys. JETP 3 , 920 (1957).
- [4] A. Shekhter, B. J. Ramshaw, R. Liang, W. N. Hardy, D. A. Bonn, F. F. Balakirev, R. D. McDonald, J. B. Betts, S. C. Riggs, and A. Migliori, Bounding the pseudogap with a line of phase transitions in YBa 2 Cu O 3 6+ δ ,

the critical strengths of V for entering the p -wave pairing phase can be determined analytically at quarter filling in the U → ∞ limit [30, 53, 54], in which the Luttinger parameter K ρ = 1 [2 + (4 / /π )sin -1 ( V/ 2)] exceeding the unity indicates that it has entered the p -wave pairing phase. In a 2D system, similar mechanism is probable to happen that the strong NN electron attraction V drives the system towards the p -wave pairing phase.

In conclusion, we adopted CPQMC calculations of the 2D extended Hubbard model with the NN attraction V to study the impact of the additional attractive V on the pairing phase and density wave state. The quantum phase diagram containing d -wave, d xy -wave and p -wave pairing phases at different doping is explored. We find that the p -wave pairing phase is completely induced by V , and the p -wave pairing strength is enhanced as V increases. However, in the intermediate coupling regime, the d -wave pairing shows insignificant response to the increase of V . The spin density wave order, which is observed near half-filling in the particle-hole channel, is also almost unaffected by the NN electron attraction V , showing a similar behavior to the d -wave pairing with V . In a way, this supports that the spin fluctuations act as the dominant microscopic mechanism to drive the unconventional d -wave pairing phase. As δ increases, the dominant region of p -wave pairing expands, further suppressing the presence of d -wave pairing. Our work suggests the p -wave pairing region can be induced and further broadened by the NN electron attraction V , which will not only deepen the understanding of the highT c pairing microscopical mechanism, but also offering a feasible mechanism to realize p -wave triplet SC [55, 56] in realistic doped cuprates.

We acknowledge useful discussion with Ji Liu. This work is supported by the National Natural Science Foundation of China (Grant No. 12204130), Shenzhen StartUp Research Funds (Grant No. HA11409065), HITSZ Start-Up Funds (Grant No. X2022000), Shenzhen Key Laboratory of Advanced Functional CarbonMaterials Research and Comprehensive Application (Grant No. ZDSYS20220527171407017). T.Y. acknowledges supports from Natural Science Foundation of Heilongjiang Province (No. YQ2023A004).

## Nature 498 , 75 (2013).

- [5] J. Xia, E. Schemm, G. Deutscher, S. A. Kivelson, D. A. Bonn, W. N. Hardy, R. Liang, W. Siemons, G. Koster, M. M. Fejer, and A. Kapitulnik, Polar Kerr-effect measurements of the high-temperature YBa Cu O 2 3 6+ x superconductor: evidence for broken symmetry near the pseudogap temperature, Physical review letters 100 , 127002 (2008).
- [6] R. Comin, A. Frano, M. M. Yee, Y. Yoshida, H. Eisaki, E. Schierle, E. Weschke, R. Sutarto, F. He, A. Soumyanarayanan, Y. He, M. Le Tacon, I. S. Elfimov, J. E. Hoffman, G. A. Sawatzky, B. Keimer, and A. Dama-

- scelli, Charge order driven by Fermi-arc instability in Bi Sr 3 2 -x La CuO x 6+ δ , Science 343 , 390 (2014).
- [7] E. G. Moon and S. Sachdev, Competition between spin density wave order and superconductivity in the underdoped cuprates, Physical review. B, Condensed matter 80 , 035117 (2009).
- [8] Z. Shi, P. G. Baity, J. Terzic, T. Sasagawa, and D. Popovi´ c, Pair density wave at high magnetic fields in cuprates with charge and spin orders, Nature communications 11 , 3323 (2020).
- [9] Y. Sato, S. Kasahara, H. Murayama, Y. Kasahara, E.G. Moon, T. Nishizaki, T. Loew, J. Porras, B. Keimer, T. Shibauchi, and Y. Matsuda, Thermodynamic evidence for a nematic phase transition at the onset of the pseudogap in YBa Cu O , Nature physics 2 3 y 13 , 1074 (2017).
- [10] B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, and J. Zaanen, From quantum matter to high-temperature superconductivity in copper oxides, Nature 518 , 179 (2015).
- [11] E. Fradkin, S. A. Kivelson, and J. M. Tranquada, Colloquium: Theory of intertwined orders in high temperature superconductors, Reviews of modern physics 87 , 457 (2015).
- [12] P. Mai, G. Balduzzi, S. Johnston, and T. A. Maier, Orbital structure of the effective pairing interaction in the high-temperature superconducting cuprates, npj Quantum Materials 6 , 1 (2021).
- [13] F. C. Zhang and T. M. Rice, Effective Hamiltonian for the superconducting Cu oxides, Physical review. B, Condensed matter 37 , 3759 (1988).
- [14] D. P. Arovas, E. Berg, S. A. Kivelson, and S. Raghu, The Hubbard model, Annual review of condensed matter physics 13 , 239 (2022).
- [15] The Hubbard model at half a century, Nature physics 9 , 523 (2013).
- [16] M. Qin, T. Sch¨ afer, S. Andergassen, P. Corboz, and E. Gull, The hubbard model: A computational perspective, Annual Review of Condensed Matter Physics 13 , 275 (2022).
- [17] Simons Collaboration on the Many-Electron Problem, J. P. F. LeBlanc, A. E. Antipov, F. Becca, I. W. Bulik, G. K.-L. Chan, C.-M. Chung, Y. Deng, M. Ferrero, T. M. Henderson, C. A. Jim´ enez-Hoyos, E. Kozik, X.W. Liu, A. J. Millis, N. V. Prokof'ev, M. Qin, G. E. Scuseria, H. Shi, B. V. Svistunov, L. F. Tocchio, I. S. Tupitsyn, S. R. White, S. Zhang, B.-X. Zheng, Z. Zhu, and E. Gull, Solutions of the Two-Dimensional Hubbard Model: Benchmarks and Results from a Wide Range of Numerical Algorithms, Physical Review X 5 , 041041 (2015).
- [18] D. J. Scalapino, A common thread: The pairing interaction for unconventional superconductors, Reviews of modern physics 84 , 1383 (2012).
- [19] B.-X. Zheng, C.-M. Chung, P. Corboz, G. Ehlers, M.-P. Qin, R. M. Noack, H. Shi, S. R. White, S. Zhang, and G. K.-L. Chan, Stripe order in the underdoped region of the two-dimensional Hubbard model, Science 358 , 1155 (2017).
- [20] Simons Collaboration on the Many-Electron Problem, M. Qin, C.-M. Chung, H. Shi, E. Vitali, C. Hubig, U. Schollw¨ck, S. R. White, and S. Zhang, Absence of Suo perconductivity in the Pure Two-Dimensional Hubbard Model, Physical Review X 10 , 031016 (2020).
- [21] H. Xu, C.-M. Chung, M. Qin, U. Schollw¨ck, S. R. White, o
- and S. Zhang, Coexistence of superconductivity with partially filled stripes in the Hubbard model, Science 384 , eadh7691 (2024).
- [22] A. Himeda, T. Kato, and M. Ogata, Stripe states with spatially oscillating d-wave superconductivity in the twodimensional t-t'-J model, Physical review letters 88 , 117001 (2002).
- [23] H.-C. Jiang and T. P. Devereaux, Superconductivity in the doped Hubbard model and its interplay with nextnearest hopping t', Science 365 , 1424 (2019).
- [24] P. Corboz, T. M. Rice, and M. Troyer, Competing states in the t-J model: uniform D-wave state versus stripe state, Physical review letters 113 , 046402 (2014).
- [25] B.-X. Zheng and G. K.-L. Chan, Ground-state phase diagram of the square lattice Hubbard model from density matrix embedding theory, Physical review. B, Condensed matter 93 , 035126 (2016).
- [26] Z. Chen, Y. Wang, S. N. Rebec, T. Jia, M. Hashimoto, D. Lu, B. Moritz, R. G. Moore, T. P. Devereaux, and Z.X. Shen, Anomalously strong near-neighbor attraction in doped 1D cuprate chains, Science 373 , 1235 (2021).
- [27] Y. Wang, Z. Chen, T. Shi, B. Moritz, Z.-X. Shen, and T. P. Devereaux, Phonon-Mediated Long-Range Attractive Interaction in One-Dimensional Cuprates, Physical review letters 127 , 197003 (2021).
- [28] T. Tang, B. Moritz, C. Peng, Z.-X. Shen, and T. P. Devereaux, Traces of electron-phonon coupling in onedimensional cuprates, Nature communications 14 , 3129 (2023).
- [29] C. Peng, Y. Wang, J. Wen, Y. S. Lee, T. P. Devereaux, and H.-C. Jiang, Enhanced superconductivity by near-neighbor attraction in the doped extended Hubbard model, Physical review. B, Condensed matter 107 , L201102 (2023).
- [30] D.-W. Qu, B.-B. Chen, H.-C. Jiang, Y. Wang, and W. Li, Spin-triplet pairing induced by near-neighbor attraction in the extended Hubbard model for cuprate chain, Communications Physics 5 , 1 (2022).
- [31] W.-C. Chen, Y. Wang, and C.-C. Chen, Superconducting phases of the square-lattice extended Hubbard model, Physical review. B, Condensed matter 108 , 064514 (2023).
- [32] J. Gukelberger, E. Kozik, L. Pollet, N. Prokof'ev, M. Sigrist, B. Svistunov, and M. Troyer, p-Wave superfluidity by spin-nematic Fermi surface deformation, Physical review letters 113 , 195301 (2014).
- [33] C. Kallin, Chiral p-wave order in Sr2RuO4, Reports on Progress in Physics 75 , 042501 (2012).
- [34] P. A. Lee and X.-G. Wen, Spin-triplet p -wave pairing in a three-orbital model for iron pnictide superconductors, Physical review. B, Condensed matter 78 , 144517 (2008).
- [35] Y. Guo and H. Tajima, Competition between pairing and tripling in one-dimensional fermions with coexistent s -and p -wave interactions, Physical review. B, Condensed matter 107 , 024511 (2023).
- [36] A. P. Mackenzie and Y. Maeno, The superconductivity of sr 2 ruo 4 and the physics of spin-triplet pairing, Reviews of Modern Physics 75 , 657 (2003).
- [37] K. Kinjo, M. Manago, S. Kitagawa, Z. Q. Mao, S. Yonezawa, Y. Maeno, and K. Ishida, Superconducting spin smecticity evidencing the Fulde-Ferrell-LarkinOvchinnikov state in Sr 2 RuO , Science 4 376 , 397 (2022).
- [38] J. A. Sauls, The order parameter for the superconducting phases of UPt , Advances in Physics 3 43 , 113 (1994).

- [39] L. Jiao, S. Howard, S. Ran, Z. Wang, J. O. Rodriguez, M. Sigrist, Z. Wang, N. P. Butch, and V. Madhavan, Chiral superconductivity in heavy-fermion metal UTe , 2 Nature 579 , 523 (2020).
- [40] Y. Jia, P. Wang, C.-L. Chiu, Z. Song, G. Yu, B. J¨ck, a S. Lei, S. Klemenz, F. A. Cevallos, M. Onyszczak, N. Fishchenko, X. Liu, G. Farahi, F. Xie, Y. Xu, K. Watanabe, T. Taniguchi, B. A. Bernevig, R. J. Cava, L. M. Schoop, A. Yazdani, and S. Wu, Evidence for a monolayer excitonic insulator, Nature physics 18 , 87 (2021).
- [41] K. Matano, M. Kriener, K. Segawa, Y. Ando, and G.-q. Zheng, Spin-rotation symmetry breaking in the superconducting state of cu x bi 2 se 3 , Nature Physics 12 , 852 (2016).
- [42] Y. Liu, T. Wei, G. He, Y. Zhang, Z. Wang, and J. Wang, Pair density wave state in a monolayer high-Tc ironbased superconductor, Nature 618 , 934 (2023).
- [43] S. D. Sarma, M. Freedman, and C. Nayak, Majorana zero modes and topological quantum computation, npj Quantum Information 1 , 1 (2015).
- [44] V. Cr´ epel and L. Fu, Spin-triplet superconductivity from excitonic effect in doped insulators, Proceedings of the National Academy of Sciences of the United States of America 119 , e2117735119 (2022).
- [45] S. Zhang, J. Carlson, and J. E. Gubernatis, Constrained path quantum Monte Carlo method for fermion ground states, Physical review letters 74 , 3652 (1995).
- [46] S. Zhang, J. Carlson, and J. E. Gubernatis, Constrained path Monte Carlo method for fermion ground states, Physical review. B, Condensed matter 55 , 7464 (1997).
- [47] See Supplemental Material, which includes Ref. [57-59], for details on the method and additional simulation results.
- [48] O. K. Andersen, A. I. Liechtenstein, O. Jepsen, and F. Paulsen, LDA energy bands, low-energy hamiltonians, t ′ , t ′′ , t ⊥ ( k ), and J ⊥ , The Journal of physics and chemistry of solids 56 , 1573 (1995).
- [49] M. Hirayama, Y. Yamaji, T. Misawa, and M. Imada, Ab initio effective Hamiltonians for cuprate superconductors,

Physical review. B, Condensed matter 98 , 134501 (2018).

- [50] T. Dahm, V. Hinkov, S. V. Borisenko, A. A. Kordyuk, V. B. Zabolotnyy, J. Fink, B. B¨ uchner, D. J. Scalapino, W. Hanke, and B. Keimer, Strength of the spin-fluctuation-mediated pairing interaction in a high-temperature superconductor, Nature physics 5 , 217 (2009).
- [51] X. Dong, L. Del Re, A. Toschi, and E. Gull, Mechanism of superconductivity in the Hubbard model at intermediate interaction strength, Proceedings of the National Academy of Sciences of the United States of America 119 , e2205048119 (2022).
- [52] Y. Otsuka, K. Seki, S. Sorella, and S. Yunoki, Dirac electrons in the square-lattice Hubbard model with a d -wave pairing field: The chiral Heisenberg universality class revisited, Physical review. B 102 , 235105 (2020).
- [53] A. Luther and I. Peschel, Calculation of critical exponents in two dimensions from quantum field theory in one dimension, Physical review b 12 , 3908 (1975).
- [54] H. J. Schulz, Correlation exponents and the metalinsulator transition in the one-dimensional Hubbard model, Physical review letters 64 , 2831 (1990).
- [55] T. Das, Nodeless superconducting gap induced by odd parity pair density wave in underdoped cuprates, Annals of physics 420 , 168251 (2020).
- [56] Y.-M. Lu, T. Xiang, and D.-H. Lee, Underdoped superconducting cuprates as topological superconductors, Nature physics 10 , 634 (2014).
- [57] K. Liu, S. Yang, W. Li, T. Ying, J. Yang, X. Sun, and X. Li, The pairing symmetries in the two-dimensional Hubbard model, Physics letters. A 392 , 127153 (2021).
- [58] A. T. Rømer, T. A. Maier, A. Kreisel, I. Eremin, P. Hirschfeld, and B. M. Andersen, Pairing in the twodimensional hubbard model from weak to strong coupling, Physical Review Research 2 , 013108 (2020).
- [59] E. Y. Loh, J. E. Gubernatis, R. T. Scalettar, S. R. White, D. J. Scalapino, and R. L. Sugar, Sign problem in the numerical simulation of many-electron systems, Physical review. B, Condensed matter 41 , 9301 (1990).

## Supplemental Material for 'Dominant p -wave pairing induced by near-neighbor attraction in the square-lattice extended Hubbard model'

Zhangkai Cao, 1, ∗ Jianyu Li, 1, 2, ∗ Jiahao Su, 1, 2 Tao Ying, 3, † and Ho-Kin Tang 1, 2, ‡

1 School of Science, Harbin Institute of Technology, Shenzhen, 518055, China 2 Shenzhen Key Laboratory of Advanced Functional Carbon Materials Research and Comprehensive Application, Shenzhen 518055, China. 3 School of Physics, Harbin Institute of Technology, Harbin 150001, China (Dated: January 20, 2025)

In this Supplemental Material, we present additional computational results. In Sec. S1, we have defined various correlation functions. In Sec. S2, we provide the effective d-wave and p-wave pairing correlation function in momentum space with the change of U , V and δ . In Sec. S3, we give a phase diagram of the dominant pairing channel in U -n space when V =0. In Sec. S4, we give a brief introduction to constraint path quantum Monte Carlo (CPQMC).

## S1. Correlation function

We have defined the effective pair momentum distribution function where ζ = s, d x 2 -y 2 , d xy , p . The effective real-space correlation of the on-site s -wave pairing operator C eff s -pair ( i, j ) = 〈 ∆ ( )∆ ( ) † s i s j 〉 -G ↑ i,j G ↓ i,j , where ∆ ( ) = † s i c † i ↑ c † i ↓ is the on-site s -wave pairing operator and G ↑ i,j G ↓ i,j is the uncorrelated pairing structure factor. G σ i,j = 〈 c iσ c † jσ 〉 is the Green's function. The effective real-space correlation of the NN sites d x 2 -y 2 -wave pairing operator C eff d x 2 -y 2 -pair ( i, j ) = ∑ δ ζ ,δ ′ ζ 〈 ∆ † d x 2 -y 2 ( i, i + δ ζ )∆ d x 2 -y 2 ( j, j + δ ′ ζ ) 〉-G ↑ i,j G ↓ i + δ ζ ,j + δ ′ ζ ), where ∆ † d x 2 -y 2 ( ) = i c † i ↑ ( c † i + x ↓ -c † i + y ↓ + c † i - ↓ x -c † i - ↓ y ) and δ ( ) ′ ζ are NN sites in the d x 2 -y 2 -wave pairing. The effective realspace correlation of the NNN sites d xy -wave pairing operator C eff d xy -pair ( i, j ) = ∑ δ ζ ,δ ′ ζ 〈 ∆ † d xy ( i, i + δ ζ )∆ d xy ( j, j + δ ′ ζ ) 〉 -G ↑ i,j G ↓ i + δ ζ ,j + δ ′ ζ ), where ∆ † d xy ( ) = i c † i ↑ ( c † i + + x y ↓ -c † i -x + y ↓ + c † i - - ↓ x y -c † i + x - ↓ y ) and δ ( ) ′ ζ are the NNN sites in the d xy -wave pairing. The effective real-space correlation of the NN sites p -wave pairing operator C eff p -pair ( i, j ) = ∑ δ ζ ,δ ′ ζ 〈 ∆ ( † p i, i + δ ζ )∆ ( p j, j + δ ′ ζ ) 〉 -G σ i,j G σ i + δ ζ ,j + δ ′ ζ ), in which δ ( ) ′ ζ are the NN sites in the p -wave pairing. ∆ ( ) = † p i c † iσ ( c † i + lσ -c † i -lσ ), with spin σ = ↑ ↓ , representing the ↑↑ and ↓↓ pairing, and l = x, y corresponds to the symmetries of p x and p y , respectively. Considering the structure of square lattice, these different pairing symmetries have the following form factor, as shown in Fig. S1. The effective on-site s -wave pairing exists only in the weak U and large V regions. The values of the effective s -wave are small and do not dominate the system, so it was not included it in the main text.

$$\ x l \, & \text{the effective pair momentum distribution function} \\ N _ { \zeta - p a i r } ( k ) = ( 1 / N ) \sum _ { i, j } \exp [ i k ( r _ { i } - r _ { j } ) ] C _ { \zeta - p a i r } ^ { \text{eff} } ( i, j ), & & ( S 1 ) \\ \ x \, & \ d \quad n \, \text{ The effective real sense correlation of the entire wave rainbow or color } \, C _ { \text{eff} } \quad ( i \, \ i ) =$$

FIG. S1. (Color online) Phase of the d x 2 -y 2 , d xy , p ( p x or p y ) pairing symmetry.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-eb9ddc1d-4ba2-4dc0-b292-d5420630a8c9.jpg)

∗ These authors contributed equally.

‡ denghaojian@hit.edu.cn

† taoying86@hit.edu.cn

## S2. The effective d -wave and p -wave pairing in momentum space evolve with the change of U , V and δ

This section we analyze the characteristic of effective d -wave and p -wave pairing correlation function in momentum space evolving with the change of U , V and δ . In Fig. S2(a), we show the effective d -wave pairing correlation function N eff d x 2 -y 2 -pair ( k ) at δ = 0.125 with U = 1 2 3 4 and , , , V = 0 0 . , -0 2 . , -0 4 . , -0 6 . , -0 8 . , -1 0. . We can compare the phase diagram (Fig. 1(b)) to the dominant region of the d x 2 -y 2 -wave pairing, d -wave pairing correlation function exhibits a singular nonzero condensation structure, and its peak occurs near nonzero momentum point Q = (2 π/ , 3 π ), which demonstrate the presence of the d -wave PDW phase. However, in the dominant region of the p -wave pairing, d -wave pairing gradually evolves into a Q = (0 , 0) condensation, exhibiting a significantly weaker strength compared to the p -wave pairing. In Fig. S2(b), we show the effective p -wave pairing correlation function N eff p -pair ( k ) in momentum space at δ = 0.125 with U = 1 2 3 4 and , , , V = 0 0 . , -0 2 . , -0 4 . , -0 6 . , -0 8 . , -1 0. . It clearly that the p -wave pairing strength increases with the increase of V and suppressed by the increase of U . Noteworthy, when the maximum value in momentum space is less than 0, it indicates the absence of p -wave pairing. In the dominant region of the p -wave pairing, the p -wave pairing always condenses at point Q = (0 , 0).

We also present how the effective d -wave and p -wave pairings evolve with doping. As shown in Fig. S3(a)-(c), with doping increases, we can see the condensation point becomes slightly dispersed and the intensity of the d -wave pairing decreases, which may be related to the suppression of the d -wave by doping. At the same time, the p -wave pairing has always been zero momentum condensation, and its intensity has remained basically unchanged when evolve with the change of δ . This indicates that as δ increases, the expansion of the p -wave pairing dominant region is mainly due to the decrease of d -wave pairing.

In the dominant region of d -wave pairing, the position of its condensation point will not change significantly, but with the change of U , V and δ , the sharpness of its nonzero condensation point will change slightly. In the dominant region of the p -wave pairing, it always condenses at point Q = (0 , 0), which supporting the formation of Cooper-pair triplets with zero center-of-mass momentum and p -wave symmetry.

FIG. S2. (Color online) (a) CPQMC simulation result of effective d -wave pair momentum distribution function N eff d x 2 -y 2 -pair ( k ) at δ = 0.125 with U = 1 2 3 4 , , , and V = 0 0 . , -0 2 . , -0 4 . , -0 6 . , -0 8 . , -1 0 . on 12 × 12 lattice. (b) CPQMC simulation result of effective p -wave pair momentum distribution function N eff p -pair ( k ) at δ = 0.125 with U = 1 2 3 4 , , , and V = 0 0 . , -0 2 . , -0 4 . , -0 6 . , -0 8 . , -1 0 on 12 . × 12 lattice.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-62679007-fd8b-4854-9e0d-f0537f046380.jpg)

FIG. S3. (Color online) (a)-(c) CPQMC simulation result of effective d -wave pair momentum distribution function N eff d x 2 -y 2 -pair ( k ) at δ = 0.014, 0.125, 0.208 with U = 3 and V = -0 3 on 12 . × 12 lattice. (d)-(f) CPQMC simulation result of effective p -wave pair momentum distribution function N eff p -pair ( k ) at δ = 0.014, 0.125, 0.208 with U = 2 and V = -0 9 . on 12 × 12 lattice.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-1a6d7913-5f71-46b7-a927-f77e9bf9a443.jpg)

FIG. S4. (Color online) The zero-temperature schematic phase diagram of the dominant pairing channel in U -n space when V =0.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01113v2-figure-bce51901-e1cc-4e55-bb6c-9b9b180d44c3.jpg)

## S3. U -n phase diagram

In this section, we mainly explore the dominant pairing channels in different interaction strength and doping regimes when V = 0. At zero temperature, we summarize the phase diagram of the square-lattice extended Hubbard model in Fig. S4, where includes the d x 2 -y 2 -wave and d xy -wave pairing phases in U -n space when V = 0. At n = 0 681 . -1, the d x 2 -y 2 -wave pairing phase dominates the system when there is no attraction V , while for small electron densities n &lt; 0 681, . the dominant pairing channel is the d xy -wave. This is a schematic phase diagram drawn based on 12 × 12 lattice, which will become more accurate as the system size increase. Our results agree qualitatively with the dynamic cluster approximation (DCA) and determinant quantum Monte Carlo (DQMC) studies [S1, S2], but there are differences in details. The phase boundary between the d x 2 -y 2 -wave and d xy -wave pairing is n ∼ 0 6 in their study, . but in our study, this boundary is n ∼ 0 7. . This is mainly because our CPQMC method is study the ground state of the system at zero temperature. Furthermore, we did not observe the presence of p -wave pairing at low electron density when V = 0.

## S4. Constraint path quantum Monte Carlo

The CPQMC method is a quantum Monte Carlo method with a constraint path approximation [S3, S4]. CPQMC method prevents the infamous sign problem [S5] encountered in the DQMC method when dealing with systems that are far from half filling. In CPQMC, the ground state wave function ψ ( n ) is represented by a finite ensemble of Slater determinants, i.e., ∣ ∣ ψ ( n ) 〉 ∝ ∑ ∣ ∣ k ∣ φ ( n ) k 〉 , where k is the index of the Slater determinants, and n is the number of iteration. The overall normalization factor of the wave function has been omitted here. The propagation of the Slater determinants dominates the computational time, as follows

$$\begin{array} { c } \colon \text{computational} \, \text{in}, \, \text{as} \, \text{nows} \\ & \left | \phi _ { k } ^ { ( n + 1 ) } \right \rangle \leftarrow \int d \vec { x } P ( \vec { x } ) B ( \vec { x } ) \left | \phi _ { k } ^ { ( n ) } \right \rangle. & ( S 2 ) \\ \\ \text{ld configuration, that we select according to the probability distribution function } P ( \vec { x } ). \\ \text{he matrix multiplication of the non-actor} \, R ( \vec { x } ) \, \text{and} \, \vec { a } ^ { ( n ) } \text{ after a series of equilibrium} \end{array}$$

where /vector x is the auxiliary-field configuration, that we select according to the probability distribution function P /vector ( x ). The propagation includes the matrix multiplication of the propagator B /vector ( x ) and φ ( n ) k . After a series of equilibrium steps, the walkers are the Monte Carlo samples of the ground state wave function φ (0) and ground-state properties can be measured.

The random walk formulation suffers from the sign problem because of the fundamental symmetry existing between the fermion ground state | ψ 0 〉 and its negative -| ψ 0 〉 . In more general cases, walkers can cross N in their propagation by e -∆ τH whose bounding surface N is defined by 〈 ψ 0 | φ 〉 = 0 and is in general unknown . Once a random walker reaches N , it will make no further contribution to the representation of the ground state since

Paths that result from such a walker have equal probability of being in either half of the Slater-determinant space. Computed analytically, they would cancel, but without any knowledge of N , they continue to be sampled in the random walk and become Monte Carlo noise.

$$\text{$me no name conunation to one repositions on one group above one} \\ \langle \psi _ { 0 } \left | \phi \right \rangle = 0 \Rightarrow \left \langle \psi _ { 0 } \left | e ^ { - \tau H } \right | \phi \right \rangle = 0 \quad \text{for any $\tau$}. \\ \text{rom such a walker have equal probability of being in either half of the Slater-determinant space}. \text{ly, they would cancel, but without any knowledge of $\mathcal{N}$, they continue to be sampled in the}$$

The decay of the signal-to-noise ratio, i.e., the decay of the average sign of 〈 ψ T | φ 〉 , occurs at an exponential rate with imaginary time. To eliminate the decay of the signal-to-noise ratio, we impose the constrained path approximation. It requires that each random walker at each step has a positive overlap with the trial wave function | ψ T 〉 :

$$\left \langle \psi _ { T } \, | \, \phi _ { k } ^ { ( n ) } \right \rangle & > 0. \\ \intertext { to the ground-state wave function. } | \psi _ { K } \right \rangle & = \Sigma _ { A } | \phi \rangle. \text{ in which all Slater deter-}$$

The details of hyper-parameters in CPQMC calculations are given as follows. The number of walkers is 1000, the number of blocks for relaxation is 10, the number of blocks for growth estimate is 3, the number of blocks after relaxation is 10, the number of steps in each block is 640, the Trotter step size is 0.01, the growth-control energy estimate is -50. To ensure the stability of the data, we intentionally increased the number of steps in each block is 640, and reduced the Trotter step size is 0.01.

This yields an approximate solution to the ground-state wave function, | ψ c 0 〉 = Σ φ | φ 〉 , in which all Slater determinants | φ 〉 satisfy Eq. S3. From Eq. S4, it follows that this approximation becomes exact for an exact trial wave function | ψ T 〉 = | ψ 0 〉 .

[S1] K. Liu, S. Yang, W. Li, T. Ying, J. Yang, X. Sun, and X. Li, The pairing symmetries in the two-dimensional Hubbard model, Physics letters. A 392 , 127153 (2021).

[S2] A. T. Rømer, T. A. Maier, A. Kreisel, I. Eremin, P. Hirschfeld, and B. M. Andersen, Pairing in the two-dimensional hubbard model from weak to strong coupling, Physical Review Research 2 , 013108 (2020).

[S3] S. Zhang, J. Carlson, and J. E. Gubernatis, Constrained path quantum Monte Carlo method for fermion ground states, Physical review letters 74 , 3652 (1995).

[S4] S. Zhang, J. Carlson, and J. E. Gubernatis, Constrained path Monte Carlo method for fermion ground states, Physical review. B, Condensed matter 55 , 7464 (1997).

[S5] E. Y. Loh, J. E. Gubernatis, R. T. Scalettar, S. R. White, D. J. Scalapino, and R. L. Sugar, Sign problem in the numerical simulation of many-electron systems, Physical review. B, Condensed matter 41 , 9301 (1990).
Zhangkai Cao Jianyu Li Jiahao Su Tao Ying Ho-Kin Tang
p-wave superconductivity induced by nearest-neighbor attraction in the square-lattice extended Hubbard model
The two-dimensional (2D) Hubbard model is widely believed to contain the key
ingredients of high-temperature superconductivity in cuprate materials. Here,
we report a constrained path quantum Monte Carlo (CPQMC) study of the
square-lattice extended Hubbard model with on-site Coulomb repulsion U and
nearest-neighbor (NN) electron attraction V. Upon doping $\delta$= 0.125, we
find that the NN electron attraction V can notably drive an exotic spin-triplet
(p-wave) superconducting (SC) phase, and enhance the p-wave SC correlations
with the increase of V. But in the intermediate coupling regime, the
$d_{x^2-y^2}$-wave (d-wave) does not significantly increase with the increase
of V, indicating that the d-wave is not affected by V in strongly correlated
system. Besides the pairing phase, a spin density wave (SDW) only exists near
the half-filling in the particle-hole channel, and doping disrupts the
formation of SDW order. Especially, the NN electron attraction V has no
significant effect on SDW, reflecting the consistent relationship between
d-wave SC and spin correlation. Moreover, as doping increases, the dominant
region of p-wave also expands, further suppressing the presence of d-wave,
which may help explain the disappearance of d-wave SC in overdoped cuprate
superconductors. We also find the d-wave exhibits a singular nonzero (near
point ($2{\pi}/3, {\pi}$)) condensation structure in momentum space, resulting
in a different staggered behavior in the x and y direction with distance in
real space. On the contrary, the p-wave condensed at zero momentum, and the
p-wave correlation exhibits exponential decay in real space. Our work suggests
the p-wave SC region can be induced and further broadened by the NN electron
attraction V, offering a feasible mechanism to realize p-wave superconductivity
in realistic cuprate materials.