text:
# Oblique frozen modes in periodic layered media

## I Introduction

Electromagnetic properties of periodic stratified media have been a subject of extensive research for decades (see, for example, Strat1 ; Strat2 ; Strat3  and references therein). Of particular interest has been the case of periodic stacks ($`1D`$ photonic crystals) made up of lossless dielectric components with different refractive indices. Photonic crystals with one-dimensional periodicity had been widely used in optics long before the term ”photonic crystals” was invented.

Let us look at the classical problem of a plane electromagnetic wave incident on the surface of semi-infinite plane-parallel periodic array, as shown in Fig. 1.

The well known effects of the slab periodicity are: (i) the possibility of omnidirectional reflectance when the incident radiation is reflected by the slab, regardless of the angle of incidence; (ii) the possibility of negative refraction , when the tangential component of the energy flux $`\stackrel{}{S}_T`$ of transmitted wave is antiparallel to that of the incident wave; (iii) dramatic slowdown of the transmitted wave near photonic band edge frequency, where the normal component of the transmitted wave group velocity $`\stackrel{}{u}`$ vanishes along with the respective energy flux $`\stackrel{}{S}_T`$. The extensive discussion on the subject and numerous references can be found in Omni1 ; Omni2 ; Omni3 ; NegRefr1 ; NegRefr2 ; Joann1 ; Scalora1 ; Agua01 ; Scalora2 ; Scalora97 . All the above effects can occur even in the simplest case of a semi-infinite periodic array of two isotropic dielectric materials with different refractive indices, for example, glass and air. The majority of known photonic crystals fall into this category. The introduction of *dielectric anisotropy*, however, can bring qualitatively new features to electromagnetic properties of periodic stratified media and open up new opportunities for practical applications (see, for example, a recent publication Manda03 ). One of such phenomena is the subject of this work.

### I.1 The Axially Frozen Mode (AFM)

Consider a semi-infinite periodic stack with at least one of the constituents being an anisotropic dielectric material with oblique orientation of anisotropic axis. A simple example of such an array is presented in Fig. 2.

We will show that under certain physical conditions, a monochromatic plane wave incident on the semi-infinite slab is converted into an abnormal electromagnetic mode with huge amplitude and nearly tangential energy density flux, as illustrated in Fig. 3.

Such a wave will be referred to as the *Axially Frozen Mode* (AFM). The use of this term is justified because the normal (axial) component $`u_z`$ of the respective group velocity becomes vanishingly small, while the amplitude of the AFM can exceed the amplitude of the incident plane wave by several orders of magnitude.

The group velocity $`\stackrel{}{u}`$ of the AFM is parallel to the semi-infinite slab boundary and, therefore, the magnitude of the tangential component $`\left(\stackrel{}{S}_T\right)_{}`$ of the respective energy density flux $`\stackrel{}{S}_T`$ is overwhelmingly larger than the magnitude of the normal component $`\left(\stackrel{}{S}_T\right)_z`$. But, although $`\left(\stackrel{}{S}_T\right)_z\left(\stackrel{}{S}_T\right)_{}`$, the normal component $`\left(\stackrel{}{S}_T\right)_z`$ of the energy density flux inside the slab is still comparable with that of the incident plane wave in vacuum. This property persists even if the normal component $`u_z`$ of the wave group velocity inside the slab vanishes, i.e.,

$$\left(\stackrel{}{S}_T\right)_z>0\text{,}\text{if}u_z=0.$$
(1)
The qualitative explanation for this is that the infinitesimally small value of $`u_z`$ is offset by huge magnitude of the energy density $`W`$ in the AFM. As the result, the product $`u_zW`$, which determines the normal component $`\left(\stackrel{}{S}_T\right)_z`$ of the energy flux, remains finite. The above behavior is totally different from what happens in the vicinity of a photonic band edge, where the normal component $`u_z`$ of the wave group velocity vanishes too. Indeed, let us introduce the transmittance $`\left(\tau \right)`$ and the reflectance $`\left(\rho \right)`$ of a lossless semi-infinite slab

$$\tau =1\rho =\frac{\left(\stackrel{}{S}_T\right)_z}{\left(\stackrel{}{S}_I\right)_z},\rho =\frac{\left(\stackrel{}{S}_R\right)_z}{\left(\stackrel{}{S}_I\right)_z}.$$
(2)
In line with Eq. (1), in the AFM regime the transmittance $`\tau `$ remains significant and can be even close to unity, as shown in an example in Fig. 4(a). In other words, the incident plane wave enters the slab with little reflectance, where it turns into an abnormal AFM with infinitesimally small normal component of the group velocity, huge amplitude, and huge tangential component of the energy density flux. By contrast, in the vicinity of a photonic band edge (at frequencies near $`\omega =\omega _b`$ in Fig. 4(a)), the transmittance of semi-infinite slab always vanishes, along with the normal component $`u_z`$ of the wave group velocity.

It turns out that at a given frequency $`\omega _0`$ the AFM regime can occur only for a special direction $`\stackrel{}{n}_0`$ of the incident plane wave propagation

$$\stackrel{}{n}_0=\stackrel{}{n}_0\left(\omega _0\right).$$
(3)
This special direction of incidence always makes an oblique angle with the normal $`z`$ to the layers. To find $`\stackrel{}{n}_0`$ for a given $`\omega _0`$ or, conversely, to find $`\omega _0`$ for a given $`\stackrel{}{n}_0`$, one has to solve the Maxwell equations in the periodic stratified medium. This problem will be addressed in Section 3. In Section 2 we consider the relation between the AFM regime and the singularity of the electromagnetic dispersion relation responsible for such a peculiar behavior. If the frequency $`\omega `$ and the direction of incidence $`\stackrel{}{n}`$ do not match explicitly as prescribed by Eq. (3), the AFM regime will be somewhat smeared.

### I.2 The vicinity of the AFM regime

Let $`\mathrm{\Psi }_T\left(z\right)`$ be the transmitted electromagnetic field inside the semi-infinite slab (the explicit definition of $`\mathrm{\Psi }_T\left(z\right)`$ is given in Eqs. (38) and (109)). It turns out that in the vicinity of the AFM regime, $`\mathrm{\Psi }_T\left(z\right)`$ is a superposition of the extended and evanescent Bloch eigenmodes

$$\mathrm{\Psi }_T\left(z\right)=\mathrm{\Psi }_{ex}\left(z\right)+\mathrm{\Psi }_{ev}\left(z\right),z>0,$$
(4)
where $`\mathrm{\Psi }_{ex}\left(z\right)`$ is an extended mode with $`u_z>0`$, and $`\mathrm{\Psi }_{ev}\left(z\right)`$ is an evanescent mode with $`\mathrm{}k_z>0`$. As shown in an example in Fig. 5, both the contributions to $`\mathrm{\Psi }_T\left(z\right)`$ have huge and nearly equal and opposite values near the slab boundary, so that their superposition (4) at $`z=0`$ is small enough to satisfy the boundary condition (111). As the distance $`z`$ from the slab boundary increases, the evanescent component $`\mathrm{\Psi }_{ev}\left(z\right)`$ decays exponentially, while the amplitude of the extended component $`\mathrm{\Psi }_{ex}\left(z\right)`$ remains constant and huge. As the result, the field amplitude $`\left|\mathrm{\Psi }_T\left(z\right)\right|^2`$ reaches its huge saturation value $`\left|\mathrm{\Psi }_{ex}\right|^2`$ at a certain distance from the slab boundary (see Eqs. (120), (121) and (122)).

When the direction of incidence $`\stackrel{}{n}`$ tends to its critical value $`\stackrel{}{n}_0`$ for a given frequency $`\omega _0`$, the respective saturation value $`\left|\mathrm{\Psi }_{ex}\right|^2`$ of the AFM amplitude $`\left|\mathrm{\Psi }_T\left(z\right)\right|^2`$ diverges as $`\left|\stackrel{}{n}\stackrel{}{n}_0\right|^{2/3}`$. Conversely, when the frequency $`\omega `$ tends to its critical value $`\omega _0`$ for a given direction of incidence $`\stackrel{}{n}_0`$, the saturation value of the AFM amplitude diverges as $`\left|\omega \omega _0\right|^{2/3}`$. In the real situation, of course, the AFM amplitude will be limited by such physical factors as: (i) nonlinear effects, (ii) electromagnetic losses, (iii) structural imperfections of the periodic array, (iv) finiteness of the slab dimensions, (v) deviation of the incident radiation from a perfect plane monochromatic wave.

Fig. 6 gives a good qualitative picture of what really happens in the vicinity of the AFM regime.

Consider a wide monochromatic beam of frequency $`\omega `$ incident on the surface of semi-infinite photonic slab. The direction of incidence $`\stackrel{}{n}_0\stackrel{}{S}_I`$ is chosen so that the condition (3) of the AFM regime is satisfied at $`\omega =\omega _0`$. As frequency $`\omega `$ tends to $`\omega _0`$ from either direction, the normal component $`u_z`$ of the transmitted wave group velocity approaches zero, while the tangential component $`\stackrel{}{u}_{}`$ remains finite

$$u_z\left|\omega \omega _0\right|^{2/3}0,\stackrel{}{u}_{}\stackrel{}{u}_0\text{ as }\omega \omega _0.$$
(5)
This relation together with the equality

$$\frac{\pi }{2}\theta _T=\mathrm{arctan}\frac{u_z}{u_{}}$$
(6)
involving the refraction angle $`\theta _T`$, yield

$$\frac{\pi }{2}\theta _T\left|\omega \omega _0\right|^{2/3}0\text{ as }\omega \omega _0.$$
(7)
Hence, in the vicinity of the AFM regime, the transmitted (refracted) electromagnetic wave can be viewed as a *grazing mode*. The most important and unique feature of this grazing mode directly relates to the fact that the transmittance $`\tau `$ of the semi-infinite slab remains finite even at $`\omega =\omega _0`$ (see, for example, Fig. 4(a)). Indeed, let $`A_I`$ and $`A_T`$ be the cross-section areas of the incident and transmitted (refracted) beams, respectively. Obliviously,

$$\frac{A_T}{A_I}=\frac{\mathrm{cos}\theta _T}{\mathrm{cos}\theta _I}$$
(8)
Let us also introduce the quantities

$$U_I=A_IS_I,U_T=A_TS_T,$$
(9)
where $`S_I`$ and $`S_T`$ are the energy density fluxes of the incident and transmitted waves. $`U_I`$ and $`U_T`$ are the total power transmitted by the incident and transmitted (refracted) beams, respectively. The expressions (8) and (9) imply that

$$\frac{U_T}{U_I}=\frac{S_T\mathrm{cos}\theta _T}{S_I\mathrm{cos}\theta _I}=\frac{\left(S_T\right)_z}{\left(S_I\right)_z}=\tau $$
(10)
which is nothing more than a manifestation of the energy conservation law. Finally, Eq. (10), together with the formula (7), yield

$$S_T=\tau S_I\frac{\mathrm{cos}\theta _I}{\mathrm{cos}\theta _T}\left|\omega \omega _0\right|^{2/3}\mathrm{}\text{ as }\omega \omega _0.$$
(11)
where we have taken into account that $`\tau S_I\mathrm{cos}\theta _I`$ is limited (of the order of magnitude of unity) as $`\omega \omega _0`$. By contrast, in the vicinity of the photonic band edge the transmittance $`\tau `$ of the semi-infinite slab vanishes along with the energy density flux $`S_T`$ of the transmitted (refracted) wave.

The expressions (7) and (11) show that in the vicinity of the AFM regime, the transmitted wave behaves like a grazing mode with huge and nearly tangential energy density flux $`S_T`$ and very small (compared to that of the incident beam) cross-section area $`A_T`$, so that the total power $`U_T=A_TS_T`$ associated with the transmitted wave cannot exceed the total power $`U_I`$ of the incident wave: $`U_T=\tau U_IU_I`$.

The above qualitative consideration is only valid on the scales exceeding the size $`L`$ of the unit cell (which is of the order of magnitude of $`c/\omega `$) and more importantly, exceeding the *transitional distance* $`l=\left(\mathrm{}k_{ev}\right)^1`$ from the slab boundary where the evanescent mode contribution to the resulting electromagnetic field $`\mathrm{\Psi }_T\left(z\right)`$ is still significant. The latter means that the width of both the incident and the refracted beams must be much larger than $`l`$. If the above condition is not met, we cannot treat the transmitted wave as a beam, and the expressions (7) through (11) do not apply. Instead, we would have to use the explicit electrodynamic expressions for $`\mathrm{\Psi }_T\left(z\right)`$, such as the asymptotic formula (122). Note that if the direction $`\stackrel{}{n}`$ of the incident wave propagation and the frequency $`\omega `$ *exactly match* the condition (3) for the AFM regime, the transmitted wave $`\mathrm{\Psi }_T\left(z\right)`$ does not reduce to a superposition (4) of canonical Bloch eigenmodes. Instead, the AFM is described by a general Floquet eigenmode $`\mathrm{\Psi }_{01}\left(z\right)`$ from Eq. (84), which diverges inside the slab as $`z`$, until the nonlinear effects or other limiting factors come into play. The related mathematical analysis is provided in Sections 3 and 4.

In some respects, the remarkable behavior of the AFM, is similar to that of the *frozen mode* related to the phenomenon of *electromagnetic unidirectionality* in nonreciprocal magnetic photonic crystals PRE01 ; PRB03 . In a unidirectional photonic crystal, electromagnetic radiation of a certain frequency $`\omega _0`$ can propagate with finite group velocity $`\stackrel{}{u}z`$ only in one of the two opposite directions, say, from right to left. The problem with the electromagnetic unidirectionality, though, is that it essentially requires the presence of magnetic materials with strong circular birefringence (Faraday rotation) and low losses at the frequency range of interest. Such materials are readily available at the microwave frequencies, but at the infrared and optical frequency ranges, finding appropriate magnetic materials is highly problematic. Thus, at frequencies above $`10^{12}`$ Hz, the electromagnetic unidirectionality along with the respective nonreciprocal magnetic mechanism of the frozen mode formation may prove to be impractical. *By contrast, the occurrence of AFM does not require the presence of magnetic or any other essentially dispersive components in the periodic stack. Therefore, the AFM regime can be realized at any frequencies, including the infrared, optical, and even ultraviolet frequency ranges*. The only essential physical requirement is the presence of anisotropic dielectric layers with proper orientation of the anisotropy axes. An example of such an array is shown in Fig. 2.

In Section 2 we establish the relation between the phenomenon of AFM and the electromagnetic dispersion relation of the periodic layered medium. This allows us to formulate strict and simple symmetry conditions for such a phenomenon to occur, as well as to find out what kind of periodic stratified media can exhibit the effect. Relevant theoretical analysis based on the Maxwell equations in stratified media is carried out in Sections 3 and 4. Finally, in Section 5 we discuss some practical aspects of the phenomenon.

## II Dispersion relation with the AFM

Now we establish the connection between the phenomenon of AFM and the electromagnetic dispersion relation $`\omega \left(\stackrel{}{k}\right),\stackrel{}{k}=(k_x,k_y,k_z)`$ of the periodic stratified medium. In a plane-parallel stratified slab, the tangential components $`(k_x,k_y)`$ of the Bloch wave vector $`\stackrel{}{k}`$ always coincide with those of the incident plane wave in Figs. 1, 3, and 6 while the normal component $`k_z`$ is different from that of the incident wave. To avoid confusion, in further consideration, the $`z`$ component of the Bloch wave vector $`\stackrel{}{k}`$ inside the periodic slab will be denoted as $`k`$ without the subscript $`z`$, namely

$$\text{Inside periodic stack: }\stackrel{}{k}=(k_x,k_y,k).$$

The value of $`k`$ is found by solving the Maxwell equations in the periodic stratified medium for given $`\omega `$ and $`(k_x,k_y)`$; $`k`$ is defined up to a multiple of $`2\pi /L`$, where $`L`$ is the period of the layered structure.

Consider now the frequency $`\omega `$ as function of $`k`$ for fixed $`(k_x,k_y)`$. A typical example of such a dependence is shown in Fig. 7(a). A large gap at the lowest frequencies is determined by the value of the fixed tangential components $`(k_x,k_y)`$ of the quasimomentum $`\stackrel{}{k}`$. This gap vanishes in the case of normal incidence, when $`k_x=k_y=0`$. An alternative and more convenient representation for the dispersion relation is presented in Fig. 7(b), where the plot of $`\omega (k)`$ is obtained for fixed $`(n_x,n_y)`$ based on

$$(n_x,n_y)=(ck_x/\omega ,ck_y/\omega ).$$
(12)
The pair of values $`(n_x,n_y)`$ coincide with the tangential components of the unit vector $`\stackrel{}{n}`$ defining the direction of the incident plane wave propagation. The dependence $`\omega (k)`$ for fixed $`(n_x,n_y)`$ or for fixed $`(k_x,k_y)`$ will be referred to as the *axial dispersion relation*.

Suppose that for $`\stackrel{}{k}=\stackrel{}{k}_0`$and $`\omega =\omega _0=\omega \left(\stackrel{}{k}_0\right)`$, one of the spectral branches $`\omega \left(k\right)`$ develops a stationary inflection point for given $`(k_x,k_y)=(k_{0x},k_{0y})`$, i.e.,

$$\left(\frac{\omega }{k}\right)_{k_x,k_y}|_{\stackrel{}{k}=\stackrel{}{k}_0}=0;\left(\frac{^2\omega }{k^2}\right)_{k_x,k_y}|_{\stackrel{}{k}=\stackrel{}{k}_0}=0;\left(\frac{^3\omega }{k^3}\right)_{k_x,k_y}|_{\stackrel{}{k}=\stackrel{}{k}_0}0,$$
(13)
The value

$$u_z=\left(\frac{\omega }{k}\right)_{k_x,k_y}$$
(14)
in Eq. (13) is the axial component of the group velocity, which vanishes at $`\stackrel{}{k}=\stackrel{}{k}_0`$. Observe that

$$u_x=\left(\frac{\omega }{k_x}\right)_{k,k_y}\text{ and }u_y=\left(\frac{\omega }{k_y}\right)_{k,k_x},$$
(15)
representing the tangential components of the group velocity, may not be zeros at $`\stackrel{}{k}=\stackrel{}{k}_0`$.

Notice that instead of (13), one can use another definition of the stationary inflection point

$$\left(\frac{\omega }{k}\right)_{n_x,n_y}|_{\stackrel{}{k}=\stackrel{}{k}_0}=0,\left(\frac{^2\omega }{k^2}\right)_{n_x,n_y}|_{\stackrel{}{k}=\stackrel{}{k}_0}=0,\left(\frac{^3\omega }{k^3}\right)_{n_x,n_y}|_{\stackrel{}{k}=\stackrel{}{k}_0}0.$$
(16)
The partial derivatives in Eqs. (16) are taken at constant $`(n_x,n_y)`$, rather than at constant $`(k_x,k_y)`$. Observe that the definitions (13) and (16) are equivalent, and we will use both of them.

In Fig. 4(b) we reproduced an enlarged fragment of the upper spectral branch of the axial dispersion relation in Fig. 7(b). For the chosen $`(n_x,n_y)`$, this branch develops a stationary inflection point (16) at $`\omega =\omega _0`$ and $`k=k_0`$. The extended Bloch eigenmode with $`\omega =\omega _0`$ and $`\stackrel{}{k}=\stackrel{}{k}_0`$, associated with the stationary inflection point, turns out to be directly related to the *axially frozen mode* (AFM).

In Sections 3 and 4, based on the Maxwell equations, we prove that the singularity (16) (or, equivalently, (13)) indeed leads to the very distinct AFM regime in the semi-infinite periodic stack. We also show that a necessary condition for such a singularity and, therefore, a necessary condition for the AFM existence is the following property of the axial dispersion relation of the periodic stack

$$\omega (k_x,k_y,k)\omega (k_x,k_y,k)\text{ or, equivalently, }\omega (n_x,n_y,k)\omega (n_x,n_y,k)$$
(17)
This property will be referred to as the *axial spectral asymmetry*. Evidently, the axial dispersion relations presented in Fig. 7, satisfy this criterion. Leaving the proof of the above statements to Section 3, let us look at the constraints imposed by the criterion (17) on the geometry and composition of the periodic stack.

### II.1 Conditions for the axial spectral asymmetry

First of all, notice that a periodic array would definitely have an *axially symmetric* dispersion relation

$$\omega (k_x,k_y,k)=\omega (k_x,k_y,k)\text{ or, equivalently, }\omega (n_x,n_y,k)=\omega (n_x,n_y,k)$$
(18)
if the symmetry group $`G`$ of the periodic stratified medium includes any of the following two symmetry operations

$$m_z,\mathrm{\hspace{0.33em}2}_z^{}=2_z\times R,$$
(19)
where $`m_z`$ is the mirror plane parallel to the layers, $`2_z`$ is the 2-fold rotation about the $`z`$ axis, and $`R`$ is the time reversal operation. Indeed, since $`2_z(k_x,k_y,k)=(k_x,k_y,k)`$ and $`R(k_x,k_y,k)=(k_x,k_y,k)`$, we have

$$2_z^{}(k_x,k_y,k)=(k_x,k_y,k),$$

which implies the relation (18) for arbitrary $`(k_x,k_y)`$. The same is true for the mirror plane $`m_z`$

$$m_z(k_x,k_y,k)=(k_x,k_y,k).$$

Consequently, a necessary condition for the axial spectral asymmetry (17) of a periodic stack is the absence of the symmetry operations (19), i.e.,

$$m_zG\text{ and }2_z^{}G.$$
(20)
In reciprocal (nonmagnetic) media, where by definition, $`RG`$, instead of Eq. (20) one can use the following requirement

$$m_zG\text{ and }2_zG.$$
(21)
Note, that the *axial spectral symmetry* (18) is different from the *bulk spectral symmetry*

$$\omega (k_x,k_y,k)=\omega (k_x,k_y,k)$$
(22)
For example, the space inversion $`I`$ and/or the time reversal $`R`$, if present in $`G`$, ensure the bulk spectral symmetry (22), but neither $`I`$ nor $`R`$ ensures the axial spectral symmetry (18).

#### II.1.1 Application of the criterion (21) to deferent periodic stacks.

The condition (21) for the axial spectral asymmetry imposes certain restrictions on the geometry and composition of the periodic stratified medium, as well as on the direction of the incident wave propagation.

##### Restrictions on the geometry and composition of periodic stack.

First of all, observe that a common periodic stack made up of *isotropic* dielectric components with different refractive indices, always has axially symmetric dispersion relation (18), no matter how complicated the periodic array is or how many different isotropic materials are involved. To prove this, it suffices to note that such a stack always supports the symmetry operation $`2_z`$.

In fact, the symmetry operation $`2_z`$ holds in the more general case when all the layers are either isotropic, or have a purely *in-plane* anisotropy

$$\widehat{\epsilon }=\left[\begin{array}{ccc}\epsilon _{xx}& \epsilon _{xy}& 0\\ \epsilon _{xy}& \epsilon _{yy}& 0\\ 0& 0& \epsilon _{zz}\end{array}\right]$$
(23)
Obviously, the in-plane anisotropy (23) does not remove the symmetry operation $`2_z`$ and, therefore, the property (18) of the axial spectral symmetry holds in this case. Thus, we can state that in order to display the axial spectral asymmetry, the periodic stack must include at least one anisotropic component, either uniaxial or biaxial. In addition, one of the principle axes of the respective dielectric permittivity tensor $`\widehat{\epsilon }`$ must make an oblique angle with the normal to the layers, which means that *at least* *one of the two components* $`\epsilon _{xz}`$ *and* $`\epsilon _{yz}`$ *of the respective dielectric tensor must be nonzero*.

The above requirement gives us a simple and useful idea on what kind of periodic stratified media can support the axial spectral asymmetry and the AFM regime. But this is not a substitute for the stronger symmetry criterion (20) or (21). For example, although the periodic stack in Fig. 8

includes the $`A`$ layers identical to those in Fig. 2, this stack does not meet the criterion (20) for the axial spectral asymmetry. Indeed, the stack in Fig. 8 supports the mirror plane $`m_z`$, which, according to the expression (19), ensures the axial spectral symmetry.

##### Restriction on the direction of incident wave propagation

Consider now an important particular case $`k_x=k_y=0`$ of the normal incidence. The criterion (17) reduces now to the simple requirement

$$\omega \left(\stackrel{}{k}\right)\omega \left(\stackrel{}{k}\right),\text{where }\stackrel{}{k}=(0,0,k),$$
(24)
of the bulk spectral asymmetry, which is prohibited in nonmagnetic photonic crystals due to the time reversal symmetry. Therefore, in the nonmagnetic case, we have the following additional condition for the axial spectral asymmetry

$$k_{}=\sqrt{k_x^2+k_y^2}0,$$
(25)
implying that the AFM cannot be excited in a nonmagnetic semi-infinite stack by a normally incident plane wave, i.e., *the incident angle must be oblique*.

Conditions (21) and (25) may not be necessary in the case of nonreciprocal magnetic stacks (see the details in PRE01 ). But as we mentioned earlier, at frequencies above $`10^{12}`$ Hz, the nonreciprocal effects in common nonconducting materials are negligible. Therefore, in order to have a robust AFM regime in the infrared or optical frequency range, we must satisfy both requirements (21) and (25), regardless of whether or not nonreciprocal magnetic materials are involved.

As soon as the above conditions are met, one can always achieve the AFM regime at any desirable frequency $`\omega `$ within certain frequency range $`\mathrm{\Delta }\omega `$. The frequency range $`\mathrm{\Delta }\omega `$ is determined by the stack geometry and the dielectric materials used, while a specific value of $`\omega `$ within the range can be selected by the direction $`\stackrel{}{n}`$ of the light incidence.

### II.2 Periodic stack with two layers in unit cell

The simplest and the most practical example of a periodic stack supporting the axial spectral asymmetry (17) and, thereby, the AFM regime, is shown in Fig. 2. It is made up of anisotropic $`A`$ layers alternated with isotropic $`B`$ layers. The respective dielectric permittivity tensors are

$$\widehat{\epsilon }_A=\left[\begin{array}{ccc}\epsilon _{xx}& 0& \epsilon _{xz}\\ 0& \epsilon _{yy}& 0\\ \epsilon _{xz}& 0& \epsilon _{zz}\end{array}\right],\widehat{\epsilon }_B=\left[\begin{array}{ccc}\epsilon _B& 0& 0\\ 0& \epsilon _B& 0\\ 0& 0& \epsilon _B\end{array}\right].$$
(26)
For simplicity, we assume

$$\widehat{\mu }_A=\widehat{\mu }_B=\widehat{I}.$$
(27)
The stack in Fig. 2 has the monoclinic symmetry

$$2_y/m_y$$
(28)
with the mirror plane $`m_y`$ normal to the $`y`$ \- axis. Such a symmetry is compatible with the necessary condition (21) for the AFM existence. But as we will see below, the symmetry (28) imposes additional constraints on the direction $`\stackrel{}{n}`$ of the incident wave propagation.

In Fig. 9 we show the axial dispersion relation $`\omega \left(k\right)`$ of this periodic array, computed for four different directions $`(n_x,n_y)`$ of incident wave propagation.

These four cases cover all the possibilities, different in terms of symmetry.

In the case (*a*) of normal incidence, when $`n_x=n_y=0`$, the dispersion relation is axially symmetric, as must be the case with any reciprocal periodic stratified medium (see the explanation after Eq. (24)).

In the case (*b*), when $`n_x=0`$ and $`n_y0`$, the two necessary conditions (21) and (25) for the axial spectral asymmetry are met. Yet, those conditions prove not to be sufficient. Indeed, if $`n_x=0`$, either of the symmetry operations

$$2_y\text{ and }m_y^{}m_y\times R$$
(29)
imposes the relation

$$\omega (0,k_y,k)=\omega (0,k_y,k)$$
(30)
which implies the axial spectral symmetry. Neither stationary inflection point, nor AFM can occur in this case.

In the case (*c*), when $`n_x0`$ and $`n_y=0`$, the situation is more complicated. The quasimomentum $`\stackrel{}{k}`$ lies now in the $`xz`$ plane, which coincides with the mirror plane $`m_y`$. Therefore, every Bloch eigenmode $`\mathrm{\Psi }_\stackrel{}{k}\left(z\right)`$ can be classified as a pure TE or pure TM mode, depending on the $`\mathrm{\Psi }_\stackrel{}{k}\left(z\right)`$ parity with respect to the mirror reflection $`m_y`$

$$\text{for TE mode }m_y\mathrm{\Psi }_\stackrel{}{k}\left(z\right)=\mathrm{\Psi }_\stackrel{}{k}\left(z\right)\text{;   for TM mode: }m_y\mathrm{\Psi }_\stackrel{}{k}\left(z\right)=\mathrm{\Psi }_\stackrel{}{k}\left(z\right).$$
(31)
The TE modes have axially symmetric dispersion relation

$$\text{for TE modes: }\omega (k_x,0,k)=\omega (k_x,0,k).$$
(32)
Indeed, the component $`\epsilon _{xz}`$ of the dielectric tensor $`\widehat{\epsilon }_A`$ does not affect the TE modes, because in this case the electric component $`𝐄(𝐫,t)`$ of the electromagnetic field is parallel to the $`y`$ axis. As a consequence, the axial dispersion relation of the TE spectral branches is similar to that of the isotropic case with $`\epsilon _{xz}=0`$, where it is always symmetric. By contrast, for the TM modes we have $`𝐄(𝐫,t)y`$. Therefore, the TM modes are affected by $`\epsilon _{xz}`$ and display axially asymmetric dispersion relation

$$\text{for TM modes: }\omega (k_x,0,k)\omega (k_x,0,k),$$
(33)
as seen in Fig. 9(c). We wish to remark, though, that the equality (32) cannot be derived from symmetry arguments only. The axial spectral symmetry of the TE modes is not exact and relies on the approximation (27) for the magnetic permeability of the $`A`$ layers. On the other hand, the fact that the spectral branches have different parity (31) with respect to the symmetry operation$`m_y`$, implies that none of the branches can develop a stationary inflection point (see Eq. (87) and explanations thereafter). Thus, in the case $`n_y=0`$, in spite of the axial spectral asymmetry, the AFM regime cannot occur either.

Finally, in the general case (*d*), when $`n_x0`$ and $`n_y0`$, all the spectral branches display the property (17) of the axial spectral asymmetry. In addition, the Bloch eigenmodes now are of the same symmetry (i. e., belong to the same irreducible representation of the wave vector symmetry group) and are neither TE, nor TM. This is exactly the case when the AFM regime can be achieved at some frequencies by proper choice of the incident angle. For instance, if we impose the equality $`n_x=`$ $`n_y`$ and change the incident angle only, it turns out that every single spectral branch at some point develops a stationary inflection point (16) and, thereby, displays the AFM at the respective frequency. If we want the AFM at a specified frequency $`\omega _0`$, then we will have to adjust both $`n_x`$ and $`n_y`$.

## III Electrodynamics of the axially frozen mode

### III.1 Reduced Maxwell equations

We start with the classical Maxwell equations for time-harmonic fields in nonconducting media

$$\times 𝐄\left(\stackrel{}{r}\right)=i\frac{\omega }{c}𝐁\left(\stackrel{}{r}\right),\times 𝐇\left(\stackrel{}{r}\right)=i\frac{\omega }{c}𝐃\left(\stackrel{}{r}\right),$$
(34)
where

$$𝐃\left(\stackrel{}{r}\right)=\widehat{\epsilon }\left(\stackrel{}{r}\right)𝐄\left(\stackrel{}{r}\right),𝐁\left(\stackrel{}{r}\right)=\widehat{\mu }\left(\stackrel{}{r}\right)𝐇\left(\stackrel{}{r}\right),$$
(35)
In a lossless dielectric medium, the material tensors $`\widehat{\epsilon }\left(\stackrel{}{r}\right)`$ and $`\widehat{\mu }\left(\stackrel{}{r}\right)`$ are Hermitian. In a stratified medium, the tensors $`\widehat{\epsilon }\left(\stackrel{}{r}\right)`$ and $`\widehat{\mu }\left(\stackrel{}{r}\right)`$ depend on a single Cartesian coordinate $`z`$, and the Maxwell equations (34) can be recast as

$$\times 𝐄\left(\stackrel{}{r}\right)=i\frac{\omega }{c}\widehat{\mu }\left(z\right)𝐇\left(\stackrel{}{r}\right),\times 𝐇\left(\stackrel{}{r}\right)=i\frac{\omega }{c}\widehat{\epsilon }\left(z\right)𝐄\left(\stackrel{}{r}\right).$$
(36)
Solutions for Eq. (36) are sought in the following form

$$𝐄\left(\stackrel{}{r}\right)=e^{i\left(k_xx+k_yy\right)}\stackrel{}{E}\left(z\right),𝐇\left(\stackrel{}{r}\right)=e^{i\left(k_xx+k_yy\right)}\stackrel{}{H}\left(z\right).$$
(37)
The substitution (37) transforms the system of six linear equation (36) into a system of four linear differential equations

$$_z\mathrm{\Psi }\left(z\right)=i\frac{\omega }{c}M\left(z\right)\mathrm{\Psi }\left(z\right),\mathrm{\Psi }\left(z\right)=\left[\begin{array}{c}E_x\left(z\right)\\ E_y\left(z\right)\\ H_x\left(z\right)\\ H_y\left(z\right)\end{array}\right]$$
(38)
The explicit expression for the Maxwell operator $`M\left(z\right)`$ is

$$M\left(z\right)=\left[\begin{array}{cc}M_{11}& M_{12}\\ M_{21}& M_{22}\end{array}\right]$$
(39)
where

$`M_{11}`$ $`=`$ $`\left[\begin{array}{cc}\frac{\epsilon _{xz}^{}}{\epsilon _{zz}}n_x\frac{\mu _{yz}}{\mu _{zz}}n_y& \left(\frac{\epsilon _{yz}^{}}{\epsilon _{zz}}+\frac{\mu _{yz}}{\mu _{zz}}\right)n_x\\ \left(\frac{\epsilon _{xz}^{}}{\epsilon _{zz}}\frac{\mu _{xz}}{\mu _{zz}}\right)n_y& \frac{\epsilon _{yz}^{}}{\epsilon _{zz}}n_y\frac{\mu _{xz}}{\mu _{zz}}n_x\end{array}\right],`$
$`M_{22}`$ $`=`$ $`\left[\begin{array}{cc}\frac{\epsilon _{yz}}{\epsilon _{zz}}n_y\frac{\mu _{xz}^{}}{\mu _{zz}}n_x& \left(\frac{\epsilon _{yz}}{\epsilon _{zz}}\frac{\mu _{yz}^{}}{\mu _{zz}}\right)n_x\\ \left(\frac{\epsilon _{xz}}{\epsilon _{zz}}\frac{\mu _{xz}^{}}{\mu _{zz}}\right)n_y& \frac{\epsilon _{xz}}{\epsilon _{zz}}n_x\frac{\mu _{yz}^{}}{\mu _{zz}}n_y\end{array}\right],`$
$`M_{12}`$ $`=`$ $`\left[\begin{array}{cc}\mu _{xy}^{}\frac{\mu _{xz}^{}\mu _{yz}}{\mu _{zz}}+\frac{n_xn_y}{\epsilon _{zz}}& \mu _{yy}\frac{\mu _{yz}\mu _{yz}^{}}{\mu _{zz}}\frac{n_x^2}{\epsilon _{zz}}\\ \mu _{xx}+\frac{\mu _{xz}\mu _{xz}^{}}{\mu _{zz}}+\frac{n_y^2}{\epsilon _{zz}}& \mu _{xy}+\frac{\mu _{xz}\mu _{yz}^{}}{\mu _{zz}}\frac{n_xn_y}{\epsilon _{zz}}\end{array}\right],`$
$`M_{21}`$ $`=`$ $`\left[\begin{array}{cc}\epsilon _{xy}^{}+\frac{\epsilon _{xz}^{}\epsilon _{yz}}{\epsilon _{zz}}\frac{n_xn_y}{\mu _{zz}}& \epsilon _{yy}+\frac{\epsilon _{yz}\epsilon _{yz}^{}}{\epsilon _{zz}}+\frac{n_x^2}{\mu _{zz}}\\ \epsilon _{xx}\frac{\epsilon _{xz}\epsilon _{xz}^{}}{\epsilon _{zz}}\frac{n_y^2}{\mu _{zz}}& \epsilon _{xy}\frac{\epsilon _{xz}\epsilon _{yz}^{}}{\epsilon _{zz}}+\frac{n_xn_y}{\mu _{zz}}\end{array}\right].`$
The Cartesian components of the material tensors $`\widehat{\epsilon }`$ and $`\widehat{\mu }`$ are functions of $`z`$ and (in dispersive media) $`\omega `$. The reduced Maxwell equation (38) should be complemented with the following expressions for the $`z`$ components of the fields

$$\begin{array}{c}E_z=\left(n_xH_y+n_yH_x\epsilon _{13}^{}E_x\epsilon _{23}^{}E_y\right)\epsilon _{zz}^1\\ H_z=\left(n_xE_yn_yE_x\mu _{13}^{}H_x\mu _{23}^{}H_y\right)\mu _{zz}^1\end{array}$$
(44)
where $`(n_x,n_y)`$ are defined in Eq. (12).

Notice that in the case of normal incidence, the Maxwell operator is drastically simplified

$$\text{ }M_{11}=M_{22}=0,\text{for }n_x=n_y=0\text{.}$$
(45)
This is the case we dealt with in PRB03  when considering the phenomenon of electromagnetic unidirectionality in nonreciprocal magnetic photonic crystals. By contrast, the objective of this Section is to show how the terms $`M_{11}`$ and $`M_{22}`$, occurring only in the case of oblique incidence, can lead to the phenomenon of AFM, regardless of whether or not the nonreciprocal effects are present. Note that $`M_{11}`$ and $`M_{22}`$ are also nonzero in materials with linear magnetoelectric effect (see, for example, Ref. 15 and references therein), but we are not considering here such an exotic situation.

Importantly, the $`4\times 4`$ matrix $`M\left(z\right)`$ in Eq. (39) has the property of $`J`$ \- Hermitivity defined as

$$\left(JM\right)^{}=JM$$
(46)
where

$$J=J^1=\left[\begin{array}{cccc}0& 0& 0& 1\\ 0& 0& 1& 0\\ 0& 1& 0& 0\\ 1& 0& 0& 0\end{array}\right]$$
(47)
Different versions of the reduced Maxwell equation (38) can be found in the extensive literature on electrodynamics of stratified media (see, for example, Tmatrix ; Abdul00 ; Abdul99 , and references therein). For more detailed studies of $`J`$ \- Hermitian and $`J`$ \- unitary operators see JHerm .

### III.2 The transfer matrix

The Cauchy problem

$$_z\mathrm{\Psi }\left(z\right)=i\frac{\omega }{c}M\left(z\right)\mathrm{\Psi }\left(z\right),\mathrm{\Psi }\left(z_0\right)=\mathrm{\Psi }_0$$
(48)
for the reduced Maxwell equation (38) has a unique solution

$$\mathrm{\Psi }\left(z\right)=T(z,z_0)\mathrm{\Psi }\left(z_0\right)$$
(49)
where the $`4\times 4`$ matrix $`T(z,z_0)`$ is so-called *transfer matrix*. From the definition (49) of the transfer matrix it follows that

$$T(z,z_0)=T(z,z^{})T(z^{},z_0),T(z,z_0)=T^1(z_0,z),T(z,z)=I.$$
(50)
The matrix $`T(z,z_0)`$ is uniquely defined by the following Cauchy problem

$$_zT(z,z_0)=i\frac{\omega }{c}M\left(z\right)T(z,z_0),T(z,z)=I.$$
(51)
The equation (51), together with $`J`$ \- Hermitivity (46) of the Maxwell operator $`M\left(z\right)`$, imply that the matrix $`T(z,z_0)`$ is $`J`$ \- *unitarity*, i.e.,

$$T^{}(z,z_0)=JT^1(z,z_0)J.$$
(52)
(see the proof in Appendix 1). The $`J`$ \- unitarity (52) of the transfer matrix imposes strong constraints on its eigenvalues (see Eq. (65)). It also implies that

$$\left|detT(z,z_0)\right|=1.$$
(53)
The transfer matrix $`T_S`$ of a stack of layers is a sequential product of the transfer matrices $`T_m`$ of the constitutive layers

$$T_S=\underset{m}{}T_m$$
(54)
If the individual layers are homogeneous, the corresponding single-layer transfer matrices $`T_m`$ are explicitly expressed in terms of the respective Maxwell operators $`M_m`$

$$T_m=\mathrm{exp}\left(iz_mM_m\right)$$
(55)
where $`z_m`$ is the thickness of the $`m`$-th layer. The explicit expression for $`M_m`$ is given by (39). Thus, formula (54), together with (55) and (39), gives us an explicit expression for the transfer matrix $`T_S`$ of an arbitrary stack of anisotropic dielectric layers. $`T_S`$ is a function of (i) the material tensors $`\widehat{\epsilon }`$ and $`\widehat{\mu }`$ in each layer of the stack, (ii) the layer thicknesses, (iii) the frequency $`\omega `$, and (iv) the tangential components $`(k_x,k_y)=(n_x\omega /c,n_y\omega /c)`$ of the wave vector.

Consider the important particular case of normal wave propagation. Using Eq. (55) and the explicit expression (39) for the Maxwell operator, one can prove that

$$det\left(T_S\right)=1,\text{for }n_x=n_y=0.$$
(56)
Additional information related to the transfer matrix formalism can be found in Tmatrix ; Abdul00 ; Abdul99  and references therein.

### III.3 Periodic arrays. Bloch eigenmodes.

In a periodic layered structure, all material tensors, along with the $`J`$ \- Hermitian matrix $`M(z)`$ in Eq. (38), are periodic functions of $`z`$

$$M\left(z+L\right)=M\left(z\right)$$
(57)
where $`L`$ is the length of a primitive cell of the periodic stack. By definition, Bloch solutions $`\mathrm{\Psi }_k\left(z\right)`$ of the reduced Maxwell equation (38) with the periodic operator $`M(z)`$ satisfy

$$\mathrm{\Psi }_k\left(z+L\right)=e^{ikL}\mathrm{\Psi }_k\left(z\right)$$
(58)
The definition (49) of the $`T`$ \- matrix together with Eq. (58) give

$$\mathrm{\Psi }_k\left(z+L\right)=T(z+L,z)\mathrm{\Psi }_k\left(z\right)=e^{ikL}\mathrm{\Psi }_k\left(z\right).$$
(59)
Introducing the transfer matrix of a primitive cell

$$T_L=T(L,0)$$
(60)
we have from Eq. (59)

$$T_L\mathrm{\Phi }_k=e^{ikL}\mathrm{\Phi }_k,\text{ where }\mathrm{\Phi }_k=\mathrm{\Psi }_k\left(0\right).$$
(61)
Thus, the eigenvectors of the transfer matrix $`T_L`$ of the unit cell are uniquely related to the eigenmodes of the reduced Maxwell equation (38) through the relations

$$\mathrm{\Phi }_{k_1}=\mathrm{\Psi }_{k_1}\left(0\right),\mathrm{\Phi }_{k_2}=\mathrm{\Psi }_{k_2}\left(0\right),\mathrm{\Phi }_{k_3}=\mathrm{\Psi }_{k_3}\left(0\right),\mathrm{\Phi }_{k_4}=\mathrm{\Psi }_{k_4}\left(0\right)$$
(62)
The respective four eigenvalues

$$X_i=e^{ik_iL},i=1,2,3,4$$
(63)
of $`T_L`$ are the roots of the characteristic equation

$$F\left(X\right)=0,\text{ where }F\left(X\right)=det\left(T_LX\widehat{I}\right)=X^4+P_3X^3+P_2X^2+P_1X+1.$$
(64)
For any given $`\omega `$ and $`(k_x,k_y)`$, the characteristic equation defines a set of four values $`\{X_1,X_2,X_3,X_4\}`$, or equivalently, $`\{k_1,k_2,k_3,k_4\}`$. Real $`k`$ correspond to propagating Bloch waves (extended modes), while complex $`k`$ correspond to evanescent modes. Evanescent modes are relevant near photonic crystal boundaries and other structural irregularities.

The $`J`$-unitarity (52) of $`T_L`$ imposes the following restriction on the eigenvalues (63) for any given $`\omega `$ and $`(k_x,k_y)`$

$$\{k_i\}\{k_i^{}\},i=1,2,3,4.$$
(65)
In view of the relation (65), one has to consider three different situation. The first possibility

$$k_1k_1^{},k_2k_2^{},k_3k_3^{},k_4k_4^{}$$
(66)
relates to the case of all four Bloch eigenmodes being extended. The second possibility

$$k_1=k_1^{},k_2=k_2^{},k_3=k_4^{},\text{where }k_3k_3^{},k_4k_4^{},$$
(67)
relates to the case of two extended and two evanescent modes. The last possibility

$$k_1=k_2^{},k_3=k_4^{},\text{where }k_1k_1^{},k_2k_2^{},k_3k_3^{},k_4k_4^{}$$
(68)
relates the case of a frequency gap, when all four Bloch eigenmodes are evanescent.

Observe that the relation

$$k_1+k_2+k_3+k_40$$

valid in the case of normal incidence (see Refs. PRE01 ; PRB03 ), may not apply now.

#### III.3.1 Axial spectral symmetry

Assume that the transfer matrix $`T_L`$ is similar to its inverse

$$T_L=U^1T_L^1U$$
(69)
where $`U`$ is an invertible $`4\times 4`$ matrix. This assumption together with the property (52) of $`J`$-unitarity, imply the similarity of $`T_L`$ and $`T_L^{}`$

$$T_L=V^1T_L^{}V\text{,  where }V=JU.$$
(70)
This relation imposes additional restrictions on the eigenvalues (63) for a given frequency $`\omega `$ and given $`(k_x,k_y)`$

$$\{k_i\}\{k_i\},i=1,2,3,4.$$
(71)
The relation (71) is referred to as the *axial spectral symmetry*, because in terms of the corresponding axial dispersion relation, it implies the equality (18) for every spectral branch.

If the sufficient condition (69) for the axial spectral symmetry is not in place, then we can have for a given $`\omega `$ and $`(k_x,k_y)`$

$$\{k_i\}\{k_i\},i=1,2,3,4$$
(72)
which implies the *axial spectral asymmetry* (17).

### III.4 Stationary inflection point

The coefficients of the characteristic polynomial $`F\left(X\right)`$ in Eq. (64) are functions of $`\omega `$ and $`(k_x,k_y)`$. Let $`F_0\left(X\right)`$ be the characteristic polynomial at the stationary inflection point (16), where $`\omega =\omega _0`$ and $`(k_x,k_y)=(k_{0x},k_{0y})`$. The stationary inflection point (16) can also be defined as follows

$$F_0\left(X\right)=0,F_0^{}\left(X\right)=0,F_0^{\prime \prime }\left(X\right)=0,F_0^{\prime \prime \prime }\left(X\right)0.$$
(73)
This relation requires the respective value of $`X_0=\mathrm{exp}\left(ik_0L\right)`$ to be a triple root of the characteristic polynomial $`F_0\left(X\right)`$ implying

$$F_0\left(X\right)=\left(XX_1\right)\left(XX_0\right)^3=0.$$
(74)
A small deviation of the frequency $`\omega `$ from its critical value $`\omega _0`$ changes the coefficients of the characteristic polynomial and removes the triple degeneracy of the solution $`X_0`$

$$XX_06^{1/3}\left(\frac{F_0/\omega }{^3F_0/X^3}\right)^{1/3}\left(\omega \omega _0\right)^{1/3}\xi ,\xi =1,e^{2\pi i/3},e^{2\pi i/3}.$$
(75)
or, in terms of the axial quasimomentum $`k`$

$$kk_06^{1/3}\left(\frac{\omega \omega _0}{\omega _0^{\prime \prime \prime }}\right)^{1/3}\xi ,\xi =1,e^{2\pi i/3},e^{2\pi i/3}$$
(76)
where

$$\omega _0^{\prime \prime \prime }=\left(\frac{^3\omega }{k^3}\right)_{k_x,k_y}|_{\stackrel{}{k}=\stackrel{}{k}_0}>0.$$
(77)
The three solutions (76) can also be rearranged as

$$\{\begin{array}{c}k_{ex}k_0+6^{1/3}\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left(\omega \omega _0\right)^{1/3},\\ k_{ev}k_0+\frac{1}{2}\left(6\right)^{1/3}\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left(\omega \omega _0\right)^{1/3}+i\frac{\sqrt{3}}{2}6^{1/3}\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left|\omega \omega _0\right|^{1/3},\\ k_{EV}k_0+\frac{1}{2}\left(6\right)^{1/3}\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left(\omega \omega _0\right)^{1/3}i\frac{\sqrt{3}}{2}6^{1/3}\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left|\omega \omega _0\right|^{1/3}.\end{array}$$
(78)
The real $`k_{ex}`$ in (78) relates to the extended mode $`\mathrm{\Psi }_{ex}\left(z\right)`$, with $`u_z=0`$ at $`\omega =\omega _0`$. The other two solutions, $`k_{ev}`$ and $`k_{EV}=k_{ev}^{}`$, correspond to a pair of evanescent modes $`\mathrm{\Psi }_{ev}\left(z\right)`$ and $`\mathrm{\Psi }_{EV}\left(z\right)`$ with positive and negative infinitesimally small imaginary parts, respectively. Those modes are truly evanescent (i.e., have $`\mathrm{}k0`$) only if $`\omega \omega _0`$, but it does not mean that at $`\omega =\omega _0`$, the eigenmodes $`\mathrm{\Psi }_{ev}\left(z\right)`$ and $`\mathrm{\Psi }_{EV}\left(z\right)`$ become extended. In what follows we will take a closer look at this problem.

#### III.4.1 Eigenmodes at the frequency of AFM

Consider the vicinity of stationary inflection point (13). As long as $`\omega \omega _0`$, the four eigenvectors (62) of the transfer matrix $`T_L`$ comprise two extended and two evanescent Bloch solutions. One of the extended modes (say, $`\mathrm{\Phi }_{k_1}`$) corresponds to the non-degenerate real root $`X_1=e^{ik_1L}`$ of the characteristic equation (64). This mode has negative axial group velocity $`u_z\left(k_1\right)<0`$ and, therefor, is of no interest for us. The other three eigenvectors of $`T_L`$ correspond to three nearly degenerate roots (75). As $`\omega `$ approaches $`\omega _0`$, these three eigenvalues become degenerate, while the respective three eigenvectors $`\mathrm{\Phi }_{k_2},\mathrm{\Phi }_{k_3},`$ and $`\mathrm{\Phi }_{k_4}`$ become collinear

$$\mathrm{\Phi }_{k_2}\alpha _1\mathrm{\Phi }_{k_0},\mathrm{\Phi }_{k_3}\alpha _2\mathrm{\Phi }_{k_0},\mathrm{\Phi }_{k_4}\alpha _3\mathrm{\Phi }_{k_0},\text{as }\omega \omega _0.$$
(79)
The latter important feature relates to the fact that at $`\omega =\omega _0`$, the matrix $`T_L`$has a nontrivial Jordan canonical form

$$U^1T_LU=\left[\begin{array}{cccc}X_1& 0& 0& 0\\ 0& X_0& 1& 0\\ 0& 0& X_0& 1\\ 0& 0& 0& X_0\end{array}\right],\text{at }\omega =\omega _0$$
(80)
and, therefore, cannot be diagonalized. It is shown rigorously in PRB03 , that the very fact that the $`T_L`$ eigenvalues display the singularity (75), implies that at $`\omega =\omega _0`$, the matrix $`T_L`$ has the canonical form (80). In line with (79), the matrix $`T_L`$ from Eq. (80) has only two (not four!) eigenvectors:

1. $`\mathrm{\Phi }_{k_1}=\mathrm{\Psi }_{k_1}\left(0\right)`$, corresponding to the non-degenerate root $`X_1`$ and relating to the extended mode with $`u_z<0`$;
2. $`\mathrm{\Phi }_{k_0}=\mathrm{\Psi }_{k_0}\left(0\right)`$, corresponding to the triple root $`X_0`$ and related to the AFM.

The other two solutions of the Maxwell equation (38) at $`\omega =\omega _0`$ are general Floquet eigenmodes, which do not reduce to the canonical Bloch form (58). Yet, they can be related to $`\mathrm{\Psi }_{k_0}\left(z\right)`$. Indeed, following the standard procedure (see, for example, Linalg ; LODE ), consider an extended Bloch solution $`\mathrm{\Psi }_k\left(z\right)`$ of the reduced Maxwell equation (38)

$$𝐋\mathrm{\Psi }_k\left(z\right)=0,\text{where}𝐋=_zi\frac{\omega }{c}M\left(z\right)$$
(81)
where both operators $`M\left(z\right)`$ and $`𝐋\left(z\right)`$ are functions of $`\omega `$ and $`(k_x,k_y)`$. Assume now that the axial dispersion relation $`\omega \left(k\right)`$ has a stationary inflection point (13) at $`k=k_0`$. Differentiating Eq. (81) with respect to $`k`$ at constant $`(k_x,k_y)`$ gives, with consideration for Eq. (13),

$$𝐋_k\mathrm{\Psi }_k\left(z\right)=0,𝐋_{kk}^2\mathrm{\Psi }_k\left(z\right)=0,\text{at }k=k_0.$$

This implies that at $`k=k_0`$, both functions

$$\mathrm{\Psi }_{01}\left(z\right)=_k\mathrm{\Psi }_k\left(z\right)|_{k=k_0},\text{and}\mathrm{\Psi }_{02}\left(z\right)=_{kk}^2\mathrm{\Psi }_k\left(z\right)|_{k=k_0}$$
(82)
are also eigenmodes of the reduced Maxwell equation at $`\omega =\omega _0`$. Representing $`\mathrm{\Psi }_k\left(z\right)`$ in the form

$$\mathrm{\Psi }_k\left(z\right)=\psi _k\left(z\right)e^{ikz},\text{where }\psi _k\left(z+L\right)=\psi _k\left(L\right),\mathrm{}k=0$$
(83)
and substituting Eq. (83) into (82) we get

$`\mathrm{\Psi }_{01}\left(z\right)`$ $`=`$ $`\overline{\mathrm{\Psi }}_{k_0}\left(z\right)+iz\mathrm{\Psi }_{k_0}\left(z\right),`$ (84)
$`\mathrm{\Psi }_{02}\left(z\right)`$ $`=`$ $`\overline{\mathrm{\Psi }}_{k_0}^{}\left(z\right)+iz\overline{\mathrm{\Psi }}_{k_0}\left(z\right)z^2\mathrm{\Psi }_{k_0}\left(z\right),`$ (85)
where

$$\overline{\mathrm{\Psi }}_{k_0}\left(z\right)=\left(_k\psi _k\left(z\right)\right)_{k=k_0}e^{ik_0z}\text{and }\overline{\mathrm{\Psi }}_{k_0}^{}\left(z\right)=\left(_{kk}^2\psi _k\left(z\right)\right)_{k=k_0}e^{ik_0z}$$

are auxiliary Bloch functions (not eigenmodes).

To summarize, at the frequency $`\omega _0`$ of AFM, there are four solutions for the reduced Maxwell equation (38)

$$\mathrm{\Psi }_{k_1}\left(z\right),\mathrm{\Psi }_{k_0}\left(z\right),\mathrm{\Psi }_{01}\left(z\right),\mathrm{\Psi }_{02}\left(z\right)$$
(86)
The first two solutions from (86) are extended Bloch eigenmodes with $`u_z<0`$ and $`u_z=0`$, respectively. The other two solutions diverges as the first and the second power of $`z`$, respectively, they are referred to as general (non-Bloch) Floquet modes.

Deviation of the frequency $`\omega `$ from $`\omega _0`$ removes the triple degeneracy (80) of the matrix $`T_L`$, as seen from Eq. (75). The modified matrix $`T_L`$ can now be reduced to a diagonal form with the set (62) of four eigenvectors comprising two extended and two evanescent Bloch solutions.

### III.5 Symmetry considerations

In Section 2, we discussed the relation between the symmetry of the axial dispersion relation of a periodic stack, and the phenomenon of AFM. At this point we can prove that indeed, the axial spectral asymmetry (17) is a necessary condition for the occurrence of the stationary inflection point and for the AFM associated with such a point. As we have seen earlier in this Section, the stationary inflection point relates to a triple root of the characteristic polynomial $`F(X)`$ from Eq. (64). Since $`F(X)`$ is a polynomial of the fourth degree, it cannot have a symmetric pair of triple roots, that would have been the case for axially symmetric dispersion relation. Hence, only asymmetric axial dispersion relation $`\omega \left(k\right)`$ can display a stationary inflection point (13) or, equivalently, (16), as shown in Fig. 4(b). In this respect, the situation with the AFM is somewhat similar to that of the frozen mode in unidirectional magnetic photonic crystals PRB03 . The difference lies in the physical nature of the phenomenon. The bulk spectral asymmetry (24) leading to the effect of electromagnetic unidirectionality, essentially requires the presence of nonreciprocal magnetic materials. By contrast, the axial spectral asymmetry (17) along with the AFM regime can be realized in perfectly reciprocal periodic dielectric stacks with symmetric bulk dispersion relation (22). On the other hand, the axial spectral asymmetry essentially requires an oblique light incidence, which is not needed for the bulk spectral asymmetry.

Another important symmetry consideration is that in the vicinity of the stationary inflection point (13), all four Bloch eigenmodes (62) must have the same symmetry, which means that all of them *must belong to the same one-dimensional irreducible representation of the Bloch wave vector group*. This condition is certainly met when the direction defined by $`(n_x,n_y)`$ is not special in terms of symmetry. Let us see what happens if the above condition is not in place. Consider the situation (31), when at any given frequency $`\omega `$ and fixed $`(n_x,n_y)=(n_x,0)`$, two of the Bloch eigenmodes are TE modes and the other two are TM modes. Note, that TE and TM modes belong to *different* one-dimensional representations of the Bloch wave vector group. In such a case, the transfer matrix $`T_L`$ can be reduced to the block-diagonal form

$$T_L=\left[\begin{array}{cccc}T_{11}& T_{12}& 0& 0\\ T_{21}& T_{22}& 0& 0\\ 0& 0& T_{33}& T_{34}\\ 0& 0& T_{43}& T_{44}\end{array}\right]$$

The respective characteristic polynomial $`F(X)`$ degenerates into

$$F(X)=F_{TE}(X)F_{TM}(X)$$
(87)
where $`F_{TE}(X)`$ and $`F_{TM}(X)`$ are independent second degree polynomials describing the TE and TM spectral branches, respectively. Obviously, in such a situation, the transfer matrix cannot have the nontrivial canonical form Eq. (80), and the respective axial dispersion relation cannot develop a stationary inflection point (73), regardless of whether or not the axial spectral asymmetry is in place.

## IV The AFM regime in a semi-infinite stack

### IV.1 Boundary conditions

In vacuum (to the left of semi-infinite slab in Fig. 1) the electromagnetic field $`\mathrm{\Psi }_V\left(z\right)`$ is a superposition of the incident and reflected waves

$$\mathrm{\Psi }_V\left(z\right)=\mathrm{\Psi }_I\left(z\right)+\mathrm{\Psi }_R\left(z\right),\text{ at }z<0$$
(88)
At the slab boundary we have

$$\mathrm{\Psi }_V\left(0\right)=\mathrm{\Psi }_I\left(0\right)+\mathrm{\Psi }_R\left(0\right)=\mathrm{\Phi }_I+\mathrm{\Phi }_R.$$
(89)
where

$`\mathrm{\Phi }_I`$ $`=`$ $`\left[\begin{array}{c}E_{I,x}\\ E_{I,y}\\ H_{I,x}\\ H_{I,y}\end{array}\right]=\left[\begin{array}{c}E_{I,x}\\ E_{I,y}\\ E_{I,x}n_xn_yn_z^1E_{I,y}\left(1n_x^2\right)n_z^1\\ E_{I,x}\left(1n_y^2\right)n_z^1+E_{I,y}n_xn_yn_z^1\end{array}\right],`$ (98)
$`\mathrm{\Phi }_R`$ $`=`$ $`\left[\begin{array}{c}E_{R,x}\\ E_{R,y}\\ H_{R,x}\\ H_{R,y}\end{array}\right]=\left[\begin{array}{c}E_{R,x}\\ E_{R,y}\\ E_{R,x}n_xn_yn_z^1+E_{R,y}\left(1n_x^2\right)n_z^1\\ E_{R,x}\left(1n_y^2\right)n_z^1E_{R,y}n_xn_yn_z^1\end{array}\right].`$ (107)
The complex vectors $`\stackrel{}{E}_I,\stackrel{}{H}_I`$ and $`\stackrel{}{E}_R,\stackrel{}{H}_R`$ are related to the actual electromagnetic field components $`𝐄_I,𝐇_I`$ and $`𝐄_R,𝐇_R`$ as

$`𝐄_I`$ $`=`$ $`e^{i\frac{\omega }{c}\left(n_xx+n_yy\right)}\stackrel{}{E}_I\left(z\right),𝐇_I=e^{i\frac{\omega }{c}\left(n_xx+n_yy\right)}\stackrel{}{H}_I,`$
$`𝐄_R`$ $`=`$ $`e^{i\frac{\omega }{c}\left(n_xx+n_yy\right)}\stackrel{}{E}_R\left(z\right),𝐇_R=e^{i\frac{\omega }{c}\left(n_xx+n_yy\right)}\stackrel{}{H}_R.`$
The transmitted wave $`\mathrm{\Psi }_T\left(z\right)`$ inside the semi-infinite slab is a superposition of two Bloch eigenmodes

$$\mathrm{\Psi }_T\left(z\right)=\mathrm{\Psi }_1\left(z\right)+\mathrm{\Psi }_2\left(z\right),\text{ at }z>0.$$
(108)
(in the case of a finite slab, all four eigenmodes (62) would contribute to $`\mathrm{\Psi }_T\left(z\right)`$). The eigenmodes $`\mathrm{\Psi }_1\left(z\right)`$and $`\mathrm{\Psi }_2\left(z\right)`$ can be both extended (with $`u_x>0`$), one extended and one evanescent (with $`u_x>0`$ and $`\mathrm{}k>0,`$ respectively), or both evanescent (with $`\mathrm{}k>0`$), depending on which of the three cases (66), (67), or (68) we are dealing with. In particular, in the vicinity of the AFM (e.g., the vicinity of $`\omega _0`$ in Fig. 4(b)), we always have the situation (67). Therefore, in the vicinity of AFM, $`\mathrm{\Psi }_T\left(z\right)`$ is a superposition of the extended eigenmode $`\mathrm{\Psi }_{ex}\left(z\right)`$ with the group velocity $`u_z>0`$, and the evanescent mode $`\mathrm{\Psi }_{ev}\left(z\right)`$ with $`\mathrm{}k>0`$

$$\mathrm{\Psi }_T\left(z\right)=\mathrm{\Psi }_{ex}\left(z\right)+\mathrm{\Psi }_{ev}\left(z\right),\text{ at }z>0.$$
(109)
The asymptotic expressions for the respective wave vectors $`k_{ex}`$ and $`k_{ev}`$ in the vicinity of AFM are given in Eq. (78).

When the frequency $`\omega `$ exactly coincides with the frequency $`\omega _0`$ of the AFM, the representation (109) for $`\mathrm{\Psi }_T\left(z\right)`$ is not valid. In such a case, according to Eq. (78), there is no evanescent modes at all. It turns out that at $`\omega =\omega _0`$, the electromagnetic field inside the slab is a superposition of the extended mode $`\mathrm{\Psi }_{k_0}\left(z\right)`$ and the (non-Bloch) Floquet eigenmode $`\mathrm{\Psi }_{01}\left(z\right)`$ from Eq. (84)

$$\mathrm{\Psi }_T\left(z\right)=\mathrm{\Psi }_{k_0}\left(z\right)+\mathrm{\Psi }_{01}\left(z\right),\text{at }\omega =\omega _0\text{ and }z>0.$$
(110)
Since the extended eigenmode $`\mathrm{\Psi }_{k_0}\left(z\right)`$ has zero axial group velocity $`u_z`$, it does not contribute to the axial energy flux $`S_z`$. By contrast, the divergent non-Bloch contribution $`\mathrm{\Psi }_{01}\left(z\right)`$ is associated with the finite axial energy flux $`S_z>0`$, although the notion of group velocity does not apply here. The detailed analysis is carried out in the next subsection.

Knowing the eigenmodes inside the slab and using the standard electromagnetic boundary conditions

$$\mathrm{\Phi }_T=\mathrm{\Phi }_I+\mathrm{\Phi }_R\text{}\text{where }\mathrm{\Phi }=\mathrm{\Psi }\left(0\right),$$
(111)
one can express the amplitude and composition of the transmitted wave $`\mathrm{\Psi }_T`$ and reflected wave $`\mathrm{\Psi }_R`$, in terms of the amplitude and polarization of the incident wave $`\mathrm{\Psi }_I`$. This gives us the transmittance and reflectance coefficients (2) of the semi-infinite slab, as well as the electromagnetic field distribution $`\mathrm{\Psi }_T\left(z\right)`$ inside the slab, as functions of the incident wave polarization, the direction $`\stackrel{}{n}`$ of incidence, and the frequency $`\omega `$.

### IV.2 Field amplitude inside semi-infinite slab

In what follows we assume that $`\omega `$ can be arbitrarily close but not equal to $`\omega _0`$, unless otherwise is explicitly stated. This will allow us to treat the transmitted wave $`\mathrm{\Psi }_T\left(z\right)`$ as a superposition (109) of one extended and one evanescent mode. Since evanescent modes do not transfer energy in the $`z`$ direction, the extended mode is solely responsible for the axial energy flux $`S_z`$

$$S_z\left(\mathrm{\Psi }_T\right)=S_z\left(\mathrm{\Psi }_{ex}\right).$$
(112)
According to Eq. (130), $`S_z`$ does not depend on $`z`$ and can be expressed in terms of the semi-infinite slab transmittance $`\tau `$ from Eq. (2)

$$S_z=\tau \left(\stackrel{}{S}_I\right)_z=\tau S_I,$$
(113)
where $`S_I=\left(\stackrel{}{S}_I\right)_z`$ is the axial energy flux of the incident wave, which is set to be unity.

The energy density $`W_{ex}`$ associated with the extended mode $`\mathrm{\Psi }_{ex}\left(z\right)`$ can be expressed in terms of the axial component $`u_z`$ of its group velocity and the axial component $`S_z\left(\mathrm{\Psi }_{ex}\right)`$ of the respective energy density flux

$$W_{ex}=u_z^1S_z\left(\mathrm{\Psi }_{ex}\right)=\left(\frac{\omega }{k}\right)_{k_x,k_y}^1\tau S_I.$$
(114)
In close proximity of the AFM frequency $`\omega _0`$, we have according to Eq. (13)

$$\omega \omega _0\frac{1}{6}\omega _0^{\prime \prime \prime }\left(kk_0\right)^3,$$
(115)
where $`\omega _0^{\prime \prime \prime }`$ is defined in Eq. (77). Differentiating Eq. (115) with respect to $`k`$

$$\left(\frac{\omega }{k}\right)_{k_x,k_y}\frac{1}{2}\omega _0^{\prime \prime \prime }\left(kk_0\right)^2\frac{6^{2/3}}{2}\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left(\omega \omega _0\right)^{2/3},$$
(116)
and plugging Eq. (116) into (114) yields

$$W_{ex}\frac{2}{6^{2/3}}\tau S_I\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left(\omega \omega _0\right)^{2/3},$$
(117)
where the transmittance $`\tau `$ depends on the incident wave polarization, the frequency $`\omega `$, and the direction of incidence $`(n_x,n_y)=(ck_x/\omega ,ck_y/\omega )`$. Formula (117) implies that the energy density $`W_{ex}`$ and, therefore, the amplitude $`\left|\mathrm{\Psi }_{ex}\left(z\right)\right|=\left|\mathrm{\Phi }_{ex}\right|`$ of the extended mode inside the stack diverge in the vicinity of the AFM regime

$$\left|\mathrm{\Phi }_{ex}\right|\sqrt{W_{ex}}\sqrt{\tau S_I}\left(\omega _0^{\prime \prime \prime }\right)^{1/6}\left|\omega \omega _0\right|^{1/3}\text{as}\omega \omega _0.$$
(118)
The divergence of the extended mode amplitude $`\left|\mathrm{\Phi }_{ex}\right|`$ imposes the similar behavior on the amplitude $`\left|\mathrm{\Psi }_{ev}\left(0\right)\right|=\left|\mathrm{\Phi }_{ev}\right|`$ of the evanescent mode at the slab boundary. Indeed, the boundary condition (111) requires that the resulting field $`\mathrm{\Phi }_T=\mathrm{\Phi }_{ex}+\mathrm{\Phi }_{ev}`$ remains limited to match the sum $`\mathrm{\Phi }_I+\mathrm{\Phi }_R`$ of the incident and reflected waves. The relation (111) together with (118) imply that there is a destructive interference of the extended $`\mathrm{\Phi }_{ex}`$ and evanescent $`\mathrm{\Phi }_{ev}`$ modes at the stack boundary

$$\mathrm{\Phi }_{ex}\mathrm{\Phi }_{ev}K\sqrt{\tau S_I}\left(\omega _0^{\prime \prime \prime }\right)^{1/6}\left(\omega \omega _0\right)^{1/3}\mathrm{\Phi }_{k_0}\text{as}\omega \omega _0,$$
(119)
Here $`\mathrm{\Phi }_{k_0}`$ is the normalized eigenvector of $`T_L`$ in Eq. (80); $`K`$ is a dimensionless parameter. The expression (119) is in compliance with the earlier made statement (79) that the column-vectors $`\mathrm{\Phi }_{ex}`$ and $`\mathrm{\Phi }_{ev}`$ become collinear as $`\omega \omega _0`$.

#### IV.2.1 Space distribution of electromagnetic field in the AFM regime

The amplitude $`\left|\mathrm{\Psi }_{ex}\left(z\right)\right|`$ of the extended Bloch eigenmode remains constant and equal to $`\left|\mathrm{\Phi }_{ex}\right|`$ from (118), while the amplitude of the evanescent contribution to the resulting field decays as

$$\left|\mathrm{\Psi }_{ev}\left(z\right)\right|=\left|\mathrm{\Phi }_{ev}\right|e^{z\mathrm{}k_{ev}},\text{ where }\mathrm{}k_{ev}\frac{\sqrt{3}}{2}6^{1/3}\left(\omega _0^{\prime \prime \prime }\right)^{1/3}\left|\omega \omega _0\right|^{1/3}.$$
(120)
At $`z\left(\mathrm{}k_{ev}\right)^1`$, the destructive interference (119) of the extended and evanescent modes becomes ineffective, and the only remaining contribution to $`\mathrm{\Psi }_T\left(z\right)`$ is the extended mode $`\mathrm{\Psi }_{ex}\left(z\right)`$ with huge and independent of $`z`$ amplitude (118). This situation is graphically demonstrated in Fig. 5.

Let us see what happens when the frequency $`\omega `$ tends to its critical value $`\omega _0`$. According to Eqs. (61) and (109), at $`z=NL,`$ $`N=0,1,2,\mathrm{}`$, the resulting field $`\mathrm{\Psi }_T\left(z\right)`$ inside the slab can be represented as

$$\mathrm{\Psi }_T\left(z\right)=\mathrm{\Phi }_{ex}e^{izk_{ex}}+\mathrm{\Phi }_{ev}e^{izk_{ev}},$$
(121)
Substituting $`k_{ex}`$ and $`k_{ev}`$ from Eq. (78) in (121), and taking into account the asymptotic relation (119), we have

$$\mathrm{\Psi }_T\left(z\right)\left(\mathrm{\Phi }_T+zK\sqrt{\frac{\tau S_I}{\omega _0^{\prime \prime \prime }}}6^{1/3}\left(\frac{i}{2}+\frac{\sqrt{3}}{2}\frac{\omega \omega _0}{\left|\omega \omega _0\right|}\right)\mathrm{\Phi }_{k_0}\right)e^{izk_0}\text{as}\omega \omega _0.$$
(122)
Although this asymptotic formula is valid only for $`z=NL,`$ $`N=0,1,2,\mathrm{}`$, it is obviously consistent with the expression (84) for the non-Bloch solution $`\mathrm{\Psi }_{10}\left(z\right)`$ of the Maxwell equation (38) at $`\omega =\omega _0`$.

#### IV.2.2 The role of the incident wave polarization

The incident wave polarization affects the relative contributions of the extended and evanescent components to the resulting field $`\mathrm{\Psi }_T\left(z\right)`$ in Eq. (109). In addition, it also affects the overall transmittance (2). The situation here is similar to that of the normal incidence, considered in PRB03 . There are two special cases, merging into a single one as $`\omega \omega _0`$. The first one occurs when the elliptic polarization of the incident wave is chosen so that it produces a single extended eigenmode $`\mathrm{\Psi }_{ex}\left(z\right)`$ inside the slab (no evanescent contribution to $`\mathrm{\Psi }_T\left(z\right)`$). In this case, $`\mathrm{\Psi }_T\left(z\right)`$ reduces to $`\mathrm{\Psi }_{ex}\left(z\right)`$, and its amplitude $`\left|\mathrm{\Psi }_T\left(z\right)\right|`$ remains limited and independent of $`z`$. As $`\omega `$ approaches $`\omega _0`$, the respective transmittance $`\tau `$ vanishes in this case, and there is no AFM regime. The second special case is when the elliptic polarization of the incident wave is chosen so that it produces a single evanescent eigenmode $`\mathrm{\Psi }_{ev}\left(z\right)`$ inside the slab (no extended contribution to $`\mathrm{\Psi }_T\left(z\right)`$). In such a case, $`\mathrm{\Psi }_T\left(z\right)`$ reduces to $`\mathrm{\Psi }_{ev}\left(z\right)`$, and the amplitude $`\left|\mathrm{\Psi }_T\left(z\right)\right|`$ decays exponentially with $`z`$ in accordance with Eq. (120). The respective transmittance $`\tau `$ in this latter case is zero regardless of the frequency $`\omega `$, because evanescent modes do not transfer energy. Importantly, as $`\omega `$ approaches $`\omega _0`$, the polarizations of the incident wave that produce either a sole extended or a sole evanescent mode become indistinguishable, in accordance with Eq. (79). In the vicinity of the AFM regime, the maximal transmittance $`\tau `$ is achieved for the incident wave polarization orthogonal to that exciting a single extended or evanescent eigenmode inside the semi-infinite stack.

### IV.3 Tangential energy flux

So far we have been focusing on the axial electromagnetic field distribution, as well as the axial energy flux $`S_z`$ inside semi-infinite slab. At the same time, in and near the AFM regime, the overwhelmingly stronger energy flux occurs in the tangential direction. Let us take a closer look at this problem.

The axial energy flux $`S_z`$ is exclusively provided by the extended contribution $`\mathrm{\Psi }_{ex}\left(z\right)`$ to the resulting field $`\mathrm{\Psi }_T\left(z\right)`$, because the evanescent mode $`\mathrm{\Psi }_{ev}\left(z\right)`$ does not contribute to $`S_z`$. Neither $`\left|\mathrm{\Psi }_{ex}\right|`$ nor $`S_z`$ depends on $`z`$ (see Eq. (130)). By contrast, both the extended and the evanescent modes determine the tangential energy flux $`\stackrel{}{S}_\tau \left(z\right)`$. Besides, according to Eq. (130), the tangential energy flux depends on $`z`$. Far from the AFM regime, the role of the evanescent mode is insignificant, because $`\mathrm{\Psi }_{ev}\left(z\right)`$ is appreciable only in a narrow region close to the slab boundary. But the situation appears quite different near the AFM frequency. Indeed, according to Eq. (120), the imaginary part of the respective Bloch wave vector $`k_{ev}`$ becomes infinitesimally small near the critical point. As a consequence, the evanescent mode extends deep inside the slab, so does its role in formation of $`\stackrel{}{S}_\tau \left(z\right)`$. The tangential energy flux $`\stackrel{}{S}_\tau \left(z\right)`$ as function of $`z`$ can be directly obtained using formula (132) and the explicit expression for $`\mathrm{\Psi }_T\left(z\right)=\mathrm{\Psi }_{ex}\left(z\right)+\mathrm{\Psi }_{ev}\left(z\right)`$. Although the explicit expression for $`\stackrel{}{S}_\tau \left(z\right)`$ is rather complicated and cumbersome, it has very simple and transparent structure. Indeed, the tangential energy flux can be represented in the following form

$$\stackrel{}{S}_\tau \left(z\right)=\stackrel{}{u}_\tau W\left(z\right),$$

where the tangential group velocity $`\stackrel{}{u}_\tau `$ behaves regularly at $`\omega =\omega _0`$. Therefore, the magnitude and the space distribution of the tangential energy flux $`\stackrel{}{S}_\tau \left(z\right)`$ in and near the AFM regime literally coincides with that of the electromagnetic energy density $`W\left(z\right)`$, which is proportional to $`\left|\mathrm{\Psi }_T\left(z\right)\right|^2`$. A typical picture of that is shown in Fig. 5(a).

## V Summary

As we have seen in the previous Section, a distinctive characteristic of the AFM regime is that the incident monochromatic radiation turns into a very unusual grazing wave inside the slab, as shown schematically in Figs. 3 and 6. Such a grazing wave is significantly different from that occurring in the vicinity of the total internal reflection regime, where the transmitted wave also propagates along the interface. The most obvious difference is that near the regime of total internal reflection, the reflectivity approaches unity, which implies that the intensity of the transmitted (refracted) wave vanishes. By contrast, in the case of AFM the light reflection from the interface can be small, as shown in an example in Fig. 4(a). Thus, in the AFM case, a significant portion of the incident light gets converted into the grazing wave (the AFM) with huge amplitude, compared to that of the incident wave. For this reason, the AFM regime can be of great utility in many applications.

Another distinctive feature of the AFM regime relates to the field distribution inside the periodic medium. The electromagnetic field of the AFM can be approximated by a divergent Floquet eigenmode $`\mathrm{\Psi }_{10}\left(z\right)`$ from (84), whose magnitude $`\left|\mathrm{\Psi }_{10}\left(z\right)\right|^2`$ increases as $`z^2`$, until nonlinear effects or other limiting factors come into play. In fact, the field amplitude inside the slab can exceed the amplitude of the incident plane wave by several orders of magnitude, depending on the quality of the periodic array, the actual number of the layers, and the width of the incident light beam.

Looking at the $`z`$ component of light group velocity and energy flux, we see a dramatic slowdown of light in the vicinity of the AFM regime, with all possible practical applications extensively discussed in the literature (see, for example, Joann1 ; Scalora1 ; Agua01 ; Scalora2 , and references therein). In principle, there can be a situation when the tangential components $`(u_x,u_y)`$ of the group velocity also vanish in the AFM regime, along with the axial component $`u_z`$. Although we did not try to achieve such a situation in our numerical experiments, it is not prohibited and might occur if the physical parameters of the periodic array are chosen properly. In such a case, the AFM regime reduces to its particular case – the frozen mode regime with $`\stackrel{}{u}=0`$ inside the periodic medium. This regime would be similar to that considered in PRB03 , with one important difference: it is not related to the magnetic unidirectionality and, hence, there is no need to incorporate nonreciprocal magnetic layers in the periodic array. The latter circumstance allows to realize the frozen mode regime at the infrared, optical, and even UV frequency range.

Acknowledgment and Disclaimer. The efforts of A. Figotin and I. Vitebskiy are sponsored by the Air Force Office of Scientific Research, Air Force Materials Command, USAF, under grant number F49620-01-1-0567. The US Government is authorized to reproduce and distribute reprints for governmental purposes notwithstanding any copyright notation thereon. The views and conclusions contained herein are those of the authors and should not be interpreted as necessarily representing the official policies or endorsements, either expressed or implied, of the Air Force Office of Scientific Research or the US Government.

## VI Appendix 1. $`J`$-unitarity of the transfer matrix

Let $`n\times n`$ matrix $`T\left(z\right)`$ satisfies the following Cauchy problem

$$_zT\left(z\right)=iJA\left(z\right)T\left(z\right),T\left(0\right)=I.$$
(123)
where $`A\left(z\right)=JM\left(z\right)`$ is a Hermitian matrix. Let us prove that the unique solution $`T\left(z\right)`$ for Eq. (123) is a $`J`$-unitary operator

$$T^{}\left(z\right)=JT^1\left(z\right)J^1.$$
(124)
To prove it, notice that Eq. (123) implies

$$_xT^{}\left(z\right)=T^{}\left(x\right)iA\left(z\right)J,T^{}\left(0\right)=I.$$
(125)
Now, let us find the respective Cauchy problem for $`T^1\left(z\right)`$. Since

$$_z\left[T\left(z\right)T^1\left(z\right)\right]=0=T\left(z\right)\left[_zT^1\left(z\right)\right]+\left[_zT\left(z\right)\right]T^1\left(z\right),$$

we have

$$_zT^1\left(z\right)=T^1\left(z\right)\left[_zT\left(z\right)\right]T^1\left(z\right),$$

which in a combination with Eq. (123) yields

$$_zT^1\left(z\right)=T^1\left(z\right)iJA\left(z\right),T^1\left(0\right)=I.$$
(126)
Finally, multiplying both sides of the equality (126) by $`J`$ and using the fact that $`J^2=I`$, we get the following Cauchy problem for $`JT^1\left(z\right)J`$

$$_z\left[JT^1\left(z\right)J\right]=\left[JT^1\left(z\right)J\right]iA\left(z\right)J,JT^1\left(0\right)J=I.$$
(127)
which is identical to that for $`T^{}\left(z\right)`$ from Eq. (125). Since both Cauchy problems (125) and (127) have unique solutions, their similarity implies the relation (124) of $`J`$-unitarity.

## VII Appendix 2. Energy density flux

The real-valued Poynting vector is defined by

$$𝐒\left(\stackrel{}{r}\right)=\frac{c}{8\pi }\mathrm{}\left[𝐄^{}\left(\stackrel{}{r}\right)\times 𝐇\left(\stackrel{}{r}\right)\right].$$
(128)
Substituting the representation (37) for $`𝐄\left(\stackrel{}{r}\right)`$ and $`𝐇\left(\stackrel{}{r}\right)`$ in Eq. (128) yields

$$𝐒\left(\stackrel{}{r}\right)=𝐒\left(z\right)=\frac{c}{8\pi }\mathrm{}\left[\stackrel{}{E}^{}\left(z\right)\times \stackrel{}{H}\left(z\right)\right]$$
(129)
implying that none of the three Cartesian components of the energy density flux $`𝐒`$ depends on the transverse coordinates $`x`$ and $`y`$. Energy conservation argument implies that the component $`S_z`$ of the energy flux does not depend on the coordinate $`z`$ either, while the transverse components $`S_x`$ and $`S_y`$ may depend on $`z`$. Indeed, in the case of steady-state oscillations in a lossless medium we have, with consideration for Eq. (129)

$$𝐒=_zS_z\left(z\right)=0$$

which together with Eq. (129) gives

$$S_z\left(\stackrel{}{r}\right)=S_z=const,S_x\left(\stackrel{}{r}\right)=S_x\left(z\right),S_y\left(\stackrel{}{r}\right)=S_y\left(z\right).$$
(130)
The explicit expression for the $`z`$ component of the energy flux (129) is

$$S_z=\frac{1}{2}\left[E_x^{}H_yE_y^{}H_x+E_xH_y^{}E_yH_x^{}\right]=\frac{1}{2}(\mathrm{\Psi },J\mathrm{\Psi }).$$
(131)
The tangential components of the energy flux can also be expressed in terms of the column vector $`\mathrm{\Psi }\left(z\right)`$ from Eq. (38). Using the expressions (44) for $`E_z`$ and $`H_z`$ and eliminating these field components from $`𝐒\left(z\right)`$ in Eq. (129) yields

$$S_x=\frac{1}{2}(\mathrm{\Psi },\widehat{G}_x\mathrm{\Psi }),S_y=\frac{1}{2}(\mathrm{\Psi },\widehat{G}_y\mathrm{\Psi }),$$
(132)
where $`G_x`$ and $`G_y`$ are Hermitian matrices

$$G_x=\left[\begin{array}{cccc}0& \frac{n_y}{\mu _{33}}& 0& \frac{\epsilon _{13}}{\epsilon _{33}}\\ \frac{n_y}{\mu _{33}}& 2\frac{n_x}{\mu _{33}}& \frac{\mu _{13}^{}}{\mu _{33}}& \frac{\epsilon _{23}}{\epsilon _{33}}\frac{\mu _{23}^{}}{\mu _{33}}\\ 0& \frac{\mu _{13}}{\mu _{33}}& 0& \frac{n_y}{\epsilon _{33}}\\ \frac{\epsilon _{13}^{}}{\epsilon _{33}}& \frac{\epsilon _{23}^{}}{\epsilon _{33}}\frac{\mu _{23}}{\mu _{33}}& \frac{n_y}{\epsilon _{33}}& 2\frac{n_x}{\epsilon _{33}}\end{array}\right],G_y=\left[\begin{array}{cccc}2\frac{n_y}{\mu _{33}}& \frac{n_x}{\mu _{33}}& \frac{\epsilon _{13}}{\epsilon _{33}}+\frac{\mu _{13}^{}}{\mu _{33}}& \frac{\mu _{23}^{}}{\mu _{33}}\\ \frac{n_x}{\mu _{33}}& 0& \frac{\epsilon _{23}}{\epsilon _{33}}& 0\\ \frac{\epsilon _{13}^{}}{\epsilon _{33}}+\frac{\mu _{13}}{\mu _{33}}& \frac{\epsilon _{23}^{}}{\epsilon _{33}}& 2\frac{n_y}{\epsilon _{33}}& \frac{n_x}{\epsilon _{33}}\\ \frac{\mu _{23}}{\mu _{33}}& 0& \frac{n_x}{\epsilon _{33}}& 0\end{array}\right].$$

Both $`G_x`$ and $`G_y`$ are functions of the Cartesian coordinate $`z`$, frequency $`\omega `$, and the direction $`\stackrel{}{n}`$ of incident wave propagation.