Demonstration of complex refractive index of graphene waveguide by microfiber-based Mach– Zehnder interferometer Baicheng Yao,1 Yu Wu,1,3 Zegao Wang,2 Yang Cheng,1 Yunjiang Rao,1,* Yuan Gong,1 Yuanfu Chen,2 and Yanrong Li2 1
2
Key Laboratory of Optical Fiber Sensing and Communications (Education Ministry of China), University of Electronic Science and Technology of China, Chengdu 610054, China State Key Laboratory of Electronic Thin Films and Integrated Devices, University of Electronic Science and Technology of China, Chengdu 610054, China 3 Department of Electronic Engineering, City University of Hong Kong, Kowloon, Hong Kong *
[email protected] Abstract: The complex refractive index (CRI) of graphene waveguide (GW) is of great importance for modeling and developing graphene-based photonic or optoelectronic devices. In this paper, the CRI of the GW is investigated theoretically and experimentally, it is found that the CRI of the GW will modulate the intensity and phase of transmitting light. The phase alterations are obtained spectrally by a Microfiber-based Mach–Zehnder interferometer (MMZI), experimental results demonstrate that the CRIs of the GW vary from 2.91-i13.92 to 3.81-i14.64 for transmitting wavelengths ranging from 1510 to 1590 nm. This method cannot only be used to determine the CRI of the GW optically and provide one of the fundamental parameters for designing graphene-based optic devices for communication and sensing applications, but also is adoptable in graphene-based transformation optics for determination of the CRI of the GW at other wavelengths. © 2013 Optical Society of America OCIS codes: (060.5060) Phase modulation; (160.4670) Optical materials; (310.2790) Guided waves.
References and links 1. 2.
A. K. Geim and K. S. Novoselov, “The rise of graphene,” Nat. Mater. 6(3), 183–191 (2007). R. R. Nair, P. Blake, A. N. Grigorenko, K. S. Novoselov, T. J. Booth, T. Stauber, N. M. R. Peres, and A. K. Geim, “Fine structure constant defines visual transparency of graphene,” Science 320(5881), 1308 (2008). 3. S. Thongrattanasiri, F. H. L. Koppens, and F. J. García de Abajo, “Complete optical absorption in periodically patterned graphene,” Phys. Rev. Lett. 108(4), 047401 (2012). 4. 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(7), 532–535 (2008). 5. M. Breusing, S. Kuehn, T. Winzer, E. Malić, F. Milde, N. Severin, J. P. Rabe, C. Ropers, A. Knorr, and T. Elsaesser, “Ultrafast nonequilibrium carrier dynamics in a single graphene layer,” Phys. Rev. B 83(15), 153410 (2011). 6. B. V. Cunning, C. L. Brown, and D. Kielpinski, “Low-loss flake-graphene saturable absorber mirror for laser mode-locking at sub-200-fs pulse duration,” Appl. Phys. Lett. 99(26), 261109 (2011). 7. G. K. Lim, Z. L. Chen, J. Clark, R. G. S. Goh, W. H. Ng, H. W. Tan, R. H. Friend, P. K. H. Ho, and L. L. Chua, “Giant broadband nonlinear optical absorption response in dispersed graphene single sheets,” Nat. Photonics 5(9), 554–560 (2011). 8. T. H. Han, Y. Lee, M. R. Choi, S. H. Woo, S. H. Bae, B. H. Hong, J. H. Ahn, and T. W. Lee, “Extremely efficient flexible organic light-emitting diodes with modified graphene anode,” Nat. Photonics 6(2), 105–110 (2012). 9. J. T. Kim and S.-Y. Choi, “Graphene-based plasmonic waveguides for photonic integrated circuits,” Opt. Express 19(24), 24557–24562 (2011). 10. M. Liu, X. Yin, E. Ulin-Avila, B. Geng, T. Zentgraf, L. Ju, F. Wang, and X. Zhang, “A graphene-based broadband optical modulator,” Nature 474(7349), 64–67 (2011). 11. Q. Bao, H. Zhang, B. Wang, Z. Ni, C. H. Y. X. Lim, Y. Wang, D. Y. Tang, and K. P. Loh, “Broadband graphene
#199279 - $15.00 USD Received 10 Oct 2013; revised 14 Nov 2013; accepted 15 Nov 2013; published 25 Nov 2013 (C) 2013 OSA 2 December 2013 | Vol. 21, No. 24 | DOI:10.1364/OE.21.029818 | OPTICS EXPRESS 29818
polarizer,” Nat. Photonics 5(7), 411–415 (2011). 12. Y. Wu, B. Yao, Y. Cheng, X. Liu, Y. Gong, and Y. Rao, “Hybrid graphene-microfiber waveguide for chemical gas sensing,” J. Sel. Top. Quantum Electron. 20(1), 4400206 (2013). 13. M. Bruna and S. Borini, “Optical constants of graphene layers in the visible range,” Appl. Phys. Lett. 94(3), 031901 (2009). 14. X. Wang, Y. P. Chen, and D. D. Nolte, “Strong anomalous optical dispersion of graphene: complex refractive index measured by Picometrology,” Opt. Express 16(26), 22105–22112 (2008). 15. A. Vakil and N. Engheta, “Transformation optics using graphene,” Science 332(6035), 1291–1294 (2011). 16. A. N. Grigorenko, M. Polini, and K. S. Novoselov, “Graphene plasmonics,” Nat. Photonics 6(11), 749–758 (2012). 17. L. Yuan, “Recent progress of in-fiber integrated interferometers,” Photonic Sens. 1(1), 1–5 (2011). 18. M. Oya, H. Kishikawa, N. Goto, and S. Yanagiya, “All-optical switch consisting of two-stage interferometers controlled by using saturable absorption of monolayer graphene,” Opt. Express 20(24), 27322–27330 (2012). 19. G. W. Hanson, “Dyadic Green’s functions and guided surface waves for a surface conductivity model of graphene,” J. Appl. Phys. 103(6), 064302 (2008). 20. http://refractiveindex.info/?group=CRYSTALS&material=MgF2. 21. X. He, Z. B. Liu, and D. N. Wang, “Wavelength-tunable, passively mode-locked fiber laser based on graphene and chirped fiber Bragg grating,” Opt. Lett. 37(12), 2394–2396 (2012). 22. B. Yao, Y. Wu, L. Jia, Y. Rao, Y. Gong, and C. Jiang, “Mode field distribution of optical transmission along micro fiber affected by CNTs films with complex refraction index,” J. Opt. Soc. Am. B 29(5), 891–895 (2012). 23. N. J. Horing, “Coupling of graphene and surface plasmons,” Phys. Rev. B 80(19), 193401 (2009). 24. E. H. Hwang and S. D. Sarma, “Dielectric function, screening, and plasmons in two-dimensional grapheme,” Phys. Rev. B 75(20), 205418 (2007). 25. V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, “Unusual microwave response of Dirac quasiparticles in graphene,” Phys. Rev. Lett. 96(25), 256802 (2006). 26. L. M. Tong, R. R. Gattass, J. B. Ashcom, S. L. He, J. Y. Lou, M. Y. Shen, I. Maxwell, and E. Mazur, “Subwavelength-diameter silica wires for low-loss optical wave guiding,” Nature 426(6968), 816–819 (2003). 27. Z. G. Wang, Y. F. Chen, P. J. Li, X. Hao, J. B. Liu, R. Huang, and Y. R. Li, “Flexible graphene-based electroluminescent devices,” ACS Nano 5(9), 7149–7154 (2011). 28. S. A. Mikhailov and K. Ziegler, “New electromagnetic mode in graphene,” Phys. Rev. Lett. 99(1), 016803 (2007). 29. H. Li, Y. Anugrah, S. J. Koester, and M. Li, “Optical absorption in graphene integrated on silicon waveguides,” Appl. Phys. Lett. 101(11), 111110 (2012).
1. Introduction Graphene with complex permittivity has been extensively investigated in recent years due to its attractive physical properties, such as strong conductivity [1], good transparency [2] and notable medium nonlinearity [3], etc. Recently, many physics phenomena, such as graphenebased Fermi effect [4] and saturable absorptions [5], have been discovered and a number of graphene-based photonic devices, such as mode locked laser [6], optical switch [7] and LED [8], have been reported. Especially, when graphene acts as an optical waveguide, it has great potential for applications to photonic integrated circuits [9], optical fiber communication and sensing, such as intensity modulator [10], polarizer [11], and gas sensor [12]. The determination of the complex refractive index (CRI) of the graphene-waveguide (GW) becomes essential for designing of novel graphene-based waveguide devices. However, in many reports, graphene was only regarded as boundary conditions, because it is still a challenge to determinate the complex refractive index (CRI) of the GW effectively. So far, the CRI of G-films rather than GWs could be determined by using an approximate theory [13] or the Picometrology method proposed for visible wavelengths, e.g. for 532nm and 633nm, the CRI of graphene film is 2.4-i1.0 and 3.0-i1.4, respectively [14]. As the permittivity of graphene is anisotropic and coupling light into and out the GW is still a challenge, it is very difficult to measure its CRI optically. This remains as a fundamental problem to be solved in graphene-based transformation optics and plasmonics [15,16]. In this paper, microfiber is used as an effective and convenient means to launch and collect evanescent light from the GW, and offers the flexibility of changing the contact length in a graphene/microfiber hybrid waveguide (GMHW). By this means, graphene induced optical phase alterations and attenuations are measured by microfiber-based Mach–Zehnder interferometer (MMZI) [17,18], thus the CRI of the GW in the wavelength window of
#199279 - $15.00 USD Received 10 Oct 2013; revised 14 Nov 2013; accepted 15 Nov 2013; published 25 Nov 2013 (C) 2013 OSA 2 December 2013 | Vol. 21, No. 24 | DOI:10.1364/OE.21.029818 | OPTICS EXPRESS 29819
1510nm-1590nm are calculated from the measured effective refractive index of the GMHW. It is summarized when the light propagates along the GW, the CRI of the single layer GW ng ranges from 2.91-i13.92 to 3.81-i14.64 at room temperature for this wavelength window. Such a method provides a powerful way to investigate the optical properties of graphene or other photonic thin films by spectral measurement. 2. Modeling and theoretical analysis Figure 1(a) schematically shows the configuration of the GMHW. A microfiber drawn from a single mode fiber (SMF) is attached onto the graphene film tightly. The graphene is deposited on a MgF2 substrate. The evanescent wave is launched to the GMHW from one taper of the microfiber and collected by the other taper at the output port. Here, the effective length of the GMHW would be adjusted conveniently via controlling the interaction distance between the microfiber and the graphene film. The geometry of the cross-section of the GMHW is shown in Fig. 1(b). The microfiber (blue) with diameter of 1.1μm and refractive index (RI) of n1 = 1.45 is located onto the graphene film (black), the permittivity of the GW is related to its conductivity σ as εg,eq = -σg,i/ω + iσg,r/ω [15]. Here σg,i and σg,r are the imaginary and real part of σ, ω is the light’s angular frequency, is the thickness of graphene. Here, the graphene is regarded as an ultrathin bulk waveguide with = 0.5nm (single layer in physics) in this model rather than a boundary condition when only considering losses [11]. Considering the pristine graphene at room temperature, the CRI of the GW is calculated in Fig. 1(c), by applying ω2μ0ε = n2k02. Here μ0 is the permeability, ε is the permittivity, and k0 is the wave number in free space [15,19]. The MgF2 substrate has a 3mm thickness and its RI is equal to 1.37 at ~1550nm [20]. By adopting the Finite Element Method (FEM) to multilayer waveguides, the effective RI of the GMHW for the fundamental mode is numerically calculated and given in Fig. 1(d). For the wavelength range from 1300nm to 1600nm, the real part of the effective index decreases from 1.405 to 1.339, while the imaginary part increases from 3.1 × 10−4 to 3.8 × 10−4. For the microfiber (D = 1.1μm) without GW attached, the effective RI for the wavelength band 1510nm-1590nm is a constant 1.426 [21]. We simulate the sectional field distribution of the transmitting light along the GMHW with λ = 1550nm, shown in Fig. 1(e). According to Fig. 1(e), we observe that the GW attracts a large portion of evanescent waves to transmit along it. Figure 1(f) shows the intensity distribution for y-axis according to Fig. 1(e). The curve is discontinuous at the graphene layer, which is caused by the inherent optical absorption of graphene [2]. It is clear that once the neff of the GMHW could be measured, the CRI of the pure GW can be numerically calculated from the experimental results. First, the process of evanescent wave couples from microfiber to the GW is investigated. The evanescent wave transmitting along microfiber is determined by its eigenequations [22,23], once the microfiber coupled the intrinsic graphene, the transmitting mode is not only HE11 any longer [24,25]. Considering boundary condition n × (H1-H2) = σE and adopting the model shown in Fig. 1(a), we write the light just influenced by the GW as E(x,y,z) = (exx + eyy + ezz)ejßzejωt, where we assume that the GW is located at the y-z flat and the wave vector is in the direction of z, here β = neff(ω)(ω/c). Referring neff(ω) = neffRE(ω) - jneffIM(ω), the transmitting wave could be written as E(x,y,z) = (exx + eyy + ezz) exp[ω(-neffIM × z/c)] exp[jω(neffRE × z/c) + t]. For transmitting waves, exp[ω(-neffIM × z/c)] is the attenuation coefficient and exp[jω(neffRE × z/c) + t] is the phase coefficient. Here ω = 2πc/λ. The neffRE induced phase change and the neffIM induced attenuation change (z = 20μm, λ = 1550nm) are simulated in Figs. 2(a) and 2(b). Figures 2(c)–2(e) show the simulated spatial distributions of electric fields transmitting along microfiber, microfiber attached MgF2 and the GMHW respectively, for z = 20μm and λ = 1550nm as well. Larger attenuation, less compressed phase distribution in the GMHW can be obtained in the simulations. Moreover, in Fig. 2(e), TE polarized surface plasmon on the GW is obvious, which has been predicted in [24].
#199279 - $15.00 USD Received 10 Oct 2013; revised 14 Nov 2013; accepted 15 Nov 2013; published 25 Nov 2013 (C) 2013 OSA 2 December 2013 | Vol. 21, No. 24 | DOI:10.1364/OE.21.029818 | OPTICS EXPRESS 29820
Fig. 1. (a) Schematic diagram of the structure for the light propagating along the GMFW (Red: SMF, White: Monolayer graphene, Cyan: MgF2 substrate). The orange arrows show the transmitting direction of the evanescent waves. (b) Geometry of the cross-section of the GMHW. (c) The theoretical CRI of the GW. (d) The effective CRI of the GMHW (Blue curve: real part, red dashed: imaginary part). (e) Simulated field intensity distribution of the GMHW. (f) The field intensity of the Microfiber/GW distributed along y-axis according to (e).
Fig. 2. For λ = 1550nm: (a) Correlation between neffRE and transmitting phase at z = 20μm. (b) Correlation between neffIM and attenuation at z = 20μm. (c) (d) (e): The 3-D distributes of the electric field intensity along the microfiber, the microfiber/MgF2 and the GMHW, respectively.
#199279 - $15.00 USD Received 10 Oct 2013; revised 14 Nov 2013; accepted 15 Nov 2013; published 25 Nov 2013 (C) 2013 OSA 2 December 2013 | Vol. 21, No. 24 | DOI:10.1364/OE.21.029818 | OPTICS EXPRESS 29821
3. Experimental verification based on specially designed MMZI In order to measure the neff of the GMHW, a MMZI with one arm with the GMHW is constructed, shown in Fig. 3(a). A SMF with a length of ~50cm is used as the lead-in fiber, and a polarization controller is adopted to optimize the incident polarization mode to be TEpolarized for the graphene layer. Through a 50:50 Coupler, light is divided into 2 interfering arms. For the arm 1, the light is guided by SMF, here a tunable attenuator is used to match the intensities of the arm1 and arm 2. The length difference of the two arms L2-L1 is ~5cm (arm 2 is longer). For the arm 2, the SMF is drawn into a microfiber with a diameter of ~1.1μm and a length of ~17cm [26], the waist part of the microfiber is attached onto the graphene film supported by the MgF2 substrate (40mm × 40mm), shown in Fig. 5. The output light is collected by the SMF at the other end of the microfiber. Finally the arm 2 and arm 1 are combined via a 50:50 coupler and sent into an optical spectrum analyzer (OSA, Micron Si 720, USA). When broadband light lunches into the MMZI, the output is a resonant spectrum with interference fringes. The GMHW in experiment is shown in Fig. 3(b) schematically. The tapers of the microfiber are fixed on the docks and the waist of the microfiber is pendulous naturally. The graphene deposited on the MgF2 substrate (G/MgF2) is fixed on a translation stage. When the microfiber is close to graphene, under the effects of Van der Waals force and electrostatic force, they can touch tightly, forming a stable structure. Figure 3(c) shows the SEM photograph of the microfiber attached onto the GW. In our experiment, graphene films were grown on Cu foils (Alfa Aesar, No. 13382) by chemical vapor deposition (CVD) in the way described in [27]. To transfer the graphene onto the objective MgF2 substrate, PMMA was spin-coated on the surface of graphene/Cu foil and then the underlying Cu foil was etched with 1M FeCl3 solution. Subsequently, the PMMA/graphene was washed in DI water several times and transferred onto the MgF2 substrate which had been ultrasonically cleaned in sequence by acetone, ethanol and DI water. Then it was allowed to dry at room temperature for overnight and baked at 180°C for 10 min. Finally, the PMMA was removed by acetone. The Raman spectrum of the GW is shown in Fig. 3(d). A very weak D peak, a ~0.23 G-to-2D intensity ratio, and a 2D-peak with a full width at half maximum of ~36.3 cm–1, suggest that the graphene adopted in this experiment is monolayer with high quality.
Fig. 3. (a) Setup of the GMHW-based MMZI. (b) The experimental details of the GMHW (c) SEM of the GMHW. (d) Raman spectrum of the GW on the MgF2 substrate, the blue curve indicates the G/MgF2 and the red line indicates the MgF2 substrate only.
The phase alteration of the microfiber by the GW in one arm of the MMZI would cause the spectral shift. For the MMZI with only microfiber in arm 2 (without graphene attached), the location of any resonance dip λd is shown in Eq. (1). Here L1 is the length of arm 1, L2 is the length of arm 2, N is a natural number related to neffRE and (L2-L1), vf is the phase velocity
#199279 - $15.00 USD Received 10 Oct 2013; revised 14 Nov 2013; accepted 15 Nov 2013; published 25 Nov 2013 (C) 2013 OSA 2 December 2013 | Vol. 21, No. 24 | DOI:10.1364/OE.21.029818 | OPTICS EXPRESS 29822
of light in fiber. Once the GW is attached, the phase velocity in the GMHW zone is vG rather than vf, so that the location of the spectral dip would shift to λd’, as shown in Eq. (2). As the phase velocity of electromagnetic waves is v = c/nRE, when the neffRE decreases, the dip location λd has a blue shift. The shift Λ = λd’-λd is shown in Eq. (3). Here nRE = neffREGMHW neffREMF