## Abstract

We employ the numerically exact superposition *T*-matrix method to perform extensive computations of electromagnetic scattering by a volume of discrete random medium densely filled with increasingly absorbing as well as non-absorbing particles. Our numerical data demonstrate that increasing absorption diminishes and nearly extinguishes certain optical effects such as depolarization and coherent backscattering and increases the angular width of coherent backscattering patterns. This result corroborates the multiple-scattering origin of such effects and further demonstrates the heuristic value of the concept of multiple scattering even in application to densely packed particulate media.

© 2007 Optical Society of America

## 1. Introduction

Several recent publications [1–5] have marked the emergence of an accurate quantitative approach to the problem of light scattering by macroscopic media composed of randomly positioned particles based on numerically exact solutions of the Maxwell equations. In particular, in [4] we simulated the onset and evolution of various multiple-scattering effects by increasing the number of particles, *N*, in a statistically homogeneous volume of discrete random medium. We argued that if specific scattering features either intensify or weaken with *N* in an expected way [6, 7] then these features can be attributed to the increasing effect of multiple scattering. The particles in [4] were assumed to be non-absorbing by assigning to them a real-valued relative refractive index. However, it is well known that absorption can be an efficient suppressant of multiple scattering [6–8] and reduces the average length of multiple-scattering paths [9]. It is, therefore, important to verify whether the multiple-scattering interpretation of the various optical displays in [4] is consistent with the way in which these features should be expected to change upon increased absorption. This verification is the main objective of the present paper.

## 2. Concept of multiple scattering

Before describing and discussing the results of numerical computations, it is instructive to explain the primordial origin of the concept of multiple scattering. The majority of theoretical techniques based on a direct solution of the frequency-domain differential Maxwell equations or their integral counterparts are applicable to an arbitrary fixed finite object, be it a single physical body or a cluster consisting of several distinct components [10, 11]. These techniques are based on treating the object as a single scatterer and yield the cumulative scattered electric field. However, if the object is a multi-particle group, such as a cloud of water droplets, then it is often convenient to represent the scattered field at a point **r** as a superposition of partial fields scattered by the individual particles, i.e.,

where *N* is the number of particles in the group, **E**(**r**) is the total field, **E**
^{inc}(**r**) is the incident field, and **E**
_{sca}
_{i}(**r**) is the *i*th partial scattered field. The partial scattered fields can be found by solving the vector Foldy–Lax equations (FLEs) [7]. Specifically, the *i*th partial scattered field is given by

where *V _{i}* is the volume occupied by the

*i*th particle, $\overleftrightarrow{G}$ (

**r**,

**r′**) is the free-space dyadic Green’s function, and

**E**

_{i}(

**r**′′) is the electric field “exciting” particle

*i*. The

*N*dyadics ${\overleftrightarrow{T}}_{i}$ can be found by solving the following equation for each

*i*separately:

where *k*
_{1} is the wave number in the host medium, *m _{i}*(

**r**) is the relative refractive index, and $\overleftrightarrow{I}$ is the identity dyadic. The ${\overleftrightarrow{T}}_{i}$ is the dyadic transition operator of particle

*i*with respect to the fixed laboratory coordinate system computed in the absence of all the other particles. In other words, the

*N*dyadic transition operators are totally independent of each other. However, the exciting fields are interdependent and must be found by solving the following system of

*N*linear integral equations:

Let us rewrite Eqs. (1), (2), and (4) in the following compact operator form:

where

Iterating Eq. (6) yields

while the substitution of Eq. (8) in Eq. (5) gives what can be called an order-of-scattering expansion of the total electric field:

Indeed, *Ĝ*
*T̂ _{i}*

*E*

_{inc}can be interpreted as the partial scattered field at the observation point generated by particle

*i*in response to the excitation by the incident field only,

*Ĝ*

*T̂*

_{i}*Ĝ*

*T̂*

_{j}*E*

_{inc}is the partial field generated by the same particle in response to the excitation caused by particle

*j*in response to the excitation by the incident field, etc. This order-of-scattering interpretation of Eq. (9) becomes even more transparent when the particles are widely separated, and single-scattering dyadics replace the dyadic transition operators as complete electromagnetic descriptors of the individual particles [7].

It is important to recognize that besides being an interpretation and visualization tool, the concept of multiple scattering does not represent a physical process per se in the framework of frequency-domain electromagnetics. Indeed, it follows from Eq. (4) that all mutual particle-particle excitations occur simultaneously and are not temporally discrete and ordered events. Nevertheless, Eq. (9) constitutes a very fruitful way of re-writing the original FLEs, while the “multiple scattering” terminology is a convenient and compact way of illustrating and interpreting their solutions [4].

## 3. Numerical results

It is well known from the theory of radiative transfer in homogeneous plane-parallel atmospheres [6–8] that the contribution of light scattered *n* times to the total specific intensity is proportional to the *n*th power of the single-scattering albedo. The single-scattering albedo is equal to unity for Im(*m*) = 0 but can decrease significantly as Im(*m*) increases, thereby suppressing various manifestations of multiple scattering. The radiative transfer theory and the concept of the single-scattering albedo may not be applicable directly to densely packed particles. Still, the various effects identified in [4] as being caused by multiple scattering should be expected to diminish and ultimately vanish with increasing absorption. Indeed, if absorption diminishes the response of a particle to a certain excitation, this effect should have a stronger influence on multiply scattered radiation than on singly scattered ration: the *n*-th order of scattering would suffer *n* times in succession from such a diminished response, Eq. (9).

To verify whether this is indeed the case, we have performed numerically exact computations of far-field scattering by a spherical volume of discrete random medium filled with 160 identical spherical particles (Fig. 1). The size parameters of the volume and the particles are fixed at *k*
_{1}
*R* = 40 and *k*
_{1}
*r* = 4, respectively, which means that the particle volume concentration is equal to 16%. The real part of the particle refractive index is fixed at 1.32, while the imaginary part is varied from 0 to 0.3. To perform averaging over particle positions, we use only one randomly configured 160-particle group and average over all possible orientations of this configuration with respect to the laboratory coordinate system [4]. This approach yields an infinite continuous set of random realizations of the 160-particle group and allows us to employ the highly efficient orientation averaging technique afforded by the superposition *T*-matrix method [10, 12, 13].

We assume that the statistically random particulate volume is illuminated by a plane electromagnetic wave or a parallel quasi-monochromatic beam of light propagating in the direction **n**̂^{inc} (Fig. 1). The observation direction is specified by the unit vector **n**̂^{sca}. Since all scattering properties of the volume are averaged over the uniform orientation distribution of the multi-particle group, we can simplify the discussion by using the scattering plane for defining the Stokes parameters of the incident and scattered light. The transformation of the Stokes parameters in the far-field zone of the volume is then written in terms of the normalized Stokes scattering matrix [4, 6, 7]:

where Θ is the scattering angle (Fig. 1). The phase function is normalized according to

The most relevant numerical results are shown in Fig. 2. For comparison, the yellow curves depict the results computed for a single sphere with *k*
_{1}
*r* = 4 and *m* = 1.32 + i0.3.

Note that increasing Im(*m*) above a certain threshold can cause decreasing rather than increasing absorption since the scatterer eventually behaves as a metallic object. Therefore, we have verified specifically that the single-scattering albedo of the entire scattering volume decreases and the total absorption cross section increases as Im(*m*) is increased from 0 to 0.3.

## 4. Discussion

The ratio *a*
_{2}(Θ)*a*
_{1}(Θ) is identically equal to unity for scattering by a single sphere. Therefore, the large deviation of this ratio from 100% for a scattering volume comprising non-absorbing particles was interpreted in [4] as being a direct consequence of the strengthened depolarizing effect of multiple scattering. The upper left panel of Fig. 2 obviously reinforces this interpretation. Indeed, the deviation of *a*
_{2}(Θ)/*a*
_{1}(Θ) from 100% decreases with increasing Im(*m*) quite significantly, even though it does not completely vanish even for Im(*m*) = 0.3 thereby revealing a residual influence of multiple scattering. The deviation of the ratio *a*
_{3}(180°)*a*
_{1}(180°) from –100% behaves quite similarly. With increasing absorption, the ratios *a*
_{2}(Θ)/*a*
_{1}(Θ), *a*
_{3}(Θ)/*a*
_{1}(Θ), and *b*
_{2}(Θ)/*a*
_{1}(Θ) for the entire scattering volume should be expected to approach those for a single constituent sphere. The left-hand panels of Fig. 2 clearly exhibit this tendency, although the residual multiple-scattering effect diminishes noticeably the large amplitude of the interference features typical of monodisperse spheres. The increasing amplitude of oscillations of the ratio *b*
_{2}(Θ)/*a*
_{1}(Θ) with increasing Im(*m*) and the similar behavior of the ratio –*b*
_{1}(Θ)*a*
_{1}(Θ) (not shown) obviously represent what is known in astrophysics of planetary surfaces as the Umov effect.

An important, inherently multiple-scattering effect is coherent backscattering (CB) [4, 7, 14–16]. The grey curves in the right-hand panels of Fig. 2 demonstrate the manifestations of CB in the same-helicity scattered intensity

as well as in the linear and circular polarization ratios defined, respectively, by

and

It is clearly seen that increasing absorption diminishes and then essentially extinguishes these features, which corroborates their multiple-scattering origin.

Further evidence comes from the examination of the dependence of the angular width of CB patterns on absorption. It is well known that absorption terminates long scattering paths and thereby can be expected to reduce the range of angles affected by CB [9]. This is clearly seen indeed, especially in the panels depicting *μ*
_{L} and *μ*
_{C}.

## 5. Conclusion

The concept of multiple scattering is native to the theories of radiative transfer and CB describing electromagnetic energy transport in low-density particulate media. In fact, the very derivation of the radiative transfer equation from the Maxwell equations is explicitly based on the order-of-scattering expansion of the far-field FLEs [7]. One might argue that the relevance of this concept to densely packed particulate media is not obvious and may not be easy to justify. Our previous results [3, 4] have demonstrated that the evolution of certain scattering patterns with increasing the number of particles in a volume of discrete random medium is qualitatively consistent with the expected result of intensifying multiple scattering. Our new numerical data show that increasing absorption diminishes and nearly extinguishes these scattering patterns and increases the angular width of the CB patterns, thereby further corroborating their multiple-scattering origin. This result provides an additional illustration of the heuristic value of the concept of multiple scattering, even when this concept is applied to densely packed particulate media.

## Acknowledgments

The authors thank Daniel Mackowski, Michiel Min, and Gorden Videen for numerous useful discussions. This research was supported by the NASA Radiation Sciences Program managed by Hal Maring and by the NASA Glory Mission project.

## References and links

**1. **S. H. Tseng, A. Taflove, D. Maitland, and V. Backman, “Pseudospectral time domain simulations of multiple light scattering in three-dimensional macroscopic random media,” Radio Sci. **41**, RS4009 (2006). [CrossRef]

**2. **D. W. Mackowski, “Direct simulation of scattering and absorption by particle deposits,” Proc. IMECE **2006**, 14615 (2006).

**3. **M. I. Mishchenko and L. Liu, “Weak localization of electromagnetic waves by densely packed many-particle groups: exact 3D results,” J. Quant. Spectrosc. Radiat. Transfer **106**, 616–621 (2007). [CrossRef]

**4. **M. I. Mishchenko, L. Liu, D. W. Mackowski, B. Cairns, and G. Videen, “Multiple scattering by random particulate media: exact 3D results,” Opt. Express **15**, 2822–2836 (2007). [CrossRef] [PubMed]

**5. **A. Penttilä and K. Lumme, “Coherent backscattering effects with Discrete Dipole Approximation method,” in *Peer-Reviewed Abstracts of the Tenth Conference on Electromagnetic & Light Scattering*,
G. Videen, M. Mishchenko, M. P. Mengüç, and N. Zakharova, eds. (http://www.giss.nasa.gov/~crmim/, 2007), pp. 157–160.

**6. **J. W. Hovenier, C. van der Mee, and H. Domke, Transfer of Polarized Light in Planetary Atmospheres – Basic Concepts and Practical Methods (Springer, Berlin, 2004).

**7. **M. I. Mishchenko, L. D. Travis, and A. A. Lacis, *Multiple Scattering of Light by Particles: Radiative Transfer and Coherent Backscattering* (Cambridge U. Press, Cambridge, UK, 2006).

**8. **H. C. van de Hulst, Multiple Light Scattering. Tables, Formulas, and Applications (Academic Press, San Diego, 1980).

**9. **S. Etemad, R. Thompson, M. J. Andrejco, S. John, and F. C. MacKintosh, “Weak localization of photons: termination of coherent random walks by absorption and confined geometry,” Phys. Rev. Lett. **59**, 1420–1423 (1987). [CrossRef] [PubMed]

**10. **M. I. Mishchenko, L. D. Travis, and A. A. Lacis, *Scattering, Absorption, and Emission of Light by Small Particles* (Cambridge U. Press, Cambridge, UK, 2002). http://www.giss.nasa.gov/~crmim/books.html.

**11. **F. M. Kahnert, “Numerical methods in electromagnetic scattering theory,” J. Quant. Spectrosc. Radiat. Transfer **79–80**, 775–824 (2003). [CrossRef]

**12. **D. W. Mackowski and M. I. Mishchenko, “Calculation of the *T* matrix and the scattering matrix for ensembles of spheres,” J. Opt. Soc. Am. A **13**, 2266–2278 (1996). [CrossRef]

**13. **D. Mackowski, K. Fuller, and M. Mishchenko, “Codes for calculation of scattering by clusters of spheres.” ftp://ftp.eng.auburn.edu/pub/dmckwski/scatcodes/index.html.

**14. **Y. Kuga and A. Ishimaru, “Retroreflectance from a dense distribution of spherical particles,” J. Opt. Soc. Am. A **1**, 831–835 (1984). [CrossRef]

**15. **L. Tsang and A. Ishimaru, “Backscattering enhancement of random discrete scatterers,” J. Opt. Soc. Am. A **1**, 836–839 (1984). [CrossRef]

**16. **V. L. Kuz’min and V. P. Romanov, “Coherent phenomena in light scattering from disordered systems,” Phys.-Uspekhi **39**, 231–260 (1996). [CrossRef]