The Origin and Non-quasiparticle Nature of Fermi Arcs in Bi2Sr2CaCu2O8+δ

T.  J.  Reber ,  N.  C.  Plumb ,  Z.  Sun ,  Y.  Cao ,  Q.  Wang ,  K.  McElroy ,  H.  Iwasawa ,  M. 1 1 1 1 1 1 2 Arita , J. S. Wen , Z. J. Xu , G. Gu , Y. Yoshida , H. Eisaki , Y. Aiura , and D. S.  Dessau 2 3 3 3 4 4 4 1

1 Dept. of Physics, University of Colorado, Boulder, 80309-0390, USA

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

3 Condensed  Matter  Physics  and  Materials  Science  Department,  Brookhaven  National  Labs, Upton, New York, 11973 USA

4 AIST Tsukuba Central 2, 1-1-1 Umezono, Tsukuba, Ibaraki 3058568, Japan

A Fermi arc 1,2 is  a  disconnected segment of a Fermi surface observed in the pseudogap phase 3,4 of cuprate  superconductors.  This  simple  description  belies  the  fundamental inconsistency  in  the  physics  of  Fermi  arcs,  specifically  that  such  segments  violate  the topological integrity of the band . Efforts to resolve this contradiction of experiment and 5 theory  have  focused  on  connecting  the  ends  of  the  Fermi  arc  back  on  itself  to  form  a pocket,  with  limited  and  controversial  success 6,7,8,9 .  Here  we  show  the  Fermi  arc,  while composed of real spectral weight, lacks the quasiparticles to be a true Fermi surface . To 5 reach  this  conclusion  we  developed  a  new  photoemission-based  technique  that  directly probes  the  interplay  of  pair-forming  and  pair-breaking  processes  with  unprecedented precision. We find the spectral weight composing the Fermi arc is shifted from the gap edge to the Fermi energy by pair-breaking processes 10 .  While real, this weight does not form a true Fermi surface, because the quasiparticles, though significantly broadened, remain at the  gap  edge.  This  non-quasiparticle  weight  may  account  for  much  of  the  unexplained behavior of the pseudogap phase of the cuprates.

In  a  solid  the  behavior  of the electrons is most fully described in terms of the electron Green's function, the poles of which map the energy vs. momentum dependence of the electronic quasiparticles (the dressed electronic states) . The locus of poles at the Fermi energy, EF, defines 5 the Fermi surface of the material, from which almost all of the electronic properties of a material emanate.    For  a  single  continuous  band,  the  Fermi  surface  should  form  a  continuous  loop. Consequently, the broken segments of Fermi surface known as Fermi arcs apparently require a major rethinking of some of the basic tenets of condensed matter physics.

One approach to resolving this problem is completely discarding the notion of electron quasiparticles, and with it almost all of the understanding of solids built up from generations of condensed matter physicists.  Much support for this line of reasoning came from angle resolved photoemission  spectroscopy  (ARPES),  which  is  unique  in  its  ability  to  directly  probe  the electronic  excitations  as  a  function  of  energy  and  momentum,  i.e.  the  quasiparticles.    Earlier ARPES studies found that the ARPES peaks were either anomalously broad or vanishingly weak 11,12,13 especially  in  the  underdoped  'pseudogap'  regime  of  the  cuprates  -  aspects  that  were widely taken as evidence for the lack of electron-quasiparticles. The recent introduction of laser

and low energy ARPES 14 made tremendous advances in the peak sharpness and spectral weight at EF, handicapping this line of argument. We for example now see sharp nodal quasiparticle-like peaks for all doping  levels of Bi2Sr2CaCu2O8+δ (Bi2212) studied (down  to moderately underdoped x=.10 samples - see Fig. 3b black), in contrast to recent studies which were unable to observe peaks below x=0.20 13 .  Therefore other ideas are needed to understand the origin of the unusual non-Fermi liquid behavior in the cuprates.  Here we show that regardless of whether the peaks are sharp enough to be considered quasiparticles, these states do not reach the Fermi energy  and  so  cannot  dominate  the  transport  and  thermodynamics  the  way  that  regular quasiparticle states normally would. This finding drops the relevance of the question of whether these states are sharp enough to be quasiparticles.  Instead, we argue that a more relevant issue for understanding the nature of the non-Fermi liquid behavior is whether these excitations reach the Fermi energy at all.

In  addition  to  its  ability  to  measure  quasiparticle  peaks,  the  momentum-resolution  of ARPES makes it a preferred tool to study gaps and Fermi surfaces, and therefore Fermi arcs. Despite  its  great  power  and  directness,  ARPES  has  lacked  appropriate  quantitative  analysis techniques, which can also significantly affect the qualitative picture that emerges. For example, gap  values  have  been  predominantly  measured  using  the  approximate  'midpoint  of  leading edge' 15 or  'peak  separation  of  symmetrized spectra' 1,2 techniques, each of which is known to fail  in  many  limits.    Phenomenological  models 16 have  also  been  used  to  fit  gapped  energy distribution curves (EDCs) 17 , but since the EDC lineshape is not yet understood 18 , these fittings are not obviously better than the approximate measures.

To  avoid  the  difficulties  inherent  in  the  traditional  techniques  of  analyzing  ARPES spectra,  we  have  developed  and  employed  a  new  method,  the  tomographic  density  of  states (TDoS),  the  creation  of  which  is  illustrated  in  Fig  1,  with  more  details  in  the  supplementary materials.  Briefly,  the  TDoS  is  a  one-dimensional  momentum-sum  of  the  coherent  electronic spectral weight (black curve in D), which is then normalized to a similar but ungapped reference momentum-sum along the node (black curve in B),  with  some  similarities  to  the  proposal  by Vehkter &amp; Varma 19 . The resulting TDoS (panel F) represents the density of states and is in many ways equivalent to a typical Giaever tunneling curve of an s-wave superconductor 20 , except that it is localized to a single slice through the band structure (hence the name tomographic, meaning

sliced  or  sectioned).  This  unique  advantage  is  indispensible  when  electron  interactions  are strongly momentum dependent (e.g. in a d-wave superconductor).

To analyze the TDoS, we use the formula first proposed by Dynes to explain tunneling spectra from strongly coupled s-wave superconductors:

$$I _ { T D o S } \left ( \omega \right ) = \rho _ { _ { D y n e s } } \left ( \omega \right ) = \text{Re} \frac { \omega - i \Gamma _ { _ { T D o S } } } { \sqrt { \left ( \omega - i \Gamma _ { _ { T D o S } } \right ) ^ { 2 } - \Delta ^ { 2 } } }$$

(1)

where  ω  is  the  energy  relative  to  EF,  Γ TDoS  is  the  pair-breaking  scattering  rate  and  Δ  is  the superconducting gap 21,22 . This formula is essentially a simple BCS density of states with gap  broadened by the parameter  ΓTDoS, which is interpreted as the rate of pair-breaking. This formula  has been extensively tested on conventional superconductors 23 , and while initially phenomenological, it has been derived from the Eliashberg theory 22 . Since each TDoS is specific to a single location on the Fermi surface the gap magnitude is single valued, allowing us to use Dynes' original form even in d-wave superconductors.  An example fit of a TDoS to (1) is shown in  panel  F,  with  the  extracted  superconducting  gap  value  =5.9±0.1  meV  and  scattering  rate  TDoS=2.6±0.1 meV. The  TDoS values extracted here are consistent with what has been measured by  tunneling 24 and  optics 25 , though  previous  ARPES-based  methods  to  extract  electronic scattering  rates  relied  on  the  EDC  or  MDC  widths  that  give  values  roughly  an  order  of magnitude larger 26 -  see for example the EDCs of panel A that have widths of order 20 meV. This  difference  arises  because  TDoS  is  only  sensitive  to  pair-breaking  interactions  while  the EDCs and MDCs are sensitive to all electronic scattering processes (see Fig SM3).

Taking  advantage  of  the  momentum  selectivity  of  the  TDoS,  Figure  2  details  the evolution of the TDoS along the Fermi surface in the near-nodal regime from an optimally doped (Tc=91K) Bi2212 sample. By fitting to the Dynes formula (eqn 1), we find the gap is linear and symmetric about the node as expected for a d-wave  gap (Fig.  2E).  In  this  near-nodal  regime,  TDoS  is  essentially  isotropic  and  of  the  scale  of  2-3  meV,  which  is  counter  to  earlier  ARPES observations in which  was large (of the order of 20-30 meV) and grew rapidly away from the node 26 .

In  Fig.  2F,  we  compare  the  TDoS  gap  measurements  with  the  EDC  leading  edge 15 technique  (data  in  panel  C)  and  symmetrized  EDC 1,2 technique  (panel  D),  all  of  which  came from the identical ARPES data sets. We show ΔTDoS is roughly the average of the 'standard' techniques away from the node. However the other methods both fail in proximity to the node. Specifically, both the EDC leading edge method and the symmetrized EDC methods have a short 'arc' near the node in which the extracted gap value is either zero or negative (which is often artificially  set  to  zero).  Directly  observing  the  sharp  cusp  of  the  d-wave  gap  so  near  the  node confirms the findings of other probes such as thermal conductivity 27 . By comparing Figs. 2e and 2f  we  see  that  the  threshold  for  the  artificial  zero  gap  in  the  conventional  techniques  is  when ΓTDoS≈Δ.  This zero gap regime is what previously would have been considered a Fermi arc, i.e. the zero gap regime would have indicated the presence of a real portion of Fermi surface . 2

To study the Fermi arc in  greater  detail  we  move  to  an  underdoped  sample  (TC=65K) with finely gridded momentum maps in the superconducting state (T=10K) and the normal (i.e. pseudogap and/or prepairing) state (T=75K) (Fig. 3). The formation of the Fermi arc is evident in the Fermi surface maps (Fig. 3A).  The symmetrized EDCs show a small (2°) Fermi arc at 10K but a much larger one at 75K, with a gap from this method (determined by the depression of weight at zero energy) only definitively present at 10° and beyond (Fig. 3B).  However the TDoS paint a different picture, showing a smooth evolution of the gap rather than a discrete change at any one angle (Fig. 3C). Our observation of a finite near-nodal gap above TC is only possible if there are pre-formed superconducting pairs in the pseudogap state 28,29 .

Assuming a simple d-wave Δ and an isotropic  ΓTDoS, we can fit the entire momentum dependence of the TDoS at once (Fig. 3D).  Both the superconducting and pseudogap states are well fit with this form.  While ΔMax has shrunk by 30% from the superconducting state to the normal state, ΓTDoS has more than doubled, completely filling in the gap for states close to the node.  These filled-in states are the source of previous observations of a Fermi arc and should be considered a manifestation of superconductive pre-pairing in the pseudogap state.  The extremely short arc in the superconducting state is most likely due to in-plane impurity scattering (e.g. Cu vacancies). 30,31 However, the growth of the arc with temperature cannot be explained by static disorder. Instead, a dynamical process must lead to the increased pair-breaking rate that causes the arc to grow as was first proposed by Norman et al. 10

A more complete temperature dependence of the Fermi arc is shown in Fig.  4A, for a similarly underdoped sample (TC=67K).  The growth of the Fermi arc with temperature  can be 2 accurately  simulated  with  an  electron  Green's  function  which  utilizes  a  tight-binding  band structure obtained from experiment, a simple d-wave gap (ΔMax=38meV), and a scattering rate  TDoS  which  is  constant  across  the  Fermi  surface  (Fig.  2E)  but  which  varies  strongly  with temperature. The values of the electron scattering rate are directly obtained from the Dynes fits up to 110K.  At higher temperatures the fits become less reliable, so we have extrapolated those values. While the gap is seen to close slightly with increasing temperature (Fig. 3D), it was held constant  in  Fig.  4  for  simplicity.  The  observed  growth  of  the  Fermi  arc  with  increased temperature can be understood with the schematic of Fig. 4B. As the dominant effect of raising temperature is to increase ΓTDoS, the threshold of ΓTDoS=Δ moves away from the node increasing the arc length.

The  notion  of discontinuous Fermi  arcs put forward  from  the previous ARPES experiments is unphysical within the context of standard condensed matter theory, because Fermi surfaces  must  be  topologically  connected.    Therefore,  significant  effort  has  been  expended  to observe  if  and  how  these  arcs  close  to  form  a  Fermi  pocket.  For  example,  recent  ARPES experiments 6,7 have  presented  evidence  for  arc  closures,  though  many  of  these  have  been  met with skepticism for reasons including potential contamination from shadow bands, superstructure bands , or extended extrapolations. With the understanding presented here, there is no physical 9 reason why the arc would need to be closed because it is not a real locus of quasiparticle states. In this light, the difficulty in observing the 'back' side of the arc is completely natural - that is, the arc is actually arc-like, as opposed to only being the front side of a small hole pocket.

High magnetic field quantum oscillation studies in the cuprates have found evidence for a small Fermi surface pocket , with these pockets potentially related to the pockets discussed in 8 some of the ARPES experiments. However, Hall effect measurements indicate that the pockets hold  electrons  instead  of  holes ,  making  it  clear  that  the  pockets  observed  in  the  quantum 8 oscillation  experiments  are  inconsistent  with  the  ARPES  Fermi  arc  with  back-side  closure. Additionally,  new  evidence  from  NMR  indicates  that  the  high-magnetic  fields  used  in  the quantum oscillation experiments cause a reconstruction of the Fermi surface that can create the small pockets observed in those experiments 32 .

Within  the  conventional  analysis  method  of  symmetrized  EDCs,  the  Fermi  arcs  are considered loci of quasiparticle states at the Fermi level that grow with increasing temperature. Instead, the above simulations and data argue that the arcs are a result of the interplay between Δ(k,T) and  TDoS (k,T) , which we find is dominated by Δ(k) and  TDoS (T). Importantly, this result indicates  that  the  arcs  are  not  composed  of  true  quasiparticles  because,  even  though  the quasiparticles exist, they are not at at E  - rather the pole locations are the gap values, Δ(k) (Fig F 4C).    Rather  than  a  true  Fermi  surface  of  quasiparticle  poles,  the  arcs  are  regions  where  real spectral weight has been scattered inside the d-wave gap, with this incoherent non-quasiparticle weight  varying  with  Fermi  surface  angle  and  temperature  (Fig.  4C).  Aside  from  the  single quasiparticle state at the node, this non-quasiparticle weight comprises the only states available to contribute to the low energy transport and thermodynamic properties.  Therefore, these states are strong candidates to explain the unusual transport, heat capacity, and other thermodynamic properties in the pseudogap state of the cuprates 33 .

Acknowledgements. We thank G. Arnold, A. Balatsky, I. Mazin, S. Golubov and M. Hermele for valuable conversations and D. H. Lu and R. G. Moore for help at SSRL. SSRL is operated by the DOE, Office of Basic Energy Sciences. ARPES experiments at the Hiroshima Synchrotron Radiation Center were performed under proposal 09-A-48.  Funding for this research was provided by DOE Grant No. DE-FG02-03ER46066 (Colorado) and DE-AC02-98CH10886 (Brookhaven) with partial support from the NSF EUV ERC and from Kakenhi (10015981 and 19340105).

## Contributions

T. R., N. P., Z. S., Y. C., Q. W., M. A. and H. I. performed the ARPES measurements.  J. W., Z. X., G. G., Y. Y. , H. E. and Y. A. grew and prepared  the samples.  T.R. analysed the ARPES data.   T.R and D.D. developed the TDoS technique and wrote the paper with suggestions and comments by N.P., Z. S. and K. M. ARPES simulations were done by T. R. D.D. is responsible for project direction, planning and infrastructure.

.

Correspondence and requests for materials should be addressed to D. D. (Dessau@Colorado.edu)

Fig 1: Creation and Fitting of the Tomographic density of states (TDoS). a) Nodal spectrum of optimally doped (TC =92K) Bi2212 at T=50K. b) (colored) 7 out of a total of 170 individual energy distribution curves (EDCs) taken along the horizontal colored lines of panel a. The green EDC corresponds to k=kF. The sum of all 170 individual EDCs gives the spectral weight curve (black). c) Slightly (6°) off nodal spectra taken along the red cut in the inset. d) EDCs and the corresponding spectral weight for the data of panel C showing the superconducting gap at EF as well as additional spectral weight above EF due to the Bogoliubov quasiparticles. e) Off-nodal spectra normalized to nodal spectral weight (color scale data) and the spectral weight curves from panels B and D (orange and red curves, respectively).  The red curve is normalized to the orange curve to create the tomographic density of states or TDoS curve (Yellow). f) Fit of the TDoS (red circles) to the Dynes tunneling formula (green line).

Fig 2: Angular dependence of the TDOS and comparison to conventional ARPES analysis techniques from an optimal TC=91K Bi2212 sample at T=50K. a) Angle dependence of ARPES spectra. b) TDoS curves vs. angle. c) EDC's at kF vs. angle. d) Symmetrized EDC's at kF vs. angle, all from the same data set. e) ΔTDoS and ΓTDoS for angles very near the node. Error bars are ±σ returned from fits to the data of panel B. f) ΔTDoS near the node compared to Δ's measured by the other two methods.  These other methods show a short but artificial 'arc' of zero gapped states near the node, while ΔTDoS extrapolates to 0 meV only at the node.

Fig 3: Comparing the TDoS and Symmetrized EDCs when determining a Fermi Arc, for an underdoped (TC=65K ) sample.  a) Spectral weight at EF both in the superconducting (T=10K, top) and pseudogap (T=75K, bottom) states. b) Symmetrized EDCs from which one would determine a finite range of gapless states at EF that grows with increasing temperature. c) Measured TDoS spectra as a function of angle away from the node. d) Two parameter fits (ΓTDoS and ΔMax) for each entire set of spectra, with parameters listed on the panels.  Δmax is the max of the d-wave gap, occurring at the antinode, with the k-dependence of the gap forced to maintain the simple d-wave form. ΓTDoS is held to a constant value throughout k-space, consistent with the results of fig 2E.

Fig 4: Simulating the observed temperature dependence of the Fermi arc from a Tc=67K underdoped Bi2Sr2CaCu2O8 sample. a) The unusual temperature-dependent arc length is seen in both the data and simulations.  The main inputs to the simulations are an experimentally determined tight binding band structure, a simple d wave gap (both below and above Tc) with independently measured d-wave gap maximum 40 meV, and a temperature dependent scattering rate ΓTDoS (shown) with units of meV.  The majority of these ΓTDoS values (&lt;110K) are also independently experimentally determined. b) Simple picture for the temperature dependence of the Fermi arc length.  The arc is approximately the region where ΓTDoS &gt; Δ, which grows rapidly with temperature as ΓTDoS grows. c) Cartoon showing change in spectral weight (red) along the Fermi surface due to increasing ΓTDOS between the superconducting and normal states as determined by the TDoS (green)

Figure 1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-5677f273-ec9c-48a8-8725-0aaed15cc28c.jpg)

Figure 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-5ef0d94e-d512-4bb0-816e-8d683ed02f00.jpg)

Figure 3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-079f0cd4-37ec-40b6-8ba2-b55ae70b5647.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-d6395474-c84d-488f-ba3e-46f4642c9d58.jpg)

## Supplemental Methods

Samples  and  experimental  equipment. The  samples  studied  in  this  letter  were  high quality Bi2Sr2CaCu2O8+  single crystals cleaved in ultra-high vacuum.  All data were taken with a Scienta R4000 analyzer using 7 eV photons.  The Scienta electron detection systems have a noticeable nonlinearity as a function of counting rate that can be up to a 30% increase 34,35 . To our knowledge this effect has so far been ignored in the analysis of ARPES spectra though  it can play  a  significant  role,  especially  near  the  Fermi  edge  where  spectral  weight  is  rapidly decreasing. We  have  calibrated and corrected for this effect. All other methods of experimentation are standard for the field.  TDoS's do not require very high angular resolution, though high energy resolution is important because we will be studying superconducting gaps of the  order  of  a  few  meV.    For  this  reason  we  use  low  energy  photons  (7  eV),  giving  an experimental resolution of approximately 4 meV FWHM.

## Data treatment- background removal and spectral weight determination.

We find it  is  very  helpful  to  isolate  the  weight  of  the  dispersive  states  from  that  of  the background  non-dispersive  states.    We  have  done  that  here  by  analyzing  each  spectrum  via momentum  distribution  curves  (MDCs).    We  fit  each  MDC  to  a  Lorentzian  with  a  constant offset, taking the area of the Lorentzian to be the dispersive spectral weight and the offset to be the background. Figure SI1a shows a sample MDC from a nodal cut at 50K (red), its fit (teal), the  background  weight  (blue)  and  the  band  weight  (green).    The  energy  dependence  of  these various spectral weights is shown in Fig SI1b, while the effect of the background removal on the overall spectrum is shown in Fig SI1c.  The BG-removed spectrum looks essentially identical to the  original  spectrum,  except  that  some  of  the  haze  around  the  dispersive  band-like  states  is missing.

## Including Energy Resolution in Fitting the TDoS,

Though small in low energy ARPES, the energy resolution is finite and for the most accurate fits we need to include it.  To properly include the energy resolution we can't simply convolve the Dynes formula with a Gaussian with the width of the energy resolution.  That method fails to consider  the  effect  of  the  rapidly  decreasing  weight  from  the  Fermi  Edge.    Instead  we  must consider how the experimental resolution affects the data prior to the normalization procedure. Consequently, we use the formula:

$$\ F i t ( \omega ) = \frac { \left ( F ( \omega ) \, \rho _ { _ { D y n e s } } ( \omega ) \right ) \otimes R ( \omega ) } { F ( \omega ) \otimes R ( \omega ) }$$

Where  F(ω)  is  the  Fermi  Function,  ρDynes(ω)  is  the  Dynes  formula,  R(ω)  is  the  experimental resolution.  The main effect of including the resolution is the deviation from the expected 'U' shape inside the gap to the 'V' shape of our TDoS.  Once the energy resolution is included the quality of the fits is excellent as shown in Fig. SI2, where we have reprinted the TDoS from Fig 2b (Fig. SI2).

## Comparing ΓTDoS to ΓMDC

The standard ARPES method for determining the scattering rate, momentum distribution curve (MDC) widths 36 , gives static values ΓMDC(  =0,T~0) nearly an order of magnitude larger than the associated values for ΓTDOS.  Note that this data is from low photon energy ARPES, which as first shown  by Koralek et al. returns ΓMDC  significantly lower than standard ARPES experiments[ 37 ].  Furthermore, while the TDoS method returns almost identical scattering rates from many samples (Fig. SI3a), the ΓMDC's from the same exact samples vary by more than a factor of 2 (Fig. SI3b).  This wide variation in MDC widths without a similar variation in the TC indicates  that  the  MDCs  (and  EDCs)  are  subject  to  scattering  processes  that  do  not  affect  the TDoS's (or the pairing). These results show that the physical processes that are responsible for

the MDC broadening ΓMDC are different from (or in addition to) those that are responsible for the TDoS broadening ΓTDoS.

To understand the difference between the ΓMDC and ΓTDoS, we plot the effect of increasing ΓTDoS on density of states using Dynes formula, in Fig. SI3c.  The dominant effect of increasing Γ is to shift weight from the peak at Δ up to EF.  Since we are restoring weight to the Fermi energy, we argue ΓTDoS represents only the scattering processes that can break pairs while ΓMDC represents all scattering processes.   For example, In a d-wave superconductor we expect pairbreaking scattering to be those scattering events with a q=k'-k which change the phase of the order parameter, whereas non-pair-breaking scattering would not change the phase and are likely related to highly prevalent forward scattering events 38 . Figure SI3d schematically illustrates both a  pair-breaking  (blue)  and  non-pair-breaking  scattering  event  (red).    We  imagine  that  out-ofplane disorder would be well screened and would predominantly cause forward scattering and not break the pairs, while in-plane disorder (e.g. Cu vacancies) would give the large q scattering necessary for pair breaking 39,40 .

## Simulating Nodal and Off-Nodal Spectra

Using  the  electron  Green's  function,  it  is  simple,  in  principle,  to  recreate  our  observed ARPES  spectra.  Including  the  Fermi  Edge,  energy  dependent  self-energy  and  experimental resolution, the ARPES spectra can be written:

$$I _ { \text{ARPES} } \left ( \omega, k \right ) = \left ( F ( \omega ) \times - \frac { 1 } { \pi } \text{Im} G ( \omega, k ) \right ) \otimes R ( \omega, k )$$

Where F(ω) is  the  Fermi  function R(ω,k) is  the  experimental  resolution  and  the  Green's function is written:

$$G ( \omega, k ) = \frac { \omega + \Sigma - \varepsilon _ { k } } { \left ( \omega + \Sigma \right ) ^ { 2 } - \left ( \varepsilon _ { k } \right ) ^ { 2 } - \Delta _ { k } ^ { 2 } }$$

where Σ is the complex self energy and εk is the bare band and Δ is the gap magnitude.   For these simulations, we used a reasonable interpretation that Im(Σ) linear and symmetric about EF with a step at 70 meV and the minimum value as determined from the MDC width (20 meV). We took the Kramers-Kronig transform of Im(Σ) to extract the Re(Σ).  We assumed a linear bare band with vBB = 2.9 eV/Å and kF = .42 1/Å.  Finally, we set T to 50K and the energy resolution at  4  meV.  For  the  nodal  (Δ=0)  case  this  simulation  does  an  excellent  job  of  recreating  the observed spectrum as shown in Fig SI4. To simulate the slightly off-nodal  spectra  we  should only have to add the presence of a small but finite Δ.  However, the addition of a finite Δ has no observable effect on the on the simulation, even though the gap is clearly resolvable in the data. We find that the minimum Im(Σ) is just too large to achieve a significant depreciation in weight at  the  Fermi  energy  for  the  off-nodal  spectrum.    We  can  achieve  the  expected  reduction  in spectral weight by drastically reducing the minimum Im(Σ) by an order of magnitude. However the  band's  width  would  be  dramatically  too  narrow  if  we  only  used  ΓTDoS.  Instead,  we  can compensate by broadening the spectrum as if the MDC width was dominated by extrinsic final state  effects  by  convolving  a  Lorentzian  in  momentum  with  the  spectrum.    The  resulting simulation  matches  the  observed  spectra  quite  well  for  both  the  gapped  and  ungapped  cases. This result further supports our conclusion that ΓTDoS more accurately describes the superconducting state scattering rates than ΓMDC.

## Creating and Simulating the Fermi Surface Maps

To create a Fermi surface map that is representative of the coherent spectral weight at EF and thus the Fermi arc, we found it necessary to remove the incoherent background and correct

for  the  photon  energy  dependent  matrix  elements.    We  removed  the  background  as  described previously for each individual angle.  For 7 eV photons, the suppression of spectral weight due to matrix  elements  is  dramatic  beyond  10°  from  the  node.    If  left  uncorrected  this  effect  would artificially shrink the observed Fermi arcs.  To compensate, we normalize the coherent spectral weight for each angle over the range -30 to -70 meV, well below the region of the gap.

To simulate the Fermi surface maps in Fig. 4, rather than the spectral function, we started

## with electron Green's function of the form

$$G ( \omega, k _ { x }, k _ { y } ) = \frac { \omega + i \Gamma _ { T D o S } - \varepsilon _ { k } } { \left ( \omega + i \Gamma _ { T D o S } \right ) ^ { 2 } - \left ( \varepsilon _ { k } \right ) ^ { 2 } - \Delta _ { k } ^ { 2 } }$$

Where ΓTDoS is the pair-breaking scattering rate,  Δk is a d-wave gap of the form

$$\Delta _ { k } = \frac { \Delta _ { 0 } } { 2 } \left ( \cos k _ { x } - \cos k _ { y } \right )$$

And εk is the band structure, a tight binding model of the form

$$\varepsilon _ { k } = u - 2 t _ { 0 } \left ( \cos k _ { x } + \cos k _ { y } \right ) - 4 t _ { 1 } \left ( \cos k _ { x } \cos k _ { y } \right ) - 2 t _ { 2 } \left ( \cos 2 k _ { x } + \cos 2 k _ { y } \right ) \\ - 4 t _ { 3 } \left ( \cos 2 k _ { x } \cos k _ { y } + \cos k _ { x } \cos 2 k _ { y } \right ) - 4 t _ { 4 } \left ( \cos 2 k _ { x } \cos 2 k _ { y } \right )$$

Extracting  just  the  imaginary  part,  we  recover  the  spectral  function  A(ω,kx,ky)  which ARPES actually measures.  To account for the large MDC widths, we broaden the spectra as previously discussed for the 2-D simulations, but along both kx and ky for the 3-D case.

## Supplemental Discussion

## Mathematical Comparison of TDoS to Normal ARPES and Tunneling

As  the  TDoS  is  a  merger  of  ARPES  data  with  tunneling  analysis,  we  find  illustrative  to directly  compare  how  ARPES,  tunneling  and  the  TDoS  measure  the  spectral  function A(k,ω). Rather than the usual kx and k y ,  we define two orthogonal momenta, kll and k  ,  where kll is  the

momentum parallel to the Fermi surface while k  is normal to the Fermi surface, though still in the xy plane ( kz would be the out-of plane direction). A(k,  ) is the spectral function, M 2 (k,  ) is the ARPES or tunneling matrix element, and f(  ) is the Fermi function.

The intensity measured by ARPES (ignoring experimental resolution) is:

$$I ( k, \omega ) _ { \text{ABS} } = A ( k, \omega ) M ^ { 2 } ( k, \omega ) f ( \omega )$$

While  (k,  ) ARPES contains  the  most  detail,  the  pairing  processes  can  be  swamped  by  nonthermodynamic  scattering  processes  like  forward  scattering.    Tunneling  experiments  integrate over all momenta such that:

$$I ( \omega ) _ { \iota u m e l i n g } = \int _ { k _ { \perp } } \int _ { k _ { l l } } A ( k, \omega ) M ^ { 2 } ( k, \omega ) d ^ { 2 } k$$

While  tunneling  is  sensitive  to  the  pairing  processes , I(  ) tunneling is  complicated  by  strong variation along kll such as to the d-wave superconducting gap, the van-Hove singularity near the antinode, etc.  However to create the TDoS we integrate over the momentum, k  :

$$I ( k _ { i l, } \omega ) _ { T D o S } = \int _ { k _ { \perp } } A ( k, \omega ) M ^ { 2 } ( k \, \omega ) f ( \omega ) d k _ { \perp }$$

By normalizing the spectra properly we obtain the simple functional as shown below,

$$\cong \int _ { k _ { \perp } } A ( k, \omega ) d k _ { \perp }$$

This form allows us to probe the thermodynamic pair-forming and pair-breaking processes with the  use  of  well  developed  tunneling  analysis  techniques  for  a  simple  (i.e.  s-wave)  gap 41 ,  even though the gap varies as a function of kll due to the d-wave nature of the order parameter.

Fig SI1:  Removing the Background, a) Sample MDC (red) 100 meV from EF, its fit (teal), background  weight  (blue),  and  band  weight  (Green).  The  MDC  at  EF  (dashed)  is  shown  for comparison. b) Energy  dependence  of  spectral  weights  for  a  nodal  spectrum  at  50K. c) Comparison of the original nodal spectrum and with the background removed.

Fig SI2: Experimental TDoS spectra and Fits. TDoS and fits as a function of angle near the node used in figure 2.

Fig SI3: Broadening mechanisms for MDC curves not present in TDoS spectra.  a) TDoS spectra at θ=6° from 3 different samples, including extracted parameters Δ and ΓTDoS. b) Nodal MDCs at EF from the same three samples as studied in panel a, showing the significant deviation in MDC widths from sample to sample. The extracted broadening parameters ΓMDC are shown in a. c) Illustration  of  how  increasing  Γ TDoS shifts  weight  from the peak to EF d) Diagram of the Fermi  surface  detailing  the  difference  between  non-pair-breaking  forward  scattering  (red)  and pair-breaking unitary scattering (blue).

## Figure SI4: Simulations of Nodal and Off-Nodal Spectra from the Electron Green's

Function. The first column shows experimental data in the superconducting state for the node (top and slightly off-node (bottom).  The second column shows simulations using Eq. S2 and S3 with a Γ=20 meV (obtained from MDC analysis).  While nodal case closely matches the experimental data, the off-nodal simulation cannot recreate the gap. Instead we must lower Γ by an order of magnitude to achieve the observed depletion of weight at EF (third column).  To recover the experimental MDC widths, we then broaden the spectra by convolving a Lorentzian with a width, Γk, in momentum (fourth column).

Figure SI1

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-a1e7421a-2d1b-40ab-9afc-86a51a3846f3.jpg)

Figure SI2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-1cc7a0b8-03e3-4020-b942-928d2dc9381f.jpg)

a

Figure SI3

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-8bc72931-011e-4067-9ebe-d300452c4389.jpg)

1

## Data

## Simulations

Figure SI4 Momentum

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0441v1-figure-84845010-4ce7-49f9-b7fe-e86a7ed6bfee.jpg)

1 Norman, M. et al. Destruction of the Fermi surface in underdoped high-TC superconductors. Nature 392 , 157(1998)

- 2 Kanigel, A. et al. Evolution of the pseudogap from Fermi arcs to the nodal liquid. Nature Physics 2 , 447 (2006).
- 3 Loeser, A. G. et al., Excitation gap in the normal state of underdoped Bi Sr CaCu2O8+δ  Science 2 2 273 , 325 (1996).
- 4 Ding, H. et al., Spectroscopic evidence for a pseudogap in the normal state of underdoped high-TC superconductors Nature 382 , 51 (1996)
- 5 Pines, P.  &amp; Nozieres, P., The Theory of Quantum Liquids , W.A. Benjamin Inc., New York (1966)
- 6 Yang, H.-B., et al.  Emergence of preformed Cooper pairs from the doped Mott insulating state in Bi2Sr2CaCu2O8+δ. Nature 456 ,77 (2008)
- 7 Meng, J. et al.  Coexistence of Fermi arcs and Fermi pockets in a highT c copper oxide superconductor. Nature (London) 462 , 335 (2009)
- 8 Doiron-Leyraud,  N.  et  al.    Quantum  oscillations  and  the  Fermi  surface  in  an  underdoped  highT c superconductor. Nature 447 , 565 (2007)
- 9 King, P.D.C.  et al, Structural origin of apparent Fermi surface pockets in angle-resolved photoemission of Bi2 Sr2-x LaxCuO6+d. Phys. Rev. Lett. 106 , 127005 (2011).
- 10 Norman, M. R., Kanigel, A. Randeria M., Chatterjee, U. &amp; Campuzano, J. C. Modeling the Fermi arc in underdoped cuprates Phys. Rev. B 76 , 174501 (2007)
- 11 Damascelli,  A.  Hussain,  Z.  &amp;  Shen,  Z-X  Angle-resolved  photoemission  studies  of  the  cuprate superconductors. Review of Modern Physics 75 , 473 (2003).
- 12 Laughlin, R.B.  Evidence for Quasipartile Decay from Photoemission in Underdoped Cuprates, Phys. Rev. Lett 79, 1726 (1997)
- 13 Fournier, D. et al, Loss of nodal quasiparticle integrity in underdoped YBa2Cu3O6+ x Nature Physics 6 , 905-911 (2010)
- 14 Koralek, J.D. et al, Laser Based Angle-Resolved Photoemission, the Sudden Approximation, and Quasiparticle-Like Spectral Peaks in Bi2Sr2CaCu2O8+d Phys. Rev. Lett. 96, 017005 (2006)
- 15 Shen, Z.-X. et al. Anomalously Large Gap Anisotropy in the a-b plane of Bi 2 Sr 2 CaCu 2 O 8+ Phys. Rev. Lett. 70 , 1553 (1993).
- 16 Norman, M. R., Randeria, M., Ding, H. &amp; Campuzano, J. C. Phenomenology of the low-energy spectral function in high-TC superconductors. Phys. Rev. B 57 , R11093-R11096 (1998).
- 17 Lee, W.S. et al.,  Abrupt onset of a second energy gap at the superconducting transition of underdoped Bi2212. Nature 450 , 81-84 (2007).
- 18 Casey, P. A., Koralek, J. D., Plumb, N. C., Dessau, D.S. &amp; Anderson, P. W.   Accurate theoretical fits to laser-excited photoemission spectra in the normal phase of high temperature superconductors. Nature Physics 4 , 210 (2008).
- 19 Vehkter, I. &amp; Varma C. M. Proposal to Determine the Spectrum of Pairing Glue in High-Temperature Superconductors, Phys. Rev. Lett. 90 , 237003 (2003)
- 20 Giaever, I., Hart, H. R. , &amp; Mergele, K.  Tunneling into superconductors at temperatures below 1 degree K Phys. Rev. 126 , 941-948 (1962).
- 21 Dynes,  R.  C.,  Narayanamurti,  V.,  Garno,  J.P.,  Direct  Measurement  of  Quasiparticle-Lifetime Broadening in a Strong-Coupled Superconductor. Phys. Rev. Lett. 41 ,1509-1512 (1978)
- 22 Mikhailovsky, A. A., Shulga, S. V., Karakozov, A. E., Dolgov, O. V., &amp; Maksimov, E. G.   Thermal pair-breaking in superconductors with strong electron-phonon interaction. Solid State Communications 80 ,7 (1991)
- 23 Wolf, E.L.  Principles of Electron Tunneling Spectroscopy, Oxford University Press, USA (1989)
- 24 Pasupathy, A. et al. Electronic Origin of the Inhomogeneous Pairing Interaction in the High-Tc Superconductor Bi2Sr2CaCu2O8+δ. Science 320 , 196 (2008).

25 Hwang, J., Timusk, T. &amp; Gu G. D.  High-transition-temperature superconductivity in the absence of the magnetic-resonance mode. Nature 427 , 6976 (2004)

- 26 Valla, T. et al. Temperature Dependent Scattering Rates at the Fermi Surface of Optimally Doped Bi2Sr2CaCu2 O 8+ δ . Phys. Rev. Lett. 85 , 828 (2000).
- 27 Sutherland,  M.  et  al.  Thermal  conductivity  across  the  phase  diagram  of  cuprates:  Low-energy quasiparticles and doping dependence of the superconducting gap. Phys. Rev. B 67 , 174520 (2003).
- 28 Emery,  V.J.  &amp;  Kivelson,  S.A.  Importance  of  phase  fluctuations  in  superconductors  with  small superfluid density. Nature 374 , 434(1995)
- 29 Y. Wang, L. Li, and N. P. Ong, Nernst Effect in High TC Superconductors Phys. Rev. B 73 , 024510 (2006)
- 30 Fehrenbacher, R. Effect of nonmagnetic impurities on the gap of a dx^2-y^2 superconductor as seen by angle-resolved photo-emission Phys. Rev. B 54 6632 (1996)
- 31 Haas, S., Balatsky, A. V., Sigrist, M. &amp; Rice, T. M. Extended gapless regions in disordered dx^2-y^2 wave superconductors Phys. Rev. B 56 5108 (1997)
- 32 Wu, T. et al., Magnetic-field-induced charge-stripe order in the high-temperature superconductor YBa2Cu3Oy.  Nature 477 , 191 (2011)
- 33 Varma, C.M., Littlewood, P. B., Schmitt-Rink, S., Abrahams, E.  &amp; Ruckenstein, A. E. Phenomenology of the normal state of Cu-O high-temperature superconductors. Phys. Rev Lett. 63 , 1996 (1989)
- 34 Kay, A. W. et al. Multiatom resonant photoemission Phys. Rev. B 63 , 115119 (2001)
- 35 Manella, N. et al. Correction of Non-linearity Effects in Detectors for Electron Spectroscopy Jour. of Elec. Spec. and Rel. Phen. 141 , 45 (2004)
- 36 T. Valla et al., Evidence for quantum critical behavior in the optimally doped cuprate Bi Sr CaCu2O8+δ 2 2 Science 285 , 2110 (1999)
- 37 Koralek, J.D. et al, Laser Based Angle-Resolved Photoemission, the Sudden Approximation, and Quasiparticle-Like Spectral Peaks in Bi2Sr2CaCu2O8+d Phys. Rev. Lett. 96, 017005 (2006)
- 38 Dahm, T., Hirschfeld, P.J., Scalapino, D. J., &amp; Zhu, L. Nodal quasiparticles in the cuprate superconductors Phys. Rev. B 72 , 214512 (2005)
- 39 Haas, S., Balatsky, A. V., Sigrist, M. &amp; Rice, T. M. Extended gapless regions in disordered dx^2-y^2 wave superconductors Phys. Rev. B 56 5108 (1997)
- 40 Fehrenbacher, R. Effect of nonmagnetic impurities on the gap of a dx^2-y^2 superconductor as seen by angle-resolved photo-emission Phys. Rev. B 54 6632 (1996)
- 41 Dynes, R. C., Narayanamurti, V., and Garno, J. P. Direct Measurement of Quasiparticle-Lifetime Broadening in a Strong-Coupled Superconductor. Phys. Rev. Lett. 41 1509 (1978)
T. J. Reber N. C. Plumb Z. Sun Y. Cao Q. Wang K. McElroy H. Iwasawa M. Arita J. S. Wen Z. J. Xu G. Gu Y. Yoshida H. Eisaki Y. Aiura D. S. Dessau
The Origin and Non-quasiparticle Nature of Fermi Arcs in Bi$_2$Sr$_2$CaCu$_2$O$_{8+δ}$
A Fermi arc is a disconnected segment of a Fermi surface observed in the
pseudogap phase of cuprate superconductors. This simple description belies the
fundamental inconsistency in the physics of Fermi arcs, specifically that such
segments violate the topological integrity of the band. Efforts to resolve this
contradiction of experiment and theory have focused on connecting the ends of
the Fermi arc back on itself to form a pocket, with limited and controversial
success. Here we show the Fermi arc, while composed of real spectral weight,
lacks the quasiparticles to be a true Fermi surface. To reach this conclusion
we developed a new photoemission-based technique that directly probes the
interplay of pair-forming and pair-breaking processes with unprecedented
precision. We find the spectral weight composing the Fermi arc is shifted from
the gap edge to the Fermi energy by pair-breaking processes. While real, this
weight does not form a true Fermi surface, because the quasiparticles, though
significantly broadened, remain at the gap edge. This non-quasiparticle weight
may account for much of the unexplained behavior of the pseudogap phase of the
cuprates.