## 1.

## INTRODUCTION

Metamaterials (sometimes termed left-handed materials (LHMs)) are materials whose permittivity *ε* and permeability *μ* are both negative and consequently have negative index of refraction. These materials are artificial and theoretically discussed first by Veselago [1] over 40 years ago. The first realization of such materials, consisting of split-ring resenators (SRRs) and continuous wires, was first introduced by Pendry [2, 3]. Regular materials are materials whose *ε* and *μ* are both positive and termed right handed materials (RHMs). R. A. Shelby et al [4] have studied negative refraction in LHMs. I. V. Shadrivov [5] has investigated nonlinear guided waves in LHMs. N. Garcia et al [6] have shown that LHMs don’t make a perfect lens. Kong [7] has provided a general formulation for the electromagnetic wave interaction with stratified metamaterial structures. M. M. Shabat et al [8] have discussed Nonlinear TE surface waves in a left-handed material and magnetic super lattice waveguide structure. I. Kourakis et al [9] have investigated a nonlinear propagation of electromagnetic waves in negative-refraction index LHM. M. F. Ubeid et al [10] have demonstrated the effect of negative permittivity and permeability on the transmission of electromagnetic waves through a structure containing left-handed material. H. Cory et al [11] and C. Sabah et al [12] have estimated high reflection coatings of multilayered structure. Oraizi et al [13] have obtained a zero reflection from multilayered metamaterial structures. In this paper we consider a structure consisting of LHM and dielectric slabs inserted in vacuum. A plane polarized wave is obliquely incident on it. Tow cases of the LHM are considered, loss-less case and loss case. Maxwell’s equations are used to determine the electric and magnetic fields in each region. Then, Snell’s law is applied and boundary conditions are imposed at each interface to obtain a number of equations with unknown parameters. The equations are solved for the unknown parameters by a MAPLE software to calculate the reflection and transmission coefficients. These coefficients are used to determine the reflected, transmitted and loss powers of the structure. The effect of many parameters like frequency, angle of incidence etc. on the mentioned powers are studied in details by changing the dissipation factor of the LHM. Throughout the computations, the frequency dependence of permittivity and permeability of the LHM are taken into account. The numerical results are in agreement with the law of conservation of energy given by [11, 14, 15]. It is also noticed that the numerical results of Figure 4 is similar to Figure 8.6 obtained by [14], this is another evidence for validity of the performed computations. The propagation of electromagnetic waves through a structure containing LHM and dielectric is considered by many previous studies. Our work is concentrated on the powers of the structure and the role of properties of the dissipation factor of LHM layer. Thus, the behavior of the powers for the incident waves against the incident angle, the frequency and LHM thickness is computed and presented in numerical results with the emphasis on the dissipation factor of the LHM.

## 2.

## THEORY

We consider four regions each with permittivity ε_{ℓ} and permeability μ_{ℓ}, where ℓ represents the region order. Region 1 and 4 are vacuums (ε_{o}, μ_{o}), Region 2 is a regular dielectric (ε_{2}, μ_{2}), Region 3 is a metamaterial (ε_{3}(ω), μ_{3}(ω)). A polarized plane wave in Region 1 incident on the plane z = 0 at some angle *θ* relative to the normal to the boundary (see Figure 1).

The electric field in each region is [7, 11]:

We use Maxwell’s equation as done by [16] to find corresponding magnetic field :

Where *A _{ℓ}* and

*B*are the amplitudes of forward and backward traveling waves. is the wave vector inside the material and

_{ℓ}*n*is the refractive index of it. Matching the boundary conditions for and fields at each layer interface, that is at z = 0,

_{ℓ}*E*

_{1}=

*E*

_{2}and

*H*

_{1}=

*H*

_{2}and so on. This yields six equations with six unknown parameters [11, 13, 16]:

Where *k*_{1x} = *k*_{2x} = *k*_{3x} = *k*_{4x} = Snell’s law and *k _{ℓz}* is given by:

Letting *A*_{1} = 1 and solving the obtained equations for the unknown parameters enables us to calculate the reflection and transmission coefficients *B*_{1} and *A*_{4} [11, 16]. The reflected power R and the transmitted power T are given by [11, 16]:

Where *B*_{1}* and *A*_{4}* are the complex conjugate of *B*_{1} and *A*_{4} respectively.

The law of conservation of energy is given by [11, 14, 15]:

Where, P_{loss} is the loss power due to losses in LHM.

## 3.

## THE EFFECTIVE PARAMETERS

For the LHM in region 3 we employ a dispersive one with *ε*_{3} and *μ*_{3} appeared in [2, 3, 16]:

where *ω _{ep}* and

*ω*are the electric and magnetic plasma frequencies,

_{mp}*ω*and

_{eo}*ω*are the electric and magnetic resonance frequencies.

_{mo}*F*and

_{e}*F*are the scaling filling parameters.

_{m}*γ*

_{e}and

*γ*

_{m}are the electric and magnetic dissipation factors. We have used the following parameters appearing in [16]:

*ω*= 2

_{mp}*π*10.95 GHz,

*ω*= 2

_{mo}*π*10.1 GHz,

*F*= .26,

_{m}*ω*= 2

_{ep}*π*13.3 GHz,

*ω*= 2

_{eo}*π*10.3 GHz,

*F*= .37. The electric and dissipation factors are selected to be equal to each other (

_{e}*γ*=

_{e}*γ*=

_{m}*γ*) and two cases of the LHM are considered, loss-less case (

*γ*= 0) and loss_case (

*γ*≠ 0). The operating frequency is assumed to be f

_{o}= 11 GHz. This frequency is chosen by an arbitrary decision, but it must be in the frequency band where the permittivity and permeability of the LHM are both simultaneously negative. Region 1, 2 and 4 in Fig. 1 are assumed to be loss-less and the permeability’s of them are equal to the permeability of free space. The refractive index of the dielectric in region 4 is selected to be 1.883. The thickness of each of LHM and dielectric slabs is equal to one half-wavelength long at the operating frequency.

## 4.

## CONSTRUCTION OF THE PROGRAMME

In this section the MAPLE programme is constructed as follows (note that, only the blue and red items appear in the programme):

• Input

The thickness of the slabs (a, d), the velocity of light (c = .3 Gm/sec) as well as the constant *π*:

The frequency and the angle of incidence:

The constants *γ*_{m}, *γ*_{e}, F_{m},ω_{mp} and ω_{mo}:

The permeability μ_{3}:

The permeabilities μ_{1}, μ_{2} and μ_{4:}

The constants F_{e}, ω_{ep} and ω_{eo}:

The permittivity ε_{3}:

The refractive index n_{3}:

To choose the negative sign of n_{3} if ε_{3} and μ_{3} are both negative:

The refractive indices n_{1}, n_{2} and n_{4}:

The wave vector components in the z direction:

The equations 3-8 are inserted:

• Output

## 5.

## NUMERICAL RESULTS

In figure 2 we show a calculation of reflected, transmitted and loss powers of the considered structure as a function of frequency for three values of dissipation factor of the LHM (*γ* = 0, .1 and .2 GHz) and for 30^{o} angle of incidence. The frequency is changed between 9 GHz and 12 GHz, because the simultaneously negative permittivity and permeability can be realized in this range in all cases, according to Eq (12) and (13). As confirmed from Fig. 2, the powers show increasing, decreasing and oscillatory behaviors in different frequency ranges. The interpretation of this is as follows: the transmission is null above 11.5 GHz frequency, where all incident radiation is being reflected. In this band *ε*_{3}(*ω*) and *μ*_{3}(*ω*) are different in signs, and then n_{3} is imaginary, this is because . In 10.3-11.5 GHz band the transmission is very good. In this band *ε*_{3} (*ω*) and *μ*_{3} (*ω*) are both negative, and then n_{3} is real. At 10.1 and 10.2 GHz frequencies the transmitted power is zero because *μ*_{3} (*ω*) is positive and *ε*_{3}(*ω*) is negative and then n_{3} is imaginary. In 9-10 GHz band the transmitted power is not zero because both *ε*_{3}(*ω*) and *μ*_{3}(*ω*) are positive and then n_{3} is real. This means that electromagnetic waves will only propagate in a medium that has a real index of refraction [16]. This interpretation is presented for loss-less case and can be applied on the loss case of LHM.

Figure 3 illustrates the variation of reflected, transmitted and loss powers with the angle of incidence. The angle of incidence is changed between 0° and 90° to realize all possible angles of incidence. Clearly the reflected power increases while the transmitted and loss powers decrease with the angle of incidence. At 90° the reflected power is maximum while the transmitted and loss powers are minimum at that angle for any value of dissipation factor. The role of the dissipation factor is clear at angles below 90°. The reflected and transmitted powers decrease while the loss power increases with the dissipation factor for any angle below 90°.

Figure 4 demonstrates the effect of thickness of LHM on the reflected, transmitted and loss powers respectively at 30° angle of incidence. The slab thickness is changed from zero to 30 mm It is noticed from the figure that the reflected power changes periodically for all dissipation factors. The transmitted power shows the same property for *γ* = 0 GHz, it shows oscillatory decreasing when *γ* = .1 and .2 GHz. On the other hand the loss power shows oscillatory increasing behavior with the increasing values of the dissipation factor.

## 6.

## CONCLUSIONS

A MAPLE program is designed to study the transmission and reflection of electromagnetic waves by a structure containing LHM with effect of the dissipation factor. The followed method is based on Maxwell’s equations and matching the boundary conditions for the electric and magnetic fields at each layer interface. The frequency dependence of *ε* and *μ* of the LHM is taken into account. The dependence of the reflected, transmitted and loss powers of the considered structure on various parameters have been investigated to observe the effect of the dissipation factor. As can be seen from the numerical results, if the dissipation factor changes, the behaviors of the powers will be affected from this change. Consequently the dissipation factor has an important role in the variations of the powers. The law of conservation of energy given by [11, 14, 15] is satisfied by our results. The discussed problem is useful for applications which require controlling of reflected and transmitted powers like antenna radome, microwave, millimeter wave and optical devices