Technopedia

The technologist’s encyclopedia

Phase Transitions

Phase transitions are abrupt changes in the macroscopic state of a system as a control parameter (temperature, pressure, field) crosses a critical value. They are relevant whenever collective behavior emerges: water freezing, magnets ordering, superconductors forming, fluids separating, or any system where a small parameter change produces a qualitative structural shift. They are the domain where statistical mechanics meets materials science, cosmology, and critical phenomena.

Authoritative sources

Order parameter

A phase transition is characterized by an order parameter ϕ\phi that is zero in the symmetric (disordered) phase and nonzero in the broken-symmetry (ordered) phase:

ϕ=0(T>Tc),ϕ≠0(T<Tc)\phi = 0 \quad (T > T_c), \qquad \phi \neq 0 \quad (T < T_c)

Examples: magnetization MM (ferromagnet), density difference ρl−ρg\rho_l - \rho_g (liquid–gas), superconducting gap Δ\Delta.

Details

Global idea. The order parameter is the macroscopic variable that distinguishes the two phases: it measures how much symmetry the disordered phase has that the ordered phase lacks. It transforms under the symmetry group that is broken. A ferromagnet breaks Z2\mathbb{Z}_2 spin-flip symmetry; a liquid–gas transition breaks the density-translation symmetry. The choice of order parameter is not unique — any quantity that vanishes in the disordered phase and is nonzero in the ordered phase works.

The Landau picture. Write the free energy as an analytic function of the order parameter, constrained by symmetry: F(ϕ)=a(T−Tc)ϕ2+bϕ4+⋯F(\phi) = a(T-T_c)\phi^2 + b\phi^4 + \cdots with b>0b > 0. Above TcT_c the coefficient of ϕ2\phi^2 is positive, so the minimum sits at ϕ=0\phi = 0 (disordered). Below TcT_c it is negative, so two symmetric minima appear at ϕ=±a(Tc−T)/(2b)\phi = \pm\sqrt{a(T_c-T)/(2b)} (ordered), giving the mean-field prediction ϕ∼(Tc−T)1/2\phi \sim (T_c - T)^{1/2}.

Example. For a ferromagnet, ϕ\phi is the net magnetization. Above the Curie temperature, thermal fluctuations point the atomic spins in random directions and M=0M = 0. Below it, the spin-spin coupling wins and the system spontaneously picks one of the two minima — a symmetry breaking choice that the Hamiltonian itself does not make.

T > Tc T = Tc T < Tc T < T c −φ₀ +φ₀ order parameter φ

The double well below TcT_c is the whole story of spontaneous symmetry breaking: the system must choose one minimum, and an infinite system cannot tunnel between them.

Authoritative sources. Landau’s original 1937 paper (republished, English translation); Cardy, Scaling and Renormalization in Statistical Physics, ch. 4.

Example problem

For the Landau free energy F(ϕ)=a(T−Tc)ϕ2+bϕ4F(\phi) = a(T-T_c)\phi^2 + b\phi^4 with a,b>0a, b > 0, find the equilibrium order parameter at T=0.96 TcT = 0.96\,T_c in terms of aa and bb, and verify it vanishes as (Tc−T)1/2(T_c - T)^{1/2}.

Solution

Minimize: dF/dϕ=2a(T−Tc)ϕ+4bϕ3=0dF/d\phi = 2a(T-T_c)\phi + 4b\phi^3 = 0. Nonzero solutions: ϕ2=a(Tc−T)/(2b)\phi^2 = a(T_c - T)/(2b), so ϕ=±a⋅0.04 Tc / (2b)=±0.02 aTc/b\phi = \pm\sqrt{a \cdot 0.04\,T_c\,/\,(2b)} = \pm\sqrt{0.02\,a T_c / b}. The square-root dependence on (Tc−T)(T_c - T) is the mean-field exponent β=1/2\beta = 1/2 in action.

Critical exponents

Near the critical point TcT_c, thermodynamic quantities vanish or diverge as power laws of the reduced temperature t=(T−Tc)/Tct = (T - T_c)/T_c:

ϕ∼∣t∣β,C∼∣t∣−α,χ∼∣t∣−γ,ξ∼∣t∣−ν\phi \sim |t|^{\beta}, \quad C \sim |t|^{-\alpha}, \quad \chi \sim |t|^{-\gamma}, \quad \xi \sim |t|^{-\nu}

where CC is the specific heat, χ\chi the susceptibility and ξ\xi the correlation length. Only ϕ\phi vanishes at TcT_c (positive exponent β\beta); the response functions CC, χ\chi and ξ\xi all diverge (negative exponents). The exponents are universal: they depend only on the spatial dimension dd and the symmetry of the order parameter, not on microscopic details.

ExponentQuantityMean-field (Landau)2D Ising (exact)
αSpecific heat C ∝ |t|−α00 (logarithmic divergence)
βOrder parameter φ ∝ |t|β1/21/8
γSusceptibility χ ∝ |t|−γ17/4
νCorrelation length ξ ∝ |t|−ν1/21
Details

Global idea. Power laws are the fingerprint of scale invariance. The correlation length ξ\xi is the distance over which fluctuations of the order parameter are correlated. At the critical point, ξ→∞\xi \to \infty, and the system becomes scale-invariant: fluctuations exist on all length scales, which is exactly why a power law (the only scale-free functional form) describes the singularities. This is why critical phenomena are described by the renormalization group (RG): integrating out short-distance modes reveals the same physics at longer scales.

Example. Cool a ferromagnet toward its Curie point from above. Domains of aligned spins grow from a few lattice spacings to macroscopic size (ξ→∞\xi \to \infty); the susceptibility χ\chi — the response to an external field — diverges because a whole correlated region flips coherently; the specific heat CC peaks (in 2D Ising it diverges only logarithmically, C∼ln⁡∣t∣C \sim \ln|t|, corresponding to α=0\alpha = 0).

Scaling relations. The exponents are not independent. The hyperscaling relation α+2β+γ=dν\alpha + 2\beta + \gamma = d\nu connects them, and it holds below the upper critical dimension. Mean-field theory (Landau) gives α=0,β=1/2,γ=1,ν=1/2\alpha = 0, \beta = 1/2, \gamma = 1, \nu = 1/2 and is exact only for d≥4d \geq 4 (the upper critical dimension). For the 2D Ising model, exact results (Onsager, Yang) give α=0,β=1/8,γ=7/4,ν=1\alpha = 0, \beta = 1/8, \gamma = 7/4, \nu = 1.

Authoritative sources. Onsager, Crystal Statistics I (1944); Goldenfeld, Lectures on Phase Transitions and the Renormalization Group, ch. 3.

Example problem

Using the hyperscaling relation α+2β+γ=dν\alpha + 2\beta + \gamma = d\nu, check the exact 2D Ising exponents α=0,β=1/8,γ=7/4,ν=1\alpha = 0, \beta = 1/8, \gamma = 7/4, \nu = 1 and the mean-field exponents α=0,β=1/2,γ=1,ν=1/2\alpha = 0, \beta = 1/2, \gamma = 1, \nu = 1/2 at d=4d = 4.

Solution

2D Ising: 0+2(1/8)+7/4=1/4+7/4=2=dν=2×10 + 2(1/8) + 7/4 = 1/4 + 7/4 = 2 = d\nu = 2 \times 1. ✓

Mean-field at d=4d = 4: 0+2(1/2)+1=2=dν=4×1/20 + 2(1/2) + 1 = 2 = d\nu = 4 \times 1/2. ✓ Both saturate the relation — the first because 2D Ising is below its upper critical dimension, the second because d=4d = 4 is the upper critical dimension of mean-field theory.

Types of transitions

TypeOrder parameterFree energyExample
Second-order (continuous)Continuous onset: φ ∝ |T−Tc|βSmooth, but derivatives divergeFerromagnet, liquid–gas (at critical point)
First-order (discontinuous)Jumps: φ finite at TcTwo minima coexist; latent heatWater freezing, solid–liquid
QuantumDriven by non-thermal parameter (pressure, field)At T = 0Superconductor–insulator, BCS gap
TopologicalNo local order parameter; invariant changesNo symmetry breakingAnyon statistics, topological insulator

The liquid–gas pressure–temperature phase diagram makes the geometry concrete: the coexistence curve separates two phases, ends at a critical point, and meets a second transition line at the triple point.

T P Solid Liquid Gas supercritical fluid triple point critical point coexistence curve (boiling) melting line sublimation
Details

Global idea. The type of transition is determined by the behavior of the free energy FF at the transition. If FF itself is continuous but a first derivative (entropy S=−∂F/∂TS = -\partial F/\partial T, volume V=∂F/∂PV = \partial F/\partial P) jumps, the transition is first-order: latent heat L=TcΔSL = T_c \Delta S is released, and two phases coexist in equilibrium. If the first derivatives are continuous but a second derivative (specific heat, compressibility, susceptibility) diverges, the transition is second-order (continuous).

Example. In the phase diagram above, follow the boiling curve from the triple point upward: crossing it at any pressure below the critical pressure is a first-order transition (water jumps from density ρg\rho_g to ρl\rho_l, absorbing latent heat). But the curve ends at the critical point — above it there is no distinction between liquid and gas, and approaching that endpoint along the curve gives a continuous transition with diverging compressibility and critical exponents. The melting line, by contrast, has no known critical point: solid and liquid differ by a symmetry (positional order), not just by degree.

Quantum and topological transitions. Quantum phase transitions occur at T=0T = 0 and are controlled by a non-thermal parameter; at finite TT they are rounded over a scale T∼ξ−zT \sim \xi^{-z} where zz is the dynamic critical exponent. Topological transitions (e.g., the Kosterlitz–Thouless transition in 2D superfluids) involve the unbinding of topological defects rather than symmetry breaking — the 2016 Nobel prize recognized this work.

Authoritative sources. Kosterlitz & Thouless, Ordering, metastability and phase transitions in two-dimensional systems (1973); Nobel Prize 2016: topological phase transitions.

Example problem

Water boils at 100 °C at 1 atm with latent heat L=2.26L = 2.26 MJ/kg. Explain why boiling is first-order while the critical-point endpoint of the boiling curve is second-order, and what happens to the latent heat as the critical point is approached.

Solution

At 100 °C, liquid and gas have different densities and entropies: SS (a first derivative of FF) jumps, so the transition is first-order and L=TΔS≠0L = T\Delta S \neq 0. As the critical point is approached along the coexistence curve, the density difference ρl−ρg\rho_l - \rho_g (the order parameter) vanishes as ∣T−Tc∣β|T - T_c|^{\beta}, so ΔS→0\Delta S \to 0 and the latent heat smoothly goes to zero — the discontinuity shrinks to nothing and only the diverging derivative (compressibility, specific heat) remains, which is the second-order signature.

Simulation: the 2D Ising model

The Ising model — spins si=±1s_i = \pm 1 on a lattice with energy E=−J∑⟨ij⟩sisjE = -J \sum_{\langle i j\rangle} s_i s_j — is the hydrogen atom of phase transitions. The simulation below runs Metropolis dynamics at temperature TT (in units of J/kBJ/k_B, with Tc=2.269T_c = 2.269): drag the slider through TcT_c and watch domains of aligned spins form, coarsen, and finally lock into a single ordered phase, while the magnetization ∣m∣|m| (the order parameter) turns on continuously.

Renormalization group

The RG is the systematic framework for critical phenomena. The idea: coarse-grain the system by integrating out short-distance degrees of freedom, then rescale. The Hamiltonian flows in parameter space:

Hℓ+1=R[Hℓ]H_{\ell+1} = \mathcal{R}[H_\ell]

Fixed points of the RG flow correspond to critical points. The relevant and irrelevant directions at a fixed point determine the critical exponents and the universality class.

critical surface (t = 0) critical fixed point ← relevant direction (reduced temperature t) → irrelevant directions flows off the critical surface run away (no criticality)
Details

Global idea. Coarse-graining replaces blocks of spins by a single effective spin and rescales lengths back. Repeating this generates a flow in the space of Hamiltonians. Criticality means scale invariance, and scale invariance means the flow must reach a fixed point (R[H∗]=H∗\mathcal{R}[H^*] = H^*): the system looks the same at every scale. Near a fixed point, perturbations decompose into relevant directions (they grow under coarse-graining and push the system away from criticality — the reduced temperature tt is one) and irrelevant directions (they shrink and die out). The rate of growth/decay is the critical exponent: ν=1/yt\nu = 1/y_t where yty_t is the relevant eigenvalue.

Example. An Ising ferromagnet, a liquid–gas system, and a lattice gas look nothing alike microscopically, yet all three flow to the same fixed point under coarse-graining — hence they share the exact same exponents (β=1/8\beta = 1/8 etc. in 2D). That is universality, and it is why studying one lattice model teaches us about real fluids.

Computing exponents. The epsilon expansion (ϵ=4−d\epsilon = 4 - d) computes exponents perturbatively around the upper critical dimension; the Monte Carlo RG and tensor network methods provide non-perturbative results in low dimensions.

Authoritative sources. Wilson & Kogut, The Renormalization Group and the ε Expansion (1974); Cardy, Scaling and Renormalization in Statistical Physics, ch. 3.

Example problem

A system sits exactly on its critical surface (t=0t = 0) but slightly off in an irrelevant direction with eigenvalue y<0y < 0. What happens to it under repeated coarse-graining, and why does this not spoil criticality?

Solution

The irrelevant perturbation shrinks as bynb^{y n} after nn coarse-graining steps of scale factor b>1b > 1, since y<0y < 0. The flow converges to the fixed point anyway, so the long-distance physics is still scale-invariant and critical. This is precisely why microscopic details do not affect critical exponents: they couple to irrelevant directions. Only a relevant direction (like t≠0t \neq 0) would grow and carry the flow away, cutting off the divergence of ξ\xi at a finite scale ξ∼∣t∣−ν\xi \sim |t|^{-\nu}.

Quiz

Quiz: correlation length near the critical point

A 2D Ising ferromagnet has critical temperature Tc=2.269 J/kBT_c = 2.269\, J/k_B. Estimate the correlation length at T=0.99 TcT = 0.99\, T_c using ν=1\nu = 1 (exact for 2D Ising). Then estimate it at T=0.999 TcT = 0.999\, T_c and explain what the trend means physically.

Solution

The correlation length scales as ξ∼∣T−Tc∣−ν\xi \sim |T - T_c|^{-\nu}. With ν=1\nu = 1:

∣T−Tc∣/Tc=0.01|T - T_c|/T_c = 0.01 so ξ∼(0.01)−1=100\xi \sim (0.01)^{-1} = 100 lattice spacings.

At T=0.999 TcT = 0.999\, T_c: ξ∼1000\xi \sim 1000 lattice spacings. Each factor of 10 closer to TcT_c multiplies the correlated region by 10. The system is approaching the critical point where ξ→∞\xi \to \infty and the correlation length becomes macroscopic: domains span the whole sample, fluctuations occur on all scales, and the power-law (scale-free) critical behavior takes over. You can watch this domain growth directly in the simulation above by sliding TT toward 2.2692.269.

How this fits together

Phase transitions are where the ensembles meet collective behavior. The grand canonical ensemble’s occupation function f(ε0)f(\varepsilon_0) diverges at the Bose condensation transition. The canonical ensemble’s free energy FF develops non-analyticities at critical points. The microcanonical ensemble’s entropy S(E)S(E) shows a double-hump for first-order transitions. The RG provides the unifying framework: all critical phenomena are fixed points of the coarse-graining flow. See GrandCanonicalEnsemble for the occupation-number picture, CanonicalEnsemble for the free-energy picture, and NonequilibriumStatisticalMechanics for driven and time-dependent transitions.