Skip to content

Phase Transitions

A phase transition is a discontinuity in a thermodynamic quantity or its derivative as a function of a state variable. Phase transitions are classified by Ehrenfest according to which derivative of the Gibbs free energy is discontinuous.

10.1 Classification of Phase Transitions

OrderDefinitionExample
First orderGG continuous; G/T\partial G/\partial T or G/P\partial G/\partial P discontinuousBoiling of water
Second orderFirst derivatives continuous; second derivatives discontinuousSuperconducting transition
Lambda (λ\lambda)Divergent second derivativesHelium-4 superfluid transition

For a first-order transition at temperature TcT_cThe latent heat is:

L=TcΔS=Tc(Sphase2Sphase1)L = T_c \Delta S = T_c \left(S_{\text{phase} 2} - S_{\text{phase} 1}\right)

The Clausius—Clapeyron equation governs the slope of the coexistence curve:

dPdT=LTcΔv\frac{dP}{dT} = \frac{L}{T_c \Delta v}

Where Δv=v2v1\Delta v = v_2 - v_1 is the change in specific volume.

10.2 Van der Waals Equation and Critical Phenomena

The van der Waals equation of state modifies the ideal gas law to account for intermolecular forces:

(P+av2)(vb)=kBT\left(P + \frac{a}{v^2}\right)(v - b) = k_B T

Where aa accounts for attractive interactions and bb for the finite molecular volume. The critical point (Tc,Pc,vc)(T_c, P_c, v_c) satisfies:

PvTc=0,2Pv2Tc=0\frac{\partial P}{\partial v}\bigg|_{T_c} = 0, \qquad \frac{\partial^2 P}{\partial v^2}\bigg|_{T_c} = 0

Solving gives:

Tc=8a27bkB,Pc=a27b2,vc=3bT_c = \frac{8a}{27bk_B}, \qquad P_c = \frac{a}{27b^2}, \qquad v_c = 3b

Near the critical point, define the reduced variables T~=T/Tc\tilde{T} = T/T_c, P~=P/Pc\tilde{P} = P/P_c, v~=v/vc\tilde{v} = v/v_c to obtain the universal form:

(P~+3v~2)(3v~1)=8T~\left(\tilde{P} + \frac{3}{\tilde{v}^2}\right)(3\tilde{v} - 1) = 8\tilde{T}

The order parameter ϕ=(vgasvliquid)/(vc)\phi = (v_{\text{gas} - v_{\text{liquid})/(v_c)}} vanishes as:

ϕ(TcT)β\phi \propto (T_c - T)^{\beta}

Where β=1/2\beta = 1/2 is the mean-field critical exponent (van der Waals prediction).

10.3 Critical Exponents

Near a second-order phase transition, thermodynamic quantities follow power laws characterized by critical exponents:

ExponentDefinitionMean-field2D Ising3D Ising (numerical)
α\alphaCtαC \propto \|t\|^{-\alpha}0 (jump)0 (log)
β\betaϕ(t)β\phi \propto (-t)^\beta1/21/21/81/80.326\approx 0.326
γ\gammaχtγ\chi \propto \|t\|^{-\gamma}117/47/4
δ\deltaϕh1/δ\phi \propto h^{1/\delta} at t=0t=03315154.789\approx 4.789

Where t=(TTc)/Tct = (T - T_c)/T_c is the reduced temperature and hh is the conjugate field.

Worked Example 10.1: Clausius--Clapeyron for Water

For the water—steam transition at 1 atm, Tc=373.15T_c = 373.15 K, L=2260L = 2260 kJ/kg, vsteam=1.673v_{\text{steam} = 1.673} m3^3/kg, vwater=1.043×103v_{\text{water} = 1.043 \times 10^{-3}} m3^3/kg.

dPdT=LTΔv=2.26×106373.15×1.673=2.26×106624.33620 Pa/K0.0357 atm/K\frac{dP}{dT} = \frac{L}{T \Delta v} = \frac{2.26 \times 10^6}{373.15 \times 1.673} = \frac{2.26 \times 10^6}{624.3} \approx 3620 \text{ Pa/K} \approx 0.0357 \text{ atm/K}

This means increasing the boiling temperature by 1 K requires increasing the pressure by about 0.036 atm.

Worked Example 10.2: Critical Parameters of CO$_2$

For CO2_2, a=0.364a = 0.364 Pa\cdotM6^6/mol2^2, b=4.27×105b = 4.27 \times 10^{-5} m3^3/mol. Using the critical point formulas:

Tc=8a27Rb=8×0.36427×8.314×4.27×105=2.9129.585×103303.7 KT_c = \frac{8a}{27Rb} = \frac{8 \times 0.364}{27 \times 8.314 \times 4.27 \times 10^{-5}} = \frac{2.912}{9.585 \times 10^{-3}} \approx 303.7 \text{ K}

Pc=a27b2=0.36427×(4.27×105)2=0.3644.923×1087.40×106 Pa=74.0 atmP_c = \frac{a}{27b^2} = \frac{0.364}{27 \times (4.27 \times 10^{-5})^2} = \frac{0.364}{4.923 \times 10^{-8}} \approx 7.40 \times 10^6 \text{ Pa} = 74.0 \text{ atm}

The experimental values are Tc=304.3T_c = 304.3 K and Pc=73.8P_c = 73.8 atm, showing good agreement.