6

Click here to load reader

Tekever drone português espionagem industrial 1.pdf

Embed Size (px)

Citation preview

Page 1: Tekever drone português espionagem industrial 1.pdf

8/17/2019 Tekever drone português espionagem industrial 1.pdf

http://slidepdf.com/reader/full/tekever-drone-portugues-espionagem-industrial-1pdf 1/6

 Fig. 1 A chiral material slab with thickness d is illuminated by a normalincident circularly polarized plane wave. A perfect electric conductorscreen is laced at a distance l from the chiral slab.

Optical Tractor Beam with Chiral LightDavid E. Fernandes, Mário G. Silveirinha 

University of Coimbra, Department of Electrical Engineering – Instituto de Telecomunicações, 3030-290, Coimbra, Portugal,

[email protected]  , [email protected]  

 Abstract—We propose a novel mechanism to promoteoptomechanical interactions that can induce the motion of achiral material body towards the optical source. Our solution isbased on the interference between a chiral light beam and itsreflection on an opaque mirror. Surprisingly, we theoreticallyshow that it is possible to tailor the response of the chiralmetamaterial such that independent of its specific location withthe respect to the mirror, it is always pushed downstream againstthe photon flow associated with the incoming wave. Moreover,we prove that by controlling the handedness of the incomingwave it may be possible to harness the sign of the optical force,and switch from a pushing force to a pulling force

 Index Terms—tractor beam, chirality, metamaterials. 

I.  I NTRODUCTION 

Light transports both energy and momentum and thus theradiation field and matter can exchange momentum [1-3]. This

 property is manifested macroscopically in terms of a radiation pressure that may be used to control the motion and position ofmicrometer-sized particles and neutral atoms [4]. In particular,optical forces have had a tremendous impact on thedevelopment of atom cooling and trapping technologies (e.g.[5]). Moreover, focused laser beams (optical tweezers [6]) arecommonly used in molecule manipulations for biological and

 physical chemistry applications [7-8]. These advances were

made possible by the advent of lasers, which provide themeans to effectively boost the usually weak optomechanicalinteractions.

Recently, there has been a great interest in new techniquesfor optical manipulation of nanoparticles based on non-structured (gradientless) light. Intuitively, one might expectthat when a uniform light beam impinges on a polarizable

 body, the translational force acting on the body should always push it towards the direction of propagation. Notably, severalgroups demonstrated theoretically, and in some cases alsoexperimentally, that tailored travelling light beams may allow

 pulling micron-sized particles towards the light source [9-24].

Such solutions are known as optical “tractor beams”, whichmay be defined as engineered light beams that exert a negativescattering force on a polarizable body forcing it to move alonga direction counter-propagating with the photon flow. A

 pulling force can be obtained with interfering non-diffractive beams [10], based on optical-conveyor belts [9], with

structured media supporting backward waves [12], with thehelp non-paraxial Bessel beams [13, 14, 16, 24], and usinggain media [22, 23].

In this article, we put forward a different paradigm for atractor beam based on unstructured light and on suitablydesigned chiral metamaterials. It is shown that by controllingthe light polarization state helicity it is possible to switch froma pushing force (upstream motion) to a pulling force(downstream motion). The emergence of pulling forcesresulting from the interaction of chiral particles with lightBessel beams has been recently discussed in [24]. Moreover, ithas been shown that light can exert a lateral optical force onchiral particles above a dielectric substrate [25]. As detailed inwhat follows, our proposal to harness the sign of the opticalforce is based on entirely different physical mechanisms, anddoes not require complex Bessel beams.For the sake of clarity, we consider a simple canonicalgeometry wherein the body of interest is a planar slab withthickness d  (Fig. 1). Nevertheless, the ideas that follow can bereadily generalized to other geometries, e.g. to smallnanoparticles. As a starting point, we calculate the force thatlight exerts on the polarizable body when it stands alone in avacuum.

In a stationary time-harmonic regime

(   , Re i t t e    E r E r ) the time-averaged force V F   is

given by ˆ ReV c

dV 

F n T   where cT   is the (complex-

valued) Maxwell stress tensor

This work was funded by Fundação para Ciência e a Tecnologia under 

 project PTDC/EEI-TEL/2764/2012. D. E. Fernandes acknowledges support

 by Fundação para a Ciência e a Tecnologia, Programa Operacional Potencial

Humano/POPH, and the cofinancing of Fundo Social Europeu under the

fellowship SFRH/BD/70893/2010.

Page 2: Tekever drone português espionagem industrial 1.pdf

8/17/2019 Tekever drone português espionagem industrial 1.pdf

http://slidepdf.com/reader/full/tekever-drone-portugues-espionagem-industrial-1pdf 2/6

* * * *0 0

0 0

1 1 1 1

2 2 2c      

T E E E E 1 B B B B 1 , 1  

is the identity tensor, V    is the boundary surface of the bodywith volume V , and n̂   is the outward unit normal vector [26,27]. Suppose that a plane electromagnetic wave illuminates the

 planar slab along the normal direction, creating in this mannerreflected and transmitted transverse plane waves characterized

 by the reflection and transmission coefficients  R  and T . Astraightforward analysis shows that in such a scenario thetranslational force (along the z-direction) acting on the body is

2 2

,av av/ 1inc

 zF A W R T   , where  A  is the area of the slab

cross-section and   2 2

av 0 0

1

4inc inc incW      E H   is the time-

averaged energy density of the incident plane wave. It isinteresting to note that ,av z

F   is equal to the difference between

the time rates of the electromagnetic momentum entering theslab front face and the momentum exiting the back face.Clearly, for passive materials the conservation of energy

implies that the optical force is nonnegative, ,av 0 zF    , suchthat the body is dragged upstream by the light wave. Notably,for active materials the transmission coefficient can be greater

than unity and hence it is possible to have2 2

1 0 R T   

and a pulling (negative) force [22, 23].Here, however, we are interested in exploring a differentsolution based on the interference of the upstream wave with adownstream wave created by reflection on an opaque mirror(Fig. 1). It is well known that the interference of the two wavescreates a standing wave, and that usually a material body is

 pushed to a spot wherein the electric field intensity reaches alocal maximum. Surprisingly, we prove in what follows that it

is possible to engineer the material response such thatindependent of the location  of the material body it is always

 pulled downstream towards the light source. Thus, in such asystem the light beam behaves as an optical tractor beam.As seen previously, a light wave propagating upstream(downstream) tends to drag the material body along the up(down) direction. Clearly, in the absence of a mirror, theoptical force is minimal ( ,av 0 zF    ) when the material body is

transparent to the incoming wave ( 0 R     and 1T   ). This

suggests that to have a negative optical force the upstreamwave should desirably be perfectly transmitted through thematerial, so that no momentum is imparted on the material

 body by this wave. On the other hand, the downstream wave isrequired to be either absorbed or reflected by the material bodyso that it can transfer a negative momentum to the body, andinduce the negative optical force. Thus, the material body mustscatter the waves propagating upstream and downstream in anasymmetric  manner. Clearly, this cannot be achieved with astandard dielectric slab because by symmetry thecorresponding scattering response is the same for the twodirections. Moreover, a stack of conventional dielectrics is alsonot a valid option because if 1T     for the upstream wave,

then, because of the Lorentz reciprocity theorem [27], T    is

also required to be equal to unity for the downstream wave andhence the optical force vanishes.

 Notably, chiral metamaterials provide a way out of this bottleneck, and enable us to engineer the scattering propertiessuch that the transmission of the upstream and downstreamwaves is strongly asymmetric. Chiral media are a special class

of isotropic media whose electromagnetic response depends onthe handedness of its structural unities [28], and ischaracterized by a magneto-electric coupling that breaks thedegeneracy between two circularly polarized plane waves sothat the refractive index increases for one circular polarizationand decreases for the other [28]. Chiral materials endow uswith the means to achieve polarization rotation (opticalactivity) [29]-[31], polarization conversion [32], and negativerefraction [33]-[36]. The constitutive relations of a chiral

material are 0 0 0i   D E H   and

0 0 0i    B E H   where      is the dimensionless

chirality parameter,    is the relative permittivity, and    is the

relative permeability. The photonic modes in an unboundedchiral medium are right circular polarized (RCP) waves or leftcircular polarized (LCP) waves characterized by distinct

refractive indices n      , where n   stands for the

refractive index of RCP waves and n  for the LCP waves [28].

It is well known that a chiral slab scatters RCP and LCP wavesdifferently, such that the transmission and reflectioncoefficients for an incident wave are [28], [37]:

02 2 2

4

1 1

ik d 

ink d  

 zeT 

 z z e

 

,

0

0

22

2 2 2

1 1

1 1

ink d  

ink d  

 z e R R

 z z e

. (1)

In the above, d   is the thickness of the chiral slab,

0k k n      , 0 /k c  , n     , and / / z n    

is the normalized impedance. Hence, the reflection coefficientis identical for the two polarization states but, crucially, thetransmission coefficient is polarization dependent. Now, thekey observation is that in presence of an opaque mirror, e.g. fora metallic mirror, the wave polarization state is switched uponreflection in the mirror, such that if the upstream wave has

RCP (LCP) polarization then the downstream wave has LCP(RCP) polarization. This implies that the scattering of theupstream and downstream waves can be asymmetric becausethe transmission coefficient for RCP waves ( T  ) differs from

the transmission coefficient for LCP waves ( T  ). It must be

highlighted that the asymmetric transmission ofelectromagnetic waves through chiral metamaterials wasdiscussed by different groups in other contexts [38, 39, 40].Let us denote 21T    as the transmission coefficient for a

circularly polarized wave propagating upstream that impinges

Page 3: Tekever drone português espionagem industrial 1.pdf

8/17/2019 Tekever drone português espionagem industrial 1.pdf

http://slidepdf.com/reader/full/tekever-drone-portugues-espionagem-industrial-1pdf 3/6

on a chiral slab standing alone in free-space ( 21T T    for an

incident RCP/LCP wave; the subscript “21” indicates that theincident wave propagates from the region 1 to the region 2),and 12T   as the transmission coefficient for the corresponding

downstream wave ( 12T T  

). One important remark is that the

conservation of energy and the fact that  R  is polarizationindependent imply that

12 21

T T    in the absence of material

absorption. Thus, a strong asymmetric transmission requiressome material loss.To put the previous ideas on a firm ground, we calculated theoptical force acting on a chiral material slab when it stands at adistance l from the opaque mirror. For simplicity, the mirror isassumed to be a perfect electric conductor (PEC), but similarresults can be obtained with other mirror types (e.g. with a

 photonic crystal wall). For the geometry of Fig. 1, theelectromagnetic fields are of the form:

 z C z C z

E e ,

0  z C z C z i  

H e   (2)

with ˆ ˆ / 2i  e x y . The signal   determines if the

incoming (upstream) plane wave has RCP (+) or LCP polarization (). The function C z  (   C z ) represents the

complex amplitude of the upstream (downstream) wave ineach point of space. Using the Maxwell-stress tensor, it can bereadily found that the translational force acting on the chiralslab is simply

  2 2

.av 0/ / 4out 

in

 z z

 z z z

F A C z C z C z C z 

, where in z z   ( out  z z ) are the coordinates of the input(output) interfaces of the slab. The optical force can also bewritten as:

  2 2 2 2

.av 0

2 z

out out in in

F C C C C  

 A

      (3)

where in inC C z  and out out  

C C z . Note that inC   is the

complex amplitude of the incident electric field at the chiralslab input interface. The coefficients in

C   and out C    are linked

 by 12in in out  C RC T C    and 21out in out  C T C RC   , where 12T  ,

21T  , and  R   are the transmission/reflection coefficients

introduced before for a slab standing alone in free-space.Moreover, for a PEC mirror the coefficients out 

C   are bound to

obey 02i k l

out out  C e C  . Putting all these results together, it is

 possible to express inC   and out 

C   as a function of the incident

wave amplitude inC   as follows:

0

0

221 12

21

i k l

in ini k l

e T T C R C 

 Re

,

0

2121out ini k l

T C C 

 Re

,

0

0

221

21

i k l

out ini k l

e T C C 

 Re

. (4)

Substituting this result into Eq. (3), we finally find that the

optical force is:0

0 0

2 22.av 21 12 21

av 2 21 2

1 1

i k l

inc z

i k l i k l

F  e T T T  W R

 A  Re Re

, (5)

where2

0av 2inc

inW C       is the time-averaged energy density

of the incident plane wave.As previously mentioned, in order to minimize the

upstream light-matter interactions, it is desirable that 0 R    

and 21 1T      so that the slab is transparent to the upstream

wave. In this ideal case, the optical force becomes 2

.av av 12 1 0inc

 zF W T  , i.e. consistent with our intuition the

optomechanical interactions force the material slab to movedownstream towards the light source, independent of itslocation with respect to the mirror. Thus, the proposed regimeis fundamentally different from the conventional light-matterinteractions in presence of standing waves, which tend to draga polarizable object to a field intensity maximum.Interestingly, the amplitude of the negative force has amaximum value when 12 0T    . A simple picture of this

scenario is that first the upstream beam overtakes the chiral

slab without any momentum transfer, and then, after reflectionin the mirror, it is totally absorbed by the chiral slab,transferring in this way its electromagnetic momentum to theslab. Clearly, this regime requires a strongly asymmetricscattering by the material slab (ideally, 0 R   , 21 1T    ,

12 0T    ). Given that the incoming wave has a positive

momentum, one may wonder how the material slab can gain anegative momentum. The justification is simple: the negativemomentum is indirectly provided by the mirror wall whichupon reflection flips the sign of the incident radiationmomentum.

A remarkable property of our system is that if the

handedness of the upstream light is reversed (e.g. from LCP toRCP) then 21T   and 12T   are interchanged. Hence, if for a certain

light polarization the upstream light wave originates a negativeoptical force (e.g. 0 R   , 21 1T    , 12 0T    ), then an

upstream light wave with opposite handedness originates a positive optical force (e.g. 0 R   , 21 0T    , 12 1T    ). Indeed,

in the latter situation the upstream light is completely absorbed by the chiral slab never reaching the opaque mirror, and thus itimparts a positive momentum to the material body. Therefore,

 by controlling the handedness of the upstream light it may be

Page 4: Tekever drone português espionagem industrial 1.pdf

8/17/2019 Tekever drone português espionagem industrial 1.pdf

http://slidepdf.com/reader/full/tekever-drone-portugues-espionagem-industrial-1pdf 4/6

 

Fig. 2 (a) Geometry of the conjugated-gammadions planar chiral metamaterial. The length of the central arm of the gammadions is denoted by l  and the width by w . The meta-atom is arranged in a periodic square lattice with period a. (b) Side view of the meta-atom: the conjugated-gammadions are made of silver

with thicknessmd   and are separated by a distance

gd  . The two gammadions are embedded in a polyimide material with thickness 2l  and electric permittivity

0 6.25 1 0.03 p i    . Two layers of silicon, with permittivity 0 11.9 1 0.004s i     and thickness 1l , are placed on each side of the polyimide material.

(c) Amplitude of: i) Transmission coefficient for incident LCP waves, ii) Transmission coefficient for incident RCP waves, iii) Reflection coefficient forincident LCP waves, iv) Reflection coefficient for incident RCP waves (d) Similar to (c) but for the phase of the transmission and reflection coefficients. (e)Circular dichroism.

 possible to transport a material particle at will, either upstreamor downstream, similar to an optical conveyer belt.

The asymmetry in the response is determined by2 2

T T  ,

which is a measure of the circular dichroism. Therefore, achiral material that enables the light source to behave as atractor beam must have a strong circular dichroism.To demonstrate these ideas, we designed a chiral metamaterial

such that the meta-atoms are made of conjugated-gammadionsas illustrated in Fig. 2a-b. Previous studies have shown thatthis metamaterial is characterized by a strong circulardichroism and large optical activity in the microwave and low-THz regime [30],[31]. Evidently the response of thismetamaterial is not isotropic, but the key requirement to have a

 pulling force is the strong circular dichroism, rather than theisotropy.To take advantage of high-power laser sources, the meta-atomwas engineered to operate around the wavelength

0 1.55μm   . The length of the central arm of the

gammadions is 151.81nml  , the width is 12.19 nmw  andits thickness is 25nmmd   . The conjugated-gammadions are

made of silver, described by a Drude dispersion model

21m p i   , with parameters

2 2175 z p     and 2 4.35 z  , consistent with

experimental data reported in the literature [41]. Thegammadions are separated by a distance 11.58 nmg

d   and

embedded in a polyimide material with thickness 2 88.4nml    

and electric permittivity 0 6.25 1 0.03 p i    [30]. We also

added two layers of silicon (   0 11.9 1 0.004s i    , one on

each side of the polyimide material, as shown in Figure 2b,with thickness 1 50.6 nml   , so that the total thickness of the

meta-atom is 1 2 02 189.6 nm 0.122d l l     . The meta-

Page 5: Tekever drone português espionagem industrial 1.pdf

8/17/2019 Tekever drone português espionagem industrial 1.pdf

http://slidepdf.com/reader/full/tekever-drone-portugues-espionagem-industrial-1pdf 5/6

 Fig. 3 Normalized force .av av

inc

 zF AW    as a function of the frequency for

two fixed positions of the chiral slab: 0 1.5k l    (solid curves) and

0 2.0k l  (dashed curves) for i) RCP upstream light and ii) LCP upstream

light.

atoms are arranged in a periodic square lattice with a deeplysubwavelength period 0202.41nm 0.131a     .

Using the full-wave electromagnetic simulator CST-MWS [42]we calculated the transmission and reflection matrices thatrelate the fields scattered by a metamaterial slab standing alonein free-space with the incident fields as

t inc xx xy x x

t inc yx yy y y

T T  E E 

T T  E E 

  and

r inc xx xy x x

r inc yx yy y y

 R R E E 

 R R E E 

.

Here, inc

i E  , t 

i E    and r 

i E    are the incident, transmitted and

reflected ,i x y   field components, respectively. The

calculation of the reflection and transmission matrices requiresilluminating the metamaterial slab with two independentexcitations, in our case with x- and y- linearly polarized waves.The rotation symmetry of the conjugated-gammadions causesthe cross polarized components of the matrices to besymmetric, so that  xy yxT T   and  xy yx R R , whereas the co-

 polarized components are equal, i.e.  xx yyT T    and  xx yy R R .

As a consequence, the eigenstates of the scattering problem are

circularly polarized waves, with eigenvectors e . Moreover,the eigenvalues of the transmission and reflection matrices are

 xx xyT T iT     and  xx xy R R iR  , respectively [35].

The numerically calculated amplitude and phase of thetransmission and reflection coefficients for incident circularly polarized light are shown in Figure 2c-d. The resultsdemonstrate that, as expected, the reflection coefficients forRCP and LCP waves are the same. Most important, thesimulations show that the chiral metamaterial has a giantoptical activity near 160THz f   , 196THz f   and

245THz f   , where the scattering of the RCP and LCP

waves is strongly asymmetric. This is highlighted in Fig. 2ewhich depicts the circular dichroism of the structure. Notably,near 196THz f   there is a spectral region wherein

2 2T T    is nearly unity while the reflection coefficient is

near zero. These are exactly the conditions required to have a pulling force with an LCP excitation. Figure 3 depicts thenormalized optical force .av av

inc

 zF AW    calculated for two

 particular positions of the slab, 0 1.5k l    (solid curves) and

0 2.0k l  (dashed curves) for an upstream wave with either

RCP or LCP polarization. The results confirm that only in afrequency range around 196THz f   , the force exerted by an

LCP excitation is negative for the two positions of the PECscreen. Moreover we see that in the same frequency range anupstream RCP wave exerts a positive force on the chiralmetamaterial, in agreement with an earlier theoreticaldiscussion.

Figure 4a and 4c depict the optical forces created by RCPand LCP upstream excitations at two fixedfrequencies 190THz f   and 196THz f   , respectively. As

seen, while in the former case –due to the interference of theupstream and downstream waves– the optical force signdepends on the specific position of the slab relative to the

mirror, in the latter case an LCP excitation always exerts anegative force on the slab whereas an RCP excitation always pushes the slab towards the mirror, independent of the distancel . Thus, the LCP excitation mimics, indeed, an optical tractor

 beam. Here, it is interesting to note that the force varies periodically with the distance between the chiral slab and themetallic screen, with a period 0 / 2l     , where 0    is the

wavelength of light in vacuum.

Because the optical force only depends on the position of theslab relative to the mirror, we can define a potential energy as

,av zV F dz , which is depicted in Fig. 4b and 4d for

190THz f   and 196THz f   , respectively. As seen, in the

former case there are potential wells wherein the material body

will inevitably become “trapped”. Quite differently, at196THz the potential energy has a monotonic dependencewith l and does not have stationary points. Hence, dependingon the light polarization, the material body is steadily pushedeither upstream (RCP polarization) or downstream (LCP polarization). Interestingly, this optical “conveyer belt”operation may occur in a spectral range from

194.419 THz f   to 198.57THz f   , which corresponds to a

 bandwidth of 4.151THz . Finally, we note that a typical laseroperating at 0 1.55μm   generates a light beam with a

Poynting vector 289.3MW mincS    [43]. In this case, an

optical force with amplitude .av av 1.0inc

 zF AW     corresponds to

an optical pressure of the order .av / 0.3 Painc

 zF A S c ,

which is a significant value in the nanoscale.In summary, it was theoretically demonstrated that by usingcircularly polarized light and a planar chiral metamaterial witha tailored response it may be possible to mimic an opticaltractor beam. By controlling the incident light polarization onemay switch from a pushing force to a pulling force,independent of the specific position of chiral slab. We envisionthat these findings may lead into new inroads in opticalmanipulation of micro and nanoparticles.

Page 6: Tekever drone português espionagem industrial 1.pdf

8/17/2019 Tekever drone português espionagem industrial 1.pdf

http://slidepdf.com/reader/full/tekever-drone-portugues-espionagem-industrial-1pdf 6/6

 

Fig. 4 (a) Normalized force .av avinc

 zF AW   and (b) Potential (in arbitrary units) as a function of the normalized distance 0k l  between the dielectric slab and the

PEC screen at a fixed frequency 190 THz f   for i) RCP upstream light and ii) LCP upstream light. (c) and (d) Similar to (a) and (b), respectively, but for afixed fre uenc 196 THz f   .

R EFERENCES [1] J. C. Maxwell, Treatise on Electricity and Magnetism, DoverPublications, 3rd Ed., 1954.

[2] P. Lebedev, Ann. der Physik, 6, 433 (1901). 

[3] E. F. Nichols and G. F. Hull, Phys. Rev. 13, 307 (1901).

[4] A. Ashkin, Phys. Rev. Lett. 24, 156 (1970); A. Ashkin, Phys. Rev. Lett. 25,1321 (1970).

[5] S. Chu. Rev. Mod. Phys. 70, 685 (1998).

[6] A. Ashkin, J. M. Dziedzic, J. E. Bjorkholm, and S. Chu, Opt. Lett. 11, 288(1986).

[7] O. M. Maragò, P. H. Jones, P. G. Gucciardi, G. Volpe and A. C. Ferrari, Nat. Nanotechnol. 8, 807 (2013).

[8] K. Dholakia, T. Čižmár, Nature Photon. 5, 335 (2011).

[9] T. Čižmár, V. Garcés-Chávez, K. Dhokalia, and P. Zemánek,  Appl. Phys. Lett. 86, 174101 (2005).

[10] S. Sukhov and A. Dogariu, Opt. Lett. 35, 3847 (2010).

[11] S.-H. Lee, Y. Roichman and D. G. Grier, Opt. Expr. 18, 6988 (2010).

[12] A. Salandrino, D. N. Christodoulides, Opt. Lett. 36, 3103, (2011).

[13] J. Chen, J. Ng, Z. Lin and C. T. Chan, Nat. Phot. 5, 531 (2011).

[14] A. Novitsky, C.-. Qiu, H. Wang, Phys. Rev. Lett., 107, 203601 (2011).

[15] S. Sukhov and A. Dogariu, Phys. Rev. Lett. 107, 203602 (2011). 

[16] A. Novitsky, C.-W. Qiu, A. Lavrinenko, Phys. Rev. Lett., 109, 023902,(2012).

[17] D. B. Ruffner and D. G. Grier, Phys. Rev. Lett. 109, 163903 (2012).

[18] O. Brzobohatý, V. Karásek, M. Šiler, L. Chvátal, T. Čižmár and P.Zemánek, Nat. Phot. 7, 123 (2013).

[19] V. Kajorndejnukul, W. Ding, S. Sukhov, C.-W. Qiu, A. Dogariu,  Nat.Phot. 7, 787 (2013). 

[20] N. Wang, W. Lu, J. Ng, Z. Lin, Opt. Lett., 39, 2399, (2014).

[21] C. E. M. Démoré, P. M. Dahl, Z. Yang, P. Glynne-Jones, A. Melzer, S.Cochran, M. P. MacDonald, G. C. Spalding, Phys. Rev. Lett.  112, 174302(2014).

[22] A. Mizrahi, Y. Fainman, Opt. Lett. 35, 3405 (2010).

[23] K. J. Webb, Shivanand, Phys. Rev. E  84, 057602 (2011).

[24] K. Ding,  J. Ng,  L. Zhou  and C. T. Chan,  Phys. Rev. A.  89, 063825(2014).

[25] S.B. Wang, C.T. Chan, Nat. Commun., 5, 3307, (2014).

[26] L. D. Landau, E. M. Lifshitz, and L. P. Pitaevskii,  Electrodynamics ofContinuous Media, Butterworth-Heinemann, 1985 (Oxford).

[27]  J. D. Jackson, Classical Electrodynamics (3rd  Ed.), John Wiley, 1999(New York, NY).

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

[29] M. Decker, M. Ruther, C. E. Kriegler, J. Zhou, C. M. Soukoulis, S.Linden, M. Wegener, Opt. Lett. 34, 2501 (2009).

[30] R. Zhao, J. Zhou, T. Koschny, E. N. Economou and C. M. Soukoulis,Phys. Rev. B 81, 235126 (2010).

[31] R. Zhao, L. Zhang, J. Zhou, T. Koschny, and C. M. Soukoulis, Phys. Rev. B 83, 035105 (2011).

[32] M. G. Silveirinha,  IEEE Trans. on Antennas and Propag., 56, 390,(2008). 

[33] J. B. Pendry, Science 306, 1353 (2004).

[34] S. A. Tretyakov, A. V. Sihvola and L. Jylhä, Phot. Nano. Fund. Appl. 3,107 (2005).

[35] J. Zhou, J. Dong, B. Wang, T. Koschny, M. Kafesaki, C. M. Soukoulis,Phys. Rev. B 79, 121104(R) (2009).

[36] T. A. Morgado, S. I. Maslovski, M. G. Silveirinha,  New J. of Phys. , 14,063002, (2012).

[37] R. Zhao, T. Koschny, C. M. Soukoulis, Opt. Express, 18, 14553, (2010).

[38] V. A. Fedotov, P. L. Mladyonov, S. L. Prosvirnin, A. V. Rogacheva, Y.Chen, and N. I. Zheludev, Phys. Rev. Lett., 97, 167401, (2006).

[39] R. Singh, E. Plum, C. Menzel, C. Rockstuhl, A. K. Azad, R. A. Cheville,F. Lederer, W. Zhang, N. I. Zheludev, Phys. Rev. B, 80, 153104, (2009).

[40] C. Menzel, C. Helgert, C. Rockstuhl, E.-B. Kley, A. Tünnermann, T.Pertsch, and F. Lederer, Phys. Rev. Lett., 104, 253902, (2010). 

[41] M. A. Ordal, R. J. Bell, R. W: Alexander, Jr., L. L. Long, and M. R.Query, Appl. Opt. 24, 4493 (1985). 

[42] CST GmbH 2013 CST Microwave Studio http://www.cst.com.

[43] G. P. Agrawal,  Nonlinear Fiber Optics, Academic Press, 2001(SanDiego, Calif.), (Chap. 4.)