Steady-State Excited Population and Fluorescence

Shine a resonant laser on a single atom and it never simply "turns on and stays on." The field drives coherent Rabi flopping between ground and excited states, while spontaneous emission — the one-way channel this chapter's aq-optics-spontaneous-emission and aq-optics-einstein-a-b-coefficients articles established — randomly kicks the atom back down at the rate $\Gamma$ . These two processes compete, and after a few excited-state lifetimes the atom forgets how it started and settles into a stationary balance: excitation up exactly matches relaxation down. That balance is the steady state, and it is the single most reused calculation in atomic, molecular and optical physics. The right language is the density matrix $\rho$ of a two-level atom evolving under the optical Bloch equations — the coherent Hamiltonian drive (Rabi frequency $\Omega$ , detuning $\delta$ ) plus dissipative decay of population and coherence at rates fixed by the natural linewidth $\Gamma$ (recalled from aq-optics-radiative-lifetime-natural-linewidth). Set every time derivative to zero and the algebra closes into a compact expression for the excited-state population $\rho_{ee}$ and hence the photon-scattering rate. The punchline is saturation. However hard you drive, $\rho_{ee}$ climbs to $1/2$ and stops — a closed, coherently driven two-level atom cannot invert, because the same field that pumps it up also stimulates it down. Only spontaneous emission breaks the symmetry, capping the fluorescence at $\Gamma/2$ photons per lifetime, and scan the detuning and that fluorescence traces a Lorentzian broadening as $\Gamma\sqrt{1+s_0}$ . This one result underwrites laser cooling, atomic clocks, single-ion readout, Doppler-free spectroscopy and saturable absorbers.

⚡ Key Takeaways

A two-level atom driven by a monochromatic field and damped by spontaneous emission relaxes to a steady state whose excited population $\rho_{ee}$ never exceeds $1/2$ — no inversion, no gain, ever.
Solving the optical Bloch equations at $\dot\rho=0$ gives $\rho_{ee}=\tfrac12\,s/(1+s)$ with saturation parameter $s=s_0/[1+(2\delta/\Gamma)^2]$ and $s_0=2\Omega^2/\Gamma^2=I/I_{sat}$ .
The photon-scattering (fluorescence) rate is $R_{sc}=\Gamma\rho_{ee}$ , a power-broadened Lorentzian that saturates at $\Gamma/2$ for strong driving.
The saturation intensity $I_{sat}=\pi hc\Gamma/3\lambda^3$ fixes the scale: at $I=I_{sat}$ (so $s_0=1$ , on resonance) the excited population is exactly $1/4$ . For sodium $I_{sat}\approx6.3$ mW cm⁻².
The fluorescence line broadens as $\Delta\omega_{FWHM}=\Gamma\sqrt{1+s_0}$ while its peak height stays pinned — the width grows, the ceiling does not.
Key idea: in steady state a closed, coherently driven two-level atom obeys $\rho_{ee}=\tfrac12\,s/(1+s)$ with $s=s_0/[1+(2\delta/\Gamma)^2]$ ; the population saturates at $1/2$ and the fluorescence rate at $\Gamma/2$ .

💡 The Core Insight

In steady state a closed, coherently driven two-level atom obeys $\rho_{ee}=\tfrac12\,s/(1+s)$ with $s=s_0/[1+(2\delta/\Gamma)^2]$ ; the excited population saturates at $1/2$ and the fluorescence rate at $\Gamma/2$ , so a two-level atom can never invert or amplify.

The physical picture

Picture the atom as a switch that the laser keeps trying to flip up while spontaneous emission keeps knocking it back down. The coherent field would, on its own, drive full Rabi oscillations — ground to excited and back, endlessly. Spontaneous emission interrupts that clean flopping at random moments, resetting the phase and dumping the atom into the ground state with a scattered photon. When the up-rate and the down-rate come into balance the populations stop changing on average, and the atom sits in a fixed mixture of ground and excited. Why does the excited population climb to $1/2$ and stop, rather than reaching $1$ ? Because the drive is symmetric. The same field amplitude that stimulates absorption ( $g\to e$ ) stimulates emission ( $e\to g$ ) at the identical rate — the Einstein $B$ coefficients are equal — so stimulated processes alone can only equalise the two populations, never pile more atoms upstairs than down. The one asymmetric channel is spontaneous emission, which always points downward. So the ceiling is a tie, $\rho_{ee}=\rho_{gg}=1/2$ , and the leftover spontaneous decay is exactly the steady stream of fluorescence photons you detect.

Figure

Saturation of the excited-state population on resonance ρ_ee = ½ s₀/(1+s₀). The curve rises linearly at weak drive, passes through exactly ¼ at the saturation intensity s₀ = 1, and asymptotes to the hard ceiling ½ — a driven two-level atom never inverts

Detuning weakens the effective drive. Off resonance the atom is out of step with the field, its coherence is small, and the scattering rate collapses as $1/\delta^2$ — the Lorentzian wings. Scan a probe laser across the transition and you map out that Lorentzian directly as fluorescence versus frequency. Turn the intensity up and something instructive happens: the peak height stops rising once $\rho_{ee}$ hits its ceiling, but the width keeps growing, because atoms that were formerly too far detuned to respond are now driven hard enough to scatter. That is power broadening, and it is the second figure.

Figure

Power broadening of the fluorescence line R_sc/Γ versus detuning δ/Γ. As the drive s₀ grows the Lorentzian widens — FWHM = Γ√(1+s₀), so s₀ = 8 triples the width to 3Γ — while the peak creeps toward but never passes the hard ceiling R_sc = Γ/2

At weak drive the whole picture collapses back to Einstein rate equations: absorption balances spontaneous decay, $\rho_{ee}$ grows linearly with intensity, and there is no saturation at all. Saturation is precisely the breakdown of that linear story once stimulated emission becomes a competitive return path — and the density-matrix treatment contains the rate-equation answer as its low-intensity limit while going beyond it to predict the ceiling, the power broadening, and (at very strong drive) the Mollow triplet.

The core equations

The steady state is what you get by demanding that nothing changes in time. Solving the optical Bloch equations at $\dot\rho=0$ (Derivation 1) delivers the excited-state population directly in terms of the drive $\Omega$ , the detuning $\delta$ , and the decay rate $\Gamma$ . This is the master result from which everything else is a rearrangement or a multiplication.

\rho_{ee}^{ss}=\frac{\Omega^2/4}{\delta^2+\Omega^2/2+\Gamma^2/4} \qquad (1)

where: $\rho_{ee}^{ss}$ — steady-state excited-state population (a probability, $0\le\rho_{ee}\le\tfrac12$ ); $\Omega$ — resonant Rabi frequency, $\Omega=-dE_0/\hbar$ ; $\delta=\omega_L-\omega_0$ — laser detuning from resonance; $\Gamma=1/\tau$ — natural linewidth / spontaneous decay rate. Physical reading: the numerator is the coherent drive, the denominator is drive plus detuning plus decay — push $\Omega\to\infty$ and the ratio tends to $1/2$ , never more. The same expression becomes far more transparent in terms of the saturation parameter, which folds the Rabi frequency and the transition's radiative strength into a single dimensionless "how hard am I driving relative to the natural scale." Dividing through by $\Gamma^2/4$ (Derivation 1's coda) produces the canonical form.

\rho_{ee}^{ss}=\frac{1}{2}\,\frac{s}{1+s},\qquad s=\frac{s_0}{1+(2\delta/\Gamma)^2},\qquad s_0=\frac{2\Omega^2}{\Gamma^2}=\frac{I}{I_{sat}} \qquad (2)

where: $s$ — detuning-dependent saturation parameter; $s_0$ — on-resonance saturation parameter, equal to the intensity in units of $I_{sat}$ ; $I$ — laser intensity; $I_{sat}$ — saturation intensity. Physical reading: $s$ measures the drive against the (detuning-widened) linewidth; $\rho_{ee}=\tfrac12\,s/(1+s)$ is a smooth rise from $0$ to the ceiling $\tfrac12$ , hitting $\tfrac14$ exactly when $s=1$ . Multiply the population by the decay rate $\Gamma$ and you get the observable: the rate at which the atom scatters photons out of the beam. This is what a photodetector counts, and it is a power-broadened Lorentzian in the detuning.

R_{sc}=\Gamma\,\rho_{ee}^{ss}=\frac{\Gamma}{2}\,\frac{s_0}{1+s_0+(2\delta/\Gamma)^2} \qquad (3)

where: $R_{sc}$ — steady-state photon-scattering (fluorescence) rate per atom. Physical reading: on resonance and at full saturation $R_{sc}\to\Gamma/2$ , the hard maximum; every scattered photon is one spontaneous decay, so the atom fluoresces at most once per two lifetimes. The saturation intensity is not a free parameter — it is fixed by the transition itself. Averaging the dipole coupling over polarisation and orientation for an ideal two-level atom gives the standard isotropic result, which sets the absolute intensity scale for the whole problem.

I_{sat}=\frac{\pi h c\,\Gamma}{3\lambda^3}=\frac{\pi h c}{3\lambda^3\tau} \qquad (4)

where: $\lambda$ — transition wavelength; $h$ — Planck constant; $c$ — speed of light; $\tau=1/\Gamma$ — radiative lifetime. Physical reading: strong (short-wavelength) or broad (large- $\Gamma$ ) transitions saturate only at high intensity; a narrow line saturates with almost no light. For the sodium D₂ line this evaluates to about $6.3$ mW cm⁻². Scanning the detuning at fixed drive traces out a Lorentzian whose width is larger than the natural linewidth — the extra width is entirely due to the drive itself. Derivation 2 extracts it from equation (3).

\Delta\omega_{FWHM}=\Gamma\sqrt{1+s_0} \qquad (5)

where: $\Delta\omega_{FWHM}$ — full width at half maximum of the scattering Lorentzian. Physical reading: the line broadens as $\sqrt{I}$ ; at $s_0=1$ it is $\Gamma\sqrt2$ , at $s_0=3$ it doubles to $2\Gamma$ , at $s_0=8$ it triples to $3\Gamma$ — while the peak height stays pinned at its saturated value. All of this descends from the optical Bloch equations: the coupled equations of motion for the excited population and the optical coherence in the frame rotating at the laser frequency, with the rotating-wave approximation applied.

\dot\rho_{ee}=-\Gamma\rho_{ee}+\frac{i\Omega}{2}(\tilde\rho_{eg}-\tilde\rho_{ge}),\qquad \dot{\tilde\rho}_{ge}=-\left(\frac{\Gamma}{2}+i\delta\right)\tilde\rho_{ge}+\frac{i\Omega}{2}(\rho_{ee}-\rho_{gg}) \qquad (6)

where: $\rho_{gg}$ — ground-state population, with $\rho_{gg}+\rho_{ee}=1$ ; $\tilde\rho_{ge}$ — slowly varying optical coherence in the rotating frame. Physical reading: population decays at $\Gamma$ , coherence at $\Gamma/2$ (the factor of two is why line shapes come out Lorentzian), and the drive $\Omega$ shuttles amplitude between the two. Finally, because each scattered photon carries momentum $\hbar k$ , the saturated fluorescence rate translates directly into a mechanical force — the engine of laser cooling.

R_{sc}^{max}=\frac{\Gamma}{2},\qquad F_{max}=\hbar k\,\frac{\Gamma}{2} \qquad (7)

where: $F_{max}$ — maximum spontaneous scattering (radiation-pressure) force; $k=2\pi/\lambda$ — laser wavenumber; $\hbar k$ — photon recoil momentum. Physical reading: the force saturates because the scattering rate does — you cannot push an atom harder than one recoil kick per two lifetimes, no matter how bright the laser.

Derivation ASteady state of the optical Bloch equations

The goal is to take the coupled optical Bloch equations for the excited population and the optical coherence, impose the steady-state condition that every time derivative vanishes, and solve the resulting algebra to obtain the closed-form excited-state population $\rho_{ee}^{ss}=(\Omega^2/4)/(\delta^2+\Omega^2/2+\Gamma^2/4)$ . We work in the frame rotating at the laser frequency under the rotating-wave approximation, use population conservation $\rho_{gg}=1-\rho_{ee}$ , and eliminate the coherence to leave a single equation for the population — which is where the ceiling of $1/2$ becomes visible.

Step 1. Set the coherence derivative to zero: $0=-(\tfrac{\Gamma}{2}+i\delta)\tilde\rho_{ge}+\tfrac{i\Omega}{2}(\rho_{ee}-\rho_{gg})$ , and abbreviate the population difference as $w=\rho_{ee}-\rho_{gg}$ . Why: in steady state the coherence is enslaved to the populations — it stops evolving and simply tracks whatever the population difference is, so we can solve for it algebraically.

Step 2. Rearrange to $\tilde\rho_{ge}=\dfrac{i\Omega/2}{\Gamma/2+i\delta}\,w$ . Why: this exposes that the coherence is proportional both to the drive $\Omega$ and to the population difference $w$ ; a fully saturated atom with $w=0$ carries no coherence, and hence radiates only incoherently.

Step 3. Insert into the steady population equation $0=-\Gamma\rho_{ee}+\tfrac{i\Omega}{2}(\tilde\rho_{eg}-\tilde\rho_{ge})$ , using $\tilde\rho_{eg}=\tilde\rho_{ge}^{*}$ so that $i(\tilde\rho_{eg}-\tilde\rho_{ge})=2\,\mathrm{Im}\,\tilde\rho_{ge}$ . Why: only the imaginary (out-of-phase) part of the coherence feeds population between the levels — the in-phase part is the light shift, and it is absorption that drives the population.

Step 4. Take the imaginary part: $\mathrm{Im}\,\tilde\rho_{ge}=\dfrac{\Omega\Gamma/2}{2(\Gamma^2/4+\delta^2)}\,w$ , giving $0=-\Gamma\rho_{ee}+\dfrac{\Omega^2\Gamma/4}{\Gamma^2/4+\delta^2}\,w$ with $w=2\rho_{ee}-1$ . Why: substituting $w=2\rho_{ee}-1$ collapses the two unknowns to a single linear equation in $\rho_{ee}$ .

Step 5. Collect the $\rho_{ee}$ terms: $\rho_{ee}\!\left[\Gamma+\dfrac{\Omega^2\Gamma/2}{\Gamma^2/4+\delta^2}\right]=\dfrac{\Omega^2\Gamma/4}{\Gamma^2/4+\delta^2}$ , and divide through to reach $\rho_{ee}^{ss}=\dfrac{\Omega^2/4}{\delta^2+\Omega^2/2+\Gamma^2/4}$ . Why: the $\Omega^2/2$ sitting in the denominator is the saturation term — it grows without bound as the drive increases, which is exactly what forces $\rho_{ee}\to1/2$ rather than $1$ .

Step 6. Divide numerator and denominator by $\Gamma^2/4$ and identify $s_0=2\Omega^2/\Gamma^2$ and $s=s_0/[1+(2\delta/\Gamma)^2]$ to land on the compact form $\rho_{ee}^{ss}=\tfrac12\,s/(1+s)$ . Why: this repackaging is what makes the ceiling and the $I_{sat}$ scale manifest, and it is the form used everywhere in laser cooling. *Verified: on resonance ( $\delta=0$ ), $\rho_{ee}=\tfrac12\,\Omega^2/(\Omega^2+\Gamma^2/2)\to1/2$ for $\Omega\gg\Gamma$ and $\to\Omega^2/\Gamma^2$ for $\Omega\ll\Gamma$ , matching the Einstein rate-equation limit. At $s_0=1$ , $\rho_{ee}=1/4$ . See Foot Ch. 7 and Loudon Ch. 2.*

Derivation BPower broadening of the fluorescence Lorentzian

The goal is to show that the scattering rate as a function of detuning is a Lorentzian whose full width at half maximum grows with the drive as $\Gamma\sqrt{1+s_0}$ , not the bare natural width $\Gamma$ . We start from the observable rate $R_{sc}=\Gamma\rho_{ee}^{ss}$ , rewrite the denominator so the detuning dependence stands alone, read off the half-width where the profile drops to half its peak, and double it. The physical content is that the drive itself broadens the line — a purely intensity-dependent effect distinct from Doppler or collisional broadening.

Step 1. Multiply the boxed population by $\Gamma$ : $R_{sc}=\dfrac{\Gamma}{2}\,\dfrac{\Omega^2/2}{\delta^2+\Omega^2/2+\Gamma^2/4}$ . Why: the scattering rate is just the population drained at the spontaneous rate $\Gamma$ , so every feature of $\rho_{ee}(\delta)$ is inherited directly by the measured fluorescence.

Step 2. Use $\Omega^2/2=s_0\Gamma^2/4$ and factor $(\Gamma/2)^2$ out of the denominator: $R_{sc}=\dfrac{\Gamma}{2}\,\dfrac{s_0(\Gamma/2)^2}{\delta^2+(\Gamma/2)^2(1+s_0)}$ . Why: writing the denominator as $\delta^2+(\text{const})$ makes it a textbook Lorentzian in $\delta$ , whose half-width is the square root of that constant.

Step 3. Read off the half-width at half maximum: the profile falls to half its peak when $\delta^2=(\Gamma/2)^2(1+s_0)$ , i.e. $\delta_{HWHM}=\tfrac{\Gamma}{2}\sqrt{1+s_0}$ . Why: for a Lorentzian $1/(\delta^2+a^2)$ the half-max point is $\delta=a$ , so the driven constant is the entire source of the extra width.

Step 4. Double it: $\Delta\omega_{FWHM}=2\delta_{HWHM}=\Gamma\sqrt{1+s_0}$ . Why: this is the power-broadened linewidth — at $s_0=3$ it doubles to $2\Gamma$ , at $s_0=8$ it triples to $3\Gamma$ , while the peak height saturates, so the profile fattens without getting taller. *Verified: the peak value at $\delta=0$ tends to $\Gamma/2$ as $s_0\to\infty$ ; the area grows only sublinearly because height and width both saturate. Extrapolating a measured width to $s_0\to0$ recovers the natural linewidth. See Demtröder, Laser Spectroscopy, Vol. 1.*

Limiting cases

Case	Result	Meaning
Weak drive, on resonance ( $s_0\ll1$ , $\delta=0$ )	$\rho_{ee}\approx\tfrac12 s_0=\dfrac{\Omega^2}{\Gamma^2}$ , $\;R_{sc}\approx\dfrac{\Gamma}{2}s_0\propto I$	Linear Einstein regime — population and fluorescence rise proportionally to intensity; no saturation
Strong drive, on resonance ( $s_0\gg1$ , $\delta=0$ )	$\rho_{ee}\to\tfrac12$ , $\;R_{sc}\to\dfrac{\Gamma}{2}$	Full saturation — equal populations, maximum fluorescence $\Gamma/2$ , no inversion
At the saturation intensity ( $s_0=1$ , $\delta=0$ )	$\rho_{ee}=\tfrac14$ , $\;R_{sc}=\dfrac{\Gamma}{4}$	Definition of $I_{sat}$ — the excited population is exactly one quarter
Far detuning ( $\vert\delta\vert\gg\Gamma,\Omega$ )	$\rho_{ee}\approx\dfrac{\Omega^2}{4\delta^2}$ , $\;R_{sc}\approx\dfrac{\Gamma\Omega^2}{8\delta^2}$	Lorentzian wings — scattering falls as $1/\delta^2$ ; the low-scatter regime of far-detuned dipole traps

⚠️ Strong driving inverts a two-level atom and makes it lase

It does not, and cannot. For a closed two-level system $\rho_{ee}\le1/2$ no matter how hard you drive, because the same field that pumps the atom up stimulates it back down at the identical rate — stimulated absorption and emission share one Einstein $B$ coefficient. The drive can at best equalise the populations; only spontaneous emission is one-way, and it always favours the ground state, capping $\rho_{ee}$ at $1/2$ . Since gain requires $\rho_{ee}>\rho_{gg}$ , a driven two-level atom never inverts and never amplifies — real lasers must use three- or four-level schemes where the pumped and lasing transitions are different. A second trap lives in the symbols: $s_0=I/I_{sat}$ is the on-resonance parameter, while $s=s_0/[1+(2\delta/\Gamma)^2]$ already carries the detuning. Use $s_0$ in the power-broadened width $\Gamma\sqrt{1+s_0}$ and in $I_{sat}$ ; use $s$ only in $\rho_{ee}=\tfrac12 s/(1+s)$ . Swapping them misplaces the saturated linewidth.

Worked examples

1 · Saturating the sodium D₂ line (confidence builder)

The sodium D₂ line has $\lambda=589$ nm and $\Gamma=2\pi\times9.79$ MHz $=6.15\times10^{7}$ s⁻¹. Find its saturation intensity, the excited fraction and scattering rate at $I=I_{sat}$ on resonance, and the maximum spontaneous force and deceleration on an atom of mass $23$ u.

Step 1. Saturation intensity from equation (4): $I_{sat}=\pi hc\Gamma/3\lambda^3$ . Numerically $\pi\times(6.626\times10^{-34})\times(3.00\times10^{8})\times(6.15\times10^{7})/[3\times(589\times10^{-9})^3]=62.6$ W m⁻². Converting, $62.6$ W m⁻² $=$ 6.3 mW cm⁻², the standard book value for the sodium D₂ line.

Step 2. Drive exactly at saturation, $I=I_{sat}$ , so $s_0=1$ and $\delta=0$ . From equation (2), $\rho_{ee}=\tfrac12\cdot\tfrac{1}{1+1}=$ $1/4$ . The scattering rate is $R_{sc}=\Gamma\rho_{ee}=\Gamma/4=1.54\times10^{7}$ photons s⁻¹ — about 15 million photons per second, well below the ceiling because we are only at $s_0=1$ .

Step 3. The absolute ceiling is at $s_0\to\infty$ : $R_{sc}^{max}=\Gamma/2=3.08\times10^{7}$ s⁻¹. With recoil momentum $\hbar k=h/\lambda=1.12\times10^{-27}$ kg m s⁻¹, the maximum force is $F_{max}=\hbar k\,\Gamma/2=(1.12\times10^{-27})(3.08\times10^{7})=$ $3.46\times10^{-20}$ N.

Step 4. With $m=23\times1.66\times10^{-27}=3.82\times10^{-26}$ kg, the peak deceleration is $a_{max}=F_{max}/m=$ $9.1\times10^{5}$ m s⁻², roughly $9\times10^{4}$ times gravity — what lets a Zeeman slower stop a $600$ m s⁻¹ thermal sodium beam in under a metre, the concrete payoff of the $\Gamma/2$ ceiling.

2 · Power broadening and a detuned scattering rate

A cooling laser drives the same sodium transition at $I=8\,I_{sat}$ , so $s_0=8$ . (a) On resonance, find $\rho_{ee}$ , $R_{sc}$ , and the linewidth. (b) Now detune to $\delta=\Gamma$ at the same intensity — find the effective $s$ and the new $\rho_{ee}$ .

Step 1. On resonance with $s_0=8$ : $\rho_{ee}=\tfrac12\cdot\tfrac{8}{1+8}=\tfrac{4}{9}=$ $0.44$ , close to but under the ceiling $1/2$ . The scattering rate is $R_{sc}=\Gamma\rho_{ee}=0.44\,\Gamma=2.7\times10^{7}$ photons s⁻¹.

Step 2. The power-broadened width from equation (5): $\Delta\omega_{FWHM}=\Gamma\sqrt{1+8}=3\Gamma$ , i.e. $3\times9.79=$ $29.4$ MHz FWHM, three times the natural $9.79$ MHz — the line has fattened threefold while the peak barely moved.

Step 3. Detune to $\delta=\Gamma$ . The effective saturation parameter is $s=s_0/[1+(2\delta/\Gamma)^2]=8/[1+4]=$ $1.6$ : the detuning has cut the effective drive from $8$ to $1.6$ .

Step 4. The detuned excited population is $\rho_{ee}=\tfrac12\cdot\tfrac{1.6}{1+1.6}=\tfrac12\times0.615=$ $0.31$ , down from $0.44$ on resonance. Equivalently $R_{sc}=0.31\,\Gamma=1.9\times10^{7}$ s⁻¹ — the atom still scatters strongly a full linewidth off resonance precisely because the line has been power-broadened.

Physics of the portrait

Saturation is one of those results that sits at the crossing of three lineages: Einstein's rate coefficients, Bloch's damped precessing magnetization, and the quantum-optics theory of resonance fluorescence. The number $\rho_{ee}=\tfrac12\,s/(1+s)$ is deceptively small, but the ceiling of $1/2$ it enforces is what makes laser cooling, single-atom detection and Doppler-free spectroscopy possible — and what forbids a two-level laser.

saturation

the ρ_ee = ½ s/(1+s) ceiling of a driven two-level atom · atomic physics and quantum optics

Saturation is where the classical rate-equation picture of Einstein's 1917 coefficients meets the coherent physics of the density matrix. Einstein already held the seed: because stimulated absorption and emission share one $B$ coefficient, no amount of driving can pump more atoms up than down, so a two-level system tends toward equal populations. What his rate equations could not capture — coherence, Rabi flopping, the light shift, the exact line shape — arrived with the optical Bloch equations, direct descendants of the phenomenological equations Felix Bloch wrote in 1946 for a damped, driven nuclear magnetization. Feynman, Vernon and Hellwarth showed in 1957 that any two-level system maps onto a precessing, decaying Bloch vector, handing magnetic-resonance machinery to optics wholesale; set that vector's motion to zero and saturation, power broadening and the scattering rate drop out in a few lines. The picture reached its sharpest form in resonance fluorescence: at strong drive the scattered light is not a simple Lorentzian but a three-peaked spectrum, the Mollow triplet predicted by Benjamin Mollow in 1969 and seen in the mid-1970s. Saturation is the time-averaged face of that same driven atom — the ceiling that survives once the coherent oscillation is averaged away.

⏳ History

1917

— Einstein derives absorption and spontaneous plus stimulated emission rates, the A and B coefficients

1946

— Bloch writes the equations of motion for a damped, driven magnetization in nuclear magnetic resonance

1957

— Feynman, Vernon and Hellwarth recast any two-level system as a precessing, decaying Bloch vector

1969

— Mollow predicts the power spectrum of resonance fluorescence, the three-peaked Mollow triplet

1985

— Chu and co-workers demonstrate optical molasses, cooling atoms with the saturated scattering force

Where this lives today

Every magneto-optical trap and atomic-beam slower runs on this one equation. A MOT's capture velocity, temperature and loading rate are all set by $R_{sc}=\Gamma\rho_{ee}$ and its $\Gamma/2$ ceiling — designers deliberately run near $s_0\sim1$ because pushing brighter buys almost no extra force (the rate has saturated) while power-broadening the line and widening the velocity capture range. Trapped-ion quantum computers read out a qubit by asking whether the ion fluoresces: a "bright" state scatters up to $\Gamma/2\approx$ tens of millions of photons per second, a "dark" state scatters none, and the two are distinguished with well above $99.9\%$ fidelity in under a millisecond — readout speed is literally bounded by the saturated scattering rate. The same ceiling shapes precision spectroscopy and photonics. Saturated-absorption spectroscopy bleaches a velocity class so a counter-propagating probe sees a Lamb dip of power-broadened width $\Gamma\sqrt{1+s_0}$ ; extrapolating to zero intensity recovers the natural linewidth and locks a laser to an atomic line — the backbone of optical clocks. In bulk media the absorption falls as $\alpha=\alpha_0/(1+I/I_{sat})$ because the excited fraction cannot pass $1/2$ , exactly the bleaching a saturable-absorber mirror exploits to mode-lock a femtosecond laser. Wherever light pushes a resonant transition as hard as it will go, the answer is this steady state.

Going further

Texts. Foot, Atomic Physics (2005), Ch. 7 (optical Bloch equations, saturation) · Loudon, The Quantum Theory of Light, 3rd ed., Ch. 2 and 8 (resonance fluorescence) · Metcalf & van der Straten, Laser Cooling and Trapping (saturation intensity, scattering rate, force) · Demtröder, Laser Spectroscopy, Vol. 1 (saturation, power broadening) · Saleh & Teich, Fundamentals of Photonics, 2nd ed., Ch. 14 (rate equations, saturation).

Provenance. The results are the standard steady-state solutions of the optical Bloch equations. A verifier should: (1) reconstruct the rotating-frame Bloch equations and set the derivatives to zero to reproduce $\rho_{ee}^{ss}=(\Omega^2/4)/(\delta^2+\Omega^2/2+\Gamma^2/4)$ (Foot Ch. 7; Loudon 3rd ed. Ch. 2); (2) confirm the saturation-parameter form with $s_0=2\Omega^2/\Gamma^2=I/I_{sat}$ and $I_{sat}=\pi hc\Gamma/3\lambda^3$ (Metcalf & van der Straten Ch. 2), checking $I_{sat}\approx6.3$ mW cm⁻² for the Na D₂ line; (3) verify the power-broadened width $\Gamma\sqrt{1+s_0}$ (Demtröder Vol. 1); (4) check the limiting behaviours against the Einstein rate-equation results, and the worked deceleration $a_{max}=\hbar k\Gamma/2m\approx9\times10^{5}$ m s⁻² for sodium. Atomic parameters (Na D line, ion linewidths) are cited as standard tabulated values.

Prerequisites. Einstein A and B Coefficients · Spontaneous Emission · Stimulated Emission · Radiative Lifetime and Natural Linewidth · Transition Dipole Moment · Lorentzian Line Shape. Leads to. Resonance Fluorescence and the Mollow Triplet · Optical Molasses and the Scattering Force · Saturation Spectroscopy · Rabi Oscillations.

Ask the tutor: Want to see how this steady state turns into the friction force of optical molasses, or how strong driving splits the fluorescence into the Mollow triplet — just ask.

Why can't I just keep raising the intensity to excite the atom more?

Because absorption and stimulated emission are perfectly symmetric — they run off the same Einstein $B$ coefficient. Every photon the field uses to promote the atom, it can equally use to stimulate it back down, so the coherent drive, however strong, only pushes the two populations toward equality, $\rho_{ee}=\rho_{gg}=1/2$ . The single asymmetric process is spontaneous emission, which always returns population to the ground state, so in balance the excited fraction can never quite reach $1/2$ and never exceeds it. Raising the intensity past $I_{sat}$ mostly just power-broadens the line; it barely raises the excited fraction.

What is the difference between s and s₀?

$s_0=I/I_{sat}=2\Omega^2/\Gamma^2$ is the on-resonance saturation parameter — it measures how hard you drive relative to the natural scale, ignoring detuning. $s=s_0/[1+(2\delta/\Gamma)^2]$ is the effective saturation parameter that folds the detuning in: as you tune away from resonance the denominator grows and the effective drive shrinks. Use $s$ in the population formula $\rho_{ee}=\tfrac12\,s/(1+s)$ , because that captures the detuned response. Use $s_0$ in the saturation intensity ( $s_0=I/I_{sat}$ ) and in the power-broadened width $\Gamma\sqrt{1+s_0}$ , both of which are defined at line centre. Mixing them up is the single most common algebra slip in this topic.

Where do the Einstein rate equations fit into this?

They are the weak-drive limit. When $s_0\ll1$ the atom rarely gets excited, coherence effects are negligible, and $\rho_{ee}\approx\tfrac12 s_0\propto I$ — exactly the linear absorption-versus-spontaneous-decay balance Einstein wrote in 1917, and what you recover if you keep only the incoherent (population) terms of the Bloch equations. Saturation is precisely the failure of that linear picture once the drive is strong enough that stimulated emission competes with absorption. The density-matrix treatment contains the rate equations as a limit and adds coherence, the light shift, power broadening and, at very strong drive, the Mollow triplet.

Why does the fluorescence rate saturate but the linewidth keep growing?

Two different quantities are doing two different things. The peak height of the fluorescence line is set by the on-resonance $\rho_{ee}$ , which is capped at $1/2$ , so the peak scattering rate cannot exceed $\Gamma/2$ — it saturates. The width, on the other hand, reflects how far off resonance an atom can still be appreciably driven. As you crank up the intensity, atoms that were previously too detuned to respond now see an effective drive strong enough to scatter, so the wings fill in and the line fattens as $\Gamma\sqrt{1+s_0}$ . Height is bounded by the population ceiling; width is bounded by nothing, so the profile grows sideways while its top stays flat.

Is the scattered light coherent or incoherent?

It depends on the drive. At weak drive ( $s_0\ll1$ ) the atom scatters mostly coherently — light at the laser frequency, phase-locked to it, essentially elastic Rayleigh scattering, a component scaling as $\Omega^2$ . As the drive grows the coherent fraction shrinks: a fully saturated atom carries no steady coherence ( $w=\rho_{ee}-\rho_{gg}\to0$ , so $\tilde\rho_{ge}\to0$ ), and almost all the scattered light is incoherent — spontaneously emitted photons of natural or power-broadened width, which at very strong drive split into the Mollow triplet. So the atom passes from a coherent mirror to an incoherent fluorescent source as you turn up the light.

How does this steady state give rise to a force on the atom?

Each absorbed photon delivers a momentum kick $\hbar k$ along the beam, while the subsequent spontaneous emission recoils in a random direction, so over many cycles the emission kicks average to zero and only the directed absorption kicks survive. The net force is $F=\hbar k\,R_{sc}=\hbar k\,\Gamma\rho_{ee}$ , and because $R_{sc}$ saturates at $\Gamma/2$ the force saturates at $F_{max}=\hbar k\Gamma/2$ — one recoil per two lifetimes, about $3.5\times10^{-20}$ N for sodium and a deceleration near $10^6$ m s⁻². This scattering force is the engine of Zeeman slowers and optical molasses; its velocity dependence through the Doppler shift in $\delta$ provides the cooling friction.