Blackbody radiation is the equilibrium electromagnetic field that arises when photons and matter interact frequently enough to share a common temperature. In such environments, photons scatter repeatedly and exchange energy with particles until the entire system reaches thermal equilibrium. Once equilibrium is established, the radiation field depends only on temperature, not on the material composition of the emitting region. Stellar interiors, deep layers of stellar atmospheres, and the early Universe before photon decoupling all approximate these conditions.

When matter and radiation are in equilibrium, the average photon energy matches the characteristic kinetic energy of nonrelativistic particles. This condition can be expressed as

$$ \frac{3}{2} kT = \left( \frac{1}{2} mv^2 \right)_{av} = h \nu_{\mathrm{av}}. $$

A hypothetical container with walls and gas at temperature \( T \) would produce this equilibrium spectrum exactly. Although astrophysical environments are more complex than an ideal cavity, sufficiently high density and optical thickness allow photons to reach the same statistical equilibrium state.

The distribution of blackbody radiation follows from the Bose–Einstein statistics of photons. The photon distribution function describes how many photons occupy a phase-space cell of volume \( h^3 \) and is given by

$$ f = \frac{2}{h^3} \frac{1}{e^{h\nu/(kT)} - 1}. $$

From this distribution one obtains the specific intensity,

$$ I_\nu(T) = \frac{2 h \nu^3}{c^2} \frac{1}{e^{h\nu/(kT)} - 1}, $$

which is known as the Planck function. It gives the radiant energy per unit area, time, solid angle, and frequency interval. The spectrum rises as \( \nu^3 \) at low frequencies, reaches a well-defined maximum, and then declines rapidly at high frequencies. Expressed in wavelength units, the Planck function becomes

$$ I_\lambda(T) = \frac{2 h c^2}{\lambda^5} \frac{1}{e^{hc/(\lambda kT)} - 1}. $$

While these two forms describe the same radiation field, their peaks occur at different numerical locations because frequency and wavelength depend on each other nonlinearly.

At low frequencies where \( h\nu \ll kT \), the exponential in the denominator can be expanded using a Taylor series:

$$ e^{h\nu/(kT)} \approx 1 + \frac{h\nu}{kT}. $$

Substituting this into the Planck function yields the Rayleigh–Jeans approximation,

$$ I_\nu(T) \approx \frac{2 \nu^2 kT}{c^2}, $$

which grows quadratically with frequency and linearly with temperature. This form applies to radio and microwave wavelengths where photon energies are much smaller than the thermal energy scale.

At high frequencies where \( h\nu \gg kT \), the exponential term dominates and Planck’s law reduces to the Wien approximation,

$$ I_\nu(T) \approx \frac{2 h \nu^3}{c^2} e^{-h\nu/(kT)}. $$

This approximation describes the exponential falloff at short wavelengths. The location of the maximum of \( I_\nu(T) \) is determined by setting the derivative with respect to frequency equal to zero, which gives

$$ \nu_{\mathrm{peak}} = 5.88 \times 10^{10} T. $$

Thus the peak frequency scales directly with temperature.

Because \( I_\nu \) and \( I_\lambda \) have different functional forms, the peak of the wavelength spectrum does not occur at \( \lambda = c/\nu_{\mathrm{peak}} \). Differentiating the wavelength form of the Planck function yields the Wien displacement law,

$$ T \lambda_{\mathrm{peak}} = 2.898 \times 10^{-3} \ \mathrm{m\,K}. $$

Increasing the temperature therefore shifts the peak emission to shorter wavelengths. Cold sources of a few kelvin peak in the millimeter regime, while hotter bodies with temperatures of thousands of kelvin peak in the visible or ultraviolet.

Integrating the Planck function over all frequencies gives the total radiation energy density. The number density of photons varies as \( T^3 \), while the radiation energy density varies as \( T^4 \). The total flux emerging from a unit surface area is given by the Stefan–Boltzmann law,

\begin{equation} F = \sigma T^4. \end{equation}

A spherical star of radius \( R \) radiates a luminosity

$$ L = 4\pi R^2 \sigma T^4. $$

This connects the effective temperature of a radiating surface to its total output of energy and allows estimates of stellar radii and luminosities when the flux and temperature are known.

When blackbody spectra are plotted as functions of frequency, \(I(\nu)\), several consistent trends appear. At low frequencies the curves rise approximately as \(\nu^2\), following the Rayleigh–Jeans form that emerges from the classical limit of the Planck law. At high frequencies the intensity drops exponentially, producing the Wien tail whose slope steepens rapidly with increasing \(\nu\). As the temperature increases, the entire curve shifts upward and the peak moves systematically to higher frequencies, illustrating the strong temperature dependence of the radiation field. The total emitted power scales as \(T^4\), so hotter blackbodies produce far greater integrated intensity even if their spectral shapes remain similar.

In the wavelength representation, \(I(\lambda)\), the same underlying physics appears but with an important difference: the peak occurs at a different location from the frequency plot. Because frequency and wavelength are related by \(\nu=c/\lambda\), equal intervals in frequency do not correspond to equal intervals in wavelength. As a result, the wavelength-domain curve rises steeply at long wavelengths and falls sharply at short wavelengths, but the maximum occurs at a longer wavelength than one might infer from the frequency-domain peak. This difference directly reflects the Jacobian factor between the two representations, since \(I(\nu)\,d\nu \neq I(\lambda)\,d\lambda\). Wien’s displacement law, \(\lambda_{\max}T \approx 2.9\times10^{-3}\,\mathrm{m\,K}\), becomes particularly clear in the wavelength plot: increasing the temperature shifts the peak toward shorter wavelengths across the infrared, visible, and ultraviolet regions. Together, the two plots illustrate how the choice of independent variable changes the appearance of the same physical spectrum.

• Blackbody radiation arises when photons and matter reach full thermal equilibrium, producing a spectrum that depends only on temperature.
• The Planck function determines the intensity at each frequency and follows from the photon distribution in phase space.
• The Rayleigh–Jeans and Wien limits describe the low- and high-frequency behavior of the spectrum and explain its characteristic shape.
• The Wien displacement law shows that higher temperatures shift the peak emission to shorter wavelengths, linking temperature to observed colors.
• The Stefan–Boltzmann law connects temperature to total radiative flux, enabling luminosities and radii of stars to be inferred.
• Dense, optically thick astrophysical regions such as stellar atmospheres and the early Universe closely approximate blackbody conditions.
• Blackbody sources serve as calibration standards because their emission depends solely on temperature and follows universal laws.

• How does the equilibrium distribution of photons lead to the mathematical form of the Planck function?
• Why do the Rayleigh–Jeans and Wien approximations follow from the limiting behavior of the exponential term in Planck’s law?
• How does expressing the spectrum in frequency versus wavelength naturally lead to different peak locations?