Calculating equilibrium concentrations of a multi-reaction-system
ChemistryThis post is old. Codes, opinions and facts could be outdated.
Help to fix potential errors: GitHub
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}\]