
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39261549
71910
10.1038/s41598-024-71910-5
Article
Polarization insensitive non-interleaved frequency multiplexed dual-band Terahertz coding metasurface for independent control of reflected waves
Iqbal Shahid 1
Noor Ahsan 2
Ullah Naeem 3
Saifullah Yasir 1
Ahmed Shehzad 3
Nisar Muhammad Shemyal 4
Wong Sai-Wai wongsaiwai@ieee.org

1
1 grid.263488.3 0000 0001 0472 9649 State Key Laboratory of Radio Frequency Heterogenous Integration, Guangdong Engineering Research Centre of Base Station Antennas, Shenzhen Key Laboratory of Antennas and Propagation, College of Electronics and Information Engineering, Shenzhen University, Shenzhen, 518060 China
2 https://ror.org/049tv2d57 grid.263817.9 0000 0004 1773 1790 Electrical and Electronics Engineering, Southern University of Science and Technology, Shenzhen, 518055 China
3 https://ror.org/01vy4gh70 grid.263488.3 0000 0001 0472 9649 College of Physics and Optoelectronics engineering, Shenzhen University, Shenzhen, 518055 China
4 https://ror.org/00ay9v204 grid.267139.8 0000 0000 9188 055X Sino-British College, University of Shanghai for Science and Technology, Shanghai, 200093 China
11 9 2024
11 9 2024
2024
14 2119926 2 2024
2 9 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Independent control of electromagnetic (EM) waves by metasurfaces for multiple tasks are highly desired and is the recent hot topic of research. In this work we contribute a polarization insensitive frequency multiplexed 2-bit coding metasurface to control the Terahertz (THz) waves in the two operating bands independently. In this regard, as a first step a cascaded meta-atom composed of square rings and/or square metallic patches separated by two polyimide substrates is designed and optimized that provides sixteen independent distinct discrete phases in the reflection geometry. These meta-atoms are then distributed with distinct coding sequences in the two-dimensional spatial plane to realize various bi-functional metasurfaces. As a proof of the concept various full structures are designed and simulated to realize a series of bi-functionalities including anomalous reflection/beam shaping, beam shaping/anomalous reflection, beam deflection/Orbital angular momentum (OAM) beam generation with distinct modes and propagating wave to surface wave (PW–SW) conversion/PW beam manipulation in the lower and higher THz bands, respectively. All the simulation results are in excellent agreement with their theoretical equivalents. We envision that the proposed meta-designs have potential applications for the multi-spectral control of EM waves in THz band. The idea can be further extended to design frequency dependent tri-functional and multi-functional THz meta-devices.

Keywords

Digital coding metasurface
Frequency multiplexed
Dual-band
Polarization insensitive
Propagating waves controlling
Multi-mode OAM beams
Surface waves generation and controlling
Subject terms

Optics and photonics
Optical materials and structures
Metamaterials
http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 62171289 6225041036 Ullah Naeem Shenzhen Science and Technology ProgramsJCYJ20190808145411289 issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Metasurfaces (MTSs)1–3, allow for the exotic control and manipulation of electromagnetic (EM) waves, which have led to the demonstration of phenomena like flat lenses3, highly efficient and switchable surface plasmon couplers4–6, high-resolution three-dimensional (3D) MTS holograms7–9, and different functional interfaces10. Technologically, MTSs have bridged the gap between the engineered structures and high technological industry, showing great potential to realize novel operational devices for high resolution imaging, wireless power transfer and modern communication11–14. The working principle of these MTS based devices is to modulate the phase, amplitude and polarization of an impinging wave in a continuous manner1–3,11. From the perspective of circuit theory these designs can be considered as analog circuits15,16. On the contrary, Cui et al.17, introduced a new approach to these engineered structures by designing and defining a new class of MTSs called the digital/coding and programmable MTSs. In this class of MTSs, the degrees of freedom of EM waves (phase, amplitude, and polarization) are mapped upon binary numbers 0, 1, 00, 01 ..., and the meta devices are studied from the perspective of digital circuits, for potential applications in the realm of information sciences17–23.

In the past decade, variety of digital/coding and programmable MTSs, in distinct frequency regimes from Microwaves (MWs) to Terahertz (THz), have been proposed17–22,24–28 to demonstrate various optical functionalities. Among the distinct coding MTSs designs17–22,24–28, a digital/coding MTS design of particular interest is the re-configurable/programmable MTSs17,19,20,22,24,29,30 realized in the microwave regimes. In these designs, the functionalities of the digital/coding MTSs are controlled and programmed in real-time by embedding a lumped element such as positive-intrinsic-negative (PIN) diode or varactor diode at each unit cell. Indeed, these re-configurable coding MTSs designs have demonstrated many mature algorithms in the information and imaging theory17,18,29–31. Another interesting class of digital/coding MTS is the bi/multi-functional MTS. Such designs, realized in both the MWs and THz regimes of EM spectrum, demonstrate independent and multiple functionalities from a single MTS design26,32–34. Operationally, independent bi/multiple functions in the digital/coding MTS designs are encoded on either the polarization of the incident waves34, or frequency of operation32,33, or both25,35. Notably, digital/coding MTS bi-functional designs are advantageous over the analog MTSs working on the principle of continuous modulation phase/amplitude counterpart36,37, with regard to their relatively easier design processes and their ability to perform digital algorithms18–21,38.

THz radiations, which lie between the MWs and infrared regimes of the EM wave spectrum, have the features of strong penetration power, low photon energy, low attenuation, and wider bandwidth B.W. With these advantages, control and manipulation of the THz radiations possess tremendous potentials in various applications including high speed communication, non-invasive detection and imaging39–41 are to mention a few. The conventional THz devices, in this regard are mostly bulky, which hinders their integration to chip scales42,43, consequently, limiting their practical applications. MTSs, however, provided a readily available plateform39 to manipulate the THz waves for various functions. Likewise, distinct digital/coding MTS designs, for THz manipulation has also been proposed21,44–48. At this point, it is worth mentioning that despite the challenges in the availability of tuning elements and mechanisms in this part of EM wave spectrum, there have been attempts to realize THz re-configurable and programmable MTSs45–48, by utilizing some novel tuneable materials such as Vanadium dioxide (VO2), graphene and liquid crystals, showing the potential of programmable meta-devices at higher frequencies. Similarly, the THz digital/coding meta-devices, demonstrating multiple tasks from the same shared aperture shall enable highly miniaturized and integrated devices21,34,49. In the past few years, some studies have been published to demonstrate bi-functional THz digital/coding MTSs based on either the polarization of the incident waves (polarization multiplexed)34,50 or the frequencies of operation (frequency multiplexed MTSs)33. Meanwhile, these designs have significant limitations, especially in the case of frequency dependent coding MTSs. For example, the frequency multiplexed structures previously reported can operate for a single polarization, have limited phase resolution, narrow band of operation and angularly unstable. Moreover, the dual-band coding MTSs33,51 have been reported so far have limited phase resolution which consequently limits their number of possible functionalities.

In a bid to cover this gap, this work demonstrates a complete 2-bit dual-band polarization insensitive coding MTS for the THz regime. The proposed MTS design with its 2-bit coding resolution in each operating band provides more flexibility in manipulating the reflected waves independently in different frequency bands , much worthier from the previously published frequency dependent THz works, as those were performing limited tasks. Functionalities demonstrated here include independent deflection of the normally incident waves to different anomalously reflected angles and in different directions, Orbital angular momentum (OAM) beam generation with distinct topological charges and deflected to distinct directions, anomalous reflection and beam shaping, conversion and controlling of normally incident propagating waves to surface waves (PW–SW) in the two operating bands independently. Moreover, based on the spatial convolution theorem a compensation scheme has been used with which, regardless of the frequency band and the direction in which the anomalously reflected wave is to be deflected, the angle can be engineered desirably with negligible distortions, mitigating the strong wavelength dependence of such kind of multi-band devices. Furthermore, besides the demonstrated functionalities, all other frequency dependent bi-functionalities can also be achieved by the proper selection and distribution of the designed coding meta-atoms. Some of the dual functionalities of the proposed work are shown in the schematic illustration in Fig. 1.Fig. 1 Schematic illustration. The schematic illustrates different bi-functionalities of the proposed design. Each bi-functionality is achieved with a specific coding sequence ‘CSn’. MTSs with ‘CS1’, ‘CS2’ and ‘CS3’ illustrates beam shaping/ beam deflection, OAM beam generation with l1=m=+1/ OAM beam generation with l2=m=+2 (deflected to different planes), and PW–SW conversion/PW beam manipulation in lower and higher bands (fl/fh), respectively.

Design and simulation

Dual-band coding meta-atom design and simulation

To achieve an independent 2-bit phase response in two different frequency bands, a total of sixteen unit cells are required. Here, a three-metallic-layer structure shown in Fig. 2a is used as the coding meta-atom for a 2-bit dual-band coding MTS. The lattice constant of the proposed meta-atom is p=70um which in terms of operating wavelength (λo) of the incident waves in free space is p≈λo/7 at ‘fl’ while at ‘fh’ p≈λo/3.5. All the metallic layers used are of lossy Aluminum, where a fully metallic bottom layer acting as a ground is separated by a polyimide spacer of thickness d1=20um from the inner layer. The inner, and outer metallic layers are separated by the same polyimide with thickness d2=10um, and dielectric properties of the polyimide are ε=3.0, and loss tangent δ=0.03. A square loop-patch structure with parameters wl(i,o), wp(i,o), w(i, o), as shown in Fig. 2b and square patch with parameters wpt(i,o), as shown in Fig. 2c, have been used to load inner and outer metallic layers. Note that the subscripts ‘i and ‘o, here suggests that a coding meta-atom, may either have both layers (inner and outer) loaded with a square loop-patch or a simple square patch, or it can be a combination of both. For design and simulation of the coding meta-atom, the commercially available software computer simulation technology (CST) Microwave studio is used. The unit cell is simulated using the unit cell boundary conditions, in the x–y plane with the same period p along the x and y-axis and Floquet ports are set at a distance of 400um along the z-axis as the excitation source. The surface of the structure is selected as the reference plane, and above that region upto 400um the phases and amplitudes of reflection are extracted. The simulation results obtained are plotted in Fig. 3, where Fig. 3a and Fig. 3b show the amplitude and phase states in the lower band (0.53–0.73 THz), having a B.W ≈0.2, while Fig. 3c and Fig. 3d show the corresponding phase and amplitude spectra in higher band (1.25–1.45 THz) with a B.W≈ 0.25. Note that the amplitudes of coding meta-atoms in each working band remain in the rage of 0.7–1 in linear scale, however at central frequencies for each row this difference reduces to ≈0.15, consequently not causing any significant changes in the overall results of the full structures. Similarly, the phase differences between the consecutive digital phase states also remain approximately to 90∘±10∘ which fulfills the definition of binary digital phase states of a 2-bit coding meta-atom. Moreover, the separation between the two operating bands is ≈0.5 THz that is large enough to avoid any cross-talk between the two functions obtained from the same design. Figure 2d provides locations of the 16 optimized coding meta-atoms in the reflection phase diagram, with horizontal and vertical axes representing the reflection phases under the lower and higher frequency, respectively. It is apparent that phase response at a lower frequency in each row remain the same and are independent of the coding states in the higher frequency band (the number before the slash). The similar kind of trend appears when seen from each column, corresponding to independent digital states at higher frequency band (the number after the slash), allowing us to encode dual functionalities using a single MTS.Fig. 2 Unit cell design. (a) Perspective transparent view of the meta-atom. (b, c) Top layer and middle layer view with square loop and patch. (d) Dual-band digital phase states.

Fig. 3 Phase and amplitude spectra. (a, b) Reflection amplitude and phase at the lower band. (c, d) Reflection amplitude and phase at the higher band.

It is also very important to mention that the phase and amplitude spectra obtained here represent the overall response of the meta atom. Any changes in the dimensions of metallic patches or the substrates sizes can potentially alters this over all response. To further investigate a tolerance analysis on the meta-atom performance is carried out (see Supplementary information, Section S1) to study how the variations in the meta-atom’s parameters affect its behavior.

Full structures design and simulation

After realizing the 2-bit dual-band coding meta-atoms, the next stage is to distribute them with different coding sequences in the x–y spatial plane for designing full structures. In the following text different bi-functional structures are discussed where each structure is encoded with a specific coding sequence to perform distinct tasks independently. In the design process of full structures the idea of super unit cell is adopted, which facilitates to minimize the unwanted inter element coupling between the coding meta-atoms, and to flexibly manipulate the scattered waves. For the simulation purposes of the full structure CST MWS is used but with time domain solver and open add space boundary condition in the ± z direction. Linearly polarized plane waves are impinged to excite the structures. Furthermore, other necessary simulation conditions and requirements are explained for each design in the following text.

Results and discussion

Bi-functional full structures for beam shaping/steering and flexible beam shaping

Beam forming/shaping that includes beam splitting and quad beam generation, can be realized by employing coding MTS with a periodic ‘00/10’ coding sequence, chess board like coding pattern, respectively. Although, it has already been established that beam shaping can be achieved using coding MTSs with 1-bit phase resolution17,21,25, however, joint independent dual-band beam shaping, and beam shaping at one frequency along other tasks at another frequency band have rarely been reported. In the proposed case since independent phase modulation with 2-bit resolution can be performed thus it is feasible to achieve flexible beam shaping in the two operating bands and also beam shaping in one band along another useful functionality such as beam deflection in the other band. For this purpose, firstly we selected the coding meta-atoms from the matrix M2-bit as shown in Eq. (1), which have digital phase states with 1-bit resolution in the lower and 2-bit resolution in the higher band, i.e., ‘00/00, 01/00, 10/10, 11/10’.1 M2-bit=00/0000/0100/1000/1101/0001/0101/1001/1110/0010/0110/1010/1111/0011/0111/1011/11

A full structure labeled as S1 is designed with period of eight coding meta-atoms and overall size of x×y=(70×72)×(70×72)um2. Full wave simulations are performed for both lower and higher bands under the normal illumination of x-polarized incident waves to measure the far field patterns. The simulation results are shown in Fig. 4a–d, where Fig. 4a and c shows the 3D far field pattern, while Fig. 4b and d show their 2D equivalents, respectively. The results clearly manifest that a splitted and a deflected beam is achieved at the lower and higher band, respectively. The angle at which each beam is deflected have the numerical values which are almost same as predicted by the Generalized Snell’s law1 described in Eq. (2).2 θd=arcsinλΓ

In addition, we designed a coding MTS of the same size labeled as ‘S2’ which can deflect the reflected beam in the lower band and split the reflected beam in to two equal parts in higher band, i.e., opposite function as that of the first structure ‘S1’. This task is accomplished by selecting the coding meta-atoms which have 2-bit phase resolution in the lower and 1-bit in the higher band. The coding sequence that can perform this function is, ‘00/00, 01/00, 10/10, 11/10’. The results obtained for this case are depicted in Fig. 4e–h, where Fig. 4e and Fig. 4g shows the far field patterns while Fig. 4f and Fig. 4h shows their corresponding 2D equivalents. The results clearly shows that beam deflection and beam splitting in the lower and higher band, respectively, has been achieved efficiently. Furthermore, a MTS that can shape the beam in two distinct patterns in the two operating bands is designed. In this case the full structure labeled as ‘S3’ can shape the reflected beam into a quad beam pattern in the lower band and a splitting pattern in the higher band. This task can be achieved by selecting the coding meta-atoms which have 1-bit phase resolution in both bands but a chess board like distribution in the lower band and a periodic distribution in the higher band, i.e., ‘00/00, 10/10, 10/00, 00/10’. The results obtained in this case are depicted in Fig. 4i–l, where Fig. 4i and Fig. 4k shows the 3D far field pattern while Fig. 4j and Fig. 4l shows the corresponding 2D patterns. Here the angle at which each beam in the quad beam pattern is deflected can be calculated using the Eqs. (3) and (4)19,21.3 ϕ=±π4andϕ=±3π4

and4 θ=arcsinλ·2Γ2i.e.,θd=arcsin2λΓ

The results shown in Fig. 4i–j have very close agreement with the theoretically predicted results of Eqs. (3) and (4).Fig. 4 Beam shaping and deflection. (a–d) 3D and 2D far field pattern of beam splitting and beam deflection in the lower and higher bands,respectively. (e–h) 3D and 2D far field pattern of beam deflection and beam splitting in the lower and higher bands, respectively. (i–l) 3D and 2D far field pattern of quad beam shaping and beam splitting in the lower and higher bands, respectively.

OAM vortex beam generation and control

OAM beams are the conical vortex beams which have the feature of non-diffractiveness and mode orthogonality52–54. Each OAM vortex beam has a mode number ‘±l’ which is orthogonal to any other OAM mode. This feature of the OAM vortex beams make them strong agents to use them as independent channels for transferring information. Theoretically, the OAM provides another degree of freedom to EM waves which can enhance the channel capacity26,55–57.

The phase profile of an OAM vortex beam is described by ‘eilφ’ where ‘l’ is the topological charge (mode number) and ‘φ’ is the azimuth angle54,55. To realize OAM vortex beam from a MTS composed of m×n elements, the required phase distribution at the mnth element should satisfy the following relation with the azimuthal angle5 ϕ(x,y)=l·arctanyx

where l=0,±1,±2,±3… is the OAM mode number and ‘x’, ‘y’ are the Cartesian coordinates in the x–y plane. The relation shown in Eq. (5), is the necessary and sufficient condition to generate an OAM vortex beam of mode ‘±l’. In order to generate OAM beam from the proposed MTS firstly we designed a full structure by selecting the off diagonal elements of the matrix M2-bit and distributed them with the pattern labeled as ‘S4’ depicted in Fig. 5a. This phase pattern fulfills the relation described in Eq. (5). Moreover the gradient of the phase pattern in lower and higher bands are opposite to each other which adds another feature to the design that is the generation of OAM beam in the two operating bands with mode number opposite to each other, i.e., ‘+l’ at ‘fl’ and ‘-l’ at ‘fh’. In this regard, a full structure composed of 40×40 elements is designed and simulated under the normal incidence of y-polarized plane waves.Fig. 5 OAM beam generation with distinct modes. (a) Coding pattern of the MTS to generate OAM with opposite modes ‘l=+1’ and ‘l=-1’ in lower nad higher bands, respectively. (b, c) Far-field pattern and the corresponding phase of OAM in lower band. (d, e) Far-field pattern and the corresponding phase of OAM in higher band. (f) Coding pattern of the MTS to generate OAM with distinct modes ‘l=+1’ and ‘l=+2’ in lower and higher bands, respectively. (g, h) Far-field pattern and the corresponding phase of OAM in lower band. (i, j) Far-field pattern and the corresponding phase of OAM in higher band.

The simulation results shown in Fig. 5b and Fig. 5d show the 2D far-field pattern in lower and higher bands, respectively, where a donut-like pattern with a null in the center can be observed which satisfies the requirements of an OAM beam. Figure 5c and Fig. 5e show the corresponding OAM phase patterns in the lower and higher bands, respectively. The phase patterns clearly show that the OAM beam generated has mode number ‘l=+1’ and ‘l=-1’ in the two operating bands, respectively. Consequently, two independent channels are established using the same design. The idea is further applied to design a full structure to generate OAM beams with two distinct integral values of topological charge l. For this purpose, we selected the coding meta-atoms which has a phase distribution shown in Fig. 5f, i.e., a single repetition of the digital phase state in the lower band and a double repeating pattern of the same digital phases in the higher band. The meta-atoms which fulfill this condition have the digital phase states of ‘00/00, 00/01, 01/10, 01/11, 10/00, 10/01, 11/10, 11/11’, which have decimal equivalents ‘0/0, 0/1, 1/2, 1/3, 2/0, 2/1, 3/2, 3/3’. A full structure labeled as ‘S5’, composed of 64×64 coding meta-atoms are distributed with the pattern as shown in Fig. 5f. The simulation results are shown in Fig. 5g–j, where Fig. 5g and Fig. 5i show the 2D far-field patterns in the lower and higher bands, respectively. Figure 4h and j shows that the generated OAM has mode numbers ‘l=1’ and ‘l=2’, respectively. Furthermore, the idea is extended to design a full structure that can steer the OAM beams in different directions independently. For this purpose, convolution operation on coding MTSs is employed to accomplish the task. Note that convolution operation on coding MTSs allows us to simply make binary addition of two or more different coding patterns that result in a new useful coding pattern where the resultant pattern can simultaneously perform all the tasks of all the individual coding patterns. For this purpose, we selected the phase states lying on the off-diagonal of the Matrix M2-bit and distributed them in the two patterns as shown in Fig. 6a and Fig. 6b, both having the same size of 64×64 coding meta-atoms. As mentioned earlier the coding pattern shown in Fig. 6a is responsible to generate the OAM beam in the normal direction with mode ‘l=+1’ and ‘l=-1’ in the two operating bands, respectively, while the coding pattern shown in Fig. 6b can be employed to deflect the reflected beam to opposite half-planes for the two frequencies of operation. Now performing the binary addition of these two coding patterns we get the new coding pattern as shown in Fig. 6c, which can perform both the tasks of OAM beam generation and deflection, simultaneously. A full structure labeled as ‘S6’ is designed according to this new coding pattern and simulated in both operating bands under the normal illumination of x-polarized plane waves. The simulation results are shown in Fig. 6d–i, where Fig. 6d and Fig. 6g show the 3D far-field patterns of the steered OAM beam, while Fig. 6e and h shows their corresponding 2D far-field patterns in the two operating bands, respectively. The deflection angle has almost the same value as predicted by the Eq. (2), which in other words verify that the simulated results have excellent agreement with theoretical results. Similarly, the phase patterns in the two operating bands for the steered OAM beams are depicted in Fig. 6f and Fig. 6i, respectively, showing that the OAM beam of mode ‘l=+1’ and ‘l=-1’ have been generated.Fig. 6 OAM beam generation with steering features. (a) Coding pattern of the MTS to generate OAM in the normal direction with opposite modes ‘l=1’ and ‘l=-1’ in lower and higher bands, respectively. (b) Coding pattern of the MTS to steer the reflected beams to opposite half planes. (c) Coding pattern of the MTS to generate and steer the OAM beams to opposite half planes. (d, e) 3D and 2D far field pattern of the generated OAM beam in the lower band, respectively (f) The phase pattern of the generated OAM beam in the lower band. (g, h) 3D and 2D far field pattern of the generated OAM beam in the higher band, respectively. (i) The phase pattern of the generated OAM beam in the higher band.

OAM mode purity

An important feature of OAM beams generated by a MTS is their mode purity, which is desired in most practical scenarios so that information can be transferred without distortion. To analyze the performance of the designed coding MTS in generation of OAM beams with different modes, the purity of different OAM modes can be extracted based on the spectral analysis of Fourier transform58. For this purpose, we analyzed the performance of the structures ‘S4–S6’ that generate OAM beams with ‘l=+1’,‘l=-1’ (design S4), ‘l=+1’,‘l=+2’(design S5) in the specular direction and ‘l=+1’,‘l=-1’ (design S6) deflected in the lower and higher band, respectively. By sampling the fields along a circle with a specific radius, ξ(ϕ), the Fourier relation can be written as6 ξ(ϕ)=∑nane-inϕ

where the amplitude ‘an’ can be calculated as7 an=12π∫02πξ(ϕ)e-inϕdϕ

In Eq. (7), ‘einϕ’ indicates the spiral harmonics of each mode. Utilizing the above two equations the normalized spectral intensities of different OAM harmonics are calculated and displayed in Fig. 7. In the lower operating band the normalized spectral intensities of the OAM modes ‘l=+1’ and ‘l=-1’, are shown in Fig. 7a and b,respectively, while Fig. 7c and d shows the same for modes ‘l=+1’ and ‘l=+2’, in the lower and higher band, respectively. It should be noted that compared to lower band the intensities of the unwanted higher order modes in the high frequency band are relatively higher which are due to the increased losses induced in the structures at the higher frequency. Moreover, the OAM mode purity of the deflected beams generated by ‘S6’ are relatively lower as compared to the other two structures which is due to the increased unwanted coupling between the adjacent unit cells resulted from the addition of two different coding sequences. The results in this case for ‘l=+1’ and ‘l=-1’, are displayed in Fig. 7e and f, in the lower and higher band, respectively. From these results it can be seen that the proposed designs have generated OAM beams with high purity since the dominant OAM modes carry ‘>87%’ of the total energy for all the modes ‘l=+1’ and ‘l=-1’ of (S4) ‘>85%’ for ‘l=+1’ and ‘l=+2’ of (S5) and ‘>80%’ of the total energy for all the modes ‘l=+1’ and ‘l=-1’ of the structure (S6). The generation of OAM beams with high purity shows their applicability for mode division multiplexing communication meta-devices.Fig. 7 OAM beam mode purity. (a, b) Normalized intensity of the OAM beams with mode ‘l=+1’ and ‘l=-1’ in lower and higher band, respectively. (c, d) Normalized intensity of the OAM beams with mode ‘l=+1’ and ‘l=+2’ in lower and higher band, respectively. (e, f) Normalized intensity of the deflected OAM beams with mode ‘l=+1’ and ‘l=-1’ in lower and higher band, respectively.

PW–SW conversion and PW beam manipulation

It has been established previously that a gradient MTS and digital/coding MTS as well can convert the PW to SW when it fulfills the required conditions5,59,60. The necessary condition for phase modulated coding MTS to convert the PW–SW is to select a gradient coding sequence that has a period length equal to or less than the operating wavelength, i.e., ‘λ≥Γ’. The proposed design which has the coding meta-atoms with the feature of deep sub-wavelength size, provides the flexibility to convert PW–SW at least at lower band operating frequencies and performs PW beam manipulation in higher band,as illustrated in Fig. 8a. Moreover, the direction of the flow of the SWs can also be controlled by changing the coding sequence of the design. In this regard, various structures are designed to convert the PW–SW in the lower band and perform other useful tasks at the higher band independently. Firstly, the coding meta-atoms which have the digital phase states lie on the off-diagonal of the Matrix M2-bit, i.e., ‘00/11, 01/10, 10/01, 11/00’, with a period of Γ=4×70μm, are selected. The length of this period is less than the operating wavelength (λo>Γ) in the lower band. The operating wavelength here we consider is the wavelength corresponding to 0.63 THz. A full structure labeled ‘S7’ is designed with an overall size of 64×64 coding elements and simulated by adjusting the waveguide port with an aperture size of 3200×2800μm2, 1000μm above the design as an excitation source so that the design can receive the plane waves. Both near and far field monitors are set at the selected frequencies in both lower and higher bands to observe the near field distributions and far-field patterns.Fig. 8 Simultaneous manipulation of surface waves and propagating waves. (a) Schematic illustration of PW–SW conversion in lower band and PW beam manipulation in higher band, (b, c) SWs received from PW–SW converter and beam deflection in lower and higher bands, respectively, employing off diagonal elements of the matrix M2-bit. (d, e) SWs received from PW–SW converter and beam deflection in lower and higher band, respectively using diagonal elements of the matrix M2-bit. (f, g) SWs received from PW–SW converter flowing in both ±x directions and beam deflection in and lower and higher band, respectively using elements, ‘00/00, 00/01, 10/10, 10/11’. (h, i) SWs received from PW–SW converter and normal reflection in lower and higher band, respectively using first row elements of the matrix M2-bit.

To record the SWs quartz slabs with permittivity, ε=4.4 and loss tangent δ=0.0004, having a size of 4480×2240×80um3 is placed at each side (‘±x’ and ‘±y’) directions of the PW–SW converter. As it is obvious that the design ‘S7’ has a period length of the coding sequence (Γ<λo) at the lower band, and (Γ>λo) at the higher band, consequently, the design should convert the PW–SW at lower band and steer PW beam at higher band. The simulation results displayed in Fig. 8b and Fig. 8c show the above predictions where in the lower band the design has efficiently converted PW–SW which flows to-x direction on the quartz substrate (see Fig. 8b), while in the higher band (Fig. 8c) the structure deflected the beam to (x-z plane) and the angle, ‘θd≈26∘’ has a close agreement with the value calculated by Eq. (2). In the second case, we designed a full structure labeled as ‘S8’ to convert PW–SW in the lower band and scatter the PW reflected beam to the same plane. This function can be achieved using coding elements lying on the diagonal of Matrix M2-bit, i.e., 00/00, 01/01, 10/10, and11/11. The results are shown in Fig. 8d and Fig. 8e which shows SW flowing on the dielectric substrates and PW deflected to an angle ‘θd≈-26∘’ at lower and higher band, respectively. To further explore the PW–SW conversion of the proposed MTS. In the third case, we designed a full structure labeled ‘S9’ that can convert the PW–SW in the lower band which flows in both ±x direction and steer the propagation waves in the higher band. For this purpose, we designed a full structure that has a beam-splitting coding sequence in the lower band and beam deflection (‘00/00, 00/01, 10/10, 10/11’) in the higher band. Here the period length is Γ=4×70um and the full structure size again is 64×64 meta-atoms. The corresponding results are shown in Fig. 8f and Fig. 8g for lower and higher bands, respectively. In the last case, another full structure labeled as ‘S10’ is constructed to scatter the reflected beam in the normal direction and to convert PW–SW in the lower band. The coding sequence for this purpose can be used is any row of the matrix M2-bit, however, for simplicity we selected the first row with phase states, ‘00/00, 01/00, 10/00, 11/00’. Clearly, this sequence has a gradient phase pattern in the lower and a constant phase pattern in the higher band which fulfills the required condition for this case. The simulation results are plotted in Fig. 8h and Fig. 8i for lower and higher bands, respectively, clearly showing the SW propagating to-x direction and PW reflected normally.

To prove the novelty and strength of the proposed design compared to the previous works, a brief comparison has been made that can be seen in Table 1. The comparison clearly shows that the proposed design is superior to previously demonstrated THz works in many aspects, including phase resolution in both operating bands, polarization insensitivity and the number of achievable bi-functionalities.Table 1 Comparison of the proposed work with previously demonstrated THz works.

Refs	Operating bands	Phase resolution	Polarization insensitivity	Band width in (THz)	Polarization	Bifunctionalities	
34	1	2-bit	No	0.8–1.3	Linear	Beam Shaping and deflection	
51	2	1-bit in L.B, 2-bit in H.B	No	0.75–0.9 in L.B and 1–1.2 in H.B	Linear	Beam splitting, deflection	
59	1	3-bit	No	0.6–1.4	Linear	Beam deflection, PW–SW conversion	
61	1	1-bit	No	0.38–43	Linear	Beam shaping, OAM beam generation	
62	1	3-bit	No	2.4–4.8	Circular	Beam shaping, OAM beam generation	
63	1	3-bit	No	0.8–1.8	Circular	Beam shaping/deflection, RCS reduction	
This work	2	2-bit in L.B and 2-bit in H.B	Yes	0.53–0.73 in L.B and 1.2–1.45 in H.B	Linear	Beam shaping/deflection, OAM, PW–SW conversion and all others 2-bit bi-functionalities	
L.B, H.B stands for lower band and higher band, respectively

Efficiency and angular stability analysis

The efficiency of the MTS for a specific functionality at a particular frequency can be calculated as:8 η=PoutPin×100%

In Eq. (8), Pout and Pin shows the scattered output and input power, respectively. For a reflection type MTS the efficiency can be estimated as the ratio of power reflected by MTS to the power reflected by a bare metal of the same size, i.e.,9 ηr=PrMTSPrMetal×100%

where PrMTS is the power reflected in the presence of MTS while PrMetal the power reflected by the bare metallic structure.

To numerically calculate the efficiency of the proposed design, for simplicity we have simulated the structure shown in Fig. 4a and Fig. 4c which has the functionality of beam splitting and beam deflection in the lower and higher band, respectively, as shown in Supplementary Fig. S4(a) and Fig. S4(c). Next, a bare metallic structure of the same size is also simulated and the reflected power at the same central frequencies are recorded and plotted (see Supplementary Fig. S4b,d), which shows the specular reflection in both bands. Furthermore, to compare the reflected power and calculate the efficiency, the reflected power are plotted in lower and higher band, for both MTS and bare metal case that can be seen in (Supplementary Fig. S1e,f). The comparison of the two plots shows that the proposed design for this functionality has the efficiency ηr=34% in lower band and ηr=85% in the higher band. The lower efficiency in lower band is due to the fact that the energy is distributed into two equal parts.

Furthermore for angular stability analysis of the proposed work, we designed a full structure composed of the off diagonal elements ‘00/11, 01/10, 10/01, 11/00’ of the Matrix M2-bit and simulated it for plane waves incidence with different oblique angles. The structure at the normal incidence ‘θinc=0∘’, deflects the reflected beam to the opposite half planes as shown in (Supplementary Fig. S5a,b). The same structure is further simulated for two distinct incidence angles including ‘θinc=15∘’ and ‘θinc=30∘’, as shown (Supplementary Fig. S5c–f) which shows the reflected E-field pattern in the lower and higher band for ‘θinc=15∘’ and ‘θinc=30∘’, respectively, confirming that the proposed design still works efficiently for oblique incidence, especially from the functionality point of view. Increasing the oblique incidence further the performance of the MTS degraded, due to the growing unwanted coupling between the meta-atoms and also possibly due to the in destructive interference between the incidence and reflected waves.

Conclusion

In summary, we have presented a bi-functional polarization insensitive dual-band digital/coding MTS to flexibly manipulate the reflected THz waves independently in the two operating bands. Various full structures are designed and encoded with different coding sequences to perform two independent tasks from the same shared aperture. A series of two independent tasks reported here from the proposed frequency-dependent meta-devices, including beam shaping and beam deflection, beam deflection to opposite half-planes, high mode purity OAM beam generation and deflecting to different directions with distinct mode numbers and lastly PW–SW conversion at the lower band and PW beam manipulation at the higher frequency band. All the simulated trends match well with their theoretical predictions. The proposed design is superior compared to the previously frequency dependent MTSs as many other bi-functionalities could also be achieved by the proper distribution of the designed meta-atoms. We believe that the proposed idea has the potential to design polarization-insensitive multi-spectral multipurpose meta-devices in the THz regime, which could find application in multi-color holography and frequency beam manipulators for the spatial separation of signals modulated with different frequencies.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-71910-5.

Author contributions

S.I.: Writing—review and editing, writing—original draft, visualization, validation, software, project administration, investigation, formal analysis, data curation, conceptualization. A.N.: Writing review and editing, software, formal analysis, data curation. N.U.: Writing—review and editing, methodology, formal analysis, data curation. Y.S.: Writing—review and editing, methodology, formal analysis, data curation. S.A.: Writing—review and editing, formal analysis, data curation. M.S.N.: Writing—review and editing, formal analysis, data curation. S.-W.W.: Writing—review and editing, supervision, project administration, methodology, funding acquisition, formal analysis, data curation.

Funding

This work is partly supported by, National Natural Science Foundation of China under grant number. 62171289 and partly supported by, State Key Laboratory of Radio Frequency Heterogeneous Integration, Shenzhen university, Shenzhen China (Independent Scientific Research Program number. 2024012).

Data availability

All the authors agreed that they have no competing interests. The authors declare that the data supporting the findings of this study are available within the paper. Any other simulation data sets are available from the corresponding author on a reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Yu N Light propagation with phase discontinuities: generalized laws of reflection and refraction Science 2011 334 333 337 10.1126/science.1210713 21885733
Yu, N. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science 334, 333–337 (2011).21885733 10.1126/science.1210713
2. Yu N A broadband, background-free quarter-wave plate based on plasmonic metasurfaces Nano Lett. 2012 12 6328 6333 10.1021/nl303445u 23130979
Yu, N. et al. A broadband, background-free quarter-wave plate based on plasmonic metasurfaces. Nano Lett. 12, 6328–6333 (2012).23130979 10.1021/nl303445u
3. Aieta F Aberration-free ultrathin flat lenses and axicons at telecom wavelengths based on plasmonic metasurfaces Nano Lett. 2012 12 4932 4936 10.1021/nl302516v 22894542
Aieta, F. et al. Aberration-free ultrathin flat lenses and axicons at telecom wavelengths based on plasmonic metasurfaces. Nano Lett. 12, 4932–4936 (2012).22894542 10.1021/nl302516v
4. Grady NK Terahertz metamaterials for linear polarization conversion and anomalous refraction Science 2013 340 1304 1307 10.1126/science.1235399 23686344
Grady, N. K. et al. Terahertz metamaterials for linear polarization conversion and anomalous refraction. Science 340, 1304–1307 (2013).23686344 10.1126/science.1235399
5. Sun S Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves Nat. Mater. 2012 11 426 431 10.1038/nmat3292 22466746
Sun, S. et al. Gradient-index meta-surfaces as a bridge linking propagating waves and surface waves. Nat. Mater. 11, 426–431 (2012).22466746 10.1038/nmat3292
6. Lin J Polarization-controlled tunable directional coupling of surface plasmon polaritons Science 2013 340 331 334 10.1126/science.1233746 23599488
Lin, J. et al. Polarization-controlled tunable directional coupling of surface plasmon polaritons. Science 340, 331–334 (2013).23599488 10.1126/science.1233746
7. Huang L Three-dimensional optical holography using a plasmonic metasurface Nat. Commun. 2013 4 2808 10.1038/ncomms3808
Huang, L. et al. Three-dimensional optical holography using a plasmonic metasurface. Nat. Commun. 4, 2808 (2013).10.1038/ncomms3808
8. Chen WT High-efficiency broadband meta-hologram with polarization-controlled dual images Nano Lett. 2014 14 225 230 10.1021/nl403811d 24329425
Chen, W. T. et al. High-efficiency broadband meta-hologram with polarization-controlled dual images. Nano Lett. 14, 225–230 (2014).24329425 10.1021/nl403811d
9. Ni X Kildishev AV Shalaev VM Metasurface holograms for visible light Nat. Commun. 2013 4 2807 10.1038/ncomms3807
Ni, X., Kildishev, A. V. & Shalaev, V. M. Metasurface holograms for visible light. Nat. Commun. 4, 2807 (2013).10.1038/ncomms3807
10. Farmahini-Farahani M Cheng J Mosallaei H Metasurfaces nanoantennas for light processing JOSA B 2013 30 2365 2370 10.1364/JOSAB.30.002365
Farmahini-Farahani, M., Cheng, J. & Mosallaei, H. Metasurfaces nanoantennas for light processing. JOSA B 30, 2365–2370 (2013).10.1364/JOSAB.30.002365
11. Chen H-T Taylor AJ Yu N A review of metasurfaces: physics and applications Rep. Prog. Phys. 2016 79 076401 10.1088/0034-4885/79/7/076401 27308726
Chen, H.-T., Taylor, A. J. & Yu, N. A review of metasurfaces: physics and applications. Rep. Prog. Phys. 79, 076401 (2016).27308726 10.1088/0034-4885/79/7/076401
12. Song M Smart table based on a metasurface for wireless power transfer Phys. Rev. Appl. 2019 11 054046 10.1103/PhysRevApplied.11.054046
Song, M. et al. Smart table based on a metasurface for wireless power transfer. Phys. Rev. Appl. 11, 054046 (2019).10.1103/PhysRevApplied.11.054046
13. Lee D Gwak J Badloe T Palomba S Rho J Metasurfaces-based imaging and applications: From miniaturized optical components to functional imaging platforms Nanoscale Adv. 2020 2 605 625 10.1039/C9NA00751B 36133253
Lee, D., Gwak, J., Badloe, T., Palomba, S. & Rho, J. Metasurfaces-based imaging and applications: From miniaturized optical components to functional imaging platforms. Nanoscale Adv. 2, 605–625 (2020).36133253 10.1039/C9NA00751B
14. Huang CX Zhang J Cheng Q Cui TJ Polarization modulation for wireless communications based on metasurfaces Adv. Funct. Mater. 2021 31 2103379 10.1002/adfm.202103379
Huang, C. X., Zhang, J., Cheng, Q. & Cui, T. J. Polarization modulation for wireless communications based on metasurfaces. Adv. Funct. Mater. 31, 2103379 (2021).10.1002/adfm.202103379
15. Chen H Ran L Huangfu J Grzegorczyk TM Kong JA Equivalent circuit model for left-handed metamaterials J. Appl. Phys. 2006 100 024915 10.1063/1.2219986
Chen, H., Ran, L., Huangfu, J., Grzegorczyk, T. M. & Kong, J. A. Equivalent circuit model for left-handed metamaterials. J. Appl. Phys. 100, 024915 (2006).10.1063/1.2219986
16. Bilotti F Equivalent-circuit models for the design of metamaterials based on artificial magnetic inclusions IEEE Trans. Microw. Theory Tech. 2007 55 2865 2873 10.1109/TMTT.2007.909611
Bilotti, F. et al. Equivalent-circuit models for the design of metamaterials based on artificial magnetic inclusions. IEEE Trans. Microw. Theory Tech. 55, 2865–2873 (2007).10.1109/TMTT.2007.909611
17. Cui TJ Qi MQ Wan X Zhao J Cheng Q Coding metamaterials, digital metamaterials and programmable metamaterials Light: Sci. Appl. 2014 3 e218 e218 10.1038/lsa.2014.99
Cui, T. J., Qi, M. Q., Wan, X., Zhao, J. & Cheng, Q. Coding metamaterials, digital metamaterials and programmable metamaterials. Light: Sci. Appl. 3, e218–e218 (2014).10.1038/lsa.2014.99
18. Ma Q Cui TJ Information metamaterials: Bridging the physical world and digital world PhotoniX 2020 1 1 32 10.1186/s43074-020-00006-w
Ma, Q. & Cui, T. J. Information metamaterials: Bridging the physical world and digital world. PhotoniX 1, 1–32 (2020).10.1186/s43074-020-00006-w
19. Liu S Cui TJ Concepts, working principles, and applications of coding and programmable metamaterials Adv. Opt. Mater. 2017 5 1700624 10.1002/adom.201700624
Liu, S. & Cui, T. J. Concepts, working principles, and applications of coding and programmable metamaterials. Adv. Opt. Mater. 5, 1700624 (2017).10.1002/adom.201700624
20. Bao L Cui TJ Tunable, reconfigurable, and programmable metamaterials Microw. Opt. Technol. Lett. 2020 62 9 32 10.1002/mop.32164
Bao, L. & Cui, T. J. Tunable, reconfigurable, and programmable metamaterials. Microw. Opt. Technol. Lett. 62, 9–32 (2020).10.1002/mop.32164
21. Liu S Convolution operations on coding metasurface to reach flexible and continuous controls of terahertz beams Adv. Sci. 2016 3 1600156 10.1002/advs.201600156
Liu, S. et al. Convolution operations on coding metasurface to reach flexible and continuous controls of terahertz beams. Adv. Sci. 3, 1600156 (2016).10.1002/advs.201600156
22. Cui T-J Liu S Li L-L Information entropy of coding metasurface Light: Sci. Appl. 2016 5 e16172 e16172 10.1038/lsa.2016.172 30167131
Cui, T.-J., Liu, S. & Li, L.-L. Information entropy of coding metasurface. Light: Sci. Appl. 5, e16172–e16172 (2016).30167131 10.1038/lsa.2016.172
23. Iqbal S Qian M Ullah N Noor A Saifullah Y Saleem A Wong S-W Dual-band and dual-polarization coding metasurface for independent controls of phases and amplitudes with cross-polarized transmissions in two bands Results Phys. 2024 58 107452 10.1016/j.rinp.2024.107452
Iqbal, S. et al. Dual-band and dual-polarization coding metasurface for independent controls of phases and amplitudes with cross-polarized transmissions in two bands. Results Phys. 58, 107452 (2024).10.1016/j.rinp.2024.107452
24. Liu C A programmable diffractive deep neural network based on a digital-coding metasurface array Nat. Electron. 2022 5 113 122 10.1038/s41928-022-00719-9
Liu, C. et al. A programmable diffractive deep neural network based on a digital-coding metasurface array. Nat. Electron. 5, 113–122 (2022).10.1038/s41928-022-00719-9
25. Iqbal S Polarization-selective dual-band digital coding metasurface for controls of transmitted waves J. Phys. D Appl. Phys. 2018 51 285103 10.1088/1361-6463/aacad5
Iqbal, S. et al. Polarization-selective dual-band digital coding metasurface for controls of transmitted waves. J. Phys. D Appl. Phys. 51, 285103 (2018).10.1088/1361-6463/aacad5
26. Iqbal S Madni HA Liu S Zhang L Cui TJ Full controls of OAM vortex beam and realization of retro and negative reflections at oblique incidence using dual-band 2-bit coding metasurface Mater. Res. Express 2019 6 125804 10.1088/2053-1591/ab5ca9
Iqbal, S., Madni, H. A., Liu, S., Zhang, L. & Cui, T. J. Full controls of OAM vortex beam and realization of retro and negative reflections at oblique incidence using dual-band 2-bit coding metasurface. Mater. Res. Express 6, 125804 (2019).10.1088/2053-1591/ab5ca9
27. Iqbal S Power modulation of vortex beams using phase/amplitude adjustable transmissive coding metasurfaces J. Phys. D Appl. Phys. 2020 54 035305 10.1088/1361-6463/abbcfc
Iqbal, S. et al. Power modulation of vortex beams using phase/amplitude adjustable transmissive coding metasurfaces. J. Phys. D Appl. Phys. 54, 035305 (2020).10.1088/1361-6463/abbcfc
28. Zhang K Phase-engineered metalenses to generate converging and non-diffractive vortex beam carrying orbital angular momentum in microwave region Opt. Express 2018 26 1351 1360 10.1364/OE.26.001351 29402010
Zhang, K. et al. Phase-engineered metalenses to generate converging and non-diffractive vortex beam carrying orbital angular momentum in microwave region. Opt. Express 26, 1351–1360 (2018).29402010 10.1364/OE.26.001351
29. Zhang L Space-time-coding digital metasurfaces Nat. Commun. 2018 9 4334 10.1038/s41467-018-06802-0 30337522
Zhang, L. et al. Space-time-coding digital metasurfaces. Nat. Commun. 9, 4334 (2018).30337522 10.1038/s41467-018-06802-0
30. Li L Electromagnetic reprogrammable coding-metasurface holograms Nat. Commun. 2017 8 197 10.1038/s41467-017-00164-9 28775295
Li, L. et al. Electromagnetic reprogrammable coding-metasurface holograms. Nat. Commun. 8, 197 (2017).28775295 10.1038/s41467-017-00164-9
31. Ma Q Smart metasurface with self-adaptively reprogrammable functions Light: Sci. App. 2019 8 98 10.1038/s41377-019-0205-3
Ma, Q. et al. Smart metasurface with self-adaptively reprogrammable functions. Light: Sci. App. 8, 98 (2019).10.1038/s41377-019-0205-3
32. Iqbal S Dual-band 2-bit coding metasurface for multifunctional control of both spatial waves and surface waves JOSA B 2019 36 293 299 10.1364/JOSAB.36.000293
Iqbal, S. et al. Dual-band 2-bit coding metasurface for multifunctional control of both spatial waves and surface waves. JOSA B 36, 293–299 (2019).10.1364/JOSAB.36.000293
33. Liu S Frequency-dependent dual-functional coding metasurfaces at terahertz frequencies Adv. Opt. Mater. 2016 4 1965 1973 10.1002/adom.201600471
Liu, S. et al. Frequency-dependent dual-functional coding metasurfaces at terahertz frequencies. Adv. Opt. Mater. 4, 1965–1973 (2016).10.1002/adom.201600471
34. Liu S Anisotropic coding metamaterials and their powerful manipulation of differently polarized terahertz waves Light: Sci. Appl. 2016 5 e16076 e16076 10.1038/lsa.2016.76 30167164
Liu, S. et al. Anisotropic coding metamaterials and their powerful manipulation of differently polarized terahertz waves. Light: Sci. Appl. 5, e16076–e16076 (2016).30167164 10.1038/lsa.2016.76
35. Wu LW Multitask bidirectional digital coding metasurface for independent controls of multiband and full-space electromagnetic waves Nanophotonics 2022 11 2977 2987 10.1515/nanoph-2022-0190
Wu, L. W. et al. Multitask bidirectional digital coding metasurface for independent controls of multiband and full-space electromagnetic waves. Nanophotonics 11, 2977–2987 (2022).10.1515/nanoph-2022-0190
36. Xu H-X Chirality-assisted high-efficiency metasurfaces with independent control of phase, amplitude, and polarization Adv. Opt. Mater. 2019 7 1801479 10.1002/adom.201801479
Xu, H.-X. et al. Chirality-assisted high-efficiency metasurfaces with independent control of phase, amplitude, and polarization. Adv. Opt. Mater. 7, 1801479 (2019).10.1002/adom.201801479
37. Liu L Broadband metasurfaces with simultaneous control of phase and amplitude Adv. Mater. 2014 26 5031 5036 10.1002/adma.201401484 24863731
Liu, L. et al. Broadband metasurfaces with simultaneous control of phase and amplitude. Adv. Mater. 26, 5031–5036 (2014).24863731 10.1002/adma.201401484
38. Wu RY Shi CB Liu S Wu W Cui TJ Addition theorem for digital coding metamaterials Adv. Opt. Mater. 2018 6 1701236 10.1002/adom.201701236
Wu, R. Y., Shi, C. B., Liu, S., Wu, W. & Cui, T. J. Addition theorem for digital coding metamaterials. Adv. Opt. Mater. 6, 1701236 (2018).10.1002/adom.201701236
39. He J Dong T Chi B Zhang Y Metasurfaces for terahertz wavefront modulation: A review J. Infrared Millim. Terahertz Waves 2020 41 607 631 10.1007/s10762-020-00677-3
He, J., Dong, T., Chi, B. & Zhang, Y. Metasurfaces for terahertz wavefront modulation: A review. J. Infrared Millim. Terahertz Waves 41, 607–631 (2020).10.1007/s10762-020-00677-3
40. Tian HW Terahertz metasurfaces: Toward multifunctional and programmable wave manipulation Front. Phys. 2020 8 584077 10.3389/fphy.2020.584077
Tian, H. W. et al. Terahertz metasurfaces: Toward multifunctional and programmable wave manipulation. Front. Phys. 8, 584077 (2020).10.3389/fphy.2020.584077
41. Al-Naib I Withayachumnankul W Recent progress in terahertz metasurfaces J. Infrared Millim. Terahertz Waves 2017 38 1067 1084 10.1007/s10762-017-0381-2
Al-Naib, I. & Withayachumnankul, W. Recent progress in terahertz metasurfaces. J. Infrared Millim. Terahertz Waves 38, 1067–1084 (2017).10.1007/s10762-017-0381-2
42. Scherger B Compression molded terahertz transmission blaze-grating IEEE Trans. Terahertz Sci. Technol. 2012 2 556 561 10.1109/TTHZ.2012.2210892
Scherger, B. et al. Compression molded terahertz transmission blaze-grating. IEEE Trans. Terahertz Sci. Technol. 2, 556–561 (2012).10.1109/TTHZ.2012.2210892
43. Gospodaric J 3D-printed phase waveplates for THZ beam shaping Appl. Phys. Lett. 2018 112 221104 10.1063/1.5027179
Gospodaric, J. et al. 3D-printed phase waveplates for THZ beam shaping. Appl. Phys. Lett. 112, 221104 (2018).10.1063/1.5027179
44. Liu S Anomalous refraction and nondiffractive Bessel-beam generation of terahertz waves through transmission-type coding metasurfaces ACS Photonics 2016 3 1968 1977 10.1021/acsphotonics.6b00515
Liu, S. et al. Anomalous refraction and nondiffractive Bessel-beam generation of terahertz waves through transmission-type coding metasurfaces. ACS Photonics 3, 1968–1977 (2016).10.1021/acsphotonics.6b00515
45. Rouhi K Rajabalipanah H Abdolali A Multi-bit graphene-based bias-encoded metasurfaces for real-time terahertz wavefront shaping: From controllable orbital angular momentum generation toward arbitrary beam tailoring Carbon 2019 149 125 138 10.1016/j.carbon.2019.04.034
Rouhi, K., Rajabalipanah, H. & Abdolali, A. Multi-bit graphene-based bias-encoded metasurfaces for real-time terahertz wavefront shaping: From controllable orbital angular momentum generation toward arbitrary beam tailoring. Carbon 149, 125–138 (2019).10.1016/j.carbon.2019.04.034
46. Xu J Liu W Song Z Terahertz dynamic beam steering based on graphene coding metasurfaces IEEE Photonics J. 2021 13 1 9
Xu, J., Liu, W. & Song, Z. Terahertz dynamic beam steering based on graphene coding metasurfaces. IEEE Photonics J. 13, 1–9 (2021).
47. Manjappa M Reconfigurable mems fano metasurfaces with multiple-input–output states for logic operations at terahertz frequencies Nat. Commun. 2018 9 4056 10.1038/s41467-018-06360-5 30283070
Manjappa, M. et al. Reconfigurable mems fano metasurfaces with multiple-input–output states for logic operations at terahertz frequencies. Nat. Commun. 9, 4056 (2018).30283070 10.1038/s41467-018-06360-5
48. Cong, L., Pitchappa, P., Wang, N. & Singh, R. Electrically programmable terahertz diatomic metamolecules for chiral optical control. Research (2019).
49. Fu X Yang F Liu C Wu X Cui TJ Terahertz beam steering technologies: from phased arrays to field-programmable metasurfaces Adv. Opt. Mater. 2020 8 1900628 10.1002/adom.201900628
Fu, X., Yang, F., Liu, C., Wu, X. & Cui, T. J. Terahertz beam steering technologies: from phased arrays to field-programmable metasurfaces. Adv. Opt. Mater. 8, 1900628 (2020).10.1002/adom.201900628
50. Yang W Full-space dual-helicity decoupled metasurface for a high-efficiency multi-folded reflective antenna Opt. Express 2022 30 33613 33626 10.1364/OE.471942 36242392
Yang, W. et al. Full-space dual-helicity decoupled metasurface for a high-efficiency multi-folded reflective antenna. Opt. Express 30, 33613–33626 (2022).36242392 10.1364/OE.471942
51. Li S-H Li J-S Frequency coding metasurface for multiple directions manipulation of terahertz energy radiation AIP Adv. 2019 9 035146 10.1063/1.5089617
Li, S.-H. & Li, J.-S. Frequency coding metasurface for multiple directions manipulation of terahertz energy radiation. AIP Adv. 9, 035146 (2019).10.1063/1.5089617
52. Allen L Padgett M Babiker M IV the orbital angular momentum of light Prog. Opt. 1999 39 291 372 10.1016/S0079-6638(08)70391-3
Allen, L., Padgett, M. & Babiker, M. IV the orbital angular momentum of light. Prog. Opt. 39, 291–372 (1999).10.1016/S0079-6638(08)70391-3
53. Iqbal S Controls of transmitted electromagnetic waves for diverse functionalities using polarization-selective dual-band 2 bit coding metasurface J. Opt. 2019 22 015104 10.1088/2040-8986/ab5e18
Iqbal, S. et al. Controls of transmitted electromagnetic waves for diverse functionalities using polarization-selective dual-band 2 bit coding metasurface. J. Opt. 22, 015104. 10.1088/2040-8986/ab5e18 (2019).10.1088/2040-8986/ab5e18
54. Yao A Padgett M Orbital angular momentum: Origins, behavior and applications Adv. Opt. Photonics 2011 3 2 161 204 10.1364/AOP.3.000161
Yao, A. & Padgett, M. Orbital angular momentum: Origins, behavior and applications. Adv. Opt. Photonics 3(2), 161–204 (2011).10.1364/AOP.3.000161
55. He J Generation and evolution of the terahertz vortex beam Opt. Express 2013 21 20230 20239 10.1364/OE.21.020230 24105568
He, J. et al. Generation and evolution of the terahertz vortex beam. Opt. Express 21, 20230–20239 (2013).24105568 10.1364/OE.21.020230
56. Tang S High-efficiency broadband vortex beam generator based on transmissive metasurface Opt. Express 2019 27 4281 4291 10.1364/OE.27.004281 30876045
Tang, S. et al. High-efficiency broadband vortex beam generator based on transmissive metasurface. Opt. Express 27, 4281–4291 (2019).30876045 10.1364/OE.27.004281
57. Iqbal S Broadband and high-efficiency manipulation of transmitted vortex beams via ultra-thin multi-bit transmission type coding metasurfaces IEEE Access 2020 8 197982 197991 10.1109/ACCESS.2020.3033568
Iqbal, S. et al. Broadband and high-efficiency manipulation of transmitted vortex beams via ultra-thin multi-bit transmission type coding metasurfaces. IEEE Access 8, 197982–197991 (2020).10.1109/ACCESS.2020.3033568
58. Shi Y Zhang Y Generation of wideband tunable orbital angular momentum vortex waves using graphene metamaterial reflectarray IEEE Access 2017 6 5341 5347 10.1109/ACCESS.2017.2740323
Shi, Y. & Zhang, Y. Generation of wideband tunable orbital angular momentum vortex waves using graphene metamaterial reflectarray. IEEE Access 6, 5341–5347 (2017).10.1109/ACCESS.2017.2740323
59. Liu S Full-state controls of terahertz waves using tensor coding metasurfaces ACS Appl. Mater. Interfaces 2017 9 21503 21514 10.1021/acsami.7b02789 28580778
Liu, S. et al. Full-state controls of terahertz waves using tensor coding metasurfaces. ACS Appl. Mater. Interfaces 9, 21503–21514 (2017).28580778 10.1021/acsami.7b02789
60. Liu S Negative reflection and negative surface wave conversion from obliquely incident electromagnetic waves Light: Sci. Appl. 2018 7 18008 18008 10.1038/lsa.2018.8 30839515
Liu, S. et al. Negative reflection and negative surface wave conversion from obliquely incident electromagnetic waves. Light: Sci. Appl. 7, 18008–18008 (2018).30839515 10.1038/lsa.2018.8
61. Chen Xi Liu Fei Yang Xiao Jian Fu Jun Wei Wu Lei Zhang Jun Yang Tie Jun Cui Programmable manipulations of terahertz beams by transmissive digital coding metasurfaces based on liquid crystals Adv. Opt. Mater. 2021 9 2100932 10.1002/adom.202100932
Chen Xi, Liu et al. Programmable manipulations of terahertz beams by transmissive digital coding metasurfaces based on liquid crystals. Adv. Opt. Mater. 9, 2100932 (2021).10.1002/adom.202100932
62. Yunping Qi Baohe Zhang Jinghui Ding Ting Zhang Xiangxian Wang Zao Yi Efficient manipulation of terahertz waves by multi-bit coding metasurfaces and further applications of such metasurfaces Chin. Phys. B 2021 30 024211
Yunping, Qi. et al. Efficient manipulation of terahertz waves by multi-bit coding metasurfaces and further applications of such metasurfaces. Chin. Phys. B 30, 024211 (2021).
63. Jiu-sheng Li Jian quan Yao Manipulation of terahertz wave using coding Pancharatnam–Berry phase metasurface IEEE Photonics J. 2018 10 1 12 10.1109/JPHOT.2018.2812238
Jiu-sheng, Li. & Jian quan, Yao. Manipulation of terahertz wave using coding Pancharatnam–Berry phase metasurface. IEEE Photonics J. 10, 1–12 (2018).10.1109/JPHOT.2018.2812238
