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  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  have studied negative refraction in LHMs. I. V. Shadrivov  has investigated nonlinear guided waves in LHMs. N. Garcia et al  have shown that LHMs don’t make a perfect lens. Kong  has provided a general formulation for the electromagnetic wave interaction with stratified metamaterial structures. M. M. Shabat et al  have discussed Nonlinear TE surface waves in a left-handed material and magnetic super lattice waveguide structure. I. Kourakis et al  have investigated a nonlinear propagation of electromagnetic waves in negative-refraction index LHM. M. F. Ubeid et al  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  and C. Sabah et al  have estimated high reflection coatings of multilayered structure. Oraizi et al  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 , 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.
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).
We use Maxwell’s equation as done by  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, E1 = E2 and H1 = H2 and so on. This yields six equations with six unknown parameters [11, 13, 16]:
Where k1x = k2x = k3x = k4x = Snell’s law and kℓz is given by:
Letting A1 = 1 and solving the obtained equations for the unknown parameters enables us to calculate the reflection and transmission coefficients B1 and A4 [11, 16]. The reflected power R and the transmitted power T are given by [11, 16]:
Where B1* and A4* are the complex conjugate of B1 and A4 respectively.
Where, Ploss is the loss power due to losses in LHM.
THE EFFECTIVE PARAMETERS
where ωep and ωmp are the electric and magnetic plasma frequencies, ωeo and ωmo are the electric and magnetic resonance frequencies. Fe and Fm are the scaling filling parameters. γe and γm are the electric and magnetic dissipation factors. We have used the following parameters appearing in : ωmp = 2π 10.95 GHz, ωmo = 2π 10.1 GHz, Fm = .26, ωep = 2π 13.3 GHz, ωeo = 2π 10.3 GHz, Fe = .37. The electric and dissipation factors are selected to be equal to each other (γ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 fo = 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.
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):
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, Fm,ωmp and ωmo:
The permeability μ3:
The permeabilities μ1, μ2 and μ4:
The constants Fe, ωep and ωeo:
The permittivity ε3:
The refractive index n3:
To choose the negative sign of n3 if ε3 and μ3 are both negative:
The refractive indices n1, n2 and n4:
The wave vector components in the z direction:
The equations 3-8 are inserted:
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 30o 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 n3 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 n3 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 n3 is imaginary. In 9-10 GHz band the transmitted power is not zero because both ε3(ω) and μ3(ω) are positive and then n3 is real. This means that electromagnetic waves will only propagate in a medium that has a real index of refraction . 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.
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