Voir la source de Activité chimique

{{Voir homonymes|Activité}}

En [[chimie physique]], et plus particulièrement en [[thermodynamique]], l''''activité chimique''', ou '''activité''', d'une [[espèce chimique]] exprime l'écart entre les propriétés de cette espèce pure ou dans un mélange réel et les propriétés de cette même espèce dans un [[état standard]] à la même température. La notion d'activité chimique est surtout employée pour les [[Phase (thermodynamique)|phases]] liquide et solide. Elle permet notamment le calcul des [[diagramme de phase|équilibres de phases]] et des [[équilibre chimique|équilibres chimiques]].

De nombreux modèles de '''coefficients d'activité''' ont été développés pour les phases liquides, avec comme état standard les liquides purs. Les '''grandeurs d'excès''', calculées à partir des coefficients d'activité, expriment l'écart entre les [[Extensivité et intensivité (physique)#extensivité|grandeurs extensives]] d'un mélange réel et celles d'une [[solution idéale]].

== Historique ==

À la suite de l'apparition de l'[[équation d'état de van der Waals]] en 1873, il devenait possible de calculer les équilibres de phases gaz-liquide. En effet, cette équation novatrice pour son époque permettait de calculer aussi bien les propriétés d'un gaz que celles d'un liquide. [[Willard Gibbs]] et [[Pierre Duhem]] introduisirent peu après (1875 à 1878) le [[potentiel chimique]], important dans les équilibres de phases et les [[équilibre chimique|équilibres chimiques]]. Toutefois, le potentiel chimique se révéla difficile à manipuler, car il ne peut être calculé qu'à une constante additive près et non de façon absolue, de plus le potentiel chimique de toute espèce tend vers l'infini négatif à dilution infinie. En 1900<ref>{{article |langue=en |volume=36 |numéro=9 |pages=145-168 |auteur=[[Gilbert Lewis|Gilbert Newton Lewis]] |titre=A New Conception of Thermal Pressure and a Theory of Solutions |journal=Proceedings of the American Academy of Arts and Sciences |date=octobre 1900 |doi=10.2307/20020988}}.</ref> et 1901<ref name=Lewis1901>{{article |langue=en |volume=37 |numéro=3 |pages=49–69 |auteur=[[Gilbert Lewis|Gilbert Newton Lewis]] |titre=The Law of Physico-Chemical Change |journal=Proceedings of the American Academy of Arts and Sciences |date=juin 1901 |doi=10.2307/20021635}}.</ref>, [[Gilbert Lewis]] introduisit la [[fugacité]], qui décrit l'écart de comportement d'un corps réel, pur ou en mélange, par rapport au même corps à l'état de gaz parfait pur. Cette notion se révéla efficace dans son application pour les gaz, mais les équations d'état représentant assez mal les phases liquides, Lewis introduisit en 1923<ref name="Lewis1923">{{Ouvrage |langue=en |auteur1=[[Gilbert Lewis|Gilbert Newton Lewis]] |auteur2=Merle Randall |titre=Thermodynamics and the free energy of chemical substances |sous-titre= |éditeur=McGraw-Hill Book Company Inc. |collection= |lieu= |année=1923 |numéro d'édition= |volume= |tome= |pages totales= |passage= |isbn= |lire en ligne= |id= }}.</ref> l'activité chimique, plus spécialement appliquée aux phases condensées (liquide ou solide). L'activité chimique est surtout utilisée dans l'expression des [[Cinétique chimique|vitesses de réaction]] et des [[constante d'équilibre|constantes d'équilibre]] dans l'étude des [[réaction chimique|réactions]] et équilibres chimiques, et dans le calcul des [[coefficient de partage|coefficients de partage]] dans l'étude des [[diagramme de phase|équilibres de phases]].

== Définitions ==
=== Activité ===

La [[fugacité]] <math>f_i</math> d'un corps <math>i</math> a été définie par Lewis par la variation isotherme du [[potentiel chimique]] <math>\mu_i</math> de ce corps :

:<math>\mathrm{d} \mu_i = RT \, \mathrm{d} \ln f_i</math>

En intégrant cette relation à température constante, on obtient la relation entre les potentiels chimiques <math>\mu_i^1</math> et <math>\mu_i^2</math> d'un corps <math>i</math> dans deux états différents 1 et 2 :

:<math>\mu_i^2 - \mu_i^1 = RT \int_{\ln f_i^1}^{\ln f_i^2 } \, \mathrm{d} \ln f_i = RT \ln \! \left( \frac{f_i^2}{f_i^1} \right)</math>

<math>f_i^1</math> et <math>f_i^2</math> étant les fugacités dans les états respectifs. Ces deux états peuvent être :
* à des pressions différentes : <math>P^1</math> et <math>P^2</math> ;
* de compositions différentes : <math>x^1=\left[x_1^1, x_2^1, \cdots, x_N^1\right]</math> et <math>x^2=\left[x_1^2, x_2^2, \cdots, x_N^2\right]</math>, avec <math>x_i^k</math> la [[fraction molaire]] de l'espèce <math>i</math> dans l'état <math>k</math> ; par ex. :
** l'état 1 peut être le corps <math>i</math> pur (<math>x_1^1 = 0, x_2^1 = 0, \cdots, x_i^1 = 1, \cdots, x_N^1 = 0</math>) ;
** l'état 2 un mélange contenant le corps <math>i</math> (<math>x_1^2 \neq 0, x_2^2 \neq 0, \cdots, x_i^2 \neq 0, \cdots, x_N^2 \neq 0</math>) ;
* dans des phases différentes : par ex. gaz parfait pour l'état 1 et liquide pour l'état 2.
Néanmoins, '''les deux états sont à la même température <math>T</math>''' :

:<math>\mu_i^2 \! \left( P^2,T,x^2 \right) - \mu_i^1 \! \left( P^1,T,x^1 \right) 
= RT \ln \! \left( \frac{f_i^2 \! \left( P^2,T,x^2 \right)}{f_i^1 \! \left( P^1,T,x^1 \right)} \right)</math>
<br>
{{boîte déroulante/début|titre = Note sur la fugacité et l'état standard gaz parfait pur.}}

La fugacité est définie par rapport à l'état standard gaz parfait pur. La relation ci-dessus est tout à fait compatible avec cette définition. En effet, on a les relations suivantes entre potentiels chimiques du même corps <math>i</math> dans divers états, ''à la même température'' :

* <math>\mu_i^1 \! \left( P^1,T,x^1 \right) - \mu_i^{\bullet,*} \! \left( P^1,T \right)
= RT \ln \! \left( \frac{f_i^1 \! \left( P^1,T,x^1 \right)}{P^1} \right)</math> pour le passage isotherme du gaz parfait pur dans les conditions de pression <math>P^1</math> et température <math>T</math> à l'état 1 dans les mêmes conditions ;
* <math>\mu_i^{\bullet,*} \! \left( P^2,T \right) - \mu_i^{\bullet,*} \! \left( P^1,T \right) 
= RT \ln \! \left( \frac{P^2}{P^1} \right)</math> pour le changement de pression isotherme d'un gaz parfait pur ;
* <math>\mu_i^2 \! \left( P^2,T,x^2 \right) - \mu_i^{\bullet,*} \! \left( P^2,T \right) 
= RT \ln \! \left( \frac{f_i^2 \! \left( P^2,T,x^2 \right)}{P^2} \right)</math> pour le passage isotherme du gaz parfait pur dans les conditions de pression <math>P^2</math> et température <math>T</math> à l'état 2 dans les mêmes conditions.

On a donc la relation :

:<math>\mu_i^2 \! \left( P^2,T,x^2 \right) - \mu_i^1 \! \left( P^1,T,x^1 \right)</math>
:<math>= \left[ \mu_i^2 \! \left( P^2,T,x^2 \right) - \mu_i^{\bullet,*} \! \left( P^2,T \right) \right] + \left[ \mu_i^{\bullet,*} \! \left( P^2,T \right) - \mu_i^{\bullet,*} \! \left( P^1,T \right) \right] + \left[ \mu_i^{\bullet,*} \! \left( P^1,T \right) - \mu_i^1 \! \left( P^1,T,x^1 \right) \right]</math>
:<math>= RT \ln \! \left( \frac{f_i^2 \! \left( P^2,T,x^2 \right)}{P^2} \frac{P^2}{P^1} \frac{P^1}{f_i^1 \! \left( P^1,T,x^1 \right)} \right)</math>
:<math>= RT \ln \! \left( \frac{f_i^2 \! \left( P^2,T,x^2 \right)}{f_i^1 \! \left( P^1,T,x^1 \right)} \right)</math>

{{boîte déroulante/fin}}

Dans la pratique, l'état 2 est l'état réel du mélange et l'état 1 un '''[[état standard]]''' dans lequel il est facile de calculer les propriétés du corps <math>i</math> ''à la même température'' que l'état réel. On note alors :

:<math>\mu_i - \mu_i^\circ = RT \ln \! \left( \frac{f_i}{f_i^\circ} \right) = RT \ln a_i</math>

avec :
* <math>T</math> la [[température]] du mélange réel ;
* <math>\mu_i</math> le [[potentiel chimique]] du corps <math>i</math> dans le mélange réel ;
* <math>\mu_i^\circ</math> le potentiel chimique du corps <math>i</math> dans l'[[état standard]] à la même <math>T</math> que le mélange réel ;
* <math>f_i</math> la [[fugacité]] du corps <math>i</math> dans le mélange réel ;
* <math>f_i^\circ</math> la fugacité du corps <math>i</math> dans l'état standard à la même <math>T</math> que le mélange réel ;
* <math>a_i</math> l'activité du corps <math>i</math> dans le mélange réel par rapport à l'état standard à la même <math>T</math> que le mélange réel ;
* <math>R</math> la [[constante universelle des gaz parfaits]].

L''''activité chimique''' est définie par le rapport des fugacités réelle et dans l'état standard ''à la même température'' que l'état réel :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Activité chimique : '''<math>a_i = \frac{f_i}{f_i^\circ}</math>
|}

soit également :

<center>'''Activité chimique : '''<math>\mu_i - \mu_i^\circ = RT \ln a_i</math></center>

L'activité est [[Grandeur sans dimension|adimensionnelle]].

=== Coefficient d'activité ===

Lorsque, pour le calcul des propriétés d'un mélange réel, l'[[état standard]] de chaque corps <math>i</math> est pris comme étant le corps <math>i</math> pur aux mêmes pression <math>P</math>, température <math>T</math> et dans la même phase (gaz, liquide ou solide) que le mélange réel, les potentiels chimiques réel <math>\mu_{i}</math> (fugacité <math>f_i</math>) et du corps pur <math>\mu_{i}^*</math> (fugacité <math>f_i^*</math>) sont liés par la relation :

:<math>\mu_{i} - \mu_{i}^* = RT \ln \! \left( \frac{f_i}{f_i^*} \right) = RT \ln a_i</math>

Par définition, le potentiel chimique <math>\mu_{i}^\text{id}</math> du corps <math>i</math> dans une [[solution idéale]] et le potentiel chimique <math>\mu_{i}^*</math> du corps <math>i</math> pur, tous deux définis aux mêmes pression <math>P</math> et température <math>T</math>, et pour la même phase, sont liés par la relation :

:<math>\mu_{i}^\text{id} = \mu_{i}^* + RT \ln x_i</math>

avec <math>x_i</math> la [[fraction molaire]] du corps <math>i</math> dans le mélange idéal.

On peut alors écrire, en considérant une solution idéale ''aux mêmes pression, température et composition, et dans la même phase'' (gaz, liquide ou solide) que le mélange réel :

:<math>\mu_{i} - \mu_{i}^\text{id} = RT \ln \! \left( \frac{f_i}{f_i^\text{id}} \right) = RT \ln \! \left( \frac{a_i}{x_i} \right)</math>

avec <math>f_{i}^\text{id} = x_i f_i^*</math> la fugacité en solution idéale. Le rapport entre la fugacité réelle et la fugacité en solution idéale est appelé '''coefficient d'activité''', il est noté <math>\gamma_i</math> :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Coefficient d'activité : '''<math>\gamma_i = \frac{f_i}{f_i^\text{id}} = \frac{f_i}{x_i f_i^*} = \frac{a_i}{x_i}</math>
|}

Le coefficient d'activité est [[Grandeur sans dimension|adimensionnel]].

Le coefficient d'activité exprime l'écart entre le potentiel chimique <math>\mu_i</math> d'un corps dans un mélange réel et le potentiel chimique <math>\mu_{i}^\text{id}</math> de ce corps dans la solution idéale aux mêmes pression <math>P</math>, température <math>T</math> et composition (fractions molaires <math>x_i</math>), dans la même phase (gaz, liquide ou solide) :

<center>'''Coefficient d'activité : '''<math>\mu_{i} - \mu_{i}^\text{id} = RT \ln \gamma_i</math></center>

Le coefficient d'activité exprime également l'écart entre le potentiel chimique <math>\mu_i</math> d'un corps dans un mélange réel et le potentiel chimique <math>\mu^*_i</math> de ce [[corps pur]] dans la même phase, aux mêmes pression et température :

:<math>\mu_i - \mu^*_i = RT \ln \! \left( x_i \gamma_i \right)</math>

Par définition, le coefficient d'activité d'un corps dans une solution idéale est donc égal à 1 :

<center>'''Dans une solution idéale : '''<math>\gamma_i^\text{id} = 1</math> et <math>a_i^\text{id} = x_i</math></center>

À fortiori pour un corps pur :

<center>'''Pour un corps pur : '''<math>\gamma_i^* = 1</math> et <math>a_i^* = 1</math></center>

La solution idéale correspond au cas idéal dans lequel les molécules du mélange s'attirent ou se repoussent entre elles, toutes espèces confondues, de la même façon que les molécules de chacun des constituants considéré pur s'attirent ou se repoussent entre elles.

Lorsque les molécules d'un mélange réel s'attirent plus entre elles que dans la solution idéale, les activités de tous les composants sont inférieures aux fractions molaires : dans ce cas <math>\gamma_i < 1</math> et <math>a_i < x_i</math>.

Inversement, lorsque les molécules d'un mélange réel se repoussent plus entre elles que dans la solution idéale, les activités de tous les composants sont supérieures aux fractions molaires : dans ce cas <math>\gamma_i > 1</math> et <math>a_i > x_i</math>.

;Note sur la relation entre coefficient de fugacité et coefficient d'activité

Avec les [[Fugacité|coefficients de fugacité]] respectifs du corps <math>i</math> dans le mélange réel <math>f_{i} = x_i \phi_i P</math> et du corps <math>i</math> pur <math>f_{i}^* = \phi_i^* P</math>, nous avons la relation entre coefficient d'activité et coefficients de fugacité :

<center>'''Coefficient d'activité : '''<math>\gamma_i = \frac{\phi_i}{\phi_i^*}</math></center>

== Calcul de l'activité et du coefficient d'activité ==
=== Expressions usuelles de l'activité ===

Pour le calcul d'un [[équilibre chimique]], certaines approximations sont souvent faites dans l'expression de l'activité chimique.

'''Pour un mélange de [[gaz parfait]]s''', l'état de référence pris pour chacun des constituants est le constituant pur à l'état de gaz parfait à la même température que le mélange, mais à la pression <math>P^\circ</math>, prise égale à 1, exprimée dans la même unité que la pression <math>P</math> : bar, Pa, atm, mmHg, etc. Il est important de vérifier dans toute littérature donnant une constante d'équilibre l'unité dans laquelle les [[Pression partielle|pressions partielles]] sont exprimées. Si l'unité conseillée est le bar, et donc <math>P^\circ</math> = 1 bar, cette convention est loin d'être utilisée systématiquement.

L'expression de l'activité :

:<math>a_i = \frac{x_i P}{P^\circ}</math>

est alors réduite à la pression partielle du corps dans le mélange :

:<math>a_i = x_i P</math>

Rappelons que l'activité n'a pas de dimension, le fait de ne pas noter <math>P^\circ</math> n'est qu'une simplification d'écriture.

'''Pour une [[Solution (chimie)|solution liquide]]''' :

*Un '''[[solvant]]''', dont la concentration est largement supérieure à celle des [[soluté]]s, peut être assimilé à un corps pur et son activité est prise égale à 1 (par exemple l'eau dans une [[réaction d'oxydoréduction]] ou une [[réaction acido-basique]] en phase aqueuse).
*Lorsque la [[concentration molaire|concentration ou molarité]] d'un '''[[soluté]]''' est faible, son activité est prise égale au rapport de sa concentration <math>\left[ {\rm C}_i \right]</math> sur une concentration de référence <math>\left[ {\rm C} \right]^\circ</math> prise égale à 1, exprimée dans la même unité que la concentration réelle : par ex. mol.l<sup>−1</sup>. L'expression de l'activité :

:<math>a_i = \frac{\left[ {\rm C}_i \right]}{\left[ {\rm C} \right]^\circ}</math>

est alors réduite à la concentration du corps dans le mélange :

:<math>a_i = \left[ {\rm C}_i \right]</math>

L'activité peut également être exprimée à partir de la [[molalité]] <math>m_i</math> avec une molalité de référence <math>m^\circ</math> de 1 exprimée dans la même unité : par ex. mol.kg<sup>−1</sup>. L'expression de l'activité :

:<math>a_i = \frac{m_i}{m^\circ}</math>

est alors réduite à la molalité du corps dans le mélange :

:<math>a_i = m_i</math>

Ceci est notamment utilisé dans les calculs d'équilibres ioniques (Cf. [[Constante de dissociation]]) ou acido-basiques (Cf. [[Potentiel hydrogène]]) en phase aqueuse.

Rappelons que l'activité n'a pas de dimension, le fait de ne pas noter <math>\left[ {\rm C} \right]^\circ</math> ou <math>m^\circ</math> n'est qu'une simplification d'écriture.

'''Pour un [[État solide|solide]]''' l'activité est prise égale à 1 dans les constantes de [[solubilité]] lors des calculs d'équilibres d'espèces ioniques avec [[Précipité|précipitation]].

=== Calcul du coefficient d'activité ===

Contrairement au potentiel chimique qui n'est toujours défini qu'à une constante additive près, le coefficient d'activité, et l'activité, peut être calculé de façon absolue, ce qui le rend particulièrement intéressant dans le calcul des propriétés des mélanges réels.

Les modèles d'activité sont surtout utilisés dans le calcul des équilibres de phases.

De nombreux modèles de coefficients d'activité pour les phases liquides ont été développés, avec comme état standard les corps liquides purs :
* [[Modèle d'activité de Margules|Margules]] ;
* {{lien|lang=en|trad=Van Laar equation|texte=Van Laar}} ;
* [http://www.lassc.ulg.ac.be/webCheng00/chim007/Thermo07Activ.pdf Wilson] ;
* {{lien|lang=en|trad=Non-random two-liquid model|texte=NRTL}} ;
* [[UNIQUAC]] ;
* [[Unifac|UNIFAC]] ;
* [[Cosmospace|COSMOSPACE]].

Hormis les modèles de Margules et Van Laar, formes totalement empiriques à régresser sur des résultats expérimentaux, les autres modèles sont basés sur des théories d'interaction entre les divers corps du mélange et peuvent également être ajustés sur des résultats expérimentaux.

== Équilibre de phases ==
{{article détaillé|Équilibre de phases}}

=== Approche γ - ϕ ===

Lorsque deux phases ou plus (gaz, liquide, solide) sont en équilibre, les potentiels chimiques de chacun des corps en présence sont homogènes sur l'ensemble des phases. Autrement dit, s'il existe à l'équilibre les phases <math>\alpha, \beta, \gamma</math>, les [[potentiel chimique#Équilibre de phases|potentiels chimiques]] de tout corps <math>i</math> présent répondent à l'égalité :

<center>'''Équilibre chimique : '''<math>\mu_i^\alpha = \mu_i^\beta = \mu_i^\gamma</math></center>

Cette égalité est équivalente à celle des [[fugacité#Équilibre de phases|fugacités]] du corps <math>i</math> dans les différentes phases : <math>f_i^\alpha = f_i^\beta = f_i^\gamma</math>.

Nous supposerons par la suite un équilibre entre une phase vapeur et une phase liquide, en introduisant les fractions molaires respectives <math>x_i^\text{g}, x_i^\text{l}</math> de chacun des corps <math>i</math>.

Pour la phase vapeur, l'approche par coefficient de fugacité est préférée car il existe de nombreuses équations d'état pouvant la représenter correctement. L'état de référence pour le calcul du potentiel chimique <math>\mu_i^\text{g}</math> est le gaz parfait pur aux mêmes pression et température que le mélange réel. On a alors, en introduisant le [[coefficient de fugacité]] <math>\phi_i^\text{g}</math> en phase vapeur, l'expression de la fugacité en phase vapeur :

<center>'''Fugacité en phase vapeur : '''<math>f_i^\text{g} = x_i^\text{g} \phi_i^\text{g} P</math></center>

Pour la phase liquide, l'approche par coefficient d'activité est préférable à l'approche par coefficient de fugacité car il existe de nombreux modèles d'enthalpie libre d'excès qui la représentent bien mieux que les équations d'état. L'état de référence pour le calcul du potentiel chimique est le corps pur liquide aux mêmes pression et température que le mélange réel. On a alors, en introduisant le coefficient d'activité <math>\gamma_i^\text{l}</math> en phase liquide et la fugacité <math>f_i^\text{l,*}</math> du corps pur en phase liquide, l'expression de la fugacité en phase liquide :

<center>'''Fugacité en phase liquide : '''<math>f_i^\text{l} = x_i^\text{l} \gamma_i^\text{l} f_i^\text{l,*}</math></center>
<br>
{{boîte déroulante/début|titre = Note sur l'état standard de la phase liquide.}}

Pour la phase liquide, l'état standard choisi est celui du liquide pur. Ceci n'est pas incompatible avec la définition de la fugacité par rapport à l'état de gaz parfait pur. En effet, pour un même corps <math>i</math>, on a les relations successives des potentiels chimiques dans divers états aux mêmes <math>P</math> et <math>T</math> :

* <math>\mu_i^\text{l,*} - \mu_i^{\bullet,*} = RT \ln \! \left( {f_i^\text{l,*} \over P} \right)</math> pour le passage du gaz parfait pur au liquide pur ;
* <math>\mu_i^\text{l,id} - \mu_i^\text{l,*} = RT \ln x_i^\text{l}</math> pour le passage du liquide pur à la solution liquide idéale ;
* <math>\mu_i^\text{l} - \mu_i^\text{l,id} = RT \ln \gamma_i^\text{l}</math> pour le passage de la solution liquide idéale à la solution liquide réelle.

On a donc bien la relation liant le potentiel chimique du corps <math>i</math> dans la solution liquide réelle à celui du même corps <math>i</math> à l'état de gaz parfait pur aux mêmes <math>P</math> et <math>T</math> :

:<math>\mu_i^\text{l} - \mu_i^{\bullet,*} = RT \ln \! \left( {f_i^\text{l} \over P} \right) = \left[ \mu_i^\text{l} - \mu_i^\text{l,id} \right] + \left[ \mu_i^\text{l,id} - \mu_i^\text{l,*} \right] + \left[ \mu_i^\text{l,*} - \mu_i^{\bullet,*} \right] = RT \ln \! \left( {x_i^\text{l} \gamma_i^\text{l} f_i^\text{l,*} \over P} \right)</math>

{{boîte déroulante/fin}}

L'égalité des fugacités à l'équilibre des phases <math>f_i^\text{g} = f_i^\text{l}</math> induit l'expression de l''''approche <math>\gamma - \phi</math>''' (gamma-phi) :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Approche <math>\gamma - \phi</math> : ''' <math>x_i^\text{g} \phi_i^\text{g} P = x_i^\text{l} \gamma_i^\text{l} f_i^\text{l,*}</math>
|}

On déduit également de cette expression le '''[[coefficient de partage]]''' du corps <math>i</math> entre les deux phases, noté <math>K_i</math> :

<center>'''Coefficient de partage : ''' <math>K_i = \frac{x_i^\text{g}}{x_i^\text{l}} = \frac{\gamma_i^\text{l}}{\phi_i^\text{g}} \frac{f_i^\text{l,*}}{P}</math></center>

Cette approche par états de référence différents pour les deux phases est appelée « '''approche <math>\gamma - \phi</math>''' » (gamma-phi). L'approche par [[Fugacité|coefficient de fugacité]] employé pour les deux phases est appelée [[Fugacité#Équilibre de phases|« '''approche <math>\phi - \phi</math>''' »]] (phi-phi). Dans cette approche <math>f_i^\text{l} = x_i^\text{l} \phi_i^\text{l} P</math> ; on a donc : <math>\phi_i^\text{l} P = \gamma_i^\text{l} f_i^\text{l,*}</math>.

L'approche <math>\gamma - \phi</math> nécessite deux modèles distincts : une équation d'état pour la phase vapeur et un modèle de coefficient d'activité pour la phase liquide ; elle est recommandée pour des pressions inférieures à 10 bar. Au-delà, l'approche <math>\phi - \phi</math> est recommandée car elle ne nécessite qu'une seule et même équation d'état pour représenter les deux phases, ce qui rend le calcul cohérent à l'approche du point critique où les deux phases doivent se rejoindre. Les deux modèles de l'approche <math>\gamma - \phi</math> posent des problèmes de convergence numérique autour du point critique.

En ce qui concerne l'état standard, le corps <math>i</math> pur à la même température que le mélange réel n'existe pas toujours à l'état liquide, aussi faut-il discerner deux cas de figure :
* le corps <math>i</math> existe à l'état liquide pur à la température <math>T</math> du mélange, c'est-à-dire qu'il existe une pression de vapeur saturante <math>P_i^\text{sat}</math> du corps <math>i</math> pur à <math>T</math> : la fugacité <math>f_i^\text{l,*}</math> est calculée selon la « '''convention symétrique''' » ;
* le corps <math>i</math> est [[Fluide supercritique|supercritique]] à la température <math>T</math> du mélange, c'est-à-dire qu'il n'existe pas de pression de vapeur saturante <math>P_i^\text{sat}</math> du corps <math>i</math> pur à <math>T</math> car <math>T > T_{c,i}</math> la [[température critique]] du corps <math>i</math> : la fugacité <math>f_i^\text{l,*}</math> (fictive) est calculée selon la « '''convention dissymétrique''' », le corps <math>i</math> est alors un gaz dissout.

=== Convention symétrique ===

Lorsque tous les constituants de la phase liquide existent à l'état liquide à la température du mélange (c'est-à-dire que la température du mélange <math>T < T_{c,i}</math> la [[température critique]] du corps <math>i</math>, pour chacun des constituants), l'état de référence liquide pur est donc réel pour chacun des corps <math>i</math>. La fugacité de tous les corps <math>i</math> dans cet état est alors établie comme suit. Cette convention étant applicable à tous les constituants du mélange, elle est appelée « '''convention symétrique''' ».

{{boîte déroulante/début|titre = Démonstration.}}

On considère d'abord le liquide pur à la température <math>T</math> du mélange. La variation isotherme du potentiel chimique d'un corps <math>i</math> pur en phase liquide en fonction de la pression est donnée par la relation :

:<math>\left( \frac{\partial \mu_i^\text{l,*}}{\partial P} \right)_{T,n} = \bar V_i^\text{l,*}</math>

avec <math>\bar V_i^\text{l,*}</math> le [[grandeur molaire|volume molaire]] du corps pur à l'état liquide.

En intégrant l'expression ci-dessus entre <math>P_i^\text{sat}</math> et <math>P</math> la pression du mélange, on obtient :

:<math>\mu_i^\text{l,*} \! \left( P, T \right) - \mu_i^\text{l,*} \! \left( P_i^\text{sat}, T \right) = \mu_i^\text{l,*} - \mu_i^\text{l,*,sat} = \int_{P_i^\text{sat}}^{P} \bar V_i^\text{l,*} \, \mathrm{d} P</math>

Nous notons la correction de Poynting :

:<math>\mathcal{P}_i = \exp \! \left( \frac{ \int_{P_i^\text{sat}}^{P} \bar V_i^\text{l,*} \, \mathrm{d} P }{RT} \right)</math>
:{{equarefa|1}} <math>\mu_i^\text{l,*} - \mu_i^\text{l,*,sat} = RT \ln \! \left( \mathcal{P}_i \right)</math>

On considère ensuite le corps <math>i</math> pur à saturation à la température <math>T</math>, la pression est alors égale à la [[pression de vapeur saturante]] <math>P_i^\text{sat}</math> correspondante. On introduit le potentiel chimique <math>\mu_i^{\bullet,*,\text{sat}}</math> du corps <math>i</math> pur à l'état de gaz parfait, les expressions du potentiel chimique et du coefficient de fugacité du corps pur (soit <math>x_i^\text{g}=x_i^\text{l}=1</math> les fractions molaires respectives en phases vapeur et liquide) à saturation pour chacune des deux phases :

:<math>\mu_i^{\bullet,*} \! \left( P_i^\text{sat}, T \right) = \mu_i^{\bullet,*,\text{sat}}</math>
:<math>\mu_i^\text{g,*} \! \left( P_i^\text{sat}, T \right) = \mu_i^\text{g,*,sat} = \mu_i^{\bullet,*,\text{sat}} + RT \ln \phi_i^\text{g,*,sat}</math>
:<math>\phi_i^\text{g,*} \! \left( P_i^\text{sat}, T \right) = \phi_i^\text{g,*,sat}</math>
:<math>\mu_i^\text{l,*} \! \left( P_i^\text{sat}, T \right) = \mu_i^\text{l,*,sat} = \mu_i^{\bullet,*,\text{sat}} + RT \ln \phi_i^\text{l,*,sat}</math>
:<math>\phi_i^\text{l,*} \! \left( P_i^\text{sat}, T \right) = \phi_i^\text{l,*,sat}</math>

Puisque, par définition, à saturation le corps <math>i</math> liquide pur est en équilibre avec sa vapeur, le potentiel chimique du liquide pur est égal à celui du gaz pur :

:<math>\mu_i^\text{l,*,sat} = \mu_i^\text{g,*,sat}</math>

nous pouvons établir les égalités :

:<math>\phi_i^\text{g,*,sat} = \phi_i^\text{l,*,sat}</math>
:{{equarefa|2}} <math>\mu_i^\text{l,*,sat} - \mu_i^{\bullet,*,\text{sat}} = RT \ln \phi_i^\text{g,*,sat}</math>

'''Note''' : bien que les coefficients de fugacité à saturation des deux phases soient en théorie égaux, il est préférable de conserver le coefficient de fugacité de la phase vapeur dans l'expression des potentiels de la phase liquide. En pratique en effet les [[équation d'état|équations d'état]] représentent mieux les phases gaz que les phases liquides, aussi vaut-il mieux, pour la précision du résultat, calculer un coefficient de fugacité d'une phase vapeur que d'une phase liquide.

Enfin, en considérant la variation du potentiel chimique d'un gaz parfait pur en fonction de la pression à température constante :

:<math>\left( \frac{\partial \mu_i^{\bullet,*}}{\partial P} \right)_{T,n} = \bar V_i^{\bullet,*} = \frac{RT}{P}</math>
:<math>\mu_i^{\bullet,*} \! \left( P_i^\text{sat}, T \right) - \mu_i^{\bullet,*} \! \left( P, T \right) = \int_{P}^{P_i^\text{sat}} \frac{RT}{P} \, \mathrm{d} P</math>
:{{equarefa|3}} <math>\mu_i^{\bullet,*,\text{sat}} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{P_i^\text{sat}}{P} \right)</math>

Nous avons donc les relations :
:{{equarefl|1}} <math>\mu_i^\text{l,*} - \mu_i^\text{l,*,sat} = RT \ln \mathcal{P}_i</math>
:{{equarefl|2}} <math>\mu_i^\text{l,*,sat} - \mu_i^{\bullet,*,\text{sat}} = RT \ln \phi_i^\text{g,*,sat}</math>
:{{equarefl|3}} <math>\mu_i^{\bullet,*,\text{sat}} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{P_i^\text{sat}}{P} \right)</math>

Nous obtenons :

:<math>\mu_i^\text{l,*} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{f_i^\text{l,*}}{P} \right) = RT \ln \! \left( \phi_i^\text{g,*,sat} \frac{P_i^\text{sat}}{P} \mathcal{P}_i \right)</math>

{{boîte déroulante/fin}}

La fugacité du corps <math>i</math> pur à l'état de liquide à <math>P</math> et <math>T</math> vaut :

<center>'''Fugacité du corps <math>i</math> liquide pur : '''<math>f_i^\text{l,*} = \phi_i^\text{g,*,sat} P_i^\text{sat} \mathcal{P}_i</math></center>

L'égalité des potentiels chimiques en phase vapeur et en phase liquide conduit selon la '''convention symétrique''' à :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Pour un solvant : ''' <math>x_i^\text{g} \phi_i^\text{g} P = x_i^\text{l} \gamma_i^\text{l} \phi_i^\text{g,*,sat} P_i^\text{sat} \mathcal{P}_i</math>
|}

On déduit également de cette expression le '''[[coefficient de partage]]''' du corps <math>i</math> entre les deux phases, noté <math>K_i</math> :

<center>'''Coefficient de partage : ''' <math>K_i = \frac{x_i^\text{g}}{x_i^\text{l}} = \gamma_i^\text{l} \frac{\phi_i^\text{g,*,sat}}{\phi_i^\text{g}} \frac{P_i^\text{sat}}{P} \mathcal{P}_i</math></center>

avec :
* <math>T</math> la température d'équilibre ;
* <math>P</math> la pression d'équilibre ;
* <math>P_i^\text{sat}</math> la [[pression de vapeur saturante]] du corps <math>i</math> pur à <math>T</math> ;
* <math>x_i^\text{g}</math> la fraction molaire du corps <math>i</math> en phase vapeur ;
* <math>x_i^\text{l}</math> la fraction molaire du corps <math>i</math> en phase liquide ;
* <math>\phi_i^\text{g}</math> le coefficient de fugacité du corps <math>i</math> en phase vapeur, à <math>P</math>, <math>T</math> et composition du mélange gazeux ;
* <math>\phi_i^\text{g,*,sat}</math> le coefficient de fugacité du corps <math>i</math> pur à saturation en phase vapeur, à <math>P_i^\text{sat}</math> et <math>T</math> ;
* <math>\gamma_i^\text{l}</math> le coefficient d'activité du corps <math>i</math> en phase liquide, à <math>P</math>, <math>T</math> et composition du mélange liquide ;
* <math>\bar V_i^\text{l,*}</math> le volume molaire liquide du corps <math>i</math> pur à <math>T</math> ;
* <math>\mathcal{P}_i = \exp \! \left( \frac{ \int_{P_i^\text{sat}}^{P} \bar V_i^\text{l,*} \, \mathrm{d} P }{RT} \right)</math> la correction de Poynting.

Les liquides étant peu compressibles, nous pouvons supposer que le volume molaire du liquide pur ne dépend que de la température, aussi obtient-on :

:<math>\int_{P_i^\text{sat}}^{P} \bar V_i^\text{l,*} \, \mathrm{d} P \approx \bar V_i^\text{l,*} \left( P - P_i^\text{sat} \right)</math>
:<math>\mathcal{P}_i \approx \exp \! \left( \frac{ \bar V_i^\text{l,*} \left( P - P_i^\text{sat} \right) }{RT} \right)</math>

Pour des pressions de l'ordre de grandeur de la [[pression atmosphérique]], la correction de Poynting est négligeable : <math>\mathcal{P}_i \approx 1</math>. Si l'on considère la phase vapeur comme un [[mélange de gaz parfaits]], alors <math>\phi_i^\text{g} = 1</math> et <math>\phi_i^\text{g,*,sat} = 1</math>. De même, si l'on considère la phase liquide comme une solution idéale, alors <math>\gamma_i^\text{l} = 1</math>. On obtient la '''[[loi de Raoult]]''' :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Loi de Raoult : ''' <math>x_i^\text{g} P = x_i^\text{l} P_i^\text{sat}</math>
|}

=== Convention dissymétrique ===
<!-- "dissout" et non "dissous", réforme de l'orthographe de 1990 -->
{{article détaillé|Loi de Henry}}

Lorsqu'un constituant <math>i</math> dans la phase liquide est [[Fluide supercritique|supercritique]] (c'est-à-dire que la température du mélange <math>T > T_{c,i}</math> la [[température critique]] du constituant <math>i</math>), le corps <math>i</math> est un gaz dissout dans un solvant liquide. L'état de référence liquide pur du constituant <math>i</math> est alors fictif : on se base par conséquent sur un autre état de référence, celui du corps <math>i</math> infiniment dilué dans le solvant aux mêmes <math>P</math> et <math>T</math>.

Pour les autres corps <math>k</math> tels que <math>T < T_{c,k}</math>, ces corps sont considérés comme des solvants et on leur applique la référence liquide pur de la convention symétrique.

Tous les corps dans le mélange liquide n'étant pas représentés de la même façon, puisque l'on distingue les gaz dissouts des solvants, cette convention est appelée '''convention dissymétrique'''.

Avant de développer le formalisme mathématique de la convention dissymétrique, il est nécessaire de définir la '''constante de Henry'''.

==== Constante de Henry ====
<!-- "dissout" et non "dissous", réforme de l'orthographe de 1990 -->

Soit <math>f_{i,s}^\infty</math> la fugacité du corps <math>i</math> à dilution infinie dans un solvant <math>s</math> :

:<math>\mu_{i,s}^{\infty} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{f_{i,s}^\infty}{P} \right)</math>

À dilution infinie d'un corps <math>i</math> dans un solvant <math>s</math>, le potentiel chimique <math>\mu_{i,s}^\infty</math> de <math>i</math> est [[Potentiel chimique|infiniment négatif]] :

:<math>\mu_{i,s}^\infty = -\infty</math>

la fugacité à dilution infinie ne peut donc être que nulle :

:<math>f_{i,s}^\infty = \lim_{x_i^\text{l} \to 0} f_i^\text{l} = 0</math>

Si l'on considère que, pour des fractions molaires de <math>i</math> proches de 0, la fugacité se comporte idéalement de façon linéaire, telle que :

:<math>f_i^\text{l} = x_i^\text{l} k_{\text{H},i,s}</math>

bien que <math>\lim_{x_i^\text{l} \to 0} f_i^\text{l} = 0</math> et <math>\lim_{x_i^\text{l} \to 0} x_i^\text{l} = 0</math>, en vertu de la [[règle de L'Hôpital]] on définit la '''constante de Henry''' <math>k_{\text{H},i,s}</math> par la limite :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Constante de Henry : ''' <math>k_{\text{H},i,s} = \lim_{x_i^\text{l} \to 0} \frac{f_i^\text{l}}{x_i^\text{l}}</math>
|}

La constante de Henry, comme la fugacité, a la dimension d'une pression.

La constante de Henry ''n'est pas'' la fugacité du corps <math>i</math> à dilution infinie : <math>k_{\text{H},i,s} \neq f_{i,s}^\infty</math>. Comme vu plus haut celle-ci est nulle : <math>f_{i,s}^\infty = 0</math>.

La constante de Henry dépend du solvant <math>s</math> dans lequel le corps <math>i</math> est dissout, ce solvant pouvant être aussi bien un corps pur qu'un mélange.

En considérant la solution du corps <math>i</math> gazeux dissout dans le solvant liquide <math>s</math> comme une [[solution idéale]], on écrit :

:<math>\mu_{i,s}^{\infty,\text{id}} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{x_i^\text{l} k_{\text{H},i,s}}{P} \right)</math>

avec <math>\mu_{i,s}^{\infty,\text{id}}</math> le potentiel chimique du corps <math>i</math> en solution idéale dans le solvant <math>s</math>. Les potentiels chimiques de la relation précédente sont définis aux mêmes <math>P</math> et <math>T</math>.

La constante de Henry est déterminée expérimentalement en fonction d'une pression de référence <math>P^{\text{réf}}</math> et de la température <math>T</math>. La pression de référence <math>P^{\text{réf}}</math> est souvent prise égale à la pression de vapeur saturante du solvant <math>s</math>, <math>P_s^\text{sat}</math>, à <math>T</math> donnée ; dans ce cas, la pression de référence varie avec la température. La pression de référence peut également être fixe (le plus souvent 1 atm).

La constante de Henry peut être calculée à une autre pression <math>P</math>, à la même température <math>T</math>, selon la relation :

:<math>\left( \frac{\partial RT \ln \! \left( k_{\text{H},i,s} \! \left( P,T \right) \right)}{\partial P} \right)_{T,n} = \bar V_{i,s}^{\infty}</math>
:<math>RT \ln \! \left( k_{\text{H},i,s} \! \left( P,T \right) \right) - RT \ln \! \left( k_{\text{H},i,s} \! \left( P^{\text{réf}},T \right) \right) = \int_{P^{\text{réf}}}^{P} \bar V_{i,s}^{\infty} \, \mathrm{d} P</math>

avec <math>\bar V_{i,s}^{\infty}</math> le [[grandeur molaire partielle|volume molaire partiel]] du corps <math>i</math> infiniment dilué dans le solvant <math>s</math>. Cette grandeur, non nulle, est déterminée expérimentalement par extrapolation du volume molaire partiel du corps <math>i</math> dans le mélange liquide :

:<math>\bar V_{i,s}^{\infty} = \lim_{x_i^\text{l} \to 0} \bar V_i^\text{l}</math>

Nous notons la correction de Poynting :

:<math>\mathcal{P}_{i,s}^{\infty} = \exp \! \left( \frac{ \int_{P^{\text{réf}}}^{P} \bar V_{i,s}^{\infty} \, \mathrm{d} P }{RT} \right)</math>

D'où l'évolution de la constante de Henry entre deux pressions :

:<math>k_{\text{H},i,s} \! \left( P,T \right) = k_{\text{H},i,s} \! \left( P^{\text{réf}},T \right) \mathcal{P}_{i,s}^{\infty}</math>

<center>'''Variation de la constante de Henry avec la pression : '''<math>k_{\text{H},i,s} = k_{\text{H},i,s}^{\text{réf}} \mathcal{P}_{i,s}^{\infty}</math></center>

En supposant que le volume molaire <math>\bar V_{i,s}^{\infty}</math> ne dépend pas de la pression, cette relation devient l''''équation de Krichevsky–Kasarnovsky''' :

<center>'''Équation de Krichevsky–Kasarnovsky : '''<math>\ln k_{\text{H},i,s} \! \left( P,T \right) = \ln k_{\text{H},i,s} \! \left( P^\text{réf},T \right) + {\bar V_{i,s}^{\infty} \left( P - P^\text{réf} \right) \over RT}</math></center>

==== Formalisme ====
<!-- "dissout" et non "dissous", réforme de l'orthographe de 1990 -->

{{boîte déroulante/début|titre = Démonstration.}}

L'état de référence considéré est celui du corps <math>i</math> infiniment dilué dans le solvant <math>s</math>. Nous introduisons <math>\gamma_i^{\oslash,\text{l}}</math> le coefficient d'activité du corps <math>i</math> en solution dans le solvant <math>s</math> à <math>P</math> et <math>T</math> dans la convention dissymétrique, correspondant à l'écart entre le potentiel chimique <math>\mu_i^\text{l}</math> du corps <math>i</math> dans le mélange liquide réel et le potentiel chimique <math>\mu_{i,s}^{\infty,\text{id}}</math> du corps <math>i</math> en solution idéale dans le solvant <math>s</math>.

Pour un corps <math>i</math> dans la convention dissymétrique, la fugacité en phase liquide s'écrit :

:{{equarefa|a}} <math>\mu_i^\text{l} - \mu_{i,s}^{\infty,\text{id}} = RT \ln \gamma_i^{\oslash,\text{l}}</math>
:<math>\mu_{i,s}^{\infty,\text{id}} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{ x_i^\text{l} k_{\text{H},i,s}}{P} \right)</math>
:<math>\mu_i^\text{l} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{x_i^\text{l} \gamma_i^{\oslash,\text{l}} k_{\text{H},i,s}}{P} \right)</math>
:{{equarefa|b}} <math>f_i^\text{l} = x_i^\text{l} \gamma_i^{\oslash,\text{l}} k_{\text{H},i,s}</math>

Dans cette convention, par définition de la constante de Henry :

:<math>\lim_{ x_i^\text{l} \to 0 } \frac{f_i^\text{l}}{x_i^\text{l}} = k_{\text{H},i,s}</math>

le passage à la limite de la dilution infinie de la relation {{equarefl|b}} donne :

:{{equarefa|c}} <math>\lim_{ x_i^\text{l} \to 0 } \gamma_i^{\oslash,\text{l}} = 1</math>

Pour un corps <math>i</math> dans la convention symétrique, la fugacité en phase liquide s'écrit :

:{{equarefa|d}} <math>\mu_i^\text{l} - \mu_i^\text{l,id} = RT \ln \gamma_i^\text{l}</math>
:<math>\mu_i^\text{l,id} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{x_i^\text{l} f_i^\text{l,*}}{P} \right)</math>
:<math>\mu_i^\text{l} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{x_i^\text{l} \gamma_i f_i^\text{l,*}}{P} \right)</math>
:{{equarefa|e}} <math>f_i^\text{l} = x_i^\text{l} \gamma_i^\text{l} f_i^\text{l,*}</math>

Dans cette convention, par définition de la fugacité du corps <math>i</math> pur en phase liquide :

:<math>\lim_{ x_i^\text{l} \to 1 } f_i^\text{l} = f_i^\text{l,*}</math>

le passage à la limite du corps pur de la relation {{equarefl|e}} donne :

:{{equarefa|f}} <math>\lim_{ x_i^\text{l} \to 1 } \gamma_i^\text{l} = 1</math>

Les limites {{equarefl|c}} et {{equarefl|f}} différentes sont dues aux états de référence différents. Les modèles de coefficients d'activité exposés au paragraphe [[#Calcul du coefficient d'activité|calcul du coefficient d'activité]] sont développés pour l'état standard corps liquide pur et <math>\gamma_i^\text{l}</math>. Néanmoins, ils sont également utilisés pour calculer le coefficient d'activité <math>\gamma_i^{\oslash,\text{l}}</math> comme démontré plus bas.

L'égalité des expressions {{equarefl|b}} et {{equarefl|e}} de la fugacité en phase liquide donne :

:<math>f_i^\text{l} = x_i^\text{l} \gamma_i^\text{l} f_i^\text{l,*} = x_i^\text{l} \gamma_i^{\oslash,\text{l}} k_{\text{H},i,s}</math>
:{{equarefa|g}} <math>\gamma_i^\text{l} f_i^\text{l,*} = \gamma_i^{\oslash,\text{l}} k_{\text{H},i,s}</math>

Pour un gaz dissout, la fugacité à l'état de liquide pur <math>f_i^\text{l,*}</math> est purement fictive.

En passant l'égalité {{equarefl|g}} à la limite de la dilution infinie :

:<math>\lim_{ x_i^\text{l} \to 0 } \left[ \gamma_i^\text{l} f_i^\text{l,*} \right] = \lim_{ x_i^\text{l} \to 0 } \left[ \gamma_i^{\oslash,\text{l}} k_{\text{H},i,s} \right]</math>
:<math>\left[ \lim_{ x_i^\text{l} \to 0 } \gamma_i^\text{l} \right] f_i^\text{l,*} = \left[ \lim_{ x_i^\text{l} \to 0 } \gamma_i^{\oslash,\text{l}} \right] k_{\text{H},i,s}</math>

car <math>f_i^\text{l,*}</math> définie à <math>x_i^\text{l} = 1</math> et <math>k_{\text{H},i,s}</math> définie à <math>x_i^\text{l} = 0</math> ne dépendent pas de <math>x_i^\text{l}</math>.

En introduisant le coefficient d'activité du corps <math>i</math> à dilution infinie défini par :

:<math>\gamma_{i,s}^{\infty} = \lim_{ x_i^\text{l} \to 0 } \gamma_i^\text{l}</math>

on obtient, avec {{equarefl|c}} :

:{{equarefa|h}} <math>\gamma_{i,s}^{\infty} f_i^\text{l,*} = k_{\text{H},i,s}</math>

En introduisant cette relation {{equarefl|h}} dans l'égalité {{equarefl|g}}, on obtient finalement :

:{{equarefa|i}} <math>\gamma_i^{\oslash,\text{l}} = \frac{\gamma_i^\text{l}}{\gamma_{i,s}^{\infty}}</math>

Les coefficients d'activité <math>\gamma_i^\text{l}</math> et <math>\gamma_{i,s}^{\infty}</math> sont calculés au moyen des modèles classiques exposés au paragraphe [[#Calcul du coefficient d'activité|calcul du coefficient d'activité]], il n'y a donc pas besoin de développer des modèles spécifiques de coefficient d'activité pour la convention dissymétrique et <math>\gamma_i^{\oslash,\text{l}}</math>. On vérifie en effet la relation {{equarefl|c}} :

:<math>\lim_{ x_i^\text{l} \to 0 } \gamma_i^{\oslash,\text{l}} = \frac{ 1 }{\gamma_{i,s}^{\infty}} . \lim_{ x_i^\text{l} \to 0 } \gamma_i^\text{l} = 1</math>

De plus, avec les relations {{equarefl|a}}, {{equarefl|d}} et {{equarefl|i}}, nous avons :

:<math>\mu_{i,s}^{\infty,\text{id}} - \mu_i^\text{l,id} = RT \ln \gamma_{i,s}^{\infty}</math>

Le coefficient d'activité à dilution infinie représente donc l'écart entre les deux solutions idéales de référence.

La relation {{equarefl|h}} donne :

:<math>\mu_i^\text{l,*} - \mu_i^{\bullet,*} = RT \ln \! \left( \frac{f_i^\text{l,*}}{P} \right) = RT \ln \! \left( \frac{1}{\gamma_{i,s}^{\infty}} \frac{k_{\text{H},i,s}}{P} \right)</math>

{{boîte déroulante/fin}}

La fugacité fictive du corps <math>i</math> à l'état de liquide pur à <math>P</math> et <math>T</math> vaut :

<center>'''Fugacité du corps <math>i</math> liquide pur : '''<math>f_i^\text{l,*} = \frac{1}{\gamma_{i,s}^{\infty}} k_{\text{H},i,s} = \frac{1}{\gamma_{i,s}^{\infty}} k_{\text{H},i,s}^{\text{réf}} \mathcal{P}_{i,s}^{\infty}</math></center>

L'égalité des potentiels chimiques en phase vapeur et en phase liquide conduit selon la '''convention dissymétrique''' à :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Pour un gaz dissout :''' <math>x_i^\text{g} \phi_i^\text{g} P = x_i^\text{l} \frac{\gamma_i^\text{l}}{\gamma_{i,s}^{\infty}} k_{\text{H},i,s}^{\text{réf}} \mathcal{P}_{i,s}^{\infty}</math></center>
|}

On déduit également de cette expression le '''[[coefficient de partage]]''' du corps <math>i</math> entre les deux phases, noté <math>K_i</math> :

<center>'''Coefficient de partage :''' <math>K_i = \frac{x_i^\text{g}}{x_i^\text{l}} = \frac{1}{\phi_i^\text{g}} \frac{\gamma_i^\text{l}}{\gamma_{i,s}^{\infty}} \frac{k_{\text{H},i,s}^{\text{réf}}}{P} \mathcal{P}_{i,s}^{\infty}</math></center>

avec :
* <math>T</math> la température d'équilibre ;
* <math>P</math> la pression d'équilibre ;
* <math>P^{\text{réf}}</math> la pression à laquelle la constante de Henry a été expérimentalement déterminée à <math>T</math> ;
* <math>k_{\text{H},i,s}^{\text{réf}}</math> la [[#Constante de Henry|constante de Henry]] du corps <math>i</math> à <math>P^{\text{réf}}</math> et <math>T</math> dans le solvant <math>s</math> ;
* <math>x_i^\text{g}</math> la fraction molaire du corps <math>i</math> en phase vapeur ;
* <math>x_i^\text{l}</math> la fraction molaire du corps <math>i</math> en phase liquide ;
* <math>\phi_i^\text{g}</math> le coefficient de fugacité du corps <math>i</math> en phase vapeur, à <math>P</math>, <math>T</math> et composition du mélange gazeux ;
* <math>\gamma_i^\text{l}</math> le coefficient d'activité du corps <math>i</math> en phase liquide, à <math>P</math>, <math>T</math> et composition du mélange liquide ;
* <math>\gamma_{i,s}^{\infty} = \lim_{ x_i^\text{l} \to 0 } \gamma_i^\text{l}</math> le coefficient d'activité du corps <math>i</math> en phase liquide, à <math>P</math>, <math>T</math> et dilution infinie dans le solvant <math>s</math> ;
* <math>\bar V_{i,s}^{\infty} = \lim_{x_i^\text{l} \to 0} \bar V_i^\text{l}</math> le volume molaire partiel du corps <math>i</math> en phase liquide, à <math>P</math>, <math>T</math> et dilution infinie dans le solvant <math>s</math> ;
* <math>\mathcal{P}_{i,s}^{\infty} = \exp \! \left( \frac{ \int_{P^{\text{réf}}}^{P} \bar V_{i,s}^{\infty} \, \mathrm{d} P }{RT} \right)</math> la correction de Poynting.

Les liquides étant peu compressibles, nous supposons que le volume molaire ne dépend que de la température, aussi obtient-on :

:<math>\int_{P^{\text{réf}}}^{P} \bar V_{i,s}^{\infty} \, \mathrm{d} P \approx \bar V_{i,s}^{\infty} \left( P - P^{\text{réf}} \right)</math>
:<math>\mathcal{P}_{i,s}^{\infty} \approx \exp \! \left( \frac{ \bar V_{i,s}^{\infty} \left( P - P^{\text{réf}} \right) }{RT} \right)</math>

Pour des pressions de l'ordre de grandeur de la pression atmosphérique, la correction de Poynting est négligeable : <math>\mathcal{P}_{i,s}^{\infty} \approx 1</math>. Si l'on considère la phase vapeur comme un mélange de gaz parfaits, alors <math>\phi_i^\text{g} = 1</math>. De même, si l'on considère la phase liquide comme une solution idéale, alors <math>\frac{\gamma_i^\text{l}}{\gamma_{i,s}^{\infty}} = 1</math>. On obtient la '''[[loi de Henry]]''' :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Loi de Henry :''' <math>x_i^\text{g} P = x_i^\text{l} k_{\text{H},i,s}^{\text{réf}}</math>
|}

== Grandeurs extensives des solutions réelles ==
=== Grandeurs d'excès ===

Une '''grandeur d'excès''' <math>X^\text{E}</math> exprime l'écart entre une grandeur thermodynamique [[Extensivité et intensivité (physique)|extensive]] <math>X^\text{l}</math> d'un mélange liquide réel et la même grandeur thermodynamique extensive <math>X^\text{l,id}</math> d'un [[solution idéale|mélange liquide idéal]] aux mêmes pression, température et composition :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Grandeur d'excès : ''' <math>X^\text{E} = X^\text{l} - X^\text{l,id}</math>
|}

Une grandeur d'excès est par construction également une grandeur [[Extensivité et intensivité (physique)|extensive]], il est donc possible de définir une [[grandeur molaire]] d'excès <math>\bar X^\text{E}</math> pour le mélange et une [[grandeur molaire partielle]] d'excès <math>\bar X_i^\text{E}</math> pour chacun des corps <math>i</math> présents dans le mélange. Les relations suivantes sont également vraies :

<center>'''Grandeur molaire d'excès : '''<math>\bar X^\text{E} = {X^\text{E} \over n} = \bar X - \bar X^\text{l,id}</math></center>
<center>'''Grandeur molaire partielle d'excès : '''<math>\bar X_i^\text{E}
= \left( {\partial X^\text{E} \over \partial n_i} \right)_{P,T,n_{j \neq i}} = \bar X_i - \bar X_i^\text{l,id}</math></center>

avec, toutes ces grandeurs étant définies aux mêmes pression, température et composition :
* <math>n</math> la [[quantité de matière]] totale dans le mélange ;
* <math>n_i</math> la quantité du corps <math>i</math> dans le mélange ;
* <math>\bar X</math> la grandeur molaire du mélange réel ;
* <math>\bar X_i</math> la grandeur molaire partielle du corps <math>i</math> dans le mélange réel ;
* <math>\bar X^\text{E}</math> la grandeur molaire d'excès du mélange ;
* <math>\bar X_i^\text{E}</math> la grandeur molaire partielle d'excès du corps <math>i</math> ;
* <math>\bar X^\text{l,id}</math> la grandeur molaire du mélange liquide idéal ;
* <math>\bar X_i^\text{l,id}</math> la grandeur molaire partielle du corps <math>i</math> dans le mélange liquide idéal.

''Note : la même démarche est en théorie applicable à toutes les phases, gaz, liquide ou solide. Néanmoins, elle est surtout employée pour les phases condensées liquide et solide. Dans ce dernier cas il faut bien entendu considérer une solution idéale solide, constituée à partir des solides purs, et dans ce qui suit un modèle de coefficient d'activité spécifique aux solides.''

=== Calcul des grandeurs d'excès ===

Pour l'[[enthalpie libre]] <math>G</math>, l'identité des [[potentiel chimique|potentiels chimiques]] <math>\mu_i</math> et des [[grandeur molaire partielle|enthalpies libres molaires partielles]] <math>\bar G_i</math> permet d'écrire, pour tout corps <math>i</math> dans un mélange liquide quelconque :

:<math>\bar G_i^\text{l} = \bar G_i^\text{l,id} + RT \ln \gamma_i^\text{l} = \mu_i^\text{l} = \mu_i^\text{l,id} + RT \ln \gamma_i^\text{l}</math>

Le terme du coefficient d'activité est donc l''''enthalpie libre molaire partielle d'excès''', ou '''potentiel chimique d'excès''' :

<center>'''Enthalpie libre molaire partielle d'excès : '''<math>\bar G_i^\text{E} = \mu_i^\text{E} = RT \ln \gamma_i^\text{l}</math></center>

À partir de l'enthalpie libre d'excès, toutes les autres grandeurs d'excès peuvent être calculées :
* <math>\bar S_i^\text{E} = -\left( \frac{\partial \bar G_i^\text{E}}{\partial T} \right)_{P,x}</math> l'[[Entropie (thermodynamique)|entropie]] molaire partielle d'excès, calculée par l'une des [[Relations de Maxwell#Équations d'état|équations d'état]] ;
* <math>\bar H_i^\text{E} = \bar G_i^\text{E} + T \bar S_i^\text{E}</math> l'[[enthalpie]] molaire partielle d'excès ;
* <math>\bar V_i^\text{E} = \left(\frac{\partial \bar G_i^\text{E}}{\partial P} \right)_{T,x}</math> le volume molaire partiel d'excès, calculé par l'une des [[Relations de Maxwell#Équations d'état|équations d'état]] ;
* <math>\bar U_i^\text{E} = \bar H_i^\text{E} -P \bar V_i^\text{E}</math> l'[[énergie interne]] molaire partielle d'excès ;
* <math>\bar F_i^\text{E} = \bar G_i^\text{E} -P \bar V_i^\text{E}</math> l'[[énergie libre]] molaire partielle d'excès.

Le [[Théorème d'Euler (fonctions de plusieurs variables)|théorème d'Euler sur les fonctions homogènes du premier ordre]] s'applique aux grandeurs d'excès puisque ce sont des grandeurs extensives, avec <math>n_i^\text{l}</math> [[quantité de matière|nombre de moles]] et <math>x_i^\text{l}</math> [[fraction molaire]] de chacun des <math>N</math> corps <math>i</math> présents dans le mélange liquide :

<center>'''Grandeur d'excès : '''<math>X^\text{E} = \sum_{i=1}^{N} n_i^\text{l} \bar X_i^\text{E}</math></center>
<center>'''Grandeur molaire d'excès : '''<math>\bar X^\text{E} = \sum_{i=1}^{N} x_i^\text{l} \bar X_i^\text{E}</math></center>

L'enthalpie libre d'excès a donc une importance particulière dans le calcul des propriétés des phases condensées :

<center>'''Enthalpie libre d'excès : '''<math>G^\text{E} = \sum_{i=1}^{N} n_i^\text{l} \bar G_i^\text{E} = RT \sum_{i=1}^{N} n_i^\text{l} \ln \gamma_i^\text{l}</math></center>
<center>'''Enthalpie libre molaire d'excès : '''<math>\bar G^\text{E} = \sum_{i=1}^{N} x_i^\text{l} \bar G_i^\text{E} = RT \sum_{i=1}^{N} x_i^\text{l} \ln \gamma_i^\text{l}</math></center>

'''Note sur la dépendance de <math>G^\text{E}</math> à la pression'''

:Si on mélange 1 litre d'eau avec 1 litre d'[[éthanol#Miscibilité avec l'eau|éthanol]], on obtient un volume total <math>V</math> d'environ {{unité|1.92|litre}}<ref>[https://www.inrs.fr/dms/inrs/FicheToxicologique/TI-FT-48/ft48.pdf Fiche INRS de l'éthanol].</ref>. Le volume idéal <math>V^\text{id}</math> est de {{unité|1|litre}} + {{unité|1|litre}} = {{unité|2|litres}} et le volume d'excès <math>V^\text{E}</math> est de {{unité|1,92|litre}} - {{unité|2|litres}} = {{unité|-0,08|litre}}, il y a contraction du volume. Le mélange eau-éthanol n'est donc pas une [[solution idéale]], il présente d'ailleurs un [[azéotrope]] que la [[loi de Raoult]] ne peut pas représenter.

:Dans les phases condensées, le volume d'excès <math>V^\text{E} = \left(\frac{\partial G^\text{E}}{\partial P} \right)_{T,x}</math> est souvent négligeable et négligé devant le volume de la solution idéale <math>V^\text{id}</math>. En pratique, aucun des modèles de coefficient d'activité (ou d'enthalpie libre molaire d'excès <math>\bar G^\text{E}</math>) listés au paragraphe ''[[#Calcul du coefficient d'activité|Calcul du coefficient d'activité]]'' ne dépend de <math>P</math>, et <math>V^\text{E}</math> est souvent considéré comme nul.

=== Calcul des grandeurs extensives réelles ===

Les grandeurs extensives <math>X^\text{l,id}</math> du mélange liquide idéal peuvent être calculées facilement à l'aide des propriétés des corps liquides purs si l'on considère une [[solution idéale]] aux mêmes pression, température et composition que le mélange réel.

Les grandeurs extensives du mélange liquide réel se calculent donc en appliquant la définition des grandeurs d'excès :

{|border="1" cellpadding="5" style="text-align:center;" align="center"
|'''Propriétés d'un mélange liquide réel : ''' <math>X^\text{l} = X^\text{l,id} + X^\text{E}</math>
|}

Ceci est appliqué surtout dans le calcul des propriétés des phases condensées, liquide ou solide (pour peu que l'on dispose d'un modèle de coefficients d'activité des solides). Pour les gaz, la solution idéale choisie est le mélange de [[gaz parfait]]s aux mêmes pression, température et composition que le mélange réel, et les propriétés sont calculées à l'aide d'une [[équation d'état]] et de [[grandeur résiduelle|grandeurs résiduelles]].

== Notes et références ==
=== Notes ===
{{Références}}

=== Bibliographie ===

;Articles
* {{Article |langue=fr |auteur1=E. Darmois |titre=La thermodynamique des solutions |périodique=J. Phys. Radium |volume=4 |numéro=7 |date=1943 |pages=129-142 |issn=|e-issn=|lire en ligne=https://hal.archives-ouvertes.fr/jpa-00233838/document |consulté le=13 juillet 2022 |id=RefDarmois}}.

;Ouvrages
* {{ouvrage |langue=fr |auteur1=Christophe Coquelet |auteur2=Dominique Richon |titre=Propriétés thermodynamiques |sous-titre=Détermination pour les fluides purs |éditeur=[[Éditions techniques de l'ingénieur |Techniques de l'ingénieur]] |collection=base documentaire ''Thermodynamique et énergétique'', pack ''Physique énergétique'', univers : ''Énergies'' |volume=BE 8030 |passage=1-8 |année=2007}}.
* {{Ouvrage |langue=fr |auteur1=Christophe Coquelet |auteur2=Dominique Richon |titre=Propriétés thermodynamiques |sous-titre=Détermination pour les mélanges |volume=BE 8 031 |éditeur=[[Éditions techniques de l'ingénieur |Techniques de l'ingénieur]] |année=2007 |pages totales=12 |passage=6 |isbn=|lire en ligne=https://books.google.fr/books?id=z9_vNuYiOiQC&pg=PA6 |id=Coquelet2007}}.
* {{Ouvrage |langue=fr |auteur1=Jean-Pierre Corriou |titre=Thermodynamique chimique |sous-titre=Définitions et relations fondamentales |éditeur=[[Éditions techniques de l'ingénieur |Techniques de l'ingénieur]] |collection=base documentaire : ''Thermodynamique et cinétique chimique'', pack ''Opérations unitaires. Génie de la réaction chimique'', univers ''Procédés chimie - bio - agro'' |lieu=|année=1984 |volume=J 1025 |tome=|pages totales=|passage=1-19 |isbn=|lire en ligne=https://books.google.fr/books?id=SIQEilUn97sC&pg=PP1 |id=Corriou1984}}.
* {{Ouvrage |langue=fr |auteur1=Jean-Pierre Corriou |titre=Thermodynamique chimique |sous-titre=Diagrammes thermodynamiques |volume=J 1026 |éditeur=[[Éditions techniques de l'ingénieur |Techniques de l'ingénieur]] |collection=base documentaire ''Thermodynamique et cinétique chimique'', pack ''Opérations unitaires. Génie de la réaction chimique'', univers ''Procédés chimie - bio - agro'' |année=1985 |passage=1-30 |isbn=}}.
* {{Ouvrage |langue=en |auteur1=Robert C. Reid |auteur2=John M. Prausnitz |auteur3=Bruce E. Poling |titre=The properties of gases and liquids |lieu=New York |éditeur=McGraw-Hill |année=1987 |numéro d'édition=4 |pages totales=741 |isbn=978-0-07-051799-8}}.
* {{Ouvrage |langue=fr |prénom1=Jean |nom1=Vidal |titre=Thermodynamique : application au génie chimique et à l'industrie pétrolière |éditeur=[[Technip |Éditions Technip]] |collection=Publications de l'[[IFP Énergies nouvelles |Institut français du pétrole]]. |lieu=Paris |année=1997 |pages totales=500 |isbn=978-2-7108-0715-5 |oclc=300489419 |id=Vidal}}.

== Voir aussi ==
{{autres projets
 | wikiversity = Thermodynamique des mélanges
 | wikiversity titre = Thermodynamique des mélanges
}}

=== Articles connexes ===

* [[Constante d'équilibre]]
* [[Équilibre chimique]]
* [[Fugacité]] 
* [[Loi de Henry]]
* [[Loi de Raoult]]
* [[Potentiel chimique]]
* [[Solution idéale]]
* Modèles de coefficient d'activité :
** [[Cosmospace|COSMOSPACE]]
** {{lien|lang=en|trad=Non-random two-liquid model|texte=NRTL}}
** [[Unifac|UNIFAC]]
** [[UNIQUAC]]

=== Liens externes ===
{{Liens}}

{{Portail|chimie|physique}}

[[Catégorie:Grandeur thermodynamique]]
[[Catégorie:Grandeur en thermochimie]]
[[Catégorie:Chimie générale]]
[[Catégorie:Grandeur sans dimension]]
[[Catégorie:Chimie des équilibres]]