Electron–Phonon Coupling Constants

• Electron–phonon coupling constants are defined via the Eliashberg spectral function, linking electron interactions with phonon density of states in periodic systems.
• Experimental techniques like ARPES, time-resolved spectroscopy, and HAS extract precise λ values, elucidating phenomena such as superconductivity and charge-density-wave formation.
• Theoretical and computational models underscore λ's dependence on band structure, phonon spectrum, and quantum-geometric effects, guiding the design of advanced materials.

Electron–phonon coupling constants are fundamental parameters that characterize the strength and nature of interactions between conduction electrons and lattice vibrations (phonons) in solid-state systems. These constants play a central role in a vast array of phenomena, including conventional and unconventional superconductivity, charge-density-wave formation, electrical resistivity, and spectral renormalization effects observed in both equilibrium and time-resolved spectroscopies.

1. Fundamental Definitions and Physical Basis

The dimensionless electron–phonon coupling constant λ is rigorously defined within the Migdal–Eliashberg theory via the Eliashberg spectral function α²F(Ω). For a general, periodically ordered system,

$$\lambda = 2 \int_0^{\infty} \frac{\alpha^2 F(\Omega)}{\Omega} d\Omega,$$

where α²F(Ω) encodes the momentum- and mode-averaged electron–phonon matrix elements and the phonon density of states. Explicitly,

$$\alpha^2 F(\Omega) = \frac{1}{N(E_F)} \sum_{k,k',n,m,l} |G_{nml}(k,k')|^2\,\delta(E_F-E_{n}(k))\,\delta(E_F-E_{m}(k'))\,\delta(\Omega-\Omega_{l}(k'-k)),$$

with $G_{nml}(k,k')$ the electron–phonon matrix elements connecting bands n, m via phonon branch l. λ controls the electronic mass renormalization $m^* = m (1+\lambda)$ and imposes a velocity or “kink” renormalization in the quasiparticle dispersion. For local and band-resolved contexts, λ can acquire explicit dependence on momentum, band indices, or even energy (Faeth et al., 2021, Yu et al., 2023, Mazzola et al., 2016, Mahatha et al., 2018).

2. Experimental Extraction Methodologies

Extraction of λ exploits various physical observables linked to its underlying definition, across multiple experimental platforms:

Angle-Resolved Photoemission Spectroscopy (ARPES):

In systems such as monolayer FeSe/SrTiO₃, ARPES directly images main and replica bands induced by strong forward-scattering electron–phonon coupling. By analyzing the ratio of first replica to main quasiparticle spectral weight $I_1/I_0$, blue shifts in replica energy relative to the phonon threshold, and linewidth broadening consistent with self-energy calculations, λ can be determined (e.g., λ = 0.19 ± 0.02 for FeSe/SrTiO₃) (Faeth et al., 2021). In organic semiconductors like pentacene, ARPES temperature broadening of the HOMO band yields λ = 0.36 ± 0.05, with the dominant coupling attributed to soft intermolecular phonons (0908.4258).

Time-Resolved Optical and Photoemission Spectroscopy:

Nonequilibrium optical pump–probe methods (such as broadband transient reflectivity or time-resolved two-photon photoemission) provide time-domain access to thermalization rates between electrons and specific phonon branches, which—when modeled by multi-temperature (2TM/3TM) models—permit determination of selective or total λ values. In MgB₂, ultrafast electron relaxation directly tracks to selective coupling with the E₂g phonon, yielding λ_{E₂g} ≈ 0.56, about half the total λ_total ≈ 1.1 (Mor et al., 4 Mar 2025).

Debye–Waller Analysis of Helium Atom Scattering (HAS):

On metal surfaces and 2D materials, the thermal attenuation (Debye–Waller factor) of specular He diffraction quantifies the magnitude of λ, both in ultrathin films and single- or multi-layer systems. The theoretical framework connects the Debye–Waller exponent linearly to the mass enhancement parameter λ (Manson et al., 2016, Benedek et al., 2020, Benedek et al., 2017).

Raman and Transport Methods:

Electronic Raman scattering, by comparison of temperature-dependent continuum and relaxation rates with calculations incorporating band-structure and electron–phonon scattering, provides λ values consistent with other one-particle and optical measurements, notably for elemental metals and metallic compounds (Ponosov et al., 2012).

3. Model-Specific and Material-Specific Considerations

Electron–phonon coupling is strongly dependent not just on the electronic density of states, but also on the details of the phonon spectrum, geometric/topological band properties, and the spatial symmetry of the electron–phonon vertex:

• Forward-Scattering and Replica Bands: In low-density 2D systems strongly coupled to high-energy polar optical phonons, forward scattering can generate well-defined replica bands in $A(k, \omega)$ and provides a direct probe of λ distinct from traditional mass renormalization (Faeth et al., 2021).
• Band, Valley, and Spin Selectivity: In multiband systems such as MgB₂ and monolayer MoS₂, λ becomes channel-selective—e.g., λ ≈ 0.05 for the upper spin–orbit band at K in MoS₂ versus λ ≈ 0.32 for the lower branch due to allowed intervalley scattering (Mahatha et al., 2018, Mor et al., 4 Mar 2025). For MgB₂, selective coupling to σ- versus π- derived bands yields sharply different λ components (Mor et al., 4 Mar 2025, Yu et al., 2023).
• Quantum-Geometry Contributions: The total λ can acquire substantial quantum-geometric contributions, quantified by the Fubini–Study metric on the Fermi surface. In MgB₂, 90% of λ can be geometric in origin, whereas in graphene near charge neutrality, λ_{geo}/λ → 0.5 (Yu et al., 2023).
• Adiabatic and Nonadiabatic Regimes: Generalized Eliashberg–McMillan approaches allow for λ to be properly defined in both adiabatic (Ω₀ ≪ E_F) and antiadiabatic (Ω₀ ≫ E_F) regimes. The relevant small parameters controlling the theory and the physical interpretation of λ and the mass renormalization ($\tilde{\lambda}$) change accordingly (Sadovskii, 2018).

4. Numerical Values and Empirical Ranges

A cross-section of experimentally extracted and theoretically calculated λ values is summarized below:

Material/System	λ Value(s)	Principal Measurement	Reference
FeSe/SrTiO₃	0.19 ± 0.02	ARPES replica band analysis	(Faeth et al., 2021)
Pentacene (crystalline)	0.36 ± 0.05	ARPES thermal broadening	(0908.4258)
Double-wall CNTs	(5.4 ± 0.9) ×10⁻⁴	Ultrafast 2PPE, TTM modeling	(Chatzakis, 2013)
MgB₂ total	1.10 ± 0.10	ARPES, reflectivity, 3TM optical fits	(Mor et al., 4 Mar 2025)
MgB₂ (E₂g phonon)	0.56 ± 0.05	Time-resolved optical, 3TM decomposition	(Mor et al., 4 Mar 2025)
Graphene σ-band	0.6(1)–0.9	DFT, ARPES, self-energy Kramers–Kronig	(Mazzola et al., 2016)
SGL MoS₂ upper VB @K	0.05 ± 0.01	ARPES linewidth vs T	(Mahatha et al., 2018)
SGL MoS₂ lower VB @K	0.32 ± 0.01	ARPES linewidth vs T	(Mahatha et al., 2018)
Alkali-metal films	0.16–0.54	He atom scattering (HAS)	(Benedek et al., 2017)
Pb/Cu(111) films	0.9–1.2	HAS, ARPES	(Benedek et al., 2017)
IrGe	~1.4	Specific heat, McMillan inversion	(Hirai et al., 2018)

Values for elemental metals span from λ ≈ 0.13 (W) to ≳1 (Nb, Pb) (Ponosov et al., 2012). For carbon-based nanostructures, absolute values of λ can decrease by several orders of magnitude due to geometrical or screening effects (Chatzakis, 2013).

5. Theoretical Frameworks and Computation

First-principles calculations, predominantly based on density-functional perturbation theory (DFPT), provide a direct route to α²F(Ω), λ, and related superconducting parameters, with efficient implementations leveraging Fermi-surface averaging and special integration schemes to accelerate convergence (Koretsune et al., 2016). Generalized methodologies allow for extraction of λ from lattice-dynamical models, tight-binding molecular dynamics (far above equilibrium temperatures), or via explicit evaluation of self-energy derivatives (Medvedev, 2023, Yu et al., 2023):

• The canonical mass enhancement parameter arises as the low-frequency slope of the real part of the electron–phonon self-energy, λ = −∂ReΣ(ω)/∂ω|_{ω=0} (Carbotte et al., 2012).
• For strong-coupling or quantum-geometrical regimes, λ splits naturally into an energetic/kinetic and a geometric/topological part (Yu et al., 2023).

6. Surface, Interface, and Reduced-Dimensionality Effects

The extension of λ to surfaces (e.g., ultrathin overlayers, 2D systems, and topological insulators) is experimentally realized via HAS Debye–Waller analysis, which can resolve the evolution of λ with layer thickness, boundary conditions, and substrate binding. For example, λ extrapolates to 0.89 (high-T model) or 0.32 (Einstein model) for free-standing graphene under cyclic boundary conditions (Benedek et al., 2020). On the surface of topological insulators such as Bi₂Te₃, optical phonon bands of narrow energy width (3–7 meV) dominate, leading to large, sharp λ ≈ 2.0, with agreement between electron and phonon spectroscopy perspectives (Howard et al., 2013).

7. Non-Equilibrium and High-Temperature Phenomena

In high-excitation, non-equilibrium, and nano-confinement regimes, the effective electron–phonon coupling constant can deviate from equilibrium values by up to an order of magnitude, primarily due to temporal and spatial non-equilibrium between electron and various phonon branches. Dynamical, nonperturbative models indicate that such effects must be included to correctly model ultrafast thermalization and transport in micro- and nano-scale systems (Miao et al., 2020, Medvedev, 2023). In semiconductors at high electronic temperature, band-resolved and nonequilibrium methods converge to within ~35% for practical purposes (Medvedev, 2023).

---

References:

(0908.4258, Ponosov et al., 2012, Carbotte et al., 2012, Chatzakis, 2013, Howard et al., 2013, Manson et al., 2016, Mazzola et al., 2016, Koretsune et al., 2016, Benedek et al., 2017, Hirai et al., 2018, Sadovskii, 2018, Mahatha et al., 2018, Benedek et al., 2020, Miao et al., 2020, Faeth et al., 2021, Yu et al., 2023, Medvedev, 2023, Mor et al., 4 Mar 2025)

The electron–phonon coupling constant λ remains a central, experimentally and theoretically tractable parameter for quantifying electron–lattice interactions, subject to refinement via new probes, materials, and theoretical insights into geometric and nonadiabatic effects.