🌡️ AMO Toolkit · Laser Cooling Simulator

Laser Cooling Simulator

From room temperature to quantum ground state. Explore Doppler cooling, sub-Doppler gray molasses, and resolved-sideband cooling: the three stages that bring atoms to the quantum limit.

Laser cooling exploits radiation pressure from near-resonant light to slow and cool atoms. Doppler cooling, the first stage, provides cooling whenever an atom moving toward a beam sees a red-shifted photon closer to resonance, leading to a velocity-dependent damping force with limit k_BT_D = ℏΓ/2. Sub-Doppler mechanisms (Sisyphus cooling, gray molasses) exploit light shifts, dark states, and optical pumping to cool below the Doppler limit; the single-photon recoil scale E_r/k_B = ℏ²k²/(2mk_B) is the natural lower energy scale. This simulator shows momentum-space dynamics in real time so you can see the transition from broad thermal distribution to cold peak.

Live Doppler Cooling Simulator

Drag the sliders — watch the velocity distribution narrow as atoms cool

Detuning δ/Γ −0.5

Saturation s₀ 0.5

T = — μK ▏= Doppler limit ← hot | velocity | hot →

Temperature Ladder

From room temperature to motional ground state, 9 orders of magnitude

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 the balance between velocity damping and momentum diffusion:

T_D = ℏΓ / (2k_B)

Force Curve F(v), interactive

Detuning δ/Γ −0.5

Saturation s₀ 0.5

α (damping) = —

T_Doppler vs Detuning

Minimum T_D = ℏΓ/(2k_B) at δ = −Γ/2

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 single-photon recoil scale E_r/k_B = ℏ²k²/(2mk_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.

Doppler Cooling Parameters by Species

⁶Li (671 nm)

MOT → cMOT → GM

Γ/2π = 5.87 MHz. T_D = 141 μK. Unresolved excited hyperfine structure, GM essential for sub-Doppler cooling. Single recoil E_r/k_B = 3.5 μK; one scatter cycle is roughly 7.1 μK without cooling.

⁸⁷Rb (780 nm)

MOT → PGC → evaporation

Γ/2π = 6.065 MHz. T_D = 146 μK. Well-resolved F levels. Sub-Doppler polarisation-gradient cooling (PGC) to ~4 μK.

¹³³Cs (852 nm)

MOT → cMOT → GM → RSB

Γ/2π = 5.234 MHz. T_D = 125 μK. Strong Sisyphus cooling. Tweezer arrays with RSB cooling to ⟨n⟩ < 0.1.

²³Na (589 nm)

Zeeman slower → MOT → PGC

Γ/2π = 9.795 MHz. T_D = 235 μK. Yellow D2 line. First atom laser (MIT 1997). High evaporation efficiency.

⁸⁸Sr (689 nm)

Blue MOT → Red MOT → magic tweezer

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)

Blue MOT → Green MOT → lattice

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)

The single-photon recoil energy E_r = ℏ²k²/(2m) sets the microscopic scale for random kicks. For Cs D2, E_r/k_B ≈ 99 nK. In practice, polarization-gradient cooling reaches μK-scale temperatures, often tens of recoil energies rather than a universal 2–4 recoil limit.

Λ-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)

Sublevel energies U±(x) and optical pumping at crests

Λ System, Gray Molasses Level Diagram

D1 Λ-system: dark state protected from scattering

⁶Li Gray Molasses: Experimental Result

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 + Γ_heat

Steady-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/4

EIT is especially useful for motional frequencies ω_trap < Γ (unresolved sideband regime), enabling ground-state cooling in shallower traps.

Sideband Spectrum, interactive

⟨n⟩ 5

η (Lamb-Dicke) 0.15

RSB/BSB asymmetry ≈ ⟨n⟩/(⟨n⟩+1), thermometry probe

⟨n⟩(t) Cooling Evolution

Initial ⟨n⟩ 20

Γ_cool/ω 0.05

⟨n⟩(t) = ⟨n⟩_ss + (n₀ − ⟨n⟩_ss)e^{−Γ_cool t}

All Cooling Methods at a Glance

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

Single recoil: E_r/k_B = ℏ²k²/(2mk_B). One absorption plus spontaneous emission gives a two-recoil heating scale ≈2E_r/k_B. For Cs D2: E_r/k_B ≈ 99 nK and 2E_r/k_B ≈ 198 nK.

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.

import pylcp
import numpy as np

# Define laser beams (3D molasses)
laser_beams = pylcp.laserBeams([
  {'kvec': np.array([1,0,0]), 's': 0.5, 'delta': -0.5},
  {'kvec': np.array([-1,0,0]), 's': 0.5, 'delta': -0.5},
  # ... ±y, ±z beams
], beam_type=pylcp.infinitePlaneWaveBeam)

# Define atom (e.g., Rb F=2 → F'=3)
atom = pylcp.hamiltonians.singleF(
    F=2, Fp=3, gamma=1, dipole_matrix_element=1
)

# Integrate force equation
eqn = pylcp.rateeq(laser_beams, magField, atom)
eqn.set_initial_pop(np.array([1/5]*5))
sol = eqn.equilibrium_populations(r=np.zeros(3), v=np.array([0.1,0,0]))

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.

from arc import Caesium

atom = Caesium()

# Get natural linewidth of D2 line
gamma = atom.getTransitionRate(6,0,0.5, 6,1,1.5)
print(f"Γ/2π = {gamma/2/np.pi/1e6:.3f} MHz")

# Rydberg C6 coefficient for 70S1/2
C6 = atom.getC6term(70,0,0.5, 70,0,0.5)
R_b = (abs(C6) / (hbar * 2*np.pi*1e6)) ** (1/6)
print(f"Blockade radius = {R_b*1e6:.1f} μm")

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, ρ})

import qutip as qt
N = 30  # Fock space truncation
a = qt.destroy(N)
H = eta * Omega * (a + a.dag())  # RSB drive
c_ops = [np.sqrt(gamma) * a]     # decay
result = qt.mesolve(H, rho0, tlist, c_ops,
                    e_ops=[a.dag()*a])

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