# A universal explanation of tunneling conductance in exotic superconductors

A longstanding mystery in understanding cuprate superconductors is the inconsistency between the experimental data measured by scanning tunneling spectroscopy (STS) and angle-resolved photoemission spectroscopy (ARPES). In particular, the gap between prominent side peaks observed in STS is much bigger than the superconducting gap observed by ARPES measurements. Here, we reconcile the two experimental techniques by generalising a theory which was previously applied to zero-dimensional mesoscopic Kondo systems to strongly correlated two-dimensional (2D) exotic superconductors. We show that the side peaks observed in tunneling conductance measurements in all these materials have a universal origin: They are formed by coherence-mediated tunneling under bias and do not directly reflect the underlying density of states (DOS) of the sample. We obtain theoretical predictions of the tunneling conductance and the density of states of the sample simultaneously and show that for cuprate and pnictide superconductors, the extracted sample DOS is consistent with the superconducting gap measured by ARPES.

Tunneling conductance of non-strongly-correlated materials (where the electron–electron interactions do not have a fundamental effect on the physics) is simply understood because the as a function of (where is the voltage difference between the sample and the tip and is the measured current) is known to correspond to the DOS of the sample. However, for strongly correlated materials (SCMs) this is not true because the strong on-site interactions fundamentally change the tunneling mechanism by which electrons move from the current source to the sample itself. Nevertheless, this interpretation is frequently used when discussing SCMs, which can lead to confusing inconsistencies. For example, a sharp peak indicating the superconducting gap of cuprate superconductors is observed in the ARPES data near the nodal region Lee-Nature450 (); Tanaka-Science314 (), but in the STS data, two prominent side peaks occur at much higher energy Lawler-Nature466 (); Pushp-Science324 (); Fischer-RMP79 (). Also, STS data for strongly correlated pnictide superconductors display two prominent side peaks which are much bigger than the energy scale of the superconducting transition temperature Chi-PRL109 (). Therefore, explaining the structure of tunneling conductance of SCMs is one of the most challenging and urgent subjects in condensed matter physics. This situation will remain until a comprehensive theory taking into account the non-equilibrium nature of the STS experiment and the strong electron–electron interactions has been developed.

The main purpose of this study is to prove that the two prominent side peaks in the STS data are produced by the interplay between strong electron correlations in the sample and the non-equilibrium situation imposed by the experimental setup. An additional highlight of this study, just as important as the interpretation of the side peaks, is that we can effectively find the DOS of the specific 2D SCM under study within the fitting procedure in the theory. We find that the DOS predicted for cuprate and pnictide superconductors are consistent with the data given by ARPES Tanaka-Science314 (); Umezawa-PRL108 ().

The non-equilibrium calculation of for SCMs is achieved by generalizing a theory for the non-equilibrium Kondo effect in zero-dimensional mesoscopic systems Hong-JPCM23a (); Hong-JPCM23b () to the extended 2D situation. Some mathematical details are presented in the Methods section and in the Supplementary Materials, but in summary, the non-equilibrium tunneling dynamics are encoded within the non-equilibrium many-body Green’s function for the mediating site (MS, indicated by the subscript “”), , which is derived from the Liouville operator. Once this is known, the tunneling conductance can be computed from the formula in equation (2). The advantage of using the Liouvillian approach instead of the Hamiltonian approach is the availability of a complete set of basis vectors, which is not possible in the latter case. The difficulty with the Liouvillian approach arises in the calculation of Liouville matrix elements for a specific model system. This problem can be resolved by making phenomenological assumptions or treating the unknowns as fitting parameters. However, for the Hamiltonian approach, one encounters a fundamental difficulty in obtaining a complete set of basis vectors at the beginning of the calculation and therefore progress is blocked. The DOS in the leads are captured by the lead DOS functions and which enter via and in equation (2). Hence, these functions are inputs to the theory, and can be chosen to represent various different scenarios.

We now describe the coherent tunneling processes which produce the dominant features of the experimental curves. Our central claim is that many-body singlet states form between electrons on the MS and in both leads (as shown by the ellipses in figures 1, 2, and 3). When a bias voltage is applied, the coherent singlet cotunnels unidirectionally through the MS. In this context, coherence implies that the singlet may perform spin exchange or change its partner in the leads without any energy cost. A mesoscopic Kondo system is sketched in figure 1a and consists of two leads with constant DOS and a central zero-dimensional MS. This system can be described by a two-reservoir Anderson impurity model Anderson-PhysRev124 (). To understand the coherent tunneling mechanism of the entangled state, we must examine the elements of the Liouville operator shown in equation (1). These are illustrated schematically in figure 2 and their full operator expressions are given in equation (S6) in the Supplementary Information. To put this discussion in context, figure 1b shows an archetypal plot of the tunneling conductance as a function of bias voltage for a quantum point contact (QPC) with symmetric coupling to the leads on both sides Sarkozy-PRB79 (). Notice there are three coherent peaks: one “zero bias peak”, and two “side peaks” which cannot be Coulomb peaks because of the small energy scale. We shall demonstrate which of the elements contribute to the formation of these peaks. We first discuss the formation of the side peaks. Figure 2a and its mirror image 2b show coherent current formed solely from singlet hopping. In this process, a left singlet performs a singlet hop to the right (second sketch) before partner exchange on the left (third sketch) enables a second singlet hop to the right which returns the MS to its original state (fourth sketch). This coherent process allows two electrons (highlighted in yellow in the sketch) to tunnel from the left lead to the right lead without paying the on-site Coulomb repulsion energy cost . Matrix elements and in equation (1) are given by symmetric and antisymmetric combinations and . It is clear that at equilibrium (i.e., when no bias is applied across the system), and will have the same magnitude because the singlet hops will occur with equal probability in each direction, so that and , indicating that contributes to the zero bias peak. However, out of equilibrium, either or will be suppressed because it must work against the external bias in the system so that and as a result, . This contributes an additional channel for current at finite bias, and this manifests as the side peaks in away from equilibrium. Therefore, it is the inclusion of left-right antisymmetric coherent superpositions which allow us to explain the existence of the side peaks in the tunneling conductance. Additionally, , where (figure 2c) represents the Kondo coupling strength on the left side, and (figure 2d) describes the current flow forming the zero bias peak, in which spin exchange is involved. However, the zero bias peak is suppressed in the 2D SCMs because of the interaction-induced gap at the Fermi level and so we do not focus on it in this study. We briefly mention that the elements are not strongly temperature dependent unless the temperature gets so high that thermal fluctuations destroy the coherence. In fact, the main thermal effects will be driven by the changes in the Fermi-Dirac distribution of the leads, and changes in the sample DOS itself.

In real STS experiments for cuprate superconductors, every copper atom has strong on-site repulsion (shown in figure 3a). We set up our model for 2D cuprate superconductors by stating that the tunneling current coming from the STS tip only enters through one copper atom. Then, the tunneling mechanism depends on the strong at that site only so that the rest of the sample can then be modeled as a system of non-interacting Bogoliubov quasiparticles Lee-Nature450 (); Yang-Nature456 (); Matsui-PRL90 () excited by applied bias, as sketched in figure 3b. Here, an “entangled state” comprising a linear combination of two singlets (solid and dashed purple loops) accompanied by other coherent spins in the tip (blue and red dots) and the sample is formed. The Hamiltonian for the left lead is chosen to be appropriate for the metallic tip so that is constant, the Hamiltonian for the right lead is taken to be a system of Bogoliubov quasiparticles Lee-Nature450 (); Yang-Nature456 (); Matsui-PRL90 () whose frequency-dependent DOS characterizes the superconductor. Therefore, the DOS of the superconductor becomes an input to the theory [ in equations (1) and (2)] and hence we can use it as a fitting function to distinguish between various theoretical propoals for the superconductor DOS. This allows us to separate the physics of the tunneling mechanism by coherent co-tunneling of singlet states through the MS from the physics of the superconductor in the bulk of the sample.

Now, we apply the theory described above to correlated superconductors to show that it gives accurate predictions for published STS data and that the features in the sample DOS lead to the fine structure in the low-bias region. We require a phenomenological frequency-dependent sample DOS to use as . We find that a DOS characteristic of a -wave superconducting order parameter for the under-doped (UD) cuprate and a DOS composed of an -wave gap and a sharp DOS barrier for the pnictide give a well-fitting tunneling conductance. In both cases, we use flat DOS in the high frequency region to emphasize that the side peaks are not created by features in the sample DOS. These DOS functions are shown as the green lines (right-hand axis) in figures 4a and 4b. Choosing the parameters given in the second and third row of Tab. 1 gives the fit (blue line) to the experimental data (red line). The quantitative agreement up to the side peaks is almost exact, and we reproduce all the low-bias features. The peaks of the DOS are located at (figure 4a) and (figure 4b), which are consistent with the superconducting gap reported in recent ARPES data Tanaka-Science314 (); Umezawa-PRL108 () and correspond to the low-energy shoulders in the . The result is natural since the Bogoliubov quasiparticles introduced in figure 3b describe the coherent superconducting state Lee-Nature450 (). We also show the and the DOS for optimally doped (OP) and highly overdoped (HOD) cases in the insets of figure 4a. It is noteworthy that the two side peaks in the UD and OP cases are formed in the region of flat DOS. This demonstrates that the side peaks are not connected to the superconducting gap, and they must be reinterpreted as a coherent peak of the non-equilibrium transport mechanism. Their sustained presence when signifies the existence of pre-formed Cooper pairs Yang-Nature456 () up to of cuprate superconductors. Since a big gap observed in the anti-nodal region in the ARPES measurements is not a coherent gap Lee-Nature450 (); Tanaka-Science314 (), the anti-nodal gap is not included in this study. A big difference occurs in the HOD case where strong correlation no longer persists. We obtain the HOD line shape using vanishing elements for a very weak effective Coulomb interaction. Note that two peaks occur within the -wave region.

In conclusion, we have generalised a non-equilibrium quantum transport theory which accurately models tunneling conductance measurements for extended 2D strongly correlated systems. This theory is therefore crucial in understanding experimental data in a wide variety of different contexts. A particularly interesting point is that our theoretical tunneling conductance and sample DOS for cuprate and pnictide superconductors are consistent with the results obtained by two different experimental techniques. This fact eliminates the temptation to think that the curve is simply a reflection of the sample DOS (as it would be in the weakly interacting case), and therefore to interpret the two side peaks as evidence of a correlated band gap in the sample DOS. We demonstrate that the correct interpretation is that the side peaks are evidence of the coherent non-equilibrium tunneling mechanism.

Figure | |||||||||
---|---|---|---|---|---|---|---|---|---|

1b | 0.53 | 0.53 | 0.68 | 1.0 | 1.17 | 1.17 | 0 | 0.8 | 0.99 |

4a | 0.01 | 0.5 | 0.5 | 1.75 | 0.35 | 2.45 | 0.1 | 0.24 | 118 |

4b | 0.01 | 0.5 | 0.5 | 3.5 | 0.7 | 3.85 | 0 | 0.24 | 13.4 |

,
,
,

,
and .

## Methods.

After the appropriate modelling described in the text, the STS experiment may be described by the following Hamiltonian:

where and correspond to the left and right leads, which are taken to be a metallic tip and a system of non-interacting Bogoliubov quasiparticles (described in detail in Section III of the Supplementary Information) excited by the applied bias at , respectively. The Hamiltonian for the MS is , and describes the tunneling between the MS and the leads. As described in the Supplementary Information, sections IB and IC, the Liouville operator describing the dynamics of the MS is

(1) |

where denotes the self-energy matrix with elements . The coefficients come from the process of matrix reduction which generates equation (1). Each matrix element has its own unique role in determining . The middle row and column of the matrix are the incoherent couplings of the MS to the leads which represent double occupancy at the MS via the fluctuations of incoherent spins from the tip or the sample. Therefore, these elements represent the effective Coulomb interaction. The doping effect is given by the imaginary part of these elements shown in equation (S7). The subscript indicates the current operator on side , whereas indicates the sum of the left and right movements.

Within this formalism, the tunneling conductance at zero temperature can be obtained as discussed in the Supplementary Information, sections IB and IE:

(2) |

where is the effective coupling of the left lead, the MS, and the right lead and is given by

(3) |

for 2D SCMs. The quantity , where the matrix is given by the elements , is the non-equilibrium local DOS at the MS (which is distinct from the sample DOS) and depends on the strong interactions. It is computed from the Liouville matrix in equation (1) as described in the Supplementary Information.

## Acknowledgements.

This work was supported by Project Code (IBS-R024-D1), the Basic Science Research Program via the NRF Korea (2012R1A1A2005220), Nordita, ERC project DM-321031, and partially supported by a KIAS grant funded by MSIP. JH thanks Piers Coleman and Peter Fulde for helpful discussions.

## Author contributions

J.H. designed the study, performed the calculations, and wrote the manuscript. D.S.L.A designed the study and wrote the manuscript.

## Competing financial interests

The authors declare no competing financial interests.

### References

- Lee, W. S. et al. Abrupt onset of a second energy gap at the superconducting transition of underdoped Bi2212. Nature 450, 81–84; DOI:10.1038/nature06219 (2007).
- Tanaka, K. et al. Distinct Fermi-momentum-dependent energy gaps in deeply underdoped Bi2212. Science 314, 1910–1913; DOI:10.1126/science.1133411 (2006).
- Lawler, M. J. et al. Intra-unit-cell electronic nematicity of the high- copper-oxide pseudogap states. Nature 466, 347–351; DOI:10.1038/nature09169 (2010).
- Pushp, A. et al. Extending universal nodal excitations optimizes superconductivity in BiSrCaCuO. Science 324, 1689–1693; DOI:10.1126/science.1174338 (2009).
- Fischer, O., Kugler, M., Maggio-Aprile, I., Berthod, C. & Renner, C. Scanning tunneling spectroscopy of high-temperature superconductors. Rev. Mod. Phys. 79, 353–419; DOI:10.1103/RevModPhys.79.353 (2007).
- Chi, S. et al. Scanning tunneling spectroscopy of superconducting LiFeAs single crystals: Evidence for two nodeless energy gaps and coupling to a bosonic mode. Phys. Rev. Lett. 109, 087002; DOI:10.1103/PhysRevLett.109.087002 (2012).
- Umezawa, K. et al. Unconventional anisotropic -wave superconducting gaps of the LiFeAs iron-pnictide superconductor. Phys. Rev. Lett. 108, 037002; DOI:10.1103/PhysRevLett.108.037002 (2012).
- Sarkozy, S. et al. Zero-bias anomaly in quantum wires. Phys. Rev. B 79, 161307; DOI:10.1103/PhysRevB.79.161307 (2009).
- Hong, J. A complete set of basis vectors of the Anderson model and its Kondo dynamics. J. Phys: Conden. Matt. 23, 225601; DOI:10.1088/0953-8984/23/22/225601 (2011).
- Hong, J. Kondo dynamics of quasiparticle tunneling in a two-reservoir Anderson model. J. Phys: Conden. Matt. 23, 275602; DOI:10.1088/0953-8984/23/27/275602 (2011).
- Anderson, P. W. Localized magnetic states in metals. Phys. Rev. 124, 41–53; DOI:10.1103/PhysRev.124.41 (1961).
- Yang, H.-B. et al. Emergence of preformed Cooper pairs from the doped Mott insulating state in BiSrCaCuO. Nature 456, 77-80; DOI:10.1038/nature07400 (2008).
- Matsui, H. et al. BCS-like Bogoliubov quasiparticles in high- superconductors observed by angle-resolved photoemission spectroscopy. Phys. Rev. Lett. 90, 217002; DOI:10.1103/PhysRevLett.90.217002 (2003).