Nonlinear Effects of Electromagnetic TM Wave Propagation in Anisotropic Layer with 0 Kerr Nonlinearity

The problem of electromagnetic TM wave propagation through a layer with Kerr nonlinearity is considered. The layer is located between two half-spaces with constant permittivities. This electromagnetic problem is reduced to the nonlinear boundary eigenvalue problem for ordinary differential equations. It is necessary to find eigenvalues of the problem propagation constants of an electromagnetic wave . The dispersion equation DE for the eigenvalues is derived. The DE is applied to nonlinear metamaterial as well. Comparison with a linear case is also made. In the nonlinear problem there are new eigenvalues and new eigenwaves. Numerical results are presented.


Introduction
Problems of electromagnetic wave propagation in nonlinear waveguide structures are intensively investigated during several decades.First known studies about nonlinear optics' problems are given in the monographs 1, 2 .Propagation of electromagnetic wave in a layer and a circle cylindrical waveguide are among such problems.Phenomena of electromagnetic wave propagation in nonlinear media have original importance and also find a lot of applications, for example, in plasma physics, microelectronics, optics, and laser technology.There are a lot of different nonlinear phenomena in media when an electromagnetic wave propagates, such as self-focusing, defocusing, and self-channeling 1-5 .
Investigation of nonlinear phenomena leads us to solve nonlinear differential equations.In some cases it is necessary to solve nonlinear boundary eigenvalue problems NBEPs , which rarely can be solved analytically.One of the important nonlinear phenomenan is the case when the permittivity of the sample depends on electric field intensity.And one of the simplest nonlinearities is a Kerr nonlinearity 4, 6, 7 .When we Advances in Mathematical Physics speak about NBEPs we mean that differential equations and boundary conditions nonlinearly depend on the spectral parameter and also the differential equations nonlinearly depend on the sought-for functions.These facts do not allow to apply well-known methods of spectral problems' investigation.
Here we consider electromagnetic TM wave propagation in a layer with Kerr nonlinearity.Perhaps, the papers 6, 7 were the first studies where some problems of electromagnetic wave propagation are considered in a strong electromagnetic statement.Propagation of polarized electromagnetic waves in a layer and in a circle cylindrical waveguide with Kerr nonlinearity is considered in this paper.When one says that the permittivity ε is described by Kerr law this means that for an isotropic material ε ε const α|E| 2 , ε const is the constant part of the permittivity ε; α is the nonlinearity coefficient; , where E E x , E y , E z is an electric field.Below we consider an anisotropic case some results were presented in 8 .The first approximation for eigenvalues of the problem is presented in 9 .
Problems of electromagnetic wave propagation in a linear layer with constant permittivity and in a linear circle cylindrical waveguide were deeply studied many years ago, see, for example, 10 .Such problems are formulated as boundary eigenvalue problems for ordinary differential equations.Indeed, the main interest in this problem is the value of the spectral parameter eigenvalues which corresponds to the propagating wave.If an eigenvalue is known it is easy to solve differential equations numerically.Otherwise numerical methods cannot be successfully applied.However, in nonlinear cases it is often paid more attention to solve the differential equations see, e.g., 11-13 .Though the first problem is to find eigenvalues therefore to find Des, from the mathematical standpoint the DE is an equation with respect to the spectral parameter.Analysis of this equation allows us to make conclusions about problem's solvability, eigenvalues' localization, and so forth.In most cases the equations of the problem cannot be integrated in an explicit form.Of course, if one has the explicit solutions of the differential equations it is easy to derive the DE.Therefore, when the equations cannot be integrated things do not go to a DE.However, in some cases the DE can be found in an explicit form and it is not necessary to have explicit solutions of differential equations.
Let us discuss in detail the case of Kerr nonlinearity.The work 4 contains a wide range of details of third-order nonlinear electromagnetic TE and TM guided waves.Problems of surface wave propagation along the interface between two semi-infinite linear or/and nonlinear media were studied completely see the results in 4 .At the same time we should notice that problems of wave propagation in a nonlinear layer that is located between two semi-infinite linear or/and nonlinear media are much more difficult than and cannot be reduced to the problems where surface waves are considered only at the interface between two semi-infinite linear or/and nonlinear media.Propagation of TE waves was more studied.The work 14 is devoted to the problem of electromagnetic wave propagation in a nonlinear dielectric layer with absorption and the case of Kerr nonlinearity is considered separately.One of the most interesting works about propagation of TE waves in a layered structure with Kerr nonlinearity is the paper 15 .Also the reader can see the work 16 , where a layer with Kerr nonlinearity without absorption is considered.
The case of TM wave propagation in a nonlinear medium is more complicated.This is due to the fact that two components of the electric field make the analysis much harder 17 .
In the work 18 a linear dielectric layer is considered.The layer is located between two half-spaces.The half-spaces are filled by nonlinear medium with Kerr nonlinearity.This problem for TE waves is solved analytically 19, 20 .For the TM case in 18 obtained DE is an algebraic equation.It should be noticed that in 18 authors simplify the problem.Earlier in 21 the DE is obtained with other simplifying assumption authors take into account only one component E x of the electric field .Later in 22 it is proved that the dominating nonlinear contribution in the permittivity is proportional to the transversal component E z .In the works 11 propagation of TM waves in a nonlinear half-space with Kerr nonlinearity is considered.Formal solutions of differential equations in quadratures are obtained.In the paper 11 DEs are presented for isotropic and anisotropic media in a half-space with nonlinear permittivity.The DEs are rational functions with respect to the value of field's components at the interface.Authors found the first integral of the system of differential equations so called a conservation law .This is also very interesting work to study, another way to simplify the problem pointed out in 23 .
In the case of TE wave you can see the papers 24-26 .Propagation of TM wave in terms of the magnetic component is studied in 12, 13 .The paper 21 is devoted to the question from physical standpoint why it is possible to take into account only one component of the electric field in the expression for permittivity in the case of TM waves in a nonlinear layer.The results are compared with the case of TE waves.
The most important results about TM wave propagation in a layer with Kerr nonlinearity system of differential equations, first integral and a circle cylindrical waveguide system of differential equations were obtained in 6, 7 .In some papers e.g., 12 polarized wave propagation in a layer with arbitrary nonlinearity is considered.However, DEs were not obtained and no results about solvability of the boundary eigenvalue problem were obtained as well.The problem of TM wave propagation in a layer with Kerr nonlinearity is solved at first for a thin layer and then for a layer of arbitrary thickness 27-29 .Theorems of existence and localization of eigenvalues are proved in 30, 31 .Some numerical results are shown in 8, 9 .
In this paper the DE is an equation with additional conditions.Only for linear media when permittivity is a constant in a layer or in a circle cylindrical waveguide the DEs are sufficiently simple but even for these cases the DEs are transcendental equations .For a nonlinear layer the DE is quite complicated nonlinear integral equation, where the integrand is defined by implicit algebraic function.It should be stressed that in spite of the fact that the DE is complicated it can be rather easily solved numerically.
This DE allows to study both nonlinear materials and nonlinear metamaterials.It should be noticed that in this paper materials with nonlinear permittivity and constant positive permeability are studied.But it is not difficult to take into account the sign of the permeability.
Problems of propagation of TE and TM waves in a nonlinear circle cylindrical waveguide are also close to the problem considered here.These problems are more complicated in comparison with corresponding problems in nonlinear layers.And even in the case of Kerr nonlinearity the results are not so complete as in the case in layers 30, 32, 33 .

Statement of the Problem
Let us consider electromagnetic wave propagation through a homogeneous anisotropic nonmagnetic dielectric layer.The layer is located between two half-spaces: x < 0 and x > h in Cartesian coordinate system Oxyz.The half-spaces are filled with isotropic nonmagnetic media without any sources and characterized by permittivities ε 1 ≥ ε 0 and ε 3 ≥ ε 0 , respectively, where ε 0 is the permittivity of free space.Assume that everywhere μ μ 0 is the permeability of free space see Figure 1 .It should be noticed that conditions ε 1 ≥ ε 0 , ε 3 ≥ ε 0 are not necessary.They are not used for derivation of DEs, but they are useful for DEs' solvability analysis.

Advances in Mathematical Physics
The electromagnetic field depends on time harmonically 6 : where ω is the circular frequency; E, E , E − , H, H , H − are real functions.Everywhere below the time multipliers are omitted.Form complex amplitudes of the electromagnetic field where the continuity condition for the tangential field components on the media interfaces x 0, x h, and the radiation condition at infinity: the electromagnetic field exponentially decays as |x| → ∞ in the domains x < 0 and x > h.
The permittivity inside the layer is described by the diagonal tensor where and a, b, ε 2 > max ε 1 , ε 3 are positive constants below the solutions are sought under more general conditions .It does not matter what a form ε yy has.Since ε yy is not contained in the equations below for the TM case, it should be noticed that ε describes tensor Kerr nonlinearity.
When a b we obtain scalar Kerr nonlinearity.Moreover, chosen nonlinearity satisfies the condition ∂εxx/∂E 2 z ∂εzz/∂E 2 x .This equation is satisfied by almost every known nonlinear Kerr mechanism, such as electronic distortion, molecular orientation, electrostriction, and Kerr nonlinearities described within the uniaxial approximation mentioned in the paper 11 .The case when ε 2x ε 2z is studied in 8 .Pay heed to the fact that the problem considered here is not studied in 31 .
The solutions to the Maxwell equations are sought in the entire space.

TM Waves
Let us consider TM waves E x , E z , H y are functions of three spatial variables.It is easy to show that the components of the fields do not depend on y.Waves propagating along medium interface z depend on z harmonically.This means that the fields components have the form where γ is the spectral parameter of the problem.So we obtain from system 2.3 6

3.3
where • ≡ d/dx.The following equation can be easily derived from the previous system: Differentiating 3.4 and using the second and the third equations of system 3.3 we obtain

3.5
Let us denote k 2 0 : ω 2 μ 0 ε 0 and perform the normalization according to the formulas Advances in Mathematical Physics Z x : E z , X x : iE x and omitting the tilde symbol, from system 3.5 we obtain

3.6
It is necessary to find eigenvalues γ of the problem that correspond to surface waves propagating along boundaries of the layer 0 < x < h.We seek the real values of the spectral parameter γ such that real solutions X x and Z x to system 3.6 exist.Indeed, in this case is a function on z.In this case the components E x and E z depend on z, but it contradicts to the choice of E x x and E z x .So we can consider only real values of γ.
We consider that

3.7
Also we assume that max ε 1 , ε 3 < γ 2 < min ε 2x , ε 2z .This two-sided inequality naturally appears for the problem in a layer with a constant permittivity tensor.
Functions X, Z are supposed to be sufficiently smooth due to physical nature of the problem It is clear that system 3.6 is an autonomous one.System 3.6 can be rewritten in a normal form.This system in the normal form can be considered as a dynamical system with analytical with respect to X and Z right-hand sides.Of course, in the domain where these right-hand sides are analytical with respect to X and Z, it is well known see, e.g., 34 that the solutions X and Z of such a system are analytical functions with respect to the independent variable as well.This is an important fact for DEs' derivation.
In accordance with the radiation condition we obtain

4.1
We assume that γ 2 − ε 1 > 0; otherwise it will be impossible to satisfy the radiation condition.
In the domain x > h we have ε ε 3 .From system 3.6 we obtain X γZ, Z γ −1 γ 2 − ε 3 X.In accordance with the radiation condition we obtain

4.2
Here for the same reason as above we consider that γ 2 − ε 3 > 0.
Constants A and B in 4.1 and 4.2 are defined by transmission conditions and initial conditions.
Inside the layer 0 < x < h system 3.6 takes the form

4.3
Differentiating the second equation and substituting its right-hand side instead of lefthand side into the first equation we can rewritten system 4.3 in the following form:

4.4
Now system 4.4 is written in a normal form.If the right-hand sides are analytic functions with respect to X and Z, then the solutions are analytic functions with respect to its independent variable.
Dividing the first equation in system 4.4 to the second one we obtain the ordinary differential equation Equation 4.5 can be transformed into a total differential equation.
Its solution first integral of system 4.4 can be easily found and be written in the following form: where C is a constant of integration.

Transmission Conditions and the Transmission Problem
Tangential components of an electromagnetic field are known to be continuous at media interfaces.In this case the tangential components are H y and E z .Hence, we obtain From the continuity conditions for the tangential components of the fields E and H and using 3.4 we obtain

5.2
where The constant E h z : E z h 0 is supposed to be known initial condition .Let us denote X 0 : X 0 , X h : X h , Z 0 : Z 0 , and Then from conditions 5.2 we obtain

5.3
In accordance with 3.6 , 3.7 inside the layer Then for x h, using 5.2 , we obtain from 5.4

5.5
From 5.5 we obtain the equation with respect to X h : Under taken assumptions in regard to ε 2 and a the value a −1 ε 2 aZ 2 h > 0. Hence, this equation has at least one real root, which is considered the root can be find explicitly by using Cardanus-Ferrari formula 35 .
Using first integral 4.6 at x h, we find the value C X h : C| x h from the equation

5.7
In order to find the values X 0 and Z 0 it is necessary to solve the following system this system is obtained using formula 5.4 at x 0 and the first integral at the same point :

5.8
It is easy to see from the second equation of this system that the values X 0 and Z 0 can have arbitrary signs.At the same time from the first equation of this system we can see that X 0 and Z 0 must be positive or negative simultaneously.
Normal components of electromagnetic field are known to be discontinues at media interfaces.And it is the discontinuity of the first kind.In this case the normal component is E x .It is also known that the value εE x is continuous at media interfaces.From the above and from the continuity of the tangential component E z it follows that the transmission conditions for the functions εX and Z are εX x 0 0, εX x h 0, Z x 0 0, Z x h 0, 5.9 where f x x 0 lim x → x 0 −0 f x − lim x → x 0 0 f x denotes a jump of the function f at the interface.
We also suppose that functions X x and Z x satisfy the condition 5.10

Dispersion Equation
Introduce the new variables Using new variables rewrite system 4.4 , and 4.6 where constant C is equal to the constant C in 4.6 .
In order to obtain the DE for the propagation constants it is necessary to find the values η 0 , η h .
It is clear that η 0 γ X 0 /Z 0 τ 0 , η h γ X h /Z h τ h .Taking into account that γ 2 X x τ x εX x and using formulas 5.2 , 5.3 , it is easy to obtain that It is easy to see that the right-hand side of the second equation of system VI is strictly positive.This means that the function η x monotonically increases on interval 0, h .Taking into account 6.4 we obtain that the function η x cannot be differentiable on the entire interval 0, h .This means that the function η x has a break point.Let x * ∈ 0, h be the break point.From 6.3 it is obvious that x * is such that τ * τ x * is a root of the equation It is natural to suppose that the function η x on interval 0, h has several break points x 0 , x 1 , . . ., x N .The properties of function η x imply where i 0, N.
Taking into account our hypothesis we will seek the solutions on each interval 0, x 0 , x 0 , x 1 , . . ., x N , h :
Substituting x 0, x x i 1 , and x x N into equations in 6.7 into the first, the second, and the third, resp., and taking into account 6.5 , we find constants c 1 , c 2 , . . ., c N 1 : where i 0, N − 1.
Using 6.8 we can rewrite 6.7 in the following form: where 0 ≤ x ≤ x 0 , x i ≤ x ≤ x i 1 , x N ≤ x ≤ h, respectively, and i 0, N − 1.
Introduce the notation T : ∞ −∞ wdη.It follows from formula 6.9 that 0 < x i 1 − x i T < h, where i 0, N − 1.This implies the convergence of the improper integral it will be proved in other way below .Now consider x in 6.9 such that all the integrals on the left side vanish i.e., x x 0 , x x i , and x x N , and sum all equations in 6.9 .We obtain 6.10 Finally we obtain where η 0 , η h are defined by formulas 6.4 .Expression 6.11 is the DE, which holds for any finite h.Let γ be a solution of DE 6.11 and an eigenvalue of the problem.Then, there are eigenfunctions X and Z, which correspond to the eigenvalue γ.The eigenfunction Z has N 1 zeros on the interval 0, h .
Notice that improper integrals in DE 6.11 converge.Indeed, function τ τ η is bounded as η → ∞ since τ γ −2 ε 2x bX 2 aZ 2 , and X, Z are bounded. Then where α > 0, β > 0 are constants.It is obvious that improper integral ∞ −∞ dη/ αη 2 β converges.Convergence of the improper integrals in 6.11 in inner points results from the requirement that the right-hand side of the second equation of system VI is positive.
The first equation of system VI jointly with the first integral can be integrated in hyperelliptic functions.The solution is expressed in implicit form by means of hyperelliptic integrals.This is the simple example of Abelian integrals.The inversion of these integrals is hyperelliptic functions and they are solutions of system VI .Hyperelliptic functions are Abelian functions, which are meromorphic and periodic functions.Since function η is expressed algebraically through τ, therefore, η is a meromorphic periodic function.This means that the break point x * is a pole of function η.

Generalized Dispersion Equation
Here we derive the generalized DE, which holds for any real values ε 2 .In addition the sign of the right-hand side of the second equation in system VI and condition max ε 1 , ε 3 < γ 2 < ε 2 are not taken into account.These conditions appear in the case of a linear layer and are used for derivation of DE 6.11 .Though on the nonlinear case it is not necessary to limit the value γ 2 from the right side, at the same time it is clear that γ is limited from the left side, since this limit appears from the solutions in the half-spaces.

Advances in Mathematical Physics 13
Now we assume that γ satisfies the following two-sided inequality: 7.1 Using first integral 6.3 it is possible to integrate formally any of the equations of system VI .As earlier we integrate the second equation, we cannot obtain the solution on the entire interval 0, h , since function η x can have break points, which belong to 0, h .It is known that function η x has break points only of the second kind η is an analytical function .
Assume that function η x on interval 0, h has N 1 break points x 0 , x 1 , . . ., x N .
It should be noticed that where i 0, N, and signs ± are independent and unknown.
Taking into account the previous, solutions are sought on each interval 0, x 0 , x 0 , x 1 , . . ., x N , h : where 0 ≤ x ≤ x 0 , x i ≤ x ≤ x i 1 , and x N ≤ x ≤ h, respectively, and i 0, N − 1.
From 7.3 , substituting x 0, x x i 1 , and x x N into the first, the second, and the third equations in 7.3 , respectively, we find required constants c 1 , c 2 , . . ., c N 1 : where 0 ≤ x ≤ x 0 , x i ≤ x ≤ x i 1 , and x N ≤ x ≤ h, respectively, and i 0, N − 1.
From formulas 7.5 we obtain that Expressions 0 < x i 1 − x i < h < ∞ imply that under the assumption about the break point existence the integral on the right side converges and wdη > 0. In the same way, from the first and the last equations of 7.5 we obtain that x 0 These considerations yield that the function w η has no nonintegrable singularities for η ∈ −∞, ∞ .And also this proves that the assumption about finite number break points is true.Now, setting x x 0 , x x i , and x x N into the first, the second, and the third equations in 7.5 , respectively, we have that all the integrals on the left sides vanish.We add all the equations in 7.5 to obtain

7.9
Advances in Mathematical Physics 15 From 7.9 we obtain It follows from formulas 7.6 that and it is necessary to choose the infinities of different signs.
Thus we obtain that Hence Now we can rewrite 7.10 in the following form: Let T ≡ ∞ −∞ wdη; then we finally obtain where η 0 , η h are defined by formulas 6.4 .Expression 7.14 is the DE, which holds for any finite h.Let γ be a solution of DE 7.14 and an eigenvalue of the problem.Then, there are eigenfunctions X and Z, which correspond to the eigenvalue γ.The eigenfunction Z has N 1 zeros on the interval 0, h .It should be noticed that for every number N 1 it is necessary to solve two DEs: for N 1 and for − N 1 .

Note 1.
If there is a certain value γ 2 * , such that some of the integrals in DEs 6.11 or 7.14 diverge at certain inner points this simply means that the value γ 2 * is not a solution of chosen DE and the value γ 2 * is not an eigenvalue of the problem.
Note 2. It is necessary to emphasize that this boundary eigenvalue transmission problem essentially depends on the initial condition Z h .The transmission problem for a linear layer does not depend on the initial condition.If the nonlinearity function is a specific one, then in some cases it will be possible to normalize the Maxwell equations in such a way that the transmission problem does not depend on initial condition Z h explicitly it is possible e.g., for Kerr nonlinearity in a layer and in a circle cylindrical waveguide .Once more we stress the fact that the opportunity of such normalization is an exceptional case.What is more, in spite of the fact that in certain cases this normalization is possible it does not mean that the normalized transmission problem is independent of the initial condition.In this case one of the problem's parameter depends on the initial condition.

Passage to the Limit in the Generalized Dispersion Equation
In this section we assume that ε 2x ε 2z ε 2 and b a.Now consider the passage to the limit as a → 0. The value a 0 corresponds to the case of a linear medium in the layer.Two cases are possible: a ε 2 > 0, b ε 2 < 0 metamaterial case .Let us examine case a .The DE for a linear case is well known 10 and has the form Using passage to the limit as a → 0 we obtain the function f 1 from the function f.We seek bounded solutions X x and Z x .This implies that the denominator of the function f 1 cannot vanish.What is more, the function f as a → 0 tends to the function f 1 uniformly on x ∈ 0, h .It is possible to pass to the limit under integral sign as a → 0 in 7.14 using results of classical analysis where η 0 , η h are defined by formulas 6.4 .
The integrals in 8.3 are calculated analytically.Calculating these integrals we obtain Expression 8.4 can be easily transformed into expression 8.1 .
Let us examine b case.We have ε 2 < 0 metamaterial and the DE for the linear case has the form 31 where γ 2 − ε 1 > 0, γ 2 − ε 2 > 0, and γ 2 − ε 3 > 0.  In the same way as above, passing to the limit in the function f as a → 0 we obtain Passing to the limit in equation 7.14 as a → 0 and integrating the function f 2 , after simple calculations, we obtain formula 8.5 The results here show that it is possible to pass to the limit as a → 0. DE in 7.14 for the nonlinear case turns into 8.1 or 8.5 for the linear case as a → 0.

Numerical Results
The way of solution to the DE in 7.14 is the following: we choose the segment on γ then cut this segment into p pieces with nods γ i , i 0, 1, . . ., p. Then for each γ i we can calculate all necessary values in 7.14 .The integral in 7.14 can be calculated using any method of numerical integration.In order to calculate the value of the integrand at the point γ i it is necessary to use first integral 6.3 .For each value γ i we calculate value h i , so we obtain the grid {γ i , h i }, i 0, 1, . . ., p. Choosing reasonably dense grid on γ we can plot dependence γ on h, as it is done below.
Dispersion curves DC calculated from 7.14 , 8.1 , and 8.5 are shown in Figures 2  and 4. Eigenmodes for eigenvalues indicated in Figures 3 and 5 are shown in Figures 3 and 5.
As it is known and it is shown in Figure 2, the line γ 3 is an asymptote for DCs in the linear case.It should be noticed that in the linear case there are no DCs in the region γ 2 ≥ ε 2 .It can be proved that function h ≡ h γ defined from equation 7.14 is continuous at the neighborhood γ 2 ε 2 when a / 0 see Figure 2 .This is the important distinction between linear and nonlinear cases.Further, it can be proved that function h ≡ h γ defined from equation 7.14 when a / 0 has the following property: lim γ 2 → ∞ h γ 0. 9.1

Advances in Mathematical Physics
In Figure 2 for h 3.206 in the case of a linear layer there are 3 eigenvalues black dots where the line h 3.206 intersects with DCs .These eigenvalues correspond to 3 eigenmodes.In the case of a nonlinear layer in Figure 2 are shown 5 eigenvalues uncolored dots .These eigenvalues correspond to 5 eigenmodes.Taking into account the last paragraph's statement it is clear that in this case there is infinite number of eigenvalues.
It is easy to see in Figure 6 that the less nonlinearity coefficient a is the more stretched DCs in the nonlinear case.The maximum points of the curves h γ in Figure 6 they are marked by asterisks move to the right.The parts of the DCs that locate below the maximum points asymptotically tend to the DCs for the linear case as a → 0.
It should be noticed that in the case of Kerr nonlinearity in a layer and TE waves there are strong constraints on the value a depending on the value ε 2 for details see 16 .It is natural to suppose that there are constraints of the kind in the case under consideration.
As far as we know experiments to observe the new nonlinear eigenmodes were not carried out.So the question if the modes corresponding to the new eigenvalues exist in an experiment stays open!We should like to emphasize that it is interesting to observe purely nonlinear modes that do not arise in a linear limiting case.If to see Figure 2 then points b and c correspond to these purely nonlinear modes.
and • T denotes the operation of transposition, and each component of the fields is a function of three spatial variables.Electromagnetic field E, H satisfies the Maxwell equations rot H −iωεE, rot E iωμH, 2.3

Figure 2 :
Figure 2: Plot of γ h .The first few dispersion curves are shown.Solid curves for the nonlinear case solutions of 7.14 ; dashed curves for the linear case solutions of 8.1 .The following parameters are used for both cases: ε 1 1.44, ε 2 9, ε 3 1, and for the nonlinear case a 0.1, and E h z 1. Dashed lines are described by formulas: h * 3.206 thickness of the layer , γ 1.2 lower bound for γ , and γ 3 upper bound for γ in the case of linear medium in the layer .

Figure 3 :
Figure 3: Eigenfunctions fields for the nonlinear problem are shown.Solid curves for X; dashed curves for Z.The same parameters as in Figure 2 are used.For a , γ 2.994; for b , γ 3.892: for c , γ 8.657, and h 3.206 is used for all three cases.The eigenvalues are marked in Figure 2.

Figure 4 :
Figure 4: Plot of γ h .The first few dispersion curves are shown.Solid curves for the nonlinear case solutions of 7.14 ; dashed curve for the linear case solutions of 8.5 .The following parameters are used for both cases: ε 1 1, ε 2 −1.5, ε 3 1, and for the nonlinear case a 5.2 and E h z 1. Dashed lines are described by formulas: h * 1.71 thickness of the layer , γ 1 lower bound for γ and γ 3 upper bound for γ in the case of linear medium in the layer .

Figure 5 :
Figure 5: Eigenfunctions fields for the nonlinear problem are shown.Solid curves for X; dashed curves for Z.The same parameters as in Figure 4 are used.For a , γ 2.620, and a 0 see 8.5 ; for b , γ 1.565; for c , γ 3.481, and h 1.71 is used for all three cases.The eigenvalues are marked in Figure 4.

Figure 6 : 1 .
Figure 6: Plot of γ h for the different values of a: 1 -a 1; 2 -a 0.1; 3 -a 0.01; 4 -a 0.001; 5a 0.0001; 6 -a 0 linear case .The following parameters are used for both cases: ε 1 4, ε 2 9, ε 3 1, and for the nonlinear case E h z 1. Dashed lines are described by formulas: γ 2 lower bound for γ , γ 3 upper bound for γ in the case of linear medium in the layer .Curves solid 1-5 are solutions of 7.14 , and curve 6 dashed is solutions of 8.1 .