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

Contents XIII

3.8 Complex Particles . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 242 3.9 Particle on or Near a Plane Surface . . . . . . . . . . . . . . . . . . . . . . . . 245 3.10 E ective Medium Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 246

A Spherical Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 253

A.1 Spherical Bessel Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 254

A.2 Legendre Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 256

B Wave Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 261

B.1 Scalar Wave Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 261

B.2 Vector Wave Functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 265

B.3 Rotations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 270

B.4 Translations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 276

References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 303

Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 317

1

Basic Theory of Electromagnetic Scattering

This chapter is devoted to present the fundamentals of the electromagnetic scattering theory which are relevant in the analysis of the null-field method. We begin with a brief discussion on the physical background of Maxwell’s equations and establish vector spherical wave expansions for the incident field. We then derive new systems of vector functions for internal field approximations by analyzing wave propagation in isotropic, anisotropic and chiral media, and present the T -matrix formulation for electromagnetic scattering. We decided to leave out some technical details in the presentation. Therefore, the integral and orthogonality relations, the addition theorems and the basic properties of the scalar and vector spherical wave functions are reviewed in Appendices A and B.

1.1 Maxwell’s Equations and Constitutive Relations

In this section, we formulate the Maxwell equations that govern the behavior of the electromagnetic fields. We present the fundamental laws of electromagnetism, derive the boundary conditions and describe the properties of isotropic, anisotropic and chiral media by constitutive relations. Our presentation follows the treatment of Kong [122] and Mishchenko et al. [169]. Other excellent textbooks on classical electrodynamics and optics have been given by Stratton [215], Tsang et al. [228], Jackson [110], van de Hulst [105], Kerker [115], Bohren and Hu man [17], and Born and Wolf [19].

The behavior of the macroscopic field at interior points in material media is governed by Maxwell’s equations:

21 Basic Theory of Electromagnetic Scattering

where t is time, E the electric field, H the magnetic field, B the magnetic induction, D the electric displacement and ρ and J the electric charge density and current density, respectively. The first three equations in Maxwell’s theory are independent, because the Gauss magnetic field law can be obtained from Faraday’s law by taking the divergence and by setting the integration constant with respect to time equal to zero. Analogously, taking the divergence of Maxwell–Ampere law and using the Gauss electric field law we obtain the continuity equation:

which expresses the conservation of electric charge. The Gauss magnetic field law and the continuity equation should be treated as auxiliary or dependent equations in the entire system of equations (1.1)–(1.5). The charge and current densities are associated with the so-called “free” charges, and for a sourcefree medium, J = 0 and ρ = 0. In this case, the Gauss electric field law can be obtained from Maxwell–Ampere law and only the first two equations in Maxwell’s theory are independent.

In our analysis we will assume that all fields and sources are time harmonic.

√

With ω being the angular frequency and j = −1, we write

E(r, t) = Re E(r)e−jωt

and similarly for other field quantities. The vector field E(r) in the frequency domain is a complex quantity, while E(r, t) in the time domain is real. As a result of the Fourier component Ansatz, the Maxwell equations in the frequency domain become

× E = j ωB ,

× H = J − j ωD ,

· D = ρ ,

· B = 0 .

Taking into account the continuity equation in the frequency domain · J −j ωρ = 0, we may express the Maxwell–Ampere law and the Gauss electric field law as

× H = −j ωDt ,

· Dt = 0 ,

where

j

Dt = D + ω J

is the total electric displacement.

41 Basic Theory of Electromagnetic Scattering

density ρs such that ρs = limh→0 hρ, we see that the tangential component of E is continuous:

n1 × (E2 − E1) = 0 ,

the tangential component of H is discontinuous:

n1 × (H2 − H1) = J s ,

the normal component of B is continuous:

n1 · (B2 − B1) = 0 ,

and the normal component of D is discontinuous:

n1 · (D2 − D1) = ρs .

Energy conservation follows from Maxwell’s equations. The vector identity

· (a × b) = b · ( × a) − a · ( × b)

yields the Poynting theorem in the time domain:

· (E × H) + H · ∂∂tB + E · ∂∂tD = −E · J ,

and the Poynting vector defined as

S = E × H

is interpreted as the power flow density. Integrating over a finite domain D with boundary S, and using the Gauss theorem, yields

where as before, n is the outward normal unit vector to the surface S. The above equation states that the power supplied by the sources within a volume is equal to the sum of the increase in electromagnetic energy and the Poynting’s power flowing out through the volume boundary. Poynting’s theorem can also be derived in the frequency domain:

· (E × H ) = jω (B · H − E · D ) − E · J ,

where the asterisk denotes a complex-conjugate value. The complex Poynting vector is defined as S = E × H and the term − D E · J dV is interpreted as the complex power supplied by the source.

In practice, the angular frequency ω is such high that a measuring instrument is not capable of following the rapid oscillations of the power flow but

rather responds to a time average power flow. Considering the time-harmonic vector fields a and b,

a(r, t) = 12 a(r)e−jωt + a (r)ejωt , b(r, t) = 12 b(r)e−jωt + b (r)ejωt ,

we express the dot product of the vectors as

c(r, t) = a(r, t) · b(r, t)

= 12 Re a(r) · b (r) + a(r) · b(r)e−2jωt .

where T is a time interval, we derive

c(r) = 12 Re {a(r) · b (r)} ,

while for the cross product of the vectors

c(r, t) = a(r, t) × b(r, t) ,

we similarly obtain

c(r) = 12 Re {a(r) × b (r)} .

Thus, the time average of the dot or cross product of two time-harmonic complex quantities is equal to half of the real part of the respective product of one quantity and the complex conjugate of the other. In this regard, the time-averaged Poynting vector is given by

S = 12 Re {E × H } .

The three independent vector equations (1.1)–(1.3) are equivalent to seven scalar di erential equations, while the number of unknown scalar functions is 16. Obviously, the three independent equations are not su cient to form a complete systems of equations to solve for the unknown functions, and for this reason, the equations given by (1.1)–(1.4) are known as the indefinite form of the Maxwell equations. Note that for a free-source medium, we have six scalar di erential equations with 12 unknown scalar functions. To make the Maxwell equations definite we need more information and this additional

61 Basic Theory of Electromagnetic Scattering

information is given by the constitutive relations. The constitutive relations provide a description of media and give functional dependence among vector fields. For isotropic media, the constitutive relations read as

D = εE ,

B = µH ,

where ε is the electric permittivity, µ is the magnetic permeability and σ is the electric conductivity. The above equations provide nine scalar relations that make the number of unknowns and the number of equations compatible, while for a source-free medium, the first two constitutive relations guarantee this compatibility. When the constitutive relations between the vector fields are specified, Maxwell equations become definite. In free space ε0 = 8.85 × 10−12 F m−1 and µ0 = 4π × 10−7 H m−1, while in a material medium, the permittivity and permeability are determined by the electrical and magnetic properties of the medium. A dielectric material can be characterized by a free-space part and a part depending on the polarization vector P such that

D = ε0E + P .

The polarization P symbolizes the average electric dipole moment per unit volume and is given by

P = ε0χeE ,

where χe is the electric susceptibility. A magnetic material can also be characterized by a free-space part and a part depending on the magnetization vector M ,

B = µ0H+µ0M ,

where M symbolizes the average magnetic dipole moment per unit volume,

M = χmH ,

and χm is the magnetic susceptibility. A medium is diamagnetic if µ < µ0 and paramagnetic if µ > µ0, while for a nonmagnetic medium we have µ = µ0. The permittivity and permeability of isotropic media can be written as

ε = ε0εr = ε0 (1 + χe) ,

µ = µ0µr = µ0 (1 + χm) ,

where εr and µr stand for the corresponding relative quantities. The constitutive relation for the total electric displacement is

Dt = εtE ,