# 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 que$^,$ :
   « À 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. »
En pratique, elle ne s'applique qu'aux faibles concentrations du soluté (fraction molaire inférieure à 0,05) et à des pressions de moins de 10 bar (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 soluté $\sigma$ dissout dans un solvant $s$. 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 selon$^,$ :
* Loi de Henry pression partielle du soluté $\sigma$ dans le solvant $s$ : $P_{\sigma }=x_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}k_{{\text{H}},\sigma ,s}$
avec :
* $P$ la pression totale du mélange ;
* $P_{\sigma }$ la pression partielle du soluté $\sigma$, par définition $P_{\sigma }=x_{\sigma }^{\text{g}}P$ ;
* $k_{{\text{H}},\sigma ,s}$ la constante de Henry du soluté $\sigma$ dans le solvant $s$, 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$, ${\mathcal {H}}$ voire $K$ dans la littérature ;
* $x_{\sigma }^{\text{g}}$ la fraction molaire du soluté $\sigma$ dans la phase vapeur ;
* $x_{\sigma }^{\text{l}}$ la fraction molaire du soluté $\sigma$ 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={1 \over 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 que$^,$ :
   « À 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. »
et s'écrit sous la forme :
$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 à 20 °C l'oxygène et l'azote, de $T_{\text{c}}$ respectives −118,67 °C et −147,1 °C), 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 à 20 °C et 1 bar le propane de $P^{\text{sat}}$ = 8,327 bar). 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 à 20 °C et 1 bar l'éthanol de $P^{\text{sat}}$ = 5,8 × 10$^{−2}$ bar), 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 soluté $\sigma$ 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 solvant $s$ 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érature, 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 soluté $\sigma$ dans le solvant $s$ 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 soluté $\sigma$ et un unique solvant $s$.
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 soluté $\sigma$ et du solvant $s$ ; ceci implique qu'elle doit être déterminée pour chaque couple « soluté $\sigma$ - solvant $s$ » 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 soluté $\sigma$ avec un solvant autre que le solvant $s$).
En pratique, la constante de Henry est déterminée expérimentalement.

#### Définition thermodynamique

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 corps $\sigma$ dans le mélange :
* à dilution infinie : $\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}f_{\sigma }^{\text{l}}=0$ ;
* pour le corps pur : $\lim _{x_{\sigma }^{\text{l}}\to 1 \atop x_{s}^{\text{l}}\to 0}f_{\sigma }^{\text{l}}=f_{\sigma }^{{\text{l}},*}$ ;
avec $f_{\sigma }^{{\text{l}},*}$ la fugacité du corps $\sigma$ à l'état de liquide pur. Cette fugacité $f_{\sigma }^{{\text{l}},*}$ peut être fictive si le corps $\sigma$ 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éaires$^,$ :
* 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 soluté $\sigma$ à dilution infinie dans le solvant $s$. 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 soluté $\sigma$ dissout en phase liquide s'annule$^,$$^,$$^,$ :
* Constante de Henry : $k_{{\text{H}},\sigma ,s}\equiv \lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}{f_{\sigma }^{\text{l}} \over x_{\sigma }^{\text{l}}}$ à pression et température constantes.
avec :
* $f_{\sigma }^{\text{l}}$ la fugacité du soluté $\sigma$ dans le mélange liquide ;
* $k_{{\text{H}},\sigma ,s}$ la constante de Henry du soluté $\sigma$ dans le solvant $s$, aux pression $P$ et température $T$ du mélange ;
* $x_{\sigma }^{\text{l}}$ la fraction molaire du soluté $\sigma$ dans le mélange liquide ;
* $x_{s}^{\text{l}}$ la fraction molaire du solvant $s$ 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 par$^,$$^,$ :
$k_{{\text{H}},\sigma ,s}=\left({\partial f_{\sigma }^{\text{l}} \over \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\leq x_{\sigma }^{\text{l}}\leq 1$ du soluté $\sigma$, 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édemment$^,$ :
$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é$^,$ :
* à dilution infinie : $\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\gamma _{{\text{H}},\sigma }^{\text{l}}=1$ ;
* pour le corps pur : $\lim _{x_{\sigma }^{\text{l}}\to 1 \atop x_{s}^{\text{l}}\to 0}\gamma _{\sigma }^{\text{l}}=1$.
On pose à dilution infinie$^,$ :
$\gamma _{\sigma ,s}^{{\text{l}},\infty }\equiv \lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\gamma _{\sigma }^{\text{l}}$
Par conséquent, à dilution infinie on a :
    $\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}{f_{\sigma }^{\text{l}} \over x_{\sigma }^{\text{l}}}=\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\left(\gamma _{{\text{H}},\sigma }^{\text{l}}k_{{\text{H}},\sigma ,s}\right)=\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\left(\gamma _{\sigma }^{\text{l}}f_{\sigma }^{{\text{l}},*}\right)$
d'où la relation$^,$ (puisque $k_{{\text{H}},\sigma ,s}$ et $f_{\sigma }^{{\text{l}},*}$ ne dépendent pas de la composition) :
$k_{{\text{H}},\sigma ,s}=\gamma _{\sigma ,s}^{{\text{l}},\infty }f_{\sigma }^{{\text{l}},*}$
En injectant la relation précédente dans les expressions de $f_{\sigma }^{\text{l}}$ on obtient$^,$ :
$\gamma _{{\text{H}},\sigma }^{\text{l}}={\gamma _{\sigma }^{\text{l}} \over \gamma _{\sigma ,s}^{{\text{l}},\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 soluté $\sigma$ 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, Van Laar (en), Wilson, NRTL (en), UNIQUAC, UNIFAC, COSMOSPACE, etc.).

#### Dépendance à la pression

La fugacité $f_{\sigma }^{\text{l}}$ du soluté $\sigma$ dans le mélange liquide varie en fonction de la pression selon :
    $\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n^{\text{l}}}={{\bar {V}}_{\sigma }^{\text{l}} \over RT}$
avec :
* $V^{\text{l}}$ le volume de la phase liquide ;
* ${\bar {V}}_{\sigma }^{\text{l}}=\left({\partial V^{\text{l}} \over \partial n_{\sigma }^{\text{l}}}\right)_{P,T,n_{s}^{\text{l}}}$ le volume molaire partiel du soluté $\sigma$ dans le mélange liquide ;
* $n_{\sigma }^{\text{l}}$ la quantité du soluté $\sigma$ dans le mélange liquide ;
* $n_{s}^{\text{l}}$ la quantité du solvant $s$ dans le mélange liquide.
Quelle que soit la fraction molaire $x_{\sigma }^{\text{l}}=n_{\sigma }^{\text{l}}/\left(n_{\sigma }^{\text{l}}+n_{s}^{\text{l}}\right)$ du soluté $\sigma$, la dérivée partielle étant effectuée à composition constante, on peut écrire :
    $\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n^{\text{l}}}=\left({\partial \left[\ln f_{\sigma }^{\text{l}}-\ln x_{\sigma }^{\text{l}}+\ln x_{\sigma }^{\text{l}}\right] \over \partial P}\right)_{T,n^{\text{l}}}=\left({\partial \ln {f_{\sigma }^{\text{l}} \over x_{\sigma }^{\text{l}}} \over \partial P}\right)_{T,n^{\text{l}}}+\underbrace {\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n^{\text{l}}}} _{=\,0\,{\text{à composition constante}}}$
En passant à la limite de la dilution infinie :
    $\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n^{\text{l}}}=\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\left({\partial \ln {f_{\sigma }^{\text{l}} \over x_{\sigma }^{\text{l}}} \over \partial P}\right)_{T,n^{\text{l}}}=\left({\partial \ln k_{{\text{H}},\sigma ,s} \over \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 :
${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\equiv \lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}{\bar {V}}_{\sigma }^{\text{l}}$
La constante de Henry dépend par conséquent de la pression selon$^,$ :
* Dépendance de la constante de Henry à la pression $\left({\partial \ln k_{{\text{H}},\sigma ,s} \over \partial P}\right)_{T}={{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}$
avec :
* $P$ la pression ;
* $T$ la température ;
* ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$ le volume molaire partiel du soluté $\sigma$ à dilution infinie dans le solvant $s$ ;
* $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}\!\left(P,T\right)-\ln k_{{\text{H}},\sigma ,s}\!\left(P^{\circ },T\right)={\int _{P^{\circ }}^{P}{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\,\mathrm {d} P \over RT}$
La pression de référence $P^{\circ }$ est le plus souvent prise égale à la pression de vapeur saturante du solvant $s$ à 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}\!\left(P^{\circ },T\right)=k_{{\text{H}},\sigma ,s}\!\left(P_{s}^{\text{sat}}\!\left(T\right),T\right)=k_{{\text{H}},\sigma ,s}^{\text{sat}}\!\left(T\right)$. La constante de Henry est alors exprimée sous la forme$^,$$^,$ :
$k_{{\text{H}},\sigma ,s}\!\left(P,T\right)=k_{{\text{H}},\sigma ,s}^{\text{sat}}\!\left(T\right)\,{\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }$
avec le facteur de Poynting$^,$ :
${\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }\equiv \exp \!\left({\int _{P_{s}^{\text{sat}}}^{P}{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\,\mathrm {d} P \over RT}\right)$
Le volume molaire partiel ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$ représente la variation de volume de la solution liquide due à la dissolution d'une mole de soluté $\sigma$ dans une quantité infinie de solvant $s$. Il peut être déterminé expérimentalement par extrapolation de ${\bar {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 ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$ peut être considéré comme ne dépendant pas de la pression, soit ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\approx {\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\!\left(T\right)$, on obtient :
    $\ln k_{{\text{H}},\sigma ,s}=\ln k_{{\text{H}},\sigma ,s}^{\text{sat}}+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over 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 ${\bar {V}}_{\sigma ,s}^{{\text{l}},\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 soluté $\sigma$ dans le mélange liquide varie en fonction de la température selon :
    $\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial {1 \over T}}\right)_{P,n^{\text{l}}}={{\bar {H}}_{\sigma }^{\text{l}}-{\bar {H}}_{\sigma }^{\bullet ,*} \over R}$
avec :
* $H^{\text{l}}$ l'enthalpie de la phase liquide ;
* ${\bar {H}}_{\sigma }^{\text{l}}=\left({\partial H^{\text{l}} \over \partial n_{\sigma }^{\text{l}}}\right)_{P,T,n_{s}^{\text{l}}}$ l'enthalpie molaire partielle du soluté $\sigma$ dans le mélange liquide ;
* ${\bar {H}}_{\sigma }^{\bullet ,*}$ l'enthalpie molaire du soluté $\sigma$ à l'état de gaz parfait pur à $T$ ;
* $n_{\sigma }^{\text{l}}$ la quantité du soluté $\sigma$ dans le mélange liquide ;
* $n_{s}^{\text{l}}$ la quantité du solvant $s$ dans le mélange liquide.
Quelle que soit la fraction molaire $x_{\sigma }^{\text{l}}=n_{\sigma }^{\text{l}}/\left(n_{\sigma }^{\text{l}}+n_{s}^{\text{l}}\right)$ du soluté $\sigma$, la dérivée partielle étant effectuée à composition constante, on peut écrire :
    $\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial {1 \over T}}\right)_{P,n^{\text{l}}}=\left({\partial \left[\ln f_{\sigma }^{\text{l}}-\ln x_{\sigma }^{\text{l}}+\ln x_{\sigma }^{\text{l}}\right] \over \partial {1 \over T}}\right)_{P,n^{\text{l}}}=\left({\partial \ln {f_{\sigma }^{\text{l}} \over x_{\sigma }^{\text{l}}} \over \partial {1 \over T}}\right)_{P,n^{\text{l}}}+\underbrace {\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial {1 \over T}}\right)_{P,n^{\text{l}}}} _{=\,0\,{\text{à composition constante}}}$
En passant à la limite de la dilution infinie :
    $\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial {1 \over T}}\right)_{P,n^{\text{l}}}=\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\left({\partial \ln {f_{\sigma }^{\text{l}} \over x_{\sigma }^{\text{l}}} \over \partial {1 \over T}}\right)_{P,n^{\text{l}}}=\left({\partial \ln k_{{\text{H}},\sigma ,s} \over \partial {1 \over 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 partielle$^,$ :
${\bar {H}}_{\sigma ,s}^{{\text{l}},\infty }\equiv \lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}{\bar {H}}_{\sigma }^{\text{l}}$
c'est-à-dire l'enthalpie molaire partielle du soluté $\sigma$ à dilution infinie dans le solvant $s$ liquide. L'enthalpie molaire du gaz parfait pur ${\bar {H}}_{\sigma }^{\bullet ,*}$ ne dépendant pas de la composition, elle reste inchangée lors du passage à la limite et on obtient :
    $\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\left({\bar {H}}_{\sigma }^{\text{l}}-{\bar {H}}_{\sigma }^{\bullet ,*}\right)={\bar {H}}_{\sigma ,s}^{{\text{l}},\infty }-{\bar {H}}_{\sigma }^{\bullet ,*}$
On pose pour les enthalpies molaires partielles$^,$ :
$\Delta _{\text{sol}}H_{\sigma ,s}\equiv {\bar {H}}_{\sigma ,s}^{{\text{l}},\infty }-{\bar {H}}_{\sigma }^{\bullet ,*}$
avec :
* ${\bar {H}}_{\sigma ,s}^{{\text{l}},\infty }$ l'enthalpie molaire partielle du soluté $\sigma$ à dilution infinie dans le solvant $s$ liquide à $T$ ;
* ${\bar {H}}_{\sigma }^{\bullet ,*}$ l'enthalpie molaire du soluté $\sigma$ à l'état de gaz parfait pur à $T$.
La constante de Henry dépend par conséquent de la température selon$^,$ :
* Dépendance de la constante de Henry à la température $\left({\partial \ln k_{{\text{H}},\sigma ,s} \over \partial {1 \over T}}\right)_{P}={\Delta _{\text{sol}}H_{\sigma ,s} \over R}$
avec :
* $P$ la pression ;
* $T$ la température ;
* $\Delta _{\text{sol}}H_{\sigma ,s}$ l'enthalpie de dissolution du soluté $\sigma$ dans le solvant $s$ à $T$ ;
* $R$ la constante universelle des gaz parfaits.
Si l'on considère l'enthalpie de dissolution $\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}\!\left(P,T\right)-\ln k_{{\text{H}},\sigma ,s}\!\left(P,T^{\circ }\right)={\Delta _{\text{sol}}H_{\sigma ,s} \over R}\left({1 \over T}-{1 \over T^{\circ }}\right)$
    $k_{{\text{H}},\sigma ,s}\!\left(P,T\right)=k_{{\text{H}},\sigma ,s}\!\left(P,T^{\circ }\right)\cdot \exp \!\left({\Delta _{\text{sol}}H_{\sigma ,s} \over R}\left({1 \over T}-{1 \over T^{\circ }}\right)\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 « soluté $\sigma$ - solvant $s$ » :
    $\ln k_{{\text{H}},\sigma ,s}\!\left(P,T\right)=A+{B \over T}$
La littérature utilise parfois l'inverse de la constante de Henry définie précédemment, $H_{\sigma ,s}={1 \over k_{{\text{H}},\sigma ,s}}$ (cette notation prête à confusion avec celle de l'enthalpie de dissolution $\Delta _{\text{sol}}H_{\sigma ,s}$), aussi trouve-t-on également les relations :
    $\left({\partial \ln H_{\sigma ,s} \over \partial {1 \over T}}\right)_{P}=-{\Delta _{\text{sol}}H_{\sigma ,s} \over R}$
    $\ln H_{\sigma ,s}\!\left(P,T\right)-\ln H_{\sigma ,s}\!\left(P,T^{\circ }\right)=-{\Delta _{\text{sol}}H_{\sigma ,s} \over R}\left({1 \over T}-{1 \over T^{\circ }}\right)$
    $H_{\sigma ,s}\!\left(P,T\right)=H_{\sigma ,s}\!\left(P,T^{\circ }\right)\cdot \exp \!\left(-{\Delta _{\text{sol}}H_{\sigma ,s} \over R}\left({1 \over T}-{1 \over T^{\circ }}\right)\right)$
L'enthalpie de dissolution $\Delta _{\text{sol}}H_{\sigma ,s}$ est la chaleur produite par la dissolution d'une mole de soluté $\sigma$ à l'état de gaz parfait pur dans une quantité infinie de solvant $s$ à 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 ${\bar {H}}_{\sigma ,s}^{{\text{l}},\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 ${\bar {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 : $\Delta _{\text{sol}}H_{\sigma ,s}=\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 $\Delta _{\text{sol}}H_{\sigma ,s}<0$, par conséquent $\ln k_{{\text{H}},\sigma ,s}$ augmente avec une diminution de ${1 \over 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é $\phi _{\sigma }^{\text{l}}$ du soluté $\sigma$ en phase liquide est défini par la relation avec la fugacité $f_{\sigma }^{\text{l}}$ :
$\phi _{\sigma }^{\text{l}}\equiv {f_{\sigma }^{\text{l}} \over x_{\sigma }^{\text{l}}P}$
Par définition de la constante de Henry, on a donc$^,$ :
${k_{{\text{H}},\sigma ,s} \over P}=\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\phi _{\sigma }^{\text{l}}$
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\!\left(V^{\text{l}},T,n^{\text{l}}\right)$, le coefficient de fugacité $\phi _{\sigma }^{\text{l}}$ se calcule selon :
    $RT\,\ln \phi _{\sigma }^{\text{l}}=-\int _{+\infty }^{V^{\text{l}}}\left[\left({\partial P \over \partial n_{\sigma }^{\text{l}}}\right)_{V,T,n_{s}^{\text{l}}}-{RT \over V}\right]\,\mathrm {d} V-RT\,\ln \!\left({PV^{\text{l}} \over \left(n_{\sigma }^{\text{l}}+n_{s}^{\text{l}}\right)RT}\right)$
avec $n_{\sigma }^{\text{l}}$ et $n_{s}^{\text{l}}$ les quantités respectives du soluté $\sigma$ et du solvant $s$ 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 solvant $s$ et du soluté $\sigma$, 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={RT \over {\bar {V}}^{\text{l}}-b_{m}^{\text{l}}}-{a_{m}^{\text{l}} \over {{\bar {V}}^{\text{l}}}^{2}}$
    $RT\,\ln \phi _{\sigma }^{\text{l}}=-{a_{m}^{\text{l}} \over {\bar {V}}^{\text{l}}}\left[\delta _{\sigma }^{\text{l}}-{b_{\sigma } \over b_{m}^{\text{l}}}\right]+{b_{\sigma } \over b_{m}^{\text{l}}}\left(P{\bar {V}}^{\text{l}}-RT\right)-RT\,\ln \!\left({P\!\cdot \!\left({\bar {V}}^{\text{l}}-b_{m}^{\text{l}}\right) \over RT}\right)$
avec :
    $a_{m}^{\text{l}}={x_{\sigma }^{\text{l}}}^{2}a_{\sigma }+2\,x_{\sigma }^{\text{l}}x_{s}^{\text{l}}{\sqrt {a_{\sigma }a_{s}}}\left(1-k_{\sigma ,s}\right)+{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{{\sqrt {a_{\sigma }}} \over a_{m}^{\text{l}}}\left(x_{\sigma }^{\text{l}}{\sqrt {a_{\sigma }}}+x_{s}^{\text{l}}{\sqrt {a_{s}}}\left(1-k_{\sigma ,s}\right)\right)$
À 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}^{\infty }=2\,{\sqrt {a_{\sigma } \over a_{s}}}\left(1-k_{\sigma ,s}\right)$
Le volume molaire ${\bar {V}}^{\text{l}}$ de la solution liquide tend vers celui ${\bar {V}}_{s}^{\text{l,*}}$ du solvant $s$ liquide pur à $P$ et $T$, calculable par l'équation d'état : ${\bar {V}}^{\text{l}}\to {\bar {V}}_{s}^{\text{l,*}}$. On obtient $k_{{\text{H}},\sigma ,s}$, la constante de Henry à $P$ et $T$ :
$RT\,\ln {k_{{\text{H}},\sigma ,s} \over P}=-{a_{s} \over {\bar {V}}_{s}^{\text{l,*}}}\left[\delta _{\sigma ,s}^{\infty }-{b_{\sigma } \over b_{s}}\right]+{b_{\sigma } \over b_{s}}\left(P{\bar {V}}_{s}^{\text{l,*}}-RT\right)-RT\,\ln \!\left({P\!\cdot \!\left({\bar {V}}_{s}^{\text{l,*}}-b_{s}\right) \over RT}\right)$
Lorsque la pression $P$ tend vers la pression de vapeur saturante $P_{s}^{\text{sat}}$ du solvant $s$ à la température $T$, l'équilibre liquide-vapeur tend vers celui du solvant $s$ pur (état de saturation) ; le volume molaire ${\bar {V}}_{s}^{\text{l,*}}$ tend vers celui ${\bar {V}}_{s}^{\text{l,*,sat}}$ du solvant $s$ 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 solvant $s$ :
$RT\,\ln {k_{{\text{H}},\sigma ,s}^{\text{sat}} \over P_{s}^{\text{sat}}}=-{a_{s} \over {\bar {V}}_{s}^{\text{l,*,sat}}}\left[\delta _{\sigma ,s}^{\infty }-{b_{\sigma } \over b_{s}}\right]+{b_{\sigma } \over b_{s}}\left(P_{s}^{\text{sat}}{\bar {V}}_{s}^{\text{l,*,sat}}-RT\right)-RT\,\ln \!\left({P_{s}^{\text{sat}}\!\cdot \!\left({\bar {V}}_{s}^{\text{l,*,sat}}-b_{s}\right) \over RT}\right)$
Les liquides étant peu compressibles, le volume molaire ${\bar {V}}_{s}^{\text{l,*}}$ du solvant $s$ liquide pur à $P$ est assimilable à celui ${\bar {V}}_{s}^{\text{l,*,sat}}$ du solvant $s$ liquide pur à saturation : ${\bar {V}}_{s}^{\text{l,*}}\approx {\bar {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}}+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over RT}$
on identifie ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$, le volume molaire partiel du soluté $\sigma$ à dilution infinie dans le solvant $s$ :
${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\approx {b_{\sigma } \over b_{s}}{\bar {V}}_{s}^{\text{l,*,sat}}$
Avec :
* $a_{\sigma }={27R^{2}{T_{{\text{c}},\sigma }}^{2} \over 64P_{{\text{c}},\sigma }}$ et $a_{s}={27R^{2}{T_{{\text{c}},s}}^{2} \over 64P_{{\text{c}},s}}$ ;
* $b_{\sigma }={R{T_{{\text{c}},\sigma }} \over 8P_{{\text{c}},\sigma }}$ et $b_{s}={R{T_{{\text{c}},s}} \over 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 soluté $\sigma$ et du solvant $s$ ;
* $P_{s}^{\text{sat}}$ la pression de vapeur saturante du solvant $s$ à la température $T$ ;
* $T_{{\text{c}},\sigma }$ et $T_{{\text{c}},s}$ les températures critiques respectives du soluté $\sigma$ et du solvant $s$ ;
* ${\bar {V}}^{\text{l}}$ le volume molaire de la solution liquide à la pression $P$ et à la température $T$ ;
* ${\bar {V}}_{s}^{\text{l,*}}$ le volume molaire du solvant $s$ liquide pur à $P$ et $T$ ;
* ${\bar {V}}_{s}^{\text{l,*,sat}}$ le volume molaire du solvant $s$ 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 solvant $s$ et du soluté $\sigma$, et aux interactions entre le soluté $\sigma$ et le solvant $s$ par l'intermédiaire du coefficient d'interaction binaire $k_{\sigma ,s}$. La constante de Henry est donc spécifique du couple « soluté $\sigma$ - solvant $s$ » et ne peut être utilisée pour d'autres mélanges binaires ; par exemple, elle n'est pas valable si le soluté $\sigma$ 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 :
$k_{{\text{H}},\sigma ,s}\!\left(P,T\right)=k_{{\text{H}},\sigma ,s}^{\text{sat}}\!\left(T\right)\cdot \exp \!\left({{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over 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 : ${\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }\approx 1$. La constante de Henry peut alors être considérée comme indépendante de la pression et être approchée par :
$k_{{\text{H}},\sigma ,s}\!\left(P,T\right)\approx k_{{\text{H}},\sigma ,s}^{\text{sat}}\!\left(T\right)$
Les formes suivantes sont souvent utilisées pour la dépendance à la température$^,$ :
    $\ln k_{{\text{H}},\sigma ,s}^{\text{sat}}\!\left(T\right)=A+{B \over T}+C\,\ln T+D\,T$
    $\ln k_{{\text{H}},\sigma ,s}^{\text{sat}}\!\left(T\right)=A+{B \over T}+{C \over T^{2}}+{D \over T^{3}}$
avec $A$, $B$, $C$ et $D$ des constantes empiriques spécifiques du couple « soluté $\sigma$ - solvant $s$ ». L'enthalpie de dissolution est alors exprimée sous les formes respectives :
    ${\Delta _{\text{sol}}H_{\sigma ,s} \over R}=B-C\,T-D\,T^{2}$
    ${\Delta _{\text{sol}}H_{\sigma ,s} \over R}=B+{2\,C \over T}+{3\,D \over T^{2}}$

### Démonstration de la loi de Henry

Lorsque l'équilibre liquide-vapeur est atteint, les fugacités du soluté $\sigma$ sont homogènes entre les deux phases :
$f_{\sigma }^{\text{g}}=f_{\sigma }^{\text{l}}$
avec :
* $f_{\sigma }^{\text{g}}$ la fugacité du soluté $\sigma$ en phase gaz (vapeur) ;
* $f_{\sigma }^{\text{l}}$ la fugacité du soluté $\sigma$ en phase liquide.
Aux basses pressions (moins de 10 bar), le gaz se comporte comme un mélange de gaz parfaits, et la fugacité du soluté $\sigma$ 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 soluté $\sigma$ dans le solvant $s$ 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 soluté $\sigma$ 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 10 bar), 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 Setchenov$^,$$^,$ :
$\ln {k_{{\text{H}},\sigma ,s} \over k_{{\text{H}},\sigma ,s}^{*}}=kI$
avec :
* $k_{{\text{H}},\sigma ,s}$ la constante de Henry du soluté $\sigma$ dans le solvant $s$ en solution avec tous les solutés ;
* $k_{{\text{H}},\sigma ,s}^{*}$ la constante de Henry du soluté $\sigma$ dans le solvant $s$ 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 CO$_2$, le dioxyde de soufre SO$_2$ et l'ammoniac NH$_3$ 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 (CO$_2$), qui forme partiellement, par réaction avec l'eau, de l'acide carbonique (H$_2$CO$_3$), qui lui-même, en fonction du pH de l'eau, forme les ions hydrogénocarbonate (HCO$_3$$^−$) et carbonate (CO$_3$$^{2−}$) selon les réactions :
* ${\rm {CO_{2}(g)\rightleftarrows CO_{2}(aq)}}$ (dissolution)
    ${\rm {[CO_{2}]}}=K_{1}\,P_{\rm {CO_{2}}}$ (loi de Henry)
* ${\rm {CO_{2}(aq)+H_{2}O\rightleftarrows H_{2}CO_{3}(aq)}}$
    ${\rm {[H_{2}CO_{3}]}}=K_{2}\,{\rm {[CO_{2}]}}$
* ${\rm {H_{2}CO_{3}(aq)+H_{2}O\rightleftarrows {HCO_{3}}^{-}(aq)+{H_{3}O}^{+}(aq)}}$
    ${\rm {[{HCO_{3}}^{-}]}}\,{\rm {[{H_{3}O}^{+}]}}=K_{3}\,{\rm {[H_{2}CO_{3}]}}$
* ${\rm {{HCO_{3}}^{-}(aq)+H_{2}O\rightleftarrows {CO_{3}}^{2-}(aq)+{H_{3}O}^{+}(aq)}}$
    ${\rm {[{CO_{3}}^{2-}]}}\,{\rm {[{H_{3}O}^{+}]}}=K_{4}\,{\rm {[{HCO_{3}}^{-}]}}$
avec :
* $K_{i}$ la constante de l'équilibre $i$ ; à 298 K : $K_{1}$ = 3,39 × 10$^{−2}$, $K_{2}$ = 1,7 × 10$^{−3}$, $K_{3}$ = 2,34 × 10$^{−4}$, $K_{4}$ = 5,01 × 10$^{−11}$ ;
* $P_{\rm {CO_{2}}}$ la pression partielle du CO$_2$ en phase gaz (en bar) ;
* ${\rm {[X]}}$ la concentration molaire de l'espèce ${\rm {X}}$ en solution aqueuse (en mol/l).
La loi de Henry à proprement parler ne s'applique qu'à la première des quatre réactions ci-dessus, celle de dissolution du CO$_2$ gazeux en CO$_2$ aqueux. Si l'on considère le CO$_2$ dissout sous toutes ses formes, le bilan de matière donne :
    ${\rm {[CO_{2}]_{\text{total}}}}={\rm {[CO_{2}]}}+{\rm {[H_{2}CO_{3}]}}+{\rm {[{HCO_{3}}^{-}]}}+{\rm {[{CO_{3}}^{2-}]}}$
La loi de Henry peut ainsi être modifiée selon :
    ${\rm {[CO_{2}]_{\text{total}}}}=\left(K_{1}+K_{1}K_{2}+{K_{1}K_{2}K_{3} \over {\rm {[{H_{3}O}^{+}]}}}+{K_{1}K_{2}K_{3}K_{4} \over {\rm {[{H_{3}O}^{+}]}}^{2}}\right)\,P_{\rm {CO_{2}}}$
À pression $P_{\rm {CO_{2}}}$ constante, lorsque ${\rm {[{H_{3}O}^{+}]}}$ diminue, ${\rm {[CO_{2}]_{\text{tot}}}}$ augmente ; donc plus la solution est basique, plus l'on peut dissoudre de CO$_2$ ; inversement, plus la solution est acide, moins l'on peut dissoudre de CO$_2$.

##### Remarque

    Puisque $K_{2}$ = 1,7 × 10$^{−3}$, on a ${\rm {[H_{2}CO_{3}]\ll [CO_{2}]}}$ quel que soit le pH. Les deuxième et troisième réactions sont regroupées en une seule réaction :
       ${\rm {CO_{2}(aq)+2\,H_{2}O\rightleftarrows {HCO_{3}}^{-}(aq)+{H_{3}O}^{+}(aq)}}$
    avec ${\rm {[{HCO_{3}}^{-}]}}\,{\rm {[{H_{3}O}^{+}]}}=K_{5}\,{\rm {[CO_{2}]}}$ et $K_{5}=K_{3}\,K_{2}$ = 3,98 × 10$^{−7}$. Le bilan de matière donne, en négligeant H$_2$CO$_3$ :
       ${\rm {[CO_{2}]_{\text{total}}}}=\left(K_{1}+{K_{1}K_{5} \over {\rm {[{H_{3}O}^{+}]}}}+{K_{1}K_{4}K_{5} \over {\rm {[{H_{3}O}^{+}]}}^{2}}\right)\,P_{\rm {CO_{2}}}$

#### Cas des solvants, loi de Raoult

Un solvant $s$ est un corps présent dans une solution liquide ayant une fraction molaire très supérieure à celle d'un soluté $\sigma$, 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 :
$P_{s}=x_{s}^{\text{g}}P=x_{s}^{\text{l}}P_{s}^{\text{sat}}{\mathcal {P}}_{s}^{{\text{l}},*}$
avec :
* $P$ la pression totale du mélange ;
* $P_{s}$ la pression partielle du solvant $s$, 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 ;
* ${\mathcal {P}}_{s}^{{\text{l}},*}=\exp \!\left({\int _{P_{s}^{\text{sat}}}^{P}{\bar {V}}_{s}^{{\text{l}},*}\,\mathrm {d} P \over RT}\right)$ le facteur de Poynting appliqué au solvant $s$ ;
* ${\bar {V}}_{s}^{{\text{l}},*}$ le volume molaire du solvant $s$ liquide pur ;
* $x_{s}^{\text{g}}$ la fraction molaire du solvant $s$ dans la phase vapeur ;
* $x_{s}^{\text{l}}$ la fraction molaire du solvant $s$ 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 ${\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }={\mathcal {P}}_{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 corps $i$, 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é $\phi _{i}^{\text{g}}$ corrigeant la loi des gaz parfaits :
    $f_{i}^{\text{g}}=x_{i}^{\text{g}}\phi _{i}^{\text{g}}P$
La fugacité $f_{i}^{\text{l}}$ en phase liquide s'écrit à l'aide d'un coefficient d'activité $\Gamma _{i}^{\text{l}}$ corrigeant la loi de Lewis et Randall :
    $f_{i}^{\text{l}}=x_{i}^{\text{l}}\Gamma _{i}^{\text{l}}f_{i}^{\text{l,*}}$
Ce coefficient est décomposé en deux termes :
    $\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({\partial G^{\text{E}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq 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 $\phi _{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 10 bar), le gaz se comporte comme un mélange de gaz parfaits, soit $\phi _{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 corps $i$ :
    é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 soluté $\sigma$ et un unique solvant $s$.

##### Lois des équilibres binaires

La fugacité (fictive si le soluté est un gaz) $f_{\sigma }^{\text{l,*}}$ du soluté $\sigma$ à 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,*}}={k_{{\text{H}},\sigma ,s} \over \gamma _{\sigma ,s}^{{\text{l}},\infty }}={k_{{\text{H}},\sigma ,s}^{\text{sat}} \over \gamma _{\sigma ,s}^{{\text{l}},\infty }}{\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }$
Pour le soluté $\sigma$, 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 soluté $\sigma$ : $x_{\sigma }^{\text{g}}\phi _{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}\gamma _{\sigma }^{\text{l}}{k_{{\text{H}},\sigma ,s}^{\text{sat}} \over \gamma _{\sigma ,s}^{{\text{l}},\infty }}{\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }$
La fugacité $f_{s}^{\text{l,*}}$ du solvant $s$ à 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,*}}=\phi _{s}^{\text{g,*,sat}}P_{s}^{\text{sat}}{\mathcal {P}}_{s}^{{\text{l}},*}$
Pour le solvant $s$, 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 solvant $s$ : $x_{s}^{\text{g}}\phi _{s}^{\text{g}}P=x_{s}^{\text{l}}\gamma _{s}^{\text{l}}\phi _{s}^{\text{g,*,sat}}P_{s}^{\text{sat}}{\mathcal {P}}_{s}^{{\text{l}},*}$
En application du théorème d'Euler, le volume molaire ${\bar {V}}^{\text{l}}$ de la phase liquide vaut :
    ${\begin{aligned}{{\bar {V}}^{\text{l}} \over RT}&={x_{\sigma }^{\text{l}}{\bar {V}}_{\sigma }^{\text{l}}+x_{s}^{\text{l}}{\bar {V}}_{s}^{\text{l}} \over RT}=x_{\sigma }^{\text{l}}\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}+x_{s}^{\text{l}}\left({\partial \ln f_{s}^{\text{l}} \over \partial P}\right)_{T,n}\\&=x_{\sigma }^{\text{l}}\left[{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}-\left({\partial \ln \gamma _{\sigma ,s}^{{\text{l}},\infty } \over \partial P}\right)_{T,n}+\left({\partial \ln \gamma _{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}\right]+x_{s}^{\text{l}}\left[{{\bar {V}}_{s}^{{\text{l}},*} \over RT}+\left({\partial \ln \gamma _{s}^{\text{l}} \over \partial P}\right)_{T,n}\right]\end{aligned}}$
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 :
$\sigma$
$x_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}k_{{\text{H}},\sigma ,s}^{\text{sat}}{\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }$
et :
$s$
$x_{s}^{\text{g}}P=x_{s}^{\text{l}}P_{s}^{\text{sat}}{\mathcal {P}}_{s}^{{\text{l}},*}$
Dans une solution idéale, les lois de Henry et de Raoult sont applicables sur l'ensemble des fractions molaires $0\leq x_{i}^{\text{l}}\leq 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 ${\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }\approx {\mathcal {P}}_{s}^{{\text{l}},*}\approx 1$.
Le volume molaire ${\bar {V}}^{\text{l,id}}$ de la phase liquide idéale vaut :
    ${\bar {V}}^{\text{l,id}}=x_{\sigma }^{\text{l}}{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }+x_{s}^{\text{l}}{\bar {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 ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$ dans le facteur de Poynting ${\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }$ ne dépend pas de la pression. On obtient, pour tout soluté $\sigma$ :
    $\ln {f_{\sigma }^{\text{g}} \over x_{\sigma }^{\text{l}}}=\ln \!\left(\gamma _{\sigma }^{\text{l}}{k_{{\text{H}},\sigma ,s}^{\text{sat}} \over \gamma _{\sigma ,s}^{{\text{l}},\infty }}{\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }\right)=\ln k_{{\text{H}},\sigma ,s}^{\text{sat}}+\ln \gamma _{\sigma }^{\text{l}}-\ln \gamma _{\sigma ,s}^{{\text{l}},\infty }+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over 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}},\infty }=\lim _{x_{\sigma }^{\text{l}}\to 0 \atop x_{s}^{\text{l}}\to 1}\ln \gamma _{\sigma }^{\text{l}}=A$
avec $x_{s}^{\text{l}}$ la fraction molaire du solvant $s$ dans la phase liquide (avec $x_{s}^{\text{l}}+x_{\sigma }^{\text{l}}=1$). On obtient l'équation de Krichevsky-Ilinskaya$^,$ :
$\ln {f_{\sigma }^{\text{g}} \over x_{\sigma }^{\text{l}}}=\ln k_{{\text{H}},\sigma ,s}^{\text{sat}}+A\cdot \left({x_{s}^{\text{l}}}^{2}-1\right)+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over 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 forme$^,$$^,$$^,$ :
$\ln {f_{\sigma }^{\text{g}} \over x_{\sigma }^{\text{l}}}=\ln k_{{\text{H}},\sigma ,s}^{\text{sat}}+{A \over RT}\cdot \left({x_{s}^{\text{l}}}^{2}-1\right)+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over RT}$
Si le mélange liquide est idéal, soit $\gamma _{\sigma }^{\text{l}}=\gamma _{\sigma ,s}^{{\text{l}},\infty }=1$ (d'où $A=0$), on obtient l'équation de Krichevsky–Kasarnovsky$^,$$^,$$^,$$^,$ :
$\ln {f_{\sigma }^{\text{g}} \over x_{\sigma }^{\text{l}}}=\ln k_{{\text{H}},\sigma ,s}^{\text{sat}}+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over 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 10 bar) le gaz se comporte comme un gaz parfait, la fugacité du soluté $\sigma$ 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 soluté $\sigma$ en phase gaz est calculée à l'aide d'un coefficient de fugacité : $f_{\sigma }^{\text{g}}=x_{\sigma }^{\text{g}}\phi _{\sigma }^{\text{g}}P$. Les équations de Krichevsky-Ilinskaya et Krichevsky-Kasarnovsky sont employées pour calculer des solubilités à haute pression, jusqu'à 1 000 bar environ$^,$.
Paramètres de l'équation de Krichevsky-Ilinskaya pour l'hydrogène H$_2$ dans divers solvants. $\ln {f_{\rm {H_{2}}}^{\text{g}} \over x_{\rm {H_{2}}}^{\text{l}}}=\ln k_{{\text{H}},{\rm {H_{2}}},s}^{\text{sat}}+{A \over RT}\cdot \left({x_{s}^{\text{l}}}^{2}-1\right)+{{\bar {V}}_{{\rm {H_{2}}},s}^{{\text{l}},\infty }\cdot \left(P-P_{s}^{\text{sat}}\right) \over RT}$ :
| Solvant $s$ | Température $T$ (K) | Constante de Henry $k_{{\text{H}},{\rm {H_{2}}},s}^{\text{sat}}$ (bar) (à $P_{s}^{\text{sat}}\!\left(T\right)$) | Paramètre de Margules $A$ (J mol$^{−1}$) | Volume molaire à dilution infinie ${\bar {V}}_{{\rm {H_{2}}},s}^{{\text{l}},\infty }$ (cm$^3$ mol$^{−1}$) |
| Monoxyde de carbone CO | 68 | 648 | 704±69 | 31,2 |
| Monoxyde de carbone CO | 78 | 476 | 704±69 | 32,6 |
| Monoxyde de carbone CO | 88 | 405 | 704±69 | 34,4 |
| Azote N$_2$ | 68 | 547 | 704±69 | 30,4 |
| Azote N$_2$ | 79 | 456 | 704±69 | 31,5 |
| Azote N$_2$ | 95 | 345 | 704±69 | 34,4 |
| Méthane CH$_4$ | 90 | 1848 | 1486±297 | 29,7 |
| Méthane CH$_4$ | 110 | 1050 | 1486±297 | 31,0 |
| Méthane CH$_4$ | 144 | 638 | 1486±297 | 36,0 |
| Éthane C$_2$H$_6$ | 144 | 2634 | 2478±198 | 37,9 |
| Éthane C$_2$H$_6$ | 200 | 1672 | 2478±198 | 44,2 |
| Éthane C$_2$H$_6$ | 228 | 1226 | 2478±198 | 54,3 |
| Propane C$_3$H$_8$ | 228 | 1692 | 2478±198 | 50 |
| Propane C$_3$H$_8$ | 255 | 1317 | 2478±198 | 51 |
| Propane C$_3$H$_8$ | 282 | 1044 | 2478±198 | 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 « soluté $\sigma$ - solvant $s$ ».
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 corps $i$, soluté ou solvant, le facteur :
$i$
$y_{i}^{\text{l}}={x_{i}^{\text{l}} \over \sum _{s}x_{s}^{\text{l}}}$
avec $x_{i}^{\text{l}}$ la fraction molaire du corps $i$ dans le mélange multicomposant liquide. $y_{\sigma }^{\text{l}}$ est la quantité de soluté $\sigma$ rapportée à la quantité totale de solvants dans le mélange liquide. $y_{s}^{\text{l}}$ est la fraction molaire du solvant $s$ dans le mélange liquide des solvants seuls, autrement dit la limite de la fraction molaire $x_{s}^{\text{l}}$ du solvant $s$ à 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 soluté $\sigma$ à l'état de liquide pur calculée à partir des propriétés du mélange binaire « soluté $\sigma$ - solvant $s$ » :
$\sigma$
$s$
$f_{\sigma ,s}^{\text{l,*}}={k_{{\text{H}},\sigma ,s} \over \gamma _{\sigma ,s}^{{\text{l}},\infty }}$
Si le soluté $\sigma$ 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 solvant $s$ à l'autre. La fugacité $f_{\sigma }^{\text{l,*}}$ de tout soluté $\sigma$ à l'état de liquide pur est construite selon l'hypothèse de la solution idéale et peut être fictive :
$\sigma$
$\ln f_{\sigma }^{\text{l,*}}=\sum _{s}y_{s}^{\text{l}}\,\ln f_{\sigma ,s}^{\text{l,*}}$
Pour tout soluté $\sigma$, 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 soluté $\sigma$ : $x_{\sigma }^{\text{g}}\phi _{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}\gamma _{\sigma }^{\text{l}}f_{\sigma }^{\text{l,*}}$
La fugacité $f_{s}^{\text{l,*}}$ d'un solvant $s$ à l'état de liquide pur dans les conditions de pression et température du mélange est donnée par la relation :
$s$
$f_{s}^{\text{l,*}}=\phi _{s}^{\text{g,*,sat}}P_{s}^{\text{sat}}{\mathcal {P}}_{s}^{{\text{l}},*}$
Pour tout solvant $s$ on pose le facteur correctif ${\hat {\gamma }}_{s}^{\text{l}}$ :
$s$
$\ln {\hat {\gamma }}_{s}^{\text{l}}=\sum _{\sigma }y_{\sigma }^{\text{l}}\left(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}}\right)$
En l'absence de tout soluté, soit $y_{\sigma }^{\text{l}}=0$ pour tout soluté $\sigma$, ou en présence d'un unique solvant $s$, soit $y_{s}^{\text{l}}=1$ pour le solvant et $f_{\sigma ,s}^{\text{l,*}}=f_{\sigma }^{\text{l,*}}$ pour tout soluté $\sigma$, on a ${\hat {\gamma }}_{s}^{\text{l}}=1$. Ceci est à fortiori vrai pour les mélanges binaires « soluté $\sigma$ - solvant $s$ ».
Pour tout solvant $s$, 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 solvant $s$ : $x_{s}^{\text{g}}\phi _{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 ${\bar {V}}^{\text{l}}$ de la phase liquide vaut :
    ${\begin{aligned}{{\bar {V}}^{\text{l}} \over RT}&={\sum _{\sigma }x_{\sigma }^{\text{l}}{\bar {V}}_{\sigma }^{\text{l}}+\sum _{s}x_{s}^{\text{l}}{\bar {V}}_{s}^{\text{l}} \over RT}=\sum _{\sigma }x_{\sigma }^{\text{l}}\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}+\sum _{s}x_{s}^{\text{l}}\left({\partial \ln f_{s}^{\text{l}} \over \partial P}\right)_{T,n}\\&=\sum _{\sigma }x_{\sigma }^{\text{l}}\left[\sum _{s}y_{s}^{\text{l}}\left({{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}-\left({\partial \ln \gamma _{\sigma ,s}^{{\text{l}},\infty } \over \partial P}\right)_{T,n}\right)+\left({\partial \ln \gamma _{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}\right]+\sum _{s}x_{s}^{\text{l}}\left[{{\bar {V}}_{s}^{{\text{l}},*} \over RT}+\left({\partial \ln \gamma _{s}^{\text{l}} \over \partial P}\right)_{T,n}\right]\end{aligned}}$
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 ${\bar {V}}_{\sigma }^{\text{l}}$ de tout soluté $\sigma$ vaut :
    ${\begin{aligned}{{\bar {V}}_{\sigma }^{\text{l}} \over RT}&=\left({\partial \ln f_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}=\underbrace {\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}} _{{\text{à composition constante}}\,=\,0}+\left({\partial \ln f_{\sigma }^{{\text{l}},*} \over \partial P}\right)_{T,n}+\left({\partial \ln \gamma _{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}\\&=\sum _{s}y_{s}^{\text{l}}\left({\partial \ln f_{\sigma ,s}^{{\text{l}},*} \over \partial P}\right)_{T,n}+\left({\partial \ln \gamma _{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}=\sum _{s}y_{s}^{\text{l}}\left(\left({\partial \ln k_{{\text{H}},\sigma ,s} \over \partial P}\right)_{T,n}-\left({\partial \gamma _{\sigma ,s}^{{\text{l}},\infty } \over \partial P}\right)_{T,n}\right)+\left({\partial \ln \gamma _{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}\\&=\sum _{s}y_{s}^{\text{l}}\left({{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}-\left({\partial \ln \gamma _{\sigma ,s}^{{\text{l}},\infty } \over \partial P}\right)_{T,n}\right)+\left({\partial \ln \gamma _{\sigma }^{\text{l}} \over \partial P}\right)_{T,n}\end{aligned}}$
Le volume molaire partiel ${\bar {V}}_{s}^{\text{l}}$ de tout soluté $s$ vaut :
    ${\begin{aligned}{{\bar {V}}_{s}^{\text{l}} \over RT}&=\left({\partial \ln f_{s}^{\text{l}} \over \partial P}\right)_{T,n}=\underbrace {\left({\partial \ln x_{s}^{\text{l}} \over \partial P}\right)_{T,n}} _{{\text{à composition constante}}\,=\,0}+\underbrace {\left({\partial \ln \!\left(\phi _{s}^{\text{g,*,sat}}P_{s}^{\text{sat}}\right) \over \partial P}\right)_{T,n}} _{{\text{à température constante}}\,=\,0}+\left({\partial \ln {\mathcal {P}}_{s}^{{\text{l}},*} \over \partial P}\right)_{T,n}+\left({\partial \ln \gamma _{s}^{\text{l}} \over \partial P}\right)_{T,n}+\left({\partial \ln {\hat {\gamma }}_{s}^{\text{l}} \over \partial P}\right)_{T,n}\\&={{\bar {V}}_{s}^{{\text{l}},*} \over RT}+\left({\partial \ln \gamma _{s}^{\text{l}} \over \partial P}\right)_{T,n}+\left({\partial \ln {\hat {\gamma }}_{s}^{\text{l}} \over \partial P}\right)_{T,n}\end{aligned}}$
avec :
    $\left({\partial \ln {\hat {\gamma }}_{s}^{\text{l}} \over \partial P}\right)_{T,n}=\sum _{\sigma }y_{\sigma }^{\text{l}}\left({{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}-\left({\partial \ln \gamma _{\sigma ,s}^{{\text{l}},\infty } \over \partial P}\right)_{T,n}-\sum _{r}y_{r}^{\text{l}}\left({{\bar {V}}_{\sigma ,r}^{{\text{l}},\infty } \over RT}-\left({\partial \ln \gamma _{\sigma ,r}^{{\text{l}},\infty } \over \partial P}\right)_{T,n}\right)\right)$
Note : dans le calcul du volume molaire ${\bar {V}}^{\text{l}}$, il apparait dans la somme $\sum _{s}x_{s}^{\text{l}}{\bar {V}}_{s}^{\text{l}}$ le terme :
    $\sum _{s}x_{s}^{\text{l}}\left({\partial \ln {\hat {\gamma }}_{s}^{\text{l}} \over \partial P}\right)_{T,n}=\left({\partial \sum _{s}x_{s}^{\text{l}}\ln {\hat {\gamma }}_{s} \over \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 corps $i$ ;
* $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 :
$i$
$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 corps $i$ on a :
    $y_{i}^{\text{l}}={n_{i}^{\text{l}} \over \sum _{s}n_{s}^{\text{l}}}={n_{i}^{\text{l}}/n^{\text{l}} \over \sum _{s}\left(n_{s}^{\text{l}}/n^{\text{l}}\right)}$
d'où :
$y_{i}^{\text{l}}={x_{i}^{\text{l}} \over \sum _{s}x_{s}^{\text{l}}}$
Pour tout solvant $s$, $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 corps $i$ on a :
    $n_{j}^{\text{l}}\left({\partial \ln x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{n_{\neq i}^{\text{l}}}={n_{j}^{\text{l}} \over x_{j}^{\text{l}}}\left({\partial x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{n_{\neq i}^{\text{l}}}=n^{\text{l}}\left({\partial x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{n_{\neq i}^{\text{l}}}$
    $\left({\partial x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{n_{\neq i}^{\text{l}}}={\begin{cases}{1 \over n^{\text{l}}}-{x_{i}^{\text{l}} \over n^{\text{l}}}&{\text{si}}\,j=i\\-{x_{j}^{\text{l}} \over n^{\text{l}}}&{\text{si}}\,j\neq i\end{cases}}$
    $\sum _{j}n_{j}^{\text{l}}\left({\partial \ln x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{n_{\neq i}^{\text{l}}}=1-\sum _{j}x_{j}$
d'où :
$\sum _{j}n_{j}^{\text{l}}\left({\partial \ln x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{n_{\neq 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 :
$G^{\text{l}}=\sum _{i}n_{i}^{\text{l}}\mu _{i}^{\text{l}}$
dans lequel, pour tout corps $i$, le potentiel chimique $\mu _{i}^{\text{l}}$ est défini comme l'enthalpie libre molaire partielle :
$\mu _{i}^{\text{l}}\equiv \left({\partial G^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}$
Par définition de la fugacité $f_{i}^{\text{l}}$, on a, en mélange :
$\mu _{i}^{\text{l}}=\mu _{i}^{\bullet ,*}+RT\,\ln {f_{i}^{\text{l}} \over P}$
et pour le corps pur :
    $\mu _{i}^{\text{l,*}}=\lim _{x_{i}^{\text{l}}\to 1 \atop \sum _{j\neq i}x_{j}^{\text{l}}\to 0}\mu _{i}^{\text{l}}=\mu _{i}^{\bullet ,*}+RT\,\ln {f_{i}^{\text{l,*}} \over P}$
avec :
* $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 :
$\sum _{i}n_{i}^{\text{l}}\,\mathrm {d} \mu _{i}^{\text{l}}=0$
d'où :
$i$
$\sum _{j}n_{j}^{\text{l}}\left({\partial \ln f_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq 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 :
${G^{\text{E}} \over RT}\equiv \sum _{i}n_{i}^{\text{l}}\,\ln \!\left({f_{i}^{\text{l}} \over x_{i}^{\text{l}}f_{i}^{\text{l,*}}}\right)=\sum _{i}n_{i}^{\text{l}}\,\ln \Gamma _{i}^{\text{l}}$
Le coefficient d'activité $\Gamma _{i}^{\text{l}}$ est défini par :
$\Gamma _{i}^{\text{l}}\equiv {f_{i}^{\text{l}} \over x_{i}^{\text{l}}f_{i}^{\text{l,*}}}$
Le coefficient $\gamma _{i}^{\text{l}}$ est défini par :
$\ln \gamma _{i}^{\text{l}}\equiv {\mu _{i}^{\text{E}} \over RT}\equiv \left({\partial \over \partial n_{i}^{\text{l}}}\left[{G^{\text{E}} \over RT}\right]\right)_{P,T,n_{\neq i}^{\text{l}}}$
avec $\mu _{i}^{\text{E}}$ le potentiel chimique d'excès. Le théorème d'Euler induit :
${G^{\text{E}} \over RT}=\sum _{i}n_{i}^{\text{l}}{\mu _{i}^{\text{E}} \over RT}=\sum _{i}n_{i}^{\text{l}}\,\ln \gamma _{i}^{\text{l}}$
Les coefficients $\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 :
$i$
$\sum _{j}n_{j}^{\text{l}}\left({\partial \ln \gamma _{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}=0$

###### Fugacité d'un soluté pur

Pour un mélange binaire « soluté $\sigma$ - solvant $s$ », 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,*}}=\lim _{x_{\sigma }^{\text{l}}\to 1 \atop x_{s}^{\text{l}}\to 0}\mu _{\sigma }^{\text{l}}=\mu _{\sigma }^{\bullet ,*}+RT\,\ln {f_{\sigma ,s}^{\text{l,*}} \over 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,*}}=\lim _{x_{\sigma }^{\text{l}}\to 1 \atop \sum _{s}x_{s}^{\text{l}}\to 0}\mu _{\sigma }^{\text{l}}=\mu _{\sigma }^{\bullet ,*}+RT\,\ln {f_{\sigma }^{\text{l,*}} \over P}$
Si dans le mélange de solvants le solvant $s$ 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,*}}=\underbrace {\left(\sum _{s}y_{s}^{\text{l}}\right)} _{=\,1}\mu _{\sigma }^{\bullet ,*}+RT\,\left(\sum _{s}y_{s}^{\text{l}}\,\ln f_{\sigma ,s}^{\text{l,*}}\right)-\underbrace {\left(\sum _{s}y_{s}^{\text{l}}\right)} _{=\,1}RT\,\ln P=\mu _{\sigma }^{\bullet ,*}+RT\,\left(\sum _{s}y_{s}^{\text{l}}\,\ln f_{\sigma ,s}^{\text{l,*}}\right)-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 :
$\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 corps $i$ on a :
    ${\begin{aligned}\left({\partial \over \partial n_{i}^{\text{l}}}\left[{G^{\text{E}} \over RT}\right]\right)_{P,T,n_{\neq i}^{\text{l}}}&=\left({\partial \over \partial n_{i}^{\text{l}}}\left[\sum _{j}n_{j}^{\text{l}}\,\ln \!\left({f_{j}^{\text{l}} \over x_{j}^{\text{l}}f_{j}^{\text{l,*}}}\right)\right]\right)_{P,T,n_{\neq i}^{\text{l}}}\\&=\ln \!\left({f_{i}^{\text{l}} \over x_{i}^{\text{l}}f_{i}^{\text{l,*}}}\right)+\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln f_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{{\text{Gibbs-Duhem}}\,=\,0}-\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{=\,0}-\sum _{j}n_{j}^{\text{l}}\left({\partial \ln f_{j}^{\text{l,*}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}\end{aligned}}$
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({\partial \ln f_{j}^{\text{l,*}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}$
On peut donc écrire :
    $\ln \gamma _{i}^{\text{l}}=\ln \!\left({f_{i}^{\text{l}} \over x_{i}^{\text{l}}f_{i}^{\text{l,*}}}\right)-\ln {\hat {\gamma }}_{i}^{\text{l}}$
On obtient :
$i$
$f_{i}^{\text{l}}=x_{i}^{\text{l}}\gamma _{i}^{\text{l}}{\hat {\gamma }}_{i}^{\text{l}}f_{i}^{\text{l,*}}$
soit :
$\Gamma _{i}^{\text{l}}={f_{i}^{\text{l}} \over x_{i}^{\text{l}}f_{i}^{\text{l,*}}}=\gamma _{i}^{\text{l}}{\hat {\gamma }}_{i}^{\text{l}}$
On a $\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 soluté $\sigma$, 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({\partial \ln f_{j}^{\text{l,*}} \over \partial n_{\sigma }^{\text{l}}}\right)_{P,T,n_{\neq \sigma }^{\text{l}}}=0$
On obtient :
$\sigma$
${\hat {\gamma }}_{\sigma }^{\text{l}}=1$
Pour un solvant $s$, 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({\partial \ln f_{j}^{\text{l,*}} \over \partial n_{s}^{\text{l}}}\right)_{P,T,n_{\neq s}^{\text{l}}}=\sum _{\sigma }n_{\sigma }^{\text{l}}\left({\partial \ln f_{\sigma }^{\text{l,*}} \over n_{s}^{\text{l}}}\right)_{P,T,n_{\neq 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}{n_{r}^{\text{l}} \over 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({\partial \ln f_{\sigma }^{\text{l,*}} \over \partial n_{s}^{\text{l}}}\right)_{P,T,n_{\neq s}^{\text{l}}}={1 \over n_{S}^{\text{l}}}\,\ln f_{\sigma ,s}^{\text{l,*}}-\sum _{r}{n_{r} \over {n_{S}^{\text{l}}}^{2}}\,\ln f_{\sigma ,r}^{\text{l,*}}={1 \over n_{S}^{\text{l}}}\left(\ln f_{\sigma ,s}^{\text{l,*}}-\sum _{r}y_{r}^{\text{l}}\,\ln f_{\sigma ,r}^{\text{l,*}}\right)={1 \over n_{S}^{\text{l}}}\left(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}}\right)$
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 :
$s$
$\ln {\hat {\gamma }}_{s}^{\text{l}}=\sum _{\sigma }y_{\sigma }^{\text{l}}\left(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}}\right)$
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}}\left(\sum _{r}y_{r}^{\text{l}}\ln f_{\sigma ,r}^{\text{l,*}}\right)$
On a donc également :
$\ln {\hat {\gamma }}_{s}^{\text{l}}=\sum _{\sigma }\sum _{r}y_{\sigma }^{\text{l}}y_{r}^{\text{l}}\left(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma ,r}^{\text{l,*}}\right)$

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

Pour tout corps $i$ on a $\Gamma _{i}^{\text{l}}=\gamma _{i}^{\text{l}}{\hat {\gamma }}_{i}^{\text{l}}$. Par définition l'enthalpie libre d'excès donne :
    ${G^{\text{E}} \over RT}=\sum _{i}n_{i}^{\text{l}}\ln \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 :
    ${G^{\text{E}} \over RT}=\sum _{i}n_{i}^{\text{l}}{\mu _{i}^{\text{E}} \over 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}}\left(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}}\right)$
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}}\left(\ln f_{\sigma ,s}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}}\right)=\sum _{\sigma }n_{\sigma }^{\text{l}}\underbrace {\left(\sum _{s}y_{s}^{\text{l}}\ln f_{\sigma ,s}^{{\text{l}},*}\right)} _{=\,\ln f_{\sigma }^{\text{l,*}}}-\underbrace {\left(\sum _{s}y_{s}^{\text{l}}\right)} _{=\,1}\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 corps $i$ 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({\partial \ln f_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}=\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{=\,0}+\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln \gamma _{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{{\text{Gibbs-Duhem}}\,=\,0}+\sum _{j}n_{j}^{\text{l}}\left({\partial \ln {\hat {\gamma }}_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}+\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln f_{j}^{\text{l,*}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{{\text{définition}}\,=\,\ln {\hat {\gamma }}_{i}^{\text{l}}}$
On a :
    $\sum _{j}n_{j}^{\text{l}}\left({\partial \ln {\hat {\gamma }}_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}=\underbrace {\left({\partial \over \partial n_{i}^{\text{l}}}\left[\sum _{j}n_{j}^{\text{l}}\,\ln {\hat {\gamma }}_{j}^{\text{l}}\right]\right)_{P,T,n_{\neq i}^{\text{l}}}} _{=\,0}-\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 corps $i$ :
    $\sum _{j}n_{j}^{\text{l}}\left({\partial \ln f_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq 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 soluté $\sigma$ à l'état de liquide pur est calculée selon la relation :
$f_{\sigma }^{\text{l,*}}={k_{{\text{H}},\sigma ,m} \over \gamma _{\sigma ,m}^{{\text{l}},\infty }}$
avec :
* $k_{{\text{H}},\sigma ,m}$ la constante de Henry du soluté $\sigma$ dans le mélange de solvants seuls ;
* $\gamma _{\sigma ,m}^{{\text{l}},\infty }$ le coefficient d'activité du soluté $\sigma$ à 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 soluté $\sigma$ la relation idéale de Krichevsky :
$\sigma$
$\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 :
$\sigma$
$x_{\sigma }^{\text{g}}P=x_{\sigma }^{\text{l}}k_{{\text{H}},\sigma ,m}$
On a par ailleurs pour tout solvant $s$ :
$s$
$\ln {\hat {\gamma }}_{s}^{\text{l}}=\sum _{\sigma }y_{\sigma }^{\text{l}}\left(\ln k_{{\text{H}},\sigma ,s}-\ln k_{{\text{H}},\sigma ,m}\right)$
Pour les basses pressions et les phases liquides idéales on a par conséquent la loi idéale pour les solvants :
$s$
$x_{s}^{\text{g}}P=x_{s}^{\text{l}}{\hat {\gamma }}_{s}^{\text{l}}P_{s}^{\text{sat}}{\mathcal {P}}_{s}^{{\text{l}},*}$
Le volume molaire ${\bar {V}}^{\text{l,id}}$ de la phase liquide idéale vaut :
    ${\bar {V}}^{\text{l,id}}=\sum _{\sigma }\sum _{s}x_{\sigma }^{\text{l}}y_{s}^{\text{l}}{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }+\sum _{s}x_{s}^{\text{l}}{\bar {V}}_{s}^{{\text{l}},*}$
L'enthalpie libre d'excès $G^{\text{E}}$ vaut :
    ${G^{\text{E}} \over RT}=\sum _{i}n_{i}^{\text{l}}\,\ln \!\left({f_{i}^{\text{l}} \over x_{i}^{\text{l}}f_{i}^{\text{l,*}}}\right)=\sum _{i}n_{i}^{\text{l}}\,\ln \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}}\left(\mu _{i}^{\bullet ,*}+RT\,\ln {f_{i}^{\text{l}} \over P}\right)$
en conséquence :
    $G^{\text{l,id}}=G^{\text{l}}-G^{\text{E}}=\sum _{i}n_{i}^{\text{l}}\left(\mu _{i}^{\bullet ,*}+RT\,\ln \!\left({x_{i}^{\text{l}}f_{i}^{\text{l,*}} \over P}\right)\right)$
On calcule le potentiel chimique idéal $\mu _{i}^{\text{l,id}}$ :
    ${\mu _{i}^{\text{l,id}} \over RT}=\left({\partial \over \partial n_{i}^{\text{l}}}\left[{G^{\text{l,id}} \over RT}\right]\right)_{P,T,n_{\neq i}^{\text{l}}}={\mu _{i}^{\bullet ,*} \over RT}+\ln \!\left({x_{i}^{\text{l}}f_{i}^{\text{l,*}} \over P}\right)+\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln x_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{=\,0}+\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln f_{j}^{\text{l,*}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{{\text{définition}}\,=\,\ln {\hat {\gamma }}_{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 {f_{i}^{\text{l,id}} \over P}$
    $f_{i}^{\text{l,id}}=x_{i}^{\text{l}}{\hat {\gamma }}_{i}^{\text{l}}f_{i}^{\text{l,*}}$
    $\Gamma _{i}^{\text{l,id}}={f_{i}^{\text{l,id}} \over 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}}\left(\mu _{i}^{\bullet ,*}+RT\,\ln \!\left({x_{i}^{\text{l}}f_{i}^{\text{l,*}} \over P}\right)\right)+RT\,\underbrace {\sum _{i}n_{i}^{\text{l}}\,\ln {\hat {\gamma }}_{i}^{\text{l}}} _{=\,0}=G^{\text{l,id}}$
et pour tout corps $i$, 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({\partial \ln f_{j}^{\text{l,id}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}=\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln f_{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{{\text{Gibbs-Duhem}}\,=\,0}-\underbrace {\sum _{j}n_{j}^{\text{l}}\left({\partial \ln \gamma _{j}^{\text{l}} \over \partial n_{i}^{\text{l}}}\right)_{P,T,n_{\neq i}^{\text{l}}}} _{{\text{Gibbs-Duhem}}\,=\,0}=0$

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

Soit un mélange liquide ternaire constitué d'un soluté $\sigma$ 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 ${\bar {G}}^{\text{E}}$ du mélange ternaire est calculable selon l'extension du modèle de Margules :
    ${{\bar {G}}^{\text{E}} \over 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 :
$\ln \gamma _{\sigma }^{\text{l}}=A_{\sigma ,s1}\,x_{s1}^{\text{l}}+A_{\sigma ,s2}\,x_{s2}^{\text{l}}-{\bar {G}}^{\text{E}}/RT$
$\ln \gamma _{s1}^{\text{l}}=A_{\sigma ,s1}\,x_{\sigma }^{\text{l}}+A_{s1,s2}\,x_{s2}^{\text{l}}-{\bar {G}}^{\text{E}}/RT$
$\ln \gamma _{s2}^{\text{l}}=A_{\sigma ,s2}\,x_{\sigma }^{\text{l}}+A_{s1,s2}\,x_{s1}^{\text{l}}-{\bar {G}}^{\text{E}}/RT$
Les coefficients d'activité du soluté $\sigma$ à dilution infinie dans chacun des deux solvants $s1$ et $s2$ sont donnés par :
    $\ln \gamma _{\sigma ,s1}^{{\text{l}},\infty }=\lim _{x_{\sigma }^{\text{l}}+x_{s2}^{\text{l}}\to 0 \atop x_{s1}^{\text{l}}\to 1}\ln \gamma _{\sigma }^{\text{l}}=A_{\sigma ,s1}$
    $\ln \gamma _{\sigma ,s2}^{{\text{l}},\infty }=\lim _{x_{\sigma }^{\text{l}}+x_{s1}^{\text{l}}\to 0 \atop x_{s2}^{\text{l}}\to 1}\ln \gamma _{\sigma }^{\text{l}}=A_{\sigma ,s2}$
On pose :
$y_{\sigma }^{\text{l}}={x_{\sigma }^{\text{l}} \over x_{s1}^{\text{l}}+x_{s2}^{\text{l}}}\quad ;\quad y_{s1}^{\text{l}}={x_{s1}^{\text{l}} \over x_{s1}^{\text{l}}+x_{s2}^{\text{l}}}\quad ;\quad y_{s2}^{\text{l}}={x_{s2}^{\text{l}} \over x_{s1}^{\text{l}}+x_{s2}^{\text{l}}}$
La fugacité $f_{\sigma }^{\text{l,*}}$ du soluté $\sigma$ à 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 {k_{{\text{H}},\sigma ,s1} \over \gamma _{\sigma ,s1}^{{\text{l}},\infty }}+y_{s2}^{\text{l}}\,\ln {k_{{\text{H}},\sigma ,s2} \over \gamma _{\sigma ,s2}^{{\text{l}},\infty }}$
d'où :
$\sigma$
$\ln f_{\sigma }^{\text{l,*}}=y_{s1}^{\text{l}}\left(\ln k_{{\text{H}},\sigma ,s1}-A_{\sigma ,s1}\right)+y_{s2}^{\text{l}}\left(\ln k_{{\text{H}},\sigma ,s2}-A_{\sigma ,s2}\right)$
Pour les deux solvants $s1$ et $s2$ on a :
    $\ln {\hat {\gamma }}_{s1}^{\text{l}}=y_{\sigma }^{\text{l}}\left(\ln f_{\sigma ,s1}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}}\right)=y_{\sigma }^{\text{l}}y_{s2}^{\text{l}}\left(\ln \!\left({k_{{\text{H}},\sigma ,s1} \over \gamma _{\sigma ,s1}^{{\text{l}},\infty }}\right)-\ln \!\left({k_{{\text{H}},\sigma ,s2} \over \gamma _{\sigma ,s2}^{{\text{l}},\infty }}\right)\right)$
    $\ln {\hat {\gamma }}_{s2}^{\text{l}}=y_{\sigma }^{\text{l}}\left(\ln f_{\sigma ,s2}^{\text{l,*}}-\ln f_{\sigma }^{\text{l,*}}\right)=y_{\sigma }^{\text{l}}y_{s1}^{\text{l}}\left(\ln \!\left({k_{{\text{H}},\sigma ,s2} \over \gamma _{\sigma ,s2}^{{\text{l}},\infty }}\right)-\ln \!\left({k_{{\text{H}},\sigma ,s1} \over \gamma _{\sigma ,s1}^{{\text{l}},\infty }}\right)\right)$
d'où :
$s1$
$s2$
$\ln {\hat {\gamma }}_{s1}^{\text{l}}=y_{\sigma }^{\text{l}}y_{s2}^{\text{l}}\,\left(\ln k_{{\text{H}},\sigma ,s1}-\ln k_{{\text{H}},\sigma ,s2}+A_{\sigma ,s2}-A_{\sigma ,s1}\right)$
$\ln {\hat {\gamma }}_{s2}^{\text{l}}=y_{\sigma }^{\text{l}}y_{s1}^{\text{l}}\,\left(\ln k_{{\text{H}},\sigma ,s2}-\ln k_{{\text{H}},\sigma ,s1}+A_{\sigma ,s1}-A_{\sigma ,s2}\right)$
Pour une solution idéale, soit ${\bar {G}}^{\text{E}}=0$, on a $\gamma _{\sigma }^{\text{l}}=\gamma _{s1}^{\text{l}}=\gamma _{s2}^{\text{l}}=\gamma _{\sigma ,s1}^{{\text{l}},\infty }=\gamma _{\sigma ,s2}^{{\text{l}},\infty }=\gamma _{\sigma ,m}^{{\text{l}},\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}}\,\left(\ln k_{{\text{H}},\sigma ,s1}-\ln k_{{\text{H}},\sigma ,s2}\right)$
    $\ln {\hat {\gamma }}_{s2}^{\text{l}}=y_{\sigma }^{\text{l}}y_{s1}^{\text{l}}\,\left(\ln k_{{\text{H}},\sigma ,s2}-\ln k_{{\text{H}},\sigma ,s1}\right)$
L'ensemble de ces expressions peut être étendu à des mélanges contenant plus de solutés et de solvants, avec :
    ${{\bar {G}}^{\text{E}} \over 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 :
    ${{\bar {G}}^{\text{E}} \over 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 solvant $s$ est très peu volatil, voire que le soluté $\sigma$ 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({\partial \ln P_{\sigma } \over \partial P}\right)_{T}=\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T}+\left({\partial \ln k_{{\text{H}},\sigma ,s} \over \partial P}\right)_{T}=\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T}+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}$
avec ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$ le volume molaire partiel du soluté $\sigma$ à dilution infinie dans le solvant $s$ (voir le paragraphe Constante de Henry - Dépendance à la pression). Le solvant $s$ étant négligeable en phase gaz, on a $P_{\sigma }=P$ et donc :
    $\left({\partial \ln P_{\sigma } \over \partial P}\right)_{T}={1 \over P}$
On obtient :
    $\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T}={1 \over P}-{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}$
En application de la loi des gaz parfaits, aux basses pressions le volume molaire de la phase gaz vaut : ${\bar {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 ${\bar {V}}^{\bullet }\gg {\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$, d'où :
$\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T}\approx {{\bar {V}}^{\bullet } \over RT}={1 \over P}$
Avec un point de solubilité connu $\left(x_{\sigma }^{{\text{l}},\circ },P^{\circ },T\right)$ on intègre à température constante :
${x_{\sigma }^{\text{l}} \over x_{\sigma }^{{\text{l}},\circ }}={P \over P^{\circ }}$
Aux basses pressions, la solubilité $x_{\sigma }^{\text{l}}$ augmente linéairement avec la pression totale $P$ lorsque le soluté $\sigma$ est seul en phase gaz.
Le volume molaire d'un gaz diminue avec une augmentation de pression ; par contre on peut considérer que ${\bar {V}}_{\sigma ,s}^{{\text{l}},\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 solvant $s$ étant négligé en phase gaz, la fugacité du soluté $\sigma$ dans cette phase varie en fonction de la pression selon :
    $\left({\partial \ln f_{\sigma }^{\text{g}} \over \partial P}\right)_{T}={{\bar {V}}_{\sigma }^{{\text{g}},*} \over RT}=\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T}+{{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}$
avec ${\bar {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équent$^,$ :
$\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial P}\right)_{T}={{\bar {V}}_{\sigma }^{{\text{g}},*}-{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty } \over RT}$
Lorsque ${\bar {V}}_{\sigma }^{{\text{g}},*}>{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$, aux basses pressions, la solubilité augmente avec la pression ; elle décroît lorsque ${\bar {V}}_{\sigma }^{{\text{g}},*}<{\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$, aux très hautes pressions (lorsque le volume molaire de la phase gaz est faible). Ainsi, lorsque ${\bar {V}}_{\sigma }^{{\text{g}},*}={\bar {V}}_{\sigma ,s}^{{\text{l}},\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 à 18 °C un maximum à environ 3 000 atm$^,$.

#### En fonction de la température

En dérivant, à pression constante, l'expression de la loi de Henry par rapport à la température, on obtient :
    $\left({\partial \ln P_{\sigma } \over \partial {1 \over T}}\right)_{P}=\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial {1 \over T}}\right)_{P}+\left({\partial \ln k_{{\text{H}},\sigma ,s} \over \partial {1 \over T}}\right)_{P}=\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial {1 \over T}}\right)_{P}+{\Delta _{\text{sol}}H_{\sigma ,s} \over R}$
avec $\Delta _{\text{sol}}H_{\sigma ,s}$ l'enthalpie de dissolution du soluté $\sigma$ dans le solvant $s$ (voir le paragraphe Constante de Henry - Dépendance à la température). Le solvant $s$ étant négligeable en phase gaz, on a $P_{\sigma }=P$ et donc :
    $\left({\partial \ln P_{\sigma } \over \partial {1 \over T}}\right)_{P}=0$
On obtient :
$\left({\partial \ln x_{\sigma }^{\text{l}} \over \partial {1 \over T}}\right)_{P}=-{\Delta _{\text{sol}}H_{\sigma ,s} \over R}$
Diverses formes de la solubilité en fonction de la température sont communément employées$^,$ :
    $\ln x_{\sigma }^{\text{l}}=a+{b \over T}+c\,\ln T+d\,T$
    $\ln x_{\sigma }^{\text{l}}=a+{b \over T}+{c \over T^{2}}+{d \over T^{3}}$
avec $a$, $b$, $c$ et $d$ des constantes empiriques spécifiques du couple « soluté $\sigma$ - solvant $s$ ». Les expressions des enthalpies de dissolution respectives sont :
    $-{\Delta _{\text{sol}}H_{\sigma ,s} \over R}=b-c\,T-d\,T^{2}$
    $-{\Delta _{\text{sol}}H_{\sigma ,s} \over R}=b+{2\,c \over T}+{3\,d \over T^{2}}$
Pour la plupart des gaz, la dissolution est exothermique aux basses pressions et températures, soit $\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 augmente$^,$. Ainsi, aux basses pressions, le minimum de solubilité de l'hélium dans l'eau se situe à environ 30 °C, ceux de l'argon, de l'oxygène et de l'azote se situent entre 92 et 93 °C et celui du xénon à environ 114 °C.
Solubilité de l'oxygène O$_2$ dans l'eau pour une pression partielle d'oxygène $P_{\rm {O_{2}}}$ de 1 atm. :
| Température $T$ (K) | Solubilité 10$^5$ × $x_{\rm {O_{2}}}^{\text{l}}$ | Enthalpie de dissolution $\Delta _{\text{sol}}H_{{\rm {O_{2}}},{\text{eau}}}$ (kJ mol$^{−1}$) | Constante de Henry 10$^{−9}$ × $k_{{\text{H}},{\rm {O_{2}}},{\text{eau}}}^{\text{sat}}$ (Pa) (à $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

Les boissons gazeuses sont maintenues sous une pression importante de dioxyde de carbone (CO$_2$). On suppose qu'il n'y a que du CO$_2$ 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 CO$_2$ 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 CO$_2$ dissout se désorbe en formant des bulles dans le liquide$^,$.

### 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 O$_2$ et d'autres gaz dont l'azote N$_2$ très majoritaire, soit les fractions molaires $x_{{\text{O}}_{2}}^{\text{g}}$ = 0,209 et $x_{{\text{N}}_{2}}^{\text{g}}$ ≈ 0,791. 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 à 10 m 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 plongeur$^,$.

### 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 à 25 °C (298,15 K).
Valeurs de la constante de Henry pour des gaz dissouts dans l'eau à 298,15 K. :
| Équation | $k_{\mathrm {H} }={P \over c_{\text{aq}}}$ | $H={c_{\text{aq}} \over P}$ | $k_{\mathrm {H} }={P \over x}$ | $H={c_{\text{aq}} \over c_{\text{gas}}}$ |
| Unité | atm l mol$^{−1}$ | mol l$^{−1}$ atm$^{−1}$ | atm | sans dimension |
| Oxygène (O$_2$) | 769,23 | 1,3 × 10$^{−3}$ | 4,259 × 10$^{−4}$ | 3,181 × 10$^{−2}$ |
| Hydrogène (H$_2$) | 1 282,05 | 7,8 × 10$^{−4}$ | 7,099 × 10$^{−4}$ | 1,907 × 10$^{−2}$ |
| Dioxyde de carbone (CO$_2$) | 29,41 | 3,4 × 10$^{−2}$ | 0,163 × 10$^{−4}$ | 0,831 7 |
| Azote (N$_2$) | 1 639,34 | 6,1 × 10$^{−4}$ | 9,077 × 10$^{−4}$ | 1,492 × 10$^{−2}$ |
| Hélium (He) | 2 702,7 | 3,7 × 10$^{−4}$ | 14,97 × 10$^{−4}$ | 9,051 × 10$^{−3}$ |
| Néon (Ne) | 2 222,22 | 4,5 × 10$^{−4}$ | 12,30 × 10$^{−4}$ | 1,101 × 10$^{−2}$ |
| Argon (Ar) | 714,28 | 1,4 × 10$^{−3}$ | 3,955 × 10$^{−4}$ | 3,425 × 10$^{−2}$ |
| Monoxyde de carbone (CO) | 1 052,63 | 9,5 × 10$^{−4}$ | 5,828 × 10$^{−4}$ | 2,324 × 10$^{−2}$ |
Avec :
* $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 corps $i$ en phase vapeur ;
* $f_{i}^{\text{l}}$ la fugacité du corps $i$ en phase liquide ;
* $f_{i}^{{\text{l}},*}$ la fugacité du corps $i$ liquide pur ;
* ${\bar {H}}_{\sigma ,s}^{{\text{l}},\infty }$ l'enthalpie molaire partielle du soluté $\sigma$ à dilution infinie dans le solvant $s$ ;
* ${\bar {H}}_{\sigma }^{\bullet ,*}$ l'enthalpie molaire du soluté $\sigma$ à 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 corps $i$ en phase liquide ;
* $k_{{\text{H}},\sigma ,s}$ la constante de Henry du soluté $\sigma$ dans le solvant $s$ ;
* $k_{{\text{H}},\sigma ,s}^{\text{sat}}$ la constante de Henry du soluté $\sigma$ dans le solvant $s$ à $P_{s}^{\text{sat}}$ ;
* $k_{{\text{H}},\sigma ,m}$ la constante de Henry du soluté $\sigma$ dans un mélange de solvants liquide ;
* $P$ la pression totale ;
* $P_{i}$ la pression partielle du corps $i$ ;
* $P_{s}^{\text{sat}}$ la pression de vapeur saturante du solvant $s$ ;
* ${\mathcal {P}}_{\sigma ,s}^{{\text{l}},\infty }$ le facteur de Poynting du soluté $\sigma$ à dilution infinie dans le solvant $s$ liquide ;
* ${\mathcal {P}}_{s}^{{\text{l}},*}$ le facteur de Poynting du solvant $s$ liquide pur ;
* $R$ la constante universelle des gaz parfaits ;
* $T$ la température thermodynamique ;
* ${\bar {V}}_{\sigma }^{{\text{g}},*}$ le volume molaire du soluté $\sigma$ gazeux pur ;
* $V^{\text{l}}$ le volume de la solution liquide ;
* ${\bar {V}}^{\text{l}}$ le volume molaire de la solution liquide ;
* ${\bar {V}}_{i}^{\text{l}}$ le volume molaire partiel du corps $i$ dans le mélange liquide ;
* ${\bar {V}}_{\sigma ,s}^{{\text{l}},\infty }$ le volume molaire partiel du soluté $\sigma$ à dilution infinie dans le solvant $s$ liquide ;
* ${\bar {V}}_{s}^{{\text{l}},*}$ le volume molaire du solvant $s$ liquide pur ;
* $x_{i}^{\text{g}}$ la fraction molaire du corps $i$ en phase gaz ;
* $x_{i}^{\text{l}}$ la fraction molaire du corps $i$ en phase liquide ;
* $y_{i}^{\text{l}}$ la quantité du corps $i$ rapportée à la quantité totale de solvants dans un mélange liquide multicomposant.

### Alphabet grec

* $\gamma _{i}^{\text{l}}$ le coefficient d'activité du corps $i$ en phase liquide ;
* $\gamma _{\sigma ,m}^{{\text{l}},\infty }$ le coefficient d'activité du soluté $\sigma$ à dilution infinie dans un mélange de solvants liquide ;
* $\gamma _{\sigma ,s}^{{\text{l}},\infty }$ le coefficient d'activité du soluté $\sigma$ à dilution infinie dans le solvant $s$ liquide ;
* ${\hat {\gamma }}_{s}^{\text{l}}$ le coefficient d'ajustement de la loi de Raoult pour le solvant $s$ dans un mélange multicomposant liquide ;
* $\Delta _{\text{sol}}H_{\sigma ,s}$ l'enthalpie de dissolution du soluté $\sigma$ dans le solvant $s$ ;
* $\phi _{i}^{\text{g}}$ le coefficient de fugacité du corps $i$ en phase gaz ;
* $\phi _{\sigma }^{\text{l}}$ le coefficient de fugacité du soluté $\sigma$ en phase liquide ;
* $\phi _{s}^{\text{g,*,sat}}$ le coefficient de fugacité du solvant $s$ pur en phase gaz à $P_{s}^{\text{sat}}$.