Lifetime, Linewidth, and Natural Broadening

Tune a laser across the sodium D line and you discover something an energy-level diagram never told you: the "line" is not a line. It is a little hill, about ten megahertz wide, with sharp Gaussian shoulders if the vapor is warm and long, lazy wings that refuse to die away. Cool the atoms, cancel the Doppler shifts, and the hill does not collapse to a spike. A floor remains: the natural linewidth, and no experimental cleverness gets under it. This article is about that floor. The physics is one Fourier fact wearing several costumes: an excited state that survives only a time τ\tau radiates a wave train that rings down in a time τ\tau, and a truncated ring-down cannot have a single frequency. We derive the exact consequence — exponential decay transforms into a Lorentzian profile of full width Γ=1/τ\Gamma = 1/\tau — then push the same integral further than the textbooks usually do: we compute what a spectrometer sees when it watches for only half a lifetime, two, ten, and watch the sinc-ridged profile relax onto the Lorentzian. That crossover is the daily bread of atomic clocks. Along the way we assemble the rate bookkeeping: why a level's width sums the rates of all its exit channels, why a line carries the widths of both its levels, and how the numbers ladder from Lyman-α (100 MHz) through the Mössbauer standard (nano-eV) to millihertz strontium clock lines. The rigorous decay theory behind the exponential is owned elsewhere and cited; here you learn to read every linewidth as a stopwatch.

⚡ Key Takeaways
  • A state decaying as ∣ce(t)∣2=e−Γt|c_e(t)|^2 = e^{-\Gamma t} has lifetime τ=1/Γ\tau = 1/\Gamma; its emitted amplitude is the damped phasor e−iω0te−Γt/2e^{-i\omega_0 t}e^{-\Gamma t/2}.
  • Its power spectrum is exactly Lorentzian: FWHM Δω=Γ\Delta\omega = \Gamma, Δν=1/2πτ\Delta\nu = 1/2\pi\tau, ΔE=ℏ/τ\Delta E = \hbar/\tau.
  • A level's width sums all its exit rates, Γ=∑fAf\Gamma = \sum_f A_f; a line's width sums both its levels, Γe+Γg\Gamma_e + \Gamma_g.
  • A finite watch TT gives a Fourier-limited profile: FWHM =11.2 Γ= 11.2\,\Gamma at T=τ/2T=\tau/2, 2.87 Γ2.87\,\Gamma at 2τ2\tau, 1.01 Γ1.01\,\Gamma at 10τ10\tau — computed exactly below.
  • Doppler broadening buries the natural width (×305 for room-temperature Lyman-α) but can be engineered away; the natural floor cannot.
  • Key idea: exponential decay in time and a Lorentzian line in frequency are one object seen through the two ends of a Fourier transform — the lifetime is the linewidth.
💡 The Core Insight

A state with lifetime τ\tau radiates a wave train that rings down in time τ\tau, and its spectrum is therefore an exact Lorentzian of width Γ=1/τ\Gamma = 1/\tau: decay rate and linewidth are the same number read in two Fourier-conjugate domains.

The physical picture

Strike a bell and listen. The pitch is definite while the tone lasts, but the tone does not last. A Fourier analyzer reports a distribution of frequencies, centred on the pitch, whose width is the reciprocal of the ring-down time: a short clang is broad, a sustained hum narrow. An excited atom is a bell struck once — its dipole oscillates at ω0\omega_0 while the excited amplitude drains at rate Γ/2\Gamma/2, and the emitted light carries exactly that damped waveform.

Figure
A struck bell in light: the damped wave train and its Lorentzian spectrum are one object
time domain: e⁻ᴳᵗ/²cos ω₀t frequency domain: Lorentzian envelope e⁻ᴳᵗ/² lifetime τ = 1/Γ sets the ring-down FWHM = Γ ω₀ heavy 1/(ω−ω₀)² wings — not Gaussian ⇌ Fourier pair

The shape matters as much as the width. A single exponential decay produces a Lorentzian — algebraic 1/(ω−ω0)21/(\omega-\omega_0)^2 wings that outlive any Gaussian: far from line centre, where the Doppler core has vanished, the damping wings of strong stellar lines are Lorentzian physics. And the shape warns against the folklore "the energy is uncertain by ℏ/2τ\hbar/2\tau": the Lorentzian is so heavy-tailed that its standard deviation diverges — the honest statements are the FWHM dictionary derived below and the orthogonalization bounds owned by The Energy-Time Uncertainty Relation and Mandelstam-Tamm Bound and Evolution Timescales. There is a second, more modern picture: the linewidth as an interrogation budget. The Lorentzian is what you record after watching forever; a finite window TT gives a ridged, broad profile narrowing as 1/T1/T until, around ten lifetimes, it settles onto the floor. The figure below is computed exactly: the same emitter looks eleven natural widths wide at T=τ/2T=\tau/2 and within one percent of the floor at T=10τT=10\tau. Clock builders live on this curve — a millihertz line is useless without the seconds of coherent watching that resolve it.

Figure
The interrogation crossover: finite watching time T morphs a broad Fourier-limited ridge into the natural Lorentzian floor
peak-normalized emission spectra vs detuning (ω−ω₀)/Γ, watched for a window T 0 −4 +4 +8 −8 T = τ/2 → FWHM 11.2 Γ T = 2τ → 2.87 Γ T = 10τ → 1.01 Γ (within 1% of the Lorentzian floor) Fourier side-ridges — real, not noise every curve computed from I(ω;T) = |∫₀ᵀ e^{i(ω−ω₀)t−Γt/2} dt|² — Fourier limit ≈ 5.57/T for T ≪ τ, natural floor Γ for T ≫ τ

Finally, the ladder of magnitudes. Natural widths are fractionally tiny — parts in 10710^7 for allowed optical lines — and plummet for forbidden transitions, because Γ∝ω03∣d∣2\Gamma \propto \omega_0^3 |\mathbf d|^2 is merciless to small matrix elements and low frequencies. That is a resource: precision measurement is the art of choosing states that refuse to decay.

Figure
The quality-factor ladder: four real lines spanning ten orders of magnitude in Q = ν₀/Δν
how sharp is sharp: Q = ν₀/Δν on a log axis 10⁶ 10⁸ 10¹⁰ 10¹² 10¹⁴ 10¹⁶ 10¹⁸ H Lyman-α Q = 2.5×10⁷ · Δν = 99.8 MHz · τ = 1.6 ns Na D₂ Q = 5.2×10⁷ · Δν = 9.79 MHz · τ = 16.25 ns ⁵⁷Fe Mössbauer, 14.4 keV Q = 3.1×10¹² · Γ = 4.7 neV · τ = 141 ns ⁸⁷Sr clock line, 698 nm Q ≈ 4×10¹⁷ · Δν ≈ 1 mHz — doubly forbidden recoil-free γ emission buys 10⁵ forbidding the decay buys another 10⁵

The working objects: rates, lifetimes, and conversions

Everything hangs on one empirical law and the amplitude behind it. An undisturbed excited population decays exponentially — what a delayed-coincidence photon counter actually records — with a rate built from the spontaneous rates of every open channel. The deep reason the decay is exponential (and when it is not) belongs to Survival Amplitude and Decay of Unstable States; here we spend the working law.

∣ce(t)∣2=e−Γt,τ=1Γ,Γ=∑fAe→f(1)|c_e(t)|^2 = e^{-\Gamma t}, \qquad \tau = \frac{1}{\Gamma}, \qquad \Gamma = \sum_f A_{e\to f} \qquad (1)

where: ce(t)c_e(t) — amplitude to still find the atom excited; Γ\Gamma — total decay rate of the level; Ae→fA_{e\to f} — the Einstein spontaneous rate into final state ff, computed in Spontaneous Emission and the Einstein A Coefficient. Read the sum physically: the level does not know which line you are watching — every open exit door drains the same reservoir, so a level with a fast ultraviolet escape route is broad even in its feeble red line. The radiating amplitude carries the oscillation and half the decay exponent — half, because the probability ∣ce∣2|c_e|^2 must decay at the full rate Γ\Gamma. This factor of two is the most common slip in linewidth arithmetic:

ce(t)=ce(0) e−iω0t e−Γt/2 θ(t)(2)c_e(t) = c_e(0)\, e^{-i\omega_0 t}\, e^{-\Gamma t/2}\,\theta(t) \qquad (2)

where: ω0=(Ee−Eg)/ℏ\omega_0 = (E_e - E_g)/\hbar — the Bohr frequency of the transition being watched; θ(t)\theta(t) — the decay starts when the state is prepared. This damped phasor is the Weisskopf-Wigner ansatz: in the quantum theory of radiation it emerges as the resummed solution of the atom-field coupling, repairing the naive first-order probability that grows linearly forever. Its rigorous credentials — when a spectrum yields exact exponential decay, and where the quadratic and power-law eras hide — are Survival Amplitude and Decay of Unstable States' property; we use the ansatz as every laser lab does. Experimentalists quote linewidths in hertz — and 2π2\pi slips are endemic — so keep the dictionary at hand:

ΔωFWHM=Γ,ΔνFWHM=Γ2π=12πτ,ΔEFWHM=ℏ Γ=ℏτ(3)\Delta\omega_{\text{FWHM}} = \Gamma, \qquad \Delta\nu_{\text{FWHM}} = \frac{\Gamma}{2\pi} = \frac{1}{2\pi\tau}, \qquad \Delta E_{\text{FWHM}} = \hbar\,\Gamma = \frac{\hbar}{\tau} \qquad (3)

where: the three widths are one Lorentzian measured in angular frequency, ordinary frequency, and energy. A 1.6 ns lifetime is Γ=6.3×108 s−1\Gamma = 6.3\times10^8\ \mathrm{s^{-1}} but Δν≈100\Delta\nu \approx 100 MHz — a factor-6.28 slip that has survived referees. The energy version, ΔE τ=ℏ\Delta E\,\tau = \hbar, is the honest "energy-time uncertainty" of a decaying state: an equality about a distribution's FWHM, not a bound involving a time operator (there is none — The Energy-Time Uncertainty Relation owns that argument, and Pauli's proof of it).

Derivation AThe Lorentzian lineshape and the finite-observation crossover

We want the spectrum a spectrometer records from a decaying emitter, twice over. First the classic result: Fourier-transform the damped phasor over an infinite window and obtain an exact Lorentzian whose FWHM is the decay rate Γ\Gamma — the lifetime-linewidth theorem. Then the part this article owns: redo the same integral with a finite watching window TT and extract how the measured width relaxes from the Fourier-limited regime onto the natural floor. The strategy is one convergent integral and one modulus squared; the thing to watch is where Γ/2\Gamma/2 versus Γ\Gamma enters — the amplitude decays at half the population rate.

Step 1. Write the emitted field as the causal damped phasor E(t)=E0 e−iω0te−Γt/2E(t) = E_0\, e^{-i\omega_0 t} e^{-\Gamma t/2} for t≥0t \ge 0 and zero before, and define the detuning δ=ω−ω0\delta = \omega - \omega_0. Its Fourier transform is elementary,

E~(ω)=E0∫0∞eiδt e−Γt/2 dt=E0Γ/2−i δ,\tilde E(\omega) = E_0 \int_0^{\infty} e^{i\delta t}\, e^{-\Gamma t/2}\, \mathrm{d}t = \frac{E_0}{\Gamma/2 - i\,\delta},

convergent precisely because the state decays. Why the lower limit is zero. The atom is prepared at t=0t=0; there is no field before. Causality truncates the integral — and that truncation is what forbids a delta-function spectrum.

Step 2. The spectrometer records the modulus squared:

I(ω)=∣E~(ω)∣2=∣E0∣2(ω−ω0)2+(Γ/2)2.I(\omega) = |\tilde E(\omega)|^2 = \frac{|E_0|^2}{(\omega-\omega_0)^2 + (\Gamma/2)^2}.

Phase information cancels; a Lorentzian centred on the Bohr frequency survives.

Step 3. Half-maximum points: II falls to half its peak at ω=ω0±Γ/2\omega = \omega_0 \pm \Gamma/2, so

ΔωFWHM=Γ=1τ.\Delta\omega_{\text{FWHM}} = \Gamma = \frac{1}{\tau}.

Why Lorentzian and not Gaussian. A Gaussian needs a spread of centre frequencies (Doppler); one emitter with one decay constant transforms to the Lorentzian with its algebraic wings.

Step 4. Now watch for only a window TT — cut the integral off. The same arithmetic gives the closed form

I(ω;T)=∣∫0Teiδt−Γt/2 dt∣2=1+e−ΓT−2 e−ΓT/2cos⁡(δ T)(ω−ω0)2+(Γ/2)2.I(\omega; T) = \left|\int_0^{T} e^{i\delta t - \Gamma t/2}\,\mathrm{d}t\right|^2 = \frac{1 + e^{-\Gamma T} - 2\,e^{-\Gamma T/2}\cos(\delta\,T)}{(\omega-\omega_0)^2 + (\Gamma/2)^2}.

The Lorentzian denominator survives, but the numerator ripples in δ\delta with period 2π/T2\pi/T — the side-ridges of the sharp cutoff, the same cos⁡(δT)\cos(\delta T) arithmetic that gives first-order perturbation theory its sinc2\mathrm{sinc}^2 profile in First-Order Transition Amplitudes and Probabilities.

Step 5. Extract the measured FWHM numerically. Short window: decay irrelevant, Fourier-limited FWHM →5.566/T\to 5.566/T (the sinc2\mathrm{sinc}^2 constant). Long window: the ripples die under e−ΓT/2e^{-\Gamma T/2} and the Lorentzian survives. The crossover, computed: FWHM =11.2 Γ= 11.2\,\Gamma at T=τ/2T=\tau/2, 2.87 Γ2.87\,\Gamma at 2τ2\tau, 1.01 Γ1.01\,\Gamma at 10τ10\tau — ten lifetimes of coherent watching recover the floor to one percent. Verified: the infinite-window transform, the half-maximum solution, the finite-window closed form, the 5.566/T Fourier constant, and all three crossover widths are re-derived symbolically and numerically in verify_aq-qm-lifetime-linewidth-broadening.py (34 checks green).

The energy reading: what the line "knows" about the state

A spectral line is also an energy histogram: expand the decaying state over the exact stationary states of atom-plus-field, and the survival amplitude is the Fourier transform of the energy distribution — the Krylov-Fock master formula that Survival Amplitude and Decay of Unstable States owns — so the observed decay determines the distribution:

P(E)=ℏΓ/2π(E−Ee)2+(ℏΓ/2)2,ΔEFWHM τ=ℏ(4)P(E) = \frac{\hbar\Gamma/2\pi}{(E - E_e)^2 + (\hbar\Gamma/2)^2}, \qquad \Delta E_{\text{FWHM}}\,\tau = \hbar \qquad (4)

where: P(E)P(E) — normalized energy distribution of the decaying state; EeE_e — the nominal level position; ℏΓ\hbar\Gamma — its FWHM. An unstable "level" is a Lorentzian band of true stationary states, and the emitted line photographs the band.

Derivation BThe energy distribution of a decaying state

We want to show that exponential decay forces the Lorentzian energy distribution and the exact product ΔEFWHM⋅τ=ℏ\Delta E_{\text{FWHM}}\cdot\tau = \hbar — and to expose what this does and does not say about "energy-time uncertainty." The route inverts the lineshape derivation: start from the survival amplitude as the Fourier transform of an unknown energy distribution, and invert. The subtlety worth watching: the resulting distribution's wings are so heavy that its variance diverges, so standard-deviation statements must be handled with care, not folklore.

Step 1. The survival amplitude is the Fourier transform of the energy distribution, a(t)=∫P(E) e−iEt/ℏ dEa(t) = \int P(E)\, e^{-iEt/\hbar}\,\mathrm dE — the spectral theorem the decay article owns. Imposing ∣a(t)∣2=e−Γt|a(t)|^2 = e^{-\Gamma t} with a pure phase fixes a(t)=e−iEet/ℏe−Γ∣t∣/2a(t) = e^{-iE_e t/\hbar} e^{-\Gamma|t|/2}.

Step 2. Invert the transform. The Fourier transform of a two-sided decaying exponential is again elementary:

P(E)=12πℏ∫−∞∞a(t) eiEt/ℏ dt=ℏΓ/2π(E−Ee)2+(ℏΓ/2)2.P(E) = \frac{1}{2\pi\hbar}\int_{-\infty}^{\infty} a(t)\, e^{iEt/\hbar}\,\mathrm dt = \frac{\hbar\Gamma/2\pi}{(E-E_e)^2 + (\hbar\Gamma/2)^2}.

Why this is the same Lorentzian as the lineshape. The emitted photon inherits the energy bookkeeping of the state that made it: line profile and level energy-distribution are Fourier images of the same decay.

Step 3. Read off the exact product. FWHM =ℏΓ= \hbar\Gamma and τ=1/Γ\tau = 1/\Gamma, so ΔEFWHM τ=ℏ\Delta E_{\text{FWHM}}\,\tau = \hbar — an equality, with the full width, not ℏ/2\hbar/2 with a standard deviation.

Step 4. Try to compute the standard deviation and fail informatively: ∫(E−Ee)2P(E) dE\int (E-E_e)^2 P(E)\,\mathrm dE diverges (checked symbolically). The Lorentzian has no σE\sigma_E, so the Robertson-style bounds of Mandelstam-Tamm Bound and Evolution Timescales — which itself flags the Lorentzian's divergent second moment — apply vacuously here. The FWHM dictionary of Step 3 is the statement with content. Verified: the inversion integral, the unit normalization, the ΔE·τ = ħ product, and the divergence of the second moment are all checked symbolically in the verification script.

Every open door widens the line

Two pieces of bookkeeping make the theorem laboratory-grade. Levels: independent channel rates add — equation (1); the branching ratio says how often each door is taken, but every observed line reflects the sum. Lines: if the lower level is itself unstable, the dipole coherence decays at the sum of the two half-rates.

Δωline=Γe+Γg(5)\Delta\omega_{\text{line}} = \Gamma_e + \Gamma_g \qquad (5)

where: Γe\Gamma_e, Γg\Gamma_g — total decay rates of upper and lower level. A transition between two broad levels is broader than either. This is why the strontium clock line pairs a ground state with a doubly forbidden excited state — both widths near zero — and why Forbidden Transitions: Magnetic Dipole and Electric Quadrupole is this article's natural sequel.

Derivation CWhy the line width adds the widths of both levels

The claim: a transition between an upper level decaying at Γe\Gamma_e and a lower level decaying at Γg\Gamma_g radiates a Lorentzian of full width Γe+Γg\Gamma_e + \Gamma_g — not Γe2+Γg2\sqrt{\Gamma_e^2+\Gamma_g^2}, not the larger of the two. The strategy is to identify the object whose transform is the line: not a population but the dipole coherence between the two amplitudes, each carrying its own half-rate damping. The pitfall is squaring too early — the rates add in the amplitude product first.

Step 1. The radiated field is proportional to the dipole coherence d(t)∝ce(t) cg∗(t)d(t) \propto c_e(t)\,c_g^{*}(t); each amplitude damps at half its level's rate, ce∼e−iωete−Γet/2c_e \sim e^{-i\omega_e t}e^{-\Gamma_e t/2}, cg∼e−iωgte−Γgt/2c_g \sim e^{-i\omega_g t}e^{-\Gamma_g t/2}.

Step 2. Multiply: d(t)∝e−iω0t e−(Γe+Γg)t/2d(t) \propto e^{-i\omega_0 t}\, e^{-(\Gamma_e + \Gamma_g)t/2}. The half-rates add in the exponent — coherence dies if either partner dies.

Step 3. This is exactly the damped phasor of the lineshape derivation with Γ→Γe+Γg\Gamma \to \Gamma_e + \Gamma_g, so the spectrum is Lorentzian with

ΔωFWHM=Γe+Γg.\Delta\omega_{\text{FWHM}} = \Gamma_e + \Gamma_g.

Why a ground state contributes zero. No decay channel means Γg=0\Gamma_g = 0: the line width collapses to the upper-level width — the default case in tables. Verified: the Fourier transform of the product coherence and the additive FWHM are checked symbolically in the verification script; the hydrogen 3p worked example below runs the arithmetic on real rates.

Limiting cases

caseresultmeaning
isolated atom, watched foreverLorentzian, FWHM Γ=1/τ\Gamma = 1/\tauthe natural floor; no experiment gets below it
finite window T≪τT \ll \tauFWHM ≈5.57/T\approx 5.57/T, sinc-ridgedFourier-limited: the apparatus sets the width
warm vaporGaussian, width ω08ln⁡2 kBT/mc2\omega_0\sqrt{8\ln 2\, k_BT/mc^2}Doppler; ≈ 305 × natural for H Lyman-α at 300 K
natural + DopplerVoigt profileGaussian core, Lorentzian wings — wings are always natural/collisional
buffer gasLorentzian, width ≈2/τcoll\approx 2/\tau_{\text{coll}}phase-interrupting collisions, same Fourier logic
lower level also decaysFWHM =Γe+Γg= \Gamma_e + \Gamma_gcoherence dies with either partner
Γ→0\Gamma \to 0 (metastability)line sharpens toward a delta functionthe stationary-state limit — the clock builder's grail
⚠️ A 1.6 ns lifetime means the emitted photon's energy is smeared — the atom's energy levels are not blurry

The Lorentzian width belongs to the decaying state as prepared — the joint atom-field dynamics — not to any fuzziness of the Hamiltonian's spectrum. The true stationary states of atom + field are perfectly sharp; the short-lived "level" is a wave packet over them, and P(E)P(E) in equation (4) is its composition. Two exam-favourite corollaries: (i) no time operator and no commutator stands behind ΔE τ=ℏ\Delta E\,\tau = \hbar — it is Fourier arithmetic on a distribution whose variance does not even exist; (ii) different lines from the same upper level carry the same total width Γ=∑fAf\Gamma = \sum_f A_f.

Worked examples

Example 1 — Lyman-α end to end: from a rate to a wavelength smear. The hydrogen 2p2p level decays only to 1s1s, at A(2p→1s)=6.27×108 s−1A(2p{\to}1s) = 6.27\times10^{8}\ \mathrm{s^{-1}} — the flagship number of Spontaneous Emission and the Einstein A Coefficient, re-derived from the radial integral in this article's verification script. Find the lifetime, linewidth, wavelength width, and the Doppler comparison.

Step 1. Lifetime: τ=1/A=1.60 ns\tau = 1/A = \textbf{1.60 ns} — set entirely by ω03 ∣⟨1s∣r∣2p⟩∣2\omega_0^3\,|\langle 1s|\mathbf r|2p\rangle|^2.

Step 2. Natural width: Δν=A/2π=99.8 MHz\Delta\nu = A/2\pi = \textbf{99.8 MHz}. Against ν0=c/λ=2.47×1015\nu_0 = c/\lambda = 2.47\times10^{15} Hz that is Q=ν0/Δν=2.5×107Q = \nu_0/\Delta\nu = \textbf{2.5}\times\textbf{10}^{7}.

Step 3. Wavelength width at λ=121.567\lambda = 121.567 nm: Δλ=λ2Δν/c=4.9×10−5 A˚\Delta\lambda = \lambda^2\Delta\nu/c = \textbf{4.9}\times\textbf{10}^{-5}\ \textbf{Å}. No grating spectrograph resolves that.

Step 4. Doppler at 300 K: ΔνD=ν08ln⁡2 kBT/mHc2=3.0×1010 Hz\Delta\nu_D = \nu_0\sqrt{8\ln 2\,k_BT/m_Hc^2} = \textbf{3.0}\times\textbf{10}^{10}\ \textbf{Hz} — a factor ≈ 305 above natural. The natural line hides under the thermal Gaussian until you cool, trap, or Doppler-cancel.

Example 2 — the hydrogen 3p ledger: one level, two lines, one width. The 3p3p level has two E1 exit doors: 3p→1s3p \to 1s (Lyman-β) at A1=1.672×108 s−1A_1 = 1.672\times10^{8}\ \mathrm{s^{-1}} and 3p→2s3p \to 2s (an Hα component) at A2=2.245×107 s−1A_2 = 2.245\times10^{7}\ \mathrm{s^{-1}} — both computed from first-principles radial integrals in the verification script (matching NIST to 0.1%). Find the level width, branching, and both line widths, given that 2s2s decays only by two-photon emission at 8.23 s−18.23\ \mathrm{s^{-1}}.

Step 1. Total rate: Γ3p=A1+A2=1.897×108 s−1\Gamma_{3p} = A_1 + A_2 = 1.897\times10^{8}\ \mathrm{s^{-1}}, so τ3p=5.27 ns\tau_{3p} = \textbf{5.27 ns}.

Step 2. Branching: 88.2%\textbf{88.2\%} of decays take the Lyman-β door, 11.8%\textbf{11.8\%} the Hα door.

Step 3. Lyman-β width: 1s1s is stable, so Δν=Γ3p/2π=30.2 MHz\Delta\nu = \Gamma_{3p}/2\pi = \textbf{30.2 MHz}.

Step 4. Hα-component width: 2s2s adds 8.23/2π=1.3 Hz8.23/2\pi = \textbf{1.3 Hz} — seven orders down — so this line is also 30.2 MHz wide. The minority line is exactly as broad as the majority line: the width belongs to the level, not the branch. And the freakish narrowness of the 2s2s door (0.12 s lifetime) is forbidden-transition physics — why 2s2s cannot take a one-photon exit at all is owned by Forbidden Transitions: Magnetic Dipole and Electric Quadrupole.

NL
Natural Linewidth
1930 · Göttingen and Berlin

The natural linewidth taught quantum theory to talk about decay honestly. Lorentz's classical damped oscillator already gave a Lorentzian, but only in 1930 did Victor Weisskopf, then a doctoral student, and Eugene Wigner compute the shape from Dirac's new radiation theory — feeding the decay back into the amplitude so a linearly growing perturbative probability became a clean exponential with a Lorentzian spectrum of width Γ=1/τ\Gamma = 1/\tau. The calculation quietly standardized a complex self-energy, a level shift alongside the width (the seed of the Lamb-shift industry seventeen years later), and the reading of an unstable state as a resonance band of true stationary states. Every linewidth in every atomic-data table since is a Weisskopf-Wigner number.

⏳ History
1916
— Einstein introduces the spontaneous-emission coefficient A, tying level lifetimes to radiation
1927
— Dirac quantizes the radiation field, deriving spontaneous emission from first principles
1930
— Weisskopf and Wigner derive the natural lineshape: exponential decay, Lorentzian profile, FWHM Γ = 1/τ
1945
— Mandelstam and Tamm prove the rigorous energy-time bound for state evolution in Moscow
1958
— Mössbauer discovers recoil-free gamma emission, exposing nuclear natural widths of nano-eV sharpness
2015
— optical lattice clocks interrogating millihertz-wide lines reach fractional uncertainties of a few parts in 10¹⁸

Where this lives today

The natural linewidth is the currency of precision measurement. Optical atomic clocks run on transitions with pathologically small Γ\Gamma: the 87^{87}Sr clock line at 698 nm has a natural width near one millihertz (Q≈4×1017Q \approx 4\times10^{17}), and lattice clocks resolving fractions of it reach 10−1810^{-18} fractional uncertainty. The interrogation crossover derived above is their daily engineering reality: seconds of coherent probing are needed before the measured line approaches the floor — which is why clock lasers must themselves be stabilized to hertz-level linewidths. The same physics runs both directions in the quantum-optics lab: the natural width sets the bandwidth and temporal shape of single photons from atoms and quantum dots, and a cavity changes Γ\Gamma (the Purcell effect) — the linewidth belongs to the atom plus its electromagnetic environment. In astrophysics the Lorentzian damping wings of Lyman-α shape galaxy absorption spectra; in nuclear metrology the 4.7 neV Mössbauer width of ⁵⁷Fe registers a building's worth of gravitational redshift.

Check yourself

Going further

References. Griffiths & Schroeter, Introduction to Quantum Mechanics, 3rd ed., Ch. 11 · Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., Ch. 5 · Loudon, The Quantum Theory of Light, 3rd ed., Ch. 2 · Cohen-Tannoudji, Diu & Laloë, Vol. II, Ch. XIII · Weisskopf & Wigner, Z. Phys. 63, 54 (1930) · Foot, Atomic Physics, Ch. 7-8.

Provenance. Lineshape and FWHM derivations follow Loudon (Quantum Theory of Light, Ch. 2) and Cohen-Tannoudji Vol. II; decay-rate formula Γ=ω03∣deg∣2/3πε0ℏc3\Gamma=\omega_0^3|d_{eg}|^2/3\pi\varepsilon_0\hbar c^3 is the standard Einstein-AA/Weisskopf–Wigner result (verify against Sakurai Ch. 5, Griffiths Ch. 11). A verifier should recompute: (1) the Fourier transform giving the Lorentzian and FWHM =Γ=\Gamma; (2) the energy distribution P(E)P(E) and ΔE τ=ℏ\Delta E\,\tau=\hbar; (3) the numerical estimate Γ∼108 s−1\Gamma\sim10^8\,\mathrm{s^{-1}}, Δν∼10 MHz\Delta\nu\sim10\,\mathrm{MHz} for λ≈600 nm\lambda\approx600\,\mathrm{nm}, ∣deg∣≈ea0|d_{eg}|\approx ea_0, and compare with the standard sodium-D lifetime ≈16 ns\approx16\,\mathrm{ns}.

Builds on: Fermi's Golden Rule · First-Order Transition Amplitudes and Probabilities · Absorption and Stimulated Emission · The Dipole Approximation. Leads to: Semiclassical Theory of Atom–Radiation Interaction · Harmonic (Sinusoidal) Perturbations and Resonance.

Does the natural width ever dominate in a real experiment, given Doppler is hundreds of times larger?

Yes, routinely — that is what "Doppler-free" methods are for. Saturated-absorption spectroscopy, two-photon spectroscopy with counter-propagating beams, atomic beams observed at right angles, and laser-cooled samples all cancel or remove the thermal shifts, and the recorded profile collapses to the natural Lorentzian (plus residual collisional and transit-time effects). Every hyperfine interval quoted to kHz was measured under the Doppler floor.

Why does the amplitude decay at Γ/2 but the probability at Γ?

Probability is the amplitude squared: |e−Γt/2|² = e−Γt. But it is the amplitude that interferes and radiates — which is why the coherence between two decaying levels damps at (Γ_e+Γ_g)/2 and the line width comes out Γ_e+Γ_g.

Is exponential decay exact? The derivation seems to assume it.

An approximation — extraordinarily good in the middle. The exact survival probability starts quadratically (the Zeno era), runs exponentially for many lifetimes, then crosses to a power law when the spectrum's lower edge asserts itself (Khalfin's theorem). Survival Amplitude and Decay of Unstable States owns that three-era anatomy and the statement that a Lorentzian band gives exact exponential decay while a physical, bounded-below spectrum is Lorentzian only in the middle. For linewidths the middle is all that matters.

What exactly does a "1 mHz natural linewidth" mean for a clock transition nobody can watch for 1000 seconds?

It means the floor is 1 mHz; the measured line is Fourier-limited by the interrogation time, exactly as in the crossover figure. A 1 s probe gives a roughly hertz-wide feature; longer coherent probes narrow it toward the floor. Clock development is therefore as much about ultra-stable lasers (keeping the probe coherent) as about the atom.

Do collisions and Doppler broadening also come from a Fourier transform?

Yes — same theorem, different truncation. Phase-interrupting collisions chop the wave train into random segments of mean length τ_coll, giving a Lorentzian of width ≈ 2/τ_coll. Doppler broadening is different in kind: each atom emits its narrow line at a thermally shifted centre, giving a Gaussian inhomogeneous profile. Convolving the two gives the Voigt profile: Gaussian core, Lorentzian wings.

Why does the same Γ show up in the decay curve, the linewidth, and the Breit-Wigner resonance shape?

All three are Fourier views of one pole. A quasi-bound state at complex energy E_r − iΓ/2 decays as e−Γt, emits a Lorentzian of FWHM Γ, and scatters with the Breit-Wigner profile T(E) ≈ (Γ/2)²/((E−E_r)²+(Γ/2)²) — Resonant Transmission and Tunnelling Resonances derives that form and reads its width as the trapped state's lifetime τ = ħ/Γ, the dictionary Conservation of Normalization plants as early as the optical-potential model.

Before this

Where this leads