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+δEd3Nq d3NpN!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=∂S∂E\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):

⟨∑ipi⋅∂H∂pi⟩=3NkBT\left\langle \sum_i \mathbf{p}_i \cdot \frac{\partial \mathcal{H}}{\partial \mathbf{p}_i} \right\rangle = 3Nk_B T

⟨∑iri⋅∂H∂ri⟩=−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 U∝rnU \propto r^nThis gives:

⟨K⟩=n2⟨U⟩\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 (N→∞N \to \infty, V→∞V \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/E∼1/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(β)=−β−1ln⁡Z(β)F(\beta) = -\beta^{-1} \ln Z(\beta) are related by the Legendre transform:

F(β)=inf⁡E[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 T→0T \to 0, the entropy of a system approaches a constant (typically zero for a non-degenerate ground state):

lim⁡T→0S(E,V,N)=kBln⁡g0\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=(N↑−N↓)μB=(2N↑−N)μBE = (N_\uparrow - N_\downarrow)\mu B = (2N_\uparrow - N)\mu B, the number of microstates with N↑N_\uparrow up-spins is:

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

Using Stirling’s approximation:

S=kBln⁡Ω=kB[Nln⁡N−N↑ln⁡N↑−(N−N↑)ln⁡(N−N↑)]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μBln⁡N−N↑N↑\frac{1}{T} = \frac{k_B}{2\mu B} \ln\frac{N - N_\uparrow}{N_\uparrow}

Solving for the magnetisation M=(N↑−N↓)μ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)=EN−1(N−1)! (ℏω)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+N−1N−1)=(n+N−1)!n!(N−1)!\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)−nln⁡n−Nln⁡N]S = k_B\left[(n+N)\ln(n+N) - n\ln n - N\ln N\right]

1T=∂S∂E=kBℏω[ln⁡(n+N)−ln⁡n]=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 (n≫Nn \gg N): E≈NkBTE \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)NEN−1(N−1)! δ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⟩/⟨E⟩∼1/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+δEd3Nq d3Np/(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=T ∂S/∂VP = T\,\partial S/\partial V
Chemical potentialμ=−T ∂S/∂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}]