Quantum Mechanics, Volume 3. Claude Cohen-Tannoudji
alt="image"/>
with, when the state of the system is given by the density operator (1):
This trace can be computed in the Fock state basis |n1, ..,ni..,nj,..〉 associated with the eigenstates basis {|uk〉} of
For i = j, this expression yields the average particle number in the individual state |ui〉.
2-a. Fermion distribution function
As the occupation number only takes the values 0 and 1, the first bracket in expression (17) is equal to [e–β(ei – μ)]; as for the other modes (k ≠ i) contribution, in the second bracket, it has already been computed when we determined the partition function. We therefore obtain:
(18)
Multiplying both the numerator and denominator by 1 + e–β (ei – μ) allows reconstructing the function Z in the numerator, and, after simplification by Z, we get:
We find again the Fermi-Dirac distribution function
(20)
This distribution function gives the average population of each individual state |ui〉 with energy e; its value is always less than 1, as expected for fermions.
The average value at thermal equilibrium of any one-particle operator is now readily computed by using (19) in relation (15).
2-b. Boson distribution function
The mode j = i contribution can be expressed as:
(21)
We then get:
(22)
which, using (11), amounts to:
where the Bose-Einstein distribution function
This distribution function gives the average population of the individual state |ui〉 with energy e. The only constraint of this population, for bosons, is to be positive. The chemical potential is always less than the lowest individual energy ek. In case this energy is zero, μ must always be negative. This avoids any divergence of the function
Hence for bosons, the average value of any one-particle operator is obtained by inserting (23) into relation (15).
2-c. Common expression
We define the function fβ as equal to either the function
where the number η is defined as:
(26)
2-d. Characteristics of Fermi-Dirac and Bose-Einstein distributions
We already gave in Complement CXIV (Figure 3) the form of the Fermi-Dirac distribution. Figure 1 shows both the variations of this distribution and the Bose-Einstein distribution. For the sake of comparison, it also includes the variations of the classical Boltzmann distribution:
which takes on intermediate values between the two quantum distributions. For a non-interacting gas contained in a box with periodic boundary conditions, the lowest possible energy e is zero and all the others are positive. Exponential eβ (ei – μ) is therefore always greater than e–βμ. We are now going to distinguish several cases, starting