# Phase coexistence and general properties of phase transitions

In order to understand a little bit more about phase transitions and the coexistence of different phases, let us consider a ${\displaystyle PVT}$ system in the liquid-gaseous transition, and let us focus on the region of liquid-gas coexistence; physically, what happens (and what can also be understood from the figure) is that the system exhibits both a liquid and a gaseous phase, and it is possible to change the volume of the system without changing its pressure at fixed temperature. This means that the work we do on the system is used only to change the proportion between the phases (breaking or forming molecular bonds, on a microscopic level).

Let us now consider a slightly more general situation: suppose we have an isolated single-component thermodynamic system described by the internal energy ${\displaystyle U}$ and extensive variables ${\displaystyle x_{i}}$ (generalized displacements); if the system is subject to reversible processes then from the expression of ${\displaystyle dU}$ we get:

${\displaystyle dS={\frac {1}{T}}dU-\sum _{i}{\frac {X_{i}}{T}}dx_{i}}$
where ${\displaystyle X_{i}}$ are the generalized forces relative to the generalized displacements ${\displaystyle x_{i}}$. Let us now suppose that two phases ${\displaystyle \alpha }$ and ${\displaystyle \beta }$ coexist and that they can exchange ${\displaystyle U}$ and ${\displaystyle x}$ (namely, they are systems in contact with each other). At equilibrium the total entropy ${\displaystyle S=S_{\alpha }+S_{\beta }}$ must be maximized, and since the whole system is in equilibrium we also have ${\displaystyle dU_{\alpha }+dU_{\beta }=0}$ and ${\displaystyle dx_{i,\alpha }+dx_{i,\beta }=0}$. Substituting in ${\textstyle dS}$ we get:
${\displaystyle dS=\left({\frac {1}{T_{\alpha }}}-{\frac {1}{T_{\beta }}}\right)dU_{\alpha }-\sum _{i}\left({\frac {X_{i,\alpha }}{T_{\alpha }}}-{\frac {X_{i,\beta }}{T_{\beta }}}\right)dx_{i,\alpha }=0}$
Therefore, since ${\displaystyle dU_{\alpha }}$ and ${\displaystyle dx_{i,\alpha }}$ are arbitrary positive quantities, we must have:
${\displaystyle T_{\alpha }=T_{\beta }\quad \qquad X_{i,\alpha }=X_{i,\beta }}$

Considering our ${\displaystyle PVT}$ system again, since the generalized displacements that are needed in order to describe it are the volume ${\displaystyle V}$ and the number of particles ${\displaystyle N}$, when two phases coexist we have:

${\displaystyle T_{\alpha }=T_{\beta }\quad \quad P_{\alpha }=P_{\beta }\quad \quad \mu _{\alpha }=\mu _{\beta }}$
Since, as we have shown in Thermodynamic potentials, the Gibbs free energy is ${\displaystyle G(T,P)=\mu N}$ (for a single-component ${\displaystyle PVT}$ system) then ${\displaystyle g(T,P)=\mu }$ and so when two phases coexist we also have:
${\displaystyle g_{\alpha }(T,P)=g_{\beta }(T,P)}$
This equality must hold along the whole coexistence line in ${\displaystyle (T,P)}$ space and so if we know a point in this space where the two phases coexist we can, at least locally, "reconstruct" the coexistence line. In fact we have ${\displaystyle dg_{\alpha }=dg_{\beta }}$, namely:
{\displaystyle {\begin{aligned}dg_{\alpha }=-s_{\alpha }dT+v_{\alpha }dP=dg_{\beta }=-s_{\beta }dT+v_{\beta }dP\quad \Rightarrow \\\Rightarrow \quad {\frac {dP}{dT}}_{|{\text{coex.}}}={\frac {s_{\alpha }-s_{\beta }}{v_{\alpha }-v_{\beta }}}={\frac {L_{\alpha ,\beta }}{T(v_{\alpha }-v_{\beta })}}\end{aligned}}}
where ${\displaystyle s=S/N}$ and ${\displaystyle v=V/N}$, and by definition ${\displaystyle L_{\alpha ,\beta }}$ is the molar latent heat needed to bring the system from phase ${\displaystyle \beta }$ to phase ${\displaystyle \alpha }$. This is known as Clausius-Clapeyron equation.

From the expression of ${\displaystyle dg_{\alpha }=dg_{\beta }}$ we can also understand some very general properties of phase transitions. In fact, from the expressions of ${\textstyle dg_{\alpha }}$ and ${\textstyle dg_{\beta }}$ at phase coexistence we have:

${\displaystyle {\frac {\partial g_{\beta }}{\partial T}}_{|P}-{\frac {\partial g_{\alpha }}{\partial T}}_{|P}=\Delta s>0\quad \quad {\frac {\partial g_{\beta }}{\partial P}}_{|T}-{\frac {\partial g_{\alpha }}{\partial P}}_{|T}=\Delta v>0}$
This means that when the system undergoes a phase transition its volume and its entropy have a jump discontinuity; as we will later see in more detail, since in this transition the first derivatives of a thermodynamic potential have a jump discontinuity, we call it a first order transition. A similar behaviour can be encountered in magnetic systems, where the magnetization ${\displaystyle M}$ has a jump at ${\displaystyle H=0}$ for temperatures lower that the critical one; in this case since ${\displaystyle M=-\partial F/\partial H}$ we see that the first derivative of the free energy ${\displaystyle F}$ with respect to ${\displaystyle H}$ has a jump discontinuity.

Furthermore, if we consider our system at the critical point we see from the ${\displaystyle (V,P)}$ projection of the phase diagram that the isothermal compressibility:

${\displaystyle K_{T}=-{\frac {1}{V}}{\frac {\partial V}{\partial P}}_{|T}}$
diverges when ${\displaystyle V=V_{c}}$. Similarly, the magnetic susceptibility ${\displaystyle \chi _{T}}$ of a magnetic system at its critical temperature diverges when ${\displaystyle H=0}$ and ${\displaystyle M=0}$. As we will see further on in Long range correlations, the divergence of response functions is a typical behaviour of thermodynamic systems in the neighbourhood of critical points and has important consequences.