This article is part of the series Low dimensional electric and photonic materials and devices.

Open Access Nano Express

Mode manipulation and near-THz absorptions in binary grating-graphene layer structures

Haiming Yuan12, Hanning Yang12, Pengzi Liu12, Xiangqian Jiang12 and Xiudong Sun12*

Author Affiliations

1 Department of Physics, Harbin Institute of Technology, Harbin 150001, China

2 Key Lab of Micro-Optics and Photonic Technology of Heilongjiang Province, Harbin 150001, China

For all author emails, please log on.

Nanoscale Research Letters 2014, 9:90  doi:10.1186/1556-276X-9-90


The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/9/1/90


Received:13 December 2013
Accepted:28 January 2014
Published:21 February 2014

© 2014 Yuan et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.

Abstract

The excitation and absorption properties of grating coupled graphene surface plasmons were studied. It was found that whether a mode can be excited is mainly determined by the frequency of incident light and the duty ratio of gratings. In the structure consisting graphene bilayer, a blueshift of the excitation frequency existed when the distance between neighbor graphene layer were decreased gradually. In graphene-grating multilayer structures, a strong absorption (approximately 90% at maximum) was found in near-THz range.

Keywords:
Graphene surface plasmons; Excitation condition; Multilayer structure absorption

Background

Recently, a lot of work has been done based on graphene due to its unique properties in electric, magnetic, thermal, etc. [1-3]. Graphene is carbon atoms arranged in a two-dimensional honeycomb lattice, in which the electrons behave like massless Dirac fermions with linear dispersion [4,5]. Graphene has strong plasmonic effects which can be modified by gating, by doping, and so on [2]. A controllable optical absorption was also found in structured graphene [6,7].

Up to date, the graphene is modeled usually to be an extremely thin film with a conductivity σ, which consists of both intraband and interband from Kubo formula [7-9]. The intraband conductivity with Drude type plays a leading role when ℏω/μc was small [10]. Both transverse electric (TE) and transverse magnetic (TM) have the dispersion relations at monolayer graphene with dielectric materials on two sides [10-12]. In other words, the charge carriers coupling to electromagnetic waves will produce a new surface wave, namely graphene surface plasmons (GSPs).

In the previous works, many numerical approaches were used to study the structured graphene, for example the finite element method (FEM) [13], finite difference time domain (FDTD) [14], and others [6,15]. A strong plasmonic response of graphene has been demonstrated in a square-wave grating with a flat graphene on top [15]. In which, the graphene-based plasmon response lead to a 45% optical absorption. In a periodic array of graphene ribbons, remarkably large GSPs result in prominent optical absorption peaks [13]. In multilayer graphene, the absorption spectrum can be decomposed into subcomponents [6], which is helpful in understanding the behavior of GSP coupling.

In this paper, we studied the binary grating bounded by graphene on both sides. The rigorous coupled wave analysis (RCWA) [16,17] was used the first time as we know to characterize the graphene-containing periodic structures. The excitation condition and excitation intensity seemed to be influenced by the grating constant, duty ratio and the distance between the graphene layers. When introducing more graphene layers into the structure periodically, a strong absorption band was found in the near-THz range.

Methods

Electromagnetic mode of binary grating-graphene

Previous research has shown that the conductivity of graphene came from the contribution of intraband and interband [18-22]. The interband conductivity tends to be ignorable when ℏωμc (see [10]). Then the intraband conductivity can be expressed as [23]

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M1">View MathML</a>

(1)

where μc is the chemical potential, relating to the electron density. Equation 1 became a Drude type when μc/kBT ≫ 1, i.e.

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M2">View MathML</a>

(2)

First, we studied graphene layer positioned at the interface of two half planes with permittivity ε1 and ε2, respectively. The electric induced current on graphene layer resulted in a magnetic field difference, which led to the coupled GSP on graphene layer. Using Maxwell equation and boundary condition, GSP modes were proved to existed for both TE and TM polarization [12,23-25]. For TE mode, the dispersion relation was as follows:

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M3">View MathML</a>

(3)

and for TM mode it became

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M4">View MathML</a>

(4)

Because the imaginary part of conductivity (2) was positive, no solution of Equation 3 was found in real, which meant the TE mode GSP could not be excited.

For TM mode, put Equation 2 into Equation 4, we found

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M5">View MathML</a>

(5)

Here, we defined neff = β/k0 = βc/ω as the effective index of GSP. After making a transformation of (ω, neff) → (ω, β), the dispersion relations were obtained and plotted in Figure  1. The wave vector was normalized by kΛ0 = 2π/λ0, λ0 = 1 μm. As a local mode, GSP modes were same as the surface plasmon polaritons (SPPs). They cannot be excited directly from the air. And in our work, gratings were used to provide an external wave vector to match the phase condition.

thumbnailFigure 1. Dispersion relations of graphene surface plasmons (GSPs) on monolayer graphene with different material on two sides. Here, we use the graphene parameters of μc = 0.2 eV, τ-1 = 1 meV.

Rigorous coupled wave analysis in graphene-containing structures

In Figure  2a, we used h to be the depth of grating (thickness of gratings). The h was also the distance between two graphene layers. In multilayer structures of Figure  2b, 2 h was the longitudinal period. The structures were designed to only contain two kinds of interfaces.

thumbnailFigure 2. Binary grating graphene structures. (a) The bilayer graphene structure. (b) The multilayer graphene structure. h is the grating layer thickness. Λ is the period of grating. L1 is the width of dielectric with ε1. L(L2) is the width of dielectric with ε2. The duty ratio is f2 = L/Λ, and f1 = 1 - f2. In this paper, we simply set ε1 = 1 and ε2 = 4.

In common, the conventional RCWA based on the Floquet's theorem [26] was unable to be used for the graphene-containing structures as the electric field will induce a current with current density J = σE, while graphene was included.

In RCWA, the field was expanded into the form of

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M6">View MathML</a>

(6)

So the current density J can also be expended to the sum of spatial harmonics with different wave vector components. To obtain the reflection, transmission, absorption, field distribution, and other optical properties of such structures as shown in Figure  2, a nonzero item must be included in the boundary condition of Hy field considering the induced current,

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M7">View MathML</a>

(7)

According to the principle of superposition, Hy will also be continuous at the interface if each spatial harmonics subcomponent satisfied the boundary conditions independently,

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M8">View MathML</a>

(8)

in which n was the order, ± in subscripts represented approaching to y0 from two different directions. After the modification on the RCWA program, we can utilize it to deal with the graphene-containing structures.

Results and discussion

Phase matching condition

For a structure with a binary grating bounded by graphene layers on two sides shown in Figure  2a, the attenuated total reflection spectrum is plotted in Figure  3 using the modified RCWA method when it was illuminated normally. A set of absorption peaks each corresponding to a GSP mode was shown in blue solid line. From left to right, each peak corresponding to a GSP mode ordered with 1, 2 …

thumbnailFigure 3. Attenuated total reflection of the structure in Figure2a. Λ = 11 μm, L = 10 μm, so the duty ratio is 10/11. Each of the absorption peaks (on blue solid line) corresponds to a GSP mode.

In the structure shown in Figure  2a, there exist two kinds of interfaces, i.e., ε1-graphene-ε1 and ε1-graphene-ε2. When GSP is propagating along the interfaces, the phase shifts were φ1 and φ2 for the two kinds of interfaces, respectively. δ was the total phase loss considering two abrupt phase changes when GSP propagates across the joints between the two kinds of interfaces in a grating period. At the excitation frequency, the phase change in a grating period should satisfy the relation

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M9">View MathML</a>

(9)

which was known to be the phase matching conditions [27,28]. In Equation 9, N is the integer and can be rewritten as

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M10">View MathML</a>

(10)

where f2 = L/Λ and f1 = 1 - f2, β1 and β2 were the wave vectors of GSP on two kinds of interfaces, respectively. When N was a nonnegative integer, the GSP mode could be excited, and N can be defined as the order of surface modes. The resonant frequencies can be obtained both from absorption spectrum in Figure  3 and theoretically from Equation 10 (δ = 0). They were given in Table  1 and agreed well for high order modes. But for low order modes, some deviations existed between numerical and theoretical caused by the coupling of GSPs on two graphene layers.

Table 1. The resonant frequency of different orders

The field distributions of orders 1 and 2 of the structure in Figure  2a were given in Figure  4. It was indicated that the GSP field distributions had nodes as standing wave because the GSP modes propagating in two directions were excited simultaneously.

thumbnailFigure 4. Field distributions of |Ey| of first (11.9 meV, left) and second (16.7 emV, middle) order GSP modes. The last figure was real part of Ey of second order.

Duty ratio and stand wave interference

By using the modified RCWA, the absorption spectrum was obtained in Figure  5 when varying f, where f = φ1/(φ1 + φ2), φ1, φ2 had the same meaning as Equation 9. From the discussion above, when the phase match conditions were satisfied, GSPs could be coupled and absorption peaks should appear. But some break points appeared as shown in Figure  5, which meant that GSPs cannot be excited at this frequency and this duty ratio and it was reasonable to believe that there must be other process related to f of controlling the excitation of GSP.

thumbnailFigure 5. Absorption spectra for duty ratio vs the frequency fixing the light path of grating period.

From the field distributions in Figure  4, each corner of the grating was a singular point of field and these scatting points became the sources of surface wave, as Figure  6 shown. In periodic, we only need to consider the scatting in one period, i.e., A and B. Each scatting point will couple to two GSP modes propagating in two directions. So, the field can be expressed in four terms, which is [28,29]

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M11">View MathML</a>

(11)

thumbnailFigure 6. Corners of grating will become the scatting points of the incident light which was the source of GSPs. These scatting points can be divided into two kinds due to the geometric symmetry, which is A and B. Each scatting point will scatter into two GSP modes propagating in two directions (blue and green).

First two terms were GSP excited by one set of points (A in Figure  6) with two propagating directions (blue and green) and the last two terms were that from another set of points (B in Figure  6), where x0 is the distance of A and B in the form of light path (k0x0 = L1β1 = φ1 = (φ1 +φ2)f = 2πNf). Because in real space, different interfaces (ε1/ε1 and ε1/ε2) had different propagating constants, the expression might be complex. Here, the light path of x was used. It is found that scatting points A and B had a phase difference of π. This was caused by the different geometric symmetries. From Equation 11, when sin(k0x0/2) = 0, i.e., f = m/N ( m = 0, 1, …, N), field amplitude F would always be 0, which meant that the field cannot be excited. It was a cancelation process of two sets of standing waves that are coherent. So, for GSP mode of N, N + 1 of none absorption points appeared.

Coupling of GSPs on different graphene layers and resonant frequency shift

From Table  1, we can see that for higher order modes, the consistency between the theory and the numerical results from RCWA was better than that of the lower order modes. It was because the structure consists of bilayer of graphene and there could be interaction between GSP modes on neighbor graphene layers determined by the depth of the grating. In order to understand the behavior of GSPs coupling, in Figure  7, the absorption spectra were given as a function of the grating deepness h. A blueshift of absorption peaks was found when the grating became thin. The oscillator model is used to describe this phenomenon of spectrum blueshift [30,31].

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M12">View MathML</a>

(12)

thumbnailFigure 7. The absorption spectrum for various grating thickness.

In Equation 12, κ(n, h, ∆θ) was the coupling coefficient and n, h, and ∆θ were order of GSP mode, thickness of grating, and phase difference of GSPs on two graphene layers, respectively. Essentially, the GSPs were surface waves so they interact with each other via evanescent interactions, and the coupling intensity decayed exponentially with h increasing. For fixed h, the lower order modes had larger skin depth (stronger coupling intensity) than the higher orders; then, the stronger coupling resulted in a large spectra shift. The phase difference of ∆θ also had affection to the absorption frequencies. However, in our case, the wavelength (15 meV ~ 82.8 μm) was much larger than the thickness of grating layer (h = 10 μm), it is reasonable to assume ∆θ is approximately 0. This can also be obtained clearly from the field distribution in Figure  4 that the electric fields on upper and lower graphene layers oscillated synchronously. This conclusion can still hold in multilayer graphene-grating structures. Finally, κ(n, h, ∆θ) ∝ e-hq(n), where <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M13">View MathML</a>.

Suppose the solution of having the form of xup = xdown = x0e-iωt (no phase difference between GSP on neighbor layers), it is found that the resonant frequency became

<a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M14">View MathML</a>

(13)

When h was small (h < 4 μm), the larger κ(n, h, ∆θ) ∝ e-h was the larger shift of resonant frequency would be. And obviously, κ(n, h, ∆θ) was approaching 0 rapidly when h was large enough, which meant that the resonant frequency became a stable value of <a onClick="popup('http://www.nanoscalereslett.com/content/9/1/90/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://www.nanoscalereslett.com/content/9/1/90/mathml/M15">View MathML</a>. Otherwise, κ(n, h, ∆θ) was also related to the order of GSP. The high order mode had a small skin deep with weak coupling intensity and less blueshift. When h tends to be 0, the grating became too thin to excite the surface mode. This was why the absorption disappeared when h = 0 in Figure  7.

Strong absorption in grating-graphene multilayers

Moreover, the behavior of multilayer structures shown in Figure  2b was also investigated using the modified RCWA and the absorption and reflection spectra were given in Figure  8. When increasing the number of graphene layers, it can be seen that the resonant frequencies do not change but for several lower order modes. Though the reflections were always weak within the resonant range, it is obvious that the more graphene layers included, the stronger the absorption is (almost 90% when it contained 26 graphene layers).

thumbnailFigure 8. The absorption spectrums of grating-graphene periodic multilayer structure. ‘Layers’, number of graphene layers, which is the odd number between 2 and 26. The frequency ranges from 0 to 60 meV (approximately 14.5 THz). The figure inset is the reflections.

The field distributions of Figure  9 also give the same conclusion that the stand waves on each graphene layer were almost oscillated synchronously. The energy was mainly located and absorbed by the graphene layer as we expected.

thumbnailFigure 9. Field distributions. The real part (a) and (b) and magnitude (c) of Ey in multilayer structure of different orders. (a) Excitation at the frequency of 24.6 meV. (b) and (c) Excitation at the frequency of 28.4 meV.

Conclusions

On conclusion, the rigorous coupled wave analysis was modified to compute the excitation of graphene surface plasmon in graphene-containing binary gratings structures. Under the phase matching conditions, the excitation of the graphene surface plasmonics was determined by the distance between graphene layers and duty ratio of gratings, and the mode suppression can be realized by modifying the grating constant and duty ratio. A blueshift of the excitation frequency was obtained for enhanced coupling between GSP of neighbor graphene layers. Increasing the number of graphene layers had almost no effect on the excitation frequency of GSP but would lead to a high absorption with negligible reflection in near-THz range. Finally, the resonant frequency and absorptions can be easily modified by manipulating the structure parameter, including grating constant, duty ratio, and distance between the graphene layers and number of grating, and graphene-containing grating might become potential applications of THz region, such as optical absorption devices, optical nonlinear, optical enhancement, and so on.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HYU, XJ, and XS conceived the idea. HYU, PL, HYA, and XS wrote the codes, calculated the results, and made the conclusions. HYU, XS, and PL contributed to the preparation and revision of the manuscript. All the authors read and approved the final manuscript.

Acknowledgements

This project was supported by the National Basic Research Program of China (no. 2013CB328702) and by the National Natural Science Foundation of China (no. 11374074).

References

  1. Geim AK, Novoselov KS: The rise of graphene.

    Nat Mater 2007, 6:183-191. PubMed Abstract | Publisher Full Text OpenURL

  2. Grigorenko A, Polini M, Novoselov K: Graphene plasmonics.

    Nat Photonics 2012, 6:749-758. Publisher Full Text OpenURL

  3. Bonaccorso F, Sun Z, Hasan T, Ferrari A: Graphene photonics and optoelectronics.

    Nat Photonics 2010, 4:611-622. Publisher Full Text OpenURL

  4. Novoselov K, Geim AK, Morozov S, Jiang D, Grigorieva MKI, Dubonos S, Firsov A: Two-dimensional gas of massless Dirac fermions in graphene.

    Nature 2005, 438:197-200. PubMed Abstract | Publisher Full Text OpenURL

  5. Ju L, Geng B, Horng J, Girit C, Martin M, Hao Z, Bechtel HA, Liang X, Zettl A, Shen YR: Graphene plasmonics for tunable terahertz metamaterials.

    Nat Nanotechnol 2011, 6:630-634. PubMed Abstract | Publisher Full Text OpenURL

  6. Koshino M, Ando T: Magneto-optical properties of multilayer graphene.

    Phys Rev B 2008, 77:115313. OpenURL

  7. Gusynin V, Sharapov S, Carbotte J: Magneto-optical conductivity in graphene.

    J Phys Condens Matter 2007, 19:026222. Publisher Full Text OpenURL

  8. Dressel M: Electrodynamics of Solids: Optical Properties of Electrons in Matter. Cambridge: Cambridge University Press; 2002. OpenURL

  9. Falkovsky L, Pershoguba S: Optical far-infrared properties of a graphene monolayer and multilayer.

    Phys Rev B 2007, 76:153410. OpenURL

  10. Mikhailov SA, Ziegler K: New electromagnetic mode in graphene.

    Phys Rev Lett 2007, 99:016803. PubMed Abstract | Publisher Full Text OpenURL

  11. Stern F: Polarizability of a two-dimensional electron gas.

    Phys Rev Lett 1967, 18:546-548. Publisher Full Text OpenURL

  12. Jablan M, Buljan H, Soljačić M: Plasmonics in graphene at infrared frequencies.

    Phys Rev B 2009, 80:245435. OpenURL

  13. Nikitin AY, Guinea F, Garcia-Vidal FJ, Martin-Moreno L: Surface plasmon enhanced absorption and suppressed transmission in periodic arrays of graphene ribbons.

    Phys Rev B 2012, 85:081405. OpenURL

  14. Nayyeri V, Soleimani M, Ramahi OM: Modeling graphene in the finite-difference time-domain method using a surface boundary condition. In IEEE Transactions on Antennas and Propagation. Piscataway: IEEE; 2013. OpenURL

  15. Peres N, Bludov YV, Ferreira A, Vasilevskiy MI: Exact solution for square-wave grating covered with graphene: surface plasmon-polaritons in the THz range.

    J. Phys. Condens. Matter 25:125303.

    arXiv preprint arXiv:12116358 2012

    OpenURL

  16. Moharam M, Grann EB, Pommet DA, Gaylord T: Formulation for stable and efficient implementation of the rigorous coupled-wave analysis of binary gratings.

    JOSA A 1995, 12:1068-1076. Publisher Full Text OpenURL

  17. Moharam M, Pommet DA, Grann EB, Gaylord T: Stable implementation of the rigorous coupled-wave analysis for surface-relief gratings: enhanced transmittance matrix approach.

    JOSA A 1995, 12:1077-1086. Publisher Full Text OpenURL

  18. Neto AC, Guinea F, Peres N, Novoselov KS, Geim AK: The electronic properties of graphene.

    Rev Mod Phys 2009, 81:109. Publisher Full Text OpenURL

  19. Ziegler K: Robust transport properties in graphene.

    Phys Rev Lett 2006, 97:266802. PubMed Abstract | Publisher Full Text OpenURL

  20. Gusynin V, Sharapov S, Carbotte J: Unusual microwave response of Dirac quasiparticles in graphene.

    Phys Rev Lett 2006, 96:256802. PubMed Abstract | Publisher Full Text OpenURL

  21. Falkovsky L, Varlamov A: Space-time dispersion of graphene conductivity.

    Eur Phys J B 2007, 56:281-284. Publisher Full Text OpenURL

  22. Falkovsky L: Optical properties of graphene.

    Phys. Conf. Ser 2008, 129(1):012004. OpenURL

  23. Hanson GW: Quasi-transverse electromagnetic modes supported by a graphene parallel-plate waveguide.

    J Appl Phys 2008, 104(8):084314-084314-5. Publisher Full Text OpenURL

  24. Mikhailov S, Ziegler K: A new electromagnetic mode in graphene.

    Phys. Rev. Lett. 2007, 99:016803. PubMed Abstract | Publisher Full Text OpenURL

  25. Economou EN: Surface plasmons in thin films.

    Phys Rev 1969, 182:539. Publisher Full Text OpenURL

  26. Petit R: Electromagnetic Theory of Gratings. Heidelberg: Springer Berlin; 1980. OpenURL

  27. Liu H, Lalanne P: Microscopic theory of the extraordinary optical transmission.

    Nature 2008, 452:728-731. PubMed Abstract | Publisher Full Text OpenURL

  28. van Beijnum F, Rétif C, Smiet CB, Liu H, Lalanne P, van Exter MP: Quasi-cylindrical wave contribution in experiments on extraordinary optical transmission.

    Nature 2012, 492:411-414. PubMed Abstract | Publisher Full Text OpenURL

  29. Lalanne P, Hugonin J, Rodier J: Theory of surface plasmon generation at nanoslit apertures.

    Phys Rev Lett 2005, 95:263902. PubMed Abstract | Publisher Full Text OpenURL

  30. Liu N, Langguth L, Weiss T, Kästel J, Fleischhauer M, Pfau T, Giessen H: Plasmonic analogue of electromagnetically induced transparency at the Drude damping limit.

    Nat Mater 2009, 8:758-762. PubMed Abstract | Publisher Full Text OpenURL

  31. Haus HA: Waves and Fields in Optoelectronics. Englewood Cliffs, NJ: Prentice-Hall; 1984. OpenURL