Grounded Graphene-Based Nano-Waveguide with Chiral Cover and Substrate: New Theoretical Investigation
11 Grounded Graphene-Based Nano-Waveguide with Chiral Cover and Substrate: New Theoretical Investigation
Mohammad Bagher Heydari , Mohammad Hashem Vadjed Samiei School of Electrical Engineering, Iran University of Science and Technology (IUST), Tehran, Iran * Corresponding author: [email protected] ; [email protected]
Abstract:
In this paper, we propose a new analytical model for grounded chiral slab waveguides covered with graphene sheets. The general waveguide is constructed of a graphene sheet sandwiched between two different chiral layers (as substrate and cover layers), each one has the permittivity, the permeability and the chirality of ɛ 𝑐 , 𝜇 𝑐 , 𝛾 𝑐 , respectively. The substrate is supposed to be a perfect electromagnetic conductor (PEMC), which is able to easily convert to a perfect electric conductor (PEC) or a perfect magnetic conductor (PMC). A novel matrix representation is obtained for the general structure to find its dispersion relation and other propagating parameters. To show the richness of the proposed waveguide, two new structures have been introduced and investigated. It has been shown that the modal properties of these exemplary structures are tunable by changing the chemical potential of the graphene and the chirality. The proposed general structure and its analytical model can be utilized for designing tunable plasmonic devices in THz frequencies. Key-words:
Analytical model, graphene sheet, chiral medium, plasmon, PEMC, effective index Introduction
Nowadays, graphene plasmonics is a new emerging science, which studies the excitation of surface plasmon polaritons (SPPs) on the graphene [1-3]. Due to some fascinating features of these SPPs, such as high mode confinement, a large variety of graphene-based devices have been reported at THz frequencies, such as waveguides [4-14], isolator [15], circulator [16, 17], coupler [18], resonator [19], antennas [20-22], filter [23], Radar Cross-Section (RCS) reduction-based devices [24-26], and graphene-based medical components [27-33]. It should be noted that noble metals support SPPs at the near-infrared and visible frequencies [19, 34-43]. Chiral materials, one of the interesting materials in chemistry science with many remarkable features such as optical activity and circular dichroism, have opened a new area for researchers to design and report innovative devices such as cloak [44], polarization rotator [45] and microstrip antenna [46]. The hybridization of graphene and chiral medium in a THz device is very fascinating because the plasmonic features of the proposed device can be controlled and adjusted simultaneously by changing the chirality and the chemical potential of the graphene [5, 47, 48]. For instance, the modal properties of a graphene-chiral interface have been discussed and studied in [47]. The authors in [5] have been investigated the effect of various factors on the plasmonic modes of a graphene-chiral-graphene slab waveguide. In this paper, we aim to analyze and consider the plasmonic features of grounded chiral slab waveguides covered with graphene sheets. To the authors’ knowledge, the general proposed structure and its analytical model has not been published in any research article. Our presented waveguide is a generalization of all types of hybrid graphene-chiral grounded slab waveguides. In our general structure, a graphene sheet has been sandwiched between two chiral layers, as cover and substrate layers. The ground is a PEMC boundary condition in the general form, which can be easily transformed into PEC or PMC boundaries. In this generalized structure, one side of the proposed waveguide is always grounded by a PEMC layer and the other side is open . It should be noticed that our presented study considers only chiral materials with real-valued chirality and the study of generalized bi-anisotropic materials with the complex-valued tensors (i.e. epsilon, mu, and gamma) is complicated and is outside the scope of this article.
The general waveguide supports hybrid plasmonic modes, which are split into two kinds of modes, called “Higher modes” and “Lower modes” in this paper. Indeed, chiral materials are responsible for generating these modes.
The paper is organized as follows. Section 2 introduces the general waveguide and proposes its analytical model. New matrix representation is obtained in this section that its determinant achieves the propagation constant of the general structure. Section 3 explains how the dispersion relation can be computed for a specific case of the general structure and also review some numerical methods for solving nonlinear equations. In section 4, two new graphene-based waveguides with chiral substrates are studied to show the richness of our proposed structure. The first waveguide is composed of a graphene sheet located on the grounded chiral slab. The ground plane in this structure is supposed to be PEC. We have achieved high mode confinement ( 𝑛 𝑒𝑓𝑓 = 28 ) at the frequency of 20 THz for it. The second case is a graphene sheet placed on the chiral slab backed by a PEMC layer. The hybridization of graphene with chiral materials allows the designer to tune and control SPPs via the chirality and the chemical potential. Finally, section 5 concludes the article. The Proposed Structure and its Analytical Model
In Fig. 1, the configuration of the proposed waveguide has been illustrated. In this structure, a graphene sheet has been sandwiched between two chiral layers (as cover and substrate layers) and it has isotropic conductivity, which is described by Kubo’s formula [49]:
22 2 g B g gBg Bg
K T je K TjeT Ln Ln eK T (1) Fig. 1.
The configuration of the proposed structure: (a) the 3D schematic, (b) the cross-section. The cover and substrate are assumed to be different chiral layers. In the above relation, ℎ is the reduced Planck’s constant, 𝐾 𝐵 is Boltzmann’s constant, ω is radian frequency, 𝑒 is the electron charge, 𝛤 is the phenomenological electron scattering rate ( 𝛤 = 1 𝜏⁄ , where 𝜏 is the relaxation time), 𝑇 is the temperature, and 𝜇 𝑔 is the chemical potential which can be altered by chemical doping or electrostatic bias [49]. Moreover, each chiral layer satisfies the following relations [50]: 𝑃𝐸𝑀𝐶 𝑡 𝑧 𝑦 𝑥 (𝒂) ɛ 𝑐,1 , 𝜇 𝑐,1 , 𝛾 𝑐,1 𝑃𝐸𝑀𝐶 𝑧 𝑥 𝑦 𝜎 (𝒃) 𝑡 ɛ 𝑐,1 , 𝜇 𝑐,1 , 𝛾 𝑐,1 ɛ 𝑐,2 , 𝜇 𝑐,2 , 𝛾 𝑐,2 ɛ 𝑐,2 , 𝜇 𝑐,2 , 𝛾 𝑐,2 𝜎 , , N c N N c N N j D E B (2) , ,
NN c N N c N j BH E (3)
Where ɛ 𝑐,𝑁 , 𝜇 𝑐,𝑁 , 𝛾 𝑐,𝑁 are the permittivity, the permeability, and the chirality of that layer ( N shows the index of each layer, N=1, 2 ), respectively. Furthermore, ɛ , 𝜇 are the permittivity and permeability of free space, respectively. In our proposed structure, the ground plane is a PEMC layer with the following boundary condition [51]: M n H E (4) In relation (4), M is a PEMC admittance that can be transformed into PEC ( 𝑀 → ±∞ ) and PMC (
𝑀 = 0 ) boundaries. To obtain the analytical model for our waveguide, we write the electromagnetic fields inside the chiral medium as the superposition of left-handed circularly polarized (LCP) and right-handed circularly polarized (RCP) waves [50]. Here, we assume that the electromagnetic waves propagate in the z-direction ( 𝑒 𝑗𝛽𝑧−𝑖𝜔𝑡 ) and they have no variations in the y-direction ( 𝜕𝜕𝑦 = 0 ). Therefore, the Helmholtz’s equation for the N -th layer ( N=1, 2 ) is [52]: R N R N R N d qdx (5) L N L N L N d qdx (6)
Where
R N R N k q (7) L N L N k q (8) In (7)-(8), 𝑘 𝑅,𝑁 , 𝑘
𝐿,𝑁 are the propagation constants of RCP and LCP waves [50]: , , ,
R N c N c N k n (9) , , , L N c N c N k n (10) Where , , , c N c N c N n (11) is the refractive index of the N -th chiral layer [50]. It should be mentioned that the impedance of the N -th chiral layer has been defined as follows [50]: ,, , c Nc N c N (12) Hence, the electromagnetic fields inside the chiral medium can be expressed as [53]: , , , , R N z N c N z N
E j H (13) , , , , L N z N c N z N
E j H (14)
Which results in [53]: , L, , R N N z N E (15) , L, , , R N N z N c N
H j (16) The transverse components of electric and magnetic fields are obtained by utilizing Maxwell’s equations: ,, x N E NNy N E QE x (17) ,, , x N H NNy N c N H QH x (18) Where the following definitions have been used in the relations (17)-(18):
EN R L j jQ k k (19)
HN R L
Q j k j k (20) ,2,,2,
R NR NN L NL N qq (21) Let us suppose that the propagation constants of RCP and LCP waves are: ,1 ,2 RR R x tj x t qq q (22) ,1 ,2 L L L x tj x t qq q (23)
Now, the Helmholtz’s relations are written as: ,1 ,1,2 cos sin 0exp
R R R RR R R
A x B x x tC x x t q qq (24) ,1 ,1,2 cos sin 0exp
L L L LL L L
A x B x x tC x x t q qq (25)
Where 𝐴 𝑅 , 𝐵 𝑅 , 𝐶 𝑅 , 𝐴 𝐿 , 𝐵 𝐿 , 𝐶 𝐿 are the unknown coefficients and should be determined by applying boundary conditions. By using relation (4), boundary conditions at 𝑥 = 0 are obtained as [51]: ,1 ,1 ,1 ,1 y y z z H M E H M E (26) At 𝑥 = 𝑡 , the boundary conditions are written as [54]: ,1 ,2 ,1 ,2 , z z z y y y E E E E E E (27) ,2 ,1 ,2 ,1 , z z y y y z H H E H H E (28)
By applying boundary conditions expressed in (26)-(28), the final matrix representation is achieved:
RLRLRL
AABS BCC (29)
In (29), the matrix 𝑆̿ is: ,1 ,1,1 ,1 ,1 ,1,1 ,1,1 ,1 ,1 ,1 ,2 ,2,1 ,1 ,1 ,1 ,1 ,1 ,1,1 ,1 ,1 R LR c L cc cR L R L R LR R L L R R L LR L R S k kj jM Mq qM Mj jq t q t q t q t q t q tk q t k q t k q t k qq q q ,1 ,2 ,2 ,2 ,2,1 ,2 ,2,2 ,2,1 ,1 ,1 ,1 ,2 ,2 ,2 ,2,1 ,1 ,1 ,1 ,2 ,2,1 ,1 ,1 ,1,1 ,1 exp exp1 1cos cos sin sin exp expsin sin R R L LL R LR LR L R L c R c Lc c c c R LR R L Lc R t k q t k q tq q qk kq t q t q t q t j q j qj j j j q t q tk q t k q tj jq ,1 ,1 ,1 ,1 ,2 ,2 ,2 ,2,1 ,1 ,1 ,1 ,1 ,1 ,2 ,2 cos cos exp exp
R R L L c R c Lc L c R c L R L j jk q t k q t q qj jq q q q t q t (30)
Now, by setting det(𝑆̿) = 0 , the dispersion relation (or the propagation constant, 𝛽 ) is found. Then, achieving plasmonic parameters such as the effective index ( 𝑛 𝑒𝑓𝑓 = 𝑅𝑒[𝛽] 𝑘 ⁄ ) and the propagation loss ( 𝐿 𝑃𝑟𝑜𝑝 = 1 (2𝐼𝑚[𝛽])⁄ ) is straightforward. Methods
In this section, we briefly study the calculation of the dispersion relations and also review numerical techniques for solving them. First, we will compute the dispersion relation for a specific case and second, we will consider the numerical solution of the relation (30) in MATLAB software.
Fig. 2.
The graphene-based waveguide with a chiral substrate backed by a PEC layer: (a)
The
3D schematic, (b)
The cross-section of the structure in the z-x plane. 𝑡 𝑃𝐸𝐶 ɛ , 𝜇 𝜎 𝑃𝐸𝐶 𝑡 𝑧 𝑦 𝑥 (𝒂) 𝑧 ɛ , 𝜇 (𝒃) ɛ 𝑐 , 𝜇 𝑐 , 𝛾 ɛ 𝑐 , 𝜇 𝑐 , 𝛾 𝑥 𝜎 Consider the waveguide depicted in Fig. 2, where a graphene sheet has been placed on the grounded chiral layer. The ground plane is supposed to be PEC. Hence, we will utilize the mathematical procedures of the previous section for
𝑀 → ∞ . Consider the propagation constants of RCP and LCP waves for this structure: RR x tj x t qq q (31) L L x tj x t qq q (32) Now, the Helmholtz’s relations are written as: cos sin 0exp R R R RR R
A x B x x tC x x t q qq (33) cos sin 0exp L L L LL L
A x B x x tC x x t q qq (34) Where 𝐴 𝑅 , 𝐵 𝑅 , 𝐶 𝑅 , 𝐴 𝐿 , 𝐵 𝐿 , 𝐶 𝐿 are the unknown coefficients and should be determined by applying boundary conditions. By using the relations (26)-(28) for 𝑀 → ∞ , the closed-form of the dispersion relation for the structure is derived: cos cos .sin cos .cos sinexp 2 1sin .sin cos .sin sin
R R R R L RR L L R L
P q t P q t q t P q t q t P q tjq t P q t q t P q t q t P q t (35) Where c R R L L R L R L R L L c R R L L
P k k q k q k k k q q q k q k q q k q k q (36) c R R L L R L R L R L L c R R L LR L R R L
P k k q k q j k k k q q q k q k q q k q k qj k q q j k q k q (37) c R R L L R L R L L R R c R R L L
P k k q k q k k k q q q k q k q q k q k q (38) c R R L L R L R L R L L c R R L L P k k q k q k k k q q q k q k q q k q k q (39) c L R c R L R L R L L R R c L RR L L c R L P k k q k k q q q k k k k q q k k q q k q qk k q q k q q (40) .2 R R L L L L R L R R L R LR L R c R R L L
P k k q k q j q k q j k k q q j k k k q q qk k q q k q q k q (41) c R R L L R L R L L R R c R R L L P k k q k q k k k q q q k q k q q k q k q (42) As seen in the above relations, the closed-form relation is so complicated and it is better to use numerical methods for solving the general matrix relation of (30). Some numerical software such as MATLAB can calculate the determinant of a matrix with an arbitrary size, by using “Cofactor expansion” [55]. Here, to obtain the determinant of matrix (30), one should define the frequency and the propagation constant as symbolic (by using command “syms” in MATLAB) and then utilize the command “det” to find the determinant of general matrix (defined in (30)) as a function of frequency and the propagation constant. Let us assume that the outcome of MATLAB software for the determinant of the relation (30) is: det , S f (43) Now, to solve the nonlinear equation of 𝑑𝑒𝑡(𝑆̿) = 0 , numerical methods in mathematics such as Newton-Raphson [56] can be used, which exist in MATLAB software. Hence, one can utilize MATLAB software to solve the nonlinear equation of 𝑑𝑒𝑡(𝑆̿) = 0. Special Cases of the Proposed Structure: Results and Discussions
This section introduces and studies two new graphene-based waveguides incorporating chiral substrates. Without loss of generality, the cover layer in both structures is supposed to be air. The substrate is a chiral layer in both waveguides and it has different parameters in them. Also, the ground plane is PEC and PEMC boundary conditions in the first and second structures, respectively. Both waveguides support hybrid plasmonic modes, which are tunable via the chirality and the chemical potential of the graphene sheet. In chiral-based waveguides, hybrid modes are split into two “bifurcated modes”, which is one of the interesting properties of these structures [50, 57, 58]. In this paper, we call them “Higher modes” and “Lower modes” that are relevant to high and low plasmon resonance frequencies, respectively. Bifurcated modes are hybrid modes with the same cut-off frequency and different propagation constants. In our studied structures, chiral materials are responsible for generating these modes. In the literature, mode bifurcation has been reported in some structures such as parallel plate waveguides [50], circular waveguides [58], rectangular waveguides [57] and diagonal anisotropic chiral waveguides [59]. The nature of bifurcating modes depends on the studied chiral structure. For instance, in circular waveguides, modes with 𝑒 −𝑖𝑚𝜑 and 𝑒 𝑖𝑚𝜑 ( 𝑚 ≠ 0 ) azimuth variations have different phase constants (above the cut-off frequency) and the same cut-off frequency, which form bifurcated modes [58]. Our proposed, general structure and its specific cases (the first and the second structures) are a kind of parallel plate waveguides. In these structures, odd and the next even mode have the same cut-off frequency and thus they are a pair of bifurcated modes [50, 58]. The mathematical explanations of emerging bifurcated modes in parallel plate chiro-waveguides have been reported in [50]. In all the following analytical results, the temperature is 𝑇 = 300 𝐾 and the relaxation time of the graphene is 𝜏 =0.4 𝑝𝑠 . The chemical potential of the graphene is assumed to be 𝜇 𝑔 = 0.3 𝑒𝑉 unless otherwise stated. It is worthwhile to be mentioned that the relaxation time of the graphene depends on the phenomenological scattering rate and the quality of graphene [1, 60-62]. In what follows, we will investigate the modal properties of these waveguides. The configuration of the first structure was introduced in the previous section (see Fig. 2) and a closed-form equation was obtained for its dispersion relation. Here, we will consider and investigate the numerical results. In our simulations, the chiral substrate has a thickness of 𝑡 = 600 𝑛𝑚 and the effective index of 𝑛 𝑐 = 1.3 . To indicate the tunability of plasmonic waves via chirality, two values of chirality (i.e. 𝛾 = 0.0025,0.0035 𝛺 −1 ) have been considered here . These values for the chirality and the index of the chiral medium depend upon the nature and structure of the organic, inorganic and biochemical molecules, which are reported in [63-69]. Fig. 3 represents the analytical results of the effective index and the propagation length. As mentioned before, two plasmonic modes (Higher and Lower modes) are appeared here due to the existence of the chiral substrate. One can observe from Fig. 3(a) that the effective index increases as the frequency increases. While the propagation length decreases with the frequency increment, as seen in Fig. 3(b ). It happens because the imaginary part of the conductivity increases with the frequency increment and thus the propagation length decreases. Furthermore, it is seen from Fig. 3(b) that the propagation length of higher modes has larger values compared to the lower modes. A large value of the effective index, amounting to 28 e.g., is seen for the chirality of 0.0035 Ω -1 at the frequency of 20 THz for the lower mode. In Fig. 4, we depict the effective index of higher mode for various chiral thicknesses ( t = 600, 650, 700 nm ). In this diagram, the chirality of the chiral substrate is 0.0025 Ω -1 . This figure clearly indicates that the effective index increases as the thickness of the chiral slab decreases. It happens due to the extreme penetration of the electromagnetic fields inside the waveguide for thinner chiral thickness. Fig. 5 illustrates the normalized profile distributions of 𝑅𝑒 [𝐸 𝑥 ] , 𝑅𝑒 [𝐸 𝑦 ] and 𝑅𝑒 [𝐸 𝑧 ] . It can be seen from this figure that higher and lower modes have similar field distributions of 𝑅𝑒 [𝐸 𝑧 ] and 𝑅𝑒 [𝐸 𝑥 ] . The field distribution of 𝑅𝑒 [𝐸 𝑦 ] is the main difference between higher and lower modes , which can distinguish the different types of modes. Fig. 3.
Analytical results of modal properties for the first structure for various chirality values: (a) the effective index as a function of frequency for lower and higher hybrid modes, (b) the propagation length as a function of frequency for lower and higher hybrid modes. The chemical potential of the graphene sheet is
Fig. 4.
The effective index of higher mode for the first structure for various chiral thicknesses ( t = 600, 650, 700 nm ). The chemical potential of the graphene sheet is 0.3 eV and the chirality of the chiral medium is 0.0025 Ω -1 . Fig. 5.
Normalized field distributions of higher and lower modes: (a) Re [E z ], (b) Re [E x ], (c) Re [E y ], with μ g = 0.3 eV, γ = 0.0035 Ω -1 and f = 14 THz. The second structure, as seen in Fig. 6, is a graphene-based waveguide with a chiral substrate, where the ground plane is supposed to be a PEMC. In the second structure, the chiral slab has the thickness of 𝑡 = 750 𝑛𝑚 and its effective index is 𝑛 𝑐 = 1.4 . The admittance parameter of PEMC is considered 𝑀 = 0.5 unless otherwise stated. Similar to the first structure, two values of the chirality (here 𝛾 = 0.003,0.004 𝛺 −1 ) are considered for the chiral substrate to study the effect of chirality on the propagation properties. Fig. 6.
The graphene-based waveguide with a chiral substrate backed by a PEMC layer: (a)
The
3D schematic, (b)
The cross-section of the structure in the z-x plane. 𝑡 𝑃𝐸𝑀𝐶 ɛ , 𝜇 𝜎 𝑃𝐸𝑀𝐶 𝜎 𝑡 𝑧 𝑦 𝑥 (𝒂) 𝑧 (𝒃) ɛ , 𝜇 𝑥 ɛ 𝑐 , 𝜇 𝑐 , 𝛾 ɛ 𝑐 , 𝜇 𝑐 , 𝛾 Fig. 7 illustrates the plasmonic features of the second structure as a function of the frequency. As seen in this figure, the effective index increases as the frequency increases while the propagation length decreases. Similar to the first structure, the propagation length reduces as the frequency increases because the imaginary part of the graphene conductivity (or losses) increases with the frequency increment. Moreover, the lower modes have larger values of the effective index compared to the higher modes. Fig. 8 represents the dependence of the mode confinement on the chemical potential of the graphene. This figure has some remarkable features. First, it allows one to design the proposed structure at the desirable effective index. For instance, to achieve 𝑛 𝑒𝑓𝑓 = 8 , one should design the waveguide with the following parameters: 𝛾 = 0.003 𝛺 −1 , 𝜇 𝑔 =0.2 𝑒𝑉 or 𝛾 = 0.004 𝛺 −1 , 𝜇 𝑔 = 0.35 𝑒𝑉 . Second, the diagram determines the effective index of higher and lower modes at the specific chemical potential. Consider 𝜇 𝑔 = 0.5 𝑒𝑉 . In this chemical doping, the effective index of higher and lower modes for the chirality of 𝛾 = 0.004 𝛺 −1 are 𝑛 𝑒𝑓𝑓 = 3, 6 , respectively. The last feature of this diagram is its ability for determining the propagation range of the plasmonic waves. For instance, hybrid modes do not propagate inside the waveguide for the chemical range of 𝜇 𝑔 > 0.4 𝑒𝑉 for the chirality of 𝛾 = 0.003 𝛺 −1 . As seen in Fig. 8, the chirality loses its effect when the chemical potential of graphene increases, and thus the structure supports only one plasmonic mode (two propagating modes are converged). We called these points “cut-off points”, which happens for the chemical potential of 0.4 eV (for the chirality of 0.003 Ω -1 ) and 0.7 eV (for the chirality of 0.004 Ω -1 ) for this studied structure at the frequency of 10 THz. At these points, two plasmonic modes are converted to radiation modes, i.e. the two lines collapse to a 𝑛 𝑒𝑓𝑓 = 1 line. The cut-off region for the structure of Fig. 6 can be obtained by supposing 𝑘 𝑅,1 > 𝑘
𝐿,1 and setting 𝛽 = 𝑘
𝑅,1 , 𝑞
𝑅,1 =0, 𝑞
𝐿,1 = √𝑘
𝑅,12 − 𝑘
𝐿,12 . Fig. 9 shows the chiral cut-off value as a function of frequency for various chemical potentials. At a specific frequency (for instance, consider 𝑓 = 20 𝑇𝐻𝑧 ), higher cut-off values are obtainable for low values of chemical potential. As a final point, the effective indices of higher and lower modes have been depicted for various values of PEMC admittance (M). As seen in Fig. 10, the mode confinement increases with the admittance increment. Therefore, one way for achieving a desirable effective index is the design of the structure in the appropriate PEMC admittance.
Fig. 7.
Analytical results of modal properties for the second structure for various chirality values: (a) the effective index as a function of frequency for lower and higher hybrid modes, (b) the propagation length as a function of frequency for lower and higher hybrid modes. The chemical potential of the graphene sheet is Fig. 8.
The effective index of hybrid modes as a function of the chemical potential for the second waveguide at the frequency 10 THz.
Fig. 9.
Chiral cutoff values as a function of the frequency for various values of chemical potential.
Fig. 10.
The effective indices of higher and lower modes for the second structure for various admittances of PEMC (
M = 0.5, 0.3, 0.1 ). The chemical potential of the graphene sheet is 0.3 eV and the chirality of the chiral medium is 0.003 Ω -1 . The dashed and solid lines represent lower and higher modes, respectively. Conclusion
A new analytical model is presented for the proposed graphene-based waveguides in this article. The modal properties of plasmonic waves, propagating inside these waveguides, are adjustable by tuning the chirality and the 2 chemical potential of the graphene. As special examples of the general rich structure, two novel waveguides have been proposed and considered. The first waveguide is composed of a graphene sheet deposited on the grounded chiral slab. The ground plane in this structure is supposed to be PEC. A large value of the effective index, amounting to 28 e.g., is obtained at the frequency of 20 THz. The second case is a graphene sheet placed on the chiral slab backed by a PEMC layer. It has been shown that the effective index of this structure can be changed by altering the admittance of the PEMC layer. The hybridization of graphene with chiral materials allows the designer to tune and control SPPs via the chirality and the chemical potential of the graphene. Our presented study can be utilized for designing new kinds of THz devices such as cloaks.
References [1] P. A. D. Gonçalves and N. M. Peres,
An introduction to graphene plasmonics : World Scientific, 2016. [2] F. J. Garcia de Abajo, "Graphene plasmonics: challenges and opportunities,"
Acs Photonics, vol. 1, pp. 135-152, 2014. [3] S. Huang, C. Song, G. Zhang, and H. Yan, "Graphene plasmonics: physics and potential applications,"
Nanophotonics, vol. 6, pp. 1191-1204, 2016. [4] J. Cui, Y. Sun, L. Wang, and P. Ma, "Graphene plasmonic waveguide based on a high-index dielectric wedge for compact photonic integration,"
Optik, vol. 127, pp. 152-155, 2016. [5] I. Toqeer, A. Ghaffar, M. Y. Naz, and B. Sultana, "Characteristics of dispersion modes supported by Graphene Chiral Graphene waveguide,"
Optik, vol. 186, pp. 28-33, 2019. [6] M. B. Heydari and M. H. V. Samiei, "Plasmonic graphene waveguides: A literature review," arXiv preprint arXiv:1809.09937,
IEEE Transactions on Magnetics, vol. 56, pp. 1-10, 2020. [8] M. B. Heydari and M. H. V. Samiei, "Analytical Study of Chiral Multi-Layer Structures Containing Graphene Sheets for THz Applications,"
IEEE Transactions on Nanotechnology, vol. 19, pp. 653-660, 2020. [9] M. B. Heydari and M. H. V. Samiei, "Novel analytical model of anisotropic multi-layer cylindrical waveguides incorporating graphene layers,"
Photonics and Nanostructures-Fundamentals and Applications, vol. 42, p. 100834, 2020. [10] M. B. Heydari and M. H. V. Samiei, "New analytical investigation of anisotropic graphene nano-waveguides with bi-gyrotropic cover and substrate backed by a PEMC layer,"
Optical and Quantum Electronics, vol. 52, pp. 1-16, 2020. [11] M. B. Heydari and M. H. Vadjed Samiei, "An Analytical Study of Magneto-Plasmons in Anisotropic Multi-layer Structures Containing Magnetically Biased Graphene Sheets,"
Plasmonics, vol. 15, pp. 1183-1198, 2020/08/01 2020. [12] M. B. Heydari and M. H. V. Samiei, "Analytical Study of TM-Polarized Surface Plasmon Polaritons in Nonlinear Multi-Layer Graphene-Based Waveguides,"
Plasmonics, arXiv preprint arXiv:2101.02536, arXiv preprint arXiv:2102.10135,
Optik, vol. 134, pp. 60-65, 2017. [16] M. B. Heydari and M. H. V. Samiei, "A Novel Graphene-Based Circulator with Multi-layer Triangular Post for THz Region," arXiv preprint arXiv:2102.02683,
Optik, p. 166457, 2021, https://doi.org/10.1016/j.ijleo.2021.166457 [18] M. B. Heydari and M. H. V. Samiei, "Graphene-Based Couplers: A Brief Review," arXiv preprint arXiv:2010.09462, arXiv preprint arXiv:2011.14767,
Optik, vol. 202, p. 163700, 2020. [21] M. Shalini, "Performance predictions of slotted graphene patch antenna for multi-band operation in terahertz regime,"
Optik, vol. 204, p. 164223, 2020. [22] G. Varshney, S. Debnath, and A. K. Sharma, "Tunable Circularly Polarized Graphene Antenna for THz Applications," Optik, p. 165412, 2020. [23] M. B. Heydari and M. H. V. Samiei, "A Short Review on Graphene-Based Filters: Perspectives and Challenges," arXiv preprint arXiv:2010.07176,
Aip Advances, vol. 4, p. 117103, 2014. [25] C. Gao, Y. Jiang, J. Zhang, W. Cao, X. Gao, X. Yu , et al. , "Graphene-based wideband absorbing screen with radar cross section reduction," in , 2016, pp. 686-689. [26] M. B. Heydari, M. Asgari, M. Zolfaghari, and N. Jafari, "Comparison of Various Full-wave Softwares in Calculating the RCS of Simple Objects,"
International Journal of Mechatronics, Electrical and Computer Technology (IJMEC), vol. 8, pp. 4002-4008, 2018. [27] F. M. Tonelli, V. A. Goulart, K. N. Gomes, M. S. Ladeira, A. K. Santos, E. Lorençon , et al. , "Graphene-based nanomaterials: biological and medical applications and toxicity,"
Nanomedicine, vol. 10, pp. 2423-2450, 2015. [28] H. Xie, T. Cao, F. J. Rodríguez-Lozano, E. K. Luong-Van, and V. Rosa, "Graphene for the development of the next-generation of biocomposites for dental and medical applications,"
Dental Materials, vol. 33, pp. 765-774, 2017. [29] Y. Li, H. Dong, Y. Li, and D. Shi, "Graphene-based nanovehicles for photodynamic medical therapy,"
International journal of nanomedicine, vol. 10, p. 2451, 2015. [30] M. B. Heydari, M. Asgari, N. Jafari, M. Mohtavipour, and M. Zolfaghari, "Various Numerical and Analytical Methods for Studying the Electromagnetic Fields Effects on the Human Body: A Classical Review Article,"
Paramedical Sciences and Military Health, vol. 13, pp. 58-72, 2018. [31] N. Jafari, M. B. Heydari, and M. Asgari, "Simulation of Mobile Radiations in Vicinity of Adult and Child Head,"
Paramedical Sciences and Military Health, vol. 13, pp. 1-7, 2018. [32] N. Jafari and M. B. Heydari, "A Brief Review of New Advances in the Treatment of Spinal Cord Injuries,"
Paramedical Sciences and Military Health, vol. 13, pp. 63-77, 2019. [33] H. Jafari, M. B. Heydari, N. Jafari, and H. Mirhosseini, "Design and Fabrication of Novel Digital Transcranial Electrical Stimulator for Medical and Psychiatry Applications," arXiv preprint arXiv:2009.04411,
IEEE Transactions on Plasma Science, vol. 47, pp. 3320-3326, 2019. [35] M. Heydari, M. Asgari, L. Gharib, A. Keshtkar, N. Jafari, and M. Zolfaghari, "A Novel Augmented Railgun Using Permanent Magnets,"
Advanced Electromagnetics, vol. 8, pp. 99-105, 2019. [36] M. B. Heydari and A. Ahmadvand, "A novel analytical model for a circularly-polarized, ferrite-based slot antenna by solving an integral equation for the electric field on the circular slot,"
Waves in Random and Complex Media, pp. 1-20, 2020. [37] A. Ahmadvand, M. B. Heydari, B. Rejaei, and J. Ghalibafan, "Analytical calculation for input impedance of a circularly ‐ polarized ferrite disk antenna including higher order modes," International Journal of RF and Microwave Computer ‐ Aided Engineering, vol. 29, p. e21783, 2019. [38] D.-G. Zhang and E. K.-N. Yung, "Investigation of waveguide Y-junction circulator with a triangular ferrite post,"
International Journal of Infrared and Millimeter Waves, vol. 21, pp. 123-131, 2000. [39] J. Helszajn, "Circulators Using Triangular and Irregular Hexagonal Planar Resonators," 2008. [40] Y. Akaiwa, "Mode classification of a triangular ferrite post for Y-circulator operation (short papers),"
IEEE Transactions on Microwave Theory and Techniques, vol. 25, pp. 59-61, 1977. [41] M. E. Hines, "Reciprocal and nonreciprocal modes of propagation in ferrite stripline and microstrip devices,"
IEEE Transactions on Microwave Theory and Techniques, vol. 19, pp. 442-451, 1971. [42] M. B. Heydari, M. Zolfaghari, M. Asgari, and N. Jafari, "Analysis of Two modified goubau waveguides at THz frequencies: Conical and elliptical structures,"
Optik, vol. 189, pp. 148-158, 2019. [43] M. B. Heydari, M. Asgari, and N. Jafari, "Novel analytical model for nano-coupler between metal–insulator–metal plasmonic and dielectric slab waveguides,"
Optical and Quantum Electronics, vol. 50, p. 432, 2018. [44] L. Zeng and R. Song, "Quantized chiral anomaly materials cloak,"
Scientific reports, vol. 7, pp. 1-8, 2017. [45] Y. Ye and S. He, "90° polarization rotator using a bilayered chiral metamaterial with giant optical activity,"
Applied Physics Letters, vol. 96, p. 203501, 2010. [46] D. M. Pozar, "Microstrip antennas and arrays on chiral substrates,"
IEEE transactions on antennas and propagation, vol. 40, pp. 1260-1263, 1992. [47] M. Z. Yaqoob, A. Ghaffar, M. Alkanhal, and S. U. Rehman, "Characteristics of light–plasmon coupling on chiral–graphene interface," JOSA B, vol. 36, pp. 90-95, 2019. [48] M. Yaqoob, A. Ghaffar, M. Alkanhal, S. ur Rehman, and F. Razzaz, "Hybrid Surface Plasmon Polariton Wave Generation and Modulation by Chiral-Graphene-Metal (CGM) Structure,"
Scientific reports, vol. 8, pp. 1-9, 2018. [49] V. Gusynin, S. Sharapov, and J. Carbotte, "Magneto-optical conductivity in graphene,"
Journal of Physics: Condensed Matter, vol. 19, p. 026222, 2006. [50] P. Pelet and N. Engheta, "The theory of chirowaveguides,"
IEEE transactions on antennas and propagation, vol. 38, pp. 90-98, 1990. [51] I. V. Lindell, A. Sihvola, S. Tretyakov, and A. Viitanen, "Electromagnetic waves in chiral and bi-isotropic media," 1994. [52] R. Oussaid and B. Haraoubia, "Behavior of a chiral material in terms of a guided wave propagation,"
International Journal of Applied Electromagnetics and Mechanics, vol. 19, pp. 631-635, 2004. [53] J.-F. Dong and J. Li, "Characteristics of guided modes in uniaxial chiral circular waveguides,"
Progress In Electromagnetics Research, vol. 124, pp. 331-345, 2012. [54] C. A. Balanis,
Advanced engineering electromagnetics : John Wiley & Sons, 2012. [55] D. Poole,
Linear algebra: A modern introduction : Cengage Learning, 2014. [56] E. Süli and D. F. Mayers,
An introduction to numerical analysis : Cambridge university press, 2003. [57] A. R. Samant and K. Whites, "Notable features of hybrid modes in a chiral-filled rectangular waveguide,"
IEEE microwave and guided wave letters, vol. 5, pp. 144-146, 1995. [58] S. Mahmoud, "On mode bifurcation in chirowaveguides with perfect electric walls,"
Journal of electromagnetic waves and applications, vol. 6, pp. 1381-1392, 1992. [59] S. Aib, F. Benabdelaziz, C. Zebiri, and D. Sayad, "Propagation in diagonal anisotropic chirowaveguides,"
Advances in OptoElectronics, vol. 2017, 2017. [60] H. Idzuchi, A. Fert, and Y. Otani, "Revisiting the measurement of the spin relaxation time in graphene-based devices,"
Physical Review B, vol. 91, p. 241407, 2015. [61] W. Choi, H. Nishiyama, Y. Ogawa, Y. Ueno, K. Furukawa, T. Takeuchi , et al. , "Relaxation of Plasma Carriers in Graphene: An Approach by Frequency ‐ Dependent Optical Conductivity Measurement,"
Advanced Optical Materials, vol. 6, p. 1701402, 2018. [62] M. T. Mihnev, F. Kadi, C. J. Divin, T. Winzer, S. Lee, C.-H. Liu , et al. , "Microscopic origins of the terahertz carrier relaxation and cooling dynamics in graphene,"
Nature communications, vol. 7, pp. 1-11, 2016. [63] N. Berova, P. L. Polavarapu, K. Nakanishi, and R. W. Woody,
Comprehensive Chiroptical Spectroscopy: Instrumentation, Methodologies, and Theoretical Simulations vol. 1: John Wiley & Sons, 2011. [64] N. Berova, P. L. Polavarapu, K. Nakanishi, and R. W. Woody,
Comprehensive chiroptical spectroscopy: applications in stereochemical analysis of synthetic compounds, natural products, and biomolecules vol. 2: John Wiley & Sons, 2012. [65] M. Hentschel, M. Schäferling, X. Duan, H. Giessen, and N. Liu, "Chiral plasmonics,"
Science advances, vol. 3, p. e1602735, 2017. [66] T. Wu, J. Kessler, and P. Bouř, "Chiral sensing of amino acids and proteins chelating with EuIII complexes by Raman optical activity spectroscopy,"
Physical Chemistry Chemical Physics, vol. 18, pp. 23803-23811, 2016. [67] B. Wang, J. Zhou, T. Koschny, M. Kafesaki, and C. M. Soukoulis, "Chiral metamaterials: simulations and experiments,"
Journal of Optics A: Pure and Applied Optics, vol. 11, p. 114003, 2009. [68] X. Ma, M. Pu, X. Li, Y. Guo, P. Gao, and X. Luo, "Meta-chirality: Fundamentals, construction and applications,"
Nanomaterials, vol. 7, p. 116, 2017. [69] S. Yoo and Q.-H. Park, "Metamaterials and chiral sensing: a review of fundamentals and applications,"
Nanophotonics, vol. 8, pp. 249-261, 2019.vol. 8, pp. 249-261, 2019.