Why Cold Atoms for Quantum Computing?
Gate fidelity for a two-qubit Rydberg entangling gate scales as ℱ ≈ 1 − α⟨n⟩ − β/Ω_R where ⟨n⟩ is the mean vibrational quantum number and Ω_R is the Rabi frequency. Each cooling stage suppresses thermal motion, reducing state-preparation errors from >10% to <0.5%. Resolved-sideband cooling to ⟨n⟩ < 0.1 is now routine in tweezer arrays.
Radiation Pressure Force
A two-level atom moving with velocity v along a laser beam of detuning δ = ω_L − ω_0 experiences a Doppler-shifted resonance at ω_0 − kv. The scattering rate
Γ_sc = (Γ/2) · s₀ / [1 + s₀ + (2δ'/Γ)²]where δ' = δ + kv is the effective detuning and s₀ = I/I_sat. With counter-propagating beams the net force becomes velocity-dependent:
F(v) ≈ −αv, α = ℏk² · 4|δ|s₀Γ / (1 + s₀ + (2δ/Γ)²)²This viscous force (α > 0 for red detuning δ < 0) cools atoms. The minimum Doppler temperature is set by photon recoil:
T_D = ℏΓ / (2k_B)Force Curve F(v) — interactive
T_Doppler vs Detuning
Doppler Limit in Practice
For Cs D2 line (Γ/2π = 5.234 MHz): T_D = 125 μK. For ⁶Li D2 (Γ/2π = 5.87 MHz): T_D = 141 μK. Sub-Doppler effects can break this limit, reaching the photon recoil temperature T_r = ℏ²k²/(mk_B).
A 3D molasses (6 beams ±x, ±y, ±z) provides isotropic cooling. Adding a magnetic field gradient creates a magneto-optical trap (MOT) that both cools and traps atoms.
⁶Li (671 nm)
Γ/2π = 5.87 MHz. T_D = 141 μK. Unresolved excited hyperfine structure — GM essential for sub-Doppler cooling. T_r = 3.5 μK.
⁸⁷Rb (780 nm)
Γ/2π = 6.065 MHz. T_D = 146 μK. Well-resolved F levels. Sub-Doppler polarisation-gradient cooling (PGC) to ~4 μK.
¹³³Cs (852 nm)
Γ/2π = 5.234 MHz. T_D = 125 μK. Strong Sisyphus cooling. Tweezer arrays with RSB cooling to ⟨n⟩ < 0.1.
²³Na (589 nm)
Γ/2π = 9.795 MHz. T_D = 235 μK. Yellow D2 line. First atom laser (MIT 1997). High evaporation efficiency.
⁸⁸Sr (689 nm)
Narrow-line intercombination: Γ/2π = 7.5 kHz. T_D = 180 nK. Direct sub-μK MOT. Magic wavelength at 813 nm for optical lattice clock.
¹⁷⁴Yb (556 nm)
Narrow ¹S₀→³P₁: Γ/2π = 182 kHz. T_D = 4.4 μK. Nuclear spin I=0. Ideal for SU(N) fermions (¹⁷¹Yb) and optical clocks.
Sisyphus Cooling (Polarisation-Gradient)
In a standing wave with spatially varying polarisation (lin⊥lin), the light-shifted sublevel energies vary with position on the scale of λ/4. An atom in |m_F = +F⟩ moving in the +x direction climbs a potential hill, losing kinetic energy. At the top it optically pumps to |m_F = −F⟩ (lowest potential at that location), and the cycle repeats — Sisyphus loses kinetic energy each half-cycle.
T_Sisyphus ∝ U₀²/(E_r) ∝ I²/(δ² E_r)Temperature is limited by the photon recoil energy E_r = ℏ²k²/(2m). For Cs: T_r ≈ 200 nK. PGC can reach 2–4 T_r ≈ a few μK.
Λ-Enhanced Gray Molasses (GM) for ⁶Li
⁶Li has an unresolved excited-state hyperfine structure (Δ_HF < Γ), making Sisyphus PGC inefficient. Instead, a Λ-system gray molasses (D1 line at 671 nm) exploits a dark state between |F=1/2⟩ and |F=3/2⟩:
|dark⟩ = cos θ |F=1/2⟩ − sin θ |F=3/2⟩Dark-state atoms are velocity-selected: slow atoms stay dark (low scatter), fast atoms couple to the bright state and get cooled. GM achieves T ≈ 40–60 μK for ⁶Li, well below T_D = 141 μK.
Our experiment: GM on the D1 line at δ/Γ ≈ +2, reaching T ≈ 50 μK before loading into the optical tweezer array.
Sisyphus Potential (lin⊥lin)
Λ System — Gray Molasses Level Diagram
Thesis Result — ⁶Li Gray Molasses
By applying a D1 gray molasses (Ω/2π ≈ 3 Γ, δ/2π ≈ +2 Γ) for 3 ms after the MOT switch-off, we achieved T ≈ 50 μK from T_MOT ≈ 200 μK — a 4× improvement in temperature that directly reduces the motional-state occupation after loading into the tweezer: ⟨n⟩_initial ∝ T/ω_trap.
Harmonic Trap & Sideband Structure
An atom in a harmonic potential with trap frequency ω_trap/2π has quantised motional states |n⟩. The motional sidebands appear at:
ω_carrier ± n·ω_trap/2πIn the Lamb-Dicke regime (η² = E_r/ℏω ≪ 1), transition rates on the red sideband (RSB, ω − ω_trap) go as η²n, while the blue sideband (BSB) goes as η²(n+1).
η = k√(ℏ/2mω_trap) = √(E_r/ℏω_trap)Resolved Sideband Cooling Cycle
Drive the RSB: |g,n⟩ → |e,n−1⟩. Spontaneous emission returns mostly to |g,n−1⟩ (Lamb-Dicke suppression). Each cycle removes one vibrational quantum:
⟨ṅ⟩ = −Γ_cool · n + Γ_heatSteady-state: ⟨n⟩_ss = (Γ_heat/Γ_cool) ≈ η²Γ²/(4ω²) for weak sideband. In practice ⟨n⟩_ss < 0.1 is achievable.
EIT-Assisted Cooling
Electromagnetically Induced Transparency (EIT) cooling uses a probe + coupling beam to engineer a narrow absorption feature at the RSB frequency. The dark state suppresses carrier scattering while the RSB absorption is enhanced, achieving faster cooling rates than conventional RSB cooling.
⟨n⟩_EIT ≈ (γ_dark/ω_trap)² · 1/4EIT is especially useful for motional frequencies ω_trap < Γ (unresolved sideband regime), enabling ground-state cooling in shallower traps.
Sideband Spectrum — interactive
⟨n⟩(t) Cooling Evolution
| Method | Stage | T_min | Limit | Typical atoms | Notes |
|---|---|---|---|---|---|
| Zeeman Slower | Pre-cooling | ~10 mK | Doppler | Most alkalis, Sr, Yb | Slows hot atomic beam before MOT |
| 3D MOT | Stage 1 | ~100–200 μK | Doppler | All | Combines cooling + 3D spatial trapping |
| Compressed MOT (cMOT) | Stage 1b | ~50 μK | Doppler / PGC | Cs, Rb, Li | Ramp B-field + detuning before sub-Doppler |
| Polarisation-Gradient (PGC) | Stage 2 | 2–10 μK | Photon recoil | Rb, Cs, K, Na | lin⊥lin or σ+σ− standing waves, no B-field |
| D1 Gray Molasses (GM) | Stage 2 | ~40–60 μK | Photon recoil | ⁶Li, K (unresolved HF) | Λ-system dark state, works on D1 line |
| Narrow-line MOT (red) | Stage 2 | ~1 μK | Photon recoil | Sr (689 nm), Yb (556 nm) | Single-photon recoil kicks visible; direct sub-μK |
| Resolved Sideband (RSB) | Stage 3 | ⟨n⟩ < 0.1 | Lamb-Dicke, η | Cs, Rb, Ca⁺, Mg⁺ in tweezers/lattices | Requires resolved ω_trap > Γ; closes on RSB |
| EIT Cooling | Stage 3 | ⟨n⟩ < 0.1 | Dark-state linewidth | Ca⁺, Mg⁺, Sr, neutral atoms | Dark resonance engineered at RSB; faster rate |
Temperature Scale Comparison (log₁₀ T / K)
Doppler Limit
T_D = ℏΓ/(2k_B). Set by the balance between laser cooling force and random recoil kicks from spontaneous emission. Scales with linewidth.
Recoil Limit
T_r = ℏ²k²/(mk_B). Every scattered photon imparts ℏk momentum. For Cs: T_r ≈ 198 nK. Fundamental barrier for free-space laser cooling.
Lamb-Dicke Limit
In a trap: ⟨n⟩_ss → η²Γ²/(4ω²). Deep into Lamb-Dicke regime (η → 0) and resolved sidebands (ω ≫ Γ), ⟨n⟩ → 0.
pylcp — Laser Cooling Physics
pylcp (Python Laser Cooling Physics) is an open-source package for
simulating laser-atom interactions. It supports rate-equation, optical Bloch equation,
and Hund's case (a) Hamiltonian approaches.
arc — Alkali Rydberg Calculator
arc provides atomic structure data, dipole matrix elements, and
Rydberg state properties for alkali atoms. Useful for computing polarizabilities,
C₆ coefficients, and Rydberg blockade radii.
QuTiP — Quantum Toolbox
For master equation simulation of the density matrix ρ under sideband cooling (Lindblad form):
ρ̇ = −i[H_eff, ρ] + Σ_k (L_k ρ L_k† − ½{L_k†L_k, ρ})Additional Libraries
AtomicUnits.jl (Julia) — unit conversions for atomic physics. QuantumOptics.jl — fast master-equation solvers in Julia. MOLSCAT — molecular scattering for Feshbach resonances. COMSOL — FEM for magnetic trap geometry and electrode design.
- [1] T. W. Hänsch and A. L. Schawlow, "Cooling of gases by laser radiation," Opt. Commun. 13, 68 (1975).
- [2] D. J. Wineland and H. Dehmelt, "Proposed 10¹⁴ δν < ν laser fluorescence spectroscopy on Tl⁺ mono-ion oscillator," Bull. Am. Phys. Soc. 20, 637 (1975).
- [3] S. Chu et al., "Three-dimensional viscous confinement and cooling of atoms by resonance radiation pressure," Phys. Rev. Lett. 55, 48 (1985). Nobel Prize 1997.
- [4] P. D. Lett et al., "Observation of atoms laser cooled below the Doppler limit," Phys. Rev. Lett. 61, 169 (1988).
- [5] J. Dalibard and C. Cohen-Tannoudji, "Laser cooling below the Doppler limit by polarization gradients: simple theoretical models," J. Opt. Soc. Am. B 6, 2023 (1989).
- [6] A. Aspect et al., "Laser cooling below single-photon recoil by velocity-selective coherent population trapping," Phys. Rev. Lett. 61, 826 (1988).
- [7] F. Sievers et al., "Simultaneous sub-Doppler laser cooling of fermionic ⁶Li and ⁴⁰K on the D1 line," Phys. Rev. A 91, 023426 (2015).
- [8] I. Diedrich et al., "Laser cooling to the zero-point energy of motion," Phys. Rev. Lett. 62, 403 (1989). First sideband cooling.
- [9] A. M. Kaufman et al., "Cooling a single atom in an optical tweezer to its quantum ground state," Phys. Rev. X 2, 041014 (2012).
- [10] J. D. Thompson et al., "Coherence and Raman sideband cooling of a single atom in an optical tweezer," Phys. Rev. Lett. 110, 133001 (2013).
- [11] A. J. Park et al., "Cavity-enhanced optical lattice clock," Phys. Rev. X 10, 021078 (2020). EIT cooling in neutral atoms.
- [12] B. J. Lester et al., "Rapid production of uniformly filled atomic arrays," Phys. Rev. Lett. 115, 073003 (2015).