Optical Dipole Traps

An optical dipole trap confines a polarizable particle through the position-dependent AC Stark shift produced by an inhomogeneous light field. For an atom in internal state $|a\rangle$, the center-of-mass Hamiltonian is approximately

$$H_{\mathrm{cm},a} = \frac{\mathbf p^2}{2m} + U_a(\mathbf r),$$

where $U_a(\mathbf r)$ is the light-shift energy of the adiabatically followed internal state. This leading model neglects geometric scalar/vector potentials from spatially changing internal eigenvectors; slow internal following alone does not always make those terms negligible. The mechanical force is

$$\mathbf F_a(\mathbf r) = -\boldsymbol{\nabla}U_a(\mathbf r).$$

An ideal dipole trap confines without resonant photon scattering and provides no intrinsic damping. Its time-independent force is conservative. Real dipole traps are only approximately conservative: off-resonant scattering, technical noise, background collisions, internal-state changes, and finite trap depth all produce heating or loss.

Depth, curvature, and coherence are distinct: a deep trap can scatter too strongly for a coherent experiment, and a tight local minimum can have a small escape barrier.

Required background. AC Stark shifts explain the internal-state energy change; dynamic polarizability provides the multilevel response; magneto-optical traps describe the dissipative preparation commonly used before loading.

The light-induced potential

Write a monochromatic electric field as

$$\mathbf E(\mathbf r,t) = \operatorname{Re} \left[ E_0(\mathbf r) \boldsymbol{\epsilon}(\mathbf r) e^{-i\omega_Lt} \right],$$

where $E_0$ is the peak field amplitude and $\boldsymbol{\epsilon}^{\,*}\mathbin{\cdot}\boldsymbol{\epsilon}=1$. The cycle-averaged intensity in vacuum is

$$I(\mathbf r) = \frac{1}{2} \epsilon_0c |E_0(\mathbf r)|^2.$$

For a passive atom, define the effective polarizability in state $|a\rangle$ by

$$\alpha_a^{\mathrm{eff}} = \boldsymbol{\epsilon}^{\,*} \mathbin{\cdot} \boldsymbol{\alpha}_a(\omega_L) \mathbin{\cdot} \boldsymbol{\epsilon}.$$

This quantity may depend on position if the polarization, quantization axis, or internal eigenstate varies.

For a two-level comparison, use atom-minus-laser detuning

$$\Delta = \omega_0-\omega_L.$$

Thus:

• red detuning has $\Delta>0$;
• blue detuning has $\Delta<0$.

This sign convention is opposite to the commonly used $\Delta_L=\omega_L-\omega_0$. A formula imported from another source must be translated before its sign is interpreted.

Hierarchy of approximations

The simplest dipole-potential model assumes:

$$|\Delta| \gg \Gamma, \qquad |\Delta| \gg |\Omega|, \qquad s_{\mathrm{eff}} \ll 1,$$

where $\Gamma$ is the natural linewidth, $\Omega$ is the local resonant Rabi frequency, and $s_{\mathrm{eff}}$ is the detuning-reduced saturation parameter. These conditions suppress real excited-state population and power broadening.

Additional checks are separate:

• Atomic structure: a two-level reduction is unreliable when several fine-structure, hyperfine, molecular, or continuum channels contribute.
• Adiabatic internal following: motion through polarization or intensity gradients must not drive transitions between local dressed states.
• Slow envelope: the optical cycle is much faster than the center-of-mass motion, allowing a cycle-averaged potential.
• Coherent beams: fields at the same frequency and polarization may interfere. Their electric fields must be added before converting to intensity.
• Incoherent or rapidly dephased beams: their time-averaged potentials can be added for the mean potential when the interference averages much faster than motion. Beating near motional resonances or noisy interference can still heat, and two beams can drive resonant internal Raman transitions.

The term far detuned is therefore relative to all nearby resonances that matter, not only to one selected line.

From polarizability to mechanical force

With the field convention above, the conservative light shift is

$$U_a(\mathbf r) = - \frac{1}{4} \operatorname{Re} \alpha_a^{\mathrm{eff}}(\omega_L) |E_0(\mathbf r)|^2.$$

Equivalently,

$$U_a(\mathbf r) = - \frac{ \operatorname{Re} \alpha_a^{\mathrm{eff}}(\omega_L) }{ 2\epsilon_0c } I(\mathbf r).$$

For spatially uniform polarizability,

$$\mathbf F_a = \frac{ \operatorname{Re}\alpha_a^{\mathrm{eff}} }{ 2\epsilon_0c } \boldsymbol{\nabla}I.$$

A positive real polarizability pulls the atom toward higher intensity. A negative real polarizability pushes it toward lower intensity. If the local polarization changes, the force must instead be calculated from the gradient of the complete eigenenergy $U_a(\mathbf r)$; differentiating intensity alone can miss vector and tensor light-shift forces.

Far-detuned two-level form

For an isolated ground-state transition in the rotating-wave and weak-excitation limits ($\Gamma,|\Omega|\ll|\Delta|\ll\omega_0$),

$$U_g(\mathbf r) \simeq - \frac{ \hbar|\Omega(\mathbf r)|^2 }{ 4\Delta }.$$

The excited-state admixture is of order

$$P_e \simeq \frac{ |\Omega|^2 }{ 4\Delta^2 }.$$

Consequently, the conservative force survives even when the instantaneous excited-state population is very small. The atom is not being pulled because it repeatedly absorbs photons from the high-intensity side. The force is the gradient of a dressed-state energy.

For an alkali atom, a single detuning is often inadequate. The D1 and D2 lines can contribute with different weights, and the vector and tensor parts depend on hyperfine state, polarization, and magnetic field. Near a cancellation between contributions, the potential may be small while the scattering rate remains appreciable. The full dynamic-polarizability calculation is then mandatory.

Red and Blue Detuning

In the ideal two-level ground-state model,

$$\begin{array}{c|c|c} \text{detuning} & U_g & \text{preferred region} \\ \hline \Delta>0\ \text{(red)} & U_g<0 & \text{high intensity} \\ \Delta<0\ \text{(blue)} & U_g>0 & \text{low intensity} \end{array}$$

Red-detuned bright traps

A focused red-detuned beam produces a negative potential whose minimum is at the intensity maximum. If the potential tends to zero far from the beam and the minimum is $-U_0$, then $U_0>0$ is the nominal trap depth in the absence of gravity or other fields.

Bright traps are simple and provide strong curvature for a given optical geometry. Their main disadvantage is that trapped particles reside where the intensity, AC Stark shift, and scattering rate are largest.

Blue-detuned dark traps

Blue-detuned light produces positive energy where the intensity is high. Particles seek an intensity minimum. A single focused blue-detuned Gaussian beam is a barrier, not a three-dimensional trap. Confinement requires a dark region surrounded by bright walls, for example a hollow beam, a bottle beam, crossed light sheets, or an interference pattern.

A dark trap can reduce the average scattering and differential light shift because a cold particle spends most of its time near low intensity. Nevertheless, the barrier height, tunneling or escape channels, residual intensity at the minimum, and nonadiabatic polarization structure must all be checked.

(a) Radial and axial sections of one red-detuned Gaussian trap, plotted against the same distance unit $w_0$, with $z_R/w_0=5$. The axial well is much broader. (b) At fixed depth and waist, the isolated-transition model requires $P/P_*=|\Delta|/|\Delta_*|$ while $\Gamma_{\rm sc}/\Gamma_*=|\Delta_*|/|\Delta|$; the vertical scale is logarithmic. Every plotted detuning must remain far from the linewidth but small compared with the optical frequency; these are scaling curves, not an extrapolation through resonance.

A Single Gaussian-Beam Trap

Let $\rho^2=x^2+y^2$ measure distance from the propagation axis $z$. A TEM$_{00}$ Gaussian beam of power $P$, waist $w_0$, and wavelength $\lambda_L$ has

$$I(\rho,z) = \frac{ 2P }{ \pi w^2(z) } \exp \left[ - \frac{ 2\rho^2 }{ w^2(z) } \right],$$

with

$$w(z) = w_0 \sqrt{ 1+\frac{z^2}{z_R^2} }, \qquad z_R = \frac{ \pi w_0^2 }{ \lambda_L }.$$

For a red-detuned state, define the positive on-axis depth

$$U_0 = \frac{ \operatorname{Re}\alpha_a^{\mathrm{eff}} }{ 2\epsilon_0c } \frac{ 2P }{ \pi w_0^2 }.$$

Then

$$U(\rho,z) = - U_0 \frac{ \exp\left[-2\rho^2/w^2(z)\right] }{ 1+z^2/z_R^2 }.$$

This compact expression contains both confinement directions. Radial confinement comes from the transverse intensity profile; axial confinement comes from diffraction.

Harmonic expansion

For $\rho\ll w_0$ and $|z|\ll z_R$,

$$U(\rho,z) \simeq - U_0 + \frac{ 2U_0 }{ w_0^2 } \rho^2 + \frac{ U_0 }{ z_R^2 } z^2.$$

Matching to

$$U \simeq - U_0 + \frac12m\omega_r^2\rho^2 + \frac12m\omega_z^2z^2$$

gives

$$\omega_r = \sqrt{ \frac{ 4U_0 }{ mw_0^2 } }, \qquad \omega_z = \sqrt{ \frac{ 2U_0 }{ mz_R^2 } }.$$

Their aspect ratio is

$$\frac{ \omega_r }{ \omega_z } = \sqrt{2} \frac{ z_R }{ w_0 }.$$

Because $z_R/w_0=\pi w_0/\lambda_L$ is usually large, a single beam is typically a tight radial trap and a weak axial trap.

General normal modes

For an arbitrary collection of beams and external potentials, first find the equilibrium position $\mathbf r_0$ from

$$\left. \boldsymbol{\nabla}U_{\mathrm{tot}} \right|_{\mathbf r_0} = 0.$$

Then form the Hessian

$$K_{ij} = \left. \frac{ \partial^2U_{\mathrm{tot}} }{ \partial x_i\partial x_j } \right|_{\mathbf r_0}.$$

For a single particle of isotropic mass $m$, the eigenvalues of $K/m$ are $\omega_i^2$. Off-diagonal entries rotate the principal axes away from the laboratory axes. For particles with an effective mass tensor, diagonalize $M^{-1/2}KM^{-1/2}$ instead.

A positive-definite Hessian establishes a strict local harmonic minimum, not a sufficient escape barrier. A zero eigenvalue requires higher-order terms: a quartic well can confine even with zero quadratic curvature. Check the barrier along every escape direction as well as local stability.

Depth Is Not Curvature

The trap depth is the minimum energy required to reach an escape channel from the local minimum:

$$U_{\mathrm{depth}} = \min_{\mathcal P} \left[ \max_{\mathbf r\in\mathcal P} U_{\mathrm{tot}}(\mathbf r) - U_{\mathrm{tot}}(\mathbf r_0) \right],$$

where $\mathcal P$ ranges over paths from the minimum to the exterior; use a supremum if escape occurs asymptotically. This is a classical potential barrier. Quantum zero-point energy reduces the energy available for binding, and a finite barrier can permit tunneling. This definition matters for crossed beams, tilted traps, apertures, and multiwell landscapes. The shallowest saddle, not the highest local barrier, sets the actual depth.

The harmonic frequencies depend on local second derivatives. Two traps can have the same depth and very different curvature, or the same curvature and very different depth. Quoting only $U_0/k_{\mathrm B}$ therefore does not specify the trap.

Anharmonicity

The Gaussian potential softens away from its center. Oscillation frequency therefore decreases with energy, and the density of states departs from the harmonic result as the ensemble approaches the lip. Consequences include:

• broadened and shifted parametric resonances;
• energy-dependent dephasing of collective motion;
• non-Boltzmann truncation near the escape threshold;
• modified evaporation efficiency;
• failure of a temperature fit that assumes an infinite harmonic well.

The harmonic approximation is controlled when the occupied energy scale is small compared with the depth:

$$k_{\mathrm B}T \ll U_{\mathrm{depth}}$$

for a classical thermal cloud, with an analogous comparison using chemical potential and excitation energy for a quantum-degenerate gas.

Photon Scattering

The imaginary polarizability gives the photon-removal rate from the trapping mode:

$$\Gamma_{\mathrm{sc},a}(\mathbf r) = \frac{ \operatorname{Im} \alpha_a^{\mathrm{eff}}(\omega_L) }{ \hbar\epsilon_0c } I(\mathbf r).$$

For a closed radiative transition this is the spontaneous scattering rate, provided the response includes radiative damping. For atoms or molecules with photoionization, photodissociation, quenching, or other absorptive channels, the imaginary response must be resolved into the physically relevant loss channels.

Combining the conservative and dissipative parts gives

$$\frac{ \hbar\Gamma_{\mathrm{sc}} }{ |U| } = \frac{ 2\operatorname{Im}\alpha^{\mathrm{eff}} }{ |\operatorname{Re}\alpha^{\mathrm{eff}}| }.$$

In the far-detuned two-level limit,

$$\Gamma_{\mathrm{sc}} \simeq \Gamma \frac{ |\Omega|^2 }{ 4\Delta^2 },$$

so that

$$\frac{ \hbar\Gamma_{\mathrm{sc}} }{ |U| } \simeq \frac{ \Gamma }{ |\Delta| }.$$

At fixed laser intensity,

$$|U| \propto \frac{I}{|\Delta|}, \qquad \Gamma_{\mathrm{sc}} \propto \frac{I}{\Delta^2}.$$

At fixed trap depth, the required intensity scales as $I\propto|\Delta|$, and therefore

$$\Gamma_{\mathrm{sc}} \propto \frac{1}{|\Delta|}.$$

Moving farther from resonance buys lower scattering only by demanding more optical power or tighter focusing. The simple scaling eventually fails when other resonances, counter-rotating terms, molecular channels, or laser damage become important.

Elastic and Raman scattering

Off-resonant scattering is not a single error channel:

• Rayleigh scattering returns the particle to the same internal state but changes momentum and can reveal state information through the photon.
• Raman scattering changes hyperfine, Zeeman, vibrational, rotational, or electronic state.
• A Raman event can move the particle into an anti-trapped state or outside the addressed qubit manifold.
• Even state-preserving scattering can dephase a superposition when the emitted photon distinguishes its components.

A total scattering rate is therefore insufficient for a coherence budget. Branching ratios and state-resolved polarizabilities are required.

Recoil Heating

For wave number $k_L=2\pi/\lambda_L$, the single-photon recoil energy is

$$E_R = \frac{ \hbar^2k_L^2 }{ 2m }.$$

For traveling-wave absorption, zero mean momentum along the beam, and uncorrelated spontaneous-emission directions, absorption and emission each add one recoil energy on average. For an approximately closed transition with comparable absorbed and emitted wavelengths,

$$\frac{d\langle E\rangle}{dt} \simeq 2E_R \left\langle \Gamma_{\mathrm{sc}} \right\rangle.$$

This is a total three-dimensional energy rate. Its division among axes depends on beam geometry and emission pattern. It is not directly $k_{\rm B} \dot T$: an equilibrated three-dimensional classical harmonic trap has $E=3k_{\rm B}T$ per atom. In the Lamb–Dicke regime use recoil matrix elements and sideband-dependent rates; the mean recoil sum can still give $2E_R$ when all relevant final motional states are included and scattering rates are effectively independent of their energy.

Photon scattering can cause loss long before recoil heating raises the mean energy to the nominal trap depth. A Raman event, optical pumping into a dark anti-trapped state, or a collision following internal excitation may remove a particle immediately.

Gravity and Other Tilts

Take $z$ upward. Gravity adds

$$U_g(z) = mgz.$$

Near a harmonic minimum,

$$U_{\mathrm{tot}}(z) \simeq \frac12m\omega_z^2z^2 + mgz,$$

so the equilibrium displacement is

$$z_{\mathrm{sag}} = - \frac{ g }{ \omega_z^2 }.$$

This local result is valid only if a true equilibrium exists. The full optical restoring force must be able to exceed $mg$ somewhere. For a red-detuned Gaussian profile along a radial vertical coordinate,

$$|F_{\mathrm{opt}}|_{\max} = \frac{ 2U_0 }{ w_0\sqrt e }.$$

Thus radial support requires

$$mg < \frac{ 2U_0 }{ w_0\sqrt e }.$$

Along the propagation axis of a single Gaussian beam,

$$|F_{\mathrm{opt}}|_{\max} = \frac{ 9 }{ 8\sqrt3 } \frac{ U_0 }{ z_R }.$$

The weak axial direction may therefore fail to support the same atom even when its harmonic sag appears moderate. Beam orientation is part of the trap design.

Magnetic gradients, electric fields, surface forces, and optical radiation-pressure offsets can tilt the potential in the same way. The effective depth must be recomputed from the lowest saddle after all such terms are included.

Loading a Conservative Trap

Hamiltonian evolution preserves fine-grained phase-space density. A slow decompression can reduce temperature by removing work without increasing phase-space density or lowering entropy. To obtain a dense, cold initial ensemble, practical loading commonly starts with dissipative preparation such as a MOT, sub-Doppler molasses, buffer-gas cooling, or laser-cooled molecules.

After cooling light and time-dependent ramps are removed, a sufficient classical condition for confinement below every escape saddle is

$$\frac{ p^2 }{ 2m } + U_{\mathrm{tot}}(\mathbf r) < U_{\mathrm{saddle}}$$

for the final static potential. Particles above the saddle can sometimes remain for long times because of dynamical restrictions; the inequality is a sufficient energetic criterion, not a complete prediction of finite-time capture. During cooling and ramps, energy is not conserved. Loading depends on:

1. spatial overlap between the cold cloud and optical well;
2. velocity distribution relative to the trap depth;
3. internal-state polarizabilities and optical pumping;
4. timing and adiabaticity of the trap ramp;
5. light-assisted collisions while near-resonant cooling light remains on;
6. gravity and other escape channels.

Turning on the trap suddenly changes particle energies according to their instantaneous positions. A sufficiently slow ramp approximately preserves motional actions rather than energy. Neither procedure increases fine-grained phase-space density without dissipation or particle loss.

Further cooling through energy-selective escape and collisional rethermalization is treated in evaporative cooling.

Technical Heating and Loss

Even when photon scattering is negligible, a noisy conservative trap can heat rapidly.

Intensity noise

Fractional intensity fluctuations modulate both depth and curvature:

$$U_0(t) = \bar U_0 \left[ 1+\epsilon(t) \right].$$

Because $\omega_i^2\propto U_0$, noise near twice a trap frequency drives parametric excitation:

$$\text{intensity-noise sensitivity} \ \text{peaks near}\ 2\omega_i.$$

The resulting heating is often exponential in mean motional energy. Its coefficient is proportional to $\omega_i^2S_\epsilon(2\omega_i)$, with numerical factors set by whether the power spectral density is one-sided or two-sided and whether frequency is measured in hertz or radians per second.

Pointing noise

Beam-position fluctuations move the trap center. Noise near $\omega_i$ resonantly drives sloshing:

$$\text{pointing-noise sensitivity} \ \text{peaks near}\ \omega_i.$$

The heating scale grows strongly with confinement, approximately as $m\omega_i^4S_{x_i}(\omega_i)$ for a fixed displacement-noise spectrum. A stiffer trap is therefore more sensitive to the same pointing jitter. With gravity or other tilts, changing intensity also moves the equilibrium, so intensity noise can excite both curvature and center-motion resonances. The two mechanisms cannot always be assigned to separate hardware sources.

Other loss channels

A complete lifetime model may include

$$\dot N = - \Gamma_1N - K_2 \int n^2(\mathbf r)\,d^3r - K_3 \int n^3(\mathbf r)\,d^3r.$$

Here $\Gamma_1$ includes one-body loss (which can be state dependent), while $K_2$ and $K_3$ are atom-loss coefficients for two- and three-body processes; multiplicities per event are included in their definitions. Measured decay need not be exponential when density-dependent loss or evaporation is important.

Additional technical channels include:

• laser-frequency noise, especially near a polarizability resonance;
• polarization noise and vector-light-shift fluctuations;
• interference fringes from unintended reflections;
• relative phase noise in crossed or standing-wave traps;
• slow waist, focus, or alignment drift;
• servo transients during power ramps.

State-Dependent Potentials

Two internal states $|a\rangle$ and $|b\rangle$ generally experience different potentials:

$$U_a(\mathbf r) \neq U_b(\mathbf r).$$

The local transition shift is

$$\delta\omega_{ba}(\mathbf r) = \frac{ U_b(\mathbf r)-U_a(\mathbf r) }{ \hbar }.$$

Spatial variation of this differential shift produces inhomogeneous broadening and dephasing. Motion through the trap can convert it into a time-dependent frequency modulation.

Magic conditions

A magic trapping condition for a chosen transition satisfies

$$\operatorname{Re}\alpha_b^{\mathrm{eff}} = \operatorname{Re}\alpha_a^{\mathrm{eff}}$$

to the required order and for the specified polarization and quantization geometry. Then the leading differential light shift vanishes. A magic wavelength is not automatically magic against:

• polarization errors;
• magnetic-field changes;
• hyperpolarizability;
• multipolar light shifts;
• motion-dependent sampling;
• photon-scattering differences.

Tune-out is different

A tune-out frequency for state $|a\rangle$ satisfies

$$\operatorname{Re}\alpha_a^{\mathrm{eff}} = 0.$$

It is a zero of one state’s conservative response, not an equality between two states. The imaginary polarizability need not vanish there, so zero dipole potential does not imply zero scattering.

Vector light shifts act like polarization-dependent effective magnetic fields, while tensor shifts depend on the alignment of angular momentum. For spin coherence, polarization purity and the direction of the bias field are therefore part of the trap calibration.

Measuring Trap Parameters

No single measurement determines depth, frequency, temperature, and scattering. Useful diagnostics are deliberately redundant.

Trap frequencies

• Center modulation: shake the trap position and identify resonant response near $\omega_i$.
• Parametric modulation: modulate intensity and identify the principal loss or heating resonance near $2\omega_i$.
• Sloshing: displace the equilibrium position, release it, and fit the center-of-mass oscillation.
• Resolved sidebands: in sufficiently tight, cold traps, use red and blue motional sidebands to determine $\omega_i$ and occupation.

Large disagreement among these methods usually signals anharmonicity, mode coupling, or a calibration error.

Depth and temperature

Trap depth can be inferred by controlled spilling, a calibrated static tilt, release-and-recapture modeling, or measured beam parameters combined with a validated polarizability. Each method has model dependence.

Temperature measurements include time of flight, release and recapture, in-situ size with calibrated curvature, and sideband thermometry. A truncated distribution in a shallow trap should not be fit blindly to an unbounded Maxwell–Boltzmann distribution.

Scattering and coherence

Scattering can be constrained through state-changing rates, recoil heating, fluorescence, trap lifetime, or decoherence. These observables measure different combinations of Rayleigh and Raman channels. A lifetime longer than the experimental sequence does not by itself establish a small coherence error.

A 1064 nm rubidium trap

Consider a single red-detuned Gaussian beam with assumed depth

$$\frac{ U_0 }{ k_{\mathrm B} } = 1.0\,\mathrm{mK},$$

waist

$$w_0 = 50\,\mu\mathrm m,$$

wavelength

$$\lambda_L = 1064\,\mathrm{nm},$$

and an $^{87}\mathrm{Rb}$ atom of mass

$$m = 1.443\times10^{-25}\,\mathrm{kg}.$$

This example assumes the depth rather than deriving it from laser power; the latter requires a state- and polarization-resolved polarizability.

The Rayleigh range is

$$z_R = \frac{ \pi w_0^2 }{ \lambda_L } \simeq 7.38\,\mathrm{mm}.$$

The harmonic frequencies are

$$\frac{ \omega_r }{ 2\pi } \simeq 1.97\,\mathrm{kHz}, \qquad \frac{ \omega_z }{ 2\pi } \simeq 9.43\,\mathrm{Hz}.$$

Thus the trap is highly anisotropic:

$$\frac{ \omega_r }{ \omega_z } \simeq 209.$$

If gravity acts along a radial direction, the harmonic sag is only

$$|z_{\mathrm{sag}}| = \frac{ g }{ \omega_r^2 } \simeq 64\,\mathrm{nm}.$$

The maximum radial optical acceleration is about

$$a_{r,\max} \simeq 2.3\times10^3\,\mathrm{m\,s^{-2}},$$

well above $g$. If the beam instead points vertically so that gravity acts along its weak axial direction, the maximum axial optical acceleration is only

$$a_{z,\max} \simeq 8.4\,\mathrm{m\,s^{-2}},$$

which is smaller than $g$. The full potential then has no supported axial equilibrium even though inserting $\omega_z$ into the harmonic sag formula would return a finite number. This is why the force-maximum test must precede the sag calculation.

The 1064 nm recoil temperature is

$$\frac{ E_R }{ k_{\mathrm B} } \simeq 97\,\mathrm{nK}.$$

If a separately calculated scattering rate were $\Gamma_{\mathrm{sc}}=0.5\,\mathrm{s^{-1}}$, the recoil-heating estimate would be

$$\frac{1}{k_{\mathrm B}} \frac{d\langle E\rangle}{dt} \simeq 97\,\mathrm{nK\,s^{-1}}.$$

That rate is not a universal property of 1064 nm light. It depends on species, state, polarization, and intensity.

From trapping to coherent control

Tight focusing produces optical tweezers; beam interference produces optical lattices. Their localization and tunneling scales follow from the same potential, but single-particle imaging and many-site dynamics require additional models. In precision spectroscopy, equal real polarizabilities cancel the leading differential shift while leaving scattering and higher-order shifts possible. A long particle lifetime therefore does not establish long internal coherence.

Exercises

1. Determine the trapping region

For the convention $\Delta=\omega_0-\omega_L$, use

$$U_g = - \frac{ \hbar|\Omega|^2 }{ 4\Delta }$$

to determine whether the ground state seeks high or low intensity for red and blue detuning. Why does a single blue-detuned Gaussian beam fail to trap in three dimensions?

Solution

For red detuning, $\Delta>0$, so $U_g<0$ and becomes more negative as $|\Omega|^2\propto I$ increases. The atom seeks high intensity.

For blue detuning, $\Delta<0$, so $U_g>0$ and increases with intensity. The atom seeks low intensity. A focused Gaussian beam has a maximum, not a minimum, at its focus. Blue-detuned atoms are repelled from that point and can escape around the beam. A dark trap requires a low-intensity region surrounded in all directions by higher-intensity barriers.

2. Cooling without increasing phase-space density

An isolated classical gas occupies an isotropic harmonic trap of frequency $\omega$. Slowly reduce the frequency to $q\omega$, with $0<q<1$, while preserving the motional actions and neglecting loss. Find the final temperature, cloud width, and peak phase-space density.

Solution

For each mode $J=E/\omega$ is adiabatically conserved, so $E'=qE$ and $T'=qT$. The rms width obeys $\sigma'^2=k_{\rm B}T'/(m\omega'^2)=\sigma^2/q$. At fixed particle number, $n_0'=q^{3/2}n_0$, while $\lambda_{\rm dB}'=\lambda_{\rm dB}/\sqrt q$. Their product $n_0'\lambda_{\rm dB}'^3=n_0\lambda_{\rm dB}^3$ is unchanged. The temperature falls by adiabatic work extraction without entropy removal or increased degeneracy. A finite-depth trap additionally needs negligible spilling throughout this ramp.

3. Maximum force of a Gaussian trap

Find the largest restoring forces of $U_r(x)=-U_0e^{-2x^2/w_0^2}$ and $U_z(z)=-U_0/(1+z^2/z_R^2)$. At what distances do they occur?

Solution

The positive force magnitudes at positive coordinates are

$$|F_r|=\frac{4U_0x}{w_0^2}e^{-2x^2/w_0^2},\qquad |F_z|=\frac{2U_0z/z_R^2}{(1+z^2/z_R^2)^2}.$$

Their derivatives vanish at $x=w_0/2$ and $z=z_R/\sqrt3$. Substitution gives $|F_r|_{\max}=2U_0/(w_0\sqrt e)$ and $|F_z|_{\max}=9U_0/(8\sqrt3z_R)$. These finite maxima determine whether gravity can be balanced; the curvature at the origin alone cannot do so.

4. Fixed-depth detuning scaling

In a two-level model,

$$|U| \propto \frac{I}{|\Delta|}, \qquad \Gamma_{\mathrm{sc}} \propto \frac{I}{\Delta^2}.$$

If $|\Delta|$ is increased by a factor of $5$ while trap depth is held fixed, how must intensity change, and what happens to the scattering rate?

Solution

Fixed depth requires

$$\frac{I'}{5|\Delta|} = \frac{I}{|\Delta|},$$

so $I'=5I$. The new scattering rate is

$$\frac{ \Gamma'_{\mathrm{sc}} }{ \Gamma_{\mathrm{sc}} } = \frac{ 5I }{ (5\Delta)^2 } \frac{ \Delta^2 }{ I } = \frac15.$$

Thus the required optical power rises by a factor of $5$, while scattering falls by a factor of $5$. This conclusion assumes the same two-level far-detuned regime remains valid.

5. Recoil-heating estimate

An atom has $E_R/k_{\mathrm B}=100\,\mathrm{nK}$ and scatters trap photons at $0.20\,\mathrm{s^{-1}}$. Estimate the three-dimensional recoil-heating rate in temperature units.

Solution

Using

$$\frac{ d\langle E\rangle }{ dt } \simeq 2E_R\Gamma_{\mathrm{sc}},$$

one obtains

$$\frac{1}{k_{\mathrm B}} \frac{ d\langle E\rangle }{ dt } \simeq 2 \left( 100\,\mathrm{nK} \right) \left( 0.20\,\mathrm{s^{-1}} \right) = 40\,\mathrm{nK\,s^{-1}}.$$

This is an energy-growth scale. If an equilibrated classical three-dimensional harmonic gas retains all particles, $E=3k_{\rm B}T$ gives $\dot T\simeq13.3\,\mathrm{nK\,s^{-1}}$. A nonthermal cloud has no unique thermodynamic temperature, and selective loss changes this conversion.

6. Sag versus support

For a vertical harmonic frequency $\omega/(2\pi)=20\,\mathrm{Hz}$, calculate the harmonic sag. Then explain why this number alone cannot establish that the trap supports the atom.

Solution

The angular frequency is

$$\omega = 2\pi \left( 20\,\mathrm{Hz} \right) \simeq 126\,\mathrm{s^{-1}}.$$

Therefore

$$|z_{\mathrm{sag}}| = \frac{ 9.81 }{ (126)^2 } \simeq 6.2\times10^{-4}\,\mathrm m = 0.62\,\mathrm{mm}.$$

The harmonic approximation extends indefinitely and therefore always returns a stationary point. A real optical well has a bounded restoring force and a finite escape barrier. One must verify $|F_{\mathrm{opt}}|_{\max}>mg$ using the full potential; otherwise no equilibrium exists and the sag result is physically irrelevant.

7. Identify the noise channel

An experiment measures trap frequencies

$$\frac{\omega_x}{2\pi} = 1.2\,\mathrm{kHz}, \qquad \frac{\omega_z}{2\pi} = 90\,\mathrm{Hz}.$$

Heating peaks appear when intensity noise is injected at $2.4\,\mathrm{kHz}$ and when beam pointing is modulated at $90\,\mathrm{Hz}$. Interpret both peaks for an approximately harmonic trap whose equilibrium does not shift with intensity.

Solution

Intensity modulation changes the curvature and drives parametric heating near twice a normal-mode frequency:

$$2 \frac{ \omega_x }{ 2\pi } = 2.4\,\mathrm{kHz}.$$

The first peak therefore identifies the $x$ mode.

Pointing modulation moves the trap center and directly drives motion near the normal-mode frequency. The $90\,\mathrm{Hz}$ peak identifies the $z$ sloshing mode. The different resonance conditions help separate intensity noise from position noise.

8. Magic, tune-out, and scattering

At one wavelength,

$$\operatorname{Re}\alpha_a = \operatorname{Re}\alpha_b \neq 0.$$

At a second wavelength,

$$\operatorname{Re}\alpha_a = 0, \qquad \operatorname{Re}\alpha_b \neq 0.$$

Classify the two wavelengths for the $a\leftrightarrow b$ transition. Does either equation guarantee zero photon scattering?

Solution

The first condition is magic for the $a\leftrightarrow b$ transition at leading order because

$$U_b-U_a \propto - \left( \operatorname{Re}\alpha_b - \operatorname{Re}\alpha_a \right) I = 0.$$

The second is a tune-out condition for state $|a\rangle$, but it is not magic for the transition because state $|b\rangle$ still shifts.

Neither equation constrains the imaginary polarizabilities. Since $\Gamma_{\mathrm{sc}}\propto\operatorname{Im}\alpha$, photon scattering may remain nonzero at either wavelength.

References

• B. Arora, M. S. Safronova, and C. W. Clark, “Tune-out wavelengths of alkali-metal atoms and their applications,” Physical Review A 84, 043401 (2011), doi:10.1103/PhysRevA.84.043401.
• A. Ashkin, “Trapping of atoms by resonance radiation pressure,” Physical Review Letters 40, 729–732 (1978), doi:10.1103/PhysRevLett.40.729.
• M. D. Barrett, J. A. Sauer, and M. S. Chapman, “All-optical formation of an atomic Bose–Einstein condensate,” Physical Review Letters 87, 010404 (2001), doi:10.1103/PhysRevLett.87.010404.
• C. J. Foot, Atomic Physics (Oxford University Press, 2005), doi:10.1093/oso/9780198506959.001.0001.
• M. E. Gehm, K. M. O’Hara, T. A. Savard, and J. E. Thomas, “Dynamics of noise-induced heating in atom traps,” Physical Review A 58, 3914–3921 (1998), doi:10.1103/PhysRevA.58.3914.
• R. Grimm, M. Weidemüller, and Y. B. Ovchinnikov, “Optical dipole traps for neutral atoms,” Advances in Atomic, Molecular, and Optical Physics 42, 95–170 (2000), doi:10.1016/S1049-250X(08)60186-X.
• H. Katori, M. Takamoto, V. G. Pal’chikov, and V. D. Ovsiannikov, “Ultrastable optical clock with neutral atoms in an engineered light shift trap,” Physical Review Letters 91, 173005 (2003), doi:10.1103/PhysRevLett.91.173005.
• S. J. M. Kuppens, K. L. Corwin, K. W. Miller, T. E. Chupp, and C. E. Wieman, “Loading an optical dipole trap,” Physical Review A 62, 013406 (2000), doi:10.1103/PhysRevA.62.013406.
• H. J. Metcalf and P. van der Straten, Laser Cooling and Trapping (Springer, 1999), doi:10.1007/978-1-4612-1470-0.
• J. D. Miller, R. A. Cline, and D. J. Heinzen, “Far-off-resonance optical trapping of atoms,” Physical Review A 47, R4567–R4570 (1993), doi:10.1103/PhysRevA.47.R4567.
• T. A. Savard, K. M. O’Hara, and J. E. Thomas, “Laser-noise-induced heating in far-off resonance optical traps,” Physical Review A 56, R1095–R1098 (1997), doi:10.1103/PhysRevA.56.R1095.