Graphene Metamaterials Based Plasmon-Induced Terahertz Modulator for High-Performance Multiband Filtering and Slow Light Applications

We proposed multilayered graphene (Gr)-based surface plasmon resonance-induced high-performance terahertz (THz) modulators with tunable resonance frequencies. Several THz plasmonic modulators based on Gr metamaterials were previously reported; however, these modulators had small group delay, low extinction ratio (ER), and difficult-to-tune resonant frequency without structural parameters in the THz range. A comprehensive investigation employing the finite-difference time-domain (FDTD) simulation technique revealed high group delay, broad tunability independent of structural parameters, and large ER for our proposed quadband and pentaband plasmonic modulators. We obtained tunable group delays with a maximum of 1.02 ps and 1.41 ps for our proposed quadband and pentaband plasmonic modulators, respectively, which are substantially greater compared to previously reported Gr-based metamaterial structures. The maximum ER of 22.3 dB was obtained which was substantially high compared to previous reports. Our proposed modulators were sensitive to the polarization angle of incident light; therefore, the transmittance at resonant frequencies was increased while the polarization angle varied from 0 to 180 degree. These high-performance plasmonic modulators have emerging potential for the design of optical buffers, slow light devices, multistop band filters, integrated photonic circuits, and various optoelectronic systems.


I. INTRODUCTION
T HE invention of metamaterials has garnered many aca- demic and technological interests due to their extraordinary electromagnetic (EM) properties, which are not easily obtained in the existing photonics and optics world.Plasmonic metamaterials utilize surface plasmons (SPs) which exhibit the most compelling EM properties among the ramifications of the metamaterial [1].Structures comprising noble metals, such as gold and silver, can overcome the limitations of classical diffraction limit and realize a photonic loop in the subwavelength range of light waves; however, significant energy loss limits practical applicability due to the poor binding effect on incoming light and scattering effect [2].Gr, a twodimensional material with a hexagonal lattice, has similar properties as noble metals in the THz and infrared regime [3].Most importantly, Gr can excite SPs with strong localization ability, low transmission loss, high tunability, and high group delay.These properties make Gr metamaterial-based structures advantageous in applications including optical cloaking [4], perfect lens [5], modulators [6], photodetectors [7], polarizers [8], bandstop filtering [9], and slow light devices [10].
Bandstop filters are essential components of information processing systems.Recently, plasmonic bandstop filters have received considerable attention from researchers.Conventional metal-based bandstop filters operate well in the visible and near-infrared ranges; however, their applications are limited by the poor SP confinement in the THz frequency range.Contrarily, Gr-based bandstop filters exhibit strong SP confinement providing excellent performance in the THz spectral band.Tunability of resonant frequency without varying structural parameters is another distinguishing attribute for Grbased structures, whereas tuning characteristics are not easily obtained in metal-based structures where liquid crystals or liquid metals with variable permittivity have been proposed.However, tunning is achievable by simply changing the chemical potential of Gr with an external gate voltage due to Gr's unique energy band structure.This allows fine-tuning of resonant frequencies and group velocities in Gr-based bandstop filters.Graphene nanoribbon (GNR)-Gr disk [11], multilayer Gr metamaterial [12], continuously patterned Gr monolayer [13], [14], H-shape Gr-metamaterial [15], split-ring resonator [16], and asynchronous crisscross structures [17] were reported; however, extinction ratio and quality factor, two performance parameters of bandstop plasmonic filter, were quite small.Therefore, there is a scope for better Gr metamaterial-based design to improve the performance of the bandstop optical filters.
Gr-based plasmonic structures play a vital role in slow light applications, such as optical buffers and quantum memory, by the tuning properties of SPs originated by the interaction of Gr and incoming light.A slow light device can slow down the speed of light pulses propagating through it.The concept behind slow light is based on the phenomenon of light propagation in a medium with a refractive index that can be manipulated by increasing the group delay of light [18], [19].Group delay, a performance parameter of slow light devices, is a measure of how much the light in a medium is delayed compared to the speed of light in the vacuum.Multilayered Gr and different shapes of GNRs exhibited slow light properties due to light-matter interactions; however, the group delay and delay-bandwidth product of the previous reports were low [12], [13], [16].Zhao et al. proposed a ring metamaterial structure comprising noble metals, which provided a high group delay; however, tunability can not be possible without changing structural parameters [20].Thus, there is a scope for significant improvement in slow light device performance with wide tunability for THz interconnects, quantum memory, and optical buffer applications.
In this work, we proposed two multiband THz modulator structures comprising GNR, Gr nanoholes, and polymethyl methacrylate (PMMA).We employed the FDTD simulation method to optimize the structures and evaluated the performance parameters of the devices.These structures exhibited dips in their transmittance spectra that arose from the interaction between the incoming light and Gr-based metamaterials.Four-and five-patterned Gr layers on PMMA substrate were utilized for the quadband and pentaband plasmonic modulators, respectively, where each Gr layer acted as an attenuating medium.Our proposed modulators achieved high ER and group delay at resonant frequencies due to strong light-matter interactions making it promising for slow light applications.Moreover, the plasmonic modulators demonstrated high tunability by varying the chemical potentials of Gr layers, which acted as a perfect multi-bandstop filter.We analyzed the polarization sensitivity of our proposed modulators when they were used as bandstop filters.Additionally, we provided a comprehensive comparative study of performance parameters among previously reported Gr-based and our proposed THz modulator structures.

II. METHODOLOGY
A unit cell of the Gr-based plasmonic quadband modulator is illustrated in Fig. 1(a).Contrary to the conventional modulator designs, the proposed structure used Gr's excellent SPP and metal-like properties to enhance the modulator's performance.The quadband modulator was designed by comprising four Gr monolayers, which were made up of GNR and Gr nanoholes on PMMA substrate as illustrated in Fig. 1(a).Here, the substrate height was selected to be 2 µm to avoid plasmonic coupling and signal interference.We optimized the GNR arrangement of the modulator's top Gr layer, as depicted in Fig. S1 of Supplementary Material.Twelve periodic airstrips or eleven GNRs on the Gr monolayer were adopted for our proposed modulators.The chemical vapor deposition (CVD) technique can grow a such Gr monolayer on PMMA [21].Fig. 1(b) presents the GNR layer of the proposed structures with length, q of 2.5 µm, width, w of 10 nm, and pitch, d of 20 nm.Inset shows the zoomed view of the GNR layer.The period, p, of the proposed structures, was set to be 3 µm.The widths of the nanoholes, i = j = k of 150 nm in the rest of the Gr layers, were adopted with length, q of 2.5 µm.The helium ion beam lithography technique can be utilized to precisely create GNR because of its ability to fabricate GNR below five nm [22].To design a pentaband plasmonic modulator, we utilized an additional Gr layer with a nanohole below, where the width and length of the Gr nanohole were the same as before.Fig. 1(c) illustrates the schematic representation of a unit cell of a pentaband plasmonic modulator.The surface current, J, of the Gr layer was modeled by employing the equation of Ohm's law: J = σ g E where σ g and E are the Gr's complex surface conductivity and electric field, respectively.We modeled the σ g of Gr by using Kubo's formulation (see Section 2 of Supplementary Material) [23].This formulation was modeled specifically for the σ g of Gr monolayer; however, it can be conducive to explaining multilayer Gr's σ g [24].The interband conductivity term can be neglected due to E F ≫ k B T in the THz frequency regime where E F , k B , and T denote the Fermi energy level, Boltzmann constant, and temperature, respectively.The σ g of Gr is given by, .
Here, e, , and ω are the electron charge, reduced Plunck's constant, and angular frequency, respectively.τ is the relaxation time of the Gr layer, which can be expressed by τ = µE F /eV F 2 .The Fermi velocity, V F , and mobility, µ were set to be 10 6 m/s, and 0.4 m 2 V −1 s −1 , respectively [25].The real and imaginary surface conductivity of Gr as a function of E F and frequency are shown in Figs.2(a) and (b).The imaginary part of the Gr conductivity rises while the real part stays constant as E F increases.Therefore, Gr exhibits more metal-like properties with the increase of E F producing SPPs.Figs.2(c) and (d) illustrate the real and imaginary surface conductivity of the Gr layer at E F = 0.6 eV.The refractive index data of PMMA was adapted from Lee et al. [26].We conducted the study utilizing the three-dimensional (3D) FDTD simulation technique.The 3D FDTD simulation range was set to be 3 µm × 3 µm × 35 µm.We utilized the periodic boundary conditions in the x-and y-directions, and 24-steep angle perfectly matched layers (PMLs) in the zdirection to absorb light after the simulation region without back reflections.We employed a transverse magnetic (TM) polarized plane-wave light source on the top of the proposed modulators along the z-direction with a center frequency of 3 THz and a frequency span of 5.7 THz.To consider the room temperature for our study, the temperature of the simulation environment was set to 300 K. FDTD is a time-varying and time-consuming technique: larger mesh sizes shorten the computational time which results in low accuracy, and smaller mesh sizes provided high accuracy, increasing computational time and memory size.In our work, smaller mesh sizes were utilized in the region of interest to resolve tiny features, whereas bigger mesh sizes were used in other regions.The non-uniform mesh was utilized to complete simulations in a reasonable amount of time with maintaining the memory requirements.We utilized a mesh size of 8 nm in the GNR and Gr nanoholes for better accuracy in our simulations.Nearfield monitors were utilized to obtain the spatial electric field profiles.Transmittance and reflectance spectra were recorded by using power monitors.

III. RESULTS AND DISCUSSION
Interaction between the incident plane wave's electric field and free electrons or plasmons of the Gr layer resulted in SP resonance on the surface of the Gr layers of the proposed modulators.As a result, the modulators attenuated light by scattering that created dips in transmittance spectra at resonant frequencies.Based on this phenomenon, we utilized each Gr layer as an attenuating medium at the resonant frequency.Most importantly, this phenomenon will be beneficial for designing tunable plasmonic bandstop filters and slow light devices.
A. Bandstop filtering applications 1) Quadband filter: We utilized the Gr layer's lightblocking phenomenon to obtain the filtering mechanism by enhancing surface light-matter interactions.This blockage of light was the origin of the bandstop plasmonic filter.Based on this mechanism, we utilized four engineered Gr layers with different chemical potentials, providing four transmittance dips due to strong light-matter interaction at resonant frequencies.Fig. 3(a) depicts the transmittance spectra of the proposed quadband plasmonic filter structure.Here, chemical potentials of Gr layers µ c1 = 0.28, µ c2 = 0.20, µ c3 = 0.40, µ c4 = 0.52 eV were adopted.Four transmittance dips formed at the resonant frequencies of 1.746, 2.487, 4.197, and 5.451 THz, respectively.Fig. 3(b) represents the transmitted THz wave in the time domain.The structure experienced a huge deviation from the reference source THz wave because of attenuation by scattering.The transmittance, reflectance, and absorptance spectra are given in Fig. S2(a) of Supplementary Material.To get a better understanding, we presented the spatial electric field distributions of our proposed plasmonic quadband modulator in Figs.3(c)-(f).Each Gr layer primarily localized the electric light field at a resonant frequency.First, second, third, and fourth transmittance dips originated from the second, third, first, and fourth Gr layers.Scattering parameters of the quadband plasmonic modulator exhibited four dips at resonant frequencies, which are provided in Fig. S3(a) of Supplementary Material.2) Pentaband filter: A pentaband plasmonic filter was modeled utilizing multilayered Gr layers with different chemical potentials.Five engineered Gr layers with different chemical potentials were employed to block light by enhancing surface light-matter interactions, as seen in Fig. 1(c).The chemical potentials of five Gr layers were set to be µ c1 = 0.28, µ c2 = 0.18, µ c3 = 0.45, µ c4 = 0.65, and µ c5 = 0.50 eV.Five transmittance dips emerged at the resonant frequencies of 1.404, 2.145, 3.171, 4.197, and 5.337 THz, as shown in Fig. 4. Pentaband filter structure experienced a larger attenuation with prolonged time than quadband filter structure due to scattering among Gr layers and localized light field.The transmittance, reflectance, and absorptance spectra are provided in Fig. S2(b) of Supplementary Material.We studied the spatial electric field distributions of our proposed structure to elevate our knowledge of the localized electric fields which were depicted at resonant frequencies in Figs.4(c)-(g).Scattering parameters of the pentaband plasmonic modulator exhibited five dips at resonant frequencies, as shown in Fig. S3(b) of Supplementary Material.The electric field of incoming light was localized at resonant frequencies by considering each Gr layer as an attenuating medium.First, second, third, fourth, and fifth transmittance dips emerged from the second, third, fourth, first, and fifth Gr layers.3) Plasmon tuning: High tunability without structural modifications for our proposed plasmonic modulators can be achieved by changing surface conductivity via applying V gn .Meanwhile, the surface conductivity varied the µ c of Gr layers [6], [8], [27].The relationship between V gn and ∆µ c is given by [28], Here, ǫ 0 , e, ǫ PMMA , ∆E F , and V Dirac represent the permittivity of vacuum, the charge of an electron, PMMA layer's permit-tivity, change of Fermi level, and the Dirac voltage of Gr layer, respectively.Voltage source varied the change of the Gr layer's chemical potential, as can be seen in Fig. S4 of Supplementary Material.Here, one assumption was made that the impedance of metal contacts is negligible; therefore, it is not considered in our study.A color map of transmittance for our proposed plasmonic modulators displays the structure's wide tunability evolution at different chemical potentials, shown in Fig. 5.
The third and the fourth resonant modes of our proposed quadband and pentaband modulators, respectively, were highly tunable which were consistent with previous report [6].These plasmonic modulators can be operated in the wide range from 1.25 to 5.85 THz by varying the µ c of Gr layers from 0.12 to 0.80 eV.It can be inferred from Fig. 5 that the dips in transmittance are blue-shifted when chemical potentials rise.

4) Impact of polarization angles:
We analyzed the effect of the polarization angles of the incident light on the transmittance spectra of our proposed modulators, as illustrated in Fig. 6.Here, θ represents the polarization angle of incoming light.The electric field of TM polarized incident light (θ = 0 • ) was strongly coupled with the proposed modulators that exhibited transmittance minima at resonant frequencies.Transmittance gradually elevated as θ increased and reached a maximum at θ = 90 • shown in Fig. 6 5) Performance analysis: We enumerated the quality factor (Q) for our proposed quad and penta bandstop plasmonic structures.Q provides the relation between full-width halfmaximum (FWHM) and resonant frequency, f c that is given by,  for the proposed quadband plasmonic modulator.Moreover, we achieved Q of 3.7, 5.38, 5.06, 9.2, and 7.8 at resonant frequencies of 1.404, 2.145, 3.171, 4.197, and 5.337 THz for the proposed pentaband plasmonic modulator.The top GNR layer provided the highest Q compared to Gr nanoholes.Q gives the idea of the filter's selectivity.As Q increases, our filter becomes more selective with decreasing FWHM.Therefore, it can be inferred that the resonance modes at the frequencies of 4.197 and 4.197 THz exhibited high selectiveness for our proposed quadband and pentaband plasmonic filters, respectively.Moreover, we calculated the ER, a performance parameter of the filter, which was solved by, Here, T T M and T T E are the transmittances of TM and TE polarized incident light, respectively.We obtained a high difference between the transmittance under TM and TE polarized light for our proposed quadband and pentaband plasmonic modulators shown in Fig. S5 of Supplementary Material.ERs of 4.6, 12.8, 14.4, and 22.2 dB were achieved at resonant frequencies for the quadband plasmonic filter.Similarly, we found ERs of -0.4,10.0, 16.1, 15.2, and 22.3 dB at resonant frequencies for the pentaband plasmonic filter.Furthermore, we calculated the attenuation of the bandstop filters, a performance parameter of the plasmonic filter, by solving −10log 10 (T T M ).Attenuations of 14.8, 19.1, 16.6, and 22.5 dB were obtained at resonant frequencies for the quadband plasmonic filter.Similarly, we achieved attenuations of 14.8, 22.1, 24.6, 20.3, and 24.5 dB at resonant frequencies for our proposed pentaband plasmonic filter.

B. Slow light applications
We investigated the emergence of slow-light properties, such as group delay and phase, for our proposed THz modulator structures.Fig. 5 demonstrates the evolution of plasmonic resonances with various cases of chemical potentials at different Gr layers.Black dotted lines of Fig. 5 indicated the various cases of quadband and pentaband plasmonic modulators.The chemical potentials of different Gr layers for these cases are provided in Table S1 and S2 of Supplementary Material.The transmittance peaks had a strong destructive interference effect, generating a strong dispersion with a positive group delay.Meanwhile, the phase of our proposed modulators experienced a sharp change while the frequency approached transmission peaks.This phenomenon is depicted in Figs.7  and 8.The group delay, τ g was calculated by, Here, ω and θ p are the angular frequency of light and the phase of the transmittance coefficient, respectively.θ p was enumerated by solving θ p = arg(t) where t denotes the modulator's transmittance coefficient.Fig. 7 illustrates the impact of changing the Gr's chemical potential on slow light properties for our proposed quadband modulator.We can change the τ g and θ p of our proposed quadband plasmonic modulator by varying the µ cn of Gr layers.τ g and θ p were blue-shifted by increasing µ cn of Gr layers, as shown in Fig. 7.
We utilized the µ cn of Gr layers from 0.12 to 0.57 eV.The chemical potentials of case 1 and case 4 were set to be µ c1 = 0.16, µ c2 = 0.12, µ c3 = 0.35, µ c4 = 0.37 and µ c1 = 0.32, µ c2 = 0.28, µ c3 = 0.55, µ c4 = 0.57 eV, respectively.The chemical potentials of the other two cases varied between the chemical potentials of case 1 and case 4. We obtained a maximum group delay of 1.02 ps for the proposed quadband plasmonic modulator.Additionally, we calculated delay-bandwidth product (DBP) by solving DBP = τ g ×∆f .Here, ∆f is the bandwidth of slow light delay spectra.We require slow light devices as a buffer in integrated systems where DBP is an important performance parameter for buffer components.Slow light performance analysis of our proposed quadband plasmonic modulator is enlisted in Table I.We adopted cases from Fig. 5  properties for our proposed pentaband plasmonic modulator was depicted in Fig. 8. Here, we varied the µ cn of Gr layers from 0.12 to 0.80 eV in our study.Increasing µ c of Gr layers provided blue-shifted τ g and θ p , as can be seen in Fig. 8.The chemical potentials of case 1 and case 5 were considered to be µ c1 = 0.24, µ c2 = 0.12, µ c3 = 0.40, µ c4 = 0.60, µ c5 = 0.45 and µ c1 = 0.36, µ c2 = 0.20, µ c3 = 0.60, µ c4 = 0.80, µ c5 = 0.57 eV, respectively.The µ cn of the other three cases were adopted between the µ cn of these two cases.We achieved a maximum group delay of 1.41 ps for the proposed pentaband plasmonic modulator, which was the Gr modulator's highest reported group delay.The group delay and phase were blue-shifted by increasing the Gr layers' µ cn .Thus, we can alter the slow-light properties of our proposed modulators by varying the µ cn of Gr layers, which allows us to create external voltage-controlled tunable slow-light devices.performance analysis of our proposed pentaband plasmonic modulator.We adopted cases from Fig. 5(b) where dot lines indicated the five cases of our proposed pentaband plasmonic modulator.A calculated maximum DBP of 1.725 was achieved for our proposed pentaband plasmonic modulator.

IV. COMPARATIVE ANALYSIS
Table III shows a detailed comparative analysis of the performance parameters among the previously reported studies and our proposed modulators.Our obtained τ g outperformed In comparison to existing Gr-based devices, our suggested quadband and pentaband plasmonic modulator structures offer higher ER with a τ g of 1.41 ps which is quite large.Dual-frequency on-off modulator 2 -7 11.55 0.239 [30] Patterned monolayer graphene metamaterial 4 -7 -0.255 [31] Single-layer graphene ribbon and strips 2.5 -5.75 11.75 0.623 [13] Rectangular defect in graphene 2.5 -5.5 -0.363 [32] Dual-Mode On-to-Off

V. CONCLUSION
In this paper, we proposed two multilayered Gr-based quadband and pentaband plasmonic modulators and utilized the FDTD technique to conduct the performance analysis of these modulators.The Gr layers strongly interacted with the incident light's electric field, resulting in SP resonance at resonant frequencies due to scattering.Varying applied voltage changed the Gr layer's chemical potential, allowing the modulators to tune the resonant frequencies without modifying structural parameters.Meanwhile, the change in chemical potentials shifted the Gr's surface conductivity.The proposed modulators exhibited high tunability over a wide operating range of THz spectrum.We obtained a noteworthy improvement in performance for the pentaband plasmonic modulator compared to the quadband plasmonic modulator.Our designed pentaband plasmonic modulator had more than a quarter group delay compared to the quadband plasmonic modulator; however, the quadband plasmonic modulator provided a 25% higher Q-factor.We achieved a maximum optical group delay of 1.41 ps and ER of 22.3 dB for our proposed pentaband plasmonic modulator, which is significantly higher than previously reported works.Therefore, the proposed highperformance modulators will benefit optical filtering and slow light applications, specifically in optical buffers.

Fig. 1 .
Fig. 1.(a) Illustration of a quadband plasmonic modulator unit cell.Here, four Gr monolayers were employed, and the structural parameters were: i = j = k = 150 nm, h = 2 µm.The output transmission signal attenuated and deviated from the source signal due to scattering.(b) Top Gr layer of the proposed structures where the period of the structure, p = 3 µm, length of GNR, q = 2.5 µm.Inset shows the magnified view of GNR.Here, the pitch distance of GNR, d is 20 nm, width, w = 10 nm.(c) Schematic view of a pentaband plasmonic modulator unit cell with geometrical parameters: i = j = k = l = 150 nm.Vg n is the applied voltage source of the n-th Gr layer via a metal contact.The Vg n varied the µc n of Gr layers and σg .

Fig. 2 .
Fig. 2. Evolution of (a) Re (σg ) and (b) Im (σg ) of Gr layer at different Fermi levels.(c) Re (σg ) and (d) Im (σg ) parts of the Gr layer at the E F of 0.6 eV.Gr exhibited more metallic behavior with increasing Fermi energy levels.

Fig. 3 .
Fig. 3. (a) Transmittance spectra of the proposed quadband plasmonic filter.Here, µc 1 = 0.28, µc 2 = 0.20, µc 3 = 0.40, and µc 4 = 0.52 eV were employed.(b) The E-field time-domain waveforms of the proposed quadband plasmonic filter at the normal incidence for TM polarized incident THz light.The transmitted waveform attenuated and deviated from the source THz waveform.The spatial E-field distribution of Gr layers at resonant frequencies of (c) 1.746,(d) 2.487, (e) 4.197, and (f) 5.451 THz for our proposed quadband plasmonic filter.The top, second, third, and bottom Gr layers showed localized E-fields at the resonant frequencies of 4.197, 1.746, 2.487, and 3.385 THz, respectively.The logarithmic scale was utilized to represent the color map of the E-field.

Fig. 5 .
Fig. 5. Evolution of plasmonic (a) quadband, and (b) pentaband transmittance dips at different chemical potentials of Gr layers.We applied a TM-polarized normal incidence plane wave THz light on our proposed modulators.Transmittance dips and peaks were blue-shifted by increasing the chemical potentials of Gr layers.Dot lines represent the different case studies of our proposed modulators.
. The symmetric nature of the proposed modulators provided the same transmittance spectra at θ = 0 • = 180 • and θ = 90 • = 270 • .Lower-frequency resonances were less sensitive to polarization angles compared to high-frequency resonances.It can be observed that highfrequency resonances were confined to the polarization angle range from 45 • to 315 • and 135 • to 225 • .Two high-frequency resonances existed around 4.197, 5.451 THz, and 4.197, 5.337 THz for our proposed quadband and pentaband plasmonic modulators.

Fig. 7 .
Fig. 7. Group delay and phase spectra of the quadband plasmonic modulator at (a) case 1, (b) case 2, (c) case 3, and (d) case 4. Black dotted lines of Fig. 5(a) indicate the case studies of our proposed modulator.
(a) where black dotted lines indicated the four cases of our proposed quadband plasmonic modulator.DBP can be adjusted without changing structural parameters by varying µ c via an external voltage source.A maximum DBP of 1.268 was obtained for our proposed quadband plasmonic modulator.The impact of varying the chemical potential of Gr layers on slow light

TABLE I CALCULATED
τg AND DBP FOR VARIOUS CASE STUDIES OF QUADBAND

TABLE II CALCULATED
[13]ND DBP FOR VARIOUS CASE STUDIES OF PENTABAND Xiong et al. reported a continuous staircase-shaped structure on the Gr monolayer[14].Their proposed structure provided a high ER of 16.56 dB with a τ g of 0.488 ps.However, a staircase shape in the Gr layer makes the fabrication procedure complex.Gao et al. proposed a modulator comprising four GNRs on silicon substrate[13].Low ER was the limiting factor of their modulator structure.
[29]et al. demonstrated a similar kind of GNR modulator which had similar disadvantages as previous work including low ER and small operating range[29].Hu et al. reported an asymmetrical Gr-based structure; however, the operating range of the structure covered a small region of the THz spectrum.

TABLE III COMPARATIVE
PERFORMANCE ANALYSIS OF GR STRUCTURES