Petersen matrix

The Petersen matrix is a comprehensive description of systems of biochemical reactions used to model reactors for pollution control (engineered decomposition) as well as in environmental systems. It has as many columns as the number of relevant involved components (chemicals, pollutants, biomasses, gases) and as many rows as the number of involved processes (biochemical reactions and physical degradation). One further column is added to host the description of the kinetics of each transformation (rate equation).

Matrix structure

The mass conservation principle for each process is expressed in the rows of the matrix. If all components are included (none omitted) then the mass conservation principle states that, for each process:

${\text{for all process }}i:\sum _{j=1}^{n}a_{ij}{\dot {\rho _{j}}}=0\;,$

where ${\dot {\rho _{j}}}$ is the density rate of each component. This can also be seen as the process stoichiometric relation.

Moreover, the rate of variation of each component for all processes simultaneous effect can be easily assessed by summing the columns:

${\text{for all component }}j:{\frac {\partial C_{j}}{\partial t}}=\sum _{i=1}^{m}a_{ij}r_{i}\;,$

where $r_{i}$ are the reaction rates of each process.

Example

A system of a third order reaction followed by a Michaelis–Menten enzyme reaction.

${\ce {{A}+ 2B -> S}}$

${\ce {{E}+S<=>[k_{f}][k_{r}]ES->[k_{\mathrm {cat} }]{E}+P}}$

where the reagents A and B combine forming the substrate S (S = AB₂), which with the help of enzyme E is transformed into the product P. Production rates for each substance is:

${\begin{aligned}{\frac {d[{\ce {A}}]}{dt}}&=-k_{1}[{\ce {A}}][{\ce {B}}]^{2}\\[6pt]{\frac {d[{\ce {B}}]}{dt}}&=-2k_{1}[{\ce {A}}][{\ce {B}}]^{2}\\[6pt]{\frac {d[{\ce {S}}]}{dt}}&=k_{1}[{\ce {A}}][{\ce {B}}]^{2}-k_{f}[{\ce {E}}][{\ce {S}}]+k_{r}[{\ce {ES}}]\\[6pt]{\frac {d[{\ce {E}}]}{dt}}&=-k_{f}[{\ce {E}}][{\ce {S}}]+k_{r}[{\ce {ES}}]+k_{\ce {cat}}[{\ce {E}}S]\\[6pt]{\frac {d[{\ce {E}}S]}{dt}}&=k_{f}[{\ce {E}}][{\ce {S}}]-k_{r}[{\ce {ES}}]-k_{\ce {cat}}[{\ce {ES}}]\\[6pt]{\frac {d[{\ce {P}}]}{dt}}&=k_{\ce {cat}}[{\ce {E}}S]\end{aligned}}$

Therefore, the Petersen matrix reads as

Components (kmol/m³) Process	A	B	S	E	ES	P	Reaction rate
P1: 2nd order formation of S from A and B	−1	−2	+1	0	0	0	$k_{1}[{\ce {A}}][{\ce {B}}]^{2}$
P2: Formation of ES from E and S	0	0	−1	−1	+1	0	$k_{f}[{\ce {E}}][{\ce {S}}]$
P3: Back decomposition of ES into E and S	0	0	+1	+1	−1	0	$k_{r}[{\ce {ES}}]$
P4: Forward decomposition of ES into E and P	0	0	0	+1	−1	+1	$k_{{\ce {cat}}}[{\ce {ES}}]$

The Petersen matrix can be used to write the system's rate equation

${\begin{pmatrix}{\frac {d}{dt}}[{\ce {A}}]\\{\frac {d}{dt}}[{\ce {B}}]\\{\frac {d}{dt}}[{\ce {S}}]\\{\frac {d}{dt}}[{\ce {E}}]\\{\frac {d}{dt}}[{\ce {ES}}]\\{\frac {d}{dt}}[{\ce {P}}]\end{pmatrix}}={\begin{bmatrix}-1&0&0&0\\-2&0&0&0\\+1&-1&+1&0\\0&-1&+1&+1\\0&+1&-1&-1\\0&0&0&+1\end{bmatrix}}{\begin{pmatrix}k_{1}[{\ce {A}}][{\ce {B}}]^{2}\\k_{f}[{\ce {E}}][{\ce {S}}]\\k_{r}[{\ce {ES}}]\\k_{\ce {cat}}[{\ce {ES}}]\end{pmatrix}}$