Plasmonic Drag in a Flowing Fermi Liquid

Haoyang Gao, Zhiyu Dong, Leonid Levitov Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA

Collective modes in two-dimensional electron fluids show an interesting response to a background carrier flow. Surface plasmons propagating on top of a flowing Fermi liquid acquire a non-reciprocal character manifest in a ±k asymmetry of mode dispersion. The nonreciprocity arises due to Fermi surface polarization by the flow. The flow-induced interactions between quasiparticles make col- lective modes of the system uniquely sensitive to subtle “motional” Fermi-liquid effects. The flow- induced Doppler-type frequency shift of plasmon resonances, arising due to electron interactions, can strongly deviate from the classical value. This opens a possibility to directly probe motional Fermi-liquid effects in plasmonic near-field imaging experiments.

Plasmonic drag, also known as the plasmonic motional effects. The quantity α describes the Doppler effect, is a motional effect that describes curvature of the band dispersion, such that α = a change in the dispersion of collective charge os- 0 for the linearly-dispersing carriers and α = 1 cillations induced by an electric current driven for parabolic dispersion. This result is valid at through the system. As a simplest case of mo- relatively weak interactions; a more complicated tional coupling between two different collective behavior is found for stronger interactions using flows, the collective oscillations and the DC cur- a relativistic Landau Fermi-liquid framework9 rent, plasmonic drag is of interest for the quest As a quick sanity check, taking F1 = G1 = 0 for new effects due to the electron-electron inter- yields the classical Doppler shift δω = ku when actions, and new transport phenomena due to band dispersion is parabolic [α = 1]. In con- such effects. Graphene plasmonics1–4, in partic- trast, for linear dispersion [α = 0] Eq.(1) pre- 10 1 ular the near-field imaging techniques developed dicts a nonclassical Doppler shift δω = 4 ku. recently5,6, provide a platform in which the plas- This behavior of plasmonic drag displays an monic drag effects can be realized and explored. interesting analogy with the seminal results Here we investigate plasmonic drag in a flow- on motional effects in superconducting Fermi ing Fermi liquid. The Fermi-liquid interactions liquids11,12. The current-current correlation are known to be unimportant for plasmons in function, which determines the response of su- systems with parabolic electron band dispersion, percurrent to vector potential, was found to be where the collective center-of-mass motion of strongly renormalized by the Fermi-liquid inter- charges can be separated from their relative mo- actions. However, these renormalization effects, tion due to the Galilean symmetry7,8. However, while nominally big, feature a cancellation for as we will see, a very different situation occurs systems with parabolic bands. for electron systems with a nonparabolic band To emphasize the sensitivity of the Doppler dispersion such as that of graphene. shift to fundamental symmetries of Bloch elec- In this case, as we will see, the Fermi-liquid in- trons, such as Galilean symmetry for parabolic teractions do renormalize the Doppler shift. Our bands and Lorentz symmetry for Dirac bands, analysis, which fully accounts for the interaction it is instructive to make comparison with light effects, predicts the change in the plasmonic fre- drag in optics. Known as Fizeau drag13, it arises quency in the presence of the flow: due to the speed of light dependence on the ve-

arXiv:1912.13409v1 [cond-mat.mes-hall] 31 Dec 2019 locity of a transparent, moving medium. For a 1 3 G + α δω = ku + 1 + O(u2). (1) slowly moving medium, Fizeau drag is a ±k-odd 4 4 1 + F1 effect first-order in the medium velocity u,

Here u is the drift velocity, F1 is the m = 1  1  δω = uk 1 − √ , (2) harmonic of the Landau interaction and G1 is n its radial derivative defined below. The quantity G1, as we will see, is uniquely sensitive to the where n is the medium refraction index. The re- 2 duction of the frequency shift compared to the ions or charge on the gates. The last two classical Doppler shift δω = uk is a distinct sig- terms in Eq.(3) represent the potential energy nature originating from the symmetry of space- of a particle due to a change in the distribu- time and special relativity. tion of other particles, those far away and those For the plasmonic Doppler effect in graphene, nearby. Distant particles contribute the long- our analysis predicts a similar suppression as range Coulomb potential which arises due to a compared to the classical Doppler effect, a dis- change in the net density of charge at a remote P 0 0 tinct behavior arising due to relativistic carrier point. The term p0 f(p, p )n(p ) is the angle- dispersion in graphene. At the same time, the dependent spatially-local interaction of the Lan- effect is considerably stronger than the light dau Fermi-liquid theory. drag, on the order δω/ω ∼ u/vp, where vp We note parenthetically that the apparent sin- is plasmon velocity. Our analysis also demon- gularity at x0 = x is an artifact of our decompo- strates that the Doppler shift is further renor- sition of the potential into a sum of the remote malized, and enhanced, by interactions in the Coulomb part and the local Fermi-liquid part, flowing Fermi liquid. where ‘local’ and ‘remote’ is defined relative to The dependence of the Doppler shift on the the Fermi wavelength. While this decomposition band curvature α and the electron interac- is somewhat ambiguous, it will be seen that the tions (F1 and G1) can provide a way to tune expression above is mathematically sound and the Doppler shift. If band curvature is large well behaved: It is free from divergences arising and positive, the Doppler shift is enhanced, at x0 ≈ x and provides a correct description in whereas when curvature is large and negative, the long-wavelength limit of interest. the Doppler shift sign is reversed. The interac- We will write the particle distribution as a tions F1 and G1 renormalize the Doppler shift and push it away from the free-particle value. sum of the parts describing a steady-state equi- Measuring plasmonic Doppler effect can there- librium in the presence of a flow and a pertur- fore be used to directly probe motional Fermi- bation describing collective charge oscillations: liquid effects. Comparison of the effects for dif- ferent electron band dispersion can shed light on n(p, x, t) = nu(p) + δn(p, x, t), (5) subtle aspects of Bloch electron dynamics 1 In our analysis, we will focus on a two- nu(p) = P 0 0 . β(0(p)+ 0 f(p,p )nu(p )−up−µ) dimensional Fermi liquid in the collisionless e p + 1 regime ω  γee, where γee is the carrier colli- sion rate. In this case, while the effects of colli- Here the subscript u indicates that the momen- sions are negligible, the effects of ee interactions tum distribution is altered by the flow. are not negligible because carriers are subject Since nu(p) appears under the Fermi function to the short-range Landau interactions in com- that defines nu(p), it may seem that the de- bination with long-range Coulomb interactions. pendence of current on the flow velocity u must This system is described by the single-particle take a nonclassical form. Yet, this dependence Hamiltonian takes a completely conventional form. This can X 0 0 H = 0(p) + eφ(x) + f(p, p )n(p , x), (3) be seen by starting with the expression for cur- p0 rent that accounts for a change in velocity due to Fermi-liquid interactions with a u-dependent where 0(p) is particle dispersion and φ(x) is the particle distribution: electrostatic potential Z n(x0) − n¯   φ(x) = d2x0e , (4) |x − x0| X X 0 0 j = e∇p 0(p) + f(p, p )nu(p , x) nu(p) R d2p P p p0 Here n(x) = 2 n(p, x) 0 ... is the den- (2π) p X sity of distant charges; the quantity −n¯ de- = eunu(p) = en¯u. (6) notes compensating background charge due to p 3

Here we integrated by parts using the identity where vF = d(p)/dp at p = pF , and we intro- duced angular harmonics of the Landau interac- ∇p ln(1 − nu(p)) (7) tion defined in the standard way:   X −im(θ 0 −θp) 0 X Fm = e p f(p, p )δ((p) − µ). (12) = β v (p) − u + ∇ f(p, p0)n (p0) n (p),  0 p u  u p0 p0 From Eq.(11) we find the relation where v (p) = ∇  (p). m u 0 p 0 ∆p = ∗ (13) The result in Eq.(6) identifies the quantity u, 1 + F1 introduced above as a convenient parameteriza- tion of the flowing carrier distribution, with the where we defined m∗ = pF /vF the quasiparticle drift velocity defined in a conventional way as effective mass. The dynamics of our system is described by j = envd. Below we study collective charge os- cillations in the presence of the flow and deter- classical equations of motion mine the plasmonic Doppler shift by carrying ∂tn + {H, n} = 0 (14) out perturbation theory in u. The relation in Eq.(6) can then be used to express the Doppler where {A, B} = ∇pA∇xB − ∇xA∇pB are clas- shift through the actual electric current. sical Poisson brackets. We linearize the Hamil- A nonclassical relation that does arise is the tonian in the carrier distribution perturbed one for the Fermi surface displacement induced away from equilibrium as given in Eq.(5), ar- by the flow. Working at small u and assuming riving at X a change in particle distribution due to current H = ε(p) + eδφ(x) + f˜(p, p0)δn(p0, x, t) that happens only near the Fermi level, we can p0 represent the distribution as a displaced Fermi (15) surface X ε(p) ≡ (p) + f(p, p0)∆n(p0), (16) p(θ) = pF + ∆p cos(θ). (8) p0 where δφ is the potential of a distant charge The amplitude of the displacement ∆p can be perturbation, δφ(x) = R d2x0 e δn(x0). The found from the relation defining the Fermi sur- |x−x0| face, quantities ε(p) and (p) are the quasiparticle energy in the presence and absence of u, re- X ˜ 0  (p) + f(p, p0)n (p0) − up = µ, (9) spectively; f(p, p ) is the Landau function for 0 u ˜ p0 a shifted Fermi surface. The relation between f and f will be discussed below. through rewriting it in terms of the change of To proceed with the analysis, we defineε ˜(p) = the distribution due to the flow ε(p) − p · u the quasiparticle energy in the pres- ence of flow viewed in the comoving frame, then ∆n(p) = nu(p) − n0(p). (10) the steady-state distribution describing current flow can be written as As always in the Fermi-liquid theory, it will be n (p) = θ(µ − ε˜(p)) (17) convenient to absorb the contribution of a non- u moving Fermi sea in the quasiparticle energy, Using this notation and the Hamiltonian in P 0 0 (p) = 0(p) + p0 f(p, p )n0(p ). Combining Eq.(15), we linearize equations of motion, with Eq.(8), we can describe the displaced Fermi Eq.(14), to obtain surface as   X ˜ 0 0 X 0 0 ∂tδn = −eE + f(p, p )∇xδn  ∇pnu(p) 0 = vF ∆p cos θ + f(p, p )∆n(p ) − upF cos θ p0 p0 − v ∇ δn (18) = vF ∆p(1 + F1) − upF cos θ, (11) p x 4

Since the Fermi surface, after being shifted, is still approximately circular (at lowest order in u/vF ), we find it convenient to reparameterize all quantities with q ≡ (q, ϕ), denoting the mo- mentum and angle measured from center of the shifted Fermi sea (See Fig.1). Using the new coordinate system, all the quantities can be written explicitly. The shifted Fermi sea at zero temperature is simply

nu(p) = θ (q(p) − pF ) . (20)

The perturbed distribution δn(p) can be ex- pressed through Fermi surface normal displace- FIG. 1. The dashed red and the black dotted lines ment vs. polar angle ϕ on the shifted Fermi mark the Fermi surface in the presence and absence surface: of the flow, respectively. Shown are the coordi- nates (p, θ) for the Fermi surface at rest, and the coordinates (q, ϕ) for the Fermi surface describing a h2 flowing Fermi liquid. Also shown is the vector v˜ = δn(p) = δ(q − pF )δn(ϕ) (21) pF ∇pε˜(p) normal to the shifted Fermi surface (dashed red line), which is the contour ofε ˜(p) = ε(p) − p · u (see text). This relation allows us to convert any integral over p involving δn(p) into an integral over ϕ: where δn and δn0 is a shorthand for δn(p, x, t) 0 2 and δn(p , x, t), respectively; E = −∇xδφ and Z dp Z p dq δn(p)U(p) = F δn(p)U(p) we defined the velocity in the lab frame vp ≡ h2 h2 ∇pε(p). Z To describe collective modes, we consider per- = dϕδn(ϕ)U(ϕ, q = pF ) (22) turbations of a plane-wave form, δn(p)eikx−iωt. Writing the field of distant charges as eE = P 0 2πe2 −ikV (k) p0 δn(p ), V (k) = k , and substi- where U(p) can be an arbitrary function. tuting in Eq.(18), gives an integral equation for Another useful property of the coordinates the collective mode: (q, ϕ) is that the ϕ directly labels the direction of X v˜, because being the gradient ofε ˜(p), the v˜ has (k · v − ω) δn(p) + k · ∇ n (p) f˜(p, p0)δn(p0) p p u to be perpendicular to the shifted circular Fermi 0 p surface, which is the contour ofε ˜(p). This fact X 0 = −k · ∇pnu(p)V (k) δn(p ). (19) will be useful later when we evaluate the velocity x p0 component v .

1 B. Wunsch, T. Stauber, F. Sols, F. Guinea, Dy- graphene. Phys. Rev. B 75, 205418 (2007). namical polarization of graphene at finite doping. 3 M. Polini, A. H. MacDonald, and G. Vignale, New J. Phys. 8, 318318 (2006). arXiv:0901.4528; S. H. Abedinpour, et al. Phys. 2 E. H. Hwang, S. Das Sarma, Dielectric func- Rev. B 84, 045429 (2011). tion, screening, and plasmons in two-dimensional 4 F. H. L. Koppens, D. E. Chang, F. J. Garca 5

de Abajo, Graphene Plasmonics: A Platform for 10 D. S. Borgnia, T. V. Phan, L. S. Levitov, Quasi- Strong LightMatter Interactions Nano Lett. 11, Relativistic Doppler Effect and Non-Reciprocal 8, 3370-3377 (2011) Plasmons in Graphene, arXiv:1512.09044 5 J. Chen et al., Optical nano-imaging of gate- 11 A. I. Larkin and A. B. Migdal, Theory of Super- tunable graphene plasmons, Nature 487, 77-81 fluid Fermi Liquid. Application to the Nucleus, (2012). Zh. Eksperim. i Teor. Fiz. 44, 1703 (1963) [En- 6 Z. Fei et al., Gate-tuning of graphene plasmons glish transl: Soviet Phys. JETP 17, 1146 (1963)]. revealed by infrared nano-imaging, Nature 487, 12 A. J. Leggett, Theory of a Superfluid Fermi Liq- 82-85 (2012). uid. I. General Formalism and Static Properties, 7 G. F. Giuliani, G. Vignale, Quantum Theory Phys. Rev. 140, A1869 (1965). of the Electron Liquid (Cambridge Univ. Press, 13 L. D. Landau, E. M. Lifshitz, Electrodynamics of 2005). Continuous Media (Pergamon Press, 1960) 8 T. N. Theis, Surf. Sci. 98, 515-532 (1980). 14 L. S. Levitov, A. V. Shtyk, and M. V. Feigel- 9 G. Baym, S. A. Chin, Landau theory of relativis- man, Electron-electron interactions and plasmon tic Fermi-liquids, Nuclear Physics A262 (1976) dispersion in graphene, Phys. Rev. B 88, 235403 527-538 (2013)