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(TTc)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ϕ4hϕ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ϕ3h=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(TcT)/b\phi = \pm\sqrt{-a/b} = \pm\sqrt{a_0(T_c - T)/b}

The order parameter grows as:

ϕ={0T>Tc±a0(TcT)/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 χ=ϕ/hh=0\chi = \partial\phi/\partial h|_{h=0} is obtained by expanding ϕ(h)=ϕ0+χh+\phi(h) = \phi_0 + \chi h + \cdots:

aϕ+bϕ3h=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(TTc)]\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(TcT)]\chi = 1/(-2a) = 1/[2a_0(T_c - T)]Giving γ"=1\gamma" = 1.

The free energy at equilibrium is:

feq={f0T>Tcf0a2/(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:

CTcCTc+=Tc2T2(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(T100)ϕ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±b24ac2c\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 TT^*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 TTc>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/(TTc)\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(T400)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(Tc400)=00.5(T_c - 400) = 0 giving Tc=400T_c = 400 K.

At T=300T = 300 K: a(300)=0.5(300400)=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(T400)]\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.