Mechanism of High-Temperature Superconductivity in Correlated-Electron Systems

It is very important to elucidate the mechanism of superconductivity for achieving room temperature superconductivity. This paper is a short review article on the mechanism of high-temperature superconductivity. In the first half of this paper, we give a brief review on mechanisms of superconductivity in many-electron systems. We believe that high-temperature superconductivity may occur in a system with interaction of large-energy scale. Empirically, this is true for superconductors that have been found so far. In the second half of this paper, we discuss cuprate high-temperature superconductors. We argue that superconductivity of high temperature cuprates is induced by the strong on-site Coulomb interaction, that is, the origin of high-temperature superconductivity is the strong electron correlation. We show the results on the ground state of electronic models for high temperature cuprates on the basis of the optimization variational Monte Carlo method. A high-temperature superconducting phase will exist in the strongly correlated region.


Introduction
][3] For this purpose, it is important to clarify the ground state and phase diagram of electronic models with strong correlation because high temperature cuprates are strongly correlated electron systems.
Most superconductors induced by the electron-phonon interaction have s-wave pairing symmetry.][6] The critical temperature T c of most of electron-phonon superconductors is very low except for exceptional compounds.Many unconventional superconductors that cannot be understood by the BCS theory have been discovered.They are, for example, heavy fermion superconductors, organic superconductors and cuprate superconductors for which the pairing mechanism is different from the electron-phonon interaction.In particular, cuprate superconductors exhibit relatively high T c and have become of great interest.A common feature in both electronphonon systems and correlated electron systems is that critical temperature may have a strong correlation with the energy scale of the interaction that induces electron pairing.
This paper has two parts.In the first part, we give a review on mechanisms of superconductivity in the electron-phonon system and in the correlated electron system.In the second part, we mainly discuss the mechanism of high-temperature cuprates.
[9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24] It is certainly a very difficult task to elucidate the phase diagram of the d-p model.Simplified models are also used to investigate the mechanism of superconductivity, for example the two-dimensional (2D) single-band Hubbard model  and ladder model [51][52][53][54][55][56] . The Hubard model was introduced to understand the metal-insulator transition 57) and was employed to understand various magnetic phenomena.][70][71] It was also expected that the Hubbard model can account for high temperature superconductivity.72) A variational Monte Carlo method is used to examine the ground state properties of strongly correlated electron systems, where we calculate the expectation values exactly using a numerical method.[28][29][30][31][36][37][38][39][40][41] We introduced the wave function of exp(−S )-type in the study of superconductivity in the Hubbard model.[73][74][75] This wave function is very excellent in the sense that the energy expectation value is lower than that of any other wave functions. 50) The papr is organized as follows.In Section 2.1, we discuss the phonon mechanism of superconductivity.In Section 2.2, we discuss the electron mechanism of superconductivity.Section 2.3 is devoted to a discussion on superconductivity in correlated electron systems. In Section 3. we show the model for high temperature cuprates.We present the optimization variational Monte Carlo method (OVMC) in Section 3.2. We ow the results on superconductivity based on the OVMC in Section 3.3.We discuss the stability of antiferromagnetic state in Section 3.4.We show the phase diagram when the hole doping rate is changed in Section 3.5.We give a summary in Section 4.

Possibility of High-T c Superconductivity
][6] There are many superconductors with s-wave pairing symmetry and most of them are due to the electron-phonon interaction.The BCS theory was successful to explain physical properties of these superconductors.
In the strong-coupling theory based on the Green function formulation, 76,77) the critical temperature T c was estimated as, 78) where λ is the electron-phonon coupling constant, θ D = ω D /k B is the Debye temperature and µ * is the renormalized Coulomb parameter defined by for µ = U/ǫ F where U is the strength of the Coulomb interaction and ǫ F is the Fermi energy.µ * is the phenomenological parameter being approximately 0.1.The electron-phonon coupling constant λ is expressed as where α(ω) is the averaged electron-phonon coupling over the Fermi surface and F(ω) indicates the product of the spectral function of phonon and the density of states.This is approximately written as where ρ(ǫ F ) is the density of states at the Fermi surface and M is the mass of an atom.McMillan predicted that T c would have a limit being of the order of 30 K from the analysis for this formula. 78)he McMillan formula was modified by replacing θ D /1.45 by logarithmic Debye frequency ω ln where 79) It was predicted that high critical temperature would be possible for large λ since T c ∝ √ λ for λ ≫ 1.If ω ln is large, λ is also large, and the crystal is stable, high T c would be realized.It was predicted that high T c would be realized in hydrogen solid with high Debye temperature. 80)[83] It is important to consider multi-band superconductors in the search for high temperature superconductors.In fact, MgB 2 and iron based superconductors are multi-band superconductors.An important role of Lifshitz transition in iron based superconductors and MgB 2 multi-band superconductors has been predicted. 84)An interesting point is that the possibility of high-T c superconductivity in materials where tuning the chemical potential shows a quasi-1D Fermi surface topology as in organics and hydrides. 85)A layered superconductor such as cuprate superconductor can be regarded as a multiband superconductor due to interlayer couplings.15][116][117][118][119] When we have multiple order parameters, there appear multiple Nambu-Goldstone bosons and Higgs bosons.104,[120][121][122][123][124][125][126][127][128][129] This will result in significant excitation modes that are unique in multiband superconductors.
It is important to include in a theoretical picture the presence of multiple electronic components with anomalous normal state properties in the charge and spin sector, e.g., the well known Fermi arcs and charge pseudogap phenomenology.The "shape resonance" scenario of multigap BCS-BEC crossover has been proposed. 130,131) 134]

Electron Correlation and Superconductivity
We discuss the electron correlation due to the Coulomb interaction between electrons.The on-site Coulomb interaction is important in the study of the metal insulator transition and magnetic properties of materials.The Hubbard model is written as 25) where t i j indicates the transfer integral and the second term denotes the Coulomb interaction with the strength U. t i j are chosen as follows.t i j = −t when i and j are nearest-neighbor pairs i j and t i j = −t ′ when i and j are next-nearest neighbor pairs.In the following, N is the number of lattices, and N e denotes the number of electrons.When two electrons spin up and down at the same site, the energy becomes higher by U where U denotes the on-site Coulomb energy.In the case of half-filling, the Mott transition occurs when U(> 0) is as large as the bandwidth and the ground state is an insulator.The effective Hamiltonian is derived in the limit of large U/t, [135][136][137] based on the canonical transformation H eff = e iS He −iS . I the limit U/t → ∞, the double occupancy is not allowed.The effective Hamiltonian is written as We write the Hamiltonian as H = H0 + H 1 where Here, we defined is the electron operator without double occupancy.We choose S to satisfy i[S , H0 ] + H 1 = 0, so that H eff reads in the subspace of no double occupancy, When we consider only the nearest-neighbor transfer t i j = −t, the effective Hamiltonian reads where j + µ and j + µ ′ denote the nearest-neighbor sites in the µ and µ ′ directions, respectively.The second term being proportional to t 2 /U contains the nearest-neighbor exchange interaction and also three-site terms when µ µ ′ .The threesite terms are of the same order as the exchange interaction.When we neglect the three-site terms, the effective Hamiltonian reduces to the t-J model given by where J = 4t 2 /U and ñi = ñi↑ + ñi↓ with ñiσ = a † iσ a iσ .High-temperature cuprates and heavy fermion systems are typical correlated electron systems and many superconductors have been reported.Most of superconductors in these systems have nodes in the superconducting gap, namely, the Cooper pair is anisotropic.This indicates that superconductivity is unconventional and does not conform to the conventional BCS theory.The mechanism of superconductivity is certainly nonphonon mechanism.We show several characteristic properties of cuprate high-temperature superconductors: 1.The Cooper pair has d-wave symmetry.
2. The superconducting phase exists near the antiferromagnetic phase and parent materials are a Mott insulator.

The CuO 2 plane is commonly contained and the on-site
Coulomb repulsive interaction works between d electrons.4. The size of Cooper pair is very small being of order of 2Å. 5.The CuO 2 plane is high anisotropic and there is a weak Josephson coupling between two layers.][140][141] A plausible non-phonon mechanism is due to the Coulomb interaction on the same atom.Because the energy scale of the Coulomb interaction is very large, which is of the order of eV, we can expect superconductivity with high critical temperature T c .The critical temperature of heavy fermion materials is, however, very low, although superconductivity occurs due to strong Coulomb interaction between f electrons.This is because the effective mass of f electrons is very large in heavy fermion systems owing to the large self-energy correction.The effective mass enhancement of heavy fermion materials becomes as large as 100-1000, which means that the effective cutoff becomes very small.As a result, the characteristic energy scale is reduced considerably and the critical temperature T c becomes very low begin of the order of 1 K.In heavy fermion systems, the characteristic energy scale is given by the Kondo temperature T K .The ratio of the effective mass m * to the band mass m 0 is approximately given as m * /m 0 ≃ D/T K for the bandwidth D and T K .Thus, the effective bandwidth for heavy fermions is given by the Kondo temperature T K ≃ D/(m * /m 0 ).Empirically, T c is lowered as the effective mass increases.This is expressed as follows: where t denotes the transfer integral proportional to the bandwidth.The estimated values of the transfer t, the ratio m * /m 0 and T c for several compounds are shown in Table I.For cuprates, the transfer t is estimated as t ∼ 0.51 eV.The bandwidth for iron pnictides is about five times smaller than that for cuprates.A list of typical superconductors in correlated electron systems is shown in Table II.

Superconductivity in Strongly Correlated Electron Systems
The possibility of superconductivity in strongly correlated electron systems has been discussed intensively.][178] There were, however, the results by quantum Monte Carlo methods, which did not support the existence of high-temperature superconductivity in the two-dimensional Hubbard model. 32,33,45) I quantum Monte Carlo calculations, the strength of the Coulomb interaction U is not large enough because the range of accessible U is very restricted.It is now certain that there is a superconducting phase in the strongly correlated region. 50)The simplest wave function of superconducting state with strong electron correlation is the Gutzwiller-projected BCS wave function: where u k and v k are BCS parameters and P G is the Gutzwiller operator to control the on-site electron correlation.P G is written as where g is a variational parameter in the range of 0 ≤ g ≤ 1.
The ratio of u k and v k is given as where ξ k denotes the electron dispersion relation measured from the Fermi energy and ∆ k is the gap function.We use the following form for the gap function in the two-dimensional case: ∆ is a constant and is treated as a variational parameter.The wave function ψ BCS −G is just the wave function that Anderson proposed as a wave function of the resonate-valence-bond (RVB) state. 72)t has been shown that the ground-state energy has a minimum at finite ∆ for the BCS-Gutzwiller wave function with d-wave symmetry in the two-dimensional Hubbard model by using the variational Monte Carlo method. 37)The superconducting condensation energy E cond per site in the limit of large system size was estimated as 37,38) E cond /N ≃ 0.2 meV, (20)   where the transfer integral t is set at 0.5 eV.The similar result was obtained for the three-band d-p model. 18)Thus, the condensation energy per atom is of the order of 10 −4 eV.The superconducting condensation energy E cond for cuprate high-temperature superconductors was evaluated by using the result of specific heat measurement for YBCO as 0.17-0.26meV per Cu atom. 37,179) he estimation of E cond from the data of critical magnetic field gives the similar result. 180)The obtained results by theoretical calculations and experimental measurements are very close each other.This agreement is very remarkable.Thus, this value indicates the characteristic energy for cuprate high-temperature superconductors.This result may support that the superconductivity in cuprate high temperature superconductors is caused by the strong electron correlation and the 2D Hubbard model includes essential ingredients.

Part II. Mechanism of Superconductivity in Cuprates
We discuss the mechanism of superconductivity in this part.We show numerical results obtained by using the optimized wave functions.

Model for High-T c Cuprates
The Hamiltonian of the d-p model for high-T c cuprates is  1).This value is determined from the sign of the transfer integral between next nearest-neighbor d orbitals.U d indicates the strength of the on-site Coulomb repulsion between d holes and U p is that between p holes.
][183][184][185] We show an example: U d = 10.5, U p = 4.0 and U d p = 1.2 in eV. 182)Here, U d p is the nearest-neighbor Coulomb interaction between holes on adjacent Cu and O orbitals and is small compared to U d .U d p is neglected in this paper.We write ∆ d p = ǫ p − ǫ d .The number of sites is denoted as N, and the energy is measured in units of t d p .

Optimization Variational Monte Carlo Method 3.2.1 Off-Diagonal Wave Function
The Gutzwiller wave function is where ψ 0 is a one-particle state.Our purpose is to improve the Gutzwiller function.We multiply the Gutzwiller function by an exponential-type operator.The wave function is given as as 50,73,[186][187][188][189][190] where K denotes the kinetic part of the Hamiltonian.λ is a newly introduced real variational parameter. 41,73,187,191) The are other methods to improve the Gutzwiller function. 43,192) he following Jastrow operator is used, 43) where d j is the operator for the doubly-occupied site given as d j = n j↑ n j↓ and e j is that for the empty site given by e j = (1 − n j↑ )(1 − n j↓ ).η is the variational parameter in the range of 0 ≤ η ≤ 1.The wave function is In this paper, we use the wave function of exponential type in Equation ( 23) because the energy is further lowered when we use this wave function. 50)The wave function for the dp model is formulated similarly.An initial state ψ 0 contains many variational parameters (t d p , tpp , t′ d , and ǫp − ǫd ): We use tdp = t d p as the energy unit.We consider the following wave function that is improved from the Gutzwiller wave function: 50,73,[186][187][188][189][190] The expectation values are evaluated by using the auxiliary field method. 73,191) he kinetic part K also contains the band parameters t pp , t ′ d and ǫ p − ǫ d as variational parameters: We take tpp = tpp , t′ d = t′ d and ǫp − ǫd = ǫp − ǫd , for simplicity.Thus, we have g, tpp , t′ d , ǫp − ǫd = ǫp − ǫd , and λ as variational parameters.The expectation values for this type of wave function are calculated on the basis of the variational Monte Carlo method.One can evaluate the expectation value correctly within statistical errors.

Antiferromagnetic Wave Function
The AF one-particle state ψ AF is formulated by the eigenfunction of the AF trial Hamiltonian: where ∆ AF is the AF order parameter and (x i , y i ) represents the coordinates of the site i.With ψ AF , the wave function is given as

Superconducting Wave Function
We start from the BCS wave function with coefficients u k and v k satisfying where P N e indicates the operator to extract the state with N e electrons.The exponential-BCS wave function is given by In this wave function, we perform the electron-hole transformation for down-spin electrons: and not for up-spin electrons: c k = c k↑ .The electron pair operator c † k↑ c † −k↓ denotes the hybridization operator c † k d k in this formulation.

Correlated Superconductivity
We first discuss the superconducting (SC) state in the two-dimensional Hubbard model.In the optimization Monte Carlo method, the SC state becomes indeed stable when the Coulomb interaction U is large to be of the order of the bandwidth.We show the ground-state energy as a function of the superconducting order parameter ∆ in Figure 2 (left).The simple Gutzwiller-projected BCS wave function predicted the possibility of superconductivity in the Hubbard model, and the improved wave function also shows a stability of the SC state.
We show the SC and antiferromagnetic (AF) order parameters as a function of U in Figure 3.The AF order parameter has a peak when U/t ∼ 10, which is of the order of the bandwidth, and the SC one also has a peak at U c that is greater than the bandwidth.This indicates that there is the possibility of high-temperature superconductivity in the strongly correlated region.
The AF correlation is maximized at U ∼ U c and decreases when U is larger than U c .We show schematic pictures in Figure 4, where the SC condensation energy as a function of U is shown in the left panel, and the AF and SC gap functions are shown in the right panel.There is a crossover from weakly correlated region to the strongly correlated region.The superconducting state is most favorable when the AF correlation is gradually suppressed in the strongly correlated region.Thus,  (Lower) The ground-state energy as a function of the electron number where ∆ is fixed for each line. 50)gh temperature superconductivity is highly promising in the strongly correlated region where U is as large as the bandwidth D or larger than D.

Stability of Antiferromagnetic State 3.4.1 Hubbard Model
Let us examine the stability of AF state.There are two parameters U and t ′ , and there is the AF region in the parameter space.High temperature superconductivity is expected near the boundary between the AF phase and the paramagnetic phase.We show the AF condensation energy ∆E AF as a function of 1 − n e in Figure 5a for t ′ = 0 and Figure 5b for t ′ = −0.2t.The AF region becomes larger as |t ′ | increases.When t ′ = −0.2t, the AF region extends up to about 20% doping.From the competition between superconductivity and AF order, t ′ = 0 is most favorable for superconductivity.

Three-Band d-p Model
In general, in the three-band d-p model, the AF correlation is very strong and the AF state is more stable than in the single-band Hubbard model.This is because d electrons are localized and easily form magnetic order. 13)To investigate the possibility of high temperature superconductivity in the d-p model, it is necessary to reveal regions with weak AF order.There are many parameters in the d-p model to control the strength of the AF correlation.Among them, the Coulomb repulsion between d electrons U d , the level difference ∆ d p = ǫ p − ǫ d , and the hole density x are important.The AF region is shown in Figure 6 where U d and ǫ p − ǫ d are varied, and the hole density is fixed at 0.1875.The AF region increases when the hole density decreases.We expect that high temperature superconductivity will occur near the boundary between AFM and PM regions.This boundary exists in the region when ∆ d p is small.High temperature superconductivity is likely occur when ∆ d p is small.There is a "onsite attractive region" when ∆ d p is large where two d electrons prefer to occupy the same site.In this region, a charge-density wave or an s-wave superconducting state will be realized.We proposed to introduce the transfer integral t ′ d to control the strength of AF correlation. 74)We show the AF region at half-filling in the t pp − t ′ d plane in Figure 7.As −t ′ d increases, there is a phase transition from the AF insulator to the paramagnetic insulator (PMI).We expect that t ′ d and t pp will play an important role to suppress AF correlation when holes are doped in the d-p model.

Phase Diagram for the Hubbard Model
We discuss the phase diagram when carrier holes are doped in the CuO 2 plane.We evaluate the energy lowering when we include the order parameter ∆.We define where E(∆) takes a minimum at ∆ = ∆ min .We show ∆E as a function of the hole doping rate x in Figure 8 where we put U/t = 18 and t ′ = 0.This phase diagram contains several interesting features.There are three phases: antiferromagnetic insulator (AFI), coexistent state (AFSC) and superconducting phase (SC).When the hole doping rate x is large, e.g., x > 0.09, the pure d-wave stat is stable.There is the possibility of high (and room) temperature superconductivity in this phase.In the underdoped region, approximately 0.06 < x < x dS C with 0.08 < x dS C < 0.09, there is the coexistent state of antiferromagnetism and superconductivity.This is the mixed phase of AF and SC.x dS C could not be deter- mined precisely.There is the possibility that both the AFSC and SC states are found for x dS C < x < 0.09, but the SC solution will have lower energy.There is the AFSC-SC transition at x = x dS C .The AFI state exists near half-filling for about x < 0.06, where doped holes form clusters and localize.
The existence of AFI phase is closely related to the phase separation 140,141) when the hole density is very small.In the phase-separated phase, the doped holes are localized and cannot be conductive.The existence of AFI phase is determined by the quantity (36)   where E(N e ) is the ground-state energy with N e electrons.δ 2 E(N e ) is approximately the second derivative of the energy E(N e ) and is proportional to the charge susceptibility.When δ 2 E(N e ) is negative, the phase separation occurs.As shown in Figure 8, the phase separation occurs for x < 0.06.Concerning the phase separation, the parameter t ′ is important because the phase separation region decreases as −t ′ increases.Thus, the AFI phase will decrease as −t ′ increases.The phase separation disappears for t ′ = −0.2t.

Summary
We have discussed the possibility of high temperature superconductivity in many-electron systems.The critical temperature T c may increase as the characteristic energy of the interaction increases.Empirically, T c is proportional to the inverse of the effective mass of electrons.T c is low when the effective mass is very heavy.A candidate of high (room) temperature superconductivity may be in materials with strong electron correlation and with small effective mass enhancement.From this view point, the repulsive Coulomb interaction can be a candidate of the origin of high temperature superconductivity.
We have shown phase diagrams for the 2D Hubbard model and the three-band d-p model.The diagram in Figure 8 exhibits the characteristic property of cuprate superconductors.This supports that the origin of high temperature superconductivity is the strong correlation between electrons.That is, the mechanism of high-T c superconductivity is the electronpair formation due to the strong on-site repulsive Coulomb interaction.The competition between antiferromagnetism and superconductivity is important in realizing high temperature superconductivity.High-T c superconductivity is expected in the region near the boundary between AF phase and paramagnetic phase.In the phase diagram for the Hubbard model, the SC phase exists near the AF phase, and AF order and superconductivity coexist where the doping rate is approximately 0.05∼0.06< x < x dS C and 0.08 < x dS C < 0.09.We expect that this coexistence may be related to anomalous metallic behavior in the underdoped region.The AF phase near half-filling is insulating, which is approximately for x < 0.06.There is the pure d-wave phase for x > x dS C .
In the d-p model, the AF region exists in the multi- The energy unit is given by t d p .AFM and PM denote the antiferromagnetic metal and paramagnetic metal, respectively.There is a "negative-U" region when the level difference is large where two d electrons prefer to occupy the same site.The ground state may be a charge-density wave state or an s-wave superconducting state.This is not clear yet.dimensional parameter space.The AF-PM boundary is a multi-dimensional region in this space.Since we expect that superconductivity occurs near the boundary, high temperature superconductivity is more likely to occur in the d-p model.There is the AF-PM boundary when the level difference ∆ d p is small.Thus, T c of high temperature cuprates will be high when ∆ d p is small.This tendency is consistent with experimental T c of cuprates.
We give a comment on the crossover between weakly correlated region and strongly correlated region.We expect that this crossover is universal in the sense that similar phenomena occur in nature.There may be a universal class.It will include the Kondo effect, [194][195][196] QCD, 197) BCS-BEC crossover, 198) sine-Gordon model, [199][200][201][202] and Gross-Neveu model.203)

2 Fig. 1 .
Fig. 1. (Upper) The transfer integral t ′ d in the CuO 2 plane.t d p and t pp are conventionally defined.(Lower) Fermi surface of the d-p model with the hole density 0.13. 74)We put t pp = 0.4t d p and ǫ p − ǫ d = 2t d p for t ′ d = 0, −0.2t d p and −0.3t d p .

Fig. 2 .
Fig. 2. (Upper) The ground-state energy as a function of the superconducting gap ∆ for the optimized wave function ψ λ for the Hubbard model on a 10 × 10 lattice with U/t = 18 and t ′ = 0.The electron density is n e = 0.88.(Lower)The ground-state energy as a function of the electron number where ∆ is fixed for each line.50)

Fig. 3 .
Fig.3.AF and SC order parameters as a function of U/t when N e = 88 for the 2D Hubbard model on a 10 × 10 lattice.The periodic boundary conditions are periodic in one direction and antiperiodic in the other direction.50)AF(G) indicates the result obtained for the simple Gutzwiller function.

Fig. 4 .
Fig. 4. (Upper) A schematic picture of the superconducting condensation energy as function of U for the 2D Hubbard model.(Lower) A schematic picture of the gap function of AF and AC states as a function of U for the 2D Hubbard model.

Fig. 6 .
Fig.6.Antiferromagnetic and paramagnetic regions in the plane of U d and ∆ d p = ǫ p − ǫ d for the d-p model.We put t pp = 0.4 and t ′ d = 0.There are 76 holes on a 8 × 8 lattice with 192 atoms in total.The energy unit is given by t d p .AFM and PM denote the antiferromagnetic metal and paramagnetic metal, respectively.There is a "negative-U" region when the level difference is large where two d electrons prefer to occupy the same site.The ground state may be a charge-density wave state or an s-wave superconducting state.This is not clear yet.

Fig. 8 .
Fig. 8.The condensation energy per site for the two-dimensional Hubbard model as a function of the hole doping rate x.The calculation was carried out on a 10 × 10 lattice.AFI indicates the antiferromagnetic insulating state and SC denotes the d-wave SC phase.There is the coexistent state indicated as AF-SC between these states.Parameters are t ′ = 0 and U/t = 18.

Table I .
The transfer integral t, effective mass m * and critical temperature T c in correlated electron systems, where m 0 denotes the band mass.The orders of these quantities are shown in the table.For heavy fermions, t/(m * /m 0 ) corresponds to the Kondo temperature T K .For Hydrides, the Debye frequency ω ln is shown.

Table II .
Superconducting materials.
Since we use the hole picture in this paper, d iσ and d † iσ represent the operators for the d hole.p i± x/2σ and p † i± x/2σ denote the operators for the p holes at the site R i± x/2 , and in a similar way p i±ŷ/2σ and p † i±ŷ/2σ are defined.n p i+ x/2σ and n p i+ŷ/2σ are the number operators of p holes at R i+ x/2 and R i+ŷ/2 , respectively.t d p is the transfer integral between adjacent Cu and O orbitals and t pp is that between nearest p orbitals.t ′ d indicates that between d orbitals where i j denotes a next nearest-neighbor pair of copper sites.ǫ i j takes the values ±1 (see Figure 2σ p i+ x/2σ + p † i+ŷ/2σ p i+ŷ/2σ )