# Distinguishing the origin of the superconducting state from the pseudogap of high-temperature superconductors

## Abstract

We consider an electronic phase separation process that generates regions of different charge densities, or local dopings, as the origin of the inhomogeneous charge density of high superconductors. We show that it gives rise to a phase boundary potential between such doping disordered regions or grains. The Bogliubov-deGennes self-consistent calculations in this disordered medium yield position dependent superconducting gaps which are, for all dopings, smaller than those derived from the local density of states with a pseudogap behavior. Studying these two sets of gaps for different temperatures and dopings, we are able to reproduce many many observed properties of superconducting cuprates. This scenario is consistent with a resistivity transition driven by Josephson coupling among the superconducting grains.

###### pacs:

74.20.-z, 74.25.Dw, 74.62.En, 74.72.KfThe origin of the superconducting gap associated with the superconducting state and the relation to the pseudogap above the transition temperature () remains one of the central questions in high- research. In this letter, we show that the solution of this important problem is connected to an electronic granular structure derived from an electronic phase separation (EPS) transition.

There are steadily accumulating evidences that the charge distribution is microscopically inhomogeneous in the planes of high temperature superconductors (HTSC)(1); (2); (3); (4); (5); (6); (7); (8). Low temperature Scanning Tunneling Microscopy (STM) has demonstrated non-uniform energy gaps that vary on the length scale of nm(2) over the whole surface of these materials. These gaps have two types of shapes(3), and some remain well above the superconducting critical temperature (4); (5). More recently, the existence of two types of energy gaps have been observed on electronic Raman scattering experiments(9), STM data(6); (7); (8), Angle Resolved Photon Emission (ARPES)(10) and combined STM-ARPES(11). However, the origin and even the existence of this two-gap picture is still a matter of debate(12); (13).

A possible reason to this complex behavior is an EPS transition driven by the minimization of the free energy(14). This minimization is due to the formation of anti-ferromagnetic (AF) regions of almost zero local doping () with lower free energies than the regions with large local doping level (), where is the average hole doping of the sample(14). As increases, the Coulomb repulsion in the large local doping regions generates a high energy cost to the phase separation process and the EPS ceases in the overdoped region, in agreement with the disappearance of the local AF fluctuations (15). This cooperative phenomenon gives rise to (grain) boundary potentials in the planes, single particle bound states and intra-grain superconductivity at low temperature. These superconducting regions develop Josephson couplings among them and the resistivity transition occurs when the Josephson energy is equal to .

We used the theory of Cahn-Hilliard (CH)(16); (17); (18) to binary alloys to describe the EPS transition. The transition order parameter is the difference between the local and the average charge density . In Fig.(1), we show a typical density map with the two (hole-rich and hole-poor) main solutions given by different colors.

In the inset of Fig.(1) we plot the usual potential from the Ginzburg-Landau (GL) free energy expansion used in the CH equation(14) with its two minima and the potential barrier between them. The barriers generate shallow potential wells in the planes that always have an exponential small bound state(19). These bound states lower the kinetic energy to allow, under favorable conditions, local hole pair formation in the antiadiabatic limit(20).

The calculated (or ) density map, as shown in Fig.(1), is used as the initial input and it is maintained fixed throughout the self-consistent Bogoliubov-deGennes (BdG) calculations. Notice that the charge inhomogeneity and are correlated. In what follows, the parameters involved in are chosen to match the average local density of states (LDOS) gaps measured by low temperature STM on Bi2212 compounds(3). All others parameters are similar to values previously used(21); (22).

The effect of the temperature in the potential is taken into account using , as demonstrated by Cahn and Hilliard(16). Thus we were able to obtain the intra-grain or local superconducting temperatures , i.e., the onset temperature for the d-wave superconducting gap at a given location . The largest value of all determines the temperature that marks the onset of local superconductivity of the sample.

Similar to granular superconductors(23), the superconducting transition occurs in two steps: first by the appearing of intra-grain superconductivity and than by Josephson coupling with phase locking at a lower temperature. This approach provides a clear interpretation to the superconducting amplitude and the measured quasiparticles dispersion above (12); (13).

By using the theory of granular superconductors(24) to these electronic grains and the calculated average superconducting amplitudes , where N is the total number of sites, we can estimate the values of .

(1) |

Where is the average of the local superconducting gaps on a (, and ) square lattice which are plotted in the top panel of Fig.(2). The inset shows the schematically single particle levels at a shallow puddle whose walls are proportional to . is the normal resistance of a given compound, which is proportional to the planar resistivity measurements(25) on the series. In the low panel of Fig.(2), the Josephson coupling is plotted together with the thermal energy whose intersection yields the critical temperature , as shown in the inset. The values are in reasonable agreement with the Bi2212 , as expected, since was chosen to match the Bi2212 low temperature LDOS(3).

In the BdG approach, the symmetric local density of states (LDOS) is proportional to the spectral function(26) and may be written as

(2) | |||||

The prime is the derivative with respect to the argument. and are respectively the eigenvectors and eigenvalues of the BdG matrix equation(18); (21), is the Fermi function and is the applied voltage. is proportional to the tunneling conductance , and we probe the effects of the inhomogeneous field by examining the ratio . contains the inhomogeneous charge distribution and vanishes around the Fermi energy due to the superconducting gap and the single particle levels in the grains. This LDOS ratio yields well-defined peaks and converges to the unity at large bias, in close agreement to similar LDOS ratios calculated from STM measurements(5).

We calculated the LDOS at many values of and to compare with the STM(2); (3); (4); (5); (6); (7) and ARPES(10); (11); (12); (13) data, but here we present detailed results only for an underdoped () and an overdoped () sample. Thus, in Fig.(3) we show the results at two representative places (shown as white dots in the middle of Fig.(1)) of the sample at various temperatures. In this case, the low temperature LDOS display very few well defined (coherent) peaks but the gaps are similar to the STM data of McElroy et al(3). The low temperature superconducting gap produces small anomalies that are marked by arrows, as reported by some STM data(6); (7); (8). This behavior was already noticed by our previous work(14). Then the gaps derived from the LDOS peaks are identified with the local pseudogap because they are larger than the local and because they vanish at a higher temperature than .

Let us investigate in detail the relationship between and . For the compound with K at a typical hole-poor place with , we obtain a low temperature meV and meV. In a representative hole-rich location with there is an interchange of intensity: meV and the superconducting gap meV. We observe that this inversion in the values of and is common to all compounds. For this sample, at all sites, and vanish with increasing temperature almost together near K. At the hole poor grain, and vanish near K, and at the hole rich grain, they close at and K. We see that, for underdoped samples, and they remain much above the resistivity transition K, i.e., both gaps are present in the pseudogap phase up to .

The difference between and increases continuously with increasing doping. Calculations with show that in some grains K while in some other locations . By increasing the doping , this difference increases as shown in Fig.(4) for . For this compound with K, most of the LDOS present well defined coherent peaks as observed by McElroy et al(3). On the top panel, we show the LDOS at a hole-rich site ; we see that meV and meV (marked by arrows in the plots) and they both vanish at . At a hole-poor site of , meV and meV. Following the temperature evolution we see that that is also the critical temperature , and vanishes at much larger value . Thus, the pseudogap phase of overdoped samples is composed mainly by .

The above results lead us to many conclusions concerning the HTSC measured properties: i) The charge inhomogeneities and the charge segregation potential occur due to the formation of (almost zero doping) AF regions. Consequently, its effect is more intense at low doping compounds, leading to large superconducting amplitudes and in underdoped samples. On the other hand, at far overdoped region, the large Coulomb repulsion in the hole-rich grains destroys the EPS transition and the superconductivity. ii) As in granular superconductors, the resistivity transition occurs by Josephson coupling among the intragrain superconducting regions. iii) In this scenario, the pseudogap is due to the weakly localized energy levels in the two dimensional puddles and it is not directly related to the intragrain superconductivity, although both phenomena are originated by the segregation potential . This different process was demonstrated by the distinct behavior of the two signal with applied magnetic fields and also with the temperature in tunneling experiments(27). iv) The observed difference in the LDOS form, called ”coherent” and ”zero temperature pseudogap”(3), is due to the LDOS at hole-rich and hole-poor locations respectively. This distinction is clearly seen in the overdoped calculations (Fig.(4)) where most LDOS presents coherent peaks in opposition to the underdoped case (as in Fig.(3)) where the rounded and ill-defined peaks are more abundant. v) The hole-poor regions have basically one electron per unit cell and the large Coulomb repulsion generates a large asymmetry between electron extraction and injection as measured by the STM experiments(2); (3); (4); (5); (6); (7); (8). The LDOS at hole-rich regions, like a homogeneous system, are expected to be more symmetric. vi) The observed anomaly or kink measured in the LDOS at very low bias(6); (7); (8) was demonstrated to be due to the low temperature superconducting gap . Fig.(3) shows clearly a small kink marked by the arrows. In general, it is more visible in underdoped samples because the is much larger than . vii) For underdoped samples, as in Fig.(3), and remain finite much above and vary very little around K. The presence of the superconducting gap and its quasiparticle dispersion above was measured in weakly underdoped Bi2212 by ARPES experiments(12) that also showed that the gap () almost does not change around . viii) The ARPES experiment of Lee et al(10) has measured increasing gapless Fermi arcs along the nodal () region with increasing doping above . As we have shown, the d-wave remains above in the underdoped regions but tends to vanish close to for overdoped compounds. Assuming that the weakly bound states () are due to the random phase boundary potential , and that they occur mainly along the (antinode) direction, together with the d-wave behavior of is agreement with the increase of the gapless Fermi arcs with above (10). ix) The calculations show that, in general, for any compound, the pseudogap is smaller but with well defined peaks at the hole-rich regions than in the hole-poor ones. Consequently, the measured larger gap values located at the low density grains with local (AF) insulator behavior have lower local conductivity () than the smaller gaps (and higher densities) as verified by Pasupathy et al(5). x) The STM measured relation (4) is reproduced closely by the optimal and the sample, considering and its associated .

In summary, we have proposed an EPS transition to describe the inhomogeneous charge distribution of HTSC generated by the lower free energy of the (undoped) AF domains. This approach yields potential wells with shallow bound states that reduce the kinetic energy and favor the superconducting pairing. The calculations obtain the inverted bell shape critical line , distinguish clearly the LDOS pseudogap from the intragrain superconducting gap energy and provide simple physical interpretations to many different STM and ARPES results.

We gratefully acknowledge partial financial aid from Brazilian agency CNPq.

### References

- J.M.Tranquada, et al Nature (London),375, 561 (1995).
- S. H. Pan et al, Nature (London), 413, 282-285 (2001).
- K. McElroy, et al cond-mat/0404005 an Phys. Rev. Lett. 94, 197005 (2005).
- Kenjiro K. Gomes et al, Nature 447, 569 (2007).
- Abhay N. Pasupathy et al, Science 320, 196 (2008).
- Takuya Kato et al, J. Phys. Soc. Jpn., 77, 054710 (2008).
- Aakash Pushp et al, Science 324 1689 (2009).
- T. Kato et al, J. Supercond. Nov. Magn. 23, 771 (2010).
- M. Le Tacon, et al, Nature Phys. 2, 537 (2006).
- W. S. Lee, et al, Nature 450, 81 (2007)
- J. H. Ma, et al, Phys. Rev. Lett. 101, 207002 (2008).
- A.Kanigel et al, Phys. Rev. Lett. 101, 137002 (2008).
- U. Chatterjee, et al Nature Phys. 6, 99 (2010).
- E.V.L. de Mello, et al J. Phys.: Condens. Matter 21, 235701 (2009).
- S. Wakimoto at al Phys. Rev. Lett. 98, 247003 (2007)
- J.W. Cahn and J.E. Hilliard, J. Chem. Phys, 28, 258 (1958).
- E.V.L de Mello et al, Physica A 347, 429 (2005).
- E.V.L. de Mello et al, Phys. Rev. B70, 224517 (2004).
- L.D. Landau, and E.M. Lifshitz, “Quantum Mechanics”, Pergamon Press, New York, 1977.
- E. V. L. de Mello, et al, Phys. Rev. B58, 9098 (1998).
- D. N. Dias et al, Phys. C468, 480 (2008).
- E.V.L. de Mello et al, Physica B404, 3119 (2009).
- L. Merchant et al, Phys. Rev. B63, 134508 (2001).
- V. Ambeogakar, and A. Baratoff, Phys. Rev. Lett. 10, 486 (1963).
- H. Takagi et al, Phys. Rev. Lett. 69, 2975 (1992).
- François Gygi, and Michael Schlüter, Phys. Rev. B43, 7609 (1991).
- V. M. Krasnov et al, Phys. Rev. Lett. 86, 2657 (2001).