Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.astro.spbu.ru/DOP/1-DEFS/intro3e.ps Äàòà èçìåíåíèÿ: Fri Nov 19 12:06:51 2010 Äàòà èíäåêñèðîâàíèÿ: Mon Oct 1 21:26:36 2012 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: annular solar eclipse

On optical characteristics of non­spherical scatterers
Abstract. In these brief notes we consider expressions defining some basic optical characteristics of non­spherical
scatterers and point out the differences between these expressions and those for a more simple case of spherical particles.
1. Introduction
Electromagnetic radiation outside and inside any macroscopic body can be described by the vectors ~
E; ~
H
being the strengths of the electric and magnetic fields, respectively.
Usually one assumes that the radiation from a very distant source is incident on a scatterer and then this
radiation can be represented by a plane wave
~
E in = (E in
`
~ i ` +E in
'
~ i ' ) e ik~n in ~r ; (1)
where k = 2ú=– is the wavenumber, – the radiation wavelength in the medium, the vector ~n in shows the direction
of propagation of the wave, i is the complex unity, ~r the radius­vector. Here we introduce the spherical coordinate
system (r; `; '), whose origin coincides with the origin of a Cartesian coordinate system that will be used below,
~ i ` ; ~ i ' are the unit vectors.
Any plane wave can be represented by a superposition of two totally polarized waves -- the TE and TM modes.
For the former, the vector of the electric field is parallel to a selected plane, for the latter, it is perpendicular.
Besides that, one can consider only time­harmonic fields since any field can be expended in Fourier series
whose components can be treated independently (see [1] for more details).
2. Scattering and Maxwell's equations
Scattering of electromagnetic radiation by an isolated macroscopic body is fully described by Maxwell's
equations with the corresponding boundary conditions. In the general case
r \Delta ~
D = ae; r \Theta ~
E + @ ~
B
@t = 0; r \Delta ~
B = 0; r \Theta ~
H = ~
J + @ ~
D
@t (2)
~
D = '' 0 ~
E + ~
P ; ~
H = 1
¯ 0
~
B \Gamma ~
M ; (3)
where D;B; ae; J; M; t are the electric and magnetic inductions, the density of free changes, their current density,
magnetization, and time, '' 0 ; ¯ 0 the dielectric permittivity and magnetic permeability of the medium surrounding
the scatterer (for vacuum '' 0 = ¯ 0 = 1). One can also introduce the material equations
~
J = oe ~
E; ~
B = ¯ ~
H; ~
P = '' 0 ü ~
E; (4)
where oe; ü are the conductivity and electric susceptibility.
The boundary conditions usually look as follows:
( ~
E sca + ~
E in ) \Theta ~n = ~
E int \Theta ~n; ( ~
H sca + ~
H in ) \Theta ~n = ~
H int \Theta ~n; (5)
where all points on the scatterer surface S are considered, ~n is the external normal to S, ~
E in ; ~
H in mean the
incident radiation, ~
E int ; ~
H int the field inside the scatterer, ~
E sca ; ~
H sca the additional (scattered) field appearing
outside the scatterer.
For time­harmonic fields, from Eqs. (2)--(4) one can easily obtain the wave equation. For a spherically
symmetric scatterer and a proper choice of the spherical coordinate system one can separate the variables both
in the equation and the boundary conditions and obtain the solution in the explicit form (so called Mie theory --
see, e.g., [1]). For a non­spherical scatterer of finite size such a separation is impossible in any coordinate system
and solution of the scattering problem becomes more complicated (see, e.g., [2]). Here in contrast with the case
of spheres, the proper choice of the directions of the axes of the Cartesian coordinate systems (there are usually
two required systems -- the laboratory one and that connected with the scatterer) plays in particular important
role.
1

3. Amplitude matrix
If one considers the scattered field only in the far zone (kr AE 1), it can be represented by a spherical wave
(see [1] for more details)
~
E sca =
i
E sca
`;0
~ i ` +E sca
';0
~ i '
j e ikr
r ; (6)
where E sca
`;0 ; E sca
';0 depend only on the direction of propagation of the scattered radiation.
For a scatterer of any shape, the linearity of Maxwell's equations leads to a linear dependence of the fields:
the incident and internal ones, and the incident and scattered ones. In the second case, the relationship can be
written using the amplitude matrix “
S
`
E sca
`
E sca
'
'
= e ikr
r
`
S 11 S 12
S 21 S 22
'`
E in
`
E in
'
'
: (7)
Here to come to the standard form of Eq. (7), we include the known radial dependence of the scattered field in
the components E sca
` = E sca
`;0 e ikr =r and so on.
The elements of the amplitude matrix “
S on the directions of propagation of the incident and scattered
radiations, the size, morphology and chemical composition of the scatterer as well as its orientation in the
laboratory system.
In the coordinate system connected with the particle, the choice of the TE or TM mode (respectively, S 11 ; S 21
or S 12 ; S 22 are utilized) leads to two solutions for the scattered radiation ~
E sca
TE and ~
E sca
TM .
4. Scattering matrix
Another way to describe the electromagnetic fields is to use the Stokes vector
~
I =
0
B B
@
I
Q
U
V
1
C C
A ; (8)
connected with the field strengths as follows:
I = E ` E \Lambda
` +E'E \Lambda
' ; Q = E ` E \Lambda
` \Gamma E'E \Lambda
' ; U = \GammaE ` E \Lambda
' \Gamma E'E \Lambda
` ; V = i(E' E \Lambda
` \Gamma E ` E \Lambda
' ): (9)
The process of scattering described above by the amplitude matrix can be also rewritten using the Stokes
vectors
~
I sca = 1
r 2
“
Z ~ I in ; (10)
where the elements of the scattering matrix “
Z can be explicitly expressed via the elements of the amplitude
matrix “
S (see [1] for more details).
For spheres, 10 out of 16 elements of the scattering matrix are equal to zero, whereas for non­spherical
scatterers it is generally not the case, which can be important in applications (see, e.g., [3]).
5. Cross­sections and albedo
There are several characteristics of scattered radiation considered in the far zone which are suitable and
often used in applications. For example, the extinction cross­section C ext shows what fraction of the incident
flux is absorbed and scattered by a particle. The product of the incident energy flux and this cross­section gives
the power removed from the incident beam as a result of its interaction with the scatterer. This cross­section
can be defined as follows:
C ext = 4ú
k Im
h
~
E sca ; ~ i
i
\Theta=0
; (11)
where \Theta is the scattering angle, i.e. the angle between the directions of the incident and scattered radiation, ~ i
the unit vector parallel to ~
E in .
2

In the case of polarized incident radiation (see, e.g., [4])
C ext = 1
2 (C TM
ext + C TE
ext ) + 1
2 (C TM
ext \Gamma C TE
ext ) cos 2\Psi Q in
I in + 1
2 (C TM
ext \Gamma C TE
ext ) sin 2\Psi U in
I in : (12)
Here (I in ; Q in ; U in ; V in ) T is the Stokes vector of the incident radiation ~ I in , \Psi the positional angle of its polar­
ization.
For spheres, obviously
C ext = 1
2 (C TM
ext + C TE
ext ); (13)
and C TM
ext = C TE
ext . For the scatterer of an arbitrary shape but non­polarized incident radiation, the expression
(13) is also true.
The scattering cross­section shows what fraction of radiation is scattered by a particle in all directions (the
product of the incident energy flux and this cross­section is equal to the total power of scattered radiation)
C sca = 1
I in
Z
4ú
I sca
d\Omega ; (14)
where I in ; I sca = j ~
E sca j 2 are the intensities of the incident and scattered radiation, respectively.
From the energy conservation law, it follows that the absorption cross­section is (see [1] for more details)
C abs = C ext \Gamma C sca : (15)
The important parameter is the single scattering albedo \Lambda defined as the fraction of the scattered radiation
in the total radiation excluded from the incident beam
\Lambda = C sca
C ext
; (16)
for non­polarized incident radiation
\Lambda = C TM
sca + C TE
sca
C TM
ext + C TE
ext
: (17)
6. Asymmetry factor and radiation pressure
The intensity of radiation scattered in different directions is proportional to the phase function. For spheres,
the phase function has the rotational symmetry relative to the direction of propagation of the incident radiation.
For non­spherical particles, such a symmetry is generally absent. Therefore, while for spheres the phase function
asymmetry factor g (or hcos \Thetai) is a scalar, for non­spherical particles it is a vector with the components (see
also [5])
g x;y;z = 1
C sca
Z
4ú
~n sca ~ i x;y;z I sca
d\Omega ; (18)
where ~n sca is the direction of propagation of the scattered radiation, ~ i x;y;z are the unit vectors of a Cartesian
system.
The absence of symmetry in the spatial distribution of the scattered radiation leads to the fact that the total
recoil of scattered photons is not directed along the direction of propagation of the incident radiation. Therefore,
in contrast to the spheres, for non­spherical particles the cross­section of radiation pressure is a vector too
~
C pr = C ext
~
R
R \Gamma C sca ~g; (19)
where R is the distance between the scatterer and the radiation source. For spheres,
C pr = C ext \Gamma gC sca = C abs + (1 \Gamma g)C sca : (20)
3

7. Absorption and re­emission
The absorption process can be described as follows:
~ I after = “
R \Gamma1 “
\Lambda part “
R ~ I before ; (21)
where ~
I before ( ~ I after ) is the Stokes vector of radiation before (after) absorption.
Here the matrix “
R
“
R =
0
B
B
@
1 0 0 0
0 cos 2\Psi sin 2\Psi 0
0 \Gamma sin 2\Psi cos 2\Psi 0
0 0 0 1
1
C
C
A (22)
makes the transition from the laboratory system to the system connected with the scatterer, and “
\Lambda part is the
albedo matrix in the particle system
“
\Lambda part =
0
B
B
@
~ l 1 ~ l 2 0 0
~ l 2 ~ l 1 0 0
0 0 ~ l 1 0
0 0 0 ~ l 1
1
C
C
A : (23)
Here ~ l 1 and ~ l 2 are
~ l 1 = \Lambda TM
1 + C TM
ext =C TE
ext
+ \Lambda TE
1 + C TE
ext =C TM
ext
;
~ l 2 = \Lambda TM
1 + C TM
ext =C TE
ext
\Gamma \Lambda TE
1 + C TE
ext =C TM
ext
; (24)
where \Lambda TE;TM = C TE;TM
sca =C TE;TM
ext is the albedo of totally polarized radiation. In the case of spherical particles
instead of “
\Lambda one has multiplication by a number -- the value of albedo.
In the case of re­emission of radiation, the Stokes vector can be written in the particle system as follows:
~ I emis /
0
B
B
@
(C TM
abs + C TE
abs )=2
\Gamma(C TM
abs \Gamma C TE
abs )=2
0
0
1
C
C
A B – (T d ); (25)
where B – (T d ) is the Plank function, T d the temperature of dust grains.
References
[1] Bohren C.F., Huffman D.R. (1983) Absorption and Scattering of Light by Small Particles. J.Wiley &
Sons, New York.
[2] Mishchenko M.I., Hovenier J.W., Travis L.D. (2000) Light Scattering by Nonspherical Particles. Academic
Press, San Diego.
[3] Gledhill T., McCall A. (2000) Mon. Not. Roy. Astr. Soc. 314, 123.
[4] Martin P.G. (1974) Astrophys.J. 187, 461.
[5] Draine B.T., Flatau (1997) Princeton Obs. Preprint POPe­695.
4