We investigate the performance of bends and splitters in graphene nanoribbon waveguides. Although the graphene waveguides are lossy themselves, we show that bends and splitters do not induce any additional loss provided that the nanoribbon width is sub-wavelength. We use transmission line theory to qualitatively interpret the behavior observed in our simulation. Our results pave a promising way to realize ultra-compact devices operating in the terahertz region.
© 2013 OSA
Graphene has attracted considerable attention due to its unique electronic and optical properties [1–5]. The high carrier mobility approaching 200,000 cm2V−1s−1 offers the possibility of building high-speed optoelectronic devices, such as photodetectors  and optical modulators . Despite being a single atomic layer, graphene can couple strongly to light, enabling a constant absorbance of πa = 2.3% in the infrared-to-visible spectral range , where a = e2/ħc is the fine-structure constant. More importantly, the optical properties of graphene can be substantially modified by changing the Fermi energy EF[9, 10]. For a slightly doped graphene (EF < ħω/2), the conductivity of graphene is dominated by the interband contributions. In the terahertz range when EF > ħω/2, the intraband contribution becomes relevant, thus allowing graphene to support graphene-plasmon polaritons (GPPs), i.e., the collective oscillation of carriers. Compared to conventional plasmonic materials, e.g. gold, doped-graphene plasmon has promising properties including tunability, extreme field confinement and low propagation loss.
Graphene-plasmon polaritons supported by doped graphene have been widely investigated [11–15]. It has been shown that remarkable optical absorption enhancement [12, 13, 16] in the graphene can be achieved by patterning the single sheet of doped graphene. Moreover, the GPP can also be confined laterally, e.g., in graphene ribbons. Recently, much attention has been focused on investigating propagation properties of GPP in graphene nanoribbons [17, 18], which can be used as a promising candidate to construct future high-speed compact-devices working in the terahertz region. Obviously it is indispensable to have graphene nanoribbon bends and splitters as the basic structures in these compact devices. However, there is a lack of characterization of the GPP in these geometries. In this paper, we investigate the performance of bends and power splitters in graphene nanoribbon waveguides, particularly focusing on the question whether these bends and splitters will induce reflection or excess scattering loss on top of the propagation loss in the nanoribbon waveguides themselves. Our analysis in this paper relies on a classical electromagnetic description where graphene is described by a local conductivity, and the results presented here do not include the microscopic details of the edges of graphene nanoribbons .
2. Mode properties of GPP in a graphene nanoribbon
Let us first investigate mode properties of GPP in a freestanding graphene nanoribbon, as shown in the inset of Fig. 1(a), with an ultra-thin width d<50 nm. The inset in Fig. 1(a) illustrates a cross-section of the nanoribbon and its corresponding axis. Due to the translational symmetry in the z direction, the electric solution of the nanoribbon waveguide has the form E⃗(r⃗, ω, t) = E⃗(x,y)exp(ikz(ω)z−iωt), where kz(ω) is the wave vector in the propagation direction z. For a waveguide in which the material is lossy, kz(ω) is complex. An effective refractive index neff of the waveguide mode is defined by neff = Re(kz(ω))/k0, and the propagation length L of the mode is given by L = 1/[2Im(kz(ω))], where k0 = 2π/λ0 (wavelength in vacuum λ0). In the simulations, the graphene ribbon is modeled as a thin layer of thickness t = 0.5 nm and characterized by a dielectric function 1 + iσ/(ωtε0), where σ is the surface conductivity of graphene, ε0 is the vacuum permittivity, and ω is the light frequency. Note that the simulation results do not rely on the thickness of the film as long as the meshing is sufficiently fine. In this paper we consider the case where the photon energy is always less than 2EF, and under this circumstance the intraband contribution to the conductivity dominates. The intraband conductivity σintra of the graphene film can be taken from the standard random-phase approximation [19, 20] In this work, we consider Fermi energy of 0.3 eV for all calculations, and the mode calculations have been performed by using a commercial available finite-element-method software package (COMSOL).
Figure 1(a) shows the effective refractive indices neff of the GPP supported by freestanding nanoribbons for two different widths (d = 20, 30 nm). In the frequency range of our interest, the nanoribbons with the widths of 20 nm and 30 nm only support the fundamental waveguide mode. It has been shown that a wide graphene ribbon may support many modes including edge modes and ”bulk” modes [17, 18]. With the decrease of the ribbon width, there exists cut-offs for the high-order modes. Eventually, the nanoribbon only has the fundamental mode with even symmetry (based on the value of Ex at y = 0), similar to the dispersion of the surface-plasmon in a metal strip . The effective refractive indices shown in Fig. 1(a) of the fundamental GPP of the nanoribbons have large values, e.g. 48 for the case of λ0 = 15 μm, and d = 30 nm, indicating that the GPP mode is extremely confined near the nanoribbon. This is also confirmed by the mode profile shown in the inset of Fig. 1(b). The effective refractive index becomes larger when decreasing the nanoribbon width from 30 nm to 20 nm. Compared to the nobel metals, graphene generally has a relatively large conductivity, thus leading to a lower propagation loss. Propagation lengths L as a function of the wavelength are plotted in Fig. 1(b). As an example, when λ0 = 10 μm, and d = 30 nm, L reaches 0.12λ0, equivalent to 1.2 μm. The propagation loss of the GPP becomes weaker when the wavelength decreases, which is mainly attributed to the optical property of the graphene material. On the other hand, the propagation length becomes shorter by shrinking the nanoribbon width for a specific working wavelength. It can be well understood by the fact that when the effective refractive index of the GPP mode becomes larger, there will be a larger field overlap with the graphene, thus leading to a higher propagation loss. Similar to surface plasmons supported by a metallic waveguide, there is also a competition in the graphene ribbon waveguide between the mode confinement and propagation length [22–25].
3. Bends in graphene nanoribbon waveguides
We now turn to the question whether waveguide bends (breaking translational invariance) lead to reflection or excess absorption? As mentioned above, the graphene nanoribbon waveguide is generally lossy. In order to evaluate the performance of the bends, we calculate the transmission coefficient and normalize it to that of a straight graphene ribbon waveguide with the same propagation length. Propagation properties of the graphene nanoribbon waveguide are analyzed by a finite-integral-method software package (CST microwave studio). In all cases, the width d of the graphene nanoribbons is chosen to be much smaller than the wavelength, so that only the fundamental waveguide mode exists, see the discussion above.
Figure 2 presents the normalized transmission through a 90° sharp graphene nanoribbon bend, therefore the additional loss associated with reflection and excess scattering loss for the sharp graphene waveguide can be obtained. One can see in Fig. 2 that there is no additional loss at longer wavelengths. Note that a similar high transmission through sharp bends can also be realized in the photonic crystal systems [26–28]. For the case of d = 30 nm, the limiting wavelength at which the transmission coefficient decreases below 99% is 14 μm. The limiting wavelength is blue-shifted to 12 μm when decreasing the nanoribbon width from 30 nm to 20 nm. To evaluate the transmittance through the sharp graphene ribbon, we can also employ the analogy between the subwavelength waveguiding in photonics and electronics and use a common quasistatic approximation [29, 30]. With this analogy and quasistatic approximation (valid as long as the effective width of a waveguide is much smaller than the effective wavelength governed by λ/neff), we can consider the sharp bend equivalent to a junction between two transmission lines with the same characteristic impedance. Within the quasistatic approximation no additional loss emerges for the bend. It is naturally understood that the operating wavelength range for having additional bend loss widens with the decrease of d. For a shorter wavelength, e.g. λ = 6 μm, the transmission is relatively low. The effective refractive index of the GPP at the shorter wavelength shown in Fig. 1(a) increases dramatically, and the effective wavelength of the GPP becomes comparable with the size of the bend, therefore leading to strong reflection. The inset in Fig. 2 shows the top-view magnetic field amplitude distribution at y = 0 when light passes through the sharp corner for the case of d = 30 nm at λ0 = 15 μm. To illustrate this behavior clearly, we neglect the material loss of the graphene. One can see that light can be perfectly turned by 90° without any noticeable back-scattering or radiation.
4. Splitters in graphene nanoribbon waveguides
Next, we focus on the propagation property of graphene nanoribbon splitters, which is another key element. The calculated normalized transmission coefficients as a function of wavelength for graphene T-shape nanoribbon waveguides are plotted in Fig. 3 for the cases of d = 20, 30 nm. The inset shows a top-view of the structure. Similar to the result for the nanoribbon bend, the transmission at the output graphene waveguides becomes larger as the wavelength increases, and it saturates at the longer wavelengths, reaching 42%. We find that the limiting wavelengths are almost the same as those for the graphene bends shown in Fig. 2. Note that the widths for the input and output graphene waveguides have the same values as d. Due to the symmetry of the structure, the transmitted power is equally distributed between the two output waveguide branches, so that the total transmission is close to 84%. With the quasistatic approximation, we can replace the splitter configuration by an equivalent network consisting of a junction of three transmission lines with the same characteristic impedance Z0. The combination of the two output transmission lines are connected to the input transmission line, thus resulting in the equivalent load impedance with a value of ZL = 2Z0. The transmission coefficient at each output governed by 1 − [(ZL − Z0)/(ZL + Z0)]2 is close to 88%. The result obtained by the transmission lines has a small deviation when compared with our simulated result. When using the transmission line theory, one has to assume that the energy of the guided mode is mainly confined in the waveguide and that the distribution of the transverse electromagnetic field is nearly uniform within the waveguide. This is only qualitatively correct for the GPP supported by the graphene nanoribbon waveguide. However, we can still use the transmission line theory to qualitatively interpret the behavior observed in our simulation. In order to improve the transmission coefficient of the graphene nanoribbon splitter, we can adjust the characteristic impedance of the input waveguide Zin by changing the width din of the input waveguide. Figure 4 illustrates the calculated reflection coefficient of the graphene T-shaped splitter at λ0 = 12 μm as a function of din/dout, where dout, the thickness of the output waveguide branches, is fixed at 20 nm. An optimal structure with nearly-zero reflection coefficient is observed at din/dout = 2.25. The lower-inset in Fig. 4 shows the top-view magnetic field amplitude distribution at y = 0 when light passes through the T-shape configuration for the case of din/dout = 2.25 at λ0 = 12 μm.
5. Summary and discussion
We have investigated the performance of bends and splitters in graphene nanoribbon waveguides, where the nanoribbon width is much smaller than the working wavelength. We have shown that sharp bends and splitters with no additional loss can be designed for longer wavelengths. The effective characteristic impendence model has been used to qualitatively interpret the simulation results. As discussed above, the GPP in the graphene nanoribbon has an extreme confinement, therefore making it difficult to couple light into the nanoribbon by use of the butt-coupling technology. Similar to the surface plasmon supported by the noble metals, a grating coupler could be used to excite the GPP in the nanoribbons. Recently, the propagating and localized graphene-plasmon polaritons have been visualized in real space by scattering-type scanning near-field optical microscopy [31, 32]. Our results further advance the potential of such devices by providing building blocks to construct graphene circuits for future ultra-compact terahertz devices.
This work is partly supported by the Center for Nanostructured Graphene, sponsored by the Danish National Research Foundation.
References and links
1. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, “Electric field effect in atomically thin carbon films,” Science 306, 666–669 (2004). [CrossRef] [PubMed]
2. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, “Two-dimensional gas of massless dirac fermions in graphene,” Nature 438, 197–200 (2005). [CrossRef] [PubMed]
4. A. N. Grigorenko, M. Polini, and K. S. Novoselov, “Graphene plasmonics,” Nat. Photonics 6, 749–758 (2012). [CrossRef]
9. Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, P. Kim, H. L. Stormer, and D. N. Basov, “Dirac charge dynamics in graphene by infrared spectroscopy,” Nat. Phys. 4, 532–535 (2008). [CrossRef]
10. F. Wang, Y. Zhang, C. Tian, C. Girit, A. Zettl, M. Crommie, and Y. R. Shen, “Gate-variable optical transitions in graphene,” Science 206, 206–209 (2008). [CrossRef]
11. M. Jablan, H. Buljan, and Soljacic, “Plasmonics in graphene at infrared frequencies,” Phys. Rev. B 80, 245435 (2009). [CrossRef]
14. L. Ju, B. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, H. A. Bechtel, X. Liang, A. Zettl, Y. R. Shen, and F. Wang, “Graphene plasmonics for tunable terahertz metamaterials,” Nat. Nanotechnol. 6, 630–634 (2011). [CrossRef] [PubMed]
15. F. Yang, J. R. Sambles, and G. W. Bradberry, “Tunable infrared plasmonic devices using graphene/insulator stacks,” Nature Nanotechnol. 7, 330–334 (2012). [CrossRef]
16. T. R. Zhan, F. Y. Zhao, X. H. Hu, X. H. Liu, and J. Zi, “Band structure of plasmons and optical absorption enhancement in graphene on subwavelength dielectric gratings at infrared frequencies,” Phys. Rev. B 86, 165416 (2012). [CrossRef]
17. A. Y. Nikitin, F. Guinea, F. Garcia-Vidal, and L. Martin-Moreno, “Edge and waveguide terahertz surface plasmon modes in graphene microribbons,” Phys. Rev. B 84, 161407 (R) (2011). [CrossRef]
18. J. Christensen, A. Manjavacas, S. Thongrattanasiri, F. Koppens, and F. J. Garcia de Abajo, “Graphene plasmon waveguiding and hybridization in individual and paired nanoribbons,” ACS Nano 6, 431–440 (2012). [CrossRef]
19. E. H. Hwand and S. Das Sarma, “Dielectric function, screening, and plasmons in two-dimensional graphene,” Phys. Rev. B 75, 205418 (2007). [CrossRef]
20. L. A. Falkovsky and A. A. Varlamov, “Space-time dispersion of graphene conductivity,” Eur. Phys. J. B 56, 281–284 (2007). [CrossRef]
21. S. A. Maier, Plasmonics: Fundamentals and Applications, 1st ed. (Springer, 2007).
22. S. A. Maier, P. G. Kik, H. A. Atwater, S. Meltzer, E. Harel, B. E. Koel, and A. A. G. Requicha, “Local detection of electromagnetic energy transport below the diffraction limit in metal nanoparticle plasmon waveguides,” Nat. Mater. 2, 229–232 (2003). [CrossRef] [PubMed]
23. S. I. Bozhevolnyi, V. S. Volkov, E. Devaux, J.-Y. Laluet, and T. W. Ebbesen, “Channel plasmon subwavelength waveguide components including interferometers and ring resonators,” Nature 440, 508–511 (2006). [CrossRef] [PubMed]
25. S. Xiao, J. Zhang, L. Peng, C. Jeppesen, R. Malureanu, A. Kristensen, and N. A. Mortensen, “Nealy-zero transmission through periodically modulated ultrathin metal films,” Appl. Phys. Lett. 97, 071116 (2010). [CrossRef]
26. A. Mekis, J. C. Chen, I. Kurland, S. Fan, P. R. Villeneuve, and J. D. Joannopoulos, “High transmission through sharp bends in photonic crystal waveguides,” Phys. Rev. Lett. 77, 3787–3790 (1996). [CrossRef] [PubMed]
27. J. D. Joannopoulos, R. D. Meade, and J. Winn, Photonic Crystals: Modling the Flow of Light, 1st. ed. (Princeton Univ., 1995).
28. S. S. Xiao and M. Qiu, “Study of transmission properties for waveguide bends by use of a circular photonic crystal,” Phys. Lett. A 340, 474–479 (2005). [CrossRef]
29. G. Veronis and S. H. Fan, “Bends and splitters in metal-dielectric-metal subwavelength plasmonic waveguides,” Appl. Phys. Lett. 87, 131102 (2005). [CrossRef]
31. J. Chen, M. Badioli, P. Alonso-Gonzalez, S. Thongrattanasiri, F. Huth, J. Osmond, M. Spasenovic, A. Centeno, A. Pesquera, P. Godignon, A. Z. Elorza, N. Camara, F. J. Garcia de Abajo, R. Hillenbrand, and F. H. L. Koppens, “Optical nano-imaging of gate-tunable graphene plasmons,” Nature 487, 77–81 (2012). [PubMed]
32. Z. Fei, A. S. Rodin, G. O. Andreev, W. Bao, A. S. McLeod, M. Wagner, L. M. Zhang, Z. Zhao, G. Thiemens, M. Dominguez, M. M. Fogler, A. H. Castro Neto, C. N. Lau, F. K., and D. N. Basov, “Gate-tuning of graphene plasmons revealed by infrared nano-imaging,” Nature 487, 82–85 (2012). [PubMed]