A1: Deriving Planck's Distribution Law

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\id}{\mathrm{id}}$$ $$\newcommand{\Span}{\mathrm{span}}$$

( \newcommand{\kernel}{\mathrm{null}\,}\) $$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$ $$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$ $$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\id}{\mathrm{id}}$$

$$\newcommand{\Span}{\mathrm{span}}$$

$$\newcommand{\kernel}{\mathrm{null}\,}$$

$$\newcommand{\range}{\mathrm{range}\,}$$

$$\newcommand{\RealPart}{\mathrm{Re}}$$

$$\newcommand{\ImaginaryPart}{\mathrm{Im}}$$

$$\newcommand{\Argument}{\mathrm{Arg}}$$

$$\newcommand{\norm}[1]{\| #1 \|}$$

$$\newcommand{\inner}[2]{\langle #1, #2 \rangle}$$

$$\newcommand{\Span}{\mathrm{span}}$$ $$\newcommand{\AA}{\unicode[.8,0]{x212B}}$$

$$\newcommand{\vectorA}[1]{\vec{#1}} % arrow$$

$$\newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow$$

$$\newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vectorC}[1]{\textbf{#1}}$$

$$\newcommand{\vectorD}[1]{\overrightarrow{#1}}$$

$$\newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}}$$

$$\newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}}$$

$$\newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} }$$

$$\newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}}$$

$$\newcommand{\avec}{\mathbf a}$$ $$\newcommand{\bvec}{\mathbf b}$$ $$\newcommand{\cvec}{\mathbf c}$$ $$\newcommand{\dvec}{\mathbf d}$$ $$\newcommand{\dtil}{\widetilde{\mathbf d}}$$ $$\newcommand{\evec}{\mathbf e}$$ $$\newcommand{\fvec}{\mathbf f}$$ $$\newcommand{\nvec}{\mathbf n}$$ $$\newcommand{\pvec}{\mathbf p}$$ $$\newcommand{\qvec}{\mathbf q}$$ $$\newcommand{\svec}{\mathbf s}$$ $$\newcommand{\tvec}{\mathbf t}$$ $$\newcommand{\uvec}{\mathbf u}$$ $$\newcommand{\vvec}{\mathbf v}$$ $$\newcommand{\wvec}{\mathbf w}$$ $$\newcommand{\xvec}{\mathbf x}$$ $$\newcommand{\yvec}{\mathbf y}$$ $$\newcommand{\zvec}{\mathbf z}$$ $$\newcommand{\rvec}{\mathbf r}$$ $$\newcommand{\mvec}{\mathbf m}$$ $$\newcommand{\zerovec}{\mathbf 0}$$ $$\newcommand{\onevec}{\mathbf 1}$$ $$\newcommand{\real}{\mathbb R}$$ $$\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}$$ $$\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}$$ $$\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}$$ $$\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}$$ $$\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}$$ $$\newcommand{\laspan}[1]{\text{Span}\{#1\}}$$ $$\newcommand{\bcal}{\cal B}$$ $$\newcommand{\ccal}{\cal C}$$ $$\newcommand{\scal}{\cal S}$$ $$\newcommand{\wcal}{\cal W}$$ $$\newcommand{\ecal}{\cal E}$$ $$\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}$$ $$\newcommand{\gray}[1]{\color{gray}{#1}}$$ $$\newcommand{\lgray}[1]{\color{lightgray}{#1}}$$ $$\newcommand{\rank}{\operatorname{rank}}$$ $$\newcommand{\row}{\text{Row}}$$ $$\newcommand{\col}{\text{Col}}$$ $$\renewcommand{\row}{\text{Row}}$$ $$\newcommand{\nul}{\text{Nul}}$$ $$\newcommand{\var}{\text{Var}}$$ $$\newcommand{\corr}{\text{corr}}$$ $$\newcommand{\len}[1]{\left|#1\right|}$$ $$\newcommand{\bbar}{\overline{\bvec}}$$ $$\newcommand{\bhat}{\widehat{\bvec}}$$ $$\newcommand{\bperp}{\bvec^\perp}$$ $$\newcommand{\xhat}{\widehat{\xvec}}$$ $$\newcommand{\vhat}{\widehat{\vvec}}$$ $$\newcommand{\uhat}{\widehat{\uvec}}$$ $$\newcommand{\what}{\widehat{\wvec}}$$ $$\newcommand{\Sighat}{\widehat{\Sigma}}$$ $$\newcommand{\lt}{<}$$ $$\newcommand{\gt}{>}$$ $$\newcommand{\amp}{&}$$ $$\definecolor{fillinmathshade}{gray}{0.9}$$

Albert Einstein developed a simple but effective analysis of induced emission and absorption of radiation along with spontaneous emission that can be used to derive the Planck formula for thermal radiation.

Consider two energy levels for the molecules in a material. The lower of the two is denoted as $$E_1$$ and the higher as $$E_2$$. The probability of a transition from level 1 up to level 2 through induced absorption is assumed to be proportional to the energy density per unit frequency interval, ($$du/d \nu$$). Likewise the probability of an induced transition from level 2 down to level 1 is assumed also to be proportional to ($$du/d\nu$$). These two probabilities are taken to be $$B_{12}(du/d\nu)$$ and $$B_{21}(du/d\nu)$$, respectively, where $$B_{12}$$ and $$B_{21}$$ are constants. The probability of a spontaneous emission is assumed to be a constant $$A_{21}$$.

Let $$N_1$$ and $$N_2$$ be the number of molecules in energy states 1 and 2, respectively. For equilibrium the number of transitions from 1 to 2 has to be equal to the number from 2 to 1; i.e.,

$\underbrace{N_1\left[B_{12}\left(\dfrac{du}{d\nu}\right)\right]}_{\text{flow up}} = \underbrace{N_2\left[B_{21}\left(\dfrac{du}{d\nu} \right)+ A_{21} \right]}_{\text{flow down}} \nonumber$

This means that the ratio of the occupancies of the energy levels must be

$\dfrac{N_2}{N_1} = \dfrac{B_{12}\left(\dfrac{du}{d\nu}\right)}{ B_{21}\left(\dfrac{du}{d\nu}\right) + A_{21}} \label{einstein2}$

But the occupancies are given by the Boltzmann distribution as

$N_1 = N_0 \exp \left(− \dfrac{E_1}{kT} \right) \nonumber$

and

$N_2 = N_0 \exp \left(−\dfrac{E_2}{kT} \right) \nonumber$

where $$k$$ is Boltzmann's constant and $$T$$ is absolute temperature. $$N_0$$ is just a constant that is irrelevant for the rest of the analysis.

Thus according to the Boltzmann distribution

$\dfrac{N_2}{N_1} = \exp \left(−\dfrac{E_2−E_1}{kT} \right) \label{boltz2}$

Therefore for radiative equilibrium, Equations \ref{boltz2} and \ref{einstein2} can be set to each other and

$\exp\left(−\dfrac{E_2−E_1}{kT} \right) = \dfrac{B_{12}\left(\dfrac{du}{d\nu}\right)}{ B_{21}\left(\dfrac{du}{d\nu}\right) + A_{21}} \nonumber$

This condition can be solved for $$(du/dν)$$; i.e.,

$\dfrac{du}{d\nu} = \dfrac{A_{21}}{B_{12}\exp \left( \dfrac{E_2−E_1}{kT} \right)−B_{21}} \nonumber$

Consider what happens to the above expression for as $$T \rightarrow \infty$$. It goes to

$\lim _ {T \rightarrow \infty} \dfrac{du}{d\nu} = \dfrac{A_{21}}{B_{12}−B_{21}} \nonumber$

Einstein maintained that $$(du/dν)$$ must go to infinity as $$T$$ goes to infinity. This requires that $$B_{12}$$ be equal to $$B_{21}$$.

Thus

$\dfrac{du}{d\nu} = \dfrac{A_{21}/B_{21}}{\exp \left(\dfrac{E_2−E_1}{kT}\right)−1} \label{eq10}$

Now Planck's assumption is introduced:

$E_2−E_1 = hν \nonumber$

Thus Equation \ref{eq10} becomes

$\dfrac{du}{d\nu} = \dfrac{A_{21}/B_{21}}{\exp \left(\dfrac{hv}{kT}\right)−1} \label{eq11}$

$\dfrac{du}{d\nu}= \dfrac{8πkTν^2}{c^3} \label{RJ}$

The Planck formula must coincide with the Rayleigh-Jeans Law for sufficiently small $$ν$$. Note that the exponent in the denominator of Equation \ref{eq11} can be expanded (via a Taylor expansion):

$\exp\left(\dfrac{hν}{kT}\right) \approx 1 + \dfrac{hν}{kT} \nonumber$

for sufficiently small $$ν$$.

This means that Equation \ref{eq11} simplifies to

$\dfrac{du}{d\nu} = \dfrac{A_{21}/B_{21}}{1 + (hν/kT) −1} = \dfrac{A_{21}/B_{21}}{hν/kT} \nonumber$

and hence

$\dfrac{du}{d\nu}= \left(\dfrac{A_{21}}{B_{21}} \right) \left( \dfrac{kT}{hν}\right) \label{eq20}$

Equating Equations \ref{RJ} and \ref{eq20} for $$(du/dν)$$ gives

$\left(\dfrac{A_{21}}{B_{21}}\right) \left(\dfrac{kT}{hν}\right) = \dfrac{8πkTν^2}{c^3} \nonumber$

which reduces to

$\dfrac{A_{21}}{B_{21}} = \dfrac{8πhν^3}{c^3} \nonumber$

Thus

$\dfrac{du}{d\nu} = \dfrac{8πhν^³}{c^³} \dfrac{1}{\exp(hν/kT)−1} \nonumber$

This is Planck's formula in terms of frequency.

Reference

1. K.D. Möller, Optics, University Science Books, Mill Valley, California, 1988.

This page titled A1: Deriving Planck's Distribution Law is shared under a CC BY 4.0 license and was authored, remixed, and/or curated by Delmar Larsen.