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.

The total power leaving \(1\,\mathrm{m}^2\) of a surface—known as the energy flux density \(\mathscr{F}(T)\) in units of W m\(^{-2}\)—is obtained by integrating the specific intensity \(I_\nu(T)\) over all frequencies and over all solid angles of the upper hemisphere. For a surface element \(dA\), the outward-directed flux is

$$ \mathscr{F} = \int_0^\infty \int_{\phi=0}^{2\pi} \int_{\theta=0}^{\pi/2} I_\nu(T)\, \cos\theta \, \sin\theta \, d\theta \, d\phi \, d\nu . $$

The integral over angle separates cleanly. Using

$$ \int_{\phi=0}^{2\pi} d\phi = 2\pi, \qquad \int_{\theta=0}^{\pi/2} \cos\theta \sin\theta \, d\theta = \frac{1}{2}, $$

the angular integration yields \( \pi \), so the flux becomes

$$ \mathscr{F} = \pi \int_0^\infty I_\nu(T)\, d\nu . $$

Substituting the Planck function,

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

gives

$$ \mathscr{F} = \int_0^\infty \frac{2\pi h \nu^3}{c^2} \frac{d\nu}{e^{h\nu/(kT)} - 1}. $$

Introducing the variable \( x = h\nu/(kT) \), the integral becomes proportional to

$$ \int_0^\infty \frac{x^3}{e^x - 1}\, dx , $$

whose value is \( \pi^4/15 \). Evaluating all constants produces the Stefan–Boltzmann law,

$$ \mathscr{F} = \sigma T^4, $$

where \( \sigma \) is the Stefan–Boltzmann constant,

$$ \sigma = \frac{2\pi^5 k^4}{15 c^2 h^3} = 5.670 \times 10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}} . $$

Thus the flux radiated per unit area of any ideal blackbody depends only on the fourth power of its surface temperature.

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

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

This expression gives the total radiant power emitted into space. In stellar applications, it connects the effective temperature of a star’s surface with its total luminosity. However, real stellar atmospheres may deviate from strict thermodynamic equilibrium, and their local temperature decreases outward. As a result, the observed flux often reflects a temperature slightly lower than the deeper, true thermal layers of the star. Absorption lines also form in cooler layers above the continuum-forming region, further influencing the emergent spectrum. Nevertheless, the Stefan–Boltzmann law remains the fundamental relation linking stellar radius, temperature, and luminosity.

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. At high frequencies the intensity falls exponentially, producing the Wien tail. As temperature increases, the entire curve shifts upward and the peak moves toward higher frequencies, reflecting both Wien’s law and the \(T^{4}\) scaling of the total emitted power.

In the wavelength representation, \(I(\lambda)\), the peak appears at a different location than in the frequency plot because frequency and wavelength are related by \(\nu = c/\lambda\). The difference arises from the Jacobian relating the two forms of the Planck function, since \(I(\nu)\,d\nu \neq I(\lambda)\,d\lambda\). Increasing the temperature moves the wavelength peak toward shorter wavelengths in accordance with Wien’s displacement law, \(\lambda_{\max}T \approx 2.9\times10^{-3}\,\mathrm{m\,K}\).

• 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?