Abstract
Tunable multi-function metasurfaces have become the latest research frontiers in planar optics. In this study, a dynamically tunable plasmon-induced transparency (PIT) structure based on a graphene split-ring resonator and graphene ribbon is proposed. The influences of the structural parameters and graphene Fermi energy on the PIT response were investigated both analytically and numerically simulations. The inclusion of an additional vanadium dioxide (VO2) substrate layer enables the metasurface to achieve dynamic switching between PIT and perfect absorption using the phase change property of VO2. The new metasurface device exhibits the PIT effect when the VO2 layer is in an insulating state and acts as a perfect absorber when it is in a metallic state. Moreover, the response of the two functions can be easily adjusted dynamically by changing the Fermi energy of graphene. In addition, both functions were highly sensitive to changes in the ambient refractive index. The results of this work have potential applications in slow-light devices, optical switches, modulators, perfect absorbers, highly sensitive sensors, and multifunctional devices.
© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
Metasurfaces are artificial array structures comprising a large number of sub-wavelength units arranged periodically or non-periodically on a two-dimensional plane, which can control electromagnetic waves flexibly [1]. Over the past few decades, metasurfaces have rapidly developed into powerful electromagnetic field regulation tools to realize various novel optical phenomena, such as electromagnetically induced transparency (EIT) [2–4] and perfect absorption [5,6], which have important application value in many fields because they possess a sub-wavelength ultrathin structure and are highly amenable to design. EIT refers to a phenomenon in which a narrow transparency window is formed within a broad absorption spectrum. It was first discovered in a three-atom system under harsh experimental conditions, including extremely low temperature, high-intensity lasers, and appropriate atomic energy level systems [7]. In contrast, plasmon-induced transparency (PIT), an analog of EIT, has no strict implementation conditions and can be simulated using metamaterials. Therefore, metasurface-based PIT has attracted great attention from researchers and has been widely used for realizing optical switches [8], modulators [9], and slow optical devices [10]. Perfect absorption refers to the complete absorption of incident light of a certain wavelength by the material. A metasurface-based perfect absorber can be achieved when the impedance of the absorber matches that of free space, which can be realized by designing metasurface cells of the appropriate structure and shape. Until now, metasurface-based perfect absorbers have been widely used in sensors [11,12] and thermal radiation detectors [13] for various wavebands.
The latest research frontier in the fields of PIT/EIT and perfect absorption is designing tunable metasurfaces to achieve tunable PIT/EIT devices [14,15] and perfect absorbers [16,17], which can significantly improve the practicability and compatibility of the devices in practical applications. This can be achieved by including materials with adjustable properties such as phase change materials [18], Dirac semimetals [19–21], and graphene [22,23]. Significant efforts have been made by researchers in this regard. For example, Wei et al. realized a dual- polarization PIT metasurface design based on graphene and black phosphorus, and both polarizations can be tuned [15]. Yao et al. studied the broadband absorption of two-layer graphene under electrical modulation and obtained an absorption bandwidth that is 1.7 times higher than that of single-layer graphene [24]. Song et al. proposed a VO2 based dual-function device, which can switch from EIT to broadband absorption when VO2 transitions from an insulator to a metal [25]. However, some of these studies focused only on the tunability of a single function [26–29], while the others studied only the switching of dual functions [30–32]. Therefore, designing devices that can simultaneously switch between PIT and perfect absorption and adjust the response frequency of the two functions separately is the next important step in this field, which motivated our work.
We propose a simple metasurface structure comprising a monolayer graphene split-ring resonator (GSRR) and graphene ribbon (GR) to realize PIT. The mechanism of PIT and the influences of the structural parameters and the Fermi energy of graphene on the PIT are analyzed and discussed. The inclusion of an additional VO2 substrate layer enables the metasurface to dynamically switch between PIT and perfect absorption by using the phase change property of VO2. When the conductivity of VO2 is low, the device can realize the PIT. When the conductivity of VO2 is high enough to block the transmission, the device acts as a perfect absorber. The response frequency of both functions can be adjusted dynamically by changing the Fermi energy of graphene. Furthermore, since the proposed bifunctional metasurface is highly sensitive to changes of the ambient refractive index, it can be used for sensing. Our work achieves switching between different functions and the dynamic modulation of each function. This may positively influence research of multifunctional devices and be used in slow-light devices, sensing, and perfect absorption.
2. Material models
The monolayer graphene was modeled as an equivalent two-dimensional impedance layer, whose impedance is Z(ω) = 1/σ(ω). According to the Pauli Exclusion Principle, the complex surface conductivity of graphene, in the THz band and at room temperature, is given by [33]
The permittivity of VO2 in the THz band can be obtained using the Drude model [34]:
3. Graphene-based PIT structures
In this section, we describe the study of the graphene-based PIT structure. The schematic of the proposed PIT structure is shown in Fig. 1. GSRR and GR are deposited on a silicon dioxide (SiO2) substrate whose thickness and permittivity are 0.2 µm and 3.9 [28,37], respectively. A layer of ion-gel (whose thickness and permittivity are 0.1 µm and 1.82, respectively) on the surface of graphene is used to deposit electrodes (Au) to control the carrier density of graphene [29,38]. The structural parameters are R = 2.4 µm, l = 3 µm, W1 = 0.8 µm, W2 = 0.9 µm, and P = 6 µm. The distance between GSRR and GR is d, and the opening angle of the GSRR is α. CST Microwave Studio, a full-wave electromagnetic simulation software based on a finite integration technique, was used to simulate the proposed structure. In the simulation, the boundary conditions in the x and y directions are provided by the unit cells, and the z direction is open (add space). The incident wave is plane wave and polarized along the x-axis throughout this study. The transmission can be obtained from T=|S21|2, where S21 represents transmission coefficient. Unless emphasized, the incident electromagnetic waves are incident normally.
The simulated transmission curves of a single GSRR (dashed red line), single GR (dashed blue line) and their combination (solid black line) corresponding to EF, α, and d of 0.6 eV, 90°, 0.7 µm, respectively, are shown in Fig. 2(a). Compared to those for the single GSRR and GR structures, the frequency of dip1 in the combined structure red shifts, while that of dip2 blue shifts owing to the coupling between GSRR and GR. A high transmission window appears between dip1 and dip2 in the combined structure. To explore the physical mechanism of this phenomenon, we simulated the z-components of the electric field (Ez) distributions at the frequencies corresponding to dip1 (3.6 THz), dip2 (4.42 THz), and the transmission peak (3.96 THz), as shown in Figs. 2(b)–(d). At 3.6THz, a strong electric dipole resonance occurs at both arms of GSRR, and GR has a weak response to THz wave, as shown in Fig. 2(b). By contrast, at 4.42 THz, GR exhibits a strong electric dipole response, and GSRR has a weak response, as shown in Fig. 2(d). At 3.96THz, due to resonance detuning, GSRR and GR have a very weak response at the same time and their induce currents are out of phase, which is the characteristic of the electromagnetic-trapped mode [14]. Due to the opposite current oscillations of GSRR and GR, the scattered field at 3.96THz is very weak. Therefore, the high transmission peak appears at 3.96 THz, as shown in Fig. 2(c). Figures 2(c) and (f) are the electric field distributions of a single GSRR and GR at 3.96THz, respectively. GSRR has a strong response while the GR response is not obvious. However, the response of their combination is suppressed, as shown in Fig. 2(c). This proves our conclusion: the anti-phase current oscillations of GSRR and GR at 3.96 cause the electric dipole moments of the two structures to cancel each other.
The two-particle model which can be used to describe the coupling between GSRR and GR treats them as two particles interacting with the incident wave. The interaction between them can be expressed as [37,40]:
In the fitting, the transmittance is defined as T = 1-Im(χ) (using the Kramer–Kronig relations). γ1, γ2 can be obtained from the linewidths of the separate GSRR and GR transmission curves. The coupling strength κ is obtained by solving the equations ${\omega _1} = \sqrt {\omega _0^2 - {\kappa ^2}} {\; }$and ${\omega _2} = \sqrt {\omega _0^2 + {\kappa ^2}} $, where ω1 and ω2 represent the band-stop frequencies. Finally, as shown in Fig. 3(a), the transmission curves of the analytic model with EF of 0.4 eV, 0.6 eV, and 0.8 eV are obtained. The curve obtained using the analytic model is evidently consistent with the simulated curve. As EF increases, the transmission window undergoes a blue shift and widens. This phenomenon is more evident in Fig. 3(b). The fitting parameters, such as the coupling strength and damping rates, are plotted in Fig. 3(c). As EF increases, γ1 and γ2 remain almost unchanged, while κ shows an obvious increase. These imply that increasing the Fermi energy strengthens the coupling between GSRR and GR.
We then discuss the influences of the structural parameters d and α on the transmission curve. Changing d essentially changes the coupling strength between GSRR and GR. When compared to the single GSRR and GR structures, the transmission curve of the combined structure exhibits the effect of the coupling between GSRR and GR via red shifts and blue shifts in the frequencies of dip1 and dip2, respectively. Therefore, when the coupling strength is enhanced by reducing d, the frequencies of dip1 and dip2 will be further red shifted and blue shifted, respectively. However, when the coupling strength is reduced by increasing d, the frequencies of dip1 and dip2 will undergo blue shifts and red shifts, respectively. A change in α affects the resonance frequency of GSRR. According to the formula of the plasmon resonance frequency for graphene f = e/(2πћ)(EF/πD)1/2 [41], when α increases, the effective length D of graphene parallel to the polarization direction decreases, leading to an increase in f, that is, the dip1 blue shifts when α increases. The increase in α does not affect the coupling between GSRR and GR; therefore, dip2 is expected to remain unchanged when α changes. To verify this, we simulated the transmission curves using different parameters. The transmission curves corresponding to different values of d at α= 90° are shown in Fig. 4(a). The figure shows that when d decreases from 1.1 µm to 0.3 µm, dip1 exhibits a red shift and its transmission increases, while dip2 exhibits a blue shift and its transmission decreases. This observation is consistent with our expectations. The transmission curves corresponding to different values of α at d= 0.7 µm are shown in Fig. 4(b), where an increase in α is observed to produce a blue shift in dip1. During this process, dip2 remains unchanged. These phenomena are also consistent with our analysis.
Since electromagnetic waves are not always normally incident into the structure, we simulated the transmission curve at different incident angles θ (in the y-z plane) assigning the values of 0.6 eV, 90°, and 0.7 µm to EF, α, and d, respectively. As shown in Figs. 5(a) and (b), the frequencies of the transmission peak and transmission valleys remain unchanged as θ increases, while the transmission decreases. This can be attributed to two factors. First, the change in the incident angle does not change the electric field component in the x-direction, so it does not affect the resonance frequencies and coupling strength of GSRR and GR, and the frequencies of the transmission peak and the transmission valleys remain unchanged. Second, an increase in the incident angle increases the reflection of the incident wave resulting in decreased transmission. Moreover, when the incident angle is less than 30°, the maximum change in transmission (3%) indicates insensitivity to small incident angles, which is beneficial for practical applications.
It is worth noting that PIT is accompanied by strong dispersion and slow light effect, which makes it important for applications including slow optical devices and optical routing. The group delay (τd) is used to quantify the slow-light effect. It is defined as the derivative of the transmission phase φ(ω) as a function of angular frequency ω, that is,
Figure 6(a) clearly demonstrates that as EF increases from 0.4 eV to 0.8 eV, the peak value of τd increases from 0.28 ps to 0.32 ps. When EF is fixed at 0.6 eV, τd increases monotonously from 0.30ps to 1.02ps with the incident angle θ increases from 0° to 85°, as shown in Fig. 6(b). Because a larger incident angle leads to a longer propagation distance in the structure resulting in a larger group delay. In sum, active tuning of the group delay can be achieved by controlling the EF of graphene and incident angle. Moreover, when θ increases from 0° to 30°, the group delay only increases by 1.2%, showing insensitivity to small incident angles. These tunable features are essential for developing integrated devices.
4. Tunable and switchable bifunctional metasurface for PIT and perfect absorption
In this section, we describe the effect of adding an additional VO2 substrate layer to realize a switchable PIT and perfect absorption bifunctional device. The optimized thickness of VO2 is 0.8 µm while the SiO2 layer is 4 µm thick. The other parameters remain unchanged. The new structure is shown in Fig. 7. The conductivity of VO2 can be adjusted independently by the applied voltage [27,31].
In this structure, three tunable parameters, the graphene Fermi energy, incident angle, and conductivity of VO2, can be used to actively adjust the transmission curve. The transmission curves and group delay curves under different graphene Fermi energies and incident angles obtained when VO2 is in the insulating state (σ =10 S/m), are shown in Figs. 8(a)–(d). In this case, the device exhibits PIT. The transmission peak blue shifts and decreases, the group delay of the transmission peak increases as EF increases. This behavior can be explained in two factors. First, the resonant frequency f positively related to the EF of graphene. Secondly, the transmission of THz waves pass through VO2 gradually decreases with the frequency increases. The transmission of the transmission peak decreases, and the group delay of the transmission peak increases with an increase in θ under EF = 0.3 eV, as shown in Figs. 8(c) and (d). This can be attributed to the increase in the incident angle increases the propagation distance of THz waves in the structure, which is equivalent to increasing the thickness of the structure, so the transmission decreases and group delay increases. The transmission curves and group delay curves under different σ when EF = 0.3 eV and θ=0° are shown in Figs. 8(e)–(f). It can be seen that the transmission and group delay decrease as σ increases, but the frequency of the peak remains unchanged. When σ=200000 S/m, VO2 becomes metallic, and blocks the transmission of electromagnetic waves.
As the transmission of electromagnetic waves is blocked when VO2 is in the metallic state, we simulated the absorption curves at EF = 1.2 eV. The absorption is defined as A = 1-|S11|2-|S21|2 = 1-|S11|2, where S11 and S21 represent complex reflection and transmission coefficients, respectively (|S21|=0 when VO2 is in the metallic state). Two nearly perfect absorption peaks at 3.94 THz and 5.04 THz can be clearly observed. We introduce the coupled mode theory (CMT) and compare it with the simulation results [28,43].
where a is the amplitude of resonance with the resonance frequency ${\omega _0}$. $\gamma $ and $\delta $ are the radiation loss and intrinsic loss rates of the resonance, respectively. ${S_ + }$ and ${S_ - }$ are amplitudes of the input and output THz waves, respectively. Finally, the absorption can be expressed as:Next, we studied the tunability of the absorber. When EF is fixed at 1.2 eV, the absorption in the frequency range 2-6 THz increases significantly as the conductivity of VO2 decreases, except for the two peaks, as shown in Figs. 10(a) and (b). We attribute this to the decrease in conductivity, which weakens the metallicity of VO2, resulting in the increased absorption of electromagnetic waves by VO2. As shown in Fig. 10(c), the absorption peak blue shifts with the increase in EF and absorption gradually increases because the resonance frequency is positively correlated with the graphene Fermi energy. In addition, the tolerance to the incident angle θ is discussed, as shown in Fig. 10(d). When EF= 1.2 eV is fixed and the incident angle θ is gradually increased, it is surprising to note that the absorption peaks remain unchanged at 0°-80° and maintain more than 90% absorption between 0°-58°. This demonstrates that our structure has a high tolerance to the incidence angle, which is very important for practical applications.
It is well known that the transmission peak of PIT is strongly dispersive and very sensitive to the refractive index. Therefore, we explored the influence of the ambient refractive index on PIT and perfect absorption. The response frequencies of the two functions vary with the refractive index n, as shown in Fig. 11. The sensing performance is quantified by calculating the sensitivity S, which is defined as the ratio of the variation of the wavelength to the refractive index ($S = {\Delta }\lambda /{\Delta }n$). We calculated the sensitivity S of the PIT as S= 18219 nm/RIU, while those of the two absorption peaks are S1= 20541 nm/RIU and S2= 23166 nm/RIU. The subscripts 1 and 2 correspond to peaks 1 and 2, respectively. For comparison, we list similar published studies in Table 1, which indicate that our proposed structure can switch between the two functions and also exhibit better sensing performance than other graphene sensors. This can be attributed to two factors. First, our structure is composed of graphene split-ring resonator and graphene ribbon. The response of an isolated structure is intense than that in other references. Secondly, our graphene structure occupies a larger proportion of the entire unit structure. Combined with our previous analysis, it can be found that our structure is highly sensitive to the refractive index and insensitive to the incident angle, and we believe that it will have good application potential in refractive index sensors.
5. Conclusion
In summary, a graphene-based PIT structure in the THz range is proposed. The influences of the structural parameters, graphene Fermi energy, and incident angle on the PIT response and group delay are investigated in detail. The proposed structure is then turned into a bifunctional device by including a VO2 substrate. The new structure can switch between PIT and perfect absorber using the phase change property of VO2. PIT occurs when the VO2 is in an insulating state, and two perfect absorption peaks occur when the VO2 is in the metallic state. The responses of the two functions can be adjusted dynamically by changing the Fermi energy of graphene. Moreover, the responses of the two functions are insensitive to the incident angle, which is suitable for practical applications. The influence of the ambient refractive index on PIT and absorber was also explored. It offers potential applications for refractive index sensing. We expect our work to promote research on multifunctional metasurfaces and their applications.
Funding
National Natural Science Foundation of China (62075047, 61965006, 61975038, 61964005,62065006); Natural Science Foundation of Guangxi Province (2020GXNSFDA297019, 2020GXNSFAA238040, 2021GXNSFAA075012); Science and Technology Project of Guangxi (AD19245064); Guangxi Key Laboratory of Automatic Detecting Technology and Instruments (YQ20107); and the Innovation Project of GUET Graduate Education (2020YCXS089).
Disclosures
The authors declare no conflicts of interest.
Data availability
Data on the relaxation time of graphene underlying the results presented in this paper are available in Ref. [26]. Data on the conductivity of VO2 underlying the results presented in this paper are available in Ref. [28,29]. The data of the permittivity of SiO2 underlying the results presented in this paper are available in Ref. [21,30].
References
1. Q. He, S. Sun, S. Xiao, and L. Zhou, “High-efficiency metasurfaces: Principles, realizations, and applications,” Adv. Opt. Mater. 6(19), 1800415 (2018). [CrossRef]
2. C. L. Garrido Alzar, M. A. G. Martinez, and P. Nussenzveig, “Classical analog of electromagnetically induced transparency,” Am. J. Phys. 70(1), 37–41 (2002). [CrossRef]
3. V. A. Fedotov, M. Rose, S. L. Prosvirnin, N. Papasimakis, and N. I. Zheludev, “Sharp trapped-mode resonances in planar metamaterials with a broken structural symmetry,” Phys. Rev. Lett. 99(14), 147401 (2007). [CrossRef]
4. S. Zhang, D. A. Genov, Y. Wang, M. Liu, and X. Zhang, “Plasmon-induced transparency in metamaterials,” Phys. Rev. Lett. 101(4), 047401 (2008). [CrossRef]
5. N. I. Landy, S. Sajuyigbe, J. J. Mock, D. R. Smith, and W. J. Padilla, “Perfect metamaterial absorber,” Phys. Rev. Lett. 100(20), 207402 (2008). [CrossRef]
6. J. Grant, Y. Ma, S. Saha, L. B. Lok, A. Khalid, and D. R. S. Cumming, “Polarization insensitive terahertz metamaterial absorber,” Opt. Lett. 36(8), 1524–1526 (2011). [CrossRef]
7. L. Mao, Y. Li, G. Li, S. Zhang, and T. Cao, “Reversible switching of electromagnetically induced transparency in phase change metasurfaces,” Adv. Photonics 2(5), 056004 (2020). [CrossRef]
8. J. Liu, I. Papakonstantinou, H. Hu, and X. Shao, “Dynamically configurable, successively switchable multispectral plasmon-induced transparency,” Opt. Lett. 44(15), 3829–3832 (2019). [CrossRef]
9. W. Wang, L. Du, J. Li, M. Hu, C. Sun, Y. Zhong, G. Zhao, Z. Li, L. Zhu, J. Yao, and F. Ling, “Active control of terahertz waves based on p-Si hybrid PIT metasurface device under avalanche breakdown,” Opt. Express 29(8), 12712–12722 (2021). [CrossRef]
10. T. Zentgraf, S. Zhang, R. F. Oulton, and X. Zhang, “Ultranarrow coupling-induced transparency bands in hybrid plasmonic systems,” Phys. Rev. B 80(19), 195415 (2009). [CrossRef]
11. N. Liu, M. Mesch, T. Weiss, M. Hentschel, and H. Giessen, “Infrared perfect absorber and its application as plasmonic sensor,” Nano Lett. 10(7), 2342–2348 (2010). [CrossRef]
12. J. Nong, L. Tang, G. Lan, P. Luo, Z. Li, D. Huang, J. Shen, and W. Wei, “Combined Visible Plasmons of Ag Nanoparticles and Infrared Plasmons of Graphene Nanoribbons for High-Performance Surface-Enhanced Raman and Infrared Spectroscopies,” Small 17(1), 2004640 (2021). [CrossRef]
13. F. Niesler, J. Gansel, S. Fischbach, and M. Wegener, “Metamaterial metal-based bolometers,” Appl. Phys. Lett. 100(20), 203508 (2012). [CrossRef]
14. Z. Jia, L. Huang, J. Su, and B. Tang, “Tunable electromagnetically induced transparency-like in graphene metasurfaces and its application as a refractive index sensor,” J. Lightwave Technol. 39(5), 1544–1549 (2021). [CrossRef]
15. P. Luo, W. Wei, G. Lan, X. Wei, L. Meng, Y. Liu, J. Yi, and G. Han, “Dynamical manipulation of a dual-polarization plasmon-induced transparency employing an anisotropic graphene-black phosphorus heterostructure,” Opt. Express 29(19), 29690–29703 (2021). [CrossRef]
16. J. Nong, L. Tang, G. Lan, P. Luo, C. Guo, J. Yi, and W. Wei, “Wideband tunable perfect absorption of graphene plasmons via attenuated total reflection in Otto prism configuration,” Nanophotonics 9(3), 645–655 (2020). [CrossRef]
17. G. Lan, W. Wei, P. J. Yi, Z. Shang, and T. Xu, “Dynamically tunable coherent perfect absorption in topological insulators at oblique incidence,” Opt. Express 29(18), 28652–28663 (2021). [CrossRef]
18. J. Huang, J. Li, Y. Yang, J. Li, J. Li, Y. Zhang, and J. Yao, “Broadband terahertz absorber with a flexible, reconfigurable performance based on hybrid-patterned vanadium dioxide metasurfaces,” Opt. Express 28(12), 17832–17840 (2020). [CrossRef]
19. F. Xia, H. Wang, D. Xiao, M. Dubey, and A. Ramasubramaniam, “Two-dimensional material nanophotonics,” Nat. Photonics 8(12), 899–907 (2014). [CrossRef]
20. M. Liu, H. Y. Hwang, H. Tao, A. C. Strikwerda, K. Fan, G. R. Keiser, A. J. Sternbach, K. G. West, S. Kittiwatanakul, J. Lu, S. A. Wolf, F. G. Omenetto, X. Zhang, K. A. Nelson, and R. D. Averitt, “Terahertz-field-induced insulator-to-metal transition in vanadium dioxide metamaterial,” Nature 487(7407), 345–348 (2012). [CrossRef]
21. O. V. Kotov and Y. E. Lozovik, “Dielectric response and novel electromagnetic modes in three-dimensional Dirac semimetal films,” Phys. Rev. B 93(23), 235417 (2016). [CrossRef]
22. J. Nong, L. Tang, G. Lan, P. Luo, Z. Li, D. Huang, J. Yi, H. Shi, and W. Wei, “Enhanced Graphene Plasmonic Mode Energy for Highly Sensitive Molecular Fingerprint Retrieval,” Laser Photonics Rev. 15(1), 2000300 (2021). [CrossRef]
23. W. Yao, L. Tang, J. Wang, Y. Jiang, and X. Wei, “Anomalous redshift of graphene absorption induced by plasmon-cavity competition,” Opt. Express 28(25), 38410–38418 (2020). [CrossRef]
24. W. Yao, L. Tang, J. Nong, J. Wang, J. Yang, Y. Jiang, H. Shi, and X. Wei, “Electrically tunable graphene metamaterial with strong broadband absorption,” Nanotechnology 32(7), 075703 (2021). [CrossRef]
25. Z. Song and B. Zhang, “Controlling wideband absorption and electromagnetically induced transparency via a phase change material,” Europhys. Lett. 129(5), 57003 (2020). [CrossRef]
26. T. T. Kim, H. D. Kim, R. Zhao, S. S. Oh, T. Ha, D. S. Chung, Y. H. Lee, B. Min, and S. Zhang, “Electrically tunable slow light using graphene metamaterials,” ACS Photonics 5(5), 1800–1807 (2018). [CrossRef]
27. C. Zhang, G. Zhou, J. Wu, Y. Tang, Q. Wen, S. Li, J. Han, B. Jin, J. Chen, and P. Wu, “Active control of terahertz waves using vanadium-dioxide-embedded metamaterials,” Phys. Rev. Appl. 11(5), 054016 (2019). [CrossRef]
28. L. Qi, C. Liu, and S. M. A. Shah, “A broad dual-band switchable graphene-based terahertz metamaterial absorber,” Carbon 153, 179–188 (2019). [CrossRef]
29. M. Chen, Z. Xiao, X. Lu, F. Lv, and Y. Zhou, “Simulation of dynamically tunable and switchable electromagnetically induced transparency analogue based on metal-graphene hybrid metamaterial,” Carbon 159, 273–282 (2020). [CrossRef]
30. H. Li and J. Yu, “Bifunctional terahertz absorber with a tunable and switchable property between broadband and dual-band,” Opt. Express 28(17), 25225–25237 (2020). [CrossRef]
31. L. Liu, L. Kang, T. S. Mayer, and D. H. Werner, “Hybrid metamaterials for electrically triggered multifunctional control,” Nat. Commun. 7(1), 13236 (2016). [CrossRef]
32. F. Ding, S. Zhong, and S. I. Bozhevolnyi, “Vanadium dioxide integrated metasurfaces with switchable functionalities at terahertz frequencies,” Adv. Opt. Mater. 6(9), 1701204 (2018). [CrossRef]
33. G. Yao, F. Ling, J. Yue, Q. Luo, and J. Yao, “Dynamically tunable graphene plasmon-induced transparency in the terahertz region,” J. Lightwave Technol. 34(16), 3937–3942 (2016). [CrossRef]
34. S. Wang, L. Kang, and D. Werner, “Hybrid resonators and highly tunable terahertz metamaterials enabled by vanadium dioxide (VO2),” Sci. Rep. 7(1), 4326 (2017). [CrossRef]
35. T. Wang, Y. Zhang, H. Zhang, and M. Cao, “Dual-controlled switchable broadband terahertz absorber based on a graphene-vanadium dioxide metamaterial,” Opt. Mater. Express 10(2), 369–386 (2020). [CrossRef]
36. Q. Wen, H. Zhang, Q. Yang, Y. Xie, K. Chen, and Y. Liu, “Terahertz metamaterials with VO2 cut-wires for thermal tunability,” Appl. Phys. Lett. 97(2), 021111 (2010). [CrossRef]
37. X. He, X. Yang, G. Lu, W. Yang, F. Wu, Z. Yu, and J. Jiang, “Implementation of selective controlling electromagnetically induced transparency in terahertz graphene metamaterial,” Carbon 123, 668–675 (2017). [CrossRef]
38. M. Liu, W. Cheng, Y. Zhang, H. Zhang, Y. Zhang, and D. Li, “Multi-controlled broadband terahertz absorber engineered with VO2-integrated borophene metamaterials,” Opt. Mater. Express 11(8), 2627–2638 (2021). [CrossRef]
39. B. Xiao, S. Tong, A. Fyffe, and Z. Shi, “Tunable electromagnetically induced transparency based on graphene metamaterials,” Opt. Express 28(3), 4048–4057 (2020). [CrossRef]
40. R. Yahiaoui, M. Manjappa, Y. K. Srivastava, and R. Singh, “Active control and switching of broadband electromagnetically induced transparency in symmetric metadevices,” Appl. Phys. Lett. 111(2), 021101 (2017). [CrossRef]
41. F. Garcia de Abajo, “Graphene plasmonics: challenges and opportunities,” ACS Photonics 1(3), 135–152 (2014). [CrossRef]
42. P. Tang, J. Li, L. Du, Q. Liu, Q. Peng, J. Zhao, B. Zhu, Z. Li, and L. Zhu, “Ultrasensitive specific terahertz sensor based on tunable plasmon induced transparency of a graphene micro-ribbon array structure,” Opt. Express 26(23), 30655–30666 (2018). [CrossRef]
43. Y. Cai, Y. Guo, Y. Zhou, X. Huang, G. Yang, and J. Zhu, “Tunable dual-band terahertz absorber with all-dielectric configuration based on graphene,” Opt. Express 28(21), 31524–31534 (2020). [CrossRef]