Reconfigurable Multifunctional Metasurface Hybridized with Vanadium Dioxide at Terahertz Frequencies

Driven by the continuous demand for system integration and device miniaturization, integrating multiple diversified functions into a single metasurface hybridized with the tunable metaparticle is highly demanding at terahertz (THz) range. However, up to now, because of the limitation of the tunable metaparticle at terahertz range, most of the metasurfaces feature a single function only or process similar functionalities at a single frequency. A reconfigurable multifunctional metasurface which can realize the switch of transmission and reflection and manipulate the linearized polarization state of electromagnetic waves simultaneously over a controllable terahertz frequency range based on the vanadium dioxide was designed for the first time in the paper. The numerical result demonstrates the validity of the appropriately designed metasurface. Simulation results show that the reconfigurable and multifunctional performance of this metasurface can be acquired over 1.59 THz to 1.74 THz without re-optimizing or re-fabricating structures, which effectively extends the operating frequencies. The proposed metasurface holds potential for electromagnetic wave manipulation and this study can motivate the realization of the wideband multifunctional metasurface and the software-driven reconfigurable metasurface at terahertz frequencies.

Previous research was mainly focused on the multifunctional metasurface without tunable metaparticles; however, the multifunctionality can only be acquired at different polarization states, incident angles or wavelengths [20][21][22][23][24][25]. Recent efforts have started to be devoted to the design of the reconfigurable multifunctional metasurface by hybridizing metasurfaces with tunable metaparticles,  Figure 2 shows the illustration of the designed metasurface and expected multiple electromagnetic functionalities. By independently controlling the conductivity of the top and bottom VO2, the metasurface was expected to manipulate the state of incident waves between transmission and reflection, the polarization of the reflected or transmitted waves and dynamically tune the resonant frequency. As shown in Figure 2, a linearly polarized terahertz plane wave normally illuminates on the metasurface along the negative direction of the z-axis. (1) When the VO2 inserted at the gap of bottom metallic gratings was in the insulating state with σ , = 2 × 10 2 S/m, the incident THz wave was rotated by 90° to its orthogonal polarization state after transmitting through the metasurface. (2) When the bottom VO2 was in its fully metallic state with σ , = 2 × 10 5 S/m, the linearly polarized incident wave would be reflected by the metasurface with the same linear polarization. Simultaneously, by varying the conductivity of VO2, which is involved into each metallic resonator on the top of this metasurface between σ , = 2 × 10 2 S/m and 2 × 10 5 S/m, the resonant frequency could be dynamically tuned and thus the working frequency range can be extended. The performance of the proposed metasurface will be discussed and demonstrated in Section 4.  Figure 2 shows the illustration of the designed metasurface and expected multiple electromagnetic functionalities. By independently controlling the conductivity of the top and bottom VO 2 , the metasurface was expected to manipulate the state of incident waves between transmission and reflection, the polarization of the reflected or transmitted waves and dynamically tune the resonant frequency. As shown in Figure 2, a linearly polarized terahertz plane wave normally illuminates on the metasurface along the negative direction of the z-axis. (1) When the VO 2 inserted at the gap of bottom metallic gratings was in the insulating state with σ VO 2 ,bottom = 2 × 10 2 S/m, the incident THz wave was rotated by 90 • to its orthogonal polarization state after transmitting through the metasurface.
(2) When the bottom VO 2 was in its fully metallic state with σ VO 2 ,bottom = 2 × 10 5 S/m, the linearly polarized incident wave would be reflected by the metasurface with the same linear polarization. Simultaneously, by varying the conductivity of VO 2 , which is involved into each metallic resonator on the top of this metasurface between σ VO 2 ,top = 2 × 10 2 S/m and 2 × 10 5 S/m, the resonant frequency could be dynamically tuned and thus the working frequency range can be extended. The performance of the proposed metasurface will be discussed and demonstrated in Section 4.
In order to understand the contribution of each VO 2 and demonstrate the performance of the proposed metasurface, we used air and gold to characterize the insulating state and all-metallic state of VO 2 , respectively, in the simulation. Three other types of unit cells were designed, as shown in Figure 3. Dimensions of these unit cells were the same as that in Figure 1, except the length l in Figure 3c. One could observe that in contrast to the unit cell in Figure 1, the bottom metallic grating with VO 2 was replaced with the metallic grating without VO 2 and metallic substrate, respectively, in Figure 3a,b; the VO 2 at the top layer is removed in Figure 3c. The conductivity of VO2 on the top layer varies between 2 × 10 2 S/m and 2 × 10 5 S/m. By varying the σ , , the resonant frequency is expected to be dynamically tuned. (1) When σ , = 2 × 10 2 S/m, the incident THz wave is rotated by 90° to its orthogonal polarization state after transmitting through the metasurface. (2) When σ , = 2 × 10 5 S/m, the linearly polarized incident wave will be reflected by the metasurface with the same linear polarization.
In order to understand the contribution of each VO2 and demonstrate the performance of the proposed metasurface, we used air and gold to characterize the insulating state and all-metallic state of VO2, respectively, in the simulation. Three other types of unit cells were designed, as shown in Figure 3. Dimensions of these unit cells were the same as that in Figure 1, except the length l in Figure  3c. One could observe that in contrast to the unit cell in Figure 1, the bottom metallic grating with VO2 was replaced with the metallic grating without VO2 and metallic substrate, respectively, in Figure 3a,b; the VO2 at the top layer is removed in Figure 3c.
Based on the above description, it should be clear that VO2 at the bottom layer was used to switch the state of waves between transmission and reflection, and VO2 at the top layer was used to dynamically tune the resonant frequency. Therefore, the contribution of VO2 at the bottom layer can be understood by setting the conductivity of VO2 at the top layer to a certain value, while VO2 at the bottom layer was in the insulating state or its all-metallic state, respectively. Simultaneously, by contrasting the performance of the unit cell in Figure 1 with VO2 at the bottom layer in the insulating state, or its all-metallic state to unit cells in Figure 3a,b, respectively, the performance of transmissionreflection switching of the proposed metasurface could be verified. The contribution of VO2 at the top layer can be understood by comparing the operating frequency of the designed metasurface when the conductivity of VO2 at the top layer is in its insulating or metallic state, respectively. By contrasting the performance of the unit cell in Figure 1 to the unit cell in Figure 3c with l = 70 μm or 77 μm, the operating frequency tunability of the proposed metasurface could be verified. The conductivity of VO 2 on the top layer varies between 2 × 10 2 S/m and 2 × 10 5 S/m. By varying the σ VO 2 ,top , the resonant frequency is expected to be dynamically tuned. (1) When σ VO 2 ,bottom = 2 × 10 2 S/m, the incident THz wave is rotated by 90 • to its orthogonal polarization state after transmitting through the metasurface. (2) When σ VO 2 ,bottom = 2 × 10 5 S/m, the linearly polarized incident wave will be reflected by the metasurface with the same linear polarization.

Mechanism of the Reconfigurable Multifunctional Metasurface
Compared with many complex unit cells of reconfigurable multifunctional metasurfaces, the unit cell of the proposed metasurface is maybe the simplest and most effective structure which can greatly reduce the manufacturing complexity and effectively acquire multifunctionalities. The EM interference model [58,59] was used to illustrate the underlying mechanism of the three-layered reconfigurable multifunctional metasurface and the analytical technique is based on tracking the various Fabry-Perot-like scattering processes within the structures. From Figure 4a, we can see that, when the incident wave i E  illuminates onto the top layer, a portion is reflected by the top layer and is converted to an y-polarized wave or remain an x-polarized wave, while the other portion can propagate through the top layer and illuminates onto the bottom layer. Then a portion of these waves can penetrate the bottom layer with the x-polarized and the y-polarized wave, whereas the other portion is reflected by the bottom layer. A portion of these reflected waves will propagate through the top layer with x-polarized and y-polarized components, and the other portion will be reflected by the top layer and go back to interact with the bottom layer again. It is clear that this process is in an infinite loop.
To verify the working mechanism of the proposed reconfigurable multifunctional terahertz metasurface, as shown in Figure 4b, the unit cell of this metasurface is decomposed into two components: (1) the top layer with a 0.5t-thick lossy polyimide substrate; (2) the bottom layer with a Based on the above description, it should be clear that VO 2 at the bottom layer was used to switch the state of waves between transmission and reflection, and VO 2 at the top layer was used to dynamically tune the resonant frequency. Therefore, the contribution of VO 2 at the bottom layer can be understood by setting the conductivity of VO 2 at the top layer to a certain value, while VO 2 at the bottom layer was in the insulating state or its all-metallic state, respectively. Simultaneously, by contrasting the performance of the unit cell in Figure 1 with VO 2 at the bottom layer in the insulating state, or its all-metallic state to unit cells in Figure 3a,b, respectively, the performance of transmission-reflection switching of the proposed metasurface could be verified. The contribution of VO 2 at the top layer can be understood by comparing the operating frequency of the designed metasurface when the conductivity of VO 2 at the top layer is in its insulating or metallic state, respectively. By contrasting the performance of the unit cell in Figure 1 to the unit cell in Figure 3c with l = 70 µm or 77 µm, the operating frequency tunability of the proposed metasurface could be verified.

Mechanism of the Reconfigurable Multifunctional Metasurface
Compared with many complex unit cells of reconfigurable multifunctional metasurfaces, the unit cell of the proposed metasurface is maybe the simplest and most effective structure which can greatly reduce the manufacturing complexity and effectively acquire multifunctionalities. The EM interference model [58,59] was used to illustrate the underlying mechanism of the three-layered reconfigurable multifunctional metasurface and the analytical technique is based on tracking the various Fabry-Perot-like scattering processes within the structures. From Figure 4a, we can see that, when the incident wave → E i illuminates onto the top layer, a portion is reflected by the top layer and is converted to an y-polarized wave or remain an x-polarized wave, while the other portion can propagate through the top layer and illuminates onto the bottom layer. Then a portion of these waves can penetrate the bottom layer with the x-polarized and the y-polarized wave, whereas the other portion is reflected by the bottom layer. A portion of these reflected waves will propagate through the top layer with x-polarized and y-polarized components, and the other portion will be reflected by the top layer and go back to interact with the bottom layer again. It is clear that this process is in an infinite loop.
To verify the working mechanism of the proposed reconfigurable multifunctional terahertz metasurface, as shown in Figure 4b, the unit cell of this metasurface is decomposed into two components: (1) the top layer with a 0.5t-thick lossy polyimide substrate; (2) the bottom layer with a 0.5t-thick lossy polyimide substrate. Dimensional parameters of this unit were same as those in Figure 1c. Their respective electromagnetic properties and the interaction in between were investigated based on the simulated scattering parameters. In the simulation, periodic boundary conditions were set in the x-direction and y-direction and were open in the z-direction under the condition of free space and the x-polarized and y-polarized incident wave was applied in the unit cell. Their respective 4 × 4 scattering matrices were as follows: where subscripts of transmission coefficients t ji and reflection coefficients r ji indicate the incident polarization i and the transmitted or reflected polarization j. The complex scattering coefficients carry the phase information about the top layer, bottom layer and the polyimide substrate. The arrows atop indicates the incidents propagation along the −z ( It is worth noting that, when the unit cell of this metasurface decomposed into two components, the coupling between the top and bottom layer was not considered, which will cause some error between the simulation and calculation results. Increasing the thickness t of the polyimide substrate would weaken the coupling effect between the top and bottom layer and the error would also be smaller.
It is clear that the process described above is an infinite loop. In addition, a 4 × 4 scattering matrices of two components in Figure 4b can be obtained by simulation. According to the simulation result, we can observe that all of the scattering coefficients used to calculate the overall It is clear that the process described above is an infinite loop. In addition, a 4 × 4 scattering matrices of two components in Figure 4b can be obtained by simulation. According to the simulation result, we can observe that all of the scattering coefficients used to calculate the overall    It is clear that the process described above is an infinite loop. In addition, a 4 × 4 scattering matrices of two components in Figure 4b can be obtained by simulation. According to the simulation result, we can observe that all of the scattering coefficients used to calculate the overall

Results and Discussion
In the simulation, periodic boundary conditions were set in the x and y directions and were open in the z-direction under the condition of free space. The Floquet excitation port was applied  Figure 2. First, the property of transmission-reflection switching and polarization control of this multifunctional metasurface was examined by simulating the unit cells mentioned in Figures 1 and  3a,b. The conductivity of VO 2 at the top layer was assumed as 2 × 10 5 S/m and the VO 2 inserted at the gap of the bottom metallic gratings was in the insulating state and its all-metallic state 2 × 10 2 S/m and 2 × 10 5 S/m, respectively. Figure 6 presents the amplitude of the transmitted coefficient T yx and the reflected coefficient R xx of these three unit cells. As seen from Figure 6a,b, the highest amplitude of T yx was 0.84, while the amplitude of R xx was 0.2 at 1.61 THz with a bottom metallic grating layer. The amplitude of T yx was 0, while the amplitude of R xx was 0.91 at 1.61 THz with a metallic substrate. As shown in Figure 6c,d at the resonant frequency 1.59 THz, when VO 2 at the bottom layer was in the insulating state with σ VO 2 ,bottom = 2 × 10 2 S/m, the amplitude of T yx was 0.82 while the amplitude of R xx was 0.18. When VO 2 at the bottom layer was in the all-metallic state with σ VO 2 ,bottom = 2 × 10 5 S/m, the amplitude of T yx was 0.02 while the amplitude of R xx was 0.89. Thus, it clearly shows that the performance of the VO 2 at the bottom layer with σ VO 2 ,bottom = 2 × 10 2 S/m and 2 × 10 5 S/m was almost the same as the ideal insulating and metallic state.

Results and Discussion
In the simulation, periodic boundary conditions were set in the x and y directions and were open in the z-direction under the condition of free space. The Floquet excitation port was applied in the unit cells, based on the commercial software ANSYS HFSS (2016.1, ANSYS, Canonsburg, PA, USA). The unit cell was normally illuminated by an x-polarized incident light, as illustrated in Figure 2. First, the property of transmission-reflection switching and polarization control of this multifunctional metasurface was examined by simulating the unit cells mentioned in Figures 1 and  3a,b. The conductivity of VO2 at the top layer was assumed as 2 × 10 5 S/m and the VO2 inserted at the gap of the bottom metallic gratings was in the insulating state and its all-metallic state 2 × 10 2 S/m and 2 × 10 5 S/m, respectively. Figure 6 presents the amplitude of the transmitted coefficient Tyx and the reflected coefficient Rxx of these three unit cells. As seen from Figure 6a,b, the highest amplitude of Tyx was 0.84, while the amplitude of Rxx was 0.2 at 1.61 THz with a bottom metallic grating layer. The amplitude of Tyx was 0, while the amplitude of Rxx was 0.91 at 1.61 THz with a metallic substrate. As shown in Figure 6c,d at the resonant frequency 1.59 THz, when VO2 at the bottom layer was in the insulating state with σ , = 2 × 10 2 S/m, the amplitude of Tyx was 0.82 while the amplitude of Rxx was 0.18. When VO2 at the bottom layer was in the all-metallic state with σ , = 2 × 10 5 S/m, the amplitude of Tyx was 0.02 while the amplitude of Rxx was 0.89. Thus, it clearly shows that the performance of the VO2 at the bottom layer with σ , = 2 × 10 2 S/m and 2 × 10 5 S/m was almost the same as the ideal insulating and metallic state. The VO 2 at the bottom layer is in its all-metallic state with σ VO 2 ,bottom = 2 × 10 5 S/m. Therefore, the proposed multifunctional terahertz metasurface cannot only switch the state of waves between transmission and reflection, but also change the polarization of EM waves with high efficiency at the designed frequency.
Then, the operating frequency tunability of the proposed metasurface was investigated by simulating unit cells mentioned in Figures 1 and 3c. The bottom VO 2 underwent the insulator-to-metal transition. The conductivity of the top VO 2 in Figure 1 varied between 2 × 10 2 S/m and 2 × 10 5 S/m. The length of the gold resonator in Figure 3c is 70 µm and 77 µm, respectively.
The simulated T yx and R xx are plotted in Figure 7. It can be found from Figure 7a,b that as the l decreased from 77 µm to 70 µm, the working frequency increased from 1.61 THz to 1.74 THz. Simultaneously, the transmissivity of cross-polarization and reflectivity of co-polarization were greater than 80% or less than 15% in the operating frequency. Figure 7c,d present the T yx and R xx of our designed frequency-tunable reconfigurable multifunctional terahertz metasurface. It can be seen that when σ VO 2 ,bottom = 2 × 10 5 S/m and 2 × 10 2 S/m, the resonant frequency was 1.59 THz and 1.73 THz, respectively. The transmissivity of cross-polarization and reflectivity of co-polarization were greater than 80% or less than 20% in the operating frequency. Figure 6. The simulated amplitude of the transmitted coefficient Tyx (red-solid line) and the reflected coefficient Rxx (blue-dashed line) of unit cells. All of the conductivity of VO2 at the top layer was assumed as σ , = 2 × 10 5 S/m. (a) Metallic gratings without VO2; (b) Metallic substrate; (c) The unit cell of designed multifunctional terahertz metasurface. The VO2 at the bottom layer is in the insulating state with σ , = 2 × 10 2 S/m; (d) The unit cell of designed multifunctional terahertz metasurface. The VO2 at the bottom layer is in its all-metallic state with σ , = 2 × 10 5 S/m. Therefore, the proposed multifunctional terahertz metasurface cannot only switch the state of waves between transmission and reflection, but also change the polarization of EM waves with high efficiency at the designed frequency.
Then, the operating frequency tunability of the proposed metasurface was investigated by simulating unit cells mentioned in Figures 1 and 3c. The bottom VO2 underwent the insulator-tometal transition. The conductivity of the top VO2 in Figure 1 varied between 2 × 10 2 S/m and 2 × 10 5 S/m. The length of the gold resonator in Figure 3c is 70 μm and 77 μm, respectively.
The simulated Tyx and Rxx are plotted in Figure 7. It can be found from Figure 7a,b that as the l decreased from 77 μm to 70 μm, the working frequency increased from 1.61 THz to 1.74 THz. Simultaneously, the transmissivity of cross-polarization and reflectivity of co-polarization were greater than 80% or less than 15% in the operating frequency. Figure 7c,d present the Tyx and Rxx of our designed frequency-tunable reconfigurable multifunctional terahertz metasurface. It can be seen that when σ , = 2 × 10 5 S/m and 2 × 10 2 S/m, the resonant frequency was 1.59 THz and 1.73 THz, respectively. The transmissivity of cross-polarization and reflectivity of co-polarization were greater than 80% or less than 20% in the operating frequency.    Figure 8 plots the T yx , R xx and working frequency for the different σ VO 2 ,top . It can be observed in Figure 8a-c that when σ VO 2 ,top ranges from 2 × 10 5 S/m to 2 × 10 2 S/m with σ VO 2 ,bottom = 2 × 10 2 S/m, T yx changes between 0.52 and 0.82 while R xx varies between 0.11 and 0.45 and the mean value of T yx and R xx is about 0.65 and 0.31 at the working frequency, respectively. As σ VO 2 ,top ranged from 2 × 10 5 S/m to 2 × 10 2 S/m with σ VO 2 ,bottom = 2 × 10 5 S/m at the working frequency, T yx and R xx remained essentially unchanged; the average value was about 0.11 and 0.81, respectively. In Figure 8d, the working frequency of the designed metasurface varies between 1.59 THz and 1.74 THz, while σ VO 2 ,top decreased from 2 × 10 5 S/m to 2 × 10 2 S/m. Therefore, it is clear that the working frequency can be effectively tuned between 1.59 THz and 1.74 THz by varying σ VO 2 ,top without re-optimizing and re-fabricating the structures of the metasurface. Simultaneously, the high performance of diversified functionalities with transmission-reflection switching and polarization control can be maintained across 1.59-1.62 THz and 1.71-1.74 THz.
and Rxx is about 0.65 and 0.31 at the working frequency, respectively. As σ , ranged from 2 × 10 5 S/m to 2 × 10 2 S/m with σ , = 2 × 10 5 S/m at the working frequency, Tyx and Rxx remained essentially unchanged; the average value was about 0.11 and 0.81, respectively. In Figure 8d, the working frequency of the designed metasurface varies between 1.59 THz and 1.74 THz, while σ , decreased from 2 × 10 5 S/m to 2 × 10 2 S/m. Therefore, it is clear that the working frequency can be effectively tuned between 1.59 THz and 1.74 THz by varying σ , without re-optimizing and re-fabricating the structures of the metasurface. Simultaneously, the high performance of diversified functionalities with transmissionreflection switching and polarization control can be maintained across 1.59-1.62 THz and 1.71-1.74 THz.     xx . Increasing the thickness t of the polyimide substrate will weaken the coupling effect between the top and bottom layer, which will also decrease the error.
Thus, it is clear that two diversified functions were integrated into the appropriately designed metasurface, and the working frequency can be dynamically tuned which effectively extend operating frequencies. The multifunctionality and operating frequency tunability can be realized by simply varying the conductivity of VO 2 without re-optimizing or re-fabricating structures of the metasurface. The unit cell of the proposed metasurface consists of three simple layers. Compared with many complex unit cells of reconfigurable multifunctional metasurfaces, the proposed unit cell is maybe the simplest and most effective structure that can greatly reduce the manufacturing complexity and effectively acquire multifunctionalities. The proposed metasurface holds great potential for EM wave manipulation and can motivate the realization of the wideband multifunctional metasurface and the software-driven reconfigurable metasurface, which has the prospect to conveniently realize a complex system integration and device miniaturization with low costs at THz frequencies. This study can pave the way to many practical applications such as telecommunications, sensing and diagnostics, nanoelectronics, antennas and automotive. t of the polyimide substrate will weaken the coupling effect between the top and bottom layer, which will also decrease the error. Thus, it is clear that two diversified functions were integrated into the appropriately designed metasurface, and the working frequency can be dynamically tuned which effectively extend operating frequencies. The multifunctionality and operating frequency tunability can be realized by simply varying the conductivity of VO2 without re-optimizing or re-fabricating structures of the metasurface. The unit cell of the proposed metasurface consists of three simple layers. Compared with many complex unit cells of reconfigurable multifunctional metasurfaces, the proposed unit cell is maybe the simplest and most effective structure that can greatly reduce the manufacturing complexity and effectively acquire multifunctionalities. The proposed metasurface holds great potential for EM wave manipulation and can motivate the realization of the wideband multifunctional metasurface and the software-driven reconfigurable metasurface, which has the prospect to conveniently realize a complex system integration and device miniaturization with low costs at THz frequencies. This study can pave the way to many practical applications such as telecommunications, sensing and diagnostics, nanoelectronics, antennas and automotive.

Conclusions
In conclusion, we propose a reconfigurable multifunctional terahertz metasurface based on VO2. The designed metasurface can manipulate the linearized polarization state of EM waves and simultaneously realize the switch of transmission and reflection in the designed frequency range by utilizing the insulator-to-metal transition in VO2 inserted at the gap of the bottom metallic gratings; that is, this metasurface can convert incident waves into cross-polarized transmitted waves and copolarized reflected waves in the designed frequency range. In addition, the operating frequency of this metasurface can be effectively tuned in the frequency range of 1.59 THz to 1.74 THz by varying the conductivity of the VO2 loaded on the top gold resonator from 2 × 10 5 S/m to 2 × 10 2 S/m, without re-optimizing and re-fabricating structures of the metasurface, which effectively extends the operating frequencies. The proposed metasurface holds great potential for EM wave manipulation and this study can motivate the realization of the wideband multifunctional metasurface and the

Conclusions
In conclusion, we propose a reconfigurable multifunctional terahertz metasurface based on VO 2 . The designed metasurface can manipulate the linearized polarization state of EM waves and simultaneously realize the switch of transmission and reflection in the designed frequency range by utilizing the insulator-to-metal transition in VO 2 inserted at the gap of the bottom metallic gratings; that is, this metasurface can convert incident waves into cross-polarized transmitted waves and co-polarized reflected waves in the designed frequency range. In addition, the operating frequency of this metasurface can be effectively tuned in the frequency range of 1.59 THz to 1.74 THz by varying the conductivity of the VO 2 loaded on the top gold resonator from 2 × 10 5 S/m to 2 × 10 2 S/m, without re-optimizing and re-fabricating structures of the metasurface, which effectively extends the operating frequencies. The proposed metasurface holds great potential for EM wave manipulation and this study can motivate the realization of the wideband multifunctional metasurface and the software-driven reconfigurable metasurface at THz frequencies, which has the prospect to conveniently realize a complex system integration and device miniaturization with low costs.
The designed metasurface integrates only two expected functions; additionally, the performance still need to be optimized. Therefore, our future work may be focused on optimizing and experimentally verifying the performance of the designed metasurface. Photolithography could be used to fabricate the proposed metasurface, and the resistive heater or the external CW laser could be used to control the temperature of VO 2 . Moreover, an infrared camera could be employed to monitor the temperature. As an emerging research area, integrating multiple diversified functions into a single metasurface over a multiwavelength or wide wavelength range based on tunable metaparticles still requires dealing with formidable challenges at terahertz frequencies.