Quasi- Perfect Vortices Generated by Pancharatnam-Berry Phase Metasurfaces for Optical Spanners and OAM Communication

doi:10.21203/rs.3.rs-997592/v1

Download PDF

Research Article

Quasi- Perfect Vortices Generated by Pancharatnam-Berry Phase Metasurfaces for Optical Spanners and OAM Communication

https://doi.org/10.21203/rs.3.rs-997592/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Optical vortex (OV) can be used in the fields of optical manipulation and optical communication because of its inherent orbital angular momentum (OAM). The size of the OV ring increases with the correlated topological charge (TC), making the OV with large TC not suitable for optical rotation and short-distance communication. Perfect vortex (PV) has attracted much attention due to that its optical transmission profile is almost independent of TC. In this manuscript, we proposed a method to generate quasi- perfect vortices (Q-PVs) by Pancharatnam–Berry (PB) phase metasurfaces, the so-called Q-PV can be regarded as an annularly focused optical vortex with spiral phase countermarked on the ring. It has a similar property to PV in that its light profile hardly changes with TC in the focal plane. We demonstrated that the Q-PV can be used for optical spanners that particles are trapped and rotated on the specific orbit. Non-coaxial and coaxial Q-PV arrays were further generated for OAM communication applications. We believe that the proposed Q-PVs has potential applications in optical manipulation and optical communication.

Scientific Communication

Electronic Materials and Devices

Electrical Engineering

Optical vortex (OV)

orbital angular momentum (OAM)

ancharatnam-Berry Phase Metasurfaces

Optical Spanners

Optical vortex (OV) as a kind of Laguerre-Gaussian (LG) mode beam with spiral phase and orbital angular momentum (OAM), has main applications in the fields of particle trapping and manipulation [1–3] and optical communication [4–6]. OV has singularity with a donut shape, leading to destructive interference at the central light field. Due to the oblique feature of the Poynting vector in the OAM field, the trapped particles will rotate along the OAM orbit field [7, 8]. In the past several years, OV has been developed as an optical spanner and received wide attention [9–13]. Various kinds of optical spanners such as focused optical spanners [14], plasmonic spanners have been developed [15]. The complex amplitude of OV has Hilbert factor exp(ilθ), where \(l\) is the number of topological charge (TC), and \(\theta\) is the azimuth angle [16]. TC is the key parameter of OV, and it is proportional to the OAM and the radius of light field. Because of the changeability of OAM, OV has a wide range applications in the field of optical communication, such as quantum information coding [17, 18], wavelength division multiplexing [19], Dammann grating splitting [20, 21]. However, the ring radius of the vortex intensity profile strongly depends on the TC, thus limiting its application on fiber coupling and multiplexing. OAM sources with changeable TCs as well as fixed rings are highly desirable for optical manipulation and communications [22].

Perfect vortex (PV), whose ring profile and divergence angle keep the same size under different TCs [23–25], could provide a perfect solution. Traditional PVs were generated by the Fourier transform of a Bessel-Gaussian beam [26] by using several bulky optical elements including conical lens, spiral phase plate or SLM, and lens [24, 27, 28]. Alternatively, PV has also been generated by using conical lens to generate annular beam and spatial light modulator (SLM) to mark vortex phase, and it has been used to achieve fixed orbit OAM manipulation in the tight focal field [27]. In 2014, a similar beam was generated by SLM imprinting the phases of the diffraction ring lens and spiral phase plate [29]. However, the generation of these beams requires complex and bulky optical elements and precise light path adjustment.

Recently, the study of metasurfaces makes optical devices tend to be flattened and integrated. Dielectric metasurfaces consist of nanofins in subwavelength scale, providing a powerful and efficient platform for light field regulation in cases of amplitude, phase, and polarization [30, 31]. PB phase optical elements use the idea of geometric phase with polarization transmission, which recently has been used in the applications of OV generation [32, 33] and OV spanners [33, 34]. There have appeared some excellent works for PV generation by using PB phase metalens. In 2017, Liu proposed the generation of PV by the combination of three metalenses with convex lens phase, Bessel phase, and vortex phase [22]. In 2018, Yang combined these phases on one single plasmonic metasurface based on PB phase to obtain a 3D broadband PV [35].

In this manuscript, the dielectric metasurfaces based on PB phase were employed to combine the annular focusing phase and OV phase to generate quasi- perfect vortex (Q-PV). It is expected that the generated Q-PVs with different TCs have fixed light field profiles in the focal planes. The generated Q-PV has technically eliminated the Bessel phase, which generally exists in traditional PV due to the Bessel series expansion [22, 35]. With the generated Q-PV, 3D optical trapping can be achieved since the size and position of the focal ring are adjustable, while only 2D trapping can be achieved with traditional PV due to the non-diffraction propagation property. Besides, we can generate a new kind of optical spanner in that particles are rotated on a fixed orbit with different torques. The Q-PV can be used in OAM communication for integrated photonic devices. Since the metasurface can easily form an array, either coaxial or non-coaxial communication can be achieved. All in all, the proposed Q-PV may broaden the applications of optical information coding and storage, photon computing, micromechanics, and lab-on-a-chip.

Theoretical description of Q-PV.

Our proposed Q-PV is defined by

\(E\left(r,\theta ,f\right)=\delta \left(r-{r}_{0}\right)\bullet \text{e}\text{x}\text{p}\left(il\theta \right), \left(l=1, 2, 3\dots \right)\). \(\left(1\right)\)

where f is the designed focal length of Q-PV, r, θ are the polar coordinates in the beam cross-section in the focal plane, δ(*) is the Dirac δ-function, l is the TC and r₀ is the ring radius of the light field. The part after the equal sign keeps exactly the same form as that of PV [25]. It means the Q-PV has the same property with PV in the focal plane that the maximum intensity ring at r₀ has noting with TC.

Function δ(*) can be expanded based on the principle of generalized Snell's law [36], and the transmission function of Q-PV can be expressed as

\(T\left(r,\theta \right)=\text{e}\text{x}\text{p}\left(ik\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\right)\bullet \text{e}\text{x}\text{p}\left(il\theta \right)\), \(\left(2\right)\)

where \(k=2\pi /\lambda\) is the wave vector at the working wavelength λ, θ is the azimuth angle, l is the number of TC, respectively. Fresnel diffraction integral theory was used to accurately calculate the light field distribution. The field strength of Q-PV on the z-transmission plane can be expressed as

\(E\left(r,\rho ,z\right)=\frac{k}{z}{\int }_{0}^{{\rho }_{0}}\text{e}\text{x}\text{p}\left(ik\sqrt{{r}^{2}+{z}^{2}}-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\bullet {J}_{l}\left(\frac{k\rho r}{z}\right)rdr\), \(\left(3\right)\)

where ρ is the integral radius of the aperture, ρ₀ is the lens size, and J_l(*) is the l^th Bessel function, which can be considered as the combination of the circular aperture diffraction of the PB phase metasurface and the l^th OV light.

The metasurface design.

We used an all-dielectric PB phase metasurface consisting of elliptic nanoposts [36] array on a SiO₂ substrate, as shown in figure 1(a). A right-handed circularly polarized (RCP) Gaussian light normally illuminated the SiO₂ substrate from the bottom up, passing through the nanoposts array and form a Q-PV with LCP polarization state [37]. The beam’s working wavelength was set with 532 nm and the waist radius was twice the radius of the metasurface. To uniformly produce a Q-PV, the nanoposts on the SiO₂ substrate were arranged in ring shape. According to PB phase principle, the rotation angle of each nanopost should be \(\alpha =\phi ⁄2\). Extract the phase term of equation 2 and accept the phase-angle conversion relationship, we obtain the following expressions for the phase and rotation angle distributions, respectively:

\(\phi \left(r,l,\theta \right)=k\left(\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\right)+l\bullet \theta\) , \(\left(4\right)\)

\(\alpha \left(r,l,\theta \right)=\frac{k}{2}\left(\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\right)+\frac{l\bullet \theta }{2}\) . \(\left(5\right)\)

figure 1(b) shows the local arrangement of the metasurface. The distance between the adjacent ring is 330 nm, and the distance between each nanopost in the ring is 330 nm, too. The toroidal arrangement of the metasurface could make the light field of Q-PV more even. Figure 1(c) shows the structure of a single nanopost with a height of H = 600 nm, the long axis is L = 250 nm, and the short axis is W = 95 nm. The TiO₂ elliptic nanopost is an artificial birefringent material with high refractive index at visible light wavelength. It can accurately control the emission phase of the partially transmitted light within 2π range [34]. For high dielectric index metasurfaces, due to a strong local waveguide effect, the interaction among adjacent nanoposts can be ignored. [34]. At the same time, the TiO₂ nanopost has a transmittance of up to 90% in the visible light band.

The generated Q-PVs by metasurfaces.

To generate Q-PV, we firstly chose metasurfaces with aperture diameter D of 20 µm, l of 1, 5, 10, f of 15 µm, r₀ of 5 µm to perform numerical simulations and made a comparison with the light field calculated by Fresnel integral diffraction theory. figures 2(a), (e), and (i) show the simulation results of the designed phase distributions of different TCs when l equals 1, 5, and 10, respectively. Phase values were obtained by extracting the transmission electrical field components and calculating with \(arctan\left[Im\right(Ex)/Re(Ex\left)\right]\). The simulated phase distributions agree with the theoretical ones. This proves that the designed metasurface can well accomplish phase control. figures 2(b), (f), and (j) respectively show the light field distributions of the Q-PV field in the xOz plane when l equals 1, 5, and 10. We can see that as the TC increases, the actual focal ring gradually deviates from the position of the designed focal length of f = 15 µm, which can be attributed to the Poynting vector’s oblique nature of the light field. However, the ring-shaped distribution in the designed focal plane remains unchanged. figures 2(c), (g), and (k) show the simulated light field in the focal plane with different l, their ring-shape distribution hardly changed with TC. The insets of figures 2(c), (g), and (k) are the phase distribution of the corresponding focal planes. The number of the mutation in the phase pattern corresponds to the number of TC. figures 2 (d), (h), and (i) show the theoretical light fields of Q-PVs calculated with equation 3, we can see that the simulation results were consistent with theoretical results. Therefore, the designed PB phase metasurface can successfully generate Q-PVs whose ring-shaped light field distributions hardly changed with TC.

As we know, the radius of the focal OV ring is positively correlated with TC. In order to study the divergence of the Q-PV, we quantitatively calculated the relationship between the measured radius (r) and TC by replacing the lens phase with the annular focusing phase. When designed radii (r₀) were different, as shown in figure 3(a). We took focal OV as a comparison, as shown by the black line, the slope of the fitted line is 0.904. When r₀ increases to 5 µm, we found that r still increased with TC, but the trend of the increase significantly reduced, appearing as the slope of the fitted straight line decreases to 0.299. When r₀ was set to 10 µm, the increasing slope of the r was further limited to 0.221. This kind of divergence was also found in reference [35] as the Q-PV has a lower slope for the relationship between TC and r when r₀ was larger. Meanwhile, when the focal length was set from 2 to 20 µm, we can see that r hardly changes with different focal lengths, as shown in figure 3(b). Therefore, the divergence characteristics of the Q-PV are not affected by the designed focal length.

Optical spanners with Q-PVs.

Particles can be trapped and rotated in the focal OV. Since Q-PV has fixed radius of the focal ring, we will explore the new possibility of Q-PV being used in the applications of optical tweezers and spanners. We firstly simulated the Q-PV to trap SiO₂ Mie-particles by choosing metasurfaces with aperture D of 20 µm, l of 1, 3, 5, f of 15 µm, r₀ of 5 µm. Meanwhile, the amplitude of the incident light was set as 300 hν (i.e., the power of the light source was set to 0.00093 mW). The position of the SiO₂ dielectric particle was set to the maximum intensity point of the focal ring (z = 14 µm), and along the x-axis from −10 µm to 10 µm, every 0.2 µm a step. We calculated the horizontal component of optical force exerted on the SiO₂ particles, as shown in figure 4(a). Since the force distribution has symmetry at x = 5 µm and −5 µm, x = −5 µm was taken as an example. The force at x = −6 µm shows that when the dielectric particle was close to the outer side of the focal ring, it firstly received an increasing optical force pointing to the center of the focal ring, and then the optical force gradually decreased to 0, which was just at the center of the focal ring of Q-PV (i.e., x = −5 µm). As the particle continue to move towards the inner side of the focal ring, it was pulled by increasing and decreasing optical pulling forces in turn, suggesting that the dielectric particle would be trapped in the center of the Q-PV focal ring. By using equation 11, we calculated the trapping potential for SiO₂ particle that integrated from x = −10 µm, as shown in figure 4(b). The trapping potential under each TC reached the lowest point of potential energy when the radius equals 5, and the depth of the potential wells reach below −200 k_BT. In general, when the depth of trapping potential reaches −1 k_BT, the particles can be stably trapped in the light field. Therefore, the dielectric particles can be stably trapped in the center of the orbit of the Q-PV.

Then we detected the OAM of Q-PV by detecting the angular optical force of metallic particles. The using of metallic particles is because they have high absorption and scattering that maximizes the transfer of OAM from the Q-PV [38]. Using Ag particles with a radius of 0.2 µm, we tested the angular force of the micro-sized particles every 30° on the focal ring, and the obtained results are shown in figure 4(c). As TC increases, the angular optical force exerts on the metallic particles gradually increases. The fluctuation of the force distribution was due to the uneven intensity of the light field, which was caused by the discrete distribution of nanoposts of the metasurface. We concluded that Q-PV not only can stably trap particles on fixed orbits, but also it can provide optical angular force that increases with TC for optical spanner use.

Q-PV arrays.

In order to design a single metasurface with multi-channel OAM states, coaxial and non-coaxial Q-PV arrays are generated by stacking the complex amplitudes of Q-PV that carrying different TCs. By superimposing complex amplitudes of different single Q-PVs, the summation complex function is calculated, and the phase profile of the Q-PV array is then determined.

By introducing a linear phase gradient along the θ direction onto the phase of Q-PV: \(\text{e}\text{x}\text{p}\left[ikrS\left(\text{s}\text{i}\text{n}\theta \bullet \text{s}\text{i}\text{n}{\theta }_{n}+\text{c}\text{o}\text{s}\theta \bullet \text{c}\text{o}\text{s}{\theta }_{n}\right)/f\right]\) [35], where (\(r,\theta\)) is the polar coordinate parameter of each sampling point on the metasurfaces, (S, θ_n) is the designed position of the Q-PV in polar coordinates in the focal plane, the Q-PV will propagate away from the \(z\)-direction at a certain diffraction angle. The single phase of Q-PV in the array could be expressed as

\({{\Phi }}_{n}=krS\left(\text{s}\text{i}\text{n}\theta \bullet \text{s}\text{i}\text{n}{\theta }_{n}+\text{c}\text{o}\text{s}\theta \bullet \text{c}\text{o}\text{s}{\theta }_{n}\right)/f+\left\{k\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)+{l}_{n}\bullet \theta \right\}\). \(\left(6\right)\)

where n = 1, 2, 3, 4. By using the phase packing formula, we finally got the combinate phase of the non-coaxial Q-PV array [39]:

\({{\Phi }}_{array}=\text{a}\text{r}\text{g}\left\{\sum _{\text{i}=1}^{n}exp\left(i{{\Phi }}_{\text{i}}\right)\right\}\). \(\left(7\right)\)

Note that the italic i in the formula is the imaginary unit.

As shown in figure 5, a non-coaxial Q-PV array composed of four single Q-PVs is demonstrated here. The size of metalens is 40 µm × 40 µm, focal length f = 50 µm, the displacement distance S_n = 10 µm, the parameters of every single Q-PV are listed as follows, Q-PV₁ (l₁ = 1, θ₁ = 90°), Q-PV₂ (l₂ = 3, θ₂ = 180°), Q-PV₃ (l₃ = 5, θ₃ = 0°) and Q-PV₄ (l₄ = 8, θ₄ = 270°). figure 5(a) shows the PB phase distribution that produced the Q-PV array. The intensity distribution of the Q-PV array at the propagation distance z = 50 µm is shown in figure 5(b). The results show that the focal ring strengths of four different TCs have the same radius. The center positions of the four Q-PVs are located at (10 µm, 90°), (10 µm, 180°), (10 µm, 0°), and (10 µm, 270°), respectively, expressed in polar coordinates, and are consistent with the designed position displacement (S, θ) of each Q-PV. The phase distribution in the focal plane is shown in figures 5(c) and (d). The blue rectangles box out the position of each Q-PV in figure 5(c), figure 5(d) shows the partially enlarged phase pattern. The interference fork fringe pattern can be observed near the Q-PV center, which is due to the fact that Q-PV exudes oblique upward at a certain angle after adding the linear phase gradient. The effect characteristics similar to isoclinic interference occur. The number of fork fringe reflects the TC of each Q-PV, which is consistent with the designed TC. The generated non-coaxial PV array has potential applications in the field of optical communication.

When several concentric rings that carrying different OAM information are combined to form a single metasurface array, the position of each OAM information will become certain, which is because the radius of the focal ring can maintain the designed ring radius \({r}_{0}\), especially when TC is in the range of 1 to 10. By combining equations 2 & 7, the phase formula of the coaxial Q-PV metasurface array is obtained, following as

\({{\Phi }}_{coaxial}=\text{a}\text{r}\text{g}\left(\sum _{\text{i}=1}^{n}exp\left(ik\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{\text{i}}\right)}^{2}}\right)+{l}_{\text{i}}\bullet \theta \right)\right)\). \(\left(8\right)\)

As is shown in figure 6, the size of metalens is 40 µm × 40 µm, focal length f = 15 µm, Δr = 2 µm, n = 7, the parameters of every single Q-PV are listed as follows, Q-PV₁ (r₁ = 3 µm, l₁ = 1), Q-PV₂ (r₂ = 5 µm, l₂ = 2), Q-PV₃ (r₃ = 7 µm, l₃ = 3), Q-PV₄ (r₄ = 9 µm, l₄ = 4), Q-PV₅ (r₅ = 11 µm, l₅ = 5), Q-PV₆ (r₆ = 13 µm, l₆ = 6), Q-PV₇ (r₇ = 15 µm, l₇ = 7). The phase distribution of the 7-ring coaxial array is shown in figure 6(a). The intensity distribution of Q-PV at z = 15 µm is shown in figure 6(b), and the measured radii of each Q-PV ring are r₁' = 2.8 µm, r₂' = 5 µm, r₃' = 7 µm, r₄' = 9 µm, r₅' = 11 µm, r₆' = 13 µm, r₇' = 15 µm. They are consistent with the designed radii. The corresponding phase distribution is shown in figure 6(c). In this case, the phase distributions at the position of each ring are less recognizable. figure 6(d) shows the extraction of the phase information of each Q-PV ring, and the OAM information of each ring can be clearly distinguished, which is consistent with the number of TCs at the designed position of each Q-PV ring. The generated coaxial PV array has potential applications in the fields of ultra-high-density information storage and photon computing.

In summary, as an alternative candidate to PV, Q-PV was produced based on PB phase metasurface, it can be regarded as an annular light field with vortex phase gradient. Compared with the PVs mentioned by the predecessors, the Q-PV was generated by using a metasurface with combined phases, which can more accurately and conveniently define the focal light field. With the increase of the designed radius, the Q-PV can better get rid of the influence of TC on the radius of the maximum light intensity ring, which is not affected by focal length as well. Through the calculation of optical forces, particles can be trapped on the focal plane and rotated around the Q-PV ring. The Q-PV can be used for optical spanners that provide different angular optical forces with different TCs. Non-coaxial and coaxial arrays with Q-PVs were successfully generated, demonstrating the possibilities in the applications of optical information coding, photonic computing, and ultra-high-density information storage.

Simulation method.

The full-wave simulation results shown in figures 2, 5, and 6 were performed using Lumerical FDTD solutions. The optical forces were calculated by the optical force calculation toolbox based on Maxwell stress tensor (MST) available in commercial software Lumerical FDTD solutions.

Calculation methods.

The force of micro-/nano- sized particles in the light field can be calculated by integrating the MST on the particle surface. The time-average force (including gradient force and scattering force) exerted on the particle can be expressed as

\(〈F〉=\int \left\{\frac{\epsilon }{2}Re\left[E\bullet n\right]{E}^{}-\frac{\epsilon }{4}\left(E\bullet {E}^{}\right)n+\frac{\mu }{2}Re\left[\mu \left(H\bullet n\right){H}^{}\right]-\frac{\mu }{4}\left(H\bullet {H}^{}\right)n\right\}ds\), \(\left(9\right)\)

where \(\epsilon\) and \(\mu\) are the relative permittivity and relative permeability of the medium among the particle under test, and \(n\) is the normal unit perpendicular to the integration area \(ds\). For further calculation of the trapping potential for a particle at ρ, it is necessary to integrate the optical force of particle at position from ∞ to ρ₀. The trapping potential can be calculated with integral of continuous or discretized forces as follows

\(\text{U}\left({{\rho }}_{0}\right)=-{\int }_{\infty }^{{{\rho }}_{0}}F\left(\rho \right)d\rho\). \(\left(10\right)\)

\(U\left({\rho }_{n}\right)=\sum _{j=1}^{n}F\left({\rho }_{j}\right){\Delta }\rho\), \(\left(11\right)\)

where ρ_n is the position of particle n on the x-axis in the focal plane, Δρ is the distance between adjacent measurement points (particle n and particle n+1), and F(ρ_j) is the force on the particle at the position j. Generally, the trapping potential of the optical trap can usually be unified in the units of k_BT, where k_B is Boltzmann’s constant, and T is normal temperature.

Data availability

The data that support the findings of this study are available within the paper.

Competing interests

The authors declare no competing interests.

Acknowledgements

This work was supported by National Natural Science Foundation of China (61805119), Natural Science Foundation of Jiangsu Province (BK20180469; BK20180468), Fundamental Research Funds for the Central Universities (30919011275). The authors wish to thank the anonymous reviewers for their valuable suggestions.

Author contributions

Z.S was responsible for the original research concept and the accurate physical interpretations. Z.Y.X was responsible for the numerical simulations and the manuscript draft. All authors reviewed the manuscript.

Albaladejo, S., Marques, M. I., Laroche, M. & Saenz, J. J. Scattering forces from the curl of the spin angular momentum of a light field. Phys Rev Lett, 102, 113602 (2009).
Ng, J., Lin, Z. & Chan, C. T. Theory of optical trapping by an optical vortex beam. Phys Rev Lett, 104, 103601 (2010).
Grier, D. G. A revolution in optical manipulation. Nature, 424, 810–816 (2003).
Gregg, P., Kristensen, P. & Ramachandran, S. Conservation of orbital angular momentum in air-core optical fibers.Optica2(2015).
Liu, E. et al. Robust transmission of orbital angular momentum mode based on a dual-cladding photonic quasi-crystal fiber.Journal of Physics D: Applied Physics52(2019).
Gibson, G. et al. Free-space information transfer using light beams carrying orbital angular momentum. Opt. Express, 12, 5448–5456 (2004).
Garces-Chavez, V. et al. Observation of the transfer of the local angular momentum density of a multiringed light beam to an optically trapped particle. Phys Rev Lett, 91, 093602 (2003).
Curtis, J. E. & Grier, D. G. Structure of optical vortices. Phys Rev Lett, 90, 133901 (2003).
Simpson, N. B., Dholakia, K., Allen, L. & Padgett, M. J. Mechanical equivalence of spin and orbital angular momentum of light: an optical spanner. Opt. Lett, 22, 52–54 (1997).
Ng, J., Lin, Z. & Chan, C. T. Theory of optical trapping by an optical vortex beam. Phys. Rev. Lett, 104, 103601 (2010).
Paterson et al. Controlled rotation of optically trapped microscopic particles.Science(2001).
MacDonald & M., P. Creation and manipulation of three-dimensional optically trapped structures., 296, 1101–1103 (2002).
Grier, D. G. A revolution in optical manipulation. Nature, 424, 810–816 (2003).
Shen, Z., Su, L., Yuan, X. C. & Shen, Y. C. Trapping and rotating of a metallic particle trimer with optical vortex.Applied Physics Letters109(2016).
Zhang, Y. et al. A plasmonic spanner for metal particle manipulation. Sci Rep, 5, 15446 (2015).
Allen, L., Beijersbergen, M. W., Spreeuw, R. J. C. & Woerdman, J. P. Orbital angular momentum of light and transformation of Laguerre Gaussian Laser modes. Phys. Rev. A, 45, 8185–8189 (1992).
Molina-Terriza, G., Torres, J. P. & Torner, L. Twisted photons. Nature Physics, 3, 305–310 (2007).
Bozinovic, N. et al. Terabit-Scale Orbital angular momentum mode division multiplexing in fibers., 340, 1545–1548 (2013).
Mehmood, M. Q. et al. Visible-frequency metasurface for structuring and spatially multiplexing optical vortices. Adv Mater, 28, 2533–2539 (2016).
Lei, T. et al. Massive individual orbital angular momentum channels for multiplexing enabled by Dammann gratings. Light Science & Applications, 4, e257 (2015).
Moreno, I., Davis, J. A., Cottrell, D. M., Zhang, N. & Yuan, X. C. Encoding generalized phase functions on Dammann gratings. Opt. Lett, 35, 1536 (2010).
Liu, Y. et al. Generation of perfect vortex and vector beams based on Pancharatnam-Berry phase elements. Sci Rep, 7, 44096 (2017).
Ostrovsky, A. S., Rickenstorff-Parrao, C. & Arrizón, V. Generation of the “perfect” optical vortex using a liquid-crystal spatial light modulator. Opt. Lett, 38, 534–536 (2013).
Garcia-Garcia, J., Rickenstorff-Parrao, C., Ramos-Garcia, R., Arrizon, V. & Ostrovsky, A. S. Simple technique for generating the perfect optical vortex. Opt Lett, 39, 5305–5308 (2014).
Ostrovsky, A. S., Rickenstorff-Parrao, C. & Arrizon, V. Generation of the "perfect" optical vortex using a liquid-crystal spatial light modulator. Opt. Lett, 38, 534–536 (2013).
García-García, J., Rickenstorff-Parrao, C., Ramos-García, R., Arrizón, V. & Ostrovsky, A. S. Simple technique for generating the perfect optical vortex. Opt. Lett, 39, 5305–5308 (2014).
Chen, M., Mazilu, M., Arita, Y., Wright, E. M. & Dholakia, K. Dynamics of microparticles trapped in a perfect vortex beam. Opt Lett, 38, 4919–4922 (2013).
Vaity, P. & Rusch, L. Perfect vortex beam: Fourier transformation of a Bessel beam. Opt Lett, 40, 597–600 (2015).
Xin, J. et al. Generation of Bessel beams from a diffractive ring lens. Opt. Commun, 310, 25–30 (2014).
Kildishev, A. V., Boltasseva, A. & Shalaev, V. M. Planar photonics with metasurfaces.Science, 1289(2013).
Yu, N. & Capasso, F. Flat optics with designer metasurfaces. Nature Materials, 13, 139–150 (2014).
Shen, Z. et al. Generation of optical vortices with polarization-insensitive metasurfaces. IEEE Photonics Journal, 12, 1–10 (2020).
Shen, Z., Xiang, Z. Y., Wang, Z. Y., Shen, Y. C. & Zhang, B. F. Optical spanner for nanoparticle rotation with focused optical vortex generated through a Pancharatnam-Berry phase metalens. Appl Optics, 60, 4820–4826 (2021).
Ma, Y., Rui, G., Gu, B. & Cui, Y. Trapping and manipulation of nanoparticles using multifocal optical vortex metalens. Sci Rep, 7, 14611 (2017).
Zhang, Y., Liu, W., Gao, J. & Yang, X. Generating focused 3D perfect vortex beams by plasmonic metasurfaces.Advanced Optical Materials6(2018).
Lin, R. & Li, X. Multifocal metalens based on multilayer Pancharatnam–Berry phase elements architecture.Optics Letters44(2019).
Chen, X. et al. Dual-polarity plasmonic metalens for visible light. Nat Commun, 3, 1198 (2012).
Zhao, Y., Edgar, J. S., Jeffries, G., McGloin, D. & Chiu, D. Spin-to-orbital angular momentum conversion in a strongly focused optical beam. Physical review letters, 99, 073901 (2007).
Yue, F. Y. et al. Multichannel polarization-controllable superpositions of orbital angular momentum states.Adv Mater29(2017).

No competing interests reported.

Download PDF

Editorial decision: Major revision
10 Dec, 2021
Reviews received at journal
09 Nov, 2021
Reviewers agreed at journal
09 Nov, 2021
Reviewers invited by journal
22 Oct, 2021
Editor assigned by journal
22 Oct, 2021
Editor invited by journal
22 Oct, 2021
Submission checks completed at journal
22 Oct, 2021
First submitted to journal
19 Oct, 2021

You are reading this latest preprint version

Quasi- Perfect Vortices Generated by Pancharatnam-Berry Phase Metasurfaces for Optical Spanners and OAM Communication

Status:

Version 1

Abstract

Figures

Introduction

Results

Theoretical description of Q-PV.

\(E\left(r,\theta ,f\right)=\delta \left(r-{r}_{0}\right)\bullet \text{e}\text{x}\text{p}\left(il\theta \right), \left(l=1, 2, 3\dots \right)\). \(\left(1\right)\)

\(T\left(r,\theta \right)=\text{e}\text{x}\text{p}\left(ik\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\right)\bullet \text{e}\text{x}\text{p}\left(il\theta \right)\), \(\left(2\right)\)

\(E\left(r,\rho ,z\right)=\frac{k}{z}{\int }_{0}^{{\rho }_{0}}\text{e}\text{x}\text{p}\left(ik\sqrt{{r}^{2}+{z}^{2}}-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\bullet {J}_{l}\left(\frac{k\rho r}{z}\right)rdr\), \(\left(3\right)\)

The metasurface design.

The generated Q-PVs by metasurfaces.

Optical spanners with Q-PVs.

Q-PV arrays.

\({{\Phi }}_{array}=\text{a}\text{r}\text{g}\left\{\sum _{\text{i}=1}^{n}exp\left(i{{\Phi }}_{\text{i}}\right)\right\}\). \(\left(7\right)\)

\({{\Phi }}_{coaxial}=\text{a}\text{r}\text{g}\left(\sum _{\text{i}=1}^{n}exp\left(ik\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{\text{i}}\right)}^{2}}\right)+{l}_{\text{i}}\bullet \theta \right)\right)\). \(\left(8\right)\)

Conclusions

Methods

Simulation method.

Calculation methods.

\(〈F〉=\int \left\{\frac{\epsilon }{2}Re\left[E\bullet n\right]{E}^{}-\frac{\epsilon }{4}\left(E\bullet {E}^{}\right)n+\frac{\mu }{2}Re\left[\mu \left(H\bullet n\right){H}^{}\right]-\frac{\mu }{4}\left(H\bullet {H}^{}\right)n\right\}ds\), \(\left(9\right)\)

\(\text{U}\left({{\rho }}_{0}\right)=-{\int }_{\infty }^{{{\rho }}_{0}}F\left(\rho \right)d\rho\). \(\left(10\right)\)

\(U\left({\rho }_{n}\right)=\sum _{j=1}^{n}F\left({\rho }_{j}\right){\Delta }\rho\), \(\left(11\right)\)

Declarations

Data availability

Competing interests

Acknowledgements

Author contributions

References

Additional Declarations

Status:

Version 1

Quasi- Perfect Vortices Generated by Pancharatnam-Berry Phase Metasurfaces for Optical Spanners and OAM Communication

Status:

Version 1

Abstract

Figures

Introduction

Results

Theoretical description of Q-PV.

\(E\left(r,\theta ,f\right)=\delta \left(r-{r}_{0}\right)\bullet \text{e}\text{x}\text{p}\left(il\theta \right), \left(l=1, 2, 3\dots \right)\). \(\left(1\right)\)

\(T\left(r,\theta \right)=\text{e}\text{x}\text{p}\left(ik\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\right)\bullet \text{e}\text{x}\text{p}\left(il\theta \right)\), \(\left(2\right)\)

\(E\left(r,\rho ,z\right)=\frac{k}{z}{\int }_{0}^{{\rho }_{0}}\text{e}\text{x}\text{p}\left(ik\sqrt{{r}^{2}+{z}^{2}}-\sqrt{{f}^{2}+{\left(r-{r}_{0}\right)}^{2}}\right)\bullet {J}_{l}\left(\frac{k\rho r}{z}\right)rdr\), \(\left(3\right)\)

The metasurface design.

The generated Q-PVs by metasurfaces.

Optical spanners with Q-PVs.

Q-PV arrays.

\({{\Phi }}_{array}=\text{a}\text{r}\text{g}\left\{\sum _{\text{i}=1}^{n}exp\left(i{{\Phi }}_{\text{i}}\right)\right\}\). \(\left(7\right)\)

\({{\Phi }}_{coaxial}=\text{a}\text{r}\text{g}\left(\sum _{\text{i}=1}^{n}exp\left(ik\left(f-\sqrt{{f}^{2}+{\left(r-{r}_{\text{i}}\right)}^{2}}\right)+{l}_{\text{i}}\bullet \theta \right)\right)\). \(\left(8\right)\)

Conclusions

Methods

Simulation method.

Calculation methods.

\(〈F〉=\int \left\{\frac{\epsilon }{2}Re\left[E\bullet n\right]{E}^{*}-\frac{\epsilon }{4}\left(E\bullet {E}^{*}\right)n+\frac{\mu }{2}Re\left[\mu \left(H\bullet n\right){H}^{*}\right]-\frac{\mu }{4}\left(H\bullet {H}^{*}\right)n\right\}ds\), \(\left(9\right)\)

\(\text{U}\left({{\rho }}_{0}\right)=-{\int }_{\infty }^{{{\rho }}_{0}}F\left(\rho \right)d\rho\). \(\left(10\right)\)

\(U\left({\rho }_{n}\right)=\sum _{j=1}^{n}F\left({\rho }_{j}\right){\Delta }\rho\), \(\left(11\right)\)

Declarations

Data availability

Competing interests

Acknowledgements

Author contributions

References

Additional Declarations

Status:

Version 1

\(〈F〉=\int \left\{\frac{\epsilon }{2}Re\left[E\bullet n\right]{E}^{}-\frac{\epsilon }{4}\left(E\bullet {E}^{}\right)n+\frac{\mu }{2}Re\left[\mu \left(H\bullet n\right){H}^{}\right]-\frac{\mu }{4}\left(H\bullet {H}^{}\right)n\right\}ds\), \(\left(9\right)\)