1.2 Blackbody Radiation and Planck’s Law
Historically, Planck’s deciphering of the spectra of incandescent heat and light sources played a key role for the development of quantum mechanics, because it included the first proposal of energy quanta, and it implied that line spectra are a manifestation of energy quantization in atoms and molecules. Planck’s radiation law is also extremely important in astrophysics and in the technology of heat and light sources.
Generically, the heat radiation from an incandescent source is contaminated with radiation reflected from the source. Pure heat radiation can therefore only be observed from a non-reflecting, i.e. perfectly black body. Hence the name blackbody radiation for pure heat radiation. Physicists in the late nineteenth century recognized that the best experimental realization of a black body is a hole in a cavity wall. If the cavity is kept at temperatureT, the hole will emit perfect heat radiation without contamination from any reflected radiation.
Suppose we have a heat radiation source (or thermal emitter) at temperatureT. The power per area radiated from a thermal emitter at temperatureT is denoted as itsexitance(oremittance)e(T ). In the blackbody experimentse(T )ãAis the energy per time leaking through a hole of areaAin a cavity wall.
To calculatee(T )as a function of the temperatureT, as a first step we need to find out how it is related to the densityu(T )of energy stored in the heat radiation.
One half of the radiation will have a velocity component towards the hole, because all the radiation which moves under an angleϑ ≤ π/2 relative to the axis going through the hole will have a velocity componentv(ϑ)=ccosϑin the direction of the hole. To find out the average speedvof the radiation in the direction of the hole, we have to averageccosϑ over the solid angle=2πsr of the forward direction 0≤ϕ≤2π, 0≤ϑ≤π/2:
v= c 2π
2π
0
dϕ π/2
0
dϑ sinϑcosϑ= c
2. (1.3)
The effective energy current density towards the hole is energy density moving in forward direction×average speed in forward direction:
u(T ) 2
c
2 =u(T )c
4, (1.4)
and during the timetan amount of energy
E=u(T )c
4t A (1.5)
will escape through the hole. Therefore the emitted power per areaE/(t A)=e(T ) is
e(T )=u(T )c
4. (1.6)
However, Planck’s radiation law is concerned with thespectral exitancee(f, T ), which is defined in such a way that
e[f1,f2](T )= f2
f1
df e(f, T ) (1.7)
is the power per area emitted in radiation with frequencies f1 ≤ f ≤ f2. In particular, the total exitance is
e(T )=e[0,∞](T )= ∞
0
df e(f, T ). (1.8)
Operationally, the spectral exitance is the power per area emitted with frequencies f ≤f≤f +f, and normalized by the widthf of the frequency interval,
e(f, T )= lim
f→0
e[f,f+f](T )
f = lim
f→0
e[0,f+f]−e[0,f](T )
f = ∂
∂fe[0,f](T ).
The spectral exitancee(f, T )can also be denoted as theemitted power per area and per unit of frequencyor as thespectral exitance in the frequency scale.
The spectral energy densityu(f, T )is defined in the same way. If we measure the energy densityu[f,f+f](T )in radiation with frequency betweenf andf +f, then the energy per volume and per unit of frequency (i.e. the spectral energy density in the frequency scale) is
u(f, T )= lim
f→0
u[f,f+f](T )
f = ∂
∂fu[0,f](T ), (1.9) and the total energy density in radiation is
u(T )= ∞
0
df u(f, T ). (1.10)
The equatione(T ) = u(T )c/4 also applies separately in each frequency interval [f, f+f], and therefore must also hold for the corresponding spectral densities,
e(f, T )=u(f, T )c
4. (1.11)
1.2 Blackbody Radiation and Planck’s Law 5 The following facts were known before Planck’s work in 1900.
• The prediction from classical thermodynamics for the spectral exitancee(f, T ) (Rayleigh-Jeans law) was wrong, and actually non-sensible!
• The exitancee(T )satisfies Stefan’s law1(Stefan, 1879; Boltzmann, 1884)
e(T )=σ T4, (1.12)
with the Stefan–Boltzmann constant
σ =5.6704×10−8 W
m2K4. (1.13)
• The spectral exitancee(λ, T )=e(f, T )
f=c/λãc/λ2per unit of wavelength (i.e.
the spectral exitance in the wavelength scale) has a maximum at a wavelength
λmaxãT =2.898ì10−3mãK=2898μmãK. (1.14) This is Wien’s displacement law (Wien, 1893).
The puzzle was to explain the observed curves e(f, T ) and to explain why classical thermodynamics had failed. We will explore these questions through a calculation of the spectral energy densityu(f, T ). Equation (1.11) then also yields e(f, T ).
The key observation for the calculation ofu(f, T )is to realize thatu(f, T )can be split into two factors. If we want to know the radiation energy densityu[f,f+df] = u(f, T )dfin the small frequency interval[f, f+df], then we can first ask ourselves how many different electromagnetic oscillation modes per volume,(f )df, exist in that frequency interval. Each oscillation mode will then contribute an energy E(f, T )to the radiation energy density, whereE(f, T )is the expectation value of energy in an electromagnetic oscillation mode of frequencyf at temperatureT,
u(f, T )df =(f )dfE(f, T ). (1.15) The spectal energy densityu(f, T )can therefore be calculated in two steps:
1. Calculate the number (f ) of oscillation modes per volume and per unit of frequency (“counting of oscillation modes”).
2. Calculate the mean energy E(f, T ) in an oscillation of frequency f at temperatureT.
The results can then be combined to yield the spectral energy densityu(f, T )= (f )E(f, T ).
1References are enclosed in square brackets, e.g. [1]. For historical context, I have also included parenthetical remarks with the names of scientists and years referring to events preceding the development of quantum mechanics and for the emergence of particular applications.
The number of electromagnetic oscillation modes per volume and per unit of frequency is an important quantity in quantum mechanics and will be calculated explicitly in Chap.12, with the result
(f )= 8πf2
c3 . (1.16)
The corresponding density of oscillation modes in the wavelength scale is (λ)=(f )
f=c/λã c λ2 = 8π
λ4. (1.17)
Statistical physics predicts that the probabilityPT(E)to find an oscillation of energy Ein a system at temperatureT should be exponentially suppressed,
PT(E)= 1 kBT exp
− E kBT
. (1.18)
The possible values of E are not restricted in classical physics, but can vary continuously between 0 ≤ E < ∞. For example, for any classical oscillation with fixed frequencyf, continuously increasing the amplitude yields a continuous increase in energy. The mean energy of an oscillation at temperatureT according to classical thermodynamics is therefore
E
classical= ∞
0
dE EPT(E)= ∞
0
dE E kBT exp
− E kBT
=kBT .
Therefore the spectral energy density in blackbody radiation and the corresponding spectral exitance according to classical thermodynamics should be
u(f, t )=(f )kBT = 8πf2
c3 kBT , e(f, T )=u(f, T )c
4 = 2πf2 c2 kBT ,
(1.19) but this is obviously nonsensical: it would predict that every heat source should emit a diverging amount of energy at high frequencies/short wavelengths! This is theultraviolet catastropheof the Rayleigh-Jeans law.
Max Planck observed in 1900 that he could derive an equation which matches the spectra of heat sources perfectly if he assumes that the energy in electromagnetic waves of frequencyf is quantized in multiples of the frequency,
E=nhf =nhc
λ , n∈N. (1.20)
The exponential suppression of high energy oscillations then reads PT(E)=PT(n)∝exp
−nhf kBT
, (1.21)
1.2 Blackbody Radiation and Planck’s Law 7 but due to the discreteness of theenergy quantahf, the normalized probabilities are now
PT(E)=PT(n)=
1−exp
− hf kBT
exp
−nhf kBT
=exp
−n hf kBT
−exp
−(n+1) hf kBT
, (1.22)
such that∞
n=0PT(n)=1.
The resulting mean energy per oscillation mode is E = ∞
n=0
nhf PT(n)
= ∞
n=0
nhfexp
−n hf kBT
− ∞
n=0
nhfexp
−(n+1) hf kBT
=
∞
n=0
nhfexp
−n hf kBT
−
∞
n=0
(n+1)hfexp
−(n+1) hf kBT
+hf
∞
n=0
exp
−(n+1) hf kBT
. (1.23)
The first two sums cancel, and the last term yields the mean energy in an electromagnetic wave of frequencyf at temperatureT as
E(f, T )=hf exp
−khfBT
1−exp
−khfBT = hf exp
hf kBT
−1
. (1.24)
Combination with(f ) from Eq. (1.16) yields Planck’s formulas for the spectral energy density and spectral exitance in heat radiation,
u(f, T )= 8π hf3 c3
1 exp
hf kBT
−1
, e(f, T )=2π hf3 c2
1 exp
hf kBT
−1
. (1.25)
These functions fitted the observed spectra perfectly! The spectrume(f, T )and the emitted powere[0,f](T )with maximal frequencyf are displayed forT =5780 K in Figs.1.1and1.2.
Fig. 1.1 The spectral emittancee(f, T )for a heat source of temperatureT =5780 K
Fig. 1.2 The emittancee[0,f](T )=f
0dfe(f, T )(i.e. emitted power per area in radiation with maximal frequencyf) for a heat source of temperatureT =5780 K. The asymptote forf → ∞ ise[0,∞](T )≡e(T )=σ T4=6.33×107W/m2for the temperatureT =5780 K