Dark resonance spectra of trapped ions under the influence of micromotion

We study the influence of micromotion on the spectrum of trapped ions with a lambda-type level scheme, leading to dark resonances due to coherent population trapping. We work with calcium ions trapped in a ring-shaped Paul trap, in which one can compensate excess micromotion for only one ion of the crystal. We observe that micromotion affects the shapes of the dark resonances and causes the appearance of “echoes” separated by intervals given by the drive frequency. We present a theoretical model that provides good fits to the measurements and can be used to estimate the amplitude of the micromotion modulation of the atomic motion. We estimate an effective temperature of the ions from the spectra and observe clear micromotion heating as well as impaired cooling for sufficiently large excess micromotion.

1 Introduction

Trapped-ion systems are versatile platforms which in the last decades led to the development of quantum technologies such as quantum simulators (Porras and Cirac, 2004; Friedenauer et al., 2008), optical clocks (Rosenband et al., 2008; Chwalla et al., 2009; Chou et al., 2010; Huntemann et al., 2012), and small-scale prototype quantum computers (Cirac and Zoller, 1995; Blatt and Wineland, 2008). They have also served as a platform for the study of fundamental physics at the quantum level and for the search of new physics via precision measurements (Berengut et al., 2018; Safronova et al., 2018; Shaniv et al., 2018; Keller et al., 2019; Counts et al., 2020).

Paul traps, and more generally radiofrequency (RF) traps, employ electric RF fields to confine ions in vacuum (Neuhauser et al., 1980; Paul, 1990). The dynamics of the ions within these traps can generally be understood in terms of two contributions: macro or secular motion, and micromotion. For appropriate trap parameters, the oscillating RF field creates an effective static harmonic potential with a characteristic frequency that is much slower than that of the RF field; the movement of the ions due to this effective potential is called secular motion. This motion can be externally controlled by cooling, coherent excitation and preparation of quantum motional states (Leibfried et al., 2003). On the other hand, the oscillating character of the trapping field creates a driven motion at the RF drive frequency, named micromotion. The micromotion amplitude is generally small and, in lowest-order approximation, proportional to the displacement of the ions with respect to the center of the trap, where the RF field vanishes.

As has been shown in the last years, it is possible to use micromotion as a resource to enhance the speed and fidelity of entangling gates in trapped ions (Bermudez et al., 2017; Ratcliffe et al., 2020; Gaebler et al., 2021). Nevertheless, it remains generally desirable to have low micromotion, for example, to achieve efficient ion cooling (DeVoe et al., 1989; Cirac et al., 1994; Sikorsky et al., 2017), reduce heating (Turchette et al., 2000; Brouard and Plata, 2001), and reduce Doppler related shifts in trapped ion clocks (Keller et al., 2015; Ludlow et al., 2015). With this goals in mind, several techniques were developed that allow for the reduction of micromotion, usually referred to as micromotion compensation (Berkeland et al., 1998; Ibaraki et al., 2011; Higgins et al., 2021).

Micromotion-related spurious effects are particularly problematic in ion crystals. In this case, the micromotion will affect the ions differently depending on their locations within the trap (van Mourik et al., 2022). The geometry of the trap and the potentials used then determine the possible ion crystal shapes and the properties of the micromotion for each ion. For instance, the linear Paul trap is specially well-suited to work with linear ion chains since ions can be made to lie along an axis of vanishing oscillating field. On the other hand, in three-dimensional (3D) traps like ring traps the RF field vanishes at only one point in space, so excess micromotion can at best be fully compensated for a single ion.

Micromotion has important modulation effects in the atomic spectra of ions, which are different depending on the width of the relevant spectral features. Considering that the typical drive frequency in RF traps is usually in the MHz regime, micromotion will be reflected either as broadening of the spectral lines, or echoes at the drive frequency, depending on whether the width of the spectral lines is wider or narrower than the drive frequency, respectively. The first case is typical of dipole transitions, and its effects are commonly used to compensate excess micromotion. The procedure consists in driving a fluorescence transition with a detuning equal to the drive frequency and minimizing the fluorescence when varying the voltage of several electrodes. This reduces the micromotion, which in turn sharpens the spectrum (Pruttivarasin et al., 2014). On the other hand, for narrow transitions like quadrupolar or stimulated Raman transitions, the spectral lines are replicated as smaller echoes separated by the drive frequency (Peik et al., 1999; Goham and Britton, 2022).

A mixed behavior between these two regimes can be found in dark resonance spectra, which we study in detail in this article. The atomic spectra of atoms with lambda-type level schemes present narrow negative peaks, known as dark resonances, when the detunings of the lasers that drive the transitions from the lower states to the excited one are similar. This is due to coherent population trapping (CPT) phenomena (Arimondo and Orriols, 1976; Gray et al., 1978; Vanier, 2005). The widths of the dark resonances are typically much smaller than those of dipolar transitions.

For multi-level systems, many dark resonances can emerge depending on the polarizations of the lasers and the geometry of the optics. A typical dark resonance spectrum is shown in Figure 1F, where the polarizations are such that four dark resonances emerge. The widths are related to the product of the two Rabi frequencies and the dephasing mechanisms at work, and reach values of up to a few MHz. For these spectra, the micromotion will manifest as a mixture of further broadening of the already broad main line together with echoes of the dark resonances which appear as dips within the main line.

Figure 1

Figure 1. (A) Frontal view photograph of the ring trap¹, where the electrodes are indicated. At the trap center, fluorescence from trapped ions can be seen. (B) Top view schematic of the trap. The lasers direction, polarization vectors $\widehat{ϵ}$ and the light acquisition system is indicated. (C) Pulsed sequence for measuring dark resonance spectra. An initial Doppler cooling of 100 μs is followed by a measurement of 50 μs. After 20 repeats, a long Doppler cooling of 1 ms is performed to avoid excess heating. The whole sequence is repeated 100 times to increase the collected ion light and thus reduce statistical uncertainties. (D) Relevant level scheme for the Doppler cooling process in ⁴⁰Ca⁺ ion with the two lasers, labelled UV and IR, that excite each of the transitions. (E) CPT allowed transitions under the polarization configuration σ₊ + σ₋ for each laser. This configuration is set using linear and perpendicular polarization to the magnetic field that is used to obtain a Zeeman splitting of the sublevels in each manifold. (F) A typical dark resonance spectrum with the four dark resonances that emerge for the allowed transitions. Each peak correspond to one of the lambda-type systems indicated in (E).

In this paper we present a detailed study of the effects of micromotion in dark resonance spectra in trapped ions. We work with ⁴⁰Ca⁺ ions confined in a 3D ring-shaped Paul trap. By displacing a single ion within the trap we measure its dark resonance spectrum under different micromotion conditions. We show how to retrieve the modulation factor of the ion using a theoretical model to fit the spectra. We also estimate an effective temperature of the ion, evidencing the presence of an additional heating mechanism due to micromotion known as RF heating (Blümel et al., 1989; Kalincev et al., 2021). Moreover, we measure the spectra of more than one ion at the same time and retrieve the modulation factor of each ion through a fit to the model for the full measurement.

The paper is organized as follows. In Section 2 we present the experimental setup and the measurement algorithm as well as the compensation scheme. In Section 3 we show how to calculate dark resonance spectra with and without micromotion. In Section 4 we present results for micromotion-modulated spectra with a single ion and analyze how it alters the final temperature of the ion. Finally, in Section 5 we present results on dark resonance spectra for two ions and study how to model the sum of their spectra.

2 Experimental setup and compensation scheme

We perform our experiments with singly-ionized ⁴⁰Ca atoms trapped in a 3D ring-shaped Paul trap, as shown in Figures 1A, B. The ions sit at the center of the ring-shaped electrode, set to a voltage of ∼ 300 V_pp oscillating at 22.1 MHz which is responsible for the radial trapping. The axial confinement, on the other hand, is produced by two endcaps, top and bottom, with ∼ − 0.5 V each. One extra electrode located below the trap is used to compensate the position of the ion. The oven, which is the source of neutral ⁴⁰Ca atoms, is also used as an additional electrode for compensation. In this work, we will vary the voltage of the top endcap electrode to controllably add micromotion to the ions.

The calcium ion of interest has a multi-lambda type electronic level system, schematized in Figure 1D. The two ground 4²S_1/2 states are connected through a dipole transition near 397 nm (which we label UV for “ultraviolet”) to two excited 4²P_1/2 states, with a lifetime of $\sim 7$ ns (Hettrich et al., 2015). From this P level, the ion has a probability of 0.94 to decay back to the S states, and a probability of 0.06 to decay to one of the four metastable 3²D_3/2 states (Ramm et al., 2013; Barreto et al., 2023a). These metastable states are connected to the excited ones through another dipolar transition near 866 nm (labelled IR for “infrared”). A closed fluorescence cycle, which avoids pumping to dark metastable states, is established through the UV and IR lasers together with an additional magnetic field (Barreto et al., 2023b). The polarizations of the lasers are set to be linear and orthogonal to the magnetic field, as depicted in Figure 1B, which has a magnitude of B ∼ 4 G. Then both lasers have σ₊ + σ₋ polarization in the atomic (magnetic) eigenbasis of the ion. In such a configuration, the atomic structure enables the formation of four lambda-type systems, illustrated in Figure 1E, each of which is responsible for one dark-resonance peak in the spectrum, as shown in Figure 1F.

We measure atomic spectra by collecting scattered ultraviolet light of the ion for different values of IR detuning Δ_IR while keeping the UV detuning Δ_UV fixed. To ensure that the temperature of the ion is the same for all data points, we perform a pulsed algorithm where we interleave cooling and detection (Reiß et al., 2002; Barreto et al., 2022). A scheme of the algorithm is represented in Figure 1C. We start by performing a Doppler cooling stage for 100 μs with the UV laser red-detuned around Δ_UV = −15 MHz and the IR laser blue-detuned around Δ_IR = +20 MHz. With these cooling parameters, which we use for all measurements, we prepare the ion in a thermal state at a temperature close to the Doppler limit for a single well-compensated ion. After this, we switch Δ_UV to −25 MHz and Δ_IR to the measurement value and collect the scattered light of the ion for 50 μs. This time interval is chosen so that the scattering does not affect significantly the temperature of the ion. We repeat this sequence 20 times, and then perform a longer Doppler cooling stage of 1 ms to avoid excess heating of the ion and to remedy heating from other extreme events such as background gas collisions. Each sub-loop is repeated 100 times to reduce statistical uncertainties. Each datapoint then consists of 2000 measurements.

In order to compensate excess micromotion for a single ion we use the DC endcaps, one compensation electrode and the atom oven, located below the trap. For the compensation we use the following algorithm. First, we set Δ_IR = +20 MHz and Δ_UV ∼ − Ω_RF, where Ω_RF is the drive frequency. Then, we manually sweep the voltage of each compensation electrode and look for a minimum of the fluorescence. This reduces the power in the first-order modulation sideband of the UV spectrum, hence finding the position of the ion with minimum micromotion along the propagation direction of the UV laser (Keller et al., 2015). In order to do a full 3D compensation, this procedure should be iterated with three different lasers coming from from three different directions, ideally orthogonal between them, finding a global minimum. For the experiments presented in this paper, we only compensate the ion in two directions.

3 Dark resonance spectra with micromotion

In this section we present measurements of dark resonance spectra including micromotion effects. We then introduce a theoretical model that fits the measurements and subsequently describe how to retrieve relevant experimental parameters with it.

3.1 Experimental spectra

We measure the dark resonance spectra of a single ion for different positions in the ion trap where it is subjected to varying degrees of micromotion. As we vary the potential applied to the top endcap electrode, V_ec, the ion is displaced to different positions with respect to the center of the trap following a straight line, as shown at the sketch in Figure 2G. When the ion is closer to the center we expect it to experience less micromotion than when it is further away. For each position a spectrum was measured and the modulation factor was determined by a fit of the model. Some representative spectra are shown in Figures 2A–F. The one with the least micromotion, when the ion is closest to the center of the trap, corresponds to Figure 2C. There, the four fundamental dark resonances can be identified. Further away from the center, when micromotion is more pronounced, echoes of the four fundamental resonances appear in the spectrum, as can be seen in the other subplots.

Figure 2

Figure 2. Dark resonance spectra with micromotion induced by (de)compensation of the position of the ion by means of an endcap voltage. Subfigures (A–F) depict the acquired spectra (dark-green dots) and the corresponding fitted model (black lines) for selected voltage values. The parameter β of the regression model serves as a measure of micromotion modulation and is expected to be proportional to the distance from the trap center. The ion position for each measurement varies proportionally to the applied voltage on the endcap, along a linear trajectory as illustrated in (G). The values of the factor β as a function of the applied voltage are plotted in (H) as red dots, accompanied by the expected hyperbolic relation fitted to the data (blue line). Residuals of the fit are presented in (I).

For all spectra in Figure 2, the data are plotted along with fits to a model which considers micromotion. In the following subsection we show how to build such a formalism, which adds a modulation at the drive frequency Ω_RF in the atomic spectra. This allows us to retrieve a modulation factor that quantifies the amount of micromotion in the ion, as well as the temperature of the ion.

3.2 Numerical calculation of spectra including micromotion

To calculate the spectra including micromotion effects we closely follow (Barreto et al., 2022) and (Oberst, 1999), in which similar theoretical approaches were developed. To model the atomic dynamics of the full eight-level system interacting with the two lasers, we resort to a master equation of the form

$\frac{d\rho }{dt}=-\frac{i}{\hslash }\left[H,\rho \right]+\mathcal{L}\left(\rho \right).(1)$

Here, ρ is the density matrix of the system, H is the Hamiltonian of the eight-level system interacting with the two lasers under the rotating wave approximation in the frame rotating with the lasers, and $\mathcal{L}$ is a Lindblad superoperator that models the spontaneous emission from the upper levels into the two lower manifolds, as well as finite laser linewidths and thermal motion as described in (46).

To include micromotion in the dynamics, we consider the Doppler shift that the ion experiences due to its periodic motion with respect to the beams. This can be treated by adding the corresponding Doppler shift to the detuning of each laser for the ion at rest, ${\mathrm{\Delta }}_{i_0}$, as

${\mathrm{\Delta }}_i={\mathrm{\Delta }}_{i_0}+{\vec{k}}_i\cdot {\vec{v}}_0\cos \left({\mathrm{\Omega }}_{\text{RF}}t\right),(2)$

where i = {IR,UV}, ${\vec{k}}_i$ is the wave vector of the corresponding laser, Ω_RF is the drive frequency and ${\vec{v}}_0\cos ({\mathrm{\Omega }}_{\text{RF}}t)$ is the instantaneous velocity of the atom, considering only the contribution due to micromotion. The magnitude of ${\vec{v}}_0$ is proportional to the displacement of the ion from the trap center. Motion due to thermal effects is introduced in the following section.

We define a dimensionless modulation factor β for each beam as

${\beta }_i{\mathrm{\Omega }}_{\text{RF}}={\vec{k}}_i\cdot {\vec{v}}_0,(3)$

which quantifies the strength of micromotion effects (Berkeland et al., 1998). It is important to note that the value of β_i depends on the product between the wavevector of the corresponding beam and the driven ion motion. As in our case both beams are co-propagating, we have that 866 nm × β_IR = 397 nm × β_UV. Therefore only one of the β_i needs to be specified as a free parameter in the mathematical model. For the sake of simplicity we refer to β ≡ β_UV and use it as the free parameter that quantifies the micromotion.

The explicit time-dependence of the dynamics is treated through a Floquet-like approach, where we map the master equation to a Liouvillian equation of the form

$\frac{d\vec{\rho }}{dt}=L\left(t\right)\vec{\rho },(4)$

where $\vec{\rho }$ is the vectorized form of ρ and L(t) is a Liouvillian superoperator in the corresponding matrix form. This can in turn be written as

$L\left(t\right)=L_0+\beta {\mathrm{\Omega }}_{\text{RF}}\mathrm{\Delta }L\cos \left({\mathrm{\Omega }}_{\text{RF}}t\right),(5)$

where L₀ is the Liouvillian for the system without modulation, and ΔL is a 8² × 8² matrix with constant elements (Oberst, 1999).

For the trivial case where β = 0, i.e., without micromotion, the steady state can be found by inverting L₀ with a normalization constraint over $\vec{\rho }$. We calculate the atomic spectrum by plotting, as a function of the detuning Δ_IR, the sum of the populations of the two excited P states ${\rho }_{{\text{P}}_{+1/2}{\text{P}}_{+1/2}}+{\rho }_{{\text{P}}_{-1/2}{\text{P}}_{-1/2}}$, which is proportional to the fluorescence of the ion. The plot of Figure 1F was calculated in this way, and exhibits the four expected dark resonances.

If β ≠ 0, the atomic populations will oscillate around a mean value with frequency components that are multiples of Ω_RF. This will be reflected in the spectrum as modulation sidebands displaced from the original dips by integer multiples of Ω_RF. To obtain the time-averaged value of the populations, we propose a solution for times longer than all transient times of the system (Barreto et al., 2023a) with the form

$\vec{\rho }\left(t\right)=\sum _{m=-\infty }^{+\infty }{\vec{\rho }}_m{e}^{im{\mathrm{\Omega }}_{\text{RF}}t}.(6)$

When introducing this ansatz into Eq. 4, we get the recursive relation

$\left(L_0-in{\mathrm{\Omega }}_{\text{RF}}\right){\vec{\rho }}_n+\frac{\beta {\mathrm{\Omega }}_{\text{RF}}}{2}\mathrm{\Delta }L\left({\vec{\rho }}_{n+1}+{\vec{\rho }}_{n-1}\right)=0.(7)$

We are interested in obtaining ${\vec{\rho }}_0$, i.e., the mean value of $\vec{\rho }(t)$. To do so, we define two operators of the form

$S_{n-1}^{+}=-{\left(L_0-in{\mathrm{\Omega }}_{\text{RF}}+\frac{\beta {\mathrm{\Omega }}_{\text{RF}}}{2}\mathrm{\Delta }L{S}_n^{+}\right)}^{-1}\frac{\beta {\mathrm{\Omega }}_{\text{RF}}}{2}\mathrm{\Delta }L,(8)$

$S_{n+1}^{-}=-{\left(L_0-in{\mathrm{\Omega }}_{\text{RF}}+\frac{\beta {\mathrm{\Omega }}_{\text{RF}}}{2}\mathrm{\Delta }L{S}_n^{-}\right)}^{-1}\frac{\beta {\mathrm{\Omega }}_{\text{RF}}}{2}\mathrm{\Delta }L.(9)$

It is straightforward to see that $S_n^{+}{\vec{\rho }}_n={\vec{\rho }}_{n+1}$ and $S_n^{-}{\vec{\rho }}_n={\vec{\rho }}_{n-1}$, so the defined operators act as raising and lowering operators respectively, in the sense that they relate components that oscillate at neighbouring frequencies.

In terms of these operators, we can evaluate Eq. 7 for n = 0 and obtain

$\left(L_0+\frac{\beta {\mathrm{\Omega }}_{\text{RF}}}{2}\mathrm{\Delta }L\left(S_0^{+}+S_0^{-}\right)\right){\vec{\rho }}_0=0.(10)$

Finally, to calculate $S_0^{+/-}$ we truncate the sum in Eq. 6 such that ${\vec{\rho }}_{\pm N_{\mathit{max}}}=0$. This is equivalent to saying that we will include only up to N_max − 1 modulation sidebands in the spectrum. The appropriate value of N_max depends on the value of β: for larger micromotion, a higher number of sidebands will appear, and a higher N_max will be required to properly describe the atomic spectra.

Micromotion can be seen as a frequency modulation of the spectrum, which allows us to use Bessel functions to estimate how many orders to use in the approximation. A frequency modulated signal with modulation factor β will have terms at NΩ_RF with their amplitude given by the Bessel function J_N(β). Therefore we could choose N_max by imposing that the Bessel function J_N(β) is negligible for N > N_max. Nevertheless, we see that setting N_max ∼ β is sufficient for all practical considerations. For all fits in this work, we set N_max = 5, since β never surpassed that value in our experiments.

3.3 Estimation of the modulation parameter

In Figures 2A–F we plot, in solid line, a fit to the model for the six measured spectra, which are shown in green dots along with their uncertainties. The fits are in very good agreement with the measurements, evidencing the robustness of the model to predict how micromotion affects the atomic spectra. From them, we retrieve a value for the modulation factor β for each voltage, as well as the temperature T of the ion as explained in the next section.

In Figure 2H we show how β depends on the endcap voltage V_ec. There is an optimal point at −0.55 (1) V in which β is minimum, but not zero. This is due to extra stray fields not compensated in directions that are orthogonal to the one of the endcaps.

Considering that the ion moves along a trajectory that misses the center of the trap while sweeping V_ec by an amount of y, that its position varies linearly with it, as it is shown in Figure 2G, following (V_ec − V₀)α_V and the proportionality of β with the distance r to the center (β = α_β r), the data can be modeled by a hyperbola function of the form

$\beta \left(V\right)={\alpha }_{\beta }\sqrt{{\alpha }_V^2{\left(V_{\text{ec}}-V_0\right)}^2+y^2}.(11)$

Bringing the α_β inside the square root, the four parameters in the equation above can be grouped and reduced to three in the regression model. The fitted model is plotted in a solid blue line over the experimental data in Figure 2H, showing good agreement with the data. From this fit we extract the minimum value of β which is β_min = 0.30 (5). Furthermore, in Figure 2I we plot the residuals of the fit, showing a random dispersion around zero.

4 Micromotion and thermal effects

4.1 Inclusion of thermal effects in the spectrum

We incorporate the temperature T of the ion in the model following the approximation presented in (46), where the effect is accounted for as an additional dephasing term proportional to the square root of T, which is added as an extra broadening of the laser linewidths

${\mathrm{\Gamma }}_D=\left|{\vec{k}}_{\text{IR}}-{\vec{k}}_{\text{UV}}\right|\sqrt{\frac{k_{B}T}{2m}}.(12)$

One can understand this approximation the following way. The Doppler effect generates a shift proportional to kv. As we are analyzing a two photon transition, we replace k for the difference $\left|{\vec{k}}_{\text{UV}}-{\vec{k}}_{\text{IR}}\right|$. Because we are considering a thermal state, v fluctuates with a standard deviation $\sqrt{k_{B}T/m}$. This is modelled through stochastic shifts in the frequency of the lasers.

This method implies then two important approximations, the first by only considering two-photon processes, and the second by approximating a Gaussian distribution by a Lorentzian one. The $\sqrt{2}$ in the denominator of Eq. 12 was chosen because it was found to provide the best approximation. We note, however, that the proportionality factor leading to the best results is generally dependent on the temperature regime of interest as well as on the relevant linewidths. A more detailed analysis of this is beyond the scope of this paper and will the subject of a future publication. Finally, to include Γ_D in the model, we multiply it by k_i/(k_UV + k_IR), where i = {UV, IR}, and add it to the corresponding linewidth of the laser.

Following this procedure, the ion temperature T becomes a parameter that we can retrieve from the fits. The temperatures estimated from the spectra of Figure 2 are plotted as a function of the modulation factor β in Figure 3, along with their uncertainties derived from the standard error of the fitted parameters in Figure 2. The plot shows a temperature that is approximately constant and close to the Doppler limit for β < 1. For larger micromotion, the temperature exhibits a significant boost upon increasing β. The temperature uncertainty is also bigger for higher T since the dark resonance thermometry method becomes less sensitive to temperature changes (Roßnagel et al., 2015).

Figure 3

Figure 3. Relationship between temperature and micromotion modulation factor β. The data inferred from the spectra in Figure 2 are depicted by red dots along with uncertainty intervals. Two models have been tested to describe this relationship: one excluding RF heating (orange dashed line) and another incorporating it (blue straight line). Further details about each model can be found in the main text. The estimated temperatures for low micromotion are consistent with the Doppler limit for detuning of around −10 MHz, represented by the gray shaded area in the graph.

To understand this behavior, we now characterize the RF heating that the ion experiences and add this effect to the equation describing the thermalization of the ion.

4.2 Thermalization in the presence of micromotion

When the ion is close to the Doppler temperature, the equation determining the energy balance under laser cooling is

$\dot{E}={\dot{E}}_{\text{Dop}}+{\dot{E}}_{\text{rec}}+{\dot{E}}_{\text{RF}},(13)$

where ${\dot{E}}_{\text{Dop}}$ represents the Doppler cooling mechanism due to the interaction of the atom with the lasers, ${\dot{E}}_{\text{rec}}$ is the recoil heating due to the random emission of one photon of the atom, and ${\dot{E}}_{\text{RF}}$ is the RF heating. In the following, we show how to calculate the three terms separately.

4.2.1 Doppler cooling

Doppler cooling works as a dissipation channel (Yan et al., 2019). It acts on the ion through a force of the form

${\vec{F}}_d=\hslash \vec{k}\mathrm{\Gamma }{\rho }_{ee}\left(\mathrm{\Delta },\beta \right),(14)$

where $\vec{k}$ is the wave vector of the Doppler cooling laser, Γ is the linewidth of the cooling transition, and ρ_ee is the excited state population, that is modulated by the micromotion. We will approximate the dynamics by a two-level system corresponding to the S−P transition which is the most relevant for the Doppler cooling. This implies that the temperature of the ion will be mainly determined by the UV laser, so all laser parameters in this subsection will correspond to it.

In the low saturation limit we can express the excited state population using Bessel functions J_n as (DeVoe et al., 1989)

${\rho }_{ee}\left(\mathrm{\Delta },\beta \right)={\left(\frac{\mu E_0}{2\hslash }\right)}^2\sum _n\frac{J_n^2\left(\beta \right)}{{\left(\mathrm{\Delta }+n{\mathrm{\Omega }}_{\text{RF}}\right)}^2+{\left(\mathrm{\Gamma }/2\right)}^2},(15)$

where μ is the electric dipole moment of the transition and E₀ is the electric field of the laser. To include the motion of the atom, we will consider low velocities $\vec{v}$ such that $|\vec{k}\cdot \vec{v}|\ll |\mathrm{\Delta }|$, where $\vec{k}$ is the wave vector of the cooling laser. We can then approximate

${\rho }_{ee}\left(\mathrm{\Delta }+\vec{k}\cdot \vec{v},\beta \right)\sim {\rho }_{ee}\left(\mathrm{\Delta },\beta \right)+\frac{d{\rho }_{ee}\left(\mathrm{\Delta },\beta \right)}{d\mathrm{\Delta }}\left(\vec{k}\cdot \vec{v}\right).(16)$

Considering that $⟨\vec{v}⟩=0$ since the atom is trapped in a harmonic potential, the cooling power turns to be

${\dot{E}}_{\text{Dop}}=\langle {\vec{F}}_d\cdot \vec{v}\rangle =-\hslash |\vec{k}{|}^2\mathrm{\Gamma }\frac{d{\rho }_{ee}\left(\mathrm{\Delta },\beta \right)}{d\mathrm{\Delta }}\langle v_{\widehat{k}}^2\rangle ,(17)$

where $v_{\widehat{k}}$ is the velocity of the ion projected in the direction of $\vec{k}$. Finally, assuming a thermal state we relate the mean value of the squared of the velocity with the temperature as $\frac{1}{2}m⟨v_{\widehat{k}}^2⟩=\frac{1}{2}{k}_{B}T$. In this form, we can write ${\dot{E}}_{\text{Dop}}$ in terms of the temperature of the ion.

4.2.2 Recoil heating

Recoil heating is a mechanism that balances Doppler cooling leading to a non-vanishing temperature limit. It is due to the spontaneous emission of a photon when the atom is excited, which causes a recoil kick of the atom in a random direction. The rate of heating is proportional to the energy of the emitted photon, the decay rate Γ and the emission probability linked to the excited state population ρ_ee, resulting in (Stenholm, 1986)

${\dot{E}}_{\text{rec}}=\frac{{\left(\hslash |\vec{k}|\right)}^2}{m}\mathrm{\Gamma }{\rho }_{ee}\left(\mathrm{\Delta },\beta \right).(18)$

4.2.3 RF heating

The RF heating mechanism is given by electric noise at frequency components of the driving field in Ω_RF ± ω, where ω is a secular frequency of the ion. This causes a periodic ponderomotive noise force at ω that adds extra heating. In (38), the resulting RF heating rate was shown to scale as

${\dot{E}}_{\text{RF}}\sim {\left[\vec{\nabla }\left({\vec{E}}^2\left(\vec{r}\right)\right)\right]}^2\sim {\beta }^2.(19)$

Here, we are using that the electric field $\vec{E}$ is approximately linear in the ion displacement from the trap center, and that the modulation factor β is also proportional to this displacement. We can then express the rate for RF heating as

${\dot{E}}_{\text{RF}}=C_{\text{RF}}{\beta }^2,(20)$

where C_RF is a factor that depends on experimental parameters like the amount of electric noise. Therefore, it is a quantity that must be estimated from the measurements.

4.2.4 Stationary temperature

Considering all the above contributions one can find the stationary value of T by setting $\dot{E}=0$. The expression for the temperature then becomes

$T_0=\frac{\hslash }{k_B}\frac{{\rho }_{ee}\left(\beta ,\mathrm{\Delta }\right)}{\left(\frac{d{\rho }_{ee}\left(\beta ,\mathrm{\Delta }\right)}{d\mathrm{\Delta }}\right)}+\frac{C_{\text{RF}}}{\left(\frac{\hslash |\vec{k}{|}^2\mathrm{\Gamma }{k}_B}{m}\right)}\frac{{\beta }^2}{\left(\frac{d{\rho }_{ee}\left(\beta ,\mathrm{\Delta }\right)}{d\mathrm{\Delta }}\right)e}.(21)$

If we set β = 0, the minimum temperature is obtained for Δ = −Γ/2, and is $T_0=\frac{\hslash \mathrm{\Gamma }}{2k_B}$, which corresponds to the well-known Doppler limit temperature.

In Figure 3 we show a fit to Eq. 21 in solid blue line, leaving C_RF and Δ as fitting parameters. For comparison, we also include in orange dashed line a fit to the model without RF heating, i.e., setting C_RF = 0 and leaving Δ as the only fitting parameter. From this last fit, we obtain Δ = 2π × − 5.1 (8) MHz. However, this model does not properly fit the measurements for β > 1. In particular, neglecting RF heating leads to the prediction of a temperature decrease when β becomes larger than $\sim 2.4$, in clear contradiction with the experiments.

We can see that the inclusion of the RF heating term captures the increasing behavior of T with increasing β, in good agreement with the measurements. From this fit we retrieve Δ = 2π × (−11.7) (Bermudez et al., 2017) MHz, consistent with the detuning used to Doppler cool the ion before the measurements.

On the other hand, we obtain C_RF = 8.4 (32) × 10^–20 J/s. We can rewrite it in terms of a phonon heating rate of one motional mode of the ion (Leibfried et al., 2003). Considering one of the radial modes of the ion, with ω = 2π × 1 MHz, the heating rate becomes

${\dot{n}}_{\text{RF}}=0.20\left(7\right)\times {\beta }^2\mathrm{p}\mathrm{h}/\mathrm{s}.(22)$

This result is in good agreement with other similar measurements of RF heating effects (Keller et al., 2015).

5 Multi-ion spectra

In a final phase of this work, we measure spectra for a system of two ions. The setup used collects the light of all the trapped ions globally on a photo-multiplier tube, not allowing to distinguish each one individually, but allowing fast measurements with low background. The spectra are therefore composed of the light emitted by all ions. As each ion is in a different position in the trap, they will experience different micromotion, and therefore contribute differently to the summed spectrum. We aimed at extracting the modulation factor and temperature of each ion from the global signal.

In Figure 4 we show one typical spectrum for two ions and different model fits. We find that a two-ion model, with an independent modulation factor β for each one, adequately fits the data, but the retrieved parameters are very sensitive in their initial guess. Moreover, we see that a one-ion model, with only one β, can also produce reasonable fits. Also, we studied the behaviour of the fitted modulation factors as a function of the position of the ion crystal in the trap, and found that the results obtained were also very sensitive to the initialization of the parameters.

Figure 4

Figure 4. Dark resonance spectrum acquired by collecting light from two trapped ions. The data are compared with the fitted model for one (red line) and two (dashed-black line) ions. The β obtained in each regression are reported in their respective colors.

We conclude that extracting independent modulation parameters or temperatures for multiple ions from a single signal is an ill-defined problem, with too many free parameters. For an appropriate measurement, independent collection of each ion’s fluorescence would be necessary. This has to be combined with the short collection times (50 μs), which ensures that the temperatures of the ions are marginally affected by the measurement. Such a procedure is rather challenging with electron multiplying charge-coupled devices (EMCCDs) or complementary metal–oxide–semiconductor (CMOS) cameras, but multi-pixel photomultiplier tubes, avalanche photodiode detectors or photon counting cameras could provide useful alternatives.

6 Conclusion

We have analyzed the effects of micromotion in dark-resonance atomic spectra. We presented experimental results obtained with a single trapped calcium ion for which the amplitude of the RF-driven motion was varied applying a controllable voltage to one trap electrode. We introduced a theoretical model that accurately describes the corresponding spectra, and used it to fit our measurements. From these fits, we were able to retrieve the micromotion modulation factor of each curve.

The same fits were used to estimate the temperature of the ion under different micromotion amplitudes. Our results evidence the presence of a micromotion-induced heating mechanism, also known as RF heating. We include this phenomenon in the standard description of Doppler cooling and show that it leads to good predictions of the final temperature of the ion depending on the value of the micromotion modulation parameter.

These methods could be used to study heat transport in trapped ion systems where excess micromotion is unavoidable due to crystal geometry, for example, 2D crystals in linear traps. However, we note that a proper description of the dynamics of several ions requires single-ion detection. Indeed, our studies of the spectra obtained by collecting the light emitted by two ions together indicate that global light collection is not sufficient for a good characterization of the motional state.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

NN: Writing–review and editing, Writing–original draft, Validation, Supervision, Software, Methodology, Investigation, Formal Analysis, Data curation, Conceptualization. MB: Writing–review and editing, Validation, Software, Methodology, Investigation, Formal Analysis. ML: Writing–review and editing, Software, Methodology, Investigation, Formal Analysis, Data curation. CC: Writing–review and editing, Supervision, Software, Methodology, Formal Analysis. CS: Writing–review and editing, Validation, Supervision, Resources, Project administration, Methodology, Investigation, Funding acquisition, Conceptualization.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. NN, MB, ML, and CS acknowledge support for grants PICT 2018-03350, PICT 2019-04349 and PICT2021-I-A-01288 from ANPCyT (Argentina) and grant UBACYT 2018 Mod I - 20020170100616BA from Universidad de Buenos Aires. CC acknowledges funding from grant PICT 2020-SERIEA-00959 from ANPCyT (Argentina).

Acknowledgments

We acknowledge the contribution of M. Drechsler in the building of the ion trap and early experiments. We thank F. Schmidt-Kaler for his invaluable support and generosity, as well as J. P. Paz for his unconditional help in the setting up of the laboratory. Finally, we thank F. Meconi for taking the photo of Figure 1A.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Footnotes

¹Photo taken by Franco Meconi, https://www.instagram.com/terrazaalcosmos/

References

Arimondo, E., and Orriols, G. (1976). Nonabsorbing atomic coherences by coherent two-photon transitions in a three-level optical pumping. Nuovo Cimento Lett. 17, 333–338. doi:10.1007/bf02746514

CrossRef Full Text | Google Scholar

Barreto, N. A. N., Cormick, C., and Schmiegelow, C. T. (2023b). Polarization vs magnetic field: competing eigenbases in laser-driven atoms. arXiv preprint arXiv:2310.18525.

Google Scholar

Barreto, N. A. N., Drechsler, M., and Schmiegelow, C. T. (2022). Three-laser coherent population trapping in a multi-Λ system: theory, experiment, and applications. Phys. Rev. A 106, 053708. doi:10.1103/PhysRevA.106.053708

CrossRef Full Text | Google Scholar

Barreto, N. A. N., Giardino, L. T., Crucianelli, C. J., Bonetto, M., Drechsler, M., and Schmiegelow, C. T. (2023a). Transient fluorescence with a single trapped ion. JOSA B 40, C48. doi:10.1364/JOSAB.482639

CrossRef Full Text | Google Scholar

Berengut, J. C., Budker, D., Delaunay, C., Flambaum, V. V., Frugiuele, C., Fuchs, E., et al. (2018). Probing new long-range interactions by isotope shift spectroscopy. Phys. Rev. Lett. 120, 091801. doi:10.1103/PhysRevLett.120.091801

PubMed Abstract | CrossRef Full Text | Google Scholar

Berkeland, D., Miller, J., Bergquist, J. C., Itano, W. M., and Wineland, D. J. (1998). Minimization of ion micromotion in a Paul trap. J. Appl. Phys. 83, 5025–5033. doi:10.1063/1.367318

CrossRef Full Text | Google Scholar

Bermudez, A., Schindler, P., Monz, T., Blatt, R., and Müller, M. (2017). Micromotion-enabled improvement of quantum logic gates with trapped ions. New J. Phys. 19, 113038. doi:10.1088/1367-2630/aa86eb

CrossRef Full Text | Google Scholar

Blatt, R., and Wineland, D. (2008). Entangled states of trapped atomic ions. Nature 453, 1008–1015. doi:10.1038/nature07125

PubMed Abstract | CrossRef Full Text | Google Scholar

Blümel, R., Kappler, C., Quint, W., and Walther, H. (1989). Chaos and order of laser-cooled ions in a Paul trap. Phys. Rev. A 40, 808–823. doi:10.1103/PhysRevA.40.808

CrossRef Full Text | Google Scholar

Brouard, S., and Plata, J. (2001). Heating of a trapped ion by random fields: the influence of the micromotion. Phys. Rev. A 63, 043402. doi:10.1103/PhysRevA.63.043402

CrossRef Full Text | Google Scholar

Chou, C., Hume, D., Koelemeij, J., Wineland, D. J., and Rosenband, T. (2010). Frequency comparison of two high-accuracy Al+ optical clocks. Phys. Rev. Lett. 104, 070802. doi:10.1103/PhysRevLett.104.070802

PubMed Abstract | CrossRef Full Text | Google Scholar

Chwalla, M., Benhelm, J., Kim, K., Kirchmair, G., Monz, T., Riebe, M., et al. (2009). Absolute frequency measurement of the Ca+ 40 4 s S 1/2 2- 3 d D 5/2 2 clock transition. Phys. Rev. Lett. 102, 023002. doi:10.1103/PhysRevLett.102.023002

PubMed Abstract | CrossRef Full Text | Google Scholar

Cirac, J., Garay, L., Blatt, R., Parkins, A., and Zoller, P. (1994). Laser cooling of trapped ions: the influence of micromotion. Phys. Rev. A 49, 421–432. doi:10.1103/PhysRevA.49.421

PubMed Abstract | CrossRef Full Text | Google Scholar

Cirac, J. I., and Zoller, P. (1995). Quantum computations with cold trapped ions. Phys. Rev. Lett. 74, 4091–4094. doi:10.1103/PhysRevLett.74.4091

PubMed Abstract | CrossRef Full Text | Google Scholar

Counts, I., Hur, J., Craik, D. P. A., Jeon, H., Leung, C., Berengut, J. C., et al. (2020). Evidence for nonlinear isotope shift in Yb+ search for new boson. Phys. Rev. Lett. 125, 123002. doi:10.1103/PhysRevLett.125.123002

PubMed Abstract | CrossRef Full Text | Google Scholar

DeVoe, R., Hoffnagle, J., and Brewer, R. (1989). Role of laser damping in trapped ion crystals. Phys. Rev. A 39, 4362–4365. doi:10.1103/PhysRevA.39.4362

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedenauer, A., Schmitz, H., Glueckert, J. T., Porras, D., and Schätz, T. (2008). Simulating a quantum magnet with trapped ions. Nat. Phys. 4, 757–761. doi:10.1038/nphys1032

CrossRef Full Text | Google Scholar

Gaebler, J., Baldwin, C., Moses, S., Dreiling, J., Figgatt, C., Foss-Feig, M., et al. (2021). Suppression of midcircuit measurement crosstalk errors with micromotion. Phys. Rev. A 104, 062440. doi:10.1103/PhysRevA.104.062440

CrossRef Full Text | Google Scholar

Goham, C. J., and Britton, J. W. (2022). Resolved-sideband micromotion sensing in Yb+ on the 935 nm repump transition. AIP Adv. 12. doi:10.1063/5.0128113

CrossRef Full Text | Google Scholar

Gray, H., Whitley, R., and Stroud, C. (1978). Coherent trapping of atomic populations. Opt. Lett. 3, 218–220. doi:10.1364/OL.3.000218

PubMed Abstract | CrossRef Full Text | Google Scholar

Hettrich, M., Ruster, T., Kaufmann, H., Roos, C., Schmiegelow, C., Schmidt-Kaler, F., et al. (2015). Measurement of dipole matrix elements with a single trapped ion. Phys. Rev. Lett. 115, 143003. doi:10.1103/PhysRevLett.115.143003

PubMed Abstract | CrossRef Full Text | Google Scholar

Higgins, G., Salim, S., Zhang, C., Parke, H., Pokorny, F., and Hennrich, M. (2021). Micromotion minimization using Ramsey interferometry. New J. Phys. 23, 123028. doi:10.1088/1367-2630/ac3db6

CrossRef Full Text | Google Scholar

Huntemann, N., Okhapkin, M., Lipphardt, B., Weyers, S., Tamm, C., and Peik, E. (2012). High-accuracy optical clock based on the octupole transition in Yb+ 171. Phys. Rev. Lett. 108, 090801. doi:10.1103/PhysRevLett.108.090801

PubMed Abstract | CrossRef Full Text | Google Scholar

Ibaraki, Y., Tanaka, U., and Urabe, S. (2011). Detection of parametric resonance of trapped ions for micromotion compensation. Appl. Phys. B 105, 219–223. doi:10.1007/s00340-011-4463-x

CrossRef Full Text | Google Scholar

Kalincev, D., Dreissen, L., Kulosa, A., Yeh, C., Fürst, H., and Mehlstäubler, T. (2021). Motional heating of spatially extended ion crystals. Quantum Sci. Technol. 6, 034003. doi:10.1088/2058-9565/abee99

CrossRef Full Text | Google Scholar

Keller, J., Kalincev, D., Burgermeister, T., Kulosa, A., Didier, A., Nordmann, T., et al. (2019). Probing time dilation in Coulomb crystals in a high-precision ion trap. Phys. Rev. Appl. 11, 011002. doi:10.1103/PhysRevApplied.11.011002

CrossRef Full Text | Google Scholar

Keller, J., Partner, H. L., Burgermeister, T., and Mehlstäubler, T. (2015). Precise determination of micromotion for trapped-ion optical clocks. J. Appl. Phys. 118. doi:10.1063/1.4930037

CrossRef Full Text | Google Scholar

Leibfried, D., Blatt, R., Monroe, C., and Wineland, D. (2003). Quantum dynamics of single trapped ions. Rev. Mod. Phys. 75, 281–324. doi:10.1103/RevModPhys.75.281

CrossRef Full Text | Google Scholar

Ludlow, A. D., Boyd, M. M., Ye, J., Peik, E., and Schmidt, P. O. (2015). Optical atomic clocks. Rev. Mod. Phys. 87, 637–701. doi:10.1103/RevModPhys.87.637

CrossRef Full Text | Google Scholar

Neuhauser, W., Hohenstatt, M., Toschek, P., and Dehmelt, H. (1980). Localized visible Ba+ mono-ion oscillator. Phys. Rev. A 22, 1137–1140. doi:10.1103/PhysRevA.22.1137

CrossRef Full Text | Google Scholar

Oberst, H. (1999). Resonance fluorescence of single barium ions (master’s thesis).

Google Scholar

Paul, W. (1990). Electromagnetic traps for charged and neutral particles. Rev. Mod. Phys. 62, 531–540. doi:10.1103/RevModPhys.62.531

CrossRef Full Text | Google Scholar

Peik, E., Abel, J., Becker, T., Von Zanthier, J., and Walther, H. (1999). Sideband cooling of ions in radio-frequency traps. Phys. Rev. A 60, 439–449. doi:10.1103/PhysRevA.60.439

CrossRef Full Text | Google Scholar

Porras, D., and Cirac, J. I. (2004). Effective quantum spin systems with trapped ions. Phys. Rev. Lett. 92, 207901. doi:10.1103/physrevlett.92.207901

PubMed Abstract | CrossRef Full Text | Google Scholar

Pruttivarasin, T., Ramm, M., and Häffner, H. (2014). Direct spectroscopy of the 2S1/2- 2P1/2 and 2D3/2- 2P1/2 transitions and observation of micromotion modulated spectra in trapped 40Ca+. J. Phys. B Atomic, Mol. Opt. Phys. 47, 135002. doi:10.1088/0953-4075/47/13/135002

CrossRef Full Text | Google Scholar

Ramm, M., Pruttivarasin, T., Kokish, M., Talukdar, I., and Häffner, H. (2013). Precision measurement method for branching fractions of excited P 1/2 states applied to Ca+ 40. Phys. Rev. Lett. 111, 023004. doi:10.1103/PhysRevLett.111.023004

PubMed Abstract | CrossRef Full Text | Google Scholar

Ratcliffe, A. K., Oberg, L. M., and Hope, J. J. (2020). Micromotion-enhanced fast entangling gates for trapped-ion quantum computing. Phys. Rev. A 101, 052332. doi:10.1103/PhysRevA.101.052332

CrossRef Full Text | Google Scholar

Reiß, D., Abich, K., Neuhauser, W., Wunderlich, C., and Toschek, P. (2002). Raman cooling and heating of two trapped Ba+ ions. Phys. Rev. A 65, 053401. doi:10.1103/PhysRevA.65.053401

CrossRef Full Text | Google Scholar

Rosenband, T., Hume, D., Schmidt, P., Chou, C. W., Brusch, A., Lorini, L., et al. (2008). Frequency ratio of Al+ and Hg+ single-ion optical clocks; metrology at the 17th decimal place. Science 319, 1808–1812. doi:10.1126/science.1154622

PubMed Abstract | CrossRef Full Text | Google Scholar

Roßnagel, J., Tolazzi, K. N., Schmidt-Kaler, F., and Singer, K. (2015). Fast thermometry for trapped ions using dark resonances. New J. Phys. 17, 045004. doi:10.1088/1367-2630/17/4/045004

CrossRef Full Text | Google Scholar

Safronova, M., Budker, D., DeMille, D., Kimball, D. F. J., Derevianko, A., and Clark, C. W. (2018). Search for new physics with atoms and molecules. Rev. Mod. Phys. 90, 025008. doi:10.1103/RevModPhys.90.025008

CrossRef Full Text | Google Scholar

Shaniv, R., Ozeri, R., Safronova, M., Porsev, S., Dzuba, V., Flambaum, V., et al. (2018). New methods for testing Lorentz invariance with atomic systems. Phys. Rev. Lett. 120, 103202. doi:10.1103/PhysRevLett.120.103202

PubMed Abstract | CrossRef Full Text | Google Scholar

Sikorsky, T., Meir, Z., Akerman, N., Ben-Shlomi, R., and Ozeri, R. (2017). Doppler cooling thermometry of a multilevel ion in the presence of micromotion. Phys. Rev. A 96, 012519. doi:10.1103/PhysRevA.96.012519

CrossRef Full Text | Google Scholar

Stenholm, S. (1986). The semiclassical theory of laser cooling. Rev. Mod. Phys. 58, 699–739. doi:10.1103/revmodphys.58.699

CrossRef Full Text | Google Scholar

Turchette, Q. A., King, B., Leibfried, D., Meekhof, D., Myatt, C., Rowe, M., et al. (2000). Heating of trapped ions from the quantum ground state. Phys. Rev. A 61, 063418. doi:10.1103/PhysRevA.61.063418

CrossRef Full Text | Google Scholar

Vanier, J. (2005). Atomic clocks based on coherent population trapping: a review. Appl. Phys. B 81, 421–442. doi:10.1007/s00340-005-1905-3

CrossRef Full Text | Google Scholar

van Mourik, M. W., Hrmo, P., Gerster, L., Wilhelm, B., Blatt, R., Schindler, P., et al. (2022). rf-induced heating dynamics of noncrystallized trapped ions. Phys. Rev. A 105, 033101. doi:10.1103/PhysRevA.105.033101

CrossRef Full Text | Google Scholar

Yan, L. L., Su, S. L., and Feng, M. (2019). Analytical investigation of one-dimensional Doppler cooling of trapped ions with Λ-type configuration. Phys. Rev. A 100, 033418. doi:10.1103/PhysRevA.100.033418

CrossRef Full Text | Google Scholar