This article is part of the series EMN Meeting.

Open Access Nano Express

Doping dependence of low-energy quasiparticle excitations in superconducting Bi2212

Akihiro Ino1*, Hiroaki Anzai2, Masashi Arita3, Hirofumi Namatame3, Masaki Taniguchi13, Motoyuki Ishikado456, Kazuhiro Fujita47, Shigeyuki Ishida45 and Shinichi Uchida4

Author Affiliations

1 Graduate School of Science, Hiroshima University, Higashi-Hiroshima 739-8526, Japan

2 Graduate School of Engineering, Osaka Prefecture University, Sakai 599-8531, Japan

3 Hiroshima Synchrotron Radiation Center, Hiroshima University, Higashi-Hiroshima 739-0046, Japan

4 Department of Physics, University of Tokyo, Tokyo 113-0033, Japan

5 National Institute of Advanced Industrial Science and Technology, Tsukuba 305-8568, Japan

6 Research Center for Neutron Science and Technology, CROSS, Tokai, Ibaraki 319-1106, Japan

7 Laboratory for Atomic and Solid State Physics, Department of Physics, Cornell University, Ithaca, New York 14853, USA

For all author emails, please log on.

Nanoscale Research Letters 2013, 8:515  doi:10.1186/1556-276X-8-515


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/8/1/515


Received:17 July 2013
Accepted:20 November 2013
Published:5 December 2013

© 2013 Ino et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The doping-dependent evolution of the d-wave superconducting state is studied from the perspective of the angle-resolved photoemission spectra of a high-Tc cuprate, Bi2Sr2CaCu2 O8+δ (Bi2212). The anisotropic evolution of the energy gap for Bogoliubov quasiparticles is parametrized by critical temperature and superfluid density. The renormalization of nodal quasiparticles is evaluated in terms of mass enhancement spectra. These quantities shed light on the strong coupling nature of electron pairing and the impact of forward elastic or inelastic scatterings. We suggest that the quasiparticle excitations in the superconducting cuprates are profoundly affected by doping-dependent screening.

Keywords:
High-Tc cuprate; Bi2212; ARPES; Superconducting gap; Effective mass; Coupling strength; 74.25.Jb; 74.72.-h; 79.60.i

Background

Electronic excitations dressed by the interaction with the medium are called quasiparticles. They serve as a direct probe of the anisotropic order parameter of a superconducting phase and also as a clue to the electron-pairing glue responsible for the superconductivity. In fact, the major unresolved issues on the mechanism of high-Tc superconductivity depend on the low-energy quasiparticle excitations. The superconducting order parameter, which is typified by the particle-hole mixing and gives rise to Bogoliubov quasiparticles (BQPs), manifests itself as an energy gap in quasiparticle excitation spectra. In cuprate superconductors, however, the energy gap increases against the decrease in critical temperature Tc with underdoping and is open even at some temperatures above Tc[1-3]. In the direction where the d-wave order parameter disappears, renormalization features have been extracted quantitatively from the gapless continuous dispersion of nodal quasiparticles (NQPs), suggesting strong coupling with some collective modes [4]. Nevertheless, the origins of these features remain controversial [4,5].

In this paper, we address the doping dependence of BQP and NQP of a high-Tc cuprate superconductor, Bi2Sr2CaCu2O8+δ (Bi2212), on the basis of our recent angle-resolved photoemission (ARPES) data [6-8]. The use of low-energy synchrotron radiation brought about improvement in energy and momentum resolution and allowed us to optimize the excitation photon energy. After a brief description of BQP and NQP spectral functions, we survey the superconducting gap anisotropy on BQPs and the renormalization features in NQPs. In light of them, we discuss possible effects of doping-dependent electronic screening on the BQP, NQP, and high-Tc superconductivity.

Methods

High-quality single crystals of Bi2212 were prepared by a traveling-solvent floating-zone method, and hole concentration was regulated by a post-annealing procedure. In this paper, the samples are labeled by the Tc value in kelvin, together with the doping-level prefix, i.e. underdoped (UD), optimally doped (OP), or overdoped (OD). ARPES experiments were performed at HiSOR BL9A in Hiroshima Synchrotron Radiation Center. The ARPES data presented here were taken with excitation-photon energies of hν = 8.5 and 8.1 eV for the BQP and NQP studies, respectively, and at a low temperature of T = 9 - 10 K in the superconducting state. Further details of the experiments have been described elsewhere [7-9].

The relation between a bare electron and a renormalized quasiparticle is described in terms of self-energy Σk(t), which can be regarded as a factor of feedback on the wave function from past to present through the surrounding medium. Incorporating a feedback term into the Schrödinger equation, we obtain

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M1">View MathML</a>

(1)

where ψk(t) and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M2">View MathML</a> denote a wave function and a bare-electron energy, respectively. It is obvious from Equation 1 that the self-energy is a linear response function. Therefore, its frequency representation, Σk(ω), obeys the Kramers-Kronig relation. As the solution of Equation 1, we obtain the form of dressed Green’s function,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M3">View MathML</a>

(2)

The spectral function given by Ak(ω) = - Im Gk(ω)/π is directly observed by ARPES experiments. The extensive treatments of the ARPES data in terms of Green’s function are given elsewhere [10].

Results

Superconducting gap anisotropy

In the superconducting state, the condensate of electron pairs allows the particle-like and hole-like excitations to turn into each other. Hence, the wave function of a hole-like excitation also comes back by way of the particle-like excitation. Such a feedback has a sign-reversed eigenenergy, <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M4">View MathML</a>, and is expressed by <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M5">View MathML</a>, where Θ(t), Δk, γk and Γk denote the step function, the particle-hole off-diagonal element, and the scattering rates of the intermediate and bare-particle states, respectively. The Fourier transform of Σk(t) gives the frequency representation of the self-energy of the BQPs,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M6">View MathML</a>

(3)

Figure 1 shows the ARPES spectra of BQPs for underdoped and overdoped Bi2212 samples with Tc = 66 and 80 K (UD66 and OD80, respectively) [8]. As shown in Figure 1b,c, an energy distribution curve was extracted from the minimum gap locus for each off-node angle θ and symmetrized with respect to the Fermi energy ω = 0. These spectra were well fitted with a phenomenological function,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M7">View MathML</a>

(4)

except for a featureless background. Equation 4 is deduced from Equation 3 and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M8">View MathML</a>, neglecting γk after Norman et al. [11]. Figure 1b,c exemplifies that the superconducting gap energy Δ at each θ is definitely determined by sharp spectral peaks. In Figure 1d,e, the obtained gap energies (small yellow circles) are plotted over the image of spectral intensity as a function of sin 2θ, so that the deviation from a d-wave gap is readily seen with reference to a straight line. While the superconducting gap of the overdoped sample almost follows the d-wave line, that of the underdoped sample is deeply curved against sin 2θ. Furthermore, Figure 1d indicates that the deviation from the d-wave gap penetrates into the close vicinity of the node and that it is difficult to define the pure d-wave region near the node. Therefore, the next-order harmonic term, sin 6θ, has been introduced, so that the smooth experimental gap profile is properly parametrized [12-14]. The next-order high-harmonic function is also expressed as Δ(θ) = ΔN sin 2θ + (ΔN)(3 sin 2θ- sin 6θ)/4, where the antinodal and nodal gap energies are defined as Δ = Δ(θ)|θ=45° and <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M9">View MathML</a>, respectively, so that ΔN = 1 is satisfied for a pure d-wave gap.

thumbnailFigure 1. Superconducting gap manifested in BQP spectra. The data are for underdoped and overdoped Bi2212 samples with Tc = 66 and 80 K (UD66 and OD80, respectively) [8]. (a) Momentum-space diagram for an off-node angle, θ, and a bonding-band (BB) Fermi surface along which the ARPES spectra were taken. (b, c) Symmetrized energy distribution curves (colored circles) and their fits (black curves). (d, e) ARPES spectral images as a function of energy ω and sin 2θ. Superimposed are the gap energies (yellow circles) and high-harmonic fit (yellow curve) as functions of sin 2θ.

The doping dependences of the superconducting gap parameters are summarized in Figure 2. One can see from Figure 2a that as hole concentration decreases with underdoping, the nodal gap energy 2ΔN closely follows the downward curve of 8.5kBTc in contrast to the monotonic increase in the antinodal gap energy 2Δ. It seems reasonable that Tc primarily depends on ΔN rather than Δ for the underdoped Bi2212, because it follows from 2Δ ≫ 4kBTc that the thermal quasiparticle excitations concentrate in the vicinity of the node and hardly occur around the antinode. The relevance of the nodal excitations has also been suggested by various experiments [15-19]. Then, the problem with Tc is that the nodal gap ΔN is suppressed relative to the antinodal gap Δ. This behavior can be associated with low superfluid density ρs[20]. Figure 2b,c shows that the doping dependence of the nodal-to-antinodal gap ratio ΔN is quite similar to that of the square-root superfluid density <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M10">View MathML</a>[8,21,22]. The normalized gap plot in Figure 2d indicates that what occurs with underdoping is analogous to the nodal gap suppression observed with increasing temperature [17] in terms of the decrease in ρs. It is notable that the square-root dependence on ρs is a typical behavior of the order parameter as expected from the Ginzburg-Landau theory [23]. These findings can be written down in a simple relational formula,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M11">View MathML</a>

(5)

where <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M12">View MathML</a>, for a wide hole-concentration range of Bi2212.

thumbnailFigure 2. Doping dependences of superconducting gap parameters. (a) Nodal gap energy 2ΔN (blue circles) and antinodal gap energy 2Δ (red squares) [8]. The solid curve denotes an energy of 8.5kBTc. (b) Square of nodal-to-antinodal gap ratio (ΔN)2 determined from ARPES [8]. (c) Superfluid density ρs determined from magnetic penetration depth (triangles) [21] and from heat capacity (crosses) [22]. (d) Superconducting gap profiles normalized to the antinodal gap for underdoped and optimally doped samples with Tc = 42, 66, and 91 K (UD42, UD66, and OP91, respectively). (e) Correlation between 2ΔN/kBTc and 2Δ/kBTc ratios. The insets illustrate the occurrence of incoherent electron pairs in strong coupling superconductivity.

As presented in Figure 2e, the correlation between the nodal and antinodal gaps provides a perspective of crossover for our empirical formula (Equation 5). It is deduced from the conventional Bardeen-Cooper-Schrieffer (BCS) theory that 2Δ/kBTc = 4.3 in the weak coupling limit for the d-wave superconducting gap [23]. However, the critical temperature Tc is often lower than that expected from the weak coupling constant and a given Δ as an effect of strong coupling. Thus, the gap-to- Tc ratio is widely regarded as an indicator for the coupling strength of electron pairing and adopted for the coordinate axes in Figure 2e. As hole concentration decreases from overdoped to underdoped Bi2212, the experimental data point moves apart from the weak coupling point toward the strong coupling side, and a crossover occurs at 8.5, which is about twice the weak coupling constant. It appears that the evolution of ΔN is confined by two lines as ΔN ≤ 0.87Δ and 2ΔN ≤ 8.5kBTc. As illustrated in the insets of Figure 2e, the strong coupling allows the electrons to remain paired with incoherent excitations. As a result, the superconducting order parameter is reduced with respect to the pairing energy. Indeed, it has been shown that the reduction factor due to the incoherent pair excitations has a simple theoretical expression <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M13">View MathML</a> and that the nodal and antinodal spectra are peaked at the order parameter and at the pairing energy, respectively, taking into account a realistic lifetime effect [24,25]. Therefore, the latter part of Equation 5 is consistent with the strong coupling scenario, and furthermore, the two distinct lines in Figure 2e are naturally interpreted as the energies of the condensation and formation of the electron pairs.

Renormalization features in dispersion

In the nodal direction where the order parameter disappears, one can investigate the fine renormalization features in dispersion. They reflect the intermediate-state energy in coupling between an electron and other excitations, and thus provide important clues to the pairing interaction. As for the electron-boson coupling, the intermediate state consists of a dressed electronic excitation and an additional bosonic excitation (Figure 3a). Averaging the momentum dependence for simplicity, the energy distribution of the intermediate state is expressed by A(ω - Ω) Θ(ω - Ω)+A(ω + Ω) Θ(-ω - Ω) for a given boson energy Ω and for zero temperature, owing to the Pauli exclusion principle. Therefore, taking into account the effective energy distribution of the coupled boson, α2F(Ω), the self-energy is written down as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M14">View MathML</a>

(6)

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M15">View MathML</a>

(7)

where 0+ denotes a positive infinitesimal.

thumbnailFigure 3. Simulation for a single coupling mode at Ω = 40 meV. Dotted and solid curves denote those with and without a d-wave gap of Δ = 30 meV, respectively. (a) Diagram of electron-boson interaction. (b) Eliashberg coupling function α2F(-ω), dispersion k(ω) = [ω + ReΣ(ω)]/v0, and momentum width Δk(ω) = -ImΣ(ω)/v0. (c) Real and imaginary parts of 1 + λ(ω).

In ARPES spectra, the real and imaginary parts of self-energy manifest themselves as the shift and width of spectral peak, respectively. Specifically, provided that the momentum dependence of Σk(ω) along the cut is negligible, and introducing bare electron velocity v0 by <a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M16">View MathML</a>, it follows from Equation 2 that the momentum distribution curve for a given quasiparticle energy ω is peaked at k(ω) = [ω-ReΣ(ω)]/v0 and has a natural half width of Δk(ω) = - ImΣ(ω)/v0.

We argue that the mass enhancement function defined as the energy derivative of the self-energy, λ(ω) ≡ -(d/dω)Σ(ω), is useful for the analysis of NQP [7,26]. The real and imaginary parts of λ(ω) are directly obtained from the ARPES data as the inverse of group velocity, vg(ω), and as the differential scattering rate, respectively.

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M17">View MathML</a>

(8)

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M18">View MathML</a>

(9)

We note that -Imλ(ω) represents the energy distribution of the impact of coupling with other excitations and can be taken as a kind of coupling spectrum. However, it should be emphasized that -Imλ(ω) is expressed as a function of quasiparticle energy ω, whereas the widely used Eliashberg coupling function α2F(Ω) is expressed as a function of boson energy Ω. For example, a simulation of λ(ω) using Equations 7 to 9 is presented in Figure 3b,c, where a single coupling mode is given at Ω = 40 meV. One can see that the peak of α2F(-ω) is reproduced by -Imλ(ω), provided that A(ω) is gapless and approximated by a constant. As an energy gap of Δ opens in A(ω), the peak in -Imλ(ω) is shifted from Ω into Ω + Δ. Nevertheless, irrespective of A(ω), the causality of Σ(ω) is inherited by λ(ω), so that Reλ(ω) and Imλ(ω) are mutually convertible through the Kramers-Kronig transform (KKT). The directness and causality of λ(ω) enable us to decompose the quasiparticle effective mass without tackling the integral inversion problem in Equation 7.

Figure 4 shows the ARPES spectra along the nodal cut perpendicular to the Fermi surface for the superconducting Bi2212 [7]. Although the splitting due to the CuO2 bilayer is minimum at the nodes, it has clearly been observed by using some specific low-energy photons [6-8]. A prominent kink in the NQP dispersion is observed at 65 meV for all the doping level, as has been reported since early years [4]. In addition to this, another small kink at 15 meV is discernible in the raw spectral image of the underdoped sample (UD66) [7,27].

thumbnailFigure 4. Dispersion kinks manifested in NQP spectra. The ARPES spectra were taken in the superconducting state for Bi2212 [7]. (a) Underdoped sample with Tc = 66 K (UD66). (b) Optimally doped sample with Tc = 91 K (OP91). (c) Overdoped sample with Tc = 80 K (OD80).

The fine renormalization features in the NQP dispersion were determined by fitting the momentum distribution curves with double Lorentzian. Figure 5a,d shows the real and imaginary parts of λ(ω)/v0 experimentally obtained as the energy derivatives of the peak position and width, respectively. The KKT of Reλ(ω)/v0 in Figure 5a is shown in Figure 5b as Imλ(ω)/v0, which is comparable with the data in Figure 5d. A step-like mass enhancement in Figure 5a and a peak-like coupling weight in Figure 5b,d are consistently observed at 65 meV. This is a typical behavior of the mode coupling, as shown by the simulation in Figure 3. It is also found that an additional feature around 15 meV is dramatically enhanced with underdoping. In order to deduce the partial coupling constant, we express the mass enhancement factor λ as the form of KKT,

<a onClick="popup('http://www.nanoscalereslett.com/content/8/1/515/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/8/1/515/mathml/M19">View MathML</a>

(10)

thumbnailFigure 5. Doping dependences of NQP properties. The real and imaginary parts of mass enhancement spectra were directly deduced from the APRES data shown in Figure 4[7]. (a) Inverse group velocity, 1/vg(ω) = [1 + Re λ(ω)]/v0, determined from (d/dω) k(ω). (b) Differential scattering rate -Im λ(ω)/v0, deduced from the Kramers-Kronig transform (KKT) of (a). (c) Partial coupling constants, λLE (red circles) and λIE (blue triangles), deduced from (b). Also shown are the inverse group velocities at ω = 0 (black circles) and at ω = -40 meV (black triangles). (d) Differential scattering rate -Im λ(ω)/v0, directly determined from -(d/dω) Δk(ω). (e) Partial coupling constants, λLE (red circles), λIE (blue triangles), and λHE (green diamonds), deduced from (d).

Dividing the energy range of the integral in Equation 10, one can quantify the contribution from a particular energy part. We refer to the KKT integrals of Im λ(ω)/v0 for the low-energy (LE; 4 < |ω| < 40 meV), intermediate-energy (IE; 40 < |ω| < 130 meV), and high-energy (HE; 130 < |ω| < 250 meV) parts as λLE/v0 (red circles), λIE/v0 (blue triangles), and λHE/v0 (green diamonds), respectively. Those obtained from the data in Figure 5b,d are plotted in Figure 5c,e, respectively. Also shown in Figure 5c are the inverse group velocities at ω = 0 meV (black circles) and at ω = -40 meV (black triangles). Figure 5c and Figure 5e consistently indicate that as hole concentration decreases, the contribution of the low-energy part rapidly increases and becomes dominant over the other parts.

Possible origins of the low-energy kink are considered from the energy of 15 meV and the evolution with underdoping. The quasiparticles that can be involved in the intermediate states are limited within the energy range of |ω| ≤ 15 meV, and the irrelevance of the antinodal states is deduced from the simulation in Figure 3c. Therefore, the low-energy kink is due to the near-nodal scatterings with small momentum transfer. The candidates for bosonic forward scatterers are the low-frequency phonons, such as the acoustic phonons and the c-axis optical phonons involving heavy cations [7,28-31]. On the other hand, it has also been argued that the elastic forward scattering by off-plane impurities may give rise to the low-energy kink for the d-wave superconductors [7,32]. In usual metal, both the potentials of the low-frequency phonons and the static impurities are strongly screened by the rapid response of electronic excitations. Therefore, the enhancement of the low-energy kink suggests the breakdown of electronic screening at low hole concentrations [7,28].

The dispersion kink at 65 meV has been ascribed to an intermediate state consisting of an antinodal quasiparticle and the B1g buckling phonon of Ω ∼ 35 meV [33]. However, the mass enhancement spectra in Figure 5a,b,d are suggestive of the presence of multiple components in the intermediate-energy range.

Discussion

We found that both the superconducting gap anisotropy and the renormalized dispersion show the striking evolution with underdoping. These behaviors are considered to be dependent on the extent of the screening. In association with the forward elastic or inelastic scatterings, the screening breakdown would enhance the low-energy kink. From the aspect of the impact of off-plane impurities, the inadequacy of static screening would inevitably lead to the nanoscale inhomogeneities, as observed by scanning tunneling microscopy experiments [34]. The forward scatterings by the remaining potential would generate additional incoherent pair excitations, as expected from the nodal gap suppression at low temperatures [8,25]. From the aspect of the electron-phonon coupling, the inhomogeneous depletion of the electrons for screening may considerably increase the coupling strength, providing an account for the unexpectedly strong dispersion kink [35] and a candidate for the strong pairing interaction [8]. The former and latter aspects have negative and positive effects, respectively, on the superconductivity. Thus, we speculate that the doping dependence of Tc is eventually determined by the balance between these effects.

Conclusions

Summarizing, the evolution of a d-wave high-Tc superconducting state with hole concentration has been depicted on the basis of the high-resolution ARPES spectra of the quasiparticles and discussed in relation to the screening by electronic excitations. The divergence between the nodal and antinodal gaps can be interpreted as an effect of the incoherent pair excitations inherent in the strong coupling superconductivity. The low-energy kink, which rapidly increases with underdoping, should be caused by the forward elastic or inelastic scatterings, although it remains as an open question which scattering is more dominant. The quantitative simulation of the doping-dependent effect will be helpful for resolving this problem.

Abbreviations

AB: Antibonding band; ARPES: Angle-resolved photoemission spectroscopy; BB: Bonding band; BCS: Bardeen-Cooper-Schrieffer; Bi2212: Bi2Sr2CaCu2 O8+δ; BQP: Bogoliubov quasiparticle; HE: High energy; IE: Intermediate energy; KKT: Kramers-Kronig transform; LE: Low energy; NQP: Nodal quasiparticle; OD: Overdoped; OP: Optimally doped; UD: Underdoped.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AI wrote the manuscript. HA and AI designed the experiment and analyzed the data with support from MT. HA acquired the ARPES data with support from AI, MA, and HN. High-quality single-crystalline samples were grown by MI, KF, SI, and SU. All authors discussed the results and commented on the manuscript. All authors read and approved the final manuscript.

Acknowledgements

We thank Z.-X. Shen and A. Fujimori for useful discussions and K. Ichiki, Y. Nakashima, and T. Fujita for their help with the experimental study. The ARPES experiments were performed under the approval of HRSC (Proposal No. 07-A-2, 09-A-11 and 10-A-24).

References

  1. Miyakawa N, Guptasarma P, Zasadzinski JF, Hinks DG, Gray KE: Strong dependence of the superconducting gap on oxygen doping from tunneling measurements on Bi2Sr2CaCu2O8-δ.

    Phys Rev Lett 1998, 80:157-160. Publisher Full Text OpenURL

  2. Campuzano JC, Ding H, Norman MR, Fretwell HM, Randeria M, Kaminski A, Mesot J, Takeuchi T, Sato T, Yokoya T, Takahashi T, Mochiku T, Kadowaki K, Guptasarma P, Hinks DG, Konstantinovic Z, Li ZZ, Raffy H: Electronic spectra and their relation to the (π,π) collective mode in high-Tc superconductors.

    Phys Rev Lett 1999, 83:3709-3712. Publisher Full Text OpenURL

  3. Loeser AG, Shen ZX, Dessau DS, Marshall DS, Park CH, Fournier P, Kapitulnik A: Excitation gap in the normal state of underdoped Bi2Sr2CaCu2O8+δ.

    Science 1996, 273:325-329. PubMed Abstract | Publisher Full Text OpenURL

  4. Lanzara A, Bogdanov PV, Zhou XJ, Kellar SA, Feng DL, Lu ED, Yoshida T, Eisaki H, Fujimori A, Kishio K, Shimoyama JI, Noda T, Uchida S, Hussain Z, Shen ZX: Evidence for ubiquitous strong electron-phonon coupling in high-temperature superconductors.

    Nature 2001, 412:510. PubMed Abstract | Publisher Full Text OpenURL

  5. Johnson PD, Valla T, Fedorov AV, Yusof Z, Wells BO, Li Q, Moodenbaugh AR, Gu GD, Koshizuka N, Kendziora C, Jian S, Hinks DG: Doping and temperature dependence of the mass enhancement observed in the cuprate Bi2Sr2CaCu2O8+δ.

    Phys Rev Lett 2001, 87(17):177007. PubMed Abstract | Publisher Full Text OpenURL

  6. Yamasaki T, Yamazaki K, Ino A, Arita M, Namatame H, Taniguchi M, Fujimori A, Shen ZX, Ishikado M, Uchida S: Unmasking the nodal quasiparticle dynamics in cuprate superconductors using low-energy photoemission.

    Phys Rev B 2007, 75:140513. OpenURL

  7. Anzai H, Ino A, Kamo T, Fujita T, Arita M, Namatame H, Taniguchi M, Fujimori A, Shen ZX, Ishikado M, Uchida S: Energy-dependent enhancement of the electron-coupling spectrum of the underdoped Bi2Sr2CaCu2O8+δ superconductor.

    Phys Rev Lett 2010, 105:227002. PubMed Abstract | Publisher Full Text OpenURL

  8. Anzai H, Ino A, Arita M, Namatame H, Taniguchi M, Ishikado M, Fujita K, Ishida S, Uchida S: Relation between the nodal and antinodal gap and critical temperature in superconducting Bi2212.

    Nat Commun 1815, 4:2013. OpenURL

  9. Hobou H, Ishida S, Fujita K, Ishikado M, Kojima KM, Eisaki H, Uchida S: Enhancement of the superconducting critical temperature in Bi2Sr2CaCu2O8+δ by controlling disorder outside CuO2 planes.

    Phys Rev B 2009, 79:064507. OpenURL

  10. Campuzano JC, Norman MR, Randeria M: Photoemission in the High Tc Superconductors. In The Physics of Superconductors. Edited by Bennemann KH, Ketterson JB. Berlin: Springer; 2004:167-273.

    [ArXiv/0209476]

    OpenURL

  11. Norman MR, Randeria M, Ding H, Campuzano JC: Phenomenology of the low-energy spectral function in high-Tc superconductors.

    Phys Rev B 1998, 57:11093-11096. Publisher Full Text OpenURL

  12. Mesot J, Norman MR, Ding H, Randeria M, Campuzano JC, Paramekanti A, Fretwell HM, Kaminski A, Takeuchi T, Yokoya T, Sato T, Takahashi T, Mochiku T, Kadowaki K: Superconducting gap anisotropy and quasiparticle interactions: a doping dependent photoemission study.

    Phys Rev Lett 1999, 83(4):840. Publisher Full Text OpenURL

  13. Angilella GGN, Sudbø A, Pucci R: Extended-wave superconductivity. Flat nodes in the gap and the low-temperature asymptotic properties of high-superconductors.

    Eur Phys J B 2000, 15(2):269-275. OpenURL

  14. Angilella GGN, Pucci R, Siringo F, Sudbø A: Sharp k-space features in the order parameter within the interlayer pair-tunneling mechanism of high-Tc superconductivity.

    Phys Rev B 1999, 59:1339-1353. Publisher Full Text OpenURL

  15. Tacon ML, Sacuto A, Georges A, Kotliar G, Gallais Y, Colson D, Forget A: Two energy scales and two distinct quasiparticle dynamics in the superconducting state of underdoped cuprates.

    Nat Phys 2006, 2:537-543. Publisher Full Text OpenURL

  16. Alldredge JW, Lee J, McElroy K, Wang M, Fujita K, Kohsaka Y, Taylor C, Eisaki H, Uchida S, Hirschfeld PJ, Davis JC: Evolution of the electronic excitation spectrum with strongly diminishing hole density in superconducting Bi2Sr2CaCu2O8+δ.

    Nat Phys 2008, 4:319-326. Publisher Full Text OpenURL

  17. Lee WS, Vishik IM, Tanaka K, Lu DH, Sasagawa T, Nagaosa N, Devereaux TP, Hussain Z, Shen ZX: Abrupt onset of a second energy gap at the superconducting transition of underdoped Bi2212.

    Nature 2007, 450:81-84. PubMed Abstract | Publisher Full Text OpenURL

  18. Pushp A, Parker CV, Pasupathy AN, Gomes KK, Ono S, Wen J, Xu Z, Gu G, Yazdani A: Extending universal nodal excitations optimizes superconductivity in Bi2Sr2CaCu2O8+δ.

    Science 2009, 324(5935):1689-1693. PubMed Abstract | Publisher Full Text OpenURL

  19. Tanaka K, Lee WS, Lu DH, Fujimori A, Fujii T, Risdiana , Terasaki I, Scalapino DJ, Devereaux TP, Hussain Z, Shen ZX: Distinct Fermi-momentumdependent energy gaps in deeply underdoped Bi2212.

    Science 2006, 314(5807):1910-1913. PubMed Abstract | Publisher Full Text OpenURL

  20. Emery VJ, Kivelson SA: Importance of phase fluctuations in superconductors with small superfluid density.

    Nature 1995, 374(6521):434-437. Publisher Full Text OpenURL

  21. Anukool W, Barakat S, Panagopoulos C, Cooper JR: Effect of hole doping on the London penetration depth in Bi2Sr1.85CaCu2O8+δ and Bi2.1Sr1.9Ca0.85Y0.15Cu2O8+δ.

    Phys Rev B 2009, 80:024516. OpenURL

  22. Tallon JL, Loram JW, Cooper JR, Panagopoulos C, Bernhard C: Superfluid density in cuprate high-Tc superconductors: a new paradigm.

    Phys Rev B 2003, 68:180501. OpenURL

  23. Schrieffer JR: Theory of Superconductivity. New York: Addison-Wesley; 1964. OpenURL

  24. Chen Q, Kosztin I, Boldizsár J, Levin K: Pairing fluctuation theory of superconducting properties in underdoped to overdoped cuprates.

    Phys Rev Lett 1998, 81:4708-4711. Publisher Full Text OpenURL

  25. Chien CC, He Y, Chen Q, Levin K: Two-energy-gap preformed-pair scenario for cuprate superconductors: implications for angle-resolved photoemission spectroscopy.

    Phys Rev B 2009, 79:214527. OpenURL

  26. Mahan GD: Many-Particle Physics. New York: Plenum; 1981. OpenURL

  27. Vishik IM, Lee WS, Schmitt F, Moritz B, Sasagawa T, Uchida S, Fujita K, Ishida S, Zhang C, Devereaux TP, Shen ZX: Doping-dependent nodal Fermi velocity of the high-temperature superconductor Bi2Sr2CaCu2O8+δ revealed using high-resolution angle-resolved photoemission spectroscopy.

    Phys Rev Lett 2010, 104:207002. PubMed Abstract | Publisher Full Text OpenURL

  28. Johnston S, Vishik IM, Lee WS, Schmitt F, Uchida S, Fujita K, Ishida S, Nagaosa N, Shen ZX, Devereaux TP: Evidence for the importance of extended Coulomb interactions and forward scattering in cuprate superconductors.

    Phys Rev Lett 2012, 108:166404. PubMed Abstract | Publisher Full Text OpenURL

  29. Rameau JD, Yang HB, Gu GD, Johnson PD: Coupling of low-energy electrons in the optimally doped Bi2Sr2CaCu2O8+δ superconductor to an optical phonon mode.

    Phys Rev B 2009, 80:184513. OpenURL

  30. Kovaleva NN, Boris AV, Holden T, Ulrich C, Liang B, Lin CT, Keimer B, Bernhard C, Tallon JL, Munzar D, Stoneham AM: c-axis lattice dynamics in Bi-based cuprate superconductors.

    Phys Rev B 2004, 69:054511. OpenURL

  31. Kulić ML: Interplay of electron-phonon interaction and strong correlations: the possible way to high-temperature superconductivity.

    Phys Rep 2000, 338:1-264. Publisher Full Text OpenURL

  32. Hong SH, Bok JM, Choi HY, Zhang W, He J, Zhou XJ: Low energy kink induced by off-plane impurities in BSCCO superconductors.

    arXiv/1306.3731 OpenURL

  33. Devereaux TP, Cuk T, Shen ZX, Nagaosa N: Anisotropic electron-phonon interaction in the cuprates.

    Phys Rev Lett 2004, 93:117004. PubMed Abstract | Publisher Full Text OpenURL

  34. McElroy K, Lee D, Hoffman J, Lang K, Lee J, Hudson E, Eisaki H, Uchida S, Davis J: Coincidence of checkerboard charge order and antinodal state decoherence in strongly underdoped superconducting Bi2Sr2CaCu2O8+δ.

    Phys Rev Lett 2005, 94(19):197005. PubMed Abstract | Publisher Full Text OpenURL

  35. Giustino F, Cohen ML, Louie SG: Small phonon contribution to the photoemission kink in the copper oxide superconductors.

    Nature 2008, 452(7190):975-978. PubMed Abstract | Publisher Full Text OpenURL