Codex @codex  0

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.

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].

An optical lattice is a periodic optical trap produced by interfering light fields. Its periodic potential gives atomic Bloch states and energy bands. For the displayed one-dimensional lattice, the recoil energy is $E_R={\hslash }^2{\pi }^2/(2m{d}^2)$ and $s$ measures the depth in recoil units.

Take the initial ground-state chemical atom at rest in the laboratory. The incident photon has energy $h\nu$ and momentum $h\nu /c$. The excited chemical atom has rest mass $M+\Delta M$ and must carry that momentum. Conservation of energy and conservation of momentum giveSquaring cancels the photon-energy square:ConsequentlyThe extra positive term is the photon absorption recoil correction and supplies atomic recoil energy. The formula is exact for the initially stationary isolated chemical atom; it does not assume $\Delta M/M\ll 1$. Initial and final chemical atoms cannot both remain at rest during single-photon absorption without another body taking momentum.

Laser cooling removes atomic kinetic energy through a velocity-dependent optical force, with scattered photons carrying away energy and entropy. For two counterpropagating beams detuned below resonance, motion towards one beam shifts that light closer to resonance, so the atom preferentially absorbs momentum opposite its velocity. At small velocity the net force is $F\simeq -{\gamma }_{v}v$ with ${\gamma }_v>0$. Three orthogonal beam pairs form optical molasses. Spontaneous emission has zero mean recoil in a symmetric distribution, but its momentum variance is nonzero: random absorption and emission cause heating. The cooling and heating balance sets a finite temperature.

For a weakly saturated two-level transition, write $\delta ={\omega }_L-{\omega }_0<0$ and let $\Gamma$ be the excited-state population decay rate. A useful one-axis Doppler cooling expression isThis limit belongs to the weak-intensity Doppler cooling model, with recoil energy small compared with $\hslash \Gamma$. Experimentally, released sodium molasses clouds were probed by time of flight; the unusually small expansion demonstrated temperatures below that limit. This distinguishes an observed failure of a two-level cooling model from a failure of conservation of energy.^[ref]

Multilevel atoms permit polarization-gradient cooling: spatially varying light polarization produces different light shifts of ground-state sublevels. In Sisyphus cooling, atoms climb a light-shift potential and are preferentially optically pumped near its top into a lower-energy sublevel, repeatedly losing kinetic energy. The optical pumping and emitted light supply dissipation. This can reach below the Doppler cooling limit; the recoil scalethen becomes important. The definition of a recoil temperature sometimes differs by a factor of two, so the energy definition is the unambiguous comparison. Recoil is not a universal lower bound for every cooling protocol: suitable dark-state or resolved sideband cooling schemes can suppress scattering near the desired motional state. Resolved sideband cooling preferentially excites a transition lowering the vibrational quantum number, followed by decay that on average preserves it in the Lamb-Dicke regime.^[ref]

Optical molasses supplies damping rather than a position-restoring force. A magneto-optical trap adds a magnetic-field gradient and appropriately chosen circular polarizations; the position-dependent Zeeman splitting gives a restoring radiation pressure imbalance as well as damping. In real alkali atoms, repumping light returns population that would otherwise accumulate in an uncoupled hyperfine level. Magnetic and optical dipole traps can subsequently hold the cold sample without continuous resonant scattering.

Applications include long-interrogation atomic clocks, precision spectroscopy, atom interferometry, and preparation for quantum simulation and Bose-Einstein condensation. Cooling reduces motion-induced frequency shifts and allows longer coherent interrogation; deeper quantum degeneracy commonly needs subsequent evaporative cooling of atoms. Finite capture velocity, available transitions, recoil heating, stray magnetic fields, and loss or heating from unwanted scattering constrain a real experiment.

For a free-space electric dipole transition, the spontaneous emission rate in SI units isHere $|\mathbf{d}{|}^2$ includes the relevant dipole components; a specified single polarization requires its corresponding matrix element.

A complete scattering cycle returns the internal state to the ground state and changes the atomic momentum to $\mathbf{P}+\mathbf{p}-{\mathbf{p}}^{\prime }$. ThusFor photon-scattering recoil heating with isotropic emission of fixed momentum magnitude, $⟨{\mathbf{p}}^{\prime }⟩=0$, givingThe first term is the velocity-dependent work of absorption and can cool an atom if the photon opposes its motion. The last two are positive absorption and emission recoil energies. If the actual absorbed directions are also isotropic and independent of $\mathbf{P}$, the first term averages to zero. For approximately equal incoming and outgoing photon momenta $\hslash k$,Isotropic incident illumination alone would not guarantee isotropic absorption if Doppler-dependent excitation selected directions.

Source qualification for the trapped calculation. Set $m=M$ for the same atom. The printed Gaussian $\exp [-x^2/(2l_0{)}^2]$ has the wrong width for the specified harmonic oscillator: the printed ground and first-excited functions have squared norms $\sqrt{2}$ and $2\sqrt{2}$, respectively. The normalized oscillator states fixed by the given $l_0=\sqrt{\hslash /(m{\omega }_{\mathrm{HO}})}$ instead useAlso, one common complex phase multiplying both off-diagonal terms does not make the Hamiltonian a Hermitian operator. Use an absorption matrix element $(\hslash \Omega /2)e^{-ikx}$ and its Hermitian conjugate for emission; choosing the opposite propagation convention reverses the imaginary sideband phase but not a transition probability. The physical results below use these explicitly stated corrections. Parts (i) and (ii) also give the literal integrals of the printed functions.

Conditional motional heating and recoil averaging. With $\eta =k{l}_0/\sqrt{2}$, an oscillator eigenstate has $⟨K{⟩}_n=⟨V{⟩}_n=(n+1/2)\hslash {\omega }_{\mathrm{HO}}/2$. Therefore the absorption outcomes haveThe same statements hold for the specified emission outcomes $|e,0⟩\to |g,0⟩$ and $|e,0⟩\to |g,1⟩$. Internal excitation energy goes into, or comes from, the photon and is distinct from this external energy change. Emission into the first motional state makes the emitted photon correspondingly less energetic.

A sudden recoil from $|0⟩$ gives probabilities $P_n=e^{-{\eta }^2}{\eta }^{2n}/n!$, so $P_1\simeq {\eta }^2\ll 1$ while $P_0\simeq 1$. The rare first sideband transfers a whole trap quantum, much larger than $E_R$; the mean is neverthelessThe emission average is the same for a one-dimensional recoil of the same magnitude. Projecting isotropic three-dimensional emission onto one axis instead gives a factor $⟨{\cos }^2\theta ⟩=1/3$ for that axis.

A sudden unobserved kick initially raises kinetic energy by $E_R$ and leaves the position distribution, hence potential energy, unchanged. It creates a coherent superposition of motional states. The equal kinetic/potential increments above concern energy-eigenstate-conditioned outcomes, or time-averaged/dephased motion; these are different measurement conditions. There is no contradiction.

With $\Gamma >{\omega }_{\mathrm{HO}}$, motional sidebands are not cleanly resolved, and a short-lived excited atom has little time to move before decay. In the stronger sudden limit $\Gamma \gg {\omega }_{\mathrm{HO}}$ the recoil-kick picture is direct. For only the stated inequality, finite linewidth weights can change numerical factors; it does not fix an exact continuous-wave scattering rate. The average heating remains of recoil order per absorption or emission, and of order $2E_R$ when both absorption and emission kicks are counted at full magnitude. For the axial component of isotropic three-dimensional emission, the emission contribution instead carries the factor $1/3$ noted above. Strong confinement suppresses the chance of changing $n$, but does not make its averaged recoil energy disappear. This agrees with free-space recoil heating while replacing small continuous energy changes by rare quantized motional changes.

For an ideal weakly driven two-level atom, opposing red-detuned beams give velocity-dependent absorption and a friction force near zero velocity. Random photon recoils also cause momentum diffusion. Their balance gives the minimum Doppler temperatureat detuning $\Delta =-\Gamma /2$ in the low-saturation model, where $\Gamma$ is the spontaneous population-decay rate, equivalently the natural full linewidth expressed as an angular frequency. This limit assumes the atom has no useful ground-state internal structure.

Real atoms possess magnetic sublevels. In a suitable multilevel transition, position-dependent polarization gives different optical couplings and different AC Stark energy shifts for those sublevels. Counterpropagating beams with orthogonal linear polarizations have nearly uniform total intensity but a periodically varying polarization. Thus spatially varying light shifts can exist even without an intensity standing wave.

In Sisyphus cooling, an atom on one ground-state branch climbs a light-shift hill and loses kinetic energy. Near the upper part of the hill, optical pumping tends to transfer population to a lower branch. spontaneous emission makes this transfer dissipative and randomizes the emitted photon direction. The atom can then climb another hill. Finite pumping time is essential: instantaneous local equilibration would remove the delayed internal-state response that produces friction. The original multilevel analysis relates these optical-pumping cycles to the spatially displaced light shifts.^[ref]

For large detuning compared with the natural linewidth, useful scales areup to angular-coupling and Rabi-frequency convention factors. Here $\Omega$ is the optical Rabi angular frequency. The slow-atom capture range is of order $kv\lesssim {\Gamma }_{\mathrm{op}}$. Typical energies in the Sisyphus regime track the light-shift scale rather than only $\hslash \Gamma$, so reducing intensity can lower temperature below the Doppler limit. Capture becomes narrower, however, and recoil heating, imperfect polarization, magnetic fields and scattered light eventually matter. Define recoil energy $E_r=(\hslash k{)}^2/(2m)$; this is an important scale, not a claim that every sub-Doppler method has the same universal final temperature.^[ref]

The experimental test must establish a narrow velocity distribution independently of fluorescence brightness. William D. Phillips's account describes moving a time-of-flight probe above and below the cloud, deliberately heating the atoms to check a previously absent signal, and release-and-recapture comparisons that included gravity and finite beam geometry. These tests separated very cold atoms from an insensitive detector or a faulty thermometer.^[ref] In a modern implementation, field cancellation and polarization checks test the mechanism, while independent expansion thermometry checks the inferred temperature. sub-Doppler cooling exploits internal-state dynamics omitted from the two-level Doppler model; it does not contradict that model within its assumptions.

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.

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]