• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Theory of unconventional superconductivity in nickelate-based materials?

    2021-10-28 07:03:24MingZhang張銘YuZhang張渝HuaimingGuo郭懷明andFanYang楊帆
    Chinese Physics B 2021年10期
    關(guān)鍵詞:楊帆

    Ming Zhang(張銘) Yu Zhang(張渝) Huaiming Guo(郭懷明) and Fan Yang(楊帆)

    1School of Physics,Beijing Institute of Technology,Beijing 100081,China

    2Institute for Quantum Science and Engineering and Department of Physics,Southern University of Science and Technology,Shenzhen 518055,China

    3School of Physics,Beihang University,Beijing 100191,China

    Keywords: nickelate superconductor,pairing symmetries

    1. Introduction

    The discovery of superconductivity(SC)with highestTcup to 15 K in the Sr-doped infinite-layer nickelate NdNiO2[1,2]has attracted great research interests both experimentally[3–7]and theoretically.[4,8–26]This superconductor was thought to be a Ni-based analogy of the hight-temperature cuprates superconductors:[2]the 3d9electronic configuration of the Ni+in the NdNiO2is similar to that of the Cu2+in the cuprates; both compounds have isostructural CuO2/NiO2plane withD4symmetry. Density functional theory (DFT)based calculations[12]suggest that the electron–phonon coupling mediated pairing mechanism cannot account for theTcas high as 15 K, and the SC in this material is more likely driven by electron–electron interactions. More remarkably, recent doping study over the parent compound reveals a dome-shaped dependence of the superconductingTcover the doping concentration,[5,6]which is similar to that of the cuprates. These similarities between the Sr-doped NdNiO2and the cuprates superconductors imply that the study on the former is helpful for the clarification of the pairing mechanism of the high-temperature cuprate superconductors.[27]

    There are, however, significant differences between the cuprates and the infinite-layer nickelate.[2]One difference lies in the presence of the small Fermi-pockets around the Aand possibly theΓ-points in the Brillouin zone in the parent compound NdNiO2, which are mainly contributed from the 5d electrons of the Nd.[8,11–14,19,22–26]These Nd 5d-stateoriginated Fermi pockets lead to self doping into the shouldbe Mott-insulating NiO2plane with extra holes. Note that the charge-transfer gap between the Ni 3d states and the O 2p states in the NdNiO2is much larger than that in the cuprates,[2,7,9,21,24]which indicts that the extra holes in the NiO2plane should locate on the Ni 3d orbitals instead of the O 2p ones. Therefore, the Zhang–Rice singlet[28]present in the cuprates is absent here. These conducting Nd 5d electrons and the extra holes in the NiO2plane make the NdNiO2metallic. The other obvious difference lies in that, unlike the cuprates, the parent compound NdNiO2is not antiferromagnetic (AFM) ordered.[2,29,30]One possible reason lies in that the AFM superexchange interaction in the NdNiO2is much weaker than that in the cuprates.[2,7]The other reason lies in that the Kondo coupling between the conduction Nd 5d electrons and the localized Ni 3d electrons leads to formation of the Kondo singlet at low temperature,[2,13,14,16,23,24]which suppresses the AFM order. Therefore,the main difference between the cuprates and the infinite-layer nickelate might lies in the presence of the low-lying Nd 5d degree of freedom in the parent compound. Such difference will be largely suppressed with the enhancement of the hole doping,[11,12,15,25]as the participation of the Nd 5d electrons in the low-energy physics simply takes the form of electron pockets, which will vanish or shrink upon hole doping. Therefore,as long as the pairing mechanism and pairing symmetry in the Sr-doped compound are considered,the Ni 3d electrons should take the main role.

    The band structure of the NdNiO2has been studied previously.[31–33]To make the succeeding calculations considering electron–electron interactions realizable,a tight-binding(TB) model with finite-range hopping integrals is required.In the TB model introduced in Ref. [22], all the orbitals involved on the FS are considered, which leads to a seven- or six-orbital model. Such a model can provide a perfect description of the band structure of NdNiO2. Then a fluctuationexchange study on this model yields the dx2?y2pairing symmetry. Here to make the physics more simple and clear, we start from a two-orbital TB model provided in Ref.[4]. In this model, only the Nd 5d orbital and the Ni 3dx2?y2orbital are considered. As the hybridization between the two orbitals is weak,[4]we perform two comparison studies on the model: in one study we adopt the realistic band structure with weak coupling between the two orbitals, in the other one we turn off the coupling and study the imaginary comparison model. Our random-phase-approximation (RPA) based calculations yield very similar results for the two studies: both studies yield the dx2?y2pairing symmetry consistent with that obtained in other studies,[11,22,24]and the spin correlations obtained in the two studies are very similar: both studies reveal strong inplane spin fluctuations,consistent with Refs.[17,20]. In both studies, either the Cooper pairing or the spin fluctuations are mainly carried out by the Ni 3dx2?y2orbital. Our results suggest that the Sr-doped NdNiO2superconductor is essentially a single Ni-3dx2?y2-orbital system, and the only role of the Nd 5d orbital is to slightly tune the doping concentration of the charge carrier. More interestingly, our study can well reproduce the dome-shaped dependence of theTcover the Sr-doping concentration revealed in the recent experiment,[5,6]implying the relevance of our study to the real material.

    The rest of this paper is organized as follows. In Section 2, we describe the adopted TB model. In Section 3, we introduce the details of the RPA approach engaged in our calculations. In Section 4,we provide the results of our numerical calculations:we study both cases of realisticVR,Ni/=0 and imaginaryVR,Ni=0 for comparison. HereVR,Ni(R=La,Nd)represents the hybridization between the Ni-3dx2?y2orbital and theR-5d orbital. It will be seen that the two cases yield very similar results, and that the Ni-3dx2?y2orbital takes the main role in determining the low-energy physics of the system.Finally,in Section 5,conclusions are reached with further discussions.

    2. The model

    TheRNiO2(R=La,Nd)has a body-centered cubic structure with the lattice constantsa=b/=c. In the cubic unit cell,R(Ni) is at the vertex (center) of the cube, and O is in the middle of each face. The band structures of LaNiO2and NdNiO2are very similar, in which only the 3dx2?y2orbital of Ni and the 5d orbitals of La or Nd play the main role near the Fermi energy.[4,11]Here we start from a minimal twoband tight binding model[4]to describe the low-energy physics dominated by the rare-earthRand Ni bands,which writes as

    We take the relevant parameters from Ref.[4],and show the band structure of the Hamiltonian(1)in Fig.1. Note that although the provided TB parameters are for LaNiO2,the obtained band structure can also well describe the band structure of NdNiO2, as the band structures of the two materials are very similar. Here only the Ni 3dx2?y2and La 5d3z2?r2orbitals are considered for simplicity. In this system, although the Ni 3d orbital and the La 5d orbital are coupled through hybridization, such coupling is on the order of O(10?2) eV,which is negligibly small compared to the intra-layer hopping amplitudes. Therefore, each of the two resulting bands has a dominating orbital component. The Ni 3d-orbital dominated band is almost dispersionless alongΓ–Z(thez-axis in the Brillouin zone), and hence it is of two-dimensional (2D)nature. In contrast,the La 5d-orbital dominated band is quite extended and three-dimensional (3D). Specially two small Fermi surface pockets are formed at theΓand theApoints,based on this 3D band. For comparison purpose,we also calculate the band structure with the inter-orbital hybridization turning off, i.e.,Vkk′=0, in which the two bands are completely decoupled. As shown in Fig. 1(b), the resulting band structure is almost identical to that in Fig.1(a),and all features are remained.

    Fig. 1. Band structure of the Hamiltonian in Eq. (1): (a) the realistic case with weak inter-orbital coupling; (b) the comparison case with the interorbital coupling turning off. In(a)((b)),the red one is mainly(completely)composed of R 5d, and the blue one is mainly (completely) composed of Ni 3dx2?y2. (c)The high symmetry points in the tetragonal Brillouin zone.The two dashed lines are the chemical potentials corresponding to doping 0 and doping ?0.2,respectively.

    The density of states (DOS) of the system is shown in Fig. 2(a), with the contributions from the two orbital components marked separately. From Fig.2(a),it is clear that for the undoped parent compound, the most part of the DOS is contributed by the Ni 3dx2?y2orbital,with theR5d orbital slightly contributing to the DOS.Further more,when the compound is hole-doped by substituting someR3+ions with Sr2+atoms,the part contributed by theR5d orbital continues decreasing and vanishes promptly. On the contrary, the part of the DOS contributed by the Ni 3dx2?y2orbital increases promptly upon the hole-doping, and experiences a 3D VHS (introduced below)at the dopingδc=?0.1. After that the DOS decreases,leading to a dome-shaped curve for the DOS contributed by the Ni 3dx2?y2orbital,as shown in Fig.2(a). When the parent compound is electron-doped, the situation becomes very different: the part of the DOS contributed by the Ni 3dx2?y2orbital decreases promptly,while that contributed by theR5d orbital orbital increases. The DOS for the comparison model with the inter-orbital coupling turning off is shown in Fig.2(b),which can not be obviously distinguished from Fig.2(a).

    Fig. 2. DOS of the Hamiltonian in Eq. (1): (a) the realistic case; (b) the comparison case. The contribution to the DOS from each orbital component is also marked,with the red one marking the La 5d component and the blue one marking the Ni 3dx2?y2 one.The grey vertical line marks the VH doping level,which is equal to ?0.1 for(a)and ?0.11 for(b).

    The 3D Fermi surface (FS) of the model is shown in Fig. 3(a) for the undoped parent compound and Fig. 3(b) for the hole-doped case with doping levelδ=?0.2 where theTcof the SC is the highest. For the zero doping, the FS includes a 2D columnar sheet centering around theM–Aline whose main orbital component is the Ni 3d orbital, and two 3D spherical Fermi pockets centering around theΓand theApoints respectively whose main orbital component is theR5d orbital. While the Ni-3d-orbital derived 2D FS column is hole pocket,the twoR-5d-orbital derived pockets are electron pockets. Forδ=?0.2, the two electron pockets contributed mainly from theR-5d-orbital greatly shrink,with theΓpocket nearly invisible. In the meantime,the columnar“hole pocket”formed mainly by the Ni-3d-orbital expands obviously with hole-doping until its top part near thekz=πplane experiences a Lifshitz transition and changes its topology,as is more clear in the 2D FS shown below.

    The 3D FSs for the comparison band structure with the inter-orbital coupling between the Ni 3d andR-5d-orbitals turning off are shown in Fig. 3(c) for the undoped case and Fig. 3(d) forδ=?0.2, which cannot be obviously distinguished from the uncoupled case shown in Figs.3(a)and 3(b).

    Fig. 3. The 3D FSs for the realistic band structure at the doping levels (a)δ =0 and(b)δ =?0.2;and those for the comparison band structure at the doping levels(c)δ =0 and(d)δ =?0.2.

    The 2D FS for the zero doping is shown in Fig. 4(a) forkz= 0 and 4(b) forkz=π, while that for the doping levelδ=?0.2 is shown in Figs. 4(c) and 4(d). Figures 4(a) and 4(b) reveal the two separateR5d-orbital derived 3D electron pockets centering at theΓandApoints,and the Ni 3d-orbital derived 2D columnar FS surrounding theM–Aline which can be regarded as a hole pocket. Clearly, the cut of this Ni 3d pocket on thekz=πplane hosts a larger area than that on thekz=0 plane. When the doping is shifted toδ=?0.2,the twoR5d-orbital derived 3D electron pockets shrink a lot,with theΓ-pocket almost invisible. In the meantime,the Ni 3d pocket expands a lot,as shown by Figs.4(c)and 4(d).Such expansion is more clear on thekz=πplane, where the different sheets of the hole pocket which used to surround theApoint atδ=0 now touch with each other and the FS on this plane has now changed to an electron pocket surrounding theZpoint. Such a FS-topology change on thekz=πplane suggests that the system has experienced a Lifshitz transition,which takes place at the doping levelδc=?0.1. Such a Lifshitz transition makes a van Hove singularity(VHS),with the VH momenta locating at theR-points. Such a VHS on the FS has enhanced the DOS,leading to a broad peak in the DOS curve shown in Fig.2. The 2D FS and its evolution with doping for the comparison band structure with the inter-orbital coupling turning off(un-shown)cannot be obviously distinguished from Fig.4.

    Fig.4. The 2D FSs for the doping levels δ =0 on the kz=0 plane(a)and kz =π plane (b), and those for the doping levels δ =?0.2 on the kz =0 plane(c)and kz=π plane(d). Only the results for the realistic band structure are shown.

    3. The RPA formulation

    Generally, repulsive Hubbard interactions suppress the charge susceptibility,but enhance the spin susceptibility.[34–46]There is a critical interaction strengthUc,where the spin susceptibility diverges, indicating that the long-range magnetic order forms,and the RPA treatment is no longer valid.[47–50,52]When the interaction strength is belowUc, short-ranged spin and charge fluctuations exist in the system. Electrons may form Cooper pairs by exchanging these fluctuations,resulting in exotic superconducting states in the system.

    According to theD4hpoint group[53]of LaNiO2, the possible superconducting pairing symmetries include A1g(s)-,Eu(p)-,B1g(dx2?y2)-,B2g(dxy)-,A2g(dx2?y2×dxy)-,A2u(pz)-,Eg(px×pz/py×pz)-, B2u(dx2?y2×pz)-, B1u(dxy×pz)- and A1u(dx2?y2×dxy×pz)-wave ones. The eigenvector(s)Δα(k)for each eigenvalueλobtained from gap equation(17)as the basis function(s)forms(form)an irreducible representation of theD4hpoint group. While the Eu(p), Eg(px×pz/py×pz)symmetries each form a 2D representation with degenerate pairing eigenvalues,other symmetries each forms a 1D representation of the point group with nondegenerate pairing eigenvalues. The character table for each representation of theD4hpoint group is shown in Table 1.

    Table 1. The D4h character table.

    4. Numerical result

    In this section, we present our numerical results for the spin correlations and pairing symmetries in this system. We compare the results obtained for two cases: one case is for the realistic band structure with weak inter-orbital coupling between the Ni 3d and La 5d orbitals, and the other one is for the imaginary comparison band structure with no coupling between the Ni 3d and La 5d orbitals. Our results suggest that the spin correlations obtained for the two cases are similar,and the pairing symmetries obtained are both dx2?y2which is robust against doping. Our results suggest that the pairing mechanism and pairing symmetry of this system are mainly determined by the 3d orbital of Ni,and the main role of La is only to adjust the effective doping concentration. In the following,we provide our results for both cases simultaneously.For the doping levels, we provide two cases for comparison:one is zero doping, and the other is?0.2 which represents a typical doping level with obviousTcfor the SC. The 3D and 2D FSs for the two doping levels are shown in Figs.3 and 4,respectively.

    In Fig. 5, we show thekdependence of the susceptibilityχ(k) in the Brillouin zone. Hereχrepresents the largest eigenvalue of the susceptibility matrix (χ). The elementχlm(k)of this matrix is defined as

    Theχ(k)for two typical doping levels,i.e.,0 and?0.2,are shown in Fig. 5. For the zero doping, such momentum dependence is shown along the high-symmetry lines in the Brillouin zone in (a), on thekz=0 plane in (b), and on thekz=πplane in (c). Clearly, there are two peaks for the distribution ofχ(q), with one near theM=(π,π,0) point, and the other near theA=(π,π,π) point. The in-plane component of both momenta is(π,π),reflecting the strong in-plane antiferromagnetic(AFM)spin correlations. The two peaks are nearly degenerate, reflecting the competition between interlayer AFM and FM correlations. When the doping concentration is changed to?0.2,such momentum dependence is shown in Figs.5(d)–5(f). In this doping level,there are still two susceptibility peaks near theMandApoints. However, theApeak is slightly higher than theM-peak, reflecting that while the in-plane spin correlation is still AFM, the inter-layer one is also AFM.Similarkdependence ofχ(k)for the comparison band structure with the inter-orbital coupling turning off is shown in Fig.6. A comparison between Figs.5 and 6 reveals that the two figures are nearly the same,suggesting that the influence of the inter-orbital coupling between the Ni 3d orbital and theR5d orbital on the spin correlations in this system can be ignored.

    Fig.5. The k-space distribution of the largest eigenvalue χ(k)of the susceptibility matrix for the realistic band structure.Of which,(a),(b),(c)represent the doping 0 and(d),(e),(f)represent the doping ?0.2. (a)and(d)are along the high-symmetry lines in the Brillouin zone, (b)and(e)are on the kz =0 plane,and(c)and(f)are on the kz=π plane.

    Fig.6. The k-space distribution of the largest eigenvalue χ(k)of the susceptibility matrix for the comparison band structure with the inter-orbital coupling turned off. The meaning of each figure is similar to that in Fig.5.

    From the eigenvectorξ(km)corresponding to the largest eigenvalueχ(km) of the susceptibility matrix atkm, the spin fluctuation in the system takes the pattern ofsiμ=ξμ(km)eikm·ri. The spin-fluctuation pattern thus obtained for the doping levelδ=?0.2 is shown from top view in Fig.7(a)and from side view in Fig. 7(b). As thekmfor this doping level is atA,the obtained spin-fluctuation pattern is the Neelordered AFM pattern both intra-layer and inter-layer. What’s more, the distribution ofsiμbetween the two orbitalsμ=1 andμ=2 is determined byξμ(km).From our calculation,it is obtained asξ1(km)/ξ2(km)=10?4, which is very small. As a result,the spin fluctuation is mainly carried by the Ni 3d orbital,as shown in Figs.7(a)and 7(b). Similar spin-fluctuation patterns thus obtained for the comparison case with the coupling between the Ni 3d and theR5d orbitals turning off are shown in Figs.7(c)and 7(d)for the top and side views,respectively. In Figs. 7(c) and 7(d), the spin pattern is completely distributed on the Ni 3d orbital asξ1(km)=0,ξ2(km)=1,suggesting that if we omitted the inter-orbital coupling the spin fluctuation should be completely carried by the Ni 3d orbital. Comparing the case with realistic weak inter-orbital coupling between the Ni 3d andR5d orbitals and the one without such inter-orbital coupling, it is known that the spin fluctuation in this system is mainly carried by the Ni 3d orbital and theR5d orbital only participates into the spin fluctuation through the tiny hybridization with Ni 3d orbital.

    Fig. 7. The spin-fluctuation pattern, in which blue color represents Ni and the red one represents La. (a),(b)are top view and side view for the realistic band structure with weak inter-orbital coupling and(c),(d)are top view and side view for the comparison band structure with the inter-orbital coupling turned off. The doping for both cases is ?0.2. The relative ratio between the lengths of the arrows representing the spins of the two orbitals has been exaggerated.

    When the interaction strength is below the critical one,the above illustrated spin fluctuation is short-ranged. Through exchanging these short-ranged spin fluctuations,a pair of electrons with opposite momenta can acquire effective attraction and pair, which leads to the SC. In our calculations, we takeU=0.9 eV,under which the renormalized spin and charge susceptibilities do not diverge,which suggest that this interaction strength is below the critical one. The pairing phase diagram thus obtained is shown in Fig. 8. Figure 8 suggests that the phase diagram (a) for the realistic band structure with weak inter-orbital coupling between Ni 3d andR5d and that(b)for the comparison band structure with the inter-orbital coupling turning off are very similar. Figure 8 reveals three experimentally relevant characters in the doping dependence of the pairing phase diagram.

    Fig. 8. Doping δ dependence of the largest pairing eigenvalue λ for each pairing symmetry for the realistic band structure (a) and the comparison band structure(b).

    Firstly, the pairingTcfor the hole-doped case is much larger than that for the electron-doped case. Particularly, theTcfor the electron-doped case decreases promptly with doping, which might possibly be invisible in experiment. This character is different from that of the cuprates superconductors which host high-temperature SC for both hole-doped and electron-doped cases and is consistent with the experiment of the Sr-doped NdNiO2superconductors wherein the Sr doping introduced holes into the material. The physical reason in our study for this character lies in that,on the one hand the Ni 3d orbital-component contributes less and less to the DOS of the electron-doped case,and on the other hand it is the Ni 3d orbital-component which contributes to the SC of the material,as more clearly illustrated below. Our study thus predicts that no obvious SC would be detected for the electron-doped case of the NdNiO2.

    Secondly, the pairingTcfor the hole-doped case exhibits a dome shape, which peaks in between the doping levels?0.1 and?0.2. Such a dome-shaped curve of the doping-dependentTcis qualitatively consistent with recent experiments[5,6]wherein the peak of theTcis near the doping?0.2. The physical reason of this character lies in that,the Ni 3d-orbital component of the DOS, which contributes mainly to the SC of the system,peaks near this doping range,as shown in the DOS curve in Fig. 2. This DOS peak originates from the VHS caused by the Lifshitz transition of the FS,as shown in Fig.9. Figures 9(a)and 9(b)show the 2D FS on thekz=0 andkz=πplanes for the VH dopingδc=?0.1 for the realistic band structure. This 2D FS can be compared with those forδ=0 andδ=?0.2 shown in Fig. 4.Clearly, while the 2D FSs for the three cases on thekz=0 plane have the same topology,those on thekz=πplane illustrate a FS-topology change. While the FS on thekz=πplane forδ=0 shown in Fig.4(b)comprises aR5d-orbital derived electron pocket and a Ni 3d-orbital derived hole pocket centering around theApoint, that forδ=?0.2 shown in Fig. 4(d)comprises two electron pockets,with the Ni 3d-orbital derived FS sheet now changing topology to an electron pocket centering around theZpoint. Clearly, the Ni 3d-orbital derived FS sheet has experienced a Lifshitz transition fromδ=0 toδ=?0.2. Figure 9(b) just illustrates the Lifshitz transition point atδc=?0.1,on which the separate FS sheets appearing in Fig. 4(b) now touch with each other at theR-points, making these points as the VH momenta. This 3D VHS will not lead to a divergence of the DOS, but it causes a broad peak in the DOS in Fig.2(a),which just leads to the broadTc-peak in Fig. 8(a). The situation for the comparison band structure(δc=?0.11)shown in Figs.9(c)and 9(d)is similar,which is related to the broadTc-peak in Fig.8(b).

    Fig.9.The 2D FS of the VH doping(δc=?0.1)for the realistic band structure on the kz=0 plane(a)and kz=π plane(b),and that(δc=?0.11)for the comparison band structure on the kz=0 plane(c)and kz=π plane(d).

    Thirdly, while all the ten possible pairing symmetries listed in Table 1 are possible, the dx2?y2-wave pairing dominates over other pairing symmetries in the whole doping range where obvious SC can be detected. This robust dx2?y2-wave pairing obtained for our simple model is consistent with those obtained for more complicated models.[11,22]The distribution of the gap function for the obtained dx2?y2-wave pairing on the FS is shown in Fig. 10(a) on thekz=0 plane and Fig. 10(b)on thekz=πplane for the doping levelδ=?0.2. Obviously, this gap function is symmetric about thexandyaxes and is antisymmetric about thex=±yaxes. Four symmetryprotected line nodes are presented in the 3D Brillouin zone whose projections on thexyplane locate on thex=±yaxes.It is also obvious that the pairing gap is mainly distributed on the Ni 3d-orbital derived FS sheets,with the component on theR5d-orbital derived FS sheets almost vanished,suggesting the dominating role of the Ni 3d orbital in determining the SC in this system. The situation for the comparison band structure is similar,as shown in Figs.10(c)and 10(d).

    Fig. 10. Distributions of the gap functions on the FS for the doping level?0.2 for the realistic band structure(a),(b)and the comparison band structure with the inter-orbital coupling turned off(c),(d). While(a)and(c)are the distribution on the kz=0 plane,(b)and(d)represent that on the kz=π plane.

    5. Discussion and conclusion

    Due to the strong on-site Coulomb interaction between the Ni 3d electrons, the nickelate superconductors might belong to strongly-correlated system, or at least intermediatelycorrelated one. Here, to avoid the difficulty in the systematic treatment of strongly- or intermediately-correlated electronic systems,we start from the weak-coupling RPA approach,with the adopted Hubbard-U(0.9 eV)several times weaker than the realistic value. Two reasons lie behind this choice ofU:[34–43]on the one hand, the RPA as a perturbational approach does not apply to the case with strongU; on the other hand, the conditionU

    Note that the weak-coupling RPA approach adopted here might have missed the Kondo physics at the small doping levels, where the Kondo coupling between the Ni-3d and theR-5d orbitals might lead to the d+is-wave pairing.[13]However,when the hole-doping levels are large, due to the electronpocket character of theR-5d dominant FSs, the DOS contributed by theR-5d orbital component is largely reduced,and correspondingly the Kondo physics should be suppressed.Therefore, the conclusion arrived at in our study should only be tenable for large hole-doping concentrations.

    In conclusion, we have performed a RPA based calculation for the Sr-doped RNiO2superconductor,adopting a realistic two-band band structure. Through comparison between the results obtained for the realistic band structure and those obtained for the imaginary band structure with the inter-orbital coupling between the Ni 3d andR5d orbitals turned off,it is revealed that the Ni 3d orbital dominates over theR5d orbital in carrying out the AFM spin fluctuation and the superconducting pairing. Our results suggest that this system is essentially a one-orbital(i.e.,Ni 3d)system,with the only role of theR5d orbital to slightly tuning the doping concentration in the hole-doped superconducting regime relevant to experiments.Therefore,this material can be well captured by a single-band Hubbard model on the square lattice.

    Further more,our results reveal three experimentally relevant doping-dependent behaviors of the pairing state. Firstly,the superconducting pairing is strongly particle–hole asymmetric,with obviousTconly detectable on the hole-doped side.Secondly, a dome-shaped doping dependence of the SCTcis obtained, which peaks within the regime of (?0.2,?0.1),qualitatively consistent with recent experiments. These two characters are determined by the doping dependence of the DOS contributed by the Ni 3d orbital, which is particle–hole asymmetric and peaks near the doping regime of(?0.2,?0.1),caused by the 3D VHS.Thirdly, the robust dx2?y2-wave pairing symmetry dominates over other pairing symmetries in the large doping regime,consistent with the results obtained from other more complicated models.

    Acknowledgement

    We are grateful to the stimulating discussions with Congjun Wu,Chen Lu,Lun-Hui Hu and Zhan Wang.

    Fig.A1. The hoping between particles in real space,and the unit of related parameters is eV.

    Appendix B: The effective interaction

    The Feynman’s diagram of RPA is shown in Fig. B1.Note that when the interaction strengthUis higher thanUc,the denominator matrixI ?χ(0)(k,iwn)Usin Eq.(8)will have zero eigenvalues for someqand the renormalized spin susceptibility diverges there, which invalidates the RPA treatment.WhenU

    Fig.B1.Feynman’s diagram for the renormalized susceptibilities in the RPA level.

    Fig. B2. Three processes which contribute the renormalized effective vertex considered in the RPA, with (a) the bare interaction vertex and (b), (c)the two second order perturbation processes during which spin or charge fluctuations are exchanged between a Cooper pair.

    This equation can be looked upon as an eigenvalue problem, where the normalized eigenvectorδα(k) represents the relative gap function on theα-th FS patches nearTc, and the eigenvalueλdeterminesTcviaTc=cutoff energy×e?1/λ,where the cutoff energy is of order bandwidth. The leading pairing symmetry is determined by the largest eigenvalueλof Eq.(B6).

    猜你喜歡
    楊帆
    Band structures of strained kagome lattices
    Effect of short-term plasticity on working memory
    《魚(yú)與蓮》
    Analytical formula describing the non-saturating linear magnetoresistance in inhomogeneous conductors
    Two regularization methods for identifying the source term problem on the time-fractional diffusion equation with a hyper-Bessel operator
    Three-dimensional spatial multi-point uniform light focusing through scattering media based on feedback wavefront shaping?
    劫后華夏再楊帆(弋陽(yáng)腔)
    影劇新作(2020年2期)2020-09-23 03:22:12
    THE QUASI-BOUNDARY VALUE METHOD FOR IDENTIFYING THE INITIAL VALUE OF THE SPACE-TIME FRACTIONAL DIFFUSION EQUATION ?
    假期(劇本)
    我學(xué)會(huì)了
    精品一区二区免费观看| 午夜激情福利司机影院| 狂野欧美激情性bbbbbb| 黄色一级大片看看| 国产大屁股一区二区在线视频| 欧美3d第一页| 日韩电影二区| 国产欧美日韩精品一区二区| 国产午夜精品一二区理论片| 国产黄片美女视频| 日韩av不卡免费在线播放| 日韩一区二区视频免费看| 国产毛片a区久久久久| 99精国产麻豆久久婷婷| 最近中文字幕2019免费版| 两个人的视频大全免费| 国产 一区精品| 97人妻精品一区二区三区麻豆| 亚洲精品一二三| 一区二区三区四区激情视频| 日本黄大片高清| 中文字幕久久专区| 特级一级黄色大片| 亚洲天堂av无毛| 少妇人妻久久综合中文| 丝袜美腿在线中文| 亚洲精品日本国产第一区| 少妇裸体淫交视频免费看高清| av在线蜜桃| 亚洲精品日本国产第一区| 丰满乱子伦码专区| 日日摸夜夜添夜夜爱| 国产人妻一区二区三区在| 久久99热6这里只有精品| 成人鲁丝片一二三区免费| 另类亚洲欧美激情| 国产有黄有色有爽视频| 丰满人妻一区二区三区视频av| 久久久久久久午夜电影| 性色avwww在线观看| 少妇人妻精品综合一区二区| 久久久久久久国产电影| a级一级毛片免费在线观看| 国产男女超爽视频在线观看| 午夜日本视频在线| av免费在线看不卡| av国产久精品久网站免费入址| 亚洲国产精品成人久久小说| 亚洲精品456在线播放app| 亚洲精华国产精华液的使用体验| 网址你懂的国产日韩在线| 大话2 男鬼变身卡| 国产一区二区三区综合在线观看 | 一本色道久久久久久精品综合| 久久人人爽av亚洲精品天堂 | 毛片一级片免费看久久久久| 国产亚洲av嫩草精品影院| 99热6这里只有精品| 最近手机中文字幕大全| 国产精品伦人一区二区| 午夜精品一区二区三区免费看| 欧美+日韩+精品| 欧美区成人在线视频| 成人免费观看视频高清| 亚洲精品中文字幕在线视频 | 亚洲国产欧美在线一区| 黄色怎么调成土黄色| 亚洲最大成人手机在线| 久久精品熟女亚洲av麻豆精品| 亚洲在久久综合| 成人亚洲精品一区在线观看 | 亚洲欧美成人精品一区二区| 十八禁网站网址无遮挡 | 免费看av在线观看网站| 久久久久久久久久久丰满| 免费大片黄手机在线观看| 寂寞人妻少妇视频99o| 97超视频在线观看视频| 欧美激情在线99| 一级毛片我不卡| 97在线人人人人妻| 麻豆成人av视频| 国产精品福利在线免费观看| 国产成年人精品一区二区| 国产精品一及| 亚洲国产精品成人久久小说| 老师上课跳d突然被开到最大视频| 国产精品.久久久| 国内少妇人妻偷人精品xxx网站| 久久精品久久久久久噜噜老黄| 听说在线观看完整版免费高清| 欧美日韩国产mv在线观看视频 | 99热这里只有是精品在线观看| 亚洲精品乱码久久久v下载方式| 亚洲经典国产精华液单| 一区二区三区精品91| 一级毛片电影观看| 国产老妇伦熟女老妇高清| 美女cb高潮喷水在线观看| 国产乱人偷精品视频| 极品教师在线视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产成人一区二区在线| 久久人人爽人人爽人人片va| 午夜福利网站1000一区二区三区| 久久鲁丝午夜福利片| 亚洲av成人精品一区久久| 国产午夜精品一二区理论片| 亚洲色图av天堂| 国产精品av视频在线免费观看| 亚洲最大成人av| 少妇熟女欧美另类| 国产淫语在线视频| 国产精品福利在线免费观看| 97在线视频观看| 丰满乱子伦码专区| 性插视频无遮挡在线免费观看| 国产成年人精品一区二区| 中国美白少妇内射xxxbb| 久久久精品94久久精品| 久久久久久国产a免费观看| xxx大片免费视频| 五月玫瑰六月丁香| 99热全是精品| 久久这里有精品视频免费| 久久久午夜欧美精品| 99久久中文字幕三级久久日本| 97热精品久久久久久| 午夜老司机福利剧场| 91久久精品国产一区二区三区| 成人毛片a级毛片在线播放| 久久久久久久亚洲中文字幕| 中文字幕亚洲精品专区| 久久久精品免费免费高清| 亚洲av日韩在线播放| 丰满人妻一区二区三区视频av| 丝袜脚勾引网站| 又黄又爽又刺激的免费视频.| 成人欧美大片| 精品熟女少妇av免费看| 在线观看三级黄色| 又爽又黄无遮挡网站| 欧美变态另类bdsm刘玥| 国产免费一级a男人的天堂| 搡女人真爽免费视频火全软件| 欧美另类一区| 国产老妇伦熟女老妇高清| 涩涩av久久男人的天堂| 日韩人妻高清精品专区| 狂野欧美激情性xxxx在线观看| 亚洲精品影视一区二区三区av| 简卡轻食公司| 午夜免费观看性视频| 大香蕉久久网| av又黄又爽大尺度在线免费看| 丰满少妇做爰视频| 18+在线观看网站| 国产爱豆传媒在线观看| 嫩草影院入口| 亚洲av在线观看美女高潮| 纵有疾风起免费观看全集完整版| 777米奇影视久久| 舔av片在线| 亚洲不卡免费看| 国产探花极品一区二区| 97超视频在线观看视频| 日韩欧美精品免费久久| 精品久久久久久久人妻蜜臀av| av.在线天堂| 五月伊人婷婷丁香| 欧美日韩国产mv在线观看视频 | 少妇 在线观看| 日韩一区二区视频免费看| 波野结衣二区三区在线| 精品一区在线观看国产| 99热全是精品| 亚洲电影在线观看av| 边亲边吃奶的免费视频| 亚洲欧美精品自产自拍| 国产成人精品婷婷| 欧美3d第一页| www.av在线官网国产| 自拍偷自拍亚洲精品老妇| 色吧在线观看| 亚洲av免费高清在线观看| 亚洲国产精品国产精品| av网站免费在线观看视频| 777米奇影视久久| 久久久精品94久久精品| 99热这里只有精品一区| 大香蕉久久网| 国产成人精品婷婷| 亚洲人与动物交配视频| 亚洲色图av天堂| 久久久久精品久久久久真实原创| 深爱激情五月婷婷| av在线app专区| 九草在线视频观看| 久久精品综合一区二区三区| 99热6这里只有精品| 人人妻人人澡人人爽人人夜夜| 精品国产乱码久久久久久小说| 久久久久九九精品影院| 欧美一级a爱片免费观看看| 久久女婷五月综合色啪小说 | 我要看日韩黄色一级片| 99热这里只有是精品在线观看| 中文字幕久久专区| 国产精品成人在线| 国产成人a∨麻豆精品| 亚洲精品一区蜜桃| 欧美三级亚洲精品| 亚洲精品乱久久久久久| 亚洲国产成人一精品久久久| 精品酒店卫生间| 欧美精品国产亚洲| 精品一区二区三区视频在线| www.av在线官网国产| 男女下面进入的视频免费午夜| 国产精品熟女久久久久浪| 免费人成在线观看视频色| 永久免费av网站大全| 久久6这里有精品| 香蕉精品网在线| 在线天堂最新版资源| 欧美日韩视频精品一区| 国产永久视频网站| 日韩欧美精品免费久久| 国产色婷婷99| 亚洲成人中文字幕在线播放| 六月丁香七月| 亚洲av中文字字幕乱码综合| 亚洲精品日韩在线中文字幕| 欧美少妇被猛烈插入视频| 九九爱精品视频在线观看| 亚洲经典国产精华液单| 日韩成人av中文字幕在线观看| 干丝袜人妻中文字幕| 欧美区成人在线视频| av免费在线看不卡| 中文字幕制服av| 国产老妇女一区| 久久精品国产自在天天线| 亚洲欧美日韩卡通动漫| 一级爰片在线观看| 亚洲欧美成人综合另类久久久| 欧美成人一区二区免费高清观看| 国产伦理片在线播放av一区| 成人黄色视频免费在线看| 99久国产av精品国产电影| 成人欧美大片| 久久精品久久久久久噜噜老黄| 最后的刺客免费高清国语| 成人亚洲精品av一区二区| 边亲边吃奶的免费视频| 黄色怎么调成土黄色| av专区在线播放| 蜜臀久久99精品久久宅男| 国语对白做爰xxxⅹ性视频网站| 在线观看人妻少妇| 欧美日韩在线观看h| 亚洲国产精品成人综合色| 国产精品嫩草影院av在线观看| 女人久久www免费人成看片| 亚洲av成人精品一二三区| 国产精品一二三区在线看| 特级一级黄色大片| 欧美少妇被猛烈插入视频| 国产亚洲一区二区精品| 午夜福利高清视频| 精品少妇久久久久久888优播| 亚洲三级黄色毛片| 自拍偷自拍亚洲精品老妇| 在现免费观看毛片| 国产乱人偷精品视频| 亚洲第一区二区三区不卡| 91在线精品国自产拍蜜月| av线在线观看网站| 久久精品国产自在天天线| 亚洲国产最新在线播放| 一个人看视频在线观看www免费| 国产精品av视频在线免费观看| 久久精品久久久久久久性| 久久女婷五月综合色啪小说 | 国产一区二区亚洲精品在线观看| 国产黄色视频一区二区在线观看| 免费观看av网站的网址| 我要看日韩黄色一级片| 伊人久久国产一区二区| 国产伦精品一区二区三区视频9| 日韩av在线免费看完整版不卡| 少妇裸体淫交视频免费看高清| 国产爽快片一区二区三区| 2022亚洲国产成人精品| 97热精品久久久久久| 在线精品无人区一区二区三 | 看非洲黑人一级黄片| 午夜爱爱视频在线播放| 日韩成人伦理影院| 亚洲精品日韩在线中文字幕| 亚洲国产精品成人久久小说| 日本熟妇午夜| 欧美性感艳星| av黄色大香蕉| 成人国产麻豆网| 国产亚洲av嫩草精品影院| 看十八女毛片水多多多| 精品一区在线观看国产| 久久精品国产a三级三级三级| 男女啪啪激烈高潮av片| 听说在线观看完整版免费高清| 久久精品熟女亚洲av麻豆精品| 神马国产精品三级电影在线观看| 久久精品夜色国产| 精品酒店卫生间| 视频中文字幕在线观看| 国产欧美日韩一区二区三区在线 | 三级国产精品欧美在线观看| 欧美成人一区二区免费高清观看| 亚洲三级黄色毛片| 自拍欧美九色日韩亚洲蝌蚪91 | 久久热精品热| 丝袜美腿在线中文| 国产白丝娇喘喷水9色精品| 中文天堂在线官网| 午夜老司机福利剧场| 免费人成在线观看视频色| 精品久久久久久久久亚洲| 各种免费的搞黄视频| 日韩国内少妇激情av| 看十八女毛片水多多多| 插逼视频在线观看| www.色视频.com| 男插女下体视频免费在线播放| 又黄又爽又刺激的免费视频.| 成人亚洲欧美一区二区av| 日本免费在线观看一区| 精品久久久久久久久av| av.在线天堂| 亚洲色图av天堂| 男人添女人高潮全过程视频| 熟妇人妻不卡中文字幕| 亚洲精品中文字幕在线视频 | 亚洲精品中文字幕在线视频 | 欧美成人午夜免费资源| 九九爱精品视频在线观看| 岛国毛片在线播放| 国产精品成人在线| 欧美另类一区| 亚洲国产精品成人综合色| 国产淫语在线视频| 日产精品乱码卡一卡2卡三| 亚洲精品久久午夜乱码| 国产精品精品国产色婷婷| 一本一本综合久久| 久久久久久久国产电影| 看黄色毛片网站| 在线精品无人区一区二区三 | 嘟嘟电影网在线观看| 亚洲在线观看片| 国产精品99久久久久久久久| 日本黄大片高清| 亚洲内射少妇av| 国产精品久久久久久久久免| 免费大片黄手机在线观看| av天堂中文字幕网| 最新中文字幕久久久久| 校园人妻丝袜中文字幕| 在线 av 中文字幕| 亚洲精品国产av蜜桃| 夜夜爽夜夜爽视频| www.av在线官网国产| 国产精品一及| 美女内射精品一级片tv| 欧美丝袜亚洲另类| 国产在线一区二区三区精| 五月开心婷婷网| 午夜福利在线在线| 九色成人免费人妻av| 综合色丁香网| 国产一区二区亚洲精品在线观看| 在线观看美女被高潮喷水网站| 日本爱情动作片www.在线观看| av一本久久久久| 男人爽女人下面视频在线观看| 成人午夜精彩视频在线观看| av在线蜜桃| 国产有黄有色有爽视频| 我的老师免费观看完整版| 一级二级三级毛片免费看| 国产在视频线精品| 亚洲天堂国产精品一区在线| 特级一级黄色大片| 国产成人精品婷婷| 丝袜脚勾引网站| 赤兔流量卡办理| 高清日韩中文字幕在线| 精品久久久久久久末码| 乱码一卡2卡4卡精品| 内射极品少妇av片p| 一级片'在线观看视频| 亚洲av成人精品一区久久| xxx大片免费视频| 男人添女人高潮全过程视频| 97在线视频观看| 麻豆精品久久久久久蜜桃| freevideosex欧美| 日韩免费高清中文字幕av| 人妻制服诱惑在线中文字幕| 精品久久久久久久久av| 涩涩av久久男人的天堂| 亚洲成人精品中文字幕电影| 亚洲国产色片| 婷婷色麻豆天堂久久| 性插视频无遮挡在线免费观看| 亚洲成人av在线免费| 国产亚洲最大av| 最近最新中文字幕大全电影3| 22中文网久久字幕| 免费观看的影片在线观看| 欧美日韩综合久久久久久| 亚洲色图综合在线观看| 小蜜桃在线观看免费完整版高清| 国产亚洲午夜精品一区二区久久 | 交换朋友夫妻互换小说| 日韩免费高清中文字幕av| 我的老师免费观看完整版| 18禁在线无遮挡免费观看视频| 色网站视频免费| 一级二级三级毛片免费看| 亚洲国产日韩一区二区| 成人漫画全彩无遮挡| 内射极品少妇av片p| 国产一区二区亚洲精品在线观看| 麻豆成人av视频| 国产av不卡久久| 国产精品人妻久久久久久| 中文精品一卡2卡3卡4更新| 男女国产视频网站| 亚洲精品一区蜜桃| 欧美bdsm另类| 久久精品国产鲁丝片午夜精品| 欧美亚洲 丝袜 人妻 在线| 国产v大片淫在线免费观看| 亚洲真实伦在线观看| 亚洲精品国产色婷婷电影| 中文在线观看免费www的网站| 亚洲精品成人久久久久久| 久久久久久久久久久丰满| 少妇人妻一区二区三区视频| 又爽又黄a免费视频| 久久精品国产自在天天线| 乱码一卡2卡4卡精品| 国产一区二区在线观看日韩| 真实男女啪啪啪动态图| 亚洲综合色惰| 国产在线男女| 啦啦啦在线观看免费高清www| 午夜视频国产福利| 国产国拍精品亚洲av在线观看| 亚洲精品一二三| 国产黄a三级三级三级人| 国产色婷婷99| 亚洲av.av天堂| 免费黄网站久久成人精品| 国产精品99久久久久久久久| 丰满乱子伦码专区| 中文资源天堂在线| 国产亚洲最大av| 黄色配什么色好看| 亚洲最大成人av| 色综合色国产| 国产亚洲最大av| 国产精品一区二区三区四区免费观看| 国产亚洲午夜精品一区二区久久 | 五月玫瑰六月丁香| kizo精华| 一级a做视频免费观看| 内射极品少妇av片p| 亚洲最大成人av| 国产精品.久久久| 天美传媒精品一区二区| 色网站视频免费| 亚洲色图综合在线观看| 国产免费福利视频在线观看| 97热精品久久久久久| 国产永久视频网站| 人妻少妇偷人精品九色| 一区二区三区精品91| 亚洲人成网站在线播| 精品国产乱码久久久久久小说| 免费黄频网站在线观看国产| 免费不卡的大黄色大毛片视频在线观看| 少妇的逼水好多| 大片电影免费在线观看免费| 十八禁网站网址无遮挡 | 国产精品国产三级专区第一集| 18禁裸乳无遮挡免费网站照片| 成人毛片60女人毛片免费| 日韩欧美精品免费久久| 国产黄色免费在线视频| 高清欧美精品videossex| 成年女人看的毛片在线观看| 久久久久性生活片| 久久99蜜桃精品久久| 极品少妇高潮喷水抽搐| 成人毛片60女人毛片免费| 我要看日韩黄色一级片| 日韩在线高清观看一区二区三区| 国产成年人精品一区二区| 尾随美女入室| 噜噜噜噜噜久久久久久91| 国产精品一区二区在线观看99| 在线观看免费高清a一片| 边亲边吃奶的免费视频| 亚洲国产精品成人久久小说| 美女主播在线视频| 建设人人有责人人尽责人人享有的 | 亚洲av不卡在线观看| 三级男女做爰猛烈吃奶摸视频| 亚洲aⅴ乱码一区二区在线播放| 在线看a的网站| 老司机影院成人| 中文欧美无线码| 国产乱来视频区| 草草在线视频免费看| 51国产日韩欧美| 婷婷色麻豆天堂久久| 亚洲va在线va天堂va国产| videos熟女内射| 九九久久精品国产亚洲av麻豆| 一本一本综合久久| 五月开心婷婷网| 国产精品久久久久久av不卡| 可以在线观看毛片的网站| eeuss影院久久| 亚洲av不卡在线观看| 国产黄频视频在线观看| 秋霞在线观看毛片| 亚洲精品视频女| 水蜜桃什么品种好| 纵有疾风起免费观看全集完整版| 热99国产精品久久久久久7| 99视频精品全部免费 在线| 97热精品久久久久久| 国产男女超爽视频在线观看| 六月丁香七月| 亚洲经典国产精华液单| 久久韩国三级中文字幕| 黄色日韩在线| 综合色丁香网| 高清视频免费观看一区二区| 欧美一区二区亚洲| 色综合色国产| 97热精品久久久久久| 亚洲成色77777| 欧美性猛交╳xxx乱大交人| 日韩成人伦理影院| 色播亚洲综合网| 人妻制服诱惑在线中文字幕| 亚洲内射少妇av| 天天一区二区日本电影三级| 超碰av人人做人人爽久久| 九九在线视频观看精品| 大片免费播放器 马上看| 80岁老熟妇乱子伦牲交| 18禁在线无遮挡免费观看视频| 联通29元200g的流量卡| 少妇裸体淫交视频免费看高清| 九九在线视频观看精品| 欧美日韩视频精品一区| 99视频精品全部免费 在线| 99久久精品一区二区三区| 在线精品无人区一区二区三 | 在线观看免费高清a一片| 免费在线观看成人毛片| 亚洲性久久影院| 日韩一本色道免费dvd| 亚洲综合色惰| 99久国产av精品国产电影| 人人妻人人爽人人添夜夜欢视频 | 99精国产麻豆久久婷婷| 18禁在线播放成人免费| 精品国产一区二区三区久久久樱花 | 欧美高清成人免费视频www| 国产欧美另类精品又又久久亚洲欧美| 亚洲av不卡在线观看| 亚洲电影在线观看av| 在线免费十八禁| 久久久久久久久久人人人人人人| 免费观看性生交大片5| 欧美丝袜亚洲另类| 夫妻性生交免费视频一级片| 91在线精品国自产拍蜜月| 国产v大片淫在线免费观看| 久久99热这里只有精品18| 亚洲av不卡在线观看| 2021天堂中文幕一二区在线观| 精品久久久精品久久久| 国产精品国产三级专区第一集| 韩国av在线不卡| 日韩人妻高清精品专区| 亚洲av不卡在线观看| 99热国产这里只有精品6| 99热6这里只有精品| 精品久久久久久久久亚洲| 熟妇人妻不卡中文字幕| 国产黄频视频在线观看| 熟妇人妻不卡中文字幕| 美女国产视频在线观看| 蜜桃久久精品国产亚洲av| 精品久久国产蜜桃| 国产成人一区二区在线| 深爱激情五月婷婷| 国产免费福利视频在线观看|