The anomalous attenuation of plasmons in strange metals and holography

The plasmon is a ubiquitous collective mode in charged liquids. Due to the long range Coulomb interaction, the massless {zero} sound mode of the neutral system acquires a finite plasmon frequency in the long wavelength limit. In the zero temperature state of conventional metals -- the Fermi liquid -- the plasmon lives infinitely long at long wavelength when the system is (effectively) translationally invariant. In contrast, we will show that in strongly entangled strange metals the protection of zero sound fails at finite frequency and plasmons are always short-lived regardless of their wavelength. Computing the explicit plasmon response in holographic strange metals as an example, we show that decay into the quantum critical continuum replaces Landau damping and this happens for any wavelength.


I. INTRODUCTION
The plasmon is a ubiquitous propagating mode in electromagnetically charged systems.
It is governed by a simple universal principle. Any system, which in the absence of electromagnetic interactions, carries a propagating gapless density mode with dispersion ω = v s |p|, will, upon switching on the static Coulomb interaction, see this dispersion change to ω = ω 2 p + v s p 2 with ω 2 p = ρ 2 /( + p) the plasma frequency; here ρ is the (charge) density, -the energy density, and p -the pressure. In non-relativistic systems this becomes ω 2 p = ne 2 /m with e and m the charge and the mass of the charged particles at a density n. The plasmon's ubiquity ranges from crystals formed from charged particles in a homogeneous background, where the longitudinal phonon becomes the plasmon [1], to superconductors.
Here the plasmon originates in the Nambu-Goldstone boson through the Anderson-Higgs mechanism and takes the role of the longitudinal photon polarization. 1 Here we focus on the plasmons found in the zero temperature metallic state formed by electrons in solids. Normal metals are described by the Fermi-liquid (FL) theory coupled to electromagnetism and the explanation of plasmons is part of its success. The zero temperature uncharged Fermi-liquid with repulsive s-wave interactions (the Landau parameter F s 0 > 0) is characterized by a propagating longitudinal density mode with dispersion ω = v s |p|: the zero sound. Similar to the usual "first sound" of classical fluids, it is hydrodynamically protected in the sense that at low energies and long wavelengths momentumand particle number conservation inhibit its decay. Upon including the long range Coulomb interaction, this zero sound is promoted to the plasmon which requires a finite frequency ω p to excite.
At this finite frequency, the universality of hydrodynamical-like protection no longer applies and one would expect the plasmon to decay. Nevertheless, the Fermi-liquid plasmon is stable. This is well understood thanks to the fact that the Fermi-liquid renormalizes into a system of non-interacting quasi-particles. Using this microscopic description one can compute that, besides the propagating plasmon, the Fermi-liquid can also react to electromagnetic probes by "incoherent" particle-hole excitations. The phase space for these excitations is limited by energy and momentum conservation. This phase space -the Lind- 1 In relativistic systems the longitidunal and transversal polarizations combine in order to form massive vector bosons with mass equal to the plasmon frequency. In non-relativistic superconductors characterized by a material velocity v F c the mass for the transversal vector polarizations, which sets the London penetration depth, is a factor c/v F smaller then ω p .
hard continuum -shrinks to a small region near zero frequency in the long wavelength limit (see Section II). Therefore at low momenta the plasmon is completely isolated from these excitations; it has no channels to decay and inherits thereby the hydrodynamical protection of the zero sound. At larger momenta the plasmon dispersion will cross the upper bound of the continuum, at which point it will rapidly decay in uncorrelated particle-hole excitations: the classic Landau damping (see Fig. 1).
In any real metal translational symmetry is broken by the atomic lattice and momentum is not strictly conserved. Hydrodynamical protection in the strictest sense therefore does not apply, neither for zero sound nor its inheritance by the plasmon. However, in weakly interacting Fermi liquids as realized in simple metals, the resulting Umklapp scattering becomes noticeable in the zero momentum charge response only upon reaching the energy of interband transitions. When ω p is not coincident with this energy, the lattice continues to be irrelevant and long wavelength plasmons continue to be long lived. This physics is at the basis of a flourishing engineering field that makes use of the beneficial properties of the plasmon as "material photon": plasmonics.
However, these well established wisdoms are violated in a spectacular fashion in the non-Fermi liquid strange metals realized in e.g. the cuprate high T c superconductors. In the p → 0 limit the plasmons can be studied via optical conductivity and are found to be almost overdamped: their inverse life time is of order of the plasmon frequency [2,3]. This is surely not due to an energy independent elastic scattering since the momentum relaxation rate observed in the DC conductivity at low temperatures is smaller by orders of magnitude than the observed plasmon attenuation. The charge density response at finite p has been measured as well a while ago with a limited energy resolution using high energy transmission energy electron loss spectroscopy. These data reveal a dispersing plasmon becoming completely overdamped at a momentum of approximately 1/3 of the Brillouin zone [4][5][6].
Very recently the data obtained by low energy reflection electron loss spectroscopy were presented [7]. Although these data suggest that the plasmon is disappearing even faster for increasing momentum, these also appear to reveal directly the nature of the incoherent excitations. Instead of the Lindhard continuum it was found that the spectrum is independent of momentum and energy up to a high energy cut-off at 1 eV.
So much is clear that the atomic lattice effects and the Umklapp scattering are less "hidden" in these strongly interacting electron systems as compared to the simple metals. This is already seen in exact diagonalization studies on small Hubbard model clusters [8] and in the recent quantum Monte Carlo results [9]. For strange metal models described by a holographic framework, which we will be considering here, it is also established that lattice effects significantly modify the finite frequency response in optical conductivity [10][11][12]. Undoubtedly, these effects have to be dealt with in order to address what is going on in experiment.
This we will not pursue. We will address instead a matter of principle of such a generality that it may well play a crucial role in any attempt to explain experiment. Consider a translationally invariant system, eliminating the complicating lattice effects. What is the fate of the plasmon in a metal that does not renormalize to the free fermion fixed point controlling the Fermi liquid?
We use the holographic strange metals as a model system which we claim is in the present context fully representative for such non-Fermi-liquids. A defining characteristic of these systems is that their charge response is characterized by two "sectors": the zero sound as in the neutral Fermi liquid, and in addition a "quantum critical" (QC) continuum [13][14][15].
The additional significance of this quantum critical sector is that it cannot be understood in terms of particle-like excitations and it therefore has a unique imprint in transport and collective responses. This sector should not be confused, however, with the fluctuations associated with a quantum critical point. The holographic strange metal shows scaling behaviors characterized by hyperscaling violation (θ) and dynamical critical (z) exponents, which are set by the deep infrared "quantum critical" sector of the dual model but that are unrelated to the familiar "bosonic" quantum criticality.
This very presence of the QC sector in holographic models, seemingly unrelated to anything that is familiar in conventional Condensed Matter, has been a puzzle. Most of the holographic strange metal studies are focused on the zero momentum collective observables (rare exceptions include [16][17][18]). Collective responses in the Fermi liquid at zero momentum, on the other hand, lack any kind of a "second sector" contribution. However, this situation changes once one turns to finite momentum. Then an algebraic second sector arises in Fermi liquid as well: the Lindhard continuum. Clearly, it is very different from the scaling infrared continuum of the strange metal, as we will discuss in more detail below. But effectively, the holographic strange metals may be indeed viewed as "strongly interacting" (in the critical theory sense) generalizations of the Fermi liquid, with two-sector collective  In conventional Fermi Liquids the width of the plasmon (blue lines) at large momenta is due to decay into particle-hole (Lindhard) continuum (left). In strange metals there is a quantum critical continuum at all frequencies and momenta, and a plasmon in a strange metal can therefore decay at any momentum (right).
responses where the QC sector supplants the Lindhard continuum.
From this perspective, the vanishing of the Lindhard continuum at zero momentum is a singular feature of the free fermion fixed point. In a non-Fermi liquid an incoherent continuum should be present at any finite frequency and momentum. The ramification is that, as we show below, even in the perfect Galilean continuum the long wavelength limit plasmon will decay in finite time (see Fig. 1). The infinitely long lifetime of the Fermi liquid plasmon is not rooted in general principle, but instead in the singularly special nature of its free particle IR fixed point.
In the remainder of this paper we will provide the evidence for this case. In Sec.II we will remind the reader of the results of electromagnetic linear response theory in metals, tying together the electrical conductivity with the charge density response and highlighting the universality of the plasmon. We will also specify how this can be married with the nominally neutral matter computed by holography. In Sec.III we introduce the holographic models for strange metals and show how the quantum critical sector is encoded in the geometry. We also review how their zero sound response is encoded in quasi-normal mode fluctuations of this dynamical space-time geometry. In Sec.IV we show how the long-range Coulomb interaction can be encoded in the holographic dictionary as a so-called double trace deformation, a point also made in the recent article [19]. This allows us to compute the corresponding charged density-density response function and study the corresponding modified quasi-normal mode spectrum. Our results are consistent with earlier studies of the plasmon in holographic strange metals where the charge response is computed in a different manner [20][21][22]. This central section also contains the main results: as a consequence of the omnipresent QC sector into which the holographic strange metal plasmon can decay, the width of the plasmon is significantly enhanced and finite even at zero momentum. We devote Sec.V to the study of the dependence of the features of plasmon mode on the characteristics of the quantum critical sector by tuning its scaling dimensions. The surprise is that the plasmon width at zero momentum is quite insensitive to these vast changes in the deep IR.
Appendix B is devoted to how to go beyond the static Coulomb interaction and include dynamical electromagnetism. Appendix C generalizes our 3D results to an effective 2D Coulomb interaction V p = e −λ|p| /|p|. We are envisioning an electrically charged liquid. Let us first collect central results from the standard theory describing the linear response of such a system to external electromagnetic fields [23,24]. The central quantity enumerating the response of the charged medium is the Polarization propagator Π(ω, p). To make contact with condensed matter physics we consider non-relativistic systems characterized by a material velocity scale that is much smaller than the velocity of light. In the transversal channel Π(ω, p) is just the photon self-energy associated with the interaction with matter. Our focus is here entirely on the longitudinal channel, which encodes the response to longitudinal electric fields. It is captured by the dielectric function, which relates the displacement to the electrical field D(ω, p) = ε(ω, p) E (assuming an isotropic medium). The dielectric function ε(ω, p), the frequency and momentum dependent conductivity σ(ω, p) and the charge density response function χ(ω, p) are determined in terms of the polarization propagator and Coulomb interaction V p = e 2 /p 2 as One notices that the optical conductivity is most closely related to Π(ω, p); this is a consequence of the continuity equation associated with the conservation of electrical charge that links the density and longitudinal currents through ∂ t n + ∇ · J = 0. The practical difficulty with the optical conductivity is that it can only be measured at p → 0 because of the smallness of the material velocity compared to the velocity of light. There is much to learn at finite momenta, but this is only accessible by mobilizing electrons with their small wavelengths through electron loss spectroscopy; these measure the loss function Im 1 ε(ω,p) or the charge susceptibility χ(ω, p) directly.
In a continuum charged Fermi liquid with a repulsive F 0 s Landau parameter the polarization propagator can be parametrized at zero temperature, and not too large momenta as This low momenta response of the polarization operator Π(ω, p) has the same structure as the (particle number) density response χ 0 (ω, p) = n(ω, p)n(−ω, −p)) of a neutral Fermi liquid like 3 He. It is characterized by a zero sound mode with velocity v s and a hydrodynamical-like damping rate Γ(p) = Dp 2 where D is the diffusivity. In the nonrelativistic Fermi liquid its pole strength equals A = n /m ( n and m being the carrier density and electron mass, respectively). In addition, it contains the incoherent "second sector" denoted with Ξ(ω, p). This is the Lindhard continuum. In the 3+1 D Fermi gas it has the well known branch cut form, This continuum is bounded from above by energy and momentum conservation for a quasiparticle-hole pair near the Fermi surface: ω max (p) = p F |p| m + p 2 2m . The bounding value ω max goes to zero in the long wavelength limit. One observes that continuum obeys a z = 1 dynamical critical scaling while at finite momentum its imaginary part increases linearly at low frequency in d ≥ 2 space dimensions (see e.g. [26]).
It follows from Eq. (1) that for small p the conductivity behaves as or in the strict p → 0 limit The second term is well known in holography as the "incoherent conductivity" [27][28][29][30][31], arising from the "quantum critical" sector of the model. In a Fermi-liquid this contribution vanishes at p = 0, but would be seen at finite momentum, being associated in this case with the Lindhard continuum. However, because Ξ(ω, p) vanishes at p = 0, all what remains in the Fermi-liquid is the infinitely long lived sound showing up in the conductivity as the diamagnetic delta function at zero frequency ∼ iA/ω. It is a peculiarity of the free fermion fixed point that the "second sector" contribution disappears in the long wavelength limit.
To finalize our brief review of electromagnetic response in charged matter, it also immediately follows from Eq. (1) that the charge susceptibility χ(ω, p) (or the loss function) has a pole when Π(ω, p) = 1/V p . Substituting in Eq. (2) one finds a propagating gapped excitation with dispersion ω = (v s p) 2 + ω 2 p , damped by Γ. The zero sound gets "promoted" to the plasmon when one observes the density response. Eventually the plasmon becomes overdamped when it "dives" in the Lindhard continuum at large momenta: the conventional Landau damping.
We wish to address the electromagnetic response in strange metals described by holography in the next section, but to do so we face a problem of principle: the standard holographic strange metals are electrically neutral fluids. The intrinsic strong interactions, which form a holographic strange metal, are in principle not of electromagnetic origin. There is however a loophole. In the treatment presented above the Coulomb interactions were seperated into two parts. The long range Coulomb force is accounted for perturbatively via the factors V p . On the other hand the short ranged part of the screened Coulomb interactions gives rise to the polarization propagator Π(ω, p) of the charged electron system. A typical representative of this kind is the Hubbard U -potential that may well be of crucial importance for the emergence of the real strange metals. As we noted below Eq. (2), the polarization propagator in charged electron Fermi-liquids, however, shares the same zero-sound structure as the density response in neutral Fermi-liquids such as 3 He. Computing the density response in a neutral Fermi-system can therefore be used as an analogy for computing the short-range electromagnetic response in a charged Fermi-system. Similarly, we can thus consider the strong intrinsic interactions of the neutral holographic model to correspond to short range Coulomb forces and regard the holographic strange metals as a model system that captures generic features of the charged strange metals. Subsequently the long range part of the Coulomb interaction can be regarded weakly coupled and treated in exactly the same fashion as above.
To make this more precise, let us specify some formalism. The model system should be characterized by a conserved current J µ with n(p) ≡ J 0 (p), p = (ω, p). This current is externally coupled to dynamical electromagnetism with a strength e -in relativistic notation, This equation of motion together with the dynamics of the strange metal is conveniently encoded in the action for the electromagnetic potential A µ Fixing a Lorentz gauge ∂ µ A µ = 0, the electromagnetic potential can be integrated out. In the absence of external photon sources the effective action becomes For non-relativistic condensed matter systems characterized by velocities v F c and η 00 η ii the electromagnetic interaction reduces to the static 3D Coulomb interaction, If the Coulomb interaction is sufficiently weak, it can be fully accounted for by time dependent mean field [25,32]: the Random Phase approximation (RPA) amounting to a simple "bubble resummation", see Fig. 2. The result is the dynamical charge susceptibility χ(ω, p), This has the same form as the electromagnetic charge density response with the difference that the density-density propagator of the neutral model system χ 0 (ω, p) has replaced the polarization propagator of the charged system. If the neutral system has a zero sound response, χ 0 has the same structure as Eq. (2). As we will see, this is also in the case of the holographic strange metals with the difference that in that case Ξ(ω, p) never vanishes. Let us therefore insert Eq. (2) in the expression for the charge susceptibility (10). This yields Dyson resummation for the "dressed" charge density two-point correlation function χ(ω, p) (10). The diagrams included in the RPA approximation are shown as bubbles, where χ 0 (ω, p) is the density-density correlator of the neutral Fermi liquid (2), and V p is the Coulomb interaction (9).
One recognizes the plasmon excitation, but its mass, attenuation and pole strength receive contributions from the incoherent "second sector", With this parametrization, we can now isolate the signatures of the specific structure of the Ξ(ω, p) in the observable density response and compare strange metals to well-known Fermi-liquids. Inspecting Eq. (12) we see that: 1. The Coulomb interaction promotes the p → 0 massless zero sound mode to the plasmon frequency ω 2 p in the longitudinal charge response. In the Fermi liquid this plasmon frequency is set by the Drude weight ω 2 p = AV p p 2 alone, as Ξ vanishes at p = 0. 2 In the two-sector strange metals this is no longer the case. There it also receives a contribution from the finite weight Re Ξ(ω, p = 0) of the underlying quantum critical continuum.
2. The inverse lifetime of the plasmonΓ is determined both by the "hydrodynamical" zero sound attenuation factor Γ p ∼ p 2 and a factor associated with the spectral function ImΞ(ω, p). In the Fermi-liquid, this is responsible for the Landau damping of the plasmon that only occurs at large momenta. In the two-sector strange metals this contribution may cause damping at all momenta. This includes the long wavelength limit all the way to p = 0.
We will now verify these predictions in a general class of holographic strange metals. 2 In real metals, interband transitions can shift the plasmon frequency.

METALS.
We first show how to compute the density response in holographic strange metals. We consider systems which at low temperature have a strongly coupled quantum critical sector determined in terms of a hyperscaling violation (θ) and dynamical critical (z) exponents.
These are readily modeled in terms of anti-de-Sitter(AdS) Einstein-Maxwell-Dilaton gravity -see [13] for a review. The U (1) gauge field theory in the bulk encodes for a global U (1) symmetry in the quantum critical matter system of the boundary, associated with the particle number density. Conventional holography does not take into account the (external) electromagnetic force in the boundary but we learned already in the previous Section how to deal with that (see also next Section). Similar to the neutral Fermi liquid, the finite density translationally invariant holographic strange metals are invariably characterized by a zero sound mode [33][34][35] [36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51][52][53][54][55]. This fact was originally part of the motivation to identify these holographic models as strange metals. Here we are interested in the density-density correlation function. As mentioned above, the number density operator n(p) is dual to temporal component A 0 (p) of the U (1) gauge field in the bulk, i.e. its boundary value A 0 (p) acts as the source for n(p). In terms of the partition function of the "boundary" strongly coupled theory it reads Applying holography, this can be equated with the (exponential of the) on-shell gravitational action The (connected) two point function is thus obtained by taking the second variation of the on-shell action (14) with respect to the sources A 0 : Doing so one arrives at the simple result that the two-point function can be read off from the near-boundary behavior of the solution to the equations of motion for A 0 (p; u). For a U (1) gauge field, this is of the form with α > 0 being the exponent of the subleading mode near the boundary, dictated by the equations of motion. As mentioned above, in standard holography the perturbative source just sets the Dirichlet boundary condition for the bulk field A 0 The on-shell action S AdS reduces to the boundary term therefore the expectation value δS/δA induced by the source is associated with the subleading coefficient b(p). Below we will have the case with α = 1 in Eq. (16) meaning that the expectation value is given by a first derivative on the boundary The holographically computed two point function is therefore (c.f. (15)) For a unique solution one does need additional boundary conditions in the interior. For a black hole spacetime as we shall consider here, choosing infalling boundary conditions at the black hole horizon results in the retarded thermal Green's function in Eq. (20) [56].
In what follows we compute the density-density correlator in the Einstein-Maxwell-Dilaton backgrounds [43,57] encoding for translationally invariant holographic strange metals with a quantum critical sector that can be tuned in a flexible manner. These have the gravitational action Without loss of generality we choose the bulk couplings for gravity 2κ 2 = 1 and for the gauge field q = 1. Moreover, we choose the following potentials We then solve the equations of motion numerically for D = 4 dimensions subject to the assumption of translational invariance and without a probe source for the density operator The resulting spacetime geometry is dual to a 2 + 1 dimensional strange metal state and the presence of a quantum critical sector is exhibited through an emergent scaling IR geometry at zero temperature of the form [58]: In what follows we will use particular values of z and θ that ensure a vanishing ground state entropy [58]. We often consider the case of large z in order to make connection with previous results in the literature (numerically the largest value we use is z = 6401). We will simply denote this value as z → ∞.
On top of such backgrounds we then study the linearized equations of motion of momentum-dependent perturbations induced by the source A 0 (p). This is a system of eleven coupled ordinary differential equations which we solve again numerically; details are also discussed in Appendix A. The result for the two-point function of the density operator is then extracted from the solution following Eq. (20). Its spectral function Imχ 0 (ω, p) is plotted in Fig. 3 for the simplest such model: locally quantum critical strange metal with z = ∞, θ = 0. One clearly sees the zero sound peak with characteristic linear dispersion and zero width in long-wavelength limit.
There is a useful technique to analyze these coherent peaks more precisely. We study the Sturm-Liouville problem, which we further discuss in Appendix A). The frequencies that satisfy this condition specify the poles of the two-point function (20) and correspond to the quasi-normal modes of the dual black hole. For a given momentum, these poles will lie in the lower half complex ω plane. The pole at (complex) ω s (p) with lowest frequency is the one corresponding to the zero sound mode with the characteristics In the bottom right plot of Fig. 3, we show the zero sound pole and its motion as a function of momentum |p|. We also show how this pole is reflected in the spectral function Imχ 0 (ω, p) which is plotted in the top right plot of the same figure. We clearly see that the coherent response is well described by a single pole approximation. The sum of the two fits (continous yellow) describes reasonably well the full density response. Here we take |p|/µ = T /µ = 0.02.
At the same time, the coherent zero sound peak is not the only contribution that can be seen in the density-density correlator χ 0 (ω, p). Fig. 4 shows that there is an additional incoherent contribution Ξ(ω, p) to the density-density correlation function, as we anticipated in Eq. (2). The "quantum critical sector", which gives this incoherent contribution, is omnipresent in holographic strange metals and reflects the "deep infrared scaling geometry" near the black hole horizon in the bulk. The effects of this second parallel sector have been studied at length in the holographic literature, especially in connection to the temperature dependence of the DC electrical conductivity, at zero frequency and momentum where it becomes possible to compute it analytically. Here we observe its manifestation in the density response, which is related to the conductivity by the continuity equation, see Eq. (1).
The observation that the optical response in holographic strange metals is governed by two seemingly independent systems -one of them responsible for the Drude peak, directly associated with zero sound in our case; the other one due to the quantum critical sectorwas already made early on [27][28][29][30][31]. At p = 0 it is easy to deduce the temperature scaling of this QC-contribution to the DC-conductivity at zero momentum using general scaling arguments [59,60]: Here d equals the spatial dimensions and z, θ are the dynamical critical and hyperscaling violation exponent determined by the deep interior geometry [61,62]. Φ encodes an anomalous scaling of the charge density and depends both on the geometry and on the gauge field part of the gravitational action [59]. In our EMD model we have Φ = z always (for the absolute value, see Sec. 8.6 of [57]); for d = 2 this gives The same scaling arguments then lead to the frequency scaling at zero temperature and small frequencies: which together with Eq. (1) suggests that At finite momentum, however, scaling only fixes Ξ up to an arbitrary function, In specific holographic models its asymptotic behavior near ω = 0 can be computed, e.g. for the class of 3+1 dimensional EMD models where θ → −∞, z → ∞ with η = − θ z fixed, one has for infinitesimal ω [17] σ QC (ω, p) = iωΞ(ω, p) ∼ ω 2ν 0 + ω 2ν + + ω 2ν − . . .
with V 0 proportional to the chemical potential µ.
Given, however, that the plasmon frequency can be quite high in general (of order of the chemical potential) it is not clear whether these powerful deep infrared scaling forms valid near ω = 0 are still appropriate. Pragmatically we find that a simple power law form, fits the numerical results for the charge susceptibility at fixed p quite well.
All along the physical meaning of this holographic "quantum critical sector" has been a mystery. In the holographic community the focus has been entirely on zero momentum given the technical simplifications and the focus on measurable transport properties. It is however quite instructive to consider finite momentum. It is a property of the Fermi liquid as well. As we discussed in Section II, its optical and density responses at finite momentum are characterized by a "Drude sector" (zero sound), augmented by the incoherent contributions from the Lindhard continuum. As we noted below Eq. The real singular aspect of the Fermi-liquid is the upper bound of the Lindhard contin-uum, ω max = p F |p| m + p 2 2m , responsible for its disappearance at p = 0. This is intimately linked to the free quasiparticle kinematics of the FL fixed point. Dealing with the "unparticle" physics associated with the "strongly interacting" holographic strange metals such an upper bound is completely unnatural. Recognizing this unique particle feature of the free fermion fixed point, it is obvious that in anything else (the "non Fermi liquid") the continuum is expected to survive at zero momentum and finite frequency. This has the the ramification implied by the universal linear response theory of Section II that the long wavelength plasmon will have a finite lifetime. This we will now show by explicit computation in holographic strange metals as a proof of principle in expectation that two-sector quantum critical strange metals of this type are realized in nature.

METALS.
We first address how to incorporate the effects of the long range Coulomb interaction V p in the nominally neutral holographic strange metal. We do so in exactly the same way as reviewed in Sec. II, by coupling the density as part of a U (1) current to the electromagnetic field in the boundary, integrate out the Maxwell potential followed by the non-relativistic limit to obtain the partition function (c.f. (9)), This is known in the context of holography as a double-trace deformation -it is quadratic in the "single trace" operator n[Ψ] [63]. 3 The holographic counterpart, the on-shell action, must be deformed correspondingly. Keeping in mind that the subleading coefficient b(p) in (16) is responsible for the expectation value of n(p) , one adds the "double-trace" boundary term [64] in terms of b(p), This extra term modifies the boundary conditions for the gauge field A 0 (p; u), which for α = 1 take the mixed form [63] (Note that for V p = 0 it reduces to (17)) The perturbative source in this case is not directly related to the leading coefficient a(p) as in Eq. (17) and therefore the variation with respect to A 0 (15) brings down extra contributions from "double-trace" boundary term (33). Keeping in mind that for a given horizon boundary conditions b(p)/a(p) = χ 0 (ω, p) (20), one arrives at the result for the "dressed" charge density two-point function [13,63,64] which coincides exactly with the RPA formula Eq.(10). In Appendix B we discuss the generalization of this expression to the fully relativistic case.
A. The quantum critical sector and the plasmon.
To elucidate the gross features of plasmons in holographic strange metals let us consider for the first example the most familiar case in the literature, the z → ∞ locally quantum critical metal, dual to the simple charged black hole of the Einstein-Maxwell theory in the bulk. Although the strict limit is pathological as it has a finite ground state entropy (see Sec. III), this has been studied in particular because of its comparatively simple bulk structure.
With the number density correlator χ 0 (ω, p) computed numerically in the previous Section in hand, we evaluate Eq. (35) and the result for the imaginary part of χ(ω, p) is shown in Fig. 5. We have chosen a value of the electromagnetic coupling e 2 = 2.5 in the Coulomb potential V p = e 2 /p 2 such that the plasmon frequency is a reasonable fraction of the chemical potential, of the kind that is typically encountered in condensed matter systems (e.g., in the cuprates ω p 1 eV for a bare Fermi energy 2 eV). As explained, we indeed recover that the sound pole of the neutral system (Fig. 3) gets promoted to a finite energy (plasmon frequency) in the long wavelength limit. It should now not be surprising that in contrast to the sound pole the plasmon is damped substantially already in the long wavelength limit: this is the effect of the "quantum critical damping" of the plasmon, decaying into the second quantum critical infrared sector.
To highlight the role of the quantum critical infrared sector, we can artificially exclude its continuum contribution from the holographic neutral density correlator χ 0 (ω, p). As we discussed in the previous section, the pole structure of χ 0 (ω, p) can be obtained from the study of the quasinormal modes shown on Fig. 3, right panel. The dominant contribution comes from the pole corresponding to zero sound and we can use this to construct a "zero sound pole-only" density response with no continuum contribution Ξ pole-only (ω, p) = 0.
Substituting the pole-only density response into the RPA formula (35), and computing the spectral weight of the corresponding charge density correlator, one gets the result shown on Fig. (5), right panel. One immediately sees both our predictions verified. Removing the quantum critical continuum (a) shifts the gap at zero momentum (similar to "undoing" the effect of interband transitions in conventional systems), and (b) reduces the decay width of the plasmon dramatically and eventually kills it completely at zero momentum. At p → 0 • e 2 =0.5 • e 2 =1.25 only the quantum critical continuum Ξ QC (ω, p) contributes to the width of the plasmon.
Similar to the case of zero sound, discussed above, one can understand the plasmon response χ(ω, p) in terms of the pole structure of the "RPA-dressed" correlator in the complexω plane. These "plasmon QNM modes" are different from the formal quasi-normal modes of the black hole as one uses different boundary conditions (34). Nevertheless, they are closely related and we use similar technique to determine them; see Appendix A for more details. In Fig. 6 we show the momentum dependence of several sets of plasmon QNM modes, obtained for different values of the electric charge e 2 . Again, it is immediately seen that, unlike the sound poles, the plasmon poles do not approach the origin as momentum tends to zero.
Instead they halt at points with finite real and imaginary parts, which correspond to the plasmon gap and width at zero momentum, respectively. The absolute values of both real and imaginary parts grow with e 2 , which is again expected from the RPA result (11). -Im( χ) There is one more way that the quantum critical sector reflects itself in the plasmon response. Given the presence of other degrees of freedom than sound at finite frequency, the single pole approximation with its hydrodynamical damping ∼ Γ p falls short of capturing the physics well. This reflects itself also in the line shape of the plasmon: it is no longer a Lorentzian indicative of a long-lived mode with a finite lifetime, see Eq. (11). The effect of the finite imaginary part of Ξ(ω, p = 0) on the line shape is illustrated in Fig. 7. The dashed lines show the Lorentzian peaks corresponding to the position of the leading plasmon mode ω p (p) in the complex-ω plane whereas the actual response is shown in the continuous line. For large electromagnetic coupling e 2 = 2.5 the effect becomes especially pronounced. In particular the peak is now clearly asymmetric and has a long heavy tail towards the UV.

V. PLASMONS AND THE FAMILY OF HOLOGRAPHIC STRANGE METALS.
As we discussed in Section III, by tuning the potentials in the bulk EMD theory we can study a large variety of holographic strange metals, characterized by vastly different scaling behaviors in terms of the dynamical critical exponent z and the hyperscaling violating parameter θ (28). In this section we explore this wide parameter space and study the observable effects in the charge density response as the features of the quantum critical sector change. In principle we can compute the density response in holographic strange metals for any finite values of (z, θ). Two subfamilies are of particular interest, however. In the family with θ = 0 one can study the clean dependence on the dynamical critical exponent z. Despite its limited applicability, we can already infer from the scaling analysis Eq. (28) that decreasing z will decrease the contribution from the QC sector (recall that ω µ.).
Studying the dependence on z is therefore a more fine-grained test of the QC contribution to the density response. The other subfamily of interest is the set where z → ∞, θ → −∞ with θ/z = −η fixed. This subfamily describes locally quantum critical theories (z → ∞) with tunable scaling of the entropy s ∼ T (d−θ)/z ∼ T that in this case the contribution of quantum critical continuum Ξ is small and the plasmon decay is essentially completely governed by the hydrodynamical-like damping Γ(p). Indeed we observe in Fig. 8 that for all these families the width of the plasmon is almost zero in long-wavelength limit and grows with momentum. We do see, however, that the value of Γ(p) does depend on the scaling features of our model. It can be understood in the following way. At low temperatures decreasing z from z = ∞ or increasing θ not only suppresses the contribution from the quantum critical sector to Ξ(ω, p) but also to Γ(p). In holographic models this is computationally clear as the diffusivity is controlled by a black hole horizon in the scaling geometry. The computed physical response at finite p reflects this as the plasmon width clearly increases as z grows (Fig. 8-(a)). On the other hand in the family with z, θ → ∞ (Fig. 8-(c)) the plasmon decay rate increases with reducing the absolute value of η = −θ/z. In essence for larger z and smaller θ, the plasmon has more channels into which it can decay and its width grows.
Next we examine a larger coupling e 2 = 1.25, see right column of Fig. 8. Here the dressed response acquires a sizable contribution from the quantum critical continuum Ξ(ω, p), see Eq. (12). This is particularly pronounced at zero momentum, whereΓ → 0 and the width of the plasmon is completely dominated by Ξ(ω, p). In this long wavelength case we can study the continuum contribution very cleanly. As discussed above in Sec. III we expect Ξ(ω, p) to be strongly dependent on the z and θ exponents (28). However, rather surprisingly, we observe that the widths of the peaks at small momenta are not very sensitive to our tuning parameters. There are only small differences which can be seen by superposing the peaks, see Fig. 9.
From a more careful examination we understand why this sensitivity is not present in the response. First we note that the expression for the width of the plasmon at zero momentum, Γ(0) ∼ −AIm(Ξ)/ω, receives contributions from both the value of the QC continuum at plasmon frequency ω p and the residue of the sound pole A. In order to figure out the interplay between these contributions we extract these parameters from a fit to the neutral "undressed" density susceptibility χ 0 (ω, p) in the different models under consideration. We use the same ansatz (2) and obtain the fits as shown on Fig. 4. The values we obtain for A and the QC power law exponent ν at small momentum p = 0.02µ (see ansatz (31)) are shown on Fig. 9. When we substitute these fit-obtained values in the expression for the dressed correlator χ(ω, p) (11) we can extract the position and the width of the plasmon pole by finding the zero of the denominator numerically. These we show as Re(ω * p ) and Im(ω * p ) in the table of Fig. 9. We see that these values describe rather well the observed width of the peaks in the directly computed full correlator χ(ω, p), shown on the right panel of Fig. 9, but also that they barely change as we change z or θ. This confirms that our ansatz correctly takes into account this low sensitivity issue.
Inspecting more closely one sees that it is true that, as expected from Eq.(28), the continuum power law exponent ν depends strongly on the background scaling features. 4 However, the sound pole residue A also changes. It contributes oppositely in such a way that it stabilizes the plasmon width. This is surprising, though in a qualitative sense one should know that A and ν are not entirely independent. There is a nontrivial connection between A and ν arising from the spectral sum rule for the number density correlator (where we use the ansatz (11) Even though Eq. (38) does not prescribe the exact relation between A and ν, it correlates the changes in Ξ(ω, p) ∼ ω ν to changes of A. Nonetheless, that they counteract each other so closely is surprising.

VI. CONCLUSIONS
As more often, holography has a remarkable capacity to suggest unusual general physics questions of special relevance to condensed matter [13] that in hindsight have simple answers.
The present story is an illuminating illustration.
The impeccable, long-lived nature of the plasmon "particle" is at the heart of the large engineering pursuit of plasmonics. It should be clear by now that this is actually rooted in the special, singular free-particle nature of the Fermi-liquid fixed point. The existence of sound is a different matter -it is a consequence of conservation laws and although there is no mathematical proof that it has to exist at zero temperature it may well be ubiquitous.
By the same token, given the universal nature of the electromagnetic linear response theory, this massless sound will be promoted to the plasma frequency at zero momentum. However, 4 Note, however, that the exact value of the fitted power law which we find doesn't coincide with the prediction of Eq. (28). The reason is that the former is fitted in the range of frequencies of order ω p (see Fig. 4), while the latter is only valid in the limit ω → 0.  width⌃ from Eq. (??). In the tables shown in Fig. 9 we see that indeed the width of the plasmon, estimated from Im ! ⇤ p , is fairly independent of z, ✓. The small di↵erences observed for di↵erent anomalous exponents are more easily seen by collapsing the plasmon peaks as shown on the right plots of Fig. 9. for the plasmon to have an infinite lifetime at this finite energy, it has to be completely isolated from all other energy eigenstates that couple to charge density. There is no general, symmetry based argument insisting that this is a natural condition. The very special nature of the Fermi liquid plays the primary role here. Given its free fermion IR fixed point the only states available are formed from the perfectly free fermionic quasiparticles and according to textbook kinematics such states have vanishing energy in the long wavelength limit. It follows immediately that, in absence of this special kinematical protection, any strongly interacting IR fixed point of fermion matter (the non-Fermi liquid, by definition) will give rise to a plasmon that will in principle decay in the long wavelength limit already.
This can be seen explicitly in computations in holographic non-Fermi liquid strange met-als. In these holographic strange metals the other "non-hydrodynamical" charge carrying states are organized in a particular elegant fashion, ruled by scaling principle. But the basic motive is completely general. Do note that, as we stressed in the Introduction, the "quantum critical" sector has a very different meaning in the holographic strange metal context than in the case of quantum criticality originating in a quantum phase transition involving the disappearance of an order parameter. We expect, however, that in these latter cases the plasmon will also decay. A case in point is the marginal Fermi liquid [65][66][67] and variations thereof. There a polarization propagator is envisioned which is not dissimilar from the holographic one, but instead it originates in the quantum critical fluctuations associated with the order parameter. It would be quite interesting to study the fate of the plasmon in this Herz-Millis framework, especially for the case of momentum conserving systems.
Zooming in on the particulars of the holographic plasmons, a quite surprising outcome is that the plasmon (as well as zero sound) is remarkably robust up to very large momenta.
We find that in all cases we looked at, it is still underdamped up to momenta |p| ∼ 1 2 µ -(half) the Fermi momentum in a condensed matter interpretation. The plasmon at such momenta still dominates to a degree that the Ξ(ω, p) continuum is concealed from direct observation. In weakly interacting Fermi liquids the plasmon enters the Lindhard continuum at momenta around this scale and is rapidly obliterated by Landau damping. There are basic scenarios where sound is less prominent, e.g. a generalization from static Coulomb interactions to relativistic electrodynamics case can make it more overdamped as we show in Appendix B. However, taking our lesson at heart, that the quantum critical sector is the generalization of the Lindhard continuum, the plasmon should be overdamped from the start. To understand why it is not, is hard to answer. Why this "hydrodynamic response" works far better than it should -an observation also made in other holographic studiesis a UV sensitive affair: the density operator of holography is associated with large central charge CFT physics, having a-priori nothing to do with the microscopic electrons and atoms of the condensed matter systems. Having the latter in mind, it is a matter of high priority to look for holographic set-ups that are more flexible in this regard. A closely related issue is the incomplete understanding of how the quantum critical sector evolves at finite momenta.
A starting point is [16,17,68], but these are pressing issues for future holographic research.
Last but not least, what is the relevance of this affair for the condensed matter laboratory? We alluded to the timely nature of the subject in the Introduction, while presently major investments are made in condensed matter laboratories to improve the electron-loss instruments. To close the gap with experiment the highest priority is to do justice in holography to the circumstance that the non-Fermi liquid electron systems in the laboratory are invariably characterized by very strong lattice potentials. Undoubtedly, the strong Umklapp scattering is a key factor in shaping the charge reponse. Although this presents quite a technical challenge, this hurdle has to be taken in order to arrive at a meaningful dialogue with experiment. The equations of motion for the background fields following from the action (21) are [58]: where we set 2κ = q 2 = 1. The Maxwell and scalar equations are: These equations may be solved numerically using the shooting method in a particularly simple gauge [58,69]. However, here we choose to use the DeTruck method to solve the equations as a boundary value problem [10,[70][71][72][73], which will facilitate our analysis of the linear fluctuations. We use the metric ansatz for d = 3: The DeTurck reference metric is just a Reissner-Nordström black hole which one gets when setting all the ansatz functions to unity. The temperature is set by the parameterμ and we can setμ = µ without loss of generality. We consider T /µ = 0.02 throughout the paper, which is small enough and accessible for our numerics. We solve the resulting nonlinear elliptic equations of motion via Newton's method using the pseudospectral collocation approximation with 60 points on the pseudospectral grid. More details on the algorithm can be found in [74].

Linear response equations
In order to evaluate the density-density correlator we first introduce the linear perturbations of all the fields with finite frequency and momentum: The analysis of the linearized equations shows that there are 2 decoupled composite modes: longitudinal : {δg tt , δg tx , δg tz , δg xx , δg xz , δg yy , δg zz , δA t , δA x , δA z , δφ} transverse : {δg ty , δg xy , δg yz , δA y } (A7) We use the DeDonder gauge in δg µν and Lorentz gauge in δA µ [10,72]. This brings all the equations of motion for the 11 fields in the longitudinal channel to the convenient second order form.
Using the pseudospectral collocation approximation on the same grid as the background solutions, we recast the linearized equations as a linear algebraic problem where the right hand side corresponds to the source in (17). This equation is efficiently solved by the LinearSolve routine of Wolfram Mathematica [75], giving the vector representation of the linearized solutions, which we use to read off the two-point function (20).
The quasinormal modes are the solutions to the Sturm-Liuville problem and are simply the eigenvalues of the matrix M(ω) (This is similar to the method used in [76]).
We also get the plasmon modes in a similar way, but one should note that in this case we use the different boundary conditions (34), and therefore the matrix of the boundary value problem gets modified [74]. The plasmon modes are therefore the eigenvalues of the deformed matrix M Vp (ω).

Appendix B: From Coulomb interaction to full dynamical electromagnetism
The non-relativistic approximation and reduction of full dynamical electromagnetism to the static Coulomb potential only (9) clearly breaks Lorentz invariance. This is justified because the Lorentz invariance, implicit in the general relativistic description of holographic strange metals, is set by the emergent speed v F in the linear dispersion relation of the gapless modes at low energies and not the fundamental speed of light in vacuum c. Nevertheless, it is straightforward to include dynamical electromagnetism by either the proper photon propagator V µν = iD µν = η µν /p 2 or a generalized one that takes into account the difference between v F = 1 and c: V 00 = −e 2 c 2 /(c 2 p 2 − ω 2 ), V xx = e 2 c 2 p 2 −ω 2 . The full electromagnetic interaction J µ V µν J ν will couple both charge densities n ≡ J 0 via electric Coulomb force, but also the spacial components of the currents J i via the magnetic Lorentz force. The RPA formula (10) becomes matrix valued and gives   χ 00 χ 0x where Π µν = J µ J ν e 2 =0 and detΠ = (Π 00 Π xx − Π x0 Π 0x ). In the non-relativistic limit c → ∞ this reduces to the previous results as should be.
For completeness we give an example of the plasmon response with dynamical electromagnetism with c 2 /v 2 F = 10 3 , 200 and e 2 = 0.5 in Fig. 10. At p → 0, V xx has a similar effect to a stronger V 00 dressing, namely the gap increases and the plasmon peak is damped more significantly. There is however, a qualitatively different feature: the on-shell photon state, which is refelcted as a zero of the density-density response at ω = c|p|. As momentum increases there is another qualitative difference: the gap and peak damping initially decrease as momentum increases until the effect of V xx is negligible and the response is controlled again by V 00 alone. c 2 p 2 −ω 2 are considered. Left column: θ = 0, z → ∞, and right column θ = 0, z = 2. The main differences observed with respect to the case of V xx 0, discussed in the main text, are the strong damping of the plasmon for larger momenta (green lines) and the zero of the spectral density, associated to the on-shell photon at ω = c|p|. We take the electromagnetic coupling e 2 = 0.5 and T = 0.02µ.

Appendix C: 2D-Coulomb interaction
In layered systems the Coulomb interaction is effectively two dimensional: V (|p|) = e −λ|p| /|p|. This interaction leads to the standard square-root dispersion relation, see Fig. 11.