# Properties of the Band Gaps in 1D Ternary Lossy Photonic Crystal Containing Double-Negative Materials

Authors: Alireza Aghajamali; Maryam Akbarimoosavi; Mahmood Barati
Journal: Advances in Optical Technologies
 (2014)

---

## Abstract

Theoretically, the characteristics matrix method is employed to investigate and compare the properties of the band gaps of the one-dimensional ternary and binary lossy photonic crystals which are composed of double-negative and double-positive materials. This study shows that by varying the angle of incidence, the band gaps for TM and TE waves behave differently in both ternary and binary lossy structures. The results demonstrate that, by increasing the angle of incidence for the TE wave, the width and the depth of zero-n¯, zero-μ, and Bragg gap increase in both ternary and binary structures. On the other hand, the enhancement of the angle of incidence for the TM wave contributes to reduction of the width and the depth of the zero-n¯ and Bragg gaps, and they finally disappear for incidence angles greater than 50° and 60° for the binary structure and 40° and 45° for the ternary structures, respectively. In addition, the details of the edges of the band gaps variations as a function of incidence angle for both structures are studied.


---

## Body

## 1. Introduction

Photonic crystals (PCs) are artificial dielectric or metallic structures in which the refractive index changes periodically. The periodic structure of the PCs contributes to photonic band gap (forbidden range of frequencies). Interference of the Bragg scattering is considered as a cause of these phenomena, the so-called Bragg gap. PCs have attracted much interest due to their novel electromagnetic waves characteristics and important scientific and engineering applications and have received attentions by many researchers in recent decades [1–8]. With possibility of producing metamaterials, PCs containing metamaterials, namely, “Metamaterial Photonic Crystals” (MPCs), have been made. The technological relevance of the MPCs is of great importance. These structures have attracted considerable attention for their various applications. Designing dielectric mirrors, effective waveguides, filters, and perfect lenses are clear-cut examples. The metamaterials are utilized to design one-dimensional (1D) binary and ternary MPCs with layers of double-negative (DNG) and double-positive (DPS) materials. In these structures, various band gaps can be seen in the transmission spectrum. When the average refractive index of the MPC is equal to zero, the new band gap which is called zero-n- gap appears in the transmission [9], which possesses some unique properties that distinguish itself from an ordinary Bragg gap. The properties and characteristics of the zero-n- gap of 1D binary PC structure consisting of DNG and DPS materials were investigated in recent years, such as bandwidth, the depth, and the central frequency of the gap [9–15].In 1D binary MPC, in addition to the zero-n- and Bragg gaps, new gaps appear in transmission spectra for both TE and TM polarized waves. In the TE mode, this gap appears at the frequency where permeability of DNG material vanishes, called the zero-μ gap, but in the TM mode, this gap arises at the frequency where permittivity of DNG material goes to zero, called the zero-ε gap. These two new types of gaps are found for the oblique incidence angles but are absent for normal incidence [16–18].The novel idea of this research is to include the loss factor which is an unavoidable consequence of double-negative materials to investigate optical properties of 1D binary and ternary MPCs. In this regard, the transmission of electromagnetic waves through a 1D ternary lossy PC consisting of layers with double-negative and double-positive materials is studied. The effects of the incidence angle and polarization on the characteristics of the zero-n-, zero-μ, zero-ε, and Bragg gaps are investigated. Furthermore, the results are compared with the similar results in the 1D binary lossy PC. The paper is organized as follows: the MPC structure’s design, the permittivity and permeability of the DNG layer, and theoretical formulation (characteristic matrix method) are described in Section 2, the numerical results and discussions are presented in Section 3, and the paper is concluded in Section 4.


## 2. Structures Design and Characteristic Matrix Method

The 1D MPC structure under study which is located in air is constituted by alternative layers of DNG and DPS materials, where the DNG material is dispersive and dissipative. We consider the 1D binary MPC with periodic structure of(AB)N and the 1D ternary MPC with periodic structure of (ABC)N, where A, B, and C denote three different materials. N is number of the lattice periods, and di, εi, μi, and ni (i=A,B,C) are thickness, permittivity, permeability, and refractive index of the layers, respectively.The permittivity and permeability of the DNG (layerA) in the microwave region are complex and are given as [9, 10, 19–21]
(1)εA(f)=1+520.92-f2-ifγ+10211.52-f2-ifγ,μA(f)=1+320.9022-f2-ifγ,
where f and γ are frequency and damping frequency, respectively, given in GHz. Various details of the real parts of the permittivity, permeability, and refractive index of layers A, εA′, μA′, and nA′ versus frequency for γ=0.2×10-3 GHz are shown in Figure 1. As seen from the figure, there exist three different regions. In the first region, f<3.13 GHz, εA′ and μA′ simultaneously are negative (DNG material) and nA′<0. For 3.13<f<3.78 GHz, in the second region, εA′<0, but μA′>0 (single-negative material extending to the epsilon-negative material), and similar to the first region, nA′<0 again. Finally, in the third region, f>3.78 GHz, both εA′ and μA′ are positive (DPS material) and nA′>0.Figure 1
Real parts of permittivity (dotted line), permeability (dashed line), and refractive index (solid line) of layersA, εA′, μA′, and nA′ versus frequency for a small loss factor.The calculation is based on the characteristic matrix method [22], which is most effective to analyze the transmission properties of PCs. The characteristic matrix M[d] for TE wave at incidence angle θ0 from vacuum to a 1D PC structure is given by [22]
(2)M[d]=∏i=1k[cos⁡γi-ipisinγi-ipisinγicos⁡γi],
where k=2 and 3 for binary and ternary PCs, respectively. γi=(ω/c)nidicos⁡θi, where  c is speed of light in vacuum, θi is the angle of refraction inside the layer with refractive index ni, and pi is given as
(3)pi=εiμicos⁡θi,
where
(4)cos⁡θi=1-n02sin2θ0ni2.The characteristic matrix forN period structure is therefore [M(d)]N. The transmission coefficient of the multilayer is given by
(5)t=2p0(m11+m12ps)p0+(m21+m22ps),
where mij(i,j=1,2) are the matrix elements of [M(d)]N and
(6)p0=n0cos⁡θi,ps=nscos⁡θs.The transmissivity of the multilayer is given by(7)T=psp0|t|2.The above formulations can be applied to TM wave by simple replacements ofpi, p0, and ps as follows:
(8)pi=μiεicos⁡θi,p0=cos⁡θ0n0,ps=cos⁡θsns.


## 3. Numerical Results and Discussion

Based on the theoretical model described in the previous section, the transmission spectrum of the presented lossy MPC structures was calculated. In the study of the 1D binary MPC, consisting of DNG (layerA) and DPS (layer B) materials, equations (1) and (2) are used for the permittivity, εA, and the permeability, μA, of layer A. Layer B is assumed to be a vacuum layer with εB=μB=1 (nB=1). The thickness of layers A and B is chosen as dA=6 mm and dB=12 mm, respectively, and the total number of lattice periods is set as N=16 [10].The transmission spectra of TE and TM polarized waves for the binary structure at various angles of incidence and forγ=0.2×10-3 GHz are shown in Figures 2(a) and 2(b), respectively. For oblique incidence, the zero-μ and zero-ε gaps appear in the transmission spectrum for TE and TM waves, respectively, as reported in [17, 18]. The gaps appear at the frequencies where the sign of the permeability or permittivity of DNG materials changes. As it is seen in Figure 2, the width of the gaps increases as the incidence angle increases, while the left edge of the gaps remains nearly unchanged. In addition, for TE waves, the width and the depth of the zero-n- and the Bragg gap increase as the incidence angle increases, as reported in [9–14]. On the contrary, for TM waves, the width and the depth decrease as the incidence angle increases (Figure 2(b)). Moreover, the figures show that the central frequencies of the zero-n-, as reported in our previous work [10], and the Bragg gap shift to the higher frequencies as the angle of incidence increases.Transmission spectra for the 1D binary MPC structure for different angles of incidence withγ=0.2×10-3 GHz: (a) TE and (b) TM waves.
(a)
(b)In the next part, the band gaps of the 1D ternary MPC are investigated. The binary structure used before is modified by introducing a third layer ofSiO2 (layer C) with refractive index of nC=1.46 and thickness of dC=6 mm in each lattice period. The other parameters are kept the same as in the binary structure. The transmission spectra of TE and TM polarized waves for the ternary structure at various angles of incidence and for γ=0.2×10-3 GHz are shown in Figures 3(a) and 3(b), respectively. As it is seen, the width of the zero-ε and zero-μ gaps increases by increasing the incidence angle. Although, for TE wave, the width and the depth of the zero-n- and the Bragg gaps increase as the incidence angles increases. For TM waves, the behaviour is completely different such that the width and the depth of the band gaps decrease when the angle of incidence increases. However, the central frequency of the zero-n- and Bragg gaps behaves the same as in the binary structure where they shift to the higher frequencies as the angle of incidence increases.Transmission spectra for the 1D ternary MPC structure for different angles of incidence withγ=0.2×10-3 GHz: (a) TE and (b) TM waves.
(a)
(b)In this part, the band gaps in binary and ternary structures are compared. As it is seen from Figures2 and 3, the width of the band gaps decreases, when another dielectric layer is added to the lattice period of binary structure. Furthermore, the zero-n- and Bragg gaps frequencies of the ternary structure are lower than the corresponding frequencies in the binary structures, while the zero-ε and zero-μ gaps appear in the same frequencies in both structures.The lower (fL) and higher (fH) frequencies of the zero-n- gap in the binary and ternary structures, as a function of incidence angle for TE and TM waves, are shown in Figures 4(a) and 4(b), respectively. As it is clearly seen, the lower frequency of the gap for the ternary structure is smaller than that in the binary structure for both TE and TM modes for different incidence angles. It is also interesting to note that the width of the gap in the ternary structure for TE wave is sensitive to the angle of incidence and increases as the incidence angle increases, while the left edge of the band gap is nearly invariant. Moreover, the width of the zero-n- gap of both binary and ternary structures is very sensitive to the incidence angle for TM wave, and the gap disappears for incidence angles greater than 50° in the binary structure as reported in [10, 12]. The gap disappears for incidence angles greater than 40° for the ternary structure.The lower and higher frequencies of the zero-n- gap as a function of angle of incidence for both binary and ternary structures: (a) TE and (b) TM waves.
(a)
(b)The lower and the higher frequencies of the zero-μ and zero-ε gaps as a function of incidence angle, for both binary and ternary structures, are shown in Figures 5(a) and 5(b), respectively. The width of the zero-μ and zero-ε gaps increases faster for the binary structure than the ternary ones. As it is clearly seen, the gaps appear at the same frequency for both binary and ternary structures; in addition, the figures indicate the upper edge of the gaps, so the width of the gaps is very sensitive to the incidence angle, while the lower edge of the gaps is nearly insensitive.The lower and higher frequencies of the (a) zero-μ and (b) zero-ε gaps as a function of angle of incidence for both binary and ternary structures.
(a)
(b)ThefL and fH frequencies of the Bragg gap of the binary and ternary structures, as a function of incidence angle for TE and TM waves, are shown in Figures 6(a) and 6(b), respectively. It is seen that the Bragg gap like the zero-n- gap appears in the lower frequencies for the ternary structure. Moreover, fL and fH shift to the higher frequencies as the incidence angle increases, and as it is seen from Figure 6(b), the gap for TM polarized wave disappears for incidence angles greater than 60° for the binary and 45° for the ternary structures.The lower and higher frequencies of the Bragg gap as a function of angle of incidence for both binary and ternary structures: (a) TE and (b) TM waves.
(a)
(b)


## 4. Conclusion

The numerical results show that including the loss factor in the permittivity and permeability of the DNG layer the changes in transmission spectrum of TE wave in 1D ternary and binary MPCs behave similar to different incidence angles. The zero-n-, zero-μ, and the Bragg gaps become wider and deeper as the angle of incidence increases. In addition, by enhancing the angle of incidence, the lower frequencies of the zero-n- and zero-μ gaps are nearly invariant in both ternary and binary MPCs, but the Bragg gap increases. Moreover, the upper edge of all three band gaps increases. Such behaviors were not clearly observed in the zero-n- gap. The transmission spectrum for TM polarized wave is rather different. It is found that, in both structures, as the angle of incidence increases, the width and the depth of the zero-ε gap increase, while for the zero-n- gap and for the angles of incidence greater than 40° for the ternary and 50° for the binary structures, the zero-n- gap disappears. Similar behavior is observed for the Bragg gap for the angles of incidence greater than 45° and 60° for binary and ternary structures, respectively. For TM waves, the lower frequency of the zero-n- and Bragg gaps increases, but the zero-ε gap decreases as the angle of incidence increases in both binary and ternary structures. The higher edge of the zero-ε and Bragg gaps increases, while it decreases for the zero-n- gap. It is interesting to note that the difference between the lower and higher frequencies of all the gaps for TM wave in the binary MPC is larger than that in the ternary MPC as all the figures show. To recapitulate what was said before, (i) the width of the band gaps decreases as a dielectric layer is added to the binary structure, (ii) the zero-n- and Bragg gaps appear at the lower frequencies in the ternary structure, and (iii) the zero-ε and zero-μ gaps appear at the same frequencies in both structures. Finally, the results of the present study could be employed in designing new edge filters, waveguides, and other optical devises in microwave engineering. In addition, detailed analysis of the band gaps in ternary metamaterial photonic crystal will certainly provide useful information for designing new types of transmission filter. This filter has a superior feature of ternary structure, which is not seen in the usual binary photonic crystals.



---

*Source: 129568-2014-02-13.xml*