We investigate the Bloch mode conversion of surface plasmon polaritons in a periodic array of graphene pairs with each consisting of two separated parallel graphene sheets. The employment of graphene pair as a unit cell in the array yields two Bloch modes belonging to different bands. By periodically modulating the permittivity of dielectrics between graphene along the propagation direction, the interband transitions occur and the modes will alternatively couple to each other, similar to traditional Rabi oscillations in quantum systems. The indirect Rabi oscillations can also be observed through introducing transverse modulation momentum. The period of Rabi oscillations can be optimized by taking advantage of the flexible tunability of graphene. The study suggests that the structure have applications in optical switches and mode converters operating on deep-subwavelength scale.
© 2015 Optical Society of America
Rabi oscillations usually refer to periodic state transitions between two atomic energy levels driven by an electromagnetic wave with the frequency tuned to the energy gap, leading to the periodic emission and reabsorption of photons [1, 2 ]. Such phenomena of alternate conversion of states coupled to external periodic perturbations have also been achieved in semiconductors and condensed matter systems [3, 4 ]. In the past decade, optical Rabi oscillations have been demonstrated in multimode waveguides, photonic crystals and waveguide arrays [5–7 ]. The optical Rabi oscillation, analogous to its initial conception, refers to the alternate conversion between distinct guided modes as the waveguides undergo a longitudinal periodic modulation of permittivity. As for waveguide arrays, the Bloch modes belonging to different bands will convert to each other and experience an alternate energy exchange during propagation. The effect may find applications in optical switches and spatial mode converters [8–10 ]. Recently, efficient Rabi oscillations of surface plasmon polaritons (SPPs) have been realized in metal-dielectric waveguides within a subwavelength oscillation period . Similar to metals in visible ranges, graphene can support SPPs in terahertz and far infrared regimes. Compared to the SPPs in metal-dielectric waveguides, the SPPs in graphene possess huge mode localization, low propagation loss and flexible tunability [12–15 ]. By taking advantage of the unique properties, the graphene constructed waveguides and arrays could provide a new platform to investigate Rabi oscillations.
In this work, we propose the graphene-pair arrays (GPAs) composed of periodically arranged graphene pairs to realize Rabi oscillations of SPPs. A graphene pair can support the symmetric and antisymmetric SPP modes [16, 17 ]. Consequently, the GPAs allow two kinds of collective Bloch modes belonging to different bands of the dispersion relation, depending on the symmetry of SPP modes in each individual graphene pair. As the permittivity of the dielectric between graphene undergoes a periodic modulation in the propagation direction, the interband transition could occur. When the permittivity modulation carries a proper transverse momentum, indirect interband transition can be stimulated as well. Here we realize efficient Rabi oscillations with a deep subwavelength period of a few hundred nanometers. The period of Rabi oscillations can also be controlled by tuning the incident Bloch momentums, operation wavelength and chemical potential of graphene.
2. Bloch modes in GPAs
We start by investigating the eigen Bloch modes in the GPAs without modulation. As shown in Fig. 1(a) , the GPAs are composed of periodically arranged graphene pairs separated by dielectrics. The period of the GPAs is denoted by d = d 1 + d 2 with d 1 being the spacing of graphene in each pair. The surface conductivity of graphene σg(λ, μc, τ, T) is modeled by the Kubo formula [18, 19 ], which is determined by the wavelength (λ), chemical potential of graphene (μc), relaxation time (τ), and temperature (T). We initially choose λ = 10 μm and μc = 0.15 eV. The relaxation time is τ = 1 ps at room temperature T = 300 K . The relative permittivity of the dielectric medium is assumed as εd = 2.13 .
We only consider transverse magnetic (TM) polarized SPPs propagate along z direction. According to Maxwell’s equations and the Bloch theorem [22, 23 ], we can obtain the dispersion relation of the Bloch modes
Figure 1(b) shows the band structure of the Bloch modes in the GPAs. Note that the band structure here denotes the diffraction relation between the propagation constant or effective index of the Bloch modes and the Bloch momentum [24, 25 ]. There are two bands in the diagram, corresponding to two kinds of Bloch modes in the array. The mode wavelength and propagation distance of the Bloch modes are given by λp = λ/Re(n eff) and L p = [2k 0∙Im(n eff)]−1 [26, 27 ]. As shown in Figs. 1(b) and 1(c), the modes of band 1 have shorter mode wavelengths and longer propagation distances by comparing with the modes of band 2. Particularly for band 1, the modes in the Brillouin zone center have smaller λp and larger L p than those at the edges. Concerning band 2, in contrast, the modes at the Brillouin zone edges have shorter mode wavelengths and larger propagation distances. As the period d of the GPAs increases, the coupling of SPPs between adjacent graphene pairs becomes weaker. Consequently, both the real and imaginary parts of the effective indices become flatten (dotted curves), which are more remarkable for band 2.
The normalized magnetic field distributions of the Bloch modes are illustrated in Fig. 2 . We denote the mode as Φm ,φ with m and φ being the band index and Bloch momentum. As each graphene pair can support two kinds of SPP mode, i.e., the symmetric and antisymmetric modes. The Bloch modes of band 1 are formed by the coupling of the symmetric SPP mode in each graphene pair with a phase difference of φ between adjacent periods. As shown in Figs. 2(a)-2(c), the mode profiles between adjacent periods are in phase as φ = 0 and out of phase as φ = π. As φ = π/2, the modes are out of phase in every two periods. The Bloch modes in band 2 are formed by the coupling of antisymmetric SPP mode in each graphene pair. Figures 2(d)-2(f) show the mode profiles in band 2 as φ = 0, π/2 and π, respectively. The symmetry of the modes in band 2 is opposite to that of the corresponding modes of band 1. The symmetry of the modes is tightly related with the conversion efficiency of Rabi oscillations, which will be discussed in the following.
3. Rabi oscillations of SPP Bloch modes
Now we consider the Rabi oscillations between the Bloch modes in bands 1 and 2 stimulated by the permittivity perturbation. The modulation of permittivity applied to the permittivity dielectrics between graphene is given by28]29, 30 ]28]. As the transverse modulation momentum is ϕ = 0, only direct interband transition between the modes with φ2 = φ1 will take place. While for nonzero modulation momentum (ϕ ≠ 0), indirect interband transition between the Bloch modes with φ2 − φ1 = ϕ can be stimulated.
Figure 3 illustrates the coupling coefficient |M12| versus the incident Bloch momentum φ1 and momentum difference Δφ = φ2 – φ1. The incident Bloch momentum can be introduced by controlling the phase difference of SPPs in adjacent graphene pairs . The momentum difference between the modes of different bands is generated by the transverse modulation momentum by choosing ϕ = Δφ . In Fig. 3(a), the transverse modulation is odd and reads Δε(x) = sin(2πx/d). In respect of the direct interband transition as Δφ = 0, M12 has the largest modulus at the Brillouin zone center φ1 = 0 and decreases from the center to the edges φ1 = ± π. While for indirect interband transition Δφ = π, the maximum of the coupling coefficient locates at the Brillouin zone edges. In Fig. 3(b), the transverse modulation is even and reads Δε(x) = cos(2πx/d), the coupling coefficients for both direct (Δφ = 0) and indirect (Δφ = π) interband transitions reach their maxima near φ1 = ± π/2. As φ1 = 0 and ± π, the coupling coefficient |M12| = 0, leading to the forbiddance of Rabi oscillations. Since the modes possess opposite symmetries for φ1 = 0 and ± π as shown in Fig. 2, the overlap integration of M12 in Eq. (4) will vanish under an even transverse modulation, and Rabi oscillations are forbidden.
For other incident Bloch momentum 0 < |φ1| < π, due to the broken of mode symmetry, the coupling coefficient |M12| is nonzero. Thus Rabi oscillations can be yielded under both odd and even transverse modulations. Figures 3(c) and 3(d) show the coupling coefficients as the transverse modulation is odd for Δε(x) = sin(2πx/d) and even for cos(2πx/d). Both situations for Δφ = ± π/2 are considered. Note that |M12| is larger under the odd modulation than the even by more than one order of magnitude. The odd transverse modulation is preferred over the even to reduce the period of Rabi oscillations. As Δφ = π/2 shown in Fig. 3(c), |M12| is largest near φ1 = − π/2 and smallest near φ1 = π/2. While in Fig. 3(d), the coupling coefficient is largest near φ1 = 0 and approaches zero near φ1 = ± π/2. Thus by introducing certain transverse modulations, the Rabi oscillations can be controlled arbitrarily.
To validate the theoretical analysis, we also perform numerical simulations by using COMSOL Multiphysics. The results are shown in Fig. 4 . In the calculation, graphene is modeled as the surface current boundary condition [31–33 ]. Figures 4(a) and 4(b) illustrate the analytical and numerical field evolution of the direct interband transition between the Bloch modes Φ1,0 and Φ2,0 under the odd modulation of Δε(x) = sin(2πx/d). They agree fairly with each other. The analytical and numerical periods of Rabi oscillations (denoted by the horizontal dotted lines) are 0.385 μm and 0.380 μm, respectively. In Fig. 4(c), the mode Φ1,0 is injected from the end of the GPAs under the even modulation of Δε(x) = cos(2πx/d). It shows that Rabi oscillations do not occur under even modulation for φ1 = 0, which confirms the analysis in Fig. 3(b). Figures 4(d) and 4(e) show the field evolution of the indirect interband transitions of Φ1,0 ↔ Φ2,π and Φ1,π ↔ Φ2,0 under the odd modulation of Δε(x) = sin(2πx/d). The periods of Rabi oscillations are 0.48 μm and 0.45 μm, which are larger than those of direct interband transitions shown in Fig. 4(b). It should be mentioned that the propagation loss of SPPs might hinder the observation of Rabi oscillations. In order to implement complete mode conversions, one can reduce the period of Rabi oscillations to combat with negative influence of the propagation loss.
Figure 5 illustrates the influence of other parameters on the period of Rabi oscillations as δε = 0.3 and Δε(x) = sin(2πx/d). In Fig. 5(a), the oscillation period as a function of the spacing of each graphene pair and the period of the GPAs is investigated. As d 1 is fixed, the period of Rabi oscillations increases as d increases. The coupling of SPPs between adjacent graphene pairs becomes weaker as d increases, making the oscillation period increase. For fixed d, the period of Rabi oscillations decreases as d 1 increases. It shows that the GPAs with close d 1 and d 2 can benefit the occurrence of Rabi oscillations. In Fig. 5(b), we investigate the influence of incident Bloch momentum and momentum difference on the period of Rabi oscillations as d 1 = 20 nm and d = 80 nm. Due to the symmetry of the band structure, the diagram is central symmetric, that is, L c(φ1, Δφ) = L c(− φ1, − Δφ). As φ1 is fixed, L c increases as Δφ increases. The oscillation period is smallest for direct transition as φ1 = 0, Δφ = 0 and largest for indirect transition as φ1 = 0, Δφ = ± π. Figure 5(c) illustrates the influence of incident wavelength and chemical potential of graphene on the period of Rabi oscillations. The oscillation period increases as λ or μc increases. In terahertz range, the intraband transition of electrons in graphene dominates, leading to a Drude-like surface conductivity of σg = ie 2μc/[πħ 2(ω + iτ−1)] . The transverse decay constant of SPPs in a single-layer graphene is κsp = 2iεdk 0/(σgη0) , thus we have Re(κsp) = 8cεdħ 2π3/(η0μcλ2 e 2) with c, e, ħ being the speed of light in vacuum, the electron charge and reduced Planck constant, respectively. Re(κsp) decreases as λ or μc increases, making the SPP mode more weakly confined to the surface of graphene. Consuquently, the contribution of the mode products in a period to the overlap integration of M12 in Eq. (4) becomes less, resulting in the increase of oscillation period.
In terms of the experimental implementions, the graphene-pair arrays can be fabricated as follows . Firstly, a dielectric film with a thickness of nanometers is coated on the substrate by using PLD or CVD method. Then a graphene sheet, usually grown on the cooper foil, can be transferred onto the dielectric film in the Fe(NO3)3 solution. After the graphene is adhered on the first dielectric film, we can repeat the processes of coating dielectric films and transferring graphene until the desirable layers are achieved. The thickness of dielectric is controllable during the coating process in order to form the pair array structure. The longitudinal permittivity modulation could be realized by holographic recording technique with a two-beam interference setup. As the dielectric material is photorefractive, the required variation of the permittivity can be controlled by the interference intensity which is recorded by the photorefractive dielectrics [7, 36 ]. On the other hand, the transverse modulation of the permittivity is realizable by coating multiple dielectric layers with distinct permittivity according to the desirable permittivity profile between the graphene sheets.
In conclusion, we investigate the Rabi oscillations of SPP Bloch modes in GPAs. The modes belonging to different bands are formed by the coupling of symmetric and antisymmetric SPP modes in individual graphene pairs. The interband transitions can be stimulated as the dielectric permittivity undergoes a longitudinal periodic modulation, leading to the occurrence of Rabi oscillations. The Bloch modes will convert to each other alternatively during propagation. As the modulation carries a transverse momentum, Rabi oscillations will occur as well but accompanied with indirect interband transitions. Due to the opposite symmetries of the Bloch modes at the center and edges of the Brillouin zone, Rabi oscillations are forbidden under even transverse modulation. The period of Rabi oscillations can be squeezed into a deep subwavelength scale of a few hundred nanometers. The oscillation period with direct transition is smaller than that with indirect transition. It can be further reduced by introducing proper modulation functions or optimizing the geometry of GPAs. To reduce the influence of propagation loss, we can reduce the oscillation period by choosing shorter operation wavelength and smaller chemical potential of graphene. The study may find applications in optical switches, couplers, and spatial mode converters.
This work is supported by the 973 Program (No. 2014CB921301), the National Natural Science Foundation of China (NSFC) (Nos. 11304108 and 11104095), and the Specialized Research Fund for the Doctoral Program of Higher Education of China (No. 20130142120091).
References and links
1. I. I. Rabi, “On the process of space quantization,” Phys. Rev. 49(4), 324–328 (1936). [CrossRef]
3. K. Schuh, J. Seebeck, M. Lorke, and F. Jahnke, “Rabi oscillations in semiconductor quantum dots revisited: influence of LO-phonon collisions,” Appl. Phys. Lett. 94(20), 201108 (2009). [CrossRef]
4. M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. Hall, M. J. Holland, J. E. Williams, C. E. Wieman, and E. A. Cornell, “Watching a superfluid untwist itself: recurrence of Rabi oscillations in a Bose-Einstein condensate,” Phys. Rev. Lett. 83(17), 3358–3361 (1999). [CrossRef]
7. K. Shandarova, C. E. Rüter, D. Kip, K. G. Makris, D. N. Christodoulides, O. Peleg, and M. Segev, “Experimental observation of Rabi oscillations in photonic lattices,” Phys. Rev. Lett. 102(12), 123905 (2009). [CrossRef] [PubMed]
8. C. N. Alexeyev, T. A. Fadeyeva, B. P. Lapin, and M. A. Yavorsky, “Generation of optical vortices in layered helical waveguides,” Phys. Rev. A 83(6), 063820 (2011). [CrossRef]
9. G. K. L. Wong, M. S. Kang, H. W. Lee, F. Biancalana, C. Conti, T. Weiss, and P. St. J. Russell, “Excitation of orbital angular momentum resonances in helically twisted photonic crystal fiber,” Science 337(6093), 446–449 (2012). [CrossRef] [PubMed]
10. H. Hu, K. Wang, H. Long, W. Liu, B. Wang, and P. Lu, “Precise determination of the crystallographic orientations in single ZnS nanowires by second-harmonic generation microscopy,” Nano Lett. 15(5), 3351–3357 (2015). [CrossRef] [PubMed]
14. A. N. Grigorenko, M. Polini, and K. S. Novoselov, “Graphene plasmonics,” Nat. Photonics 6(11), 749–758 (2012). [CrossRef]
15. C. Qin, B. Wang, H. Long, K. Wang, and P. Lu, “Bloch mode engineering in graphene modulated periodic waveguides and cavities,” J. Opt. Soc. Am. B 32(8), 1748–1753 (2015). [CrossRef]
16. B. Wang, X. Zhang, X. Yuan, and J. Teng, “Optical coupling of surface plasmons between graphene sheets,” Appl. Phys. Lett. 100(13), 131111 (2012). [CrossRef]
17. A. Auditore, C. de Angelis, A. Locatelli, and A. B. Aceves, “Tuning of surface plasmon polaritons beat length in graphene directional couplers,” Opt. Lett. 38(20), 4228–4231 (2013). [CrossRef] [PubMed]
18. N. M. R. Peres, “Colloquium: the transport properties of graphene: an introduction,” Rev. Mod. Phys. 82(3), 2673–2700 (2010). [CrossRef]
20. C. H. Gan, “Analysis of surface plasmon excitation at terahertz frequencies with highly doped graphene sheets via attenuated total reflection,” Appl. Phys. Lett. 101(11), 111609 (2012). [CrossRef]
21. Y. Sun, Z. Zheng, J. Cheng, J. Liu, J. Liu, and S. Li, “The un-symmetric hybridization of graphene surface plasmons incorporating graphene sheets and nano-ribbons,” Appl. Phys. Lett. 103(24), 241116 (2013). [CrossRef]
22. A. Yariv and P. Yeh, Photonics: Optical Electronics in Modern Communication (Academic, 2006).
23. B. Wang, X. Zhang, F. J. García-Vidal, X. Yuan, and J. Teng, “Strong coupling of surface plasmon polaritons in monolayer graphene sheet arrays,” Phys. Rev. Lett. 109(7), 073901 (2012). [CrossRef] [PubMed]
26. S. A. Maier, Plasmonics: Fundamentals and Applications (Springer Verlag, 2007).
29. K. Okamoto, Fundamentals and Applications of Optical Waveguides (Academic, 2006).
30. Z. Yu and S. Fan, “Complete optical isolation created by indirect interband photonic transitions,” Nat. Photonics 3(2), 91–94 (2009). [CrossRef]
31. A. Y. Nikitin, F. Guinea, and L. Martin-Moreno, “Resonant plasmonic effects in periodic graphene antidot arrays,” Appl. Phys. Lett. 101(15), 151119 (2012). [CrossRef]
36. V. Gericke, P. Hertel, E. Krätzig, J. P. Nisius, and R. Sommerfeldt, “Light-induced refractive index changes in LiNbO3: Ti waveguides,” Appl. Phys. B 44(3), 155–162 (1987). [CrossRef]