Huang–Rhys Factor: Electron–Phonon Coupling

• Huang–Rhys Factor is a dimensionless measure that quantifies electron–phonon coupling by estimating the average number of phonons involved in an electronic transition.
• It directly affects spectral line shapes and energy transfer mechanisms in systems ranging from light-harvesting complexes to quantum dots.
• Both advanced computational techniques and experimental methods, such as force-based approximations and spectral decomposition, are used to calculate and extract its value.

The Huang–Rhys factor is a dimensionless parameter that quantifies the strength of electron–phonon coupling in a wide variety of condensed-matter and molecular systems. It expresses the degree to which an electronic excitation, transition, or transport event is accompanied by vibrational (phononic) excitations due to the reorganization of nuclear coordinates. Its value underlies fundamental phenomena ranging from photoluminescence line shapes, transport blockades in quantum dots, nonradiative defect capture, to energy transfer in light-harvesting complexes.

1. Mathematical Definition and Physical Meaning

The Huang–Rhys factor, often denoted as $S$, arises naturally in the shifted-harmonic-oscillator (Franck–Condon) model. For a system in which an electronic transition shifts the equilibrium position of a vibrational mode of frequency $\omega$, the coupling is characterized as: $S = \frac{M \omega}{2\hbar} (\Delta x)^2$ where $M$ is the effective mass, $\hbar$ Planck’s constant, and $\Delta x$ the shift in equilibrium coordinate between ground and excited electronic states (Wang et al., 2015). Equivalently, for electron–phonon coupling constant $g$,

$S = \left( \frac{g}{\hbar\omega} \right)^2$

For systems with multiple modes, $S$ is additive: $$S = \sum_k S_k, \quad S_k = \frac{\omega_k q_k^2}{2\hbar}$$ where $q_k$ is the mass-weighted coordinate displacement in mode $k$ (Schulze et al., 2017, Mackoit-Sinkeviciene et al., 2019, Turiansky et al., 13 Jun 2025).

Physically, $S$ is the average number of phonons emitted or absorbed during an electronic transition. For $S \ll 1$, transitions are mostly “zero-phonon”; for $S \gtrsim 1$, multiphonon processes dominate.

2. Manifestation in Spectral Densities and Line Shapes

The Huang–Rhys factor determines the Franck–Condon progression in both absorption and emission spectra. The intensities of vibrational sidebands follow a Poisson distribution: $I_n \propto e^{-S} \frac{S^n}{n!}$ with $I_0$ corresponding to the zero-phonon line (ZPL) and $I_{n > 0}$ to phonon sidebands. The Debye–Waller factor, $w_{\mathrm{ZPL}} = e^{-S}$, gives the ZPL fraction (Mackoit-Sinkeviciene et al., 2019, Bhunia et al., 2023). A large $S$ shifts oscillator strength from the ZPL to the sidebands, broadening the emission and reducing photon indistinguishability, which is crucial in single-photon emitters (SPEs) (Bhunia et al., 2023, Hassanzada et al., 2020).

In models with continuous phonon baths, $S$ is related to the spectral density $J(\omega)$ via: $$S = \int_0^\infty \frac{J(\omega)}{\omega^2} d\omega$$ This form connects $S$ directly to the total area under the normalized spectral density, independent of its detailed shape (Schulze et al., 2017, Shen et al., 2023).

3. Methodologies for Calculation and Experimental Extraction

Calculation of $S$ requires knowledge of excited- and ground-state geometry and vibrational mode structure. Standard ab initio approaches—$\Delta$SCF, constrained-DFT, or vibrational analysis—yield the mass-weighted displacements $\Delta Q_k$ and frequencies $\omega_k$, entering the sum for $S$ (Mackoit-Sinkeviciene et al., 2019, Turiansky et al., 13 Jun 2025, Habis et al., 28 Oct 2025).

To avoid the computational cost of full excited-state relaxations and phonon calculations:

• Force-based approximations: $S$ can be estimated from excited-state forces at the ground-state geometry, using the force-mode or accepting-mode approaches (Turiansky et al., 13 Jun 2025).
• Orbital-based descriptors: The change in chemical bonding character, quantified via crystal orbital Hamilton populations (COHP), provides a rapid estimate of $S$ and can be used in computational screening pipelines (Habis et al., 28 Oct 2025).

Experimental extraction commonly employs:

• Spectral decomposition: From photoluminescence, $S = -\ln(\text{ZPL}/[\text{ZPL}+\text{PSB}])$ (Bhunia et al., 2023, Hassanzada et al., 2020).
• Tunneling spectroscopy: In quantum dots, the width of the Franck–Condon blockade in electron tunneling gives $S \approx \sqrt{e V_{\mathrm{gap}}/\hbar \omega}$ (Wang et al., 2015).
• STM-induced luminescence: The peak area ratios in $d^2I/dV^2$ match the Poisson form, allowing $S$ to be fitted or directly extracted as $S \approx I_1^{(2)}/I_0^{(2)}$ (Wen et al., 2023).
• Parity-forbidden systems: For transitions dominated by Herzberg–Teller vibronic coupling, $S$ is extracted not from the omitted ZPL but from the ratio of third- to first-order sideband intensities as $S = \sqrt{2I_3/(9I_1)}$ (Wang et al., 1 Sep 2025).

4. Role in Energy Transfer, Carrier Capture, and Photo-Physics

The Huang–Rhys factor is central to dynamics in light-harvesting complexes, quantum transport, and nonradiative processes:

• Exciton Transport: In pigment–protein complexes such as the FMO complex, the total $S$ determines the timescale of vibrational relaxation and exciton trapping. Large $S$ accelerates energy transfer to the lowest-energy site and damps coherence more rapidly; the detailed mode structure modulates vibrational and vibronic population pathways (Schulze et al., 2017, Calderón et al., 2024).
• Carrier Trap Centers: In semiconductor defects, $S$ sets the magnitude of the nonradiative capture coefficient. Very large values (e.g., $S > 300$ for iodine interstitials in halide perovskites) indicate strong multi-phonon emission and fast nonradiative recombination, directly impacting device performance (Whalley et al., 2021).
• Franck–Condon Blockade: In QDs and molecular junctions, strong electron–phonon coupling ($S \gtrsim 1$) produces a bias threshold below which conductance is suppressed; the gap scales as $S^2\hbar\omega$ (Wang et al., 2015).

5. Variants: Polaritonic Huang–Rhys Factor and Extensions

Beyond electron–phonon coupling, the analogy with vibrational physics has motivated the definition of polaritonic Huang–Rhys factors. Here, the displacement is in the field coordinate due to light–matter coupling and permanent dipoles: $$S_{\mathrm{pol}} = \frac{1}{\hbar\pi\epsilon_0 c^2}\int_0^\infty d\omega~\Delta\mu\cdot \mathrm{Im}~\overline{G}(r_M, r_M, \omega)\cdot \Delta\mu$$ This quantity determines the reduction in effective light–matter coupling and the appearance of multipolariton transitions, analogous to multi-phonon sidebands (Wei et al., 2022, López et al., 16 Sep 2025). In cavity QED, $S_{\mathrm{pol}}$ may become significant in ultrastrong regimes, modifying Rabi splittings and leading to phenomena such as light–matter decoupling, strong enhancement of two-photon processes, and non-radiative progression (Wei et al., 2022).

In parity-forbidden transitions and systems that require tensorial transition operators, the standard theory fails, and Herzberg–Teller (linear vibronic) terms dominate. The modified analytic structure of vibrational intensities necessitates extraction formulas for $S$ based on the relative intensities of odd sidebands (Wang et al., 1 Sep 2025).

6. Representative Values Across Systems

System / Material	S	Method / Source
PbS quantum dots	1.7–2.5	Franck–Condon blockade (Wang et al., 2015)
FMO complex (biological pigment sites)	0.42–0.96	Exciton modeling (Schulze et al., 2017)
h-BN single-photon emitter (C₂C_N)	0.6	Debye–Waller ratio (Bhunia et al., 2023)
h-BN C₂ dimer defect	~2	First-principles (Mackoit-Sinkeviciene et al., 2019)
2D-SiC Stone–Wales defects	0.7–1.7	Ab initio (Hassanzada et al., 2020)
MAPbI₃: Iodine interstitial	> 300	First-principles (Whalley et al., 2021)
TMDs (MoS₂, WSe₂, etc.)	0.1–1	MBO fitting (Shen et al., 2023)
Mn$^{4+}$ parity-forbidden transition	10⁻³–10⁻²	Sideband ratio (Wang et al., 1 Sep 2025)

These values illustrate the broad range of $S$ encountered, with strong implications for quantum efficiency, spectral purity, and dynamic timescales.

7. Significance and Design Implications

The value of the Huang–Rhys factor is a key design criterion for light-emitting devices, quantum emitters, and defect-engineered spin qubits. Low $S$ ($\lesssim 1$) indicates weak coupling, supporting a sharp ZPL and high photon indistinguishability—a requirement for quantum networks and cavity integration (Bhunia et al., 2023, Hassanzada et al., 2020). High $S$ underpins efficient nonradiative trapping, energy dissipation, and transport control, but at the cost of emission broadening or transport blockage (Whalley et al., 2021).

Novel computational techniques leveraging chemical bonding changes (Habis et al., 28 Oct 2025), local phonon mode projections (Turiansky et al., 13 Jun 2025), and rapid screening approaches have emerged to rationalize and predict $S$ across large materials spaces. These developments support high-throughput identification of materials with tailored electron–phonon coupling, advancing defect engineering, photonics, and quantum information science.

In summary, the Huang–Rhys factor provides a quantitative unifying framework for describing, measuring, and engineering vibronic and polaritonic effects in solids, molecules, and hybrid systems. Its theoretical calculation, experimental determination, and manipulation are central to contemporary condensed-matter, nano-optics, and quantum technologies.