Abstract

A multiple plasmon-induced transparency (PIT) device operated in the mid-infrared region has been proposed. The designed model is comprised of one graphene ribbon as main waveguide and two narrow graphene sheets resonators. The phase coupling between two graphene resonators has been investigated. The multimode PIT resonances have been found in both cases and can be dynamically tuned via varying the chemical potential of graphene resonators without optimizing its geometric parameters. In addition, this structure can get multiple PIT effect by equipping extra two sheets on the symmetric positions of graphene waveguide. The simulation results based on finite element method (FEM) are in good agreement with the resonance theory. This work may pave new way for graphene-based thermal plasmonic devices applications.

1. Introduction

Surface plasmon polaritons (SPPs), coupled modes of plasmons and photons, propagate along the metal-dielectric interface [1, 2]. SPPs provide an intriguing approach to localize and guide light in subwavelength metallic structures. In recent years, a great variety of efforts have been devoted to applications based on SPPs, which promotes the development for high density integration of photonic circuits. Electromagnetically induced transparency (EIT) is a sharp transparency window within a broad absorption spectrum in the medium, which is due to the coherent interaction between atomic upper levels and applied optical field [3, 4]. Recently, a novel analogy phenomenon to EIT, known as plasmonically induced transparency (PIT), has been manifested in various systems based on metamaterial structures [512]. The PIT effect has been applied in various applications such as nonlinear optical processes, ultrafast switching, slow light, and signal processing [13, 14]. However, the most of PIT effects are based on MIM (metal-insulator-metal) structure, and the adjustment of PIT effects must rely on changing the geometrical parameter of structures. Recently, multimode PIT effects have been achieved by dual coupled graphene-integrated ring resonators [15, 16]. The MIM was replaced by graphene structure to achieve PIT effects without modifying parameter of structure. The dynamically tunable PIT planar hybrid metamaterials have been proposed, in which the VO2 stripes are filled in the grooves of the device surface. The simulation results indicate that the PIT effects can be achieved by changing the temperature of VO2 stripes [17]. To compare with [18], the structure of our proposal is simpler and more convenient in modulating PIT effects.

In recent years, graphene has aroused extensive attention for its unique and fantastic properties in optics and electronics [19]. A monolayer graphene has been widely researched as a promising platform for plasmonics for well supporting SPPs on graphene [20, 21]. Most importantly, the surface conductivity of graphene could be dynamically tuned by chemical potential via gate voltage, electric field, magnetic field, and chemical doping [22], which increasingly promotes the development of active plasmonic devices including polarizers [23], high-speed optical modulators [24], absorbers [25], and transformation optical devices [26]. Zeng et al. have studied a graphene-integrated Fabry-Perot (F-P) microcavity to efficiently modulate plasmon-induced transparency effect [27]. Shi et al. numerically demonstrated a plasmonic analog to EIT based on graphene nanostructures [28]. Therefore, the monolayer graphene is a promising platform to design a dynamically wavelength tunable PIT planar device.

In this paper, we have proposed and investigated a PIT device operating in the mid-infrared region. The designed model consisted of one graphene ribbon as main waveguide and two narrow graphene sheets as resonators. Firstly, a graphene waveguide side-coupled with single graphene resonator has been studied as a filter. Next, the phase coupling between two graphene resonators side-coupled to the main graphene waveguide has been investigated. In order to modulate multipeak PIT, four sheets were inserted into the symmetric position of graphene waveguide. The multimode and multipeak PIT resonances can be actively controlled via changing the chemical potential (Fermi energy) of graphene. Our designs have been well verified by numerical simulation based on finite element method (FEM). Hence, the proposed structures may have great potential in ultra-compact graphene thermal plasmonic devices in the infrared range.

2. Method

The structure of proposed device is schematically depicted in Figure 1. It consists of one graphene ribbon as main waveguide and two uniform graphene sheets as resonators. For experimental consideration, the designed graphene system is assumed in the SiO2 medium. It should be noted that the SiO2 medium is replaced with air in simulation for simplicity. The Kubo formula has governed the optical conductivity of graphene including the interband and intraband transition contributions [29]. It depends on the momentum relaxation time , temperature , incident wavelength (angular frequency ), and chemical potential (Fermi energy ) . The surface conductivity of graphene follows the equation

In our simulation, the employed incident light is in the mid-infrared region where the intraband transition contribution dominates in monolayer graphene [30]. Under this condition, the optical conductivity is given byHere, is the electron charge and the carrier relaxation time relates to the carrier mobility , chemical potential , and Fermi velocity  m/s in graphene. The equivalent permittivity of graphene is given by the equation [26] where is wave number in vacuum and  Ω is the intrinsic impedance of air. The carrier mobility is reasonably chosen to be  cm2 V−1 s−1 from experiment results [25, 31]. The thickness of graphene is assumed as  nm in order to take the numerical simulation, but is not the real thickness of the graphene [19, 21]. TM- (Transverse Magnetic-) polarized SPPs on graphene are only in consideration for the investigation. The dispersion relation of this TM SPP surface wave follows the equationwhere is the propagation constant of graphene SPPs in air (). The effective refractive index of graphene SPPs shows the ability to confine SPPs on graphene, which is defined as . The propagation length that indicates the SPP propagation loss in graphene. The dependences of and on the Fermi level and incident light wavelength are shown in Figure 2. Obviously, from Figure 2(a), increases as the Fermi level decreases for a fixed wavelength, which means that SPPs are better confined when Fermi level is lower. Nevertheless, the tendency in Figure 2(b) is evidently opposite to that in Figure 2(a), indicating that a lower Fermi level gives a smaller propagation length. Thus, both these important factors should be taken into consideration in the design of PIT systems. Interestingly, varies greatly when the Fermi level is changed slightly, which provides a way to actively control our configuration.

3. Result and Discussion

At first, we focus on one single graphene resonator side-coupled to graphene waveguide as a filter. As shown in Figure 1, there are two graphene sheet resonators whereas we set the right resonator inactive by fixing its at 0.05 eV. Under this condition SPPs will not be supported on graphene because the interband transition of electrons occurs at this Fermi energy [11]. The Fermi level of graphene waveguide is fixed at 0.4 eV. When SPPs on graphene meets the resonant condition, the graphene resonator will influence the transmission of main waveguide by forming standing wave. The resonance wavelength is determined by the width and Fermi level of the graphene resonator, which can be deduced from quasi-static analysis and Drude-like model [19, 32]. It follows the following equation:Here, is the velocity of light in vacuum and parameter is introduced for a freestanding graphene sheet () [33]. The fitting dimensionless parameter can be derived from numerical simulation [19, 31, 34] and is the width of graphene sheet. The resonance condition is given by [34]. Here, is an arbitrary positive integer, resonance mode number, and is an additional phase shift because of the edge ends of the graphene sheets. Actually, the influence of the edge ends is only effective when the structure is in extremely small size, so we can neglect the influence in most cases [33, 35]. Besides, another equation for propagation constant of SPPs on monolayer graphene has been derived by combining the Drude-like model with the dispersion relation [36]:

By using the FEM, the transmission spectra have been calculated for different Fermi level of left graphene resonator ranging from 3.6 μm to 5 μm, shown in Figure 3(a). The above discussion indicates that SPPs are better confined when Fermi level is lower. So the Fermi level is selected to be 0.42 eV, 0.43 eV, 0.44 eV, 0.45 eV, and 0.46 eV, respectively. The corresponding transmission spectra are denoted by blue, pink, green, purple, and red lines, respectively. As the Fermi level decreases, an obvious red shift of the transmission dip can be found. Obviously, the larger Fermi levels correspond to the higher transmission dips at the same resonance mode, shown in Figure 3(a). It is because that the loss of propagation mode determined by the imaginary part of is inversely proportional to the Fermi level. As shown in Figures 3(b)3(d), when SPPs on graphene meets the resonant condition, the energy can be confined in graphene sheet by forming a standing wave mode, which leads transmission dips. According to the above analysis, this structure can achieve filter effects by controlling Fermi level without optimizing geometric parameters.

Next, the influence of interval has been investigated. The transmission spectra have been calculated for different values of (45 nm, 50 nm, 55 nm, 60 nm, and 65 nm), shown in Figure 4. The transmission spectra of selected to be 45, 50, 55, 60, and 65 nm are represented by blue, pink, green, purple, and red lines, respectively. Obviously, as the interval increases, the coupling process becomes weaker. In order to investigate the coupling between main waveguide and graphene sheet concretely, the coupled mode theory is introduced to interpret the simulation results. By using coupled mode theory, the transmission can be expressed as [37]. Here, and are the incident frequency and the resonance frequency, respectively. is the decay rate resulting from the intrinsic loss of graphene sheet resonator, and is the decay rate due to the power escaping from the graphene waveguide. With the growing of interval , the decay rate decreases while the decay rate remains the same. In this way, the energy confined in graphene resonator will become less for an increasing interval , which leads to a growing value of the transmission dip. At the same time, the spectra width is proportional to the decay rate . That explains that the spectrum width becomes narrower with a growing interval at one resonance mode. This structure can achieve best filter effect when interval is fixed at 50 nm.

In this section, the discussion focuses on the phase coupling between two graphene sheet resonators. To avoid the direct interference between two graphene sheets, the separation has been chosen to be large enough. The other geometric parameters are selected as follows:  nm and  nm. We know that the resonance wavelength can be tuned via gate voltage on graphene sheets through the resonant condition. The difference of resonant wavelength between two graphene resonators is defined as (, , for each mode). The of left graphene sheet is fixed at 0.42 eV, which presents three resonance dips (4.655 μm, 4.165 μm, 3.8 μm) corresponding to first-order, second-order, and third-order resonance modes depicted in Figure 3(a). The gradually increasing leads to a broader gap between two resonance dips for three resonance modes as of the right sheet grows from 0.425 eV to 0.44 eV. In Figures 5(a), 5(b), and 5(c), for each resonance mode, all PIT windows show obvious off-to-on feature and the transparency bands become broader with the growing. The magnetic distributions corresponding to the transparency peaks of 1st, 2nd, and 3rd resonance modes are depicted in Figures 5(e), 5(f), and 5(g). In this way, the tunable PIT effects based on graphene sheet resonators are achieved.

In this section, the phase coupling between the four graphene sheets structure has been investigated by contrasting with the two grapheme sheets structure. The extra graphene sheets are placed to graphene waveguide symmetric position on the basis of two sheets structure, which is shown in insert picture of Figure 6(d). The transmission spectra of second-order mode and third-order mode are depicted in Figures 6(a), 6(b), 6(c), and 6(d); the modulation of multipeak PIT can be achieved by controlling the Fermi level of four graphene sheets. By adjusting Fermi level gradually can result in tunable off-to-on multipeak PIT feature. Figure 7 shows the transmission spectra in the mid-infrared region of two peaks PIT based on the four graphene sheets structure. The upper and lower picture of Figure 7 show the magnetic distributions of the second-order mode under μm, 4.12 μm, and 4.14 μm, which correspond to the standing wave modes originating from the splitting of the second-order resonance. In addition, Figure 8 shows the transmission spectra of three peaks PIT based on the four graphene sheets structure. The upper and lower picture of Figure 8 show the magnetic distributions of the second-order mode under μm, 4.1 μm, and 4.12 μm. Due to the great difference of resonant wavelengths in first-order mode from Figure 3(a), the first-order mode of multipeak PIT has disorderly phenomenon because of the interference between two random resonance wavelengths.

4. Conclusion

In this paper, we have numerically investigated a multimode plasmonically induced transparency (PIT) device operating in the mid-infrared region. The proposed model consisted of one graphene ribbon as main waveguide and two graphene sheets as resonators. In phase coupling systems, the simulation results evidently demonstrate the multimode PIT phenomenon. On the other hand, the multipeak PIT spectral responses emerge in four graphene sheets system. For phase coupling, the off-to-on PIT feature could be effectively achieved by varying the Fermi level of one graphene sheet resonator. Besides, by turning off one graphene resonator, this structure can achieve filter effects. The filter wavelength can be adjusted via varying the Fermi energy of graphene resonator without optimizing geometric parameters. All of the observed PIT spectral features are of great importance to provide platforms for the design of the dynamical control PIT systems. The work may open up new ways for the development of compact thermal plasmonic devices such as tunable sensors, spectral splitter, and switchers in mid-infrared region.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the National Basic Research Program of China (Grant no. 2013CBA01702), the National Natural Science Foundation of China (Grant no. 61377016, 61308017, and 11504139), and the Natural Science Foundation of Jiangsu Province (Grant no. BK20140167).