Petersen matrix

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:

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

where ${\displaystyle {\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:

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

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

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

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

${\displaystyle {\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:

${\displaystyle {\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	${\displaystyle k_{1}[{\ce {A}}][{\ce {B}}]^{2}}$
P2: Formation of ES from E and S	0	0	−1	−1	+1	0	${\displaystyle k_{f}[{\ce {E}}][{\ce {S}}]}$
P3: Back decomposition of ES into E and S	0	0	+1	+1	−1	0	${\displaystyle k_{r}[{\ce {ES}}]}$
P4: Forward decomposition of ES into E and P	0	0	0	+1	−1	+1	${\displaystyle k_{{\ce {cat}}}[{\ce {ES}}]}$

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

${\displaystyle {\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}}}$