Activated complex theory: the Eyring equation and the Arrhenius prefactor

Intuition Beginner

Every chemical reaction requires the reactants to climb over an energy barrier. The arrangement of atoms at the top of this barrier is called the transition state or activated complex. It is not a stable molecule: it is the highest-energy configuration along the pathway from reactants to products, and it exists for roughly one molecular vibration period before collapsing forward to products or backward to reactants.

The Arrhenius equation $k=A{e}^{-E_a/RT}$ fits experimental rate constants with two parameters: the activation energy $E_a$ and the pre-exponential factor $A$. The exponential term accounts for the fraction of molecules energetic enough to surmount the barrier. But what determines $A$? Transition-state theory answers this question from first principles using statistical mechanics.

The central result is the Eyring equation:

$$k=\frac{k_{B}T}{h}{e}^{\Delta S^{‡}/R}{e}^{-\Delta H^{‡}/RT}$$

The universal prefactor $k_{B}T/h\approx 6.2\times 10^{12}{\text{s}}^{-1}$ at 298 K replaces the empirical $A$. The activation enthalpy $\Delta H^{‡}$ replaces $E_a$, and the activation entropy $\Delta S^{‡}$ encodes how ordered or disordered the transition state is compared with the reactants. A tight, constrained transition state has $\Delta S^{‡}<0$, reducing $A$; a loose, floppy transition state has $\Delta S^{‡}>0$, increasing $A$.

Visual Beginner

  Potential energy surface for A + BC -> AB + C

  Energy
    |
    |           (AB-C)ddagger
    |            /  \
    |           /    \
    |          /      \
    |   A+BC /        \ AB+C
    |        /          \
    |_______/            \________
    |
    +----- Reaction coordinate ----->

    Delta-H-ddagger = height of saddle above reactants
    Delta-S-ddagger = entropy change at the saddle
    k = (kB T / h) exp(DSdd/R) exp(-DHdd/RT)

  Arrhenius vs Eyring:
    Arrhenius:  k = A exp(-Ea / RT)       [empirical]
    Eyring:     k = (kBT/h) exp(DSdd/R) exp(-DHdd / RT)  [derived]

    A = (e * kBT / h) exp(DSdd / R)      [for unimolecular reactions]

The transition state sits at the saddle point of the potential energy surface: a maximum along the reaction coordinate but a minimum in all other directions. Think of a mountain pass. You climb up to it (the reaction coordinate), but the terrain drops away on either side (the bound vibrations). The partition function of the transition state counts the states available in this pass.

Worked example Beginner

Problem: The gas-phase decomposition of cyclopropane to propene has $k=2.7\times 10^{13}{e}^{-272000/(RT)}{\text{s}}^{-1}$. Use the Eyring equation to estimate the entropy of activation at 500 K.

Solution:

The Arrhenius form gives $A=2.7\times 10^{13}{\text{s}}^{-1}$. For a unimolecular gas-phase reaction, the Eyring equation gives:

$$A=e\cdot \frac{k_{B}T}{h}{e}^{\Delta S^{‡}/R}$$

At $T=500$ K: $k_{B}T/h=1.044\times 10^{13}{\text{s}}^{-1}$.

$$2.7\times 10^{13}=e\times 1.044\times 10^{13}\times e^{\Delta S^{‡}/R}$$ $$\frac{2.7\times 10^{13}}{2.718\times 1.044\times 10^{13}}=e^{\Delta S^{‡}/R}$$ $$0.951=e^{\Delta S^{‡}/R}$$ $$\Delta S^{‡}/R=\ln (0.951)=-0.050$$ $$\Delta S^{‡}=-0.050\times 8.314=-0.4\text{J/(mol K)}$$

The near-zero entropy of activation is consistent with a unimolecular isomerisation: the transition state has roughly the same number of accessible configurations as the reactant. The ring-opening does not dramatically constrain or liberate the molecule.

Check your understanding Beginner

Formal definition Intermediate+

Potential energy surfaces and the transition state

Consider a reaction $A+B\to P$. The electronic energy of the combined system of nuclei and electrons, computed at the Born-Oppenheimer level, is a function $E(\mathbf{R})$ of the $3N$ nuclear coordinates $\mathbf{R}$. This function is the potential energy surface (PES). For a nonlinear system with $N$ nuclei, the PES is a surface in $3N-6$ dimensions (subtracting translation and rotation).

A saddle point on the PES is a stationary point ($\nabla E=0$) where the Hessian matrix of second derivatives has exactly one negative eigenvalue. The eigenvector corresponding to this negative eigenvalue is the reaction coordinate; the remaining $3N-7$ positive eigenvalues correspond to bound vibrations. The configuration at the saddle is the transition state $‡$.

The Eyring equation from partition functions

Transition-state theory (TST) models the reaction as a quasi-equilibrium between reactants and the activated complex:

$$A+B\rightleftharpoons (AB{)}^{‡}\to P$$

The equilibrium constant for forming the activated complex is

$$K^{‡}=\frac{q^{‡}/V}{(q_A/V)(q_B/V)}{e}^{-\Delta E_0/k_{B}T},$$

where $q^{‡}$, $q_A$, $q_B$ are molecular partition functions and $\Delta E_0$ is the zero-point-corrected energy difference between the transition state and the reactants. The partition function of the transition state is factorised by separating the reaction-coordinate vibration. The activated complex has $3N-6$ vibrational modes (for a nonlinear system), but one of these is the unstable mode along the reaction coordinate. Its frequency ${\nu }^{‡}$ is imaginary (the force constant is negative), but formally it contributes a vibrational partition function factor

$$q_{\text{r.c.}}^{‡}=\frac{1}{1-e^{-h{\nu }^{‡}/k_{B}T}}\approx \frac{k_{B}T}{h{\nu }^{‡}},$$

where the last step takes $h{\nu }^{‡}\ll k_{B}T$ (the barrier is broad compared to the quantum level spacing along the reaction coordinate). Define $q^{‡\prime }$ as the partition function of the transition state with the reaction-coordinate mode removed:

$$q^{‡}=q^{‡\prime }\cdot \frac{k_{B}T}{h{\nu }^{‡}}.$$

The rate of passage over the barrier is the frequency ${\nu }^{‡}$ multiplied by the concentration of activated complexes:

$$k={\nu }^{‡}\cdot K^{‡}={\nu }^{‡}\cdot \frac{q^{‡\prime }}{q_A{q}_B/V}\cdot \frac{k_{B}T}{h{\nu }^{‡}}\cdot e^{-\Delta E_0/k_{B}T}.$$

The unknown frequency ${\nu }^{‡}$ cancels:

$$\boxed{k=\frac{k_{B}T}{h}\cdot \frac{q^{‡\prime }/V}{(q_A/V)(q_B/V)}\cdot e^{-\Delta E_0/k_{B}T}.}$$

This is the Eyring equation in its partition-function form. The universal prefactor $k_{B}T/h$ replaces the empirical Arrhenius $A$ with a quantity computable from fundamental constants and temperature alone.

Thermodynamic formulation

Writing $K^{‡\prime }=\exp (-\Delta G^{‡}/RT)$ where $\Delta G^{‡}=\Delta H^{‡}-T\Delta S^{‡}$ is the Gibbs energy of activation (excluding the reaction-coordinate mode), the Eyring equation becomes

$$k=\frac{k_{B}T}{h}\exp \left(\frac{\Delta S^{‡}}{R}\right)\exp \left(-\frac{\Delta H^{‡}}{RT}\right).$$

A transmission coefficient $\kappa$ (close to 1 for direct reactions) may be prefixed:

$$k=\kappa \cdot \frac{k_{B}T}{h}\exp \left(\frac{\Delta S^{‡}}{R}\right)\exp \left(-\frac{\Delta H^{‡}}{RT}\right).$$

Plotting $\ln (k/T)$ vs. $1/T$ (an Eyring plot) gives a straight line:

$$\ln \frac{k}{T}=\ln \frac{k_B\kappa }{h}+\frac{\Delta S^{‡}}{R}-\frac{\Delta H^{‡}}{R}\cdot \frac{1}{T}.$$

The slope yields $\Delta H^{‡}$ and the intercept yields $\Delta S^{‡}$.

Connection to the Arrhenius equation

Comparing the Eyring and Arrhenius forms $k=A{e}^{-E_a/RT}$:

$$E_a=\Delta H^{‡}+nRT,A=e^n\cdot \frac{k_{B}T}{h}\exp \left(\frac{\Delta S^{‡}}{R}\right),$$

where $n=1$ for unimolecular gas-phase reactions and $n=2$ for bimolecular gas-phase reactions. For solution-phase reactions, $n=0$ at constant pressure because the $PV$ work of bringing molecules together is absorbed into the enthalpy. The factor $e^n$ accounts for the difference between $\Delta H^{‡}$ and the temperature-dependent part of the Eyring prefactor.

Key derivation Intermediate+

Partition-function analysis of the activated complex

The partition function of the activated complex $q^{‡\prime }$ factorises in the same way as for a stable molecule, but with one fewer vibrational mode:

$$q^{‡\prime }=q_{\text{trans}}^{‡}\cdot q_{\text{rot}}^{‡}\cdot q_{\text{vib}}^{‡\prime }\cdot q_{\text{elec}}^{‡}.$$

Here $q_{\text{vib}}^{‡\prime }$ is the product over the $3N-7$ bound vibrational modes of the activated complex (for a nonlinear system; $3N-6$ for a linear $\to$ nonlinear transition). The structure and vibrational frequencies of the transition state are not directly measurable but can be computed from quantum chemistry (saddle-point optimisation followed by frequency analysis). The transition state has all real vibrational frequencies except one imaginary frequency, which corresponds to the reaction coordinate.

For a bimolecular gas-phase reaction $A+B\to P$, the rate constant is:

$$k_2=\frac{k_{B}T}{h}\frac{(q_{\text{trans}}^{‡}/V)(q_{\text{rot}}^{‡})(q_{\text{vib}}^{‡\prime })(q_{\text{elec}}^{‡})}{(q_{\text{trans}}^A/V)(q_{\text{rot}}^A)(q_{\text{vib}}^A)(q_{\text{elec}}^A)(q_{\text{trans}}^B/V)(q_{\text{rot}}^B)(q_{\text{vib}}^B)(q_{\text{elec}}^B)}{e}^{-\Delta E_0/k_{B}T}.$$

The ratio of translational partition functions per unit volume gives $q_{\text{trans}}^{‡}/V$ divided by $(q_{\text{trans}}^A/V)(q_{\text{trans}}^B/V)$. Since the activated complex has mass $m_{AB}=m_A+m_B$ and the translational partition function scales as $m^{3/2}$, this ratio introduces a factor proportional to $(2\pi \mu k_{B}T/h^2{)}^{-3/2}$ where $\mu$ is the reduced mass. The result is that $k_2$ has units of $\text{volume}\cdot {\text{time}}^{-1}$ (e.g. ${\text{M}}^{-1}{\text{s}}^{-1}$), consistent with a second-order rate constant.

Entropy of activation and mechanistic diagnosis

The entropy of activation is the single most diagnostic quantity in mechanism analysis:

Reaction type	$\Delta S^{‡}$ (J/(mol K))	$A$ (typical)
Unimolecular dissociation	+10 to +50	$10^{14}$--$10^{16}{\text{s}}^{-1}$
Unimolecular isomerisation	-5 to +5	$10^{12}$--$10^{13}{\text{s}}^{-1}$
Bimolecular association (gas)	-60 to -120	$10^7$--$10^{11}{\text{M}}^{-1}{\text{s}}^{-1}$
Bimolecular (tight cyclic TS)	-120 to -170	$10^3$--$10^7{\text{M}}^{-1}{\text{s}}^{-1}$

A dissociation increases the number of accessible configurations, giving positive $\Delta S^{‡}$. An association decreases translational entropy (two molecules become one), giving large negative $\Delta S^{‡}$. A tight cyclic transition state (e.g., a concerted pericyclic reaction) further restricts internal rotations, driving $\Delta S^{‡}$ even more negative.

The isothermal-isobaric TST formulation

For reactions in solution, the constant-pressure formulation is more natural. The Gibbs energy of activation is $\Delta G^{‡}=-RT\ln K^{‡\prime }$, and the rate constant is:

$$k=\kappa \frac{k_{B}T}{h}\exp \left(-\frac{\Delta G^{‡}}{RT}\right)=\kappa \frac{k_{B}T}{h}\exp \left(\frac{\Delta S^{‡}}{R}\right)\exp \left(-\frac{\Delta H^{‡}}{RT}\right).$$

In solution, $\Delta H^{‡}$ includes contributions from solvation changes during activation. Solvent reorganisation around the activated complex can add or subtract from the gas-phase $\Delta H^{‡}$ by tens of kJ/mol. This is the physical basis for solvent effects on reaction rates: polar solvents stabilise polar transition states, lowering $\Delta H^{‡}$ and accelerating polar reactions relative to nonpolar ones.

Exercises Intermediate

Tunneling corrections, kinetic isotope effects, and limitations of conventional TST Master

Quantum tunneling through the barrier

TST assumes that trajectories cross the dividing surface only once and do not recross. It also treats the reaction coordinate classically. Both assumptions fail for reactions involving light atoms. The reaction coordinate is a vibrational mode with an imaginary frequency, and the quantum-mechanical treatment of motion along this coordinate allows tunneling through the barrier when the de Broglie wavelength of the transferring particle is comparable to the barrier width.

The Wigner tunneling correction (Wigner 1932, Z. Phys. Chem. B 19) is the leading-order quantum correction to the TST rate:

$${\kappa }_{\text{Wigner}}=1+\frac{1}{24}{\left(\frac{h|{\nu }^{‡}|}{k_{B}T}\right)}^2,$$

where ${\nu }^{‡}$ is the magnitude of the imaginary frequency at the saddle point. This correction is small at high temperature but becomes significant when $h|{\nu }^{‡}|\gtrsim k_{B}T$. For hydrogen-atom transfer reactions with $|{\nu }^{‡}|\approx 2000i{\text{cm}}^{-1}$, $h|{\nu }^{‡}|/k_{B}T\approx 10$ at 300 K, giving ${\kappa }_{\text{Wigner}}\approx 5$.

More accurate tunneling corrections include:

Eckart barrier. The potential along the reaction coordinate is modelled by the Eckart function $V(x)=V_0{\text{sech}}^2(x/a)$, for which the quantum transmission probability has an analytic form. The Eckart barrier is parameterised by the forward and reverse barrier heights and the imaginary frequency, and the correction ${\kappa }_{\text{Eckart}}$ is computed by numerically integrating the transmission probability weighted by the Boltzmann distribution of reactant energies. The Eckart correction typically gives $\kappa =2$--$20$ for H-transfer reactions at room temperature.

Small-curvature tunneling (SCT) and large-curvature tunneling (LCT). These methods, developed by Truhlar and coworkers, account for the curvature of the reaction path. The tunneling particle does not follow the classical minimum-energy path but takes a shortcut through the corner-cutting region of the PES. SCT assumes the tunneling follows the concave side of the reaction path; LCT assumes a straight-line path through the PES between reactant and product valleys. The multi-dimensional tunneling corrections (microcanonical optimized tunneling, $\mu$OMT) combine both and are the standard in modern variational TST calculations.

Instanton theory. The semiclassical instanton approach (Langer 1969; Miller 1975) treats tunneling as motion along a periodic orbit in imaginary time (the "instanton"). The instanton is a saddle point of the Euclidean action and corresponds to the most probable tunneling path. The instanton rate is:

$$k_{\text{inst}}=\frac{k_{B}T}{h}\sqrt{\frac{S_0}{2\pi \hslash }}{e}^{-S_0/\hslash },$$

where $S_0$ is the Euclidean action of the instanton orbit. The instanton automatically finds the optimal tunneling path, including corner-cutting and multi-dimensional effects. It is the preferred method for deep tunneling ($\kappa >100$) and becomes exact in the low-temperature limit. At high temperature, the instanton collapses to the classical saddle point and reproduces conventional TST.

Kinetic isotope effects

The most direct experimental test of TST (and tunneling) is the primary kinetic isotope effect (KIE): the ratio of rate constants for isotopologues. Replacing hydrogen with deuterium in a bond-breaking reaction changes the rate through three mechanisms:

1. Zero-point energy effect. The C-H stretching frequency ($\sim 3000{\text{cm}}^{-1}$) drops to $\sim 2200{\text{cm}}^{-1}$ for C-D because $m_{\text{D}}=2m_{\text{H}}$. The zero-point energy $\frac{1}{2}h\nu$ is lower for C-D. If the C-H(D) stretch is the reaction coordinate, the zero-point energy disappears at the transition state and does not affect the barrier. But the other vibrational modes of the transition state involving the transferring H/D retain their isotope dependence. The net effect is that the activation energy for the deuterated reaction is typically 1--5 kJ/mol higher, giving $k_{\text{H}}/k_{\text{D}}\approx 3$--$8$ at room temperature from zero-point effects alone.

2. Tunneling contribution. Because the de Broglie wavelength scales as $\lambda \propto m^{-1/2}$, deuterium tunnels less effectively than hydrogen. The tunneling contribution to the KIE adds a multiplicative factor to the classical ratio. For reactions with significant tunneling, $k_{\text{H}}/k_{\text{D}}$ can reach 50--100 at room temperature, far exceeding the maximum classical prediction of about 10. The Swain-Schaad relationship $k_{\text{H}}/k_{\text{T}}=(k_{\text{H}}/k_{\text{D}}{)}^{1.44}$ provides a diagnostic for tunneling: deviations from this exponent indicate non-classical behavior.

3. Partition-function ratio. The full TST expression for the KIE is

$$\frac{k_{\text{H}}}{k_{\text{D}}}=\frac{(q^{‡\prime }/q_{\text{react}}{)}_{\text{H}}}{(q^{‡\prime }/q_{\text{react}}{)}_{\text{D}}}\cdot e^{-(\Delta E_0^{\text{H}}-\Delta E_0^{\text{D}})/k_{B}T}\cdot \frac{{\kappa }_{\text{H}}}{{\kappa }_{\text{D}}},$$

where the first term is the ratio of partition-function products, the second is the zero-point-energy correction, and the third is the tunneling correction ratio. Each term can be computed independently if the transition-state structure and frequencies are known from quantum chemistry.

Variational transition-state theory

Conventional TST places the dividing surface at the saddle point of the PES. This choice is optimal only if every trajectory that crosses the surface proceeds to products without recrossing. In reality, some trajectories recross the surface, leading TST to overestimate the rate. Variational TST (VTST, Garrett and Truhlar 1979) places the dividing surface at the point along the reaction coordinate that minimises the TST rate, thereby providing an upper bound that is tighter than the conventional TST result.

For a one-dimensional reaction coordinate $s$ (measured as the distance along the minimum-energy path from reactants), the VTST rate is:

$$k^{\text{VTST}}=\underset{s}{\min }\left[\frac{k_{B}T}{h}\frac{Q^{‡}(s)}{Q_{\text{react}}}{e}^{-V(s)/k_{B}T}\right],$$

where $Q^{‡}(s)$ is the partition function evaluated at the dividing surface orthogonal to $s$ and $V(s)$ is the potential energy along the path. The minimisation over $s$ is performed numerically. In the classical limit, VTST with a good reaction coordinate reduces recrossing errors to a few percent. Combined with tunneling corrections (VTST/MT), the method achieves chemical accuracy ($\sim 1\text{kJ/mol}$) for gas-phase reactions when the electronic structure method provides accurate barrier heights.

Solvent effects and Kramers theory

In solution, the reacting molecules are surrounded by solvent that exerts friction. Henry Eyring's original formulation treats the solvent implicitly through $\Delta G^{‡}$, but the dynamics of barrier crossing in a viscous medium differ fundamentally from the gas phase. Kramers (1940, Physica 7) showed that the rate constant for barrier crossing in a medium with friction coefficient $\gamma$ is:

$$k_{\text{Kramers}}=\frac{{\omega }_0}{2\pi }\left(\sqrt{1+\frac{{\gamma }^2}{4{\omega }_b^2}}-\frac{\gamma }{2{\omega }_b}\right)e^{-\Delta G^{‡}/k_{B}T},$$

where ${\omega }_0$ is the frequency at the bottom of the reactant well and ${\omega }_b$ is the magnitude of the imaginary frequency at the barrier top. In the low-friction limit ($\gamma \to 0$), $k\to ({\omega }_0/2\pi )\exp (-\Delta G^{‡}/k_{B}T)$, recovering the energy-diffusion-controlled regime where the rate is limited by the rate at which the solvent activates the reactant. In the high-friction limit ($\gamma \gg {\omega }_b$), $k\to ({\omega }_b/2\pi )({\omega }_b/\gamma )\exp (-\Delta G^{‡}/k_{B}T)$, the spatial-diffusion-controlled regime where the rate is limited by the solvent viscosity. The TST result corresponds to the turnover region $\gamma \sim {\omega }_b$ where the prefactor is maximised.

The Kramers turnover formula shows that TST provides an upper bound to the true rate in solution, and the transmission coefficient $\kappa <1$ accounts for the friction-induced reduction. Grote and Hynes (1981, J. Chem. Phys. 84) generalised Kramers theory to frequency-dependent friction (memory effects), producing the Grote-Hynes correction that is standard in modern condensed-phase TST.

Limitations of conventional TST

Several assumptions underlie the Eyring equation, each with known failure modes:

Quasi-equilibrium. TST assumes a Boltzmann distribution of reactant energies and that the transition-state population is in equilibrium with reactants. This fails for reactions faster than vibrational relaxation (e.g., fast exothermic radical reactions at low pressure) and for reactions with early barriers where the energy distribution at the transition state deviates from Boltzmann.

No recrossing. The fundamental postulate that trajectories cross the dividing surface only once is violated in viscous solvents, for reactions with late barriers (where the product vibrational energy can feed back into the reaction coordinate), and for systems with multiple saddle points connected by flat regions of the PES.

Classical reaction coordinate. The reaction coordinate is treated as a classical degree of freedom. This fails for reactions involving H-atom transfer at low temperature, where tunneling dominates.

Separability. TST assumes the reaction coordinate is separable from the other $3N-7$ bound modes. This is an approximation: the reaction coordinate couples to the transverse vibrations (curvature of the reaction path), and the dividing surface is not uniquely defined for systems with strongly curved or bifurcating reaction paths.

Born-Oppenheimer PES. The entire TST framework assumes that the reaction occurs on a single electronic potential energy surface. Reactions involving spin-orbit coupling, conical intersections, or non-adiabatic transitions (e.g., electron transfer) require surface-hopping or multi-state TST methods.

Despite these limitations, conventional TST with the Eyring equation remains the workhorse of chemical kinetics because it provides a quantitative molecular interpretation of the Arrhenius parameters at modest computational cost. For most thermal reactions above 200 K, the errors from the TST assumptions are smaller than the uncertainty in the computed barrier height, and the Eyring equation's prediction of $A$ from partition functions is accurate to within an order of magnitude without any empirical input.

Connections Master

• Partition functions for chemical systems 14.07.02. The Eyring equation is a ratio of partition functions: $q^{‡\prime }/(q_A{q}_B/V)$ multiplied by $k_{B}T/h$ and the Boltzmann factor. Every factor in this ratio uses the same factorisation (translational, rotational, vibrational, electronic) developed in 14.07.02. The transition state has one fewer vibrational mode than a stable molecule with the same atoms, because the reaction coordinate is removed.

• Chemical kinetics: rate laws and the Arrhenius equation 14.08.01. The Eyring equation provides the molecular interpretation of the Arrhenius parameters $A$ and $E_a$. The activation energy is $E_a=\Delta H^{‡}+nRT$, and the pre-exponential factor encodes the entropy of activation. The Eyring plot ($\ln (k/T)$ vs. $1/T$) is the statistical-mechanical analogue of the Arrhenius plot ($\ln k$ vs. $1/T$).

• Chemical equilibrium 14.06.04. The equilibrium constant for the activation step $K^{‡}=\exp (-\Delta G^{‡}/RT)$ uses the same thermodynamic machinery as equilibrium constants for overall reactions. The Hammond postulate connects the position of the transition state along the reaction coordinate to the overall thermodynamics.

• Quantum chemistry and computational chemistry 14.04.04. The transition-state structure, vibrational frequencies, and energy are computed using quantum chemistry (DFT, coupled-cluster). Saddle-point optimisation algorithms and intrinsic reaction coordinate (IRC) calculations locate the transition state on the PES. The accuracy of the computed rate constant is limited by the accuracy of the electronic structure method for the barrier height.

• Spectroscopy 14.12.01. The partition functions of reactants are evaluated from spectroscopic constants (rotational constants, vibrational frequencies). Transition-state frequencies come from quantum chemistry rather than spectroscopy, because the activated complex cannot be observed directly.

• Enzyme kinetics. Enzyme-catalysed reactions have lower $\Delta G^{‡}$ than the uncatalysed reaction. The enzyme stabilises the transition state through specific interactions (hydrogen bonds, electrostatic complementarity, strain), lowering $\Delta H^{‡}$ by 40--80 kJ/mol. The entropy of activation is typically more negative for the enzymatic reaction (the substrate loses conformational freedom upon binding), but the enthalpy saving dominates.

• Atmospheric chemistry. The rate constants for key atmospheric reactions ($\text{OH}+{\text{CH}}_4$, $\text{Cl}+{\text{O}}_3$, ${\text{NO}}_2$ photolysis) are computed using TST with tunneling corrections. The NASA/JPL kinetic data evaluation uses TST as the primary theoretical framework for temperature extrapolation of rate constants.

Historical context Master

Henry Eyring developed transition-state theory in 1934--1935 while at Princeton, stimulated by the problem of computing chemical reaction rates from molecular properties. His 1935 paper "The activated complex in chemical reactions" (J. Chem. Phys. 3, 107) is one of the most cited papers in physical chemistry. Eyring recognised that the potential energy surface -- computable, in principle, from the Schrodinger equation -- contains all the information needed to predict a reaction rate, and that the saddle-point geometry provides a natural dividing surface between reactants and products. The key mathematical insight is the cancellation of the unknown reaction-coordinate frequency ${\nu }^{‡}$ against the translational partition function factor $k_{B}T/h{\nu }^{‡}$, leaving the universal prefactor $k_{B}T/h$.

Eyring worked independently of Michael Polanyi and Meredith Evans at Manchester, who published a similar formulation in the same year (Evans and Polanyi 1935, Trans. Faraday Soc. 31, 875). The two groups arrived at equivalent results from different starting points: Eyring from partition functions and statistical mechanics, Evans and Polanyi from an empirical energy-barrier model. The Eyring formulation is more general and is the basis for all modern TST implementations.

Eugene Wigner, who was at Princeton during the same period, contributed the tunneling correction (Wigner 1932, Z. Phys. Chem. B 19) and later provided a rigorous dynamical foundation for TST through phase-space dividing-surface arguments (Wigner 1938, Trans. Faraday Soc. 34, 29). Wigner showed that TST gives an upper bound to the true rate constant because it counts all trajectories that cross the dividing surface in the product direction, including those that recross. The variational principle for TST -- minimising the rate over the choice of dividing surface -- was implicit in Wigner's work and made explicit by Keck (1960, J. Chem. Phys. 32) and later by Truhlar and Garrett (1979, J. Phys. Chem. 83).

The thermodynamic formulation of TST ($\Delta G^{‡}$, $\Delta H^{‡}$, $\Delta S^{‡}$) was developed by Evans and Polanyi (1935) and widely adopted because it connects directly to measurable Arrhenius parameters. The Eyring plot ($\ln (k/T)$ vs. $1/T$) became the standard tool for extracting activation thermodynamics from rate data, replacing the Arrhenius plot when mechanistic interpretation was desired.

The computational implementation of TST was revolutionised by the development of quantum chemistry methods for locating saddle points on potential energy surfaces. The Berny optimisation algorithm (Schlegel 1982) and the synchronous transit methods (Peng and Schlegel 1993) made saddle-point searches routine. Combined with DFT (which provides barrier heights accurate to $\sim 5\text{kJ/mol}$ for many reactions) and higher-level methods (CCSD(T) with large basis sets, accurate to $\sim 1\text{kJ/mol}$), TST became a predictive tool. Modern computational kinetics (Polyrate, GaussRate, KiSThelP) automates the entire workflow: PES computation, saddle-point location, frequency analysis, partition-function evaluation, tunneling correction, and VTST rate-constant calculation.

The Kramers theory of barrier crossing in solution (Kramers 1940) extended TST to condensed phases by introducing friction as a dynamical parameter. The Grote-Hynes generalisation (1981) incorporated frequency-dependent friction and is the standard framework for condensed-phase TST. Transition-state theory in solution differs from gas-phase TST because the solvent provides both a frictional damping and a fluctuating force that can assist or oppose barrier crossing.

Bibliography Master

@article{Eyring1935,
  author = {Eyring, H.},
  title = {The Activated Complex in Chemical Reactions},
  journal = {Journal of Chemical Physics},
  volume = {3},
  year = {1935},
  pages = {107--115}
}

@article{Wigner1938,
  author = {Wigner, E.},
  title = {The Transition State Method},
  journal = {Transactions of the Faraday Society},
  volume = {34},
  year = {1938},
  pages = {29--41}
}

@article{EvansPolanyi1935,
  author = {Evans, M. G. and Polanyi, M.},
  title = {Some Applications of the Transition State Method to the Calculation of Reaction Velocities, Especially in Solution},
  journal = {Transactions of the Faraday Society},
  volume = {31},
  year = {1935},
  pages = {875--894}
}

@article{Kramers1940,
  author = {Kramers, H. A.},
  title = {Brownian Motion in a Field of Force and the Diffusion Model of Chemical Reactions},
  journal = {Physica},
  volume = {7},
  year = {1940},
  pages = {284--304}
}

@article{GroteHynes1981,
  author = {Grote, R. F. and Hynes, J. T.},
  title = {The Stable States Picture of Chemical Reactions. {II.} Rate Constants for Condensed and Gas Phase Reaction Models},
  journal = {Journal of Chemical Physics},
  volume = {74},
  year = {1981},
  pages = {4465--4475}
}

@article{TruhlarGarrett1979,
  author = {Truhlar, D. G. and Garrett, B. C.},
  title = {Variational Transition-State Theory},
  journal = {Accounts of Chemical Research},
  volume = {13},
  year = {1980},
  pages = {440--448}
}

@book{Laidler1987,
  author = {Laidler, K. J.},
  title = {Chemical Kinetics},
  edition = {3},
  publisher = {Harper \& Row},
  year = {1987}
}

@book{SteinfeldFranciscoHase1999,
  author = {Steinfeld, J. I. and Francisco, J. S. and Hase, W. L.},
  title = {Chemical Kinetics and Dynamics},
  edition = {2},
  publisher = {Prentice Hall},
  year = {1999}
}

@book{AtkinsPaula2023,
  author = {Atkins, P. and de Paula, J.},
  title = {Physical Chemistry},
  edition = {12},
  publisher = {Oxford University Press},
  year = {2023}
}

@article{Hammond1955,
  author = {Hammond, G. S.},
  title = {A Correlation of Reaction Rates},
  journal = {Journal of the American Chemical Society},
  volume = {77},
  year = {1955},
  pages = {334--338}
}

@article{Bell1933,
  author = {Bell, R. P.},
  title = {The Application of Quantum Mechanics to Chemical Kinetics},
  journal = {Proceedings of the Royal Society of London A},
  volume = {139},
  year = {1933},
  pages = {466--474}
}

@book{Johnston1966,
  author = {Johnston, H. S.},
  title = {Gas Phase Reaction Rate Theory},
  publisher = {Ronald Press},
  year = {1966}
}

@book{Forst2003,
  author = {Forst, W.},
  title = {Unimolecular Reactions: A Concise Introduction},
  publisher = {Cambridge University Press},
  year = {2003}
}