Open access peer-reviewed chapter

Dispersion Properties of TM and TE Modes of Gyrotropic Magnetophotonic Crystals

By Alexander A. Shmat’ko, Viktoria N. Mizernik, Eugene N. Odarenko and Viktor T. Lysytsya

Submitted: May 25th 2017Reviewed: September 26th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.71273

Downloaded: 418

Abstract

This chapter discusses the propagation of TM and TE waves in the one-dimensional gyrotropic magnetophotonic crystals with ferrite and plasma-like layers. Elements of the transfer matrix are calculated in closed analytical form on the base of electrodynamic problem rigorous solution for arbitrary location of the gyrotropic elements on the structure period. Dispersion equation of the layered periodic structure with gyrotropic elements is obtained. Dispersion properties of the structure for TE and TM modes are analyzed for different configurations of magnetophotonic crystals (ferrite and plasma-like layers). Existence areas of transmission bands for surface and bulk waves are obtained. The effect of problem parameters on the dispersion properties of magnetophotonic crystals for TM and TE modes is investigated. Regimes of complete transmission of wave through limited magnetophotonic crystal are analyzed for bulk and surface waves.

Keywords

  • magnetophotonic crystal
  • gyrotropic media
  • dispersion diagrams
  • TE and TM modes
  • bulk and surface waves

1. Introduction

Photonic crystals (PCs) are artificial periodic structures with spatially modulated refractive index in one or more coordinates [1, 2]. Their outstanding optical properties are due to the existence of frequency band gaps where the propagation of electromagnetic waves is impossible. Application of these structures became very attractive for modern optoelectronics which uses the various waveguides, resonators, sensors, and other devices on the basis of PC [3, 4]. Moreover, the control of the PC structure characteristics is the important problem that is usually solved using external electric or magnetic fields. These methods of providing controllability are based on the variation of refractive index of special materials such as liquid crystals and magnetic materials [5, 6]. Since these sensitive materials are anisotropic, then theoretical analysis of their properties is more complicated.

When at least one of the PCs’ unit cell components is a magnetically sensitive (gyrotropic) material, they exhibit unique magneto-optical properties and identified as magnetophotonic crystals (MPCs). Investigations of the MPCs are begun for simplest one-dimensional structures [7, 8]. However, one-dimensional MPCs are the basis elements for various active field-controlling applications so far [9, 10]. Changing of the permeability by external magnetic field is one of the main phenomena that allow developing electronically tuned devices in different frequency bands: filters, circulators and so on [11, 12, 13].

Along with the properties inherent in conventional PCs, these structures have additional optical and magneto-optical properties which considerably expand their functionality. Kerr effect, Faraday rotation and optical nonlinearity can be enhanced in MPC due to light localization within magnetic multilayer. Magneto-optical system with large Faraday or Kerr rotation can be used for effective optical isolators [14, 15], spatial light-phase modulators [16] and magnetic field and current sensor [17] development. Furthermore, one can obtain stronger enhancement of the magneto-optical phenomena due to resonant effects in the MPCs [18], which characterized by specific polarization properties. Using PCs with magneto-optical layers provides possibility of control of optical bistability threshold in structure based on graphene layer [19]. It should be noted that not only magnetic materials are suitable for MPC. Namely, one-dimensional PC with plasma layers can be tunable by external magnetic field [20].

A number of applications of the MPCs are inspired by their nonreciprocal properties. For example, special spatial structure of the MPC layers provides the asymmetry of dispersion characteristics and, as a result, the effect of unidirectional wave propagation [21]. This phenomenon allows enhancing field amplitude in the MPC without any periodicity defects. In this case, the so-called frozen mode regime occurs instead the defect mode one.

One of the unique properties of gyrotropic materials is the possibility of negative values of material parameters under the certain conditions. Usually, these are so-called single-negative media that are divided into epsilon-negative media (plasma) and mu-negative ones (gyrotropic magnetic materials). The term “double-negative media” or “left-handed materials” is used for media with negative values of both permittivity and permeability and often replaced by term “metamaterials.” Application of metamaterials in one-dimensional PC systems results in unusual regularities of bulk and surface wave propagation and is the subject of experimental and theoretical research [22, 23].

Theoretical description of the various types of one-dimensional PCs is usually based on the transfer-matrix method of Abelès [24] that was applied by Yeh et al. to periodic layered media [25]. This method cannot be applied in general case for anisotropic multilayer structures because of mode coupling. However, this is possible in special cases, namely, in two-dimensional model of wave propagation in periodic layered media [26]. This case is considered in this chapter. Such an approach makes it possible to simplify significantly the analysis of physical phenomena in complex layered media with various combinations of gyrotropic and isotropic elements. Moreover using well-known permutation duality principle of Maxwell’s equations results in a reduction of unique combinations number. In turn, this allows better understanding of regularities of bulk and surface wave propagation in one-dimensional MPC and finding new modes for applications in modern microwave, terahertz and optical devices.

2. Formulation and solution of the problem for modes of gyrotropic periodic structures

2.1. Basic relationships

We study electromagnetic wave propagation in periodic structure in general case with bigyrotropic layers (one-dimensional MPC) (Figure 1). Each of two layers on the structure period L = a + b is an anisotropic medium (plasma or ferrite or their combinations). Their permittivity and permeability are characterized by tensor values of standard form [26]:

εj=εjiεaj0iεajεj000εj,μj=μjiμaj0iμajμj000μj.E1

Figure 1:

Schematic of the periodic structure.

For plasma media, the value of the permittivity is tensor, whereas the value of the permeability is scalar. Such media are called electrically gyrotropic. It is opposite for the ferrite media; the permeability is tensor, whereas permittivity is scalar. Such media are usually called magnetically gyrotropic. If the permittivity and permeability are simultaneously described by the tensors (Eq. (1)), such media are called gyrotropic or bigyrotropic. The material parameters included in tensors εjand μjare defined by the value of the external bias magnetic field H0=z0H0, which is directed along Oz axis.

The study of the general case of gyrotropic media with material parameters of form (Eq. (1)) is reasonable, primarily because it allows using the permutation duality principle when obtaining main equations for fields and characteristic Eqs. [26]. According to this principle generalized for gyrotropic media, namely, when simultaneously the substitution of fields EHand material parameters εμis done, the receiving of general equations, from which the equations for magneto-gyrotropic media (ferrite), electro-gyrotropic media (plasma), and gyrotropic media (bianisotropic media) can be obtained easily, turns to be more simple than in each of the mentioned particular cases separately.

Indeed, it follows from Maxwell equations:

rotE=1cBt,rotH=1cDt,E2

where D=εEis inductance vector of electric field and B=μHis inductance vector of magnetic field. Components of these vectors can be written as:

Dx=εEx=εExiεaEy,Dy=εEy=iεaEx+εEy,Dz=εEz=εEz,I8
Bx=μHx=μHxiμaHy,By=μHy=iμaHx+μHy,Bz=μHz=μHz.I9

In general case, one can obtain two connected differential equations for longitudinal components of electromagnetic fields Ez and Hz (along Oz axis and the direction of bias magnetic field H0=z0H0). In two-dimensional case (∂/∂z = 0), these equations can be broken up into two independent Helmholtz equations with respect to the selected longitudinal field components Ez and Hz (TM and TE waves) [26]. Indeed, we can show it for two-dimensional case and harmonic dependence exp(−iωt) of the fields on time t.

Using the relation between field components Ez and Hz via transverse field components, one can obtain the Helmholtz equation for two polarizations Ez and Hz, respectively:

2Ezx2+2Ezy2+k2μεEz=0,2Hzx2+2Hzy2+k2εμHz=0.E3

Here:

μj=μj1μaj2μj2,εj=εj1εaj2εj2.

Eq. (3) describes TM waves (Hx, Hy, Ez), Ez-polarization (s-polarization) and TE waves (Ex, Ey, Hz), Hz-polarization (p-polarization). Therefore, the vector of the electric field for s-polarization is directed perpendicular to the xy plane and the vector of the electric field for the p-polarization is parallel to this plane. It is necessary to use the boundary conditions for the tangential components of the electric and magnetic fields at the media interfaces to solve the boundary electrodynamic problem for the eigenvalues and the eigenfunctions of the Laplace operator. We use the conditions of continuity of the components Ez and Hy for TM waves and Hz and Ey for TE waves.

This analysis shows the important conclusion that in the general case of a gyrotropic medium, the fields of TM and TE waves with respect to the direction of the bias magnetic field H0in the two-dimensional case are divided into two independent solutions of the Maxwell equations. Moreover, Eq. (3) and the expressions for the tangential components of fields Ey and Hy yield the permutation duality principle for the TM and TE waves. These equations are transformed each other after replacing the component of field Ez by Hz, and simultaneously, the effective permeability μ should be replaced by the effective permittivity (−ε) and ε by (−μ). The boundary conditions for the components Ey and Hy also satisfy the permutation duality principle. Hence, it is possible to simplify this general electrodynamic problem, and we can consider only one type of TM or TE waves for any kind of media in order to obtain a solution for another type of wave. Below, we consider the propagation of Ez-polarized (TM) waves in a gyrotropic PC structure.

To determine the eigenvalues and the corresponding eigenfunctions of the two-layer gyromagnetic MPC, we consider Helmholtz equation for the Ez-polarization with the corresponding boundary conditions for the tangential components of the fields Ez and Hy at the interfaces of the periodic structure layers. For the case of Hz-polarization (TE waves), it is necessary to use the permutation duality principle in the solution for TM waves.

The solution of the Helmholtz equation for TM waves for both tangential components of the fields Ez and Hy can be written in the following form:

Ez1xy=an1eiξ1xn1L+bn1eiξ1xn1Leiβy,b<xn1L<LHy1xy=ξ1kμ1an1g1+eiξ1xn1Lbn1g1eiξ1xn1Leiβy,E4
Ez2xy=cneiξ2xnL+dneiξ2xnLeiβy,0<xnL<bHy2xy=ξ2kμ2cng2+eiξ2xnLdng2eiξ2xnLeiβy.E5

Here, n = 1, 2, … is a number of the period; ξj=k2μjεjβ2are the transverse wave numbers in gyrotropic layers in the direction of the Ox axis; β is the longitudinal wave number of the MPC; an, bn, cn, and dn are the wave amplitudes in the layers; and gj±=1±iμajμjξjβ.

Let us note one feature in the expressions for the electromagnetic fields (Eqs. (4) and (5)). The presence of gyrotropy in the layers (μaj ≠ 0) leads to the fact that the distributions of the amplitude of the fields in the layers differ by factors gj±for the forward and backward waves, which propagating in the layers along the direction of periodicity.

To find the dispersion equation that relates the longitudinal wave number β with the structure parameters for a given frequency ω, it is necessary to use the boundary conditions on the interfaces of the layers and the Floquet-Bloch theorem [27, 28] for the periodic structure. Using the conditions of continuity of the Ez and Hy components of the fields at the interfaces x − (n − 1)L = 0 and x = nL, we get matrix equations:

an1bn1=a11a12a21a22cndn,cndn=b11b12b21b22anbn.E6

where

a11=12g1+ξ2ξ1μ1μ2g2+eiξ2L,a12=12g1ξ2ξ1μ1μ2g2eiξ2L,E7
a21=12g1+ξ2ξ1μ1μ2g2+eiξ2L,a22=12g1++ξ2ξ1μ1μ2g2eiξ2L.E8
b11=12g2+ξ1ξ2μ2μ1g1+eiξ2aeiξ1a,b12=12g2ξ1ξ2μ2μ1g1eiξ2aeiξ1a,E9
b21=12g2+ξ1ξ2μ2μ1g1+eiξ2aeiξ1a,b22=12g2++ξ1ξ2μ2μ1g1eiξ2aeiξ1a.E10

Eliminating the coefficients cn, dn in the matrix equations (Eq. (6)), we obtain the relation for the coefficients in identical layers for two neighboring periods of the structure:

an1bn1=a11a12a21a22cndn=a11a12a21a22b11b12b21b22anbn=ABCDanbn.E11

The elements of the ABCD matrix are calculated by the rule of multiplying two matrices. Using the Eqs. (7)(10), we find the elements of the given transfer matrix, namely:

A=cosξ2bi12ξ1ξ2μ2μ1+ξ2ξ1μ1μ2+β2ξ1ξ2μ2μ1μa1μ1μ1μ2μa2μ22sinξ2beiξ1a,E12
D=cosξ2b+i12ξ1ξ2μ2μ1+ξ2ξ1μ1μ2+β2ξ1ξ2μ2μ1μa1μ1μ1μ2μa2μ22sinξ2beiξ1a,E13
B=i12sinξ2bξ2ξ1μ1μ2+ξ1ξ2μ2μ11iβξ1μa1μ1μa2μ2μ1μ22eiξ1a,E14
C=i12sinξ2bξ2ξ1μ1μ2+ξ1ξ2μ2μ11+iβξ1μa1μ1μa2μ2μ1μ22eiξ1a.E15

An important property of the ABCD matrix is the unimodularity property, when the ratio between the elements of the matrix is fulfilled: AD − BC = 1. Using the expressions for the elements of the matrix ABCD, one can show that this condition is satisfied. We note that when ξ1 is a real number, then the elements of the matrix A = D и B = C are pairwise conjugate.

The resulting matrix equation (Eq. (11)), which determines the relationship of the unknown coefficients in two identical layers of different periods of the periodic structure and the Floquet-Bloch theorem, allows us to find the characteristic (dispersion) equation for determining the previously introduced unknown longitudinal wave number β for a wave propagating along gyrotropic layers (along the Oy axis) and the Floquet-Bloch wave number K. According to the Floquet-Bloch theorem in its matrix formulation, one can obtain

ABCDanbn=eiKLanbn.E16

The phase factor eiKL is the eigenvalue of the transfer-matrix ABCD, which is determined from the characteristic equation:

eiKL=12A+D±i112A+D2.E17

The unknown real values of the roots of the characteristic equation have the form:

KTMβ=1Larccoscosξ2bcosξ1a12ξ1ξ2μ2μ1+ξ2ξ1μ1μ2++β2ξ1ξ2μ2μ1μa1μ1μ1μ2μa2μ22sinξ2bsinξ1a.E18

It is easy to show that this expression is transformed to the well-known solution of the dispersion equation for the case of two magnetodielectric layers (μa1 = μa2 = 0) [25].

Note that the longitudinal wave number enters into equation as β squared. This indicates that its absolute value is the same for opposite directions of wave propagation along the Oy axis. That is, the dispersion is the same for forward and backward waves propagating along the layers. However, the field distributions for the case of gyrotropic media in the direction of periodicity for forward and backward waves are different.

Using the permutation duality principle, we found the solutions of the electrodynamic problem for TE wave propagation in the gyrotropic MPC. For this case, we change the material parameters according to the rule μεin the transfer-matrix elements (Eqs. (12)(15)), in the dispersion relation, and in the solution (Eq. (18)). Then, we obtain

ATE=cosξ2bi12ξ1ξ2ε2ε1+ξ2ξ1ε1ε2+β2ξ1ξ2ε2ε1εa1ε1ε1ε2εa2ε22sinξ2beiξ1a,E19
DTE=cosξ2b+i12ξ1ξ2ε2ε1+ξ2ξ1ε1ε2+β2ξ1ξ2ε2ε1εa1ε1ε1ε2εa2ε22sinξ2beiξ1a,E20
BTE=i12sinξ2bξ2ξ1ε1ε2+ξ1ξ2ε2ε11iβξ1εa1ε1εa2ε2ε1ε22eiξ1a,E21
CTE=i12sinξ2bξ2ξ1ε1ε2+ξ1ξ2ε2ε11+iβξ1εa1ε1εa2ε2ε1ε22eiξ1a.E22

Here ξ1=k2ε1μ1β2and ξ2=k2ε2μ2β2.

It is apparent that in this case, one can write

KTEβ=1Larccoscosξ2bcosξ1a12ξ1ξ2ε2ε1+ξ2ξ1ε1ε2++β2ξ1ξ2ε2ε1εa1ε1ε1ε2εa2ε22sinξ2bsinξ1a.E23

In the absence of gyrotropy (εaj = 0), this solution is transformed into the well-known expression for the Bloch wave number for TE waves for the magnetodielectric PC [25].

The given elements of the transfer-matrix ABCD and the solutions of the dispersion equations (Eqs. (18) and (23)) for both s- and p-polarizations (TM and TE waves) are suitable for analysis of wide variety of MPCs with different material parameters: two isotropic layers on the crystal period (dielectric, magnetic, magnetodielectric), one isotropic layer with another anisotropic layer, two gyroelectric layers or gyromagnetic ones or a combination of them, and two gyrotropic layers. Such an abundance of variants makes Eq. (17) universal in terms of analyzing the dispersion characteristics of TE and TM waves and establishing the features of their propagation in various one-dimensional MPCs.

2.2. Eigen regimes of MPCs

We perform the analysis of the features for the propagation of the electromagnetic waves in different MPCs for various eigenmodes. We identified 10 variants of such regimes.

1. The crystal contains two layers of magnetodielectric: εa1 = εa2 = μa1 = μa2 = 0.

2. The crystal contains magnetodielectric layer, εa1 = μa1 = 0, and the layer of a semiconductor plasma: εa2 ≠ 0 and μa2 = 0.

3. The crystal contains magnetodielectric layer, εa1 = μa1 = 0, and the ferrite layer: μa2 ≠ 0 and εa2 = 0.

Taking into account the permutation duality principle, we obtain equations analogous to variant 2.

4. The crystal contains magnetodielectric layer and gyrotropic one: εa1 = μa1 = 0, εa2 ≠ 0, and μa2 ≠ 0.

5. The crystal contains layer of semiconductor plasma: εa1 ≠ 0 and μa1 = 0 and the ferrite layer: εa2 = 0 and μa2 ≠ 0.

6. The crystal contains two layers of semiconductor plasma: εa1 ≠ 0, εa2 ≠ 0 and μa1 = 0, μa2 = 0.

7. The crystal contains two ferrite layers: μa1 ≠ 0, μa2 ≠ 0 and εa1 = 0, εa2 = 0.

Taking into account the permutation duality principle, we obtain equations analogous to variant 6.

8. The crystal contains the plasma layer and gyrotropic layer: εa1 ≠ 0, εa2 ≠ 0 and μa1 = 0, μa2 ≠ 0.

9. The crystal contains the ferrite layer and gyrotropic layer: εa1 = 0, εa2 ≠ 0 and μa1 ≠ 0, μa2 ≠ 0.

Taking into account the permutation duality principle, we obtain equations analogous to variant 8.

10. The crystal contains two gyrotropic layers μa1 ≠ 0, μa2 ≠ 0 and εa1 ≠ 0, εa2 ≠ 0

In this case, one can use solutions Eqs. (18) and (23) for TM and TE waves, respectively.

3. Analysis of the propagation of TE and TM waves in MPCs

3.1. General aspects of wave propagation in MPCs

Let us do a physical analysis of the obtained results and determine general rules of electromagnetic wave propagation in MPC. If the wave number K(β) is real, then the electromagnetic wave propagates in a MPC without attenuation. These particular Floquet-Bloch wave numbers K(β) correspond to the transmission bands and satisfy the condition |cosK(β)L| ≤ 1. On the other hand, a different behavior is observed when the wave number is complex: K(β) = K(β) + iK(β). Such waves cannot propagate in a MPC and thus decay along the direction of periodicity. As a result, the forbidden zones are formed that satisfied the condition |cosK(β)L| > 1. This condition allows easy determination of the imaginary part K(β) of the wave number. When K(β) ≠ 0, then sinK(β)L = 0 and K(β)L = πm (where m = 0, 1, 2… are forbidden zone numbers). As a result, we can obtain the equation for K(β): chKβL=1m12A+D.

The value m = 0 corresponds to the zeroth forbidden zone where the wave number is purely imaginary. It is important to distinguish two cases: the absence (εaj = 0, μaj = 0) and presence (εaj ≠ 0, μaj ≠ 0) of gyrotropy in a MPC. Besides of that, if in jth layer the conditions k2εjμj − β2 < 0 or k2μjεj − β2 < 0 are satisfied, and in one layer ε⊥1 < 0 (orμ⊥1 < 0) and in another layer ε⊥2 > 0(μ⊥2 > 0), then the wave in such layer would decay in amplitude along the x-axis. This wave is a surface wave and delayed in respect to the speed of light.

3.2. Analysis of dispersion characteristics

Let us consider first the case of a MPC in the absence of gyrotropy (εaj = 0 and μaj = 0). This structure is a periodic sequence of magnetodielectric layers. Dispersion equations for these structures were previously considered in a simplified version by many authors. However, the first investigation of such characteristic (dispersion) equations was carried out as early as the nineteenth century by Rayleigh [29] in the solution of the one-dimensional Hill equation. If the medium is periodic, then the latter equation becomes the traditional one-dimensional Helmholtz equation. Further investigations were performed by various researchers, for example, Brillouin [30], Yeh et al. [25], and Bass [31].

The solutions of the dispersion equations for the magnetodielectric periodic structures written out in this section include all possible combinations of the signs and magnitudes of the material parameters. As an example, let us consider dispersion characteristics of one-dimensional magnetodielectric PCs for several combinations of the parameters.

Figure 2 shows dispersion curves for the value of the longitudinal wave number β = 0.9. In this case, the periodic structure consists of two dielectrics (μ1 = μ2 = μ = 1) with positive values of permittivity ε1 = ε∥1 = 2 and ε2 = ε∥2 = 9. The normalized width of the layers is a/L = 0.8 and b/L = 0.2. Solid and dotted curves denote the real and imaginary parts of the Floquet-Bloch wave number, respectively. The imaginary parts of the Floquet-Bloch wave number characterize the degree of wave decay in the forbidden bands. The wave decay is different in forbidden zones. Moreover, the maximum attenuation is observed for the zeroth zone.

Figure 2:

Dispersion characteristics of the one-dimensional photonic crystal.

If we consider a MPC consisting of two different magnetodielectrics, then there are no fundamental differences.

Figure 3 shows the dispersion diagrams for PC with two dielectric layers (μ1 = μ2 = 1) on the structure period and with different signs of the permittivities (ε1 > 0 and ε2 < 0), namely, ε1 = ε∥1 = 2 and ε2 = ε∥2 =  − 6. First the problem of wave propagation at the boundary of two half-spaces from dielectrics with opposite signs of permittivity values was considered by Sommerfeld [32].

Figure 3.

Dispersion diagrams for both polarizations.

The existence of surface waves at the boundaries of PC layers is illustrated in the dispersion diagram (Figure 3) for TE waves (p-polarization). The region of these waves existence is located below the light line k = β(ε1)−1/2 (the solid line on the diagram). Taking into account the identical physical nature of the surface wave existence in PC and Zenneck-Sommerfeld wave for two media [33], this regime can be classified as a modified surface Zenneck-Sommerfeld wave. For this wave, the condition β=kε1ε2ε1+ε2is satisfied approximately, as follows from the dispersion equation. It should be noted that if we consider PC with magnetic layers, then a surface wave will exist for s-polarization.

With the advent of new artificial media (metamaterials) for which the permittivity and permeability are simultaneously negative, optoelectronics devices with new functionalities are developed. In this connection, it is expedient to consider a MPC, one of whose layers on the structure period is a metamaterial (e.g., ε2 < 0, μ2 < 0).

In Figure 4, dispersion diagrams are calculated for the one-dimensional PC with alternative layers of magnetodielectric and metamaterial (ε2 < 0, μ2 < 0). The following parameters were selected: ε1 = 4, μ1 = 3, ε2 =  − 2, μ2 =  − 6, and a/L = 0.5. Figure 5a and b shows the dispersion characteristics for such MPC at the values β = 0 and β = 2π/L, respectively. On the presented graphs, identical dispersion diagrams for both polarizations are seen. Therefore, it is a case of polarization indifference for the forbidden zones and transmission ones. For the presented case, the refractive indices of the layers are the same in absolute value and have phase difference of π, namely, n1=ε1μ1=12and n2=ε2eμ2e=e12. For this case, the complete Floquet-Bloch wave transmission in PC for two polarizations is possible.

Figure 4.

Dispersion diagram of the photonic crystal with metamaterial layer.

Figure 5.

Dispersion characteristics for different values of longitudinal wavenumber: (a) β = 0; (b) β = 2π/L.

It is clear that there are transmission bands only for surface slow waves (modified Zenneck-Sommerfeld waves) that are located below the light line ω12=βc. For bulk waves (k2εjμj > β2), the transmission bands degenerate into curves. These curves are described by equation (2πm)2 = (k2|εjμj| − β2)L2 (m = 1, 2, …) for both polarizations.

Dispersion curves for real K and imaginary K values (Figure 5a and b) show that the wave decay for both polarizations is the same. The propagation of waves is observed at discrete points on which the conditions K = 0 and |cosK(β)L| = 1 are satisfied. These points are located outside the transmission bands of surface waves. The phenomenon of polarization indifference makes it possible to develop devices with a finite number of periods of a MPC in which a complete narrow-band propagation of the wave is possible simultaneously for a whole spectrum of discrete frequencies and both polarizations [34].

We now turn to an analysis of the propagation of waves in MPC consisting of a magnetodielectric and a semiconductor plasma layer [35]. It is advisable to consider two cases in the presence of gyrotropy of the medium (εa2 ≠ 0), when ε⊥2 > 0 and ε⊥2 < 0 for TE waves. We note that the dispersion characteristics for TM waves are the same as for a conventional magnetodielectric. Let us first consider the features of TE wave propagation in MPC when the plasma layer gyrotropy is such that condition ε⊥2 > 0 (εa2 < ε2) is fulfilled. From an analysis of the solutions of the dispersion equation, it follows that in the periodic structure at such a value of the effective permittivity of a plasma medium, there can exist two modes of TE wave propagation. The first mode is bulk waves whose domain of existence is determined by the simultaneous fulfillment of two conditions: k2ε⊥1μ∥1 − β2 > 0 and k2ε⊥2μ∥2 − β2 > 0. The second mode is the propagation mode of gyrotropic surface waves. This mode is observed when three conditions are fulfilled simultaneously: k2ε⊥1μ∥1 − β2 < 0, k2ε⊥2μ∥2 − β2 < 0, and ξ1ξ2ε2ε1+ξ2ξ1ε1ε2+β2ξ1ξ2ε1ε2εa2ε22<0.

If the above conditions are satisfied, then the existence of real values of the Floquet-Bloch wave number other than zero is possible when the condition |cosK(β)L| ≤ 1 is fulfilled for pure imaginary values of transverse wave numbers ξ1 and ξ2.

The results of calculations of the dispersion diagram are shown in Figure 6. The calculation was carried out with the following parameters: a/L = 0.8, μ1 = μ2 = μ∥1 = μ∥2 = 1, ε1 = 2, ε2 = 4, εa2 = 3.2, and ε⊥2 = 1.44. The solid line shows the light line for the first layer on the period of the structure. It can be seen from the figure that when the conditionskε2μ2<β, kε1μ1<β, and ε2ε2β2ξ1ξ2ε1ε2εa2ε22ξ1ξ2ε2ε1+ξ2ξ1ε1ε2>0are satisfied, there exists a region of real values of the Floquet-Bloch wave number for which the gyrotropic surface wave regime is observed (more shaded area of the transmission zone).

Figure 6.

Band structure for TE polarization.

In the case ε⊥2 < 0, as for PC with dielectric layers, the regime of a modified Zenneck-Sommerfeld surface wave can exist in MPC. However, in addition to that wave, there is also a gyrotropic surface wave for other parameters of the problem. Indeed, such a wave exists when both transverse wave numbers in two layers are pure imaginary ξ12<0ξ22<0. Also under these conditions and ε⊥2 < 0, the condition β2ξ1ξ2ε1ε2εa2ε22<ξ1ξ2ε2ε1+ξ2ξ1ε1ε2must be fulfilled, that is, opposite to the condition for positive values of the second-layer effective permittivity.

Figure 7 shows the dispersion diagrams for both polarizations in the case ε⊥2 < 0. Parameters of the problem were chosen as follows: a/L = 0.8, μ1 = μ2 = μ∥1 = μ∥2 = 1, ε1 = 2, ε2 = 4, εa2 = 6.7, and ε⊥2 =  − 7.2. The solid line k = β(ε⊥1μ∥1)−1/2 separates the areas of the existence of bulk waves and surface waves for p-polarization. One area relates to the modified Zenneck-Sommerfeld waves, and the other one relates to the gyrotropic surface wave which also exists for positive values of the effective permittivity, as shown above. We note that in the transmission bands of TE waves in dispersion diagrams, there are two regions in which surface waves propagate.

Figure 7.

Dispersion diagrams of the magnetophotonic crystal containing plasma layers.

The case of MPC with a magnetodielectric and ferrite layer is analogous to that considered earlier by virtue of the permutation duality principle. Let us now consider the following case, when both layers on the structure period are ferrite with different material parameters. In this case, modes of bulk wave’s existence in the transmission bands can be observed when two conditions ξ12=k2μ1ε1β2>0and ξ22=k2μ2ε2β2>0are fulfilled. Regime of surface waves is realized when opposite conditions (ξ12<0, ξ22<0) and also the additional condition are fulfilled:

μ1μ1μ2μ2ξ1ξ2μ2μ1+ξ2ξ1μ1μ2β2ξ1ξ2μ2μ1μa1μ1μ1μ2μa2μ22<0E24

Note that this condition is equally suitable for both positive and negative values of the effective magnetic permeability of ferrite μj. Here, as in the case of the plasma semiconductor layer, the existence of gyrotropic surface waves is possible [36].

In Figure 8, we represent the dispersion diagrams of TM waves for the considered above MPC with two ferrite layers at the period of the structure. In the calculation, the following task parameters were chosen: a/L = 0.5, μ1 = μ∥1 = 2, μ2 = μ∥2 = 3, ε1 = 2, ε2 = 4, μa2 = 2.9,μ⊥2 = 0.197μa1 = 0.1,and μ⊥1 = 1.995.

Figure 8.

Band structure of the ferrite magnetophotonic crystal.

The solid line in the figure, below which the solutions of the dispersion equation are in the regime of surface waves, is determined by the equation k = β(μ⊥1ε∥1)−1/2. Analysis of the dispersion equation solutions and the conditions for the existence of surface waves show that for the chosen values of the material parameters in the ferrite MPC, there exists a range of parameters on the dispersion diagram that corresponds to the regime of propagation of the surface wave at the boundaries of the layers (in this figure it is marked as more shaded area in transmission zone). We emphasize that the regime of the surface wave is realized with positive values of the both magnetic permeabilities of the ferrite layers.

Figure 9 illustrates the case of the existence of a modified Zenneck-Sommerfeld wave for the case when the effective magnetic permeability of one of the layers is negative. Here, the dispersion diagrams are calculated with the following parameters of the problem: a/L = 0.85, μ1 = μ∥1 = 2, μ2 = μ∥2 = 3, ε1 = 1.5, ε2 = 1.8, μa2 = 7.9, μ⊥2 =  − 17.8, μa1 = 0.7, and μ⊥1 = 1.755.

Figure 9.

Dispersion diagrams of the ferrite magnetophotonic crystal for both polarizations.

Figure 10 illustrates the evolution of dispersion diagrams at change of ferrite effective magnetic permeability μ⊥2(μa2 = 5.6; 5.8; 6.1) for a/L = 0.83, μ1 = μ∥1 = 2, μ2 = μ∥2 = 3, ε1 = 4, ε2 = 2, and μa1 = 0.4. A darker shade shows the transmission bands in which the conditions of existence of surface waves are fulfilled. The change of the bias magnetic field value leads to a significant change of the parameters of these areas. An increase in the value of the effective magnetic permeability μ⊥2 owing to the magnitude μa2 results in a change in the width and location of the transmission band of the modified Zenneck-Sommerfeld surface wave.

Figure 10.

Dispersion diagrams for different values of bias magnetic field: (a) μa2 = 5.6; (b) μa2 = 5.8; (c) μa2 = 6.1.

Figure 11 shows dispersion diagrams for the modified Zenneck-Sommerfeld surface wave of MPC with two ferrite layers on the structure period at change of width of the layer b (b/L = 0.3; 0.5; 0.7; 0.8) with μ⊥2 < 0 and for the following parameters of the problem: μ1 = μ∥1 = 2, μ2 = μ∥2 = 3, ε1 = 4, ε2 = 2, μa1 = 0.4, μ⊥1 = 1.92, μa1 = 5.6, and μ⊥2 =  − 7.45.

Figure 11.

Dispersion diagrams for different values of second layer thickness: (a) b/L = 0.3; (b) b/L = 0.5; (c) b/L = 0.7; (d) b/L = 0.8.

Increase of the second-layer width b leads to an expansion of existence area of surface waves with a simultaneous shift of the bandwidth toward the value β = 0 (Figure 11a and b). The largest existence area of surface waves is realized in dispersion diagram for approximately equal values of the layer thicknesses. Further increase of the parameter b is accompanied by a displacement of this area to the opposite direction (Figure 11c and d). Thus the surface wave modes in a ferrite MPC are determined both by the material parameters of the system and by the width of the layers on the period of the structure.

Let us move on the dispersion diagrams for two bigyrotropic layers on the structure period. Taking into account the complete symmetry of the dispersion for TE and TM waves, the case of polarization indifference can be realized in this case. We will show this by example.

Figure 12 shows the dispersion diagrams for TE and TM waves for MPC consisting of two bigyrotropic layers on the period of the structure. Figure 12a corresponds to the following values of the parameters: a/L = 0.85, μj = μj = 2, εj = εj = 2, εa1 = μa1 = 0.5, εa2 = μa2 = 5.2, μ⊥1 = 1.875, and μ⊥2 =  − 11.52. In Figure 12b, only the values of the effective permittivity and permeability of the second layer differ: εa2 = μa2 = 6.8.

Figure 12.

Dispersion diagrams of the bigyrotropic magnetophotonic crystal: (a) εa2 = μa2 = 5.2; (b) εa2 = μa2 = 6.8.

The solid lines in the figures distinguish the area of fast (upper part of figures) and slow wave (the lower part of figures). The complete identity of the dispersion diagrams for the transmission bands of both surface and bulk waves follows from the figures and formulas (18) and (23). By changing the bias magnetic field, it is possible to control the width and location of transmission bands for both polarizations.

Dispersion diagrams in Figure 13 correspond to the case when only the surface wave transmission bands for both polarizations are realized. The parameters of the problem were chosen as follows: a/L = 0.28, μ1 = μ∥1 = 1.5, μ2 = μ∥2 = 2, ε1 = ε∥1 = 1.2, ε2 = ε∥2 = 1.8, εa1 = 4.95, εa2 = 1.7, μa1 = 0.7, and μa2 = 5. The polarization sensitivity of the MPC is realized in this case. Only the surface waves with certain polarization can propagate through the periodic structure for defined values of parameters k and β.

Figure 13.

Transmission bands for surface waves in the bigyrotropic magnetophotonic crystal.

Therefore, we have considered the main features of the propagation of TE and TM waves in various magnetophotonic gyrotropic crystals. Important application of this one-dimensional PC theory is the problem of the electromagnetic waves scattering by a structures with limited number of periods.

4. Scattering of a plane wave by a MPC

In this section, the scattering of the plane wave on gyromagnetic MPC with N periods is considered. When a p-polarized plane wave is scattered on MPC with a limited number of periods, the problem is divided into three stages. At the first stage, the problem of scattering of a plane wave on the first gyromagnetic layer of MPC is solved. In the second stage, the coupling between the field coefficients of the first and last layers of the MPC is used in the problem for the MPC modes. And, finally, in the third stage, the problem of wave transmission from the last layer of the structure to the surrounding area is considered. Following [37], we write out the final expression for the reflection and transmission coefficients for gyrotropic MPC:

R=n21M11+n22M21k11+n21M12+n22M22k21n11M11+n12M21k11+n11M12+n12M22k21E25
T=1n11M11+n12M21k11+n11M12+n12M22k21E26

Here, matrix M is the Nth power of an ABCD matrix. Other notations correspond to Ref. [37]. Figure 14a shows the dependences of the transmission coefficient modulus on the normalized frequency in the case of normal wave incidence (β = 0) on MPC with 20 periods in the regime of bulk waves. Calculation parameters are the same as in Figure 7.

Figure 14.

The transmittance vs frequencies for the case of normal incidence of wave: (a) spectral characteristic; (b) fine structure of spectral characteristic in second transmission band; (c) fine structure of spectral characteristic in third transmission band.

There are three transmission bands in the frequency range under consideration. Each of these bands contains resonances which observed with respect to the frequency of the complete transmission of the wave (Figure 14b and c). Frequency resonances correspond to different modes of the periodic structure. The number of modes is determined by the number of periods of the structure (N − 1).

Figure 15 depicts the frequency dependences of the transmission coefficient modulus in the regime of the surface waves. In this case the incident angle of wave is greater than the angle of total internal reflection. We can see one transmission band in this case. Inset in Figure 15 shows enlarged frequency dependence within this band. Complete propagation is observed for each resonant frequency of modes of the limited periodic structure.

Figure 15.

The transmittance vs. frequencies for the case of surface wave mode.

5. Conclusions

The electrodynamic problem is solved for the proper TE and TM waves of a MPC with two gyrotropic layers. The elements of the transmission matrix, the dispersion equation, and its solution are obtained analytically. An analysis of the dispersion properties of TE and TM waves for MPC is carried out, and features of the existence of fast and slow waves are revealed. Different regimes of gyrotropic surface waves are found. The conditions for the existence of surface waves are established for positive and negative values of the permittivity and permeability. Analytic expressions for the reflection and transmission coefficients for a limited MPC are obtained, and their analysis is performed for the regime of bulk and surface waves. Complete transmission of the wave through this structure is realized at resonant frequencies that correspond to different spatial distributions of the mode field in limited MPC.

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Alexander A. Shmat’ko, Viktoria N. Mizernik, Eugene N. Odarenko and Viktor T. Lysytsya (December 20th 2017). Dispersion Properties of TM and TE Modes of Gyrotropic Magnetophotonic Crystals, Theoretical Foundations and Application of Photonic Crystals, Alexander Vakhrushev, IntechOpen, DOI: 10.5772/intechopen.71273. Available from:

chapter statistics

418total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Structure-Induced Ultratransparency in Photonic Crystals

By Jie Luo and Yun Lai

Related Book

First chapter

Introductory Chapter: Molecular Dynamics: Basic Tool of Nanotechnology Simulations for “Production 4.0” Revolution

By Alexander V. Vakhrushev

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us