Skip to content

Common Pitfalls

  • Confusing the microcanonical, canonical, and grand canonical ensembles. The microcanonical ensemble describes an isolated system with fixed E,V,NE, V, N. The canonical ensemble describes a system in contact with a heat bath at fixed T,V,NT, V, N. The grand canonical ensemble describes a system exchanging both energy and particles, at fixed μ,V,T\mu, V, T.

  • Forgetting the 1/N!1/N! for indistinguishable particles. Without this factor, entropy is not extensive and the Gibbs paradox arises. This is essential for all quantum statistical mechanics.

  • Applying the equipartition theorem to quantum systems. At temperatures below the characteristic energy spacing (kBTΔEk_BT \ll \Delta E), the relevant degrees of freedom are “frozen out” and do not contribute to CVC_V.

  • Assuming the classical limit always applies. Electrons in metals are degenerate (TTFT \ll T_F) and must be treated with Fermi-Dirac statistics. Helium-4 at low temperatures exhibits Bose-Einstein condensation and superfluidity. The classical limit nλth31n\lambda_{\mathrm{th}}^3 \ll 1 is violated in these cases.

  • Confusing μ=0\mu = 0 for bosons with μ\mu for fermions. For bosons, με0\mu \leq \varepsilon_0 and μ0\mu \to 0 at BEC. For fermions, μεF\mu \approx \varepsilon_F at low temperatures and can be much larger than ε0\varepsilon_0.

  • Using mean field critical exponents in 2D. Mean field theory gives β=1/2\beta = 1/2 everywhere, but the exact 2D Ising result is β=1/8\beta = 1/8. Mean field theory is qualitatively wrong in low dimensions.

Ensemble selection. Identify which quantities are fixed in the physical setup: isolated system \to microcanonical; in contact with heat bath \to canonical; open to particle exchange \to grand canonical. Always specify the control parameters before choosing the ensemble.

The 1/N!1/N! factor. For NN indistinguishable particles, the partition function is ZN=zN/N!Z_N = z^N/N! where z=Z1z = Z_1 is the single-particle partition function. This ensures the entropy S=kBlnZNS = k_B \ln Z_N scales linearly with NN and resolves the Gibbs paradox.

Equipartition in the classical limit. Equipartition holds only when kBTk_B T is much larger than the spacing between energy levels. For a harmonic oscillator: CV=kBC_V = k_B at high TT, but CV0C_V \to 0 as T0T \to 0 due to quantization. The Einstein model captures this for lattice vibrations.

Recognizing degeneracy. The condition nλth31n\lambda_{\mathrm{th}}^3 \ll 1 must be checked before applying Maxwell-Boltzmann statistics. Here λth=h/2πmkBT\lambda_{\mathrm{th}} = h/\sqrt{2\pi m k_B T} is the thermal de Broglie wavelength. If this condition fails, quantum statistics are required.

Problem 1 (Ensemble Selection). A container of gas is placed in a heat bath at temperature TT, has fixed volume VV, but the walls are permeable to particles. Which ensemble should be used?

Solution. Since TT and μ\mu are fixed (particle exchange with bath), and VV is fixed, use the grand canonical ensemble. The grand partition function Ξ=NeβμNZN\Xi = \sum_N e^{\beta\mu N} Z_N describes the system, where ZNZ_N is the canonical partition function. \blacksquare

Problem 2 (Gibbs Paradox). Two identical ideal gases at the same temperature and pressure are separated by a partition. The partition is removed. Compute the entropy change with and without the 1/N!1/N! factor.

Solution. Without the 1/N!1/N! factor: Si=2NkB(lnV+const)S_i = 2Nk_B(\ln V + \text{const}), Sf=2NkB(ln(2V)+const)S_f = 2Nk_B(\ln(2V) + \text{const}) so ΔS=2NkBln2\Delta S = 2Nk_B\ln 2. But the gases are identical, so mixing should produce no entropy change. With the 1/N!1/N! factor: the factor of N!N! in the denominator cancels this spurious increase, giving ΔS=0\Delta S = 0. \blacksquare

Problem 3 (Low-Temperature Heat Capacity). Estimate the heat capacity of electrons in a metal at TTFT \ll T_F, where TF=εF/kBT_F = \varepsilon_F/k_B is the Fermi temperature.

Solution. Only electrons within kBTk_B T of the Fermi surface are excited. The fraction of excited electrons is T/TF\sim T/T_F, and each gains energy kBT\sim k_B T. So the electronic heat capacity is CVγTC_V \approx \gamma T where γ=(π2/2)(kB2/εF)\gamma = (\pi^2/2)(k_B^2/\varepsilon_F). This linear TT dependence is a key signature of a Fermi liquid. \blacksquare

  • Confusing heat and temperature. Heat QQ is energy transferred due to temperature difference; temperature TT is the thermodynamic potential that determines the direction of heat flow. A system can have high temperature but low heat content (e.g., a spark at 10001000 K has little energy).

  • Misapplying CVC_V and CPC_P. CVC_V applies at constant volume, CPC_P at constant pressure. For an ideal gas, CPCV=nRC_P - C_V = nR. For condensed phases, the difference is much smaller. Using the wrong one leads to errors in enthalpy and entropy calculations.

  • Forgetting the Maxwell relations. The four Maxwell relations are derived from the equality of mixed partial derivatives of thermodynamic potentials. Neglecting them can lead to inconsistent thermodynamic cycles.

  • Assuming reversible processes everywhere. Real processes have irreversibilities (friction, uncontrolled expansion, mixing). Entropy is not conserved in irreversible processes. Always compute dSδQ/TdS \geq \delta Q/T.

  • Phase transition classification errors. First-order transitions have latent heat and discontinuous order parameter. Second-order transitions have continuous order parameter but divergent susceptibility. The order of transition must be determined from the free energy, not from intuition.

PotentialSymbolNatural VariablesDifferential
Internal energyUUS,V,NS, V, NdU=TdSPdV+μdNdU = TdS - PdV + \mu dN
EnthalpyHHS,P,NS, P, NdH=TdS+VdP+μdNdH = TdS + VdP + \mu dN
Helmholtz free energyFFT,V,NT, V, NdF=SdTPdV+μdNdF = -SdT - PdV + \mu dN
Gibbs free energyGGT,P,NT, P, NdG=SdT+VdP+μdNdG = -SdT + VdP + \mu dN
Grand potentialΦ\PhiT,V,μT, V, \mudΦ=SdTPdVNdμd\Phi = -SdT - PdV - Nd\mu
  • Confusing CPC_P and CVC_V derivatives. CPCV=T(P/T)V(V/T)PC_P - C_V = T(\partial P/\partial T)_V (\partial V/\partial T)_P. For water near 44^\circC, (V/T)P=0(\partial V/\partial T)_P = 0 so CP=CVC_P = C_V, which is unusual. Always verify the relation for the specific substance.
  • Neglecting the chemical potential in open systems. In grand canonical ensembles, μ\mu is fixed by the reservoir. For photon gases, μ=0\mu = 0 because photon number is not conserved. For electrons, μ\mu is determined by the density and temperature via the Fermi-Dirac distribution.
  • Confusing thermodynamic potentials and their natural variables. U(S,V,N)U(S,V,N), F(T,V,N)F(T,V,N), G(T,P,N)G(T,P,N), and H(S,P,N)H(S,P,N) are the correct pairings. Using U(T,V,N)U(T,V,N) instead of F(T,V,N)F(T,V,N) leads to incorrect Maxwell relations and thermodynamic identities. The natural variables determine which potential is most convenient for a given process.

  • Forgetting the factor of kBk_B in entropy. The statistical definition S=kBlnΩS = k_B \ln \Omega is dimensionally correct only with Boltzmann’s constant. Without kBk_B, the numerical value of entropy has units of information (bits/nats), not energy per temperature.

  • Assuming all phase transitions are of first or second order. The 2D Ising model has a continuous phase transition with a logarithmic divergence of specific heat (α=0\alpha = 0), which does not fit neatly into the Ehrenfest classification. The modern classification uses order parameter behaviour and symmetry breaking rather than derivatives of free energy.

  • Confusing thermodynamic and kinetic stability. A system can be thermodynamically unstable (negative curvature of free energy) yet kinetically stable (metastable) due to energy barriers. Diamond at room temperature is metastable — it is kinetically stable but thermodynamically unstable relative to graphite. The phase diagram shows equilibrium, not kinetics.

  • Misapplying the Sackur-Tetrode equation. The Sackur-Tetrode equation for the entropy of an ideal gas assumes translational degrees of freedom only. For diatomic gases, rotational and vibrational contributions must be included at sufficiently high temperatures. The formula S=NkB[ln(V/Nλ3)+5/2]S = Nk_B[\ln(V/N\lambda^3) + 5/2] is valid only for monatomic ideal gases.

  • Neglecting the temperature dependence of the chemical potential. For fermions at low temperatures, μεF[1(π2/12)(kBT/εF)2]\mu \approx \varepsilon_F[1 - (\pi^2/12)(k_BT/\varepsilon_F)^2]. The Sommerfeld expansion shows that μ\mu decreases quadratically with TT, not linearly. Assuming μ=εF\mu = \varepsilon_F at all temperatures leads to errors in the electronic heat capacity calculation.