Loi de Henry

En physique, et plus particulièrement en thermodynamique, la loi de Henry, établie empiriquement par le physicien britannique William Henry en 1803, énonce queModèle:, :

Modèle:Début citationÀ température constante et à saturation, la pression partielle dans la phase vapeur d'un soluté volatil est proportionnelle à la fraction molaire de ce corps dans la solution liquide.Modèle:Fin citation

En pratique, elle ne s'applique qu'aux faibles concentrations du soluté (fraction molaire inférieure à Modèle:Unité) et à des pressions de moins de Modèle:Unité (domaine d'application de la loi des gaz parfaits). Le soluté peut être un gaz dissout ou plus généralement tout corps volatil très faiblement soluble ou très dilué. Elle n'est également appliquable qu'à des mélanges binaires, ne contenant qu'un seul soluté et un seul solvant. Par extension à l'aide de coefficients de fugacité et d'activité, elle peut être appliquée à des mélanges multicomposants réels. Le pendant de la loi de Henry pour les solvants est la loi de Raoult.

Elle est utilisée dans de nombreux domaines de la chimie, de la physique et de la météorologie.

Énoncé, définitions et démonstration

Énoncé de la loi de Henry

On considère une solution liquide constituée d'un Modèle:Nobr dissout dans un Modèle:Nobr. La loi de Henry relie la pression partielle $P_{\sigma }$ du soluté en phase gazeuse à sa fraction molaire $x_{\sigma }^{\text{l}}$ en phase liquide à l'équilibre liquide-vapeur selonModèle:, :

Loi de Henry pression partielle du Modèle:Nobr dans le Modèle:Nobr : $P_{\sigma }=x_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}$

avec : Modèle:Début de colonnes

• $P$ la pression totale du mélange ;
• $P_{\sigma }$ la pression partielle du Modèle:Nobr, par définition $P_{\sigma }=x_{\sigma }^{\text{g}}P$ ;
• $k_{\text{H},\sigma ,s}$ la constante de Henry du Modèle:Nobr dans le Modèle:Nobr, aux pression $P$ et température $T$ du mélange ; la constante de Henry a la dimension d'une pression ; $k_{\text{H}}$ est la notation recommandée par le Green Book de l'Union internationale de chimie pure et appliquée (IUPAC), on trouve également $H$, $\mathscr{H}$ voire $K$ dans la littérature ;
• $x_{\sigma }^{\text{g}}$ la fraction molaire du Modèle:Nobr dans la phase vapeur ;
• $x_{\sigma }^{\text{l}}$ la fraction molaire du Modèle:Nobr dans la phase liquide.

La littérature utilise parfois l'inverse de la constante de Henry $k_{\text{H}}$ définie précédemment, $H=\frac{1}{k_{\text{H}}}$, et l'appelle également constante de Henry. Sa dimension est alors l'inverse de celle d'une pression, et la loi de Henry énonce queModèle:, :

Modèle:Début citationÀ température constante et à saturation, la quantité de gaz dissout dans un liquide est proportionnelle à la pression partielle qu'exerce ce gaz sur le liquide.Modèle:Fin citation

et s'écrit sous la forme :

Loi de Henry : $x_{\sigma }^{\text{l}}=P_{\sigma }{H}_{\sigma ,s}$

D'autres formes de la loi de Henry sont écrites en fonction de la concentration molaire ou de la molalité du soluté. Les lecteurs de littérature spécialisée doivent être attentifs à noter quelle version de l'équation de la loi de Henry est utilisée (voir le paragraphe Constantes de Henry pour divers gaz dissouts dans l'eau).

Un soluté répondant à la loi de Henry est une espèce chimique capable de passer en phase gaz dans les conditions de pression et de température considérées, c'est-à-dire un corps volatil. La loi de Henry ne s'applique donc pas aux solutés solides tels les sels. Dans les conditions de température $T$ et de pression $P$ de la solution, un soluté répondant à la loi de Henry est typiquement un fluide supercritique, pour lequel il n'existe pas de pression de vapeur saturante $P^{\text{sat}}$ à la température $T$, c'est-à-dire pour lequel $T>T_{\text{c}}$ sa température critique (par exemple à Modèle:Unité l'oxygène et l'azote, de $T_{\text{c}}$ respectives Modèle:Unité et Modèle:Unité), ou un fluide subcritique gazeux, pour lequel à la température $T$ il existe une pression de vapeur saturante telle que $P^{\text{sat}}>P$ (par exemple à Modèle:Unité et Modèle:Unité le propane de Modèle:Nobr). Le soluté peut également être un fluide subcritique liquide, pour lequel à la température $T$ il existe une pression de vapeur saturante telle que $P^{\text{sat}}<P$ (par exemple à Modèle:Unité et Modèle:Unité l'éthanol de Modèle:Nobr), présent en faible quantité dans la solution, soit parce qu'il est peu soluble, soit parce qu'il est fortement dilué. De façon générale, un Modèle:Nobr répondant à la loi de Henry est un corps volatil dans les conditions du mélange et dont la fraction molaire en phase liquide est faible, soit $x_{\sigma }^{\text{l}}\approx 0$. Le Modèle:Nobr est alors un corps se comportant quasiment comme un corps pur, soit $x_s^{\text{l}}\approx 1$. La relation de Duhem-Margules impose que si l'équilibre liquide-vapeur du soluté répond à la loi de Henry, alors celui du solvant répond à la loi de Raoult, et réciproquement.

L'équilibre liquide-vapeur déterminé par la loi de Henry est un état stable, appelé état de saturation du solvant par le soluté. Dans les conditions de pression et température données, le solvant peut contenir plus de soluté que la quantité déterminée par la loi de Henry, mais il s'agit alors d'un état d'équilibre instable dit de sursaturation. Dans ce cas la moindre perturbation (choc sur le récipient contenant le liquide, introduction d'une poussière formant un site de nucléation pour les bulles de gaz, fluctuation de pression ou de températureModèle:Etc.) peut provoquer le dégazage de l'excès de soluté dissout jusqu'à l'établissement de l'état stable dicté par la loi. De même, la quantité du soluté dissout peut être inférieure à celle déterminée par la loi de Henry : il y a sous-saturation. Dans ce cas, si le soluté est présent en phase gaz, la phase liquide absorbe du soluté gazeux jusqu'à atteindre l'équilibre stable. La fraction $x_{\sigma }^{\text{l}}$ déterminée par la loi de Henry est donc la fraction molaire maximale de soluté que peut contenir la phase liquide de façon stable : la fraction $x_{\sigma }^{\text{l}}$ est la solubilité du Modèle:Nobr dans le Modèle:Nobr dans les conditions de pression et de température données.

Constante de Henry

Les définitions et formules suivantes ne sont valables que pour un mélange binaire comprenant un unique Modèle:Nobr et un unique Modèle:Nobr.

La constante de Henry est définie rigoureusement en thermodynamique à partir de la fugacité.

Contrairement à ce que peut laisser entendre le terme de constante, la constante de Henry dépend de la pression et de la température. En revanche, elle ne dépend pas de la composition du mélange. La constante de Henry $k_{\text{H},\sigma ,s}$ dépend également de la nature du Modèle:Nobr et du Modèle:Nobr ; ceci implique qu'elle doit être déterminée pour chaque couple « Modèle:Nobr - Modèle:Nobr » et n'est pas valable si l'un de ces deux corps est considéré dans un mélange binaire autre que celui pour lequel elle a été déterminée (par exemple le Modèle:Nobr avec un solvant autre que le Modèle:Nobr).

En pratique, la constante de Henry est déterminée expérimentalement.

Définition thermodynamique

Évolution de la fugacité en fonction de la fraction molaire à pression et température constantesModèle:,.

En thermodynamique, à pression et température constantes, la fugacité $f_{\sigma }^{\text{l}}$ d'une espèce chimique $\sigma$ (soluté) en phase liquide, en présence d'une deuxième espèce $s$ (solvant), possède deux limites, avec $x_{\sigma }^{\text{l}}$ la fraction molaire du Modèle:Nobr dans le mélange :

• à dilution infinie : $\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{f}_{\sigma }^{\text{l}}=0$ ;
• pour le corps pur : $\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 1}{x_s^{\text{l}}\to 0}}{\lim }{f}_{\sigma }^{\text{l}}=f_{\sigma }^{\text{l},*}$ ;

avec $f_{\sigma }^{\text{l},*}$ la fugacité du Modèle:Nobr à l'état de liquide pur. Cette fugacité $f_{\sigma }^{\text{l},*}$ peut être fictive si le Modèle:Nobr est gazeux à l'état pur dans les conditions de pression et température données.

L'évolution de la fugacité en fonction de la composition est encadrée par deux lois linéairesModèle:, :

Loi de Henry : $f_{\sigma }^{\text{l}}\approx x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}$ aux faibles concentrations. Loi de Lewis et Randall : $f_{\sigma }^{\text{l}}\approx x_{\sigma }^{\text{l}}{f}_{\sigma }^{\text{l},*}$ aux fortes concentrations.

La constante de Henry $k_{\text{H},\sigma ,s}$ n'est pas la fugacité du Modèle:Nobr à dilution infinie dans le Modèle:Nobr. La fugacité $f_{\sigma }^{\text{l}}$ tend vers zéro lorsque $x_{\sigma }^{\text{l}}$ tend vers zéro. Aussi la constante de Henry est-elle définie comme étant la limite lorsque la quantité de Modèle:Nobr dissout en phase liquide s'annuleModèle:,Modèle:,Modèle:, :

Constante de Henry : $k_{\text{H},\sigma ,s}\equiv \underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }\frac{f_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l}}}$ à pression et température constantes.

avec : Modèle:Début de colonnes

• $f_{\sigma }^{\text{l}}$ la fugacité du Modèle:Nobr dans le mélange liquide ;
• $k_{\text{H},\sigma ,s}$ la constante de Henry du Modèle:Nobr dans le Modèle:Nobr, aux pression $P$ et température $T$ du mélange ;
• $x_{\sigma }^{\text{l}}$ la fraction molaire du Modèle:Nobr dans le mélange liquide ;
• $x_s^{\text{l}}$ la fraction molaire du Modèle:Nobr dans le mélange liquide ($x_{\sigma }^{\text{l}}+x_s^{\text{l}}=1$).

En application de la règle de L'Hôpital, la constante de Henry peut également être définie parModèle:,Modèle:, :

Constante de Henry : $k_{\text{H},\sigma ,s}={\left(\frac{\partial f_{\sigma }^{\text{l}}}{\partial x_{\sigma }^{\text{l}}}\right)}_{P,T,x_{\sigma }^{\text{l}}=0}$

La constante de Henry est donc la pente de la fugacité à dilution infinie.

Quelle que soit la concentration $0\le x_{\sigma }^{\text{l}}\le 1$ du Modèle:Nobr, sa fugacité réelle $f_{\sigma }^{\text{l}}$ peut être exprimée en fonction d'un coefficient d'activité à partir des deux lois linéaires idéales définies précédemmentModèle:, :

$f_{\sigma }^{\text{l}}=x_{\sigma }^{\text{l}}{\gamma }_{\text{H},\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}=x_{\sigma }^{\text{l}}{\gamma }_{\sigma }^{\text{l}}{f}_{\sigma }^{\text{l},*}$

en posant :

• ${\gamma }_{\text{H},\sigma }^{\text{l}}$ le coefficient d'activité défini par rapport à la loi de Henry ;
• ${\gamma }_{\sigma }^{\text{l}}$ le coefficient d'activité défini par rapport à la loi de Lewis et Randall.

Puisque les deux limites de la fugacité $f_{\sigma }^{\text{l}}$ sont définies, on a les limites des coefficients d'activitéModèle:, :

• à dilution infinie : $\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\gamma }_{\text{H},\sigma }^{\text{l}}=1$ ;
• pour le corps pur : $\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 1}{x_s^{\text{l}}\to 0}}{\lim }{\gamma }_{\sigma }^{\text{l}}=1$.

On pose à dilution infinieModèle:, :

Coefficient d'activité à dilution infinie : ${\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}\equiv \underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\gamma }_{\sigma }^{\text{l}}$

Par conséquent, à dilution infinie on a :

$\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }\frac{f_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l}}}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }\left({\gamma }_{\text{H},\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}\right)=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }\left({\gamma }_{\sigma }^{\text{l}}{f}_{\sigma }^{\text{l},*}\right)$

d'où la relationModèle:, (puisque $k_{\text{H},\sigma ,s}$ et $f_{\sigma }^{\text{l},*}$ ne dépendent pas de la composition) :

Constante de Henry : $k_{\text{H},\sigma ,s}={\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}{f}_{\sigma }^{\text{l},*}$

En injectant la relation précédente dans les expressions de $f_{\sigma }^{\text{l}}$ on obtientModèle:, :

${\gamma }_{\text{H},\sigma }^{\text{l}}=\frac{{\gamma }_{\sigma }^{\text{l}}}{{\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}$

La première relation permet de déterminer $k_{\text{H},\sigma ,s}$ si l'on connait $f_{\sigma }^{\text{l},*}$. Inversement, connaissant $k_{\text{H},\sigma ,s}$ on peut extrapoler $f_{\sigma }^{\text{l},*}$ si le Modèle:Nobr n'existe pas à l'état de liquide pur dans les conditions de pression et température données. La deuxième relation montre que les deux coefficients d'activité ${\gamma }_{\text{H},\sigma }^{\text{l}}$ et ${\gamma }_{\sigma }^{\text{l}}$ ne sont pas indépendants, bien que liés à des états de référence différents. La loi de Henry, quelle que soit sa forme, peut ainsi être employée avec les modèles classiques développés pour la loi de Lewis et Randall (Margules, Modèle:Lien, Wilson, Modèle:Lien, UNIQUAC, UNIFAC, COSMOSPACEModèle:Etc).

Dépendance à la pression

La fugacité $f_{\sigma }^{\text{l}}$ du Modèle:Nobr dans le mélange liquide varie en fonction de la pression selon :

${\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n^{\text{l}}}=\frac{{\overline{V}}_{\sigma }^{\text{l}}}{RT}$

avec : Modèle:Début de colonnes

• $V^{\text{l}}$ le volume de la phase liquide ;
• ${\overline{V}}_{\sigma }^{\text{l}}={\left(\frac{\partial V^{\text{l}}}{\partial n_{\sigma }^{\text{l}}}\right)}_{P,T,n_s^{\text{l}}}$ le volume molaire partiel du Modèle:Nobr dans le mélange liquide ;
• $n_{\sigma }^{\text{l}}$ la quantité du Modèle:Nobr dans le mélange liquide ;
• $n_s^{\text{l}}$ la quantité du Modèle:Nobr dans le mélange liquide.

Quelle que soit la fraction molaire $x_{\sigma }^{\text{l}}=n_{\sigma }^{\text{l}}/(n_{\sigma }^{\text{l}}+n_s^{\text{l}})$ du Modèle:Nobr, la dérivée partielle étant effectuée à composition constante, on peut écrire :

${\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n^{\text{l}}}={\left(\frac{\partial [\ln f_{\sigma }^{\text{l}}-\ln x_{\sigma }^{\text{l}}+\ln x_{\sigma }^{\text{l}}]}{\partial P}\right)}_{T,n^{\text{l}}}={\left(\frac{\partial \ln \frac{f_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l}}}}{\partial P}\right)}_{T,n^{\text{l}}}+\underset{=0\text{à composition constante}}{\underbrace{{\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n^{\text{l}}}}}$

En passant à la limite de la dilution infinie :

$\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n^{\text{l}}}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\left(\frac{\partial \ln \frac{f_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l}}}}{\partial P}\right)}_{T,n^{\text{l}}}={\left(\frac{\partial \ln k_{\text{H},\sigma ,s}}{\partial P}\right)}_T$

La référence à la composition constante disparait dans la dérivée partielle de la constante de Henry, puisque celle-ci ne dépend pas de la composition. On pose pour le volume molaire partiel :

Volume molaire partiel du soluté à dilution infinie : ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\equiv \underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\overline{V}}_{\sigma }^{\text{l}}$

La constante de Henry dépend par conséquent de la pression selonModèle:, :

Dépendance de la constante de Henry à la pression ${\left(\frac{\partial \ln k_{\text{H},\sigma ,s}}{\partial P}\right)}_T=\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}$

avec : Modèle:Début de colonnes

• $P$ la pression ;
• $T$ la température ;
• ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ le volume molaire partiel du Modèle:Nobr à dilution infinie dans le Modèle:Nobr ;
• $R$ la constante universelle des gaz parfaits.

En intégrant cette relation entre une pression de référence $P^{\circ }$ et la pression $P$ :

$\ln k_{\text{H},\sigma ,s}(P,T)-\ln k_{\text{H},\sigma ,s}(P^{\circ },T)=\frac{\underset{P^{\circ }}{\overset{P}{\int }}{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\mathrm{d}P}{RT}$

La pression de référence $P^{\circ }$ est le plus souvent prise égale à la pression de vapeur saturante du Modèle:Nobr à la température $T$ du mélange : $P^{\circ }=P_s^{\text{sat}}\left(T\right)$. En conséquence, on peut réduire la constante d'intégration à une fonction de la température seule : $k_{\text{H},\sigma ,s}(P^{\circ },T)=k_{\text{H},\sigma ,s}(P_s^{\text{sat}}\left(T\right),T)=k_{\text{H},\sigma ,s}^{\text{sat}}\left(T\right)$. La constante de Henry est alors exprimée sous la formeModèle:,Modèle:, :

$k_{\text{H},\sigma ,s}(P,T)=k_{\text{H},\sigma ,s}^{\text{sat}}\left(T\right){𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$

avec le facteur de PoyntingModèle:, :

Facteur de Poynting : ${𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\equiv \exp \left(\frac{\underset{P_s^{\text{sat}}}{\overset{P}{\int }}{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\mathrm{d}P}{RT}\right)$

Le volume molaire partiel ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ représente la variation de volume de la solution liquide due à la dissolution d'une mole de Modèle:Nobr dans une quantité infinie de Modèle:Nobr. Il peut être déterminé expérimentalement par extrapolation de ${\overline{V}}_{\sigma }^{\text{l}}$ établi pour plusieurs concentrations de soluté dans le mélange liquide ; il existe également des corrélations telles que celle de Brelvi-O'Connell. Les liquides étant peu compressibles, le volume molaire partiel ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ peut être considéré comme ne dépendant pas de la pression, soit ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\approx {\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\left(T\right)$, on obtient :

$\ln k_{\text{H},\sigma ,s}=\ln k_{\text{H},\sigma ,s}^{\text{sat}}+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}$

Il peut être aussi bien positif (la dissolution du gaz provoque une dilatation du liquide) que négatif (la dissolution du gaz provoque une contraction du liquide). Si le volume molaire partiel ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ est positif alors la constante de Henry $k_{\text{H},\sigma ,s}$ augmente avec la pression $P$.

Dépendance à la température

La fugacité $f_{\sigma }^{\text{l}}$ du Modèle:Nobr dans le mélange liquide varie en fonction de la température selon :

${\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial \frac{1}{T}}\right)}_{P,n^{\text{l}}}=\frac{{\overline{H}}_{\sigma }^{\text{l}}-{\overline{H}}_{\sigma }^{\bullet ,*}}{R}$

avec : Modèle:Début de colonnes

• $H^{\text{l}}$ l'enthalpie de la phase liquide ;
• ${\overline{H}}_{\sigma }^{\text{l}}={\left(\frac{\partial H^{\text{l}}}{\partial n_{\sigma }^{\text{l}}}\right)}_{P,T,n_s^{\text{l}}}$ l'enthalpie molaire partielle du Modèle:Nobr dans le mélange liquide ;
• ${\overline{H}}_{\sigma }^{\bullet ,*}$ l'enthalpie molaire du Modèle:Nobr à l'état de gaz parfait pur à $T$ ;
• $n_{\sigma }^{\text{l}}$ la quantité du Modèle:Nobr dans le mélange liquide ;
• $n_s^{\text{l}}$ la quantité du Modèle:Nobr dans le mélange liquide.

Quelle que soit la fraction molaire $x_{\sigma }^{\text{l}}=n_{\sigma }^{\text{l}}/(n_{\sigma }^{\text{l}}+n_s^{\text{l}})$ du Modèle:Nobr, la dérivée partielle étant effectuée à composition constante, on peut écrire :

${\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial \frac{1}{T}}\right)}_{P,n^{\text{l}}}={\left(\frac{\partial [\ln f_{\sigma }^{\text{l}}-\ln x_{\sigma }^{\text{l}}+\ln x_{\sigma }^{\text{l}}]}{\partial \frac{1}{T}}\right)}_{P,n^{\text{l}}}={\left(\frac{\partial \ln \frac{f_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l}}}}{\partial \frac{1}{T}}\right)}_{P,n^{\text{l}}}+\underset{=0\text{à composition constante}}{\underbrace{{\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial \frac{1}{T}}\right)}_{P,n^{\text{l}}}}}$

En passant à la limite de la dilution infinie :

$\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial \frac{1}{T}}\right)}_{P,n^{\text{l}}}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\left(\frac{\partial \ln \frac{f_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l}}}}{\partial \frac{1}{T}}\right)}_{P,n^{\text{l}}}={\left(\frac{\partial \ln k_{\text{H},\sigma ,s}}{\partial \frac{1}{T}}\right)}_P$

La référence à la composition constante disparait dans la dérivée partielle de la constante de Henry, puisque celle-ci ne dépend pas de la composition. On pose pour l'enthalpie molaire partielleModèle:, :

Enthalpie molaire partielle à dilution infinie du Modèle:Nobr : ${\overline{H}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\equiv \underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\overline{H}}_{\sigma }^{\text{l}}$

c'est-à-dire l'enthalpie molaire partielle du Modèle:Nobr à dilution infinie dans le Modèle:Nobr liquide. L'enthalpie molaire du gaz parfait pur ${\overline{H}}_{\sigma }^{\bullet ,*}$ ne dépendant pas de la composition, elle reste inchangée lors du passage à la limite et on obtient :

$\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }({\overline{H}}_{\sigma }^{\text{l}}-{\overline{H}}_{\sigma }^{\bullet ,*})={\overline{H}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}-{\overline{H}}_{\sigma }^{\bullet ,*}$

On pose pour les enthalpies molaires partiellesModèle:, :

Enthalpie de dissolution : ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}\equiv {\overline{H}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}-{\overline{H}}_{\sigma }^{\bullet ,*}$

avec :

• ${\overline{H}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ l'enthalpie molaire partielle du Modèle:Nobr à dilution infinie dans le Modèle:Nobr liquide à $T$ ;
• ${\overline{H}}_{\sigma }^{\bullet ,*}$ l'enthalpie molaire du Modèle:Nobr à l'état de gaz parfait pur à $T$.

La constante de Henry dépend par conséquent de la température selonModèle:, :

Dépendance de la constante de Henry à la température ${\left(\frac{\partial \ln k_{\text{H},\sigma ,s}}{\partial \frac{1}{T}}\right)}_P=\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}$

avec : Modèle:Début de colonnes

• $P$ la pression ;
• $T$ la température ;
• ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}$ l'enthalpie de dissolution du Modèle:Nobr dans le Modèle:Nobr à $T$ ;
• $R$ la constante universelle des gaz parfaits.

Si l'on considère l'enthalpie de dissolution ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}$ comme constante, alors, en intégrant cette relation entre une température de référence $T^{\circ }$ et la température $T$ :

$\ln k_{\text{H},\sigma ,s}(P,T)-\ln k_{\text{H},\sigma ,s}(P,T^{\circ })=\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}(\frac{1}{T}-\frac{1}{T^{\circ }})$

$k_{\text{H},\sigma ,s}(P,T)=k_{\text{H},\sigma ,s}(P,T^{\circ })\cdot \exp \left(\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}(\frac{1}{T}-\frac{1}{T^{\circ }})\right)$

Cette forme n'est applicable que sur des plages de température relativement étroites. Elle est généralisée au moyen de deux constantes $A$ et $B$ empiriques spécifiques du couple « Modèle:Nobr - Modèle:Nobr » :

$\ln k_{\text{H},\sigma ,s}(P,T)=A+\frac{B}{T}$

La littérature utilise parfois l'inverse de la constante de Henry définie précédemment, $H_{\sigma ,s}=\frac{1}{k_{\text{H},\sigma ,s}}$ (cette notation prête à confusion avec celle de l'enthalpie de dissolution ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}$), aussi trouve-t-on également les relations :

${\left(\frac{\partial \ln H_{\sigma ,s}}{\partial \frac{1}{T}}\right)}_P=-\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}$

$\ln H_{\sigma ,s}(P,T)-\ln H_{\sigma ,s}(P,T^{\circ })=-\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}(\frac{1}{T}-\frac{1}{T^{\circ }})$

$H_{\sigma ,s}(P,T)=H_{\sigma ,s}(P,T^{\circ })\cdot \exp (-\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}(\frac{1}{T}-\frac{1}{T^{\circ }}))$

L'enthalpie de dissolution ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}$ est la chaleur produite par la dissolution d'une mole de Modèle:Nobr à l'état de gaz parfait pur dans une quantité infinie de Modèle:Nobr à l'état liquide. Elle est déterminée expérimentalement par calorimétrie en extrapolant la chaleur de dissolution d'une mole de soluté dans plusieurs quantités de solvant. On peut considérer l'enthalpie molaire partielle d'un corps dans un mélange liquide, ici ${\overline{H}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$, comme indépendante de la pression, les liquides étant peu compressibles. De même, en vertu de la deuxième loi de Joule, l'enthalpie molaire d'un gaz parfait, ici ${\overline{H}}_{\sigma }^{\bullet ,*}$, ne dépend pas de la pression. Ainsi, il peut être considéré que l'enthalpie de dissolution ne dépend que de la température : ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}={\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}\left(T\right)$. Elle peut être négative (dissolution exothermique, l'opération de dissolution dégage de la chaleur), positive (dissolution endothermique, l'opération de dissolution absorbe de la chaleur) ou nulle (dissolution athermique). Pour la plupart des gaz à température ambiante la dissolution est exothermique, soit ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}<0$, par conséquent $\ln k_{\text{H},\sigma ,s}$ augmente avec une diminution de $\frac{1}{T}$ et la constante de Henry $k_{\text{H},\sigma ,s}$ augmente avec la température $T$.

Calcul par une équation d'état

Le coefficient de fugacité ${\varphi }_{\sigma }^{\text{l}}$ du Modèle:Nobr en phase liquide est défini par la relation avec la fugacité $f_{\sigma }^{\text{l}}$ :

Coefficient de fugacité : ${\varphi }_{\sigma }^{\text{l}}\equiv \frac{f_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l}}P}$

Par définition de la constante de Henry, on a doncModèle:, :

$\frac{k_{\text{H},\sigma ,s}}{P}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }{\varphi }_{\sigma }^{\text{l}}$ à pression et température constantes

Si l'on dispose d'une équation d'état d'une phase liquide donnant la pression $P$ en fonction du volume $V^{\text{l}}$, de la température $T$ et de la composition $n^{\text{l}}$, soit $P=P(V^{\text{l}},T,n^{\text{l}})$, le coefficient de fugacité ${\varphi }_{\sigma }^{\text{l}}$ se calcule selon :

$RT\ln {\varphi }_{\sigma }^{\text{l}}=-\underset{+\mathrm{\infty }}{\overset{V^{\text{l}}}{\int }}[{\left(\frac{\partial P}{\partial n_{\sigma }^{\text{l}}}\right)}_{V,T,n_s^{\text{l}}}-\frac{RT}{V}]\mathrm{d}V-RT\ln \left(\frac{P{V}^{\text{l}}}{(n_{\sigma }^{\text{l}}+n_s^{\text{l}})RT}\right)$

avec $n_{\sigma }^{\text{l}}$ et $n_s^{\text{l}}$ les quantités respectives du Modèle:Nobr et du Modèle:Nobr dans la solution liquide. Il est donc possible de calculer la constante de Henry $k_{\text{H},\sigma ,s}$ à partir d'une équation d'état de la phase liquide. Toutefois, les équations d'état telles que les équations d'état cubiques sont généralement développées pour représenter des phases gazeuses et représentent assez mal les phases liquides. Cette démarche reste donc théorique ; en pratique la constante de Henry est plutôt déterminée expérimentalement sous des formes empiriques présentées au paragraphe Formes usuelles. La relation établie ci-dessus et l'exemple ci-dessous montrent cependant la dépendance de la constante de Henry $k_{\text{H},\sigma ,s}$ aux propriétés du Modèle:Nobr et du Modèle:Nobr, et aux interactions entre les deux constituants.

L'équation d'état de van der Waals donne, à pression $P$ et température $T$, pour un mélange binaire liquide de composition $x_{\sigma }$ et $x_s$ :

$P=\frac{RT}{{\overline{V}}^{\text{l}}-b_m^{\text{l}}}-\frac{a_m^{\text{l}}}{{{\overline{V}}^{\text{l}}}^2}$

$RT\ln {\varphi }_{\sigma }^{\text{l}}=-\frac{a_m^{\text{l}}}{{\overline{V}}^{\text{l}}}[{\delta }_{\sigma }^{\text{l}}-\frac{b_{\sigma }}{b_m^{\text{l}}}]+\frac{b_{\sigma }}{b_m^{\text{l}}}(P{\overline{V}}^{\text{l}}-RT)-RT\ln \left(\frac{P\cdot ({\overline{V}}^{\text{l}}-b_m^{\text{l}})}{RT}\right)$

avec :

$a_m^{\text{l}}={x_{\sigma }^{\text{l}}}^2{a}_{\sigma }+2x_{\sigma }^{\text{l}}{x}_s^{\text{l}}\sqrt{a_{\sigma }{a}_s}(1-k_{\sigma ,s})+{x_s^{\text{l}}}^2{a}_s$

$b_m^{\text{l}}=x_{\sigma }^{\text{l}}{b}_{\sigma }+x_s^{\text{l}}{b}_s$

${\delta }_{\sigma }^{\text{l}}=2\frac{\sqrt{a_{\sigma }}}{a_m^{\text{l}}}(x_{\sigma }^{\text{l}}\sqrt{a_{\sigma }}+x_s^{\text{l}}\sqrt{a_s}(1-k_{\sigma ,s}))$

À dilution infinie ($x_{\sigma }^{\text{l}}\to 0$ et $x_s^{\text{l}}\to 1$ à $P$ et $T$ constantes) on a :

$a_m^{\text{l}}\to a_s$

$b_m^{\text{l}}\to b_s$

${\delta }_{\sigma }^{\text{l}}\to {\delta }_{\sigma ,s}^{\mathrm{\infty }}=2\sqrt{\frac{a_{\sigma }}{a_s}}(1-k_{\sigma ,s})$

Le volume molaire ${\overline{V}}^{\text{l}}$ de la solution liquide tend vers celui ${\overline{V}}_s^{\text{l,*}}$ du Modèle:Nobr liquide pur à $P$ et $T$, calculable par l'équation d'état : ${\overline{V}}^{\text{l}}\to {\overline{V}}_s^{\text{l,*}}$. On obtient $k_{\text{H},\sigma ,s}$, la constante de Henry à $P$ et $T$ :

$RT\ln \frac{k_{\text{H},\sigma ,s}}{P}=-\frac{a_s}{{\overline{V}}_s^{\text{l,*}}}[{\delta }_{\sigma ,s}^{\mathrm{\infty }}-\frac{b_{\sigma }}{b_s}]+\frac{b_{\sigma }}{b_s}(P{\overline{V}}_s^{\text{l,*}}-RT)-RT\ln \left(\frac{P\cdot ({\overline{V}}_s^{\text{l,*}}-b_s)}{RT}\right)$

Lorsque la pression $P$ tend vers la pression de vapeur saturante $P_s^{\text{sat}}$ du Modèle:Nobr à la température $T$, l'équilibre liquide-vapeur tend vers celui du Modèle:Nobr pur (état de saturation) ; le volume molaire ${\overline{V}}_s^{\text{l,*}}$ tend vers celui ${\overline{V}}_s^{\text{l,*,sat}}$ du Modèle:Nobr liquide pur à saturation à $T$. On obtient $k_{\text{H},\sigma ,s}^{\text{sat}}$, la constante de Henry à $P_s^{\text{sat}}$ et $T$, c'est-à-dire dans les conditions de saturation du Modèle:Nobr :

$RT\ln \frac{k_{\text{H},\sigma ,s}^{\text{sat}}}{P_s^{\text{sat}}}=-\frac{a_s}{{\overline{V}}_s^{\text{l,*,sat}}}[{\delta }_{\sigma ,s}^{\mathrm{\infty }}-\frac{b_{\sigma }}{b_s}]+\frac{b_{\sigma }}{b_s}(P_s^{\text{sat}}{\overline{V}}_s^{\text{l,*,sat}}-RT)-RT\ln \left(\frac{P_s^{\text{sat}}\cdot ({\overline{V}}_s^{\text{l,*,sat}}-b_s)}{RT}\right)$

Les liquides étant peu compressibles, le volume molaire ${\overline{V}}_s^{\text{l,*}}$ du Modèle:Nobr liquide pur à $P$ est assimilable à celui ${\overline{V}}_s^{\text{l,*,sat}}$ du Modèle:Nobr liquide pur à saturation : ${\overline{V}}_s^{\text{l,*}}\approx {\overline{V}}_s^{\text{l,*,sat}}$. En considérant la relation avec le facteur de Poynting :

$\ln k_{\text{H},\sigma ,s}\approx \ln k_{\text{H},\sigma ,s}^{\text{sat}}+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}$

on identifie ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$, le volume molaire partiel du Modèle:Nobr à dilution infinie dans le Modèle:Nobr :

${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\approx \frac{b_{\sigma }}{b_s}{\overline{V}}_s^{\text{l,*,sat}}$

Avec : Modèle:Début de colonnes

• $a_{\sigma }=\frac{27R^2{T_{\text{c},\sigma }}^2}{64P_{\text{c},\sigma }}$ et $a_s=\frac{27R^2{T_{\text{c},s}}^2}{64P_{\text{c},s}}$ ;
• $b_{\sigma }=\frac{R{T}_{\text{c},\sigma }}{8P_{\text{c},\sigma }}$ et $b_s=\frac{R{T}_{\text{c},s}}{8P_{\text{c},s}}$ ;
• $k_{\sigma ,s}$ un coefficient d'interaction binaire ;
• $P_{\text{c},\sigma }$ et $P_{\text{c},s}$ les pressions critiques respectives du Modèle:Nobr et du Modèle:Nobr ;
• $P_s^{\text{sat}}$ la pression de vapeur saturante du Modèle:Nobr à la température $T$ ;
• $T_{\text{c},\sigma }$ et $T_{\text{c},s}$ les températures critiques respectives du Modèle:Nobr et du Modèle:Nobr ;
• ${\overline{V}}^{\text{l}}$ le volume molaire de la solution liquide à la pression $P$ et à la température $T$ ;
• ${\overline{V}}_s^{\text{l,*}}$ le volume molaire du Modèle:Nobr liquide pur à $P$ et $T$ ;
• ${\overline{V}}_s^{\text{l,*,sat}}$ le volume molaire du Modèle:Nobr liquide pur à saturation, soit à $P_s^{\text{sat}}$ et $T$.

Cet exemple montre la dépendance de la constante de Henry $k_{\text{H},\sigma ,s}$ aux propriétés du Modèle:Nobr et du Modèle:Nobr, et aux interactions entre le Modèle:Nobr et le Modèle:Nobr par l'intermédiaire du coefficient d'interaction binaire $k_{\sigma ,s}$. La constante de Henry est donc spécifique du couple « Modèle:Nobr - Modèle:Nobr » et ne peut être utilisée pour d'autres mélanges binaires ; par exemple, elle n'est pas valable si le Modèle:Nobr est dissout dans un solvant autre que $s$.

Formes usuelles

La constante de Henry est souvent utilisée sous la forme obtenue par intégration par rapport à la pression :

Forme usuelle de la constante de Henry $k_{\text{H},\sigma ,s}(P,T)=k_{\text{H},\sigma ,s}^{\text{sat}}\left(T\right)\cdot \exp \left(\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}\right)$

La correction de Poynting ne devient significative qu'aux hautes pressions. Pour des pressions de l'ordre de grandeur de la pression atmosphérique, le facteur de Poynting est négligeable : ${𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\approx 1$. La constante de Henry peut alors être considérée comme indépendante de la pression et être approchée par :

Aux basses pressions : $k_{\text{H},\sigma ,s}(P,T)\approx k_{\text{H},\sigma ,s}^{\text{sat}}\left(T\right)$

Les formes suivantes sont souvent utilisées pour la dépendance à la températureModèle:, :

$\ln k_{\text{H},\sigma ,s}^{\text{sat}}\left(T\right)=A+\frac{B}{T}+C\ln T+DT$

$\ln k_{\text{H},\sigma ,s}^{\text{sat}}\left(T\right)=A+\frac{B}{T}+\frac{C}{T^2}+\frac{D}{T^3}$

avec $A$, $B$, $C$ et $D$ des constantes empiriques spécifiques du couple « Modèle:Nobr - Modèle:Nobr ». L'enthalpie de dissolution est alors exprimée sous les formes respectives :

$\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}=B-CT-D{T}^2$

$\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}=B+\frac{2C}{T}+\frac{3D}{T^2}$

Démonstration de la loi de Henry

Lorsque l'équilibre liquide-vapeur est atteint, les fugacités du Modèle:Nobr sont homogènes entre les deux phases :

$f_{\sigma }^{\text{g}}=f_{\sigma }^{\text{l}}$

avec :

• $f_{\sigma }^{\text{g}}$ la fugacité du Modèle:Nobr en phase gaz (vapeur) ;
• $f_{\sigma }^{\text{l}}$ la fugacité du Modèle:Nobr en phase liquide.

Aux basses pressions (moins de Modèle:Unité), le gaz se comporte comme un mélange de gaz parfaits, et la fugacité du Modèle:Nobr en phase gaz peut être assimilée à sa pression partielle :

$f_{\sigma }^{\text{g}}\approx x_{\sigma }^{\text{g}}P$

D'autre part, par définition, aux faibles concentrations la fugacité du Modèle:Nobr dans le Modèle:Nobr en phase liquide suit approximativement la loi linéaire :

$f_{\sigma }^{\text{l}}\approx x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}$

Ainsi, aux basses pressions et aux faibles concentrations, l'équilibre liquide-vapeur du Modèle:Nobr est approché par la relation :

$x_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}$

qui est la loi de Henry. Aux fortes concentrations, la fugacité en phase liquide suit approximativement la loi de Lewis et Randall : aux basses pressions et aux fortes concentrations ceci conduit à la loi de Raoult qui s'applique aux solvants.

Limites et extensions de la loi de Henry

Limites de la loi idéale

Conditions de pression

La loi de Henry n'est valable que si la phase gaz peut être considérée comme un mélange de gaz parfaits. Autrement dit, elle ne s'applique qu'à des pressions partielles de soluté de l'ordre de la pression atmosphérique (moins de Modèle:Unité), dans le domaine d'application de la loi des gaz parfaits.

Composition de la phase liquide

Concentration du soluté

La loi de Henry est une loi limitante qui ne s'applique qu'aux solutions suffisamment diluées ayant un comportement idéal. Dans ces solutions, le soluté et le solvant ont des comportements similaires. Typiquement, la loi de Henry ne s'applique que si la fraction molaire $x_{\sigma }^{\text{l}}$ du soluté est inférieure à 0,03 ou 0,05.

Présence d'autres solutés

La loi de Henry est établie pour un soluté unique dissout dans un solvant unique. Si le solvant contient plusieurs solutés la constante de Henry est modifiée et dépend de la composition. Ainsi la solubilité d'un gaz dans l'eau de mer est-elle inférieure à celle dans l'eau douce en raison de la compétition entre le gaz dissout et les sels dissouts. La constante de Henry pourra être corrigée selon l'équation empirique de SetchenovModèle:,Modèle:, :

Équation de Setchenov : $\ln \frac{k_{\text{H},\sigma ,s}}{k_{\text{H},\sigma ,s}^{*}}=kI$

avec : Modèle:Début de colonnes

• $k_{\text{H},\sigma ,s}$ la constante de Henry du Modèle:Nobr dans le Modèle:Nobr en solution avec tous les solutés ;
• $k_{\text{H},\sigma ,s}^{*}$ la constante de Henry du Modèle:Nobr dans le Modèle:Nobr en solution avec $\sigma$ seul soluté ;
• $k$ le coefficient de Setchenov, qui dépend des solutés et du solvant ;
• $I$ la force ionique ; on trouve aussi cette équation exprimée en fonction de la molalité des sels dissouts.

Dissolution réactive

Un soluté peut réagir avec le solvant dans lequel il est dissout, produisant d'autres espèces chimiques. La loi de Henry ne s'applique alors qu'à la réaction de dissolution du soluté dans le solvant et non aux équilibres chimiques impliquant cette espèce dissoute et les produits de réaction. Si l'on tient compte des produits de réaction, alors la solubilité d'un soluté réactif peut être très différente de ce que prédit la loi de Henry en l'absence de réaction. En solution aqueuse, le dioxyde de carbone Modèle:Fchim, le dioxyde de soufre Modèle:Fchim et l'ammoniac Modèle:Fchim sont des exemples de solutés réagissant avec le solvant (l'eau) et produisant des ions en fonction du pH.

Un exemple usuel de gaz réagissant avec le solvant est le dioxyde de carbone (Modèle:Fchim), qui forme partiellement, par réaction avec l'eau, de l'acide carbonique (Modèle:Fchim), qui lui-même, en fonction du pH de l'eau, forme les ions hydrogénocarbonate (Modèle:Fchim) et carbonate (Modèle:Fchim) selon les réactions :

• $\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{(}\mathrm{g}\mathrm{)}\mathrm{\rightleftarrows }\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}$ (dissolution)

  $\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{]}=K_1{P}_{\mathrm{C}{\mathrm{O}}_{\mathrm{2}}}$ (loi de Henry)

• $\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}\mathrm{+}{\mathrm{H}}_{\mathrm{2}}\mathrm{O}\mathrm{\rightleftarrows }{\mathrm{H}}_{\mathrm{2}}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}$

  $\mathrm{[}{\mathrm{H}}_{\mathrm{2}}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}\mathrm{]}=K_2\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{]}$

• ${\mathrm{H}}_{\mathrm{2}}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}\mathrm{+}{\mathrm{H}}_{\mathrm{2}}\mathrm{O}\mathrm{\rightleftarrows }{\mathrm{H}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{-}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}\mathrm{+}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}$

  $\mathrm{[}{\mathrm{H}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{-}}\mathrm{]}\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}=K_3\mathrm{[}{\mathrm{H}}_{\mathrm{2}}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}\mathrm{]}$

• ${\mathrm{H}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{-}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}\mathrm{+}{\mathrm{H}}_{\mathrm{2}}\mathrm{O}\mathrm{\rightleftarrows }{\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{2}\mathrm{-}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}\mathrm{+}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}$

  $\mathrm{[}{\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{2}\mathrm{-}}\mathrm{]}\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}=K_4\mathrm{[}{\mathrm{H}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{-}}\mathrm{]}$

avec :

• $K_i$ la constante de l'équilibre $i$ ; à Modèle:Unité : Modèle:Nobr, Modèle:Nobr, Modèle:Nobr, Modèle:Nobr ;
• $P_{\mathrm{C}{\mathrm{O}}_{\mathrm{2}}}$ la pression partielle du Modèle:Fchim en phase gaz (en Modèle:Unité) ;
• $\mathrm{[}\mathrm{X}\mathrm{]}$ la concentration molaire de l'espèce $\mathrm{X}$ en solution aqueuse (en Modèle:Unité).

La loi de Henry à proprement parler ne s'applique qu'à la première des quatre réactions ci-dessus, celle de dissolution du Modèle:Fchim gazeux en Modèle:Fchim aqueux. Si l'on considère le Modèle:Fchim dissout sous toutes ses formes, le bilan de matière donne :

$\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}{\mathrm{]}}_{\text{total}}=\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{]}+\mathrm{[}{\mathrm{H}}_{\mathrm{2}}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}\mathrm{]}+\mathrm{[}{\mathrm{H}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{-}}\mathrm{]}+\mathrm{[}{\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{2}\mathrm{-}}\mathrm{]}$

La loi de Henry peut ainsi être modifiée selon :

$\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}{\mathrm{]}}_{\text{total}}=(K_1+K_1{K}_2+\frac{K_1{K}_2{K}_3}{\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}}+\frac{K_1{K}_2{K}_3{K}_4}{{\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}}^2})P_{\mathrm{C}{\mathrm{O}}_{\mathrm{2}}}$

À pression $P_{\mathrm{C}{\mathrm{O}}_{\mathrm{2}}}$ constante, lorsque $\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}$ diminue, $\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}{\mathrm{]}}_{\text{tot}}$ augmente ; donc plus la solution est basique, plus l'on peut dissoudre de Modèle:Fchim ; inversement, plus la solution est acide, moins l'on peut dissoudre de Modèle:Fchim.

Remarque

Puisque Modèle:Nobr, on a $\mathrm{[}{\mathrm{H}}_{\mathrm{2}}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}\mathrm{]}\mathrm{\ll }\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{]}$ quel que soit le pH. Les deuxième et troisième réactions sont regroupées en une seule réaction :

$\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}\mathrm{+}\mathrm{2}{\mathrm{H}}_{\mathrm{2}}\mathrm{O}\mathrm{\rightleftarrows }{\mathrm{H}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{-}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}\mathrm{+}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{(}\mathrm{a}\mathrm{q}\mathrm{)}$

avec $\mathrm{[}{\mathrm{H}\mathrm{C}{\mathrm{O}}_{\mathrm{3}}}^{\mathrm{-}}\mathrm{]}\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}=K_5\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}\mathrm{]}$ et Modèle:Nobr. Le bilan de matière donne, en négligeant Modèle:Fchim :

$\mathrm{[}\mathrm{C}{\mathrm{O}}_{\mathrm{2}}{\mathrm{]}}_{\text{total}}=(K_1+\frac{K_1{K}_5}{\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}}+\frac{K_1{K}_4{K}_5}{{\mathrm{[}{{\mathrm{H}}_{\mathrm{3}}\mathrm{O}}^{\mathrm{+}}\mathrm{]}}^2})P_{\mathrm{C}{\mathrm{O}}_{\mathrm{2}}}$

Cas des solvants, loi de Raoult

Un Modèle:Nobr est un corps présent dans une solution liquide ayant une fraction molaire très supérieure à celle d'un Modèle:Nobr, soit $x_s^{\text{l}}\gg x_{\sigma }^{\text{l}}$. Ce corps peut quasiment être considéré comme pur, soit $x_s^{\text{l}}\approx 1$.

La relation de Duhem-Margules implique que si un soluté suit la loi de Henry, alors le solvant suit la loi de Raoult qui relie sa pression partielle $P_s$ en phase gazeuse à sa fraction molaire $x_s^{\text{l}}$ en phase liquide à l'équilibre liquide-vapeur selon :

Loi de Raoult : $P_s=x_s^{\text{g}}P=x_s^{\text{l}}{P}_s^{\text{sat}}{𝒫}_s^{\text{l},*}$

avec : Modèle:Début de colonnes

• $P$ la pression totale du mélange ;
• $P_s$ la pression partielle du Modèle:Nobr, par définition $P_s=x_s^{\text{g}}P$ ;
• $P_s^{\text{sat}}$ la pression de vapeur saturante du composé $s$ à la température $T$ du mélange ;
• ${𝒫}_s^{\text{l},*}=\exp \left(\frac{\underset{P_s^{\text{sat}}}{\overset{P}{\int }}{\overline{V}}_s^{\text{l},*}\mathrm{d}P}{RT}\right)$ le facteur de Poynting appliqué au Modèle:Nobr ;
• ${\overline{V}}_s^{\text{l},*}$ le volume molaire du Modèle:Nobr liquide pur ;
• $x_s^{\text{g}}$ la fraction molaire du Modèle:Nobr dans la phase vapeur ;
• $x_s^{\text{l}}$ la fraction molaire du Modèle:Nobr dans la phase liquide.

La relation de Duhem-Margules induit également que si l'on néglige la correction de Poynting pour le soluté, alors elle est également négligeable pour le solvant, soit ${𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}={𝒫}_s^{\text{l},*}=1$.

Extension aux mélanges réels

Coefficients de fugacité et d'activité, équilibres idéaux

À l'équilibre liquide-vapeur, on a pour tout Modèle:Nobr, soluté ou solvant, l'égalité des fugacités en phase vapeur et en phase liquide, soit :

$f_i^{\text{g}}=f_i^{\text{l}}$

La fugacité $f_i^{\text{g}}$ en phase gaz s'écrit à l'aide d'un coefficient de fugacité ${\varphi }_i^{\text{g}}$ corrigeant la loi des gaz parfaits :

$f_i^{\text{g}}=x_i^{\text{g}}{\varphi }_i^{\text{g}}P$

La fugacité $f_i^{\text{l}}$ en phase liquide s'écrit à l'aide d'un coefficient d'activité ${\mathrm{\Gamma }}_i^{\text{l}}$ corrigeant la loi de Lewis et Randall :

$f_i^{\text{l}}=x_i^{\text{l}}{\mathrm{\Gamma }}_i^{\text{l}}{f}_i^{\text{l,*}}$

Ce coefficient est décomposé en deux termes :

${\mathrm{\Gamma }}_i^{\text{l}}={\gamma }_i^{\text{l}}{\hat{\gamma }}_i^{\text{l}}$

Le coefficient ${\gamma }_i^{\text{l}}$ est défini à partir de l'enthalpie libre d'excès $G^{\text{E}}$ :

$RT\ln {\gamma }_i^{\text{l}}\equiv {\left(\frac{\partial G^{\text{E}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}$

Le coefficient ${\hat{\gamma }}_i^{\text{l}}$ permet le calcul du coefficient d'activité dans un mélange multicomposant à partir des propriétés des mélanges binaires « soluté - solvant » ; il vaut 1 dans un mélange binaire.

L'extension des lois de Henry et de Raoult aux mélanges réels nécessite donc une équation d'état pour le calcul des coefficients ${\varphi }_i^{\text{g}}$ et un modèle d'enthalpie libre d'excès pour le calcul des coefficients ${\gamma }_i^{\text{l}}$. Pour des pressions proches de la pression atmosphérique (moins de Modèle:Unité), le gaz se comporte comme un mélange de gaz parfaits, soit ${\varphi }_i^{\text{g}}=1$ pour tout constituant. Dans une solution idéale liquide, on a ${\gamma }_i^{\text{l}}=1$ pour tout constituant. Dans les équilibres idéaux, on a donc, pour tout Modèle:Nobr :

équilibre idéal : $x_i^{\text{g}}P=x_i^{\text{l}}{\hat{\gamma }}_i^{\text{l}}{f}_i^{\text{l,*}}$

Mélanges binaires

On considère un mélange binaire ne comprenant qu'un unique Modèle:Nobr et un unique Modèle:Nobr.

Lois des équilibres binaires

La fugacité (fictive si le soluté est un gaz) $f_{\sigma }^{\text{l,*}}$ du Modèle:Nobr à l'état de liquide pur dans les conditions de pression et température du mélange est donnée par la relation :

$f_{\sigma }^{\text{l,*}}=\frac{k_{\text{H},\sigma ,s}}{{\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}=\frac{k_{\text{H},\sigma ,s}^{\text{sat}}}{{\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$

Pour le Modèle:Nobr, l'équilibre liquide-vapeur est calculé par l'extension de la loi de Henry :

Extension de la loi de Henry aux mélanges binaires réels pour le Modèle:Nobr : $x_{\sigma }^{\text{g}}{\varphi }_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}{\gamma }_{\sigma }^{\text{l}}\frac{k_{\text{H},\sigma ,s}^{\text{sat}}}{{\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$

La fugacité $f_s^{\text{l,*}}$ du Modèle:Nobr à l'état de liquide pur dans les conditions de pression et température du mélange est donnée par la relation :

$f_s^{\text{l,*}}={\varphi }_s^{\text{g,*,sat}}{P}_s^{\text{sat}}{𝒫}_s^{\text{l},*}$

Pour le Modèle:Nobr, l'équilibre liquide-vapeur est calculé par l'extension de la loi de Raoult :

Extension de la loi de Raoult aux mélanges binaires réels pour le Modèle:Nobr : $x_s^{\text{g}}{\varphi }_s^{\text{g}}P=x_s^{\text{l}}{\gamma }_s^{\text{l}}{\varphi }_s^{\text{g,*,sat}}{P}_s^{\text{sat}}{𝒫}_s^{\text{l},*}$

En application du théorème d'Euler, le volume molaire ${\overline{V}}^{\text{l}}$ de la phase liquide vaut :

$\begin{array}{cc}\frac{{\overline{V}}^{\text{l}}}{RT}& =\frac{x_{\sigma }^{\text{l}}{\overline{V}}_{\sigma }^{\text{l}}+x_s^{\text{l}}{\overline{V}}_s^{\text{l}}}{RT}=x_{\sigma }^{\text{l}}{\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}+x_s^{\text{l}}{\left(\frac{\partial \ln f_s^{\text{l}}}{\partial P}\right)}_{T,n}\\ & =x_{\sigma }^{\text{l}}[\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}-{\left(\frac{\partial \ln {\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{\partial P}\right)}_{T,n}+{\left(\frac{\partial \ln {\gamma }_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}]+x_s^{\text{l}}[\frac{{\overline{V}}_s^{\text{l},*}}{RT}+{\left(\frac{\partial \ln {\gamma }_s^{\text{l}}}{\partial P}\right)}_{T,n}]\end{array}$

Le plus souvent, les modèles d'activité ${\gamma }_i^{\text{l}}$ ne dépendent pas de la pression, le volume molaire est alors calculé selon le modèle idéal.

Lois des équilibres binaires idéaux

Pour un mélange liquide binaire idéal aux basses pressions on retrouve les lois idéales :

Loi de Henry pour le Modèle:Nobr : $x_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}^{\text{sat}}{𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$

et :

Loi de Raoult pour le Modèle:Nobr : $x_s^{\text{g}}P=x_s^{\text{l}}{P}_s^{\text{sat}}{𝒫}_s^{\text{l},*}$

Dans une solution idéale, les lois de Henry et de Raoult sont applicables sur l'ensemble des fractions molaires $0\le x_i^{\text{l}}\le 1$. Pour le soluté pur, soit $x_{\sigma }^{\text{l}}=1$, on a $f_{\sigma }^{\text{l,id,*}}=k_{\text{H},\sigma ,s}$. Dans une solution idéale binaire, la constante de Henry est donc à la fois la pente de la fugacité du soluté, qui varie linéairement sur toute la plage de composition, et la fugacité du soluté pur. Aux basses pressions, on a approximativement ${𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\approx {𝒫}_s^{\text{l},*}\approx 1$.

Le volume molaire ${\overline{V}}^{\text{l,id}}$ de la phase liquide idéale vaut :

${\overline{V}}^{\text{l,id}}=x_{\sigma }^{\text{l}}{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}+x_s^{\text{l}}{\overline{V}}_s^{\text{l},*}$

Équations de Krichevsky-Kasarnovsky et Krichevsky-Ilinskaya

Pour un mélange réel, on développe l'égalité des fugacités $f_{\sigma }^{\text{g}}=f_{\sigma }^{\text{l}}$ en supposant que le volume molaire ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ dans le facteur de Poynting ${𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ ne dépend pas de la pression. On obtient, pour tout Modèle:Nobr :

$\ln \frac{f_{\sigma }^{\text{g}}}{x_{\sigma }^{\text{l}}}=\ln \left({\gamma }_{\sigma }^{\text{l}}\frac{k_{\text{H},\sigma ,s}^{\text{sat}}}{{\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\right)=\ln k_{\text{H},\sigma ,s}^{\text{sat}}+\ln {\gamma }_{\sigma }^{\text{l}}-\ln {\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}$

On suppose que le coefficient d'activité suit le modèle de Margules à un paramètre :

$\ln {\gamma }_{\sigma }^{\text{l}}=A{x_s^{\text{l}}}^2$

$\ln {\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 0}{x_s^{\text{l}}\to 1}}{\lim }\ln {\gamma }_{\sigma }^{\text{l}}=A$

avec $x_s^{\text{l}}$ la fraction molaire du Modèle:Nobr dans la phase liquide (avec $x_s^{\text{l}}+x_{\sigma }^{\text{l}}=1$). On obtient l'équation de Krichevsky-IlinskayaModèle:, :

Équation de Krichevsky-Ilinskaya $\ln \frac{f_{\sigma }^{\text{g}}}{x_{\sigma }^{\text{l}}}=\ln k_{\text{H},\sigma ,s}^{\text{sat}}+A\cdot ({x_s^{\text{l}}}^2-1)+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}$

que l'on trouve aussi, en écrivant le modèle de coefficient d'activité selon $RT\ln {\gamma }_{\sigma }^{\text{l}}=A{x_s^{\text{l}}}^2$, sous la formeModèle:,Modèle:,Modèle:, :

Équation de Krichevsky-Ilinskaya $\ln \frac{f_{\sigma }^{\text{g}}}{x_{\sigma }^{\text{l}}}=\ln k_{\text{H},\sigma ,s}^{\text{sat}}+\frac{A}{RT}\cdot ({x_s^{\text{l}}}^2-1)+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}$

Si le mélange liquide est idéal, soit ${\gamma }_{\sigma }^{\text{l}}={\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}=1$ (d'où $A=0$), on obtient l'équation de Krichevsky–KasarnovskyModèle:,Modèle:,Modèle:,Modèle:, :

Équation de Krichevsky–Kasarnovsky $\ln \frac{f_{\sigma }^{\text{g}}}{x_{\sigma }^{\text{l}}}=\ln k_{\text{H},\sigma ,s}^{\text{sat}}+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}$

L'équation de Krichevsky-Kasarnovsky ne s'emploie que pour de faibles concentrations de soluté (solutions liquides idéales), l'équation de Krichevsky-Ilinskaya est valable pour des concentrations plus fortes. Pour des pressions proches de la pression atmosphérique (moins de Modèle:Unité) le gaz se comporte comme un gaz parfait, la fugacité du Modèle:Nobr en phase vapeur est alors égale à sa pression partielle : $f_{\sigma }^{\text{g}}=x_{\sigma }^{\text{g}}P$. Pour des pressions plus importantes, la fugacité du Modèle:Nobr en phase gaz est calculée à l'aide d'un coefficient de fugacité : $f_{\sigma }^{\text{g}}=x_{\sigma }^{\text{g}}{\varphi }_{\sigma }^{\text{g}}P$. Les équations de Krichevsky-Ilinskaya et Krichevsky-Kasarnovsky sont employées pour calculer des solubilités à haute pression, jusqu'à Modèle:Unité environModèle:,.

Paramètres de l'équation de Krichevsky-Ilinskaya pour l'hydrogène Modèle:Fchim dans divers solvants.
$\ln \frac{f_{{\mathrm{H}}_{\mathrm{2}}}^{\text{g}}}{x_{{\mathrm{H}}_{\mathrm{2}}}^{\text{l}}}=\ln k_{\text{H},{\mathrm{H}}_{\mathrm{2}},s}^{\text{sat}}+\frac{A}{RT}\cdot ({x_s^{\text{l}}}^2-1)+\frac{{\overline{V}}_{{\mathrm{H}}_{\mathrm{2}},s}^{\text{l},\mathrm{\infty }}\cdot (P-P_s^{\text{sat}})}{RT}$

Modèle:Nobr	Température $T$ (Modèle:Unité)	Constante de Henry $k_{\text{H},{\mathrm{H}}_{\mathrm{2}},s}^{\text{sat}}$ (Modèle:Unité) (à $P_s^{\text{sat}}\left(T\right)$)	Paramètre de Margules $A$ (Modèle:Unité)	Volume molaire à dilution infinie ${\overline{V}}_{{\mathrm{H}}_{\mathrm{2}},s}^{\text{l},\mathrm{\infty }}$ (Modèle:Unité)
Monoxyde de carbone Modèle:Fchim	68	648	704±69	31,2
	78	476		32,6
	88	405		34,4
Azote Modèle:Fchim	68	547	704±69	30,4
	79	456		31,5
	95	345		34,4
Méthane Modèle:Fchim	90	1848	1486±297	29,7
	110	1050		31,0
	144	638		36,0
Éthane Modèle:Fchim	144	2634	2478±198	37,9
	200	1672		44,2
	228	1226		54,3
Propane Modèle:Fchim	228	1692	2478±198	50
	255	1317		51
	282	1044		63

Mélanges multicomposants

On considère un mélange liquide composé de plusieurs solutés, notés $\sigma$, et plusieurs solvants, notés $s$ ou $r$, aux pression $P$ et température $T$. Il est possible de calculer l'équilibre liquide-vapeur de ce mélange à partir des données d'équilibre de chacun des couples « Modèle:Nobr - Modèle:Nobr ».

Pour les mélanges multicomposants, l'état de référence est le mélange idéal des solvants seuls, en l'absence de tout soluté, autrement dit l'état de dilution infinie de l'ensemble des solutés simultanément. On définit les grandeurs $y_i^{\text{l}}$, $f_{\sigma }^{\text{l,*}}$ et ${\hat{\gamma }}_s^{\text{l}}$ relatives à cet état. Une somme $\sum _s$ ou $\sum _r$ est effectuée sur l'ensemble des solvants du mélange, une somme $\sum _{\sigma }$ sur l'ensemble des solutés du mélange.

Lois des équilibres multicomposants

On définit, pour tout Modèle:Nobr, soluté ou solvant, le facteur :

Pour tout Modèle:Nobr, soluté ou solvant : $y_i^{\text{l}}=\frac{x_i^{\text{l}}}{\sum _s{x}_s^{\text{l}}}$

avec $x_i^{\text{l}}$ la fraction molaire du Modèle:Nobr dans le mélange multicomposant liquide. $y_{\sigma }^{\text{l}}$ est la quantité de Modèle:Nobr rapportée à la quantité totale de solvants dans le mélange liquide. $y_s^{\text{l}}$ est la fraction molaire du Modèle:Nobr dans le mélange liquide des solvants seuls, autrement dit la limite de la fraction molaire $x_s^{\text{l}}$ du Modèle:Nobr à dilution infinie de tous les solutés du mélange liquide multicomposant ; on a $\sum _s{y}_s^{\text{l}}=1$.

On note $f_{\sigma ,s}^{\text{l,*}}$ la fugacité du Modèle:Nobr à l'état de liquide pur calculée à partir des propriétés du mélange binaire « Modèle:Nobr - Modèle:Nobr » :

Pour tout couple Modèle:Nobr - Modèle:Nobr : $f_{\sigma ,s}^{\text{l,*}}=\frac{k_{\text{H},\sigma ,s}}{{\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}$

Si le Modèle:Nobr n'existe pas à l'état de liquide pur dans les conditions de pression $P$ et température $T$ données du mélange, cette fugacité est fictive et sa valeur peut être différente d'un Modèle:Nobr à l'autre. La fugacité $f_{\sigma }^{\text{l,*}}$ de tout Modèle:Nobr à l'état de liquide pur est construite selon l'hypothèse de la solution idéale et peut être fictive :

Pour tout Modèle:Nobr : $\ln f_{\sigma }^{\text{l,*}}=\sum _s{y}_s^{\text{l}}\ln f_{\sigma ,s}^{\text{l,*}}$

Pour tout Modèle:Nobr, l'équilibre liquide-vapeur est calculé par l'extension de la loi de Henry :

Extension de la loi de Henry aux mélanges multicomposants réels pour tout Modèle:Nobr : $x_{\sigma }^{\text{g}}{\varphi }_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}{\gamma }_{\sigma }^{\text{l}}{f}_{\sigma }^{\text{l,*}}$

La fugacité $f_s^{\text{l,*}}$ d'un Modèle:Nobr à l'état de liquide pur dans les conditions de pression et température du mélange est donnée par la relation :

Pour tout Modèle:Nobr : $f_s^{\text{l,*}}={\varphi }_s^{\text{g,*,sat}}{P}_s^{\text{sat}}{𝒫}_s^{\text{l},*}$

Pour tout Modèle:Nobr on pose le facteur correctif ${\hat{\gamma }}_s^{\text{l}}$ :

Pour tout Modèle:Nobr : $\ln {\hat{\gamma }}_s^{\text{l}}=\sum _{\sigma }{y}_{\sigma }^{\text{l}}(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}})$

En l'absence de tout soluté, soit $y_{\sigma }^{\text{l}}=0$ pour tout Modèle:Nobr, ou en présence d'un unique Modèle:Nobr, soit $y_s^{\text{l}}=1$ pour le solvant et $f_{\sigma ,s}^{\text{l,*}}=f_{\sigma }^{\text{l,*}}$ pour tout Modèle:Nobr, on a ${\hat{\gamma }}_s^{\text{l}}=1$. Ceci est à fortiori vrai pour les mélanges binaires « Modèle:Nobr - Modèle:Nobr ».

Pour tout Modèle:Nobr, l'équilibre liquide-vapeur est calculé par l'extension de la loi de Raoult :

Extension de la loi de Raoult aux mélanges multicomposants réels pour tout Modèle:Nobr : $x_s^{\text{g}}{\varphi }_s^{\text{g}}P=x_s^{\text{l}}{\gamma }_s^{\text{l}}{\hat{\gamma }}_s^{\text{l}}{f}_s^{\text{l,*}}$

En application du théorème d'Euler, le volume molaire ${\overline{V}}^{\text{l}}$ de la phase liquide vaut :

$\begin{array}{cc}\frac{{\overline{V}}^{\text{l}}}{RT}& =\frac{\sum _{\sigma }{x}_{\sigma }^{\text{l}}{\overline{V}}_{\sigma }^{\text{l}}+\sum _s{x}_s^{\text{l}}{\overline{V}}_s^{\text{l}}}{RT}=\sum _{\sigma }{x}_{\sigma }^{\text{l}}{\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}+\sum _s{x}_s^{\text{l}}{\left(\frac{\partial \ln f_s^{\text{l}}}{\partial P}\right)}_{T,n}\\ & =\sum _{\sigma }{x}_{\sigma }^{\text{l}}[\sum _s{y}_s^{\text{l}}(\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}-{\left(\frac{\partial \ln {\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{\partial P}\right)}_{T,n})+{\left(\frac{\partial \ln {\gamma }_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}]+\sum _s{x}_s^{\text{l}}[\frac{{\overline{V}}_s^{\text{l},*}}{RT}+{\left(\frac{\partial \ln {\gamma }_s^{\text{l}}}{\partial P}\right)}_{T,n}]\end{array}$

Le plus souvent, les modèles de coefficient d'activité ${\gamma }_i^{\text{l}}$ ne dépendent pas de la pression, le volume molaire est alors calculé selon le modèle idéal.

Le volume molaire partiel ${\overline{V}}_{\sigma }^{\text{l}}$ de tout Modèle:Nobr vaut :

$\begin{array}{cc}\frac{{\overline{V}}_{\sigma }^{\text{l}}}{RT}& ={\left(\frac{\partial \ln f_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}=\underset{\text{à composition constante}=0}{\underbrace{{\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}}}+{\left(\frac{\partial \ln f_{\sigma }^{\text{l},*}}{\partial P}\right)}_{T,n}+{\left(\frac{\partial \ln {\gamma }_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}\\ & =\sum _s{y}_s^{\text{l}}{\left(\frac{\partial \ln f_{\sigma ,s}^{\text{l},*}}{\partial P}\right)}_{T,n}+{\left(\frac{\partial \ln {\gamma }_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}=\sum _s{y}_s^{\text{l}}({\left(\frac{\partial \ln k_{\text{H},\sigma ,s}}{\partial P}\right)}_{T,n}-{\left(\frac{\partial {\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{\partial P}\right)}_{T,n})+{\left(\frac{\partial \ln {\gamma }_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}\\ & =\sum _s{y}_s^{\text{l}}(\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}-{\left(\frac{\partial \ln {\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{\partial P}\right)}_{T,n})+{\left(\frac{\partial \ln {\gamma }_{\sigma }^{\text{l}}}{\partial P}\right)}_{T,n}\end{array}$

Le volume molaire partiel ${\overline{V}}_s^{\text{l}}$ de tout soluté $s$ vaut :

$\begin{array}{cc}\frac{{\overline{V}}_s^{\text{l}}}{RT}& ={\left(\frac{\partial \ln f_s^{\text{l}}}{\partial P}\right)}_{T,n}=\underset{\text{à composition constante}=0}{\underbrace{{\left(\frac{\partial \ln x_s^{\text{l}}}{\partial P}\right)}_{T,n}}}+\underset{\text{à température constante}=0}{\underbrace{{\left(\frac{\partial \ln \left({\varphi }_s^{\text{g,*,sat}}{P}_s^{\text{sat}}\right)}{\partial P}\right)}_{T,n}}}+{\left(\frac{\partial \ln {𝒫}_s^{\text{l},*}}{\partial P}\right)}_{T,n}+{\left(\frac{\partial \ln {\gamma }_s^{\text{l}}}{\partial P}\right)}_{T,n}+{\left(\frac{\partial \ln {\hat{\gamma }}_s^{\text{l}}}{\partial P}\right)}_{T,n}\\ & =\frac{{\overline{V}}_s^{\text{l},*}}{RT}+{\left(\frac{\partial \ln {\gamma }_s^{\text{l}}}{\partial P}\right)}_{T,n}+{\left(\frac{\partial \ln {\hat{\gamma }}_s^{\text{l}}}{\partial P}\right)}_{T,n}\end{array}$

avec :

${\left(\frac{\partial \ln {\hat{\gamma }}_s^{\text{l}}}{\partial P}\right)}_{T,n}=\sum _{\sigma }{y}_{\sigma }^{\text{l}}(\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}-{\left(\frac{\partial \ln {\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{\partial P}\right)}_{T,n}-\sum _r{y}_r^{\text{l}}(\frac{{\overline{V}}_{\sigma ,r}^{\text{l},\mathrm{\infty }}}{RT}-{\left(\frac{\partial \ln {\gamma }_{\sigma ,r}^{\text{l},\mathrm{\infty }}}{\partial P}\right)}_{T,n}))$

Note : dans le calcul du volume molaire ${\overline{V}}^{\text{l}}$, il apparait dans la somme $\sum _s{x}_s^{\text{l}}{\overline{V}}_s^{\text{l}}$ le terme :

$\sum _s{x}_s^{\text{l}}{\left(\frac{\partial \ln {\hat{\gamma }}_s^{\text{l}}}{\partial P}\right)}_{T,n}={\left(\frac{\partial \sum _s{x}_s^{\text{l}}\ln {\hat{\gamma }}_s}{\partial P}\right)}_{T,n}=0$

Composition

Dans le mélange liquide multicomposant, on note $i$ ou $j$ n'importe quel constituant (soluté ou solvant), $\sigma$ un soluté, et $s$ ou $r$ un solvant, avec les quantités :

• $n_i^{\text{l}}$ du Modèle:Nobr ;
• $n^{\text{l}}=\sum _i{n}_i^{\text{l}}=\sum _{\sigma }{n}_{\sigma }^{\text{l}}+\sum _s{n}_s^{\text{l}}$ totale ;
• $n_S^{\text{l}}=\sum _s{n}_s^{\text{l}}$ totale de solvants.

On définit :

Pour tout Modèle:Nobr fraction molaire : $x_i^{\text{l}}\equiv n_i^{\text{l}}/n^{\text{l}}$ $y_i^{\text{l}}\equiv n_i^{\text{l}}/n_S^{\text{l}}$

Pour tout Modèle:Nobr on a :

$y_i^{\text{l}}=\frac{n_i^{\text{l}}}{\sum _s{n}_s^{\text{l}}}=\frac{n_i^{\text{l}}/n^{\text{l}}}{\sum _s(n_s^{\text{l}}/n^{\text{l}})}$

d'où :

$y_i^{\text{l}}=\frac{x_i^{\text{l}}}{\sum _s{x}_s^{\text{l}}}$

Pour tout Modèle:Nobr, $y_s^{\text{l}}=n_s^{\text{l}}/n_S^{\text{l}}$ est la fraction molaire dans le mélange de solvants seuls. On a :

$\sum _s{y}_s^{\text{l}}=\left(\sum _s{n}_s^{\text{l}}\right)/n_S^{\text{l}}$

d'où :

$\sum _s{y}_s^{\text{l}}=1$

Pour tout Modèle:Nobr on a :

$n_j^{\text{l}}{\left(\frac{\partial \ln x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{n_{\ne i}^{\text{l}}}=\frac{n_j^{\text{l}}}{x_j^{\text{l}}}{\left(\frac{\partial x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{n_{\ne i}^{\text{l}}}=n^{\text{l}}{\left(\frac{\partial x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{n_{\ne i}^{\text{l}}}$

${\left(\frac{\partial x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{n_{\ne i}^{\text{l}}}=\{\begin{array}{cc}\frac{1}{n^{\text{l}}}-\frac{x_i^{\text{l}}}{n^{\text{l}}}& \text{si}j=i\\ -\frac{x_j^{\text{l}}}{n^{\text{l}}}& \text{si}j\ne i\end{array}$

$\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{n_{\ne i}^{\text{l}}}=1-\sum _j{x}_j$

d'où :

$\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{n_{\ne i}^{\text{l}}}=0$

Enthalpie libre, potentiel chimique, fugacité

L'enthalpie libre $G^{\text{l}}$ du mélange multicomposant est donnée par le théorème d'Euler :

Théorème d'Euler : $G^{\text{l}}=\sum _i{n}_i^{\text{l}}{\mu }_i^{\text{l}}$

dans lequel, pour tout Modèle:Nobr, le potentiel chimique ${\mu }_i^{\text{l}}$ est défini comme l'enthalpie libre molaire partielle :

Potentiel chimique : ${\mu }_i^{\text{l}}\equiv {\left(\frac{\partial G^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}$

Par définition de la fugacité $f_i^{\text{l}}$, on a, en mélange :

Fugacité : ${\mu }_i^{\text{l}}={\mu }_i^{\bullet ,*}+RT\ln \frac{f_i^{\text{l}}}{P}$

et pour le corps pur :

${\mu }_i^{\text{l,*}}=\underset{\genfrac{}{}{0ex}{}{x_i^{\text{l}}\to 1}{\sum _{j\ne i}{x}_j^{\text{l}}\to 0}}{\lim }{\mu }_i^{\text{l}}={\mu }_i^{\bullet ,*}+RT\ln \frac{f_i^{\text{l,*}}}{P}$

avec : Modèle:Début de colonnes

• $f_i^{\text{l}}$ la fugacité dans le mélange ;
• $f_i^{\text{l,*}}$ la fugacité du corps pur ;
• ${\mu }_i^{\text{l}}$ le potentiel chimique dans le mélange ;
• ${\mu }_i^{\text{l,*}}$ le potentiel chimique du corps pur ;
• ${\mu }_i^{\bullet ,*}$ le potentiel chimique à l'état de gaz parfait pur.

D'autre part, la relation de Gibbs-Duhem induit :

Relation de Gibbs-Duhem à pression et température constantes $\sum _i{n}_i^{\text{l}}\mathrm{d}{\mu }_i^{\text{l}}=0$

d'où :

Pour tout Modèle:Nobr : $\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}=0$

Enthalpie libre d'excès, coefficient d'activité

L'enthalpie libre d'excès $G^{\text{E}}$ est définie par :

Enthalpie libre d'excès : $\frac{G^{\text{E}}}{RT}\equiv \sum _i{n}_i^{\text{l}}\ln \left(\frac{f_i^{\text{l}}}{x_i^{\text{l}}{f}_i^{\text{l,*}}}\right)=\sum _i{n}_i^{\text{l}}\ln {\mathrm{\Gamma }}_i^{\text{l}}$

Le coefficient d'activité ${\mathrm{\Gamma }}_i^{\text{l}}$ est défini par :

Coefficient d'activité : ${\mathrm{\Gamma }}_i^{\text{l}}\equiv \frac{f_i^{\text{l}}}{x_i^{\text{l}}{f}_i^{\text{l,*}}}$

Le coefficient ${\gamma }_i^{\text{l}}$ est défini par :

$\ln {\gamma }_i^{\text{l}}\equiv \frac{{\mu }_i^{\text{E}}}{RT}\equiv {\left(\frac{\partial }{\partial n_i^{\text{l}}}\left[\frac{G^{\text{E}}}{RT}\right]\right)}_{P,T,n_{\ne i}^{\text{l}}}$

avec ${\mu }_i^{\text{E}}$ le potentiel chimique d'excès. Le théorème d'Euler induit :

$\frac{G^{\text{E}}}{RT}=\sum _i{n}_i^{\text{l}}\frac{{\mu }_i^{\text{E}}}{RT}=\sum _i{n}_i^{\text{l}}\ln {\gamma }_i^{\text{l}}$

Les coefficients ${\mathrm{\Gamma }}_i^{\text{l}}$ et ${\gamma }_i^{\text{l}}$ sont définis indépendamment l'un de l'autre. Les deux relations précédentes sur $G^{\text{E}}$ ne permettent pas d'établir d'équivalence entre eux.

D'autre part, la relation de Gibbs-Duhem induit :

Pour tout Modèle:Nobr : $\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln {\gamma }_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}=0$

Fugacité d'un soluté pur

Pour un mélange binaire « Modèle:Nobr - Modèle:Nobr », on définit le potentiel chimique ${\mu }_{\sigma ,s}^{\text{l,*}}$ et la fugacité $f_{\sigma ,s}^{\text{l,*}}$ du soluté pur par la relation :

${\mu }_{\sigma ,s}^{\text{l,*}}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 1}{x_s^{\text{l}}\to 0}}{\lim }{\mu }_{\sigma }^{\text{l}}={\mu }_{\sigma }^{\bullet ,*}+RT\ln \frac{f_{\sigma ,s}^{\text{l,*}}}{P}$

Le potentiel ${\mu }_{\sigma }^{\bullet ,*}$ du soluté à l'état de gaz parfait pur est défini indépendamment du solvant considéré. Pour le même soluté dans un mélange de solvants, on définit sa fugacité à l'état de corps pur par :

${\mu }_{\sigma }^{\text{l,*}}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}\to 1}{\sum _s{x}_s^{\text{l}}\to 0}}{\lim }{\mu }_{\sigma }^{\text{l}}={\mu }_{\sigma }^{\bullet ,*}+RT\ln \frac{f_{\sigma }^{\text{l,*}}}{P}$

Si dans le mélange de solvants le Modèle:Nobr est représenté par la fraction molaire $y_s^{\text{l}}$, on a :

${\mu }_{\sigma }^{\text{l,*}}=\sum _s{y}_s^{\text{l}}{\mu }_{\sigma ,s}^{\text{l,*}}=\underset{=1}{\underbrace{\left(\sum _s{y}_s^{\text{l}}\right)}}{\mu }_{\sigma }^{\bullet ,*}+RT(\sum _s{y}_s^{\text{l}}\ln f_{\sigma ,s}^{\text{l,*}})-\underset{=1}{\underbrace{\left(\sum _s{y}_s^{\text{l}}\right)}}RT\ln P={\mu }_{\sigma }^{\bullet ,*}+RT(\sum _s{y}_s^{\text{l}}\ln f_{\sigma ,s}^{\text{l,*}})-RT\ln P$

On obtient ainsi la fugacité $f_{\sigma }^{\text{l,*}}$ pour le mélange multicomposant à partir des propriétés des mélanges binaires :

Pour tout soluté $\sigma$ : $\ln f_{\sigma }^{\text{l,*}}=\sum _s{y}_s^{\text{l}}\ln f_{\sigma ,s}^{\text{l,*}}$

La fugacité $f_{\sigma ,s}^{\text{l,*}}$ d'un soluté pur par rapport à un solvant ne dépend pas de la composition. La fugacité $f_{\sigma }^{\text{l,*}}$ d'un soluté pur par rapport à un mélange de solvants dépend des quantités des solvants. La fugacité $f_s^{\text{l,*}}$ d'un solvant pur ne dépend pas de la composition.

Fugacité d'un corps en mélange

Pour tout Modèle:Nobr on a :

$\begin{array}{cc}{\left(\frac{\partial }{\partial n_i^{\text{l}}}\left[\frac{G^{\text{E}}}{RT}\right]\right)}_{P,T,n_{\ne i}^{\text{l}}}& ={\left(\frac{\partial }{\partial n_i^{\text{l}}}[\sum _j{n}_j^{\text{l}}\ln \left(\frac{f_j^{\text{l}}}{x_j^{\text{l}}{f}_j^{\text{l,*}}}\right)]\right)}_{P,T,n_{\ne i}^{\text{l}}}\\ & =\ln \left(\frac{f_i^{\text{l}}}{x_i^{\text{l}}{f}_i^{\text{l,*}}}\right)+\underset{\text{Gibbs-Duhem}=0}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}-\underset{=0}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}-\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l,*}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}\end{array}$

On définit le coefficient ${\hat{\gamma }}_i^{\text{l}}$ par :

$\ln {\hat{\gamma }}_i^{\text{l}}\equiv \sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l,*}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}$

On peut donc écrire :

$\ln {\gamma }_i^{\text{l}}=\ln \left(\frac{f_i^{\text{l}}}{x_i^{\text{l}}{f}_i^{\text{l,*}}}\right)-\ln {\hat{\gamma }}_i^{\text{l}}$

On obtient :

Pour tout Modèle:Nobr : $f_i^{\text{l}}=x_i^{\text{l}}{\gamma }_i^{\text{l}}{\hat{\gamma }}_i^{\text{l}}{f}_i^{\text{l,*}}$

soit :

${\mathrm{\Gamma }}_i^{\text{l}}=\frac{f_i^{\text{l}}}{x_i^{\text{l}}{f}_i^{\text{l,*}}}={\gamma }_i^{\text{l}}{\hat{\gamma }}_i^{\text{l}}$

On a ${\mathrm{\Gamma }}_i^{\text{l}}={\gamma }_i^{\text{l}}$ si ${\hat{\gamma }}_i^{\text{l}}=1$, c'est-à-dire si les $f_j^{\text{l,*}}$ ne dépendent pas de $n_i^{\text{l}}$.

Coefficient ${\hat{\gamma }}_i^{\text{l}}$

Pour un Modèle:Nobr, les fugacités $f_j^{\text{l,*}}$ ne dépendent pas de sa quantité $n_{\sigma }^{\text{l}}$. En conséquence :

$\ln {\hat{\gamma }}_{\sigma }^{\text{l}}=\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l,*}}}{\partial n_{\sigma }^{\text{l}}}\right)}_{P,T,n_{\ne \sigma }^{\text{l}}}=0$

On obtient :

Pour tout Modèle:Nobr : ${\hat{\gamma }}_{\sigma }^{\text{l}}=1$

Pour un Modèle:Nobr, seules les fugacités $f_{\sigma }^{\text{l,*}}$ dépendent de sa quantité $n_s^{\text{l}}$. En conséquence :

$\ln {\hat{\gamma }}_s^{\text{l}}=\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l,*}}}{\partial n_s^{\text{l}}}\right)}_{P,T,n_{\ne s}^{\text{l}}}=\sum _{\sigma }{n}_{\sigma }^{\text{l}}{\left(\frac{\partial \ln f_{\sigma }^{\text{l,*}}}{n_s^{\text{l}}}\right)}_{P,T,n_{\ne s}^{\text{l}}}$

Par construction on a :

$\ln f_{\sigma }^{\text{l,*}}=\sum _r{y}_r^{\text{l}}\ln f_{\sigma ,r}^{\text{l,*}}=\sum _r\frac{n_r^{\text{l}}}{n_S^{\text{l}}}\ln f_{\sigma ,r}^{\text{l,*}}$

Les $f_{\sigma ,r}^{\text{l,*}}$ ne dépendant pas de la composition, on a :

${\left(\frac{\partial \ln f_{\sigma }^{\text{l,*}}}{\partial n_s^{\text{l}}}\right)}_{P,T,n_{\ne s}^{\text{l}}}=\frac{1}{n_S^{\text{l}}}\ln f_{\sigma ,s}^{\text{l,*}}-\sum _r\frac{n_r}{{n_S^{\text{l}}}^2}\ln f_{\sigma ,r}^{\text{l,*}}=\frac{1}{n_S^{\text{l}}}(\ln f_{\sigma ,s}^{\text{l,*}}-\sum _r{y}_r^{\text{l}}\ln f_{\sigma ,r}^{\text{l,*}})=\frac{1}{n_S^{\text{l}}}(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}})$

La somme pondérée par les quantités $n_{\sigma }^{\text{l}}$, avec $y_{\sigma }^{\text{l}}=n_{\sigma }^{\text{l}}/n_S^{\text{l}}$, donne :

Pour tout Modèle:Nobr : $\ln {\hat{\gamma }}_s^{\text{l}}=\sum _{\sigma }{y}_{\sigma }^{\text{l}}(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}})$

On introduit $\sum _r{y}_r^{\text{l}}=1$ et $\ln f_{\sigma }^{\text{l,*}}=\sum _s{y}_r^{\text{l}}\ln f_{\sigma ,r}^{\text{l,*}}$ :

$\ln {\hat{\gamma }}_s^{\text{l}}=\sum _{\sigma }{y}_{\sigma }^{\text{l}}\left(\sum _r{y}_r^{\text{l}}\right)\ln f_{\sigma ,s}^{\text{l,*}}-\sum _{\sigma }{y}_{\sigma }^{\text{l}}(\sum _r{y}_r^{\text{l}}\ln f_{\sigma ,r}^{\text{l,*}})$

On a donc également :

$\ln {\hat{\gamma }}_s^{\text{l}}=\sum _{\sigma }\sum _r{y}_{\sigma }^{\text{l}}{y}_r^{\text{l}}(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma ,r}^{\text{l,*}})$

Vérification du théorème d'Euler pour l'enthalpie libre d'excès

Pour tout Modèle:Nobr on a ${\mathrm{\Gamma }}_i^{\text{l}}={\gamma }_i^{\text{l}}{\hat{\gamma }}_i^{\text{l}}$. Par définition l'enthalpie libre d'excès donne :

$\frac{G^{\text{E}}}{RT}=\sum _i{n}_i^{\text{l}}\ln {\mathrm{\Gamma }}_i^{\text{l}}=\sum _i{n}_i^{\text{l}}\ln {\gamma }_i^{\text{l}}+\sum _i{n}_i^{\text{l}}\ln {\hat{\gamma }}_i^{\text{l}}$

Puisque par définition $RT\ln {\gamma }_i^{\text{l}}\equiv {\mu }_i^{\text{E}}$, elle doit aussi vérifier le théorème d'Euler :

$\frac{G^{\text{E}}}{RT}=\sum _i{n}_i^{\text{l}}\frac{{\mu }_i^{\text{E}}}{RT}=\sum _i{n}_i^{\text{l}}\ln {\gamma }_i^{\text{l}}$

Ceci impose de vérifier que :

$\sum _i{n}_i^{\text{l}}\ln {\hat{\gamma }}_i^{\text{l}}=0$

La somme pondérée par les quantités $n_s^{\text{l}}$ des coefficients $\ln {\hat{\gamma }}_s^{\text{l}}$ trouvés précédemment donne :

$\sum _s{n}_s^{\text{l}}\ln {\hat{\gamma }}_s^{\text{l}}=\sum _s\sum _{\sigma }{n}_s^{\text{l}}{y}_{\sigma }^{\text{l}}(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}})$

et, avec $n_s^{\text{l}}{y}_{\sigma }^{\text{l}}=y_s^{\text{l}}{n}_{\sigma }^{\text{l}}$ :

$\sum _s{n}_s^{\text{l}}\ln {\hat{\gamma }}_s^{\text{l}}=\sum _s\sum _{\sigma }{y}_s^{\text{l}}{n}_{\sigma }^{\text{l}}(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}})=\sum _{\sigma }{n}_{\sigma }^{\text{l}}\underset{=\ln f_{\sigma }^{\text{l,*}}}{\underbrace{(\sum _s{y}_s^{\text{l}}\ln f_{\sigma ,s}^{\text{l},*})}}-\underset{=1}{\underbrace{\left(\sum _s{y}_s^{\text{l}}\right)}}\sum _{\sigma }{n}_{\sigma }^{\text{l}}\ln f_{\sigma }^{\text{l,*}}=0$

Puisque ${\hat{\gamma }}_{\sigma }^{\text{l}}=1$ pour tout soluté, on peut généraliser :

$\sum _i{n}_i^{\text{l}}\ln {\hat{\gamma }}_i^{\text{l}}=0$

Le théorème d'Euler est donc vérifié pour l'enthalpie libre d'excès.

Vérification de la relation de Gibbs-Duhem sur les fugacités

Pour tout Modèle:Nobr on a $f_i^{\text{l}}=x_i^{\text{l}}{\gamma }_i^{\text{l}}{\hat{\gamma }}_i^{\text{l}}{f}_i^{\text{l,*}}$. Par conséquent :

$\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}=\underset{=0}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}+\underset{\text{Gibbs-Duhem}=0}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln {\gamma }_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}+\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln {\hat{\gamma }}_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}+\underset{\text{définition}=\ln {\hat{\gamma }}_i^{\text{l}}}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l,*}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}$

On a :

$\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln {\hat{\gamma }}_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}=\underset{=0}{\underbrace{{\left(\frac{\partial }{\partial n_i^{\text{l}}}[\sum _j{n}_j^{\text{l}}\ln {\hat{\gamma }}_j^{\text{l}}]\right)}_{P,T,n_{\ne i}^{\text{l}}}}}-\ln {\hat{\gamma }}_i^{\text{l}}=-\ln {\hat{\gamma }}_i^{\text{l}}$

On vérifie en conséquence la relation de Gibbs-Duhem sur les fugacités pour tout Modèle:Nobr :

$\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}=0$

Lois des équilibres multicomposants idéaux

Comme pour un mélange binaire, la fugacité (fictive pour les gaz) $f_{\sigma }^{\text{l,*}}$ du Modèle:Nobr à l'état de liquide pur est calculée selon la relation :

$f_{\sigma }^{\text{l,*}}=\frac{k_{\text{H},\sigma ,m}}{{\gamma }_{\sigma ,m}^{\text{l},\mathrm{\infty }}}$

avec :

• $k_{\text{H},\sigma ,m}$ la constante de Henry du Modèle:Nobr dans le mélange de solvants seuls ;
• ${\gamma }_{\sigma ,m}^{\text{l},\mathrm{\infty }}$ le coefficient d'activité du Modèle:Nobr à dilution infinie dans le mélange de solvants seuls.

La relation $\ln f_{\sigma }^{\text{l,*}}=\sum _s{y}_s^{\text{l}}\ln f_{\sigma ,s}^{\text{l,*}}$ donne pour tout Modèle:Nobr la relation idéale de Krichevsky :

Équation de Krichevsky pour tout Modèle:Nobr : $\ln k_{\text{H},\sigma ,m}=\sum _s{y}_s^{\text{l}}\ln k_{\text{H},\sigma ,s}$

Pour les basses pressions et les phases liquides idéales on a par conséquent la loi idéale pour les solutés :

Loi de Henry étendue aux mélanges multicomposants idéaux pour tout Modèle:Nobr : $x_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,m}$

On a par ailleurs pour tout Modèle:Nobr :

pour tout Modèle:Nobr : $\ln {\hat{\gamma }}_s^{\text{l}}=\sum _{\sigma }{y}_{\sigma }^{\text{l}}(\ln k_{\text{H},\sigma ,s}-\ln k_{\text{H},\sigma ,m})$

Pour les basses pressions et les phases liquides idéales on a par conséquent la loi idéale pour les solvants :

Loi de Raoult étendue aux mélanges multicomposants idéaux pour tout Modèle:Nobr : $x_s^{\text{g}}P=x_s^{\text{l}}{\hat{\gamma }}_s^{\text{l}}{P}_s^{\text{sat}}{𝒫}_s^{\text{l},*}$

Le volume molaire ${\overline{V}}^{\text{l,id}}$ de la phase liquide idéale vaut :

${\overline{V}}^{\text{l,id}}=\sum _{\sigma }\sum _s{x}_{\sigma }^{\text{l}}{y}_s^{\text{l}}{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}+\sum _s{x}_s^{\text{l}}{\overline{V}}_s^{\text{l},*}$

L'enthalpie libre d'excès $G^{\text{E}}$ vaut :

$\frac{G^{\text{E}}}{RT}=\sum _i{n}_i^{\text{l}}\ln \left(\frac{f_i^{\text{l}}}{x_i^{\text{l}}{f}_i^{\text{l,*}}}\right)=\sum _i{n}_i^{\text{l}}\ln {\mathrm{\Gamma }}_i^{\text{l}}=\sum _i{n}_i^{\text{l}}\ln {\gamma }_i^{\text{l}}$

Une solution est idéale si $G^{\text{E}}=0$, avec ${\gamma }_i^{\text{l}}=1$ pour tous les constituants. L'enthalpie libre d'excès $G^{\text{E}}$ est par définition l'écart entre l'enthalpie libre $G^{\text{l}}$ du mélange réel et l'enthalpie libre $G^{\text{l,id}}$ de la solution idéale correspondante :

$G^{\text{E}}\equiv G^{\text{l}}-G^{\text{l,id}}$

Le théorème d'Euler donne :

$G^{\text{l}}=\sum _i{n}_i^{\text{l}}{\mu }_i^{\text{l}}=\sum _i{n}_i^{\text{l}}({\mu }_i^{\bullet ,*}+RT\ln \frac{f_i^{\text{l}}}{P})$

en conséquence :

$G^{\text{l,id}}=G^{\text{l}}-G^{\text{E}}=\sum _i{n}_i^{\text{l}}({\mu }_i^{\bullet ,*}+RT\ln \left(\frac{x_i^{\text{l}}{f}_i^{\text{l,*}}}{P}\right))$

On calcule le potentiel chimique idéal ${\mu }_i^{\text{l,id}}$ :

$\frac{{\mu }_i^{\text{l,id}}}{RT}={\left(\frac{\partial }{\partial n_i^{\text{l}}}\left[\frac{G^{\text{l,id}}}{RT}\right]\right)}_{P,T,n_{\ne i}^{\text{l}}}=\frac{{\mu }_i^{\bullet ,*}}{RT}+\ln \left(\frac{x_i^{\text{l}}{f}_i^{\text{l,*}}}{P}\right)+\underset{=0}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln x_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}+\underset{\text{définition}=\ln {\hat{\gamma }}_i^{\text{l}}}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l,*}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}$

Par conséquent, on a le potentiel chimique, la fugacité et le coefficient d'activité idéaux :

${\mu }_i^{\text{l,id}}={\mu }_i^{\bullet ,*}+RT\ln \frac{f_i^{\text{l,id}}}{P}$

$f_i^{\text{l,id}}=x_i^{\text{l}}{\hat{\gamma }}_i^{\text{l}}{f}_i^{\text{l,*}}$

${\mathrm{\Gamma }}_i^{\text{l,id}}=\frac{f_i^{\text{l,id}}}{x_i^{\text{l}}{f}_i^{\text{l,*}}}={\hat{\gamma }}_i^{\text{l}}$

On vérifie le théorème d'Euler pour l'enthalpie libre idéale :

$\sum _i{n}_i^{\text{l}}{\mu }_i^{\text{l,id}}=\sum _i{n}_i^{\text{l}}({\mu }_i^{\bullet ,*}+RT\ln \left(\frac{x_i^{\text{l}}{f}_i^{\text{l,*}}}{P}\right))+RT\underset{=0}{\underbrace{\sum _i{n}_i^{\text{l}}\ln {\hat{\gamma }}_i^{\text{l}}}}=G^{\text{l,id}}$

et pour tout Modèle:Nobr, puisque $f_i^{\text{l}}={\gamma }_i^{\text{l}}{f}_i^{\text{l,id}}$, la relation de Gibbs-Duhem :

$\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l,id}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}=\underset{\text{Gibbs-Duhem}=0}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln f_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}-\underset{\text{Gibbs-Duhem}=0}{\underbrace{\sum _j{n}_j^{\text{l}}{\left(\frac{\partial \ln {\gamma }_j^{\text{l}}}{\partial n_i^{\text{l}}}\right)}_{P,T,n_{\ne i}^{\text{l}}}}}=0$

Application à un mélange ternaire : un soluté et deux solvants

Soit un mélange liquide ternaire constitué d'un Modèle:Nobr et de deux solvants, notés $s1$ et $s2$. On suppose que l'on dispose des constantes de Henry $k_{\text{H},\sigma ,s1}$ et $k_{\text{H},\sigma ,s2}$ dans les solutions binaires, et que l'enthalpie libre d'excès molaire ${\overline{G}}^{\text{E}}$ du mélange ternaire est calculable selon l'extension du modèle de Margules :

$\frac{{\overline{G}}^{\text{E}}}{RT}=A_{\sigma ,s1}{x}_{\sigma }^{\text{l}}{x}_{s1}^{\text{l}}+A_{\sigma ,s2}{x}_{\sigma }^{\text{l}}{x}_{s2}^{\text{l}}+A_{s1,s2}{x}_{s1}^{\text{l}}{x}_{s2}^{\text{l}}$

Les coefficients d'activité des trois corps sont donnés par :

coefficients d'activité $\ln {\gamma }_{\sigma }^{\text{l}}=A_{\sigma ,s1}{x}_{s1}^{\text{l}}+A_{\sigma ,s2}{x}_{s2}^{\text{l}}-{\overline{G}}^{\text{E}}/RT$ $\ln {\gamma }_{s1}^{\text{l}}=A_{\sigma ,s1}{x}_{\sigma }^{\text{l}}+A_{s1,s2}{x}_{s2}^{\text{l}}-{\overline{G}}^{\text{E}}/RT$ $\ln {\gamma }_{s2}^{\text{l}}=A_{\sigma ,s2}{x}_{\sigma }^{\text{l}}+A_{s1,s2}{x}_{s1}^{\text{l}}-{\overline{G}}^{\text{E}}/RT$

Les coefficients d'activité du Modèle:Nobr à dilution infinie dans chacun des deux solvants $s1$ et $s2$ sont donnés par :

$\ln {\gamma }_{\sigma ,s1}^{\text{l},\mathrm{\infty }}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}+x_{s2}^{\text{l}}\to 0}{x_{s1}^{\text{l}}\to 1}}{\lim }\ln {\gamma }_{\sigma }^{\text{l}}=A_{\sigma ,s1}$

$\ln {\gamma }_{\sigma ,s2}^{\text{l},\mathrm{\infty }}=\underset{\genfrac{}{}{0ex}{}{x_{\sigma }^{\text{l}}+x_{s1}^{\text{l}}\to 0}{x_{s2}^{\text{l}}\to 1}}{\lim }\ln {\gamma }_{\sigma }^{\text{l}}=A_{\sigma ,s2}$

On pose :

$y_{\sigma }^{\text{l}}=\frac{x_{\sigma }^{\text{l}}}{x_{s1}^{\text{l}}+x_{s2}^{\text{l}}};y_{s1}^{\text{l}}=\frac{x_{s1}^{\text{l}}}{x_{s1}^{\text{l}}+x_{s2}^{\text{l}}};y_{s2}^{\text{l}}=\frac{x_{s2}^{\text{l}}}{x_{s1}^{\text{l}}+x_{s2}^{\text{l}}}$

La fugacité $f_{\sigma }^{\text{l,*}}$ du Modèle:Nobr à l'état de liquide pur est calculée selon :

$\ln f_{\sigma }^{\text{l,*}}=y_{s1}^{\text{l}}\ln f_{\sigma ,s1}^{\text{l,*}}+y_{s2}^{\text{l}}\ln f_{\sigma ,s2}^{\text{l,*}}=y_{s1}^{\text{l}}\ln \frac{k_{\text{H},\sigma ,s1}}{{\gamma }_{\sigma ,s1}^{\text{l},\mathrm{\infty }}}+y_{s2}^{\text{l}}\ln \frac{k_{\text{H},\sigma ,s2}}{{\gamma }_{\sigma ,s2}^{\text{l},\mathrm{\infty }}}$

d'où :

pour le Modèle:Nobr $\ln f_{\sigma }^{\text{l,*}}=y_{s1}^{\text{l}}(\ln k_{\text{H},\sigma ,s1}-A_{\sigma ,s1})+y_{s2}^{\text{l}}(\ln k_{\text{H},\sigma ,s2}-A_{\sigma ,s2})$

Pour les deux Modèle:Nobr on a :

$\ln {\hat{\gamma }}_{s1}^{\text{l}}=y_{\sigma }^{\text{l}}(\ln f_{\sigma ,s1}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}})=y_{\sigma }^{\text{l}}{y}_{s2}^{\text{l}}(\ln \left(\frac{k_{\text{H},\sigma ,s1}}{{\gamma }_{\sigma ,s1}^{\text{l},\mathrm{\infty }}}\right)-\ln \left(\frac{k_{\text{H},\sigma ,s2}}{{\gamma }_{\sigma ,s2}^{\text{l},\mathrm{\infty }}}\right))$

$\ln {\hat{\gamma }}_{s2}^{\text{l}}=y_{\sigma }^{\text{l}}(\ln f_{\sigma ,s2}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}})=y_{\sigma }^{\text{l}}{y}_{s1}^{\text{l}}(\ln \left(\frac{k_{\text{H},\sigma ,s2}}{{\gamma }_{\sigma ,s2}^{\text{l},\mathrm{\infty }}}\right)-\ln \left(\frac{k_{\text{H},\sigma ,s1}}{{\gamma }_{\sigma ,s1}^{\text{l},\mathrm{\infty }}}\right))$

d'où :

pour les deux solvants $s1$ et $s2$ $\ln {\hat{\gamma }}_{s1}^{\text{l}}=y_{\sigma }^{\text{l}}{y}_{s2}^{\text{l}}(\ln k_{\text{H},\sigma ,s1}-\ln k_{\text{H},\sigma ,s2}+A_{\sigma ,s2}-A_{\sigma ,s1})$ $\ln {\hat{\gamma }}_{s2}^{\text{l}}=y_{\sigma }^{\text{l}}{y}_{s1}^{\text{l}}(\ln k_{\text{H},\sigma ,s2}-\ln k_{\text{H},\sigma ,s1}+A_{\sigma ,s1}-A_{\sigma ,s2})$

Pour une solution idéale, soit ${\overline{G}}^{\text{E}}=0$, on a ${\gamma }_{\sigma }^{\text{l}}={\gamma }_{s1}^{\text{l}}={\gamma }_{s2}^{\text{l}}={\gamma }_{\sigma ,s1}^{\text{l},\mathrm{\infty }}={\gamma }_{\sigma ,s2}^{\text{l},\mathrm{\infty }}={\gamma }_{\sigma ,m}^{\text{l},\mathrm{\infty }}=1$ avec $A_{\sigma ,s1}=A_{\sigma ,s2}=A_{s1,s2}=0$. On obtient :

$\ln k_{\text{H},\sigma ,m}=y_{s1}^{\text{l}}\ln k_{\text{H},\sigma ,s1}+y_{s2}^{\text{l}}\ln k_{\text{H},\sigma ,s2}$ (équation de Krichevsky)

$\ln {\hat{\gamma }}_{s1}^{\text{l}}=y_{\sigma }^{\text{l}}{y}_{s2}^{\text{l}}(\ln k_{\text{H},\sigma ,s1}-\ln k_{\text{H},\sigma ,s2})$

$\ln {\hat{\gamma }}_{s2}^{\text{l}}=y_{\sigma }^{\text{l}}{y}_{s1}^{\text{l}}(\ln k_{\text{H},\sigma ,s2}-\ln k_{\text{H},\sigma ,s1})$

L'ensemble de ces expressions peut être étendu à des mélanges contenant plus de solutés et de solvants, avec :

$\frac{{\overline{G}}^{\text{E}}}{RT}=\sum _i\sum _{j>i}{A}_{i,j}{x}_i^{\text{l}}{x}_j^{\text{l}}$

Il est alors pratique de négliger les interactions entre solutés si ceux-ci ne sont que faiblement concentrés, soit $A_{i,j}=0$ si $i$ et $j$ sont deux solutés :

$\frac{{\overline{G}}^{\text{E}}}{RT}=\sum _{\sigma }\sum _s{A}_{\sigma ,s}{x}_{\sigma }^{\text{l}}{x}_s^{\text{l}}+\sum _s\sum _{r>s}{A}_{s,r}{x}_s^{\text{l}}{x}_r^{\text{l}}$

Applications

Étude de la solubilité

La loi de Henry, sous la forme $P_{\sigma }=x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}$, permet d'établir l'évolution de la solubilité $x_{\sigma }^{\text{l}}$ en fonction de la pression et de la température dans son domaine d'application, aux basses pressions (loi des gaz parfaits) et faibles solubilités (solution idéale). On suppose également que le Modèle:Nobr est très peu volatil, voire que le Modèle:Nobr est seul en phase gaz, soit, pour la pression partielle du soluté, $P_{\sigma }\approx P$, la pression totale.

En fonction de la pression

En dérivant, à température constante, l'expression de la loi de Henry par rapport à la pression, on obtient :

${\left(\frac{\partial \ln P_{\sigma }}{\partial P}\right)}_T={\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_T+{\left(\frac{\partial \ln k_{\text{H},\sigma ,s}}{\partial P}\right)}_T={\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_T+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}$

avec ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ le volume molaire partiel du Modèle:Nobr à dilution infinie dans le Modèle:Nobr (voir le paragraphe Constante de Henry - Dépendance à la pression). Le Modèle:Nobr étant négligeable en phase gaz, on a $P_{\sigma }=P$ et donc :

${\left(\frac{\partial \ln P_{\sigma }}{\partial P}\right)}_T=\frac{1}{P}$

On obtient :

${\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_T=\frac{1}{P}-\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}$

En application de la loi des gaz parfaits, aux basses pressions le volume molaire de la phase gaz vaut : ${\overline{V}}^{\bullet }=RT/P$. D'autre part, toujours aux basses pressions, le volume molaire de la phase gaz est très supérieur au volume molaire partiel en phase liquide, soit ${\overline{V}}^{\bullet }\gg {\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$, d'où :

Aux basses pressions : ${\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_T\approx \frac{{\overline{V}}^{\bullet }}{RT}=\frac{1}{P}$

Avec un point de solubilité connu $(x_{\sigma }^{\text{l},\circ },P^{\circ },T)$ on intègre à température constante :

$\frac{x_{\sigma }^{\text{l}}}{x_{\sigma }^{\text{l},\circ }}=\frac{P}{P^{\circ }}$

Aux basses pressions, la solubilité $x_{\sigma }^{\text{l}}$ augmente linéairement avec la pression totale $P$ lorsque le Modèle:Nobr est seul en phase gaz.

Le volume molaire d'un gaz diminue avec une augmentation de pression ; par contre on peut considérer que ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ est indépendant de la pression, les liquides étant très peu compressibles. Aux pressions élevées, le modèle des gaz parfaits ne s'applique plus, la loi de Henry est appliquée sous la forme $f_{\sigma }^{\text{g}}=x_{\sigma }^{\text{l}}{k}_{\text{H},\sigma ,s}$. Le Modèle:Nobr étant négligé en phase gaz, la fugacité du Modèle:Nobr dans cette phase varie en fonction de la pression selon :

${\left(\frac{\partial \ln f_{\sigma }^{\text{g}}}{\partial P}\right)}_T=\frac{{\overline{V}}_{\sigma }^{\text{g},*}}{RT}={\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_T+\frac{{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}$

avec ${\overline{V}}_{\sigma }^{\text{g},*}$ le volume molaire du soluté pur à l'état gazeux dans les conditions de pression et température du mélange. On obtient par conséquentModèle:, :

Aux hautes pressions : ${\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial P}\right)}_T=\frac{{\overline{V}}_{\sigma }^{\text{g},*}-{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}}{RT}$

Lorsque ${\overline{V}}_{\sigma }^{\text{g},*}>{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$, aux basses pressions, la solubilité augmente avec la pression ; elle décroît lorsque ${\overline{V}}_{\sigma }^{\text{g},*}<{\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$, aux très hautes pressions (lorsque le volume molaire de la phase gaz est faible). Ainsi, lorsque ${\overline{V}}_{\sigma }^{\text{g},*}={\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$, la solubilité $x_{\sigma }^{\text{l}}$ atteint un maximum en fonction de la pression. Ceci a été vérifié expérimentalement pour la solubilité de l'azote dans l'eau, qui atteint à Modèle:Unité un maximum à environ Modèle:UnitéModèle:,.

En fonction de la température

Solubilité du dioxyde de carbone dans l'eau en fonction de la température à pression atmosphérique.

En dérivant, à pression constante, l'expression de la loi de Henry par rapport à la température, on obtient :

${\left(\frac{\partial \ln P_{\sigma }}{\partial \frac{1}{T}}\right)}_P={\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial \frac{1}{T}}\right)}_P+{\left(\frac{\partial \ln k_{\text{H},\sigma ,s}}{\partial \frac{1}{T}}\right)}_P={\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial \frac{1}{T}}\right)}_P+\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}$

avec ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}$ l'enthalpie de dissolution du Modèle:Nobr dans le Modèle:Nobr (voir le paragraphe Constante de Henry - Dépendance à la température). Le Modèle:Nobr étant négligeable en phase gaz, on a $P_{\sigma }=P$ et donc :

${\left(\frac{\partial \ln P_{\sigma }}{\partial \frac{1}{T}}\right)}_P=0$

On obtient :

${\left(\frac{\partial \ln x_{\sigma }^{\text{l}}}{\partial \frac{1}{T}}\right)}_P=-\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}$

Diverses formes de la solubilité en fonction de la température sont communément employéesModèle:, :

$\ln x_{\sigma }^{\text{l}}=a+\frac{b}{T}+c\ln T+dT$

$\ln x_{\sigma }^{\text{l}}=a+\frac{b}{T}+\frac{c}{T^2}+\frac{d}{T^3}$

avec $a$, $b$, $c$ et $d$ des constantes empiriques spécifiques du couple « Modèle:Nobr - Modèle:Nobr ». Les expressions des enthalpies de dissolution respectives sont :

$-\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}=b-cT-d{T}^2$

$-\frac{{\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}}{R}=b+\frac{2c}{T}+\frac{3d}{T^2}$

Pour la plupart des gaz, la dissolution est exothermique aux basses pressions et températures, soit ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}<0$, par conséquent la solubilité diminue lorsque la température augmente. De nombreux gaz présentent un minimum de solubilité, la solubilité augmentant après avoir diminué lorsque la température augmenteModèle:,. Ainsi, aux basses pressions, le minimum de solubilité de l'hélium dans l'eau se situe à environ Modèle:Unité, ceux de l'argon, de l'oxygène et de l'azote se situent entre 92 et Modèle:Unité et celui du xénon à environ Modèle:Unité.

Solubilité de l'oxygène Modèle:Fchim dans l'eau
pour une pression partielle d'oxygène $P_{{\mathrm{O}}_{\mathrm{2}}}$ de Modèle:Unité.

Température $T$ (Modèle:Unité)	Solubilité Modèle:Nobr	Enthalpie de dissolution ${\mathrm{\Delta }}_{\text{sol}}{H}_{{\mathrm{O}}_{\mathrm{2}},\text{eau}}$ (Modèle:Unité)	Constante de Henry Modèle:Nobr (Modèle:Unité) (à $P_{\text{eau}}^{\text{sat}}\left(T\right)$)
273,15	3,9594	-17,43	2,5591
278,15	3,4651	-16,26	2,9242
283,15	3,0736	-15,13	3,2966
288,15	2,7603	-14,04	3,6708
293,15	2,5071	-12,99	4,0415
298,15	2,3009	-11,97	4,4038
303,15	2,1317	-10,98	4,7533
308,15	1,9923	-10,03	5,0859
313,15	1,8769	-9,11	5,3985
318,15	1,7813	-8,21	5,6882
323,15	1,7021	-7,35	5,9531
328,15	1,6365	-6,51	6,1917

Pressurisation des boissons gazeuses

En ouvrant une bouteille de champagne la pression chute brutalement dans la bouteille et, selon la loi de Henry, la solubilité du gaz diminue. Le dioxyde de carbone (Modèle:CO2) dissout dans le vin se désorbe et forme des bulles dans le liquide.

Les boissons gazeuses sont maintenues sous une pression importante de dioxyde de carbone (Modèle:CO2). On suppose qu'il n'y a que du Modèle:CO2 dans le gaz, soit $x_{{\text{CO}}_2}^{\text{g}}=1$, la présence d'eau en phase gaz étant négligeable, soit $x_{{\text{H}}_2\text{O}}^{\text{g}}\approx 0$. La pression partielle du Modèle:CO2 vaut alors la pression totale : $P_{{\text{CO}}_2}=x_{{\text{CO}}_2}^{\text{g}}P=P$. Cette pression permet de dissoudre une grande quantité de gaz dans le liquide en vertu de la loi de Henry :

$P_{{\text{CO}}_2}=x_{{\text{CO}}_2}^{\text{l}}{k}_{{\text{H,CO}}_2,\text{boisson}}$

En ouvrant la bouteille la pression $P_{{\text{CO}}_2}$ chute brusquement. Puisque $k_{{\text{H,CO}}_2,\text{boisson}}$ est quasi constante, alors $x_{{\text{CO}}_2}^{\text{l}}$ ne peut que diminuer : la solubilité du dioxyde de carbone chute. En conséquence, le Modèle:CO2 dissout se désorbe en formant des bulles dans le liquideModèle:,.

Solubilité des gaz atmosphériques dans le sang

La composition de l'atmosphère est considérée comme constante : elle contient environ 20,9 % molaire d'oxygène Modèle:Fchim et d'autres gaz dont l'azote Modèle:Fchim très majoritaire, soit les fractions molaires Modèle:Nobr et Modèle:Nobr. La loi de Henry donne les relations entre la pression atmosphérique totale, $P$, les pressions partielles d'oxygène et d'azote, $P_{{\text{O}}_2}$ et $P_{{\text{N}}_2}$, et les teneurs en oxygène et azote dissouts dans le sang, $x_{{\text{O}}_2}^{\text{l}}$ et $x_{{\text{N}}_2}^{\text{l}}$ :

$x_{{\text{O}}_2}^{\text{g}}P=P_{{\text{O}}_2}=x_{{\text{O}}_2}^{\text{l}}{k}_{{\text{H,O}}_2,\text{sang}}$

$x_{{\text{N}}_2}^{\text{g}}P=P_{{\text{N}}_2}=x_{{\text{N}}_2}^{\text{l}}{k}_{{\text{H,N}}_2,\text{sang}}$

En altitude, par exemple en montagne, la pression atmosphérique $P$ est plus basse qu'au niveau de la mer. La pression partielle en oxygène $P_{{\text{O}}_2}$ y est donc plus faible. En conséquence de la loi de Henry, la teneur en oxygène dissout dans le sang, $x_{{\text{O}}_2}^{\text{l}}$, est plus basse en montagne qu'à des altitudes moins élevées. Cet état est appelé hypoxie et peut provoquer le mal aigu des montagnes si ces altitudes sont atteintes trop rapidement, sans acclimatation progressive (par exemple à l'atterrissage en montagne d'un avion parti du niveau de la mer).

Au contraire, sous l'eau la pression est plus importante qu'au niveau de la surface. En un point où la pression $P$ est le double de la pression atmosphérique standard (c'est le cas à Modèle:Unité de profondeur), les solubilités des gaz, $x_{{\text{O}}_2}^{\text{l}}$ et $x_{{\text{N}}_2}^{\text{l}}$, sont doublées par rapport à la surface. Un plongeur consomme l'oxygène mais stocke l'azote de l'air dissout dans son organisme. Lorsque le plongeur remonte vers la surface, l'azote se désorbe en raison de la baisse de pression. Un accident de décompression survient si la remontée est trop rapide : le plongeur ne peut évacuer ce gaz par la respiration et l'azote forme des bulles dans le sang. Les bulles ainsi créées se dilatent dans les vaisseaux sanguins, toujours en raison de la baisse de pression (loi de Boyle-Mariotte), et peuvent provoquer une embolie gazeuse et le décès du plongeurModèle:,.

Constantes de Henry pour divers gaz dissouts dans l'eau

Le tableau suivant donne quelques valeurs de la constante de Henry pour des gaz dissouts dans l'eau à Modèle:Tmp (Modèle:Unité).

Valeurs de la constante de Henry
pour des gaz dissouts dans l'eau à Modèle:Unité.

Équation	$k_{\mathrm{H}}=\frac{P}{c_{\text{aq}}}$	$H=\frac{c_{\text{aq}}}{P}$	$k_{\mathrm{H}}=\frac{P}{x}$	$H=\frac{c_{\text{aq}}}{c_{\text{gas}}}$
Unité	Modèle:Unité	Modèle:Unité	Modèle:Unité	sans dimension
Oxygène (Modèle:Fchim)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre
Hydrogène (Modèle:Fchim)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre
Dioxyde de carbone (Modèle:Fchim)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre
Azote (Modèle:Fchim)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre
Hélium (He)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre
Néon (Ne)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre
Argon (Ar)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre
Monoxyde de carbone (CO)	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre	Modèle:Nombre

Avec : Modèle:Début de colonnes

• $c_{\text{aq}}$ la concentration du gaz en solution aqueuse (mol/l) ;
• $c_{\text{gaz}}$ la concentration du gaz en phase vapeur (mol/l) ;
• $P$ la pression partielle du gaz en phase vapeur (atm) ;
• $x$ la fraction molaire du gaz en solution aqueuse (sans dimension).

Notations

Alphabet latin

• $a,b,c,d,A,B,C,D$ les paramètres de formules empiriques ;
• $f_i^{\text{g}}$ la fugacité du Modèle:Nobr en phase vapeur ;
• $f_i^{\text{l}}$ la fugacité du Modèle:Nobr en phase liquide ;
• $f_i^{\text{l},*}$ la fugacité du Modèle:Nobr liquide pur ;
• ${\overline{H}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ l'enthalpie molaire partielle du Modèle:Nobr à dilution infinie dans le Modèle:Nobr ;
• ${\overline{H}}_{\sigma }^{\bullet ,*}$ l'enthalpie molaire du Modèle:Nobr à l'état de gaz parfait pur ;
• $n^{\text{l}}$ la quantité de matière totale en phase liquide ;
• $n_i^{\text{l}}$ la quantité du Modèle:Nobr en phase liquide ;
• $k_{\text{H},\sigma ,s}$ la constante de Henry du Modèle:Nobr dans le Modèle:Nobr ;
• $k_{\text{H},\sigma ,s}^{\text{sat}}$ la constante de Henry du Modèle:Nobr dans le Modèle:Nobr à $P_s^{\text{sat}}$ ;
• $k_{\text{H},\sigma ,m}$ la constante de Henry du Modèle:Nobr dans un mélange de solvants liquide ;
• $P$ la pression totale ;
• $P_i$ la pression partielle du Modèle:Nobr ;
• $P_s^{\text{sat}}$ la pression de vapeur saturante du Modèle:Nobr ;
• ${𝒫}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ le facteur de Poynting du Modèle:Nobr à dilution infinie dans le Modèle:Nobr liquide ;
• ${𝒫}_s^{\text{l},*}$ le facteur de Poynting du Modèle:Nobr liquide pur ;
• $R$ la constante universelle des gaz parfaits ;
• $T$ la température thermodynamique ;
• ${\overline{V}}_{\sigma }^{\text{g},*}$ le volume molaire du Modèle:Nobr gazeux pur ;
• $V^{\text{l}}$ le volume de la solution liquide ;
• ${\overline{V}}^{\text{l}}$ le volume molaire de la solution liquide ;
• ${\overline{V}}_i^{\text{l}}$ le volume molaire partiel du Modèle:Nobr dans le mélange liquide ;
• ${\overline{V}}_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ le volume molaire partiel du Modèle:Nobr à dilution infinie dans le Modèle:Nobr liquide ;
• ${\overline{V}}_s^{\text{l},*}$ le volume molaire du Modèle:Nobr liquide pur ;
• $x_i^{\text{g}}$ la fraction molaire du Modèle:Nobr en phase gaz ;
• $x_i^{\text{l}}$ la fraction molaire du Modèle:Nobr en phase liquide ;
• $y_i^{\text{l}}$ la quantité du Modèle:Nobr rapportée à la quantité totale de solvants dans un mélange liquide multicomposant.

Alphabet grec

• ${\gamma }_i^{\text{l}}$ le coefficient d'activité du Modèle:Nobr en phase liquide ;
• ${\gamma }_{\sigma ,m}^{\text{l},\mathrm{\infty }}$ le coefficient d'activité du Modèle:Nobr à dilution infinie dans un mélange de solvants liquide ;
• ${\gamma }_{\sigma ,s}^{\text{l},\mathrm{\infty }}$ le coefficient d'activité du Modèle:Nobr à dilution infinie dans le Modèle:Nobr liquide ;
• ${\hat{\gamma }}_s^{\text{l}}$ le coefficient d'ajustement de la loi de Raoult pour le Modèle:Nobr dans un mélange multicomposant liquide ;
• ${\mathrm{\Delta }}_{\text{sol}}{H}_{\sigma ,s}$ l'enthalpie de dissolution du Modèle:Nobr dans le Modèle:Nobr ;
• ${\varphi }_i^{\text{g}}$ le coefficient de fugacité du Modèle:Nobr en phase gaz ;
• ${\varphi }_{\sigma }^{\text{l}}$ le coefficient de fugacité du Modèle:Nobr en phase liquide ;
• ${\varphi }_s^{\text{g,*,sat}}$ le coefficient de fugacité du Modèle:Nobr pur en phase gaz à $P_s^{\text{sat}}$.

Notes et références

Bibliographie

Articles

• Modèle:Article.
• Modèle:Article.
• Modèle:Article.

Ouvrages

• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.

Articles connexes

• Activité chimique
• Équilibre de phases
• Fugacité
• Relation de Duhem-Margules
• Solution idéale