## Abstract

The study of light propagating and scattering for various particles has always been important in many practical applications, such as optical diagnostics for combustion, monitoring of atmospheric pollution, analysis of the structure and pathological changes of the biological cell, laser Doppler technology, and so on. This chapter discusses propagation and scattering through particles. The description of the solution methods, numerical results, and potential application of the light scattering by typical particles is introduced. The generalized Lorenz-Mie theory (GLMT) for solving the problem of Gaussian laser beam scattering by typical particles with regular shapes, including spherical particles, spheroidal particles, and cylindrical particles, is described. The numerical methods for the scattering of Gaussian laser beam by complex particles with arbitrarily shape and structure, as well as random discrete particles are introduced. The essential formulations of numerical methods are outlined, and the numerical results for some complex particles are also presented.

### Keywords

- light scattering
- small particles
- Gaussian laser beam
- generalized Loren-Mie theory
- numerical method

## 1. Introduction

The investigation of light propagation and scattering by various complex particles is of great importance in a wide range of scientific fields, and it has lots of practical applications, such as detection of atmospheric pollution, optical diagnostics for aerosols, remote sensing of disasters [1–3]. Over the past few decades, some theories and numerical methods have been developed to study the light wave propagation and scattering through various particles. For the particles with special shape, such as spheres, spheroids, and cylinders, the generalized Lorenz-Mie theory (GLMT) [4–15] can obtain an analytic solution in terms of a limited linear system of equations using the method of separation of variables to solve the Helmholtz equation in the corresponding coordinate system. For the complex particles of arbitrary shapes and structure, some numerical methods, such as the discrete dipole approximation (DDA), the method of moments (MOM), the finite element method (FEM), and the finite-difference time-domain (FDTD), have been utilized. For random media composed of many discrete particles, the T-matrix method, the sparse-matrix canonical-grid (SMCG) method, and the characteristic basis function method (CBFM) can be applied to obtain simulation results.

This chapter discusses the light propagation and scattering through particles. Without loss of generality, the incident light is assumed to be Gaussian laser beam, which can be reduced to conventional plane wave. The detailed description of the solution methods, numerical results, and potential application of the light scattering by systems of particles is introduced.

## 2. Light scattering by regular particles

### 2.1. Light scattering by a homogeneous sphere

The geometry of light scattering of a Gaussian beam by a homogeneous sphere is illustrated in Figure 1. As shown in Figure 1, two Cartesian coordinates

Within the framework of the GLMT, the electromagnetic field components of the illuminating beam are described by partial wave expansions over a set of basic functions, e.g., vector spherical wave functions in spherical coordinates, vector spheroidal wave functions in spheroidal coordinates, and vector cylindrical wave functions in cylindrical coordinates. The expansion coefficients or sub-coefficients are named as beam-shape coefficients (BSCs) are denoted as * X* is TE, transverse electric, or TM, transverse magnetic, with

where the superscript “

The

the index corresponds to the spherical Bessel functions of the first, second, third, or fourth kind (

The electric component of the internal field and the scattered field can be expanded in terms of vector spherical wave functions in the particle coordinate system

To solve the scattering problem, light scattering, the scattering coefficients

where

where

Figure 2 presents the normalized DSCS for the scattering of a Gaussian beam by a homogeneous spherical dielectric particle. The radius of the spherical particle is

### 2.2. Light scattering by a spheroidal particle

Light scattering by a spheroid has been of great interest to many researchers in the past several decades since it provides an appropriate model in many practical situations. For example, during the atomization processes, the shape of fuel droplets departs from sphere to spheroid when it impinges on the wall and breaks. Due to the inertial force, the raindrop also departs from the spherical particle to the near-spheroidal one. A rigorous solution to the scattering problem concerning a homogeneous spheroid illuminated by a plane wave was first derived by Asano and Yamamoto [16]. It was extended later to the cases of shaped beam illumination [17], a layered spheroid [18], and a spheroid with an embedded source [19]. Nevertheless, only parallel incident of the shaped beam, including on-axis and off-axis Gaussian beam scattered by a spheroid was studied, that is to say, the propagation direction of the incident beam is assumed to be parallel to the symmetry axis of the spheroid. An extension of shaped beam scattering with arbitrary incidence was developed within the framework of GLMT by Han et al. [20–22] and Xu et al. [23, 24].

To deal with the shaped beam scattering of a spheroidal particle within the framework of GLMT, the incident Gaussian beam is required to be expanded in terms of the vector spheroidal wave functions in spheroidal coordinates, which can be achieved using the relationship between the vector spherical wave functions and the spherical wave functions. The geometry of shaped beam scattering by a prolate spheroid is illustrated in Figure 3. According to the expansion of shaped beam in unrotated spherical coordinates in Eq. (1), we can rewrite it as

where we have

where

Considering the vector spheroidal wave functions in the spheroidal coordinates, whose explicit expressions are the same as the ones used in Refs. [14, 25], the relationship between the vector spherical wave functions and vector spheroidal wave functions is given as:

From Eq. (15), we can obtain the expansion of Gaussian beam in spheroidal coordinates. Accordingly, the electric component of the internal field and the scattered field can be expanded in terms of vector spheroidal wave functions. The unknown scattered coefficients can be determined by applying the boundary conditions of continuity of the tangential electromagnetic fields over the surface of the particle. Thus, the solution of scattering for Gaussian beam by a spheroidal particle can be obtained.

For the purpose of demonstration, Figure 4 shows angular distributions of the DSCS for a spheroid with a semimajor axis and a semiminor axis being

### 2.3. Light scattering by a circular cylindrical particle

The geometry of shaped beam scattering by a circular cylinder is illustrated in Figure 5. Similarly to a spheroidal particle, due to the lack of spherical symmetry, the arbitrary orientation is also compulsory in the case of GLMTs for cylinders. The expansion of the case of an arbitrary-shaped beam propagating in an arbitrary direction, based on which an approach to expand the shaped beam in terms of cylindrical vector wave functions natural to an infinite cylinder of arbitrary orientation is given below.

The vector cylindrical wave functions in the cylindrical coordinates

where

The relationship between the vector spherical wave functions and the vector cylindrical wave functions is defined as

where

Based on Eqs. (17) and (1), we can obtain the expansion of the incident shaped beam in terms of the cylindrical vector wave functions in cylindrical coordinates as

where

Accordingly, the electric component of the internal field and the scattered field can be expanded in terms of cylindrical vector wave functions. The unknown scattered coefficients can be determined by applying the boundary conditions of continuity of the tangential electromagnetic fields over the surface of the particle. Thus, the solution of scattering for Gaussian beam by a cylindrical particle can be obtained.

## 3. Light scattering by complex particles of arbitrary shapes and structure

### 3.1. Surface integral equation method

Many particles encountered in nature or produced in industrial processes, such as raindrops, ice crystals, biological cells, dust grains, daily cosmetics, and aerosols in the atmosphere, not only have irregular shapes but also have complex structures. The study of light scattering by these complex particles is essential in a wide range of scientific fields with many practical applications, including optical manipulation, particle detection and discrimination, design of new optics devices, etc. Here, we introduce the surface integral equation method (SIEM) [26–28] to simulate the light scattering by arbitrarily shaped particles with multiple internal dielectric inclusions of arbitrary shape, which can be reduced to the case of arbitrarily shaped homogeneous dielectric particles. For SIEM, the incident Gaussian beam can be described using the method of combing Davis-Barton fifth-order approximation [29] in combination with rotation Euler angles [30].

Now, let us consider the problem of Gaussian beam scattering by an arbitrarily shaped particle with multiple dielectric inclusions of arbitrary shape. As illustrated in Figure 6, let

where the integral operators

in which the subscript “0” represents the medium in which the scattered fields are computed and the superscript “

Also based on the surface equivalence principle, the fields

where

By enforcing the continuity of the tangential electromagnetic fields across each surface, the following integral equations may be established

where the subscripts “

The resultant matrix Eq. (33) can also be solved iteratively by employing the multilevel fast multipole algorithm (MLFMA). Once obtained the unknown equivalent electromagnetic currents, the far-zone scattered fields and DSCS can be calculated.

### 3.2. Numerical results

First, we consider the reduced case of arbitrarily shaped homogeneous dielectric particles. To illustrate the validity of the proposed method, the scattering of a focused Gaussian beam by a homogeneous spherical dielectric particle is considered. The radius of the spherical particle is

To illustrate the validity of the proposed method for composite particles with inclusions, we consider the scattering a Gaussian beam by a spheroidal particle with a spherical inclusion at the center, as shown in Figure 8. The semimajor axis and the semiminor axis of the host spheroid are

Finally, the scattering of an obliquely incident Gaussian beam by a cubic particle containing 27 randomly distributed spherical inclusions is considered to illustrate the capabilities of the proposed method. The center of the host cube is located at the origin of the particle system and the side length of the cube is

## 4. Light scattering by random discrete particles

Due to the wide range of possible applications in academic research and industry, the problem of light scattering by random media composed of many discrete particles is a subject of broad interest. Over the past few decades, some theories and numerical methods have been developed to study the light scattering by random discrete particles [32–48]. In this section, we introduce a hybrid finite element-boundary integral-characteristic basis function method (FE-BI-CBFM) to simulate the light scattering by random discrete particles [49]. In this hybrid technique, the finite element method (FEM) is used to obtain the solution of the vector wave equation inside each particle and the boundary integral equation (BIE) is applied on the surfaces of all the particles as a global boundary condition. To reduce computational burdens, the characteristic basis function method (CBFM) is introduced to solve the resultant FE-BI matrix equation. The incident light is assumed to be Gaussian laser beam.

### 4.1. FE-BI-CBFM for random discrete particles

Now, let us consider the scattering of an arbitrarily incident focused Gaussian beam by multiple discrete particles with a random distribution, as depicted in Figure 12. For simplicity, the background region, which is considered to be free space, is denoted as

where

where

Since Eq. (35) is independent of the excitation, we can remove the interior unknowns to derive a matrix equation that only includes the unknowns on

where

For the convenience of description, we write the relation between

where

It is worth to notice that the calculations of Eqs. (37) and (39) in each particle are independent and can be completely parallelized. Furthermore, since the particles are uniform, the coefficient matrices are the same for each particle. This implies that only one particle needs to be dealt with to obtain all the matrices

To formulate the field in region

Enforcing boundary condition on

and a magnetic field integral equation (MFIE)

where the subscript “

where the subscript

The expressions of the elements for matrices

The above equation can be written in a more compact form as

where

### 4.2. Numerical results and discussion

In what follows, some numerical results are presented. First, we consider the scattering of Gaussian beam by 125 randomly distributed conducting spherical particles with a radius of

We then consider the multiple scattering of an obliquely incident Gaussian beam by 512 randomly distributed inhomogeneous spherical particles with which the fractional volume is

Finally, we use the present numerical method to simulate the multiple scattering of Gaussian beam by 1000 randomly distributed homogeneous dielectric spherical particles with which the fractional volume is

## References

- 1.
Gouesbet G and Onofri FRA. Preface: laser-light and interactions with particles (LIP), 2014. Journal of Quantitative Spectroscopy & Radiative Transfer. 2015;162:1–7. - 2.
Mishechenko MI, Travis LD, and Lacis AA. Scattering, absorption, and emission of light by small particles. 2002, Cambridge: Cambridge University Press. - 3.
Doicu A, Wriedt T, and Eremin YA. Light scattering by systems of particles null-field method with discrete sources: theory and programs. 2006, Berlin: Springer-Verlag. - 4.
Gouesbet G and Gréhan G. Generalized Lorenz-Mie theories. 2011, Berlin: Springer. - 5.
Lock JA and Gouesbet G. Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. I. On-axis beams. The Journal of the Optical Society of America A. 1994;11(9):2503–2515. - 6.
Gouesbet G and Lock JA. Rigorous justification of the localized approximation to the beam-shape coefficients in generalized Lorenz-Mie theory. II. Off-axis beams. The Journal of the Optical Society of America A. 1994;11(9):2516–2525. - 7.
Gouesbet G, Maheu B, and Grehan G. Light scattering from a sphere arbitrarily located in a Gaussian beam, using a Bromwich formulation. The Journal of the Optical Society of America A., 1988;5:1427–1443. - 8.
Wang JJ, Gouesbet G, Han YP, and Grehan G. Study of scattering from a sphere with an eccentrically located spherical inclusion by generalized Lorenz-Mie theory: internal and external field distribution. The Journal of the Optical Society of America A. 2011;28:24–39. - 9.
Han L, Han YP, Wang JJ, and Cui ZW. Internal and near-surface electromagnetic fields for a dielectric spheroid illuminated by a zero-order Bessel beam. Journal of the Optical Society of America a-Optics Image Science and Vision. 2014;31(9):1946–1955. - 10.
Onofri F, Grehan G, and Gouesbet G. Electromagnetic scattering from a multilayered sphere located in an arbitrary beam. Applied Optics. 1995;30:7113–7124. - 11.
Gouesbet G, Wang JJ, and Han YP. Transformations of spherical beam shape coefficients in generalized Lorenz-Mie theories through rotations of coordinate systems. I. General formulation. Optics Communication. 2010;283:3218–3225. - 12.
Wang JJ, Gouesbet G, and Han YP. Transformations of spherical beam shape coefficients in generalized Lorenz-Mie theories through rotations of coordinate systems. II. Axisymmetric beams. Optics Communication. 2010;283:3226–3234. - 13.
Gouesbet G, Lock JA, Wang JJ, and Grehan G. Transformations of spherical beam shape coefficients in generalized Lorenz-Mie theories through rotations of coordinate systems. V. Localized beam models. Optics Communication. 2011;284:411–417. - 14.
Han YP, Zhang Y, Zhang HY, and Han GX. Scattering of typical particles by beam shape in oblique illumination. Journal of Quantitative Spectroscopy & Radiative Transfer. 2009;110:1375–1381. - 15.
Edmonds AR. Angular momentum in quantum mechanics. 1957, Princeton: Princeton University Press. - 16.
Asano S and Yamamoto G. Light scattering by a spheroid particle. Applied Optics. 1975;14:29–49. - 17.
Barton JP. Internal and near-surface electromagnetic fields for a spheroidal particle with arbitrary illumination. Applied Optics. 1995;34:5542–5551. - 18.
Barton JP. Internal, near-surface, and scattered electromagnetic fields for a layered spheroid with arbitrary illumination. Applied Optics. 2001;40(21):3598–3607. - 19.
Barton JP. Electromagnetic fields for a spheroidal particle with an arbitrary embedded sources. The Journal of the Optical Society of America A. 2000;17:458–464. - 20.
Yiping H and Zhensen W. The expansion coefficients of a spheroidal particle illuminated by Gaussian beam. IEEE Transactions on Antennas and Propagation. 2001;49(4):615–620. - 21.
Han YP, Grehan G, and Gouesbet G. Generalized Lorenz-Mie theory for a spheroidal particle with off-axis Gaussian-beam illumination. Applied Optics. 2003;42:6621–6629. - 22.
Han YP, Mees L, Ren KF, Grehan G, Wu ZS, and Gouesbet G. Far scattered field from a spheroid under a femtosecond pulsed illumination in a generalized Lorenz-Mie theory framework. Optics Communication. 2004;231:71–77. - 23.
Xu F, Ren KF, Gouesbet G, Grehan, and Cai X. Generalized Lorenz-Mie theory for an arbitrary oriented, located, and shaped beam scattered by homogeneous spheroid. The Journal of the Optical Society of America A. 2007;24:119–131. - 24.
Xu F, Ren K, and Cai X. Expansion of an arbitrarily oriented, located, and shaped beam in spheroidal coordinates. The Journal of the Optical Society of America A. 2007;24(1):109–118. - 25.
Flammer C. Spheroidal wave functions. 1957, California: Stanford U.P. - 26.
Cui ZW, Han YP, and Zhang HY. Scattering of an arbitrarily incident focused Gaussian beam by arbitrarily shaped dielectric particles. The Journal of the Optical Society of America B. 2011;28:2625–2632. - 27.
Han YP, Cui ZW, and Gouesbet G. Numerical simulation of Gaussian beam scattering by complex particles of arbitrary shape and structure. Journal of Quantitative Spectroscopy & Radiative Transfer. 2012;113:1719–1727. - 28.
Han YP, Cui ZW, and Zhao WJ. Scattering of Gaussian beam by arbitrarily shaped particles with multiple internal inclusions. Optics Express. 2012;20:718–731. - 29.
Barton JP and Alexander DR. Fifth-order corrected electromagnetic fields components for a fundamental Gaussian beam. Journal of Applied Physics. 1989;66:2800–2802. - 30.
Edmonds AR. Angular momentum in quantum mechanics. 1957, Princeton: Princeton University Press. - 31.
Tsang L, Kong JA, Ding KH, and Ao CO. Scattering of electromagnetic waves, numerical simulations. 2001, New York: Wiley. - 32.
Ishimaru A. Wave propagation and scattering in random media. 1978, New York: Academic. - 33.
Foldy LL. The multiple scattering of waves. I. General theory of isotropic scattering by randomly distributed scatterers. Physical Review Letters. 1945;67:107–119. - 34.
Lax M. Multiple scattering of waves. Reviews of Modern Physics. 1951;23:287–310. - 35.
Varadan VV and Varadan VK. Multiple scattering of electromagnetic waves by randomly distributed and oriented dielectric scatters. Physical Review. 1980;21:388–394. - 36.
Varadan VK, Bringi VN, Varadan VV, and Ishimaru A. Multiple scattering theory for waves in discrete random media and comparison with experiments. Radio Science. 1983;18:321–327. - 37.
Furutsu K. Multiple scattering of waves in a medium of randomly distributed particles and derivation of the transport equation. Radio Science. 1975;10:29–44. - 38.
Tsang L, Kong JA, and Ding KH. Scattering of electromagnetic waves, theories and applications. 2000, New York : Wiley. - 39.
Tishkovets VP and Jockers K. Multiple scattering of light by densely packed random media of spherical particles: dense media vector radiative transfer equation. Journal of Quantitative Spectroscopy & Radiative Transfer. 2006;101:54–72. - 40.
Lu CC, Chew WC, and Tsang L. The application of recursive aggregate T-matrix algorithm in the Monte Carlo simulations of the extinction rate of random distribution of particles. Radio Science. 1995;30:25–28. - 41.
Chew WC, Lin JH, and Yang XG. An FFT T-matrix method for 3D microwave scattering solution from random discrete scatterers. Microwave and Optical Technology Letters. 1995;9:194–196. - 42.
Siqueira PR and Sarabandi K. T-matrix determination of effective permittivity for three-dimensional dense random media. IEEE Transactions on Antennas and Propagation. 2000;48:317–327. - 43.
Mishchenko MI, Liu L, Mackowski DW, Cairns B, and Videen G. Multiple scattering by random particulate media: exact 3D results. Optics Express. 2007;15:2822–2836. - 44.
Chart CH and Tsang L. A sparse-matrix canonical-grid method for scattering by many scatterers. Microwave and Optical Technology Letters. 1995;8:114–118. - 45.
Barrowes BE, Ao CO, Teixeira FL, and Kong JA. Sparse matrix/canonical grid method applied to 3-D dense medium simulations. IEEE Transactions on Antennas and Propagation. 2003;51:48–58. - 46.
Cui ZW, Han YP, and Xu Q. Numerical simulation of multiple scattering by random discrete particles illuminated by Gaussian beams. The Journal of the Optical Society of America A. 2011;28:2200–2208. - 47.
Sun YF, Chan CH, Mittra R, and Tsang L. Characteristic basis function method for solving large problem arising in dense medium scattering. IEEE Antennas and Propagation Society International Symposium. 2003;2:1068–1071. - 48.
Mackowski DW and Mishchenko MI. Direct simulation of multiple scattering by discrete random media illuminated by Gaussian beams. Physical Review Letters A. 2011;83:013804. - 49.
Cui ZW, Han YP, and Li CY. Simulation of electromagnetic scattering by random discrete particles using a hybrid FE-BI-CBFM technique. Waves Random Complex Media. 2012;22:207–221. - 50.
Jin JM. The finite element method in electromagnetics. 2002, New York: Wiley.