From a Reaction Step to a Rate
A QSP model represents biological mechanisms as changes among molecular species. To describe those changes quantitatively, it needs more than a list of species: it needs to specify which species participate in each reaction event and how frequently those events occur. A rate law connects the concentrations of participating reactants to the event rate. That connection is particularly direct for an elementary reaction step, because the step represents one mechanistic event rather than a summary of several hidden events.
This tutorial focuses on translating an elementary step into its mass-action rate law. The key restriction is that the step must genuinely be elementary. A reaction equation that summarizes an overall process does not, by itself, establish the rate law for that process.
For example, \(A+B\rightarrow C\) says that an event consumes one molecule of \(A\) and one molecule of \(B\), and forms one molecule of \(C\). It does not mean that one molecule of \(A\) is somehow converted into one molecule of \(B\); the arrow separates the reactants from the products. The coefficients specify the amounts participating in each event. If no coefficient is written, it is understood to be one.
In a concentration-based model, \([A]\) denotes the concentration of species \(A\), commonly in mol L\(^{-1}\), abbreviated M. A rate law predicts the event rate from the reactant concentrations and a rate constant. It describes how fast events occur under specified conditions; it is not itself a count of how much material has accumulated.
The Mass-Action Rule for an Elementary Step
For an elementary step, the mass-action rule assigns one concentration factor for each reactant molecule participating in an event. In general, for the step \(aA+bB\rightarrow\) products, the rate law is proportional to \([A]^a[B]^b\). The proportionality constant \(k\) is the rate constant for that step.
The exponents in this rule come from the molecular participation in the elementary event. In a single event involving two \(A\) molecules, for instance, the rate depends on the concentration of \(A\) twice, so the concentration factor is \([A]^2\). The rate constant collects the proportionality between this concentration product and the event rate; its units therefore depend on the concentration factors in the law.
A brief probabilistic justification helps explain the product form. In a well-mixed system, the chance of finding a specified reactant molecule in a small volume is proportional to its concentration. For an event requiring one \(A\) and one \(B\), the chance of a suitable encounter is proportional to the chance of finding \(A\) times the chance of finding \(B\), giving a factor proportional to \([A][B]\). If the event requires two \(A\) molecules, the chance of finding a suitable pair scales with \([A]^2\). For a finite population of identical molecules, the number of distinct pairs is proportional to \(N(N-1)/2\); when the population is large, this scales approximately with \(N^2\). Fixed volume converts these molecule-number dependencies into concentration dependencies, while constant factors are included in \(k\). This justification relies on the step being elementary and on mixing and conditions being sufficiently uniform for concentrations to represent encounter opportunities.
The rule does not license reading a rate law from every reaction equation. A net equation can conceal a sequence of events, and the concentrations of intermediate species or the slowest step can affect the observed rate. The mass-action exponents follow directly from the reactant participation only when the equation describes an elementary step.
Worked Example: One Reactant in an Elementary Step
Suppose a hypothetical elementary step is \(A\rightarrow P\), with rate constant \(k=0.12\ \text{s}^{-1}\). At a particular time, \([A]=0.050\ \text{M}\). Find the event rate.
The step has one \(A\) molecule as a reactant, so the mass-action rate law is
Substitute the stated values:
Thus, under the stated conditions, the hypothetical event rate is \(0.0060\ \text{mol L}^{-1}\text{s}^{-1}\). The units are concentration per time because \(\text{s}^{-1}\) multiplied by concentration gives concentration per time.
Worked Example: Two Different Reactants
A hypothetical elementary step is \(A+B\rightarrow C\). Its rate constant is \(k=0.40\ \text{L mol}^{-1}\text{s}^{-1}\). At a specified time, \([A]=0.030\ \text{M}\) and \([B]=0.020\ \text{M}\). Calculate the event rate.
There is one molecule of each reactant in an event, so the rate law is \(r=k[A][B]\). Substitution gives
The rate is \(2.4\times 10^{-4}\ \text{mol L}^{-1}\text{s}^{-1}\). The units confirm the result: \((\text{L mol}^{-1}\text{s}^{-1})(\text{mol L}^{-1})^2=\text{mol L}^{-1}\text{s}^{-1}\). If either reactant concentration were zero, this law would give zero event rate because the specified event requires both reactants.
Units, Conditions, and Scope
A useful check is to make the units of the rate law agree with the units of the event rate. If concentration is measured in M and time in seconds, then \(r\) must have units M s\(^{-1}\). The units of \(k\) must supply whatever is needed to turn the concentration factors into M s\(^{-1}\). They are not arbitrary labels: they are part of a dimensionally consistent model.
The rate law also has a domain of validity. It assumes that the stated step is elementary, that concentrations adequately describe the relevant reactant availability, and that the rate constant applies under the conditions being modeled. Temperature, solvent, pH, or other conditions can alter a rate constant; treating \(k\) as fixed means those conditions are being held fixed or their variation is being neglected. In a spatially nonuniform system, a single well-mixed concentration may not represent local encounters. In that case, a spatially resolved description or a different effective model may be needed.
Worked Example: Two Identical Reactants and a Unit Check
Consider the hypothetical elementary step \(2A\rightarrow D\), with \(k=0.50\ \text{L mol}^{-1}\text{s}^{-1}\) and \([A]=0.040\ \text{M}\). Calculate the event rate and check its units.
Two \(A\) molecules participate in each event, so the mass-action law is \(r=k[A]^2\). The calculation is
The predicted event rate is \(8.0\times 10^{-4}\ \text{mol L}^{-1}\text{s}^{-1}\). For the units, \((\text{L mol}^{-1}\text{s}^{-1})(\text{mol L}^{-1})^2=\text{mol L}^{-1}\text{s}^{-1}\), as required. This calculation concerns the rate of events, not a separately written equation for each species concentration.
How to Check a Proposed Rate Law
Before using a mass-action law in a QSP model, check that its meaning and assumptions match the mechanism. The following sequence is useful in model construction and review.
Write the elementary step and identify which species are consumed and formed in one event.
For an elementary step, use each reactant's participation amount as the power on its concentration.
Multiply the rate constant by the corresponding reactant concentration factors.
Confirm that the result has concentration-per-time units and that the assumptions behind the elementary-step law are defensible.
Several observations can signal that a proposed law or its use is going wrong. A rate with units of concentration alone is not an event rate. A law that predicts a nonzero rate when a required reactant is absent contradicts the specified event. A negative concentration supplied to a positive-integer power is outside the physical domain of this concentration model. Finally, a rate law that fits measurements poorly may indicate incorrect parameters, missing mechanisms, changing conditions, or an incorrect assumption that the proposed step is elementary. A unit check can catch algebraic errors, but it cannot establish that the mechanism itself is correct.
Mass action is not the only way to represent reaction rates. An empirical rate law can be fitted to observed rates without asserting that its powers arise from a single elementary event. Such a law may be appropriate when the microscopic mechanism is unresolved, but it requires evidence over the concentration and condition range where it will be used. For an established elementary step, mass action provides a mechanistic link between event participation and concentration dependence. For a net process with hidden steps, imposing that link directly may be unjustified.
Common Mistake and Key Limitation
A frequent mistake is to take an overall reaction equation and assign its coefficients as rate-law exponents without checking whether the equation represents one elementary event. For example, the equation \(A+B\rightarrow C\) might describe a single event, or it might summarize several intermediate steps. Only in the first case does the elementary-step mass-action rule directly give \(r=k[A][B]\). In the second case, the observed rate depends on the underlying mechanism, so the overall equation alone is not enough to determine the law.
Another error is to omit a participating reactant or to use a concentration unit inconsistent with the fitted value of \(k\). The correction is to start from the stated elementary event, write one concentration factor for each participating reactant molecule, and check dimensions explicitly. Even then, the result is a model under assumptions, not proof that the proposed mechanism is biologically complete.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions below.
- For the hypothetical elementary step \(X+Y\rightarrow Z\), write the mass-action rate law.
- A proposed rate law is \(r=k[A][B]\), where \([A]\) and \([B]\) are in M and \(r\) is in M s\(^{-1}\). What units must \(k\) have?
- For \(2X\rightarrow P\), with \(k=0.10\ \text{L mol}^{-1}\text{s}^{-1}\) and \([X]=0.030\ \text{M}\), calculate the event rate.
- Why is it not generally valid to infer rate-law exponents from the coefficients in an overall reaction equation?
Solutions. (1) The step has one \(X\) and one \(Y\) per event, so \(r=k[X][Y]\). (2) The product \([A][B]\) has units M\(^{2}\); therefore \(k\) must have units M\(^{-1}\)s\(^{-1}\), equivalently L mol\(^{-1}\)s\(^{-1}\), so that \(r\) has units M s\(^{-1}\). (3) The elementary step gives \(r=k[X]^2=(0.10)(0.030)^2=0.000090\ \text{mol L}^{-1}\text{s}^{-1}\), or \(9.0\times10^{-5}\ \text{mol L}^{-1}\text{s}^{-1}\). (4) An overall equation can summarize multiple unresolved events. Its coefficients describe net stoichiometry, but only the reactant participation in a single elementary event directly supplies the exponents in the mass-action law for that step.