Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Generation and dynamics of optical beams with polarization singularities

Open Access Open Access

Abstract

We present a convenient method to generate vector beams of light having polarization singularities on their axis, via partial spin-to-orbital angular momentum conversion in a suitably patterned liquid crystal cell. The resulting polarization patterns exhibit a C-point on the beam axis and an L-line loop around it, and may have different geometrical structures such as “lemon”, “star”, and “spiral”. Our generation method allows us to control the radius of L-line loop around the central C-point. Moreover, we investigate the free-air propagation of these fields across a Rayleigh range.

© 2013 Optical Society of America

1. Introduction

It has been widely studied and directly observed that both longitudinal and transverse waves may have structures rich of singular points [1, 2]. Optical light beams with singular scalar and vectorial features received a particular attention [2]. In the scalar case, these singularities in the beam transverse plane are known as “wave dislocations” or “optical vortices” [1, 3]. They correspond to points where the optical field amplitude vanishes and the optical phase becomes undefined. Near the singular point, the optical phase changes continuously with a specific topological charge, given by the phase change in units of 2π along a closed path surrounding the singularity [3]. Beams with optical vortices carry orbital angular momentum (OAM), which is simply proportional to the topological charge when the vortex is located on the beam axis [4].

When also the vectorial nature of the optical field is taken into account, the variety of possible singular points in the transverse plane is increased [511]. In the paraxial approximation, the optical field is purely transverse, i.e. orthogonal to the propagation direction. In each point of the transverse plane, the direction of the electric field changes in time tracing an ellipse, whose ellipticity and orientation define the local “polarization” state of the field. The polarization can be also parametrized by introducing the three reduced Stokes parameters, s1, s2 and s3, or the polar and azimuthal angles on the Poincaré sphere. The vector-field singularities form points or lines in the transverse plane where one of the polarization parameters is undefined. Examples are the so-called C-points and L-lines, where the orientation and the handedness of the polarization ellipse are undefined, respectively [6,7]. Optical beams exhibiting vector singularities of this kind have been given various names in the literature, such as for example (full)-Poincaré beams [12, 13]. Here, we simply call them “ polarization-singular beams” (PSB). The PSB may have different transverse polarization patterns, characterized by a different integer or half-integer topological index η as defined by the variation of the polarization ellipse orientation for a closed path around a C-point, such as “star”, “monstar”, “lemon”, “spiral” (with a central “source/sink” singularity), “hyperbolic” (with a central “saddle” singularity), etc. Some examples are shown in Fig. 1. Optical PSB possessing high values of positive and negative topological charge form so-called “polarization flowers” and “hyperbolic webs”, respectively [8].

 figure: Fig. 1

Fig. 1 Polarization patterns of three different polarization singular beams: “star” (index η = −1/2), “lemon” (η = 1/2) and “spiral” (η = 1). The central point at the origin indicates the C-point, while dashed blue lines show the L-line circles with radius ρ0. The handednesses of the polarization ellipses in the outer region, ρ > ρ0, and in the inner one, ρ < ρ0, are opposite.

Download Full Size | PDF

Polarization singularities are typically generated either by interferometric methods (see, e.g., [14]), which is usually very delicate and unstable, or by propagation through suitable birefringent solid crystals (see, e.g., [15] and references therein), a method which does not lend itself to a flexible control of the output. In this paper, we study the generation and propagation dynamics of vectorial polarization singular beams generated by “q-plates” [16], i.e. birefringent liquid crystal cells having a singular pattern of the optical axis. In a recent work, we have reported the generation of singular vortex beams with uniform polarization ellipticity and a central vortex singularity with vanishing intensity, as generated by “tuned” q-plates [17]. Here, by exploiting “untuned” q-plates, we conveniently generate and control beams having a space-variant polarization ellipticity, with a central C-point singularity, and a lemon, star or spiral ellipse-orientation pattern. Finally, we analyze the varying polarization structures in different longitudinal positions along the beam, so as to investigate the effect of propagation.

2. Polarization singular beams

A class of polarization singular beams can be generated by the superposition of two optical fields bearing orthogonal polarization states and different OAM values. Here, we consider in particular a beam described by the following expression:

|PSB=cos(δ2)LG0,0(ρ,ϕ,z)|L+eiαsin(δ2)LG0,m(ρ,ϕ,z)|R
where ρ, ϕ, z are cylindrical coordinates, δ is a tuning parameter, |L〉 and |R〉 stand for left- and right-circular polarizations, respectively, α is a phase, and LGp,m denotes a Laguerre-Gauss mode of radial index p and OAM index m, with a given waist location (z = 0) and radius w0.

The LG modes with m ≠ 0 have a doughnut shape, with vanishing intensity at the center. The PSB resulting from Eq. (1) then have a polarization pattern, shown in Fig. 1, characterized by the following four elements: (i) at the center (that is the beam axis), the polarization is left-circular (C-point), since the last term in Eq. (1) vanishes here; (ii) as the radius increases, the second term contributes and changes the circular polarizations into left-hand elliptical; (iii) at a certain radius ρ = ρ0(δ), where the two waves in Eq. (1) have equal amplitude, the polarization becomes exactly linear (L-lines); (iv) at larger radii, that is for ρ > ρ0, the second term in Eq. (1) prevails and the polarization state changes into right-handed elliptical. It is worth noting that, in this case, the polarization pattern has a topological index η = m/2. The initial angle at azimuthal location ϕ = 0 of the polarization ellipse major axis is fixed by the relative phase α in Eq. (1). Of course, the |L〉 and |R〉 polarizations in Eq. (1) can also be swapped, leading to an inverted polarization ellipticitiy everywhere.

We generate PSB by exploiting partial spin-to-orbital angular momentum conversion in q-plates. The optical action of q-plates is described in [16, 1820]. Depending on its birefringent retardation δ, the q-plate converts a varying fraction of the spin angular momentum change suffered by the light passing through the device into OAM with m = ±2q, where q is the q-plate topological charge and the sign is determined by the input circular polarization handedness. The conversion is full when the q-plate optical phase retardation is δ = π (or any odd multiple of π). For other values of δ, the q-plate is said to be untuned and the spin-to-orbital conversion is only partial. In the untuned q-plate a fraction of the incoming beam passes through unchanged. The field at the exit plane of an untuned q-plate for a Gaussian left-circular polarized input beam is given by the following expression:

U^|L=HyGG|q|2,q2(ρ,d/n¯)[cos(δ2)|L+eiαsin(δ2)|Re2iqϕ],
where HyGG|q|2,q2 denotes the amplitude profile of a hypergeometric-Gaussian beam [21], describing with good approximation the output mode profile of a q-plate of charge q, thickness d and average refractive index for a Gaussian input (see Ref. [22] for more details), and the phase α in the second term depends on the initial angle (on the x axis) of the q-plate optical-axis pattern [16,18]. The converted and non-converted terms in Eq. (2), that is in the near field, have a common amplitude profile, which can therefore be factorized. The far-field profiles produced by the converted and non-converted parts of the beam are different, however, because the different OAM values affect the beam propagation. The HyGG term appearing in Eq. (2), after multiplication by the phase factor exp(2iqϕ), can be expanded in a series of LG modes [21]. Truncating this expansion to the first useful order (which is a good approximation in the far field), Eq. (2) is reduced to Eq. (1) with m = 2q. The optical retardation δ of the q-plate, that is the q-plate tuning, can be adjusted continuously by applying an appropriate voltage [20]. Different polarization topologies can be obtained by using q-plates with different topological charges. In this work, we used q = −1/2, +1/2, and +1, generating correspondingly equal indices of the polarization patterns.

3. Experiment

The experiment layout is shown in Fig. 2. The generated PSB polarization distribution was studied both in the near and far fields generated in the image and focal planes of a lens. The polarization distribution was analyzed by polarization tomography: the local Stokes parameter maps have been calculated from the intensity images taken by a CCD camera (Sony AS-638CL), after projection onto horizontal (H), vertical (V), anti-diagonal (A), diagonal (D), L and R polarization states, respectively. The polarization pattern was thus entirely reconstructed, as shown in form of little ellipses in Fig. 2 and subsequent figures. In order to reduce noise and experimental errors, the average intensity was recorded in a grid of 20×20 pixel square area [17]. The missing ellipses in the experimentally reconstructed polarization patterns (e.g., in Fig. 3) correspond to regions with low intensity, so that the polarization reconstruction was too affected by stray light. The superposition coefficients in Eq. (2) have been adjusted by changing the q-plate optical retardation δ with the applied voltage. The radius where the L-line forms in the polarization pattern depends on the ratio between the absolute value of the two coefficients in Eq. (2) and, hence, on the voltage applied to the q-plate. In particular, for optical retardation δ = 0 and δ = π the beam profile is Gaussian fully left-circular and doughnut-shaped fully right-circular, respectively. In the intermediate case, the pattern is a superposition and a circular L-line appears in the polarization pattern, whose radius ρ0 depends on δ. From Eq. (1) one finds ρ0=w0/(2tanδ/2). Figure (3) shows the intensity and reconstructed polarization pattern for different q-plate retardations in steps Δδ = π/4, and the corresponding measured L-lines radii. As shown in the figure, also the intensity profile changes from a Gaussian (left-circular) to a doughnut (right-circular) shape. In Fig. 3 (c) and (g), the optical retardations are π/2 and 3π/2, respectively; the π difference in the relative phases results into a 180° rotation of the polarization patterns. The topological index of the polarization pattern, however, depends only on the q-value of the q-plate, and in the case of Fig. 3 we have η = q = 1/2.

 figure: Fig. 2

Fig. 2 Experimental setup. A He-Ne laser beam is first spatial-mode filtered by focusing via a microscope objective (4X) into a 50 μm pinhole and recollimating by a lens (f). Its polarization is then prepared in a left (right) circular state by a polarizing beam-splitter (PBS) and a quarter-wave plate (Q). An untuned q-plate then transforms the beam into a PSB. The PSB polarization pattern was analyzed by a quarter-wave plate, a half-wave-plate (H) and another polarizer, followed by imaging on a CCD camera. The δ parameter of the PSB was adjusted by electrically tuning the q-plate. In order to study the propagation dynamics of the PSB, the CCD camera was mounted on a translation stage and moved around the focal plane of an output lens. Upper insets show the intensity and reconstructed polarization patterns, in the near field, of the lemon, star, and spiral beams generated by q-plates with q = 1/2, −1/2, 1, respectively.

Download Full Size | PDF

 figure: Fig. 3

Fig. 3 Intensity distribution and reconstructed polarization patterns of beam generated by a q-plate with q = 1/2 for seven different optical retardations: (a) δ = 0 (or 2π), (b) δ = π/4, (c) δ = π/2, (d) δ = 3π/4, (e) δ = π, (f) δ = 5π/4, and (g) δ = 3π/2. The corresponding L-line radii relative to the beam waist w0 are the following: (a) undefined (b) ρ0 = 2.0w0, (c) ρ0 = 1.4w0, (d) ρ0 = 1.1w0, (e) undefined, (f) ρ0 = 1.0w0, and (g) ρ0 = 1.5w0.

Download Full Size | PDF

In a second experiment, we studied the dynamics of different PSBs with η = −1/2, +1/2 and +1 topologies under free-air propagation for fixed parameter δ = π/2. To this purpose, we moved the CCD camera along the propagation axis around the imaging-lens focal plane. We recorded the beam polarization patterns at six different planes in the range −zRzzR, where zR is the lens Rayleigh parameter and z = 0 corresponds to the beam waist location. The experimental results are shown in Fig. 4. As it can be seen, the polarization pattern evolves during propagation. Indeed LG00 and LG0,m modes of Eq. (1) have different z dependences of their Gouy phases; the relative Gouy phase between them is given by ψ = |m|arctan(z/z0), which in the explored region varies in the range |m|[−π/4, π/4]. In the case m = ±1 (i.e. η = 1/2) the phase evolution leads to a rigid rotation of the whole polarization structure by an angle equal to ψ [as shown in Fig. 4 (a) and (b)], while when m = 2 (i.e. η = 1), this dephasing leads to a more complex pattern dynamics, in which the polarization structure changes from radial to spiral and then to azimuthal (Fig. 4 (c)).

 figure: Fig. 4

Fig. 4 Reconstructed experimental polarization patterns of different PSB beams. (a) η = m/2 = +1/2, (b) η = m/2 = −1/2 and (c) η = m/2 = +1. Patterns have been reconstructed by measuring the maps of reduced Stokes parameters in six different longitudinal planes within the beam Rayleigh range, from −zR to +zR. The corresponding rotation of polarization patterns for (a) and (b) are 90° and 88° = (30 + 30 + 28) °, respectively.

Download Full Size | PDF

4. Conclusion

We reported a non-interferometric technique to generate beams carrying a polarization singularity on their axis, based on a partial spin-to-orbital angular momentum conversion in a detuned q-plate device. By this approach, we generated beams with C-points on their optical axis, surrounded by a circular L-line whose radius can be controlled by the voltage applied to the q-plate. We demonstrated different patterns such as stars, lemons, and spirals, corresponding to vector topological indices of −1/2, 1/2 and +1, respectively. We have also studied the pattern evolution in free-space propagation. Such polarization patterns could find practical application in lithographic optical techniques, as for the examples discussed in Refs. [23, 24].

Acknowledgments

We thank Sergei Slussarenko for preparing the q-plates. We acknowledge financial support of the Future and Emerging Technologies FET-Open programme, within the 7th Framework Programme for Research of the European Commission, under grant number 255914 - PHOR-BITECH.

References and links

1. J. F. Nye and M. V. Berry, “Dislocations in Wave Trains,” Proc. R. Soc. Lond. A 336, 165–190 (1974). [CrossRef]  

2. M. R. Dennis, K. O’Holleran, and M. J. Padgett, “Optical Vortices and Polarization Singularities,” Prog. Opt. 53, 293–363 (2009). [CrossRef]  

3. M. S. Soskin and M. V. Vasnetsov, (Ed. E. Wolf), “Singular optics,” Prog. Opt.42, 219–276 (2001). [CrossRef]  

4. S. Franke-Arnold, L. Allen, and M. J. Padgett, “Advances in optical angular momentum,” Laser Photonics Rev. 2, 299–313 (2008). [CrossRef]  

5. J. F. Nye, “Polarization effects in the diffraction of electromagnetic waves: the role of disclinations,” Proc. Roy. Soc. Lond. A 387, 105–132 (1983). [CrossRef]  

6. J. F. Nye, “Lines of circular polarization in electromagnetic wave fields,” Proc. Roy. Soc. Lond. A 389, 279–290 (1983). [CrossRef]  

7. M. R. Dennis, “Polarization singularities in paraxial vector fields: morphology and statistics,” Opt. Commun. 213, 201–221 (2002). [CrossRef]  

8. I. Freund, “Polarization flowers,” Opt. Commun. 199, 47–63 (2001). [CrossRef]  

9. O. Angelsky, A. Mokhun, I. Mokhun, and M. Soskin, “The relationship between topological characteristics of component vortices and polarization singularities,” Opt. Commun. 207, 57–65 (2002). [CrossRef]  

10. M. V. Berry, M. R. Dennis, and R. L. Lee Jr., “Polarization singularities in the clear sky,” New J. Phys. 6, 162 (2004). [CrossRef]  

11. M. S. Soskin, V. Denisenko, and I. Freund, “Optical polarization singularities and elliptic stationary points,” Opt. Lett. 28, 1475–1477 (2003). [CrossRef]   [PubMed]  

12. A. M. Beckley, T. G. Brown, and M. A. Alonso, “Full Poincare beams,” Opt. Express 18, 10777–10785 (2010). [CrossRef]   [PubMed]  

13. E. J. Galvez, S. Khadka, W. H. Schubert, and S. Nomoto, “Poincaré-beam patterns produced by nonseparable superpositions of Laguerre-Gauss and polarization modes of light,” Appl. Opt. 51, 2925–2934 (2012). [CrossRef]   [PubMed]  

14. P. Kurzynowski, W. A. WoŸniak, M. Zdunek, and M. Borwińska, “Singularities of interference of three waves with different polarization states,” Opt. Express 20, 26755–26765 (2012). [CrossRef]   [PubMed]  

15. T. A. Fadeyeva, C. N. Alexeyev, P. M. Anischenko, and A. V. Volyar, “Engineering of the space-variant linear polarization of vortex-beams in biaxially induced crystals,” Appl. Opt. 51, C224–C230 (2012). [CrossRef]   [PubMed]  

16. L. Marrucci, C. Manzo, and D. Paparo, “Optical Spin-to-Orbital Angular Momentum Conversion in Inhomogeneous Anisotropic Media,” Phys. Rev. Lett. 96, 163905 (2006). [CrossRef]   [PubMed]  

17. F. Cardano, E. Karimi, S. Slussarenko, L. Marrucci, C. de Lisio, and E. Santamato, “Polarization pattern of vector vortex beams generated by q-plates with different topological charges,” Appl. Opt. 51, C1–C6 (2012). [CrossRef]   [PubMed]  

18. L. Marrucci, E. Karimi, S. Slussarenko, B. Piccirillo, E. Santamato, E. Nagali, and F. Sciarrino, “Spin-to-orbital conversion of the angular momentum of light and its classical and quantum applications,” J. Opt. 13, 064001 (2011). [CrossRef]  

19. E. Karimi, B. Piccirillo, E. Nagali, L. Marrucci, and E. Santamato, “Efficient generation and sorting of orbital angular momentum eigenmodes of light by thermally tuned q-plates,” Appl. Phys. Lett. 94, 231124 (2009). [CrossRef]  

20. S. Slussarenko, A. Murauski, T. Du, V. Chigrinov, L. Marrucci, and E. Santamato, “Tunable liquid crystal q-plates with arbitrary topological charge,” Opt. Express 19, 4085–4090 (2011). [CrossRef]   [PubMed]  

21. E. Karimi, G. Zito, B. Piccirillo, L. Marrucci, and E. Santamato, “Hypergeometric-Gaussian modes,” Opt. Lett. 32, 3053–3055 (2007). [CrossRef]   [PubMed]  

22. E. Karimi, B. Piccirillo, L. Marrucci, and E. Santamato, “Light propagation in a birefringent plate with topological charge,” Opt. Lett. 34, 1225–1227 (2009). [CrossRef]   [PubMed]  

23. C. Hnatovsky, V. Shvedov, W. Krolikowski, and A. Rode, “Revealing Local Field Structure of Focused Ultrashort Pulses,” Phys. Rev. Lett. 106, 123901 (2011). [CrossRef]   [PubMed]  

24. A. Ambrosio, L. Marrucci, F. Borbone, A. Roviello, and P. Maddalena, “Light-induced spiral mass transport in azo-polymer films under vortex-beam illumination,” Nat. Commun. 3, 989 (2012). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (4)

Fig. 1
Fig. 1 Polarization patterns of three different polarization singular beams: “star” (index η = −1/2), “lemon” (η = 1/2) and “spiral” (η = 1). The central point at the origin indicates the C-point, while dashed blue lines show the L-line circles with radius ρ0. The handednesses of the polarization ellipses in the outer region, ρ > ρ0, and in the inner one, ρ < ρ0, are opposite.
Fig. 2
Fig. 2 Experimental setup. A He-Ne laser beam is first spatial-mode filtered by focusing via a microscope objective (4X) into a 50 μm pinhole and recollimating by a lens (f). Its polarization is then prepared in a left (right) circular state by a polarizing beam-splitter (PBS) and a quarter-wave plate (Q). An untuned q-plate then transforms the beam into a PSB. The PSB polarization pattern was analyzed by a quarter-wave plate, a half-wave-plate (H) and another polarizer, followed by imaging on a CCD camera. The δ parameter of the PSB was adjusted by electrically tuning the q-plate. In order to study the propagation dynamics of the PSB, the CCD camera was mounted on a translation stage and moved around the focal plane of an output lens. Upper insets show the intensity and reconstructed polarization patterns, in the near field, of the lemon, star, and spiral beams generated by q-plates with q = 1/2, −1/2, 1, respectively.
Fig. 3
Fig. 3 Intensity distribution and reconstructed polarization patterns of beam generated by a q-plate with q = 1/2 for seven different optical retardations: (a) δ = 0 (or 2π), (b) δ = π/4, (c) δ = π/2, (d) δ = 3π/4, (e) δ = π, (f) δ = 5π/4, and (g) δ = 3π/2. The corresponding L-line radii relative to the beam waist w0 are the following: (a) undefined (b) ρ0 = 2.0w0, (c) ρ0 = 1.4w0, (d) ρ0 = 1.1w0, (e) undefined, (f) ρ0 = 1.0w0, and (g) ρ0 = 1.5w0.
Fig. 4
Fig. 4 Reconstructed experimental polarization patterns of different PSB beams. (a) η = m/2 = +1/2, (b) η = m/2 = −1/2 and (c) η = m/2 = +1. Patterns have been reconstructed by measuring the maps of reduced Stokes parameters in six different longitudinal planes within the beam Rayleigh range, from −zR to +zR. The corresponding rotation of polarization patterns for (a) and (b) are 90° and 88° = (30 + 30 + 28) °, respectively.

Equations (2)

Equations on this page are rendered with MathJax. Learn more.

| PSB = cos ( δ 2 ) LG 0 , 0 ( ρ , ϕ , z ) | L + e i α sin ( δ 2 ) LG 0 , m ( ρ , ϕ , z ) | R
U ^ | L = HyGG | q | 2 , q 2 ( ρ , d / n ¯ ) [ cos ( δ 2 ) | L + e i α sin ( δ 2 ) | R e 2 i q ϕ ] ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.