Stoichiometry in the mathematical modeling of chemical reactions is the manner in which the quantities of reactants and products involved in chemical reactions are tracked. Of particular interest is the number of molecules require to catalyze a reaction and the number produced by a reaction. For physical reasons, the stoichiometric quantities are typically assumed to be whole numbers.
In a general reaction of the form , a vector is assigned to the reactant and product complexes according to and . For example, in the reaction network
the following vectors are assigned
to indicate that the reactant complex of the first reaction consists of one element of the first species and second species, the product complex of the first reaction consists of one element of the third species, and so on.
The complex matrix is the matrix with columns given by the stoichiometrically distinct complexes. For the reaction network given above we have
The stoichiometric subspace is given by the span of the reaction vectors , . In other words, we define
The dimension of the stoichiometric subspace is commonly denoted . For the network above we have
In other words, there is a net loss of one molecule of species one and two and a net gain of one molecule of species three as a result of reaction one, and a loss of two molecules of species three and gain of one molecule of species four as a result of reaction two. It follows that
and the dimension is .
The stoichiometric matrix is the matrix with the column given by the reaction vector . For the network above we have
Stoichiometric compatibility classes
The stoichiometric compatibility classes of a chemical reaction network are the sets . They are affine translations of the stoichiometric subspaces and their interiors are locally homeomorphic to the Euclidean space of dimension .
Stoichiometric compatibility classes arise from the general kinetic form of a chemical reaction network. Integrating directly gives
for all . When , it follows that solutions may not wander freely about the positive orthant ; rather, they are restricted to a translation (determined by the initial condition) of stoichiometric subspace. With modest assumptions of the form of the rate functions , it follows that for all . For the above network, even though the system is five-dimensional, trajectories remain in two-dimensional translations of .
- Aizik I. Vol'pert and Sergei I. Hudjaev, Analysis in Classes of Discontinuous Functions and Equations of Mathematical Physics, Martinus Nijhoff Publishers, Dordrecht, Netherlands, 1985.