Polaritonic Huang–Rhys Factor

• Polaritonic Huang–Rhys factor is a dimensionless metric that measures the effective displacement of electromagnetic modes resulting from molecular dipole changes during excitation.
• It modifies light–matter coupling by exponentially suppressing radiative transitions, thereby controlling multiphoton emission and spectral progression in various dielectric environments.
• Its analytical and numerical predictions within macroscopic QED frameworks enable precise experimental design in cavity QED and polaritonic chemistry applications.

The polaritonic Huang–Rhys factor ($S_p$) generalizes the well-established vibrational Huang–Rhys parameter to systems in which molecular excitations interact with quantized electromagnetic (EM) fields. $S_p$ quantifies the effective displacement of EM modes induced by changes in a molecule’s permanent dipole moment upon electronic or vibronic excitation, encapsulating both light–matter coupling and self-energy renormalization effects. This dimensionless factor governs exponential suppression of photonic coupling (“Franck–Condon” reduction), determines the non-perturbative modification of cavity QED coupling strengths, and enables a unified framework to describe spectral progressions, multiphoton emission, and light–matter decoupling phenomena in arbitrary dielectric environments (Wei et al., 2022).

1. Theoretical Definition and Physical Interpretation

In macroscopic quantum electrodynamics (QED), the polaritonic Huang–Rhys factor $S_p$ is analogous to its vibrational counterpart $S_v$, but replaces nuclear coordinate displacements with EM field mode displacements resulting from a change in the permanent molecular dipole ($\Delta \mu = \mu_g - \mu_e$) during an electronic transition. The photon field effectively undergoes a multimode shift in response, analogous to a “polaron” displacement (Wei et al., 2022).

The explicit definition of $S_p$ is:

$$S_p = \frac{1}{\hbar\,\pi\,\epsilon_0\,c^2} \int_0^\infty d\omega\, \Delta\mu \cdot \operatorname{Im} \overline{G}(r_m,r_m,\omega)\cdot\Delta\mu$$

where $\overline{G}(r_m,r_m,\omega)$ is the dyadic Green's function at the molecular position, embodying all effects of geometry, absorption, and dispersion.

Physically, $S_p$ measures the extent to which the EM field vacuum is reorganized by the molecular transition. As $S_p$ increases, the overlap between the original and displaced EM modes decreases, leading to an exponential Franck–Condon–like suppression $S_p$0 in radiative transitions and catalyzing higher-order (“polaritonic progression”) processes.

2. Fundamental Equations and Operational Framework

The core macroscopic QED results defining how $S_p$1 alters light–matter coupling are summarized as follows (Wei et al., 2022):

• Polaritonic coupling shift:

$S_p$2

where $S_p$3 is the transition dipole.

• Modified light–matter coupling:

$S_p$4

$S_p$5 is the bare vacuum-induced coupling:

$S_p$6

These equations are parameter-free once $S_p$7 is computed for the environment. The quantities $S_p$8, $S_p$9, and $S_p$0 can be determined analytically (for simple geometries) or numerically (via FDTD or BEM solvers).

3. Regimes of Polaritonic Progression and Physical Consequences

The value of $S_p$1 controls the nature of photonic spectral features and polaritonic state formation (Wei et al., 2022). The transition probabilities for photon emission follow a Poisson-like distribution:

$S_p$2

key regimes are:

• $S_p$3: Single-photon transitions dominate; $S_p$4.
• $S_p$5: Multipolaritonic transitions (i.e., $S_p$6 emission), enhanced non-radiative transitions.
• Large $S_p$7: Suppression of $S_p$8 transitions (“light–matter decoupling”); strong multiphoton processes analogous to Purcell factor breakdown in deep-strong-coupling regimes.

These progressions directly correspond to the polaron physics and the familiar vibrational Franck–Condon progressions, but now in the photonic sector.

4. Polaritonic Huang–Rhys Factor in Model Hamiltonians

The Holstein–quantum–Rabi (HQR) model and its extensions incorporate the Huang–Rhys factor in the polariton eigenproblem, revealing rich photonic and vibronic structure as a function of $S_p$9 (López et al., 16 Sep 2025). The model Hamiltonian typically contains:

• Molecular harmonic potentials (diabatic states): $S_v$0, $S_v$1 with displacement parameter $S_v$2 (Huang–Rhys parameter), with $S_v$3.
• Photon–molecule coupling: Either a constant or coordinate-dependent term.
• Hybrid basis: The displaced oscillator states $S_v$4 infuse $S_v$5 into polaritonic eigenstates and matrix elements.

Key results include:

• Polariton Rabi splittings at resonance are modulated by $S_v$6 via the factor $S_v$7 (where $S_v$8 is the Laguerre polynomial).
• Dramatic restructuring of the energy spectrum with increasing $S_v$9, including a proliferation of avoided crossings and altered multiphoton dynamics.
• Maximum photon yield after a vibronic–photonic cycle is sensitively dependent on both the magnitude and sign of $\Delta \mu = \mu_g - \mu_e$0.

Non-adiabatic and counter-rotating terms modify the sequence of light-induced crossings, enabling multiple photon generation pathways for positive $\Delta \mu = \mu_g - \mu_e$1 and suppressing them for negative $\Delta \mu = \mu_g - \mu_e$2 (López et al., 16 Sep 2025).

5. Extensions: Cavity and Polariton Branch–Specific Factors

When a vibrating molecule is resonantly coupled to a single cavity mode, the vibrational Huang–Rhys parameter divides into distinct polaritonic factors ($\Delta \mu = \mu_g - \mu_e$3, $\Delta \mu = \mu_g - \mu_e$4) for the upper and lower branches (Rokaj et al., 2023). These are given by:

$\Delta \mu = \mu_g - \mu_e$5

• At resonance ($\Delta \mu = \mu_g - \mu_e$6), $\Delta \mu = \mu_g - \mu_e$7.
• Off-resonance, $\Delta \mu = \mu_g - \mu_e$8 approaches $\Delta \mu = \mu_g - \mu_e$9 while $S_p$0 (or vice versa depending on detuning).
• This tunability enables “Franck–Condon transfer,” with vibrational overlap and transition probabilities distributed and controlled among polaritonic branches.

This framework generalizes to more complex environments and predicts the efficient redistribution of Franck–Condon factors between light–matter hybridized states. Enhanced photochemical processes such as molecular photoassociation have been experimentally linked to such polariton-enhanced factors (Rokaj et al., 2023).

6. Experimental Significance and Quantitative Benchmarks

Application of the macroscopic QED formalism, utilizing experimentally determined or simulated Green’s functions, shows that:

• In typical nanoplasmonic and Fabry–Pérot environments, $S_p$1, leading to minimal Franck–Condon suppression and negligible ($S_p$2–$S_p$3) coupling shift.
• Calculated modified coupling $S_p$4 is within a factor of $S_p$5–$S_p$6 of observed Rabi splittings, indicating quantitative robustness (Wei et al., 2022).

These results confirm the polaritonic Huang–Rhys factor as a quantitative, predictive tool for interpreting and designing cavity QED and polaritonic chemistry experiments, particularly in contexts where permanent dipole changes and field-induced self-energies are non-negligible.

7. Broader Impact and Connections

The polaritonic Huang–Rhys factor unifies multiple strands of quantum optics, molecular photodynamics, and polariton chemistry by offering a rigorous, parameter-free metric for the strength and consequences of permanent-dipole EM field dressing. It enables precise control and prediction of multiphoton emission, non-radiative channels, and light–matter decoupling phenomena, spanning regimes from ultrafast photodynamics to ultracold, cavity-modified chemistry. The control parameters—the Green's function (environment), molecular dipole changes, and light-matter detuning—offer experimentalists and theorists handles to shape photonic and vibronic phenomena at the quantum level (Wei et al., 2022, López et al., 16 Sep 2025, Rokaj et al., 2023).