Transition State Theory - Physical Chemistry 2

SECTION 1

Historical Context & Motivation

The quest to understand why chemical reactions proceed at particular rates stretches back to the earliest kinetic studies of the nineteenth century. The Arrhenius equation, proposed in 1889, provided an empirical relationship between the rate constant k and temperature through an activation energy parameter, but it offered no microscopic explanation of what the activation energy actually represented. Meanwhile, collision theory — developed primarily by Max Trautz and William Lewis — modeled reactions as hard-sphere collisions with an energy threshold, yet consistently overestimated rate constants for complex molecules because it ignored the internal degrees of freedom and the geometry of molecular encounters.

The deficiency was clear: chemists needed a theory that could bridge molecular structure with reaction kinetics. Specifically, a framework was required that could account for the rearrangement of bonds during a reaction, the role of vibrational and rotational modes, and the statistical likelihood that a system could reach the critical configuration necessary for transformation. This intellectual gap drove the development of Transition State Theory in the 1930s.

1889 • Arrhenius Equation

Svante Arrhenius publishes his empirical rate-temperature relationship k = A·exp(−E_a/RT), introducing the concept of an activation energy but leaving its physical origin unexplained.

1918 • Collision Theory Formalized

Trautz and Lewis develop collision theory, treating molecules as hard spheres. The theory succeeds for atom–atom reactions but fails for polyatomic systems due to its neglect of molecular geometry and internal energy distributions.

1931 • Potential Energy Surfaces

Henry Eyring and Michael Polanyi construct the first quantum-mechanical potential energy surface (PES) for the H + H₂ reaction, providing the conceptual landscape on which transition state theory would be built.

1935 • Birth of TST

Henry Eyring, and independently Meredith Gwynne Evans and Michael Polanyi, publish the formal Transition State Theory (also called Activated Complex Theory or Absolute Rate Theory), providing a statistical-mechanical expression for rate constants.

1960s–present • Extensions & Refinements

Variational TST, quantum tunneling corrections (Wigner, Eckart), and computational chemistry (DFT-based PES scans) extend TST's accuracy and applicability to complex systems including enzyme catalysis and atmospheric reactions.

The central question that Transition State Theory answered was deceptively simple: Can we predict a reaction's rate constant from first principles — using only the molecular properties of reactants and the fleeting activated complex at the reaction's energetic summit? Eyring's affirmative answer transformed chemical kinetics from an empirical discipline into one rooted in statistical mechanics and quantum chemistry.

SECTION 2

Core Principles & Foundational Assumptions

Transition State Theory rests on a set of well-defined assumptions that enable a tractable statistical-mechanical treatment of reaction rates. Understanding these assumptions is essential not only for applying the theory but also for recognizing its limitations. The theory envisions a reaction as a journey across a potential energy surface, where reactants start in one valley (energy minimum), pass through a saddle point — the transition state — and descend into the product valley. The saddle-point configuration is called the activated complex, denoted AB^‡.

Quasi-Equilibrium Hypothesis

The activated complex AB^‡ is in thermodynamic quasi-equilibrium with the reactants. This allows us to write an equilibrium constant K^‡ between reactants and the transition state, applying all standard statistical-mechanical partition function machinery.

Separability of the Reaction Coordinate

One vibrational mode of the activated complex — the reaction coordinate — is separable from all other degrees of freedom and is treated as a translational motion along the PES ridge. This mode has an imaginary frequency, corresponding to the bond-breaking/forming motion.

No Recrossing

Once the system crosses the transition state dividing surface in the forward direction, it proceeds to products without returning. This assumption leads to an upper-bound estimate of the true rate constant, since real trajectories may recross the barrier.

Classical Nuclear Motion

Nuclear motion over the barrier is treated classically — quantum tunneling through the barrier is neglected in the basic theory. Corrections for tunneling (e.g., the Wigner tunneling correction) can be added but are not part of the canonical formulation.

Born–Oppenheimer Separation

Electronic and nuclear motions are separable, so the system moves on a single, well-defined potential energy surface throughout the reaction. Non-adiabatic transitions between electronic states are not considered.

✦

KEY TAKEAWAY

Think of TST like modeling freeway traffic through a mountain pass. The pass (transition state) is the narrowest, highest point. TST assumes that the density of cars approaching the pass is in a steady state (quasi-equilibrium), that every car that crests the summit continues downhill without making a U-turn (no recrossing), and that no car can tunnel through the mountain (classical motion). These simplifications let you calculate traffic flow from the geometry and altitude of the pass alone — analogous to computing a rate constant from the properties of the activated complex.

SECTION 3

The Potential Energy Surface & Transition State

The conceptual heart of TST is the potential energy profile along the reaction coordinate. This one-dimensional slice through the multidimensional PES plots energy as a function of the progress of the reaction, from reactants through the transition state to products. The diagram below illustrates the essential features: the reactant well, the activation energy barrier ΔG^‡, the activated complex at the saddle point, and the product well. The difference in free energy between the two wells is ΔG^° for the overall reaction, while the height of the barrier above the reactant level determines how fast the reaction proceeds.

The transition state (AB‡) sits at the saddle-point maximum along the reaction coordinate. The height ΔG‡ is the Gibbs free energy of activation, which governs the rate. The overall thermodynamic driving force is ΔG°, the standard free energy change of the reaction.

Several critical features of this diagram warrant emphasis. First, the transition state is not a stable species — it exists for roughly one vibrational period (≈ 10⁻¹³ s) and cannot be isolated or directly observed under normal circumstances. Second, the reaction coordinate is a composite variable that may represent the simultaneous stretching and compression of multiple bonds. Third, the barrier height ΔG^‡ is a free-energy quantity — it includes both the enthalpic cost of distorting bonds and the entropic cost of organizing the reactants into the transition-state geometry.

SECTION 4

Mathematical Framework

The derivation of the TST rate constant begins with the quasi-equilibrium assumption. Consider a bimolecular elementary reaction A + B → Products. TST posits a quasi-equilibrium between reactants and the activated complex: A + B ⇌ AB^‡ → Products. The equilibrium concentration of the activated complex is related to the reactant concentrations through the equilibrium constant K^‡, which can be expressed in terms of molecular partition functions.

The Eyring Equation

EYRING EQUATION (THERMODYNAMIC FORM)

k = (k_B T / h) × exp(−ΔG‡ / RT)

Where k = rate constant, k_B = Boltzmann constant (1.381 × 10⁻²³ J·K⁻¹), T = absolute temperature (K), h = Planck's constant (6.626 × 10⁻³⁴ J·s), ΔG‡ = Gibbs free energy of activation (J·mol⁻¹), R = gas constant (8.314 J·mol⁻¹·K⁻¹). The prefactor k_BT/h has units of s⁻¹ and represents the universal frequency at which activated complexes decompose to products.

The elegance of the Eyring equation lies in the prefactor k_BT/h. It arises from treating the reaction coordinate as a one-dimensional translational mode. At 298 K, this factor equals approximately 6.21 × 10¹² s⁻¹, meaning the activated complex vibrates at roughly 10¹³ Hz along the decomposition coordinate. The exponential term provides the fraction of systems that have sufficient energy to reach the transition state.

Enthalpy–Entropy Decomposition

ENTHALPY–ENTROPY FORM

k = (k_B T / h) × exp(ΔS‡ / R) × exp(−ΔH‡ / RT)

Expanding ΔG‡ = ΔH‡ − TΔS‡ separates the rate constant into enthalpic and entropic contributions. ΔH‡ (enthalpy of activation) reflects the energy required to distort bonds to the transition-state geometry. ΔS‡ (entropy of activation) captures the change in order — a negative ΔS‡ indicates a more constrained transition state relative to the reactants.

Statistical-Mechanical Form

PARTITION FUNCTION FORM

k = (k_B T / h) × (q‡ / q_A · q_B) × exp(−ΔE₀ / RT)

Here q‡ is the partition function of the activated complex (with the reaction-coordinate mode removed), q_A and q_B are the partition functions of reactants A and B, and ΔE₀ is the zero-point-corrected barrier height. All partition functions are evaluated per unit volume.

ℹ

Connecting to Arrhenius

Comparing the Eyring equation with the Arrhenius expression k = A·exp(−E_a/RT), one can show that E_a = ΔH‡ + RT for a unimolecular reaction and E_a = ΔH‡ + 2RT for a bimolecular reaction in the gas phase. Thus, TST provides a microscopic interpretation of the empirical Arrhenius parameters.

SECTION 5

The Eyring Plot & Extracting Activation Parameters

Just as the Arrhenius plot (ln k vs. 1/T) yields the activation energy from its slope, the Eyring plot provides both the enthalpy and entropy of activation. Taking the natural logarithm of the Eyring equation and rearranging gives a linear relationship: ln(k/T) = −ΔH‡/R × (1/T) + ln(k_B/h) + ΔS‡/R. By plotting ln(k/T) against 1/T, the slope yields −ΔH‡/R and the y-intercept yields ln(k_B/h) + ΔS‡/R.

An Eyring plot linearizes the Eyring equation. The slope gives −ΔH‡/R, and the intercept gives ln(k_B/h) + ΔS‡/R. Measuring rate constants at several temperatures and performing linear regression provides both activation parameters from a single experiment.

LINEARIZED EYRING EQUATION

ln(k/T) = −(ΔH‡/R) × (1/T) + ln(k_B/h) + ΔS‡/R

At 298 K, ln(k_B/h) = ln(1.381 × 10⁻²³ / 6.626 × 10⁻³⁴) = 23.76. This constant is useful when computing ΔS‡ from the intercept.

The Eyring plot is particularly powerful because it separates the enthalpic and entropic contributions to the activation barrier. For instance, a reaction with a small ΔH‡ but a very negative ΔS‡ is entropically disfavored — the transition state is much more ordered than the reactants, as is typical for bimolecular association reactions where translational degrees of freedom are lost. Conversely, unimolecular dissociations often have positive ΔS‡ because the transition state is looser than the reactant molecule.

SECTION 6

Worked Example

Let us apply the Eyring equation to compute the activation parameters for a hypothetical reaction using temperature-dependent rate-constant data. This exercise demonstrates the complete workflow from experimental data to ΔH‡ and ΔS‡.

Determining ΔH‡ and ΔS‡ from an Eyring Plot

1
Step 1 — State the Problem

The rate constant for a unimolecular isomerization reaction was measured at five temperatures. At T₁ = 300 K, k₁ = 2.50 × 10⁻⁴ s⁻¹; at T₂ = 320 K, k₂ = 1.45 × 10⁻³ s⁻¹; at T₃ = 340 K, k₃ = 7.20 × 10⁻³ s⁻¹; at T₄ = 360 K, k₄ = 3.10 × 10⁻² s⁻¹; at T₅ = 380 K, k₅ = 1.15 × 10⁻¹ s⁻¹. Determine ΔH‡ and ΔS‡ using an Eyring plot.
2
Step 2 — Compute ln(k/T) and 1/T

For each data point, calculate 1/T (in K⁻¹) and ln(k/T). For T₁ = 300 K: 1/T = 3.333 × 10⁻³ K⁻¹, k/T = 2.50 × 10⁻⁴ / 300 = 8.33 × 10⁻⁷, so ln(k/T) = −14.00. For T₂ = 320 K: 1/T = 3.125 × 10⁻³ K⁻¹, ln(k/T) = ln(4.53 × 10⁻⁶) = −12.30. Repeating for T₃: ln(k/T) = −10.76; T₄: ln(k/T) = −9.37; T₅: ln(k/T) = −8.10.

Five data pairs (1/T, ln(k/T)) ready for linear regression.
3
Step 3 — Perform Linear Regression

Plotting ln(k/T) vs. 1/T and performing a least-squares fit yields slope m = −8,530 K and intercept b = 14.43. The slope has units of kelvin because 1/T is in K⁻¹ and ln(k/T) is dimensionless.

Slope = −8,530 K, Intercept = 14.43
4
Step 4 — Extract ΔH‡

From the Eyring equation, slope = −ΔH‡/R. Therefore: ΔH‡ = −slope × R = 8,530 K × 8.314 J·mol⁻¹·K⁻¹ = 70,920 J·mol⁻¹ ≈ 70.9 kJ·mol⁻¹.

ΔH‡ = 70.9 kJ·mol⁻¹
5
Step 5 — Extract ΔS‡

The intercept equals ln(k_B/h) + ΔS‡/R = 23.76 + ΔS‡/R. Setting this equal to the observed intercept: 14.43 = 23.76 + ΔS‡/R, so ΔS‡/R = −9.33, giving ΔS‡ = −9.33 × 8.314 = −77.6 J·mol⁻¹·K⁻¹. The negative value indicates that the transition state is more ordered than the reactant, consistent with a cyclic transition state in an isomerization.

ΔS‡ = −77.6 J·mol⁻¹·K⁻¹
6
Step 6 — Verify with ΔG‡ at 300 K

ΔG‡ = ΔH‡ − TΔS‡ = 70,900 − 300 × (−77.6) = 70,900 + 23,280 = 94,180 J·mol⁻¹ = 94.2 kJ·mol⁻¹. Plugging back into the Eyring equation: k = (1.381 × 10⁻²³ × 300) / (6.626 × 10⁻³⁴) × exp(−94,180 / (8.314 × 300)) = 6.25 × 10¹² × exp(−37.81) = 6.25 × 10¹² × 3.93 × 10⁻¹⁷ = 2.46 × 10⁻⁴ s⁻¹, which agrees with the experimental value of 2.50 × 10⁻⁴ s⁻¹ to within rounding error.

Calculated k(300 K) = 2.46 × 10⁻⁴ s⁻¹ ✓ — consistent with experimental data.

SECTION 7

Strengths, Limitations, and Comparison with Collision Theory

No theory is universally applicable, and understanding the boundaries of TST is as important as mastering its equations. The table below contrasts the strengths and limitations of Transition State Theory, and compares it to the simpler collision theory to highlight where TST provides genuinely superior predictions.

Feature	Collision Theory	Transition State Theory
Molecular model	Hard-sphere collisions; no internal structure	Full molecular structure via partition functions; vibrational, rotational, and translational modes included
Steric factor	Introduced empirically as an adjustable parameter P (0 < P ≤ 1)	Emerges naturally from the geometry and symmetry of the activated complex
Activation parameters	Predicts E_a only	Provides ΔH‡, ΔS‡, and ΔG‡ separately
Solvent effects	Not easily incorporated	Naturally handled through ΔS‡ (solvation reorganization changes entropy)
Tunneling	Neglected	Neglected in basic theory; corrections available (Wigner, Eckart)
Recrossing	Not applicable	Assumed absent; k_TST is therefore an upper bound
Computational cost	Minimal — requires only molecular masses and collision cross-sections	Moderate to high — requires knowledge of the transition-state structure and its vibrational frequencies

Comparison of Collision Theory and Transition State Theory

✦

KEY TAKEAWAY

TST's greatest strength is also the source of its primary limitation. By assuming quasi-equilibrium and no recrossing, the theory provides a clean, closed-form expression for the rate constant — but it inherently overestimates the true rate. Think of it as predicting traffic flow by assuming every car heading toward a bridge will cross it: the prediction sets an upper bound, and the actual flow may be lower if some drivers turn around before crossing. Variational TST addresses this by optimizing the location of the dividing surface to minimize the calculated flux, providing tighter bounds.

SECTION 8

Connections to Advanced Theory

Classical TST serves as the foundation for several more sophisticated theoretical frameworks. As computational power has grown, so has the ability to refine the approximations inherent in the basic theory. Three major extensions deserve mention: Variational Transition State Theory (VTST), quantum tunneling corrections, and the Marcus theory of electron transfer, which applies TST-like reasoning to redox reactions.

Extension	What It Addresses	Key Improvement
Variational TST (VTST)	Recrossing — conventional TST places the dividing surface at the saddle point, which may not minimize the reactive flux	Optimizes the dividing-surface position along the reaction coordinate to locate the free-energy maximum, providing a tighter upper bound on k
Tunneling Corrections	Quantum-mechanical barrier penetration, especially important for light atoms (H, D) and low-temperature reactions	Multiplicative correction κ (transmission coefficient): k_corrected = κ × k_TST. Wigner: κ ≈ 1 + (1/24)(hν‡/k_BT)²
Marcus Theory	Electron-transfer reactions where the nuclear rearrangement is dominated by solvent reorganization rather than bond breaking	Introduces reorganization energy λ and predicts the inverted region: beyond −ΔG° = λ, rates decrease with increasing driving force
Kramers Theory	Solution-phase reactions where solvent friction affects barrier crossing dynamics	Couples TST with Langevin dynamics; predicts rate turnover in high-friction solvents (diffusion-controlled limit)

Major extensions of classical Transition State Theory

These extensions illustrate that TST is not a dead-end theory but a living framework continually refined by advances in quantum chemistry and molecular dynamics. Modern computational chemistry routinely uses density functional theory (DFT) to locate transition states on multidimensional potential energy surfaces, compute vibrational frequencies of the activated complex, and evaluate rate constants via the partition-function form of the Eyring equation. Enzyme kinetics, atmospheric chemistry, combustion modeling, and materials degradation pathways all rely on TST-based computational predictions.

💡

Looking Forward

In advanced kinetics courses, you will encounter reactive flux methods and ring-polymer molecular dynamics (RPMD), which go beyond TST entirely by computing exact classical or quantum-mechanical rate constants from molecular dynamics simulations. These methods use TST as a starting point — the initial dividing surface — and then correct for recrossing and tunneling dynamically.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain why the TST rate constant is considered an upper bound to the true rate constant. What physical phenomenon does the theory neglect that causes this overestimation, and how does Variational TST attempt to correct it?

Show answer

TST assumes that every trajectory crossing the transition-state dividing surface in the forward direction proceeds to products. In reality, some trajectories recross the dividing surface — they temporarily enter the product region but are reflected back to reactants. Because TST counts these recrossing events as successful reactive events, it overestimates the rate constant. Variational TST addresses this by optimizing the location of the dividing surface along the reaction coordinate. Rather than fixing it at the saddle point, VTST searches for the position that minimizes the calculated flux (equivalently, maximizes the free energy of activation), thereby reducing the overcounting of recrossing trajectories and providing a tighter upper bound.

Question 2: BASIC CALCULATION

Calculate the rate constant at 298 K for a reaction with ΔG‡ = 85.0 kJ·mol⁻¹ using the Eyring equation. Use k_B = 1.381 × 10⁻²³ J·K⁻¹, h = 6.626 × 10⁻³⁴ J·s, and R = 8.314 J·mol⁻¹·K⁻¹.

Show answer

Apply the Eyring equation: k = (k_BT/h) × exp(−ΔG‡/RT). First compute the prefactor: k_BT/h = (1.381 × 10⁻²³ × 298) / (6.626 × 10⁻³⁴) = 6.21 × 10¹² s⁻¹. Next compute the exponent: −ΔG‡/RT = −85,000 / (8.314 × 298) = −34.31, so exp(−34.31) = 1.22 × 10⁻¹⁵. Therefore k = 6.21 × 10¹² × 1.22 × 10⁻¹⁵ = 7.58 × 10⁻³ s⁻¹. This corresponds to a reaction with a half-life of about 91 seconds — moderately fast at room temperature.

Question 3: INTERMEDIATE

An Eyring plot for a bimolecular reaction gives a slope of −12,400 K and a y-intercept of 11.82. Determine ΔH‡ and ΔS‡. Then calculate the Arrhenius activation energy E_a at 300 K, recalling that for a bimolecular gas-phase reaction E_a = ΔH‡ + 2RT.

Show answer

From the slope: ΔH‡ = −slope × R = 12,400 × 8.314 = 103,100 J·mol⁻¹ = 103.1 kJ·mol⁻¹. From the intercept: ΔS‡/R = intercept − ln(k_B/h) = 11.82 − 23.76 = −11.94, so ΔS‡ = −11.94 × 8.314 = −99.3 J·mol⁻¹·K⁻¹. For the Arrhenius activation energy at 300 K: E_a = ΔH‡ + 2RT = 103,100 + 2(8.314)(300) = 103,100 + 4,988 = 108.1 kJ·mol⁻¹. The large negative ΔS‡ indicates significant loss of translational freedom as two molecules combine into a single activated complex.

Question 4: APPLIED

An enzyme-catalyzed reaction has ΔG‡ = 55 kJ·mol⁻¹, while the corresponding uncatalyzed reaction has ΔG‡ = 100 kJ·mol⁻¹ at 310 K (body temperature). Calculate the ratio of the catalyzed to uncatalyzed rate constants (k_cat/k_uncat) and comment on the biological significance.

Show answer

Since both reactions have the same prefactor k_BT/h at the same temperature, the ratio simplifies to k_cat/k_uncat = exp[(ΔG‡_uncat − ΔG‡_cat)/RT] = exp[(100,000 − 55,000)/(8.314 × 310)] = exp[45,000/2,577] = exp(17.46) = 3.83 × 10⁷. The enzyme accelerates the reaction by a factor of nearly 40 million. This enormous rate enhancement arises from a reduction of only 45 kJ·mol⁻¹ in the activation barrier, demonstrating the exponential sensitivity of reaction rates to barrier heights — a mere 17% reduction in ΔG‡ produces a seven-order-of-magnitude speedup, which is the fundamental reason enzymes are essential for life.

Question 5: CRITICAL THINKING

A colleague measures the kinetic isotope effect (KIE) for a C−H bond-breaking reaction and finds k_H/k_D = 12 at 298 K. The maximum primary KIE predicted by TST (considering only zero-point energy differences) is approximately 6.9 at this temperature. Propose an explanation for why the observed KIE exceeds the semiclassical TST prediction, and discuss what this implies about the assumptions of conventional TST.

Show answer

A k_H/k_D value of 12, significantly exceeding the semiclassical TST maximum of ~6.9, strongly suggests that quantum mechanical tunneling contributes significantly to the reaction rate. Hydrogen, being lighter than deuterium, tunnels through the barrier more effectively, inflating k_H relative to k_D beyond what zero-point energy differences alone can explain. This observation directly violates the classical nuclear motion assumption of conventional TST (Principle 4). To reconcile theory with experiment, one must apply tunneling corrections (e.g., the Wigner or Eckart correction, or multidimensional tunneling approaches like the small-curvature tunneling approximation). The Wigner correction κ ≈ 1 + (1/24)(hν‡/k_BT)² is isotope-dependent because ν‡ scales as m⁻¹/², increasing the effective KIE. Large KIE values are commonly observed in enzymatic reactions involving hydrogen atom or proton transfers, where the narrow, high barrier provides favorable tunneling conditions.

SUMMARY

Summary

Transition State Theory provides a powerful statistical-mechanical framework for computing reaction rate constants from the molecular properties of reactants and the activated complex. The theory rests on the quasi-equilibrium hypothesis (the transition state is in equilibrium with reactants), the no-recrossing assumption (every trajectory crossing the dividing surface proceeds to products), and classical nuclear motion (no quantum tunneling). The central result is the Eyring equation: k = (k_BT/h) × exp(−ΔG‡/RT), which decomposes the activation barrier into enthalpic (ΔH‡) and entropic (ΔS‡) contributions, both extractable from the Eyring plot of ln(k/T) vs. 1/T.

While TST yields an upper-bound rate constant due to the no-recrossing assumption, extensions such as Variational TST and tunneling corrections refine the theory's accuracy. TST connects microscopic molecular properties — partition functions, vibrational frequencies, and potential energy surfaces — to macroscopic observables like rate constants and activation energies, making it an indispensable tool in physical chemistry, enzyme kinetics, atmospheric science, and computational chemistry.