## Three-dimensional trapping of individual Rydberg atoms in ponderomotive bottle beam traps

D. Barredo, V. Lienhard, P. Scholl, S. de L´ es´ eleuc, ∗ T. Boulier, A. Browaeys, and T. Lahaye Laboratoire Charles Fabry, Institut d'Optique Graduate School, CNRS, Universit´ e Paris-Saclay, 91127 Palaiseau Cedex, France

We demonstrate three-dimensional trapping of individual Rydberg atoms in holographic optical bottle beam traps. Starting with cold, ground-state 87 Rb atoms held in standard optical tweezers, we excite them to nS 1 / 2 , nP 1 / 2 , or nD 3 / 2 Rydberg states and transfer them to a hollow trap at 850 nm. For principal quantum numbers 60 /lessorequalslant n /lessorequalslant 90 , the measured trapping time coincides with the Rydberg state lifetime in a 300 K environment. We show that these traps are compatible with quantum information and simulation tasks by performing single qubit microwave Rabi flopping, as well as by measuring the interaction-induced, coherent spin-exchange dynamics between two trapped Rydberg atoms separated by 40 µ m. These results will find applications in the realization of high-fidelity quantum simulations and quantum logic operations with Rydberg atoms.

Neutral atoms excited to Rydberg states are an attractive platform for large-scale quantum simulation and computation [1, 2]. The strong, controllable interactions between these excited states can be used to implement high-fidelity quantum logic gates, or to engineer various types of spin Hamiltonians that are difficult to study on classical computers [3]. These ideas have been intensively explored in the last years and several important milestones have been achieved [4]. Prominent examples of this progress are the demonstration of strong optical non-linearities [5], single-photon sources [6], conditional phase shifters [7], single-photon transistors [8, 9], the experimental realizations of two-qubit gates [10-13], or the first quantum simulations of spin models with tens of particles in optical lattices [14-16] and in arrays of optical tweezers [1720].

far only in one dimension. Here we go beyond those initial demonstrations to show efficient three-dimensional trapping of cold individual Rydberg atoms in micron-size optical potentials [34, 35]. We use holography to create bottle beam (BoB) traps [36, 37] which are deterministically loaded with single Rydberg atoms. We characterize the depth and trapping frequencies of these traps and observe that the trapping time for principal quantum numbers in the range 60 /lessorequalslant n /lessorequalslant 90 is mainly limited by the lifetime of the Rydberg states in the presence of blackbody radiation at 300 K. Finally, we illustrate the compatibility of these traps with quantum simulation by driving Rabi flopping between different Rydberg states,

In none of the above experiments were the Rydberg atoms trapped. However, control over the motion of Rydberg atoms during gate operation and in quantum simulations is advantageous, since finite atom temperatures and mechanical forces due to the strong interactions between the particles ultimately limit quantum state fidelities [4, 21] and the available time for coherent dynamics [18, 22, 23]. Rydberg trapping is also a prerequisite for precision measurements of fundamental constants using circular Rydberg states [24, 25] or positronium [26].

To date, clouds of Rydberg atoms have been confined to three-dimensional, millimeter-size regions using static magnetic [27-29] or electric fields [30, 31]. In an inhomogeneous AC electric field that oscillates faster than any internal frequency of the Rydberg atom, the weakly-bound Rydberg electron experiences an oscillating force that can be used for laser trapping of the Rydberg atom [32]. The so-called ponderomotive potential, which is the time-averaged kinetic energy of the nearly free Rydberg electron oscillating in the laser field, is proportional to the light intensity. Therefore, to obtain a 3D trap, one must create a dark region surrounded by light in all directions; since the atom trapping arises mainly from the ponderomotive potential experienced by the electron, such traps can be used to confine Rydberg states whatever their n, l, j, m j quantum numbers.

Rydberg atoms have been efficiently confined in the tight potentials of ponderomotive optical lattices [32, 33], but so

FIG. 1. Main components of the experimental setup. We use a spatial light modulator (SLM) and a high NA aspherical lens to generate a BoB trap for Rydberg atoms (blue beam, see text). The zoom shows two-dimensional cuts of the reconstructed light intensity distribution near the focal plane, measured using an electrically-tunable lens (not shown) to scan the images along the optical axis. To load the BoB trap, we use regular optical tweezers (red beam) created by light reflected from a second SLM and superimposed onto the BoB trap beam using a polarizing beam splitter (PBS).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-1c19e889-96eb-4e16-aa10-d4a4369aa0aa.jpg)

and by observing the coherent exchange of internal states induced by the dipole-dipole interaction for two atoms confined in BoB traps separated by 40 µ m.

Far from resonances and in a Born-Oppenheimer-like approximation, the Rydberg atom trapping potential in a light field of intensity I and angular frequency ω L can be described as [40]:

The experimental setup is based on the one described in Ref. [38]. Briefly, we use a spatial light modulator (SLM) to create arbitrary arrays of micron-size optical tweezers (Fig. 1). To generate the BoB traps that will host Rydberg atoms, we follow the procedure of Refs. [37, 39]. We use a second SLM to imprint on another light beam at 850 nm a phase pattern that is composed of two terms added modulo 2 π : (i) the standard hologram to create the desired number and positions of point traps, and (ii) a centered disk of radius r 0 where the phase is shifted by π . The value of r 0 is adjusted such that, for each trap, destructive interference occurs between the central and outer parts of the Gaussian beam at the focal point, thus creating the needed dark region surrounded by light in all directions. The local maximum of light intensity encountered when moving away from the origin is smallest ( ∼ 10% of the maximal intensity) on two 'escape cones' originating from the trap center. By changing the diameter of the beam impinging on the SLM and adjusting r 0 accordingly, the size of the trapping region can be tuned. An example of BoB trap created using this method is shown in the inset of Fig. 1, where we used an electrically-tunable lens in the imaging path to record the 3D light intensity distribution near the focal plane of the aspherical lens. To allow for an efficient transfer of the atoms between the two different traps (for ground-state and Rydberg atoms), the light from the BoB SLM is then overlapped with the standard optical tweezers beam using a polarizing beam splitter (PBS), ensuring that both traps are centered at the same position.

$$& \text{as } \lfloor + \upsilon \rfloor. & \text{atom} \\ & U _ { n l j m _ { j } } ( R ) = \int d ^ { 3 } r \, V _ { P } ( R + r ) \, | \psi _ { n l j m _ { j } } ( r ) | ^ { 2 } \,. \quad ( 1 ) \quad \text{optic} \\ & \text{deper} \\ & \text{at} \, \upsilon _ { - } \, \upsilon ^ { 2 } \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upseta \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsery \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsery \, \upsilon \, \upsilon \, \upsilon \, \upsery \, \upsilon \, \upsery \, \upsilon \, \upsery \, \upsilon \, \upsery \, \upsilon \, \upsery \, \upsilon \, \upsery \, \upsery \, \upsilon \, \upsery \, \upsery \, \upsilon \, \upsery \, \upsery \, \upsilon \, \upsery \, \upsilon \, \upsery \, \upsery \, \upsilon \, \upsery \, \upsery \, \upsery \, \upsery \, \upsilon \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upseries \, \text{at} \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsigma \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsigma \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsigma \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upset \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upsilon \, \upser \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsry \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsERY \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsEvery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \ups EVERY \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsory \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsevery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsvery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upstery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsrying \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsrys \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \upsery \, \ups$$

Here, V P ( r ) = e I 2 ( r ) / (2 m /epsilon1 e 0 cω 2 L ) is the repulsive ponderomotive shift experienced by the nearly free electron of mass m e and charge e , ψ is the Rydberg wave function, R is the center of mass coordinate of the atom, and r is the relative position of the Rydberg electron.

An example of such calculation for the 84 S 1 / 2 Rydberg state is shown in the axial cut displayed in Fig. 2b, calculated from the measured light intensity displayed on Fig. 2a. The trapping potential is harmonic in the axial direction z , whereas in the transverse directions it has approximately a quartic form close to the center, and maxima separated by ∼ 2 5 . µ m. From the measured three-dimensional light intensity distribution and a total power of 400 mW (Fig. 2a), we obtain minimum trap barriers (saddle point along the escape cones) of around 0.6 mK (Fig. 2b). When the variation of the field starts to be substantial over a length scale comparable to the atom size, the trapping potential U nljm j ( R ) depends strongly on the specific Rydberg state. This is illustrated in

FIG. 2. (a) Example of measured light intensity distribution near the focal plane of the aspherical lens for a total beam power of 400 mW. (b) Trapping potential for the 84 S 1 / 2 Rydberg state calculated from Eq. (1) and the intensity distribution shown in (a). (c) Onedimensional cuts of the trapping potential along z = 0 (dashed line in (b)) for different nS 1 / 2 Rydberg states. A schematic of the 60 S and 120 S Rydberg orbitals (to scale) is shown as an inset.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-022ba2e3-de0b-41a9-b842-b43aae25f9a2.jpg)

the one-dimensional cuts along the transverse direction for different Rydberg states displayed in Fig. 2c. When the orbital radius ∼ n a 2 0 of the atom becomes comparable to the trap size (here, for n ∼ 120 ), the potential does not have a local minimum any longer, and Rydberg atoms are then repelled from the BoB trap.

The experimental sequence is shown in Fig. 3. First, 87 Rb atoms from a magneto-optical trap (MOT) are loaded in the optical tweezers (with temperatures in the range 3 -130 µ K depending on the experiment). We detect the occupancy of each trap by collecting the fluorescence of the atoms at 780 nm with an electron-multiplying CCD camera in 20 ms. After successful detection of a filled trap, atoms are optically pumped in the ground state | g 〉 = 5 ∣ ∣ S 1 / 2 , F = 2 , m F = 2 〉 , the optical tweezers are switched off, and a two-photon STImulated Raman Adiabatic Passage (STIRAP) prepares the target Rydberg state | r 〉 with ∼ 90% efficiency. We then switch on the BoB trap for a variable duration τ . Finally, a Rydberg de-excitation pulse, resonant with the | r 〉 -5 P 1 / 2 transition, brings the population in | r 〉 back to the ground state | g 〉 via spontaneous decay from 5 P 1 / 2 , and the atom is then imaged by fluorescence. The result of one of these experiments is plotted in Fig. 3. Without the BoB trap (red curve), this release and recapture experiment can be used to measure single atom temperatures [41]. The recapture probability for atoms in 84 S 1 / 2 is reduced by a factor ten after 30 µ s. From this curve we extract a single atom temperature of ∼ 130 µ K. When the BoB trap is applied (blue curve), the recapture probability is drastically enhanced, and only 30 % of the atoms are

FIG. 3. (a) Time sequence of the experiment. Ground-state atoms are initially loaded in optical tweezers. The tweezers are then switched off for Rydberg excitation. After 2 µ s the BoB trap is applied for a variable duration τ . Following Rydberg de-excitation in free-flight, the tweezers are switched on again and the atoms are imaged. As the BoB trap is strongly repulsive for ground state atoms, only atoms that were in the Rydberg state before de-excitation are recaptured. (b) Recapture probability for atoms excited to the 84 S 1 / 2 state as a function of the trapping time τ (disks). The recapture curve in the absence of BoB trap (squares) is shown for comparison. Error bars are the standard error of the mean.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-cb979c14-e8a8-4087-949b-5bec4d572865.jpg)

lost in the same amount of time, consistent with the expected Rydberg lifetime (see below). We do not observe any appreciable heating of the atoms after they are transferred back to the ground-state tweezers.

Following a procedure similar to the one in Ref. [42], we determined the typical radial oscillation frequency of the Rydberg atom in the trap. After the atoms are loaded into the BoB trap, we excite breathing oscillations by switching off the trap for three microseconds. During this trap-off time the atom leaves the bottom of the potential. We then switch on the trap again for a variable time T , during which the atom oscillates in the trap with a typical frequency ω BoB . The recapture probability after a fixed switch off time oscillates with a frequency 2 ω BoB . For atoms excited to the 84 S 1 / 2 state we observe oscillations in the recapture probability that are well reproduced by Monte-Carlo simulations taking into account the experimental parameters, where the sampled atom trajectories inside the trapping field are computed classically [43]. For a trap depth of ∼ 0 6 . mK the atoms oscillate in the traps with a measured frequency of ω BoB / (2 π ) = 59 ± 3 kHz, in fair agreement with the simulations, which predict 59 kHz for the same parameters. The necessary power for efficient trapping depends on the trap size, the initial temperature of the atoms, and the target Rydberg state. For our smallest trap (with maxima separated by ∼ 2 µ m in the radial direction), and single atom temperatures of ∼ 3 µ K, we need only 20 mW of power to reach our higher trapping efficiencies for the 60 S 1 / 2 Rydberg state, which shows that it will be possible to scale this technique up to large arrays of BoB traps.

FIG. 4. (a) Recapture probability of a 84 S 1 / 2 Rydberg atom, as a function of the time spent in the BoB trap, showing a roughly exponential decay with lifetime 222 ± 3 µ s (dashed line). The solid line is the result of a simulation without any adjustable parameter (see text). (b) Measured lifetimes of various nS 1 / 2 states; the dashed line is the theoretical value [44].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-497c4068-8072-4597-9045-6e8dbf32b5f4.jpg)

To evaluate the quality of the BoB potential as a trap for Rydberg atoms, we measured the trap-decay times for different Rydberg states. The recapture probability as a function of the time spent in the BoB trap is shown for the 84 S 1 / 2 state in Figure 4a. An exponential fit to the experimental data (dashed line) gives a 1 /e decay time of 222 ± 3 µ s, in excellent agreement with the calculated lifetime 228 µ s of this state at 300 K [44-46]. The solid line is the result of a Monte-Carlo simulation of the atomic trajectories in the BoB trap potential (calculated using Eq. (1) and the measured light intensity distribution) and taking into account not only radiative decay, but also the possibility for some energetic atoms to spill over the trap barrier. All the parameters entering the simulation were given their experimentally measured values [43]. Figure 4b shows the results of the same experiment for other nS 1 / 2 states; the measured decay times match almost perfectly the expected lifetime (dashed line) up to n ∼ 90 . This indicates that atom losses are mainly due to Rydberg decay, and that other processes, such as photoionization [21, 34, 47], are negligible. We also performed the same measurement for the 84 P 1 / 2 and 82 D 3 / 2 states (not shown) and obtained lifetimes ∼ 20% shorter than their theoretical values.

Beyond n = 90 , the extension of the radial wave function becomes comparable to the trap size and trapping is less efficient. As a consequence, some atoms can escape the trap, and the loss rate is significantly higher than expected from the Rydberg lifetime alone. Beyond n ∼ 100 the potential does not have a local minimum and the Rydberg atoms are quickly expelled from the BoB traps. Monte-Carlo simulations [43] show that for n = 60 , atoms can be lost only through the narrow escape cones of the BoB trap, while for n = 100 , we observe much more pronounced losses over the radial barrier due to the small depth of the trap.

We now investigate the performance of our BoB trap for quantum simulation tasks by studying its compatibility with one- and two-qubit operations. For low principal quantum numbers ( n &lt; 60 ), the trapping potential corresponds approximately to the one of a free electron, and it is there-

FIG. 5. (a) Rabi oscillation between the states |↑〉 = ∣ ∣ 82 D 3 / 2 , m j = 3 2 / 〉 and |↓〉 = ∣ ∣ 83 P 1 / 2 , m j = 1 2 / 〉 . The total trapping time is 50 µ s. (b) Excitation-hopping oscillations in the population of the pair states |↑↓〉 and |↓↑〉 , mediated by the dipole-dipole interaction between the Rydberg states |↑〉 and |↓〉 at a distance of 40 µ m. The temperature of the atoms is ∼ 3 µ K. Solid lines are damped sine fits to the data. Error bars represent the standard error of the mean and are most often smaller than the symbol size.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-7437b0b4-1338-42dc-a8fc-c4f4a3d36ed2.jpg)

fore almost independent of n or the atom's angular momenta. For higher Rydberg states the wave function of the electron must be considered in the convolution (1). Most quantum simulation experiments, however, involve only couplings to adjacent Rydberg levels, for which changes in the squared radial part of the wave function play only a minor role in the effective potential. The angular part, in turn, is identical for states with different quantum number l , but same j , m j (e.g., for ∣ ∣ nS 1 / 2 , m j = 1 2 / 〉 and ∣ ∣ nP 1 / 2 , m j = 1 2 / 〉 ). This 'quasi-magic' trapping condition strongly suppresses differential light shifts in Rabi oscillations between different trapped Rydberg states. In a first experiment, we analyze the coherence in the spin manipulation of a single trapped atom. We apply the BoB trap for a total time of 50 µ s. We let the atoms move in the trapping potential during ∼ 35 µ s and then we drive microwave Rabi oscillations between the states |↑〉 = 82 ∣ ∣ D 3 / 2 , m j = 3 2 / 〉 and |↓〉 = 83 ∣ ∣ P 1 / 2 , m j = 1 2 / 〉 , while the atoms are still confined (Fig. 5a). We observe spin flip oscillations without appreciable damping, despite the two states having different angular wave functions. The constant finite contrast of the oscillations is mainly due to the excitation efficiency and the limited lifetime of the involved Rydberg states. We measured similar Rabi flopping curves even for atom temperatures as high as 130 µ K, where the atoms explore a large volume of the BoB trap.

As a final illustration of the usefulness of our trapping scheme for quantum simulation, we measure spin-exchange dynamics driven by dipole-dipole interaction between two atoms. We create two traps separated by a distance of R =

40 µ m(Fig. 5b). Immediately after loading the Rydberg atoms in the BoB traps, we use a resonant microwave field and local addressing [48] to prepare the first atom in |↑〉 and the second one in |↓〉 . In these two Rydberg states, the atoms are resonantly coupled by a dipole-dipole interaction of strength U dd /h = C /R 3 3 ∼ 0 4 . MHz, and the dynamics is governed by an XY-spin Hamiltonian H = U dd ( σ + 1 σ -2 + σ -1 σ + 2 ), where σ ± i are the Pauli matrices for atom i = 1 2 , . We observe coherent spin exchange between the pair states |↑↓〉 and |↓↑〉 . This experiment demonstrates the compatibility of our Rydberg trapping scheme with quantum simulations using current experimental setups, with levels of contrast and damping comparable or better that previously achieved [22, 48].

In conclusion, we have demonstrated laser trapping of individual Rydberg atoms in microscopic potentials and shown the suitability of such traps for quantum information tasks. This trapping scheme can be extended to larger arrays with moderate laser power for atoms which can be cooled to the motional ground state before Rydberg excitation [49, 50]. In addition, the ponderomotive traps can, in principle, be used for circular Rydberg atoms in a cryogenic environment, opening the door to unprecedented trapping times [51]. The results presented here are suited for experiments that occur entirely in the Rydberg manifold. For quantum logic gates or for quantum simulation involving also the ground state in addition to Rydberg states, it would be possible, by using shorter wavelength trapping light, to realize BoB traps that satisfy the ground-Rydberg 'magic' condition and would minimize heating and decoherence rates [34, 52].

Note: Two-dimensional trapping of circular Rydberg atoms in the ponderomotive potential of a hollow laser beam has been observed recently at Coll` ege de France [53].

This work benefited from financial support by the EU (FETFlag 817482, PASQUANS), by the R´ egion ˆ Ile-de-France in the framework of DIM SIRTEQ (project CARAQUES), and by the IXCORE-Fondation pour la Recherche. T. B. acknowledges the support of the European Marie Skłodowska-Curie Actions (H2020-MSCA-IF-2015 Grant No. 701034).

∗ Present address: Institute for Molecular Science, National Institutes of Natural Sciences, Myodaiji, Okazaki 444-8585, Japan.

- [2] M. Saffman, T. G. Walker, and K. Mølmer, Rev. Mod. Phys. 82 , 2313 (2010).
- [1] D. Jaksch, J. I. Cirac, P. Zoller, S. L. Rolston, R. Cˆ ot´, e and M. D. Lukin, Phys. Rev. Lett. 85 , 2208 (2000).
- [3] H. Weimer, M. Muller, I. Lesanovsky, P. Zoller, and H. P. B¨ uchler, Nature Phys. 6 , 382 (2010).
- [5] J. D. Pritchard, D. Maxwell, A. Gauguet, K. J. Weatherill, M. P. A. Jones, and C. S. Adams, Phys. Rev. Lett. 105 , 193603 (2010).
- [4] M. Saffman, J. Phys. B: At. Mol. Opt. Phys. 49 , 202001 (2016).
- [6] Y. O. Dudin and A. Kuzmich, Science 336 , 887 (2012).
- [8] D. Tiarks, S. Baur, K. Schneider, S. D¨rr, u and G. Rempe, Phys.
- [7] O. Firstenberg, T. Peyronel, Q.-Y. Liang, A. V. Gorshkov, M. D. Lukin, and V. Vuleti´ c, Nature 502 , 71 (2013).

- Rev. Lett. 113 , 053602 (2014).
- [10] T. Wilk, A. Ga¨ etan, C. Evellin, J. Wolters, Y . Miroshnychenko, P. Grangier, and A. Browaeys, Phys. Rev. Lett. 104 , 010502 (2010).
- [9] H. Gorniaczyk, C. Tresp, J. Schmidt, H. Fedder, and S. Hofferberth, Phys. Rev. Lett. 113 , 053601 (2014).
- [11] L. Isenhower, E. Urban, X. L. Zhang, A. T. Gill, T. Henage, T. A. Johnson, T. G. Walker, and M. Saffman, Phys. Rev. Lett. 104 , 010503 (2010).
- [13] H. Levine, A. Keesling, A. Omran, H. Bernien, S. Schwartz, A. S. Zibrov, M. Endres, M. Greiner, V. Vuleti´ c, and M. D. Lukin, Phys. Rev. Lett. 121 , 123603 (2018).
- [12] Y.-Y. Jau, A. M. Hankin, T. Keating, I. H. Deutsch, and G. W. Biedermann, Nature Physics 12 , 71 (2016).
- [14] P. Schauß, M. Cheneau, M. Endres, T. Fukuhara, S. Hild, A. Omran, T. Pohl, C. Gross, S. Kuhr, and I. Bloch, Nature 491 , 87 (2012).
- [16] J. Zeiher, J.-y. Choi, A. Rubio-Abadal, T. Pohl, R. van Bijnen, I. Bloch, and C. Gross, Phys. Rev. X 7 , 041063 (2017).
- [15] P. Schauß, J. Zeiher, T. Fukuhara, S. Hild, M. Cheneau, T. Macr` ı, T. Pohl, I. Bloch, and C. Gross, Science 347 , 1455 (2015).
- [17] H. Labuhn, D. Barredo, S. Ravets, S. de L´ es´ eleuc, T. Macr` ı, T. Lahaye, and A. Browaeys, Nature 534 , 667 (2016).
- [19] H. Kim, Y. Park, K. Kim, H.-S. Sim, and J. Ahn, Phys. Rev. Lett. 120 , 180502 (2018).
- [18] H. Bernien, S. Schwartz, A. Keesling, H. Levine, A. Omran, H. Pichler, S. Choi, A. S. Zibrov, M. Endres, M. Greiner, V. Vuleti´ c, and M. D. Lukin, Nature 551 , 579 (2017).
- [20] S. de L´ es´leuc, V . Lienhard, P. Scholl, D. Barredo, S. Weber, e N. Lang, H. P. B¨ uchler, T. Lahaye, and A. Browaeys, Science (2019), 10.1126/science.aav9105.
- [22] D. Barredo, H. Labuhn, S. Ravets, T. Lahaye, A. Browaeys, and C. S. Adams, Phys. Rev. Lett. 114 , 113002 (2015).
- [21] M. Saffman and T. G. Walker, Phys. Rev. A 72 , 022347 (2005).
- [23] S. de L´ es´leuc, e D. Barredo, V. Lienhard, A. Browaeys, and T. Lahaye, Phys. Rev. A 97 , 053803 (2018).
- [25] U. D. Jentschura, P. J. Mohr, J. N. Tan, and B. J. Wundt, Phys. Rev. Lett. 100 , 160404 (2008).
- [24] A. Ramos, K. Moore, and G. Raithel, Phys. Rev. A 96 , 032513 (2017).
- [26] D. B. Cassidy, Eur. Phys. J. D 72 , 53 (2018).
- [28] D. A. Anderson, A. Schwarzkopf, R. E. Sapiro, and G. Raithel, Phys. Rev. A 88 , 031401 (2013).
- [27] J.-H. Choi, J. R. Guest, A. P. Povilus, E. Hansis, and G. Raithel, Phys. Rev. Lett. 95 , 243001 (2005).
- [29] A. D. Bounds, N. C. Jackson, R. K. Hanley, R. Faoro, E. M.
- Bridge, P. Huillery, and M. P. A. Jones, Phys. Rev. Lett. 120 , 183401 (2018).
- [31] S. D. Hogan and F. Merkt, Phys. Rev. Lett. 100 , 043001 (2008).
- [30] E. Vliegen, S. D. Hogan, H. Schmutz, and F. Merkt, Phys. Rev. A 76 , 023405 (2007).
- [32] S. E. Anderson, K. C. Younge, and G. Raithel, Phys. Rev. Lett. 107 , 263001 (2011).
- [34] S. Zhang, F. Robicheaux, and M. Saffman, Phys. Rev. A 84 , 043408 (2011).

[33]

L. Li, Y. O. Dudin, and A. Kuzmich, Nature

498

,

466 (2013).

- [35] S. Zhang, PhD thesis, University of Wisconsin-Madison (2012).

R1750 (1999).

[36]

R. Ozeri, L. Khaykovich, and N. Davidson, Phys. Rev. A

59

,

- [37] P. Xu, X. He, J. Wang, and M. Zhan, Opt. Lett. 35 , 2164 (2010).
- [39] J. L. Chaloupka, Y. Fisher, T. J. Kessler, and D. D. Meyerhofer, Opt. Lett. 22 , 1021 (1997).
- [38] D. Barredo, V. Lienhard, S. de L´ es´leuc, e T. Lahaye, and A. Browaeys, Nature 561 , 79 (2018).
- [40] S. K. Dutta, J. R. Guest, D. Feldbaum, A. Walz-Flannigan, and G. Raithel, Phys. Rev. Lett. 85 , 5551 (2000).
- [42] Y. R. P. Sortais, H. Marion, C. Tuchendler, A. M. Lance, M. Lamare, P. Fournet, C. Armellin, R. Mercier, G. Messin, A. Browaeys, and P. Grangier, Phys. Rev. A 75 , 013406 (2007).
- [41] C. Tuchendler, A. M. Lance, A. Browaeys, Y. R. P. Sortais, and P. Grangier, Phys. Rev. A 78 , 033425 (2008).
- [43] See the Supplemental Material for more details.
- [45] N. Sibali´ c, ˇ J. Pritchard, C. Adams, and K. Weatherill, Computer Physics Communications 220 , 319 (2017).
- [44] I. I. Beterov, I. I. Ryabtsev, D. B. Tretyakov, and V. M. Entin, Phys. Rev. A 79 , 052504 (2009).
- [46] M. Archimi, C. Simonelli, L. Di Virgilio, A. Greco, M. Ceccanti, E. Arimondo, D. Ciampini, I. Ryabtsev, I. Beterov, and O. Morsch, arXiv:1907.01254 (2019).
- [48] S. de L´ es´leuc, e D. Barredo, V. Lienhard, A. Browaeys, and T. Lahaye, Phys. Rev. Lett. 119 , 053202 (2017).
- [47] R. M. Potvliege and C. S. Adams, New J. Phys. 8 , 163 (2006).
- [49] A. M. Kaufman, B. J. Lester, and C. A. Regal, Phys. Rev. X 2 , 041014 (2012).
- [51] T. L. Nguyen, J. M. Raimond, C. Sayrin, R. Corti˜as, T. Cantatn Moltrecht, F. Assemat, I. Dotsenko, S. Gleyzes, S. Haroche, G. Roux, T. Jolicoeur, and M. Brune, Phys. Rev. X 8 , 011032 (2018).
- [50] J. D. Thompson, T. G. Tiecke, A. S. Zibrov, V. Vuleti´ c, and M. D. Lukin, Phys. Rev. Lett. 110 , 133001 (2013).
- [52] M. J. Piotrowicz, M. Lichtman, K. Maller, G. Li, S. Zhang, L. Isenhower, and M. Saffman, Phys. Rev. A 88 , 013420 (2013).
- [53] Michel Brune, private communication.

## Supplemental Material: Three-dimensional trapping of individual Rydberg atoms in ponderomotive bottle beam traps

D. Barredo, V. Lienhard, P. Scholl, S. de L´s´leuc, T. Boulier, A. Browaeys and T. Lahaye e e Laboratoire Charles Fabry, Institut d'Optique Graduate School,

CNRS, Universit´ e Paris-Saclay, 91127 Palaiseau Cedex, France

## S.1. NUMERICAL SIMULATION OF THE ATOMIC MOTION IN THE BOB TRAP

To assess the quality of the BoB traps, we performed classical Monte-Carlo simulations of the atomic motion in the trapping potential. This allowed us to compute the recapture probability obtained in a release and recapture experiment (see Fig. 3a of the main text for the sequence) and compare it with experimental data. We simulate the full experimental sequence, including the two free-flight stages during Rydberg excitation and de-excitation, and the atomic motion in the BoB trap for a time τ . To determine if an atom is recaptured at the end of the sequence, we compare its kinetic energy to the potential energy of the ground state trap at the final position of the atom.

We account for the finite Rydberg lifetime at 300 K by introducing, for each Rydberg state, a single exponential probability distribution function in the dynamics. This approximation is justified because the two processes contributing to the Rydberg lifetime (decay to low-lying states and blackbody radiation-induced emission to other Rydberg states) produce the same outcome in this experiment: a loss of the atom in the final fluorescence image. Atoms that are not successfully excited by the STIRAP pulses, or Rydberg atoms that spontaneously decay to the ground state, are quickly repelled from the trapping region by the BoB potential and not recaptured at the end of the sequence. On the contrary, Rydberg atoms transferred to other Rydberg states by black-body radiation can still remain trapped in the BoB, but are not projected to the ground state by the de-excitation pulse. As a consequence, they are repelled by the repulsive potential of the standard optical tweezers before imaging, and finally lost.

For each simulation we use the same parameters as in the experiment (namely, a total power of 400 mW, the recorded three-dimensional trap intensity distribution shown in Fig. 2a of the main text, and the measured STIRAP excitation efficiency). The final result of the simulation is obtained by averaging over ∼ 10 4 realizations, starting from random atom positions (with rms values σ r = √ k B T/mω 2 r ∼ 230 nm; σ z = √ k B T/mω 2 z ∼ 1 15 . µ m) and velocities (rms value σ v = √ k B T/m ∼ 110 nm /µ s), according to a thermal atomic distribution in the standard optical tweezers at T = 130 µ K, with radial and longitudinal trapping frequencies ω / r (2 π ) = 75 kHz, and ω / z (2 π ) = 15 kHz, respectively.

Figure S1 shows the results of the simulation for two different Rydberg states, 75 S 1 / 2 , and 84 S 1 / 2 . The agreement between the experimental data and the simulation with no free parameters is very good over the full time span. At short times ( τ &lt; 50 µ s) we observe a fast decrease in the recapture probability (see inset) both in the experimental data and in the simulation. In the simulation, this initial drop has two contributions: (i) about 10 -15% of this decrease (depending on the dataset) corresponds to the loss of atoms that remain in the ground state after the STIRAP excitation pulses; (ii) the other 10% is due to energetic atoms having enough kinetic energy to overcome the potential barrier of the BoB trap.

Figure S1: Simulated recapture probability as a function of the trapping time τ (solid lines) for the Rydberg states 75 S 1 / 2 (a) and 84 S 1 / 2 (b). The experimental data (disks) are shown for comparison. Insets are zooms on the short-time behavior (dashed lines) showing, in linear scale, a fast atom loss at short times.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-caa918a9-2983-4253-82b2-4fd913f7c5dd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-92c4c480-9f16-45ec-8dd1-a44b6f9e3eab.jpg)

Figure S2: (a) Simulated recapture probability as a function of the trapping time for different Rydberg states assuming an infinite Rydberg lifetime. (b) Classical trajectories of Rydberg atoms evolving in the BoB potential. For the 60 S 1 / 2 state mechanical losses occur close to the saddle points ('escape cones'), whereas for 100 S 1 / 2 atoms leave the trap uniformly in the radial direction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-007c5607-e848-4bb3-ba1a-5e6da4e525b1.jpg)

## S.2. ANALYSIS OF THE MECHANICAL LOSSES IN THE BOB TRAP

In order to investigate the (small) initial atom loss and to distinguish 'mechanical losses' due to inefficient or unstable trapping from atom losses due to Rydberg state decay, we performed the same simulation as described in the previous section, but now assuming an infinite Rydberg lifetime. The results are shown in Fig. S2a.

For n ≤ 100, we observe again the fast initial decrease in the recapture probability in a time scale of ∼ 50 µ s, before reaching an almost stationary value at around τ = 1 ms. For times longer than τ ∼ 50 µ s the decay rate due to mechanical losses is comparatively much lower than the inverse of the Rydberg state lifetime. The fact that extra losses of atoms occur mainly at short times allows us to measure Rydberg trapping times close to the theoretical Rydberg state lifetimes. For n = 120, the potential does not have a local minimum (see Fig. 2 of the main text), and Rydberg atoms are quickly repelled from the trapping region.

The atom trajectories shown in Fig. S2b give us information about the escape regions in the trapping potential. For the 60 S 1 / 2 atoms escape mainly though the saddle points. For higher Rydberg states, the effect of the convolution is to smooth the saddle cones at the expense of decreasing the effective potential barrier, and atoms escape the trap in the radial direction in a more homogeneous manner.

## S.3. OSCILLATION FREQUENCY OF THE ATOMS IN THE BOB TRAP

We present in this section an experiment to measure the typical radial trapping frequency in the BoB trap. The experimental sequence is shown in Fig. S3a. We first transfer the Rydberg atom to the BoB trap. We then wait for 30 µ s. During this preliminary waiting time, the hottest atoms escape, as explained in the previous section. Then, we follow the same approach as in [S1], replacing the standard optical tweezers by our BoB trap. We switch off the trap during 4 µ s to excite the radial breathing mode. Then, atoms oscillate in the BoB trap for a varying duration τ . We

Figure S3: (a) Sketch of the experimental sequence to measure the oscillation frequency of the atoms in the BoB trap. (b) Experimental data (disks) and classical simulation (solid curve) of the 'release and recapture' experiment performed using the sequence in (a) for a trap depth of 0.6 mK.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00853v1-figure-50bcd81d-8d67-4093-904a-9a68bd40f5ec.jpg)

release the atoms for 10 µ s and finally de-excite and recapture them in the standard optical tweezers. For a harmonic trap, the recapture probability is expected to oscillate at 2 ω , with ω the trapping frequency.

Fig. S3b shows the measured recapture probability for a trap depth of 0.6 mK. From the observed oscillations we extract a trapping frequency, ω BoB / (2 π ) = 59 ± 3 kHz. The observed damping of the oscillations is most likely due to the finite atom temperature and the anharmonicity of the trapping potential. In the simulation, we follow the same procedure described above to numerically compute the recapture probability with no adjustable parameters. The result (solid line) is in reasonable agreement with the experimental data, as we find a trapping frequency ω BoB = 2 π × 59 kHz.

[S1] C. Tuchendler et al. , Phys. Rev. A 78 , 033425 (2008).
Daniel Barredo Vincent Lienhard Pascal Scholl Sylvain de Léséleuc Thomas Boulier Antoine Browaeys Thierry Lahaye
Three-dimensional trapping of individual Rydberg atoms in ponderomotive bottle beam traps
We demonstrate three-dimensional trapping of individual Rydberg atoms in
holographic optical bottle beam traps. Starting with cold, ground-state
$^{87}$Rb atoms held in standard optical tweezers, we excite them to
$nS_{1/2}$, $nP_{1/2}$, or $nD_{3/2}$ Rydberg states and transfer them to a
hollow trap at 850 nm. For principal quantum numbers $60 \leqslant n \leqslant
90$, the measured trapping time coincides with the Rydberg state lifetime in a
300~K environment. We show that these traps are compatible with quantum
information and simulation tasks by performing single qubit microwave Rabi
flopping, as well as by measuring the interaction-induced, coherent
spin-exchange dynamics between two trapped Rydberg atoms separated by 40
$\mu$m. These results will find applications in the realization of
high-fidelity quantum simulations and quantum logic operations with Rydberg
atoms.