Gas-Phase Hydrogen Exchange Reaction (HER)

• Gas-phase HER is a fundamental isotopic exchange process with a defined activation barrier on the H₃ potential energy surface.
• Quantum scattering and DFT analyses reveal detailed state-to-state dynamics, including tunneling and resonance phenomena, that govern the reaction rate.
• Cavity-induced modifications enable hybridization between transition-state and dynamical-barrier resonances, yielding measurable shifts in reactive decay pathways.

Gas-phase hydrogen exchange reaction (HER) commonly denotes the elementary isotopic exchange $D + H_2 \rightarrow HD + H$, a prototypical reactive process on the $H_3$ potential energy surface (PES) and its isotopologues. In the treatment centered on the DHH isotopologue, the system functions as a canonical benchmark for quantum reaction dynamics because it is the smallest neutral reactive system with an energy barrier, admits highly accurate PES and quantum scattering descriptions, and exhibits tunneling, shape/Feshbach resonances, rich state-to-state dynamics, and well-defined $S$-matrix poles. In free space the process follows the conventional exchange channel, whereas in a resonant dark cavity it can additionally proceed through the radiative pathway $D + H_2 \rightarrow HD + H + E_{photon}$ with $E_{photon}=\hbar\omega_{cav}$, so that vacuum-field hybridization alters the decay structure of reactive resonances (Berenstein et al., 13 Jul 2025).

1. Benchmark character and elementary physical content

The gas-phase HER studied in the cited cavity-QED work is the isotopic hydrogen exchange reaction $D + H_2 \rightarrow HD + H$. It is nearly thermoneutral, with small zero-point-energy differences that can make the reaction slightly exo- or endothermic depending on rovibrational state, but the dominant dynamical feature is the activation barrier on the $H_3$ PES near the saddle point. The reaction is therefore a stringent test case for any framework that aims to connect microscopic resonance structure with measurable reaction rates (Berenstein et al., 13 Jul 2025).

Its benchmark status follows from several convergent features. The system has minimal nuclear complexity yet nontrivial reactivity; it supports highly accurate quantum-scattering treatments; and it displays phenomena that are central to modern reaction dynamics, including tunneling, resonance structure, and state-resolved channel competition. Within the cavity-control program, these properties make HER unusually suitable for a falsifiable comparison between theory and experiment, because the cavity modifies a reaction whose underlying free-space dynamics are already sharply constrained.

A central conceptual point is that the reactive flux is not treated as a featureless barrier crossing. Instead, it is organized by resonant states of the scattering matrix. In free space, the top-of-barrier resonance acts as the transition-state resonance (TSR), and its decay width provides the microcanonical rate input. Inside a cavity, this resonance can hybridize with a lower-lying resonance, opening a radiative decay pathway and thereby changing the reaction rate without external pumping.

2. Free-space dynamics and reaction-path formulation

The free-space treatment constructs the DHH PES for the collinear reaction by DFT at the mPW2PLYP/6-311++G(d,p) level using Gaussian 16, scans the Intrinsic Reaction Coordinate (IRC), and evaluates the mass-weighted Hessian to obtain normal-mode frequencies along the reaction coordinate $X$ (Berenstein et al., 13 Jul 2025). On that basis, the adiabatic reaction-path Hamiltonian is reduced to a one-dimensional potential,

$V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$

where $V_{rc}(X)$ is the electronic PES along the reaction coordinate and $H_3$0 is the sum of the bound perpendicular-mode frequencies, comprising the symmetric stretch and two degenerate bends.

The resulting $H_3$1 is given a smooth analytic representation by the Schlessinger continued-fraction method. This step is technically important because the resonance analysis requires a stable analytic continuation of the barrier region and its flanking wells. The framework separates motion along $H_3$2 from the bound perpendicular modes, but does not treat those perpendicular modes as dynamically irrelevant: their variation with $H_3$3 is precisely what generates the dynamical-barrier structure that underlies the cavity-active resonance pair.

Resonances are obtained by uniform complex scaling. In one-dimensional form, the resonance states satisfy

$H_3$4

with

$H_3$5

After the transformation $H_3$6, the outgoing-wave resonance eigenfunctions become square-integrable, and the resonance poles are exposed as $H_3$7-independent for sufficiently large scaling angle; in the D + $H_3$8 case, $H_3$9 was used. The top-of-barrier pole is identified as the TSR, and its width $S$0 sets the free-space microcanonical rate.

3. Resonance picture of cross sections and rates

The free-space reaction rate is framed through the TSR width,

$S$1

This is the microcanonical quantity that anchors the later cavity analysis. Standard thermal rate constants are then obtained by the conventional energy average of reactivity over a Maxwell-Boltzmann distribution. For state-to-state scattering, the integral cross section is

$S$2

where $S$3 and $S$4 are $S$5-matrix probabilities. The thermal rate constant is

$S$6

Within the non-Hermitian resonance description, $S$7 and hence $S$8 are dominated near the barrier top by the TSR. The free-space theory in the cited work emphasizes relative cavity-induced changes rather than tabulating an exhaustive set of thermal rates, but the formal bridge to experiment is explicit: the microscopic decay width from the TSR supplies the resonance input for the thermal observable.

A second resonance family is equally important for the cavity problem. Because the perpendicular normal-mode frequencies vary significantly along the reaction coordinate, the adiabatic potential develops a dynamical barrier and a ladder of lower-lying dynamical-barrier (DB) resonances. These DB states are not merely spectators. They provide the lower-energy partner that can couple resonantly to the TSR when the cavity mode is tuned to the corresponding energy gap.

4. Dark-cavity formulation and polaritonic mechanism

In the cavity-QED treatment, “dark cavity” means that no external driving field is applied; all rate modifications arise from the vacuum field and light-matter hybridization. The full Hamiltonian is written schematically as

$S$9

with $D + H_2 \rightarrow HD + H + E_{photon}$0 and $D + H_2 \rightarrow HD + H + E_{photon}$1 the annihilation and creation operators of the selected cavity mode. In the length gauge,

$D + H_2 \rightarrow HD + H + E_{photon}$2

whereas in the acceleration gauge convenient for scattering,

$D + H_2 \rightarrow HD + H + E_{photon}$3

The relevant effective matrix element is

$D + H_2 \rightarrow HD + H + E_{photon}$4

Restricting the problem to the dominant resonance pair $D + H_2 \rightarrow HD + H + E_{photon}$5 and the lowest photon-number subspace $D + H_2 \rightarrow HD + H + E_{photon}$6 yields a $D + H_2 \rightarrow HD + H + E_{photon}$7 non-Hermitian polaritonic Hamiltonian,

$D + H_2 \rightarrow HD + H + E_{photon}$8

On resonance, the vacuum-emission condition is

$D + H_2 \rightarrow HD + H + E_{photon}$9

Under that condition, the initially populated $E_{photon}=\hbar\omega_{cav}$0 state mixes strongly with $E_{photon}=\hbar\omega_{cav}$1, and the cavity opens the additional radiative reaction channel

$E_{photon}=\hbar\omega_{cav}$2

The physical mechanism has two indispensable ingredients. First, the reaction path must possess a single dominant saddle region so that a well-defined TSR is initially populated. Second, the perpendicular-mode frequencies must vary enough along $E_{photon}=\hbar\omega_{cav}$3 to generate the lower-lying DB resonances. The cavity does not create reactivity from nothing; it reroutes reactive decay by hybridizing two specific $E_{photon}=\hbar\omega_{cav}$4-matrix poles.

5. Quantitative predictions, collective coupling, and isotopic symmetry breaking

For an initial $E_{photon}=\hbar\omega_{cav}$5 preparation, diagonalization of $E_{photon}=\hbar\omega_{cav}$6 gives polaritonic eigenvalues $E_{photon}=\hbar\omega_{cav}$7 and the cavity-modified rate

$E_{photon}=\hbar\omega_{cav}$8

with polaritons

$E_{photon}=\hbar\omega_{cav}$9

Near the barrier, the enhancement factor is

$D + H_2 \rightarrow HD + H$0

The single-photon field amplitude scales as $D + H_2 \rightarrow HD + H$1, the effective coupling obeys $D + H_2 \rightarrow HD + H$2, and for $D + H_2 \rightarrow HD + H$3 identical noninteracting molecules coupled to the same mode the collective coupling follows the Tavis-Cummings scaling $D + H_2 \rightarrow HD + H$4 (Berenstein et al., 13 Jul 2025).

The quantitative predictions are strongly parameter dependent:

Regime	Coupling condition	Predicted $D + H_2 \rightarrow HD + H$5
Many-molecule limit $D + H_2 \rightarrow HD + H$6	$D + H_2 \rightarrow HD + H$7–$D + H_2 \rightarrow HD + H$8 a.u.	typical $D + H_2 \rightarrow HD + H$9
Many-molecule limit $H_3$0	$H_3$1 a.u.	$H_3$2
Many-molecule limit $H_3$3	$H_3$4 a.u.	$H_3$5
Single/several molecules $H_3$6–$H_3$7	same $H_3$8 range	enhancement without attenuation

These results establish that cavity control is not generically accelerative. In the many-molecule limit, the same dark-cavity mechanism can enhance or suppress the rate, depending on the coupling value. The strongest changes occur when $H_3$9 is tuned to the TS-DB gap and the coupling is large enough to compete with the resonance-width mismatch and the cavity loss rate $X$0, but not so large that the system moves beyond the optimal mixing regime typical of non-Hermitian two-level problems. A plausible implication is that cavity design must target a particular resonance pair rather than maximize coupling indiscriminately.

The microscopic origin of light-matter coupling is isotopic symmetry breaking. Homonuclear $X$1 has no permanent dipole, but along the D + $X$2 reaction path the transient DHH complex is mass-asymmetric and vibrationally anharmonic, yielding a nonzero resonance-to-resonance transition amplitude,

$X$3

or equivalently, in the acceleration gauge, $X$4. The small H/D mass imbalance is therefore sufficient to break the symmetry that would suppress coupling in a strictly homonuclear setting.

6. Experimental observables, methodological limits, and nomenclature

The experimentally decisive observables are the HD product yield or state distribution in and out of the cavity, together with correlated photon emission at $X$5. The proposed radiative pathway $X$6 is described as the “smoking gun” for cavity-induced dynamics. Suggested platforms include merged molecular beams with velocity-map imaging and cryogenic buffer-gas cells integrated with infrared cavities. The latter are described as capable of delivering internally cold molecules at densities $X$7 and low velocities, with demonstrated trapping potentials $X$8 using $X$9 mode waists and $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$0 IR lasers; larger mode volumes with moderate power scaling are stated to be feasible. Expected populations exceed $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$1 trapped species in the cavity mode volume, which would permit exploration of collective coupling through $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$2 (Berenstein et al., 13 Jul 2025).

The proposed measurement protocol is correspondingly concrete: tune $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$3 to the computed TS-DB gap, compare HD formation with and without the cavity, and count cavity photons simultaneously. Photon detection can then trigger time-of-flight mass spectrometry to isolate the cavity-induced pathway from background reactions. The relevant cavity design targets are a sufficiently small mode volume to ensure noticeable $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$4 and a sufficiently high $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$5 factor, equivalently small $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$6, so that hybridization competes with $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$7 and $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$8.

The theoretical framework is deliberately restrictive. It assumes a single-mode cavity and the long-wavelength dipole approximation; employs a nuclear adiabatic approximation along the reaction coordinate; constructs $V_{ad}(X)=V_{rc}(X)+\frac{\hbar}{2}\Omega(X),$9 from a collinear DFT-based PES rather than the most elaborate ab initio surface; uses a non-Hermitian complex-scaling treatment to extract rates from resonance widths; and treats molecules as noninteracting except for their common coupling to the cavity mode. These assumptions do not negate the benchmark value of the system, but they define the scope within which the predicted enhancement and attenuation factors should be interpreted.

A recurrent source of confusion is nomenclature. In the gas-phase reaction-dynamics literature discussed here, HER denotes hydrogen exchange reaction. In electrochemistry, HER commonly denotes hydrogen evolution reaction, as in studies of bubble trapping and transport in porous electrodes during water-splitting conditions (Ferguson et al., 17 Dec 2025). The two usages are unrelated in mechanism: the gas-phase HER concerns exchange dynamics on the $V_{rc}(X)$0 PES and cavity-induced modification of reactive resonances, whereas the electrochemical usage concerns interfacial proton reduction, gas evolution, and multiphase transport.