Multistep Reaction Energy Profile

SECTION 1

Historical Context & Motivation

The notion that chemical reactions proceed through a single concerted step dominated early mechanistic thinking. By the late nineteenth century, however, observations of complex reaction kinetics—rate laws that could not be explained by a single elementary step—demanded a more nuanced picture. The development of multistep reaction energy profiles arose from the confluence of thermodynamics, transition state theory, and experimental kinetics, providing chemists with a visual framework to map the energetic landscape of reactions that pass through one or more reactive intermediates before arriving at products.

1889

Arrhenius Equation

Svante Arrhenius proposed his celebrated equation relating rate constants to temperature, introducing the concept of activation energy (E_a) as an energetic barrier that reactants must overcome.

1935

Transition State Theory

Henry Eyring, Meredith Gwynne Evans, and Michael Polanyi independently formulated transition state theory (TST), placing the activated complex on a potential energy surface and enabling rigorous construction of energy profiles for elementary steps.

1944

LFER & Hammond's Postulate Foundations

The development of linear free-energy relationships (LFER) connected thermodynamic stability of intermediates to kinetic barrier heights, enabling predictive energy profile sketching for multistep mechanisms.

1955

Hammond's Postulate

George Hammond articulated his postulate: for an exothermic elementary step the transition state resembles the reactants, while for an endothermic step it resembles the products. This principle became central to interpreting the shapes of multistep energy diagrams.

1968

Computational Potential Energy Surfaces

Advances in quantum chemistry enabled Roald Hoffmann and others to compute full potential energy surfaces, confirming experimentally inferred multistep profiles and elucidating previously undetectable intermediates.

The central question that multistep energy profiles address is deceptively simple: when a reaction proceeds through several elementary steps, which step controls the overall rate, and how do the energies of intermediates relate to the transition states that connect them? Answering this question is essential for catalyst design, pharmaceutical synthesis, and understanding biological enzyme mechanisms.

SECTION 2

Core Principles & Definitions

A multistep energy profile plots the free energy (or potential energy) of a reacting system along the reaction coordinate, which is an abstract axis representing the progress from reactants through intermediates to products. Unlike a single-step profile with one maximum, a multistep profile features multiple peaks and valleys—each peak corresponding to a transition state and each valley to a reactive intermediate. Understanding these features requires internalizing several foundational concepts.

Transition State (TS)

A saddle point on the potential energy surface representing the highest-energy configuration along the minimum-energy pathway for a single elementary step. It cannot be isolated and exists for roughly 10⁻¹³ s. Denoted by the ‡ (double-dagger) symbol.

Reactive Intermediate

A species that lies in a local energy minimum between two transition states. Unlike a TS, an intermediate has a finite (though often short) lifetime and can sometimes be detected spectroscopically. Examples include carbocations, carbanions, free radicals, and enzyme–substrate complexes.

Activation Energy (Eₐ)

The energy difference between the reactant (or intermediate) and the immediately following transition state for a given elementary step. Each step in a multistep mechanism has its own Eₐ, and the step with the largest Eₐ is typically the rate-determining step.

Rate-Determining Step (RDS)

The elementary step with the highest-energy transition state relative to the starting materials (or the highest individual activation energy). It acts as the kinetic bottleneck, governing the overall reaction rate and determining the form of the experimental rate law.

Overall ΔG (or ΔH)

The free energy difference between the final products and the initial reactants. This thermodynamic quantity is path-independent—it does not depend on the number of steps—and determines whether the overall process is spontaneous (ΔG < 0) or nonspontaneous (ΔG > 0).

✦ KEY TAKEAWAY

Think of a multistep energy profile like a mountain hike with several peaks and valleys. The transition states are the mountain passes you must cross, and the intermediates are the rest stops in valleys between passes. The rate-determining step is the highest pass—no matter how easy the other passes are, the speed at which your group completes the hike is limited by the time it takes to cross that tallest ridge. Meanwhile, the overall elevation change from trailhead to destination tells you the thermodynamic favorability (ΔG), which is independent of the route taken.

SECTION 3

Visual Explanation — The Two-Step Energy Profile

The following diagram illustrates a generic two-step exothermic reaction in which reactants (A + B) form an intermediate (I) that then converts to products (C + D). Two transition states (TS₁ and TS₂) separate the three energy minima. Pay close attention to the relative heights of the two barriers: the taller barrier identifies the rate-determining step.

Figure 1. A two-step exothermic energy profile. The reactants (A + B) climb over TS₁ to form intermediate I, which then crosses TS₂ to yield products (C + D). The dashed vertical arrows indicate the activation energies for each step (E_a,1 and E_a,2). Because TS₂ is the global maximum, step 2 is the rate-determining step.

Several features deserve attention. First, intermediate I sits in a local minimum—it is a real, albeit transient, species with a definite geometry and electronic structure. Second, the rate-determining step is identified not merely by comparing E_a,1 and E_a,2 directly, but by determining which transition state has the highest absolute free energy relative to the starting materials. Third, the overall thermodynamic driving force (ΔG < 0) is independent of the barrier heights—a reaction can be highly exergonic yet exceedingly slow if its activation barriers are large.

SECTION 4

Mathematical Framework

Quantifying a multistep energy profile requires connecting the macroscopic rate constant to the microscopic barrier heights through transition state theory and the steady-state or pre-equilibrium approximations. The Eyring equation provides the bridge between the free energy of activation for each step and the experimentally measured rate constant.

EYRING EQUATION

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

where k is the rate constant for an elementary step, k_B is Boltzmann's constant (1.381 × 10⁻²³ J K⁻¹), T is temperature in kelvin, h is Planck's constant (6.626 × 10⁻³⁴ J s), ΔG‡ is the Gibbs free energy of activation, and R is the gas constant (8.314 J mol⁻¹ K⁻¹).

For a two-step mechanism A → I → P, if the first step is fast and reversible (the pre-equilibrium approximation), the equilibrium between A and I is established quickly relative to the slower second step. The equilibrium constant K_eq for the first step is related to the free energy change of that step.

PRE-EQUILIBRIUM RELATION

K_eq = k₁ / k₋₁ = exp(−ΔG₁° / RT)

K_eq is the equilibrium constant for step 1, k₁ and k₋₁ are the forward and reverse rate constants, and ΔG₁° is the standard free energy change for step 1.

OVERALL OBSERVED RATE CONSTANT

k_obs = K_eq × k₂ = (k₁ / k₋₁) × k₂

The observed rate constant k_obs combines the pre-equilibrium constant with the rate constant of the slow (rate-determining) second step. This shows why the energy of both the intermediate and the highest transition state matter for the overall rate.

ARRHENIUS EQUATION

k = A × exp(−Eₐ / RT)

The Arrhenius form relates the rate constant to the activation energy E_a and the pre-exponential factor A. An Arrhenius plot (ln k vs. 1/T) yields a straight line whose slope equals −E_a/R, providing experimental access to individual step activation energies when rate constants for each step can be isolated.

⚗️ Steady-State vs. Pre-Equilibrium

When neither step is overwhelmingly fast or slow, the steady-state approximation (d[I]/dt ≈ 0) is preferred over the pre-equilibrium approach. Both approximations ultimately identify the rate-determining step—the one whose transition state sits highest on the energy profile—but they differ in their mathematical treatment of the intermediate concentration.

SECTION 5

Catalyzed vs. Uncatalyzed Energy Profiles

One of the most powerful applications of multistep energy profiles is visualizing the effect of a catalyst. A catalyst provides an alternative reaction pathway with a lower overall activation energy barrier. Crucially, the catalyst does not alter the overall ΔG of the reaction—it lowers the kinetic barrier without changing the thermodynamic outcome. The diagram below contrasts the single-step uncatalyzed pathway with a two-step catalyzed pathway for the same overall transformation.

Figure 2. Comparison of uncatalyzed (dashed red) and catalyzed (solid green) pathways. The catalyst introduces an intermediate and splits the single high barrier into two lower barriers (TS₁' and TS₂'). The overall ΔG remains unchanged.

Notice that even though the catalyzed pathway has two transition states, both TS₁' and TS₂' lie below the single uncatalyzed TS. The overall activation energy for the catalyzed process (measured from reactants to the highest catalyzed TS) is substantially lower, resulting in a dramatically larger rate constant. Enzymes, heterogeneous catalysts, and organocatalysts all operate on this principle—they stabilize transition states differentially, not products. This is why the energy profile diagram is the single most useful tool for understanding catalytic rate enhancement.

⚠️ Common Misconception

A catalyst does not lower the energy of the products or change ΔG. It also does not shift the equilibrium position. It accelerates both the forward and reverse reactions equally by lowering the activation barrier of the rate-determining step (or all steps along the pathway).

SECTION 6

Worked Example — S_N1 Hydrolysis of tert-Butyl Bromide

Consider the S_N1 hydrolysis of tert-butyl bromide ((CH₃)₃CBr) in aqueous solution. The mechanism proceeds in two steps: (1) rate-limiting ionization to form a carbocation intermediate, and (2) fast nucleophilic attack by water. We will construct the energy profile and estimate the overall activation energy.

Constructing the Energy Profile for S_N1 Hydrolysis of (CH₃)₃CBr

Step 1 — Identify the Elementary StepsStep 1 (slow): (CH₃)₃CBr → (CH₃)₃C⁺ + Br⁻. This ionization is endothermic and has a large activation energy (E_a,1 ≈ 105 kJ mol⁻¹). Step 2 (fast): (CH₃)₃C⁺ + H₂O → (CH₃)₃COH₂⁺ → (CH₃)₃COH + H⁺. This step has a much smaller barrier (E_a,2 ≈ 15 kJ mol⁻¹).

Two elementary steps identified: ionization (RDS) and nucleophilic capture.

Step 2 — Place Reactants and Products on the Energy AxisThe reactant (CH₃)₃CBr is assigned a reference energy of 0 kJ mol⁻¹. The overall reaction is exergonic with ΔG° ≈ −30 kJ mol⁻¹, so the product (CH₃)₃COH sits 30 kJ mol⁻¹ below the reactant level.

Reactant at 0 kJ mol⁻¹; product at −30 kJ mol⁻¹.

Step 3 — Place the IntermediateThe carbocation intermediate (CH₃)₃C⁺ is higher in energy than the reactant because ionization is endothermic. Typical ΔG° for this ionization in water is approximately +70 kJ mol⁻¹, placing the intermediate at +70 kJ mol⁻¹.

Intermediate at +70 kJ mol⁻¹.

Step 4 — Place the Transition StatesTS₁ lies E_a,1 = 105 kJ mol⁻¹ above the reactant, so TS₁ is at +105 kJ mol⁻¹. TS₂ lies E_a,2 = 15 kJ mol⁻¹ above the intermediate, so TS₂ is at 70 + 15 = +85 kJ mol⁻¹.

TS₁ at +105 kJ mol⁻¹ (global max → RDS); TS₂ at +85 kJ mol⁻¹.

Step 5 — Identify the Rate-Determining StepTS₁ (+105 kJ mol⁻¹) is the highest point on the entire profile. Therefore, step 1 (ionization) is the rate-determining step. The overall activation energy for the reaction equals E_a,1 = 105 kJ mol⁻¹. The rate law predicted is: Rate = k[(CH₃)₃CBr], which is first order in substrate and zero order in nucleophile—consistent with experimental observations.

Overall Eₐ = 105 kJ mol⁻¹; Rate = k[(CH₃)₃CBr]

SECTION 7

Strengths, Limitations & Common Pitfalls

Strengths and limitations of multistep energy profiles

Feature	Strengths	Limitations
Visual Clarity	Makes abstract kinetic and thermodynamic concepts tangible; clearly shows which step controls the rate.	Reduces a multidimensional potential energy surface to a one-dimensional projection—nuances of molecular motion are lost.
RDS Identification	Directly reveals the rate-determining step by the height of the global maximum.	In some multistep mechanisms, multiple steps may have similar barrier heights, making a single RDS assignment an oversimplification.
Catalyst Design	Guides catalyst optimization by showing precisely which TS to stabilize.	Does not reveal the structural details of the TS—computational chemistry or isotope-effect studies are needed for that.
Intermediate Stability	Predicts whether intermediates might accumulate and be detectable.	A shallow minimum may not correspond to a kinetically significant intermediate in solution; solvent effects are often approximated.
Quantitative Use	Barrier heights can be extracted from Arrhenius or Eyring plots.	Accuracy depends on the assumption that TST is valid (breaks down for tunneling, recrossing, or very low barriers).

✦ KEY TAKEAWAY

A multistep energy profile is like a topographical cross-section of a hiking trail: it faithfully shows elevation changes along the path you've chosen, but it does not reveal alternative routes over adjacent ridges or the full three-dimensional terrain. For simple to moderately complex mechanisms, the one-dimensional projection is extraordinarily useful; for complex enzyme catalysis or atmospheric reactions with dozens of pathways, a full potential energy surface computed ab initio is required.

SECTION 8

Connection to Advanced Theory

The one-dimensional energy profile taught in introductory kinetics is a slice through a much richer mathematical object—the multidimensional potential energy surface (PES). In advanced physical chemistry and computational chemistry courses, you will encounter the full PES, intrinsic reaction coordinates (IRC), and Marcus theory for electron transfer reactions. The table below maps familiar concepts from this lesson to their advanced counterparts.

Introductory concepts and their advanced counterparts

Introductory Concept	Advanced Extension
Reaction coordinate (1-D)	Intrinsic Reaction Coordinate (IRC) on a multidimensional PES — the minimum-energy path connecting reactants, TS, and products in 3N−6 dimensional space.
Transition state as a single point	Transition state as a first-order saddle point: a maximum in the reaction coordinate direction and minimum in all other directions.
Activation energy (Arrhenius Eₐ)	ΔG‡ decomposed into ΔH‡ and −TΔS‡ via Eyring analysis; temperature dependence of Eₐ exposed.
Rate-determining step (single bottleneck)	Energetic span model (Kozuch & Shaik): defines the TOF-determining intermediate (TDI) and TOF-determining transition state (TDTS) for catalytic cycles.
Hammond's Postulate (qualitative)	Marcus theory: quantitative parabolic intersection model relating ΔG° to ΔG‡ through the reorganization energy λ.

As you progress to graduate-level kinetics, you will encounter variational transition state theory (VTST), which allows the dividing surface to be optimized rather than fixed at the saddle point, and Rice–Ramsperger–Kassel–Marcus (RRKM) theory for unimolecular reactions, which explicitly treats the flow of energy among internal degrees of freedom. These theories enrich and sometimes correct the simplified one-dimensional picture, but the intuition built from multistep energy profiles remains the indispensable starting point for mechanistic reasoning.

SECTION 9

Practice Problems

PROBLEM 1 — CONCEPTUAL

A reaction proceeds through two elementary steps. Step 1 has E_a,1 = 60 kJ mol⁻¹ and ΔG₁° = +40 kJ mol⁻¹. Step 2 has E_a,2 = 50 kJ mol⁻¹. Which step is the rate-determining step, and why can you not answer this question by simply comparing E_a,1 and E_a,2 alone?

PROBLEM 2 — BASIC CALCULATION

Given a two-step mechanism where the first step has ΔG‡₁ = 85 kJ mol⁻¹ at 298 K, calculate the rate constant k₁ using the Eyring equation. (k_B = 1.381 × 10⁻²³ J K⁻¹, h = 6.626 × 10⁻³⁴ J s, R = 8.314 J mol⁻¹ K⁻¹)

PROBLEM 3 — INTERMEDIATE

A three-step mechanism has the following energies (kJ mol⁻¹) relative to reactants: TS₁ = +75, Int₁ = +30, TS₂ = +95, Int₂ = +10, TS₃ = +55, Products = −40. (a) Identify the RDS. (b) Calculate the activation energy for each elementary step. (c) Determine the overall ΔG for the reaction.

PROBLEM 4 — APPLIED

An enzyme catalyzes a reaction by lowering the activation energy of the rate-determining step from 90 kJ mol⁻¹ to 45 kJ mol⁻¹ at 310 K (body temperature). By what factor does the rate constant of the RDS increase? Use the Arrhenius equation and assume the pre-exponential factor A remains unchanged.

PROBLEM 5 — CRITICAL THINKING

Consider a two-step catalytic cycle where the catalyst is regenerated. The energetic span model defines the turnover frequency (TOF) as being determined by the energy difference between the TOF-determining transition state (TDTS) and the TOF-determining intermediate (TDI). For a cycle with Int₁ at −10 kJ mol⁻¹, TS₁ at +60 kJ mol⁻¹, Int₂ at −25 kJ mol⁻¹, and TS₂ at +45 kJ mol⁻¹ (all relative to the catalyst resting state), with an overall ΔG_rxn = −50 kJ mol⁻¹, identify the TDTS and TDI and calculate the energetic span δE.

SUMMARY

Summary

A multistep reaction energy profile plots free energy against the reaction coordinate for mechanisms involving two or more elementary steps. Each energy maximum represents a transition state (‡), while each local minimum between maxima corresponds to a reactive intermediate. The rate-determining step is the one whose transition state has the highest absolute energy relative to the starting materials—not simply the step with the largest individual activation energy. The overall thermodynamic driving force (ΔG) is path-independent and equals the energy difference between products and reactants.

Quantitatively, each elementary rate constant is related to its barrier height through the Eyring equation or the Arrhenius equation. A catalyst provides an alternative pathway with a lower overall barrier but does not change ΔG. The pre-equilibrium and steady-state approximations connect the energy profile to the observed rate law. In advanced work, the one-dimensional energy profile extends to the full multidimensional potential energy surface analyzed by computational methods such as IRC calculations, Marcus theory, and the energetic span model for catalytic cycles.