This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Also at ]Istituto di Scienze e Tecnologie per l’Energia e la Mobilità Sostenibili, Consiglio Nazionale delle Ricerche, P.le Tecchio 80, 80125 Napoli, Italy.Also at]CNR-ISASI, Institute of Applied Science and Intelligent Systems, Via Campi Flegrei 34, 80078 Pozzuoli (NA), Italy Also at]INFN, Sez. di Napoli, Complesso Universitario di Monte Sant’Angelo, via Cinthia, 80126 Napoli, Italy

Engineering monstar polarization disclination through geometric phase

Verónica Vicuña-Hernández    Pegah Darvehi [    Lorenzo Marrucci [    Bruno Piccirillo bruno.piccirillo@unina.it [ Dipartimento di Fisica “Ettore Pancini”, Università degli Studi di Napoli Federico II, Complesso Universitario di Monte Sant’Angelo, Via Cintia, 80126 Napoli, Italy
Abstract

We present a method to generate monstar singularities via Pancharatnam-Berry phase by the coherent collinear superposition of two Free-Form Dark Hollow beams, FFDHqmFFDH^{m}_{q}, of topological charge qq, and order of symmetry mm. FFDH beams are generated with the geometrical parameters of a closed curve exploited to obtain a nonuniform rotation rate of the local polarization azimuth and generate monstar patterns. We report space-variant polarization patterns: radial- and azimuthal-like, lemon- and star-like for symmetric disclinations, and the asymmetric monstar disclination. We present theory and measurements, and find excellent agreement between the two.

preprint: APS/123-QED

I Introduction

The generation of classical and quantum states of light is the result of the ability to give structure to the degrees of freedom of light. It has served and impacted in a wide variety of applications, in domains ranging from basic science to manufacturing to communications [1]. Recently, the study of topological features of light has been increased giving rise to discoveries and descriptions about wavefront dislocations [2, 3, 4], and phase and polarization singularities [5]. In particular the inhomogeneous space-variant polarization structures in paraxial and non-paraxial regimens that take place in anisotropic media have been investigated in the last years due to their riches in spin-orbit optical interactions [6]. Poincaré beams (PBs) are beams of light with inhomogeneous space-variant polarization distribution containing rotational dislocations in directional order [7]. They compass different spatial regions on the Poincaré sphere surface over their transverse plane [8]. The singular point of the dislocation is a V- or C-point. Polarization singularities arise when one of the parameters that defines the state of polarization of light is undefined. V-point singularities are related to vectors fields in which the linearly polarized optical field has both orientation and handedness undefined. Poincare-Hopf index η\eta defines the order of the V-point singularity [9]. C-point singularity occurs in ellipse fields and is associated with circular polarization in which the azimuth of the state of polarization is undefined and the orientation of the surrounded ellipses rotate clockwise or counterclockwise about it [9]. Such rotation is represented by an integer or half-integer index ICI_{C}, indicating that the azimuthal coordinate on the Poincaré sphere rotates 2|IC|2|I_{C}| times per turn about the singularity. Basic polarization disclinations are classified as lemons (IC=1/2I_{C}=1/2, 1 radial line), stars (IC=1/2I_{C}=-1/2, 3 radial lines) and monstars (IC=1/2I_{C}=1/2, 3 radial lines). Another characteristic of the disclination is the number of radial lines, L-lines. These are straight lines that originate at the singularity. The polarization orientation is radial at all points along a radial line. Lemons with IC=1/2I_{C}=1/2 and stars with IC=1/2I_{C}=-1/2 have symmetric patterns with one and three radial lines respectively. Monstar disclination is asymmetric and occurs when the azimuth rotation rate about a C-point exceeds a threshold value and produces more radial lines than lemons or stars with the same ICI_{C} index [10, 11, 12]. Different approaches have been taken in order to generate the three types of basic disclinations patterns in optical beams. In the majority of cases, the methods adopted are based on the coherent superposition of two distinct spatial modes in two orthogonal polarization states. However, distinct approaches have been exploited depending on the actual tools used for molding wavefronts: sub-wavelength structures [13, 14], stress birefringence [8], interferometer-based superpositions [15, 7, 16, 17, 18, 19], q-plates [20, 21, 12]. In particular the generation of monstar singularities requires the rupture of symmetry of the circular rotation rate of the polarization ellipses semi-major axis. The superpositions of two modes carrying orbital angular momentum (OAM) in orthogonal states of circular polarization breaks off the symmetry by varying the asymmetry of one the optical vortices [22]. Other methods of monstar generation have used geometric phase optical elements as elliptically-symmetric q-plates [12] and Spatially Varying Axis retardation waveplates (SVAPs) [23, 24]. Q-plates are nematic liquid-crystal cells where a disclination of order qq is patterned by the alignment direction of the liquid crystal molecules called ’director’. The disclination pattern orientation is a function of the angular coordinate ϕ\phi. Symmetry dislocations as lemons and stars have been generated with q-plates devices with circular symmetry in the orientation of the director pattern. Elliptically-symmetric q-plates and SVAPs have non-circularly symmetry disclination patterns and their director pattern follows a mm-fold rotational symmetry encoded in the form of the Pancharatnam-Berry phase. The rotational symmetry of order mm and other geometrical properties encoded into the Pancharatnam-Berry phase are able to add radial lines to the polarization cross-section structure of the light beam. If θ\theta is the orientation of polarization at a particular point, and ϕ\phi is the angular coordinate, the radial orientation is defined as the orientation of the disclination relative to the radial direction [11, 25]:

θr=θϕ.\theta_{r}=\theta-\phi\,. (1)

Radial lines occur when θr=nπ\theta_{r}=n\pi, with nn integer. Radial lines mark the boundaries of angular sectors. The patterns contained within the angular sectors have three types of lines classified as elliptic, hyperbolic, and parabolic. The radial orientation in these sectors change from line to another. Elliptic regions have lines that begin and end at the C-point and the change of radial orientation, δθr\delta\theta_{r}, is +π+\pi. Lines in hyperbolic regions avoid the C-point and δθr=π\delta\theta_{r}=-\pi. Parabolic regions appear as a characteristic in monstar patterns. In theses regions one end of the lines connects with the C-point. The number of radial lines for lemons and stars follows the equation [11, 25]:

N=|2(IC1)|.N=|2(I_{C}-1)|. (2)

Lemons have the rate of rotation of the semi-major axis constant along a circular path centered about the C-point. They have only one direction where the semi-major axis is radial: the semi-major axis rotates at half the rate of circulation about the C-point. Stars rotate in the sense opposite to the path, the axis is radial along three lines. Monstar pattern is more general than lemon and star because in these patterns the semi-major axis has not a constant rate of rotation and the axis may be radial along three or more lines. It contain parabolic sectors delineated by the separatrices, where all the lines of curvature have the C-point as an end point. Monstar patterns have a number of radial lines that do not follow Eq. 2, and contain parabolic sectors in addition to one or more sectors that can be either parabolic or elliptical. Monstars patterns with positive-, negative- and zero-indexes have been reported generated with non-separable superpositions of three spatial modes carrying orbital angular momentum (OAM) in opposite states of circular polarization [7, 26, 22]. The criterion of monstar generation has being retaining the index of lemons and stars while adding more radial lines and sectors to invalidate the Eq. 2. In this article we show that the structural features of monstar dislocations can be generated by the Spatially Varying Axis Waveplate (SVAP) device. Monstar singularities require to manipulate the rotation rate of the polarization azimuth around the C- or V -points. SVAP is a device in which the function of the disclination pattern orientation draws the background of a beam carrying OAM with topologial charge qq in superposition with a phase modulator factor Φ(ϕ)\Phi(\phi) which is the unitary normal vector orientation of a closed curve [23, 24]. This method is able to engineer and set on demand the rotation rate and the structure of the local polarization azimuth around the singularity. The azimuth orientation information is actually encoded in the form of Pancharatnam-Berry phase. SVAP device was introduced in the design and generation of the Free Form Dark Hollow (FFDHmFFDH^{m}_{\ell}) beams of topological charge \ell and symmetry order mm. It was used to the generating of exotic polarization-structured beams [23, 24]. Non-separable superpositions of FFDH in orthogonal states of polarization produce beams of light that have spatially variable polarization. They map the Poincaré sphere onto its transverse mode and carry disclinations by design and have been used to produce high-order disclinations. FFDH can easily produce symmetric lemons and stars. In this work we demostrated the production of monstars with non-separable superposition of FFDH by increasing the number of radial lines (L-lines) and angular sectors in lemon- and star-like patterns while the index IC=1/2I_{C}=-1/2 is maintained with the same value. We illustrate the operation of SVAPs encoding geometric phase distribution with a non-uniform azimuthal dependence, explain how we use them to create various polarization patterns and how to modulate its parameters in order to get monstar singularities. we introduce some examples of such beams and demonstrate their properties on experimental grounds.

I.1 Operation of SVAPs

The operation principle of SVAP devices has been already described elsewhere [27, 21, 28, 29]. In particular, the optical effect of a generic Free-Form Azimuthal SVAP can be described by the operator Q^(q,m,ξ)\hat{Q}(q,m,\xi) acting on the local polarization state of the field, defined as

Q^(q,m,ξ)\displaystyle\hat{Q}(q,m,\xi) =\displaystyle= cosδ2(|LL|+|RR|)\displaystyle\cos{\frac{\delta}{2}}\left(|L\rangle\,\langle L|+|R\rangle\,\langle R|\right) (3)
+isinδ2(ei2Ψ(ϕ;q,m,ξ)|RL|)\displaystyle+{\rm i}\sin{\frac{\delta}{2}}\left({\rm e}^{{\rm i}2\Psi(\phi;q,m,\xi)}|R\rangle\langle L|\right)
+isinδ2(ei2Ψ(ϕ;q,m,ξ)|LR|),\displaystyle+{\rm i}\sin{\frac{\delta}{2}}\left({\rm e}^{-{\rm i}2\Psi(\phi;q,m,\xi)}|L\rangle\langle R|\right),

in which the phase factor ei 2Ψ(ϕ;q,m,ξ){\rm e}^{{\rm i}\,2\,\Psi(\phi;q,m,\xi)} corresponds to Pancharatnam-Berry phase function of the parameters: topological charge 2q2q, phase modulation ϕ\phi, mm-fold symmetry of a closed curve ρ\rho, and ξ\xi the geometrical properties of the closed curve ρ\rho. All these parameters are directly mapped to the generated Free-Form Dark Hollow, FFDHmFFDH^{m}_{\ell} beam. Our device SVAP allows the electrical control of the phase retardation (birefringence) δ\delta, which can be set at δ=π\delta=\pi called half wave plate configuration HWP or at δ=π/2\delta=\pi/2 the quarter wave plate QWP configuration. Let us consider the two operation modes of SVAP device and linear and circular input polarization to obtain four polarization structures: radial-like, azimuthal-like, lemon-like, and star-like. In the first case, the SVAP retardation is set to δ=π/2\delta=\pi/2 (quarter-wave plate operation) and the input light is a circularly polarized Gaussian beam, the output beam is a coherent superposition of a Gaussian beam and a beam carrying the geometric phase with orthogonal polarization respect to the input beam. These beams are so-called Modulate Poincaré Beams due to their diversity in polarization singularity landscapes [24]. In the near field, the resulting MPB has a linear polarization map just point-by-point coincident with the SVAP optic axis distribution, with a central C-point of index IC=qI_{C}=q. A second structure of the polarization map can then be obtained by using a haf-wave plate after the SVAP to exchange the handedness of the circular components of the MPB. SVAP device reverses the circular polarization direction, introduce Orbital Angular Momentum (OAM) by 2qϕ2q\,\phi, and the modulation phase factor giving by qΦ(ϕ)q\,\Phi(\phi). In the third configuration mode, the SVAP retardation is set to δ=π\delta=\pi (half-wave plate operation) and the input polarization is linear. SVAP transforms linearly polarized light into a coherent superposition of right- and left-handed polarized light with opposite OAM and an additional mm-fold rotation-symmetric phase values, which add at every point in the beam transverse plane leading to a vector field with a central V-point with IV=2qI_{V}=2q. The MPB vector field is of course also dependent on the input linear polarization orientation, which introduces a further local rotation of the polarization by a uniform angle giving as output either radial- or azimuthal-like polarizations.

Refer to caption
Figure 1: Modulated vs unmodulated helical wavefronts: (a) uniform helical wavefront for =2\ell=2 (γ\gamma_{\infty}); (b) modulated helical wavefront for =2\ell=2 and γ5(a=b,n1=1/2;n2=n3=4/3)\gamma_{5}(a=b,n_{1}=1/2;n_{2}=n_{3}=4/3).

An FFH spatial mode of mm-order and topological charge =2q\ell=2q is generated by adding an additional mm-fold rotation-symmetric phase, ψ¯(ϕ)\bar{\psi}(\phi) to a beam carrying OAM with topological charge \ell as we can see the pictured example in Fig. 1. This additional phase gives a nonuniform distribution of the phase with an mm-fold symmetry structure. The modulation azimuthal phase factor, eiqψ¯(ϕ)e^{{\rm i}\,q\,\bar{\psi}(\phi)}  [23], corresponds to the optic axis distribution, which is giving by the orientation of the normal vector 𝐧^\mathbf{\hat{n}} along the closed curve ρ(ϕ)\rho(\phi).
The Pancharatnam–Berry phase for the SVAP device is:

ei2Ψ(ϕ;q,m,ξ)=eiϑ(ϕ),e^{{\rm i}2\Psi(\phi;q,m,\xi)}=e^{{\rm i}\ell\vartheta(\phi)}, (4)

where ϑ\vartheta is defined as the polar angle of the unit vector 𝒏^=(nx,ny)\hat{\bm{n}}=(n_{x},n_{y}) which reads:

nx(ϕ)=cosϑ=ρcosϕ+ρ˙sinϕ(ρ˙2+ρ2)1/2,\displaystyle n_{x}(\phi)=\cos\vartheta=\frac{\rho\cos\phi+\dot{\rho}\sin{\phi}}{\left(\dot{\rho}^{2}+\rho^{2}\right)^{1/2}},
ny(ϕ)=sinϑ=ρsinϕρ˙cosϕ(ρ˙2+ρ2)1/2,\displaystyle n_{y}(\phi)=\sin\vartheta=\frac{\rho\sin\phi-\dot{\rho}\cos{\phi}}{\left(\dot{\rho}^{2}+\rho^{2}\right)^{1/2}}, (5)

the closed curve and its first derivative read:

ρ(ϕ)\displaystyle\rho(\phi) =\displaystyle= (|cosmϕ4a|n2+|sinmϕ4b|n3)1n1,\displaystyle\left(\left|\frac{\cos{\frac{m\phi}{4}}}{a}\right|^{n_{2}}+\left|\frac{\sin{\frac{m\phi}{4}}}{b}\right|^{n_{3}}\right)^{-\frac{1}{n_{1}}},
ρ˙\displaystyle\dot{\rho} =\displaystyle= dρdϕ.\displaystyle\frac{d\rho}{d\phi}. (6)

The Eq. (I.1) represents the superformula in polar coordinates for multiple classes of plane curves γm(a,b,n1,n2,n3)\gamma_{m}(a,b,n_{1},n_{2},n_{3}) of the most diverse kinds [30]. Specifically ρ\rho is the distance of a point of the curve γm\gamma_{m} from the origin of the coordinate system as a function of the azimuthal angle ϕ\phi, mm is an integer number, n1n_{1}, n2n_{2} and n3n_{3} are three integers controlling its local radius of curvature and, finally, the positive real numbers aa and bb parameterize the radii of the circumferences respectively inscribed and circumscribed to the curve γm\gamma_{m}. It is worth noticing that, for m=4m=4, a=ba=b and n2=n3>2n_{2}=n_{3}>2, the superformula simply returns the superellipses first introduced by G. Lamé in 1818 [31].
The phase of an FFH mode is hence proportional to the angle ϑ(ϕ)\vartheta(\phi) formed by the unit vector 𝒏^\hat{\bm{n}} normal to the curve γm(a,b,n1,n2,n3)\gamma_{m}(a,b,n_{1},n_{2},n_{3}) with respect to a specified reference axis xx.

An \ell-charge mm-order FFH mode obviously does not carry a well-defined OAM and its spectrum includes multiple indices reflecting the mm-fold rotational symmetry of the generating curve. In order to determine such indices, it is convenient to express the phase Ψ(ϕ)\Psi(\phi) as follows

ei2Ψ(ϕ)=(nx+iny)=[ρ(ϕ)iρ˙(ϕ)ρ(ϕ)+iρ˙(ϕ)]qei 2qϕ=eiqψ¯(ϕ)ei 2qϕ,e^{{\rm i}2\Psi(\phi)}=\left(n_{x}+{\rm i}\,n_{y}\right)^{\ell}=\left[\frac{\rho(\phi)-{\rm i}\,\dot{\rho}(\phi)}{\rho(\phi)+{\rm i}\,\dot{\rho}(\phi)}\right]^{q}e^{{\rm i}\,2q\,\phi}=e^{{\rm i}\,q\,\bar{\psi}(\phi)}\,e^{{\rm i}\,2q\,\phi}, (7)

where we set q=/2q=\ell/2. qq is hence integer or half-integer depending on whether \ell is even or odd, respectively. The FFH phase factor therefore can be separated into the product of the helical background phase factor eiϕ=ei 2qϕe^{{\rm i}\ell\phi}=e^{{\rm i}\,2q\,\phi} times the modulation phase factor eiqψ¯(ϕ)e^{{\rm i}\,q\,\bar{\psi}(\phi)}. The latter is a periodic function of ϕ\phi with period 2π/m2\pi/m, i.e.

eiqψ¯(ϕ+2πm)=eiqψ¯(ϕ),e^{{\rm i}\,q\,\bar{\psi}(\phi+\frac{2\pi}{m})}=e^{{\rm i}\,q\,\bar{\psi}(\phi)}, (8)

and hence can be expanded in a Fourier series

eiqψ¯(ϕ)\displaystyle e^{{\rm i}\,q\,\bar{\psi}(\phi)} =\displaystyle= hχheihmϕ,\displaystyle\sum_{h}{\chi_{h}}e^{{\rm i}\,h\,m\,\phi},
χh\displaystyle{\chi_{h}} =\displaystyle= m2π02πm[ρ(ϕ)iρ˙(ϕ)ρ(ϕ)+iρ˙(ϕ)]qeihmϕ𝑑ϕ.\displaystyle\frac{m}{2\pi}\int_{0}^{\frac{2\pi}{m}}{\left[\frac{\rho(\phi)-{\rm i}\,\dot{\rho}(\phi)}{\rho(\phi)+{\rm i}\,\dot{\rho}(\phi)}\right]^{q}\,e^{-{\rm i}\,h\,m\,\phi}\,d\phi}. (9)

Then, the OAM spectrum of an \ell-charge mm-order FFH beam turns out to include only the components with indices ±hm\ell\pm h\,m, hh being any integer and =2q\ell=2q the OAM index corresponding to the background helical mode. The azimuthal phase factor in Eq. (10) can be finally expressed as the following helical mode expansion

ei2Ψ(ϕ)=hχhei(2q+hmϕ)=lcleilϕ,l=2q+hm.e^{{\rm i}2\Psi(\phi)}=\sum_{h}\chi_{h}e^{{\rm i}\,(2q+h\,m\,\phi)}=\sum_{l}c_{l}e^{{\rm i}\,l\,\phi},\hskip 28.45274ptl=2q+h\,m. (10)

The mean value of the OAM in an mm-order FFH mode can be easily calculated from Eqs. (4) and (I.1),

Lz\displaystyle\langle L_{z}\rangle =\displaystyle= i2π02πei2Ψ(ϕ)ϕei2Ψ(ϕ)\displaystyle-\frac{{\rm i}\hbar}{2\pi}\int_{0}^{2\pi}{e^{-{\rm i}2\Psi(\phi)}\,\frac{\partial}{\partial\phi}\,e^{{\rm i}2\Psi(\phi)}} (11)
=\displaystyle= πΨ(ϕ)|02π=πqϑ(ϕ)|02π=2q.\displaystyle\frac{\hbar}{\pi}\,\left.\Psi(\phi)\right|_{0}^{2\pi}=\frac{\hbar}{\pi}q\left.\vartheta(\phi)\right|_{0}^{2\pi}=2\hbar q.

This result indicates that the mean OAM carried by an \ell-charge FFH mode is always =2q\hbar\ell=2\hbar q per photon, independently of mm. It is therefore coincident with the OAM per photon of the background helical mode and it is proportional to the overall topological charge of the transverse wavefront. The phase modulation factor eiqψ¯(ϕ)e^{{\rm i}\,q\,\bar{\psi}(\phi)} defined in Eq. (7) has a zero mean OAM. However, the specific value of the fraction |cl|2|c_{l}|^{2} of the total power of the optical field carrying an OAM proportional to ll depends on the geometric details of the generating curve γm(a,b,n1,n2,n3)\gamma_{m}(a,b,n_{1},n_{2},n_{3}) and ultimately on the specific values of the parameters in Eq. (I.1). Uniform helical modes can be easily obtained as a special case of FFH modes when the generating curve γm\gamma_{m} degenerates into a circumference (γmγ\gamma_{m}\rightarrow\gamma_{\infty}). Free space propagation obviously alters the shape of FFH modes, but it does not influence their symmetry properties.

In this work we introduce three examples of SVAP device with parameters m=6m=6, n2=2.3n_{2}=2.3, n3=2.3n_{3}=2.3, a=ba=b and distinct parameters n1n_{1} for each distinct SVAP profile. With these examples, we illustrate the ability of the SVAP device to produce monstar singularities by adding radial lines and angular sectors to the polarization structure that originally performs as a lemon-like or as a star-like . The three profiles samples maintained the same mm-fold rotational symmetry and also the parameters a=ba=b, n2=n3n_{2}=n_{3}, only the parameter n1n_{1} was set distinct for each profile. Parameter n1n_{1} is able to modulate the radius of curvature of γm(a,b,n1,n2,n3)\gamma_{m}(a,b,n_{1},n_{2},n_{3}), while 1/n11/n_{1} increases in value, then the radius of curvature of γm(a,b,n1,n2,n3)\gamma_{m}(a,b,n_{1},n_{2},n_{3}) decreases respect to the origin of the coordinate system. It seems as if the beam were compressed towards the center. In Fig. 2 can be observed the transverse intensity of the beams generated with different parameter 1/n11/n_{1}. The beam is more compressed as the value of 1/n11/n_{1} increases and also, the number of radial lines L-lines increases as well. We show this in the section of results and discussion.

Refer to caption
Figure 2: Transverse intensity of the beams generated by the SVAP devices designed with a 6-fold symmetry, n2=n3=2.3n_{2}=n_{3}=2.3, a=ba=b and different number of L-lines modulated by the parameter n1n_{1}. a) Profile 1 has n1=2/3n_{1}=2/3. b) Profile 2 has n1=3/10n_{1}=3/10. c) Profile 3 has n1=1/6n_{1}=1/6.

II Experiment

Refer to caption
Figure 3: Experimental setup for the measurements in near field.

Our liquid crystals SVAPs have been fabricated adopting a “direct-write approach” to pattern the LC anchoring to the ITO-coated glass walls of the cell. Electrical control of the retardation δ\delta – or electrical tunability – and non-diffractive operation are just some of the main advantages of this technology. The pattern of each profile was generated by following the orientation respect to the xx-axis of the normal vector of the 6-fold symmetry curve as it was described in section I.1.

To generate the beams with spatially varying polarization we use a beam from a Helium-Neon laser centred at λ=633\lambda=633nm with a beam diameter of  0.2mm. The Fig. 3 shows the experimental setup. The beam is sent through the lenses L1 and L2 with focal lengths of f1=f_{1}=5 cm and f2=f_{2}=20 cm and separated by 25 cm forming a telescope T1 to magnifying the beam by a factor of  4×\times. A linear polarizer (P1) was placed to orientate the beam into linear polarization. We used a quarter wave plate (QWP1) or a half wave plate (HWP1) to obtain either circular or linear input polarization. The polarized beam impinges on the SVAP device. SVAP turns the Gaussian beam into a structured beam carrying the Pancharatnam-Berry phase. To generate the spatially-variable polarization singularities, the retardation δ\delta must be put in half wave configuration δ=π\delta=\pi (HWP) or quarter wave plate configuration δ=π/2\delta=\pi/2 (QWP) to get the coherent superposition of two modes with perpendicular polarization. Our SVAP devices produce spatially varying polarization singular patterns shaped by the symmetry properties of the curve encoded in the Pancharatnam-Berry phase. We measured the transverse polarization distribution in the near field. The light at the plane of SVAP (near field) was re-imaged onto a second plane by two lenses, L3L_{3} and L4L_{4}, with focal length f=15f=15cm in 4-f configuration. The plane of near field was located at the focal distance of L4, in which a CCD camera was placed. In order to obtain the state of polarization of the beams, we measured the Stokes parameters with the half wave plate (HWP4) and with the QWP2 which were placed before of CCD plane and the analyser (P2).

III Results and Discussion

In this section we present the transverse structures of polarization that were obtained by using the three designed SVAP devices. We have called these patterns as radial-like, azimuthal-like, lemon-like, and star-like due to they resemble the basic singularities patterns with the additional m-symmetry phase factor which draws the shape of symmetry on the pattern. We show how to generate monstar singularities by breaking the circular symmetry of the rotation rate of the axis of orientation of the polarization of lemon- and the star-like patterns so that we obtain monstar with positive and negative index ICI_{C}. The intensity distributions showed in Fig. 2 corresponds to the far field, HWP operation mode δ=π\delta=\pi, and circular input polarization. The parameter n1n_{1} is associated to the number of radial lines L-lines and angular sectors present in the polarization structure. It appears as the fractional power 1n1\frac{1}{n_{1}} in the expression for ρ(ϕ)\rho(\phi) shown in Eq. I.1 having numerical values for our designed SVAPs; 32\frac{3}{2}, 103\frac{10}{3}, and 6, for profile 1, profile 2, and profile 3 respectively. We can see that by increasing the parameter 1n1\frac{1}{n_{1}} the structures of the transverse intensity get narrower, in particular, the radius of curvature measured with respect to the phase singularity is reduced as if the beam were compressed towards the center shaping narrower edges. We obtained the polarization structure in near field for each profile and showed the evolution of the basic-like structures turning into monstar singularities when L-lines and angular sectors are added by increasing the parameter 1n1\frac{1}{n_{1}} while the value of index IC=|1/2|I_{C}=|1/2| is maintained.
We show theory and measurements for the four configurations of input polarization and retardance value δ\delta. Azimuthal- and radial-like patterns are performed with retardance δ=π\delta=\,\pi and horizontal and vertical input polarization respectively. Lemon- and star-like configurations are performed with retardance δ=π/2\delta=\pi/2 and right circular input polarization to obtain star-like pattern. A half wave plate is placed after the SVAP, showed as HWP3 in Fig. 3 to flip the handedness of circular polarization and obtain lemon-like pattern. Figures 4, 5 and 6 show polarization structure in near field obtained with SVAP of: 1n1=32\frac{1}{n_{1}}=\frac{3}{2} (profile 1), 1n1=103\frac{1}{n_{1}}=\frac{10}{3} (profile 2), and 1n1=6\frac{1}{n_{1}}=6 (profile 3) respectively, the first row shows the theoretical patterns and the second row shows the experimental measurements. The patterns appear in the following order: azimuthal-like, radial-like, star-like, and lemon-like. We can observe the evolution of the angle of polarization orientation in the polarizations structures for the three samples with different value of parameter 1/n11/n_{1}. Modulation phase factor ψ¯(ϕ)\bar{\psi}(\phi) gives structure to the transverse intensity and polarization. In the azimuthal-like configurations in figures 4a), 5a), 6a) we observe a structure like flower petals inherited from the 6 order rotational symmetry of the curve ρ\rho, and that is stronger as parameter 1n1\frac{1}{n_{1}} increases. In particular, the star-like configurations show more radial lines and angular sectors as the parameter 1/n11/n_{1} increases turning the usually basic star structure of profile 1 into a more complex structures for profiles 2 and 3 which have six and nine radial lines and angular sectors respectively. These last two are monstar patterns.

Radial lines are indicated with a yellow background colour in figures 7 and 8. It shows the patterns star-like and lemon-like of the profiles 1, 2, and 3 with a background colour corresponding to the orientation of polarization relative to the radial direction θr=θϕ\theta_{r}=\theta-\phi. Radial lines are indicated in yellow colour and angles are in radians. Monstars distinguishes from lemons and stars because their parabolic sectors, delimited by radial lines where the lines radiate from the singularity. A disclination pattern in which the number of radial lines does not follow the equation 2 is defined as a monstar disclination [4]. We can see the transition from the star- and lemon-like patterns to monstar patterns by adding sectors and radial lines by the modulation of the parameter 1n1\frac{1}{n_{1}} in each sample of SVAP device. In the star-like Figure 7b) is shown the star-like pattern of profile 1 which has 1n1=32\frac{1}{n_{1}}=\frac{3}{2}, all sectors are hyperbolic and it has three radial lines as the classical star pattern. In Fig. 7d) is shown the star-like pattern of profile 2 which has 1n1=103\frac{1}{n_{1}}=\frac{10}{3}, there are six radial lines, it seems as if each radial line of the basic star pattern were split into two. Within these radial lines, there are parabolic sectors that change their radial orientation from one line to another by δrπ/8\delta_{r}\leq\pi/8. Finally, Fig. 7f) shows the star-like pattern of profile 3 which has 1n1=6\frac{1}{n_{1}}=6, in this pattern the parabolic sectors and radial lines split into three and the radial orientation change form one to another by δrπ/8\delta_{r}\leq\pi/8. The structure has nine lines and parabolic sectors from one line to another. Figures 7d) and 7f) are monstar because the number of radial lines does not follow the equation 2. These monstars patterns have an index IC=1/2I_{C}=-1/2.

The generation of monstar patterns with the SVAP device allows adding sectors and radial lines by the modulation of the Pancharatnam-Berry phase while maintaining the value of orbital angular momentum (OAM) which is directly related to the index ICI_{C}. The transition from lemon patterns to monstar patterns is illustrated in Fig. 8. Radial lines pointed out with the yellow colour background indicating that relative angle of the polarization orientation and the radial direction is zero. False colour encodes the orientation of the polarization relative to the radial direction θr\theta_{r}. A remarkable characteristic about the monstars in Figures 7 and 8 is that all the lines of curvature in the sectors delineated by the separatrices, have the C-point as an end point. In lemon-like patterns we can observe how the radial lines and parabolic sectors increase in the patterns of profiles 2 and 3. Profile 2 has for example, two parabolic regions localized at the angular positions: 3636^{\circ} to 5454^{\circ}, and 26-26^{\circ} to 49-49^{\circ}, and the usual radial line for lemons localized at zero rads. Profile 3 also has two paraboloc regions localized at the angular positions: 3636^{\circ} to 4949^{\circ}, and 36-36^{\circ} to 57-57^{\circ} and the usual radial line for lemons localized at zero rads. Patterns shown in Figures 8d) and 8f) are monstar because the number of radial lines does not follow the equation 2. These monstars patterns have an index IC=1/2I_{C}=1/2. The images have a resolution of 60 ×\times 60 pixels in size.

Refer to caption
Figure 4: Theoretical and experimental measurements of the polarisation structure in near field for the profile 1. First row shows the theoretical results: a) azimuthal-like pattern, b) radial-like pattern, c) star-like pattern, and d) lemon-like pattern. Second row shows the experimental measurements: e) azimuthal-like pattern, f) radial-like pattern, g) star-like pattern, and h) lemon-like pattern.
Refer to caption
Figure 5: Theoretical and experimental measurements of the polarisation structure in near field for the profile 2. First row shows the theoretical results: a) azimuthal-like pattern, b) radial-like pattern, c) star-like pattern, and d) lemon-like pattern. Second row shows the experimental measurements: e) azimuthal-like pattern, f) radial-like pattern, g) star-like pattern, and h) lemon-like pattern.
Refer to caption
Figure 6: Theoretical and experimental measurements of the polarisation structure in near field for the profile 3. First row shows the theoretical results: a) azimuthal-like pattern, b) radial-like pattern, c) star-like pattern, and d) lemon-like pattern. Second row shows the experimental measurements: e) azimuthal-like pattern, f) radial-like pattern, g) star-like pattern, and h) lemon-like pattern.
Refer to caption
Figure 7: Numerical simulations and experimental measurements of the star patterns. False color encodes the orientation of the polarization relative to the radial direction: θr=θϕ\theta_{r}=\theta-\phi. L-lines are along the yellow colour.
Refer to caption
Figure 8: Numerical simulations and experimental measurements of the star patterns. False color encodes the orientation of the polarization relative to the radial direction: θr=θϕ\theta_{r}=\theta-\phi. L-lines are along the yellow colour

IV Conclusions

We presented the generation of beams with spatially varying polarization emphasizing monstars patterns through the non- separable coherent superpositions of FFDH beams in orthogonal states of polarization produced by Spatially Varying Axis waveplates SAVP devices that were designed with parameters of geometric properties of a closed curve. We showed how the star- and lemon-like patterns evolve from the basic structures and number of radial lines and angular sectors to more complex structures, and radial lines, and angular sectors which are monstar patterns. It is remarkable that these monstars patterns were generated without changing the index ICI_{C} to add angular sectors to the basic structures and the rupture of symmetry was done on the rate rotation of the optic axis of the devices which means also a non-symmetric Panchratnam-Berry phase. The rupture of symmetry can also be seen in the rotation rate of the local polarization azimuth around the singularity. We showed theoretical and experimental measurements with excellent agreement.

Acknowledgements.
This work was supported by the European Union (EU) within Horizon 2020—European Research Council Advanced Grant No. 694683 (PHOSPhOR) and by INFN – Naples, ETHIOPIA (CSN 5).

References