Skip to content

Statistical Mechanics

2.1 Microstates and Macrostates

A microstate is a complete specification of the state of a system (positions and momenta of all particles). A macrostate is specified by macroscopic variables (energy, volume, particle number).

The fundamental postulate of statistical mechanics states that for an isolated system in equilibrium, every accessible microstate is equally probable.

Definition. The multiplicity Ω(E,V,N)\Omega(E, V, N) is the number of microstates consistent with the macrostate (E,V,N)(E, V, N). The statistical entropy is

S=kBlnΩS = k_B \ln \Omega

Proposition 2.1. This definition of entropy agrees with the thermodynamic entropy: ΔS=kBln(Ωf/Ωi)\Delta S = k_B \ln(\Omega_f / \Omega_i).

2.2 The Boltzmann Distribution

Theorem 2.1 (Canonical Ensemble). For a system in thermal equilibrium with a heat bath at temperature TTThe probability of the system being in microstate ii with energy EiE_i is

Pi=1ZeEi/(kBT)P_i = \frac{1}{Z} e^{-E_i / (k_B T)}

Where the partition function is

Z=ieEi/(kBT)Z = \sum_{i} e^{-E_i / (k_B T)}

Proof. Consider the combined system (system + reservoir) with total energy EtotE_{\mathrm{tot}}. The probability of the system being in state ii is proportional to the number of reservoir microstates compatible with it, which is ΩR(EtotEi)\Omega_R(E_{\mathrm{tot} - E_i)}. Using SR=kBlnΩRS_R = k_B \ln \Omega_R:

P_i \propto \Omega_R(E_{\mathrm{tot} - E_i) = \exp\left(\frac{S_R(E_{\mathrm{tot} - E_i)}{k_B}\right)}}

Expanding SRS_R around EtotE_{\mathrm{tot}}: SR(EtotEi)SR(Etot)Ei(SRE)=SR(Etot)EiTS_R(E_{\mathrm{tot} - E_i) \approx S_R(E_{\mathrm{tot}) - E_i \left(\frac{\partial S_R}{\partial E}\right) = S_R(E_{\mathrm{tot}) - \frac{E_i}{T}}}}

Since (SR/E)=1/T(\partial S_R / \partial E) = 1/T. Therefore PieEi/(kBT)P_i \propto e^{-E_i / (k_BT)}And normalising gives the result. \blacksquare

2.3 Thermodynamic Quantities from the Partition Function

Theorem 2.2. The partition function determines all thermodynamic quantities:

E=lnZβ,F=kBTlnZ,S=EFT\langle E \rangle = -\frac{\partial \ln Z}{\partial \beta}, \quad F = -k_B T \ln Z, \quad S = \frac{\langle E \rangle - F}{T}

Where β=1/(kBT)\beta = 1/(k_BT).

Proof. E=iEiPi=1ZiEieβEi=1ZZβ=lnZβ\langle E \rangle = \sum_i E_i P_i = \frac{1}{Z}\sum_i E_i e^{-\beta E_i} = -\frac{1}{Z}\frac{\partial Z}{\partial \beta} = -\frac{\partial \ln Z}{\partial \beta}.

F=kBTlnZF = -k_BT \ln Z follows from F=UTS=ETSF = U - TS = \langle E \rangle - TS and the identification Z=eβFZ = e^{-\beta F}. \blacksquare

2.4 Ideal Gas

Theorem 2.3 (Partition Function of an Ideal Gas). For NN indistinguishable particles in a 3D box of volume VV:

Z_N = \frac{1}{N!}\left(\frac{V}{\lambda_{\mathrm{th}^3}\right)^N, \quad \lambda_{\mathrm{th} = \frac{h}{\sqrt{2\pi m k_B T}}}}

Where λth\lambda_{\mathrm{th}} is the thermal de Broglie wavelength.

Proof. The single-particle energy levels in a 3D box of side LL (V=L3V = L^3) are:

εnx,ny,nz=h28mL2(nx2+ny2+nz2)\varepsilon_{n_x, n_y, n_z} = \frac{h^2}{8mL^2}(n_x^2 + n_y^2 + n_z^2)

The single-particle partition function is:

z=nx,ny,nz=0eβh2(nx2+ny2+nz2)/(8mL2)(0eβh2n2/(8mL2)dn)3=(Lh2πmβ)3=Vλth3z = \sum_{n_x, n_y, n_z = 0}^{\infty} e^{-\beta h^2(n_x^2 + n_y^2 + n_z^2)/(8mL^2)} \approx \left(\int_0^{\infty} e^{-\beta h^2 n^2/(8mL^2)} dn\right)^3 = \left(\frac{L}{h}\sqrt{\frac{2\pi m}{\beta}}\right)^3 = \frac{V}{\lambda_{\mathrm{th}^3}}

For NN indistinguishable particles (correct Boltzmann counting): ZN=zN/N!Z_N = z^N / N!. \blacksquare

Corollary 2.4. From ZNZ_NWe recover the ideal gas law:

F = -k_BT \ln Z_N = -k_BT\left[N\ln\left(\frac{V}{\lambda_{\mathrm{th}^3}\right) - \ln N!\right]}

P=(FV)T=NkBTVP = -\left(\frac{\partial F}{\partial V}\right)_T = \frac{Nk_BT}{V}

Giving PV=NkBTPV = Nk_BT.

2.5 The Equipartition Theorem

Theorem 2.5 (Equipartition). For a classical system in thermal equilibrium, each quadratic degree of freedom in the Hamiltonian contributes 12kBT\frac{1}{2}k_BT to the average energy.

Proof. Consider a single degree of freedom with Hamiltonian H=ap2H = ap^2 (or bq2bq^2). The average energy is:

ap2=ap2eβap2dpeβap2dp=βln(eβap2dp)=βln(πaβ)=12β=kBT2\langle ap^2 \rangle = \frac{\int_{-\infty}^{\infty} ap^2 e^{-\beta ap^2}\, dp}{\int_{-\infty}^{\infty} e^{-\beta ap^2}\, dp} = -\frac{\partial}{\partial \beta}\ln\left(\int_{-\infty}^{\infty} e^{-\beta ap^2}\, dp\right) = -\frac{\partial}{\partial \beta}\ln\left(\sqrt{\frac{\pi}{a\beta}}\right) = \frac{1}{2\beta} = \frac{k_BT}{2}

The same calculation for bq2bq^2 gives another kBT/2k_BT/2. \blacksquare

Application. A monatomic ideal gas has 3 translational degrees of freedom: U=32NkBTU = \frac{3}{2}Nk_BT and CV=32NkBC_V = \frac{3}{2}Nk_B. A diatomic gas also has 2 rotational degrees of freedom: U=52NkBTU = \frac{5}{2}Nk_BT and CV=52NkBC_V = \frac{5}{2}Nk_B (at temperatures where vibration is frozen out).

2.6 Quantum Statistical Distributions

Fermi—Dirac Statistics (for fermions, particles with half-integer spin):

ni=1e(Eiμ)/(kBT)+1\langle n_i \rangle = \frac{1}{e^{(E_i - \mu)/(k_BT)} + 1}

Where μ\mu is the chemical potential.

Bose—Einstein Statistics (for bosons, particles with integer spin):

ni=1e(Eiμ)/(kBT)1\langle n_i \rangle = \frac{1}{e^{(E_i - \mu)/(k_BT)} - 1}

Maxwell—Boltzmann Statistics (classical limit, μ\mu very negative):

ni=e(Eiμ)/(kBT)\langle n_i \rangle = e^{-(E_i - \mu)/(k_BT)}

The classical limit applies when the thermal de Broglie wavelength is much smaller than the inter-particle spacing: λth3V/N\lambda_{\mathrm{th}^3 \ll V/N}.

2.7 The Fermi Gas

Definition. The Fermi energy εF\varepsilon_F is the chemical potential at T=0T = 0:

εF=22m(3π2n)2/3\varepsilon_F = \frac{\hbar^2}{2m}\left(3\pi^2 n\right)^{2/3}

Where n=N/Vn = N/V is the particle number density.

Proposition 2.6. At T=0T = 0All states with EεFE \leq \varepsilon_F are occupied and all states with E>εFE > \varepsilon_F are empty. The ground-state energy of a 3D Fermi gas is:

U0=35NεFU_0 = \frac{3}{5}N\varepsilon_F

Proof. U0=0εFEg(E)dEU_0 = \int_0^{\varepsilon_F} E \cdot g(E)\, dE where g(E)=V2π2(2m2)3/2Eg(E) = \frac{V}{2\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2}\sqrt{E} is the density of states. Evaluating: U0=V2π2(2m2)3/225εF5/2=35NεFU_0 = \frac{V}{2\pi^2}\left(\frac{2m}{\hbar^2}\right)^{3/2} \cdot \frac{2}{5}\varepsilon_F^{5/2} = \frac{3}{5}N\varepsilon_F. \blacksquare

2.8 Blackbody Radiation

Planck”s Law gives the spectral energy density of blackbody radiation:

u(ν,T)=8πhν3c31ehν/(kBT)1u(\nu, T) = \frac{8\pi h \nu^3}{c^3} \cdot \frac{1}{e^{h\nu/(k_BT)} - 1}

Stefan—Boltzmann Law: The total radiated power per unit area:

j=σT4,σ=π2kB4603c2j = \sigma T^4, \quad \sigma = \frac{\pi^2 k_B^4}{60 \hbar^3 c^2}

Wien’s Displacement Law: The peak frequency satisfies νmax/T=const\nu_{\mathrm{max} / T = \mathrm{const}}.

2.9 Worked Examples

Problem. Calculate the Fermi energy and Fermi temperature for copper. Given: electron density n8.5×1028m3n \approx 8.5 \times 10^{28}\,\mathrm{m}^{-3}, me=9.109×1031m_e = 9.109 \times 10^{-31} kg.

Solution

εF=22me(3π2n)2/3\varepsilon_F = \frac{\hbar^2}{2m_e}(3\pi^2 n)^{2/3}

=(1.055×1034)22×9.109×1031×(3π2×8.5×1028)2/3= \frac{(1.055 \times 10^{-34})^2}{2 \times 9.109 \times 10^{-31}} \times (3\pi^2 \times 8.5 \times 10^{28})^{2/3}

(3π2×8.5×1028)1/3=(2.52×1030)1/31.36×1010(3\pi^2 \times 8.5 \times 10^{28})^{1/3} = (2.52 \times 10^{30})^{1/3} \approx 1.36 \times 10^{10}

(3π2n)2/3=(1.36×1010)2=1.85×1020(3\pi^2 n)^{2/3} = (1.36 \times 10^{10})^2 = 1.85 \times 10^{20}

εF=1.113×10681.822×1030×1.85×10201.13×1018J7.0eV\varepsilon_F = \frac{1.113 \times 10^{-68}}{1.822 \times 10^{-30}} \times 1.85 \times 10^{20} \approx 1.13 \times 10^{-18}\,\mathrm{J} \approx 7.0\,\mathrm{eV}

TF=εF/kB=1.13×1018/1.381×102381800KT_F = \varepsilon_F / k_B = 1.13 \times 10^{-18} / 1.381 \times 10^{-23} \approx 81800\,\mathrm{K}

The Fermi temperature is much larger than room temperature, confirming that copper electrons are in the degenerate regime. \blacksquare

Worked Example: Entropy of Mixing

Solution. Two ideal gases of NN particles each, initially separated by a partition, are allowed to mix. Calculate the entropy change.

Before mixing: the total entropy is 2×NkB[ln(VNλ3)+52]2 \times Nk_B\left[\ln\left(\frac{V}{N\lambda^3}\right) + \frac{5}{2}\right] (for a monatomic gas).

After mixing: each gas occupies volume 2V2VSo the total entropy is:

Sf=2×NkB[ln(2VNλ3)+52]S_f = 2 \times Nk_B\left[\ln\left(\frac{2V}{N\lambda^3}\right) + \frac{5}{2}\right]

ΔSmix=SfSi=2NkBln(2VNλ3)2NkBln(VNλ3)=2NkBln2\Delta S_{\mathrm{mix} = S_f - S_i = 2Nk_B\ln\left(\frac{2V}{N\lambda^3}\right) - 2Nk_B\ln\left(\frac{V}{N\lambda^3}\right) = 2Nk_B\ln 2}

For 1 mole of each gas: ΔSmix=2Rln211.5J/K\Delta S_{\mathrm{mix} = 2R\ln 2 \approx 11.5\,\mathrm{J}/K}.

Gibbs paradox. If the two gases are identical, the entropy of mixing is zero (no physical change). The resolution is that identical particles are indistinguishable, and the correct counting already accounts for this via the 1/N!1/N! factor in the partition function. \blacksquare

2.10 Common Pitfalls

  • The classical limit does not always apply. When λth3V/N\lambda_{\mathrm{th}^3 \gtrsim V/N}Quantum …/4-statistics-and-probability/2_statistics (Fermi-Dirac or Bose-Einstein) must be used. This is critical for electrons in metals and for helium-4 at low temperatures.
  • The Boltzmann distribution applies to systems in contact with a heat bath, not isolated systems. For isolated systems, use the microcanonical ensemble (all accessible microstates equally probable).
  • The partition function must account for indistinguishability. The 1/N!1/N! factor in ZNZ_N is essential for obtaining the correct entropy (otherwise the entropy is not extensive and the Gibbs paradox arises).