Steady-State Approximation

A short-lived intermediate can control an entire reaction, yet its concentration is almost impossible to measure directly. Chemists solve this puzzle with a simple trick. They assume the intermediate's concentration barely changes once the reaction gets going. This idea, the steady-state approximation, converts an intractable set of coupled differential equations into a solvable algebra problem. It is the single most useful tool for turning a proposed multistep mechanism into a testable, closed-form rate law. This article builds the approximation from the mechanism up. You will see the mathematical justification. You will see two full worked derivations. Then you will see where the method fails and what it sets up for enzyme kinetics and chain reactions.

Key takeaways
  • The steady-state approximation sets d[I]/dt≈0d[\text{I}]/dt \approx 0 for a reactive intermediate I after a brief induction period.
  • It converts coupled rate equations into simple algebraic ones, letting you solve for [I][\text{I}] in terms of reactants.
  • It is valid when the intermediate is highly reactive, so it is consumed almost as fast as it forms.
  • It underlies the derivations of Michaelis–Menten kinetics, the Lindemann mechanism, and radical chain reactions.
  • It gives the same rate law as the more restrictive pre-equilibrium approximation when the first step is fast and reversible.
  • Key idea: treat a fleeting intermediate's concentration as constant, not zero, because its formation and destruction rates nearly balance.
The core insight

When an intermediate is consumed almost as quickly as it is made, its concentration rises briefly, then stays low and nearly flat, so we may set its net rate of change to zero and solve the resulting algebraic equations for the observable rate law.

The physical picture

Consider a mechanism that proceeds through an intermediate I, formed from reactant A and destroyed to give product P:

A→k1I→k2P\text{A} \xrightarrow{k_1} \text{I} \xrightarrow{k_2} \text{P}

If k2≫k1k_2 \gg k_1, then as soon as a molecule of I is born it is converted onward almost instantly. The population of I builds up briefly during an induction period, then levels off at a small, nearly constant value while A is depleted and P accumulates. During this long "steady" stretch, the rate at which I is created by step 1 is essentially equal to the rate at which it is destroyed by step 2. The net rate of change of [I][\text{I}] is therefore very small compared with the rates of change of [A][\text{A}] and [P][\text{P}], even though [I][\text{I}] itself is not zero. This is fundamentally different from assuming [I]=0[\text{I}] = 0. We keep [I][\text{I}] as a small, dynamically-maintained quantity and use the condition d[I]/dt≈0d[\text{I}]/dt \approx 0 as an algebraic constraint that lets us eliminate [I][\text{I}] from the rate law entirely. The result is an expression for the reaction rate written purely in terms of species we can actually monitor by experiment — reactants, products, and sometimes a stable catalyst. The approximation is a close cousin of the pre-equilibrium approximation encountered in consecutive reactions, but it is more general: pre-equilibrium requires the first step to equilibrate rapidly and reversibly, while the steady-state method only requires the intermediate to be short-lived, which is true far more often.

Figure
Concentration profiles for A, the fleeting intermediate I, and product P versus time
time, t concentration A I (small, nearly flat plateau) P induction period ends

The core relations

d[I]dt≈0\frac{d[\text{I}]}{dt} \approx 0

where: this holds only after the brief induction period, while I remains at a small, slowly-drifting concentration.

rateformation of I≈rateconsumption of I\text{rate}_{\text{formation of I}} \approx \text{rate}_{\text{consumption of I}}

where: this balance is what "steady state" means physically — it is not that nothing happens to I, but that its net change is negligible.

[I]ss=k1[A]k2[\text{I}]_{ss} = \frac{k_1[\text{A}]}{k_2}

where: for the simple two-step chain above, this is the algebraic solution obtained by setting d[I]/dt=0d[\text{I}]/dt = 0.

rate=k2[I]ss=k1[A]\text{rate} = k_2[\text{I}]_{ss} = k_1[\text{A}]

where: the observable rate becomes first order in A alone, with the slow step (formation of I) controlling the overall rate.

Derivation ASteady-state rate law for A → I → P

Goal. Derive the overall rate law for the mechanism

A→k1I→k2P\text{A} \xrightarrow{k_1} \text{I} \xrightarrow{k_2} \text{P}

without assuming either step is at equilibrium, using only the steady-state condition on I.

Step 1. Write the exact rate equations for every species.

d[A]dt=−k1[A]\frac{d[\text{A}]}{dt} = -k_1[\text{A}]
d[I]dt=k1[A]−k2[I]\frac{d[\text{I}]}{dt} = k_1[\text{A}] - k_2[\text{I}]
d[P]dt=k2[I]\frac{d[\text{P}]}{dt} = k_2[\text{I}]

These three coupled first-order ODEs are exactly solvable here (it is a classic consecutive-reaction problem), but in general multistep mechanisms give ODEs with no closed form. The steady-state method sidesteps that difficulty.

Step 2. Apply the steady-state condition to the intermediate. Assume k2≫k1k_2 \gg k_1, so I is consumed almost as fast as it forms. After a short induction period,

d[I]dt≈0\frac{d[\text{I}]}{dt} \approx 0

This converts the differential equation for I into an algebraic one:

0=k1[A]−k2[I]ss0 = k_1[\text{A}] - k_2[\text{I}]_{ss}

Step 3. Solve algebraically for the steady-state intermediate concentration.

[I]ss=k1k2[A][\text{I}]_{ss} = \frac{k_1}{k_2}[\text{A}]

Because k2≫k1k_2 \gg k_1, the ratio k1/k2≪1k_1/k_2 \ll 1, confirming self-consistently that [I]ss[\text{I}]_{ss} stays small relative to [A][\text{A}] throughout the reaction — this is the internal consistency check every steady-state solution should pass.

Step 4. Substitute back into the rate of product formation.

d[P]dt=k2[I]ss=k2(k1k2[A])=k1[A]\frac{d[\text{P}]}{dt} = k_2[\text{I}]_{ss} = k_2 \left(\frac{k_1}{k_2}[\text{A}]\right) = k_1[\text{A}]

Step 5. Compare with the exact solution. The exact solution of the linear ODE system gives

[P](t)=[A]0(1−k2e−k1t−k1e−k2tk2−k1)[\text{P}](t) = [\text{A}]_0\left(1 - \frac{k_2 e^{-k_1 t} - k_1 e^{-k_2 t}}{k_2 - k_1}\right)

When k2≫k1k_2 \gg k_1, the term e−k2te^{-k_2 t} decays essentially instantly (that is the induction period), leaving

[P](t)≈[A]0(1−e−k1t)[\text{P}](t) \approx [\text{A}]_0\left(1 - e^{-k_1 t}\right)

which is exactly what integrating d[P]/dt=k1[A]d[\text{P}]/dt = k_1[\text{A}] predicts. The steady-state shortcut reproduces the exact long-time behavior.

Boxed result.

rate=k1[A],[I]ss=k1k2[A]\boxed{\text{rate} = k_1[\text{A}], \qquad [\text{I}]_{ss} = \frac{k_1}{k_2}[\text{A}]}

The overall rate law depends only on the slow rate constant k1k_1 — the fast second step never appears in the final observable rate, even though it is mechanistically essential.

Worked examples

Example 1 — Foundational: gas-phase decomposition via an energetic intermediate. Consider the mechanism for the thermal decomposition of a hypothetical gas A₂B into products, going through intermediate I:

A2B→k1I+X,I→k2Y+Z\text{A}_2\text{B} \xrightarrow{k_1} \text{I} + \text{X}, \qquad \text{I} \xrightarrow{k_2} \text{Y} + \text{Z}

with k1=2.5×10−3 s−1k_1 = 2.5 \times 10^{-3}\ \text{s}^{-1} and k2=4.8×104 s−1k_2 = 4.8 \times 10^{4}\ \text{s}^{-1} at 500 K, starting from [A2B]0=0.100 mol dm−3[\text{A}_2\text{B}]_0 = 0.100\ \text{mol dm}^{-3}. Step 1: check the steady-state condition is reasonable.

k1k2=2.5×10−34.8×104=5.2×10−8\frac{k_1}{k_2} = \frac{2.5\times10^{-3}}{4.8\times10^{4}} = 5.2\times10^{-8}

Since k2/k1≈1.9×107k_2/k_1 \approx 1.9\times10^7, I is consumed roughly 19 million times faster than it forms — an excellent steady-state candidate. Step 2: steady-state intermediate concentration at t=0t=0 (when [A2B]≈0.100 mol dm−3[\text{A}_2\text{B}] \approx 0.100\ \text{mol dm}^{-3}).

[I]ss=k1k2[A2B]=(5.2×10−8)(0.100 mol dm−3)=5.2×10−9 mol dm−3[\text{I}]_{ss} = \frac{k_1}{k_2}[\text{A}_2\text{B}] = (5.2\times10^{-8})(0.100\ \text{mol dm}^{-3}) = 5.2\times10^{-9}\ \text{mol dm}^{-3}

This is indeed negligible compared with [A2B][\text{A}_2\text{B}], confirming self-consistency. Step 3: the overall rate.

rate=k1[A2B]=(2.5×10−3 s−1)(0.100 mol dm−3)=2.5×10−4 mol dm−3s−1\text{rate} = k_1[\text{A}_2\text{B}] = (2.5\times10^{-3}\ \text{s}^{-1})(0.100\ \text{mol dm}^{-3}) = 2.5\times10^{-4}\ \text{mol dm}^{-3}\text{s}^{-1}

Sanity check: the reaction is first order overall, with a half-life

t1/2=ln⁡2k1=0.6932.5×10−3 s−1=277 st_{1/2} = \frac{\ln 2}{k_1} = \frac{0.693}{2.5\times10^{-3}\ \text{s}^{-1}} = 277\ \text{s}

comfortably longer than the induction time for I, which scales as 1/k2≈2.1×10−5 s1/k_2 \approx 2.1\times10^{-5}\ \text{s} — four orders of magnitude shorter, so the steady-state plateau is established almost immediately on the timescale of the overall reaction.

Example 2 — Applied: NO-catalyzed ozone decomposition (three-step radical mechanism). A classic atmospheric mechanism for ozone loss catalyzed by NO is:

NO+O3→k1NO2+O2(slow-ish)\text{NO} + \text{O}_3 \xrightarrow{k_1} \text{NO}_2 + \text{O}_2 \qquad (\text{slow-ish})
NO2+O→k2NO+O2(fast)\text{NO}_2 + \text{O} \xrightarrow{k_2} \text{NO} + \text{O}_2 \qquad (\text{fast})

with atomic oxygen O supplied by a background photolysis step

O3→k3O2+O\text{O}_3 \xrightarrow{k_3} \text{O}_2 + \text{O}

Here NO₂ is the reactive intermediate (NO and O are catalyst/co-reactant, regenerated or supplied externally). We treat NO₂ with the steady-state approximation. Step 1: write d[NO2]/dtd[\text{NO}_2]/dt.

d[NO2]dt=k1[NO][O3]−k2[NO2][O]\frac{d[\text{NO}_2]}{dt} = k_1[\text{NO}][\text{O}_3] - k_2[\text{NO}_2][\text{O}]

Step 2: apply steady state.

0=k1[NO][O3]−k2[NO2]ss[O]0 = k_1[\text{NO}][\text{O}_3] - k_2[\text{NO}_2]_{ss}[\text{O}]
[NO2]ss=k1[NO][O3]k2[O][\text{NO}_2]_{ss} = \frac{k_1[\text{NO}][\text{O}_3]}{k_2[\text{O}]}

Step 3: rate of ozone loss. Ozone is consumed in step 1 directly, and O is consumed in step 2, which regenerates NO. The rate-determining chemistry for O₃ loss is

−d[O3]dt=k1[NO][O3]-\frac{d[\text{O}_3]}{dt} = k_1[\text{NO}][\text{O}_3]

Because the NO consumed in step 1 is exactly regenerated in step 2 (it is a catalyst), [NO][\text{NO}] stays constant, and the observed rate is first order in both NO and O₃. Numeric plug-in. At 220 K (stratospheric conditions), literature-typical values are k1=1.8×10−14 cm3molecule−1s−1k_1 = 1.8\times10^{-14}\ \text{cm}^3\text{molecule}^{-1}\text{s}^{-1}, with typical stratospheric concentrations [NO]=5.0×108 molecule cm−3[\text{NO}] = 5.0\times10^{8}\ \text{molecule cm}^{-3} and [O3]=3.0×1012 molecule cm−3[\text{O}_3] = 3.0\times10^{12}\ \text{molecule cm}^{-3}.

−d[O3]dt=(1.8×10−14)(5.0×108)(3.0×1012) molecule cm−3s−1-\frac{d[\text{O}_3]}{dt} = (1.8\times10^{-14})(5.0\times10^{8})(3.0\times10^{12})\ \text{molecule cm}^{-3}\text{s}^{-1}
=2.7×107 molecule cm−3s−1= 2.7\times10^{7}\ \text{molecule cm}^{-3}\text{s}^{-1}

Sanity check: this is a bimolecular rate law, exactly the form measured experimentally for this catalytic cycle, and the derivation never required knowing the absolute concentration of the NO₂ intermediate — only its steady-state ratio to NO and O₃, which is the entire point of the method.

Limiting cases & common traps

SituationWhat to do
k2≫k1k_2 \gg k_1 (intermediate very reactive)Steady-state approximation is excellent; induction period is negligible on the reaction timescale
k2∼k1k_2 \sim k_1 (comparable rates)Steady-state approximation breaks down; must solve exact coupled ODEs
First step reversible and fast, second step slowPre-equilibrium approximation applies and typically agrees with steady-state result
Multiple intermediates in seriesApply d[Ii]/dt≈0d[\text{I}_i]/dt \approx 0 separately to each fast-decaying intermediate
Very early times (within the induction period)Steady-state result is invalid; [I][\text{I}] is still rising, not yet plateaued

The chemist behind it

The steady-state approximation grew out of early 20th-century efforts to explain gas-phase reaction orders that didn't match naive collision theory — most famously the puzzling kinetics of unimolecular gas reactions. Frederick Alexander Lindemann proposed in 1922 that an energized molecule must first be produced by collision before it can react unimolecularly, introducing exactly the kind of transient intermediate this approximation was built to handle. Around the same time, Irving Langmuir and others applied similar reasoning to catalytic surface reactions, while chemists studying free-radical chain reactions — such as the hydrogen–bromine reaction — found that assuming steady concentrations of radical intermediates was the only way to reproduce the observed complicated rate law. The method became a standard toolkit item once Michaelis and Menten's 1913 enzyme kinetics work was reinterpreted through this lens by Briggs and Haldane in 1925, who explicitly invoked the steady-state condition on the enzyme–substrate complex rather than Michaelis and Menten's original equilibrium assumption.

Origins · 1 minGeorge Edward Briggs and John Burdon Sanderson Haldane and how this was discovered
GE
George Edward Briggs and John Burdon Sanderson Haldane
1893–1985 (Briggs); 1892–1964 (Haldane) · reformulated Michaelis–Menten enzyme kinetics using the steady-state approximation in 1925

:::quiz[In the derivation of A → I → P, what does setting d[I]/dt ≈ 0 allow you to do?]

  • Ignore the second step entirely
  • Convert a differential equation for [I] into an algebraic equation, solvable for [I]_ss ✓
  • Prove that k1 equals k2
  • Replace the rate law with the equilibrium constant expression

:::faq[Why can we set d[I]/dt to zero if the intermediate's concentration is changing?] It isn't truly zero — it's small compared with the rates of change of reactants and products. After a brief induction period, formation and consumption of I nearly balance, so [I][\text{I}] drifts very slowly compared with the timescale of the overall reaction. Treating d[I]/dtd[\text{I}]/dt as zero is an approximation valid on the timescale we care about, not a claim that nothing happens to I at the microscopic level.

Ask the tutor: Want to see how this exact steady-state logic turns into the Michaelis–Menten equation for enzyme catalysis?

How do I know if the steady-state approximation applies to a given mechanism?

Check whether the intermediate is consumed by a much faster step than the one that forms it — ideally kconsumption≫kformationk_{\text{consumption}} \gg k_{\text{formation}} by at least one or two orders of magnitude. Reactive species like free radicals, excited states, or unstable complexes are usually excellent candidates. If the rate constants are comparable, you must solve the exact coupled differential equations instead.

Does the steady-state approximation always give the same rate law as the pre-equilibrium approximation?

When the first step really is fast and reversible with a much slower second step, both methods converge on the same rate law, because pre-equilibrium is a special case of steady state. But the steady-state approximation is more broadly applicable — it works even when the first step is not reversible, as long as the intermediate is short-lived, so it is the more general and more commonly used tool.

Where does the steady-state approximation show up later in kinetics?

It is the mathematical backbone of Michaelis–Menten enzyme kinetics, where the enzyme–substrate complex is the steady-state intermediate. It also underlies the Lindemann–Hinshelwood mechanism for unimolecular gas reactions, and it is essential for deriving rate laws in radical chain reactions, where multiple reactive intermediates must all be treated as steady-state species simultaneously.

Before this

Where this leads