Skip to content

Landau Theory of Phase Transitions

Landau theory provides a phenomenological framework for continuous (second-order) phase transitions by expanding the free energy in powers of an order parameter ϕ\phi.

The Landau free energy density (in the absence of external fields) is:

f(ϕ,T)=f0(T)+12a(T)ϕ2+14bϕ4+16cϕ6+⋯f(\phi, T) = f_0(T) + \frac{1}{2}a(T)\phi^2 + \frac{1}{4}b\phi^4 + \frac{1}{6}c\phi^6 + \cdots

Assumptions:

  • ff is analytic in ϕ\phi near the transition
  • Symmetry ϕ→−ϕ\phi \to -\phi (e.g., Ising systems) eliminates odd powers
  • b>0b > 0 for stability
  • a(T)=a0(T−Tc)a(T) = a_0(T - T_c) changes sign at TcT_c

With an external field hh conjugate to ϕ\phiAdd −hϕ-h\phi:

f(ϕ,T)=f0+12a(T)ϕ2+14bϕ4−hϕf(\phi, T) = f_0 + \frac{1}{2}a(T)\phi^2 + \frac{1}{4}b\phi^4 - h\phi

The equilibrium order parameter minimizes ff:

∂f∂ϕ=aϕ+bϕ3−h=0\frac{\partial f}{\partial \phi} = a\phi + b\phi^3 - h = 0

For h=0h = 0:

  • T>TcT > T_c (a>0a > 0): minimum at ϕ=0\phi = 0 (disordered phase)
  • T<TcT < T_c (a<0a < 0): minima at ϕ=±−a/b=±a0(Tc−T)/b\phi = \pm\sqrt{-a/b} = \pm\sqrt{a_0(T_c - T)/b}

The order parameter grows as:

ϕ={0T>Tc±a0(Tc−T)/bT<Tc\phi = \begin{cases} 0 & T > T_c \\ \pm\sqrt{a_0(T_c - T)/b} & T < T_c \end{cases}

This yields the mean-field critical exponent β=1/2\beta = 1/2.

The susceptibility χ=∂ϕ/∂h∣h=0\chi = \partial\phi/\partial h|_{h=0} is obtained by expanding ϕ(h)=ϕ0+χh+⋯\phi(h) = \phi_0 + \chi h + \cdots:

aϕ+bϕ3−h=0  ⟹  (a+3bϕ02)χ=1a\phi + b\phi^3 - h = 0 \implies (a + 3b\phi_0^2)\chi = 1

  • T>TcT > T_c: ϕ0=0\phi_0 = 0So χ=1/a=1/[a0(T−Tc)]\chi = 1/a = 1/[a_0(T - T_c)]Giving γ=1\gamma = 1.
  • T<TcT < T_c: ϕ02=−a/b\phi_0^2 = -a/bSo χ=1/(−2a)=1/[2a0(Tc−T)]\chi = 1/(-2a) = 1/[2a_0(T_c - T)]Giving γ"=1\gamma" = 1.

The free energy at equilibrium is:

feq={f0T>Tcf0−a2/(4b)T<Tcf_{\text{eq} = \begin{cases} f_0 & T > T_c \\ f_0 - a^2/(4b) & T < T_c \end{cases}}

The specific heat discontinuity is:

CTc−−CTc+=−Tc∂2∂T2(−a24b)∣Tc=Tca022bC_{T_c^-} - C_{T_c^+} = -T_c \frac{\partial^2}{\partial T^2}\left(\frac{-a^2}{4b}\right)\bigg|_{T_c} = \frac{T_c a_0^2}{2b}

This is a finite jump (α=0\alpha = 0 in mean-field theory).

Worked Example 11.1: Landau Free Energy Minimum

Consider f=12(T−100)ϕ2+14ϕ4f = \frac{1}{2}(T - 100)\phi^2 + \frac{1}{4}\phi^4 (in arbitrary units where a0=b=1a_0 = b = 1).

At T=50T = 50 (a=−50a = -50): f=−25ϕ2+14ϕ4f = -25\phi^2 + \frac{1}{4}\phi^4.

∂f∂ϕ=−50ϕ+ϕ3=0  ⟹  ϕ=0 (max)orϕ=±50=±7.07 (min)\frac{\partial f}{\partial \phi} = -50\phi + \phi^3 = 0 \implies \phi = 0 \text{ (max)} or \phi = \pm\sqrt{50} = \pm 7.07 \text{ (min)}

fmin=−25(50)+14(2500)=−1250+625=−625f_{\text{min} = -25(50) + \frac{1}{4}(2500) = -1250 + 625 = -625}

At T=150T = 150 (a=50a = 50): f=25ϕ2+14ϕ4f = 25\phi^2 + \frac{1}{4}\phi^4.

∂f∂ϕ=50ϕ+ϕ3=0  ⟹  ϕ=0 (min)\frac{\partial f}{\partial \phi} = 50\phi + \phi^3 = 0 \implies \phi = 0 \text{ (min)}

fmin=0f_{\text{min} = 0}

The free energy drops by 625 units when going below Tc=100T_c = 100Driving the transition.

Worked Example 11.2: First-Order Transition in Landau Theory

When b<0b < 0 (which can happen in systems with first-order transitions), we must include the ϕ6\phi^6 term with c>0c > 0:

f=12a(T)ϕ2+14bϕ4+16cϕ6f = \frac{1}{2}a(T)\phi^2 + \frac{1}{4}b\phi^4 + \frac{1}{6}c\phi^6

The equilibrium condition ∂f/∂ϕ=0\partial f/\partial \phi = 0 gives:

ϕ(a+bϕ2+cϕ4)=0\phi(a + b\phi^2 + c\phi^4) = 0

The quartic factor has solutions when:

ϕ2=−b±b2−4ac2c\phi^2 = \frac{-b \pm \sqrt{b^2 - 4ac}}{2c}

This requires b2>4acb^2 > 4acWhich occurs when TT is below some temperature T∗>TcT^* > T_c. Between TcT_c and T∗T^*The system undergoes a first-order transition because the order parameter jumps discontinuously from zero to a finite value.

  • Critical exponents (mean-field): β=1/2\beta = 1/2 (order parameter), γ=1\gamma = 1 (susceptibility), α=0\alpha = 0 (specific heat jump), δ=3\delta = 3 (critical isotherm: h∝ϕ3h \propto \phi^3 at T=TcT = T_c).
  • Universality: Systems with the same symmetry and dimensionality share the same critical exponents, regardless of microscopic details. Landau theory gives mean-field exponents, which are exact only above the upper critical dimension (d>4d > 4 for short-range interactions).
  • Clausius-Clapeyron analogue: At a first-order transition (when b<0b < 0), the discontinuity in the order parameter gives a latent heat L=TcΔsL = T_c \Delta s where Δs=−∂f/∂T\Delta s = -\partial f/\partial T is the entropy jump.
  • Ginzburg criterion: Mean-field theory is valid when fluctuations are small, i.e., when ∣T−Tc∣>Tc(a02kB2)/(32π2b2ξ0d)|T - T_c| > T_c(a_0^2 k_B^2)/(32\pi^2 b^2 \xi_0^d). For d<4d < 4, this fails very close to TcT_c.
  • Assuming Landau theory is always valid: It is a mean-field theory. Near TcT_c in low dimensions, critical fluctuations dominate and renormalisation group methods are required.
  • Forgetting that bb can be negative: If b<0b < 0, the ϕ6\phi^6 term must be included to ensure stability. The transition becomes first-order, and the simple ϕ=±−a/b\phi = \pm\sqrt{-a/b} solution does not apply.
  • Neglecting the role of symmetry: The form of the Landau expansion depends on the symmetry of the order parameter. A vector order parameter (e.g., in the XY model) requires a different expansion than a scalar.
  • Confusing the order parameter with a physical observable: The order parameter ϕ\phi is an abstract quantity. For a ferromagnet it is the magnetisation; for a superfluid it is the condensate wavefunction; for a liquid-gas transition it is the density difference.
  • Ferromagnetic transitions: The Landau theory with ϕ=M\phi = M (magnetisation) predicts the Curie temperature and the Curie-Weiss law χ∝1/(T−Tc)\chi \propto 1/(T - T_c) for the susceptibility above TcT_c.
  • Superfluid helium: The order parameter is the complex condensate wavefunction ψ\psi. The Landau-Ginzburg expansion includes ∣ψ∣2|\psi|^2 terms and gradient terms, leading to the Ginzburg-Landau theory of superconductivity.
  • Binary alloys: The order parameter describes the degree of chemical ordering (e.g., Cu-Zn ordering in brass). The Landau theory predicts the order-disorder transition temperature.
  • Liquid crystals: Nematic-isotropic transitions can be described by a tensor order parameter QijQ_{ij}. The Landau expansion includes both scalar and tensor invariants.

11.8 Worked Example: Finding the Transition Temperature

Section titled “11.8 Worked Example: Finding the Transition Temperature”

A magnetic system has Landau coefficients a(T)=0.5(T−400)a(T) = 0.5(T - 400) K and b=2.0b = 2.0 (arbitrary units). Find TcT_c and the magnetisation at T=300T = 300 K.

At TcT_c, the coefficient a(Tc)=0a(T_c) = 0, so 0.5(Tc−400)=00.5(T_c - 400) = 0 giving Tc=400T_c = 400 K.

At T=300T = 300 K: a(300)=0.5(300−400)=−50a(300) = 0.5(300 - 400) = -50. The equilibrium magnetisation is:

M=−a/b=50/2.0=25=5.0M = \sqrt{-a/b} = \sqrt{50/2.0} = \sqrt{25} = 5.0

The free energy at equilibrium: feq=−a2/(4b)=−2500/(8)=−312.5f_{\text{eq}} = -a^2/(4b) = -2500/(8) = -312.5.

The susceptibility above TcT_c: χ=1/a=1/[0.5(T−400)]\chi = 1/a = 1/[0.5(T - 400)]. At T=500T = 500 K, χ=1/50=0.02\chi = 1/50 = 0.02.

The specific heat jump at TcT_c: ΔC=Tca02/(2b)=400×0.25/(4)=25\Delta C = T_c a_0^2/(2b) = 400 \times 0.25/(4) = 25.