Generation of a lattice of spin-orbit beams via coherent averaging

Generation of a lattice of spin-orbit beams via coherent averaging

D. Sarenac Department of Physics, University of Waterloo, Waterloo, ON, Canada, N2L3G1 Institute for Quantum Computing, University of Waterloo, Waterloo, ON, Canada, N2L3G1    D. G. Cory Institute for Quantum Computing, University of Waterloo, Waterloo, ON, Canada, N2L3G1 Department of Chemistry, University of Waterloo, Waterloo, ON, Canada, N2L3G1 Perimeter Institute for Theoretical Physics, Waterloo, ON, Canada, N2L2Y5 Canadian Institute for Advanced Research, Toronto, Ontario, Canada, M5G 1Z8    J. Nsofini Department of Physics, University of Waterloo, Waterloo, ON, Canada, N2L3G1 Institute for Quantum Computing, University of Waterloo, Waterloo, ON, Canada, N2L3G1    I. Hincks Institute for Quantum Computing, University of Waterloo, Waterloo, ON, Canada, N2L3G1 Department of Applied Math, University of Waterloo, Waterloo, ON, Canada, N2L3G1    P. Miguel Institute for Quantum Computing, University of Waterloo, Waterloo, ON, Canada, N2L3G1 Department of Chemistry, University of Waterloo, Waterloo, ON, Canada, N2L3G1    M. Arif National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA    Charles W. Clark National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA Joint Quantum Institute, National Institute of Standards and Technology and University of Maryland, College Park, Maryland 20742, USA    M. G. Huber National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA    D. A. Pushin Department of Physics, University of Waterloo, Waterloo, ON, Canada, N2L3G1 Institute for Quantum Computing, University of Waterloo, Waterloo, ON, Canada, N2L3G1
July 7, 2019

We describe a highly robust method, applicable to both electromagnetic and matter-wave beams, that can produce a beam consisting of a lattice of orbital angular momentum (OAM) states coupled to a two-level system. We also define efficient protocols for controlling and manipulating the lattice characteristics. These protocols are applied in an experimental realization of a lattice of optical spin-orbit beams. The novel passive devices we demonstrate here are also a natural alternative to existing methods for producing single-axis OAM and spin-orbit beams. Our techniques provide new tools for investigations of chiral and topological materials with light and particle beams.

Valid PACS appear here
preprint: APS/123-QED

Since their experimental demonstrations a quarter-century ago Bazhenov et al. (1990); Allen et al. (1992), there has been great progress in generation, detection, and applications of “structured waves” of light and quantum particles, where the wavefront is patterned to attain nontrivial propagation characteristics such as orbital angular momentum (OAM), non-diffraction, and self-healing Harris et al. (2015); Garces-Chavez et al. (2002); Siviloglou et al. (2007); Barnett et al. (2017); Molina-Terriza et al. (2007). Beams of light Allen et al. (1992), neutrons Clark et al. (2015); Sarenac et al. (2016) and electrons Uchida and Tonomura (2010); McMorran et al. (2011) can carry orbital angular momentum parallel to their propagation axis. Furthermore, lattices of optical OAM beams have been produced and studied Courtial et al. (2006); Vyas and Senthilkumaran (2007); Wei et al. (2009); Masajada et al. (2007). The structured OAM waves have demonstrated a number of applications in microscopy, encoding and multiplexing of communications, and manipulation of matter Padgett and Bowman (2011); Mair et al. (2001); Wang et al. (2012).

Of particular interest are “spin-orbit” beams where the orbital degree of freedom is coupled to a two-level system such as polarization for light or spin for electrons and neutrons. These beams have found applications in high resolution optical imaging, high-bandwidth communication, and optical metrology Rubinsztein-Dunlop et al. (2016); Marrucci et al. (2011); Milione et al. (2015). Spin-orbit states of light beams may be prepared by an interferometric method using a spatial light modulator Maurer et al. (2007), or via q-plates Marrucci et al. (2006). The latter method is similar to preparing spin-orbit states via a space-variant Wien filter for electrons Karimi et al. (2012) or via a quadrupole magnetic field for neutrons Nsofini et al. (2016).

The utility of the spin-orbit beams may be enhanced by producing a periodic lattice of such states, the lattice constants of which are matched to characteristic length scales of target materials. Here, we describe a universal parallel multiplexing technique that can produce a beam consisting of a lattice of OAM states coupled to a two-level system. Our protocols use coherent averaging and spatial control methods borrowed from nuclear magnetic resonance Zhang et al. (1995); Cory et al. (1989); Sodickson and Cory (1998); Levitt (2007) to prepare a general pulse sequence for producing the lattices. Spin and polarization enter here as natural manifestations of the two degrees of freedom of light and spin particles. The approach could be extended to systems with more degrees of internal freedom, such as atoms with higher spin.

Figure 1: a) The isomorphism between the Bloch sphere representing the spin states of fermions and that of the Poincaré sphere representing the polarization states of light . The corresponding eigenvectors are chosen as shown to ensure that are the transverse coordinates of the incoming beam. b) The lattices of optical spin-orbit beams are produced by passing a circularly polarized light beam through N sets of Lattice of Optical Vortices (LOV) prism pairs. A LOV prism pair consists of two perpendicular optical birefringent prisms where one prism has the optical axis along the prism incline and the second prism has the optical axis offset by . The lattice constant is given by , where is the prism incline angle, and are the extraordinary and ordinary refractive indices, and is the wavelength of the incoming light.

To describe the protocols for creating and optimizing the lattices of spin-orbit beams, we first analyze a single spin-orbit state. It is convenient to consider a (light or particle) wavepacket traveling along the -direction with momentum and expectation values of momentum in the transverse plane equal to zero. When the transverse coherence lengths are equal , where , and are the and spreads of the wavepacket’s transverse momentum distributions, the eigenstates in cylindrical coordinates can be expressed as


where is the normalization constant, is the dimensionless radial coordinate, is the azimuthal coordinate, is the radial quantum number, is the azimuthal quantum number, are the associated Laguerre polynomials, is the longitudinal wavefunction often approximated by a Gaussian wavepacket, and is the polarization state of light (or as per Fig. 1a we may use in the case of spin particles). Applying the OAM operator shows that the wavepacket carries an OAM of . The coherence length is important when dealing with particle beams where the beam is generally an incoherent mixture of coherent wavepackets, whereas for light one may simply consider the beam waist and the Laguerre-Gaussian modes. However, although the polarization-orbit beam can cleanly be described via Laguerre-Gaussian modes, the beam carrying a lattice of polarization-orbit states can not due to the translational symmetry.

When considering beams carrying OAM of major importance is the one fixed axis in space about which the OAM is quantized. In the case of beams carrying a lattice of OAM states there is a two-dimensional array of such axes and we are interested in what happens locally within each cell. Particularly, when this beam interacts with a material then the region around the local OAM axes becomes important.

To prepare states with coupled polarization and OAM we can start with circularly polarized light,


and apply a coupling operator of the form Nsofini et al. (2016)


Here, are the raising and lowering OAM operators, and are the Pauli operators, and . The length is defined as the smallest radial distance at which the polarization degree of freedom undergoes a -rotation. At radii different than , other rotation angles will occur producing the spin-orbit state


where we have set . describes a vector vortex beam where the OAM is induced via Pancharatnam-Berry geometrical phase Pancharatnam (1956); Berry (1987). The polarization distribution and the intensity post-selected on the right circularly polarized light is depicted on the zoomed-in plot of Fig. 1b. It has been shown that there is a correlation between the two degrees of freedom whereby post-selecting on one degree of freedom determines the value for the other, and that this correlation is maximized under the condition Nsofini et al. (2016).

Figure 2: Intensity profiles post-selected on a particular polarization state of the lattices of optical spin-orbit beams. Top are the simulated profiles and the bottom are the observed images. The lattice constant specified by Eq. 8 for nm light and our quartz LOV prism sets is mm; the measured lattice constant at the camera being slightly larger due to beam divergence. If desired, the lattice constant can easily be pushed into the m-range by fabricating prisms with a larger incline angle out of a high birefringent material such as TiO.

Our proposed procedure for producing a beam with a lattice of optical spin-orbit states consists of a sequence of linear birefringent gradients that are equal in magnitude and perpendicular to each other and the polarization axis of the incoming light. This procedure may be motivated by applying the Suzuki-Trotter expansion to Eq. 3a, i.e.


where we have switched from radial to Cartesian coordinates, and . Examining and truncating the right hand side of this relation, we see that it can be interpreted as a sequence of perpendicular linear gradients. Generalizing to put the origin of the gradients at and choosing that the gradients be independent of N, we define the operators


In the case of photons one way to produce the operators is via optical birefringent prisms as shown in Fig. 1b. Placing one prism with an optical axis along the prism incline and a second prism with an optical axis offset by results in the product operation . We term such a set a “Lattice of Optical Vortices (LOV) prism pair”. Eq. 6 shows that a physical shift by a distance, , of a prism along its incline direction (x or y) results in a simple phase shift of around the corresponding axis. A sequence of sets of LOV prism pairs generates a lattice of optical spin-orbit beams, calculated as

Figure 3: Phase imaging of the lattice of optical spin-orbit beams where we postselect on the polarization carrying the OAM. The sets of LOV prism pairs are placed in one path of the interferometer and a linear phase gradient is applied in the other path ( rad/mm) by tilting a mirror in order to pronounce the well known fork structure holograms in the lattice, which indicate the presence of OAM beams.

This process is shown in Fig. 1b for . The spin-orbit states in these lattices form a two-dimensional array with a lattice constant of


where and are the birefringence and the incline angle of the LOV prism pairs. The OAM structure of the resulting beam can be analyzed by looking at the phase profile of the polarization state which is correlated with the OAM:


By analyzing Eq. 9 it can be observed that the lattice cells are centered on a phase structure, while the lattice cell corners are on a structure. Although the number () of LOV prism pairs does not affect the phase profile, in any lattice cell the number of well defined intensity rings is equal to . Therefore, N provides control over the mean radial quantum number in a lattice cell, and even expansions of Eq. 5 should be used. In the case both polarization states are similarly coupled to the OAM, and both and phase structures are illuminated. Similar vortex-antivortex structures can also be obtained via Wollaston prisms Kurzynowski et al. (2006, 2010).

The simulated and observed polarization profiles for and are plotted in Fig. 2, and are in a good agreement. For our LOV prism pairs the lattice constant given by Eq. 8 comes out to be mm, though it was measured to be slightly larger due to beam divergence.

The period of the lattice can span a large range. LOV prism pairs fabricated from TiO (birefringence of ) with an incline angle of would produce a lattice period of m for a light wavelength of 532 nm. Furthermore, if birefringent materials which exhibit the Pockel’s effect are used then with the addition of external electric field control a variable period may be obtained via the electro-optic effect.

The doughnut structure in the spin-orbit states shown in Fig. 2 is indicative of the polarization profile of the polarization-orbit state (Eq. 4) and not due to the OAM structure. To show that there is a lattice of OAM states we measure the phase profile of the beam using an interferometer. The schematic of the setup is shown in Fig. 3 where a linear phase gradient in one path has been introduced to observe the characteristic fork structure hologram indicative of OAM. A lattice of fork structures can clearly be seen, indicating an at each lattice center.

Lattices of spin-orbit states may be obtained by orienting the first prism of the LOV prism pairs along the negative -direction. While various sequences of LOV prism pairs and polarization filters may be used to achieve higher order OAM structures in the outgoing beam. For example, to increment the OAM values to which the polarization states are coupled to by an integer “m”, the following sequence may be used:


where is the operator for a polarization filter along the direction. Lastly, using LOV prism pairs which produce different lattice constants results in a “superlattice” which has an overlay of the distinct lattice constants.

The described protocols provide a two-dimensional control of the characteristic length scale of the single spin-orbit features. It may be possible to create a lattice of ring-shaped optical atomic traps, individual instance of which have figured prominently in recent studies of atomtronic circuits Eckel et al. (2014). One can also envisage vortex pinning in Bose-Einstein condensates via these beams Tung et al. (2006). Lattices of polarization coupled optical vortices may also be fruitful in microscopy or basic studies of the interaction of structured light Andersen et al. (2006); He et al. (1995); Friese et al. (1996); Brullot et al. (2016) with individual atoms or molecules. This is because OAM is defined with respect to a single axis perpendicular to the wavefront. Thus, in studies using a single OAM axis, only atoms or molecules in the region of a fraction of a wavelength about that axis are subject to the OAM selection rules Afanasev et al. (2016); Schmiegelow et al. (2016). This technique extends those rules across a region proportional to the area of the fully-structured wavefront.

Our technique is particularly useful for matter-wave beams where the beam is generally an incoherent mixture of coherent wavepackets. In the case of spin particles, to create a lattice of spin-orbit states one requires a magnetic prism set with the magnetic field along the direction of the prism incline, and where the prisms are perpendicular to each other and the spin state of the incoming particles. Matter-wave lattices of spin-orbit beams may thus be generated where the OAM axis is specified along the coherent wavepacket rather than the beam axis. This opens the door for new types of studies of chiral and topological materials via particle beams.

I Methods

A laser of wavelength 532 nm was used, along with standard polarizers, wave-plates, and optical components. The LOV prism pairs were circular quartz wedges (birefringence of ) with a wedge angle of and diameter of cm. One wedge had the optical axis aligned with wedge angle while the other wedge had the optical axis aligned to wedge angle.

For images shown in Fig. 2 the setup consisted of a laser, a linear polarization filter, a quarter-wave plate, LOV prism pairs, a quarter-wave plate, a linear polarization filter, and a CMOS camera. For beam phase imaging shown in Fig. 3, a four-mirror interferometer was used because it allowed for compensation of the beam deviation due to the LOV prism pairs. An alternative method would have been to add a non-birefringent prism after each prism of the LOV prism pair in order to compensate for the beam deviation. A linear phase gradient in Fig. 3 was introduced to obtain the fork structure holograms by tilting the mirror of the interferometer path which did not contain the LOV prism pairs.

Ii Acknowledgements

This work was supported by the Canadian Excellence Research Chairs (CERC) program, the Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery program, Collaborative Research and Training Experience (CREATE) program, the Canada First Research Excellence Fund (CFREF), and the National Institute of Standards and Technology (NIST) Quantum Information Program. The authors thank Connor Kapahi for his contributions to the figures.


  • Bazhenov et al. (1990) V. I. Bazhenov, M. V. Vasnetsov,  and M. S. Soskin, Pisma v Zhurnal Eksperimentalnoi i Teoreticheskoi Fiziki 52, 1037 (1990).
  • Allen et al. (1992) L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw,  and J. P. Woerdman, Phys. Rev. A 45, 8185 (1992).
  • Harris et al. (2015) J. Harris, V. Grillo, E. Mafakheri, G. C. Gazzadi, S. Frabboni, R. W. Boyd,  and E. Karimi, Nature Physics 11, 629 (2015).
  • Garces-Chavez et al. (2002) V. Garces-Chavez, D. McGloin, H. Melville, W. Sibbett,  and K. Dholakia, Nature 419, 145 (2002).
  • Siviloglou et al. (2007) G. Siviloglou, J. Broky, A. Dogariu,  and D. Christodoulides, Physical Review Letters 99, 213901 (2007).
  • Barnett et al. (2017) S. M. Barnett, M. Babiker,  and M. J. Padgett, Philosophical Transactions of the Royal Society A 375 (2017).
  • Molina-Terriza et al. (2007) G. Molina-Terriza, J. P. Torres,  and L. Torner, Nature Physics 3, 305 (2007).
  • Clark et al. (2015) C. W. Clark, R. Barankov, M. G. Huber, D. G. Cory,  and D. A. Pushin, Nature 525, 504 (2015).
  • Sarenac et al. (2016) D. Sarenac, M. G. Huber, B. Heacock, M. Arif, C. W. Clark, D. G. Cory, C. B. Shahi,  and D. A. Pushin, Optics Express 24, 22528 (2016).
  • Uchida and Tonomura (2010) M. Uchida and A. Tonomura, nature 464, 737 (2010).
  • McMorran et al. (2011) B. J. McMorran, A. Agrawal, I. M. Anderson, A. A. Herzing, H. J. Lezec, J. J. McClelland,  and J. Unguris, science 331, 192 (2011).
  • Courtial et al. (2006) J. Courtial, R. Zambrini, M. R. Dennis,  and M. Vasnetsov, Optics express 14, 938 (2006).
  • Vyas and Senthilkumaran (2007) S. Vyas and P. Senthilkumaran, Applied optics 46, 2893 (2007).
  • Wei et al. (2009) G.-X. Wei, L.-L. Lu,  and C.-S. Guo, Optics Communications 282, 2665 (2009).
  • Masajada et al. (2007) J. Masajada, A. Popiolek-Masajada,  and M. Leniec, Optics express 15, 5196 (2007).
  • Padgett and Bowman (2011) M. Padgett and R. Bowman, Nature Photonics 5, 343 (2011).
  • Mair et al. (2001) A. Mair, A. Vaziri, G. Weihs,  and A. Zeilinger, Nature 412, 313 (2001).
  • Wang et al. (2012) J. Wang, J.-Y. Yang, I. M. Fazal, N. Ahmed, Y. Yan, H. Huang, Y. Ren, Y. Yue, S. Dolinar, M. Tur, et al., Nature Photonics 6, 488 (2012).
  • Rubinsztein-Dunlop et al. (2016) H. Rubinsztein-Dunlop, A. Forbes, M. Berry, M. Dennis, D. L. Andrews, M. Mansuripur, C. Denz, C. Alpmann, P. Banzer, T. Bauer, et al., Journal of Optics 19, 013001 (2016).
  • Marrucci et al. (2011) L. Marrucci, E. Karimi, S. Slussarenko, B. Piccirillo, E. Santamato, E. Nagali,  and F. Sciarrino, Journal of Optics 13, 064001 (2011).
  • Milione et al. (2015) G. Milione, M. P. Lavery, H. Huang, Y. Ren, G. Xie, T. A. Nguyen, E. Karimi, L. Marrucci, D. A. Nolan, R. R. Alfano, et al., Optics letters 40, 1980 (2015).
  • Maurer et al. (2007) C. Maurer, A. Jesacher, S. Fürhapter, S. Bernet,  and M. Ritsch-Marte, New Journal of Physics 9, 78 (2007).
  • Marrucci et al. (2006) L. Marrucci, C. Manzo,  and D. Paparo, Phys. Rev. Lett. 96, 163905 (2006).
  • Karimi et al. (2012) E. Karimi, L. Marrucci, V. Grillo,  and E. Santamato, Physical review letters 108, 044801 (2012).
  • Nsofini et al. (2016) J. Nsofini, D. Sarenac, C. J. Wood, D. G. Cory, M. Arif, C. W. Clark, M. G. Huber,  and D. A. Pushin, Phys. Rev. A 94, 013605 (2016).
  • Zhang et al. (1995) Y. Zhang, W. Maas,  and D. Cory, Molecular Physics 86, 347 (1995).
  • Cory et al. (1989) D. Cory, J. Miller, A. Garroway,  and W. Veeman, Journal of Magnetic Resonance (1969) 85, 219 (1989).
  • Sodickson and Cory (1998) A. Sodickson and D. G. Cory, Progress in nuclear magnetic resonance spectroscopy 33, 77 (1998).
  • Levitt (2007) M. H. Levitt, Composite pulses (Wiley Online Library, 2007).
  • Pancharatnam (1956) S. Pancharatnam, in Proceedings of the Indian Academy of Sciences-Section A, Vol. 44 (Springer, 1956) pp. 247–262.
  • Berry (1987) M. V. Berry, Journal of Modern Optics 34, 1401 (1987).
  • Kurzynowski et al. (2006) P. Kurzynowski, W. A. Wozniak,  and E. Fraczek, Applied optics 45, 7898 (2006).
  • Kurzynowski et al. (2010) P. Kurzynowski, W. Woźniak,  and M. Borwińska, Journal of Optics 12, 035406 (2010).
  • Eckel et al. (2014) S. Eckel, J. G. Lee, F. Jendrzejewski, N. Murray, C. W. Clark, C. J. Lobb, W. D. Phillips, M. Edwards,  and G. K. Campbell, Nature 506, 200 (2014).
  • Tung et al. (2006) S. Tung, V. Schweikhard,  and E. A. Cornell, Physical review letters 97, 240402 (2006).
  • Andersen et al. (2006) M. F. Andersen, C. Ryu, P. Cladé, V. Natarajan, A. Vaziri, K. Helmerson,  and W. D. Phillips, Physical Review Letters 97, 170406 (2006).
  • He et al. (1995) H. He, M. Friese, N. Heckenberg,  and H. Rubinsztein-Dunlop, Physical Review Letters 75, 826 (1995).
  • Friese et al. (1996) M. Friese, J. Enger, H. Rubinsztein-Dunlop,  and N. R. Heckenberg, Physical Review A 54, 1593 (1996).
  • Brullot et al. (2016) W. Brullot, M. K. Vanbel, T. Swusten,  and T. Verbiest, Science advances 2, e1501349 (2016).
  • Afanasev et al. (2016) A. Afanasev, C. E. Carlson,  and A. Mukherjee, Journal of Optics 18, 074013 (2016).
  • Schmiegelow et al. (2016) C. T. Schmiegelow, J. Schulz, H. Kaufmann, T. Ruster, U. G. Poschinger,  and F. Schmidt-Kaler, Nature communications 7 (2016).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description