Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Dielectric approximation media to reproduce dispersion for field transformation

Open Access Open Access

Abstract

Field transformation is complementary to transformation optics as it converts polarization without reflection for any incidence, despite not being impedance-matched to vacuum. However, to ease manufacture, a dielectric approximation may degrade the transmission performance due to the approximate dispersion relation for the exact reciprocal medium. In this paper, we instead propose a different dielectric equivalent medium to guarantee the exact dispersion. We also discuss the birefringence phenomenon in such dielectric media, and their unavoidable induced reflections.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. INTRODUCTION

It is well known that a planar wave propagates in a homogenous medium, but it is generally analytically elusive what kind of waves can move in inhomogeneous media. Whatsoever, reflection and refraction should generally occur at the boundary of two distinct media. An exception comes around from the mathematical method of transformation optics (TO) [13], which reveals a general recipe to design an inhomogeneous medium aiming for predefined wave behaviors of null reflection, privileged by the state-of-the-art to manufacture artificial electromagnetic materials. The surprising point of reflectionlessness for an incidence waves from air to traversing the transformation media is attributed to the invariance of wave impedance. However, experimental constraints yet require a further dielectric approximation to remove the demand of magnetic materials [4]. Moreover, other reflectionless media are also available by making use of the spatial Kronig–Kramer method and phase-amplitude equations when the impedance match condition is relaxed [58]. The merit of this line of reflectionless media is not only to give fundamental gain to knowledge of wave equations, but also to pave the path towards engineering design to control waves in practical efficiency.

Nevertheless, the field transformation (FT) method also gives reflectionless media a reciprocal material design, which was proposed as complementary to TO [9]. Starting from Tretyakov’s proposal [10], it exploits the local unitary gauge symmetry for electric and magnetic fields simultaneously and provides an additional tunability for their polarization status. FT is intimately connected to the concept of near-to-far FT (Stratton–Chu formula) [11], which has yet not been exploited throughly. It formally packs the electromagnetic (EM) fields into a spinor form as a Riemann–Silberstein vector [1214], and recasts Maxwell equations into a Dirac-like equation. Specifically, it can be adopted in two-dimensional (2D) fabric to keep energy density and the in-plane Poynting vector invariant [9,15,16]. In the 2D form, all fields are reduced to two-component spinor ${({E_z},{\textit{iH}_z})^{T}}$, from which their transverse components can be determined accordingly. This is peculiar as the ideal FT medium remains reflectionless despite not being exactly impedance-matched to vacuum. However, a dielectric approximation is still called on to ease the realization of reciprocal materials, considering the reproduction dispersion relation approximately, which will compromise performance to some extent in the case of polarization converters [17,18].

The dielectric approximation medium [17] cannot fully reproduce perfect transmission and zero reflection as the ideal (exact) FT one does. One could speculate that such imperfection comes from the unfaithful dispersion from the dielectric approximation, which is not the case as our work will show. Hereby, our paper provides a different dielectric medium by reproducing the exact dispersion, instead of dropping a small term in previous work [17]. In Section 2, our dielectric medium will be shown as a biaxial anisotropic medium, which is manufacturable with current metamaterial designs. We also discuss in Section 3 the birefringent behaviors for the light transversing the exact medium and both of its dielectric media and explain how such birefringent rays converge to single transmission after transmitting back to vacuum. The intrinsic reflective property of dielectric approximation media is also discussed to reveal the peculiar non-reflectionlessness for exact FT media designed from the FT approach. Our simulation results indicate that the effort of making the exact dispersion still gives comparable polarization conversion. It seems, in principle, that the underperformance for dielectric media is inevitable because only the exact (magnetic) FT medium is perfect in the discussed polarization transferring scenario. We note that other dielectric mediums can also be made up from tilt layered structures [1921], which is beyond the scope of this paper.

2. EXACT DISPERSION

A. Preliminary

In this section, we start from previous derivation of FT [9,16] to set up preliminaries of our dielectric approximation method. Suppose that our system remains invariant along the $z$ direction, and we separate all EM vectors into transverse ($ T $) and longitudinal ($ z $) components so that a general anisotropic constitutive relation is written as

$$\left({\begin{array}{c}{D_T}\\{D_z}\\{{\textit{iB}_T}}\\{{\textit{iB}_z}}\end{array}} \right) = \left({\begin{array}{cccc}{{\varepsilon _{{TT}}}}&{{\varepsilon _{{Tz}}}}&0&0\\{{\varepsilon _{{zT}}}}&{{\varepsilon _{{zz}}}}&0&0\\0&0&{{\mu _{{TT}}}}&{{\mu _{{Tz}}}}\\0&0&{{\mu _{{zT}}}}&{{\mu _{{zz}}}}\end{array}} \right)\left({\begin{array}{c}{E_T}\\{E_z}\\{{\textit{iH}_T}}\\{{\textit{iH}_z}}\end{array}} \right).$$

Maxwell equations in Heaviside–Lorentz units (H. L. U.) with time-harmonic terms ${e^{- i\omega t}}$ of free charge and current,

$${-}i\omega {\bf D} = c\nabla \times {\bf H},$$
$$i\omega {\bf B} = c\nabla \times {\bf E},$$
can be recast as the propagation of only two coupled $ z $ polarizations ${E_z}$ and ${H_z}$, which also determine the $ T $ components as below:
$$\begin{split}&\frac{1}{{{k_0}}}(\nabla \cdot \hat z \times)\left({\begin{array}{c}{{\textit{iH}_T}}\\{E_T}\end{array}} \right) + \left({\begin{array}{cc}{{\varepsilon _{{zT}}}}&0\\0&{{\mu _{{zT}}}}\end{array}} \right)\left({\begin{array}{c}{E_T}\\{{\textit{iH}_T}}\end{array}} \right) \\&\qquad + \left({\begin{array}{cc}{{\varepsilon _{{zz}}}}&0\\0&{{\mu _{{zz}}}}\end{array}} \right)\left({\begin{array}{c}{E_z}\\{{\textit{iH}_z}}\end{array}} \right) = 0,\end{split}$$
$$\begin{split}\left({\begin{array}{cc}{{\varepsilon _{{TT}}}}&0\\0&{{\mu _{{TT}}}}\end{array}} \right)\left({\begin{array}{c}{E_T}\\{{\textit{iH}_T}}\end{array}} \right) & = - \frac{{\hat z \times \nabla}}{{{k_0}}}\left({\begin{array}{c}{{\textit{iH}_z}}\\{E_z}\end{array}} \right)\\&\quad - \left({\begin{array}{cc}{{\varepsilon _{{Tz}}}}&0\\0&{{\mu _{{Tz}}}}\end{array}} \right)\left({\begin{array}{c}{E_z}\\{{\textit{iH}_z}}\end{array}} \right),\end{split}$$
where the $\hat z$ operator is defined as
$$\hat z \times : = \left({\begin{array}{cc}0&{- 1}\\1&0\end{array}} \right)$$
with its columns and rows for $ x $, $ y $, respectively, $c$ represents the light velocity in vacuum, and ${k_0}$ is the vacuum wave-number ${k_0} = \omega /c$. Equations (4) and (5) are combined into a general wave equation for spinor ${(E_z^{(0)},{\textit{iH}}_z^{(0)})^{T}}$:
$$[\nabla \cdot ({\bf C} \cdot \nabla + \alpha) - {\bf A} - \beta \cdot \nabla]\left({\begin{array}{c}{E_z}\\{{\textit{iH}_z}}\end{array}} \right) = {\bf 0},$$
where
$${\bf C}: = - \left({\begin{array}{cc}0&{\hat z \times}\\{\hat z \times}&0\end{array}} \right){\left({\begin{array}{cc}{{\varepsilon _{{TT}}}}&0\\0&{{\mu _{{TT}}}}\end{array}} \right)^{- 1}}\left({\begin{array}{cc}0&{\hat z \times}\\{\hat z \times}&0\end{array}} \right),$$
$${\boldsymbol \alpha}: = - {k_0}{\bf C}\left({\begin{array}{cc}0&{\hat z \times}\\{\hat z \times}&0\end{array}} \right)\left({\begin{array}{cc}{{\varepsilon _{{Tz}}}}&0\\0&{{\mu _{{Tz}}}}\end{array}} \right),$$
$${\boldsymbol \beta}: = {k_0}\left({\begin{array}{cc}{{\varepsilon _{{zT}}}}&0\\0&{{\mu _{{zT}}}}\end{array}} \right)\left({\begin{array}{cc}0&{\hat z \times}\\{\hat z \times}&0\end{array}} \right){\bf C},$$
$${\bf A}: = - k_0^2\left({\begin{array}{cc}{{\varepsilon _{{zz}}}}&0\\0&{{\mu _{{zz}}}}\end{array}} \right) + \beta \cdot {{\bf C}^{- 1}} \cdot {\boldsymbol \alpha}.$$

The components ${E_z}$ and ${H_z}$ or, equivalently, three components of electric vectors are tangled in the general wave of Eq. (7),

$$\sum\limits_{k = x,y,z} {E_{z,k,k}} = \sum\limits_{m = x,y,z} \frac{{{\varepsilon _{{zm}}}}}{{{c^2}}}\partial _t^2{E_m},$$
which one cannot make an exact geometry [22] out of. The dispersion relation is obtained by demanding ${E_z}$ and ${H_z}$ to take non-trivial solutions, so that the determinant of the coefficient matrix for eigenvector ${({E_z}, \pm {\textit{iH}_z})^{T}}$ with eigenvalues $\pm1$, respectively, stands as
$$\det [{\bf K}(- {\bf C}{{\bf K}^{T}} - {\boldsymbol \alpha}) + {\boldsymbol \beta} \cdot {\bf K} + {\bf A}] = 0,$$
where the symbol matrices are defined as
$${\bf K}: = \left({\begin{array}{cccc}{i{k_x}}&{i{k_y}}&0&0\\0&0&{i{k_x}}&{i{k_y}}\end{array}} \right).$$

In this paper, we take a polarization converter in Fig. 1 as an example to illustrate our dielectric approximation. Throughout such a converter, an incidence wave ${(E_z^{(0)},{\textit{iH}}_z^{(0)})^{T}}$ transverses the FT block while experiencing a unitary local phase rotation $\phi (x,y)$ and converts into a transmission wave ${({E_z},{\textit{iH}_z})^{T}}$, which is massaged into a compact matrix form,

$$\left({\begin{array}{c}{E_z}\\{{\textit{iH}_z}}\end{array}} \right) = \left({\begin{array}{cc}{\cos \phi}&{- \sin \phi}\\{\sin \phi}&{\cos \phi}\end{array}} \right)\left({\begin{array}{c}{E_z^{(0)}}\\{{\textit{iH}}_z^{(0)}}\end{array}} \right).$$
 figure: Fig. 1.

Fig. 1. Schematic of field-transformation where the blue block serves as the FT medium to convert incidence wave (black arrow) to transmission wave (blue arrow) and reflection wave (yellow dashed arrow).

Download Full Size | PDF

Note that the local gauge symmetry above for Maxwell’s equations is reminiscent of the duality transformation, which leaves the in-plane Poynting vector and the spinor modulus invariant [23]. In the case of our polarization converter, the transverse electric (TE) mode converts completely into the transverse magnetic (TM) mode, i.e., ${(E_z^{(0)},{\textit{iH}}_z^{(0)})^{T}}={ (1,0)^{T}},({E_z},{\textit{iH}_z}{)^{T}} = (0, - {1)^{T}}$ if the gauge phase is defined as $\phi = - \pi /2$. This induces cross terms between transverse and longitudinal components in permittivity and permeability matrices, explicitly as

$$\bar {\bar \varepsilon} : = \left({\begin{array}{cc}{{\varepsilon _{{TT}}}}&{{\varepsilon _{{Tz}}}}\\{{\varepsilon _{{zT}}}}&{{\varepsilon _{{zz}}}}\end{array}} \right) = \left({\begin{array}{ccc}n&0&{A_y}\\0&n&{- {A_x}}\\{A_y}&{- {A_x}}&n\end{array}} \right),$$
$$\bar {\bar \mu}: = \left({\begin{array}{cc}{{\mu _{{TT}}}}&{{\mu _{{Tz}}}}\\{{\mu _{{zT}}}}&{{\mu _{{zz}}}}\end{array}} \right) = \left({\begin{array}{ccc}n&0&{- {A_y}}\\0&n&{A_x}\\{- {A_y}}&{A_x}&n\end{array}} \right).$$

Symmetric material matrices are assumed for simplicity, although asymmetric ones are also extensible, for example, in gyromagnetic materials. In order to achieve reflectionlessness, the diagonal term is taken to be unity $n = 1$ to match with vacuum and the gauge vector ${\bf A} = ({A_x},{A_y}{)^{T}} = \nabla \phi /{k_0}$ parallel to the interface. This is the trivial case of null pseudo-magnetic force $\nabla \times {\bf A} = 0$. Resulting from FT, the anisotropic material of Eqs. (16) and (17) is remarkable because it guarantees perfect non-reflection for the incidence wave ${(E_z^{(0)},{\textit{iH}}_z^{(0)})^{T}}$, although it is not impedance-matched to vacuum [cf. Figs. 2(a) and 2(b)].

 figure: Fig. 2.

Fig. 2. (a) and (b) Dispersion contours for ideal FT media in blue solid curves, Li’s approximation medium in purple dashed curves, our dielectric media 1 and 2 in red solid curves, and vacuum in a blue dotted curve. (b) Blue arrow for incidence wave, green arrow for reflected wave, red arrows for refracted wave vectors, and purple arrows for refracted rays. Three-dimensional normal surfaces for (c) ideal FT media and (d) our dielectric media 1 along with its three principal axes marked in arrows; chained red circles for 2D dispersion contours as in (a) and (b). Parameters: $\theta ={ 45^ \circ}$, ${A_x} = 0$, ${A_y} = - 0.25$, $n = 1$.

Download Full Size | PDF

However, such magnetic materials are always difficult to manufacture, and one often resorts to a more realizable scheme of dielectric approximation, which results from enforcing a non-magnetic material while keeping the original wave equation solution [4] or its approximate dispersion [9]. Our dielectric approximation is achieved via reproducing the exact dispersion with dielectrics, as will be detailed below. We note that when separately treating TE and TM modes to keep their dispersions we will obtain an estranged dispersion off of the original one.

B. Li’s Dielectric Approximation

It is straightforward to obtain the explicit dispersion curve for exact medium in Eq. (13), and they are quartic closed curves [24,25], shown as blue chained circles in Fig. 2(a). The convectional dielectric approximation (Li’s approximation) derives from enforcing non-magnetic materials,

$$\bar {\bar \varepsilon} = \left({\begin{array}{ccc}{{n_1}}&0&{b{A_y}}\\0&{{n_1}}&{- b{A_x}}\\{b{A_y}}&{- b{A_x}}&{{n_1}}\end{array}} \right),$$
$$\bar {\bar \mu} = {1_3},$$
which will be substituted into Eq. (13) to give the explicit dispersion
$$\begin{split}&\frac{{k_x^4}}{{{n^2}}} + \frac{{k_y^4}}{{{n^2}}} + \frac{2}{{{n^2}}}k_x^2k_y^2 + \left({- 2k_0^2 - \frac{{2A_x^2k_0^2}}{{{n^2}}} + \frac{{2A_y^2k_0^2}}{{{n^2}}}} \right)k_x^2 \\&\quad + \left({- 2k_0^2 + \frac{{2A_x^2k_0^2}}{{{n^2}}} - \frac{{2A_y^2k_0^2}}{{{n^2}}}} \right)k_y^2 \\ &\quad -\frac{{2{A_x}{A_y}k_0^2}}{{{n^2}}}{k_x}{k_y} - 2k_0^4(A_x^2 + A_y^2) + \frac{{k_0^4}}{{{n^2}}}{(A_x^2 + A_y^2)^2} = 0.\end{split}$$

The solutions for Li’s dielectric approximation are then solved by neglecting quadratic terms of wave-numbers contrasted with the quartic of them, which approximates well except for the grazing angles

$$\begin{split}\frac{{{n^2}}}{{{n_1}}} &= \frac{{A_y^2{b^2}{n^2} - 2{n^2}n_1^2}}{{2(A_y^2 - A_x^2)n_1^2 - 2n_1^2{n^2}}} \\&= \frac{{A_x^2{b^2}{n^2} - 2{n^2}n_1^2}}{{2(A_x^2 - A_y^2)n_1^2 - 2n_1^2{n^2}}} = \frac{{{n^2}{b^2}}}{{4n_1^2}}\end{split}$$
to obtain
$${n_1} = \sqrt {A_x^2 + A_y^2 + {n^2}} ,$$
$$b = 2\sqrt {A_x^2 + A_y^2 + {n^2}} ,$$
and the dispersion curve for an exemplary case ${A_x} = 0$, ${A_y} = - 0.25$, $n = 1$ is shown in dashed purple in Fig. 2(a) for comparison. Note that $b$ can also take opposite values of ${-}2\sqrt {{n^2} - A_x^2 - A_y^2}$, but we dismiss these values since they give the same wave behaviors as the counterpart in Eq. (23).
 figure: Fig. 3.

Fig. 3. (a)–(f) Field plots from wave simulation for incidence angle $\theta ={ 45^ \circ}$ when the beam transverses the designed material for a thickness of $h = \lambda$ demarcated by two black lines. (a) and (b) Ideal FT medium; (c) and (d) Li’s approximation medium; (e) and (f) our dielectric medium 1. (a), (c), and (e) for ${E_z}$ and (b), (d), and (f) for ${H_z}$.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. Conversion ratio $\eta (\theta)$ from TE to TM for a plane wave with incidence angle $\theta$, as in Eq. (24), for four media as listed in the lengend. It samples a line-segment $l$ between the two points $(1.5,1)\lambda$ and $(2.5,1)\lambda$.

Download Full Size | PDF

The wave amplitudes for both (a), (b) exact and (b)–(d) in Li’s approximation media are shown in Fig. 3, where TE converts to the TM mode for an incidence angle of 45° after transversing a thickness of $h = \lambda$, respectively. For Li’s approximation medium, the conversion ratio $\eta$ from TE to TM remains rather close to that of the exact medium until the incidence angle gets more than 45° [cf. red dashed curve in Fig. 4], where the conversion ratio is defined to sample only line-segment $l$ between the two points $(1.5,1)\lambda$ and $(2.5,1)\lambda$ [26]:

$$\eta = \frac{{\int_l \text{d}l|{H_z}{|^2}}}{{\int_l \text{d}l(|{E_z}{|^2} + |{H_z}{|^2})}}.$$

Nevertheless, the conversion ratio of Li’s approximation still outperforms that of the exact FT medium in black in Fig. 4 when approaching the grazing angle of 79.2°.

C. Our Dielectric Media 1 and 2 of Perfect Dispersion

Li’s approximation medium requires three constant diagonal terms for simplicity, but we note that a further degree of freedom can be exploited to remove the dispersion deviation in Li’s approximation medium in the hope of achieving a different approximation. Hence, our dielectric media is designed as follows. In order to achieve the dispersion for exact medium, we allow the diagonal terms of permittivity to differ to realize the perfect dispersion rather than approximately:

$$\bar {\bar \varepsilon} = \left({\begin{array}{ccc}{{n_1}}&0&{b{A_y}}\\0&{{n_1}}&{- l{A_x}}\\{b{A_y}}&{- l{A_x}}&{{n_3}}\end{array}} \right),$$
$$\bar {\bar \mu} = {1_3}.$$

Our strategy taken here is to keep all coefficients for wave-numbers ${k_x}$, ${k_y}$, different from Li’s approximation in Subsection 2.B, which neglects the quadratic terms. We equate every coefficient of dispersion curves in Eq. (20) of our dielectric media and the exact FT medium, both from Eq. (13). This gives

$$\begin{split}\frac{1}{{{n^2}{n_1}}}& = \frac{{- 2 - 2A_x^2/{n^2} + 2A_y^2/{n^2}}}{{A_y^2{b^2} - {n_1}{n_2} - {n_1}{n_3}}}\\& = \frac{{- 2 + 2A_x^2/{n^2} - 2A_y^2/{n^2}}}{{A_x^2{l^2} - {n_1}{n_2} - {n_1}{n_3}}} = \frac{{8{A_x}{A_y}}}{{2{n^2}{A_x}{A_y}bl}},\quad {n_1} \ne 0,\end{split}$$
which boils down to two meaningful solutions
$$\left({\begin{array}{l}{{n_1} = {n^2} - A_x^2 - A_y^2}\\{{n_3} = 3A_x^2 + 3A_y^2 + {n^2}}\\{b = - 2\sqrt {{n^2} - A_x^2 - A_y^2}}\\{l = - 2\sqrt {{n^2} - A_x^2 - A_y^2}}\end{array}} \right.,\quad \left({\begin{array}{l}{{n_1} = {n^2} + A_x^2 + A_y^2 - 2\sqrt {{n^2}(A_x^2 + A_y^2)}}\\{{n_3} = {n^2} - 3A_x^2 - 3A_y^2 + 2n\sqrt {A_x^2 + A_y^2}}\\{b = - 2\frac{{|{A_x}|}}{A_y}\sqrt {+ 2\sqrt {{n^2}(A_x^2 + A_y^2)} - {n^2} - A_x^2 - A_y^2}}\\{l = 2\frac{A_y}{{|{A_x}|}}\sqrt {+ 2\sqrt {{n^2}(A_x^2 + A_y^2)} - {n^2} - A_x^2 - A_y^2}}.\end{array}} \right.$$

Note that the opposite coefficients ${-}b$, ${-}l$ in the cross terms can also be effectively the same for wave propagation and omitted due to their unrealistic values. We shall name the two realistic solutions as mediums 1 and 2, respectively, in the rest of our paper. For the exemplary parameters in Fig. 2, the permittivities for media 1 and 2 are explicitly

$$\begin{split}{{\bar {\bar \varepsilon}} _1} &= \left({\begin{array}{ccc}{0.9375}&0&{0.4841}\\0&{0.9375}&0\\{0.4841}&0&{1.1875}\end{array}} \right),\\ {{\bar {\bar \varepsilon}} _2} &= \left({\begin{array}{ccc}{0.4179}&0&{+ i0.3232}\\0&{0.4179}&{- i0.3232}\\{+ i0.3232}&{- i0.3232}&{1.3321}\end{array}} \right),\end{split}$$
of which, for ${{\bar {\bar \varepsilon}} _2}$, only the solution with positive diagonal terms is taken, whereas the imaginary off-diagonal terms indicate the lossy medium instead of the phase-shifting one.

Our dielectric media 1 and 2 shall guarantee that the dispersion (red) is perfectly aligned with the original one (blue), so that in Fig. 2(b) only the former reveals itself in the picture. However, a perfect dispersion does not assure that our dielectric media is reflectionless as the exact medium. A relook into the FT approach indicates that only the exact medium perfectly converts polarization as designed [9], and both Li’s and our dielectric media have to induce intrinsic reflection for a non-perpendicular incidence as Figs. 3(b)–3(f) demonstrate. Dielectric media always cause reflection in TM modes [panels 3(d) and 3(f)], as they encode less so as not to be completely impedance-matched to the air when enforcing a non-magnetic medium, as will be investigated further in Sec. 3.

The conversion ratio of plane waves for our dielectric medium 1 resembles Li’s approximation within the interval [0, 58.6]° and even beats Li’s around grazing angle $\theta \ge {79.2^ \circ}$, although generally is lower, as seen in Fig. 4. This is the central result of this paper. For our dielectric medium 2 (lossy medium), its conversion ratio diminishes severely for incidence angles greater than 45° [blue dotdashed] due to the complex terms in ${{\bar {\bar \varepsilon}} _2}$, which is very inefficient and thus not in our favor.

3. DISCUSSION

A. Not Perfect Conversion of Polarization

The dielectric media presented above is not an exact equivalent medium. After all, it would seem too good to be true that one could achieve the perfect polarization conversion of the ideal FT medium only via a dielectric one. All dielectric media cannot be exactly equivalent to a magnetic medium. A signature of that is the reflections, which are spotted in Figs. 3(c)–3(f), while they are absent in Figs. 3(a) and 3(b). Only the exact medium indeed perfectly guarantees no reflection, as shown in Figs. 3(a) and 3(b), which appears similar to the impedance-matched medium such as transformation media [16,2729]. Note that although the dispersion curves align exactly, the actual normal surfaces differ when viewed from the same angle in Figs. 2(c) and 2(d), and henceforth we confirm that the dielectric medium cannot fully reproduce the exact medium for waves.

B. Indeed Birefringence

It is natural and familiar in the language of crystal optics to consider birefringence in anisotropic media [30]. It would, at first glance, seem irrelevant for birefringence in previous work on FT, as only one polarization is transferred into another. However, the three-dimensional dispersion for non-trivial electric field ${\bf E}$ in the exact reciprocal medium reveals that it is indeed the case. Considering a planar electric wave ${\bf E}$ in a general magnetic anisotropic medium,

 figure: Fig. 5.

Fig. 5. (a) Dispersion contours for the ideal FT media in blue solid curves, Li’s approximation medium in purple dashed curves, our dielectric media 1 and 2 in red solid curves, and vacuum in a blue dotted curve. Field plots from wave simulation for incidence angle $\theta ={ 80^ \circ}$ for our dielectric medium 1 (b) for ${E_z}$ and (c) for ${H_z}$. Other parameters are the same as Figs. 2(b), 3(e), and 3(f).

Download Full Size | PDF

$${\bf k} \times {{\bar {\bar \mu}}^{- 1}} \cdot ({\bf k} \times {\bf E}) + {\omega ^2}\bar{\bar \epsilon} \cdot {\bf E} = 0,$$
and demanding the determinant of the coefficient matrix before the three-vector electric field vanishes gives the three-dimensional wave normal surface for the electric field in the ideal medium, which is plotted in Fig. 2(c). One also sees that the projection red curve in Fig. 2(c) on the $xy$ plane aligns the blue chained circles in Figs. 2(a) and 2(b). For our dielectric medium 1, its normal surface including its projection dispersion curve are presented in Fig. 2(d), aligning with red circles in Fig. 2(b). As to the polarization states, each of the two rays of beams are one of the two such as ${\psi _ +} = ({E_z},{\textit{iH}_z}{)^{T}}$ and ${\psi _ -} = ({E_z}, - {\textit{iH}_z}{)^{T}}$ in Eq. (15). Then, it is clearly demonstrated that both the exact medium and our dielectric medium 1 are birefringent for an incidence wave with any incidence direction.

Moreover, both dielectric media are biaxial anisotropic with tilting optical axes. With off-diagonal terms in their permittivities, their principal axes are positioned tilted. For our dielectric medium 1, one can calculate the Euler angles [31] and figure out the three principal axes $({k_{{x^\prime}}},{k_{{y^\prime}}},{k_{{z^\prime}}})$ with three principal values {1.5625, 0.9375, 0.5625} marked in Fig. 2(d), so that any three vectors in $({k_x},{k_y},{k_z})$ space are rotated according to

$$\left({\begin{array}{c}{{k_{{x^\prime}}}}\\{{k_{{y^\prime}}}}\\{{k_{{z^\prime}}}}\end{array}} \right) = \left({\begin{array}{ccc}{0.6124}&0&{0.7906}\\0&{- 1}&0\\{- 0.7906}&0&{0.6124}\end{array}} \right)\left({\begin{array}{c}{k_x}\\{k_y}\\{k_z}\end{array}} \right).$$

To implement such a bi-axial dielectric medium, one can build a periodic lattice of printed board circuits for radio frequency waves, which lose their rotational symmetry in the transverse plane [3235]. It is also of interest to see that the two intersection points of dispersion curves on the $({k_x},{k_y})$ plane happen to be two of the four dimples, a general hallmark of a biaxial medium, where the wave group velocity is undefined in the doubly degenerate wave surfaces.

One may also wonder why no bifurcating waves are seen in the transmission beams in Figs. 3(b)–3(f). To answer this question, it is helpful to pinpoint the divergent wave routes through the dispersion curve only. For a TE wave with the incidence angle of 45° [blue arrow in Fig. 2(b)], two distinct wave vectors [red arrows] and two rays [purple arrows] appear after refraction upon the bottom surface. Note that purple arrows are parallel, and hence no birefringence is observed for the transmission beams in Figs. 3(b) and 3(f). The birefringence still appears for the grazing incidence angle. For a near-grazing angle of 80°, the double-refracted wavevectors are marked in two red arrows Fig. 5(a). It is peculiar to note that a partially downward propagation, in fact, derives from refraction instead of reflection. This is confirmed in the reflection part of magnetic field (TM mode), which accounts to the severe reflection in the lower-right region in Fig. 5(c). Therefore, both the exact medium and its dielectric medium 1 hide birefringence beneath the facade of reflectionlessness.

 figure: Fig. 6.

Fig. 6. Amplitude ratio $\zeta (\phi)$ and the phase ${\tau _x}(\phi ,\theta)$ (inset) with the full range of $\phi \in [- \pi /2,0]$ and incidence angle $\theta \in [0,\pi /2]$.

Download Full Size | PDF

For the wave fields, the two separate beams with distinct polarizations are visible from the bifurcating reflection beams in Figs. 3(d) and 3(f). They are not due to double reflection, as they might appear to at first glance, because reflection law requires a single reflection wavevector only ${k_r} = - {k_i}$, as the green arrow indicates in Fig. 2(b). Instead, another reflection beam reflects off of the top surface after transversing the dielectric medium for a short distance.

C. Remark on Polarization Status for Transmitted Waves

Our proposed method is also applicable to other FT devices in the sense to convert linear polarization to other polarization circular ones, similar to previous works by our coauthors [17,18]. Moreover, it is even possible to convert linear polarization to other elliptic polarization when one tunes the gauge phase $\phi$ in the range $[- \pi /2,0]$. We add a remark in this subsection to address the generality of polarization conversion for the ideal FT medium [9], considering the ideal FT medium in Eq. (16) and the incidence wave of TE mode ${(E_z^{(0)},{\textit{iH}}_z^{(0)})^{T}}={ (1,0)^{T}}$. From Eq. (15),

$${E_z} = E_z^{(0)}\cos \phi ,$$
$${\textit{iH}_z} = E_z^{(0)}\sin \phi .$$

According to Eq. (5),

$${E_T} = - \frac{1}{{{k_0}}}\varepsilon _{{TT}}^{- 1}\hat z \times \nabla \left({\begin{array}{c}{E_T}\\{{\textit{iH}_T}}\end{array}} \right) - \varepsilon _{{TT}}^{- 1}\left({\begin{array}{c}{{\varepsilon _{{Tz}}}{E_z}}\\{{\mu _{{Tz}}}{\textit{iH}_z}}\end{array}} \right).$$

For Eq. (16), ${A_x} = 0$, ${A_y} = \phi /({k_0}h)$, with ${k_x} = {k_0}\sin \theta$, ${k_y} = {k_0}\cos \theta$ and planar wave approximation ${E_z} \approx {e^{i({k_x}x + {k_y}y - \omega t)}}$, it is straightforward to derive

$${E_T} = \left({\begin{array}{c}{\sqrt {\frac{{k_y^2\mathop {\sin}\nolimits^2 \phi}}{{k_0^2}} + \frac{{{\phi ^2}\mathop {\cos}\nolimits^2 \phi}}{{4{\pi ^2}}}} {e^{i{\tau _x}}}}\\{\sqrt {\frac{{k_x^2\mathop {\sin}\nolimits^2 \phi}}{{k_0^2}}} {e^{i\frac{\pi}{2}}}}\end{array}} \right),$$
where
$${\tau _x} = \text{ArcTan}2\left[{\frac{\phi}{{2\pi}},\frac{k_y}{{{k_0}}}\sin \phi} \right].$$

Then, it is important to look at the amplitude ratio $\zeta$ and the phase ${\tau _x}$,

$$\zeta = \frac{{|{E_T}|}}{{|{E_z}|}} = \sqrt {\mathop {\tan}\nolimits^2 \phi + \frac{{{\phi ^2}}}{{4{\pi ^2}}}} .$$

As shown in Fig. 6 and its inset, one sees that it is not completely tunable to achieve any polarization status $\zeta$, ${\tau _x}$ for transmitted light via tuning gauge phase $\phi$, because incidence angle $\theta$ however also affects the phase ${\tau _x}$ (see the inset). This may also account for the deterioration of the conversion ratio for the grazing incidence angle in Fig. 4. Therefore, we remark that with the full range of $[- \pi /2,0]$ for the gauge phase one cannot achieve arbitrary polarization generally.

4. CONCLUSION

In this article, to exploit the extra freedom of dielectric approximation media for polarization conversion purposes, we design different biaxial dielectrics to reproduce the exact dispersion for the FT medium. We also further analyze the birefringence phenomenon prevailing in our biaxial anisotropic medium and point out the inevitable reflections for all dielectric equivalent media. Hence, other types of reflectionless anisotropic media are worth future investigations.

Funding

Department of Science and Technology, Hubei Provincial People’s Government (2018CFB148); National Natural Science Foundation of China (NSFC11804087); Hubei University (030-090105).

Acknowledgment

L. Y. acknowledges useful discussion with Xu Donghui, Xu Yadong, and Xu Lin, and also the hospitality received during the single-day visit to Xi Lei’s lab.

Disclosures

The authors declare no conflicts of interest.

REFERENCES

1. U. Leonhardt, “Optical conformal mapping,” Science 312, 1777–1780 (2006). [CrossRef]  

2. J. B. Pendry, “Controlling electromagnetic fields,” Science 312, 1780–1782 (2006). [CrossRef]  

3. U. Leonhardt, “Notes on conformal invisibility devices,” New J. Phys. 8, 118 (2006). [CrossRef]  

4. D. Schurig, J. B. Pendry, and D. R. Smith, “Calculation of material properties and ray tracing in transformation media,” Opt. Express 14, 9794–9804 (2006). [CrossRef]  

5. T. G. Philbin, “Making geometrical optics exact,” J. Mod. Opt. 61, 552–557 (2014). [CrossRef]  

6. B. Vial, Y. Liu, S. A. R. Horsley, T. G. Philbin, and Y. Hao, “A class of invisible inhomogeneous media and the control of electromagnetic waves,” Phys. Rev. B 94, 245119 (2016). [CrossRef]  

7. Y. Liu, B. Vial, S. A. R. Horsley, T. G. Philbin, and Y. Hao, “Direct manipulation of wave amplitude and phase through inverse design of isotropic media,” New J. Phys. 19, 073010 (2017). [CrossRef]  

8. C. King, S. Horsley, and T. Philbin, “Designing scattering-free isotropic index profiles using phase-amplitude equations,” Phys. Rev. A 97, 053818 (2018). [CrossRef]  

9. F. Liu, Z. Liang, and J. Li, “Manipulating polarization and impedance signature: a reciprocal field transformation approach,” Phys. Rev. Lett. 111, 033901 (2013). [CrossRef]  

10. S. A. Tretyakov, I. S. Nefedov, and P. Alitalo, “Generalized field-transforming metamaterials,” New J. Phys. 10, 115028 (2008). [CrossRef]  

11. J. A. Kong, Electromagnetic Wave Theory, 2nd ed. (Wiley, 1990).

12. A. Saa, “Local electromagnetic duality and gauge invariance,” Classical Quantum Gravity 28, 127002 (2011). [CrossRef]  

13. R. P. Cameron, “On the ‘second potential’ in electrodynamics,” J. Opt. 16, 015708 (2014). [CrossRef]  

14. S. M. Barnett, “Optical Dirac equation,” New J. Phys. 16, 093008 (2014). [CrossRef]  

15. A. B. Khanikaev, S. H. Mousavi, W. K. Tse, M. Kargarian, A. H. MacDonald, and G. Shvets, “Photonic topological insulators,” Nat. Mater. 12, 233–239 (2013). [CrossRef]  

16. F. Liu and J. Li, “Gauge field optics with anisotropic media,” Phys. Rev. Lett. 114, 103902 (2015). [CrossRef]  

17. J. Zhao, L. Zhang, J. Li, Y. Feng, A. Dyke, S. Haq, and Y. Hao, “A wide-angle multi-octave broadband waveplate based on field transformation approach,” Sci. Rep. 5, 17532 (2015). [CrossRef]  

18. H. Shi and Y. Hao, “Wide-angle optical half-wave plate from the field transformation approach and form-birefringence theory,” Opt. Express 26, 20132–20144 (2018). [CrossRef]  

19. S. Xi, H. Chen, B.-I. Wu, and J.-A. Kong, “One-directional perfect cloak created with homogeneous material,” IEEE Microw. Wireless Compon. Lett. 19, 131–133 (2009). [CrossRef]  

20. Y. Huang and L. Gao, “Effective negative refraction in anisotropic layered composites,” J. Appl. Phys. 105, 013532 (2009). [CrossRef]  

21. K. Dolgaleva and R. W. Boyd, “Local-field effects in nanostructured photonic materials,” Adv. Opt. Photon. 4, 1–77 (2012). [CrossRef]  

22. U. Leonhardt and T. Philbin, Geometry and Light : The Science of Invisibility, Dover Books on Physics (Dover, 2010).

23. J. D. Jackson, Classical Electrodynamics, 3rd ed. (Wiley, 1999).

24. M. Born and E. Wolf, Principles of Optics (Pergamon, 2009).

25. A. Yariv and P. Yeh, Optical Waves in Crystals: Propagation and Control of Laser Radiation (Wiley, 2003).

26. J. Hao, Y. Yuan, L. Ran, T. Jiang, J. Kong, C. Chan, and L. Zhou, “Manipulating electromagnetic wave polarizations by anisotropic metamaterials,” Phys. Rev. Lett. 99, 063908 (2007). [CrossRef]  

27. M. W. McCall, P. Kinsler, and R. D. M. Topf, “The refractive index of reciprocal electromagnetic media,” J. Opt. 18, 044017 (2016). [CrossRef]  

28. A. Favaro and L. Bergamin, “The non-birefringent limit of all linear, skewonless media and its unique light-cone structure,” Ann. Phys. 523, 383–401 (2011). [CrossRef]  

29. G. Gok and A. Grbic, “Tailoring the phase and power flow of electromagnetic fields,” Phys. Rev. Lett. 111, 233904 (2013). [CrossRef]  

30. E. Hecht, Optics, 5th ed. (Pearson, 2017).

31. H. Goldstein, C. Pool, and J. Safko, Classical Mechanics, 3rd ed. (Higher Education, 2005).

32. Y. Chen, R. Y. Zhang, Z. Xiong, Z. H. Hang, J. Li, J. Q. Shen, and C. T. Chan, “Non-Abelian gauge field optics,” Nat. Commun. 10, 3125 (2019). [CrossRef]  

33. W. Ji, J. Luo, and Y. Lai, “Extremely anisotropic epsilon-near-zero media in waveguide metamaterials,” Opt. Express 27, 19463–19473 (2019). [CrossRef]  

34. F. Liu, T. Xu, S. Wang, Z. H. Hang, and J. Li, “Polarization beam splitting with gauge field metamaterials,” Adv. Opt. Mater. 7, 1801582 (2019). [CrossRef]  

35. W. Cai, U. K. Chettiar, A. V. Kildishev, and V. M. Shalaev, “Optical cloaking with metamaterials,” Nat. Photonics 1, 224–227 (2007). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1.
Fig. 1. Schematic of field-transformation where the blue block serves as the FT medium to convert incidence wave (black arrow) to transmission wave (blue arrow) and reflection wave (yellow dashed arrow).
Fig. 2.
Fig. 2. (a) and (b) Dispersion contours for ideal FT media in blue solid curves, Li’s approximation medium in purple dashed curves, our dielectric media 1 and 2 in red solid curves, and vacuum in a blue dotted curve. (b) Blue arrow for incidence wave, green arrow for reflected wave, red arrows for refracted wave vectors, and purple arrows for refracted rays. Three-dimensional normal surfaces for (c) ideal FT media and (d) our dielectric media 1 along with its three principal axes marked in arrows; chained red circles for 2D dispersion contours as in (a) and (b). Parameters: $\theta ={ 45^ \circ}$ , ${A_x} = 0$ , ${A_y} = - 0.25$ , $n = 1$ .
Fig. 3.
Fig. 3. (a)–(f) Field plots from wave simulation for incidence angle $\theta ={ 45^ \circ}$ when the beam transverses the designed material for a thickness of $h = \lambda$ demarcated by two black lines. (a) and (b) Ideal FT medium; (c) and (d) Li’s approximation medium; (e) and (f) our dielectric medium 1. (a), (c), and (e) for ${E_z}$ and (b), (d), and (f) for ${H_z}$ .
Fig. 4.
Fig. 4. Conversion ratio $\eta (\theta)$ from TE to TM for a plane wave with incidence angle $\theta$ , as in Eq. (24), for four media as listed in the lengend. It samples a line-segment $l$ between the two points $(1.5,1)\lambda$ and $(2.5,1)\lambda$ .
Fig. 5.
Fig. 5. (a) Dispersion contours for the ideal FT media in blue solid curves, Li’s approximation medium in purple dashed curves, our dielectric media 1 and 2 in red solid curves, and vacuum in a blue dotted curve. Field plots from wave simulation for incidence angle $\theta ={ 80^ \circ}$ for our dielectric medium 1 (b) for ${E_z}$ and (c) for ${H_z}$ . Other parameters are the same as Figs. 2(b), 3(e), and 3(f).
Fig. 6.
Fig. 6. Amplitude ratio $\zeta (\phi)$ and the phase ${\tau _x}(\phi ,\theta)$ (inset) with the full range of $\phi \in [- \pi /2,0]$ and incidence angle $\theta \in [0,\pi /2]$ .

Equations (37)

Equations on this page are rendered with MathJax. Learn more.

( D T D z iB T iB z ) = ( ε T T ε T z 0 0 ε z T ε z z 0 0 0 0 μ T T μ T z 0 0 μ z T μ z z ) ( E T E z iH T iH z ) .
i ω D = c × H ,
i ω B = c × E ,
1 k 0 ( z ^ × ) ( iH T E T ) + ( ε z T 0 0 μ z T ) ( E T iH T ) + ( ε z z 0 0 μ z z ) ( E z iH z ) = 0 ,
( ε T T 0 0 μ T T ) ( E T iH T ) = z ^ × k 0 ( iH z E z ) ( ε T z 0 0 μ T z ) ( E z iH z ) ,
z ^ × := ( 0 1 1 0 )
[ ( C + α ) A β ] ( E z iH z ) = 0 ,
C := ( 0 z ^ × z ^ × 0 ) ( ε T T 0 0 μ T T ) 1 ( 0 z ^ × z ^ × 0 ) ,
α := k 0 C ( 0 z ^ × z ^ × 0 ) ( ε T z 0 0 μ T z ) ,
β := k 0 ( ε z T 0 0 μ z T ) ( 0 z ^ × z ^ × 0 ) C ,
A := k 0 2 ( ε z z 0 0 μ z z ) + β C 1 α .
k = x , y , z E z , k , k = m = x , y , z ε z m c 2 t 2 E m ,
det [ K ( C K T α ) + β K + A ] = 0 ,
K := ( i k x i k y 0 0 0 0 i k x i k y ) .
( E z iH z ) = ( cos ϕ sin ϕ sin ϕ cos ϕ ) ( E z ( 0 ) iH z ( 0 ) ) .
ε ¯ ¯ := ( ε T T ε T z ε z T ε z z ) = ( n 0 A y 0 n A x A y A x n ) ,
μ ¯ ¯ := ( μ T T μ T z μ z T μ z z ) = ( n 0 A y 0 n A x A y A x n ) .
ε ¯ ¯ = ( n 1 0 b A y 0 n 1 b A x b A y b A x n 1 ) ,
μ ¯ ¯ = 1 3 ,
k x 4 n 2 + k y 4 n 2 + 2 n 2 k x 2 k y 2 + ( 2 k 0 2 2 A x 2 k 0 2 n 2 + 2 A y 2 k 0 2 n 2 ) k x 2 + ( 2 k 0 2 + 2 A x 2 k 0 2 n 2 2 A y 2 k 0 2 n 2 ) k y 2 2 A x A y k 0 2 n 2 k x k y 2 k 0 4 ( A x 2 + A y 2 ) + k 0 4 n 2 ( A x 2 + A y 2 ) 2 = 0.
n 2 n 1 = A y 2 b 2 n 2 2 n 2 n 1 2 2 ( A y 2 A x 2 ) n 1 2 2 n 1 2 n 2 = A x 2 b 2 n 2 2 n 2 n 1 2 2 ( A x 2 A y 2 ) n 1 2 2 n 1 2 n 2 = n 2 b 2 4 n 1 2
n 1 = A x 2 + A y 2 + n 2 ,
b = 2 A x 2 + A y 2 + n 2 ,
η = l d l | H z | 2 l d l ( | E z | 2 + | H z | 2 ) .
ε ¯ ¯ = ( n 1 0 b A y 0 n 1 l A x b A y l A x n 3 ) ,
μ ¯ ¯ = 1 3 .
1 n 2 n 1 = 2 2 A x 2 / n 2 + 2 A y 2 / n 2 A y 2 b 2 n 1 n 2 n 1 n 3 = 2 + 2 A x 2 / n 2 2 A y 2 / n 2 A x 2 l 2 n 1 n 2 n 1 n 3 = 8 A x A y 2 n 2 A x A y b l , n 1 0 ,
( n 1 = n 2 A x 2 A y 2 n 3 = 3 A x 2 + 3 A y 2 + n 2 b = 2 n 2 A x 2 A y 2 l = 2 n 2 A x 2 A y 2 , ( n 1 = n 2 + A x 2 + A y 2 2 n 2 ( A x 2 + A y 2 ) n 3 = n 2 3 A x 2 3 A y 2 + 2 n A x 2 + A y 2 b = 2 | A x | A y + 2 n 2 ( A x 2 + A y 2 ) n 2 A x 2 A y 2 l = 2 A y | A x | + 2 n 2 ( A x 2 + A y 2 ) n 2 A x 2 A y 2 .
ε ¯ ¯ 1 = ( 0.9375 0 0.4841 0 0.9375 0 0.4841 0 1.1875 ) , ε ¯ ¯ 2 = ( 0.4179 0 + i 0.3232 0 0.4179 i 0.3232 + i 0.3232 i 0.3232 1.3321 ) ,
k × μ ¯ ¯ 1 ( k × E ) + ω 2 ϵ ¯ ¯ E = 0 ,
( k x k y k z ) = ( 0.6124 0 0.7906 0 1 0 0.7906 0 0.6124 ) ( k x k y k z ) .
E z = E z ( 0 ) cos ϕ ,
iH z = E z ( 0 ) sin ϕ .
E T = 1 k 0 ε T T 1 z ^ × ( E T iH T ) ε T T 1 ( ε T z E z μ T z iH z ) .
E T = ( k y 2 sin 2 ϕ k 0 2 + ϕ 2 cos 2 ϕ 4 π 2 e i τ x k x 2 sin 2 ϕ k 0 2 e i π 2 ) ,
τ x = ArcTan 2 [ ϕ 2 π , k y k 0 sin ϕ ] .
ζ = | E T | | E z | = tan 2 ϕ + ϕ 2 4 π 2 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.