Modèle d'activité de Margules

Le modèle d'activité de Margules est un modèle thermodynamique simple de coefficient d'activité applicable à des espèces chimiques dans un mélange liquide. Il est présenté pour la première fois par Max Margules en 1895 comme une solution possible à la relation de Duhem-Margules.

L'enthalpie libre d'excès calculée par ce modèle présente, sous certaines conditions, des extrémums (minimums ou maximums) en fonction de la composition du mélange, ce qui permet de représenter les azéotropes des équilibres liquide-vapeur.

Description du modèle

Enthalpie libre d'excès

On considère un mélange de deux molécules en phase liquide notées $1$ et $2$.

Avec les enthalpies libres molaires ${\bar {G}}_{1}^{*}$ et ${\bar {G}}_{2}^{*}$ des corps $1$ et $2$ purs dans les mêmes conditions de température et de pression, ainsi qu'à l'état liquide comme le mélange, l'enthalpie libre molaire ${\bar {G}}$ du mélange vaut :

${\bar {G}}=x_{1}{\bar {G}}_{1}^{*}+x_{2}{\bar {G}}_{2}^{*}+{\bar {G}}^{\text{mix}}$

Dans un mélange idéal, on considère qu'il n'y a aucune interaction entre les molécules. La création d'enthalpie libre pour un mélange idéal, ou enthalpie libre molaire de mélange idéale ${\bar {G}}^{\text{mix,id}}$, est donnée par :

${{\bar {G}}^{\text{mix,id}} \over RT}=x_{1}\,\ln x_{1}+x_{2}\,\ln x_{2}$

Cette valeur est toujours négative puisque les fractions molaires sont inférieures à 1, et le mélange des deux composants est donc toujours possible.

Dans un mélange réel on doit prendre en compte les interactions entre molécules : la répulsion entre molécules $1$ et molécules $2$ a tendance à empêcher le mélange, alors que l'attraction entre ces molécules a tendance à le favoriser. Il résulte de ces interactions une enthalpie libre molaire d'excès ${\bar {G}}^{\text{E}}$ par rapport au mélange idéal. L'enthalpie libre molaire de mélange réelle ${\bar {G}}^{\text{mix}}$ vaut ainsi :

${\bar {G}}^{\text{mix}}={\bar {G}}^{\text{mix,id}}+{\bar {G}}^{\text{E}}$

Margules a exprimé l'enthalpie libre molaire d'excès en fonction des fractions molaires des composants selon :

${{\bar {G}}^{\text{E}} \over RT}=x_{1}x_{2}\left(A_{2,1}\,x_{1}+A_{1,2}\,x_{2}\right)+{x_{1}}^{2}{x_{2}}^{2}\left(B_{2,1}\,x_{1}+B_{1,2}\,x_{2}\right)+\cdots$

avec :

• $x_{1}$ et $x_{2}$ les fractions molaires des deux composants du mélange ;
• $A_{2,1}$, $A_{1,2}$, $B_{2,1}$, $B_{1,2}$ et suivants des coefficients qui décrivent l'interaction entre $1$ et $2$ dans le milieu.

Dans la grande majorité des cas on utilise l'ordre 1, ce qui permet d'obtenir l'équation à deux paramètres :

${{\bar {G}}^{\text{E}} \over RT}=x_{1}x_{2}\left(A_{2,1}\,x_{1}+A_{1,2}\,x_{2}\right)$

Cette équation peut encore être simplifiée en supposant que la fonction est symétrique, c'est-à-dire que $A_{2,1}=A_{1,2}=A$, ce qui donne :

${{\bar {G}}^{\text{E}} \over RT}=A\,x_{1}x_{2}$

puisque $x_{1}+x_{2}=1$ pour un mélange binaire.

Coefficients d'activité

Les coefficients d'activité sont obtenus en calculant la dérivée de l'enthalpie d'excès par rapport à la quantité de chaque composé, respectivement $n_{1}$ et $n_{2}$. Si on utilise l'équation d'ordre 1 à deux coefficients :

${{\bar {G}}^{\text{E}} \over RT}=x_{1}x_{2}\left(A_{2,1}\,x_{1}+A_{1,2}\,x_{2}\right)$

on obtient :

$\ln \gamma _{1}=\left({\partial \over \partial n_{1}}\left[{\left(n_{1}+n_{2}\right)\,{\bar {G}}^{\text{E}} \over RT}\right]\right)_{n_{2}}=\left(A_{1,2}+2\left(A_{2,1}-A_{1,2}\right)x_{1}\right){x_{2}}^{2}$

$\ln \gamma _{2}=\left({\partial \over \partial n_{2}}\left[{\left(n_{1}+n_{2}\right)\,{\bar {G}}^{\text{E}} \over RT}\right]\right)_{n_{1}}=\left(A_{2,1}+2\left(A_{1,2}-A_{2,1}\right)x_{2}\right){x_{1}}^{2}$

Les coefficients d'activité à dilution infinie sont obtenus selon :

$\ln \gamma _{1,2}^{\infty }=\lim _{x_{1}\to 0 \atop x_{2}\to 1}\ln \gamma _{1,2}=A_{1,2}$

$\ln \gamma _{2,1}^{\infty }=\lim _{x_{2}\to 0 \atop x_{1}\to 1}\ln \gamma _{2,1}=A_{2,1}$

Si on utilise l'équation simplifiée à 1 paramètre, l'expression des coefficients d'activité se simplifie en :

$\ln \gamma _{1}=A\,{x_{2}}^{2}$

$\ln \gamma _{2}=A\,{x_{1}}^{2}$

$\ln \gamma _{1,2}^{\infty }=A$

$\ln \gamma _{2,1}^{\infty }=A$

Extrémum

En dérivant $\ln \gamma _{1}$ par rapport à la fraction molaire d'un composant du mélange, par exemple $x_{1}$, on peut calculer la position d'un extrémum du coefficient d'activité à la position $x_{1,{\text{max}}}$ :

$x_{1,{\text{max}}}={1-2\,A_{1,2}/A_{2,1} \over 3\left(1-A_{1,2}/A_{2,1}\right)}$

Pour que cet extrémum soit compris entre 0 et 1, il faut que $A_{1,2}/A_{2,1}<1/2$ ou $A_{1,2}/A_{2,1}>2$.

Quelques exemples

Le tableau suivant regroupe quelques valeurs pour les coefficients de Margules trouvées dans la littérature scientifique et technique.

corps 1	corps 2	$A_{1,2}$	$A_{2,1}$
acétone	chloroforme	-0,8404	-0,5610
acétone	méthanol	0,6184	0,5788
acétone	eau	2,0400	1,5461
tétrachlorure de carbone	benzène	0,0948	0,0922
chloroforme	méthanol	0,8320	1,7365
éthanol	benzène	1,8362	1,4717

Voir aussi

Notes et références

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Margules activity model » (voir la liste des auteurs).

Liens internes

• Activité chimique
• Relation de Duhem-Margules