Electron-momentum dependence of electron-phonon coupling underlies dramatic phonon renormalization in YNi2B2C.

Electron-phonon coupling, i.e., the scattering of lattice vibrations by electrons and vice versa, is ubiquitous in solids and can lead to emergent ground states such as superconductivity and charge-density wave order. A broad spectral phonon line shape is often interpreted as a marker of strong electron-phonon coupling associated with Fermi surface nesting, i.e., parallel sections of the Fermi surface connected by the phonon momentum. Alternatively broad phonons are known to arise from strong atomic lattice anharmonicity. Here, we show that strong phonon broadening can occur in the absence of both Fermi surface nesting and lattice anharmonicity, if electron-phonon coupling is strongly enhanced for specific values of electron-momentum, k. We use inelastic neutron scattering, soft x-ray angle-resolved photoemission spectroscopy measurements and ab-initio lattice dynamical and electronic band structure calculations to demonstrate this scenario in the highly anisotropic tetragonal electron-phonon superconductor YNi2B2C. This new scenario likely applies to a wide range of compounds.

Introduction

Interacting degrees of freedom in solids underlie new emergent ground states and competing phases with potential for new functionalities. Vibrations of the atomic lattice, i.e. phonons, can couple to electrons[1], magnetic[2], or orbital degrees of freedom[3]. In particular, electron–phonon coupling (EPC) received a lot of attention as a microscopically understood origin of superconductivity. Furthermore, EPC has recently been in the focus of investigations of materials with competing phases, such as cuprates[4-8] and layered transition-metal dichalcogenides[9-13].

Reduced phonon lifetimes are typically related to nesting where large sections of the Fermi surface (FS) are connected by a single phonon wave vector[14-16]. The hallmark of nesting is a singularity or at least a strong peak in the electronic susceptibility for a particular phonon wave vector q connecting parallel sections of the FS[17,18].

Such phonon anomalies received a lot of attention recently in the context of charge density wave (CDW) formation, which is often driven by a soft phonon mode triggering a structural distortion at the CDW transition. In some CDW compounds, such as 1D conductors, this behavior is indeed related to the nesting of the FS[19-21]. Yet in others, such as 2H-NbSe2, nesting is absent[22-26] and cannot explain the soft-mode properties. While the CDW in 2H-NbSe2 originates from EPC[22,27], it has been proposed that the phonon softening and broadening on cooling towards the CDW transition temperature TCDW can be explained only by taking into account lattice fluctuations[25], and anharmonic effects may play an important role[24]. Phonon anomalies in cuprates remain enigmatic though a strong response to the onset of superconductivity is evident[4,7]. In other cases, an expected phonon broadening related to nesting is absent[28] or incomplete[29]. In principle, the phonon momentum, q, dependence of the EPC itself, expressed in the EPC matrix elements , can determine the wavevector of the phonon broadening in the absence of FS nesting[26,27]. Yet, changes of the phonon linewidth because of the temperature dependence of up to room temperature are expected to be small.

Here, we propose a scenario in which the electron momentum, k, dependence of the EPC matrix elements comes into play. In such a scenario, the , is particularly large for electrons on certain parts of the FS. We show that such k-selective EPC can be strong enough to significantly broaden phonons even when k-integrated quantities like the electronic susceptibility χ lack particular features at the phonon momentum q. In this case, the broadening can sensitively depend on the temperature-induced changes of electronic states at the Fermi energy. This scenario can explain large temperature-dependent phonon line widths in the absence of both nesting and anharmonicity.

We demonstrate this scenario on the electron–phonon superconductor YNi2B2C (Tc = 15.2 K). YNi2B2C is known for unusual low-temperature phonon lineshapes reflecting the anisotropic superconducting energy gap[30-32], EPC distributed over 160 meV in phonon energy[33] and phonons with unusual eigenvectors mediating superconductivity[34]. This work focuses on the microscopic origin of strongly increased line widths of certain phonons observed over a large range of wavevectors in the normal state at low temperatures, i.e. T = 20 K.

We first argue that previously proposed explanations for the phonon properties in YNi2B2C, i.e., FS nesting[35], can be ruled out. Then we demonstrate that results of our comprehensive inelastic neutron scattering (INS) and soft x-ray angle-resolved photoemission spectroscopy (SX-ARPES) measurements agree well with ab-initio lattice dynamical and electronic band structure calculations, which allows us to use these calculations to gain insights into the microscopic origin of EPC. The calculations highlight the importance of 2D electronic joint density of states (2D-eJDOS). It is defined as the usual electronic joint density of states (eJDOS) but evaluated for 2D slices of the reciprocal space at a specific component of qk, where the z-direction is defined to be parallel to the crystallographic c-axis (Fig. 1a, b). The results highlight the decisive role of the interplay between the k-dependence of the EPC combined with a strongly k-dependent 2D-eJDOS, and explain strongly temperature-dependent phonon broadening in the absence of FS nesting and lattice anharmonicity.

Fig. 1

Electron–phonon coupling.

a Crystal structure of YNi2B2C (, , ) and Feynman diagram of the EPC process including the EPC matrix element . Wave vectors and refer to electron and phonon momentum space, respectively. b Calculated partial FS of YNi2B2C. Colors denote different bands crossing the Fermi energy. c Calculated phonon dispersion featuring an anomaly at . Note nearly parallel sections of the Fermi surface at kz = 0.5 where 2D-eJDOS is enhanced.

Electron–phonon coupling renormalizes the phonon energy and reduces the phonon lifetime through emission or absorption of electron–hole pairs. Quantitatively, phonon energy is affected by the real part of the dependent electronic susceptibility χ, whereas the lifetime is determined by the imaginary part (see Feynman diagram in Fig. 1a). The former depends on electronic states near as well as far from the Fermi surface, whereas only electronic states near the FS contribute to the latter.

Thus knowledge of both electrons and phonons as well as the coupling between them is necessary to understand phonon self-energy in metals. Electronic states near the FS can be accurately measured by angle-resolved photoemission (ARPES) and the phonon linewidths can be measured by neutron scattering. On the other hand, measurements of electronic states far from the FS that need to be included to obtain phonon peak shift due to interaction with electrons are a lot more challenging especially for the states that are unoccupied. Thus it is easier to compare theory with experiments for the phonon lifetimes than for the peak positions.

Phonon lifetimes are reduced when phonons are absorbed or emitted by electronic excitations across the Fermi surface. Since the energy and momentum are conserved, only occupied electronic states in a small energy window can contribute to this process, where is the energy of a phonon mode with wave vector q in the dispersion branch λ (Fig. 1b). Reduced lifetimes translate into large phonon spectral line widths.

Two quantities are relevant for determining phonon line broadening due to the EPC: (1) The likelihood that a particular phonon with is scattered/absorbed by a particular electron in a state with wave vector k and energy , which is expressed as the EPC matrix element . (2) The imaginary part of the electronic susceptibility reflects the number of electronic states at the FS which are connected by q, i.e. the eJDOS, which is defined as and equivalent to the so-called nesting function [36]. (In the following, we will use the term eJDOS. The equivalent term nesting function could cause confusion since conventional FS nesting plays absolutely no role in the reported effects.) Both and χ are in general q dependent and, therefore, their interplay in q-space determines the q-dependence of phonon linewidths (Fig. 1c). Eq. (1) describes the contribution to the phonon linewidth because of EPC:

Apart from their intrinsic q dependence, the EPC matrix elements can amplify or suppress contributions to the phonon line width for electronic states at different k. This k-selectivity has not yet been considered in the analysis of EPC. We show that it can have a profound impact on the lattice dynamical properties as exemplified by the behavior of the phonon broadening in YNi2B2C.

Results

Lattice dynamics—previous results

YNi2B2C offers unique insights into the interplay between electronic and lattice degrees of freedom via superconductivity-induced phonon anomalies[30,32,37]. Phase competition of superconductivity with CDW order[35], the role of FS nesting[38-40] as well as its lattice dynamics have been investigated extensively[30,41] including in our earlier work[32-34]. We first demonstrate that temperature and dependence of phonon renormalization in YNi2B2C is inconsistent with either standard mechanism: FS nesting or anharmonicity.

Overall, ab-initio lattice dynamical calculations of YNi2B2C (I4/mmm, a = b = 3.51 Å, c = 10.53 Å) agree with phonon spectroscopy results at energies up to 160 meV with regard to both phonon energies and corresponding phonon line widths due to the EPC[32-34]. The strongest phonon anomalies are observed in c-axis polarized transverse acoustic (TA) phonon modes[34]. One such mode at , the point of the Brillouin zone (BZ), is measured at , i.e., in the BZ adjacent to (Fig. 2a). Another strong coupling TA mode is observed at , i.e., about halfway between the and points of the BZ. For simplicity, we call it here the anomaly. We measured it at , i.e. close to the same BZ center as the point phonon, i.e. (Fig. 2b). The anomaly is the only part of the phonon dispersion, which is not very well captured by ab-initio calculations in that the observed anomaly is stronger than predicted[34].

Fig. 2

Phonon broadening in YNi2B2C.

a, b Raw data from inelastic neutron scattering (INS) revealing the strong softening and broadening of c-axis polarized acoustic phonons at a , b , c and d corresponding to the reduced phonon wave vectors , , and , respectively, within the Brillouin zone adjacent to . Horizontal bars indicate the full-width at half-maximum of the peak observed at . Error bars represent s.d. e Calculated total electronic joint density-of-states (eJDOS) along (solid blue line, bottom axis) and (dashed green line, top axis) in space. Circles indicate the phonon wave vectors corresponding to those shown in (a) and (b).

Key enigmatic results are (1) that strong phonon broadening appears over a wide range of phonon wave vectors , , i.e. , but quickly vanishes going away from along the [110] direction and (2) that broadening of the TA mode at is much stronger than at the point. Both TA phonons display a pronounced broadening upon cooling from room temperature to easily identified in high-energy resolution INS data (see the “Methods” section) (Fig. 2a–d). The new data allows for a detailed study of the momentum dependence of phonon broadening (see the “Results” section and Figs. 5, S7–S9).

Fig. 5

Momentum-dependent phonon broadening (DFPT + INS).

a Illustration of the scattering geometry in the horizontal [110]−[001] scattering plane. Measurements were also performed at out-of-plane wave vectors offset along the vertical direction by (red dashed plane) and (blue dashed plane). The , and points of the Brillouin zone are indicated (dots). b–d Color-coded plot of the calculated electronic contribution to the phonon line width (see color bars for corresponding values in units of [meV]) for the transverse acoustic phonon mode investigated close to . Results are shown in the (H,H,L) scattering plane for (b) , (c) , and (d) . Symbols refer to results from INS measurements and indicate phonon renormalization as indicated in the legend of panel (d). (zero structure factor) indicates that the phonon intensity was too weak to be observed—in agreement with DFPT structure factor calculations. Size of symbols (dots, circles) scale with the strength of the phonon renormalization within each panel. The white arrow in (c) indicates corresponding to a reduced wave vector .

The observed broadening at is much stronger than at the point (Fig. 2a, b). In contrast, the calculated eJDOS[34] along the [100] and [110] directions displays a peak only for (Fig. 2e). Yet, the calculations correctly predict the broad momentum range along the [100] direction over which phonon broadening is observed[34]. Therefore, nesting is unlikely to be the origin of the observed strong phonon renormalization. On the other hand, EPC matrix elements [see Eq. (1)] can underlie a dependence not reflected by the FS geometry. However, the matrix elements are expected to change only for temperatures on the order of an electron-volt corresponding to much higher values than those in Fig. 2.

Anharmonic broadening of phonons is typically attributed to phonon–phonon scattering which is typically strong at elevated temperatures.[42,43] At low temperatures anharmonic broadening can occur near or at structural phase transitions as it is argued for the case of the CDW transition in 2H-NbSe2 ().[24] In fact, it was argued that YNi2B2C is close to a structural phase transition related to the soft mode at only precluded by the onset of superconductivity.[35] Hence, the large phonon linewidth just above could still be anharmonic. However, previously reported superconductivity-induced changes of the line shape of the anomalous TA phonons to rule out a non-electronic origin.[30] This redistribution of phonon spectral weight[32] directly reflects the opening of the superconducting gap 2Δ (for more details see Fig. S1 and Supplementary Note 1). A model[37] that accurately describes these observations is purely EPC-based, which is possible only if phonon broadening in YNi2B2C occurs only because of electronic scattering, i.e., is EPC in nature.

In the following, we present evidence that the strong TA phonon broadening in YNi2B2C originates from EPC that is strongly enhanced on parts of the Fermi surface, which are connected by the correct values of the phonon momentum. This scenario explains the results from phonon spectroscopy (see Fig. 2a, b) in the absence of FS nesting or lattice anharmonicity.

Band structure from electron spectroscopy

There has been considerable discussion about FS nesting for in (Lu/Y)Ni2B2C including theoretical[35,38,44] and experimental work[39,40,45,46]. Yet, the clearest signature of electron–phonon coupling, the -induced phonon effect, appears over a large momentum region along the [100] direction[32], which is inconsistent with nesting. In order to clarify the origin of this EPC, we performed SX-ARPES measurements on samples cut from the large YNi2B2C single crystal, which we used for our current and also previous INS studies[32-34].

SX-ARPES measurements at photon energies have a resolution, which is sharp enough to reliably resolve the -dispersion effects in YNi2B2C. (The half out-of-plane periodicity of the BZ, i.e. the distance [see Fig. 3f], is [], which corresponds to a difference in the incident photon energy of ). We compared the experimentally observed band structure with density functional theory (DFT) calculations carried out in the framework of the mixed basis pseudopotential method[47] using the local density approximation (LDA) (for details see ref. [34]).

Fig. 3

Fermi surface from soft x-ray angle-resolved photoemission spectroscopy.

Comparison of calculated Fermi surface (lines) and SX-ARPES intensities (white/black: high/low intensities) observed in the Г − X/Z − M plane of the Brillouin zone (BZ) with (a) k = 23 r.l.u. = 13.7 Å−1 achieved at , (b) at and (c) at . d, e Calculated Fermi surface (lines) and SX-ARPES intensities for the plane at Measurements shown in a–d were performed with p polarization except for the lower right quadrant in (a) which was obtained with s polarized light as data in (e). f Sketch of the BZ of YNi2B2C. Note that X and Z points are equivalent. Wave vectors are given in reciprocal lattice units (r.l.u.) of .

Selected cuts of the FS in the plane in Fig. 3a–c and the plane in Fig. 3d, e overlaid with corresponding DFT calculations (solid blue lines), show a good agreement (see also Figs. S2 and S3 in SI). We focus on the vicinity of the point, i.e. at since the DFT predicts a square-shaped FS around it for certain values prone to nesting[34]. Our results demonstrate the significant evolution of the FS at : An elliptical FS stretches between the four visible points at (Fig. 3a, ), a nearly square-shaped feature is found at (Fig. 3b, ), and a large nearly circular ellipse is observed at (Fig. 3c, ) with the long axis rotated by 90° relative to the orientation at . Note that due to the particular symmetry of the BZ, the positions of and points are reversed upon moving from to (see Fig. 3f).

The experimentally observed spectral weight with a square-like shape for is relatively diffuse, in contrast to the sharp FS expected from the calculated Fermi contour (Fig. 3b). This originates from two features of the band structure in the vicinity of the M point: (1) The dispersion of the inner band is very strong (Figs. 3d, e and S4c) and (2) both the inner and outer bands forming the square have a shallow in-plane energy dispersion of <0.5 eV near (see Fig. S3c). Hence, due to the uncertainty of the inner potential V0 used in rendering photon energy into and to the limited and energy resolution of the experiment, we pick up intensity inside the square.

We used both linear vertical, i.e. , and linear horizontal, i.e. polarization of the light in order to exploit photoemission matrix element effects. polarization is sensitive to the inner (outer) band for odd (even) values of , whereas polarization reveals inner and outer bands for both odd and even . This is evident in Fig. 3a, c where the spectral weight is located at the center and at the edges of the ellipse at the point, respectively, with , but uniformly distributed with (bottom right in Fig. 3a). This is also visible on the dispersion for r.l.u. (Fig. 3d, e). With polarization (Fig. 3d), the experimental periodicity of the band structure looks like it is doubled in comparison to the calculations due to these matrix elements effects. In fact, we recover the correct periodicity using polarization of the light (Fig. 3e). Note that the data of Fig. 3d cut the point in the plane between the first and second BZ, while the data of Fig. 3e cut the point between the second and third BZ. However, Fig. 3a–c indicate that variations of the photoemission matrix elements across different BZs are smaller than variations due to light polarization.

Electronic joint density of states from density functional theory

The SX-ARPES study detailed above shows that the band structure of YNi2B2C is well described by our DFT calculations. Therefore, we rely on the analysis of the calculated FS (see Fig. S4) to assess the nesting properties of YNi2B2C via the calculated eJDOS. The eJDOS along the [100] and [110] directions in (Figs. 2c and 4a) lacks a peak at the position of the most pronounced phonon broadening. However, our calculations reveal a peculiar momentum dependence when we look at the joint electronic density of states analyzed for fixed values of , which we call a 2D-eJDOS. Detailed -dependent calculations for phonon wave vectors in the entire plane reveal that the 2D-eJDOS is low for (Fig. 4b) but increases strongly to (Figs. 4c and S5 and S6).

Fig. 4

k-dependence of electronic joint density-of-states (eJDOS).

a eJDOS in the plane with averaged over the full 3D FS. b, c 2D-eJDOS() considering only the FS states with (b) and (c) . All results are given in states/eV² and the same color-code applies to (a–c). d–f Detailed dependences of the -averaged 2D electronic joint density of states averaged over three different ranges in : (d) full plane [large blue square in (b)], (e) [red bar in (b)] and (f) [green bar in (b)], respectively. The black dashed horizontal line in (d) denotes the -averaged value of deduced from the data shown in (a).

We report the dependence of the 2D-eJDOS in more detail for three ranges defined in Fig. 4b: (1) with [large blue square], (2) [red bar] and (3) [green bar]. Starting with range (1), the corresponding -averaged 2D-eJDOS, , has a clear peak at (Fig. 4d). Corresponding data for ranges (2) (Fig. 4e) and (3) (Fig. 4f) show that the peak of the at is even more pronounced at the wave vectors of the strong phonon broadening. Overall, our analysis demonstrates that the 2D-eJDOS is extremely sensitive to and depends less on the phonon wave vector .

Momentum-dependent phonon broadening from inelastic neutron scattering

So far, we focused on the mechanism explaining the strong broadening of the two TA modes at and (Fig. 2a, b). While these effects are not sharply localized in space, they do depend on the phonon wave vector[34] and only a rigorous experimental test enables an overall verification of the validity of the DFT calculations. Therefore, we performed neutron scattering experiments with improved energy resolution at phonon wave vectors along as well as off the high-symmetry directions investigated previously[30,32,34]. We compare our results to ab-initio calculations of (see Eq. (1)), which were calculated using the linear response technique or density functional perturbation theory (DFPT)[48] in combination with the mixed-basis pseudopotential method[49].

For the present study, we measured phonons at wave vectors close to and the point at room temperature and and looked for peak broadening. We observed broad but clear momentum dependences, e.g., in measurements going away from the point (see Fig. 2a, c, d): Already at close-by wavevectors, the broadening is clearly reduced (Fig. 2c) and completely absent further away in momentum space (Fig. 2d). More INS data taken along different directions for both anomalous modes are shown in Fig. S7 of SI. Figure 5 shows a comparison of the observed (symbols) and calculated (color code) phonon broadening in absolute wave vectors to highlight the clear momentum dependences in 3D (see also Fig. S8 and Supplementary Note 2). Overall, we find good agreement.

Finally, we note that the -averaged eJDOS is about 50% higher for the point than for (Figs. 2e and 4a), whereas the (Figs. 4c, e, f and  S6) and (Fig. 5b, c) have similar values at these two wave vectors. This is further evidence that the electronic states at the Fermi level with are selected by EPC matrix elements and mediate the coupling responsible for the broadening of the TA phonons in YNi2B2C.

Discussion

EPC has been studied for a long time because it can stabilize emergent ground states as well as lead to phase competition. The particular role/impact of EPC matrix elements was discussed already in the 1970s[50,51], though experimental evidence for the decisive role of the momentum dependence of the EPC matrix elements was only reported in the last decade[24,26,27,52-54]. Still, only the dependence of the EPC matrix element was scrutinized. Our work takes the full and momentum dependence of into account. We can explain not only phonon broadening due to EPC as a function of but also its pronounced weakening at increased temperatures in the absence of FS nesting and anharmonicity. The peak of the as function of (Fig. 4d–f) accounts for a high sensitivity to temperature via the electronic band structure if the contribution of the electronic states with to EPC is boosted by -selective matrix elements. The peak of the (k) (Fig. 4d–f) is superficially reminiscent of classic FS nesting scenario where the peak in the -averaged eJDOS() associated with the nesting condition is responsible for a phonon anomaly. At increased temperatures, smearing of the Fermi function quickly suppresses this peak in the eJDOS(q) and reduces the phonon broadening[18]. In our case, the peaks in the (k) (Fig. 4d–f) are expected to also be strongly broadened and reduced in height with increasing temperatures. This is exactly what is observed in experiments. Our results demonstrate that selectivity of the EPC can boost the contribution of only a small part of the Fermi surface to the phonon lifetime. This new point of view has to be taken into account generally when assessing phonon anomalies in metallic systems.

This finding makes it necessary to reconsider our understanding of some well-known materials. NbSe2 is a model system for various fundamental aspects (phase competition as function of dimensionality[12]; quantum phase transitions[55,56]; 2D quantum metal state[57]; Higgs mode in condensed matter physics[58,59]). The correct view of EPC is the indispensable ingredient for explaining many of the observed effects. In particular, NbSe2 is considered to be a showcase for -dependent EPC matrix elements featuring concomitant charge-density-wave (CDW) order () and phonon-mediated superconductivity ()[60]. Calculations[25] and experiments[26] have shown that FS nesting is absent and concluded that the periodicity of the CDW is determined by the momentum dependence of the EPC matrix elements. Thus, the phonon broadening and softening close to is due to the EPC. However, lattice fluctuations overwhelm EPC away from and explain the strongly reduced phonon renormalization at elevated temperatures[24,25]. We find that the calculated eJDOS across the CDW ordering wavevector in 2H-NbSe2 is nearly featureless whereas the calculated line width due to EPC, , shows a pronounced broad maximum around in agreement with experiment[27]. The situation is, in fact, very similar to that observed for the anomaly in YNi2B2C (see Figs. 2e and S8a). In analogy to the results presented here, the strong temperature dependence of the phonon line width in 2H-NbSe2 could also originate from an interplay of -selective EPC and the electronic band structure. In fact, Flicker et al. [25] already considered an orbital-dependent EPC matrix element. Yet, it remains unclear whether the model calculations also support a strongly dependent 2D-eJDOS.

Recently, phonon anomalies related to CDW order in cuprates have attracted large scientific interest. In particular, a strong phonon broadening was observed at and below the onset of CDW fluctuations competing with superconductivity[4,6,8]. The abrupt decrease of the phonon line width in YBa2Cu3O6.6 on entering the superconducting state[4] indicates a high-sensitivity to the opening of the superconducting energy gap on the FS mediated by the EPC. Yet, the mechanism determining the periodicity of the CDW remains unclear. The reported sharp dependence of the phonon broadening is reminiscent of a FS nesting-type origin but the nesting could not be identified in the electronic band structure. selective EPC in concert with large 2D-eJDOS on a small part of the FS is a novel approach for an improved understanding of these observations.

Methods

Inelastic neutron scattering

Experiments were performed at the 1T and IN8[61] triple-axis spectrometers located at Laboratoire Léon Brillouin (LLB), CEA Saclay, and Institute Laue-Langevin (ILL), Grenoble, respectively. Double focusing pyrolytic graphite monochromators and analyzers were employed on the 1T spectrometer. We used double-focusing copper (Cu200) monochromators and analyzers on IN8. A fixed analyzer energy of 14.7 meV allowed us to use a graphite filter in the scattered beam to suppress higher orders. The phonon scattering wave vector is expressed in reciprocal lattice units (r.l.u.) with the lattice constants and of the tetragonal unit cell. The single-crystal sample was mounted in a standard orange cryostat at ILL and in a closed-cycle refrigerator at LLB, allowing measurements down to T = 2 K.

Soft x-ray angle-resolved photoemission spectroscopy

experiments were performed on in-situ cleaved single crystals of YNi2B2C (001 surface) at the SX-ARPES endstation[62] at ADRESS beamline[63] of the Swiss Light Source (SLS) with incident photon energies in the range (for more details see[62]). The increase of the photoelectron mean-free path related to the increased kinetic energy in the soft x-ray regime translates into a proportional improvement of the intrinsic resolution of the photoemission experiment, which is indispensable for our investigation. The measurements were performed at low sample temperature in order to suppress relaxation of -selectivity because of thermal motion[64]. For consistency with the results from neutron scattering, the electron scattering wave vector is expressed in r.l.u. .

Density functional theory

calculations were carried out in the framework of the mixed basis pseudopotential method[47] using the local density approximation (LDA) (for details see ref. [34]) with a grid of points in space. Plots of the electronic band structure and its properties are based on an interpolated grid to enhance the in-plane accuracy. All results presented here were obtained for the experimental lattice constants of the tetragonal structure (, ) with internally relaxed atomic coordinates. Phonon properties and EPC matrices were calculated for a grid in phonon, i.e. momentum space using the linear response technique or density functional perturbation theory (DFPT)[48] in combination with the mixed-basis pseudopotential method[49]. While the eJDOS is defined as an integral over the 3D Brillouin zone (see main text), the 2D-eJDOS is evaluated as a 2D integral for a fixed value of k, i.e., it samples the electronic phase-space over slices parallel to the k − k plane. δ-functions are approximated by Gaussians of width σ = 0.05 eV.