Stoichiometric Matrices

Introduction

A network diagram becomes easier to calculate with when its reaction coefficients are placed in a matrix. Each row represents a species, each column a reaction and each entry records how many units of that species are created or consumed by one occurrence of the reaction. This stoichiometric matrix separates chemical bookkeeping from kinetic assumptions and provides a compact route to mass balances.

Core explanation

Choose a fixed species order and reaction order. For A → B, the column is (−1, +1)ᵀ for rows A and B. For 2B → C, entries are 0 for A, −2 for B and +1 for C. A negative entry denotes consumption, a positive entry production, and zero means the species does not participate in that step. The matrix depends on how reactions are written; reversing a reaction reverses its column sign. Its entries have stoichiometric meaning, not kinetic-order meaning.

Let N be the stoichiometric matrix, v the vector of reaction fluxes in concentration per time, and c the species concentration vector. In a well-mixed closed, constant-volume system, dc/dt = Nv. Row i adds contributions from every reaction to the net rate of species i . If reaction extents ξ are measured in concentration units, a finite composition change is Δc = NΔξ. Open reactors add inlet and outlet terms. An ACS chemical-network treatment uses this matrix–flux structure to connect network topology with dynamics.

The same N can be paired with different kinetic laws. The step A → B could have v = k[A] if elementary and first order, or a more complex expression if the written arrow lumps multiple steps. Matrix entries do not reveal which law is correct. They do allow immediate consistency checks: an element-composition row vector times N should be zero for each closed-system reaction. Charge conservation can be checked similarly if electrons and counterions are represented properly.

Matrix construction requires careful species definitions. H⁺, H₂O and OH⁻ may need explicit rows for acid–base chemistry. An adsorbed species A and vacant site are distinct rows if a catalyst site balance is important. If a reaction is reversible, one can represent it as one net flux that may be positive or negative, or as separate forward and reverse columns with nonnegative fluxes. Both conventions are valid if applied consistently.

Step-by-step reasoning

1. Write each reaction balanced with all tracked species. 2. Fix row order for species and column order for reactions. 3. Enter products as positive and reactants as negative coefficients. 4. Multiply each column by its flux to obtain its contribution to species rates. 5. Add columns and verify conservation laws before fitting kinetics.

Visual explanation

Draw a table with rows A, B, C and columns R₁: A → B and R₂: 2B → C. Fill R₁ with −1, +1, 0; fill R₂ with 0, −2, +1. Put a flux vector (v₁, v₂) beside it and show the resulting row equations: d[A]/dt = −v₁, d[B]/dt = v₁ − 2v₂ and d[C]/dt = v₂.

Real-world analogy

A warehouse ledger lists each transaction as a column: remove one part A, add one part B, or remove two B and add one C. The transaction matrix shows inventory changes regardless of how often transactions occur. Chemical stoichiometry works similarly, but conservation of atoms imposes extra constraints on which ledger transactions are physically possible.

Real-world example

A kinetic modeller studies A → B → C and A → D. Writing N immediately reveals that A disappears through two columns and B is made and consumed. When measurements show C rising, the model can attribute it to the B → C flux. If a proposed side reaction produces D without accounting for an oxygen atom, an element-matrix check detects the missing species before parameter fitting begins.

Why?

Why separate N from v? Stoichiometry is fixed by the chosen reactions, while flux changes with concentrations, temperature and catalysts. The same balanced network can be tested with different kinetic hypotheses without rebuilding the conservation bookkeeping each time.

Common misconception

“A coefficient of 2 in N proves second-order kinetics” is false for an overall or lumped reaction. “Every nonzero entry is positive” ignores consumption. “Matrix order changes chemistry” is false if row and column labels are tracked, though misordered multiplication can create errors. “A valid matrix proves the mechanism” is false; it only encodes the proposed reactions consistently.

Worked example

For rows A, B, C and columns R₁: A → B, R₂: 2B → C, use N = [[−1,0],[+1,−2],[0,+1]]. If v₁ = 4 mmol L⁻¹ min⁻¹ and v₂ = 1 mmol L⁻¹ min⁻¹, then Nv = (−4, 4 − 2, 1) = (−4, +2, +1) mmol L⁻¹ min⁻¹. A falls by 4 concentration units per minute, B rises by 2 and C rises by 1. If B contains the same atoms as A and C represents a dimer of B, the conserved count is [A] + [B] + 2[C]; its derivative is −4 + 2 + 2(1) = 0. A naive total [A] + [B] + [C] would not be conserved because C contains two units.

Quick check

1. What does a negative stoichiometric matrix entry mean? Answer: That species is consumed by the reaction represented by that column.

Exam focus

Construct N for a two- or three-reaction network with labelled row and column order. Multiply by a flux vector to calculate net species rates. Distinguish stoichiometric coefficients from kinetic orders and verify a conserved atom or moiety count.

Advanced insight

Linear algebra reveals more than rates. A left-null vector l satisfying lᵀN = 0 defines a conserved quantity in a closed system. A right-null vector Nv = 0 describes a set of reaction fluxes that creates no net change in tracked species. Interpretation depends on the network boundary and whether all relevant species are included.

Summary

A stoichiometric matrix encodes which species each reaction consumes or produces. Multiplying it by reaction fluxes gives species rates, while its null spaces expose conservation and possible steady flux patterns. Kinetic laws remain separate hypotheses.

Practice questions

1. What are N entries for A → 2B in rows A, B? Answer: −1 for A and +2 for B. 2. If v₁ = 3 for A → B and v₂ = 1 for B → C, what is d[B]/dt? Answer: v₁ − v₂ = 2 in the flux units. 3. Does N alone tell whether a reaction is first order? Answer: No. It contains stoichiometry, not an experimentally established rate law. 4. What does lᵀN = 0 represent? Answer: A weighted species combination conserved by every encoded reaction in a closed system.