Transmissive terahertz metasurfaces with vanadium dioxide split-rings and grids for switchable asymmetric polarization manipulation

Metasurfaces containing arrays of thermally tunable metal-free (double-)split-ring meta-atoms and metal-free grids made of vanadium dioxide (VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2), a phase-change material can deliver switching between (1) polarization manipulation in transmission mode as well as related asymmetric transmission and (2) other functionalities in the terahertz regime, especially when operation in the transmission mode is needed to be conserved for both phases of VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2. As the meta-atom arrays function as arrays of metallic subwavelength resonators for the metallic phase of VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2, but as transmissive phase screens for the insulator phase of VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2, numerical simulations of double- and triple-array metasurfaces strongly indicate extreme scenarios of functionality switching also when the resulting structure comprises only VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 meta-atoms and VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 grids. More switching scenarios are achievable when only one meta-atom array or one grid is made of VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 components. They are enabled by the efficient coupling of the geometrically identical resonator arrays/grids that are made of the materials that strongly differ in terms of conductivity, i.e. Cu and VO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_2$$\end{document}2 in the metallic phase.

Metasurface research has been flourishing during the past decade, as quasiplanar designs enabling manipulation of the polarization state and wavefront are being investigated 1 . In particular, tunable metasurfaces have attracted a lot of attention in recent years. Various materials and tunability mechanisms have been proposed for use in metasurfaces designed to operate as spectral filters, polarization convertors, wave-front manipulating devices, etc. 4,5 . Transition from an insulator (I) phase to a metallic (M) phase, or vice versa, under the influence of a control parameter is a general feature of many natural materials 2,3 that allows control of the permittivity tensor in the complex plane at a fixed frequency. It has particular significance for metasurfaces with switchable functionality.
Vanadium dioxide (VO 2 ) is an attractive phase-change material displaying a hysteretic insulator-metal-insulator transition [6][7][8] . The underlying monoclinic-to-tetragonal phase change occurs as the temperature T is raised from a value slightly lower than 58 • C to a value slightly above 72 • C , the transition being reversible on cooling [9][10][11] . The values of the (complex) relative permittivity of VO 2 at any value of the free-space wavelength are significantly different in its two crystallographic phases. As VO 2 is a dissipative insulator (or semiconductor) (I phase) when monoclinic but metallic (M phase) when tetragonal provided that the free-space wavelength 1300 nm, it can be used in reconfigurable metasurfaces in the infrared 12,13 and terahertz 14,15 spectral regimes. The capability in switching of functionality in structures comprising tunable materials is often considered from the multifunctionality perspective [16][17][18] . Besides metasurfaces and metamaterials, applications such as smart windows 19 should be mentioned. Engineering of phase transition of VO 2 promises many new applications 20,21 .
The functionality of a structure comprising tunable materials such as VO 2 depends on the manner of their deployment in that structure. Commonly, a metasurface is an array of meta-atoms fabricated atop a slab that may be stratified in the thickness direction but is transversely homogeneous. The transverse dimensions of meta-atoms , geometrical features of VO 2 meta-atoms and grids proposed here, and specifics of their design that enable transmission and polarization manipulation properties as well as AT [57][58][59] , which are changed in a desired manner when VO 2 thermally transits from the M/I to I/M phase. A particular goal is to investigate the effects of coupling of meta-atom arrays or meta-atom grids that are made of different materials (i.e., VO 2 and a metal such as Cu), and their use for functionality switching.
In line with the aims of this study, the meta-atom is chosen to be much thinner than its transverse dimensions, so that resonances appear when VO 2 is in the M phase but not in the I phase. Thus, a VO 2 metasurface is simply a phase screen 56 when VO 2 is in the I phase, and the same is true for a VO 2 grid of proper thickness. This feature differs from our previous work 60 , in which the height of thermally tunable microrods is large enough to enable resonances for the I phase while a phase-screen regime corresponds to an intermediate phase; however, resonances in the I phase are undesirable for our present purpose.
Our preference here is that the metasurfaces operate in the transmission mode while thermal switching of functionality may include polarization manipulation and AT for one of the VO 2 phases. However, the scenarios of functionality switching that change the transmission mode into the reflection mode, or vice versa, are also studied. The proposed meta-atoms differ from the hybrid meta-atoms because they never combine VO 2 and a metal in one meta-atom. The concept is introduced and numerically validated for metasurfaces representing cascades of meta-atom arrays, but we emphasize that design optimization is beyond the scope of this paper. In contrast with the earlier works on polarization manipulation with thin VO 2 meta-atoms 44,45,47 , we consider polarization manipulation only in the transmission mode. Also, in contrast with Ref. 32 , where some of AT regimes were studied, we achieve AT without hybrid VO 2 -metal meta-atoms. In contrast with Ref. 27 , VO 2 grids in the proposed designs directly contribute to polarization manipulation and AT. The structures belonging to the two large classes of the polarization-converting metasurfaces capable of AT 59,61 are examined. Notably, although a single-array metasurface may produce cross-polarized components, at least two coupled arrays are needed to obtain well pronounced AT for linearly polarized (LP) waves.

Generic metasurface with thin and arbitrary shaped meta-atoms
A generic single-array metasurface comprising VO 2 meta-atoms is illustrated in Fig. 1. It is assumed that the meta-atoms may have rather arbitrary geometry, except that the thickness (here, along the z axis) of the VO 2 meta-atom should be chosen so that the VO 2 array can function as a phase screen for the I phase. For the M phase, the achievable functionality is determined by the geometrical and material properties of the metaatoms. They are arranged in a biperiodic array in the xy plane and are illuminated normally by a LP plane wave. CST Studio Suite 62 was used to calculate the transmission coefficients τ Q nm and the reflection coefficients r Q nm , where m ∈ x, y identifies the direction of the electric field of the incident plane wave, n ∈ x, y identifies the direction of the far-zone transmitted electric field, and the superscript Q = I when VO 2 is in the I phase but www.nature.com/scientificreports/ Q = M when VO 2 is in the M phase. The Cartesian unit vectors are denoted by û x , û y , and û z . The coefficients τ Q nm for front-face and back-face illumination, i.e., for propagation along the −z and +z directions, are additionally superscripted by → and ← . Since all materials used are Lorentz reciprocal 63 , |τ → yy | = |τ ← yy | , |τ → xx | = |τ ← xx | , |τ → xy | = |τ ← yx | , and |τ → yx | = |τ ← xy | . Additionally, some kind of geometric symmetry is needed to obtain |τ → xx | = |τ → yy | and |τ ← xx | = |τ ← yy | . More details can be found in Ref. 57 . The strongest asymmetry in transmission is expected to be achieved when, say, a forward-propagating incident wave of a specific LP state is fully converted to a transmitted wave of the orthogonal LP state, and the backward-propagating wave of the same LP state is completely reflected; i.e., either or The relative permittivity of VO 2 is given by ε VO 2 the angular frequency and f is the frequency, ε ∞ = 9 , ω VO 2 p = (σ VO 2 /ε 0τ VO 2 ) 1/2 is the plasma angular frequency, ε 0 is the permittivity of free space, τ VO 2 = 2.27 fs is the relaxation time, and the conductivity σ VO 2 = 40 S m −1 for the I phase ( T = 300 K) and 5 × 10 5 S m −1 for the M phase ( T = 400 K) 64 . For all data reported here, the VO 2 thickness h was chosen to be realizable with the use of available fabrication techniques 34,35 . In order to enable the phase-screen regime at 1.1 THz with VO 2 in the I phase, h = 8.75 µm so that (2πh/ ) ε VO 2 r 1/2 = 0.55 + 0.022i . The conductivity of Cu was fixed as σ Cu = 5.96 × 10 7 S m −1 .

Results and discussion
Switchable metasurfaces with one array of double split rings. Let us begin with the simplest illustration of switching that exploits the nearly complete electromagnetic disappearance of VO 2 in the I phase to enable switching for spectral filtering 65 . In Fig. 2, the meta-atom of metasurface 1A comprises two concentric rings of VO 2 , each with a small air gap g aligned parallel to û y , on top of a teflon substrate of relative permittivity equal to 2.1 and thickness denoted by d. Whereas the meta-atom is of transverse dimensions a × a , the thickness of the double split-ring (DSR) structure along the z axis is denoted by h. The inner split ring has mean radius r 1 and width w. The outer split ring has outer diameter b < a , width w, and mean radius r 2 = (b − w)/2. Representative spectra of the magnitudes of the co-polarized transmission coefficients |τ M xx | , |τ M yy | , |τ I xx | , and |τ I yy | are presented in Fig. 2b,c. Whereas frequency-selective transmission occurs when VO 2 is in the M phase, nearly frequency-independent and high transmission occurs for f ∈ [0.45, 1.35] THz when VO 2 is in the I phase. In the latter case, ε VO 2 r = 8.99 + 0.72i at f = 1 THz so that VO 2 cannot be said to be in the vacuum state (i.e., Re ε r = 1 on contrary to Ref. 60 ), but metasurface 1A functions as a phase screen with a relatively slight attenuation.
When VO 2 is in the M phase, transmission through the metasurface 1A depends on the polarization state of the incident plane wave. As can be observed in Fig. 2b, the minimums and maximums of |τ mm | have different spectral locations for E inc �û x (i.e., m = x ) than for E inc �û y (i.e., m = y ). These spectral locations can be shifted by changing h. Whereas max |τ M xx | > 0.8 but min |τ M yy | < 0.2 in some spectral regimes, max |τ M yy | > 0.8 but min |τ M xx | < 0.2 in other spectral regimes. In contrast, |τ I xx | ≃ |τ I yy | > 0.87 for f ∈ [0.45, 1.35] THz. Phase change may lead to a strong difference in transmission strength ( |τ I | 2 /|τ M | 2 > 50 ) in some spectral regimes, but a weak difference ( |τ I | 2 /|τ M | 2 ≤ 1.2 ) in others.
Rotation of a single-array metasurface, such as metasurface 1A, is a simple way to obtain significant crosspolarized transmission when the incident plane wave is linearly polarized 66 . To demonstrate the capability of VO 2 DSR structures for cross-polarized transmission, we rotated metasurface 1A about the z axis by 45 • to www.nature.com/scientificreports/ obtain metasurface 1B. Figure 2c presents spectra of the transmission coefficients |τ xy | = |τ yx | of metasurface 1B. The results are presented for the both M and I phases of VO 2 . Switching between nearly zero (OFF state, I phase) and moderately efficient (ON state, M phase) cross-polarized transmission is achieved. Manipulation of LP state by a metasurface has been often considered in connection with AT 57,59 , when transmission coefficients depend on which one of the faces of the metasurface is being illuminated. The single-array metasurfaces such as 1B are not capable of exhibiting AT in the case of LP incidence, regardless of the crystallographic phase of VO 2 .
Switchable metasurfaces with coupled arrays of double split rings. To demonstrate polarization manipulation and, in particular, asymmetry of the cross-polarized components (i.e., |τ M→ xy | � = |τ M← xy | ) due to metasurfaces comprising metal-free VO 2 meta-atoms and metal-free VO 2 grids, we first consider the ones with both front and back faces having DSRs in their unit cells; see Fig. 3. Each meta-atom in the double-array metasurface 2A has a VO 2 DSR on the front face and its rotated version on the back face, with a teflon ( ε r = 2.1 ) spacer between the front and back DSRs. Geometrically similar metasurfaces but with metallic resonators have been previously investigated for diodelike AT 57,58 . For metasurface 2A, cross-polarized transmission is expected to appear due to the coupling of two metal-free arrays in the M-phase of VO 2 .
Representative spectra of |τ Q→ xy | = |τ Q← yx | and |τ Q→ yx | = |τ Q← xy | are shown in Fig. 3d. Here, we obtain www.nature.com/scientificreports/  with (a,b); dark green and yellow colors stand for VO 2 and Cu components, respectively, and blue color stands for the teflon spacers. www.nature.com/scientificreports/ while |τ M→ xx | = |τ M→ yy | = 0.44 and 0.3, respectively. Hence, there is quite a good coupling between the M-phase VO 2 DSRs, although metallic components are not used. Cross-polarized transmission does not occur in the I phase of VO 2 , because the VO 2 arrays function as phase screens. Thus, cross-polarized transmission can be switched ON/OFF, i.e., the VO 2 phase change delivers a change between (i) a spectrally selective transmission regime with cross-polarized components and (ii) highly efficient and frequency-independent transmission without such components.
Intuitively, the efficiency can be increased by decreasing the overall losses and/or increasing the overall conduction capability. This is demonstrated by metasurface 2B, which differs from 2A only in that the DSRs on the back face are made of Cu. Therefore, in contrast with many polarization-manipulating metasurfaces 57-59 , metasurface 2B contains two dissimilar arrays of DSRs. Despite the significant difference between the DSR arrays in terms of conductivity, they are well coupled and deliver moderately high cross-polarized transmission, when VO 2 is in the M phase. The transmission spectra presented in Fig. 3e  , the ratio of the overall forward-to-backward transmission intensities is about 4.3 when E inc �û x and 3.7 when E inc �û y . This indicates a relatively low contrast for AT, so that a different design is needed to suppress co-polarized transmission, per conditions set by Eqs. (1) and (2). In terms of functionality switching, metasurface 2B is similar to 2A, because the transmission mode is conserved for both the M and I phases. However, transmission in the I phase now depends on the frequency.
To ensure suppression of both co-polarized and one of the two cross-polarized components, in line with Eqs. (1) and (2), we designed the triple-array metasurface 2C to have VO 2 DSRs on the front face and rotated VO 2 DSRs on the back face, with a VO 2 grid inserted in the middle of the spacer, i.e., similarly to the Cu grid in Ref. 59 . Representative transmission spectra are shown in Fig. 3f. Now we have |τ M→ xy | ≫ |τ M→ yx | and |τ M→ xy | ≫ max(|τ M→ xx |, |τ M→ yy |) , as desired. Metasurface 2C works as an AT device in the M phase and somewhat as a phase screen in the I phase. However, the peak values of |τ M→ xy | are relatively low, because of high ohmic losses and quite strong reflection.
In order to enhance cross-polarized transmission, we replaced the VO 2 grid in metasurface 2C by a geometrically identical Cu grid. The transmission spectra of the resulting metasurface 2D are presented in Fig. 3g. Now, we have |τ M→ xy | = 0.21 instead of 0.16 at f = 0.72 THz and |τ M→ xy | = 0.27 instead of 0.185 near f = 1.14 THz. Metasurface 2D exhibits AT when VO 2 is in the M phase but functions as an almost perfect reflector when VO 2 is in the I phase. Hence, in terms of functionality, it is not identical to metasurface 2C. As an example of AT efficiency in the M phase, |τ M→ | 2 /|τ M← | 2 ≈ 270 and 28 at f = 0.72 THz and 1.14 THz, respectively.
The triple-array metasurface 2E was designed to further enhance AT efficiency while preserving the capability of ON/OFF switching. It has VO 2 DSRs on the front face, rotated Cu DSRs on the back face, and a Cu grid at the middle. Transmission spectra are presented in Fig. 3h. Compared to metasurface 2D, |τ M→ xy | is now as twice as large. An efficient coupling of two different arrays of DSRs in metasurface 2E for VO 2 in the M phase is evident. For VO 2 in the I phase, metasurface 2E shows high reflection for a larger part of the considered frequency range but has windows of small-to-moderate co-polarized transmission. Figure 3i presents some results for metasurface 2F, in which the grid is made of VO 2 whereas all DSRs are made of Cu. When VO 2 is in the M phase, AT is strong with the unwanted transmission components suppressed. In particular, |τ M→ xy | = 0.465 at f = 0.74 THz and |τ M→ xy | = 0.44 at f = 1.179 THz, but |τ M→ xx | ≤ 0.06 and |τ M yx →| ≤ 0.01 . When VO 2 is in the I phase, the grid functions as a phase screen and both |τ I→ xx | and |τ I→ xy | are significant; e.g., |τ I→ xx | = |τ I→ yy | = 0.37 and |τ I→ xy | = 0.52 at f = 0.65 THz while |τ I→ yx | ≤ 0.05 . To compare, |τ I→ xx | = |τ I→ yy | = 0.22 and |τ I→ xy | = 0.52 at f = 1.065 THz, whereas |τ I→ | 2 /|τ I← | 2 = 5.58 (the case of a relatively weak asymmetry in transmisson). Roughly speaking, whereas metasurface 2F for the M phase is similar to 2E for the M phase, 2F for the I phase is similar to 2B for the M phase. In terms of functionality, the phase change of VO 2 results here in switching between AT for the M phase and spectrally selective transmission with significant co-and cross-polarized components for the I phase.
The obtained results show that high-contrast switchable AT can be achieved (at least with moderate efficiency) in the structures comprising two coupled VO 2 DSR arrays, whereas transmission mode is conserved for both the M and I phases of VO 2 . Using only one VO 2 array and coupling of VO 2 and Cu arrays, we can achieve higher transmittance in the AT regime when VO 2 is the M phase, but at the price of abandoning the transmission mode for the I phase. It is worth noting that using two different arrays in one structure (such as for the metasurfaces 2B and 2E) generally leads to |τ M→ xx | � = |τ M→ yy | and |τ M← xx | � = |τ M← yy | . However, for the proposed metasurfaces, selection of one or the other LP state for the incident wave does not lead to the principal difference in the functionality-switching scenarios.
Switchable few-layer metasurfaces with one split-ring array. Next, we present a set of few-array metasurfaces designed for switchable AT enabled by polarization conversion. In contrast with "Switchable metasurfaces with coupled arrays of double split rings", the metasurfaces studied here comprise only one array of subwavelength split rings. Notably, few-array non-switchable metasurfaces have been studied earlier that are either capable or incapable of polarization manipulation and AT 61 www.nature.com/scientificreports/ broadband polarization conversion and related AT in a switching scenario. The polarization-converting (but non-switchable) structures from Refs. 61,73 were taken as the pre-prototypes. Each unit cell of the triple-array metasurfaces 3A-3D shown in Fig. 4 has a front grid of strips parallel to the y axis, a single split ring (SSR) in the middle with a gap aligned parallel at 45 • to both the x and y axes, and a back grid of strips parallel to the x axis. Any two consecutive arrays are separated by a silica ( ε r = 2.25 ) spacer. The strips on the front face are made of VO 2 (3A, 3C, 3D) or Cu (3B), the SSR is made of VO 2 (3B, 3C) or Cu (3A, 3D), and the strips on the back face are made of VO 2 (3A, 3C) or Cu (3B, 3D). Transmission spectra are presented in Fig. 4, first of all, to demonstrate the capability of exhibiting AT.
As shown in Fig. 4b for metasurface 3A, the M phase of VO 2 delivers strong polarization conversion and AT, since all three arrays work as metallic arrays. Both grids operate as highly transmissive phase screens when VO 2 is in the I phase, so that the Cu SSR array by itself enables the cross-polarized components but without AT. For f ∈ [0.8, 1] THz, switching between the M and I phases of VO 2 leads to the functionality switching between Interestingly, at f = 1.3 THz, |τ → yx | is weakly sensitive to a phase change, but |τ → xy | is strongly sensitive. Comparable co-and cross-polarized components occur for the I phase at 1.15 THz. Thus, diverse scenarios can be achieved with metasurface 3A functioning in the transmission mode for both phases of VO 2 .
VO 2 and Cu are interchanged in metasurface 3B. Figure 4c shows that strong polarization conversion and AT are achieved when VO 2 is in the M phase, so that |τ M→ yx | ≫ |τ M← yx | for the entire frequency range, and |τ M→ yx | ≫ |τ M→ xx | for most of it. The VO 2 SSR array functions as a transmissive phase screen in the I phase, i.e., without AT. In terms of functionality, we obtain switching between transmission-mode polarization conversion and AT for the M phase and reflection in the I phase. This metasurface is well suited to function as an ON/OFF switchable diodelike AT device, but the transmission mode is not conserved while switching.
Both strip grids and the SSR arrays in metasurface 3C function as highly transmissive phase screens when VO 2 is in the I phase, but as metal arrays when it is in the M phase. As a result, we obtain strong polarization conversion and related AT for the M phase and the almost perfect and f-independent transmission without polarization conversion for the I phase for 0.65 < f < 1.05 THz in Fig. 4d. Thus, broadband ON/OFF switching of AT can be obtained in the transmission mode. For example, when VO 2 is in the M phase, |τ M→ | 2 /|τ M← | 2 = 2.85 × 10 3 and |τ M→ yx | = 0.64 at f = 0.88 THz. Only the strip grid on the front face is made of VO 2 in metasurface 3D, the other two arrays being made of Cu. Similarly to Fig. 3g, there are two geometrically identical arrays (strip grids, here), which are rotated with respect to each other and made of different materials. The VO 2 strip grid is a transmissive phase screen when VO 2 is in the I phase, but a metallic grid otherwise. In Fig. 4e THz. At f = 1.24 THz and 1.32 THz, |τ → yx | is insensitive to the crystallographic phase of VO 2 . The obtained switching scenarios strongly depend on the polarization state of the incident wave. When E inc �û x , switching between two AT regimes with strong cross-polarized and weak co-polarized transmission components can be achieved by thermally changing the crystallographic phase of VO 2 . When E inc �û y , switching between strong co-polarized transmission for the I phase of VO 2 and strong reflection for the M phase of VO 2 is obtained. However, in the contrast with metasurfaces 3A and 3C, nearly perfect frequency-and polarization independent co-polarized transmission cannot be achieved for 3D within the used frequency range.
(1) and (2). Since the switching of AT band is enabled by one of the cross-polarized transmission components, it is quantified here in terms of the difference (in dB) between |τ M→ yx | and |τ I→ yx |. A remarkable result observed in Fig. 4 is that the highly efficient, switchable, broadband polarization conversion and AT can be obtained in the transmission mode for both the M and I phases of VO 2 . In contrast with Refs. [29][30][31] , placing VO 2 pads/inserts into subwavelength metallic resonators is unnecessary. Moreover, the SSRs can be made of Cu, while only the grids are made of VO 2 , to enable switching while operating in the transmission mode. One strip grid can be made of VO 2 in order to obtain a specific switching between two directionally selective regimes by thermally changing the crystallographic phase of VO 2 . Note that the VO 2 grids in M phase directly contribute to polarization manipulation for metasurfaces 3A, 3C, and 3D, so that their role is not restricted to tunability.

Discussion
To summarize, we have proposed and numerically validated that few-layer metasurfaces comprising electrically thin metal-free meta-atoms and grids made of VO 2 are capable of polarization manipulation and asymmetric transmission. The advantages of using a thermally tunable material, such VO 2 , include that: (a) the temperaturerange required for extreme ON/OFF switching is relatively narrow and (b) no biasing circuit is needed. Only a simple heater is needed on-site. We have demonstrated that the thermal transition from the M/I to the I/M phase of VO 2 may result in extreme functional reconfiguration, while either conserving the transmission mode of operation or changing it for the reflection mode.
A key feature is that the arrays of meta-atoms and/or grids made of VO 2 function as highly transmissive phase screens for VO 2 in the I phase. Metasurfaces comprising designed meta-atoms and grids are capable in diverse scenarios of switchable polarization conversion and AT. In particular, triple-array metasurfaces comprising one SSR array and two strip grids, or two DSR arrays and one grid, may enable switchable polarization conversion and AT, provided that VO 2 components are properly incorporated in the device. To operate in the transmission mode for both phases of VO 2 and have AT for one of the two phases, (a) all SSR/DSR arrays and grids, or (b) the central grid, or (c) grids on both faces should be made of VO 2 . The two DSR arrays can be well coupled even if one of them is made of VO 2 and the other of a metal, i.e., when the difference in DSR conductivities is significant; indeed, the smaller conductivity of VO 2 in the M phase, as compared to that of the metal, does not lead to the disappearance of switchable polarization conversion and AT. The same is true for triple-array metasurfaces with one face grid made of VO 2 and the other face grid made of Cu. Significant dependence of the co-polarized transmission on the polarization state of the incident wave may enable an additional degree of freedom in functionality switching. The proposed concept and sample designs constitute a perfect platform to design new switchable, multifunctional, polarization-manipulating, and diodelike THz devices.
Although considerable theoretical literature on metasurfaces containing VO 2 exists, idealized constitutive parameters are used therein. As stated earlier, it is challenging to deposit pure VO 2 because it is difficult to avoid the formation of polycrystalline vanadium oxide (i.e., VO ζ ) 34,35 . The value of ζ depends on the deposition process and the process parameters. Most noticeably, the thermal hysteretic insulation-metal-insulator transition is affected by polycrystallinity 74 . As our proposed designs avoid the temperature-regime of hysteresis, we expect that their performance is going to be affected only to a limited degree by |ζ − 2| ≤ 0.1 . Therefore, with stochiometry engineering 8 , vanadium-based materials 74,75 have a high potential in novel physical scenarios and related applications, which still need to be explored.

Methods
An intuitive design approach, based on the qualitative analysis of the dominant physics and simple estimates, was applied. CST Studio Suite 62 , a commercial software was used for numerical simulations. It is based on the finite integration method with controllable convergence and accuracy, and is particularly appropriate for unitary metasurfaces comprising unit cells of complex geometry. The frequency-domain solver and unit-cell (Floquet-Bloch) boundary conditions were adopted along with a tetrahedral mesh for simulations. www.nature.com/scientificreports/