Sub-Doppler Cooling and Sisyphus Mechanism

For a two-level atom the story of laser cooling ends at the Doppler limit TD=ℏΓ/2kBT_D=\hbar\Gamma/2k_B — roughly 235 μ235\,\muK for sodium — where the viscous drag of red-detuned light balances the momentum diffusion of random photon scattering. In 1988 the NIST group released sodium optical molasses and clocked its ballistic expansion at 43 μ43\,\muK — five times colder than a two-level atom is allowed to be, and impossible in Doppler theory. The missing ingredient was the internal structure quantum mechanics hands every real alkali — a manifold of ground-state Zeeman sublevels — together with two effects this article takes as given from the atom–light chapter: the light shift (AC Stark shift) a laser imposes on each sublevel, and optical pumping, the scattering that shuffles population between them. Standing waves supply the third ingredient, a polarization gradient: two counter-propagating beams with orthogonal polarizations build a field whose polarization state, not just its intensity, cycles on the scale of the wavelength. Each ground sublevel then feels a light shift keyed to the local polarization, so the sublevel energies become spatially periodic potentials — corrugated hills and valleys locked to the standing wave, with position-dependent optical pumping tying the atom's internal state to where it sits. The result is the Sisyphus mechanism: a moving atom is repeatedly caught climbing a hill, then pumped back to a well bottom, converting kinetic energy into light-shift energy that departs with a blue-shifted photon. Temperatures plunge toward a small multiple of the recoil energy, and the microkelvin regime — the doorway to Bose–Einstein condensation — swings open.

Key takeaways
  • Sub-Doppler cooling reaches temperatures far below TD=ℏΓ/2kBT_D=\hbar\Gamma/2k_B by exploiting the atom's ground-state Zeeman sublevels, which acquire spatially modulated light shifts inside a polarization gradient.
  • In the lin⟂lin (Sisyphus) geometry two orthogonally polarized counter-propagating beams make the polarization cycle σ+→\sigma^+\to lin →σ−\to\sigma^- every λ/4\lambda/4, modulating the mgm_g sublevel potentials a quarter wavelength out of phase.
  • A moving atom climbs a hill in one sublevel, is optically pumped near the top into a valley of the other, and climbs again — always uphill — dumping kinetic energy into blue-shifted light.
  • The equilibrium temperature scales as kBT∝U0∝s0/∣δ∣k_BT\propto U_0\propto s_0/|\delta| — it falls linearly with intensity, opposite to the intensity-independent Doppler limit, and that scaling was the experimental smoking gun.
  • A distinct σ+σ−\sigma^+\sigma^- mechanism cools with uniform light shifts through a velocity-induced population orientation, giving friction with no Sisyphus hills.
  • Key idea: when an atom has several ground sublevels, a spatially varying polarization modulates their light-shift potentials; optical pumping keeps a moving atom mostly climbing hills, continuously converting kinetic energy into light-shift potential energy that leaves with scattered photons — Sisyphus cooling.
The core insight

When an atom has several ground sublevels, a spatially varying polarization makes their light-shift potentials spatially modulated; a moving atom is optically pumped so that it spends most of its time climbing hills, continuously converting kinetic energy into light-shift potential energy that departs with scattered photons — Sisyphus cooling.

Where this lives today

Every laboratory that makes ultracold atoms uses this mechanism, usually without ceremony. In a magneto-optical trap the magnetic gradient and scattering force capture a cloud, but the final temperature is set by a brief molasses stage: the field is ramped to zero and the beams detuned further red, and polarization-gradient cooling drops the cloud to tens of microkelvin in milliseconds before it is handed to a magnetic or optical dipole trap. That low starting temperature and high phase-space density are what make the subsequent evaporative cooling to Bose–Einstein condensation feasible within a finite trap lifetime. For species where lin⟂lin cooling is inefficient — lithium and potassium, with their compressed excited states — gray molasses on the D1D_1 line uses dark states to combine Sisyphus dissipation with suppressed reabsorption, reaching a few microkelvin. The same physics sets the floor in the caesium and rubidium fountain clocks that define the SI second, and in the ultracold sources feeding atom-interferometer gravimeters and inertial sensors.

The physical picture

Picture the two ground sublevels as interleaved corrugated landscapes, shifted by half a period so every valley of one sits under a hilltop of the other. Optical pumping is a one-way trapdoor: whenever an atom drifts near the top of a hill in the landscape it occupies, the local polarization is exactly the one that pumps it out of that sublevel and into a valley of the other. Having fallen, a moving atom can only climb again. Averaged over many cycles it is almost always going uphill, paying kinetic energy at every rise — the punishment of Sisyphus, run in reverse. The energy budget is transparent: the atom absorbs a photon near a potential maximum and emits one that is, on average, blue-shifted, carrying away U0U_0 plus a recoil. The light field pockets a sliver of the atom's mechanical energy every turn.

Figure
The Sisyphus landscape the two ground sublevels m = +1/2 and −1/2 see light-shift potentials modulated a quarter wavelength out of phase (blue cos²kz, red sin²kz). A moving atom climbs a blue hill losing kinetic energy, is optically pumped near the top (green arrow) into a red valley, then climbs again — always uphill. Depth here U₀ ≈ 12 μK for rubidium at s₀ = 1, δ = −3Γ; polarization cycles σ⁺ → lin → σ⁻ every λ/4
Two sublevel potentials, π/2 out of phase — optical pumping keeps the atom climbing optical pumping U₀ ≈ 12 μK m = +1/2 (cos²kz) m = −1/2 (sin²kz) λ/4 σ⁺ lin σ⁻ lin σ⁺ lin σ⁻ σ⁺

Why so cold? Because the well depth U0∼ℏΩ2/(4∣δ∣)U_0\sim\hbar\Omega^2/(4|\delta|) is a knob we control — large detuning and modest intensity make it small — whereas the Doppler limit ℏΓ/2\hbar\Gamma/2 is fixed by nature. Cooling continues until the kinetic energy is comparable to a well depth, below which the atom can no longer climb out and the engine stalls, so kBT∼U0k_BT\sim U_0, tunable toward the recoil energy. That is a qualitatively new prediction: lower the intensity and the temperature should fall, the fingerprint the second figure shows. The σ+σ−\sigma^+\sigma^- configuration cools by a subtler route with no wells at all, but the lin⟂lin picture is the one to hold in your head.

Figure
The experimental fingerprint sub-Doppler temperature falls linearly with laser intensity, k_B T ∝ U₀ ∝ s₀ — the opposite of the Doppler limit, which is fixed at T_D = 235 μK for sodium no matter the intensity (far above the top of this frame). Cooling floors a few times above the single-photon recoil temperature T_r ≈ 2.4 μK. Points sketch rubidium-like molasses at detuning δ = −3Γ
The smoking gun: T rises linearly with intensity (Doppler theory predicts a constant) 0 10 20 30 40 50 0 1 2 3 saturation parameter s₀ (intensity per beam) temperature T (μK) recoil floor T_r ≈ 2.4 μK Sisyphus: k_B T ∝ U₀ ∝ s₀ ↑ Doppler limit T_D = 235 μK — a flat line, ×6 above this frame

The core equations

Everything begins with the benchmark that sub-Doppler cooling beats: Doppler theory balances the drag of a red-detuned two-level atom against the diffusion of spontaneous emission and finds a floor set only by the linewidth.

TD=ℏΓ2kB(1)T_D = \frac{\hbar\Gamma}{2k_B} \qquad (1)

where: Γ\Gamma — natural linewidth (spontaneous decay rate); ℏ\hbar — reduced Planck constant; kBk_B — Boltzmann constant. The limit depends on nothing tunable — one number per atom (235 μ235\,\muK Na, 146 μ146\,\muK Rb, 125 μ125\,\muK Cs) that two-level theory cannot beat. The heart of Sisyphus cooling is that the sublevels acquire periodic potentials. In the Jg=1/2J_g=1/2 model the two sublevels trace cos⁡2\cos^2 and sin⁡2\sin^2 along the standing wave, offset by λ/4\lambda/4: mg=+1/2m_g=+1/2 shifts most where the field is σ+\sigma^+, mg=−1/2m_g=-1/2 where it is σ−\sigma^-.

U+(z)=U0cos⁡2(kz),U−(z)=U0sin⁡2(kz)(2)U_+(z)=U_0\cos^2(kz),\qquad U_-(z)=U_0\sin^2(kz) \qquad (2)

where: U±U_\pm — light-shift potential of the mg=±1/2m_g=\pm1/2 sublevel; U0U_0 — peak-to-valley depth; k=2π/λk=2\pi/\lambda. A valley of one landscape sits beneath a hilltop of the other — exactly the geometry the trapdoor of optical pumping needs. The depth U0U_0 is the ground-state AC Stark shift. Its exact form saturates logarithmically, but far-detuned and weakly saturated it collapses to a clean form proportional to intensity and inverse detuning.

U0=ℏδ2 ln⁡ ⁣(1+s01+(2δ/Γ)2)  →∣δ∣≫Γ,  s0≪1  ℏΩ24δ(3)U_0 = \frac{\hbar\delta}{2}\,\ln\!\left(1+\frac{s_0}{1+(2\delta/\Gamma)^2}\right)\;\xrightarrow{|\delta|\gg\Gamma,\;s_0\ll1}\;\frac{\hbar\Omega^2}{4\delta} \qquad (3)

where: δ=ωL−ω0\delta=\omega_L-\omega_0 — detuning (negative for cooling); s0=I/Isats_0=I/I_{\rm sat} — saturation parameter per beam; Ω=Γs0/2\Omega=\Gamma\sqrt{s_0/2} — Rabi frequency. Since U0∝s0/δU_0\propto s_0/\delta, the wells shrink at lower power or larger detuning — the mechanism by which temperature falls. The Sisyphus cycle is clocked by optical pumping, which far-detuned is just the photon scattering rate and falls off more steeply with detuning than U0U_0.

Γ′≈Γ s01+(2δ/Γ)2  →∣δ∣≫Γ  ΓΩ24δ2(4)\Gamma' \approx \Gamma\,\frac{s_0}{1+(2\delta/\Gamma)^2}\;\xrightarrow{|\delta|\gg\Gamma}\;\frac{\Gamma\Omega^2}{4\delta^2} \qquad (4)

where: Γ′\Gamma' — pumping/scattering rate, inverse of the cycle time τp\tau_p. The key ratio U0/ℏΓ′∼∣δ∣/ΓU_0/\hbar\Gamma'\sim|\delta|/\Gamma means each pump event carries away a light-shift energy large compared with a recoil — so the cooling is photon-efficient. Balancing the climb against the reset gives a friction linear in velocity whose coefficient, remarkably, does not vanish as intensity is lowered.

F=−αv,α≃− ℏk2 δΓ(∣δ∣≫Γ)(5)F=-\alpha v,\qquad \alpha\simeq-\,\frac{\hbar k^2\,\delta}{\Gamma}\quad(|\delta|\gg\Gamma) \qquad (5)

where: α\alpha — friction coefficient; kk — laser wavenumber. Because α∝ℏk2δ/Γ\alpha\propto\hbar k^2\delta/\Gamma is independent of s0s_0, the drag stays strong even as U0→0U_0\to0 — unlike Doppler friction. The price is a small capture velocity: this strong drag grips only already-slow atoms. The steady-state temperature is a fluctuation–dissipation balance of recoil heating against Sisyphus drag, fixing kBTk_BT at a small multiple of the well depth.

kBT≃Dα  ∝  ℏΩ2∣δ∣  ∝  U0(6)k_BT\simeq\frac{D}{\alpha}\;\propto\;\frac{\hbar\Omega^2}{|\delta|}\;\propto\;U_0 \qquad (6)

where: D∼ℏ2k2Γ′D\sim\hbar^2k^2\Gamma' — momentum-diffusion coefficient; TT — steady-state temperature. Thus kBT∝U0∝s0/∣δ∣k_BT\propto U_0\propto s_0/|\delta| — temperature falls linearly with intensity, the plotted signature and the break from Doppler physics. Cooling cannot run forever: once the wells are shallower than one photon recoil, a single emission kicks the atom out. This sets the floor, a few times the recoil temperature.

Er=ℏ2k22m,Tr=ℏ2k2mkB=2ErkB(7)E_r=\frac{\hbar^2k^2}{2m},\qquad T_r=\frac{\hbar^2k^2}{mk_B}=\frac{2E_r}{k_B} \qquad (7)

where: mm — atomic mass; ErE_r — single-photon recoil energy. Numerically Tr≈2.4 μT_r\approx2.4\,\muK (Na), 0.36 μ0.36\,\muK (Rb), 0.20 μ0.20\,\muK (Cs) — the wall below which sub-recoil methods must take over.

Derivation AThe polarization gradient of the lin-perp-lin standing wave

The goal is to prove the geometric fact everything rests on: two counter-propagating beams with orthogonal linear polarizations build not an intensity standing wave but a field of nearly constant magnitude whose polarization state cycles σ+→\sigma^+\to lin →σ−\to\sigma^- with period λ/2\lambda/2, its pure circular points a quarter wavelength apart — the modulation that creates the Sisyphus hills, here as a consequence of superposition alone.

Step 1. Superpose an x^\hat x-wave in +z+z and a y^\hat y-wave in −z-z of equal amplitude: the polarization is ∝x^ eikz+y^ e−ikz\propto\hat x\,e^{ikz}+\hat y\,e^{-ikz} once e−iωte^{-i\omega t} is stripped off. Why: reducing the field to one zz-dependent complex vector exposes the local ellipse directly.

Step 2. At kz=0kz=0 the vector is x^+y^\hat x+\hat y — equal, in-phase components — linear at +45∘+45^\circ. Why: this reference axis fixes the phase convention.

Step 3. At kz=π/4kz=\pi/4 it becomes x^ eiπ/4+y^ e−iπ/4∝x^−iy^\hat x\,e^{i\pi/4}+\hat y\,e^{-i\pi/4}\propto\hat x-i\hat y, a π/2\pi/2 relative phase — pure circular, a σ−\sigma^- point. Why: a quarter-wave of retardation turns equal linear components into circular light.

Step 4. At kz=π/2kz=\pi/2 the vector is ∝x^−y^\propto\hat x-\hat y (linear at −45∘-45^\circ); at kz=3π/4kz=3\pi/4 it is ∝x^+iy^\propto\hat x+i\hat y, the opposite handedness — a σ+\sigma^+ point. Why: two more steps complete the cycle and show the σ−\sigma^- and σ+\sigma^+ points are λ/4\lambda/4 apart.

Step 5. So the polarization runs σ− (λ/4) \sigma^-\,(\lambda/4)\,lin (λ/4) σ+\,(\lambda/4)\,\sigma^+ at constant magnitude. Since σ+\sigma^+ shifts mg=+1/2m_g=+1/2 most and σ−\sigma^- shifts mg=−1/2m_g=-1/2 most, the two potentials are modulated π/2\pi/2 out of phase. Why: this is the origin of the offset cos⁡2\cos^2 and sin⁡2\sin^2 landscapes — the geometry of the light, not any dynamics, plants the corrugation.

Go deeper · Derivation B · 1 minWhy the Sisyphus friction does not vanish at low intensity
Derivation BWhy the Sisyphus friction does not vanish at low intensity

The aim is to build the friction of equation (5) from an energy-loss argument and expose its most surprising property — it stays large as power drops, contradicting Doppler friction. The strategy: follow one atom through a pumping cycle, find the energy lost climbing before the reset, convert it to a force, then track the intensity scaling so the cancellation shows.

Step 1. A moving atom travels ∼vτp\sim v\tau_p between pumps (τp=1/Γ′\tau_p=1/\Gamma'); optimal operation is kvτp∼1kv\tau_p\sim1, where it climbs of order a quarter well — a height ∼U0\sim U_0 — before being reset. Why: matching travel distance to the well period guarantees the atom is caught near a hilltop, where pumping harvests its energy.

Step 2. The mechanical loss rate is energy per cycle times cycle rate, weighted by the fraction of a well climbed: E˙∼U0k2v2/Γ′\dot E\sim U_0k^2v^2/\Gamma'. Why: "energy per event × events per second" is the cleanest route to a rate; the two velocity factors encode climb height and timing.

Step 3. A loss quadratic in vv is a force linear in vv: from E˙=Fv\dot E=Fv, α∼U0k2/Γ′\alpha\sim U_0k^2/\Gamma'. Why: power equals force times velocity — the definition turning dissipation into a drag coefficient.

Step 4. Insert the scalings U0∝ℏΩ2/δU_0\propto\hbar\Omega^2/\delta and Γ′∝ΓΩ2/δ2\Gamma'\propto\Gamma\Omega^2/\delta^2: the intensity factor Ω2\Omega^2 cancels, leaving α∝ℏk2δ/Γ\alpha\propto\hbar k^2\delta/\Gamma, independent of s0s_0, while the capture velocity vc∼Γ′/kv_c\sim\Gamma'/k shrinks with power. Why: that cancellation is the whole point — Sisyphus friction is strong but narrow, cooling only atoms Doppler cooling has already slowed, so the two act in series.

Go deeper · Derivation C · 1 minThe equilibrium temperature from fluctuation and dissipation
Derivation CThe equilibrium temperature from fluctuation and dissipation

This derivation delivers the central law kBT∝U0∝s0/∣δ∣k_BT\propto U_0\propto s_0/|\delta| by treating the atom as a Brownian particle whose momentum obeys a Langevin equation. The plan: write that equation, extract the steady-state variance by equipartition, insert the diffusion and friction coefficients already built, and combine — so the temperature comes out proportional to the well depth.

Step 1. Model the momentum by mv˙=−αv+ξ(t)m\dot v=-\alpha v+\xi(t) with ⟨ξ(t)ξ(t′)⟩=2D δ(t−t′)\langle\xi(t)\xi(t')\rangle=2D\,\delta(t-t'): deterministic drag plus white-noise recoil kicks. Why: laser cooling is exactly a particle in a viscous bath, so the Langevin equation is the correct minimal description.

Step 2. Multiply by vv and average in steady state (ddt⟨v2⟩=0\frac{d}{dt}\langle v^2\rangle=0): α⟨v2⟩=D/m\alpha\langle v^2\rangle=D/m, so 12m⟨v2⟩=12kBT\tfrac12 m\langle v^2\rangle=\tfrac12 k_BT with kBT=D/αk_BT=D/\alpha. Why: equipartition converts the mechanical variance into the temperature a time-of-flight actually measures.

Step 3. Estimate the diffusion: each scattered photon shifts momentum by ∼ℏk\sim\hbar k at rate Γ′\Gamma', so D∼ℏ2k2Γ′D\sim\hbar^2k^2\Gamma'. Why: recoil diffusion is the irreducible heating any scheme must overcome; it fixes the numerator of TT.

Step 4. Insert D∼ℏ2k2Γ′D\sim\hbar^2k^2\Gamma' and α∼U0k2/Γ′\alpha\sim U_0k^2/\Gamma' into kBT=D/αk_BT=D/\alpha to get ℏ2(Γ′)2/U0\hbar^2(\Gamma')^2/U_0; then use ℏΓ′∼U0Γ/∣δ∣\hbar\Gamma'\sim U_0\Gamma/|\delta| to eliminate Γ′\Gamma', giving kBT≃C U0k_BT\simeq C\,U_0 with C=O(1)C=O(1). Why: the ratio U0/ℏΓ′∼∣δ∣/ΓU_0/\hbar\Gamma'\sim|\delta|/\Gamma collapses two intensity-dependent quantities into one clean result — the falsifiable prediction that halving intensity halves temperature.

Result.

kBT=Dα≃C ℏΩ24∣δ∣=C U0,C=O(1)(8)k_BT=\frac{D}{\alpha}\simeq C\,\frac{\hbar\Omega^2}{4|\delta|}=C\,U_0,\qquad C=O(1) \qquad (8)

Limiting cases

The character of the cooling — and the temperature it reaches — depends on which regime of velocity, intensity, and detuning the atom sits in.

CaseResultMeaning
Doppler regime (fast atoms, v≫Γ/kv\gg\Gamma/k)kBTD=ℏΓ2k_BT_D=\dfrac{\hbar\Gamma}{2}Sub-Doppler friction has a small capture velocity, so fast atoms are cooled only by ordinary Doppler molasses, which pre-cools them to TDT_D before Sisyphus takes over.
Sisyphus regime (slow atoms, low s0s_0)kBT≃C U0∝s0∣δ∣k_BT\simeq C\,U_0\propto \dfrac{s_0}{|\delta|}Temperature falls linearly with intensity, lowered by reducing s0s_0 or raising ∣δ∣|\delta| at fixed U0U_0 (CC order unity) — the sub-Doppler window.
Recoil floorkBTr=ℏ2k2m=2Erk_BT_r=\dfrac{\hbar^2k^2}{m}=2E_rCooling saturates a few times above the recoil temperature (≈ 2.4 μK Na, 0.36 μK Rb, 0.20 μK Cs); a recoil kick ejects atoms from wells shallower than ErE_r.
High-intensity breakdownU0≳ℏΓ⇒T→TDU_0\gtrsim\hbar\Gamma\Rightarrow T\to T_DWhen wells approach the Doppler energy scale, saturation and multi-photon effects wash out the advantage and TT climbs back toward TDT_D.

Worked examples

1 · The benchmark and the anomaly (confidence builder)

Sodium's cooling transition has linewidth Γ/2π=9.79 \Gamma/2\pi=9.79\,MHz and wavelength 589 589\,nm (m=23 m=23\,u). Find its Doppler limit and recoil temperature, and place the 1988 NIST measurement of 43 μ43\,\muK among them.

Step 1. With Γ=2π(9.79×106)=6.15×107 \Gamma=2\pi(9.79\times10^6)=6.15\times10^7\,s⁻¹, the Doppler limit is TD=ℏΓ/2kB=(1.055×10−34)(6.15×107)/[2(1.381×10−23)]=235 μKT_D=\hbar\Gamma/2k_B=(1.055\times10^{-34})(6.15\times10^7)/[2(1.381\times10^{-23})]=\mathbf{235\,\mu K}.

Step 2. With k=2π/λ=1.07×107 k=2\pi/\lambda=1.07\times10^7\,m⁻¹ and m=3.82×10−26 m=3.82\times10^{-26}\,kg, the recoil temperature is Tr=ℏ2k2/(mkB)=2.4 μKT_r=\hbar^2k^2/(mk_B)=\mathbf{2.4\,\mu K} (the single-photon recoil is Er/kB=1.2 μE_r/k_B=1.2\,\muK).

Step 3. The measured 43 μ43\,\muK sits a factor 5.55.5 below TDT_D and ∼18\sim18 above TrT_r — in the forbidden zone for two-level theory yet above the recoil floor. The 1988 result is impossible in Doppler theory (43 μK vs a 235 μK limit) and about 18× above recoil — the fingerprint of a sub-Doppler mechanism.

2 · Well depth and temperature of a rubidium molasses

Rubidium-87 has Γ/2π=6.07 \Gamma/2\pi=6.07\,MHz. Run a molasses at s0=1s_0=1 per beam, δ=−3Γ\delta=-3\Gamma. Estimate the Rabi frequency, well depth U0U_0, pumping time, and expected temperature.

Step 1. Ω=Γs0/2=Γ/2=0.707 Γ\Omega=\Gamma\sqrt{s_0/2}=\Gamma/\sqrt2=0.707\,\Gamma, with Γ=3.81×107 \Gamma=3.81\times10^7\,s⁻¹, so Ω=2.70×107 \Omega=2.70\times10^7\,s⁻¹.

Step 2. From equation (3), U0=ℏΩ2/(4∣δ∣)=ℏΩ2/(12Γ)=ℏΓ/24U_0=\hbar\Omega^2/(4|\delta|)=\hbar\Omega^2/(12\Gamma)=\hbar\Gamma/24 (using Ω2=Γ2/2\Omega^2=\Gamma^2/2). Since ℏΓ/kB=2TD=291 μ\hbar\Gamma/k_B=2T_D=291\,\muK, U0/kB=291/24=12 μKU_0/k_B=291/24=\mathbf{12\,\mu K} — the depth in the first figure.

Step 3. The pumping rate is Γ′≃ΓΩ2/(4δ2)=Γ/72=5.3×105 \Gamma'\simeq\Gamma\Omega^2/(4\delta^2)=\Gamma/72=5.3\times10^5\,s⁻¹, so τp=1.9 μs\tau_p=\mathbf{1.9\,\mu s}. The capture velocity is vc∼Γ′/k=6.6 cm/sv_c\sim\Gamma'/k=\mathbf{6.6\,cm/s} — tiny beside the Doppler capture speed Γ/k≈4.7 \Gamma/k\approx4.7\,m/s.

Step 4. With kBT≃C U0k_BT\simeq C\,U_0 (CC of order unity), the predicted TT is a few times 12 μ12\,\muK — tens of microkelvin, below the 146 μ146\,\muK Rb Doppler limit and above the 0.36 μ0.36\,\muK recoil floor, matching real Rb molasses. Halving the intensity would halve U0U_0 and halve TT — the linear fall a constant TDT_D cannot produce. U0/kB≈12 μU_0/k_B\approx12\,\muK, τp≈1.9 μ\tau_p\approx1.9\,\mus, vc≈6.6 v_c\approx6.6\,cm/s; T∝s0T\propto s_0.

Physics of the portrait

Sub-Doppler cooling is where a puzzle became a Nobel Prize: a 1988 anomaly met, within a year, a theory built from the humblest atomic-physics ingredients — and the person most identified with that synthesis had spent his career turning atom–photon interactions into a language precise enough to predict it.

Origins · 1 minClaude Cohen-Tannoudji and how this was discovered
CC
Claude Cohen-Tannoudji
b. 1933 · polarization-gradient cooling and atom–photon interactions

Claude Cohen-Tannoudji, born in Constantine, Algeria, in 1933, built much of the modern language of atom–light interaction — the dressed-atom picture, the theory of light shifts, and the master-equation treatment of optical pumping developed at the École Normale Supérieure in the lineage of Alfred Kastler and Jean Brossel. When the NIST group reported sodium molasses far below the Doppler limit in 1988, two-level theory forbade the result: a genuine crisis. Within a year Cohen-Tannoudji and his student Jean Dalibard published the polarization-gradient (Sisyphus) mechanism in the Journal of the Optical Society of America B (1989), showing that spatially varying light shifts plus optical pumping among ground sublevels reproduce the observed temperatures and their T∝s0/∣δ∣T\propto s_0/|\delta| scaling; Steven Chu's group reached the same understanding independently that year. For laser cooling and trapping he shared the 1997 Nobel Prize with Chu and William Phillips. His dressed-atom methods underlie the whole toolkit of ultracold physics — optical lattices, dipole traps, and the sub-recoil schemes that push below the very floor this article ends at.

History
1975
— Hänsch and Schawlow propose Doppler cooling of neutral atoms; Wineland and Dehmelt propose it for trapped ions
1985
— Chu, Hollberg, Bjorkholm, Cable and Ashkin demonstrate three-dimensional optical molasses at Bell Labs
1988
— Lett, Watts, Westbrook and Phillips at NIST measure sodium molasses at 43 μK, far below the 235 μK Doppler limit
1989
— Dalibard and Cohen-Tannoudji, and independently Ungar, Weiss, Riis and Chu, explain polarization-gradient (Sisyphus) cooling
1997
— Chu, Cohen-Tannoudji and Phillips share the Nobel Prize in Physics for laser cooling and trapping

Going further

Texts. Foot, Atomic Physics (2005), ch. 9 · Metcalf and van der Straten, Laser Cooling and Trapping (1999), ch. 8 · Cohen-Tannoudji, Dupont-Roc and Grynberg, Atom–Photon Interactions (1992).

Primary sources. J. Dalibard and C. Cohen-Tannoudji, J. Opt. Soc. Am. B 6, 2023 (1989) — the polarization-gradient theory · P. J. Ungar, D. S. Weiss, E. Riis and S. Chu, J. Opt. Soc. Am. B 6, 2058 (1989) · P. D. Lett et al., Phys. Rev. Lett. 61, 169 (1988) — the sodium sub-Doppler observation.

Provenance. All equations follow the standard treatment (Dalibard and Cohen-Tannoudji 1989; Ungar, Weiss, Riis and Chu 1989; Foot ch. 9; Metcalf and van der Straten). A verifier should confirm the lin⟂lin cycle with σ±\sigma^\pm points λ/4\lambda/4 apart (Derivation 1), the scaling kBT∝U0∝s0/∣δ∣k_BT\propto U_0\propto s_0/|\delta| (Derivation 3), the Doppler limits 235 μK (Na), 146 μK (Rb), 125 μK (Cs) and recoil temperatures 2.4 μK (Na), 0.36 μK (Rb), 0.20 μK (Cs) from the linewidths and wavelengths, and the well depth U0/kB=ℏΓ/24kB≈12 μU_0/k_B=\hbar\Gamma/24k_B\approx12\,\muK. Order-unity constants are configuration-dependent and omitted rather than fabricated. (Some summaries mislabel sodium's Doppler limit as 146 μK, which is actually rubidium's.)

Prerequisites. Doppler Cooling and the Doppler Limit · Optical Molasses · The Dipole Force and the Optical Dipole Potential · Radiation Pressure and the Scattering Force. Leads to. The Recoil Limit and Sub-Recoil Cooling · The Magneto-Optical Trap · Evaporative Cooling · Bose–Einstein Condensation in Dilute Gases.

Ask the tutor: Ask me to walk through why the Sisyphus friction survives at vanishing laser power while Doppler friction dies, or to work the well depth and predicted temperature for your own atom, detuning, and intensity.

Why does the atom always seem to be climbing rather than descending?

That asymmetry is the whole mechanism. Optical pumping is not random in space: it is fastest where the light shift is largest, and delivers the atom to the sublevel whose well is deepest there. So an atom is preferentially pumped near a hilltop of its current potential — having climbed there — and lands in a valley of the other, whereupon it must climb again. The correlation between position, internal state, and pumping rate, engineered by the polarization gradient, biases the atom uphill on average. Remove the spatial modulation and the bias vanishes.

Where does the dissipated energy actually go?

Into the light field, as a tiny blue shift on the scattered photons. In each cycle the atom absorbs a photon near a potential maximum, where its light-shift energy is high, and emits one as it is pumped down; on average the emitted photon carries slightly more energy than the absorbed one — the difference being the U0U_0 the atom gave up climbing, plus a recoil. Summed over the cloud the fluorescence is imperceptibly blue-shifted relative to the drive: that spectral asymmetry is the energy budget balancing, mechanical energy leaving with the light.

If Sisyphus friction is so strong, why not skip Doppler cooling?

Because the strong friction grips only atoms that are already slow. Its capture velocity is one well width per pumping time — a few cm/s in the rubidium example, versus several m/s for the Doppler range. Atoms straight from an oven or a Zeeman slower move far too fast; they sail through the wells before pumping catches them near a hilltop. Doppler cooling supplies the wide velocity capture to decelerate them into the narrow window where polarization-gradient cooling then takes over — broad-but-weak, then narrow-but-strong.

The σ⁺σ⁻ configuration has uniform light shifts — how can it cool?

Through a velocity-induced orientation rather than a Sisyphus hill. Two counter-propagating σ+\sigma^+ and σ−\sigma^- beams make a field of constant magnitude whose linear polarization twists helically along zz. The light shifts are position-independent, so a stationary atom settles into a symmetric distribution with no net force. But a moving atom sees the polarization axis rotating in its frame; if pumping cannot keep the populations adiabatically symmetric, a net orientation ⟨m⟩∝v\langle m\rangle\propto v builds up, making it scatter preferentially from the beam it moves against. The result is damping with no potential wells — a distinct mechanism needing a Jg=1J_g=1 or higher ground state.

Why can't we reach zero temperature by lowering the intensity?

Two things break at low intensity. First, the recoil floor: as U0U_0 falls below the recoil energy ErE_r, one spontaneous emission kicks the atom clean out of a well, so cooling stalls at a few times TrT_r — about 0.36 μ0.36\,\muK for rubidium. Second, the capture velocity vc∼Γ′/kv_c\sim\Gamma'/k also shrinks, so fewer atoms are slow enough to be gripped and residual scattering heats them. The linear law T∝s0T\propto s_0 is only the ideal trend; going lower needs sub-recoil schemes such as velocity-selective coherent population trapping or Raman sideband cooling.

How deep should the wells be for cooling to work well?

There is a comfortable window, Er≲U0≲ℏΓE_r\lesssim U_0\lesssim\hbar\Gamma. Shallower than the recoil energy and a single recoil ejects the atom — recoil-limited. Deeper than the Doppler scale (U0≳ℏΓU_0\gtrsim\hbar\Gamma) and saturation and multi-photon processes take over, the two-sublevel picture fails, and TT climbs back toward TDT_D. In practice one sits in the middle: modest saturation and a detuning of a few linewidths, giving wells tens of times the recoil energy but well below ℏΓ\hbar\Gamma — the s0∼1s_0\sim1, δ∼−3Γ\delta\sim-3\Gamma regime of the worked example.

Before this

Where this leads