Spin-orbit interactions in optically active materials

Spin-orbit interactions in optically active materials

Chandroth P. Jisha Centro de Física do Porto, Faculdade de Ciências, Universidade do Porto, R. Campo Alegre 687, Porto 4169-007, Portugal.    Alessandro Alberucci Optics Laboratory, Tampere University of Technology, FI-33101 Tampere, Finland alessandro.alberucci@gmail.com
July 15, 2019

We investigate the inherent influence of light polarization on the intensity distribution in anisotropic media undergoing a local inhomogeneous rotation of the principal axes. Whereas in general such configuration implies a complicated interaction between geometric and dynamic phase, we show that, in a medium showing an inhomogeneous circular birefringence, the geometric phase vanishes. Due to the spin-orbit interaction, the two circular polarizations perceive reversed spatial distribution of the dynamic phase. Based upon this effect, polarization-selective lens, waveguides and beam deflectors are proposed.

The simplest approaches in wave optics model light as a scalar wave, an approximation valid only for paraxial beams propagating in isotropic homogeneous media Lax:1975 (). The electromagnetic nature of light implies that photons have spin, appearing in the Maxwell’s equations as the field polarization. Nonetheless, in case of paraxial waves propagating in isotropic materials, the spatial degree of freedom (i.e., the field distribution) is independent of the polarization. This is not rigorously true due to the vectorial nature of the Maxwell’s equation, leading to a so-called classical entanglement Spreeuw:1998 (); Berg-Johansen:2015 (). The discrepancies between the two approaches are significant, for example, in case of non-paraxial beams leading to polarization-dependent focusing Dorn:2003 (), or for polarization-dependent trajectory of light in inhomogeneous materials, the so called optical Magnus effect, aka spin Hall effect (SHE) of light Bliokh:2004 (); Bliokh:2007 (); Bliokh:2008 (); Zhang:2016 ().

In this context, the investigation of spin-orbit interaction is rapidly becoming a central topic in optics Bliokh:2015 (). Photons have been demonstrated to be a unique tool for the investigation of basic quantum field theory Neugenbauer:2015 (); Bliokh:2015_1 (). On the other hand, spin-orbit interaction paves the way to a new family of ultra-thin photonic devices, including gratings Bomzon:2002 (), lenses Lin:2014 (), polarimeters Todorov:1992 (); BalthasarMueller:2016 (), beam shaper Marrucci:2006 (); Arbabi:2015 (), sensors Petersen:2014 (), deflectors Li:2013 () and so on.

In essence, spin-orbit interactions in isotropic materials are connected with the fact that singular plane waves composing the beam are transverse electromagnetic in different frameworks. Practical examples include beams bended by gradient index Bliokh:2008 (), interaction with interfaces Hosten:2008 () and even free space Dorn:2003 (). These effects can be interpreted in terms of a geometric phase called the Rytov-Vladimirskii-Berry phase.

In anisotropic media another type of geometric phase arises, the Pancharatnam-Berry phase (PBP). PBP appears in the presence of a rotation of the beam polarization for a fixed wave-vector. Noteworthy, in anisotropic materials spin-orbit effects, such as birefringence or spatial walk-off, are observed even for plane waves due to the dependency on polarization of the light-matter interaction Saleh:2007 (). Nonetheless, the presence of PBP in anisotropic media widens the spectrum of the observable spin-orbit effects Bomzon:2002 (); Marrucci:2006 (); Lin:2014 (), including the nonlinear case Li:2015 (); Nookala:2016 (). The PBP affects the wavefront, and thus light propagation, when the anisotropic medium is inhomogeneous and showing a point-wise rotation of the principal axes across the intensity profile. In this Letter we first analyze light propagation in inhomogeneously rotated materials considering the trade-off between the PBP and diffraction. We show that PBP vanishes if the material is optically active, the field propagation then being affected by dynamic phase alone. Unlike Refs. Bliokh:2004 (); Bliokh:2007 (); Zhang:2016 () where geometric optics is used, here the wave behavior of the electromagnetic radiation is accounted for. Differently from the case of linear birefringence Alberucci:2016 (), we demonstrate that in optically active media the two circular polarizations see a photonic potential reversed in sign, yielding polarization-dependent guiding and SHE.
Let us consider the monochromatic propagation of an electromagnetic field [amplitude ] in an inhomogeneous anisotropic material. The average wavevector is parallel to the axis . We consider a non-magnetic material (). Supposing that the medium permittivity does not vary along , the electric field in the paraxial limit obeys


where the anisotropy in the diffraction coefficients has been neglected and is the vacuum wavenumber. The paraxial hypothesis allows us to neglect the longitudinal field in (1). Finally, in writing (1) we have assumed that the input beam, polarized on the plane , does not excite any appreciable polarization parallel to , that is, the linear walk-off vanishes.
Generally speaking, the two-dimensional tensor will change on the transverse plane either due to a local change in the values of its eigenvalues and/or due to a local rotation of the principal axes. Variation of eigenvalues is associated with changes in the local refractive indices. Principal axes distribution, in our case, can be modeled like a rotation by an angle around the propagation direction , . Thus, for a generic vector it is where the subscript indicates the framework of the local principal axes. To analyze the effect of the local rotation , let us introduce the field in the rotated framework . After setting , (1) provides Alberucci:2016 ()


A straightforward computation yields


where represent the Pauli matrices.
Let us now discuss how the solutions of (Spin-orbit interactions in optically active materials) depend on the features of . From a general point of view, transparent anisotropic materials can be classified into two categories depending on the polarization of the plane waves which are supported as eigenmodes by the medium: we talk of linear or of circular birefringence when the eigenmodes are linearly or circularly polarized waves, respectively Saleh:2007 (). The last family can be easily generalized to the case of elliptic polarizations. Mathematically, the anisotropic material shows linear birefringence if the components of are symmetric, whereas circular birefringence is related with an anti-symmetric . Media fulfilling the latter condition are called gyrotropic due to the appearance of optical activity, and includes both chiral and magneto-optical materials, where the optical activity is spontaneous or induced by an external magnetic field, respectively Condon:1937 (); Saleh:2007 (). For a lossless medium, the components should be also complex conjugate for the tensors to be Hermitian, that is, .
In the presence of linear birefringence, is diagonal with two distinct eigenvalues, let us call them and . Then we can always set , where is the optical anisotropy. Substituting back into (Spin-orbit interactions in optically active materials), we find terms proportional to two different Pauli matrices, thus (Spin-orbit interactions in optically active materials) cannot be diagonalized in any basis Weinberg:2012 (). This means that the polarization of the field will always vary while evolving along , in turn leading to the occurrence of a transversely-varying PBP, whenever is not uniform. In accordance with the former statement, in the case of linear birefringence it is well known that the wavefront is strongly affected by the PBP, both for short Bomzon:2002 (); Marrucci:2006 (); Li:2013 () or long propagation Slussarenko:2016 (); Alberucci:2016 () with respect to the Rayleigh distance.
Light propagation changes drastically in the presence of circular birefringence, that is, , where is the modulus of the gyration vector. As a matter of fact, if then is proportional to , and (Spin-orbit interactions in optically active materials) can be diagonalized in the circular basis composed by left circular polarization (LCP) and right circular polarization (RCP) . Complex unit vectors and can also be rewritten as , where is the photon helicity, with and . Hence, no exchange of power between the two circular polarizations takes place versus . This behavior agrees with physical intuition: a rotation around the propagation distance does not change the polarization state for circular polarizations, the latter being eigenvectors of the rotation matrix .
Above discussion holds valid even when varies across (variations must be small on the wavelength scale, otherwise (1) is no more valid). Formally, let us introduce the slowly varying envelopes , where . Then (Spin-orbit interactions in optically active materials) can be diagonalized using the space-independent transformation between linear and circular polarization basis,


Substitution of (4) to (Spin-orbit interactions in optically active materials) yields the two uncoupled equations


Eventually, application of the polarization-dependent gauge transformation transforms Eqs. (Spin-orbit interactions in optically active materials-6) into two paraxial Helmholtz equations . There is no exchange of power between the LCP and the RCP, meaning no contribution from the PBP: only a polarization-dependent and inhomogeneous dynamic phase is acting on the electromagnetic field. Summarizing, light propagates as if subject to a polarization-dependent distribution of refractive index given by


We tested the theoretical predictions by means of in-house BPM (Beam Propagation Method) codes in (1+1)D and (2+1)D, simulating directly (1) in the paraxial approximation (second derivative along is neglected) and in the presence of a transversely-varying circular birefringence. Two types of inhomogeneity for will be considered: with reference to the (1+1)D geometry, an even-symmetric (in two dimensions cylindrically symmetric) and an anti-symmetric profile (in two dimensions mirror symmetry will be imposed). According to Eqs. (7-8), in the first case a focusing/defocusing behavior is expected according to the beam helicity, whereas in the second case a spin-dependent deflection should occur.

Figure 1: Polarization-dependent confinement. Top row: light undergoes either guiding or anti-guiding behavior depending on the input polarization of the beam; for m and RCP (left) shows focusing and LCP (right) shows defocusing. Bottom row: beam width versus the propagation distance for varying strength of optical activity and RCP input, as marked in each panel. Different curves correspond to different width of the optical activity profile: solid blue lines correspond to m, dashed red lines correspond to m and black dotted lines to m. Input beam is Gaussian with a waist of m placed in .

For simplicity, we begin with the one-dimensional case and later on extend the results to the two-dimensional case. For the even case, optical activity is supposed to be Gaussian, i.e., , where is the maximum of the optical activity and is the width of the transverse distribution. Hereafter the input beam is a fundamental Gaussian beam and a Laguerre-Gaussian (LG) beam , in the (1+1)D and (2+1)D cases respectively. We also take a wavelength of m and : the results can be easily generalized to other cases using standard normalizations. BPM results confirm the analytical results, i.e., the two polarizations evolve independently, seeing an effective index distribution provided by Eqs. (7-8) (see top panels of Fig. 1). The bottom panels of Fig. 1 depict the evolution of the beam width of the confined polarization for different values of the parameters, and , of the graded waveguide. The confining effect becomes stronger, e.g. the number of bounded modes raises up, either increasing and/or the width of the optical activity profile , as witnessed by the oscillations in the beam width. Power is almost completely coupled into the fundamental mode when m, m and (first panel in Fig. 1).

Figure 2: Even symmetric : Profiles of the input LG beam with width m at different propagation planes as marked (normalized with power at each plane). The confined beam corresponds to the RCP and the diffracting part corresponds to the LCP component of the beam. The last row shows the phase and the Stokes parameter for the beam at m. The parameters of the medium are and m.

Analogous behavior occurs in two-dimensional geometries. To achieve a complete description of light with one single simulation for each geometry, in (2+1)D input beams are taken linearly polarized, that is, a superposition with the same weight of RCP and LCP waves. The RCP and LCP will be either focused and defocused, according to the sign of . To stress out the differences with respect to birefringent media Marrucci:2006 (), we will consider beams carrying angular momentum. Accordingly, input beams are Laguerre-Gaussian beam carrying a unitary charge of orbital angular momentum (OAM), whereas for a cylindrically-symmetric Gaussian profile is set. Results are plotted in Fig. 2. The RCP component sees a focusing potential, coupling almost all its energy to guided modes of the waveguide. On the other hand, the LCP component perceives a repelling potential and spreads out faster than in a homogeneous medium. The transverse distribution of the beam phase, plotted in Fig. 2, shows that the diffracting and the confined beam carry the same OAM of the input beam. In contrast, in an anisotropic medium with transversely varying optic axis (a so-called q-plate Marrucci:2006 ()), the spin angular momentum (SAM) and the OAM are inextricably connected due to the PBP. The polarization behavior is addressed plotting the component of the Stokes vector : for a purely RCP beam , whereas for a LCP beam , with all the other Stokes components being equal for both the circular polarizations. Polarization distribution on the beam cross section , plotted in Fig. 2, shows that the confined part is RCP and the diffracting halo is LCP, fully confirming the predictions.

Figure 3: Anti-symmetric in (1+1)D: Evolution of the Stokes parameter showing the splitting of an input linearly polarized beam with m. The beam in the upper and the lower half plane are LCP and RCP, respectively. Each panel corresponds to different longitudinal length (in the figure expressed in unit of the Rayleigh length , equal to 11.8m in this case) of the medium, the latter extending from to the white lines (after free space is assumed). Here and m.

We change now the spatial symmetry of the optical activity, i.e., we consider an anti-symmetric profile for , . Depending on the input polarization, the beam will get deflected either to the left or to the right, i.e., optical SHE Zhang:2016 (). Figure 3 shows the evolution of the Stokes parameter in the (1+1)D case for a linearly polarized input. The RCP and LCP components separate spatially and are deflected symmetrically to either sides of the input direction. With increasing , the strength of deflection is increased (not shown here). Deflection angle also depends on the length of the optically active medium, see the different panels of Fig. 3. Finally, for longer propagation lengths in the optically active medium, the beam begins to develop fringes due to self-interference effects (last panel).

Figure 4: Anti-symmetric in (2+1)D: Profiles of the input LG beam with width m at different propagation planes as marked. The right moving beam corresponds to the RCP and the left moving part corresponds to the LCP component of the beam. The last row shows the phase of the beam at m and the Stokes parameter for the beam at m. The white rectangles show the position of the two vortices, each of them carrying the original topological charge. The parameters of the medium are and m.

Next, we want to check the polarization-dependent deflection for beams carrying OAM. Input is a linearly polarized LG beam with unitary OAM charge. Like the SHE in the 1D case, the vortex beam splits into two specular halves, the LCP and the RCP. The two circular polarizations, deflected by the same angle in absolute value, broaden due to diffraction as well (Fig. 4). The OAM of the two separated beams is same as that of the input, i.e., there is no conversion of SAM into OAM Marrucci:2006 ().
In conclusion, we demonstrated that, in locally twisted anisotropic materials, the Pancharatnam-Berry phase vanishes in the presence of circular birefringence. Unlike with linear birefringence, a stronger spin-orbit effect occurs due to the monotonic accumulation of the polarization-dependent effects in propagation. Although in natural media the smallness of optical activity hinders the observation of the proposed effects, man-made metamaterials with enhanced circular birefringence are ideal candidates to verify our theoretical results Plum:2007 (); Gansel:2009 (); Schaferling:2012 (); Kaschke:2015 (), with possible extensions to the nonlinear case as well Ren:2012 ().

Funding. Fundação para a Ciência e a Tecnologia, POPH-QREN and FSE (FCT, Portugal), fellowship SFRH/BPD/77524/2011. Academy of Finland FiDiPro no. 282858.


  • (1) M. Lax, W. H. Louisell, and W. B. McKnight, Phys. Rev. A 11, 1365 (1975).
  • (2) R. J. C. SpreeuwFoundations of Physics 28, 361 (1998).
  • (3) S. Berg-Johansen, F. Töppel, B. Stiller, P. Banzer, M. Ornigotti, E. Giacobino, G. Leuchs, A. Aiello, and C. Marquardt, Optica 2, 864 (2015).
  • (4) R. Dorn, S. Quabis, and G. Leuchs, Phys. Rev. Lett. 91, 233901 (2003).
  • (5) K. Y. Bliokh and Y. P. Bliokh, Phys. Rev. E 70, 026605 (2004).
  • (6) K. Y. Bliokh, D. Y. Frolov, and Y. A. Kravtsov, Phys. Rev. A 75, 053821 (2007).
  • (7) K. Y. Bliokh, A. Niv, V. Kleiner, and E. Hasman, Nat. Photon. 2, 748 (2008).
  • (8) Y. Zhang, L. Shi, and C. Xie, Opt. Lett. 41, 3070 (2016).
  • (9) K. Y. Bliokh, F. J. Rodriguez-Fortuno, F. Nori, and A. V. Zayats, Nat. Photon. 9, 796 (2015).
  • (10) M. Neugebauer, T. Bauer, A. Aiello, and P. Banzer, Phys. Rev. Lett. 114, 063901 (2015).
  • (11) A. Y. Bekshaev, K. Y. Bliokh, and F. Nori, Phys. Rev. X 5, 011039 (2015).
  • (12) Z. Bomzon, G. Biener, V. Kleiner, and E. Hasman, Opt. Lett. 27, 1141 (2002).
  • (13) D. Lin, P. Fan, E. Hasman, and M. L. Brongersma, Science 345, 298 (2014).
  • (14) T. Todorov and L. Nikolova, Opt. Lett. 17, 358 (1992).
  • (15) J. P. B. Mueller, K. Leosson, and F. Capasso, Optica 3, 42 (2016).
  • (16) L. Marrucci, C. Manzo, and D. Paparo, Phys. Rev. Lett. 96, 163905 (2006).
  • (17) A. Arbabi, Y. Horie, M. Bagheri, and A. Faraon, Nat. Nanotech. 10, 937 (2015).
  • (18) J. Petersen, J. Volz, and A. Rauschenbeutel, Science 346, 67 (2014).
  • (19) G. Li, M. Kang, S. Chen, S. Zhang, E. Y.-B. Pun, K. W. Cheah, and J. Li, Nano Lett. 13, 4148 (2013). PMID: 23965168.
  • (20) O. Hosten and P. Kwiat, Science 319, 787 (2008).
  • (21) B. E. A. Saleh and M. C. Teich, Fundamentals of Photonics (Wiley Publishing, 2007).
  • (22) G. Li, S. Chen, N. Pholchai, B. Reineke, P. W. H. Wong, E. Y. B. Pun, K. W. Cheah, T. Zentgraf, and S. Zhang, Nat. Mater. 14, 607 (2015).
  • (23) N. Nookala, J. Lee, M. Tymchenko, J. S. Gomez-Diaz, F. Demmerle, G. Boehm, K. Lai, G. Shvets, M.-C. Amann, A. Alu, and M. Belkin, Optica 3, 283 (2016).
  • (24) A. Alberucci, C. P. Jisha, L. Marrucci, and G. Assanto, ACS Photonics 3, 2249 (2016).
  • (25) E. U. Condon, Rev. Mod. Phys. 9, 432 (1937).
  • (26) S. Weinberg, ed., Lectures on Quantum Mechanics (Cambridge University Press, Cambridge, 2012).
  • (27) S. Slussarenko, A. Alberucci, C. P. Jisha, B. Piccirillo, E. Santamato, G. Assanto, and L. Marrucci, Nat. Photon. 10, 571 (2016).
  • (28) E. Plum, V. A. Fedotov, A. S. Schwanecke, N. I. Zheludev, and Y. Chen, Appl. Phys. Lett. 90, 223113 (2007).
  • (29) J. K. Gansel, M. Thiel, M. S. Rill, M. Decker, K. Bade, V. Saile, G. von Freymann, S. Linden, and M. Wegener, Science 325, 1513 (2009).
  • (30) M. Schäferling, D. Dregely, M. Hentschel, and H. Giessen, Phys. Rev. X 2, 031010 (2012).
  • (31) J. Kaschke and M. Wegener, Opt. Lett. 40, 3986 (2015).
  • (32) M. Ren, E. Plum, J. Xu, and N. I. Zheludev, Nat. Commun. 3, 833 (2012).
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