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
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, , of topological charge , and order of symmetry . 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.
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 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 , indicating that the azimuthal coordinate on the Poincaré sphere rotates times per turn about the singularity. Basic polarization disclinations are classified as lemons (, 1 radial line), stars (, 3 radial lines) and monstars (, 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 and stars with 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 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 is patterned by the alignment direction of the liquid crystal molecules called ’director’. The disclination pattern orientation is a function of the angular coordinate . 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 -fold rotational symmetry encoded in the form of the Pancharatnam-Berry phase. The rotational symmetry of order 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 is the orientation of polarization at a particular point, and is the angular coordinate, the radial orientation is defined as the orientation of the disclination relative to the radial direction [11, 25]:
(1) |
Radial lines occur when , with 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, , is . Lines in hyperbolic regions avoid the C-point and . 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]:
(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 in superposition with a phase modulator factor 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 () beams of topological charge and symmetry order . 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 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 acting on the local polarization state of the field, defined as
(3) | |||||
in which the phase factor corresponds to Pancharatnam-Berry phase function of the parameters: topological charge , phase modulation , -fold symmetry of a closed curve , and the geometrical properties of the closed curve . All these parameters are directly mapped to the generated Free-Form Dark Hollow, beam. Our device SVAP allows the electrical control of the phase retardation (birefringence) , which can be set at called half wave plate configuration HWP or at 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 (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 . 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 , and the modulation phase factor giving by . In the third configuration mode, the SVAP retardation is set to (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 -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 . 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.

An FFH spatial mode of -order and topological charge is generated by adding an additional -fold rotation-symmetric phase, to a beam carrying OAM with topological charge as we can see the pictured example in Fig. 1. This additional phase gives a nonuniform distribution of the phase with an -fold symmetry structure. The modulation azimuthal phase factor, [23], corresponds to the optic axis distribution, which is giving by the orientation of the normal vector along the closed curve .
The Pancharatnam–Berry phase for the SVAP device is:
(4) |
where is defined as the polar angle of the unit vector which reads:
(5) |
the closed curve and its first derivative read:
(6) |
The Eq. (I.1) represents the superformula in polar coordinates for multiple classes of plane curves of the most diverse kinds [30]. Specifically is the distance of a point of the curve from the origin of the coordinate system as a function of the azimuthal angle , is an integer number, , and are three integers controlling its local radius of curvature and, finally, the positive real numbers and parameterize the radii of the circumferences respectively inscribed and circumscribed to the curve . It is worth noticing that, for , and , 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 formed by the unit vector normal to the curve with respect to a specified reference axis .
An -charge -order FFH mode obviously does not carry a well-defined OAM and its spectrum includes multiple indices reflecting the -fold rotational symmetry of the generating curve. In order to determine such indices, it is convenient to express the phase as follows
(7) |
where we set . is hence integer or half-integer depending on whether is even or odd, respectively. The FFH phase factor therefore can be separated into the product of the helical background phase factor times the modulation phase factor . The latter is a periodic function of with period , i.e.
(8) |
and hence can be expanded in a Fourier series
(9) |
Then, the OAM spectrum of an -charge -order FFH beam turns out to include only the components with indices , being any integer and 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
(10) |
The mean value of the OAM in an -order FFH mode can be easily calculated from Eqs. (4) and (I.1),
(11) | |||||
This result indicates that the mean OAM carried by an -charge FFH mode is always per photon, independently of . 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 defined in Eq. (7) has a zero mean OAM. However, the specific value of the fraction of the total power of the optical field carrying an OAM proportional to depends on the geometric details of the generating curve 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 degenerates into a circumference (). 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 , , , and distinct parameters 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 -fold rotational symmetry and also the parameters , , only the parameter was set distinct for each profile. Parameter is able to modulate the radius of curvature of , while increases in value, then the radius of curvature of 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 . The beam is more compressed as the value of increases and also, the number of radial lines L-lines increases as well. We show this in the section of results and discussion.

II Experiment

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 – 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 -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 nm 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 5 cm and 20 cm and separated by 25 cm forming a telescope T1 to magnifying the beam by a factor of 4. 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 must be put in half wave configuration (HWP) or quarter wave plate configuration (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, and , with focal length cm 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 . The intensity distributions showed in Fig. 2 corresponds to the far field, HWP operation mode , and circular input polarization. The parameter is associated to the number of radial lines L-lines and angular sectors present in the polarization structure. It appears as the fractional power in the expression for shown in Eq. I.1 having numerical values for our designed SVAPs; , , and 6, for profile 1, profile 2, and profile 3 respectively. We can see that by increasing the parameter 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 while the value of index is maintained.
We show theory and measurements for the four configurations of input polarization and retardance value . Azimuthal- and radial-like patterns are performed with retardance and horizontal and vertical input polarization respectively. Lemon- and star-like configurations are performed with retardance 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: (profile 1), (profile 2), and (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 . Modulation phase factor 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 , and that is stronger as parameter increases. In particular, the star-like configurations show more radial lines and angular sectors as the parameter 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 . 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 in each sample of SVAP device. In the star-like Figure 7b) is shown the star-like pattern of profile 1 which has , 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 , 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 . Finally, Fig. 7f) shows the star-like pattern of profile 3 which has , in this pattern the parabolic sectors and radial lines split into three and the radial orientation change form one to another by . 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 .
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 . 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 . 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: to , and to , and the usual radial line for lemons localized at zero rads. Profile 3 also has two paraboloc regions localized at the angular positions: to , and to 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 . The images have a resolution of 60 60 pixels in size.





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 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
- Forbes et al. [2021] A. Forbes, M. de Oliveira, and M. R. Dennis, Structured light, Nature Photonics 15, 253 (2021).
- M.V and M.R [2001] B. M.V and D. M.R, Polarization singularities in isotropic random vector waves, Proc. R. Soc. Lond. A 457, 141?155 (2001).
- Gbur [2015] G. Gbur, Singular Optics. The Optics Encyclopedia (Wiley-VCH Verlag Gmbh & Co. KGaA, Weinheim, 2015).
- Galvez et al. [2014] E. J. Galvez, B. L. Rojec, V. Kumar, and N. K. Viswanathan, Generation of isolated asymmetric umbilics in light’s polarization, Phys. Rev. A 89, 031801 (2014).
- Torres and Torner [2011] J. Torres and L. Torner, Twisted Photons: Applications of Light with Orbital Angular Momentum (Wiley, 2011).
- Cardano and Marrucci [2015] F. Cardano and L. Marrucci, Spin-orbit photonics, Nature Photonics 9, 776 (2015).
- Galvez et al. [2012] 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 (2012).
- Beckley et al. [2010] A. M. Beckley, T. G. Brown, and M. A. Alonso, Full poincaré beams, Opt. Express 18, 10777 (2010).
- Senthilkumaran [2018] P. Senthilkumaran, Singularities in Physics and Engineering, 2053-2563 (IOP Publishing, 2018).
- Freund [2011] I. Freund, Mbius strips and twisted ribbons in intersecting gauss-laguerre beams, Optics Communications 284, 3816 (2011).
- Khajavi and Galvez [2016a] B. Khajavi and E. J. Galvez, High-order disclinations in space-variant polarization, Journal of Optics 18, 084003 (2016a).
- Cvarch et al. [2017] B. A. Cvarch, B. Khajavi, J. A. Jones, B. Piccirillo, L. Marrucci, and E. J. Galvez, Monstar polarization singularities with elliptically-symmetric q-plates, Opt. Express 25, 14935 (2017).
- Niv et al. [2005] A. Niv, G. Biener, V. Kleiner, and E. Hasman, Rotating vectorial vortices produced by space-variant subwavelength gratings, Opt. Lett. 30, 2933 (2005).
- Beresna et al. [2011] M. Beresna, M. Gecevičius, P. G. Kazansky, and T. Gertus, Radially polarized optical vortex converter created by femtosecond laser nanostructuring of glass, Applied Physics Letters 98, 201101 (2011), https://doi.org/10.1063/1.3590716 .
- Maurer et al. [2007] C. Maurer, A. Jesacher, S. Frhapter, S. Bernet, and M. Ritsch-Marte, Tailoring of arbitrary optical vector beams, New Journal of Physics 9, 78 (2007).
- Kumar and Viswanathan [2013] V. Kumar and N. K. Viswanathan, Topological structures in the poynting vector field: an experimental realization, Opt. Lett. 38, 3886 (2013).
- Vyas et al. [2013] S. Vyas, Y. Kozawa, and S. Sato, Polarization singularities in superposition of vector beams, Opt. Express 21, 8972 (2013).
- Moreno et al. [2014] I. Moreno, J. A. Davis, D. M. Cottrell, and R. Donoso, Encoding high-order cylindrically polarized light beams, Appl. Opt. 53, 5493 (2014).
- Otte et al. [2016] E. Otte, C. Alpmann, and C. Denz, Higher-order polarization singularitites in tailored vector beams, Journal of Optics 18, 074012 (2016).
- Cardano et al. [2013] F. Cardano, E. Karimi, L. Marrucci, C. de Lisio, and E. Santamato, Generation and dynamics of optical beams with polarization singularities, Opt. Express 21, 8815 (2013).
- Piccirillo et al. [2013] B. Piccirillo, S. Slussarenko, L. Marrucci, and E. Santamato, The orbital angular momentum of light: Genesis and evolution of the concept and of the associated photonic technology, La Rivista del Nuovo Cimento 36, 501 (2013).
- Khajavi and Galvez [2016b] B. Khajavi and E. J. Galvez, High-order disclinations in space-variant polarization, Journal of Optics 18, 084003 (2016b).
- Piccirillo et al. [2020] B. Piccirillo, E. Piedipalumbo, and E. Santamato, Geometric-phase waveplates for free-form dark hollow beams, Frontiers in Physics 8, 94 (2020).
- Darvehi et al. [2021] P. Darvehi, V. Vicuña-Hernández, L. Marrucci, E. Piedipalumbo, E. Santamato, and B. Piccirillo, Increasing the topological diversity of light with modulated poincaré beams, Journal of Optics 23, 054002 (2021).
- Galvez and Khajavi [2017] E. J. Galvez and B. Khajavi, Monstar disclinations in the polarization of singular optical beams, J. Opt. Soc. Am. A 34, 568 (2017).
- Galvez et al. [2013] E. J. Galvez, B. L. Rojec, and X. Cheng, Polarization singularities in poincare optical beams, in The Rochester Conferences on Coherence and Quantum Optics and the Quantum Information and Measurement meeting (Optical Society of America, 2013) p. M2A.3.
- Piccirillo et al. [2010] B. Piccirillo, V. D’Ambrosio, S. Slussarenko, L. Marrucci, and E. Santamato, Photon spin-to-orbital angular momentum conversion via an electrically tunable q-plate, Applied Physics Letters 97, 241104 (2010), https://doi.org/10.1063/1.3527083 .
- Alemán-Castaneda et al. [2019] L. A. Alemán-Castaneda, B. Piccirillo, E. Santamato, L. Marrucci, and M. A. Alonso, Shearing interferometry via geometric phase, Optica 6, 396 (2019).
- Rubano et al. [2019] A. Rubano, F. Cardano, B. Piccirillo, and L. Marrucci, Q-plate technology: a progress review (invited), J. Opt. Soc. Am. B 36, D70 (2019).
- Gielis [2003] J. Gielis, A generic geometric transformation that unifies a wide range of natural and abstract shapes, American Journal of Botany 90, 333 (2003).
- Sokolov [2001] D. Sokolov, Encyclopedia of mathematics, in Lamé curves, Encyclopedia of Mathematics, edited by M. Hazewinkel (Springer Science+Business Media B.V. / Kluwer Academic Publishers, 2001).