Electron-Phonon Coupling (EPC) Overview

• Electron-Phonon Coupling (EPC) is the interaction between electrons and lattice vibrations that governs phenomena like superconductivity, polaron formation, and charge-density waves.
• First-principles methods such as DFPT and spectroscopic techniques like ARPES, Raman, and RIXS quantitatively resolve mode-specific EPC and its impact on material properties.
• EPC insights enable the design of quantum materials with tailored electrical and thermal behavior, while challenges remain in strongly correlated and topologically complex systems.

Electron-Phonon Coupling (EPC) is the fundamental interaction between electronic excitations and lattice vibrations in crystalline, molecular, and low-dimensional systems. This coupling is one of the central paradigms in condensed matter physics, underpinning phenomena such as conventional superconductivity, charge-density-wave formation, electrical and thermal transport, band structure renormalization, polaron formation, and the nonequilibrium modification of materials by optical driving. In systems with strong correlations, spin-orbit coupling, topological bandstructure, or complex moiré superlattices, the precise characterization and control of EPC are essential for identifying emergent collective states. The quantitative and mode- and momentum-resolved determination of EPC remains a forefront challenge, requiring combined theoretical and experimental advances.

1. Formal Definitions and Theoretical Framework

A general Hamiltonian describing electrons, phonons, and their coupling is

$$H = H_{\text{el}} + H_{\text{ph}} + H_{\text{el-ph}}$$

where

• $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$ describes the electronic spectrum,
• $$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$ is the phonon Hamiltonian, and
• $H_{\text{el-ph}}$ encodes the coupling:

$$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$

The EPC matrix element

$$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$

quantifies the amplitude for an electron to scatter from $|n,\mathbf{k}\rangle$ to $|m,\mathbf{k+q}\rangle$ by emitting or absorbing a phonon with branch index $\nu$ and wavevector $\mathbf{q}$; $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$0 is the change in the self-consistent potential under the phonon displacement.

The central quantity characterizing the total strength of EPC at the Fermi surface is the Eliashberg spectral function:

$$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$1

and the dimensionless mass enhancement parameter ("Eliashberg $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$2"):

$$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$3

This $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$4 parameter controls the electron mass renormalization, $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$5, and enters Migdal-Eliashberg theory for phonon-mediated pairing.

For conventional superconductors, the transition temperature $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$6 is given semi-quantitatively by the McMillan-Allen-Dynes equation:

$$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$7

where $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$8 is the logarithmic average phonon frequency, and $$H_{\text{el}} = \sum_{n\mathbf{k}} \epsilon_{n\mathbf{k}}\, c^\dagger_{n\mathbf{k}} c_{n\mathbf{k}}$$9 is the Coulomb pseudopotential (Wang et al., 2024, Zhong et al., 2022).

2. First-Principles Approaches to EPC Calculations

First-principles EPC calculations are commonly performed in the density-functional perturbation theory (DFPT) framework, which enables evaluation of phonon frequencies and $$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$0 matrix elements over the Brillouin zone (Ouyang et al., 5 Nov 2025, Zhan et al., 2024, Wang et al., 2024).

There are two mathematically equivalent but technically distinct electronic-structure approaches:

• Derivative-of-Hamiltonian (dH): Computes $$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$1 as the derivative of the Hamiltonian with respect to phonon normal coordinate, acting on unperturbed wavefunctions.
• Derivative-of-States (d$$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$2): Perturbs the Kohn-Sham states and eigenvalues via first-order perturbation theory and computes $$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$3 from overlaps between perturbed and unperturbed wavefunctions.

Rigorous benchmarking across CP2K (Gaussian basis) and VASP+PAW+Wannier90 (plane-wave/MLWF) implementations demonstrates that the dH approach is numerically robust across codes, whereas the d$$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$4 method is more sensitive to level crossing and eigenstate matching, particularly for low-frequency modes and in cases of near-degeneracy (Merkel et al., 28 Jul 2025). For most molecular and crystalline systems, agreement of $$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$5 values is better than $$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$6 in dimensionless units across codes and methods.

For high-throughput applications and reduction of computational cost, E(3)-equivariant neural networks have been applied to predict the Hamiltonian and its gradients, enabling evaluation of the EPC matrix within milliseconds at DFT accuracy, as demonstrated on H$$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$7O and MoS$$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$8 (Zhong et al., 2023).

3. Experimental Probes and Mode-Resolved Analyses

EPC manifests in a broad array of experimental signatures:

• Angle-Resolved Photoemission Spectroscopy (ARPES): "Kinks" in the quasiparticle dispersion due to EPC yield direct information on the real part of the electronic self-energy, from which the mode-resolved $$H_{\text{ph}} = \sum_{\mathbf{q}\nu} \hbar\omega_{\mathbf{q}\nu}(b^\dagger_{\mathbf{q}\nu} b_{\mathbf{q}\nu} + 1/2)$$9 and $H_{\text{el-ph}}$0 can be extracted via inversion methods (Zhong et al., 2022, Zhu et al., 2013).
• Raman and Infrared Spectroscopies: Phonon frequency shifts, linewidth broadening, and Fano lineshape asymmetry (in the presence of electronic continuum) reveal the impact of EPC on zone-center modes. Mode-resolved coupling constants can be estimated via Allen’s formula:

$H_{\text{el-ph}}$1

where $H_{\text{el-ph}}$2 is the measured phonon linewidth (Zhang et al., 2013).

• Resonant Inelastic X-ray Scattering (RIXS): RIXS quantifies mode- and momentum-resolved $H_{\text{el-ph}}$3 parameters via multi-phonon overtone intensities and detuning analysis, allowing extraction of dimensionless $H_{\text{el-ph}}$4 for specific phonon branches (Braicovich et al., 2019, Peng et al., 2021).
• Two-Dimensional EPC Spectroscopy: Recent developments allow for direct mapping of the $H_{\text{el-ph}}$5- and mode-resolved EPC matrix via ultrafast, coherently-initiated dynamics, distinguishing local (Holstein) and nonlocal (SSH-type) couplings based on the energy dependence and vanishing at specific $H_{\text{el-ph}}$6 (Qu et al., 2023).

4. EPC in Model Systems: Materials-Specific Insights

Kagome Metals and 1D Carbon Systems

In kagome metal CsV$H_{\text{el-ph}}$7Sb$H_{\text{el-ph}}$8, ARPES measurements find an intermediate $H_{\text{el-ph}}$9 (Sb 5p and V 3d bands) supporting conventional $$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$0 estimates up to $$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$1 K, demonstrating that BCS pairing is viable, though intertwined with charge/spin/nematic orders (Zhong et al., 2022). Ab initio studies of (3,0) carbon nanotubes report exceptionally strong mode-resolved EPC—dominated by three breathing/stretching modes—enabling a $$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$2 of $$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$3 K at ambient pressure, the highest among elemental 1D systems (Ouyang et al., 5 Nov 2025).

Transition Metal Oxides and Nickelates

In pressurized La$$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$4Ni$$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$5O$$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$6, unique out-of-plane A$$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$7 breathing phonons couple selectively to Ni $$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$8 orbitals, but the total $$H_{\text{el-ph}} = \sum_{mn,\mathbf{k},\mathbf{q},\nu} g_{mn,\nu}(\mathbf{k},\mathbf{q})\, c^\dagger_{m,\mathbf{k}+\mathbf{q}} c_{n,\mathbf{k}} (b^\dagger_{\mathbf{q}\nu} + b_{-\mathbf{q}\nu}).$$9 remains below the superconducting threshold. Nevertheless, cooperation between EPC and frustrated electronic $$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$0 pairing (revealed by functional renormalization group) exponentially boosts $$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$1 by releasing pairing frustration (Zhan et al., 2024). In infinite-layer LaNiO$$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$2, antiferromagnetic ordering strongly enhances low-frequency EPC (from $$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$3 to $$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$4), producing a clear kink at 15 meV in the spectral function—a spectroscopic hallmark of magnetism-enhanced EPC (Zhang et al., 17 Apr 2025).

Cuprates, Dichalcogenides, and Layered Superconductors

In cuprates, RIXS and Raman measurements consistently show strong, mode- and momentum-dependent EPC, especially in bond-buckling and bond-stretching phonons ($$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$5–$$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$6 meV, $$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$7–$$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$8), with the small-$$g_{mn,\nu}(\mathbf{k},\mathbf{q}) = \left \langle \psi_{m,\mathbf{k}+\mathbf{q}} \left | \Delta_{\mathbf{q},\nu} V_{\rm KS}\right | \psi_{n,\mathbf{k}} \right \rangle$$9 "forward-scattering" part robust across the phase diagram and enhanced near the pseudogap critical doping. This EPC contributes additively to $|n,\mathbf{k}\rangle$0-wave pairing, potentially raising $|n,\mathbf{k}\rangle$1 by $|n,\mathbf{k}\rangle$2–$|n,\mathbf{k}\rangle$3\% in concert with spin fluctuations (Peng et al., 2021, Zhang et al., 2013, Braicovich et al., 2019).

Analysis of 1H-NbSe$|n,\mathbf{k}\rangle$4 and 1T-VTe$|n,\mathbf{k}\rangle$5 shows that it is the $|n,\mathbf{k}\rangle$6-resolved, mode-specific EPC—and not just Fermi surface nesting—that drives the lattice instabilities and CDW formation, with phonon softening and gap opening localized by the structure of $|n,\mathbf{k}\rangle$7 in the Brillouin zone (Wang et al., 2022).

5. Advanced Computational and Spectroscopic Developments

• Exchange-Correlation Functionals: The r2scan meta-GGA functional substantially improves EPC predictions for $|n,\mathbf{k}\rangle$8/$|n,\mathbf{k}\rangle$9-electron oxides and non-polar materials, capturing Fröhlich polar coupling, dielectric properties, and phonon stability without requiring empirical $|m,\mathbf{k+q}\rangle$0 corrections, as validated for CoO, NiO, and MgB$|m,\mathbf{k+q}\rangle$1 (Wang et al., 2024).
• Dynamical Mean-Field Theory (DMFT): Inclusion of local dynamical electron correlation via DFT+DMFT modifies the EPC vertex $|m,\mathbf{k+q}\rangle$2, leading to substantial frequency-dependent corrections absent in static $|m,\mathbf{k+q}\rangle$3 or hybrid-DFT. In SrVO$|m,\mathbf{k+q}\rangle$4, this increases the Jahn-Teller coupling by nearly $|m,\mathbf{k+q}\rangle$5 and shifts phonon lifetimes and electron scattering rates, providing direct evidence for the necessity of dynamical many-body treatments in correlated metals (Abramovitch et al., 6 May 2025).
• Twisted Bilayer Graphene (tBLG): Near the magic angle, EPC is strongly enhanced due to the resonance of the narrow bandwidth and dominant low-frequency phonons (5–13 meV, especially layer-breathing and shearing modes at $|m,\mathbf{k+q}\rangle$6), supporting $|m,\mathbf{k+q}\rangle$7 K in the appropriate $|m,\mathbf{k+q}\rangle$8 window. The magnitude and profile of $|m,\mathbf{k+q}\rangle$9 directly follow the strength with which phonons modulate the moiré potential (Zhu et al., 2024, Wang et al., 2024).

6. Functional Consequences and Ultrafast Control

Beyond equilibrium states, EPC provides a channel for coherent manipulation of band structure:

• Tr-HHG (time-resolved high-harmonic generation) reveals the direct action of coherent phonons on bandgap modulation, with both energy and intensity domains encoding absolute EPC strength $\nu$0. The phase and strength of EPC can be coherently manipulated by varying pump polarization, demonstrating full optical control on $\nu$1 fs timescales (Zhang et al., 2024).
• Two-dimensional EPC spectroscopy resolves mode-resolved, electronic-energy-dependent EPC, enabling direct discrimination between nonlocal SSH and local Holstein mechanisms via the $\nu$2 dependence of the matrix element $\nu$3. This paves the way for targeted ultrafast phonon control and the design of mode-selective, photoinduced collective phenomena (Qu et al., 2023).

7. Open Challenges and Outlook

Current frontiers involve accurate treatment of EPC in strongly correlated and topological systems, quantitative mapping of mode- and momentum-dependent EPC by advanced spectroscopy, and development of machine learning–assisted EPC evaluations for high-throughput materials discovery (Merkel et al., 28 Jul 2025, Zhong et al., 2023). In unconventional superconductors, the precise quantification of EPC—and its interplay with electronic correlations—remains pivotal for distinguishing the role of phonons in the mechanism of high-$\nu$4 pairing.

The combination of first-principles calculations, spectroscopic benchmarking, and ultrafast control schemes is rapidly refining our microscopic understanding of EPC, connecting fundamental many-body theory with practical routes to engineer quantum materials for next-generation device functionality.