Relation de Gibbs-Duhem

La relation de Gibbs-Duhem, du nom de Willard Gibbs et Pierre Duhem, est une relation de thermodynamique qui s'applique à tout système thermodynamique à l'équilibre contenant plusieurs espèces chimiques en une seule phase. Lorsque l'on parle de la relation de Gibbs-Duhem, sans autre précision, il s'agit de la relation liée à l'enthalpie libre qui s'écrit sous la forme :

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

avec :

  • n i {\displaystyle n_{i}} la quantité, ou nombre de moles, de l'espèce i {\displaystyle i} au sein du mélange ;
  • μ i {\displaystyle \mu _{i}} le potentiel chimique de l'espèce i {\displaystyle i} au sein du mélange ;
  • V {\displaystyle V} le volume du mélange ;
  • P {\displaystyle P} la pression du mélange ;
  • S {\displaystyle S} l'entropie du mélange ;
  • T {\displaystyle T} la température du mélange ;
  • N {\displaystyle N} le nombre d'espèces chimiques dans le mélange.

À pression et température constantes, cette relation devient :

À pression et température constantes : i = 1 N n i d μ i = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}=0}

De ces premières relations basées sur l'enthalpie libre découlent d'autres relations impliquant les fugacités, coefficients de fugacité, activités chimiques, coefficients d'activité d'un mélange. La relation de Gibbs-Duhem peut être généralisée à toute grandeur extensive autre que l'enthalpie libre.

Démonstration

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

μ i = G ¯ i = ( G n i ) P , T , n j i {\displaystyle \mu _{i}={\bar {G}}_{i}=\left({\frac {\partial G}{\partial n_{i}}}\right)_{P,T,n_{j\neq i}}}

Selon le théorème d'Euler sur les fonctions homogènes du premier ordre, on peut écrire pour l'enthalpie libre, qui est une grandeur extensive :

G = i = 1 N n i G ¯ i {\displaystyle G=\sum _{i=1}^{N}n_{i}{\bar {G}}_{i}}

et donc :

G = i = 1 N n i μ i {\displaystyle G=\sum _{i=1}^{N}n_{i}\mu _{i}}

En différentiant l'identité d'Euler pour l'enthalpie libre, nous obtenons :

(a) d G = i = 1 N n i d μ i + i = 1 N μ i d n i {\displaystyle \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 :

(b) d G = V d P S d T + i = 1 N μ i d n i {\displaystyle \mathrm {d} G=V\,\mathrm {d} P-S\,\mathrm {d} T+\sum _{i=1}^{N}\mu _{i}\,\mathrm {d} n_{i}}

Nous pouvons identifier les différents termes des équations (a) et (b), on obtient la relation de Gibbs-Duhem :

V d P S d T = i = 1 N n i d μ i {\displaystyle V\,\mathrm {d} P-S\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}}

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

V ¯ d P S ¯ d T = i = 1 N x i d μ i {\displaystyle {\bar {V}}\,\mathrm {d} P-{\bar {S}}\,\mathrm {d} T=\sum _{i=1}^{N}x_{i}\,\mathrm {d} \mu _{i}}

avec :

  • V ¯ = V / n {\displaystyle {\bar {V}}=V/n} le volume molaire du mélange ;
  • S ¯ = S / n {\displaystyle {\bar {S}}=S/n} l'entropie molaire du mélange ;
  • x i = n i / n {\displaystyle x_{i}=n_{i}/n} la fraction molaire de l'espèce i {\displaystyle i} dans le mélange.

Relations dérivées

Représentation énergie et représentation entropie

La relation définie précédemment est appelée relation de Gibbs-Duhem en représentation énergie :

Relation de Gibbs-Duhem en représentation énergie : S d T V d P + i = 1 N n i d μ i = 0 {\displaystyle S\,\mathrm {d} T-V\,\mathrm {d} P+\sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}=0}

Il existe une seconde relation de Gibbs-Duhem basée sur l'entropie, que l'on obtient selon :

G = U + P V T S {\displaystyle G=U+PV-TS}
S = U T + P V T G T {\displaystyle S={U \over T}+{PV \over T}-{G \over T}}

avec G = i = 1 N μ i n i {\displaystyle G=\sum _{i=1}^{N}\mu _{i}n_{i}} , en différentiant :

(c) d S = d U T + U d ( 1 T ) + P T d V + V d ( P T ) i = 1 N μ i T d n i i = 1 N n i d ( μ i T ) {\displaystyle \mathrm {d} S={\mathrm {d} U \over T}+U\,\mathrm {d} \!\left({1 \over T}\right)+{P \over T}\,\mathrm {d} V+V\,\mathrm {d} \!\left({P \over T}\right)-\sum _{i=1}^{N}{\mu _{i} \over T}\,\mathrm {d} n_{i}-\sum _{i=1}^{N}n_{i}\,\mathrm {d} \!\left({\mu _{i} \over T}\right)}

or la différentielle de l'énergie interne U {\displaystyle U} donne :

d U = P d V + T d S + i = 1 N μ i d n i {\displaystyle \mathrm {d} U=-P\,\mathrm {d} V+T\,\mathrm {d} S+\sum _{i=1}^{N}\mu _{i}\,\mathrm {d} n_{i}}
(d) d S = d U T + P T d V i = 1 N μ i T d n i {\displaystyle \mathrm {d} S={\mathrm {d} U \over T}+{P \over T}\,\mathrm {d} V-\sum _{i=1}^{N}{\mu _{i} \over T}\,\mathrm {d} n_{i}}

Nous pouvons identifier les différents termes des équations (c) et (d), on obtient la relation de Gibbs-Duhem en représentation entropie[1] :

Relation de Gibbs-Duhem en représentation entropie : U d ( 1 T ) + V d ( P T ) i = 1 N n i d ( μ i T ) = 0 {\displaystyle U\,\mathrm {d} \!\left({1 \over T}\right)+V\,\mathrm {d} \!\left({P \over T}\right)-\sum _{i=1}^{N}n_{i}\,\mathrm {d} \!\left({\mu _{i} \over T}\right)=0}

ou encore :

H d ( 1 T ) + V T d P i = 1 N n i d ( μ i T ) = 0 {\displaystyle H\,\mathrm {d} \!\left({1 \over T}\right)+{V \over T}\,\mathrm {d} \!P-\sum _{i=1}^{N}n_{i}\,\mathrm {d} \!\left({\mu _{i} \over T}\right)=0}

avec H = U + P V {\displaystyle H=U+PV} l'enthalpie.

Relation de Gibbs-Duhem générale

Toute grandeur extensive X {\displaystyle X} peut être exprimée en fonction de la pression P {\displaystyle P} , de la température T {\displaystyle T} et des quantités de matière n i {\displaystyle n_{i}}  : X = X ( P , T , n ) {\displaystyle X=X\!\left(P,T,n\right)} , même si P {\displaystyle P} et T {\displaystyle T} ne sont pas ses variables naturelles. Aussi peut-on écrire la différentielle totale de toute grandeur extensive sous la forme :

d X = ( X P ) T , n d P + ( X T ) P , n d T + i = 1 N X ¯ i d n i {\displaystyle \mathrm {d} X=\left({\partial X \over \partial P}\right)_{T,n}\,\mathrm {d} P+\left({\partial X \over \partial T}\right)_{P,n}\,\mathrm {d} T+\sum _{i=1}^{N}{\bar {X}}_{i}\,\mathrm {d} n_{i}}

avec X ¯ i = ( X n i ) P , T , n j i {\displaystyle {\bar {X}}_{i}=\left({\partial X \over \partial n_{i}}\right)_{P,T,n_{j\neq i}}} la grandeur molaire partielle de l'espèce i {\displaystyle i} .

Le théorème d'Euler sur les fonctions homogènes du premier ordre lie toute grandeur extensive X {\displaystyle X} aux grandeurs molaires partielles X ¯ i {\displaystyle {\bar {X}}_{i}} , définies aux mêmes P {\displaystyle P} , T {\displaystyle T} et composition, selon :

X = i = 1 N n i X ¯ i {\displaystyle X=\sum _{i=1}^{N}n_{i}{\bar {X}}_{i}}

En différentiant cette expression il vient :

d X = i = 1 N n i d X ¯ i + i = 1 N X ¯ i d n i {\displaystyle \mathrm {d} X=\sum _{i=1}^{N}n_{i}\,\mathrm {d} {\bar {X}}_{i}+\sum _{i=1}^{N}{\bar {X}}_{i}\,\mathrm {d} n_{i}}

En identifiant terme à terme les deux expressions de d X {\displaystyle \mathrm {d} X} , nous obtenons la relation de Gibbs-Duhem générale :

Relation de Gibbs-Duhem générale : ( X P ) T , n d P + ( X T ) P , n d T = i = 1 N n i d X ¯ i {\displaystyle \left({\partial X \over \partial P}\right)_{T,n}\,\mathrm {d} P+\left({\partial X \over \partial T}\right)_{P,n}\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} {\bar {X}}_{i}}
Exemple
Pour le volume V {\displaystyle V} d'un mélange on pourra écrire :
( V P ) T , n d P + ( V T ) P , n d T = i = 1 N n i d V ¯ i {\displaystyle \left({\partial V \over \partial P}\right)_{T,n}\,\mathrm {d} P+\left({\partial V \over \partial T}\right)_{P,n}\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} {\bar {V}}_{i}}
Avec :
  • α = 1 V ( V T ) P , n {\displaystyle \alpha ={1 \over V}\left({\partial V \over \partial T}\right)_{P,n}} le coefficient de dilatation isobare ;
  • χ T = 1 V ( V P ) T , n {\displaystyle \chi _{T}=-{1 \over V}\left({\partial V \over \partial P}\right)_{T,n}} la compressibilité isotherme ;
  • V ¯ i = ( V n i ) P , T , n j i {\displaystyle {\bar {V}}_{i}=\left({\partial V \over \partial n_{i}}\right)_{P,T,n_{j\neq i}}} le volume molaire partiel de l'espèce i {\displaystyle i}  ;
on obtient :
V χ T d P + V α d T = i = 1 N n i d V ¯ i {\displaystyle -V\chi _{T}\,\mathrm {d} P+V\alpha \,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} {\bar {V}}_{i}}

Solutions idéales et mélanges de gaz parfaits

Déclinée au cas d'une solution idéale, la relation de Gibbs-Duhem donne :

Relation de Gibbs-Duhem pour une solution idéale : V id d P S id d T = i = 1 N n i d μ i id {\displaystyle V^{\text{id}}\,\mathrm {d} P-S^{\text{id}}\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{\text{id}}}

en particulier, un mélange de gaz parfaits étant une solution idéale selon le théorème de Gibbs :

Relation de Gibbs-Duhem pour un mélange de gaz parfaits : V d P S d T = i = 1 N n i d μ i {\displaystyle V^{\bullet }\,\mathrm {d} P-S^{\bullet }\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{\bullet }}

avec :

  • V id {\displaystyle V^{\text{id}}} le volume de la solution idéale ;
  • V {\displaystyle V^{\bullet }} le volume du mélange de gaz parfaits ;
  • S id {\displaystyle S^{\text{id}}} l'entropie de la solution idéale ;
  • S {\displaystyle S^{\bullet }} l'entropie du mélange de gaz parfaits ;
  • μ i id {\displaystyle \mu _{i}^{\text{id}}} le potentiel chimique du corps i {\displaystyle i} dans la solution idéale ;
  • μ i {\displaystyle \mu _{i}^{\bullet }} le potentiel chimique du corps i {\displaystyle i} dans le mélange de gaz parfaits.

Dans une solution idéale, le potentiel chimique μ i {\displaystyle \mu _{i}} de toute espèce chimique i {\displaystyle i} est liée au potentiel chimique μ i {\displaystyle \mu _{i}^{*}} de cette même espèce pure aux mêmes pression et température et dans la même phase par la fraction molaire x i {\displaystyle x_{i}} de cette espèce selon :

μ i id = μ i + R T ln x i {\displaystyle \mu _{i}^{\text{id}}=\mu _{i}^{*}+RT\ln x_{i}}

Le potentiel chimique μ i {\displaystyle \mu _{i}^{*}} de toute espèce pure étant identique à l'enthalpie libre molaire G ¯ i {\displaystyle {\bar {G}}_{i}^{*}} , il ne dépend que de la pression et de la température :

μ i = G ¯ i = μ i ( P , T ) {\displaystyle \mu _{i}^{*}={\bar {G}}_{i}^{*}=\mu _{i}^{*}\left(P,T\right)}

à pression et température constantes :

d μ i = 0 {\displaystyle \mathrm {d} \mu _{i}^{*}=0}

La relation de Gibbs-Duhem à pression et température constantes donne :

i = 1 N n i d μ i id = i = 1 N n i d μ i = 0  à P et T constantes + i = 1 N n i R T d ln x i = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{\text{id}}=\underbrace {\sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{*}} _{=0{\text{ à P et T}} \atop {\text{constantes}}}+\sum _{i=1}^{N}n_{i}RT\,\mathrm {d} \ln x_{i}=0}
i = 1 N n i R T d ln x i = 0 {\displaystyle \sum _{i=1}^{N}n_{i}RT\,\mathrm {d} \ln x_{i}=0}

En divisant par n R T {\displaystyle nRT} on obtient la relation entre fractions molaires :

À pression et température constantes : i = 1 N x i d ln x i = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln x_{i}=0}

soit :

i = 1 N x i d ln x i = i = 1 N x i d x i x i = i = 1 N d x i = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln x_{i}=\sum _{i=1}^{N}x_{i}\,{\mathrm {d} x_{i} \over x_{i}}=\sum _{i=1}^{N}\,\mathrm {d} x_{i}=0}

Cette relation découle directement de la relation entre fractions molaires :

i = 1 N x i = 1 {\displaystyle \sum _{i=1}^{N}x_{i}=1}

Elle implique également que, pour tout corps j {\displaystyle j} du mélange :

i = 1 N n i ( ln x i n j ) P , T , n k j = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\left({\frac {\partial \ln x_{i}}{\partial n_{j}}}\right)_{P,T,n_{k\neq j}}=0}

Activités chimiques

En soustrayant la relation des solutions idéales à la relation des mélanges réels, toutes les grandeurs étant définies aux mêmes pression, température et composition que le mélange réel, on a :

( V V id ) d P ( S S id ) d T = i = 1 N n i d ( μ i μ i id ) {\displaystyle \left(V-V^{\text{id}}\right)\,\mathrm {d} P-\left(S-S^{\text{id}}\right)\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} \!\left(\mu _{i}-\mu _{i}^{\text{id}}\right)}

en identifiant les grandeurs d'excès on obtient :

Relation de Gibbs-Duhem pour les grandeurs d'excès : V E d P S E d T = i = 1 N n i d μ i E {\displaystyle V^{\text{E}}\,\mathrm {d} P-S^{\text{E}}\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{\text{E}}}

avec :

  • V E = V V id {\displaystyle V^{\text{E}}=V-V^{\text{id}}} le volume d'excès ;
  • S E = S S id {\displaystyle S^{\text{E}}=S-S^{\text{id}}} l'entropie d'excès ;
  • μ i E = μ i μ i id {\displaystyle \mu _{i}^{\text{E}}=\mu _{i}-\mu _{i}^{\text{id}}} le potentiel chimique d'excès du corps i {\displaystyle i} .

En introduisant le coefficient d'activité γ i {\displaystyle \gamma _{i}} de chaque espèce, défini par :

μ i E = μ i μ i id = R T ln γ i {\displaystyle \mu _{i}^{\text{E}}=\mu _{i}-\mu _{i}^{\text{id}}=RT\ln \gamma _{i}}

la relation précédente devient, à pression et température constantes :

i = 1 N n i d μ i E = i = 1 N n i R T d ln γ i = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{\text{E}}=\sum _{i=1}^{N}n_{i}RT\,\mathrm {d} \ln \gamma _{i}=0}

En divisant par n R T {\displaystyle nRT} on obtient la relation entre coefficients d'activité :

À pression et température constantes : i = 1 N x i d ln γ i = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln \gamma _{i}=0}

Cette relation implique que, pour tout corps j {\displaystyle j} du mélange :

i = 1 N n i ( ln γ i n j ) P , T , n k j = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\left({\frac {\partial \ln \gamma _{i}}{\partial n_{j}}}\right)_{P,T,n_{k\neq j}}=0}

On a également, en considérant l'activité chimique a i = x i γ i {\displaystyle a_{i}=x_{i}\gamma _{i}}  :

i = 1 N x i d ln γ i = i = 1 N x i d ln ( a i x i ) = i = 1 N x i d ln a i i = 1 N x i d ln x i = 0 = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln \gamma _{i}=\sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln \!\left({a_{i} \over x_{i}}\right)=\sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln a_{i}-\underbrace {\sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln x_{i}} _{=0}=0}

d'où la relation entre activités :

À pression et température constantes : i = 1 N x i d ln a i = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln a_{i}=0}

Cette relation implique que, pour tout corps j {\displaystyle j} du mélange :

i = 1 N n i ( ln a i n j ) P , T , n k j = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\left({\frac {\partial \ln a_{i}}{\partial n_{j}}}\right)_{P,T,n_{k\neq j}}=0}

Fugacités

En prenant comme solution idéale de référence le mélange de gaz parfaits aux mêmes pression, température et composition que le mélange réel, on a :

( V V ) d P ( S S ) d T = i = 1 N n i d ( μ i μ i ) {\displaystyle \left(V-V^{\bullet }\right)\,\mathrm {d} P-\left(S-S^{\bullet }\right)\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} \!\left(\mu _{i}-\mu _{i}^{\bullet }\right)}

en identifiant les grandeurs résiduelles on obtient :

Relation de Gibbs-Duhem pour les grandeurs résiduelles : V RES d P S RES d T = i = 1 N n i d μ i RES {\displaystyle V^{\text{RES}}\,\mathrm {d} P-S^{\text{RES}}\,\mathrm {d} T=\sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{\text{RES}}}

avec :

  • V RES = V V {\displaystyle V^{\text{RES}}=V-V^{\bullet }} le volume résiduel ;
  • S RES = S S {\displaystyle S^{\text{RES}}=S-S^{\bullet }} l'entropie résiduelle ;
  • μ i RES = μ i μ i {\displaystyle \mu _{i}^{\text{RES}}=\mu _{i}-\mu _{i}^{\bullet }} le potentiel chimique résiduel du corps i {\displaystyle i} .

En introduisant le coefficient de fugacité ϕ i {\displaystyle \phi _{i}} de chaque espèce, défini par :

μ i RES = μ i μ i = R T ln ϕ i {\displaystyle \mu _{i}^{\text{RES}}=\mu _{i}-\mu _{i}^{\bullet }=RT\,\ln \phi _{i}}

la relation précédente devient, à pression et température constantes :

i = 1 N n i d μ i RES = i = 1 N n i R T d ln ϕ i = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\,\mathrm {d} \mu _{i}^{\text{RES}}=\sum _{i=1}^{N}n_{i}RT\,\mathrm {d} \ln \phi _{i}=0}

En divisant par n R T {\displaystyle nRT} on obtient la relation entre coefficients de fugacité :

À pression et température constantes : i = 1 N x i d ln ϕ i = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln \phi _{i}=0}

Cette relation implique que, pour tout corps j {\displaystyle j} du mélange :

i = 1 N n i ( ln ϕ i n j ) P , T , n k j = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\left({\frac {\partial \ln \phi _{i}}{\partial n_{j}}}\right)_{P,T,n_{k\neq j}}=0}

On a également, en considérant la fugacité f i = x i ϕ i P {\displaystyle f_{i}=x_{i}\phi _{i}P}  :

i = 1 N x i d ln ϕ i = i = 1 N x i d ln ( f i x i P ) = i = 1 N x i d ln f i i = 1 N x i d ln x i = 0 i = 1 N x i d ln P = 0 à pression constante = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln \phi _{i}=\sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln \!\left({f_{i} \over x_{i}P}\right)=\sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln f_{i}-\underbrace {\sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln x_{i}} _{=0}-\underbrace {\sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln P} _{=0 \atop {\text{à pression constante}}}=0}

On obtient la relation entre fugacités :

À pression et température constantes : i = 1 N x i d ln f i = 0 {\displaystyle \sum _{i=1}^{N}x_{i}\,\mathrm {d} \ln f_{i}=0}

Cette relation implique que, pour tout corps j {\displaystyle j} du mélange :

i = 1 N n i ( ln f i n j ) P , T , n k j = 0 {\displaystyle \sum _{i=1}^{N}n_{i}\left({\frac {\partial \ln f_{i}}{\partial n_{j}}}\right)_{P,T,n_{k\neq j}}=0}

Exemples d'application

Lois des équilibres de phases

La relation de Gibbs-Duhem intervient dans l'établissement théorique de nombreuses formules décrivant les équilibres de phases, elle entre au nombre des relations prises en compte pour établir la règle des phases.

La formule de Clapeyron décrit l'évolution de la pression de changement d'état d'un corps pur en fonction de la température. Cette formule peut être étendue aux mélanges.

Plusieurs lois décrivent les propriétés colligatives d'une solution liquide, c'est-à-dire la différence entre une propriété donnée d'un solvant pur liquide et la même propriété de ce solvant en présence d'un soluté très dilué :

La loi de Raoult concerne les équilibres liquide-vapeur idéaux. La relation de Duhem-Margules permet de démontrer les règles de Konovalov concernant la dépendance des pressions partielles d'une phase vapeur à la composition du liquide en équilibre à température constante.

L'équation de Schröder-van Laar concerne les équilibres liquide-solide.

Le théorème de Gibbs-Konovalov concerne les azéotropes et les points de fusion congruents.

Détermination des grandeurs molaires partielles, méthode de Roozeboom

Article détaillé : Méthode de Roozeboom.

Détermination graphique des potentiels chimiques

Pour un mélange binaire, ne contenant que deux espèces chimiques, la relation de Gibbs-Duhem à pression et température constantes permet d'établir la variation de l'enthalpie libre molaire G ¯ {\displaystyle {\bar {G}}} en fonction de la fraction molaire x 1 {\displaystyle x_{1}} de l'espèce 1 à P {\displaystyle P} et T {\displaystyle T} constantes :

Variation de l'enthalpie libre molaire d'un mélange binaire avec la composition : ( G ¯ x 1 ) P , T = μ 1 μ 2 {\displaystyle \left({\partial {\bar {G}} \over \partial x_{1}}\right)_{P,T}=\mu _{1}-\mu _{2}}
Démonstration

Pour un mélange binaire la relation de Gibbs-Duhem à pression et température constantes se réduit à :

n 1 d μ 1 + n 2 d μ 2 = 0 {\displaystyle n_{1}\,\mathrm {d} \mu _{1}+n_{2}\,\mathrm {d} \mu _{2}=0}

avec :

  • n 1 {\displaystyle n_{1}} et n 2 {\displaystyle n_{2}} les quantités respectives des espèces 1 et 2 ;
  • μ 1 {\displaystyle \mu _{1}} et μ 2 {\displaystyle \mu _{2}} les potentiels chimiques respectifs des espèces 1 et 2.

Soit n = n 1 + n 2 {\displaystyle n=n_{1}+n_{2}} la quantité totale de matière dans le mélange. On note x 1 = n 1 n {\displaystyle x_{1}={n_{1} \over n}} et x 2 = n 2 n {\displaystyle x_{2}={n_{2} \over n}} les fractions molaires respectives des espèces 1 et 2. En divisant la relation de Gibbs-Duhem par n {\displaystyle n} , on a :

x 1 d μ 1 + x 2 d μ 2 = 0 {\displaystyle x_{1}\,\mathrm {d} \mu _{1}+x_{2}\,\mathrm {d} \mu _{2}=0}

Puisque les fractions molaires sont liées par la relation :

(1) x 1 + x 2 = 1 {\displaystyle x_{1}+x_{2}=1}

on réécrit :

(2) x 1 d μ 1 + ( 1 x 1 ) d μ 2 = 0 {\displaystyle x_{1}\,\mathrm {d} \mu _{1}+\left(1-x_{1}\right)\,\mathrm {d} \mu _{2}=0}

Les potentiels chimiques dépendent de la pression, de la température et des fractions molaires des espèces présentes dans le mélange :

μ 1 = μ 1 ( P , T , x 1 , x 2 ) {\displaystyle \mu _{1}=\mu _{1}\left(P,T,x_{1},x_{2}\right)}
μ 2 = μ 2 ( P , T , x 1 , x 2 ) {\displaystyle \mu _{2}=\mu _{2}\left(P,T,x_{1},x_{2}\right)}

Étant donné la relation (1), les potentiels chimiques d'un mélange binaire ne dépendent que d'une seule variable de composition, nous choisissons x 1 {\displaystyle x_{1}} . Nous avons donc à pression et température constantes :

d μ 1 = ( μ 1 x 1 ) P , T d x 1 {\displaystyle \mathrm {d} \mu _{1}=\left({\partial \mu _{1} \over \partial x_{1}}\right)_{P,T}\,\mathrm {d} x_{1}}
d μ 2 = ( μ 2 x 1 ) P , T d x 1 {\displaystyle \mathrm {d} \mu _{2}=\left({\partial \mu _{2} \over \partial x_{1}}\right)_{P,T}\,\mathrm {d} x_{1}}

En substituant ces relations dans l'équation (2), on obtient :

( x 1 ( μ 1 x 1 ) P , T + ( 1 x 1 ) ( μ 2 x 1 ) P , T ) d x 1 = 0 {\displaystyle \left(x_{1}\left({\partial \mu _{1} \over \partial x_{1}}\right)_{P,T}+\left(1-x_{1}\right)\left({\partial \mu _{2} \over \partial x_{1}}\right)_{P,T}\right)\,\mathrm {d} x_{1}=0}

La variation d x 1 {\displaystyle \mathrm {d} x_{1}} étant quelconque, ceci implique que :

(3) x 1 ( μ 1 x 1 ) P , T + ( 1 x 1 ) ( μ 2 x 1 ) P , T = 0 {\displaystyle x_{1}\left({\partial \mu _{1} \over \partial x_{1}}\right)_{P,T}+\left(1-x_{1}\right)\left({\partial \mu _{2} \over \partial x_{1}}\right)_{P,T}=0}

D'autre part, selon le théorème d'Euler sur les fonctions homogènes du premier ordre, on peut écrire, étant donné l'égalité des potentiels chimiques et des enthalpies libres molaires partielles :

G = n 1 μ 1 + n 2 μ 2 {\displaystyle G=n_{1}\mu _{1}+n_{2}\mu _{2}}

En divisant par la quantité totale de matière n {\displaystyle n} , on obtient, en introduisant l'enthalpie libre molaire G ¯ {\displaystyle {\bar {G}}}  :

G ¯ = G n = x 1 μ 1 + x 2 μ 2 {\displaystyle {\bar {G}}={G \over n}=x_{1}\mu _{1}+x_{2}\mu _{2}}

En considérant la relation (1), on réécrit :

(4) G ¯ = x 1 μ 1 + ( 1 x 1 ) μ 2 {\displaystyle {\bar {G}}=x_{1}\mu _{1}+\left(1-x_{1}\right)\mu _{2}}

De même que les potentiels chimiques, l'enthalpie libre molaire dépend de la pression, de la température et des fractions molaires des espèces présentes dans le mélange :

G ¯ = G ¯ ( P , T , x 1 , x 2 ) {\displaystyle {\bar {G}}={\bar {G}}\left(P,T,x_{1},x_{2}\right)}

Étant donné la relation (1), l'enthalpie libre d'un mélange binaire ne dépend que d'une seule variable de composition, nous choisissons x 1 {\displaystyle x_{1}} . On a, en dérivant la relation (4) :

( G ¯ x 1 ) P , T = μ 1 + x 1 ( μ 1 x 1 ) P , T μ 2 + ( 1 x 1 ) ( μ 2 x 1 ) P , T {\displaystyle \left({\partial {\bar {G}} \over \partial x_{1}}\right)_{P,T}=\mu _{1}+x_{1}\left({\partial \mu _{1} \over \partial x_{1}}\right)_{P,T}-\mu _{2}+\left(1-x_{1}\right)\left({\partial \mu _{2} \over \partial x_{1}}\right)_{P,T}}

En considérant la relation (3), il vient finalement :

( G ¯ x 1 ) P , T = μ 1 μ 2 {\displaystyle \left({\partial {\bar {G}} \over \partial x_{1}}\right)_{P,T}=\mu _{1}-\mu _{2}}

Si l'on parvient à déterminer expérimentalement l'enthalpie libre d'un mélange binaire à pression et température données en fonction de la composition, pour déterminer les potentiels chimiques des deux espèces on trace la courbe de l'enthalpie libre molaire G ¯ {\displaystyle {\bar {G}}} en fonction de la fraction molaire x 1 {\displaystyle x_{1}} . Cette courbe admet deux limites selon la relation (4) :

  • pour x 1 = 0 {\displaystyle x_{1}=0} l'espèce 2 est pure et G ¯ = μ 2 {\displaystyle {\bar {G}}=\mu _{2}^{*}} le potentiel chimique de l'espèce 2 pure à P {\displaystyle P} et T {\displaystyle T}  ;
  • pour x 1 = 1 {\displaystyle x_{1}=1} l'espèce 1 est pure et G ¯ = μ 1 {\displaystyle {\bar {G}}=\mu _{1}^{*}} le potentiel chimique de l'espèce 1 pure à P {\displaystyle P} et T {\displaystyle T} .

La pente de la tangente à la courbe ainsi obtenue en un point quelconque de composition x 1 {\displaystyle x_{1}} vaut μ 1 μ 2 {\displaystyle \mu _{1}-\mu _{2}} (les deux potentiels dépendant de x 1 {\displaystyle x_{1}} ). Cette tangente a pour équation :

g ( z ) = ( μ 1 μ 2 ) ( z x 1 ) + G ¯ ( x 1 ) {\displaystyle g\left(z\right)=\left(\mu _{1}-\mu _{2}\right)\left(z-x_{1}\right)+{\bar {G}}\left(x_{1}\right)}

soit, étant donné la relation (4) :

g ( z ) = ( μ 1 μ 2 ) ( z x 1 ) + x 1 μ 1 + ( 1 x 1 ) μ 2 {\displaystyle g\left(z\right)=\left(\mu _{1}-\mu _{2}\right)\left(z-x_{1}\right)+x_{1}\mu _{1}+\left(1-x_{1}\right)\mu _{2}}
g ( z ) = ( μ 1 μ 2 ) z + μ 2 {\displaystyle g\left(z\right)=\left(\mu _{1}-\mu _{2}\right)z+\mu _{2}}
 
Détermination des potentiels chimiques μ 1 {\displaystyle \mu _{1}} et μ 2 {\displaystyle \mu _{2}} et des activités a 1 {\displaystyle a_{1}} et a 1 {\displaystyle a_{1}} d'un mélange binaire par étude de l'évolution de l'enthalpie libre molaire G ¯ {\displaystyle {\bar {G}}} en fonction de la fraction molaire x 1 {\displaystyle x_{1}} à pression P {\displaystyle P} et température T {\displaystyle T} données.

L'équation de la tangente en x 1 {\displaystyle x_{1}} à la courbe G ¯ {\displaystyle {\bar {G}}} à P {\displaystyle P} et T {\displaystyle T} données est :

Équation de la tangente à la courbe G ¯ {\displaystyle {\bar {G}}}  : g ( z ) = ( μ 1 μ 2 ) z + μ 2 {\displaystyle g\left(z\right)=\left(\mu _{1}-\mu _{2}\right)z+\mu _{2}}

Cette tangente intercepte[2] :

  • l'axe x 1 = 0 {\displaystyle x_{1}=0} en μ 2 ( P , T , x 1 ) {\displaystyle \mu _{2}\!\left(P,T,x_{1}\right)} ,
  • l'axe x 1 = 1 {\displaystyle x_{1}=1} en μ 1 ( P , T , x 1 ) {\displaystyle \mu _{1}\!\left(P,T,x_{1}\right)} .

De plus, si l'on connaît μ 1 {\displaystyle \mu _{1}^{*}} et μ 2 {\displaystyle \mu _{2}^{*}} les potentiels chimiques respectifs des deux espèces pures à P {\displaystyle P} et T {\displaystyle T} , alors cette courbe permet également de déterminer a 1 {\displaystyle a_{1}} et a 2 {\displaystyle a_{2}} les activités chimiques respectives des deux espèces, puisque :

μ 1 = μ 1 + R T ln a 1 {\displaystyle \mu _{1}=\mu _{1}^{*}+RT\ln a_{1}}
μ 2 = μ 2 + R T ln a 2 {\displaystyle \mu _{2}=\mu _{2}^{*}+RT\ln a_{2}}

Les écarts :

G ¯ 1 mix = μ 1 mix = μ 1 μ 1 = R T ln a 1 {\displaystyle {\bar {G}}_{1}^{\text{mix}}=\mu _{1}^{\text{mix}}=\mu _{1}-\mu _{1}^{*}=RT\ln a_{1}}
G ¯ 2 mix = μ 2 mix = μ 2 μ 2 = R T ln a 2 {\displaystyle {\bar {G}}_{2}^{\text{mix}}=\mu _{2}^{\text{mix}}=\mu _{2}-\mu _{2}^{*}=RT\ln a_{2}}

sont les enthalpies libres molaires partielles de mélange, ou potentiels chimiques de mélange.

Rappelons également que dans le cas particulier du potentiel chimique :

  • pour x 1 = 0 {\displaystyle x_{1}=0} l'espèce 2 est pure :
    • μ 1 = {\displaystyle \mu _{1}=-\infty } et a 1 = 0 {\displaystyle a_{1}=0}  ;
    • μ 2 = μ 2 {\displaystyle \mu _{2}=\mu _{2}^{*}} et a 2 = 1 {\displaystyle a_{2}=1}  ;
  • pour x 1 = 1 {\displaystyle x_{1}=1} l'espèce 1 est pure :
    • μ 1 = μ 1 {\displaystyle \mu _{1}=\mu _{1}^{*}} et a 1 = 1 {\displaystyle a_{1}=1}  ;
    • μ 2 = {\displaystyle \mu _{2}=-\infty } et a 2 = 0 {\displaystyle a_{2}=0} .

Généralisation de la méthode

Étant donné que la relation de Gibbs-Duhem est généralisable à toute grandeur extensive X {\displaystyle X} , la méthode graphique employée ci-dessus pour déterminer les potentiels chimiques est déclinable pour la détermination des grandeurs molaires partielles X ¯ i {\displaystyle {\bar {X}}_{i}} d'un mélange binaire.

Il suffit donc de connaitre, à pression et température données, l'évolution de la grandeur molaire partielle X ¯ {\displaystyle {\bar {X}}} d'un mélange binaire sur toute l'étendue de composition entre le corps 1 pur et le corps 2 pur. La courbe de X ¯ {\displaystyle {\bar {X}}} est tracée en fonction de la composition, par ex. en fonction de x 1 {\displaystyle x_{1}} la fraction molaire du corps 1. Pour une valeur donnée de x 1 {\displaystyle x_{1}} , les grandeurs molaires partielles X ¯ 1 {\displaystyle {\bar {X}}_{1}} du corps 1 et X ¯ 2 {\displaystyle {\bar {X}}_{2}} du corps 2 sont déterminées par les intersections de la tangente à la courbe X ¯ {\displaystyle {\bar {X}}} en ce point avec les droites respectivement x 1 = 1 {\displaystyle x_{1}=1} et x 1 = 0 {\displaystyle x_{1}=0} . Cette méthode est appelée méthode de Roozeboom[a],[3].

Le tableau ci-dessous résume les équivalences entre grandeurs sur les exemples de l'enthalpie libre et du volume. Nous rappelons que toutes ces grandeurs sont déterminées aux mêmes pression et température.

Détermination des grandeurs molaires partielles
Grandeur extensive du mélange
X {\displaystyle X}
Enthalpie libre
G {\displaystyle G}
Volume
V {\displaystyle V}
Grandeur molaire du mélange
X ¯ {\displaystyle {\bar {X}}}
G ¯ {\displaystyle {\bar {G}}} V ¯ {\displaystyle {\bar {V}}}
Grandeur molaire partielle du corps i {\displaystyle i}
X ¯ i {\displaystyle {\bar {X}}_{i}}
G ¯ i = μ i {\displaystyle {\bar {G}}_{i}=\mu _{i}} V ¯ i {\displaystyle {\bar {V}}_{i}}
Grandeur molaire du corps i {\displaystyle i} pur
X ¯ i {\displaystyle {\bar {X}}_{i}^{*}}
G ¯ i = μ i {\displaystyle {\bar {G}}_{i}^{*}=\mu _{i}^{*}} V ¯ i {\displaystyle {\bar {V}}_{i}^{*}}
Grandeur molaire partielle de mélange du corps i {\displaystyle i}
X ¯ i mix = X ¯ i X ¯ i {\displaystyle {\bar {X}}_{i}^{\text{mix}}={\bar {X}}_{i}-{\bar {X}}_{i}^{*}}
G ¯ i mix = μ i mix {\displaystyle {\bar {G}}_{i}^{\text{mix}}=\mu _{i}^{\text{mix}}}
= G ¯ i G ¯ i {\displaystyle ={\bar {G}}_{i}-{\bar {G}}_{i}^{*}}
= R T ln a i {\displaystyle =RT\ln a_{i}}
V ¯ i mix {\displaystyle {\bar {V}}_{i}^{\text{mix}}}
= V ¯ i V ¯ i {\displaystyle ={\bar {V}}_{i}-{\bar {V}}_{i}^{*}}

Test de cohérence

Si l'on dispose, pour un système binaire à pression et température données constantes, de mesures expérimentales des coefficients d'activité γ 1 {\displaystyle \gamma _{1}} et γ 2 {\displaystyle \gamma _{2}} des deux composants, représentés respectivement par les fractions molaires x 1 {\displaystyle x_{1}} et x 2 {\displaystyle x_{2}} , sur une large plage de composition du mélange binaire (du corps 1 pur au corps 2 pur), la relation de Gibbs-Duhem permet d'établir un test de cohérence de ces mesures. On a les relations suivantes :

  • enthalpie libre molaire d'excès : G ¯ E R T = x 1 ln γ 1 + x 2 ln γ 2 {\displaystyle {{\bar {G}}^{\text{E}} \over RT}=x_{1}\ln \gamma _{1}+x_{2}\ln \gamma _{2}}
  • relation de Gibbs-Duhem à pression et température constantes : x 1 d ln γ 1 + x 2 d ln γ 2 = 0 {\displaystyle x_{1}\,\mathrm {d} \ln \gamma _{1}+x_{2}\,\mathrm {d} \ln \gamma _{2}=0}
  • contrainte sur les fractions molaires : x 1 + x 2 = 1 {\displaystyle x_{1}+x_{2}=1}

Rappelons que par définition du coefficient d'activité a i = x i γ i {\displaystyle a_{i}=x_{i}\gamma _{i}} , d'où :

  • pour x 1 = 0 {\displaystyle x_{1}=0} l'espèce 2 est pure :
    • a 1 = 0 {\displaystyle a_{1}=0} et γ 1 = γ 1 , 2 0 {\displaystyle \gamma _{1}=\gamma _{1,2}^{\infty }\neq 0} , coefficient d'activité de 1 infiniment diluée dans 2, de valeur non nulle ;
    • a 2 = 1 {\displaystyle a_{2}=1} et γ 2 = 1 {\displaystyle \gamma _{2}=1}  ;
  • pour x 1 = 1 {\displaystyle x_{1}=1} l'espèce 1 est pure :
    • a 1 = 1 {\displaystyle a_{1}=1} et γ 1 = 1 {\displaystyle \gamma _{1}=1}  ;
    • a 2 = 0 {\displaystyle a_{2}=0} et γ 2 = γ 2 , 1 0 {\displaystyle \gamma _{2}=\gamma _{2,1}^{\infty }\neq 0} , coefficient d'activité de 2 infiniment diluée dans 1, de valeur non nulle.

En différentiant l'expression de l'enthalpie libre d'excès :

d [ G ¯ E R T ] = ln γ 1 d x 1 + ln γ 2 d x 2 + x 1 d ln γ 1 + x 2 d ln γ 2 = 0 relation de Gibbs-Duhem {\displaystyle \mathrm {d} \left[{{\bar {G}}^{\text{E}} \over RT}\right]=\ln \gamma _{1}\,\mathrm {d} x_{1}+\ln \gamma _{2}\,\mathrm {d} x_{2}+\underbrace {x_{1}\,\mathrm {d} \ln \gamma _{1}+x_{2}\,\mathrm {d} \ln \gamma _{2}} _{=0 \atop {\text{relation de Gibbs-Duhem}}}}

En différentiant la relation entre fractions molaires :

d x 1 + d x 2 = 0 {\displaystyle \mathrm {d} x_{1}+\,\mathrm {d} x_{2}=0}

En retenant arbitrairement x 1 {\displaystyle x_{1}} comme variable, on obtient :

d [ G ¯ E R T ] = ln γ 1 d x 1 ln γ 2 d x 1 = ln ( γ 1 γ 2 ) d x 1 {\displaystyle \mathrm {d} \left[{{\bar {G}}^{\text{E}} \over RT}\right]=\ln \gamma _{1}\,\mathrm {d} x_{1}-\ln \gamma _{2}\,\mathrm {d} x_{1}=\ln \!\left({\gamma _{1} \over \gamma _{2}}\right)\,\mathrm {d} x_{1}}

Si l'on intègre cette relation sur toute l'étendue de la plage de variation de x 1 {\displaystyle x_{1}}  :

x 1 = 0 x 1 = 1 d [ G ¯ E R T ] = G ¯ E R T | x 1 = 1 G ¯ E R T | x 1 = 0 = 0 1 ln ( γ 1 γ 2 ) d x 1 {\displaystyle \int _{x_{1}=0}^{x_{1}=1}\,\mathrm {d} \left[{{\bar {G}}^{\text{E}} \over RT}\right]=\left.{{\bar {G}}^{\text{E}} \over RT}\right|_{x_{1}=1}-\left.{{\bar {G}}^{\text{E}} \over RT}\right|_{x_{1}=0}=\int _{0}^{1}\ln \!\left({\gamma _{1} \over \gamma _{2}}\right)\,\mathrm {d} x_{1}}

Or par définition pour x 1 = 1 {\displaystyle x_{1}=1} le corps 1 est pur et G ¯ E ( x 1 = 1 ) = 0 {\displaystyle {\bar {G}}^{\text{E}}\!\left(x_{1}=1\right)=0} . De même pour x 1 = 0 {\displaystyle x_{1}=0} le corps 2 est pur et G ¯ E ( x 1 = 0 ) = 0 {\displaystyle {\bar {G}}^{\text{E}}\!\left(x_{1}=0\right)=0} . On obtient un test de cohérence sur les coefficients d'activité du système binaire :

Test de cohérence : 0 1 ln ( γ 1 γ 2 ) d x 1 = 0 {\displaystyle \int _{0}^{1}\ln \!\left({\gamma _{1} \over \gamma _{2}}\right)\,\mathrm {d} x_{1}=0}

Dans une représentation graphique de ln ( γ 1 γ 2 ) {\displaystyle \ln \!\left({\gamma _{1} \over \gamma _{2}}\right)} en fonction de x 1 {\displaystyle x_{1}} , l'aire située entre l'axe des abscisses et les valeurs positives de la fonction doit être égale à l'aire entre l'axe des abscisses et les valeurs négatives de la fonction ; si ce n'est pas le cas, cela indique que les mesures sont biaisées.

Ce test ne garantit pas la justesse des coefficients d'activité mesurés : il indique seulement que toutes les mesures ont été effectuées de la même façon, mais pas que le protocole de mesure est correct.

Notes et références

Notes

  1. Hendrik Willem Bakhuis Roozeboom (1854-1907), chimiste et professeur d'université néerlandais, membre de l'Académie royale néerlandaise des arts et des sciences.

Références

  1. Réseau Numéliphy, Relations de Gibbs-Duhem.
  2. Jacques Schwartzentruber, École nationale supérieure des mines d'Albi-Carmaux, « Construction géométrique du potentiel chimique », (consulté le ).
  3. Bernard Claudel, Propriétés thermodynamiques des fluides, vol. B 8020, Techniques de l'ingénieur, , 46 p. (lire en ligne), p. 6.

Bibliographie

  • Jean-Pierre Corriou, Thermodynamique chimique : Définitions et relations fondamentales, vol. J 1025, Techniques de l'ingénieur, coll. « base documentaire : Thermodynamique et cinétique chimique, pack Opérations unitaires. Génie de la réaction chimique, univers Procédés chimie - bio - agro », (lire en ligne), p. 1-19.
  • Jean-Pierre Corriou, Thermodynamique chimique : Diagrammes thermodynamiques, vol. J 1026, coll. « base documentaire : Thermodynamique et cinétique chimique, pack : Opérations unitaires. Génie de la réaction chimique, univers : Procédés chimie - bio - agro », (lire en ligne), p. 1-30.
  • Jean Vidal, Thermodynamique : application au génie chimique et à l'industrie pétrolière, Paris, Éditions Technip, coll. « Publications de l'Institut français du pétrole. », , 500 p. (ISBN 978-2-710-80715-5, OCLC 300489419, lire en ligne), p. 153-154.

Liens externes

Articles connexes

  • icône décorative Portail de la physique
  • icône décorative Portail de la chimie