arxiv: v1 [physics.optics] 8 May 2018

Transcription

1 Orbital Angular Momentum Induced by Optical Elements Dengke Zhang, Xue Feng, and Yidong Huang Tsinghua National Laboratory for Information Science and Technology, Department of Electronic Engineering, Tsinghua University, Beijing 18, China arxiv:185.37v1 [physics.optics] 8 May 18 INTRODUCTION Aim to manipulate orbital angular momentum OAM of light, there is a great interest in designing various optical elements and the operating waveband is ranging from microwave to deep ultra-violent light. While the light beam passing through optical elements, the carried OAM can be varied accordingly. Usually, there are two contributions from dynamic and geometric phases, respectively. And the dynamic contribution by optical elements has been clearly recognized. However, the contribution of geometric phase still seems obscure, especially when considering the vector vortex beams. In this work, an explicit formula is deduced to describe the OAM variation after the light beam passed the optical elements. It is found that the geometric contribution can be further divided into two parts: one is directly related to optical elements and the other one explicitly only relies on the vortices themselves before and after the transformations. Based on it, the same OAM variation can be achieved with different combinations of controlling the dynamic and/or geometric contributions from optical elements. Through some numerical simulations, we have demonstrated that transformation of the optical vortices can be fully and flexibly designed with a family of optical elements. Light can carry spin and orbital angular momentum SAM and OAM, which are corresponding to the polarization and spatial degrees of freedom within both classical and quantum context [1 3]. In the limit of paraxial approximation, spin and orbital components of the angular momentum can be decomposed within isotropic homogeneous media []. Each photon possesses ± the reduced Planck s constant of SAM depending on the handedness of the circular polarization, but OAM of optical vortex becomes intriguing even under paraxial approximation. For a scalar vortex beam, OAM per photon can be read as l for the spiral wavefront of the field with an azimuthal phase expilφ, where l could be any integer [5]. However, for vector vortex beams, there is not an explicit form for the OAM charge due to geometric contribution induced by a space-variant state of polarization SOP. To address it, several approaches have been proposed to extract the geometric contribution with highorder Poincar spheres [6, 7] or introducing topological Pancharatnam charge [8, 9]. In modern optics, booming materials and structures have provided more and more possibility and options to manipulate light beams or photons so that the spin-orbit interactions SOIs is inevitable and have been observed while light propagating [1, 11], focusing [1], etc. The SOIs have evoked the various investigations of physical phenomena such as the spin-hall effect [13 15], extraordinary momentum states [16] and even extended to cavity-quantum electrodynamics CQEDs [17]. Thus, it is a topic of resurgent interest to investigate the interactions between SAM and OAM in inhomogeneous or anisotropic media. Beyond dynamic evolution, the concept of geometric phases necessarily underpins the SOI of light and should be seriously considered [18, 19]. Actually, tailoring OAM of light with SOI have to be performed by various optical elements to meet the demand of various applications. Some pioneering works have been reported that a scalar vortex can acquire OAM via optical elements e.g. spiral phase plates [, 1], q-plates [, 3], J-plates [], where the spiral wavefront are originated from the SOI contributing to dynamic phase and/or geometric phase. But the mechanism and process of SOI would be much more complicated while considering vector vortex beams passing through optical elements. Consequently, a question, which has not been addressed yet to the best of our knowledge, is whether the OAM variation is distinguishable for contributions from dynamic or geometric phases for arbitrary vortex beams. And does the design of optical elements to tailor OAM have variety options as steering SAM? These issues are much more appealing for both the theoretical understanding and practical application. In this work, we have tackled both issues. Specifically, we provide an intuitive picture of OAM variation introduced by optical elements, which offers much deeper insights into how the dynamic and geometric phases evolve therein. As a concrete example, we have further developed the design rule for transformations from one scalar vortex to another to demonstrate the flexibility to achieve the same OAM variation. Then a case of design for transformations on vector vortex is presented. To that end, we have discussed the function of reported optical elements within our theoretical approach and give a summary. RESULT Theoretical principle Under the paraxial approximation, an electric field of a fully polarized paraxial vector vortex beam with the angular frequency ω and propagating along z direction can be written as [5]

2 E x, y = iω αˆx + βŷ + i α k x + β ẑ e ikz, 1 y where αβ represent the normalized complex amplitude of xy component of electric field, as functions of x, y omitted for simplicity. Then the SOP of this light at each site can be described by a 1 Jones vector a = a x, a y T, where a x a y = αβ/ I E represent the normalized complex amplitude and I E = α + β is the electric intensity of light. The corresponding Stokes vector S =S 1, S, S 3 T is defined by S j = a σ j a j = 1,, 3, where σ j are the Pauli matrices and S = a σ a, where σ equals to identity matrix [6]. Thus, S = 1 for fully polarized light and S 3 = ±1 for left/right circularly polarized field e ± = 1/ 1, ±i T. Then plotting S directly in three-dimensional Cartesian coordinates, the Poincar sphere is constructed and the corresponding azimuth ψ S and ellipticity χ S angles of SOP can be resolved by tanψ S = S /S 1 and sinχ S = S 3 /S [7]. With the aforementioned notation, the average OAM charge for a fully polarized paraxial vector vortex beam can be calculated by OAM density j o z as [3, 8] l = j o z rdrdφ ωɛ IE S rdrdφ. Next, similar as our previous work [9], the OAM density can be expressed by introducing the Pancharatnam connection between two different SOPs [9]. Here, circularly polarized fields are adopted as reference fields and then the phase difference for any field a = a x, a y T can be written as ψ P± = arg e ± a. According to Ref. [9], the OAM density can be obtained as j o z = S 3 ψ P+ ωɛ I E + 1 S 3 ψ P. 3 The detailed deductions of equations and 3 also can be found in section S1 of Supplementary Materials. In equation 3, the derivative of ψ P± is known as topological Pancharatnam charge [8, 9]. With equations and 3, the average OAM charge carried by light beam can be expressed with the SAM S 3 and topological Pancharatnam charge ψ P± /, which can depict the OAM of vortices on a single Poincar sphere as Ref. [9, 3]. Thus, the corresponding geometric phase for any transformations can be conveniently identified on the same Poincar spheres. As shown in the following, our approach provides an explicit and convenient way to design the proper optical element for manipulating the OAM charge and investigate the evolution of light beam propagating in an optical system. Figure 1. Manipulating both SAM and OAM. For the same transformation from a to b, different linear operations i.e. with different J i, where i = 1, can induce distinct OAM charge since there are different combinations of dynamic phases ψ i D and geometric phases Ω i G/ i = 1,. The magnitude of Ω i G depends on the eigen-polarization q i 1 and birefringent phase difference ψ i B of J i. The inset shows that one linear operations can transform a pair of orthogonal SOP scalar vortices a ± to another pair of SOP scalar vortices b ± with different OAM variations at the same time. Here, let us consider a scenario that a light beam passes through a non-absorbing optical element. The SOPs of input and output field are denoted as a = a x, a y T and b = b x, b y T, respectively, and the optical element is characterized by a unitary Jones matrix J, i.e., J = J 1. Thus, the light field transformation can be described as b = J a see Fig. 1. Mathematically, the eigenvalues and eigenstates of J are µ 1 and q 1, respectively, and the corresponding Stokes vectors for eigenstates can be calculated as S J = S q1 = S q = S1, J S, J S3 J T, where Sj J = q 1 σ j q 1. With these notations, the variation of OAM density can be deduced according to equation 3. For the transformation from a to b, beyond SAM variation from S3 a to S3, b there is also variation for ψp± a to ψp± b, where the superscript ab refers to the parameters related to state a b. According to Refs. [6, 31], the phase difference ψp± a b = ψb P± ψa P± can be rewritten as ψ a b P± = ψ D + ΩJ abb a ΩC abe ±, where ψ D = argµ 1 µ / presents the dynamic phase as the light beam / propagating through the optical element, Ω C abe ± is the induced geometric phase due to varied SOP between the output and input states, which corresponds to parallel transport of the state around a closed loop a b e ± a on the Poincar sphere see Fig. A, Ω J abb a is a spherical quadrangle corresponding to the closed trajectory a b b J a J a, as blue area shown in Fig. B also

3 3 see fig. S1, where a J b J holds the Stokes vector S a J b J = S ab / S ab S J S J. It can be found that the term Ω J abb a is the geometric phase explicitly related with the optical / element J. It / should be noticed that although Ω C abe ± and Ω J abb a are terms related to the geometric phase, they would affect the final OAM density with different manner. To distinct expressing the contribution of optical element and the impact of varied SOP between the output and input states, the variation of OAM density can be deduced with equation as followed form j o z = ψ D ωɛ I E + + Ω J abb a [ S3 a ψs a ψ Sb S b 3 Ω C abb a ], 5 where Ω C abb a is a spherical quadrangle defined by states a, b, b C and a C as green area shown in Fig. B, where a C b C holds the Stokes vector Sa C b C = S ab S ab S e+ S e+. It is easy to find Ω C abb a = Ω C abe + +Ω C abe see section S of Supplementary Materials for details. According to equation 5, the variation of OAM charge can be attributed to three terms. The first term C d = ψd is dynamic contribution and presents the OAM variation induced by the dynamic phase delay that only depends on ψ D of J for the optical elements, regardless the SOP of input beam. The rest two terms present the geometric contributions C g that rely on not only optical element but also the SOP of light beams. Specifically, the second term Cg J = 1 Ω J abb a is related to eigen-polarization S J i.e. q 1 and birefringent phase difference ψ B equals argµ 1µ of the considered optical element J see corresponding spherical quadrangle Ω i G in Fig. 1 or ΩJ abb a in Fig. B. The third rest term Cg V presents the geometric contribution stemming from the different SOP of input and output vortices, which explicitly depends on the inputoutput fields themselves. Namely Cg V can be fully determined by S a and S b for Ω C abb a, see Fig. B. Note that there is C g Cg J + Cg V. It means that, for an input vortex and the targeted output vortex with a given OAM variation, the geometric contribution Cg V is also determined. Thus, we can design the optical element with either or both of dynamic C d and geometric Cg J contributions to satisfy the desired OAM variation. Actually, equation 5 provides great flexibility in designing the optical element to transform one vortex beam to another. However, it has not been fully perceived and discovered at present. To demonstrate the mentioned above, some simulations have been carried out for both the scalar and vector vortex beams. Transformations on scalar vortices A Figure. Solid angle associated with geometric phase in the transformation. A Geodesic triangle Ω C abe + reddish and Ω C abe buff on the Poincar sphere. B Spherical quadrangle Ω J abb a blue on the Poincar sphere, which is a portion of lune of dihedral angle defined by states a, b and q 1. The state a J b J is a reflection of state a b referring to mirror plane of a great circle, which is perpendicular to the axis joining the states q 1 and q. Similarly, spherical quadrangle Ω C abb a green is a portion of lune of dihedral angle defined by states a, b and e ±. The state a C b C is a reflection of state a b referring to mirror plane of the equator. Firstly, we demonstrate the point-to-point PP transformation on Poincar sphere for a scalar vortex beam with the same OAM variation but different design, as sketched in Fig. 1. As a scalar vortex, the input light can be expressed as SOP of {ψ S, χ S } and OAM charge of l. For sake of simplicity but without loss of generality, ψ S = is settled since the absolute azimuth angle is not relevant due to the rotation symmetry of the coordinate. Thus, the input scalar vortex can be expressed as a + = e ilφ cosχ S, i sinχ S T. The output vortex b + is considered as b + = e il+ lφ cosχ S, i sinχ S T with a flipped handedness and OAM variation l after a optical element, which is described by Jones matrix of J see the inset of Fig. 1. It is easy to find that, for such a transformation b + = J a +, the term e ilφ can be canceled, which means the transformation is independent of OAM charge of input beam. Generally, J would linearly birefringent without considering chiral or magnetooptic materials. Thus the eigenstates of J are two linear and orthogonal eigen-polarizations. Correspondingly, the eigenvalues of J are given by e iψd {e iψb/, e iψb/ } with dynamic phase delay ψ D and birefringent phase difference ψ B after considering the unitary nature of J. And the orthogonal eigen-polarizations can be written as Rψ R 1, T and Rψ R, 1 T, where R is the standard rotation matrix and ψ R is the orientation angle of linear eigen-polarizations. Thus, for the considered optical elements, there are three parameters of {ψ D φ, ψ B φ, ψ R φ} to determine the transformation matrix J for each φ. It should be noticed that there are only two equations b + = J a + to define the relations of such three parameters. Thus, there would be various strategies to set the J for the same transformation. In B

4 D Output.P B Output Input C Input Output.P1 A Figure 3. Transformations on scalar vortices. A The design parameters {ψd, ψb, ψr } of Jones matrices for optical plates P1P from the top down. B The contributions for OAM variations from dynamic term Cd and geometric term Cg for input SOP {ψs, χs } = {, 5 }. Note that we set IE φ = 1. C SOP trace on Poincar sphere top panel and spatial distribution of SOP middle and bottom panels for input left panel and output right panel fields with a designed PP transformation. For the same input SOP but different OAM charges l = middle panel or l = 1 bottom panel, they can be transferred to both the same SOP and OAM variation by any of P1P. D SOP trace on Poincar sphere left panel and spatial distribution of SOP left panel for input fields not with designed SOP {ψs, χs } = {, 8 } top panel. The transferred fields are not the same due to different attributes of dynamic P1, middle panel and geometric P, bottom panel contributions. other words, once the dynamic phase delay ψd φ is set, a combination of {ψb φ, ψr φ} can always be found. Obviously, it is a family of optical plates to perform the PP transformation that only relies on input SOP but independent on the carried OAM charge details are discussed in section S3 of Supplementary Materials. To demonstrate this unique feature, some simulations have been performed. In Fig. 3, four different optical plates denoted as P1P are shown to transform left-handed elliptical polarization {ψs, χs } = {, 5 } vortex to right-handed elliptical polarization {ψs, χs } = {, 5 } vortex with l = 1. It can be found that CgV = for this PP transformation. Figure 3A gives the parameters of J for each optical plate and Fig. 3B shows the corresponding dynamic term Cd and geometric term CgJ that would induce the OAM variations CgV =. For P1 and P, the OAM variation is induced only by dynamic or geometric contributions, respectively. For P1, both ψb and ψr keep constant while ψd keeps constant for P. As a comparison, both dynamic and geometric terms would contribute to the OAM variation for P3 and P. For P3, the dynamic and geometric contributions are designed as homogeneous distribution along φ while inhomogeneous distribution for P. Obviously, P provides an example of more general and flexible design. Moreover, the corresponding SOP on Poincar sphere and the electric field distribution have been calculated for each optical plate to verify the equivalence of transformation. As expected, the final results are the same for P1P as shown in Figure 3C. It also can be found that the same OAM variation can be obtained with input of l = or l = 1 for P1P. It coincides with that the OAM variation is only related to SOP of input light and independent on the initial OAM charge. It should be mentioned that if the SOP of input light changes, P1P will introduce different transformations since the geometric contribution is heavily dependent on the SOP of the input beam. Thus, different design of optical plates will bring a diverse variation of SAM and OAM when the SOP of input beam does not match the designed SOP. Figure 3D shows the transformations for input plane-wave with polarization {ψs, χs } = {, 8 } by P1 and P. For P1, the output will also be a scalar vortex with l = 1 since there is the only pure geometric contribution as shown in Fig. 3D. However, the average OAM variation would be l =.87 for P and the output would be a vector vortex field as shown in Fig. 3D. The reason is that both two geometrical terms would contribute to the OAM variation see fig. S. It should be mentioned that, for the orthogonal input SOPs antipodal points on Poincar sphere, there are equal but opposite geometric contributions since they have opposite evolution direction on Poincar sphere. Thus, for the orthogonal input SOPs, the same dynamic term Cd and opposite geometric term Cg would be introduced so that

5 5 A Input Output A B SPP-I SPP-II J-plate q-plate B C 1 1 Figure. Transformations on vector vortex. A SOP trace on Poincar sphere and spatial distribution of SOP for input left panel and output right panel vector vortices with a designed transformation by any of optical plates P5P7. B The design parameters {ψ D, ψ B, ψ R} of Jones matrices for optical plates P5P7 from the top down. C The contributions for OAM variation from dynamic term C d, geometric term Cg J and Cg V for input vector vortex shown in A from the top down. The contribution from Cg V are the same for any of P5P7. Note that we set I Eφ = 1. l C d ± C g. Transformations on vector vortices More generally, equation 5 can be applied to transformations on vector vortices. As shown in Fig. A, the input beam is a cylindrical vortex with χ S = 3 and l = while the output beam is another cylindrical vortex with χ S = 3 and l = 1 see fig. S3 for {ψ S, χ S } of the input and output beams. For such a transformation, the linearly birefringent unitary J is employed to design the optical plates. By numerically solving the J, the design parameters for three different optical plates named P5P7 are displayed in Fig. B. The corresponding dynamic and geometric contributions are given in Fig. C. It can be found that although the terms of Cg V are the same but not equal to zero, proper designs of optical elements are needed. Similar to PP transformations, we have designed P5 with only dynamic contribution C d and P6 with only geometric contribution Cg J. Moreover, both two terms are designed for P7. It should be noticed that there is no theoretical limitation for applying equation 5 on designing optical element. However, in reality, it is not easy to achieve arbitrary transformation on vector vortices since the physically implemented Jones matrices would be restricted by the available materials and structures. Figure 5. Features of spiral phase plates SPPs, J- plates, and q-plates. A Design parameters {ψ D, ψ B, ψ R} of Jones matrices from the top down. For SPP-I and SPP- II, the desired OAM variation is l and the portion from dynamic contribution is set as t for SPP-II. For J-plate, the desired OAM variation is l 1/ for input light with left/righthanded circular polarization. For qplate, the desired OAM variation is ± l for input light with left/right-handed circular polarization. B The corresponding dynamic upper panel and geometric lower panel contributions for OAM variation. For J-plate and q-plate, the opposite geometric contributions will be induced according to handedness of input light. DISCUSSION So far, there are three common optical platesspiral phase plates, q-plates, and J-plates, which have been reported to generate and manipulate OAM beams [,, 3, 33]. Actually, for all of them, the operation mechanism can be more clearly understood within our theoretical approach. Here, we will give some discussions and comments about them. For spiral phase plates SPPs, there are two types according to equation 5. The first one SPP-I is fabricated with homogeneous materials and the required phase delay is induced by the spiral design of plate []. Thus, there is only the pure dynamic contribution so that the same dynamic phase, as well as the same OAM variation, can be obtained in spite of the SOP of input beams. For the second type SPP-II, the desired OAM variation is achieved only for a specific SOP of the input beam and the output beam would possess the same SOP [1]. Thus, SPP-II as a particular case of PP transformation, the SOP of the input beam is just coincided with one of eigen-polarizations of the plate to remain the same on SOP. But it just works on a specific SOP, which means it is not a pure dynamic contribution based optical plate. Thus, for SPP-II, both the dynamic and geometric contributions have to be taken into account for the OAM variations case by case. Obviously, the operation mechanism of SPP-II is totally different with SPP-I since not only pure dynamic contribution

6 6 would introduce the OAM variation. For the J-plates [], it could be considered as a special type of PP transformation. The special constraint is that the orthogonal SOPs of input beam should be transferred to a flipped handedness with a different OAM variation at the same time. As shown in the inset of Fig. 1, for a J-plate, the OAM variation l 1 and l should be obtained for b + = J a + and b = J a, respectively and simultaneously. According to our framework, it means that the dynamic contribution is always the same but the geometric contribution is opposite for orthogonal SOPs. Thus, the J-plate can be designed as that the dynamic OAM variation is l 1 + l / and opposite geometric OAM variation should be ± l 1 l / according to the handedness of SOP of the input beam. Then the combinations of three parameters {ψ D, ψ B, ψ R } can be readily obtained. Particularly, if the input field is circular polarization χ S = ±π/, three parameters would hold simple forms as shown in Fig. 5A. This kind of J-plate is a half-plate and can flip circular polarizations with different OAM variations. Specifically, if l 1 = l = l, it is, in fact, the well-known q-plate, in which only the pure geometric contribution is introduced [, 3, 33]. The full parameters of J for these mentioned optical plates and the corresponding contributions of each term are summarized and presented in Fig. 5. In this work, an explicit form to evaluate the OAM variation due to the optical elements have been presented with both dynamic and geometric phases. With the help of the topological Pancharatnam charge, the geometric phases can be further separated into two parts. One is directly related to optical elements and the other one explicitly only relies on SOPs of the input and output lights. This treatment is not just a mathematical trick but would introduce a new viewpoint to fully understand the operation mechanism and explore the designing flexibility of an optical element. Based on this, targeted applications could be developed, for instance, pure dynamic contribution based optical plates can achieve identical OAM variations in spite of the SOP of input beam while pure geometric contribution based optical plates can serve as a mode sorter for both SAM and OAM in modern optical communication systems. Moreover, our theoretical approach can be employed for other optical systems to analyze influences due to the dynamic and geometric phases. For the simulation results, only the case of linear orthogonal eigen-polarizations of Jones matrix is considered since it is the common response of most materials and structures. It should be mentioned that our theoretical approach is not limited by this constraint. Actually, if the eigen-polarizations of the Jones matrix could be arbitrary, more complicated function can be achieved for various potential applications. Additionally, there are several assumptions in our theoretical deduction such as unitary Jones matrix, paraxial beam, and fully polarized fields. Actually, breaking either of them would introduce some more interesting investigations, e.g. considering inhomogeneous Jones matrix [3], non-hermitian including PT -symmetry systems [35], or non-reciprocal systems [36]. Furthermore, only classic beams are considered in our work, but we believe that the similar work about quantum counterpart would bring more things of new physics and our work could evoke some fundamental research about spin-orbit interaction and the related fields. MATERIALS AND METHODS Vortex beams in numerical simulations For general vector beams, such as cylindrical vortices, the field can be expressed as E φ = 1 π cos χ S ˆx + iŷe illφ ψs + 1 π sin χ S ˆx iŷe ilrφ+ψs, 6 where l L and l R are topological charges of field components with left and right-handed circular polarization, respectively. With equation 6, there is I E φ = 1. For a scalar vortex beam with topological charge l, it is easy to be obtained by setting l L = l R = l. For PP transformation shown in Fig. 3C, we set the input light with SOP of {ψ S, χ S } = {, 5 } and topological charge of l L = l R = or 1. For the transformation shown in Fig. 3D, we set the input light with {ψ S, χ S } = {, 8 } and l L = l R =. And for the transformation on vector vortex shown in Fig., we set the input light with SOP of {ψ S, χ S } = {π φ, π/6} and topological charge of {l L, l R } = {1, 3}, and the detailed SOPs of input and output vector vortices are presented in fig. S3. This work was supported by the National Key Research and Development Program of China 17YFA337, the National Natural Science Foundation of China Grant No dkzhang@outlook.com x-feng@tsinghua.edu.cn [1] Padgett, M., Courtial, J. & Allen, L. Light s orbital angular momentum. Phys. Today 57, 35. [] Franke-Arnold, S., Allen, L. & Padgett, M. Advances in optical angular momentum. Laser Photon. Rev., [3] Barnett, S. M. Optical angular-momentum flux. J. Opt. B Quantum Semiclassical Opt., S7. [] Yao, A. M. & Padgett, M. J. Orbital angular momentum: Origins, behavior and applications. Adv. Opt. Photonics 3, [5] Allen, L., Beijersbergen, M. W., Spreeuw, R. J. C. & Woerdman, J. P. Orbital angular momentum of light and the

7 7 transformation of Laguerre-Gaussian laser modes. Phys. Rev. A 5, [6] Milione, G., Sztul, H. I., Nolan, D. A. & Alfano, R. R. Higher-order Poincaré sphere, Stokes parameters, and the angular momentum of light. Phys. Rev. Lett. 17, [7] Milione, G., Evans, S., Nolan, D. A. & Alfano, R. R. Higher order Pancharatnam-Berry phase and the angular momentum of light. Phys. Rev. Lett. 18, [8] Niv, A., Biener, G., Kleiner, V. & Hasman, E. Manipulation of the Pancharatnam phase in vectorial vortices. Opt. Express 1, 8 6. [9] Zhang, D., Feng, X., Cui, K., Liu, F. & Huang, Y. Identifying orbital angular momentum of vectorial vortices with Pancharatnam phase and Stokes parameters. Sci. Rep. 5, [1] Bliokh, K. Y., Niv, A., Kleiner, V. & Hasman, E. Geometrodynamics of spinning light. Nat. Photonics, [11] Bliokh, K. Y., Bekshaev, A. Y. & Nori, F. Optical momentum, spin, and angular momentum in dispersive media. Phys. Rev. Lett. 119, [1] Zhao, Y., Edgar, J., Jeffries, G., McGloin, D. & Chiu, D. Spin-to-orbital angular momentum conversion in a strongly focused optical beam. Phys. Rev. Lett. 99, [13] Onoda, M., Murakami, S. & Nagaosa, N. Hall effect of light. Phys. Rev. Lett. 93, [1] Hosten, O. & Kwiat, P. Observation of the spin Hall effect of light via weak measurements. Science 319, [15] Bliokh, K. Y., Smirnova, D. & Nori, F. Quantum spin Hall effect of light. Science 38, [16] Bliokh, K. Y., Bekshaev, A. Y. & Nori, F. Extraordinary momentum and spin in evanescent waves. Nat. Commun. 5, [17] Deng, Y., Cheng, J., Jing, H. & Yi, S. Bose-Einstein condensates with cavity-mediated spin-orbit coupling. Phys. Rev. Lett. 11, [18] Souza, C. E. R., Huguenin, J. A. O., Milman, P. & Khoury, A. Z. Topological phase for spin-orbit transformations on a laser beam. Phys. Rev. Lett. 99, [19] Ma, L. B. et al. Spin-orbit coupling of light in asymmetric microcavities. Nat. Commun. 7, [] Oemrawsingh, S. S. R. et al. Production and characterization of spiral phase plates for optical wavelengths. Appl. Opt. 3, [1] Yu, N. et al. Light propagation with phase discontinuities: Generalized laws of reflection and refraction. Science 33, [] Marrucci, L., Manzo, C. & Paparo, D. Optical spin-toorbital angular momentum conversion in inhomogeneous anisotropic media. Phys. Rev. Lett. 96, [3] Bomzon, Z., Biener, G., Kleiner, V. & Hasman, E. Spacevariant Pancharatnam-Berry phase optical elements with computer-generated subwavelength gratings. Opt. Lett. 7, [] Devlin, R. C., Ambrosio, A., Rubin, N. A., Mueller, J. P. B. & Capasso, F. Arbitrary spin-to-orbital angular momentum conversion of light. Science 358, [5] Torres, J. P. & Torner, L. Twisted Photons: Applications of Light with Orbital Angular Momentum John Wiley & Sons, 11. [6] Gutiérrez-Vega, J. C. Pancharatnam-Berry phase of optical systems. Opt. Lett. 36, [7] Born, M. & Wolf, E. Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light Elsevier, 13. [8] Allen, L. & Padgett, M. J. The Poynting vector in Laguerre-Gaussian beams and the interpretation of their angular momentum density. Opt. Commun. 18, [9] Berry, M. The adiabatic phase and Pancharatnam s phase for polarized light. J. Mod. Opt. 3, [3] Zhao, P. et al. Identifying the tilt angle and correcting the orbital angular momentum spectrum dispersion of misaligned light beam. Sci. Rep. 7, [31] Martínez-Fuentes, J. L., Albero, J. & Moreno, I. Analysis of optical polarization modulation systems through the Pancharatnam connection. Opt. Commun. 85, [3] Biener, G., Niv, A., Kleiner, V. & Hasman, E. Formation of helical beams by use of Pancharatnam-Berry phase optical elements. Opt. Lett. 7, [33] Karimi, E. et al. Generating optical orbital angular momentum at visible wavelengths using a plasmonic metasurface. Light Sci. Appl. 3, e [3] Cerjan, A. & Fan, S. Achieving arbitrary control over pairs of polarization states using complex birefringent metamaterials. Phys. Rev. Lett. 118, [35] Lawrence, M. et al. Manifestation of PT symmetry breaking in polarization space with terahertz metasurfaces. Phys. Rev. Lett. 113, [36] Mahmoud, A. M., Davoyan, A. R. & Engheta, N. Allpassive nonreciprocal metastructure. Nat. Commun. 6,

8 1 SUPPLEMENTARY MATERIALS for Orbital Angular Momentum Induced by Optical Elements ORBITAL ANGULAR MOMENTUM OF AN OPTICAL VORTEX Under the paraxial approximation, the electric and magnetic fields of a fully polarized vectorial vortex beam of angular frequency ω propagate along z direction can be written as [5] E x, y = iω B x, y = ik αˆx + βŷ + i α k x + β y βˆx + aŷ + i β k x + α y ẑ e ikz, e ikz, ẑ S1a S1b where α and β represent the complex amplitude of of x and y component of electric field, respectively. They can be written as αx, y = I E x, ya x x, y = I E x, ya x x, ye iδxx,y, βx, y = I E x, ya y x, y = I E x, ya y x, ye iδyx,y, Sa Sb where a x a y = αβ/ I E are normalized electric field components with electric intensity of I E = α + β, A xy and δ xy are real numbers and represent amplitude and phase, respectively. Meanwhile, we define δ s = δ y δ x is the phase difference between x and y components of electric field. So the polarization state of this light at each site can be described by a 1 Jones vector a = a x, a y T. Then, Stokes vector S = S 1, S, S 3 T is defined by S j = a σ j a j = 1,, 3, where σ j are the Pauli matrices [6], σ 1 = 1, σ 1 = 1, σ 1 3 = i. i Meanwhile S = a σ a, where σ equals to identity matrix. Thus. for left/right circularly polarized light e ± = 1/ 1, ±i T i.e. δ s = ±π/, there is S 3 = ±1. Then plotting Stokes vector S on three-dimensional Cartesian coordinates, the Poincar sphere could be constructed and the corresponding azimuth ψ S and ellipticity χ S angles are resolved by [7], respectively tanψ S = S /S 1, sinχ S = S 3 /S. S3a S3b The linear momentum density, which is defined as of p = ε E B, can be expressed and divided into transverse and longitudinal components p = i ωε [α α + β β α α β β + α β β α ẑ], Sa p z = ωkε α + β = ωkε I E S. Sb Meanwhile, the energy density of such a beam is w = cp z = ε ω α + β = ε ω I E S. S5 Then, the cross product of line momentum density with r radius vector gives the angular momentum density, and z component of angular momentum density is

9 j z = r p z = i ωε Further, j z can be divided into spin and orbital parts as [ α α + β β α α β β + r ] r α β β α. S6 j s z = iωε r r α β β α = ωε r I ES 3, S7a r j o z = i ωε α α + β β α α β β = ωε I E A δ x x + A y δ y. S7b With the ratio of angular momentum over energy that is examined by Allen [8], the average SAM charge and OAM charge can be calculated as j s s = ω z rdrdφ IE S 3 rdrdφ = wrdrdφ IE S rdrdφ, j o l = ω z rdrdφ j o = z rdrdφ wrdrdφ ωɛ IE S rdrdφ. S8a S8b Then, introducing the phase difference for two different state of polarization SOP of e A and e B through the Pancharatnam connection, which is defined by [9] ψ P = arg e A e B. S9 Here, using left/right circularly polarized fields as reference fields, the phase difference for any field a = a x, a y T can be written as ψ P± = arg e ± a. S1 According to Ref. [9], we can obtain j o z ψ P± = S ωɛ I E ± S S 3 ψ S, S11 then using equation S11, we can get a relation then equation S11 can be rewritten as j o z ωɛ I E = 1 ψ S = 1 ψp+ ψ P, S1 S + S 3 ψ P+ + S S 3 ψ P. S13 Substituting equations S11 or S13 into equation S8b, we can calculate the average OAM charge for any vortex beams. In equation S13, the derivative of ψ P± is known as topological Pancharatnam charge. OAM VARIATION INDUCED BY OPTICAL ELEMENTS Here, we consider a scenario that the light beam passes a non-absorbing optical element and investigate the OAM variation due an optical element, which is characterized by a unitary Jones matrix J, i.e., J = J 1 with the eigenvalues of µ 1, eigenstates of q 1, and the corresponding Stokes vectors of S J = S q1 = S q = S J 1, S J, S J 3 T S J j = q 1 σ j q 1 [6]. When light beam denoted as a passes through optical element J, the output beam is expressed as b = J a. From equation S8b, it could be known that the variation of OAM only depends the variation of j o z due to non-absorbing nature I a E = Ib E = I E. Then with equations S1 and S13, the variation of j o z can be deduced as

10 j o z ωɛ I E = 1 S b 3 S a 3 ψa P+ + Sb + S3 b ψa b P+ + S a 3 S b 3 ψa P 3 + Sb S3 b ψa b P, S1 where the superscript ab refers to the parameters related to state a b and ψp± a b = ψb P± ψa P± is the difference of phases defined by the Pancharatnam connection between a and b. According to Refs. [6, 31], ψp± a b can be written as ψ a b P± = ψ D ΩC abe ± + ΩJ abb a, S15 where/ ψ D = argµ 1 µ / is dynamic phase gained by the beam when it propagates through the optical element, Ω C abe ± is the geometric phase, related to the referenced circularly polarized field, which corresponds to parallel transport of/ the state around a closed loop a b e ± a on the Poincar sphere see Fig. A in main text, and Ω J abb a is the geometric phase introduced by the optical element J. For the third term, Ω J abb a is a spherical quadrangle corresponding to the closed trajectory a b b J a J a, as shown in Fig. B in main text, where a J b J holds the Stokes vector Sa J b J = S ab S ab S J S J. Further, using equations S1 and S15, the equation S1 can be rewritten as j o z = S b ψ D ωɛ I E + Ω J Sb abb a + [ S3 a ψs a ψ Sb S b 3 Ω C Sb abb a ], S16 where Ω C abb a is a spherical quadrangle defined by states a, b, b C and a C as shown in Fig. b in main text, where a C b C holds the Stokes vector Sa C b C = S ab S ab S e+ S e+. Thus, using equation S8b the variation of OAM charge can be solved l = IE j o zrdrdφ ωɛ IE S a rdrdφ. S17 Note: For any spherical triangle defined by states a, b and c, whose Stokes vectors are S a, S b and S c, respectively, the triangular area Ω abc is [ S a S b S c ] Ω abc = arctan 1 + S a S b + S b S c + S c S a. S18 From equation S18, it can found that clockwise and anticlockwise walk on the sphere surface will induce opposite value of solid angles. As shown in Fig. B in main text, the spherical lune formed by two geodesics connecting the antipodal states q 1 and q passing through a and b and forming a dihedral angle ψ B = argµ argµ 1 = argµ 1µ, which is introduced by birefringent of optical element also see fig. S1. The corresponding lune area equals to ψ B = Ω abq1 Ω abq. In particular, when q 1 coincides with e ±, ψ B also equals to ψs b ψa S. It is easy to find that there is a relation Ω C abb a = Ω abq1 + Ω abq and Ω J abb a can be solved by similar approach. PP TRANSFORMATION ON SCALAR VORTEX The input scalar vortex a + is set as polarization azimuth of ψ S, ellipticity of χ S and carrying OAM of l a + = e ilφ cosχs Rψ S, S19 i sinχ S where R is the standard rotation matrix. Then using an optical element with Jones matrix J transferring it to the output vortex b + with a flipped handedness and OAM charge m as b + = e imφ cosχs Rψ S, S i sinχ S then we can obtain the relation e imφ Rψ S cosχs = e ilφ cosχs JRψ i sinχ S S. S1 i sinχ S

11 For a scalar vortex, due to the rotation symmetry of the coordinate choice for polarization azimuth, the ψ S = is settled for simplicity but without loss of generality, thus equation S1 can be rewritten as e i l1φ cosχs cosχs = J, S i sinχ S i sinχ S where l 1 = m l is the variation of OAM from a + to b +. Generally, without considering chiral or magneto-optic materials, J is linearly birefringent and the eigenvalues will correspond to linearly polarized orthogonal eigen-polarizations. Being unitary, eigenvalues of J are given by complex exponentials of e iψd {e iψb/, e iψb/ } for dynamic phase delay ψ D and birefringent induced phase difference ψ B. And such two orthogonal eigen-polarizations can be written as Rψ R 1, T and Rψ R, 1 T, where ψ R is the orientation angle of eigen-polarizations. Then, for such kinds of optical elements, the Jones matrix J can be expressed as J = J1 J = e iψd J 3 J cos ψb i sin i sin ψb ψb cos ψ R i sin ψb sin ψ R sin ψ R cos ψb + i sin ψb cos ψ R From equation S3, it can be found that there are there parameters {ψ D φ, ψ B φ, ψ R φ} to determine the Jones matrix for each φ since J = J 3, while there are just two equations to define their relations by equation S. The result is that we have infinite choices to construct J to achieve the same transformation. That is once a contribution from dynamic phase delay ψ D φ is set, we always can find a selection of {ψ B φ, ψ R φ}. Combining equations S and S3, we can rewrite J J = e iψd cos l1 φ ψ D cosχ S + i sin l 1 φ ψ D i cos l 1 φ ψ D sinχ S i cos l 1 φ ψ D sinχ S cos l 1 φ ψ D cosχ S i sin l 1 φ ψ D. S3. S With both S3 and S, for any input light with SOP of χ S to obtain a desired OAM variation l 1, the design of three parameters {ψ D φ, ψ B φ, ψ R φ} can be solved for J of optical element. To present the effect from {ψ D, ψ B } of the optical element with eigen-polarizations q 1, there are some comments and discussions on the considered transformation. As shown in fig. S1, for the input state a, we can use a geodesic arc join q 1, a and q and let the arc go a rotation of ψ B around the axis defined by its eigen-polarizations q 1. Then the final state b can be obtained on the corresponding location as shown in fig. S1. From such transformation, the introduced dynamic phase is always equal to ψ D, while the introduced geometric phase is Ω G / = Ω J abb a /. It is easy to find that the geometric phase equals to ψ B / if a = q 1 and ψ B / if a = q. Moreover, if the input SOP is orthogonal to a, i.e. the antipodal point on the Poincar sphere, the evolution will encircle on an opposite direction but the same area using the same optical element, this means that the introduced geometric phase is Ω G / for input light with orthogonal SOP. J-plates and q-plates For the J-plates [], there is a special constrain on PP transformation. That is, the input field holding with an orthogonal SOP and OAM charge k gives a = e ikφ i sinχs Rψ S, S5 cosχ S and the transferred output field also with a flipped handedness and OAM charge n, yielding b = e inφ i sinχs Rψ S. S6 cosχ S With the similar approach and setting ψ S =, we can find another relation for J as e i lφ i sinχs i sinχs = J, S7 cosχ S cosχ S where l = n k is the variation of OAM from a to b. Then, combining equations S and S7, we can get the design parameters of J. For a special and simple case of input field with circular polarization χ S = ±π/, the combination of equations S and S7 can reduce J to

12 J 1 = J = 1 e i l 1φ e i lφ J = J 3 = i. S8 e i l 1φ + e i lφ From equation S8, we can find the eigenvalues as i l1 + l φ µ 1 = i exp, S9 i l1 + l φ µ = i exp 5 and eigen-polarizations as l1 l φ π 1 q 1 = R, S3 l1 l φ π q = R 1 thus the parameters {ψ D φ, ψ B φ, ψ R φ} for J can be found as ψ D φ = l 1 + l φ ψ B φ = π ψ R φ = l 1 l φ π. S31 It can be found that this kind of J-plates is a half-plate but possesses spatial variable dynamic phase delays and orientation angles and can transfer scalar vortices with circular polarizations. Further, q-plates would be as a special case of J-plates with opposite variation of OAM, i.e. l 1 = l. Thus, for the q-plate, there is only the contribution from geometric phases while none from dynamic phases. It should be noted that, for J-plates and q-plates, the geometric contributions actually come from the third term of equation S16. And to calculate the OAM variation, the ψ S should be calculated using equation S1 while not equation S3a due to there is a singularity for circularly polarized fields. And this singularity also induces the geometric contribution is transformed to Cg V while not Cg J i.e., not same with the demonstration of PP transformation in main text. However, this phenomenon occurs just because the circularly polarized fields is set as the reference fields and it can be resolved if another pair of reference fields is adopted. Spiral phase plates SPPs Type I: This type of spiral phase plates SPP-I is fabricated with homogeneous materials and can provide required phase delay by designing path length. In the transformation, only dynamic phases ψ D φ should be considered since ψ B φ =. From equation S3, the corresponding Jones matrix for this kind of spiral phase plates can be written 1 J = e iψd. S3 1 So they can not change the SOP but always introduce the same dynamic phase for any SOP. That is to say that the same OAM variation can be achieved for any input fields. Type II: There is another kind of spiral phase plates SPP-II can provide a desired OAM variation only for a specific SOP of input light beam and the output also possesses the same SOP. Actually, the input SOP is just coincided with one of eigen-polarizations so that the dynamic or geometric or both of them would contribute to the OAM variations. For simplicity, we set the input SOP as {ψ S, χ S } = {, } and then ψ R φ =. Thus, with equation S3, we can obtain the Jones matrix e iψ B/ J = e iψd e iψb/. S33

13 It can found that if ψ D φ = and all the contribution comes from geometric phase. For this case, the orthogonal input SOP will get an opposite OAM variation. There is an extreme case of ψ B φ =, which is exactly the Type I plates. Except for such the two extreme cases, this type of plates actually is a common PP transformation plate with same input and output SOPs. 6

14 7 Figure S1. For the transformation from a to b, there will be the introduced dynamic phase ψ D and geometric phase Ω G/ by an optical element of {ψ D, ψ B} and eigen-polarizations q 1. Ω G = Ω J abb a is a spherical quadrangle corresponding to the closed trajectory a b b J a J a, where a J b J holds the Stokes vector Sa J b J = S ab S ab S J S J E in.ψ S E in.χ S E out.ψ S E out.χ S C d C J g C V g C d + C g φ/π φ/π Figure S. Transformation with P for input light with SOP of {ψ S, χ S} = {, 8 }. A SOPs for input and output fields. B Detailed contributions for OAM variarion from each term. E in E out E in E out ψs/π χs φ/π φ/π Figure S3. SOPs for input and output vector vortices shown in Fig. in main text. A Azimuth angle ψ S of the vector vortices. B Ellipticity angle χ S of the vector vortices.