 Solution

Each electromagnetic normal mode is a quantum harmonic oscillator with $E_n=(n+1/2)h\nu$. For $x=h\nu /(k_{B}T)$ its canonical partition function and occupation probabilities areThe geometric distribution is normalized and gives thermal mode energy $h\nu /(e^x-1)$. The zero-point term $h\nu /2$ is not part of the thermally emitted black-body intensity. The equilibrium photon chemical potential is zero because photon number is not conserved.

Counting two transverse polarizations gives $8\pi {\nu }^2/c^3$ modes per volume per frequency interval. Thus the spectral energy density is $u_{\nu }=8\pi h{\nu }^3/[c^3(e^x-1)]$. In an isotropic radiation field, spectral radiance is $B_{\nu }=c{u}_{\nu }/(4\pi )$. Transforming the spectral measure, $B_{\lambda }=B_{\nu }|d\nu /d\lambda |$, givesHere intensity means radiance per wavelength and solid angle; the hemispheric spectral exitance of a black surface is $M_{\lambda }=\pi B_{\lambda }$. The extra Jacobian is essential: the wavelength peak is not found by simply converting the frequency peak using $\lambda =c/\nu$.

At short wavelength the exponential suppresses radiation. At long wavelength the Rayleigh-Jeans limit is $B_{\lambda }\simeq 2c{k}_{B}T/{\lambda }^4$; using classical equipartition at every frequency instead would produce the ultraviolet catastrophe. The wavelength maximum solves $5(1-e^{-x})=x$, with nonzero root $x=4.9651$, yielding Wien's displacement law. For the specified temperature,Integrating the hemispheric flux gives the Stefan-Boltzmann law, $M=\sigma T^4\simeq 3.01\times 10^{-6}\mathrm{W}{\mathrm{m}}^{-2}$. The plot uses the specified ideal temperature, rather than substituting a later measured value. The physical black-body description of the cosmic microwave background is tested by comparing the sky spectrum with an absolute reference blackbody and separating instrumental and Galactic contributions.^[ref]