Lecture 26.

Blackbody Radiation
(Ch. 7)

Two types of bosons:

(a) Composite particles which contain an even

number of fermions. These number of these
particles is conserved if the energy does not
exceed the dissociation energy (~ MeV in the
case of the nucleus).

(b) particles associated with a field, of which the

most important example is the photon. These
particles are not conserved: if the total
energy of the field changes, particles appear
and disappear. We’ll see that the chemical
potential of such particles is zero in
equilibrium, regardless of density.
 Radiation in Equilibrium with Matter
Typically, radiation emitted by a hot body, or from a laser is not in equilibrium: energy
is flowing outwards and must be replenished from some source. The first step towards
understanding of radiation being in equilibrium with matter was made by Kirchhoff,
who considered a cavity filled with radiation, the walls can be regarded as a heat
bath for radiation.
The walls emit and absorb e.-m. waves. In equilibrium, the walls and radiation must
have the same temperature T. The energy of radiation is spread over a range of
frequencies, and we define uS (,T) d as the energy density (per unit volume) of the
radiation with frequencies between  and +d. uS(,T) is the spectral energy
density. The internal energy of the photon gas: 
u T    uS  , T  d
0

In equilibrium, uS (,T) is the same everywhere in the cavity, and is a function of

frequency and temperature only. If the cavity volume increases at T=const, the
internal energy U = u (T) V also increases. The essential difference between the
photon gas and the ideal gas of molecules: for an ideal gas, an isothermal expansion
would conserve the gas energy, whereas for the photon gas, it is the energy density
which is unchanged, the number of photons is not conserved, but proportional to
volume in an isothermal change.
A real surface absorbs only a fraction of the radiation falling on it. The absorptivity 
is a function of  and T; a surface for which ( ) =1 for all frequencies is called a
black body.
 Photons
The electromagnetic field has an infinite number of modes (standing
waves) in the cavity. The black-body radiation field is a superposition of T
plane waves of different frequencies. The characteristic feature of the
radiation is that a mode may be excited only in units of the quantum
of energy h (similar to a harmonic oscillators) :
  ni  1 / 2 h
i

This fact leads to the concept of photons as quanta of the electromagnetic field. The
state of the el.-mag. field is specified by the number n for each of the modes, or, in other
words, by enumerating the number of photons with each frequency.

According to the quantum theory of radiation, photons are massless E ph  h

bosons of spin 1 (in units ħ). They move with the speed of light :
E ph  cp ph
The linearity of Maxwell equations implies that the photons do not E ph 
p ph  h
interact with each other. (Non-linear optical phenomena are c c
observed when a large-intensity radiation interacts with matter).

The mechanism of establishing equilibrium in a photon gas is absorption and emission

of photons by matter. Presence of a small amount of matter is essential for establishing
equilibrium in the photon gas. We’ll treat a system of photons as an ideal photon gas,
and, in particular, we’ll apply the BE statistics to this system.
 Chemical Potential of Photons = 0
The mechanism of establishing equilibrium in a photon gas is absorption F 
and emission of photons by matter. The textbook suggests that N can be    0
  N T ,V
found from the equilibrium condition:

F  Thus, in equilibrium, the chemical

On the other hand,     ph  ph  0
  N T ,V
potential for a photon gas is zero:

However, we cannot use the usual expression for the chemical potential, because one
cannot increase N (i.e., add photons to the system) at constant volume and at the same
time keep the temperature constant: F 
  - does not exist for the photon gas
  N T ,V

G  N  F  F T ,V 
Instead, we can use G  F  PV P    
 V T V
- by increasing the volume at T=const, we proportionally scale F

F - the Gibbs free energy of an G

Thus, GF V 0  ph  0
V equilibrium photon gas is 0 ! N
For  = 0, the BE distribution reduces to the Planck’s distribution:

n ph  f ph  , T  
1 1 Planck’s distribution provides the average

    h 
exp    1 exp    1 number of photons in a single mode of
 B 
k T  k BT  frequency  = /h.
 h
  n h 
The average energy in the mode:  h 
exp    1
 B 
k T

In the classical (h << kBT) limit:   k BT

In order to calculate the average number of photons per small energy interval d, the
average energy of photons per small energy interval d, etc., as well as the total
average number of photons in a photon gas and its total energy, we need to know the
density of states for photons as a function of photon energy.

kz Density of States for Photons

1 4 / 3 k 3 k 3 volume k3
N k    G k  
8    6 2 6 2
kx
Lx Ly Lz

ky dG  3 2
g      cp  ck G   g 3D
  
d 6  2 c  2  2 c 
3 ph 3

extra factor of 2 due

g    g  
3D 3D d 
h 2
h 
2
8 2
 3 3D 8 2
g    3
to two polarizations: d  c 3
ph ph ph
c c
Spectrum of Blackbody Radiation h g   n   d  uS  , T d
The average energy of photons with frequency photon average number
between  and +d (per unit volume): energy of photons

h g( ) n( )
n( )

h g( )
8 2
g ph    3
3D

c
 3
 =
  
- the spectral density of the
8 h 3
us  , T   h g   f    3 black-body radiation
c exp  h   1 (the Plank’s radiation law)

d
u as a function of the energy: u  , T d  u  , T d u  , T   u  , T   u h , T  h
d

8 3
u  , T  
hc 3 exp 
  
  1
 B 
k T

u(,T) - the energy density per unit photon

energy for a photon gas in equilibrium with
a blackbody at temperature T.
 Classical Limit (small f, large ), Rayleigh-Jeans Law

At low frequencies or high temperatures:  h  1 exp  h  1   h

8 h 3 8  2
us  , T   3
- purely classical result (no h), can be
 3 k BT
c exp  h   1 c obtained directly from equipartition

Rayleigh-Jeans Law

This equation predicts the so-called

ultraviolet catastrophe – an infinite
amount of energy being radiated at
high frequencies or short
wavelengths.
 Rayleigh-Jeans Law (cont.)
u as a function of the wavelength:
3
 c
h 
 d hc  8    hc  8 hc
u  , T d  u  , T d   1
 d   u  , T    
2  hc 3 exp  hc   1  2  5 exp  hc   1
k T  k T 
 B   B 

In the classical limit of large :

8 k BT
u  , T  large   1
4
4
 High  limit, Wien’s Displacement Law
At high frequencies/low temperatures:  h  1 exp  h  1  exp  h 

8 h 3
u s  , T    exp   h  Nobel 1911
c3

The maximum of u() shifts toward higher frequencies with increasing temperature. The
position of maximum:  3

 h 
   
du d   B  k T   3x 2 x 3e x 
 const   const   x  2
0
d  h  
d 
 h  
  exp    1  e  1 x

e  1 
 max  2.8
k BT  k BT    k BT  
h
3  x e x  3  x  2.8

h max Wien’s displacement law

 2 .8 - discovered experimentally
k BT
by Wilhelm Wien
u(,T)

- the “most likely” frequency of a photon in a

blackbody radiation with temperature T
Numerous applications
(e.g., non-contact radiation thermometry)

 max  max
h max
u , T  u , T   2.8 - does this mean that
k BT
hc
 2.8 ? No!
 max max k BT max
3
 c
h 
 d hc  8    hc  8 hc
u  , T d  u  , T d   1
 d   u  , T    
2  hc 3 exp  hc   1  2  5 exp  hc   1
k T  k T 
 B   B 
du d 
 const   5
1  
 const   6
5 
 5

 x 2 exp 1 / x  
0
 2
df dx  x exp 1 / x   1  x exp 1 / x   1 x exp 1 / x   1 

5xexp 1/ x 1  exp 1/ x

hc
 max 
5 k BT

T = 300 K  max  10 m

“night vision” devices

 Solar Radiation
The surface temperature of the Sun - 5,800K.

hc
max   0.5 m
5 k BT

As a function of energy, the spectrum of umax  h max  2.8k BT  1.4 eV

sunlight peaks at a photon energy of
- close to the energy gap in Si, ~1.1 eV, which has been so far the best material for solar cells

Spectral sensitivity of human eye:

 Stefan-Boltzmann Law of Radiation
The total number of photons per unit volume:
  3
8 2 8
3
 k BT  x 2 dx  kB  3
n    n  g  d  3
N
V 0 c 0  h 
d 
c3

 h 
 0 e x  1  hc  T  2.4
 8
exp    1
 k BT  - increases as T 3

The total energy of photons per unit volume :   g  


8 5 k BT 
4
uT    
U
d 
(the energy density of a photon gas) V 0 exp    1 15hc 
3

2 5 k B
4
the Stefan-Boltzmann 4 4 the Stefan-
 constant u T   T Boltzmann Law
15h 3c 2 c

uT  8 5 k BT  hc  4
4 3
The average energy per photon:    k BT  2.7 k BT
15hc  8 k BT   2.4 15  2.4
3 3
N
(just slightly less than the “most” probable energy)

Some numbers: The value of the Stefan-Boltzmann constant:   5.76 10 W / K m

8 4 2
 
Consider a black body at 310K. The power emitted by the body:  T 4  500 W / m2
While the emissivity of skin is considerably less than 1, it still emits a considerable
power in the infrared range. For example, this radiation is easily detectable by modern
techniques (night vision).
 Power Emitted by a Black Body

For the “uni-directional” motion, the flux of energy per unit area  c u T

energy density u

1m2
c  1s

Integration over all angles 1

power emitted by unit area  cu
provides a factor of ¼: 4
(the hole size must be >> the wavelength)

Thus, the power emitted by a unit-area c 4 4

power  u T   
c
surface at temperature T in all directions:
T  T4
4 4 c

The total power emitted by a black-body sphere of radius R:  4R 2 T 4

 Sun’s Mass Loss
The spectrum of the Sun radiation is close to the black body spectrum with the
maximum at a wavelength  = 0.5 m. Find the mass loss for the Sun in one second.
How long it takes for the Sun to loose 1% of its mass due to radiation? Radius of the
Sun: 7·108 m, mass - 2 ·1030 kg.
hc hc  6.6 1034  3 108 
max = 0.5 m  max   T   K  5,740 K 
5 k BT 5 k B max  23
 5 1.38 10  0.5 10
6


2 5 k B
4

P power emitted by a sphere   4R 2 T 4

W
  5.7  10 8

15h 3c 2 m2K 4

This result is consistent with the flux of the solar radiation energy received by the Earth
(1370 W/m2) being multiplied by the area of a sphere with radius 1.5·1011 m (Sun-Earth
distance).
4
 
P  4 RSun   
2 hc
 2

  4 7 108 m  5.7 10 8 2 4  5,740K 4  3.8 10 26 W
W
 2.8 k B max  mK

dm P 3.8 1026 W
the mass loss per one second  2   4.2  10 9
kg/s
dt c 3 10 m
8
2

0.01M 2 1028 kg
1% of Sun’s mass will be lost in t    4.7 1018 s  1.5 1011 yr
dm / dt 4.2 10 kg/s
9
Dewar
Radiative Energy Transfer
Liquid nitrogen and helium are stored in a vacuum or Dewar flask, a
container surrounded by a thin evacuated jacket. While the thermal
conductivity of gas at very low pressure is small, energy can still be
transferred by radiation. Both surfaces, cold and warm, radiate at a rate:

J rad  1  r  Ti
4 i=a for the outer (hot) wall, i=b for the inner (cold) wall,
W / m2
r – the coefficient of reflection, (1-r) – the coefficient of emission

Let the total ingoing flux be J,

J  1  r  Ta4  rJ  J   1  r  Tb4  rJ
and the total outgoing flux be J’:

1 r
The net ingoing flux: J  J 
1 r

 Ta4  Tb4 
If r=0.98 (walls are covered with silver mirror), the net flux is reduced to 1% of the
value it would have if the surfaces were black bodies (r=0).
 Superinsulation
Two parallel black planes are at the temperatures T1 and T2 respectively. The
energy flux between these planes in vacuum is due to the blackbody
radiation. A third black plane is inserted between the other two and is allowed
to come to an equilibrium temperature T3. Find T3 , and show that the energy
flux between planes 1 and 2 is cut in half because of the presence of the third
plane. T1 T3 T2

Without the third plane, the energy flux per unit area is: 
J 0   T1  T2
4 4

The equilibrium temperature of the third plane can be found from the energy balance:
1/ 4

  
 T14  T3 4   T3 4  T2 4  T1  T2  2T3
4 4 4  T14  T2 4 
T3   

 2 
The energy flux between the 1st and 2nd planes in the presence of the third plane:


J   T1  T3
4 4
  4 T14  T2 4  1
   T1   1
   T14  T2 4  J 0
 2
- cut in half
 2  2
Superinsulation: many layers of aluminized Mylar foil loosely
wrapped around the helium bath (in a vacuum space between the
walls of a LHe cryostat). The energy flux reduction for N heat
shields: J0
JN 
N 1
 The Greenhouse Effect
Absorption:
 
2
4 R 
Power in    RE  TSun   Sun 
2

 Rorbit 
the flux of the solar radiation energy
received by the Earth ~ 1370 W/m2
Power out  4 RE  TE
2 4
Emission:

1/ 4
  R  2 
TE    Sun   TSun
 4  Rorbit  

Rorbit = 1.5·1011 m RSun = 7·108 m

Transmittance of the Earth atmosphere

 = 1 – TEarth = 280K
In reality  = 0.7 – TEarth = 256K
To maintain a comfortable temperature on the
Earth, we need the Greenhouse Effect !

However, too much of the greenhouse effect

leads to global warming:
 Thermodynamic Functions of Blackbody Radiation
 uT  16
The heat capacity of a photon gas at constant volume: cV     VT 3
 T V c

This equation holds for all T (it agrees with c T  16V 16
T T
S T    V
c 0
the Nernst theorem), and we can integrate it dT   T  2
dT   VT 3

to get the entropy of a photon gas: 0

T 3c

Now we can derive all thermodynamic functions of blackbody radiation:

4 16 4
the Helmholtz free energy: F  U  TS  VT 4  VT 4   VT 4
c 3c 3c
the Gibbs free energy: G  U  TS  PV  F  PV  0  N 
 F  4 4 1 1
the pressure of a photon gas P     T  u PV  U
(radiation pressure)  V T , N 3c 3 3

For comparison, for a non-relativistic monatomic gas – PV = (2/3)U. The difference

– because the energy-momentum relationship for photons is ultra-relativistic, and
the number of photon depends on T.

1 1 1
In terms of the average
PV  U  V  n    V  n  2.7k BT  0.9 N k BT
density of phonons: 3 3 3
 Radiation in the Universe
Approximately 98% of all the photons
emitted since the Big Bang are observed
now in the submillimeter/THz range.

The dependence of the radiated energy

In the spectrum of the Milky Way galaxy, versus wavelength illustrates the main
at least one-half of the luminous power is sources of the THz radiation: the
emitted at sub-mm wavelengths interstellar dust, emission from light and
heavy molecules, and the 2.7-K cosmic
background radiation.
Cosmic Microwave Background

A. Penzias
R. Wilson
Nobel 1978

In the standard Big Bang model, the radiation

is decoupled from the matter in the Universe
about 300,000 years after the Big Bang,
when the temperature dropped to the point
where neutral atoms form (T~3000K). At this
moment, the Universe became transparent
for the “primordial” photons. The further
expansion of the Universe can be
considered as quasistatic adiabatic
(isentropic) for the radiation:

16
S T   VT 3  const
3c
Since V  R3, the isentropic expansion
leads to :
T  R 1
 “… for their
CMBR (cont.) discovery of the
blackbody form and
anisotropy of the
CMBR”.
Mather, Smoot, Nobel 2006

At present, the temperature of the Planck’s

distribution for the CMBR photons is 2.735
K. The radiation is coming from all
directions and is quite distinct from the
radiation from stars and galaxies.

Alternatively, the later evolution of the radiation temperature may be considered as a

result of the red (Doppler) shift (z). Since the CMBR photons were radiated at T~3000K,
the red shift z~1000.
 Problem 2006 (blackbody radiation)
The cosmic microwave background radiation (CMBR) has a temperature of
approximately 2.7 K.
(a) (5) What wavelength λmax (in m) corresponds to the maximum spectral density
u(λ,T) of the cosmic background radiation?
(b) (5) What frequency max (in Hz) corresponds to the maximum spectral density
u(,T) of the cosmic background radiation?
(c) (5) Do the maxima u(λ,T) and u(,T) correspond to the same photon energy? If
not, why?

hc 6.6 1034  3 108 hc

 1.1 meV
(a) max    1.1103 m  1.1mm
 23
5 k BT 5 1.38 10  2.7 max

k BT 2.8 1.38 1023  2.7 h max  0.65 meV

(b)  max  2.8  34
 1.58 1011 Hz
h 6.6 10

(c) the maxima u(λ,T) and u(,T) do not correspond to the same photon energy. The
reason of that is

 d c 8 hc
u  , T d  u  , T d   1
 d   u  , T 
2  5  hc 
exp    1

 B 
k T
 Problem 2006 (blackbody radiation)
(d) (15) What is approximately the number of CMBR photons hitting the earth per
second per square meter [i.e. photons/(s·m2)]?

(d)
4
 
J   TCMBR  5.7 108 W / K 4  m2  2.7 K 4  3 106W / m2
4

uT  8 5 k BT  hc  4
4 3
   k BT  2.7 k BT
15hc  8 k BT   2.4 15  2.4
3 3
N

W 
J 2 
 photons   m  3 106 16 photons
N     3  10
 s  m2   J  2.7 1.38 1023  2.7 s  m2