distribution functions of Fermi-Dirac (for fermions, lower curve) and of Bose-Einstein (for bosons, upper curve) as a function of the dimensionless variable β(e – μ); the dashed line intermediate curve represents the classical Boltzmann distribution e–β (e – μ). In the right-hand side of the figure, corresponding to large negative values of μ, the particle number is small (the low density region) and the two distributions practically join the Boltzmann distribution. The system is said to be non-degenerate, or classical. As μ increases, we reach the central and left hand side of the figure, and the distributions become more and more different, reflecting the increasing gas degeneracy. For bosons, μ cannot be larger than the one-particle ground state energy, assumed to be zero in this case. The divergence observed for μ = 0 corresponds to Bose-Einstein condensation. For fermions, the chemical potential μ can increase without limit, and for all the energy values, the distribution function tends towards 1 (but never exceeding 1 due to the Pauli exclusion principle).
1 (iii) For a boson system, the chemical potential cannot be larger than the lowest e0 individual energy value, which we assumed to be zero. As μ tends towards zero through negative values and —kBT ≪ μ < 0, the distribution function denominator becomes very small leading to very large populations of the corresponding states. The boson gas is then said to be “degenerate”. On the other hand, for energies of the order or larger than μ, and as was the case for fermions, the boson distribution becomes practically equal to the Boltzmann distribution.
2 (iv) Finally, for situations intermediate between the extreme cases described above, the gas is said to be “partially degenerate”.
3. Two-particle operators
For a two-particle symmetric operator
with:
(29)
As the exponential operator in the trace is diagonal in the Fock basis states |n1, .., ni,.., nj,..〉, this trace will be non-zero on the double condition that the states i and j associated with the creation operator be exactly the same as the states k and l associated with the annihilation operators, whatever the order. In other words, to get a non-zero trace, we must have either i = l and j = k, or i = k and j = l, or both.
3-a. Fermions
As two fermions cannot occupy the same quantum state, the product
and the case i = k and j = l yields:
We begin with term (30). As i and j are different, operators
(32)
(33)
As for the second term (31), it is just the opposite of the first one. Consequently, we finally get:
The first term on the right-hand side is called the direct term. The second one is the exchange term, and has a minus sign, as expected for fermions.
3-b. Bosons
For bosons, the operators a commute with each other.
α. Average value calculation
If i ≠ j, a calculation, similar to the one we just did, yields:
which differs in two ways from (34): the result now involves the Bose-Einstein distribution, and the exchange term is positive.
If i = j, only one individual state comes into a new calculation, which we now perform. Using for ρeq expression (5) we get, after summing as in (11) a geometric series:
The sum appearing in this equation can be written as:
(37)
The first order derivative term yields:
(38)
and the second order derivative term is:
(39)
Summing these two
