Laser cooling of trapped ions

Ion traps and quantum information

Thus far, we have focused our attention on neutral atoms, and their interactions with the electromagnetic field. Charged atoms are an important physical system which are also important to study, for a variety of reasons. They can be confined in with much deeper potential wells, using electric fields, than is possible with light forces alone. Ions can also be readily trapped in small numbers, allowing single atoms or small, discrete numbers of atoms to be individually manipulated. This capability has made trapped ions an excellent physical system for creation of exotic quantum states, such as Schr\"odinger cat superpositions, and highly entangled multi-atom states. Trapped ions have also proven to be an excellent platform for implementation of many quantum information protocols, including teleportation, superdense coding, and quantum error correction; simple quantum algorithms such as the Deutsch-Jozsa and Grover quantum search algorithms have also been demonstrated. Efforts are underway in the community to realize trapped ion systems with not just handfuls, but rather, hundreds to thousands of individually controlled ion qubits, in highly integrated trapped ion chips. This chapter begins with a description of laser cooling in the trapped ion physical system. We then describe how arbitrary motional states of ions can be engineered, in the second section. The third section describes quantum computation with trapped ions, and the fourth section covers concepts in quantum simulation with ions.

Laser cooling of trapped ions

A simple comparison of potential depths provides useful perspective for understanding the experimental convenience of trapped ions versus neutral atoms. Recall that the depth of a dipole force trap is $\hbar \Omega ^{2}/4\delta$; for a $1$ Watt laser at $1.06$ $\mu$m focused to a $10$ $\mu$m waist used to trap atomic sodium, the trap depth is $260$ $\mu$K; this is not much more than the Doppler cooling limit. A magneto-optical trap is much deeper; for typical laboratory magnetic field gradients, depths of $\sim 10^{-2}$ K can be obtained. In contrast, a singly charged ion in an electric field of $\sim 100$ V/mm, which is easily obtained in the laboratory, the depth of an RF ion trap can be over $\sim 10^{4}$ K. This allows hot ions to be readily trapped, and subsequently laser cooled. A unique aspect of trapped ions, versus neutral atoms, is the importance of the role of the trapping potential in the laser cooling process. As the motion of ions is reduced to the characteristic size of their confining potential, quantum aspects of their motion become accessible. This allows a form of laser cooling known as {\em resolved sideband cooling} to be performed, through which ions can be cooled to their motional ground states.

Physical system

Two kinds of ion traps are widely used. The Penning trap employs static magnetic and electric fields to confine ions. The Paul trap employs just oscillating electric fields, and is the configuration we shall study here. Consider this electrode configuration:

\noindent The inner surfaces are hyperboloids, the top and bottom electrodes are tied together, and the drive voltage oscillates sinusoidally at frequency $\Omega$, such that near the center of the electrodes, the potential is

$V(t)={\frac {V_{0}\cos \Omega t}{2r_{0}^{2}}}(x^{2}+y^{2}-2z^{2})\,.$

For a singly charged ion of mass $M$ located in this trap, the solutions to the equations of motion balancing the force of the electric field against ion motion are

${\ddot {x_{k}}}+2q_{k}{\frac {\Omega ^{2}}{4}}x_{k}\cos \Omega t=0\,.$

This is a Mathieu equation, where

$q_{x}=q_{y}={\frac {2eV_{0}}{Mr_{0}^{2}\Omega ^{2}}}$

and $q_{z}=-2q_{x}$ is a non-dimensional parameter. When $q\leq 0.9$, stable solutions exist. These solutions describe a fast oscillation known as micromotion, superposed on top of a slow harmonic motion with secular frequency

$\nu _{x}={\frac {q_{x}\Omega }{2{\sqrt {2}}}}\,.$

The amplitude of micromotion, in a well designed experiment, is very small, and thus can be neglected. The dominant physics of a trapped ion system are thus described by an atom coupled to a simple harmonic oscillator. Typically, RF frequencies are $\Omega \sim 10$ MHz or higher, and secular frequencies are $\nu \sim 1$ MHz. These may be compared with typical spontaneous emission rates of $\Gamma \sim 1$ MHz, and the recoil energy from a single photon emission of

$\omega _{R}={\frac {\hbar k^{2}}{2M}}\sim 100~{\rm {kHz}}\,,$

where $\omega _{L}=ck$ is the frequency of the laser. The three systems involved in this scenario are thus the trapped ion, with its atomic levels, the harmonic oscillator, and incident electromagnetic radiation on the atom. That is, we have a system of atom + photon + phonon, where the phonon represents quantized vibrational modes of the trapped atom.

Hamiltonian: classical field + 2-level atom

Let us model the trapped ion system initially using a classical field, and neglect spontaneous emission. Let us also approximate the atom to be a two-level system. The Hamiltonian for such a system has the interaction

${\begin{array}{rcl}H_{I}&=&{\vec {d}}\cdot {\vec {E}}\\&=&\left[{d{\vec {S}}}\right]\cdot \left[{E{\hat {x}}\cos(kz-\omega t)}\right]\,,\end{array}}$

where we use ${\vec {S}}$ to denote the spin of the atom, and assume the electric field is along the ${\hat {x}}$ direction. $z$ is the position of the ion in the field; we quantize this degree of freedom by representing it as motion in a harmonic oscillator, letting

$z=z_{0}(a+a^{\dagger })\,,$

where $z_{0}$ is the characteristic length scale of the harmonic motion. The interaction Hamiltonian can thus be written as

$H_{I}=dES_{x}\cos \left[{kz_{0}(a+a^{\dagger })-\omega t}\right]\,.$

The Lamb-Dicke parameter

When the ions are well-confined, a natural small parameter arises, in which this interaction can be expanded. This parameter describes the extent to which the ion is localized in the trap, relative to the incident light. We define this Lamb-Dicke parameter as

$\eta =kz_{0}={\frac {2\pi z_{0}}{\lambda }}={\frac {\omega _{R}}{\nu }}\,.$

It is the ratio of the size of the ground state wavefunction of the motion in the harmonic oscillator, to the incident laser wavelength:

\noindent Note that $\eta$ can also be understood as being the ratio of the recoil frequency to the vibrational frequency.

Expansion in the Lamb-Dicke parameter

Let us now expand $H_{I}$ to leading order in $\eta$. Recall that $S_{x}=(S_{+}+S_{-})/2$, and let $\Omega =(d\times E)/2\hbar$ be the Rabi frequency. In terms of these,

${\begin{array}{rcl}H_{I}&=&\hbar \Omega \left[{S_{+}+S_{-}}\right]\cos \left[{\eta (a+a^{\dagger })-\omega t}\right]\\&=&{\frac {\hbar \Omega }{2}}\left[{S_{+}+S_{-}}\right]\left[{e^{i\eta (a+a^{\dagger })-\omega t}+e^{-i\eta (a+a^{\dagger })+\omega t}}\right]\,.\end{array}}$

The exponentials can be expanded to leading order in $\eta$, resulting in terms of the form $\left[{1+i\eta (a+a^{\dagger })}\right]e^{-i\omega t}$, such that

${\begin{array}{rcl}H_{I}&=&{\frac {\hbar \Omega }{2}}\left[{S_{+}e^{-i\omega t}+S_{-}e^{i\omega t}}\right]+{\frac {i\eta \hbar \Omega }{2}}(S_{+}+S_{-})\left[{a+a^{\dagger }}\right]\left({e^{-i\omega t}-e^{i\omega t}}\right)\,.\end{array}}$

In the above expansion, it is assumed that the laser's frequency is close to that of the atomic transition, such that to good approximation, terms oscillating as $\sim 2\omega$ can be dropped (the rotating wave approximation). The first term in this expression describes the carrier transition, in which the light changes the internal atomic state, and the second term describes sideband transitions, in which the light changes both the internal atomic state as well as its motional state. These transitions have a frequency spectrum which looks like:

Moving into the rotating frame defined by the atom's internal and motional states simplifies this Hamiltonian, providing a time-independent form. Let this frame be defined by the system Hamiltonian

$H_{0}=\hbar \omega _{0}S_{z}+\hbar \nu a^{\dagger }a\,,$

such that in the rotating frame, operators are transformed according to

${\begin{array}{rcl}S_{\pm }&\rightarrow &S_{\pm }e^{\pm i\omega _{0}t}\\a&\rightarrow &ae^{-i\nu t}\\a^{\dagger }&\rightarrow &a^{\dagger }e^{i\nu t}\,.\end{array}}$

The interaction Hamiltonian in this frame is approximated by

${\begin{array}{rcl}H'_{I}&=&e^{iH_{0}t/\hbar }H_{I}e^{-iH_{0}t/\hbar }\\&=&\left\{{\begin{array}{lr}{\frac {\hbar \Omega }{2}}\left[{S_{+}+S_{-}}\right]&{\omega \approx \omega _{0}~~{\rm {carrier}}}\\{\frac {i\eta \hbar \Omega }{2}}(S_{+}a^{\dagger }-S_{-}a)&{\rm {\omega \approx \omega _{0}+\nu ~~{\rm {blue~sideband}}}}\\{\frac {i\eta \hbar \Omega }{2}}(S_{+}a-S_{-}a^{\dagger })&{\rm {\omega \approx \omega _{0}-\nu ~~{\rm {red~sideband}}}}\end{array}}\right.\,.\end{array}}$

The first few energy levels are diagrammed here:

Resolved sideband cooling

When the ion is hot, doppler shifts will dominate, and the sidebands will be unresolvable. However, for many ions, standard doppler cooling is sufficient to reach the point at which the motional sidebands become resolvable. In that regime, a different cooling method can be applied; this method is known as {\em resolved sideband cooling}. The basic idea of this method is analogous to optical pumping: the laser is detuned to cause transitions on the red sideband of the ion, removing one quantum of motion for each photon absorbed. Just as for laser cooling of neutral atoms, let us calculate the cooling limit of such a procedure. Let $|n{\rangle }$ be an eigenstate of the harmonic motion of the ion, with motional quantum number $n$, and recall that $a|n{\rangle }={\sqrt {n}}|n-1{\rangle }$ and $a^{\dagger }|n{\rangle }={\sqrt {n+1}}|n+1{\rangle }$. The transition amplitude between motional states $|n{\rangle }$ and $|n\pm 1{\rangle }$ is given by this matrix element of the interaction Hamiltonian:

${\begin{array}{rcl}\langle n\pm 1|H'_{I}|n\rangle \approx \langle n\pm 1|i\eta {\frac {\hbar \Omega }{2}}\left[{{\sqrt {n}}|n-1\rangle \langle n|-{\sqrt {n+1}}|n+1\rangle \langle n|}\right]|n{\rangle }\,,\end{array}}$

in the limit that $\eta$ is small, and we're interested only in terms which couple $n$ and $n+1$. Let us define

$W(\Delta )={\frac {1}{\left({\frac {2\Delta }{\Gamma }}\right)^{2}+1}}$

as the lineshape function; it captures the frequency dependent response of the atom to a laser detuned by $\Delta =\omega _{L}-\omega _{0}$ from the atomic resonance. The incident laser light can connect three basic transitions: the carrier, the blue sideband, and the red sideband. How do those transitions contribute to changing the motional quantum number $n$? Excitations of the carrier transition couple $|g,n{\rangle }$ (using $|{\rm {atom}},{\rm {motion}}{\rangle }$) to the state $|e,n{\rangle }$, with rate proportional to $W(\Delta )$. The excited state can then decay to $|g,n-1{\rangle }$ or $|g,n+1{\rangle }$, changing the motional quantum number; these transitions occur with rates proportional to $\Gamma \eta ^{2}n\alpha$ and $\Gamma \eta ^{2}(n+1)\alpha$, respectively, where $\alpha$ is a geometric factor describing the probability for spontaneous emission to change $n$. For dipole emission into free space, $\alpha \approx 2/5$. These transitions are diagrammed as follows:

On the blue and red sidebands, absorption happens at a rate proportional to $n\eta ^{2}W(\Delta +\nu )$ and $(n+1)\eta ^{2}W(\Delta -\nu )$, while spontaneous emission happens at rate proportional to $\Gamma$, all-together connecting $|g,n{\rangle }$ with $|g,n\pm 1{\rangle }$. These transitions are diagrammed as follows:

To obtain the cooling limits of resolved sideband cooling, we need to write down the rate equations for these phonon-number changing carrier and sideband excitation processes, and solve for the steady state average phonon number. For the transitions

\noindent where $P_{n}$ is the probability of being in the $n$ phonon state, the transition rate coefficients are

$A_{\pm }={\frac {\Omega ^{2}}{\Gamma }}\eta ^{2}\left[{W(\Delta \mp \nu )+\alpha W(\Delta )}\right]\,.$

The first term in the brackets comes from the sideband excitations, and the second from the carrier. The rate equations for the populations are thus

${\dot {P}}_{n}=\left[{nP_{n-1}-(n+1)P_{n}}\right]A_{+}+\left[{(n+1)P_{n+1}-nP_{n}}\right]A_{-}\,.$

The average phonon number ${\bar {n}}$ is

${\bar {n}}=\sum _{n=0}^{\infty }nP_{n}\,,$

which evolves as

${\begin{array}{rcl}{\frac {d}{dt}}{\bar {n}}&=&A_{-}\left[{\sum _{n}n(n+1)P_{n+1}-n^{2}P_{n}}\right]+A_{+}\left[{\sum _{n}n^{2}P_{n-1}-n(n+1)P_{n}}\right]\\&=&-{\bar {n}}A_{-}+({\bar {n}}+1)A_{+}\,.\end{array}}$

This differential equation has a solution in which ${\bar {n}}(t)$ decays exponentially as $e^{(A_{-}-A_{+})t}$ to an equilibrium average phonon number ${\bar {n}}_{eq}$,

${\bar {n}}_{eq}={\frac {A_{+}}{A_{-}-A_{+}}}\,.$

The corresponding equilibrium temperature is given by

${\frac {P_{n+1}}{P_{n}}}=e^{-\hbar \nu /k_{B}T}={\frac {A_{+}}{A_{-}}}\,.$

It is insightful to evaluate these expressions in two limits. When the ion is weakly confined, such that $\Gamma \gg \nu$, then the sidebands are poorly resolved:

\noindent The optimum detuning is $\Delta =\omega _{L}-\omega _{0}=-\Gamma /2$, and cooling corresponds to standard Doppler cooling of free particles, giving

${\bar {n}}={\frac {\Gamma }{2\nu }}\,,$

with a final temperature of around

$k_{B}T=\hbar \nu \left({{\bar {n}}+{\frac {1}{2}}}\right)\approx {\frac {\hbar \Gamma }{2}}\,.$

Including geometric factors reduces this slightly, by $(1+\alpha )/2$. The strong confinement limit is reached when $\Gamma \ll \nu$, so that individual sidebands are well resolved:

\noindent This is the case of resolved sideband cooling, as mentioned above. For this case, the laser is detuned to $\omega _{L}=\omega _{0}-\nu$, such that $A_{+}\sim \alpha W(-\nu )$, and $A_{-}\sim W(0)=1$, such that the equilibrium average phonon number, given by

${\bar {n}}={\frac {1}{{\frac {A_{-}}{A_{+}}}-1}}\,,$

is approximately ${\bar {n}}\approx (\alpha \Gamma /2\nu )^{2}$, which is much less than one. This corresponds to a temperature of

$k_{B}T={\frac {\hbar \nu }{-\log {\frac {A_{+}}{A_{-}}}}}={\frac {\hbar \nu }{2\log \left({\frac {2\nu }{\alpha \Gamma }}\right)}}\,.$

Essentially, the final temperature is set by zero-point motion of the ion in the harmonic trap.