# Multiple Fractional Dual Solutions in PEMC Backed Chiral Nihility Slab

Authors: Kanwer Aftab; M. Naveed; Aqeel A. Syed; Shakeel Ahmed; Q. A. Naqvi
Journal: Mathematical Problems in Engineering
 (2011)

---

## Abstract

Inside a chiral nihility slab which is backed by perfect electromagnetic conductor, all valid fractional dual solutions are determined using the fractional curl operator. Fractional dual solutions may be regarded as intermediate step between the original and dual to the original solution. Corresponding fields outside the nihility slab for each case are also determined and analyzed. It is noted that only for one out of three cases the corresponding field outside chiral nihility slab is also fractional dual, while for other two cases, corresponding fields are not fractional dual. It is noted that time average power density vanishes inside the PEC-backed chiral nihility slab, while it exists outside the nihility slab.


---

## Body

## 1. Introduction

Chiral nihility is a special kind of chiral medium, for which the real part of permittivity and permeability are simultaneously zero or refractive index become zero at certain frequency known as nihility frequency. The realization of the chiral nihility medium has been discussed in [1–3]. Tretyakov et al. [1], showed that this material can be realized in a similar way as mixture with “ordinary” nihility by addig chiral inclusions all of the same handedness. At the frequency where the real parts of both permittivity and permeability become zeros, the chirality parameter is nonzero, and, the imaginary parts of all the parameters can be rather small compared to the chirality parameter. It has been proved that at the frequency equal to 6.38 GHz, both the permittivity and permeability of the material have zero real parts as discussed by Tretyakov et al. [2]. Also, Cheng et al. [3] proposed the realization of chiral nihility medium inside waveguide at frequency 0.5 GHz. Chiral nihility is extension of nihility concept, introduced by Lakhtakia [4] and studied by researchers [5–7], for ordinary dielectric medium. In chiral nihility, the two nihility, the two eigenwaves are still circularly polarized but one of them is a backward wave. For backward waves, phase velocity is antiparallel to the corresponding Poynting vector. Phenomena of negative refraction occurs when a plane wave enters from vacuum to chiral nihility medium.Due to its interesting and extraordinary characteristics, chiral nihility has attracted many researchers [8–21]. Constitutive relations for isotropic chiral nihility medium are [2] (1.1)D=-iκϵ0μ0H,B=iκϵ0μ0E.Taking time dependence asexp(iωt), Maxwell equations for chiral nihility medium can be written as (1.2)∇×E=k0κE,∇×H=k0κH,
where k0=ω/c is the wavenumber in vacuum. It must be noted that electric and magnetic fields in chiral nihility medium are not independent but are related through wave impedance (1.3)η=η0lim⁡ϵ→0,μ→0μϵ.It may be noted that, in the present paper, impedance of chiral nihility medium has been considered close to impedance of free spaceη0.According to the duality principle for ordinary isotropic medium, if(E,η0H) is one solution to Maxwell equations, then (η0H,-E) is another solution to Maxwell equations. Solution (η0H,-E) is termed as dual solution to the original solution (E,η0H), because Maxwell equations remain unchanged for both solutions. It is obvious, from above given Maxwell equations for chiral nihility medium, that if (E,ηH) is one solution to Maxwell equations, then other solutions for which Maxwell equations remain unchanged are (∓ηH,±E), (-E,-ηH), and (±ηH,±E). Set of solutions which may be discerned as dual solutions to the original solution (E,ηH) is {(ηH,-E),(-ηH,E),(-E,-ηH)}, as discussed by Naqvi [11]. This means that simultaneous rotation of both electric and magnetic fields in clockwise or counterclockwise direction by amount of integer multiple of π/2 is allowed to get the dual solution. Unallowed dual solution means that one field of original solution is being rotated in clockwise direction, while the other is being rotated counterclockwise and vice versa. That is, (ηH,E) and (-ηH,-E) are invalid dual solutions. Now, if an achiral-chiral-nihility interface is excited by a uniform plane wave, there exists more than one dual to original solution in chiral nihility medium while only one dual to original solution in ordinary medium.For an isotropic and homogeneous medium, the solutions which may be regarded as intermediate step between the original and dual to the original solutions may be obtained using the following relations:(1.4)Efd=1(ik0)α(∇×)αE,η0Hfd=1(ik0)α(∇×)αη0H,
where α is fractional parameter and (∇×)α means fractional curl operator. k0=ωμ0ϵ0 is the wavenumber and η0=μ0/ϵ0 is impedance of the medium. The intermediate solutions between the above two sets of solutions may be obtained by varying parameter α between zero and one. A number of contributions in this direction are available in the published literature [22–40].The other important material which is to be discussed here in this paper is the perfect electromagnetic (PEMC) material. PEMC material is the generalization of PEC and PMC materials. Perfect electromagnetic conductor (PEMC) was introduced by Lindell and Sihvola [41]. PEMC has attracted many researchers to contribute in the field of the field of electromagnetics using this material [42–54].In this paper, we have determined all valid fractional dual solutions inside a chiral nihility slab which is backed by PEMC surface, using the fractional curl operator. Fractional dual solutions inside the PEMC-backed chiral nihility slab and corresponding fields for region outside the slab are determined and analyzed. It is noted that only for one out of three cases corresponding field outside chiral nihility slab is also fractional dual, while for other two cases, corresponding fields are not fractional dual. It is shown that time average power density vanishes inside the PEC-backed chiral nihility slab while it exists outside the nihility slab. It is seen that duality theorem provides a way to obtain solution of a problem without going into analytical details. In ordinary medium, only one solution/problem (known as dual) can be treated if original solution is known, whereas chirality nihility provides opportunity to address more than one such situations. Present study incorporates all possible (three possible dual solutions) combinations of original and dual to original solutions.


## 2. Chiral Nihility Slab Backed by PEMC

Consider a slab of chiral nihility metamaterial. The slab is of infinite length and is backed by PEMC. Front face of the chiral nihility slab is located atz=d1, while perfect electromagnetic conductor is located at location z=d2, where d2>d1.A linearly polarized uniform plane wave, with time dependency time harmonicexp(-jωt), is incident on the chiral nihility slab backed by PEMC. For the sake of simplicity, it is assumed that direction of incidence is parallel to the normal at chiral nihility slab. The electric and magnetic fields inside and outside the grounded chiral nihility slab may be written in terms of unknown coefficients as (2.1)E0=x̂exp⁡(ik0zz)+A-NR-exp⁡(-ik0zz)+B-NL-exp⁡(-ik0zz),z<d1,E1=E+MR+exp⁡(ikz+z)+F+ML+exp⁡(ikz-z)+E-MR-exp⁡(-ikz+z)+F-ML-exp⁡(-ikz-z),d1<z<d2,H0=1η0ŷexp⁡(ik0zz)-iη0{A-NR-exp⁡(-ik0zz)-B-NL-exp⁡(-ik0zz)},z<d1,H1=-iη[E+MR+exp⁡(ikz+z)-F+ML+exp⁡(ikz-z)+E-MR-exp⁡(-ikz+z)-F-ML-exp⁡(-ikz-z)],d1<z<d2,
where (2.2)NR±=x̂±iŷ,NL±=x̂∓iŷ,MR±=x̂±iŷ,ML±=x̂∓iŷ.Superscript ± in (2.2) represents the eigenwaves propagating in the ±z direction. The subscript R and L refer to the RCP and LCP eigenwaves satisfying the dispersion relations (2.3)kz±=±ωκ,
at the nihility frequency. In the above equations, k0z=ωμ0ϵ0 and η=μ/ϵ. It may be noted that relation kz+=-kz- holds for all modes propagating inside the slab.Unknown coefficients in field expressions (2.1) may be obtained by using the related boundary conditions. At z=d2, tangential components of field quantity (ME1+H1) must be [40].Tangential components of electric and magnetic fields across the dielectric interface located at z=d1 must be continuous. Application of these boundary conditions yield unknown coefficients
(2.4)E-=12(Mη+ii-Mη)exp⁡⁡(ik0zd1+ikz+d1),E+=exp⁡⁡(ik0zd1-ikz+d1)2,A-=(Mη+ii-Mη)exp⁡⁡(2ik0zd1)2,B-=(i-MηMη+i)exp⁡(2ik0zd1)2.It may be noted that all coefficients are independent of parameterd2.


## 3. Fractional Dual Solutions in Chiral Nihility Slab

Using the concept of fractional curl operator(∇×)α [14], we can find fractional dual solutions to the Maxwell equations in chiral nihility slab. Discussion has been divided into three cases as:Case 1.
Fractional dual fields between(E,ηH) and (ηH,-E) for region inside and outside the slab may be obtained from Maxwell equations as
(3.1)E0fd±=(k̂i×)αE0±,η0H0fd±=(k̂i×)αη0H0±,E1fd±=(k̂i×)αE1±,ηH1fd=(k̂i×)αηH1±.
Above fields must satisfy the Maxwell equations. Subscript fd stands for fractional dual. Sign ± are for propagation in positive and negative z-direction. In order to deal with above equations, eigenvalues and eigenvector of operators (k̂i×)=±ẑ× are required. So, the eigenvalues and eigenvectors of these cross-product operators are given below. Eigenvalues and eigenvectors of operator k̂1×=ẑ×(3.2)A1=12[x̂+iŷ],a1=-i,A2=12[x̂-iŷ],a2=+i,A3=iẑ,a3=0.
Eigenvalues and eigenvectors of operator k̂2×=-ẑ× are
(3.3)A1=12[x̂+iŷ],a1=+i,A2=12[x̂-iŷ],a2=-i,A3=-iẑ,a3=0.
Fractional dual fields are obtained by fractionalizing the eigenvalues of corresponding linear operator; that is,(3.4)E0fd=[(-i)α12A1+(+i)α12A2]exp⁡⁡(ik0zz)+[(-i)αA-NR-+(i)αB-NL-]exp⁡(-ik0zz),E1fd=[(-i)αE+MR+exp⁡⁡(ikz+z)+(i)αF+ML+exp⁡⁡(ikz-z)+(-i)αE-MR-exp⁡(-ikz+z)+(i)αF-ML-exp⁡(-ikz-z)],η0H0fd=[(-i)α+112A1+(+i)α+112A2]exp⁡⁡(ik0zz)+[(-i)α+1A-NR-+(i)α+1B-NL-]exp⁡(-ik0zz),ηH1fd=[(-i)α+1E+MR+exp⁡(ikz+z)+(i)α+1F+ML+exp⁡⁡(ikz-z)+(-i)α+1E-MR-exp⁡⁡(-ikz+z)+(i)α+1F-ML-exp⁡⁡(-ikz-z)].
Expressing kz- as (-kz+), above equations yields the following:
(3.5)E0fd=[cos⁡(απ2)x̂+sin⁡(απ2)ŷ]exp⁡(ik0zz)-1(Mη)2+1[{((Mη)2-1)cos⁡(απ2)+2Mηsin⁡(απ2)}x̂-{((Mη)2-1)sin⁡(απ2)-2Mηcos⁡(απ2)}ŷ]exp⁡(-ik0z(z-2d1)),E1fd=(-ix̂+ŷ)Mη+i[Mηsin⁡(απ2)-cos⁡(απ2)]exp⁡(ik0zd1+ikz+(z-d1))+(-ix̂-ŷ)Mη+i[Mηsin⁡(απ2)-cos⁡(απ2)][Mη+ii-Mη]×exp⁡(ik0zd1-ikz+(z-d1)),η0H0fd=[cos⁡((α+1)π2)x̂+sin⁡((α+1)π2)ŷ]exp⁡(ik0zz)-1(Mη)2+1[{((Mη)2-1)cos⁡(απ2)+2Mηsin⁡((α+1)π2)}x̂-{((Mη)2-1)sin⁡((α+1)π2)-2Mηcos⁡((α+1)π2)}ŷ]×exp⁡(-ik0z(z-2d1)),ηH1fd=(-ix̂+ŷ)Mη+i[Mηsin⁡((α+1)π2)-cos⁡((α+1)π2)]×exp⁡(ik0zd1+ikz+(z-d1))+(-ix̂-ŷ)Mη+i[Mηsin⁡((α+1)π2)-cos⁡((α+1)π2)][Mη+ii-Mη]×exp⁡(ik0zd1-ikz+(z-d1)).
Boundary conditions can be easily verified at two interfaces. Changing values ofα between zero and one, one can find behavior of fields inside intermediate geometries. For α=0 and α=1, one reproduces the original and dual to original solutions in chiral nihility slab. It may be noted that corresponding field outside the slab is also fractional dual between the original and dual to the original field.Case 2.
Inside the chiral nihility slab, solutions which may be regarded as intermediate step between the original solution(E,ηH) and dual to original solution (-ηH,E) are obtained as
(3.6)E1fd=-ŷ+ix̂Mη+1[Mηsin⁡⁡(απ2)+cos⁡⁡(απ2)]exp⁡(ik0zd1+ikz+(z-d1))+ŷ+ix̂i-Mη[Mηsin⁡⁡(απ2)+cos⁡⁡(απ2)]exp⁡(ik0zd1-ikz+(z+d1)),η0H1fd=(ŷ-ix̂)sin⁡⁡((α+1)π2)exp⁡(ik0zd1+ikz+(z-d1))-(ŷ+ix̂)Mη+ii-Mηsin⁡⁡((α+1)π2)exp⁡(ik0zd1-ikz+(z-d1)).
Corresponding fields outside the nihility slab is given below
(3.7)E0fd=[cos⁡⁡(απ2)x̂-sin⁡⁡(απ2)ŷ]exp⁡⁡(ik0zz)-11+(Mη)2[(((Mη)2-1)cos⁡⁡(απ2)-2Mηsin⁡⁡(απ2))x̂+(((Mη)2-1)sin⁡⁡(απ2)+2Mηcos⁡⁡(απ2))ŷ]×exp⁡(-ik0z(z-2d1)),η0H0fd=[cos⁡⁡((α+1)π2)x̂-sin⁡⁡((α+1)π2)ŷ]exp⁡⁡(ik0zz)-11+(Mη)2×[(((Mη)2-1)cos⁡⁡((α+1)π2)-2Mηsin⁡⁡(α+1π2))x̂+(((Mη)2-1)×sin⁡((α+1)π2)+2Mηcos⁡⁡((α+1)π2))ŷ]×exp⁡(-ik0z(z-2d1)).
Equivalent cross-product operator for the corresponding field outside the nihility slab is (k̂i×)=∓ẑ×. It may be noted that outside fields are not fractional dual solutions.Case 3.
Solutions which may be regarded as intermediate step between the original solution(E,ηH) and (-E,-ηH) are obtained as
(3.8)E1fd=-ŷ+ix̂Mη+1[cos⁡⁡(απ)-Mηsin⁡⁡(απ)]exp⁡(ik0zd1+ikz+(z-d1))+ŷ+ix̂i-Mη[cos⁡⁡(απ)-Mηsin⁡⁡(απ)]exp⁡(ik0zd1-ikz+(z-d1)),η0H1fd=(ŷ-ix̂)sin⁡⁡((2α+1)π2)exp⁡(ik0zd1+ikz+(z-d1))-(ŷ+ix̂)Mη+ii-Mηsin⁡⁡((2α+1)π2)exp⁡(ik0zd1-ikz+(z-d1)).
Corresponding fields outside the nihility slab are(3.9)E0fd=[cos⁡(απ)x̂+sin⁡(απ)ŷ]exp⁡⁡(ik0zz)-11+(Mη)2[((Mη)2-1)cos⁡⁡(απ)+2Mηsin⁡⁡(απ)]x̂-[((Mη)2-1)sin⁡⁡(απ)-2Mηcos⁡⁡απ]ŷexp⁡⁡(-ik0z(z-2d1)),η0H0fd=[cos⁡((2α+1)π2)x̂+sin⁡((2α+1)π2)ŷ]exp⁡⁡(ik0zz)-11+(Mη)2[((Mη)2-1)cos⁡((2α+1)π2)+2Mηcos⁡⁡((2α+1)π2)]x̂-[((Mη)2-1)sin⁡((2α+1)π2)-2Mηsin⁡⁡((2α+1)π2)]ŷ×exp⁡(-ik0z(z-2d1)).Equivalent cross-product operator for the corresponding field outside the nihility slab is (k̂i×)=±iẑ×. It may be noted that outside fields do not construct fractional dual solution.
The average power densityPav(z) can be obtained using the relation
(3.10)Pav(z)=12R(E(z)×H*(z)).
It is noted that average power density is zero inside the chiral nihility slab, whereas it has nonzero value outside the slab.


## 4. Conclusions

It is concluded that inside the chiral nihility slab, three dual to original solution are possible, while outside the chiral nihility slab, only one dual to original solution exist. Out of three dual solutions for region inside the chiral nihility slab, only one dual solution has correspondence with the dual solution for region outside the slab. Fractional dual solutions inside the PEMC-backed chiral nihility slab and corresponding fields for region outside the slab are determined and analyzed. Average power density outside the slab is nonzero, whereas it yields zero for region inside the slab.



---

*Source: 402943-2011-10-20.xml*