Jiaxin Zhouab,
Yueke Wang*ab,
Mengjia Luab,
Jian Dingab and
Lei Zhouab
aOptical Information Science and Technology Department, Jiangnan University, Wuxi, Jiangsu 214122, China. E-mail: ykwang@jiangnan.edu.cn
bOptoelectronic Engineering and Technology Research Center, Jiangnan University, Wuxi, Jiangsu 214122, China
First published on 21st October 2019
In this letter, we propose a structure based on double-layer graphene-based planar chiral metasurface with a J-shaped pattern to generate asymmetric transmission for circularly polarized waves in the mid-infrared region. Asymmetric transmission of the double-layer structure can reach to 16.64%, which is much larger than that of the monolayer. The mechanism of asymmetric transmission is attributed to enantiomerically sensitive graphene's surface plasmons. Besides, asymmetric transmission can be dynamically tuned by changing the Fermi energy and is affected by intrinsic relaxation time. All simulations are conducted by the finite element method. Our findings provide a feasibility of realizing photonic devices in tunable polarization-dependent operation, such as asymmetric wave splitters and circulators.
As a kind of promising electrically tunable plasmonic material, periodical graphene metasurfaces can excite graphene surface plasmons in the mid-infrared and terahertz (THZ) region.20–22 And for graphene metasurface, novel optical properties are dynamically tunable by changing the Fermi energy through voltage control or the changing intrinsic relaxation time through chemical doping.23–26 Compared with metallic surface plasmons, graphene surface plasmons have the advantages of low loss, and high localization.27 Thus, graphene will revolutionize metamaterials and metadevices, and promote the development of nano-optics.28 There are plenty of researches on metallic chiral metamaterials,29 but they are generally not tunable, which can be easily realized in graphene metasurface. Recently, people propose kinds of monolayer graphene chiral metasurface,30,31 and study the asymmetric transmission for circularly polarized waves. The researches of graphene planar chiral structure will deepen the understanding of the interaction between light and matter in planar chiral structure, which can be applied to tunable polarization sensitive devices, and circular polarizers. But the asymmetric transmission of monolayer graphene chiral metasurface is only up to 5%.
In this letter, we propose a double-layer graphene-based planar chiral metasurface structure to achieve giant enhancement of asymmetric transmission for circularly polarized waves incidence in mid-infrared region. In the monolayer metasurface which is composed of J-shaped hollow graphene patterns, the asymmetric transmission ΔT can reach a small value of 2.05%. When another J-shaped pattern, which is designed to be rotated by 90 degrees after being mirror symmetric operation with respect to the previous layer, is added, the asymmetric transmission ΔT can achieve 16.64%. The mechanism of asymmetric transmission can be explained with the enantiomerically sensitive graphene surface plasmons, which induce the charge-field excitations. Also, the peaks of asymmetric transmission blue shift with increasing Fermi level, and the values of peaks increase with increasing electron scattering time. Finite Element Method (FEM) simulations are conducted to verify our findings.
According to the Kubo formula, the conductivity of graphene consists of two parts: intraband and interband elements:32
(1) |
(2) |
(3) |
(4) |
The pattern of the J-shaped graphene metasurface has a character of chirality, which can reveal asymmetric transmission when the incident circularly polarized waves have different handedness or propagation directions. Fig. 1(c) and (d) demonstrate the different incident directions of LCP waves, which are coming perpendicularly from the front and the back of monolayer graphene-based planar chiral metasurface. In the complex circular transmission matrices Ei = tijE0j, the indices i, j represents right-handed circularly polarized (RCP) or left-handed circularly polarized (LCP) components. The transmission matrix elements are defined as Tij = |tij|2. The circular transmission T can be applied to any dispersive optical system under the condition of coherent light incidence. The transmission matrix elements are the same for circularly polarized waves of opposite handedness propagating in opposite directions. It means that the LCP(RCP) waves incident forward equal to the RCP(LCP) waves incident backward. Thus, we can define the transmittances of the LCP waves with opposite propagation directions as
Tf = |tLL|2 + |tRL|2 | (5) |
Tb = |tRR|2 + |tLR|2 | (6) |
ΔT = ΔLCPcirc = Tf − Tb = −ΔRCPcirc | (7) |
Thus, the asymmetric transmission is only decided by circular polarization conversion transmission.
Fig. 2(a) and (b) show the calculated results of the four transmission matrix elements when the Fermi energy Ef is fixed to be 0.80 eV, for the monolayer J-shaped graphene metasurface. Black and red lines represent LCP incidence along the forward and backward directions, respectively. It is found that the direct transmissions (TLL and TRR) always keeps the same with each other with the wavelength changing, and reaches minimum when the incident wavelength is 14.83 μm, shown as Fig. 2(a). It proves that graphene-based planar chiral metasurface is lack of circular dichroism. In Fig. 2(b), TRL is always a little larger than TLR with the wavelength ranging from 13 to 17 μm, and the maxima of TRL and TLR are 15.07% and 13.02% at 14.83 μm, respectively. Thus, the asymmetric transmission for LCP waves reaches the maximum value of 2.05% at the wavelength of 14.83 μm, as shown in Fig. 2(c). It is believed that the graphene surface plasmon is resonantly excited along the edges of the hollow J-shaped pattern at 14.83 μm, and in charge of the asymmetric transmission. For RCP incident waves, the asymmetric transmission is quite the contrary to that of LCP incidence based on eqn (7), shown as Fig. 2(d).
Although the asymmetric transmission has been obtained, the ΔT is small, compared with that of metallic planar chiral metamaterial structure.34,35 To realize larger asymmetric transmission, a double-layer structure is proposed based on the monolayer graphene-based J-shaped metasurface structure. Fig. 1(e) shows the structure composed of two graphene-based planar metasurface layers with J-shaped patterns. The pattern of the bottom layer is shown as Fig. 1(f), which is the same with the monolayer structure of Fig. 1(b); while the upper layer shown as Fig. 1(g) is the pattern of the bottom J-shaped pattern being clockwise rotated by 90 degrees after mirror symmetric transformation about y-axis. The medium between the two graphene layers is quartz, and the thickness of the quartz is 125 nm. The LCP waves are still coming perpendicularly from the front and the back of graphene-based planar chiral metasurfaces, respectively.
Fig. 3(a) and (b) show the calculated results of the four transmission matrix elements, when the Fermi energy Ef is fixed to be 0.80 eV for the double-layer structure shown as Fig. 1(e). Fig. 3(a) shows that TLL and TRR are still identical, which proves that asymmetric transmission is independent of direct transmission. In Fig. 3(b), the transmission matrix element TRL is slightly higher than that of the monolayer chiral metasurface shown in Fig. 2(b), while TLR goes down sharply to almost zero at the wavelength of 14.90 μm. Compared with the monolayer J-shaped metasurface, the resonant wavelength of graphene surface plasmon redshifts a little. Due to the existence of the near to zero value of TLR, the asymmetric transmission ΔT can be increased a lot. As shown in the Fig. 3(c), the maximum of ΔT can reach 16.64% at the wavelength of 14.90 μm, which is much larger than 2.05% of the monolayer graphene-based planar chiral metasurface. Similarly, the asymmetric transmission peak of RCP incident waves reaches −16.64%, just the opposite of LCP. The giant enhancement of circular asymmetric transmission suggests that the structure of double-layer J-shaped graphene metasurface can greatly improve chirality.
In order to understand the physical mechanisms behind the asymmetric transmission, the induced electric field in the hollow area of the bottom layer of graphene metasurface is studied in Fig. 4 and the induced electric vector field is indicated by red arrows. At the resonance wavelength, a standing wave of current oscillations is excited by electromagnetic radiation. This coupled oscillation of the electron density and electromagnetic field is a localized plasmon. Plasmonic excitation of chiral structure is different from that of ordinary metamaterials since it is enantiomerically sensitive.36 The excitation level is decided by the direction and polarization state of the incident waves. Here, graphene localized surface plasmon is excited along the J-shaped edge because the electric field is enhanced around the J-shaped edge shown in Fig. 4. Here, the asymmetric transmission is attributed to the enantiomerically sensitive graphene surface plasmons. When the LCP (or RCP) at the resonance wavelength is incident on the graphene metasurface from the forward (or the backward) direction, the enantiomerically graphene surface plasmon is excited, enhances the local electric field, and excites the induced electric field in the substrate, which determines the radiation (direct transmission, and conversion transmission). The striking difference between plasmonic excitations of opposite handedness leads to the different induced electric field. With greater chirality, the induced electric field is more different and the difference of the circular polarization conversion between the opposite handedness will be stronger. The components Ey of the induced electric field for LCP(RCP) incidences can reflect the excitation level of localized plasmon, and explain the difference between TLR and TRL. Here, the yellow and green arrows was painted to demonstrate the components Ey along y-directions as shown in Fig. 4. Fig. 4(a) and (b) show the induced electric field distribution under the monolayer graphene metasurface in the substrate. The green arrows for the RCP incidence are stronger than that of the LCP incidence for the monolayer metasurface structure, which indicates the suppression of circular polarization conversion is stronger. Thus, TLR is smaller than TRL, and an asymmetric transmission (ΔT = TRL − TLR) of 2.05% is achieved for the monolayer metasurface structure. Fig. 4(c) and (d) show the induced electric field distribution under the double-layer graphene metasurface in the substrate when the incident wavelength is 14.90 μm. The induced electric field components Ey for the LCP incidence indicated by yellow arrows in Fig. 4(c) are stronger than those in Fig. 4(a), which indicates that the suppression of circular polarization conversion is smaller. Thus, TRL of double-layer increases to 17.13%, compared with 15.07% of the monolayer. As shown in the Fig. 4(d), more electric fields indicated by the green arrows and less electric fields indicated by the yellow arrows than those of Fig. 4(b) means that the electric field is weakly coupling to free space, scattering is low and suppression of circular polarization conversion is much stronger than that of monolayer metasurface.17,36 Thus, TLR is near to zero, and the asymmetric transmission of the double-layer metasurface can be greatly increased and ΔT can reach 16.64%.
To study the dynamical tunability of asymmetric transmission, we changed the Fermi energy and intrinsic relaxation time of graphene layers. When increasing the Fermi energy from 0.76 to 1.00 eV, the peak of asymmetric transmission has a blue shift changing from 15.30 to 13.30 μm, as shown in Fig. 5(a). It suggests that graphene-based planar chiral metasurface layer structure has a feasibility of realizing polarization sensitive devices in a wide wavelength range. This blue shift behaviour can be interpreted by the following: the excitation of graphene surface plasmonic wave vector is satisfied with the equation Kspp = ħω2/2α0Efc ∝ 1/Lg37 where α0 = e2/ħc is the fine structure constant. The resonant wavelength λ can be written as:
(8) |
Fig. 5 (a) The relation between the asymmetric transmission ΔT and wavelength under different Fermi energy. (b) The relation between ΔT and wavelength under different scattering time τ, when Ef is 0.8 eV. Other parameters are the same with that of Fig. 1(e). |
Other double-layer structures of J-shaped pattern have been also studied, and circular polarization conversion are calculated, when Ef is fixed to be 0.80 eV. Fig. 6(a), (d) and (g) show the upper patterns of the three typical kinds of double-layer structures, which are the same with the bottom layer pattern, the bottom layer pattern being rotated by 90 degrees and the bottom layer pattern being mirror rotated about y-axis, respectively. Fig. 6(b), (e) and (h) show the circular polarization conversion of the three kinds of two-layer structures, it is found that the difference between TLR and TRL varies with wavelength. For the structures of Fig. 6(a) and (g), the asymmetric transmission can reach maximum 5.19% and 7%, when λ locates at 15.1 and 14.8 μm, shown as Fig. 6(c) and (i). And for the structure of 90 degrees rotation, the asymmetric transmission almost disappears. Thus, the largest asymmetric transmission can be obtained for the double-layer structure of Fig. 1(e).
This journal is © The Royal Society of Chemistry 2019 |