Module 7 · Quantum Effects
Quantum Tunneling & Kinetic Isotope Effects
Hydrogen-transfer reactions in enzymes and constrained organic molecules are governed not only by classical over-barrier passage but by quantum mechanical tunneling through the barrier. The signature of tunneling is the kinetic isotope effect (KIE): the ratio of the rate constants for protium (H) and deuterium (D) transfer.
7.1 Transition State Theory and Tunneling Corrections
The classical Eyring rate constant is:
\[ k_{\text{TST}}(T) = \frac{k_BT}{h} \exp\!\left(-\frac{\Delta G^\ddagger}{k_BT}\right) \tag{7.1}\]
Tunneling corrections multiply the classical rate by a factor \(\kappa(T) \geq 1\):
\[ k(T) = \kappa(T) \cdot k_{\text{TST}}(T) \tag{7.2}\]
The Bell correction (the simplest quantum correction) treats the barrier as an inverted parabola with imaginary frequency \(\omega^\ddagger\) at the transition state:
\[ \kappa_{\text{Bell}}(T) = \frac{u/2}{\sin(u/2)}, \qquad u = \frac{\hbar |\omega^\ddagger|}{k_BT} \tag{7.3}\]
The Wigner correction is the leading-order quantum correction in \(\hbar^2\):
\[ \kappa_{\text{Wigner}}(T) = 1 + \frac{1}{24}\!\left(\frac{\hbar\omega^\ddagger}{k_BT}\right)^2 + O(\hbar^4) \tag{7.4}\]
7.2 The KIE and Its Physical Content
\[ \text{KIE}(T) = \frac{k_H}{k_D} = \frac{\kappa_H}{\kappa_D} \cdot \exp\!\left(-\frac{\Delta G^\ddagger_H - \Delta G^\ddagger_D}{k_BT}\right) \tag{7.5}\]
The \(\Delta\Delta G^\ddagger\) term reflects the difference in zero-point energies (ZPE) between H and D in the reactant and transition states. For a classical transfer, \(\text{KIE}_{\text{classical}} \approx \exp(\Delta ZPE_{H-D}/k_BT) \approx 2\)–3 at 300 K. Larger KIEs (4–10) indicate significant tunneling contributions through \(\kappa_H / \kappa_D > 1\).
7.3 Tunneling Suppression by Geometric Confinement
Tunneling probability depends exponentially on the barrier width sampled by the transferring nucleus. In systems where the donor–acceptor distance \(d_{D\cdots A}\) fluctuates freely, the transfer occasionally samples compressed geometries where the effective barrier width is thin — enhancing tunneling. In conformationally rigid scaffolds, \(\langle(\delta d_{D\cdots A})^2\rangle\) is small: the system rarely accesses the compressed geometry, and the effective tunneling path is longer on average.
\[ P_{\text{tunnel}} \propto \exp\!\left(-\frac{2}{\hbar}\int_{x_1}^{x_2}\sqrt{2m\bigl(V(x) - E\bigr)}\,dx\right) \tag{7.6}\]
The integral is the WKB tunneling integral over the classically forbidden region \([x_1, x_2]\). Scaffold rigidity increases the effective width \(x_2 - x_1\) (by reducing sampling of compressed \(d_{D\cdots A}\) geometries) and therefore suppresses \(P_{\text{tunnel}}\) — driving \(\kappa \to 1\) and \(\text{KIE} \to \text{KIE}_{\text{classical}} \approx 2\)–3.
Foundational Experimental Observation
This mechanism — suppression of tunneling by geometric confinement in constrained bicyclic scaffolds, despite thermodynamic favourability of the transfer — was first observed in the study of hydrogen transfer in conformationally restricted bicyclic systems (IFP, 1985). The observation anticipated modern quantum biology by nearly two decades and is the experimental anchor for the theoretical framework developed in this course.
7.4 Benchmark KIE Values Across Scaffold Geometries
| System | \(\langle\delta d^2\rangle^{1/2}\) (Å) | KIE (RPMD/DFT) | KIE (MAML-NequIP) | KIE (Classical) | Tunneling regime |
|---|---|---|---|---|---|
| M0 — Malonaldehyde | 0.18 | 4.5 | 4.3 | 2.1 | strong |
| C1 — Cyclohexadienyl | 0.12 | 3.8 | 3.6 | 1.9 | moderate |
| B2 — Norbornyl | 0.08 | 2.9 | 2.8 | 1.6 | partial |
| B3 — Bicyclo[2.2.2] | 0.05 | 2.1 | 2.0 | 1.4 | weak |
| IFP — Constrained bicyclic | 0.03 | 1.5 | 1.4 | 1.1 | suppressed |