Ultrafast reprogrammable multifunctional vanadium-dioxide-assisted metasurface for dynamic THz wavefront engineering

In this paper, for the first time, a new generation of ultrafast reprogrammable multi-mission bias encoded metasurface is proposed for dynamic terahertz wavefront engineering by employing VO2 reversible and fast monoclinic to tetragonal phase transition. The multi-functionality of our designed VO2 based coding metasurface (VBCM) was guaranteed by elaborately designed meta-atom comprising three-patterned VO2 thin films whose operational statuses can be dynamically tuned among four states of “00”-“11” by merely changing the biasing voltage controlled by an external Field-programmable gate array platform. Capitalizing on such meta-atom design and by driving VBCM with different spiral-like and spiral-parabola-like coding sequences, single vortex beam and focused vortex beam with interchangeable orbital angular momentum modes were satisfactorily generated respectively. Additionally, by adopting superposition theorem and convolution operation, symmetric/asymmetric multiple beams and arbitrarily-oriented multiple vortex beams in pre-demined directions with different topological charges are realized. Several illustrative examples successfully have clarified that the proposed VBCM is a promising candidate for solving crucial terahertz challenges such as high data rate wireless communication where ultrafast switching between several missions is required.

simultaneous manipulations of EM waves in both space and frequency domains which substantially expand the application scope of digital metamaterials [33][34][35][36] . Until the present, some papers have brought forward electronically tunable coding metasurfaces in microwave frequency by using pin diodes in each coding element 24,37,38 . For example, Huang et al. presented a method to design a 2-bit digitally-controlled coding metasurface in order to realize different far-field patterns. By integrating two pin diodes in each coding particles and switching their operating states, producing four phases responses was realized 39 . In light of the complexity and lack of commercial access of active elements (pin or varactor diodes) at high frequencies, scaling these devices to THz regime is very difficult if not impossible 40 .
For real-time manipulating THz waves, an efficient technique is to integrate standard metasurface with phase-change materials (PCM), for instance, graphene [41][42][43] , liquid crystal [44][45][46][47] , GST 48,49 and vanadium dioxide (VO2) [50][51][52][53] . To implement such a platform, we have benefited from VO2 exotic properties. VO2 is a smart material, that undergoes an ultrafast and brutal reversible first-order phase transition from insulating monoclinic (P c 2 / 1 ) to the metallic tetragonal (P mnm 4 / 2 ) phase, above critical temperature =°T 68 C c 53 . Since the critical temperature is a function of V-V distance in crystal structure, decreasing in critical temperature would be feasible by adopting doping technique 54,55 . This metal-insulator transition (MIT) can be provoked by thermal 56,57 , optical [58][59][60] or electrical (charge injection or Joule heating) 61,62 stimuli. The origin of such an insulator to metal transition is still under debate [63][64][65] . This transition in VO2 can occur within an order of several nanoseconds or even in picoseconds range for optical activation 66 . The VO2 phase alternation that was studied in ref. 67 and ref. 68 through time-resolved x-ray diffraction and time-resolved optical transmission respectively, revealed that only a short time (<500 fs) was required for taking place the MIT. Moreover, Cavalleri's study indicated that the photo-induced transition time for VO2 thin film (50 nm) could reach 80 fs 69 .
Electrical and optical property of VO2 has dramatic changes (4 to 5 order of magnitude change of the electrical conductivity) across the two phases 70 . Owing to the ultrafast transition duration, almost near room critical temperature and fruitful structural transition behaviors, VO2 has became a striking material in tunable metamaterial devices at GHz 71,72 , optics [73][74][75] and has numerous fantastic applications in THz frequencies such as reconfigurable THz filters 76,77 , polarization converter [78][79][80][81][82] , reflection/transmission THz waves modulator [83][84][85] , tunable THz absorbers 86,87 and reconfigurable antennas 88,89 . In 2016, a simple VO2-assisted digital metasurface was proposed to dynamically control the near-infrared light 90 . By allowing voltage to be locally applied to the VO2 and distributing unit cells in one direction, switchable beam splitters with only limited splitting angles have been realized. Regarding that the reflection phase of the proposed structure is either 0° or 180° relying on the state of VO2, several basic functionalities and fundamental THz challenges such as anomalous reflection, multi-beams and vortex beam generation with adjustable properties are not accessible with this structure.
In departure from the abovementioned work, wherein the mission of utilizing VO2 in the structure is limited to only a single tunable function, we present, to the best of author's knowledge, the first VO2-based coding metasurface (VBCM) that can be reprogrammed for realizing multi-type functionalities from vortex beams with different topological charges toward emitting multiple arbitrarily-oriented pencil beams. Although, some articles have reported graphene-based coding metasurface for real-time manipulation of THz waves recently, but the graphene is of great loss compared to VO2. Graphene consists only of a monolayer of carbon atoms compared to volumetric VO2 and its switching time is thousands of times slower than VO2. Disadvantages such as difficult fabrication process and very low switching time, made graphene to be impractical in real-time wave manipulation meta-devices 90,91 .
Although today's applications in THz regime as high data rate wireless communication and ultra-massive MIMO communication are in dire need of ultrafast THz wavefront manipulation, we are not aware of any reports of ultrafast real-time THz wavefront engineering and this field is still largely unexplored and we believe that our proposed ultrafast VBCM structure has the great potentials to fill this gap.
In this paper, an ultrafast versatile 2-bit VO2-based coding metasurface with programmable meta atoms is designed whose operational statuses can be dynamically switched between four states of "00", "01", "10", and "11". By applying external bias voltage controlled by an FPGA platform to the elaborately designed meta-atom, the capability of manipulating THz wavefront in a real-time manner has been achieved. Each constitutive unit-cell of the VBCM integrates three-patterned VO2 layers fed by two biasing voltages (ON/OFF), therefore it can alterably possess four reflection phase responses of 0, π/2, π, and π 3 /2 without changing the geometrical parameters. In order to expose the capacities of our structure in multifunctional wavefront engineering, we have arranged various coding samples to accomplish different interchangeable functions from vortex beams generation with different OAM modes toward emitting multiple arbitrarily-oriented pencil beams without re-optimizing or re-fabricating the structure. By adopting convolution and superposition theorems on the far-field pattern, several THz fundamental challenges such as symmetric/asymmetric multiple pencil beams and arbitrarily-oriented multiple vortex beams with different OAM modes, with ultrafast switching time are solved by our proposed VBCM structure. The presented concepts have been confirmed through a number of numerical simulations and theoretical predictions. To the best of our insight, this is the first ultrafast re-programmable coding metasurface based on VO2 phase transition that armed a platform for realizing some rich interchangeable missions at THz frequencies.
The authors believe that the proposed VBCM paves the way for ultrafast multifunctional THz wavefront engineering and future practical applications.

Results
Design of the switchable multifunctional VBCM. Figure 1 represents the basic meta-atom of the VO2 integrated metasurface composing of three layers, which, from top to bottom, are the VO2 layers, dielectric substrate, and a gold plane to impede the transmission energy into the back of the VBCM. At the top layer, three patterned VO2 layers deposited on the sapphire as the suitable substrate (ε δ = . . = . 9 98 to reactive electron-beam evaporation 99 , but it is preferable to use sapphire substrate to reach high-quality VO2 thin films due to the beneficial lattice matching effect 100 .
The complex dielectric properties of the VO2 can be characterized by the Bruggeman effective-medium theory in the THz range, wherein, ε d and ε m denotes the dielectric constant of the semiconductor and metallic regions respectively and V represents the volume fraction of metallic regions 78 .
At room temperature the dielectric constant of VO2 is about 9 in the insulting state 79,82,101,102 and by applying external bias voltage directly to VO2 thin film, the structural transformation occurs and VO2 turns into the rutile phase. Typical VO2 films that are grown on c-type or r-type sapphire substrate, display electrical conductivity in the range of 10~100 S/m in insulating state and as high as an order of 10 5 S/m in the metallic state 76 . Frequency independent conductivity of VO2 is set to be σ = S m 10 / (OFF state) 92,103 and σ = × S m 5 10 / 5 (ON state) 62,80,104 in the insulating and metallic phase corresponding to T c = 300 K and T h = 400 K respectively. The periodicity of our subwavelength meta-atoms is μ = P m 100 . The other geometrical parameters are W 1 = 90 μm, W 2 = 62 μm, W 3 = 14 μm, a 1 = 39 μm, a 2 = 26 μm, b 1 = 90 μm, b 2 = 52 μm, b 3 = 19 μm, c 1 = 14 μm, c 2 = 26 μm respectively. The thickness of VO2 films is assumed to be μ = t m 1 . All of these geometrical dimensions are extracted from an extensive simulation and optimization process to obtain four reflection phases of 0, π/2, π, and π 3 /2 (phase step of 90°) to mimic four digital states of "00", "01", "10", "11" and any change in dimensions will worsen this phase differences. Contrary to geometrically encoded metasurface, the employed meta-atoms of the VBCM structure in this study have the same dimensions but with different properly bias voltages leading to call this structure as bias-encoded metasurface. All the numerical simulations are carried out by means of Commercial software CST Microwave Studio. For evaluating the reflection characteristics for the infinite array of digital elements, open (add space) boundary condition is applied along z-axis whilst periodic boundary conditions are utilized along the xand y-directions to incorporate the mutual coupling effect between neighboring elements. The simulated reflection spectra are depicted in Fig. 2(a,b) for different sets of external bias voltage. It is clearly observed from Fig. 2(b) that with appropriately biasing three patterned VO2 layers from top to bottom as follows: ("OFF/OFF/ OFF", "ON/ON/ON", "ON/OFF/ON" and "OFF/ON/OFF"), our elaborately designed meta-particle has been successful in providing a phase step of 90° at 0.44 THz. Also, the reflection amplitudes are above 0.82 which means that the meta-atom non-absorptive behavior has been achieved. We have deliberately selected these two temperature points since VO2 at these two points is in the dielectric or metallic steady-state phases 90 . Furthermore, the Ohmic loss of the meta-atom is maximum at intermediate temperatures which leads to a sharp drop in the reflection amplitude which makes device efficiency to be minimized. It is worth mentioning that to maintain the symmetry of the simulation results, both top and bottom VO2 layers (red colors in Fig. 1) must be geometrically identical and switched "ON" and "OFF" simultaneously. Finally considering the abovementioned constraints and avoiding VO2 intermediate phases, employing three-patterned VO2 layers in the meta-particle that leads to the complexity of the biasing system is forced to users. It can be found from Fig. 2(b) that the perfect 2-bit characteristic occurs at 0.44 THz which means that the approximate size of our subwavelength meta-atoms equals λ/7. Figure 1. Sketch representation of reprogrammable VO2 based coding metasurface controlled by an FPGA platform. Different spatial coding patterns can be encoded onto the structure simultaneously through computer-programmed biasing network at distinct moments of t 1 , t 2 , t 3 ….
Meanwhile, a maximum phase range of 260° is attained here that is wide enough for acceptable functioning of our ultrafast multifunctional VBCM for real-time THz wavefront engineering.
Before beginning the next section, it is worthwhile mentioning that the multi-functionality of the VBCM structure can be realized when external dc bias voltage/current is applied to each individual unit-cell independently through the Au electrodes deposited on the top of the patterned VO2 layers to form a good Ohmic contact. The main mechanism for this voltage/current driven MIT process remains controversial as it may originate from Joule heating 61,105 or electric field effects 106,107 . Recently the viewpoint of Zimmers et al. 61 , was that the local Joule heating plays a predominant role in the dc voltage (or DC current) induced MIT. Wu et al. 106 claim that the Joule heating effect was negligible and the electric field alone is sufficient to induce MIT. Beyond these two competing claims, some papers were commented that these two main mechanisms are considered to be mixed and were difficult to disentangle these two effects in the VO2 voltage-driven MIT process 108,109 .
To minimize the EM coupling between the adjacent meta-particles, the VBCM structure comprises of × N N array of × M M identical subwavelength unit-cells that constructs the so-called super-unit cell or lattices. The length of one lattice is MP and the length of the whole VBCM is equaled to NMP. The metasurface is assumed to be illuminated by a y-polarized normal incident plane wave throughout the paper. From the well-known antenna theory, the far-field scattering pattern function of the VBCM can be expressed by 11 : in the above equations, θ and ϕ are the elevation and azimuth angles of the desired direction, P demonstrates the period of lattice along both x-and y-directions, a mn and ϕ mn are the demonstrator of reflection amplitude and phase of each lattices respectively and π λ = k 2 / where λ is a working wavelength. Establishing a 2D inverse fast Fourier transform (2D-IFFT) to accelerate the calculation make this formula beneficial for the prediction of scattering patterns caused by different coding sequences. The VBCM structure occupied with × 8 8 array of × 6 6 identical unit-cells throughout this study. In the two next layouts, only digital elements of "00" and "10" with 180° phase difference are adopted. By distributing these two digital elements in alternate rows or columns, the reflected wave is split into two symmetrically oriented scattered beams (θ = 34°, ϕ = 0°, 180°) governed by generalized Snell's law 110 as shown in Fig. 2(d). By arranging the above two digital elements in a chessboard configuration (see Fig. 2(c)), the incoming energy redirected into four scattered main beams of (θ = 53°, ϕ = 45°, 135°, 225°, 315°) at 0.44 THz that have an excellent conformity with the theoretical predictions. According to generalized Snell's law, with increasing the observation frequency, the reflection beams come closer to the boresight direction 111 .
Obviously in anomalous reflection phenomena, by purposefully arranging the coding sequences, the reflected beam can be directed toward a pre-determined azimuth angle in each of four quadrants of (ϕ r = 45°, 135°, 225°, 315°) and at the same frequency, the generated beam elevation angle can be shifted by changing the dimension of lattices of the VBCM structure and beam steering functionality could be envisioned. To validate the concept, the VBCM employed with different size lattices of M = 3, 6, 7 are analyzed by the full-wave simulation (see Fig. 3(c-e)) and a good agreement between numerical simulations and theoretical predictions (see Fig. 3(f)) confirm the validity of the presented beam steering analysis in a real-time ultrafast VBCM structure. These examples clearly illustrate that the instantaneous access to an anomalous reflection with outgoing directions is simply realized with our designed VBCM structure. Eventually, by varying the biasing voltage controlled by an FPGA platform, a single beam that can be dynamically altered in pre-determined angles with ultrafast switching time is provided with our presented structure which has fascinating functionality to implement in various applications such as ultrafast THz wireless communications and tracking systems.

multi-beam generation.
Recently, it has been revealed that when two different coding patterns are added together through the superposition theorem 112 , a combined coding pattern will perform both functionalities simultaneously aid to reach a metasurface with several missions such as multi-beam generation. Two different layouts will be anticipated here to picture the multi-mission capability of the VBCM device. Let us consider two metasurfaces driven by gradient coding sequences along different directions to generate two single beams toward (θ r = 25°, ϕ r = 45°) and (θ r = 25°, ϕ r = 315°) respectively. We adopted the superposition theorem here to design a www.nature.com/scientificreports www.nature.com/scientificreports/ 2-bit VBCM that redirect the incident THz wavefront from the normal direction into two reflected beams with (θ r = 25°, ϕ r = 45°, 315°) driven by a coding sequence generated by: j j j in which ϕ 1 and ϕ 2 are arguments of the two primary complex codes and ϕ 0 represents the argument of superimposed complex code. The 2D phase map of two independent phase gradient metasurface and the final required phase profile of the mixed metasurface to generate proposed multi-beams is depicted in Fig. 4(a) and the 3D scattering pattern of superimposed VBCM is shown in Fig. 4(b). The next example is dedicated to adopting the superposition theorem to add two different gradient coding patterns with distinct elevation angles. Before we delve into the full-wave simulations, it should be noted that when multiple independent pencil beams with different elevation angles are added together, the superimposed coding metasurface generates multiple beams with asymmetric power ratio levels. More recently, by revisiting the addition theorem in the metasurface, our team provided a generalized version of the superposition theorem to estimate the exact amount of power ratio of the multiple beams 113 . Considering the cosine function as the element factor of the metasurface particles, the amount of power distribution of two asymmetric pencil beams can be estimated as follows: which means that in our designed VBCM structure with two asymmetrically oriented scattering beams, the arbitrarily oriented pencil beam with a higher elevation angle (θ) carrying lower power intensity. Let us consider, two metasurfaces with phase gradient coding sequences one of which makes the incident wave reflect at oblique angles of (θ r = 17°, ϕ r = 180°) while the other generates single anomalously scattered beam toward (θ r = 25°, ϕ r = 45°). Exploiting the same design approach results in a VBCM structure that has incorporated both aforesaid gradient codes satisfactorily generates two asymmetrically oriented reflected beams along the pre-determined directions as shown in Fig. 4(d). The initial phase distribution of the metasurfaces before additional operation and the final phase distribution of the mixed metasurface are also provided in Fig. 4(c). Inspired by the addition theorem, these examples proved the correct functioning of our proposed VBCM structure to create symmetric/ asymmetric multiple beams in pre-determined directions which is dynamically interchangeable.

Spiral coding Sequence: vortex beam generation. Single vortex wavefront carrying different OAM
modes. Since being discovered in 1992 114 , vortex beams carrying orbital angular momentum (OAM) has experienced increasing levels of attention for its potential opportunities in high speed communication 115 , fast imaging 116 , optical manipulation 116 , etc. Advantages such as improving the channel capacity without increasing the bandwidth and orthogonality of different topological charges 117 have prompted researchers to present various novel methods for generating vortex beams from spiral phase plates 118 to antenna arrays 119 . However, there are still critical bottlenecks that are waiting for the solutions so that only a few reports involve THz vortex beam  www.nature.com/scientificreports www.nature.com/scientificreports/ realization. Furthermore, future progress in classical and quantum systems requires rapid switch between different OAM modes and at present, they are suffering from lack of ultrafast reconfigurable THz device in order to switch between different OAM modes in a real time-manner.
In order to address the aforesaid restrictions and to further clarify the versatility of our proposed ultrafast tunable VBCM structure, four THz vortex wavefront with different OAM modes (l = 2, −2, +1, −1) are envisaged. In each case, the coding metasurface is encoded by a spiral coding sequences generated by: a rctan (8) In another word, after dividing coding metasurface into N equal segments with phase differences ϕ ∆ of the neighboring segments, the relation between OAM topological charges and the number of segments for our 2-bit coding metasurface yields to: Toward this aim, in order to generate OAM beams with = ± l 1 and = ± l 2, we divided the VBCM structure into four and eight segments with phase shift ranging from 0 to 2π and 0 to 4π respectively. The instant electric field intensities of these four configurations when the observation plane is set as 170 μm away from the center of VBCM structure with an area of μ μ × m m 4800 4800 , shown in Fig. 5(c-f). 3D far-field scattering patterns (see Fig. 5(a,b)), demonstrate that the VBCM can successfully generate vortex beams carrying OAM (l = 2, −2) at 0.44 THz. According to Fig. 5, a typical doughnut-like intensity profile with an amplitude null in the center (15 dB lower than the annular high-intensity region) satisfies the far-field feature of OAM beams. The capability of ultrafast switching between different topological charges in a real-time manner is provided by our VBCM structure that can empower dramatic advances in wideband OAM based multi-user system where the beams topological charges identify the routing.
Arbitrarily-oriented multiple vortex beams. As a fascinating application in THz OAM-based MIMO systems, generating oriented vortex beam at a predetermined direction which can be dynamically tuned can solve the crucial challenges of this research area. Convolution operation as a simple yet wonderful solution has armed a platform to produce obliquely directed OAM vortex wavefront. Regarding the Fourier transform relation between the coding patterns and its far-field scattering patterns, the scattering pattern shift functionality can be achieved by the means of convolution theorem 120 by adding spiral phase distribution with gradient coding sequences. We start with a simple gradient coding sequence of [00, 01, 10, 11 …/00, 01, 10, 11 …] which generates single anomalously scattered beam toward (θ r = 17°, ϕ r = 180°). In this example, the OAM-generating spiral-type coding pattern (l = 1) is multiplied by the above gradient coding sequence, yielding a tilted vortex scattering pattern along predetermined direction. (see Fig. 6(b,d)). As the next layout, a phase gradient coding sequences of M 1 = [00, 01, 10, 11 …/01, 10, 11, 00 …/10, 11, 00, 01 …/11, 00, 01, 10 …] has added to elaborately encoded four 2 ). The 2D reflection phases map and 3D far-field scattering pattern of such encoded metasurface (M 3 ) is depicted in Fig. 6(a,c). As can be deduced from Fig. 6, the rotation angle of generated arbitrarily-oriented OAM beam with l = 1, which is dictated by the phase gradient coding sequence has good conformity with our theoretical predictions based on Eqs. 4 and 5.
focused vortex beam. The versatility of our designed VBCM also has equipped a platform to focus the incident THz wavefront into a pre-determined point. The focal length (Z focal ) can be dynamically altered by suitably changed the biasing system of VBCM driven by a parabola phase distribution along the radial direction. Furthermore, our proposed structure provides the ability to obtain several focused vortex beams with ultrafast switching time between different topological charges and focal lengths. To engineer this feature of the work, the phase profile of such encoded VBCM must involve both spiral and parabolic phase distributions simultaneously which can be expressed by 121 : To validate the concept, four focused wavefronts carrying OAM with diverse topological charges and focal l e n g t h s o f l Z m ( , ) ( 1, 800 ) have been designed and exemplified in Fig. 8. In each case, the required 2D spiral-parabola phase map for realizing the corresponding digital states has been pictured in Fig. 8(e-h) respectively. Eventually, illuminating by a y-polarized normal incident plane wave, the focused-vortex-generating VBCM structures are build-up and the simulated normalized electric near-field intensities in the sampling planes at the corresponding focal lengths (see Fig. 8(a-d)) demonstrate the capability and flexibility of the designed VBCM to generate THz focused vortex beams with different OAM modes and focal lengths which is dynamically switchable. potential fabrication procedure of the designed structure. In this section, we furnish a brief presentation on the current fabrication technologies 20,122 for fabricating our proposed VO2-based meta atom structure. This process can follow the steps below (see Fig. 9): (a) photoresist is spin-coated and deposited on a substrate using photolithography to form the pattered VO2 bricks; (b) a desired sample pattern is transferred to photoresist by writing the pattern with an electron beam lithography; (c) 1 μm thick VO2 layer is then prepared on 38 μm c-type sapphire substrate using magnetron sputtering technique; (d) lift-off process is adopted in order to dissolve the photoresist and leave behind the film only in the patterned area. Remaining patterned VO2 bricks is annealed at ~450 °C; (e) finally gold layer is deposited on the backside of the substrate as a ground layer. Following the above steps, one can envision a practical fabrication within the scope of the current fabrication technologies. www.nature.com/scientificreports www.nature.com/scientificreports/

Discussion
The necessity of integrating multiple diversified functionalities in a single structure at THz frequencies sparks our curiosity to design a re-programmable multi-mission coding metasurface incorporated with phase-changes materials. In this paper, for the first time, we have proposed a novel reconfigurable metasurface to dynamically manipulate THz wavefront by utilizing the insulator to metal transition of VO2. VO2 thin film exhibits ultrafast switching time which is of great practical significance in ultrafast THz communication. The versatility of our proposed VBCM structure was successfully clarified with different illustrative examples from anomalous reflection to focused vortex beam generation. Furthermore, symmetric/asymmetric multiple beams along the pre-determined directions have also been realized by applying superposition operation. By encoding the VBCM by spiral-like and spiral-parabola-like coding sequences, single vortex beam carrying OAM and focused vortex beam were satisfactorily generated respectively. Additionally, by adopting convolution operation, arbitrarily-oriented multiple vortex beams in pre-determined directions with different topological charges in a real-time manner was realized with Figure 9. Potential fabrication procedure of the designed structure: (a) 38 μm thick c-cut sapphire is prepared; (b) Photoresist is spin coated and deposited on a substrate using photolithography; (c) Photoresist is patterned with an electron beam lithography; (d) 1 μm thick VO2 layer is deposited using magnetron sputtered technique; (e) lift-off process in order to dissolve the photoresist and final VBCM structure. Remaining patterned VO2 bricks is annealed at 450 °C. Note: VO2 is shown in two different colors (red and blue) just for better visualization. www.nature.com/scientificreports www.nature.com/scientificreports/ our elaborately designed VBCM. Our ultrafast reconfigurable THz device meets well the THz future industrial demands which require rapidly switch between different OAM modes. In the two last sections, we demonstrated the feasibility of fabricating the VBCM within the realm of the current fabrication technologies and also provided a method for biasing the structure. This simple yet fruitful structure holds great potential for dynamically THz wavefront engineering and can enable advanced applications such as high data rate wireless communication and ultra-massive MIMO communication.

Methods
We used Matlab software version 2016a (https://www.mathworks.com/products/matlab.html) to generate required phase on the surface. All the full-wave numerical simulations are carried out by the means of CST Microwave Studio version 2016 (https://www.3ds.com/products-services/simulia/products/cst-studio-suite/).
The biasing mechanism of the proposed VBCM structure can be realized by photo-induced phase transition for very fast switching time or employing resistive heater electrode for Joule heating actuation. For the sake of briefness, the latter method is only investigated in this section. A resistive heater electrode is placed at the perimeter of the layers to trigger the insulater to metal transition independently for each patterned VO2 film 56 (See Fig. 10(a)). Two thin Au pads are deposited on each patterned VO2 layer which serves as the local bias elements. To apply current to each heating electrode, three metallic via holes are drilled through the top substrate to connect one of the two electrical contacts for each patterned VO2 with the ground layer as a negative electrode. The other three metallics via holes are drilled through the top and bottom substrates to attach other electrical contacts with three pieces of separated patches on the bottom metallic layer which acts as positive electrodes to apply the DC bias current. To ensure that the positive metallic via hole is electrically isolated from the negative electrode (ground plane), a metal ring can be subtracted from the second metallic layer. Figure 10(b,c) shows the two-dimensional view of the top metallic layer (ground plane) as a negative electrode and three pieces of separated patches on the bottom metallic layer which act as positive electrodes to apply the DC bias current. It is important to note that the reflection spectra of the unit-cell is not influenced by adding a new substrate and metallic patches underneath the ground plane. Moreover, very thin heater electrodes and biasing lines with respect to the operating wavelength (analogous to this paper), will not perturb the far-field pattern and therefore have a negligible impact on the macroscopic behavior of our designed VBCM and the individual response of meta-atoms 56,123,124 .