For bosons the occupation number is . Summing the geometric series for the single-mode grand canonical ensemble gives the Bose-Einstein distribution
The Bose-Einstein distribution is
Positivity requires not to exceed the ground-state energy, here zero. For a finite system above the condensation limit one has ; it approaches zero from below at Bose-Einstein condensation.
In a collapsing spacetime, positive-frequency in-modes are defined at past null infinity and outgoing modes at future null infinity. Tracing a late outgoing wave packet backward toward the forming event horizon produces an exponential blueshift. Its rapidly varying phase has both positive- and negative-frequency parts relative to the past time coordinate, so the Bogoliubov coefficient is nonzero. The resulting bosonic occupation numbers have the Bose-Einstein distribution at
where is the surface gravity. The outgoing quanta form Hawking radiation, while their correlated partners cross the horizon.
For fixed momentum put
In the contour formula, has poles at . Deforming the contour onto those poles and using the oddness of gives the standard Bosonic Matsubara sum
Equivalently,
where is the Bose-Einstein distribution. The first term is the zero-point fluctuation and the second is its thermal occupation.
For an ideal gas of identical bosons, the Bose-Einstein distribution gives the mean occupation of a discrete one-particle state as
where is the chemical potential. If the spectrum is sufficiently dense to use a density of states , the discrete sum becomes
For , the preceding density of states has
When , the Bose-Einstein distribution has the classical approximation
Each of the modes therefore contributes to the energy, so the Dulong-Petit law gives
When , the upper integration limit tends to infinity. The Bose integral
then yields
Hence the heat capacity at constant volume is
so . The high-temperature constant is the classical equipartition theorem result, while the low-temperature power law makes as required by the Third law of thermodynamics.
Subtract the zero-temperature term by using . In three dimensions the thermal part is
where the last integral uses the Bose-Einstein distribution. Gaussian phase fluctuations reduce the order parameter by . Taking symmetry restoration to occur when the thermal variance is of order one gives
The effective action in part c has , and hence . This is an order-of-magnitude estimate because near restoration the phase-only effective field theory omits large amplitude fluctuations and sensitivity to its ultraviolet cutoff.