Buffer gas cooling of ions in radio-frequency traps using ultracold atoms

Reaching ultracold temperatures within hybrid atom–ion systems is a major limiting factor for control and exploration of the atom–ion interaction in the quantum regime. In this work, we present results on numerical simulations of trapped ion buffer gas cooling using an ultracold atomic gas in a large number of experimentally realistic scenarios. We explore the suppression of micromotion-induced heating effects through optimization of trap parameters for various radio-frequency (rf) traps and rf driving schemes including linear and octupole traps, digital Paul traps, rotating traps and hybrid optical/rf traps. We find that very similar ion energies can be reached in all of them even when considering experimental imperfections that cause so-called excess micromotion. Moreover we look into a quantum description of the system and show that quantum mechanics cannot save the ion from micromotion-induced heating in an atom–ion collision. The results suggest that buffer gas cooling can be used to reach close to the ion’s groundstate of motion and is even competitive when compared to some sub-Doppler cooling techniques such as Sisyphus cooling. Thus, buffer gas cooling is a viable alternative for ions that are not amenable to laser cooling, a result that may be of interest for studies into cold controlled quantum chemistry and charged impurity physics.


I. INTRODUCTION
By immersing a trapped ion in a cold atomic bath, one creates a system that benefits both from the spatial localization and addressability of the ion, as well as from the scalability and ultracold temperatures of the atom cloud [1][2][3].The result is a unique combination of two well-controllable platforms linked together by the intermediate range ion-atom interaction, with applications in impurity physics [4][5][6][7][8][9][10][11][12], cold chemistry [13][14][15][16] and quantum technology [17,18].Buffer gas cooling of trapped ions has long been used to obtain cold samples [19][20][21][22].However, in efforts to reach quantized ion motion it has been outperformed by sub-Doppler laser cooling techniques, which have allowed small ion crystals to be cooled to their ground state of motion [23,24].Meanwhile, developments in laser and evaporative cooling techniques have opened up the possibility to prepare long-lived, ultracold atomic gases below about 1 µK, making them excellent buffer gases for cooling single ions.This progress has allowed buffer gas cooling to become a promising route towards ultracold atom-ion mixtures and exploring the realm of quantum effects within them.
Reaching the lowest possible temperatures in a hybrid atom-ion system is crucial to taking full advantage of the research opportunities it has to offer.In particular, an ion cooled to the regime of quantized motion would make an excellent probe to monitor interactions within the atomic bath and study impurity physics.Furthermore, cooling of both the atomic bath and the ion to atom-ion collision energies below the quantum limit would allow tuning of their interaction via s-wave Feshbach resonances [25,26], similar to those widely used in neutral atom systems [27] and recently observed in atom-ion systems [28].In previous work, we demonstrated that buffer gas cooling can be used to cool an atom-ion system to the edge of the s-wave collision energy, where the ion reached an average of 3.7 motional quanta [29].
However, efforts to reach deeper into the quantum regime or to buffer gas cool an ion to its motional ground state have remained elusive.
A major obstacle limiting buffer gas cooling of an ion in a radiofrequency (rf) trap is the heating effect of the time-dependent potential.In such a trap, energy stored in the ion's micromotion can be released during an atom-ion collision, heating the system [30].
This competes with the desired cooling effects of the buffer gas collisions.As a result, the temperature of the ion after thermalization is always higher than that of the buffer gas in which it is immersed.Reduction of atom-ion collision energies to the quantum regime are therefore predominantly limited by the kinetic energy of the ion, when trapping in a rf ion trap.For this reason, minimizing the ion's excess micromotion (EMM) through compensation of experimental imperfections is critical to buffer gas cooling experiments [31].
Furthermore, choosing an atom-ion species combination with a large mass ratio greatly suppresses the effects of micromotion-induced heating [30,32].Therefore, the coldest atomion collision energies and subsequent observation of quantum effects have been reported for Li/Yb + [29] and Li/Ba + [28] systems.Attempts to avoid these micromotion heating effects altogether by trapping ions in other types of trap, for example optical traps, are ongoing, yet offer different challenges [33].For an ion in a Paul trap, the question remains: what are the lowest temperatures achievable in a rf-trapped ion-atom system?
In this paper, we numerically simulate buffer gas cooling in time-dependent traps, showing that the optimization of trap parameters is crucial to reaching the lowest possible energies.
We focus on the Li/Yb + species combination, and find that parameter optimization can allow cooling to collision energies approximately two times smaller than what has previously been achieved experimentally [29].For these settings, we expect an occupancy of about one motional quanta on average, putting buffer gas cooling on par with sub-Doppler cooling techniques such as Sisyphus cooling [34].We study the collision energies that can be achieved by buffer gas cooling in alternative rf traps including hybrid rf-optical traps, digital Paul traps, rotating rf-traps and hexapole and octupole traps.We find that all investigated traps can be used to achieve the same low collision energies as in a linear Paul trap.The only investigated traps that offer a clear advantage over the linear Paul trap are the rf-optical hybrid trap and multipole traps, in which even lower temperatures can be reached.However these depend on the use of ambitious experimental parameters and, in the case of the hybrid trap, suppression of photon scattering effects which could otherwise lead to heating of the atom cloud.All simulations show that within these time-dependent traps, we can reach the energy regime where a quantum mechanical treatment of the atom-ion interaction becomes necessary.Therefore, we also discuss how to describe the atom-ion collision within a quantum description and show that the micromotion-induced heating survives.
The paper is organised as follows.In section II we outline the method with which we classically simulate the buffer gas cooling process, indicating how we extract and compare the ion kinetic energies and atom-ion collision energies reached in these simulations.In section III we detail the time-dependent potential used to confine an ion in a linear Paul trap, and show simulation results for the ion energies that can be reached in such a trap with variable parameters.In section IV, we consider other types of time-dependent ion traps including digital, rotating and multipole traps, comparing buffer gas cooling simulations in each of these to the results of the linear Paul trap.We make a similar comparison in section V, now focusing on a hybrid optical-rf trap.Finally, in section VI, we present a quantum description of the atom-ion interaction to support our results also in the regime of low kinetic energies of the ion.

II. METHOD
For our investigation into buffer gas cooling, we simulate repeated collisions between a trapped ion and atoms from an ultracold cloud.During these simulations, there are two effects that dictate the motion of the colliding particles as they are propagated towards each other in time.These are the atom-ion interaction V ai (r), and the time-dependent potential used to trap the ion V trap (t, r).The latter is adjusted to describe the rf trap we wish to investigate.Our simulations follow the method outlined in ref. [31], which has been shown to be in good agreement with experimental results [29].We adapt this process to allow for a broader choice of V trap (t, r), corresponding to a range of rf ion traps.We study the role played by V trap (t, r), and its defining parameters, in dictating the properties of the simulated ion after buffer gas cooling.
The atom-ion interaction consists of a long-range induced-dipole interaction, in which the charge of the ion polarizes the atom, giving rise to an attraction between them, as well as a short-range repulsion.The interaction is described in our simulation as with r the separation of the colliding particles.Here, C 4 and C 6 are the long-range attraction and short-range repulsion interaction coefficients, respectively.The choice of C 6 dictates the point at which the atom-ion interaction becomes repulsive.It was shown in [31] that its exact value is unimportant as long as the r −6 term starts to dominate at a short enough length scale of ∼1 nm.Meanwhile, the value of C 4 depends on the polarizability of the species of buffer gas atoms.To describe buffer gas cooling of an ion to the mK range, the short-range repulsion is enough [20].However for cooling beyond this, the induced-dipole interaction has been shown to play a significant role [30], as it is the mechanism by which the ion can be pulled from the trap centre into a region where micromotion heating limits its final temperature.Both terms are therefore crucial to our simulation of atom-ion collisions at the µK energy level.
During each simulated atom-ion collision, we use a fourth-order step-size-adaptive Runge-Kutta algorithm to propagate one atom towards the confined ion, under the influence of V trap (t, r) and V ai (r).We assume the atom cloud is dilute enough to ignore non-binary collisions between the ion and the atoms.The ion's initial kinetic energy is set to 0 K.The atom position is chosen at random on a sphere centred around the ion at the start of each collision.The sphere size is set to be large enough to accommodate large ion motion if and when the ion is heated.The atom's velocity in the radial (angular) direction is also selected at random from a Weibull (Gaussian) distribution corresponding to a chosen buffer gas temperature in the µK range.The two particles are then propagated until the atom leaves the sphere.After each simulated collision we record the ion's average kinetic energy (E kin ) and motional properties.These are not reset after each collision.We investigate the final values achieved after the ion has thermalized, usually after some hundreds of collisions.
To stay within a reasonable simulation time, the temperature of the buffer gas is set to 2 µK, which corresponds to the experimental values of ref. [29].
An example of the typical ion thermalization data obtained in our simulation is shown in Figure 1.For each simulation, we run at least 200 repeats of the thermalization process, averaging over them to find the mean ion energy reached after the nth collision (green line).
The blue error bars correspond to the standard deviation in these energies.In general, higher ion energies are accompanied by larger error bars.This is because, dependent on the mass ratio of the atom-ion species combination, the ion energy statistics after multiple collisions follow either a Tsallis or Maxwell-Boltzmann distribution [29,35,36].In both cases, the higher the mean ion energy, the larger the statistical spread of ion energies about this mean, as reflected in the error bars.To extract the final steady-state ion energy, we perform a weighted least-squares fit to our data.The weighting of the fit is determined by the uncertainty at each data point.We fit according to the equation Here, E kin,0 is the initial kinetic energy of the ion, Γ indicates the rate of thermalization and E kin,f is the kinetic energy reached by the ion at steady state.The latter two are fitting parameters.In figure 1, the fit is shown by the black curve, and the steady-state kinetic energy is indicated by the grey, dashed line.Note that n in this plot includes both Langevin and glancing collisions, meaning that the value of Γ extracted from the fit depends strongly on the initial atom-ion separation.However, within the range used in our simulation, this separation does not influence E kin,f as shown in ref. [31].
We focus our simulations on the Li/Yb + species combination, due to its proved effectiveness in reaching s-wave collision energies in buffer gas cooling experiments [29].The species choice has implications for the atom-ion interaction strength.Choosing 6 Li/ 171 Yb + corresponds to setting C 4 = 1.121 × 10 −56 Jm 4 .For the short-range repulsion we use C 6 = 2.804 × 10 −75 Jm 6 .The induced-dipole interaction has a characteristic range R * = µC 4 2 , with a corresponding energy scale given by E * = 2 2µ(R * ) 2 , which also represents the s-wave energy limit for the system.In a 6 Li/ 171 Yb + system, these take the values 69.75 nm and 8.58 µK, respectively.To reach deep into the quantum regime with buffer gas cooling, the kinetic energy of the ion must be reduced such that the collision energy of the atom-ion system E col = µ m i E kin,i + µ ma E kin,a is significantly lower than E * .Here, µ is the reduced mass.

III. THE LINEAR PAUL TRAP
In an ideal linear Paul trap, the time-dependent potential can be written as with α x = α y ≈ −α z /2 = −0.5 and α x = −α y = 1, α z = 0 [31,37].Here, v dc and v rf correspond to the direct-current (dc) and radio-frequency (rf) voltages applied to the trap electrodes.Meanwhile, R dc and R rf describe the fixed geometry of the trap and are approximately half the separation between opposing dc and rf electrodes, respectively.The frequency at which the rf voltage oscillates is given by Ω rf .The two terms of equation (3) are commonly referred to as V dc trap (r) and V rf trap (t, r), respectively.The resulting ion motion is described by the Mathieu equations, rj = −(a j + 2q j cos (Ω rf t)) Ω 2 rf 4 r j , with a j , q j the stability parameters of the ion trap for all three directions (j = x, y, z).The stability parameters can be defined in terms of the trap geometry and voltage signals as a x,y ≈ −a z /2 = 4ev dc /(m i R 2 dc Ω 2 rf ), and q ≡ q x = −q y = 2ev rf /(m i R 2 rf Ω 2 rf ), q z =0, with e the charge of the ion.The stability parameters combined with Ω rf determine the frequency ω sec,j ≈ Ω rf 2 a j + q 2 j 2 of the trapped ion's secular motion.
In an ideal trap, the ion's energy is mainly given by the kinetic energy, which is averaged over one period of the ion's secular motion.It is defined as and depends on the secular frequency and amplitude A j of the ion motion, and the trap parameters q and Ω rf .The latter two can be easily changed in experiments by changing the frequency and amplitude of the rf signal applied to the trap electrodes, and will be studied in our work.Typical numbers used in our Li/Yb + experiment are q = 0.5, Ω rf = 2π × 2 MHz and v dc up to 6.3 V, corresponding to ω sec,x and ω sec,y up to 2π×350 kHz, and ω sec,z up to 2π×120 kHz.The size of our trap is R dc = 5 mm and R rf = 1.5 mm, as outlined in ref. [38].
The ion's kinetic energy is connected to the number of excited motional quanta n j of the ion in each direction j by the following equations and T sec,j + T sec,x + T sec,y .
Here T sec,j , represents the secular temperature in direction j and k B is the Boltzmann constant.The kinetic energy contains both secular motion (first term) and intrinsic micromotion (second and third terms) contributions [39].Intrinsic micromotion (IMM) is driven by the rf field, and so depends on the radial temperature of the ion only.In our simulations, we focus on the average number of excited motional quanta in the radial direction (n = (n x + n y )/2), as the axial motion of the ion is largely unaffected by variations in our parameters of interest, q and Ω rf .
In addition to q and Ω rf , the ion's kinetic energy is also influenced by trap imperfections, which can lead to excess micromotion (EMM).Although trap-dependent, in a realistic simulation of the linear Paul trap these effects should be included [29,31].There are three forms of EMM which all scale differently with q and Ω rf .Axial EMM arises from the pickup of the rf signal on the dc electrodes, and contributes a kinetic energy (E ax ) that scales with q 2 and Ω −2 rf .Phase EMM, E ph , results from a phase difference between the rf voltage signals on opposing electrodes, and scales with q 2 and Ω 4 rf .Finally, E rad scales with q −2 and Ω −2 rf , and arises when stray electric fields in the radial direction of the trap displace the ion from the trap centre.The ion's total kinetic energy is therefore Here, b, c, d are fitting parameters indicating the contributions of the experiment-specific energies E rad , E ax and E ph .The fitting parameter g is added when an overall offset of E kin needs to be taken into account.In experiments, we can use ion position tracking or Ramsey microwave spectroscopy to measure stray fields to an accuracy of ∼0.01 Vm −1 [29].
It is feasible that radial stray fields can be compensated down to 0.05 Vm −1 .Note that the finite size of the trap can lead to q z = 0, which can give rise to additional EMM.However, since q z is much smaller than q and its EMM effects can be minimised by selection of an appropriate ion position along the trap axis, we assume it is small enough to be ignored in our simulations.
We investigate the values that E kin and n can reach after buffer gas cooling by simulating collisions between 6 Li atoms and an 171 Yb + ion confined in a linear Paul trap with varying levels of EMM.The temperature of the Li cloud is set to 2 µK.We vary q and Ω rf independently, to explore their impact, setting Ω rf = 2π × 2 MHz while q is varied and q = 0.219 while Ω rf is varied.The simulated trap shares its geometry with that of ref. [38] and has v dc = 0.787 V such that ω sec,z = 2π × 4.24 kHz.We include experimentally feasible levels of axial and phase EMM corresponding to an axial rf pickup of 2.02 Vm −1 (E ax = 42.4 µK) and a phase of 0.041 millirad (E ph = 9.18 µK) when q = 0.219 and Ω rf = 2π × 2 MHz [31,38,40].
We also vary the strength of the radial stray field around experimentally realistic values.
Our simulations show that the optimal values of q and Ω rf for buffer gas cooling depend on the radial stray fields present in our system.This is shown in figure 2. Blue, orange and green data corresponds to a trap with a radial stray field of strength 0 Vm −1 , 0.05 Vm −1 and 0.1 Vm −1 , respectively.Figures 2a and b show how the ion energy E kin varies with q and Ω rf .In general, the optimal trap parameters for minimising the ion energy take non-zero values, and depend on the strength of the radial stray field.During atom-ion collisions, the trap's rf field does work to heat the ion mediated by the long-range attractive atomion interaction [30].This heating competes with buffer gas cooling and leads to a trade-off between two effects.On the one hand, high values of q and Ω rf provide tight ion confinement, meaning that the pull from the atoms does not significantly alter its motion.On the other hand, low values of q and Ω rf will decrease the energy of the rf electric field, lowering the energy that it can donate to the ion during an atom-ion collision.The optimal trap parameters therefore provide a trade-off between these effects.Moreover, EMM provides additional trap-parameter-dependent heating.The higher the field-strength of stray radial electric fields that additionally pull the ion from the Paul trap centre, the higher the favoured optimal values of q and Ω rf .Appropriate selection of q is particularly important for optimal buffer gas cooling.This can be seen by the steeper slope of the fit lines in The lowest achievable kinetic energy depends significantly on the possibility to reduce electric stray fields and optimal choice of q and Ω rf .For E kin we find optimal values by fitting the data points using equation (7).The fit values can be found in Appendix A. For the chosen parameters, the lowest ion kinetic energy achievable by buffer gas cooling in our Paul trap is 64(1) µK, which can be achieved using trap parameters q =0.147 and Ω rf = 2.28 MHz.In a 2 µK atom cloud, this corresponds to E col = 0.59×E * , significantly below the s-wave limit.It is also two times lower than the lowest collision energy achieved in previous experiments [29], which is indicated by the red data point in a.This experiment yielded collision energies 1.15 × E * by buffer gas cooling an ion to kinetic energies of 193 µK.The lowest achievable ion kinetic energies indicated by our simulation are therefore roughly 3 times colder than what has been achieved in experiments.Note that the discrepancy between the experimental data point and the simulation predictions can be explained by an additional cooling step of adiabatic decompression that was used in the plotted experiment, but is not included in our simulation.Our results suggest that with careful parameter optimization, this final step is not necessary for reaching low temperatures.The value 64 µK is a lower limit on what can be achieved for the ion in this setup with an atomic bath of 2 µK.Experiments will be limited to higher energies due to background heating from technical noise [38], which varies across experiments and is not included in our simulations.Quantum effects may also play a role in reducing the rate of buffer gas cooling at very low temperatures, as the rate of atom-ion collisions is expected to be lower once E col reaches below the s-wave energy [29].
More details on a quantum treatment of the atom-ion interaction can be found in section VI.
Depending on the stray fields, n is minimized for different values of the trap parameters to E kin .This can be seen when comparing figure 2a and b to c and d.Each of these variables is affected differently by sources of EMM.Therefore, to optimize buffer gas cooling experiments, one should consider whether to prioritize reaching quantized ion motion or swave collision energies to inform the choice of trap parameters.The lowest average number of excited motional quanta to which the ion can be cooled by an atom cloud at 2 µK in our setup is 0.92 (7), at a stability parameter q = 0.26.This is 4 times smaller than what has been achieved in buffer gas cooling experiments, as indicated by the red data point in c.It allows us to reach quantized ion motion, making buffer gas cooling competitive with some sub-Doppler laser cooling techniques, such as Sisyphus cooling [34].However, it is not capable of reaching n 1, as is possible using, for example, resolved sideband cooling [24].

IV. OTHER TYPES OF RF ION TRAPS
Since the ion's final kinetic energy is limited by the trap parameters and excess micromotion in a linear Paul trap, this raises the question of whether other trap configurations are better suited for buffer gas cooling.Therefore, in this section, we simulate a digital ion trap (DIT), multipole ion trap (MIT) and rotating ion trap (RIT).We investigate how the dependence of E kin on the trap parameters may differ for different traps and, where relevant, consider the roles played by new parameters that are unique to certain traps.
The DIT is a modification of the previously-discussed linear Paul trap, whereby the rf voltage is driven digitally (stepwise constant) instead of analogue (sinusoidally).The switching of the voltage leads to times where the ion dynamics are that of a free particle as the rf voltage is zero and this might lead to more efficient buffer gas cooling in these periods.
However, collisions occurring during a voltage switch will experience a violent change in potential, which may result in larger ion heating than in the traditional analogue Paul trap.
Similar to equation ( 3), the time-dependent potential of a DIT can be described by Here, the variables are equivalent to those in equation ( 3) apart from the shape of the rf signal P DIT (t), which is defined as Here, T = 2π/Ω rf is the period of the signal, t r is the remainder t = κ × T + t r [41], with κ an integer, and τ is the fractional pulse width.The digital voltage driving is shown in Figure 3a.For a DIT, τ adds an additional trap parameter that affects the stability of the trap [42].Approximating the digital signal with the first term of the corresponding Fourier series [43], the q parameter for a DIT becomes q x,y,DIT = 16ev rf sin (τ π)/(πm i R 2 rf Ω 2 rf ).We investigate how E kin varies when buffer gas cooling an ion in a DIT, dependent on the temperature of the atomic bath T a , and the trap parameters q and τ .The results can be seen in figure 3b and c.In these simulations, we use the same geometry, trap parameters and atomic bath temperature as in the linear Paul trap simulations, except where otherwise stated.In figure 3b, we choose τ = 0.285 and v rf = 34.467V such that we have q = 0.219, to match the previous simulations.In figure 3c, we vary q and τ while keeping the trap geometry constant.For a fixed τ , the variation in q corresponds only to a change in the voltage of the digital pulses v rf .The black lines show the trends as given by an analogue Paul trap with the same q, Ω rf and EMM parameters.Radial stray field strengths are set to 0 or 0.1 Vm −1 in these simulations, while sources of axial and phase MM are not included.
Our simulation shows that buffer gas cooling using a digital rf driving results in equal ion temperatures to those in an analogue linear Paul trap.In figure 3b, DIT results are in excellent agreement with the analogue Paul trap trend line.Therefore, for the chosen τ and q of this plot, the energies reached in the DIT are identical to those of an analogue trap, irrespective of the atomic bath temperature.In figure 3c, we see that in a DIT of variable q, the data points corresponding to the lowest E kin also lie on the trend line of the analogue trap.The range of q for which the DIT and analogue trends match is affected by the fractional pulse width τ .While τ provides an extra degree of freedom for tuning the trap stability, its tunability does not lead to lower ion kinetic energies than are achievable with an analogue Paul trap.The fact that the DIT results are so similar to those of the analogue Paul trap suggests that the time scale of an atom-ion collision is longer than the period of the rf trap, such that the shape of the rf voltage signal during the collision period is unimportant [30].

B. Rotating ion trap
In a RIT, ion confinement is provided by a rotating potential, instead of the oscillating potential typical for linear Paul traps.The ion motion is described by four compoments, instead of the infinite sum of Fourier components needed for a linear Paul trap.Since fast MM plays an important role in atom-ion collisions, we may therefore expect a different outcome for buffer gas cooling when comparing the two traps.
A RIT time-dependent trapping potential can be written as [44,45] Here, variables are equivalent to those in equation ( 3), yet the rf driving part has changed.
The resulting ion motion is then rx = −(a x r x +2q x (r x cos Ω rf t−r y sin Ω rf t)) 4 .An important consequence of this is that the ion's motion is now circular, instead of oscillating, and consists of four Fourier components: two secular motions and two IMM components.The frequency of the secular motion in the RIT can be written to first order in q as Here a = a x = a y , and q = q x = −q y = 2ev rf /(m i R 2 rf Ω 2 rf ).Note that if a is very small, the secular frequency differs from that of a linear ion trap by a factor √ 2. Therefore, the ion confinement in a RIT is twice as tight as in a linear trap with the same trap parameters.
We simulate buffer gas cooling of an ion in a RIT of variable q and Ω rf , and explore how the dependence of E kin on these parameters compares to the case of the linear Paul trap.
Results can be seen in figure 4. We set the same ion trap and atom cloud parameters as in our linear trap simulations.We simulate radial stray field strengths of 0 and 0.1 Vm −1 and do not include sources of axial and phase EMM.We set the stability parameter q set,RIT in terms of the trap geometry and driving signal as q set,RIT = q lin = 2ev rf /(m i R 2 rf Ω 2 rf ), which is the same definition as used in the linear Paul trap.However, due to the rotation of the potential, its effect on the ion energy differs from that of q lin .We simulate buffer gas cooling in a RIT with and without a stray electric field of 0.1 Vm −1 , as shown by the blue circles and triangles in figure 4. The data is compared to the equivalent linear trap trends, represented by the black dotted and dashed lines.
We find good agreement between the RIT data and the linear Paul trap trend at low q when using q eff,RIT = √ 2q lin to rescale the latter.Here, q eff,RIT refers to the effective stability parameter of the RIT and q lin is that of the linear trap.In figure 4a the x-axis of the black trend line is rescaled by √ 2 to obtain the match with the RIT data at low q, to within 3 µK.The black line in b is also in good agreement with the RIT data when we use the same scaling of q lin .The E kin dependence on Ω rf and q after this adjustment is thus similar for both traps, and significantly lower energies cannot be reached in the RIT.In figure 4a, the minimum E kin is achieved in the low-q region.Only at high q do the linear and RIT data differ significantly, and in this limit E kin is higher in the RIT.Note that the energies reached in this graph cannot be compared to those in figure 2, as the latter includes sources of axial and phase MM, making the system inherently hotter.We conclude that although the q set,RIT we set in the simulation is defined identically in v rf , Ω rf and R rf to q lin , the effect of the rotating potential is that of a linear trap that is a factor √ 2 higher in q.In a RIT, we can therefore use a lower-amplitude voltage signal v rf to achieve similar ion temperatures as in a linear Paul trap.

C. Multipole ion trap
The benefits of using a MIT for buffer gas cooling of molecular and atomic ions on the ∼ K-mK level has been extensively explored [46][47][48].It has been shown that an increased number of trap poles minimises rf heating effects, allowing more efficient buffer gas cooling than would otherwise be achievable, as well as stable cooling of systems with a smaller atomion mass ratio than is possible in quadrupole traps [49].Here we investigate whether the same benefits extend to cooling of the high mass-ratio Li/Yb + system, beyond the lowest possible energies that can be reached in a quadrupole trap.
In a MIT, the four electrodes of the Paul trap are replaced by a higher number of electrodes, which are used to create an rf potential with N poles.This requires N ground and N rf electrodes.For a Paul trap N = 2, while a hexapole trap uses N = 3 and an octupole trap has N = 4.The resulting potential can be written as [50] V trap (t, r) = V dc trap (r) + Here we use cylindrical coordinates r ⊥ = r 2 x + r 2 y and φ = tan −1 (r y /r x ).The resulting motion in each direction j = x, y, z will depend on N and can be expressed as rx = −(a x + N (N −1) η x cos Ω rf t) Here, all parameters are defined as in equation (3), except now the rfdependent motion is written in terms of a generalized stability parameter η j .We define this in terms of q as such that η x = −η y = q when N = 2. Here, q = q x = −q y = 2ev rf /(m i R 2 rf Ω 2 rf ) and q z = 0.
For N > 2, η j has a r ⊥ -dependence that is not present in a quadrupole trap.Depending on its energy, the ion oscillates around the trap axis with an amplitude A ⊥ .Therefore, in a MIT used for buffer gas cooling, the stability parameter η j changes throughout the cooling process as the ion's energy decreases.In experiments, to ensure that η(A ⊥ ) is kept roughly constant, changes to A ⊥ can be compensated by ramping the other trap parameters, for instance v rf can be easily ramped.
We simulate buffer gas cooling of an ion in a hexapole (N = 3) and octupole (N = 4) trap of variable stability.We use the same trap parameters as in the Paul trap from section III, except for v rf , which we use to tune the trap stability η j .We set v rf to define η j (r ⊥ ) = η at a chosen ion position r ⊥ = 100 nm.We repeat the simulations with and without a radial stray electric field of 0.05 Vm −1 .Attempts were made to simulate the process at stronger fields, however for N > 2 this frequently resulted in a very unstable trap, as the ion was pulled to high r ⊥ .Propagating the ion under these conditions is beyond the capability of our simulation.Similarly, for simulations with N = 4 and η > 0.2, frequent stability problems were encountered, making it difficult to repeat the simulations.Note that these simulations do not include sources of axial and phase EMM.
The greater the number of poles of the MIT, the more efficient the buffer gas cooling process.Figure 5 shows the variation of E kin with η for each trap, with (b) and without (a) stray electric fields.We vary η between 0.1 and 0.5.Black lines show results for a quadrupole trap (N = 2), while red and olive data points show results for the hexapole and octupole traps, respectively.In a, the trends of the hexapole and octupole data follow similar shapes to the lines of the quadrupole trap, but with a shallower gradient and an overall lower E kin .
Meanwhile in b we see indications that the optimal value of η in an octupole trap is higher than in a quadrupole trap, and gives access to lower energies.The gradients of the MIT trends in b are also shallower than in the quadrupole trap.The lowest temperatures are therefore reached in traps of higher N , and the consequences of a sub-optimal choice of η are most severe for traps with fewer poles.
To achieve the lowest possible ion energies in a MIT of high N requires well-compensated stray electric fields and a high-amplitude voltage signal.This makes the use of high-order traps for buffer gas cooling practically challenging.When cooled to the µK range, an ion in the simulated trap will oscillate with an amplitude on the scale of 10 2) , traps of increasing N require increasing v rf to maintain a reasonable The dotted and dash-dotted lines are the fit to the data from the corresponding quadrupole trap without and with the radial stray field, respectively.Errors are standard errors and are sometimes smaller than the point size.
stability.For example, an octupole trap with η = 0.2 requires v rf = 4.72 GV, while the same stability parameter can be achieved in a quadrupole trap with v rf = 62.95 V.This also has severe implications for axial EMM in the trap, as these voltages are likely to lead to significant rf pickup on the dc electrodes.Even with a suitable choice of v rf , the N -dependent scaling of η with r ⊥ means that traps with high N are more likely to become unstable when stray electric fields are present.These practical limitations should be considered when selecting how many poles would make a MIT most suitable for buffer gas cooling.

V. HYBRID OPTICAL-PAUL TRAP
Another possible extension of the Paul trap is to combine it with an optical tweezer for the ion.This hybrid optical-Paul trap would give an additional confinement for the ion which could benefit the attainable ion kinetic energies and suppress the impact of stray fields.Several approaches of using optical ion trapping in atom-ion setups have been experimentally explored [51][52][53], as an optical trap avoids the ion-heating effects present in time-dependent trapping potentials.However, the optical trap for an ion is typically less deep than an rf trap, meaning ions that gain energy due to inelastic collisions are more easily lost.Here, we therefore focus on a hybrid approach, designed to combine the benefits of rf and optical traps.Unlike the traps discussed so far, a deep optical-rf hybrid trap will also affect the density and temperature of the atom cloud.This leads to a tightly-compressed buffer gas and off-resonant heating of the atoms.However, here we focus on the effect on the ion alone, and assume in our simulations that the atoms are unaffected by the optical potential.
The modeled hybrid optical-Paul trap consists of a quadrupole linear Paul trap with an optical tweezer directed along its axis for additional ion confinement.The setup is shown in the inset of Figure 6.This corresponds to using the rf potential of ( 3) and adding to it a Gaussian potential Here r ⊥ = r 2 x + r 2 y , w 0 is the beam waist of our optical tweezer, z R = πw 2 0 λ is the Rayleigh length, and λ is the frequency of the tweezer beam.The effective trap depth U is determined by the beam power P and the polarizability α of the species in the trap.It is defined as , with 0 the permittivity of a vacuum and c the speed of light.The effect of the optical potential on the ion motion is to give an additional position-dependent force towards the trap centre Here ω opt,j is the trapping frequency of the optical trap in each direction ω opt,x,y = . The confinement of the tweezer in the axial direction is sufficiently weak for us to assume that ω opt,z plays no significant role here.
We simulate buffer gas cooling of an ion in this hybrid optical-Paul trap, as a function of T a and q for varying optical confinements.The contribution from the linear Paul trap is taken with identical trap parameters and geometry to those used in section III.We choose a stray radial electric field of either 0 or 0.1 Vm −1 .For the former, no other EMM effects are included.The optical tweezer has a wavelength of 1070 nm, a variable power set to 0.423, 4.23 and 2.12 W and a variable beam waist set to 2, 2 and 1 µm, respectively.The Stark shift of a 171 Yb + ion in such a tweezer gives it a dynamic polarizability α of approximately 66 au (1.09×10 −39 Jm 2 V −2 ) [55].The chosen tweezer parameters therefore correspond to respective trapping frequencies ω opt,x,y /ω sec,x,y = 0.230, 0.728 and 2.05 for the 171 Yb + ion.
Here, trapping frequencies are given as a ratio of the linear Paul trap's secular frequency, as defined in section III.FIG. 6. Ion energy after buffer gas cooling in an optical-rf hybrid trap.The tweezer has a radial frequency ω opt /ω sec of 0.23 (brown), 0.728 (orange) and 2.05 (pink).a Variation with trap parameter q.The black trend lines are simulations for a Paul trap alone with varying radial stray fields of 0.1 (black, dotted), 0.05 (black, dash-dotted) and 0 Vm −1 (black, dashed).In this plot, the stray electric field for the data points is 0.1 Vm −1 .Inset: illustration of a hybrid optical-rf trap.b Varying buffer gas temperature.The black, dashed line represents a Paul trap without radial stray fields, while the grey dash-dotted line represents a unity gradient.Data points are taken for a setup with no radial stray field.Errors are standard errors and are sometimes smaller than the point size.
Our simulations show that adding an optical tweezer aids in decreasing the effects of radial stray fields on the ion cooling.Therefore at low q, the ion can be cooled to temperatures significantly lower than those accessible in a Paul trap alone.Figure 6a shows the simulation data of E kin for variable q in a hybrid trap with stray electric fields of 0.1Vm −1 .The data is compared to the trend lines from a pure Paul trap with varying levels of stray fields.All lines and data points in this plot also include E ax = 42.4 µK and E ph = 9.18 µK at q = 0.219 and Ω rf = 2π × 2 MHz.For a hybrid trap with increasing ω opt , the data points tend increasingly towards the trend line of an rf trap with no stray electric fields.In other words, the tighter the optical tweezer, the smaller the heating effects of radial EMM, allowing the ion to reach much lower temperatures than in an equivalent regular Paul trap.The same suppression of EMM effects is seen when analysing simulations of a hybrid trap with variable Ω rf .
A tighter confinement of the optical tweezers leads to a more efficient buffer gas cooling process.Figure 6b shows the linear relationship between the ion energies reached in a hybrid trap of variable optical depth, and the temperature T a of the buffer gas in which it is submerged.We also plot the same relationship for a Paul trap without an optical tweezer as well as a line with unity gradient, for comparison.The latter corresponds to a static trap.We find that the hybrid trap results follow lines of increasingly shallow gradient for an increasingly tight optical tweezer.This reflects the time-independent tweezer potential becoming increasingly dominant over rf contributions in determining the energy reached by the ion.The value of E kin,i at T a = 0 K also decreases for an increasingly tight tweezer.
To benefit from the optical tweezer, it must be perfectly aligned with the Paul trap axis.
Figure 7 shows the simulation of buffer gas cooling in a hybrid trap in which the optical tweezer runs parallel to the Paul trap axis, with a varying radial offset ∆ between them.
For comparison, we also indicate the ion energy reached in an equivalent Paul trap with no optical tweezer.Offsets of 50 nm or greater lead to ion energies above this value, and the greater the offset, the higher the ion energy.For misalignment, the effect of the tweezer is to pull the ion from the Paul trap axis, providing an additional source of excess micromotion, akin to a radial stray field.
Another important limitation for practical implementation of the optical-rf hybrid trap is the effect of the optical potential on the atoms.At the high trapping frequencies required to observe a significant benefit of the tweezer on the ion, heating and confinement of the atoms become significant.If we wish to trap an Yb + ion in a tweezer with λ = 1070 nm and ω opt,x,y = 2.05 × ω sec,x,y , a cloud of Li atoms in the path of the beam will experience a trap depth of 82 mK [54].This will lead to a tightly-compressed atom cloud.Combined with the high photon-intensity of the tweezer, these atom densities give a possibility for unwanted chemistry [16,56].The same tweezer would also lead to heating of the Li atoms due to photon scattering at a rate of 128 µKs −1 , significantly reducing the efficiency of buffer gas cooling.To address the problem of trap depth, it may be possible that for some atom/ion species combinations one could confine the ion in an optical tweezer with a wavelength equal to a tune-out wavelength of the atoms, as has been proposed for multi-species atom experiments [57].However for our Li/Yb + species, this is not possible, as the only Li tuneout wavelength that is red-detuned for Yb + is too close to the D1 and D2 lines of Li [58].
Suitable combinations of bichromatic optical dipole traps can also be used for all-optical trapping of both atoms and ions [33].We consider this as an alternative method of reducing the effect of the tweezer on the atom cloud in an optical-rf hybrid trap.However, the cooling in a bichromatic trap is highly sensitive to relative position fluctuations between the two trapping beams.Furthermore, for our species and desired trap depth, the wavelength and intensity requirements cannot be met simultaneously with the requirement for a large detuning of each beam from the atom and ion transition lines, to avoid off-resonant heating.
In fact, the heating problem would be intensified by the higher-frequency beam.
A way to avoid the negative effects of the tweezer on the atoms may be to use a two-ion crystal, with one ion submerged in the atomic bath and the other in the tweezer.Using the ion chain's collective modes of oscillation, the pinned ion could then be used to sympathetically enhance the cooling of the ion in the buffer gas.We simulate this option and show the results in figure 8.It shows the energy reached by the untweezed ion in the crystal, dependent on the trapping frequency of the optical tweezer.We also show for comparison E kin for a single ion (dashed line) and a two-ion crystal (solid line) after buffer gas cooling in a Paul trap without an optical tweezer.The results show that the higher the trapping frequency of the optical tweezer, the lower the achievable ion energies.However, these energies tend towards a minimum equal to the energy of a single ion in a Paul trap, for large tweezer depths.The effect of the tweezer on the two-ion crystal is to decouple the motion of the two ions, thereby giving no additional cooling benefits.

VI. THE QUANTUM DYNAMICS OF MICROMOTION-INDUCED HEATING
Up until now, we have used classical physics to describe the dynamics of the trapped ion in the buffer gas.However, at the simulated temperatures, quantum effects may play a role.
In particular, the angular momentum in the atom-ion collision and the energy of the ion in its trap will both be quantized.Furthermore, the dispersion of the atomic wavepacket may become much larger than the micromotion amplitude of the ion such that it cannot be resolved at all in a collision.
Specifically, when looking into micromotion-induced heating due to the time-dependence of the trap, the classical and quantum mechanical picture might differ.In the classical simulations of colliding atoms and trapped ions, the heating stems from a temporary binding between the ion and atom during which multiple collisions can occur [30].At some point, and depending on the phase of the rf trap, the ion and atom can violently separate again.
During the whole process, energy can be absorbed from the rf field, leading to heating of the system.One may wonder whether this behaviour survives in a quantum treatment, in which binding can only occur in discrete molecular states with energy level spacings that may be larger than the collision energy.In this section we answer the question: Can quantum mechanics save us from micromotion-induced heating in an atom-ion collision?
One approach is to consider the same case that we studied in the classical regime, but this time launching an atomic wavepacket at an ionic wavepacket that is held in a (quantized) rf Paul trap.Floquet theory allows us to solve the Schrödinger equation, similar to previous theoretical investigations into mechanically-controlled interactions between trapped atoms and ions [59][60][61][62][63] and individually-trapped atoms [64].However, even if we restrict ourselves to one dimension to reduce the complexity, the required Floquet Hilbert space is very large and convergence is very slow.Instead, we follow a qualitative analysis of the problem, showing that quantum mechanics cannot eliminate micromotion-induced heating.
A. Quantum description of an ion in a Paul trap We start by studying the ion motion quantum mechanically, when the ion is confined by an rf-potential Vrf (t) = where ω sec,i takes the same value as in our classical equations (see section III).Here, the ion intrinsic micromotion Hamiltonian Ĥmm (t) is given by Ĥmm Furthermore, g = [2(1+2a/q 2 )] −1/2 and {., .}denotes the anti-commutator.Note that excess micromotion is not included in our quantum study.
Typically in analysis of the trapped ion, one would proceed with applying the secular approximation.The basis for the secular approximation is to neglect Ĥmm (t) such that the dynamics of the ion is given by a static harmonic oscillator alone.Note that the relative size of the terms in Ĥmm (t) are of a similar size to those of the static harmonic oscillator and in fact, the mathematical basis for the secular approximation is quite delicate.The reason it works at all is the structure of the eigenstates of the harmonic potential.The time-dependent part Ĥmm (t) contains terms ∝ x2 i .Taking the Fock states of the static harmonic oscillator as a convenient basis in which to study the Hamiltonian Ĥi (t), we see that strict selection rules apply, i.e. n|x 2 i |m = 0 except for n = m ± 2 and n = m, with |n denoting the n-th Fock state.Since Ω rf ω sec,i such transitions cannot be resonant.
Therefore, the only effect of Ĥmm (t) here is to cause a small increase in the effective secular trapfrequency that becomes more prominent for large q.
However, any anharmonicity in the trapping potential can cause a serious breakdown of the secular approximation.This is because transitions |n → |m close to the resonance condition |ω sec,i (n − m)| = Ω rf are allowed in anharmonic traps.The observation of such parametric resonances in rf traps, in which ions show rapid heating, dates back a long time [66][67][68].Thus, for anharmonic potentials, the contribution of Ĥmm (t) can not be ignored.A similar effect is expected when describing an ion in the presence of a buffer gas, which adds an additional atom-ion interaction potential.

B. Quantum description of atom-ion interactions
Now we add the atom to our quantum description.Then, the Hamiltonian is given by where Vai = −C 4 /r 4 is the atom-ion interaction potential, r is the separation between the atom and the ion and C 4 takes the same value as that in our classical simulations (see section II).We transform to the relative (r) and centre of mass (R) coordinates, to obtain the atom-ion dynamics [62].In these coordinates, eq. ( 18) can be written as with ω R = ω sec,i m i /M and ω r = ω sec,i m a /M denoting the centre of mass and relative frequencies of the system, respectively.The micromotion Hamiltonians are Ĥrr mm (t) = −µg 2 ω 2 r r2 cos(2Ωt) − gω r {r, p} sin(Ω rf t) , Here P and p are the centre of mass and relative momenta, respectively.
To understand where heating can arise in this system, we must look at each of the above equations individually.Let us begin with the centre of mass Hamiltonian eq. ( 19).Excluding micromotion effects (i.e.ĤRR mm (t) = 0), the centre of mass motion can be treated as that of a harmonic oscillator, with solutions that are Fock states |n R and energies n R ω R .If we include micromotion, ĤR (t) resembles the Hamiltonian Ĥi (t) of a single ion in a Paul trap, as described in section VI A, now with m i → M and ω sec,i → ω R .Since this Hamiltonian is without anharmonicity, it describes a stably trapped particle for which we may employ the secular approximation.As a result, looking at the centre of mass motion alone, it would seem that no significant heating is expected.
The relative Hamiltonian, eq. ( 20), takes a similar form to that of the centre of mass, except it includes the additional static potential Vai .Provided that there is no micromotion (i.e.Ĥrr mm (t) = 0), this static term does not give rise to heating effects.However, it does mean that, unlike the centre of mass Hamiltonian, the solutions to the unperturbed relative Hamiltonian are not Fock states.Instead the dynamics must be obtained by solving the Schrödinger equation by means of quantum defect theory [62].In the limit where r → 0, the energy is dominated by the atom-ion interaction, and we obtain even and odd solutions given by [59,60] ψe (r) ∝ |r| sin ψo (r) ∝ r sin where we set the even and odd short range phases as φ e = π/4 and φ o = −π/4.We can use q = 0.5, similar to the experimental parameters of ref. [29].The results are shown in Fig. 9.
We see that the couplings due to micromotion can be up to ∼ 0.1 ω sec,i for an energy difference δE ∼ Ω rf .This energy scale is comparable to the fractional change of kinetic energy K for a classical hardcore collision.Assuming the atom is at rest initially and the initial ion kinetic energy is K 1 , we get for the final ion kinetic energy K 2 = 1 − 2β (1+β) 2 K 1 = 0.93K 1 , with β = m a /m i .For a cold ion K 1 ∼ ω sec,i , giving K 1 − K 2 ∼ 0.07 ω sec,i .We therefore conclude that buffer gas cooling to the ion ground state is not possible, even in the quantum picture, as micromotion-induced heating will impose a limit on the ion temperature, similar in scale to that seen in classical analysis.

VII. CONCLUSION AND OUTLOOK
In this work, we demonstrate that careful optimization of ion trap parameters is crucial for reaching cold temperatures when buffer gas cooling ions with ultracold atoms in timedependent traps.We show that, under appropriate conditions, buffer gas cooling should allow one to reach the regime of quantized ion motion.For the Li/Yb + system, parameter optimization can allow buffer gas cooling to ion kinetic energies three times colder than the previous best experimental values, corresponding to collision energies approximately 0.6 times the s-wave limit.This makes buffer gas cooling capable of outperforming laser Doppler cooling in some situations, and it is even competitive when compared to some sub-Doppler cooling techniques such as Sishyphus cooling.These low energies are achievable by sup-pression of micromotion-induced heating effects, through optimization of q and Ω rf , and by compensation of trap imperfections to low, but experimentally realistic, levels.Additionally, we benefit from a large ion-atom mass ratio.For equal mass systems, similar trap parameter dependencies were found for a linear Paul trap [69], however the achievable temperatures for such systems remain far from the s-wave regime.
We further demonstrate that the lowest possible energies achievable in a linear Paul trap can also be reached through similar parameter optimization in a range of rf traps.
In the case of the multipole and hybrid optical-rf traps, even lower temperatures can be reached, albeit with the requirement for challenging experimental parameters and, for the hybrid trap, a method of avoiding photon-scattering heating effects.In addition, we provide a quantum description of the atom-ion interaction, and argue that micromotion-induced heating continues to play a significant role in this regime.
Our study shows a route towards reaching deeper into the quantum regime with timedependent traps, and presents buffer gas cooling as a feasible alternative to laser cooling techniques.Especially for studies of quantum chemistry and precision spectroscopy, cold controlled collisions in hybrid atom-ion systems are important.Moreover, in the quantum regime one can benefit from the tunability of the atom-ion interaction through Feshbach resonances [25,26,28] and use the ground state cooled-ion for high-precision studies of charged impurity physics [8,12].

VIII. ACKNOWLEDGEMENTS
Fit parameter Trend in q Trend in Ω rf A x 2.483×10 −8 1.335×10   2, but also to data from similar simulations that excluded one or more forms of EMM.
Fig 2a, compared to b.

FIG. 2 .
FIG.2.Calculated ion energy and motion after buffer gas cooling in a linear Paul trap.The ion's steady-state kinetic energy E kin and number of excited motional quanta n are shown for cooling inside a trap of variable stability parameter q and radio frequency Ω rf .Blue, orange and green data corresponds to a trap with a radial stray field of strength 0 Vm −1 , 0.05 Vm −1 and 0.1 Vm −1 , respectively.Red data points in a and b show experimental results from ref.[29], while the red dashed lines in c and d indicate the value of Ω rf used in the same experiment.The lines in a and b correspond to equation(7), fit to the data using fitting parameters b, c, d, g and A x,y,z .The fit values of these parameters can be found in Appendix A. Errors are standard errors and are sometimes smaller than the point size.

FIG. 3 .
FIG.3.Ion energy after buffer gas cooling in a DIT. a Analogue (black, dashed) versus digital (purple and red) driving of the rf voltage.The fractional pulse width τ defines the driving.Digital signals are shown for τ = 0.285 (purple) and 0.5 (red).b The kinetic energy reached for variable buffer gas temperatures, for q = 0.219 and τ = 0.285.c The kinetic energy reached for variable trap parameter q for τ = 0.285 (purple), 0.4 (green) and 0.5 (red).The black dotted lines represent the trend in an analogue linear Paul trap.Errors are standard errors and are sometimes smaller than the point size.

FIG. 4 .
FIG.4.Ion energy after buffer gas cooling in a RIT.Results are shown for traps of various q (a) and Ω rf (b), with (circles) and without (triangles) a radial stray field of 0.1V/m.The black dotted and dashed lines represent the trend in a linear Paul trap with similar trap parameters.Note that the black lines have been adjusted according to the relation q eff,RIT = √ 2q lin (see text).Errors are standard errors and are sometimes smaller than the point size.

FIG. 5 .
FIG.5.Ion energy after buffer gas cooling in a multipole ion trap.Results are shown for an ion in a hexapole (red) and octupole (olive) trap, without (a) and with (b) a radial stray field of 0.05 V/m.The dotted and dash-dotted lines are the fit to the data from the corresponding quadrupole trap without and with the radial stray field, respectively.Errors are standard errors and are sometimes smaller than the point size.

FIG. 7 .
FIG.7.Ion energy after buffer gas cooling in a displaced rf-optical hybrid trap.Orange data points mark the temperatures reached for a variable radial offset (∆) between the Paul trap axis and the optical beam.The black dashed line shows the energy achieved when simulating the same process with the dipole beam switched off.Errors are standard errors and are sometimes smaller than the point size.

1 |EFIG. 9 .
FIG. 9. Matrixelements D nr,n r = | n R | r|H mm,1 |r |n R | versus the energy difference δE = |E nr − E n r |.Data is shown for a system of6 Li and 171 Yb + , confined in a Paul trap with parameters Ω rf = 2π × 2 MHz, a = 0 and q = 0.5. ) FIG.8.Energy of one ion in a two-ion crystal after buffer gas cooling, with the second ion confined by an optical tweezer outside of the atom bath.Data points show results for a tweezer of varying depth.For comparison, the black lines show the temperature reached by an ion in a Paul trap without a tweezer for the single ion case (black, dashed) and two-ion case (black, solid).Errors are standard errors.The grey shaded regions indicate standard errors on the lines.
[65] i Ω rf (a + 2q cos (Ω rf t)) x2i .Its motion can be described by the Hamiltonian Ĥi (t) = p2i /2m i + Vrf (t), with wavefunction |Ψ i (t) and momentum pi .Using the Cook and Shankland transformation the wavefunction of the ion can be written as|Ψ CS (t) = exp [−i/(4 )m i qΩ rf x2 sin (Ω rf t)] |Ψ i (t)[65]. Te fast dynamics is now contained in the exponential factor and the slow-varying part of the ion wavefunction evolves under the Hamiltonian Ĥi (t).The rf-potential Vrf (t) can be separated into a secular and intrinsic

TABLE I .
Fit parameters for the ion energy and motion after buffer gas cooling in a linear Paul trap.These parameters together with equation (7) are used to obtain the lines in Fig2a and bfor varying radial micromotion contributions.Note that these results were obtained not only by fitting to the data shown in figure