Case Study 2: WKB and Quantum Tunneling in Chemical Reactions

Overview

Quantum tunneling is not confined to nuclear physics. In chemistry, the transfer of light particles — protons, hydrogen atoms, and even heavier groups — through potential energy barriers is a measurable phenomenon with profound consequences for reaction rates, enzyme catalysis, and the chemistry of the very early universe. The WKB approximation provides the theoretical framework for understanding these effects, bridging quantum mechanics and chemical kinetics.

This case study explores how the WKB tunneling formalism applies to chemical reactions, from simple hydrogen transfer to the catalytic machinery of enzymes, and examines the experimental signatures that distinguish quantum tunneling from classical barrier crossing.

Part 1: Classical Transition State Theory and Its Failures

The Arrhenius Framework

Classical chemical kinetics describes reaction rates through the Arrhenius equation:

$$k(T) = A\,\exp\!\left(-\frac{E_a}{k_BT}\right)$$

where $k(T)$ is the rate constant, $A$ is the pre-exponential (frequency) factor, and $E_a$ is the activation energy — the height of the potential energy barrier separating reactants from products. In transition state theory (TST), the rate is:

$$k_{\text{TST}}(T) = \frac{k_BT}{h}\,\exp\!\left(-\frac{\Delta G^\ddagger}{k_BT}\right)$$

where $\Delta G^\ddagger$ is the free energy of activation.

Both expressions assume that the reacting system must acquire enough thermal energy to surmount the barrier. Particles with energy less than $E_a$ are assumed to be reflected completely. This is a purely classical assumption.

Where Classical Theory Fails

By the mid-20th century, several experimental observations defied classical transition state theory:

1. Anomalous isotope effects. When hydrogen (H) is replaced by deuterium (D) in a reaction, classical theory predicts a rate decrease due to the larger reduced mass affecting the zero-point energy. The predicted "kinetic isotope effect" (KIE) for C-H vs. C-D bond breaking is typically $k_H/k_D \approx 2$–$7$. But some reactions show KIEs of 50 or even 100 — far exceeding the classical prediction.

2. Temperature-independent rates. Some reactions continue at measurable rates even at cryogenic temperatures where $k_BT \ll E_a$. If the only mechanism is thermal activation over the barrier, the rate should be immeasurably small.

3. Curved Arrhenius plots. A plot of $\ln k$ vs. $1/T$ should be a straight line (constant activation energy). Instead, some reactions show downward curvature at low temperatures — the rate decreases more slowly than expected, as though an alternative, temperature-independent pathway opens up.

All three observations are signatures of quantum tunneling.

Part 2: WKB Analysis of Chemical Tunneling

The Chemical Barrier

In a typical chemical reaction, the potential energy surface along the reaction coordinate $x$ has the form:

$$V(x) = E_a \cdot f(x/w)$$

where $E_a$ is the barrier height and $w$ is the barrier width. For hydrogen transfer reactions, typical values are $E_a \sim 0.2$–$1.0$ eV and $w \sim 0.5$–$1.5$ \AA.

The WKB tunneling probability for a particle of mass $m$ and energy $E < E_a$ is:

$$T(E) = \exp\!\left(-\frac{2}{\hbar}\int_{x_1}^{x_2}\sqrt{2m[V(x) - E]}\,dx\right)$$

The Bell Tunneling Correction

In 1959, Ronald Bell introduced the simplest model for incorporating tunneling into chemical kinetics. For a parabolic barrier:

$$V(x) = E_a\left(1 - \frac{x^2}{w^2}\right)$$

the WKB tunneling integral can be evaluated analytically:

$$T(E) = \exp\!\left(-\frac{\pi w}{\hbar}\sqrt{2m(E_a - E)}\right) = \exp\!\left(-\frac{\pi w\sqrt{2mE_a}}{\hbar}\cdot\sqrt{1 - E/E_a}\right)$$

The thermal tunneling correction factor $\kappa$ modifies the classical rate:

$$k = \kappa \cdot k_{\text{TST}}$$

where:

$$\kappa = \frac{\int_0^{\infty}T(E)\,e^{-E/k_BT}\,dE}{\int_{E_a}^{\infty}e^{-E/k_BT}\,dE}$$

For a symmetric parabolic barrier, the Bell tunnel correction factor is:

$$\kappa = \frac{u/2}{\sin(u/2)}, \quad u = \frac{h\nu^\ddagger}{k_BT}$$

where $\nu^\ddagger = (1/2\pi)\sqrt{2E_a/(mw^2)}$ is the imaginary frequency of the inverted barrier.

Estimating Tunneling Importance

The dimensionless parameter that controls the importance of tunneling is:

$$\xi = \frac{w\sqrt{2mE_a}}{\hbar}$$

This is essentially the Gamow exponent evaluated at $E = 0$. For $\xi \gg 1$, tunneling is negligible and classical TST applies. For $\xi \sim 1$, tunneling is significant.

The critical observation: tunneling is most important for the lightest particles. Electrons tunnel readily (which is the basis of all of solid-state physics). Protons tunnel significantly. Heavier atoms tunnel negligibly under normal chemical conditions.

Part 3: The Kinetic Isotope Effect as a Tunneling Diagnostic

Classical Kinetic Isotope Effect

When a hydrogen atom (mass $m_H$) in a reaction is replaced by deuterium ($m_D \approx 2m_H$), the classical kinetic isotope effect (KIE) arises from the difference in zero-point energies:

$$\text{KIE}_{\text{classical}} = \frac{k_H}{k_D} = \exp\!\left(\frac{\Delta\text{ZPE}}{k_BT}\right)$$

For a C-H stretch ($\nu \approx 3000$ cm$^{-1}$) vs. C-D ($\nu \approx 2200$ cm$^{-1}$), $\Delta\text{ZPE} \approx 0.05$ eV, giving $\text{KIE} \approx 7$ at room temperature.

Tunneling-Enhanced KIE

When tunneling is important, the KIE is dramatically amplified. Since $\gamma \propto \sqrt{m}$, the tunneling probability ratio is:

$$\frac{T_H}{T_D} = \exp\!\left(2\gamma_D - 2\gamma_H\right) = \exp\!\left(2\gamma_H(\sqrt{2} - 1)\right) \approx \exp(0.83 \cdot 2\gamma_H)$$

For $2\gamma_H = 10$ (a typical value for proton transfer through a chemical barrier), this gives:

$$\frac{T_H}{T_D} \approx e^{8.3} \approx 4000$$

This is an enormous enhancement — the tunneling contribution to the KIE alone is $\sim 4000$, far exceeding the classical maximum of $\sim 7$.

Experimental Evidence

Alcohol dehydrogenase (ADH) catalyzes the transfer of a hydride (H$^-$) from ethanol to NAD$^+$. The measured KIE at 25$^\circ$C is $k_H/k_D \approx 3.4$ (close to classical). But at $-50^\circ$C (in a cryosolvent), the KIE rises to $\sim 50$, indicating that tunneling contributes significantly at low temperature.

Methylamine dehydrogenase (MADH) shows a primary KIE of $k_H/k_D \approx 17$ at room temperature — well above the classical limit of $\sim 7$. The Arrhenius plot shows curvature, and the activation energies for H and D transfer differ by more than the expected zero-point energy difference. All three signatures point to tunneling.

Soybean lipoxygenase (SLO) shows the most dramatic evidence: $k_H/k_D \approx 80$ at 35$^\circ$C, with an Arrhenius prefactor ratio $A_H/A_D \gg 1$ (classically, this ratio should be $\sim 1$). Theoretical analysis using WKB-based models confirms that proton tunneling dominates this reaction at room temperature.

Part 4: Tunneling in Enzyme Catalysis

The Marcus-Like Model for Proton-Coupled Electron Transfer

Modern theories of enzyme-catalyzed hydrogen transfer combine Marcus theory (for the heavy-atom reorganization) with WKB tunneling (for the hydrogen transfer):

$$k = \int_0^{\infty} P(R)\, T(R, E)\, e^{-E/k_BT}\, dE\, dR$$

Here $R$ is the donor-acceptor distance, $P(R)$ is the probability of finding the donor and acceptor at distance $R$ (determined by the protein dynamics), and $T(R, E)$ is the WKB tunneling probability through the barrier at that distance and energy.

The key insight is that the enzyme controls the donor-acceptor distance through its protein structure and dynamics. By positioning the donor and acceptor close together (reducing $R$, which narrows the barrier), the enzyme exponentially enhances the tunneling rate. This is tunneling-promoting vibration — the protein breathes, occasionally bringing the donor and acceptor close enough for efficient tunneling.

The WKB Barrier Width Is the Control Parameter

For a given barrier height $E_a$, the WKB tunneling exponent scales linearly with barrier width $w$:

$$2\gamma \propto w\sqrt{2mE_a}$$

Reducing $w$ by just 0.1 \AA (a tiny structural change, well within the range of protein flexibility) changes $2\gamma$ by:

$$\Delta(2\gamma) \approx \frac{\sqrt{2m_H \cdot E_a}}{\hbar}\cdot 0.1\,\text{\AA}$$

For $E_a = 0.5$ eV: $\Delta(2\gamma) \approx 1.6$, giving a tunneling enhancement of $e^{1.6} \approx 5$. A 0.1 \AA change in donor-acceptor distance changes the tunneling rate by a factor of 5!

This extreme sensitivity explains why enzyme active sites are so precisely structured: small changes in geometry produce large changes in catalytic rate through the tunneling channel.

Part 5: Tunneling at Extreme Conditions

Interstellar Chemistry

In the cold molecular clouds of interstellar space ($T \sim 10$–$50$ K), thermal energies are negligible compared to most chemical activation barriers. Yet complex molecules — including water, methanol, formaldehyde, and even amino acid precursors — are synthesized on the surfaces of dust grains.

At these temperatures, classical rates for H-atom addition reactions would be immeasurably slow. The reactions proceed by quantum tunneling. The WKB tunneling rate for hydrogen atom addition to CO on an ice surface, for example, has been measured in the laboratory and modeled computationally:

$$k_{\text{tunnel}} \approx \nu_0\, \exp\!\left(-\frac{2}{\hbar}\int\sqrt{2m_H[V(x) - E]}\,dx\right)$$

where $\nu_0 \sim 10^{12}$ Hz is the attempt frequency and the barrier is the activation energy for the H + CO $\to$ HCO reaction ($E_a \approx 0.08$ eV, $w \approx 0.9$ \AA). The calculated rate at 10 K is $k \sim 10^4$ s$^{-1}$ — fast enough to be astrophysically relevant on the $\sim 10^5$-year timescales of molecular cloud lifetimes.

Tunneling in Solid-State Chemistry at 0 K

In the 1970s, Goldanskii predicted that certain chemical reactions should proceed at measurable rates even at absolute zero, driven entirely by tunneling. This was confirmed experimentally by observing the polymerization of formaldehyde (H$_2$C=O) in the solid state at 4 K. The reaction rate was nonzero, temperature-independent, and consistent with WKB tunneling predictions.

Muon-Catalyzed Fusion

A muon ($\mu^-$) is 207 times heavier than an electron but still much lighter than a proton. When a muon replaces an electron in a hydrogen molecule, the bond length shrinks by a factor of $\sim 207$, bringing the two nuclei close enough for the Coulomb barrier to become thin enough for tunneling. Muon-catalyzed fusion occurs at room temperature — a stunning demonstration of WKB tunneling physics.

The tunneling integral scales as $\sqrt{m_\mu}\,R$, where $R$ is the internuclear distance. Since $R \propto 1/m_\mu$ (the muonic "Bohr radius"), the Gamow exponent scales as $\gamma \propto \sqrt{m_\mu}/m_\mu = 1/\sqrt{m_\mu}$ — smaller than for electronic hydrogen, enabling fusion.

Part 6: Computational Methods for Chemical Tunneling

The Wigner Tunneling Correction

The simplest estimate of tunneling in chemical reactions is the Wigner correction (1932):

$$\kappa_W = 1 + \frac{1}{24}\left(\frac{h\nu^\ddagger}{k_BT}\right)^2$$

where $\nu^\ddagger$ is the imaginary frequency of the transition state. This is a second-order expansion of the exact Bell factor and is valid when $h\nu^\ddagger \ll k_BT$ (high temperature limit).

The Eckart Barrier

For more accurate WKB calculations, the Eckart barrier provides an analytically tractable model:

$$V(x) = \frac{V_f\, e^{x/a}}{(1 + e^{x/a})^2} + \frac{(V_f - V_r)e^{x/a}}{1 + e^{x/a}}$$

where $V_f$ and $V_r$ are the forward and reverse barrier heights and $a$ is the width parameter. The exact transmission coefficient for the Eckart barrier is known analytically, providing a benchmark for WKB calculations.

Instanton Theory

The most sophisticated semiclassical approach to chemical tunneling is instanton theory, where the tunneling path is found by solving Newton's equations in imaginary time (i.e., in the inverted potential). The instanton path minimizes the WKB action integral:

$$S_{\text{inst}} = \int \sqrt{2m[V(\mathbf{x}) - E]}\,|d\mathbf{x}|$$

along the optimal tunneling path through the multi-dimensional potential energy surface. This path is not necessarily the minimum energy path — the tunneling path can "cut corners" through the barrier, especially for light particles.

Modern computational chemistry codes (e.g., DL-FIND, Polyrate, i-PI) implement instanton methods to compute tunneling corrections for reactions in full dimensionality, including all atomic degrees of freedom.

Discussion Questions

1. Why is quantum tunneling more important for hydrogen transfer than for heavier atom transfer (e.g., carbon or nitrogen)? Quantify your argument using the WKB tunneling exponent.

2. The kinetic isotope effect is often used as evidence for tunneling. But a large KIE can also arise from classical zero-point energy effects. How can experimentalists distinguish between these two mechanisms? What additional measurements (beyond the KIE) are diagnostic of tunneling?

3. Enzymes have evolved over billions of years to optimize their catalytic rates. If tunneling contributes significantly to proton transfer reactions, should we expect natural selection to have optimized the donor-acceptor distance for tunneling? What structural features of enzyme active sites might reflect this optimization?

4. In interstellar chemistry, tunneling enables reactions at 10 K that would be impossible classically. Does this mean that the chemistry of molecular clouds is fundamentally quantum mechanical? Or is tunneling just a correction to an otherwise classical picture?

5. Muon-catalyzed fusion demonstrates that tunneling can enable nuclear reactions at room temperature. Why is muon-catalyzed fusion not a practical energy source? What limits its efficiency, and how does the WKB framework help us understand those limits?

Quantitative Exercises

E1. For the reaction H + H$_2$ $\to$ H$_2$ + H (hydrogen atom exchange), the barrier height is $E_a = 0.42$ eV and the barrier width is approximately $w = 0.87$ \AA.

(a) Calculate the Gamow exponent $2\gamma$ for a proton at $E = 0$ (the maximum tunneling barrier).

(b) Calculate $2\gamma$ for deuterium under the same conditions.

(c) Compute the tunneling-contribution to the KIE at low temperature (where tunneling dominates): $k_H/k_D \approx T_H/T_D$.

(d) At what temperature does the classical KIE equal 7? At what temperature does tunneling begin to dominate over classical barrier crossing?

E2. A proton transfer reaction in an enzyme has barrier height $E_a = 0.25$ eV and donor-acceptor distance $R = 2.7$ \AA (so the barrier width is $w \approx R - 2r_{\text{vdW}} \approx 0.5$ \AA, where $r_{\text{vdW}} \approx 1.1$ \AA for O or N).

(a) Estimate the tunneling probability for a proton at the zero-point level.

(b) If protein breathing reduces $R$ by 0.2 \AA (to 2.5 \AA), by what factor does the tunneling probability increase?

(c) Estimate the enzymatic rate enhancement due to tunneling at the optimal donor-acceptor distance, compared to a "frozen" enzyme with $R = 3.0$ \AA.

E3. Calculate the WKB tunneling rate for the reaction H + CO $\to$ HCO on an interstellar dust grain at $T = 10$ K. Use $E_a = 0.08$ eV, $w = 0.9$ \AA, $m = m_H$, and attempt frequency $\nu_0 = 10^{12}$ Hz. Compare with the classical Arrhenius rate $k_{\text{class}} = \nu_0 \exp(-E_a/k_BT)$ at the same temperature.