The grand canonical ensemble allows particle number to fluctuate. The chemical potential μ controls the average particle number ⟨N⟩ via the grand partition function Ξ = Σ z^N Z_N where z = exp(μ/kT) is the fugacity. Fluctuations grow as √⟨N⟩.