Abstract
Focused higher-order Poincaré (HOP) beams are of particular interest because they facilitate understanding the exotic properties of structured light and their applications in classical physics and quantum information. However, generating focused HOP beams using metasurfaces is challenging. In this study, we proposed a metasurface design comprising two sets of metal nanoslits for generating coaxially focused HOP beams. The nanoslits were interleaved on equispaced alternating rings. The initial rings started at the two adjacent Fresnel zones to provide opposite propagation phases for overall elimination of the co-polarization components. With the designed hyperbolic and helical profiles of the geometric phases, the two vortices of the opposite cross-circular-polarizations were formed and selectively focused, realizing HOP beams of improved quality. Simulations and experimental results demonstrated the feasibility of the proposed metasurface design. This study is of significance in the integration of miniaturized optical devices and enriches the application areas of metasurfaces.
© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
Vector beams (VBs) have garnered significant interest in recent decades owing to their unique characteristics such as spatially variant polarization [1]. They are applied extensively in fields, such as optical trapping [2], high-resolution imaging [3], ultrasensitive metrology [4], and particle acceleration [5]. Considering physical essentials of VBs, spin angular momentum (SAM) and orbital angular momentum (OAM) are associated with circular polarizations (CPs) and optical vortices, respectively, and they can be coupled into the total angular momentum (TAM) states of vector vortices [6]. Two conjugate TAM states constitute higher-order Poincaré (HOP) beams [7] and full-Poincaré beams [8], which are a fundamental category of VBs. The non-separability of these beams in space and polarization degrees of freedom is similar to the local entanglement in a bipartite system [9] that enables a range of methods and applications in quantum systems [10], such as superdense teleportations [11], asymmetric optical quantum networks [12], and quantum key distribution [13]. Accordingly, VBs provide a novel resource of protocols for encoding rotationally invariant qubits in alignment-free communication over a particular distance [14–16]. The focus of most studies is on the generation of VBs for which various methods have been developed using different interferometers, such as Michelson, Mach–Zehnder, and Sagnac [17–19] with optical elements, including Dammann gratings [20], q-plates [17], spatial light modulators [21], and spiral-phase plates [22]. However, cascaded bulky devices limit the flexibility of VB manipulation and are not suitable for integration into compact optical systems.
Metasurfaces are composed of metallic or dielectric subwavelength meta-atoms with variable geometric parameters (size, shape, and orientation) and provide a powerful and flexible platform for manipulating the light field by modifying multiple degrees of freedom, including amplitudes, phases, and polarizations. These metasurfaces have been widely applied in miniaturized and multifunctional optical devices and systems, such as broadband achromatic metalenses [23,24], OAM multiplexed holograms [25–27], VB generators [28–31], and polarization beam splitters [32,33]. They have also been used to implement multichannel transforms and develop elements for entangled photons [34,35]. Metasurfaces are powerful tools for manipulating VBs in subwavelength dimension owing to the strong spin–orbit interactions of light with the constituent meta-atoms. Particularly, VBs in frequency spectrum of microwaves [36–38], terahertz waves [39,40] and light-waves have been generated using metasurfaces with meta-atoms of sophisticated multilayer structures [36–38] and anisotropic resonators [41]. In these categories, metasurfaces for generating light-wave VBs, with the advantages of miniaturization and integration over the conventional optical systems, have attracted extensive interest and have been applied in wide range of applications such as quantum entanglement [42,43], on-chip manipulation of VBs [44], vector fiber laser [45], and full-color vectorized hologram [46]. Furthermore, in the entities of light waves, metasurfaces have been used to generate versatile VBs with radial and azimuthal polarizations [47–49], and complicated HOP [50,51] and hybrid HOP beams [31,52]. Initially, VB generation using metasurfaces were realized by directly controlling the polarization and phase of the output wave of the meta-atom [48,53–55]. Subsequently, Yue et al. [29] introduced the method of superposing two optical vortices of opposite CPs [7] into the metasurface design, and HOP beams in multiple channels were generated [56]. Devlin et al. [52] manipulated arbitrary HOP beams using the spin–orbit conversion of meta-atoms to create two circularly polarized vortices. The superposition of vortices of cross-circular polarizations (CPCs) has been used to generate enriched vector fields, including multichannel and perfect HOP beams [31,57–60]. Recently, Wang et al. proposed the solid Poincaré sphere and designed the metasurface to manipulate degree of polarization of light waves [61]. However, the VBs generated by these methods are unfocused and usually have comparatively large beam sizes; even for perfect VBs, the radii of the doughnuts are as large as several tens of micrometers, which hinders certain specific but important applications.
With the beam size reduced to the subwavelength scale [62], tightly focused VBs exhibit fascinating properties and unusual abilities, including strong manipulations [63] and super-resolutions [3]. Furthermore, they have also contributed to the discovery of exotic optical phenomena in focal regions, such as topological knots [64], Möbius strips [65], and lasing emission chirality [66]. In recent years, the generation of focused VBs with metasurfaces has attracted considerable research attention; however, most studies have been implemented on simple VBs with radial and azimuthal polarizations using either simulation designs [57,62] or experimental realizations [39,40,66]. In addition, multichannel-focused un-coaxial HOP beams have been investigated [67,68] in which the cores of multiple VBs deviating from the optical axis. However, few studies have focused on the generation of coaxially focused HOP beams [50]. The difficulties encountered in this research field are due to the fact that the dark cores of the VBs would arrive at the optical axis and coincide with the probable maximum spot of the remaining co-polarization components owing to the focusing effect of the metasurface, thus resulting in deteriorated beam quality. Even in common metadevice performances, the remaining co-polarization components may also cause opposite effects. Essentially, the components are originated from such factors as errors of the meta-atoms deviating from the nano-half-waveplate in metasurface fabrication, misalignments in optical setup and mismatches of illuminating wavelengths. To decrease the influences of the co-polarization components, different schemes have been developed in recent years. Firstly, the components can be filtered out using the polarization filter composed of quarter-wave plate and polarizer [69,70]; however, this method is inapplicable in VB generation because two CPs in illuminating light are used simultaniously [56]. Then, components are also eliminated by using the nano-polarizers as the meta-atoms of the metasurfaces based on Malus law [71,72]. Besides, they can be eliminated by the destructive interference of the light waves transmitting meta-atom-pairs, which was demonstrated in the focused VB generations [50]. Yet, the last two schemes require two meta-atoms to constitute atom-pairs, leading to inconvenience for metasurface designs and fabrications.
In this study, we propose a metal metasurface comprising two sets of nanoslits with identical dimensions to generate HOP beams by overall elimination of co-polarization components based on the principle of Fresnel zones. The initial rings of the two sets of nanoslits start at two adjacent Fresnel zones [41,50,73] to provide opposite propagation phases for their wavefields. By alternately interleaving the subsequent equispaced rings and optimizing the number of rings, the overall superposition of the light field of the two nanoslits effectively eliminate the co-polarization components. The orientation angles of the nanoslits are rotated in opposite directions along the radius and identically in the azimuthal direction to form the overall hyperbolic and helical profiles of the geometric phases for the corresponding CPs; subsequently, two vortices of the opposite CPCs can be generated and selectively focused. The initial phases and amplitudes of the two CPC vortices are controlled by the initial orientation angles of the two-set nanoslits and elliptical polarization states of the incident light. Consequently, the HOP beams are realized at different points on the HOP sphere. In the practical performace, the finite-difference time-domain (FDTD) method was used to simulate HOP beams of different orders and optimize the metasurfaces. Experimentally, the samples producing different HOP beams were fabricated, and the measurements of the generated beams were conducted. The results demonstrated the feasibility of the proposed method for generating focused HOP beams of satisfactory quality. This study is of significance for compact metasurface devices and integrated optical systems. It has widespread applications in both classical and quantum physics.
2. Principles of metasurface design
2.1 Outline of the design
Figure 1 shows a schematic of the design of a gold-film metasurface for the generation of focused coaxial HOP beams. The metasurface comprises two sets of nanoslits interleaved on alternate rings with equal increments in radius. In Fig. 1(a), the metasurface on plane OXY is illuminated by an incident light of elliptical polarization |Ein > . It is the linear combination of left-handed circular polarization (LCP) |L > and right-handed circular polarization (RCP) |R > with weights ε1 and ε2, respectively; the focused VB is observed on focal plane oxy with focal length f behind the object plane. The two sets of nanoslits in the metasurface are specified as sets A and B, respectively, and the radii of the two adjacent rings in each set have the same increment δR as demonstrated in Fig. 1(b). A nanoslit of the metasurface with length L and width W is shown in green-squared inset in Fig. 1(a). It acts as a linear polarizer [50,70] with a transmitted light field perpendicular to the longer side and converts the illuminating light of the CPs into the co-polarization and cross-polarization components. With φ’ denoting the angle of the longer side of the slit with respect to x-axis, the slit orientation angles are defined as φ = φ’ + π/2. They are designed to form the geometrical phase profiles of the helical and hyperbolic phase profiles ϕh (θ) and ϕf (R), which vary with the azimuth θ and radius R to realize vortices and focusing, respectively. Furthermore, as shown in Fig. 1(a), under illuminating RCP and LCP lights, the two helical wavefronts of the corresponding cross- LCP and RCP are focused for the two sets of slits, respectively. To reduce the total co-polarization components in the paraxial area by overall destructive interference, the annulus area of nanoslit set B is increased by a Fresnel zone in the innermost and outermost radii with respect to set A as demonstrated by the schematics and bottom-line slits in Fig. 1(b). The color backgrounds depict the Fresnel zones. By controlling the elliptical polarization states of the incident light in match with metasurface design, the VBs corresponding to a point with certain latitude and longitude on the HOP sphere are obtained as shown in purple-squared inset in Fig. 1(a).
Suppose that the metasurface contains nanoslits of 2N rings and that each set A and B occupies N rings. The first rings of sets A and B are located at the Fresnel zones of first even and odd numbers n1o and n1e = n1o + 1, and the last rings at the last zones nlo and nle = n1o + 1, respectively. Thus, the areas of sets A and B are overlapped except for the innermost and outermost zones. The metasurface covers Nfz = nle- n1o Fresnel zones altogether, while each set covers Nfz-1 zones. For a nanoslit sl at point (R, θ) on the metasurface, its orientation angle is set as φ (R, θ) = φh (θ) + φf (R), and the geometric phase for the CPC is ϕ (R, θ) = 2σφ (R, θ) = ϕh (θ) + ϕf (R). Meanwhile, ϕh (θ) and ϕf (R) are given as
2.2 Transmitted field of nanoslits
As depicted in Fig. 1(c), the unit vectors parallel and perpendicular to the longer side of the nanoslit sla are represented by $\hat{e}_{ra}^s = \textrm{cos}\varphi {^{\prime}_a}{\hat{{\boldsymbol e}}_x} + \textrm{sin}\varphi {^{\prime}_a}{\hat{{\boldsymbol e}}_y}$ and $\hat{e}_{\varphi a}^s ={-} \textrm{sin}\varphi {^{\prime}_a}{\hat{{\boldsymbol e}}_x} + \textrm{cos}\varphi {^{\prime}_a}{\hat{{\boldsymbol e}}_y}$, respectively. Here, the subscript a specifies nanoslit in set A; and ${\hat{{\boldsymbol e}}_x}$ and ${\hat{{\boldsymbol e}}_y}$ are unit vectors in the x and y directions, respectively. In addition, $\varphi {^{\prime}_a}{\boldsymbol = }{\varphi _a} - \pi /2$ is the angle of the longer side of sla with respect to x-axis as defined earlier for φ. In polar coordinates, the radial and azimuthal unit vectors are ${\hat{{\boldsymbol e}}_r} = \textrm{cos}{\theta _a}{\hat{{\boldsymbol e}}_x} + \textrm{sin}{\theta _a}{\hat{{\boldsymbol e}}_y}$ and ${\hat{{\boldsymbol e}}_\theta } ={-} \textrm{sin}{\theta _a}{\hat{{\boldsymbol e}}_x} + \textrm{cos}{\theta _a}{\hat{{\boldsymbol e}}_y}$, respectively. Then, the unit vectors $\hat{e}_{ra}^s$ and $\hat{e}_{\varphi a}^s$ are expressed as follows:
When the illuminating light is circularly polarized with chirality σ, it is represented as
From Eqs. (3) and (5), we obtain ${\hat{{\boldsymbol e}}_r} \cdot \hat{\boldsymbol e}_\varphi ^s = \textrm{cos(}{\varphi _a}\textrm{ - }{\theta _a}\textrm{)}$ and ${\hat{{\boldsymbol e}}_\theta } \cdot \hat{\boldsymbol e}_\varphi ^s = \textrm{sin(}{\varphi _a}\textrm{ - }{\theta _a}\textrm{)}$, and it follows that ${\boldsymbol E}_{in}^\sigma \cdot \hat{\boldsymbol e}_\varphi ^s = {e^{i\sigma {\varphi _a}}}/\sqrt 2$. Substituting Eqs. (3) and (4) in (5), the transmitted field of sla is expressed as
2.3 Destructive and constructive interferences of co-polarization components and CPCs under CP illuminations
When the illuminating light is circularly polarized with chirality σ, the light field U (r, α; σ) generated by the metasurface at point q (r, α) on the observation plane are expressed as ${{\boldsymbol U}_{a\textrm{, }b}}\textrm{(}r\textrm{,}\alpha \textrm{;}\sigma \textrm{) = }{\boldsymbol U}_{a\textrm{, }b}^\sigma \textrm{(}r\textrm{,}\alpha \textrm{) + }{\boldsymbol U}_{a\textrm{, }b}^{ - \sigma }\textrm{(}r\textrm{,}\alpha \textrm{)}$ with subscripts a and b referring to sets A and B, and the superscripts σ and -σ represent the co-polarization component and the CPC, respectively. The light field U (r, α; σ) produced by the metasurface is expressed as
We primarily analyzed the light field Ua (r, α; σ) of set A. Based on the Huygens–Fresnel principle for surface plasmon polariton [74] and considering that the element area occupied by the slit sla at point p (Ra, θ) is dA = RadRdθ = RaδRdθ, Ua (r, α; σ) in Eq. (7) can be written as:
The co-polarization component ${\boldsymbol U}_a^\sigma \textrm{(}r\textrm{,}\alpha \textrm{)}$ in Ua (r, α; σ) as given in Eq. (11) is expressed as follows:
To elucidate the elimination of the co-polarization component by destructive interference, we divide the summation in Eq. (12) into two layers. The inner layer is implemented on the rings included in two adjacent Fresnel zones as a group, and the outer layer is conducted on groups of the adjacent zones of set A. As the nanoslits of set A start at zone n1o and cover Nfz-1 zones, we assumed that the number of nanoslit rings of set A in two adjacent Fresnel zones of 2 M + 1 and 2 M + 2 (M = 0, 1, ……., (Nfz - 1)/2) is 2NMf for convenience, with NMf being an odd number. Then, the radius of the ja-th ring in set A is ${R_{ja}} = {R_M} + \textrm{(}ja - {M_f}\textrm{)}{\delta _R}$, where RM is the radius of the first slit ring in the (2 M + 1)-th zone with its ring number in set A denoted by Mf, and ja = Mf, …, Mf + 2NMf - 1 within the two adjacent zones. Equation (12) is then written as follows:
The analysis was similar for the co-polarization component ${\boldsymbol U}_b^\sigma \textrm{(}r\textrm{,}\alpha \textrm{)}$ of light field Ub (r, α) of set B. The nanoslits started at n1o + 1 zone and also covered Nfz-1 zones. In the adjacent two Fresnel zones, 2NMf rings are also included, and the radius of each ring is ${R_{jb}} = {R_M} + \textrm{(}jb - {M_f} + {N_{Mf}} + 1/2\textrm{)}{\delta _R}$. Thus, ${\boldsymbol U}_b^\sigma \textrm{(}r\textrm{,}\alpha \textrm{)}$ can be expressed as follows:
Although the summations over M in Eqs. (13) and (12) appear to be the same, the increased value of Rjb compared to Rja indicates that the summation in Eq. (12) is implemented on the annular area enlarged by a Fresnel zone. By calculating the integrals over θ in Eqs. (13) and (14), the co-polarized fields for sets A and B were obtained as follows:
When ja = jb in Eqs. (13) and (14) for sets A and B, respectively, the difference in radius between Rja and Rjb is $\delta R_{ja}^{\textrm{(}b - a\textrm{)}} = {\textrm{(}{R_{jb}} - {R_{ja}}\textrm{)}_{jb = ja}} = \textrm{(}{N_{Mf}} + 1/2\textrm{)}{\delta _R}$, which is also the difference in radius between two adjacent Fresnel zones. Correspondingly, the optical path difference $\delta s_{ja}^{\textrm{(}b - a\textrm{)}} = {\textrm{(}{s_{jb}} - {s_{ja}}\textrm{)}_{jb = ja}}$ between the rings of ja and jb to the central point o of observation plane is $\delta s_{ja}^{\textrm{(}b - a\textrm{)}} = \lambda /2$, and resultantly, ${s_{jb}} = {s_{ja}} + \delta s_{ja}^{\textrm{(}b - a\textrm{)}} = {s_{ja}} + \lambda /2$ for ja = jb. In Eqs. (15) and (16), the diffraction spots of ${J_0}\textrm{(}kr{R_{ja}}/f\textrm{) }$ and ${J_0}\textrm{(}kr{R_{jb}}/f\textrm{) }$ are produced by the slit rings with ${R_{ja}}$ and ${R_{jb}}$, respectively. For the rings within the two adjacent zones, the spot does not vary significantly in dimension, and the approximation ${R_{ja}}{J_0}\textrm{(}kr{R_{ja}}/f\textrm{) } \approx {R_{jb}}{J_0}\textrm{(}kr{R_{jb}}/f\textrm{) }$ is adopted. Therefore, when the metasurface is illuminated by CP light, the co-polarized field ${{\boldsymbol U}^\sigma }\textrm{(}r\textrm{,}\alpha \textrm{)}{ = }{\boldsymbol U}_a^\sigma \textrm{(}r\textrm{,}\alpha \textrm{)} + {\boldsymbol U}_b^\sigma \textrm{(}r\textrm{,}\alpha \textrm{)}$ is calculated as follows:
The derivation of Eq. (17) indicates that the difference in radius between the nanoslit ring sets A and B was the width of a Fresnel zone. The derivation of the above equation indicates that the difference in radius of a nanoslit ring in set A with the corresponding ring in set B for ja = jb was the width of a Fresnel zone. Accordingly, the destructive interference approximately eliminated the overall co-polarization component ${{\boldsymbol U}^\sigma }\textrm{(}r\textrm{,}\alpha \textrm{)}$ near the center point of the observation plane.
We now consider the light fields of CPCs on the observation plane. For the orientation angle φa (Ra, θ) = φh (θ) + φfa (R) of slit sla in set A, the transmitted CPC was imposed with the geometric phase ϕa (R, θ) given in Eqs. (1) and (2). Notably, φfa (R) takes plus sign in Eq. (2) for set A, and the corresponding illumination is RCP. By substituting φa (Ra, θ) in Eq. (11), the CPC of the light field of set A can be obtained as follows:
In the calculations, the geometric phase ϕfa (R) for RCP (σ = -1) illumination cancels the hyperbolic phase profile implicated in ${e^{ik{s_{ja}}}}$ in Eq. (11); therefore, the CPC of set A is focused as expressed by Eq. (18). Although this equation appears simple, the summation therein is practically the same as in Eq. (14). In addition, the integral over θ in Eq. (18) is calculated as ${J_l}\textrm{(}k{R_{ja}}r/f\textrm{)}{e^{ - il\alpha }}$, which is the focused Bessel vortex beam with topological charge of -l = 2σm = -2 m. The superposition of the Bessel function ${J_l}\textrm{(}k{R_{ja}}r/f\textrm{)}$ over Rja is evaluated as an integral over the radial coordinate R. Then, the ${\boldsymbol U}_a^{ - \sigma }\textrm{(}r\textrm{,}\alpha \textrm{)}$ is derived as follows:
This equation indicates that the focused circularly polarized vortex beam of CPC has an initial phase 2σφ0 and that its doughnut profile of amplitude is expressed as [75]
While the LCP light with σ = 1 illuminates set A, the geometric phase is reversed by the opposite chirality σ, resulting in an opposite hyperbolic phase profile that diverges rather than focuses the light field. The light field of the CPC is insignificant and negligible owing to divergence; as a result,
Similarly, when the light of LCP with σ = 1 illuminates the nanoslits in set B, owing to the radial variation of orientation φbr (Rb) opposite to φar (Ra) for set A, the vortex of CPC in RCP is focused; therefore,
When the light of RCP (σ = -1) is used to illuminate nanoslits in set B for the same reason as in the calculation of Eq. (21), we obtain
Thus, selective control of the CPCs in LCP and RCP was realized by the focusing effect using the geometric phases.
2.4 Realization of HOP beams under illumination of elliptical polarizations
We now consider a general case in which the metasurface is illuminated by an elliptically polarized light uin represented as follows:
Because each CP component in uin produces the corresponding co-polarization component and CPC in the transmitted field, the co-polarization components on the observation plane are ${{\boldsymbol U}^{{\sigma _1}}}{ = }{\boldsymbol U}_a^{{\sigma _1}} + {\boldsymbol U}_b^{{\sigma _1}} \approx \textrm{0}$ and ${{\boldsymbol U}^{{\sigma _2}}}{ = }{\boldsymbol U}_a^{{\sigma _2}} + {\boldsymbol U}_b^{{\sigma _2}} \approx \textrm{0}$ based on Eq. (18), and according to Eqs. (21) and (23), the divergent CPC terms are ${\boldsymbol U}_a^{ - {\sigma _1}}\textrm{(}r\textrm{,}\alpha \textrm{) } \approx \textrm{0}$ and ${\boldsymbol U}_b^{ - {\sigma _2}}\textrm{(}r\textrm{,}\alpha \textrm{) } \approx \textrm{0}$. These results, along with those obtained using Eqs. (19) and (22) yield the focused vector field on the observation plane as follows:
3. Simulations
In this study, we first performed the 3D FDTD simulations to demonstrate VB generation, with the mesh size of 0.25 nm. The metasurface was designed using rectangular nanoslits with optimized dimensions of 250 × 80 nm, which were etched into a 200 nm-thick gold film. Identical nanoslits were arranged in a period of 500 nm on the perimeters of the circular rings and the increment in radius of the rings was constant. Each metasurface was designed to contain a total of 24 rings, which were for the nanoslits of sets A and B to be interleaved alternately. A perfectly matched layer was employed as the absorption boundary condition to avoid the influence of adjacent periodic elements. A planar light wave of wavelength 632.8 nm was incident on the 400 nm-thick substrate side of the metasurface. The monitors were set on the air side and focused vector fields were generated on the plane at z = f = 20 μm where their components and total intensities were observed and analyzed.
In Fig. 2(a), Panel (i) shows the image of the designed sample S1. The upper-right inset is the locally enlarged view of the sample. The sample diameter is 24 μm. Nanoslit sets A and B of the sample have the same rotation order of m = 0.5; they focus the vortices of the cross- LCP and RCP with topological charges of -1 and 1, respectively, and their superposition generates VBs of order l = 1. The lower-right inset is the designed phase map φa (R, θ) for the focused VB. Therein the dotted circle denotes the initial ring of set A, and the area is filled with faded colors, instead of the original blank space, for convenient comparison with the familiar focused vortex phase map. Panel (ii) in Fig. 2(a) shows the x-z cross-section image of the total intensity I = |Ex|2 +|Ey|2 +|Ez|2 of the radially polarized VB generated by S1 under illumination of horizontal linear polarization (LP). The unneglectable z-component intensity |Ez|2 resulted in a nonzero distribution in the image around the z-axis, which is the basic characteristic of focusing radially polarized VBs near the focus [50]. Panel (iii) in Fig. 2(a) shows the curves of the total and component intensities versus x. From the results, we see that an obvious longitudinal component field Ez was obtained, and its curve of Iz = |Ez|2 exhibited a focused profile with full width at half-maximum (FWHM) of 0.91λ. Meanwhile, the intensity curve of the radial component Ir has a hollow-core profile with FWHM 0.70λ. Thus, the intensity I = |Ex|2+ |Ey|2+ |Ez|2 exhibited a double-peak curve. Though it does not have the single sharp peak as the tightly focused radially polarized VBs, whose Iz curve usually has a peak value several times higher than that of the Ir curve, the peak value of its Iz curve reaches half that of Ir curve, demonstrating that the generated radially polarized VB is well focused. Besides, the tightness for focusing the VB may be reasonably enhanced by increase the size of the metasurface. Panels (iv–vii) in Fig. 2(a) from left to right indicate the simulated images of IA = |Ex|2 + |Ey|2, Ix = |Ex|2 and Iy = |Ey|2 for the four VBs generated by S1 when illuminating light are horizontally, 45°, vertically, and 135°-linearly polarized from top to bottom, respectively. The purple double arrows in the title column indicate the direction of polarization of incident light, and the black arrows overlaid on the intensity doughnuts in the first column indicate the polarization state of the VBs. The light intensity images of x and y-component fields exhibit the familiar petal-like distributions, and the two lobes in the x-component images of different VBs are oriented and rotate in the same direction as the corresponding incident polarization.
Panel (i) in Fig. 2(b) is the image of sample S2 in which the two sets of nanoslits have the same rotation order m = -0.5. The topological charges of the corresponding focused vortices of RCP and LCP are l = -1 and 1, respectively, thus generating the π-phased VBs with l = -1. Panel (ii) in Fig. 2(b) shows the x-z cross-sectional image of the total intensity of the π-phased VB of radial polarization under the illumination of horizontally linear polarization. The hollow intensity distribution indicates the zero value of |Ez|2 ∼ 0 near the vortex core for this VB. As shown in panel (iii) in Fig. 2(b), for the π-radially polarized VB, the intensity curve of Iφ has a profile with a hollow core with FWHM of 0.79λ, and Iz curve exhibits the hollow core and does not have focused central spot. This demonstrating the basic characteristics of the π-radially polarized VBs. Panels (iv–vii) in Fig. 2(b) show the intensity images of VBs of order l = -1; i.e., the π-radial, π-135°, π-azimuthal, and π-45° polarizations, respectively in the focal plane. Noticeably, the orientations of the two lobes in the x-component images were observed to rotate in the direction opposite to the corresponding incident polarizations, which differs from the previously discussed VBs. In addition to the numerical simulation of the VBs of order l = 1 and -1 (i.e., π-phase VBs) generated by the two samples S1 and S2, the light fields of all the samples in this study were simulated and optimized by FDTD before sample fabrications, although no simulation results are presented.
4. Experiment
4.1 Optical measurement setup and sample fabrications
Figure 3(a) shows a schematic of the experimental setup. A He-Ne laser with a wavelength of 632.8 nm is used as the light source. The incident light is passed through an attenuator (A) for power adjustment, which illuminates the sample from the side of the fused silica substrate, and the focused VBs are observed on the plane at a distance of 20 μm from the sample surface. The focused VBs are enlarged and imaged using a microscope objective (MO) and the enlarged images are recorded on a scientific complementary metal–oxide semiconductor (sCMOS). The analyzing polarizer (P) is placed in front of the sCMOS to obtain x-, 45°-, y-, and 135°- component intensity images. A half-wave plate (HP) and quarter-wave plate (QWP) are placed in front of the sample and combined to adjust the major-axis direction and ellipticity of the incident light of elliptical polarization corresponding to an arbitrary point on the traditional Poincaré sphere.
In the sample fabrication, a 200 nm-thick gold film was evaporated on the 1 mm-thick fused silica substrate using magnetron sputtering and nanoslits were etched on the gold film using the focused ion beam (FIB) method. Seven samples were fabricated and denoted as samples S1–S7. The parameters of the fabricated samples are presented in Table 1. Figure 3(b) shows images of sample S1 generated through scanning electron microscopy (SEM), whose rotation order is m = 0.5 for generating VBs of order l = 1. Inset shows a locally enlarged view. The rotation orders of samples S2, S3, and S4 are m = -0.5, 1, and 1.5, respectively and their initial orientation angles φ0 = 0. Correspondingly, they generate VBs of orders l = -1, 2, and 3 on the observation plane. Samples S5, S6, and S7 have the same rotation order as sample S1, i.e., m = 0.5; however, they have different initial orientation angles φ0. Under the illuminating light of certain polarizations, they can generate VBs of the order l = l with different polarization states.
4.2 VBs of different orders with linear polarized states
First, using the illumination of LP and changing the polarization direction (i.e., the polarization angle Φ/2 of the LP), the VBs were generated on the equator of the HOP sphere. Figure 4(a) shows a Poincaré sphere. For simplicity and convenience, we used this sphere to discuss the illuminating light on the traditional Poincaré sphere and the generated VBs on the HOP spheres of different orders. Figures 4(b)–(e) present the images of VBs of orders l = 1, -1, 2, and 3 on the equators of the HOP spheres generated by samples S1-S4, respectively. The purple double arrows to the left of the images indicate the polarizations of the incident light corresponding to the four equidistant points A, B, C, and D moving right-handedly on the equator of the conventional Poincaré sphere. The green double arrows in the dotted circles in the top title row of the figures indicate the transmission direction of the analyzed polarizer. The images in the first column are theoretical doughnuts, the overlaid line segments are schematics of the polarization distributions of the corresponding VBs, and images of the component intensities of the VBs generated by each sample are presented in the second to fifth columns.
Comparing the component images of samples S1 (m = 0.5) in Fig. 4(b) and S2 (m = -0.5) in Fig. 4(c), we observed that the intensity lobes of the π-phase VBs in the latter rotate with the analyzing polarizer opposite to the ordinary VBs in the former. In addition, the polarization distributions in Figs. 4(b) and (c) indicate that, for both the VBs of orders l = 1 and l = -1, the four points A’, B’, C’, and D’ representing the polarization states on their HOP spheres, respectively, move oppositely to the right-hand direction of points A, D, C, and B on the conventional Poincaré sphere. Particularly, the intensity images in the first and third rows in both Fig. 4(b) and (c) are the radial, azimuthal, π-radial, and π-azimuthal VBs of the Bell-like states [76], respectively, which correspond to the simulation results shown in Figs. 2(a) and (b). The measured results for all VBs were in good consistency.
Figures 4(d) and (e) show the focused VBs of l = 2 and 3 generated by samples S3 and S4 with rotation orders m = 1 and 1.5, respectively. The intensity images from rows 2 to 5 in each figure for the VBs of radial, 135°-slanted, azimuthal, and 45°-slanted polarizations also corresponded to points A’, B’, C’, and D’ on the equators of the corresponding HOP spheres, respectively. They exhibit properties similar to those of the VBs produced by sample S1 and are consistent with the familiar theoretical results.
4.3 HOP beams of different orders in different longitudes
For demonstrating the focused HOP beams by superposing the two CP vortices of non-equal weights, samples S1 and S4 were used as examples to generate VBs at nine points on the longitudes Φ'/2 = 135°, and Φ'/2 = 45°, on the corresponding HOP spheres of l = 1 and l = 3 as depicted in the left panels of Figs. 5(a) and (b), respectively. The corresponding theoretical, simulation, and experimental results are presented in the right panels of Figs. 5(a) and (b). On the right panel of each figure, the elliptical arrows in top title row denote the elliptical polarizations of the illuminating light, corresponding to the latitude Θ with sin Θ varying from 1 to -1 by an increment of -0.25. The images in the first and third rows are doughnuts of the theoretical and experimental total intensities of IA, respectively, and the schematic of the polarization states is overlaid on the first row. The second and fourth rows show the simulated and experimental images of x-component intensities, respectively, and the corresponding values of latitudes on the HOP spheres are shown in the second row.
These images show the evolution of the HOP beams at points from the south-pole point I to the north-pole point A on the longitudes of the corresponding HOP spheres. Starting from the LCP vortex beam at south-pole point I, the VBs successively underwent elliptical polarizations of left-handedness in the southern hemisphere, LP at the equator point, elliptical polarizations of right-handedness in the northern hemisphere, and ended at the north-pole point A of the RCP vortex beam. According to the transformation relation (Θ, Φ) → (Θ’ = π-Θ, Φ’ = 2φ0-Φ) in the above process, incident light transited in a direction opposite to RCP at the north pole to LCP at the south pole along the longitudes Θ = π-Θ’ on the conventional Poincaré sphere. As shown in the experimental images in Fig. 5(a), the lobes in the x-component images have the maximum contrast at the equator point E; however, they gradually blurred with decreasing contrast when the polarization state points moved toward the two pole points A and I. Similarly, the component intensity images of the HOP beams with l = 3 in Fig. 5(b) evolve in the same manner, but the x-component intensity images of the VBs contain 2 l = 6 lobes.
4.4 VBs generated by controlling initial orientations of nanoslits
To demonstrate manipulating polarization states of the VBs by the initial orientation angle of the nanoslits, we designed and fabricated samples S1, S5, S6, and S7 as a group with initial angles of 0, π/8, π/4, and 3π/8. According to Eq. (26), when the incident light kept horizontally polarized, the four samples generated the VBs of radial, 135°-slanted, azimuthal, and 45°-slanted polarizations, respectively. Figure 6 shows the SEM images of the four samples and experimental intensity images of the corresponding VBs. The upper-left panels in Figs. 6(a)–(d) show the SEM images and enlarged views of the four samples, respectively, in which the initial orientation of the first nanoslit in set A is marked. The upper-right panels of Figs. 6(a)–(d) depict images of the total intensity IA with overlaid schematics of the polarization states. The lower-left and right panels of Figs. 6(a)–(d) depict images of the component intensities Ix and Iy, respectively. These images show that samples S1, S5, S6, and S7 generated radial, 135°-slanted, azimuthal, and 45°-slanted VBs, respectively, which are consistent with the theoretical results. These images demonstrate that the initial orientation angles of the samples can control the polarization states of the generated VBs. Thus, using the metasurfaces designed in this study, focused VBs were experimentally generated with satisfactory beam quality.
5. Discussion
In the metasurface design, we assumed that in Eq. (17), in a pair of Fresnel zones, the sum of the co-polarization components of the nanoslits in sets A and B is approximately zero. This was realized by the destructive interference arising from the increase of optical path of nanoslit set B by a half-wavelength. Strictly speaking, the increment in radius of the Fresnel zones varied with the radius such that the number of rings 2NMf of sets A and B in certain adjacent zones may not be equal. However, their co-polarization components can be considered along with those of the other two similar adjacent zones. Thus, based on the overall summation, the co-polarization components in Eq. (17) can be eliminated. In the practical design of metasurfaces, this was realized by using FDTD simulations to determine the ring numbers of sets A and B and the total ring number 2N of the metasurfaces. The diameter of the metasurfaces is 24 μm, and each metasurface include 2352 nanoslits. Additionally, other factors that adversely affect VB generation, including practical errors in the amplitude and small defocusing of the two vortices, were finely optimized by the simulation. The metasurface parameters were finally determined as R1a = 3.5 μm, R1b = 4.95 μm, δR = 0.6 μm, NA = NB = 12. It is noted that under the preset focal length of 20 μm, the radii of the first and second Fresnel zone were calculated as 3.572 μm and 5.071 μm, respectively, and they were the theoretical values we set for R1a and R1b. In fact, the nanoslit would impose an inherent phase retardation on the transmitted waves [77]. To compensate this retardation, R1a and R1b were finely adjusted, and the optimization were performed using the FDTD simulations of the VBs. Finally, to further check the consistency of the experimental results with the metasurface designs, we measured the focal length of the metasurfaces. The measurement was performed using a piezoelectric nano-translation stage (PI E-516) with precision of 10 nm, which could finely move the metasurface longitudinally. The focal length of the sample S1 was measured 18.52 ± 0.21 μm, demonstrating a satisfactory consistency.
6. Conclusion
This study demonstrated the generation of focused coaxial HOP beams based on the proposed metasurface design with theoretical and experimental implementation. Metasurfaces were fabricated on gold films etched with nanoslits on equidistant circular rings. Using the geometric phases of two interleaved sets of nanoslits on a metasurface, two vortices with opposite CPs were generated and the focusing were realized simultaneously. By controlling the initial radii of the two sets of nanoslits to produce an optical path difference of half-wavelength, the opposite phase of the overall co-polarization fields of the two nanoslit sets was achieved. By the optimized designs with FDTD simulations, the co-polarization components were eliminated and focused coaxial HOP beams were generated. The proposed metasurface design provided a flexible and feasible approach for the manipulation of focused HOP beams. This study has significance in integrated and miniaturized optical devices and has potential applications in versatile areas in classical and quantum sciences.
Funding
National Natural Science Foundation of China (12174226, 62175134, 62375159, 12004215); Natural Science Foundation of Shandong Province (ZR2022MF248).
Disclosures
The authors declare no conflicts of interest.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
References
1. C. Rosales-Guzmán, B. Ndagano, and A. Forbes, “A review of complex vector light fields and their applications,” J. Opt. 20(12), 123001 (2018). [CrossRef]
2. Q. Zhan, “Trapping metallic Rayleigh particles with radial polarization,” Opt. Express 12(15), 3377–3382 (2004). [CrossRef]
3. X. Xie, Y. Chen, K. Yang, et al., “Harnessing the point-spread function for high-resolution far-field optical microscopy,” Phys. Rev. Lett. 113(26), 263901 (2014). [CrossRef]
4. V. D’Ambrosio, N. Spagnolo, L. Del Re, et al., “Photonic polarization gears for ultra-sensitive angular measurements,” Nat. Commun. 4(1), 2432 (2013). [CrossRef]
5. L. J. Wong and F. X. Kärtner, “Direct acceleration of an electron in infinite vacuum by a pulsed radially-polarized laser beam,” Opt. Express 18(24), 25035–25051 (2010). [CrossRef]
6. V. D’Ambrosio, G. Carvacho, F. Graffitti, et al., “Entangled vector vortex beams,” Phys. Rev. A 94(3), 030304 (2016). [CrossRef]
7. G. Milione, H. I. Sztul, D. A. Nolan, et al., “Higher-order Poincaré sphere, Stokes parameters, and the angular momentum of light,” Phys. Rev. Lett. 107(5), 053601 (2011). [CrossRef]
8. A. Beckley, T. Brown, and M. Alonso, “Full Poincaré beams,” Opt. Express 18(10), 10777–10785 (2010). [CrossRef]
9. B. Zhao, X. B. Hu, V. Rodríguez-Fajardo, et al., “Determining the non-separability of vector modes with digital micromirror devices,” Appl. Phys. Lett. 116(9), 091101 (2020). [CrossRef]
10. Z. Wan, Y. Shen, Q. Liu, et al., “Multipartite classically entangled scalar beams,” Opt. Lett. 47(8), 2052–2055 (2022). [CrossRef]
11. T. M. Graham, H. J. Bernstein, T. C. Wei, et al., “Superdense teleportation using hyperentangled photons,” Nat. Commun. 6(1), 7185 (2015). [CrossRef]
12. J. Alicea and P. Fendley, “Topological phases with parafermions: theory and blueprints,” Annu. Rev. Condens. Matter Phys. 7(1), 119–139 (2016). [CrossRef]
13. H. K. Lo and H. F. Chau, “Unconditional security of quantum key distribution over arbitrarily long distances,” Science 283(5410), 2050–2056 (1999). [CrossRef]
14. V. D’Ambrosio, E. Nagali, S. P. Walborn, et al., “Complete experimental toolbox for alignment-free quantum communication,” Nat. Commun. 3(1), 961 (2012). [CrossRef]
15. G. Vallone, V. D’Ambrosio, A. Sponselli, et al., “Free-space quantum key distribution by rotation-invariant twisted photons,” Phys. Rev. Lett. 113(6), 060503 (2014). [CrossRef]
16. V. Parigi, V. D’Ambrosio, C. Arnold, et al., “Storage and retrieval of vector beams of light in a multiple-degree-of-freedom quantum memory,” Nat. Commun. 6(1), 7706 (2015). [CrossRef]
17. D. Naidoo, F. S. Roux, A. Dudley, et al., “Controlled generation of higher-order Poincaré sphere beams from a laser,” Nat. Photon. 10(5), 327–332 (2016). [CrossRef]
18. C. Chen, Y. Zhang, L. Ma, et al., “Flexible generation of higher-order Poincaré beams with high efficiency by manipulating the two eigenstates of polarized optical vortices,” Opt. Express 28(7), 10618–10632 (2020). [CrossRef]
19. P. Li, Y. Zhang, S. Liu, et al., “Generation of perfect vectorial vortex beams,” Opt. Lett. 41(10), 2205–2208 (2016). [CrossRef]
20. P. Chen, S. J. Ge, W. Duan, et al., “Digitalized geometric phases for parallel optical spin and orbital angular momentum encoding,” ACS Photonics 4(6), 1333–1338 (2017). [CrossRef]
21. X. L. Wang, J. Ding, W. J. Ni, et al., “Generation of arbitrary vector beams with a spatial light modulator and a common path interferometric arrangement,” Opt. Lett. 32(24), 3549–3551 (2007). [CrossRef]
22. Z. Y. Rong, Y. J. Han, S. Z. Wang, et al., “Generation of arbitrary vector beams with cascaded liquid crystal spatial light modulators,” Opt. Express 22(2), 1636–1644 (2014). [CrossRef]
23. W. T. Chen, A. Y. Zhu, V. Sanjeev, et al., “A broadband achromatic metalens for focusing and imaging in the visible,” Nat. Nanotechnol. 13(3), 220–226 (2018). [CrossRef]
24. W. Zang, Q. Yuan, R. Chen, et al., “Chromatic dispersion manipulation based on metalenses,” Adv. Mater. 32(27), e1904935 (2020). [CrossRef]
25. J. P. Balthasar Mueller, N. A. Rubin, R. C. Devlin, et al., “Metasurface polarization optics: independent phase control of arbitrary orthogonal states of polarization,” Phys. Rev. Lett. 118(11), 113901 (2017). [CrossRef]
26. M. Liu, W. Zhu, P. Huo, et al., “Multifunctional metasurfaces enabled by simultaneous and independent control of phase and amplitude for orthogonal polarization states,” Light: Sci. Appl. 10(1), 107 (2021). [CrossRef]
27. Q. Fan, M. Liu, C. Zhang, et al., “Independent amplitude control of arbitrary orthogonal states of polarization via dielectric metasurfaces,” Phys. Rev. Lett. 125(26), 267402 (2020). [CrossRef]
28. E. Maguid, I. Yulevich, M. Yannai, et al., “Multifunctional interleaved geometric-phase dielectric metasurfaces,” Light: Sci. Appl. 6(8), e17027 (2017). [CrossRef]
29. F. Yue, D. Wen, J. Xin, et al., “Vector vortex beam generation with a single plasmonic metasurface,” ACS Photonics 3(9), 1558–1563 (2016). [CrossRef]
30. F. Ding, Y. Chen, Y. Yang, et al., “Multifunctional metamirrors for broadband focused vector-beam generation,” Adv. Opt. Mater. 7(22), 1900724 (2019). [CrossRef]
31. M. Liu, P. Huo, W. Zhu, et al., “Broadband generation of perfect Poincaré beams via dielectric spin-multiplexed metasurface,” Nat. Commun. 12(1), 2230 (2021). [CrossRef]
32. S. Gao, C. S. Park, S. S. Lee, et al., “A highly efficient bifunctional dielectric metasurface enabling polarization-tuned focusing and deflection for visible light,” Adv. Opt. Mater. 7(9), 1801337 (2019). [CrossRef]
33. S. Li, X. Li, G. Wang, et al., “Multidimensional manipulation of photonic spin Hall effect with a single-layer dielectric metasurface,” Adv. Opt. Mater. 7(5), 1801365 (2019). [CrossRef]
34. Y. J. Gao, Z. Wang, Y. Jiang, et al., “Multichannel Distribution and Transformation of Entangled Photons with Dielectric Metasurfaces,” Phys. Rev. Lett. 129(2), 023601 (2022). [CrossRef]
35. J. Zhang, J. Ma, M. Parry, et al., “Spatially entangled photon pairs from lithium niobate nonlocal metasurfaces,” Sci. Adv. 8(30), eabq4240 (2022). [CrossRef]
36. Y. Yuan, K. Zhang, B. Ratni, et al., “Independent phase modulation for quadruplex polarization channels enabled by chirality-assisted geometric-phase metasurfaces,” Nat. Commun. 11(1), 4186 (2020). [CrossRef]
37. P. Xu, H. Liu, R. Li, et al., “Dual-Band Spin-Decoupled Metasurface for Generating Multiple Coaxial OAM Beams,” IEEE Trans. Antennas Propag. 70(11), 10678–10690 (2022). [CrossRef]
38. J. Zhu, Y. Yang, J. Lai, et al., “Additively Manufactured Polarization Insensitive Broadband Transmissive Metasurfaces for Arbitrary Polarization Conversion and Wavefront Shaping,” Adv. Opt. Mater. 10(21), 2200928 (2022). [CrossRef]
39. C. Zheng, J. Li, J. Liu, et al., “Creating longitudinally varying vector vortex beams with an all-dielectric metasurface,” Laser Photonics Rev. 16(10), 2200236 (2022). [CrossRef]
40. J. Li, J. Li, Z. Yue, et al., “Structured vector field manipulation of terahertz wave along the propagation direction based on dielectric metasurfaces,” Laser Photonics Rev. 16(12), 2200325 (2022). [CrossRef]
41. Z. Yue, J. Liu, J. Li, et al., “Vector beam generation based on spin-decoupling metasurface zone plate,” Appl. Phys. Lett. 120(19), 191704 (2022). [CrossRef]
42. T. Stav, A. Faerman, E. Maguid, et al., “Quantum entanglement of the spin and orbital angular momentum of photons using metamaterials,” Science 361(6407), 1101–1104 (2018). [CrossRef]
43. M. Wang, L. Chen, D. Y. Choi, et al., “Characterization of Orbital Angular Momentum Quantum States Empowered by Metasurfaces,” Nano Lett. 23(9), 3921–3928 (2023). [CrossRef]
44. J. Ji, Z. Wang, J. Sun, et al., “Metasurface-Enabled On-Chip Manipulation of Higher-Order Poincaré Sphere Beams,” Nano Lett. 23(7), 2750–2757 (2023). [CrossRef]
45. H. Hua, C. Zeng, Z. He, et al., “Plasmonic metafiber for all-fiber Q-switched cylindrical vector lasers,” Nanophotonics 12(4), 725–735 (2023). [CrossRef]
46. Z. L. Deng, M. Jin, X. Ye, et al., “Full-Color Complex-Amplitude Vectorial Holograms Based on Multi-Freedom Metasurfaces,” Adv. Funct. Mater. 30(21), 1910610 (2020). [CrossRef]
47. F. Ding, Y. Chen, and S. I. Bozhevolnyi, “Gap-surface plasmon metasurfaces for linear-polarization conversion, focusing, and beam splitting,” Photonics Res. 8(5), 707–714 (2020). [CrossRef]
48. J. T. Heiden, F. Ding, J. Linnet, et al., “Gap-surface plasmon metasurfaces for broadband circular-to-linear polarization conversion and vector vortex beam generation,” Adv. Opt. Mater. 7(9), 1801414 (2019). [CrossRef]
49. Y. Deng, C. Wu, C. Meng, et al., “Functional metasurface quarter-wave plates for simultaneous polarization conversion and beam steering,” ACS Nano 15(11), 18532–18540 (2021). [CrossRef]
50. X. Zeng, Y. Zhang, M. Gu, et al., “Arbitrary manipulations of focused higher-order Poincaré beams by a Fresnel zone metasurface with alternate binary geometric and propagation phases,” Photonics Res. 10(4), 1117–1126 (2022). [CrossRef]
51. Y. Bao, J. Ni, and C.-W. Qiu, “A minimalist single-layer metasurface for arbitrary and full control of vector vortex beams,” Adv. Mater. 32(6), e1905659 (2020). [CrossRef]
52. R. C. Devlin, A. Ambrosio, N. A. Rubin, et al., “Arbitrary spin-to–orbital angular momentum conversion of light,” Science 358(6365), 896–901 (2017). [CrossRef]
53. A. Arbabi, Y. Horie, M. Bagheri, et al., “Dielectric metasurfaces for complete control of phase and polarization with subwavelength spatial resolution and high transmission,” Nat. Nanotechnol. 10(11), 937–943 (2015). [CrossRef]
54. P. Yu, S. Chen, J. Li, et al., “Generation of vector beams with arbitrary spatial variation of phase and linear polarization using plasmonic metasurfaces,” Opt. Lett. 40(14), 3229–3232 (2015). [CrossRef]
55. Z. Wu, F. Dong, S. Zhang, et al., “Broadband dielectric metalens for polarization manipulating and superoscillation focusing of visible light,” ACS Photonics 7(1), 180–189 (2020). [CrossRef]
56. F. Yue, D. Wen, C. Zhang, et al., “Multichannel polarization-controllable superpositions of orbital angular momentum states,” Adv. Mater. 29(15), 1603838 (2017). [CrossRef]
57. T. Li, X. Li, S. Yan, et al., “Generation and conversion dynamics of dual Bessel beams with a photonic spin-dependent dielectric metasurface,” Phys. Rev. Appl. 15(1), 014059 (2021). [CrossRef]
58. J. Yang, T. K. Hakala, and A. T. Friberg, “Generation of arbitrary vector Bessel beams on higher-order Poincaré spheres with an all-dielectric metasurface,” Phys. Rev. A 106(2), 023520 (2022). [CrossRef]
59. Z. H. Jiang, L. Kang, T. Yue, et al., “A single noninterleaved metasurface for high-capacity and flexible mode multiplexing of higher-order Poincaré sphere beams,” Adv. Mater. 32(6), e1903983 (2020). [CrossRef]
60. D. Wen, K. Pan, J. Meng, et al., “Broadband multichannel cylindrical vector beam generation by a single metasurface,” Laser Photonics Rev. 16(10), 2200206 (2022). [CrossRef]
61. S. Wang, S. Wen, Z.-L. Deng, et al., “Metasurface-Based Solid Poincaré Sphere Polarizer,” Phys. Rev. Lett. 130(12), 123801 (2023). [CrossRef]
62. R. Zuo, W. Liu, H. Cheng, et al., “Breaking the diffraction limit with radially polarized light based on dielectric metalenses,” Adv. Opt. Mater. 6(21), 1800795 (2018). [CrossRef]
63. J. Yang, S. Gurung, S. Bej, et al., “Active optical metasurfaces: comprehensive review on physics, mechanisms, and prospective applications,” Rep. Prog. Phys. 85(3), 036101 (2022). [CrossRef]
64. E. Pisanty, G. J. Machado, V. Vicuña-Hernández, et al., “Knotting fractional-order knots with the polarization state of light,” Nat. Photonics 13(8), 569–574 (2019). [CrossRef]
65. H. Larocque, D. Sugic, D. Mortimer, et al., “Reconstructing the topology of optical polarization knots,” Nat. Phys. 14(11), 1079–1082 (2018). [CrossRef]
66. G. Makey, Ö. Yavuz, D. K. Kesim, et al., “Breaking crosstalk limits to dynamic holography using orthogonality of high-dimensional random vectors,” Nat. Photonics 13(4), 251–256 (2019). [CrossRef]
67. E. Wang, L. Shi, J. Niu, et al., “Multichannel spatially nonhomogeneous focused vector vortex beams for quantum experiments,” Adv. Opt. Mater. 7(8), 1801415 (2019). [CrossRef]
68. C. Zheng, G. Wang, J. Li, et al., “All-dielectric metasurface for manipulating the superpositions of orbital angular momentum via spin-decoupling,” Adv. Opt. Mater. 9(10), 2002007 (2021). [CrossRef]
69. M. Khorasaninejad, W. Chen, R. C. Devlin, et al., “Metalenses at visible wavelengths: Diffraction-limited focusing and subwavelength resolution imaging,” Science 352(6290), 1190–1194 (2016). [CrossRef]
70. R. Zhang, M. Gu, R. Sun, et al., “Plasmonic metasurfaces manipulating the two spin components from spin–orbit interactions of light with lattice field generations,” Nanophotonics 11(2), 391–404 (2022). [CrossRef]
71. R. Fu, L. Deng, Z. Guan, et al., “Zero-order-free meta-holograms in a broadband visible range,” Photonics Res. 8(5), 723–728 (2020). [CrossRef]
72. G. Zheng, R. Fu, L. Deng, et al., “On-axis three-dimensional meta-holography enabled with continuous-amplitude modulation of light,” Opt. Express 29(4), 6147–6157 (2021). [CrossRef]
73. Y. Hu, X. Liu, M. Jin, et al., “Dielectric metasurface zone plate for the generation of focusing vortex beams,” PhotoniX 2(1), 1–9 (2021). [CrossRef]
74. T. V. Teperik, A. Archambault, F. Marquier, et al., “Huygens-Fresnel principle for surface plasmons,” Opt. Express 17(20), 17483–17490 (2009). [CrossRef]
75. V. V. Kotlyar, S. N. Khonina, A. A. Kovalev, et al., “Diffraction of a plane, finite-radius wave by a spiral phase plate,” Opt. Lett. 31(11), 1597–1599 (2006). [CrossRef]
76. L. J. Pereira, A. Z. Khoury, and K. Dechoum, “Quantum and classical separability of spin-orbit laser modes,” Phys. Rev. A 90(5), 053842 (2014). [CrossRef]
77. X. Li, R. Zhang, Y. Zhang, et al., “Slit width oriented polarized wavefields transition involving plasmonic and photonic modes,” New J. Phys. 20(6), 063037 (2018). [CrossRef]