Tutorials › Quantitative Systems Pharmacology (QSP) › Chemical Equilibrium Versus Steady State

Reaction Kinetics and Mass Action · Tutorial 38 of 799

Chemical Equilibrium Versus Steady State

Learn how to identify equilibrium and steady state from reaction fluxes and concentration changes, and why constant concentration does not necessarily mean that reactions have stopped.

Beginner 10 min read

What You'll Learn

  • Define equilibrium and steady state using reaction fluxes and concentration changes.
  • Derive the equilibrium condition for a reversible first-order reaction.
  • Calculate a driven steady-state concentration and its continuing flux.
  • Recognize why a momentary zero rate of change is not necessarily steady state.
  • Compare equilibrium with steady state and identify common diagnostic errors.

Two Different Questions About a Reaction System

A reaction system can have unchanging concentrations while reactions continue. For example, an intermediate may be produced at the same rate that it is consumed. Conversely, a system may have no net reaction flux only after it has reached a particular composition. These statements describe different features: one concerns reaction flux, and the other concerns how concentrations change over time.

This distinction matters when interpreting a QSP model. A flat concentration-time profile does not by itself show that the system is at chemical equilibrium. It may instead reflect a balance of ongoing processes, maintained by feeds, sinks, or other reactions. To decide which interpretation is appropriate, inspect both the concentration derivatives and the reaction fluxes, as well as the conditions imposed on the system.

Definition: A system is at steady state with respect to a set of species when the concentrations of those species are constant in time under the specified conditions: \(d[X_i]/dt=0\) for each species \(X_i\) in the set. A system is at chemical equilibrium when it is in a thermodynamic equilibrium state, with no sustained net reaction progress or other driving force. For an elementary reversible reaction at equilibrium, its forward and reverse fluxes are equal, so its net reaction flux is zero.

A reaction flux is a rate of reaction progress, commonly expressed as concentration per unit time. For a reversible step, the net flux is the forward flux minus the reverse flux. Equal forward and reverse fluxes do not mean that molecular events stop: molecules continue to react in both directions, but the two directions balance macroscopically. A steady state, by contrast, is defined by the net rates of change of concentrations. Those derivatives can be zero even when individual reaction fluxes are positive.

When Equilibrium and Steady State Coincide

Consider a closed, well-mixed system of constant volume containing a single reversible first-order reaction, \(A \rightleftharpoons B\). Suppose its forward and reverse fluxes follow mass action, with rate constants \(k_f\) and \(k_r\). There are no external feeds or removals, and the total concentration of \(A\) and \(B\) is conserved. The net flux in the forward direction is the forward flux minus the reverse flux.

$$ v_{\mathrm{net}} = k_f[A] - k_r[B], \qquad \frac{d[B]}{dt}=v_{\mathrm{net}}, \qquad \frac{d[A]}{dt}=-v_{\mathrm{net}}. $$

In this simple system, setting either concentration derivative to zero sets the net flux to zero. Thus, the stationary composition is also the equilibrium composition. This coincidence depends on the setup: one reversible step, a closed system, and no other processes that can balance its flux. It should not be generalized to open or multi-step systems without checking their fluxes and boundary conditions.

Worked Example: Equilibrium in a Closed Reversible Reaction

Hypothetical data. A closed, well-mixed system has \(A \rightleftharpoons B\), with \(k_f=0.40\ \mathrm{min}^{-1}\), \(k_r=0.10\ \mathrm{min}^{-1}\), and conserved total concentration \([A]+[B]=10\ \mathrm{\mu mol\,L^{-1}}\). Assume both directions follow first-order mass-action kinetics and that no material enters or leaves.

At equilibrium, the net flux is zero, so the forward and reverse fluxes must match:

$$ k_f[A]_{\mathrm{eq}}=k_r[B]_{\mathrm{eq}}, \qquad \frac{[B]_{\mathrm{eq}}}{[A]_{\mathrm{eq}}} =\frac{k_f}{k_r} =\frac{0.40\ \mathrm{min}^{-1}}{0.10\ \mathrm{min}^{-1}} =4. $$

Let \([A]_{\mathrm{eq}}=a\). Then \([B]_{\mathrm{eq}}=4a\), and conservation gives \(a+4a=10\ \mathrm{\mu mol\,L^{-1}}\). Therefore, \([A]_{\mathrm{eq}}=2.00\ \mathrm{\mu mol\,L^{-1}}\) and \([B]_{\mathrm{eq}}=8.00\ \mathrm{\mu mol\,L^{-1}}\).

The forward flux is \((0.40\ \mathrm{min}^{-1})(2.00\ \mathrm{\mu mol\,L^{-1}})=0.800\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\). The reverse flux is \((0.10\ \mathrm{min}^{-1})(8.00\ \mathrm{\mu mol\,L^{-1}})=0.800\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\). Their difference is zero, even though each one-way flux is nonzero. In this closed one-step model, both concentrations are constant because the opposing fluxes balance.

How a Driven System Can Be at Steady State

In an open system, external conditions can maintain a persistent flow of material. A concentration may remain constant because its rate of production equals its rate of removal—not because each reaction has stopped or because the system is at thermodynamic equilibrium. One convenient mathematical statement is that the net rate of change of each steady-state species is zero, even though the fluxes contributing to that rate are nonzero.

$$ \frac{d[B]}{dt} = \text{rate of formation of }B - \text{rate of removal of }B =0. $$

The zero derivative is a balance condition. It does not require each term on the right-hand side to be zero. A fixed concentration can therefore coexist with ongoing chemical conversion, especially when boundary reservoirs or other external processes continually supply and remove material.

Worked Example: An Intermediate Sustained by Continuous Throughput

Hypothetical data. A well-mixed compartment receives \(A\) from a maintained reservoir. The formation flux of intermediate \(B\) is constant at \(3.00\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\), and \(B\) is converted to product \(C\) with first-order rate constant \(k=0.50\ \mathrm{min}^{-1}\). Product \(C\) is removed, so it does not accumulate in the compartment. Assume constant volume, and treat the reservoir conditions and formation flux as fixed.

The formation flux is \(3.00\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\); the removal flux for \(B\) is \(k[B]\). The balance equation is

$$ \frac{d[B]}{dt} = 3.00\ \mathrm{\mu mol\,L^{-1}\,min^{-1}} - (0.50\ \mathrm{min}^{-1})[B]. $$

At steady state, set the derivative to zero and solve:

$$ 0=3.00\ \mathrm{\mu mol\,L^{-1}\,min^{-1}} -(0.50\ \mathrm{min}^{-1})[B]_{\mathrm{ss}}, \qquad [B]_{\mathrm{ss}}=6.00\ \mathrm{\mu mol\,L^{-1}}. $$

At this concentration, the removal flux is \((0.50\ \mathrm{min}^{-1})(6.00\ \mathrm{\mu mol\,L^{-1}})=3.00\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\), matching the formation flux. Thus \(B\) remains constant while \(3.00\ \mathrm{\mu mol\,L^{-1}}\) of \(B\) is converted to \(C\) per minute. The maintained input and product removal drive ongoing throughput. The constant concentration of \(B\) demonstrates steady state for \(B\); it does not establish thermodynamic equilibrium.

Why Zero Concentration Change Is Not Enough

A zero derivative at one instant is weaker than steady state. Steady state describes concentrations that remain constant under the specified conditions, not merely a momentary balance that is about to change. If the forcing conditions change, a concentration can have zero derivative at the instant of the change and then move away from its previous value.

Worked Example: A Momentary Balance After the Input Changes

Hypothetical data. Let \(B\) be formed at an input rate \(u(t)\) and removed with first-order rate constant \(k=2.0\ \mathrm{min}^{-1}\). Immediately before \(t=0\), the input is \(6.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\) and \([B]=3.0\ \mathrm{\mu mol\,L^{-1}}\). At \(t=0\), the input drops to \(4.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\) and remains there. Assume the removal law and volume remain unchanged.

At the instant of the change, the derivative calculated using the new input is

$$ \left.\frac{d[B]}{dt}\right|_{t=0^+} = 4.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}} - (2.0\ \mathrm{min}^{-1})(3.0\ \mathrm{\mu mol\,L^{-1}}) =-2.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}. $$

This shows that the concentration begins to decrease immediately after the input changes; the pre-change balance does not persist. More generally, if the input had momentarily equaled \(6.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\) at \([B]=3.0\ \mathrm{\mu mol\,L^{-1}}\), the instantaneous derivative would have been zero, but a subsequent input change would break that balance. For the stated post-change conditions, the trajectory is

$$ [B](t)=2.0\ \mathrm{\mu mol\,L^{-1}} + 1.0\ \mathrm{\mu mol\,L^{-1}}\,e^{-(2.0\ \mathrm{min}^{-1})t}, \qquad t\geq 0. $$

At \(t=0\), this gives \(3.0\ \mathrm{\mu mol\,L^{-1}}\); as time passes, \([B]\) approaches \(2.0\ \mathrm{\mu mol\,L^{-1}}\), the steady-state concentration for the new input. This example emphasizes that steady-state claims require the conditions to be specified and maintained.

Comparing the Concepts and Diagnosing a Model

QuestionChemical equilibriumSteady state
What is being described?A thermodynamic condition with no sustained net reaction progress.Constant concentrations for specified species under specified conditions.
Must reaction fluxes be zero?There is no sustained net reaction flux at equilibrium; for an elementary reversible step, its forward and reverse fluxes balance.No. Formation and removal fluxes can be nonzero and equal.
Can external driving be important?A persistent external drive is incompatible with thermodynamic equilibrium.Yes. Continuous input, removal, or other maintained conditions can sustain steady state.

These descriptions are related but not interchangeable. An equilibrium state is stationary under its equilibrium conditions, so its concentrations do not change. But stationarity alone does not prove equilibrium: an open system can maintain constant concentrations through ongoing throughput. The earlier material on reversible reactions and net flux provides the flux language used here; the key additional question is whether the observed balance reflects thermodynamic equilibrium or maintained input and output.

When assessing a model or data, use the following checks:

  • Write the scope of the claim. State which species are claimed to be at steady state and whether feeds, sinks, or reservoir concentrations are held fixed.
  • Check the derivatives. Verify that the model gives \(d[X_i]/dt=0\) for each claimed steady-state species under the stated conditions. A flat plot over a short interval is evidence to inspect, not by itself proof of an exact steady state.
  • Inspect the component fluxes. If formation and removal fluxes are both positive, the system may be at a driven steady state. Do not infer zero reaction flux from a zero concentration derivative.
  • Check for changing conditions. A changing input, parameter, or boundary condition can invalidate a steady-state calculation even when the balance held before the change.
  • Use equilibrium only when justified. For equilibrium claims, check that the system is not being maintained away from equilibrium by external driving, and that the relevant reversible reactions have no sustained net progress.

A common mistake is to reason, “The concentration is constant, so the reaction must be at equilibrium.” The correction is to separate the two tests: constant concentration means the net rate of change is zero for the species being considered; equilibrium additionally requires the thermodynamic equilibrium condition. In a driven system, nonzero formation and removal fluxes can cancel exactly.

These distinctions also have limits. A numerical model may show very small derivatives because of solver tolerances or a slowly changing transient; “approximately constant over the observed interval” is not automatically an exact steady state. Likewise, a steady-state calculation for one intermediate does not imply that every species in the larger system is stationary. State which species and conditions the conclusion covers.

Key takeaway: Equilibrium concerns the absence of sustained net reaction progress; steady state concerns unchanging concentrations. In a closed, single-step reversible reaction, the two can coincide. In a continuously driven system, a steady concentration can be maintained by equal, nonzero formation and removal fluxes.

Check Your Understanding

Answer each question using the distinction between net flux and concentration change. Solutions are given below.

  1. For a closed reaction \(A \rightleftharpoons B\), the forward and reverse fluxes are each \(0.20\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\). What is the net flux, and do the nonzero one-way fluxes rule out equilibrium?
  2. An intermediate is formed at \(5.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\) and removed at \(5.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\). What can be concluded about its concentration derivative? Can equilibrium be concluded from this information alone?
  3. A species has a zero concentration derivative at one recorded instant. What additional question should be asked before calling it steady state?
  4. For a first-order removal process with \(k=0.25\ \mathrm{min}^{-1}\), a constant input of \(2.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\), and no other formation or removal processes, calculate the steady-state concentration.

Solutions. (1) The net flux is \(0.20-0.20=0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\). Nonzero one-way fluxes do not rule out equilibrium; they balance for the reversible step. (2) Its derivative is \(5.0-5.0=0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}\), so its concentration is instantaneously balanced. Equilibrium cannot be inferred: ongoing nonzero fluxes may be maintained by external driving. (3) Ask whether the derivative remains zero under maintained conditions, rather than being zero only momentarily; also identify which species and conditions are covered by the claim. (4) Set input equal to removal: \(0=2.0\ \mathrm{\mu mol\,L^{-1}\,min^{-1}}-(0.25\ \mathrm{min}^{-1})[B]_{\mathrm{ss}}\). Thus \([B]_{\mathrm{ss}}=2.0/0.25=8.0\ \mathrm{\mu mol\,L^{-1}}\).

Clinical Biostats