Abstract
Manipulating the circular polarization of light is of great importance in chemistry and biology, as chiral molecules exhibit different physiological properties when exposed to different circularly polarized waves. Here we suggest a graphene/dielectric-stacked structure, which has both the properties of an epsilon-near-zero material and the high Hall conductivity of graphene. The proposed sub-wavelength structure demonstrates efficient manipulation of circular polarization properties of light. In a quite broad frequency range and at a large oblique incidence angle, the present magnetically active structure is transparent for one circularly polarized wave, and opaque for another. Such an effect can be further tuned by changing the magnitude of the applied magnetic field and chemical potential of graphene.
© 2014 Optical Society of America
1. Introduction
The capability to manipulate the circular polarization of light is one of the essential optical applications of molecular biology, medical science and analytical chemistry. Most biomolecules, including the basic block of life, DNA, are chiral, which are sensitive to optical stimuli and thus behave differently when exposed to left-handed circular (LCP) waves and right-handed circular (RCP) waves. The most commonly used method to manipulate circular polarization of light is to use half- or quarter-wave plates with a retardation effect in an anisotropic crystal. However, the small differences in permittivity between the corresponding crystallographic directions require a thick retardation distance, resulting in dozens of wavelength thicknesses. On the other hand, progress with chiral heterostructrues [1–3] provides an unprecedented opportunity to manipulate the circular polarization of light by tailoring the refractive index of the material. Nevertheless, most chiral metastructures are based on Bragg reflection and localized plasmons, which lead to a narrow frequency band. Although the spectrum can be broadened by stacking multiple polarizers or introducing a gradient helical pitch [4], the relative strength of different circularly polarized waves and the working spectrum cannot be tuned without changing the physical structure of the chiral material.
Graphene nanophotonics has attracted considerable interest due to its unique optical properties [5–9] and many related applications, ranging from graphene-based optical modulation [10], transformation optics [11], field-effect transistors [12] to numerous other devices [13–16]. Apart from being the thinnest material, graphene is attractive for its physical flexibility, high electron mobility and the possibility of controlling its carrier concentration via external gate voltages or chemical doping. Another promising application of graphene is in magneto-optics. The ability to create a giant Faraday rotation angle has been indentified experimentally in both single-layer and multilayered continuous graphene [17]. The large Hall conductivity [18], which produces extremely significant Faraday rotation originating from the cyclotron effect in the classical region and the inter-laudau-level transitions in the quantum regime, inevitably affects RCP and LCP waves differently, both in amplitude and phase. Epsilon-near-zero (ENZ) material also exhibits highly unusual and intriguing optical properties. One of these properties is the near-zero phase delays while a wave propagates though such an ENZ anisotropic material, which has been proposed for extraordinary transmission in the bend channel [19]. ENZ materials were also employed for radiation pattern control [20, 21], and perfect absorption [22]. Combined with magneto-optical materials, ENZ materials may provide a mechanism for enhancing the effects of nonreciprocity and time-reversal symmetry breaking [23, 24].
In this paper, by employing the properties of ENZ materials and the high Hall conductivity of graphene, we propose a fast tunable magneto-optical sub-wavelength structure that demonstrates efficient manipulation of circular polarization properties of light. The magnetically active structure is transparent for one circularly polarized wave, and opaque for another. In a fairly broad frequency range and at a large oblique incidence angle, such a polarizing effect can be realized and tuned by changing the magnitude of the applied magnetic field or the chemical potential of graphene.
2. Approach for the design of graphene-based circular polarization beam splitter
To model one-atom-thick graphene in a macroscopic electromagnetic description, the graphene can be treated as an ultrathin film with thickness [11]. Asumming exp(-iωt) time harmonic variation, the fourth Maxwell equation can be rearranged as Denoting the equivalent complex permittivity of Δ thick graphene, we obtain ξg = iσ/(ξ0ωΔ) + 1. As , due to the large value of the first part, the second part could be neglected. The equivalent parallel complex permittivity of the Δ-thick graphene layer is ξg = iσ/(ξ0ωΔ). In the presence of a static magnetic field B applied perpendicular to its surface, graphene is gyrotropic and its anisotropic permittivity can be described by a tensor:
where ξgd = iσxx/(ξ0ωΔ) and ξgod = iσyx/(ξ0ωΔ). σxx and σyx are the longitudinal conductivity (along the E-direction) and Hall conductivity (along the direction perpendicular to E), respectively.Consider a graphene layer which is biased by a static magnetic field perpendicular to the surface, and its electronic property can be characterized by the surface conductivity, which can be modeled by the Kubo formula [25] through a quantum mechanical analysis, as
where νf ≈106 m/s is the Fermi velocity, Mn = are the Laudau energy levels, and is the Laudau energy scale. fd(ξ) = [exp(ξ-µc)/(kBT) + 1]−1 is the Fermi-Dirac distribution, µc is the chemical potential, kB is the Boltzmann constant, T is the environment temperature, ћ = h/(2π), h is the Planck constant, and Γ is the phenomenological scattering rate. The meaning of arrow in Eq. (2) is “replacing Mn with -Mn in the expression of items which have the same priority in calculation”.The total conductivity depends on two parts of electron transition, namely, the intraband transition and interband transition [7]. The former transition involves energy levels in the same band (both in the conduction band or the valence band), and the latter transition involve levels in different bands. As a result, interband transitions occur essentially at frequencies ћω2µc. The spectral range we study in this paper is consistent with the restriction ћω<2µc. Therefore the results of this paper are derived by including the intraband term in (1). Combining this with the fact that the transitions between the levels around μc are the strongest ones, it shows that the conductivities follow the Drude model form [26]
and
where
In this semiclassical expression for the conductivity, τ = 1/(2Γ) is the scattering time and ωc is the cyclotron frequency, which corresponds to the difference between the neighbor Laudau energy around the Fermi level.
A DC mobility of high-quality suspended graphene is expected to be 200,000 cm2V−1s−1 [27, 28], and μ > 60,000 cm2V−1s−1 in graphene on a hexagonal boron nitride (h-BN) substrate was experimentally achieved [29]. In a graphene rinbbon, the relaxation time τ = 20 fs was reported [30]. For a multilayered epitaxial graphene on 4H-SiC substrate, the relaxation time τ = 1 ps was determined through carrier mobility measurements [31]. As graphene sheets with a representative value = 0.2 eV in our suggested structure, it is possible to consider two values for the relaxation time, namely, τ = 2 ps, which is expected to have high mobility 100,000 cm2V−1s−1 by using [32] and result in larger resonance q-factor in grpahene magnetoplasmons, and τ = 0.2 ps, corresponding to lower mobility 10,000 cm2V−1s−1, which can be possibly achieved in some experimental study of graphene.The schematic diagram of the suggested structure is shown in Fig. 1, which consists of periodically alternating graphene and dielectric layers. The dielectric layer has permittivity ξd = 3. A static magnetic field B is applied perpendicularly to the interface of such two materials. The graphene/dielectric stack can be described as an anisotropic metamaterial with effective constitutive parameters, since the thickness of the dielectric layer td is much smaller than the wavelength. There are three effective permittivity parameters, the parallel diagonal and off-diagonal elements ξp and g for the electric fields along the graphene/dielectric interface, and ξt for the electric field perpendicular to the surface.
According to the effective medium theory, we have ξp = f∙ξgd + (1-f)∙ξd, g = f∙i∙ξgod and 1/ξt = f/ξz + (1-f)/ξd, where f = ∆/(∆ + td) is the filling factor of graphene. When ∆ is infinitely small, the three effective permittivity parameters can be reduced to ξp = f∙ξgd + ξd, g = f∙i∙ξgod and ξt = ξd. Since ξg = iσ/(ξ0ω∆) and f = ∆/td, we can easily get ξp = ξd + iσxx/(ξ0ωtd) and g = -σyx/(ξ0ωtd). Given a frequency in the THz range, without properly tuning the magnitude of the applied static magnetic field B or by properly setting the chemical potential µc, which relates to electrostatic biasing and chemical doping, we can achieve near-zero permittivity for ξp. The dependence of the frequency of ENZ with both magnetic field B and chemical potential µc is rendered in Fig. 2(a) and 2(b). Throughout the paper, the two possible values for the dielectric layer thickness are assumed to be 50 nm and 300 nm, which are far below the wavelength corresponding to dozens of THz (the frequency of ENZ is shown in Fig. 2(a), 2(b) and 3(a)), and the effective medium method used here is valid.Optical properties of such a graphene/dielectric stack magnetized along the z-axis are characterized by the permittivity matrix (5). Let us consider the incident electromagnetic wave with the form:
where and correspond to the right-handed and left-handed polarization waves, respectively. θ is the angle of incidence, and the transverse wave vector γ = k0sin(θ).
In the magnetized graphene/dielectric layer:
Note that the electromagnetic field in the magnetic layer consists of four plane waves with different polarization states, which could contain both LH and RH polarization components, especially in oblique incidence. With the Fresnel equation, and in the magnetic media can be calculated as: For the plane wave propagating along the magnetization direction, we have γ = 0, and in the forward direction Eq. (8a) reduces to a well known relation for the corresponding wavevectors . The plus sign is for a forward propagating RCP wave, and the minus sign is for a forward propagating LCP wave. As the convention for the polarization handedness, the backward propagating waves have the wave vectors, . The plus sign is for a LCP wave, and the minus sign is for a RCP wave. The eigenmodes in such a multi-layered structure are two circularly polarized wave , corresponding to wavevector . Assuming a lossless condition, the value of goes to zero, and the wave vectors are real. Therefore the medium is transparent to forward RCP and backward LCP. The other wave vectors are purely imaginary. This means that the medium is opaque for the forward LCP and backward RCP.In the graphene/dielectric alternately stacked material, one can achieve a real part of the diagonal element ξp in matrix (5) near zero, while its imaginary part remains quite small. As magnetically biased graphene possessing giant gyrotropic properties at THz frequencies, we can also achieve a large off-diagonal element g compared to the conventional material. In Fig. 2(c) and 2(d) we show off-diagonal element g as a function of both the applied static magnetic field B and the chemical potential μc. The value of g has a large real part in the frequency range of dozens of THz, which becomes more significant when it is applied by a stronger magnetic field. According to Eq. (4), we can see that the cyclotron frequency has an inverse relationship with the chemical potential, and in Eq. (3) the most significant magneto-optic effect occurs around the cyclotron frequency. If we increase the chemical potential, the frequency of the ENZ will be further away from the cyclotron frequency, resulting in a decrease in the value of g. In Fig. 2(a)-2(d), here are some variations for both the frequency of ENZ and the real value of g, when the relaxation time decreases from 2 ps to 0.2 ps by an order. That is because the introduction of smaller relaxation time will inevitably suppress the strength of gyrotropic resonance.
3. Reflection and transmission for different handed polarization waves by the graphene/dielectric-stacked material
With different values of applied magnetic field and chemical potential, the transmission amplitudes for both circularly polarized waves are shown in Fig. 2(e)-2(h). In Fig. 2(e) and 2(g) (the left panel), as the distance increases, we observe that the RCP forward propagating wave almost preserves its amplitude and has only a little attenuation due to the introduction of an imaginary part in the permittivity. In contrast, the LCP propagating wave is almost completely reduced (most of the energy is reflected) and only a small fraction of it tunnels trough the stack. Such a phenomenon of selected circular polarization transmission will become more significant when we increase the applied static magnetic field and decrease the chemical potential. In the case τ = 0.2 ps (the right panel), both circular polarization waves decay faster than those in the case τ = 2 ps. This is because the imaginary part of the diagonal element ξp in Eq. (5) has a reverse relationship with relaxation time. With the premise at ENZ frequency, the reduced form of Eq. (8a), can be further simplified to and . It can be clearly seen that the loss introduced by reducing the relaxation time will finally bring down the transmission strength for both circular polarization waves. The physics behind such an effect is magnetoplasmons damping, which is related to the relaxation time τ, caused by inter- or intra-band Laudau damping as well as inelastic scattering with phonons and charged impurities. The latter process could be a dominant factor in the situation when the graphene has small a relaxation time.
An incident electromagnetic wave has a linear polarization, which can be decomposed into two opposite circular polarizations, namely, the RCP and LCP propagating waves. Here we assume both circularly polarized waves have unit amplitude and impinge on a graphene/dielectric-stacked layer 15 In order to achieve large value of the off-diagnol element g at the frequency of ENZ, we increase the thickness of dielectric layer to 300 nm. The relationships are shown in Fig. 3(a) and 3(b). With the same strength of magnetic field and chemical potential, the frequency of ENZ has a red shift and g is almost twice as large as the one in Fig. 2. The increase in value of g is due to the approaching of the frequency of ENZ to the cyclotron frequency. The transmission and reflection spectrum of the RCP component are shown in Fig. 3(c). In the lower frequency, most of the wave is reflected. When it comes to the ENZ region, the structure starts to become almost transparent. When we increase the magnitude of the applied magnetic field, the curves have a red shift due to the stronger magneto-optic effect. The spectrum of the LCP component is shown in Fig. 3(d), which is similar to the one for the RCP component: at a lower frequency most of the wave is reflected, and the medium becomes transparent when the frequency increases. However, the positions of the turning points are located at higher frequencies as compared with the RCP components. When we increase the magnitude of the applied magnetic field, the spectrum has a blue shift, which has an opposite moving direction to the RCP component. It is obvious that there is an overlapping region between such two groups of spectra, where the RCP component of the wave can transmit through the structure while the LCP component is reflected. Now we define Pt as the product of the transmission intensity of the RCP component and the quality degree of circular polarization, which is defined as . In the lower frequency, the transmission light has a high quality of circular polarization but a low intensity, while in the higher frequency, the transmission light has a high intensity but low quality of circular polarization, resulting in a near-zero value of Pt for both sides. In the region of ENZ, the suggested structure has both high transmission intensity and high quality of circular polarization. Such a phenomenon is shown in Fig. 3(e). The full width at half maximum (FWHM) of such a peak represents the tolerance of the polarizing effect near the frequency point of ENZ, and has a positive correlation to the magnitude of the applied magnetic field. The relationship is shown in Fig. 3(f). The width of the polarizing band has an almost linear relationship with the magnitude of the magnetic field. However, when we decrease the chemical potential while increasing the magnetic field, the polarizing band will become too close to the cyclotron frequency, where the permittivity of the medium shows some dramatic fluctuation, which results in the disorder at the end of the red line in Fig. 3(f).
Figure 4(a) and 4(b) illustrate the situation of oblique incidence. Here we consider point in Fig. 3(a), with following parameters: the applied magnetic field 3 T, graphene relaxation time 2 ps, the chemical potential 0.2 eV and the corresponding ENZ frequency 9.43 THz. The thickness of the composite layer is 15 μm. The RCP wave and LCP wave incidence are represented in Fig. 4(a) and 4(b), respectively. According to Eq. (8), there are four waves mixing in the magneto-optic material. When we increase the incident angle, it is more likely to excite all four waves in the medium, and thus produce its opposite circular polarization in the transmission spectrum. Both figures show that the graphene/dielectric-stacking polarizer maintains its efficiency at the ENZ region in quite a large incident angle. As the figures show, when the incident angle is less than 45 degrees, the amplitudes of light remain above 80% for both the transmission of RCP incident wave and the reflection of LCP incident wave.
4. Conclusion
We theoretically analyzed the transmission and reflection coefficients of both the RCP and LCP waves at the ENZ frequency range of the graphene/dielectric-stacked material, and these results have been further validated by FDTD calculations, which are much more time-consuming. The composite material shows efficient manipulation of the polarization of the electromagnetic wave in a fairly large range of the incidence angle. Both the transmission ratio of the two circularly polarized waves and the working frequency range can be tuned by changing the magnitude of the applied magnetic field or the chemical potential. The suggested structure with the possibility of easy and fast tuning, together with a strong magneto-optic effect (which is not present in other known materials), may have an impact in a variety of novel devices and applications, beyond circular polarization manipulation and control.
Acknowledgment
The authors are grateful to the partial supports from NSFCs 61178062, 60990322 and 91233208, the National High Technology Research and Development Program (863 Program) of China (No. 2012AA030402), Swedish VRgrant (# 621-2011-4620) and SOARD.
References and links
1. A. Papakostas, A. Potts, D. M. Bagnall, S. L. Prosvirnin, H. J. Coles, and N. I. Zheludev, “Optical manifestations of planar chirality,” Phys. Rev. Lett. 90(10), 107404 (2003). [CrossRef] [PubMed]
2. P. Biagioni, J. S. Huang, L. Duò, M. Finazzi, and B. Hecht, “Cross resonant optical antenna,” Phys. Rev. Lett. 102(25), 256801 (2009). [CrossRef] [PubMed]
3. M. Hentschel, M. Schäferling, T. Weiss, N. Liu, and H. Giessen, “Three-dimensional chiral plasmonic oligomers,” Nano Lett. 12(5), 2542–2547 (2012). [CrossRef] [PubMed]
4. J. K. Gansel, M. Thiel, M. S. Rill, M. Decker, K. Bade, V. Saile, G. von Freymann, S. Linden, and M. Wegener, “Gold helix photonic metamaterial as broadband circular polarizer,” Science 325(5947), 1513–1515 (2009). [CrossRef] [PubMed]
5. A. K. Geim and K. S. Novoselov, “The rise of graphene,” Nat. Mater. 6(3), 183–191 (2007). [CrossRef] [PubMed]
6. A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, “The electronic properties of graphene,” Rev. Mod. Phys. 81(1), 109–162 (2009). [CrossRef]
7. F. H. L. Koppens, D. E. Chang, and F. J. García de Abajo, “Graphene plasmonics: a platform for strong light-matter interactions,” Nano Lett. 11(8), 3370–3377 (2011). [CrossRef] [PubMed]
8. A. N. Grigorenko, M. Polini, and K. S. Novoselov, “Graphene plasmonics,” Nat. Photonics 6(11), 749–758 (2012). [CrossRef]
9. L. Ju, B. Geng, J. Horng, C. Girit, M. Martin, Z. Hao, H. A. Bechtel, X. G. Liang, A. Zettl, Y. R. Shen, and F. Wang, “Graphene plasmonics for tunable terahertz metamaterials,” Nat. Nanotechnol. 6(10), 630–634 (2011). [CrossRef] [PubMed]
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). [CrossRef] [PubMed]
11. A. Vakil and N. Engheta, “Transformation optics using graphene,” Science 332(6035), 1291–1294 (2011). [CrossRef] [PubMed]
12. F. N. Xia, D. B. Farmer, Y. M. Lin, and P. Avouris, “Graphene field-effect transistors with high on/off current ratio and large transport band gap at room temperature,” Nano Lett. 10(2), 715–718 (2010). [CrossRef] [PubMed]
13. R. Alaee, M. Farhat, C. Rockstuhl, and F. Lederer, “A perfect absorber made of a graphene micro-ribbon metamaterial,” Opt. Express 20(27), 28017–28024 (2012). [CrossRef] [PubMed]
14. Q. Bao, H. Zhang, B. Wang, Z. Ni, C. H. Y. X. Lim, Y. Wang, D. Y. Tang, and K. P. Loh, “Broadband graphene polarizer,” Nat. Photonics 5(7), 411–415 (2011). [CrossRef]
15. S. He, X. Z. Zhang, and Y. R. He, “Graphene nano-ribbon waveguides of record-small mode area and ultra-high effective refractive indices for future VLSI,” Opt. Express 21(25), 30664–30673 (2013). [CrossRef] [PubMed]
16. S. He and T. Chen, “Broadband THz absorbers with graphene-based anisotropic metamaterial films,” IEEE Trans. Terahertz Sci. Technol. 3(6), 757–763 (2013). [CrossRef]
17. I. Crassee, J. Levallois, A. L. Walter, M. Ostler, A. Bostwick, E. Rotenberg, T. Seyller, D. van der Marel, and A. B. Kuzmenko, “Giant faraday rotation in single- and multilayer graphene,” Nat. Phys. 7(1), 48–51 (2011). [CrossRef]
18. D. L. Sounas and C. Caloz, “Electromagnetic nonreciprocity and gyrotropy of graphene,” Appl. Phys. Lett. 98(2), 021911 (2011). [CrossRef]
19. M. G. Silveirinha and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using epsilon-near-zero materials,” Phys. Rev. Lett. 97(15), 157403 (2006). [CrossRef] [PubMed]
20. S. Enoch, G. Tayeb, P. Sabouroux, N. Guérin, and P. A. Vincent, “A metamaterial for directive emission,” Phys. Rev. Lett. 89(21), 213902 (2002). [CrossRef] [PubMed]
21. A. Alu, M. G. Silveirinha, A. Salandrino, and N. Engheta, “Epsilon-near-zero metamaterials and electromagnetic sources: Tailoring the radiation phase pattern,” Phys. Rev. B 75(15), 155410 (2007). [CrossRef]
22. S. M. Zhong, Y. Ma, and S. He, “Perfect absorption in ultrathin anisotropic ε-near-zero metamaterials,” Appl. Phys. Lett. 105(2), 023504 (2014). [CrossRef]
23. A. R. Davoyan, A. M. Mahmoud, and N. Engheta, “Optical isolation with epsilon-near-zero metamaterials,” Opt. Express 21(3), 3279–3286 (2013). [CrossRef] [PubMed]
24. A. R. Davoyan and N. Engheta, “Theory of wave propagation in magnetized near-zero-epsilon metamaterials: evidence for one-way photonic states and magnetically switched transparency and opacity,” Phys. Rev. Lett. 111(25), 257401 (2013). [CrossRef] [PubMed]
25. R. Kubo, “Statistical-mechanical theory of irreversible processes. i.general theory and simple applications to magnetic and conduction problems,” J. Phys. Soc. Jpn. 12(6), 570–586 (1957). [CrossRef]
26. V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, “On the universal ac optical background in graphene,” New J. Phys. 11(9), 095013 (2009). [CrossRef]
27. K. I. Bolotin, K. J. Sikes, Z. Jiang, M. Klima, G. Fudenberg, J. Hone, P. Kim, and H. L. Stormer, “Ultrahigh electron mobility in suspended graphene,” Solid State Commun. 146(9-10), 351–355 (2008). [CrossRef]
28. S. V. Morozov, K. S. Novoselov, M. I. Katsnelson, F. Schedin, D. C. Elias, J. A. Jaszczak, and A. K. Geim, “Giant intrinsic carrier mobilities in graphene and its bilayer,” Phys. Rev. Lett. 100(1), 016602 (2008). [CrossRef] [PubMed]
29. C. R. Dean, A. F. Young, I. Meric, C. Lee, L. Wang, S. Sorgenfrei, K. Watanabe, T. Taniguchi, P. Kim, K. L. Shepard, and J. Hone, “Boron nitride substrates for high-quality graphene electronics,” Nat. Nanotechnol. 5(10), 722–726 (2010). [CrossRef] [PubMed]
30. H. G. Yan, T. Low, W. J. Zhu, Y. Q. Wu, M. Freitag, X. S. Li, F. Guinea, P. Avouris, and F. N. Xia, “Damping pathways of mid-infrared plasmons in graphene nanostructures,” Nat. Photonics 7(5), 394–399 (2013). [CrossRef]
31. M. Orlita, C. Faugeras, P. Plochocka, P. Neugebauer, G. Martinez, D. K. Maude, A. L. Barra, M. Sprinkle, C. Berger, W. A. de Heer, and M. Potemski, “Approaching the dirac point in high-mobility multilayer epitaxial graphene,” Phys. Rev. Lett. 101(26), 267601 (2008). [CrossRef] [PubMed]
32. J. Christensen, A. Manjavacas, S. Thongrattanasiri, F. H. Koppens, and F. J. G. de Abajo, “Graphene plasmon waveguiding and hybridization in individual and paired nanoribbons,” ACS Nano 6(1), 431–440 (2012). [CrossRef] [PubMed]