Spontaneous Emission and Einstein Coefficients

Fill a box with two-level atoms and balance spontaneous emission (A), absorption (B_ab ρ) and stimulated emission (B_ba ρ). At thermal equilibrium the Boltzmann ratio N_a/N_b = e^(ħω₀/k_BT) gives ρ(ω₀) = A/(e^(ħω₀/k_BT)B_ab − B_ba); matching Planck's blackbody formula forces B_ab = B_ba and A = (ħω³/π²c³)B_ba, and with B_ba = π|p|²/(3ε₀ħ²) the spontaneous emission rate comes out as A = ω³|p|²/(3πε₀ħc³).

As I mentioned before, spontaneous emission can occur even with no incoming photons: an excited atom couples to the quantized electromagnetic field and releases some of its own excitation energy as a photon. And the scientist who succeeded in writing spontaneous emission down as equations is… ugh… Einstein, apparently.

Alright, let’s get started.

Einstein A & B coefficients

A box packed with atoms, each drawn as a circle with two energy levels

Let’s say a box is packed full of atoms, each one a two-level system. We take the two levels to be nondegenerate (or equally degenerate), with an isotropic thermal radiation field. Here $\rho(\omega)$ is energy density per unit angular frequency.

Two-level diagram: from ψ_b down to ψ_a go spontaneous emission A (red) and stimulated emission B_ba × ρ(ω₀); from ψ_a up to ψ_b goes absorption B_ab × ρ(ω₀)

Inside the box, let the number of atoms in the ψ_a state be N_a, and the number of atoms in the ψ_b state be N_b.

And if I write out what $A$ and those other numbers in the figure above mean,

A: spontaneous emission rate; B_ab ρ(ω₀): absorption rate; B_ba ρ(ω₀): stimulated emission rate. Clearly, absorption or stimulated emission will be proportional to the incoming light’s energy density, the ω₀ energy density at that!

If we use these numbers to write the rate of change of $N_b$,

Rate of change of N_b: dN_b/dt = −N_b A − N_b B_ba ρ(ω₀) + N_a B_ab ρ(ω₀)

That time-dependent population change is measurable and useful too. Here, though, let’s focus on thermal equilibrium.

What we’re curious about is the equilibrium at $t \to \infty$, so the rate of change there will be 0!!!! That’s what equilibrium means!!!!!

So, like, the number $N_b$ will look constant. (It means the system has reached thermal equilibrium at temperature $T$.)

Solving dN_b/dt = 0 gives ρ(ω₀) = A/((N_a/N_b)B_ab − B_ba)

And the energy of the atoms in the box is either $E_a$ or $E_b$, right???? In a system at temperature $T$, an atom’s energy is $E_a$ or $E_b$, and how many atoms sit at each one follows the Boltzmann distribution!!!!

(You can check my Thermal & Statistical Mechanics series for that~)

That is,

Boltzmann distribution: N_a and N_b are proportional to e^(−E_a/k_BT) and e^(−E_b/k_BT), so the number ratio is N_a/N_b = e^(−E_a/k_BT)/e^(−E_b/k_BT) = e^((E_b − E_a)/k_BT) = e^(ħω₀/k_BT)

Plugging that into ρ(ω₀) = A/((N_a/N_b)B_ab − B_ba) gives ρ(ω₀) = A/(e^(ħω₀/k_BT)B_ab − B_ba)

This, apparently, can be connected to Planck’s blackbody radiation formula.

So what’s blackbody radiation?

(You can also check The Rayleigh-Jeans Formula and Blackbody Radiation (Modern Physics I Studied #7) for this.) To put it simply, when we heat something up, light comes out, right???? If you study that light,

Blackbody spectrum from the textbook: energy density against frequency (in 10¹⁴ Hz) at 2000 K, 4000 K and 6000 K, with the visible region marked

it comes out like this, apparently.

The classical Rayleigh–Jeans formula works in the low-frequency limit $\hbar\omega\ll k_BT$, but fails at high frequencies. Planck’s formula covers both limits, ta-da~!

Planck’s formula: ρ(ω) = ħω³/(π²c³(e^(ħω/k_BT) − 1))

That’s what he fixed it to, apparently.

(Thanks to that, the curtain rose on quantum mechanics, too.)

So we can compare this with the blackbody radiation formula!!!!

At thermal equilibrium, the atoms emit and absorb radiation at balanced rates — this is not a directed beam. What we did was write the spectral energy density of that ω₀ light in terms of the Einstein coefficients, right!!!! Let’s compare this formula with Planck’s formula, the one that matches the experiments!!!!!

Planck’s formula set equal to the Einstein-coefficient result: ħω³/(π²c³(e^(ħω/k_BT) − 1)) = A/(e^(ħω/k_BT)B_ab − B_ba)

As you can see from the comparison, we find B_ab = B_ba.

We proved P_a→b = P_b→a before, didn’t we?????? So you might even think it’s obvious!!!! hehe

And $A$ is

A = (ħω³/(π²c³)) B_ab

Earlier, we derived

The transition rate from before (written P_b→a): (π/(3ε₀ħ²))|p|²ρ(ω₀)

didn’t we?? That is,

B_ba = π|p|²/(3ε₀ħ²)

and therefore $A$ is

A = (ħω³/(π²c³))(π|p|²/(3ε₀ħ²)) = ω³|p|²/(3πε₀ħc³)

(What is this, an all-star reunion of the universe’s physical constants? lol)

And that’s how we get the spontaneous emission rate $A$~


Heads-up. The figures here are the original notes, slips and all. A few things to watch for:

  • In $\rho(\omega_0)$, the ratio in the denominator is $\frac{N_a}{N_b}$, not $\frac{N_a}{N_a}$: setting $\frac{dN_b}{dt} = -N_bA - N_bB_{ba}\rho(\omega_0) + N_aB_{ab}\rho(\omega_0) = 0$ gives $\rho(\omega_0) = \frac{A}{\frac{N_a}{N_b}B_{ab} - B_{ba}}$ (the same ratio in the line where $e^{\hbar\omega_0/k_BT}$ is plugged in).
  • In the Boltzmann ratio, $\frac{N_a}{N_b} = e^{\frac{1}{k_BT}(E_b - E_a)}$, with no minus sign in front (the next lines, ending in $e^{\hbar\omega_0/k_BT}$, are right), and each count carries the normalization: $N_a = \frac{N}{Z}e^{-E_a/k_BT}$ and $N_b = \frac{N}{Z}e^{-E_b/k_BT}$, where $Z$ cancels in the ratio.
  • The formula taken from the last post, $\frac{\pi}{3\varepsilon_0\hbar^2}|p|^2\rho(\omega_0)$, is the transition rate $R_{b\to a}$ (per second), not a probability $P_{b\to a}$; $R_{b\to a} = B_{ba}\rho(\omega_0)$ then gives $B_{ba}$.

Comments

Discussion happens via GitHub Discussions. You'll need a GitHub account to comment.