# Directly characterizing the relative strength and momentum dependence of electron-phonon coupling using resonant inelastic x-ray scattering

###### Abstract

The coupling between lattice and charge degrees of freedom in condensed matter materials is ubiquitous and can often result in interesting properties and ordered phases, including conventional superconductivity, charge density wave order, and metal-insulator transitions. Angle-resolved photoemission spectroscopy and both neutron and non-resonant x-ray scattering serve as effective probes for determining the behavior of appropriate, individual degrees of freedom – the electronic structure and lattice excitation, or phonon dispersion, respectively. However, each provides less direct information about the mutual coupling between the degrees of freedom, usual through self-energy effects, which tend to renormalize and broaden spectral features precisely where the coupling is strong, impacting ones ability to quantitively characterize the coupling. Here we demonstrate that resonant inelastic x-ray scattering, or RIXS, can be an effective tool to directly determine the relative strength and momentum dependence of the electron-phonon coupling in condensed matter systems. Using a diagrammatic approach for an 8-band model of copper oxides, we study the contributions from the lowest order diagrams to the full RIXS intensity for a realistic scattering geometry, accounting for matrix element effects in the scattering cross-section as well as the momentum dependence of the electron-phonon coupling vertex. A detailed examination of these maps offers a unique perspective into the characteristics of electron-phonon coupling, which complements both neutron and non-resonant x-ray scattering, as well as Raman and infrared conductivity.

## I Introduction

The ability to characterize excited states and fundamental excitations in solids remains one of the forefront challenges in condensed matter physics. This is particular evident in materials that display emergent phases, such as unconventional superconductivity, that stem from intertwined lattice, charge, spin and orbital degrees of freedom Fradkin2015 (). Typically simple pictures for excitations derived from well-controlled perturbation expansions do not capture the richness of these materials due to the lack of a small parameter around which to expand. This is also evident particularly for excitations around a quantum critical point in correlated materials, where the low energy degrees of freedom are presumably intertwined. The lack of knowledge of how to truly quantify the strength of couplings for all momenta throughout the Brilloiun zone (BZ) has been one of the major roadblocks to understanding the complex phase diagrams that typically emerges near phase boundaries between various electronic, magnetic, or lattice instabilities.

In the last decade resonant inelastic x-ray scattering (RIXS) has enabled seminal progress in the understanding of fundamental excitations in correlated materials Ament (). Exploiting the role of strong spin-orbit coupling in core levels in Cu -edge RIXS, it has been shown recently that spin flip excitations can be probed across a wide variety of materials, complementing RIXS studies of charge transfer, orbital, and bi-magnon excitations down to a now available resolution better than 100 meV GG (); AS (). Perhaps more intriguingly from the perspective of a single tool that can provide complete characterization of charge, orbital, spin, and lattice exciations, excitations off of the elastic line oxygen -edge that occur at multiples of optical phonon energies have become visible in one dimensional edge-shared CuO systems LeePhononsI (), while phonon side bands have emerged off of excitations measured at the Cu -edge LeePhononsII (). This progress have been made with a resolution of about 120 meV at the Cu -edge obtained already ten years ago at the Swiss Light Source (SLS) Ghiringhelli2006 (); AS (). The possibilities expanded further less than one year ago by the new standard of 35 meV at the Cu -edge at the European Synchrotron Radiation Facility (ESRF) Brookes2016 (). Since the phenomena of emergence involves necessarily the intertwined nature of different degrees of freedom in the excitation spectra, knowledge of the changes of excitation spectra is a goal that is well-addressed in present and future RIXS experimental configurations.

Specifically, one of the outstanding issues in condensed matter physics is to determine the strength of coupling of electrons to lattice excitations across the entire BZ. This has important ramifications to our understanding of emergent phenomena such as superconductivity and density wave order. For conventional isotropic superconductors, the momentum dependence of the electron-phonon coupling is rather irrelevant to pairing as all modes contribute to the Cooper instability. However, this is not the case in unconventional (sign-changing) anisotropic superconductors Bulut (); OKA (); Johnston (). While electron-phonon coupling at large momentum transfers that connects momentum points with different signs of the superconducting gap leads to a suppression of superconductivity, scattering that involves small momentum transfers always enhances superconducting pairing JJLee2014 (). Likewise, it has now become accepted that for dimensions greater than 1, the momentum dependence of the electron-phonon coupling largely controls the ordering wavevector CDW (). Therefore a determination of is highly desirable on rather general grounds.

While Raman and optical measurements can extract the coupling at zone center buckling (); RamanRMP (); neutron (), a reliable method for extracting the strength of lattice couplings for all momenta to date has been inferred indirectly from angle-resolved photoemission (ARPES) experiments, and more directly from inelastic neutron or high-energy non-resonant x-ray scattering, which couple directly to lattice vibrations Johnston (); SpecRev (); AREV (); IXS (); Reznik (); BiPhonons (); SJJacoustic (). While ARPES probes electrons at well defined electron momenta in the BZ, a sum over phonons at all bosonic momenta contributes to renormalization effects detected as ”kinks” in dispersion, or abrupt changes of spectral linewidths at the energies of the lattice modes. This precludes as a practical manner a momentum-resolved way to detect throughout the BZ.

On the other hand momentum dependent coupling can be determined from lineshape fitting of x-ray or neutron scattering. However, these techniques present a fundamental challenge to the extraction of lattice coupling itself. Due to the scattering process which couples to the atomic degrees of freedom, lattice vibrations are particular sharp and easy to detect when the lattice coupling to electrons is weak. In contrast, when the coupling is large, as seen for example in bond-stretching modes in the cuprates or near charge density wave instabilities, phonon lineshapes become broad in both energy and momentum which interfere with other phononsReznik (). Thus the region of energy and momentum where the lattice coupling is strongest tends to be also the region where a direct extraction of the coupling is most difficult. Therefore one desires a complementary method that can be used to determine coupling in this region.

As suggested in earlier studies LeePhononsI (); exprixs (); EPL (); Yao (), RIXS offers a unique insight into the momentum dependence and magnitude of electron-phonon coupling. In this paper we examine this selectivity of the light scattering process to determine how and which phonons can be resonantly excited during a resonant x-ray process, and the resulting information that can be obtained about the relative strength and momentum dependence of electron-phonon coupling. We primarily consider an -like resonant x-ray scattering where a core electron is excited into the conduction band, and assume that the most interesting phonons are those that couple to the conduction/valence electrons/holes. We will neglect effects from direct phonon coupling to core levels, which would primarily only affect the intermediate state (absorption) profile. We particularly highlight the role of polarization and selection rules of the light-electron-phonon coupling, which are already well-characterized only at zone center in Raman measurements. Specifically we find that RIXS can provide detailed clues in regions of momentum space with strong electron-phonon coupling with atomic specificity, to further complement neutron and IXS measurements.

## Ii Expressions for phonon contributions to RIXS

Since we are interested in the generic ways in which phonons can be coupled in RIXS, we consider weak electron-phonon coupling and examine the one-phonon contribution to RIXS (Figs. 1 and 2). While other considerations involving the ratio of main intensities to satellites offer a way to quantify electron-phonon coupling for generic model Hamiltonians LeePhononsI (); EPL (), our approach follows an extension of resonant Raman scattering from phonons, extended to arbitrary momentum transfers and different x-ray resonant edge processes. We examine the interplay of light polarization configurations with specific momentum dependent electron-phonon coupling for different classes of coupling (for example, deformation, piezoelectric, and electrostatic), for some oxygen modes common to perovskites. Here we neglect the role of Coulomb interactions, both among and between the valence electrons and the core hole. While the role of Coulomb interactions is crucial for indirect RIXS processes Ament (), and for couplings at BZ center, for larger momentum processes it plays only a minor role Johnston (); tpd (); multiband ().

The phonon branches are characterized by a frequency for each mode as well as a generalized coupling to electrons . Here denotes the coupling of phonon mode to an electron in band carrying momentum scattering into band with momentum . The total RIXS response contains two contributions from these two diagrams each of which depend upon where and .

For the bare loop contribution (Fig. 1) we have

(1) |

where denotes the imaginary part and denote the conduction, core electron, respectively, appropriate for a direct RIXS transition. The matrix elements depend generally on the incident and scattered photon momentum , energy and polarizations , respectively, including the light dipole couplings involving the resonant photoexcitation of a core hole with momentum into the band with momentum .

As a concrete example, our target will be an understanding of lattice coupling in the cuprates, to which end we consider an Cu edge scattering process where a Cu core electron is photoexcited into the valence band at a resonant energy around 931 eV (Cu -edge). However, the formalism will work equally well for other resonant excitations involving core electrons and valence states, such as the oxygen -edge.

We only consider particle-hole excitations in the Cu-O hybridized band which cuts through the Fermi level. This allows us to write with the eigenenergies, the Cu projected character of the band, and sign. Expressions for these functions can be determined by Wannier downfolding in density functional based approaches, or via tight-binding parameterizations for free electrons, or exact diagonalization studies for correlated models. In addition expressions for the momentum dependent electron-phonon coupling constants can be determined via largely the same methods.

The localized core Cu orbitals are represented by having an energy measured between the bottom of the Cu level and the Cu level, and a phenomenological core-hole lifetime . One can show that such a resonance approximation for the propagator (which neglects the contributions from other edges) works well for RIXS at either the Cu or edge Marra2012 (); Jia2015 ().

The limit leads to the collapse of the frequency integral. After substitution of the expressions for the Green functions one obtains

(2) |

By taking the limit we neglect the multi-particle interactions within the electrons, focusing attention on a simple understanding of the role of lattice coupling. For fully interacting electrons, one still can proceed by taking Eq. (1) in place of (2).

As the main focus of this study, the phonon contribution (represented diagrammatically in Fig. 2), can be cast as

(3) |

with indices running over phonon modes. In accordance with Fermi’s golden rule, this contribution involves the product of two matrix elements and the phonon density of states represented by the phonon propagator . In the limit of large core-hole energy , the matrix elements can be written as

(4) |

and

(5) |

Using the previous Green functions definitions these equations become

(6) |

and

(7) |

Here we do not assume free electrons, and keep finite.

These matrix elements carry the resonant enhancement of phonon scattering via coupling to electrons. Importantly, they show that the resonance process itself imparts a projection of the electron-phonon coupling to the portion of the coupling involving the photo-excited valence electron, in this case copper. If we instead focused on the oxygen -edge, the matrix elements would involve the projection onto the oxygen contribution to the hybridized valence band. This is one important facet that already distinguishes RIXS from other measurements, and allows RIXS from phonons to have an element selectivity to not only the particle-hole excitations, but also the momentum-dependent phonon coupling to them.

Finally, for the phonon contribution one obtains

(8) |

It can easily be seen that there can be several sources of enhancement from the resonant process. First, if both the and states are occupied, no scattering can occur since a core electron is Pauli blocked. Therefore at least one or both of those states must be unoccupied. However the largest contribution results when only one of the momentum states is unoccupied so that an occupied electron state can refill the core hole.

Since we have assumed that the phonons are eigenstates indexed by , the loss peaks are determined by the phonon mode energies and their dispersion with momentum, multiplied by overall factors that represent the resonant excitation process. This demonstrates how the electron-hole pairs created by the process couple to those phonons in a symmetry-dependent way that is sampled by the incident and scattered polarization orientations. This is analogous to the way in which photon polarization selection rules can be used to probe unconventional superconductivity RamanRMP (); tpd ().

## Iii General behavior of RIXS from Phonons

We can examine a few insights from the momentum-dependent electron-phonon coupling matrix element. From Eqs (6-8), it is clear that the overall intensity of the phonon contribution to RIXS is governed by the strength and momentum dependence of the electron-phonon coupling . We first consider general couplings that determine its momentum dependence.

For electron-phonon coupling of the deformation type, which involves lattice motion that modifies the kinetic energy of electrons via orbital overlaps on adjacent ions, the electron-phonon coupling depends strongly on transfered momentum , being largest for large momentum transfers and vanishing in the limit of Mahan (). These phonons, such as the longitudinal acoustic modes or the bond-stretching modes in the cuprates, fail to appear in Raman or optical conductivity and likewise would not contribute to RIXS for small momentum transfers via Eq. (6-7). Moreover, the largest RIXS intensity from those phonons will occur for large momentum at the BZ boundaries where deformation electron-phonon coupling is largest.

On the other hand, for couplings that are driven by modulations of the electrostatic energies via atomic motion, the momentum dependence is biased towards small momentum transfers. These include out-of-plane -axis oxygen vibrations in transition metal oxides, which show up in Raman and optical measurements and likewise in RIXS at small wavevector transfers Johnston ().

Also any phonons which do not transform according the full irreducible point group symmetry of the lattice will likewise vanish at zone center provided that the dipole matrix elements do not depend upon the magnitude of the momenta involved in scattering. Since the dipole matrix element involves photo excitations of core levels and is thus quite local, this would be a very good approximation. As a consequence, scattering would be suppressed from zone center phonons having or symmetries in systems. This is quite different from resonant optical Raman scattering where electrons can be non-locally excited out of the valence band into fully itinerant bands and thus depend on the momenta involved in the scattering process.

We can show that this is the case for Cu edge RIXS via the following consideration. Following Ref. Ament, the RIXS vertex in the dipole approximation (here for the process from the initial to intermediate state of RIXS) is

(9) |

where () is the incoming photon momentum (polarization), is the lattice site, and is the position operator of electron at site .

This can be written in the second quantization form (for the transition metal ion edge) as

(10) |

where () are the quantum numbers of the electron in the () orbital. Here the RIXS dipole operator matrix element does not depend on the momentum of the valence or core electrons, since it follows from the expression

(11) |

(where and are quantum numbers of the electrons in the or orbitals, are lattice sites, and are Wannier states centered on atom ) and the fact that

(12) |

This is due to the localized nature of the core hole. This property is the one that distinguishes RIXS from optical Raman scattering.

In this paper where we are interested only in RIXS in the energy range of single phonons, we keep only the orbital and neglect the other Cu 3 orbitals which would appear if we consider larger energy transfers at the excitation lines. Moreover, we are solely interested in the RIXS transitions for which the spin of the electrons is conserved: this is possible when choosing a particular combination of the polarization vectors and Groot1998 (); Veenendaal2006 (); AmentPolarization (); Haverkort2010 (); Marra2012 (). Thus, assuming that the latter indeed takes place, we can neglect the spin quantum number carried by the electron and we are left with the dipole matrix elements which depend solely on the orbital quantum number: where . We note that in RIXS experiment for cuprates, different momentum transfers are usually realized by changing the relative angle between the incoming photons and the sample. Thus, the incoming and outgoing photon polarizations may effectively depend on the momentum transfer AmentPolarization (). We do not consider this effect in the current manuscript.

## Iv An example: 8-Band Cu-O model

We can evaluate RIXS for a simple model as given in Refs. Johnston, and SJJacoustic, which derived momentum dependent electron-phonon couplings for Cu-O planar longitudinal acoustic, bond stretching and bond bending modes, as well as axis apical modes. Besides accounting for the phonon eigenmodes, the model consists of 2 Cu (3 and ) orbitals octahedrally coordinated and hybridized with O orbitals that are -bonded to the Cu 3 orbital, plus the apical O orbital hybridized with Cu . Oxygen hopping is taken into account. We finally add the 3 Cu core levels from which electrons are resonantly excited.

For the phonon excitations, we specifically focus on the oxygen modes common to all perovskite oxides. These phonons typically couple to CuO electrons either electrostatically in the case of axis phonon modes or via bond deformation for in-plane stretching modes. The forms for the electron-phonon coupling depend strongly on which process is dominant, and further can depend strongly on either fermion momentum or phonon momentum transfer . In Ref. Johnston, these modes were examined and expressions for their couplings were derived from general symmetry grounds. Their forms can be expressed as:

(13) |

with and the lattice constant set to 1. Here denotes the longitudinal acoustic modes, the Cu-O in-plane bond stretching modes, the -axis oxygen modes involving charge transfer between Cu-O and O-O, respectively, and the apical oxygen modes. The prefactors for mode are rather weakly dependent on momenta and are generally set by the orbital content of the bands that are coupled to the lattice modes. Here we limit focus on them to solely setting an overall magnitude relative to each other. For more details, the reader is directed to Ref. Johnston, . Finally we take the apical and modes to be dispersionless phonons of frequencies 85, 40 and 35 meV, respectively, and the dispersions of the acoustic and breathing modes in eV as and to match the dispersions in Bi-based bi-layer cuprates BiPhonons ().

We neglect any momentum or angle dependence of the polarization vectors and take them to be constants such that Eqs. (1) and (8) can be simply integrated numerically. Here we have taken a simple downfolded three-band Cu-O tight-binding model with Cu-O hopping eV and charge transfer 0.9 eV to mimic the anti-bonding band structure of the cuprates, doped for an overall filling neutron (); buckling (). Additionally we have set since the resonance is quite broad due to the short core-hole lifetime ( eV). The width of the all phonons is set to 5 meV, and eV is taken as the width of the -function in Eqs. (1) and (8) to phenomenologically mimic inelastic relaxation in the particle-hole continuum. The resulting RIXS response is shown in the Figs. 3-9 for momentum cuts along the BZ -axis. Since the phonons enter additively, we plot them individually as their intensities can be set from considerations given in the references cited above.

Fig. 3 shows the bare contribution from particle-hole excitations to RIXS from Eq. (1). The particle-hole continuum is resonantly enhanced, with the leading edge rapidly dispersing away from BZ center with the Fermi velocity. The overall continuum expands with increasing energy loss and momentum transfer due to kinematic conditions of particle-hole creation, with the overall strength and sharpness of the continuum being set by the inelastic relaxation . We note that the continuum dips down close to zero energy transfer around a momentum , reflecting the weak nesting condition of the band structure and Fermi surface for this momentum transfer.

Figs. 4-8 show the single phonon contribution to RIXS for the individual phonons listed above. For the acoustic phonon branch (Fig. 4) it is clear that the dispersion for RIXS follows the phonon dispersion, having an intensity which follows the momentum dependence of the electron-phonon coupling . The relative strength of the phonon contribution to RIXS compared to the bare response (Fig. 4) is set both by the overall electron-phonon coupling as well as the particle-hole relaxation . Therefore while the overall magnitude of the electron-phonon coupling cannot be straightforwardly determined solely from the single phonon RIXS signal, the momentum dependence of the coupling as well as the dispersion can be clearly read off directly from the RIXS intensity plots. In contrast, the magnitude of the coupling can be determine by the ratio of the phonon side bands when multiple phonon inelastic events are considered LeePhononsI (); LeePhononsII ().

Consideration of the interplay of electron-phonon coupling to the continuum as well as the overall symmetry between the electron-phonon coupling and the edge RIXS process can be seen clearly in Fig. 5 which shows the contribution to RIXS from the phonons (charge transfer between O and O in the Cu-O unit cell). Both the overall dispersion as well as the intensity vary nonmonotonically across the BZ. We first consider symmetry effects. As already noted, this phonon branch at zone center transforms as a lower representation of the point group, and therefore, due to the locality of the RIXS process, the phonon branch does not contribute to the RIXS intensity at zone center and grows for increasing momentum transfer. This is in strong contrast to either resonant optical Raman scattering or infrared conductivity for which the intensities are strong for this phonon and show clear Fano interference with the particle-hole continuum buckling ().

Moreover, a clear interference effect is seen when the particle-hole continuum crosses the phonon line. This occurs weakly near zone center and much more clearly around the weak nesting momentum . At both of these points the phonon contribution to RIXS inherits contributions from the particle-hole continuum and shows a strong momentum-dependent Fano lineshape, distorting along the lines of the soft continuum. The strength of this Fano effect is like that in Raman scattering, determined both by the strength of the electron-hole continuum at the frequency of the phonon and the magnitude of electron-phonon coupling, but in addition the range of momentum transfers available through RIXS gives a clear qualitative way of understand electron-lattice coupling nearinstabilities, allowing for the particle-hole continuum to interfere with the phonon in a narrow range of mode momenta.

This Fano effect is indicated also clearly in the phonon contribution to RIXS, shown in Fig. 6. Here the symmetry of the modes (involving charge transfer between Cu-O) allow for finite coupling at zone center, where the overall bare electron-phonon coupling is largest. However a clear drop in intensity is observed where the electron-hole continuum crosses the phonon energy. This intensity would add to the intensity of the bare continuum shown in Fig. 3. The momentum-dependent coupling falls off for large momentum transfers and gives vanishing intensity at the BZ boundary.

In Fig. 7, the Cu-O bond stretching phonons disperse downward away from zone center, and due to the strong momentum dependence of the bare electron-phonon coupling, the intensity grows as to be largest at the BZ boundary. The coupling’s momentum dependence minimizes the Fano effect near the BZ zone center, yet allows for a weak coupling around the weak nesting feature near where the intensity extends from the phonon line down to lower energy transfers.

In Fig. 8 the dispersionless apical oxygen contribution to RIXS is shown. The bare electron-phonon coupling is quite anisotropic in momentum space due to the charge transfer pathway between planar Cu and the apical oxygens, and reaches its strongest coupling for zone center. However, as in the case of phonons, the Fano effect is clearly observed at the intersection of the apical phonon with the main peak in the electron-hole continuum. The bare intensity would normally fall smoothly for larger momentum transfers, yet in contrast to the phonons and more like the phonons, the intensity remains large as the phonon line disperses across the weak nesting momentum before abruptly falling for larger momentum transfers.

As the electron-phonon coupling is only considered at lowest order, these contributions add to the overall RIXS signal, as shown in Fig. 9. One can clearly see that different phonon contributions are illuminated depending upon the overall magnitude of the electron-phonon coupling as well as the intersection of the phonon line with the electron-hole continuum.

Finally, Fig. 10 plots the total RIXS intensity map from both the bare electron-hole continuum and the contributions from all phonon branches considered. The figure demonstrates the rich overall spectra having a strong momentum dependence. A detailed examination of the structure of the RIXS maps offers an unique perspective into the characteristics of lattice-electron coupling in a way which complements neutron and non-resonant x-ray scattering, as well as Raman and infrared conductivity. Thus RIXS at low energies offers a direct insight into couplings to lattice modes.

We add a remark concerning the effect of screening via the long-range Coulomb interaction, which has not been included in these considerations. It is well known that Coulomb screening mediates charge backflow that enforces particle number conservation and can lead to strong modifications of the electron-phonon coupling for small momentum transfers provided the symmetry of the coupling coincides with the full symmetry of the lattice tpd (); multiband (); RamanRMP (). In addition in layer materials a screening enhancement can occur when the phonon energy crosses the effective 3D plasmon Johnston (). Both of these effects become prominent for small momentum transfers and do not effect the overall considerations given here for larger momentum transfers.

## V Summary and Connection to experiments

With the recent development of x-ray Synchrotron Light Source and facilities, two significant improvements have made RIXS a possible approach for directly addressing phonon excitations. First, RIXS for the transition-metal -edge has shown energy resolution around 100meV, enabling RIXS to study not only the specific phonon modes on top of the excitations or at the oxygen -edge directly off the elastic line, but also the momentum dependent one-phonon contributions from different phonon modes directly. Second, although previously polarization discrimination had been limited to only the incoming photon, the recently constructed ESRF RIXS facility has made fully distinguishable polarization measurement for both incoming and outgoing photons possible LeTaconPolarization (). Other end-stations also would allow for the full photon polarization discrimination in the near future. Full polarization discrimination measurements can disentangle spin-flip, or single (para)magnon, excitations when the incoming and outgoing photon polarizations are perpendicular to one another (or cross-polarization) and the non-spin-flip or electron-hole excitations when the incoming and outgoing photon polarizations are parallel with each other (or parallel-polarization).

With the recent development of x-ray Synchrotron Light Source and facilities, two significant improvements have made RIXS a possible approach for directly approaching of the phonon excitations. First, RIXS for the transition-metal - edge has reached at the ESRF an energy resolution down to 35 meV in the most difficult case (Cu-) among , enabling RIXS to study not only the specific phonon modes on top of the excitations or at the oxygen K-edge, but also the momentum dependent one-phonon contributions from different phonon modes directly esrf (). Second, although previously the phonon polarization discrimination is only for the incoming photon, recently constructed ESRF RIXS facility has made the fully distinguishable polarization measurement for both incoming and outgoing photons possible Ghiringhelli2006 () with the first application LeTaconPolarization (). With the polarimeter the present standard is 90 meV with a reduction of the counting rate by a factor of 10.

To summarize, we have derived and calculated the electron-hole and one-phonon contributions to direct RIXS, specifically at the Cu -edge, based on a simple diagrammatic approach. We have examined an 8-band Cu-O model with oxygen phonon modes, including Cu-O planar longitudinal acoustic, Cu-O bond-stretching and apical O modes. This model reveals that the momentum dependence of the electron-phonon coupling can be extracted from the intensity of RIXS excitations in the low energy regime, touting RIXS as a potential method to directly characterize the momentum dependent electron-phonon coupling. Our theoretical approach provides an intuitive picture of how different phonon modes contribute to RIXS and interacts with charge densities, making RIXS a potential new technique for electron-phonon coupling measurement, complementary to ARPES, inelastic neutron scattering and non-resonant inelastic x-ray scattering.

###### Acknowledgements.

This research was supported by the U.S. Department of Energy (DOE), Office of Basic Energy Sciences, Division of Materials Sciences and Engineering, under Contract No. DE-AC02-76SF00515, SLAC National Accelerator Laboratory (SLAC), Stanford Institute for Materials and Energy Sciences. K.W. acknowledges support from the Polish National Science Center (NCN) under Project No. 2012/04/A/ST3/00331.## References

- (1) Eduardo Fradkin, Steven A. Kivelson, and John M. Tranquada, Rev. Mod. Phys. 87, 457 (2015).
- (2) L. Ament, M. van Veenendaal, T. Devereaux, J. Hill and J. van den Brink, Rev. Mod. Phys. 83, 705 (2011).
- (3) L. Braicovich, J. van den Brink, V. Bisogni, M. Moretti Sala, L. J. P. Ament, N. B. Brookes, G. M. De Luca, M. Salluzzo, T. Schmitt, V. N. Strocov, and G. Ghiringhelli, Phys. Rev. Lett. 104, 077002 (2010); M. Dean et al., Nat. Mater. 12, 1019 (2013); M. Le Tacon et al., Nat. Phys. 7, 725 (2011); K.-J. Zhou et al., Nat. Comm. 4, 1470 (2013).
- (4) V. N. Strocov et al., J. Sync. Rad. 17, 631 (2010).
- (5) W. S. Lee et al., Phys. Rev. Lett. 110, 265502 (2013).
- (6) J. J. Lee et al., Phys. Rev. B 89, 041104(R) (2014).
- (7) G. Ghiringhelli, A. Piazzalunga, C. Dallera, G. Trezzi, L. Braicovich, T. Schmitt, V. N. Strocov, R. Betemps, L. Patthey, X. Wang and M. Grioni, Rev. Sci. Instrum. 77, 113108 (2006).
- (8) N. Brookes, L. Braicovich, and G. Ghiringhelli et al, to be published.
- (9) N. Bulut and D. J. Scalapino, Phys. Rev. B 54, 14971 (1996).
- (10) S. Y. Savrasov and O. K. Andersen, Phys. Rev. Lett. 77, 4430 (1996).
- (11) S. Johnston, F. Vernay, B. Moritz, Z.-X. Shen, N. Nagaosa, J. Zaanen, and T. P. Devereaux, Phys. Rev. B 82, 064513 (2010).
- (12) J. J. Lee, F. T. Schmitt, R. G. Moore, Y. T. Cui, W. Li, M. Yi, Z. K. Liu, M. Hashimoto, Y. Zhang, D. H. Lu, T. P. Devereaux, D. -H. Lee, and Z.-X. Shen, Nature 515, 245 (2014).
- (13) M. D. Johannes and I. I. Mazin, Phys. Rev. B 77, 165135 (2008); H.-M. Eiter, M. Lavagnini, R. Hackl, E.A. Nowadnick, A.F. Kemper, T.P. Devereaux, J.-H. Chu, J. G. Analytis, I. R. Fisher, and L. Degiorgi, Proc. Nat. Acad. Sciences 110, 64 (2013).
- (14) M. Opel, R. Hackl, T. P. Devereaux, A. Virosztek, A. Zawadowski, A. Erb, E.Walker, H. Berger, and L. Forró, Phys. Rev. B 60, 9836 (1999).
- (15) T. P. Devereaux and R. Hackl, Rev. Mod. Phys. 79, 175 (2007); T. P. Devereaux and A. P. Kampf, Int. J. Mod. Phys. B 11, 2093 (1997).
- (16) T. P. Devereaux, A. Virosztek, and A. Zawadowski, Phys. Rev. B 59, 14618 (1999); ibid 51, 505 (1995).
- (17) A. Damascelli, Z. Hussain, and Z.-X. Shen, Rev. Mod. Phys. 75, 473 (2003); L. Pintschovius, Phys. Stat. Sol. (b) 242, No. 1, 30 (2005).
- (18) T. Cuk, D. H. Lu, X. J. Zhou, Z.-X. Shen, T. P. Devereaux, and N. Nagaosa, Phys. Stat. Sol. (b), 242, 11 (2005).
- (19) M. Maschek, S. Rosenkranz, R. Heid, A. H. Said, P. Giraldo-Gallo, I. R. Fisher, and F. Weber, Phys. Rev. B 91, 235146 (2015).
- (20) D. Reznik, L. Pintschovius, M. Ito, S. Iikubo, M. Sato, H. Goka, M. Fujita, K. Yamada, G.D. Gu, and J.M. Tranquada, Nature 440, 1170 (2006); S. R. Park, T. Fukuda, A. Hamann, D. Lamago, L. Pintschovius, M. Fujita, K. Yamada, and D. Reznik, Phys. Rev. B 89, 020506 (R) (2014); M. Le Tacon, A. Bosak, S. M. Souliou, G. Dellea, T. Loew, R. Heid, K.-P. Bohnen, G. Ghiringhelli, M. Krisch, B. Keimer, Nat. Phys. 10, 52 (2014); C. J. Bonnoit, D. R. Gardner, R. Chisnell, A. H. Said, Y. Okada, T. Kondo, T. Takeuchi, H. Ikuta, D. E. Moncton, Y. S. Lee, arXiv:1202.4994.
- (21) M. Braden, L. Pintschovius, T. Uefuji, and K. Yamada, Phys. Rev. B 72, 184517 (2005); J. Graf, M. dÕAstuto, C. Jozwiak, D. R. Garcia, N. L. Saini, M. Krisch, K. Ikeuchi, A. Q. R. Baron, H. Eisaki, and A. Lanzara, Phys. Rev. Lett. 100, 227002 (2008).
- (22) S. Johnston, I. M. Vishik, W. S. Lee, F. Schmitt, S. Uchida, K. Fujita, S. Ishida, N. Nagaosa, Z. X. Shen, and T. P. Devereaux, Phys. Rev. Lett. 108, 166404 (2012).
- (23) S. Moser, S. Fatale, P. Krüger, H. Berger, P. Bugnon, A. Magrez, H. Niwa, J. Miyawaki, Y. Harada, and Grioni, Phys. Rev. Lett. 115, 096404 (2015); Y. Y. Peng et al., Phys. Rev. B 92 064517 (2015).
- (24) L. J. P. Ament, M. van Veenendaal, and J. van den Brink, Euro. Phys. Lett. 95, 2 (2011).
- (25) Y. Wang, B. Moritz, C.-C. Chen, C. Jia, M. van Veenendaal, and T. P. Devereaux, Phys. Rev. Lett. 116, 086401 (2016).
- (26) T. P. Devereaux, D. Einzel, B. Stadlober, and R. Hackl, Phys. Rev. Lett. 72, 396 (1994); ibid., 3291 (1994).
- (27) T. P. Devereaux, A. Virosztek, and A. Zawadowski, Phys. Rev. B 54, 12523 (1996).
- (28) P. Marra, K. Wohlfeld, and J. van den Brink, Phys. Rev. Lett. 109, 117401 (2012).
- (29) C. Jia, K. Wohlfeld, Y. Wang, B. Moritz, T. P. Devereaux, ArXiv:1510.05068, to appear in Phys. Rev. X.
- (30) G. Mahan, Many-Particle Physics, Springer 3rd ed. (2000).
- (31) F. M. F. de Groot, P. Kuiper, and G. A. Sawatzky, Phys. Rev. B 57, 14584 (1998).
- (32) M. van Veenendaal, Phys. Rev. Lett. 96, 117404 (2006).
- (33) Luuk J. P. Ament, Giacomo Ghiringhelli, Marco Moretti Sala, Lucio Braicovich, and Jeroen van den Brink, Phys. Rev. Lett. 103, 117003 (2009).
- (34) M. W. Haverkort, Phys. Rev. Lett. 105, 167404 (2010).
- (35) M. Minola et al., Phys. Rev. Lett. 114, 217003 (2015); L. Braicovich et al., Rev. Sci. Inst. 85, 115104 (2014).
- (36) http://www.esrf.eu/UsersAndScience/Experiments/EMD/ID32