Sub-Doppler Cooling

Sub-Doppler cooling is any laser-cooling mechanism whose stationary motional distribution lies below the two-level Doppler scale

$$k_{\mathrm B}T_D = \frac{\hbar\Gamma}{2}.$$

The name describes an achieved regime, not a unique force. Several distinct pieces of physics can evade the assumptions behind the Doppler result:

• polarization gradients make degenerate ground-state sublevels experience different position-dependent light shifts;
• optical pumping can reset an atom preferentially after it has climbed a light-shift potential;
• coherent superpositions can be dark only for selected velocity classes;
• Raman transitions can address kinetic-energy classes without relying on the optical linewidth;
• confinement can resolve discrete motional sidebands and make the vibrational ground state dark to a cooling transition.

These mechanisms do not violate recoil heating or detailed balance. They change where, when, and for which states scattering occurs.

Canonical Scope

This page owns the conceptual and quantitative bridge from two-level Doppler cooling to four families of sub-Doppler physics:

1. polarization gradients in multilevel atoms;
2. Sisyphus cooling by light shifts and optical pumping;
3. dark-state and gray-molasses cooling;
4. free-particle Raman cooling and trapped-particle Raman sideband cooling.

Doppler Cooling owns the friction–diffusion derivation of $T_D$, including its detuning and saturation dependence. AC Stark Shift owns the general derivation of conservative light shifts. Electromagnetically Induced Transparency owns the full susceptibility and dressed-state treatment of a driven $\Lambda$ system. Raman Spectroscopy owns two-photon resonance conventions and effective Raman couplings.

Optical Tweezers specializes Raman-sideband and imaging-compatible cooling to individual particles and arrays. Later pages treat optical lattices, ion traps, and platform-specific control. Here they appear only far enough to explain why a motional ladder can replace velocity as the resolved cooling coordinate.

What the Doppler model leaves out

The standard Doppler limit assumes a broad, closed, nondegenerate two-level transition and a force obtained from the velocity dependence of its scattering rate. Real alkali ground manifolds instead contain Zeeman and hyperfine sublevels. The coupling strength depends on local polarization, and repeated scattering optically pumps population among those sublevels.

Three omissions are decisive.

Internal degeneracy

For total angular momentum $F>0$, the ground manifold contains $2F+1$ states $|F,m_F\rangle$. A local spherical polarization component $q=0,\pm1$ couples them with Clebsch–Gordan amplitude

$$\Omega_{m_F,q}(\mathbf r) \propto \mathcal E_q(\mathbf r) \langle F',m_F+q|d_q|F,m_F\rangle.$$

Different magnetic sublevels therefore see different light shifts and different optical-pumping rates.

Spatial polarization structure

Counterpropagating beams can have constant total intensity while their local polarization changes on the scale of a wavelength. An atom moving through that field samples a position-dependent internal Hamiltonian even when the scalar intensity is uniform.

Coherence and quantized motion

The internal steady state need not be an incoherent population. A coherent superposition can decouple from the light, and confinement can quantize motion into states $|n\rangle$. Neither effect is represented by a two-level rate equation for a continuous velocity.

Three different ways to suppress the Doppler floor. In polarization-gradient cooling, optical pumping occurs after an atom climbs a light-shift potential. Dark-state cooling suppresses scattering in a narrow velocity class. Raman sideband cooling removes one trapped motional quantum per cycle and leaves $|\downarrow,0\rangle$ dark.

Scale hierarchy

The relevant energy and frequency scales should be written down before a mechanism is named:

$$\begin{aligned} E_r &= \frac{(\hbar k)^2}{2m}, & T_r &= \frac{E_r}{k_{\mathrm B}}, \\ E_D &= k_{\mathrm B}T_D = \frac{\hbar\Gamma}{2}, & U_{\mathrm{ls}} &\sim \frac{\hbar|\Omega|^2}{4|\Delta|}, \\ \hbar\gamma_p &\sim \hbar\Gamma \frac{|\Omega|^2}{4\Delta^2}, & \hbar\omega_t &= \text{trap level spacing}. \end{aligned}$$

Here $E_r$ is the one-photon recoil energy, $U_{\mathrm{ls}}$ is a representative differential light-shift depth, $\gamma_p$ is an optical-pumping rate in the far-detuned scaling regime, and $\omega_t$ is a trap frequency. Numerical factors and Clebsch–Gordan coefficients are transition dependent.

For a broad optical transition,

$$E_r \ll \hbar\Gamma.$$

This separation leaves a large interval between the recoil and Doppler scales. Polarization-gradient cooling commonly occupies that interval. When the distribution approaches a few recoils, momentum changes can no longer be treated as continuous diffusion; dark-state or motional-state descriptions become essential.

Polarization gradients

Constant intensity does not imply a uniform coupling

Consider one-dimensional lin-perpendicular-lin molasses. Two equal-frequency beams propagate along $\pm z$ with orthogonal linear polarizations:

$$\mathbf E^{(+)}(z,t) = E_0 \left( \mathbf e_x e^{ikz} + \mathbf e_y e^{-ikz} \right) e^{-i\omega_L t}.$$

Because $\mathbf e_x\cdot\mathbf e_y=0$, the cycle-averaged intensity is independent of $z$:

$$I(z) \propto \mathbf E^{(-)}\cdot\mathbf E^{(+)} = 2|E_0|^2.$$

The normalized Stokes parameters, however, rotate:

$$s_2(z) = \cos(2kz), \qquad s_3(z) = \sin(2kz).$$

Equivalently, up to the convention for handedness,

$$I_{\sigma^\pm}(z) = \frac{I}{2} \left[ 1 \pm \sin(2kz) \right].$$

The polarization changes from linear to circular, back to linear, and then to the opposite circular polarization over a distance $\lambda/2$. The intensity is flat while the couplings between magnetic sublevels are not.

Position-dependent light shifts

For atom-minus-laser detuning

$$\Delta = \omega_0-\omega_L,$$

red detuning has $\Delta>0$. In a far-detuned, low-saturation treatment, the ground-state light shift of sublevel $m$ is schematically

$$U_m(z) \simeq - \frac{\hbar}{4\Delta} \sum_q \left| \Omega_{m,q}(z) \right|^2.$$

The sum includes polarization amplitudes, Clebsch–Gordan coefficients, and, when necessary, several excited hyperfine manifolds. The common scalar part does not cool by itself. Cooling comes from differential vector and tensor shifts that make the $U_m(z)$ curves differ.

The same off-resonant coupling produces optical pumping. A representative far-detuned rate scales as

$$\gamma_{m\rightarrow m'} \sim \Gamma \sum_q \frac{ |\Omega_{m,q}|^2 }{ 4\Delta^2 } B_{m,q\rightarrow m'},$$

where $B_{m,q\rightarrow m'}$ is a branching factor. Near resonance, the denominator must retain $\Gamma^2/4$, and interference among excited hyperfine levels may matter.

Sisyphus cooling

The energy ledger

A minimal Sisyphus cycle has four steps:

1. the atom occupies a magnetic sublevel whose light shift rises along its trajectory;
2. kinetic energy is converted into light-shift potential energy;
3. near the crest, optical pumping transfers the atom to a sublevel with a lower local light shift;
4. spontaneous emission exports the energy difference, and the cycle repeats.

If the atom climbs through a differential depth $U_0$ before pumping, the mean mechanical energy removed in a successful cycle is of order

$$\Delta E_{\mathrm{cool}} \sim U_0.$$

This is why the achievable scale is governed by light shifts rather than directly by $\hbar\Gamma$. Recoil from the pumping photon still heats the motion. The net change per cycle is only favorable when

$$\langle \Delta E_{\mathrm{cool}}\rangle > \langle \Delta E_{\mathrm{recoil}}\rangle + \langle \Delta E_{\mathrm{technical}}\rangle.$$

The spontaneous photon is not an incidental nuisance: it is the irreversible reset that makes repeated downhill transfers possible.

Why timing matters

Let the characteristic polarization period be of order $1/k$. An atom moving at speed $v$ traverses it at rate $kv$. Optical pumping acts at rate $\gamma_p$. Efficient Sisyphus cycles require comparable scales:

$$kv \sim \gamma_p.$$

If $kv\gg\gamma_p$, the atom crosses many hills before its internal state changes. If $kv\ll\gamma_p$, it is pumped before climbing appreciably. This gives the rough capture scale

$$v_c \sim \frac{\gamma_p}{k},$$

supplemented by the energetic condition that the kinetic energy not greatly exceed the useful potential depth:

$$\frac{1}{2}mv^2 \lesssim U_0.$$

These are scaling tests, not universal equalities.

Friction and diffusion

At sufficiently low velocity, a polarization-gradient force can again be written

$$F(v) \simeq -\alpha_{\mathrm{PG}}v.$$

But neither $\alpha_{\mathrm{PG}}$ nor the momentum-diffusion coefficient is given by the two-level Doppler expressions. Both must be derived from the multilevel optical Bloch or master equation, including position-dependent couplings and recoil.

In a broad semiclassical Sisyphus regime, a common scaling is

$$k_{\mathrm B}T \sim C\,U_0,$$

where $C$ depends on angular momenta, detuning, polarization geometry, and dimensionality. This scaling explains the often-observed approximately linear dependence of temperature on intensity at low saturation. It should not be promoted to a universal temperature formula.

Lin-perpendicular-lin and circular configurations

The potential-hill picture is clearest for lin-perpendicular-lin light. In a counterpropagating circular-polarization geometry, the local polarization may have constant ellipticity while its orientation or internal-state basis varies. Cooling then involves motion-induced coherences and population transfer among Zeeman sublevels. It is still called polarization-gradient cooling, but a literal scalar hill-and-reset cartoon can be misleading.

The robust statement is:

Motion through a spatially varying multilevel coupling creates a lag between internal-state adaptation and the local dressed-state energy.

That lag can produce friction below the two-level Doppler scale.

Magnetic-field sensitivity

Sub-Doppler coherences and optical pumping operate on scales far smaller than the optical linewidth. A residual magnetic field introduces a Larmor frequency

$$\omega_B = \frac{|g_F|\mu_{\mathrm B}B}{\hbar}.$$

Cooling is strongly perturbed when $\omega_B$ becomes comparable to the optical-pumping rate or differential light-shift frequency:

$$\omega_B \gtrsim \min \left( \gamma_p, \frac{U_0}{\hbar} \right).$$

This criterion explains why a field tiny compared with the field required to Zeeman-shift an optical transition by $\Gamma$ can still destroy sub-Doppler cooling.

Dark-state cooling

The Lambda-system dark state

Consider two ground states coupled to a common excited state:

$$|g_1\rangle \xleftrightarrow{\ \Omega_1\ } |e\rangle \xleftrightarrow{\ \Omega_2\ } |g_2\rangle.$$

In the rotating-wave approximation, the interaction is

$$H_{\mathrm{int}} = \frac{\hbar}{2} \left( \Omega_1|e\rangle\langle g_1| + \Omega_2|e\rangle\langle g_2| + \text{h.c.} \right).$$

At two-photon resonance, define

$$\Omega = \sqrt{ |\Omega_1|^2 + |\Omega_2|^2 }.$$

The superposition

$$|D\rangle = \frac{ \Omega_2|g_1\rangle - \Omega_1|g_2\rangle }{ \Omega }$$

obeys

$$H_{\mathrm{int}}|D\rangle = 0.$$

It is dark because the two excitation amplitudes cancel, not because either beam is weak. The orthogonal bright state couples to $|e\rangle$ and scatters.

Velocity selection

For wave vectors $\mathbf k_1$ and $\mathbf k_2$, motion changes the two-photon detuning:

$$\delta_2(\mathbf v) = \delta_2(\mathbf 0) - \Delta\mathbf k\cdot\mathbf v, \qquad \Delta\mathbf k = \mathbf k_1-\mathbf k_2.$$

Only a selected velocity class satisfies the dark resonance. If its effective half-width is $\gamma_D$, the corresponding velocity width scales as

$$\Delta v_D \sim \frac{\gamma_D}{|\Delta\mathbf k|}.$$

A useful schematic scattering profile is

$$R_{\mathrm{sc}}(v) \simeq R_b \frac{ \delta_2(v)^2 }{ \delta_2(v)^2+\gamma_D^2 },$$

where $R_b$ is the bright-state scattering scale. Atoms outside the dark window scatter and receive random recoil; some are transferred into the dark velocity class, where scattering is suppressed. Cooling is therefore accumulation in a low-scattering state rather than a simple linear drag.

Ground-state decoherence, differential AC Stark shifts, magnetic-field inhomogeneity, laser phase noise, and collisions all broaden $\gamma_D$ and raise the final kinetic scale.

Velocity-selective coherent population trapping

Velocity-selective coherent population trapping, or VSCPT, treats internal and translational states together. In an ideal symmetric one-dimensional model, the dark state near zero mean momentum has the form

$$|D_0\rangle \propto |g_1,-\hbar k\rangle - |g_2,+\hbar k\rangle.$$

The two components absorb into the same excited momentum state, so their amplitudes cancel. Momentum classes away from resonance are bright and are recycled by spontaneous emission.

The resulting distribution can have narrow peaks and long nonthermal tails. Near or below the recoil scale, quoting a single temperature may conceal the relevant physics. Report the momentum distribution, rms width, kinetic energy, and dark-state population.

Gray molasses

Gray molasses combines dark-state protection with polarization-dependent bright-state light shifts. Spatially varying bright and dark eigenstates make moving atoms nonadiabatically enter a bright state, climb a repulsive light-shift landscape, and optically pump back into a dark state. Slow atoms remain dark for longer and scatter less.

In a $\Lambda$-enhanced implementation, two optical frequencies address ground hyperfine states through a common excited manifold. Performance is sharply sensitive to the two-photon Raman detuning. The method is especially valuable when unresolved excited structure or small hyperfine splittings make ordinary polarization-gradient cooling inefficient.

Gray molasses is not “Sisyphus cooling without spontaneous emission.” Scattering remains the dissipative reset; the adjective gray means that the low-energy state is weakly coupled rather than uniformly bright.

Raman cooling of free particles

Two-photon Raman transitions can select velocity much more narrowly than an optical linewidth. For a transition that changes momentum by $\hbar\Delta\mathbf k$, energy conservation gives

$$\omega_1-\omega_2 = \omega_{21} + \Delta\mathbf k\cdot\mathbf v + \frac{ \hbar|\Delta\mathbf k|^2 }{ 2m }.$$

The last term is the two-photon recoil shift. A cooling sequence uses:

1. a velocity-selective Raman pulse that transfers atoms in a chosen moving class to another internal state and changes their momentum toward zero;
2. optical pumping that returns the internal state;
3. a pulse schedule whose spectral response leaves a narrow dark interval around zero velocity.

The coherent step can remove more directed kinetic energy than the optical pumping step adds on average. By narrowing the unaddressed velocity window over successive pulses, free-particle Raman cooling can produce subrecoil, strongly non-Gaussian distributions.

This is distinct from Raman sideband cooling. Free Raman cooling selects a continuum of velocities. Raman sideband cooling selects discrete vibrational quantum numbers in a trap.

Raman sideband cooling

Resolved motional quanta

Approximate a trapped direction as a harmonic oscillator:

$$H_{\mathrm{mot}} = \hbar\omega_t \left( a^\dagger a + \frac{1}{2} \right).$$

Two long-lived internal states form $|\downarrow,n\rangle$ and $|\uparrow,n\rangle$. A Raman pair with effective wave-vector difference $\Delta\mathbf k$ has Lamb–Dicke parameter

$$\eta = \Delta k\,x_0, \qquad x_0 = \sqrt{ \frac{\hbar}{2m\omega_t} }.$$

In the Lamb–Dicke regime,

$$\eta\sqrt{2n+1} \ll 1,$$

recoil rarely changes the motional state by more than one quantum.

Cooling cycle

Tune the Raman difference frequency to the red motional sideband:

$$|\downarrow,n\rangle \longrightarrow |\uparrow,n-1\rangle.$$

For small $\eta$, its Rabi frequency scales as

$$\Omega_{n,n-1} \simeq \eta\Omega_0\sqrt n.$$

Optical pumping then resets the internal state,

$$|\uparrow,n-1\rangle \longrightarrow |\downarrow,n-1\rangle,$$

while approximately preserving $n$ in the Lamb–Dicke regime. Each ideal cycle removes $\hbar\omega_t$. Because

$$\Omega_{0,-1} = 0,$$

$|\downarrow,0\rangle$ is dark to the red sideband.

The sidebands must be spectrally resolved:

$$\omega_t \gg \Gamma_{\mathrm{eff}}, \Omega_{\mathrm{sb}},$$

where $\Gamma_{\mathrm{eff}}$ is the effective linewidth of the Raman or repumping cycle. Otherwise the carrier and blue sideband are driven and heating competes directly with cooling.

Rate model and endpoint

Let $A_-$ and $A_+$ be the single-quantum cooling and heating coefficients. For a thermal-like motional state,

$$\frac{d\bar n}{dt} = - \left( A_--A_+ \right) \bar n + A_+.$$

If $A_->A_+$, the stationary occupation is

$$\bar n_{\mathrm{ss}} = \frac{A_+}{A_--A_+}.$$

The endpoint is therefore set by off-resonant carrier and blue-sideband excitation, repumping recoil, trap-frequency noise, anharmonicity, and background heating. “Ground-state cooled” should be supported by a measured $\bar n$ or ground-state probability, not inferred merely from a narrow spatial image.

Sideband thermometry

For a thermal oscillator,

$$\bar n = \frac{1}{ \exp( \hbar\omega_t/k_{\mathrm B}T ) -1 }, \qquad P_0 = \frac{1}{1+\bar n}.$$

In the weak-probe limit, the integrated red- and blue-sideband strengths obey

$$\frac{I_{\mathrm{red}}}{I_{\mathrm{blue}}} = \frac{\bar n}{\bar n+1} = \exp \left( - \frac{\hbar\omega_t}{k_{\mathrm B}T} \right).$$

This asymmetry is a direct motional thermometer when state preparation, probe saturation, line overlap, and detection biases are controlled.

Worked scale audit for rubidium

For the $^{87}\mathrm{Rb}$ D$_2$ line, take

$$\lambda \simeq 780\,\mathrm{nm}, \qquad \frac{\Gamma}{2\pi} \simeq 6.07\,\mathrm{MHz}, \qquad m \simeq 1.443\times10^{-25}\,\mathrm{kg}.$$

The two-level Doppler temperature is

$$T_D = \frac{\hbar\Gamma}{2k_{\mathrm B}} \simeq 146\,\mu\mathrm K.$$

The one-photon recoil scales are

$$\begin{aligned} v_r &= \frac{\hbar k}{m} \simeq 5.9\,\mathrm{mm\,s^{-1}}, \\ \frac{E_r}{h} &\simeq 3.77\,\mathrm{kHz}, \\ T_r &= \frac{E_r}{k_{\mathrm B}} \simeq 181\,\mathrm{nK}. \end{aligned}$$

Thus

$$\frac{T_D}{T_r} \simeq 8.0\times10^2.$$

A measured temperature of $10\,\mu\mathrm K$ is deeply sub-Doppler while still about $55$ recoil temperatures. A semiclassical polarization-gradient description may be appropriate there. A distribution with width of order one recoil demands a quantum momentum treatment, and a trapped sample with $\bar n<1$ demands resolved motional-state language. “Sub-Doppler,” “subrecoil,” and “motional ground state” are not synonyms.

Mechanism comparison

Mechanism	Selected coordinate	Dissipative reset	Characteristic endpoint	Necessary signature
two-level Doppler	velocity through optical detuning	spontaneous emission	$\hbar\Gamma$ scale	damping plus recoil diffusion
polarization-gradient Sisyphus	position and magnetic sublevel	optical pumping	differential light-shift scale	polarization and field dependence
gray molasses	local bright/dark dressed state	weak bright-state scattering	dark-resonance width and recoil	sharp Raman-detuning feature
VSCPT	internal–momentum superposition	recycling of bright momentum classes	subrecoil width possible	non-Gaussian dark momentum peak
free Raman cooling	continuum velocity class	optical pumping between coherent pulses	pulse-defined dark window	velocity-selective Raman spectrum
Raman sideband cooling	trapped vibrational number	internal-state repumping	$\bar n\ll1$ possible	resolved red/blue sideband asymmetry

The table is a diagnostic map, not a ranking. A practical experiment often uses several rows in sequence: broad-line Doppler capture, polarization- gradient or gray-molasses cooling, trap loading, and then sideband cooling.

Experimental validation

Demonstrate that the sample is actually colder

Use at least one calibrated motional observable:

• time-of-flight expansion;
• release–recapture;
• Doppler or Raman spectroscopy;
• velocity-selective recoil measurements;
• sideband asymmetry;
• spatial width in a calibrated harmonic trap.

At recoil scale, fit the full distribution rather than forcing a Gaussian temperature onto narrow peaks and broad wings.

Vary the mechanism-specific controls

A sub-Doppler attribution should survive discriminating tests:

• reverse or remove the polarization gradient;
• scan intensity and test light-shift scaling;
• apply a controlled magnetic field;
• scan the two-photon Raman detuning;
• vary the dark time or coherent pulse width;
• change trap frequency and verify resolved-sideband scaling;
• interrupt repumping and test whether cooling stops;
• measure heating after the cooling light is removed.

Account for particle loss and selection

Cooling can be mimicked by preferential loss of hot particles. Report atom number, phase-space density, capture fraction, and the initial and final distributions. A narrower surviving cloud is not by itself evidence of entropy removal.

Failure modes

Calling every low temperature Sisyphus cooling

A temperature below $T_D$ establishes that the two-level Doppler model is insufficient. It does not identify the replacement mechanism. Test polarization, magnetic-field, intensity, and Raman-detuning dependence.

Treating the Doppler limit as fundamental

$\hbar\Gamma/2$ is the equilibrium scale of a declared two-level friction–diffusion model. Multilevel light shifts and dark resonances remove its assumptions.

Treating the recoil temperature as an absolute floor

Continual bright-state scattering tends to produce recoil-scale heating, but coherent dark states can suppress scattering for the coldest atoms. Subrecoil distributions are possible. Their quantum and nonthermal nature must be described explicitly.

Ignoring excited hyperfine structure

If the detuning is not small compared with excited-state hyperfine splittings, several pathways interfere. A single isolated-$F'$ model can predict the wrong light shifts, dark states, or cooling sign.

Using a scalar AC Stark shift

Polarization-gradient cooling depends on differential vector and tensor light shifts and optical pumping. A scalar potential alone cannot encode the internal-state cycle.

Forgetting the magnetic-field scale

Comparing a Zeeman shift only with $\Gamma$ misses the relevant sensitivity. Compare the Larmor frequency with $\gamma_p$, $U_0/\hbar$, and the dark-state linewidth.

Confusing Raman cooling with Raman sideband cooling

Free Raman cooling selects velocity classes. Sideband cooling requires a resolved trapped ladder and a Lamb–Dicke parameter. The word Raman names the coherent coupling, not the motional regime.

Inferring ground-state cooling from temperature alone

For a quantized trap, the decisive variables are $\bar n$, $P_0$, and the mode-resolved distribution. A classical temperature can be ambiguous in a nonthermal or anisotropic state.

Analysis workflow

1. State the internal manifold. List $F$, $F'$, Zeeman states, branching, and repump pathways.
2. Write the optical geometry. Include wave vectors, phases, polarizations, and intensity imbalance.
3. Fix detuning conventions. Distinguish one-photon and two-photon detunings.
4. Compute local dressed states. Identify light shifts, bright states, and dark states.
5. Compare motion with pumping. Evaluate $kv$, $\gamma_p$, $U_0/\hbar$, and magnetic precession.
6. Choose the motional description. Use semiclassical diffusion only when many recoils occupy the distribution; otherwise retain momentum or trap states.
7. Include all heating channels. Recoil, off-resonant excitation, technical noise, collisions, and trap heating belong in the endpoint.
8. Predict a discriminating scan. A mechanism claim should imply a measurable dependence that a competing mechanism does not.

References

1. P. D. Lett, R. N. Watts, C. I. Westbrook, W. D. Phillips, P. L. Gould, and H. J. Metcalf, “Observation of atoms laser cooled below the Doppler limit,” Physical Review Letters 61, 169–172 (1988), doi:10.1103/PhysRevLett.61.169.
2. J. Dalibard and C. Cohen-Tannoudji, “Laser cooling below the Doppler limit by polarization gradients: simple theoretical models,” Journal of the Optical Society of America B 6, 2023–2045 (1989), doi:10.1364/JOSAB.6.002023.
3. Y. Castin, H. Wallis, and J. Dalibard, “Limit of Doppler cooling,” Journal of the Optical Society of America B 6, 2046–2057 (1989), doi:10.1364/JOSAB.6.002046.
4. P. J. Ungar, D. S. Weiss, E. Riis, and S. Chu, “Optical molasses and multilevel atoms: theory,” Journal of the Optical Society of America B 6, 2058–2071 (1989), doi:10.1364/JOSAB.6.002058.
5. A. Aspect, E. Arimondo, R. Kaiser, N. Vansteenkiste, and C. Cohen-Tannoudji, “Laser cooling below the one-photon recoil energy by velocity-selective coherent population trapping,” Physical Review Letters 61, 826–829 (1988), doi:10.1103/PhysRevLett.61.826.
6. M. Kasevich and S. Chu, “Laser cooling below a photon recoil with three-level atoms,” Physical Review Letters 69, 1741–1744 (1992), doi:10.1103/PhysRevLett.69.1741.
7. J. Reichel, F. Bardou, M. Ben Dahan, E. Peik, S. Rand, C. Salomon, and C. Cohen-Tannoudji, “Raman cooling of cesium below 3 nK,” Physical Review Letters 75, 4575–4578 (1995), doi:10.1103/PhysRevLett.75.4575.
8. C. Monroe, D. M. Meekhof, B. E. King, S. R. Jefferts, W. M. Itano, D. J. Wineland, and P. Gould, “Resolved-sideband Raman cooling of a bound atom to the 3D zero-point energy,” Physical Review Letters 75, 4011–4014 (1995), doi:10.1103/PhysRevLett.75.4011.
9. S. E. Hamann, D. L. Haycock, G. Klose, P. H. Pax, I. H. Deutsch, and P. S. Jessen, “Resolved-sideband Raman cooling to the ground state of an optical lattice,” Physical Review Letters 80, 4149–4152 (1998), doi:10.1103/PhysRevLett.80.4149.
10. A. T. Grier, I. Ferrier-Barbut, B. S. Rem, M. Delehaye, L. Khaykovich, F. Chevy, and C. Salomon, “$\Lambda$-enhanced sub-Doppler cooling of lithium atoms in D$_1$ gray molasses,” Physical Review A 87, 063411 (2013), doi:10.1103/PhysRevA.87.063411.
11. C. N. Cohen-Tannoudji, “Nobel Lecture: Manipulating atoms with photons,” Reviews of Modern Physics 70, 707–719 (1998), doi:10.1103/RevModPhys.70.707.
12. W. D. Phillips, “Nobel Lecture: Laser cooling and trapping of neutral atoms,” Reviews of Modern Physics 70, 721–741 (1998), doi:10.1103/RevModPhys.70.721.
13. H. J. Metcalf and P. van der Straten, Laser Cooling and Trapping (Springer, 1999), doi:10.1007/978-1-4612-1470-0.
14. C. Cohen-Tannoudji and D. Guéry-Odelin, Advances in Atomic Physics: An Overview (World Scientific, 2011), doi:10.1142/6631.

Exercises

1. Polarization without an intensity gradient

For

$$\mathbf E^{(+)}(z) = E_0 \left( \mathbf e_x e^{ikz} + \mathbf e_y e^{-ikz} \right),$$

show that the intensity is constant and derive the circular-component intensities. Identify positions of pure circular polarization.

Solution

Orthogonality eliminates the interference term in the scalar intensity:

$$\begin{aligned} |\mathbf E^{(+)}|^2 &= |E_0|^2 \left( |\mathbf e_x|^2 + |\mathbf e_y|^2 \right) \\ &= 2|E_0|^2. \end{aligned}$$

Using circular basis vectors, the normalized component intensities are, up to which handedness is called plus,

$$\frac{I_{\sigma^\pm}}{I} = \frac{1}{2} \left[ 1 \pm \sin(2kz) \right].$$

One component vanishes and the other carries all the intensity when $\sin(2kz)=\pm1$, namely

$$z = \frac{\lambda}{8} + n\frac{\lambda}{4}.$$

The handedness alternates between neighboring positions.

2. Energy removed in an ideal Sisyphus cycle

Model two shifted potentials as

$$U_\pm(z) = \frac{U_0}{2} \left[ 1 \pm \cos(2kz) \right].$$

An atom begins at a minimum of $U_-$, climbs to its maximum, and is pumped to $U_+$. Neglect recoil. How much mechanical energy is removed?

Solution

At a minimum,

$$U_-^{\min} = 0.$$

At the corresponding crest,

$$U_-^{\max} = U_0, \qquad U_+^{\min} = 0.$$

The atom converts $U_0$ of kinetic energy into potential energy while climbing. Optical pumping changes the internal dressed state at fixed position and releases the potential difference to the light field and spontaneously emitted photon. The ideal mechanical energy loss is therefore

$$\Delta E_{\mathrm{mech}} = -U_0.$$

Real cycles pump over a distribution of positions and add recoil, so the mean magnitude is smaller.

3. Doppler and recoil scales

Using $\lambda=780\,\mathrm{nm}$, $\Gamma/2\pi=6.07\,\mathrm{MHz}$, and $m=1.443\times10^{-25}\,\mathrm{kg}$, calculate $T_D$, $T_r$, and their ratio for $^{87}\mathrm{Rb}$.

Solution

The wave number is

$$k = \frac{2\pi}{\lambda} \simeq 8.06\times10^6\,\mathrm{m^{-1}}.$$

Then

$$T_D = \frac{\hbar\Gamma}{2k_{\mathrm B}} \simeq 1.46\times10^{-4}\,\mathrm K,$$

and

$$T_r = \frac{\hbar^2k^2}{2mk_{\mathrm B}} \simeq 1.81\times10^{-7}\,\mathrm K.$$

Therefore

$$\frac{T_D}{T_r} \simeq 8.0\times10^2.$$

The broad transition leaves almost three orders of magnitude between the two-level Doppler and one-photon recoil temperatures.

4. Verify the dark state

For the $\Lambda$ interaction Hamiltonian in the text, verify that

$$|D\rangle = \frac{ \Omega_2|g_1\rangle-\Omega_1|g_2\rangle }{ \sqrt{|\Omega_1|^2+|\Omega_2|^2} }$$

has no excited-state coupling.

Solution

Only the excitation terms contribute:

$$\begin{aligned} H_{\mathrm{int}}|D\rangle &= \frac{\hbar}{2\Omega} |e\rangle \left( \Omega_1\Omega_2 - \Omega_2\Omega_1 \right) \\ &= 0. \end{aligned}$$

The cancellation is coherent. Population loss, dephasing, or unequal two-photon phases can spoil it even when both one-photon couplings remain strong.

5. Dark-resonance velocity width

Counterpropagating Raman beams have nearly equal wave-number magnitude $k$, so $|\Delta\mathbf k|\simeq2k$. If the dark resonance has angular-frequency half-width $\gamma_D/2\pi=20\,\mathrm{kHz}$ at $\lambda=780\,\mathrm{nm}$, estimate the one-dimensional velocity half-width.

Solution

The velocity width is

$$\Delta v_D \sim \frac{\gamma_D}{2k}.$$

Using

$$\gamma_D = 2\pi(20\times10^3)\,\mathrm{s^{-1}}, \qquad k = \frac{2\pi}{780\times10^{-9}\,\mathrm m},$$

gives

$$\Delta v_D \simeq 7.8\times10^{-3}\,\mathrm{m\,s^{-1}}.$$

This is comparable to a single-photon recoil velocity for rubidium, so a continuous classical diffusion model is already becoming questionable.

6. Sideband cooling step

Show that the Raman transition

$$|\downarrow,n\rangle \rightarrow |\uparrow,n-1\rangle$$

followed by recoil-free repumping removes $\hbar\omega_t$ of motional energy. Why does the cycle stop at $n=0$?

Solution

The motional energy changes from

$$E_n = \hbar\omega_t \left( n+\frac{1}{2} \right)$$

to

$$E_{n-1} = \hbar\omega_t \left( n-\frac{1}{2} \right).$$

Their difference is

$$E_{n-1}-E_n = -\hbar\omega_t.$$

Ideal repumping changes only the internal state, so this motional reduction survives the reset. The red-sideband matrix element is proportional to $\sqrt n$ and vanishes at $n=0$; there is no state $|-1\rangle$.

7. Sideband thermometry

A weak probe measures

$$\frac{I_{\mathrm{red}}}{I_{\mathrm{blue}}} = 0.12.$$

Assuming a thermal mode, find $\bar n$ and $P_0$.

Solution

Let $r=0.12$. Since

$$r = \frac{\bar n}{\bar n+1},$$

the mean occupation is

$$\bar n = \frac{r}{1-r} = \frac{0.12}{0.88} \simeq 0.136.$$

The thermal ground-state probability is

$$P_0 = \frac{1}{1+\bar n} = 1-r = 0.88.$$

8. Diagnose a mechanism claim

An alkali cloud reaches $8\,\mu\mathrm K$, below its Doppler temperature. The temperature rises sharply with a residual magnetic field, scales nearly linearly with cooling intensity, and changes little when a two-photon Raman detuning is scanned. Which mechanism is best supported, and what additional test would you perform?

Solution

The magnetic sensitivity and approximate light-shift scaling support polarization-gradient Sisyphus cooling. The absence of a narrow two-photon-detuning feature argues against a $\Lambda$ dark-state mechanism as the dominant effect.

A discriminating test is to change the relative beam polarizations while holding total intensity and one-photon detuning fixed. Suppressing the polarization gradient should weaken the cooling. One should simultaneously measure atom number to rule out velocity-selective loss.