# Focusing, power tunneling and rejection from chiral and/or chiral nihility/nihility metamaterials layers.

1. INTRODUCTIONChiral media have been well known in the optical frequency range because of its special properties. Of these properties, optical activity and circular dichroism are of particular interest [1-6]. Chiral medium can be thought as composed of numerous randomly oriented chiral objects which can never be brought into congruence with their mirror images by any translation or rotation [7-10]. Many objects found in nature can be regarded as chiral: such as irregular tetrahedrons, sugar molecules and wire helices. In fact the word "chirality" is derived from the Greek word "cheir" having meaning of hand; an object possessing above mentioned property of incongruence with its mirror image. In science particular word "chirality" instead of dissymmetry was first introduced by Lord Kelvin, Professor of natural philosophy in the University of Glasgow. Now this word is used to describe the media microscopically composed of the chiral objects.

When linearly polarized wave falls on a slab of chiral medium, it splits into two circularly polarized waves: one left circularly polarized and the other right circularly polarized [8-10]. After passing through the slab the two waves combine to yield a linearly polarized wave whose plane of polarization is rotated with respect to that of the plane of polarization of the incident wave. Effect of chirality of the medium which rotates the plane of linearly polarized wave passing through it is studied by Arago [1], Biot [2] and Fresnel [3]. This phenomenon is named as optical activity of the medium. The amount of rotation depends upon the distance traveled by the wave in the slab. This phenomenon is named as circular birefringence. The phenomena circular dichroism alludes to the fact that different amount of field absorbtion for the left and right circularly polarized waves occurs as they pass through the chiral media. The constitutive relations for chiral metamaterial [8] are given below,

D = [epsilon]E + i[kappa]H (1)

B = [mu]H - i[kappa]E (2)

The behavior of propagation and radiation of the field in chiral medium has been investigated by several authors [8-10].

Left-handed or double negative (DNG) materials are metamaterials in which electric, magnetic and the wave vectors follow the so called left hand rule. In fact the behavior of electromagnetic waves in "left handed media" goes back a long way, e.g., Lamb [11], Schuster [12], Pocklington [13] and Malyuzhinets [14] contributed on this topic. The actual study of waves in media with simultaneously negative permittivity and permeability was first done by Sivukhin [15]. Veselago in detail studied the interesting property of negative refraction in left handed materials [16]. In chiral medium, negative real parts of the permittivity and permeability lead an isotropic chiral medium to exhibit circular dichroism that is reverse with respect to that exhibited by an identical medium but with positive real parts of permittivity and permeability. Left handed materials exhibit some other interesting properties, namely, backward wave and double negative parameters [17]. Negative refraction can also be achieved by using the chiral metamaterials [1821]. In chiral matematerials strong enough chirality produces negative refractive index for one of the eigenwaves [22,23].

In electromagnetics, Lakhtakia introduced the concept of nihility metamaterial for a material whose permittivity and permeability approach to zero [24]. Constitutive relations for nihility metamaterial are,

D = 0 (3)

B = 0 (4)

Lakhtakia showed that propagation is not possible in nihility metamaterial. Later, Tretyakov et al. [25] extended this concept of nihility for the isotropic chiral medium. A chiral nihility metamaterial is one for which, at certain value of frequency, real parts of permittivity and permeability both are simultaneously zero. That is, [epsilon] [right arrow] 0, [mu] [right arrow] 0 and [kappa] [not equal to] 0 at nihility frequency. Constitutive relations for chiral nihility metamaterials are,

D = i[kappa]H (5)

B = -i[kappa]E (6)

Chiral nihility metamaterial has two wavenumbers of equal magnitude but opposite signs. When a dielectric-chiral nihility half space is excited by an oblique incident linearly polarized plane wave, two circularly polarized plane waves having opposite handedness are produced. One of them is forward wave while the other is backward wave which gives rise to the phenomenon of negative refraction in chiral nihility metamaterial.

If a perfect electric conductor (PEC) interface placed in chiral nihility metamaterial is excited by a forward plane wave. Only one reflected wave parallel to the incident wave is produced, i.e., negative reflection happens. The reflected wave is backward wave which cancels the incident forward wave to produce zero power propagation. Interestingly, if a chiral nihility slab backed by PEC interface is excited by linearly polarized plane wave, effect of each eigenwave in chiral nihility after reflection from PEC interface is canceled leading to a situation as if front face of the geometry is PEC and there is no chiral slab [26]. A PEC waveguide coated with chiral nihility metamaterial confines propagation of power within un-coated region of the guide [27].

In electromagnetics, Brewster angle is defined as angle of incidence for which reflection power is zero. In all introductory text books on electromagnetic theory [28,29], reflection of vertically polarized (parallel polarized or transverse magnetic (TM)) plane wave and horizontally polarized (perpendicular polarized or transverse electric (TE)) plane wave from a dielectric-dielectric half space boundary is usually presented to explain the concept of Brewster angle. It is mentioned that vertically polarized wave experiences zero reflection at one incident angle whereas, reflection is non-zero for horizontal polarization for all angles. This means that if dielectric-dielectric half space is excited by a wave, having both the vertical and horizontal polarizations, at the Brewster angle the reflected wave will be linearly polarized with horizontal polarization only. The angle of incidence that allows total reflection of power from a planar dielectric-dielectric interface is known as critical angle. Critical angle exists only for perpendicular polarization if the wave propagates from a denser dielectric medium to a rare dielectric.

In the same context, reflection of plane wave from an interface of an achiral dielectric and a chiral/chiral nihility metamaterial was discussed by Qiu et al. [30]. Their study reveals that for certain values of constitutive parameters, the results opposite to that of conventional dielectric-dielectric interface can also be achieved. That is, existence of zero reflection for only parallel polarization in the dielectric-chiral case is possible. They also showed that, for certain values it is possible to have no Brewster angle for both polarizations and total reflection can be achieved for a wide range of incident angles. Power corresponding to the electromagnetic waves associated with a planar interface of two chiral and/or chiral nihility metamaterials excited by a plane wave was studied by Ahmad et al. [31]. Their results suggest interesting characteristics such as complete power transmission/rejeection and band pass/band reject filter of the interface.

2. GO FIELD AND MASLOV'S METHOD

Geometrical optics (GO) is a powerful tool for the study of wave motion at high frequencies [32-34]; however, it fails at the caustics. A caustic is a region where area of a ray tube is zero and hence introduces a mathematical singularity when ever intensity per unit area is needed to be counted, though physically this is not the case. Physically the field is always finite at caustics and caustics are of practical interest in many applications including defense and medical sciences. Maslov proposed an alternative method to find the fields in the caustic region [35]. Maslov's method combines the simplicity of asymptotic ray theory and the generality of the Fourier transform method. This is achieved by representing the GO fields in terms of mixed coordinates consisting of wave vector coordinates and spatial coordinates. Maslov's method have been applied to study the field near the caustics of the focusing systems by many authors [36-45]. Focusing from a PEC cylindrical and spherical reflectors using Maslov's method was first treated by Ji and Hongo [38]. Faryad and Naqvi extended this work for reflector coated with chiral metamaterial [44]. Illahi and Naqvi studied the focusing from chiral nihility coated PEC and perfect electromagnetic conductor (PEMC) cylindrical reflector [46]. In the present work, one interest is also to unify all these cases so that these become special cases of the geometry, under investigation.

Consider a three dimensional vector wave equation, in cartesian coordinates r = (x, y, z), describing the field in medium having wavenumber [k.sub.0]

[[nabla].sup.2]U(r) + [k.sup.2.sub.0]U(r) = 0

Expressing the solution U(r) of wave equation in terms of well known Luneberg-Kline series yields Eikonal equation for phase s(r) and transport equation for amplitude. For homogeneous and lossless medium, Eikonal equation reduces to Hamiltonian equation as

H(r, p) = (p x p - 1)/2 = 0, p = [nabla]s (7)

where p is the wave vector. The solution of Hamiltonian equation is given below

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where ([xi], [eta], [zeta]) and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are the initial values of cartesian coordinates (x, y, z) and wave vector coordinates ([P.sub.x], [P.sub.y], [P.sub.z]) respectively and t is parameter along the ray. Finite field around caustic may be obtained using following expression obtained from Maslov's method [34, 40]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where U(r) denotes the electric field strength. [D(t)/D(0)] is the Jacobian for transformation from ray coordinates to cartesian coordinate and so represents the initial phase. Above equation provides uniform solution around the caustic.

In this paper, transmitted and reflected powers from two parallel layers of chiral and/or chiral nihility metamaterials with changing angle of incidence and chiralities of the layers are analyzed. Limiting cases of permittivity, with value of chirality fixed to zero, are also taken into account and their interpretation given. Focusing from a large size paraboloidal reflector, composed of two layers of chiral and/or chiral nihility metamaterials, is studied by utilizing the Maslov's method to give the remedy of GO which fails at caustic.

3. FORMULATION

Consider a geometry consisting of two infinite parallel layers of chiral metamaterials placed in air as shown in Figure 1. Physical width of the two layers is taken as [d.sub.1] and [d.sub.2], respectively. Constitutive parameters for the two layers are denoted by ([[epsilon].sub.1], [[mu].sub.1], [[kappa].sub.1]) and ([[epsilon].sub.2], [[mu].sup.2], [[kappa].sub.2]), respectively. Air in which layers are placed have constitutive parameters [[epsilon].sub.0] and [[mu].sub.0].

The above mentioned geometry is excited by an oblique incident linearly polarized electromagnetic plane wave. [[theta].sub.i], [[theta].sub.r] and [[theta].sub.t] are the angles of incidence, reflection and transmission with respect to the walls of layers, respectively. As chiral medium supports two circularly polarized eigen waves having different wavenumbers so fields in each layer are written as linear combination of left circularly polarized (LCP) and right circularly polarized (RCP) plane waves. Both LCP and RCP waves propagate in forward as well as backward direction in each layer. The wavenumbers of both eigen waves in a chiral media are different. Refractive indices, wavenumbers and intrinsic impedance of the medium for first chiral layer are

[n.sub.1(R,L)] = [square root of ([[epsilon].sub.r1][[mu].sub.r1])] [+ or -] [[kappa].sub.1] (8)

[k.sub.1(R,L)] = [omega] ([square root of ([[epsilon].sub.1][[mu].sub.1])] [+ or -] [[kappa].sub.1]) (9)

[[eta].sub.1] = [square root of ([[[mu].sub.1]/[[epsilon].sub.1])] (10)

where [n.sub.1R] and [k.sub.1R] are refractive index and wavenumber for RCP wave, respectively. Quantity [[epsilon].sub.r1] = [[epsilon].sub.1]/[[epsilon].sub.0] is the relative permittivity of the medium occupying first layer.

For second chiral layer, these quantities are

[n.sub.2(R,L)] = [square root of ([[epsilon].sub.r2][[mu].sub.r2])] [+ or -] [[kappa].sup.2] (11)

[k.sub.2(R,L)] = [omega] ([square root of ([[epsilon].sub.2][[mu].sub.2])] [+ or -] [[kappa].sub.2]) (12)

[[eta].sub.2] = [square root of ([[mu].sub.2]/[[epsilon].sub.2])]. (13)

Wavenumber and intrinsic impedance of the host medium are given as,

[[eta].sub.0] = [square root of ([[mu].sub.0]/[[epsilon].sub.0])] (14)

[k.sub.0] = [omega] [square root of ([[mu].sub.0][[epsilon].sub.0])]. (15)

Expressions for the incident and reflected fields and the field transmitted after passing through both the layers are given in the following:

[E.sub.inc] = [[E.sub.i[parallel]] ([a.sub.z] sin [[theta].sub.i] + [a.sub.x] cos [[theta].sub.i]) + [E.sub.i] [perpendicular to][a.sub.y]] exp[-i[k.sub.0](z cos [[theta].sub.i] - x sin [[theta].sub.i])} (16)

[H.sub.inc] = 1/[[eta].sub.0][- [E.sub.i] [perpendicular to] ([a.sub.z] sin [[theta].sub.i] + [a.sub.x] cos [[theta].sub.i]) + [E.sub.i[parallel]][a.sub.y]] exp[-i[k.sub.0](z cos [[theta].sub.i] - x sin [[theta].sub.i])] (17)

[E.sub.ref] = [[E.sub.r[parallel]] (-[a.sub.z] sin [[theta].sub.r] + [a.sub.x] cos [[theta].sub.r]) + [E.sub.r][perpendicular to][a.sub.y]] exp[i[k.sub.0](z cos [[theta].sub.r] + x sin [[theta].sub.r])] (18)

[H.sub.ref] = 1/[[eta].sub.0][[E.sub.r][perpendicular to] (-[a.sub.z] sin [[theta].sub.r] + [a.sub.x] cos [[theta].sub.r]) - [E.sub.r[parallel]][a.sub.y]] exp[i[k.sub.0](z cos [[theta].sub.r] + x sin [[theta].sub.r])] (19)

[E.sub.tra] = [[E.sub.t[parallel]] ([a.sub.z] sin [[theta].sub.t] + [a.sub.x] cos [[theta].sub.t]) + [E.sub.t][perpendicular to][a.sub.y]] exp[-i[k.sub.0](z cos [[theta].sub.t] - x sin [[theta].sub.t])] (20)

[H.sub.tra] = -1/[[eta].sub.0] [[E.sub.t][perpendicular to] ([a.sub.z] sin [[theta].sub.t] + [a.sub.x] cos [[theta].sub.t]) - [E.sub.t[parallel]] [a.sub.y]] exp[-i[k.sub.0](z cos [[theta].sub.t] - x sin [[theta].sub.t])] (21)

In the above equations, [E.sub.r[parallel]], [E.sub.r[perpendicular to]], [E.sub.t[parallel]], and [E.sub.t[perpendicular to]] are unknown coefficients corresponding to parallel and perpendicular components of the respective fields.

Expressions of fields inside each chiral layer can be written as,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

where [k.sub.pL] and [k.sub.pR] represent the wavenumbers for LCP and RCP waves, respectively. Value of p, either 1 or 2, describes fields for first or the second layer, whereas, subscripts b and f are used to represent forward and backward waves, respectively. In the above equations, [E.sub.fL], [E.sub.fR], [E.sub.bL], and [E.sub.bR] are unknown coefficients. Tangential components of fields are continuous at interfaces located at z = 0, [d.sub.1], [d.sub.1] + [d.sub.2] and these boundary conditions are used to find the unknown coefficients. Snell's law is used to determine the relation among the angles of incidence, reflection and transmission by using following relations:

[[theta].sub.r] = [[theta].sub.i]

[k.sub.0] sin [[theta].sub.i] = [k.sub.1(R,L)] sin [[theta].sub.1(R,L)] = [k.sub.2(R,L)] sin [[theta].sub.2(R,L)] = [k.sub.t] sin [[theta].sub.t],

where [k.sub.0] and [k.sub.t] denote the wavenumbers for medium before and after the layers and are equal when these two layers are placed in air. Using transfer matrix method (TMM) discussed in [47}, three matching matrices each relating the fields at an interface are obtained. Moreover, two 4 x 4 propagation matrices, linking fields between two consecutive interfaces are also obtained. Numerical results for powers corresponding to parallel and perpendicular components of transmitted and reflected fields are obtained in the following. In Section 2, it is assumed that optical width of each layer is [[lambda].sub.0]/4, where [[lambda].sub.0] is wave length in air corresponding to the operating frequency.

4. RESULTS AND DISCUSSION

In this section, numerical results for powers corresponding to parallel and perpendicular components of the reflected and transmitted fields are presented and analyzed. Four different cases, summarized below, are considered for this purpose.

i. Both layers are of chiral metamaterials (c-c).

ii. First layer is of chiral nihility and second is of chiral metamaterial (cn-c).

iii. First layer is of chiral whereas second is of chiral nihility metamaterial (c-cn).

iv. Both layers are of chiral nihility metamaterials (cn-cn).

Results corresponding to above cases are shown by Figure 3 to Figure 9.

In addition to the above, the limiting cases of permittivity of second layer while keeping [[kappa].sub.2] = 0 are also considered for the discussion. In these situations, the above four cases essentially present following scenario.

a. First layer is of chiral metamaterial whereas permittivity of second layer is taken very large with chirality [[kappa].sub.2] equal to zero making it essentially a PEC boundary (c-PEC).

b. First layer is of chiral nihility metamaterial whereas permittivity of second layer is taken very large with chirality [[kappa].sub.2] equal to zero (cn-PEC).

c. First layer is of chiral. metamaterial whereas permittivity of second layer is taken very small with chirality [[kappa].sub.2] equal to zero making it essentially a nihility interface (c-n).

d. First layer is of chiral nihility metamaterial whereas permittivity of second layer is taken very small with chirality [[kappa].sub.2] equal to zero (cn-n).

Results corresponding to above situations are depicted through Figure 10 to Figure 13. [parallel] and [perpendicular to] signs are used to represent power for parallel and perpendicular components, respectively. For example: [c-c.sub.[parallel]] and [c-c.sub.[perpendicular to]] represent parallel and perpendicular components of power for chiral-chiral cases, respectively.

In Figure 2, power corresponding to parallel and perpendicular components of reflected fields versus angle of incidence is shown. It is assumed that the intrinsic impedances of metamaterials filling both layers are same. Figure contains plots for four different cases, i.e., c-c, cn-c, c-cn, and cn-cn. Brewster angle is defined as the angle of incidence for which reflection of power, both for parallel and perpendicular components, is zero. Reflected powers for c-c case are zero for [[theta].sub.i] < 23[degrees] and this gives a range of brewster angles. In the case of cn-cn, a range of Brewster angles exists for [[theta].sub.i] < 17.5[degrees]. After this range, power for parallel component increases whereas power for perpendicular component remains negligible for the entire range of incident angles. In the case of cn-c, parallel component is nonzero everywhere and near grazing incidence almost total reflection of power happens whereas power for perpendicular component remains negligible through the entire range of incident angles. It may be noted that range of Brewster angles exist only for c-c and cn-cn cases and range of angles yielding almost total reflection exists only for parallel component of cn-c case. For c-cn case, power only for perpendicular component is zero for [[theta].sub.i] < 27.5[degrees] and corresponding parallel component is nonzero everywhere.

Figure 3 describes behavior of the parallel and perpendicular components of transmitted powers for impedance matching of both layers. For c-c case, power for parallel and perpendicular components have same initial amplitude, i.e., 0.5. Considering cn-c case, power of both parallel and perpendicular components of transmitted power follow each other and decrease gradually with incidence angle. In c-cn case, perpendicular component is dominant over parallel component whereas for cn-cn case, it is observed that power of parallel component is dominant over the perpendicular component.

Figure 4 shows the behavior of reflected power for layers having different intrinsic impedances. For c-c and cn-c cases, no Brewster angle exists and power for parallel component is dominant over corresponding perpendicular component. For c-c case, reflected power only for perpendicular component is zero for [[theta].sub.i] < 34.3[degrees] and [[theta].sub.i] > 85[degrees]. For cn-c case, when [[theta].sub.i] > 22[degrees], power for parallel component of reflected field starts increasing and after [[theta].sub.i] = 80[degrees], almost total reflection is observed because perpendicular component has insignificant reflected power. Perpendicular component has very low value between [[theta].sub.i] = 34.3[degrees] and [[theta].sub.i] = 78[degrees] and is zero elsewhere. Both c-cn and cn-cn cases have same range of Brewster angles, i.e., [[theta].sub.i] [less than or equal to] 11.46[degrees].

Figure 5 is about transmitted power in case of impedance mismatch. For c-c case, both components of power have equal contribution at [[theta].sub.i] = 31.5[degrees]. For cn-c and cn-cn cases, power for parallel component has higher values and follows the power for corresponding perpendicular component.

Figure 6 shows plots of power for impedance matching case versus chirality of the first layer. For the purpose of analysis chirality values are considered between 0 and 5. For cn-c case, almost total power reflection in terms of parallel component occurs at very low values of chirality and perpendicular component is almost zero for considered range of chirality. For c-cn case, zero reflection of parallel power component is obatined at chirality value [[kappa].sub.1] = 1.9 whereas at [[kappa].sub.1] = 1.8, 2.2, 4 perpendicular component of reflected power is zero. In cn-cn case, oscillating behavior is observed for both components with periodic behavior for cn-[cn.sub.[perpendicular to]]. At [[kappa].sub.1] = 0, 2, 4 perpendicular component of power becomes zero and at [[kappa].sub.1] = 1, 3, 5 parallel component of power is zero. At [[kappa].sub.1] = 0, total reflection occurs. Zeros of parallel and perpendicular components of power at specific values of chirality can be used to get or avoid particular component of polarization.

Figure 7 gives transmission behavior of power. It is noted that both components for each case contain either zeros or values approaching to zero. These findings may be used to get or avoid particular component of polarization. In Figure 8, no Brewster angle exist for all four cases except for cn-c case where both components are negligibly small at [[kappa].sub.2] = 1.9. It is also noted that perpendicular components of power for cn-cn and cn-c are negligibly small. Conclusions similar to Figure 7 may be drawn from Figure 9.

Figure 10 shows the plots for reflected power assuming very high value of permittivity of second layer. This situation may be considered as a chiral layer backed by perfect electric conductor interface. In general reflected power has both components, i.e., both c-PE[C.sub.[parallel]] and c-PE[C.sub.[perpendicular to]] are nonzero. But for certain values of chirality of first layer, reflected field contains only co-polarized component. Contribution of parallel component is very large as compared to the perpendicular component of reflected power. For cn-PEC case, complete reflection in terms of parallel component is observed. Figure 11 shows zero transmission of power for c-PEC and cn-PEC cases.

Another interesting case is shown in Figure 12 where second layer is assumed of nihility metamaterial. This has been achieved by taking [[kappa].sub.2] = 0 and setting very low value of permittivity of second layer. This situation may be seen as a chiral layer backed by nihility interface. For c-n case, almost total reflection in terms of parallel component of the reflected power is observed. For the case of cn-n, perpendicular component of reflected power has also significant contribution for certain range of chirality. It is noted that for certain range of chirality, overlap of cn-n and c-n cases happen. Figure 13 shows plots of transmitted power for c-n and cn-n cases and no transmission of power is obtained. It may be noted that conclusions drawn from Figure 10 to Figure 13 agree with published work.

Focusing of electromagnetic waves from a paraboloidal reflector composed of chiral and/or chiral nihility metamaterial using the Maslov's method is studied in next section. PEC and nihility backed chiral/chiral nihility paraboloidal reflector are also discussed.

5. EVALUATION OF FINITE FIELD AROUND FOCUS

Consider a bilayer paraboloidal reflector placed in air as shown in Figure 14. The equation for the surface of a paraboloidal reflector is given by,

[zeta] = g([xi], [eta]) = f - [[[rho].sup.2]/4f] = f - [[xi].sup.2] + [[eta].sup.2]/4f] (26)

where ([xi], [eta], [zeta]) are the Cartesian coordinate of any point on the surface of paraboloidal reflector. f is the focal length of the paraboloidal reflector and [[rho].sup.2] = [[xi].sup.2] + [[eta].sup.2].

The incident field traveling along z axis is expressed as,

[E.sub.i] = [a.sub.x] exp (-i[k.sub.0]z) (27)

The incident plane wave makes an angle [alpha] with surface normal, where surface normal is given as,

[a.sub.n] = sin [alpha] cos [gamma][a.sub.x] + sin [alpha] sin [gamma][a.sub.y] + cos [alpha][a.sub.z] (28)

where [alpha] and [gamma] are given as,

sin [alpha] = [rho]/[square root of ([[rho].sup.2] + 4[f.sup.2])], (29)

cos [alpha] = 2f/[square root of ([[rho].sup.2] + 4[f.sup.2])], (30)

tan [gamma] = [eta]/[xi] (31)

The wave reflected from paraboloidal reflector is given below,

[E.sub.r] = [E.sub.0] exp[i[k.sub.0] (x sin 2[alpha] cos [gamma] + y sin 2[alpha] sin [gamma] + z cos 2[alpha])] (32)

Rectangular components of field [E.sub.0], at the surface of reflector, are given below,

[E.sub.x0] = [B.sub.[perpendicular to]] [sin.sup.2] [gamma] - [B.sub.[parallel]] [cos.sup.2] [gamma] cos 2[alpha] (33)

[E.sub.y0] = - cos [gamma] sin [gamma]([B.sub.[parallel]] cos 2[alpha] + [B.sub.[perpendicular to]]) (34)

[E.sub.z0] = [B.sub.[parallel]] sin 2[alpha] cos [gamma] (35)

where [B.sub.[parallel]] and [B.sub.[perpendicular to]] are parallel and perpendicular components of GO reflected field from planar layers of chiral and/or chiral nihility layers. The field expression in polar coordinates which is valid around the focal point is given below [40],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (36)

The upper limit of integration with respect to a is taken as,

H = [tan.sup.-1] ([D/2f])

where D is the height of the paraboloidal reflector from horizontal axis.

6. NUMERICAL RESULTS AND DISCUSSION

The field behavior around the focal region of a paraboloidal reflector are obtained by solving the equation, obtained using Maslov's method. Numerical results are presented with the variation of different parameters, i.e., chirality, relative permittivity and thickness of the layer. Chiral nihility is introduced as limiting case of chiral medium where value of relative permittivity and relative permeability of the chiral medium is assumed as [10.sup.-5], keeping nonzero chirality whereas for nihility material, the value of chirality parameter is taken equal to zero. In all the simulations, parameters ([k.sub.0], f, H) are taken as (1, 100, [pi]/4). For all cases, permeability of each chiral layer is taken as unity. All plots deal with behavior of |U(r)| versus kz.

6.1. Chiral-chiral Paraboloidal Reflector

From Figure 15 and Figure 16, it is noted that increasing the permittivity of each layer of the paraboloidal reflector increases the strength of electric field around the focus. Comparison of Figures 15 and 16 also shows that strength of field around focus is more dependent on the permittivity of second layer.

In Figure 17 and Figure 18, effects of variation of thickness of each layer of paraboloidal reflector are presented. It is noted that increase of thickness of first or second layer also increases the strengths of field around focus. It is noticed that first layer has more impact on the field strength as compared to the variation of thickness of the second layer.

Figures 19 and 20 show the effect of change of chirality parameter of each layer while keeping all the other parameters fix. It is observed that increase in chirality of each layer also strengthen the field around focal point of the paraboloidal reflector. Change in chirality of first layer produces greater change in reflected field strength around the focal point as compared to chirality variation for second layer.

6.2. Chiral-PEC and Chiral-nihility Paraboloidal Reflector

In this section two special cases for paraboloidal reflector are discussed one by one: paraboloidal reflector composed of c-PEC and c-n cases. Figure 21 shows reflected field strength for three different values of permittivity of the first layer when second layer is considered as PEC (very large value of permittivity, i.e., [10.sup.5]). It is noted that increase in permittivity of first layer also increases the field strength at the focal point. Figure 22 shows that increase in the chirality of first layer correspondingly increases the field strength at the focus. In Figures 23 and 24, second layer is considered as nihility (by assuming very small values of permittivity and permeability, i.e., 10-5). From Figures 23 and 24, it is observed that by increasing the values of permittivity and chirality, the field strength around the focal point also increases correspondingly.

7. CONCLUSIONS

Characteristics of reflected and transmitted powers from a structure of two parallel layers filled with chiral and/or chiral-nihility metamaterials are studied. Special cases of PEC and nihility backed chiral/chiral nihility layer are also discussed. It is observed that the behavior of reflected or transmitted powers strongly depend on the angle of incidence and values of constitutive parameters describing the metamaterials. For specific values and ranges of incident angles, total reflection of power has been observed hence structure can yield range of Brewster angles. Parallel or perpendicular component of reflected/transmitted power may be selected. That is, either co-polarized or cross polarized field is produced. Similarly for some values of chirality, total reflection/transmission of power is also observed. This arrangement of parallel layers can also be used as a power divider.

Arrangement of two parallel layers mentioned above is utilized to study the behavior of field around the focal region of a large size paraboloidal reflector. Paraboloidal reflector, composed of c-c, c-PEC and c-n layers is discussed for this purpose. PEC and nihility backed cases give higher values of reflected field around the focus as compared to composed of chiral-chiral layers. Nihility material has been introduced in the design of paraboloidal antenna for the purpose of analysis. It is noted that, it provides higher strength around focus as it is a requirement for focusing systems. Use of chiral layers provides an additional parameter, i.e., "chirality", which helps in adjusting the field strength around the focal point. It is concluded that increase in the relative permittivity of each layer also increases field strength around the focus. An increase in the thickness and chirality parameter of the layer also have the same effect, i.e., increase in values of these parameters yields correspondingly higher focusing of the field. It is also noted that different parameters, mentioned above, affect the strength of the field differently. Moreover in the case of c-c paraboloidal reflector, variation of a parameter for first and second layer yields different reflected field strength when all other parameters are kept constant. In this case, change in permittivity of second layer produces more variation in the reflected field strength around the focal point as compared to the first layer, whereas for change in chirality and thickness, field strength of the reflector is more sensitive for the first layer. Study of focusing field strength is of great importance in defence and other electronic systems. This study provides an investigation to design an efficient reflector coated with materials having different material properties. This whole discussion suggests to get control over the focusing of the reflected field by bringing variations in the permittivity, thickness or the chirality of the corresponding layer.

REFERENCES

[1.] Arago, D. F., "Sur une modification remarquable qu' eprouvent les rayons lumineux dans leur passage a travers certains corps diaphanes, et sur quelques autres nouveaux phenomnnes d'optique," Mem. Inst., Vol. 1, 93-134, 1811.

[2.] Biot, J. B., "Phenomenes de polarisation successive, observs dans des fluides homogenes," Bull. Soc. Philomath., 190-192, 1815.

[3.] Fresnel, A., "Memoire sur la double refraction que les rayons lumineux eprouvent en traversant les aiguilles de cristal de roche suivant des directions paralleles a l'axe," Oeuvres, Vol. 1, 731-751, 1822.

[4.] Biot, J. B., "Memoire sur la polarisation circulaire et sur ses applications a la chimie organique," Mem. Acad. Sci., Vol. 13, 39-175, 1835.

[5.] Lindman, K. F., "Ober eine durch ein isotropes system von spiralformigen resonatoren erzeugte rotationspolarisation der elektromagnetischen wellen," Ann. Phys., Vol. 63, 621-644, 1920.

[6.] Lindman, K. F., "Uber die durch ein aktives raumgitter erzeugte rotationspolarisation der elektromagnetischen wellen," Ann. Phys., Vol. 69, 270-284, 1922.

[7.] Jaggard, D. L., A. R. Mickelson, and C. H. Papas, "On electromagnetic waves in chiral media," Appl. Phys., Vol. 18, 211-216, 1979.

[8.] Lindell, I. V., A. H. Sihvola, S. A. Tretyakov, and A. J. Viitanen, Electromagnetic Waves in Chiral and Bi-isotropic Media, Artech House, Boston, 1994.

[9.] Bassiri, S., C. H. Papas, and N. Engheta, "Electromagnetic wave propagation through a dielectric-chiral interface and through a chiral slab," J. Opt. Soc. of Am. A, Vol. 5, 1450-1459, 1988.

[10.] Lakhtakia, A., Beltrami Fields in Chiral Media, World Scientific, Singapore, 1994.

[11.] Lamb, H., "On group-velocity," Proc. London Math. Soc., Vol. 1, 473-479, 1904.

[12.] Schuster, A., An Introduction to the Theory of Optics, Edward Arnold, London, 1904.

[13.] Pocklington, H. C., "Growth of a wave-group when the group velocity is negative," Nature, Vol. 71, 607-608, 1905.

[14.] Malyuzhinets, G. D., "A note on the radiation principle," Zh. Tekh. Fiz., Vol. 21, 940-942, 1951.

[15.] Sivukhin, D. V., "The energy of electromagnetic waves in dispersive media," Opt. Spektrosk., Vol. 3, 308-312, 1957.

[16.] Veselago, V. G., "The electrodynamics of substances with simultaneously negative values of permittivity and permeability," Sov. Phys. Usp., Vol. 10, 509-514, 1968.

[17.] Qiu, C. W., H. Y. Yao, S. Zouhdi, L. W. Li, and M. S. Leong, "On the constitutive relations of G-chiral media and the possibility to realize negative-index media," Microwave Opt. Technol. Lett., Vol. 48, 2534-2538, 2006.

[18.] Tretyakov, S., A. Sihvola, and L. Jylh, "Backward-wave regime and negative refraction in chiral composites," Photonics Nanostruct. Fundam. Appl., Vol. 3, 107-115, 2005.

[19.] Pendry, J. B., A. J. Holden, W. J. Stewart, and I. Youngs, "Extremely low frequency plasmons in metallic mesostructures," Phys. Rev. Lett., Vol. 76, 4773-4776, 1996.

[20.] Pendry, J. B., "A chiral route to negative refraction," Science, Vol. 306, 1353-1355, 2004.

[21.] Qiu, C. W., H. Y. Yao, L. W. Li, S. Zouhdi, and T. S. Yeo, "Routes to left-handed materials by magnetoelectric couplings," Phys. Rev. B, Vol. 75, 245214, 2007.

[22.] Zhang, S., Y. S. Park, J. Li, X. Lu, W. Zhang, and X. Zhang, "Negative refractive index in chiral metamaterials," Phys. Rev. Lett., Vol. 102, 023901, 2009.

[23.] Mackay, T. G. and A. Lakhtakia, "Simultaneous negative- and positive-phase-velocity propagation in an isotropic chiral medium," Microwave Opt. Technol. Lett., Vol. 49, 1245-1246, 2007.

[24.] Lakhtakia, A., "An electromagnetic trinity from negative permittivity and negative permeability," Int. J. Inf. and Mil. Wav., Vol. 22, 1731-1734, 2001.

[25.] Tretyakov, S., I. Nefedov, A. H. Sihvola, S. Maslovki, and C. Simovski, "Waves and energy in chiral nihility," Journal of Electromagnetic Waves and Applications, Vol. 17, No. 5, 695-706, 2003.

[26.] Naqvi, Q. A., "Planar slab of chiral nihility metamaterial backed by fractional dual/PEMC interface," Progress In Electromagnetics Research, Vol. 85, 381-391, 2008.

[27.] Baqir, M. A., A. A. Syed, and Q. A. Naqvi, "Electromagnetic fields in a circular waveguide containing chiral nihility metamaterial," Progress In Electromagnetics Research M, Vol. 16, 85-93, 2011.

[28.] Balanis, C. A., Advanced Engineering Electromagnetics, 2nd edition, John Willey and Sons, 2012.

[29.] Cheng, D. K., Fields and Wave Electromagnetics, Addison- Wesley, New York, 1989.

[30.] Qiu, C. W., N. Burokur, S. Zouhdi, and L. W. Li, "Chiral nihility effects on energy flow in chiral materials," J. Opt. Soc. of Am., Vol. 25, 55-63, 2008.

[31.] Ahmad, F., S. N. Ali, A. A. Syed, and Q. A. Naqvi, "Chiral and/or chiral nihility interfaces: Parametric dependence, power tunneling and rejection," Progress In Electromagnetics Research M, Vol. 23, 167-180, 2012.

[32.] Felson, L. B., Hybrid Formulation of Wave Propagation and Scattering, Nato Science Series E, Martinus Nijho, Dordrecht, Netherlands, 1984.

[33.] Dechamps, G. A., "Ray techniques in electromagnetics," Proc. IEEE, Vol. 60, 1022-1035, 1972.

[34.] Chapman, C. H. and R. Drummond, "Body wave seismograms in inhomogeneous media using Maslov asymptotic theory," Bull. Seismol., Soc. Am., Vol. 72, 277-317, 1982.

[35.] Maslov, V. P., "Perturbation theory and asymptotic method," Gos. Moskov. Univ., Moscow, 1965 (in Russian), Translated into Japanese by Ouchi et al., Iwanami, Tokyo, 1976.

[36.] Ghaffar, A., Q. A. Naqvi, and K. Hongo, "Analysis of the fields in three dimensional Cassegrain system," Progress In Electromagnetics Research, Vol. 72, 215-240, 2007.

[37.] Ji, Y. and K. Hongo, "Analysis of electromagnetic waves refracted by a spherical dielectric interface by Maslov's method," J. Opt. Soc. ofAm. A, Vol. 8, 541-548, 1991.

[38.] Ji, Y. and K. Hongo, "Field in the focal region of a dielectric spherical lens by Maslov's method," J. Opt. Soc. Of Am. A, Vol. 8, 1721-1728, 1991.

[39.] Hongo, K., Y. Ji, and E. Nakajima, "High frequency expression for the field in the caustic region of a reflector using Maslov's method," Radio Sci., Vol. 21, 911-919, 1986.

[40.] Hongo, K. and Y. Ji, "High frequency expression for the field in the caustic region of a cylindrical reflector using Maslov's method," Radio Sci., Vol. 22, 357-366, 1987.

[41.] Hongo, K. and Y. Ji, "Study of the field around the focal region of spherical reflector antenna by Maslovs method," IEEE Trans. Antennas Propagat., Vol. 36, 592-598, 1988.

[42.] Ziolkowski, R. W. and G. A. Deschamps, "Asymptotic evaluation of high frequency field near a caustic: An introduction to Maslov's method," Radio Sci., Vol. 19, 1001-1025, 1984.

[43.] Faryad, M. and Q. A. Naqvi, "High frequency expression for the field in the caustic region of cylindrical reflector placed in chiral medium," Progress In Electromagnetics Research, Vol. 76, 153-182, 2007.

[44.] Faryad, M. and Q. A. Naqvi, "High frequency expression for the field in the caustic region of a parabolic reflector coated with isotropic chiral medium," Journal of Electromagnetic Waves and Applications, Vol. 22, No. 7, 965-986, 2008.

[45.] Rahim, T., M. J. Mughal, Q. A. Naqvi, and M. Faryad, "Focal region field of a paraboloidal reflector coated with isotropic chiral medium," Progress In Electromagnetics Research, Vol. 94, 351-366, 2009.

[46.] Illahi, A. and Q. A. Naqvi, "Study of focusing of electromagnetic waves reflected by a PEMC backed chiral nihility reflector using Maslov's method," Journal of Electromagnetic Waves and Applications, Vol. 23, No. 7, 863-873, 2009.

[47.] Sabah, C. and S. Uckun "Mirrors with chiral slabs," Journal of Optoelectronics and Advanced Materials, Vol. 8, 1918-1924, 2006.

Syed T. H. Shah (1), Faiz Ahmad (2), Nayyar A. Shah (1), Aqeel A. Syed (1), Qaisar A. Naqvi (1) *

(1) Department of Electronics, Quaid-i-Azam University, Islamabad 45320, Pakistan

(2) Department of Electronics, CIIT, Islamabad, Pakistan

Received 7 March 2013, Accepted 3 May 2013, Scheduled 8 May 2013

* Corresponding author: Qaisar Abbas Naqvi (qaisar@qau.edu.pk).

Printer friendly Cite/link Email Feedback | |

Author: | Shah, Syed T.H.; Ahmad, Faiz; Shah, Nayyar A.; Syed, Aqeel A.; Naqvi, Qaisar A. |
---|---|

Publication: | Progress In Electromagnetics Research B |

Geographic Code: | 1USA |

Date: | Jun 1, 2013 |

Words: | 6566 |

Previous Article: | Design of reconfigurable miniaturized UWB-BPF with tuned notched band. |

Next Article: | Improving convergence time of the electromagnetic inverse method based on genetic algorithm using the PZMI and neural network. |

Topics: |