Blackbody Radiation & Planck's Law

4 hours ~10 min read
Prerequisites

Every warm body glows. Classical physics — Maxwell's waves plus statistical mechanics — makes a clean, parameter-free prediction for the glow's spectrum, and the prediction is not just wrong but infinitely wrong. Planck's 1900 repair, restricting each radiation mode to energies En=nhνE_n = nh\nu, is the first appearance of hh in physics and the opening shot of the quantum revolution. A cat on a radiator has always known thermal physics matters; here we learn exactly how much.

Learning Objectives

After this lesson you will be able to:

  1. Explain why the cavity radiator realizes an ideal blackbody, and state the Stefan–Boltzmann and Wien displacement laws with correct constants.
  2. Derive the Rayleigh–Jeans law from standing-wave mode counting plus equipartition, and explain the ultraviolet catastrophe.
  3. Derive Planck's law from the hypothesis En=nhνE_n = nh\nu via Boltzmann-weighted geometric sums, and recover both limits.
  4. Obtain the Stefan–Boltzmann constant and Wien's constant by integrating and maximizing Planck's law.
  5. Compute blackbody spectra numerically and apply them to real sources (the Sun, the CMB).
  6. State Kirchhoff's law of thermal radiation and use it to explain why the blackbody spectrum depends on temperature alone.

Intuition

Heat a poker: first infrared, then dull red, then orange, then white. Two facts are universal — hotter bodies radiate more (total power grows steeply with TT) and radiate bluer (the peak shifts up in frequency). For an ideal absorber the entire spectrum depends on temperature alone, not on material — so it probes something fundamental about radiation itself.

The classical picture is seductive. The field in a hot cavity is a collection of standing-wave modes — exactly the modes of P.1.1 — and equipartition assigns each mode average energy kBTk_BT. But a cavity supports ever more modes at higher frequency, without bound, so classical physics predicts infinite energy density, diverging in the ultraviolet. Planck's fix: a mode of frequency ν\nu holds energy only in lumps of hνh\nu. High-frequency modes, whose lump dwarfs kBTk_BT, freeze out — the spectrum bends over and the total is finite.


Theory

Light as an electromagnetic wave

Maxwell's equations in vacuum make each field component obey the wave equation of P.1.1, ∇2E=c−2 ∂t2E\nabla^2E = c^{-2}\,\partial_t^2E with c=1/μ0ε0=2.998×108 m s−1c = 1/\sqrt{\mu_0\varepsilon_0} = 2.998\times10^8\ \mathrm{m\,s^{-1}} and λν=c\lambda\nu = c. Light is transverse, with two independent polarizations per propagation direction — remember that factor of 2. The spectrum runs from radio (ν∼106\nu\sim10^6 Hz) through visible (44–7.5×10147.5\times10^{14} Hz, 750–400 nm) to gamma rays. Accelerating charges radiate; thermally agitated charges in matter emit a continuous spectrum: thermal radiation.

Blackbodies and the cavity radiator

The spectral radiancy RT(ν) dνR_T(\nu)\,d\nu is the power emitted per unit area in [ν,ν+dν][\nu,\nu+d\nu]. Real surfaces emit a fraction ϵ≤1\epsilon\le1 (emissivity) of the ideal maximum; a blackbody is the ideal ϵ=1\epsilon=1 absorber-and-emitter. The lab realization is a cavity with a small hole: entering radiation rattles around and is absorbed before escaping, so the hole absorbs everything — it is a blackbody — and its leakage samples the equilibrium energy density u(ν,T)u(\nu,T) inside. A kinetic-theory angle average relates the two: RT(ν)=c4 u(ν,T)R_T(\nu) = \tfrac{c}{4}\,u(\nu,T).

Kirchhoff's law: why the spectrum is universal

Why should the glow depend on temperature alone and not on the material? Kirchhoff (1859) answered before anyone knew the spectrum. Let the spectral irradiance IνI_\nu be the power per unit area and unit frequency falling on a surface, and its spectral absorptance aνa_\nu (0≤aν≤10 \le a_\nu \le 1) the fraction it absorbs. Place any body inside a closed cavity whose walls are held at temperature TT, and wait. At equilibrium the body must emit, in every frequency band, exactly what it absorbs in that band: Rν=aν IνR_\nu = a_\nu\,I_\nu. The balance has to hold band by band, not just in total. Otherwise a filter passing only one band would let the body warm or cool spontaneously against an equal-temperature cavity, a violation of the second law. (The same argument works direction by direction and polarization by polarization.) The irradiance inside the cavity is the same whatever body sits in it, so Kirchhoff's law follows: the ratio Rν/aνR_\nu/a_\nu is a universal function of ν\nu and TT, identical for every material. A blackbody has aν=1a_\nu = 1, so its emission is that universal function, Rνbb(T)=Iν(T)=c4 u(ν,T)R^{\text{bb}}_\nu(T) = I_\nu(T) = \tfrac{c}{4}\,u(\nu,T). Two corollaries:

  • Good absorbers are good emitters: a surface's emissivity at each frequency equals its absorptance, ϵν=aν\epsilon_\nu = a_\nu. Polished silver, which absorbs little, glows feebly when red hot; soot glows brightly.
  • Every real body radiates ϵν\epsilon_\nu times the blackbody spectrum at its temperature, so finding u(ν,T)u(\nu,T) once describes thermal radiation from everything. That universal function is what Rayleigh, Jeans and Planck set out to derive.

The empirical laws

Two facts were nailed down before 1900. Stefan–Boltzmann (1879/84): the total radiancy is

RT=σT4,σ=5.670×10−8 W m−2 K−4. R_T = \sigma T^4, \qquad \sigma = 5.670\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}} .

Wien displacement (1893): the peak of the wavelength spectrum obeys

λmax T=2.898×10−3 m K. \lambda_{max}\,T = 2.898\times10^{-3}\ \mathrm{m\,K} .

A 5800 K surface (the Sun) peaks near 500 nm, mid-visible; your 310 K body peaks near 9.4 μm.

Rayleigh–Jeans: counting modes, feeding them equipartition

Take a cubic cavity of side aa (V=a3V = a^3) with conducting walls, so fields vanish at the walls — the 3D standing waves of P.1.1, E∝sin⁡(nxπx/a)sin⁡(nyπy/a)sin⁡(nzπz/a)E \propto \sin(n_x\pi x/a)\sin(n_y\pi y/a)\sin(n_z\pi z/a) with ki=niπ/ak_i = n_i\pi/a, ni=1,2,3,…n_i = 1,2,3,\dots The dispersion relation ∣k∣=2πν/c|\mathbf k| = 2\pi\nu/c gives

nx2+ny2+nz2=(2aνc) ⁣2≡r2. n_x^2+n_y^2+n_z^2 = \left(\frac{2a\nu}{c}\right)^{\!2} \equiv r^2 .

Modes are lattice points in the positive octant of n\mathbf n-space (sign flips give the same standing wave). Counting points below frequency ν\nu as the octant volume of a sphere of radius rr, times 2 for polarization:

N(<ν)=18⏟octant⋅43π(2aνc) ⁣3⋅2⏟polarizations=8πVν33c3⟹1VdNdν=8πν2c3. N(<\nu) = \underbrace{\tfrac18}_{\text{octant}}\cdot\tfrac43\pi\left(\frac{2a\nu}{c}\right)^{\!3}\cdot\underbrace{2}_{\text{polarizations}} = \frac{8\pi V\nu^3}{3c^3} \quad\Longrightarrow\quad \frac{1}{V}\frac{dN}{d\nu} = \frac{8\pi\nu^2}{c^3} .

Classically each mode is a harmonic degree of freedom in equilibrium, and equipartition (Boltzmann-weighting over a continuum of energies, ⟨E⟩=∫0∞Ee−E/kBTdE/∫0∞e−E/kBTdE\langle E\rangle = \int_0^\infty E e^{-E/k_BT}dE \big/ \int_0^\infty e^{-E/k_BT}dE) assigns it ⟨E⟩=kBT\langle E\rangle = k_BT independent of ν\nu, with kB=1.381×10−23 J K−1k_B = 1.381\times10^{-23}\ \mathrm{J\,K^{-1}}. Modes × energy-per-mode is the Rayleigh–Jeans law:

u(ν,T)=8πν2c3 kBT. u(\nu,T) = \frac{8\pi\nu^2}{c^3}\,k_BT .

It fits experiment at low frequency — and the total energy diverges: ∫0∞u dν∝∫0∞ν2 dν=∞\int_0^\infty u\,d\nu \propto \int_0^\infty \nu^2\,d\nu = \infty. Every oven should contain infinite energy, mostly beyond the ultraviolet: the ultraviolet catastrophe. Nothing here is sloppy — the mode counting survives into quantum theory unchanged, and equipartition is a theorem. Classical physics itself is wrong.

Planck's hypothesis and the quantized average energy

Planck (1900) kept the mode counting and changed the energy assignment: a mode of frequency ν\nu may only hold

En=n hν,n=0,1,2,…,h=6.626×10−34 J s. E_n = n\,h\nu, \qquad n = 0,1,2,\dots, \qquad h = 6.626\times10^{-34}\ \mathrm{J\,s} .

Boltzmann factors still weight each allowed energy, but integrals become sums. With x≡hν/kBTx \equiv h\nu/k_BT and q≡e−xq \equiv e^{-x}, both sums are geometric:

Z=∑n=0∞e−nx=11−e−x,∑n=0∞n e−nx=−dZdx=e−x(1−e−x)2, Z = \sum_{n=0}^{\infty} e^{-nx} = \frac{1}{1-e^{-x}}, \qquad \sum_{n=0}^{\infty} n\,e^{-nx} = -\frac{dZ}{dx} = \frac{e^{-x}}{(1-e^{-x})^2}, ⟨E⟩=∑nnhν e−nx∑ne−nx=hν e−x1−e−x=hνehν/kBT−1. \langle E\rangle = \frac{\sum_n nh\nu\,e^{-nx}}{\sum_n e^{-nx}} = h\nu\,\frac{e^{-x}}{1-e^{-x}} = \frac{h\nu}{e^{h\nu/k_BT}-1} .

Multiplying by the unchanged mode density gives Planck's law:

 u(ν,T)=8πν2c3 hνehν/kBT−1  \boxed{\,u(\nu,T) = \frac{8\pi\nu^2}{c^3}\,\frac{h\nu}{e^{h\nu/k_BT}-1}\,}

Low frequency (hν≪kBTh\nu\ll k_BT): ex−1≈xe^x - 1 \approx x, so ⟨E⟩→kBT\langle E\rangle \to k_BT — Rayleigh–Jeans recovered exactly where it worked. High frequency (hν≫kBTh\nu\gg k_BT): ⟨E⟩≈hν e−hν/kBT\langle E\rangle \approx h\nu\,e^{-h\nu/k_BT}, the exponential Wien form u∝ν3e−hν/kBTu \propto \nu^3 e^{-h\nu/k_BT}. Modes whose quantum exceeds the thermal budget are almost never excited: discreteness starves the ultraviolet.

The demo plots the same physics per unit wavelength, which moves the peak (the caution under Wien displacement below explains why), with the Rayleigh–Jeans curve as an optional overlay:

Interactive demo: Blackbody spectrum. Heat a black body and watch Planck’s curve rise and its peak slide, beside the classical ultraviolet catastrophe.

Caution. Planck quantized the energies of the cavity modes (equivalently, of the material oscillators exchanging energy with them) — he did not claim free light is made of particles. Radiation in flight was still, for Planck, a classical wave. The bolder step — the photon — is Einstein's (1905), next course: P.3.1 Light as Particles.

Stefan–Boltzmann from Planck

Integrate over all frequencies with x=hν/kBTx = h\nu/k_BT, dν=(kBT/h) dxd\nu = (k_BT/h)\,dx:

utot(T)=8πhc3(kBTh) ⁣4∫0∞ ⁣x3ex−1 dx=8π5kB415 h3c3 T4, u_{tot}(T) = \frac{8\pi h}{c^3}\left(\frac{k_BT}{h}\right)^{\!4}\int_0^\infty\!\frac{x^3}{e^x-1}\,dx = \frac{8\pi^5k_B^4}{15\,h^3c^3}\,T^4,

using ∫0∞x3/(ex−1) dx=π4/15\int_0^\infty x^3/(e^x-1)\,dx = \pi^4/15. With RT=(c/4) utotR_T = (c/4)\,u_{tot}:

RT=σT4,σ=2π5kB415 h3c2=5.670×10−8 W m−2 K−4. R_T = \sigma T^4, \qquad \sigma = \frac{2\pi^5k_B^4}{15\,h^3c^2} = 5.670\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}} .

The empirical T4T^4 law now follows from first principles. Historically the logic ran in reverse: Planck fitted his law to the data and extracted the first accurate values of both hh and kBk_B.

Wien displacement from Planck

Maximize u∝x3/(ex−1)u \propto x^3/(e^x-1) over x=hν/kBTx = h\nu/k_BT: setting the derivative to zero, 3x2(ex−1)=x3ex3x^2(e^x-1) = x^3e^x, i.e.

x=3(1−e−x), x = 3\left(1 - e^{-x}\right),

a transcendental equation with nonzero root x≈2.821x \approx 2.821. Hence νmax/T=2.821 kB/h=5.88×1010 Hz K−1\nu_{max}/T = 2.821\,k_B/h = 5.88\times10^{10}\ \mathrm{Hz\,K^{-1}}.

Caution. The peak of the wavelength distribution is not at c/νmaxc/\nu_{max}. Since u(λ)∣dλ∣=u(ν)∣dν∣u(\lambda)|d\lambda| = u(\nu)|d\nu| with Jacobian c/λ2c/\lambda^2, the λ\lambda-form is u(λ,T)=(8πhc/λ5)[ehc/λkBT−1]−1u(\lambda,T) = (8\pi hc/\lambda^5)[e^{hc/\lambda k_BT}-1]^{-1}, whose maximum solves x=5(1−e−x)x = 5(1-e^{-x}), root x≈4.965x \approx 4.965 — giving the familiar λmaxT=hc/(4.965 kB)=2.898×10−3 m K\lambda_{max}T = hc/(4.965\,k_B) = 2.898\times10^{-3}\ \mathrm{m\,K}. Same spectrum, different binning; always say which variable you maximized over.

Looking ahead: a cavity mode is a harmonic oscillator

The ladder En=nhνE_n = nh\nu is no accident: each electromagnetic mode is mathematically a harmonic oscillator, whose exact quantum levels are En=(n+12)hνE_n = (n+\tfrac12)h\nu — the ladder returns, with a zero-point offset, in P.5.3 The Harmonic Oscillator. The same quantized-mode structure, engineered in microwave resonators coupled to artificial atoms, is circuit quantum electrodynamics — the physics of superconducting quantum processors (Term 4.4).


Worked Examples

Example 1 — The Sun as a blackbody

Solar surface: T≈5800T \approx 5800 K, radius R⊙=6.96×108R_\odot = 6.96\times10^{8} m.

Peak: λmax=2.898×10−3/5800=5.00×10−7\lambda_{max} = 2.898\times10^{-3}/5800 = 5.00\times10^{-7} m =500= 500 nm — mid-visible, matching the evolution of eyes. Flux: RT=σT4=5.670×10−8×(5800)4=6.42×107 W m−2R_T = \sigma T^4 = 5.670\times10^{-8}\times(5800)^4 = 6.42\times10^{7}\ \mathrm{W\,m^{-2}}. Luminosity: L=4πR⊙2RT=4π(6.96×108)2×6.42×107=3.9×1026L = 4\pi R_\odot^2R_T = 4\pi(6.96\times10^8)^2\times6.42\times10^7 = 3.9\times10^{26} W — the accepted solar output, from two constants and a temperature.

Example 2 — Why high-frequency modes freeze out

At T=300T = 300 K, kBT=4.14×10−21k_BT = 4.14\times10^{-21} J =0.0259= 0.0259 eV.

Radio mode, ν=1\nu = 1 GHz: hν=6.63×10−25h\nu = 6.63\times10^{-25} J, so x=1.6×10−4x = 1.6\times10^{-4} and ⟨E⟩=hν/(ex−1)≈kBT(1−x/2)=0.99992 kBT\langle E\rangle = h\nu/(e^x-1) \approx k_BT(1-x/2) = 0.99992\,k_BT — equipartition holds to 0.01%; Rayleigh–Jeans is excellent here.

Visible mode, ν=5×1014\nu = 5\times10^{14} Hz: hν=3.31×10−19h\nu = 3.31\times10^{-19} J =2.07= 2.07 eV, so x=80x = 80 and ⟨E⟩≈hν e−80≈6×10−54\langle E\rangle \approx h\nu\,e^{-80} \approx 6\times10^{-54} J — below equipartition by a factor ∼1033\sim10^{33}. The mode exists, but the bath almost never scrapes together its first quantum. That suppression is the resolution of the catastrophe — and why a warm room is dark.


Hands-on (Python)

import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import quad
from scipy.optimize import brentq

h, c, kB = 6.626e-34, 2.998e8, 1.381e-23        # J s, m/s, J/K

def planck(nu, T):
    """Spectral energy density u(nu,T) in J m^-3 Hz^-1."""
    return (8 * np.pi * h * nu**3 / c**3) / np.expm1(h * nu / (kB * T))

def rayleigh_jeans(nu, T):
    return 8 * np.pi * nu**2 * kB * T / c**3

def wien_approx(nu, T):
    return (8 * np.pi * h * nu**3 / c**3) * np.exp(-h * nu / (kB * T))

# --- 1. Planck vs Rayleigh-Jeans vs Wien at T = 5800 K -------------------
T = 5800.0
nu = np.linspace(1e12, 2.5e15, 2000)
plt.plot(nu, planck(nu, T), label="Planck")
plt.plot(nu, rayleigh_jeans(nu, T), "--", label="Rayleigh–Jeans (diverges)")
plt.plot(nu, wien_approx(nu, T), ":", label="Wien approximation")
plt.ylim(0, 1.8 * planck(nu, T).max()); plt.xlabel(r"$\nu$ (Hz)")
plt.ylabel(r"$u(\nu,T)$"); plt.legend(); plt.title("T = 5800 K"); plt.show()
# Expected: RJ hugs Planck at low nu then blows up; Wien matches the high-nu tail.

# --- 2. Verify the T^4 law and extract sigma ------------------------------
for T in [300.0, 1000.0, 3000.0, 5800.0]:
    nu_max = 100 * kB * T / h              # integrand is dead beyond x = 100
    u_tot, _ = quad(planck, 1.0, nu_max, args=(T,))
    print(f"T = {T:6.0f} K   sigma = {(c/4) * u_tot / T**4:.4e}")
# Expected: sigma = 5.676e-08 W m^-2 K^-4 for EVERY T -> R_T = sigma T^4.
# (CODATA: 5.670e-08; the 0.1% offset is rounding in our 4-digit h, kB.)

# --- 3. Wien displacement via brentq, applied to the CMB ------------------
x_nu  = brentq(lambda x: x - 3 * (1 - np.exp(-x)), 0.1, 10)
x_lam = brentq(lambda x: x - 5 * (1 - np.exp(-x)), 0.1, 10)
print(x_nu, x_lam)                                   # 2.8214...  4.9651...
print(f"Wien constant = {h*c/(x_lam*kB):.4e} m K")   # 2.898e-03  (verified)

T_cmb = 2.725                                        # cosmic microwave background
print(f"CMB peak = {x_nu * kB * T_cmb / h / 1e9:.1f} GHz")   # ~160.2 GHz

The CMB — the 2.725 K afterglow of the Big Bang — is the most perfect blackbody ever measured; COBE found deviations from Planck's law below 50 parts per million.


Exercises

E1 (easy). Your body: T=310T = 310 K, area A=1.7 m2A = 1.7\ \mathrm{m^2}, ϵ≈1\epsilon \approx 1 in the infrared. Find λmax\lambda_{max} and the total power you radiate.

Solution

λmax=2.898×10−3/310=9.35 μ\lambda_{max} = 2.898\times10^{-3}/310 = 9.35\ \mum (far infrared — thermal cameras look here). P=σT4A=5.670×10−8×(310)4×1.7≈890P = \sigma T^4A = 5.670\times10^{-8}\times(310)^4\times1.7 \approx 890 W. You also absorb σTenv4A≈690\sigma T_{env}^4A \approx 690 W from a 293 K room; the net ∼\sim150–200 W is your metabolic scale.

E2 (easy). Expand Planck's ⟨E⟩\langle E\rangle for hν≪kBTh\nu \ll k_BT one order beyond equipartition and show ⟨E⟩≈kBT−12hν\langle E\rangle \approx k_BT - \tfrac12h\nu.

Solution

With x≪1x \ll 1: ex−1=x(1+x/2+O(x2))e^x-1 = x(1 + x/2 + O(x^2)), so ⟨E⟩=hνx(1+x/2)−1≈kBT(1−x/2)=kBT−12hν\langle E\rangle = \frac{h\nu}{x}(1+x/2)^{-1} \approx k_BT(1 - x/2) = k_BT - \tfrac12h\nu. Quantization always lowers the average below equipartition, more at higher frequency.

E3 (medium). Reproduce the geometric-sum derivation of ⟨E⟩=hν/(ehν/kBT−1)\langle E\rangle = h\nu/(e^{h\nu/k_BT}-1) without looking, evaluating both sums explicitly.

Solution

Let x=hν/kBTx = h\nu/k_BT, q=e−x<1q = e^{-x} < 1. Partition sum Z=∑n≥0qn=1/(1−q)Z = \sum_{n\ge0}q^n = 1/(1-q); energy sum ∑nnhν qn=hν q ddq11−q=hν q(1−q)2\sum_n nh\nu\,q^n = h\nu\,q\,\frac{d}{dq}\frac{1}{1-q} = h\nu\,\frac{q}{(1-q)^2}. Ratio: ⟨E⟩=hνq1−q=hνex−1\langle E\rangle = h\nu\frac{q}{1-q} = \frac{h\nu}{e^x-1}. Limits: kBTk_BT as x→0x\to0; hνe−xh\nu e^{-x} as x→∞x\to\infty.

E4 (medium). Compute the CMB's total energy density utot=4σcT4u_{tot} = \frac{4\sigma}{c}T^4 at T=2.725T = 2.725 K, in J m−3\mathrm{J\,m^{-3}} and eV cm−3\mathrm{eV\,cm^{-3}}.

Solution

a≡4σ/c=7.566×10−16 J m−3K−4a \equiv 4\sigma/c = 7.566\times10^{-16}\ \mathrm{J\,m^{-3}K^{-4}} and T4=55.1T^4 = 55.1, so utot=4.17×10−14 J m−3=2.60×105 eV m−3=0.26 eV cm−3u_{tot} = 4.17\times10^{-14}\ \mathrm{J\,m^{-3}} = 2.60\times10^{5}\ \mathrm{eV\,m^{-3}} = 0.26\ \mathrm{eV\,cm^{-3}}. Every cubic centimeter of the universe holds about a quarter electron-volt of primordial light.

E5 (hard). Derive the λ\lambda-form of Planck's law from the ν\nu-form, maximize it to get x=5(1−e−x)x = 5(1-e^{-x}), and compute Wien's constant. Then explain in one sentence why λmaxνmax≠c\lambda_{max}\nu_{max} \ne c.

Solution

Equal energy in corresponding bins: u(λ)=u(c/λ) ∣dν/dλ∣=u(c/λ) c/λ2=8πhcλ51ehc/λkBT−1u(\lambda) = u(c/\lambda)\,|d\nu/d\lambda| = u(c/\lambda)\,c/\lambda^2 = \frac{8\pi hc}{\lambda^5}\frac{1}{e^{hc/\lambda k_BT}-1}. With x=hc/λkBTx = hc/\lambda k_BT, maximizing x5/(ex−1)x^5/(e^x-1) gives 5x4(ex−1)=x5ex5x^4(e^x-1) = x^5e^x, i.e. x=5(1−e−x)x = 5(1-e^{-x}), root x≈4.965x \approx 4.965; then λmaxT=hc4.965 kB=(6.626×10−34)(2.998×108)4.965×1.381×10−23=2.898×10−3 m K\lambda_{max}T = \frac{hc}{4.965\,k_B} = \frac{(6.626\times10^{-34})(2.998\times10^{8})}{4.965\times1.381\times10^{-23}} = 2.898\times10^{-3}\ \mathrm{m\,K}. ✓ The two peaks differ because maximizing a density depends on the variable it is a density in — the Jacobian c/λ2c/\lambda^2 reweights the spectrum, so xλ=4.965≠xν=2.821x_\lambda = 4.965 \ne x_\nu = 2.821.

E6 (easy). Stefan's T4T^4 law (1879) began with data published by Tyndall: a platinum wire at 1200 °C radiated about 11.7 times as much as at 525 °C. (a) Show that a T4T^4 law fits, and what it would have predicted had the temperatures been used in Celsius. (b) Platinum is far from black. What does the agreement silently assume, by Kirchhoff's law?

Solution

(a) In kelvin the ratio is 1473/798=1.8461473/798 = 1.846, and 1.8464=11.61.846^4 = 11.6, within 1% of the measured 11.7. In Celsius, (1200/525)4=2.294=27(1200/525)^4 = 2.29^4 = 27, wildly off: only an absolute temperature can enter a power law. (b) A real surface radiates R=ϵσT4R = \epsilon\sigma T^4, so the measured ratio is (ϵ1473/ϵ798)(1473/798)4(\epsilon_{1473}/\epsilon_{798})(1473/798)^4. The fit assumes the platinum's emissivity was the same at both temperatures, which is not true of metals in general. The T4T^4 law was put on firm ground later: Boltzmann derived it from thermodynamics in 1884, and cavity measurements in the 1890s confirmed it for true blackbody radiation.

E7 (easy). Two stars have their spectral peaks (per unit wavelength) at 510 nm and 350 nm. (a) Find their surface temperatures. (b) Find the ratio of the power each radiates per unit area of surface. (c) Can you tell which star is more luminous?

Solution

(a) T=b/λmaxT = b/\lambda_{max}: 2.898×10−3/510×10−9=5.7×1032.898\times10^{-3}/510\times10^{-9} = 5.7\times10^3 K and 2.898×10−3/350×10−9=8.3×1032.898\times10^{-3}/350\times10^{-9} = 8.3\times10^3 K. The bluer star is hotter. (b) RT∝T4R_T \propto T^4: (8280/5682)4=1.4574=4.5(8280/5682)^4 = 1.457^4 = 4.5. (c) No. Luminosity is L=4πR2σT4L = 4\pi R^2\sigma T^4, and the radius is not known: a cool giant easily outshines a hot dwarf. Color gives temperature; it takes a distance or a size as well to get luminosity.


Checkpoint

  1. Why does a small hole in a cavity behave as a perfect blackbody?
  2. Derive the mode density 8πν2/c38\pi\nu^2/c^3, identifying where the octant factor 18\tfrac18 and the polarization factor 2 enter.
  3. Which classical ingredient causes the ultraviolet catastrophe — mode counting or equipartition — and how does Planck's hypothesis repair it?
  4. Show in two lines how Planck's law reduces to Rayleigh–Jeans for hν≪kBTh\nu\ll k_BT and to Wien's exponential for hν≫kBTh\nu\gg k_BT.
  5. What did Planck quantize — and what did he not quantize?
  6. State Kirchhoff's law. Why does it make the blackbody spectrum a universal function of ν\nu and TT?
Answers
  1. Any ray entering the hole reflects many times and is absorbed before escaping, so the hole absorbs essentially everything (ϵ=1\epsilon = 1); in equilibrium it must therefore also emit the ideal blackbody spectrum of the cavity interior.
  2. Standing waves require ki=niπ/ak_i = n_i\pi/a with ni>0n_i > 0; counting lattice points inside the sphere ∣n∣=2aν/c|\mathbf n| = 2a\nu/c uses the positive octant only (factor 18\tfrac18), and each spatial mode has two transverse polarizations (factor 2): N=8πVν3/3c3N = 8\pi V\nu^3/3c^3, whose ν\nu-derivative per volume is 8πν2/c38\pi\nu^2/c^3.
  3. Equipartition. The mode counting survives into QM; the error is giving every mode kBTk_BT. With En=nhνE_n = nh\nu, ⟨E⟩=hν/(ehν/kBT−1)→0\langle E\rangle = h\nu/(e^{h\nu/k_BT}-1) \to 0 exponentially at high ν\nu, so the total energy integral converges.
  4. hν≪kBTh\nu\ll k_BT: ehν/kBT−1≈hν/kBTe^{h\nu/k_BT}-1 \approx h\nu/k_BT, so u→(8πν2/c3)kBTu \to (8\pi\nu^2/c^3)k_BT. hν≫kBTh\nu\gg k_BT: drop the −1-1, so u→(8πhν3/c3)e−hν/kBTu \to (8\pi h\nu^3/c^3)e^{-h\nu/k_BT}.
  5. He quantized the allowed energies of cavity modes / material oscillators (En=nhνE_n = nh\nu) — energy exchange with radiation. He did not quantize free light; the photon is Einstein's step, in P.3.1.
  6. In equilibrium a body emits, band by band, what it absorbs: Rν=aνIνR_\nu = a_\nu I_\nu. The irradiance inside an isothermal cavity is the same whatever body sits there, so Rν/aνR_\nu/a_\nu is the same for all materials. A blackbody (aν=1a_\nu = 1) emits exactly that universal function, and every other body emits a fraction ϵν=aν\epsilon_\nu = a_\nu of it.

Further Reading

  • [ER] Eisberg & Resnick, §1.1–1.6 — thermal radiation, cavity mode counting, and Planck's postulate; the definitive elementary treatment.
  • [ER] Eisberg & Resnick, Appendix to Ch. 1 — the equipartition theorem and its breakdown.
  • [Gri] Griffiths & Schroeter, §5.4 — the photon gas and the blackbody spectrum, revisited with proper quantum statistics.

← Prev: Hamiltonian Mechanics · Up: Pre-Term · Next: Atomic Models & Spectral Series →

Ready to measure your state?

7 exercises · 12 checkpoint questions

Start the quiz