Steady-State Approximation in Chemical Kinetics | Intermediate Concentration

The Michaelis-Menten mechanism describes enzyme-catalyzed reactions through two steps: enzyme-substrate complex formation (E + S ⇌ ES) and product formation (ES → P + E). The steady state approximation assumes constant enzyme-substrate complex concentration, leading to the relationship k1[E][S] = (k-1 + k2)[ES]. Total enzyme concentration equals free enzyme plus complex concentration, enabling derivation of the Michaelis-Menten equation.

The Michaelis-Menten equation (v = Vmax[S]/(Km + [S])) describes enzyme kinetics by relating reaction velocity to substrate concentration, where Vmax is the maximum velocity when all enzyme active sites are saturated, and Km is the substrate concentration at half-maximal velocity; the derivation assumes substrate concentration is much higher than enzyme concentration and that the rate of product formation (k2) is much smaller than the rate constants for complex formation (k1) and dissociation (k-1), using the steady-state approximation where the rate of enzyme-substrate complex formation equals its breakdown rate.

The Michaelis-Menten equation (V = Vmax[S]/(Km + [S])) is derived by applying the steady-state approximation to the enzyme-catalyzed reaction mechanism, where the enzyme-substrate complex concentration remains constant initially, allowing the rate of product formation to be expressed in terms of substrate concentration using the equilibrium between complex formation and dissociation rates.

This comprehensive section covers the foundational theory of enzyme kinetics. It begins with the basic mechanism where enzymes accelerate reactions by forming enzyme-substrate complexes that convert to products. The steady state approximation is introduced as the key assumption that the enzyme-substrate intermediate concentration remains constant, meaning formation equals loss rates. Through systematic mathematical derivation using rate equations and total enzyme concentration relationships, the Michaelis-Menten equation is derived: V₀ = (Vmax × [S])/(Km + [S]). Vmax represents maximum reaction velocity when all enzymes are saturated, while Km is defined as the substrate concentration yielding half-maximal velocity, providing a quantitative measure of enzyme-substrate affinity.

In 1913, Michaelis and Menten proposed a model to describe enzyme kinetics. The mechanism involves: (1) Enzyme (E) binds with substrate (S) to form an enzyme-substrate complex (ES) with rate constant k1 (forward) and k-1 (reverse). (2) The ES complex then converts to product (P) and releases the enzyme with rate constant k2 (also called kcat or catalytic constant). The overall reaction is: E + S ⇌ ES → E + P. This two-step mechanism forms the basis for understanding how enzymes catalyze reactions.

Chain reactions and polymerization are complex reactions characterized by initiation, propagation, and termination steps; the quasi-steady state approximation allows derivation of kinetic rate expressions by assuming highly reactive intermediates (like radicals) maintain constant concentrations, enabling simplification of complex reaction mechanisms to observable rate laws such as R = k[Acetaldehyde]^(3/2) for acetaldehyde decomposition and explaining the unimodal molecular weight distribution in polymerization reactions.

In multi-step reactions, the overall rate depends on the slowest step (rate-determining step). For polymerization, if initiation is slow, it controls the rate; if propagation is slow, it controls the rate. Using steady state approximation for radicals, the rate of polymerization equals the rate of propagation. Substituting the expression for radical concentration derived from steady state conditions yields: Rate = kp × [M] × (ki × [I])^(1/2) / kt^(1/2). The kinetic chain length (ν) represents the average number of monomer units added per initiated chain: ν = kp × [M] / ki × [I]. These equations provide quantitative tools for predicting and controlling polymer properties based on reaction conditions.

The quasi-steady state assumption states that radical concentration remains constant, meaning initiation rate equals termination rate. This allows calculation of radical concentration as [R•] = √(fkd[I]/kt). The polymerization rate equals propagation rate: Rp = kp[M][R•]. The kinetic chain length (ν) is the ratio of propagation to initiation rate. The number-average degree of polymerization depends on termination mechanism: DPn = 2ν for recombination, DPn = ν for disproportionation, and DPn = ν(1+f) when both mechanisms occur.

The kinetic chain length ν is defined as the average number of monomer units added per active center before termination. For free radical polymerization, ν = R_p/R_i, where R_p is the rate of propagation and R_i is the rate of initiation. Given R_p = 4 × 10⁻³ M/s, R_i = k₁[I] = (2 × 10⁻³ M⁻¹s⁻¹)(10⁻³ M) = 2 × 10⁻⁶ M/s, then ν = (4 × 10⁻³)/(2 × 10⁻⁶) = 2,000. This means each active center adds an average of 2,000 monomer units before termination.

The overall polymerization rate equals the propagation rate when initiation and termination are fast compared to propagation. The propagation rate depends on radical and monomer concentrations multiplied by the propagation rate constant. The steady-state approximation assumes constant radical concentration because radicals form during initiation and destroy during termination at equal rates: 2fkp[M][R•] = kt[R•]^2. Solving gives [R•] = √(fkdkp[I][M]/kt). The average degree of polymerization (DPn) equals propagation rate divided by termination rate: DPn = Rp/Rt = √(fkdkp[I][M]/kt). When chain transfer occurs, DPn decreases according to 1/DPn = 1/DPn_0 + Cm[M] + Cs[S] + Ci[I], where C represents transfer constants.

The general solution combines homogeneous and particular solutions: u(x) = u_h(x) + u_p(x). For the example, the homogeneous solution is e^(-αx) and a particular solution is sin(x), giving u(x) = e^(-αx) + sin(x). When α is very large, the exponential decays extremely rapidly, causing the solution to quickly transition to the sinusoidal steady-state. This rapid change over a short interval creates the defining characteristic of stiff problems, where numerical methods struggle to capture the solution accurately without impractically small step sizes.

The Quasi-Steady State Assumption (QSSA) is a classical method to remove stiffness by assuming fast-evolving species have zero net production rate, converting their differential equations to algebraic equations. This reduces stiff ODEs to non-stiff differential-algebraic equations (DAEs) while preserving accuracy for slow-evolving species. Despite being over 100 years old, QSSA remains widely used in modern turbulent combustion simulations on supercomputers. The approach enables faster explicit solvers instead of expensive stiff solvers, making it computationally efficient while maintaining accuracy for the species of primary interest.

The steady-state approximation is used when an intermediate is formed and consumed rapidly. It assumes the concentration of the intermediate remains constant (d[intermediate]/dt ≈ 0). For a mechanism with fast equilibrium followed by a slow step: (1) Set rate of formation = rate of consumption for the intermediate; (2) Solve for the intermediate concentration; (3) Substitute into the rate law for the slow step. This eliminates intermediates from the rate law.

Stiff equations have terms changing at vastly different rates, requiring very small time steps for stability. The Robertson chemical reaction equations are a classic example. NDSolve's StiffnessSwitching option handles stiff equations efficiently by switching between explicit and implicit methods.

Stiff differential equations are those where explicit numerical methods (like explicit Euler) fail to produce physically meaningful results because they require impractically small step sizes to maintain stability, often causing oscillations or unphysical values; implicit methods (like implicit Euler) are preferred for stiff problems because they automatically damp fast-decaying components without requiring extremely small step sizes, preserving physical properties like positivity and monotonic decay.