Physical

Transition state theory: the Eyring equation

12 min read Updated

Arrhenius fits the temperature dependence of a rate constant and reports two numbers: A (the pre-exponential factor) and E_a (the activation energy). Both are empirical. Transition state theory (TST) — Henry Eyring’s 1935 upgrade — reframes those two numbers into molecular properties of a specific chemical species: the activated complex sitting at the top of the reaction barrier. Eyring’s equation gives ΔH‡ and ΔS‡, activation enthalpy and activation entropy, and separates “why is this reaction slow?” into two independent questions with different remedies.

Assume familiarity with Arrhenius and integrated rate laws from rate laws from mechanism.

The Eyring equation

Eyring’s derivation (McQuarrie §29) treats the activated complex — the arrangement of atoms at the top of the reaction-coordinate barrier — as being in quasi-equilibrium with the reactants. Applying statistical mechanics to that equilibrium gives the rate constant:

k=κ⋅kBThexp⁡(ΔS‡R)exp⁡(−ΔH‡RT)k = \kappa \cdot \frac{k_B T}{h} \exp\left(\frac{\Delta S^{\ddagger}}{R}\right) \exp\left(-\frac{\Delta H^{\ddagger}}{RT}\right)

with:

  • k_B = Boltzmann constant = 1.380 649 × 10⁻²³ J K⁻¹ (CODATA 2018, exact).
  • h = Planck constant = 6.626 070 15 × 10⁻³⁴ J s (CODATA 2018, exact).
  • T = temperature in K.
  • R = gas constant = 8.314 462 618 J K⁻¹ mol⁻¹.
  • κ = transmission coefficient (fraction of activated complexes that actually make it over the barrier to products; usually taken as 1 for simple reactions).
  • ΔH‡ = activation enthalpy (J/mol) — the “height” of the barrier in energy.
  • ΔS‡ = activation entropy (J K⁻¹ mol⁻¹) — the “narrowness” of the barrier in configurational space.

The prefactor k_B T/h ≈ 6.21 × 10¹² s⁻¹ at 298 K — sometimes called the “universal frequency factor.” Every rate constant, in Eyring’s framework, is that number times two exponentials.

Compare Arrhenius:

k=Aexp⁡(−Ea/RT)k = A \exp(-E_a / RT)

Correspondence: E_a ≈ ΔH‡ + RT (usually within 2-5 kJ/mol at room T), and A ≈ (k_B T / h) exp(ΔS‡/R). Arrhenius A absorbs the entropy factor; TST separates them.

What ΔH‡ tells you

ΔH‡ is the enthalpic cost of reaching the transition state — bond stretches, angle bends, partial bond breakings that the reactants must undergo. Larger ΔH‡ means the reaction is intrinsically hard, requires more energy input per event, and slows dramatically as T decreases.

Diagnostic values:

  • Diffusion-controlled reactions in water: ΔH‡ ≈ 15-20 kJ/mol (activation energy is essentially just the energy to break the solvent cage).
  • Fast enzymatic reactions: ΔH‡ = 25-40 kJ/mol.
  • Typical bench organic reactions at 25 °C: ΔH‡ = 50-90 kJ/mol.
  • Slow reactions (days at 25 °C): ΔH‡ = 90-120 kJ/mol.
  • Reactions requiring reflux (80 °C): ΔH‡ = 90-110 kJ/mol.

Reducing ΔH‡ is the traditional goal of catalysis. Enzymes achieve rate enhancements of 10⁶-10¹⁰ almost entirely by lowering ΔH‡; heterogeneous catalysts (Pd/C, Pt/C) do the same for hydrogenation and coupling reactions.

What ΔS‡ tells you

ΔS‡ is the entropic cost of assembling the transition state from reactants. If forming the TS requires two independent molecules to line up in a specific relative orientation and freeze several rotational degrees of freedom, ΔS‡ is strongly negative (order-of-magnitude −100 J K⁻¹ mol⁻¹). If the TS is looser than the reactant — for example, a unimolecular dissociation that partially breaks a bond — ΔS‡ is positive.

Diagnostic values:

  • Bimolecular associations (SN2, Diels-Alder, radical recombination): ΔS‡ ≈ −80 to −150 J K⁻¹ mol⁻¹. Very negative.
  • Unimolecular dissociations (SN1 ionization, retro-Diels-Alder): ΔS‡ ≈ +10 to +30 J K⁻¹ mol⁻¹. Slightly positive.
  • Isomerizations, conformational changes: ΔS‡ ≈ 0 to −20 J K⁻¹ mol⁻¹. Near zero.
  • Enzyme reactions with well-positioned substrates: ΔS‡ often less negative than the uncatalyzed case, because the enzyme has pre-organized the reactants — one of enzymes’ entropic tricks.

The sign of ΔS‡ is a mechanism diagnostic. A strongly negative ΔS‡ almost always means an associative mechanism; a positive ΔS‡ often means a dissociative mechanism. Steinfeld Ch. 10 has case studies of TST used this way in organic mechanistic assignment.

Extracting ΔH‡ and ΔS‡ from rate data

Take the log of the Eyring equation:

ln⁡kT=−ΔH‡RT+ln⁡kBh+ΔS‡R\ln\frac{k}{T} = -\frac{\Delta H^{\ddagger}}{R T} + \ln\frac{k_B}{h} + \frac{\Delta S^{\ddagger}}{R}

A plot of ln(k/T) vs 1/T is linear. Slope = −ΔH‡/R, intercept = ln(k_B/h) + ΔS‡/R. This is the Eyring plot — the same idea as an Arrhenius plot but with the y-axis rescaled to divide by T, and the intercept now returns ΔS‡ rather than ln A.

From two data points at temperatures T₁ and T₂ with rate constants k₁ and k₂:

ΔH‡=R⋅ln⁡(k1T2/k2T1)1/T2−1/T1\Delta H^{\ddagger} = R \cdot \frac{\ln(k_1 T_2 / k_2 T_1)}{1/T_2 - 1/T_1}

Then ΔS‡ falls out of Eyring’s equation at either temperature. Modern kinetics software does this automatically; the half-life calculator does not do Eyring plots, but any Arrhenius-based analysis converts to Eyring with the same data.

Worked example — SN2 reaction from two temperatures

Worked example — SN2 of CH3I with OH-: computing DHt and DSt
T1
298 K, k1 = 2.5e-4 M^-1 s^-1
T2
318 K, k2 = 1.7e-3 M^-1 s^-1
Constants
kB/h = 2.084e10 K^-1 s^-1 at 298 K; R = 8.314 J K^-1 mol^-1
  1. Formula
  2. Substitute
  3. Note

    ln(k1 T2 / (k2 T1)) = ln(2.5e-4 x 318 / (1.7e-3 x 298)) = ln(0.0795 / 0.5066) = ln(0.1569) = -1.852. 1/T2 - 1/T1 = 1/318 - 1/298 = 3.145e-3 - 3.356e-3 = -2.115e-4 K^-1. DHt = 8.314 x (-1.852) / (-2.115e-4) = 72,800 J/mol ~ 73 kJ/mol (close to the 70 kJ/mol textbook Laidler value for methyl iodide SN2 with hydroxide). For DSt, use Eyring at T1: 2.5e-4 = (1.381e-23 x 298 / 6.626e-34) exp(DSt/R) exp(-73000 / (8.314 x 298)). kBT/h = 6.21e12 s^-1. exp(-73000/2477) = exp(-29.47) = 1.6e-13. So 2.5e-4 = 6.21e12 x 1.6e-13 x exp(DSt/R) = 0.994 exp(DSt/R). exp(DSt/R) = 2.5e-4. DSt/R = -8.29. DSt = -8.29 x 8.314 = -69 J K^-1 mol^-1 (close to the -85 J K^-1 mol^-1 order-of-magnitude expectation for a bimolecular SN2). The mechanistic reading: a substantial enthalpic barrier (73 kJ/mol - typical of C-I bond breaking) and a strongly negative entropy of activation (bimolecular association bringing OH- and CH3I into a specific back-side approach - classic SN2 signature).

  4. Result

Why TST helps where Arrhenius does not

Two mechanisms that give the same Arrhenius fit can be distinguished by TST. Consider two hypothetical reactions:

  • Reaction A: ΔH‡ = 60 kJ/mol, ΔS‡ = −20 J K⁻¹ mol⁻¹ (mild association).
  • Reaction B: ΔH‡ = 40 kJ/mol, ΔS‡ = −150 J K⁻¹ mol⁻¹ (strong association).

At 298 K both give similar k (≈ 10⁻³ s⁻¹), so an Arrhenius plot at 25 °C alone cannot tell them apart. But the temperature dependence of the temperature dependence differs: Reaction A speeds up more with T (steeper ln k vs 1/T slope). If you have data over a wide T range, the Eyring plot cleanly separates them. Steinfeld Ch. 10 has real examples where mechanistic assignment turned on exactly this discrimination.

TST also handles isotope effects (KIE) more cleanly: the ratio k_H/k_D can be computed from the difference in zero-point vibrational energies of the C-H versus C-D bond, and the ratio’s magnitude tells you whether the H atom is being transferred at the TS. Primary KIE k_H/k_D = 6-8 → H being broken in the RDS. Secondary KIE k_H/k_D = 0.9-1.1 → H not directly involved but the bond hybridization changes.

Practice problems

Problem 1 Easy A reaction has DHt = 50 kJ/mol and DSt = -100 J K^-1 mol^-1 at 298 K. What is the rate constant?

Answer: k = (kB T/h) exp(DSt/R) exp(-DHt/RT). kBT/h at 298 K = 6.21e12 s^-1. exp(-100/8.314) = exp(-12.03) = 5.96e-6. exp(-50000/(8.314 x 298)) = exp(-20.18) = 1.72e-9. k = 6.21e12 x 5.96e-6 x 1.72e-9 = 6.4e-2 s^-1. Half-life ~ 11 s (first order).

Problem 2 Intermediate An SN1 reaction gives k = 3.2e-6 s^-1 at 25 C. If DHt = 90 kJ/mol, estimate DSt.

Answer: Rearrange Eyring: DSt/R = ln(k h / (kB T)) + DHt/RT. ln(3.2e-6 / 6.21e12) = ln(5.15e-19) = -42.11. DHt/RT = 90000 / (8.314 x 298) = 36.32. DSt/R = -42.11 + 36.32 = -5.79. DSt = -48 J K^-1 mol^-1. Somewhat negative — smaller in magnitude than the SN2 example (-85), consistent with SN1's unimolecular RDS being less entropically demanding than a bimolecular SN2 attack. Note: real SN1 reactions often show slightly positive DSt because the ionization TS has stretched, more disordered bonds; the -48 here is on the transition between the two regimes.

Problem 3 Hard At 298 K a Diels-Alder reaction has k = 1e-3 M^-1 s^-1. At 348 K, k = 3e-2 M^-1 s^-1. Compute DHt and estimate DSt.

Answer: ln(k1 T2 / (k2 T1)) = ln(1e-3 x 348 / (3e-2 x 298)) = ln(0.0389) = -3.246. 1/T2 - 1/T1 = 1/348 - 1/298 = -4.81e-4. DHt = R x (-3.246) / (-4.81e-4) = 56.1 kJ/mol. For DSt at 298 K: ln(1e-3 x 6.626e-34 / (1.381e-23 x 298)) = ln(1.61e-16) = -36.36. DSt/R = -36.36 + 56100/(8.314 x 298) = -36.36 + 22.65 = -13.7. DSt = -114 J K^-1 mol^-1. Consistent with a bimolecular cyclic-TS reaction — diene and dienophile line up in a specific relative orientation with two forming bonds and a locked ring, high entropic cost.

Where to go next

  • The next layer of theory is the Marcus equation for electron-transfer reactions, which further decomposes ΔH‡ into reorganization energy and driving-force terms. Steinfeld Ch. 10 covers Marcus theory.
  • For heterogeneous catalysis and enzyme active sites, TST extends via the Bronsted-Evans-Polanyi (BEP) relation: linear correlation between ΔH‡ and reaction ΔH across a series of similar reactions. It is the theoretical foundation of most scaling relations in catalysis research.
  • The half-life calculator uses first-order kinetics; once you have k from Eyring, half-life follows immediately.
  • The sig figs cheatsheet reminds that activation-parameter uncertainty compounds quickly through the exponential; report ΔH‡ to ±2 kJ/mol and ΔS‡ to ±5 J K⁻¹ mol⁻¹ unless your data are exceptional.

TST’s real value is diagnostic. Arrhenius gives you two numbers that fit; Eyring gives you two numbers that mean something about the transition state. When two mechanisms fit the same Arrhenius plot, ΔS‡ often names which one is running.

Sources

  • Atkins & de Paula, Physical Chemistry, 11th ed. (2018), Ch. 21 — Reaction dynamics
  • McQuarrie & Simon, Physical Chemistry: A Molecular Approach (1997), Ch. 29-30
  • Steinfeld, Francisco & Hase, Chemical Kinetics and Dynamics, 2nd ed. (1998), Ch. 10
  • H. Eyring (1935), J. Chem. Phys. 3, 107 — Original transition-state theory paper
  • K. J. Laidler & M. C. King (1983), J. Phys. Chem. 87, 2657 — TST retrospective
  • CODATA 2018 — k_B = 1.380 649 x 10^-23 J K^-1, h = 6.626 070 15 x 10^-34 J s