Файл: Doicu A., Wriedt T., Eremin Y.A. Light scattering by systems of particles (OS 124, Springer, 2006.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 28.06.2024

Просмотров: 763

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

1.1

Maxwell’s Equations and Constitutive Relations

7

where the complex permittivity εt is given by

 

 

 

εt

= ε0εrt = ε0 1 + χe +

jσ

 

 

ωε0

 

 

 

 

with εrt being the complex relative permittivity. Both the conductivity and the susceptibility contribute to the imaginary part of the permittivity, Imt} = ε0Ime} + Re{σ/ω}, and a complex value for εt means that the medium is absorbing. Usually, Ime} is associated with the “bound” charge current density and Re{σ/ω} with the “free” charge current density, and absorption is determined by the sum of these two quantities. Note that for a free-source medium, σ = 0 and εrt = εr = 1 + χe. The simplest solution to Maxwell’s equations in source-free media is the vector plane wave solution. The behavior of a vector plane wave in an isotropic medium is characterized by the dispersion relation

k = ωεµ ,

which relates the wave number k to the properties of the medium and to the angular frequency ω of the wave. The dimensionless quantity

m = cεµ

is the refractive index of the medium, where c = 1/ε0µ0 is the speed of light in vacuum, and if k0 = ωε0µ0 is the wave number in free space, we see that

m =

 

k

.

 

 

 

 

 

 

 

k0

 

The constitutive relations for anisotropic media are

 

D =

 

E ,

 

ε

 

B =

 

H ,

(1.7)

µ

where ε and µ are the permittivity and permeability tensors, respectively. In our analysis we will consider electrically anisotropic media for which the permittivity is a tensor and the permeability is a scalar. Except for amorphous materials and crystals with cubic symmetry, the permittivity is always a tensor, and in general, the permittivity tensor of a crystal is symmetric. Since there exists a coordinate transformation that transforms a symmetric matrix into a diagonal matrix, we can take this coordinate system as reference frame and we have

 

 

 

εx

0

0

 

 

 

 

= 0

εy

0

.

(1.8)

ε

 

 

 

0

0

εz

 

 

 

 

 

 

 

This reference frame is called the principal coordinate system and the three coordinate axes are known as the principal axes of the crystal. If εx = εy = εz ,


81 Basic Theory of Electromagnetic Scattering

the medium is biaxial, and if ε = εx = εy and ε = εz , the medium is uniaxial. Orthorhombic, monoclinic and triclinic crystals are biaxial, while tetragonal, hexagonal and rhombohedral crystals are uniaxial. For uniaxial crystals, the principal axis that exhibits the anisotropy is called the optic axis. The crystal is positive uniaxial if εz > ε and negative uniaxial if εz < ε.

In our analysis, we will investigate the electromagnetic response of isotropic, chiral media exposed to arbitrary external excitations. The lack of geometric symmetry between a particle and its mirror image is referred to as chirality or optical activity. A chiral medium is characterized by either a leftor a right-handedness in its microstructure, and as a result, leftand right-hand circularly polarized fields propagate through it with di ering phase velocities. For a source-free, isotropic, chiral medium, the constitutive relations read as

D = εE + βε × E ,

B = µH + βµ × H ,

where the real number β is known as the chirality parameter. The Maxwell equations can be written compactly in matrix form as

E

E

 

E

 

 

× H

= K H

,

· H

= 0 ,

(1.9)

where

1βk2 jωµ

K = 1 − β2k2 jωε βk2

and k = ωεµ.

Without loss of generality and so as to simplify our notations we make the following transformations:

1 1

E → ε0 E, H → µ0 H ,

D → ε0D, B → µ0B .

As a result, the Maxwell equations for a free-source medium become more “symmetric”:

× E = jk0B ,

× H = jk0D ,

· D = 0 ,

· B = 0 ,

(1.10)

the constitutive relations are given by (1.6) and (1.7) with ε and µ being the relative permittivity and permeability, respectively, the wave number is k = k0εµ, and the K matrix in (1.9) takes the form

1

βk2

jk0µ

(1.11)

K =

 

jk0

ε βk2 .

1 − β2k2



1.2 Incident Field

9

1.2 Incident Field

In this section, we characterize the polarization state of vector plane waves and derive vector spherical wave expansions for the incident field. The first topic is relevant in the analysis of the scattered field, while the second one plays an important role in the derivation of the transition matrix.

1.2.1 Polarization

In addition to intensity and frequency, a monochromatic (time harmonic) electromagnetic wave is characterized by its state of polarization. This concept is useful when we discuss the polarization of the scattered field since the polarization state of a beam is changed on interaction with a particle.

We consider a right-handed Cartesian coordinate system OXY Z with a fixed spatial orientation. This reference frame will be referred to as the global coordinate system or the laboratory coordinate system. The direction of propagation of the vector plane wave is specified by the unit vector ek , or equivalently, by the zenith and azimuth angles β and α, respectively (Fig. 1.2). The polarization state of the incident wave will be described in terms of the vertical polarization unit vector eα = ez ×ek /|ez ×ek | and the horizontal polarization vector eβ = eα × ek . Note that other names for vertical polarization are TM polarization, parallel polarization and p polarization, while other names for horizontal polarization are TE polarization, perpendicular polarization and s polarization.

z

ek

 

 

eα

 

 

epol

 

β

αpol

 

 

 

 

eβ

 

 

ke

O

 

y

 

 

 

α

 

x

Fig. 1.2. Wave vector in the global coordinate system

10 1 Basic Theory of Electromagnetic Scattering

In the frequency domain, a vector plane wave propagating in a medium with constant wave number ks = k0εsµs is given by

Ee(r) = Ee0ejke·r , Ee0 · ek = 0 ,

(1.12)

where k0 is the wave number in free space, ke is the wave vector, ke = ksek , Ee0 is the complex amplitude vector,

Ee0 = Ee0eβ + Ee0eα ,

and Ee0and Ee0are the complex amplitudes in the β- and α-direction, respectively. An equivalent representation for Ee0 is

Ee0 = |Ee0| epol ,

(1.13)

where epol is the complex polarization unit vector, |epol| = 1, and

1

epol = |Ee0| (Ee0eβ + Ee0eα) .

Inserting (1.13) into (1.12), gives the representation

Ee(r) = |Ee0| epolejke·r , epol · ek = 0 ,

and obviously, |Ee(r)| = |Ee0|.

There are three ways of describing the polarization state of vector plane waves.

1. Setting

Ee0= aβ ejδβ ,

 

Ee0= aαejδα ,

(1.14)

where aβ and aα are the real non-negative amplitudes, and δβ and δα are the real phases, we characterize the polarization state of a vector plane wave by aβ , aα and the phase di erence ∆δ = δβ − δα.

2. Taking into account the representation of a vector plane wave in the time domain

Ee(r, t) = Re Ee(r)ejωt = Re Ee0ej(ke·rωt) ,

where Ee(r, t) is the real electric vector, we deduce that (cf. (1.14))

Ee(r, t) = aβ cos (δβ + ke · r − ωt) ,

Ee(r, t) = aα cos (δα + ke · r − ωt) ,

where

Ee(r, t) = Ee(r, t)eβ + Ee(r, t)eα .


1.2 Incident Field

11

At any fixed point in space the endpoint of the real electric vector describes an ellipse which is also known as the vibration ellipse [17]. The vibration ellipse can be traced out in two opposite senses: clockwise and anticlockwise. If the real electric vector rotates clockwise, as viewed by an observer looking in the direction of propagation, the polarization of the ellipse is right-handed and the polarization is left-handed if the electric vector rotates anticlockwise. The two opposite senses of rotation lead to a classification of

vibration ellipses according to their handedness. In addition to its handed-

ness, a vibration ellipse is characterized by E0 = a2 + b2, where a and b are the semi-major and semi-minor axes of the ellipse, the orientation angle ψ and the ellipticity angle χ (Fig. 1.3). The orientation angle ψ is the angle between the α-axis and the major axis, and ψ [0, π). The ellipticity angle χ is usually expressed as tan χ = ±b/a, where the plus sign corresponds to right-handed elliptical polarization, and χ [−π/4, π/4].

We now proceed to relate the complex amplitudes Ee0and Ee0to the ellipsometric parameters E0, ψ and χ. Representing the semi-axes of the

vibration ellipse as

 

b = ±E0 sin χ ,

 

a = E0 cos χ ,

(1.15)

where the plus sign corresponds to right-handed polarization, and taking into account the parametric representation of the ellipse in the principal coordinate system Oα β

Ee(r, t) = ±b sin (ke · r − ωt) ,

Ee(r, t) = a cos (ke · r − ωt) ,

β

β

 

 

 

 

 

 

 

 

α

 

 

χ

 

b

a

ψ

α

 

 

 

 

 

O

 

 

Fig. 1.3. Vibration ellipse