# Équilibre thermodynamique

En thermodynamique, un équilibre thermodynamique est l'état d'un système ne subissant aucune évolution à l'échelle macroscopique. Les grandeurs intensives caractérisant ce système (notamment la pression, la température et les potentiels chimiques) sont alors homogènes dans l'espace et constantes dans le temps. Cet équilibre est stable s'il n'est pas susceptible d'évoluer sur une échelle de temps infinie, ou sous le coup d'une perturbation extérieure ou d'une fluctuation locale interne.
Tout système évolue selon les premier et deuxième principes de la thermodynamique, qui définissent les fonctions d'état énergie interne et entropie. Un équilibre thermodynamique est stable lorsque le système ne peut plus évoluer, car cela impliquerait une diminution de l'entropie contraire au deuxième principe. L'équilibre thermodynamique d'un système isolé (qui n'échange ni travail, ni chaleur, ni matière avec l'extérieur) est caractérisé par un maximum de l'entropie et un minimum de l'énergie interne. L'équilibre d'un système fermé (qui n'échange pas de matière avec l'extérieur) sous contrainte est caractérisé par le minimum d'un potentiel thermodynamique spécifique (qui est également une fonction d'état) correspondant à un maximum de l'entropie globale du système et de son extérieur.
L'entropie et les potentiels thermodynamiques possèdent en conséquence des caractéristiques de concavité ou de convexité qui déterminent le signe de certains coefficients calorimétriques et thermoélastiques, ainsi que certaines propriétés des potentiels chimiques des espèces constituant le système. D'autre part, deux systèmes en équilibre l'un avec l'autre ont atteint simultanément :
* un équilibre mécanique lorsque les deux systèmes ont la même pression ;
* un équilibre thermique lorsque les deux systèmes ont la même température ;
* un équilibre diffusif lorsque chaque espèce chimique a le même potentiel chimique dans les deux systèmes.
Les équilibres de phases et les équilibres chimiques sont des exemples d'équilibres thermodynamiques. Ces principes expliquent également la coalescence des milieux dispersés, ceux-ci tendant à minimiser leur surface.

## Définitions et principes généraux
### Éléments de définition d'un équilibre thermodynamique
#### Définition de l'équilibre thermodynamique

Les états d'équilibre sont les objets d'étude de la thermodynamique. Un processus thermodynamique est une transformation d'un système passant d'un état d'équilibre à un autre. Une transformation réversible est effectuée par l'intermédiaire d'équilibres successifs, tandis qu'une transformation irréversible est effectuée via des états hors équilibre. Dans les faits, la définition d'un état d'équilibre est bien souvent circulaire : un système à l'équilibre thermodynamique est un système qui peut être décrit par les lois de la thermodynamique, tandis qu'un échec de la thermodynamique peut être interprété comme une preuve que le système n'est pas à l'équilibre.
L'équilibre thermodynamique peut être défini ainsi$^,$$^,$ :
* Équilibre thermodynamique État que tout système isolé atteint après un temps suffisamment long et caractérisé macroscopiquement par un nombre fini de variables d'état.
Les variables, ou grandeurs, décrivant le système sont extensives : en général on considère le volume $V$, l'énergie interne $U$ et les quantités $n_{1}$, $n_{2}$, etc. des espèces chimiques constituant le système. Les variables décrivant l'équilibre sont intensives : en général on considère la pression $P$, la température $T$ et les potentiels chimiques $\mu _{1}$, $\mu _{2}$, etc. des espèces$^,$. Ces grandeurs intensives ne sont définies et mesurables qu'à l'équilibre.
On ajoute à cette définition la nécessité pour le système d'être homogène et isotrope, c'est-à-dire le fait que les grandeurs intensives sont identiques en tout point du système et ne dépendent pas de la direction selon laquelle le système est observé. Dans les faits, un système peut être à l'équilibre sans être homogène, comme dans le cas d'une colonne de liquide dans laquelle la pression dépend de l'emplacement dans la colonne. Le critère de l'homogénéité permet de simplifier la description thermodynamique du système et de son équilibre. D'autre part, les grandeurs décrivant l'équilibre sont invariantes dans le temps : le système n'évolue plus et ne peut plus évoluer.
Cette définition est issue d'observations courantes qui montrent qu'un système isolé a tendance à évoluer vers un état indépendant de son histoire. Dans la pratique, cette définition pose certaines difficultés, notamment quant à l'établissement d'un système isolé et quant à la durée d'observation nécessaire à l'obtention de l'équilibre.

#### Systèmes thermodynamiques

* Dans la colonne d'un baromètre à mercure, la pression est quasi nulle au sommet et égale à la pression atmosphérique au pied, par effet de la gravité.
* Un ménisque est un effet de surface.
* Un piston est un système fermé. Son enceinte ne laisse ni entrer ni sortir la matière. Par contre elle est déformable et laisse passer la chaleur.
* Un tube en aluminium (5 mm d'épaisseur) ayant subi une pression extérieure de 700 bar.
* Le vase de Dewar est un récipient aux parois adiabatiques.
Un système thermodynamique, ou système, est une portion de l'Univers délimitée par la pensée du reste de l'Univers, appelé système extérieur ou extérieur$^,$. La paroi, aussi appelée enceinte ou enveloppe, qui marque la frontière entre le système et l'extérieur, qui peut être virtuelle, détermine la nature du système selon les échanges qu'elle rend possibles au cours d'une transformation du système. Un système ouvert peut échanger de l'énergie (sous forme de travail et de chaleur) et de la matière avec l'extérieur. Un système fermé peut échanger de l'énergie mais pas de matière. Un système isolé n'échange ni énergie ni matière.
Ces divers types de système sont supposés être homogènes, ce qui signifie que les grandeurs intensives (pression, température, potentiels chimiques) qui décrivent un système sont identiques en tout point de ce système. Ceci implique que le système ne baigne pas dans un champ de force extérieur, notamment un champ électromagnétique ou gravitationnel. La gravité induit que la pression n'est pas homogène dans un système réel : dans une colonne de liquide, en raison des lois de l'hydrostatique, la pression est plus forte au fond de la colonne qu'en surface. Par exemple, 10 m de colonne d'eau induisent une différence de pression d'environ 1 bar, et la loi de Henry induit que les concentrations en gaz dissouts dans une telle colonne sont deux fois plus importantes au fond de la colonne qu'en surface. La gravité induit des gradients de pression et de concentration dans une colonne de liquide, qui n'est donc pas homogène. La gravité est un facteur moins influent pour les gaz, qui ont des masses volumiques plus faibles que les liquides. Le volume du système étudié doit par conséquent être suffisamment petit pour pouvoir, expérimentalement, négliger l'influence de champs de forces extérieurs, et notamment de la gravité.
Le volume doit toutefois rester suffisamment important pour éviter les effets de surface. En effet, la matière peut avoir un comportement altéré aux interfaces entre phases : entre un gaz et un liquide, entre les parois d'un récipient et son contenu. Des gradients de concentration peuvent, par exemple, s'établir à ces interfaces, des forces particulières (tension superficielle, etc.) peuvent y agir. Ces phénomènes sont cependant très localisés et négligeables sur des volumes suffisamment importants.
Un système fermé n'échange aucune matière avec l'extérieur. Sa paroi est imperméable, aucune matière ne peut la traverser. Par contre, cette paroi peut être déformable, le volume du système peut ainsi évoluer au cours de la transformation par le travail des forces de pression. De même, cette paroi peut laisser passer de la chaleur par conduction ou rayonnement. Un piston est un exemple typique de système fermé. Expérimentalement, l'obtention d'une étanchéité parfaite est le point le plus délicat d'un tel système. De plus, la paroi elle-même ne doit pas subir de dégradation au cours du temps, sous l'effet des déformations ou de la chaleur. Son matériau doit être inerte, il ne doit pas relarguer de composants vers le système, ou absorber des composants du système, ou réagir avec le contenu.
L'absence d'échange de travail avec l'extérieur suppose, outre l'absence d'un champ de force extérieur, que les forces de pression ne peuvent pas déformer le système, et donc que son volume reste constant au cours de la transformation. La paroi du système doit pouvoir résister, le cas échéant, à un vide poussé ou à de fortes pressions, tant à l'intérieur qu'à l'extérieur du système.
L'absence d'échange de chaleur du système avec son environnement est approximativement réalisable à l'aide d'isolants thermiques. Il n'existe pas d'isolant parfait, tout matériau laissant passer, à des degrés divers, de la chaleur par conduction ou rayonnement. La capacité d'un matériau à s'opposer à un flux thermique est appelée résistance thermique ; une transformation effectuée en l'absence d'échange de chaleur avec l'extérieur est dite « adiabatique » et suppose que les parois du système isolé ont une résistance thermique infinie. De tels matériaux n'existent pas, et les déperditions thermiques du système ne peuvent être négligées que sur des temps d'observation réduits.
Dans la pratique, un système isolé ne peut donc être réalisé que sur de petits volumes, afin de négliger l'influence de la gravité, mais suffisamment importants pour rendre les effets de surface négligeables. Sa paroi doit être totalement imperméable et inerte vis-à-vis du contenu et de l'extérieur du système, et résister à la différence de pression entre contenu et extérieur. Elle doit de plus empêcher toute déperdition thermique. Dans les faits, un tel système est idéal et n'existe pas, mais on peut obtenir des systèmes s'en approchant sur des temps relativement courts.

#### Échelle macroscopique, équilibre stationnaire

* À l'échelle microscopique les molécules se déforment en permanence.
* À l'échelle microscopique les molécules se déplacent en permanence (mouvement brownien).
* À l'échelle macroscopique aucun mouvement n'est observé.
En thermodynamique, la matière étudiée est le plus souvent constituée de particules telles que des molécules ou des ions. Ces structures composées d'atomes se déforment en permanence selon leurs degrés de liberté internes. Dans un fluide (liquide, gaz) ces espèces chimiques se déplacent en permanence (mouvement brownien) dans l'espace. Dans un solide les atomes vibrent en permanence autour d'une position moyenne. De même, dans les équilibres chimiques des espèces se transforment en permanence, des réactifs donnant des produits et réciproquement. Dans les équilibres de phases, des particules passent en permanence d'une phase à l'autre au niveau de l'interface entre phases. À l'échelle microscopique, les composants d'un système thermodynamique quelconque ne sont donc jamais figés dans l'espace et le temps. À cette échelle, un système thermodynamique subit en permanence des fluctuations spatiales et temporelles des grandeurs qui le décrivent$^,$.
À l'échelle macroscopique ces fluctuations et mouvements ne sont pas perceptibles directement, mais sont reflétés par des grandeurs intensives telles que la pression, la température, les potentiels chimiques des espèces chimiques, etc. Lorsque ces grandeurs sont homogènes dans l'espace et constantes dans le temps, l'état du système est décrit comme étant celui d'un équilibre$^,$. Cet équilibre est qualifié d'équilibre stationnaire ou d'équilibre dynamique, car à l'échelle microscopique les échanges entre particules et les mouvements de celles-ci ont toujours lieu (au contraire d'un équilibre statique).
Ce sont donc les moyennes des grandeurs intensives mesurées qui décrivent son équilibre : localement, ces grandeurs fluctuent dans le temps et l'espace. Ainsi, la thermodynamique décrit les états macroscopiques de la matière et énonce les lois qui régissent les échanges d'énergie entre de grands ensembles de particules. Elle ne dit rien sur le détail des interactions et mouvements à l'échelle microscopique, dont l'étude relève d'autres domaines de la physique.

#### Équilibre stable

* Exemple d'équilibre instable. Le diamant (à gauche) est une forme instable du carbone. Sur une très longue échelle de temps, le diamant se transforme en graphite (à droite).
* Exemple d'équilibre métastable. La perturbation appliquée à cette eau liquide à température négative, en état de surfusion, déclenche son passage à l'état solide.
Certaines transformations extrêmement lentes conduisent à l'observation d'états hors équilibre, dits équilibres instables, qui peuvent paraitre stables sur une échelle de temps relativement longue. Ainsi le diamant est-il une forme instable du carbone : il tend naturellement à se transformer en graphite, qui est la forme stable dans les conditions normales de température et de pression, mais cette réaction n'est observable qu'à des échelles de temps géologiques. De même, des phénomènes de diffusion de la matière peuvent être extrêmement lents et non observables sur des temps courts, l'équilibre apparent observé est également instable.
Certains états observés sont des états hors équilibre qui perdurent dans le temps jusqu'à ce qu'une perturbation déclenche une transformation souvent violente : ces états sont dits équilibres métastables. Par exemple, un mélange d'hydrogène et d'oxygène peut rester stable dans le temps, mais une simple étincelle déclenche une réaction violente qui produit de l'eau. De même, l'eau peut être maintenue liquide en dessous de 0 °C à pression atmosphérique : cet état est appelé état de surfusion, il suffit alors de l'introduction d'une poussière ou d'une fluctuation de température pour que l'eau se solidifie brutalement. Des solutions liquides peuvent contenir un soluté (gaz ou solide) dissout en état de sursaturation, c'est-à-dire à des concentrations plus importantes que celles décrites par la thermodynamique, jusqu'à ce qu'une perturbation déclenche un phénomène brutal de désaturation (respectivement un dégazage et une précipitation) permettant d'établir l'état d'équilibre stable de la solution, l'état de saturation prédit par la thermodynamique.
Un équilibre stable peut donc se définir comme étant un état non susceptible d'évoluer sur une échelle de temps infinie, ni sous le coup d'une perturbation extérieure ou d'une fluctuation locale interne$^,$$^,$.

#### Équilibre thermodynamique local

Un système thermodynamique peut parfois être décomposé, à l'échelle mésoscopique (intermédiaire entre microscopique et macroscopique$^,$$^,$), en sous-systèmes décrits individuellement par des grandeurs intensives mesurables (pression, température, potentiels chimiques, etc.) dont les valeurs diffèrent d'un sous-système à l'autre. Le système se trouve donc dans un état hors équilibre, et chacun des sous-systèmes est décrit comme étant dans un état d'équilibre thermodynamique local. Les inhomogénéités des grandeurs entre sous-systèmes induisent entre eux des phénomènes de transport de matière et d'énergie (voir la section Équilibre d'un système hétérogène). Les équilibres locaux sont par conséquent instables et peuvent évoluer ; l'équilibre stable du système correspond à l'homogénéité des grandeurs sur l'ensemble de ses sous-systèmes$^,$$^,$$^,$$^,$.
Par exemple, une barre de fer maintenue entre deux températures différentes à ses extrémités présente un gradient de température sur sa longueur, induisant un phénomène de conduction thermique. De même, à l'interface d'un glaçon et de l'eau liquide dans un verre, les deux phases sont à la même température (0 °C), mais le glaçon peut être plus froid en son cœur et l'eau liquide plus chaude loin du glaçon ; l'interface solide-liquide est en état d'équilibre thermodynamique local, instable puisque le glaçon fond. L'atmosphère est un exemple de système composé d'équilibres thermodynamiques locaux, puisque l'on peut y définir localement des pressions et des températures qui sont globalement inhomogènes sur l'ensemble de la planète.
Il existe des cas dans lesquels des systèmes hors équilibre ne peuvent être décomposés en sous-systèmes en équilibre local, certaines grandeurs, par exemple la température, ne pouvant être définies localement.

### Principes de la thermodynamique
#### Premier principe de la thermodynamique

Le premier principe de la thermodynamique est un principe de conservation de l'énergie interne, notée $U$.
* Premier principe de la thermodynamique Au cours d'une transformation quelconque, la variation de l'énergie interne $U$ d'un système fermé est égale à la somme du travail $W$ et de la chaleur $Q$ échangés par le système avec l'extérieur : $\mathrm {d} U=\delta W+\delta Q$
En conséquence, l'énergie interne d'un système isolé ne varie pas au cours de la transformation, soit :
$\mathrm {d} U=0$
Si un système fermé, d'énergie interne $U$, et son extérieur, d'énergie interne $U_{\text{ext}}$, forment ensemble un système isolé, alors l'énergie globale $U+U_{\text{ext}}$ du système isolé est invariante : $\mathrm {d} \left[U+U_{\text{ext}}\right]=\mathrm {d} U+\mathrm {d} U_{\text{ext}}=0$. Ainsi, l'une des deux énergies $U$ et $U_{\text{ext}}$ peut varier du moment que l'autre énergie compense cette variation par une variation opposée.

#### Deuxième principe de la thermodynamique

Le deuxième principe de la thermodynamique est un principe d'évolution de la transformation ayant lieu dans le système. Ce principe introduit la notion d'entropie, notée $S$.
Pour un système fermé échangeant de la chaleur $Q$ avec un autre système, l'entropie est définie par l'inégalité de Clausius :
* Entropie, inégalité de Clausius : $\mathrm {d} S\geq {\delta Q \over T}$
avec $T$ la température absolue. La transformation est :
* réversible si $\mathrm {d} S={\delta Q \over T}$ ;
* irréversible si $\mathrm {d} S>{\delta Q \over T}$.
Dans un système isolé, n'échangeant pas de chaleur avec l'extérieur, soit $\delta Q=0$, un processus réversible se traduit donc par $\mathrm {d} S=0$ et un processus irréversible par $\mathrm {d} S>0$. Le deuxième principe de la thermodynamique énonce que :
* Deuxième principe de la thermodynamique Au cours d'une transformation quelconque, l'entropie $S$ d'un système isolé ne peut qu'augmenter : $\mathrm {d} S\geq 0$
Soit un système fermé ou isolé, siège d'une transformation quelconque impliquant un échange de chaleur $\delta Q$ à température $T$. On note$^,$$^,$ :
$\mathrm {d} S=\mathrm {d} S_{\text{e}}+\mathrm {d} S_{\text{i}}$
avec :
* $\mathrm {d} S$ la variation d'entropie du système due à la transformation ;
* $\mathrm {d} S_{\text{e}}=\delta Q/T$ la variation d'entropie due exclusivement à l'échange de chaleur avec l'extérieur ; $S_{\text{e}}$ est l'entropie d'échange ;
* $\mathrm {d} S_{\text{i}}$ la génération d'entropie (aussi appelée création ou production d'entropie), liée à l'irréversibilité de la transformation ; $S_{\text{i}}$ est l'entropie d'irréversibilité.
Le terme $\mathrm {d} S_{\text{e}}$ peut être positif (le système gagne de la chaleur), nul (processus adiabatique) ou négatif (le système perd de la chaleur). Le terme $\mathrm {d} S_{\text{i}}$ ne peut être que positif (dans une transformation irréversible) ou nul (dans une transformation réversible). On note également :
$T\,\mathrm {d} S=\delta Q+\delta \Pi$
Le terme $\delta \Pi =T\,\mathrm {d} S_{\text{i}}$ est appelé travail non compensé ou chaleur non compensée. Il ne peut être que positif ou nul$^,$.
Si $\delta Q$ est négatif, l'entropie $S$ d'un système fermé peut diminuer, soit $\mathrm {d} S\leq 0$. On note $S_{\text{ext}}$ l'entropie de l'extérieur de ce système. En vertu du deuxième principe, le système et son extérieur formant un système isolé, l'entropie globale $S+S_{\text{ext}}$ ne peut qu'augmenter, soit $\mathrm {d} \left[S+S_{\text{ext}}\right]=\mathrm {d} S+\mathrm {d} S_{\text{ext}}\geq 0$. L'entropie du système peut donc diminuer si la variation d'entropie de l'extérieur compense cette diminution de telle sorte que $\mathrm {d} S_{\text{ext}}\geq -\mathrm {d} S\geq 0$. Si l'entropie $S$ du système fermé diminue, l'entropie $S_{\text{ext}}$ de l'extérieur ne peut qu'augmenter de façon que l'entropie globale $S+S_{\text{ext}}$ du système isolé augmente. Inversement, l'entropie de l'extérieur peut également diminuer si l'entropie du système fermé augmente et si l'entropie globale du système isolé augmente.

#### Potentiels thermodynamiques

Considérons un système fermé, appelé « système » par la suite, subissant une transformation quelconque. Ce système et son extérieur forment ensemble un système isolé qui possède son propre extérieur. On pose :
* $P$ la pression, en considérant que le système et son extérieur sont constamment à la même pression, celle-ci pouvant évoluer ;
* $T$ la température, en considérant que le système et son extérieur sont constamment à la même température, celle-ci pouvant évoluer ;
* $V$ le volume du système ;
* $V_{\text{ext}}$ le volume de l'extérieur du système.
Par définition, le système isolé n'échange pas de travail avec son propre extérieur. Il n'y a donc aucune force de pression pour déformer son volume $V+V_{\text{ext}}$ : le travail $\delta W_{P}=-P\,\mathrm {d} \left[V+V_{\text{ext}}\right]=0$. Le volume $V$ du système étudié peut donc évoluer au cours de la transformation en respectant la contrainte :
    $\mathrm {d} V+\mathrm {d} V_{\text{ext}}=0$
Pour le système, on considère :
* qu'il produit un travail des forces de pression égal à $\delta W=-P\,\mathrm {d} V$ ;
* qu'il est le siège d'une transformation potentiellement irréversible, soit $T\,\mathrm {d} S\geq \delta Q$ ; il peut s'agir, entre autres :
** du travail d'une force utile non récupéré, qui est alors dégradé en chaleur ;
** de travaux de forces non récupérables, telles que les frottements ou la viscosité, qui sont également dégradés en chaleur ;
** d'une réaction chimique ou d'un phénomène diffusif.
On pose : $\delta Q=T\,\mathrm {d} S+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$. Le terme $\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$ représente la variation d'énergie due à la transformation du système fermé. Ce terme doit être négatif selon le deuxième principe, il s'agit de l'opposé de la chaleur non compensée : $\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-\delta \Pi =-T\,\mathrm {d} S_{\text{i}}\leq 0$. La variation de l'énergie interne du système $\mathrm {d} U=\delta W+\delta Q$ vaut en conséquence :
$\mathrm {d} U=-P\,\mathrm {d} V+T\,\mathrm {d} S+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$
Pour l'extérieur, on considère :
* qu'il produit un travail des forces de pression égal à $\delta W_{\text{ext}}=-P\,\mathrm {d} V_{\text{ext}}=P\,\mathrm {d} V=-\delta W$ ;
* qu'il reçoit la chaleur émise par le système et que la transformation est réversible : $\delta Q_{\text{ext}}=T\,\mathrm {d} S_{\text{ext}}=-\delta Q$.
La variation de l'énergie interne de l'extérieur $\mathrm {d} U_{\text{ext}}=\delta W_{\text{ext}}+\delta Q_{\text{ext}}$ vaut en conséquence :
$\mathrm {d} U_{\text{ext}}=P\,\mathrm {d} V+T\,\mathrm {d} S_{\text{ext}}$
En application du premier principe, puisque le système et l'extérieur forment un système isolé, l'énergie globale $U+U_{\text{ext}}$ est conservée :
$\mathrm {d} U+\mathrm {d} U_{\text{ext}}=\mathrm {d} \left[U+U_{\text{ext}}\right]=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}+T\,\left[\mathrm {d} S+\mathrm {d} S_{\text{ext}}\right]=0$
En application du deuxième principe, puisque le système et l'extérieur forment un système isolé, l'entropie globale $S+S_{\text{ext}}$ ne peut que croître :
$\mathrm {d} S+\mathrm {d} S_{\text{ext}}=\mathrm {d} \left[S+S_{\text{ext}}\right]\geq 0$
À partir des relations (1), (2) et (3), on appelle potentiel thermodynamique une fonction d'état, notée de façon générique $\Phi$, telle que :
* Potentiel thermodynamique $\mathrm {d} \Phi =\mathrm {d} U+P\,\mathrm {d} V-T\,\mathrm {d} S=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]\leq 0$
On vérifie l'inégalité de Clausius : $\mathrm {d} S-{\delta Q \over T}=-{1 \over T}\left[\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}\right]\geq 0$.
Les divers potentiels thermodynamiques sont créés en fixant certaines conditions de la transformation :
* énergie interne $U$ :
    à volume $V$ et entropie $S$ constants : $\mathrm {d} \Phi =\mathrm {d} U\leq 0$ ;
* énergie libre, ou fonction de Helmholtz, $F=U-TS$ :
    à volume $V$ et température $T$ constants : $\mathrm {d} \Phi =\mathrm {d} U-T\,\mathrm {d} S=\mathrm {d} \left[U-TS\right]=\mathrm {d} F\leq 0$ ;
* enthalpie $H=U+PV$ :
    à pression $P$ et entropie $S$ constantes : $\mathrm {d} \Phi =\mathrm {d} U+P\,\mathrm {d} V=\mathrm {d} \left[U+PV\right]=\mathrm {d} H\leq 0$ ;
* enthalpie libre, ou fonction de Gibbs, $G=U+PV-TS$ :
    à pression $P$ et température $T$ constantes : $\mathrm {d} \Phi =\mathrm {d} U+P\,\mathrm {d} V-T\,\mathrm {d} S=\mathrm {d} \left[U+PV-TS\right]=\mathrm {d} G\leq 0$.
Au cours de la transformation, l'une des deux entropies $S$ et $S_{\text{ext}}$ peut décroître, du moment que l'entropie globale $S+S_{\text{ext}}$ du système isolé croît. Les potentiels thermodynamiques $\Phi$ ne peuvent que décroître.

### Convexité des fonctions thermodynamiques
#### Concavité de l'entropie

Les variables de l'entropie sont toutes extensives : volume $V$, énergie interne $U$, quantité de matière $n$, etc. On les note de façon générique $V,U,X_{1},X_{2},\cdots$ À l'équilibre, l'entropie est elle-même une fonction extensive : si l'on met en contact deux systèmes $1$ et $2$ à l'équilibre, alors l'entropie globale du système résultant $1+2$ vaut $S^{1+2}=S^{1}+S^{2}$. À contrario, si les deux systèmes fusionnés ne sont pas à l'équilibre, il y a, en vertu du deuxième principe, création d'entropie au cours du processus de fusion, et, lorsque le nouveau système a atteint l'équilibre, on a$^,$ : $S^{1+2}>S^{1}+S^{2}$.
Considérons deux systèmes thermodynamiques quelconques isolés, notés $\alpha$ et $\beta$, par exemple deux phases différentes, deux réservoirs de réactifs chimiques, deux solides à températures différentes, etc. Chacun de ces deux systèmes est à l'équilibre thermodynamique indépendamment de l'autre. Ils ont respectivement pour entropie$^,$ :
    $S^{\alpha }=S\!\left(V^{\alpha },U^{\alpha },X_{1}^{\alpha },X_{2}^{\alpha },\cdots \right)$
    $S^{\beta }=S\!\left(V^{\beta },U^{\beta },X_{1}^{\beta },X_{2}^{\beta },\cdots \right)$
On prélève une fraction quelconque $\tau$ (comprise entre 0 et 1) du système $\alpha$ et une fraction $1-\tau$ du système $\beta$. Les variables de l'entropie étant extensives, ces fractions de systèmes $\tau \alpha$ et $\left[1-\tau \right]\beta$ sont définies respectivement par des fractions des variables initiales :
    $V^{\tau \alpha }=\tau \,V^{\alpha },\,U^{\tau \alpha }=\tau \,U^{\alpha },\,X_{1}^{\tau \alpha }=\tau \,X_{1}^{\alpha },\,X_{2}^{\tau \alpha }=\tau \,X_{2}^{\alpha },\cdots$
    $V^{\left[1-\tau \right]\beta }=\left[1-\tau \right]\,V^{\beta },\,U^{\left[1-\tau \right]\beta }=\left[1-\tau \right]\,U^{\beta },\,X_{1}^{\left[1-\tau \right]\beta }=\left[1-\tau \right]\,X_{1}^{\beta },\,X_{2}^{\left[1-\tau \right]\beta }=\left[1-\tau \right]\,X_{2}^{\beta },\cdots$
L'entropie étant elle-même extensive à l'équilibre, les entropies de ces deux fractions de systèmes sont des fractions des entropies des systèmes $\alpha$ et $\beta$ initiaux et valent respectivement$^,$ :
    $S^{\tau \alpha }=S\!\left(\tau \,V^{\alpha },\tau \,U^{\alpha },\tau \,X_{1}^{\alpha },\tau \,X_{2}^{\alpha },\cdots \right)=\tau \,S\!\left(V^{\alpha },U^{\alpha },X_{1}^{\alpha },X_{2}^{\alpha },\cdots \right)=\tau \,S^{\alpha }$
    $S^{\left[1-\tau \right]\beta }=S\!\left(\left[1-\tau \right]\,V^{\beta },\left[1-\tau \right]\,U^{\beta },\left[1-\tau \right]\,X_{1}^{\beta },\left[1-\tau \right]\,X_{2}^{\beta },\cdots \right)=\left[1-\tau \right]\,S\!\left(V^{\beta },U^{\beta },X_{1}^{\beta },X_{2}^{\beta },\cdots \right)=\left[1-\tau \right]\,S^{\beta }$
Ces deux fractions sont mises en contact de façon à former un nouveau système isolé $\tau \alpha +\left[1-\tau \right]\beta$. Puisque les variables de l'entropie sont toutes extensives, elles s'additionnent. En particulier, en vertu du premier principe, le nouveau système étant isolé, son énergie interne ne peut pas évoluer et ne peut qu'être égale à $\tau \,U^{\alpha }+\left[1-\tau \right]\,U^{\beta }$. Le nouveau système est par conséquent défini par les variables :
    $V^{\tau \alpha +\left[1-\tau \right]\beta }=\tau \,V^{\alpha }+\left[1-\tau \right]\,V^{\beta },\,U^{\tau \alpha +\left[1-\tau \right]\beta }=\tau \,U^{\alpha }+\left[1-\tau \right]\,U^{\beta },\,X_{1}^{\tau \alpha +\left[1-\tau \right]\beta }=\tau \,X_{1}^{\alpha }+\left[1-\tau \right]\,X_{1}^{\beta },\,X_{2}^{\tau \alpha +\left[1-\tau \right]\beta }=\tau \,X_{2}^{\alpha }+\left[1-\tau \right]\,X_{2}^{\beta },\cdots$
et son entropie est une fonction telle que$^,$ :
    $S^{\tau \alpha +\left[1-\tau \right]\beta }=S\!\left(\tau \,V^{\alpha }+\left[1-\tau \right]\,V^{\beta },\tau \,U^{\alpha }+\left[1-\tau \right]\,U^{\beta },\tau \,X_{1}^{\alpha }+\left[1-\tau \right]\,X_{1}^{\beta },\tau \,X_{2}^{\alpha }+\left[1-\tau \right]\,X_{2}^{\beta },\cdots \right)$
Si les deux fractions sont déjà à l'équilibre thermodynamique entre elles au moment de leur mise en contact, alors le nouveau système n'évolue pas et a pour entropie $\tau \,S^{\alpha }+\left[1-\tau \right]\,S^{\beta }$. Si les deux fractions ne sont pas immédiatement à l'équilibre entre elles, alors le nouveau système est le siège d'une transformation (par exemple un transfert de chaleur, une diffusion de matière d'un sous-système à l'autre, une réaction chimique…). En vertu du deuxième principe, l'entropie du nouveau système ne peut qu'augmenter au cours de cette transformation et, lorsque celui-ci a atteint l'équilibre, ne peut qu'être supérieure à la somme des entropies des deux fractions. À l'équilibre de ce nouveau système on a par conséquent$^,$ :
$S^{\tau \alpha +\left[1-\tau \right]\beta }\geq \tau \,S^{\alpha }+\left[1-\tau \right]\,S^{\beta }$
Les systèmes $\alpha$ et $\beta$ sont quelconques, la fraction $\tau$ est également quelconque, l'entropie est donc une fonction concave par rapport à ses variables naturelles. En termes mathématiques, la matrice hessienne de l'entropie est définie négative, d'où$^,$$^,$$^,$ :
* L'entropie est une fonction concave : $\mathrm {d} ^{2}S\leq 0$

#### Convexité de l'énergie interne

On considère un système fermé siège d'une transformation quelconque, que l'on traduit par le terme $\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$ dans la différentielle du premier ordre de $U$ :
    $\mathrm {d} U=-P\,\mathrm {d} V+T\,\mathrm {d} S+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$
Les variables $X_{i}$ sont les mêmes pour l'énergie interne et l'entropie. Les variables de l'énergie interne (volume, entropie, quantité de matière, etc.) sont par conséquent toutes extensives. On a la différentielle du deuxième ordre de l'énergie interne :
    $\mathrm {d} ^{2}U=-\mathrm {d} P\,\mathrm {d} V+\mathrm {d} T\,\mathrm {d} S+\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}$
Pour l'entropie on a, par réagencement de $\mathrm {d} U$ :
    $\mathrm {d} S={P \over T}\,\mathrm {d} V+{1 \over T}\,\mathrm {d} U-\sum _{i}{U_{X_{i}} \over T}\,\mathrm {d} X_{i}$
    ${\begin{aligned}\mathrm {d} ^{2}S&=\mathrm {d} {P \over T}\,\mathrm {d} V+\mathrm {d} {1 \over T}\,\mathrm {d} U-\sum _{i}\mathrm {d} {U_{X_{i}} \over T}\,\mathrm {d} X_{i}\\&={1 \over T}\,\mathrm {d} P\,\mathrm {d} V-{P \over T^{2}}\,\mathrm {d} T\,\mathrm {d} V-{1 \over T^{2}}\,\mathrm {d} T\,\mathrm {d} U+\sum _{i}{U_{X_{i}} \over T^{2}}\,\mathrm {d} T\,\mathrm {d} X_{i}-\sum _{i}{1 \over T}\,\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\end{aligned}}$
Par substitution de $\mathrm {d} U$ on obtient :
    $\mathrm {d} ^{2}S={1 \over T}\,\mathrm {d} P\,\mathrm {d} V-{1 \over T}\,\mathrm {d} T\,\mathrm {d} S-\sum _{i}{1 \over T}\,\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}$
Par conséquent :
$\mathrm {d} ^{2}U=-T\,\mathrm {d} ^{2}S$
Puisque la température thermodynamique $T$ ne peut être que positive, alors $\mathrm {d} ^{2}S$ et $\mathrm {d} ^{2}U$ sont de signes opposés. L'entropie étant une fonction concave, l'énergie interne est une fonction convexe par rapport à ses variables naturelles. En termes mathématiques, la matrice hessienne de l'énergie interne est définie positive, d'où :
* L'énergie interne est une fonction convexe : $\mathrm {d} ^{2}U\geq 0$

#### Convexité des autres potentiels thermodynamiques

Les potentiels thermodynamiques sont des transformées de Legendre de l'énergie interne $U$. Cette transformation remplace une variable extensive de $U$ par une variable intensive du potentiel thermodynamique résultant, tout en inversant la convexité de la fonction par rapport à cette variable. En conséquence, les potentiels thermodynamiques sont concaves par rapport à leur(s) variable(s) intensive(s) et convexes par rapport à leurs variables extensives (voir la section Notations et outils mathématiques)$^,$ :
* l'énergie interne $U$ est convexe par rapport à toutes ses variables, qui sont toutes extensives : volume $V$, entropie $S$, quantité de matière $n$, etc. ;
* l'énergie libre $F$ est concave par rapport à la température $T$ et convexe par rapport au volume $V$, à la quantité de matière $n$, etc. ;
* l'enthalpie $H$ est concave par rapport à la pression $P$ et convexe par rapport à l'entropie $S$, la quantité de matière $n$, etc. ;
* l'enthalpie libre $G$ est concave par rapport à la pression $P$ et la température $T$, et convexe par rapport à la quantité de matière $n$, etc.
Les potentiels thermodynamiques $F$, $H$ et $G$ sont donc des fonctions convexes-concaves.

## Conditions d'établissement d'un équilibre stable
### Conditions générales d'évolution et d'équilibre

Soit un système fermé n'échangeant avec l'extérieur que le travail des forces de pression $\delta W=-P\,\mathrm {d} V$ et de la chaleur $\delta Q=T\,\mathrm {d} S+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$. Le terme générique $\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$ représente le travail de forces autres que les forces de pression, ou une réaction chimique, etc., ce travail non récupéré étant dégradé en chaleur. Les variables $X_{i}$ sont toutes extensives (par exemple les quantités des diverses espèces présentes). L'énergie interne $U$ du système varie donc selon :
$\mathrm {d} U=-P\,\mathrm {d} V+T\,\mathrm {d} S+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$
L'entropie $S$ du système est une fonction concave :
$\mathrm {d} ^{2}S=\mathrm {d} {P \over T}\,\mathrm {d} V+\mathrm {d} {1 \over T}\,\mathrm {d} U-\sum _{i}\mathrm {d} {U_{X_{i}} \over T}\,\mathrm {d} X_{i}\leq 0$
et son énergie interne $U$ une fonction convexe :
$\mathrm {d} ^{2}U=-\mathrm {d} P\,\mathrm {d} V+\mathrm {d} T\,\mathrm {d} S+\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} ^{2}S\geq 0$
Le système évolue selon les premier principe et deuxième principe de la thermodynamique :
* Condition générale d'évolution d'un système fermé $\mathrm {d} \Phi =\mathrm {d} U+P\,\mathrm {d} V-T\,\mathrm {d} S=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]\leq 0$
À l'équilibre l'entropie globale $S+S_{\text{ext}}$ n'évolue plus, soit $\mathrm {d} \left[S+S_{\text{ext}}\right]=0$. On a la condition d'équilibre du système fermé :
* Condition générale d'équilibre $\mathrm {d} \Phi =0$

### Équilibre d'un système isolé
#### Conditions d'évolution et d'équilibre d'un système isolé

Un système isolé n'échange ni travail, soit $\delta W=0$, ni chaleur, soit $\delta Q=0$, ni matière avec l'extérieur. Selon le premier principe de la thermodynamique, le système isolé évolue à énergie interne $U$ constante, soit $\mathrm {d} U=0$. En l'absence de tout travail d'une force de pression échangé avec l'extérieur, soit $\delta W=-P\,\mathrm {d} V=0$, le système évolue à volume constant, soit $\mathrm {d} V=0$. En l'absence de tout autre échange avec l'extérieur, l'énergie interne $U_{\text{ext}}$ et l'entropie $S_{\text{ext}}$ de celui-ci ne varient pas, soit $\mathrm {d} U_{\text{ext}}=0$ et $\mathrm {d} S_{\text{ext}}=0$. Le système isolé évolue selon les premier principe et deuxième principe de la thermodynamique :
* Condition d'évolution d'un système isolé $\mathrm {d} \Phi =-T\,\mathrm {d} S=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}\leq 0$
avec :
    $\mathrm {d} ^{2}S=-\sum _{i}\mathrm {d} {U_{X_{i}} \over T}\,\mathrm {d} X_{i}\leq 0$
    $\mathrm {d} ^{2}U=\mathrm {d} T\,\mathrm {d} S+\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\geq 0$
À l'équilibre l'entropie du système n'évolue plus, soit $\mathrm {d} S=0$. On a :
    $\mathrm {d} \Phi =\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=0$
et par conséquent :
    $\mathrm {d} ^{2}S=-{1 \over T}\,\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\leq 0$
    $\mathrm {d} ^{2}U=\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\geq 0$
On a les conditions d'équilibre du système isolé :
* Conditions générales d'équilibre stable d'un système isolé $\mathrm {d} S=0$ et $\mathrm {d} ^{2}S\leq 0$
On étudie les conséquences de la variation d'un paramètre $X_{i}$ quelconque du système à partir de cet équilibre où le volume $V$, l'énergie interne $U$ et l'entropie $S$ du système isolé sont constantes, soit $\mathrm {d} V=0$, $\mathrm {d} U=0$ et $\mathrm {d} S=0$.

#### Principe de l'entropie maximale

Tout au long de la transformation, l'entropie du système isolé évolue selon $\mathrm {d} S>0$ et $\mathrm {d} ^{2}S\leq 0$. À l'équilibre, $\mathrm {d} S=0$ et $\mathrm {d} ^{2}S\leq 0$ se traduisent, à volume $V$ et énergie interne $U$ constants, pour tout paramètre $X_{i}$, par :
    $S_{X_{i}}=-{U_{X_{i}} \over T}=\left({\partial S \over \partial X_{i}}\right)_{V,U,X_{\neq i}}=0$
    $S_{X_{i}X_{i}}=\left({\partial ^{2}S \over {\partial X_{i}}^{2}}\right)_{V,U,X_{\neq i}}\leq 0$
À l'équilibre l'entropie $S$ atteint donc un maximum (voir la figure 1). Toute variation d'un paramètre $X_{i}$ quelconque à partir de cet équilibre, à volume $V$ et énergie interne $U$ constants, ne peut que faire diminuer l'entropie, ce qu'interdit le deuxième principe de la thermodynamique. Le système ne peut plus évoluer, l'équilibre est stable, d'où le principe de l'entropie maximale$^,$$^,$ :
* Principe de l'entropie maximale L'état d'équilibre stable d'un système isolé est caractérisé par une entropie $S$ maximale : $\mathrm {d} S=0$ condition d'équilibre, $\mathrm {d} ^{2}S\leq 0$ condition de stabilité.

#### Principe de l'énergie interne minimale

Tout au long de la transformation, l'énergie interne du système isolé évolue selon $\mathrm {d} U=0$ et $\mathrm {d} ^{2}U\geq 0$. À l'équilibre, ceci se traduit, à volume $V$ et entropie $S$ constants, pour tout paramètre $X_{i}$, par :
    $U_{X_{i}}=-TS_{X_{i}}=\left({\partial U \over \partial X_{i}}\right)_{V,S,X_{\neq i}}=0$
    $U_{X_{i}X_{i}}=\left({\partial ^{2}U \over {\partial X_{i}}^{2}}\right)_{V,S,X_{\neq i}}\geq 0$
L'équilibre correspond donc à un minimum de l'énergie interne $U$ (voir la figure 2). Toute variation d'un paramètre $X_{i}$ quelconque à partir de cet équilibre, à volume $V$ et entropie $S$ constants, ne peut que faire augmenter l'énergie interne, ce qu'interdit le premier principe de la thermodynamique. Le système ne peut plus évoluer, l'équilibre est stable, d'où le principe de l'énergie minimale$^,$ :
* Principe de l'énergie interne minimale L'état d'équilibre stable d'un système isolé est caractérisé par une énergie interne $U$ minimale : $\mathrm {d} U=0$ condition d'équilibre, $\mathrm {d} ^{2}U\geq 0$ condition de stabilité.
Si $X_{i}$ est le seul paramètre variable au cours de la transformation, l'intersection de la surface bleue et du plan jaune sur la figure 1 montre le chemin représentatif de la transformation du système isolé vers l'équilibre E. La transformation s'effectue à volume $V$ constant, l'énergie interne $U$ est constante et convexe, l'entropie $S$ est concave. Selon les conditions initiales de la transformation, $X_{i}$ peut augmenter ou diminuer de telle sorte que $S$ ne fasse qu'augmenter. Lorsque l'équilibre est atteint au point E, $S$ est maximale (figure 1) et $U$ minimale (figure 2) par rapport à $X_{i}$.

### Équilibre d'un système fermé sous contrainte
#### Conditions d'évolution et d'équilibre d'un système fermé

Un système fermé échange du travail et de la chaleur avec l'extérieur, mais pas de matière. On considère qu'il forme avec l'extérieur un système isolé. L'extérieur constitue un réservoir qui permet de maintenir constantes certaines caractéristiques du système fermé étudié. Ce réservoir est supposé suffisamment grand, en regard du système étudié, pour que les échanges avec le système ne modifient pas de façon significative ses propres caractéristiques, qui sont donc considérées comme constantes au cours de la transformation. Par exemple, un réacteur (système étudié) en équilibre de pression et de température avec l'atmosphère (extérieur) reste constamment à pression et température atmosphériques, qui ne varient pas ; l'atmosphère reçoit le travail des forces de pression et la chaleur du réacteur, mais cet apport d'énergie est bien insuffisant pour modifier sa pression et sa température. L'extérieur impose donc des contraintes au système fermé étudié. Au cours de cette transformation, l'une des deux entropies $S$ et $S_{\text{ext}}$ peut décroître, du moment que l'entropie globale $S+S_{\text{ext}}$ du système isolé croît. Le système fermé évolue selon les premier principe et deuxième principe de la thermodynamique :
* énergie interne globale invariante : $\mathrm {d} \left[U+U_{\text{ext}}\right]=0$ entropie globale croissante : $\mathrm {d} \left[S+S_{\text{ext}}\right]\geq 0$
soit :
* Condition générale d'évolution d'un système fermé $\mathrm {d} \Phi =\mathrm {d} U+P\,\mathrm {d} V-T\,\mathrm {d} S=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]\leq 0$
La condition d'évolution induit que $S+S_{\text{ext}}$ et $\Phi$ varient en sens inverses. Le deuxième principe de la thermodynamique interdisant toute diminution de $S+S_{\text{ext}}$, il interdit également toute augmentation de $\Phi$. La convexité du potentiel $\Phi$ peut être déterminée par comparaison de $\mathrm {d} ^{2}\Phi$ et $\mathrm {d} ^{2}U\geq 0$ au cas par cas, selon les contraintes de la transformation.
À l'équilibre, l'entropie globale $S+S_{\text{ext}}$ atteint un maximum, et par conséquent $\Phi$ atteint un minimum. Depuis l'équilibre, toute variation d'un paramètre $X_{i}$ quelconque ne peut induire qu'une diminution de $S+S_{\text{ext}}$ et donc une augmentation de $\Phi$, ce qui est impossible. On a les conditions d'équilibre du système fermé :
* Conditions générales d'équilibre stable d'un système fermé $\mathrm {d} \Phi =0$ et $\mathrm {d} ^{2}\Phi \geq 0$
$\Phi$
Les différentielles de $S+S_{\text{ext}}$ s'écrivent :
    $\mathrm {d} \left[S+S_{\text{ext}}\right]=-\sum _{i}{U_{X_{i}} \over T}\,\mathrm {d} X_{i}$
    $\mathrm {d} ^{2}\left[S+S_{\text{ext}}\right]=-\sum _{i}\mathrm {d} {U_{X_{i}} \over T}\,\mathrm {d} X_{i}={1 \over T^{2}}\,\mathrm {d} T\,\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}-{1 \over T}\,\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\leq 0$
L'entropie globale $S+S_{\text{ext}}$, somme de deux fonctions concaves, est concave. Les différentielles de $\Phi$ s'écrivent :
    $\mathrm {d} \Phi =\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]$
    $\mathrm {d} ^{2}\Phi =\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}=-\mathrm {d} T\,\mathrm {d} \left[S+S_{\text{ext}}\right]-T\,\mathrm {d} ^{2}\left[S+S_{\text{ext}}\right]$
La convexité de $\Phi$ ne peut être déterminée sans plus d'informations sur la transformation. À l'équilibre, l'entropie globale ne varie plus, soit $\mathrm {d} \left[S+S_{\text{ext}}\right]=0$, d'où $\mathrm {d} \Phi =0$ et :
    $\mathrm {d} ^{2}\Phi =-T\,\mathrm {d} ^{2}\left[S+S_{\text{ext}}\right]\geq 0$
À l'équilibre $\Phi$ est convexe et atteint un minimum.

#### À volume et entropie constants

L'énergie interne $U$ d'un système varie de façon générale selon :
    $\mathrm {d} U=-P\,\mathrm {d} V+T\,\mathrm {d} S+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$
    $\mathrm {d} ^{2}U=-\mathrm {d} P\,\mathrm {d} V+\mathrm {d} T\,\mathrm {d} S+\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}$
Pour un système fermé évoluant à volume $V$ et entropie $S$ du système constants, soit $\mathrm {d} V=0$ et $\mathrm {d} S=0$, on a la condition d'évolution :
* Condition d'évolution d'un système fermé à volume et entropie constants $\mathrm {d} \Phi =\mathrm {d} U=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} S_{\text{ext}}\leq 0$
ainsi que :
    $\mathrm {d} ^{2}\Phi =\mathrm {d} ^{2}U=\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\geq 0$
L'entropie $S$ du système est constante, mais l'entropie $S_{\text{ext}}$ de l'extérieur augmente. Le système évolue spontanément dans le sens d'une diminution de l'énergie interne $U$, mais l'énergie interne globale $U+U_{\text{ext}}$ est constante. Dans ces conditions, l'énergie interne $U$ est convexe. Lorsque l'équilibre est atteint, $S_{\text{ext}}$ ne varie plus, soit $\mathrm {d} S_{\text{ext}}=0$. On a les conditions d'équilibre :
* Conditions d'équilibre stable $\mathrm {d} U=0$ et $\mathrm {d} ^{2}U\geq 0$
Par conséquent :
* Dans une transformation à volume $V$ et entropie $S$ constants, l'énergie interne $U$ d'un système fermé diminue. L'équilibre stable est caractérisé par un minimum de $U$.

#### À volume et température constants

On introduit l'énergie libre $F$ :
    $F=U-TS$
    $\mathrm {d} F=-P\,\mathrm {d} V-S\,\mathrm {d} T+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$
    $\mathrm {d} ^{2}F=-\mathrm {d} P\,\mathrm {d} V-\mathrm {d} S\,\mathrm {d} T+\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}$
Pour un système fermé évoluant à volume $V$ et température $T$ constants, soit $\mathrm {d} V=0$ et $\mathrm {d} T=0$, on a la condition d'évolution :
* Condition d'évolution d'un système fermé à volume et température constants $\mathrm {d} \Phi =\mathrm {d} F=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]\leq 0$
ainsi que :
    $\mathrm {d} ^{2}\Phi =\mathrm {d} ^{2}U=\mathrm {d} ^{2}F=\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\geq 0$
L'entropie globale $S+S_{\text{ext}}$ augmente. Le système évolue spontanément dans le sens d'une diminution de l'énergie libre $F$. Dans ces conditions, l'énergie libre $F$ est convexe. Lorsque l'équilibre est atteint, $S+S_{\text{ext}}$ ne varie plus, soit $\mathrm {d} \left[S+S_{\text{ext}}\right]=0$. On a les conditions d'équilibre :
* Conditions d'équilibre stable $\mathrm {d} F=0$ et $\mathrm {d} ^{2}F\geq 0$
Par conséquent :
* Dans une transformation à volume $V$ et température $T$ constants, l'énergie libre $F$ d'un système fermé diminue. L'équilibre stable est caractérisé par un minimum de $F$.

#### À pression et entropie constantes

On introduit l'enthalpie $H$ :
    $H=U+PV$
    $\mathrm {d} H=V\,\mathrm {d} P+T\,\mathrm {d} S+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$
    $\mathrm {d} ^{2}H=\mathrm {d} V\,\mathrm {d} P+\mathrm {d} T\,\mathrm {d} S+\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}$
Pour un système fermé évoluant à pression $P$ et entropie $S$ constantes, soit $\mathrm {d} P=0$ et $\mathrm {d} S=0$, on a la condition d'évolution :
* Condition d'évolution d'un système fermé à pression et entropie constantes $\mathrm {d} \Phi =\mathrm {d} H=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} S_{\text{ext}}\leq 0$
ainsi que :
    $\mathrm {d} ^{2}\Phi =\mathrm {d} ^{2}U=\mathrm {d} ^{2}H=\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\geq 0$
L'entropie $S$ du système est constante, mais l'entropie $S_{\text{ext}}$ de l'extérieur augmente. Le système évolue spontanément dans le sens d'une diminution de l'enthalpie $H$. Dans ces conditions, l'enthalpie $H$ est convexe. Lorsque l'équilibre est atteint, $S_{\text{ext}}$ ne varie plus, soit $\mathrm {d} S_{\text{ext}}=0$. On a les conditions d'équilibre :
* Conditions d'équilibre stable $\mathrm {d} H=0$ et $\mathrm {d} ^{2}H\geq 0$
Par conséquent :
* Dans une transformation à pression $P$ et entropie $S$ constantes, l'enthalpie $H$ d'un système fermé diminue. L'équilibre stable est caractérisé par un minimum de $H$.

#### À pression et température constantes

On introduit l'enthalpie libre $G$ :
    $G=U+PV-TS$
    $\mathrm {d} G=V\,\mathrm {d} P-S\,\mathrm {d} T+\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}$
    $\mathrm {d} ^{2}G=\mathrm {d} V\,\mathrm {d} P-\mathrm {d} S\,\mathrm {d} T+\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}$
Pour un système fermé évoluant à pression $P$ et température $T$ constantes, soit $\mathrm {d} P=0$ et $\mathrm {d} T=0$, on a la condition d'évolution :
* Condition d'évolution d'un système fermé à pression et température constantes $\mathrm {d} \Phi =\mathrm {d} G=\sum _{i}U_{X_{i}}\,\mathrm {d} X_{i}=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]\leq 0$
ainsi que :
    $\mathrm {d} ^{2}\Phi =\mathrm {d} ^{2}U=\mathrm {d} ^{2}G=\sum _{i}\mathrm {d} U_{X_{i}}\,\mathrm {d} X_{i}\geq 0$
L'entropie globale $S+S_{\text{ext}}$ augmente. Le système évolue spontanément dans le sens d'une diminution de l'enthalpie libre $G$. Dans ces conditions, l'enthalpie libre $G$ est convexe. Lorsque l'équilibre est atteint, $S+S_{\text{ext}}$ ne varie plus, soit $\mathrm {d} \left[S+S_{\text{ext}}\right]=0$. On a les conditions d'équilibre :
* Conditions d'équilibre stable $\mathrm {d} G=0$ et $\mathrm {d} ^{2}G\geq 0$
Par conséquent :
* Dans une transformation à pression $P$ et température $T$ constantes, l'enthalpie libre $G$ d'un système fermé diminue. L'équilibre stable est caractérisé par un minimum de $G$.

## Implications
### Conditions de stabilité
#### Expression des conditions de stabilité

La concavité de l'entropie et la convexité de l'énergie interne induisent que leurs dérivées secondes, ainsi que celles des autres potentiels thermodynamiques, doivent répondre à certaines contraintes, notamment quant à leur signe (voir la section Notations et outils mathématiques). L'équilibre d'un système dans lequel l'une de ces dérivées secondes ne vérifie pas ces contraintes ne peut être stable, car il enfreint le deuxième principe de la thermodynamique. Ainsi les critères de convexité de l'entropie et des potentiels thermodynamiques induisent-ils les conditions de stabilité d'un système. Les conditions données ci-dessous sont toutefois nécessaires, mais pas suffisantes.
Dans les expressions ci-dessous, le théorème de Schwarz permet d'écrire, pour l'entropie comme pour tout potentiel thermodynamique, et pour deux variables quelconques : $\Phi _{XY}=\Phi _{YX}$. Ainsi, les conditions données pour $\Phi _{XY}$ sont également valables pour $\Phi _{YX}$.
Les conditions de stabilité relatives aux dérivées partielles secondes de l'entropie sont :
* Conditions de stabilité selon l'entropie (csS1) : $\left({\partial ^{2}S \over {\partial V}^{2}}\right)_{U,n}\leq 0$ ; (csS2) : $\left({\partial ^{2}S \over {\partial U}^{2}}\right)_{V,n}\leq 0$ ; (csS3) : $\left({\partial ^{2}S \over {\partial n_{i}}^{2}}\right)_{V,U,n_{\neq i}}\leq 0$ (csS4) : ${\left({\partial ^{2}S \over \partial V\partial U}\right)_{n}}^{2}\leq \left({\partial ^{2}S \over {\partial V}^{2}}\right)_{U,n}\,\left({\partial ^{2}S \over {\partial U}^{2}}\right)_{V,n}$ ; (csS5) : ${\left({\partial ^{2}S \over \partial V\partial n_{i}}\right)_{U,n_{\neq i}}}^{2}\leq \left({\partial ^{2}S \over {\partial V}^{2}}\right)_{U,n}\,\left({\partial ^{2}S \over {\partial n_{i}}^{2}}\right)_{V,U,n_{\neq i}}$ (csS6) : ${\left({\partial ^{2}S \over \partial U\partial n_{i}}\right)_{V,n_{\neq i}}}^{2}\leq \left({\partial ^{2}S \over {\partial U}^{2}}\right)_{V,n}\,\left({\partial ^{2}S \over {\partial n_{i}}^{2}}\right)_{V,U,n_{\neq i}}$ ; (csS7) : ${\left({\partial ^{2}S \over \partial n_{i}\partial n_{j}}\right)_{V,U,n_{\neq i,j}}}^{2}\leq \left({\partial ^{2}S \over {\partial n_{i}}^{2}}\right)_{V,U,n_{\neq i}}\,\left({\partial ^{2}S \over {\partial n_{j}}^{2}}\right)_{V,U,n_{\neq j}}$
L'entropie est une fonction concave par rapport à toutes ses variables, qui sont toutes extensives.
Les conditions de stabilité relatives aux dérivées partielles secondes des potentiels thermodynamiques sont :
* Conditions de stabilité selon l'énergie interne (csU1) : $\left({\partial ^{2}U \over {\partial V}^{2}}\right)_{S,n}\geq 0$ ; (csU2) : $\left({\partial ^{2}U \over {\partial S}^{2}}\right)_{V,n}\geq 0$ ; (csU3) : $\left({\partial ^{2}U \over {\partial n_{i}}^{2}}\right)_{V,S,n_{\neq i}}\geq 0$ (csU4) : ${\left({\partial ^{2}U \over \partial V\partial S}\right)_{n}}^{2}\leq \left({\partial ^{2}U \over {\partial V}^{2}}\right)_{S,n}\,\left({\partial ^{2}U \over {\partial S}^{2}}\right)_{V,n}$ ; (csU5) : ${\left({\partial ^{2}U \over \partial V\partial n_{i}}\right)_{S,n_{\neq i}}}^{2}\leq \left({\partial ^{2}U \over {\partial V}^{2}}\right)_{S,n}\,\left({\partial ^{2}U \over {\partial n_{i}}^{2}}\right)_{V,S,n_{\neq i}}$ (csU6) : ${\left({\partial ^{2}U \over \partial S\partial n_{i}}\right)_{V,n_{\neq i}}}^{2}\leq \left({\partial ^{2}U \over {\partial S}^{2}}\right)_{V,n}\,\left({\partial ^{2}U \over {\partial n_{i}}^{2}}\right)_{V,S,n_{\neq i}}$ ; (csU7) : ${\left({\partial ^{2}U \over \partial n_{i}\partial n_{j}}\right)_{V,S,n_{\neq i,j}}}^{2}\leq \left({\partial ^{2}U \over {\partial n_{i}}^{2}}\right)_{V,S,n_{\neq i}}\,\left({\partial ^{2}U \over {\partial n_{j}}^{2}}\right)_{V,S,n_{\neq j}}$
* Conditions de stabilité selon l'énergie libre (csF1) : $\left({\partial ^{2}F \over {\partial V}^{2}}\right)_{T,n}\geq 0$ ; (csF2) : $\left({\partial ^{2}F \over {\partial T}^{2}}\right)_{V,n}\leq 0$ ; (csF3) : $\left({\partial ^{2}F \over {\partial n_{i}}^{2}}\right)_{V,T,n_{\neq i}}\geq 0$ (csF4) : ${\left({\partial ^{2}F \over \partial V\partial n_{i}}\right)_{T,n_{\neq i}}}^{2}\leq \left({\partial ^{2}F \over {\partial V}^{2}}\right)_{T,n}\,\left({\partial ^{2}F \over {\partial n_{i}}^{2}}\right)_{V,T,n_{\neq i}}$ ; (csF5) : ${\left({\partial ^{2}F \over \partial n_{i}\partial n_{j}}\right)_{V,T,n_{\neq i,j}}}^{2}\leq \left({\partial ^{2}F \over {\partial n_{i}}^{2}}\right)_{V,T,n_{\neq i}}\,\left({\partial ^{2}F \over {\partial n_{j}}^{2}}\right)_{V,T,n_{\neq j}}$
* Conditions de stabilité selon l'enthalpie (csH1) : $\left({\partial ^{2}H \over {\partial P}^{2}}\right)_{S,n}\leq 0$ ; (csH2) : $\left({\partial ^{2}H \over {\partial S}^{2}}\right)_{P,n}\geq 0$ ; (csH3) : $\left({\partial ^{2}H \over {\partial n_{i}}^{2}}\right)_{P,S,n_{\neq i}}\geq 0$ (csH4) : ${\left({\partial ^{2}H \over \partial S\partial n_{i}}\right)_{P,n_{\neq i}}}^{2}\leq \left({\partial ^{2}H \over {\partial S}^{2}}\right)_{P,n}\,\left({\partial ^{2}H \over {\partial n_{i}}^{2}}\right)_{P,S,n_{\neq i}}$ ; (csH5) : ${\left({\partial ^{2}H \over \partial n_{i}\partial n_{j}}\right)_{P,S,n_{\neq i,j}}}^{2}\leq \left({\partial ^{2}H \over {\partial n_{i}}^{2}}\right)_{P,S,n_{\neq i}}\,\left({\partial ^{2}H \over {\partial n_{j}}^{2}}\right)_{P,S,n_{\neq j}}$
* Conditions de stabilité selon l'enthalpie libre (csG1) : $\left({\partial ^{2}G \over {\partial P}^{2}}\right)_{T,n}\leq 0$ ; (csG2) : $\left({\partial ^{2}G \over {\partial T}^{2}}\right)_{P,n}\leq 0$ ; (csG3) : $\left({\partial ^{2}G \over {\partial n_{i}}^{2}}\right)_{P,T,n_{\neq i}}\geq 0$ (csG4) : ${\left({\partial ^{2}G \over \partial P\partial T}\right)_{n}}^{2}\leq \left({\partial ^{2}G \over {\partial P}^{2}}\right)_{T,n}\,\left({\partial ^{2}G \over {\partial T}^{2}}\right)_{P,n}$ ; (csG5) : ${\left({\partial ^{2}G \over \partial n_{i}\partial n_{j}}\right)_{P,T,n_{\neq i,j}}}^{2}\leq \left({\partial ^{2}G \over {\partial n_{i}}^{2}}\right)_{P,T,n_{\neq i}}\,\left({\partial ^{2}G \over {\partial n_{j}}^{2}}\right)_{P,T,n_{\neq j}}$

#### Conclusions générales

En résumé, pour tout potentiel thermodynamique $\Phi =\Phi \!\left(Z\right)$, $Z$ étant les variables naturelles de $\Phi$, dont $X$ et $Y$ :
* si $X$ est intensive (pression, température) :
       $\Phi _{XX}=\left({\partial ^{2}\Phi \over {\partial X}^{2}}\right)_{Z\neq X}\leq 0$
    Les potentiels thermodynamiques sont concaves par rapport à leurs variables naturelles intensives$^,$ ;
* si $X$ est extensive (volume, entropie, quantité de matière) :
       $\Phi _{XX}=\left({\partial ^{2}\Phi \over {\partial X}^{2}}\right)_{Z\neq X}\geq 0$
    Les potentiels thermodynamiques sont convexes par rapport à leurs variables naturelles extensives$^,$ ;
* si $X$ et $Y$ sont simultanément intensives ou extensives :
       ${\left({\partial ^{2}\Phi \over \partial X\partial Y}\right)_{Z\neq X,Y}}^{2}\leq \left({\partial ^{2}\Phi \over {\partial X}^{2}}\right)_{Z\neq X}\,\left({\partial ^{2}\Phi \over {\partial Y}^{2}}\right)_{Z\neq Y}$, soit ${\Phi _{XY}}^{2}\leq \Phi _{XX}\Phi _{YY}$.
    Cette contrainte donne une indication sur la valeur absolue de la dérivée croisée, mais pas sur son signe ;
* si $X$ est intensive et $Y$ extensive, il n'y a pas de condition sur la dérivée croisée $\Phi _{XY}=\left({\partial ^{2}\Phi \over \partial X\partial Y}\right)_{Z\neq X,Y}$.
Soit $\Phi _{X}$ la variable conjuguée de $X$ :
$\Phi _{X}=\left({\partial \Phi \over \partial X}\right)_{Z\neq X}$
On note $\mathrm {d} \Phi _{X}|_{X}$ la variation de $\Phi _{X}$ induite par une modification de $X$ à $Z\neq X$ constantes :
    $\left({\partial \Phi _{X} \over \partial X}\right)_{Z\neq X}=\left({\partial ^{2}\Phi \over {\partial X}^{2}}\right)_{Z\neq X}=\Phi _{XX}$
    $\mathrm {d} \Phi _{X}|_{X}=\Phi _{XX}\,\mathrm {d} X$
On a donc, pour toutes variables $Z\neq X$ constantes :
$\mathrm {d} \Phi _{X}|_{X}\,\mathrm {d} X=\Phi _{XX}\,\left(\mathrm {d} X\right)^{2}$
Les conditions de stabilité donnent $\Phi _{XX}\leq 0$ si $X$ est intensive, et $\Phi _{XX}\geq 0$ si $X$ est extensive. Par conséquent tout système évolue, à $Z\neq X$ constantes, selon :
* Si $X$ est intensive, alors $\mathrm {d} \Phi _{X}|_{X}\,\mathrm {d} X\leq 0$. Si $X$ est extensive, alors $\mathrm {d} \Phi _{X}|_{X}\,\mathrm {d} X\geq 0$.
En conclusion :
| La modification d'une variable $X$ intensive induit une modification opposée de sa variable conjuguée extensive $\Phi _{X}$. | La modification d'une variable $X$ extensive induit une modification dans le même sens de sa variable conjuguée intensive $\Phi _{X}$. |
Exemple 1 - Évolutions réciproques de la pression et du volume.
   * Selon (edeV) et l'enthalpie libre $G$, la variable conjuguée de la pression $P$ est le volume : $G_{P}=V$. La pression étant une variable intensive, une augmentation de pression, à température et composition constantes, induit une diminution du volume, variable extensive. Il s'agit de la conséquence de (csG1).
   * Selon (edeP) et l'énergie libre $F$, la variable conjuguée du volume $V$ est l'opposé de la pression : $F_{V}=-P$. Le volume étant une variable extensive, une augmentation du volume, à température et composition constantes, induit une augmentation de $-P$, variable intensive, soit une diminution de la pression. Il s'agit de la conséquence de (csF1).
Soit une variable $Y$ de $\Phi$ différente de $X$. On note $\mathrm {d} \Phi _{X}|_{Y}$ la variation de $\Phi _{X}$ induite par une modification de $Y$ à $Z\neq Y$ constantes :
    $\left({\partial \Phi _{X} \over \partial Y}\right)_{Z\neq Y}=\left({\partial ^{2}\Phi \over \partial Y\partial X}\right)_{Z\neq X,Y}=\Phi _{YX}=\Phi _{XY}$
    $\mathrm {d} \Phi _{X}|_{Y}=\Phi _{XY}\,\mathrm {d} Y$
On a donc, pour toutes variables $Z\neq Y$ constantes :
$\mathrm {d} \Phi _{X}|_{Y}\,\mathrm {d} Y=\Phi _{XY}\,\left(\mathrm {d} Y\right)^{2}$
$X$ et $Y$ étant différentes, le signe de $\Phi _{XY}$ n'est pas imposé par les conditions de stabilité, que ces deux variables soient simultanément ou non extensives ou intensives. Le signe de $\mathrm {d} \Phi _{X}|_{Y}\,\mathrm {d} Y$ n'est donc pas défini, il peut être positif ou négatif. En conclusion :
* Il est impossible de prévoir le sens de variation de $\Phi _{X}$ en fonction d'une modification d'une autre variable $Y$ du potentiel thermodynamique $\Phi$.
Exemple 2 - Évolution du volume en fonction de la quantité d'un produit.
    Le signe du volume molaire partiel ${\bar {V}}_{i}=\left({\partial V \over \partial n_{i}}\right)_{P,T,n_{\neq i}}$ d'un corps $i$ dans un mélange n'est pas imposé par les conditions de stabilité. Avec (edeV), soit $G_{P}=V$, on a ${\bar {V}}_{i}=\left({\partial ^{2}G \over \partial n_{i}\partial P}\right)_{T,n_{\neq i}}$. Il n'y a pas de contrainte de signe sur cette dérivée croisée : on ne peut prévoir la variation de $G_{P}$ en fonction de la quantité $n_{i}$ de l'espèce $i$. Le volume molaire partiel d'un corps $i$ peut être aussi bien positif (le volume d'un mélange augmente lorsque l'on ajoute $i$) que négatif (le volume d'un mélange diminue lorsque l'on ajoute $i$), ce qui est toutefois assez rare. Par exemple, lorsque l'on dissout 0,1 mole de sulfate de magnésium MgSO$_4$ dans un litre d'eau, le volume de la solution obtenue est inférieur à un litre. Le volume se contracte par addition de MgSO$_4$, son volume molaire partiel est donc négatif.
    Par ailleurs, le théorème de Schwarz donne ${\bar {V}}_{i}=\left({\partial \mu _{i} \over \partial P}\right)_{T,n}$ avec la définition (mu) $G_{n_{i}}=\mu _{i}$. Il est donc également impossible de prévoir l'évolution d'un potentiel chimique $\mu _{i}$ en fonction de la pression $P$.

#### Coefficients calorimétriques et thermoélastiques

Parmi les coefficients calorimétriques et thermoélastiques sont définis :
* la capacité thermique isochore $C_{V}$ ; étant donné les équations d'état (edeS) et (edeT) on a :
       $C_{V}=T\left({\partial S \over \partial T}\right)_{V,n}=-T\left({\partial ^{2}F \over {\partial T}^{2}}\right)_{V,n}={T \over \left({\partial T \over \partial S}\right)_{V,n}}={T \over \left({\partial ^{2}U \over {\partial S}^{2}}\right)_{V,n}}$
    Les conditions de stabilité (csF2) et (csU2) donnent$^,$ :
* (csCV) : $C_{V}\geq 0$
* le coefficient de compressibilité isentropique $\chi _{S}$ ; étant donné les équations d'état (edeV) et (edeP) on a :
       $\chi _{S}=-{1 \over V}\left({\partial V \over \partial P}\right)_{S,n}=-{1 \over V}\left({\partial ^{2}H \over {\partial P}^{2}}\right)_{S,n}=-{1 \over V\left({\partial P \over \partial V}\right)_{S,n}}={1 \over V\left({\partial ^{2}U \over {\partial V}^{2}}\right)_{S,n}}$
    Les conditions de stabilité (csH1) et (csU1) donnent$^,$ :
* (csChiS) : $\chi _{S}\geq 0$
* un coefficient sans nom $\lambda$ ; étant donné l'équation d'état (edeP) on a :
       $\lambda =T\left({\partial S \over \partial P}\right)_{V,n}={T \over \left({\partial P \over \partial S}\right)_{V,n}}=-{T \over \left({\partial ^{2}U \over \partial S\partial V}\right)_{n}}$
    La condition de stabilité (csU4) donne :
* (csLambda) : $TVC_{V}\chi _{S}\leq \lambda ^{2}$
* la capacité thermique isobare $C_{P}$ ; étant donné les équations d'état (edeS) et (edeT) on a :
       $C_{P}=T\left({\partial S \over \partial T}\right)_{P,n}=-T\left({\partial ^{2}G \over {\partial T}^{2}}\right)_{P,n}={T \over \left({\partial T \over \partial S}\right)_{P,n}}={T \over \left({\partial ^{2}H \over {\partial S}^{2}}\right)_{P,n}}$
    Les conditions de stabilité (csG2) et (csH2) donnent$^,$ :
* (csCP) : $C_{P}\geq 0$
* le coefficient de compressibilité isotherme $\chi _{T}$ ; étant donné les équations d'état (edeV) et (edeP) on a :
       $\chi _{T}=-{1 \over V}\left({\partial V \over \partial P}\right)_{T,n}=-{1 \over V}\left({\partial ^{2}G \over {\partial P}^{2}}\right)_{T,n}=-{1 \over V\left({\partial P \over \partial V}\right)_{T,n}}={1 \over V\left({\partial ^{2}F \over {\partial V}^{2}}\right)_{T,n}}$
    Les conditions de stabilité (csG1) et (csF1) donnent$^,$ :
* (csChiT) : $\chi _{T}\geq 0$
* le coefficient de dilatation isobare $\alpha$ ; étant donné l'équation d'état (edeV) on a :
       $\alpha ={1 \over V}\left({\partial V \over \partial T}\right)_{P,n}={1 \over V}\left({\partial ^{2}G \over \partial T\partial P}\right)_{n}$
    La condition de stabilité (csG4) donne :
* (csAlpha) : $TV\alpha ^{2}\leq C_{P}\chi _{T}$
Les capacités thermiques $C_{V}$ et $C_{P}$ positives correspondent aux observations communes : un corps absorbe de la chaleur lorsque sa température augmente et en restitue lorsqu'elle diminue. De même, le volume d'un corps diminue sous l'effet de la pression : les coefficients de compressibilité $\chi _{S}$ et $\chi _{T}$ sont positifs. La thermodynamique n'interdit pas que ces coefficients soient négatifs, mais un corps présentant de telles propriétés serait instable considéré seul car il diminuerait l'entropie, en contradiction avec le deuxième principe de la thermodynamique : une telle situation est donc difficilement observable. Cependant des coefficients négatifs peuvent être observés dans un contexte impliquant des phénomènes compensant cette instabilité. En physique stellaire la stabilité des étoiles est expliquée par une capacité thermique négative due à l'attraction gravitationnelle entre ses constituants. Une étoile génère par fusion nucléaire plus d'énergie qu'elle ne peut en rayonner, ce qui, avec une capacité thermique positive, induirait une telle accumulation de chaleur, et donc une telle augmentation de température, que l'étoile serait instable et mourrait rapidement. La capacité thermique négative permet d'accumuler la chaleur tout en maintenant une température stable. D'autre part, des coefficients de compressibilité négatifs ont été observés sur des mousses métalliques et des cristaux composés d'eau et de méthanol, ces phénomènes étant expliqués par l'architecture des cristaux à l'échelle moléculaire$^,$$^,$.
Le volume d'un corps augmente généralement sous l'effet d'une augmentation de la température, aussi le coefficient de dilatation isobare $\alpha$ est-il le plus souvent positif. Néanmoins, la relation (csAlpha) n'impose pas le signe de ce coefficient, qui peut donc être négatif pour un corps stable. L'eau liquide en est un exemple entre 0 °C et 4 °C sous 1 atm : une augmentation de la température provoque une contraction du volume, d'où un maximum de densité à 4 °C, constituant une anomalie dilatométrique. Si l'on considère la variation du volume en fonction de la pression, de la température et de la composition, on a :
    $\mathrm {d} V=\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}\left({\partial V \over \partial n_{i}}\right)_{P,T,n_{\neq i}}\,\mathrm {d} n_{i}=-V\chi _{T}\,\mathrm {d} P+V\alpha \,\mathrm {d} T+\sum _{i=1}^{N}{\bar {V}}_{i}\,\mathrm {d} n_{i}$
Les conditions de stabilité donnent $\chi _{T}\geq 0$ mais ne donnent pas le signe de $\alpha$ et ${\bar {V}}_{i}$, qui peuvent donc être positifs ou, plus rarement, négatifs.

#### Potentiel chimique

Les conditions de stabilité (csU3), (csF3), (csH3), (csG3), avec les définitions (mu), donnent respectivement :
* Conditions de stabilité (csMu1) : $\left({\partial \mu _{i} \over \partial n_{i}}\right)_{V,S,n_{\neq i}}\geq 0$ ; (csMu2) : $\left({\partial \mu _{i} \over \partial n_{i}}\right)_{V,T,n_{\neq i}}\geq 0$ (csMu3) : $\left({\partial \mu _{i} \over \partial n_{i}}\right)_{P,S,n_{\neq i}}\geq 0$ ; (csMu4) : $\left({\partial \mu _{i} \over \partial n_{i}}\right)_{P,T,n_{\neq i}}\geq 0$
Le potentiel chimique $\mu _{i}$ d'une espèce chimique $i$ quelconque 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.
Les conditions de stabilité (csU5), (csU6), (csF4) et (csH4), avec les définitions (mu), les équations d'état (edeP) et (edeT), et le théorème de Schwarz, donnent respectivement, pour toute espèce chimique $i$ :
* Conditions de stabilité (csMu5) : ${\left({\partial \mu _{i} \over \partial V}\right)_{S,n}}^{2}={\left({\partial P \over \partial n_{i}}\right)_{V,S,n_{\neq i}}}^{2}\leq {1 \over V\chi _{S}}\,\left({\partial \mu _{i} \over \partial n_{i}}\right)_{V,S,n_{\neq i}}$ (csMu6) : ${\left({\partial \mu _{i} \over \partial S}\right)_{V,n}}^{2}={\left({\partial T \over \partial n_{i}}\right)_{V,S,n_{\neq i}}}^{2}\leq {T \over C_{V}}\,\left({\partial \mu _{i} \over \partial n_{i}}\right)_{V,S,n_{\neq i}}$ (csMu7) : ${\left({\partial \mu _{i} \over \partial V}\right)_{T,n}}^{2}={\left({\partial P \over \partial n_{i}}\right)_{V,T,n_{\neq i}}}^{2}\leq {1 \over V\chi _{T}}\,\left({\partial \mu _{i} \over \partial n_{i}}\right)_{V,T,n_{\neq i}}$ (csMu8) : ${\left({\partial \mu _{i} \over \partial S}\right)_{P,n}}^{2}={\left({\partial T \over \partial n_{i}}\right)_{P,S,n_{\neq i}}}^{2}\leq {T \over C_{P}}\,\left({\partial \mu _{i} \over \partial n_{i}}\right)_{P,S,n_{\neq i}}$
Ces conditions montrent, entre autres, qu'il est impossible de prévoir la variation de la pression $P$ ou de la température $T$ d'un système en fonction de la variation (ajout ou extraction en l'absence de réaction) de la quantité $n_{i}$ de l'un de ses constituants $i$. Autrement dit, lors de l'ajout d'un constituant $i$ à un système, toutes choses étant égales par ailleurs, la pression et la température peuvent indifféremment augmenter ou diminuer dans un système stable.
Les conditions de stabilité (csU7), (csF5), (csH5) et (csG5), avec les définitions (mu) et le théorème de Schwarz, donnent respectivement, pour deux espèces chimiques $i$ et $j$ quelconques :
* Conditions de stabilité (csMu9) : ${\left({\partial \mu _{i} \over \partial n_{j}}\right)_{V,S,n_{\neq j}}}^{2}={\left({\partial \mu _{j} \over \partial n_{i}}\right)_{V,S,n_{\neq i}}}^{2}\leq \left({\partial \mu _{i} \over \partial n_{i}}\right)_{V,S,n_{\neq j}}\,\left({\partial \mu _{j} \over \partial n_{j}}\right)_{V,S,n_{\neq j}}$ (csMu10) : ${\left({\partial \mu _{i} \over \partial n_{j}}\right)_{V,T,n_{\neq j}}}^{2}={\left({\partial \mu _{j} \over \partial n_{i}}\right)_{V,T,n_{\neq i}}}^{2}\leq \left({\partial \mu _{i} \over \partial n_{i}}\right)_{V,T,n_{\neq i}}\,\left({\partial \mu _{j} \over \partial n_{j}}\right)_{V,T,n_{\neq j}}$ (csMu11) : ${\left({\partial \mu _{i} \over \partial n_{j}}\right)_{P,S,n_{\neq j}}}^{2}={\left({\partial \mu _{j} \over \partial n_{i}}\right)_{P,S,n_{\neq i}}}^{2}\leq \left({\partial \mu _{i} \over \partial n_{i}}\right)_{P,S,n_{\neq i}}\,\left({\partial \mu _{j} \over \partial n_{j}}\right)_{P,S,n_{\neq j}}$ (csMu12) : ${\left({\partial \mu _{i} \over \partial n_{j}}\right)_{P,T,n_{\neq j}}}^{2}={\left({\partial \mu _{j} \over \partial n_{i}}\right)_{P,T,n_{\neq i}}}^{2}\leq \left({\partial \mu _{i} \over \partial n_{i}}\right)_{P,T,n_{\neq i}}\,\left({\partial \mu _{j} \over \partial n_{j}}\right)_{P,T,n_{\neq j}}$
Ces conditions montrent qu'il est impossible de prévoir l'évolution du potentiel chimique $\mu _{i}$ d'un constituant $i$ dans un mélange en fonction de la variation de la quantité $n_{j}$ d'un autre constituant $j$, toutes choses étant égales par ailleurs.

### Équilibre d'un système hétérogène
#### Conditions générales d'équilibre

On considère un système isolé hétérogène : il peut être subdivisé en plusieurs sous-systèmes ayant chacun sa propre pression, sa propre température et ses propres potentiels chimiques (ces sous-systèmes sont en état d'équilibre thermodynamique local). Il peut s'agir par exemple de phases différentes, mais il peut également s'agir de deux zones distinctes au sein d'une même phase (par exemple une barre de métal ayant des températures différentes aux deux extrémités). Soit deux sous-systèmes notés $\alpha$ et $\beta$. Le système global isolé subit une transformation afin d'atteindre l'équilibre. On écrit pour chacun des deux sous-systèmes 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 $X^{\phi }$ la grandeur $X$ relative au sous-système $\phi$. Soit, pour le système isolé global, les fonctions globales énergie interne globale $U=U^{\alpha }+U^{\beta }$, entropie globale $S=S^{\alpha }+S^{\beta }$, volume global $V=V^{\alpha }+V^{\beta }$ et quantité globale $n_{i}=n_{i}^{\alpha }+n_{i}^{\beta }$ de chacun des $N$ corps $i$. On a les relations sur le système global isolé, en l'absence de réaction chimique :
* $\mathrm {d} U=\mathrm {d} U^{\alpha }+\mathrm {d} U^{\beta }=0$ (premier principe de la thermodynamique) ;
* $\mathrm {d} S=\mathrm {d} S^{\alpha }+\mathrm {d} S^{\beta }\geq 0$ (deuxième principe de la thermodynamique) ;
* $\mathrm {d} V=\mathrm {d} V^{\alpha }+\mathrm {d} V^{\beta }=0$ (absence de travail des forces de pression) ;
* $\mathrm {d} n_{i}=\mathrm {d} n_{i}^{\alpha }+\mathrm {d} n_{i}^{\beta }=0$ pour tout $i\in [1,\cdots ,N]$ (conservation de la matière).
Les grandeurs globales du système global isolé (de l'ensemble constitué par les deux sous-systèmes) sont constantes (hormis l'entropie), mais les grandeurs propres à chacun des deux sous-systèmes, elles, peuvent varier au cours de la transformation, leurs variations étant opposées. Par exemple, le volume global du système isolé, $V=V^{\alpha }+V^{\beta }$, est constant, mais les volumes des sous-systèmes $V^{\alpha }$ et $V^{\beta }$ peuvent varier et $\mathrm {d} V^{\beta }=-\mathrm {d} V^{\alpha }$.
On obtient les relations globales d'évolution du système isolé :
* Évolution $\mathrm {d} U=-\left(P^{\alpha }-P^{\beta }\right)\,\mathrm {d} V^{\alpha }+T^{\alpha }\,\mathrm {d} S^{\alpha }+T^{\beta }\,\mathrm {d} S^{\beta }+\sum _{i=1}^{N}\left(\mu _{i}^{\alpha }-\mu _{i}^{\beta }\right)\,\mathrm {d} n_{i}^{\alpha }=0$ $\mathrm {d} S=\mathrm {d} S^{\alpha }+\mathrm {d} S^{\beta }\geq 0$
L'équilibre est atteint lorsque l'entropie globale du système atteint un maximum. Les variations d'entropie des deux sous-systèmes 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=-\left(P^{\alpha }-P^{\beta }\right)\,\mathrm {d} V^{\alpha }+\left(T^{\alpha }-T^{\beta }\right)\,\mathrm {d} S^{\alpha }+\sum _{i=1}^{N}\left(\mu _{i}^{\alpha }-\mu _{i}^{\beta }\right)\,\mathrm {d} n_{i}^{\alpha }=0$ $\mathrm {d} S=\mathrm {d} S^{\alpha }+\mathrm {d} S^{\beta }=0$
Quand les deux sous-systèmes sont à l'équilibre, cela ne signifie pas qu'il n'y a pas d'échanges entre eux. Les deux sous-systèmes 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 est qualifié de dynamique ou stationnaire (voir la section Échelle macroscopique, équilibre 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 sous-systèmes à 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 (voir la figure 1), alors l'énergie globale du système isolé (constante) correspond à un minimum par rapport aux diverses variables du système (voir la figure 2), ce qui implique que :
* $\left({\partial U \over \partial V^{\alpha }}\right)_{S^{\alpha }\!,\,n^{\alpha }}=-\left(P^{\alpha }-P^{\beta }\right)=0$ ;
* $\left({\partial U \over \partial S^{\alpha }}\right)_{V^{\alpha }\!,\,n^{\alpha }}=T^{\alpha }-T^{\beta }=0$ ;
* $\left({\partial U \over \partial n_{i}^{\alpha }}\right)_{V^{\alpha }\!,\,S^{\alpha }\!,\,n_{\neq 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 sous-systèmes $\alpha$ et $\beta$ :
| Équilibre mécanique | homogénéité de la pression $P$ |
| Équilibre thermique | homogénéité de la température $T$ |
| Équilibre diffusif | homogénéité des potentiels chimiques pour tout corps $i$ : $\mu _{i}^{\alpha }=\mu _{i}^{\beta }$ |
Si le système contient un troisième sous-système $\gamma$, celui-ci doit être à l'équilibre avec les sous-systèmes $\alpha$ et $\beta$ : il sera par conséquent à l'équilibre aux mêmes pression, température et potentiels chimiques. Ces conditions sont donc généralisables à plus de deux sous-systèmes en équilibre (gaz-liquide-liquide ou gaz-liquide-solide par exemple, ou multiples phases liquides ou solides).

##### 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 }\neq \mu _{2}^{\alpha }\neq \cdots \neq \mu _{N}^{\alpha }$.
En conséquence, si deux sous-systèmes d'un système isolé ont des températures différentes, celles-ci tendent à s'homogénéiser par transfert thermique (transfert de chaleur). De même, si à un instant donné une espèce chimique $i$ quelconque a des potentiels chimiques différents en deux points de ce système, alors par diffusion les potentiels chimiques tendent à s'homogénéiser. Par exemple, au sein d'une même phase, des concentrations ou des fractions molaires ou massiques tendent à s'égaliser en tout point. Dans la plupart des cas, l'équilibre mécanique s'établit très rapidement, suivi par l'équilibre thermique. L'équilibre diffusif, ou chimique, peut être beaucoup plus lent à atteindre.

#### Le transfert thermique

Lorsqu'un milieu est hétérogène en température, celle-ci n'est pas identique en tout point du milieu. Spontanément, de la chaleur migre des zones de température élevée vers les zones de température plus basse : l'entropie du système est ainsi maximisée, conformément au deuxième principe de la thermodynamique. Ce phénomène est appelé transfert thermique.
On considère un système isolé divisé en deux compartiments, $\alpha$ et $\beta$, séparés par une paroi fixe diathermane$^,$. Le seul échange possible entre ces deux compartiments est donc un transfert thermique, soit un échange de chaleur. Le volume du système global ne varie pas en l'absence de travail des forces de pression avec l'extérieur, et les volumes des compartiments sont également invariants par la fixité de la paroi.
Il n'est fait aucune hypothèse quant aux pressions, compositions et phases des compartiments, qui peuvent donc être différentes d'un compartiment à l'autre. Les pressions peuvent évoluer au cours de la transformation. On suppose que les compositions sont invariantes en l'absence de réaction chimique.
La variation d'énergie de chacun des compartiments vaut, en considérant que les transformations n'impliquent que de la chaleur :
    $\mathrm {d} U^{\alpha }=\delta Q^{\alpha }=T^{\alpha }\,\mathrm {d} S^{\alpha }$
    $\mathrm {d} U^{\beta }=\delta Q^{\beta }=T^{\beta }\,\mathrm {d} S^{\beta }$
Le système étant isolé, on a les relations :
* $\mathrm {d} U=\mathrm {d} U^{\alpha }+\mathrm {d} U^{\beta }=0$ (premier principe, énergie globale $U=U^{\alpha }+U^{\beta }$) ;
* $\mathrm {d} S=\mathrm {d} S^{\alpha }+\mathrm {d} S^{\beta }\geq 0$ (deuxième principe, entropie globale $S=S^{\alpha }+S^{\beta }$).
Les températures absolues $T^{\alpha }$ et $T^{\beta }$ étant positives, la première relation induit que les entropies $S^{\alpha }$ et $S^{\beta }$ des deux compartiments varient en sens inverses : $T^{\beta }\,\mathrm {d} S^{\beta }=-T^{\alpha }\,\mathrm {d} S^{\alpha }$. En conséquence $Q^{\alpha }$ et $Q^{\beta }$ varient également en sens inverses : si l'une augmente, l'autre diminue. On obtient finalement :
* À énergie interne constante $T^{\beta }\,\mathrm {d} S=-\left(T^{\alpha }-T^{\beta }\right)\,\mathrm {d} S^{\alpha }$
On suppose que la température du compartiment $\alpha$ est plus importante que celle du compartiment $\beta$, soit $T^{\alpha }>T^{\beta }$. La température absolue $T^{\beta }$ étant 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} S^{\alpha }<0$, d'où $\mathrm {d} S^{\beta }>0$, $\delta Q^{\alpha }<0$ et $\delta Q^{\beta }>0$ : la chaleur diminue dans le compartiment $\alpha$ et augmente dans le compartiment $\beta$. En conséquence, la chaleur migre du compartiment $\alpha$, dans lequel la température est la plus élevée, vers le compartiment $\beta$, de moindre température. Cette relation ne permet pas de déterminer les variations des températures $T^{\alpha }$ et $T^{\beta }$ des deux compartiments.
Les deux compartiments évoluant à volume constant et en l'absence de réaction chimique, on a par définition des capacités thermiques isochores $C_{V}^{\alpha }$ et $C_{V}^{\beta }$ respectives des deux compartiments :
    $T^{\alpha }\,\mathrm {d} S^{\alpha }=C_{V}^{\alpha }\,\mathrm {d} T^{\alpha }$
    $T^{\beta }\,\mathrm {d} S^{\beta }=C_{V}^{\beta }\,\mathrm {d} T^{\beta }$
La condition de stabilité (csCV) donne :
    $C_{V}^{\alpha }\geq 0$
    $C_{V}^{\beta }\geq 0$
Puisque $\mathrm {d} S^{\alpha }<0$ et $\mathrm {d} S^{\beta }>0$, on a $\mathrm {d} T^{\alpha }<0$ et $\mathrm {d} T^{\beta }>0$. Dans le compartiment $\alpha$, où elle est la plus élevée, la température diminue. Elle augmente dans le compartiment $\beta$, où elle est la plus basse. Au cours du transfert thermique, l'écart $T^{\alpha }-T^{\beta }>0$ tend donc à diminuer.
Le transfert thermique cesse lorsque les températures des deux compartiments sont égales, soit $T^{\alpha }=T^{\beta }$, d'où $\mathrm {d} S=0$. L'entropie ne varie plus et le système atteint un équilibre. Si l'entropie $S^{\alpha }$ continuait à diminuer, alors $T^{\alpha }$ continuerait à diminuer et $T^{\beta }$ à augmenter, d'où $T^{\alpha }<T^{\beta }$ et $\mathrm {d} S<0$. De même, si l'entropie $S^{\alpha }$ augmentait à partir de cet équilibre, soit $\mathrm {d} S^{\alpha }>0$, alors $T^{\alpha }$ augmenterait et $T^{\beta }$ diminuerait, d'où $T^{\alpha }>T^{\beta }$ et $\mathrm {d} S<0$. L'entropie globale $S$, à partir de l'équilibre atteint, ne peut que diminuer, ce qui est contraire au deuxième principe de la thermodynamique : l'entropie a donc atteint un maximum et ne peut plus évoluer, cet équilibre est stable (voir la figure 1 avec $X_{i}=S^{\alpha }$).
Si l'on considère une transformation à entropie constante ($\mathrm {d} S=\mathrm {d} S^{\alpha }+\mathrm {d} S^{\beta }=0$), alors :
* À entropie constante $\mathrm {d} U=\left(T^{\alpha }-T^{\beta }\right)\,\mathrm {d} S^{\alpha }$
Un raisonnement similaire à celui tenu à énergie constante montre que toute évolution de $S^{\alpha }$ dans ces conditions à partir de l'équilibre ($T^{\alpha }=T^{\beta }$, d'où $\mathrm {d} U=0$) ne peut induire qu'une augmentation de l'énergie interne globale $U$, soit $\mathrm {d} U>0$, ce qui est contraire au premier principe de la thermodynamique : l'énergie interne a donc atteint un minimum, cet équilibre est stable (voir la figure 2 avec $X_{i}=S^{\alpha }$).
Pour résumer, si deux compartiments sont mis en contact de telle sorte qu'ils ne puissent échanger que de la chaleur, la chaleur passe du compartiment de plus haute température à celui de moindre température. La température la plus haute diminue et la plus basse augmente. L'équilibre thermique est atteint lorsque les deux températures s'égalisent. L'entropie du système a alors atteint un maximum et son énergie interne un minimum.

#### Le principe zéro de la thermodynamique

Le principe zéro de la thermodynamique énonce que$^,$ :
* Principe zéro de la thermodynamique Si deux systèmes A et B sont en équilibre thermique stable avec un troisième C, ils sont eux-mêmes en équilibre thermique stable.
Pour l'état initial, on a les deux équilibres, entre A et B d'une part, et B et C d'autre part :
    $T^{\boldsymbol {\mathsf {A}}}\,\mathrm {d} S^{\boldsymbol {\mathsf {A+C}}}=-\left(T^{\boldsymbol {\mathsf {C}}}-T^{\boldsymbol {\mathsf {A}}}\right)\,\mathrm {d} S^{\boldsymbol {\mathsf {C}}}=0$
    $T^{\boldsymbol {\mathsf {B}}}\,\mathrm {d} S^{\boldsymbol {\mathsf {B+C}}}=-\left(T^{\boldsymbol {\mathsf {C}}}-T^{\boldsymbol {\mathsf {B}}}\right)\,\mathrm {d} S^{\boldsymbol {\mathsf {C}}}=0$
On a alors $T^{\boldsymbol {\mathsf {A}}}=T^{\boldsymbol {\mathsf {B}}}=T^{\boldsymbol {\mathsf {C}}}$. Lorsque l'on met en contact A et B, les deux systèmes sont immédiatement à l'équilibre, puisque :
    $T^{\boldsymbol {\mathsf {B}}}\,\mathrm {d} S^{\boldsymbol {\mathsf {A+B}}}=-\left(T^{\boldsymbol {\mathsf {A}}}-T^{\boldsymbol {\mathsf {B}}}\right)\,\mathrm {d} S^{\boldsymbol {\mathsf {A}}}=0$
Comme vu dans la section Le transfert thermique, à partir de cet équilibre l'entropie $S^{\boldsymbol {\mathsf {A+B}}}$ du système A+B ne peut que diminuer, ce qu'interdit le deuxième principe de la thermodynamique. Le système A+B est donc d'emblée à l'équilibre thermique stable, sans qu'un transfert thermique soit nécessaire pour l'établir. Ainsi le principe zéro de la thermodynamique se déduit-il des premier et deuxième principes.
Ce principe permet de définir la notion de température comme une grandeur caractéristique de l'état d'équilibre thermodynamique. Tous les corps à la même température sont en équilibre thermique. En conséquence, le thermomètre est un instrument qui permet de vérifier cet équilibre.

#### La diffusion de la matière

Lorsqu'un milieu est hétérogène en composition, le potentiel chimique d'au moins une espèce chimique 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, ou migration.
On considère un système isolé divisé en deux compartiments $\alpha$ et $\beta$ séparés par une paroi. Cette paroi n'est perméable qu'à la seule espèce chimique $i$. Elle peut être mobile : le volume de chacun des deux compartiments peut varier, mais le volume global $V$ du système, en l'absence de travail des forces de pression avec l'extérieur, est constant. 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$). On suppose enfin que la pression $P$ et la température $T$ sont constantes et que le 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 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 }$
Le système étant isolé, en l'absence de réaction, on a les relations :
* $\mathrm {d} U=\mathrm {d} U^{\alpha }+\mathrm {d} U^{\beta }=0$ (premier principe, énergie globale $U=U^{\alpha }+U^{\beta }$) ;
* $\mathrm {d} S=\mathrm {d} S^{\alpha }+\mathrm {d} S^{\beta }\geq 0$ (deuxième principe, entropie globale $S=S^{\alpha }+S^{\beta }$) ;
* $\mathrm {d} V^{\alpha }+\mathrm {d} V^{\beta }=0$ (absence de travail des forces de pression) ;
* $\mathrm {d} n_{i}^{\alpha }+\mathrm {d} n_{i}^{\beta }=0$ (conservation de la matière).
Les quantités $n_{i}^{\alpha }$ et $n_{i}^{\beta }$ dans les deux compartiments varient donc en sens inverses : $\mathrm {d} n_{i}^{\beta }=-\mathrm {d} n_{i}^{\alpha }$, si l'une augmente, l'autre diminue. On obtient finalement :
* À énergie interne constante $T\,\mathrm {d} S=-\left(\mu _{i}^{\alpha }-\mu _{i}^{\beta }\right)\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$ : la quantité de $i$ diminue dans le compartiment $\alpha$ et augmente dans le compartiment $\beta$. 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. Cette relation ne permet pas de déterminer les variations des potentiels chimiques $\mu _{i}^{\alpha }$ et $\mu _{i}^{\beta }$ dans les deux compartiments.
Les deux compartiments évoluant à pression et température constantes, et en l'absence de réaction chimique, on a les variations des potentiels chimiques dans les deux compartiments :
    $\mathrm {d} \mu _{i}^{\alpha }=\left({\partial \mu _{i}^{\alpha } \over \partial n_{i}^{\alpha }}\right)_{P,T,n_{\neq i}^{\alpha }}\,\mathrm {d} n_{i}^{\alpha }$
    $\mathrm {d} \mu _{i}^{\beta }=\left({\partial \mu _{i}^{\beta } \over \partial n_{i}^{\beta }}\right)_{P,T,n_{\neq i}^{\beta }}\,\mathrm {d} n_{i}^{\beta }$
La condition de stabilité (csMu4) donne :
    $\left({\partial \mu _{i}^{\alpha } \over \partial n_{i}^{\alpha }}\right)_{P,T,n_{\neq i}^{\alpha }}\geq 0$
    $\left({\partial \mu _{i}^{\beta } \over \partial n_{i}^{\beta }}\right)_{P,T,n_{\neq i}^{\beta }}\geq 0$
Puisque $\mathrm {d} n_{i}^{\alpha }<0$ et $\mathrm {d} n_{i}^{\beta }>0$, on a $\mathrm {d} \mu _{i}^{\alpha }<0$ et $\mathrm {d} \mu _{i}^{\beta }>0$. Dans le compartiment $\alpha$ , où il est le plus élevé, le potentiel chimique diminue. Il augmente dans le compartiment $\beta$, où il est le plus bas. Au cours de la migration, l'écart $\mu _{i}^{\alpha }-\mu _{i}^{\beta }>0$ tend donc à diminuer.
La migration cesse lorsque les potentiels des deux compartiments sont égaux, soit $\mu _{i}^{\alpha }=\mu _{i}^{\beta }$, d'où $\mathrm {d} S=0$. L'entropie ne varie plus et le système atteint un équilibre. Si la quantité $n_{i}^{\alpha }$ continuait à diminuer, alors $\mu _{i}^{\alpha }$ continuerait à diminuer et $\mu _{i}^{\beta }$ à augmenter, d'où $\mu _{i}^{\alpha }<\mu _{i}^{\beta }$ et $\mathrm {d} S<0$. De même, si la quantité $n_{i}^{\alpha }$ augmentait à partir de cet équilibre, soit $\mathrm {d} n_{i}^{\alpha }>0$, alors $\mu _{i}^{\alpha }$ augmenterait et $\mu _{i}^{\beta }$ diminuerait, d'où $\mu _{i}^{\alpha }>\mu _{i}^{\beta }$ et $\mathrm {d} S<0$. L'entropie globale $S$, à partir de l'équilibre atteint, ne peut que diminuer, ce qui est contraire au deuxième principe de la thermodynamique : l'entropie a donc atteint un maximum et ne peut plus évoluer, cet équilibre est stable (voir la figure 1 avec $X_{i}=n_{i}^{\alpha }$).
Si l'on considère une transformation à entropie constante ($\mathrm {d} S=\mathrm {d} S^{\alpha }+\mathrm {d} S^{\beta }=0$), alors :
* À entropie constante $\mathrm {d} U=\left(\mu _{i}^{\alpha }-\mu _{i}^{\beta }\right)\mathrm {d} n_{i}^{\alpha }$
Un raisonnement similaire à celui tenu à énergie constante montre que toute évolution de $n_{i}^{\alpha }$ dans ces conditions à partir de l'équilibre ($\mu _{i}^{\alpha }=\mu _{i}^{\beta }$, d'où $\mathrm {d} U=0$) ne peut induire qu'une augmentation de l'énergie interne globale $U$, soit $\mathrm {d} U>0$, ce qui est contraire au premier principe de la thermodynamique : l'énergie interne a donc atteint un minimum, cet équilibre est stable (voir la figure 2 avec $X_{i}=n_{i}^{\alpha }$).
Pour résumer, si deux compartiments sont mis en contact de telle sorte qu'ils ne puissent échanger que de la matière, une espèce chimique migre du compartiment dans lequel son potentiel chimique est le plus fort à celui de moindre potentiel. Le potentiel le plus fort diminue et le plus faible augmente. L'équilibre diffusif est atteint lorsque les deux potentiels s'égalisent. L'entropie du système a alors atteint un maximum et son énergie interne un minimum.

### Équilibre de phases
#### Règle des phases

Soit un mélange quelconque composé de $N$ espèces chimiques réparties en $\phi$ phases. Les conditions d'équilibre stable d'un système hétérogène définies plus haut impliquent que :
* les $\phi$ pressions sont liées par les $\phi -1$ égalités :
    $P^{\alpha }=P^{\beta }=\cdots =P^{\phi }$
* les $\phi$ températures sont liées par les $\phi -1$ égalités :
    $T^{\alpha }=T^{\beta }=\cdots =T^{\phi }$
* pour chacun des $N$ constituants, les $\phi$ potentiels chimiques sont liés par les $\phi -1$ égalités :
    ${\begin{cases}\mu _{1}^{\alpha }=\mu _{1}^{\beta }=\cdots =\mu _{1}^{\phi }\\\mu _{2}^{\alpha }=\mu _{2}^{\beta }=\cdots =\mu _{2}^{\phi }\\\cdots \\\mu _{N}^{\alpha }=\mu _{N}^{\beta }=\cdots =\mu _{N}^{\phi }\end{cases}}$
Pour chaque phase, les diverses variables sont également liées par la relation de Gibbs-Duhem, soit $\phi$ équations supplémentaires :
    ${\begin{cases}V^{\alpha }\,\mathrm {d} P^{\alpha }-S^{\alpha }\,\mathrm {d} T^{\alpha }=\sum _{i=1}^{N}n_{i}^{\alpha }\,\mathrm {d} \mu _{i}^{\alpha }\\V^{\beta }\,\mathrm {d} P^{\beta }-S^{\beta }\,\mathrm {d} T^{\beta }=\sum _{i=1}^{N}n_{i}^{\beta }\,\mathrm {d} \mu _{i}^{\beta }\\\cdots \\V^{\phi }\,\mathrm {d} P^{\phi }-S^{\phi }\,\mathrm {d} T^{\phi }=\sum _{i=1}^{N}n_{i}^{\phi }\,\mathrm {d} \mu _{i}^{\phi }\end{cases}}$
On a donc un total de $X=\left(2+N\right)\phi$ variables et $Y=\left(2+N\right)\left(\phi -1\right)+\phi$ équations. La variance $v=X-Y$ indique le nombre de variables, c'est-à-dire le nombre de conditions opératoires (parmi la pression, la température et les compositions des phases), que l'opérateur peut fixer indépendamment l'une de l'autre en restant dans le domaine d'équilibre stable des phases, d'où la règle des phases :
$v=N+2-\phi$
Par exemple, pour un corps pur, soit $N=1$ :
* en présence d'une seule phase (par exemple un liquide seul) , soit $\phi =1$, on a $v=2$ : l'opérateur peut imposer la pression et la température indépendamment l'une de l'autre ;
* en présence de deux phases (par exemple un équilibre liquide-gaz), soit $\phi =2$, on a $v=1$ : l'opérateur peut imposer la pression mais la température d'équilibre est subie (et inversement) ;
* en présence de trois phases (par exemple le point triple de l'équilibre solide-liquide-gaz), soit $\phi =3$, on a $v=0$ : l'opérateur ne peut imposer aucune condition de l'équilibre.
Un mélange de $N$ constituants peut au plus se présenter sous $\phi =N+2$ phases à l'équilibre stable. Il ne peut y avoir au plus qu'une seule phase gaz, mais il peut y avoir plusieurs phases liquides et solides à l'équilibre.

#### L'entropie dans un changement d'état

On apporte, à pression $P$ et quantité de matière $n$ constantes, de la chaleur $Q$ à un corps pur. Son énergie interne $U$, en présence du travail des seules forces de pression, évolue selon :
    $\mathrm {d} U=-P\,\mathrm {d} V+\delta Q_{P}$
En introduisant l'enthalpie $H=U+PV$, on a, à pression $P$ constante :
    $\delta Q_{P}=\mathrm {d} U+P\,\mathrm {d} V=\mathrm {d} \left[U+PV\right]=\mathrm {d} H$
La transformation est réversible, on a par conséquent, avec l'entropie $S$ :
$\delta Q_{P}=T\,\mathrm {d} S=\mathrm {d} H$
Les différentielles de deuxième ordre des deux fonctions sont données par (d2S) et (d2H) à pression $P$ et quantité de matière $n$ constantes :
    $\mathrm {d} ^{2}S=\mathrm {d} {1 \over T}\,\mathrm {d} U+\mathrm {d} {P \over T}\,\mathrm {d} V=\mathrm {d} {1 \over T}\,\mathrm {d} \left[U+PV\right]=\mathrm {d} {1 \over T}\,\mathrm {d} H=-{1 \over T^{2}}\,\mathrm {d} T\,\mathrm {d} H$
    $\mathrm {d} ^{2}H=\mathrm {d} T\,\mathrm {d} S$
La transformation ayant lieu à pression constante, on a par définition de la capacité thermique isobare $C_{P}$ :
    $T\,\mathrm {d} S=C_{P}\,\mathrm {d} T$
La température thermodynamique $T$ est positive, et, selon (csCP), on a $C_{P}\geq 0$. En conséquence, avec (a) on obtient :
* (b) : $\mathrm {d} ^{2}S=-{1 \over C_{P}\,T^{2}}\,\left(\mathrm {d} H\right)^{2}\leq 0$ (c) : $\mathrm {d} ^{2}H={T \over C_{P}}\,\left(\mathrm {d} S\right)^{2}\geq 0$
Les équations (a), (b) et (c) montrent que, à pression et composition constantes, l'entropie $S$ est une fonction croissante concave de l'enthalpie $H$, qui est elle-même une fonction croissante convexe de l'entropie.
Sur la figure ci-contre, la courbe ABC bleue représente l'évolution, à pression et quantité de matière constantes, de l'entropie d'un corps pur dans une phase $a$, la courbe CDE orange celle du même corps pur dans une phase $b$. Les deux courbes sont concaves : si l'on joint deux points bleus, ou deux points orange, par un segment de droite, celui-ci est en-dessous de la courbe considérée. Cependant, la courbe globale ABCDE n'est pas dérivable au point C. De plus, elle est localement convexe entre B et D : si l'on joint un point de BC à un point de CD par un segment de droite, celui-ci est au-dessus de la courbe. Entre les points B et D, le corps pur ne peut être en équilibre stable s'il suit la courbe BCD. Pour conserver ses propriétés de continuité, dérivabilité et concavité, l'entropie passe par le segment de droite BD gris. Au point B le segment gris est tangent à la courbe bleue, au point D il est tangent à la courbe orange. La courbe ABDE est globalement dérivable et concave, contrairement à la courbe ABCDE. Le segment BD est la seule trajectoire qui rend l'entropie globalement dérivable et concave de A à E, tout autre point que B ou D, toute autre courbe qu'un segment de droite entre ces deux points engendre une courbe d'entropie localement non dérivable ou convexe.
Aux basses enthalpies, en amont du point B, le corps pur est exclusivement en phase $a$. En aval du point D, aux hautes enthalpies, le corps pur est exclusivement en phase $b$. Entre B et D, pour une enthalpie donnée, le segment gris représente une entropie supérieure aux entropies des phases $a$ et $b$ seules représentées respectivement par les courbes bleue et orange : cette entropie maximale dans ces conditions de pression et température représente le point de stabilité du système, mais ne représente pas une phase seule. Entre B et D, le corps pur subit un changement d'état et les deux phases $a$ et $b$, à enthalpie donnée, coexistent en un équilibre stable.
La relation $\mathrm {d} H=T\,\mathrm {d} S$ induit que sur le segment BD la température est constante : à la pression donnée, les deux phases du corps pur ne peuvent coexister qu'à une seule température, conformément à la règle des phases. Ainsi, les grandeurs $S^{\boldsymbol {\mathsf {B}}}$, $H^{\boldsymbol {\mathsf {B}}}$, $S^{\boldsymbol {\mathsf {D}}}$ et $H^{\boldsymbol {\mathsf {D}}}$ sont établies aux mêmes pression $P$, température $T_{a\to b}$ et quantité de matière $n$. La température de changement d'état est donnée par le rapport :
$T_{a\to b}={H^{\boldsymbol {\mathsf {D}}}-H^{\boldsymbol {\mathsf {B}}} \over S^{\boldsymbol {\mathsf {D}}}-S^{\boldsymbol {\mathsf {B}}}}$
L'écart $S^{\boldsymbol {\mathsf {D}}}-S^{\boldsymbol {\mathsf {B}}}$ rapporté à la quantité de matière totale $n$ mise en jeu est l'entropie de changement d'état du corps pur. L'écart $H^{\boldsymbol {\mathsf {D}}}-H^{\boldsymbol {\mathsf {B}}}$ rapporté à la quantité de matière totale $n$ mise en jeu est l'enthalpie de changement d'état du corps pur. L'enthalpie de changement d'état est la chaleur qu'il faut apporter à 1 mole de corps pur pour passer entièrement d'un état à l'autre :
* Entropie de changement d'état : $\Delta _{a\to b}S={S^{\boldsymbol {\mathsf {D}}}-S^{\boldsymbol {\mathsf {B}}} \over n}$ Enthalpie de changement d'état : $\Delta _{a\to b}H={H^{\boldsymbol {\mathsf {D}}}-H^{\boldsymbol {\mathsf {B}}} \over n}=T_{a\to b}\,\Delta _{a\to b}S$
L'équation du segment BD gris s'écrit :
    $S=\left(1-\tau ^{b}\right)\times S^{\boldsymbol {\mathsf {B}}}+\tau ^{b}\times S^{\boldsymbol {\mathsf {D}}}$
    $H=\left(1-\tau ^{b}\right)\times H^{\boldsymbol {\mathsf {B}}}+\tau ^{b}\times H^{\boldsymbol {\mathsf {D}}}$
Les points B et D représentant la même quantité $n$ de corps pur, on peut écrire, en introduisant les grandeurs molaires respectives :
    $S=\left(1-\tau ^{b}\right)n\times {\bar {S}}^{\boldsymbol {\mathsf {B}}}+\tau ^{b}\,n\times {\bar {S}}^{\boldsymbol {\mathsf {D}}}$
    $H=\left(1-\tau ^{b}\right)n\times {\bar {H}}^{\boldsymbol {\mathsf {B}}}+\tau ^{b}\,n\times {\bar {H}}^{\boldsymbol {\mathsf {D}}}$
Pour $\tau ^{b}=0$ le corps pur est au point B, entièrement en phase $a$ ; pour $\tau ^{b}=1$ il est au point D, entièrement en phase $b$. Pour $0<\tau ^{b}<1$ le corps pur est biphasique, entre B et D : une quantité $n^{a}=\left(1-\tau ^{b}\right)n$ de phase $a$ coexiste en un équilibre stable avec une quantité $n^{b}=\tau ^{b}\,n$ de phase $b$. Le paramètre $\tau ^{b}$ est le taux de changement d'état, ou fraction de phase $b$, c'est-à-dire la quantité de phase $b$ rapportée à la quantité de matière totale $n$ mise en jeu lors de la transformation :
$\tau ^{b}={n^{b} \over n^{a}+n^{b}}={n^{b} \over n}$
Autrement dit, si l'on dispose, à la pression $P$ et la température $T_{a\to b}$, d'une part d'une quantité $n$ de corps pur en phase $a$ et que l'on en prélève une fraction $1-\tau ^{b}$, et d'autre part d'une quantité $n$ de corps pur en phase $b$ et que l'on en prélève une fraction $\tau ^{b}$, si l'on met en contact ces deux fractions, alors celles-ci sont d'emblée en équilibre stable, sans nécessité d'une transformation pour atteindre cet équilibre (voir la section Concavité de l'entropie).
Les branches BC bleue et CD orange sont des branches de métastabilité des phases correspondantes. Si le corps pur évolue sur ces branches (y compris au-delà de C pour la phase $a$ et en deçà de C pour la phase $b$), la moindre perturbation peut provoquer un brusque changement d'état, amenant l'entropie sur la courbe ABDE stable, à enthalpie constante. L'état métastable du liquide non solidifié en dessous de sa température de fusion est appelé surfusion, il s'observe de façon naturelle lors des pluies verglaçantes (voir la vidéo présentée dans la section Équilibre stable). L'état métastable du liquide non vaporisé au-dessus de sa température d'ébullition est appelé surébullition ou retard à l'ébullition, il peut s'observer dans une tasse de café chauffé au microonde qui entre en ébullition brutale lorsque l'on y ajoute du sucre.
Le même raisonnement peut être tenu pour tout autre changement d'état du corps pur. Il est similaire pour les changements d'état des mélanges.

### Équilibre chimique
#### Conditions de l'équilibre chimique

Soit un système fermé siège d'une réaction chimique. Son énergie interne $U$ varie selon :
    $\mathrm {d} U=-P\,\mathrm {d} V+T\,\mathrm {d} S-{\mathcal {A}}\,\mathrm {d} \xi$
    $\mathrm {d} ^{2}U=-\mathrm {d} P\,\mathrm {d} V+\mathrm {d} T\,\mathrm {d} S-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \geq 0$
avec :
* ${\mathcal {A}}$ l'affinité chimique ; $-{\mathcal {A}}=U_{\xi }=F_{\xi }=H_{\xi }=G_{\xi }$ ;
* $\xi$ l'avancement de réaction.
Pour un système isolé, à $V$ et $U$ constants :
* la réaction évolue selon $-T\,\mathrm {d} S=-{\mathcal {A}}\,\mathrm {d} \xi \leq 0$ ;
* l'équilibre chimique est atteint lorsque :
** l'entropie $S$ a atteint un maximum : $\mathrm {d} S=0$ et $\mathrm {d} ^{2}S={1 \over T}\,\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \leq 0$ (voir la figure 1 avec $X_{i}=\xi$) ;
** l'énergie interne $U$ est à un minimum : $\mathrm {d} U=0$ et $\mathrm {d} ^{2}U=-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \geq 0$ (voir la figure 2 avec $X_{i}=\xi$).
Pour un système fermé contraint, selon les conditions opératoires :
* à volume $V$ et entropie $S$ constants :
** la réaction évolue selon $\mathrm {d} U=-{\mathcal {A}}\,\mathrm {d} \xi \leq 0$ ;
** l'équilibre chimique est atteint lorsque $\mathrm {d} U=0$ et $\mathrm {d} ^{2}U=-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \geq 0$ ; l'énergie interne $U$ est à un minimum ;
* à volume $V$ et température $T$ constants :
** la réaction évolue selon $\mathrm {d} F=-{\mathcal {A}}\,\mathrm {d} \xi \leq 0$ ;
** l'équilibre chimique est atteint lorsque $\mathrm {d} F=0$ et $\mathrm {d} ^{2}F=-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \geq 0$ ; l'énergie libre $F$ est à un minimum ;
* à pression $P$ et entropie $S$ constantes :
** la réaction évolue selon $\mathrm {d} H=-{\mathcal {A}}\,\mathrm {d} \xi \leq 0$ ;
** l'équilibre chimique est atteint lorsque $\mathrm {d} H=0$ et $\mathrm {d} ^{2}H=-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \geq 0$ ; l'enthalpie $H$ est à un minimum ;
* à pression $P$ et température $T$ constantes :
** la réaction évolue selon $\mathrm {d} G=-{\mathcal {A}}\,\mathrm {d} \xi \leq 0$ ;
** l'équilibre chimique est atteint lorsque $\mathrm {d} G=0$ et $\mathrm {d} ^{2}G=-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \geq 0$ ; l'enthalpie libre $G$ est à un minimum.
Pour un système isolé comme pour un système fermé, quelles que soient les conditions opératoires, la réaction chimique évolue donc selon :
* Condition générale d'évolution d'une réaction chimique ${\mathcal {A}}\,\mathrm {d} \xi \geq 0$
On note $X$ et $Y$ les deux conditions opératoires maintenues constantes. La seule variable de ${\mathcal {A}}$ laissée libre est $\xi$. On a par conséquent, à $X$ et $Y$ constantes :
    $\mathrm {d} {\mathcal {A}}=\left({\partial {\mathcal {A}} \over \partial \xi }\right)_{X,Y}\,\mathrm {d} \xi$
À l'équilibre, puisque dans tous les cas de figure $-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \geq 0$ :
    $\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi =\left({\partial {\mathcal {A}} \over \partial \xi }\right)_{X,Y}\,\left(\mathrm {d} \xi \right)^{2}\leq 0$
Ainsi, quelles que soient les conditions opératoires $X$ et $Y$ maintenues constantes, on a à l'équilibre$^,$ :
* Conditions générales d'équilibre chimique stable condition d'équilibre : ${\mathcal {A}}=0$ condition de stabilité : $\left({\partial {\mathcal {A}} \over \partial \xi }\right)_{X,Y}\leq 0$

#### Déplacement de l'équilibre

Soit un milieu réactionnel à l'équilibre. On modifie l'une des conditions du mélange réactionnel $X$ (pression, température, quantité de l'un des constituants, etc.) en considérant dans un premier temps que toutes les autres conditions $Y\neq X$ du milieu restent constantes. Soit $\Phi _{X}$ la variable conjuguée de $X$ :
$\Phi _{X}=\left({\partial \Phi \over \partial X}\right)_{Y\neq X}$
La fonction $\Phi =\Phi \left(Y\right)$ est le potentiel thermodynamique dont les $Y$ sont les variables naturelles, celles-ci incluant la variable $X$ modifiée et l'avancement de réaction $\xi$. En l'absence de réaction, soit à $\xi$ constant, le système évolue selon (voir la section Conclusions générales) :
* Sans la réaction Si $X$ est intensive, alors $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}\,\mathrm {d} X\leq 0$. Si $X$ est extensive, alors $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}\,\mathrm {d} X\geq 0$.
Le milieu réactionnel étant initialement à l'équilibre, son affinité chimique initiale est nulle : ${\mathcal {A}}_{0}=0$. La modification de $X$ induit la variation de l'affinité : $\mathrm {d} {\mathcal {A}}=\left(\partial {\mathcal {A}}/\partial X\right)_{Y\neq X}\,\mathrm {d} X$. L'affinité chimique est l'opposé de la dérivée partielle du potentiel thermodynamique $\Phi$ selon l'avancement de réaction $\xi$ : ${\mathcal {A}}=-\left(\partial \Phi /\partial \xi \right)_{Y\neq \xi }=-\Phi _{\xi }$. Le théorème de Schwarz donne :
    $\left({\partial {\mathcal {A}} \over \partial X}\right)_{Y\neq X}=-\left({\partial \over \partial X}\left({\partial \Phi \over \partial \xi }\right)_{Y\neq \xi }\right)_{Y\neq X}=-\left({\partial \over \partial \xi }\left({\partial \Phi \over \partial X}\right)_{Y\neq X}\right)_{Y\neq \xi }=-\left({\partial \Phi _{X} \over \partial \xi }\right)_{Y\neq \xi }$
On note la variation de $\Phi _{X}$ due à la réaction dans les conditions $Y\neq \xi$ constantes :
$\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }=\left({\partial \Phi _{X} \over \partial \xi }\right)_{Y\neq \xi }\,\mathrm {d} \xi$
On a par conséquent :
    $\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi =\left({\partial {\mathcal {A}} \over \partial X}\right)_{Y\neq X}\,\mathrm {d} X\,\mathrm {d} \xi =-\left({\partial \Phi _{X} \over \partial \xi }\right)_{Y\neq \xi }\,\mathrm {d} X\,\mathrm {d} \xi =-\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\,\mathrm {d} X$
Après modification de $X$, en l'absence de réaction, l'affinité chimique devient donc ${\mathcal {A}}={\mathcal {A}}\,_{0}+\mathrm {d} {\mathcal {A}}$. L'équilibre se déplace ensuite selon le sens dicté par le deuxième principe de la thermodynamique : $-{\mathcal {A}}\,\mathrm {d} \xi \leq 0$, d'où $-\mathrm {d} {\mathcal {A}}\,\mathrm {d} \xi \leq 0$. On obtient finalement :
* Par la réaction $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\,\mathrm {d} X\leq 0$
On en déduit que, par la réaction :
* une augmentation de $X$, soit $\mathrm {d} X>0$, induit une diminution de $\Phi _{X}$ : $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\leq 0$ ;
* une diminution de $X$, soit $\mathrm {d} X<0$, induit une augmentation de $\Phi _{X}$ : $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\geq 0$.
Par conséquent :
* La modification d'une variable quelconque $X$ déplace la réaction dans le sens impliquant une variation opposée de sa variable conjuguée $\Phi _{X}$.
On modifie la condition $X$ d'un milieu réactionnel initialement à l'équilibre. Le système évolue à la fois par cette modification externe et par le déplacement interne de l'équilibre chimique, jusqu'à atteindre un nouvel équilibre. Les conditions opératoires laissées constantes tout au long de ce processus sont les variables de $\Phi$ communes aux deux ensembles $\left[Y\neq X\right]$ et $\left[Y\neq \xi \right]$, soit $\left[Y\neq X\right]\cap \left[Y\neq \xi \right]$. Lorsque l'on modifie $X$, la variation totale de sa variable conjuguée $\Phi _{X}$ vaut par conséquent :
* Variation totale de $\Phi _{X}$ sans la réaction : $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}$, avec la réaction : $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}+\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }$.
On en déduit que$^,$$^,$ :
* si $X$ est intensive, sa variable conjuguée $\Phi _{X}$ est extensive :
** une augmentation de $X$, soit $\mathrm {d} X>0$, induit que $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}\leq 0$ et $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\leq 0$ ;
** une diminution de $X$, soit $\mathrm {d} X<0$, induit que $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}\geq 0$ et $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\geq 0$ ;
** la variation de $\Phi _{X}$ est plus importante avec la réaction que sans la réaction : ${\biggl |}\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}+\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }{\biggr |}\geq {\biggl |}\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}{\biggr |}$ ;
* si $X$ est extensive, sa variable conjuguée $\Phi _{X}$ est intensive :
** une augmentation de $X$, soit $\mathrm {d} X>0$, induit que $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}\geq 0$ et $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\leq 0$ ;
** une diminution de $X$, soit $\mathrm {d} X<0$, induit que $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}\leq 0$ et $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\geq 0$ ;
** la variation de $\Phi _{X}$ est moins importante avec la réaction que sans la réaction : ${\biggl |}\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}+\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }{\biggr |}\leq {\biggl |}\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}{\biggr |}$.
* La modification d'une variable intensive $X$ induit une variation de sa variable conjuguée extensive $\Phi _{X}$ plus grande avec la réaction que sans la réaction. Autrement dit, la modification de $X$ déplace l'équilibre dans le sens qui amplifie la modification de $\Phi _{X}$$^,$$^,$.
* La modification d'une variable extensive $X$ induit une variation de sa variable conjuguée intensive $\Phi _{X}$ plus petite avec la réaction que sans la réaction. Autrement dit, la modification de $X$ déplace l'équilibre dans le sens qui s'oppose à la modification de $\Phi _{X}$$^,$$^,$.

#### Applications courantes

Le tableau suivant résume les modifications de conditions opératoires les plus courantes d'un équilibre chimique.
Conditions d'évolution d'un système sans et avec la réaction. :
| Variable modifiée $X$ | Conditions opératoires constantes $\left[Y\neq X\right]\cap \left[Y\neq \xi \right]$ | Variable conjuguée $\Phi _{X}$ | Sans la réaction $\Phi _{XX}\,\left(\mathrm {d} X\right)^{2}=\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}\,\mathrm {d} X$ | Par la réaction ${\mathcal {A}}=-\Phi _{\xi }$ $\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }\,\mathrm {d} X\leq 0$ | Variation globale avec la réaction $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}+\mathrm {d} \Phi _{X}|_{\xi }^{Y\neq \xi }$ et sans la réaction $\mathrm {d} \Phi _{X}|_{X}^{Y\neq X}$ |
| $V$ | $T$ | (edeP) $-P=F_{V}$ | (csF1) $F_{VV}\,\left(\mathrm {d} V\right)^{2}=-\mathrm {d} P|_{V}^{T,n}\,\mathrm {d} V\geq 0$ | ${\mathcal {A}}=-F_{\xi }$ $-\mathrm {d} P|_{\xi }^{V,T}\,\mathrm {d} V\leq 0$ | $0\leq -\left(\mathrm {d} P|_{V}^{T,n}+\mathrm {d} P|_{\xi }^{V,T}\right)\,\mathrm {d} V\leq -\mathrm {d} P|_{V}^{T,n}\,\mathrm {d} V$ |
| $P$ | $T$ | (edeV) $V=G_{P}$ | (csG1) $G_{PP}\,\left(\mathrm {d} P\right)^{2}=\mathrm {d} V|_{P}^{T,n}\,\mathrm {d} P\leq 0$ | ${\mathcal {A}}=-G_{\xi }$ $\mathrm {d} V|_{\xi }^{P,T}\,\mathrm {d} P\leq 0$ | $\left(\mathrm {d} V|_{P}^{T,n}+\mathrm {d} V|_{\xi }^{P,T}\right)\,\mathrm {d} P\leq \mathrm {d} V|_{P}^{T,n}\,\mathrm {d} P\leq 0$ |
| $T$ | $V$ | (edeS) $-S=F_{T}$ | (csF2) $F_{TT}\,\left(\mathrm {d} T\right)^{2}=-\mathrm {d} S|_{T}^{V,n}\,\mathrm {d} T\leq 0$ | ${\mathcal {A}}=-F_{\xi }$ $-\mathrm {d} S|_{\xi }^{V,T}\,\mathrm {d} T\leq 0$ | $-\left(\mathrm {d} S|_{T}^{V,n}+\mathrm {d} S|_{\xi }^{V,T}\right)\,\mathrm {d} T\leq -\mathrm {d} S|_{T}^{V,n}\,\mathrm {d} T\leq 0$ |
| $T$ | $P$ | (edeS) $-S=G_{T}$ | (csG2) $G_{TT}\,\left(\mathrm {d} T\right)^{2}=-\mathrm {d} S|_{T}^{P,n}\,\mathrm {d} T\leq 0$ | ${\mathcal {A}}=-G_{\xi }$ $-\mathrm {d} S|_{\xi }^{P,T}\,\mathrm {d} T\leq 0$ | $-\left(\mathrm {d} S|_{T}^{P,n}+\mathrm {d} S|_{\xi }^{P,T}\right)\,\mathrm {d} T\leq -\mathrm {d} S|_{T}^{P,n}\,\mathrm {d} T\leq 0$ |
| $n_{i}$ | $V,T$ | (mu) $\mu _{i}=F_{n_{i}}$ | (csF3) $F_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}=\mathrm {d} \mu _{i}|_{n_{i}}^{V,T,n_{\neq i}}\,\mathrm {d} n_{i}\geq 0$ | ${\mathcal {A}}=-F_{\xi }$ $\mathrm {d} \mu _{i}|_{\xi }^{V,T}\,\mathrm {d} n_{i}\leq 0$ | $0\leq \left(\mathrm {d} \mu _{i}|_{n_{i}}^{V,T,n_{\neq i}}+\mathrm {d} \mu _{i}|_{\xi }^{V,T}\right)\,\mathrm {d} n_{i}\leq \mathrm {d} \mu _{i}|_{n_{i}}^{V,T,n_{\neq i}}\,\mathrm {d} n_{i}$ |
| $n_{i}$ | $P,T$ | (mu) $\mu _{i}=G_{n_{i}}$ | (csG3) $G_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}=\mathrm {d} \mu _{i}|_{n_{i}}^{P,T,n_{\neq i}}\,\mathrm {d} n_{i}\geq 0$ | ${\mathcal {A}}=-G_{\xi }$ $\mathrm {d} \mu _{i}|_{\xi }^{P,T}\,\mathrm {d} n_{i}\leq 0$ | $0\leq \left(\mathrm {d} \mu _{i}|_{n_{i}}^{P,T,n_{\neq i}}+\mathrm {d} \mu _{i}|_{\xi }^{P,T}\right)\,\mathrm {d} n_{i}\leq \mathrm {d} \mu _{i}|_{n_{i}}^{P,T,n_{\neq i}}\,\mathrm {d} n_{i}$ |
Exemple 1 - Modification du volume à température constante.
    Si l'on augmente le volume, soit $\mathrm {d} V>0$, à température constante, on a :
   * sans la réaction, une diminution de la pression : $\mathrm {d} P|_{V}^{T,n}\leq 0$ ;
   * par la réaction, une augmentation de la pression : $\mathrm {d} P|_{\xi }^{V,T}\geq 0$ ;
   * la variation globale de la pression : $\mathrm {d} P|_{V}^{T,n}\leq \mathrm {d} P|_{V}^{T,n}+\mathrm {d} P|_{\xi }^{V,T}\leq 0$.
    Si l'on augmente le volume à température constante, la pression diminue moins avec la réaction que sans la réaction. Inversement, si l'on diminue le volume, la pression augmente moins avec la réaction que sans la réaction. Lorsque l'on modifie le volume $V$ (grandeur extensive), la réaction s'oppose à la modification de la pression $P$ (grandeur intensive).
Exemple 2 - Modification de la pression à température constante.
    Si l'on augmente la pression à température constante, le volume diminue plus avec la réaction que sans la réaction. Inversement, si l'on diminue la pression, le volume augmente plus avec la réaction que sans la réaction. Lorsque l'on modifie la pression $P$ (grandeur intensive), la réaction amplifie la modification du volume $V$ (grandeur extensive).
Exemple 3 - Modification de la température.
    Si l'on augmente la température à pression constante ou volume constant, l'entropie du système augmente plus avec la réaction que sans la réaction. Inversement, si l'on diminue la température, l'entropie diminue plus avec la réaction que sans la réaction. Lorsque l'on modifie la température $T$ (grandeur intensive), la réaction amplifie la modification de l'entropie $S$ du milieu réactionnel (grandeur extensive), c'est-à-dire l'échange de chaleur du milieu réactionnel avec l'extérieur.
    Une réaction qui augmente l'entropie du milieu réactionnel, soit $\mathrm {d} S|_{\xi }\geq 0$, absorbe de la chaleur, elle est endothermique. Une réaction qui diminue l'entropie, soit $\mathrm {d} S|_{\xi }\leq 0$, dégage de la chaleur, elle est exothermique. Lorsque l'on chauffe un milieu réactionnel pour atteindre une certaine température, il faut donc lui fournir plus de chaleur avec la réaction que sans la réaction, car la réaction absorbe une partie de la chaleur fournie. Si l'on fournit une certaine quantité de chaleur au milieu réactionnel, la température obtenue est moins élevée avec la réaction que sans la réaction. En augmentant la température, l'équilibre se déplace dans le sens endothermique : il amplifie l'absorption de chaleur et semble ainsi s'opposer à la modification de la température. Il s'agit de la loi de van 't Hoff.
Exemple 4 - Modification de la composition.
    Si l'on augmente la quantité d'un constituant (réactif, produit ou inerte), que ce soit à volume et température constants, ou à pression et température constantes, le potentiel chimique de ce constituant augmente moins avec la réaction que sans la réaction. Inversement, si l'on diminue la quantité d'un constituant, le potentiel chimique de ce constituant diminue moins avec la réaction que sans la réaction. Lorsque l'on modifie la quantité $n_{i}$ d'un constituant $i$ quelconque du milieu réactionnel (grandeur extensive), la réaction s'oppose à la modification du potentiel chimique $\mu _{i}$ de cette espèce (grandeur intensive).
    Ceci n'implique pas que la réaction se déplace dans le sens qui consomme le constituant ajouté et diminue sa quantité. Il existe en effet des cas où la réaction se déplace dans le sens qui produit le constituant ajouté et augmente donc sa quantité, son potentiel chimique étant néanmoins diminué.

### Thermodynamique des surfaces
#### Surface minimale

On considère un système fermé séparé de l'extérieur par une interface formant une surface d'aire $A$ (par exemple une bulle de gaz dans un liquide). Lorsque l'on déforme cette interface, le travail de la force exercée sur sa surface vaut $\gamma \,\mathrm {d} A$. La tension superficielle $\gamma$ est une grandeur positive : on doit fournir un travail au système pour agrandir sa surface. Ce travail est généralement négligeable et négligé devant le travail des forces de pression. À composition constante, l'énergie interne $U$ du système varie selon :
$\mathrm {d} U=-P\,\mathrm {d} V+\gamma \,\mathrm {d} A+T\,\mathrm {d} S$
On a donc la condition d'évolution du système à composition constante :
* Condition d'évolution $\mathrm {d} \Phi =\mathrm {d} U+P\,\mathrm {d} V-T\,\mathrm {d} S=\gamma \,\mathrm {d} A=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]\leq 0$
La tension superficielle peut ainsi être définie indifféremment par :
* Tension superficielle $\gamma =\left({\partial U \over \partial A}\right)_{V,S,n}=\left({\partial F \over \partial A}\right)_{V,T,n}=\left({\partial H \over \partial A}\right)_{P,S,n}=\left({\partial G \over \partial A}\right)_{P,T,n}$
Quelles que soient les conditions opératoires considérées, le système évolue donc selon $\mathrm {d} A\leq 0$ : la surface de l'interface ne peut que décroître. À l'équilibre stable $\Phi$ atteint un minimum. Toute variation d'un paramètre de l'équilibre ne peut alors que faire croître $\Phi$, soit $\mathrm {d} \Phi \geq 0$ ; par conséquent, à partir de l'équilibre $\mathrm {d} A\geq 0$ : la surface de l'interface ne peut que croître. Quelles que soient les conditions opératoires, on a donc :
* À l'équilibre thermodynamique stable, la surface de l'interface entre le système et l'extérieur est minimale.
Une bulle de gaz dans un liquide est donc sphérique (si sa vitesse de déplacement est faible, sinon elle se déforme), de même qu'une gouttelette de liquide dans un aérosol. Un liquide dans un verre forme une interface plane avec l'air. Si l'on déforme ces interfaces, leur surface augmente. Puis, lorsque l'on supprime la perturbation, la surface tend à reprendre sa forme initiale, minimale.

#### La coalescence

Soit un système composé d'une phase dispersée dans une autre phase (l'extérieur). Ce système comporte un ensemble de $K$ sous-systèmes notés $k$, de tailles variées et à l'équilibre entre eux. Chaque sous-système $k$ est caractérisé par un volume $V^{k}$, une surface $A^{k}$, une entropie $S^{k}$ et les quantités $n_{i}^{k}$ des $N$ constituants $i$. Ces sous-systèmes sont aux mêmes pression $P$, température $T$ et composition (fractions molaires), d'où pour chaque constituant $i$ le même potentiel chimique $\mu _{i}$ dans chacun des sous-systèmes. Puisque l'énergie interne est extensive, la variation d'énergie interne globale $U$ du système est la somme des variations d'énergie interne des sous-systèmes :
    $\mathrm {d} U=\sum _{k=1}^{K}\left[-P\,\mathrm {d} V^{k}+\gamma \,\mathrm {d} A^{k}+T\,\mathrm {d} S^{k}+\sum _{i=1}^{N}\mu _{i}\,\mathrm {d} n_{i}^{k}\right]$
On pose le volume total du système $V=\sum _{k=1}^{K}V^{k}$, sa surface totale $A=\sum _{k=1}^{K}A^{k}$, son entropie $S=\sum _{k=1}^{K}S^{k}$ et l'enthalpie libre $G=U+PV-TS$. La quantité globale de chacun des constituants $i$ est invariante :
    $\sum _{k=1}^{K}\mathrm {d} n_{i}^{k}=0$
On a finalement, à pression, température et composition constantes, la condition d'évolution de la phase dispersée :
* Condition d'évolution $\mathrm {d} \Phi =\mathrm {d} G=\gamma \,\mathrm {d} A=-T\,\mathrm {d} \left[S+S_{\text{ext}}\right]\leq 0$
Note : le même raisonnement mené non pas à pression constante mais à volume constant, soit $\sum _{k=1}^{K}\mathrm {d} V^{k}=0$, conduit à l'étude de l'énergie libre $F=U-TS$ à la place de $G$. Les conclusions suivantes restent identiques.
La surface $A$ globale du système ne peut que décroître. En conséquence, lorsque deux sous-systèmes initialement séparés, de surfaces $A^{1}$ et $A^{2}$ quelconques, fusionnent pour ne plus en former qu'un, la surface $A^{1+2}$ du système résultant est, à l'équilibre, inférieure à la somme des deux surfaces initiales, soit $A^{1+2}\leq A^{1}+A^{2}$. Deux gouttes d'eau sphériques qui fusionnent donnent une nouvelle goutte sphérique ; le volume est conservé, mais la surface diminue.
Ceci explique la tendance des milieux dispersés (émulsions, aérosols, mousses, etc.) à la coalescence : ils sont instables par nature. Les sous-systèmes, bulles ou gouttes, fusionnent spontanément, ce phénomène accroît l'entropie et diminue l'énergie du milieu en diminuant sa surface. Dans un empilement de bulles de savon (une mousse), les parois internes tendent à se rompre afin de ne former qu'une seule bulle. Dans une vinaigrette (une émulsion), les gouttelettes d'huile dispersées dans la phase aqueuse tendent spontanément à fusionner jusqu'à ne plus former qu'une seule phase continue : l'huile et la phase aqueuse se séparent pour atteindre l'équilibre.

## Notations et outils mathématiques
### Notations
#### Alphabet latin

* ${\mathcal {A}}$ l'affinité chimique ;
* $A$ la surface ;
* $F$ l'énergie libre ;
* $G$ l'enthalpie libre ;
* $H$ l'enthalpie ;
* $N$ le nombre d'espèces chimiques différentes dans le système ;
* $n_{i}$ la quantité de l'espèce chimique $i$ ;
* $P$ la pression ;
* $Q$ la chaleur ;
* $S$ l'entropie du système thermodynamique étudié ;
* $S_{\text{ext}}$ l'entropie de l'extérieur ;
* $T$ la température thermodynamique ;
* $U$ l'énergie interne du système thermodynamique étudié ;
* $U_{\text{ext}}$ l'énergie interne de l'extérieur ;
* $V$ le volume du système thermodynamique étudié ;
* $V_{\text{ext}}$ le volume de l'extérieur ;
* $W$ le travail ;
* ${\bar {X}}$ une grandeur molaire ;
* ${\bar {X}}_{i}$ une grandeur molaire partielle liée à l'espèce chimique $i$ ;
* $X$, $Y$, $Z$ des variables quelconques.

#### Alphabet grec

* $\gamma$ la tension superficielle ;
* $\xi$ l'avancement de réaction ;
* $\mu _{i}$ le potentiel chimique de l'espèce chimique $i$ ;
* $\Phi$, $\Psi$ des fonctions quelconques ;
* $\Phi _{X_{i}}$ la dérivée partielle de $\Phi$ par rapport à la variable $X_{i}$ ;
* $\Phi _{X_{i}X_{j}}$ la dérivée partielle seconde de $\Phi$ par rapport aux variables $X_{i}$ et $X_{j}$.

### Différentielles et dérivées des fonctions thermodynamiques
#### Série de Taylor

Soit une fonction $\Phi$ quelconque, dont les $n$ variables sont $X_{1},X_{2},\cdots ,X_{n}$ :
$\Phi =\Phi \!\left(X_{1},X_{2},\cdots ,X_{n}\right)$
On suppose cette fonction au moins $K$ fois dérivable par rapport à chacune de ses variables. Selon le théorème de Taylor, on peut poser en tout point $X^{\circ }=(X_{1}^{\circ },X_{2}^{\circ },\cdots ,X_{n}^{\circ })$, pour tout ensemble de variations élémentaires $\mathrm {d} X=(\mathrm {d} X_{1},\mathrm {d} X_{2},\cdots ,\mathrm {d} X_{n})$ quelconques :
$\Phi (X^{\circ }+\mathrm {d} X)\approx \Phi (X^{\circ })+\sum _{k=1}^{K}{1 \over k!}\,\mathrm {d} ^{k}\Phi (X^{\circ },\mathrm {d} X)$
Dans cette expression $\mathrm {d} ^{k}\Phi (X^{\circ },\mathrm {d} X)$ est la différentielle d'ordre $k$ de $\Phi$ au point $X^{\circ }$ avec les variations élémentaires $\mathrm {d} X$ :
* Différentielle d'ordre $k$ ${\begin{aligned}\mathrm {d} ^{k}\Phi (X^{\circ },\mathrm {d} X)&=\sum _{i_{1}=1}^{n}\sum _{i_{2}=1}^{n}\cdots \sum _{i_{k}=1}^{n}\left({\partial ^{k}\Phi \over \partial X_{i_{1}}\partial X_{i_{2}}\cdots \partial X_{i_{k}}}\right)\!\left(X^{\circ }\right)\,\mathrm {d} X_{i_{1}}\,\mathrm {d} X_{i_{2}}\,\cdots \,\mathrm {d} X_{i_{k}}\\&=\left(\left[\sum _{i=1}^{n}\left({\partial \,\cdot \over \partial X_{i}}\right)\,\mathrm {d} X_{i}\right]^{k}\Phi \right)\!\left(X^{\circ }\right)\end{aligned}}$
Le terme $\left({\partial ^{k}\Phi \over \partial X_{i_{1}}\partial X_{i_{2}}\cdots \partial X_{i_{k}}}\right)\!\left(X^{\circ }\right)$ est la dérivée partielle d'ordre $k$ de la fonction $\Phi$ par rapport à ses variables $X_{i_{1}},X_{i_{2}},\cdots ,X_{i_{k}}$ calculée au point $X^{\circ }$.
On suppose que toutes les fonctions thermodynamiques, entropie $S$ et potentiels thermodynamiques $U$, $F$, $H$ et $G$, sont au moins deux fois dérivables par rapport à l'ensemble de leurs variables.

#### Différentielle du premier ordre

La différentielle du premier ordre d'une fonction $\Phi$ quelconque s'écrit$^,$ :
    ${\begin{aligned}\mathrm {d} \Phi &=\left({\partial \Phi \over \partial X_{1}}\right)_{X_{\neq 1}}\,\mathrm {d} X_{1}+\left({\partial \Phi \over \partial X_{2}}\right)_{X_{\neq 2}}\,\mathrm {d} X_{2}+\cdots +\left({\partial \Phi \over \partial X_{n}}\right)_{X_{\neq n}}\,\mathrm {d} X_{n}\\&=\sum _{i=1}^{n}\left({\partial \Phi \over \partial X_{i}}\right)_{X_{\neq i}}\,\mathrm {d} X_{i}\end{aligned}}$
On note toute dérivée partielle du premier ordre :
    $\Phi _{X_{i}}=\left({\partial \Phi \over \partial X_{i}}\right)_{X_{\neq i}}$
on obtient :
* Différentielle du premier ordre $\mathrm {d} \Phi =\sum _{i=1}^{n}\Phi _{X_{i}}\,\mathrm {d} X_{i}$
Les différentielles du premier ordre des fonctions thermodynamiques s'écrivent, dans leurs variables naturelles :
* entropie $S$ :
       $\mathrm {d} S=S_{V}\,\mathrm {d} V+S_{U}\,\mathrm {d} U+\sum _{i=1}^{N}S_{n_{i}}\,\mathrm {d} n_{i}={P \over T}\,\mathrm {d} V+{1 \over T}\,\mathrm {d} U-\sum _{i=1}^{N}{\mu _{i} \over T}\,\mathrm {d} n_{i}$
* énergie interne $U$ :
       $\mathrm {d} U=U_{V}\,\mathrm {d} V+U_{S}\,\mathrm {d} S+\sum _{i=1}^{N}U_{n_{i}}\,\mathrm {d} n_{i}=-P\,\mathrm {d} V+T\,\mathrm {d} S+\sum _{i=1}^{N}\mu _{i}\,\mathrm {d} n_{i}$
* énergie libre $F$ :
       $\mathrm {d} F=F_{V}\,\mathrm {d} V+F_{T}\,\mathrm {d} T+\sum _{i=1}^{N}F_{n_{i}}\,\mathrm {d} n_{i}=-P\,\mathrm {d} V-S\,\mathrm {d} T+\sum _{i=1}^{N}\mu _{i}\,\mathrm {d} n_{i}$
* enthalpie $H$ :
       $\mathrm {d} H=H_{P}\,\mathrm {d} P+H_{S}\,\mathrm {d} S+\sum _{i=1}^{N}H_{n_{i}}\,\mathrm {d} n_{i}=V\,\mathrm {d} P+T\,\mathrm {d} S+\sum _{i=1}^{N}\mu _{i}\,\mathrm {d} n_{i}$
* enthalpie libre $G$ :
       $\mathrm {d} G=G_{P}\,\mathrm {d} P+G_{T}\,\mathrm {d} T+\sum _{i=1}^{N}G_{n_{i}}\,\mathrm {d} n_{i}=V\,\mathrm {d} P-S\,\mathrm {d} T+\sum _{i=1}^{N}\mu _{i}\,\mathrm {d} n_{i}$

##### Notes

* La pression $P$, la température $T$ et les potentiels chimiques $\mu _{i}$ sont des grandeurs intensives. Les fonctions $S$, $U$, $F$, $H$ et $G$, le volume $V$ et les quantités $n_{i}$ sont des grandeurs extensives.
* Les variables de l'entropie $S\!\left(V,U,n_{1},\cdots ,n_{N}\right)$ et celles de l'énergie interne $U\!\left(V,S,n_{1},\cdots ,n_{N}\right)$ sont toutes extensives.
* Dans le cas de l'entropie $S$, les grandeurs $S_{X_{i}}$ sont intensives et le produit $S_{X_{i}}X_{i}$ est homogène à une entropie : $S_{X_{i}}$ et $X_{i}$ sont « conjuguées par rapport à l'entropie ». Dans le cas des potentiels thermodynamiques $U$, $F$, $H$ et $G$, les grandeurs $\Phi _{X_{i}}$ et $X_{i}$ sont « conjuguées par rapport à l'énergie » : l'une des deux est extensive et l'autre intensive, et leur produit $\Phi _{X_{i}}X_{i}$ est homogène à une énergie.

#### Transformation de Legendre

Soit une fonction $\Phi$ quelconque, dont les variables naturelles sont $Y$ et $X_{1},\cdots ,X_{n}$ :
$\Phi =\Phi \!\left(Y,X_{1},\cdots ,X_{n}\right)$
avec :
    $\mathrm {d} \Phi =\Phi _{Y}\,\mathrm {d} Y+\sum _{i=1}^{n}\Phi _{X_{i}}\,\mathrm {d} X_{i}$
On pose :
$Z=\Phi _{Y}$
Soit $\Psi$ la transformée de Legendre de $\Phi$ telle que$^,$ :
$\Psi =\Phi -ZY$
On a par conséquent :
    ${\begin{aligned}\mathrm {d} \Psi &=\mathrm {d} \Phi -Z\,\mathrm {d} Y-Y\,\mathrm {d} Z=\mathrm {d} \Phi -\Phi _{Y}\,\mathrm {d} Y-Y\,\mathrm {d} Z\\&=-Y\,\mathrm {d} Z+\sum _{i=1}^{n}\Phi _{X_{i}}\,\mathrm {d} X_{i}\end{aligned}}$
La nouvelle fonction $\Psi$ dépend de $Z$ et des $X_{1},\cdots ,X_{n}$ :
$\Psi =\Psi \!\left(Z,X_{1},\cdots ,X_{n}\right)$
Avec :
    $\mathrm {d} \Psi =\Psi _{Z}\,\mathrm {d} Z+\sum _{i=1}^{n}\Psi _{X_{i}}\,\mathrm {d} X_{i}$
on obtient :
* $\Psi _{Z}=-Y$ $\Psi _{X_{i}}=\Phi _{X_{i}}$
On pose $\Pi =-P$. On a $\mathrm {d} \Pi =-\mathrm {d} P$ et pour toute fonction $\Phi _{\Pi }=-\Phi _{P}$ et $\Phi _{\Pi }\,\mathrm {d} \Pi =\Phi _{P}\,\mathrm {d} P$. Les tableaux suivants résument les relations entre les principaux potentiels thermodynamiques$^,$.
Tableau 1 - Transformation de Legendre. :
| Potentiel $\Phi \!\left(Y,X_{1},\cdots ,X_{n}\right)$ | Potentiel $\Psi \!\left(Z,X_{1},\cdots ,X_{n}\right)=\Phi -ZY$ | Variable $Y=-\Psi _{Z}$ | Variable $Z=\Phi _{Y}$ | Variables communes $X_{1},\cdots ,X_{n}$ |
| Énergie interne $U$ | Énergie libre $F=U-TS$ | $S=-F_{T}$ | $T=U_{S}$ | $V,n_{1},\cdots ,n_{N}$ |
| Énergie interne $U$ | Enthalpie $H=U-\Pi V$ | $V=-H_{\Pi }$ | $\Pi =U_{V}$ | $S,n_{1},\cdots ,n_{N}$ |
| Énergie libre $F$ | Enthalpie libre $G=F-\Pi V$ | $V=-G_{\Pi }$ | $\Pi =F_{V}$ | $T,n_{1},\cdots ,n_{N}$ |
| Enthalpie $H$ | Enthalpie libre $G=H-TS$ | $S=-G_{T}$ | $T=H_{S}$ | $\Pi ,n_{1},\cdots ,n_{N}$ |

##### Notes

* Les transformées de Legendre de l'entropie sont les fonctions de Massieu.
* L'enthalpie libre est le résultat d'une double transformation de Legendre : $G=F-\Pi V=H-TS=U-\Pi V-TS$.
* Les variables $Y$ et $Z$ sont conjuguées : l'une est extensive, l'autre est intensive, et leur produit $ZY$ est homogène à une énergie.
* Les variables naturelles de l'énergie interne $U$ sont toutes extensives. La transformation de Legendre remplace une variable extensive $Y$ de $U$ par une variable intensive $Z$ d'un autre potentiel : remplacement du volume $V$ de $U$ par la pression $P$ de $H$ et $G$, remplacement de l'entropie $S$ de $U$ par la température $T$ de $F$ et $G$.
Tableau 2 - Relations entre dérivées premières. :
| Potentiel $\Psi =\Phi -ZY$ | Relation $\Psi _{X_{i}}=\Phi _{X_{i}}$ | Relation $\Psi _{X_{i}}=\Phi _{X_{i}}$ |
| Potentiel $\Psi =\Phi -ZY$ | Équations d'état | Potentiel chimique |
| $F=U-TS$ | (edeP) $F_{V}=U_{V}=\Pi$ | $F_{n_{i}}=U_{n_{i}}=\mu _{i}$ |
| $H=U-\Pi V$ | (edeT) $H_{S}=U_{S}=T$ | $H_{n_{i}}=U_{n_{i}}=\mu _{i}$ |
| $G=F-\Pi V$ | (edeS) $G_{T}=F_{T}=-S$ | $F_{n_{i}}=G_{n_{i}}=\mu _{i}$ |
| $G=H-TS$ | (edeV) $G_{\Pi }=H_{\Pi }=-V$ | $G_{n_{i}}=H_{n_{i}}=\mu _{i}$ |

##### Notes

* Les équations d'état sont des relations définissant la pression (edeP), la température (edeT), l'entropie (edeS) et le volume (edeV) en tant que dérivées partielles des potentiels thermodynamiques.
* Le potentiel chimique de l'espèce chimique $i$ est défini indifféremment par :
$\mu _{i}=U_{n_{i}}=F_{n_{i}}=H_{n_{i}}=G_{n_{i}}$

#### Différentielle du deuxième ordre

La différentielle du deuxième ordre d'une fonction $\Phi$ quelconque s'écrit$^,$ :
    ${\begin{aligned}\mathrm {d} ^{2}\Phi &=\left({\partial \over \partial X_{1}}\left({\partial \Phi \over \partial X_{1}}\right)_{X_{\neq 1}}\right)_{X_{\neq 1}}\,\mathrm {d} X_{1}\,\mathrm {d} X_{1}+\left({\partial \over \partial X_{2}}\left({\partial \Phi \over \partial X_{1}}\right)_{X_{\neq 1}}\right)_{X_{\neq 2}}\,\mathrm {d} X_{2}\,\mathrm {d} X_{1}+\cdots +\left({\partial \over \partial X_{n}}\left({\partial \Phi \over \partial X_{1}}\right)_{X_{\neq 1}}\right)_{X_{\neq n}}\,\mathrm {d} X_{n}\,\mathrm {d} X_{1}\\&+\left({\partial \over \partial X_{1}}\left({\partial \Phi \over \partial X_{2}}\right)_{X_{\neq 2}}\right)_{X_{\neq 1}}\,\mathrm {d} X_{1}\,\mathrm {d} X_{2}+\left({\partial \over \partial X_{2}}\left({\partial \Phi \over \partial X_{2}}\right)_{X_{\neq 2}}\right)_{X_{\neq 2}}\,\mathrm {d} X_{2}\,\mathrm {d} X_{2}+\cdots +\left({\partial \over \partial X_{n}}\left({\partial \Phi \over \partial X_{2}}\right)_{X_{\neq 2}}\right)_{X_{\neq n}}\,\mathrm {d} X_{n}\,\mathrm {d} X_{2}\\&+\cdots \\&+\left({\partial \over \partial X_{1}}\left({\partial \Phi \over \partial X_{n}}\right)_{X_{\neq n}}\right)_{X_{\neq 1}}\,\mathrm {d} X_{1}\,\mathrm {d} X_{n}+\left({\partial \over \partial X_{2}}\left({\partial \Phi \over \partial X_{n}}\right)_{X_{\neq n}}\right)_{X_{\neq 2}}\,\mathrm {d} X_{2}\,\mathrm {d} X_{n}+\cdots +\left({\partial \over \partial X_{n}}\left({\partial \Phi \over \partial X_{n}}\right)_{X_{\neq n}}\right)_{X_{\neq n}}\,\mathrm {d} X_{n}\,\mathrm {d} X_{n}\\&=\sum _{i=1}^{n}\sum _{j=1}^{n}\left({\partial \over \partial X_{i}}\left({\partial \Phi \over \partial X_{j}}\right)_{X_{\neq j}}\right)_{X_{\neq i}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}\end{aligned}}$
On note toute dérivée partielle du deuxième ordre :
    $\Phi _{X_{i}X_{j}}=\left({\partial ^{2}\Phi \over \partial X_{i}\partial X_{j}}\right)_{X_{\neq i,j}}=\left({\partial \over \partial X_{i}}\left({\partial \Phi \over \partial X_{j}}\right)_{X_{\neq j}}\right)_{X_{\neq i}}$
d'où :
    $\mathrm {d} ^{2}\Phi =\sum _{i=1}^{n}\sum _{j=1}^{n}\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}$
Le théorème de Schwarz permettant d'écrire, pour toutes variables $X_{i}$ et $X_{j}$$^,$ :
* Théorème de Schwarz $\Phi _{X_{i}X_{j}}=\Phi _{X_{j}X_{i}}$
on obtient :
* Différentielle du deuxième ordre $\mathrm {d} ^{2}\Phi =\sum _{i=1}^{n}\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}+2\,\sum _{i=1}^{n}\sum _{j>i}^{n}\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}$
On a par ailleurs la différentielle pour toute variable $X_{i}$, en fonction des variables de $\Phi$ :
    ${\begin{aligned}\mathrm {d} \Phi _{X_{i}}&=\left({\partial \Phi _{X_{i}} \over \partial X_{1}}\right)_{X_{\neq 1}}\,\mathrm {d} X_{1}+\left({\partial \Phi _{X_{i}} \over \partial X_{2}}\right)_{X_{\neq 2}}\,\mathrm {d} X_{2}+\cdots +\left({\partial \Phi _{X_{i}} \over \partial X_{n}}\right)_{X_{\neq n}}\,\mathrm {d} X_{n}\\&=\sum _{j=1}^{n}\Phi _{X_{j}X_{i}}\,\mathrm {d} X_{j}=\sum _{j=1}^{n}\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{j}\end{aligned}}$
on obtient :
* Différentielle du deuxième ordre $\mathrm {d} ^{2}\Phi =\sum _{i=1}^{n}\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}$
Les différentielles du deuxième ordre des fonctions thermodynamiques s'écrivent, dans leurs variables naturelles :
* entropie $S$ (d2S) :
       ${\begin{aligned}\mathrm {d} ^{2}S&=S_{VV}\,\left(\mathrm {d} V\right)^{2}+S_{UU}\,\left(\mathrm {d} U\right)^{2}+\sum _{i=1}^{N}S_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\,S_{VU}\,\mathrm {d} V\,\mathrm {d} U+2\sum _{i=1}^{N}S_{Vn_{i}}\,\mathrm {d} V\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}S_{Un_{i}}\,\mathrm {d} U\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}\sum _{j>i}^{N}S_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\\&=\mathrm {d} {P \over T}\,\mathrm {d} V+\mathrm {d} {1 \over T}\,\mathrm {d} U-\sum _{i=1}^{N}\mathrm {d} {\mu _{i} \over T}\,\mathrm {d} n_{i}\end{aligned}}$
* énergie interne $U$ (d2U) :
       ${\begin{aligned}\mathrm {d} ^{2}U&=U_{VV}\,\left(\mathrm {d} V\right)^{2}+U_{SS}\,\left(\mathrm {d} S\right)^{2}+\sum _{i=1}^{N}U_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\,U_{VS}\,\mathrm {d} V\,\mathrm {d} S+2\sum _{i=1}^{N}U_{Vn_{i}}\,\mathrm {d} V\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}U_{Sn_{i}}\,\mathrm {d} S\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}\sum _{j>i}^{N}U_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\\&=-\mathrm {d} P\,\mathrm {d} V+\mathrm {d} T\,\mathrm {d} S+\sum _{i=1}^{N}\mathrm {d} \mu _{i}\,\mathrm {d} n_{i}\end{aligned}}$
* énergie libre $F$ (d2F) :
       ${\begin{aligned}\mathrm {d} ^{2}F&=F_{VV}\,\left(\mathrm {d} V\right)^{2}+F_{TT}\,\left(\mathrm {d} T\right)^{2}+\sum _{i=1}^{N}F_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\,F_{VT}\,\mathrm {d} V\,\mathrm {d} T+2\sum _{i=1}^{N}F_{Vn_{i}}\,\mathrm {d} V\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}F_{Tn_{i}}\,\mathrm {d} T\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}\sum _{j>i}^{N}F_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\\&=-\mathrm {d} P\,\mathrm {d} V-\mathrm {d} S\,\mathrm {d} T+\sum _{i=1}^{N}\mathrm {d} \mu _{i}\,\mathrm {d} n_{i}\end{aligned}}$
* enthalpie $H$ (d2H) :
       ${\begin{aligned}\mathrm {d} ^{2}H&=H_{PP}\,\left(\mathrm {d} P\right)^{2}+H_{SS}\,\left(\mathrm {d} S\right)^{2}+\sum _{i=1}^{N}H_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\,H_{PS}\,\mathrm {d} P\,\mathrm {d} S+2\sum _{i=1}^{N}H_{Pn_{i}}\,\mathrm {d} P\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}H_{Sn_{i}}\,\mathrm {d} S\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}\sum _{j>i}^{N}H_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\\&=\mathrm {d} V\,\mathrm {d} P+\mathrm {d} T\,\mathrm {d} S+\sum _{i=1}^{N}\mathrm {d} \mu _{i}\,\mathrm {d} n_{i}\end{aligned}}$
* enthalpie libre $G$ (d2G) :
       ${\begin{aligned}\mathrm {d} ^{2}G&=G_{PP}\,\left(\mathrm {d} P\right)^{2}+G_{TT}\,\left(\mathrm {d} T\right)^{2}+\sum _{i=1}^{N}G_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\,G_{PT}\,\mathrm {d} P\,\mathrm {d} T+2\sum _{i=1}^{N}G_{Pn_{i}}\,\mathrm {d} P\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}G_{Tn_{i}}\,\mathrm {d} T\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}\sum _{j>i}^{N}G_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\\&=\mathrm {d} V\,\mathrm {d} P-\mathrm {d} S\,\mathrm {d} T+\sum _{i=1}^{N}\mathrm {d} \mu _{i}\,\mathrm {d} n_{i}\end{aligned}}$

##### Notes

    Avec $\Pi =-P$ et $X$ une autre variable quelconque, on a pour toute fonction $\Phi$ :
   * $\mathrm {d} \Pi =-\mathrm {d} P$ ;
   * $\Phi _{\Pi X}=-\Phi _{PX}$, d'où $\Phi _{\Pi X}\,\mathrm {d} \Pi \,\mathrm {d} X=\Phi _{PX}\,\mathrm {d} P\,\mathrm {d} X$ ; les conditions de stabilité n'imposent pas le signe de ces dérivées croisées ;
   * $\Phi _{\Pi \Pi }=\Phi _{PP}$, d'où $\Phi _{\Pi \Pi }\,\left(\mathrm {d} \Pi \right)^{2}=\Phi _{PP}\,\left(\mathrm {d} P\right)^{2}$ ; $\Phi$ a la même convexité par rapport à $\Pi$ et $P$.
    Il est donc indifférent de travailler avec $\Pi$ ou $P$ pour déterminer la convexité de la fonction $\Phi$.

### Critères de convexité - stabilité
#### Convexité et extrémums

Une fonction $\Phi$ est convexe si sa matrice hessienne est définie positive en tout point. Entre autres méthodes, il suffit de montrer que, quelles que soient les valeurs des variables $X_{i}$ et leurs variations élémentaires $\mathrm {d} X_{i}$, on a :
    ${\begin{pmatrix}\mathrm {d} X_{1}&\mathrm {d} X_{2}&\cdots &\mathrm {d} X_{n}\end{pmatrix}}{\begin{pmatrix}\Phi _{X_{1}X_{1}}&\Phi _{X_{1}X_{2}}&\cdots &\Phi _{X_{1}X_{j}}&\cdots &\Phi _{X_{1}X_{n}}\\\Phi _{X_{2}X_{1}}&\Phi _{X_{2}X_{2}}&\cdots &\Phi _{X_{2}X_{j}}&\cdots &\Phi _{X_{2}X_{n}}\\\vdots &\vdots &\ddots &\vdots &\ddots &\vdots \\\Phi _{X_{i}X_{1}}&\Phi _{X_{i}X_{2}}&\cdots &\Phi _{X_{i}X_{j}}&\cdots &\Phi _{X_{i}X_{n}}\\\vdots &\vdots &\ddots &\vdots &\ddots &\vdots \\\Phi _{X_{n}X_{1}}&\Phi _{X_{n}X_{2}}&\cdots &\Phi _{X_{n}X_{j}}&\cdots &\Phi _{X_{n}X_{n}}\end{pmatrix}}{\begin{pmatrix}\mathrm {d} X_{1}\\\mathrm {d} X_{2}\\\vdots \\\mathrm {d} X_{n}\end{pmatrix}}=\sum _{i=1}^{n}\sum _{j=1}^{n}\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}\geq 0$
Cette forme quadratique selon les $\mathrm {d} X_{i}$ est la différentielle d'ordre deux de $\Phi$, d'où :
* Fonction convexe Une fonction $\Phi$ est convexe si et seulement si $\mathrm {d} ^{2}\Phi \geq 0$.
Si $\Phi$ est concave, sa matrice hessienne est définie négative, d'où :
* Fonction concave Une fonction $\Phi$ est concave si et seulement si $\mathrm {d} ^{2}\Phi \leq 0$.
Les conditions de stabilité établies ci-après à partir de ces formes quadratiques ne sont nécessaires et suffisantes, pour démontrer que la matrice hessienne est définie positive ou négative, que dans les cas de fonctions à une ou deux variables. Pour les fonctions à plus de deux variables elles ne sont que nécessaires, elles ne sont pas suffisantes.
Lorsque l'équilibre thermodynamique est atteint, la différentielle du premier ordre $\mathrm {d} \Phi =0$. Le signe de la différentielle de deuxième ordre $\mathrm {d} ^{2}\Phi$ donne la nature de l'extrémum atteint : maximum si négatif, minimum si positif. Si $\mathrm {d} ^{2}\Phi$ est nulle, il est nécessaire d'étudier les différentielles d'ordre supérieur ; ce cas n'est pas traité ici, toutefois :
* s'il existe $n$ entier naturel non nul tel que $\mathrm {d} ^{2n}\Phi \neq 0$, toutes les différentielles d'ordre inférieur étant nulles, alors :
** si $\mathrm {d} ^{2n}\Phi <0$, alors $\Phi$ est concave, l'extrémum est un maximum ;
** si $\mathrm {d} ^{2n}\Phi >0$, alors $\Phi$ est convexe, l'extrémum est un minimum ;
* s'il existe $n$ entier naturel non nul tel que $\mathrm {d} ^{2n+1}\Phi \neq 0$, quel que soit son signe, toutes les différentielles d'ordre inférieur étant nulles, alors le point est un point d'inflexion, l'équilibre est métastable ;
* si quel que soit $n$ entier naturel non nul $\mathrm {d} ^{n}\Phi =0$, alors l'équilibre est neutre, indifférent à toute modification de ses conditions d'établissement.

#### Comparaison de convexités

Soient deux potentiels thermodynamiques $\Phi$ et $\Psi$ liés par une transformation de Legendre : $\Psi =\Phi -ZY$. On a les différentielles du deuxième ordre :
    $\mathrm {d} ^{2}\Phi =\mathrm {d} \Phi _{Y}\,\mathrm {d} Y+\sum _{i=1}^{n}\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}$
    $\mathrm {d} ^{2}\Psi =\mathrm {d} \Psi _{Z}\,\mathrm {d} Z+\sum _{i=1}^{n}\mathrm {d} \Psi _{X_{i}}\,\mathrm {d} X_{i}$
Avec les relations $\Phi _{Y}=Z$, $\Psi _{Z}=-Y$ et $\Psi _{X_{i}}=\Phi _{X_{i}}$ (voir la section Transformation de Legendre), on obtient :
* $\mathrm {d} ^{2}\Phi =\mathrm {d} Z\,\mathrm {d} Y+\sum _{i=1}^{n}\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}$ $\mathrm {d} ^{2}\Psi =-\mathrm {d} Y\,\mathrm {d} Z+\sum _{i=1}^{n}\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}$
Si l'on connait le signe de $\mathrm {d} ^{2}\Phi$, alors on connait le signe de $\mathrm {d} ^{2}\Psi$ sous certaines conditions. Pour cela, on fixe certaines variables de $\Psi$ :
* à $X_{i}$ constantes, soit $\mathrm {d} X_{i}=0$ pour tout $i$ :
       $\mathrm {d} ^{2}\Psi =-\mathrm {d} Y\,\mathrm {d} Z=-\mathrm {d} ^{2}\Phi$
    Dans ces conditions les deux différentielles du deuxième ordre sont de signes opposés. La transformation de Legendre inverse la convexité de la fonction : si $\Phi$ est convexe par rapport à sa variable $Y$, alors $\Psi$ est concave par rapport à sa variable $Z$ ;
* à $Z$ constante, soit $\mathrm {d} Z=0$ :
       $\mathrm {d} ^{2}\Psi =\sum _{i=1}^{n}\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}=\mathrm {d} ^{2}\Phi$
    Dans ces conditions les deux différentielles du deuxième ordre sont de même signe. Si $\Phi$ est convexe par rapport à ses variables $X_{i}$, alors $\Psi$ l'est aussi, ces variables étant communes aux deux potentiels.
Si $\Phi$ est convexe par rapport à toutes ses variables, alors $\Psi$ est concave par rapport à la variable transformée $Z$ et convexe par rapport aux variables communes $X_{i}$, elle est convexe-concave. Pour une variable commune $X_{i}$ quelconque, les autres variables $X_{\neq i}$ étant fixées, on a $\mathrm {d} ^{2}\Phi =\mathrm {d} Z\,\mathrm {d} Y+\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}$ et $\mathrm {d} ^{2}\Psi =-\mathrm {d} Y\,\mathrm {d} Z+\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}$. Si l'on connait le signe de $\mathrm {d} ^{2}\Phi$, on ne peut en déduire celui de $\mathrm {d} ^{2}\Psi$, qui est donc indéterminé dans ces conditions.
Les conditions de stabilité par rapport aux potentiels thermodynamiques $F$, $H$ et $G$ s'obtiennent en comparant leur différentielle du deuxième ordre à celle de l'énergie interne $U$ que l'on sait convexe par rapport à toutes ses variables naturelles, soit $\mathrm {d} ^{2}U\geq 0$. Pour un potentiel thermodynamique $\Psi$ quelconque, il suffit de trouver dans quelles conditions $\mathrm {d} ^{2}\Psi =\pm \,\mathrm {d} ^{2}U$. Les variables de l'énergie interne $U$ sont toutes extensives. La transformation de Legendre remplace une variable extensive $Y$ de $U$ par une variable intensive $Z$ de $\Psi$. En conséquence, les potentiels thermodynamiques sont tous concaves par rapport à leur(s) variable(s) intensive(s)$^,$. Les autres variables $X_{i}$ communes à $U$ et $\Psi$ sont donc extensives. En conséquence, les potentiels thermodynamiques sont tous convexes par rapport à leurs variables extensives$^,$.

#### Fonction convexe

Soit une fonction quelconque $\Phi$ convexe par rapport à sa seule variable $X_{i}$ (toutes ses autres variables étant fixées, $\mathrm {d} X_{\neq i}=0$) :
    $\mathrm {d} \Phi _{X_{i}}=\Phi _{X_{i}X_{i}}\,\mathrm {d} X_{i}$
    $\mathrm {d} ^{2}\Phi =\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}=\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}\geq 0$
On obtient immédiatement la condition de stabilité :
* Fonction convexe - condition de stabilité si $\mathrm {d} ^{2}\Phi =\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}\geq 0$ alors $\Phi _{X_{i}X_{i}}\geq 0$
Pour une fonction quelconque $\Phi$ convexe par rapport à ses deux variables $X_{i}$ et $X_{j}$ (toutes ses autres variables étant fixées, $\mathrm {d} X_{\neq i,j}=0$) :
    $\mathrm {d} \Phi _{X_{i}}=\Phi _{X_{i}X_{i}}\,\mathrm {d} X_{i}+\Phi _{X_{j}X_{i}}\,\mathrm {d} X_{j}$ ; $\mathrm {d} \Phi _{X_{j}}=\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}+\Phi _{X_{j}X_{j}}\,\mathrm {d} X_{j}$
    $\Phi _{X_{j}X_{i}}=\Phi _{X_{i}X_{j}}$ (théorème de Schwarz)
    $\mathrm {d} ^{2}\Phi =\mathrm {d} \Phi _{X_{i}}\,\mathrm {d} X_{i}+\mathrm {d} \Phi _{X_{j}}\,\mathrm {d} X_{j}=\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}+2\,\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}+\Phi _{X_{j}X_{j}}\,\left(\mathrm {d} X_{j}\right)^{2}\geq 0$
À $X_{j}$ constant on a $\mathrm {d} ^{2}\Phi =\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}\geq 0$ et par conséquent la condition de stabilité $\Phi _{X_{i}X_{i}}\geq 0$. Le même raisonnement mené à $X_{i}$ constant conduit à la condition de stabilité $\Phi _{X_{j}X_{j}}\geq 0$. Par ailleurs, quelles que soient les variations élémentaires $\mathrm {d} X_{i}$ et $\mathrm {d} X_{j}$, la différentielle $\mathrm {d} ^{2}\Phi$ ne peut pas prendre de valeur négative. On pose$^,$ :
    $x=\mathrm {d} X_{i}$ ; $a=\Phi _{X_{i}X_{i}}$ ; $b=2\,\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{j}$ ; $c=\Phi _{X_{j}X_{j}}\,\left(\mathrm {d} X_{j}\right)^{2}$
    $\mathrm {d} ^{2}\Phi =a\,x^{2}+b\,x+c\geq 0$
Le polynôme du second degré caractéristique ainsi obtenu ne peut avoir que des valeurs positives ou nulles. Il ne peut avoir par conséquent que des racines complexes (le polynôme ne s'annule jamais avec des valeurs de $x$ réelles) ou une racine réelle double (le polynôme s'annule et atteint simultanément un minimum à une racine réelle double). Son discriminant répond donc à la contrainte $b^{2}-4ac\leq 0$. On obtient la condition de stabilité : ${\Phi _{X_{i}X_{j}}}^{2}\leq \Phi _{X_{i}X_{i}}\,\Phi _{X_{j}X_{j}}$. Ce résultat est indépendant du choix de $x=\mathrm {d} X_{i}$ ou $x=\mathrm {d} X_{j}$. Cette condition ne renseigne pas sur le signe de la dérivée croisée $\Phi _{X_{i}X_{j}}$, mais elle donne une indication sur sa valeur absolue : $\left|\Phi _{X_{i}X_{j}}\right|\leq {\sqrt {\Phi _{X_{i}X_{i}}\,\Phi _{X_{j}X_{j}}}}$. Par le théorème de Schwarz cette condition est identique pour $\Phi _{X_{j}X_{i}}=\Phi _{X_{i}X_{j}}$. Pour une fonction convexe par rapport à deux variables, on a donc les conditions de stabilité (théorème de Monge)$^,$$^,$ :
* Fonction convexe - conditions de stabilité si $\mathrm {d} ^{2}\Phi =\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}+2\,\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}+\Phi _{X_{j}X_{j}}\,\left(\mathrm {d} X_{j}\right)^{2}\geq 0$ alors $\Phi _{X_{i}X_{i}}\geq 0$, $\Phi _{X_{j}X_{j}}\geq 0$ et ${\Phi _{X_{i}X_{j}}}^{2}\leq \Phi _{X_{i}X_{i}}\,\Phi _{X_{j}X_{j}}$
Pour une fonction convexe par rapport à plus de deux variables, on revient au cas à deux variables en fixant les variables surnuméraires.
Les potentiels thermodynamiques $U$, $F$, $H$ et $G$ sont tous convexes par rapport à leurs variables extensives. Les variables de l'énergie interne $U$ sont toutes extensives.

#### Fonction concave

Le raisonnement pour une fonction concave est similaire à celui pour une fonction convexe. Sa différentielle du deuxième ordre ne peut avoir que des valeurs négatives ou nulles.
Pour une fonction concave par rapport à une seule variable, on a :
* Fonction concave - condition de stabilité si $\mathrm {d} ^{2}\Phi =\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}\leq 0$ alors $\Phi _{X_{i}X_{i}}\leq 0$
Pour une fonction concave par rapport à deux variables et plus, on a (théorème de Monge)$^,$$^,$ :
* Fonction concave - conditions de stabilité si $\mathrm {d} ^{2}\Phi =\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}+2\,\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}+\Phi _{X_{j}X_{j}}\,\left(\mathrm {d} X_{j}\right)^{2}\leq 0$ alors $\Phi _{X_{i}X_{i}}\leq 0$, $\Phi _{X_{j}X_{j}}\leq 0$ et ${\Phi _{X_{i}X_{j}}}^{2}\leq \Phi _{X_{i}X_{i}}\,\Phi _{X_{j}X_{j}}$
La dernière condition sur les dérivées croisées est identique à celle de la fonction convexe : elle impose que le polynôme caractéristique ne puisse changer de signe.
L'entropie $S$ est concave par rapport à toutes ses variables, qui sont toutes extensives. Les potentiels thermodynamiques $F$, $H$ et $G$ sont concaves par rapport à leurs variables intensives (l'énergie interne $U$ n'en a pas).

#### Fonction convexe-concave

Il est toujours possible de déterminer le signe de la dérivée seconde $\Phi _{X_{i}X_{i}}$ d'une fonction thermodynamique $\Phi$ (entropie ou potentiel thermodynamique) ; on sait donc si la fonction est convexe ou concave par rapport à cette seule variable $X_{i}$.
Si $\Phi$ est convexe par rapport à $X_{i}$, soit $\Phi _{X_{i}X_{i}}\geq 0$, et concave par rapport à $X_{j}$, soit $\Phi _{X_{j}X_{j}}\leq 0$, alors la fonction est convexe-concave. Le polynôme caractéristique doit pouvoir changer de signe et son discriminant est tel que $b^{2}-4ac\geq 0$. On obtient :
* Fonction convexe-concave - condition de stabilité si $\mathrm {d} ^{2}\Phi =\Phi _{X_{i}X_{i}}\,\left(\mathrm {d} X_{i}\right)^{2}+2\,\Phi _{X_{i}X_{j}}\,\mathrm {d} X_{i}\,\mathrm {d} X_{j}+\Phi _{X_{j}X_{j}}\,\left(\mathrm {d} X_{j}\right)^{2}$ indéterminé avec $\Phi _{X_{i}X_{i}}\geq 0$ et $\Phi _{X_{j}X_{j}}\leq 0$ alors ${\Phi _{X_{i}X_{j}}}^{2}\geq \Phi _{X_{i}X_{i}}\,\Phi _{X_{j}X_{j}}$
Cette dernière condition ne renseigne ni sur le signe, ni sur la valeur absolue de $\Phi _{X_{i}X_{j}}$, puisque ${\Phi _{X_{i}X_{j}}}^{2}\geq 0$ et $\Phi _{X_{i}X_{i}}\,\Phi _{X_{j}X_{j}}\leq 0$.
Ce cas est rencontré si $X_{i}$ est extensive et $X_{j}$ intensive, pour les potentiels thermodynamiques énergie libre $F$, enthalpie $H$ et enthalpie libre $G$.

### Établissement des critères de stabilité

Les conditions établies ci-après ne sont nécessaires et suffisantes, pour établir la stabilité du système thermodynamique, que dans le cas de transformations n'impliquant qu'une seule ou deux variables thermodynamiques. Pour des transformations impliquant plus de variables, elles ne sont que nécessaires à la stabilité, elles ne sont pas suffisantes.

#### Critères selon l'entropie

La relation (d2S) donne, à composition constante, soit $\mathrm {d} n_{i}=0$ pour toute espèce chimique $i$ :
    $\mathrm {d} ^{2}S=S_{VV}\,\left(\mathrm {d} V\right)^{2}+S_{UU}\,\left(\mathrm {d} U\right)^{2}+2\,S_{VU}\,\mathrm {d} V\,\mathrm {d} U$
En conséquence, étant donné la concavité de l'entropie, soit $\mathrm {d} ^{2}S\leq 0$, on obtient immédiatement$^,$ :
* (csS1) : $S_{VV}\leq 0$ ; (csS2) : $S_{UU}\leq 0$ (csS4) : ${S_{VU}}^{2}\leq S_{VV}\,S_{UU}$
En considérant tous les autres jeux de deux variables possibles de $S$, on obtient de la même façon, pour toutes espèces chimiques $i$ et $j$ :
$S_{n_{i}n_{i}}\leq 0$
${S_{Vn_{i}}}^{2}\leq S_{VV}\,S_{n_{i}n_{i}}$
${S_{Un_{i}}}^{2}\leq S_{UU}\,S_{n_{i}n_{i}}$
${S_{n_{i}n_{j}}}^{2}\leq S_{n_{i}n_{i}}\,S_{n_{j}n_{j}}$
L'entropie $S$ est concave par rapport à ses variables $V$, $U$ et $n_{i}$, qui sont toutes extensives.

#### Critères selon l'énergie interne

La relation (d2U) donne, à composition constante, soit $\mathrm {d} n_{i}=0$ pour toute espèce chimique $i$ :
    $\mathrm {d} ^{2}U=U_{VV}\,\left(\mathrm {d} V\right)^{2}+U_{SS}\,\left(\mathrm {d} S\right)^{2}+2\,U_{VS}\,\mathrm {d} V\,\mathrm {d} S$
En conséquence, étant donné la convexité de l'énergie interne, soit $\mathrm {d} ^{2}U\geq 0$, on obtient immédiatement$^,$ :
* (csU1) : $U_{VV}\geq 0$ ; (csU2) : $U_{SS}\geq 0$ (csU4) : ${U_{VS}}^{2}\leq U_{VV}\,U_{SS}$
En considérant tous les autres jeux de deux variables possibles de $U$, on obtient de la même façon, pour toutes espèces chimiques $i$ et $j$ :
$U_{n_{i}n_{i}}\geq 0$
${U_{Vn_{i}}}^{2}\leq U_{VV}\,U_{n_{i}n_{i}}$
${U_{Sn_{i}}}^{2}\leq U_{SS}\,U_{n_{i}n_{i}}$
${U_{n_{i}n_{j}}}^{2}\leq U_{n_{i}n_{i}}\,U_{n_{j}n_{j}}$
L'énergie interne $U$ est convexe par rapport à ses variables $V$, $S$ et $n_{i}$, qui sont toutes extensives.

#### Critères selon l'énergie libre

Les relations (d2F) et (d2U) donnent, à température constante, soit $\mathrm {d} T=0$ :
    ${\begin{aligned}\mathrm {d} ^{2}F&=-\mathrm {d} P\,\mathrm {d} V+\sum _{i=1}^{N}\mathrm {d} \mu _{i}\,\mathrm {d} n_{i}=\mathrm {d} ^{2}U\\&=F_{VV}\,\left(\mathrm {d} V\right)^{2}+\sum _{i=1}^{N}F_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\sum _{i=1}^{N}\,F_{Vn_{i}}\,\mathrm {d} V\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}\sum _{j>i}^{N}\,F_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\geq 0\end{aligned}}$
On obtient immédiatement, pour toutes espèces chimiques $i$ et $j$$^,$ :
$F_{VV}\geq 0$
$F_{n_{i}n_{i}}\geq 0$
${F_{Vn_{i}}}^{2}\leq F_{VV}\,F_{n_{i}n_{i}}$
${F_{n_{i}n_{j}}}^{2}\leq F_{n_{i}n_{i}}\,F_{n_{j}n_{j}}$
À volume et composition constants, soit $\mathrm {d} V=0$ et $\mathrm {d} n_{i}=0$ pour toute espèce chimique $i$ :
    ${\begin{aligned}\mathrm {d} ^{2}F&=-\mathrm {d} S\,\mathrm {d} T=-\mathrm {d} ^{2}U\\&=F_{TT}\,\left(\mathrm {d} T\right)^{2}\leq 0\end{aligned}}$
On obtient immédiatement$^,$ :
$F_{TT}\leq 0$
L'énergie libre $F$ est convexe par rapport à ses variables extensives $V$ et $n_{i}$, et concave par rapport à sa variable intensive $T$ ; on n'a donc aucune condition de stabilité concernant $F_{VT}$ et $F_{Tn_{i}}$.

#### Critères selon l'enthalpie

Les relations (d2H) et (d2U) donnent, à pression constante, soit $\mathrm {d} P=0$ :
    ${\begin{aligned}\mathrm {d} ^{2}H&=\mathrm {d} T\,\mathrm {d} S+\sum _{i=1}^{N}\mathrm {d} \mu _{i}\,\mathrm {d} n_{i}=\mathrm {d} ^{2}U\\&=H_{SS}\,\left(\mathrm {d} S\right)^{2}+\sum _{i=1}^{N}H_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\sum _{i=1}^{N}\,H_{Sn_{i}}\,\mathrm {d} S\,\mathrm {d} n_{i}+2\sum _{i=1}^{N}\sum _{j>i}^{N}\,H_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\geq 0\end{aligned}}$
On obtient immédiatement, pour toutes espèces chimiques $i$ et $j$ :
$H_{SS}\geq 0$
$H_{n_{i}n_{i}}\geq 0$
${H_{Sn_{i}}}^{2}\leq H_{SS}\,H_{n_{i}n_{i}}$
${H_{n_{i}n_{j}}}^{2}\leq H_{n_{i}n_{i}}\,H_{n_{j}n_{j}}$
À entropie et composition constantes, soit $\mathrm {d} S=0$ et $\mathrm {d} n_{i}=0$ pour toute espèce chimique $i$ :
    ${\begin{aligned}\mathrm {d} ^{2}H&=\mathrm {d} V\,\mathrm {d} P=-\mathrm {d} ^{2}U\\&=H_{PP}\,\left(\mathrm {d} P\right)^{2}\leq 0\end{aligned}}$
On obtient immédiatement :
$H_{PP}\leq 0$
L'enthalpie $H$ est convexe par rapport à ses variables extensives $S$ et $n_{i}$, et concave par rapport à sa variable intensive $P$ ; on n'a donc aucune condition de stabilité concernant $H_{PS}$ et $H_{Pn_{i}}$.

#### Critères selon l'enthalpie libre

Les relations (d2G) et (d2U) donnent, à pression et température constantes, soit $\mathrm {d} P=0$ et $\mathrm {d} T=0$ :
    ${\begin{aligned}\mathrm {d} ^{2}G&=\sum _{i=1}^{N}\mathrm {d} \mu _{i}\,\mathrm {d} n_{i}=\mathrm {d} ^{2}U\\&=\sum _{i=1}^{N}G_{n_{i}n_{i}}\,\left(\mathrm {d} n_{i}\right)^{2}+2\sum _{i=1}^{N}\sum _{j>i}^{N}\,G_{n_{i}n_{j}}\,\mathrm {d} n_{i}\,\mathrm {d} n_{j}\geq 0\end{aligned}}$
On obtient immédiatement, pour toutes espèces chimiques $i$ et $j$ :
* (csG3) : ${G_{n_{i}n_{i}}}\geq 0$ (csG5) : ${G_{n_{i}n_{j}}}^{2}\leq G_{n_{i}n_{i}}\,G_{n_{j}n_{j}}$
À composition constante, soit $\mathrm {d} n_{i}=0$ pour toute espèce chimique $i$ :
    ${\begin{aligned}\mathrm {d} ^{2}G&=\mathrm {d} V\,\mathrm {d} P-\mathrm {d} S\,\mathrm {d} T=-\mathrm {d} ^{2}U\\&=G_{PP}\,\left(\mathrm {d} P\right)^{2}+G_{TT}\,\left(\mathrm {d} T\right)^{2}+2\,G_{PT}\,\mathrm {d} P\,\mathrm {d} T\leq 0\end{aligned}}$
On obtient immédiatement$^,$$^,$ :
* (csG1) : $G_{PP}\leq 0$ ; (csG2) : $G_{TT}\leq 0$ (csG4) : ${G_{PT}}^{2}\leq G_{PP}\,G_{TT}$
L'enthalpie libre $G$ est convexe par rapport à ses variables extensives $n_{i}$ et concave par rapport à ses variables intensives $P$ et $T$ ; on n'a donc aucune condition de stabilité concernant $G_{Pn_{i}}$ et $G_{Tn_{i}}$.