A dynamically reconfigurable Fano metamaterial through graphene tuning for switching and sensing applications

We report on a novel electrically tunable hybrid graphene-gold Fano resonator. The proposed metamaterial consists of a square graphene patch and a square gold frame. The destructive interference between the narrow- and broadband dipolar surface plasmons, which are induced respectively on the surfaces of the graphene patch and the gold frame, leads to the plasmonic equivalent of electromagnetically induced transparency (EIT). The response of the metamaterial is polarization independent due to the symmetry of the structure and its spectral features are shown to be highly controllable by changing a gate voltage applied to the graphene patch. Additionally, effective group index of the device is retrieved and is found to be very high within the EIT window suggesting its potential use in slow light applications. Potential outcomes such as high sensing ability and switching at terahertz frequencies are demonstrated through numerical simulations with realistic parameters.

We report on a novel electrically tunable hybrid graphene-gold Fano resonator. The proposed metamaterial consists of a square graphene patch and a square gold frame. The destructive interference between the narrow-and broadband dipolar surface plasmons, which are induced respectively on the surfaces of the graphene patch and the gold frame, leads to the plasmonic equivalent of electromagnetically induced transparency (EIT). The response of the metamaterial is polarization independent due to the symmetry of the structure and its spectral features are shown to be highly controllable by changing a gate voltage applied to the graphene patch. Additionally, effective group index of the device is retrieved and is found to be very high within the EIT window suggesting its potential use in slow light applications. Potential outcomes such as high sensing ability and switching at terahertz frequencies are demonstrated through numerical simulations with realistic parameters.
T he idea of electromagnetically induced transparency (EIT) and the concept of Fano resonance were originally discovered in the context of quantum mechanics 1 . EIT is generated when a narrowband discrete state destructively interferes with a broadband continuum; the resulting spectrum has the Fano line-shape 1 . Plasmonic analogues of EIT have been recently generated via coupling either antiparallel dipolar surface plasmon polaritons (SPPs) or bright dipolar SPPs with dark higher order ones 2,3 , all induced on metal-dielectric interfaces. The excitation of the originally dark higher order SPPs is often made possible by means of symmetry breaking [3][4][5] . Generation of physical phenomena equivalent to EIT at terahertz (THz) frequencies is currently one of the most exciting topics of plasmonics research 2 . Since most of the biomolecules have their vibrational modes oscillating at THz frequencies 6,7 , this means that EIT generated at these frequencies offers a viable way for biosensing. The drastic ''slowing'' of light around the narrow Fano resonance results in an increased sensitivity to changes in the medium's refractive index. Plasmonic EIT could also be used in designing efficient switches for modulating the amplitude and phase of waves transmitted through metamaterials 8 . Possibility of designing ultra-sensitive bio-detectors and efficient switches urges the investigation of plasmonic metamaterials capable of supporting Fano resonances at THz frequencies. The unique and highly tunable electrical properties of graphene 9 observed within this band of the spectrum render it an attractive candidate as a building block of such metamaterials.
Since its first practical isolation by Novoselov and Geim in 2004 9 , graphene solicited a keen interest among physicists and engineers. This could be attributed to graphene's unprecedented properties that cannot be found in any other material 10 . High electron mobility 11 , significant white light absorption 8,12,13 , ability to support SPPs [14][15][16][17] are among these characteristics that are relevant in the fields of electromagnetics, optics, and photonics 18 .
In this work, graphene's ability to support SPPs at THz frequencies is exploited to design a Fano resonator that hybridizes SPPs generated on graphene and gold surfaces. Graphene surface plasmon polaritons (GSPPs) have several advantages when compared to SPPs generated on metallic surfaces: GSPPs have higher volume confinement (exceeding 10 6 times the diffraction limit), are easier to tune (obtained via applying a gate voltage to the graphene), and propagate longer distances and have narrower spectral support (due to the lower intrinsic losses in graphene) 19,20 . Not surprisingly, these superior features of GSPPs have fueled research in several directions such as development of GSPP waveguides 21,22 , THz antennas [23][24][25] , perfect absorbers 26,27 , novel Fourier optics devices 28 , THz cloaks 29,30 , photonic crystal nano-cavities 31 , modulators 32 and sensors 33 .
The hybrid graphene-gold Fano resonator proposed in this work is a doubly periodic array of a unit cell constructed using a square graphene patch located at the center of a square gold frame. The Fano resonance is obtained from the destructive interference between the dipolar SPPs generated on the graphene patch and the gold frame at THz frequencies and its spectral shape and location can easily be tuned by applying a gate voltage to the graphene patch.

Results
Physical mechanism to generate Fano resonance. The Fano resonator proposed in this work is a doubly periodic array of a unit cell that consists of a square graphene patch located at the center of a square gold frame. The resonator is embedded in a dielectric substrate. (see Fig. 1). The proposed Fano resonator is excited with a plane wave propagating in the z-direction with electric field polarized in the xdirection [Figs. 1(b) and 1(c)]. It should be noted here that due to the symmetry of the unit cell and the same periodicity in x-and y-directions, the response of the resonator is polarization independent. When excited with this incident field separately, the gold frame and the graphene patch support dipolar SPPs. As a result, the resonator constructed using the gold frame and the graphene patch supports an SPP mode hybridized between these two dipolar modes. Since the dipolar SPP induced on the frame has a much broader spectral support than the one induced on the graphene patch due to the presence of higher intrinsic losses in gold, these two modes' destructive interference generates a Fano line-shape in the response of the resonator. It should be noted here that the plasma frequencies of gold and graphene are separated by a large offset. This means that the resonance frequencies corresponding to the dipolar SPPs induced on the gold frame and the graphene patch are expected to be separated by a large gap. To move the resonance frequencies closer to each other and have the SPPs' spectral supports overlap at THz frequencies, the geometrical dimensions of the frame and the patch are chosen to be in mm scale as shown in Fig. 1(a). Obviously, in this design, SPP of the graphene patch acts like the ''dark mode'' of the traditional Fano resonator designs made of only metals 3 . Additionally, unlike these traditional resonators, the proposed design does not require its symmetry to be ''broken'' since the narrower mode can be directly excited by the incident field even if the structure is fully symmetric. The spectral location and line-shape of the graphene's narrow dipolar SPP is determined by graphene's complex relative permittivity e V,G , which can be controlled by adjusting the chemical potential m c (see Methods Section on Material Models for details). Variation in m c can be achieved by applying a gate voltage to the graphene patch using nearly transparent electrodes without perturbing the response of the resonator as suggested recently in reference 8 .
Proof of concept via numerical experiments. The hybridization of the dipolar SPPs induced on the gold frame and the graphene patch is demonstrated by an example. For this simulation, the dimension of the gold frame S 1 5 5.5 mm, the dimension of the graphene patch S 2 5 1.6 mm, the relative permittivity of the substrate e d 5 3.5, and graphene's electron mobility m 5 10,000 cm 2 /Vs and chemical potential m c 5 1500 meV. It should be noted here that the value assigned to m is a rather conservative choice considering the latest experimental results 18 . The transmittance of three structures constructed using unit cells with only the gold frame, only the graphene patch, and both the gold frame and the graphene patch are computed [Figs. 2(a), 2(b), and 2(c), respectively]. As expected, the dipolar SPPs induced on the gold frame (marked as D 1 ) and the graphene patch (marked as D 2 ) have broad (continuum-like) and very narrow spectral supports, respectively. Figure 2(c) clearly demonstrates the asymmetric Fano-like spectral line-shape and a narrow EIT window in the response of the resonator constructed using both the gold frame and the graphene patch, which results from destructive interference of D 1 and D 2 .
The dipolar nature of the SPPs is exhibited in Fig. 3 showing the surface charge distributions and the magnetic field norms computed at several frequency points. Fields due to D 1 are clearly more dominant at point I, which is far away from the resonance frequency of D 2 . Surface charge distributions at the frequency points II, III, and IV (around the Fano resonance), clearly demonstrate that D 2 interferes   with D 1 . Additionally, comparison of charge distributions (and magnetic field norms) at points II and III reveals that the sharp roll-off of the transmittance between the EIT point at 8.35 THz and the graphene's dipolar SPP at 9 THz is due to the phase change of the field distribution at points II and III.
The effect of geometrical dimensions on the response of the Fano resonator is characterized next. For the first set of simulations, S 1 5 5.5 mm, e d 5 3.5, m 5 10,000 cm 2 /Vs, and m c 5 1500 meV, while S 2 is varied between 1 mm and 2 mm. Transmittance of the resonator is computed for each value of S 2 [ Fig. 4(a)]. Increasing S 2 (for example by setting it to 2 mm) red-shifts the resonance of D 2 as it increases the graphene patch's effective dipole length. This moves the resonant frequencies of D 1 and D 2 closer to each other and forces them to couple more strongly. As a result, a higher increase in the transmittance within the EIT band [with respect to the transmittance of only D 1 plotted in Fig. 4(a) as thin blue line] is observed. On the other hand, when S 2 is decreased (for example by setting it to 1 mm), the resonance frequencies of D 1 and D 2 get sufficiently detuned. In this case, transmittance follows very closely the transmittance of only D 1 with a much smaller relative increase in the EIT band. For the second set of simulations, S 2 5 1.6 mm, e d 5 3.5, m 5 10,000 cm 2 /Vs, and m c 5 1500 meV, while S 1 is varied between 4 mm and 6 mm. Transmittance of the resonator is computed for each value of S 1 [ Fig. 4(b)]. As expected, increasing S 1 red-shifts the resonance of D 1 as it increases the frame's effective dipole length. For S 1 5 6 mm, a higher increase in the transmittance within the EIT band [with respect to the transmittance of only D 1 plotted in Fig. 4(b) as thin blue line] is observed. For S 1 between 4.5 mm and 4 mm, the asymmetric Fano line-shape appears on the left side of the resonance frequency of D 1 . These results clearly indicate that depending on the target application, the dimensions of the graphene patch and the gold frame can be fine-tuned. For example, if one wants to design an optical modulator using the Fano resonator proposed here, a greater modulation depth can be achieved by choosing the geometry dimensions that result in a large variation in the transmittance within the EIT band.
The effect of graphene's intrinsic loss, i.e., the value of its electron mobility, m, on the response of the Fano resonator is also analyzed. For this set of simulations, S 1 5 5.5 mm, S 2 5 1.6 mm, e d 5 3.5, and m c 5 1500 meV, while m is swept from 1,000 cm 2 /Vs to 250,000 cm 2 / Vs. Figure 4(c) shows the transmittance of the resonator computed for each value of m. For lower values of m, i.e., for higher loss, the absorption channel becomes dominant and the transmission efficiency is deteriorated as expected. For realistic values of m around 10,000 cm 2 /Vs, the response of the Fano resonator is good; the transmittance efficiency reaches almost 60%.
Finally, the effect of the substrate on the response of the resonator is characterized. For this set of simulations, S 1 5 5.5 mm, S 2 5 1.6 mm, m 5 10,000 cm 2 /Vs, and m c 5 1500 meV, while e d is varied between 1 (no substrate) and 5. The transmittance of the resonator is computed for each value of e d [ Fig. 4(d)]. It is clearly shown in the figure that a gradual blue-shift is observed in the response of the resonator as e d is increased. A gradual decrease in the maximum amplitude in transmittance is also observed with an increase in e d .
Applications. Switching. The EIT window generated via the destructive interference of the dipolar SPPs induced on the gold frame and the graphene patch is tuned by controlling a gate voltage applied to the graphene patch (as discussed for example in recent studies 8 ). To simulate the effect of the variation in the gate voltage, graphene's chemical potential m c is swept between 500 meV and 1500 meV while e d 5 3.5 and m 5 10,000 cm 2 /Vs.
The transmittance of the structure constructed using unit cells with only the graphene patch with S 2 5 1.6 mm is computed for various values of m c [ Fig. 2 (b)]. The figure clearly demonstrates that increasing m c leads to a blue shift in the resonance frequency of D 2 and comparatively stronger extinction amplitude 1 2 T (where T denotes the transmittance) at the resonance frequency.
Similarly, transmittance and the phase of transmission of the structure constructed using unit cells with the gold frame with S 1 5 5.5 mm and the graphene patch with S 2 5 1.6 mm is computed for the same values of m c [Figs. 5(a) and 5(b)]. The figures clearly show that the spectral location and the shape of the EIT window can be tuned by varying m c . This easily tunable EIT window can be utilized as a mechanism to make switches at THz frequencies. For example, as shown in Fig. 6(a), the transmittance at 7.94 THz can be switched between 0.4% to 53% by simply varying the chemical potential m c by an amount of Dm c 5 (1430 2 1280) meV 5 150 meV. As a result a maximum modulation index of 0.52 is achieved through a small variation of Dm c 5 150 meV. The voltage-controlled resonator has in addition the potential to be used as a phase modulation device. The highly dispersive propagation within the EIT band has a steep spectral variation in transmission phase w [ Fig. 2(c)]. This feature allows to modulate the phase of the transmitted signal to a substantial degree. As shown in Fig. 6(b), the phase of the transmission at 8.22 THz can be changed by Dw 5 0.68 rad, by simply varying the chemical potential m c by the same amount as for the amplitude. The amplitude and phase modulation index can be further improved by optimization of spectral positions for the interference of the resonances. For example, maximum modulation index for the amplitude can be increased to 0.65 by setting S 1 5 6 mm [ Fig. 4(b)]. The commercial applications for this integrated THz modulator are in the areas of communication systems, high speed Mach-Zehnder modulators, phase array antennas and time-domain spectroscopy 35 .
Slow light and sensing. The effective refractive and group indices, n e and n g , of the resonator design with S 1 5 5.5 mm, S 2 5 1.6 mm, e d 5 3.5, m 5 10,000 cm 2 /Vs, and m c 5 1500 meV are extracted. To this end, a homogenous slab with refractive index n e , which generates the same S-parameters as the proposed design, is found using the retrieval method described in reference 36 . Then, the effective group index is computed using the relation n g (v) 5 n e (v) 1 vhn e (v)/hv 36 . The retrieved n e and n g are shown in Figs. 7(a) and 7(b), respectively; the value of n g exceeds 1,400 within the EIT window. This value is  much higher than that of the other plasmonic resonators solely made of metals 2 . High values of n g clearly demonstrate that the proposed Fano resonator design has the potential to be used in slow light applications including ultra-sensitive biomolecule detection.

Discussion
To the best of our knowledge, our work is the first to report on EIT generated at THz frequencies using a hybrid graphene-gold structure. This novel design makes use of the fact that dipolar SPP induced on graphene are much narrower than those induced on gold to induce the asymmetric Fano resonance shape in the spectrum. This concept is novel and has not been exploited before. The use of this idea equipped the proposed Fano resonator with the following properties that are superior to ''classical'' Fano resonators constructed using only noble metals: 1) Polarization Independence: The design is polarization independent since it does not require excitation of originally dark modes via symmetry breaking. 2) Dynamic Tunability: The spectral location and shape of the Fano resonance (and the EIT window) can be tuned by varying the chemical potential of the graphene patch. This can be dynamically achieved by applying a gate voltage to the graphene patch. 3) Higher Group Index: Group index extracted around the EIT window is higher than that reported in literature previously. This equips our design with a high sensitivity to the changes in the background medium's refractive index.
It should be emphasized here that because gold's intrinsic loss is higher than graphene's, one would expect the overall loss in the hybrid resonator would increase destroying the EIT. But it is observed that within the EIT region, where the destructive interference occurs, the radiation losses are suppressed. This is indeed why adding the gold frame to the system decreases the total losses within the EIT region instead of increasing them. The high transmittance values, which reach 80% as could be seen in Fig. 4(c), fully demonstrate this fact.
The proposed resonator has potential applications in designing efficient switches and ultra-sensitive bio-detectors that can be operated at THz frequencies. Our work also demonstrates the possibility of these via numerical simulations with realistic parameters.

Methods
Material models. The complex dielectric constant of gold is accurately modeled at THz frequencies using the Drude model with plasma frequency v p 5 1.37 3 10 16 rad/s and damping constant c c 5 39.47 3 10 12 rad/s 34 . The complex surface conductivity s S,G for a graphene layer is calculated from Kubo's formula 8,14 : Here, s intra and s inter represent the intra-and inter-band transitions in the graphene layer and their expressions could be found in many recent studies 19,20 . At low THz frequencies, where we have s intra ? s inter , s S,G could be approximated by a Drude where v is the angular frequency, q is the electron charge, is the reduced Planck constant, m c is the chemical potential of the doped graphene layer, and c is the damping constant. In this work, it is assumed that the thickness of the graphene layer d 5 1 nm; this choice was mainly motivated and justified by the seminal work of Vakil and Engheta 19 . Since d is much smaller than the wavelength at THz frequencies, graphene's bulk conductivity can be simplified as s V,G 5 s S,G /d. Using Ampere's law in stationary regime and Ohm's law, one can calculate complex relative bulk permittivity of graphene as    dependence suggests that m c can be varied to tune the spectral location and line-shape of the Fano resonance.
These material models are used in the finite element program COMSOL Multiphysics to carry out the simulations of electromagnetic wave interactions on the design presented in Figs. 1(a), 1(b), and 1(c). The results of these simulations are presented in the previous sections. These simulation results can also be ''predicted'' by the response of an RLC circuit as described in the next section.
RLC analytic model. The optical response of the proposed Fano resonator can be mathematically modeled using coupled oscillator equations. In this work, an RLC circuit model is used to replicate the optical response of the Fano resonator 37 . In the circuit model, each of the dipolar SPPs induced on the gold frame and the graphene patch are represented by an RLC loop as shown in Fig. 8. Note that the additional third RLC loop in the circuit is needed to take into account the coupling from higher order SPPs induced on the gold frame and graphene patch beyond 12 THz. The loops are connected through three capacitors C 12 , C 13 , and C 23 , which model the coupling between the different SPPs. This RLC circuit can be mathematically represented as a system of three coupled equations: Here, q i , i 5 1,2,3, is the charge due to the steady state current I i 5 jvq i flowing in the loop indexed with i. The term V ij 5 (L i C ij ) 20.5 , i,j 5 1,2,3, i , j, couples the three equations to each other, L i~{ v 2 zjvc i zv 2 i represents the ''self-coupling'', where 20 13 zC {1 23 . The terms Q i 5 V i /L i , i 5 1,2, on the right side of the system of equations (2) represent the excitation, i.e., the direct coupling from the incident field to the SPPs. The time average power drawn by the circuit is P 5 1/2Re[V 1 I 1 * 1 V 2 I 2 *] and this quantity should be equivalent to the extinction coefficient 1 2 T of the Fano resonator. By minimizing the difference between P and 1 2 T numerically, v i , c i , Q i , and V ij can be found. Once these coefficients are known, one can easily obtain the physical parameters, R i , L i , C i , i 5 1,2,3 and C 12 , C 13 , and C 23 , which describe the RLC circuit.
The mathematical model described by the system of equations (2) provides additional physical insight into the response of the Fano resonator. For example, Q i , i 5 1,2, represent the amount of relative power each resonant mode receives from the incident electromagnetic field. Non-zero values of Q i indicate that both modes have dipole moments along the polarization vector of the incident field and they can be directly excited. Variables V 12 , V 13 , and V 23 represent the amount of (energy) coupling between the modes. For example a high V 12 means that the coupling between D 1 and D 2 is strong, which indicates high transmittance T within the EIT region. This is demonstrated by an example as described next.
Parameters of the circuit with three loops is obtained by minimizing the difference between P and 1 2 T of the Fano resonator with S 1 5 5.5 mm, S 2 5 1.6 mm, e d 5 3.5, and m 5 10,000 cm 2 /Vs for various values of m c between 500 meV and 1500 meV. The values of the extracted parameters are provided for each value of m c in Table I. The effect of the increase in m c can be seen with a blue shift in the resonance frequency v 2 that corresponds to D 2 . Additionally, a consistent increase in V 12 , which represents the strength of coupling between D 1 and D 2 , is observed. As mentioned above, the increase in coupling strength, i.e., increase in V 12 translates as increased T within the EIT region [see Fig. 5(a)]. Note that as expected, increasing m c does not effect characteristics of D 1 . As shown in Table I, the resonant and damping frequencies, v 1 and c 1 , corresponding to D 1 remain unchanged.
Table I | Parameters of the RLC circuit model representing the Fano resonator with S 1 5 5.5 mm, S 2 5 1.6 mm, e d 5 3.5, and m 5 10,000 cm 2 /Vs for various values of m c between 500 meV and 1500 meV. Units of v 1 , v 2 , v 3 , V 12 , V 13 , and V 23 are (meV) while the units of Q 1 and Q 2 are (mV/H)