Questions: The Grand Partition Function and Grand Thermodynamic Potential

Question 1

A student claims: 'The grand partition function is basically the canonical partition function — it just adds a fugacity correction factor to each term, but the underlying ensemble is the same.' What is fundamentally wrong with this claim?

• A. The claim is correct — for large systems the grand canonical and canonical ensembles give identical results
• B. The grand partition function sums over all possible particle numbers N, not a fixed N — it describes a different ensemble where particle number is a fluctuating quantity, not a constraint
• C. The canonical partition function already accounts for particle exchange through its dependence on chemical potential
• D. The fugacity correction only applies at high temperatures, so the claim is correct in the low-temperature limit

Answer: (B)

The grand partition function Ξ = Σ_{N,i} exp[−(E_i^(N) − μN)/kT] sums over every possible particle number N and every microstate within each N-particle sector. This is not a correction to the canonical ensemble — it is a completely different ensemble where N is a random variable with its own distribution. The canonical ensemble fixes N as a constraint; the grand canonical ensemble replaces that constraint with a chemical potential μ that controls the average N. The two ensembles agree on averages in the thermodynamic limit, but the grand canonical description is essential when particle exchange is physically real (e.g., open systems, quantum gases).

Question 2

You want to derive the mean occupation number of a single fermionic energy mode with energy ε. What structural feature of the grand partition function makes this straightforward?

• A. Fermions obey Maxwell-Boltzmann statistics at the single-mode level, so the sum reduces to a simple exponential
• B. For a single fermionic mode, the sum over particle number runs only from N = 0 to N = 1 due to the Pauli exclusion principle, making the sum trivial to evaluate
• C. The grand potential equals the Helmholtz free energy for non-interacting fermions, so no new calculation is needed
• D. Chemical potential is always negative for fermions, which terminates the sum after the N = 0 term

Answer: (B)

The Pauli exclusion principle limits occupancy of any single mode to 0 or 1 particle. So the grand partition function for that mode is simply Ξ_mode = e^{β(μ−ε)·0} + e^{β(μ−ε)·1} = 1 + e^{β(μ−ε)}. The mean occupation is then ⟨n⟩ = e^{β(μ−ε)} / (1 + e^{β(μ−ε)}) = 1/(e^{β(ε−μ)} + 1) — the Fermi-Dirac distribution, obtained with almost no calculation. The analogous sum for bosons (running from N = 0 to ∞) yields the Bose-Einstein distribution. Both iconic results emerge from the grand partition function's structure without additional machinery.

Question 3

The chemical potential μ plays the same role in the grand partition function that temperature plays in the canonical partition function: it controls the driving force for particle exchange between system and reservoir, just as temperature controls energy exchange.

• T. True
• F. False

Answer: True

This analogy is exact and important. Temperature T appears in the Boltzmann factor e^{−βE} and controls how readily a system accepts additional energy from a heat bath — high T means energy fluctuations are cheap. Chemical potential μ appears in the fugacity factor e^{βμN} and controls how readily a system accepts additional particles — high μ (relative to the single-particle energies) makes large N favorable. The grand partition function encodes both: Ξ = Σ_{N,i} exp[β(μN − E_i^(N))]. The parallel structure is the conceptual foundation for understanding why μ sets the equilibrium condition for particle exchange just as T sets the condition for thermal equilibrium.

Question 4

Near a thermodynamic critical point, the particle-number variance ⟨(ΔN)²⟩ = kT ∂⟨N⟩/∂μ will decrease significantly because the system becomes more ordered and resistive to fluctuations.

• T. True
• F. False

Answer: False

The opposite is true: near a critical point, particle-number fluctuations diverge rather than shrink. The variance ⟨(ΔN)²⟩ = kT ∂⟨N⟩/∂μ is large when ⟨N⟩ changes rapidly with μ — meaning the system is highly compressible and responsive to small changes in chemical potential. At a critical point, the distinction between phases (e.g., liquid and gas) vanishes and long-range correlations develop, causing the compressibility — and hence the fluctuations — to diverge. This divergence of fluctuations is a key signature of critical phenomena, not a sign of order.

Question 5

Why are the Fermi-Dirac and Bose-Einstein distributions described as 'emerging naturally' from the grand partition function, rather than requiring separate derivations? What feature of the grand canonical framework produces both distributions?

Model answer: Both distributions emerge from applying the grand partition function to a single quantum mode and evaluating the mean occupation ⟨n⟩ = −∂Ω/∂ε (or equivalently, computing the weighted average directly). For fermions, Pauli exclusion limits N_mode to 0 or 1, so the mode's grand partition function has just two terms; the mean occupation is 1/(e^{β(ε−μ)} + 1). For bosons, N_mode is unrestricted, so the sum is a geometric series that converges to 1/(e^{β(ε−μ)} − 1). The grand canonical framework produces both results because it naturally handles variable particle number — no additional assumptions about quantum statistics are needed beyond the allowed values of N.

This is the deepest payoff of the grand canonical formalism: quantum statistics are not extra assumptions bolted onto classical statistical mechanics but fall out automatically once you allow particle number to fluctuate and respect the constraints (Pauli exclusion for fermions, no constraint for bosons). The same mathematical machinery — sum over N with fugacity weighting — gives different results depending solely on the allowed occupation numbers. This is why the grand canonical ensemble is the natural framework for all quantum gases, and why it is introduced before studying Fermi gases and Bose-Einstein condensation.