ABSTRACT

Electromagnetic wave propagation is first analyzed in a composite material mde of chiral nano-inclusions embedded in a dielectric, with the help of Maxwell-Garnett formula for permittivity and permeability and its reciprocal for chirality. Then, this composite material appears as an homo-geneous isotropic chiral medium which may be described by the Post constitutive relations. We analyze the propagation of an harmonic plane wave in such a medium and we show that two different modes can propagate. We also discuss harmonic plane wave scattering on a semi-infinite chiral composite medium. Then, still in the frame of Maxwell-Garnett theory, the propagation of TE and TM fields is investigated in a periodic material made of nano dots immersed in a dielectric. The periodic fields are solutions of a Mathieu equation and such a material behaves as a diffraction grating.

1. Introduction

Nanotechnology is blossoming with in particular the inclusion of nano-particles (nano-dots) in some specific support [1,2]. Then, to analyze electromagnetic wave propagation such as light or X rays in these composite materials, we need a theory able to calculate average, macroscopic values from their granular microscopic properties. This job is performed by the Maxwell-Garett theory [3-6].

In this work, the support is a dielectric with permittivity e₁, permeability m and we consider two situations according that nano-dots are chiral or periodically distributed along a direction of the structure.

In the first case (chiral nano-dots) the permittivity e of the composite material is according to the MaxwellGarett formula

[e – e₁][e + 2e₁]^-1 = f [e₂ – e₁][e₂ + 2e₁]^-1(1)

f is the filling factor of inclusions (their volume fraction) in the host material, the subscripts 1, 2 corresponding to host and inclusions respectively and we get from (1)

e  = e₁ (1 + 2af )(1 - af)^-1a = (e₁ - e₂)(e₁ + 2e₂)^-1,(2)

Permeability µ is assumed the same for nano-dots and dielectric.

The relation (2) has been generalized [7,8] to chirality x when both inclusions and host materials are chiral. But here, the situation is different since only inclusions have this property and the relation (1) with x, x₁, x₂ instead of e, e₁, e₂ has no meaning when x₁ = 0. To cope with this difficulty, we introduce a reciprocal Maxwell-Garnett relation obtained by applying to (1) the transformation (x₁, x₂, f) Þ (fx₂, x₁, 1/f) which gives

[x – fx₂][x + 2fx₂]^-1 = f ^-1[x₁ – fx₂][x₁ + 2fx₂]^-1 (3)

reducing for x₁ = 0 to x = -2fx₂(1 - f)(1 + 2f)^-1 (4)

From now on, we assume f << 1, e₁ > 0, e₂ < 0, a > 0 and x₂ < 0 so that the 0(f ²) approximation of (2) and (4) gives e = e₁(1 + 3af ) > 0 (5)

x = -2fx₂ > 0(6)

So, this composite material made of nano-chiral particles included in a dielectic may be hand-led as an homogeneous chiral medium with permittivity and chirality (5) and (6) and permeabiity µ > 0 assumed to be the same for inclusions and dielectric.

In the second situation (periodically distributed nanorods), the relation (1) is still valid with f changed into a periodic function f (x). Assuming f (x) = f cos(2ax), we write the permittivity e (x) in the following form reducing to (5) to the 0(f ²) order e (x) = e₁exp[3afcos(2ax)](7)

Using (5)-(7), we shall analyze harmonic plane wave propagation in both composite materials, chiral and periodic.

2. Harmonic Plane Wave Propagation in a Chiral Composite Medium

We suppose this chiral medium endowed with the Post constitutive relations in which e, x have the expressions (5) and (6) [9,10]

D = e E + ix B, H = B/µ + ix E, i = (8)

This choice is not arbitrary because the Post constitutive relations, in their general form, are covariant under the proper Lorentz group as Maxwell’s equations which guarantees a consistent theory with a simple mathematical formalism, in agreement with the statement that only covariant mathematical expressions have a physical meaning.

Plane wave scattering from a semi-infinite chiral medium was discussed some time ago by Bassiri et al [11], also using the Post constitutive relations, but we proceed differently from these authors working with the Fresnel reflection and transmission amplitudes.

2.1 Refractive Indices

We consider harmonic plane waves with amplitudes E, B, D, H

(E,B,D,H)(x,t) = (E,B,D,H)y(x,t)(9)

in which y (x,t) = exp[iw(t + nsinθ x/c + ncosθ z/c)] (10)

in which n is a refractive index to be determined.

Substituting (9) into the Maxwell equations

ÑÙE + 1/cθ _t B = 0, Ñ.B = 0

ÑÙH - 1/cθ _t D = 0, Ñ.D = 0(11)

and taking into account (10) give the equations for the amplitudes E, B, D, H

-ncosθ E_y + B_x = 0,

n(cosθ E_x - sinθ E_z) + B_y = 0,

nsinθ E_y + B_z = 0, (12)

ncosθ H_y + D_x = 0

n(cosθ H_x - sinθ H_z) - D_y = 0

nsinθ H_y - D_z = 0(13)

with the divergence equations sinq B_x + cosq B_z = 0, sinq D_x + cosq D_z = 0(14)

We get at once from (8) and (14), the divergence equation satisfied by the electric field sinq E_x + cosq E_z = 0(15)

Substituting (8) into (13) gives ncosq (B_y /m + ixE_y) + eE_x + ix B_x = 0 ncosq (B_x /m + ixE_x) - nsinq (B_z /m + ixE_z) - eE_y - ixB_y = 0 nsinq (B_y /m + ixE_y) - eE_z -ixB_z = 0 (16)

Taking into account (12), these equations become ncosq B_y /m + eE_x + 2ixB_x = 0 ncosq B_x /m - nsinq B_z/m - eE_y - 2ixB_y = 0 nsinq B_y /m - eE_z - 2ixB_z = 0(17)

Then, eliminating B between (12) and (17) gives the homogeneous system of equations in which s = 2nx

(n²/m - e)E_x - iscosq E_y = 0

(n²/m - e)E_y - is(sinq E_z - cosq E_x) = 0

(n²/m - e)E_z + issinq E_y = 0(18)

This homogeneous system has nontrivial solutions if its determinant is null and a simple cal-culation gives

(n²/m - e)[(n²/m - e)² -  s²] = 0(19)

Deleting (n²/m - e) = 0 which would correspond to an a-chiral medium, we get from (11) two modes (n_±²/m - e) = ± s in which which s = 2nx so that the refractive index depends not only on permittivity and permeability but also on chirality with the positive expressions n₊ = xm + (x²m² + em)^1/2, (20)

n_- = -xm + (x²m² + em)^1/2 (21)

Changing the square root into its opposite gives negative refractive indices.

Consequently, two modes with respectively the refractive indices n₊, n_- can propagate in the metachiral slab, they are independent as long as the medium is infinite, otherwise they become coupled at boundaries. The amplitudes of the field components in these two modes have now to be determined.

2.2 Electromagnetic Fields

1) We first suppose n₊²/m - e = s and n₊ = xm + (x²m² + em)^1/2 with e and µ > 0: fields and parameters are characterized by superscripts or subscripts + respectively.

Then, we get at once from (18) and (12) in terms of E⁺_y

E ⁺_x = icosq₊ E⁺_y, E⁺_z = - isinq₊ E⁺_yB⁺_x = n₊ cosq₊ E⁺_yB⁺_y = -i n₊ E⁺_y, B⁺_z = - n₊ sin₊ E⁺_y (22)

and substituting (22) into (8)

 D⁺_x = icosq₊ l₊E⁺_yD⁺_y = l₊E⁺_yD⁺_z = -isinq₊ l₊E⁺_y (23)

H⁺_x = cosq₊ n₊E⁺_yH⁺_y = -i n₊E⁺_yH⁺_z = -sinq ₊ n₊E⁺_y (24)

in which l₊ = e + xn₊, n₊ = n₊/m -  x = (x ² + e /m)^1/2 (25)

2) For n_-²/m - e = - s and n_- = - xm + (x ²m ² + em)^1/2, we get at once with now super-scripts and subscripts -:

    E^-_x = - i cosq_- E_y, E^-_z = i sinq_- E^-_yB^-_x = n_- cosq_- E^-_y, B^-_y = i n_-E^-_yB^-_z = - n_- sinq_- E^-_y(26)

and substituting (26) into (8)

D^-_x = -icosq_- l_-E^-_y, D^-_y = l_-E^-_y,

D^-_z = isinq_- l_-E^-_y(27)

H^-_x = cosq_- n_-E^-_y, H^-_y = i n_-E^-_y,

H^-_z = - sinq_- n_-E^-_y (28)

with l_- = e - xn_-, n_- = n_-/m +  x = (x ² + e/m)^1/2 = n₊ (29)

Then, according to (9) and (10), the electromagnetic field of the plus and minus modes, each depending on an arbitrary amplitude E⁺_y , E^-_y, is

(E^±, B^±, D^±, H^±) (x,t) = (E^±, B^±, D^±, H^±) y_±(x, t) (30)

with the amplitudes given by (22)-(24) and (26)-(28) and the phase functions y_±(x,t) = exp[iw (t + n_± sinq_± x/c + n_± cosq_± z/c)](31)

2.3 Plane Wave Scattering from a Semi-Infinite Chiral Composite Medium

We suppose that the chiral composite material fulfills the half space z < 0 on which impinges from z > 0 on the interface z = 0 an harmonic plane wave characterized by the phase factor y(q_i)

y (q_i) = exp[-iwn₀(xsinq_i + zcosq_i)](32)

n₀ is the refractive index in z > 0 and the components of the incident electromagnetic field are [12] with two amplitudes M_i, N_i:

Eⁱ_x = -cosq_i M_i y (q_i), Eⁱ_y = N_i y (q_i), Eⁱ_z = sinq_i M_i y (q_i)

Hⁱ_x = -n₀cosq_i N_i y (q_i), Hⁱ_y = -n₀M_i y (q_i)Hⁱ_z = n₀sinq_i N_i y(q_i) (33)

The reflected field in the half-space z > 0 has a similar expression with (M_i, N_i, q_i) changed into (M_r, N_r, q_r) while the refracted field in z < 0 is supplied by (30).

According to (31) and (32), also valid for the reflected wave, the continuity of the phase at z = 0 implies the Descartes-Snell relations n₀ sinq_i = n₀ sinq_r = n₊ sinq₊ = n_- sinq_- (34)

The continuity of the components E_x,y, H_x,y, at z = 0 supplies four boundary conditions to de-termine in terms of M_i, N_i the amplitudes M_r, N_r of the reflected field and those E ⁺_y, E ^-_y of the refracted field.

According to (22), (26) and (33) and taking into account (34), we get for the E_x,y components cosq _i(M_r - M_i ) = icosq ₊ E ⁺_y - icosq_- E ^-_y

N_r + N_i = E ⁺_y + E ^-_y(35)

while for H_x,y, according to (24), (28) and (33), we have since n_- = n₊ (= n )

n₀ cosq_i(N_r - N_i ) = n (cosq₊ E ⁺_y + icosq_- E ^-_y)

n₀( M_r + M_i) = n(E ⁺_y - E ^-_y )(36)

To make calculations easier, we introduce the notations M_r + M_i = M, N_r + N_i = NM_r - M_i = M’, N_r - N_i = N’(37)

and a = n₀/n (cosq₊ + cosq_-)^-1 (38)

Then, we get at once from (36)

E ⁺_y = a(cosq_i N’ + cosq_- M)

E ^-_y = a(cosq_i N’ - cosq₊ M) (39)

and, substituting (39) into (35) gives cosq_i M’ = a₁₁ N’ + a₁₂ M N = a₂₁ N’ + a₂₂ M(40)

in which a₁₁ = iacosq_i (cosq ₊ + cosq_-), a₁₂ = 2iacosq ₊ cosq_-

a₂₁ = acosq _i, a₂₂ = a(cosq _- + cosq ₊ ) (41)

Taking into account (37) the system (40) becomes

(cosq _i - a₁₂)M_r + a₁₁ N_r = (cosq _i + a₁₂) M_i - a₁₁ N_i

a₂₂M_r - (1 - a₂₁) N_r = -a₂₂M_i + (1 + a₂₁)N_i(42)

from which we easily get the amplitudes M_r, N_r of the reflected field and consequently M’, N’ according to (37) to obtain finally the amplitudes E_y^± of the refracted field from (39).

One has a simple result for a normal incidence q_i = q_r = q_± = 0 since the Equations (35) and (36) reduce to M_r - M_i = i(E ⁺_y - E ^-_y), N_r + N_i = E ⁺_y + E ^-_y

-2n₀ N_i = n (E ⁺_y + E ^-_y), n₀(M_r + M_i) = n (E ⁺_y - E ^-_y) (43)

with the solution M_r = (n + in₀) (n - in₀)^-1 M_i, N_r = -(1 + 2n₀ /n)N_i (44)

E ⁺_y = n₀(n - in₀)M_i - n₀/n N_i,

E ^-_y = -n₀(n - in₀)M_i - n₀ /n N_i (45)

Remark 1. If the angles q₊, q_-  obtained from (34) are real, the plus and minus modes propa-gate in the chiral medium. If they are both purely imaginary, we get from (34)

cos(q_±) = -i[(n₀/n_±)² sin²q_i - 1]^1/2 (46)

the negative sign in front of the square root in (46) corresponds to the physical situation: refracted waves are evanescent and, incident waves undergo a total reflection, with as consequence for beams of plane waves a GoösHanken lateral shift and a Imbert-Fedorov transverse shift [13]. Of course with a single angle pure imaginary, only one mode propagates, the other mode giving rise to an evanescent wave.

Remark 2. At the expense of more intricacy, the present formalism may be generalized to wave propagation in a chiral slab located between z = 0 and z = - d. Then, two more fields exist respectively reflected at z = -d inside the slab and refracted outside in the z < -d region, supplying four supplementary amplitudes matched by the boundary conditions at z = -d. But, instead of a 4 ´ 4 system of equations to get the amplitudes of the electromagnetic field, we have to deal with a 8 ´ 8 system  more difficult to solve.

3. Harmonic Plane Wave Propagation in a Two Dimensional Nano-Periodic Medium

With B = µH, D = e (x)E, and exp(-iwt) implicit, the Maxwell equations are for E(x,z), H(x,z)

θ_zE_y - iwm/c H_x = 0, θ_zH_y + iw e(x)/c E_x = 0

θ_zE_x – θ_xE_z + iwm/c H_y = 0θ_zH_x – θ_xH_z – iw e(x)/c E_y = 0

θ_xE_y + iwm/c H_z = 0, θ_xH_y – iwe(x)/c E_z = 0 (47)

with the divergence equations

[e’ + eθ_x]E_x + eθ_zE_z(x,z) = 0, θ_xH_x + θ_zH_z = 0(48)

giving rise to TE (E_y, H_x, H_z) and TM (H_y, E_x, E_z) waves.

3.1 TE Wave Propagation

Assuming f << 1, we work with the Maxwell-Garnett 0(f²) approximation of (7)

e (x) = e₁  +  h f cos(2ax), h = 3ae₁ (49)

The component E_y satisfies the Helmholtz equation in which ∆ = θ_x² + θ_z²

[∆ + w²me (x)/c²]E_y(x,z) = 0(50)

We look for the solutions of this equation in the form, A being an arbitrary amplitude E_y(x,z) = A exp(ik_zz) y (x) (51)

Substituting (51) into (50) and taking into account (49), gives the differential equation satisfied by y(x)

[θ_x² + k₀² + f k_e² cos(2ax)] y (x) = 0 (52)

in which k₀² = w²me₁/c² - k_z², k_e² = w²mh /c²(53)

Using the variable z = k₁ x , Equation (52) becomes a Mathieu equation [14,15]

[θ_z² + c² + f cos(2az/k_e)]y (z) = 0, c² = k₀²/k_e²(54)

with solutions in the form [14,15,16] where v has to be determined y (z) = ∑_m=-¥^¥ c_m exp([i(v + 2m) az/k_e)](55)

Substituting (55) into (54) gives the following recurrence relation [15] for the coefficients c_m

c_m + g_m(v) (c_m-1 + c_m+1) = 0 (56)

with g_m(v) = -f /2 [(2m + v )² - c²](57)

Now, the main difficulty [14,15] is to get v in terms of f and c, but f being small, the infinite determinant of the system (56) supplies v to the 0(f ³) order [15]

cos(vπ) = cos(cπ) + πf ² [4c ²(1 - c ²)^1/2] ^-1 sin(cπ) (58)

Once v known, the c_m coefficients may be obtained by numerical methods based on the recurrence relations (36) or on some variant of it. It is shown [15] how for moderate values of c and f, these relations can be transformed into convergent continued fractions R_m(n) = c_m/c_m-1, L_m(n) = c_m/c_m+1.

So, according to (51) and (55), E_y(x,z) = E_y(x + π/a,z) and E_y(x,z) = A exp(ik_zz) ∑_m=-¥^¥ c_m exp[i(v + 2m)ax)]0 ≤ x < π/a(59)

and taking into account the Maxwell Equation (47), the other two components H_x, H_z of the TE field are obtained from θ_zE_y and θ_xE_y respectively. Writing (59)

E_y(x,z) = A∑_m=-¥^¥ c_m exp(ik_zz + ik_mx)k_m = (v + 2m)a, 0 ≤ x < π/a (60)

E_y(x,z) appears as a periodic beam of plane waves propagating in the directions defined by the wave vectors with components (k_z, k_m), their amplitude being weighted by the coefficients c_m.

3.2 TM Wave Propagation

For TM waves (H_y, E_x, E_z), we start with the expression (7) of e (x). Then, according to the Maxwell Equation (47) the component H_y satisfied the equation

[∆ + w²me (x)/c² - {e’(x)/e(x)}θ_x] H_y(x,z) = 0 (61)

We look for the solutions of (61) in the form H_y(x,z) = A exp(ik_zz) y(x) (62)

y (x) = u(x) f (x), f (x) = exp[ f₁/2 cos(2ax)]f₁ = f h (63)

A simple calculation gives the first and second derivative of y (x)

y’(x) = [u’/u - af₁ sin(2ax)]y (x)

 y’’(x) = [u’’/u - 2a u’/u f₁ sin(2ax)

- 2a² f₁ cos(2ax) + a²f₁² sin² (2ax)]y (x)(64)

and since e’/e = -2a f₁ sin(2ax), we get to the 0( f₁²) order y’’ - e’/e y’ = [u’/u - 2a² f₁ cos(2ax)] y (x) + 0( f₁²) (65)

so that

[y’’- e’/e y’]H_y (x,z) = [u’/u - 2a² f₁ cos(2ax)]H_y(x,z)(66)

Then, according to (62) and (66), we get from (61), the differential equation satisfied by u(x)

[θ_x² + w²me(x)/c² - k_z² - 2a² f₁ cos(2ax)]u(x) = 0 (67)

which becomes with the Maxwell-Garnett approximation (49) of e (x)

[θ_x² + k₀² + f k_h² cos(2ax)]u(x) = 0(68)

with k₀² given by (53) while k_h² = w²me₁/c² - 2a²h(69)

The comparison of (52) and (68) shows that, to the 0(f ²) order, one has just to change k_e into k_h to go from TE to TM waves so that all the calculations of Subsection 3.1 can be repeated mutatis mutandis.

3.3 TE Wave Scattering in a Semi Infinite Nano-Periodic Material

The granular material, made of nano dots immersed in a dielectric, lies in the z < 0 half-space and we suppose that a TE harmonic plane wave (E_yⁱ, H_xⁱ, H_zⁱ) impinges from the upper half-space z > 0 with refractive index n and permeability m on the z = 0 interface.

The components E_yⁱ, E_y^r of the incident and reflected waves are E_yⁱ(x,z) = A_i exp[iwn/c (x sinq_i + z cosq_i)]

E_y^r(x,z) = A_r exp[iwn/c (x sinq_i - z cosq_i)](70)

and according to the Maxwell Equation (47), the components H_xⁱ, H_x^r involved in the boundary conditions are H_xⁱ(x,z) = n /m cosq_i E_yⁱ(x,z)H_x^r(x,z) = -n /m cosq_i E_y^r(x,z) (71)

Now, the refracted periodic field in z < 0 has the form (59)

E_y^t(x,z) = A_t exp(ik_zz) ∑_m=-¥^¥ c_m exp[i(v + 2m) ax]0 ≤ x < π/a(72)

and, still using (47)

H_x^t(x,z) = g E_y^t(x,z), g = ck_z /wm(73)

the boundary conditions impose the continuity on z = 0 of E_y and H_x, that is, according to (70)-(73)

(A_i + A_r) exp(iwn /c x sinq_i)

= A_t ∑_{m = –¥}^¥ c_m exp[i(v + 2m)ax]0 ≤ x < π/a (74)

    n /m cosq_i (A_i - A_r) exp(iwn /c x sinq_i)

= gA_t ∑_m=-¥^¥ c_m exp[i(v + 2m)ax]0 ≤ x < π/a (75)

Let us write (74)

A_i + A_r = A_t W( b p)W(b p) = ∑_m=-¥^¥ c_m exp[i(k_m - k_i) ßπ](76)

in which according to (60) and (70)

k_m = (v + 2m)a, k_i = wn /c sinq_i(77)

with 0 ≤ ß < 1 since 0 ≤ x ≤ π/a.

So, taking into account (60), the granular periodic semi-infinite material behaves as a diffraction grating: the beam of plane waves propagating in the directions defined by the wave vectors with components (k_z, k_m) have their amplitudes modulated by the coefficients c_mexp(-ik_i ßπ). And, acccording to (76), the relations (74) and (75) become A_i + A_r = A_t W(b p), n/m cosq_i (A_i - A_r) = gA_t W(b p)0 ≤ ß < 1(78)

from which we get in terms of the incident amplitude A_i

A_r = -(gm - n cosq_i) (gm + n cosq_i)^-1A_iA_t = 2n cosq_i (gm + n cosq_i)^-1 W^-1(b p) A_i(79)

So, the amplitude A_t is not constant on the interval (0, π/a).

4. Discussion

The relation (6), leads to a consistent formalism but further work is needed to prove or to amend it. In any case, two different modes of harmonic plane waves propagate in these chiral materials. The Post constitutive relations used to characterize such media, allow to get exact analytic expressions for the amplitude of the electromagnetic field in each mode, a note-worthy property due, as noticed in the introduction, to the covariance of Post’s relations under the proper Lorentz group. An excellent review of chiral nano-technology may be found in [17] with a discussion of two topics: nanoscale approaches to chiral technology and, corresponding to the situation considered here, nanotechnology that benefits from chirality. In particular, a section is devoted to chiral carbon nanotechnology and the authors conclude “possible applications of such materials in the field of biomedecine and biotechnology range from prepara-tion of novel antibacterial, cyclotonic and drug delivery agents to catalysis and materials science applications”.

Remark: The analysis of Section 2 may be performed in left-handed chiral materials with negative e, m: just change e, m into -|e|, |-m|.

Granular periodic materials are currently used in mechanical engineering and, with the ob-jective to appraise their properties, theoretical studies have been devoted to acoustic wave propagation in these structures [18]. In electrical engineering, photonic crystals [19,20] are the main illustration of periodic nanomaterials and they take an increasing importance in today technology. But, they are not composite with inclusions immersed in a dielectric structure. For instance, a one-dimensional photonic crystal with a permittivity periodic in the direction of propagation may be described by an expansion in which U is the unit step function e(z) = e₁ ∑_n [U(z - 2na) - U(z - {2n + 1}a)]

+ e₂ ∑_n [U(z - {2n + 1}a) - U(z - {2n + 2}a)] (80)

and, the solutions of Maxwell’s equations are the Bloch functions ∑_m c_k,m exp(ikz+2iπmz/a) to be compared with (59) (and (80) with (49)). Incidently, (80) has a simple expression in terms of the square-sine function e(z) = e + r sin(az) / |sin(az)| U(z)e₁ = e + r, e₂ = e - r(81)

which suggests to work with the Laplace transform of Maxwell’s equations since tanh(πp/2a) is the Laplace transform of the square-sine function [21]. People fluent with the Laplace transform, could think in terms of p instad of z as they use to do with w instead of t.

In opposite to photonic crystals, composite granular materials with a continuous filling factor have no lattice structure and, as shown in Section 3.3, they rather behave as a smooth dielectric grating [22]. Some of the restrictive assumptions on the filling factor f could be somewhat released at the expense of more intricacy:

1) It would be interesting to check what happens when a higher order approximation than 0( f ²) is used;

2) When f(x) = f cos(2ax) is changed into f (x) = ∑₀^¥ f_m(cos2max), the Mathieu equation becomes a Hill equation [14,15] with solutions similar to (55) but the recurrence relations bet-ween the coefficients c_m is more intricate;

3) Finally a generalization to a two-dimensional filling factor f (x,y), periodic in x and y would approach more closely a real physical situation.

To sum up, the application of the Maxwell-Garnett theory to nano composites deserves further research, taking into account the innocuity or not of such materials in biomedecine [23]. This theory is also used to analyze, in the frame of surface plasmon polaritons, the scattering of TE, TM light waves from a composite material made of metallic nano spherical particles immersed inside a metallic structure such as Ag particles in a Sio² matrix [24].

The 0( f ²) Maxwell-Garnett approximation of the periodic permittivity in the nanodoped medium of Section 3 implies that TE, and TM fields are solutions of the Mathieu equation as if they were diffracted from a dielectric grating [25].

REFERENCES