Laser cooling

Laser cooling uses laser light to reduce the random velocities of chemical atoms, and therefore their translational temperature. The striking ingredient is selective photon absorption: an atom moving towards a beam sees a Doppler shift that can make it absorb more strongly. Each absorbed photon transfers momentum along the beam; repeated absorption and spontaneous emission can then remove kinetic energy. Doppler cooling makes this selection with the atom's motion, while sub-Doppler cooling also exploits its internal quantum states.

Cooling makes chemical atoms move slowly enough for longer interrogation in atomic clocks and atom interferometry, and easier loading into an optical trap. The same radiation pressure that pushes an atom can cool it: with suitably tuned opposing beams, whichever direction the atom moves, it receives the stronger push against that motion.

Consider equal opposing laser beams along the $x$ axis and an atom with velocity $v$. Let ${\omega }_L$ be their angular frequency, ${\omega }_0$ the atomic resonance angular frequency, $k=2\pi /\lambda >0$ their wavenumber, and $\delta ={\omega }_L-{\omega }_0$ the laser detuning. Red detuning means $\delta <0$. To first order in $v/c$, the atom seesFor small positive $v$, the beam travelling towards negative $x$ is closer to resonance and causes more photon absorption; its momentum kicks point against $v$. The beam travelling towards positive $x$ is further from resonance. Reversing $v$ exchanges the roles of the two beams. This is a Doppler effect on the absorption rates, rather than a need to follow each atom with a moving beam.

An absorbed photon gives the atom momentum $\pm \hslash k$, where $\hslash =h/(2\pi )$ is the reduced Planck constant. Subsequent spontaneous emission gives a recoil opposite to the emitted photon. For an emission pattern symmetric under reversal of direction, that recoil has zero expected value, but nonzero variance. Its mean therefore leaves the directional absorption force intact, while its randomness adds recoil energy. Random absorption times and choices of beam add further momentum fluctuations. Zero mean recoil does not mean zero heating.

Use a closed two-level system with spontaneous emission population decay $P_e(t)=P_e(0)e^{-\Gamma t}$. Here $\Gamma$ is a decay rate in inverse seconds and also the full width at half maximum of the natural linewidth expressed in angular frequency. The corresponding full width in hertz is $\Gamma /(2\pi )$; the angular half width is $\Gamma /2$. These conventions determine every factor of two and $2\pi$ below.

Let $s$ be the on-resonance saturation parameter of each beam, with total excitation weak enough that saturation and coherent coupling between beams can be neglected. The independent-beam scattering rates and net force areThe force is the photon momentum per absorption multiplied by the difference in absorption rates. With dimensionless laser detuning $d=\delta /\Gamma$ and dimensionless velocity $u=kv/\Gamma$, writeFor the example $d=-0.5$ and $u=0.25$, the positive and negative beam factors are respectively $4/13$ and $4/5$. Thus $F(0.25)=-32/65\simeq -0.4923$: the atom is moving right and the net force points left. All values in this example are dimensionless.

Subtracting the two fractions gives an especially useful sign check:For red detuning, the force opposes nonzero $v$ throughout this ideal model; its useful magnitude eventually falls at large $|u|$. Near zero, a Taylor series gives dampingAt $d=-0.5$, $F^{\prime }(0)=-2$ and $\alpha =\hslash k^{2}s$. For blue detuning the sign reverses and the small-velocity motion is amplified; at exact resonance the opposing mean forces cancel in this symmetric model. Its derivation assumes weak excitation, a closed two-level system, equal beam intensities and nonrelativistic velocities; real multilevel atoms need additional physics.

Damping removes random kinetic energy, but the discrete photon kicks cause diffusion in momentum. Their balance gives a nonzero temperature. Define the momentum diffusion coefficient $D_p$ by the noise contribution $d⟨p_x^2⟩/dt=2D_p$. With atom mass $m$ and linear damping coefficient $\alpha$,at stationary balance, where $k_B$ is the Boltzmann constant and mean momentum is zero.

For ideal three-dimensional optical molasses with three identical independent beam pairs and isotropic spontaneous emission, the $x$-pair absorption rate at rest is $R_x=\Gamma s/(1+4d^2)$. Absorption from that pair contributes $(\hslash k{)}^2{R}_x$ to the noise growth of $⟨p_x^2⟩$. Emission from all three pairs contributes the same amount: its total rate is $3R_x$ and each emitted photon has mean squared $x$-momentum $(\hslash k{)}^2/3$. Hence $D_p=(\hslash k{)}^2{R}_x$. Using the Doppler cooling force coefficient givesPut $a=-d>0$. Then $(1+4d^2)/(-d)=a^{-1}+4a\ge 4$, with equality at $a=1/2$. The ideal weak-field Doppler cooling limit is therefore $T_D=\hslash \Gamma /(2k_B)$ at $\delta =-\Gamma /2$. It is the minimum of this specified diffusion-and-damping model, rather than a universal lower bound for laser cooling. The diffusion approximation also assumes recoil energy much smaller than $\hslash \Gamma$. In a one-pair geometry, emission patterns and transverse recoil change the diffusion bookkeeping. ^[ref].

Optical molasses surrounds slow chemical atoms with opposing red-detuned laser beams, usually three pairs along perpendicular axes. In the elementary Doppler cooling picture, each pair damps one component of velocity. The result resembles viscous drag, even though the atoms are moving through a vacuum rather than a liquid.

Uniform balanced optical molasses supplies velocity damping, without a restoring force towards a chosen position. An atom at rest has zero mean force wherever it sits within the uniform illuminated region. Without noise or external forces, $m\dot{v}=-\alpha v$ givesThe atom approaches a position determined by its initial condition, not the beam intersection. Recoil-driven diffusion and gravity eventually remove atoms from finite beams. A long residence time in optical molasses is therefore compatible with the absence of position trapping. An optical trap provides spatial confinement by a different force mechanism; a magneto-optical trap adds a position-dependent magnetic field and suitable beam electromagnetic polarization to create a restoring force.

Steven Chu's Bell Laboratories team prepared sodium atoms by pulsed-laser evaporation of a pellet. A frequency-swept opposing beam slowed them before six crossed laser beams supplied optical molasses. The vacuum chamber used a shield coated with graphite and cooled with liquid nitrogen as a pump. Fluorescence provided the readout; the initial beams were aligned as closely counterpropagating as possible. ^[ref].

After equilibration, the team switched off the cooling light for a variable interval, then restored it. Faster atoms escaped during the dark interval; slower atoms were recaptured. The recaptured fraction versus delay constrained the velocity distribution and hence temperature. The inference depended on the initial cloud shape and recapture volume: Steven Chu later acknowledged adjusting a geometry correction towards the expected Doppler cooling limit. Delay was the scanned variable, and the cooling beams themselves supplied the release/recapture control; visible fluorescence alone could not establish temperature. ^[ref].

Real atoms have multiple internal quantum states. Optical pumping and spatially varying electromagnetic polarization can remove kinetic energy by additional mechanisms, producing sub-Doppler cooling below the ideal Doppler cooling limit. In polarization-gradient cooling, state-dependent energy shifts from the Stark effect and delayed optical pumping make atoms repeatedly climb potential energy hills. This mechanism needs the internal-state structure omitted by the independent-beam two-level system above. ^[ref].

William D. Phillips's team tested unexpectedly cold sodium molasses with time-of-flight atom detection. A probe above the cloud initially saw no atoms; deliberately heating them verified that the detector worked. Moving the probe below the cloud produced arrival signals consistent with temperatures below the Doppler cooling limit. Independent release/recapture checks used sharply apertured beams and included gravity; uncertain cloud and probe geometry still limited the inferred temperature. These controls tested the thermometer as well as the cooling theory. ^[ref].

The 1988 experiment loaded sodium optical molasses from a slowed atomic beam. After stopping loading, a shutter extinguished the six cooling beams. Released atoms crossed a lower probe; a photomultiplier recorded fluorescence versus arrival time. Calculations integrated the measured initial cloud distribution and an assumed three-dimensional Maxwell-Boltzmann distribution, including gravity and detection geometry. The minimum was $43\pm 20\mu \mathrm{K}$, versus the $240\mu \mathrm{K}$ Doppler cooling limit.^[ref]

Probe separation, probe intensity, cloud density and shutter timing were varied. Fountain-height, horizontal-spreading and release/recapture measurements also gave below-limit results. Geometry remained an uncertainty; agreement among complementary thermometers made a detector artifact less plausible. These measured released atoms, rather than an in-situ velocity distribution.^[ref]

Steven Chu's account connects the explanation to optical pumping, Stark effect shifts and spatially changing electromagnetic polarization; Claude Cohen-Tannoudji's lecture separates this mechanism from subrecoil schemes. Their lectures and William D. Phillips's describe complementary experimental and theoretical routes.^[ref][ref]

A spatial image after a chosen dark interval is a different readout from continuously recording arrivals at a fixed probe. Here is an original idealised fluorescence imaging model: release a dilute cloud, let it move freely, then apply a short uniform probe pulse. Neglect collisions, residual forces, attenuation and motion during the pulse. Use each image's cloud centre to remove mean drift. Uniform gravity shifts that centre without changing the centred variance.

For each atom let $x(t)=x_0+v_{0}t$ along a horizontal imaging axis. Expanding the square about the expected values gives exactlyNo Gaussian distribution is needed for this variance identity. A jointly Gaussian distribution of initial positions and velocities does make the later density Gaussian. For a separable three-dimensional cloud viewed along $z$, integration along the viewing direction gives the column densityHere $N$ is atom number. Integrating this image along $y$ leaves a one-dimensional Gaussian distribution with standard deviation ${\sigma }_x(t)$. This is a projected density measurement, not a three-dimensional radius or a direct velocity histogram. Under the uniform, optically thin probe assumptions, expected fluorescence is proportional to this column density. A known, independent point spread function with finite variance adds its own variance to the measured width; probe duration and intensity gradients need separate treatment.

For a Maxwell-Boltzmann distribution, $T_x=m{\sigma }_v^2/k_B$, where $m$ is atomic mass and $k_B$ the Boltzmann constant. For a nonthermal distribution this defines a directional second-moment temperature, which need not describe every feature of the distribution. A covariance between position and velocity can initially narrow the cloud even though ${\sigma }_v^2$ is positive. The physical covariance matrix requires $C_{xv}^2\le {\sigma }_{x0}^2{\sigma }_v^2$.

Choose ${\sigma }_{x0}=0.60\mathrm{mm}$, ${\sigma }_v=0.060\mathrm{m}{\mathrm{s}}^{-1}$ and $C_{xv}=-0.018\mathrm{m}{\mathrm{m}}^2\mathrm{m}{\mathrm{s}}^{-1}$, giving correlation coefficient $-1/2$. With $t$ in milliseconds, the numerical variance in square millimetres is $0.36-0.036t+0.0036t^2$. Its minimum occurs at $5\mathrm{ms}$, with standard deviation $\sqrt{0.27}\mathrm{mm}$. These chosen values satisfy the covariance matrix constraint; early contraction is compatible with free motion.

Generate thirteen synthetic variance observations at equally spaced times from $0$ to $30\mathrm{ms}$, adding independent Gaussian noise of standard deviation $0.012\mathrm{m}{\mathrm{m}}^2$. Linear least squares fits all three coefficients in $a+bt+c{t}^2$. The example yields $a\simeq 0.37026\mathrm{m}{\mathrm{m}}^2$, $b\simeq -0.036413\mathrm{m}{\mathrm{m}}^2\mathrm{m}{\mathrm{s}}^{-1}$ and $c\simeq 0.0035820\mathrm{m}{\mathrm{m}}^2\mathrm{m}{\mathrm{s}}^{-2}$. With the chosen approximate sodium mass $m=3.82\times 10^{-26}\mathrm{kg}$, the true $T_x$ is $9.96\mu \mathrm{K}$ and the fitted value is $9.91\mu \mathrm{K}$. This is a fit to generated width data, not historical observations.

If the same data are fitted to $a+c{t}^2$, forcing $C_{xv}=0$, the inferred $T_x$ falls to $6.79\mu \mathrm{K}$. The omitted term is partly absorbed into both fitted coefficients. In noiseless linear least squares with a free intercept, its precise effect on the quadratic coefficient iswhere the subscript denotes equal-weight moments over the chosen observation times, using one consistent unit system. Since the ratio is positive for distinct nonnegative times, negative $C_{xv}$ biases this estimate downwards. An incorrect release-time origin can also generate a linear term, so fitting it does not by itself identify a physical initial covariance. Timing calibration, spatial resolution, residual magnetic fields and probe-induced recoil energy remain separate controls. Neither a cold-looking image nor a single long expansion time identifies all three coefficients.

Polarization-gradient cooling uses spatial variation in electromagnetic polarization together with optical pumping between internal quantum states to cool moving chemical atoms.

Jean Dalibard and Claude Cohen-Tannoudji distinguish cooling by changing elliptical polarization from cooling by rotating linear polarization. The hill picture below concerns the former.^[ref]

Finite optical pumping time lets an atom climb a light-shift hill before transfer into a lower magnetic sublevel. Spontaneous emission dissipates the gained potential energy; climbing alone is conservative.^[ref]

To make the geometry and signs explicit, take an original choice of phase origin for equal opposing waves:The real electric field is $\mathrm{Re}[\mathcal{E}(z)e^{-i{\omega }_{L}t}]$; $\lambda$ is wavelength. The fractions $w_{\pm }$ are squared projections onto the defined circular basis vectors, divided by total intensity. Their sum is one: interference changes electromagnetic polarization despite constant total intensity. At $z=0,\lambda /4,\lambda /2$ the field is respectively ${\mathbf{ϵ}}_{+},{\mathbf{ϵ}}_{-},{\mathbf{ϵ}}_{+}$; halfway between these positions it is linearly polarized. Specifying the vectors avoids an observer-dependent verbal handedness convention.

For an illustrative $J_g=1/2\to J_e=3/2$ angular momentum transition, use the squared Clebsch-Gordan coefficients $1$ and $1/3$ for the strong and weak circular couplings from each lower magnetic sublevel. Interchanging circular components interchanges these strengths. The weighted sums are $w_{+}+w_{-}/3=(2+\cos (2kz))/3$ and $w_{-}+w_{+}/3=(2-\cos (2kz))/3$. Rescaling their common negative shift givesBoth potential energies are negative relative to the unilluminated level. Each hill has height $U_0$, and opposite hills are displaced by $\lambda /4$. An atom starting at the $U_{+}$ minimum at $z=0$ loses $U_0$ of kinetic energy if it reaches $z=\lambda /4$ on that branch. A pump transfer there lowers its potential energy by $U_0$, leaving it at a $U_{-}$ minimum. The next uphill segment ends at $z=\lambda /2$, where the favoured transfer is reversed. Actual transfers are stochastic and need not occur exactly at these extrema.

For one completed climb and pump, conservation of energy gives $\Delta K=-U_0$, while the ground-state potential energy returns to its starting value. Neglecting recoil and treating the scattered photon as the sole outgoing radiation change, $\hslash ({\omega }_{\mathrm{out}}-{\omega }_L)=U_0$. A potential drop is an internal-state reset, not an automatic equal upward kick in kinetic energy. Random emission nevertheless adds momentum fluctuations and recoil energy, so useful cooling requires the dissipated energy to compete with that heating.

The recoil energy $E_R=(\hslash k{)}^2/(2m)$ sets a different scale from $k_B{T}_D=\hslash \Gamma /2$. Conventional Sisyphus cooling reaches the recoil regime when reducing the hill height also weakens cooling; it is not a claim of arbitrarily small temperature at vanishing intensity. Subrecoil techniques need further physics.^[ref]

Consider an original reduced rate equation model using the same potentials. Write $s(z)=\cos (2kz)$ and stipulate positive optical pumping rates $R_{+\to -}={\gamma }_p(1-s)$, $R_{-\to +}={\gamma }_p(1+s)$, with ${\gamma }_p>0$. They pump towards $+$ at $z=0$ and towards $-$ at $z=\lambda /4$, consistently with the diagram. These rates retain the relevant symmetry and delay, while omitting excited-state dynamics, quantum coherence and recoil diffusion.

For populations $p_{+}+p_{-}=1$, put $q=p_{+}-p_{-}$. The rate equations conserve population and give $\dot{q}=2{\gamma }_p(s-q)$. The stationary atom has $q=s$; a moving atom cannot follow it instantaneously. On a prescribed trajectory $z=vt$, the periodic solution isSubstitution into $v{q}^{\prime }=2{\gamma }_p(s-q)$ verifies both coefficients. The mean mechanical force is $F=-p_{+}{U}_{+}^{\prime }-p_{-}{U}_{-}^{\prime }=-U_{0}kq\sin (2kz)$. Averaging over a spatial period, with mean squared sine $1/2$ and mean sine times cosine zero, givesThus the force opposes either sign of velocity and vanishes at rest. A small pumping rate makes the small-velocity friction slope large while narrowing its useful velocity range. Holding $v$ fixed while sending ${\gamma }_p$ to zero makes the force vanish; the apparently divergent zero-velocity slope is not unlimited cooling. Because this model has no momentum diffusion, it predicts no equilibrium temperature. It demonstrates the lag mechanism rather than replacing a multilevel quantum calculation.

A magneto-optical trap combines laser cooling with a position-dependent magnetic field and suitable beam electromagnetic polarization to produce a restoring force. Its spatial confinement is additional to the velocity damping of uniform optical molasses.

A device that slows a beam of chemical atoms with laser electromagnetic radiation, while a spatially varying magnetic field compensates the changing Doppler effect.

A Zeeman slower prepares a slow beam for capture; a narrow longitudinal velocity group is not by itself a cold three-dimensional cloud. The historical apparatus measured slowed-beam velocities with a separate fluorescence probe scanned in frequency, just after switching off the slowing beam. This avoided simultaneous-beam complications.^[ref]

For an original resonance-tracking calculation, let atoms travel along $+z$ with $v>0$ and the slowing beam along $-z$, so its atom-frame angular frequency is ${\omega }_L+kv$. Approximate the selected cycling transition by ${\omega }_{eg}(B)={\omega }_0+\zeta B$ on one field branch; $\zeta$ is its signed angular frequency shift per unit magnetic field. A multilevel atom may require a nonlinear ${\omega }_{eg}(B)$ outside that branch. Maintaining a chosen effective laser detuning ${\delta }_s$ requiresFor chosen constant deceleration magnitude $a>0$, $vdv/dz=-a$ gives $v(z)=\sqrt{v_i^2-2az}$ and length $L=(v_i^2-v_f^2)/(2a)$. Substitution supplies a square-root magnetic field profile, plus a constant offset. The signs follow from the specified directions; reversing the selected transition's magnetic slope reverses the required field variation.

As a hypothetical design, take $v_i=300\mathrm{m}{\mathrm{s}}^{-1}$, $v_f=30\mathrm{m}{\mathrm{s}}^{-1}$ and $a=5.0\times 10^4\mathrm{m}{\mathrm{s}}^{-2}$. The required length is $0.891\mathrm{m}$ and flight time $(v_i-v_f)/a=5.4\mathrm{ms}$. For chosen $\lambda =589\mathrm{nm}$ and $m=3.82\times 10^{-26}\mathrm{kg}$, photon recoil velocity $\hslash k/m\simeq 0.0294\mathrm{m}{\mathrm{s}}^{-1}$ implies roughly $9.17\times 10^3$ directed absorption events. The mean required scattering rate is about $1.70\times 10^6{\mathrm{s}}^{-1}$. This budget must lie below the available scattering rate, with margin for capture and intensity variation. Optical pumping into unwanted hyperfine structure or magnetic sublevels can interrupt the cycle. Longitudinal slowing, transverse recoil heating and final cloud temperature are separate quantities.

Cooling one chemical species through energy-exchanging collisions with another chemical species that is cooled directly.