• Photonics Research
  • Vol. 11, Issue 10, 1723 (2023)
Zi-Wen Zhang1、2, Yu-Lu Lei1、2, Juan-Feng Zhu3, and Chao-Hai Du1、2、*
Author Affiliations
  • 1Center for Carbon-based Electronics, School of Electronics, Peking University, Beijing 100871, China
  • 2State Key Laboratory of Advanced Optical Communication Systems and Networks, School of Electronics, Peking University, Beijing 100871, China
  • 3Science, Mathematics, and Technology (SMT), Singapore University of Technology and Design, Singapore 487372, Singapore
  • show less
    DOI: 10.1364/PRJ.497533 Cite this Article Set citation alerts
    Zi-Wen Zhang, Yu-Lu Lei, Juan-Feng Zhu, Chao-Hai Du, "Plasmonic vortex beam emitter," Photonics Res. 11, 1723 (2023) Copy Citation Text show less

    Abstract

    Terahertz vortices prompt numerous advanced applications spanning classical and quantum communications, sensing, and chirality-based detection, owing to the inherent physical properties of terahertz waves and orbital angular momentum (OAM). Nonetheless, existing methodologies for generating terahertz vortices face challenges such as unalterable topological charges and intricate feed networks. To address these limitations, we propose a novel approach to generate multi-mode and tunable vortex beams based on chiral plasmons. Through eigenmode analysis, the uniform helical gratings are demonstrated to support chiral plasmons carrying OAM. By leveraging their vortex characteristics and introducing modulation into the periodic system, these chiral plasmons are alternatively diffracted into high-purity vortex radiations according to the Bragg law. To validate the theory, the vortex beam emitter is fabricated and measured in the microwave regime based on the modulated scheme. Experimental results confirm the emission of vortex beams with desirable phase distributions and radiation patterns. Our findings highlight the potential of chiral plasmons as seeds for tunable and compact vortex radiation, offering promising applications in tunable vortex sources.

    1. INTRODUCTION

    The angular momentum of photons encompasses both spin and orbital components, which are related to the polarization and wavefront distributions of electromagnetic waves. In contrast to the limited states of spin angular momentum, such as right and left rotations, orbital angular momentum (OAM) exhibits an infinite range of states, spanning from negative infinity to positive values, as demonstrated in seminal work by Allen et al. [1]. Exploiting the advantages of OAM, such as topological protection and mode orthogonality within different quantum states [2,3], optical vortex photons have found numerous applications in optical tweezers [4], communications [5], quantum entanglement [6], and nonlinear optics [7]. Extending optical vortices into the terahertz (THz) band, leveraging the broadband and unique properties of THz waves, further promotes the development of THz communications [8] and chiral spectroscopy [9]. As a result, significant efforts have been dedicated by the community toward generating high-performance THz vortex beams. However, achieving tunable vortex beams over a broad range remains an urgent need for advanced multi-mode and high-purity applications. The spiral phase plate is first proposed to apply helical phase constraint for the incident Gaussian beam, where the topological charge (TC) is structurally limited to certain values without tunability in OAM and operation band [10]. Inspired by holographic optical photonics, a computer-assisted THz hologram has been fabricated to spatially govern the in-plane phased distributions of the incident beam, whereas their tunability is still limited [11]. Through the advancement and miniaturization of metamaterials, the V-shaped slits have been tailored to function as vortex plates, precisely designed to operate at frequencies that align with the structural resonant frequencies [12]. Based on the comparable mechanism, emission of nondiffractive Bessel vortex beams at 2.13 THz has been achieved through the implementation of Fresnel spiral plates, which have been experimentally demonstrated using the Novosibirsk free-electron lasers [13,14]. Furthermore, the conversion of radially polarized beams into high-power and broadband vortex beams has been realized through the use of achromatic elements across the frequency range of 0.75 to 2 THz [15]. To summarize, these methodologies have been meticulously crafted to manifest THz vortices with a precise TC at either a singular frequency or within narrow bandwidths. However, the constrained tunability of the TC and the narrow operating bands pose limitations on their potential for advanced multi-mode applications. To surmount these obstacles, the utilization of chiral plasmonics carrying OAMs holds immense promise, offering a potential solution for achieving a tunable TC separately in multiple narrow bands.

    Spoof surface plasmons (SSPs) are artificial counterparts of optical surface plasmons that operate across the microwave to infrared bands by leveraging structural periodic systems [16]. Remarkably, these plasmons eliminate the requirement for metals with negative permittivity, presenting a novel avenue for replicating the behavior of natural plasmons [17]. Due to their dispersion relations exhibiting large longitudinal wave vectors, these mimicking modes exhibit near-field local effects and transverse surface-field enhancement. Consequently, they enable compact and intense applications in the THz bands, including integrated plasmonic sensors and waveguides [1820], subwavelength focusing [21], and far-field optical hyperlens [22]. Advancements in material topology and functional requirements have led to the evolutionary development of distinctive plasmonic modes characterized by their structural and functional properties [23]. Examples include planar SSPs [24], conformal SSPs [24], and localized SSPs [24]. Among this diverse range of SSPs, chiral ones offer inherent versatility for achieving tunable vortex radiation, as they naturally carry OAMs associated with different mode orders [25]. In efforts to generate vortex beams, a compound scheme has been devised where supplemental diffraction structures are placed around the helical grating, resulting in the emission of multi-mode vortex radiation [26]. However, this approach necessitates an additional concentric diffractive structure, which poses challenges to practical robustness. Consequently, the exploration of a scheme for diffracting localized chiral plasmons into vortex beams remains an active area of research. Promisingly, the implementation of periodic modulations holds potential in overcoming this hurdle [2729]. While vortex free-electron radiation based on chiral plasmons can generate tunable vortex beams, the constraints associated with generating symmetry-customized electron beams pose challenges in the design of electron optical systems and dynamic vacuum systems [29]. As a result, it remains a considerable distance away from practical applications and development into a high-performance vortex beam generator. Therefore, the exploration of engineering-feasible and implementable generators that enable the creation of multi-mode and tunable vortex beams remains an ongoing pursuit. Such research would undoubtedly make a substantial and valuable contribution to THz spectroscopy and communication.

    Here, a novel approach is presented for the design of a compact and tunable vortex beam emitter that leverages chiral plasmons with separate OAMs. Through eigenmode analysis, the presence of distinct TCs carried by the photons is revealed by the spiral wavefronts observed in the spatial phase distributions in real space. To achieve the desired properties, double-period depth modulation is incorporated into the proposed helical periodic systems, enabling a gentle phase-matching condition between the chiral plasmons and free space, as dictated by the fundamental principle of the Bragg law. Consequently, vortex beams with desired TCs are alternatively generated as the frequency increases. Furthermore, a gradient coupling section is designed to effectively excite these plasmons and facilitate their transformation into high-purity vortex radiations. Experimentally, a plasmonic emitter is fabricated in the microwave region to validate the operational principle in terms of radiation patterns and the phase distributions of TCs. The proposed methodology fully exploits the OAM of the chiral plasmons and Brillouin-folding phenomena, boosting the development of compact and tunable generations of THz vortex beams.

    2. RESULTS AND DISCUSSION

    A. Chiral Plasmons Carrying OAMs

    As illustrated in Fig. 1(a), the double-slots helical grating is employed here for achieving the compact and tunable vortex beam emitter since its eigenmodes are a series of chiral plasmons carrying OAMs. Therefore, the structural description and dispersion analysis are necessary for the following discussion and demonstration of the vorticity of these plasmons. The propagation and modulation periods of the helical periodic system are d and L=2d, and the slot width is a=d/2. The outer radius of the helical grating is R1=3.5d, whereas the inner radii of the two slots are R2=R3=2d for the uniform scheme. In the following, based on the multi-slots uniform helical grating, the general solution for its eigenmodes will be revealed to facilitate the understanding of the origin of the OAMs carried by chiral plasmons, where the cylindrical coordinate system is chosen for analysis. It is also supposed that all the field components are steady-state varying sinusoidally in time with monochromatic circular frequency ω, that is, time-harmonic fields. The helical grating used for the analysis contains G small slots of propagation period d within the modulation period L, thus G=L/d. Assuming the G slots have the same depth, the helical grating is thus uniform with a minimum period of d. For forward TM-polarized (electric field parallel to the z direction) modes existing in a right-handed periodic system, the E components should be Bloch states that can be labeled by the Bloch wave vector k=kzbz. Among them, kz is the propagation constant lying in the Brillouin zone, and bz=y^2π/d is the primitive reciprocal lattice vector. Thus, these eigenmodes are in the form [30]Ekz(r,φ,z)=ejkzz·ukz(r,φ,z),where (r,φ,z) are the coordinates in the cylindrical coordinate system, and ukz(r,φ,z) is a periodic function consistent with the periodicity of the helical system. According to the discrete translational symmetry along the z direction and the rotational symmetry around the z axis, the function ukz(r,φ,z) should satisfy the following periodic conditions: ukz(r,φ,z)=ukz(r,φ,z+κd),κZ,ukz(r,φ,z)=ukz(r,φ+l·2π,z),lZ.

    Schematic diagram and dispersion analysis of double-slots helical grating with or without depth modulation. Panels (a), (b), and (c) present the structural descriptions, dispersion curves, and field distributions of eigenmodes for the uniform grating. On the other hand, panels (d), (e), and (f) depict the structural descriptions, dispersion curves, and field distributions of eigenmodes for the modulated grating.

    Figure 1.Schematic diagram and dispersion analysis of double-slots helical grating with or without depth modulation. Panels (a), (b), and (c) present the structural descriptions, dispersion curves, and field distributions of eigenmodes for the uniform grating. On the other hand, panels (d), (e), and (f) depict the structural descriptions, dispersion curves, and field distributions of eigenmodes for the modulated grating.

    In satisfying the periodic conditions above, the periodic function ukz(r,φ,z) relating to φ and z can be expanded into a double Fourier series of the form ukz(r,φ,z)=R(r)l=Aφlejlφκ=Azκejκ2πdz,and the field expression Eq. (1) becomes Ekz(r,φ,z)=R(r)l=Aφlejlφκ=Azκej(kz+κ2πd)z,where R(r) is an eigenfunction of radial distance r, and Aφl and Azκ are the complex Fourier amplitudes of the electric field associated with the azimuthal lth and axial κth space harmonics, and can be evaluated by Aφl·Azκ=02π0dukz(r,φ,z)R(r)  ejlφejκ2πdzdφdz.

    The initial coordinate is assumed to be (r0, φ0, z0), where r0 is the outer radius of the helical grating. When the eigenmode described by Eq. (5) propagates an arbitrary distance Δz along the z direction, to remain invariant, its azimuthal coordinate φ should also be rotated by an angle Δφ satisfying Δφ=Δz2πd·G,where G is the number of slots in a period. This one-to-one correspondence is also requested inherently by the periodic helical structure. Hence, the periodic function at coordinates (r0, φ0, z0) and (r0, φ0+Δφ, z0+Δz) should be identically described by the following relationship: ukz(r0,φ0,z0)=ukz(r0,φ0+Δφ,z0+Δz)=ukz(r0,φ0+Δz2πd·G,z0+Δz).

    Since the electric field has the form of Fourier series, the above equation becomes R(r0)l=Aφ0lejlφ0κ=Az0κejκ2πdz0=R(r0)l=Aφ0lejlφ0κ=Az0κejκ2πdz0·ej2πd(lGκ)Δz.

    To ensure that the above equation is valid and physically achievable, the azimuthal and axial index of space harmonic should be identical, which gives {κ=lG,Aφ0l0,Az0κ0κlG,Aφ0l=0,Az0κ=0.

    Hence, the electric field in Eq. (5) must be single Fourier series with κ=l/G: Ekz(r,φ,z)=m=Fl(r)ejlφej(kz+lG2πd)z,where the Fl(r)=Aφl·Azκ·R(r) is a function that satisfies Helmholtz’s equation in cylindrical coordinates. Hence, in the region of r>r0, the function Fl(r) should be a series of modified Bessel function Kl(τlr), where τl2=k2[kz+2πl/(Gd)]2 denotes the transverse angular wavenumber and k=ω(μ0ε0)(1/2) is the wavenumber of the free space. By utilizing the method of Borgnis’s potential, the Ez field components of the slow wave in the region of r>r0 could be expressed as follows: Ez=l=τl2AlKl(τlr)ejlφej(kz+lG2πd)z,where Al is the coefficient of the Fourier series of the scalar function utilized in the method of Borgnis’s potential: BU=l=AlKl(τlr)ejlφej(kz+lG2πd)z.

    As described in Eq. (12), each eigenmode is characterized by the quantum number l, which is both the angular order and the longitudinal order of space harmonics. These eigenmodes represent space-sinusoidal waves, and their combination must satisfy the boundary condition of the helical grating. Interestingly, each mode exhibits a phase rotation factor ejlφ, which corresponds to the rotational factor in the expressions of vortex beams. Thus, the lth-order mode carries OAM with a TC of l. When the dispersion relations of these space harmonics fall below the light line, chiral SSPs with TC=l are formed. In particular, the analysis pertains to a right-handed helical system. In the case of a left-handed helical system, the OAM exhibits rotation in the opposite direction, i.e., TC=l.

    To validate the accuracy of the theoretical analysis, we employed the Eigenmode Solver in CST Studio Suite to calculate and analyze the dispersion and field distributions of the uniform double-slot helical grating, as illustrated in Figs. 1(b) and 1(c). Notably, the plasmonic modes showcased in this study exhibit positive group velocities and possess transverse even symmetry. This selection is based on the subsequent utilization of coaxial modes characterized by uniform field distributions in the transverse direction, which gradually transform into plasmonic modes. As this transformation occurs, the modes undergo equal compression within the two slots of the helical gratings. Consequently, only modes exhibiting even rotational symmetries in the transverse direction can be efficiently excited. This behavior is evident from the field distributions of plasmonic modes M0, M2, and M4. For expression Mm±, the lower and upper corners denote the absolute value of the TCs and the rotational direction of OAM. For mode M0, the first row of Fig. 1(c) presents its Ez field and phase distributions. Its electric field exhibits an almost uniform distribution along the z axis, while the phase remains nearly constant. To verify the value and direction of the OAM carried by it, the normalized spectral weight or the so-called mode purity of the TC of mode M0 is calculated by [31,32]Wm=Pml=Pl,where Pm indicates the OAM spectrum weight of the spiral harmonic basis ejmφ. To evaluate the mode purity, a circular curve, aligned with the helical grating, is employed as the sampling curve. Data points are uniformly placed along the sampling curve to capture the field properties accurately. In experimental measurements, both the amplitude and phase of the Ez field can be accurately measured, and these measurements can be conveniently converted to determine the real and imaginary parts of the field components. The mode purity can then be assessed using the discrete Fourier series provided by the following equation: Pl(ρ,z)=1Nn=1NQ(ρ,φn,z)ejlφn,where N is the number of data points and Q(ρ,φn,z) is the complex quantity consisting of real and imaginary parts of the Ez field sampled on the sampling curve from the simulation: Q(ρ,φn,z)=R(Ez)+jI(Ez).

    As illustrated in the first row of Fig. 1(c), the calculated mode purities based on Eq. (12) demonstrate that mode M0 exhibits a high-purity OAM with a TC=0, yielding a purity value of 0.9926. For mode M2, its field and phase distributions are illustrated in the second row of Fig. 1(c). When revolving around the z axis, the energy distribution remains nearly uniform, while the phase completes two full cycles from 0 to 2π. This behavior signifies its possession of OAM with a TC=2, as corroborated by its mode purity of 0.9933. Likewise, for mode M4, its phase distribution is depicted in the third row of Fig. 1(c), revealing its association with OAM characterized by a TC=4 and a remarkable purity of 0.9918. Thus, M0, M2, and M4 correspond to chiral plasmons carrying high-purity OAM with TCs of 0, 2, and 4, respectively, aligning seamlessly with the theoretical predictions provided by Eq. (12). Under double-slots helical grating configuration, the Ez field of each dispersion line ω(m,kz0) shown in Fig. 1(b) with TC=m in the region 0<(L/2π)·kz0<2 can be expressed as Ez,(m,kz0)0=τm2AmKm(τmr)ejmφejkz0z,mZ,0<L2π·kz0<2,which is a simplified form of the general solution of Eq. (12), and can be utilized to describe the vortex characteristics of the supported chiral plasmons of the helical grating. In the following, we will focus only on the dispersion lines within 0<(L/2π)·kz0<2. As depicted in Fig. 1(b), this is the smallest region within the momentum space for which the dispersion ω(m,kz0) possesses the full characteristics of the helical periodic system owing to translational symmetry. In addition, according to the periodic conditions Eqs. (2), (3), and (7), the total dispersion relation ω(l,kz) over the momentum space should also be periodic in kz and l, thereby satisfying the following condition: ω(l,kz)=ω(l+N·G,kz+N2πd)=ω(l+N·G,kz+N·G·2πL),where NZ. Since G equals 2 in the proposed M0 double-slots system, the above equation should take the following form: ω(l+2N,kz+2N·2πL)=ω(l,kz).

    Accordingly, the Ez field of an arbitrary mode M|m+2N|± within 2N<(L/2π)·kz2N<2(N+1) can be expressed by the following transformation of Eq. (17) according to Eq. (19): Ez,(m,kz2N)2N=τm2AmKm(τmr)ejlφejkz2Nz,2N<L2π·kz2N<2(N+1),where l=m+2N is defined as the integrated TC of the dispersion line for its corresponding Brillouin region.

    For the uniform grating, it can be observed from Fig. 1(b) that the dispersion curves of the supported chiral plasmons lie below the light cone and their wave vectors exceed those of free space. Consequently, these plasmons are tightly confined to the surface of the grating. To facilitate their diffraction into free space and the generation of vortex beams, a scheme employing modulation of the grating slot depth is adopted. This approach harnesses the Brillouin-zone folding phenomenon, compensating for the wave vector mismatch between these modes and the free-space modes. By exploiting band folding, the discrepancy in wave vectors is remedied, allowing for efficient diffraction of the excited chiral plasmons into vortex wave beams. Figure 1(d) presents the structural schematic of the modulation scheme, where one slot is deeper than the other, as indicated by R2=0.9R3. The field and phase distributions of modes M0, M2, and M4, computed using the CST Eigenmode Solver, are shown in Fig. 1(f). These distributions confirm that these modes are chiral plasmons carrying OAMs with TCs of 0, 2, and 4, respectively. The dispersion curves of the modulated scheme obtained through eigenmode analysis are presented in Fig. 1(e). Compared to the uniform one, these curves fold back within the boundaries of the Brillouin zone, falling into the fast-wave region (yellow shaded region). Consequently, when excited, these chiral plasmons within the fast-wave region undergo diffraction with the 1st spatial harmonic, resulting in the generation of free-space vortex beams. According to Eq. (17) and Bloch’s theorem, the diffraction wave of the mth mode M|m|+() within 0<(L/2π)·kz0<2 can be written in the following general form: Ez,(m,kz0)0=τm2AmKm(τmr)ejmφejkz0z·Z(φ,z),0<L2π·kz0<2,where Z(φ,z) is a periodic function with the same periodicity as the modulated grating. Therefore, the following three periodic conditions hold: {Z(φ,z)=Z(φ,z+n·L),nZ,Z(φ,z)=Z(φ+n·2π,z),Z(φ,z)=Z(φ+Δφ,z+Δz)=Z(φ+Δz2πL,z+Δz),where the distance Δz and the angle Δφ should satisfy Δφ=2πΔz/L. Similar to the analysis method for ukz(r,φ,z) relating to Eqs. (4) to (10), Z(φ,z) can be expressed as a Fourier series of kz: Z(φ,z)=n=Bnejnφejn2πLz,where n is the diffraction order of the space harmonic belonging to the modulated helical system, and the coefficient Bn can be calculated by Bn=02π0LZ(φ,z)ejnφejn2πdzdφdz.

    Therefore, the fundamental general form Eq. (21) of the mode within 0<(L/2π)·kz0<2 becomes Ez,(m,n,kz0)0=n=τm2AmBnKm(τmr)ejlφejkznz,n<L2π·kzn<n+2,where l indicates the integrated TC carried by the nth-order diffraction of modes M|m|± within 0<(L/2π)·kz0<2, l=m+n.

    Then, the Ez field of an arbitrary mode M|m+N|± within N<(L/2π)·kzN<N+1 can be obtained by Ez,(m,n,kzN)N=n=τm2AmBnKm(τmr)ejlφejkzn+Nz,n+N<L2π·kzn+N<n+N+2,where l indicates the integrated TC carried by the nth-order diffraction of modes M|m+N|± within N<(L/2π)·kzN<N+1, l=m+n+N,thereby enabling additive (subtractive) manipulation of the TC through diffraction. In addition to this, the nth-order diffraction of M|m+N|± also superimposes the propagation constants kzN with n2π/L, i.e., kzn+N=kz0+(n+N)2π/L. This is the key to achieving wave vector compensation between the chiral plasmons and free-space vortex beams, thus converting the modes M|m+N|± into vortex beams. Regarding Eq. (27), the compensation requires the relation 1<(L/2π)·kzn+N<+1 to be established, i.e., n+N=1. Hence, for arbitrary SSP M|m|±, the ultimate TC of the vortex beam by its diffraction is l=m1.

    Based on the above discussion, the (1+N)th-order diffraction of an arbitrary mode M|m+N|± within N<(L/2π)·kzN<N+1 yields the vortex beam with integrated TC=m1. Since the intensity of the fundamental space harmonics is the highest, only the diffraction characteristics of the SSPs M|m|± within 0<(L/2π)·kz0<2 will be discussed in the following. For the sake of meeting the compensation condition, the diffraction order should be n=1, and Eq. (25) becomes Ez,(m,1,kz0)0=n=τm2AmBnKm(τmr)ej(m1)φejkz1z,1<L2π·kz1<+1.

    It can be confirmed that the propagation constant kz0 decreases by 2π/L from Eq. (30). Therefore, as the modes within the light cone undergo 1st-order diffraction, their wave vectors remain confined within the light cone. This wavenumber compensation process enables the direct transformation of these chiral plasmons into vortex beams. Therefore, with the introduction of depth modulation in the helical grating, the modes M0, M2, and M4 carrying TCs of 0, 2, and 4, respectively, undergo diffraction and transform into vortex beams with TCs of 1, 3, and 5 according to the ultimate TC indicated by Eq. (29).

    B. Excitation of Chiral Plasmons and Generation of Vortex Beams

    To efficiently excite these chiral plasmons and validate their conversion into vortex radiation, a coupler is designed as shown in Fig. 2(a). Based on a gradient grating and a carefully designed outer conductor structure, this coupler enables the transformation of coaxial modes into plasmons carrying OAMs. The coupling section consists of an inner conductor with gradually increasing grating depth and an outer conductor with contour lines gradually moving away from the grating. The total length along the z direction of the coupling section is G=8d, corresponding to eight transmission periods. As the grating depth of the inner conductor gradually increases from 0 to h, the contour lines of the outer conductor move away from the inner conductor according to the equation y=C1eαz+C2(z1<z<z2), gradually converting coaxial modes into plasmonic modes on the helical grating, where C1=(y2y1)/(eαz2eαz1), C2=(y1eαz2y2eαz1)/(eαz2eαz1), and α=0.04. To illustrate the working principle of the coupling section more clearly, the dispersion of mode M0 as a function of grating depth, with a spacing of 0.17h, is shown in Fig. 2(b). When the grating depth is 0, the dispersion of the mode M0 coincides with that of the coaxial mode, perfectly matching the dispersion curve of the light. As the depth of the gradient grating gradually increases, the dispersion curve of the mode M0 gradually splits from the light line, with an increasing wave vector in the longitudinal direction, exhibiting the characteristics of near-field localization of slow-wave modes. Therefore, based on the depth-gradient grating, the mismatched wave vectors can be gradually compensated, allowing for an effective transformation of coaxial modes into plasmonic modes. To validate the performance of the designed coupler, the field distribution of the system at the normalized frequency fnorm=0.24 is also presented in Fig. 2(c). It can be observed that the coaxial mode excited at Port 1 gradually converts into the mode M0 on the grating with the assistance of the coupler. As the plasmonic waves crawl over the grating surface and reach another coupler, the mode M0 is efficiently converted back into the coaxial mode based on the reciprocity principle and transmitted to Port 2. Apart from mode M0, similar analysis can be applied to modes M2 and M4, which can also be effectively excited from coaxial modes with the assistance of the coupler.

    Excitation of the chiral plasmons. (a) Coupler utilized for the transformation of coaxial modes into plasmonic modes. (b) Dispersion curves of mode M0 as the slot depth increases. (c) Field distribution of plasmonic transmission.

    Figure 2.Excitation of the chiral plasmons. (a) Coupler utilized for the transformation of coaxial modes into plasmonic modes. (b) Dispersion curves of mode M0 as the slot depth increases. (c) Field distribution of plasmonic transmission.

    To further validate the excitation effectiveness, a sampling line is introduced along the propagation direction on the grating surface to examine their spatial distribution in the wavenumber domain, as shown by the red dashed line in Fig. 3(a). This sampling line allowed us to capture the spatial distribution of the excited modes at different frequencies. Subsequently, a one-dimensional spatial Fourier transform is performed on the field distribution at each frequency, resulting in the wavenumber distribution of the excited modes, represented by the red circles in Fig. 3(b). Upon analysis, it is observed that as the frequency increases, the wavenumber distribution of the excited modes closely matches the dispersion curve of the grating. This alignment indicated an effective transformation of the coaxial modes into modes M0, M2, and M4, which can be further supported by the excellent transmission characteristics from Port 1 to Port 2, as depicted in Fig. 3(c). However, it is important to note the presence of a photonic bandgap (PBG) between the dispersion curves of modes M2 and M4. The frequency range of this bandgap, marked by the red shaded region in Fig. 3(b), indicated that the grating does not support the propagation of modes within this range. As a result, the transmission performance experienced a decline in this frequency range, as demonstrated in Fig. 3(c). Consequently, when the coaxial mode is output from the coupler, it cannot be converted into plasmons propagating along the grating surface within the PBG. Instead, it directly generates vortex radiation propagating toward free space, carrying OAM akin to its corresponding mode. Figure 3(b) reveals that the wavenumber distribution of the field detected within the PBG is almost an extension of the dispersion curve of mode M4 in the fast-wave region. Hence, this mode can be considered a variant of mode M4 within the outer conductor and helical grating system, exhibiting similar transverse field distribution and carrying the same OAM as mode 4. Consequently, the vortex radiation generated within the PBG carries a TC=4.

    Verification of the excitation of the chiral plasmons. Panels (a), (b), and (c) show the schematic diagrams, sampling results, and transmission characteristics of the uniform grating, while panels (d), (e), and (f) depict the schematic diagrams, sampling results, and transmission characteristics of the modulated grating.

    Figure 3.Verification of the excitation of the chiral plasmons. Panels (a), (b), and (c) show the schematic diagrams, sampling results, and transmission characteristics of the uniform grating, while panels (d), (e), and (f) depict the schematic diagrams, sampling results, and transmission characteristics of the modulated grating.

    Then, the depth modulation scheme is introduced to convert the excited plasmons carrying OAMs into vortex radiation, as illustrated in Fig. 3(d). Building upon the theoretical analysis presented earlier, the effectiveness of this excitation method is further validated by examining the wavenumber distribution of the field excited on the modulated grating surface at different frequencies, as depicted by the red circles in Fig. 3(e). Remarkably, the distribution of the excited field in momentum space is closely aligned with the dispersion curves of the eigenmodes supported by the modulated grating, thereby substantiating the efficacy of this excitation approach. In comparison to the uniform grating, the presence of Brillouin zone folding led to a partial overlap between the dispersion curves of modes M2 and M4 with the fast-wave region, as evidenced by the orange and blue shaded regions in Fig. 3(e). These regions are identified as bands B1 and B3, respectively. Consequently, the modes M2 and M4 within the bands B1 and B3 undergo 1st-order diffraction, resulting in their conversion into vortex radiation propagating toward free space. Specifically, the vortex beam from mode M2 acquired a TC=1, while the other one exhibited a TC=3. These findings highlight the transformation of coaxial modes into chiral plasmons through the interaction with the modulated grating. In addition to these discovered radiation bands, two additional bands labeled as bands B2 and B4 emerged due to the emergence of the PBG, as indicated by the red shaded regions in Fig. 3(e). Band B4, formed within the frequency interval created by the PBG between modes M2 and M4, exhibited a transverse field distribution similar to that of mode M4. However, owing to its compatibility with free-space wave vectors, the waves within the band B4 directly transformed into vortex radiation carrying an OAM with TC=4. Similarly, band B2 originated from the PBG within the slow-wave region resulting from the splitting band of the mode M2. The transverse field distribution of the converted coaxial modes within this frequency range resembled mode M2, enabling their direct transformation into vortex radiation with the OAM TC of 2, thanks to their wave vector matching with free space. Consequently, the depth-modulated helical system facilitated the simultaneous generation of vortex beams carrying TCs of 1, 2, 3, and 4 within the frequency intervals of radiation bands B1,B2, B3, and B4, respectively. The radiation characteristics of these bands are further corroborated through the analysis of the transmission properties depicted in Fig. 3(f). As anticipated, the presence of 1st-order diffraction and the PBG resulted in a reduction in the transmission performance within the frequency intervals of radiation bands B1,B2, B3, and B4, aligning with the simulated transmission characteristics and confirming the successful conversion and transmission of chiral plasmons in the system.

    To further confirm the validity of the aforementioned analysis, we perform simulations at normalized frequencies 0.2587, 0.3227, 0.3941, and 0.4203 within the bands B1,B2, B3, and B4, respectively. These calculations focus on evaluating the distribution of the radiation field, phase distribution, far-field directivity, and mode purity associated with the OAMs carried by the beams. The results, illustrated in Fig. 4, unequivocally demonstrate the generation of vortex beams with distinct TCs of 1, 2, 3, and 4, respectively. The vortex radiation within bands B1 and B3 results from the 1st-order spiral diffraction of modes M0 and M2, transforming into vortex beams propagating in free space. According to Eq. (1), their expected TCs are 1 and 3, which harmoniously align with the results obtained from the simulation analysis in Figs. 4(a) and 4(c). Additionally, their mode purities are determined to be 0.934 and 0.768, respectively. In contrast, the vortex radiations within bands B2 and B4 are directly generated from the fast-wave modes with transverse distributions resembling modes M2 and M4. Hence, their TCs are expected to be 2 and 4, respectively, in agreement with the results of the simulation analysis in Figs. 4(b) and 4(d). Their corresponding mode purities are evaluated as 0.789 and 0.876, respectively. Thus, the depth-modulated scheme enables the simultaneous generation of vortex beams with TCs ranging from 1 to 4 within the four radiation bands. This achievement not only validates the effectiveness of depth modulation for controlled vortex plasmon emission but also underscores its potential for various applications.

    Generation of the vortex beams. (a), (b), (c), and (d) are radiation field distributions, far-field directivities, and mode purities of generated vortex beams in bands B1,B2, B3, and B4.

    Figure 4.Generation of the vortex beams. (a), (b), (c), and (d) are radiation field distributions, far-field directivities, and mode purities of generated vortex beams in bands B1,B2, B3, and B4.

    3. EXPERIMENTAL VERIFICATION

    To experimentally validate the working principle, we fabricate the plasmonic vortex beam emitter in the microwave band, as shown in Fig. 5(d). The helical grating is first printed out with the whole frame by aluminum alloy AlSi10Mg through a commercial three-dimensional printing technique and then coated with copper by the electroplating process. The structural parameters suitable for the microwave band are as follows: d=3.60  mm, L=7.20  mm, a=1.80  mm, R1=12.60  mm, R2=7.20  mm, and R3=6.48  mm. It is noteworthy that the actual manufactured helical grating has increased in size after being coated with copper, and the measured outer radius R1 is 12.62 mm instead of 12.60 mm. Thus, the machining error is approximately 0.1587%, which is within an acceptable range. The overall experimental setups for the vortex beam antenna are presented in Fig. 5(d), comprising an Agilent N5245A PNA-X network analyzer (10 MHz to 50 GHz), a high-gain horn antenna for receiving the vortex radiation signal, a turntable for signal reception from different directions, a computer for programming the vector network analyzer (VNA), and the fabricated vortex emitter. During the tests, the vortex beam emitter is connected to the VNA through an N-type connector and powered through the SubMiniature version A (SMA) port. Port 1 and Port 2 of the VNA are connected to the couplers of the plasmonic emitter, while Port 3 is connected to the horn antenna for measuring far-field radiation performance. Figure 5(c) presents the transmission characteristics at both ends of the vortex beam emitter obtained from experimental testing and simulations accounting for losses, marked by blue and red lines, respectively. The simulation results exhibit good agreement with the experimental measurements, showing a consistent decrease in transmission characteristics across bands B1,B2, B3, and B4, as predicted by the theoretical analysis. Waves that do not transmit through the grating are transformed into free-space vortex beams, which are detected by the horn antenna in the far field, as confirmed by the measured dispersion and beam-scanning characteristics shown in Figs. 5(a) and 5(b).

    Experimental setups and measured results. (a) Measured dispersion of the radiation waves from the vortex beam emitter. (b) Measured and theoretical beam scanning characteristics. (c) Measured and simulated transmission performances between two couplers. (d) Experimental setups. (e), (f), (g), and (h) are phase distributions of generated vortex beams in bands B1,B2, B3, and B4.

    Figure 5.Experimental setups and measured results. (a) Measured dispersion of the radiation waves from the vortex beam emitter. (b) Measured and theoretical beam scanning characteristics. (c) Measured and simulated transmission performances between two couplers. (d) Experimental setups. (e), (f), (g), and (h) are phase distributions of generated vortex beams in bands B1,B2, B3, and B4.

    Within bands B1 and B3, the vortex radiation results from the 1st-order diffraction of modes M0 and M2, giving rise to wide-angle beam scanning characteristics. As the frequency decreases, the radiation direction increases from 30° to 150°. The phase distributions at frequencies 0.2618 and 0.3818 within these bands are shown in Figs. 5(e) and 5(g), respectively. As the detection horn rotates a full circle, the phase of the longitudinal electric field at these frequencies decreases by 1 and 3 times, respectively, confirming their TCs of 1 and 3. In addition to these two bands, two radiation spectra are observed within the 0° to 30° radiation range, as depicted in Fig. 5(b). These spectra correspond to the PBG bands B2 and B4, whose waves are directly radiated out of the coupler without beam scanning capability. The phase distributions at frequencies 0.3314 and 0.4131 within these bands are shown in Figs. 5(f) and 5(h), respectively, decreasing by 2 and 4 times from 0 to 2π as the horn rotates, confirming their TCs of 2 and 4. Thus, through the spectral and phase distributions of the radiation waves, we have demonstrated that the designed vortex beam emitter can simultaneously generate vortex radiation with TCs ranging from 1 to 4.

    4. CONCLUSION

    In conclusion, we present a theoretical proposal and experimental demonstration of a chiral plasmons-enabled vortex beam emitter, capable of simultaneously generating multi-mode and highly pure vortex beams separate in multiple narrow bands. By incorporating depth modulation into the helical periodic system, the eigenmodes carrying OAMs undergo alternating transformations into vortex radiations, following the principles of the Bragg law. Through experimental observations of the measured radiation patterns and phase distributions, our designed plasmonic vortex beam emitter convincingly exhibits the generation of vortex beams with TCs of 1, 2, 3, and 4. This proposed scheme offers a promising avenue for the generation of multi-mode, tunable, and high-purity vortex beams in the THz frequency range, paving the way for advanced OAM-based communication and spectroscopy.

    Acknowledgment

    Acknowledgment. The authors acknowledge discussions with Fan-Hong Li.

    References

    [1] L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, J. P. Woerdman. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Phys. Rev. A, 45, 8185-8189(1992).

    [2] G. Indebetouw. Optical vortices and their propagation. J. Mod. Opt., 40, 73-87(1993).

    [3] Y. Shen, X. Wang, Z. Xie, C. Min, X. Fu, Q. Liu, M. Gong, X. Yuan. Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities. Light Sci. Appl., 8, 90(2019).

    [4] M. Padgett, R. Bowman. Tweezers with a twist. Nat. Photonics, 5, 343-348(2011).

    [5] A. E. Willner, H. Huang, Y. Yan, Y. Ren, N. Ahmed, G. Xie, C. Bao, L. Li, Y. Cao, Z. Zhao, J. Wang, M. P. J. Lavery, M. Tur, S. Ramachandran, A. F. Molisch, N. Ashrafi, S. Ashrafi. Optical communications using orbital angular momentum beams. Adv. Opt. Photonics, 7, 66-106(2015).

    [6] A. Sit, F. Bouchard, R. Fickler, J. Gagnon-Bischoff, H. Larocque, K. Heshami, D. Elser, C. Peuntinger, K. Günthner, B. Heim, C. Marquardt, G. Leuchs, R. W. Boyd, E. Karimi. High-dimensional intracity quantum cryptography with structured photons. Optica, 4, 1006-1010(2017).

    [7] J. Vieira, R. M. G. M. Trines, E. P. Alves, R. A. Fonseca, J. T. Mendonça, R. Bingham, P. Norreys, L. O. Silva. Amplification and generation of ultra-intense twisted laser pulses via stimulated Raman scattering. Nat. Commun., 7, 10371(2016).

    [8] H. Zhou, X. Su, A. Minoofar, R. Zhang, K. Zou, H. Song, K. Pang, H. Song, N. Hu, Z. Zhao, A. Almaiman, S. Zach, M. Tur, A. F. Molisch, H. Sasaki, D. Lee, A. E. Willner. Utilizing multiplexing of structured THz beams carrying orbital-angular-momentum for high-capacity communications. Opt. Express, 30, 25418-25432(2022).

    [9] A. A. Sirenko, P. Marsik, C. Bernhard, T. N. Stanislavchuk, V. Kiryukhin, S. W. Cheong. Terahertz vortex beam as a spectroscopic probe of magnetic excitations. Phys. Rev. Lett., 122, 237401(2019).

    [10] K. Miyamoto, K. Suizu, T. Akiba, T. Omatsu. Direct observation of the topological charge of a terahertz vortex beam generated by a Tsurupica spiral phase plate. Appl. Phys. Lett., 104, 261104(2014).

    [11] Z. Xie, X. Wang, J. Ye, S. Feng, W. Sun, T. Akalin, Y. Zhang. Spatial terahertz modulator. Sci. Rep., 3, 3347(2013).

    [12] J. He, X. Wang, D. Hu, J. Ye, S. Feng, Q. Kan, Y. Zhang. Generation and evolution of the terahertz vortex beam. Opt. Express, 21, 20230-20239(2013).

    [13] B. A. Knyazev, Y. Y. Choporova, M. S. Mitkov, V. S. Pavelyev, B. O. Volodkin. Generation of terahertz surface plasmon polaritons using nondiffractive Bessel beams with orbital angular momentum. Phys. Rev. Lett., 115, 163901(2015).

    [14] B. Volodkin, Y. Choporova, B. Knyazev, G. Kulipanov, V. Pavelyev, V. Soifer, N. Vinokurov. Fabrication and characterization of diffractive phase plates for forming high-power terahertz vortex beams using free electron laser radiation. Opt. Quantum Electron., 48, 223(2016).

    [15] R. Imai, N. Kanda, T. Higuchi, K. Konishi, M. Kuwata-Gonokami. Generation of broadband terahertz vortex beams. Opt. Lett., 39, 3714-3717(2014).

    [16] J. B. Pendry, L. Martín-Moreno, F. J. Garcia-Vidal. Mimicking surface plasmons with structured surfaces. Science, 305, 847-848(2004).

    [17] Z. Gao, L. Wu, F. Gao, Y. Luo, B. Zhang. Spoof plasmonics: from metamaterial concept to topological description. Adv. Mater., 30, 1706683(2018).

    [18] L. Tong, H. Wei, S. Zhang, H. Xu. Recent advances in plasmonic sensors. Sensors, 14, 7959-7973(2014).

    [19] S. Roh, T. Chung, B. Lee. Overview of the characteristics of micro- and nano-structured surface plasmon resonance sensors. Sensors, 11, 1565-1588(2011).

    [20] L. Liu, Z. Han, S. He. Novel surface plasmon waveguide for high integration. Opt. Express, 13, 6645-6650(2005).

    [21] L. Yin, V. K. Vlasko-Vlasov, J. Pearson, J. M. Hiller, J. Hua, U. Welp, D. E. Brown, C. W. Kimball. Subwavelength focusing and guiding of surface plasmons. Nano Lett., 5, 1399-1402(2005).

    [22] Z. Liu, H. Lee, Y. Xiong, C. Sun, X. Zhang. Far-field optical hyperlens magnifying sub-diffraction-limited objects. Science, 315, 1686(2007).

    [23] F. J. Garcia-Vidal, A. I. Fernández-Domínguez, L. Martin-Moreno, H. C. Zhang, W. Tang, R. Peng, T. J. Cui. Spoof surface plasmon photonics. Rev. Mod. Phys., 94, 025004(2022).

    [24] C. R. Williams, S. R. Andrews, S. A. Maier, A. I. Fernández-Domínguez, L. Martín-Moreno, F. J. García-Vidal. Highly confined guiding of terahertz surface plasmon polaritons on structured metal surfaces. Nat. Photonics, 2, 175-179(2008).

    [25] F. Rüting, A. I. Fernández-Domínguez, L. Martín-Moreno, F. J. García-Vidal. Subwavelength chiral surface plasmons that carry tuneable orbital angular momentum. Phys. Rev. B, 86, 075437(2012).

    [26] J.-F. Zhu, C.-H. Du, Z.-W. Zhang, F.-H. Li. Generating a multi-mode vortex beam based on spoof surface plasmon polaritons. Opt. Lett., 47, 4459-4462(2022).

    [27] A. Bera, R. K. Barik, M. Sattorov, O. Kwon, S.-H. Min, I.-K. Baek, S. Kim, J.-K. So, G.-S. Park. Surface-coupling of Cerenkov radiation from a modified metallic metamaterial slab via Brillouin-band folding. Opt. Express, 22, 3039-3044(2014).

    [28] J. F. Zhu, C. H. Du, F. H. Li, L. Y. Bao, P. K. Liu. Free-electron-driven multi-frequency terahertz radiation on a super-grating structure. IEEE Access, 7, 181184(2019).

    [29] Z.-W. Zhang, J.-F. Zhu, C.-H. Du, F. Gao, F.-Y. Han, P.-K. Liu. Chiral plasmons enable coherent vortex Smith–Purcell radiation. Laser Photonics Rev., 17, 2200420(2023).

    [30] K. Zhang, D. Li, K. Chang, K. Zhang, D. Li. Electromagnetic Theory for Microwaves and Optoelectronics(1998).

    [31] H. I. Sztul, R. R. Alfano. The Poynting vector and angular momentum of Airy beams. Opt. Express, 16, 9411-9416(2008).

    [32] R. Zambrini, S. M. Barnett. Quasi-intrinsic angular momentum and the measurement of its spectrum. Phys. Rev. Lett., 96, 113901(2006).

    Zi-Wen Zhang, Yu-Lu Lei, Juan-Feng Zhu, Chao-Hai Du, "Plasmonic vortex beam emitter," Photonics Res. 11, 1723 (2023)
    Download Citation