Potentiel chimique

En thermodynamique chimique, le potentiel chimique d'une espèce chimique est la variation d'énergie d'un système thermodynamique due à la variation de la quantité (nombre de moles) de cette espèce dans ce système. Étroitement lié au deuxième principe de la thermodynamique, le potentiel chimique permet d'étudier la stabilité des espèces chimiques et leur tendance à changer d'état, à réagir chimiquement ou à migrer par diffusion.

La fugacité et l'activité chimique, définies à partir du potentiel chimique, sont plus faciles à manipuler que celui-ci.

Historique

Le potentiel chimique est une notion introduite entre 1875 et 1878 par Willard Gibbs et Pierre Duhem. Il se révéla toutefois difficile à manipuler, car un potentiel chimique ne peut être calculé qu'à une constante additive près et non de façon absolue ; de plus le potentiel chimique de toute espèce tend vers l'infini négatif à dilution infinie. En 1900 et 1901 Gilbert Lewis introduit la fugacité, efficace dans son application pour les gaz, puis en 1923 l'activité chimique, plus spécialement employée pour les phases condensées (liquide ou solide).

Définition

Soit un mélange de $N$ constituants. Le potentiel chimique ${\mu }_i$ de tout constituant $i$ est défini pour chacun des quatre potentiels thermodynamiques :

• énergie interne (c'est par cette relation que Gibbs introduisit le potentiel chimique) :

$\mathrm{d}U=-P\mathrm{d}V+T\mathrm{d}S+\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i$

• énergie libre : $F=U-TS$

$\mathrm{d}F=-P\mathrm{d}V-S\mathrm{d}T+\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i$

• enthalpie : $H=U+PV$

$\mathrm{d}H=V\mathrm{d}P+T\mathrm{d}S+\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i$

• enthalpie libre : $G=H-TS$

$\mathrm{d}G=V\mathrm{d}P-S\mathrm{d}T+\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i$

avec :

• $P$ la pression ;
• $T$ la température ;
• $V$ le volume ;
• $S$ l'entropie ;
• $n_i$ la quantité (nombre de moles) du corps $i$.

Le potentiel chimique, en tant que dérivée partielle d'un potentiel thermodynamique, peut donc être défini de diverses façons, toutes équivalentes :

Potentiel chimique : ${\mu }_i={\left(\frac{\partial U}{\partial n_i}\right)}_{V,S,n_{j\ne i}}={\left(\frac{\partial F}{\partial n_i}\right)}_{V,T,n_{j\ne i}}={\left(\frac{\partial H}{\partial n_i}\right)}_{P,S,n_{j\ne i}}={\left(\frac{\partial G}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$

Le potentiel chimique a la dimension d'une énergie molaire et s'exprime en joules par mole, Modèle:Unité, dans le système international d'unités (SI).

La relation liant le potentiel chimique à l'énergie libre $F$ est particulièrement importante en calorimétrie, où les réactions chimiques sont étudiées à volume et température constants. Elle tient également une place centrale en physique statistique.

C'est toutefois la dernière de ces définitions, celle qui lie le potentiel chimique à l'enthalpie libre $G$, qui est la plus importante, car les réactions chimiques sont généralement étudiées à pression et température constantes (voir l'article Équilibre chimique).

Comme les potentiels thermodynamiques, le potentiel chimique ne peut pas être calculé de façon absolue, il est toujours établi à une constante additive près. Cette indétermination et sa limite dans le cas de la dilution infinie rendent son usage difficile. D'autres grandeurs définies à partir du potentiel chimique permettent de pallier cette difficulté : la fugacité et l'activité chimique peuvent être calculées de façon absolue et ne présentent pas de limite problématique.

Relation avec l'enthalpie libre

Enthalpie libre molaire partielle

Le potentiel chimique est lié de façon particulière à l'enthalpie libre $G$, car il s'agit du seul potentiel thermodynamique dont le potentiel chimique est la grandeur molaire partielle :

Enthalpie libre molaire partielle : ${\mu }_i={\overline{G}}_i={\left(\frac{\partial G}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$

Identité d'Euler

Selon le théorème d'Euler sur les fonctions homogènes du premier ordre on peut écrire, pour toute grandeur extensive $X$, la relation entre cette grandeur et les grandeurs molaires partielles ${\overline{X}}_i$ :

$X=\sum _{i=1}^N{n}_i{\overline{X}}_i$

Dans le cas particulier de l'enthalpie libre cela implique, étant donné l'identité des potentiels chimiques et des enthalpies libres molaires partielles :

$G=\sum _{i=1}^N{n}_i{\overline{G}}_i$

${\mu }_i={\overline{G}}_i$

Identité d'Euler pour l'enthalpie libre : $G=\sum _{i=1}^N{n}_i{\mu }_i$

En divisant par $n=\sum _{i=1}^N{n}_i$ la quantité totale de matière dans le mélange, on a également la relation :

Enthalpie libre molaire du mélange : $\overline{G}=\sum _i^N{x}_i{\mu }_i$

avec :

• $\overline{G}=\frac{G}{n}$, l'enthalpie libre molaire du mélange ;
• $x_i=\frac{n_i}{n}$, la fraction molaire du corps $i$ dans le mélange.

Relation de Gibbs-Duhem

En différentiant l'identité d'Euler pour l'enthalpie libre, on obtient :

Modèle:Equarefa $\mathrm{d}G=\sum _{i=1}^N{n}_i\mathrm{d}{\mu }_i+\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i$

La différentielle de l'enthalpie libre s'écrit, dans ses variables naturelles :

Modèle:Equarefa $\mathrm{d}G=V\mathrm{d}P-S\mathrm{d}T+\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i$

On peut identifier les différents termes des équations Modèle:Equarefl et Modèle:Equarefl, on obtient la relation de Gibbs-Duhem :

Relation de Gibbs-Duhem : $V\mathrm{d}P-S\mathrm{d}T=\sum _{i=1}^N{n}_i\mathrm{d}{\mu }_i$

Cette relation n'est valable que pour des transformations réversibles dans lesquelles le travail n'est dû qu'à des forces de pression.

Autres relations avec l'enthalpie libre molaire

Puisque ${\mu }_i={\left(\frac{\partial G}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$ et $G=n\overline{G}$ on peut écrire :

${\mu }_i=\overline{G}+n{\left(\frac{\partial \overline{G}}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$

et, puisque $\overline{G}$ est une fonction des fractions molaires $x$, en employant le théorème de dérivation des fonctions composées :

${\mu }_i=\overline{G}+{\left(\frac{\partial \overline{G}}{\partial x_i}\right)}_{P,T,x_{k\ne i}}-\sum _{j=1}^N{x}_j{\left(\frac{\partial \overline{G}}{\partial x_j}\right)}_{P,T,x_{k\ne j}}$

avec :

• $\overline{G}$ l'enthalpie libre molaire ;
• $n=\sum _{i=1}^N{n}_j$ la quantité de matière totale dans le mélange ;
• $n_i$ la quantité du corps $i$ ;
• $x_i$ la fraction molaire du corps $i$, $x_i=n_i/n$.

Autres relations

Relations avec d'autres grandeurs molaires partielles

En appliquant le théorème de Schwarz, le potentiel chimique est lié à d'autres grandeurs molaires partielles :

• ${\overline{V}}_i={\left(\frac{\partial V}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$ le volume molaire partiel, d'après l'équation d'état $V={\left(\frac{\partial G}{\partial P}\right)}_{T,n}$ :

${\overline{V}}_i={\left(\frac{\partial {\mu }_i}{\partial P}\right)}_{T,n}$

• ${\overline{S}}_i={\left(\frac{\partial S}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$ l'entropie molaire partielle, d'après l'équation d'état $S=-{\left(\frac{\partial G}{\partial T}\right)}_{P,n}$ :

${\overline{S}}_i=-{\left(\frac{\partial {\mu }_i}{\partial T}\right)}_{P,n}$

• ${\overline{H}}_i={\left(\frac{\partial H}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$ l'enthalpie molaire partielle, d'après la relation de Gibbs-Helmholtz $H={\left(\frac{\partial \frac{G}{T}}{\partial \frac{1}{T}}\right)}_{P,n}$ :

${\overline{H}}_i={\left(\frac{\partial \frac{{\mu }_i}{T}}{\partial \frac{1}{T}}\right)}_{P,n}$

• ${\overline{F}}_i={\left(\frac{\partial F}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$ l'énergie libre molaire partielle, d'après ${\overline{F}}_i={\overline{G}}_i-P{\overline{V}}_i$ :

${\overline{F}}_i={\mu }_i-P{\left(\frac{\partial {\mu }_i}{\partial P}\right)}_{T,n}$

• ${\overline{U}}_i={\left(\frac{\partial U}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$ l'énergie interne molaire partielle, d'après ${\overline{U}}_i={\overline{H}}_i-P{\overline{V}}_i$ :

${\overline{U}}_i={\left(\frac{\partial \frac{{\mu }_i}{T}}{\partial \frac{1}{T}}\right)}_{P,n}-P{\left(\frac{\partial {\mu }_i}{\partial P}\right)}_{T,n}$

Enfin, par définition de l'enthalpie libre, $G=H-TS=U+PV-TS=F+PV$, en dérivant selon $n_i$ à pression, température et quantités des autres composants constantes :

${\mu }_i={\overline{G}}_i={\overline{H}}_i-T{\overline{S}}_i={\overline{U}}_i+P{\overline{V}}_i-T{\overline{S}}_i={\overline{F}}_i+P{\overline{V}}_i$

Variation isotherme du potentiel chimique

Les équations d'état et le théorème de Schwarz permettent d'établir la variation du potentiel chimique en fonction de la pression $P$, la température $T$ ou l'entropie $S$. Les deux relations les plus importantes sont les suivantes :

• avec ${\mu }_i={\left(\frac{\partial F}{\partial n_i}\right)}_{V,T,n_{j\ne i}}$ et ${\left(\frac{\partial F}{\partial V}\right)}_{T,n}=-P$, on obtient :

${\left(\frac{\partial {\mu }_i}{\partial V}\right)}_{T,n}=-{\left(\frac{\partial P}{\partial n_i}\right)}_{V,T,n_{j\ne i}}$

• avec ${\mu }_i={\left(\frac{\partial G}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$ et ${\left(\frac{\partial G}{\partial P}\right)}_{T,n}=V$, on obtient :

${\left(\frac{\partial {\mu }_i}{\partial P}\right)}_{T,n}={\left(\frac{\partial V}{\partial n_i}\right)}_{P,T,n_{j\ne i}}$

cette dernière expression étant le volume molaire du corps $i$ pur ou le volume molaire partiel du corps $i$ dans le mélange, noté ${\overline{V}}_i$ :

${\left(\frac{\partial {\mu }_i}{\partial P}\right)}_{T,n}={\overline{V}}_i$

On a par conséquent, avec $f_i$ la fugacité du corps $i$ :

À température et composition constantes : $\mathrm{d}{\mu }_i=RT\mathrm{d}\ln f_i=-{\left(\frac{\partial P}{\partial n_i}\right)}_{V,T,n_{j\ne i}}\mathrm{d}V={\overline{V}}_i\mathrm{d}P$

Dépendance à la composition

Le potentiel chimique de tout corps $i$ à pression et température constantes, dans un mélange quelconque dans lequel il est représenté par la fraction molaire $x_i$, admet deux limites :

• pour le corps $i$ pur :

$\underset{x_i\to 1}{\lim }{\mu }_i={\mu }_i^{*}$

Cette valeur finie ne dépend que de la pression $P$, de la température $T$ et de la nature de $i$ ; par définition ${\mu }_i^{*}={\overline{G}}_i^{*}$ l'enthalpie libre molaire du corps $i$ pur à $P$ et $T$ ;

• pour le corps $i$ (soluté) infiniment dilué dans un solvant $s$Modèle:, :

$\underset{\genfrac{}{}{0ex}{}{x_i\to 0}{x_s\to 1}}{\lim }{\mu }_i=-\mathrm{\infty }$

À dilution infinie d'une espèce $i$ dans une espèce $s$, que la solution soit gazeuse, liquide ou solide, l'espèce $s$ exerce une attraction infiniment forte sur l'espèce $i$. Il est par conséquent impossible d'obtenir l'espèce $s$ infiniment pure. Une espèce purifiée a naturellement tendance à se polluer par attraction d'autres espèces présentes dans son environnement (par exemple, dans le matériau de son contenant).

Les potentiels chimiques peuvent être écrits en fonction de la pression, de la température et de la composition :

${\mu }_i={\mu }_i(P,T,n)$

$\mathrm{d}{\mu }_i={\left(\frac{\partial {\mu }_i}{\partial P}\right)}_{T,n}\mathrm{d}P+{\left(\frac{\partial {\mu }_i}{\partial T}\right)}_{P,n}\mathrm{d}T+\sum _{j=1}^N{\left(\frac{\partial {\mu }_i}{\partial n_j}\right)}_{P,T,n_{k\ne j}}\mathrm{d}{n}_j$

On a ainsi, par les relations avec le volume molaire partiel ${\overline{V}}_i$ et l'entropie molaire partielle ${\overline{S}}_i$ vues plus haut :

$\mathrm{d}{\mu }_i={\overline{V}}_i\mathrm{d}P-{\overline{S}}_i\mathrm{d}T+\sum _{j=1}^N{\left(\frac{\partial {\mu }_i}{\partial n_j}\right)}_{P,T,n_{k\ne j}}\mathrm{d}{n}_j$

Pour un corps pur $i$ à pression $P$ et température $T$ constantes, la relation précédente donne : $\mathrm{d}{\mu }_i^{*}={\left(\frac{\partial {\mu }_i^{*}}{\partial n_i}\right)}_{P,T}\mathrm{d}{n}_i$. Or dans ces conditions la relation de Gibbs-Duhem donne : $n_i\mathrm{d}{\mu }_i^{*}=0$. On en déduit que :

${\left(\frac{\partial {\mu }_i^{*}}{\partial n_i}\right)}_{P,T}=0$

À pression et température constantes, le potentiel chimique ${\mu }_i^{*}$ d'un corps pur ne dépend pas de la quantité $n_i$ de ce constituant.

Pour un mélange, les conditions de stabilité d'un équilibre thermodynamique imposent que :

${\left(\frac{\partial {\mu }_i}{\partial n_i}\right)}_{V,S,n_{k\ne i}}\ge 0$ ; ${\left(\frac{\partial {\mu }_i}{\partial n_i}\right)}_{V,T,n_{k\ne i}}\ge 0$ ; ${\left(\frac{\partial {\mu }_i}{\partial n_i}\right)}_{P,S,n_{k\ne i}}\ge 0$ ; ${\left(\frac{\partial {\mu }_i}{\partial n_i}\right)}_{P,T,n_{k\ne i}}\ge 0$

Le potentiel chimique ${\mu }_i$ d'une espèce chimique $i$ quelconque dans un mélange augmente lorsque l'on augmente la quantité $n_i$ de cette espèce, les quantités des autres constituants restant constantes (en l'absence de réaction chimique), et ce quelles que soient les autres conditions opératoires également maintenues constantes.

En appliquant le théorème de Schwarz aux relations définissant le potentiel chimique par rapport aux potentiels thermodynamiques, on obtient, pour toutes espèces chimiques $i$ et $j$ d'un mélange :

${\left(\frac{\partial {\mu }_i}{\partial n_j}\right)}_{V,S,n_{k\ne j}}={\left(\frac{\partial {\mu }_j}{\partial n_i}\right)}_{V,S,n_{k\ne i}}$ ${\left(\frac{\partial {\mu }_i}{\partial n_j}\right)}_{V,T,n_{k\ne j}}={\left(\frac{\partial {\mu }_j}{\partial n_i}\right)}_{V,T,n_{k\ne i}}$ ${\left(\frac{\partial {\mu }_i}{\partial n_j}\right)}_{P,S,n_{k\ne j}}={\left(\frac{\partial {\mu }_j}{\partial n_i}\right)}_{P,S,n_{k\ne i}}$ ${\left(\frac{\partial {\mu }_i}{\partial n_j}\right)}_{P,T,n_{k\ne j}}={\left(\frac{\partial {\mu }_j}{\partial n_i}\right)}_{P,T,n_{k\ne i}}$

À pression et températures constantes, la relation de Gibbs-Duhem donne : $\sum _{i=1}^N{n}_i\mathrm{d}{\mu }_i=0$. On en déduit que, pour toute espèce chimique $j$ du mélange :

$\sum _{i=1}^N{n}_i{\left(\frac{\partial {\mu }_i}{\partial n_j}\right)}_{P,T,n_{k\ne j}}=0$

En substituant la dernière relation obtenue plus haut par le théorème de Schwarz, on obtient :

$\sum _{i=1}^N{n}_i{\left(\frac{\partial {\mu }_j}{\partial n_i}\right)}_{P,T,n_{k\ne i}}=0$

Équilibres

Équilibre de phases

Soit un système isolé de volume $V$, contenant $N$ espèces chimiques : par définition il n'échange ni travail, ni chaleur, ni matière avec l'extérieur. On considère également que ce système n'est le siège d'aucune réaction chimique, aussi peut-on écrire :

• le système n'échange pas de travail avec l'extérieur, ses parois sont indéformables : $\mathrm{d}V=0$ ;
• le système n'échange pas de chaleur avec l'extérieur, il est isentropique : $\mathrm{d}S=0$ ;
• le système n'échange pas de matière avec l'extérieur, en l'absence de réaction chimique on peut écrire pour chaque corps $i$ : $\mathrm{d}{n}_i=0$.

On a donc, selon le premier principe de la thermodynamique, conservation de l'énergie interne $U$ pour l'ensemble du système :

$\mathrm{d}U=-P\mathrm{d}V+T\mathrm{d}S+\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i=0$

Si le système isolé contient deux phases, appelées $\alpha$ et $\beta$, à l'équilibre thermodynamique, on écrit pour chacune des deux phases la variation d'énergie interne :

$\mathrm{d}{U}^{\alpha }=-P^{\alpha }\mathrm{d}{V}^{\alpha }+T^{\alpha }\mathrm{d}{S}^{\alpha }+\sum _{i=1}^N{\mu }_i^{\alpha }\mathrm{d}{n}_i^{\alpha }$

$\mathrm{d}{U}^{\beta }=-P^{\beta }\mathrm{d}{V}^{\beta }+T^{\beta }\mathrm{d}{S}^{\beta }+\sum _{i=1}^N{\mu }_i^{\beta }\mathrm{d}{n}_i^{\beta }$

avec :

• $U^{\alpha }$ et $U^{\beta }$ les énergies internes respectives des deux phases ;
• $P^{\alpha }$ et $P^{\beta }$ les pressions respectives des deux phases ;
• $V^{\alpha }$ et $V^{\beta }$ les volumes respectifs des deux phases ;
• $T^{\alpha }$ et $T^{\beta }$ les températures respectives des deux phases ;
• $S^{\alpha }$ et $S^{\beta }$ les entropies respectives des deux phases ;
• ${\mu }_i^{\alpha }$ et ${\mu }_i^{\beta }$ les potentiels chimiques respectifs du corps $i$ dans chacune des deux phases ;
• $n_i^{\alpha }$ et $n_i^{\beta }$ les quantités respectives du corps $i$ dans chacune des deux phases.

Avec les relations sur le système global isolé, en l'absence de réaction chimique :

• énergie interne totale $U=U^{\alpha }+U^{\beta }$ constante : $\mathrm{d}U=\mathrm{d}{U}^{\alpha }+\mathrm{d}{U}^{\beta }=0$, selon le premier principe de la thermodynamique de conservation de l'énergie ;
• volume total $V=V^{\alpha }+V^{\beta }$ constant : $\mathrm{d}V=\mathrm{d}{V}^{\alpha }+\mathrm{d}{V}^{\beta }=0$ ;
• quantités de matière totales $n_i=n_i^{\alpha }+n_i^{\beta }$ constantes : $\mathrm{d}{n}_i=\mathrm{d}{n}_i^{\alpha }+\mathrm{d}{n}_i^{\beta }=0$ pour tout $i\in [1,\cdots ,N]$.

Attention : ce sont les grandeurs totales du système global isolé (de l'ensemble constitué par les deux phases) qui sont constantes, pas les grandeurs propres à chacune des deux phases qui, elles, peuvent varier au cours de la transformation, leurs variations étant opposées : par ex., les volumes des phases $V^{\alpha }$ et $V^{\beta }$ peuvent varier, mais le volume total des deux phases $V=V^{\alpha }+V^{\beta }$ est constant et $\mathrm{d}{V}^{\beta }=-\mathrm{d}{V}^{\alpha }$.

On obtient la relation globale pour l'ensemble des deux phases, selon le premier principe de la thermodynamique :

$\mathrm{d}U=-(P^{\alpha }-P^{\beta })\mathrm{d}{V}^{\alpha }+T^{\alpha }\mathrm{d}{S}^{\alpha }+T^{\beta }\mathrm{d}{S}^{\beta }+\sum _{i=1}^N({\mu }_i^{\alpha }-{\mu }_i^{\beta })\mathrm{d}{n}_i^{\alpha }=0$

Selon le deuxième principe de la thermodynamique, l'entropie totale $S=S^{\alpha }+S^{\beta }$ du système isolé ne peut qu'augmenter. L'équilibre est atteint lorsque l'entropie totale du système a atteint un maximum : les variations d'entropie des deux phases sont alors liées par : $\mathrm{d}S=\mathrm{d}{S}^{\alpha }+\mathrm{d}{S}^{\beta }=0$. On peut donc écrire pour le système global à l'équilibre :

À l'équilibre : $\mathrm{d}U=-(P^{\alpha }-P^{\beta })\mathrm{d}{V}^{\alpha }+(T^{\alpha }-T^{\beta })\mathrm{d}{S}^{\alpha }+\sum _{i=1}^N({\mu }_i^{\alpha }-{\mu }_i^{\beta })\mathrm{d}{n}_i^{\alpha }=0$

Quand les deux phases sont à l'équilibre, cela ne signifie pas qu'il n'y a pas d'échanges entre elles. Les deux phases continuent à échanger de la matière et de l'énergie, mais globalement les échanges se font dans les deux sens et se compensent : l'état d'équilibre obtenu dans ce cas peut être qualifié de dynamique ou stationnaire. Les variations $\mathrm{d}{V}^{\alpha }$, $\mathrm{d}{S}^{\alpha }$ et $\mathrm{d}{n}_i^{\alpha }$ pour tout $i\in [1,\cdots ,N]$ sont donc non nulles à l'équilibre. La relation précédente implique par conséquent l'homogénéité des potentiels mécanique (pression), thermique (température) et chimiques (potentiels chimiques) entre les deux phases à l'équilibre :

• $P^{\alpha }=P^{\beta }$ ;
• $T^{\alpha }=T^{\beta }$ ;
• ${\mu }_i^{\alpha }={\mu }_i^{\beta }$ pour tout $i\in [1,\cdots ,N]$.

Une autre façon d'obtenir ce résultat est de considérer que lorsque l'équilibre est atteint, c'est-à-dire lorsque l'entropie du système isolé est maximale, alors l'énergie totale du système (constante) correspond à un minimum par rapport aux diverses variables du système, ce qui implique que :

• ${\left(\frac{\partial U}{\partial V^{\alpha }}\right)}_{S^{\alpha },n^{\alpha }}=-(P^{\alpha }-P^{\beta })=0$ ;
• ${\left(\frac{\partial U}{\partial S^{\alpha }}\right)}_{V^{\alpha },n^{\alpha }}=T^{\alpha }-T^{\beta }=0$ ;
• ${\left(\frac{\partial U}{\partial n_i^{\alpha }}\right)}_{V^{\alpha },S^{\alpha },n_{j\ne i}^{\alpha }}={\mu }_i^{\alpha }-{\mu }_i^{\beta }=0$ pour tout $i\in [1,\cdots ,N]$.

Les conditions d'équilibre sont donc données par :

Conditions d'équilibre de deux phases $\alpha$ et $\beta$

Équilibre mécanique	homogénéité de la pression $P$
Équilibre thermique	homogénéité de la température $T$
Équilibre chimique	homogénéité des potentiels chimiques pour tout corps $i$ : ${\mu }_i^{\alpha }={\mu }_i^{\beta }$

Ces conditions sont généralisables à plus de deux phases en équilibre (gaz-liquide-liquide ou gaz-liquide-solide par exemple).

Note sur l'homogénéité des potentiels chimiques

L'homogénéité des potentiels chimiques signifie que pour le corps 1 ${\mu }_1^{\alpha }={\mu }_1^{\beta }$, pour le corps 2 ${\mu }_2^{\alpha }={\mu }_2^{\beta }$, … , pour le corps $N$ ${\mu }_N^{\alpha }={\mu }_N^{\beta }$, mais pas que ces potentiels sont égaux entre corps différents : ${\mu }_1^{\alpha }\ne {\mu }_2^{\alpha }\ne \cdots \ne {\mu }_N^{\alpha }$.

Équilibre chimique

Soit une réaction chimique quelconque notée sous la forme :

${\nu }_1{\mathrm{C}}_1+{\nu }_2{\mathrm{C}}_2+\cdots +{\nu }_N{\mathrm{C}}_N=0$

en attribuant selon la convention stœchiométrique une valeur négative aux coefficients stœchiométriques des réactifs, et positive à ceux des produits :

• ${\nu }_i<0$ pour un réactif ;
• ${\nu }_i>0$ pour un produit ;
• ${\nu }_i=0$ pour un inerte.

L'affinité chimique $𝒜$ est formellement définie par :

Affinité chimique : $𝒜=-\sum _{i=1}^N{\nu }_i{\mu }_i$

On note $\xi$ l'avancement de réaction et $n_i$ la quantité du composant $i$. On a $\mathrm{d}\xi =\frac{\mathrm{d}{n}_i}{{\nu }_i}$ pour tout $i\in [1,\cdots ,N]$. Le deuxième principe de la thermodynamique implique la condition d'évolution spontanée de toute réaction chimique :

Condition d'évolution spontanée : $𝒜\mathrm{d}\xi =-\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i>0$

$𝒜$ et $\mathrm{d}\xi$ ne peuvent donc être que de même signe :

• $𝒜>0$ et $\mathrm{d}\xi >0$, la réaction progresse : des réactifs sont consommés et des produits apparaissent,
• $𝒜<0$ et $\mathrm{d}\xi <0$, la réaction régresse : des produits sont consommés et des réactifs apparaissent.

En termes de potentiel thermodynamique, selon les conditions de la réaction, ceci va se traduire :

• par une décroissance de l'énergie libre $F$ pour une évolution à volume $V$ et température $T$ constants : $\mathrm{d}F=\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i<0$ ;
• par une décroissance de l'enthalpie libre $G$ pour une évolution à pression $P$ et température $T$ constantes : $\mathrm{d}G=\sum _{i=1}^N{\mu }_i\mathrm{d}{n}_i<0$.

Dans les deux cas, le potentiel thermodynamique décroît. Lorsqu'il atteint un minimum, ce qui correspond à un maximum d'entropie, le système atteint un équilibre.

À l'équilibre, l'entropie atteint un maximum et ne varie plus, alors $𝒜\mathrm{d}\xi =0$. Lorsqu'une réaction est équilibrée, cela ne signifie pas qu'il n'y a plus de transformation chimique dans le milieu : des réactifs continuent à se transformer en produits (réaction directe), et des produits à se transformer en réactifs (réaction inverse). L'état d'équilibre obtenu dans ce cas peut être qualifié de dynamique ou stationnaire : les réactions ont toujours lieu, mais globalement leurs effets s'annulent ; par exemple l'une absorbe totalement la chaleur dégagée par l'autre, d'où un bilan énergétique globalement nul. En conséquence, à l'équilibre $\mathrm{d}\xi \ne 0$, d'où :

À l'équilibre : $𝒜=0$

Migration

Lorsqu'un milieu est hétérogène, le potentiel chimique de chaque espèce n'est pas identique en tout point du milieu. Spontanément, chaque espèce migre vers les lieux où son potentiel chimique est le plus bas : l'entropie du système est ainsi maximisée, conformément au deuxième principe de la thermodynamique. Ce phénomène est appelé diffusion de la matière.

On considère un système thermodynamique isolé : il n'échange ni travail, ni chaleur, ni matière avec l'extérieur. Ceci induit notamment que le volume $V$ du système est constant. On suppose que ce système est composé de deux compartiments, $\alpha$ et $\beta$, séparés par une paroi perméable à la seule espèce chimique $i$. Ces deux compartiments sont supposés être en permanence à l'équilibre mécanique (même pression $P$) et à l'équilibre thermique (même température $T$). Enfin, ce système n'est le siège d'aucune réaction chimique.

Il n'est fait aucune hypothèse quant aux autres espèces chimiques du système (qui peuvent donc être différentes d'un compartiment à l'autre, ou à des concentrations différentes de part et d'autre de la paroi). Les deux compartiments peuvent être dans des phases différentes. La pression $P$ et la température $T$ peuvent varier. Les volumes des deux compartiments peuvent également varier du moment que le volume global $V$ du système reste constant (la paroi séparant les deux compartiments peut être mobile).

La variation d'énergie de chacun des compartiments vaut, en considérant que seule la quantité de l'espèce $i$ varie dans les compartiments :

$\mathrm{d}{U}^{\alpha }=-P\mathrm{d}{V}^{\alpha }+T\mathrm{d}{S}^{\alpha }+{\mu }_i^{\alpha }\mathrm{d}{n}_i^{\alpha }$

$\mathrm{d}{U}^{\beta }=-P\mathrm{d}{V}^{\beta }+T\mathrm{d}{S}^{\beta }+{\mu }_i^{\beta }\mathrm{d}{n}_i^{\beta }$

avec :

• $n_i^{\alpha }$ et $n_i^{\beta }$ les quantités respectives de l'espèce $i$ dans les deux compartiments ;
• $S^{\alpha }$ et $S^{\beta }$ les entropies respectives des deux compartiments ;
• $U^{\alpha }$ et $U^{\beta }$ les énergies internes respectives des deux compartiments ;
• $V^{\alpha }$ et $V^{\beta }$ les volumes respectifs des deux compartiments ;
• ${\mu }_i^{\alpha }$ et ${\mu }_i^{\beta }$ les potentiels chimiques respectifs de l'espèce $i$ dans les deux compartiments.

Le système étant isolé, on a, étant donné le premier principe de la thermodynamique, conservation de son énergie interne totale $U=U^{\alpha }+U^{\beta }$ :

$\mathrm{d}U=\mathrm{d}{U}^{\alpha }+\mathrm{d}{U}^{\beta }=-P(\mathrm{d}{V}^{\alpha }+\mathrm{d}{V}^{\beta })+T(\mathrm{d}{S}^{\alpha }+\mathrm{d}{S}^{\beta })+{\mu }_i^{\alpha }\mathrm{d}{n}_i^{\alpha }+{\mu }_i^{\beta }\mathrm{d}{n}_i^{\beta }=0$

Le volume total $V=V^{\alpha }+V^{\beta }$ étant constant, on a :

$\mathrm{d}V=\mathrm{d}{V}^{\alpha }+\mathrm{d}{V}^{\beta }=0$

On pose $S=S^{\alpha }+S^{\beta }$ l'entropie totale du système, on a :

$\mathrm{d}S=\mathrm{d}{S}^{\alpha }+\mathrm{d}{S}^{\beta }$

En l'absence d'échange de matière avec l'extérieur et de réaction chimique, on a conservation de la quantité totale $n_i=n_i^{\alpha }+n_i^{\beta }$ de $i$ :

$\mathrm{d}{n}_i=\mathrm{d}{n}_i^{\alpha }+\mathrm{d}{n}_i^{\beta }=0$

$\mathrm{d}{n}_i^{\beta }=-\mathrm{d}{n}_i^{\alpha }$

On obtient finalement :

$T\mathrm{d}S=-({\mu }_i^{\alpha }-{\mu }_i^{\beta })\mathrm{d}{n}_i^{\alpha }$

On suppose que le potentiel chimique de l'espèce $i$ est plus important dans le compartiment $\alpha$ que dans le compartiment $\beta$, soit ${\mu }_i^{\alpha }>{\mu }_i^{\beta }$. La température absolue $T$ étant toujours positive et, selon le deuxième principe de la thermodynamique, l'entropie du système ne pouvant qu'augmenter, soit $\mathrm{d}S>0$, alors $\mathrm{d}{n}_i^{\alpha }<0$, d'où $\mathrm{d}{n}_i^{\beta }>0$. En conséquence, le corps $i$ migre du compartiment $\alpha$, dans lequel son potentiel est le plus important, vers le compartiment $\beta$, de moindre potentiel.

La migration cesse lorsque les potentiels des deux compartiments sont égaux, soit ${\mu }_i^{\alpha }={\mu }_i^{\beta }$ : l'entropie a alors atteint un maximum, soit $\mathrm{d}S=0$.

Définitions basées sur le potentiel chimique

Solution idéale

Une solution idéale est définie par la relation :

Solution idéale : ${\mu }_i^{\phi ,\text{id}}(P,T,x)={\mu }_i^{\phi ,*}(P,T)+RT\ln x_i$

avec :

• $x_i$ la fraction molaire du constituant $i$ dans la solution ;
• ${\mu }_i^{\phi ,\text{id}}$ le potentiel chimique du composant $i$ dans la solution idéale dans la phase $\phi$ ;
• ${\mu }_i^{\phi ,*}$ le potentiel chimique du composant $i$ pur, aux mêmes $P$, $T$ et phase $\phi$ que la solution idéale ;
• $R$ la constante universelle des gaz parfaits.

Le concept de solution idéale lie les propriétés des corps purs à celles d'un mélange dans lequel les molécules des diverses espèces interagissent avec les molécules des autres espèces de la même façon que ces mêmes molécules interagissent avec les molécules de même espèce à l'état de corps purs. Un mélange de gaz parfaits est un exemple de solution idéale, selon le théorème de Gibbs.

Cette relation permet également de montrer que :

• pour le corps $i$ pur : $\underset{x_i\to 1}{\lim }{\mu }_i={\mu }_i^{*}$ ;
• pour le corps $i$ infiniment dilué dans un solvant $s$ : $\underset{x_i\to 0}{\lim }{\mu }_i=-\mathrm{\infty }$.

Fugacité et coefficient de fugacité

En 1900 et 1901, Gilbert Lewis introduisit la notion de fugacité d'un corps $i$, notée $f_i$, qu'il définit ainsi en 1923Modèle:, :

Variation isotherme du potentiel chimique : $\mathrm{d}{\mu }_i=RT\mathrm{d}\ln f_i$

La fugacité a la dimension d'une pression.

En intégrant entre l'état de gaz parfait pur et l'état réel, les deux états étant aux mêmes pression et température :

${\mu }_i-{\mu }_i^{\bullet ,*}=RT\underset{\ln f_i^{\bullet ,*}}{\overset{\ln f_i}{\int }}\mathrm{d}\ln f_i=RT\ln \left(\frac{f_i}{f_i^{\bullet ,*}}\right)$

La pression étant la fugacité du gaz parfait pur : $f_i^{\bullet ,*}=P$, la fugacité exprime l'écart entre le potentiel chimique ${\mu }_i$ d'un corps $i$ dans un mélange réel et le potentiel chimique ${\mu }_i^{\bullet ,*}$ du même corps à l'état de gaz parfait pur aux mêmes pression et température que le mélange réel :

Fugacité : ${\mu }_i-{\mu }_i^{\bullet ,*}=RT\ln \left(\frac{f_i}{P}\right)$

Un mélange de gaz parfaits étant une solution idéale selon le théorème de Gibbs, on a la relation :

${\mu }_i^{\bullet }={\mu }_i^{\bullet ,*}+RT\ln x_i$

en introduisant le potentiel chimique ${\mu }_i^{\bullet }$ du même corps $i$ dans un mélange de gaz parfaits aux mêmes pression, température et composition (fractions molaires $x_i$) que le mélange réel. On a donc, en substituant ${\mu }_i^{\bullet }$ à ${\mu }_i^{\bullet ,*}$ dans la relation obtenue précédemment pour ${\mu }_i$:

${\mu }_i-{\mu }_i^{\bullet }=RT\underset{\ln f_i^{\bullet }}{\overset{\ln f_i}{\int }}\mathrm{d}\ln f_i=RT\ln \left(\frac{f_i}{f_i^{\bullet }}\right)=RT\ln \left(\frac{f_i}{x_{i}P}\right)=RT\ln {\varphi }_i$

Le rapport entre la fugacité $f_i$ réelle et la fugacité $f_i^{\bullet }=x_{i}P$ du gaz parfait en mélange (qui est égale à sa pression partielle) est appelé coefficient de fugacité, il est noté ${\varphi }_i$ :

Coefficient de fugacité ${\mu }_i-{\mu }_i^{\bullet }=RT\ln {\varphi }_i$ ${\varphi }_i=\frac{f_i}{f_i^{\bullet }}=\frac{f_i}{x_{i}P}$

Le coefficient de fugacité est adimensionnel.

Le coefficient de fugacité peut être établi pour n'importe quelle phase (gaz, liquide, solide), pourvu que l'on dispose d'une équation d'état permettant de décrire les propriétés de cette phase. Dans la pratique, le coefficient de fugacité est surtout utilisé pour représenter des phases gazeuses.

Activité chimique et coefficient d'activité

L'activité chimique est une notion introduite par Lewis en 1923. Elle est définie en intégrant la relation qui lie la variation isotherme du potentiel chimique et celle de la fugacité entre un état standard et l'état réel d'un corps $i$ dans un mélange, les états réel et standard étant à la même température :

${\mu }_i-{\mu }_i^{\circ }=RT\underset{\ln f_i^{\circ }}{\overset{\ln f_i}{\int }}\mathrm{d}\ln f_i=RT\ln \left(\frac{f_i}{f_i^{\circ }}\right)=RT\ln a_i$

avec :

• ${\mu }_i$ le potentiel chimique du corps $i$ dans la solution réelle ;
• ${\mu }_i^{\circ }$ le potentiel chimique du corps $i$ dans l'état standard à la même température que l'état réel ;
• $f_i$ la fugacité du corps $i$ dans le mélange réel ;
• $f_i^{\circ }$ la fugacité du corps $i$ dans l'état standard à la même température que l'état réel ;
• $a_i$ l'activité chimique du corps $i$ entre le mélange réel et l'état standard à la même température que l'état réel.

L'activité chimique est définie par le rapport des fugacités réelle et dans l'état standard :

Activité ${\mu }_i-{\mu }_i^{\circ }=RT\ln a_i$ $a_i=\frac{f_i}{f_i^{\circ }}$

L'activité est adimensionnelle.

La notion d'activité est surtout employée pour exprimer l'écart entre le potentiel chimique ${\mu }_i$ d'un corps $i$ dans un mélange réel et le potentiel chimique ${\mu }_i^{\text{id}}$ du même corps dans une solution idéale aux mêmes pression, température et composition (fractions molaires $x_i$), et dans la même phase que le mélange réel :

${\mu }_i-{\mu }_i^{\text{id}}=RT\underset{\ln f_i^{\text{id}}}{\overset{\ln f_i}{\int }}\mathrm{d}\ln f_i=RT\ln \left(\frac{f_i}{f_i^{\text{id}}}\right)=RT\ln \left(\frac{a_i}{x_i}\right)=RT\ln {\gamma }_i$

avec $f_i^{\text{id}}=x_i\cdot f_i^{*}$ la fugacité en solution idéale. L'état standard est le corps $i$ pur aux mêmes pression et température, dans la même phase que la solution réelle, dont la fugacité est $f_i^{*}$. Le rapport entre la fugacité réelle et la fugacité en solution idéale est appelé coefficient d'activité, il est noté ${\gamma }_i$ :

Coefficient d'activité ${\mu }_i-{\mu }_i^{\text{id}}=RT\ln {\gamma }_i$ ${\gamma }_i=\frac{f_i}{f_i^{\text{id}}}=\frac{f_i}{x_i\cdot f_i^{*}}=\frac{a_i}{x_i}$

Le coefficient d'activité est adimensionnel.

Le coefficient d'activité peut être établi pour n'importe quelle phase (gaz, liquide, solide), pourvu que l'on dispose d'un modèle permettant de décrire les propriétés de cette phase. Dans la pratique, le coefficient d'activité est surtout utilisé pour représenter les phases condensées (liquide et solide), avec comme état de référence le corps $i$ pur aux mêmes pression et température, et dans la même phase que le mélange réel.

Notes et références

Bibliographie

• J. Willard Gibbs, traduit par Henry Le Chatelier, « Équilibre des systèmes chimiques », Modèle:Éd. G. Carré et C. Naud (Paris), 1899, Accessible sur Gallica.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.
• Modèle:Ouvrage.

Articles connexes

• Activité chimique
• Affinité chimique
• Enthalpie libre
• Équilibre chimique
• Équilibre de phases
• Équilibre thermodynamique
• Fugacité
• Potentiel thermodynamique
• Solution idéale