Blackbody radiation: Planck distribution, Stefan-Boltzmann law, Wien displacement law

Intuition Beginner

A blackbody is a perfect absorber of electromagnetic radiation: every photon that lands on it is absorbed, none reflected. Heat the blackbody to temperature $T$ and it glows — not because of its chemistry, but purely because of its temperature. The light it emits is called blackbody radiation, and its colour depends only on $T$.

Look at a hot stove element. Cold, it is dark. Warm it to a few hundred degrees and it glows a dull red. Push the temperature higher and it brightens to orange, then yellow-white. The sun, at about 5800 K, glows white with a peak in the yellow-green band of visible light. Hotter stars (Sirius, 9900 K) shimmer blue. Cooler bodies (your skin at 310 K, the cosmic microwave background at 2.725 K) radiate too — but in the infrared and microwave, invisible to the eye.

Two simple laws describe the glow. Stefan-Boltzmann says the total power per unit area radiated by a blackbody grows as $T$ to the fourth power. Double the temperature, the brightness multiplies by sixteen. Wien displacement says the wavelength of peak emission shrinks in inverse proportion to $T$: hotter bodies peak at shorter wavelength, bluer light. Both laws were discovered empirically in the late 1800s.

Classical physics could not explain them. Classical electromagnetism predicts that a hot cavity should radiate ever more energy as the frequency rises, with no cutoff — the so-called ultraviolet catastrophe. Experiments showed the opposite: the spectrum has a peak and falls off at high frequency. In 1900 Max Planck proposed a fix: electromagnetic energy at frequency $f$ comes only in indivisible chunks of size $hf$, where $h$ is a new constant of nature. High-frequency modes, requiring large energy chunks, are exponentially suppressed at moderate temperature. The catastrophe vanishes, and the observed spectrum is reproduced exactly. Quantum mechanics was born.

Visual Beginner

Three curves on the same axes, each plotting blackbody intensity as a function of wavelength at a different temperature. A cool body (about 1000 K) peaks far to the right, deep in the infrared, with a tiny red tail visible to the eye. The sun-like curve (about 5800 K) peaks in the middle of the visible band. A hot star (about 10000 K) peaks in the ultraviolet with most of its visible output appearing white-blue. As the temperature rises, the peak shifts left (Wien) and the total area under the curve grows dramatically (Stefan-Boltzmann).

A second panel shows the same three temperatures plotted against frequency. The classical Rayleigh-Jeans prediction is also drawn for the 5800 K case: a parabola that rises without bound. The Planck curve agrees with Rayleigh-Jeans at low frequency, then turns over and decays exponentially, averting the ultraviolet catastrophe.

Worked example Beginner

The cosmic microwave background (CMB) is the relic photon gas filling the universe at temperature $T=2.725\text{K}$. Compute its peak wavelength and the total power it radiates per unit area.

Step 1. Wien displacement law gives the peak wavelength: ${\lambda }_{\max }=b/T$, where $b=2.898\times 10^{-3}\text{m}\cdot \text{K}$ is the Wien displacement constant.

Step 2. Plug in $T=2.725\text{K}$:

$${\lambda }_{\max }=\frac{2.898\times 10^{-3}}{2.725}\approx 1.063\times 10^{-3}\text{m}\approx 1.06\text{mm}.$$

Step 3. Stefan-Boltzmann gives the radiated power per unit area: $j=\sigma T^4$, with $\sigma =5.670\times 10^{-8}\text{W}/({\text{m}}^2\cdot {\text{K}}^4)$. Then

$$j=5.670\times 10^{-8}\times (2.725{)}^4\approx 5.670\times 10^{-8}\times 55.15\approx 3.13\times 10^{-6}\text{W}/{\text{m}}^2.$$

What this tells us: the CMB peaks at about 1 mm wavelength — in the microwave band, which is why it was first detected as unexplained antenna noise by Penzias and Wilson in 1965. The total radiated power is tiny, only a few microwatts per square metre, because the gas is so cold. For comparison, the sun at 5800 K radiates $\sigma T^4\approx 6.3\times 10^7\text{W}/{\text{m}}^2$ — twenty trillion times brighter per unit area.

Check your understanding Beginner

Formal definition Intermediate+

Consider an isotropic cavity of volume $V$ at thermal equilibrium with temperature $T$, containing electromagnetic radiation but no other matter. The walls absorb and re-emit photons freely, so photon number is not conserved; the equilibrium chemical potential is $\mu =0$. Each cavity mode, labelled by wave vector $\mathbf{k}$ and polarisation, is a quantum harmonic oscillator of frequency $\omega =c|\mathbf{k}|$ with energy levels $n\hslash \omega$ for $n=0,1,2,\dots$.

Definition (Planck spectral energy density). The energy per unit volume per unit angular frequency interval in a blackbody cavity at temperature $T$ is

$$u(\omega ,T)=\frac{\hslash }{{\pi }^2{c}^3}\frac{{\omega }^3}{e^{\hslash \omega /k_{B}T}-1}.$$

In terms of ordinary frequency $\nu =\omega /(2\pi )$ this reads

$$u(\nu ,T)d\nu =\frac{8\pi h{\nu }^3}{c^3}\frac{d\nu }{e^{h\nu /k_{B}T}-1},$$

and in terms of wavelength $\lambda =c/\nu$,

$$u(\lambda ,T)d\lambda =\frac{8\pi hc}{{\lambda }^5}\frac{d\lambda }{e^{hc/(\lambda k_{B}T)}-1}.$$

The Planck distribution interpolates between two classical limits. At low frequency $h\nu \ll k_{B}T$, expanding the denominator $e^x-1\approx x$ gives the Rayleigh-Jeans law $u(\nu ,T)\approx 8\pi {\nu }^2{k}_{B}T/c^3$ — quadratic in $\nu$, divergent when integrated to infinity. At high frequency $h\nu \gg k_{B}T$, the exponential dominates and $u(\nu ,T)\approx (8\pi h{\nu }^3/c^3)e^{-h\nu /k_{B}T}$ — the Wien spectrum that he proposed empirically in 1896. Planck's formula was the first expression that fit experimental data across the whole frequency range.

Stefan-Boltzmann law

Integrating $u(\nu ,T)$ over all frequencies gives the total energy per unit volume

$$u(T)={\int }_0^{\infty }u(\nu ,T)d\nu =\frac{8{\pi }^5{k}_B^4}{15h^3{c}^3}{T}^4\equiv a{T}^4,$$

where $a=7.566\times 10^{-16}\text{J}/({\text{m}}^3\cdot {\text{K}}^4)$ is the radiation constant. A blackbody surface radiates the cavity energy density outward at speed $c/4$ (one quarter for an isotropic gas hitting a planar surface, including the $\cos \theta$ projection averaged over the hemisphere), giving the Stefan-Boltzmann law:

$$j(T)=\frac{c}{4}u(T)=\sigma T^4,\sigma =\frac{2{\pi }^5{k}_B^4}{15h^3{c}^2}=5.670374419\times 10^{-8}\frac{\text{W}}{{\text{m}}^2\cdot {\text{K}}^4}.$$

Wien displacement law

The peak of $u(\lambda ,T)$ in wavelength is found by setting $du/d\lambda =0$. Let $x=hc/(\lambda k_{B}T)$; differentiating $u\propto {\lambda }^{-5}/(e^x-1)$ and simplifying yields the transcendental equation

$$x=5(1-e^{-x}),$$

with numerical root $x\approx 4.96511$. Therefore

$${\lambda }_{\max }T=\frac{hc}{x{k}_B}=2.897771955\times 10^{-3}\text{m}\cdot \text{K}\equiv b.$$

This is the Wien displacement law: the wavelength of peak emission is inversely proportional to temperature. The corresponding peak in frequency satisfies $h{\nu }_{\max }/(k_{B}T)=2.82144$ (a different transcendental root because the change of variables $\nu =c/\lambda$ is nonlinear).

Counterexamples to common slips

• Confusing ${\lambda }_{\max }$ in frequency and wavelength. The frequency peak ${\nu }_{\max }=2.821k_{B}T/h$ does not correspond to the wavelength $c/{\nu }_{\max }$; that would give $\lambda \cdot T=5.10\times 10^{-3}\text{m}\cdot \text{K}$, not $2.898\times 10^{-3}$. The peak position depends on which spectral variable one differentiates, because $u(\nu )d\nu$ and $u(\lambda )d\lambda$ are different distributions related by the Jacobian $d\nu /d\lambda$.
• Assuming photon number is conserved. It is not: photons can be created or destroyed at the cavity walls, so the chemical potential must be set to zero. For a massive Bose gas with fixed particle number, $\mu <0$ and the spectrum differs.
• Treating the spectrum as classical. The factor $1/(e^{h\nu /k_{B}T}-1)$ is intrinsically quantum: it is the average occupation of a harmonic oscillator with energy levels $nh\nu$ at temperature $T$. Replacing it with the classical equipartition value $k_{B}T/(h\nu )$ recovers Rayleigh-Jeans and the ultraviolet catastrophe.

Key theorem with proof Intermediate+

Theorem (Planck spectrum from the canonical ensemble). The spectral energy density of electromagnetic radiation in thermal equilibrium at temperature $T$ in a cavity of volume $V$ is

$$u(\omega ,T)=\frac{\hslash }{{\pi }^2{c}^3}\frac{{\omega }^3}{e^{\hslash \omega /k_{B}T}-1},$$

and the total energy density is $u(T)=a{T}^4$ with $a={\pi }^2{k}_B^4/(15{\hslash }^3{c}^3)$. As a corollary, the radiated flux from a blackbody surface satisfies the Stefan-Boltzmann law $j=\sigma T^4$ with $\sigma ={\pi }^2{k}_B^4/(60{\hslash }^3{c}^2)\equiv ac/4$, and the wavelength peak satisfies Wien displacement ${\lambda }_{\max }T=b$ with $b=hc/(4.96511k_B)$.

Proof. We work in three steps: mode counting, single-mode statistics, and the spectral integral.

Step 1 (mode counting). In a cubic cavity of side $L$ with periodic boundary conditions, the allowed wave vectors are $\mathbf{k}=(2\pi /L)(n_1,n_2,n_3)$ for $n_i\in \mathbb{Z}$. Each wave vector carries two transverse polarisation states. The number of modes with $|\mathbf{k}|$ between $k$ and $k+dk$ is the volume of a spherical shell in $\mathbf{k}$-space, divided by the volume per mode $(2\pi /L{)}^3$, times two polarisations:

$$dN=2\cdot \frac{4\pi k^{2}dk}{(2\pi /L{)}^3}=\frac{V{k}^2}{{\pi }^2}dk.$$

Substituting $\omega =ck$ and $d\omega =cdk$ yields the mode density in frequency:

$$g(\omega )d\omega =\frac{V{\omega }^2}{{\pi }^2{c}^3}d\omega .$$

Step 2 (single-mode statistics). Each mode of frequency $\omega$ is a quantum harmonic oscillator with energy levels $E_n=n\hslash \omega$ (we drop the zero-point energy $\frac{1}{2}\hslash \omega$, which is a constant offset and cannot be radiated). The canonical partition function is a geometric series:

$$Z_{\omega }=\sum _{n=0}^{\infty }{e}^{-\beta n\hslash \omega }=\frac{1}{1-e^{-\beta \hslash \omega }}.$$

The mean energy per mode follows from $⟨E⟩=-\partial \ln Z/\partial \beta$:

$$⟨E(\omega )⟩=\frac{\hslash \omega }{e^{\beta \hslash \omega }-1}.$$

This is the Bose-Einstein occupation $⟨n(\omega )⟩=1/(e^{\beta \hslash \omega }-1)$ at $\mu =0$ multiplied by the photon energy $\hslash \omega$; the photon-gas chemical potential vanishes because photon number is not a conserved quantity (the cavity walls freely create and absorb photons to maintain thermal equilibrium).

Step 3 (spectral integral). Combining the mode density with the mean energy per mode gives the spectral energy density (energy per unit volume per unit angular frequency):

$$u(\omega ,T)=\frac{1}{V}g(\omega )⟨E(\omega )⟩=\frac{{\omega }^2}{{\pi }^2{c}^3}\cdot \frac{\hslash \omega }{e^{\beta \hslash \omega }-1}=\frac{\hslash }{{\pi }^2{c}^3}\frac{{\omega }^3}{e^{\beta \hslash \omega }-1}.$$

For the total energy density, substitute $x=\beta \hslash \omega$:

$$u(T)={\int }_0^{\infty }u(\omega ,T)d\omega =\frac{(k_{B}T{)}^4}{{\pi }^2{\hslash }^3{c}^3}{\int }_0^{\infty }\frac{x^3}{e^x-1}dx.$$

The integral evaluates by expanding $1/(e^x-1)={\sum }_{n=1}^{\infty }{e}^{-nx}$ and integrating term-by-term:

$${\int }_0^{\infty }\frac{x^3}{e^x-1}dx=\sum _{n=1}^{\infty }{\int }_0^{\infty }{x}^3{e}^{-nx}dx=\sum _{n=1}^{\infty }\frac{6}{n^4}=6\zeta (4)=\frac{{\pi }^4}{15}.$$

Therefore $u(T)={\pi }^2(k_{B}T{)}^4/(15{\hslash }^3{c}^3)=a{T}^4$ with $a={\pi }^2{k}_B^4/(15{\hslash }^3{c}^3)$.

For the Stefan-Boltzmann surface flux, a planar opening into the cavity emits photons with mean perpendicular speed $c⟨\cos \theta {⟩}_{\mathrm{hemisphere}}=c/2$, multiplied by the factor $1/2$ for half-space directions: the net outward flux is $j=(c/4)u(T)=\sigma T^4$ with $\sigma =ac/4={\pi }^2{k}_B^4/(60{\hslash }^3{c}^2)=5.670\times 10^{-8}\text{W}/({\text{m}}^2\cdot {\text{K}}^4)$.

For Wien displacement, write $u$ in wavelength form via $u(\nu )d\nu =u(\lambda )d\lambda$ with $\nu =c/\lambda$. The result is $u(\lambda )\propto {\lambda }^{-5}/(e^{hc/\lambda k_{B}T}-1)$. Setting $du/d\lambda =0$ and substituting $x=hc/(\lambda k_{B}T)$ gives

$$\frac{d}{d\lambda }\left[{\lambda }^{-5}{\left(e^{hc/\lambda k_{B}T}-1\right)}^{-1}\right]=0⟺5-x\frac{e^x}{e^x-1}=0⟺x=5(1-e^{-x}).$$

Numerical iteration of $x_{n+1}=5(1-e^{-x_n})$ starting from $x_0=5$ converges to $x\approx 4.96511$, so ${\lambda }_{\max }T=hc/(x{k}_B)=2.898\times 10^{-3}\text{m}\cdot \text{K}$. $\square$

Bridge. The Planck spectrum builds toward 11.05.01 Bose-Einstein distribution, of which it is the $\mu =0$ photon-gas special case. The mode-density factor $V{\omega }^2/({\pi }^2{c}^3)$ appears again in 12.13.01 Bosonic Fock space as the photon-state density of the electromagnetic field operator algebra. The foundational reason that the spectrum has a finite total is the quantum hypothesis $E=n\hslash \omega$: this is exactly the high-frequency cutoff that averts the ultraviolet catastrophe of Rayleigh-Jeans. The proof generalises directly to phonons (replacing $c$ by the sound speed and adding the Debye cutoff) and to any gas of massless bosons in $d$ spatial dimensions, where the energy density scales as $T^{d+1}$ by dimensional analysis applied to the same integral structure.

Exercises Intermediate+

Advanced results Master

Theorem 1 (Kirchhoff's law of thermal radiation). At thermal equilibrium, the ratio of emissivity $\epsilon (\nu )$ to absorptivity $\alpha (\nu )$ at frequency $\nu$ is the same for all bodies and equal to the Planck function $B(\nu ,T)$. Equivalently, $\epsilon (\nu )=\alpha (\nu )$, and a perfect absorber ($\alpha =1$) is a perfect emitter at every frequency — the blackbody is the maximum-emission bound at fixed temperature.

Kirchhoff (1859) proved this from a thought experiment with two bodies inside a closed cavity in thermal equilibrium: if either body's $\epsilon /\alpha$ ratio exceeded the other's, heat would flow from cold to hot through radiation, violating the second law. The universal function $B(\nu ,T)$ was the central unknown of 19th-century theoretical physics; Planck's 1900 derivation fixed its form.

Theorem 2 (photon-gas equation of state). The photon gas obeys $PV=U/3$, equivalently $P=a{T}^4/3$. The Helmholtz free energy is $F=-\frac{1}{3}a{T}^{4}V$, and the chemical potential vanishes: $\mu =(\partial F/\partial N{)}_{T,V}=0$ identically because photon number is not an independent thermodynamic variable.

The factor $1/3$ contrasts with the non-relativistic ideal gas's $P=(2/3)U/V$ — both follow from the general result $P=U/(dV)$ for ultra-relativistic gases and $P=(2/d)U/V$ for non-relativistic ones in $d$ spatial dimensions. The photon-gas equation of state is responsible for the radiation-pressure-dominated era of early-universe cosmology: at $z\gtrsim 3400$, ${\rho }_{\gamma }>{\rho }_m$ and the universe's expansion was driven by radiation rather than matter.

Theorem 3 (Wien displacement and adiabatic invariance). For reversible adiabatic compression or expansion of a photon gas, $T{V}^{1/3}=\text{const}$, equivalently $T\lambda =\text{const}$ for each spectral mode. The spectrum retains its Planck form throughout the process; only the temperature parameter shifts. This is Wien's 1893 derivation of the displacement law from adiabatic invariance, predating Planck.

The proof uses adiabatic invariants of the harmonic-oscillator modes: under slow (adiabatic) change of the cavity size, the action $\oint pdq=E/\omega =n\hslash$ is preserved, so each mode's quantum number $n$ stays fixed while $\omega$ scales with the cavity size. Wien obtained displacement before quantum theory existed; the modern reading is that quantum numbers are exactly the adiabatic invariants.

Theorem 4 (Einstein 1916 detailed balance). In thermal equilibrium with blackbody radiation, the ratio of Einstein coefficients for a two-level system with energy gap $\hslash {\omega }_0$ is

$$\frac{A}{B}=\frac{\hslash {\omega }_0^3}{{\pi }^2{c}^3}.$$

This is the same factor that appears in the Planck distribution; the underlying reason is mode-density. In frequency form $A/B=8\pi h{\nu }^3/c^3$.

The result is foundational for laser physics: it shows that spontaneous and stimulated emission have related strengths fixed by the photon-mode density, and that population inversion ($N_2>N_1$) is required for net amplification. Einstein 1916 anticipated the laser by 44 years through this argument alone, without knowing the underlying field-quantisation framework.

Theorem 5 (CMB blackbody and COBE-FIRAS). The cosmic microwave background spectrum is a Planck blackbody at

$$T_{\mathrm{CMB}}=2.72548\pm 0.00057\text{K},$$

with deviations from the ideal Planck form smaller than $10^{-5}$ across the FIRAS frequency range (60–600 GHz). This makes the CMB the most precisely-measured blackbody spectrum in any laboratory or natural setting.

Mather et al. (1994) reported the original FIRAS result with deviations bounded by $\Delta u/u<5\times 10^{-5}$; refined COBE analysis (Fixsen 2009) tightened the bound to $|\Delta u/u|<10^{-5}$ except for a known dipole anisotropy from solar motion through the CMB rest frame. The thermalisation of the CMB during recombination ($z\approx 1100$) is a near-perfect example of Kirchhoff's universal-function principle realised in nature.

Theorem 6 (Stefan-Boltzmann from electromagnetic thermodynamics). Boltzmann (1884) derived the $T^4$ law before Planck's spectrum was known, using thermodynamic arguments on the energy density and pressure of radiation. Starting from $u=u(T)$ and $P=u/3$ (deducible from Maxwell's stress tensor without quantisation), the first and second laws give

$${\left(\frac{\partial u}{\partial V}\right)}_T=T{\left(\frac{\partial P}{\partial T}\right)}_V-P,$$

which reduces to $u=T{u}^{\prime }(T)/3-u/3$, i.e., $T{u}^{\prime }(T)=4u$, with solution $u(T)=a{T}^4$.

Boltzmann's 1884 paper is a landmark of classical thermodynamics: a major physical law derived without statistical mechanics, using only Maxwell's stress tensor for radiation pressure and the second law for the cycle. Planck's 1900 derivation later fixed the constant $a$ in terms of the elementary constants $h,c,k_B$.

Synthesis. The blackbody calculations are the foundational reason that statistical mechanics matters for physics beyond classical gases: this is exactly the case where averaging over a continuum of harmonic-oscillator modes with quantised energy levels produces an experimentally observed law that classical physics cannot reach. The central insight is that the same Bose-Einstein occupation number $⟨n(\omega )⟩=1/(e^{\hslash \omega /k_{B}T}-1)$, applied to photons (chemical potential zero, since photon number is not conserved), reproduces Wien displacement, Stefan-Boltzmann, and the full Planck spectrum together. Putting these together with the radiation equation of state $P=u/3$, the cosmological redshift relation $T\propto a^{-1}$ follows from adiabatic invariance, and the CMB temperature today is fixed by the entropy of the universe per comoving volume. The bridge is between thermodynamics (Boltzmann 1884), wave optics (Wien 1893, 1896), and quantum statistics (Planck 1900, Bose-Einstein 1924) — three independent derivations that converge on the same $a{T}^4$ result. Einstein's 1916 detailed-balance argument identifies the spontaneous-emission rate $A$ with the stimulated rate $B$ times $8\pi h{\nu }^3/c^3$, the same mode-density factor appearing in the Planck function; this pattern recurs throughout quantum optics as the mode density of the electromagnetic vacuum, and generalises to phonons, magnons, and any gas of massless bosons via the same dispersion-relation arithmetic.

Full proof set Master

Proposition 1 (Stefan-Boltzmann from thermodynamics). For any system with $u=u(T)$ depending only on temperature and $P=u/3$, the energy density obeys $u(T)=a{T}^4$ for some constant $a$.

Proof. The thermodynamic identity (a Maxwell relation derived from $dU=TdS-PdV$) reads

$${\left(\frac{\partial U}{\partial V}\right)}_T=T{\left(\frac{\partial P}{\partial T}\right)}_V-P.$$

For radiation, $U/V=u(T)$ depends only on $T$, so $(\partial U/\partial V{)}_T=u(T)$. With $P=u/3$, the right-hand side is $(T/3)u^{\prime }(T)-u/3$. Equating:

$$u=\frac{T}{3}{u}^{\prime }(T)-\frac{u}{3}⟺T{u}^{\prime }(T)=4u(T).$$

The differential equation $du/dT=4u/T$ separates: $du/u=4dT/T$, hence $u(T)=a{T}^4$ for a constant of integration $a>0$. $\square$

Proposition 2 (Wien displacement transcendental root). The equation $x=5(1-e^{-x})$ has a unique positive root $x^ \approx 4.96511$,andthewavelengthdisplacementconstantis$b = hc/(x^* k_B) = 2.897771955 \times 10^{-3},\text{m}\cdot\text{K}$.*

Proof. Let $\varphi (x)=x-5(1-e^{-x})=x-5+5e^{-x}$. Compute $\varphi (0)=0$, ${\varphi }^{\prime }(x)=1-5e^{-x}$, ${\varphi }^{\prime }(0)=-4<0$. As $x\to \infty$, $\varphi (x)\to +\infty$. The derivative vanishes at $x_c=\ln 5\approx 1.609$, where $\varphi (x_c)=\ln 5-5+1=\ln 5-4\approx -2.39$. So $\varphi$ decreases from 0 at $x=0$ to $\varphi (x_c)\approx -2.39$, then increases monotonically toward $+\infty$. Therefore $\varphi$ has exactly one positive zero, in $(x_c,\infty )$.

Numerical bisection or fixed-point iteration $x_{n+1}=5(1-e^{-x_n})$ from $x_0=5$ converges rapidly: $x_1=5(1-e^{-5})=5(1-0.00674)\approx 4.9663$; $x_2=5(1-e^{-4.9663})\approx 4.96512$; $x_3\approx 4.96511$. With Planck's constant $h=6.62607015\times 10^{-34}\text{J}\cdot \text{s}$, speed of light $c=2.998\times 10^8\text{m}/\text{s}$, and Boltzmann constant $k_B=1.381\times 10^{-23}\text{J}/\text{K}$,

$$b=\frac{hc}{x^{\ast }{k}_B}=\frac{6.626\times 10^{-34}\times 2.998\times 10^8}{4.96511\times 1.381\times 10^{-23}}=2.8978\times 10^{-3}\text{m}\cdot \text{K}.\square$$

Proposition 3 (photon-number density formula). The number of photons per unit volume in a blackbody at temperature $T$ is

$$n_{\gamma }(T)=\frac{2\zeta (3)}{{\pi }^2}{\left(\frac{k_{B}T}{\hslash c}\right)}^3,$$

with $\zeta (3)\approx 1.20206$, the Apéry constant.

Proof. Sum the Bose-Einstein occupation against the mode density:

$$n_{\gamma }=\frac{1}{V}{\int }_0^{\infty }g(\omega )⟨n(\omega )⟩d\omega =\frac{1}{{\pi }^2{c}^3}{\int }_0^{\infty }\frac{{\omega }^2}{e^{\beta \hslash \omega }-1}d\omega .$$

Substitute $x=\beta \hslash \omega$:

$$n_{\gamma }=\frac{(k_{B}T{)}^3}{{\pi }^2{\hslash }^3{c}^3}{\int }_0^{\infty }\frac{x^2}{e^x-1}dx.$$

Expand $1/(e^x-1)={\sum }_{n=1}^{\infty }{e}^{-nx}$ and integrate:

$${\int }_0^{\infty }\frac{x^2}{e^x-1}dx=\sum _{n=1}^{\infty }\frac{2}{n^3}=2\zeta (3).$$

So $n_{\gamma }=2\zeta (3)(k_{B}T/\hslash c{)}^3/{\pi }^2$. The Apéry constant $\zeta (3)$ is irrational (Apéry 1979); its appearance in a physical observable is unusual and the photon-number formula is one of the few elementary expressions where it occurs. $\square$

Proposition 4 (photon-gas entropy density). The entropy density of the photon gas is $s(T)=(4a/3)T^3$ with $a={\pi }^2{k}_B^4/(15{\hslash }^3{c}^3)$. Equivalently, the entropy per photon is $⟨s⟩=(2{\pi }^4/45)(k_B/2\zeta (3))\approx 3.602k_B$.

Proof. Photon number is not a conserved quantity, so the chemical potential vanishes and the Gibbs-Duhem relation reduces to $U+PV=TS$ (no $\mu N$ term). Using $P=u/3$ from Proposition 2's framework,

$$s=\frac{S}{V}=\frac{u+P}{T}=\frac{u+u/3}{T}=\frac{4u}{3T}=\frac{4a{T}^3}{3}.$$

For entropy per photon, divide by the photon number density: $⟨s⟩=s/n_{\gamma }=(4a{T}^3/3)/[(2\zeta (3)/{\pi }^2)(k_{B}T/\hslash c{)}^3]$. The $T^3$ factors cancel:

$$⟨s⟩=\frac{4}{3}\cdot \frac{{\pi }^2{k}_B^4}{15{\hslash }^3{c}^3}\cdot \frac{{\pi }^2{\hslash }^3{c}^3}{2\zeta (3)k_B^3}=\frac{4{\pi }^4}{45\cdot 2\zeta (3)}{k}_B=\frac{2{\pi }^4}{45\zeta (3)}{k}_B\approx 3.602k_B.$$

Each photon in a blackbody carries about $3.6k_B\approx 5\times 10^{-23}\text{J/K}$ of entropy, independent of the radiation temperature. The vanishing of $s$ as $T\to 0$ confirms the third law: in the vacuum state with no photons, entropy is zero. $\square$

Proposition 5 (Stefan-Boltzmann constant value). The Stefan-Boltzmann constant evaluates to

$$\sigma =\frac{2{\pi }^5{k}_B^4}{15h^3{c}^2}=5.670374419\times 10^{-8}\frac{\text{W}}{{\text{m}}^2\cdot {\text{K}}^4}$$

using the 2019 SI-fixed values $h=6.62607015\times 10^{-34}\text{J}\cdot \text{s}$, $k_B=1.380649\times 10^{-23}\text{J/K}$, $c=2.99792458\times 10^8\text{m/s}$.

Proof. The radiation constant is $a={\pi }^2{k}_B^4/(15{\hslash }^3{c}^3)=8{\pi }^5{k}_B^4/(15h^3{c}^3)$ after substituting $\hslash =h/2\pi$ (the factor $(2\pi {)}^3=8{\pi }^3$ multiplies, combined with the ${\pi }^2$ to give $8{\pi }^5$). The Stefan-Boltzmann constant is $\sigma =ac/4=2{\pi }^5{k}_B^4/(15h^3{c}^2)$. Numerically:

$$\sigma =\frac{2{\pi }^5(1.380649\times 10^{-23}{)}^4}{15(6.62607015\times 10^{-34}{)}^3(2.99792458\times 10^8{)}^2}.$$

Carrying out the arithmetic: numerator $=2{\pi }^5\times 3.633\times 10^{-92}\approx 2.222\times 10^{-90}$; denominator $=15\times 2.910\times 10^{-100}\times 8.988\times 10^{16}\approx 3.922\times 10^{-83}$. The ratio is $5.670\times 10^{-8}\text{W}/({\text{m}}^2\cdot {\text{K}}^4)$, as quoted in CODATA tables. The value is now exact in the SI system because all four constituent constants are defined exactly. $\square$

Connections Master

• Bose-Einstein distribution 11.05.01. The Planck spectrum is the $\mu =0$ special case of the Bose-Einstein distribution applied to photons in a cavity. The non-conservation of photon number forces $\mu =0$, and the spectrum follows by combining the Bose-Einstein occupation with the electromagnetic mode density $g(\omega )=V{\omega }^2/({\pi }^2{c}^3)$. Massive Bose gases at fixed particle number have $\mu <0$ and exhibit Bose-Einstein condensation below a critical temperature, behaviour that has no analogue for photons.

• Bosonic Fock space and second quantisation 12.13.01. The photon-gas Hilbert space is the symmetric Fock space over single-photon modes, with creation operators $a_k^{†}$ obeying canonical commutation relations $[a_k,a_{k^{\prime }}^{†}]={\delta }_{k,k^{\prime }}$. The number operator $N_k=a_k^{†}{a}_k$ has expectation $⟨N_k⟩=1/(e^{\beta \hslash {\omega }_k}-1)$ in the thermal state, providing the second-quantised framework that underlies the canonical-ensemble derivation of the Planck spectrum.

• CMB physics (pending) 13.08.04 pending. The cosmic microwave background is the cosmological realisation of blackbody radiation: a photon gas thermalised during recombination ($z\approx 1100$) and cooled by adiabatic expansion to $T=2.725\text{K}$ today. The COBE-FIRAS spectrum confirms the Planck form to one part in $10^5$ — the most precise blackbody ever measured. Anisotropies $\Delta T/T\sim 10^{-5}$ encode the primordial density perturbations.

• Scattering: Thomson, Rayleigh, Mie (pending) 10.07.04 pending. Rayleigh scattering's ${\nu }^4$ frequency dependence comes from a driven harmonic oscillator far below resonance, related to the same mode-counting integrals that govern blackbody emission. The blue colour of the sky is the Rayleigh scattering of the Planck spectrum's blue tail by atmospheric molecules; sunsets are red because longer wavelengths penetrate the atmospheric path with less scattering.

Historical & philosophical context Master

Kirchhoff 1859 ^{[Kirchhoff1859]} established the universality of blackbody radiation by a thought experiment with two bodies enclosed in a cavity at thermal equilibrium: any deviation from $\epsilon (\nu )/\alpha (\nu )=B(\nu ,T)$ would allow heat to flow from cold to hot through radiation, violating the second law. The universal function $B(\nu ,T)$ became the central problem of late-19th-century theoretical physics.

Stefan 1879 ^[Stefan1879] inferred the $T^4$ law empirically from Tyndall's measurements of the radiant heat from a platinum wire at different temperatures, and Boltzmann 1884 ^{[Boltzmann1884]} derived the same law theoretically from classical thermodynamics and Maxwell's stress tensor for electromagnetic radiation — without any quantum hypothesis. Wien 1893 obtained the displacement law ${\lambda }_{\max }T=b$ from adiabatic-invariance arguments applied to a moving-mirror cavity, and Wien 1896 ^[Wien1896] proposed the exponential high-frequency form $u(\nu )\propto {\nu }^3{e}^{-h\nu /k_{B}T}$ that matched 1890s short-wavelength spectroscopy.

Rayleigh 1900 and Jeans 1905 derived the low-frequency limit $u(\nu )\propto {\nu }^2{k}_{B}T/c^3$ from equipartition applied to electromagnetic modes, exposing the ultraviolet catastrophe: the classical theory predicted infinite total radiated energy. The catastrophe set the stage for Planck's 1900 interpolation. Working at the Physikalisch-Technische Reichsanstalt in Berlin, Planck 1900 ^[Planck1900] proposed that energy elements of size $E=nh\nu$ are exchanged between the radiation field and oscillators in the cavity walls, with $h=6.62607015\times 10^{-34}\text{J}\cdot \text{s}$ (the modern value, now fixed as exact in the 2019 SI redefinition). Planck himself regarded the quantum hypothesis as a mathematical device and spent two decades trying to derive his formula classically.

Einstein 1905 ^{[Einstein1905]} argued that the quantum hypothesis is physically real, not just formal: light itself consists of discrete energy packets, and this interpretation explained the photoelectric effect (the cutoff frequency below which no photoelectrons are emitted, regardless of intensity). The 1921 Nobel Prize was awarded for this argument, not for relativity. Einstein 1916 ^{[Einstein1916]} returned to blackbody radiation with the detailed-balance argument, deriving the $A/B=8\pi h{\nu }^3/c^3$ relation that underlies laser physics and quantum optics. Bose's 1924 derivation of Planck's law from indistinguishable-particle counting, communicated through Einstein, completed the statistical-mechanics framework.

The cosmic microwave background was discovered accidentally by Penzias and Wilson 1965 ^{[PenziasWilson1965]} as unexplained antenna noise at 4080 MHz at the Bell Labs Holmdel horn antenna; Robert Dicke's group at Princeton had been preparing to search for exactly this signal as the relic of an early hot universe. The COBE satellite's FIRAS instrument (Mather et al. 1994 ^[Mather1994]) measured the CMB spectrum to extraordinary precision, finding agreement with the Planck form to better than one part in $10^5$ across 60–600 GHz. The CMB is the universe's most perfect blackbody, and Planck's 1900 hypothesis remains valid across $T$ scales spanning thirteen orders of magnitude, from millikelvin laboratory radiation to the inferred temperatures of Big Bang nucleosynthesis at $T\approx 10^9\text{K}$.

The conceptual lineage of the blackbody problem illustrates the pre-quantum understanding of statistical physics. Kirchhoff identified the universality of $B(\nu ,T)$ in 1859 from a thermodynamic argument that contained no statistical-mechanical content; the universal function had to exist, but its form had to be measured. Stefan extracted the $T^4$ scaling from a noisy data set covering only one decade in temperature; Boltzmann's theoretical derivation supplied the explanation in terms of radiation pressure, but neither author could compute the proportionality constant. Wien's adiabatic-invariance argument fixed the spectral-displacement law without quantisation, and his 1896 high-frequency exponential matched the visible spectrum to within experimental error, leading many physicists in the 1890s to believe the blackbody problem was essentially solved. The Lummer-Pringsheim-Rubens-Kurlbaum measurements at the PTR in 1899–1900 destroyed this confidence by demonstrating clear deviations from Wien at long wavelength; it was these deviations that Planck's interpolation was constructed to fit.

Planck's October 1900 paper proposed the spectral formula as an empirical interpolation; his December 1900 follow-up gave a derivation based on the entropy of $N$ oscillators sharing $P$ energy elements of size $h\nu$, with $h$ fitted to the spectrum. Planck regarded the discreteness as a mathematical fiction enabling the combinatorial counting and resisted any literal interpretation. The transition from Planck's reluctant introduction of $h$ to Einstein's bold assertion that light itself is quantised took five years (Einstein 1905), with the photoelectric-effect cutoff as the empirical evidence. Even Einstein recognised in 1905 that the photon hypothesis seemed to contradict Maxwell's electromagnetism's continuum description of light; it took until the 1920s for the unified picture (wave-particle complementarity, second quantisation) to emerge.

The blackbody calculation also revealed the inadequacy of classical equipartition. The Rayleigh-Jeans law follows directly from assigning $k_{B}T$ to each electromagnetic mode in the cavity; the ultraviolet divergence is a fundamental classical pathology, not a calculation error. Lord Rayleigh himself noted in 1900 that the formula must fail at high frequency for some unknown reason, and Jeans confirmed the derivation in 1905, but neither could supply the cutoff mechanism. Planck's quantum hypothesis was the missing ingredient. The historical irony is that the resolution required abandoning the very principle of equipartition that had unified statistical mechanics in the 1870s — quantum mechanics arrived not as a refinement of classical mechanics, but as a structural rupture motivated by a single experimental fact: the finite total radiation density of a hot cavity.

Modern applications and extensions Master

The blackbody framework underpins much of modern physics, astrophysics, and engineering. We outline several quantitatively important extensions.

Astrophysical thermometry. Stellar spectral classification (Hertzsprung-Russell diagram) uses the Planck-spectrum fit to assign effective temperatures: the O-type stars at $T_{\mathrm{eff}}\approx 40000\text{K}$ peak in the ultraviolet; M-type dwarfs at $T_{\mathrm{eff}}\approx 3000\text{K}$ peak in the near infrared. The Wien displacement law converts the colour index (the magnitude difference between two photometric filters) directly to a temperature. For solar-type stars the Planck fit is accurate to within a few percent across the visible band; deviations from blackbody behaviour (absorption lines, limb darkening, chromospheric continuum) are second-order corrections analysed via radiative-transfer theory. Pyrometers used in industrial furnace monitoring apply the same principle in reverse: they measure infrared flux at two wavelengths and use the Planck-spectrum ratio to infer temperature without thermometer contact, a method essential for furnaces operating above $1000{}^{\circ }$C where thermocouples melt.

Cosmic microwave background. The CMB temperature $T=2.72548\pm 0.00057\text{K}$ and its anisotropy spectrum $C_{\ell }$ provide the most stringent test of the hot Big Bang model. The dipole anisotropy at $\ell =1$ ($\Delta T/T\approx 1.2\times 10^{-3}$) measures the Earth's motion relative to the CMB rest frame ($v\approx 369\text{km/s}$). The acoustic peaks at $\ell \approx 200,540,800,\dots$ encode the baryon density, dark matter density, and curvature of the universe. The CMB photon-to-baryon ratio $\eta =n_{\gamma }/n_B\approx 6\times 10^{-10}$ is computed from Proposition 3 above (with $T=2.725\text{K}$ giving $n_{\gamma }\approx 4.1\times 10^8{\text{m}}^{-3}$) divided by the baryon density inferred from Big Bang nucleosynthesis (which depends on $\eta$); the self-consistent solution constrains the baryon density of the universe to better than 1%. The polarisation of the CMB (B-modes from primordial gravitational waves, E-modes from acoustic perturbations) extends the blackbody picture to a full electromagnetic-field characterisation of the radiation.

Greenhouse effect and planetary energy budgets. The Earth's effective radiating temperature $T_{\mathrm{eff}}=255\text{K}$ from the energy-balance Stefan-Boltzmann calculation is 33 K colder than the actual surface temperature $T_s=288\text{K}$. The difference is the greenhouse effect: the atmosphere is largely opaque to thermal infrared (the molecular vibration-rotation bands of $H_{2}O$, $C{O}_2$, $C{H}_4$), so outgoing radiation escapes from a "radiating layer" several kilometres above the surface. Convective coupling fixes the lapse rate, and the lower atmosphere sits at higher temperature than the radiating layer. Doubling atmospheric $C{O}_2$ raises the effective altitude of the radiating layer; the lapse rate then forces a surface-temperature increase of order $1.5-4\text{K}$ (the climate sensitivity), the central quantitative parameter of climate physics. The Stefan-Boltzmann law also fixes the maximum greenhouse warming: above some atmospheric opacity, the surface temperature is limited by convective instability rather than radiation.

Quantum optics and cavity QED. The Einstein A and B coefficients (Theorem 4 above) are the foundation of laser physics. Light amplification by stimulated emission of radiation requires the upper-state population $N_2$ to exceed the lower-state population $N_1$ weighted by their degeneracies — population inversion. The threshold gain is set by the cavity loss rate; once it is exceeded, the stimulated rate $Bu$ dominates the spontaneous rate $A$, producing coherent monochromatic emission. The same coefficients govern photon counting in cavity quantum electrodynamics, where single atoms in high-finesse optical cavities can be observed emitting and reabsorbing individual photons, with rates set by the cavity mode density (which differs from the free-space $8\pi h{\nu }^3/c^3$ due to the discrete cavity-mode structure — the Purcell effect, a direct experimental test of the mode-density origin of the Einstein A coefficient).

Sub-Kelvin precision thermometry. Modern low-temperature physics relies on blackbody radiation as a thermometric reference. A microwave cavity at $T=100\text{mK}$ radiates with peak frequency ${\nu }_{\max }\approx 5.9\text{GHz}$ and total power $\sigma T^4\approx 5.7\times 10^{-12}\text{W}/{\text{m}}^2$ — measurable with superconducting bolometers cooled to dilution-refrigerator temperatures. Such measurements provide an absolute primary thermometer, traceable directly to Planck's constant via the spectral density formula. The Boltzmann redefinition of the kelvin in 2019 (fixing $k_B$ exactly) made the Stefan-Boltzmann constant exact, and primary radiation thermometry is one of the realisations of the kelvin in the modern SI.

Beyond the photon gas. The blackbody-spectrum machinery extends to any gas of non-interacting massless bosons. Phonons (lattice vibrations) obey the same spectrum at low frequency, with $c$ replaced by the sound speed $c_s$ and a Debye cutoff at ${\omega }_D$ where the wavelength matches the lattice spacing. The Debye specific heat $C_V\propto T^3$ at low $T$ is the phononic analogue of Stefan-Boltzmann. Magnons (spin waves in a ferromagnet) have a dispersion $\omega =D{k}^2$ giving $C_V\propto T^{3/2}$ — the Bloch $T^{3/2}$ law. Gluons in the quark-gluon plasma formed at the early universe and reproduced at the LHC also follow approximately blackbody scaling at high temperature, with corrections from the QCD running coupling and confinement. Hawking radiation from a black hole is a blackbody at the Hawking temperature $T_H=\hslash c^3/(8\pi GM{k}_B)$: a stellar-mass black hole radiates at $T_H\approx 10^{-7}\text{K}$, an undetectably cold thermal flux, while a primordial black hole at the present epoch with $M\approx 10^{15}\text{g}$ would radiate at $T_H\approx 10^{12}\text{K}$ and emit gamma rays detectable in current surveys. The blackbody framework, originally a model for hot iron in a Berlin laboratory, now spans regimes from astrophysical thermometry to quantum gravity.

Bibliography Master

Primary literature:

• Kirchhoff, G. (1859). "Über den Zusammenhang zwischen Emission und Absorption von Licht und Wärme." Monatsbericht Akademie der Wissenschaften Berlin, October 1859, 783–787.

• Stefan, J. (1879). "Über die Beziehung zwischen der Wärmestrahlung und der Temperatur." Sitzungsberichte der Mathematisch-naturwissenschaftliche Classe der Kaiserlichen Akademie der Wissenschaften, Wien 79, 391–428.

• Boltzmann, L. (1884). "Ableitung des Stefan'schen Gesetzes betreffend die Abhängigkeit der Wärmestrahlung von der Temperatur aus der electromagnetischen Lichttheorie." Annalen der Physik und Chemie 22, 291–294.

• Wien, W. (1896). "Über die Energievertheilung im Emissionsspectrum eines schwarzen Körpers." Annalen der Physik 58, 662–669.

• Planck, M. (1900). "Zur Theorie des Gesetzes der Energieverteilung im Normalspectrum." Verhandlungen der Deutschen Physikalischen Gesellschaft 2, 237–245.

• Einstein, A. (1905). "Über einen die Erzeugung und Verwandlung des Lichtes betreffenden heuristischen Gesichtspunkt." Annalen der Physik 17, 132–148.

• Einstein, A. (1916). "Strahlungs-Emission und -Absorption nach der Quantentheorie." Verhandlungen der Deutschen Physikalischen Gesellschaft 18, 318–323.

• Bose, S. N. (1924). "Plancks Gesetz und Lichtquantenhypothese." Zeitschrift für Physik 26, 178–181.

• Penzias, A. A. & Wilson, R. W. (1965). "A measurement of excess antenna temperature at 4080 Mc/s." Astrophysical Journal 142, 419–421.

• Mather, J. C. et al. (1994). "Measurement of the cosmic microwave background spectrum by the COBE FIRAS instrument." Astrophysical Journal 420, 439–444.

• Fixsen, D. J. (2009). "The Temperature of the Cosmic Microwave Background." Astrophysical Journal 707, 916–920.

Modern textbook anchors:

• Pathria, R. K. & Beale, P. D. (2021). Statistical Mechanics, 4th ed. (Elsevier), §7.3 "The photon gas and blackbody radiation."

• Landau, L. D. & Lifshitz, E. M. (1980). Statistical Physics, Part 1, 3rd ed. (Course of Theoretical Physics Vol. 5, Pergamon), §63 "Black-body radiation."

• Mandel, L. & Wolf, E. (1995). Optical Coherence and Quantum Optics (Cambridge University Press), Ch. 13 "Coherence properties of blackbody radiation."

• Reif, F. (1965). Fundamentals of Statistical and Thermal Physics (McGraw-Hill), §9.13–9.14 "The photon gas and the Stefan-Boltzmann law."

• Griffiths, D. J. & Schroeter, D. F. (2018). Introduction to Quantum Mechanics, 3rd ed. (Cambridge University Press), §5.4 "Quantum Statistical Mechanics."

• Carroll, B. W. & Ostlie, D. A. (2007). An Introduction to Modern Astrophysics, 2nd ed. (Pearson Addison-Wesley), §3.5 "Blackbody Radiation."

• Eisberg, R. & Resnick, R. (1985). Quantum Physics of Atoms, Molecules, Solids, Nuclei, and Particles, 2nd ed. (Wiley), Ch. 1 "Thermal Radiation and Planck's Postulate."