# Penetration Effect: Exotic Behavior of a Wave in Anisotropic Media.

1. IntroductionA phenomenon of plane electromagnetic wave propagation in vacuum parallel to an interface between vacuum and an anisotropic medium has been considered in [1]. It has been shown analytically in that paper for the first time that a bulk wave can be excited within an anisotropic medium in the considered case. Also this phenomenon has been defined as "penetration effect" for the first time. It has been obtained that a real part of the Poynting vector normal component of a wave propagating in an anisotropic medium is not equal to zero when an incident wave propagates tangentially to an interface between the media. This fact gives us the possibility of talking about electromagnetic wave propagation within an anisotropic medium under an oblique angle. In addition, the authors have suggested the possible practical applications of this effect. In particular, it has been proposed to use these structures to construct frequency detectors and power dividers for microwave communication systems. However, such an important point as behavior of a reflected wave has not been considered in [1]. The hypothesis about reflected wave propagation perpendicular to an interface has been presented in [2]. but this problem has not been studied in detail. For example, in that work, a deep derivation of the expressions for the reflection and transmission coefficients has not been presented. The physical effect in [2] has not been also substantiated.

It is obvious that studying behavior of the transmitted and reflected waves when an incident wave propagates parallel to an interface between anisotropic and isotropic media has important theoretical and practical significance. Indeed this issue has not been studied theoretically despite the fact that a similar phenomenon has been described even in hydrodynamics [3]. From a practical point of view, this question is important as this phenomenon leads to additional power losses. Indeed, earlier, power losses determined by "penetration effect" have not been taken into account. Furthermore, penetration effect allows us to design new optical and microwave devices.

In this paper, we consider the ideal case of wave propagation over an infinite interface between free space and a uniaxial dielectric anisotropic medium without frequency dispersion. An orientation of the anisotropy axis is taken arbitrary in our treatment. At the same time, we do not study the problem of the possibility of exciting such a wave and the grounds for this assumption. It is important to note that this particular case cannot be described by using well-known methods [4] as the transmission coefficient would be zero if the classical theory was used. However, in accordance with the solution of Maxwell's equations, there is a transmitted wave in the anisotropic medium. Obviously, the results presented in this study require further experimental verification. However, at the same time, all expressions are obtained directly from Maxwell's equations under the classical boundary conditions for the field and wave vector components by using accurate analytical transformations only. These facts confirm the correctness of the presented results.

The behavior of both transmitting in the anisotropic medium and reflecting in vacuum waves is investigated in this paper too. It is shown analytically that a reflected wave should propagate perpendicularly to the interface to satisfying Maxwell's equations and the boundary conditions simultaneously. The expressions of reflection and transmission coefficients for a semi-infinite medium are obtained and the analysis of these coefficients is carried out.

2. Statement of the Problem

In this paper, we study the particular case of the uniaxial anisotropic dielectric medium with the scalar permeability [mu] and the permittivity dyadic in the following form:

[mathematical expression not reproducible] (1)

For example, it can be a plasma medium near the resonance region [5] when external bias field is small and therefore the off-diagonal components [[epsilon].sub.12] and [[epsilon].sub.21] in (1) can be negligible.

It is assumed that the anisotropy axis is orientated in the incidence plane under an arbitrary angle [beta] with the normal to the interface (Figure 1). Now using Euler transformation [6], we can write the electric permittivity dyadic in the xyz-system as follows:

[mathematical expression not reproducible] (2)

The geometry of the problem is presented in Figure 1. Here X' y' z' is the coordinate system connected with the permittivity dyadic (1) (z' is the anisotropy axis); xyz is the coordinate system connected with the normal to the interface, where z is a normal to the interface.

In the considered case, a plane harmonic wave with the wavenumber [k.sub.inc] = [k.sub.y] = [omega] [square root of ([[mu].sub.0][[epsilon].sub.0])] propagates above the interface, where [omega] is the frequency, [[epsilon].sub.0] is the dielectric permittivity of vacuum, and [[mu].sub.0] is the magnetic permeability of vacuum.

3. Maxwell's Equations

The solutions of Maxwell's equations [4] for the considered geometry can be taken in the exponential form exp(j[omega]t - j[k.sub.x]x - j[k.sub.y]y - j[k.sub.z]z). Thus Maxwell's equations [2] in the scalar form are written as follows:

[mathematical expression not reproducible] (3)

Here, it is taken into account that

[mathematical expression not reproducible] (4)

Since the statement of the problem is shown in Figure 1, we have [k.sub.x] = 0. Therefore,

[mathematical expression not reproducible] (5)

The analogous transformations can be presented for the electrical field.

The solutions of (3) for this case give us two independent waves: an ordinary wave containing [H.sub.y]-, [H.sub.z]-, [E.sub.x]- components and an extraordinary wave containing the [H.sub.x]-, [E.sub.y]-, [E.sub.z]-components. Note that the wavenumbers of forward and backward waves are equal [2] for the ordinary wave, and these wavenumbers are not equal for the extraordinary wave. Now let us consider the behavior of these waves in detail.

4. Ordinary Wave

The geometry of the considered problem is shown in Figure 2. Using the first, fifth, and sixth equations of (3), we obtain the relations between the components of the electromagnetic fields for an ordinary wave (a wave containing [H.sub.y]-, [H.sub.z]-, [E.sub.x]- components) in an anisotropic medium:

[H.sub.y] = 1/[[rho].sub.0] [square root of ([[epsilon].sub.xx]/[[epsilon].sub.0] - 1)][E.sub.x], [H.sub.z] = - 1/[[rho].sub.0] [E.sub.x], (6)

where [[rho].sub.0] = [square root of ([[mu].sub.0]/[[epsilon].sub.0])] is wave impedance of vacuum.

It is important that the normal magnetic field component [H.sub.z] and the tangential electric field component [E.sub.x] do not depend on the frequency and an inclination angle of the anisotropy axis. The frequency dependence of the tangential component of the magnetic field [H.sub.y] arises only in the case of frequency dispersion of the electric permittivity components in (1) [5].

Let us find the normal component of the wave vector [k.sub.z]. Taking into account the continuity conditions for the tangential components of the wave vector, we obtain [k.sub.y] = [k.sub.0]. The normal component of the wave vector in an anisotropic medium is found by substituting (4) into the first equation in (3):

[K.sub.z] = [omega] [square root of ([[epsilon].sub.0]/[[epsilon].sub.xx] - [[epsilon].sub.0])] ([[epsilon].sub.xx] [square root of ([[mu].sub.0]/[[epsilon].sub.0])] - [square root of ([[epsilon].sub.0][[mu].sub.0])]) = [omega] [square root of ([[mu].sub.0])] [square root of ([[epsilon].sub.xx] - [[epsilon].sub.0])]. (7)

It is followed from (7) that the refraction angle is

[alpha] = arctan [k.sub.y]/[k.sub.z] = arctan [square root of ([[epsilon].sub.0]/[[epsilon].sub.xx] - [[epsilon].sub.0])]. (8)

The Poynting vector components are calculated by using the well-known method [4]:

[S.sub.z] = [E.sub.x][H.sup.*.sub.y] = 1/[[rho].sub.0] [square root of ([[epsilon].sub.xx]/[[epsilon].sub.0] - 1)][H.sup.2.sub.x], (9)

[S.sub.y] = -[E.sub.x][H.sup.*.sub.z] = [square root of ([[epsilon].sub.0]/[[mu].sub.0])][H.sup.2.sub.x] = - 1/[[rho].sub.0][H.sup.2.sub.x]. (10)

Thus the ordinary wave has both the real normal component of the wave vector (7) and the real normal component of the Poynting vector (9). Consequently the bulk wave must exist in an anisotropic medium. And now the power transmission angle (see (9) and (10)) is

[[alpha].sub.s] = arctan [S.sub.y]/[S.sub.z] = arctan [square root of ([[epsilon].sub.0]/[[epsilon].sub.xx] - [[epsilon].sub.0])]. (11)

It is seen from (8) and (11) that [[alpha].sub.s] = [alpha] for the considered wave type. Therefore, both the power flow and the wave front propagate in the same direction analogously to a usual isotropic medium.

5. Extraordinary Wave

The geometry of the considered problem is presented in Figure 3. Solving the second and third equations of (3) simultaneously, we obtain the relations between the field components for the extraordinary wave ([E.sub.x], [H.sub.y], [H.sub.z]):

[mathematical expression not reproducible] (12)

The normal component of the wave vector can be found by substituting the expressions (12) and (13) in the fourth equation of (3) that is equal to

[mathematical expression not reproducible] (13)

The forward wave is described by (13) with a "plus" and the backward wave with a "minus." Obviously, the angle between wave front propagation and the normal is

[alpha] = arctan [k.sub.y]/[k.sub.z]. (14)

The Poynting vector components of the extraordinary wave are equal to

[S.sub.z] = [E.sub.x][H.sup.*.sub.y] = - [k.sub.z][[epsilon].sub.zz] + [k.sub.y][[epsilon].sub.yz]/[omega]([[epsilon].sub.yy][[epsilon].sub.zz] - [[epsilon].sub.yz][[epsilon].sub.zy]) [E.sup.2.sub.x], (15)

[S.sub.y] = -[E.sub.x][H.sup.*.sub.z] = - [k.sub.y][[epsilon].sub.yy] + [k.sub.z][[epsilon].sub.zy]/[omega]([[epsilon].sub.yy][[epsilon].sub.zz] - [[epsilon].sub.yz][[epsilon].sub.zy]) [E.sup.2.sub.x], (16)

where [H.sup.*.sub.y] is the complex conjugate of [H.sub.y]. The angle between the electromagnetic power flow propagation and a normal to the interface is as follows:

[[alpha].sub.s] = arctan [S.sub.y]/[S.sub.z] = arctan [k.sub.y][[epsilon].sub.yy] + [k.sub.z][[epsilon].sub.zy]/[k.sub.z][[epsilon].sub.zz] + [k.sub.y][[epsilon].sub.yz]. (17)

The analysis of the obtained expressions shows that the transmitted wave has the real normal (13) and tangential (16) components of the Poynting vector. Therefore the power of this wave propagates under an oblique angle in an anisotropic medium. The power propagation angle (15) does not coincide with the wave front propagation angle (14).

6. Behavior of a Reflected Wave

It has been shown earlier that a plane harmonic wave propagating parallel to an interface between vacuum and an anisotropic medium can excite a bulk wave propagating under an arbitrary angle within an anisotropic medium. Now it is necessary to describe behavior of a reflected wave in vacuum if it existed.

Let us take into account the fact that the tangential components of the wavenumbers of incident, reflected, and refracted waves must be equal one to other [k.sub.inc] = [k.sub.y refr] = [k.sub.y refl] in accordance with Snell's Law [4] and these wavenumbers are also equal to the wavenumber in vacuum [k.sub.0] = [omega][square root of ([[epsilon].sub.0][[mu].sub.0])]. Then a reflected wave must propagate along an interface or perpendicular to an interface. Otherwise the frequency of a reflected wave must be more than the frequency of an incident wave to satisfy the boundary conditions. Indeed if a reflected wave is propagating under an oblique angle and the tangential component is equal to the wavenumber in vacuum then [k.sub.refl] = [square root of ([k.sup.2.sub.y] + [k.sup.2.sub.z])] = [square root of ([k.sup.2.sub.0] + [k.sup.2.sub.z])] > [k.sub.0].

Moreover a refracted wave has the tangential component of magnetic field for an ordinary wave and the tangential component of the electric field for an extraordinary wave. Therefore such component must be in a reflected wave as this component is absent in an incident wave. But such phenomenon is impossible when a reflected wave propagates along an interface. Based on these considerations, we assume that a reflected linearly polarized wave should propagate perpendicular to an interface between vacuum and an anisotropic medium.

7. Coefficients of Propagation and Reflection of an Ordinary Wave

Let us consider the ideal case of plane harmonic wave propagation above an interface between vacuum and a uniaxial anisotropic medium. It follows from (6) that a transmitting wave has the components [H.sub.y], [H.sub.z], [E.sub.x]. Thus a reflected wave must also contain the [H.sub.y]-component as this is absent in an incident wave.

The geometry of the problem is shown in Figure 2. Here [k.sub.0] is the wavenumber in vacuum; [H.sub.inc], [E.sub.inc] are the field components of an incident wave propagating along an interface; [H.sub.refl], [E.sub.refl] are the field components of the reflected wave; [H.sub.refr], [E.sub.refr] are the field components of the transmitting wave.

Let us write the continuity conditions for the field components at the interface between vacuum and an anisotropic medium [4]:

[E.sub.inc] + [E.sub.refl] = [E.sub.x refr], - [H.sub.y refl] = [H.sub.y refr] + J, [B.sub.z inc] = [B.sub.z refr]. (18)

Take into account the fact that, in free spaces [H.sub.y] = (1/[[rho].sub.0])[E.sub.x] and [H.sub.z] = -(1/[[rho].sub.0])[H.sub.x], we can obtain from (18)

1 + R = T, - 1/[[rho].sub.0]R = [X.sub.y]T + J/[E.sub.inc], - [[epsilon].sub.0]/[[rho].sub.0] = [[epsilon].sub.0][X.sub.z]T. (19)

Here [[rho].sub.0] = [square root of ([[mu].sub.0]/[[epsilon].sub.0])] is the impedance of vacuum, J is the surface current density at the interface, R = [E.sub.x refl]/[E.sub.x inc] is the reflection coefficient, and T = [E.sub.x refr] /[E.sub.x inc] is the transmission coefficient. Now we can find the transmission and reflection coefficients, as well as the surface charge density:

T = 1, (20)

R = 0, (21)

J = -[[rho].sup.2.sub.0][E.sub.inc]. (22)

Thus, if a harmonic plane wave containing the components [H.sub.y], [H.sub.z], [E.sub.x] propagates along an interface between vacuum and an anisotropic medium, then there is no reflected wave (21), the transmission coefficient is equal to one (20), and also there are the surface currents (22) at the interface. These currents are actually bias currents and these are associated with the molecule polarization phenomenon in an electromagnetic field [4].

8. Coefficients of Transmission and Reflection of an Extraordinary Wave

The geometry of the problem is given in Figure 3. The boundary conditions for the electromagnetic field components in this case are as follows:

[H.sub.x inc] + [H.sub.x refl] = [H.sub.x refr], - [E.sub.y refl] = [E.sub.y refr], - [[epsilon].sub.0] [E.sub.z inc] = [[epsilon].sub.zy][E.sub.y refr] + [[epsilon].sub.zz][E.sub.z refr] + [sigma]. (23)

Taking into account the relations [E.sub.x] = - [[rho].sub.0][H.sub.y], [E.sub.x] = [[rho].sub.0][H.sub.z] [4] and introducing the reflection R = [H.sub.x refl]/[H.sub.x inc] and transmission T = [H.sub.x refr]/[H.sub.x inc] coefficients, we can write

1 - R = T, - [[rho].sub.0]R = [X.sub.y]T, - [[mu].sub.0][[rho].sub.0] = [[epsilon].sub.zy][X.sub.y]T + [[epsilon].sub.zz][X.sub.z]T + [sigma]/[H.sub.x inc], (24)

instead of (23). Now it can be found from (24) that

[mathematical expression not reproducible] (25)

The analysis of the obtained expression shows that the parts of an electromagnetic wave transmitting in a medium and reflected from the interface depend on the parameters of the medium and do not depend on frequency. Indeed, [k.sub.y] and [k.sub.z] include w in the first degree. Therefore w can be reduced in (25). The frequency dependence is observed only in the case of frequency dispersion of the dielectric permittivity [5]. Moreover, the coefficients in (23) depend on the inclination angle of the anisotropy axis. At the interface, there is electrical charge which is a result of the polarization of the molecules in the electromagnetic field.

9. Numerical Results

In this section, the numerical example of the extraordinary electromagnetic wave behavior in the considered particular case is presented. We consider the quartz medium with the typical parameters [[epsilon].sub.11] = 4.5, [[epsilon].sub.33] = 4.6 [7, 8] and without frequency dispersion. Only the extraordinary wave is studied because the behavior of an ordinary wave is analogous to the case of an isotropic medium. Indeed the wavenumbers of the forward and backward waves are the same as the ordinary wave.

In Figure 4(a), the dependence of the transmission angle of the wave front on the inclination angle of the anisotropy axis is presented. In Figure 4(b), the dependence of the power transmission angle on the inclination angle of the anisotropy axis is presented. It is seen that the angles have low dependence on [beta] (the reflection angle is [[alpha].sub.refr] = 27.799[degrees] / 28.12[degrees] if the inclination angle is [beta] = -180[degrees] / 180[degrees]), and the directions of the wave front propagation and power propagation are different.

It has been also verified by numerical calculations that there is no frequency dependence of these angles. Dependence can be seen only if frequency dispersion of the constitutive parameters [[epsilon].sub.11], [[epsilon].sub.33] is taken into account. For example, it can be plasma medium in a resonance domain.

The dependencies of the transmission and reflection coefficients for the medium parameters given above are presented in Figure 5. It is seen from Figure 5 that the maximum transmission corresponds to the tangential orientation of the anisotropy axis and the maximum reflection corresponds to the normal orientation of the anisotropy axis. Moreover, the reflection and transmission coefficients are symmetrical functions with response to the normal orientation of the anisotropy axis in spite of the fact the normal wave vector constant is asymmetric.

We must also note that T is not the transmission coefficient in the classical meaning [4]. Indeed a transmitted wave power also depends on the relation between the components [E.sub.x], [H.sub.y], [H.sub.z].

10. Conclusions

In this paper, we consider the behavior of the harmonic plane wave propagating parallel to the interface between vacuum and a uniaxial anisotropic medium. It is shown that the wave propagating in an anisotropic medium has both the normal component of the wave vector and the normal component of the Poynting vector.

It is also proved that in this case the reflected wave should propagate perpendicularly to an interface. At the same time, the reflected wave is absent for an ordinary wave. Moreover, in this case, the reflection and transmission coefficients do not depend on the frequency if frequency dispersion of material parameters is absent.

It is important to dwell on the practical applications of this phenomenon. In particular, we note that it must be taken into account the calculation and design of microwave and optical communication systems. Furthermore, penetration effect enables the creation of new microwave and optical devices such as a power divider [1], a frequency detector with isolator properties [9], and a duplex isolator [10,11] for telecommunication systems.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this manuscript.

https://doi.org/ 10.1155/2017/4082948

References

[1] K. Vytovtov and L. Mospan, "Penetration effect in gyrotropic slab: theory and applications," Journal of the Optical Society of America A: Optics and Image Science, and Vision, vol. 29, no. 6, pp. 877-882, 2012.

[2] K. Vytovtov and S. S. Zouhdi, "Exotic reflection of plane waves by anisotropic structures," in Proceedings of the Metamaterials '2012: The 6th International Congress on Advanced Electromagnetic Materials in Microwaves and Optics, Metamaterials, St. Petersburg, Russia, 2012.

[3] S. A. Andilevko, "A hydrodynamics model of super-deep penetration effect," Journal de Physique IV, vol. 7, pp. 729-734, 1997.

[4] M. Born and E. Wolf, Principles of Optics, 7th Edition, Cambridge, UK, 1999.

[5] T. H. Stix, Waves in Plasma, Springer Science & Business Media, Media, 1992.

[6] K. Shoemake, "Animating rotation with quaternion curves," in Proceedings of the 12th Annual Conference on Computer Graphics and Interactive Techniques SIGGRAPH '85, pp. 245-254.

[7] L. Ferrari, C. Wu, D. Lepage, X. Zhang, and Z. Liu, "Hyperbolic metamaterials and their applications," Progress in Quantum Electronics, vol. 40, pp. 1-40, 2015.

[8] D. A. Robinson, S. B. Jones, J. M. Wraith, D. Or, and S. P. Friedman, "A review of advances in dielectric and electrical conductivity measurement in soils using time domain reflectometry," Vadose Zone Journal, vol. 2, no. 4, pp. 444-475, 2003.

[9] K. Vytovtov, V. Gnatushenko, and S. Zouhdi, "Terahertz range double direction isolator based on stratified antiferromagnetic-dielectric structures. Theoretical investigation," Elektronika, vol. 2, pp. 75-78, 2015.

[10] K. Vytovtov, S. Zouhdi, R. Dubrovka, and V. Hnatushenko, "The terahertz controlled duplex isolator: physical grounds and numerical experiment," International Journal of Microwave Science and Technology, vol. 2016, Article ID 1468508, 7 pages, 2016.

[11] I. O. Barabanov, N. S. Maltseva, and E. A. Barabanova, "Switching cell for information transmission optical systems," in Proceedings of the 2016 International Conference on Actual Problems of Electron Devices Engineering APEDE 2016, pp. 343-347, September 2016.

K. Vytovtov and O. Pischin

Astrakhan State Technical University, Tatischeva 16, Astrakhan 414056, Russia

Correspondence should be addressed to K. Vytovtov; vytovtov_konstan@mail.ru

Received 15 February 2017; Revised 26 March 2017; Accepted 2 May 2017; Published 23 May 2017

Academic Editor: Mustapha C. E. Yagoub

Caption: Figure 1: Geometry of the problem.

Caption: Figure 2: Geometry of the problem of the propagation and reflection of an ordinary wave.

Caption: Figure 3: Geometry of the problem of the propagation and reflection of an extraordinary wave.

Caption: Figure 4: The dependencies of the transmission angles on the inclination angle.

Caption: Figure 5: The dependencies of the transmission and reflection coefficients on the inclination angle.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Vytovtov, K.; Pischin, O. |

Publication: | International Journal of Microwave Science and Technology |

Article Type: | Report |

Date: | Jan 1, 2017 |

Words: | 3758 |

Topics: |