Sideband Cooling

A trapped ion sits in a harmonic well, and like any harmonic oscillator its motion is quantized: the ion does not wiggle with an arbitrary amplitude but occupies discrete vibrational levels En=ℏωz(n+12)E_n=\hbar\omega_z(n+\tfrac12) separated by one phonon ℏωz\hbar\omega_z. Ordinary laser (Doppler) cooling brings that motion down to a thermal cloud of order ten phonons — cold by the standards of atomic physics, but far too hot for a quantum computer, where the shared vibration is the wire that carries information between ion qubits. Sideband cooling is the technique that empties the wire. It exploits a beautiful piece of spectroscopy: because the ion's internal transition is frequency-modulated by its own oscillation, the absorption line sprouts motional sidebands one phonon above and below the carrier, and a laser tuned to the lower (red) sideband can only be absorbed if the ion simultaneously surrenders a quantum of motion. Three facts organize everything that follows. First, the trick needs resolution: the sidebands are only distinct lines if the trap frequency beats the transition linewidth, ωz≫γ\omega_z\gg\gamma — the resolved-sideband regime. Second, it needs smallness: the ion's zero-point wavepacket must be tiny compared with the optical wavelength, measured by the Lamb–Dicke parameter η=kz0≪1\eta=k z_0\ll1, so that each photon changes the phonon number by at most one. Third, it has a natural stop: the state ∣g,0⟩\lvert g,0\rangle has no lower rung, so the red-sideband coupling ∝n\propto\sqrt{n} vanishes there — the ground state is dark, and population piles up in it. The residual heating from off-resonant scattering fixes a floor nˉss≈(γ/2ωz)2\bar n_{\mathrm{ss}}\approx(\gamma/2\omega_z)^2, routinely below 0.050.05 in modern traps. This is the article that turns a warm ion into the cold, well-defined motional bus the DiVincenzo criteria demand.

Key takeaways
  • Sideband cooling drives the red motional sideband of an internal transition, so each absorption–emission cycle removes one phonon ℏωz\hbar\omega_z and pumps a trapped ion or atom toward its motional ground state ∣n=0⟩\lvert n=0\rangle.
  • It works only in the resolved-sideband regime, ωz≫γ\omega_z\gg\gamma, where the trap frequency exceeds the transition linewidth so the red sideband is spectrally distinct from the carrier and the heating (blue) sideband.
  • The Lamb–Dicke parameter η=kz0=ℏk2/2mωz\eta=k z_0=\sqrt{\hbar k^2/2m\omega_z} controls everything: for η≪1\eta\ll1 the recoil rarely changes nn, and the sideband couplings are Ωn,n−1=Ωηn\Omega_{n,n-1}=\Omega\eta\sqrt{n}.
  • The steady-state occupation is nˉss≈(γ/2ωz)2≪1\bar n_{\mathrm{ss}}\approx(\gamma/2\omega_z)^2\ll1, set by balancing resonant cooling against off-resonant reheating — narrow transitions and stiff traps push the floor down quadratically.
  • Ground-state cooling is the mandatory first step for phonon-mediated entangling gates (Cirac–Zoller, Mølmer–Sørensen), for optical-clock state preparation, and for the DiVincenzo cold-qubit requirement.
  • Key idea: modulating an internal transition by harmonic motion produces resolved motional sidebands; resonantly driving the red sideband (ωL=ω0−ωz\omega_L=\omega_0-\omega_z) lets each absorption–emission cycle subtract one phonon, pumping the system to its dark motional ground state ∣g,0⟩\lvert g,0\rangle.

Where this lives today

Sideband cooling is the indispensable opening move of every trapped-ion experiment, and its reach keeps widening. In quantum computers built from ion chains — the machines of IonQ, Quantinuum, and academic groups at Innsbruck, Maryland, and NIST — every entangling gate begins by cooling the shared motional modes to nˉ≪1\bar n\ll1; residual phonons are a leading error source, so groups now use electromagnetically-induced-transparency (EIT) and polarization-gradient cooling to chill many modes of a long chain at once, faster than pulsed sideband cooling allows. In the world's best optical clocks — the Al+\mathrm{Al}^+ quantum-logic clock and Yb+\mathrm{Yb}^+ and Sr+\mathrm{Sr}^+ standards reaching fractional uncertainties near 10−1810^{-18} — the same cooling suppresses the motional (second-order Doppler) time-dilation shift that would otherwise dominate the error budget. And the resolved-sideband idea has jumped substrate entirely: cavity optomechanics uses a red-detuned laser to cool gram-to-nanogram mechanical oscillators toward their quantum ground state, and circuit-QED microwave drives cool superconducting and mechanical resonators. An open frontier is scalable cooling — keeping the bus cold in traps holding tens to hundreds of ions, and re-cooling spectator modes mid-algorithm without disturbing the qubits, a bottleneck fault-tolerant ion-trap architectures must solve.

The physical picture

Picture the ion as a tiny bell whose pitch ω0\omega_0 is faintly frequency-modulated by its swing in the trap at ωz\omega_z. That modulation dresses the single absorption line with a comb: a strong carrier at ω0\omega_0 and weak satellites — motional sidebands — one phonon below at ω0−ωz\omega_0-\omega_z and one above at ω0+ωz\omega_0+\omega_z. Now shine a laser tuned exactly to the lower satellite. A photon at ω0−ωz\omega_0-\omega_z is one phonon short of the internal transition, so it can only be absorbed if the ion makes up the deficit from its own motion — the ion must give up a vibrational quantum to get excited. That is the cooling event: ∣g,n⟩→∣e,n−1⟩\lvert g,n\rangle\to\lvert e,n-1\rangle. The ion then re-emits a photon spontaneously, and because the recoil is tiny (Lamb–Dicke regime), the most probable emission leaves the phonon number untouched, returning the ion to ∣g,n−1⟩\lvert g,n-1\rangle. Net result of the round trip: one phonon gone.

Figure
The cooling cycle on the ladder of dressed states. A red-sideband photon drives the ion up one internal level and down one motional rung (n to n−1), then spontaneous emission returns it at the same motional level (still n−1). Each full cycle removes one phonon ħω_z so Δn = −1, and the cycle stalls at the motional ground state because there is no lower rung.
One cooling cycle removes one phonon: |g,n⟩ → |e,n−1⟩ → |g,n−1⟩, Δn = −1 ground |g,n⟩ excited |e,n⟩ n=0n=1n=2n=3 n=0n=1n=2n=3 ħω_z red sideband ω_L = ω₀ − ω_z spontaneous emission Δn ≈ 0 net −1 |g,0⟩ dark: √n = 0, no lower rung

Each cycle is a ratchet, stepping the ion down the vibrational ladder until it reaches ∣g,0⟩\lvert g,0\rangle. There the ratchet jams: the red-sideband coupling scales as n\sqrt{n}, so it vanishes at n=0n=0 — there is no n=−1n=-1 state to fall into, the laser can no longer be absorbed, and population accumulates in the ground state, which is dark to the cooling light. The whole scheme lives or dies on spectral resolution. If the excited state is short-lived, its linewidth γ\gamma smears the levels; when γ≳ωz\gamma\gtrsim\omega_z the red and blue satellites blur into the carrier, and the laser starts driving the carrier (which does not cool) and the blue sideband (which heats). Only when the trap quantum ℏωz\hbar\omega_z comfortably exceeds the energy blur ℏγ\hbar\gamma do the sidebands stand apart as clean, separately addressable lines. That is what the second figure shows, with real numbers.

Figure
The resolved-sideband absorption spectrum for a typical trapped-ion optical qubit: three sharp lines — red sideband, carrier, blue sideband — spaced by the trap frequency ω_z/2π = 1 MHz, each of width γ/2π = 20 kHz. Because ω_z ≫ γ the satellites are cleanly resolved, and the residual cooling floor is n̄ss ≈ (γ/2ω_z)² ≈ 1×10⁻⁴.
Resolved sidebands: ω_z ≫ γ makes the red satellite separately addressable laser detuning δ = ω_L − ω₀ → ω₀ − ω_z red (cools) ω₀ carrier ω₀ + ω_z blue (heats) ω_z/2π = 1 MHz γ/2π = 20 kHz (γ ≪ ω_z) tune laser here → floor n̄ss ≈ (γ/2ω_z)² ≈ 1×10⁻⁴

The core equations

The starting point is the Hamiltonian of a two-level ion whose motion along one trap axis is quantized, coupled to a travelling laser field. The spatial phase of that field, evaluated at the ion's operator-valued position z^=z0(a+a†)\hat z=z_0(a+a^\dagger), is what mixes internal and motional degrees of freedom and generates the sidebands.

H=ℏωz a†a+ℏω02σz+ℏΩ2(σ+ei(η(a+a†)−ωLt)+h.c.)(1)H = \hbar\omega_z\, a^\dagger a + \frac{\hbar\omega_0}{2}\sigma_z + \frac{\hbar\Omega}{2}\left(\sigma_+ e^{i(\eta(a+a^\dagger)-\omega_L t)} + \mathrm{h.c.}\right) \qquad (1)

where: ωz\omega_z — trap (axial) frequency; a,a†a,a^\dagger — motional annihilation/creation operators; ω0\omega_0 — internal transition frequency; σz,σ±\sigma_z,\sigma_\pm — qubit Pauli and raising/lowering operators; Ω\Omega — bare carrier Rabi frequency; η\eta — Lamb–Dicke parameter; ωL\omega_L — laser frequency. The exponential eiη(a+a†)e^{i\eta(a+a^\dagger)} is the whole story: expand it and the sidebands fall out. The single dimensionless number that governs the coupling is the Lamb–Dicke parameter — the ratio of the ion's zero-point motion z0z_0 to the reduced laser wavelength 1/k1/k, or equivalently the square root of the photon recoil energy measured in phonons.

η=kz0=kℏ2mωz=ℏk22mωz(2)\eta = k z_0 = k\sqrt{\frac{\hbar}{2 m \omega_z}} = \sqrt{\frac{\hbar k^2}{2 m \omega_z}} \qquad (2)

where: kk — laser wavevector projected on the motional axis; z0=ℏ/2mωzz_0=\sqrt{\hbar/2m\omega_z} — ground-state position spread; mm — ion mass. Small η\eta means the wavepacket is far smaller than a wavelength, so a photon's recoil kick almost never changes nn — the physical content of the Lamb–Dicke regime. Expanding the phase factor in powers of η\eta sorts the ion–laser coupling into a carrier and a hierarchy of sidebands. Each successive term shifts the phonon number by one more unit, and in the Lamb–Dicke regime only the first two terms survive.

eiη(a+a†)≈1+iη(a+a†)−η22(a+a†)2+⋯(3)e^{i\eta(a+a^\dagger)} \approx 1 + i\eta(a+a^\dagger) - \frac{\eta^2}{2}(a+a^\dagger)^2 + \cdots \qquad (3)

where: the 11 term is the carrier (Δn=0\Delta n=0); iη(a+a†)i\eta(a+a^\dagger) gives the first red (aa, Δn=−1\Delta n=-1) and blue (a†a^\dagger, Δn=+1\Delta n=+1) sidebands; the η2\eta^2 term is a Debye–Waller correction and higher sidebands. Truncating after first order is exactly the resolved two-state picture that makes cooling analytically tractable. Reading off the matrix elements of the first-order term between Fock states gives the effective Rabi frequency on each sideband. The bosonic factors n\sqrt{n} and n+1\sqrt{n+1} are the mathematical heart of why ∣0⟩\lvert 0\rangle is dark.

Ωn,n−1=Ω ηn  (red),Ωn,n+1=Ω ηn+1  (blue)(4)\Omega_{n,n-1}=\Omega\,\eta\sqrt{n}\ \ (\text{red}),\qquad \Omega_{n,n+1}=\Omega\,\eta\sqrt{n+1}\ \ (\text{blue}) \qquad (4)

where: Ωn,n∓1\Omega_{n,n\mp1} — effective Rabi frequency coupling neighbouring Fock states; the factor n\sqrt{n} comes from a∣n⟩=n∣n−1⟩a\lvert n\rangle=\sqrt{n}\lvert n-1\rangle. At n=0n=0 the red-sideband coupling is exactly zero, so the ground state cannot absorb the cooling light — it is dark. Balancing the resonant cooling rate against off-resonant reheating fixes the lowest occupation the method can reach. In the resolved-sideband limit this floor scales as the square of the ratio of linewidth to trap frequency.

nˉss≈(γ2ωz)2(ωz≫γ)(5)\bar{n}_{\mathrm{ss}} \approx \left(\frac{\gamma}{2\omega_z}\right)^2 \qquad (\omega_z\gg\gamma) \qquad (5)

where: γ\gamma — natural linewidth (spontaneous decay rate) of the cooling transition; nˉss\bar n_{\mathrm{ss}} — steady-state mean phonon number. A narrower transition or a stiffer trap lowers the floor quadratically — the design rule behind narrow-line and Raman sideband cooling. Before sideband cooling can even begin, the ion must be pre-cooled below the point where higher sidebands proliferate. Doppler cooling on a strong dipole transition does this, reaching a mean occupation set by the linewidth-to-trap ratio.

kBTD=ℏγ2 ⇒ nˉDoppler≈kBTDℏωz=γ2ωz(6)k_B T_D=\frac{\hbar\gamma}{2}\ \Rightarrow\ \bar{n}_{\mathrm{Doppler}}\approx\frac{k_B T_D}{\hbar\omega_z}=\frac{\gamma}{2\omega_z} \qquad (6)

where: TDT_D — Doppler temperature; kBk_B — Boltzmann constant. Because Doppler cooling uses a broad transition with γ≳ωz\gamma\gtrsim\omega_z, it leaves nˉ\bar n of order ten — cold enough to enter the Lamb–Dicke regime, but not the ground state. Finally, driving the red sideband coherently produces Rabi flopping whose amplitude carries the n\sqrt{n} signature, and this is both the cooling engine and the thermometer. The vanishing of the signal at n=0n=0 is the experimental fingerprint of a cold ion.

Persb(t)=sin⁡2 ⁣(Ωηn2 t),nˉ=PrsbPbsb−Prsb(7)P_e^{\mathrm{rsb}}(t)=\sin^2\!\left(\frac{\Omega\eta\sqrt{n}}{2}\,t\right),\qquad \bar n=\frac{P_{\mathrm{rsb}}}{P_{\mathrm{bsb}}-P_{\mathrm{rsb}}} \qquad (7)

where: PersbP_e^{\mathrm{rsb}} — excited-state probability under red-sideband drive; Prsb,PbsbP_{\mathrm{rsb}},P_{\mathrm{bsb}} — red- and blue-sideband excitation after a fixed pulse. The red sideband is silent for n=0n=0, so its ratio to the blue sideband reads out nˉ\bar n with no calibration — the standard sideband thermometry.

Derivation ASideband structure from the Lamb–Dicke expansion

The goal is to show, starting from the ion–laser coupling of equation (1), exactly how the single exponential eiη(a+a†)e^{i\eta(a+a^\dagger)} splits the absorption line into a carrier and first-order red and blue sidebands, and to pin down the red-sideband Rabi frequency Ωn,n−1=Ωηn\Omega_{n,n-1}=\Omega\eta\sqrt{n} that both cools the ion and, later, entangles two of them. We move to the interaction picture, expand to first order in the small parameter η\eta, pick out the resonant term by choosing the laser detuning, and evaluate the surviving matrix element on Fock states.

Step 1. Write the interaction-picture coupling Hint=ℏΩ2 σ+ eiη(a+a†)e−iδt+h.c.H_{\mathrm{int}}=\tfrac{\hbar\Omega}{2}\,\sigma_+\,e^{i\eta(a+a^\dagger)}e^{-i\delta t}+\mathrm{h.c.}, where δ=ωL−ω0\delta=\omega_L-\omega_0 is the laser detuning from the carrier and the motional operators rotate as a→a e−iωzta\to a\,e^{-i\omega_z t}, a†→a†e+iωzta^\dagger\to a^\dagger e^{+i\omega_z t}. Why: transforming into the frame rotating at ω0\omega_0 and ωz\omega_z separates the fast optical phase from the slow motional dynamics, so that resonance conditions become visible as time-independent terms.

Step 2. Expand eiη(a+a†)≈1+iη ⁣(a e−iωzt+a†e+iωzt)e^{i\eta(a+a^\dagger)}\approx 1+i\eta\!\left(a\,e^{-i\omega_z t}+a^\dagger e^{+i\omega_z t}\right) to first order in η\eta. The three pieces now oscillate at distinct frequencies: the 11 term at δ\delta (carrier), the aa term at δ−ωz\delta-\omega_z (red), the a†a^\dagger term at δ+ωz\delta+\omega_z (blue). Why: each power of η\eta carries a definite number of motional quanta, so first order is precisely the Δn=±1\Delta n=\pm1 sidebands and nothing else — legitimate whenever ηn+1≪1\eta\sqrt{n+1}\ll1.

Step 3. Tune the laser to the red sideband, δ=−ωz\delta=-\omega_z. Then the aa term becomes time-independent while the carrier and blue terms rotate at ∓ωz\mp\omega_z and average away (rotating-wave approximation), leaving Hrsb=iℏΩη2 σ+a+h.c.H_{\mathrm{rsb}}=\tfrac{i\hbar\Omega\eta}{2}\,\sigma_+ a+\mathrm{h.c.} Why: choosing the detuning selects one and only one energy-conserving process, ∣g,n⟩→∣e,n−1⟩\lvert g,n\rangle\to\lvert e,n-1\rangle — this is the anti-Jaynes–Cummings coupling, the same interaction that powers the Cirac–Zoller gate.

Step 4. Evaluate the matrix element on Fock states: ⟨n−1∣a∣n⟩=n\langle n-1\rvert a\lvert n\rangle=\sqrt{n}, so Ωn,n−1=Ω η n\Omega_{n,n-1}=\Omega\,\eta\,\sqrt{n}, and likewise the blue sideband gives Ωn,n+1=Ωηn+1\Omega_{n,n+1}=\Omega\eta\sqrt{n+1}. Why: the bosonic n\sqrt{n} is decisive — it makes cooling slow down as nn shrinks and vanish outright at n=0n=0, which is exactly why the ground state is dark and cooling self-terminates there.

Go deeper · Derivation B · 1 minSteady-state phonon number from rate equations
Derivation BSteady-state phonon number from rate equations

The goal is to derive the resolved-sideband cooling floor nˉss≈(γ/2ωz)2\bar n_{\mathrm{ss}}\approx(\gamma/2\omega_z)^2 by treating the Fock-state populations as classical occupations evolving under scattering rates, and balancing the phonon-removing (red) channel against the phonon-adding (blue and off-resonant carrier) channels. This is the calculation that tells an experimentalist how cold the ion can get and how the answer improves with a narrower transition or a stiffer trap; it uses Fermi's golden rule with Lorentzian lineshapes and a detailed-balance steady state.

Step 1. Write the scattering rates into the Δn=∓1\Delta n=\mp1 channels as Lorentzians centred on the sidebands, A±∝η2Ω2 γ/2(δ±ωz)2+(γ/2)2A_\pm\propto\eta^2\Omega^2\,\dfrac{\gamma/2}{(\delta\pm\omega_z)^2+(\gamma/2)^2}, where A−A_- removes a phonon (cooling) and A+A_+ adds one (heating). Why: the excited state has a finite lifetime 1/γ1/\gamma, so its absorption profile is Lorentzian of width γ\gamma; the off-resonant tail of that Lorentzian is precisely the leak that limits cooling.

Step 2. Tune to the red sideband, δ=−ωz\delta=-\omega_z, so the cooling channel sits on resonance (denominator (γ/2)2(\gamma/2)^2) while the heating channel is detuned by 2ωz2\omega_z (denominator (2ωz)2+(γ/2)2(2\omega_z)^2+(\gamma/2)^2). Why: this is the whole point of resolving the sidebands — resonantly enhancing cooling while pushing the heating process far into the wing of the line.

Step 3. Take the ratio A+A−=(γ/2)2(2ωz)2+(γ/2)2≈(γ4ωz)2\dfrac{A_+}{A_-}=\dfrac{(\gamma/2)^2}{(2\omega_z)^2+(\gamma/2)^2}\approx\left(\dfrac{\gamma}{4\omega_z}\right)^2 in the resolved limit ωz≫γ\omega_z\gg\gamma. Why: the suppression of heating is set entirely by how far the blue sideband lies out on the Lorentzian tail, i.e. by the spectral mismatch 2ωz2\omega_z measured in linewidths.

Step 4. Impose detailed balance. In steady state the up- and down-flows through each rung match, giving nˉss=A+A−−A+≈A+A−\bar n_{\mathrm{ss}}=\dfrac{A_+}{A_--A_+}\approx\dfrac{A_+}{A_-}; including the recoil bookkeeping of spontaneous emission promotes the prefactor to the standard nˉss≈(γ/2ωz)2\bar n_{\mathrm{ss}}\approx(\gamma/2\omega_z)^2, which vanishes as ωz→∞\omega_z\to\infty. Why: a geometric-series steady state of a birth–death process is fixed by the single ratio A+/A−A_+/A_-, so the cooling floor is quadratic in γ/ωz\gamma/\omega_z — halving the linewidth quarters the residual motion.

Go deeper · Derivation C · 1 minSideband thermometry from the red-to-blue asymmetry
Derivation CSideband thermometry from the red-to-blue asymmetry

The goal is to derive the calibration-free thermometer nˉ=Prsb/(Pbsb−Prsb)\bar n=P_{\mathrm{rsb}}/(P_{\mathrm{bsb}}-P_{\mathrm{rsb}}) that lets an experiment read the mean phonon number straight off two pulse-excitation measurements, and to explain why it works only because the red and blue sidebands carry the asymmetric bosonic factors n\sqrt{n} and n+1\sqrt{n+1}. This is the diagnostic every ground-state-cooling experiment reports; the derivation weights the Fock-state couplings by a thermal distribution and compares the two sidebands in the weak-excitation limit.

Step 1. For a thermal motional state the Fock populations are pn=nˉn/(1+nˉ)n+1p_n=\bar n^n/(1+\bar n)^{n+1}. Under a short red-sideband pulse the excitation from rung nn is ∝Ωn,n−12=Ω2η2n\propto\Omega_{n,n-1}^2=\Omega^2\eta^2 n, and under a blue pulse ∝Ωn,n+12=Ω2η2(n+1)\propto\Omega_{n,n+1}^2=\Omega^2\eta^2(n+1). Why: in the weak-drive limit excitation probability is proportional to the squared coupling, so the whole nn-dependence reduces to the bosonic factors nn and n+1n+1.

Step 2. Average over the thermal distribution: Prsb∝∑npn n=nˉP_{\mathrm{rsb}}\propto\sum_n p_n\,n=\bar n and Pbsb∝∑npn (n+1)=nˉ+1P_{\mathrm{bsb}}\propto\sum_n p_n\,(n+1)=\bar n+1, using ⟨n⟩=nˉ\langle n\rangle=\bar n. Why: the mean of nn over a thermal state is by definition nˉ\bar n, so the two sideband signals differ by exactly the additive +1+1 of the blue channel — the spontaneous-versus-stimulated asymmetry of a boson.

Step 3. Form the ratio PrsbPbsb=nˉnˉ+1\dfrac{P_{\mathrm{rsb}}}{P_{\mathrm{bsb}}}=\dfrac{\bar n}{\bar n+1} and invert to get nˉ=PrsbPbsb−Prsb\bar n=\dfrac{P_{\mathrm{rsb}}}{P_{\mathrm{bsb}}-P_{\mathrm{rsb}}}; the common prefactor Ω2η2\Omega^2\eta^2 cancels, so no knowledge of Ω\Omega, η\eta, or pulse area is needed. Why: taking a ratio of the two sidebands divides out every experimental scale factor, which is what makes the measurement absolute — and as nˉ→0\bar n\to0 the red sideband disappears, the unmistakable signature of the ground state.

Limiting cases

The behaviour of sideband cooling is fixed by two dimensionless ratios — the sideband resolution γ/ωz\gamma/\omega_z and the Lamb–Dicke parameter η\eta — and every regime below is one corner of that plane:

CaseResultMeaning
Resolved sideband, deep Lamb–Dickeωz≫γ, η≪1\omega_z\gg\gamma,\ \eta\ll1Ideal regime: clean red-sideband cooling to nˉ≈(γ/2ωz)2≪1\bar n\approx(\gamma/2\omega_z)^2\ll1; ground state reliably reached
Unresolved sidebandsωz≲γ\omega_z\lesssim\gammaSatellites overlap the carrier; only Doppler cooling is possible, nˉ≈γ/2ωz∼\bar n\approx\gamma/2\omega_z\sim tens — used as the pre-cooling stage
Large Lamb–Dicke parameterη≳1\eta\gtrsim1Higher sidebands (±2,±3,…\pm2,\pm3,\dots) and recoil-driven Δn>1\Delta n>1 jumps activate; the simple two-state cooling picture fails, demanding a full optical-Bloch treatment
Strong drive (carrier-limited)Ωeff≳ωz\Omega_{\mathrm{eff}}\gtrsim\omega_zPower broadening off-resonantly excites the carrier and blue sideband, raising the floor; Ω\Omega is optimized to trade cooling rate against final nˉ\bar n
Stiff-trap / narrow-line limitωz/γ→∞\omega_z/\gamma\to\inftynˉss→0\bar n_{\mathrm{ss}}\to0: quadrupole or Raman transitions with tiny γ\gamma push the floor arbitrarily low, at the cost of a slow cooling rate

Worked examples

1 · Lamb–Dicke parameter and recoil for a calcium ion (confidence builder)

A single 40Ca+^{40}\mathrm{Ca}^+ ion is cooled on the narrow S1/2→D5/2S_{1/2}\to D_{5/2} quadrupole transition at λ=729 nm\lambda=729\ \mathrm{nm} in an axial trap of frequency ωz/2π=1.0 MHz\omega_z/2\pi=1.0\ \mathrm{MHz}. Find the zero-point extent z0z_0, the Lamb–Dicke parameter η\eta, and the recoil energy in units of a phonon.

Step 1. Zero-point spread. With m=40×1.6605×10−27=6.64×10−26 kgm=40\times1.6605\times10^{-27}=6.64\times10^{-26}\ \mathrm{kg} and ωz=2π×106=6.28×106 s−1\omega_z=2\pi\times10^6=6.28\times10^6\ \mathrm{s^{-1}},

z0=ℏ2mωz=1.055×10−342(6.64×10−26)(6.28×106)=1.12×10−8 m=11.2 nm.z_0=\sqrt{\frac{\hbar}{2m\omega_z}}=\sqrt{\frac{1.055\times10^{-34}}{2(6.64\times10^{-26})(6.28\times10^6)}}=\mathbf{1.12\times10^{-8}\ m}=11.2\ \mathrm{nm}.

The ion's quantum jitter is about 11 nm — some sixty times smaller than the 729 nm wavelength, the first hint that η\eta is small.

Step 2. Lamb–Dicke parameter. The wavevector is k=2π/λ=2π/(729×10−9)=8.62×106 m−1k=2\pi/\lambda=2\pi/(729\times10^{-9})=8.62\times10^6\ \mathrm{m^{-1}}, so

η=kz0=(8.62×106)(1.12×10−8)=0.097≈0.10.\eta=k z_0=(8.62\times10^6)(1.12\times10^{-8})=\mathbf{0.097}\approx0.10.

This is safely in the Lamb–Dicke regime η≪1\eta\ll1, so first-order sidebands dominate and the two-state cooling picture holds.

Step 3. Recoil in phonons. The recoil energy compared with one motional quantum is η2=(0.097)2=9.4×10−3\eta^2=(0.097)^2=\mathbf{9.4\times10^{-3}}. A scattered photon deposits under one percent of a phonon of recoil, so emission almost always leaves nn unchanged — exactly the condition equation (3) needs for the ratchet to work cleanly. z₀ = 11.2 nm, η ≈ 0.10, and recoil = 0.9% of a phonon — deep Lamb–Dicke.

2 · From Doppler pre-cooling to the ground state (cooling budget)

The same ion is first Doppler-cooled on the strong S1/2→P1/2S_{1/2}\to P_{1/2} dipole line at 397 nm397\ \mathrm{nm} (γ/2π=22.4 MHz\gamma/2\pi=22.4\ \mathrm{MHz}), then resolved-sideband cooled on the quadrupole transition quenched to an effective linewidth γeff/2π=20 kHz\gamma_{\mathrm{eff}}/2\pi=20\ \mathrm{kHz}, with ωz/2π=1.0 MHz\omega_z/2\pi=1.0\ \mathrm{MHz}. Estimate the phonon number after each stage and the resulting ground-state population.

Step 1. Doppler stage. On the broad dipole line γ≫ωz\gamma\gg\omega_z (22.4 MHz vs 1 MHz), so the sidebands are unresolved and only Doppler cooling operates, with floor

nˉDoppler≈γ2ωz=22.42(1.0)=11.2.\bar n_{\mathrm{Doppler}}\approx\frac{\gamma}{2\omega_z}=\frac{22.4}{2(1.0)}=\mathbf{11.2}.

The ion arrives at roughly eleven phonons — cold enough to sit inside the Lamb–Dicke regime, but nowhere near the ground state.

Step 2. Sideband stage. Switching to the quenched quadrupole transition makes the sidebands resolved, ωz/γeff=1000/20=50≫1\omega_z/\gamma_{\mathrm{eff}}=1000/20=50\gg1, so equation (5) gives the floor

nˉss≈(γeff2ωz)2=(202×1000)2=(0.01)2=1×10−4.\bar n_{\mathrm{ss}}\approx\left(\frac{\gamma_{\mathrm{eff}}}{2\omega_z}\right)^2=\left(\frac{20}{2\times1000}\right)^2=(0.01)^2=\mathbf{1\times10^{-4}}.

Four orders of magnitude below the Doppler floor, and well under the nˉ<0.05\bar n<0.05 benchmark for a good motional bus.

Step 3. Ground-state population. For a thermal residual state the ground fraction is P0=1/(1+nˉ)=1/(1+10−4)=0.9999P_0=1/(1+\bar n)=1/(1+10^{-4})=\mathbf{0.9999}. The ion spends 99.99% of its time in ∣n=0⟩\lvert n=0\rangle — precisely the cold, definite bus the Cirac–Zoller and Mølmer–Sørensen gates require. Doppler leaves n̄ ≈ 11; resolved-sideband cooling reaches n̄ ≈ 1×10⁻⁴, i.e. P₀ ≈ 99.99%.

3 · Reading the temperature off the sidebands (thermometry)

After cooling, a fixed-area probe pulse excites the red sideband with probability Prsb=0.05P_{\mathrm{rsb}}=0.05 and the blue sideband with Pbsb=0.55P_{\mathrm{bsb}}=0.55. What is the mean phonon number?

Step 1. Apply the calibration-free formula (7): nˉ=Prsb/(Pbsb−Prsb)\bar n=P_{\mathrm{rsb}}/(P_{\mathrm{bsb}}-P_{\mathrm{rsb}}). The unknown pulse area and η\eta cancel in the ratio, so the two numbers alone determine nˉ\bar n.

Step 2. Substitute: nˉ=0.05/(0.55−0.05)=0.05/0.50=0.10\bar n=0.05/(0.55-0.05)=0.05/0.50=\mathbf{0.10}. The strong red–blue asymmetry (a red sideband eleven times weaker than the blue) signals a nearly cold ion; had the ion been in the ground state the red sideband would have vanished entirely. n̄ = 0.10, read directly from the sideband asymmetry with no calibration.

Physics of the portrait

Sideband cooling is one man's idea carried, over four decades, from a two-page proposal to the beating heart of every trapped-ion quantum computer. David Wineland proposed laser cooling of bound ions before it had ever been done, demonstrated it, wrote down its limits, and led the group that first reached the motional ground state — work recognized with the 2012 Nobel Prize.

Origins · 1 minDavid Wineland and how this was discovered
DW
David Wineland
born 1944 · laser cooling and trapped-ion quantum control

David J. Wineland is the physicist who turned the trapped ion into the most precisely controlled object in physics. In 1975, with Hans Dehmelt, he proposed that light scattered by a bound ion could carry away its motional energy — laser cooling by driving the lower motional sideband — three years before anyone cooled an ion at all. In 1978 his group at NBS (now NIST) in Boulder demonstrated laser cooling of trapped Mg+\mathrm{Mg}^+, and the following year he and Wayne Itano laid out the theory of the cooling limits, including the resolved-sideband floor nˉ≈(γ/2ωz)2\bar n\approx(\gamma/2\omega_z)^2. The decisive experiment came in 1989, when Diedrich, Bergquist, Itano, and Wineland cooled a single Hg+\mathrm{Hg}^+ ion to its motional ground state for the first time, verifying it by the disappearance of the red sideband. When Cirac and Zoller proposed in 1995 that a shared, ground-state-cooled phonon mode could entangle ion qubits, Wineland's group implemented the first such quantum logic gate within months — welding sideband cooling permanently to quantum computing. His career, from the hydrogen maser through optical clocks to quantum information, is a sustained argument that if you can cool and control a particle's every quantum, you can measure and compute with it better than any other way. He shared the 2012 Nobel Prize in Physics with Serge Haroche for these methods of measuring and manipulating individual quantum systems.

History
1975
— Wineland and Dehmelt propose laser cooling of bound ions via the motional sideband
1978
— first laser cooling of trapped ions demonstrated at NBS Boulder and by Neuhauser–Dehmelt
1979
— Wineland and Itano work out the laser-cooling limits, including the resolved-sideband floor
1989
— Diedrich, Bergquist, Itano, Wineland cool a single Hg⁺ ion to its motional ground state
1995
— Monroe et al. achieve Raman sideband cooling of Be⁺ to the 3D ground state; Cirac–Zoller gate proposed

Going further

Texts and reviews. D. Leibfried, R. Blatt, C. Monroe, D. Wineland, "Quantum dynamics of single trapped ions", Rev. Mod. Phys. 75, 281 (2003) — the definitive review of the formalism used here · D. Wineland, W. Itano, "Laser cooling of atoms", Phys. Rev. A 20, 1521 (1979) — original cooling-limit theory · F. Diedrich, J. Bergquist, W. Itano, D. Wineland, "Laser cooling to the zero-point energy of motion", Phys. Rev. Lett. 62, 403 (1989) — first ground-state cooling · C. Monroe et al., "Resolved-sideband Raman cooling of a bound atom to the 3D zero-point energy", Phys. Rev. Lett. 75, 4011 (1995) · H. Häffner, C. Roos, R. Blatt, "Quantum computing with trapped ions", Phys. Rep. 469, 155 (2008) · Nielsen & Chuang, Quantum Computation and Quantum Information, Ch. 7 (physical implementations). The key relations in this article are: the Lamb–Dicke parameter η=kℏ/2mωz≈0.10\eta=k\sqrt{\hbar/2m\omega_z}\approx0.10 for 40Ca+^{40}\mathrm{Ca}^+ at 729 nm with ωz/2π=1\omega_z/2\pi=1 MHz; the sideband Rabi scaling Ωn,n−1=Ωηn\Omega_{n,n-1}=\Omega\eta\sqrt{n}; the resolved-sideband floor nˉss≈(γ/2ωz)2\bar n_{\mathrm{ss}}\approx(\gamma/2\omega_z)^2; the Doppler relation kBTD=ℏγ/2⇒nˉ≈γ/2ωzk_B T_D=\hbar\gamma/2\Rightarrow\bar n\approx\gamma/2\omega_z; and the thermometry formula nˉ=Prsb/(Pbsb−Prsb)\bar n=P_{\mathrm{rsb}}/(P_{\mathrm{bsb}}-P_{\mathrm{rsb}}). The cooling-budget numbers in Worked Example 2 use representative values (γ397/2π=22.4\gamma_{397}/2\pi=22.4 MHz, γeff/2π=20\gamma_{\mathrm{eff}}/2\pi=20 kHz) and are order-of-magnitude illustrations, not a specific published dataset.

Prerequisites. Quantum Harmonic Oscillator · Trapped-Ion Qubits · Two-Level Atom and Rabi Oscillations · Spontaneous Emission and Optical Bloch Equations. Leads to. Cirac–Zoller Gate · Mølmer–Sørensen Gate · DiVincenzo Criteria.

Ask the tutor: Ask me to walk through why the red-sideband coupling's n\sqrt{n} factor makes the motional ground state dark, or to work a cooling budget for your own ion species and trap frequency.

What is a phonon in a trapped ion, and why do we care about removing them?

A phonon here is one quantum of the ion's vibrational motion in the trap, energy ℏωz\hbar\omega_z; the ion's motion is a quantum harmonic oscillator, so its energy comes in these discrete lumps rather than continuously. We care because in a trapped-ion quantum computer the shared vibration of an ion chain is the bus that carries quantum information between qubits during an entangling gate. A gate maps the qubit state onto the motion and back, and if the motion is not in a single, definite state — the ground state ∣0⟩\lvert 0\rangle — the mapping is ambiguous and the gate is imperfect. Removing the phonons puts the bus in a clean, known starting state so the gate does exactly what it should.

Why does the red sideband cool while the blue sideband heats?

Tuning the laser one trap-quantum below the transition, to ω0−ωz\omega_0-\omega_z, means an absorbed photon is short of the energy needed to excite the ion. The ion makes up the deficit by borrowing a phonon from its own motion, so absorption removes one quantum of vibration: ∣g,n⟩→∣e,n−1⟩\lvert g,n\rangle\to\lvert e,n-1\rangle. The blue sideband at ω0+ωz\omega_0+\omega_z is the mirror image — the photon carries excess energy, which is dumped into the motion, adding a phonon. Since spontaneous emission on average returns the phonon number unchanged, driving red gives net cooling and driving blue gives net heating. Resolving the two sidebands is what lets you pick the cooling one and avoid the heating one.

What exactly is the Lamb–Dicke regime, and why must η be small?

The Lamb–Dicke regime is the condition that the ion's quantum position spread z0z_0 is much smaller than the laser wavelength, so the dimensionless parameter η=kz0≪1\eta=k z_0\ll1. Physically it means a photon's recoil kick is a tiny fraction of a phonon (η2≪1\eta^2\ll1), so absorbing or emitting a photon almost never changes the motional quantum number by more than one. That is what keeps the spectrum simple: only the carrier and the first red and blue sidebands matter, and the clean Δn=−1\Delta n=-1 cooling ratchet works. If η\eta approached one, higher sidebands and multi-phonon recoil events would proliferate, scrambling the motion faster than the red sideband could cool it.

How do you know when you have actually reached the ground state?

By the vanishing of the red sideband. In a thermal state the red-sideband signal is proportional to nˉ\bar n and the blue-sideband signal to nˉ+1\bar n+1; taking their ratio gives nˉ=Prsb/(Pbsb−Prsb)\bar n=P_{\mathrm{rsb}}/(P_{\mathrm{bsb}}-P_{\mathrm{rsb}}) with no calibration needed. As the ion approaches ∣0⟩\lvert 0\rangle the red sideband fades toward zero — because its coupling scales as n\sqrt{n} — while the blue sideband stays strong. A red-to-blue ratio near zero is the direct, model-free announcement that the ion is in its motional ground state. This is the measurement Diedrich and Wineland used in 1989.

Why must you Doppler-cool first instead of jumping straight to sideband cooling?

Because sideband cooling only behaves simply once the ion is already inside the Lamb–Dicke regime. A hot ion has a large motional amplitude, so ηn\eta\sqrt{n} is no longer small, higher-order sidebands are strong, and the tidy first-order picture fails. Doppler cooling, using a broad dipole transition, is fast and robust and brings the ion from room-temperature loading down to nˉ\bar n of order ten — cold enough that ηn≪1\eta\sqrt{n}\ll1. Only then does resolved-sideband cooling on a narrow transition take over to remove the last handful of phonons. The two stages are complementary: broad-and-fast first, narrow-and-deep second.

Does the same idea really work for macroscopic mechanical oscillators?

Yes — that is one of the most striking generalizations in the field. Cavity optomechanics couples a laser-driven optical cavity to a mechanical oscillator (a suspended mirror, a membrane, a nanobeam) whose motion shifts the cavity frequency, producing exactly the same red and blue motional sidebands. Driving the red sideband of the cavity scatters photons that carry away mechanical quanta, cooling the oscillator toward its ground state, with the same resolved-sideband condition (cavity linewidth much less than mechanical frequency) and the same residual floor set by the sideband ratio. Groups have cooled oscillators from billions of phonons to below one this way. The physics is identical; only the mass and frequency scales differ.

Is the cooling floor n̄ss ≈ (γ/2ω_z)² a fundamental limit?

It is the fundamental limit for a given cooling transition in the resolved-sideband, weak-drive regime — the irreducible floor from off-resonant scattering into the heating sideband. But it is not a fixed number: you lower it by choosing a narrower transition (smaller γ\gamma) or a stiffer trap (larger ωz\omega_z), and it drops quadratically. This is why experiments use narrow quadrupole lines or two-photon Raman transitions, whose effective linewidth can be made very small, and push trap frequencies into the megahertz range. In practice other effects — anomalous trap heating from noisy electrodes, laser intensity noise, and off-resonant carrier excitation at strong drive — often set the achieved nˉ\bar n before this theoretical floor does.

Before this

Where this leads