Abstract
Transverse spin, a spin component with unique characteristics, provides a new dimension for plenty of applications, such as optical trapping, imaging, and communication. Here, we analyze the pure transverse spin in the Bessel beam, which is solely present in the azimuthal direction. Based on a single layer dielectric metasurface, we efficiently generate Bessel beams with pure transverse spin in a compact optical system. As designed, the transverse spin is flexibly tunable by converting the polarization of the incident light. Furthermore, in the scattered Bessel beam, the local electromagnetic field oscillates around the transverse axis, which is perpendicular to the beam propagation. At certain positions, the local polarization ellipse degenerates into a perfect circle, resulting in a ring-periodic distribution of circularly polarized points (C points) in the beam. This suggests that the local polarization demonstrates a nontrivial periodic structure. This work deepens our understanding of spin-related physics and opens a new avenue for the study and application of transverse spins in ultracompact, flat, multifunctional nanophotonic platforms.
© 2024 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
Transverse spin (T spin), an interesting component of electromagnetic (EM) field spin, has attracted much attention in recent years [1–4]. In contrast to the conventional longitudinal spin, the helix-independent T spin is perpendicular to the direction of light propagation. In physics, the T spin is related to the inhomogeneity of momentum distribution in the EM field. The direction of the spin in the evanescent waves is locked by momentum, leading to the concept of optical spin-momentum locking [1,2,5–9]. For beams with complex momentum structures, the T spin originates from the longitudinal electric field, as explained by the theory of photonic wheels [10]. As a physically measurable quantity [11–13], T spin sparks research interest in a wide range of fascinating phenomena such as optical manipulation [9,11,13–15], imaging [16], and communications [17] due to its unique physical properties.
Bessel beams are typical nondiffractive beams, and have properties of self-healing and self-accelerating [18–21]. In the conventional methods, Bessel beams can be generated by utilizing diffractive optical elements or axicons with a high numerical aperture (NA) [22,23], which provides the possibility of applying T spin according to the theory of photonic wheels [10]. However, employing axicons result in a complex optical system with limited flexibility and substantial volume. Fortunately, metasurfaces have emerged as potential candidates for manipulating light in compact optical systems due to the outstanding advantages of their thinness, ease of integration, and design flexibility [24–26]. With the capability to flexibly manipulate the wavefront, polarization and dispersion [27,28], the metasurface allows for the generation of Bessel beams with high NA, which also implies that the metasurface can be a powerful platform for the observation and application of T spins. In addition, compared to conventional tightly focused beams (such as Gaussian beams [29]), the nondiffractive beams allow the pure T spin not restricted to the focal plane. In order to obtain the pure T spin in the whole space, we choose nondiffractive Bessel beams to study.
In this work, we analytically and numerically demonstrate the pure T spin existing in Bessel beams, which is oriented along the azimuthal direction and polarization-dependent. Based on the principle of the Pancharatnam-Berry (PB) phase [24,30,31], we designed a single layer dielectric metasurface, in which arrayed $\text {TiO}_2$ nanopillars provide the required geometrical phase distribution for the incident linearly polarized (LP) light and efficiently generate Bessel beams with pure T spin. According to the design, the different linear polarizations of the incident light affect the polarization of the scattered vector Bessel beam, which results in altering the ratio of the electric and magnetic spin components. To be specific, $x$ LP incident light generates a T spin with the pure electric component, and in contrast, $y$ LP incident light generates a pure magnetic T spin. In addition, the T spin allows the local EM field to oscillate around the transverse axis, which is perpendicular to the direction of propagation [29,32]. Upon analyzing the local polarization, we observe that the C point [29,33], symbolizing the transition of the polarization ellipse into a perfect circle, displays a periodic circular distribution within the scattered Bessel beam.
2. Method and structure description
To obtain the pure T spin in the beam, the electric field of the vector Bessel beam is analyzed, which is written in the cylindrical coordinate as [34–37],
where $(r, \phi, z)$ denote the cylindrical coordinates. The coefficients $c_1$ and $c_2$ describe the complex amplitude of the transverse electric (TE, with $c_1=1$ and $c_2=0$) and transverse magnetic (TM, with $c_1=0$ and $c_2=1$) wave components, respectively. The vector cylindrical wave functions $\mathbf {M}^{(c)}_n(r, \phi, z)$ and $\mathbf {N}^{(c)}_n(r, \phi, z)$ read [34,37,38],As a consequence, we utilize the principle of the PB phase to design the metasurface [24]. Since the adjustment of the meta-atom, the circularly polarized incident light is converted into the scattered cross-polarized field and carries a phase of $2\sigma \theta$ [39], where $\sigma$ characterizes the polarization state of the incident light, $\sigma =\pm 1$ corresponds to the right (RCP) and left (LCP) circularly polarized light, respectively. And $\theta$ is the rotation angle of the meta-atom. Actually, this additional phase is a geometric phase resulting from the evolution of the polarization state of light on the Poincare sphere. The phase imparted to the incident light from an axicon can be expressed as $\varphi _1=-kr\text {sin}\alpha$ [31,39,40], and we consider loading this phase for the meta-atom. Here, the numerical aperture of the axicon is $\text {NA}=\text {sin}\alpha$. For an LP incident light, the meta-atoms individual effect of its LCP and RCP components. As an example, for $x$ LP light $\begin {bmatrix}1 \\0\end {bmatrix}=\frac {1}{2}\left (\begin {bmatrix}1 \\i\end {bmatrix}+\begin {bmatrix}1 \\-i\end {bmatrix} \right )$, to obtain a scattered Bessel beam $\frac {1}{2}\left (\begin {bmatrix}1 \\-i\end {bmatrix}e^{ikr\text {sin}\alpha }+\begin {bmatrix}1 \\i\end {bmatrix}e^{-ikr\text {sin}\alpha } \right )=\begin {bmatrix}1 \\0\end {bmatrix}e^{-ikr\text {sin}\alpha }$, we expect that the RCP and LCP components of the LP incident light are compensated for the phase with a difference of $2m\pi$ [39], i.e., $2kr\text {sin}\alpha =2m\pi$, where $m$ is a positive integer. This implies that the nanopillars are arranged with certain conditions in the radial direction. Thereafter, in order to induce the electric field in the transverse direction as shown in Eq. (3), the phase $\varphi _2=\text {arctan}(y/x)$ is assigned to the meta-atom, which causes different circularly polarized light to acquire opposite topological charges [31,40], and it is the key to generating radial polarization. For $x$ LP incident light, the RCP and LCP components obtain the phase of $\pm \varphi _2$ respectively, and the scattered light can be written as $\frac {1}{2}\left (\begin {bmatrix} 1 \\ -i \end {bmatrix} e^{i \varphi _2}+\begin {bmatrix} 1 \\ i \end {bmatrix} e^{-i \varphi _2}\right ) e^{-i k r\text {sin}\alpha }=\begin {bmatrix} \text {cos}\phi \\ \text {sin}\phi \end {bmatrix}$, which is exactly the radially polarized light. In summary, the characterizations of the designed metasurface can be expressed as
Equation (5) characterizes the angle of rotation that the meta-atoms acquire and the positional conditions that should be satisfied to arrange them.
We designed a metasurface to produce the Bessel beam with a cone angle $\alpha =70^{\circ }$ and pure T spin when illuminated by LP incident light with a wavelength $\lambda =630$ nm. Figures 1(a) and 1(b) indicate that for $x$ and $y$ LP incidence, the metasurface generates pure electric spin $\mathbf {s}_E$ and magnetic spin $\mathbf {s}_H$, respectively. To establish the condition in Eq. (5), the meta-atoms are arranged in cylindrical coordinate $(r,\phi,z=0)$ and the interval between the meta-atoms in the radius direction is calculated to be $\delta r=335$ nm, with 40 arranged rings. In the azimuthal direction, the number of nanopillars arranged on the $m$th circle is 6$m$. As shown in Fig. 1(c), the metasurface consists of $\text {TiO}_2$ nanopillars and a glass substrate (with diameter $D=27$ $\mathrm{\mu}$m), and the relative refractive indexes are set as 2.40 and 1.45, respectively. To obtain high phase response and transmittance, the dimensions of the $\text {TiO}_2$ nanopillars are chosen to be length $W_x=330$ nm, width $W_y=90$ nm, and height $H=600$ nm (black dot in Figs. 1(d) and (e)). Figures 1(f) and (g) demonstrate the additional phases $\varphi _\text {RCP}$ and $\varphi _\text {LCP}$ produced by the metasurface for the RCP and LCP components, respectively.
3. Results and discussions
Firstly, we give the analytical calculations according to Eqs. (3) and (4). Figures 2(a) and 2(b) show the intensity distributions of the electric and magnetic fields in the Bessel beams, respectively. The operating wavelength and the cone angle are also fixed at $\lambda =630$ nm and $\alpha =70^{\circ }$. While the polarization parameters $c_1=0$, $c_2=1$ ($c_1=1$, $c_2=0$) correspond to the radially polarized electric field (magnetic field) in Fig. 2(a) (Fig. 2(b)). The phase difference between the longitudinal and radial fields (see insets of Fig. 2(a) and 2(b)) creates the azimuthal T spin, which is the pure T spin in Eq. (4). Figure 2(c) is the analytical calculation according to Eq. (4) on the $xy$ plane. The intensity of spin density $|\mathbf {s}|^2$, and its components $|s_x|^2$, $|s_y|^2$ and $|s_z|^2$ remain identical for both $c_1=0$, $c_2=1$, and $c_1=1$, $c_2=0$ in Eqs. (4b) and (4c), so we only have to give the intensity of $c_1=0$, $c_2=1$ in Fig. 2(c). The intensities of the transverse components of the spin density ($|s_x|^2$ and $|s_y|^2$) show a two-lobe structure, while the intensity of the longitudinal component ($|s_z|^2$) is strictly zero, indicating the presence of pure T spin in the beam. Simulations with the three dimensional finite difference time domain (FDTD) method (Lumerical Solutions) are performed to demonstrate the designed metasurface. The entire simulation is immersed in a vacuum and selected as a perfectly matched layer (PML) in the $x$, $y$ and $z$ directions, and the incident light is an LP plane wave propagating in the $z$ direction. The radially polarized ($c_1=0$, $c_2=1$) Bessel beam $\begin {bmatrix}\text {cos}\phi \\\text {sin}\phi \end {bmatrix} e^{-ikr\text {sin}\alpha }$ in Fig. 2(a) is induced by the $x$ LP incident plane wave $\begin {bmatrix}1 \\0\end {bmatrix}$. Figure 3 is the numerical results of the metasurface. The electric field intensities of the scattered Bessel beam are exhibited in Fig. 3(a), which correspond to the analytical result of Eq. (3), and the $y$ component of the electric field is neglected in the $xoz$ plane. Hereinafter, we fix the reference plane at $z = 3$ $\mathrm{\mu}$m. Figures 3(b) and 3(c) show the intensity distributions of the spin density, which are normalized by the maximum value of $|\mathbf {s}_E|^2$. As shown in Fig. 3(b), the intensity distributions of the electric spin density $|\mathbf {s}_{E}|^2$, and its components $|s_{E,x}|^2$, $|s_{E,y}|^2$ and $|s_{E,z}|^2$, which are similar to the analytical results in Fig. 2(c). Figure 3(c) shows the intensity distributions of the magnetic spin densities, which are much smaller than the electric spin density. Therefore, the pure T spin originating from electric fields is observed in scattered Bessel beams, which corresponds to the analytical predictions in Fig. 2(b).
Switch the incident light to $y$ LP light $\begin {bmatrix}0 \\1\end {bmatrix}=\frac {i}{2}\left (\begin {bmatrix}1 \\i\end {bmatrix}e^{i\pi }+\begin {bmatrix}1 \\-i\end {bmatrix}\right )$, the polarization of the scattered Bessel beam is changed into $c_1=1$, $c_2=0$, which leads the pure T spin originated from the magnetic spin component $\mathbf {s}_H$. This function, which tunes the polarization and spin of the Bessel beam by merely switching the incident plane wave, is impossible to achieve with traditional optical elements such as axicon. The effect of the metasurface on the $y$ LP light is $\frac {i}{2}\left (\begin {bmatrix} 1 \\ -i \end {bmatrix}e^{i(\pi +\varphi _2)}+\begin {bmatrix} 1 \\ i \end {bmatrix} e^{-i \varphi _2}\right ) e^{-i k r\text {sin}\alpha }=-\begin {bmatrix} -\text {sin}\phi \\ \text {cos}\phi \end {bmatrix} e^{-i k r\text {sin}\alpha }$, which indicates that the transverse component of electric field in the Bessel beam is azimuthally polarized according to Eq. (3), and the transverse magnetic field is radially polarized. The simulation results of $y$ LP incidence are shown in Fig. 4. Figure 4(a) indicates the magnetic field intensities of the scattered Bessel beam in the $xoz$ plane. And the magnetic field component $|H_y|^2$, which corresponds to the Eq. (3) for parameters $c_1=1$ and $c_2=0$, is negligible. The electric spin density $|\mathbf {s}_E|^2$, and the components $|s_{E,x}|^2$, $|s_{E,y}|^2$, $|s_{E,z}|^2$ are negligible in Fig. 4(b). Figure 4(c) illustrates the intensity distribution of the magnetic spin density, which is consistent with the analytical results in Fig. 2(c). Whether the incident light is $x$ or $y$ linear polarization, the designed metasurface can flexibly generate pure T spin as analytically predicted in the demonstration, which further validates that the metasurface can be a new platform for T spin applications.
Furthermore, it is worth to notice that the pure T spin in the scattered Bessel beam allows the oscillation axis of the local electric field to be perpendicular to the propagation direction [32,33], thus exhibiting a remarkable polarization structure. To illustrate this phenomenon, we use the concept of C point, which represents the local polarization ellipse of the EM field at that position degenerating into a perfect circle [32], where $\mathbf {s}_E$ is the normalized electric field spin density with the electric field energy density $w_E=\varepsilon (|E_x|^2+|E_y|^2+|E_z|^2)/4$ [5,8], and
with the ratios 0 and 1 (normalized by $1/\omega$) characterize the linear and circular polarization of the electric field, respectively. In Fig. 5(a), we obtain the intensity distribution of $|\mathbf {s}_E|/w_E$ in the Bessel beam, with the simulation of incident $x$ LP light. This density distribution exhibit a periodic circular pattern in the $xy$ plane. It indicates that the metasurface modulates an $x$ LP plane wave (with a homogeneous polarization) into a beam with an interesting polarization structure. The polarization ellipse clearly describes the oscillation of the EM field around the spin $\mathbf {s}$. As shown in Fig. 5(b), the polarization ellipse is defined by the vectors $\mathbf {A}$ (major axes) and $\mathbf {B}$ (minor axes), which satisfying $|\mathbf {A}|\geqslant |\mathbf {B}|$. The explicit formula for the elliptic axis read4. Conclusion
In conclusion, we successfully analyze and demonstrate the pure T spin along the azimuthal direction in a Bessel beam, which is related to the polarization parameter. Based on the principle of the PB phase, we have designed a single layer dielectric metasurface that can convert $x$ LP light into a Bessel beam with pure T spin. By switching the incident light to $y$ LP light, the transverse electric field of the scattered Bessel beam is changed from radial polarization to azimuthal polarization, and the ratio of the electric and magnetic spin components is tunable. Moreover, the local EM field oscillates around the transverse axis, creating the C points with the ring profile, and the local polarization also exhibits a periodic distribution associated with these C points. With the method of FDTD, we validated the function of the metasurface and compared it with the results of analytical calculations to ensure the reliability of the results. This work enhances our understanding of spin-related physics and opens a new route for spin-based applications.
Funding
Ministry of Science and Technology of the People's Republic of China (2022YFC2808203); National Natural Science Foundation of China (62375245).
Disclosures
The authors declare no conflicts of interest.
Data availability
All data in the main text are available from the corresponding author upon reasonable request.
References
1. K. Y. Bliokh, D. Smirnova, and F. Nori, “Quantum spin Hall effect of light,” Science 348(6242), 1448–1451 (2015). [CrossRef]
2. K. Y. Bliokh, F. J. Rodríguez-Fortu no, F. Nori, et al., “Spin-orbit interactions of light,” Nat. Photonics 9(12), 796–808 (2015). [CrossRef]
3. K. Y. Bliokh, H. Punzmann, H. Xia, et al., “Field theory spin and momentum in water waves,” Sci. Adv. 8(3), eabm1295 (2022). [CrossRef]
4. Y. Chen, Y. Li, B. Xu, et al., “Pure transverse spin perpendicular to wave propagation in multiple plane waves,” Phys. Rev. B 108(24), L241112 (2023). [CrossRef]
5. K. Y. Bliokh and F. Nori, “Transverse and longitudinal angular momenta of light,” Phys. Rep. 592, 1–38 (2015). [CrossRef]
6. T. Van Mechelen and Z. Jacob, “Universal spin-momentum locking of evanescent waves,” Optica 3(2), 118–126 (2016). [CrossRef]
7. P. Shi, L. Du, C. Li, et al., “Transverse spin dynamics in structured electromagnetic guided waves,” Proc. Natl. Acad. Sci. U. S. A. 118(6), e2018816118 (2021). [CrossRef]
8. P. Shi, L. Du, A. Yang, et al., “Dynamical and topological properties of the spin angular momenta in general electromagnetic fields,” Commun. Phys. 6(1), 283 (2023). [CrossRef]
9. X. Yu, Y. Li, B. Xu, et al., “Anomalous lateral optical force as a manifestation of the optical transverse spin,” Laser Photonics Rev. 17(10), 2300212 (2023). [CrossRef]
10. A. Aiello, P. Banzer, M. Neugebauer, et al., “From transverse angular momentum to photonic wheels,” Nat. Photonics 9(12), 789–795 (2015). [CrossRef]
11. K. Y. Bliokh, A. Y. Bekshaev, and F. Nori, “Extraordinary momentum and spin in evanescent waves,” Nat. Commun. 5(1), 3300 (2014). [CrossRef]
12. M. Neugebauer, T. Bauer, A. Aiello, et al., “Measuring the transverse spin density of light,” Phys. Rev. Lett. 114(6), 063901 (2015). [CrossRef]
13. M. Antognozzi, C. Bermingham, R. Harniman, et al., “Direct measurements of the extraordinary optical momentum and transverse spin-dependent force using a nano-cantilever,” Nat. Phys. 12(8), 731–735 (2016). [CrossRef]
14. G. Ha, H. Zheng, X. Yu, et al., “Polarization effect on optical manipulation in a three-beam optical lattice,” Phys. Rev. A 100(3), 033817 (2019). [CrossRef]
15. X. Yu, Q. Ye, H. Chen, et al., “Simple algorithm for partial wave expansion of plasmonic and evanescent fields,” Opt. Express 25(4), 4201–4215 (2017). [CrossRef]
16. G. Araneda, S. Walser, Y. Colombe, et al., “Wavelength-scale errors in optical localization due to spin-orbit coupling of light,” Nat. Phys. 15(1), 17–21 (2019). [CrossRef]
17. Z. Shao, J. Zhu, Y. Chen, et al., “Spin-orbit interaction of light induced by transverse spin angular momentum engineering,” Nat. Commun. 9(1), 926 (2018). [CrossRef]
18. O. C. Vicente and C. Caloz, “Bessel beams: a unified and extended perspective,” Optica 8(4), 451–457 (2021). [CrossRef]
19. J. Zhi, Y. Guo, B. Hu, et al., “Generation of polarization rotation function Bessel beams based on all-dielectric metasurfaces,” Opt. Commun. 550, 130014 (2024). [CrossRef]
20. J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58(15), 1499–1501 (1987). [CrossRef]
21. Z. Bouchal, J. Wagner, and M. Chlup, “Self-reconstruction of a distorted nondiffracting beam,” Opt. Commun. 151(4-6), 207–211 (1998). [CrossRef]
22. T. Grosjean and I. Gauthier, “Longitudinally polarized electric and magnetic optical nano-needles of ultra high lengths,” Opt. Commun. 294, 333–337 (2013). [CrossRef]
23. O. Brzobohatý, T. Čižmár, and P. Zemánek, “High quality quasi-bessel beam generated by round-tip axicon,” Opt. Express 16(17), 12688–12700 (2008). [CrossRef]
24. L. Zhang, S. Mei, K. Huang, et al., “Advances in full control of electromagnetic waves with metasurfaces,” Adv. Opt. Mater. 4(6), 818–833 (2016). [CrossRef]
25. Y. Wang, Y. Yuan, Y. Liu, et al., “Extreme diffraction management in phase-corrected gradient metasurface by fourier harmonic component engineering,” Laser Photonics Rev. 17(7), 2300152 (2023). [CrossRef]
26. J. Li, Y. Yuan, Q. Wu, et al., “Bi-isotropic huygens’ metasurface for polarization-insensitive cross-polarization conversion and wavefront manipulation,” IEEE Trans. Antennas Propag. 72(3), 2445–2454 (2024). [CrossRef]
27. J. Li, Y. Yuan, G. Yang, et al., “Hybrid dispersion engineering based on chiral metamirror,” Laser Photonics Rev. 17(3), 2200777 (2023). [CrossRef]
28. D. Zhang, X. Liu, Y. Sun, et al., “Dispersion engineering of spoof plasmonic metamaterials via interdigital capacitance structures,” Opt. Lett. 48(6), 1383–1386 (2023). [CrossRef]
29. T. Bauer, M. Neugebauer, G. Leuchs, et al., “Optical polarization Möbius strips and points of purely transverse spin density,” Phys. Rev. Lett. 117(1), 013601 (2016). [CrossRef]
30. M. V. Berry, “Quantal phase factors accompanying adiabatic changes,” Proc. R. Soc. Lond. A 392(1802), 45–57 (1984). [CrossRef]
31. J. Yang, R. Zhao, Y. Li, et al., “Generation of the Bessel beam of longitudinally varied polarization with dielectric metasurfaces,” Adv. Opt. Mater. 11(10), 2202896 (2023). [CrossRef]
32. J. F. Nye and J. V. Hajnal, “The wave structure of monochromatic electromagnetic radiation,” Proc. R. Soc. Lond. A 409(1836), 21–36 (1987). [CrossRef]
33. M. V. Berry, “Index formulae for singular lines of polarization,” J. Opt. A: Pure Appl. Opt. 6(7), 675–678 (2004). [CrossRef]
34. Z. Bouchal and M. Olivík, “Non-diffractive vector Bessel beams,” J. Mod. Opt. 42(8), 1555–1566 (1995). [CrossRef]
35. J. Chen, J. Ng, P. Wang, et al., “Analytical partial wave expansion of vector Bessel beam and its application to optical binding,” Opt. Lett. 35(10), 1674–1676 (2010). [CrossRef]
36. N. Wang, J. Chen, S. Liu, et al., “Dynamical and phase-diagram study on stable optical pulling force in Bessel beams,” Phys. Rev. A 87(6), 063812 (2013). [CrossRef]
37. G. Ha, H. Zheng, X. Yu, et al., “Analytically decomposing optical force on a spherical particle in Bessel beams into conservative and non-conservative parts,” J. Opt. Soc. Am. B 37(1), 67–73 (2020). [CrossRef]
38. X. Yu, Y. Jiang, H. Chen, et al., “Approach to fully decomposing an optical force into conservative and nonconservative components,” Phys. Rev. A 100(3), 033821 (2019). [CrossRef]
39. Q. Jiang, H. Xiang, and D. Han, “Generation of Bessel beams with tunable longitudinal electric and magnetic fields using an all-dielectric metasurface,” Opt. Lett. 48(4), 920–923 (2023). [CrossRef]
40. Q. Zhou, M. Liu, W. Zhu, et al., “Generation of perfect vortex beams by dielectric geometric metasurface for visible light,” Laser Photonics Rev. 15(12), 2100390 (2021). [CrossRef]