Next Article in Journal
The Effect of Silica Particle Size on the Mechanical Enhancement of Polymer Nanocomposites
Previous Article in Journal
Influence of Bulk Nanobubbles Generated by Acoustic Cavitation on Powder Microstructure and Rehydration Characteristics of Spray-Dried Milk Protein Concentrate Powders
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Integrated Generation and Shaping of Airy and Bessel Vortex Beams Based on All-Dielectric Metasurface

Guangdong Provincial Key Laboratory of Nanophotonic Functional Materials and Devices, School for Information and Optoelectronic Science and Engineering, South China Normal University, 378 Waihuan West Road Panyu District, Guangzhou 510006, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(6), 1094; https://doi.org/10.3390/nano13061094
Submission received: 22 February 2023 / Revised: 6 March 2023 / Accepted: 10 March 2023 / Published: 17 March 2023
(This article belongs to the Section Nanophotonics Materials and Devices)

Abstract

:
Integrating multiple independent functions into a single optical component is one of the most important topics in research on photoelectric systems. In this paper, we propose a multifunctional all-dielectric metasurface that can achieve a variety of non-diffractive beams depending on the polarization state of the incident light. Using the anisotropic TiO2 rectangular column as the unit structure, the three functions of generating polygonal Bessel vortex beams under left-handed circularly polarized incidence, Airy vortex beams under right-handed circularly polarized incidence and polygonal Airy vortex-like beams under linearly polarized incidence are realized. In addition, the number of polygonal beam sides and the position of focal plane can be adjusted. The device could facilitate further developments in scaling complex integrated optical systems and fabricating efficient multifunctional components.

1. Introduction

In 1992, Allen et al. [1] discovered for the first time that photons can carry orbital angular momentum (OAM) in addition to general spin angular momentum (SAM) [2], which laid the foundation for the study of vortex beams. Optical vortex beam (OVB) is a special beam carrying an OAM of per photon. Due to its singularity effect, OVB has a helical phase exp(i) in the wavefront structure, where l is the topological charge (TC) number which can take any integer value, φ denotes the azimuthal angle and ħ is the Planck constant. Since then, OVB has been widely used in quantum information processing [3,4], quantum entanglement [5], optical micromanipulation [6,7], optical communication [8,9,10] and nonlinear optics [11,12]. Hence, with the increasing demand for practical applications, OVB needs to be further studied: from generating [13,14,15] and measuring [16] to shaping OVB [17]. Especially in recent years, shaping OVB has a research hot spot. Current common shaping methods such as multi-OVB interaction [18,19,20] and optical oscillatory elements [21] require complex holography or high-precision diffractive elements. Until 2019, a new type of phase structure called higher-order cross-phase (HOCP) was introduced into the modulation of the beam. The function of HOCP is to introduce astigmatism into OVB in phase, so that the optical field mode changes during the propagation process, thereby changing the shape of OVB [22]. Compared with previous methods, it has higher integration and simpler implementation, which opens up a new horizon for shaping OVB.
However, in practical applications, a single OVB will be affected by many clutters in the propagation, which will deteriorate the performance. In order to further stabilize and improve the performance of OVB, a combination of OVB with other beams is a solution. In 1979, Berry and Balazs [23] first theoretically proved that the Schrödinger equation has a solution in the form of an Airy function and predicted a non-diffractive Airy wave packet. The optical version of the Airy wave packet is known as the Airy beam. It has the same properties of propagation intensity distribution invariance [24], non-diffraction [25,26] and self-healing [27] as the traditional non-diffracting beams (such as Bessel beam). However, unlike traditional non-diffracting beams, Airy beam can undergo a lateral self-accelerating process without any external force [28], which is one of its most remarkable features. Although Airy beam possesses novel propagation properties, theoretically, infinite energy is required to produce Airy beam. Therefore, it was not until 2007, after the first experimental demonstration that Airy beam can be realized in finite energy by Civiroglu et al. [29], that Airy beam started to be used in many studies, including the Airy laser [30,31,32], particle manipulation [33,34,35], space-time wave packets [36,37] and so on. Further, by superimposing the vortex phase on Airy beam, the Airy vortex beam can be generated. As a special vortex, Airy vortex beam combines the dual characteristics of Airy beam and vortex singularity, and it has good prospects in medical treatment [38]. Traditional methods of generating Airy beams and vortex beams rely on the use of spatial light modulators [39,40] or diffractive optical elements [41,42], among others. Nonetheless, the above methods are often greatly limited in practical applications due to their complex structures or low efficiency.
As a reduced-dimensional representation of metamaterials, the two-dimensional (2D) metasurface is composed of subwavelength structures that can flexibly manipulate the amplitude, phase and polarization of incident light waves, and it has developed rapidly in recent years. Compared with traditional bulk optical elements, metasurface has the advantages of small size and flexible design adjustment. Due to its unique physical properties, metasurface has been widely used in various applications such as anomalous refraction and reflection [43], spin Hall effect [44,45], stealth [46] and hologram [47,48]. Since optical metasurface is sensitive to the polarization state of light, a series of polarization-sensitive multifunctional metasurface has been extensively explored and studied by taking advantage of this feature. For example, Wen et al. realized a multifunctional metasurface which could produce a single Airy beam and arrays of Airy beams in 2020 [49]. Li et al. formed a vortex beam with multi-channel polarization conversion on a single metasurface in 2021 [50]. These two studies not only prove the enthusiasm of research on multi-channel metasurface, but also they demonstrate the needs of generating Airy beams and vortex beams in light field regulation.
Combining with HOCP, this paper proposes an all-dielectric multifunctional metasurface to realize multiple polygonal vortex beams, as shown in Figure 1. The design of the metasurface is based on the use of anisotropic TiO2 structure as a unit. To achieve an adjustable polarization sensitivity difference in the x and y directions, the phase of transmitted light is modulated by the combination of the PB phase and the propagation phase. The polarization principle of the PB phase can be understood by the Poincaré sphere. Using the conventional visible light wavelength = 532.8 nm) as incident wavelength, under the incidence of right-handed circularly polarized (RCP) light, the generation of polygonal Bessel beam can be realized based on HOCP, where the number of the sides of polygonal beam is equal to the order of HOCP. Additionally, a non-diffracting Airy beam with a controllable focus position is generated in the case of the left-handed circularly polarized (LCP) light incidence. Moreover, under the linearly polarized (LP) incidence, the amplitudes and phases of the two beams are superimposed on each other, so both the Airy vortex-like beam and polygonal shaping can be realized. We only need to change the order of the HOCP to obtain a desired polygonal vortex beam, which greatly simplifies the shaping process. Compared with other multifunctional metasurfaces, we only use the difference in polarization sensitivity to realize the generation of three kinds of non-diffracting beams, and add beam shaping function on the basis of traditional vortex, which has certain reference significance for the development of multifunctional metasurfaces. Further, since the modulated beam has a stable distribution at the far-field focus, it is of great significance to particle trapping, 3D optical tweezers, parallel laser printing and multiplex fluorescence imaging.

2. Theoretical Analysis and Design Method

2.1. Design of Multifunctional Metasurface

In this section, we will describe the physical mechanism that achieves the polarization sensitivity difference in the x and y direction in the transmission mode. We consider an arbitrary anisotropic nanopillar unit placed parallel to the z-axis, with the long and short axes along the x and y directions, respectively. When the unit is excited by the incident light with different polarization states, it produces responses with different phases and amplitudes along the x and y directions. Under LP light incidence, the Jones matrix of the transmitted field for an arbitrarily anisotropic unit rotated by an angle around the z-axis in the xoy plane is expressed as follows:
T ( x , y ) = R ( θ ) ( T x e i φ x 0 0 T y e i φ y ) R ( θ ) = ( cos θ sin θ sin θ cos θ ) ( T x e i φ x 0 0 T y e i φ y ) ( cos θ sin θ sin θ cos θ )
where Tx, Ty, φx and φy represent the transmission amplitudes and phases along the x and y axes. θ is the rotation angle of the unit relative to the reference coordinate system. R represents a 2 × 2 rotation matrix. Here, we assume that the transmission amplitude of the unit remains unity amplitude (Tx = Ty = 1) with a π phase difference (φx = φy + π) for the x and y polarization cases.
The independent phase modulation of LCP and RCP on the metasurface can realize the generation of φ1(x, y) and φ2(x, y) two different phase distribution beams. This means that when the RCP is incident, the metasurface will only be modulated by φ1(x, y). Similarly, the metasurface will only generate a beam modulated by φ2(x, y) when the LCP is incident. This modulation mechanism on metasurface can be represented by a matrix J(x, y):
J ( x , y ) | L C P = exp [ i φ 1 ( x , y ) ] | R C P
J ( x , y ) | R C P = exp [ i φ 2 ( x , y ) ] | L C P
The matrix form of | L C P corresponds to [ 1 i ] , and | R C P corresponds to [ 1 i ] . After the matrix transformation in Equations (2) and (3), the matrix J(x, y) can be expressed as:
J ( x , y ) = [ e i φ 1 ( x , y ) e i φ 2 ( x , y ) i e i φ 1 ( x , y ) i e i φ 2 ( x , y ) ] [ 1 1 i i ] 1
Then, by making the transmission matrix T(x, y) of Equation (1) equivalent to the modulation matrix J(x, y) of Equation (4), the phase distribution along the x and y directions and the unit rotation angle can be calculated, respectively.
φ x ( x , y ) = [ φ 1 ( x , y ) + φ 2 ( x , y ) ] / 2
φ y ( x , y ) = [ φ 1 ( x , y ) + φ 2 ( x , y ) ] / 2 π
θ ( x , y ) = [ φ 1 ( x , y ) φ 2 ( x , y ) ] / 4
The arrangement of units in the metasurface is determined by Equations (5)–(7). The phase design is mainly determined by the combination of the Pancharatnam–Berry (PB) phase adjustment based on Equation (5) and the propagation phase adjustment based on Equation (7). The PB phase is based on the adjustment of the phase by rotating the direction of the half-wave plate under circularly polarized light, which can be explained by the Poincaré sphere [51]. As shown in Figure 2, the north and south levels of the Poincaré sphere, respectively, represent LCP and RCP, and each point on the equator line corresponds to a different linear polarization state, while the entire spherical surface The other points represent elliptical polarization states for different polarization states. When a point on the sphere rotates an angle Ψ with the center of the sphere as the center, the phase change shown on the sphere is twice the rotation angle, which is the three-dimensional interpretation of the geometric phase. Therefore, when the rotation angle of nanopillars varies from 0 to π, the PB phase of 0 to 2π can be obtained. Since this phase is only related to the orientation and size of units, and since the different polarization states of the incident light will produce different phase responses, this feature can be used to perform differential encoding of left-handed and right-handed incidence, thereby realizing multifunctional metasurfaces.

2.2. Design of the Unit Structure

In order to better realize the combination of the transmission phase and the PB phase, the design of the units in the metasurface needs to be adjusted by the rotation angle and size. The design of a unit is shown in Figure 3. Due to the high refractive index and relatively low loss of TiO2 dielectric material at visible wavelength, TiO2 semiconductor is selected as the rectangular unit material which placed on the SiO2 substrate. The refractive index parameters (n, k) of both TiO2 and SiO2 materials are taken from [52]. To further discover the realization mechanism of phase, we performed numerical simulations using a three dimensional finite difference time domain (FDTD) solutions. We use scripts to import the above theoretical content into FDTD, and obtain corresponding results under 3D simulation to support the above theoretical derivation. The size of a mesh cell we set is 0.02 µm, and the simulation wavelength is set as λ = 532 nm. In the case where the simulation sets this value to the wavelength condition, the real part (n) of refractive indices of TiO2 and SiO2 are n T i O 2 = 2.17 and n S i O 2 = 1.48 , and the imaginary part (k) of the refractive indexes are both approximately zero. Considering both integration and transmission effects in this design, we set H = 900 nm and P = 400 nm, and L and W were varied from 50 nm to 400 nm, respectively. Figure 3b shows the normalized transmittance and phase magnitude of the 11 units selected after the simulation. It can be clearly seen that the selected units satisfy the full coverage of the 2π range well, and has a constant phase response difference of approximately π for x and y polarization. We apply periodic boundary conditions in the x and y directions and implement a perfectly matched layer (PML) in the z direction. The unit can be thought as a rectangular waveguide with birefringent properties due to the asymmetry of its geometry. Therefore, the phase delay along the x and y axes can be achieved by changing L and W.
The 11 black dots in Figure 4 correspond to the units in Figure 3b from left to right and from top to bottom. Figure 4a,c depict the relationship between the size (L, W) of units and the phase shift and transmission coefficient under x-LP light, respectively. It can also be seen that the phase retardation can well span the entire 2π range, and the transmittance is always at a high value of more than 80%. Therefore, the desired phase combination can be obtained by rationally choosing the unit size. In addition, the selected units are well satisfied that the phase difference are approximately π under x-LP and y-LP incidence, which can be clearly drawn from Figure 4b. The calculated phase shifts and transmission coefficients of individual TiO2 unit under x-LP and y-LP light incidence are presented in Figure 4d. In a word, the independent regulation of RCP and LCP can be well achieved by all units.

3. Results and Discussion

3.1. Polygonal Bessel Vortex Beams

The phase φ1(x, y) of the cross phase (CP) in Cartesian coordinates (x, y) has the form [22]:
φ 1 ( x , y ) = u ( x p cos θ y q sin θ ) ( x p sin θ + y q cos θ )
where the parameter u controls the conversion rate, the azimuth factor θ characterizes the rotation angle of the transformed beam on the Fourier plane, p and q are arbitrary positive integers, and the value of p plus q is equal to the order of CP. When CP is low-order cross phase (LOCP), that is, both p and q have a value of 1, the number of topological charges can be determined. Further, when the sum of p and q exceeds 2, CP as the high-order cross phase (HOCP) can realize adjusting the beam to an arbitrary polygon. This paper mainly exploits the ability of the CP to shape the beams, so the HOCP mode is adopted.
For computational convenience, we consider the case where the CP does not rotate (θ = 0). The phase φx(x, y) of adding HOCP to the Bessel vortex beam is:
φ x ( x , y ) = 2 π λ x 2 + y 2 sin ( β ) + l arctan ( y x ) + u x p y q + 2 η π
where λ represents the incident wavelength, β represents the angle between the transmitted wave and z-axis (set to 12.7°), l represents the topological charge, the parameter u controls the beam conversion rate and η is an arbitrary real number. It can be seen from Equation (9) that the magnitude of the value of u is affected by the order of CP. Since the designed structure is nanoscale, the value of u will be relatively large. As an example, in the following Figure 5, the generation process of the phase mask composed of the HOCP and the 2nd-order Bessel vortex beam is shown. Figure 5a1–a4 corresponds to the phase masks from the 3rd-order CP to the 6th-order CP, respectively, and Figure 5b1–b4 shows the phase of the 2nd-order Bessel vortex beam after superimposing HOCP in Figure 5a1–a4. It can be clearly seen that the Bessel vortex beam can be shaped into any polygon after adding HOCP, where the number of polygon sides is equal to the order of HOCP. Another way of saying this is that when the order of HOCP is three, four, five and six, the resulting OVB is shaped into triangular, quadrilateral, pentagonal and hexagonal polygons.
Based on the above theory, considering the influence of p and q values on beam shaping, the xoy plane electric field simulated at the focal position under RCP incidence is shown in Figure 6. Figure 6a1–d1 represents the electric field and phase distribution of the 2nd-order and 4th-order Bessel vortex beams without adding HOCP, respectively. It can be observed that Bessel vortex beams of any order can be well generated. When HOCP is added, the electric field presents a polygonal ring of corresponding order, and the phase presents a polygonal distribution. Among them, some basic properties are the same as the traditional Bessel vortex light, for example: the radius of the beam will increase with the increase in the order of the Bessel vortex beam, and there will be secondary halos of various levels outside the primary halo. From Figure 6b1–b5 and Figure 6d1–d5, it can be clearly known that the adjustment of HOCP is based on the phase. By polygonal shaping the vortex part of the phase mask, the shape of the spatial beam can be changed, which can confirm the numerical results derived in Figure 5. Since the manipulated particle will move in the direction of the energy flow, as long as the values of p and q are controlled according to the desired result, the motion trajectory of the manipulated instance can be accurately determined.

3.2. Airy Beams

Secondly, the phase of the Airy beam is:
φ y ( x , y ) = τ 3 3 ( m 3 x 3 + n 3 y 3 ) π
where τ = k/f, k = 2π/λ, λ represents the wavelength of the incident light, f is the focused wavelength, and m and n are arbitrary constants. To verify that the phase structure can realized the Airy beams, we assume the complex amplitude of the plane wave moving along the z-axis:
E ( x 0 , y 0 , z ) = A exp ( i k z )
After Fraunhofer diffraction of light field, the expression for the outgoing electric field is:
E ( x , y ) = exp ( i k z 1 ) i λ z 1 exp [ i k 2 z 1 ( x 2 + y 2 ) ] E ( x 0 , y 0 , z ) exp [ i k z 1 ( x x 1 + y y 1 ) ] exp [ i φ y ( x 1 , y 1 ) ] d x 1 d y 1
where z1 is the distance from the observation screen to the diffraction screen. Setting z1 = f and substituting Equations (10)–(12), the change in the electric field can be deduced:
E ( x , y ) = exp ( i k z 1 ) i λ z 1 exp [ i τ 2 ( x 2 + y 2 ) ] E ( x 0 , y 0 , z ) exp [ i τ ( x x 1 + y y 1 ) ] + i τ 3 3 ( m 3 x 1 3 + n 3 y 1 3 ) ] d x 1 d y 1
E ( x , y ) = 1 | m n | exp ( i k z 1 ) i λ z 1 exp [ i τ 2 ( x 2 + y 2 ) ] E ( x 0 , y 0 , z ) A i ( x m ) A i ( y n ) d x 1 d y 1
As can be seen from Equation (14), the light field of an Airy beam is affected by m and n. We set m = n = a, where a is an any constant. Corresponding to different a values under LCP incidence, the xoz and yoz electric field distributions and the intensity distributions at the focus position (the dotted line position in Figure 7a1–a5 are shown in Figure 7. Figure 7c shows the intensity distribution of the light field along the x-axis at the focus. With the increase in the value of a, the focal distance will gradually approach the metasurface and the central intensity of the focal position will progressively increase. In other words, the Airy beam becomes more concentrated as the value of a increases.

3.3. Airy Vortex-like Beams

Next, according to the focus position of the Bessel vortex beam, we correspondingly select the Airy beam with a = 16 × 10−7. Firstly, considering the case of the general Bessel vortex beam, when LP light is incident, the light fields generated by the incident LCP and RCP will act simultaneously. We slightly adjusted the output intensity of the two beams in order to achieve better fusion under the incidence of LP light. Additionally, we adjust the two light fields to focus on the same plane, and there will be an interaction between the two light fields, resulting in an Airy vortex-like beam, which can be verified by Figure 8. The light field and phase distribution of the xoy cross-section at the focal position are shown in Figure 8b,c. It can be found that after the superposition of the light fields, the central spot of the Airy beam is well integrated into the vortex part, forming a series of concentric rings similar to the Airy vortex beam, and its vortex state can be further confirmed in the phase field. The beam combines the singularity properties of vortices and the self-healing properties of Airy beams, which makes it promising for future laser applications.
Finally, under the control of adding HOCP to the Bessel vortex beam, we analyze the polygonal Airy vortex-like beam generated by introducing HOCP. Figure 9 shows the light field schematic diagrams of the xoz plane at the focus position, corresponding to the polygonal Airy vortex-like beams generated after the introduction of HOCP from the 3rd order to the 6th order, respectively. The generation of triangular, quadrilateral, pentagonal and hexagonal Airy vortex-like beams is well achieved at the focal position. It can be clearly analyzed in Figure 9 that the Airy beam has been well integrated into the polygonal Bessel beam, and the beam still retains the properties of the Airy beam. This design provides a reference value for the study of special vortex beams in the optical field.

4. Conclusions

In summary, we simulated an all-dielectric multifunctional metasurface with high transmittance by FDTD and realized three different functions of generating polygonal Bessel vortex beams, Airy beams and polygonal Airy vortex-like beams. The generation of polygons mainly depends on a new phase structure, namely, HOCP. We confirmed from theoretical calculations and numerical simulations that the metasurface can well shape general vortex beams into arbitrary polygons by controlling the order of the HOCP. Additionally, it is further discovered that the HOCP only changes the light field intensity distribution but not the OAM spectrum. All generated beams have relatively narrow beam widths and extended undiffracted propagation distances. This simple and general approach provides another feasible way to realize the design of multifunctional components. Notably, the functions generated by the metasurface do not interfere with each other, which means the system possesses potential applications in the fields of particle trapping, 3D optical tweezers, parallel laser printing and multiplex fluorescence imaging.

Author Contributions

Conceptualization, K.G.; methodology, K.G.; software, K.G.; validation, K.G. and H.L.; formal analysis, K.G., Y.L., Z.W. and H.L.; investigation, K.G. and H.L.; data curation, K.G.; writing—original draft preparation, K.G.; writing—review and editing, Y.L., Z.W. and H.L.; visualization, K.G.; supervision, H.L.; project administration, H.L.; funding acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China (62175070, 61875057); The Natural Science Foundation of Guangdong, China (2021A1515012652); and Science and Technology Program of Guangzhou (2019050001).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allen, L.; Beijersbergen, M.; Spreeuw, R.; Woerdman, J. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Phys. Rev. A 1992, 45, 8185–8189. [Google Scholar] [CrossRef] [PubMed]
  2. Beth, R. Mechanical Detection and Measurement of the Angular Momentum of Light. Phys. Rev. 1936, 50, 115–125. [Google Scholar] [CrossRef]
  3. Ge, F.; Chen, B.; Wang, Y.; Zhou, F.; Zheng, R.; Yang, X.; Qian, P.; Xu, N. A wideband balun-based microwave device for quantum information processing with nitrogenvacancy centers in diamond. J. Light. Technol. 2022, 40, 7572–7577. [Google Scholar] [CrossRef]
  4. Li, Z.; Ruan, Y.; Chen, P.; Tang, J.; Hu, W.; Xia, K.; Lu, Y. Liquid crystal devices for vector vortex beams manipulation and quantum information applications. Chin. Opt. Lett. 2021, 19, 112601. [Google Scholar] [CrossRef]
  5. Banerji, A.; Singh, R.P.; Banerjee, D.; Bandyopadhyay, A. Entanglement propagation of a quantum optical vortex state. Opt. Commun. 2016, 380, 492–498. [Google Scholar] [CrossRef] [Green Version]
  6. Qiu, S.; Wang, C.; Liu, T.; Ren, Y. Generation of spiral optical vortex with varying OAM for micro-manipulation. Opt. Commun. 2022, 524, 128767. [Google Scholar] [CrossRef]
  7. Maji, S.; Brundavanam, M.M. Topological transformation of fractional optical vortex beams using computer generated holograms. J. Opt. 2018, 20, 045607. [Google Scholar] [CrossRef]
  8. Shao, W.; Huang, S.; Liu, X.; Chen, M. Free-space optical communication with perfect optical vortex beams multiplexing. Opt. Commun. 2018, 427, 545–550. [Google Scholar] [CrossRef]
  9. Karahroudi, M.K.; Moosavi, S.A.; Mobashery, A.; Parmoon, B.; Saghafifar, H. Performance evaluation of perfect optical vortices transmission in an underwater optical communication system. Appl. Opt. 2018, 57, 9148–9154. [Google Scholar] [CrossRef]
  10. Yan, X.; Guo, L.; Cheng, M.; Li, J. Controlling abruptly autofocusing vortex beams to mitigate crosstalk and vortex splitting in free-space optical communication. Opt. Express 2018, 26, 12605–12619. [Google Scholar] [CrossRef]
  11. Engay, E.; Rodrigo, P.J. Nonlinear optical vortex coronagraph. Opt. Lett. 2020, 45, 1579–1582. [Google Scholar] [CrossRef] [PubMed]
  12. Tang, Y.; Li, K.; Zhang, X.; Deng, J.; Li, G.; Brasselet, E. Harmonic spin–orbit angular momentum cascade in nonlinear optical crystals. Nat. Photonics 2020, 14, 658–662. [Google Scholar] [CrossRef]
  13. Yang, Y.; Ye, X.; Niu, L.; Wang, K.; Yang, Z.; Liu, J. Generating terahertz perfect optical vortex beams by diffractive elements. Opt. Express 2020, 28, 1417–1425. [Google Scholar] [CrossRef] [PubMed]
  14. Wang, B.; Liu, W.; Zhao, M.; Wang, J.; Zhang, Y.; Chen, A.; Guan, F.; Liu, X.; Shi, L.; Zi, J. Generating optical vortex beams by momentum-space polarization vortices centred at bound states in the continuum. Nat. Photonics 2020, 14, 623–628. [Google Scholar] [CrossRef]
  15. Kitamura, K.; Kitazawa, M.; Noda, S. Generation of optical vortex beam by surface-processed photonic-crystal surface-emitting lasers. Opt. Express 2019, 27, 1045–1050. [Google Scholar] [CrossRef]
  16. Zhao, Q.; Dong, M.; Bai, Y.; Yang, Y. Measuring high orbital angular momentum of vortex beams with an improved multipoint interferometer. Photonics Res. 2020, 8, 745–749. [Google Scholar] [CrossRef]
  17. Ma, H.; Li, X.; Zhang, H.; Tang, J.; Li, H.; Tang, M.; Wang, J.; Cai, Y. Optical vortex shaping via a phase jump factor. Opt. Lett. 2019, 44, 1379–1382. [Google Scholar] [CrossRef]
  18. Cookson, T.; Kalinin, K.; Sigurdsson, H.; Töpfer, J.D.; Alyatkin, S.; Silva, M.; Langbein, W.; Berloff, N.G.; Lagoudakis, P.G. Geometric frustration in polygons of polariton condensates creating vortices of varying topological charge. Nat. Commun. 2021, 12, 2120. [Google Scholar] [CrossRef]
  19. Gong, L.; Wang, X.; Zhu, Z.; Lai, S.; Feng, H.; Wang, J.; Gu, B. Transversal energy flow of tightly focused off-axis circular polarized vortex beams. Appl. Opt. 2022, 61, 5076–5082. [Google Scholar] [CrossRef]
  20. Xia, T.; Cheng, S.; Tao, S. Two polygon-like beams generated by a modified interfering vortex spiral zone plate. Results Phys. 2021, 29, 104762. [Google Scholar] [CrossRef]
  21. Shen, Y.; Wan, Z.; Meng, Y.; Fu, X.; Gong, M. Polygonal vortex beams. IEEE Photonics J. 2018, 10, 1–16. [Google Scholar] [CrossRef]
  22. Ren, Y.; Wang, C.; Liu, T.; Wang, Z.; Yin, C.; Qiu, S.; Li, Z.; Wu, H. Polygonal shaping and multi-singularity manipulation of optical vortices via high-order cross-phase. Opt. Express 2020, 28, 26257–26266. [Google Scholar] [CrossRef] [PubMed]
  23. Berry, M.V.; Balazs, N.L. Nonspreading wave packets. Am. J. Phys. 1979, 47, 264–267. [Google Scholar] [CrossRef]
  24. Chen, R.P.; Ying, C.F. Beam propagation factor of an Airy beam. J. Optics 2011, 13, 085704. [Google Scholar] [CrossRef]
  25. Lin, Z.; Guo, X.; Tu, J.; Ma, Q.; Wu, J.; Zhang, D. Acoustic non-diffracting Airy beam. J. Appl. Phys. 2015, 117, 104503. [Google Scholar] [CrossRef]
  26. Minovich, A.E.; Klein, A.E.; Neshev, D.N.; Pertsch, T.; Kivshar, Y.S.; Christodoulides, D.N. Airy plasmons: Non-diffracting optical surface waves. Laser Photonics Rev. 2014, 8, 221–232. [Google Scholar] [CrossRef]
  27. Broky, J.; Siviloglou, G.A.; Dogariu, A.; Christodoulides, D.N. Self-healing properties of optical Airy beams. Opt. Express 2008, 16, 12880–12891. [Google Scholar] [CrossRef] [Green Version]
  28. Efremidis, N.K.; Chen, Z.; Segev, M.; Christodoulides, D.N. Airy beams and accelerating waves: An overview of recent advances. Optica 2019, 6, 686–701. [Google Scholar] [CrossRef] [Green Version]
  29. Siviloglou, G.A.; Christodoulides, D.N. Accelerating finite energy Airy beams. Opt. Lett. 2007, 32, 979–981. [Google Scholar]
  30. Porat, G.; Dolev, I.; Barlev, O.; Arie, A. Airy beam laser. Opt. Lett. 2011, 36, 4119–4121. [Google Scholar] [CrossRef] [Green Version]
  31. Campargue, G.; Zielinski, B.; Courvoisier, S.; Sarpe, C.; Winkler, T.; Sentfleben, A.; Bonacina, L.; Baumert, T.; Wolf, J.P. Live cells assessment of opto-poration by a single femtosecond temporal Airy laser pulse. AIP Adv. 2018, 8, 125105. [Google Scholar] [CrossRef] [Green Version]
  32. Hou, T.; Chang, Q.; Chang, H.; Liu, W.; Ma, P.; Zhou, P. Higher-Order Airy patterns and their application in tailoring orbital angular momentum beams with fiber laser arrays. J. Light. Technol. 2021, 39, 4758–4768. [Google Scholar] [CrossRef]
  33. Chen, M.; Huang, S.; Liu, X.; Chen, Y.; Shao, W. Optical trapping and rotating of micro-particles using the circular Airy vortex beams. Appl. Phys. B 2019, 125, 1–6. [Google Scholar] [CrossRef]
  34. Lu, W.; Sun, X.; Chen, H.; Liu, S.; Lin, Z. Abruptly autofocusing property and optical manipulation of circular Airy beams. Phys. Rev. A 2019, 99, 013817. [Google Scholar] [CrossRef]
  35. Shou, Q.; Kuang, W.; Liu, M.; Zhou, Z.; Chen, Z.; Hu, W.; Guo, Q. Two dimensional large-scale optical manipulation of microparticles by circular Airy beams with spherical and oblique wavefronts. Opt. Commun. 2022, 525, 128561. [Google Scholar] [CrossRef]
  36. Yessenov, M.; Free, J.; Chen, Z.; Johnson, E.G.; Lavery, M.P.; Alonso, M.A.; Abouraddy, A.F. Space-time wave packets localized in all dimensions. Nat. Commun. 2022, 13, 1–13. [Google Scholar] [CrossRef]
  37. Hall, L.A.; Abouraddy, A.F. Consequences of non-differentiable angular dispersion in optics: Tilted pulse fronts versus space-time wave packets. Opt. Express 2022, 30, 4817–4832. [Google Scholar] [CrossRef]
  38. Juhasz, T.; Loesel, F.H.; Kurtz, R.M.; Horvath, C.; Bille, J.F.; Mourou, G. Corneal refractive surgery with femtosecond lasers. IEEE J. Sel. Top. Quantum Electron. 1999, 5, 902–910. [Google Scholar] [CrossRef]
  39. Latychevskaia, T.; Schachtler, D.; Fink, H.W. Creating Airy beams employing a transmissive spatial light modulator. Appl. Opt. 2016, 55, 6095–6101. [Google Scholar] [CrossRef]
  40. Szatkowski, M.; Masajada, J.; Augustyniak, I.; Nowacka, K. Generation of composite vortex beams by independent Spatial Light Modulator pixel addressing. Opt. Commun. 2020, 463, 125341. [Google Scholar] [CrossRef]
  41. Siemion, A. The magic of optics—An overview of recent advanced terahertz diffractive optical elements. Sensors 2020, 21, 100. [Google Scholar] [CrossRef]
  42. Cao, F.; Zhao, Y.; Yao, C.; Xie, C. All diffractive optical element setup for creating and characterizing optical vortices with high topological charges: Analytical models and numerical results. Opt. Commun. 2021, 495, 127119. [Google Scholar] [CrossRef]
  43. Zhang, J.; Zhao, X.; Zheng, Y.; Chen, X. Generalized nonlinear Snell’s law at χ2 modulated nonlinear metasurfaces: Anomalous nonlinear refraction and reflection. Opt. Lett. 2019, 44, 431–434. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Li, S.; Li, X.; Wang, G.; Liu, S.; Zhang, L.; Zeng, C.; Wang, L.; Sun, Q.; Zhao, W.; Zhang, W. Multidimensional manipulation of photonic spin Hall effect with a single-layer dielectric metasurface. Adv. Opt. Mater. 2019, 7, 1801365. [Google Scholar] [CrossRef]
  45. Kim, M.; Lee, D.; Yang, Y.; Kim, Y.; Rho, J. Reaching the highest efficiency of spin Hall effect of light in the near-infrared using all-dielectric metasurfaces. Nat. Commun. 2022, 13, 1–7. [Google Scholar] [CrossRef]
  46. Wang, Z.; Wang, G.; Cui, X.; Bai, H. Hybrid metasurface-based broadband high gain stealth antenna. Opt. Express 2022, 30, 32833–32846. [Google Scholar] [CrossRef]
  47. Hermon, S.; Ma, A.; Yue, F.; Kubrom, F.; Intaravanne, Y.; Han, J.; Ma, Y.; Chen, X. Metasurface hologram for polarization measurement. Opt. Lett. 2019, 44, 4436–4438. [Google Scholar] [CrossRef]
  48. Li, Y.; Yang, Q.; Xiong, J.; Yin, K.; Wu, S.T. 3D displays in augmented and virtual realities with holographic optical elements. Opt. Express 2021, 29, 42696–42712. [Google Scholar] [CrossRef]
  49. Wen, J.; Chen, L.; Yu, B.; Nieder, J.B.; Zhuang, S.; Zhang, D.; Lei, D. All-dielectric synthetic-phase metasurfaces generating practical airy beams. ACS Nano 2021, 15, 1030–1038. [Google Scholar] [CrossRef]
  50. Li, J.; Zheng, C.; Li, J.; Wang, G.; Liu, J.; Yue, Z.; Hao, X.; Yang, Y.; Li, F.; Tang, T.; et al. Terahertz wavefront shaping with multi-channel polarization conversion based on all-dielectric metasurface. Photonics Res. 2021, 9, 1939–1947. [Google Scholar] [CrossRef]
  51. Berry, M. The adiabatic phase and Pancharatnam’s phase for polarized light. J. Mod. Opt. 1987, 34, 1401–1407. [Google Scholar] [CrossRef]
  52. Sarkar, S.; Gupta, V.; Kumar, M.; Schubert, J.; Probst, P.T.; Joseph, J.; König, T.A.F. Hybridized Guided-Mode Resonances via Colloidal Plasmonic Self-Assembled Grating. ACS Appl. Mater. Interfaces 2019, 11, 13752–13760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Schematic diagram of multifunctional metasurfaces for realizing non-diffracting beams. (a) Generation of polygonal Bessel vortex beam under RCP incidence. (b) Generation of Airy vortex-like beam under LP incidence. (c) Generation of Airy beam under LCP incidence.
Figure 1. Schematic diagram of multifunctional metasurfaces for realizing non-diffracting beams. (a) Generation of polygonal Bessel vortex beam under RCP incidence. (b) Generation of Airy vortex-like beam under LP incidence. (c) Generation of Airy beam under LCP incidence.
Nanomaterials 13 01094 g001
Figure 2. Poincaré sphere model.
Figure 2. Poincaré sphere model.
Nanomaterials 13 01094 g002
Figure 3. (a) Schematic perspective and top view of the TiO2 unit structure. (b) Corresponding normalized transmission coefficients (black) and phase shifts (blue) for selected 11 units under x-LP and y-LP incidences.
Figure 3. (a) Schematic perspective and top view of the TiO2 unit structure. (b) Corresponding normalized transmission coefficients (black) and phase shifts (blue) for selected 11 units under x-LP and y-LP incidences.
Nanomaterials 13 01094 g003
Figure 4. (a) Phase shift under x-LP incidence. (b) Absolute value of the phase offset difference between x-LP and y-LP. (c) Normalized transmission coefficient under x-LP incidence. (d) Normalized transmission coefficient under y-LP incidence. The black pentagrams represent the selected 11 units.
Figure 4. (a) Phase shift under x-LP incidence. (b) Absolute value of the phase offset difference between x-LP and y-LP. (c) Normalized transmission coefficient under x-LP incidence. (d) Normalized transmission coefficient under y-LP incidence. The black pentagrams represent the selected 11 units.
Nanomaterials 13 01094 g004
Figure 5. Calculation of phase distribution of Bessel vortex beam modulated by the HOCP. (a1a4) Phase distribution of HOCP. (b1b4) Phase distribution of the Bessel vortex beam after adding HOCP modulation. The parameters u in the 3rd-order CP, 4th-order CP, 5th-order CP and 6th-order CP are 3 × 1015, 4 × 1020, 5 × 1025 and 6 × 1030, respectively.
Figure 5. Calculation of phase distribution of Bessel vortex beam modulated by the HOCP. (a1a4) Phase distribution of HOCP. (b1b4) Phase distribution of the Bessel vortex beam after adding HOCP modulation. The parameters u in the 3rd-order CP, 4th-order CP, 5th-order CP and 6th-order CP are 3 × 1015, 4 × 1020, 5 × 1025 and 6 × 1030, respectively.
Nanomaterials 13 01094 g005
Figure 6. Simulation of the light field and phase distribution of Bessel vortex beams modulated by HOCP under RCP incidence. (a1a5) Optical field distribution of the 2nd-order Bessel vortex beam modulated by HOCP. (b1b5) Phase distribution of the 2nd-order Bessel vortex beam modulated by HOCP. (c1c5) Optical field distribution of the 4th-order Bessel vortex beam modulated by HOCP. (d1d5) Phase distribution of the 4th-order Bessel vortex beam modulated by HOCP. Corresponding to no HOCP, 3rd-HOCP, 4th-HOCP, 5th-HOCP and 6th-HOCP, respectively.
Figure 6. Simulation of the light field and phase distribution of Bessel vortex beams modulated by HOCP under RCP incidence. (a1a5) Optical field distribution of the 2nd-order Bessel vortex beam modulated by HOCP. (b1b5) Phase distribution of the 2nd-order Bessel vortex beam modulated by HOCP. (c1c5) Optical field distribution of the 4th-order Bessel vortex beam modulated by HOCP. (d1d5) Phase distribution of the 4th-order Bessel vortex beam modulated by HOCP. Corresponding to no HOCP, 3rd-HOCP, 4th-HOCP, 5th-HOCP and 6th-HOCP, respectively.
Nanomaterials 13 01094 g006
Figure 7. Simulation of Airy beam generation under the LCP incidence for different a values. (a1a5) Electric field in the xoz plane. (b1b5) Electric field in the xoy plane at the dashed position. (c) Normalized intensity contrast plot at the dotted line.
Figure 7. Simulation of Airy beam generation under the LCP incidence for different a values. (a1a5) Electric field in the xoz plane. (b1b5) Electric field in the xoy plane at the dashed position. (c) Normalized intensity contrast plot at the dotted line.
Nanomaterials 13 01094 g007
Figure 8. The (a1a3) xoz plane light field, (b1b3) xoy plane light field and (c1c3) xoy plane phase field of the beam generated by RCP incidence, LP incidence and LCP incidence were simulated without HOCP modulation.
Figure 8. The (a1a3) xoz plane light field, (b1b3) xoy plane light field and (c1c3) xoy plane phase field of the beam generated by RCP incidence, LP incidence and LCP incidence were simulated without HOCP modulation.
Nanomaterials 13 01094 g008
Figure 9. Simulation of polygonal similar Airy vortex beams modulated by adding HOCP, with (a) 3rd-HOCP, (b) 4th-HOCP, (c) 5th-HOCP, and (d) 6th-HOCP.
Figure 9. Simulation of polygonal similar Airy vortex beams modulated by adding HOCP, with (a) 3rd-HOCP, (b) 4th-HOCP, (c) 5th-HOCP, and (d) 6th-HOCP.
Nanomaterials 13 01094 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, K.; Liu, Y.; Wei, Z.; Liu, H. Numerical Simulation of Integrated Generation and Shaping of Airy and Bessel Vortex Beams Based on All-Dielectric Metasurface. Nanomaterials 2023, 13, 1094. https://doi.org/10.3390/nano13061094

AMA Style

Guo K, Liu Y, Wei Z, Liu H. Numerical Simulation of Integrated Generation and Shaping of Airy and Bessel Vortex Beams Based on All-Dielectric Metasurface. Nanomaterials. 2023; 13(6):1094. https://doi.org/10.3390/nano13061094

Chicago/Turabian Style

Guo, Kuangling, Yue Liu, Zhongchao Wei, and Hongzhan Liu. 2023. "Numerical Simulation of Integrated Generation and Shaping of Airy and Bessel Vortex Beams Based on All-Dielectric Metasurface" Nanomaterials 13, no. 6: 1094. https://doi.org/10.3390/nano13061094

APA Style

Guo, K., Liu, Y., Wei, Z., & Liu, H. (2023). Numerical Simulation of Integrated Generation and Shaping of Airy and Bessel Vortex Beams Based on All-Dielectric Metasurface. Nanomaterials, 13(6), 1094. https://doi.org/10.3390/nano13061094

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop