Skip to content

Microcanonical Ensemble

The microcanonical ensemble describes an isolated system with fixed total energy EEParticle number NNAnd volume VV.

The number of microstates with energy between EE and E+δEE + \delta E is:

Ω(E,V,N)=E<H<E+δEd3Nqd3NpN!h3N\Omega(E, V, N) = \int_{E < \mathcal{H} < E + \delta E} \frac{d^{3N}q\, d^{3N}p}{N!h^{3N}}

The entropy (Boltzmann entropy):

S(E,V,N)=kBlnΩ(E,V,N)S(E, V, N) = k_B \ln \Omega(E, V, N)

The temperature is defined via:

1T=SE\frac{1}{T} = \frac{\partial S}{\partial E}

14.2 The Ideal Gas in the Microcanonical Ensemble

Section titled “14.2 The Ideal Gas in the Microcanonical Ensemble”

For NN non-interacting particles in volume VV with total energy EE:

Ω=VNN!(2πmE)3N/2EΓ(3N/2)h3NδEE\Omega = \frac{V^N}{N!}\frac{(2\pi m E)^{3N/2}}{E\, \Gamma(3N/2)\, h^{3N}} \cdot \frac{\delta E}{E}

Using Stirling”s approximation and the large-argument expansion of the Gamma function:

S=NkB[ln ⁣(VN)+32ln ⁣(4πmE3Nh2)+52]S = Nk_B\left[\ln\!\left(\frac{V}{N}\right) + \frac{3}{2}\ln\!\left(\frac{4\pi m E}{3Nh^2}\right) + \frac{5}{2}\right]

This is the Sackur—Tetrode equation, identical to the canonical ensemble result (as expected by ensemble equivalence).

From 1/T=S/E1/T = \partial S/\partial E:

E=32NkBTE = \frac{3}{2}Nk_B T

Reproducing the equipartition theorem.

For a system with Hamiltonian H=ipi2/(2mi)+U(r1,,rN)\mathcal{H} = \sum_i p_i^2/(2m_i) + U(\mathbf{r}_1, \ldots, \mathbf{r}_N):

ipiHpi=3NkBT\left\langle \sum_i \mathbf{p}_i \cdot \frac{\partial \mathcal{H}}{\partial \mathbf{p}_i} \right\rangle = 3Nk_B T

iriHri=3NkBT\left\langle \sum_i \mathbf{r}_i \cdot \frac{\partial \mathcal{H}}{\partial \mathbf{r}_i} \right\rangle = -3Nk_B T

For a power-law potential UrnU \propto r^nThis gives:

K=n2U\langle K \rangle = \frac{n}{2}\langle U \rangle

(For the harmonic oscillator, n=2n = 2: K=U\langle K \rangle = \langle U \rangle.)

14.4 Equivalence of Ensembles in the Thermodynamic Limit

Section titled “14.4 Equivalence of Ensembles in the Thermodynamic Limit”

In the thermodynamic limit (NN \to \infty, VV \to \infty, N/VN/V fixed), the microcanonical, canonical, and grand canonical ensembles produce identical thermodynamic predictions. This is a consequence of the fact that the energy fluctuations in the canonical ensemble scale as ΔE/E1/N\Delta E / E \sim 1/\sqrt{N}, vanishing in the limit.

Proposition 14.1. For a system with Hamiltonian H\mathcal{H}, the microcanonical entropy S(E)S(E) and the canonical free energy F(β)=β1lnZ(β)F(\beta) = -\beta^{-1} \ln Z(\beta) are related by the Legendre transform:

F(β)=infE[Eβ1S(E)]F(\beta) = \inf_E [E - \beta^{-1} S(E)]

14.5 The Third Law of Thermodynamics from the Microcanonical Ensemble

Section titled “14.5 The Third Law of Thermodynamics from the Microcanonical Ensemble”

Proposition 14.2 (Nernst’s Theorem). As T0T \to 0, the entropy of a system approaches a constant (typically zero for a non-degenerate ground state):

limT0S(E,V,N)=kBlng0\lim_{T \to 0} S(E, V, N) = k_B \ln g_0

where g0g_0 is the degeneracy of the ground state.

In the microcanonical picture, at T=0T = 0 the system occupies only the ground state microstate(s). If the ground state is unique, Ω=1\Omega = 1 and S=0S = 0.

Problem. Consider NN non-interacting spin-1/2 particles in a magnetic field BB. Each spin has energy ±μB\pm \mu B. Find the microcanonical entropy and the equation of state.

Solution

For a system with total energy E=(NN)μB=(2NN)μBE = (N_\uparrow - N_\downarrow)\mu B = (2N_\uparrow - N)\mu B, the number of microstates with NN_\uparrow up-spins is:

Ω(N)=(NN)=N!N!(NN)!\Omega(N_\uparrow) = \binom{N}{N_\uparrow} = \frac{N!}{N_\uparrow! (N - N_\uparrow)!}

Using Stirling’s approximation:

S=kBlnΩ=kB[NlnNNlnN(NN)ln(NN)]S = k_B \ln \Omega = k_B[N\ln N - N_\uparrow\ln N_\uparrow - (N - N_\uparrow)\ln(N - N_\uparrow)]

From 1/T=S/E1/T = \partial S/\partial E:

1T=kB2μBlnNNN\frac{1}{T} = \frac{k_B}{2\mu B} \ln\frac{N - N_\uparrow}{N_\uparrow}

Solving for the magnetisation M=(NN)μBM = (N_\uparrow - N_\downarrow)\mu B:

M=Nμtanh(μBkBT)M = N\mu \tanh\left(\frac{\mu B}{k_B T}\right)

This is the Brillouin function for spin-1/2, matching the canonical ensemble prediction.

Worked Example 14.2: Density of States for $N$ Harmonic Oscillators

For NN independent harmonic oscillators with frequency ω\omegaTotal energy EE:

Ω(E)=EN1(N1)!(ω)N\Omega(E) = \frac{E^{N-1}}{(N-1)!\,(\hbar\omega)^N}

Proof: The number of ways to distribute E/(ω)E/(\hbar\omega) energy quanta among NN oscillators is the stars-and-bars problem:

Ω=(n+N1N1)=(n+N1)!n!(N1)!\Omega = \binom{n + N - 1}{N - 1} = \frac{(n+N-1)!}{n!(N-1)!}

Where n=E/(ω)n = E/(\hbar\omega). For large nn using Stirling’s approximation:

S=kB[(n+N)ln(n+N)nlnnNlnN]S = k_B\left[(n+N)\ln(n+N) - n\ln n - N\ln N\right]

1T=SE=kBω[ln(n+N)lnn]=kBωln ⁣(1+Nn)\frac{1}{T} = \frac{\partial S}{\partial E} = \frac{k_B}{\hbar\omega}\left[\ln(n+N) - \ln n\right] = \frac{k_B}{\hbar\omega}\ln\!\left(1 + \frac{N}{n}\right)

At high TT (nNn \gg N): ENkBTE \approx Nk_B T (equipartition, each oscillator has energy kBTk_B T).

14.7 Worked Example: Ideal Gas in Two Dimensions

Section titled “14.7 Worked Example: Ideal Gas in Two Dimensions”

Problem. Find the microcanonical entropy of an ideal gas confined to a two-dimensional area AA with NN particles and total energy EE.

Solution

The phase space volume for NN particles in 2D with energy less than EE is:

Σ(E)=ANN!h2N(2πmE)NΓ(N+1)\Sigma(E) = \frac{A^N}{N! h^{2N}} \cdot \frac{(2\pi m E)^N}{\Gamma(N+1)}

The number of states with energy between EE and E+δEE + \delta E is Ω=(Σ/E)δE\Omega = (\partial\Sigma/\partial E)\,\delta E:

Ω=ANN!h2N(2πm)NEN1(N1)!δE\Omega = \frac{A^N}{N! h^{2N}} \cdot \frac{(2\pi m)^N E^{N-1}}{(N-1)!}\,\delta E

Using Stirling’s approximation:

S=NkB[ln ⁣(AN)+ln ⁣(2πmEh2)+2]S = Nk_B\left[\ln\!\left(\frac{A}{N}\right) + \ln\!\left(\frac{2\pi m E}{h^2}\right) + 2\right]

The equation of state is PA=NkBTPA = Nk_B T (the 2D analogue of PV=NkBTPV = Nk_B T), and the internal energy is E=NkBTE = Nk_B T.

\blacksquare

In the canonical ensemble, the energy fluctuates around its mean value E\langle E \rangle. The variance is related to the heat capacity:

(ΔE)2=kBT2CV\langle (\Delta E)^2 \rangle = k_B T^2 C_V

The relative fluctuation (ΔE)2/E1/N\sqrt{\langle (\Delta E)^2 \rangle}/\langle E \rangle \sim 1/\sqrt{N}, vanishing in the thermodynamic limit. This justifies the equivalence of microcanonical and canonical ensembles for macroscopic systems.

QuantityExpression
Microcanonical partition functionΩ(E,V,N)=E<H<E+δEd3Nqd3Np/(N!h3N)\Omega(E, V, N) = \int_{E < \mathcal{H} < E + \delta E} d^{3N}q\,d^{3N}p / (N! h^{3N})
Boltzmann entropyS=kBlnΩS = k_B \ln \Omega
Temperature1/T=S/E1/T = \partial S/\partial E
PressureP=TS/VP = T\,\partial S/\partial V
Chemical potentialμ=TS/N\mu = -T\,\partial S/\partial N
Sackur-Tetrode (ideal gas)S=NkB[ln(V/N)+32ln(4πmE/(3Nh2))+52]S = Nk_B[\ln(V/N) + \frac{3}{2}\ln(4\pi m E/(3Nh^2)) + \frac{5}{2}]