Neha Patil (Editor)

Petersen matrix

Updated on
Edit
Like
Comment
Share on FacebookTweet on TwitterShare on LinkedInShare on Reddit

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 got 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).

Contents

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:

for all process  i : j = 1 n a i j ρ j ˙ = 0 ,

where ρ 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:

for all component  j : C j t = i = 1 m a i j 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.

A + 2 B S E + S k r k f ES k c a t E + P

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

d [ A ] d t = k 1 [ A ] [ B ] 2 d [ B ] d t = 2 k 1 [ A ] [ B ] 2 d [ S ] d t = k 1 [ A ] [ B ] 2 k f [ E ] [ S ] + k r [ ES ] d [ E ] d t = k f [ E ] [ S ] + k r [ ES ] + k cat [ E S ] d [ E S ] d t = k f [ E ] [ S ] k r [ ES ] k cat [ ES ] d [ P ] d t = k cat [ E S ]

Therefore, the Petersen matrix reads as

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

( d d t [ A ] d d t [ B ] d d t [ S ] d d t [ E ] d d t [ ES ] d d t [ P ] ) = [ 1 0 0 0 2 0 0 0 + 1 1 + 1 0 0 1 + 1 + 1 0 + 1 1 1 0 0 0 + 1 ] ( k 1 [ A ] [ B ] 2 k f [ E ] [ S ] k r [ ES ] k cat [ ES ] )

References

Petersen matrix Wikipedia