Derivation of Planck's Law
About this page
This page provides a brief derivation of Planck’s law from basic statistical principles. For more information, the reader is referred to the textbook by Rybicki and Lightman (Radiative Processes in Astrophysics, Wiley, 2004) http://books.google.com/books?id=LtdEjNABMlsC&dq=isbn:0471827592&ei=0KPFSOKvE4mIjwGQ2Oj3BA. The reader might also find interest in the historical development of early research in radiation physics as surveyed by Barr http://scitation.aip.org/getabs/servlet/GetabsServlet?prog=normal&id=AJPIAS000028000001000042000001&idtype=cvips&gifs=yes.
Photon gas in a box
First, consider a cubic box with each side of length L that is filled with electromagnetic (EM) radiation (a so-called ‘photon gas’) that forms standing waves whose allowable wavelengths are restricted by the size of the box. We will assume that the waves do not interact and therefore can be separated into the three orthogonal Cartesian directions such that the allowable wavelengths are:
where ni is an integer greater than zero, and i represents one of the three Cartesian directions—x, y, or z.
From quantum mechanics, the energy of a given mode (i.e., an allowable set nx,ny,nz) can be expressed as
where h is Planck’s constant (6.626×10 − 34J s). The number N represents the number of such modes, or photons, of the given energy. Importantly, unlike electrons, an unlimited number of modes, or photons, of a given energy can exist; thus, photons are governed by Bose-Einstein statistics.
Statistical mechanics of the photon gas
To derive the energy density in this photon gas, we first need to know the relative probability with which a given energy state E(N) is occupied at a given temperature. Here, we turn to statistical mechanics, which reveals this probability as
where β is the inverse of thermal energy, or β = (kBT) − 1, and Z(β) is a factor, called the partition function, that normalizes the probability as
where is the energy of a single photon, and the latter equality derives from the relationship between the wavelength λ and the ni indices of the EM waves in the box. This wavelength is related to the speed of light c and frequency ν through the familiar relation
Again from statistical mechanics (and specifically Bose-Einstein statistics), the average energy within a given mode (which is related to the average number of photons N) can be expressed as
Energy density of the photon gas
Now that we have an expression for the average energy of a given mode, we can sum (integrate) over all modes to find the total energy within the photon gas. The total energy can be expressed as an integral over all energies as
where is an important function called the density of states. This function gives the number of allowed modes per unit energy within an interval between and . This function can be derived from the allowable wavelengths and ‘n’ indices as
Thus, the energy per unit volume can be expressed as
where the integrand is the spectral energy density u. This function can be expressed in terms of energy, wavelength, or frequency through the relation such that different forms of u are commonly used. However, they are each integrands in expressions that are used to calculate the overall energy density as
The corresponding expressions for spectral energy density follow:
Blackbody emission intensity
Now assume that a small hole is cut into the box. All radiation emanating from this hole will be moving at the speed of light c. Also, the radiation will be uniformly distributed throughout the hemisphere of solid angles (2π steradians), and one half of the energy will be oriented such that it can move outward through the hole. The spectral radiation intensity is defined as the rate of energy emitted per unit area per unit solid angle and per unit wavelength. The rate of energy emitted per area is simply the product of the energy density derived above and the speed of light (i.e., the distance swept by a ray per unit of time). Therefore, the spectral intensity becomes
Similarly, the spectral intensity (per unit frequency instead of wavelength) is