Steady-State Approximation

The problem the approximation solves

Real reactions rarely happen in one step. A balanced equation such as 2 N₂O₅ → 4 NO₂ + O₂ hides a sequence of elementary steps that pass through fleeting species — atoms, radicals, ions, or excited molecules — that never accumulate in a flask. These are the reactive intermediates. Write the differential rate equations for a three-step mechanism and you get a coupled system of nonlinear ODEs with no closed-form solution. The chemist wants a single, tidy rate law in terms of species you can actually weigh out and measure. The steady-state approximation (SSA) is the trick that bridges that gap.

The physical insight is simple. An intermediate that is extremely reactive is destroyed almost the instant it appears. So after a short induction period, the rate at which it is made comes into near-balance with the rate at which it is destroyed. Its concentration climbs to a small plateau and then drifts only slowly. On the timescale of the bulk reaction, the slope d[I]/dt is negligible compared with the large, nearly equal formation and destruction fluxes that produce it. We approximate that slope as zero:

d[I]/dt = (rate of formation) − (rate of destruction) ≈ 0

This is an algebraic statement — it lets us solve for [I] without integrating anything. Note the word approximation: d[I]/dt is not truly zero, it is just tiny relative to the terms that make it up. For this reason it is more precisely called the quasi-steady-state or pseudo-steady-state approximation.

A worked mechanism, step by step

Consider the canonical two-step scheme where reactant A reversibly forms intermediate I, which then turns into product P:

A  ⇌  I   (k₁ forward, k₋₁ reverse)
I  →  P   (k₂)

The intermediate is formed by step 1 (rate k₁[A]) and destroyed two ways — by the reverse of step 1 (k₋₁[I]) and by step 2 (k₂[I]). The full balance is:

d[I]/dt = k₁[A] − k₋₁[I] − k₂[I]

Set it to zero and solve for the intermediate concentration:

[I] = k₁[A] / (k₋₁ + k₂)

The observed reaction rate is the rate of product formation, k₂[I]. Substituting:

rate = d[P]/dt = k₁k₂[A] / (k₋₁ + k₂)

The intermediate has vanished from the answer. Two limits are instructive. If the intermediate prefers to go forward to product (k₂ ≫ k₋₁), the denominator becomes k₂, the k₂ cancels, and rate = k₁[A] — the first step is rate-determining and the reaction is cleanly first-order. If the intermediate prefers to fall back to A (k₋₁ ≫ k₂), the rate becomes (k₁k₂/k₋₁)[A] = K·k₂[A], where K = k₁/k₋₁ is the equilibrium constant of step 1. This second limit is exactly the pre-equilibrium approximation — so pre-equilibrium is just the steady state in the regime where the reverse step dominates.

Lindemann: how a "unimolecular" reaction is really bimolecular

The first triumph of the SSA was explaining why gas-phase isomerizations and decompositions — like cyclopropane → propene — look first-order at high pressure but turn second-order at low pressure. Frederick Lindemann's 1922 mechanism (refined by Cyril Hinshelwood) says a molecule A is energized to A* by collision, and A* either is de-energized by another collision or reacts:

A + M → A* + M  (k₁)
A* + M → A + M  (k₋₁)
A* → P  (k₂)

Applying d[A*]/dt ≈ 0 gives [A*] = k₁[A][M]/(k₋₁[M] + k₂), and the rate becomes k₁k₂[A][M]/(k₋₁[M] + k₂). At high pressure (large [M]), the [M] terms dominate and rate = (k₁k₂/k₋₁)[A] — first-order. At low pressure (small [M]), the rate becomes k₁[A][M] — second-order, because energizing collisions become the bottleneck. The single steady-state expression captures the entire fall-off curve that bridges the two regimes. The crossover pressure for cyclopropane is around 1 to 100 Torr depending on temperature.

Enzymes: the Michaelis-Menten law

The most famous use of the SSA in biochemistry is the Briggs-Haldane (1925) derivation of enzyme kinetics. For E + S ⇌ ES → E + P, applying d[ES]/dt ≈ 0 yields [ES] = [E][S]/Kₘ, where the Michaelis constant Kₘ = (k₋₁ + k₂)/k₁. Combined with the conservation law [E]total = [E] + [ES], this gives the saturation curve:

v = Vmax[S] / (Kₘ + [S])

where Vmax = k₂[E]total. When [S] = Kₘ the velocity is exactly half-maximal — Kₘ values for real enzymes span roughly 10⁻⁶ to 10⁻² M. The turnover number k₂ (also written kcat) ranges from a few per second to the spectacular ~10⁶ s⁻¹ of carbonic anhydrase and ~4×10⁷ s⁻¹ of catalase. Note that Michaelis and Menten's original 1913 treatment assumed rapid equilibrium (k₋₁ ≫ k₂); Briggs and Haldane's steady-state version is more general and reduces to it in that limit.

Steady-state vs. pre-equilibrium vs. rate-determining step

Three approximations are routinely used to extract a rate law from a mechanism. They are not competitors so much as a nested hierarchy — the steady state is the most general, and the others are special cases or different framings of it.

Approximation	Core assumption	Gives [I] as	Valid when	Result for A ⇌ I → P
Steady-state (SSA)	d[I]/dt ≈ 0	k₁[A] / (k₋₁ + k₂)	I is highly reactive; [I] stays small after a short induction	rate = k₁k₂[A]/(k₋₁ + k₂)
Pre-equilibrium	Step 1 at equilibrium	(k₁/k₋₁)[A]	k₋₁ ≫ k₂ (intermediate falls back faster than it advances)	rate = (k₁k₂/k₋₁)[A]
Rate-determining step	One step far slower than all others	not needed — slow step sets rate	one step has the largest activation energy	rate = rate of slow step

The key relationships: pre-equilibrium is the SSA result in the limit k₋₁ ≫ k₂. The rate-determining-step idea is what the SSA reduces to when one step's rate constant dwarfs the rest. In practice you reach for the SSA whenever an intermediate is too reactive to maintain a real equilibrium, and you fall back on pre-equilibrium when the reverse of the first step is genuinely fast.

Radical chains and atmospheric chemistry

Bodenstein's original problem, the H₂ + Br₂ → 2 HBr chain, is the textbook radical example. Br atoms and H atoms are the intermediates; applying the SSA to both gives the experimentally observed rate law, rate = k[H₂][Br₂]^½ / (1 + k′[HBr]/[Br₂]) — including the curious half-order in Br₂ and the inhibition by product HBr that no single elementary step could explain. Steady-state radical concentrations in such chains are vanishingly small, often 10⁻⁸ to 10⁻¹¹ M, which is exactly why the approximation is so accurate: the intermediate cannot accumulate.

The same logic governs free-radical polymerization, where the steady-state radical concentration (set by balancing initiation against termination) gives the classic rate ∝ [monomer][initiator]^½, and the Chapman cycle for stratospheric ozone, where applying d[O]/dt ≈ 0 to oxygen atoms predicts the steady-state O₃ profile and lets you model how chlorine and NOₓ catalysts perturb it. In combustion, atmospheric modeling, and chemical-reactor design, the SSA is what makes otherwise intractable networks of hundreds of radical reactions numerically solvable — modern stiff-ODE solvers exploit it implicitly.

When the approximation breaks down

The SSA is an approximation, and it has two well-defined failure regions. During the brief induction period at the very start, [I] is still climbing from zero toward its plateau, so d[I]/dt is large and the approximation underestimates [I]. Near the end of reaction, when reactants are nearly exhausted, the formation flux collapses and the intermediate decays — again d[I]/dt is no longer negligible. Between these bookends, in the long quasi-steady regime, the SSA can be accurate to a fraction of a percent. The rule of thumb is that the intermediate's chemical lifetime (1/(k₋₁ + k₂)) must be far shorter than the timescale over which reactant concentrations change appreciably. When that separation of timescales is large, the approximation is excellent; when intermediate and reactant lifetimes are comparable, you must integrate the full system numerically.