Show Summary Details
The one-stop-shop for nano science.

# nano Online

More options …

## Search

Search publication

## Background

Surface plasmon polaritons (SPPs) are surface waves that propagate along the boundary surface between dielectric and metallic materials with fields decaying exponentially in both sides, thereby creating the subwavelength confinement of electromagnetic waves [ 1]. These are mainly electromagnetic modes resulting from the resonant interaction between light waves and the collective electron oscillations, which leads to its unique properties [ 2]. Plasmonic nanostructures offer the potential to overcome diffraction limits in dielectric structures, enabling us to miniaturize optical devices [ 3]. For example, plasmonic has been widely researched in integrated photonic circuits [ 4], photonic crystals [ 5], optical antennas [ 6, 7], nano-laser [ 8], data recording [ 9], filters [ 10], refractive index sensor [ 11], biological sensors [ 12], metalens [ 13], plasmonic lens [ 14], and so forth. Among the structures based on SPPs, the metal-insulator-metal (MIM) structure has been investigated extensively in designing plasmonic Bragg reflector. For example, periodic changes in the dielectric materials of the MIM waveguides have been proposed to design effective filtering around the Bragg frequency [ 15]; the thick-modulated and index-modulated Bragg reflectors have been reported to widen bandgap [ 16]; metal-embedded MIM structure also has been studied to improve the performance of conventional step profile MIM plasmonic Bragg reflectors (PBRs) [ 17]. However, plasmonic materials, usually noble metals, are hardly tunable and have great ohmic losses at the wavelength regimes of interest, therefore limiting their potential for some specific applications.

Graphene, a single layer of carbon atoms densely arranged into a honeycomb pattern, has been widely explored as a newly alternative to plasmonic material [ 18, 19]. Graphene plasmonics, similar to metal plasmonics at the visible region, can be easily induced in the near-infrared to terahertz (THz) regime. In particular, the surface charge density, namely chemical potential, can be actively modified by chemical doping or external gate voltage, thus giving rise to dramatic changes in the optical properties [ 20]. Additionally, SPPs bound to graphene display a strong field confinement, already verified by experiments [ 21, 22]. These remarkable and outstanding properties in turn enable a utility optical material in optoelectronic applications. In recent years, great attention has been focused on graphene-based plasmonic waveguides [ 2327]. de Abajo et al. have researched the propagation properties of graphene plasmonic waveguide constituted by individual and paired nanoribbons [ 28]. The tunable nano-modulators based on graphene plasmonic waveguide modulators have been proposed and numerically demonstrated [ 29]. Lu et al. have designed a slow-light waveguide based on graphene and silicon-graded grating [ 30]. Wang et al. have utilized a graphene waveguide achieving a tunable plasmonic Bragg reflector [ 31].

In this paper, we propose a PBR structure consisting of a single-layer graphene and silicon grating and numerically study its performance. We employ a silica spacer layer to separate the monolayer graphene and silicon grating and an external voltage gating to tune the surface conductivity of graphene. The finite element method (FEM) [ 32] has been utilized to perform the simulation work. We demonstrate that SPP modes on graphene exhibit a stopband around the Bragg wavelengths. Based on Bragg scattering condition, several modulation schemes have been used to adjust the characteristics of transmission spectrum. Furthermore, we introduce a defect into PBR, consequently realizing a resonant defect mode with a high and tunable Q factor. Based on the discussion of one defect cavity, we further study the multi-defect cases. At last, the PBR unit is designed into a multi-step pattern to reduce the rippling sidelobes and insertion loss. Such proposed designs, we believe, may help build some actively tunable modulators.

## Methods

As plotted in Fig.  1a, the proposed PBR in this work is composed of a single-layer graphene and a silicon grating substrate between which a silica layer has been embedded. And major structural parameters are labeled in Fig.  1a, b. In our investigation, the incident light at mid-infrared regime intrigues an excitation of a transverse magnetic (TM)-polarized SPP mode propagating along the graphene sheet. By solving the Maxwells equations with boundary conditions [ 33], we obtain the dispersion relation for TM modes supported on the graphene layer which is surrounded with air and silica: $\frac{{\epsilon }_{\mathrm{Air}}}{\sqrt{{{k}_{\mathrm{GSPP}}}^{2}-\frac{{\epsilon }_{\mathrm{Air}}{\omega }^{2}}{{c}^{2}}}}+\frac{{\epsilon }_{{\mathrm{S}\mathrm{i}\mathrm{O}}_{2}}}{\sqrt{{{k}_{\mathrm{GSPP}}}^{2}-\frac{{\epsilon }_{{\mathrm{S}\mathrm{i}\mathrm{O}}_{2}}{\omega }^{2}}{{c}^{2}}}}=\phantom{\rule{0.25em}{0ex}}-\frac{\sigma \left(\omega ,\phantom{\rule{0.25em}{0ex}}{k}_{\mathrm{GSPP}}\right)i}{\omega {\epsilon }_{0}}$1

Fig. 1

a, b 3D and 2D schematic illustration of PBR. A single layer of graphene and a silicon (Si) grating substrate between which a silica ( SiO 2) layer has been embedded. The Si substrate has a groove grating structure with period number N, period w 1+ w 2, groove width w 2, and groove depth d. The constant distance ( D) between graphene and Si grating is 100 nm. c, d The real part of effective refractive index ( Re( N eff)) in relation to the gate voltage ( V g) and the thickness of silica layer ( h). e, f The propagation distance ( L GSPP) in relation to the gate voltage ( V g) and the thickness of silica layer ( h)

Since we only consider the non-retarded regime ( k GSPP >> ω), the Eq. ( 1) can be simplified to ${k}_{\mathrm{GSPP}}\approx \phantom{\rule{0.25em}{0ex}}{\epsilon }_{0}\frac{{\epsilon }_{\mathrm{Air}}+{\epsilon }_{{\mathrm{S}\mathrm{i}\mathrm{O}}_{2}}}{2}\frac{2i\omega }{\sigma \left(\omega ,{k}_{\mathrm{GSPP}}\right)}$2

Here, k GSPP is the wave vector of SPPs on graphene layer, and the dielectric constants of air and SiO 2 are assumed to be 1 and 3.9, respectively. The optical conductivity of graphene is σ( ω, k GSPP) determined by Kobo formula [ 34]. At the mid-infrared frequency range, σ( ω, k GSPP) can be simplified into a Drude-like equation [ 35]: $\sigma \left(\omega \right)=\phantom{\rule{0.25em}{0ex}}\frac{{e}^{2}{\mu }_{c}}{\pi {\hslash }^{2}}\frac{i}{\omega +i{\tau }^{-1}}$3

In graphene layer, τ denoting the relaxation time can be expressed as μμ c/( ev f 2), which relates to the carrier mobility μ and Fermi velocity v f = 10 6 m/s [ 36]. The carrier mobility is reasonably chosen to be μ = 20,000 cm 2 V −1 s −1 from experiment results [ 37]. And, μ c =  ħv f( πn) 1/2 is the chemical potential where the surface charge carrier density is expressed as n =  ε 0 ε d V g/( eh) [ 18]. Here, ε 0 and ε d are the dielectric constants of free space and SiO 2, respectively. V g is the applied gate voltage, e is the electron charge, and h is the thickness of silica layer. In Fig.  1a, b, h equals to D at non-groove sections and h equals to D +  d at groove sections. This expression also indicates that the chemical potential can be induced by not only a voltage gate but also the thickness of silica layer. From the above equations, a more specific graphene plasmonic dispersion relation is obtained as follows: ${k}_{\mathrm{GSPP}}=\phantom{\rule{0.25em}{0ex}}\frac{\pi \hslash {\epsilon }_{0}\left({\epsilon }_{\mathrm{Air}}+{\epsilon }_{{\mathrm{S}\mathrm{i}\mathrm{O}}_{2}}\right)}{{e}^{2}{v}_{f}{\left(\pi \frac{{\epsilon }_{0}{\epsilon }_{Si{O}_{2}}{V}_{g}}{eh}\right)}^{\frac{1}{2}}}\left(1+\frac{i}{\tau \omega }\right){\omega }^{2}$4

Another important parameter derived from the above equation is N GSPP =  k GSPP/ k 0—the effective refractive index of GSPP, which shows the ability to confine GSPP on graphene. The propagation length is defined as L GSPP = 1/[2 k 0Im( N GSPP)] featuring the GSPP propagation loss. Throughout the paper, the influence of substrate silicon on the dispersion relation is negligible when the silica layer is above 100 nm [ 30]. The dependence of Re( N GSPP) and L GSPP on the gate voltage V g and the thickness of silica layer h are illustrated in Fig.  1c–f where the wavelength range is 6 to 9 μm. Obviously, from Fig.  1c, d, the Re( N GSPP) shows a significant increase when the gate voltage V g decreases or the thickness of silica layer h goes up. In Fig.  1e, f, however, a longer propagation distance is achieved by a growing V g or a decreasing h. Hence, the two important factors should be both taken into consideration for our designs.

Generally, for wavelength-sensitive operations, a plasmonic Bragg reflector is constructed by periodically modulating the effective refractive index of the waveguide. There are some popular accesses to achieve this such as width modulation [ 38] and refractive index modulation [ 39]. In our work, a graphene-based Bragg reflector is formed by periodically modulating N GSPP. According to the aforesaid discussion, this can be realized by alternatively varying the thickness of silica layer yielding a silicon grating substrate, shown in Fig.  1a, b. Thus, the Bragg scattering condition [ 39] in our case can be formulated as: ${w}_{1}\mathrm{R}\mathrm{e}\left({N}_{\mathrm{GSPP}1}\right)+{w}_{2}\mathrm{R}\mathrm{e}\left({N}_{\mathrm{GSPP}2}\right)=m\frac{{\lambda }_{b}}{2}$5

Here, λ b is the Bragg wavelength and m is an integer assumed to be 1 in our discussion. N GSPP1 and N GSPP2 are the effective refractive index of GSPP on differently doped areas of graphene, respectively (see Fig.  1a, b). The Bragg wavelength will be stopped when Eq. ( 5) is satisfied.

## Results and discussion

At first, we discuss the influence of period number on transmission spectra of PBR. The parameters are set as w 1 =  w 2 = 40 nm, D =  d = 100 nm, and V g = 45 V. The period number N is ranged from 4 to 14. The simulated transmission spectra for different period numbers are shown in Fig.  2a. As the period number decreases, the propagation loss is lower whereas the stopband is narrower. To balance these two trends, we choose period number to be 8 in the following discussion. Besides, these spectra all display some sidelobes outside the stopband caused by light scattering at the end of PBR. Next, we study the effects of grating groove depth d and gate voltage V g on the operating wavelengths of PBR. Figure  2b indicates a pronounced red-shift of Bragg wavelengths and a widened stopband with a growing groove depth d; nevertheless, the transmission is gradually lowered as the groove depth d increases. In Fig.  2c, we see a blue-shift of Bragg wavelengths and a growing transmission when the gate voltage goes up. The shifting effect can be attributed to the alteration of N GSPP with an increasing groove depth d or a growing gate voltage V g. And the modification of L GSPP by varying groove depth d or gate voltage V g accounts for the change of transmission above. We further study another two parameters from the Bragg condition, w 1 and w 2. As displayed in Fig.  2d, we see a noticeable red-shift of the central wavelengths while the width of stopband is almost the same when w 1 is increasing. By varying w 2, the transmission spectra show almost similar features.

Fig. 2

Simulated transmission spectra. a The PBR with the different period numbers for the groove depth d = 100 nm, the grating period w 1 +  w 2 = 80 nm ( w 1 =  w 2), and the gate voltage V g = 45 V. b The PBR with the different groove depth for the period number N = 8, the grating period w 1 +  w 2 = 80 nm ( w 1 =  w 2), and the gate voltage V g = 45 V. c The PBR with the different gate voltages for the period number N = 8, the groove depth d = 100 nm, and the grating period w 1 +  w 2 = 80 nm ( w 1 =  w 2). d The PBR with the different w 1 for the period number N = 8, the groove depth d = 100 nm, and the gate voltage V g = 45 V

In addition, we introduce a nano-cavity into the plasmonic Bragg grating causing a defect in the PBR’s periodicity. At some specific wavelengths, resonant modes will be formed around the nano-cavity, and accordingly, there will be a transmission peak inside the stopband [ 40]. In Fig.  3a, we set the nano-cavity in the center of PBR and as same as that in the study above. The thickness of silica layer below the nano-cavity graphene is set as D d = 100 nm, and firstly, we discuss the length of defect cavity w d ranging from 70 to 120 nm. In Fig.  3b, we see clear peaks at the stopping range of transmission spectra as we have expected. By lengthening the nano-cavity step by step, we find an obvious red-shifting of the resonant defect modes with the stopband unmoving. Secondly, we fix the length of nano-cavity as 90 nm and research the influence of D d on the defect modes. The simulated results plotted in Fig.  3c present similar characteristics to that in Fig.  3b. Therefore, we have an access to tune the defect modes by varying the length of nano-cavity or the thickness of silica layer below the nano-cavity graphene. Next, we employ varied gating voltage on graphene sheet, and consequently, the spectra response is shifting as a whole, exhibited in Fig.  3d. Lastly, we set the defect cavity at different locations and study the features of transmission spectra. In Fig.  3e, defect cavity is moved from one side to the center. To get into more principles behind the off-to-on effect, we plot the electric field profiles at λ = 7.105 μm for the three cases (defect at third, fifth, and seventh location sites).

Fig. 3

a 3D schematic illustration of PBR by introducing a defect cavity located at ninth site with the width w d and the distance D d between graphene and Si. b Transmission spectra of the defecting PBR with different w d, D d = 100 nm, and the gate voltage V g = 45 V. c Transmission spectra of the defecting PBR with different D d, w d = 90 nm, and the gate voltage V g = 45 V. d Transmission spectra of the defecting PBR with different V g, w d = 90 nm, and D d = 100 nm. e Transmission spectra of the defecting PBR with different defect locations, w d = 90 nm, D d = 100 nm, and V g = 45 V

As shown in Fig.  4, pronounced resonant modes form at defect regions dividing the original Bragg reflector into two new Bragg reflectors. When the reflection by the two new Bragg reflectors adds up destructively, a high transmission occurs. However, when the two new Bragg reflectors are different from each other, the reflected beams from them usually cannot perform a well-destructive interference, which explain the off-to-on effect.

Fig. 4

The electric field profiles ( E y) for the defecting PBR with different defect locations, λ = 7.105 μm, w d = 90 nm, D d = 100 nm, and V g = 45 V. The dashed lines denote the defect-mode areas

Since the defect cavity works differently when they are located at different sites, we further propose multiple defect cavities. At first, we compare the transmission spectra for three cases: one-defect, two-defect, and three-defect in PBR. All defect cavities are set at the same condition: w d = 90 nm and D d = 100 nm. Through Fig.  5a, it is demonstrated that multiple transmission peaks inside the stopband can be induced by adding the number of defect cavity in PBR. In Fig.  5b, we plot the electric field profiles at λ = 7.105 μm for these three cases. The red dashed circles denoting the resonant areas indicate that these defect modes from each nano-cavity will interfere with each other leading to these amazing features in transmission spectra. Besides, stopband can be greatly widened with multiple defect cavities. Then, we place one defect cavity in the center and place another one at different locations. Figure  6a points out that there is a notable splitting of defect-resonance mode in the two-defect case. It should be noted, however, that this splitting effect is only achieved when these two defect cavities are close enough to each other. As illustrated in Fig.  6b, the closer these two defects are together, the stronger the destructive interference between them, which gives rise to the splitting effect. Furthermore, we apply various gate voltages to the two-defect case producing a shifting effect of the multiple transmission peaks, seen in Fig.  7a. Similar shifting effect is also found in the three-defect case by varying gate voltage, shown in Fig.  7b.

Fig. 5

a Transmission spectra of the defecting PBR at three different cases: one-defect, two-defect, and three-defect. b The electric field profiles ( E y) at λ = 7.105 μm for these three cases. The red dashed circles denote the resonant areas

Fig. 6

a Transmission spectra of the PBR with double defect for three different location combinations: third–ninth, fifth–ninth, and seventh–ninth. b The electric field profiles ( E y) for these three cases at λ = 7.105 μm (pointed by the arrows in transmission spectra). The red dashed circles denote the central defect-mode region

Fig. 7

a Transmission spectra of the PBR with two defects for different gate voltages. b Transmission spectra of the PBR with three defects for different gate voltages

In the end, we design the PBR unit into a multi-step pattern to deal with the high insertion loss and severe rippling in transmission spectra resulting from the abrupt change of N GSPP in the groove depth [ 38]. A three-step and a six-step version of PBR unit are well illustrated in Fig.  8a. After comparing the transmission spectra in Fig.  8b, we find the expected enhancement on transmission spectra and rippling suppression. Additionally, the stopband is gradually narrowed when the PBR unit is changed into more steps. The multiple steps in a PBR unit actually adding multiple reflections into the Bragg reflection process make it harder to satisfy the Bragg condition, which results in the sidelobe suppression and narrowed stopband.

Fig. 8

a 2D schematic illustration of PBR unit designed into a multi-step pattern: a three-step PBR unit and a six-step PBR unit, respectively. b Transmission spectra of PBR for three different cases. The red line represents the original PBR design (groove depth d = 120 nm); the green one and the blue one denote the three-step and six-step pattern cases, respectively

## Conclusions

In conclusion, we design a PBR structure consisting of a single-layer graphene and silicon grating and numerically study its performance. We employ an external voltage gating to tune the surface conductivity of graphene. It is found that SPP modes on graphene exhibit a stopband around the Bragg wavelengths. Based on Bragg scattering condition, we discuss several modulation schemes to adjust the characteristics of transmission spectrum. Furthermore, we introduce a nano-cavity into PBR, consequently realizing a resonant defect mode. We further propose multi-defect PBR and achieve multiple peaks inside the stopband. At last, by designing the PBR unit into multi-step pattern, we lower the insertion loss and suppress the rippling in transmission spectra. We hope all the proposed designs above can help pave new ways in actively tunable modulation application.

## Acknowledgements

This work is supported by the National Natural Science Foundation of China (Grant Nos. 11504139, 11504140), the Natural Science Foundation of Jiangsu Province (Grant Nos. BK20140167, BK20140128), Jiangsu Provincial Research Center of Light Industrial Optoelectronic Engineering and Technology (Grant No. BM2014402), and the Key Laboratory Open Fund of Institute of Semiconductors of CAS (Grant No. KLSMS-1405).

JW and CS conceived the idea. CS, XX, and YL calculated the properties of the proposed structure, analyzed the data, and wrote the manuscript. JW and ZDH revised the manuscript. ZDH and JW supervised the project. All authors read and approved the final manuscript.

The authors declare that they have no competing interests.

## References

• 1.

Barnes WL, Dereux A, Ebbesen TW (2003) Surface plasmon subwavelength optics. Nature 424:824–830Google Scholar

• 2.

Ozbay E (2006) Plasmonics: merging photonics and electronics at nanoscale dimensions. Science 311:57–58Google Scholar

• 3.

Zia R, Schuller JA, Chandran A (2006) Plasmonics: the next chip-scale technology. Mater Today 9:20–27Google Scholar

• 4.

Bozhevolnyi SI, Volkov VS, Devaux E, Laluet JY, Ebbesen TW (2006) Channel plasmon subwavelength waveguide components including interferometers and ring resonators. Nature 440:508–511Google Scholar

• 5.

Zhang AJ, Guo ZY, Tao YF, Wang W, Mao XQ, Fan GH, Zhou KY, Qu SL (2015) Advanced light-trapping effect of thin-film solar cell with dual photonic crystals. Nanoscale Res Lett 10:1–10Google Scholar

• 6.

Muhlschlegel P, Eisler HJ, Martin OJF, Hecht B, Pohl DW (2006) Resonant optical antennas. Science 308:1607–1609Google Scholar

• 7.

Wang J, Wang Y, Zhang X, Yang K, Song Y, Liu S (2010) Splitting and unidirectional excitation of surface plasmon polaritons by two uniform metallic nanoslits with a nanocavity antenna. J Mod Opt 57:1630–1634Google Scholar

• 8.

Oulton RF, Sorger VJ, Zentgraf T, Ma RM, Gladden C, Dai L, Bartal G, Zhang X (2009) Plasmon lasers at deep subwavelength scale. Nature 461:629–632Google Scholar

• 9.

Zijlstra P, Chon JW, Gu M (2009) Five-dimensional optical recording mediated by surface plasmons in gold nanorods. Nature 459:410–413Google Scholar

• 10.

Tang B, Wang J, Xia X, Liang X, Song C, Qu S (2015) Plasmonic induced transparency and unidirectional control based on the waveguide structure with quadrant ring resonators. Appl Phys Express 8:032202Google Scholar

• 11.

Luo LB, Ge CW, Tao YF, Zhu L, Zheng K, Wang W, Sun YX, Shen F, Guo ZY (2016). High-efficiency refractive index sensor based on the metallic nanoslit arrays with gain-assisted materials. Nanophotonics 5:139-146Google Scholar

• 12.

Skorobogatiy M, Kabashin AV (2006) Photon crystal waveguide-based surface plasmon resonance biosensor. Appl Phys Lett 89:143518Google Scholar

• 13.

Wang W, Guo ZY, Zhou KY, Sun YX, Shen F, Li Y, Qu SL, Liu ST (2015) Polarization-independent longitudinal multi-focusing metalens. Opt Express 23:29855–29866Google Scholar

• 14.

Chen SW, Huang YH, Chao BK, Hsueh CH, Li JH (2014) Electric field enhancement and far-field radiation pattern of the nanoantenna with concentric rings. Nanoscale Res Lett 9:681Google Scholar

• 15.

Qu S, Song C, Xia X, Liang X, Tang B, Hu ZD, Wang J (2016) Detuned plasmonic Bragg grating sensor based on defect metal-insulator-metal waveguide. Sensors 16:784Google Scholar

• 16.

Liu JQ, Wang LL, He MD, Huang WQ, Wang D, Zou BS, Wen S (2008) A wide bandgap plasmonic Bragg reflector. Opt Express 16:4888–4894Google Scholar

• 17.

Tian M, Lu P, Chen L, Liu DM, Peyghambarian N (2012) Plasmonic Bragg reflectors based on metal-embedded MIM structure. Opt Commun 285:5122–5127Google Scholar

• 18.

Novoselov KS, Geim AK, Morozov SV, Jiang D, Zhang Y, Dubonos SV, Grigorieva IV, Firsov AA (2004) Electric field effect in atomically thin carbon films. Science 306:666–669Google Scholar

• 19.

Geim AK, Novoselov KS (2007) The rise of graphene. Nature Mater 6:183–191Google Scholar

• 20.

Wang F, Zhang Y, Tian C, Girit C, Zettl A, Crommie M, Shen YR (2008) Gate-variable optical transitions in graphene. Science 320:206–209Google Scholar

• 21.

Fei Z, Rodin AS, Andreev GO, Bao W, Mcleod AS, Wagner M, Zhang LM, Zhao Z, Thiemens M, Dominguez G, Fogler MM, Castro Neto AH, Lau CN, Keilmann F, Basov DN (2012) Gate-tuning of graphene plasmons revealed by infrared nano-imaging. Nature 487:82–85Google Scholar

• 22.

Brar VW, Jang MS, Sherrott M, Lopez JJ, Atwater HA (2013) Highly confined tunable mid-infrared plasmonics in graphene nanoresonators. Nano Lett 13:2541–2547Google Scholar

• 23.

Zheng J, Yu L, He S, Dai D (2015) Tunable pattern-free graphene nanoplasmonic waveguides on trenched silicon substrate. Science 5:7987–7987Google Scholar

• 24.

Kong XT, Bai B, Dai Q (2015) Graphene plasmon propagation on corrugated silicon substrates. Opt Lett 40:1–4Google Scholar

• 25.

Zhu X, Yan W, Mortensen NA, Xiao S (2013) Bends and splitters in graphene nanoribbon waveguides. Opt Express 21:3486–3491Google Scholar

• 26.

Yuan H, Yang H, Liu P, Jiang X, Sun X (2014) Mode manipulation and near-THz absorptions in binary grating-graphene layer structures. Nanoscale Res Lett 9:90Google Scholar

• 27.

Wang X, Xia X, Wang J, Zhang F, Hu ZD, Liu C (2015) Tunable plasmonically-induced transparency with unsymmetrical graphene-rings resonators. J Appl Phys 118:013101Google Scholar

• 28.

Christensen J, Manjavacas A, Thongrattanasiri S, Koppens FHL, Abajo FJGD (2012) Graphene plasmon waveguiding and hybridization in individual and paired nanoribbons. ACS Nano 6:431–440Google Scholar

• 29.

Lao J, Tao J, Wang Q, Huang X (2014) Tunable graphene-based plasmonic waveguides: nano modulators and nano attenuators. Laser Photon Rev 8:569–574Google Scholar

• 30.

Lu H, Zeng C, Zhang Q, Liu X, Hossain MM, Reineck P, Gu M (2015) Graphene-based active slow surface plasmon polaritons. Sci Rep 5:8443Google Scholar

• 31.

Tao J, Yu X, Hu B, Dubrovkin A, Wang Q (2014) Graphene-based tunable plasmonic Bragg reflector with a broad bandwidth. Opt Lett 39:271–274Google Scholar

• 32.

Jin J (2002) The finite element method in electromagnetics. New York: Wiley-IEEE PressGoogle Scholar

• 33.

Jablan M, Buljan H, Soljacic M (2009) Plasmonics in graphene at infra-red frequencies. Phys Rev B 80:308–310Google Scholar

• 34.

Fal'kovskii LA (2008) Optical properties of graphene. J Exp Theor Phys 115:496–508Google Scholar

• 35.

Hanson GW (2008) Quasi-transverse electromagnetic modes supported by a graphene parallel-plate waveguide. J Appl Phys 104:084314Google Scholar

• 36.

Gao W, Shu J, Qiu C, Xu Q (2012) Excitation of plasmonic waves in graphene by guided-mode resonances. ACS Nano 6:7806–7813Google Scholar

• 37.

Thongrattanasiri S, Koppens FHL, Abajo FJGD (2012) Complete optical absorption in periodically patterned graphene. Phys Rev Lett 108:799–802Google Scholar

• 38.

Liu Y, Liu Y, Kim J (2010) Characteristics of plasmonic Bragg reflectors with insulator width modulated in sawtooth profiles. Opt Express 18:11589–11598Google Scholar

• 39.

Hossieni A, Massoud Y (2006) A low-loss metal-insulator-metal plasmonic Bragg reflector. Opt Express 14:11318–11323Google Scholar

• 40.

Chang YJ, Chen CY (2013) Ultracompact, narrowband three-dimensional plasmonic waveguide Bragg grating in metal/multi-insulator/metal configuration. Appl Opt 52:889–896Google Scholar