# The Grand Canonical Ensemble

$$\newcommand{\vecs}{\overset { \rightharpoonup} {\mathbf{#1}} }$$ $$\newcommand{\vecd}{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}{\| #1 \|}$$ $$\newcommand{\inner}{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\kernel}{\mathrm{null}\,}$$ $$\newcommand{\range}{\mathrm{range}\,}$$ $$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$ $$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}{\| #1 \|}$$ $$\newcommand{\inner}{\langle #1, #2 \rangle}$$ $$\newcommand{\Span}{\mathrm{span}}$$$$\newcommand{\AA}{\unicode[.8,0]{x212B}}$$

To consider theories for fluctuations in the number of particles we require an ensemble that keeps V, T, and the chemical potential, m constant, a grand canonical ensemble. To construct the grand canonical ensemble, the system is enclosed in a container that is permeable both to heat and to the passage of particles. The number of particles in the system can range over all possible values. As in the canonical ensemble we have occupation numbers aNj describing the number of systems that have energy Ej and N particles. There are two indices for summation since neither the energy of the system nor the number of particles is the same in all of the systems. We can specify the state of the ensemble by specifying that aN1, aN2, aN3, of the systems are in states 1, 2, 3, …, respectively, with energies EN1, EN2, EN3,…. depending on N, the number of particles in each system. In the grand canonical ensemble the occupation numbers obey three conditions. Following the principle of a priori probabilities, we assume that every distribution of occupation number is equally probable. As we have seen for the canonical ensemble we can use the method of the most probable distribution to derive the form of the distribution function and the partition function for the grand ensemble.

The number of ways W(a) or W(aN1, aN2, aN3,….) that any particular distribution of aNj’s can be achieved is given by As before we assume that the systems are macroscopic and are distinguishable objects that can be distributed among available states. In any particular distribution aNj/A is the fraction of systems of the canonical ensemble in the jth energy state containing N particles. The overall probability PNj that a system is in the jth quantum state with N particles is obtained by averaging aNj/A over all the allowed distributions. The notation of summing over a means that the value of aj depends on the distribution and that the summations are over all distributions that satisfy the constraints. The most probable distribution is the distribution that maximizes W. The maximum W will be found by setting the derivative ( lnW/ aNi) = 0 subject to the constraints above that the aNj must sum to A, the total energy E is equal to the sum of aNjEj, and the total number of particles N is aNjN. This implies that In other words, there is no change in the total number of systems A, total energy E, and the total number of particles, N with respect to changes in the occupation numbers. The procedure followed here is analogous to that used for the canonical ensemble; we maximize W subject to the constraints. The difference is that there is one addition constraint on the number of particles, N that was not present in the canonical ensemble.

To maximize subject to constraints we use the method of LaGrange undetermined multipliers. where we have moved the summation symbol in front of the three terms. The constants a , –b , and g are the undetermined multipliers. We first carry out the derivative and then find the value of the multipliers.

We can evaluate

( lnW/ aNi) = ( lnA!/ aNi) – SN Sj ( ln aNj!/ aNi)

using Stirling's approximation: lnx! » xlnx – x.

To simplify ( lnW/ aNi) the first step is to note that

lnW = lnA! - SN Sj ln aNj! = AlnA – A - SN Sj aNj ln aNj - SN Sj aNj Since A = SN Sj aNJ the last two terms cancel to give

lnW = AlnA - SN Sj aNj ln aNj.

Note that exactly the same procedure and algebra are used to show that the derivative of ln W is equal to the most probable distribution in the canonical ensemble, and so we have: The most probable distribution is: Now we only need to find the undetermined multipliers a , b, and g . By summing both sides the indices N and j we can obtain a .

The left-hand side is equal to one and ea = 1/X where X is the grand canonical partition function or the grand partition function (for short). The Boltzmann distribution in this ensemble can be written The star indicates that this is the most probable distribution as shown above by maximizing W with respect to the occupation numbers.

The averages of mechanical properties E, P, and N are We have shown that b = 1/kBT where kB is Boltzmann’s constant. This is true for all of the ensembles. A similar approach can be used to show that g = -m /kBT. The differential of the grand partition function is Using the definitions of mechanical properties from above we have The last term is the ensemble-averaged work done by the system.

We add d( Eñ ) + d(g N) to both sides The thermodynamic equation

dE = TdS - PdV + m dN

can be rearranged to

TdS = dE - m dN + PdV Comparing these two equations by dividing through by b we find that

g = -mb = -m /kT from the second term in each equation.

Since we can replace g /b by -m we have which gives the entropy as Since

G = m N = E + PV - TS

We can use the above equation to determine that There are several ways to express the grand partition function. Starting with the definition, we can define the canonical partition function for N particles as and then insert this expression into the grand partition function  We know also that eg = em/kT . The quantity em/kT is often denoted l .

For indistinguishable particles Q = qN/N! where q is the molecular partition function. Therefore, 