Файл: Doicu A., Wriedt T., Eremin Y.A. Light scattering by systems of particles (OS 124, Springer, 2006.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 28.06.2024
Просмотров: 831
Скачиваний: 0
1.4 Scattered Field |
49 |
be the electromagnetic power received by the detector downstream from the particle, and W0 the electromagnetic power received by the detector if the particle is removed. Evidently, W0 > W and we say that the presence of the particle has resulted in extinction of the incident beam. For a nonabsorbing medium, the electromagnetic power removed from the incident beam W0 −W is accounted for by absorption in the particle and scattering by the particle.
We now consider extinction from a computational point of view. In order to simplify the notations we will use the conventional expressions of the Poynting vectors and the electromagnetic powers in terms of the transformed fields introduced in Sect. 1.1 (we will omit the multiplicative factor 1/√ε0µ0). The time-averaged Poynting vector S can be written as [17]
S = 12 Re {E × H } = Sinc + Sscat + Sext ,
where E = Es + Ee and H = Hs + He are the total electric and magnetic fields,
Sinc = 12 Re {Ee × He }
is the Poynting vector associated with the external excitation,
Sscat = 12 Re {Es × Hs }
is the Poynting vector corresponding to the scattered field and
Sext = 12 Re {Ee × Hs + Es × He }
is the Poynting vector caused by the interaction between the scattered and
incident fields.
Taking into account the boundary conditions n×Ei = n×E and n×Hi = n × H on S, we express the time-averaged power absorbed by the particle as
|
1 |
|
Wabs = − |
|
S n · Re {Ei × Hi } dS |
2 |
||
= − |
1 |
|
|
S n · Re {E × H } dS . |
|
2 |
With Sc being an auxiliary surface enclosing S (Fig. 1.11), we apply the Green second vector theorem to the vector fields E and E in the domain D bounded
by S and Sc. We obtain |
|
|
|
S n · (E × H + E × H) dS = |
Sc n · (E × H + E × H) dS , |
and further |
|
|
|
S n · Re {E × H } dS = |
Sc n · Re {E × H } dS . |
50 1 Basic Theory of Electromagnetic Scattering
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S |
|
|
|
|
|
Es |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
Ee |
|
ke |
|
|
|
|
|
|
|
|
|
|
|
|
|
E = Ee+Es |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
O |
Di |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Sc |
|
|
|
|
|
|
|
|
|
|
Fig. 1.11. Auxiliary surface Sc |
|
|
|
||||||||||||||
The time-averaged power absorbed by the particle then becomes |
|
|||||||||||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Wabs = − |
|
Sc n · Re {E × H } dS = − |
Sc n · S dS |
|
||||||||||||||||||
|
|
2 |
|
|||||||||||||||||||||
where |
|
|
= Winc − Wscat + Wext , |
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
Winc = − |
Sc n · Sinc dS = − |
|
|
|
Sc n · Re {Ee × He } dS , |
(1.80) |
||||||||||||||||||
2 |
|
|||||||||||||||||||||||
W |
|
= |
n |
S |
scat |
dS = |
1 |
|
|
|
|
n |
· |
Re |
{ |
E |
s × |
H |
} |
dS , |
(1.81) |
|||
|
|
|
|
|
||||||||||||||||||||
|
scat |
Sc |
|
· |
|
|
2 |
|
Sc |
|
|
|
s |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|||
Wext = − |
Sc n · Sext dS = − |
|
|
|
Sc n · Re {Ee × Hs + Es × He } dS . |
|||||||||||||||||||
2 |
|
(1.82)
The divergence theorem applied to the excitation field in the domain Dc bounded by Sc gives
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Dc · (Ee × He ) dV = |
|
Sc n · (Ee × He ) dS , |
|||||||||||||
whence, using |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
· |
(E |
e × |
H ) = jk |
µ |
s |
| |
H |
e| |
2 |
− |
ε |
E |
e| |
2 |
, |
|
e |
0 |
|
|
|
s | |
|
|
and
Re { · (Ee × He )} = 0 ,
yield
Winc = 0 .
1.4 Scattered Field |
51 |
Thus, Wext is the sum of the electromagnetic scattering power and the electromagnetic absorption power
Wext = Wscat + Wabs .
For a plane wave incidence, the extinction and scattering cross-sections are given by
Cext = |
|
|
Wext |
|
|
, |
(1.83) |
|
1 |
|
|
|
|
2 |
|||
|
εs |
|Ee0 |
| |
|||||
|
2 |
|
µs |
|
|
|
||
Cscat = |
|
Wscat |
|
|
, |
(1.84) |
||
1 |
|
|
|
|
2 |
|||
|
εs |
|Ee0 |
| |
|||||
|
2 |
|
µs |
|
|
|
the absorption cross-section is
Cabs = Cext − Cscat ≥ 0 ,
while the single-scattering albedo is
ω = Cscat ≤ 1 .
Cext
Essentially, Cscat and Cabs represent the electromagnetic powers removed from the incident wave as a result of scattering and absorption of the incident radiation, while Cext gives the total electromagnetic power removed from the incident wave by the combined e ect of scattering and absorption. The optical cross-sections have the dimension of area and depend on the direction and polarization state of the incident wave as well on the size, optical properties and orientation of the particle. The e ciencies (or e ciency factors) for extinction, scattering and absorption are defined as
Qext = |
Cext |
, |
Qscat = |
Cscat |
, |
Qabs = |
Cabs |
, |
|
G |
G |
G |
|||||||
|
|
|
|
|
|
where G is the particle cross-sectional area projected onto a plane perpendicular to the incident beam. In view of the definition of the normalized di erential scattering cross-section, we set G = πa2c , where ac is the area-equivalent-circle radius. From the point of view of geometrical optics we expect that the extinction e ciency of all particles would be identically equal to unity. In fact, there are many particles which can scatter and absorb more light than is geometrically incident upon them [17].
The scattering cross-section is the integral of the di erential scattering cross-section over the unit sphere. To prove this assertion, we express
Cscat as
|
1 |
|
µs |
|
||
Cscat = |
|
|
|
|
|
Sc er · Re {Es × Hs } dS , |
|Ee0|2 |
|
|
εs |
52 1 Basic Theory of Electromagnetic Scattering
where Sc is a spherical surface situated at infinity and use the far-field representation
|
1 |
|
εs |
|
|
|
|
2 |
1 |
|
|
||||
er · (Es × Hs ) = |
|
|
|
|
|
|
|
|Es∞| |
+ O |
|
|
|
, r → ∞ |
||
r2 |
|
|
µs |
|
r |
||||||||||
to obtain |
1 |
|
|
|
|
|
|
|
|
||||||
Cscat = |
|
|Es∞|2 dΩ . |
(1.85) |
||||||||||||
|
|
|
|
|
2 |
||||||||||
|
|
|
|Ee0| |
|
Ω |
|
|
|
|
|
The scattering cross-section can be expressed in terms of the elements of the phase matrix and the Stokes parameters of the incident wave. Taking into account the expressions of Ie and Is, and using (1.76) we obtain
Cscat = |
1 |
|
Is (er ) dΩ (er ) |
|
|
|
Ie |
|
|||
|
|
|
Ω |
|
|
= |
|
1 |
|
[Z11 (er , ek ) Ie + Z12 (er , ek ) Qe |
|
Ie |
|
||||
|
|
|
Ω |
|
|
|
+ Z13 (er , ek ) Ue + Z14 (er , ek ) Ve] dΩ (er ) . |
(1.86) |
The phase function is related to the di erential scattering cross-section by the relation
p (er , ek ) = |
4π |
|Es∞ (er )|2 |
Cscat |Ee0|2 |
and in view of (1.85) we see that p is dimensionless and normalized, i.e.,
1
4π Ω
p dΩ = 1 .
The mean direction of propagation of the scattered field is defined as
|
g = |
|
1 |
|
2 |
|
|
|
|
|
Ω |Es∞ (er )| |
|
er dΩ (er ) |
(1.87) |
|||
|
Cscat |Ee0|2 |
|
||||||
and obviously |
|
|
|
|
|
|
|
|
g = |
1 |
|
|
|
|
|
|
|
|
|
Ω Is (er ) er dΩ (er ) |
|
|
|
|||
|
CscatIe |
|
|
|
||||
|
1 |
|
|
|
|
|
|
|
= |
|
|
Ω [Z11 (er , ek ) Ie + Z12 (er , ek ) Qe |
|
||||
|
CscatIe |
|
+ Z13 (er , ek ) Ue + Z14 (er , ek ) Ve] er dΩ (er ) .
The asymmetry parameter cos Θ is the dot product between the vector g and the incident direction ek [17, 169],