The canonical ensemble at work

Suppose we have a system whose Hamiltonian is:

(for example, can be the gravitational potential and an interaction between the particles). Then the canonical probability density is:
Note that and don't "get mixed up"; therefore, when calculating the partition function we can separate the kinetic and configurational contributions[1]:
Now, since:
we have:
In order to simplify this expression, we define the configurational partition function (or configurational sum) as:
and the thermal wavelength as:
so that, in general:

If we now consider an ideal gas, namely we set , then:

And using Stirling's approximation (see the appendix Stirling's approximation) the free energy of the system turns out to be:
Since from thermodynamics we know that , we get:
We have therefore found the state equation of an ideal gas in the canonical ensemble! We are thus really verifying that the canonical and microcanonical ensembles are equivalent. For example, we can also compute the entropy of this ideal gas, and we get:
which is exactly what we have found in the microcanonical ensemble.

Column of gas in a gravitational potential[edit | edit source]

Just as another example of the power of the canonical ensemble, we now want to see that in a particular case we can obtain the same results that we can determine from simple kinetic considerations. In particular, we want to see how the density of a gas subjected to a uniform gravitational potential varies with the height from the ground.

We begin with the "kinetic" computation. Consider a column of gas subjected to a constant gravitational potential , where is the height from the ground. The action of gravity will in some way change the density of the particles (we expect that the gas will be more dense near the ground and become thinner and thinner as grows) so let us call and the density and the pressure of the gas at height . The situation is represented in the following figure:

Column of gas subject to gravity

Let us call is the number of particles that are between the heights and , where is the particle number density at height ; these will be subjected to the pressures and , and if we call the cross section of the column then the net force acting on the particles is . Since the system is in equilibrium, this force must be balanced by the weight of the gas itself:

Therefore, dividing by and and taking the limit :
If we then consider a volume between and we have:
and since , then:
or equivalently:
Thus, the density of the gas decreases exponentially with the height.

We now want to show that the tools of the canonical ensemble allow us to obtain the very same result. Let us consider the single-particle canonical partition function[2]:

Since is such that and , integrating we will have:
The single-particle probability density of having momentum and height is:
Therefore we can obtain the spatial probability density integrating over , and ; in the end we get:
where we have relabelled the proportionality constant. Now, the particle number density at height will be , so in the end we have:
which is exactly the result previously found .

  1. This "separation" was not possible within the microcanonical ensemble and as we will later see (for example in Statistical mechanics of phase transitions) this will simplify the computations when we will consider non ideal systems, namely if the particles actually interact; in fact since any reasonable interaction does not involve the particles' momenta but only their positions only (see below) will be affected by the presence of the interaction.
  2. We use the single-particle partition function because in the end we want to find the single-particle configurational probability density.