arxiv:cond-mat/ v1 [cond-mat.supr-con] 21 Jul 2004

Transcription

1 Influence of the Third Dimension of Quasi-Two-Dimensional Cuprate Superconductors on Angle-Resolved Photoemission Spectra arxiv:cond-mat/47555v [cond-mat.supr-con] 2 Jul 24 A. Bansil, M. Lindroos,2, S. Sahrakorpi, and R.S. Markiewicz Physics Department, Northeastern University, Boston, Massachusetts 25, USA 2 Institute of Physics, Tampere University of Technology, P.O. Box 692, 33 Tampere, Finland (Dated: December 2, 23) Angle-resolved photoemission spectroscopy (ARPES) presents significant simplications in analyzing strictly two-dimensional (2D) materials, but even the most anisotropic physical systems display some residual three-dimensionality. Here we demonstrate how this third dimension manifests itself in ARPES spectra of quasi-2d materials by considering the example of the cuprate Bi 2Sr 2CaCu 2O 8 (Bi222). The intercell, interlayer hopping, which is responsible for k z-dispersion of the bands, is found to induce an irreducible broadening to the ARPES lineshapes with a characteristic dependence on the in-plane momentum k. Our study suggests that ARPES lineshapes can provide a direct spectroscopic window for establishing the existence of coherent c-axis conductivity in a material via the detection of this new broadening mechanism, and bears on the understanding of 2D to 3D crossover and pseudogap and stripe physics in novel materials through ARPES experiments. PACS numbers: i, 7.8.+y, Hs Angle-resolved photoemission spectroscopy (ARPES) has been applied extensively in the recent years for investigating the electronic and quasiparticle properties of high-temperature superconductors. Much of the existing ARPES work on the cuprates and other quasi-2d materials implicitly assumes perfect two-dimensionality, ignoring the effects of dispersion in the third dimension. As energy resolutions of the order of a few mev s have now become possible in the state-of-the-art ARPES instrumentation, it is natural to ask the question: How will k z -dispersion play out in the ARPES spectra of the quasi-2d materials? An obvious answer is that spectral peaks will undergo shifts with photon energy albeit small much like in the 3D case. 3,4,5 Our analysis however reveals a very different scenario in that insofar as ARPES response to k z -dispersion is concerned, quasi-2d systems differ fundamentally from their 3D counterparts. We find that in the presence of typical final state dampings of the order of ev s, initial state dispersions over much smaller energy scales will not appear as energy shifts with hν, but will instead induce an irreducible linewidth in the spectral peaks with a characteristic k dependence. Our results bear on a variety of important issues in cuprate physics and give insight into a number of puzzling features of ARPES in the cuprates. Since our new broadening mechanism does not have its origin in a scattering process, it explains why the lifetimes derived from ARPES spectra are generally found to be shorter than those obtained from other experiments. Moreover, Fermi surface (FS) maps obtained via ARPES often display broad patches of spectral weight rather than well-defined FS imprints, particularly near the antinodal point; such patches occur naturally in our calculations. Indeed, much of the pseudogap phenomena in particular the lack of well-defined quasiparticles near (π, ) will need to be reevaluated in light of the present findings. In this vein, these results also impact the analysis of stripe physics through ARPES experiments. Finally, detection of the k z related linewidths in the ARPES spectra offers a unique spectroscopic window for establishing the existence of coherent c-axis conductivity and intercell coupling in a system; this important intrinsic property is not accessible directly through other techniques. Specifically, we focus in this article on the tetragonal body-centered Bi 2 Sr 2 CaCu 2 O 8 (Bi222) compound, which has been a workhorse of ARPES studies. Bi222 is a nearly 2D material with two CuO 2 -layers in the primitive unit cell. The intracell interaction between the two CuO 2 -layers, spaced a relatively short distance of 3.2 Å apart, results in the well-known bilayer splitting,6,7,8,9,. On the other hand, the intercell coupling between the bilayer slabs in different unit cells is expected to be smaller due to the larger intercell Cu-Cu distance of 2 Å resulting in weak but non-vanishing k z -dispersion. The intercell coupling will be even more pronounced in other high-tc s since Bi222 presents one of the longest c-axes in the cuprate family. Concerning computational details, the band structure for Bi222 was obtained within the local-densityapproximation (LDA) by using the well-established Green s function methodology. The crystal potential used is the same as that employed in our previous studies of Bi222 and involves 3 atoms per conventional unit cell 6,7,2, yielding good agreement with the experimentally observed Fermi surface (FS),3,4,5,6. ARPES intensities have been computed within the one-step photoemission formalism; see, Refs. 6 and 7 for details. Figure (a), which considers the familiar antibonding (A) and bonding (B) bands in Bi222, shows that the k -dispersion depends strongly on k z as the associated bands wander over the two sets of shaded areas. The k z - dispersion (i.e. the vertical width of the shaded areas) displays a striking dependence on k and nearly vanishes at the antinodal point k = (π/a, ). A clear bilayer splitting between the A and B bands is seen at all k- values, except at k = k =.2(2π/a) for k z = (solid

2 2 Binding energy [ev] BE [ev].5 k [2π/a] A B c) BE [ev] a). k [2π/c] z k [2π/a] b) k FIG. : (a): Calculated first-principles k -dispersion in Bi222 along []-direction for the B and A bands at three different k z values (in units of 2π/c): k z = (solid lines); k z =.5 (dashed), and k z = (dash-dotted). Shading denotes the regions over which the bands wander as a function of k z. (b): k z-dispersion of the B band at five different k values ranging from to π/a. (c): Same as (a), except that these results are based on a tight binding formalism. lines), where A and B bands touch. This level crossing leads to the anomalous k z -dispersion shown in Fig. (b) for the B band: The shape changes from being described approximately by the form sin 2 (k z c/4) for k =, to sin(k z c/4) at k, and finally back to sin2 (k z c/4) for k a/(2π) =.5. The antibonding band behaves similarly (not shown for brevity). The complex behavior of the bilayer splitting and the k z -dispersion shown in the first-principles computations of Figs. (a) and (b) can be modeled reasonably well by a relatively simple tightbinding (TB) Hamiltonian as seen from Fig. (c). We return below to discuss this point further. We next turn our attention to how k z -dispersion will manifest itself in ARPES spectra. The key insight in this regard is provided by the simulations 7 of Fig. 2, which consider the evolution of the lineshape of the energy distribution curves (EDCs) as the final state broadening introduced via the imaginary part of the final state self-energy 8, Σ f, is increased from a very small value in (a) to a realistic value in (c). In (a), the position of the spectral peak undergoes the familiar shift as hν varies. This shift results from the fact that in the photoexcitation process k must remain unchanged in transmission of the electron across the surface and the initial and final states can connect only at a specific value of k z in order to conserve energy. Note that the total shift in the peak position in the EDCs of (a) gives the size of FIG. 2: Simulated ARPES lineshapes (EDCs) in Bi222 for a series of photon energies (hν = ev) at a fixed k = (.34,.9)2π/a-point using three different values of the final state broadening given by the indicated imaginary parts of the self-energy, Σ f. In order to highlight the influence of k z-dispersion, the initial state broadening is chosen to be very small, Σ i =.2 mev. 7 the k z -dispersion of the initial state. This should not be confused with the change in hν needed to probe such a band. The change in hν is controlled by the final state dispersion, which is generally much larger than that of the initial state. Continuing to the intermediate case of Σ f =. ev in Fig. 2(b), the shift in the peak position from the bottom to the top of the initial state band is once again evident, but the lineshapes are quite different, even though the initial state damping Σ i in (b) is identical to that in (a). It is striking that some spectral intensity appears in (b) at all energies encompassed by the initial state band at every hν. This remarkable effect comes about because the energy uncertainty permitted by the width of the final state allows the photoelectron to couple with initial states off-the-energy-shell. The lineshape thus develops a new component with an hν-independent width equal to the initial state bandwidth, which rides on top of the energy conserving peak in (a). Notice also how the changes in the peak position could allow observing different FSs as one maps different values of k z. Finally, for the realistic final state width of Σ f = ev in (c), the initial state bandwidth component dominates the lineshape. The lineshape of the EDC curve is now virtually hν-independent. Despite the large final-state broadening, the energy spread of the EDC remains equal to the initial state bandwidth in k z -direction because outside of this interval, there are no initial state electrons capable of absorbing the photon. Fig. 3 elaborates on these points by considering some of the results of Fig. 2 over a much broader range of photon energies. Fig. 3(a) presents the spectra of Fig. 2(b)

3 3 FIG. 3: (Color) (a) Simulated ARPES intensities of Fig. 2(b) are depicted over a wide photon energy range of 9 27 ev for Σ f =. ev. Reds denote highs. (b) Final state band structure as a function of k z. (c) Same as (a) except that Σ f = ev, which corresponds to the case of Fig. 2(c). FIG. 4: Different imprints of the FS of Bi222. (a): Section of the FS at k z =. (b): Projection of the FS on to the ()-plane, where the broadening of the FS reflects the effect of k z-dispersion. (c) and (d) correspond to (a) and (b), respectively, except that the effect of adding a finite energy window of ±3 mev is included in the computations. for Σ f =. ev for hν 9 27 ev. Different colored bands here represent the same initial states being excited to different final states. Comparing Figs. 3(a) and 3(b), we see that the peaks of Fig. 3(a) more or less follow the final state band structure of Fig. 3(b) [The agreement is not expected to be perfect of course, since the initial state dispersion is not entirely negligible]. In Fig. 3(c), when a more realistic final state width of Σ f = ev is assumed, much of the structure is lost. The large variations in colors (or intensities) observed in Figs. 3(a) and 3(c) are the consequence of the well-known k and hν -dependency of the ARPES matrix element.,6,7,5,9 The preceding discussion of Figs. 2 and 3 makes it obvious that in a quasi-2d system, k z -dispersion leads to an irreducible linewidth in the ARPES spectra, which cannot be resolved by changing photon frequency. This effect will also appear in the FS maps observed by measuring emission from E F in the (k x, k y )-plane. Fig. 4 illustrates how this plays out. The standard bilayer split FS in the k z = plane, given in Fig. 4(a), is the FS usually thought to be measured in ARPES. In Fig. 4(b) the full 3D FS has been projected on to the ()-plane by collecting individual FS cuts corresponding to different k z values. The k z -dispersion is now seen to introduce an effective broadening to the FS imprint. In a perfect 2D-system, the maps in (b) and (a) would be identical. The effect of using a finite energy window of ±3 mev around E F in the computations is shown in Figs. 4(c) and 4(d). These simulations indicate the influence of a finite experimental resolution on the results. 3,2 Note that even with the window of ±3 mev in (d), the broadening effect of k z - dispersion is not washed out. In fact, the broadening is somewhat enhanced, especially near the antinodal point due to the contribution of the flat bands related to the van Hove singularity (VHS). It is also interesting to consider changes in the linewidth in momentum, k, as one moves away from the Fermi level. In general, our simulations indicate that k increases with increasing binding energy (BE), due mainly to the flattening of bands and the concomitant reduction in the band velocity. 4 In any event, the broadening as a function of BE is neither simply quadratic nor exponential. Furthermore, considering the the effect of varying Σ i at a general k -point, we have found that as the initial state damping due to intrinsic scattering mechanisms (simulated via the value of Σ i ) decreases, the linewidth E becomes increasingly dominated by the irreducible width associated with k z -dispersion. Along the antinodal direction this dispersion is negligible, and E and Σ i are related linearly.2 We return now to comment briefly on the bands of Figs. (a) and (b). Insight into how the complex k - and k z -dependencies of these bands reflect intercell as well as intracell hopping effects can be gained by modeling these bands within the TB framework. In the absence of intercell coupling, the conventional bilayer splitting possesses the form t bi = t z (c x c y ) 2, with c i = cos(k i a), i = x, y. The intercell coupling in the cuprates may be included in the one-band model with dispersion 22 ǫ k = 2t(c x + c y ) 4t c x c y T z (k, s z )[(c x c y ) 2 /4 + a ], where s z = sin(k z c/4) and a =.6 corrects for the finite bilayer splitting found at k =. The form of T z depends on the particular cuprate considered. While Ref. 22 lists seven inequivalent interlayer hopping

4 4 parameters, we find that we can describe the dispersion reasonably by including only the dominant contribution associated with hopping between Cu 4s levels by introducing T z = ± (t z t z )2 + 4t z t z s2 z, where plus (minus) sign refers to the bonding (antibonding) solution. t z is a constant associated with intracell interlayer hoppings, and t z = 4t zcos(k x a/2)cos(k y a/2) accounts for intercell hopping. The extra angular dependence in this term arises because for intercell hopping the two CuO 2 -planes are offset, so that one Cu atom sits above an empty cell. Note that for t z = one obtains the simple bilayer splitting t bi with no k z -dispersion. TB bands of Fig. (c) assume (in ev): t =.42, t =.2, t z =.25, and t z =.4. A comparison of the first principles and TB bands in Figs. (a) and (c), respectively, shows that our TB model reproduces the bilayer splitting, the collapse of the k z -dispersion around (π/a, ), and the anomalous dispersions associated with level crossing to a reasonably good degree, some minor discrepancies notwithstanding. We emphasize that t z, which provides intercell coupling, is essential for obtaining coherent c-axis conductivity. t z thus controls the intrinsic resistivity anisotropy and gives a measure for discriminating between coherent and incoherent c-axis hopping. 23 The lifetimes in the cuprates extracted from ARPES are considerably smaller than optical lifetimes or those deduced from transport or tunneling measurements 24. A part of this difference may be explained to be the result of an unresolved bilayer splitting 25. The remaining anomalous broadening is often interpreted 26 in terms of small-angle scattering, which does not contribute to transport, but its origin remains obscure 27. The present analysis indicates that part of the broadening of ARPES lines has its origin in the k z - dispersion, unrelated to any scattering mechanism. Our predicted effects of k z -dispersion on the linewidths in Bi222 should be resolvable with the power of currently available high resolution ARPES instrumentation. 28 The size of our novel line broadening mechanism will be larger in materials with greater 3D character (e.g. YBCO, LSCO, NCCO, ruthnates, manganites, etc.), although we expect the results of Figs. -4 to provide at least a qualitative handle in quasi-2d materials more generally. Notably, coherent 3D coupling has recently been reported in a Tl cuprate 29. Additionally, a similar k z -induced broadening may be expected to have an impact on other experimental probes, such as Raman scattering 3. In conclusion, we have demonstrated that residual k z - dispersion in quasi-2d materials will induce an irreducible linewidth in ARPES spectra. This intrinsic linewidth offers a new spectroscopic window for understanding inas well as out-of-plane scattering mechanisms and the nature of the 2D to 3D crossover in the cuprates. This highly anisotropic line broadening mechanism, which has not been recognized previously and is unrelated to 2D physics, indicates that the existing analysis of stripe and pseudogap physics based on ARPES spectra should be reexamined. Our study shows how ARPES can be extended to unravel the hidden third dimension of energy bands and Fermi surfaces in quasi-2d systems with wideranging consequences for understanding the nature of electronic states in many novel materials. Acknowledgments This work is supported by the US Department of Energy contract DE-AC3-76SF98, and benefited from the allocation of supercomputer time at NERSC, Northeastern University s Advanced Scientific Computation Center (ASCC), and the Institute of Advanced Computing (IAC), Tampere. One of us (S.S.) acknowledges Suomen Akatemia and Vilho, Yrjö ja Kalle Väisälän Rahasto for financial support. A. Damascelli et al., Rev. Mod. Phys. 75, 473 (23). 2 Throughout this article, z-axis is defined to be normal to the Cu-O planes. 3 J.E. Inglesfield and E.W. Plummer on S.D. Kevan (ed.), Angle Resolved Photoemission, Theory and Current Applications, Elsevier N.V. Smith et al., Phys. Rev. B 47, 5476 (993). 5 M. Lindroos et al., Phys. Rev. Lett. 77, 2985 (996). 6 A. Bansil et al., Phys. Rev. Lett. 83, 554 (999). 7 M. Lindroos et al., Phys. Rev. B 65, 5454 (22). 8 D.L. Feng et al., Phys. Rev. Lett. 86, 555 (2). 9 Y.-D. Chuang et al., Phys. Rev. Lett. 87, 72 (2). P. V. Bogdanov et al., Phys. Rev. B 64, 855(R) (2). See, e.g., A. Bansil et al., Phys. Rev. B 6, 3396 (999), and references therein. 2 The well-known Bi pockets in the vicinity of the M-point have been lifted above the Fermi energy by modifying the Bi and O potentials, consistent with experimental data. 3 It should be kept in mind that LDA calculations systematically overestimate band dispersions, typically by a factor of about A. Bansil et al., J. Phys. Chem. Solids 63, 275 (22). 5 M. C. Asensio et al., Phys. Rev. B 67, 459 (23). 6 Y.-D. Chuang et al., Phys. Rev. B 69, 9455 (24). 7 The spectra have not been convoluted with the Fermi function throughout this article. 8 The imaginary part of the self-energy Σ is proportional to the inverse of the lifetime of the (quasi)particle in that state. 9 S. Sahrakorpi et al., Phys. Rev. B 68, (23). 2 The ±3 mev window used would simulate an effective experimental resolution of about ±8 mev [see Ref. 3]. 2 The situation at the nodal point is more complicated, due to a small residual bilayer splitting. See also A.A. Kordyuk et al., cond-mat/337 (unpublished). 22 O.K. Andersen et al., J. Phys. Chem. Solids 56, 573

5 5 (995). 23 Of course, in Bi222, the anomalous temperature dependence of ρ c strongly hints that it is incoherent. 24 B.W. Hoogenboom et al., Phys. Rev. B 67, (23). 25 P.V. Bogdanov et al., Phys. Rev. Lett. 89, 672 (22) 26 E. Abrahams et al., Proc. Natl. Acad. Sci. 97, 574 (2). 27 A.J. Millis et al., Phys. Rev. B 67, 2457 (23). 28 Note from Fig. that a typical broadening is a significant fraction of the full bilayer splitting, which is experimentally found to be about mev [recall Ref. 3]. 29 N.E. Hussey et al., Nature 425, 84 (23). 3 T.P. Devereaux et al., Phys. Rev. B 59, 64 (999).