Next Article in Journal
Scanning Tunneling Microscopy of Biological Structures: An Elusive Goal for Many Years
Previous Article in Journal
Multiferroic and Magnetodielectric Effects in Multiferroic Pr2FeAlO6 Double Perovskite
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generation and Superposition of Perfect Vortex Beams in Terahertz Region via Single-Layer All-Dielectric Metasurface

1
State Key Laboratory of Transient Optics and Photonics, Xi’an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences, Xi’an 710119, China
2
School of Opto-Electronic Technology, University of Chinese Academy of Sciences, Beijing 100049, China
3
Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan 030006, China
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(17), 3010; https://doi.org/10.3390/nano12173010
Submission received: 7 August 2022 / Revised: 25 August 2022 / Accepted: 26 August 2022 / Published: 30 August 2022

Abstract

:
Terahertz (THz) orbital angular momentum (OAM) technology provides promising applications in future wireless communication with large bandwidth and high capacity. However, the ring radius of the conventional THz vortex beam is related to the topological charge, limiting the co-propagation of multiple OAM modes in the THz communication systems. Although the perfect vortex beam (PVB) based on traditional methods can solve this problem, they are usually bulky and unstable. Here, we demonstrate two PVB generators based on a single all-dielectric metasurface to obtain polarization-independent PVB and spin multiplexed PVB, respectively. The former regulates the propagation phase by using isotropic unit cells; the latter simultaneously manipulates the propagation and geometric phase to achieve the spin-decoupled phase control by arranging anisotropic unit cells. In addition, we also demonstrate the stable generation of a perfect Poincaré beam with arbitrary polarization and phase distribution on a hybrid-order Poincaré Sphere via a spin-decoupled metasurface, which is achieved by the linear superposition of two PVBs with orthogonal circular polarizations. The proposed scheme provides a compact and efficient platform for the generation and superposition of PVBs in THz region, and will speed up the progress of THz communication systems, complex light field generation, and quantum information sciences.

1. Introduction

It is well known that light has both spin angular momentum (SAM) and orbital angular momentum (OAM), which are manifested as circular polarization and azimuthal phase of a light beam, respectively [1,2]. Owing to the annular intensity distribution and the particular helical phase structure, light with OAM (namely vortex beams) has been widely used in optical communication [3], optical trapping [4], metrology [5], and quantum memories [6]. Terahertz (THz) waves located between far-infrared and microwave bands with rich frequency resources and have not been utilized extensively [7]. The related wireless THz communication has been recognized as one of the potential key technologies for high-speed communication in the future due to its great performance in high rate, broad bandwidth, and strong anti-interference [8,9]. The combination of THz waves and vortex beams will further extend communication capacity because different OAM modes are strictly orthogonal, and the multiple OAM states of vortex beams can be used to transmit different data [10,11]. However, the ring radiuses of these conventional vortex beams strongly depend on their topological charges, so it is not easy to combine multiple OAM beams into a single fiber used for multiplexed communications simultaneously [12,13].
A similar situation occurred in the vector vortex beams (VVBs) domain. Such spiral phase distribution and spatial inhomogeneous polarization state provide more interaction between light and matter, and have been utilized in high-resolution lithography [14], plasmonic couplers [15], and optical encryption [16]. Due to the inhomogeneous spatial polarization and phase distribution of VVBs, the Poincaré Sphere (PS) proposed to describe fundamental polarization states is not applicable [17]. Therefore, the higher-order Poincaré Sphere (HOPS) has been developed to visually describe the states of polarization and phase [18]. However, the HOPS model is still limited in specific cases because the OAM on the south and north poles is required to be opposite. In order to overcome this limitation, the hybrid-order Poincaré Sphere (HyOPS) has been proposed [19]. Similar to the PS, each point on a HyOPS describes a specific VVB represented by a linear superposition of the two poles with different phase and amplitude offsets. The beam on PS, HOPS, and HyOPS has been uniformly defined as the Poincaré beams (PBs), which describe most of polarized fields [20]. It has been proved that the superposition of vortex beams is an effective way to generate PBs [21]. However, this is usually complicated because the ring radius of conventional vortex beams is related to the topological charges, resulting in the polarization state being unstable and the intensity profile collapsing [22,23].
In order to overcome the above issues, Ostrovsky et al. introduced the concept of the perfect vortex beam (PVB), whose ring radius is independent of the topological charge [24,25]. Plenty of methods have been reported to obtain the PVB, including spatial light modulator [26], axicon [27], interferometer [28], and polymer-based phase plate [29]. Furthermore, the perfect Poincaré beam (PPB) theory that the ring radius is independent of the topological charge has been demonstrated using combinations of conventional optical components [30,31], such as axicons, spatial light modulators, q-plates, and Fourier transform lenses. However, these methods obtain the PVB and PPBs accompanied by complex optical path design, difficult fabrication, and bulky setup, running counter to the increasing miniaturization and integration requirements of photonic devices. Meanwhile, practical devices are absent in THz region, causing a challenge in generation and superposition of THz PVBs.
Metasurfaces, composed of subwavelength optical elements, can arbitrarily modulate light’s phase, polarization, and amplitude, providing an efficient, versatile platform for compact planar optics [32,33,34]. Despite many advances, the fixed wave-manipulation functionalities are still the main drawbacks in metasurfaces. In this case, many efforts have been devoted to achieving multifunctional and tunable metasurfaces that can dynamically control light upon external tuning. For example, active materials gain dynamic control of effective permeability, transmittance, and Fano resonances via the analytical designs of interfacial phase discontinuities [35,36,37]. Researchers have recently proposed a plasmonic metasurface in near-infrared (NIR) to generate PVBs based on geometric phase theory [38]. However, the transmission of the incident beam through the metasurface is deficient due to the large ohmic absorption caused by Au nanostructures. In 2021, Xie et al. designed a polarization-controllable dielectric metasurface in NIR to generate two different PVBs by switching the incident polarization utilizing propagation phase theory [39]. Compared to metal metasurface, dielectric metasurface displays higher work efficiency. In 2020, Bao et al. proposed a silicon metasurface to generate PPBs [40]. However, this design requires simultaneous phase and amplitude modulation as well as a complicated oblique incidence scheme, which limits its efficiency and operating bandwidth. In 2021, Zhou et al. demonstrated the efficient and broadband generation of PVBs at visible wavelengths using a single-layer dielectric metasurface [41]. However, the meta-device is single function, and the performance is sensitive to incident polarization due to the limitation of geometric phase theory, which will reduce its practicability. Meanwhile, a silicon metasurface has been designed to generate PBs in the THz range [42]. However, this work is only limited to the exceptional cases, such as radially or azimuthally polarized beams, and their annular intensity profiles are still dependent on the topological charge. Recently, Liu et al. proposed so-called quasi-perfect vortex beams (Q-PVBs), a new concept to realize vortex beams with a controllable divergence angle independent of topological charges. Moreover, the multichannel terahertz Q-PVBs generation is also theoretically and experimentally demonstrated [43]. However, this work is also limited by the geometric phase method. For example, it cannot be utilized to generate PVBs composed of arbitrarily circularly polarized vortex beams because there will always be one diverging circularly polarized output due to reversed phase sign. Therefore, to our best knowledge, a compact setup capable of generating an arbitrary PVB and further superposition to obtain PPBs in THz region is highly desirable and remains unexplored.
Here, two kinds of single-layer all-silicon metasurfaces have been designed to generate PVBs and PPBs in THz region. One generates polarization-independent PVBs by the propagation phase of cross-shaped unit cells. The other, made up of rectangular-shaped silicon pillars, can simultaneously employ the geometric and dynamic phase to generate spin-decoupled PVBs. Therefore, the radius and topological charges of the PVB can be conveniently varied by switching the incident spin states. Most importantly, two completely different PVBs can be mixed by a spin-decoupled metasurface. All orders of PPBs with different phase and polarization distributions can be generated by choosing the polarization states of incident light and designing topological charges of output PVBs. Thus, our work provides a potential compact platform for accelerating THz communication systems, complex light field generation, and quantum information sciences.

2. Materials and Methods

2.1. Polarization-Independent Metasurface

We design a metasurface that can generate polarization-independent PVBs. As shown in Figure 1a,b, the incident light with arbitrary polarization through the metasurface can produce a vortex beam such that radius is irrelevant to the topological charge. Theoretically, the complex amplitude of the PVB can be described as [24]:
E 0 ( ρ , θ ) = δ ( ρ - ρ 0 ) exp ( i l θ )
where ( ρ , θ )   represent the polar coordinates, δ ( ) is the Dirac function in the polar coordinate system, l is the topological charge, and ρ 0   is the radius of bright annular intensity. Due to the Dirac function only existing in the ideal state, this perfect vortex model is not accessible in practice. According to the reported work, PVBs can be generated by taking the Fourier transform of an ideal higher-order Bessel beam [26], which can be presented in the following form,
E B ( r , z ) = J l ( k r ρ ) exp ( i l θ + i k z z )
Here, J l   is the first kind of l-th order Bessel function; k = K r 2 + K z 2 = 2 π λ ,   ( k r , k z ) are the radial and longitudinal wave vectors, respectively, and λ is the incident wavelength. However, it is not easy to produce ideal Bessel beams experimentally. In most practical cases, the accessible Bessel model is the Bessel–Gaussian (BG) beam. Hence, the PVB beam is typically generated from the Fourier transform of the BG beam, which is usually achieved by a cascading optical setup, including a spiral phase plate, an axicon lens, and a Fourier transform lens [44]. Specifically, the generation of PVBs has two steps. First, a spiral phase plate and an axicon are used to convert the Gaussian beam to the corresponding BG beam. Second, a Fourier lens is utilized to transform the BG beam into PVBs. Due to the unprecedented light manipulating ability, a single-layer metasurface can work as a spiral phase plate, an axicon, and a Fourier transformation lens simultaneously, as shown in Figure 1c. The total phase profile can be expressed as:
φ t o t a l ( x , y ) = φ spiral ( x , y ) + φ a x i c o n ( x , y ) + φ l e n s ( x , y )
φ spiral ( x , y ) = l arctan ( x y )
φ axicon ( x , y ) = 2 π x 2 + y 2 d
φ l e n s ( x , y ) = π ( x 2 + y 2 ) λ f
where x and y represent the position coordinates of the unit cell in the metasurface concerning the center of the origin, Equation (4) represents the spiral phase profile of an optical vortex with the topological charge l. Equation (5) shows an axicon phase profile, and d is the axicon period relating to the ring radius of the PVB. Equation (6) introduces the phase profile of the Fourier transform lens with the focal length of f, and λ is the working wavelength.
In order to realize the design requirement, we used the finite-difference time-domain (FDTD) method to optimize the geometric parameters of the unit cell. Single-layer plasmonic structures, depending on metallic materials with varying geometries, have been utilized to construct metasurfaces. However, they suffer from low efficiency due to the inherent ohmic loss and weak coupling between the incident and cross-polarized fields [45]. Although multilayer plasmonic structures can increase the efficiency, they also increase fabrication difficulties [46]. In addition, plasmonic structures are sensitive to the polarization of incident light, making them unsuitable for practical applications. It has been proved that the most straightforward approach to achieve polarization-independent performance is to find the unit cell with 4-fold symmetry [47]. Therefore, as shown in Figure 2a, the cross-shaped unit cell with the lattice constant of P = 130 μm, the height of H = 300 μm, and the exact width W = 50 µm are selected at the working frequency 0.75 THz. High-resistance silicon (ρ > 5000 Ω m) has been chosen as the designed material due to the negligible absorption and very low dissipation for terahertz wave as well as mature etching techniques, and its dielectric constant is set as ε = 11.9 [48,49].
By scanning the length (L) of the unit cell from 40 μm to 125 μm, both the transmittance and phase shift have been calculated, as shown in Figure 2b. Moreover, the inset illustrates the side views of the simulated magnetic amplitude distributions when L is fixed as 120 μm. The magnetic fields are mainly localized inside the silicon pillars, indicating that the unit cell can be seen as a truncated low-loss waveguide that works as a weakly coupled low-quality factor Fabry–Perot resonator [50]. Next, as shown in Figure 2c, eight units are chosen to realize the complete phase control with high transmittance. Their phase shifts cover the 2π range, and the transmittances are maintained over 70%. Figure 2d shows phase wavefronts under the x-polarized wave incident at the pillar height H of 330 μm and the length L of 40, 82, and 120 μm, respectively. Clearly, the relative phase difference between these pillars at the output port (at 330 μm) corresponds to π, and 2π, respectively.

2.2. Spin-Decoupled Metasurface

We design a spin-decoupled metasurface to generate two completely different PVBs only by varying the incident spin state. As shown in Figure 3a, left-circularly polarized (LCP) incident light through the metasurface can be transformed into a PVB with opposite-handedness, and the same metasurface can also transform right-circularly polarized (RCP) incident light into a different PVB with LCP state. Theoretically, the device must provide two independent spatial phase profiles corresponding to LCP and RCP, which cannot be achieved by using a single geometric phase or propagation phase.
For an anisotropic unit cell, the corresponding Jones matrix is:
T ( x , y ) = R ( θ ) T R ( θ ) = ( t x cos 2 θ + t y sin 2 θ t y sin θ cos θ t x sin θ cos θ t y sin θ cos θ t x sin θ cos θ t x cos 2 θ + t y sin 2 θ )
where t x = T x e i φ x , t y = T y e i φ y   are complex amplitude, T i and φ i represent the transmission amplitude and the phase delay while the incident light is along with the i -polarized wave with i ϵ ( x , y ) , respectively. R ( θ ) = [ cos θ sin θ sin θ cos θ ]   is the rotation matrix with the rotation angle   θ . When the metasurface is illuminated by circularly polarized light as E i n = 2 2 [ 1 ± i ] , the output electric field can be expressed as:
E o u t = T ( x , y ) E i n = 2 4 ( t x + t y ) [ 1 ± i ] + 2 4 ( t x t y ) e ± i 2 θ [ 1 i ]
Here   [ 1 i ] and [ 1 i ]   represent LCP and RCP, respectively. The transmitted light can be divided into two parts: the co-polarized component without any phase regulation and the cross-polarized component with conjugated phase modulation ± 2 θ . Combining the orientation-dependent geometric phase with the dimension-dependent propagation phase effectively will address the spin-locked limitation of geometry metasurfaces [51,52]. Consider the unit cell is lossless ( T x = T y = 1 ) with π phase difference ( | φ x ( x , y ) φ y ( x , y ) = π | ) , indicating the pillars working as a half waveplate. The output electric field can then be expressed as:
[ E o u t x E o u t y ] = 2 2 e φ x e ± i 2 θ [ 1 i ]
Assuming the expected phase distributions of the spin-decoupled are φ 1 ( x , y )   and φ 2 ( x , y ) . From the Equation (9), we can derive that:
{ φ x ( x , y ) + 2 θ ( x , y ) = φ 1 ( x , y ) φ x ( x , y ) 2 θ ( x , y ) = φ 2 ( x , y )
where φ x ( x , y )   represents the propagation phase, which is influenced by the material and geometry parameters of the unit pillars; ± 2 θ ( x , y ) is the geometric phase, which is only determined by the rotation angle of the unit cell. Therefore, the expected phase distributions φ 1 ( x , y )   and φ 2 ( x , y ) can be satisfied simultaneously by designing the geometry and rotation angle of the unit structure, as shown in Figure 3b.
To construct the metasurface, the sub-wavelength rectangular-shaped silicon pillar has been chosen, as shown in Figure 4a. The simulation is based on the FDTD method to obtain the pillars’ transmittance and phase delay database under the x- and y-polarized incident waves. The length L and width W varied from 35 μm to 135 μm with the fixed lattice constant of P = 140 μm and the height of H = 330 μm at the working frequency of 0.75 THz. Figure 4b,c show the transmittance as the function of the pillar sizes in L and W, respectively. The corresponding simulated phase shifts ( φ x ,   φ y ) are illustrated in Figure 4d,e. Both φ x   and φ y cover 2π degrees in the scanning range, and all the transmittances are more than 70%. Therefore, by carefully selecting L and W from the database, a nearly arbitrary phase combination ( φ x ,   φ y ) with high transmittance can be obtained for the silicon pillar. For our design, nine different silicon pillars are chosen, and the corresponding transmittances, phase shifts, and polarization conversion efficiencies (PCE, defined as the intensity ratio between the converted spin component and the total transmitted beam) are shown in Figure 4f. Phase shifts φ x (red square) and φ y   (green square) are linearly increased with a step of 2π/9. The corresponding transmittance T x (red sphere) and T y (green sphere) are more than 70%. Moreover, the PCE (orange column) is maintained above 70%. Most importantly, the phase differences between φ x and φ y are close to π, and the transmittance T x and T y are almost equal, indicating these selected pillars are working well as the half-wave plates at 0.75 THz. The specific geometric parameters and PCE of the nine selected silicon pillars are listed in Table 1.

2.3. Superpositions of PVBs by a Spin-Decoupled Metasurface

Under the paraxial approximation, the vectorial field of a monochromatic Poincaré beam can be decomposed into the superposition of two orthogonal circularly polarized vortex beams with different coefficients ψ N l and ψ s m as described [22]:
| ψ l , m = ψ N l | R l + ψ S m | L m
where | R l = e i l φ ( e x + i e y ) / 2 , | L m = e i m φ ( e x i e y ) / 2 are the RCP and LCP vortex beams with topological charges of l and m, respectively, while ex and ey are the unit vectors along the x and y axes. φ is the azimuthal angle in the polar coordinates. In order to obtain PPBs that radius is irrelevant to topological charges, Equation (11) can be further deduced as:
| ψ l , m = ψ N l | P O V R , l + ψ S m | P O V L , m
Here, | P O V R , l = e i φ t o t a l ( e x + i e y ) / 2 and | P O V L , m = e i φ t o t a l ( e x i e y ) / 2 represents the RCP and LCP PVBs with topological charges of l and m, φ t o t a l   represent the phase profile of PVBs obtained from Equation (3).
According to Equation (12), two essential parts, orthogonal circularly polarized PVBs and coefficients ψ N   l and ψ s m , should be resolved for generating arbitrary PPBs. In order to describe the process more clearly, we use the PS to describe scalar fields of incident light and the HyOPS to describe the vectorial field of output PPBs. As shown in Figure 5a, all the possible states of fundamental polarization can be mapped onto the PS. Because the polarization state of the scalar field can be represented with orthogonal circular polarization, the coefficients ψ N l and ψ s m vary with the incident polarization. As shown in Figure 5b, the spin-decoupled metasurface can generate two completely different PVBs with orthogonal circularly polarization. If the axicon period and the focal length are identical, arbitrary PPBs can be generated by superposing the two circularly polarized PVBs on a HyOPS. As shown in Figure 5c, each point on the surface of HyOPS with spherical coordinates ( α , β ) can indicate a PPB, where α ( 0 , 2 π ) , β ( 0 , π ) . The two poles of the HyOPS represent two PVBs with opposite circular polarization and topological charges of l and m , respectively. Points on the equator of HyOPS represent a superposition of equal intensities and different phase differences of the opposite circular polarization states. For example, the points   ( 0 , π / 2 ) , ( π , π / 2 ) represent the radially and azimuthally vector beam, respectively.
For the PPBs on a HyOPS, the polarization state can be mapped by representing the Stokes parameters in the spherical Cartesian coordinates, which are defined as follows [20]:
S 0 l , m = | ψ N l | 2 + | ψ S m | 2
S 1 l , m = 2 | ψ N l | | ψ S m | 2 cos Φ
S 2 l , m = 2 | ψ N l | | ψ S m | 2 sin Φ
S 3 l , m = | ψ N l | 2 | ψ S m | 2
Here Φ = arg ( ψ N l ) ) arg ( ψ N l ) is the phase difference between the two PVBs. | ψ N l | 2 and | ψ S m | 2 are the intensity of RCP and LCP PVB, respectively. The concepts of polarization order P = ( l m ) / 2 and the topological Pancharatnam charge l P = ( l + m ) / 2 are also used to characterize PPBs [22]. This can be explained by deducting Equation (12) as:
| ψ l , m = e i l p φ ( ψ N l e i p φ | P O V R , l + ψ S m e i p φ | P O V L , m )
Equation (17) indicates that the polarization state of PPBs is only determined by the parameters within the brackets, that is, the polarization order P and the HyOPS coordinate ( α , β ) . The topological Pancharatnam charge l P characterizes the phase distribution of the PPBs, which can be represented as e i φ l P .

3. Results and Discussion

Based on the above design principles, we designed and numerically simulated several metasurface devices (diameter is 14 mm) to demonstrate the practicability of our proposed design by employing the FDTD method. The perfectly matching layers (PMLs) were used in the x, y, and z directions. Plane-wave sources were used in all simulations.

3.1. Polarization-Independent PVBs Generator

The generation of polarization-independent PVBs has been investigated. The design phase distribution has been described in Equation (3). Figure 6a displays the intensity profiles in the x-z plane of the generated PVB with the parameters of d = 4 mm and l = 2 at 0.75 THz. Figure 6b plots the corresponding cross-section images of the generated PVB at different longitudinal positions. The THz beam passing through the metasurface presents an annular ring profile near the focal plane (z = 14 mm). Moreover, the intensity profile of the ring maintains a circular shape, but the radius increases with increasing propagation distance.
Figure 7a depicts the normalized intensity and phase distribution of the generated PVBs with different topological charges from l = 1 to l = 4 at the focus point when d = 3 mm and 4 mm, respectively. The focal length of the Fourier- transformed lens is designed as f = 14 mm at 0.75 THz. Figure 7b shows the corresponding cross-sections of the normalized intensity profiles extracted along the x-direction (solid line) and y-direction (dashed line). The consistent cross-sectional intensity profiles show that the generated vortex beams are insensitive to the topological charge, and the phase distribution is consistent with the designed topological charge. The full-wave simulation generation efficiencies (defined as the ratio of the PVB’s intensity to the intensity of the incident beam) of all polarization-independent metasurfaces are more than 49%. The polarization-independent properties and the relationship between the ring radius and the axicon period are analyzed in Sections S1 and S2 of the Supplemental Materials. Section S3 of the Supplemental Materials characterizes the working bandwidth of the polarization-independent PVBs generator.

3.2. Spin-Decoupled PVBs Generator

According to theoretical analysis of a spin-decoupled PVBs generator, a single metasurface can generate two distinct PVBs. The corresponding phase distribution ( φ 1 ,   φ 2 ) can be expressed as follows:
{ φ 1 ( x , y ) = l 1 arctan ( x y ) + ( 2 π x 2 y 2 d 1 ) + π ( x 2 + y 2 ) λ f 1 φ 2 ( x , y ) = = l 2 arctan ( x y ) + ( 2 π x 2 y 2 d 2 ) + π ( x 2 + y 2 ) λ f 2
where λ is the incident wavelength, li, fi, and di represent corresponding parameters for the incident RCP (i = 1) and LCP (i = 2) light, respectively. Two metasurface devices are designed based on Equations (10) and (18) to demonstrate the above function. The first metasurface is designed with l1 = 2, l2 = 3, d1 = d2 = 4 mm, and f1 = f2 = 14 mm at 0.75 THz. As shown in Figure 8a,c, when the incident light is converted from RCP to LCP, the transmitted light is converted into a corresponding cross-polarized component, and the ring radius does not change significantly, but the topological charge l changes from 2 to 3, which is consistent with our design. In addition, we can also control the ring radius of generated PVBs by switching the incident spin state. The second device with the parameters l1 = l2 = 3, d1 = 4 mm, d2 = 6 mm, and f1 = f2 =14 mm at 0.75 THz is also designed. As shown in Figure 8b,d, the topological charge value l does not change, but the ring radius is inversely proportional to the axicon periods d when the incident light is converted from RCP to LCP. It is worth noting that the transmitted light still carries with a co-polarized component, as shown in Figure 8a,b, because the selected unit cells are not a perfect half-wave plate. The full-wave simulation generation efficiencies of the spin-decoupled metasurfaces were more than 53%, which can be improved by selecting materials with higher transmittance and a more careful phase scanning operation with smaller step size (more details on working bandwidth of the meta-device seen also in Section S4 of the Supplemental Materials). Therefore, the spin-decoupled and polarization-independent PVB meta-devices will pave the way for developing a miniaturized and convenient platform for OAM-related applications in the terahertz band, especially in improving the transmission capacity of terahertz communications.

3.3. Superpositions of Perfect Vortex Beams

According to the theoretical analysis of generating PPBs in Section 2.3, the topological charges (l, m) of orthogonal circularly polarized PVBs and the corresponding coefficients ψ N l and ψ s m can determine a PPB. Specifically, different topological charges (l, m) satisfied by spin-decoupled metasurface will determine a unique HyOPS, different coefficients ψ N l and ψ s m   will define the position of a PPB on this HyOPS by calculating Stokes parameters. In order to verify the above principle, three spin-decoupled metasurfaces with the expected parameters d = 4 mm, and f =14 mm at 0.75 THz have been designed. The relevant polarization order P and topological Pancharatnam charge l P are P 1 = ( l 1 m 1 ) / 2 = 1 , l P 1 = ( l 1 + m 1 ) / 2 =0, P 2 = ( l 2 m 2 ) / 2 = 2 , l P 2 = ( l 2 + m 2 ) / 2   = 0, and P 3 = ( l 3 m 3 ) / 2 = 1 , l P 3 = ( l 3 + m 3 ) / 2 = 2 . Five points with their coordinates on the HyOPS are selected, given by the red dots in Figure 9a, and the corresponding incident polarization states are represented by the red points located on the PS shown in Figure 9b. Figure 9c shows the intensity patterns of three designed metasurfaces in the x-y plane at the focus point. The annular intensity patterns are obtained by a horizontal linear polarizer to analyze the polarization property. Because the polarization order P determines the number of polarization rotations per round trip, the transmitted intensity patterns through a horizontal linear polarizer will have   2 | P | lobes.
The case with the same topological Pancharatnam charge l P and different polarization order P is firstly discussed. As shown in the first row and the second row in Figure 9c, the transmitted intensity patterns are split into two and four lobes, respectively, which satisfy the theoretical predictions of P 1 = 1 , P 2 = 2 . Then, the influence of the topological Pancharatnam charge l P is considered. The polarization order is fixed as P = 1 , and varied with the topological Pancharatnam charge l P = 0 ,   2 , respectively. The simulated annular intensity distributions in the first and third row show the same annular patterns, demonstrating the independence of the polarization state of PPBs from the topological Pancharatnam charge, matching well with the theoretical prediction of Equation (17). It should be noticed that the third row has a slight rotation of the measured intensity induced by the Gouy phase [53]. In principle, the Gouy phase is dependent on the topological charge, and the rotation of the intensity distribution around the beam axis occurs only when the topological charges | l | | m | . Therefore, this approach provides a potential solution to measure the Gouy phase.
To further characterize the capability of arbitrary PPBs generation, another five incident lights represented by the blue points on the Poincaré sphere are chosen, as shown in Figure 9b, where the two circular eigenstates have a phase difference of π/2. The corresponding PPBs with their coordinates on the HyOPS are expressed by the blue dots in Figure 9a, and the measured intensity patterns through a horizontal linear polarizer are shown in Figure 9d. Unsurprisingly, the introduced phase difference of π/2 between the two circular eigenstates results in the change of coefficients ψ N l and ψ s m , further causing the rotation of the intensity profiles and demonstrating the feasibility of PPBs by changing the polarization of incident light.
In addition, the Stokes parameters have also been used to verify inhomogeneous polarization states of the PPBs, which are given by [54]:
S 0 = I 0 + I 90
S 1 = I 0 I 90
S 2 = I 45 I 135
S 3 = I R I L
where I 0 , I 90 , I 45 and I 135 are the intensities of PPBs through a linear polarizer oriented at 0°, 45°, 90°, and 145°. I R and I L are the intensities of PPBs through a right and a left circular polarizer, which can be achieved by cascading a linear polarizer and a quarter waveplate. With these Stokes parameters, the detailed polarization of a PPB can be obtained. The spherical coordinates ( α , β ) of the points on the HyOPS can be determined by the Stokes parameters:
α = a r c t a n ( S 2 / S 1 )
α = a r c c o s ( S 2 / S 1 )
Figure 10a shows the calculated results of normalized Stokes parameters by choosing point III in Figure 9a for three kinds of PPBs. The first Stokes parameter S0, representing the light intensity, shows the same ring pattern consistent with our theoretical analysis. By analyzing Stokes parameters, S1 and S2, the polarization orientation angle σ can be determined by:
σ = 1 2 arctan ( S 2 S 1 )
The corresponding polarization orientations (orange double arrows) and distributions of the PPBs are shown in Figure 10b. It can be seen that the polarization orientations rotate 2π and 4π per round trip for cases of P 1 = P 3 = 1 , P 2 = 2 , respectively, which is consistent to the definition of polarization order. It should be noticed that the spherical coordinate of point III is ( 0 , π / 2 ) , which means the Stokes parameter S3 is zero theoretically. The error appearing in Figure 10a mainly comes from two reasons: the phase difference of the metasurface between the actual phase and the theoretical phase, and the act that the selected unit cells are not a perfect half-wave plate. The full-wave simulation generation efficiencies of the PPBs (defined as the ratio of the PPB’s intensity to the intensity of the incident beam) were more than 51%. These results and the analysis demonstrate that spin-decoupled metasurfaces can efficiently generate arbitrary PPBs on the HyOPS in THz region. Therefore, the terahertz PVB meta-devices designed here will also provide a potential compact platform for vector beam related applications in the terahertz band, such as reconstructing incident photons’ frequency and polarization state [55].

4. Conclusions

In conclusion, we propose a general method to generate PVBs using a single-layer dielectric metasurface in THz regime. This is achieved by integrating the function of spiral phase plate, axicon, and Fourier transformation lens into a single-layer metasurface. The ring radiuses of these generated vortex beams are independent of the topological charges and inversely proportional to the axicon periods. Meanwhile, the unit-cell cross-section geometries with 4-fold symmetry can make the generated PVBs independent of incident polarization. Moreover, to realize more freedom, a spin-decoupled metasurface is designed by simultaneously manipulating the dynamic phase and geometric phase of the unit cells. The metasurface can generate spin-decoupled PVBs with arbitrary radius or topological charges by changing the incident spin states. In addition, the generation of an arbitrary perfect Poincaré beam at any position on the surface of the HyOPS is also demonstrated by superposing two perfect vortex beams. Various perfect Poincaré beams were generated, and the results agreed well with the theory. This work will provide more opportunities and possibilities for THz communication, complex light field generation, and quantum information sciences.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12173010/s1, Figure S1: (a) The intensity of the PVBs at the focus point under different polarized incident light (depicted by the red double arrow). (b) Corresponding normalized cross-section of the annular intensity distribution along x-direction (solid line) and y-direction (dashed line). (c) The relationship between the ring radius and the axicon period with the unified parameters l = 2, f = 14 mm. Figure S2. Characterization of the bandwidth of the polarization-independent PVBs generator. (a) The intensity of the PVBs at the focus point with a frequency ranging from 0.65 THz to 0.9 THz. (b–g) Corresponding normalized cross-section of the annular intensity distribution along x-direction (blue line) and y-direction (red line). Figure S3. Analysis of the bandwidth of the spin-decoupled PVBs generator. (a) The intensity of the PVBs at the focus point with a frequency ranging from 0.5 THz to 1.0 THz. (b–g) Corresponding normalized cross-section of the annular intensity distribution along x-direction (blue line) and y-direction (red line).

Author Contributions

Conceptualization, Q.W.; methodology, Q.W.; software, Q.W.; validation, Q.W. and W.F.; formal analysis, Q.W., C.Q. and W.F.; investigation, Q.W.; resources, W.F.; data curation, Q.W.; writing—original draft preparation, Q.W.; writing—review and editing, Q.W., C.Q. and W.F.; visualization, Q.W.; supervision, W.F.; project administration, W.F.; funding acquisition, W.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (No. 61675230).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allen, L.; Beijersbergen, M.W.; Spreeuw, R.J.; Woerdman, J.P. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Phys. Rev. A 1992, 45, 8185–8189. [Google Scholar] [CrossRef]
  2. Cai, X.L.; Wang, J.W.; Strain, M.J.; Johnson-Morris, B.; Zhu, J.B.; Sorel, M.; O’Brien, J.L.; Thompson, M.G.; Yu, S.T. Integrated compact optical vortex beam emitters. Science 2012, 338, 363–366. [Google Scholar] [CrossRef]
  3. Vallone, G.; D’Ambrosio, V.; Sponselli, A.; Slussarenko, S.; Marrucci, L.; Sciarrino, F.; Villoresi, P. Free-space quantum key distribution by rotation-invariant twisted photons. Phys. Rev. Lett. 2014, 113, 060503. [Google Scholar] [CrossRef]
  4. Ng, J.; Lin, Z.; Chan, C.T. Theory of optical trapping by an optical vortex beam. Phys. Rev. Lett. 2010, 104, 103601. [Google Scholar] [CrossRef]
  5. Lavery, M.P.; Speirits, F.C.; Barnett, S.M.; Padgett, M.J. Detection of a spinning object using light’s orbital angular momentum. Science 2013, 341, 537–540. [Google Scholar] [CrossRef]
  6. Parigi, V.; D’Ambrosio, V.; Arnold, C.; Marrucci, L.; Sciarrino, F.; Laurat, J. Storage and retrieval of vector beams of light in a multiple-degree-of-freedom quantum memory. Nat. Commun. 2015, 6, 7706. [Google Scholar] [CrossRef]
  7. Ferguson, B.; Zhang, X.C. Materials for terahertz science and technology. Nat. Mater. 2002, 1, 26–33. [Google Scholar] [CrossRef]
  8. Tonouchi, M. Cutting-edge terahertz technology. Nat. Photonics 2007, 1, 97–105. [Google Scholar] [CrossRef]
  9. Nagatsuma, T.; Ducournau, G. Advances in terahertz communications accelerated by photonics. Nat. Photonics 2016, 10, 371–379. [Google Scholar] [CrossRef]
  10. Wei, X.; Zhu, L.; Zhang, Z.; Wang, K.; Liu, J.; Wang, J. Orbit Angular Momentum Multiplexing in 0.1-THz Free-Space Communication via 3D Printed Spiral Phase Plates; CLEO; Optica Publishing Group: San Jose, CA, USA, 2014; p. STu2F.2. [Google Scholar]
  11. Wu, G.B.; Chan, K.F.; Shum, K.M.; Chan, C.H. Millimeter-wave and terahertz OAM discrete-lens antennas for 5G and beyond. IEEE Commun. Mag. 2022, 60, 34–39. [Google Scholar] [CrossRef]
  12. Yan, H.W.; Zhang, E.T.; Zhao, B.Y.; Duan, K.L. Free-space propagation of guided optical vortices excited in an annular core fiber. Opt. Express 2012, 20, 17904–17915. [Google Scholar] [CrossRef]
  13. Li, H.S.; Ren, G.B.; Zhu, B.F.; Gao, Y.X.; Yin, B.; Wang, J.; Jian, S.S. Guiding terahertz orbital angular momentum beams in multimode Kagome hollow-core fibers. Opt. Lett. 2017, 42, 179–182. [Google Scholar] [CrossRef]
  14. Xie, X.; Chen, Y.; Yang, K.; Zhou, J. Harnessing the point-spread function for high-resolution far-field optical microscopy. Phys. Rev. Lett. 2014, 113, 263901. [Google Scholar] [CrossRef]
  15. Ramachandran, S.; Kristensen, P.; Yan, M.F. Generation and propagation of radially polarized beams in optical fibers. Opt. Lett. 2009, 34, 2525–2527. [Google Scholar] [CrossRef]
  16. Fang, X.Y.; Ren, H.R.; Gu, M. Orbital angular momentum holography for high-security encryption. Nat. Photonics 2020, 14, 102–108. [Google Scholar] [CrossRef]
  17. Lissberger, P.H. Ellipsometry and polarised Light. Nature 1977, 269, 270. [Google Scholar] [CrossRef]
  18. Milione, G.; Sztul, H.I.; Nolan, D.A.; Alfano, R.R. Higher-order Poincaré sphere, stokes parameters, and the angular momentum of Light. Phys. Rev. Lett. 2011, 107, 053601. [Google Scholar] [CrossRef]
  19. Yi, X.; Liu, Y.; Ling, X.; Zhou, X.; Ke, Y.; Luo, H.; Wen, S.; Fan, D. Hybrid-order Poincare sphere, Phys. Rev. A 2015, 91, 023801. [Google Scholar] [CrossRef]
  20. Galvez, E.J. Light beams with spatially variable polarization. In Photonics—Fundamentals of Photonics and Physics; Wiley: Hoboken, NJ, USA, 2015; Volume 1, pp. 61–76. [Google Scholar]
  21. Han, J.; Intaravanne, Y.; Ma, A.N.; Wang, R.X.; Li, S.T.; Li, Z.C.; Chen, S.Q.; Li, J.S.; Chen, X.Z. Optical metasurfaces for generation and superposition of optical ring vortex beams. Laser Photonics Rev. 2022, 14, 2000146. [Google Scholar] [CrossRef]
  22. Niv, A.; Biener, G.; Kleiner, V.; Hasman, E. Manipulation of the Pancharatnam phase in vectorial vortices. Opt. Express 2006, 14, 4208–4220. [Google Scholar] [CrossRef]
  23. Philip, G.M.; Kumar, V.; Milione, G.; Viswanathan, N.K. Manifestation of the Gouy phase in vector-vortex beams. Opt. Lett. 2012, 37, 2667–2669. [Google Scholar] [CrossRef] [PubMed]
  24. Ostrovsky, A.S.; Rickenstorff-Parrao, C.; Arrizón, V. Generation of the “perfect” optical vortex using a liquid-crystal spatial light modulator, Opt. Lett. 2013, 38, 534–536. [Google Scholar]
  25. Fu, S.; Gao, C.; Wang, T.; Zhang, S.; Zhai, Y. Simultaneous generation of multiple perfect polarization vortices with selective spatial states in various diffraction orders. Opt. Lett. 2016, 41, 5454–5457. [Google Scholar] [CrossRef]
  26. Vaity, P.; Rusch, L. Perfect vortex beam: Fourier transformation of a Bessel beam. Opt. Lett. 2015, 40, 597–600. [Google Scholar] [CrossRef]
  27. Jabir, M.V.; Apurv Chaitanya, N.; Aadhi, A.; Samanta, G.K. Generation of “perfect” vortex of variable size and its effect in angular spectrum of the down-converted photons. Sci. Rep. 2016, 6, 21877. [Google Scholar] [CrossRef] [PubMed]
  28. Li, P.; Zhang, Y.; Liu, S.; Ma, C.; Han, L.; Cheng, H.; Zhao, J. Generation of perfect vectorial vortex beams. Opt. Lett. 2016, 41, 2205–2208. [Google Scholar] [CrossRef]
  29. Guo, Z.; Liu, H.; Xiang, L.; Chen, L.; Yang, J.; Wen, J.; Shang, Y.; Wang, T.; Pang, F. Generation of perfect vortex beams with polymer-based phase plate. IEEE Photon. Technol. Lett. 2020, 32, 565–568. [Google Scholar] [CrossRef]
  30. Gu, Z.; Yin, D.; Gu, F.; Zhang, Y.; Nie, S.; Feng, S.; Ma, J.; Yuan, C. Generation of concentric perfect Poincaré beams. Sci. Rep. 2019, 9, 15301. [Google Scholar] [CrossRef]
  31. Li, D.; Feng, S.; Nie, S.S.; Chang, C.; Ma, J.; Yuan, C. Generation of arbitrary perfect Poincaré beams. J. Appl. Phys. 2019, 125, 073105. [Google Scholar] [CrossRef]
  32. Kildishev, A.V.; Boltasseva, A.; Shalaev, V.M. Planar photonics with metasurfaces. Science 2013, 339, 1232009. [Google Scholar] [CrossRef]
  33. Qiu, C.W.; Zhang, T.; Hu, G.; Kivshar, Y. Quo Vadis, Metasurfaces? Nano Lett. 2021, 21, 5461–5474. [Google Scholar] [CrossRef] [PubMed]
  34. Rubin, N.A.; Shi, Z.; Capasso, F. Polarization in diffractive optics and metasurfaces. Adv. Opt. Photon. 2021, 13, 836–970. [Google Scholar] [CrossRef]
  35. Němec, H.; Kuel, P.; Kadlec, F.; Kadlec, C.; Mounaix, P. Tunable terahertz metamaterials with negative permeability. Phys. Rev. B 2009, 79, 241108. [Google Scholar] [CrossRef]
  36. Kowerdziej, R.; Jaroszewicz, L.; Olifierczuk, M.; Parka, J. Experimental study on terahertz metamaterial embedded in nematic liquid crystal. Appl. Phys. Lett. 2015, 106, 092905. [Google Scholar] [CrossRef]
  37. Wang, S.; Kang, L.; Werner, D.H. Hybrid resonators and highly tunable terahertz metamaterials enabled by vanadium dioxide (VO2). Sci. Rep. 2017, 7, 4326. [Google Scholar] [CrossRef]
  38. Zhang, Y.C.; Liu, W.W.; Gao, J.; Yang, X.D. Generating focused 3D perfect vortex beams by plasmonic metasurfaces. Adv. Opt. Mate. 2018, 6, 1701228. [Google Scholar] [CrossRef]
  39. Xie, J.; Guo, H.; Zhuang, S.; Hu, J. Polarization-controllable perfect vortex beam by a dielectric metasurface. Opt. Express 2021, 29, 3081–3089. [Google Scholar] [CrossRef]
  40. Bao, Y.; Ni, J.; Qiu, C.-W. A minimalist single-layer metasurface for arbitrary and full control of vector vortex beams. Adv. Mater. 2020, 32, 1905659. [Google Scholar] [CrossRef]
  41. Zhou, Q.; Liu, M.; Zhu, W.; Chen, L.; Ren, Y.; Lezec, H.J.; Lu, Y.; Agrawal, A.; Xu, T. Generation of perfect vortex beams by dielectric geometric metasurface for visible light. Laser Photonics Rev. 2021, 15, 2100390. [Google Scholar] [CrossRef]
  42. Xu, Y.H.; Zhang, H.F.; Li, Q.; Zhang, X.Q.; Xu, Q.; Zhang, W.T.; Hu, C.; Zhang, X.X.; Han, J.G.; Zhang, W.L. Generation of terahertz vector beams using dielectric metasurfaces via spin-decoupled phase control. Nanophotonics 2020, 9, 3393–3402. [Google Scholar] [CrossRef]
  43. Liu, W.; Yang, Q.; Xu, Q.; Jiang, X.; Wu, T.; Gu, J.; Han, J.; Zhang, W. Multichannel terahertz quasi-perfect vortex beams generation enabled by multifunctional metasurfaces. Nanophotonics 2022, 11, 3631–3640. [Google Scholar] [CrossRef]
  44. Liu, Y.; Ke, Y.; Zhou, J.; Liu, Y.; Luo, H.; Wen, S.; Fan, D. Generation of perfect vortex and vector beams based on Pancharatnam-Berry phase elements. Sci. Rep. 2017, 7, 44096. [Google Scholar] [CrossRef]
  45. Wang, Q.; Zhang, X.; Xu, Y.; Tian, Z.; Gu, J.; Yue, W.; Zhang, S.; Han, J.; Zhang, W. A broadband metasurface-based terahertz flat-lens array. Adv. Opt. Mater. 2015, 3, 779–785. [Google Scholar] [CrossRef]
  46. Jia, M.; Wang, Z.; Li, H.; Wang, X.; Luo, W.; Sun, S.; Zhang, Y.; He, Q.; Zhou, L. Efficient manipulations of circularly polarized terahertz waves with transmissive metasurfaces. Light Sci. Appl. 2019, 8, 16. [Google Scholar] [CrossRef]
  47. Shalaginov, M.Y.; An, S.; Zhang, Y.; Yang, F.; Su, P.; Liberman, V.; Chou, J.B.; Roberts, C.M.; Kang, M.; Rios, C.; et al. Reconfigurable all-dielectric metalens with diffraction-limited performance. Nat. Commun. 2021, 12, 1225. [Google Scholar] [CrossRef]
  48. Zhang, H.; Zhang, X.; Xu, Q.; Wang, Q.; Xu, Y.; Wei, M.; Li, Y.; Gu, J.; Tian, Z.; Ouyang, C.; et al. Polarization-independent all-silicon dielectric metasurfaces in the terahertz regime. Photon. Res. 2018, 6, 24–29. [Google Scholar] [CrossRef]
  49. 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. Photon. Res. 2021, 9, 1939–1947. [Google Scholar] [CrossRef]
  50. Arbabi, A.; Horie, Y.; Ball, A.J.; Bagheri, M.; Faraon, A. A subwavelength-thick lenses with high numerical apertures and large efficiency based on high-contrast transmit arrays. Nat. Commun. 2015, 6, 7069. [Google Scholar] [CrossRef]
  51. Yin, X.; Ye, Z.; Rho, J.; Wang, Y.; Zhang, X. Photonic spin Hall effect at metasurfaces. Science 2013, 339, 1405–1407. [Google Scholar] [CrossRef]
  52. Balthasar Mueller, J.P.; Rubin, N.A.; Devlin, R.C.; Groever, B.; Capasso, F. Metasurface polarization optics: Independent phase control of arbitrary orthogonal states of polarization. Phys. Rev. Lett. 2017, 118, 113901. [Google Scholar] [CrossRef]
  53. Baumann, S.M.; Kalb, D.M.; MacMillan, L.H.; Galvez, E.J. Propagation dynamics of optical vortices due to Gouy phase. Opt. Express 2009, 17, 9818–9827. [Google Scholar] [CrossRef]
  54. Schaefer, B.; Collett, E.; Smyth, R.; Barrett, D.; Fraher, B. Measuring the stokes polarization parameters. Am. J. Phys. 2007, 75, 163–168. [Google Scholar] [CrossRef]
  55. Yu, F.; Chen, J.; Huang, L.; Zhao, Z.; Wang, J.; Jin, R.; Chen, J.; Wang, J.; Miroshnichenko, A.E.; Li, T.; et al. Photonic slide rule with metasurfaces. Light Sci. Appl. 2022, 11, 77. [Google Scholar] [CrossRef]
Figure 1. Generation of polarization-independent PVBs. (a) Schematic of the polarization-independent metasurface. (b) Simulated intensity distribution and phase distribution of PVBs. (c) The phase profile of the designed metasurface includes the phase of a spiral phase plate, an axicon, and a Fourier transformation lens.
Figure 1. Generation of polarization-independent PVBs. (a) Schematic of the polarization-independent metasurface. (b) Simulated intensity distribution and phase distribution of PVBs. (c) The phase profile of the designed metasurface includes the phase of a spiral phase plate, an axicon, and a Fourier transformation lens.
Nanomaterials 12 03010 g001
Figure 2. Characterization of the designed unit cell for modulating the propagation phase. (a) Schematic of the unit cell: the lattice constant P, the unit length L, the width W, and the height H. (b) Transmittance and phase shifts as a function of the unit length L. The inset illustrates the side views of the simulated magnetic amplitude distributions of the unit cell with L = 120 μm. (c) The phase shifts and transmittance of eight selected silicon units. (d) Simulated phase wavefront in three silicon units with L = 40, 82, and 120 μm, respectively.
Figure 2. Characterization of the designed unit cell for modulating the propagation phase. (a) Schematic of the unit cell: the lattice constant P, the unit length L, the width W, and the height H. (b) Transmittance and phase shifts as a function of the unit length L. The inset illustrates the side views of the simulated magnetic amplitude distributions of the unit cell with L = 120 μm. (c) The phase shifts and transmittance of eight selected silicon units. (d) Simulated phase wavefront in three silicon units with L = 40, 82, and 120 μm, respectively.
Nanomaterials 12 03010 g002
Figure 3. Generation of spin-decoupled PVBs. (a) Schematic of the spin-decoupled metasurface. (b) The designed total phase of the single-layer metasurface simultaneously contains the geometric phase and the dynamic phase, which are only related to geometry parameters and the rotation angle of the unit structure, respectively.
Figure 3. Generation of spin-decoupled PVBs. (a) Schematic of the spin-decoupled metasurface. (b) The designed total phase of the single-layer metasurface simultaneously contains the geometric phase and the dynamic phase, which are only related to geometry parameters and the rotation angle of the unit structure, respectively.
Nanomaterials 12 03010 g003
Figure 4. (a) Schematic of the rectangular-shape silicon pillar: the lattice constant P, the length L, the width W, and the height H, respectively. The angle θ represents the rotation of the unit cell. (be) Simulated transmittance and phase shift as a pillar sizes (L, W) function under x-polarized and y-polarized illumination at 0.75 THz, respectively. (f) Transmittances, phase shifts, and polarization conversion efficiencies of the nine selected unit cells.
Figure 4. (a) Schematic of the rectangular-shape silicon pillar: the lattice constant P, the length L, the width W, and the height H, respectively. The angle θ represents the rotation of the unit cell. (be) Simulated transmittance and phase shift as a pillar sizes (L, W) function under x-polarized and y-polarized illumination at 0.75 THz, respectively. (f) Transmittances, phase shifts, and polarization conversion efficiencies of the nine selected unit cells.
Nanomaterials 12 03010 g004
Figure 5. Illustration of generating PPBs through a dielectric metasurface with the aid of PS and HyOPS. (a) Schematic of Poincaré sphere. A point on the Poincaré sphere represents the polarization state of the incident light. (b) Schematic of the spin-decoupled metasurface that can superpose two completely different PVBs. (c) The output PPBs is represented by points on the HyOPS. The inset illustrates the radially and azimuthally polarized beams’ intensity and polarization distribution (depicted by white arrows) on the HyOPS with l = −1 and m = 1.
Figure 5. Illustration of generating PPBs through a dielectric metasurface with the aid of PS and HyOPS. (a) Schematic of Poincaré sphere. A point on the Poincaré sphere represents the polarization state of the incident light. (b) Schematic of the spin-decoupled metasurface that can superpose two completely different PVBs. (c) The output PPBs is represented by points on the HyOPS. The inset illustrates the radially and azimuthally polarized beams’ intensity and polarization distribution (depicted by white arrows) on the HyOPS with l = −1 and m = 1.
Nanomaterials 12 03010 g005
Figure 6. (a) The intensity profiles in the x-z plane of the generated PVB with the parameters of d = 4 mm and l = 2. (b) Corresponding 2D cross-section images of the generated PVB at different longitudinal positions.
Figure 6. (a) The intensity profiles in the x-z plane of the generated PVB with the parameters of d = 4 mm and l = 2. (b) Corresponding 2D cross-section images of the generated PVB at different longitudinal positions.
Nanomaterials 12 03010 g006
Figure 7. (a) The intensity and phase profiles of the PVBs with the topological charges varied from l = 1 to l = 4 when the axicon periods d were selected as d = 3 mm and 4 mm, respectively. (b) Corresponding normalized cross-section of the annular intensity distribution along x-direction (solid line) and y-direction (dashed line) at 0.75 THz when d = 3 mm (upper) and d = 4 mm (down), respectively.
Figure 7. (a) The intensity and phase profiles of the PVBs with the topological charges varied from l = 1 to l = 4 when the axicon periods d were selected as d = 3 mm and 4 mm, respectively. (b) Corresponding normalized cross-section of the annular intensity distribution along x-direction (solid line) and y-direction (dashed line) at 0.75 THz when d = 3 mm (upper) and d = 4 mm (down), respectively.
Nanomaterials 12 03010 g007
Figure 8. (a,b) The normalized intensity profile and phase distribution for devices I and II at focus point when only the RCP or LCP beam is incident. (c,d) Corresponding normalized cross-section of the annular intensity distribution for RCP (blue) and LCP (red) polarization along x-direction (solid line) and y-direction (dashed line).
Figure 8. (a,b) The normalized intensity profile and phase distribution for devices I and II at focus point when only the RCP or LCP beam is incident. (c,d) Corresponding normalized cross-section of the annular intensity distribution for RCP (blue) and LCP (red) polarization along x-direction (solid line) and y-direction (dashed line).
Nanomaterials 12 03010 g008
Figure 9. Characterization of the generated PPBs. (a) The selected states of PPBs are represented by red and blue points on the HyOPS with their coordinates. The phase difference between the two circular eigenstates of the selected blue dots is π/2. (b) Corresponding incident polarization states are represented by the red and blue points on the Poincaré sphere. (c,d) The simulated intensity profiles of the output PPBs correspond to the red points (c) and blue points (d) in (a) for three designed metasurfaces through a horizontal linear polarizer (depicted by the white double arrow).
Figure 9. Characterization of the generated PPBs. (a) The selected states of PPBs are represented by red and blue points on the HyOPS with their coordinates. The phase difference between the two circular eigenstates of the selected blue dots is π/2. (b) Corresponding incident polarization states are represented by the red and blue points on the Poincaré sphere. (c,d) The simulated intensity profiles of the output PPBs correspond to the red points (c) and blue points (d) in (a) for three designed metasurfaces through a horizontal linear polarizer (depicted by the white double arrow).
Nanomaterials 12 03010 g009
Figure 10. (a) The calculated normalized Stokes parameters correspond to point III in Figure 9a for three kinds of PPBs. (b) The calculated polarization orientations and distributions for three PPBs. Note that the orange double arrow represents the linear polarization orientation.
Figure 10. (a) The calculated normalized Stokes parameters correspond to point III in Figure 9a for three kinds of PPBs. (b) The calculated polarization orientations and distributions for three PPBs. Note that the orange double arrow represents the linear polarization orientation.
Nanomaterials 12 03010 g010
Table 1. Geometric parameters and polarization conversion efficiency for the selected unit cells.
Table 1. Geometric parameters and polarization conversion efficiency for the selected unit cells.
Unit123456789
L (μm)10939495560107107103103
W (μm)609810110510535445158
PCE (%)73.870.475.180.177.070.671.978.680.3
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, Q.; Fan, W.; Qin, C. Generation and Superposition of Perfect Vortex Beams in Terahertz Region via Single-Layer All-Dielectric Metasurface. Nanomaterials 2022, 12, 3010. https://doi.org/10.3390/nano12173010

AMA Style

Wu Q, Fan W, Qin C. Generation and Superposition of Perfect Vortex Beams in Terahertz Region via Single-Layer All-Dielectric Metasurface. Nanomaterials. 2022; 12(17):3010. https://doi.org/10.3390/nano12173010

Chicago/Turabian Style

Wu, Qi, Wenhui Fan, and Chong Qin. 2022. "Generation and Superposition of Perfect Vortex Beams in Terahertz Region via Single-Layer All-Dielectric Metasurface" Nanomaterials 12, no. 17: 3010. https://doi.org/10.3390/nano12173010

APA Style

Wu, Q., Fan, W., & Qin, C. (2022). Generation and Superposition of Perfect Vortex Beams in Terahertz Region via Single-Layer All-Dielectric Metasurface. Nanomaterials, 12(17), 3010. https://doi.org/10.3390/nano12173010

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