Effect of small anisotropy and absorption on metamaterial applications : “ non-ideal ” features of propagation and tunneling of electromagnetic waves

For cases of isotropic, uniaxial, and biaxial electromagnetic metamaterials (MM), a comparative analysis of the effect of small deviations of local material parameters from “ideal” values on the realization of MM applications (“zero” media, the Veselago-Pendry superlens) has been carried out. On the basis of the detailed investigation of the solutions of dispersion equations, it is established that even a very small dielectric and (or) magnetic anisotropy of a general form is the universal “non-ideal” factor determining (to a much greater extent than small losses) the operability of thoseMMapplications where the wavemisphasing in the effectivemedium is undesirable. The characteristics of wave attenuation in the absorbing isotropic and weakly anisotropic MM are mainly comparable for the applications. Limitations of the traditional approaches using the second-order curves (or surfaces) for analytic modeling of the absorbing MM dispersion equations are shown.


Introduction
The main ideas and experimental realization of the composite electromagnetic materials those macroscopic properties are defined by the elements of their inner microstructure (metamaterials (MM)) [1][2][3] have led to the rapid development of MM theory and applications during this century (see, e.g., recent reviews [4][5][6][7][8][9][10] and numerous references in these works).Among the promising and actively discussed MM applications, there are development of the "superlens" (the Veselago-Pendry lens, SL) with resolution higher than the diffraction limit, usage of MM with the unique optical properties including negative, near to zero or extremely large electromagnetic material parameters, creation of the various technologies for controlling the visibility of micro-and macroscopic objects, integral optics, information technologies, and many others [4][5][6][7][8][9][10][11][12].
The majority of the realized electromagnetic MM are characterized by optical absorption or losses, and composite MM are often effectively anisotropic (to a greater or lesser extent) in their dielectric and (or) magnetic properties.The production of isotropic magnetic MM in the optical range is a difficult technological problem [13][14][15].On the other hand, it is well known that namely a "nonidealness" of the effective electromagnetic characteristics according to various criteria is one of the main factors limiting the performance capabilities of MM applications.So, optical absorption worsens considerably the resolution of isotropic SL [3,[16][17][18], leads to wave misphasing and deterioration of the image [19].As shown in paper [20], even a very small uniaxial optical anisotropy of a general form is a more important problem for SL than absorption.It is caused by the fact that the misphasing of evanescent waves (that are responsible for "superresolution") even in the weakly anisotropic media increases rapidly with increasing the transversal component of wave vector.The opposite situation takes place for absorbing isotropic media: the evanescent wave phase shift decreases with increasing this wave vector component [19,20].Undesirable effects of the phase shift in wave propagation limit also the use of "zero" media (epsilon-mu-near-zero (EMNZ), epsilon-near-zero (ENZ)) that are promising for many integral optics applications [9,[21][22][23][24].
Moreover, the account of interconnected effects of absorption (losses) and anisotropy on propagation and attenuation of electromagnetic waves is also important for the MM applications.Such investigations were performed in many works only for some specific cases and concrete MM realizations (see, e.g., [25][26][27][28][29][30]).However, it is interesting to carry out the general comparative analysis of the "non-ideal" (from viewpoint of MM applications) features of electromagnetic wave propagation and tunneling for various absorbing MM: isotropic, uniaxial (at arbitrary orientations of the optical axis), and biaxial (at least, for the biaxial MM with most symmetric electromagnetic properties).Electromagnetic waves in such MM are inhomogeneous, and, in the general case, they are not reduced to TM and TE waves in the presence of effective anisotropy [11,20].To reveal the general patterns of the phase and amplitude effects of the radiation interaction with the weakly "non-ideal" MM, it is useful to investigate, at first, the features of solutions of the corresponding dispersion equations.
In the work, the exact solutions of the dispersion equations (that describe the complex wave vector components) are investigated analytically and numerically for a number of "non-ideal" effectively isotropic and anisotropic MM with local electromagnetic properties.The analysis demonstrates new features of proper waves (PW), backward waves, wave propagation and tunneling in MM, and the effect of absorption and anisotropy on the MM applications.The results are obtained using the invariant (independent of the coordinate system choice) calculation methods and the refraction vector formalism [31] taking into account absorption of inhomogeneous PW in the media.
The rest of the paper is organized as follows.In Section 2, the general problem statement and main assumptions used in the work are briefly considered.Features of the exact solutions of the dispersion equations for isotropic and anisotropic (uniaxial, biaxial) absorbing MM are investigated in Sections (3)- (5).By way of examples of the important MM applications (SL, EMNZ, ENZ), the general peculiarities of PW misphasing and damping in such materials are ascertained.First of all, these peculiarities are caused by very small deviations of MM effective electromagnetic parameters from the required "ideal" values.Each of Sections (3)-( 5) includes an analytical consideration of the corresponding type of MM and a detailed graphical analysis of the effects of small deviations of material parameters, absorption, and anisotropy on the MM applications.The conclusion (Sect.6) includes a brief discussion of the obtained results and summarizes the paper.The additional material excluded from the main text of the work is contained in Appendix.

Main assumptions
Let us consider interaction of plane monochromatic waves with effectively homogeneous media characterized by complex tensors of dielectric permittivity, e, and magnetic permeability, m.All or some components of tensors e, m can have the negative real parts, all the imaginary parts of the components are positive (absorbing linear MM).Phase multiplier exp[ik 0 (mr) À ivt] is used in expressions for the fields where m = (m 1 , m 2 , m 3 ) = k/k 0 is a complex refraction vector [31], k is a wave vector (vectors are denoted with boldface symbols), k 0 = v/c is a wave number for vacuum, i 2 = À 1, r = (x, y, z) is a radius-vector.Here and below (uv) and [u Â v] denote scalar and vector product of some (both real and complex) vectors u, v.The denotations a 0 = Re(a), a 00 = Im(a) are used for the real and imaginary parts of scalar or vector quantities.The invariant calculation methods [31] are applied partially in the present work.At first, the invariant relations are considered, and then the results are investigated in the chosen coordinate system where m = (m 1 , 0, m 3 ), that is, plane XZ is assumed to be the incidence plane.Meanwhile, the transversal component of refraction vector is supposed to be real (m 1 does not change at different boundaries in solutions of the corresponding boundary problems) and the longitudinal component m 3 is complex.The usage of the chosen coordinate system is caused by the necessity of obtaining the resulting expressions that include explicitly the components of tensors e, m, anisotropy and absorption parameters for the following analysis of the typical for MM cases of negative or small real parts of the tensors e, m components.
The validity of material equations D = eE, B = mH for the effective medium and Maxwell's equations for monochromatic waves [31] where e À1 , m À1 are inverse to e, m tensors, is also assumed.Equations ( 1) and the following below ones are true for both homogeneous and inhomogeneous waves (when the vectors corresponding to the real, m 0 , and imaginary, m 00 , parts of vector m can be non-parallel).Quantities m 0 and m 00 describe refraction and damping of the electromagnetic wave propagating in an absorbing medium, respectively.We will use the accurate relations considered below for the graphical analysis of dispersion dependences: equations ( 2) and (3) for isotropic, equations ( 13) and ( 14) for uniaxial, and equation ( 24) for biaxial MM.Let us investigate functions m 3 0 (m 1 ), m 3 00 (m 1 ), m 3 0 (m 3 00 ) only for the transmitted waves (m 3 00 > 0) at small deviations of effective dielectric and magnetic characteristics from the "ideal" values corresponding to the realization of two important MM applications: EMNZ (ENZ) and SL.We will consider the effects of both propagation (when |m 1 | values do not exceed several units) and tunneling of the electromagnetic waves (up to the values |m 1 | 50) in the investigated MM.Correspondingly, we will use two scales of the parameter m 1 changes in the figures below.The scales will be designated for convenience as "small" (when the graphs specify features of propagation of the waves) and "large" (when evanescent waves are essentially considered).

Isotropic media
Using the invariant dispersion equation for isotropic magnetic media, m 2 = em, one can obtain the following relations for the real and imaginary parts of the longitudinal component of refraction vector in the chosen coordinate system As usual, the value of m 1 is assumed to be given and the sign of quantity m 3 00 is chosen according to the causality condition: m 3 00 > 0 for the wave propagating into MM in the positive direction of axis Z with the amplitude decrease determined by parameter exp(À m 3 00 z) (further this wave is named "transmitted").The sign of quantity m 3 0 is determined with account of the previously chosen value of m 3 00 and the correspondence of the obtained complex value of m 3 to the dispersion equation.Obviously, parameter m 3 0 z determines the phase shift of the propagating wave.The similar order of the sign choice of the m 3 00 , m 3 0 values is used for uniaxial (Eqs.(12)(13)(14)) and biaxial (Eqs.( 22) and ( 24)) media.
According to equations ( 2) and ( 3), even functions m 3 0 (m 1 ) and m 3 00 (m 1 ) are near-hyperbolic and nearlinear, correspondingly, and we have the limiting relations |m 3 0 (m 1 ) (at (em) 00 ≠ 0) and m 3 00 (m 1 ) have the single extremum at m 1 = 0 (the presence of either maximum or minimum depends on the sign choice for these functions).One can simply show that the intersection of the graph of function m 3 0 (m 1 ) with the abscissa m 3 0 = 0 is impossible.That means the impossibility of transition of the forward wave (m 3 0 > 0) into backward one (m 3 0 < 0) and vice versa in the isotropic medium with changing parameter m 1 .We also have the relation 2m 3 0 m 3 00 = (em) 00 from the dispersion equation m 2 = em, so dependences m 3 0 (m 3 00 ), m 3 00 (m 3 0 ) are hyperbolic (at (em) 00 ≠ 0).A decrease (increase) of the wave attenuation in the isotropic absorbing medium thus leads to an increase (decrease) of the phase shift.
Equations ( 2) and (3) can be transformed to the form of dispersion equations in coordinates (m 3 0 , m 1 ), (m 3 00 , m 1 ): So, even in the simplest case of the isotropic absorbing medium, equations (4) and ( 5) are not described by the usually used second-order curves (or the second-order surfaces in the general case).Equations ( 4) and ( 5) are the fourth-order equations.
Let us consider some particular cases following from equations ( 2) and (3).The case m 3 0 = 0 (the purely evanescent wave) is realized for all values of m 1 only under the condition (em) 00 = e 0 m 00 + e 00 m 0 = 0 that can be satisfied for absorbing isotropic media (e 00 , m 00 > 0) at e 0 m 0 < 0.
In this case, we have According to equation (2), realization of the condition e 0 m 0 < 0 enables to decrease the phase shift (to 0 at e 0 m 00 + e 00 m 0 = 0) in wave propagation.It is of interest for the applications where such phase shift is undesirable (SL, "zero" MM).As follows from equations ( 2) and (3), both equality m 3 00 = 0 and the case of strictly linear dependence m 3 00 (m 1 ) (at m 1 0 or m 1 ≥ 0) are not realized in the isotropic absorbing media (for both "conventional" media and MM).
For the analysis of a "non-idealness" of the effective dielectric and magnetic properties of isotropic MM, let us assume that quantities e, m are of the form: where a, b are characteristic ("ideal") real values of e, m for some effects (particularly, a = b = 0 for EMNZ, a = b = À 1 for SL MM), De, Dm are small complex additions to a, b (|De|, |Dm| << 1) characterizing the "non-idealness" of MM properties.Deviations De, Dm from the ideal conditions can be caused by losses, MM fabrication methods and other factors [13][14][15].The numerical analysis of equations ( 2) and ( 3) shows that under the conditions, e.g., |De|, |Dm| < 0.1 even the linear Taylor expansion of quantities m 3 0 , m 3 00 in terms of parameters De, Dm ensures the high accuracy.One can obtain in this case Equations ( 7) describe essentially evanescent wave tunneling in MM ðm 2 1 > abÞ.According to equations (7), quantities m 3 0 , m 3 00 are determined by parameters a, b, m 1 and only the imaginary (for m 3 0 ) or real (for m 3 00 ) parts of deviations De, Dm.So, for evanescent waves in the isotropic absorbing medium, the phase shift is determined essentially by the absorption parameters (De 00 , Dm 00 > 0) and the wave attenuation is determined by deviations De 0 , Dm 0 .
Let us consider the form of equations (7) for some important MM applications.We have in the case of ideal EMNZ MM (a = b = 0): that is, the phase shift is absent and the wave attenuation parameter is equal in magnitude to the value of m 1 .
For ideal ENZ MM ða ¼ 0; b ¼ 1Þ where m 1 ≠ 0, we have the strictly hyperbolic dependence m 3 0 (m 1 ) and the difference of the linear and hyperbolic dependences for function m 3 00 (m 1 ).Meanwhile, parameter Dm has no effect on quantity m 3 for EMNZ and ENZ MM in the first approximation.
In the case of ideal SL MM (a = b = À 1), we have: where |m 1 | > 1, the dielectric and magnetic absorption effects amplify the wave phase shift simultaneously.It is possible to decrease the wave damping at De 0 + Dm 0 < 0 in comparison with the more "ideal" case De 0 = Dm 0 = 0.
For nonmagnetic SL (NMSL) ða ¼ À1; b ¼ 1Þ the full phase shift compensation at De 00 = Dm 00 and the wave damping decrease at De 0 > Dm 0 are possible.
In spite of the approximate character, equations ( 7)-( 11) (along with the accurate relations ( 2) and ( 3)) can be applied for optimization of the devices using EMNZ (ENZ) and SL (NMSL) ММ.That is confirmed by the results of the numerical and graphical analysis executed below.
Subject to parity of functions m 3 0 (m 1 ) and m 3 00 (m 1 ) for the isotropic media (Eqs.( 2) and ( 3)), their graphs are symmetric with respect to the ordinate m 1 = 0. So, it is enough to consider the effect of only positive values of parameter m 1 on quantities m 3 0 , m 3 00 .Figures 1 and 2 illustrate the effect of small values of quantities e 0 , m 0 (the top panels) and e 00 , m 00 (the bottom panels) on dependences m 3 0 (m 1 ), m 3 00 (m 1 ) for the weakly absorbing EMNZ and ENZ MM.At the large scale, dependences m 3 0 (m 1 ) are near-hyperbolic (Fig. 1 a, c, e and g), and the phase shift decreases quickly with the parameter m 1 growth for both EMNZ and ENZ MM.At the small scale (Fig. 1b, d, f, h), the dependences are different from hyperbolic and have the pronounced maximum at m 1 = 0 (the normal incidence of the wave on MM when the boundary problem is solved).For EMNZ MM, values of the maximum at variation of the parametrs e 0 , m 0 (Fig. 1b) are greater (approximately by a factor of 10 2 ) than ones at the comparable variation of quantities e 00 , m 00 (Fig. 1d).For ENZ MM, maxima of function m 3 0 (m 1 ) have the close values with changing parameters e 0 and e 00 (Fig. 1f and h).
The similar variations of quantities e 0 , m 0 and e 00 , m 00 affect dependences m 3 00 (m 1 ) (characterising the wave amplitude attenuation) considerably more weakly (Fig. 2).At the large scale, the graphs in Figure 2a are practically coincident and correpond to all the parameters sets given in Figure 2 (for both EMNZ and ENZ, Fig. 2b-e).Dependences m 3 00 (m 1 ) for EMNZ and ENZ MM are slightly different at small values of m 1 .In particular, the variation of quantities e, m affects parameter m 3 00 for ENZ ММ at the same values of m 1 more strongly (Fig. 2b and c in comparison with Fig. 2d and e).The growth of the values of e 0 , m 0 (e ", m 00 ) leads to a decrease (increase) of parameter m 3 00 , correspondingly, at fixed values of m 1 (Fig. 2b, d, c and  e).
The considered features are characteristic and take place at various values of the parameters, e.g., under the condition e 0 m 0 < 0 (Fig. 2b, the top curve).At e 0 ¼ m 0 ¼ À0:001; À0:01; À0:1 the graphs of function m 3 00 (m 1 ) are practically invariable in comparison with ones given in Figure 2a, and graphs m 3 0 (m 1 ) became "mirror-image" with respect to the abscissa axis (that is, the transmitted wave becames backward wave).
The graphical analysis also shows that at the relation (em) 00 ≈ 0 (the values je 0 j ¼ jm 0 j ¼ 0:001; 0:01; 0:1 , e 0 m 0 < 0 are taken) the condition m 3 0 ≈ 0 is satisfied with the very high accuracy up to distinct values of absorption parameters of the order of e 00 , m 00 ≈ 0.1.In this case, the graphs of function m 3 00 (m 1 ) correspond to ones in Figure 2a.
The results of the graphical analysis which is similar to Figures 1 and 2 data (but for the conditions of isotropic SL and NMSL realization) are given in Figures 3 and 4. According to data in Figure 3 a-d, the transmitted wave in the case of SL is backward (m 3 0 < 0), and functions m 3 0 (m 1 ) at varying parameters e 0 , m 0 and e 00 , m 00 are characterized by the close values of minima m 3 0 (0) in the diapason (À0.8, À1) (Fig. 3b and d).In the NMSL case, the wave is forward (m 3 0 > 0, Fig. 3e-h).At varying parameters e 0 and e 00 , functions m 3 0 (m 1 ) have maxima m 3 0 (0) which values differ in approximately two orders (Fig. 3f and h).For both cases (SL, NMSL) the phase shift decreases quickly according to the near-hyperbolic law with the growth of m 1 up to the large values (Fig. 3a, c, e  and g).In the case of NMSL, the smaller values of |m 3 0 | in comparison with SL at the same m 1 values (Fig. 3a, c, and e, g) are realized.In Figure 3e and f the graphs are given for the cases e 00 = 0.001, m 00 = 0 ("magnetic" absorption is absent) and e 00 = m 00 = 0.001 (the coincident parameters of "dielectric" and "magnetic" absorption that corresponds to the condition De 00 = Dm 00 in the expression for m 3 0 (Eq.( 11)).The significant decrease of the phase shift takes place in the last case, as here the condition (em) 00 ≈ 0 is satisfied with the high accuracy.
The graphs in Figure 4 for function m 3 00 (m 1 ) at the SL, NMSL regimes are qualitatively similar to the corresponding graphs in Figure 2 at small values of m 1 and are practically coincident at the large scale (Fig. 2a and Fig. 4a).However, we have the relations of the parameters describing the given effects at small coincident values of m 1 : ðm 3 00 Þ SL << ðm 3 00 Þ EMNZ (Fig. 4b, c and Fig. 2b, c), ðm 3 00 Þ NMSL > ðm 3 00 Þ ENZ (Fig. 4d, e and Fig. 2d, e).At small values of m 1 , varying quantities e 0 , m 0 affects dependences m 3 00 (m 1 ) for NMSL more strongly in comparison with SL (Fig. 4b and d).Changing e 00 , m 00 values does not affect practically the graphs of function m 3 00 (m 1 ) for NMSL (Fig. 4e) but becomes significantly apparent for SL (Fig. 4c).
The data in Figures 1-4 are represented in Figure 5 in axes (m 3 0 ,m 3 00 ) where the left (Fig. 5a and c) and right (Fig. 5b and d) panels characterize the effect of varying parameters e 0 , m 0 and e 00 , m 00 .For the graphs in Figure 5, the removal of the curves from the coordinate origin (m 3 0 = m 3 00 = 0) corresponds to amplification of the "non-idealness" of the MM effective characteristics.According to data in Figure 5, dependences m 3 0 (m 3 00 ) are hyperbolic for the isotropic MM, and the absorption parameters effect on m 3 0 (m 3 00 ) graphs determinatively (especially for ENZ, SL, NMSL MM, Fig. 5b and d).The same changes of parameters e 0 , m 0 and e 00 , m 00 lead to the practically coincident m 3 0 (m 3 00 ) graphs for EMNZ MM and strongly different ones for ENZ MM (Fig. 5a and b).The variation of the values of quantities e 0 , m 0 effects on dependences m 3 0 (m 3 00 ) for ENZ, SL, NMSL MM weakly (Fig. 5a and c).The comparison of the data in Figures 1-5 shows that the phase shift decrease is accompanied by the increase of the wave damping with the growth of m 1 (in the general case, |m 1 |) values.So, for the isotropic absorbing EMNZ (ENZ) and SL (NMSL) MM the phase and amplitude effects actuate mutually antithetically with the growth of transversal wave number k 1 = k o m 1 (that is, at the transition to "more evanescent" waves): the phase shift decrease is accompanied by the increase of wave attenuation, and vice versa.The increase of the effective medium absorption leads, as a rule, to the phase shift increase (with the exception of the case (em) 00 ≈ 0 where the "dielectric" and "magnetic" absorption effects compensate mutually their impact on the phase shift).

Uniaxial media
In the general case of the uniaxial magnetic medium, nondiagonal tensors e, m have the form [31] corresponding to TM and TE waves only in the particular case of the optical axis orientation in the incidence plane.
According to equations (12), the waves in such medium are not strictly propagating or evanescent (for the transmitted waves always m 3 00 > 0 and, in the general case, m 3 0 ≠ 0 even at the incidence on a MM layer of the purely evanescent wave for which m 3 0 = 0, m 3 00 > 0).Further we use the chosen coordinate system where m = (m 1 , 0, m 3 ), c = (cos(a), cos(b), cos(g)) and consider the first equation from equations (12) corresponding to TM waves in the particular case cos(b) = 0.The results presented below are analogous to the second equation from equations (12) (corresponding to TE waves at cos(b) = 0) at the substitutions m o ↔ e o , e e ↔ m e .Consequently, we obtain the quadratic equation in complex parameter m 3 with complex coefficients:   The extraction of the real and imaginary parts in equation ( 13) leads to a system of two complete fourthorder equations of the form f(m 3 0 , m 1 ) = 0, g(m 3 00 , m 1 ) = 0 both in quantities m 3 0 , m 3 00 and parameter m 1 .Meanwhile, the solubility conditions for these equations (availability of the real solutions for m 3 0 , m 3 00 having the physical sense) are identical.The real solutions of these equations can be obtained using the standard analytical method on the basis of the cubic resolvent: where the denotations are used (p = e e /e o , q = e e m o ): According to equations (14), quantities m 3 0 , m 3 00 are determined by the following parameters: the optical axis direction, quantity m 1 , two complex parameters dependent on the medium material properties, e e /e o and e e m o .Thus, the waves with the same refraction vector can be excited in the different uniaxial media (characterized by the different parameters e o , e e , m o ) at the identical aforesaid parameters.It is also followed from equations (14), that functions m 3 0 (m 1 ), m 3 00 (m 1 ) are near-linear, especially with the growth of the |m 1 | values corresponding to evanescent waves.So, dependences m 3 00 (m 1 ) for the isotropic and uniaxial media are qualitatively similar, and dependences m 3 0 (m 1 ) are different in essence.With increasing |m 1 | values the phase shift of the wave in the isotropic medium decreases quickly according to a near-hyperbolic law, but in the uniaxial medium it increases practically linearly.One can simply show that, in the general case, dependences m 3 0 (m 1 ), m 3 00 (m 1 ) take place in any (non-magnetic, magnetic, MM) uniaxial absorbing medium.Meanwhile, in contrast to the isotropic medium, the accurate (and even approximate, to the first order approximation on the absorption and anisotropy parameters) realization of the condition m 3 0 = 0 at all the parameter m 1 values is impossible.So, it is impossible to obtain the strictly plane dispersion surface required for many MM applications.According to equations ( 14), the intersection of the function m 3 0 (m 1 ) graph with the abscissa m 3 0 = 0 is possible.That is one more difference from the isotropic medium case, and it means the possibility of transition of the forward wave into backward one (and vice versa) in the uniaxial medium at changing parameter m 1 [32].
The important particular case of equations ( 12) is determined by the condition [m Â c] 2 = 0 corresponding to realization of the peculiar PW, the Fedorov À Petrov waves [33] considered for uniaxial MM in [20].The other interesting cases are realizations of the conditions e o = 0 or e e = 0 in equation ( 13) corresponding to the cases of "partially absorbing" ENZ MM.
At e o = 0, any values of m o and e e ≠ 0, we have from equation ( 13):  that is, an arbitrarily large MM absorption determined by parameters m o 00 e e 00 does not influence quantities m 3 0 , m 3 00 .It is interesting that, according to the relations for PW polarization in the uniaxial medium [20], vector E is "parallel" to the optical axis in this case (it means that [E Â c] = 0 with account of complexity of E).
At e e = 0, any values of m o and e o ≠ 0, one can obtain from equation ( 13): The feature of equations ( 15) and ( 16) is their linearity with respect to quantities m 1 , m 3 .Meanwhile, they determine the only possible value of m 3 0 (and in the case of Eq. ( 15) also the only value of m 3 , and, therefore, the only PW polarization state) at a given m 1 for both the transmitted (m 3 00 > 0) and reflected (m 3 00 < 0) waves characterized by equation (13).One of these waves is forward (m 3 0 > 0) and the other is backward (m 3 0 < 0) subject to the optical axis orientation.Note, that the similar case of the linear dispersion equation is impossible for the isotropic medium.It is also seen that explicit dependences of quantities m 3 0 , m 3 00 on MM material parameters are vanished in these particular cases, and quantity m 3 is determined only by the optical axis direction and m 1 value.Under the additional conditions to the optical axis orientations cos(a) = 0, cos(g) ≠ 0 (the case of Eq. ( 15)) or cos(g) = cos(b) = 0 (the case of Eq. ( 16)) such MM are "ideal" ENZ ММ: m 3 0 = m 3 00 = 0 at any values of m 1 .Thus, for the considered cases (Eqs.( 15) and ( 16)) the dispersion equations are determined only by the optical axis direction that, obviously, depends on tensor e, m components.
Let us use equation ( 13) for obtaining the approximate relations that are similar to equations ( 7) and characterize dependences m 3 0 (m 1 ), m 3 00 (m 1 ) in the first order approximation on the "non-ideal" deviations of the absorption and anisotropy parameters.As in equations ( 6), we assume that where m 2 1 > ab, Q ¼ absin 2 ðgÞ À m 2 1 ½cos 2 ðaÞ À cos 2 ðgÞ, and only the upper or only lower sign is chosen in each relation.Expressions (18) and ( 19) generalize the relations for SL MM (a = b = À 1) [20] in the case of arbitrary real values of a ≠ 0, b.
At a = 0 (EMNZ or ENZ ММ) the following expression is a highly accurate approximation for parameter m 3

See equation (20) below:
According to equation (20), parameter m 3 is directly proportional to m 1 and independent of the effective medium magnetic properties.
In the general case of anisotropic media, functions m 3 0 (m 1 ) and m 3 00 (m 1 ) are not appeared to be even (Eqs.( 14) and ( 24)).So, let us consider the effect of not only positive but also negative values of parameter m 1 in the diapason |m 1 | 50 on quantities m 3 0 , m 3 00 for the transmitted waves.For all the graphs in Figures 6-8 the values of the angles a = p/4, b = g = p/3 determining the optical axis direction in the chosen coordinate system are used.The features considered below take place also for the different orientations of the optical axis.
The graph of linear function m 3 0 (m 1 ) for the case e o = 0 (Eq.(15))isgiveninFigure6a.In thiscase, thewaveamplitude attenuation is absent: m 3 00 = 0.The very small change of parameter e o (from 0 to 10 À11 i, Fig. 6b) leads to the appearance of the strong and near-linear dependences m 3 00 (m 1 ) that are very "sensitive" to very small changes of quantity e e 0 .The ) leads to near to piecewise-linear graphs of functions m 3 0 (m 1 ) and m 3 00 (m 1 ) that are asymmetric with respect to the positive and negative values of m 1 (Fig. 6c and d).With that, it is possible to change the transmitted wave type (forward, m 3 0 > 0, backward, m 3 0 < 0) with varying parameter m 1 (Fig. 6a and Fig. 6c, curve 1).Also the different character of the changes of dependences m 3 00 (m 1 ) at m 1 < 0 and m 1 > 0 (Fig. 6d) is possible.The graphs in Figure 6 do not practically change also at the small scale of values of m 1 and under the condition m o = 1, that is, forthe case of weakly anisotropic uniaxial ENZ MM.
According to Figure 6 data, the very small changes of parameters e o , e e (of the order of |e o |, |e e | ≈ 10 À5 À 10 À4 and less) lead to the significant changes of quantities m 3 0 (Fig. 6c) and m 3 00 (Fig. 6b and d and the transition from Fig. 6a to b) for evanescent waves (at |m 1 | >> 0).Such features of the considered near to "zero" MM can be of interest for controlling these waves characteristics and sensor applications.
The data in Figures 7 and 8 characterize the effect of parameters of small birefringence (Fig. 7), dichroism (Fig. 8 a and d), birefringence and dichroism (Fig. 8b-d) on the conditions of the SL MM realization (the analogous features take place for NMSL MM).At the large scale of m 1 values, dependences m 3 0 (m 1 ), m 3 00 (m 1 ) are near-linear (Fig. 7a, c, and Fig. 8a-c), and the graphs m 3 0 (m 1 ) are qualitatively similar to the graphs for EMNZ ММ (Fig. 6c).However, the phase shifts are appeared to be in 3-4 order smaller in comparison with the EMNZ ММ regime with changing values m 1 in the case of SL (Fig. 7a and Fig. 8a-c in comparison with Fig. 6a and c).At small m 1 values, dependences m 3 0 (m 1 ) are nonlinear and symmetric with respect to axis m 1 = 0 (Fig. 7b, d, and Fig. 8d).At the large scale, the function m 3 00 (m 1 ) graphs in Figure 7c coincide and correspond to all the parameters sets given in Figures 7  and 8.The comparison of the data in Figure 8a-c shows that asymmetry of the function m 3 0 (m 1 ) graphs with respect to the positive and negative m 1 values is determined mainly by the effective medium birefringence.With that, the phase shifts due to birefringence are appeared to be in several times greater than ones due to dichroism (at the same values of birefringence, |e e 0 À e o 0 |, and dichroism, |e e 00 À e o 00 |, parameters, curves 4-9 and 1-3 in Fig. 8).
Figure 9 generalizes Figures 6-8 data in the form of dependences m 3 0 (m 3 00 ).According to Figure 9 graphs at the large scale of m 1 values (at large m 3 00 values), dependences m 3 0 (m 3 00 ) are near-linear for the weakly uniaxial MM.At small m 1 values (at small m 3 00 values), graphs m 3 0 (m 3 00 ) are near-linear only for EMNZ MM (for the data in Fig. 9a and  b).In the case of SL MM with small uniaxial anisotropy, dependences m 3 0 (m 3 00 ) can be close to piecewise-linear (Fig. 9d and f) and nonlinear (Fig. 9h).Also the data in Figure 9 illustrate clearly the different scale of parameter m 3 0 changes with the growth of the absorption parameters for the considered EMNZ and SL MM.

Biaxial media
The cases of biaxial absorbing media can be realized in modeling and experimental determination of effective electromagnetic MM characteristics.With that, MM dielectric and magnetic characteristics can demonstrate the different symmetry properties.Consequently, e.g., effective tensors e, m cannot be reduced to the diagonal form in the same coordinate system.Obviously, the properties of such MM become more complex in comparison with the cases of isotropic and uniaxial media.For such problems (as for cases of "conventional" anisotropic media), the invariant methods of phenomenological electrodynamics [31,34] are high-performance ones.Meanwhile, as a rule, MM are artificial materials from "meta-atoms" ordered spatially according to some rules, that is, the given materials are frequently characterized by pronounced symmetry properties.Therefore, let us consider as a biaxial MM model a rather general case of one from the most symmetric configurations such as absorbing magnetic medium of orthorhombic symmetry with the identical structure of effective tensors e, m.
In this case, non-diagonal tensors e, m are reduced to the diagonal form in the same real coordinate system and can be expressed in the following form [31]: where e q , m q are complex scalars (eigenvalues of tensors e, m), c q are unit real vectors (q = 1, 2, 3, c 2 q ¼ 1) directed along the second order symmetry axes which are mutually orthogonal in the considered media (c q are also perpendicular to the planes of the medium mirror symmetry), c q .cq are dyads.For example, for layered MM, one of the vectors c q can be perpendicular to the layers boundaries and two others can lie in the layers plane.
Substitution of expressions ( 21) into the well-known invariant equation of normals for magnetoanisotropic media characterized by arbitrary tensors e, m [31,34] leads to the dispersion equation: where the denotations are used: M q = (mc q ), F pr = m p /e p + m r /e r (q, p, r = 1, 2, 3).One can show that only in the particular cases M q = 0 the left part of equation ( 22) is representable in the form of the product of two multipliers that are quadratic with respect to refraction vector m.Then complex vector m "lies" in one from three principal planes of the medium (one from three scalar products (mc q ) is equal to 0).Only in these cases it is possible to separate PW in the medium into two groups corresponding to TM and TE waves at the transition to diagonal tensors e, m.
One more particular case of equation ( 22) is realized under the condition e = Km where K is a complex scalar (direct proportionality of tensors e and m, the onerefringent medium case [31]).Under this condition, equation (22) leads to the relations In this case, refraction vector depends on the direction of the wave propagation but it does not depend on the wave polarization (birefringence is absent) [31].
At the transition to the chosen coordinate system where m = (m 1 , 0, m 3 ), c q = (cos(a q ), cos(b q ), cos(g q )) equation ( 22) leads to the complete fourth-order equation in quantity m 3 (and also in m 1 ) where the expressions for coefficients A, B, C, D, E are given in Appendix.Extraction of the real and imaginary parts of equation ( 24) leads to the system of two sufficiently awkward equations and each of them is the complete fourth-order equation in quantities m 3 0 , m 3 00 , m 1 .In the general case, the analytical solution of the given system is seemed to be impossible.Obviously, knowing the coefficients of equation (24), one can obtain complex analytical solutions for quantity m 3 using the well-known formulas or applying the numerical methods to find four complex roots of the equation.Both these methods (analytical and numerical) of the solution of dispersion equation (24) are used for finding and verification of the values of m 3 .Then the values of m 3 0 , m 3 00 are chosen with account of execution of the causality conditions and correpondence of the obtained values pairs (m 3 0 , m 3 00 ) to equation (24).The results of the graphical analysis of equation (24) under the conditions that are close to realization of EMNZ (Figs. 10 and 13a-c) and SL (Figs. 11,12, 13d and e) MM are given in Figures 10-13.For visualization, the different  ) for the biaxial MM.The graphs at the large scale of m 1 values in Figure 10 (that is, besides of Fig. 10e) correspond qualitatively to the case of the weakly anisotropic EMNZ MM (Fig. 6b and c).In this case, even very small birefringence parameters (of the order of 10 À5 , Fig. 10a and b) lead to the phase shifts that exceed in 5-6 orders and greater the corresponding values for the isotropic MM (Fig. 1).The growth of the dielectric and magnetic anisotropy parameters (Fig. 10d and e) leads to more stronger m 3 0 (m 1 ) dependences and nonlinearity of these dependences at small m 1 values (Fig. 10e).It is significant that dependences m 3 00 (m 1 ) are appeared to be "steady" to MM small anisotropy (Fig. 10c, f, and 11d, the similar graphs take place also for the data in Fig. 12).
For the nonabsorbing SL MM (Fig. 11, the data in Fig. 11a at small m 3 0 and m 1 values are specified by the graphs in Fig. 11b and c), the presence of small birefringence of the order of 10 À5 leads to the qualitative changes of dependences m 3 0 (m 1 ) in comparison with the isotropic SL case (Fig. 3).In particular, the forward wave appears additionally to the backward wave in MM, with that, phase incursions of these waves rise with growing |m 1 | values (Fig. 11b).At the same parameters and m q 0 = 1 (the case of NMSL), the function m 3 00 (m 1 ) graphs (Fig. 11d) do not change practically.In this case, dependences m 3 0 (m 1 ) correspond to Fig. 11b (but without the features in the diapason m 1 ∈ (À 1, 1), the graphs are near-linear).The similar characteristic properties take place also at small deviations of the same order from the values e q 0 = À 1 to the larger values.
In the case of SL MM (Fig. 12a), the presence of dichroism leads to m 3 0 (m 1 ) dependences that are practically symmetric with respect to the axis m 1 = 0 and become stronger with the growth of quantities |e 2 00 À e 1 00 |, |e 3 00 À e 1 00 | (determined by parameter p).These dependences are near-linear at the values |m 1 | > 5 À 10.As for the uniaxial MM (Figs. 7 and 8), birefringence (characterized in this case by quantities |e 2 0 À e 1 0 |, |e 3 0 À e 1 0 |) determines the asymmetry of the function m 3 0 (m 1 ) graphs with respect to the parameter m 1 sign (Fig. 12b).Small dichroism and anisotropy (of the order of 10 À4 -10 À3 ) in the considered MM lead to the phase shift values corresponding to rather large absorption in the isotropic ММ (Fig. 1).At the small scale of |m 1 | values, the graphs in Figure 12a and b are practically coincident and correspond to the data in Figure 12c.The similar features take place also in the presence of small magnetic anisotropy of the same order and at the e q , m q values near to À1, 1 correspondingly (the NMSL case).
The Figures 10 and 12 data are illustrated by Figure 13 in axes (m 3 0 ,m 3 00 ).In the case of EMNZ MM, the presence of only small birefringence leads to near-linear dependences m 3 0 (m 3 00 ) for both the large (Fig. 13a) and small scale of |m 1 | values.With the growth of the dielectric and magnetic anisotropy parameters, dependences m 3 0 (m 3 00 ) became significantly stronger (Fig. 13b) and nonlinear at small values of m 1 (Fig. 13c).In the case of SL MM (Fig. 13d and  e), the mutual effect of birefringence and dichroism also leads to the nonlinear graphs for both small and large absorption corresponding to the small and large values of parameter |m 1 |.In the given case (as for the uniaxial MM), the different scale of parameter m 3 0 changes takes place (the difference is approximately in a three order of magnitude) for EMNZ (Fig. 13a and b) and SL (Fig. 13d and e) MM at the comparable variations of the anisotropy parameters.
The features of the dispersion dependences of uniaxial (Figs.6-9) and biaxial (Figs.10-13) MM are characteristic and take place also at other MM material parameters, directions of the optical axis for the uniaxial MM and orientations of vectors c q for the biaxial MM.

Conclusion
The executed analysis points to the crucial role of effective dielectric and magnetic anisotropy (in comparison with absorption) in forming the undesirable phase incursion of electromagnetic waves in the MM applications.The presence of even a very small anisotropy leads generally to near-linear (piecewise-linear) dependences of the phase shift on the transversal wave number (the phase incursion increases quickly with the growth of values |m 1 |).The characteristics of wave attenuation for the absorbing isotropic and weakly anisotropic MM comparable (the function m 3 00 (m 1 ) graphs at the large scale of m 1 values are weakly discriminated for all the considered cases excepting some peculiar cases, Fig. 6b).However, very small (and so hardly controllable for real MM) anisotropy (of the order of |De|, |Dm| ≈ 10 À5 À 10 À3 ) leads to the same phase shifts as ones caused by rather large absorption in the isotropic MM (of the order of e 00 , m 00 ≈ 0.1 and larger) in the range of "working values" of m 1 .That causes a quick violation of "ideal" conditions for the MM applications: the image deterioration and withdrawal from the "superlensing" regime for SL (NMSL), unwanted misphasing of the waves for EMNZ (ENZ) MM.So, one of the conclusions of paper [20] (where uniaxial MM and the realization of SL were considered) about the crucial influence of small anisotropy on the SL operability can be significantly generalized.Small electromagnetic anisotropy of a general form (dielectric, magnetic, uniaxial, biaxial) is the universal "non-ideal" factor determining (to a much greater extent than small losses) the operability of SL (NMSL), EMNZ (ENZ) MM and, apparently, also the other MM applications where the wave misphasing in the effective medium is undesirable.The presence of effective anisotropy (either required for applications or unwanted) is the feature of many realized MM.Therefore, the careful account of both losses and even a very small permittivity (permeability) anisotropy for various MM applications is necessary.Besides the comparative analysis of dispersion equations for MM with various anisotropic properties, a number of important particular cases (that can be of interest for optimization of the MM applications) has been investigated.
Obviously, the obtained results accounting mainly the dispersion equations can be considered only as initial estimations of the operability of "non-ideal" MM applications.For a more detailed analysis, e.g., of SL MM, one can determine the fields in the image region or transfer functions for the concrete SL devices [3,17,19,26,35] and consider the change of the images with increasing absorption and anisotropy parameters.However, the ascertained features of the wave misphasing and damping in MM are rather general ("from first principles").So, they will appear also in cases of a more detailed consideration of the concrete devices using MM.The results can be useful also in the investigations of interaction of electromagnetic fields of a more complicated structure (wave beams, waveguide modes, et al.) with MM for the cases when these fields are represented by monochromatic waves superpositions.
Moreover, one can note the limitations of the traditional approaches using the second-order curves (or surfaces) for analytical modeling of the dispersion curves (surfaces) of absorbing MM.The obtained analytical and numerical results illustrate clearly that even a very small electromag-   netic absorption and anisotropy lead to the cardinal changes of complex solutions of the second-order dispersion equations.Even for the simplest cases of the absorbing isotropic uniaxial media (including MM), dependences m 3 0 (m 1 ) and m 3 00 (m 1 ) are characterized by the fourth-order curves, and these curves (surfaces) have even higher than four orders in the case of absorbing biaxial media.Rise of the order of the dispersion dependences complicates appreciably MM investigations but leads to more realistic and useful models of MM and devices on the basis of MM.

Implications and influences
On the basis of the detailed analysis "from first principles", it is established that even a very small (hundredths of a percent and less) effective electromagnetic anisotropy is the universal "non-ideal" factor determining (to a much greater extent than small losses) the operability of the MM applications ("zero" media, "superlens").So, a careful account of both the material absorption and even small permittivity (permeability) anisotropy of the general form for various MM applications is necessary.Limitations of the traditional approaches using the second-order dispersion curves (or surfaces) for analytic modeling of the absorbing electromagnetic MM with local material parameters are shown.It is required to use the dispersion equations of the fourth and higher orders as more realistic and useful models of the considered MM.

Fig. 1 .
Fig. 1.The effect of small values of quantities e 0 , m 0 (the top panels) and e 00 , m 00 (the bottom panels) on dependences m 3 0 (m 1 ) for isotropic EMNZ (a-d) and ENZ (e-h) ММ.Figures b, d, f and h detail the corresponding graphs at small values of m 1 .
: e = e o + (e e À e o )c.c, m = m o + (m e À m o )c.cwhere e o (m o ) is double singular, e e (m e ) is non-singular complex eigenvalue of tensor e(m), c is a unit real vector in the optical axis direction, c.c is a dyadic product.This case is described by the invariant dispersion equations [20] e e ðm 2 À e o m o Þ ¼ ðe e À e o Þ½m Â c 2 ;

Fig. 2 .
Fig. 2. The effect of small values of quantities e 0 , m 0 (the top panels) and e 00 , m 00 (the bottom panels) on dependences m 3 00 (m 1 ) for isotropic EMNZ (a-c) and ENZ (d, e) ММ.Figure a corresponds to all the data given in Figures b-e.
Fig. 2. The effect of small values of quantities e 0 , m 0 (the top panels) and e 00 , m 00 (the bottom panels) on dependences m 3 00 (m 1 ) for isotropic EMNZ (a-c) and ENZ (d, e) ММ.Figure a corresponds to all the data given in Figures b-e.

Fig. 3 .
Fig. 3.The effect of small values of quantities e 0 , m 0 (the top panels) and e 00 , m 00 (the bottom panels) on dependences m 3 0 (m 1 ) for isotropic SL (a-d) and NMSL (e-h) ММ.Figures b, d, f and h detail the corresponding graphs at small values of m 1 .

Fig. 4 .
Fig. 4. The effect of small values of quantities e 0 , m 0 (the top panels) and e 00 , m 00 (the bottom panels) on dependences m 3 00 (m 1 ) for isotropic SL (a-c) and NMSL (d, e) MM.Figure (a) corresponds to all the data given in b-e.Fig. 5.The effect of small values of quantities e 0 , m 0 (the left panels) and e 00 , m 00 (the right panels) on dependences m 3 0 (m 3 00 ) for the isotropic EMNZ and ENZ (a, b), SL and NMSL (c, d) ММ.The corresponding data of Figures 1-4 are taken for the function m 3 0 (m 3 00 ) graphs.

5 .
Fig. 4. The effect of small values of quantities e 0 , m 0 (the top panels) and e 00 , m 00 (the bottom panels) on dependences m 3 00 (m 1 ) for isotropic SL (a-c) and NMSL (d, e) MM.Figure (a) corresponds to all the data given in b-e.Fig. 5.The effect of small values of quantities e 0 , m 0 (the left panels) and e 00 , m 00 (the right panels) on dependences m 3 0 (m 3 00 ) for the isotropic EMNZ and ENZ (a, b), SL and NMSL (c, d) ММ.The corresponding data of Figures 1-4 are taken for the function m 3 0 (m 3 00 ) graphs.
where a, b, are "ideal" (for some effects) real values of quantities e o , m o , small complex parameters De o , Dm o , De (|De o |, |Dm o |, |De| << 1) describe small deviations of the MM effective parameters from "ideal" conditions, and De characterizes small dielectric anisotropy.According to the numerical analysis, the condition |De o |, |Dm o |, |De| 0.1 ensures the high accuracy again.Consequently, we obtain

Fig. 6 .
Fig. 6.The effect of small birefringence (a, b) and mutual effect of birefringence and dichroism (c, d, values of the parameters for curves 1-3 are given in the table) on dependences m 3 0 (m 1 ) (the left panels) and m 3 00 (m 1 ) (the right panels) for the uniaxial EMNZ ММ.

Fig. 8 .
Fig. 8.The effect of small dichroism (a), mutual effect of birefringence and dichroism (b, c) (values of the parameters for curves 1-9 are given in the corresponding tables) on dependences m 3 0 (m 1 ) for the uniaxial SL ММ. Figure d corresponds to the cases (a-c) at small values of m 1 .

Fig. 7 .
Fig. 7.The effect of small birefringence on dependences m 3 0 (m 1 ) (a, b) and m 3 00 (m 1 ) (c, d) for the nonabsorbing uniaxial SL ММ.Figures b, d detail the corresponding graphs at small values of m 1 .

Fig. 9 .
Fig. 9. Dependences m 3 0 (m 3 00 ) for the data given in Figures 6-8.The data correspondence is: a À Figures 6a, b; b À Figures 6c, d; c, d À Figure 7a-d; e, f À Figure 8a; g, h À Figure 8b (the same line types are used for the corresponding data).Figures d, f and h detail the data in Figures c, e and g at small values of m 3 00 .
Figure 10a-c, Figures 11 and 13a illustrate the effect of small dielectric birefringence (here the very small e 00 , m 00 values are taken for extraction of the transmitted waves), Figure 10 d-f, c and Figures 12, 13b-e characterize the mutual effect of birefringence and dichroism on dependences m 3 0 (m 1 ), m 3 00 (m 1 ), m 3 0 (m 3 00

Fig. 10 .
Fig. 10.The effect of small birefringence (a-c), the mutual effect of birefringence and dichroism (d-f) on dependences m 3 0 (m 1 ) (a, b, d, e) and m 3 00 (m 1 ) (c, f) for the biaxial EMNZ ММ. Figure b details the data in Figure a at small values of m 3 0 , Figure e details the data in Figure d at small values of m 1 .

Fig. 13 .
Fig. 13.Dependences m 3 0 (m 3 00 ) for the data given in Figures 10, 12.The data correspondence is: a À Figures 10a-c; b, c À Figures 10d-f; d, e À Figures 12b, c (at m 1 < 0 and m 1 > 0).The same line types are used for the corresponding data in Figures 10, 12 and 13. Figure c details the data in Figure b at small values of m 3 00 .

Fig. 11 .
Fig. 13.Dependences m 3 0 (m 3 00 ) for the data given in Figures 10, 12.The data correspondence is: a À Figures 10a-c; b, c À Figures 10d-f; d, e À Figures 12b, c (at m 1 < 0 and m 1 > 0).The same line types are used for the corresponding data in Figures 10, 12 and 13. Figure c details the data in Figure b at small values of m 3 00 .

Fig. 12 .
Fig. 12.The effect of small dichroism (a), the mutual effect of birefringence and dichroism (b) on dependences m 3 0 (m 1 ) for the biaxial SL ММ. Figure c corresponds to the cases (a, b) at small values of m 1 .