Huang–Rhys Factor in Electron–Phonon Coupling

• The Huang–Rhys factor is a dimensionless parameter that quantifies electron–phonon coupling by measuring the average number of phonons involved in electronic transitions.
• It is derived from the configuration-coordinate model and computed using first-principles methods, DFPT, and bonding-descriptor approaches to determine reorganization energies.
• The factor governs spectroscopic signatures and nonradiative recombination rates, critically impacting the performance of quantum emitters and optoelectronic devices.

The Huang–Rhys factor, commonly denoted as $S$, is a central dimensionless parameter quantifying electron–phonon coupling strength in solids, molecules, and nanostructures. It rigorously measures the average number of phonons emitted or absorbed in a given optical or electronic transition, encoding the degree of nuclear reorganization induced by electronic excitation or carrier trapping. The value of $S$ not only dictates spectroscopic observables such as the zero-phonon line (ZPL) fraction, phonon sideband intensity, and line broadening, but also fundamentally controls nonradiative processes, recombination rates, and the effectiveness of quantum emitter platforms across material classes from bulk crystals to quantum dots and 2D solids.

1. Theoretical Formulation

The canonical theoretical basis for the Huang–Rhys factor arises from the configuration-coordinate (CC) model, originally introduced by Huang and Rhys for multiphonon processes in point defects and extended to electron–phonon, exciton–phonon, and polariton–phonon interactions in a broad array of systems. In a harmonic approximation, the adiabatic potential energy surfaces of two electronic or charge states are represented as parabolas along a mass-weighted collective coordinate $Q$. The dimensionless Huang–Rhys factor $S$ is defined by

$$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$

where $M$ is the effective mass, $\omega$ is the effective phonon frequency, $\Delta Q$ is the equilibrium geometry shift between the two electronic states, and $\Delta E$ is the reorganization energy—the vertical energy gap at fixed $Q$ between the two parabolic potentials (Whalley et al., 2021).

More generally, in the multidimensional case the total $S$0 is a sum over all normal modes $S$1: $S$2 where $S$3 is the projection of the geometry displacement onto the $S$4th normal mode with frequency $S$5 (Turiansky et al., 13 Jun 2025, Hassanzada et al., 2020).

$S$6 can equivalently be interpreted as the Poisson-mean of the phonon number distribution in the Franck–Condon progression. For a single mode, the rate of a multiphonon process scales as $S$7, directly controlling the spectral shape and sideband intensities.

2. Methods of Determination and Computational Approaches

Direct First-Principles Approaches

Ab initio methods (hybrid functionals, constrained-occupation DFT) are used to relax ground- and excited-state geometries, compute normal modes and frequencies, and evaluate the mass-weighted displacements required for $S$8 (Hassanzada et al., 2020, Biktagirov et al., 2024). Mode-resolved $S$9 contributions are obtained by projecting the displacement vector onto computed phonon eigenvectors via density-functional perturbation theory (DFPT).

An alternative protocol, suitable when excited-state relaxations are unstable or too expensive, projects excited-state forces at the ground-state geometry onto the ground-state normal modes, using

$Q$0

and assembles $Q$1. $Q$2 can be converged to within $Q$3 by restricting to a small local subspace of first- and second-neighbor atomic displacements, providing computationally efficient estimates (Turiansky et al., 13 Jun 2025).

An upper bound for $Q$4 is provided by the accepting-mode (1D) approximation: $Q$5 where $Q$6 is the average frequency along the total relaxation vector, and $Q$7 is the global displacement. This bound is especially useful for high-throughput screening (Turiansky et al., 13 Jun 2025).

Bonding-Descriptor Based Methods

A recently established methodology connects $Q$8 to differences in bonding character, quantified by Crystal Orbital Hamilton Population (COHP) analysis. The excited-state force at the ground-state geometry is approximated via

$Q$9

with $S$0 a scaling factor. GERD ("ground–excited reflective deformation") enables extraction of the effective displacement and hence $S$1 from ground-state calculations alone, with semi-quantitative accuracy over a wide range of $S$2 (Habis et al., 28 Oct 2025).

Spectroscopic and Experimental Determination

$S$3 is extracted from the intensity ratio of ZPL to phonon sidebands (Debye–Waller factor), the area or shape of the photoluminescence spectrum, or from temperature-dependent linewidth and Stokes shift: $S$4

or, where area fraction $S$5 in the ZPL is measured, $S$6 (Bhunia et al., 2023). In single-molecule STM-induced luminescence, $S$7 is obtained from the Poisson progression of peaks in the $S$8 tunneling spectrum (Wen et al., 2023).

In some parity-forbidden transitions, the standard Huang–Rhys theory breaks down and must be generalized. The sideband intensity distribution is no longer a simple Poisson, and $S$9 is extracted from the ratio of intensities in odd-order sidebands: $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$0 where $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$1 and $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$2 are the first and third phonon replica intensities (Wang et al., 1 Sep 2025).

3. Physical Interpretation and Regimes

The value of $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$3 encapsulates the degree of electron–phonon coupling:

• $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$4: Weak coupling regime. Most emission or absorption occurs in the ZPL; sidebands are negligible. Quantum emitters in this regime exhibit high photon indistinguishability, narrow linewidth, and long coherence (Bhunia et al., 2023, Biktagirov et al., 2024).
• $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$5: Moderate coupling. Pronounced but not dominant sidebands. The ZPL carries a significant fraction of total emission.
• $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$6: Strong coupling. The ZPL is exponentially suppressed, with emission dominated by multiphonon sidebands. This regime leads to rapid vibronic relaxation and is characteristic of soft lattices, defective perovskites, and polaronic semiconductors (Whalley et al., 2021).

The tailored design of $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$7 is crucial in optimizing defect-based single-photon sources, minimizing nonradiative losses in photonic devices, and tuning dynamics in molecular polaritonics and quantum-dot systems.

4. Representative Systems and Observed Values

The following table summarizes characteristic $$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$8 values, their physical origin, and associated Debye–Waller factors, as established in key studies:

System/Defect	$$S = \dfrac{M\omega}{2\hbar}(\Delta Q)^2 = \dfrac{\Delta E}{\hbar\omega}$$9	DW Factor	Reference
I$M$0 in CH$M$1NH$M$2PbI$M$3	$M$4350	$M$5	(Whalley et al., 2021)
C$M$6C$M$7 in h-BN (single-photon)	0.6	$M$850%	(Bhunia et al., 2023)
C$M$9–C$\omega$0 dimer in h-BN	2.0	$\omega$114%	(Mackoit-Sinkeviciene et al., 2019)
PbS quantum dots	1.7–2.5	--	(Wang et al., 2015)
2D-SiC Stone–Wales defect (SW$\omega$2)	0.94	39%	(Hassanzada et al., 2020)
(6,5) SWCNT helical Stone–Wales defect	0.74	$\omega$350%	(Biktagirov et al., 2024)
Sr$\omega$4Eu$\omega$5BP$\omega$6O$\omega$7 (Eu$\omega$8)	2.96–3.98	$\omega$92–5%	(Hou et al., 2013)
Mn$\Delta Q$0 parity-forbidden transitions	$\Delta Q$1 (from $\Delta Q$2)	--	(Wang et al., 1 Sep 2025)

In halide perovskites, $\Delta Q$3 can be exceptionally high (e.g., $\Delta Q$4 for iodine interstitial), leading to ultrafast nonradiative carrier capture and signifying a strong correlation between soft lattice dynamics and defect-induced losses (Whalley et al., 2021). In color centers in h-BN and low-dimensional materials, low $\Delta Q$5 (e.g., $\Delta Q$6) is associated with highly coherent single-photon emission (Bhunia et al., 2023, Biktagirov et al., 2024). In quantum dots, strong electron–phonon coupling (e.g., PbS QDs with $\Delta Q$7) is manifested in Franck–Condon blockade regimes and substantial sideband structure (Wang et al., 2015).

5. Implications for Material Functionality and Spectroscopy

The Huang–Rhys factor is directly linked to device-relevant photophysical properties:

• Zero-phonon fraction: The Debye–Waller factor, $\Delta Q$8, quantifies the ZPL intensity fraction—a key metric for indistinguishable single-photon sources (Bhunia et al., 2023, Hassanzada et al., 2020, Biktagirov et al., 2024).
• Carrier capture and recombination: In soft and defect-rich semiconductors (e.g., halide perovskites), large $\Delta Q$9 amplifies multiphonon-mediated nonradiative carrier capture, suppressing optoelectronic efficiency. Engineering stiffer frameworks or reducing octahedral rotations can mitigate $\Delta E$0 (Whalley et al., 2021).
• Exciton/polariton dynamics: In molecular cavity QED, $\Delta E$1 governs ultrafast population transfer, polariton formation, and photonic yield. Large $\Delta E$2 results in light–matter decoupling and multipolariton states (Wei et al., 2022, López et al., 16 Sep 2025, Bittner et al., 2021).
• Spectral line shape: The shape (width, asymmetry, progression) of optical emission/absorption directly reflects the magnitude and distribution of $\Delta E$3 across modes, enabling model-free extraction of $\Delta E$4 from lineshape fitting (Hassanzada et al., 2020, Shen et al., 2023, Hou et al., 2013).
• Selection rules: In parity-forbidden transitions, the structure of vibronic sidebands and the intensity ratios of odd/even phonon progressions provide a direct, model-independent measure of $\Delta E$5 via Herzberg–Teller extensions (Wang et al., 1 Sep 2025).

6. Advanced Generalizations: Polaritonic and Nonadiabatic Huang–Rhys Factors

Beyond standard electron–phonon coupling, the Huang–Rhys concept has been rigorously generalized:

• Polaritonic Huang–Rhys factor ($\Delta E$6): In macroscopic QED, $\Delta E$7 quantifies the cumulative displacement of photonic (polaritonic) modes induced by changes in molecular permanent dipoles upon excitation. It provides an explicit measure for light–matter decoupling in optical cavities, entering directly into corrected coupling strengths and polaritonic progression spectra (Wei et al., 2022).
• Sign and phase of $\Delta E$8 in nonadiabatic contexts: When including nonadiabatic couplings and cavity quantum electrodynamics, both the sign and magnitude of $\Delta E$9 dictate the order and efficiency of photonic energy conversion, determining the regions of the nuclear coordinate where rotating- and counter-rotating light-induced crossings occur (López et al., 16 Sep 2025).

7. Summary and Design Considerations

The Huang–Rhys factor $Q$0 is a universal gauge of electron–phonon coupling, with values spanning more than two orders of magnitude across materials and applications. It critically determines:

• The balance between coherent (ZPL) and incoherent (sidebands or nonradiative) emission.
• The efficiency of single-photon emitters and optoelectronic devices.
• The dynamical pathways for nonadiabatic energy redistribution in hybrid light–matter systems.

Analytical, spectroscopic, and computational advances now enable accurate, rapid prediction and measurement of $Q$1, including bonding-descriptor approaches suitable for high-throughput screening (Habis et al., 28 Oct 2025), force-based models requiring minimal excited-state input (Turiansky et al., 13 Jun 2025), and expanded theory for parity-forbidden transitions (Wang et al., 1 Sep 2025). The rational engineering of $Q$2, via chemical design and defect/environment selection, is thus a key actionable parameter for optimizing quantum photonic and optoelectronic systems.