Huang–Rhys Theory Overview

• Huang–Rhys theory is a displaced-oscillator framework that quantifies electron–phonon (vibronic) coupling using a dimensionless Huang–Rhys factor.
• It models electronic transitions by treating ground and excited states as displaced harmonic potential surfaces, linking Franck–Condon overlaps with vibrational sidebands.
• The theory underpins both spectroscopic analysis and computational modeling to extract reorganization energies, phonon sideband intensities, and thermal effects in diverse systems.

Searching arXiv for recent and foundationally relevant Huang–Rhys literature to support a comprehensive article. Huang–Rhys theory is the standard displaced-oscillator framework for describing how electronic transitions couple to vibrational degrees of freedom in molecules, defects, and solids. In its classical form, it treats the ground and excited electronic states as harmonic potential-energy surfaces displaced along one or more normal coordinates, so that an optical transition samples Franck–Condon overlaps rather than a purely electronic line. The central quantity is the Huang–Rhys factor $S$, a dimensionless measure of electron–phonon or vibronic coupling strength that controls the partition of spectral weight between the zero-phonon line (ZPL) and phonon sidebands, the reorganization energy, and, in many applications, the propensity for radiative broadening or nonradiative multiphonon relaxation (Hashemi et al., 2020, Libbi et al., 2021).

1. Classical framework and physical interpretation

The classical Huang–Rhys construction is the displaced harmonic oscillator model under the Franck–Condon and Condon approximations. Electronic transitions are taken to occur vertically, with nuclei fixed during the optical event, and the transition dipole is assumed independent of nuclear coordinates. In solids this is commonly phrased as two Born–Oppenheimer surfaces approximated by parallel parabolas in multidimensional nuclear configuration space, each associated with harmonic phonon modes (Hashemi et al., 2020, Libbi et al., 2021).

For a single effective mode of frequency $\omega$, the Huang–Rhys factor can be written as

$S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$

or, equivalently, as the ratio of the reorganization energy to the mode quantum,

$S=\frac{\lambda}{\hbar\omega}.$

In the multimode form used for defects in solids, the mode-resolved factor is

$$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$

with the mass-weighted projection

$$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$

In dimensionless normal coordinates $q_k$, this becomes $S_k=\tfrac12(\Delta q_k)^2$ (Hashemi et al., 2020).

The physical meaning of $S$ is twofold. At $T=0$, it equals the average number of phonons emitted in the transition. It also measures the displacement between the equilibrium nuclear configurations of the two electronic states, so small $\omega$0 implies weak structural relaxation and a strong ZPL, whereas large $\omega$1 implies strong lattice reorganization, a dominant phonon sideband, and a strongly reduced Debye–Waller factor (Hashemi et al., 2020, Bhunia et al., 2023).

Historically, the classical theory is tied to Huang and Rhys’s Poisson-like sideband law under the Condon approximation, with later Herzberg–Teller generalizations becoming necessary when parity-forbidden transitions invalidate a coordinate-independent transition dipole (Wang et al., 1 Sep 2025).

2. Core quantities, line shapes, and finite-temperature generalization

The main observables of Huang–Rhys theory are compactly summarized by a small set of relations.

Quantity	Expression	Meaning
Partial HR factor	$\omega$2	Coupling of mode $\omega$3
Total HR factor	$\omega$4	Total vibronic coupling
Debye–Waller factor	$\omega$5	ZPL fraction at $\omega$6
Reorganization energy	$\omega$7	Relaxation energy
One-mode progression	$\omega$8	$\omega$9-phonon sideband intensity

In the one-dimensional displaced-oscillator limit, the vibronic progression is Poissonian: $S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$0 As $S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$1 increases, spectral weight is transferred from the ZPL into the phonon sideband, which becomes broader and smoother. The displaced-parabola model also relates the Stokes shift to the reorganization energy by $S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$2 (Hashemi et al., 2020, Bhunia et al., 2023).

The multimode generalization is most naturally expressed through a spectral function,

$S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$3

and a generating-function or cumulant representation of the emission profile. One form used in first-principles defect calculations is

$S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$4

with $S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$5 the Fourier transform of $S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$6, and the emission intensity proportional to $S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$7. A finite-temperature cumulant form writes

$S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$8

where

$S=\frac{M\omega}{2\hbar}\,\Delta Q^2,$9

The corresponding finite-temperature ZPL weight is

$S=\frac{\lambda}{\hbar\omega}.$0

These expressions preserve the classical displaced-oscillator picture while admitting realistic multimode spectra and thermal broadening (Hashemi et al., 2020).

A related time-domain formulation appears in two-dimensional spectroscopy, where the energy-gap correlation function is expanded in Huang–Rhys factors. In that setting, $S=\frac{\lambda}{\hbar\omega}.$1 governs oscillatory modulations of rephasing and non-rephasing contributions, and the vibrational line-shape function enters the third-order response directly (Mancal et al., 2010).

3. Spectroscopic interpretation and experimental extraction

In experiment, Huang–Rhys theory is usually encountered through the decomposition of a spectrum into a ZPL and a phonon sideband. The Debye–Waller factor is defined as the fraction of emission in the ZPL,

$S=\frac{\lambda}{\hbar\omega}.$2

and in the single-mode or weak-coupling approximation one uses

$S=\frac{\lambda}{\hbar\omega}.$3

This relation is exact in the single-mode displaced-oscillator model at zero temperature and remains a common working approximation in weakly multimode situations (Bhunia et al., 2023).

A room-temperature example is provided by deterministic single-photon emitters in layered h-BN. The reported emitter shows a ZPL at $S=\frac{\lambda}{\hbar\omega}.$4 nm, a phonon sideband centered at $S=\frac{\lambda}{\hbar\omega}.$5 nm, a Debye–Waller factor $S=\frac{\lambda}{\hbar\omega}.$6, and a Huang–Rhys factor $S=\frac{\lambda}{\hbar\omega}.$7 extracted from $S=\frac{\lambda}{\hbar\omega}.$8. Raman measurements identify the $S=\frac{\lambda}{\hbar\omega}.$9 mode at $$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$0, corresponding to $$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$1–$$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$2 meV, matching the ZPL–PSB separation and supporting a dominant optical-phonon picture (Bhunia et al., 2023).

The same paper also uses a displaced-oscillator Stokes-shift estimate,

$$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$3

to obtain $$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$4, consistent with the Debye–Waller extraction. That agreement is typical when one dominant optical mode controls the vibronic structure and thermal corrections do not overwhelm the ZPL fraction (Bhunia et al., 2023).

Huang–Rhys factors can also be inferred from transport spectroscopy. In STM-induced luminescence, the inelastic electron-scattering mechanism yields step-like structure in $$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$5 and vibronically resolved peaks in $$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$6, whose relative amplitudes follow Franck–Condon factors. In the proposed low-temperature protocol,

$$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$7

and, when $$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$8 varies slowly with $$S_k=\frac{\omega_k}{2\hbar}(\Delta Q_k)^2, \qquad S=\sum_k S_k,$$9,

$$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$0

This directly links tunneling-threshold spectroscopy to Huang–Rhys parameter extraction (Wen et al., 2023).

4. Defects, luminescence, and the range of Huang–Rhys behavior in solids

Color centers in semiconductors provide some of the clearest realizations of Huang–Rhys physics. In 4H–SiC, first-principles lineshape calculations show a narrow ZPL followed by a phonon sideband that starts about $$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$1 meV below the ZPL, with widths of about $$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$2–$$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$3 meV correlated with $$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$4. Across different defects, low-energy bulk-like modes dominate the sideband at about $$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$5–$$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$6 meV, while localized modes near $$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$7–$$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$8 meV and $$\Delta Q_k=\sum_i \sqrt{m_i}\,\mathbf e_{i,k}\cdot\big(\mathbf R_i^e-\mathbf R_i^g\big).$$9–$q_k$0 meV sculpt defect-specific shoulders and high-energy features (Hashemi et al., 2020).

The quantitative spread in $q_k$1 and Debye–Waller factors across 4H–SiC is substantial. For negatively charged silicon vacancies, $q_k$2 is typically around $q_k$3 with $q_k$4, except for the $q_k$5 configuration with $q_k$6 and $q_k$7. The doubly negative silicon vacancy $q_k$8 (“W2”) has $q_k$9 and $S_k=\tfrac12(\Delta q_k)^2$0, whereas the carbon antisite–vacancy pair $S_k=\tfrac12(\Delta q_k)^2$1 (“A” transition) has $S_k=\tfrac12(\Delta q_k)^2$2 and $S_k=\tfrac12(\Delta q_k)^2$3. These values map directly onto application-relevant trade-offs: higher Debye–Waller factors favor bright ZPL emission for cavity coupling and spin–photon interfaces, while larger $S_k=\tfrac12(\Delta q_k)^2$4 broadens the PSB and suppresses the ZPL (Hashemi et al., 2020).

A parallel first-principles application appears in $S_k=\tfrac12(\Delta q_k)^2$5-SiAlON:$S_k=\tfrac12(\Delta q_k)^2$6, where the emission spectrum at $S_k=\tfrac12(\Delta q_k)^2$7 K is reproduced with $S_k=\tfrac12(\Delta q_k)^2$8. The phonon signature is robust across low-energy structural variants and consists of a dominant Eu-localized mode near $S_k=\tfrac12(\Delta q_k)^2$9 meV, a broad delocalized band between $S$0 and $S$1 meV, and a breathing-like mode near $S$2 meV. The weak-to-moderate coupling explains the persistence of resolved vibronic replicas even as the Al/O concentration $S$3 increases (Bouquiaux et al., 11 May 2026).

At the opposite extreme, nonradiative multiphonon capture in soft semiconductors can exhibit giant effective Huang–Rhys factors. For electron capture by the neutral iodine interstitial $S$4 in methylammonium lead iodide, the effective accepting mode has frequency $S$5, the effective Huang–Rhys factor is approximately $S$6, and the calculated anharmonic capture coefficient is $S$7. The low-energy accepting mode is localized around the defect with inverse participation ratio $S$8, reflecting the large lattice relaxation associated with octahedral rotations (Whalley et al., 2021).

These examples delimit the practical range of Huang–Rhys theory in solids: from low-$S$9 single-photon emitters with strong ZPL fractions to giant-$T=0$0 defect centers whose functionality is governed by multiphonon capture rather than photon emission.

5. Breakdown of the classical model and symmetry-based extensions

The classical Huang–Rhys model is predictive when the ZPL is symmetry allowed, the potential surfaces are approximately harmonic and parallel, the optical dipole is only weakly coordinate dependent, and a small number of independent modes dominate the coupling (Libbi et al., 2021). Several important systems violate those assumptions.

A central misconception is that a weak or absent ZPL necessarily implies a large Huang–Rhys factor. For the negatively charged boron vacancy $T=0$1 in monolayer hBN, many-body perturbation theory shows that the phononless transition is dark because the dominant exciton is symmetry forbidden in $T=0$2. The observed luminescence near $T=0$3 eV is activated by phonons of $T=0$4 and $T=0$5 symmetry, and static Jahn–Teller distortion alone does not recover the experimental intensity. In such a case, inferring $T=0$6 from $T=0$7 is physically misleading, because the ZPL suppression originates in symmetry rather than in large geometric reorganization (Libbi et al., 2021).

Parity-forbidden transitions require an explicit Herzberg–Teller extension. In Mn$T=0$8-doped fluoride phosphors, the conventional Condon approximation fails because the bare electric-dipole matrix element vanishes for the $T=0$9 transition. The linear Herzberg–Teller term introduces odd-parity coordinate dependence and yields a distinctive selection rule: odd-order sidebands are enhanced whereas even-order sidebands are strongly suppressed. In the weak-coupling, low-temperature, single-mode limit, the odd-order intensities obey

$\omega$00

and one obtains the practical estimator

$\omega$01

This is not a perturbative correction to classical Huang–Rhys theory but a change in the intensity law forced by symmetry (Wang et al., 1 Sep 2025).

A separate generalization transfers Huang–Rhys reasoning from nuclear coordinates to quantized electromagnetic degrees of freedom. The “polaritonic Huang–Rhys factor” treats permanent-dipole-induced displacements of macroscopic-QED field coordinates in analogy with vibronic displacements of nuclei, leading to

$\omega$02

and a modified light–matter coupling

$\omega$03

Here the displaced-oscillator analogy remains intact, but the “mode” is an electromagnetic continuum shaped by the dielectric environment rather than a lattice vibration (Wei et al., 2022).

6. Computational practice, approximations, and design use

Modern first-principles implementations of Huang–Rhys theory typically combine electronic-structure calculations with harmonic phonons and a generating-function lineshape. In 4H–SiC, for example, the workflow uses PAW-based density functional theory, PBEsol structural relaxations and phonons, HSE06 ZPL energies, constrained-occupancy excited states, force-constant fitting, and a generating-function construction of $\omega$04 from mode-resolved $\omega$05 (Hashemi et al., 2020). In $\omega$06-SiAlON:$\omega$07, an analogous strategy combines $\omega$08SCF excited states, embedded interatomic force constants in supercells up to $\omega$09 atoms, and a generating-function evaluation of the vibronic spectrum (Bouquiaux et al., 11 May 2026).

Because excited-state relaxations and full phonon calculations are often the computational bottleneck, several approximation strategies have emerged. One such method reconstructs the multidimensional coupling from excited-state forces evaluated only at the ground-state equilibrium geometry. In that framework, the vertical excitation energy is

$\omega$10

the modal displacements are obtained by projecting the excited-state force onto normal modes, and the total ZPL is approximated as

$\omega$11

A key result is that the ZPL can be approximated with just a single mode, whereas the Huang–Rhys factor converges by including displacements up to the second nearest neighbors. The same analysis proves that the accepting-mode Huang–Rhys factor is a strict upper bound on the full multidimensional Huang–Rhys factor (Turiansky et al., 13 Jun 2025).

A more chemically interpretable approach uses orbital bonding descriptors. For hBN defects and the diamond NV center, a COHP-based descriptor estimates excited-state forces from bonding-character differences between the initial and final defect orbitals, and a Ground-Excited Reflective Deformation procedure then reconstructs the excited-state geometry without an explicit excited-state relaxation. In this picture, small changes in bonding character imply small forces, small structural reorganization, and therefore small Huang–Rhys factors; conversely, bonding-to-antibonding changes tend to produce large $\omega$12 (Habis et al., 28 Oct 2025).

Huang–Rhys parameters also serve as design variables outside static defect spectroscopy. In organic microcavities, the upper-polariton to lower-polariton vibronic transfer scales with the molecular Huang–Rhys parameter, so small $\omega$13 suppresses nonadiabatic relaxation and can lock molecules near the Franck–Condon geometry for tens to thousands of picoseconds, whereas relaxed lower-polariton lifetimes remain near $\omega$14–$\omega$15 ps for the chosen cavity parameters (Bittner et al., 2021). In excitation-transport models with discrete intramolecular modes, keeping $\omega$16 fixed while increasing the acceptor vibrational frequency increases the equilibrium acceptor population, but under biologically relevant nonequilibrium steady-state conditions the enhancement becomes negligible once trapping and recombination are included (Calderón et al., 2024).

Taken together, these developments recast Huang–Rhys theory from a post hoc fitting tool into a quantitative descriptor of structure–spectrum relations, nonradiative kinetics, and strong-coupling dynamics. Its classical form remains indispensable, but its modern use is inseparable from symmetry analysis, many-body treatments where the Condon picture fails, and computational strategies that make $\omega$17, $\omega$18, and related reorganization energies accessible at scale (Libbi et al., 2021, Turiansky et al., 13 Jun 2025, Habis et al., 28 Oct 2025).