arxiv:cond-mat/ v3 3 Jan 2006

Transcription

1 Estimating the mass of vortices in the cuprate superconductors Lorenz Bartosch a,b, Leon Balents c, and Subir Sachdev a arxiv:cond-mat/ v3 3 Jan 2006 a Department of Physics, Harvard University, Cambridge MA b Institut für Theoretische Physik, Universität Frankfurt, Postfach , Frankfurt, Germany c Department of Physics, University of California, Santa Barbara, CA Abstract We explore the experimental implications of a recent theory of the quantum dynamics of vortices in superfluids proximate to Mott insulators. The theory predicts modulations in the local density of states in the regions over which the vortices execute their quantum zero point motion. We use the spatial extent of such modulations in scanning tunnelling microscopy measurements (Hoffman et al, Science 295, 466 (2002)) on the vortex lattice of Bi 2 Sr 2 CaCu 2 O 8+δ to estimate the inertial mass of a point vortex. We discuss other, more direct, experimental signatures of the vortex dynamics. 1 Introduction It is now widely accepted that superconductivity in the cuprates is described, as in the standard Bardeen-Cooper-Schrieffer (BCS) theory, by the condensation of charge 2e Cooper pairs of electrons. However, it has also been apparent that vortices in the superconducting state are not particularly well described by BCS theory. While elementary vortices do carry the BCS flux quantum of hc/2e, the local electronic density of states in the vortex core, as measured by scanning tunnelling microscopy (STM) experiments, has not been explained naturally in the BCS framework. Central to our considerations here are the remarkable STM measurements of Hoffman et al. [1] (see also Refs. [2,3]) who observed modulations in the local density of states (LDOS) with a period of approximately 4 lattice spacings in the vicinity of each vortex core of a vortex lattice in Bi 2 Sr 2 CaCu 2 O 8+δ. This paper shall present some of the physical implications of a recent theory of two-dimensional superfluids in the vicinity of a quantum phase transition Preprint submitted to Elsevier Science 3 January 2006

2 to a Mott insulator [4,5] (see also Ref. [6]). By Mott insulator we mean here an incompressible state which is pinned to the underlying crystal lattice, with an energy gap to charged excitations. In the Mott insulator, the average number of electrons per unit cell of the crystal lattice, n MI, must be a rational number. If the Mott insulator is not fractionalized and if n MI is not an even integer, then the Mott insulator must also spontaneously break the space group symmetry of the crystal lattice so that the unit cell of the Mott insulator has an even integer number of electrons. There is evidence that the hole-doped cuprates are proximate to a Mott insulator with n MI = 7/8 [7], and such an assumption will form the basis of our analysis of the STM experiments on Bi 2 Sr 2 CaCu 2 O 8+δ. The electron number density in the superfluid state, n S, need not equal n MI and will be assumed to take arbitrary real values, but not too far from n MI. A key ingredient in our analysis will be the result that the superfluid carries a subtle quantum order, which is distinct from Landau-Ginzburg order of a Cooper pair condensate. In two dimensions, vortices are point-like excitations, and are therefore bona fide quasiparticle excitations of the superfluid. The quantum order is reflected in the wavefunction needed to describe the motion of the vortex quasiparticle. For n MI not an even integer, the low energy vortices appear in multiple degenerate flavors, and the space group symmetry of the underlying lattice is realized in a projective unitary representation that acts on this flavor space. Whenever a vortex is pinned (either individually due to impurities, or collectively in a vortex lattice), the space group symmetry is locally broken, and hence the vortex necessarily chooses a preferred orientation in its flavor space. As shown in Ref. [4], this implies the presence of modulations in the LDOS in the spatial region over which the vortex executes its quantum zero point motion [8]. The short-distance structure and period of the modulations is determined by that of the Mott insulator at density n MI, while its long-distance envelope is a measure of the amplitude of the vortex wavefunction (see Fig. 1). Consequently, the size of the region where the modulations are present is determined by the inertial mass of the vortex. Here we will show how these ideas can be made quantitatively precise, and use current experiments to obtain an estimate of the vortex mass, m v. There have been a number of theoretical discussions of m v using BCS theory [9,10,11,12,13,14,15], and they lead to the order of magnitude estimate m v m e (k F ξ) 2, where m e is the electron mass, k F is the Fermi wavevector, and ξ is the BCS coherence length. 2 Vortex equations of motion We begin with a simple model computation of the vortex dynamics. Consider a system of point vortices moving in a plane at positions r j, where j is a 2

3 Fig. 1. Schematic of the modulations in the LDOS of a vortex lattice. The short distance modulations in each vortex halo are determined by the orientation of the vortex in flavor space, as discussed in Ref. [4]. The envelope of these modulations is Ψ(r j ) 2 where Ψ is the wavefunction of the vortices, and its characteristics are computed in the present paper. label identifying the vortices. We do not explicitly identify the orientation of each vortex in flavor space, because we are interested here only in the longdistance envelope of the LDOS modulations; the flavor orientation does not affect the interactions between well-separated vortices, and so plays no role in determining the wavefunction of the vortex lattice. In a Galilean-invariant superfluid, the vortices move under the influence of the Magnus force d 2 r j m v dt = hn ( S v 2 2a 2 s (r j ) dr ) j ẑ, (1) dt where t is time, h = 2π is Planck s constant, v s (r) is the superfluid velocity at the position r, and n S /a 2 is the electron number density per unit area (a 2 is the area of a unit cell of the underlying lattice). One point of view is that the force in Eq. (1) is that obtained from classical fluid mechanics after imposing the quantization of circulation of a vortex. However, Refs. [16] emphasized the robust topological nature of the Magnus force and its connection to Berry phases, and noted that it applied not only to superfluids of bosons, but quite generally to superconductors of paired electrons. Here, we need the modification of Eq. (1) by the periodic crystal potential and the proximate Mott insulator. This was implicit in the results of Ref. [4], and we present it in more physical terms. It is useful to first rewrite Eq. (1) as m v d 2 r j dt 2 = F E(j) + F B (j), (2) where F E is the first term proportional to v s and F B is the second term. Our notation here is suggestive of a dual formulation of the theory in which the vortices appear as charges, and these forces are identified as the electrical 3

4 and magnetic components. In the Galilean invariant superfluid, the values of F E and F B are tied to each other by a Galilean transformation. However, with a periodic crystal potential, this constraint no longer applies, and their values renormalize differently as we now discuss. The influence of the crystal potential on F E is simple, and replaces the number density of electrons, n S, by the superfluid density. Determining v s (r j ) as a sum of contributions from the other vortices, we obtain [17] F E (j) = 2πρ s k j r j r k r j r k 2, (3) where ρ s is the superfluid stiffness (in units of energy). It is related to the London penetration depth, λ, by where d is the interlayer spacing. ρ s = 2 c 2 d 16πe 2 λ2, (4) The modification of F B is more subtle. This term states that the vortices are charges moving in a magnetic field with n S /2 flux quanta per unit cell of the periodic crystal potential. In other words, the vortex wavefunction is obtained by diagonalizing the Hofstadter Hamiltonian which describes motion of a charged particle in the presence of a magnetic field and a periodic potential. As argued in Ref. [4], it is useful to examine this motion in terms of the deviation from the rational flux n MI /2 = p/q (p, q are relatively prime integers) associated with the proximate Mott insulator. The low energy states of the rational flux Hofstadter Hamiltonian have a q-fold degeneracy, and this constitutes the vortex flavor space noted earlier [18]. However, these vortex states describe particle motion in zero magnetic field, and only the deficit (n S n MI )/2 acts as a magnetic flux. This result is contained in the action in Eq. (2.46) of Ref. [4], which shows that the dual gauge flux fluctuates about an average flux determined by (n S n MI ). The action in Ref. [4] has a relativistic form appropriate to a system with equal numbers of vortices and anti-vortices. Here, we are interested in a system of vortices induced by an applied magnetic field, and can neglect anti-vortices; so we should work with the corresponding non-relativistic version of Eq. (2.46) of Ref. [4]. In its first-quantized version, this non-relativistic action for the vortices leads to the Lorentz force in Eq. (2) given by F B (j) = h(n S n MI ) dr j 2a 2 dt ẑ, (5) If the density of the superfluid equals the commensurate density of the Mott insulator, then F B = 0; however, F E remains non-zero because we can still have ρ s 0 in the superfluid. These distinct behaviors of F E,B constitute a 4

5 Fig. 2. Dispersion of the magnetophonon modes of the vortex lattice (here we take ω c = 0.5ω p ). key difference from Galilean-invariant superfluids. In experimental studies of vortex motion in superconductors [19], a force of the form of Eq. (5) is usually quoted in terms of a Hall drag co-efficient per unit length of the vortex line, α; Eq. (5) implies α = h(n S n MI ). (6) 2a 2 d Thus the periodic potential has significantly reduced the magnitude of α from the value nominally expected [20] by subtracting out the density of the Mott insulator. A smaller than expected α is indeed observed in the cuprates [19]. It is worth emphasizing that F B (but not F E ) is an intrinsic property of a single vortex. Moreover, we expect that, taken together, the relation Eq. (6) and the flavor degeneracy q are robust universal measures of the quantum order of a clean superconductor, independent of details of the band structure, etc. 3 Vortex lattice normal modes We solved Eqs. (2-5) in the harmonic approximation in deviations from a perfect triangular lattice of vortices at positions R j. This leads to the vortex magnetophonon modes shown in Fig 2. The computation of these modes is a generalization of other vortex oscillation modes discussed previously in superconductors [21,22], rotating superfluids [23,24,25], and, in a dual picture of charges, also to oscillations of electronic Wigner crystals in a magnetic field [26]. Quantizing these modes, we determine the mean square displacement of each vortex due to the quantum zero point motion of the vortex lattice, which 5

6 we denote u 2 rms = r j R j 2 /2. We found an excellent fit of this solution to the interpolation formula m v = a2 v 2 F(x), (7) ρ s u 4 rms where a v is the separation between nearest neighbor vortices, x α du 2 rms /, and F(x) x 2 ; (8) Eq. (8) holds only as long as the r.h.s. is positive, while F(x) = 0 for larger x (we will see below and in Fig. 3 that this apparent upper bound on x is relaxed once we allow for viscous damping). Similarly, for the frequency ω mp in Fig. 2, we obtain ω mp = ω 2 p + ω2 c with ω p = 35.45ρ su 2 rms a 2 v [F(x)]1/2 ; ω c = 27.57ρ su 2 rms xf(x). (9) a 2 v For the experiments of Ref. [1] we estimate ρ s = 12 mev [27], u rms = 20 Å [28], a = 3.83 Å, and a v = 240 Å. The overall scale for m v is determined by setting n S = n MI so that x = 0 and F(x) = 1. This yields m v 8m e and ω p 3 mev (or ν p 0.7 THz). For a more accurate determination, we need n S, for which there is considerable uncertainty e.g. for n S n MI = 0.015, we find x = 1.29, m v 3m e and ω p 5 mev. 4 Limitations We now consider the influence of a variety of effects which have been neglected in this simple computation: 4.1 Viscous drag It is conventional in models of vortex dynamics at low frequencies [29] to include a dissipative viscous drag term in the equations of motion, contributing an additional force F D (j) = ηd dr j dt, (10) to the r.h.s. of Eq. (2). There are no reliable theoretical estimates for the viscous drag co-efficient, η, for the cuprates. However, we can obtain estimates of its value from measurements of the Hall angle, θ H, which is given by [29,30] tan θ H = α/η. Harris et al. [30] observed a dramatic increase in the value tanθ H, to the value 0.85, at low T in 60 K YBa 2 Cu 3 O 6+y crystals, suggesting a small η, and weak dissipation in vortex motion. For our purposes, we 6

7 Fig. 3. Plot of the function F(x,y) which replaces F(x) in Eqs. (7,8,9) upon including viscous drag, η (y ηdu 2 rms / ). The argument x measures the Hall drag α (x α du 2 rms / ) and the Hall angle is determined by tan θ H = x/y. need the value of η for frequencies of order ω p, and not just in the d.c. limit. The very naïve expectation that η(ω) behaves like the quasiparticle microwave conductivity would suggest it decreases rapidly beyond a few tens of GHz, well below ω p [31]. Lacking solid information, we will be satisfied with an estimate of the influence of viscous drag obtained by neglecting the frequency dependence of η (a probable overestimate of its influence). The resulting corrections to Eqs. (7-9) are easily obtained (as in Ref. [11]), and can be represented by the replacement and F(x) F(x, y), where y ηdu2 rms (11) F(x, 0) = F(x). (12) The sketch of the function F(x, y) is in Fig. 3; as long as y > 0, we have F(x, y) > 0. As expected, the viscous damping decreases the estimate of the mass, and this decrease is exponential for large y, e.g. at x = 0 we have the interpolation formula F(0, y) ( y y 2 )e 3.43y. (13) 4.2 Meissner screening The interaction in Eq. (3) is screened at long distances by the supercurrents, and the intervortex coupling becomes exponentially small. This does have an important influence at small momenta in that the shear mode of the vortex lattice disperses as [21,22] k 2. However, as long as a v λ, there will not be a significant influence on u rms or ω p. 7

8 4.3 Retardation The interaction Eq. (3) is assumed to be instantaneous; in reality it is retarded by the propagation of the charged plasmon mode of the superfluid. We have estimated the corrections due to this mode in a model of superfluid layers coupled by the long-range Coulomb interaction. Briefly, we compare the energy per unit area of a phase fluctuation at the wavevector of the vortex lattice Brillouin zone boundary ( ρ s /a 2 v ) with its electrostatic energy ( 2 ωp/(e 2 2 d)); this shows that such corrections are of relative order ( 2 /m v )/(e 2 d). For the parameters above and d = 7.5 Å this is 0.009, and hence quite small. 4.4 Nodal quasiparticles We expect that nodal quasiparticles contribute to the viscous drag, and so their contribution was already included in the experimentally determined estimate of η in (i). The quasiparticle contribution to m v has been estimated [9,10] in a BCS theory of a continuum of core levels, valid in the limit k F ξ 1, and a 1/ H dependence (H is the magnetic field) was obtained. The applicability of such a model to the cuprates is questionable, but nevertheless the estimates of m v above should provisionally be assumed to be specific to the H field in the STM experiments. 4.5 Disorder We have assumed here a triangular lattice of vortices. In reality, STM experiments show significant deviations from such a structure, presumably because of an appreciable random pinning potential. This pinning potential will also modify the vortex oscillation frequencies and its mean square displacement. Both pinning and damping η tend to reduce vortex motion. For this reason, the estimates of m v above in which these effects are neglected must be regarded as upper bounds. The above considerations make it clear that new experiments on cleaner underdoped samples, along with a determination of the spatial dependence of the hole density (to specify α), are necessary to obtain a more precise value for m v ; determining the H dependence of m v will enable confrontation with theory. 8

9 5 Implications An important consequence of our theory is the emergence of ω p as a characteristic frequency of the vortex dynamics. It would therefore be valuable to have an inelastic scattering probe which can explore energy transfer on the scale of ω p, and with momentum transfer on the scale of h/a v, possibly by neutron [32] or X-ray scattering. A direct theoretical consideration of magnetoconductivity in our picture would have implications for far-infrared or THz spectroscopy, allowing comparison to existing experiments [33]; further such experiments on more underdoped samples would also be of interest. Another possibility is that the zero point motion of the vortices emerges in the spectrum of the LDOS measured by STM at an energy of order ω p. We speculate that understanding the vortex core states observed in STM studies [2,3] will require accounting for the quantum zero point motion of the vortices; it is intriguing that the measured energy of these states is quite close to our estimates of ω p. 6 Acknowledgements We thank J. Brewer, E. Demler, Ø. Fischer, M. P. A. Fisher, W. Hardy, B. Keimer, N. P. Ong, T. Senthil, G. Sawatzky, Z. Tešanović and especially J. E. Hoffman and J. C. Seamus Davis for useful discussions. This research was supported by the NSF under grants DMR (L. Balents), DMR and DMR (S.S.), the Packard Foundation (L. Balents), the Deutsche Forschungsgemeinschaft under grant BA 2263/1-1 (L. Bartosch), and the John Simon Guggenheim Memorial Foundation (S.S.). References [1] J. E. Hoffman, E. W. Hudson, K. M. Lang, V. Madhavan, S. H. Pan, H. Eisaki, S. Uchida, and J. C. Davis, Science 295, 466 (2002). [2] B. W. Hoogenboom, K. Kadowaki, B. Revaz, M. Li, Ch. Renner, and Ø. Fischer, Phys. Rev. Lett. 87, (2001); G. Levy, M. Kugler, A. A. Manuel, Ø. Fischer, and M. Li, preprint. [3] S. H. Pan, E. W. Hudson1, A. K. Gupta, K.-W. Ng, H. Eisaki, S. Uchida, and J. C. Davis, Phys. Rev. Lett. 85, 1536 (2000). [4] L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and K. Sengupta, Phys. Rev. B 71, (2005). 9

10 [5] L. Balents, L. Bartosch, A. Burkov, S. Sachdev, and K. Sengupta, Phys. Rev. B 71, (2005). [6] Z. Tešanović, Phys. Rev. Lett. 93, (2004); M. Franz, D. E. Sheehy, and Z. Tešanović, Phys. Rev. Lett. 88, (2002). [7] J. M. Tranquada, B. J. Sternlieb, J. D. Axe, Y. Nakamura, and S. Uchida, Nature 375, 561 (1995); J. M. Tranquada, H. Woo, T. G. Perring, H. Goka, G. D. Gu, G. Xu, M. Fujita, and K. Yamada, Nature 429, 534 (2004). [8] The LDOS modulations are a signal of valence-bond-solid (VBS) order that can be associated with the vortex flavors [4] (see also C. Lannert, M. P. A. Fisher, and T. Senthil, Phys. Rev. B 63, (2001)). VBS order in vortex cores was predicted in K. Park and S. Sachdev, Phys. Rev. B 64, (2001). [9] G. E. Volovik, Pisma Zh. Eksp. Teor. Fiz (1997) [JETP Lett (1997)]. [10] N. B. Kopnin, Theory of Nonequilibrium Superconductivity, Cambridge University Press, Cambridge (2001); N. B. Kopnin and V. M. Vinokur, Phys. Rev. Lett. 81, 3952 (1998). [11] G. Blatter, V. B. Geshkenbein, and V. M. Vinokur, Phys. Rev. Lett. 66, 3297 (1991); G. Blatter and B. Ivlev, Phys. Rev. Lett. 70, 2621 (1993). [12] J.-M. Duan and A. J. Leggett, Phys. Rev. Lett. 68, 1216 (1992). [13] D. M. Gaitonde and T. V. Ramakrishnan, Phys. Rev. B 56, (1997). [14] M. W. Coffey, J. Phys. A 31, 6103 (1998). [15] J. H. Han, J. S. Kim, M. J. Kim, and P. Ao, Phys. Rev. B 71, (2005). [16] F. D. M. Haldane and Y.-S. Wu, Phys. Rev. Lett. 55, 2887 (1985); P. Ao and D. J. Thouless, Phys. Rev. Lett. 70, 2158 (1993). [17] The long-range force in Eq. (3) requires that boundary conditions be chosen carefully. In the dual charge analogy, we need a uniform neutralizing background charge, and this contributes an additional force 2π 2 ρ s r j /A v to Eq. (3), where A v = 3a 2 v/2 is the area per vortex. [18] At n MI /2 = p/q, the vortex degeneracy is q in the short-range pairing models of Ref. [5]. In L. Balents and S. Sachdev, to appear, we show that a similar degeneracy can be understood in long-range pairing models of d-wave superconductors based upon the gauge theory of X.-G. Wen and P. A. Lee, Phys. Rev. Lett. 76, 503 (1996). Any distinctions between these cases have no impact on the considerations of the present paper. [19] J. M. Harris, N. P. Ong, P. Matl, R. Gagnon, L. Taillefer, T. Kimura, and K. Kitazawa, Phys. Rev. B 51, R12053 (1995). [20] N. B. Kopnin and V. E. Kravtsov, Pis ma Zh. Eksp. Teor. Fiz. 23, 631 (1976) [JETP Lett. 23, 578 (1976)]. 10

11 [21] P. G. de Gennes and J. Matricon, Rev. Mod. Phys. 36, 45 (1964). [22] A. L. Fetter, P. C. Hohenberg, and P. Pincus, Phys. Rev. 147, 140 (1966); A. L. Fetter, Phys. Rev. 163, 390 (1967). [23] A. L. Fetter, Phys. Rev. 162, 143 (1967). [24] V. K. Tkachenko, Zh. Eksp. Teor. Fiz. 49, 1875 (1965) [Sov. Phys. JETP 22, 1282 (1966)]. [25] G. Baym, Phys. Rev. B 51, (1995); Phys. Rev. Lett. 91, (2003). [26] H. Fukuyama, Sol. State. Comm. 17, 1323 (1975). [27] M. Weber et al., Phys. Rev. B 48, (1993). [28] We estimate u rms from the published data in Ref. [1]. Also, using the correspondence between the LDOS modulations and the vortex core states [2], we obtain a similar number from the measurements of Ref. [3]. [29] J. Bardeen and M. J. Stephen, Phys. Rev. 140, A1197 (1965); P. Nozières and W. F. Vinen, Philos. Mag. 14, 667 (1966). [30] J. M. Harris, Y. F. Yan, O. K. C. Tsui, Y. Matsuda, and N. P. Ong, Phys. Rev. Lett. 73, 1711 (1994). [31] P. J. Turner et al, Phys. Rev. Lett. 90, (2003). [32] B. Keimer, W. Y. Shih, R. W. Erwin, J. W. Lynn, F. Dogan, and I. A. Aksay, Phys. Rev. Lett. 73, 3459 (1994). [33] H.-T. S. Lihn, S. Wu, H. D. Drew, S. Kaplan, Qi Li, and D. B. Fenner, Phys. Rev. Lett. 76, 3810 (1996) and references therein. 11