Polarization and Incident Angle Independent Multifunctional and Multiband Tunable THz Metasurface Based on VO₂

1. Introduction

Terahertz (THz) technology and its devices have sparked curiosity in researcher's minds due to the rapid advancement in THz science and technology. Their broad spectrum of applications in non-intrusive testing, Spectroscopic detection, security screening, sensing, optical imaging, communications, etc. has intensified the allure and importance of their cutting-edge technology [1]. Traditional THz wave transmission regulation techniques rely mostly on phase accumulation. The immediate interaction of THz radiation with natural materials poses a significant challenge due to the characteristics of THz radiation [2]. Metasurfaces (MSs) have developed as an extremely effective solution to overcome this fundamental constraint. Such metasurfaces are characterized by their exceptionally thin metamaterial structures, comprised of planar electromagnetic (EM) microstructures capable of changing the polarization, amplitude, and phase of electromagnetic waves. Research has explored diverse, linearly polarized (LP) converters [3,4,5,6,7] and circularly polarized (CP) converters [8,9,10,11]. There has also been extensively investigated dual-band [12,13,14], multiband [15,16,17,18] and broadband [19,20,21] absorbers. Therefore, there is a strong emphasis on tunable multifunctional devices exploiting metasurfaces consisting of active functional materials.

Polarization refers to light oscillation, commonly explained through the electric field vector [22,23]. This fundamental property of light finds widespread utility across diverse domains, including quantum optics, imaging, optical displays, light-matter interaction, and sensing [24,25,26,27]. Notably, polarization converters demonstrate the ability to manipulate the polarization of THz waves [28]. Such converters facilitate polarization conversion through several modes, encompassing LP to LP [29,30], LP to CP [8,31,32,33], and CP to LP [34,35] conversions.

Another prominent area of research within the domain of metasurfaces is absorbers, which have emerged as significant contributors in the burgeoning field of energy-manipulating devices, primarily owing to their remarkable light-harvesting capabilities. Capitalizing on their intrinsic high-loss properties, these tailored surfaces offer promising prospects for revolutionizing diverse sectors, including solar energy collection [14], thermal imaging technology [15], and the design of ultra-sensitive photodetectors [16]. Landy et al. pioneered the ideal absorber in 2008, ingeniously merging a split-ring resonator with a metal-wire design to achieve outstanding single-band absorption [17]. Nevertheless, the ability to dynamically modulate the properties of metamaterial is essential for their practical applications. Recognizing this, researchers have begun a material exploration, investigating unconventional platforms such as graphene [18,19,20,21], indium tin oxide (ITO) [22], vanadium dioxide (${\mathrm{V}\mathrm{O}}_2$) [23,24], and photosensitive silicon [25,26]. These efforts have resulted in adaptable metasurface absorbers that are not only dual-band but also multi-band.

Previous studies predominantly concentrated on the tunability of a singular functionality, such as single-band absorption [36,37,38]. It also focused on the tunability of broadband polarization conversion for singular functionality in the THz range [39,40,41]. Currently, there is a concerted effort within the scientific community to investigate methodologies for integrating multiple electromagnetic functionalities into a singular device. This pursuit aims to mitigate the manufacturing expenses associated with metasurfaces while broadening their application domains.

So far, limited research endeavors have explored the utilization of metasurfaces to achieve multiband absorption and polarization conversion under linear and circular polarized incident waves while remaining independent of both incident and polarization angles. Moreover, metasurface functional devices face challenges, including low efficiency, restricted bandwidth, lack of tunability, and intricate control methodologies. Hence, our work focuses on designing a multifunctional and multiband metasurface characterized by a simple structure and outstanding performance to facilitate the design of THz photonic devices.

In this paper, we propose a THz multifunctional and multi-band tunable metasurface based on a ${\mathrm{V}\mathrm{O}}_2$-metal hybrid structure. This configuration enables functional switching between single-band and multi-band absorption and polarization conversion. Specifically, when the conductivity of ${\mathrm{V}\mathrm{O}}_2$ is $10 S/m$, it behaves as a dielectric, facilitating multi-band polarization conversion with a polarization conversion ratio (PCR) exceeding 100% for both LP and CP incident waves at frequencies of 0.11 THz and 0.21 THz, respectively. Upon increasing the conductivity of ${\mathrm{V}\mathrm{O}}_2$ to $4\times {10}^3 S/m$, the metasurface transitions to an absorber, achieving an absorption rate of 86% at 0.11 THz. Additionally, when the conductivity of ${\mathrm{V}\mathrm{O}}_2$ reaches $2\times {10}^{5 }S/m$, the proposed design attains absorption across an ultra-broadband range spanning from 0.27 THz to 0.36 THz, achieving a rate of 61%.

Moreover, the proposed design exhibits a stable response to incident and polarization angles up to 90°. The distributions of surface-induced currents on the unit cell and the relative impedance are investigated to elucidate the underlying physical mechanisms. This work represents a significant advancement in THz metasurface technology, showcasing potential applications in light manipulation and laying the groundwork for future developments in photonic device integration.

2. Structure Design and Simulation

To achieve multiple functionalities from a single device, it is essential to meticulously design the proposed metasurface structure to elicit distinct responses. The proposed metasurface is meticulously crafted from various layers of materials to exhibit multifunctional capabilities. By adjusting the conductivities of ${\mathrm{V}\mathrm{O}}_2$material, the topology can be configured to operate as either a multiband polarization converter or multiband absorber for LP and CP incident waves. The schematic representation of the proposed snowflake-shaped metasurface is depicted in Figure 1a. Figure 1b,c show the unit segment's top-view and three-dimensional view.

The optimized geometric measurements of unit cell’s are as follows: P=300 µm, $G_{ }=10\mathrm{\mu }\mathrm{m}$, $L_{ }=15\mathrm{\mu }\mathrm{m}$ $R_{ }=145\mathrm{\mu }\mathrm{m}$, $D_{ }=195.06\mathrm{\mu }\mathrm{m}$ and $h=100\mathrm{\mu }\mathrm{m}$, as depicted in Figure 1c. The lossless polyimide, placed between the upper and bottom layers, has a relative permittivity of 3.5 and a tangent loss is 0.02 [42]. Furthermore, the metal surface is comprised of gold $\left(Au\right)$ with a conductivity of ${4.561x10}^{7}S/m$ [43]. By employing the Drude model, the permittivity of ${\mathrm{V}\mathrm{O}}_2$ is obtain by using equation (1) [44].

$${\epsilon }^{ }\left(\omega \right)={\epsilon }_{\mathrm{\infty }}-\frac{{\omega }_p^2\left(\sigma \right)}{({\omega }^2+{i\gamma \omega )}^{ }}$$

(1)

where, ${\epsilon }_{\mathrm{\infty }}=12$ , ${\gamma }_{ }={5.75\mathrm{x}10}^{13}\mathrm{r}\mathrm{a}\mathrm{d}/\mathrm{s}$, the plasma frequency is obtained at the conductive state ($\sigma$) by using ${\omega }_p^2\left(\sigma \right)=\sigma /{\sigma }_0{\omega }_p^2\left({\sigma }_0\right)$ in which ${\sigma }_0={3x10}^{5}s/m$ and ${\omega }_p^{ }\left({\sigma }_0\right)=1.4x{10}^{15}rad/s$. At $T_C=340K$, convert the status from an insulating to conducting and the conductivity and permittivity fluctuate dramatically across the insulating to conducting state [45]. In this work, we solely investigate the conductive state of ${\mathrm{V}\mathrm{O}}_2$ for the polarization conversion and absorption at $\sigma ={10}^{ }s/m$ and $\sigma ={2x10}^{5}s/m$, respectively.

3. Metasurface Performance

3.1. Performance of the Multifunctional (SFL) Metasurface as a Polarization Conversion

The efficacy of the proposed Snowflake (SFL) shaped metasurfaces can be evaluated through the utilization of a comprehensive full-wave simulation tool, such as CST MW Studio. The unit cell exhibits periodic boundary conditions along the x and y axes, with an open boundary introduced along the z direction to facilitate wave propagation along the z-axis. This arrangement allows for incident terahertz waves with different polarizations in the frequency range of 0.01 THz to 0.4 THz.

In accordance with the principles of polarization, the Jones matrix can be used for corelating incident polarized waves with their reflected counterparts [33]. The relation is expressed as in equation (2).

$$\left(\begin{array}{c}{E}_r^x\\ E_r^y\end{array}\right)=\left(\begin{array}{cc}{R}_{xx}^{ }& R_{xy}^{ }\\ R_{yx}^{ }& R_{yy}^{ }\end{array}\right)\left(\begin{array}{c}{E}_i^y\\ E_i^y\end{array}\right)=R\left(\begin{array}{c}{E}_i^x\\ E_i^Y\end{array}\right)$$

(2)

Herein, R signifies the Cartesian Jones reflection matrix. It operates in collaboration with the incident electric field $E_i^{x\left(y\right)}$ and reflected electric field $E_r^{x\left(y\right)}$, aligned in $x\left(y\right)$ directions. The reflection matrix for circularly polarizations is attained by employing the linear polarized reflection coefficient subsequent to the conversion form Cartesian to circular base.

$$\begin{gathered} R_{CP}=\left(\begin{array}{cc}{R}_{++}& R_{\pm }\\ R_{\mp }& R_{--}\end{array}\right)={\bigwedge }^{-1}R\bigwedge = \\ \frac{1}{2}\left(\begin{array}{cc}{R}_{xx}-R_{yy}-i\left(R_{xy}+{ R}_{yx}\right)& R_{xx}+{ R}_{yy}+i\left(R_{xy}-R_{yx}\right)\\ R_{xx}+{ R}_{yy}-i\left(R_{xy}-{ R}_{yx}\right)& R_{xx}-R_{yy}+i\left(R_{xy}+R_{yx}\right)\end{array}\right) \end{gathered}$$

(3)

Herein, $\bigwedge$ denotes the coordinate transformation matrix is utilized for conversion from Cartesian to circular base, defined as $\bigwedge =\frac{1}{\sqrt{2}}\left(\begin{array}{cc}1& 1\\ i& -i\end{array}\right)$. Here, in equation (3) the sign ‘+’ and ‘-’ indicates the right-handed and left-handed circularly polarized waves, respectively. The unit cell reflection coefficients for normal $x$ and $y$ incidences or RCP and LCP polarized waves are depicted in Figure 2, where (${\frac{R_{xx}}{R_{yy}}}_{ })$ ,$({\frac{R_{++}}{R_{--}})}_{ }\mathrm{a}\mathrm{n}\mathrm{d} ({\frac{R_{yx}}{R_{xy}})}_{ }$,$({\frac{R_{-+}}{R_{+-}}}_{ })$ represent the reflection coefficients for co and cross-polarized waves, respectively.

When the conductivity of ${\mathrm{V}\mathrm{O}}_2$ is ${10}^{ }s/m$, the proposed structure act as a multiband polarization conversion for both LP and CP incident waves. In Figure 2 (a) and 2(b), when the incident wave is LP and CP, the cross-polarized reflection coefficient ($R_{yx}$, $R_{xy}$_,$R_{-+}$,$R_{+-}$) is notably observed to attain a value of 100% at 0.11THz and 0.26 THz, also achieved more than 50% values in bandwidth (0.10 THz – 0.15 THz) and (0.22 THz – 0.30 THz), while co-polarization reflection coefficient reaches ($R_{\mathrm{x}\mathrm{x}}$, $R_{yy}$_,$R_{++}$,$R_{--}$) to 0.1 at frequency 0.11 THz and 0 at 0.26 THz. Hence, the envisioned metasurface configuration adeptly transforms linear and circular polarizations into their corresponding cross-polarizations across multiple bands. Despite the absence of C4 symmetry within the structure configuration, a salient feature lies in its inherent and display mirror symmetry along u-axis, resulting in the equivalence of $R_{xx}=R_{yy}$, $R_{yx}=R_{xy}$, $R_{-+}=R_{+-}$ and $R_{++}=R_{--}$.

The analysis of the proposed metasurface, cross-polarized conversion (CPC) is further elucidated through the calculation of the polarization conversion ratio (PCR), as defined by equation (4).

$$P\mathrm{C}\mathrm{R}=\frac{{\left|R_{cross}\right|}^2}{{\left|R_{cross}\right|}^2+{\left|R_{co}\right|}^2}$$

(4)

Figure 2c,d illustrated the PCR associated with LP and CP, respectively. notably within the frequency range of (0.1 THz – 0.4 THz), the PCR attains 100% efficiency within specific frequency, namely, 0.11 THz and 0.26 THz and more than 90% at frequency intervals (0.1 THz – 0.12THz) and (0.261 THz – 0.276 THz). The polarization conversion ratio is same for both LP and CP incident waves due to unique structure design of unit cell.

In practical applications, the assessment of metasurface performance often necessitates consideration of wide-angle incidence scenarios. Figure 3 depicts the influence of incident LP and CP light at various angles of incidence and azimuth on the polarization conversion effect within metasurface structures. Figure 3a illustrates the consistent Polarization Conversion Ratio (PCR) within an incidence angle range of 0° to 85°. Additionally, Figure 3b provides an examination of diverse azimuthal incidences on PCR, indicating a stable PCR bandwidth within an azimuthal angle range of 0° to 90°. This observed angular stability is ascribed to the diminutive dielectric thickness and unit cell size. Considering the potential for incoming waves to exhibit arbitrary incidence angles in practical scenarios, the metasurface's insensitivity to azimuth and incidence angles renders it a promising candidate for a variety of applications.

Further, the physical mechanism is essential to analyzing the performance of the polarization conversion of the proposed metasurface. Therefore, considering the surface distribution at the different frequencies, namely, 0.106 THz, 0.118 THz, 0.14 THz, 0.244 THz, 0.268THz, and 0.29 THz of the top gold surface of the proposed SFLM structure and the bottom layer (ground plan) of the unit cell for y-polarized incident waves. According to Faraday’s law, a changing magnetic field between the two metals causes surface current to flow in opposite directions on the top and bottom of the metallic layers. The black-colored arrow indicates the net current shown in Figure 4. The equation (5) can be used to determine the relationship between electric and magnetic dipole moments and average electromagnetic fields:

$$\left[\begin{array}{c}p\\ m\end{array}\right]=\left[\begin{array}{cc}{\alpha }_{ee}& {\alpha }_{em}\\ {\alpha }_{me}& {\alpha }_{mm}\end{array}\right]\left[\begin{array}{c}E\\ H\end{array}\right]$$

(5)

Here, $P$=${\left[\begin{array}{cc}{p}_x& p_y\end{array}\right]}^T$,$m={\left[\begin{array}{cc}{m}_x& m_y\end{array}\right]}^T$is the electric dipole moment and magnetic dipole moment, respectively, while the $E$=${\left[\begin{array}{cc}{E}_x& E_y\end{array}\right]}^T$and $H$=${\left[\begin{array}{cc}{H}_x& H_y\end{array}\right]}^T$is the average tangential electric and magnetic fields, respectively at the metasurface. The time changing electric and magnetic surface current polarization will cause electric and magnetic surface currents on the metasurface, which is expressed in the equation:

$$\left[\begin{array}{c}J\\ M\end{array}\right]=i\omega \left[\begin{array}{cc}{\alpha }_{ee}& {\alpha }_{em}\\ {\alpha }_{me}& {\alpha }_{mm}\end{array}\right]\left[\begin{array}{c}E\\ H\end{array}\right]$$

(6)

where $J$=${\left[\begin{array}{cc}{J}_x& J_y\end{array}\right]}^T$ and $M$=${\left[\begin{array}{cc}{M}_x& M_y\end{array}\right]}^T$ electric and magnetic current densities, respectively and $\omega$ is the angular frequency of the incident electromagnetic wave. The relations between the surface current density J and radiated for fields are given by

$$E=-i\frac{\omega \mu }{4\pi }\int J\left(x,y\right)\frac{e^{-ikR}}{R}dxdy$$

(7)

Figure 4. Simulated distributions of the Surface current: (a) and (g) at f=0.106THz, (b) and (h) at f=0.118THz, (c) and (i) at f=0.140THz, (d) and (j) at f=0.268THz, (e) and (k) at f=0.244THz, (f) and (l) at f=0.290THz.

Figure 4. Simulated distributions of the Surface current: (a) and (g) at f=0.106THz, (b) and (h) at f=0.118THz, (c) and (i) at f=0.140THz, (d) and (j) at f=0.268THz, (e) and (k) at f=0.244THz, (f) and (l) at f=0.290THz.

According to Equation (7), for CPC, an electric field that is polarized along the $x$-direction will result in current flow on the metasurfaces in the $y$-direction, while an electric field polarized along the $y$-direction will induce current flow in the $x$-direction. As shown in Figure 4a,g,b,h,c,i demonstrated the cross-polarization characteristics at various operating frequencies. At 1.1 THz, 0.11 THz, and 0.14 THz, the surface current on the top metasurface and the ground plate is inversely parallel, leading to the excitation of magnetic resonance and the generation of an induced magnetic field. This induced magnetic field results in the cross-polarization effect where the reflected THz wave becomes $x$-polarized, as indicated by Equation (6). Similarly, Figure 4d,j illustrated that, at a frequency of 0.24 THz, the surface current on the top layer becomes corresponding to the ground layer, leading to the generation of electric resonance and the formation of an induced electric field. By decomposing the electric field into its orthogonal components $x$ and $y$, it is evident that the electric field along the $x$-axis can cross-couple with the incident electric field to form CPC and result in the reflected THz wave being $x$-polarized according to Equation (5). Figure 4e,k,f,l illustrated that, at frequencies 0.26 THz and 0.29 THz, the primarily distributed surface current on the top layer results in the creation of electrical resonance, forming an electric dipole, where the $x$ component of the induced electric field plays a significant role in generating cross-polarization effect. These observations further validate the cross-polarization characteristics and support the electromagnetic behavior described by the equation, showcasing the potential applications of these findings in THz wave manipulation and control.

Furthermore, the phenomena of cross-polarization conversion can be comprehended through the inherent anisotropy in the designed structure. By rotating the standard $x$-$y \mathrm{t}\mathrm{h}\mathrm{e}$ coordinate system is rotated 45° to establish a unique u-v coordinate system, as shown in Figure 5a, the anisotropic nature of the snowflake-like structure allows the incoming $y$-polarized wave can be divided orthogonally into constituents along the u and v axes. In resulting of this division enables the mathematical representation of the incident and reflected waves as follows [39]:$\overrightarrow{E_i}=\left(\overrightarrow{u}{E}_{ui}+\overrightarrow{v}{E}_{vi}\right)\mathrm{e}\mathrm{x}\mathrm{p}\left(i\phi \right)$ and $\overrightarrow{E_r}=\left(\overrightarrow{u}{R_{uu}E}_{ui}+\overrightarrow{v}{E_{vv}E}_{vi}\right)\mathrm{e}\mathrm{x}\mathrm{p}\left(i\phi \right)$. Here, $R_{uu}=\left|E_{ur}\right|/\left|E_{ui}\right|$ and $R_{vv}=\left|E_{vr}\right|/\left|E_{vi}\right|$ represent the reflection coefficients in the u- and v-directions, respectively, with subscripts $i$ and $r$ indicating incident and reflected waves.

Moreover, the phase disparity $\left(\mathit{\Delta }\mathit{\phi }\right)$ arises from structural asymmetry, leading to the relationship between reflection coefficients can be written as: $\mathit{R}\mathit{v}\mathit{v}={\mathit{R}}_{\mathit{u}\mathit{u}}\mathit{e}\mathit{x}\mathit{p}\left(\mathit{j}\mathit{\Delta }\mathit{\phi }\right).$ Figure 5b depicts the reflection coefficients of the reflected wave, while Figure 5c depicts the phase and phase difference. In the frequency range of 0.15 to 0.4 THz, the reflection coefficients $R_{uu}$ and $R_{vv}$ approaches to 1, with a phase difference of around ±180°. These results indicate that the synthesized waves of $R_{uu}$ and $R_{vv}$ deviate by 90° from the incident wave, showcasing the metasurface's capability to transform y-polarized incident waves into x-polarized reflected waves.

3.2. Performance of the Multifunctional (SFL) Metasurface as an Absorption

Besides the performance of the proposed metasurface in polarization conversion, the multifunctional (SFL) metasurface designed, to operate as an absorber owing to the inherent conductivity of vanadium dioxide (${\mathrm{V}\mathrm{O}}_2$). As depicted in Figure 1, a strategically placed ${\mathrm{V}\mathrm{O}}_2$film on the corner of the snowflake-like structure is instrumental in constituting the multiband absorber.

The absorptivity of the multifunctional (SFL) metasurface is determined through the application of Equation (6).

$$A\left(w\right)=1-R\left(w\right)-T\left(w\right)$$

(6)

Herein, $R={\left|S_{11}\right|}^2$ represent the reflectivity, $T={\left|S_{21}\right|}^2$signifies the transmittivity of the metasurface. It is noteworthy that the gold is utilized as the reflector, resulting the transmission is zero, T = 0. The absorptivity of the metasurface for LP waves is computed as an equation

$$A\left(\omega \right)=1-{\left|R_{xx}\right|}^2-{\left|R_{yx}\right|}^2$$

(7)

And for CP wave,

$$A\left(\omega \right)=1-{\left|R_{\raisebox{1ex}{$++$}\!\left/ \!\raisebox{-1ex}{$--$}\right.}\right|}^2-{\left|R_{-+/+-}\right|}^2$$

(8)

According to equations (7) and (8), when the incident wave is composed of both LP and CP components and given a conductivity of ${\mathrm{V}\mathrm{O}}_2$ as ${4x10}^{3}s/m$, an absorption rate of 86% is attained at 0.11 THz. Moreover, an absorption exceeding 50% is achieved within the frequency bandwidth ranging from 0.09 THz to 0.15 THz, as illustrated in Figure 6a,b.

When the conductivity of ${\mathrm{V}\mathrm{O}}_2$ is increased to ${2x{10}^5}^{ }s/m$, the absorption of 61% is obtained at 0.34 THz and achieved more than 50% absorption in bandwidth (0.27 THz – 0.36THz) for LP and CP incident wave as shown in Figure 6c,d.

Impedance matching theory is an essential principle to evaluate the absorption. In this study, the absorption and relative impedance is achieved as [46].

$$A=1-R=1-{\left|\frac{{Z-Z}_{air}}{{Z-Z}_{air}}\right|}^2=1-{\left|\frac{Z_r-1}{Z_r+1}\right|}^2$$

(8)

$$Z_r=\sqrt{\frac{{(1+S_{11})}^2-S_{21}^2}{{(1-S_{11})}^2-S_{21}^2}} \left(9\right)$$

(9)

where, $Z_{air}$is the impedance of the free space and $Z_{ }$is the impedance of the metasurface. The absorption of the metasurface occurs when the effective impedance of the multiband absorber aligns with that of the free space. As depicted in Figure 7a,b, the real part approaches unity and the imaginary part approaches zero. The is calculated by employing Equations (8) and (9) for the conductivity of the ${\mathrm{V}\mathrm{O}}_2$ is ${4x10}^{3}s/m$ and ${2x{10}^5}^{ }s/m$, respectively.

The necessity for tunability and angular stability is paramount across numerous prospective applications. To ascertain the tunability, an examination of absorption across various conductivities of ${\mathrm{V}\mathrm{O}}_2$ is conducted, with the findings consolidated in Figure 8. As the conductivity of the ${\mathrm{V}\mathrm{O}}_2$ increases from ${10}^{ }s/m$ to ${2x{10}^5}^{ }s/m$ , there is a corresponding rise in plasmonic absorption within the ${\mathrm{V}\mathrm{O}}_2$ layer, resulting to an enhancement in the absorption of incident light. However, at extremely high conductivities, plasmonic absorption declines due to the heavy damping of plasmons, weakening their coupling with incident light and decreasing the absorption rate. At .11THz absorption reaches to 86% with the conductivity of ${\mathrm{V}\mathrm{O}}_2$ is ${4x{10}^3}^{ }s/m$ and 61% with the conductivity of ${\mathrm{V}\mathrm{O}}_2$ is ${2x{10}^5}^{ }s/m$ (conductive state), dropping to 0% at ${10}^{ }s/m$ (insulating state), demonstrating adjustable multiband absorption from 0% to 86% as shown in Figure 8a.

Furthermore, the impact of oblique incident angles is examined when the conductivity of ${\mathrm{V}\mathrm{O}}_2$ is ${4x{10}^3}^{ }s/m$ and ${2x{10}^5}^{ }s/m$ respectively, with the result shown in Figure 8b,c. The investigation reveals the proposed structure robustness resilience across a wide incidence angle spectrum (0^o to 85^o). Meanwhile, the absorption performance is investigated for LP and CP incident waves at the conductivity of ${\mathrm{V}\mathrm{O}}_2$ is ${4x{10}^3}^{ }s/m$ and ${2x{10}^5}^{ }s/m$ respectively, which demonstrates polarization independence when the incident terahertz wave is normally incident in wide incidence angles spectrum (0^o to 85^o), owing to the unit cell’s unique rotational symmetry structure, as shown in Figure 8d,e. The stability of the large incident angles and wide polarization angles stems from the unique configuration of the unit cell. In resulting, the proposed SFL Metasurface design for the absorber has excellent characteristics of wide-angle incidences and wide polarization incidences, which are significant for practical applications.

5. Conclusions

In summary, this research article represents multiple functionalities of THz metasurface based on vanadium dioxide (${\mathrm{V}\mathrm{O}}_2$), which has the functionalities of multiband absorption and multiband polarization conversion. When ${\mathrm{V}\mathrm{O}}_2$ is an insulating phase, the designed SFLM acts as a multiband polarization converter within the frequencies range of 0.11 THz – 0.12 THz and 0.26 THz – 0.27 THz; the designed metasurface exhibits robust capability for converting linear and circular polarizations to their respective cross-polarizations When the ${\mathrm{V}\mathrm{O}}_2$ transitions into a conducting state, the proposed metasurface operates as a multiband absorber, demonstrating absorption peaks at three distinct frequencies. Specifically, absorption rates exceeding 86% at 0.11 THz, over 60% at 0.288 THz, and 61% at 0.343 THz are attainable. Notably, the investigation emphasizes the stability of polarization conversion and absorption spectrum across varying polarization angles while sustaining remarkable multiband absorption capabilities even under high incidence angles. This multifunctional metasurface thus presents promising avenues for advancing future terahertz (THz) devices.