Skip to content

The Grand Canonical Ensemble

3.1 Definition and Motivation

In many physical situations, a system exchanges both energy and particles with a reservoir. The grand canonical ensemble describes such open systems. The macroscopic variables are the chemical potential μ\muThe volume VVAnd the temperature TT.

Definition. The grand partition function is

Ξ=N=0ieβ(Ei(N)μN)\Xi = \sum_{N=0}^{\infty} \sum_{i} e^{-\beta(E_{i}^{(N)} - \mu N)}

Where the outer sum is over all possible particle numbers NN and the inner sum is over all states with NN particles.

The probability that the system is in state ii with NN particles is

Pi,N=eβ(Ei(N)μN)ΞP_{i,N} = \frac{e^{-\beta(E_{i}^{(N)} - \mu N)}}{\Xi}

3.2 Thermodynamic Relations

Theorem 3.1. The grand potential ΦG=kBTlnΞ\Phi_G = -k_BT \ln \Xi satisfies

ΦG=FμN=PV\Phi_G = F - \mu N = -PV

Proof. For a classical ideal gas, Ξ=N=0eβμNZN\Xi = \sum_{N=0}^{\infty} e^{\beta \mu N} Z_N where ZN=zN/N!Z_N = z^N/N! is the canonical partition function. Therefore:

Ξ=N=0(zeβμ)NN!=exp(zeβμ)\Xi = \sum_{N=0}^{\infty} \frac{(z e^{\beta \mu})^N}{N!} = \exp(z e^{\beta \mu})

ΦG=kBTlnΞ=kBTzeβμ=PV\Phi_G = -k_BT \ln \Xi = -k_BT \cdot z e^{\beta \mu} = -PV

The last equality follows from the ideal gas law PV=NkBTPV = Nk_BT with N=zeβμN = z e^{\beta \mu}. More generally, ΦG=PV\Phi_G = -PV holds for all systems. \blacksquare

Key relations from lnΞ\ln \Xi:

N=1βlnΞμT,V,E=lnΞβμ,V+μβlnΞμT,V\langle N \rangle = \frac{1}{\beta}\frac{\partial \ln \Xi}{\partial \mu}\bigg|_{T,V}, \quad \langle E \rangle = -\frac{\partial \ln \Xi}{\partial \beta}\bigg|_{\mu,V} + \frac{\mu}{\beta}\frac{\partial \ln \Xi}{\partial \mu}\bigg|_{T,V}

S=kB(lnΞ+βEβμN)S = k_B\left(\ln \Xi + \beta \langle E \rangle - \beta \mu \langle N \rangle\right)

3.3 Number Fluctuations

Theorem 3.2. The particle number fluctuations in the grand canonical ensemble satisfy

N2N2=kBT(Nμ)T,V\langle N^2 \rangle - \langle N \rangle^2 = k_BT \left(\frac{\partial \langle N \rangle}{\partial \mu}\right)_{T,V}

Proof. N2N2=1β22lnΞμ2=1βμ(1βlnΞμ)=1βNμ\langle N^2 \rangle - \langle N \rangle^2 = \frac{1}{\beta^2}\frac{\partial^2 \ln \Xi}{\partial \mu^2} = \frac{1}{\beta}\frac{\partial}{\partial \mu}\left(\frac{1}{\beta}\frac{\partial \ln \Xi}{\partial \mu}\right) = \frac{1}{\beta}\frac{\partial \langle N \rangle}{\partial \mu}. \blacksquare

For an ideal gas, N=zeβμ\langle N \rangle = z e^{\beta \mu}So N/μ=βN\partial \langle N \rangle / \partial \mu = \beta \langle N \rangleGiving relative fluctuations:

N2N2N2=1N\frac{\langle N^2 \rangle - \langle N \rangle^2}{\langle N \rangle^2} = \frac{1}{\langle N \rangle}

This is Poisson …/4-statistics-and-probability/2_statistics: fluctuations scale as 1/N1/\sqrt{N}Negligible for macroscopic systems.

3.4 Worked Example: Ideal Gas in the Grand Canonical Ensemble

Problem. Compute Ξ\Xi, N\langle N \rangleAnd E\langle E \rangle for a classical ideal gas in the grand canonical ensemble.

Solution

The single-particle partition function is z=V/λth3z = V/\lambda_{\mathrm{th}^3} where λth=h/2πmkBT\lambda_{\mathrm{th} = h/\sqrt{2\pi m k_BT}}. The canonical partition function for NN indistinguishable particles is ZN=zN/N!Z_N = z^N/N!. The grand partition function:

Ξ=N=0zNN!eβμN=N=0(zeβμ)NN!=ezeβμ\Xi = \sum_{N=0}^{\infty} \frac{z^N}{N!} e^{\beta \mu N} = \sum_{N=0}^{\infty} \frac{(ze^{\beta \mu})^N}{N!} = e^{ze^{\beta \mu}}

lnΞ=zeβμ=Vλth3eβμ\ln \Xi = ze^{\beta \mu} = \frac{V}{\lambda_{\mathrm{th}^3} e^{\beta \mu}}

Average particle number:

N=1βlnΞμ=Vλth3eβμ\langle N \rangle = \frac{1}{\beta}\frac{\partial \ln \Xi}{\partial \mu} = \frac{V}{\lambda_{\mathrm{th}^3} e^{\beta \mu}}

Solving for the chemical potential: μ=kBTln(Nλth3/V)\mu = k_BT \ln(\langle N \rangle \lambda_{\mathrm{th}^3 / V)}.

Average energy (using E=lnΞ/β+μN/(kBT)\langle E \rangle = -\partial \ln \Xi / \partial \beta + \mu \langle N \rangle / (k_BT)):

E=32NkBT\langle E \rangle = \frac{3}{2}\langle N \rangle k_BT

This recovers the equipartition result. \blacksquare