# Some Further Notes on Blackbody Radiation

*Michael Fowler 7/23/08*

## Thermodynamic Derivation of Stefan’s Law

In 1884, Poynting found expressions for the energy density,
energy flow and momentum density and flow in an electromagnetic field: the energy density is _{} the energy flow rate
is given by the Poynting vector, _{}, and the momentum density _{}. (We write these in
MKS, he used cgs.)

For a plane wave, *B* = *E*/*c*, (_{}) so the energy density _{} and the momentum
density _{}. This was all deduced
in the 1880’s from a fairly lengthy (but straightforward, see for example
*E* = *mc*^{2}: the energy density *u* corresponds to a mass density _{}, the momentum density for a plane wave is then _{}. If we think of a
light wave as having a finite total length, it follows that the energy-momentum
relation must be

*E* = *cp*.

(Of course, this follows immediately from special relativity for massless particles, the point here is that it was previously known—in the 1880’s—from electromagnetic theory.)

Now imagine an oven filled with radiation, energy density *u*, and consider the pressure on a
wall. The argument is the same as for a
gas of particles, except that the radiation travels at *c* and has energy-momentum relation *E* = *cp*. Taking the radiation to be perfectly
reflected, and all directions to be equally likely, one finds

_{}

(Imagine first the radiation to be all perpendicular to a
wall: in the steady state, there would be momentum density *u*/2*c* in each direction
(towards the wall and away from it). A
quantity *uc*/2*c* would bounce off the wall each second, for a total momentum
change (hence pressure!) of *u*. In fact, only the perpendicular component of
the radiation (which is going equally in all directions) counts, this gives the
1/3, exactly as for particles in the kinetic theory of gases.)

Another point: can we say the oven walls are perfectly reflecting? Didn’t we say they were black? Pippard points out that two such different ovens must, at the same temperature, have the same radiation distribution inside, or one could put a small pipe from one to the other and build a temperature difference, contradicting the Second Law.

At this point, we’re ready to apply standard thermodynamics to the radiation. Assume the oven filled with radiation is equivalent to a classical gas in a piston (so the volume can be adjusted).

Following Pippard (*Classical
Thermodynamics*), the fundamental equation

_{}

can be written in terms of energy and entropy *densities* *u* = *U*/*V*, *s* = *S*/*V*,
and remembering _{}, as

_{}

But *u* and *s* only depend on temperature, not volume,
so

_{}

and therefore

_{}

Hence

_{}

so

_{}

from which

_{}

completing the thermodynamic derivation of Stefan’s law.

## Thermodynamic Derivation of Wien’s Displacement Law

The idea here is to expand a cubical container isotropically and slowly (so thermal equilibrium is maintained), allowing no heat transfer (adiabatic expansion). Imagining the radiation inside as standing waves in a reflecting box, the wavelengths will scale upwards with the dimensions of the box, so if we track waves of a particular wavelength, the wavelength will change as

_{}.

(this argument is from Wannier’s book).

Now (from the previous section)

_{}

and in the adiabatic expansion *dS* = 0, so

_{}

from which it is easy to check that the equation of the adiabat is

_{}

and since the radiation pressure goes as *T *^{4},

_{}

Suppose now we track waves having frequency in the small
interval (*f*_{1}, *f*_{1}+ *df*_{1}) as the volume increases slowly from *V*_{1} to *V*_{2}. The
temperature will go from *T*_{1} to *T*_{2}, the frequency range
from (*f*_{1}, *f*_{1}+ *df*_{1}) to (*f*_{2}, *f*_{2}+ *df*_{2}) and from *λf* = *c*,

_{}

Now we know from the Stefan-Boltzmann law that the total radiation intensity changes as:

_{}

Wien argued that Boltzmann’s thermodynamic analysis (previous section) could be applied to just the radiation in a small frequency interval, so one can conclude that

_{}

Changing notation: *T*_{1} goes to *T*, *T*_{2} goes to 1, *f*_{1} goes to *f*, we have

_{}

so apart from the overall factor *T *^{3}, the curve has exactly the same shape as a function
of *f*/*T* at all temperatures. This
is exactly what is observed experimentally. In particular, if the radiation is most intense at frequency *f*_{max}, then we find **Wien’s Displacement Law**:

_{}

The curve can equivalently be written

_{}

probably this form inspired Wien’s Radiation Law conjecture, _{}.

## The Factor ¼ for Radiation from a Hole

If one imagines an energy density *ρ*(*f, T*) streaming down a pipe of cross-section area *A* (the size of the hole) at speed *c*, one would find radiation intensity *Ac ρ*(*f, T*). So why the factor of ¼? The energy density inside the oven has waves equally
moving in opposite directions, so only half of it is moving the right way to
get through the hole. Furthermore, it’s
not all coming directly at the hole, it’s hitting it from all inside
directions, and that coming in at an angle *θ* to the vertical sees effectively a smaller
hole, by a factor cos*θ*. Taking account of this, and adding radiation
from all angles, gives the ¼.

## How Planck Related the Energy Density in the Radiation Field to That in a Single Oscillator

He chose the simplest model possible: one linear oscillator
(in the wall of the oven) driven at frequency _{} by a monochromatic
oscillating electric field, the oscillator constantly damped by emitting energy
as radiation.

It was well known from classical electrodynamics that an
accelerating charge lost energy in this way, and it could be included in the
oscillator dynamics as an effective damping force _{}, so the equation of motion is

_{}

_{} being the usual spring
constant. (*Note*: these are Planck’s cgs units. In MKS, we should replace *e*^{2} by _{}, but we won’t bother to add that clutter, since the *e*^{2} terms cancel out in the
final result.)

This can be further simplified by the following
trick: the oscillator is driven at the external frequency _{}, one can (at least for small damping) write

_{},

say. This is
now a familiar equation, an oscillator in a viscous medium with a drag force
proportional to velocity, and is readily solved, most simply by using complex
numbers. The amplitude *A* of the driven oscillations is:

_{}

The oscillator has energy _{} For small damping,
there is a sharp maximum at the natural frequency of the oscillator, *ω*_{0}.

*To summarize*:
the physics of the above model is that the charged oscillator is being driven
by the single-frequency external electromagnetic field, which is therefore
feeding energy into the oscillator, which energy is constantly being drained
out by radiation damping, that is, emission of outgoing radiation at the same
frequency—so in the steady state the oscillator is in equilibrium with the
external field.

Planck’s next step was to place this oscillator in a
real radiation field, so the single frequency driving source is replaced by a
very large number of driving fields at all different frequencies, and with no
overall coherence. (If we take the oven
to have reflective walls, there will be one driving field for each allowed
standing wave, counting also a factor of two for polarization.) But there’s a big simplification: if we take
the damping to be small, the only *significant* driving frequencies are those close to the natural frequency of the oscillator,
where

_{}

Now, going from a single driving field to an actual
radiation field amounts to summing over the modes inside the narrow peak at *ω*_{0}. The weight can be found by integration: it’s

_{}

from which:

_{}

Now this is only the *x*-component
of the radiation field, going to the full field gets rid of the 3 in the
denominator. The energy density in
(radial) frequency in the radiation field is then found by replacing

_{},

giving

_{},

or in terms of frequency in cycles per second,

_{}

where *U*(*f*, *T*)
is the energy of the oscillator, _{}the radiation energy intensity.

## Radiation Damping and Reversibility

It is perhaps worth noting that although the radiation term
is *formally* equivalent to a viscous
damping term, it does not produce “heat” in quite the same way—the energy lost
goes back into the driving system, and the damping constant is not arbitrary,
but related to the driving coupling. Planck believed at first that this absorption of radiation and
reemission was irreversible, in line with the Second Law. But Boltzmann pointed out to him that the
electromagnetic processes were reversible, just as mechanical processes were,
so this couldn’t be quite right.

### Counting the Electromagnetic Standing Wave Modes in the Oven

To make any *quantitative* progress
in analyzing the radiation, we must have a clear picture of the degrees of freedom
of this system: how many oscillators have frequencies in a given energy
range? These oscillations are standing
electromagnetic waves. The waves are
contained in the oven, so the electric field intensity drops rapidly to zero on
approaching and going into the walls, because inside the walls the electric
energy will be rapidly dissipated by currents or polarization. In fact, then, the boundary condition at the
walls is much like that for waves on a string fixed at both ends, where the
wave amplitude goes to zero at the ends.

As a warm up exercise, consider the frequency distribution function of the
possible different modes of vibration (that is, the different degrees of
freedom) of a string stretched between two points a distance *L* apart.

Possible values of the wavelength are:

_{}

so the frequency

_{}

The allowed frequencies are thus equally spaced *c*/2*L* apart. We define the spectral density by
stating that the number of modes between *f* and *f* + Δ*f* is *N*(*f*)Δ*f* for
Δ small compared with the range of frequencies in the system, but large
compared with the spacing between successive allowed frequencies. Evidently, in this simple one-dimensional
example *N*(*f*) is a constant equal to 2*L*/*c*.

The amplitude of oscillation as a function of time has the form:

_{}

This is more conveniently written:

_{}

where

_{}

using _{}

Here *w* is of course the *radial *frequency—how many* radians* per second, rather than cycles per
second, the time dependent component of the wave moves through. Analogously, *k*, called the *wave
number*, counts the radians of spatial change in unit length (one meter in
SI, but note that practicing physicists often still use CGS).

We are now ready to tackle the more complex problem of three-dimensional standing electromagnetic waves in a cubical oven. The first guess, just generalizing the expression above, would be:

_{}

This will satisfy the wave equation

_{}

if

_{}

Also the electric field will be zero at all the walls if we choose the *k*’s
appropriately, that is, so that sin*k _{x}L* = 0, etc.

**However, this wave is not quite
correct! **

There is no electric charge in the oven away from the walls, so the
divergence of the electric field, _{}, must be identically zero away from the walls, a condition
this wave *does not satisfy*. The
divergence is:

_{}

and, for example, taking any *x* for which sin*k _{x}x* =
0, the second and third terms above are identically zero, but the first term
certainly isn’t.

This problem can be fixed by thinking more carefully about the boundary
conditions at the walls. It is true that
any component of the electric field* parallel* to the wall will be
attenuated rapidly by currents or polarization in the wall. However, there could be an electric field *perpendicular* to the wall, because there could be *surface
charge* on the wall.

This means a possible solution is:

_{}

These components are deliberately chosen to give

_{}

This is telling us something we should have realized from the beginning: in
the standing electromagnetic wave in the cavity, just as in a propagating wave, *the electric field is perpendicular to
the direction in which the wave is moving*. In other words, it is a transverse wave, and in fact there are two
independent polarizations, which we must remember to count when we find the
total number of degrees of freedom.

It remains to find the allowed frequencies of vibration—the normal modes—of
the electromagnetic radiation in the oven. The first step is to find the allowed values of the wave numbers *k _{x}*,

*k*and

_{y}*k*. These are fixed by the boundary conditions sin

_{z}*k*= 0, etc., so

_{x}L_{}

so the vector

_{}

where( *l*, *m*, *n*) are positive integers, like (1, 1, 1)
or (5, 13, 4): in other words, each possible standing electromagnetic wave in
the oven corresponds to a point in the (*k _{x}*,

*k*,

_{y}*k*) space labeled by three positive integers. These are the

_{z}*intersection points of a cubic lattice*in the space

*k*> 0,

_{x}*k*> 0,

_{y}*k*> 0.

_{z}Each such point can be associated with the little lattice cube of volume _{} for which the point in
question is the furthest corner of the little cube from the *k*-space
origin. These little cubes stack
together to fill all that part of *k*-space having *k _{x}* > 0,

*k*> 0 and

_{y}*k*> 0.

_{z}The frequency of vibration of the wave having wave numbers (*k _{x}*,

*k*,

_{y}*k*) is

_{z}_{}

Recall that the experimental measurement of black body radiation from the
oven detects the intensity of radiation in a given *frequency* range. It does* not* tell us the wave numbers
of the radiating modes! Hence, to
compare with these experimental results, we must choose a small frequency range
and find how many possible sets of wave numbers correspond to modes of
vibration having frequencies within that range.

Notice that a fixed value of the frequency corresponds to a spherical
surface in *k*-space:

_{}

Therefore, the number of possible modes of radiation in the oven having
frequencies in the range _{}is equal to the number of lattice-point wave number values _{} between two spherical
surfaces centered at the origin and having radii _{}, and only counting in the positive octant, *k _{x}* > 0,

*k*> 0 and

_{y}*k*> 0.

_{z}Now, as we argues above, each lattice point can be associated with a small cube
of volume_{}. Assuming we choose _{} there will be many of
these small cubes between the spherical surfaces, and the total number of
lattice points between the spheres in the positive octant will be just the
volume of the space between the spheres divided by the volume of one of these
cubes.

So the number of possible modes of radiation in the oven having frequencies
in the range _{} is:

_{}

Putting this in terms of the frequency *f* in cycles per second, _{} and *inserting an extra factor of two* for the
two independent polarizations of the wave (discussed above) we find the
distribution function for modes in frequency is given by:

_{}

where *V* = *L*^{3} is the volume of the oven.