An electronic solid with itinerant carriers and localized magnetic moments represents a paradigmatic strongly correlated system. The electrical transport properties associated with the itinerant carriers, as they scatter off these local moments, have been scrutinized across a number of materials. Here, we analyze the transport characteristics associated with ultraclean PdCrO 2 -a quasi-two-dimensional material consisting of alternating layers of itinerant Pd-electrons and Mott-insulating CrO 2 layers-which shows a pronounced regime of T -linear resistivity over a wide range of intermediate temperatures. By contrasting these observations to the transport properties in a closely related material PdCoO 2 , where the CoO 2 layers are band-insulators, we can rule out the traditional electron-phonon interactions as being responsible for this interesting regime. We propose a previously ignored electron-magneto-elastic interaction between the Pd-electrons, the Cr local moments and an out-of-plane phonon as the main scattering mechanism that leads to the significant enhancement of resistivity and a T -linear regime in PdCrO 2 at temperatures far in excess of the magnetic ordering temperature. We suggest a number of future experiments to confirm this picture in PdCrO 2 as well as other layered metallic/Mott-insulating materials. Recent years have witnessed a resurgence of interest in the microscopic origin of an electrical resistivity that scales linearly with temperature (1,2) and exhibiting a Planckian scattering rate, Γ = Ck B T /ℏ, where C ∼ O( 1) coefficient (3)(4)(5)(6)(7). In conventional (simple) metals at room temperature, this phenomenology is readily understood as a consequence of electrons scattering off thermally excited phonons in their equipartition regime (8). On the other hand, numerous "correlated" materials belonging to the cuprate (6,7,9,10), pnictide (11,12), ruthenates (13)(14)(15)(16)(17), rare-earth (18), and moiré bilayers (5,19) display Planckian scattering down to low temperatures, likely driven by purely electronic interactions and where a priori it is unclear whether phonons play an essential role (3,(20)(21)(22). It is challenging to disentangle the role of electron-electron and electronphonon interactions on scattering lifetimes. It is quite natural to ask whether materials with a nearly identical phonon spectrum and distinct electronic spectra can lead to a distinct temperature dependence of their respective resistivities.

The goal of this paper is to resolve a conundrum inspired by electrical transport measurements in two isostructural quasi-two-dimensional compounds with distinct electronic structures: PdCoO 2 and PdCrO 2 . Their structural motif consists of alternately stacked layers of highly conducting Pd and insulating CoO 2 /CrO 2 in a triangular lattice arrangement (24-28); see Fig. 1A. The phonon spectra for the two compounds are nearly identical; some of the differences arise from the distinct ionic masses (27)(28)(29)(30), and recent analysis has also revealed that the unit cell of PdCrO 2 is slightly enlarged due to the presence of magnetic moments (31). On the other hand, their electronic spectra are different since the CrO 2 layers are Mott-insulating with the local moments interacting via antiferromagnetic (AFM) exchange interactions (25,32), while the CoO 2 layers are nonmagnetic (24,26,33). The photoemission spectrum of PdCrO 2 contains prominent features, absent in PdCoO 2 (28,(34)(35)(36), that can be understood from an effective interlayer Kondo lattice model (28). The in-plane resistivities for the two compounds (23) are shown in Fig. 1B, respectively. The salient features are as follows: i) The magnitude of the resistivity for both compounds is small, suggesting that they are "good" metals with a long mean-free path. ii) PdCrO 2 is considerably more resistive than PdCoO of the resistivity above T ≳ 150 K (far above T N ≈ 37.5 K, the Néel temperature for 120 • -antiferromagnetism) (25,32,37,38), and with a slope that is greater than the average slope of ab (T ) in PdCoO 2 in the same temperature range (39).

The central puzzle that we address in this paper concerns the microscopic origin of the excess T -linear resistivity in PdCrO 2 relative to isostructural PdCoO 2 (Δ ≡ PdCrOfoot_1 ab -PdCoO 2 ab > 0), going beyond the conventional electron-phonon scattering mechanism, and in a temperature regime where the long-range magnetic order is lost. Given the contrast between PdCrO 2 and PdCoO 2 , it is plausible that the fluctuations of the Crlocal moments play a crucial role in the electronic transport lifetimes even at the relatively high temperatures of interest (i.e., for T ≳ T N ). However, recent work (40) has demonstrated that electrons scattering off the fluctuations of a "cooperative" paramagnet (41)(42)(43)(44)(45)(46)(47)(48) can not account for a T -linear resistivity; instead, the resistivity saturates to a temperature-independent value for T T N . Starting with a microscopic model, we will now demonstrate that the resolution to the conundrum lies in a previously ignored and nontrivial interaction term between the Pd-electrons, the Cr-spins, and phonons, as encoded in an electron-magneto-elastic (EME) coupling. Although specifically motivated by PdCrO 2 , our theory has relevance beyond this single material. There are a number of exciting new material platforms that have come to the forefront in recent years that consist of stacks of metallic and Mott insulating layers (49)(50)(51)(52)(53)(54)(55)(56)(57).

In what follows, we develop a general framework to address electrical transport in such layered material platforms, and our conclusions can be used to disentangle the various sources of interaction between electrons, local moments, and phonon degrees of freedom.

Consider a quasi-two-dimensional (2D) layered model defined on a triangular lattice, where the electronic and local-moment degrees of freedom reside on alternating layers. The effective 2D Hamiltonian is given by (see ref. 39 for a microscopic derivation of H EME ),

Here, p † k , p k denote the Pd-electron creation and annihilation operators with momentum k and spin = ±1/2. The electronic dispersion is given by k , and is the chemical potential; experiments in PdCrO 2 indicate that the Pd-electronic structure is well captured using first-and second-neighbor hoppings and the conduction band, dominantly of Pd character, is very close to half-filling (23,28). The local moments, S i , interact mutually via nearest-neighbor antiferromagnetic Heisenberg exchange (J H > 0), and with the Pd-electron spin-density via a Kondo exchange (J K > 0), respectively. The Cr-electrons form S = 3 2 local moments (28). Finally, we include three phonon fields, ( ) with = I a , I o , O, corresponding to in-plane acoustic (I a ), in-plane optical (I o ), and out-of-plane (O) lattice vibrations, respectively. We set the mass, M , to be equal for these modes for simplicity. The in-plane modes I a,o couple to the p-electron density with strengths a,o , respectively, while the out-of-plane mode, O, couples to the "interlayer" Kondo interaction with EME strength, . We set the lattice constant a = 1, unless stated otherwise. Importantly, we include couplings to both in-plane acoustic and optical modes to account for the full T -dependence of ab in the absence of magnetism (i.e., in PdCoO 2 ) (29). Henceforth, we neglect the weak momentum dependence and the form factors associated with the different interaction terms in the quasi-twodimensional setting to simplify our discussion (28). We will restore these additional complexities when considering out-ofplane transport for reasons to be made clear below.

Let us begin by considering the simpler case where the optical modes are Einstein phonons, with I o ,q = 0 and O,q = 0 , while for the acoustic mode I a ,q = cq, with a corresponding Debye frequency D . The experimental regime of interest corresponds to { 0 , J H } ≲ T F , where F is the Fermi energy for the Pd-electrons. Moreover, we shall consider the limit where J H { ℏ/M 0 , J K } and thereby ignore the feedback of both electrons and phonons on the properties of the local moments. In recent work (40), some of us analyzed the properties of a subset of the terms (H el + H S + H K ) in Eq. 1a at leading order in a small J K by approximating the local moments as O(3) vectors, but capturing their complex precessional dynamics using the Landau-Lifshitz equations (41)(42)(43)(44)(45)(46)(47)(48). While this leads to an interesting frequency dependence and momentum-dependent cross-overs in the electronic self-energy, the temperature dependence can be understood entirely based on a high-temperature expansion with uncorrelated local moments. The present manuscript will treat the local moments on the same footing but include the additional interaction effects due to (H ph + H el-ph + H EME ).

We will analyze electrical transport for the model defined above within the framework of the traditional Landau-Boltzmann paradigm (8). This is justified based on the magnitude of resistivity being much smaller than the characteristic scale of ha B /e 2 (a B ≡ Bohr radius) and k F mfp 1 over the entire temperature range of interest ( mfp being the mean free path) (23). Moreover, there is direct experimental evidence for the value of the dimensionless Kondo-coupling being small, based on recent photoemission experiments (28).

Considering the full H , we have multiple sources of scattering for the electrons. Within Boltzmann theory, the total transport scattering rate satisfies Matthiessen's rule (8):

and the in-plane resistivity is given by ab = m/(ne 2 tr ), where m is the effective mass and n is the electron density (39). Experiments with controlled amounts of irradiation shift the overall resistivity curves by a constant (and relatedly, the residual resistivity), without affecting the slope in the T -linear regime (31). We start by describing the electron-phonon contribution, 1/ el-ph ≡ 1/ el-ph,a +1/ el-ph,o . Previous works have obtained 1/ el-ph for PdCoO 2 and highlighted the importance of a highfrequency optical mode (which is not entirely in the equipartition regime at T ≲ 0 ) for the observed superlinear-scaling of ab (T ) (29,58). We have fully reproduced these results for PdCoO 2 based on the same procedure (39) within our 2D model of the Fermi surface; see Fig. 2. In PdCrO 2 , the scattering of electrons off the local moments due to the bare Kondo interaction, 1/ K , can lead to a sublinear T -dependent contribution to ab at temperatures T ≳ J H , before saturating to the T -independent value (40) (see K = m/(ne 2 K ) in the Inset of Fig. 2). We now turn to the important role of the EME term in PdCrO 2 . For 0 ≳ J H , we find that -1 EME follows closely the temperature dependence of the scattering rate of electrons interacting with an optical mode of frequency 0 with a modified dimensionless coupling, Here, 0 is the electronic density of states at the Fermi level, and scattering off the (spin-S) local-moment fluctuations via the EME interaction leads to the additional factor of [S(S + 1)] (39). Ignoring the constant offset, 1/ K at T J H , we find that

where the dimensionless coefficients are given by el-ph,o =

Consequently, at the highest temperatures, the effect of the EME term considered in this work is to enhance the slope (A) of a T -linear resistivity, ab -0 = AT . This constitutes our first important result. Importantly, even if the bare EME coupling is weak relative to the electron-acoustic phonon coupling (i.e., / a 1), the dimensionless coupling EME is not necessarily small compared to el-ph,a . Furthermore, if the out-of-plane phonon is soft ( 0 / D / a ), the presence of the EME interaction can dramatically reduce the onset of T -linear resistivity to O( 0 ).

Assembling all of the above ingredients, we can now reproduce the resistivity in PdCrO 2 , including the effect of the EME term; see Fig. 2. Note that the scattering rates due to the electronphonon interaction in PdCoO 2 (58,59) and Kondo coupling in PdCrO 2 (28) are fixed by previous experiments, which leaves two independent parameters in our theory-EME and 0 [J H is also fixed (28)]. We determine these parameters by fitting the excess resistivity Δ in Fig. 1B to the analytical form of -1 EME (39). The values obtained by this procedure are consistent with the characteristic out-of-plane lattice vibration frequency being naturally softer than the in-plane one, 0 ≳ 0 ≳ J H ; however, we note that our theory extends beyond this regime. The resulting contribution to the resistivity, EME , is shown in the Inset of Fig. 2. Overall, the prominent T -linear resistivity at intermediate T stems from i) the EME scattering rate and ii) the combined sublinear and superlinear contributions of K and el-ph , respectively (39) in their equipartition regime. For PdCrO 2 , this corresponds to T ∼ O(1,000) K due to the large value of 0 and is hence not directly related to the behavior presented in Fig. 2.

Our discussion thus far has focused on the simplified limit of an EME coupling to optical phonons. Let us now analyze the effects of an EME coupling to an acoustic phonon. For H ph in Eq. 1e, this amounts to replacing 2 0

, where q = cq in the limit of small q. Similarly, for H EME in Eq. 1g, this amounts to replacing

i . Note that, in practice, the out-of-plane vibrations that couple the layers are at a finite wavevector, namely, q ≈ c 2 q 2

x + c 2 q 2 y + c 2 z q 2 z,0 in the limit where c z q z,0 T . It is worth noting that unlike the conventional electron-phonon interaction, where scattering is mainly smallangle up to the BG temperature, the EME term induces largeangle scattering even at low T due to large momentum transfer to the local moments, which serve as a "bath." If the spin structure factor exhibits nontrivial correlations in the Brillouin zone (i.e., the spin-correlation length is finite with remnants of Bragg-like peaks), the intermediate-scale transport behavior is controlled by the Pd-electron Fermi-surface geometry. However, when the spin-correlation length is short, 1/ EME shows two distinct regimes (39). For T ≳ T BG ≡ 2 ck F , the Bloch-Grüneisen temperature, the result reduces to the case of optical phonons, 1/ EME = 2 EME T with EME = 0 2 M c 2 S(S + 1). On the other hand, for J H ≲ T ≲ T BG , we encounter an unexpected 1/ EME = 2 EME T 2 / T BG , instead of the usual ∼Tfoot_3 regime in two-dimensions for the phase-space reasons introduced above. Interestingly, this is an example of a T 2 "quasi-elastic" scattering due to the EME term (instead of the usual T 2 due to Umklapp scattering).

Given that the proposed EME interaction in PdCrO 2 originates from fluctuations of the (interlayer) Kondo coupling, applying caxis pressure is expected to enhance the slope of the T -linear resistivity for the following reason. The bare EME coupling is controlled in part by the Kondo scale, ∝ t 2 cp /U , where t cp is the interplane hybridization between the Pd and Crelectrons and is the on-site Coulomb repulsion for Cr-electrons (28,39). Upon applying c-axis strain, the interlayer distance reduces, thereby increasing t cp , which is exponentially sensitive to the deformation; the stiffening of the out-of-plane phonons is at best algebraic. Therefore, the dimensionless EME coupling, EME ∝ t 4 cp /(U 0 ) 2 is expected to show a significant increase, along with an enhancement of the Kondo coupling which affects the constant shift in the resistivity at high T . The predicted form of the in-plane transport is depicted in the Inset of Fig. 3B for a range of c-axis strain ( zz ) for bare microscopic parameters as chosen in Fig. 2.

The electrical resistivity along the c-axis provides a direct window into the interlayer nature of the magnetic interactions, which we have absorbed so far in the effective 2D model. We consider the leading contributions to the c-axis conductivity within linearresponse theory, which is given by

where coh arises from the "coherent" channel due to interlayer p to p hoppings (t pp ), and K , EME represent the "incoherent" channels due to interlayer spin-assisted and spin-phonon-assisted contributions, respectively (39). For simplicity, we ignore the contribution due to an interlayer "incoherent" phonon-assisted hopping not involving the local moments; such a term does not affect our results at a qualitative level. The leading-order Feynman diagrams corresponding to each of these contributions are depicted in Fig. 3A. We have coh = e 2 (n/m) c tr , where (n/m) c is related to the c-axis dispersion and tr is given by Eq. 2; the detailed expressions for the incoherent channels appear in (39). The coherent channel dominates c up to temperatures T ≲ T * ∼ 10 3 K [determined by the condition coh (T * ) ≈ EME (T * ) (39)], such that, in this T -regime, c and ab follow the same T -scaling as they share the same transport lifetime. The contribution of the incoherent channels becomes significant at temperatures T ≳ T * , where, due to the weak c-axis dispersion, their temperature dependence is determined by the T -scaling of the current vertices, rather than the transport lifetime. In particular, for T ≳ D , 0 , while coh ∼ 1/T as in the in-plane case, K is independent of T and EME ∝ T . As a result, the c-axis resistivity becomes sublinear at sufficiently high temperatures, as depicted in Fig. 3B.

The interplay between the different conduction mechanisms has signatures in the behavior under c-axis pressure. The in-plane conductivity is expected to decrease with c-axis compression since the EME scattering is enhanced (because it is proportional to the interplane hopping strength). In contrast, the coherent part of the c-axis conductivity increases, as a result of the enhanced interlayer hopping (39). Interestingly, this increase in coh is associated solely with the el-ph scattering term. The contributions to coh from the Kondo and EME terms are proportional to t 2 pp /t 4 cp (39); assuming that t pp ∝ t 2 cp , the ratio t 2 pp /t 4 cp is unchanged by c-axis strain. Similarly, the incoherent parts of the conductivity increase with compression. However, they do so slightly in excess of coh , which in turn reduces the cross-over scale T * , as manifested by the increase in curvature at intermediate temperatures with increasing compression; see Fig. 3B (39). There are measurements of the c-axis resistivity in the literature (60,61), but the reported values are inconsistent with each other. The reason for this discrepancy is currently unclear. To resolve these issues, more accurate measurements are needed, using, for example, the techniques described in ref. 62.

Our conjectured magneto-elastic mechanism for the enhanced T -linear resistivity in PdCrO 2 relies on quasi-elastic scattering, where the phonons are in the equipartition regime. There is experimental evidence for the Lorenz ratio satisfying the Wiedemann-Franz law in the T -linear (31), which is consistent with our mechanism. Interestingly, the extracted transport scattering rate in the same regime of T -linear resistivity is Planckian with C ≈ 0.9 (31). Within our model and in the quasi-elastic regime, this is not indicative of any fundamental principle, such as a bound associated with an inelastic scattering rate. It is, however, far from obvious why the scattering rate turns out to be Planckian.

A natural future direction is to study the effects of the EME term at low temperature. In particular, the magneto-elastic coupling might be evident in the electronic spectral function. Signatures of phonon drag observed in PdCoO 2 (58) are expected to be suppressed in PdCrO 2 due to EME-induced large-angle scattering of phonons off magnetic moments (25). Pronounced magnetic correlations may also lead to a generalized Kohn anomaly (63) associated with phonon softening at the AFM wavevectors. A detailed understanding of the low-temperature properties of this interesting system remains an open problem.

Data, Materials, and Software Availability. All study data are included in the article and/or SI Appendix. Previously published data were used for this work (31).