TEMagnetostatic SurfaceWaves in Symmetric Dielectric Negative Permittivity Material Waveguide

Nonlinear magnetostatic surface wave in a slab waveguide structure has been investigated. The design consisted of dielectric film between two thick nonlinear nonmagnetic negative permittivity material (NPM) layers. A dispersion relation for TE nonlinear Magnetostatic surface waves (NMSSWs) has been derived into the proposed structure and has been numerically investigated. Effective refractive index decreases with thickness and frequency increase have been found. Effective refractive index decrease with optical nonlinearity increase and switching to negative values of effective refractive index at a certain value of optical nonlinearity have been found. This meant that the structure behaved like a left-handed material over certain range. We found that the power flow was changing by changing the operating frequency, the dielectric film thickness, and the optical nonlinearity. Also, the effective refractive index and power flow attained constant values over certain values of dielectric constant values.


Introduction
Great interest is focused on the propagation of electromagnetic waves in artificial materials, and particularly on materials with negative index of refraction: materials which exhibited both negative permeability and permittivity over a certain range of frequencies.In those materials, there were the wave vector, the electric field, and the magnetic field form a left-handed system.Thus, they were called left-handed materials (LHMs).
A group of researchers at the University of San Diego was able to synthesize an artificial dielectric medium (metamaterials).They were able to demonstrate that those materials exhibit both negative dielectric permittivity and magnetic permeability simultaneously over a certain range of frequencies [1].By that realization, the prediction of Veselago [2] in his pioneer paper that electromagnetic propagation in an isotropic medium with negative dielectric permittivity ε(ω) < 0 and negative permeability μ(ω) < 0 that could exhibit unusual properties were realized.In such (LHMs) there appeared the electric field vector E, the magnetic field vector H, and the wave vector kform a left-hand orthogonal set.Those recent demonstrations on the existence of the LHM resulted in a wide open to unique possibilities in the design of a novel type of device based on electromagnetic wave propagation in those materials, but in a nonconventional way.
Recent studies by Assa'd et al. [6,7] had proposed a new class of materials based on Shadrivov's and others work, which was the nonlinear nonmagnetic negative permittivity materials (NPMs).The nonlinear nonmagnetic NPM behaves like a metal with a better advantage; we can control its physical characteristics and the operating frequency range.This controllability is not widely available when use certain class of metal compared to NPMs [6].
Since most communication devices (e.g., waveguides and microstrips) include dielectric materials, which invited us to investigate a slab waveguide structure of dielectric film between two thick layers of nonlinear nonmagnetic negative permittivity material (NPM).
The propagation of TE waves NMSSW in slab waveguide structure was investigated.They were dielectric film between two layers of nonlinear nonmagnetic negative permittivity material (NPM).
This paper has been organized as follows.Section 2 derived the dispersion relation of the surface waves in nonlinear nonmagnetic NPM dielectric nonlinear nonmagnetic NPM structure, and the power flow.In Section 3, we discussed the numerical results.Section 4 was solely devoted to the conclusion.

Theory
The dispersion relation for TE wave propagation in the x-axis with propagation wave constant β is represented in the form e ik0 (βx−ct) , where β = k/k 0 , k is the complex propagation constant, and k 0 is the free space wave number which equals ω/c, where c is the velocity of light, and ω is the applied angular frequency.
In Figure 1, the linear dielectric film of thickness t and dielectric constant ε is bounded by two layers of nonlinear nonmagnetic NPM in the regions z < 0 and z > t, and its dielectric constant is ε NL eff .Here we consider the Kerr-like nonlinearity of the dielectric in the composite material [8][9][10], that is, where ω p is the plasma frequency, ω is the applied frequency, and ε D = ε 1 + αE 2 y , ε 1 is the linear part of the dielectric constant, and α is the nonlinear coefficient.The permeability of the nonlinear nonmagnetic NPM is considered μ = 1 [6,11].The electromagnetic field components are (2)
where z 0 is the position of the maximum of the field component in the nonlinear cover, and The magnetic field components in the nonlinear cover at z > t are and at z < t 2.2.In Linear Dielectric Region (0 < z < t).Also, substitution of (2) into Maxwell's equation yields the following linear differential equation to satisfy in the linear dielectric film: The solution of ( 6) is given by where The magnetic field components in the linear dielectric are Advances in Condensed Matter Physics The dispersion relation can be found by matching the field components at the interface z = 0 and z = t, that is, where and u = tanh(k 1 z 0 ) is the optical nonliearaty.
Power Flow.The power flux of the TE surface waves along the direction of propagation can be found by integrating the Poynting vector E y H z dz = P 1NL + P D + P 2NL , where

Numerical Results and Discussion
The dispersion relation, (11), numerically solved to find the complex effective wave index β as a function of the angular frequency ω, and the optical nonlinearity u for different values of dielectric film thicknesses (t = 2, 4, and 6 μm).
The power flux for the structure under investigation has been investigated for the film thicknesses.Besides that, we investigated the effect of the dielectric constant of the dielectric material.The parameters of the nonlinear nonmagnetic negative permittivity NPM and linear dielectric are adjusted so that the parameter ε NL eff is negative in the same frequency range, which lies between 1.3 ∼ 2.2 GHz.The parameters that were used in carrying out the numerical calculations are [6,13]: In Figure 2, we plotted the effective refractive index, β, versus the frequency range at optical nonlinearity, u = 0.27, at different dielectric film thicknesses.It is noticed that the effective refractive index is smoothly decreasing with frequency increase and negative slope.In the middle range of frequencies the effective refractive index switches to negative values with positive slope.After that, at higher frequencies the effective refractive index switches back to positive values with negative slope again.Thus, the structure behaved like left-handed material (LHM) in the middle range frequencies.The slope of the dispersion relation represents the group velocity.When the slope has a negative gradient this means that the group velocity is negative, and the structure behaves like LHM.These properties are due to the fact that the Poynting vector in the nonlinear nonmagnetic  NPM is antiparallel to the wave number, which characterized by negative dielectric constant [14,15].
The effective refractive index β of the structure versus the nonlinearity u at different film thickness and at angular frequency ω = 1.7 GHz is shown in Figure 3.We notice that the effective refractive index, β, of the structure smoothly decreases with the nonlinearity u, despite the small region in the nonlinearity (u = 0.25 ∼ 0.28) in which drastic change happens in the values of the effective refractive index, it turns to negative values.We noticed that the effective refractive index is sensitive to the dielectric film thickness.The effective refractive index decreases with dielectric film thickness increase for all nonlinearity values.
Figure 4 shows the normalized power flow (P/P 0 ), where P 0 = 1/(2αε 0 ω), versus the operating frequency with optical nonlinearity u = 0.27.At low frequencies in the range the power flow decreases with frequency increase.Then the power flow switches to negative values in the middle range frequencies.At high frequencies in the range the power flow switches to positive values.The power flow is sensitive to dielectric medium slab thickness.The power flow decreases with dielectric film thickness increase.
In Figures 5(a) and 5(b) (for clarity), we plotted the power versus the optical nonlinearity at ω = 1.7 GHz operating frequency and dielectric constant, ε 2 = 3.The power flow decreases with optical nonlinearity increase for all dielectric film thicknesses.
In Figures 6 and 7, we explored the effect of the dielectric constant on the effective refractive index of the structure and power flow.For some operating frequency range the effective refractive index of the structure and power flow have almost constant values.The effective refractive index and power flow decreases with dielectric film thickness increase for all dielectric constant values taken in the range.

Conclusion
We analytically studied the TE surface waves in a slab waveguide.It contained dielectric film between two nonlinear nonmagnetic negative permittivity material (NPM) thick layers.A dispersion relation for TE surface waves has been derived and numerically investigated.There found out that the wave effective refractive index was changing in value and sign depending on the operating frequency and the dielectric slab thickness.The effective refractive index and power flow were sensitive to optical nonlinearity and dielectric slab thicknesses; they decreased with nonlinearity and thickness increase.We noticed that the effective refractive index and power flow attained constant values, at certain values of nonlinearity and operating frequency, versus dielectric slab dielectric constant.

3 −Figure 2 :
Figure 2: The effective nonlinear refractive index versus the angular frequency at nonlinearity u = 0.27, for different dielectric slab thicknesses.

Figure 7 :
Figure 7: The normalized power flux versus dielectric constant at operating angular frequency and optical nonlinearity, ω = 1.7 GHz, u = 0.27, for different dielectric film thicknesses.