🧲 Tool 06 · MOT & Trap Designer

MOT & Trap Designer

Full parameter design tool for magneto-optical traps and Ioffe-Pritchard magnetic traps β€” from laser detuning and gradient to trap frequencies, depth, and evaporation efficiency.

A magneto-optical trap (MOT) is the first stage of every cold-atom experiment. It simultaneously cools and confines atoms via radiation pressure from three pairs of counter-propagating beams in a magnetic quadrupole field. Understanding the damping coefficient Ξ±, the spring constant ΞΊ, and the capture velocity v_c is essential for optimizing loading rates and determining whether subsequent sub-Doppler cooling (gray molasses, PGC) is needed before loading into tweezers or a lattice. Section 02 covers the Ioffe-Pritchard magnetic trap used for evaporative cooling to quantum degeneracy.
Model validity: the MOT calculator uses a 1D, low-velocity, two-level scattering-force model. Real alkali MOTs are multilevel systems with Clebsch-Gordan factors, optical pumping, polarization gradients, beam imbalance, reabsorption, and sub-Doppler physics. Treat Ξ±, ΞΊ, and capture velocity as scaling estimates, not calibrated predictions for a loaded cloud.

01 Magneto-Optical Trap (MOT)

Laser & Field Parameters
0.5Ξ“βˆ’Ξ“ (opt. Doppler)5Ξ“
0.11 (sat.)3
⁸⁷Rb D2: Ξ“/2Ο€ = 6.065 MHz Β· Ξ» = 780 nm Β· I_sat = 1.67 mW/cmΒ²
Damping coeff. Ξ±
β€”
10⁻²⁰ kg/s (1D, 2 beams)
F = βˆ’Ξ±v
Spring constant ΞΊ
β€”
10⁻¹⁷ N/m
F = βˆ’ΞΊr (radial)
MOT trap frequency
β€”
Hz
Ο‰ = √(ΞΊ/m)
Capture velocity
β€”
m/s
v_cap β‰ˆ |Ξ”|/k
Doppler temperature
β€”
ΞΌK
ℏΓ/2k_B (limit)
Max acceleration
β€”
m/sΒ²
a_max = ℏkΞ“s/(2m(1+s))

MOT force in 1D (two-beam, low velocity limit)

Each beam exerts a scattering force ℏk Γ— R_sc on the atom. Two counter-propagating beams with the quadrupole field create both a velocity-dependent damping force and a position-dependent restoring force. In the low-saturation, low-velocity limit:

$$F_{\rm MOT} \approx -\alpha v - \kappa r$$ $$\alpha = \frac{8\hbar k^2 s |\Delta|/\Gamma}{(1 + s + (2\Delta/\Gamma)^2)^2} \quad \text{[damping]}$$ $$\kappa = \alpha \cdot \frac{\mu_{\rm eff}}{\hbar k} \cdot \frac{dB}{dz} \quad \text{[spring via B-field gradient]}$$ $\mu_{\rm eff} = g_F m_F \mu_B$ (effective magnetic moment)
Live MOT Visualization  Β·  2D cross-section (horizontal + vertical beam pair)

02 Ioffe-Pritchard Magnetic Trap

IP Trap Field Parameters
10 G/cm200 G/cm500 G/cm
1 G20 G100 G
Radial freq Ο‰_r /2Ο€
β€”
Hz
√(ΞΌ_eff Bβ€²Β² / m Bβ‚€)
Axial freq Ο‰_z /2Ο€
β€”
Hz
√(ΞΌ_eff Bβ€³ / m)
Geometric mean freq
β€”
Hz = (Ο‰_rΒ² Ο‰_z)^(1/3) / 2Ο€
relevant for BEC T_c
Trap depth Uβ‚€
β€”
ΞΌK
ΞΌ_eff Γ— Bβ‚€ / k_B
Ξ· parameter
β€”
Evaporation efficiency: Ξ· = Uβ‚€ / (k_B T)
Ξ·=1 (poor)Ξ·=5 (good)Ξ·=10+ (excellent)
β€”

03 Trap Potential Profile

U(r) = ΞΌ_eff Β· |B(r)|, radial cut through IP trap
U(r), radial potential k_B T (atom energy) Trap depth Uβ‚€ = ΞΌ_eff Bβ‚€

04 Physics Reference

Radiation pressure and position-dependent detuning

In a MOT the two counter-propagating beams create opposing radiation pressure forces. For a moving atom at position r in the quadrupole field, the effective detuning for each beam is Ξ”_eff = Ξ” Β± kv Β± ΞΌ_eff B(r)/ℏ. This creates both velocity (Doppler) and position (Zeeman) dependent forces.

Common mistake, detuning sign conventions. This page uses \(\Delta=\omega_L-\omega_0\), so red detuning is negative. Some AMO notes use \(\Delta=\omega_0-\omega_L\), which flips every sign in the damping-force algebra. Always define the convention before interpreting "positive" and "negative" spring constants.

MOT restoring force

Sentence: red-detuned light plus Zeeman shifts makes the atom scatter more from the beam pushing it home.

F β‰ˆ -Ξ±v - ΞΊr

Failure mode: increasing gradient until optical pumping or capture velocity gets worse.

Capture is not temperature

Sentence: a larger capture velocity loads more atoms but does not guarantee colder atoms.

vc ∝ sqrt(aL)

Failure mode: optimizing fluorescence only and calling it temperature optimization.

Sub-Doppler needs details

Sentence: polarization gradients, dark states, and magnetic field zeroing decide the final cloud temperature.

Tfinal can be β‰ͺ TD

Failure mode: using Doppler temperature as a hard lower bound for alkalis.

$$R_{\rm sc} = \frac{\Gamma}{2} \cdot \frac{s}{1 + s + (2\Delta_{\rm eff}/\Gamma)^2}$$ Net force (two beams, $|kv| \ll |\Delta|$): $F = F_+ - F_- \approx -\alpha v - \kappa r$
Damping (1D, 2 beams): $$\alpha = \frac{8\hbar k^2 s|\Delta|/\Gamma}{(1 + s + (2\Delta/\Gamma)^2)^2}$$ Spring constant: $$\kappa = \frac{\mu_{eff}}{\hbar k}\cdot\frac{\partial B}{\partial r}\cdot\alpha = \mu_{eff}\cdot\frac{dB}{dz}\cdot\frac{8k^2 s|\Delta|/\Gamma}{(1+s+(2\Delta/\Gamma)^2)^2}$$ Optimal damping: $\Delta = -\Gamma/(2\sqrt{3}) \approx -0.29\Gamma \Rightarrow \alpha_{max}$

Non-zero minimum prevents Majorana losses

An IP trap combines a radial quadrupole field (wires or Ioffe bars: gradient Bβ€²) with axial pinch coils that create a non-zero minimum field Bβ‚€ and axial curvature Bβ€³. The non-zero field bottom eliminates Majorana spin-flip losses that plague quadrupole traps.

Field magnitude near trap center ($\rho \ll B_0/B'$): $$|B(\rho,z)| \approx B_0 + \left(\frac{B'^2}{2B_0} - \frac{B''}{2}\right)\rho^2 + \frac{B''}{2}z^2$$ Trap frequencies: $$\omega_r = \sqrt{\frac{\mu_{\rm eff}(B'^2/B_0 - B'')}{m}} \quad \text{[radial]}, \qquad \omega_z = \sqrt{\frac{\mu_{\rm eff} B''}{m}} \quad \text{[axial]}$$ Stability: $B'^2/B_0 > B''$ (radial trap not cancelled by curvature)
Trap depth: $U_0 \approx \mu_{eff} \cdot B_0$

QUIC trap (Esslinger group) and cloverleaf trap (Ketterle group) are two common IP-type implementations. Typical: Bβ‚€ = 1–10 G, Bβ€² = 100–300 G/cm, Bβ€³ = 50–200 G/cmΒ².

Ξ· = Uβ‚€ / (k_B T), truncation parameter

Evaporative cooling removes the hottest atoms (tail of the Maxwell-Boltzmann distribution) by RF or microwave induced spin-flip at a tunable cutoff energy Ξ· Γ— k_B T. As hot atoms leave and the cloud rethermalizes via elastic collisions, the temperature drops.

$$\frac{dT}{T} = -\frac{\eta - 4}{\eta - 3 + \ldots}\cdot\frac{dN}{N}$$ "Runaway evaporation": $\gamma_{\rm el}/\gamma_{\rm inel} > 150$ (Rb-87 3-body limit)
Good: $\eta \gtrsim 5$; optimal: $\eta \approx 8$–$10$ (BEC achievable); insufficient: $\eta < 4$
RF evaporation ramp: $\nu_{\rm RF}(t)$ decreases from $\nu_0$ to $\nu_f$ $$E_{cut} = h\nu_{RF} - m_F g_F \mu_B B_0$$

Spin-flip at the zero of the field

A pure quadrupole trap has a zero-field point at the center. Atoms that pass near this point can undergo a non-adiabatic Majorana spin-flip from a weak-field-seeking to a strong-field-seeking state, ejecting them from the trap. The rate increases as the cloud cools and atoms approach the zero-field region.

$$\Gamma_M \approx \omega_r \cdot \frac{\hbar}{m\omega_r\langle r^2\rangle} \approx \omega_r \cdot \frac{T_{\rm rec}}{T}$$ $$\omega_r = \sqrt{\frac{\mu_B B'^2}{m B_0}} \quad \text{[for quadrupole: } B_0 \to 0\text{]}$$ For quadrupole ($B_0 = 0$): $\Gamma_M \to \infty$, catastrophic losses
Fix: add bias field (IP trap, TOP trap) or optical plug (Rb BEC 1995).
Lifetime in quadrupole: $\tau_M \sim m\bar{v}/(\mu_B B')$ where $\bar{v} = \sqrt{2k_{\rm B}T/m}$

05 References