# Fluorescence during Doppler cooling of a single trapped atom

###### Abstract

We investigate the temporal dynamics of Doppler cooling of an initially hot single trapped atom in the weak binding regime using a semiclassical approach. We develop an analytical model for the simplest case of a single vibrational mode for a harmonic trap, and show how this model allows us to estimate the initial energy of the trapped particle by observing the fluorescence rate during the cooling process. The experimental implementation of this temperature measurement provides a way to measure atom heating rates by observing the temperature rise in the absence of cooling. This method is technically relatively simple compared to conventional sideband detection methods, and the two methods are in reasonable agreement. We also discuss the effects of RF micromotion, relevant for a trapped atomic ion, and the effect of coupling between the vibrational modes on the cooling dynamics.

###### pacs:

32.80.Lg,32.80.Pj,42.50.VkPresent address: ]Lockheed Martin, Huntsville, AL Present address: ]Weizmann Institute of Science, Rehovot, Israel

Laser cooling of trapped neutral atoms and atomic ions is a well established technique: for example, cooling to the motional ground state Diedrich et al. (1989); Monroe et al. (1995); Perrin et al. (1998) and motional state tomography Leibfried et al. (1998) are routinely performed with resolved motional-sideband excitation techniques. Sideband techniques require the natural linewidth of the cooling transition to be small compared to the vibrational frequency of the trapped particle, in order to allow the motional sidebands to be resolved. Many experiments are, however, conducted in the “weak-binding regime”, where is larger than the oscillation frequency. Here, the cooling process is essentially the same as Doppler cooling of free atoms, because the spontaneous decay process is short compared to the atom’s oscillation period Wineland and Itano (1979). Even in experimental setups that implement sideband techniques, an initial stage of such “Doppler cooling” is often employed. The first examinations of Doppler cooling of trapped ions Wineland et al. (1978); Neuhauser et al. (1978); Wineland and Itano (1979); Javanainen and Stenholm (1980, 1981a, 1981b); Itano and Wineland (1982) did not take into account the effects of micromotion due to the trapping RF field. After cooling and heating effects related to micromotion were observed, these effects were explained theoretically Blümel et al. (1989); DeVoe et al. (1989); Cirac et al. (1994) by including the effects of micromotion.

Here, we consider Doppler cooling of a single trapped atom or ion. While most previous work has focused on the final stages of cooling, our focus will be on the temporal dynamics of the cooling, particularly in the “hot regime” where the Doppler shift due to atom motion is comparable to or much larger than . For the 1-D case we find that the cooling rate can be calculated analytically in the weak-binding regime without assuming the atom to be in the Lamb-Dicke regime. For a trapped ion, when we take RF micromotion into consideration, stable, highly excited states emerge when only one mode is considered Peik et al. (1999). When all three vibrational modes of the ion are considered, we find that couplings between the modes tend to break the stability of such points allowing cooling to reach the Doppler limit.

A practical application of our results is to estimate the initial motional energy of an atom or ion from observations of the time dependence of the fluorescence during the cooling process. As mentioned above, sideband spectroscopy is the conventional technique for characterizing motional states, and it has been used to characterize the heating rate of ions in the absence of cooling Diedrich et al. (1989); Monroe et al. (1995); Turchette et al. (2000); Seidelin et al. (2006); Deslauriers et al. (2006); Pearson et al. (2006); Epstein et al. (2007). However, it is more complicated to implement experimentally than Doppler cooling, requiring more laser beams. Currently, considerable effort is being devoted to understanding the anomalous heating observed in ion traps Turchette et al. (2000); Seidelin et al. (2006); Deslauriers et al. (2006); Pearson et al. (2006); Epstein et al. (2007). A less complicated technique for measuring temperature could simplify this work.

This paper is structured as follows: In Sec. I we present a semiclassical model of the Doppler cooling process for a bound atom in the weak-binding regime. In Secs. II and III we analyze the fluorescence predicted by the model for a single vibrational mode unaffected by micromotion. Here, we consider single cooling trajectories and average over these with a given distribution of initial motional energies. We derive expressions useful for estimating initial temperature from fluorescence observations in these sections. Sec. IV discusses how to minimize the total measurement time required to estimate the mean initial energy. In Sec. V we consider the effects of other motional modes with and without taking into account any RF micromotion experienced by such modes. Sec. VI suggests modifications to the basic experimental protocol that might provide improved sensitivity of the temperature measurements. Sec. VII concludes the paper.

## I Model

We consider a semiclassical model of Doppler cooling of a single weakly trapped atom Wineland and Itano (1979); Itano and Wineland (1982). We will initially consider only a single mode of motion, taken to be along the direction. We assume a harmonic potential with oscillation frequency . In Sec. V we consider a more detailed model that includes three dimensions and micromotion for ions.

The atom is Doppler-cooled by a single laser beam of angular frequency and wave-vector , detuned by from the resonance frequency of a two-level, or “cycling”, transition between two internal states, and , of the atom. We write the coupling Hamiltonian as

(1) |

where is the atom position, is Planck’s constant, and is the resonant Rabi frequency.

We assume the atom is weakly bound in the direction, that is, is much less than the excited state decay rate . The atom’s level populations are then approximately in steady state with respect to the instantaneous effective detuning, , including the Doppler shift , where and are the -components of the velocity and wave-vector. The excited state population is then Loudon (1973)

(2) |

Here is the saturation parameter, proportional to the cooling beam intensity, .

The excited state population is associated with the photon scattering rate by the relation . While the momentum kicks associated with photon emission are assumed to average to zero over many absorption-emission cycles, the absorbed photons will impart a velocity dependent momentum transfer due to the scattering that can be described by a velocity-dependent force

(3) |

where is the atom’s mass. This velocity-dependent force will in general change the motional energy of the atom. If the relative change in energy over a motional cycle is small, we can average the effect of over the oscillatory motion to find the evolution of :

(4) |

where the average is over one motional cycle. The average energy change per scattering event is .

In addition to , the atom will experience a stochastic force due to photon recoil that, assuming isotropic emission, will cause heating at a rate Itano and Wineland (1982)

(5) |

where is the recoil energy associated with the scattering. We will mostly ignore the effects of recoil heating in what follows since it will be important only near the cooling limit.

## Ii Analysis

We will now analyze the time-dependence of the atom fluorescence during the Doppler cooling process, as predicted by the model introduced above.

To simplify the algebra, we will scale energies by times half the power-broadened linewidth, and time by the resonant scattering rate:

As an example of typical values, we consider a trapped ion, where the cooling transition at has a natural linewidth of . At a detuning of with and , we find that , where is Boltzmann’s constant, and . The detuning and recoil parameters are and .

The maximal change in energy per scattering event at a given energy is . The energy at which the maximal Doppler shift, which in the scaled units is equal to , is equal to the power-broadened linewidth, , is of interest during the cooling process. For reference we note that this energy corresponds to

(7) |

For the typical experimental parameters considered above, is equal to or for .

For harmonic oscillations, the instantaneous Doppler shift is distributed according to the probability density

(8) |

where is the Dirac function. Since the average energy change per scattering event is , and the instantaneous scattering rate is , where , the rate of change of averaged over the secular oscillations by Eq. (4) takes the form:

(9) |

as illustrated in Fig. 1. We can evaluate the integral as detailed in Appendix A, to find that

(10a) | ||||

(10b) |

where . The asymptotic approximation (10b) corresponds to approximating by , which is reasonable in the “hot” regime, where the peaks of have small overlap with the Lorentz line profile.

The scattering rate averaged over the motion is analogous to Eq. (10) and is given by

(11) |

as illustrated in Fig. 1. In the limit of , we find , so that according to Eq. (10b) we have in this limit . This corresponds to each photon on average extracting an energy of . This can be understood by noting that in the limit of , is to a good approximation uniform over the Lorentz line profile, and so the value of averaged over the scattering events will be zero. Since each scattering event extracts an energy of , the average cooling per scattering event should indeed be .

The time-dependence of is formally found by integrating as given by Eq. (10). For the asymptotic approximation (10b) we find

(12) |

where is the energy at , as plotted in the lower part of Fig. 1. For the exact expression (10a), we must resort to numerical methods, although we do find analytically that the cooling rate is maximal for related to by

(13) |

which quantifies our previous observation that is a typical energy scale of the cooling process.

The behavior of is qualitatively different for being smaller or larger than a critical detuning, . For , has a maximum at

(14) |

For the example parameters listed below Eq. (6), corresponds to a detuning of . Closer to resonance, i.e., when , no maximum occurs, as illustrated in Fig. 2. The maximal scattering rate is reached when one of the peaks of the Doppler distribution (II) is in resonance with the cooling transition, as illustrated in the insets of Fig. 1. In the regime where a maximum exists, the maximal scattering rate is found to exceed the steady state scattering rate by a factor of

(15) |

We emphasize that the only approximations made above are the weak-binding approximation and the omission of recoil heating. In particular, the trapped particle is not assumed to be in the Lamb-Dicke regime. For the weak-binding regime, implies that the motion is well outside the Lamb-Dicke regime. To find the cooling rate predicted by (10) in the Lamb-Dicke limit, we note that to first order in we have . This corresponds to decreasing exponentially with . Except for the omission of recoil heating, the value of the decay time agrees with previous work that assumed the atom was in the Lamb-Dicke regime Wineland and Itano (1979); Wells and Cook (1990).

In the above analysis, we have ignored recoil heating as given by Eq. (5). In the limit of , the ratio of heating to cooling is seen to be , which is a small fraction for realistic parameters. For , the cooling is less efficient and the contribution from recoil becomes more significant, leading to a nonzero steady-state energy. Nevertheless, ignoring recoil heating is reasonable when considering only fluorescence, since the scattering rate has almost reached its steady-state value when the effect of recoil becomes important. We have omitted recoil in this analysis to make the only free parameter and simplify the discussion. Recoil can be included in calculations by combining Eqs. (5), (10a), and (II).

## Iii Thermal averaging

An application of the analysis presented above is to estimate the initial motional energy of a trapped atom from the fluorescence observed during the cooling process. Using this method, we can estimate the average rate of heating experienced by a trapped atom in the absence of cooling by first allowing the atom to heat up without cooling for a certain period and then observing the time dependence of the fluorescence as the atom is re-cooled. As discussed in Sec. II, we have for that the average cooling per scattering event is . The approximate total number of photons scattered during the cooling of an atom with initial motional energy can consequently be approximated by . For the example parameters given in Sec. II, this corresponds to photons for , corresponding to . With typical photon detection efficiencies of less than , very few photons are registered in a single experiment. We must therefore repeat many experimental cycles consisting of a heating period and a cooling period.

We now consider the form of the fluorescence signal when averaged over many such experimental cycles. Here, we will assume the heating is stochastic and take the distribution of the motional energies at the beginning of each cooling period to be the Maxwell-Boltzmann distribution with mean energy ,

(16) |

However, the results below hold for any form of .

The thermally averaged scattering rate is conveniently written in terms of the propagator, , of : Let denote the energy at time of an atom with initial energy . We can then write the thermally averaged scattering rate at time as

(17) |

This can be efficiently computed numerically by noting that , as detailed in Appendix B. Figure 3 shows the thermally averaged scattering rate for a few different parameters.

The fluorescence predicted by Eq. (17) has been found to be in good agreement with experimentally observed fluorescence. We show one experimental data set for comparison in Fig. 4; the experiments are more fully described in Epstein et al. (2007). Furthermore, the resulting estimated heating rates have been found to agree well with results obtained using the Raman sideband technique Seidelin et al. (2006); Epstein et al. (2007). This agreement may at first seem surprising, given that the two methods probe very different energy scales. For the measurements based on the Raman sideband technique the ion was only allowed to heat for a few milliseconds thereby gaining a few motional quanta while the measurement results presented in Fig. 4 are based on heating periods allowing the ion to gain many motional quanta. However, the results should agree if, as expected, the heating rate is constant over these energy scales.

## Iv Optimal experimental parameters

We now examine how the total measurement time required to reach a given accuracy on the heating rate estimate depends on the choice of experimental parameters. As the recoil parameter will be fixed by choice of atom, we consider only the choice of optimal values for , , and laser beam intensity.

It is clear from Fig. 2 that the relevant size of the signal, in terms of fluorescence photons emitted, for a given initial motional energy increases with decreasing detuning: the re-cooling is slower and the change in scattering rate is larger. For a given experimental setup, the optimal detuning is decided as a compromise between re-cooling signal and ability to re-cool atoms that have been highly excited by e.g., collisions.

For a given value of , the experimental signal, in terms of the number of photons scattered before steady state is reached, does not depend on the laser beam intensity. Since is the average initial energy relative to , which is proportional to the power-broadened linewidth, a lower laser beam intensity will give a larger signal for a given heating period. This suggests using the smallest feasible laser intensity, requiring a compromise with respect to robust cooling and detector dark counts. From this standpoint we want to keep the saturation parameter below, but probably close to, .

For a given detuning and laser intensity, an additional choice of the length of the heating period in each experimental cycle has to be made: Should we perform a relatively low number of cycles with long heating periods or more cycles with shorter heating periods?

To answer this question, we estimate the total measurement time, , required to reach a certain relative accuracy on the estimate of the heating rate. We assume a constant heating rate and assume that the total time is dominated by the heating periods, so that is proportional to the average initial energy, , and to the the number of runs.

We will consider a setup where the observed fluorescence is collected in sequential time-bins that are short compared to the total time required for the cooling process. In the limit where the distribution of the integrated number of counts, , in time-bin is described by a normal distribution with variance , we can estimate the uncertainty on the maximum-likelihood estimate of for a given dataset by Press and William (1989)

(18) |

It follows from Eqs. (10a), (II), and (17), that in the 1-D case the cooling dynamics can be rewritten in a form independent of by reparametrizing in terms of , , and . We will denote the reparametrized scattering rate by

(19) |

Since the relative uncertainty on the heating rate estimate is equal to and , Eqs. (18) and (19) allow us to estimate the time required to obtain a given relative uncertainty on the heating rate:

(20) |

Note that the right hand side depends only on and . Fig. 5 shows calculated for different detunings. The figure confirms that a low detuning is indeed favorable, and also shows that for a given detuning, decreases with increasing . This is not surprising, given that the time to cool by a certain amount of energy increases with atom temperature, as illustrated by Fig. 1. It is clear from Fig. 5 that the heating period should be chosen long enough to get a significant signal, i.e., , but the optimal heating period must be decided based on other experimental parameters such as trap depth and background gas collision rate.

## V Cooling in three dimensions

So far, we have considered only cooling in one dimension. In this section we will consider the effect of the vibrational modes in other directions on the cooling process. Our goal is to gain a qualitative understanding of the effects of the transverse modes on the cooling dynamics of the mode, with the intent of establishing to what extent the simple 1-D model presented above is a reasonable approximation.

### v.1 3-D cooling of neutral atoms

For a neutral atom, the confinement transverse to is not associated with micromotion, as it is for ions, and the 1-D weak-binding model extends immediately to three dimensions. Let , , denote the motional energy in mode , the Doppler shift, and the maximum Doppler shift. Although all modes are formally identical in the absence of micromotion, we will discuss the cooling dynamics with a focus on the mode.

In experiments it is typically easy to make the frequencies of the three modes incommensurate, which we will assume here. In that case, we can write the rate of change of as

(21) |

where is the effective line profile experienced by the mode, obtained by convolving the Lorentz line profile with the distribution of the combined Doppler shift due to the and “spectator” modes,

(22) |

As illustrated in Figs. 6 and 7, is peaked (diverges) at . If , the peaks are separated by more than the width of the Lorentz profile, and will be double-peaked, as illustrated in Fig. 7. It follows from Eq. (V.1) that the cooling rate in the limit of small is proportional to the slope of at , and that the rate of change of is positive if the slope is negative. If , this will result in heating of the mode, at least as long as are both inside the peaks of , i.e., while , as illustrated in Fig. 7 for and DeVoe et al. (1989); Peik et al. (1999). The figure also shows that this thermalization or energy equilibration effect is not present if , as is not double-peaked in this case. Mathematically, , , and are all conveniently expressed as convolution integrals of functions with known Fourier transforms.

The dashed lines in Fig. 8 show the cooling rates predicted by Eq. (V.1) for the case of only one excited spectator mode at different energies. When is large compared to , we see that the cooling rate is almost unaffected by the spectator mode. This can be understood by noting that in the limit where is uniform over the values of where is nonzero, the symmetry of implies that the average energy change per scattering event is , as also discussed in Sec. II. Since for , this implies that the temperature of the spectator modes will not affect the cooling rate in this limit. At lower values of , we generally see a decrease in the cooling rate in a gradual approach to the thermalization regime discussed above.

The consequences of thermalization/equilibration process are complex, when considering the full 3-D cooling problem. Consider for instance the case where only one mode is initially hot. According to the discussion above, this will result in heating of the two remaining modes, until the fastest heating mode has reached a value of similar to that of the initially hot mode. After this thermalization, the modes will be cooled simultaneously at a cooling rate significantly lower than the cooling rate for a single hot mode.

At this point, it is worth reconsidering the validity of our omission of recoil heating: The recoil heating rate as given by Eq. (5) is seen to have a maximum value of at the resonant scattering rate. It is clear from Fig. 8 that for typical values of on the order of , recoil is insignificant at high energies.

### v.2 3-D cooling of ions including the effects of micromotion

For an ion in a linear Paul trap, if we take the direction to be the axis, confinement in the transverse and directions is provided by the ponderomotive potential of an RF quadrupole field. The full 3-D cooling problem including micromotion on the transverse modes is very complex even in the Lamb-Dicke regime DeVoe et al. (1989); Walther (1993); Cirac et al. (1994). At low saturation, the effects of the micromotion caused by an RF field of frequency can be modeled by including micromotion sidebands in the line profile DeVoe et al. (1989): Micromotion with peak amplitude , where is the ion position averaged over one period of the RF field, can be described by the line profile,

(23) |

where is the effective detuning, is the micromotion modulation index, is the scaled RF frequency, and is the -th Bessel function.

In contrast to the situation in Ref. DeVoe et al. (1989), we are considering a case where changes during the secular motion. Since and depend on and , respectively, we parametrize the secular motion by the instantaneous phases, , where , and where is slowly varying and . Choosing the and axes so that the RF field is proportional to , we find that . In the limit where the transverse confinement is modified only weakly by static potentials, so that , we find in the pseudopotential approximation that , which we note to be independent of the secular frequencies. In this case we have

(24) |

where the integral is over in all dimensions. Note that since the modulation index depends only on the transverse components of the motion, the effect of excited transverse modes on the cooling of the mode can still be described in terms of an effective line profile, as in Eq. (V.1).

For the cooling of the transverse modes, the effects of micromotion on the cooling rates is pronounced, as illustrated by Fig. 8. A very clear qualitative difference from the cooling rate in the micromotion-free case is that at sufficiently high RF frequencies ( for ), stable points for the transverse mode energies develop even when the remaining modes are cold. This effect has been discussed in Ref. Peik et al. (1999) and is attributed to the heating peak of the Doppler distribution becoming resonant with a micromotion sideband, as described by Eq. (23). This might be related to the bistable behavior reported in some single ion experiments Sauter et al. (1988a, b); Walther (1993). The stability breaks down when thermalization is taken into consideration. Consider for instance the stable point indicated by Fig. 8 to exist for . Here, it is clear from the figure that when the mode has heated to , cooling of the mode will commence.

When for the transverse modes, we find that only the term of Eq. (23) contributes significantly, and the argument of Sec. V.1 that the cooling rate for the mode is not affected by excited transverse modes when also applies here, as illustrated by Fig. 8.

It is clear from the results above that we cannot ignore the transverse modes if their associated maximal Doppler shifts are comparable to that of the mode. If, however, we assume the transverse modes are cold enough to avoid the heating effects described in Figs. 7 and 8, we have seen above that the primary effect of the transverse modes will be to slow down the cooling of the mode. This would result in the 1-D model overestimating the mean initial energy of the mode. However, for many experiments that use linear RF traps, it is reasonable to assume that the transverse modes are heated significantly less than the mode. This is because most investigations of the anomalous heating in ion traps have found the results to be consistent with heating rates having a frequency dependence of with Turchette et al. (2000); Deslauriers et al. (2006); Epstein et al. (2007). Since the transverse mode frequencies are often an order of magnitude larger than , this would indeed lead to the transverse modes being significantly colder than the mode. Also, since the energy in the transverse modes only affects the cooling of the mode through the resulting Doppler shift, the effect of the transverse modes could be further reduced by aligning to have a smaller projection on the transverse modes. This would however reduce the efficiency of cooling of the transverse modes Itano and Wineland (1982).

Finally, another effect with respect to micromotion is that the presence of uncontrolled static stray fields can result in the ion experiencing micromotion even at the ion equilibrium position. At temperatures where , the first order effect according to Eq. (23) of this will be a reduction of the central spectral component by a factor of ; see for example Ref. Berkeland et al. (1998). We note that this effect can be compensated by using an effective saturation parameter based on the steady-state fluorescence observed in the trap.

### v.3 Departures from the weak-binding, low-saturation limit

In most experimental situations, we will not strictly fulfill the requirements of low saturation or weak binding. In particular, for the trap referenced in Fig. 4 the secular frequencies of the transverse modes are approximately equal to half the linewidth of the Doppler cooling transition, making the weak binding assumption only approximate. Also, the illustrated data were obtained at a saturation parameter of , outside the validity region of the line-profile model that accounts for RF micromotion (24). To validate our claim that the fluorescence signal predicted by the 1-D model is a good approximation if the heating rate is assumed to be a strongly decreasing function of , we performed a numerical Monte Carlo simulation of the fluorescence, based on integrating the optical Bloch equations through a large number of cooling trajectories. For each trajectory, we propagate the density matrix of the ion’s internal state according to the master equation

(25) |

where is the Lindblad operator for excited state decay and for the and introduced above. Coupling to the motional state is modeled by the average light force, . This model assumes neither the atoms to be weakly bound nor the cooling beam intensity to be low but does neglect recoil heating.

Figure 9 shows the result of fitting simulations with different assumptions for the frequency dependence of the heating to the dataset presented in Fig. 4. We find that if we assume the transverse modes are not heated, we obtain a temperature estimate of , in agreement with the result of fitting the 1-D model to the data, as illustrated by Fig. 4. The mode temperature of estimated from the model is close to, and slightly smaller than, this value, and agrees with the temperature estimate of based on extrapolating heating rates measured with the Raman sideband technique for the same trap configuration. This particular form of the frequency dependence of the heating rate was observed for the same trap when the Raman sideband technique Epstein et al. (2007) was used, and similar frequency dependencies have been observed in other geometries Turchette et al. (2000); Deslauriers et al. (2006). If we instead assume an dependence of the heating, the results only change slightly.

Our main conclusions from the simulation results are that the primary effect of the presence of weakly heated spectator modes will be to slow down cooling due to thermalization. If heating of the transverse modes is assumed, the 1-D model will somewhat overestimate the motional temperature of the axial mode.

## Vi Modified experimental protocols

We consider two modifications to the experimental protocol to reduce the total measurement time. Both are motivated by the fact that the size of the signal from a given amount of heating increases with increased initial energy.

One approach would be to coherently add a known amount of energy to the mode at the start of the heating period. If the added energy is enough to bring the atom into the slow-cooling regime, this will increase the signal change due to a given amount of additional heating, as illustrated in Fig. 10.

Alternatively, parametric amplification Dehmelt and Walls (1968); Caves (1982); Heinzen and Wineland (1990) could be employed after the heating cycle to modify the thermal distribution. Parametric amplification can be implemented by modulating the trap potential at , and leads to amplification of one quadrature of the motion while damping the other quadrature. For a low value of , parametric amplification would increase the fraction of experiments in which the atom is in the slow-cooling regime at the beginning of the cooling process, thus increasing the signal for a given heating period, as illustrated in Fig. 11.

## Vii Conclusion

In conclusion, we have shown that the motional energy of a trapped atom or ion can be estimated from the temporal changes in fluorescence observed when Doppler cooling is applied. Specifically, the initial energy can be estimated by fitting Eq. (17), where the mean initial motional energy is the only free variable, to the observed fluorescence. Our analysis assumes the oscillation frequency of the atoms is much smaller than the linewidth of the optical transition used for Doppler cooling and the motional energy at the start of the cooling is thermal.

Compared to Raman sideband transition methods for heating rate measurements, this method is simpler to implement experimentally but requires longer measurement duration for traps with low heating rates. On the other hand, for high heating rates, where sideband cooling is inefficient, this may be the method of choice. We have shown that in the typical situation, where the time for heating dominates, the total measurement time decreases with decreasing laser intensity, decreasing laser detuning, and increased heating period duration. We have compared the trade-off between these parameters (Fig. 5 and Eq. (20)). Finally, we show (Sec. IV) that the total measurement time can be reduced by adding additional energy to more quickly bring the ion into the low fluorescence regime.

By comparison with various models of 3-dimensional Doppler cooling, we have established that under typical experimental conditions the effects of the high-frequency modes are small, and that they will lead to temperature estimates that are somewhat higher than the actual temperature of the low-frequency mode.

Work supported by DTO and NIST. J.H.W. acknowledges support from The Danish Research Agency. R.J.E. acknowledges National Research Council Research Associateship Awards. S.S. acknowledges support from the Carlsberg Foundation. J.P.H. acknowledges support from a Lindemann Fellowship. We thank J. J. Bollinger and C. Ospelkaus for comments on the manuscript. This manuscript is a publication of NIST and is not subject to U.S. copyright.

## Appendix A Integrals

The integrals appearing in Eqs. (9), (II), and (V.1) are all convolution integrals of elements with analytical Fourier transforms and can thus be easily evaluated in Fourier space. For the 1-D integrals, the inverse Fourier transform can also be performed analytically. Here we present a more direct approach to evaluating the 1-D integrals.

For we define as

(26) |

where . Noting that

(27) |

we find that, according to Eq. (26),

Taking the branch cut discontinuity for to be along the negative real axis, we have for that , so that

(28) |

## Appendix B Numerical calculation of the averaged scattering rate

In this section we present an efficient numerical method for evaluating the averaged scattering rate given by Eq. (17).

Introducing , for , we note that . The values of are the energies along a single cooling trajectory. If the scattering rate can be considered constant on time scales of ,

(29) |

we find that the thermally averaged scattering rate, as given by Eq. (17), averaged over the same intervals can be approximated by

(30) |

Since the values of the are independent of , is easily calculated for different by list convolution.

For the Maxwell-Boltzmann distribution, a numerically stable form of the weight factors appearing in (30) is

## References

- Diedrich et al. (1989) F. Diedrich, J. C. Bergquist, W. M. Itano, and D. J. Wineland, Phys. Rev. Lett. 62, 403 (1989).
- Monroe et al. (1995) C. Monroe, D. M. Meekhof, B. E. King, S. R. Jefferts, W. M. Itano, D. J. Wineland, and P. Gould, Phys. Rev. Lett. 75, 4011 (1995).
- Perrin et al. (1998) H. Perrin, A. Kuhn, I. Bouchoule, and C. Salomon, Europhys. Lett. 42, 395 (1998).
- Leibfried et al. (1998) D. Leibfried, T. Pfau, and C. Monroe, Phys. Today 51, 22 (1998).
- Wineland and Itano (1979) D. J. Wineland and W. M. Itano, Phys. Rev. A 20, 1521 (1979).
- Wineland et al. (1978) D. J. Wineland, R. E. Drullinger, and F. L. Walls, Phys. Rev. Lett. 40, 1639 (1978).
- Neuhauser et al. (1978) W. Neuhauser, M. Hohenstatt, P. Toschek, and H. Dehmelt, Phys. Rev. Lett. 41, 233 (1978).
- Javanainen and Stenholm (1980) J. Javanainen and S. Stenholm, Appl. Phys. 21, 283 (1980).
- Javanainen and Stenholm (1981a) J. Javanainen and S. Stenholm, Appl. Phys. 24, 71 (1981a).
- Javanainen and Stenholm (1981b) J. Javanainen and S. Stenholm, Appl. Phys. 24, 151 (1981b).
- Itano and Wineland (1982) W. M. Itano and D. J. Wineland, Phys. Rev. A 25, 35 (1982).
- Blümel et al. (1989) R. Blümel, C. Kappler, W. Quint, and H. Walther, Phys. Rev. A 40, 808 (1989).
- DeVoe et al. (1989) R. G. DeVoe, J. Hoffnagle, and R. G. Brewer, Phys. Rev. A 39, 4362 (1989).
- Cirac et al. (1994) J. I. Cirac, L. J. Garay, R. Blatt, A. S. Parkins, and P. Zoller, Phys. Rev. A 49, 421 (1994).
- Peik et al. (1999) E. Peik, J. Abel, T. Becker, J. von Zanthier, and H. Walther, Phys. Rev. A 60, 439 (1999).
- Turchette et al. (2000) Q. A. Turchette, D. Kielpinski, B. E. King, D. Leibfried, D. M. Meekhof, C. J. Myatt, M. A. Rowe, C. A. Sackett, C. S. Wood, W. M. Itano, et al., Phys. Rev. A 61, 63418 (2000).
- Seidelin et al. (2006) S. Seidelin, J. Chiaverini, R. Reichle, J. J. Bollinger, D. Leibfried, J. Britton, J. H. Wesenberg, R. B. Blakestad, R. J. Epstein, D. B. Hume, et al., Phys. Rev. Lett. 96, 253003 (2006).
- Deslauriers et al. (2006) L. Deslauriers, S. Olmschenk, D. Stick, W. K. Hensinger, J. Sterk, and C. Monroe, Phys. Rev. Lett. 97, 103007 (2006).
- Pearson et al. (2006) C. E. Pearson, D. R. Leibrandt, W. S. Bakr, W. J. Mallard, K. R. Brown, and I. L. Chuang, Phys. Rev. A 73, 032307 (2006).
- Epstein et al. (2007) R. J. Epstein, D. Leibfried, S. Seidelin, J. H. Wesenberg, J. J. Bollinger, J. M. Amini, R. B. Blakestad, J. Britton, J. P. Home, W. M. Itano, et al., Simplified measurements of ion heating rates in a surface-electrode trap (2007), eprint arXiv:0707.xxxx.
- Loudon (1973) R. Loudon, The Quantum Theory of Light (Clarendon, Oxford, 1973).
- Wells and Cook (1990) A. L. Wells and R. J. Cook, Phys. Rev. A 41, 3916 (1990).
- Press and William (1989) I. Press and H. William, Numerical Recipes in Pascal (Cambridge University Press, 1989).
- Walther (1993) H. Walther, Adv. Atom. Mol. Opt. Phys. 31, 137 (1993).
- Sauter et al. (1988a) T. Sauter, H. Gilhaus, W. Neuhauser, R. Blatt, and P. E. Toschek, Europhys. Lett. 7, 317 (1988a).
- Sauter et al. (1988b) T. Sauter, H. Gilhaus, I. Siemers, R. Blatt, W. Neuhauser, and P. E. Toschek, Z. Phys. D 10, 153 (1988b).
- Berkeland et al. (1998) D. J. Berkeland, J. D. Miller, J. C. Bergquist, W. M. Itano, and D. J. Wineland, J. Appl. Phys. 83, 5025 (1998).
- Dehmelt and Walls (1968) H. G. Dehmelt and F. L. Walls, Phys. Rev. Lett. 21, 127 (1968).
- Caves (1982) C. M. Caves, Phys. Rev. D 26, 1817 (1982).
- Heinzen and Wineland (1990) D. J. Heinzen and D. J. Wineland, Phys. Rev. A 42, 2977 (1990).