**2**, 6

https://doi.org/10.1051/epjam/2015009

Research Article

## Nihility in non-reciprocal bianisotropic media

Department of Radio Science and Engineering, Aalto University, P.O. Box 13000, 00076
Aalto, Finland

^{*} e-mail: younes.radi@aalto.fi

Received:
10
August
2015

Accepted:
6
November
2015

Published online: 23 December 2015

Here we consider electromagnetic response of non-reciprocal bianisotropic materials in some extreme regimes. The magneto-electric coupling is modeled by symmetric and antisymmetric uniaxial dyadics, which correspond to the so called artificial Tellegen media and moving media, respectively. Extreme electromagnetic properties of uniaxial non-reciprocal bianisotropic materials in the limiting case of *nihility*, when both permittivity and permeability of the media tend to zero, and only the magneto-electric parameters define the material response, are studied. Among other interesting effects, we show that the moving nihility materials provide the extreme asymmetry in the phase shift of transmitted waves propagating along the opposite directions. Furthermore, we reveal a possibility to create an angular filter with extreme sensitivity to the incidence angle, also using moving nihility slabs.

Key words: Bianisotropic media / Wave propagation

*© Y. Ra’di et al., Published by EDP Sciences, 2015*

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## 1 Introduction

The metamaterial paradigm offers possibilities to design and realize electromagnetic materials with desired values of constitutive parameters (naturally, with some limitations). It is of theoretical and practical interest to explore the properties of artificial materials with extreme values of material parameters, either very large or very small. Such extreme materials show interesting, often extreme, electromagnetic properties [1]. Recently, materials with very small values of the permittivity and/or permeability attracted attention due to some interesting physical effects and application possibilities, see e.g. [2–4]. General properties of isotropic media with either very small or very large parameter values were discussed in reference [1]. In this paper we study some extreme cases of bianisotropic materials, focusing on the case of non-reciprocal materials.

It is well known that the difference of the plane-wave propagation factors *β* of the two (circularly polarized) eigenwaves in isotropic chiral media is defined by the chirality parameter *κ* in the Tellegen formalism as [5] *β* = *k*
_{0} (*N* ± *κ*), where is the free-space wavenumber and is the refractive index. In the extreme case when *N* → 0 but *κ* ≠ 0, the propagation factors of the two eigenmodes differ by sign, maximizing birefringence of medium. In particular, mixtures of optimal helices [6, 7] with *κ* ≈ 1 realize effective media with the propagation factor *β* = −*k*
_{0} for one of the circular polarizations, while the same medium is transparent for the orthogonal circular polarization [8]. This extreme-parameter medium is called *chiral nihility* [9]. The concept of chiral nihility leads to understanding of the chiral route to negative refraction and superlensing with the use of chiral structures [9, 10]. In the context of this work, we can say that the asymmetry in the propagation constants of the two eigenwaves comes to its extreme (the two values differ by sign) in the limiting case of chiral nihility. The wave impedances of the two eigenwaves are always the same as they do not depend on the chirality parameter *κ*.

More recently, it was understood that there exist optimal parameters (the optimal shapes and sizes of meta-atom inclusions) also for another fundamental class of reciprocal bianisotropic media: omega materials [11]. In contrast to chiral media where the chirality breaks the symmetry of the propagation constants of eigenwaves while the wave impedances are not affected, in omega materials the properties are dual: Magneto-electric coupling breaks the symmetry of the wave impedances of the counter-propagating modes while the propagation constants remain symmetric [12]. It was shown that *omega nihility* media have extreme properties in reflection from a planar slab: the reflection coefficient changes its sign when one or the other face of the slab is illuminated [13].

The goal of this work is to study the fundamental properties of non-reciprocal bianisotropic materials in the case when the permittivity and permeability are negligibly small as compared with the magneto-electric coupling coefficients and understand what kinds of extreme wave properties can be expected for such special materials.

## 2 General solution for plane waves propagating in uniaxial non-reciprocal media

In this work we concentrate on the properties of bianisotropic materials in which the magneto-electric coupling is due to non-reciprocal effects. Thus, we assume that the material under study is non-chiral and there is no omega coupling. In this case, the appropriate constitutive relations for non-reciprocal bianisotropic media having uniaxial symmetry (the axis direction is along the unit vector **z**
_{0}) are described by uniaxial material dyadics as [12](1)

Here, and are the material permittivity and permeability, is the transverse unit dyadic, and is the 90 degree rotator in the *x* − *y* plane. Parameters *α* and *χ*
_{(t, n)} measure the strength of magneto-electric coupling and are called the effective velocity (in artificial moving media) and the Tellegen parameter, respectively. For lossless materials, these coupling parameters are real-valued.

In this paper we will derive field equations, study the propagation factors and the wave impedances of the eigenwaves in the general case of uniaxial non-reciprocal bianistropic media, and then consider the extreme cases of nihility materials. In the following derivations we use the general method for analysing plane waves in uniaxial materials, described in [12, 14]. We look for time-harmonic plane-wave solutions of the Maxwell equations(2)in form *e*^{–jk·r}, where the wavevector **k** = *β***z**_{0} + **k**_{t} (with **z**_{0} · **k**_{t} = 0). As we are considering structures with the uniaxial symmetry, it is convenient to split the field vectors into the longitudinal and transverse parts with respect to the material axis **z**_{0}: (3)(here the vectors marked by index t are orthogonal to **z**_{0}). Substituting these fields into the Maxwell equations and performing the two-dimensional Fourier transform in the transverse plane (*x* − *y*), after some mathematical manipulations, we arrive at a system of two equations for the transverse components of electric and magnetic fields: (4)

Here is the free-space wave impedance, and **k**_{t} stands for the two-dimensional Fourier variable (the transverse component of the propagation factor). Eliminating the magnetic field from (4), we come up with a second-order equation for the transverse electric field component **E**_{t}: (5)where the coefficients read(6)with . To derive the wave propagation factors *β*, we assume a plane-wave solution in the form exp(–*jβz*) and equate the determinant of (5) to zero. We get a bi-quadratic eigenvalue equation for the propagation factor *β*: (7)where . Solutions of this equation for the wave propagation factor *β* read (8)

To get a complete description of eigenwaves, we also need to find the wave impedances for waves travelling in these materials. Dyadic wave impedances define the relations between the transverse electric and magnetic fields components in plane waves propagating in an unbounded medium. We define the relation between transverse electric and magnetic fields as (9)

In this definition, the eigenvalues of which have positive real parts correspond to waves propagating along the *z*-axis, and eigenvalues having negative real part correspond to eigenwaves propagating in the opposite direction. For the four eigensolutions, the wave impedances can be readily found from (4) by substituting the corresponding propagation factors from (8). After some dyadic algebra we find (10)where *β* is given by (8).

## 3 Tellegen and moving nihility

In the previous section, we derived the general solutions for the propagation constant and wave impedance in non-reciprocal uniaxial bianisotropic media. Next, we will study what extreme properties, we can expect if the permittivity and permeability of the materials are negligibly small (*bianisotropic nihility*).

### 3.1 Tellegen nihility

To find the wave propagation factor for the case of Tellegen nihility, we just need to substitute *α* = 0 in (8). Let us consider the Tellegen nihility case where *ϵ*
_{t}, *ϵ*
_{n}, *μ*
_{t}, and *μ*
_{n} tend to zero. In this case the wave propagation constant reads (11)

As it can be seen from this equation, for lossless media the propagation constants can be purely real, purely imaginary, or zero. For the cases where *β* ≠ 0, for arbitrary propagation directions all the components of the wave impedance dyadics tend to infinity meaning that these modes do not couple to free-space modes at the sample boundaries. For the case where *β* = 0, all the components of the wave impedance dyadics tend to infinity except the component corresponding to which tends to zero.

Another interesting extreme case is the case where only the transverse components of permittivity and permeability tend to zero *ϵ*
_{t}, *μ*
_{t} → 0 but the normal components of the permittivity and permeability dyadics may have any arbitrary non-zero values. The propagation constant for this case reads(12)

In this case, for some design parameters, the propagation constants can even take real values. All the components of the wave impedance dyadics tend to infinity except component which gets a purely real value. This interesting wave-propagation phenomenon in Tellegen media happens only in the case of non-axial propagation. For axial propagation, the propagation constants are purely imaginary and all components of the wave impedance dyadics tend to infinity.

### 3.2 Moving nihility

For the case of an artificial moving medium, by substituting χ_{t} = χ_{n} = 0 in (8), the four values of the propagation constant are given by(13)

Now let us consider the moving nihility case where *ϵ*
_{t}, *ϵ*
_{n}, *μ*
_{t}, and *μ*
_{n} tend to zero with the same rate. In the case of oblique propagation (*k*
_{t} ≠ 0), the wave propagation factors read (14)while the wave impedances are (15)however, for the case of axial propagation (*k*
_{t} = 0), the propagation factors and wave impedances simplify to (16)

These interesting results show that such a medium can support both forward and backward waves propagating along the axis. In other words, it behaves as a Veselago medium [15] for waves propagating in one direction along the axis and as vacuum for waves propagating in the opposite direction. Let us consider a material slab of thickness *d* in free space, filled with a moving nihility material under illumination of a normally incident plane wave. Such a layer is matched for waves hitting any of its sides, however, the transmission coefficient *T* is different: (17)

Similar results were already discovered for metasurfaces with moving magneto-electric coupling [16].

Another interesting example is the case where only transverse components of permittivity and permeability tend to zero. The propagation constant equals (18)

while the wave impedance becomes (19)

For the axial propagation the wave impedance equals to that of free space for both propagation directions:(20)

It is seen that for axially propagating plane waves in a moving nihility medium we have *η*
_{0} wave impedance for waves propagating in both directions, while for waves propagating in any other direction the diagonal components of the impedance dyadics tend to infinity. A slab of such material constitutes an ideal angular filter.

## 4 Conclusion

Exotic and extreme properties of non-reciprocal nihility materials have been considered. We have studied Tellegen and moving nihility media. It was found that, for specific material parameters, oblique wave propagation is possible in a Tellegen nihility medium while axial propagation is not possible. It was shown that lossless moving nihility materials exhibit extreme asymmetry in the phase propagation constants for oppositely directed waves, so that the phase shift for waves transmitted through a slab in one direction is negative with respect to the phase shift for the oppositely-bound waves. We have also revealed an effect of extreme sensitivity of the reflection coefficient on the incidence angle for moving nihility slabs. It was found that the rate in which *ϵ* and *μ* tend to zero is, for some cases, a factor which determines the electromagnetic response of the composite.

## References

- A. Sihvola, S. Tretyakov, A. de Baas, Metamaterials with extreme material parameters, J. Commun. Technol. Electron. 52 (2007) 986–990. [CrossRef] (In the text)
- R.W. Ziolkowski, Propagation in and scattering from a matched metamaterial having a zero index of refraction, Phys. Rev. E 70 (2004) 046608. [CrossRef] (In the text)
- A.M. Mahmoud, N. Engheta, Wave-matter interactions in epsilon-and-mu-near-zero structures, Nature Commun. 5 (2014) 5638. [CrossRef]
- N. Engheta, Pursuing near-zero response, Science 340 (2013) 286–287. [CrossRef] (In the text)
- I.V. Lindell, A.H. Sihvola, S.A. Tretyakov, A.J. Viitanen, Electromagnetic waves in chiral and bi-isotropic media, Artech House, Norwood, MA, 1994. (In the text)
- I.V. Semchenko, S.A. Khakhomov, A.L. Samofalov, Radiation of circularly polarized electromagnetic waves by the artificial flat lattice with two-turn helical elements, in Proc. of the 10th International Conference on Complex Media and Metamaterials, Ghent, Belgium, 2004, pp. 236–239. (In the text)
- I.V. Semchenko, S.A. Khakhomov, A.L. Samofalov, Optimal shape of spiral: equality of dielectric, magnetic and chiral properties, in Proc. of NATO ARW & META’08, Marrakesh, Morocco, 2008, pp. 71–80. (In the text)
- I.V. Semchenko, S.A. Khakhomov, S.A. Tretyakov, Chiral metamaterial with unit negative refraction index, Eur. Phys. J. Appl. Phys. 46 (2009) 32607. [CrossRef] [EDP Sciences] (In the text)
- S. Tretyakov, I. Nefedov, A. Sihvola, S. Maslovski, C. Simovski, Waves and energy in chiral nihility, J. Electromagn. Waves Appl. 17 (2003) 695–706. [CrossRef] (In the text)
- S. Tretyakov, A. Sihvola, L. Jylhä, Backward-wave regime and negative refraction in chiral composites, Photon. Nanostruct. Fund. Appl. 3 (2005) 107–115. [CrossRef] (In the text)
- S.A. Tretyakov, The optimal material for interactions with linearly-polarized electromagnetic waves, in Proc. of The Fourth International Congress on Advanced Electromagnetic Materials in Microwaves and Optics (Metamaterials’2010), Karlsruhe, Germany, 2010, pp. 65–67. (In the text)
- A.N. Serdyukov, I.V. Semchenko, S.A. Tretyakov, A. Sihvola, Electromagnetics of bi-anisotropic materials: theory and applications, Gordon and Breach Science Publishers, Amsterdam, 2001. (In the text)
- Y. Radi, S.A. Tretyakov, Electromagnetic phenomena in omega nihility media, in Proc. of The Sixth International Congress on Advanced Electromagnetic Materials in Microwaves and Optics (Metamaterials’2012), St. Petersburg, Russia, 2012, pp. 764–766. (In the text)
- S.A. Tretyakov, A.A. Sochava, Reflection and transmission of plane electromagnetic waves in uniaxial bianisotropic materials, Int. J. Infrared Millim. Waves 15 (1994) 829–856. [CrossRef] (In the text)
- V.G. Veselago, The electrodynamics of substances with simultaneously negative values of ϵ and μ, Sov. Phys. Usp. 10 (1968) 509–514. [CrossRef] (In the text)
- Y. Ra’di, V.S. Asadchy, S.A. Tretyakov, One-way transparent sheets, Phys. Rev. B 89 (2014) 075109. [CrossRef] (In the text)

**Cite this article as**: Ra’di Y, Asadchy VS & Tretyakov SA: Nihility in non-reciprocal bianisotropic media. EPJ Appl. Metamat. 2015, **2**, 6.