Lean, informalized

Nielsen & Chuang · Chapter 7, Quantum computers: physical realization · p. 316

Exercise 7.29 - Spontaneous emission of an excited atom into the modes of free space

An atom’s decay probability grows linearly in time.

Hover over a blue or green term in the text, or a node in the graph, to find it in the other. Select it to read its declaration both in Lean and in words.

The Exercise

An atom with a transition of frequency ω0\omega_0 and dipole moment μ\mu is prepared in its excited state and left alone in free space. Every mode of the electromagnetic field is coupled to it, and the atom can decay by emitting a photon into any one of them. The exercise asks for the probability P(t)P(t) that the atom has decayed by time tt, summed over all the modes, and for the rate at which this probability grows.

In the library, P(t)P(t) is assembled from the probability of decay into a single mode, and the claim gives P(t)P(t) and its time derivative in closed form.

Decay Into One Mode

Besides ω0\omega_0 and μ\mu, the squared coupling depends on the reduced Planck constant ℏ\hbar, the permittivity of the vacuum ε0\varepsilon_0 and the speed of light cc. The squared coupling g2(ω)g^2(\omega) of the atom to a field mode of frequency ω\omega is g2(ω)=ω02μ22ℏωε0c2.g^2(\omega) = \frac{\omega_0^2\,\mu^2}{2\hbar\,\omega\,\varepsilon_0\,c^2}.(1)✦

To lowest order in the coupling, the decay probability pdecay(g2,ω,ω0,t)p_{\mathrm{decay}}(g^2,\omega,\omega_0,t) of an atom of transition frequency ω0\omega_0, coupled with squared strength g2g^2 to a single mode of frequency ω\omega initially empty, is the probability that it has decayed into that mode by time tt: pdecay(g2,ω,ω0,t)=g24sin⁡2((ω−ω0)t/2)(ω−ω0)2.p_{\mathrm{decay}}(g^2,\omega,\omega_0,t) = g^2\,\frac{4\sin^2\!\big((\omega-\omega_0)\,t/2\big)}{(\omega-\omega_0)^2}.(2)

Decay Into Free Space

The total decay probability P(t)P(t) of the atom at time tt is the decay probability into one mode, with the coupling taken at resonance, integrated over all frequencies ω\omega: P(t)=1(2πc)38π3∫−∞∞ω02pdecay(g2(ω0),ω,ω0,t)dω.P(t) = \frac{1}{(2\pi c)^3}\,\frac{8\pi}{3} \int_{-\infty}^{\infty} \omega_0^2\; p_{\mathrm{decay}}\big(g^2(\omega_0),\omega,\omega_0,t\big)\,d\omega.(3)✦

Carrying out the integral gives a constant, the atom’s spontaneous emission rate γrad\gamma_{\mathrm{rad}}: γrad=ω03μ23πℏε0c5.\gamma_{\mathrm{rad}} = \frac{\omega_0^3\,\mu^2}{3\pi\,\hbar\,\varepsilon_0\,c^5}.(4)

The Claim

Theorem 1. Let ω0\omega_0, ℏ\hbar, ε0\varepsilon_0 and cc be nonzero real numbers and μ\mu any real number. At every time t≥0t \ge 0 the total decay probability is P(t)=γradt,P(t) = \gamma_{\mathrm{rad}}\,t,(5)✦ and at every time t>0t > 0 its derivative is dPdt(t)=γrad.\frac{dP}{dt}(t) = \gamma_{\mathrm{rad}}.(6)✦

What Is Checked

The claim is stated over the total decay probability PP of (3), a real function of time once the five parameters are fixed, and over the rate γrad\gamma_{\mathrm{rad}} of (4); PP in turn rests on the coupling (1) and the decay probability into one mode (2).