Photonics Research, Volume. 11, Issue 11, 1959(2023)

On-chip sorting of orbital angular momentum beams using Bloch surface wave structures

Nannan Li1、†, Qi Zou1、†, Yizhi Lan1, Yaqi Wang1, Jun Zhang1, Michael Somekh1,2,3,4、*, Changjun Min1,5、*, Fu Feng1,2,6、*, and Xiaocong Yuan1,2,7、*
Author Affiliations
  • 1Nanophotonics Research Center, Institute of Microscale Optoelectronics & State Key Laboratory of Radio Frequency Heterogeneous Integration, Shenzhen University, Shenzhen 518060, China
  • 2Research Center for Humanoid Sensing, Zhejiang Lab, Hangzhou 311100, China
  • 3Faculty of Engineering, University of Nottingham, Nottingham NG7 2RD, UK
  • 4e-mail: michael.somekh@nottingham.ac.uk
  • 5e-mail: cjmin@szu.edu.cn
  • 6e-mail: fufeng@zhejianglab.com
  • 7e-mail: xcyuan@szu.edu.cn
  • show less

    Owing to their unique optical properties and new degrees of freedom, orbital angular momentum (OAM) beams have been applied in various fields. Detection of the topological charges (TCs) of OAM beams is the key step for their applications. However, on-chip sorting of OAM beams with large TCs still remains a challenge. In this paper, Bloch surface wave (BSW) structures with five semi-ring shaped nanoslits are modeled. A spatial separation of 135 nm on the chip is obtained between two neighboring OAM states. OAM beams with TCs up to 35 can be successfully sorted by the BSW structures, which is much larger than that using metallic structures (only seven). BSW structures exhibit better OAM sorting performances than metallic structures. We systematically show how the lower attenuation of BSW structures leads to far superior separation ability compared to surface plasmons propagating on metallic structures. In addition, sorting of two OAM beams with different TCs simultaneously can be achieved in this way. Our results reveal that BSW structures should be an excellent solution for OAM sorting with large TCs, which is beneficial for applications in integrated on-chip devices and optical communications.

    1. INTRODUCTION

    Orbital angular momentum (OAM) beams possess a helical phase wavefront and carry OAM information [14]. There is a phase singularity at the center of OAM beams, giving rise to a hollow intensity distribution. The azimuthal variation of the electric field profiles of OAM beams can be expressed as exp(ilφ), where l and φ are the topological charge (TC) and the azimuthal angle, respectively [4]. TC is a crucial characteristic for OAM beams. OAM beams with different TCs are orthogonal to each other [5,6]. The orthogonal property makes it possible to use multiple OAM beams as signal channels to carry information in free space without crosstalk, which is beneficial for high-capacity data transmission and optical communication [711]. In addition, owing to their unique optical properties and new degrees of freedom, OAM beams have been applied in various fields, including optical tweezers [1214], quantum entanglement [1517], nonlinear optics [1820], nanofabrication and optical machining [2123], high-resolution imaging [2426], chemical and bio-molecule detection [2729], and astronomy [3032].

    Detection of the TCs of OAM beams, namely, OAM sorting, is a key means to investigate their properties. Sorting of OAM beams also exhibits technological importance. The traditional strategies of OAM sorting usually employ spatial light modulators (SLMs) or other phase masks to carry out phase transformation and reconstruct the target phase. These strategies include fork diffraction gratings [3336], liquid crystal q-plates [3739], Faraday rotation interferometric devices [40], and optical geometric transformation [4143]. Nevertheless, these traditional strategies require complex and bulky equipment, and the target beams need to be aligned accurately with the equipment. These drawbacks limit their applications. Recently, on-chip detection of OAM beams has emerged rapidly [44], since the detection chips can be highly compact, portable, and easy to integrate with other devices [45]. The previous reports describing on-chip sorting of OAM beams all employ metallic structures [4649]. However, surface plasmon polaritons (SPPs) using metallic structures usually suffer from severe Ohmic loss. The propagation length of SPPs is small, and this attenuation brings severe constraints in the geometry of the resulting structure. As a result, on-chip sorting of OAM beams with large TCs still remains a challenge.

    On the other hand, Bloch surface waves (BSWs) are surface electromagnetic waves that appear at the interface between a periodic multilayer of dielectric materials and the surrounding medium, which are considered as the dielectric analogue of SPPs [50,51]. BSWs can be excited under the incident light with both transverse electric (TE) and transverse magnetic (TM) polarizations [51,52]. Since BSW structures consist of a periodic multilayer of dielectric materials, absorption loss of BSWs is much lower in comparison with that of SPPs [5052]. The low-loss property endows BSWs with longer propagation lengths and sharper resonance peaks, and BSW structures have been widely applied in high-sensitivity sensing [5154], high-contrast imaging [55], metallic particle trapping [56], and diffraction-free waveguides [5759]. BSW structures therefore provide a promising solution for on-chip sorting of OAM beams with large TCs, which has not been explored before.

    In this paper, we numerically investigate the on-chip sorting of OAM beams using BSW structures. The designed BSW structures consist of three dielectric layers, and five semi-ring shaped nanoslits were introduced on the top surface of the BSW structures. BSWs can be excited and modulated by the helical phase wavefront of the incident OAM beams. There is a transverse shift for the focal point owing to the constructive interference of the BSW. A spatial interval of 135 nm on the chip between two neighboring OAM states can be obtained by using BSW structures. Furthermore, OAM sorting behavior of metallic structures with the same morphology is also analyzed. OAM beams with TCs up to 35 can be successfully sorted by the BSW structures, which is much larger than that using metallic structures (only seven). BSW structures exhibit superior OAM sorting performances than SPPs due to the lower attenuation, which is beneficial for OAM sorting with large TCs. In addition, sorting of two OAM beams with different TCs simultaneously can be achieved in this way as long as the difference between these two TCs is no less than four. To our knowledge, this is the first demonstration of on-chip sorting of OAM beams with TCs up to 35 as well as the first report of on-chip sorting of two OAM beams simultaneously. Our results reveal that BSW structures could be an excellent solution for OAM sorting with large TCs, which is beneficial for applications in integrated on-chip devices and optical communication.

    2. PRINCIPLE AND MODELING

    The designed structure for exciting BSWs is schematically sketched in Fig. 1(a). The designed BSW structure consists of three layers, which are Si (top), SiO2 (middle), and Si substrate (bottom). The refractive indexes of Si and SiO2 were set as nh=3.48+0.001i and nl=1.47+0.005i [52], respectively. The thicknesses of the top Si layer and SiO2 layer were set as dh=0.12  μm and dl=1.88  μm, respectively. This structure was illuminated by incident light from the bottom. Finite element methods were employed to numerically calculate the absorption and electric field distribution of this structure. Boundary conditions with perfectly matched layers were utilized on all sides in the simulations. Under incident light with TE polarization, the dispersion relation of the BSW structure described by absorption as a function of incidence angle is shown in Fig. 1(b). The white dashed line indicates the wavelength of 1.55 µm, which is chosen as the working wavelength in this study. The absorption spectrum of the structure as a function of incidence angle at 1.55 µm under TE polarization is also plotted [top subfigure in Fig. 1(b)]. It can be observed that the absorption reaches a maximum value of 0.985 at 24.98° at 1.55 µm. Figure 1(c) shows the electric field distribution of this structure under an incidence angle of 24.98° at 1.55 µm under TE polarization. The electric field reaches maximum at the top surface of the structure and decays gradually in air. The real part of the refractive index of the structure is also plotted for comparison [top subfigure in Fig. 1(c)]. These results show that BSWs are successfully excited on the top surface of this structure.

    Designed structure for exciting BSWs. (a) Schematic of the structure composed of three layers: Si (top), SiO2 (middle), and Si substrate (bottom). (b) Dispersion relation of the structure described by absorption as a function of incidence angle under TE polarization. The white dashed line indicates the wavelength of 1.55 µm. The absorption spectrum of the structure as a function of incidence angle at 1.55 µm is plotted in the top subfigure. (c) Electric field distribution of the structure at incidence angle of 24.98° at 1.55 µm under TE polarization. The white curve refers to the electric field. The real part of the refractive index of the structure is plotted in the top subfigure. (d) Schematic of the mechanism of the OAM sorting methods. Only one nanoslit is considered. The phase along the semi-ring nanoslit has a linear dependence on the azimuthal angle φ. (e) Schematic of the designed BSW structure with five semi-ring shaped nanoslits. (f) Schematic of the geometric parameters of the BSW structure.

    Figure 1.Designed structure for exciting BSWs. (a) Schematic of the structure composed of three layers: Si (top), SiO2 (middle), and Si substrate (bottom). (b) Dispersion relation of the structure described by absorption as a function of incidence angle under TE polarization. The white dashed line indicates the wavelength of 1.55 µm. The absorption spectrum of the structure as a function of incidence angle at 1.55 µm is plotted in the top subfigure. (c) Electric field distribution of the structure at incidence angle of 24.98° at 1.55 µm under TE polarization. The white curve refers to the electric field. The real part of the refractive index of the structure is plotted in the top subfigure. (d) Schematic of the mechanism of the OAM sorting methods. Only one nanoslit is considered. The phase along the semi-ring nanoslit has a linear dependence on the azimuthal angle φ. (e) Schematic of the designed BSW structure with five semi-ring shaped nanoslits. (f) Schematic of the geometric parameters of the BSW structure.

    In order to perform on-chip sorting of OAM beams, several semi-ring shaped nanoslits were introduced on the top surface of the BSW structure. The mechanism of the OAM sorting methods is schematically illustrated in Fig. 1(d). We first consider the case of only one nanoslit. Considering the circularly polarized incident light with a phase profile of exp(ilφ), the electric field at the surface can be expressed as [46] E(r,φ)=0π|r0r|12exp[i(l+s)φ]exp(ikBSW|r0r|)dφ,where r is the position vector of the observing point P. r0 is the position vector of an arbitrary point on the semi-ring nanoslit. l denotes the TC of the incident OAM beams. s refers to the spin angular momentum (SAM) of the circularly polarized incident light. φ is the azimuthal angle. kBSW refers to the wavevector of the BSW. Due to the helical phase of the OAM beams, the phase along the semi-ring nanoslit has a linear dependence on the azimuthal angle φ. In this case, the excited BSWs do not focus at the original point O, and there is a transverse shift for the focal region. The details of the derivation of Eq. (1) are introduced in Appendix A. The transverse peak position of the focal point can be derived as xj=sgn(j)μj2πλBSW,where j=l+s, and j is the total angular momentum. μj is the first non-zero coefficient of Jj(x) [the derivative of Jj(x)]. λBSW is the effective wavelength of BSW. According to Eq. (2), the transverse peak position of the focal point is only determined by j and λBSW. s equals one for right-handed circular polarization (RCP), while s equals 1 for left-handed circular polarization (LCP). For multiple semi-ring nanoslits with different radii, the transverse position of the focal point remains the same since Eq. (2) is independent of the radius of the nanoslit. In addition, the above analysis of OAM sorting mechanism is also valid for radially polarized and linearly polarized incident light (see Appendix A).

    The designed structure for OAM sorting is schematically sketched in Fig. 1(e). The number of the nanoslits was chosen as five (as explained in the following section). The geometric parameters of the BSW structure are schematically shown in Fig. 1(f). The width of the nanoslits w was set as 0.6 µm, and the depth of the nanoslits was set as 0.2 µm. The geometry of the nanoslits was optimized to make the electric field at the focal point as large as possible. The period of the nanoslits s (the radial distance between two neighboring slits) was set to be λBSW to enhance the electric field intensity at the focal point. λBSW can be calculated as [52] λBSW=λ0neff=λ0nsinθ,where λ0 is the wavelength of the incident light (1.55 µm), neff is the effective refractive index, n is the refractive index of the Si substrate (3.48), and θ is the incidence angle where BSW is excited (24.98°). λBSW is calculated to be 1.054 µm. In this study, the period of the nanoslits d was set as 1.05 µm. The excitation angle of the TM mode of this dielectric structure is different from that of TE mode, so the effective wavelength of the TM mode is different. The five nanoslits with a period of 1.05 µm can serve as optical gratings, which selectively excite TE mode.

    3. RESULTS AND DISCUSSION

    A schematic of the designed BSW structure with five semi-ring shaped nanoslits is shown in Fig. 2(a). The radius of the outermost semi-ring nanoslit was set to 12.2 µm, and the radius of the innermost nanoslit was set to 8 µm. This structure was illuminated vertically by OAM beams with circular polarization from the bottom, and the power distribution of the OAM beams is uniform. Circularly polarized OAM beams are employed as incident light in this study, because they carry uniformly distributed electric fields along the azimuthal direction, and they are more stable compared with radially polarized OAM beams [46]. Figure 2(b) shows the electric field intensity (|E|2) distribution of the BSW structure under an incident OAM beam with TC of one. A bright focal point appears around the center with a transverse shift. We focus on the central square region with an edge length of 3 µm [surrounded by a white dashed line in Fig. 2(b)]. Spatial phase patterns of OAM beams with TCs of two, four, six, eight, and ten (from top to bottom) are plotted in Fig. 2(c) (left), and the electric field intensity of the BSW structures under the corresponding incident OAM beams are shown in Fig. 2(c) (right). The white dashed lines indicate the positions where Y=0. Figure 2(d) shows the normalized electric field intensity along the white dashed lines under incident OAM beams with TCs of one to ten. It clearly shows that the incident OAM beams with larger TCs can lead to a larger transverse shift of the peak positions of the focal points, which is consistent with our former analysis. Figure 2(e) shows the peak positions of the focal points as a function of the TCs of incident OAM beams. A linear fitting is performed, and R2 for this linear fitting is 0.999. The result reveals that a spatial interval of 135 nm on the chip between two neighboring OAM states can be obtained in this method.

    On-chip sorting of OAM beams with different TCs using BSW structures. (a) Schematic of the designed BSW structure with five semi-ring shaped nanoslits. The structure was illuminated vertically by the OAM beam with circular polarization from the bottom. (b) Electric field intensity (|E|2) of the BSW structure under incident OAM beam with TC of one. The edge length of the square region surrounded by the white dashed line is 3 µm. (c) Spatial phase patterns of OAM beams with TCs of two, four, six, eight, and ten (from top to bottom) (left) and the electric field intensity of the BSW structures under the corresponding incident OAM beams (right). All the subfigures share the same scale bar. The white dashed line indicates the position where Y=0. The green dashed line indicates the peak position of the focal points. (d) Normalized electric field intensity along the white dashed lines in (c) under incident OAM beams with TCs of one to ten. (e) Linear fitting of the peak positions of the focal points as a function of the TCs of incident OAM beams. R2 for this linear fitting is 0.999.

    Figure 2.On-chip sorting of OAM beams with different TCs using BSW structures. (a) Schematic of the designed BSW structure with five semi-ring shaped nanoslits. The structure was illuminated vertically by the OAM beam with circular polarization from the bottom. (b) Electric field intensity (|E|2) of the BSW structure under incident OAM beam with TC of one. The edge length of the square region surrounded by the white dashed line is 3 µm. (c) Spatial phase patterns of OAM beams with TCs of two, four, six, eight, and ten (from top to bottom) (left) and the electric field intensity of the BSW structures under the corresponding incident OAM beams (right). All the subfigures share the same scale bar. The white dashed line indicates the position where Y=0. The green dashed line indicates the peak position of the focal points. (d) Normalized electric field intensity along the white dashed lines in (c) under incident OAM beams with TCs of one to ten. (e) Linear fitting of the peak positions of the focal points as a function of the TCs of incident OAM beams. R2 for this linear fitting is 0.999.

    As a comparison, on-chip sorting of OAM beams using a metallic structure is also investigated. The metallic structure for exciting SPPs is schematically sketched in Fig. 3(a). This structure consists of an Au film (top) and a SiO2 substrate (bottom). The dielectric function of Au at 1.55 µm was set as εAu=116.62+11.46i. The refractive index of SiO2 at 1.55 µm was set as 1.47+0.005i. The thickness of Au film was set as 0.2 µm. Similarly, five semi-ring shaped nanoslits were introduced on the top surface of this structure. The depth of the nanoslits was set as 0.2 µm, and the width of the nanoslits was set as 0.6 µm. The period of the nanoslits (the radial distance between two neighboring slits) was set to be the effective wavelength of the SPPs. The effective wavelength of the SPPs can be calculated as [60] λSPP=Re[2π/(ωc0εdεmεd+εm)],where ω is the angular frequency of the incident light, c0 is the speed of light in free space, and εd and εm are the dielectric function of air and the Au film, respectively. λSPP is calculated to be 1.543 µm. In this study, the period of the nanoslits on the metallic structure was set to 1.54 µm. The radius of the outermost semi-ring nanoslit on the metallic structure was set as 12.2 µm (the same as that on the BSW structure), and the radius of the innermost nanoslit was set as 6.04 µm. This structure was illuminated vertically by OAM beams with circular polarization from the bottom. Electric field intensity (|E|2) of the metallic structure under an incident OAM beam with TC of one is plotted in Fig. 3(b). Similar to the case of the BSW structure, there is a focal point around the center with a transverse shift. The difference is that the electric field intensity at the focal point is relatively low. Spatial phase patterns of OAM beams with TCs of one, two, three, four, and five (from top to bottom) are plotted in Fig. 3(c) (left), and the electric field intensity of the metallic structure under the corresponding incident OAM beams is shown in Fig. 3(c) (right). The white dashed lines indicate the positions where Y=0. Figure 3(d) shows the normalized electric field intensity along the white dashed lines under incident OAM beams with TCs of one to five. It shows that as the TC of the incident OAM beams becomes larger, the peak position of the focal point shifts to the left gradually. The peak positions of the focal points as a function of the TCs of incident OAM beams are plotted in Fig. 3(e). A linear fitting is performed, and R2 for this linear fitting is 0.997. A spatial interval of 270 nm on the chip between two neighboring OAM states can be obtained using this metallic structure.

    On-chip sorting of OAM beams with different TCs using metallic structures. (a) Schematic of the metallic structure with five semi-ring nanoslits. This structure consists of an Au film (top) and a SiO2 substrate (bottom). The structure was illuminated vertically by the OAM beam with circular polarization from the bottom. (b) Electric field intensity (|E|2) of the metallic structure under incident OAM beam with TC of one. The edge length of the square region surrounded by the white dashed line is 3 µm. (c) Spatial phase patterns of OAM beams with TCs of one, two, three, four, and five (from top to bottom) (left) and electric field intensity of the metallic structure under the corresponding incident OAM beams (right). All the subfigures share the same scale bar. The white dashed line indicates the position where Y=0. The green dashed line indicates the peak position of the focal points. (d) Normalized electric field intensity along the white dashed lines in (c) under incident OAM beams with TCs of one to five. (e) Linear fitting of the peak positions of the focal points as a function of the TCs of incident OAM beams. R2 for this linear fitting is 0.997.

    Figure 3.On-chip sorting of OAM beams with different TCs using metallic structures. (a) Schematic of the metallic structure with five semi-ring nanoslits. This structure consists of an Au film (top) and a SiO2 substrate (bottom). The structure was illuminated vertically by the OAM beam with circular polarization from the bottom. (b) Electric field intensity (|E|2) of the metallic structure under incident OAM beam with TC of one. The edge length of the square region surrounded by the white dashed line is 3 µm. (c) Spatial phase patterns of OAM beams with TCs of one, two, three, four, and five (from top to bottom) (left) and electric field intensity of the metallic structure under the corresponding incident OAM beams (right). All the subfigures share the same scale bar. The white dashed line indicates the position where Y=0. The green dashed line indicates the peak position of the focal points. (d) Normalized electric field intensity along the white dashed lines in (c) under incident OAM beams with TCs of one to five. (e) Linear fitting of the peak positions of the focal points as a function of the TCs of incident OAM beams. R2 for this linear fitting is 0.997.

    The detection limits of OAM beams using metallic and BSW structures are carefully investigated. Figure 4(a) shows the electric field intensity of the metallic structure under an incident OAM beam with TCs of one, seven, and eight. We focus on the central line where Y=0 since the peak position of focal points lies in this line (green dashed line). The electric field intensity of metallic structure under an incident OAM beam with other TCs (two to twelve) is shown in Appendix B. When the metallic structure is illuminated by an OAM beam with TC not larger than seven, the electric field intensity at the focal point is the largest peak. When the TC of an incident OAM beam is greater than seven, however, the electric field intensity at the focal point is not the largest peak. In these cases, we cannot distinguish the peak positions of the focal points by simply examining the largest peaks of the electric field intensity along the central line. In this study, we set the detection criterion as follows: the electric field intensity at the focal point is the largest peak. Therefore, the detection limit of OAM beams using a metallic structure is seven. As a comparison, Fig. 4(b) shows the electric field intensity of the BSW structure under an incident OAM beam with TCs of 1, 35, and 36. The electric field intensity of the BSW structure under an incident OAM beam with other TCs (10–40) is also presented in Appendix B. When the BSW structure is illuminated by an OAM beam with TC not larger than 35, the electric field intensity at the focal point is the largest peak. When the TC of an incident OAM beam is greater than 35, however, the electric field intensity at the focal point is not the largest peak. Using the same detection criterion, we find the detection limit of OAM beams using the BSW structure is up to 35, which is much larger than that of the metallic structure (only seven).

    Comparison of the electric field intensity (|E|2) of the metallic structure and the BSW structure under incident OAM beams. (a) Electric field intensity of the metallic structure under incident OAM beams with TCs of one, seven, and eight. Subfigures in (a) share the same color scale bar. (b) Electric field intensity of the BSW structure under incident OAM beams with TCs of 1, 35, and 36. Subfigures in (b) share the same color scale bar. The green dashed lines in (a) and (b) indicate the position where Y=0. The white curves in (a) and (b) refer to the electric field intensity along the green dashed line. (c) Electric field intensity of the focal points on metallic structure (SPP, orange) and BSW structure (cyan) under incident OAM beams with TCs of one to seven. Electric field intensity of the focal points on BSW structure over that on metallic structure is also plotted. (d) Electric field intensity of the focal points on metallic structure (SPP, black) under incident OAM beams with TCs of one to seven and electric field intensity of the focal points on BSW structure (red) under incident OAM beams with TCs of 1 to 35.

    Figure 4.Comparison of the electric field intensity (|E|2) of the metallic structure and the BSW structure under incident OAM beams. (a) Electric field intensity of the metallic structure under incident OAM beams with TCs of one, seven, and eight. Subfigures in (a) share the same color scale bar. (b) Electric field intensity of the BSW structure under incident OAM beams with TCs of 1, 35, and 36. Subfigures in (b) share the same color scale bar. The green dashed lines in (a) and (b) indicate the position where Y=0. The white curves in (a) and (b) refer to the electric field intensity along the green dashed line. (c) Electric field intensity of the focal points on metallic structure (SPP, orange) and BSW structure (cyan) under incident OAM beams with TCs of one to seven. Electric field intensity of the focal points on BSW structure over that on metallic structure is also plotted. (d) Electric field intensity of the focal points on metallic structure (SPP, black) under incident OAM beams with TCs of one to seven and electric field intensity of the focal points on BSW structure (red) under incident OAM beams with TCs of 1 to 35.

    The reason for the large detection limit of OAM beams using the BSW structure is further explored. The electric field intensity of the metallic structure and the BSW structure is quantitatively calculated by averaging the intensity at a small region (0.1  μm×0.1  μm) around the peak position of the focal point. As Fig. 4(c) shows, the electric field intensity of the focal points on the BSW structure is much higher than that on the metallic structure under an incident OAM beam with the same TC. For incident OAM beams with TCs of one to seven, the ratios of electric field intensity of the BSW structure to that of the metallic structure are 8.17, 8.22, 8.71, 8.88, 8.68, 9.29, and 9.78, respectively. The electric field intensity of metallic and BSW structures along the central line where Y=0 under an incident OAM beam with different TCs is presented in detail in Appendix B. The electric field intensity of the focal points of metallic and BSW structures is summarized in Fig. 4(d). We can see that the electric field intensity of the focal points decreases with the increase of the TC of the incident OAM beam. For the BSW structure, the electric field intensity of the focal points decreases rapidly when the TC of the incident OAM beam is small (<10), while it decreases slowly when the TC of the incident OAM beam gets large (>20).

    The electric field intensity of focal points on the BSW structure is much higher than that of the metallic structure (over eight times). The reason for this is explained below. Since the BSW structures consist of several dielectric layers, absorption loss of BSWs is much lower in comparison with that of SPPs [5052]. To examine the contribution of the absorption loss, a low-loss artificial “ideal” material is modeled. The dielectric function of this artificial material is set as εa=116.62+0.01i. As a comparison, the dielectric function of Au at 1.55 µm is εAu=116.62+11.46i. The effective wavelength of this artificial material is 1.54 µm. The electric field intensity of the artificial structure with five nanoslits under different incident OAM beams is presented in Appendix C. The radius of the outermost nanoslit is 12.2 µm. Using the same detection criterion, the detection limit of OAM beams using this artificial structure is 13. When the focal point approaches the innermost nanoslit, the energy transfers to the nanoslit and the detection limit is severely restricted. It appears that the nanoslit with a relatively small radius limits the number of TCs that can be detected, but it is appropriate to increase the radius of the nanoslits since the attenuation of the “ideal” material is so low. The size of this structure is further increased, and its electric field intensity under different incident OAM beams is also shown in Appendix C. The radius of the outermost nanoslit is 18.36 µm, and the detection limit of OAM beams increases dramatically to 26. Since the real part of the dielectric function of this artificial material is the same as that of the Au, this result demonstrates that the low-absorption-loss property of the material is the major contribution to the larger detection limit. The reason for the difference between the “ideal” material and the BSW lies in the difference in the effective wavelengths of the different wavemodes. The effective wavelength of BSWs (1.05 µm) is smaller than that of SPPs (1.54 µm). The smaller effective wavelength results in a stronger focusing ability (the ability to focus energy to a small spot), which causes a smaller spot size and higher energy density at the spot size [61,62]. Taken together, the low absorption loss and smaller effective wavelength are the reasons for the larger detection limit of OAM beams using the BSW structure (35) than that using the metallic structure (7).

    In addition, the full width at half-maximum (FWHM) of electric field intensity of the metallic structure and BSW structure under incident OAM beams with different TCs is summarized in Table 1. FWHM of electric field intensity of the BSW structure is much smaller (approximately one half) than that of the metallic structure. The reason can be ascribed to the smaller effective wavelength of BSWs (1.05 µm) than that of SPPs (1.54 µm). The smaller effective wavelength results in a stronger focusing ability, which causes a smaller spot size and smaller FWHM of electric field intensity [61,62]. FWHM of electric field intensity of the BSW structure does not get broadened with the increase of the TC of the OAM beam. Although the FWHM of electric field intensity does not affect the detection limit directly, it can provide higher resolution for the detection. The improvement of detection resolution is essential for optical multiplexing systems that exploit the orthogonality of the OAM modes, especially for sorting multiple OAM beams simultaneously [43]. Taken together, BSW structures exhibit better OAM sorting performances than SPPs in terms of higher electric field intensity and smaller FWHM, which is beneficial for OAM sorting with large TCs. BSW structures are therefore a suitable and promising alternative for OAM sorting with large TCs.

    Full Width at Half-Maximum of Electric Field Intensity Distribution of the Metallic Structure and BSW Structure under Incident OAM Beams with Different TCs

    TC12345678910
    FWHM (SPP) (μm)0.780.780.780.800.800.820.82
    FWHM (BSW) (μm)0.420.420.420.420.420.420.400.400.400.40

    The above discussion is based on the condition that only one OAM beam is illuminated on the BSW structure. On the other hand, sorting of two or several OAM beams simultaneously is of vital importance to OAM demultiplexing and optical communication. To our knowledge, however, the demonstration of on-chip sorting of two OAM beams simultaneously has never been reported. Electric field intensity of the BSW structure under an incident OAM beam with TC of 3 is presented in Fig. 5(a) for comparison. Figure 5(b) shows the electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 1 and 3. Two focal points can be barely observed in this case. The positions of the focal points are quantitatively analyzed in Figs. 5(c)–5(e). For incident OAM beams with TCs of 1 and 3, the positions of these two focal points are X=0.24  μm and 0.3 µm [Fig. 5(c)], respectively. There is a dip between these two points. Since the position of the first focal point is beyond the first-order dark ring of the second focal point, these two focal points can be barely separated according to the Rayleigh criterion. Therefore, to separate two OAM beams, a difference of at least four for TCs is required. For incident OAM beams with TCs of 3 and 3, the intensity of the dip between the two focal points approaches zero [Fig. 5(d)], so these two focal points can be completely separated. For incident OAM beams with TCs of 5 and 3, the dip between the intensity of these two focal points is rather broad. Both the positions and intensities of these two focal points are approximately the same as the positions and intensities of the focal points generated by a single OAM beam [Fig. 5(e)]. In all these cases, the positions of these two focal points coincide with that of the focal point generated by a single OAM beam. Taken together, on-chip sorting of two OAM beams simultaneously using BSW structures can be achieved as long as the difference between these two TCs is no less than four, and the positions of these two focal points are the same as that of the focal points generated by a single OAM beam. By examining the positions of the two focal points, the information of the two incident OAM beams can be easily recovered. The above discussion is based on two incident beams with the same polarization state (circular polarization). However, if the polarization states of the two incident beams are different from each other (linear polarization and circular polarization), the sorting of the two OAM beams would still work, but the sorting performance is worse than sorting two OAM beams with the same polarization state (circular polarization).

    On-chip sorting of two OAM beams simultaneously. (a) Electric field intensity of the BSW structure under the incident OAM beam with TC of −3. The edge length of the square region surrounded by the white dashed line is 2 µm. (b) Spatial phase patterns of OAM beams with TCs of 1 and −3 (from top to bottom) (left) and electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 1 and −3 (right). The edge length of the square region surrounded by the white dashed line is 2 µm. (c) Electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 1 and −3 (top), and the electric field intensity along the green dashed line (bottom). (d) Electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 3 and −3 (top), and the electric field intensity along the green dashed line (bottom). (e) Electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 5 and −3 (top), and the electric field intensity along the green dashed line (bottom); (c)–(e) share the same color scale bar. The green dashed line indicates the position where Y=0. The purple dashed lines indicate the peak positions of the two focal points.

    Figure 5.On-chip sorting of two OAM beams simultaneously. (a) Electric field intensity of the BSW structure under the incident OAM beam with TC of 3. The edge length of the square region surrounded by the white dashed line is 2 µm. (b) Spatial phase patterns of OAM beams with TCs of 1 and 3 (from top to bottom) (left) and electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 1 and 3 (right). The edge length of the square region surrounded by the white dashed line is 2 µm. (c) Electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 1 and 3 (top), and the electric field intensity along the green dashed line (bottom). (d) Electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 3 and 3 (top), and the electric field intensity along the green dashed line (bottom). (e) Electric field intensity of the BSW structure under the incidence of two OAM beams with TCs of 5 and 3 (top), and the electric field intensity along the green dashed line (bottom); (c)–(e) share the same color scale bar. The green dashed line indicates the position where Y=0. The purple dashed lines indicate the peak positions of the two focal points.

    The reason for choosing five semi-ring shaped nanoslits is explained below. As we mentioned before, the period of the nanoslits was fixed at 1.05 µm to enhance the energy at the focal point. Figure 6(a) shows the electric field intensity of the BSW structure with five semi-ring nanoslits. This structure was illuminated by a circularly polarized OAM beam with TC of one from the bottom. A bright focal point can be clearly observed around the center with a transverse shift. We focus on the central square region with edge lengths of 2 µm (surrounded by the white dashed line). The BSW structures with one to six semi-ring shaped nanoslits are schematically plotted in Fig. 6(b) (top), and the electric field intensities of the corresponding structures under incident OAM beams is shown in Fig. 6(b) (bottom). It shows that the peak positions of the focal points of the BSW structures with different nanoslits remain the same, while the electric field intensity of the focal points increases gradually. The electric field intensity is quantitatively calculated by averaging the intensity at a small region around the focal point (0.1  μm×0.1  μm), as plotted in Fig. 6(c). It clearly shows that the electric field intensity of the BSW structures becomes larger when the nanoslit number increases, but it gets saturated gradually. The electric field intensity of the BSW structure with five nanoslits is 99.3% that of the structure with six nanoslits. The scattering loss caused by the inner nanoslits is the main reason. The transmission coefficient t is defined as the percentage of the energy that can reach the focal point after scattering by the inner nanoslit. The transmission coefficient t of each nanoslit is calculated to be 0.383. Therefore, the ratio of the electric field intensity of the five-nanoslit structure to that of the six-nanoslit structure can be thought of as the ratio between two geometric series with five and six terms, respectively. The ratio is (1t5)/(1t6)=0.995, which is consistent with the result in Fig. 6(c). To make a balance between high electric field intensity and a small footprint, the BSW structures with five semi-ring shaped nanoslits were employed throughout this study.

    Comparison of BSW structures with different semi-ring shaped nanoslits. (a) Electric field intensity of the BSW structure with five semi-ring nanoslits. The structure was illuminated by a circularly polarized OAM beam with TC of one from the bottom. The edge length of the square region surrounded by the white dashed line is 2 µm. (b) Schematic of the BSW structures with one to six semi-ring shaped nanoslits (top) and the electric field intensity of the corresponding structures (bottom) under incidence of an OAM beam with TC of one. All the subfigures share the same scale bar. (c) Normalized electric field intensity of the region around the peak positions of focal points as a function of nanoslit numbers.

    Figure 6.Comparison of BSW structures with different semi-ring shaped nanoslits. (a) Electric field intensity of the BSW structure with five semi-ring nanoslits. The structure was illuminated by a circularly polarized OAM beam with TC of one from the bottom. The edge length of the square region surrounded by the white dashed line is 2 µm. (b) Schematic of the BSW structures with one to six semi-ring shaped nanoslits (top) and the electric field intensity of the corresponding structures (bottom) under incidence of an OAM beam with TC of one. All the subfigures share the same scale bar. (c) Normalized electric field intensity of the region around the peak positions of focal points as a function of nanoslit numbers.

    The contributions from different electric field components (|Ex|, |Ey|, and |Ez|) are further analyzed. Electric field components of BSW and metallic structures under incident OAM beams with TC of five are presented in Appendix D. For the BSW structure, it can be clearly observed that the magnitudes of in-plane electric field components (|Ex| and |Ey|) are much larger than that of the out-of-plane electric field component (|Ez|). It proves that the BSW is excited by the TE polarization. On the contrary, for the metallic structure, the magnitudes of |Ex| and |Ey| are much smaller than that of |Ez|, which proves that the SPP is excited by the TM polarization. Although we only perform numerical simulations in this study, we propose the possible methods to measure the positions of the focal points in experiments. Near-field scanning optical microscopy (NSOM) with spatial resolution below the diffraction limit is a powerful technique, which has been used to characterize nanostructures with complex morphologies. The typical spatial resolution of an NSOM system is around 10–20 nm [63,64], which is capable of measuring the positions of the focal points in this study. Depending on the characteristics of the probes, NSOM can be generally divided into two approaches. Aperture-type NSOM employs optical fibers as probes, and it is sensitive to the in-plane electric field component [6567]. An aperture-type NSOM system can therefore be used to measure the focal points on the BSW structures. In contrast, apertureless NSOM employs machined tips (such as metallic nanoparticles) as probes, which is more sensitive to the out-of-plane electric field component [6567]. An apertureless NSOM system can be used to measure the focal points on the metallic structures.

    4. CONCLUSION

    In summary, on-chip sorting of OAM beams using BSW structures is numerically investigated in this study. The designed BSW structures consist of three dielectric layers, which are Si (top), SiO2 (middle), and Si substrate (bottom). Several semi-ring shaped nanoslits were introduced on the top surface of the BSW structures, and the number of nanoslits was optimized to be five. BSWs can be excited and modulated by the helical phase wavefront of the incident OAM beams. There is a transverse shift for the focal point owing to the constructive interference of the BSWs. An on-chip spatial interval of 135 nm between two neighboring OAM states can be obtained by using BSW structures. Furthermore, OAM sorting behaviors of metallic structures with the same morphology are also analyzed. OAM beams with TCs up to 35 can be successfully sorted by the BSW structures. As a comparison, it is difficult to detect OAM beams with TCs larger than seven by using metallic structures. BSW structures exhibit better OAM sorting performances than SPPs in terms of higher electric field intensity and smaller FWHM, which is beneficial for OAM sorting with large TCs. In addition, sorting of two OAM beams simultaneously can be achieved in this way as long as the difference between these two TCs is no less than four. By employing an ideal lossless material, we demonstrate that the greatly reduced loss of the BSW structure accounts for a big improvement of the detection limit from TC of 7 to TC of 26; the remaining improvement accounts for the smaller effective wavelength of the BSWs, which results in a tighter and more intense focus. To our knowledge, this is the first demonstration of on-chip sorting of OAM beams with TCs up to 35 as well as the first report of on-chip sorting of two OAM beams simultaneously. Our results reveal that BSW structures can be a remarkable solution for OAM sorting with large TCs. The investigation of BSWs will undoubtably broaden the OAM sorting strategies, and will be potentially applied in optical field manipulation, integrated photonic circuits, optical communications, and on-chip device design.

    Acknowledgment

    Acknowledgment. The authors acknowledge the Photonics Center of Shenzhen University for technical support.

    APPENDIX A: MECHANISM OF OAM SORTING

    The mechanism of the orbital angular momentum (OAM) sorting methods is schematically illustrated in Fig. 1(d). We first consider the case of only one nanoslit. Assuming the incident light possesses a phase profile of exp(ilφ) with radial polarization, the electric field at the surface can be expressed as [46] E(r,φ)=0π|r0r|12exp(ilφ)exp(ikBSW|r0r|)dφ,where r is the position vector of the observing point P. r0 is the position vector of an arbitrary point on the semi-ring nanoslit. l denotes the topological charge (TC) of the incident OAM beams. φ is the azimuthal angle. kBSW refers to the wavevector of the Bloch surface wave (BSW). Due to the helical phase of the OAM beams, the phase along the semi-ring nanoslit has a linear dependence on the azimuthal angle φ. In this case, the excited BSWs do not focus at the original point O, and there is a transverse shift for the focal region. Since the size of the focal region is much smaller than the diameter of the semi-ring nanoslit, the focal region can be treated as a single point. The above equation can be further expressed as [68,69] E(x)Jl(kBSWx),where Jl is the lth-order Bessel function. The transverse peak position of the focal point can therefore be derived to xl=sgn(l)μlkBSW.

    The above equation can be further simplified to xl=sgn(l)μl2πλBSW,where μl is the first non-zero coefficient of Jl(x) [the derivative of Jl(x)] and μl=0=0; λBSW is the effective wavelength of BSW. According to Eq. (A4), the transverse peak position of the focal point is only determined by l and λBSW, since μl is fixed for a specific material.

    When the incident light is a plane wave (l=0), the excited BSW focuses at the original position O. For incident OAM beams with positive TC (l>0), the excited BSWs focus at the left side of O. Correspondingly, for incident OAM beams with negative TC (l<0), the excited BSWs focus at the right side of O due to symmetry. In addition, incident OAM beams with a larger absolute value of TC |l| result in a larger transverse shift of peak position |x|.

    For multiple semi-ring nanoslits with different radii, the transverse position of the focal point remains the same since Eq. (A4) is independent of the radius of the nanoslit. The electric field intensity at the focal point can be enhanced or weakened depending on the phase differences of the BSWs excited by different nanoslits. Therefore, the period of the nanoslits (the radial distance between two neighboring slits) was set to be λBSW to enhance the electric field intensity at the focal point.

    In addition, the above analysis of the OAM sorting mechanism can be further extended to linearly polarized and circularly polarized incident light. For linearly polarized light with fixed polarization angle ζ, the electric field vector can be decomposed into two components. One component is along the radial direction, and the other component is perpendicular to the radial direction. Since only the component along the radial direction can excite BSWs, only this component will be considered. The electric field at the surface can be expressed as E(r,φ)=0π|r0r|12cos(ζφ)exp(ilφ)exp(ikBSW|r0r|)dφ.

    In the direction where φ=ζ or φ=ζ±π, the electric field reaches the maximum value. In the direction where φ=ζ±2/π, the electric field becomes zero. The overall integral of the electric field under linear polarization is obviously smaller than that under radial polarization.

    For circularly polarized light, the phase item can be expressed as exp[i(l+s)φ], where s denotes the spin angular momentum (SAM). Both OAM and SAM can contribute to a transverse shift of the position of the focal point. In this case, the electric field at the surface can be expressed as E(r,φ)=0π|r0r|12exp[i(l+s)φ]exp(ikBSW|r0r|)dφ.Similarly, the transverse peak position of the focal point can be expressed as xj=sgn(j)μj2πλBSW,where j=l+s, and j is the total angular momentum. s equals one for right-handed circular polarization (RCP), while s equals 1 for left-handed circular polarization (LCP). Since the circular polarization state of incident light can be fixed (s=±1), the OAM can be derived from the total angular momentum (l=js=j±1). Taken together, the mechanism of OAM sorting methods is valid for radially polarized, linearly polarized, and circularly polarized incident light.

    APPENDIX B: OAM SORTING PERFORMANCE COMPARISON

    OAM sorting performances of metallic structure and BSW structure are compared in Figs. 79. The electric field intensity of metallic structure under an incident OAM beam with TCs of 2, 4, 6, 9, 10, and 12 is shown in Fig. 7. The detection limit of OAM beams using this metallic structure is 7. As a comparison, the electric field intensity of BSW structure under an incident OAM beam with TCs of 10, 20, 30, 37, 38, and 40 is shown in Fig. 8. The detection limit of OAM beams using this BSW structure is 35. The electric field intensity of metallic structure and BSW structure along the central line under incident OAM beams with different TCs is further presented in Fig. 9. The electric field intensity of focal points on BSW structure is much higher than that of metallic structure (over eight times).

    Electric field intensity of the metallic structure under incident OAM beams with TCs of 2, 4, 6, 9, 10 and 12. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar.

    Figure 7.Electric field intensity of the metallic structure under incident OAM beams with TCs of 2, 4, 6, 9, 10 and 12. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar.

    Electric field intensity of the BSW structure under incident OAM beams with TCs of 10, 20, 30, 37, 38, and 40. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar.

    Figure 8.Electric field intensity of the BSW structure under incident OAM beams with TCs of 10, 20, 30, 37, 38, and 40. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar.

    (a), (b) Electric field intensity of the metallic structure along the Y=0 line under incidence of OAM beams with TCs of 1 to 7 with a broader X-axis scale (a) and a narrower X-axis scale (b). (c), (d) Electric field intensity of the BSW structure along the Y=0 line under incidence of OAM beams with TCs of 1 to 35 with a broader X-axis scale (c) and a narrower X-axis scale (d).

    Figure 9.(a), (b) Electric field intensity of the metallic structure along the Y=0 line under incidence of OAM beams with TCs of 1 to 7 with a broader X-axis scale (a) and a narrower X-axis scale (b). (c), (d) Electric field intensity of the BSW structure along the Y=0 line under incidence of OAM beams with TCs of 1 to 35 with a broader X-axis scale (c) and a narrower X-axis scale (d).

    APPENDIX C: OAM SORTING PERFORMANCE OF THE ARTIFICIAL MATERIAL

    The dielectric function of the artificial material is set as εa=116.62+0.01i. The effective wavelength of this artificial material is 1.54 µm, and its electric field intensity under incident OAM beams with TCs of 5, 10, 12, 13, 14, and 15 is shown in Fig. 10. The radius of the outermost semi-ring nanoslit was 12.2 µm. The detection limit of OAM beams using this artificial material is 13. As a comparison, the electric field intensity of this artificial material with larger radii under incident OAM beams with TCs of 5, 15, 25, 26, 27, and 28 is shown in Fig. 11. The radius of the outermost semi-ring nanoslit was 18.36 µm. The detection limit of OAM beams using this artificial material with larger radii is 26.

    Electric field intensity of the artificial structure under incident OAM beams with TCs of 5, 10, 12, 13, 14, and 15. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar. The radius of the outermost semi-ring nanoslit is 12.2 µm.

    Figure 10.Electric field intensity of the artificial structure under incident OAM beams with TCs of 5, 10, 12, 13, 14, and 15. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar. The radius of the outermost semi-ring nanoslit is 12.2 µm.

    Electric field intensity of the artificial structure with larger radius under incident OAM beams with TCs of 5, 15, 25, 26, 27, and 28. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar. The radius of the outermost semi-ring nanoslit is 18.36 µm.

    Figure 11.Electric field intensity of the artificial structure with larger radius under incident OAM beams with TCs of 5, 15, 25, 26, 27, and 28. The green dashed line indicates the position where Y=0. The white curve refers to the electric field intensity along the green dashed line. All figures share the same color scale bar. The radius of the outermost semi-ring nanoslit is 18.36 µm.

    APPENDIX D: ELECTRIC FIELD COMPONENT ANALYSIS

    Electric field components of BSW structure under incident OAM beams with a TC of 5 are shown in Fig. 12. For BSW structure, the magnitudes of |Ex| and |Ey| are much larger than that of |Ez|. It proves that the BSW is excited by the TE polarization. Electric field components of metallic structure under incident OAM beams with a TC of 5 are shown in Fig. 13. For metallic structure, the magnitudes of |Ex| and |Ey| are much smaller than that of |Ez|, which proves that the SPP is excited by the TM polarization.

    Electric field intensity (|E|2) and electric field components (|E|, |Ex|, |Ey|, and |Ez|) of the BSW structure under incident OAM beams with TC of five. The bottom three figures share the same color scale bar.

    Figure 12.Electric field intensity (|E|2) and electric field components (|E|, |Ex|, |Ey|, and |Ez|) of the BSW structure under incident OAM beams with TC of five. The bottom three figures share the same color scale bar.

    Electric field intensity (|E|2) and electric field components (|E|, |Ex|, |Ey|, and |Ez|) of the metallic structure under incident OAM beams with TC of five. The bottom three figures share the same color scale bar.

    Figure 13.Electric field intensity (|E|2) and electric field components (|E|, |Ex|, |Ey|, and |Ez|) of the metallic structure under incident OAM beams with TC of five. The bottom three figures share the same color scale bar.

    Tools

    Get Citation

    Copy Citation Text

    Nannan Li, Qi Zou, Yizhi Lan, Yaqi Wang, Jun Zhang, Michael Somekh, Changjun Min, Fu Feng, Xiaocong Yuan. On-chip sorting of orbital angular momentum beams using Bloch surface wave structures[J]. Photonics Research, 2023, 11(11): 1959

    Download Citation

    EndNote(RIS)BibTexPlain Text
    Save article for my favorites
    Paper Information

    Category: Surface Optics and Plasmonics

    Received: Aug. 10, 2023

    Accepted: Sep. 21, 2023

    Published Online: Nov. 3, 2023

    The Author Email: Michael Somekh (michael.somekh@nottingham.ac.uk), Changjun Min (cjmin@szu.edu.cn), Fu Feng (fufeng@zhejianglab.com), Xiaocong Yuan (xcyuan@szu.edu.cn)

    DOI:10.1364/PRJ.502760

    Topics