Margules activity model explained

The Margules activity model is a simple thermodynamic model for the excess Gibbs free energy of a liquid mixture introduced in 1895 by Max Margules.[1] [2] After Lewis had introduced the concept of the activity coefficient, the model could be used to derive an expression for the activity coefficients

\gammai

of a compound i in a liquid, a measure for the deviation from ideal solubility, also known as Raoult's law.

In chemical engineering the Margules Gibbs free energy model for liquid mixtures is better known as the Margules activity or activity coefficient model. Although the model is old it has the characteristic feature to describe extrema in the activity coefficient, which modern models like NRTL and Wilson cannot.

Equations

Excess Gibbs free energy

Margules expressed the intensive excess Gibbs free energy of a binary liquid mixture as a power series of the mole fractions xi:

Gex
RT

=X1X2(A21X1+A12X2)+

2
X
1
2
X
2

(B21X1+B12X2)+...+

m
X
1
m
X
2

(M21X1+M12X2)

In here the A, B are constants, which are derived from regressing experimental phase equilibria data.Frequently the B and higher order parameters are set to zero. The leading term

X1X2

assures that the excess Gibbs energy becomes zero at x1=0 and x1=1.

Activity coefficient

The activity coefficient of component i is found by differentiation of the excess Gibbs energy towards xi. This yields, when applied only to the first term and using the Gibbs–Duhem equation,:[3]

\left\{\begin{matrix}ln \gamma1=[A12+2(A21-A12

2
)x
2 \ ln \gamma

2=[A21+2(A12-A21

2
)x
1 \end{matrix}\right.

In here A12 and A21 are constants which are equal to the logarithm of the limiting activity coefficients:

ln (

infty)
\gamma
1
and
infty
ln (\gamma
2

)

respectively.

When

A12=A21=A

, which implies molecules of same molecular size but different polarity, the equations reduce to the one-parameter Margules activity model:

\left\{\begin{matrix}

2
ln \gamma
2 \ ln \gamma
2
1 \end{matrix}\right.

In that case the activity coefficients cross at x1=0.5 and the limiting activity coefficients are equal. When A=0 the model reduces to the ideal solution, i.e. the activity of a compound is equal to its concentration (mole fraction).

Extrema

Using simple algebraic manipulation, it can be stated that

dln\gamma1/dx1

increases or decreases monotonically within all

x1

range, if

A12<0

or

A21>0

with

0.5<A12/A21<2

, respectively.When

A12<A21/2

and

A12<0

, the activity coefficient curve of component 1 shows a maximum and compound 2 minimum at:

x1=

1-2A12/A21
3(1-A12/A21)

Same expression can be used when

A12<A21/2

and

A12>0

, but in this situation the activity coefficient curve of component 1 shows a minimum and compound 2 a maximum. It is easily seen that when A12=0 and A21>0 that a maximum in the activity coefficient of compound 1 exists at x1=1/3. Obviously, the activity coefficient of compound 2 goes at this concentration through a minimum as a result of the Gibbs-Duhem rule.

The binary system Chloroform(1)-Methanol(2) is an example of a system that shows a maximum in the activity coefficient of Chloroform. The parameters for a description at 20 °C are A12=0.6298 and A21=1.9522. This gives a minimum in the activity of Chloroform at x1=0.17.

In general, for the case A=A12=A21, the larger parameter A, the more the binary systems deviates from Raoult's law; i.e. ideal solubility. When A>2 the system starts to demix in two liquids at 50/50 composition; i.e. plait point is at 50 mol%. Since:

A=ln

infty
\gamma
1

=ln

infty
\gamma
2
infty
\gamma
1

=

infty
\gamma
2

>\exp(2)7.38

For asymmetric binary systems, A12≠A21, the liquid-liquid separation always occurs for

[4]

A21+A12>4

Or equivalently:

infty
\gamma
1
infty
\gamma
2

>\exp(4)54.6

The plait point is not located at 50 mol%. It depends on the ratio of the limiting activity coefficients.

Recommended values

An extensive range of recommended values for the Margules parameters can be found in the literature.[5] [6] Selected values are provided in the table below.

SystemA12A21
Acetone(1)-Chloroform(2)-0.8404-0.5610
Acetone(1)-Methanol(2)0.61840.5788
Acetone(1)-Water(2)2.04001.5461
Carbon tetrachloride(1)-Benzene (2)0.09480.0922
Chloroform(1)-Methanol(2)0.83201.7365
Ethanol(1)-Benzene(2)1.83621.4717
Ethanol(1)-Water(2)1.60220.7947

See also

Literature

  1. Max. Margules. Über die Zusammensetzung der gesättigten Dämpfe von Misschungen. Sitzungsberichte der Kaiserliche Akadamie der Wissenschaften Wien Mathematisch-Naturwissenschaftliche Klasse II. 104. 1895. 1243–1278. https://archive.org/details/sitzungsbericht10wiengoog
  2. N.A.. Gokcen. Gibbs-Duhem-Margules Laws. Journal of Phase Equilibria. 17. 1. 1996. 50–51. 10.1007/BF02648369. 95256340.
  3. Phase Equilibria in Chemical Engineering, Stanley M. Walas, (1985) p180 Butterworth Publ.
  4. Jaime. Wisniak. Liquid—liquid phase splitting—I analytical models for critical mixing and azeotropy. Chem Eng Sci. 38. 6. 1983. 969–978. 10.1016/0009-2509(83)80017-7.
  5. Book: Chemistry Data Series, Volume I: Vapor-Liquid Equilibrium Data Collection. Gmehling. J.. Onken. U.. Arlt. W.. Grenzheuser. P.. Weidlich. U.. Kolbe. B.. Rarey. J.. Dechema. 1991–2014.
  6. Book: Perry's Chemical Engineers' Handbook. Perry. Robert H.. Green. Don W.. McGraw-Hill. 1997. 978-0-07-115982-1. 7th. New York. 13:20.

External links