Controllable motion of optical vortex arrays using electromagnetically induced transparency

The ability to control optical vortices is a desired feature in numerous applications. We demonstrate control of the collective motion of a vortex array using an electromagnetically induced transparency media. Scanning the frequency detuning between the pump and probe fields changes the susceptibility of the media, producing a unique effective diffraction of the vortex array for each detuning. We measure several experimental configurations and compare them to numerical simulations. © 2012 Optical Society of America OCIS codes: (020.1670) Coherent optical effects; (050.4865) Optical vortices; (070.4690) Morphological transformations. References 1. K. Boller, A. Imamoglu, and S. Harris, “Observation of electromagnetically induced transparency,” Physical Review Letters 66, 2593–2596 (1991). 2. M. Fleischhauer, A. Imamoglu, and J. Marangos, “Electromagnetically induced transparency: Optics in coherent media,” Reviews of Modern Physics 77, 633–673 (2005). 3. P. K. Vudyasetu, D. Starling, and J. C. Howell, “All Optical Waveguiding in a Coherent Atomic Rubidium Vapor,” Physical Review Letters 102, 123602 (2009). 4. O. Firstenberg, P. London, M. Shuker, A. Ron, and N. Davidson, “Elimination, reversal and directional bias of optical diffraction,” Nature Physics 5, 665–668 (2009). 5. O. Firstenberg, M. Shuker, N. Davidson, and A. Ron, “Elimination of the Diffraction of Arbitrary Images Imprinted on Slow Light,” Physical Review Letters 102, 043601 (2009). 6. R. Pugatch, M. Shuker, O. Firstenberg, A. Ron, and N. Davidson, “Topological Stability of Stored Optical Vortices,” Physical Review Letters 98, 203601 (2007). 7. H. He, M. Friese, N. Heckenberg, and H. Rubinsztein-Dunlop, “Direct Observation of Transfer of Angular Momentum to Absorptive Particles from a Laser Beam with a Phase Singularity,” Physical Review Letters 75, 826– 829 (1995). 8. A. Mair, A. Vaziri, G. Weihs, and A. Zeilinger, “Entanglement of the orbital angular momentum states of photons.” Nature 412, 313–316 (2001). 9. S. Franke-Arnold, L. Allen, and M. Padgett, “Advances in optical angular momentum,” Laser & Photonics Review 2, 299–313 (2008). 10. Z. Chen, J. Pu, and D. Zhao, “Tight focusing properties of linearly polarized Gaussian beam with a pair of vortices,” Physics Letters A 375, 2958–2963 (2011). 11. D. Neshev, “Motion control of ensembles of ordered optical vortices generated on finite extent background,” Optics Communications 151, 413–421 (1998). 12. G. Indebetouw, “Optical Vortices and Their Propagation,” Journal of Modern Optics 40, 73–87 (1993). 13. D. Rozas, Z. S. Sacks, and G. a. Swartzlander, Jr., “Experimental Observation of Fluidlike Motion of Optical Vortices,” Physical Review Letters 79, 3399–3402 (1997). 14. D. Rozas, C. T. Law, and G. A. Swartzlander, Jr., “Propagation dynamics of optical vortices,” Journal of the Optical Society of America B 14, 3054–3065 (1997). 15. X. Gan, J. Zhao, S. Liu, and L. Fang, “Generation and motion control of optical multi-vortex,” Chinese Optics Letters 7, 1142–1145 (2009). 16. O. Firstenberg, M. Shuker, R. Pugatch, D. Fredkin, N. Davidson, and A. Ron, “Theory of thermal motion in electromagnetically induced transparency: Effects of diffusion, Doppler broadening, and Dicke and Ramsey narrowing,” Physical Review A 77, 043830 (2008). 17. M. Chen and F. S. Roux, “Accelerating the annihilation of an optical vortex dipole in a Gaussian beam.” Journal of the Optical Society of America. A, Optics, image science, and vision 25, 1279–1286 (2008). 18. O. Firstenberg, P. London, D. Yankelev, R. Pugatch, M. Shuker, and N. Davidson, “Self-Similar Modes of Coherent Diffusion,” Physical Review Letters 105, 183602 (2010). 19. A. Y. Bekshaev, M. S. Soskin, and M. V. Vasnetsov, “Optical vortex symmetry breakdown and decomposition of the orbital angular momentum of light beams,” Journal of the Optical Society of America A 20, 1635–1643 (2003). 20. P. K. Vudyasetu, R. M. Camacho, and J. C. Howell, “Storage and Retrieval of Multimode Transverse Images in Hot Atomic Rubidium Vapor,” Physical Review Letters 100, 123903 (2008). 21. M. Shuker, O. Firstenberg, R. Pugatch, A. Ron, and N. Davidson, “Storing Images in Warm Atomic Vapor,” Physical Review Letters 100, 223601 (2008). 22. E. E. Mikhailov, I. Novikova, M. D. Havey, and F. A. Narducci, “Magnetic field imaging with atomic Rb vapor,” Optics Letters 34, 3529–3531 (2009).


Introduction
The optical susceptibility of atomic media due to a weak resonant probe field can be changed dramatically in the presence of a strong pump field. The imaginary part of the susceptibility vanishes, thus leading to a transparency in the otherwise opaque media, while the real part of the susceptibility becomes very steep, resulting a great reduction in the group velocity as compared with the speed of light. This effect is termed electromagnetically induced transparency (EIT) [1,2]. An extremely rich variety of spatially dependent optical phenomena that is derived from the complex susceptibility of EIT media was demonstrated recently, such as waveguiding [3], frequency dependent diffraction, diffractionless propagation, Negative diffraction and induced deflection [4,5]. Transmission of phase dependent patterns were also achieved and highlighted as an effective way to reduce diffusion effects in storage processes [6,7].
Here we report on a different phenomenon related to the complex susceptibility, namely, a controllable collective motion of optical vortex arrays. An optical vortex is an electromagnetic mode with zero intensity in the center and a cylindrical phase e imφ , where φ is defined as an angle variable around the vortex core and m is an integer quantizing the orbital angular momentum (OAM) of the vortex. Optical vortices have a wide range of applications such as optical trapping and tweezing [8], quantum entanglement [9,10] and contrast enhancement in light microscopy [10]. Single vortex beams were used as a probe beam in EIT experiments, and demonstrated higher stability of the vortex against diffusion effects, as opposed to a bagel shaped flat phase beam [11]. The transverse collective motion of optical vortex array patterns as they propagate in free-space was studied extensively [12][13][14] and was found to depend upon the total OAM of the array (denoted by L), the average distance between the vortices, the vortices' core size and the background profile they are embedded in. In general, vortex arrays with L =0 were shown to rotate as the beam propagates due to the collective OAM of the beam, while vortex arrays with L=0 do not rotate at all but can contract and annihilate. The movement of each vortex is determined by the phase and intensity gradients of the beam in the vortex core. Thus the movement of each vortex is driven by the phase created by the other vortices in the array, leading to a collective motion [15][16][17]. We control this collective motion of the vortex array in EIT media by changing the frequency of our pump field, varying it from resonance with the two-photon transition of the EIT (two-photon detuning).
The concept behind the manipulation of vortex arrays is the controllable diffraction of the probe image due to the tunable susceptibility. For an undepleted plane wave pump beam, the susceptibility of a spatially varying probe in EIT can be written as [18]: where k ⊥ is the transverse wavevector of the probe profile, g is the dipole interaction strength, n 0 is the density of the vapor, K(ω) is the dimensionless susceptibility in the absence of the pump field andL is the following spatio-temporal spectral function : Here Ω c is the Rabi frequency of the coupling field, Δ is the two photon detuning, Γ hom is the decoherence rate and D = v 2 th γ is the Dicke-Doppler diffusion coefficient where γ is the relaxation rate due to collisions with the buffer gas and v th = k B T m is the thermal velocity of the atoms. The dynamics of the probe pattern due to this complex susceptibility is given by where λ is the optical wavelength. This dynamics can be described as a diffusion equation with a complex diffusion constant, meaning that it undergoes a frequency dependent optical diffraction together with a diffusion due to the atomic motion [18]. The diffraction component creates a behavior similar to the vortex array motion in free-space, but with a frequency controlled effective diffraction length. Hence, as the two photon detuning is scanned we expect to see the effects that were described before in free-space, such as a rotation of vortex arrays with L =0 and a contraction in the case of arrays with L=0. The diffusion component leads to a motion also in the direction of the intensity gradient and not only perpendicular to it. This outcome is unique to the propagation in diffusive media and has no counterpart in free-space propagation.
A major advantage of this method for applications such as optical trapping and optical wrenches is the ability to move vortices without a physical change of the input probe field, which is usually a slow process, but a frequency control which is potentially much faster and more robust.

Experimental setup
Our experimental setup is described in Fig. 1. As the EIT medium we use the D1 transition of 87 Rb, where the Λ configuration is comprised of the hyperfine levels. We use a distributed feedback laser as our source of 795 nm beam. This beam is separated into pump and probe beams using a beam splitter. Both beams pass through an acousto-optic modulator, the probe for amplitude control and the pump for amplitude and frequency control using a double pass configuration. A hyperfine splitting of 6.8 GHz is created using an electro-optic modulator on the probe beam followed by an etalon used to filter out the unwanted 0 and -1 orders. The EOM is used also to compensate for the frequency mismatch of the probe due to the double pass of the pump. The spatial phase dependence of the beam is created using a liquid crystal spatial light modulator (BNS, 512 × 512 pixels). Finally, both the pump and the probe beams are combined using a Glan-Taylor polarizer and pass through a 7.5 cm quartz cell heated to ∼ 40 o C. The cell contains isotopically pure 87 Rb and 20 Torr Neon as buffer gas. After the cell, the pump and probe beams are separated using another Glan-Taylor polarizer and the probe beam is measured using a CCD camera (PCO Pixelfly qe) with short exposure times of 5 − 10μs in order to avoid distortions due to phase fluctuations. For the results described here we use a probe beam with a power of 5μW and a pump beam with a power of 150μW and a waist of 2.5 mm. Figure 2 shows the phase patterns imprinted on the spatial light modulator for a 2 × 2 vortex array in two configurations. The vortex array in Fig. 2(a) comprises two vortices having OAM of +1 and two vortices having OAM of -1, so the total OAM of the pattern is L=0, while in Fig.  2(b) all the vorices have OAM of +1, so the total OAM of the pattern is L=4. We chose these  beam with a 2 × 2 vortex array embedded in it. Figure 3(a) is produced using the phase pattern depicted in Fig. 2(a), while Fig. 4(a) is produced using the phase pattern presented in Fig. 2(b). (1) of both the situations discussed with good agreement to the experimental results. The simulation uses the measured temperature and pump intensity, while γ is calculated ab initio using a hard sphere model. g and K(ω) are the only fitting parameters (minor adjustments from calculated/measured values due to experimental uncertainties) used in order to achieve the measured EIT absorption spectrum. Once the EIT spectrum is matched, the spatial propagation is calculated for all patterns without any further fitting.

Figures 3(d) and 4(d)
show the dynamics of each vortex center as the two photon detuning is scanned. For L=4 a global rotation accompanied by vortex repulsion is observed (Fig. 4(c)). The rotation is a signature of the non-zero angular momentum of the array while the repulsion has to do with the movement of the vortices due to the gaussian envelope of the beam. This repulsion is higher when the rotation is higher thus leading to a more distorted images for higher total L patterns. The circular trajectory of each vortex as described in Fig. 4(d) is very different from free-space propagation where a linear trajectory was observed [14,15]. The circular trajectory is due to the combined effects of diffusion and diffraction. L=0 has a different signature with a global contraction (Fig. 3(c)) and no global rotation. This effect happens since each two adjacent vortices have inverse OAM thus leading to rotation in an opposite direction. Moreover, as the pattern with total OAM of L=0 contracts, the vortices start to join together leading to one large dark spot in the middle of the beam. Simulation shows that for strong enough interaction the global contraction for patterns may result in even a total annihilation of the vortices, in agreement with published free-space results of vortex dipoles [19]. As a quantitative measurement of the global rotation that arrays with L =0 experience, we measure the angle of the connecting line between two vortices with respect to the x-axis. A total rotation span of more than 0.3 radians is visible (Fig. 5(a)).
The rotation has a dispersive signature, similar to the real part of the susceptibility, where it is possible to see an extra rotation in one side of the resonance and a contra free-space rotation on the other side. The EIT media in this case enables a manipulation of the effective optical path, thus allows us to change the optical pattern as if it was taken at different positions along the beam traveling axis. Actually, free space propagation like this would have taken about twice the Rayleigh length which is 1 meter in our case. In this view it is possible to think of the medium as a very high numerical aperture lens [12]. Moreover, contra free-space rotation is analogous to a backward movement of the beam from the origin. This effect is a signature of an index of refraction which is smaller than one. Figure 5(a) also demonstrates that positive or negative total angular momentum creates an opposite rotation (blue squares and red circles curves), while zero OAM produces no rotation at all (green diamonds). Due to the imaginary part of the complex susceptibility an actual diffusion of the beam occurs. This effect was observed and measured before by Firstenberg et al. [20]. Here we show an easy way to measure it by looking at the average distance between vortices with the same optical angular momentum as depicted in Fig. 5(b). This increase in the distance between the vortices agrees well with the EIT line presented also in Fig. 5(b). Larger arrays of vortices such as 3 × 3 and 4 × 4 were also created with similar effects as shown for the 2 × 2 array case. Due to the fact that the vortices cover a large part of the Gaussian beam, the EIT peak diminishes, thus leading to slightly less pronounced effects. As the dynamics of each vortex is affected by the local field in the point of the vortex, excluding the vortex itself, using a Gaussian field offers a way to manipulate even a single off-centered vortex alone [21]. Moreover, here we exploit the linear regime of the EIT media, but higher intensity will probe the nonlinear Kerr medium with extended rotation spans, not possible in the linear regime [16].
Vortex arrays can be stored in EIT media in a similar manner to previous experiments of storage of patterns using EIT [6,11,22]. We measure stored patterns of vortex arrays for different storage times. The stored patterns are similar to the patterns produced in the slow light regime apart from atomic diffusion effects that cause the dark center of the vortices to fill up as the storage time increases. This trend is in agreement with single vortex storage results presented by Pugatch et al. [11]. A typical measurement of such storage dynamics for two vortices with total L=0 can be seen in Fig. 6.
One advantage stored patterns have over slow light patterns is an enhanced purity of the measured vortex array state. The reason for this enhancement is the fact that only the EIT controlled light is stored. In contrast, in the slow light regime, various spectral contaminations can occur in the beam profile due to minor optical pumping effects and incoherent light.

Conclusion
We spectrally control the collective motion of optical vortex arrays by the use of EIT media. This control is done by changing the two photon detuning between the pump and probe beams. We describe how to use this method for various patterns of vortex arrays and explain the phenomenological differences between arrays with L=0 and L =0, namely rotation and expansion of the array. This control may be used for fast and accurate displacement of optical vortices in applications such as optical wrenches and the transfer of optical angular momentum to particles [8]. The ability to detect small changes in the displacement of the vortices may also be used in magnetic sensing [23] and quantum measurements of multiphoton states with nontrivial spatial profiles.