Calculating equilibrium concentrations of a multi-reaction-system

Tim Bernhard

Chemistry

If you have given the energies of the states as well as the transition states, or already the equilibrium constants, you can calculate the equilibrium concentrations. The way to go is to use the general formulas for the equilibrium constant \(K=exp(-E_A/(R \cdot T))\) (with \(E_A=E_{Educt} - E_{TransitionState}\)) and \(K=C_{Product}/C_{Educt}\) and apply them.

If you have a system with more than one reaction, this can become a tedious task very fast to calculate the corresponding equilibrium concentrations. That's when you want another way to go. To be precise, you need a program which does it for you. What is has to do in general is to solve a system of equations, with the same number of equations as reactions. This can generally be done by using a matrix; here, the matrix will have as many rows as reactions and as many columns as substances which are involved in the reaction stochiometry. Apart from the values in the matrix, this is equivalent to a stochiometric matrix. For the following calculations, we will start with a stochiometric matrix.

The most important step now is to determine the values of the matrix cells. You should assert on a piece of paper if the following statements apply to you too. \(m_{ij}\) is the value of the new, non-stochiometric matrix, wheras \(n_{ij}\) is the value of the sotchiometric matrix. \(K_i\) is the equilibrium constant for reaction \(i\).

\[m_{ij} = -1 \cdot n_{ij} \cdot K_i^{(n_{ij}+1)/2}\]
This only works for stochiometry 1/-1/0 as I did not find a possibility yet to square the concentration. New columns would have to be introduced but the context of e.g. [A] to [A]2 could no be given :/. Also, this post was a draft for 3 years, no idea about the context anymore. But hey, maybe it helps someone else.
Webmentions

No webmentions found for this page.