Arxiv:1907.02313V1

Home , Cuprate superconductor

arXiv:1907.02313v1 [cond-mat.str-el] 4 Jul 2019 a temperature cal urt uecnutr xii eaieyhigh relatively exhibit superconductors cuprate urt uecnutr o hc h arn ehns is mechanism pairing di the which for superconductors cuprate a example, superconductors for organic are, superconductors, They fermion discovered. heavy by been have understood theory be BCS the cannot un- Many that compounds. superconductors exceptional conventional for except low very is tors at emil ics h ehns fhigh-temperature of second mechanism the the In cuprates. discuss system. mainly electron we correlated part, the in electron–phonon and the in system superconductivity of mechanisms on pairing. electron induces that energ interaction the the with of correlation scale crit strong that a have is may systems temperature electron cal correlated electron– both and in systems feature phonon common A interest. great of become aldtedpmdlo h he-adHbadmodel. Hubbard three-band the or model d-p the called neo nooeeu ttsrpre o urts uhas such cuprates, for reported appear states the understand inhomogeneous to of possible ance is it model, Hubbard the of transition metal–insulator the derstand model ladder a ensuiditnieyfrmr hn3 years. 30 than more it for indeed intensively and studied superconductivity, been has temperature high of anism Introduction 1. w-iesoa 2)snl-adHbadmodel the inves- Hubbard example to single-band for used (2D) superconductivity, also two-dimensional of are mechanism models the Simplified tigate model. d-p the of scranyavr di very a certainly is h Bardeen–Cooper–Schrie the oudrtn aiu antcphenomena. magnetic various understand to hsclpoete of properties physical eato have teraction correlated strongly are correlation cuprates systems. and electron strong state temperature with ground high models the because electronic clarify of to diagram important phase is it purpose, this odne atr21,4 57 4, 2019, Matter Condensed ff h oe o CuO for model The review a give we part, first the In parts. two has paper This ti hlegn eerhsbett lrf h mech- the clarify to subject research challenging a is It otsprodcosidcdb h lcrnpoo in- electron–phonon the by induced superconductors Most rn rmteeeto–hnnitrcin nparticular In interaction. electron–phonon the from erent ot al ehd ihtmeauesprodcigph tempe superconducting high high-temperature for A models method. electronic Carlo of Monte state ground super the high-temperature on of sults origin the is, tem ha that high second interaction, of the superconductivity In that far. argue We so perconductors. found been have that superconductors eprtr uecnutvt a cu nasse ihin with system a supercond in of occur mechanisms may on mechanism superconductivity review the temperature brief on a article give review we paper, short this a is paper This tivity. ainlIsiueo dacdIdsra cec n Tech and Science Industrial Advanced of Institute National ti eyipratt lcdt h ehns fsupercond of mechanism the elucidate to important very is It 51–56) ehns fHg-eprtr uecnutvt nCorre in Superconductivity High-Temperature of Mechanism s wv arn ymty ecnunderstand can We symmetry. pairing -wave T c h ubr oe a nrdcdt un- to introduced was model Hubbard The . fms feeto–hnnsuperconduc- electron–phonon of most of ffi utts oeuiaetepaediagram phase the elucidate to task cult 2 s ln ncpaesprodcosis superconductors cuprate in plane wv uecnutvt ae on based superconductivity -wave ff r(C)theory. (BCS) er 57) n a employed was and 8 59) 58, 4–6) T ntebasis the On c h criti- The 25–50) n have and aah Yanagisawa Takashi 1–3) 7–24) and For Systems nd It i- odciiyi h togeeto orlto.W hwth show We correlation. electron strong the is conductivity y - , 1 ooy111Ueoo skb,Iaai3586,Japan 305-8568, Ibaraki Tsukuba, Umezono, 1-1-1 nology eauecpae sidcdb h togo-ieCoulomb on-site strong the by induced is cuprates perature auecpae ntebsso h piiainvariationa optimization the of basis the on cuprates rature eato flreeeg cl.Eprcly hsi true is this Empirically, scale. large-energy of teraction es hr ecluaeteepcainvle xcl usi method. sy exactly numerical values electron a expectation the correlated calculate strongly we where of tems, properties state ground utr with ductors formulation, superco these of properties ductors. physical was theory explain BCS to The successful interaction. electron–phonon the to due twsas xetdta h ubr oe a con for account can model superconductivity. Hubbard temperature the high that expected also was It stripes as, ino exp( of tion ue oidc trcieitrcinbteneetosa is electrons symmetry between pairing the interaction attractive induce to sumed . osblt fHigh-T of Possibility Systems Many-Electron in 2.1 Superconductivity I. Part 2. give We 3.5. 4. Section Section in in changed summary diagram is a phase rate the doping show hole We the antifer- 3.4. when of Section stability in the state on discuss We romagnetic based 3.3. superconductivity Section on in OVMC results the the Sec- show We in (OVMC) 3.2. op- method tion the Carlo present We Monte variational cuprates. temperature timization sh we high 3.1, for Section model In the supercondu systems. on electron discussion correlated a in tivity to devoted superconduc- Sec- is of 2.3 mechanism In Section electron tivity. superconductivity. the discuss of we 2.2, mechanism tion phonon the cuss h es htteeeg xetto au slwrta tha than lower functions. wave is other value any expectation of energy the that sense the model. Hubbard where fhg-eprtr uecnutvt.I h rthl of half first the In superconductivity. high-temperature of s ileiti h togycreae region. correlated strongly the in exist will ase fo hsppr edsuscpaehg-eprtr su- high-temperature cuprate discuss we paper, this of lf ciiyi ayeeto ytm.W eiv hthigh- that believe We systems. many-electron in uctivity aitoa ot al ehdi sdt xmn the examine to used is method Carlo Monte variational A ntesrn-opigter ae nteGenfunction Green the on based theory strong-coupling the In nteBSter,teeeto–hnnitrcini as- is interaction electron–phonon the theory, BCS the In h ae sognzda olw.I eto .,w dis- we 2.1, Section In follows. as organized is paper The 78) ciiyfraheigro eprtr superconduc- temperature room achieving for uctivity 60–67) λ steeeto–hnnculn constant, coupling electron–phonon the is − T n hcebadlk est aestates. wave density checkerboard-like and 6 77) 76, c s S wv arn ymtyadms fte are them of most and symmetry pairing -wave = -yei h td fsprodciiyi the in superconductivity of study the in )-type 73–75) 1 h rtcltemperature critical the θ . 45 D 83,36–41) 28–31, hswv ucini eyecletin excellent very is function wave This exp s -wave. c

− Superconductivity lated-Electron 50) λ eitoue h aefunc- wave the introduced We − 4–6) 1 µ . 04(1 ∗ hr r aysupercon- many are There (1 72) + + 0 λ . 62 ) T c λ ) a estimated was ! re- e , Yanagisawa T for l θ 68–71) D ow (1) ng nd n- c- s- = t Condens. Matter 2019, 4, 57 T Yanagisawa

~ωD/kB is the Debye temperature and µ∗ is the renormalized e.g., the well known Fermi arcs and charge pseudogap phe- Coulomb parameter defined by nomenology. The “shape resonance” scenario of multigap µ BCS-BEC crossover has been proposed.130, 131) The study µ∗ = , (2) of the electronic structure of the cuprates superconductors 1 + µ ln(ǫF /ωD) Bi Sr CaCuO and La CuO doped by mobile oxygen in- = 2 2 8+y 2 4+y for µ U/ǫF where U is the strength of the Coulomb inter- terstitials using local probes has shown a scenario made of ǫ µ action and F is the Fermi energy. ∗ is the phenomenological two electronic components: a strongly correlated Fermi liquid parameter being approximately 0.1. which coexists with stripes made of anisotropic polarons con- The electron–phonon coupling constant λ is expressed as densed into a generalized Wigner charge density wave.132–134) α(ω)2F(ω) = ∞ λ 2 dω , (3) 2.2 Electron Correlation and Superconductivity Z0 ω We discuss the electron correlation due to the Coulomb in- where α(ω) is the averaged electron–phonon coupling over teraction between electrons. The on-site Coulomb interaction the Fermi surface and F(ω) indicates the product of the spec- is important in the study of the metal insulator transition and tral function of phonon and the density of states. This is ap- magnetic properties of materials. The Hubbard model is writ- proximately written as ten as25) 2 ρ(ǫF ) I λ h i , (4) H = ti jc† c jσ + U ni ni , (6) ≃ Mω2 iσ ↑ ↓ D Xi jσ Xi where ρ(ǫF ) is the density of states at the Fermi surface and where ti j indicates the transfer integral and the second term M is the mass of an atom. McMillan predicted that T would c denotes the Coulomb interaction with the strength U. ti j are have a limit being of the order of 30 K from the analysis for chosen as follows. ti j = t when i and j are nearest-neighbor this formula.78) − pairs i j and ti j = t′ when i and j are next-nearest neighbor The McMillan formula was modified by replacing θD/1.45 pairs.h Ini the following,− N is the number of lattices, and N 79) e by logarithmic Debye frequency ωln where denotes the number of electrons. When two electrons spin up and down at the same site, the 2 ∞ 2 ln ω ωln = exp dωα(ω) F(ω) . (5) energy becomes higher by U where U denotes the on-site λ Z0 ω ! Coulomb energy. In the case of half-filling, the Mott transi- It was predicted that high critical temperature would be tion occurs when U(> 0) is as large as the bandwidth and possible for large λ since Tc √λ for λ 1. If ω is large, ∝ ≫ ln the ground state is an insulator. The effective Hamiltonian is λ is also large, and the crystal is stable, high Tc would be derived in the limit of large U/t,135–137) based on the canon- realized. It was predicted that high Tc would be realized in H = eiS He iS U/t 80) ical transformation eff − . In the limit , hydrogen solid with high Debye temperature. In fact, high the double occupancy is not allowed. The effective Hamilto-→ ∞ temperature superconductors with Tc above 200 K were dis- nian is written as covered under extremely high pressure (160 200 GPa) in hy- 2 81–83)∼ i drogen compounds such as H3S and LaH10. Heff = H + i[S, H] + [S, [S, H]] + . (7) It is important to consider multi-band superconductors in 2 ··· the search for high temperature superconductors. In fact, We write the Hamiltonian as H = H˜0 + H1 where MgB2 and iron based superconductors are multi-band super- H˜0 = ti j(a† a jσ + d† d jσ) + U ni ni , (8) conductors. An important role of Lifshitz transition in iron iσ iσ ↑ ↓ Xi jσ Xi based superconductors and MgB2 multi-band superconductors has been predicted.84) An interesting point is that the H1 = ti j(ai†σd jσ + d†jσaiσ). (9) possibility of high-Tc superconductivity in materials where Xi jσ tuning the chemical potential shows a quasi-1D Fermi sur- 85) Here, we defined aiσ = ciσ(1 ni, σ) and diσ = ciσni, σ. face topology as in organics and hydrides. A layered su- − − − aiσ = ciσ(1 ni. σ) is the electron operator without double perconductor such as cuprate superconductor can be regarded − occupancy. We− choose S to satisfy i[S, H˜ ] + H = 0, so that as a multiband superconductor due to interlayer couplings. 0 1 H reads in the subspace of no double occupancy, A multi-band superconductivity has been investigated as a eff generalization of the BCS theory since early works on the i i2 H ff = ti ja† a jσ + [S, H ] + [S, [S, H ]] + . (10) two-band superconductivity.86–89) There will appear many in- e iσ 2 1 3 1 ··· Xi jσ teresting properties in superconductors with multiple gaps such as time-reversal symmetry breaking,90–103) the exis- When we consider only the nearest-neighbor transfer ti j = 104–109) ff 110–114) t, theeffective Hamiltonian reads tence of massless modes, unusual isotope e ect − and fractional-flux quantum vortices.115–119) When we have t2 Heff = t a† a jσ a† a† a j a j+µ multiple order parameters, there appear multiple Nambu– − iσ − U j+µ j ↓ ′↑ 104, 120–129) Xi j σ Xjµµ ↑ ↓ Goldstone bosons and Higgs bosons. This will re- h i ′ sult in significant excitation modes that are unique in multi- + a† a† a a + a† a† a a + a† a† a a , j j+µ j+µ′ j j+µ j j+µ′ j j j+µ j j+µ band superconductors. ↑ ↓ ↓ ↑ ↑ ↓ ↓ ↑ ↑ ↓ ↓ ↑ It is important to include in a theoretical picture the (11) presence of multiple electronic components with anoma- where j + µ and j + µ′ denote the nearest-neighbor sites in lous normal state properties in the charge and spin sector, the µ and µ′ directions, respectively. The second term being

2 Condens. Matter 2019, 4, 57 T Yanagisawa

Table I. The transfer integral t, effective mass m and critical tempera- 2 ∗ proportional to t /U contains the nearest-neighbor exchange ture Tc in correlated electron systems, where m0 denotes the band mass. interaction and also three-site terms when µ , µ′. The three- The orders of these quantities are shown in the table. For heavy fermions, site terms are of the same order as the exchange interaction. t/(m∗/m0) corresponds to the Kondo temperature TK . For Hydrides, the De- When we neglect the three-site terms, the effective Hamilto- bye frequency ωln is shown. nian reduces to the t-J model given by ∗ ∗ t or ωln m /m0 t/(m /m0) Tc 1 Cuprates 5000 K 5 1000 100 K t 0.51eV142) ∼ 143) Heff = (ai†σa jσ + h.c.) + J Si S j n˜in˜ j , (12) Fe pnictides 1000 K 2 500 50 K t 0.1eV − · − 4 ! ∼ ∼ Xi j σ Xi j Heavy fer 10000 K 100 1000 10 100 1 10 K 144–146) h i h i ∼ ∼ ∼ Organics 200 500 K 2 5 100 10 K 147) 2 ∼ ∼ 148) where J = 4t /U andn ˜i = n˜i + n˜i withn ˜iσ = ai†σaiσ. Hydrides 1000 K 1 1000 100 K ωln High-temperature cuprates↑ and↓ heavy fermion systems are ∼ typical correlated electron systems and many superconductors have been reported. Most of superconductorsin these systems have nodes in the superconducting gap, namely, the Cooper Table II. Superconducting materials. pair is anisotropic. This indicates that superconductivityis unconventional and does not conform to the conventional BCS Materials Tc Pair Symmetry Refs 149,150) theory. The mechanism of superconductivity is certainly non- CeCu2Si2 0.6 K s or d 151) UPt3 0.52 K p or f phonon mechanism. We show several characteristic properties 152) UBe13 0.86 K p of cuprate high-temperature superconductors: 153–155) URu2Si2 1.2 K 156) 1. The Cooper pair has d-wave symmetry. CeRu2 6.2 K s 157–159) UPd2Al3 2K d 2. The superconducting phase exists near the antiferromag- 158,160) UNi2Al3 1K p ? netic phase and parent materials are a Mott insulator. 161,162) CeCoIn5 2.3 K d 163) 3. The CuO2 plane is commonly contained and the on-site CeRhIn5 2.1 K d Coulomb repulsive interaction works between d elec- (16.3 kbar) 164) trons. CeRh2Si2 0.35 K (9 kbar) 4. The size of Cooper pair is very small being of order of 165) UGe2 0.8 K p ? 2Å. (13.5 kbar) 5. The CuO2 plane is high anisotropic and there is a weak URhGe 0.25 K p ? 166) 167) Josephson coupling between two layers. Sr2RuO4 1.5 K p or f 168) PrOs4Sb12 1.85 K line nodes ? The small size of Cooper pair also supports the non-phonon 169) 138–141) NaxCoO2 y H2O 5K p ? mechanism of cuprate superconductivity. A plausible − · 170) Ba1 xKxBiO3 30 K s − 171) non-phonon mechanism is due to the Coulomb interaction MgB2 39 K s on the same atom. Because the energy scale of the Coulomb La2 xSrxCuO4 36 K d − interaction is very large, which is of the order of eV, we YBa2Cu3O6+x 90 K d Tl2Ba2Can 1CunO2n+4 125 K d can expect superconductivity with high critical temperature − HgBa2Can 1CunO2n+2+δ 135 K d T . The critical temperature of heavy fermion materials is, − 172) c LaO1 xFxFeAs 26K − 173) however, very low, although superconductivity occurs due to NdFeAsO1 y 54 K − 81) strong Coulomb interaction between f electrons. This is be- H3S 203 K s 82,83,174,175) cause the effective mass of f electrons is very large in heavy LaH10 260 K s 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 charac- 2.3 Superconductivity in Strongly Correlated Electron Sys- teristic energy scale is reduced considerably and the critical tems temperature T becomes very low begin of the order of 1 K. c The possibility of superconductivity in strongly correlated In heavy fermion systems, the characteristic energy scale is electron systems has been discussed intensively. The per- given by the Kondo temperature T . The ratio of the effec- K turbative calculations such as the fluctuation-exchange ap- tive mass m to the band mass m is approximately given as ∗ 0 proximation (FLEX) have been performed to investigate the m /m D/T for the bandwidth D and T . Thus, the ef- ∗ 0 K K superconducting ground state.176–178) There were, however, fective bandwidth≃ for heavy fermions is given by the Kondo the results by quantum Monte Carlo methods, which did not temperature T D/(m /m ). Empirically, T is lowered as K ∗ 0 c support the existence of high-temperature superconductivity the effective mass≃ increases. This is expressed as follows: in the two-dimensional Hubbard model.32, 33, 45) In quantum kBTc 0.1t/(m∗/m0), (13) Monte Carlo calculations, the strength of the Coulomb inter- ≃ action U is not largeenough becausethe rangeof accessible U where t denotes the transfer integral proportional to the band- is veryrestricted. It is now certain that there is a superconduct- width. The estimated values of the transfer t, the ratio m /m ∗ 0 ing phase in the strongly correlated region.50) The simplest and T for several compounds are shown in Table I. For c wave function of superconducting state with strong electron cuprates, the transfer t is estimated as t 0.51 eV. The band- correlation is the Gutzwiller-projected BCS wave function: width for iron pnictides is about five times∼ smaller than that for cuprates. A list of typical superconductors in correlated ψBCS G = PG (uk + vkc† c† ) 0 , (14) − k k | i electron systems is shown in Table II. Yk ↑ − ↓

3 Condens. Matter 2019, 4, 57 T Yanagisawa

where uk and vk are BCS parameters and PG is the Gutzwiller + tdp [d† (pi+xˆ/2σ + pi+yˆ/2σ pi xˆ/2σ pi yˆ/2σ) + h.c.] iσ − − − − operator to control the on-site electron correlation. PG is writ- Xiσ ten as + tpp [pi†+yˆ/2σ pi+xˆ/2σ pi†+yˆ/2σ pi xˆ/2σ pi† yˆ/2σ pi+xˆ/2σ − − − − PG = (1 (1 g)n j n j ), (15) Xiσ − − ↑ ↓ Yj +pi† yˆ/2σ pi xˆ/2σ + h.c.] where g is a variational parameter in the range of 0 g 1. − − ≤ ≤ The ratio of uk and vk is given as + td′ ǫi j(d† d jσ + h.c.) + Ud d† di d† di iσ i ↑ i ↓ Xi j σ Xi ↑ ↓ v ∆ hh ii k = k , (16) u 2 2 1/2 + p p + p p k ξk + (ξk + ∆k) Up (ni+xˆ/2 ni+xˆ/2 ni+yˆ/2 ni+yˆ/2 ). (21) Xi ↑ ↓ ↑ ↓ where ξk denotes the electron dispersion relation measured from the Fermi energy and ∆k is the gap function. We use the Since we use the hole picture in this paper, diσ and di†σ rep- following form for the gap function in the two-dimensional resent the operators for the d hole. pi xˆ/2σ and pi† xˆ/2σ denote ± ± case: the operators for the p holes at the site Ri xˆ/2, and in a similar p ± p ∆ = ∆(cos k cos k ) d wave, (17) way pi yˆ/2σ and pi† yˆ/2σ are defined. ni+xˆ/2σ and ni+yˆ/2σ are the k x y ± ± − − number operators of p holes at Ri+xˆ/2 and Ri+yˆ/2, respectively. ∆k = ∆(cos kx + cos ky) anisotropic s wave,(18) t is the transfer integral between adjacent Cu and O orbitals − − dp and tpp is that between nearest p orbitals. t′ indicates that be- ∆k = ∆ isotropic s wave. (19) d − tween d orbitals where i j denotes a next nearest-neighbor hh ii ∆ is a constant and is treated as a variational parameter. The pair of copper sites. ǫi j takes the values 1 (see Figure 1). ± wave function ψBCS G is just the wave function that Ander- This value is determined from the sign of the transfer inte- − son proposed as a wave function of the resonate-valence-bond gral between next nearest-neighbor d orbitals. Ud indicates (RVB) state.72) the strength of the on-site Coulomb repulsion between d holes It has been shown that the ground-state energy has a min- and Up is that between p holes. imum at finite ∆ for the BCS-Gutzwiller wave function with The values of band parameters were evaluated by several 181–185) d-wave symmetry in the two-dimensional Hubbard model by works. We show an example: Ud = 10.5, Up = 4.0 37) 182) using the variational Monte Carlo method. The supercon- and Udp = 1.2 in eV. Here, Udp is the nearest-neighbor ducting condensation energy Econd per site in the limit of large Coulomb interaction between holes on adjacent Cu and O or- 37, 38) system size was estimated as bitals and is small compared to Ud. Udp is neglected in this paper. We write ∆dp = ǫp ǫd. The number of sites is denoted Econd/N 0.2 meV, (20) − ≃ as N, and the energy is measured in units of tdp. 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 con- 4 densation energy per atom is of the order of 10− eV. The superconducting condensation energy Econd for cuprate high-temperature superconductors was evaluated by using the result of specific heat measurement for YBCO as 37, 179) 0.17–0.26 meV per Cu atom. The estimation of Econd 3.2 Optimization Variational Monte Carlo Method from the data of critical magnetic field gives the similar 3.2.1 Off-Diagonal Wave Function result.180) The obtained results by theoretical calculations and The Gutzwiller wave function is experimental measurements are very close each other. This ψ = P ψ , agreement is very remarkable. Thus, this value indicates the G G 0 (22) characteristic energy for cuprate high-temperature supercon- where ψ0 is a one-particle state. Our purpose is to improve the ductors. This result may support that the superconductivity Gutzwiller function. We multiply the Gutzwiller function by in cuprate high temperature superconductors is caused by an exponential-type operator. The wave function is given as the strong electron correlation and the 2D Hubbard model as50, 73, 186–190) includes essential ingredients. ψλ = exp( λK)ψG, (23) − where K denotes the kinetic part of the Hamiltonian. λ 3. Part II. Mechanism of Superconductivity in Cuprates is a newly introduced real variational parameter.41, 73, 187, 191) We discuss the mechanism of superconductivityin this part. There are other methods to improve the Gutzwiller func- We show numerical results obtained by using the optimized tion.43, 192) The following Jastrow operator is used,43) wave functions.

PJdh = 1 (1 η) d j(1 e j+τ) + e j(1 d j+τ) , T  − − − −  3.1 Model forHigh- c Cuprates Yj  Yτ h i   The Hamiltonian of the d-p model for high-Tc cuprates is  (24) where d j is the operator for the doubly-occupied site given as Hdp = ǫd d† diσ + ǫp (p† pi+xˆ/2σ + p† pi+yˆ/2σ) iσ i+xˆ/2σ i+yˆ/2σ d j = n j n j and e j is that for the empty site given by e j = Xiσ Xiσ ↑ ↓ (1 n j )(1 n j ). η is the variational parameter in the range − ↑ − ↓

4 Condens. Matter 2019, 4, 57 T Yanagisawa

Cu 3.2.2 Antiferromagnetic Wave Function The AF one-particle state ψAF is formulated by the eigen- function of the AF trial Hamiltonian:

’ xi+yi t ! HAF = ti jc† c jσ ∆AF ( 1) σniσ, (29) d iσ − − Xi jσ Xiσ

where ∆AF is the AF order parameter and (xi, yi) represents tdp tpp ! the coordinates of the site i. With ψAF , the wave function is O given as

ψλ,AF = exp( λK)PGψAF. (30) − px 3.2.3 Superconducting Wave Function d x2 −y2 We start from the BCS wave function 1

t ' = 0 t ' = −0.2 ψBCS = (uk + vkck† c† k ) 0 , (31) d d ↑ − ↓ | i 0.8 t ' = −0.3 Yk d 2 2 with coeﬃcients uk and vk satisfying uk + vk = 1. We choose

π | | / 0.6 2 2 y uk/vk = ∆k/(ξk + ξ + ∆ ) for the gap function ∆k and ξk = k k k q ǫk µ. We assume ∆k = ∆SC(cos kx cos ky). The Gutzwiller- 0.4 projected− BCS wave function is −

0.2 ψG BCS = PN PGψBCS , (32) − e

where PNe indicates the operator to extract the state with Ne 0 0 0.2 0.4 0.6 0.8 1 electrons. The exponential-BCS wave function is given by k /π = λK x ψλ e− PGψBCS . (33) In this wave function, we perform the electron–hole trans- Fig. 1. (Upper) The transfer integral td′ in the CuO2 plane. tdp and tpp are conventionally defined. (Lower) Fermi surface of the d-p model with the hole formation for down-spin electrons: 74) density 0.13. We put tpp = 0.4tdp and ǫp ǫd = 2tdp for td′ = 0, 0.2tdp − − d = c† , d† = c ; (34) and 0.3tdp. k k k k − − ↓ − ↓ and not for up-spin electrons: ck = ck . The electron pair oper- ↑ ator ck† c† k denotes the hybridization operator ck†dk in this formulation.↑ − ↓ of 0 η 1. The wave function is ≤ ≤ 3.3 Correlated Superconductivity ψη = PJdhψG. (25) We first discuss the superconducting (SC) state in the In this paper, we use the wave function of exponential type two-dimensional Hubbard model. In the optimization Monte in Equation (23) because the energy is further lowered when Carlo method, the SC state becomes indeed stable when the 50) we use this wave function. The wave function for the d- Coulomb interaction U is large to be of the order of the band- p model is formulated similarly. An initial state ψ0 contains width. We show the ground-state energy as a function of the many variational parameters (t˜dp, t˜pp, t˜′ ,andǫ ˜p ǫ˜d): ∆ d − superconducting order parameter in Figure 2 (left). The simple Gutzwiller-projected BCS wave function predicted ψ0 = ψ0(t˜dp, t˜pp, t˜′ , ǫ˜p ǫ˜d). (26) d − the possibility of superconductivity in the Hubbard model, We use t˜dp = tdp as the energy unit. We consider the follow- and the improved wave function also shows a stability of the ing wave function that is improved from the Gutzwiller wave SC state. 50, 73, 186–190) function: We show the SC and antiferromagnetic (AF) order parameters as a function of U in Figure 3. The AF order parameter ψλ = exp( λK)ψG. (27) − has a peak when U/t 10, which is of the order of the band- ∼ The expectation values are evaluated by using the auxiliary width, and the SC one also has a peak at Uc that is greater 73, 191) field method. The kinetic part K also contains the band than the bandwidth. This indicates that there is the possibility parameters tpp, t and ǫp ǫd as variational parameters: d′ − of high-temperature superconductivity in the strongly correlated region. K = K(tˆpp, tˆd′ , ǫˆp ǫˆd). (28) − The AF correlation is maximized at U Uc and decreases We take tˆpp = t˜pp, tˆ = t˜ andǫ ˆp ǫˆd = ǫ˜p ǫ˜d, for sim- ∼ d′ d′ − − when U is larger than Uc. We show schematic pictures in Fig- plicity. Thus, we have g, t˜pp, t˜ ,ǫ ˆp ǫˆd = ǫ˜p ǫ˜d, and λ as d′ − − ure 4, where the SC condensation energy as a function of U variational parameters. The expectation values for this type of is shown in the left panel, and the AF and SC gap functions wave function are calculated on the basis of the variational are shown in the right panel. There is a crossover from weakly Monte Carlo method. One can evaluate the expectation value correlated region to the strongly correlated region. The super- correctly within statistical errors. conducting state is most favorable when the AF correlation is gradually suppressed in the strongly correlated region. Thus,

5 Condens. Matter 2019, 4, 57 T Yanagisawa

-0.5 0.4 U = 18 t -0.505 n = 0.88 e 0.3 AF(G) ∆ -0.51 E/N 0.2

-0.515 AF

0.1 SC -0.52 0 0.02 0.04 0.06 0.08 0.1 0 0 5 10 15 20 25 -0.5 ∆ U/t ∆ = 0.02 ∆ = 0.04 -0.505 Fig. 3. AF and SC order parameters as a function of U/t when Ne = 88 for ∆ = 0.05 the 2D Hubbard model on a 10 10 lattice. The periodic boundary conditions = 0.06 × ∆ are periodic in one direction and antiperiodic in the other direction.50) AF(G) ∆ = 0.08 -0.51 indicates the result obtained for the simple Gutzwiller function. E/N

-0.515 to control the strength of the AF correlation. Among them, -0.52 the Coulomb repulsion between d electrons Ud, the level difference ∆dp = ǫp ǫd, and the hole density x are important. − The AF region is shown in Figure 6 where Ud and ǫp ǫd -0.525 − 43.8 43.9 44 44.1 44.2 44.3 are varied, and the hole density is ﬁxed at 0.1875. The AF Ne/2 region increases when the hole density decreases. We expect that high temperature superconductivity will occur near the Fig. 2. (Upper) The ground-state energy as a function of the supercon- boundary between AFM and PM regions. This boundary ex- ducting gap ∆ for the optimized wave function ψ for the Hubbard model on λ ists in the region when ∆ is small. High temperature super- a 10 10 lattice with U/t = 18 and t = 0. The electron density is ne = 0.88. dp × ′ (Lower) The ground-state energy as a function of the electron number where conductivity is likely occur when ∆dp is small. There is a “on- 50) ∆ is ﬁxed for each line. site attractive region” when ∆dp 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 td′ to control high temperature superconductivity is highly promising in the the strength of AF correlation.74) We show the AF region at strongly correlated region where U is as large as the band- half-ﬁlling in the tpp td′ plane in Figure 7. As td′ increases, width D or larger than D. there is a phase transition− from the AF insulator− to the para-

magnetic insulator (PMI). We expect that td′ and tpp will play 3.4 Stability of Antiferromagnetic State an important role to suppress AF correlation when holes are 3.4.1 Hubbard Model doped in the d-p model. Let us examine the stability of AF state. There are two parameters U and t′, and there is the AF region in the param- 3.5 Phase Diagram for the Hubbard Model eter space. High temperature superconductivity is expected We discuss the phase diagram when carrier holes are doped near the boundary between the AF phase and the paramag- in the CuO2 plane. We evaluate the energy lowering when we netic phase. We show the AF condensation energy ∆EAF as include the order parameter ∆. We deﬁne a function of 1 ne in Figure 5a for t′ = 0 and Figure 5b − ∆E = E(∆= 0) E(∆min), (35) for t′ = 0.2t. The AF region becomes larger as t′ increases. − − | | When t = 0.2t, the AF region extendsup to about 20% dop- where E(∆) takes a minimum at ∆ = ∆ . We show ∆E ′ − min ing. From the competition between superconductivity and AF as a function of the hole doping rate x in Figure 8 where order, t′ = 0 is most favorable for superconductivity. we put U/t = 18 and t′ = 0. This phase diagram contains several interesting features. There are three phases: antifer- 3.4.2 Three-Band d-p Model romagnetic insulator (AFI), coexistent state (AFSC) and su- In general, in the three-band d-p model, the AF corre- perconducting phase (SC). When the hole doping rate x is lation is very strong and the AF state is more stable than large, e.g., x > 0.09, the pure d-wave stat is stable. There is in the single-band Hubbard model. This is because d elec- the possibility of high (and room) temperature superconduc- trons are localized and easily form magnetic order.13) To in- tivity in this phase. In the underdoped region, approximately vestigate the possibility of high temperature superconductiv- 0.06 < x < xdSC with 0.08 < xdSC < 0.09, there is the coexis- ity in the d-p model, it is necessary to reveal regions with tent state of antiferromagnetism and superconductivity. This weak AF order. There are many parameters in the d-p model is the mixed phase of AF and SC. xdSC could not be deter-

6 Condens. Matter 2019, 4, 57 T Yanagisawa 0.04 (a)

0.03 U = 8

/N U = 12 U = 14 AF

E U = 18 !" ∆ 0.02 ∆

0.01

0 0 0.05 0.1 0.15 0.2 0 5 10 15 20 25 x U/t 0.08 (b) 0.06 U = 12 AF /N U = 14

∆ AF

E U = 18 0.04 ∆

0.02

SC 0 0 0.05 0.1 0.15 0.2 0.25 x 0 ~ D U Fig. 5. AF condensation energy ∆EAF as a function of the hole doping rate 75) x = 1 ne on a 10 10 lattice for: t = 0 (a); and t = 0.2t (b). We put − × ′ ′ − Fig. 4. (Upper) A schematic picture of the superconducting condensation U/t = 12, 14 and 18. 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.

interaction increases. Empirically, Tc is proportional to the inverse of the effective mass of electrons. Tc is low when the mined precisely. There is the possibility that both the AFSC effective mass is very heavy. A candidate of high (room) temperature superconductivity may be in materials with strong and SC states are found for xdSC < x < 0.09, but the SC solu- tion will have lower energy. There is the AFSC-SC transition electron correlation and with small effective mass enhancement. From this view point, the repulsive Coulomb interaction at x = xdSC. The AFI state exists near half-filling for about x < 0.06, where doped holes form clusters and localize. can be a candidate of the origin of high temperature supercon- The existence of AFI phase is closely related to the phase ductivity. separation140, 141) when the hole density is very small. In the We have shown phase diagrams for the 2D Hubbard model phase-separated phase, the doped holes are localized and can- and the three-band d-p model. The diagram in Figure 8 ex- not be conductive. The existence of AFI phase is determined hibits the characteristic property of cuprate superconductors. by the quantity This supports that the origin of high temperature superconductivity is the strong correlation between electrons. That is, 2 2 δ E(Ne) [E(Ne+δNe) 2E(Ne)+E(Ne δNe)]/(δNe) , (36) ≡ − − the mechanism of high-Tc superconductivity is the electron- pair formation due to the strong on-site repulsive Coulomb where E(Ne) is the ground-state energy with Ne electrons. 2 interaction. The competition between antiferromagnetism and δ E(Ne) is approximately the second derivative of the energy superconductivity is important in realizing high temperature E(Ne) and is proportional to the charge susceptibility. When 2 superconductivity. High-Tc superconductivity is expected in δ E(Ne) is negative, the phase separation occurs. As shown in Figure 8, the phase separation occurs for x < 0.06. Concern- the region near the boundary between AF phase and paramagnetic phase. In the phase diagram for the Hubbard model, ing the phase separation, the parameter t′ is important because the SC phase exists near the AF phase, and AF order and the phase separation region decreases as t′ increases. Thus, − superconductivity coexist where the doping rate is approxi- the AFI phase will decrease as t′ increases. The phase sepa- − mately 0.05 0.06 < x < xdSC and0.08 < xdSC < 0.09.We ex- ration disappears for t′ = 0.2t. − pect that this∼ coexistence may be related to anomalous metal- 4. Summary lic behavior in the underdoped region. The AF phase near We have discussed the possibility of high temperature su- half-filling is insulating, which is approximately for x < 0.06. perconductivity in many-electron systems. The critical tem- There is the pure d-wave phase for x > xdSC. In the d-p model, the AF region exists in the multi- perature Tc may increase as the characteristic energy of the

7 Condens. Matter 2019, 4, 57 T Yanagisawa

! ""#$! '"

"" AFM # Ud

!"# "

PM E/N ! ∆ !" AFI On-site Attractive region

"# # SC # ! " # !# AF ∆dp SC

" Fig. 6. Antiferromagnetic and paramagnetic regions in the plane of Ud and " " # "! "!# "" ""# ∆dp = ǫp ǫd for the d-p model. We put tpp = 0.4 and t = 0. There are 76 x − d′ holes on a 8 8 lattice with 192 atoms in total. The energy unit is given × by tdp. AFM and PM denote the antiferromagnetic metal and paramagnetic Fig. 8. The condensation energy per site for the two-dimensional Hubbard metal, respectively. There is a “negative-U” region when the level diﬀerence model as a function of the hole doping rate x. The calculation was carried out is large where two d electrons prefer to occupy the same site. The ground on a 10 10 lattice. AFI indicates the antiferromagnetic insulating state and × state may be a charge-density wave state or an s-wave superconducting state. SC denotes the d-wave SC phase. There is the coexistent state indicated as This is not clear yet. AF-SC between these states. Parameters are t′ = 0 and U/t = 18.

!"&

crossover,198) sine-Gordon model,199–202) and Gross–Neveu model.203) !"%

PMI d Acknowledgments This work was supported by a Grant- -t' in-Aid for Scientiﬁc Research from the Ministry of Education, !"# Culture, Sports, Science and Technology of Japan (Grant No. 17K05559). A part of the computations was supported by the Supercomputer Center of the Institute for Solid State Physics, !"! AFI the University of Tokyo.

!"! 1) Bednorz, J.B.; Müller, K.A. Possible high Tc superconductivity in the !"! !" !"! !"" !"# !"$ Ba-La-Cu-O system. Z. Phys. B Condens. Matter 1986, B64, 189–193. tpp 2) Keimer, B.; Kivelson, S.A.; Norman, M.R.; Uchida, S.; Zaanen, J. From quantum matter to high-temperature superconductivity in cop- Fig. 7. AF and paramagnetic insulator phases for the d-p model on a 6 6 per oxides. Nature 2015, 518, 179. 193) × lattice. Parameters are Ud = 8tdp, Up = 0, ǫp ǫd = tdp. 3) Rybicki, D.; Jurkutat, M.; Reichart, S.; Kapusta, C.; Haase, J. Perspec- − tive on the phase diagram of cuprate high-temperature superconductors. Nat. Commun. 2016, 7, 11413. 4) Cooper, L.N. Bound Electron Paira in a Degenerate Fermi Gas. Phys. Rev. 1956, 104, 1189. dimensional parameter space. The AF-PM boundary is a 5) Bardeen, J.; Cooper, L.N.; Schrieffer, J.R. Microscopic Theory of Su- multi-dimensional region in this space. Since we expect that perconductivity. Phys. Rev. 1957, 106, 162. ff superconductivity occurs near the boundary, high temperature 6) Bardeen, J.; Cooper, L.N.; Schrie er, J.R. Theory of Superconductiv- ity. Phys. Rev. 1957, 108, 1175. superconductivity is more likely to occur in the d-p model. 7) Emery, V.J. Theory of high-Tc superconductivity in oxides. Phys. Rev. There is the AF–PM boundary when the level difference ∆dp Lett. 1987, 58, 2794. is small. Thus, Tc of high temperature cuprates will be high 8) Hirsch, J.E.; Loh, E.Y.; Scalapino, D.J.; Tang, S. Pairing interaction in CuO clusters. Phys. Rev. B 1989, 39, 243. when ∆dp is small. This tendency is consistent with experimental T of cuprates. 9) Scalettar, R.T.; Scalapino, D.J.; Sugar, R.L.; White, S.R. Antiferro- c magnetic, charge-transfer, and pairing correlations in the three-band We give a comment on the crossover between weakly Hubbard model. Phys. Rev. B 1991, 44, 770. correlated region and strongly correlated region. We expect 10) Unger, P.; Fulde, P. Spectral function of holes in the Emergy model. that this crossover is universal in the sense that similar Phys. Rev. B 1993, 48, 16607. phenomena occur in nature. There may be a universal class. 11) Oguri, A.; Asahatam, T.; Maekawa, S. Gutzwiller wave function in the ff 194–196) 197) three-band Hubbard model: A variational Monte Carlo study. Phys. It will include the Kondo e ect, QCD, BCS-BEC Rev. B 1994, 49, 6880.

8 Condens. Matter 2019, 4, 57 T Yanagisawa

12) Koikegami, S.; Yamada, K. Antiferromagnetic and superconducting 42) Bulut, N. dx2 y2 superconductivity and the Hubbard model. Adv. Phys. correlations on the d-p model. J. Phys. Soc. Jpn. 2000, 69, 768. 2002, 51, 1587–1667.− 13) Yanagisawa, T.; Koike, S.; Yamaji, K. Ground state of the three-band 43) Yokoyama, H.; Tanaka, Y.; Ogata, M.; Tsuchiura, H. Crossover of su- Hubbard model. Phys. Rev. B 2001, 64, 184509. perconducting properties and kinetic-energy gain in two-dimensional 14) Koikegami, S.; Yanagisawa, T. Superconducting gap of the two- Hubbard model. J. Phys. Soc. Jpn. 2004, 73, 1119–1122. dimensional d-p model with small Ud. J. Phys. Soc. Jpn. 2001, 70, 44) Yokoyama, H.; Ogata, M.; Tanaka, Y. Mott transitions and d-wave su- 3499–3502. perconductivity in half-filled-band Hubbard model on square lattice 15) Yanagisawa, T.; Koike, S.; Yamaji, K. Lattice distortions, incommen- with geometric frustration. J. Phys. Soc. Jpn. 2006, 75, 114706. surability, and stripes in the electronic model for high-Tc cuprates. 45) Aimi, T.; Imada, M. Does simple two-dimensional Hubbard model ac- Phys. Rev. B 2003, 67, 132408. count for high-Tc superconductivity in copper oxides? J. Phys. Soc. 16) Koikegami, S.; Yanagisawa, T. Superconductivity in Sr2RuO4 medi- Jpn. 2007, 76, 113708. ated by Coulomb scattering. Phys. Rev. B 2003, 67, 134517. 46) Miyazaki, M.; Yanagisawa, T.; Yamaji, K. Diagonal stripe states in the 17) Koikegami, S.; Yanagisawa, T. Superconductivity in multilayer per- light-doping region in the two-dimensional Hubbard model. J. Phys. ovskite. J. Phys. Soc. Jpn. 2006, 75, 034715. Soc. Jpn. 2004, 73, 1643–1646. 18) Yanagisawa, T.; Miyazaki, M.; Yamaji, K. Incommensurate antiferro- 47) Yanagisawa, T. Phase diagram of the t-U2 Hamiltonian of the weak magnetism coexisting with superconductivity in two-dimensional d-p coupling Hubbard model. New J. Phys. 2008, 10, 023014. model. J. Phys. Soc. Jpn. 2009, 78, 031706. 48) Yanagisawa, T. Enhanced pair correlation functions in the two- 19) Weber, C.; Lauchi, A.; Mila, F.; Giamarchi, T. Orbital currents in dimensional Hubbard model. New J. Phys. 2013, 15, 033012. extended Hubbard model of High-Tc cuprate superconductors. Phys. 49) Yokoyama, H.; Ogata, M.; Tanaka, Y.; Kobayashi, K.; Tsuchiura, H. Rev. Lett. 2009, 102, 017005. Crossover between BCS superconductor and doped Mott insulator of 20) Lau, B.; Berciu, M.; Sawatzky, G.A. High spin polaron in lightly d-wave pairing state in two-dimensional Hubbard model. J. Phys. Soc. doped CuO2 planes. Phys. Rev. Lett. 2011, 106, 036401. Jpn. 2013, 82, 014707. 21) Weber, C.; Giamarchi, T.; Varma, C.M. Phase diagram of a three- 50) Yanagisawa, T. Crossover from wealy to strongly correlated regions orbital model for high-Tc cuprate superconductors. Phys. Rev. Lett. in the two-dimensional Hubbard model—Off-diagonal Monte Carlo 2014, 112, 117001. studies of Hubbard model II—. J. Phys. Soc. Jpn. 2016, 85, 114707. 22) Avella, A.; Mancini, F.; Paolo, F.; Plekhanov, E. Emery vs. Hubbard 51) Noack, R.M.; White, S.R.; Scalapino, D.J. The doped two-chain Hub- model for cuprate superconductors: A composite operator method bard model. EPL 1995, 30, 163. study. Eur. Phys. J. 2013, B86, 265. 52) Noack, R.M.; Bulut, N.; Scalapino, D.J.; Zacher, M.J. Enhanced 23) Ebrahimnejad, H.; Sawatzky, G.A.; Berciu, M. Differences between dx2 y2 pairing correlations in the two-leg Hubbard ladder Phys. Rev. the insulating limit quasiparticles of one-band and three-band cuprate B 1997− , 56, 7162. models. J. Phys. Condens. Matter 2016, 28, 105603. 53) Yamaji, K.; Shimoi, Y.; Yanagisawa, T. Superconductivity indications 24) Tamura, S.; Yokoyama, H. Variational study of magnetic ordered state of the two-chain Hubbard model due to the two-band effect. Physica in d-p model. Phys. Procedia 2016, 81, 5–8. C 1994, 235, 2221. 25) Hubbard, J. Electron correlations in narrow energy bands. Proc. R. Soc. 54) Yanagisawa, T.; Shimoi, Y.; Yamaji, K. Superconducting phase of a Lond. 1963, 276, 238–257. two-chain Hubbard model. Phys. Rev. B 1995, 52, R3860. 26) Hubbard, J. Electron correlations in narrow energy bands III. Proc. R. 55) Nakano, T.; Kuroki, K.: Onari, S. Superconductivity due to spin fluc- Soc. Lond. 1964, 281, 401–419. tuations originating from multiple Fermi surfaces in the double chain 27) Gutzwiller, M.C. Effect of correlation on the ferromagnetism of tran- superconductor Pr2Ba4Cu7O15 δ. Phys. Rev. B 2007, 76, 014515. sition metals. Phys. Rev. Lett. 1963, 10, 159. 56) Koike, S.; Yamaji, K.; Yanagisawa,− T. Effect of the medium-range 28) Ceperley, D.; Chester, G.V.; Kalos, K.H. Monte Carlo simulation of a transfer energies to the superconductivity in the two-chain Hubbard many-fermion study. Phys. Rev. B 1977, 16, 3081. model. J. Phys. Soc. Jpn. 1999, 68, 1657–1663. 29) Gros, C.; Joynt, R.; Rice, T.M. Antiferromagnetic correlations in 57) Mott, N.F. Metal-Insulator Transitions; Taylor and Francis Ltd.: Lon- almost-localized Fermi liquids. Phys. Rev. B 1987, 36, 381. don, UK, 1974. 30) Yokoyama, H.; Shiba, H. Variational Monte Carlo studies of Hubbard 58) Moriya, T. Spin Fluctuations in Itinerant Electron Magnetism; model I. J. Phys. Soc. Jpn. 1987, 56, 1490–1506. Springer: Berlin, Germany, 1985. 31) Giamarchi, T.; Lhuillier, C. Phase diagrams of the two-dimensional 59) Yosida, K. Theory of Magnetism; Springer: Berlin, Germany, 1996. Hubbard and t-J models by a variational Monte Carlo study. Phys. Rev. 60) Tranquada, J.; J.; Axe, J. D.; Ichikawa, N.; Nakamura, Y.; Uchida, S.; B 1991, 43, 12943. Nachumi, B. Neutron-scattering study of stripe-phase order of holes 32) Zhang, S.; Carlson, J.; Gubernatis, J.E. Constrained path Monte Carlo and spins in La1.48Nd0.4Sr0.12CuO4. Phys. Rev. B 1996, 54, 7489. method for fermion ground states. Phys. Rev. B 1997, 55, 7464. 61) Suzuki, T.; Goto, T.; Shinoda, T.; Fukase, T.; Kimura, H.; Yamada, 33) Zhang, S.; Carlson, J.; Gubernatis, J.E. Pairing correlation in the two- K.: Ohashi, M.; Yamaguchi, Y. Observation of modulated magnetic dimensional Hubbard model. Phys. Rev. Lett. 1997, 78, 4486. long-range order in La1.88Sr0.12CuO4. Phys. Rev. B 1998, 57, R3229. 34) Yanagisawa, T.; Shimoi, Y. Exact results in strongly correlated elec- 62) Yamada, K.; Lee, C. H.; Kurahashi, K.; Wada, J.; Wakimoto, S.; Ueki, trons. Int. J. Mod. Phys. 1996, B10, 3383–3450. S.; Kimura, H.; Endoh, Y.; Hosoya, S.; Shirane, G.; Birgeneau, R. J.; 35) Yanagisawa, T.; Shimoi, Y.Ground state of the Kondo-Hubbard model Greven, M.; Kastner, M. A.; Kim, Y. J. Doping dependence of the spa- at half-filling. Phys. Rev. Lett. 1995, 74, 4939. tially modulated dynamical spin correlations and the superconducting- 36) Nakanishi, T.; Yamaji, K.; Yanagisawa, T. Variational Monte Carlo transition temperature in La2 xSrxCuO4. Phys. Rev. B 1998, 57, 6165. indications of d-wave superconductivity in the two-dimensional Hub- 63) Arai, M.; Nishijima, T.; Endoh,− Y.; Egami, T.; Tajima, S.; Tomi- bard model. J. Phys. Soc. Jpn. 1997, 66, 294–297. moto, K.; Shiohara, Y.; Takahashi, M.; Garrett, A.; Bennington, S. 37) Yamaji, K.; Yanagisawa, T.; Nakanishi, T.; Koike, S. Variational M. Incommensurate spin dynamics of underdoped superconductor Monte Carlo study on the superconductivity in the two-dimensional YBa2Cu3Y6.7. Phys. Rev. Lett. 1999, 83, 608. Hubbard model. Physica C 1998, 304, 225–238. 64) Mook, H.A.; Pengcheng, D.; Dogan, F.; Hunt, R. D. One-dimensional 38) Yamaji, K.; Yanagisawa,T.; Koike, S. Bulk limit of superconducting nature of the magnetic fluctuations in YBa2Cu3O6.6. Nature 2000, condensation energy in 2D Hubbard model. Physica B 2000, 284–288, 404, 729. 415–416. 65) Wakimoto, S.; Birgeneau, R. J.; Kastner, M. A.; Lee, Y. S.; Erwin, 39) Yamaji, K.; Yanagisawa, T.; Miyazaki, M.; Kadono, R. Superconduct- R.; Gehring, P. M.; Lee, S. H.; Fujita, M.; Yamada, K.; Endoh, Y.; ing condensation energy of the two-dimensional Hubbard model in the Hirota, K.; Shirane, G. Direct observation of a one-dimensional static large-negative-t’ region. J. Phys. Soc. Jpn. 2011, 80, 083702. spin modulation in insulating La1.95Sr0.05CuO4. Phys. Rev. B 2000, 61, 40) Hardy, T.M.; Hague, P.; Samson, J.H.; Alexandrov, A.S. Superconduc- 3699. tivity in a Hubbard-Fröhlich model in cuprates. Phys. Rev. B 2009, 79, 66) Bianconi, A.; Saini, N. L.; Lanzara, A.; Missori, M.; Rossetti, T.; Oy- 212501. anagi, H.; Yamaguchi, H.; Oka, K.; Ito, T. Determination of the local 41) Yanagisawa, T.; Miyazaki, M.; Yamaji, K. Correlated-electron systems lattice distortions in the CuO2 plane of La1.85Sr0.15CuO4. Phys. Rev. and high-temperature superconductivity. J. Mod. Phys. 2013, 4, 33. Lett. 1996, 76, 3412.

9 Condens. Matter 2019, 4, 57 T Yanagisawa

67) Bianconi, A. Quantum materials: Shape resonances in superstripes. in multi-band superconductors. J. Phys. Soc. Jpn. 2012, 81, 024712. Nat. Phys. 2013, 9, 536. 95) Hu, X.; Wang, Z. Stability and Josephson effect of time-reversal- 68) Hoffman, J.E.; McElroy, K.; Lee, D.-H.; Lang, K. M.; Eisaki, symmetry-broken multicomponent superconductivity induced by frus- H.; Uchida, S.: Davis, J. C. Imaging quasiparticle interference in trated intercomponent coupling. Phys. Rev. B 2012, 85, 064516. Bi2Sr2CaCu2O8+δ. Science 2002, 295, 466. 96) Stanev, V. Model of collective modes in three-band superconductors 69) Wise, W.D.; Boyer, M. C.; Chatterjee, Kamalesh; Kondo, T.; Takeuchi, with repulsive interband interactions. Phys. Rev. B 2012, 85, 174520. T.; Ikuta, J.; Wang Yayu; Hudson, E. W. Charge-density-wave origin of 97) Platt, C.; Thomale, R.; Homerkamp, C.; Zhang, S.C.; Hanke, W. cuprate checkerboard visualized by scanning tunnelling microscopy. Mechanism for a pairing with time-reversal symmetry breaking in Nat. Phys. 2008, 4, 696. iron-based superconductors. Phys. Rev. B 2012, 85, 180502. 70) Hanaguri, T.; Lupien, C.; Kohsaka, Y.; Lee, D.-H.; Azuma, M.; 98) Maiti, S.; Chubukov, A.V. s + is state with broken time-reversal sym- Takano, M.; Takagi, H.; Davis, J. A checkerboard electronic crys- metry in Fe-based superconductors. Phys. Rev. B 2013, 87, 144511. tal state in lightly hole-doped Ca2 xNaxCuO2Cl2. Nature 2004, 430, 99) Wilson, B.J.; Das, M.P. Time-reversal-symmetry-broken state in the 1001. − BCS formalism for a multi-band superconductor. J. Phys. Condens. 71) Miyazaki, M.; Yanagisawa, T.; Yamaji, K. Checkerboard states in the Matter 2013, 25, 425702. two-dimensional Hubbard model with the Bi2212-type band. J. Phys. 100) Ganesh, R.; Baskaran, G.; van den Brink, J.; Efremov, D.V. Theoreti- Soc. Jpn. 2009, 78, 043706. cal prediction of a time-reversal broken chiral superconducting phase 72) Anderson, P.W. The resonating valence bond states in La2CuO4. Sci- driven by electronic correlations in a single TiSe2 layer. Phys. Rev. ence 1987, 253, 1196. Lett. 2014, 113, 177001. 73) Yanagisawa, T.; Koike, S.; Yamaji, K. Off-diagonal wave function 101) Yerin, Y.S.; Omelyanchouk, A.N.; Il’ichev, E. Dc SQUID based on a Monte Carlo Studies of Hubbard model I. J. Phys. Soc. Jpn. 1998, 67, three-band superocnductor with broken time-reversal symmetry. Su- 3867–3874. perconduct. Sci. Technol. 2015, 28, 095006. 74) Yanagisawa, T.; Miyazaki, M. Mott transition in cuprate high- 102) Hillier, A.D.; Quintanilla, J.; Cywinskii, R. Evidence for time- temperature superconductors. EPL 2014, 107, 27004. reversal symmetry breaking in the noncentrosymmetric superconduc- 75) Yanagisawa, T. Antiferromagnetism, superconductivity and phase di- tor LaNiC2. Phys. Rev. Lett. 2009, 102, 117007. agram in the two-dimensional Hubbard model—Off-diagonal wave 103) Hase, I.; Yanagisawa, T. Electronic structure of RNiC2 (R=La, Y, and function Monte Carlo studies of Hubbard model III—. J. Phys. Soc. Th). J. Phys. Soc. Jpn. 2009, 78, 084724. Jpn. 2019. 88, 054702. 104) Yanagisawa, T.; Hase, I. Massless modes and abelian gauge fields in 76) Eliashberg, G.M. Interactions between electrons and lattice vibrations multi-band superconductors. J. Phys. Soc. Jpn. 2013, 82, 124704. in a superconductor. Sov. Phys. JETP 1960, 11, 696–702. 105) Lin, S.Z.; Hu, X. Phase solitons in multi-band superconductors with 77) Carbotte, J.P. Properties of boson-exchange superconductors. Rev. and without time-reversal symmetry. New J. Phys. 2012, 14, 063021. Mod. Phys. 1990, 62, 1027. 106) Kobayashi, K.; Machida, M.; Ota, Y.; Nori, F. Massless collective ex- 78) McMillan, W.L. Transition temperature of strong-coupled supercon- citations in frustrated multiband superconductors. Phys. Rev. B 2013, ductors. Phys. Rev. 1968, 167, 331. 88, 224516. 79) Allen, P.B.; Dynea, R.C. Transition temperature of strong-coupled su- 107) Koyama, T. Collective modes in multiband superconductors: Rigorous perconductors reanalyzed. Phys. Rev. B 1968, 12, 905. study based on the Ward-Takahashi relations. J. Phys. Soc. Jpn. 2014, 80) Ashcroft, N.W. Metallic hydrogen: A high-temperature superconduc- 83, 074715. tor? Phys. Rev. Lett. 1968, 21, 1748. 108) Yanagisawa, T.; Tanaka, Y. Fluctuation-induced Nambu-Goldstone 81) Drozdov, A.P.; Eremets, M.I.; Troyan, I.A.; Ksenofontov, V.; Shylin, bosons in a Higs-Josephson model. New J. Phys. 2014, 16, 123014. S.I. Conventional superconductivity at 203 kelvin at high pressures in 109) Tanakai, Y.; Hase, I.; Yanagisawa, T.; Kato, G.; Nishio, T.; Arisawa, the sulfur hydride system. Nature 2015, 525, 73. S. Current-induced massless mode of the interband phase difference 82) Peng, F.; Sun, Y.; Pickard, C.J.; Needs, R.J.; Wu, Q.; Ma, Y. Hydro- in two-band superconductors. Physica C 2015, 516, 10. gen clathrate structures in rare earth hydrides at high pressures: Possi- 110) Valletta, A.; Bianconi, A.; Perali, A.; Saini, N.L. Electronic and su- ble route to room-temperature superconductivity. Phys. Rev.Lett. 2017, perconducting properties of a superlattice of quantum stripes at the 119, 107001. atomic limit. Z. Physik B Condens. Matter 1997, 104, 707. 83) Liu, H.; Naumov, I.I.; Hoffmann, R.; Schcroft, N.W.; Hemley, R.J. 111) Choi, H.Y.; Yun, J.H.; Bang, Y.; Lee, H.C. Model for the inverse iso- Potential high-Tc superconducting lanthanum and yttrium hydrides at tope effect of FeAs-based superconductors in the π-phase-shifted pair- high pressure. Proc. Natl. Acad. Sci. USA 2017, 114, 6990–6995. ing state. Phys. Rev. B 2009, 80, 052505. 84) Innocenti, D.; Poccia, N.; Ricci, A.; Valletta, A.; Caprara, S.; Perali, 112) Shirage, P.M.; Kihou, K.; Miyazawa, K.; Lee, C.-H.; Kito, H.; Eisaki, A.; Bianconi, A. Resonant and crossover phenomena in a multiband H.; Yanagisawa, T.; Tanaka, Y.; Iyo, A. Inverse iron isotope effect on superconductor: Tuning the chemical potential near a band edge. Phys. the transition temperature of the (Ba,K)Fe2As2 superconductor. Phys. Rev. B 2010, 82, 184528. Rev. Lett. 2009, 103, 257003. 85) Mazziotti, M.V.; Valletta, A.; Campi, G.; Innocenti, D.; Perali, A.; 113) Yanagisawa, T.; Odagiri, K.; Hase, I.; Yamaji, K.; Shirage, P. M.; Bianconi, A. Possible Fano resonance for high-Tc multi-gap supercon- Tanaka, Y.; Iyo, A.; Eisaki, H. Isotope effect in multi-band and multi- ductivity in p-Terphenyl doped at the Lifshitz transition. EPL 2017, channel attractive systems and inverse isotope effect in iron-based su- 118, 37003. perconductors. J. Phys. Soc. Jpn. 2009, 78, 094718. 86) Moskalenko, V.A.; Fiz. Metal and Metallored. 1959, 8, 2518. 114) Perali, A.; Innocenti, D.; Valletta, A.; Bianconi, A. Anomalous isotope 87) Suhl, H.; Mattias, B.T.; Walker, L.R. Bardeen-Cooper-Schrieffer the- effect near a 2.5 Lifshitz transition in a multi-band multi-condensates ory of superconductivity in the case of overlapping bands. Phys. Rev. superconductor made of a superlattice of stripes. Superconduct. Sci. Lett. 1959, 3, 552. Technol. 2012, 25, 124002. 88) Peretti, J. Superconductivity of transition elements. Phys. Lett. 1962, 115) Izyumov, Yu. A.; Laptev, V.M. Vortex structure in superconductors 2, 275. with a many-component order parameter. Phase Transit. 1990, 20, 95. 89) Kondo, J. Superconductivity in transition metals. Prog. Theor. Phys. 116) Volovik, G.E. The Universe in a Helium Droplet; Oxford University 1963, 29, 1–9. Press: Oxford, UK, 2009. 90) Stanev, V.; Tesanovic, Z. Three-band superconductivity and the order 117) Kuplevakhsky, S.V.; Omelyanchouk, A.N.; Yerin, Y.S. Soliton states in parameter that breaks time-reversal symmetry. Phys. Rev. B 2010, 81, mesoscopic two-band superconducting cylinders. J. Low Temp. Phys. 134522. 2011, 37, 667. 91) Tanaka, Y.; Yanagisawa, T. Chiral ground state in three-band super- 118) Tanaka, Y.; Yamamori, H.; Yanagisawa, T.; Nishio, T.; Arisawa, S. conductors. J. Phys. Soc. Jpn. 2010, 79, 114706. Experimental formation of a fractional vortex in a superconducting 92) Tanaka, Y.; Yanagisawa, T. Chiral state in three-gap superconductors. bi-layer. Physica C 2018, 548, 44. Solid State Commun. 2010, 150, 1980–1982. 119) Yanagisawa, T.; Hase, I.; Tanaka, Y. Massless and quantized modes of 93) Dias, R.G.; Marques, A.M. Frustrated multiband superconductivity. kinks in the phase space of superconducting gaps. Phys. Lett. A 2018, Superconduct. Sci. Technol. 2011, 24, 085009. 382, 3483. 94) Yanagisawa, T.; Tanaka, Y.; Hase, I.; Yamaji, K. Vortices and chirality 120) Littlewood, P.B.; Varma, C.M. Amplitude collective modes in super-

10 Condens. Matter 2019, 4, 57 T Yanagisawa

conductors and their coupling to charge-density wave. Phys. Rev. B University Press: Cambridge, UK, 1993. 1982, 26, 4883. 146) Onuki, Y. Physics of Heavy Fermions: Heavy Fermions and Strongly 121) Cea, T.; Benfatto, L. Nature and Raman signature of the Higgs am- Correlated Electron Systems; World Scientific Pub Co Inc.: Singapore, plitude modes in the coexisting superconducting and charge-density- 2018. wave mode. Phys. Rev. B 2014, 90, 224515. 147) Ishiguro, T.; Yamaji, K.; Saito, G. Organic Superconductors; Springer: 122) Pekker, D.; Varma, C.M. Amplitude/Higgs modes in condensed matter Berlin, Germany, 2012. physics. Annu. Rev. Condens. Matter Phys. 2015, 6, 269–297. 148) Akashi, R.; Kawamura, M.; Tsuneyuki, S.; Nomura, Y.; Arita, R. First- 123) Cea, T.; Castellani, C.; Seibold, G.;Benfatto, L. Nonrelativistic dy- principles study of the pressure and crystal-structure dependences of namics of the amplitude (Higgs) mode in superconductors. Phys. Rev. the superconducting transition temperature in compressed sulfur hy- Lett. 2015, 115, 157002. drides. Phys. Rev. B 2015, 91, 224513. 124) Koyama, T. Perturbative approach to the collective modes in the 149) Steglich, R.; Asrts, J.; Bredl, C. D.; Lieke, W.; Meschede, D.; Franz, TRSB phase of multiband superconductors. J. Phys. Soc. Jpn. 2016, W.; Schäfer, H. Superconductivity in the presence of strong Pauli para- 85, 064715. Nambu-Goldstone-Leggett and Higgs modes in multi- magnetism: CeCu2Si2. Phys. Rev. Lett. 1979, 43, 1892. condensate superconductors. J. Superconduct. Novel Magn. 2016, 29, 150) Kittaka, S.; Aoki, Y.; Shimura, Y.; Sakakibara, T.; Seiro, S.; Geibel, 3099. C.; Steglich, F.; Ikeda, H.; Machida, K. Multiband superconductivity 125) Yanagisawa, T. Fluctuation modes in multi-gap superconductors. with unexpected deficiency of nodal quasiparticles in CeCu2Si2. Phys. In Vortices and Nanostructured Superconductors; Crisan, A., Ed.; Rev. Lett. 2014, 112, 0607002. Springer: Berlin, Germany, 2017. 151) Stewart, G.R.; Fisk, Z.; Smith, J.L.; Willis, J.O.; Wire, M.S. New 126) Yanagisawa, T. Nambu-Goldstone bosons characterized by the order heavy-fermion system NpBe13 with a comparison to UBe13 and parameters in spontaneous symmetry breaking. J. Phys. Soc. Jpn. 2017, PuBe13. Phys. Rev. Lett. 1984, 52, 679. 86, 104711. 152) Ott, H.R.; Rudingier, H.; Delsing, P.; Fisk, Z. Magnetic ground state 127) Aitchison, I.J.R.; Ao, P.; Thouless, D.; Zhu, X.M. Effective La- of a heavy-electron system: U2Zn17. Phys. Rev. Lett. 1984, 52, 1551. grangian for BCS superconductors at T = 0. Phys. Rev. B 1995, 51, 153) Palstra, T.T.; Menovsky, A.A.; van den Berg, J. Superconducting and 6531. magnetic transitions in the heavy-fermion system URu12. Phys. Rev. 128) Murotani, Y.; Tsuji, N.; Aoki, H. Theory of light-induced resonances Lett. 1985, 55, 2727. with collective Higgs and Leggett modes in multiband superconduc- 154) Amitsuka, H.; Sato, M.; Metoki, N.; Yokoyama, M.; Kuwahara, K.; tors. Phys. Rev. B 2017, 95, 104503. Sakakibara, T.; Morimoto, H.; Kawarazaki, S.; Miyako, Y.; Mydosh, 129) Yanagisawa, T. Theory of Green’s functions of Nambu-Goldstone and J.A.Effect of pressure on tiny antiferromagnetic moment in the heavy- Higgs modes in superconductors. J. Superconduct. Novel Magn. 2019, electron compound URu2Si2. Phys. Rev. Lett. 1999, 83, 5114. doi:10.1007/s10948-018-4983-2. 155) Ohkuni, H.; nada, Y.; Tokiwa, Y.; Sakurai, K.; Settai, R.; Honma, 130) Perali, A.; Bianconi, A.; Lanzara, A.; Saini, N.L. The gap amplifi- T.; Haga, Y.; Yamamoto, E.; Onuki, Y.; Yamagami, H.; Takahashi, cation at a shape resonance in a superlattice of quantum stripes: A S.; Yanagisawa, T. Fermi surface properties and de Haas-van Alphen mechanism for high Tc. Solid State Commun. 1996, 100, 181–186. oscillation in both the normal and superconducting mixed states of 131) Bianconi, A.; Valletta, A.; Perali, A.; Saini, N.L. Superconductivity of URu2Si2. Philos. Mag. 1999, B79, 1045. s striped phase at the atomic limit. Physica C 1998, 296, 269–280. 156) Hedo, M.; Inada, Y.; Yamamoto, E.; Haga, Y.; Onuki, Y.; Aoki, Y.; 132) Kusmartsev, F.V.; Di Castro, D.; Bianconi, G.; Bianconi, A. Transfor- Matsuda, D.; Sato, H.; Takahashi, S. Superconducting properties of mation of strings into an inhomogeneous phase of stripes and itinerant CeRu2. J. Phys. Soc. Jpn. 1998, 67, 272–279. carriers. Phys. Lett. A 2000, 275, 118. 157) Geibel, G.; Schank, C.; Thies, S.; Kitazawa, H.; Bredl, C. D.; Bohm, 133) Müller, K.A.; Zhao, G.M.; Conder, K.; Keller, H. The ratio of small A.; Rau, M.; Grauel, A.; Caspary, R.; Helfrich, R.; Ahlheim, U.; We- polarons to free carriers in derived from susceptibility measurements. ber, G.; Steglich, F. Heavy-fermion superconductivity at T=2K in the J. Phys. Condens. Matter 1998, 10, L291. antiferromagnet UPd2Al3. Z. Phys. B 1991, 84, 1–2. 134) Bianconi, A. On the Fermi liquid coupled with a generalized Wigner 158) Kyogaku, M.; Kitaoka, Y.; Asayama, K.; Geibel, C.; Schank, C.; polaronic CDW giving high Tc superconductivity. Solid State Com- Steglich, F. NMR and NQR studies of magnetism and superconductiv- mun. 1994, 91, 1–5. ity in the antiferromagnetic heavy fermion superconductors UM2Al3 135) Harris, A.B.; Lange, R.V. Single-particle excitations in narrow energy (M=Ni and Pd). J. Phys. Soc. Jpn. 1993, 62, 4016–4030. bands. Phys. Rev. 1967, 157, 295. 159) Inada, Y.; Yamagami, H.; Haga, Y.; Sakurai, K.; Tokiwa, Y.; Honma, 136) Chao, K.A.; Spalek, J.; Ole´s, A.M. Kinetic exchange interaction in a T.; Yamamoto, E.; Onuki, Y.; Yanagisawa, T. Fermi surface and de narrow s-band. J. Phys. C 1977, 10, L271. Haas-van Alphen oscillation in both the normal and superconducting 137) Chao, K.A.; Spalek, J.; Ole´s, A.M. Canonical perturbation expansion mixed states of UPd2Al3. J. Phys. Soc. Jpn. 1999, 68, 3643–3654. of the Hubbard model. Phys. Rev. B 1978, 18, 3453. 160) Ishida, K.; Mukuda, H.; Kitaoka, Y.; Asayama, K.; Mao, Z. Q.; Mori, 138) Micnas, R.; Ranninger, J.; Robaszkiewicz, S. Superconductivity in Y.; Maeno, Y. Spin-triplet superconductivity in UNi2Al3 revealed by narrow-band systems with local nonretarded attractive interactions. the 27Al knight shift measurement. Phys. Rev. Lett. 2002, 89, 037002. Rev. Mod. Phys. 1990, 62, 113. 161) Petrovic, C.; Pagliuso, P. G.; Hundley, M. F.; Movshovich, R.; Sarrao, 139) Robaszkiewicz, S.; Pawlowski, G. Effects of finite pair binding energy J. L.; Thompson, J. D.; Fisk, Z.; Monthoux, P. heavy-fermion super- in a model for a superconductor with local electron pairing. Physica C conductivity in CeCoIn5 at 2.3 K. J. Phys. Condens. Matter 2001, 13, 1993, 210, 61–79. L337. 140) Arrigoni, E.; Strinati, G.C. Doping-induced incommensurate antifer- 162) Izawa, K.; Yamaguchi, H.; Matsuda, Y.; Shishido, H.; Settai, R.; romagnetism in a Mott-Hubbard insulator. Phys. Rev. B 1991, 44, Onuki, Y. Angular position of nodes in the superconducting gap of 7455. quasi-2D heavy-fermion superconductor CeCoIn5. Phys. Rev. Lett. 141) Kapcia, K.; Robaszkiewicz, S.; Micnas, R. Phase separation in a lattice 2001, 87, 057002. model of a superconductor with pair hopping. J. Phys. Condens. Matter 163) Hegger, H.; etrovic, C.; Moshopoulou, E. G.; Hundley, M. F.; Sarrao, 2012, 24, 215601. J. L.; Fisk, Z.; Thompson, J. D. Pressure-induced superconductivity in 142) Feiner, L.F.; Jefferson, J.H.; Raimondi, R. Effective single-band mod- quasi-2D CeRhIn5. Phys. Rev. Lett. 2000, 84, 4986. els for high-Tc cuprates I Coulomb interactions. Phys. Rev. B 1996, 53, 164) Movshovich, R.; Graf, T.; mandrus, D.; Thompson, J. D.; Smith, J. L.; 8751. Fisk, Z. Superconductivity in heavy-fermion CeRh2Si2. Phys. Rev. B 143) Kuroki, K.; Onari, S.; Arita, R.; Usui, H.; Tanaka, Y.; Kontani, H.; 1996, 53, 8241. Aoki, H. Unconventional Pairing Originating from the Disconnected 165) Saxena, S.S.; Agarwal, P.; Ahilan, K.; Grosche, F. M.; Haselwimmer, Fermi Surfaces of Superconducting LaFeAsO1 xFx. Phys. Rev. Lett. R. K. W.; Steiner, M. J.; Pugh, E.; Walker, I. R.; Julian, S. R.; Mon- − 2008, 101, 087004. thoux, P.; Lonzarich, G. G.; Huxley, A.; Sheikin, I.; Braithwaite, D.; 144) Grewe, N.; Steglich, F. Heavy fermions. In Handb. Phys. Chem. Of Flouquet, J. Superconductivity on the border of itinerant-electron fer- Rare Earths; Volume 14, p. 343; North-Holland; Amsterdam, Nether- romagnetism in UGe2. Nature 2000, 406, 587. lands, 1991. 166) Aoki, D.; Huxley, A.; Ressouche, E.; Braithwaite, D.; Flouquet, J.; 145) Hewson, A.C. The Kondo Problem to Heavy Fermions; Cambridge Brison, J.-P.; Lhotel, E.; Paulsen, C. Coexistence of superconductivity

11 Condens. Matter 2019, 4, 57 T Yanagisawa

and ferromagnetism in URhGe. Nature 2001, 413, 613. formed by hole doping in the high-Tc superconductors. Physica C 167) Maeno, Y.; Hashimoto, H.; Yoshida, K.; Nishizaki, S.; Fujita, T.; Bed- 1989, 160, 424–430. norz, J.G.; Lichtenberg, F. Superconductivity in a layered perovskite 184) McMahan, A.K.; Annett, J.F.; Martin, R.M. Cuprate parameters from without copper. Nature 1994, 372, 532. numerical Wannier functions. Phys. Rev. B 1990, 42, 6268. 168) Bauer, E.D.; Frederick, N. A.; Ho, P.-C.; Zapf, V. S.; Maple, M. B. 185) Eskes, H.; Sawatzky, G. Single-, triple-, or multiplel-band Hubbard Superconductivity and heavy fermion behavior in PrOs4Sb12. Phys. models. Phys. Rev. B 1991, 43, 119. Rev. B 2002, 65, R100506. 186) Otsuka, H. Variational Monte Carlo studies of the Hubbard model in 169) Takada, K.; Sakurai, H.; Takayama-Muromachi, E.; Izumi, F.; Dila- one- and two-dimensions. J. Phys. Soc. Jpn. 1992, 61, 1645–1656. nian, F. A.; Sakai, T. Superconductivity in two-dimensional CoO2 lay- 187) Yanagisawa, T.; Koike, S.; Yamaji, K. d-wave state with multiplicative ers. Nature 2003, 422, 53. correlation factors for the Hubbard model. J. Phys. Soc. Jpn. 1999, 68, 170) Chaillout, C.; Remeika, J.P.; Santoro, A.; Mareizo, M. The determina- 3608–3614. tion of the Bi valence state in BaBiO3 by neutron powder diffraction 188) Eichenberger, D.; Baeriswyl, D. Superconductivity and antiferromag- data. Solid State Commun. 1985, 56, 829–831. netism in the-dimensional Hubbard model: A variational study. Phys. 171) Nagamatsu, J.; Nakagawa, N.; Muranaka, T.; Zenitani, Y.; Akimitsu, J. Rev. B 2007, 76, 180504. Superconductivity at 39 K in magnesium diboride. Nature 2001, 410, 189) Baeriswyl, D.; Eichenberger, D.; Menteshashvii, M. Variational 63. ground states of the two-dimensional Hubbard model. New J. Phys. 172) Kamihara, Y.; Watanabe, T.; Hirano, M.; Hosono, H. Iron-based lay- 2009, 11, 075010. ered superconductor La(O1 xFx)FeAs (x = 0.05 0.12) with Tc = 26 190) Baeriswyl, D. Superconductivity in the repulsive Hubbard model. J. − − K. J. Am. Chem. Soc. 2000, 130, 3296. Superconduct. Novel Magn. 2011, 24, 1157–1159. 173) Kito, H.; Eisaki, H.; Iyo, A. Superconductivity at 54 K in F-free 191) Yanagisawa, T. Quantum Monte Carlo diagonalization for many- NdFeAsO1 y. J. Phys. Soc. Jpn. 2008, 77, 063707. fermion systems. Phys. Rev. B 2007, 75, 224503. − 174) Drozdov, A.P.; Minkov, V.S.; Besedin, S.P.; Kong, P.P.; Kuzovnikov, 192) Misawa, T.; Imada, M. Origin of high-Tc superconductivity in doped M.A.; Knyazev, D.A.; Eremets, M.I. Superconductivity at 250 K in Hubbard models and their extensions: Roles of uniform charge fluctu- lanthanum hydride at high pressures. Nature 2019, 569, 528. 175) Somayazulu, M.; Ahart, M.; Mishra, A.K.; Geballe, Z.M.; Baldini, M.; Meng, Y.; Struzhkin, V.V.; Hemley, R.J. Evidence for superconduc- ations. Phys. Rev. B 2014, 90, 115137. tivity above 260 K in lanthanum superhydride at megabar pressures. 193) Yanagisawa, T.; Miyazaki, M.; Yamaji, K. Crossover induced electron Phys. Rev. Lett. 2019, 122, 027001. pairing and superconductivity by kinetic renormalization in correlated 176) Bickers, N.E.; Scalapino, D.J.; White, S.R. Convserving Approxima- electron systems. Condens. Matter 2018, 3, 26. tions for Strongly Correlated Electron Systems: Bethe-Salpeter Equa- 194) Kondo, J. The Physics of Dilute Magnetic Alloys; Cambridge Univer- tion and Dynamics for the Two-Dimensional Hubbard Model. Phys. sity Press: Cambridge, UK, 2012. Rev. Lett. 1989, 62, 961. 195) Yanagisawa, T. Kondo effect in the presence of spin-orbit coupling. J. 177) Pao, C.-H.; Bickers, N.E. Anisotropic superconductivity in the 2D Phys. Soc. Jpn. 2012, 81, 094713. Hubbard model: Gap function and interaction weight. Phys. Rev. B 196) Yanagisawa, T. Kondo effect in Dirac systems. J. Phys. Soc. Jpn. 2015, 1994, 49, 1586. 84, 074705. 178) Monthoux, P.; Scalapino, D.J. Self-consistent d2 d2 pairing in a two- 197) Ellis, R.K.; Stirling, W.J.; Webber, B.R. QCD and Collider Physics; x − y dimensional Hubbard model. Phys. Rev. Lett. 1994, 72, 1874. Cambridge University Press: Cambridge, UK, 1996. 179) Loram, J.W.; Mirza, K.A.; Cooper, J.R.; Liang, W.Y. Electronic spe- 198) Nozieres, P.; Schmitt-Rink, S. Bose condensation in an attractive fermi cific heat of YBa2Cu3O6+x from 1.8 to 300 K. Phys. Rev. Lett. 1993, gas: From weak to strong coupling superconductivity. J. Low Temp. 71, 1740. Phys. 1985, 59, 195–211. 180) Hao, Z.; Clem, J.R.; McElfresh, M.W.; Civale, L.; Malozemoff, A.P.; 199) Rajaraman, R. Solitons and Instantons; North-Holland: Amsterdam, Holtzberg, F. Model for the reversible magnetization of high-κ type-II The Netherlands, 1989. superconductors: Application to high-Tc superconductors. Phys. Rev. 200) Solyom, J. The Fermi gas model of one-dimensional conductors. Adv. B 1991, 43, 2844. Phys. 1979, 28, 201–303. 181) Weber, C.; Haule, K.; Kotliar, G. Critical weights and waterfalls in 201) Yanagisawa, T. Chiral sine-Gordon model. EPL 2016, 113, 41001. doped charge-transfer insulators. Phys. Rev. B 2008, 78, 134519. 202) Yanagisawa, T. Renormalization group analysis of the hyperbolic sine- 182) Hybertsen, M.S.; Schlüter, M.; Christensen, N.E. Calculation of Gordon model. Prog. Theor. Exp. Phys. 2019, 2019, 023A01. Coulomb-interaction parameter for La2CuO4 using a constrained- 203) Gross, D.; Neveu, A. Dynamical Symmetry breaking in asymptotically density-functional approach. Phys. Rev. B 1989, 39, 9028. free field theories. Phys. Rev. D 1974, 10, 3235. 183) Eskes, H.; Sawatzky, G.A.; Feiner, L.F. Effective transfer for singlets