# Single-site addressing of ultracold atoms beyond the diffraction limit via position dependent adiabatic passage

###### Abstract

We propose a single-site addressing implementation based on the sub-wavelength localization via adiabatic passage (SLAP) technique. We consider a sample of ultracold neutral atoms loaded into a two-dimensional optical lattice with one atom per site. Each atom is modeled by a three-level system in interaction with a pump and a Stokes laser pulse. Using a pump field with a node in its spatial profile, the atoms at all sites are transferred from one ground state of the system to the other via stimulated Raman adiabatic passage, except the one at the position of the node that remains in the initial ground state. This technique allows for the preparation, manipulation, and detection of atoms with a spatial resolution better than the diffraction limit, which either relaxes the requirements on the optical setup used or extends the achievable spatial resolution to lattice spacings smaller than accessible to date. In comparison to techniques based on coherent population trapping, SLAP gives a higher addressing resolution and has additional advantages such as robustness against parameter variations, coherence of the transfer process, and the absence of photon induced recoil. Additionally, the advantages of our proposal with respect to adiabatic spin-flip techniques are highlighted. Analytic expressions for the achievable addressing resolution and efficiency are derived and compared to numerical simulations for Rb atoms in state-of-the-art optical lattices.

## I Introduction

Ultracold neutral atoms in an optical lattice with single-atom and single-site resolution constitute an ideal physical system to investigate strongly correlated quantum phases Bloch'11 which, in turn, has interesting applications in quantum optics Ginsburg'07 , quantum simulation Bloch'08 and quantum information processing atoms_opt_lat ; Brennen'99 ; Raussendorf'01 , among others. The first approaches towards single-site addressing considered the use of lattices with relatively large site separations inc_latt_space . However, to have access to the regime of strongly correlated systems, typical lattice spacings well below m are needed since the tunneling rate has to be comparable to the on site interactions. In this case, the diffraction limit imposes strong restrictions on the addressability of individual lattice sites. To overcome this limitation different techniques have been investigated. For instance, spatially dependent electric and magnetic fields have been used to induce position dependent energy shifts on the atom ext_fields , allowing for site-selective addressability. Alternatively, a scanning electron microscopy system to remove atoms from individual sites with a focused electron beam SEM has been reported. However, in this case, atoms need to be reloaded into the emptied sites after each detection event. More recently, high resolution fluorescence imaging techniques, that make use of an optical system with high numerical aperture, have been implemented to perform in situ single-atom and single-site imaging for strongly correlated systems HRFI . In this context, a single-site addressing (SSA) scheme based on focused laser beams inducing position-dependent energy shifts of hyperfine states has been theoretically Zhang'06 and experimentally reported Weitenberg'11 . In the experiment, an intense addressing beam is tightly focused by means of a high resolution optical system. This beam produces spatial dependent light shifts bringing the addressed atom into resonance with a chirped microwave pulse and eventually inducing a spin-flip between two different hyperfine levels of the atom.

On the other hand, during past years, several proposals based on the interaction of spatially dependent fields, e.g., standing waves, with three-level atoms in a -type configuration have been considered, not only for single-site addressing in optical lattices but, more generally, for sub-wavelength resolution and localization Agarwal'06 ; swl_dark_state ; Yavuz'07 ; SLAP . In the first approaches Agarwal'06 ; swl_dark_state ; Yavuz'07 , a spatially modulated dark state is created by means of either electromagnetically induced transparency (EIT) or coherent population trapping (CPT) ArimondoBook'96 ; CPT-EIT , which allows for a tight localization of the atomic population in one of the ground states, around the position of the nodes of the spatially dependent field. More recently, it has been shown that the resolution achieved with those CPT or EIT based techniques can be surpassed using stimulated Raman adiabatic passage (STIRAP) Bergmann'98 processes, by means of the so-called sub-wavelength localization via adiabatic passage (SLAP) technique SLAP . The SLAP technique relies on a position-dependent STIRAP of atoms between the two ground states of a three-level atomic system, and has additional advantages compared with CPT or EIT techniques such as (i) robustness against parameter variations, (ii) coherence of the transfer process, that allows for its implementation also in Bose-Einstein condensates SLAP , and (iii) the absence of photon induced recoil.

In this work, we apply the SLAP technique SLAP to ultracold atoms in an optical lattice, where single-site addressing (SSA) requires one to overcome the diffraction limit. In order to address only a single site, we use here Stokes and pump pulses with Gaussian shaped spatial distributions, with the pump presenting a node centered at the lattice site that we want to address. Assuming that all the atoms in the optical lattice are initially in the same internal ground state and applying the standard STIRAP counterintuitive temporal sequence for the light pulses Bergmann'98 , we will demonstrate that it is possible to adiabatically transfer all the atoms, except the one at the node of the pump field, to an auxiliary ground state. We will show that this process is performed with higher efficiency and yields better spatial resolution than the CPT based techniques Agarwal'06 ; swl_dark_state ; Yavuz'07 . Also, we will demonstrate that our addressing technique requires shorter times than in the adiabatic spin-flip technique discussed in Ref. Weitenberg'11 , and that larger addressing resolutions can be achieved using similar focusing of the addressing fields. Moreover, our technique has the additional advantage that it can be applied between two degenerated ground-state levels.

The paper is organized as follows. In Sec. II we introduce the physical system under consideration. In Sec. III, we present a protocol to achieve SSA, and we derive analytical expressions for the spatial resolution and addressing efficiency of our technique. In addition, a comparison with CPT based techniques is provided. Next, in Sec. IV, we perform a numerical investigation of the proposed technique for a single-occupancy optical lattice loaded with Rb atoms by integrating the corresponding atomic density-matrix equations. Finally, in Sec. V, we summarize the results and present the conclusions.

## Ii Physical model

Figure 1(a) illustrates the physical system under consideration. It consists of a sample of ultracold neutral atoms loaded into a two-dimensional (2D) square optical lattice with spatial period , placed in the plane () and illuminated by a pump and a Stokes laser pulse with Rabi frequencies and , respectively, propagating in the direction with a selectable time delay. The spatial profile of the pump pulse has a node coinciding with a particular lattice site, our target site, at which the Stokes pulse is also centered. With the spatial profiles of the pulses having revolution symmetry around the propagation axis [dashed line in Fig. 1(a)], in the following we consider only the transverse spatial dimension , without loss of generality. In our model, we assume the system to be in the Mott insulator regime with only one atom per lattice site. Each atom is considered to have only three relevant energy levels in a -type configuration, defined by the interaction with the light pulses as shown in Fig. 1(b). Here, () is the spontaneous transition rate from the excited state to the ground state () and () is the detuning of the pump (Stokes) field. We assume that all atoms are initially in state .

Our approach to achieve single site addressing is based on the SLAP technique SLAP , where, depending on their position, the atoms are transferred between two internal ground states by means of the STIRAP technique Bergmann'98 . STIRAP consists in adiabatically following one of the energy eigenstates of the system, the so-called dark state, which under the two-photon resonance condition, i.e., , has the form

(1) |

where . Starting with all the population in , it is possible to coherently transfer the atomic population to state changing adiabatically from to by means of a convenient time sequence of the fields. This time sequence corresponds to apply first the Stokes pulse and, with a certain temporal overlap, the pump pulse. Since the process involves one of the eigenstates of the system, the population transfer is robust under fluctuations of the parameter values if these are adiabatically changed and the system does not evolve near degenerate energy eigenvalues.

In the SLAP technique, the pump field has a spatial structure with nodes yielding state-selective localization at those positions where the adiabatic passage process does not occur, i.e., those atoms placed at the nodes of the pump field remain in , while those interacting with both fields, pump and Stokes, are transferred to . For our purposes, we use the SLAP technique with a pump field having a single node at the position of the target site. Therefore, at the end of the SLAP process the population of all atoms illuminated is transferred from to except for the one at the node of the pump field. The spatial and temporal profiles for pump and Stokes Rabi frequencies are given by

(2) | |||||

(3) |

where and are the peak Rabi frequencies, and are the centers of the temporal Gaussian profiles, and are the spatial widths of the node in the pump and of the Stokes field, respectively, and is the temporal width.

There exist several methods to create the required pump intensity profile with a central node: e.g. (i) re-imaging of a Gaussian beam with a dark central spot created by a circular absorption mask, using (ii) a Laguerre-Gaussian laser beam Heckenberg'92 or (iii) a “bottle beam” created by the interferometric overlap of two Gaussian beams with differing waists Isenhower'09 , or (iv) a flexible intensity pattern generated by spatial light modulators and subsequent imaging Bergamini'04 ; TUD_SLM_work .

## Iii Single-Site Addressing

In our model, we assume that the spatial wavefunctions of the individual atoms placed at the different sites, centered at (being the site index), correspond to the ground state of the trapping potential, which in first approximation can be considered harmonic. Therefore, the full atomic distribution in the lattice is given, initially, by

(4) |

where is the width of the initial atomic distribution at an individual site, is the mass of the trapped atom and is the harmonic trapping frequency. We assume that the addressed site is and their nearest neighbors are at a distance where is the wavelength of the fields that create the optical lattice.

In order to characterize our single-site addressing technique we consider that, once the SLAP technique has been applied, the final atomic population distribution in , , is given by

(5) |

where is the probability distribution that an atom remains in state after the SLAP process. Using the SLAP technique, the addressing resolution that one can obtain is related to the global adiabaticity condition Bergmann'98 at each spatial position ,

where and is a dimensionless constant that, for optimal Gaussian temporal profiles and overlapping times, takes values around 10 Bergmann'98 ; LaA . In Eq. (III), the equality gives a spatial threshold above which the adiabaticity condition is fulfilled. Assuming that the full width at half maximum (FWHM) of is and expanding Eq. (III) up to first order in one obtains

where and . Equation (III) gives the width of the addressing region, and it tends to zero as increases. Moreover, since must be real valued, we find that the inequality

(8) |

must be fulfilled. In this paper we will consider that is fixed, so can be varied through , and .

Two conditions should be satisfied for our SSA technique to work. First, the population of the atom in the addressed site must remain in state after the action of the fields and, second, the rest of the atoms of the lattice have to be transferred to level . Therefore, taking into account the overlap between and in Eq. (5), it is clear that the FWHM of the probability distribution should satisfy

(9) |

where , and is the position of the nearest-neighboring site. Using Eq. (III), it is easy to see that these conditions fix the range for to obtain SSA using the SLAP technique:

(10) |

where

(11) |

with and . Note that the upper limit for Eq. (10) is more restrictive than Eq. (8).

In order to have a quantitative description of the SSA performance, let us introduce the SSA efficiency as

(12) |

where corresponds to the probability of finding the atom at the addressed site in state , while corresponds to the probability that the atom in the neighbor site has been transferred to a different internal state. We define

(13) |

with and , whereas and have been defined in Eqs. (5) and (4), respectively. Using Eqs. (4) and (5), the explicit forms for Eqs. (13) are

(14) | |||||

(15) |

From these expressions, it can be seen that, for , the limits given by Eq. (9), i.e., and , correspond to SSA efficiencies of 0.70 and 0.94, respectively.

An alternative technique to perform atomic localization based on spatial dependent dark states is the coherent population trapping (CPT) technique Agarwal'06 . In the CPT ArimondoBook'96 technique the dark state is populated after several cycles of coherent excitation followed by spontaneous emission from to the ground states. Note that, while CPT relies on spontaneous emission, the SLAP technique is fully coherent. Moreover, the latter provides higher resolution, as shown in Ref. SLAP , and does not suffer from recoil since the localized atoms have not interacted with light. In what follows, we compare the range of parameters necessary to perform SSA considering both SLAP and CPT techniques. Note that we focus our comparative analysis in the CPT technique, although similar results are obtained considering the method proposed in Ref. Yavuz'07 , where the spatial dependent dark state, , is created via the STIRAP technique by switching off the fields before completing the transfer process.

In order to compare both techniques, we define the final population distribution in using CPT as in an analogous way as it has been done for the SLAP technique in Eq. (5). Then the FWHM of the corresponding probability function is obtained by imposing that and in Eqs. (2) and (3):

(16) |

Fixing the desired and using Eqs. (III) and (16) for SLAP and CPT, respectively, the constraints for the relevant parameters for each technique can be obtained. For simplicity, we consider that the Stokes pulse parameters and are fixed, and only the node width and pump peak Rabi frequency can be varied. Note that, for the SLAP case, we have to fix also and .

Taking , half of the site separation, the required values for and are plotted in Fig. 2(a). To simultaneously illuminate a large number of sites, we use a large Stokes beam waist of , where is the wavelength of both pump and Stokes fields Note_Lambda . The solid (dashed) line corresponds to the SLAP (CPT) case with the parameter values and LaA . As decreases, in both SLAP and CPT cases, lower values of are needed to reach the fixed resolution , since the narrower the node of the pump field, the narrower the probability distribution of atoms remaining in . In addition, for any given width of the node of the pump field , the required values of are lower in the SLAP case than in the CPT case.

It is important to realize that Eqs. (III) and (16) show the possibility to obtain, for certain parameter values, widths of the probability distribution, or , smaller than . In particular, using the SLAP technique this can be achieved with moderate values. This is shown in Fig. 2(b), where (solid line) and (dashed line) are represented as a function of for and the rest of the parameters as in Fig. 2(a). The FWHM of the atomic distribution, , is depicted with a horizontal dotted line to indicate the values where . As we stated in the discussion of Eqs. (14) and (15), this limit corresponds to a SSA efficiency of . This regime of parameters is interesting because it shows that the SLAP technique could be used for applications in site-selective imaging with a resolution down to the width of the atomic distribution at each site.

## Iv Numerical Simulations

In this section, by numerically integrating the corresponding density-matrix equations, we study the implementation of the SLAP-based SSA technique for -type three-level Rb atoms in a single-occupancy optical lattice. Numerical calculations using the CPT technique are also presented for comparison. The wavelength of the lasers that create the optical lattice, red detuned with respect to the line of Rb, is nm. The potential depth of the optical lattice is chosen as , where is the recoil energy. This corresponds to a harmonic trapping frequency of kHz Bloch'08 . Therefore, the FWHM of the atom distribution at each site due to the confining potential is nm. Pump and Stokes fields with nm are coupled to and , respectively, where , and are hyperfine energy levels of the line of Rb. The excited state has a spontaneous transition rate MHz ( MHz) to state (), and we assume no spin decoherence during the interaction time. We consider equal temporal pulse widths of s with a temporal delay , in such a way that the total SSA process time is . The Stokes pulse has a maximum Rabi frequency MHz, while the maximum Rabi frequency of the pump is varied through the parameter , since . Concerning the spatial profiles of the fields, we assume a wide Stokes profile, , and a narrow node for the pump, .

As it has been discussed in the previous section, to properly perform SSA, the population of all the atoms in the lattice, except the one in the addressed site, must be transferred from to with high probability. In what follows, those requirements for the realization of the SSA are studied by numerically evaluating the SSA efficiency.

The signatures of SSA are shown in Fig. 3, where the numerically evaluated efficiency and probabilities defined in Eqs. (12) and (13) are plotted as a function of for both the SLAP (circles) and CPT (crosses) techniques. Figures 3(a) and 3(b) show the probabilities and , respectively. For large values of the probability of finding the atom at the addressed site in is higher with the CPT than with the SLAP technique [see Fig. 3(a)], while for small values of the probability of removing the atom from the neighboring site [see Fig. 3(b)] is higher in the SLAP case. This is explained by the fact that for a given value of , as shown in Fig. 2(b). The SSA efficiency results obtained by the numerical evaluation of Eqs. (12) and (13) are shown in Fig. 3(c), together with the corresponding analytical curves (solid and dashed lines) obtained using Eqs. (14) and (15), added for comparison. A good agreement is found between numerical and analytical results. From Fig. 3(c) it is clear that the SLAP technique is more efficient for lower values of the intensity ratio of the addressing fields () than the CPT technique. Certainly, this is advantageous for the experimental implementation with limited laser power available.

An example of the final population distribution after performing SSA with the SLAP (circles) and the CPT (crosses) techniques is plotted in Fig. 4(b) for the particular case of and with the rest of the parameters as in Fig. 3. The initial population distribution in , , at the addressed site (), and the two next neighbors nm is shown as a solid line. Note that, in the SLAP case (circles), the population of state around remains almost the same after the addressing process, while in the first neighbor sites it is practically zero. On the other hand, for the CPT case (crosses), the population in the addressed site remains also nearly unchanged, but it exhibits a significant amount of population in the neighbor sites. This is in full agreement with the discussion following Figs. 3(a) and 3(b). For this example, the total efficiencies found are and according to the corresponding values shown in Fig. 3(c). In addition, as shown in Fig. 4(a), the width of the pump node required to perform the SSA method is much larger than the addressed region . In particular, for the case shown in Fig. 4, and , thus obtaining addressing resolution beyond the diffraction limit.

Finally, we carry out a comparison between our proposal and the experiment reported in Ref. Weitenberg'11 , where a focused laser beam induces position-dependent light shifts, allowing one to perform a spin flip by means of a resonant microwave pulse at the addressed site. Since the microwave field involved in the experiment has a Rabi frequency of kHz, the total spin-flip time is in the order of ms. In contrast, as our proposal makes use of only optical fields, the addressing time is three orders of magnitude below (s). Specifically, to achieve similar values of the addressing resolution in both techniques, nm, we have obtained an addressing time of 40s. This decrease of the total addressing time needed implies a reduction of the effects caused by spontaneous scattering of photons, which are a limitation for the light shifts-based proposals Zhang'06 ; Weitenberg'11 . These effects could be strongly reduced in our case. Also, we have compared the resolution obtained in both techniques. In Ref. Weitenberg'11 , they use an addressing beam with an intensity FWHM of approximately 600 nm, and obtain a spin-flip probability distribution with FWHM=330 nm. In our technique, using and a width of the pump node nm, which corresponds to the width of their addressing beam, a very similar FWHM of the addressing probability distribution is obtained: nm. Note that this value can be reduced by increasing the ratio between the pump and the Stokes intensities, e.g., implies nm and for we obtain nm.

## V Conclusions

In this work we have discussed a proposal to perform single-site addressing (SSA) in an optical lattice by using the SLAP technique. With respect to other dark state based techniques such as CPT, this method is fully coherent, robust against variations of the parameter values, and we have found that it yields higher efficiencies for smaller values of the intensity ratio between the pump and the Stokes fields. Moreover, the addressed atom does not interact with the fields, minimizing all possible decoherent processes. On the other hand, with respect to the recent experiment on SSA using adiabatic spin flips Weitenberg'11 , the present proposal allows one to use two degenerate ground levels, takes shorter times to perform the addressing process, and provides similar or even larger addressing resolutions for similar focusing of the addressing fields.

The proposed method provides an achievable addressing resolution that can be pushed well below the diffraction limit of the addressing light field and of the optical setup used for addressing or detection of atoms at closely spaced lattice sites. This relaxes the requirements on the optical setup or extends the achievable spatial resolution to lattice spacings smaller than accessible to date. Through analytical considerations, we have derived the range of parameters for which SSA is properly achieved. Moreover, we have obtained expressions to estimate the resolution and the efficiency of the SLAP-based addressing method, and we have compared them with the analogous expressions obtained using the CPT technique. Next, by integrating the density-matrix equations with realistic parameter values for state-of-the-art optical lattices loaded with Rb atoms, we have checked the validity of the analytical approach.

###### Acknowledgements.

We thank Albert Benseny for fruitful discussions and comments. We acknowledge financial support from the Spanish Ministry of Science and Innovation under Contracts No. FIS2008-02425, No. FIS2011-23719, and No. CSD2006-00019 (Consolider project Quantum Optical Information Technologies), from the Catalan Government under Contract No. SGR2009-00347, from the German Research Foundation DFG (Contract No. BI 647/3-1), and from the German Academic Exchange Service DAAD (Contract No. 0804149).## References

- (1) M. Endres, M. Cheneau, T. Fukuhara, C. Weitenberg, P. Schauß, C. Gross, L. Mazza, M. C. Banuls, L. Pollet, I. Bloch and S. Kuh, Science 334, 200 (2011).
- (2) N. S. Ginsberg, S. R. Garner, and L. V. Hau, Nature 445, 623 (2007).
- (3) I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 (2008).
- (4) G. K. Brennen, C. M. Caves, P. S. Jessen, and I. H. Deutsch, Phys. Rev. Lett. 82, 1060 (1999).
- (5) R. Raussendorf and H. J. Briegel, Phys. Rev. Lett. 86, 5188 (2001).
- (6) D. Jaksch and P. Zoller, Ann. Phys. (N.Y.) 315, 52 (2005).
- (7) R. Scheunemann, F. S. Cataliotti, T. W. Hänsch, and M. Weitz, Phys. Rev. A 62, 051801(R) (2000); R. Dumke, M. Volk, T. Müther, F. B. J. Buchkremer, G. Birkl, and W. Ertmer, Phys. Rev. Lett. 89, 097903 (2002); S. Peil, J. V. Porto, B. Laburthe Tolra, J. M. Obrecht, B. E. King, M. Subbotin, S. L. Rolston, and W. D. Phillips, Phys. Rev. A 67, 051603(R) (2003); Y. Miroshnychenko, W. Alt, I. Dotsenko, L. Förster, M. Khudaverdyan, D. Meschede, D. Schrader, and A. Rauschenbeutel, Nature 442, 151 (2006); K. Nelson, X. Li, and D. Weiss, Nature Phys. 3, 556 (2007); G. Birkl and J. Fortágh, Laser & Photon. Rev. 1, 12 (2007).
- (8) D. Schrader, I. Dotsenko, M. Khudaverdyan, Y. Miroshnychenko, A. Rauschenbeutel, and D. Meschede, Phys. Rev. Lett. 93, 150501 (2004).
- (9) P. Würtz, T. Langen, T. Gericke, A. Koglbauer, and H. Ott, Phys. Rev. Lett. 103, 080404 (2009); T. Gericke, P. Würtz, D. Reitz, T. Langen, and H. Ott, Nature Phys. 4, 949 (2008).
- (10) W. S. Bakr, J. I. Gillen, A. Peng, S. Fölling, and M. Greiner, Nature 462, 74 (2009); W. S. Bakr, A. Peng, M. E. Tai, R. Ma, J. Simon, J. I. Gillen, S. Fölling, L. Pollet, and M. Greiner, Science 329, 547 (2010); N. Gemelke, X. Zhang, C.-L. Hung, and C. Chin, Nature 460, 996 (2009). M. Karski, L. Förster, J. M. Choi, W. Alt, A. Widera, and D. Meschede, Phys. Rev. Lett. 102, 053001 (2009); J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, I. Bloch, and S. Kuhr, Nature 467, 68 (2010).
- (11) C. Zhang, S. L. Rolston, and S. Das Sarma, Phys. Rev. A 74, 042316 (2006).
- (12) C. Weitenberg, M. Endres, J. F. Sherson, M. Cheneau, P. Schauß, T. Fukuhara, I. Bloch, and S. Kuhr, Nature 471 319 (2011).
- (13) G. S. Agarwal, and K. T. Kapale, J. Phys. B: At. Mol. Opt. Phys. 39, 3437 (2006).
- (14) J. Cho, Phys. Rev. Lett. 99, 020502 (2007); A. V. Gorshkov, L. Jiang, M. Greiner, P. Zoller, and M. D. Lukin, Phys. Rev. Lett. 100, 093005 (2008); M. Holland, S. Marksteiner, P. Marte, and P. Zoller, Phys. Rev. Lett. 76, 3683 (1996); E. Paspalakis and P. L. Knight, Phys. Rev. A 63, 065802 (2001); M. Sahrai, H. Tajalli, K. T. Kapale, and M. S. Zubairy, Phys. Rev. A 72, 013820 (2005); M. Kiffner, J. Evers, and M. S. Zubairy, Phys. Rev. Lett. 100, 073602 (2008).
- (15) D. D. Yavuz and N. A. Proite, Phys. Rev. A 76, 041802(R) (2007).
- (16) J. Mompart, V. Ahufinger, and G. Birkl, Phys. Rev. A 79, 053638 (2009).
- (17) E. Arimondo, in Progress in Optics, edited by E. Wolf (Elseveir Science, Amsterdam, 1996), Vol. XXXV, p. 257.
- (18) G. Alzetta, L. Moi and G. Orriols, Nuovo Cimento 52B, 209 (1979); S. E. Harris, Phys. Today 50(7) 36, (1997); J. P. Marangos, J. Mod. Opt. 45, 471 (1998).
- (19) K. Bergmann, H. Theuer, and B. W. Shore, Rev. Mod. Phys. 70, 1003, (1998).
- (20) N.R. Heckenberg, R. McDuff, C.P. Smith, H. Rubinstein-Dunlop, and M.J. Wegener, Opt. Quantum Electron. 24, S951, 1992.
- (21) L. Isenhower, W. Williams, A. Dally, and M. Saffman, Optics Letters 34, 1159 (2009).
- (22) S. Bergamini, B. Darquie, M. Jones, L. Jacubowiez, A. Browaeys, and P. Grangier, J. Opt. Soc. B 21, 1889, 2004.
- (23) J. Kruse, C. Gierl, M. Schlosser, and G. Birkl, Phys. Rev. A 81, 060308(R), 2010.
- (24) Although usually the value is accepted, we have found that gives a more accurate description for the values of and considered along this paper.
- (25) We choose this ratio in order to be consistent with the specific numerical values used in Sec. IV.