Excess Gibbs Energy models

[siddharth]

These are some empirical excess Gibbs energy models used to calculate activity coefficients for non-ideal liquid mixtures. I didn't find many on the web, so I thought I'll put them here for future reference

All x_i refer to the mole fraction of the ith species in the solution.
To calculate the activity coefficient \gamma_i from the excess free energy, we have
ln \gamma_i=\left(\frac{\partial\left(nG^{E}/RT\right)}{\partial n_i} \right)_{P,T,n_j}

Most of these are for binary solutions, but some can be extended into multi-component systems.

1) Van Laar equation

\frac{G^{E}}{RT} = \frac{ABx_1x_2}{Ax_1+B_x_2}

So that,

ln \gamma_1 = A\left(1+\frac{Ax_1}{Bx_2}\right)^{-2}

ln \gamma_2 = B\left(1+\frac{Bx_2}{Ax_1}\right)^{-2}

2) Porters equation

\frac{G^{E}}{RT}=Ax_1x_2

So that,

ln \gamma_1 = Ax_2^2

ln \gamma_2= Ax_1^2

3) Margules equation

\frac{G^{E}}{RT}=x_1x_2\left(Ax_1+Bx_2\right)

So that,

ln \gamma_1 = x_2^2\left(A+2\left(B-A\right)x_1\right)

ln \gamma_2 = x_1^2\left(B+2\left(A-B\right)x_2\right)

4) Wilsons equation (*local composition model*)

\frac{G^{E}}{RT}=-x_1ln \left(x_1+A'x_2\right) - x_2ln \left(x_2 + B'x_1\right)

So that,

ln \gamma_1=-ln(x1+A'x_2) + x2 \left(\frac{A'}{x_1+A'x_2} - \frac{B'}{x_2+B'x_1}\right)

ln \gamma_2 = -ln(x2+B'x_1) - x1 \left(\frac{A'}{x_1+A'x_2} - \frac{B'}{x_2+B'x_1}\right)

I'll add the NRTL equation and maybe the UNIQUAC model later. Feel free to correct me if you spot any errors.