#### 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-, 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- 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.

#### 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- 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- 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- 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-, 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 and the 1D ternary MPC with periodic structure of , where , , and denote three different materials. is number of the lattice periods, and , , , and () are thickness, permittivity, permeability, and refractive index of the layers, respectively.

The permittivity and permeability of the DNG (layer ) in the microwave region are complex and are given as [9, 10, 19–21] where 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 , , , and versus frequency for GHz are shown in Figure 1. As seen from the figure, there exist three different regions. In the first region, GHz, and simultaneously are negative (DNG material) and . For GHz, in the second region, , but (single-negative material extending to the epsilon-negative material), and similar to the first region, again. Finally, in the third region, GHz, both and are positive (DPS material) and .

The calculation is based on the characteristic matrix method [22], which is most effective to analyze the transmission properties of PCs. The characteristic matrix for TE wave at incidence angle from vacuum to a 1D PC structure is given by [22] where and for binary and ternary PCs, respectively. , where is speed of light in vacuum, is the angle of refraction inside the layer with refractive index , and is given as where

The characteristic matrix for period structure is therefore . The transmission coefficient of the multilayer is given by where are the matrix elements of and

The transmissivity of the multilayer is given by

The above formulations can be applied to TM wave by simple replacements of , , and as follows:

#### 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 (layer ) and DPS (layer ) materials, equations (1) and (2) are used for the permittivity, , and the permeability, , of layer . Layer is assumed to be a vacuum layer with (). The thickness of layers and is chosen as mm and mm, respectively, and the total number of lattice periods is set as [10].

The transmission spectra of TE and TM polarized waves for the binary structure at various angles of incidence and for 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- 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-, as reported in our previous work [10], and the Bragg gap shift to the higher frequencies as the angle of incidence increases.

**(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 of (layer ) with refractive index of and thickness of 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 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- 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- and Bragg gaps behaves the same as in the binary structure where they shift to the higher frequencies as the angle of incidence increases.

**(a)**

**(b)**

In this part, the band gaps in binary and ternary structures are compared. As it is seen from Figures 2 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- 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 () and higher () frequencies of the zero- 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- 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 in the binary structure as reported in [10, 12]. The gap disappears for incidence angles greater than for the ternary structure.

**(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.

**(a)**

**(b)**

The and 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- gap appears in the lower frequencies for the ternary structure. Moreover, and 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 for the binary and for the ternary structures.

**(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-, 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- 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- 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- gap and for the angles of incidence greater than for the ternary and for the binary structures, the zero- gap disappears. Similar behavior is observed for the Bragg gap for the angles of incidence greater than and for binary and ternary structures, respectively. For TM waves, the lower frequency of the zero- 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- 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- 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.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of the paper.

#### Acknowledgment

Alireza Aghajamali would like to acknowledge his gratitude to Parisa Shams for her help and useful discussion.