Файл: Doicu A., Wriedt T., Eremin Y.A. Light scattering by systems of particles (OS 124, Springer, 2006.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 813
Скачиваний: 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, Im{εt} = ε0Im{χe} + Re{σ/ω}, and a complex value for εt means that the medium is absorbing. Usually, Im{χe} 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 = Ee0,β eβ + Ee0,αeα ,
and Ee0,β and Ee0,α are 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| (Ee0,β eβ + Ee0,αeα) .
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)e−jω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 Ee0,β and Ee0,α to 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