Next Article in Journal
Emission Enhancement in fs + ns Dual-Pulse LIBS of Cu
Previous Article in Journal
The Design of Large Curved Waveguide Based on Sunflower Graded Photonic Crystal
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extinction and Independent Scattering Criterion for Clusters of Spherical Particles Embedded in Absorbing Host Media

1
School of Energy and Power Engineering, Shandong University, Jinan 250061, China
2
Optics & Thermal Radiation Research Center, Shandong University, Qingdao 266237, China
*
Author to whom correspondence should be addressed.
Submission received: 5 June 2023 / Revised: 1 July 2023 / Accepted: 3 July 2023 / Published: 5 July 2023

Abstract

:
In practical applications, the independent scattering approximation (ISA) is widely used to analyze light transfer in nanoparticle systems. However, the traditional independent scattering criterion is obtained under the assumption that the host medium surrounding particles is nonabsorbing, and thus may be invalid in certain circumstances. In this work, to explore the applicability of the ISA for small particles in absorbing host media, we calculate the extinction efficiency of particle clusters by direct solutions of macroscopic Maxwell equations. Using the far-field and distance-independent definitions of extinction, the computational efficiency multi-sphere method is applied for particle clusters in absorbing host, and its accuracy is verified with the discrete dipole approximation method. It is well known that for small particles, the dependent scattering in transparent host always enhances the extinction of the cluster and the criterion for the ISA is nearly independent of the particle refractive index and particle size. We show, however, that when the host medium is absorbing, the dependent scattering between particles can lead to a decreased or even negative extinction, and thus the ISA criterion depends on the particle refractive index, size, and host medium absorption index. In this result, the generalized criteria for absorbing host media may differ significantly from the conventional ones for transparent host media. The results can provide guidance in solving problems related to light transfer in nanoparticle systems, particularly in the presence of absorption in the host medium.

1. Introduction

The phenomenon of radiative (light) transfer in nanoparticle systems is ubiquitous and plays a vital role in numerous applications, including climate science, ocean optics, remote sensing, biomedicine, color paints, solar energy utilization and radiative cooling [1,2,3,4,5,6,7,8,9,10,11,12]. Accurate prediction of the radiative properties of nanoparticle systems is central to such applications. An approximate method based on the radiative transfer equation (RTE) is one of the most used methods for calculating the radiative properties of nanoparticle systems [13,14]. Prior to solving the RTE, the scattering properties of particles should be determined. Most works usually calculate the scattering properties of particles based on the independent scattering approximation (ISA) [1,2,15,16]. Under the ISA, the scattering properties of an individual particle in the nanoparticle system are independent of the other particles [1,2]. While the ISA predictions are accurate when the interparticle separation is sufficiently large, they may deviate significantly from the actual properties when particles are in close proximity [2,15,17]. Therefore, the main task of the radiative property calculation is to decide whether the ISA is applicable for the specific practical problem at hand [18].
So far, plenty of works have been conducted to investigate the applicability conditions of the ISA (independent scattering criterion) [2,19,20]. Tien and coworkers [21,22,23] established a regime map for deciding whether dependent scattering can be neglected over a wide range of size parameters and volume fractions based on theoretical analysis and experimental research. Their results suggest that the ISA is satisfied if the clearance-to-wavelength ratio is larger than 0.5 (or the particle volume fraction is less than 0.006). Quirantes et al. [24] used the T-matrix method to calculate scattering cross sections of two-sphere clusters and used the results to determine the maximum distance required to produce interparticle light scattering interaction. Hohmann et al. [25] investigated the effects of multiple and dependent scattering on the Mueller matrix by comparing results of a polarization-sensitive radiative transfer solution to the Maxwell theory. Mishchenko et al. [26] studied the applicability conditions of the RTE based on comparisons of the experimentally measured Stokes reflection matrix for suspensions of microscopic latex particles in water with solutions of the vector RTE. Galy et al. [15] proposed a new criterion for the scattering cross section and the asymmetry factor based on electromagnetic scattering calculation of particle systems with up to eight spherical particles. However, previous studies typically assumed that particles are embedded in the nonabsorbing (transparent) host medium, neglecting the effect of host medium absorption on particle scattering [15,16,18,27,28].
Cases where the absorption of the host medium cannot be neglected are frequently encountered in practical applications. Examples include light scattering by cloud particles surrounded by water vapor in the atmosphere, light scattering by particles in polymers, biomedicine, and many others [4,29,30]. Neglecting host medium absorption can result in errors, because the host medium absorption nullifies the validity of many scattering quantities that are well defined for the nonabsorbing host medium [31]. To obtain appropriate scattering properties for absorbing host media, various approaches have been proposed from different perspectives. One widely used approach is based on the asymptotic form of the scattered field in the far-field zone [32,33,34]. Another widely used approach is based on the integration of the field on the surface of the particle [35,36]. Mishchenko [37,38,39] has developed an alternative approach using the vector volume integral equation, where the far-field extinction cross section and scattering matrix represent observable quantities and can be applied to the RTE. Using the above or other approaches, the optical properties of individual spherical particles [31,40,41,42,43,44,45,46,47,48] and individual aspherical [49,50,51,52] particles have been studied with consideration of the effect of host medium absorption on particle scattering. Additionally, light transfer in the nanoparticle system, which consists of monodisperse or polydisperse spheres randomly dispersed in an absorbing host medium, has also been investigated based on the RTE and the ISA [29,30,53]. However, the independent scattering criterion used in these studies is derived for nonabsorbing host media, and thus may be invalid for absorbing host media. Recently, multiple scattering of light by spherical particles embedded in absorbing host media has been investigated using the statistical approach, and several numerical examples have been presented to illustrate the spectral transmittance and reflectance of thin films containing a two-dimensional array of spherical particles [54].
Many aspects of light scattering by particles in absorbing host media have been studied; however, little work has focused on the effects of host medium absorption on the independent scattering criterion. Considering the widespread use of the ISA (and the RTE) and the presence of absorption in the host medium in many applications, it is necessary to explore the applicability conditions of the ISA for absorbing host media. In this paper, we aim at providing a comprehensive analysis of the applicability conditions of the ISA for absorbing host media. To accomplish this objective, we must obtain the scattering properties of a cluster of nanoparticles considering both the host medium absorption and dependent scattering effects. Using the far-field and distance-independent definitions of extinction [37] and the multi-sphere method (MSM) [55], we can obtain the scattering properties that take these two effects into account simultaneously. Note that the MSM is one of the most efficient numerically exact methods for solving the macroscopic Maxwell equations, which can take into account the interactions between particles. Although the MSM is typically used for nonabsorbing host media, it can be used for absorbing host media with appropriate modifications of the definition of the scattering quantities. With the help of the above methods, we first investigate the extinction properties of a cluster of particles embedded in absorbing host media considering the dependent scattering effect. Then, the effects of host medium absorption on the applicability conditions of the ISA are analyzed. Various factors, including particle refractive index, particle radius, and host medium absorption index, are discussed.

2. Methodology

In this section, we describe the MSM for the absorbing host medium and the cluster models for investigating the independent scattering criterion. The MSM has been well established, and several codes, such as the GMM code [56], the MSTM code [57], the FaSTMM code [58], and the CELES code [59], are freely available. The main idea of the MSM is to expand the scattered fields from all the spheres into a single expansion written about the origin of the cluster. This procedure is based on the addition theorem for vector spherical wave functions, which allows a scattered field from one sphere to be represented as an exciting field in another sphere [55]. It should be noted that to account for the host medium absorption effect, the far-field scattering quantities need to be redefined. In this work, the definition given by Mishchenko [37,60] is adopted because the corresponding scattering quantities can be applied to the RTE.

2.1. Expansion of the Scattered Field

We consider the model of a cluster consisting of N s nonoverlapping spheres embedded in an absorbing host medium with known complex refractive index m h = m h + i m h . Each sphere is characterized by its radius a i and complex refractive index m p , i = m p , i + i m p , i , for i = 1,2 N s . The central point of the cluster O serves as the origin of the primary coordinate system and the central point of the i th sphere O i serves as the origin of the i th coordinate system, as shown in Figure 1.
According to the MSM, the scattered field of the cluster E s can be represented by the superposition of scattered fields from each of the spheres in the cluster [55]
E s = i = 1 N s E s i .
The scattered field of the i th sphere E s i can be expanded in terms of vector spherical wave functions about its origin O i as
E s i = n = 1 L i m = n n a m n i N m n 3 k h r i , θ i , ϕ i + b m n i M m n 3 k h r i , θ i , ϕ i ,
where L i is the truncation order [61], k h = k h + i k h = 2 π m h / λ is the complex wavenumber in the host medium, λ is the wavelength in a vacuum, and a m n i and b m n i are the scattering expansion coefficients for sphere i . M m n υ and N m n υ represent the vector spherical wave functions and are given by [62]
M m n υ = 1 Γ m n × r u m n υ , N m n υ = 1 k h × M m n υ ,
where Γ m n = n n + 1 2 n + 1 n + m ! n m ! , u m n 1 ρ , θ , ϕ = j n ρ P n m cos θ e i m ϕ , and u m n 3 ρ , θ , ϕ = h n ρ P n m cos θ e i m ϕ . In addition, j n , h n and P n m denote the spherical Bessel function, the spherical Hankel function of the first kind, and the associated Legendre function, respectively [1,63,64]. In a similar fashion, the vector spherical harmonic expansion of the incident field is written about the origin of the i th sphere O i as [55]
E inc i = n = 1 L i m = n n p m n i N m n 1 k h r i , θ i , ϕ i + q m n i M m n 1 k h r i , θ i , ϕ i ,
where p m n i and q m n i are the incident expansion coefficients about the origin O i . Using the above equations and the addition theorem for vector spherical wave functions, the scattered field of the cluster as a whole can be expanded about the origin of the cluster O as [55]
E s = n = 1 L i m = n n a m n N m n 3 k h r , θ , ϕ + b m n M m n 3 k h r , θ , ϕ ,
where a m n and b m n are the scattering expansion coefficients of the cluster as a whole in the primary coordinate system. Note that the complex-valued wavevector k h does not affect the form of the formulas j n k h r and h n k h r [65]. Thus, the expressions of the fields for an absorbing medium do not differ from those for a nonabsorbing medium, except that the refractive index of the host medium is complex-valued.

2.2. Far-Field Extinction and Scattering of Sphere Clusters

As shown in Figure 2, we assume the incident light propagates along the z axis, and the plane containing the scattering direction n ^ s (or e ^ r ) and the incident direction n ^ i (or z ^ ) defines the scattering plane.
Following these conventions, the incident electric field is resolved as [1,64]
E i n c = E 0 i e i k h z = e ^ i E 0 i + e ^ i E 0 i e i k h z .
Similarly, the scattered field in the far-field region is resolved as [64]
E s = e i k h r r E 0 s = e i k h r r e ^ s E 0 s + e ^ s E 0 s ,
or as
E s = e i k h r r e ^ s θ E 0 s θ + e ^ s ϕ E 0 s ϕ ,
where e ^ s = e ^ s θ and e ^ s = e ^ s ϕ . Note that the wavenumber k h = k h + i k h is complex for the absorbing host medium. The relation between the scattered and incident fields can be represented by a 2 × 2 amplitude scattering matrix [64]
E 0 s E 0 s = s 11 n ^ s , n ^ i s 12 n ^ s , n ^ i s 21 n ^ s , n ^ i s 22 n ^ s , n ^ i E 0 i E 0 i .
With the elements of the amplitude scattering matrix known, the optical quantities of interest can be derived from them.
Using Mishchenko’s definition [37,60], the far-field and distance-independent extinction and scattering cross sections and the scattering phase function are written as
C ext = 4 π k h Im E 0 s n ^ i E 0 i E 0 i 2 ,
C sca = 2 π 0 π 1 2 s 11 2 + s 22 2 sin θ d θ ,
p θ = 4 π C s c a 1 2 s 11 2 + s 22 2 ,
where k h is the real part of the complex wavenumber in the absorbing host medium, Im represents the imaginary part of the argument, and the asterisk denotes the complex conjugate. Taking the extinction cross section as an example, the explicit expressions for incident x -polarized and y -polarized lights are written as
C ext , x = 2 π k h Re 1 k h n = 1 L i i n + 1 2 n + 1 a 1 n , x + b 1 n , x + a 1 n , x + b 1 n , x ,
C ext , y = 2 π k h Re 1 k h n = 1 L i i n + 1 i 2 n + 1 a 1 n , y + b 1 n , y + a 1 n , y b 1 n , y ,
where the subscripts x and y on a m n (or b m n ) denote the arguments calculated for incident x -polarized and y -polarized lights, respectively, and Re represents the real part of the argument. The formulas reduce to the conventional formulas in a nonabsorbing host medium [61]. For incident unpolarized light, the extinction cross section is simply given by C ext = C ext , x + C ext , y / 2 , and the corresponding extinction efficiency factor is given by Q ext = C ext / π a eff 2 , while a eff is the effective volume-equivalent-sphere radius. For monodispersed spheres, we have a eff = a i N s 1 / 3 .

2.3. Nanoparticle Cluster Model for Investigating the Independent Scattering Criterion

Similar to previous studies [15,16,18,27], the model of an imaginary spherical volume filled with N s randomly distributed identical spheres (see Figure 3) is used to explored the effect of host medium absorption on the scattering properties of a nanoparticle cluster. By keeping the particle radius fixed while scaling all particle coordinates, the clusters with different average interparticle clearances c a v e (or particle volume fraction f v ) are obtained. The average clearance c a v e is calculated based on the relative positions of particles.
If the interactions between particles can be neglected, the total extinction cross section of the N s -sphere cluster C ext is equal to N s times the corresponding single-particle extinction cross section C e x t s i n g l e . Therefore, in the independent scattering regime, the value of the nondimensional extinction cross section
η ext = C ext N s C e x t s i n g l e
must be equal to unity [27]. In practical applications, the independent scattering regime is usually considered to be reached when the extinction of the cluster deviates less than 5% from the summation of extinction by the individual particles (the value of η ext locates between 0.95 and 1.05). In this work, the ratio of average interparticle clearance to wavelength c = c a v e / λ is used to determine the ISA criterion, and the ISA is considered to be applicable when 0.95   η ext 1.05.
First, we investigate the extinction properties of sphere cluster embedded in absorbing host media with different particle radii a i and particle numbers N s for both fixed and random orientations. Then, the impact of host medium absorption on the applicability conditions of the ISA is analyzed, considering the effects of the particle refractive index and particle radius.

3. Validation

In this section, we validate the MSM by comparing it with the discrete dipole approximation (DDA) method [50,51] for absorbing host media. The extinction efficiency factor and the scattering phase function are used for validation in this work, because these two parameters represent observable quantities in the absorbing host medium. The clusters composed of N s = 2, 3 and 8 spheres are considered. The complex refractive indices of the spheres and the host medium are m p = 1.5 and m h = 1.0 + 0.02 i , respectively. The incident light propagates along the z axis, and the wavelength in a vacuum is λ = 2 π μm. Figure 4a–d shows the extinction efficiency factors Q ext of two-sphere clusters, calculated by the MSM and the DDA method, as a function of the ratio of interparticle clearance c to the radius of sphere 1, c / a 1 . N d represents the total number of discretized sub-volumes of the multiple-sphere cluster. Different numbers of discretized sub-volumes are used to ensure the accuracy of the DDA results. As shown, the extinction efficiency factors obtained by the MSM coincide very well with those obtained by the DDA method. Figure 4e,f shows the extinction efficiency factors Q ext of three-sphere and eight-sphere clusters as a function of sphere radius a i . Some minor differences are observed for large radius, but the differences become smaller with increasing the number of the discretized sub-volumes N d .
Figure 5 compares the scattering phase functions p θ calculated by the MSM and DDA method. It is seen that the MSM results are also in good agreement with the DDA results. In addition, the running times of the MSM and the DDA ( N d = 526,016) method for calculating the scattering properties of the eight-sphere cluster ( a i = 5.0 μm) are 0.463 s and 2423 s, respectively. Note that the calculations are performed on a 64-bit Windows 10 operating system with an Intel Xeon 6144 CPU using one core. Obviously, MSM is much more efficient in calculating particle clusters than the DDA method.

4. Results and Discussion

4.1. Extinction Properties of Sphere Clusters Embedded in Absorbing Host Media

In this subsection, we will investigate the effect of the host medium absorption on the extinction properties of multiple-sphere clusters. Figure 6 depicts the extinction efficiency factors Q ext of two-sphere clusters at a fixed orientation versus the ratio of interparticle gap to sphere radius c / a i with a i = 0.2, 1.0 and 2.0 μm, respectively. The incident light propagates (a)–(c) perpendicularly to and (d)–(f) along the line connecting the centers of the two spheres. The complex refractive indices of the particle and the host medium are m p = 1.59 and m h = 1.0 + i m h with m h = 0.0, 0.001, 0.01 and 0.02, respectively. The incident wavelength is λ = 2 π μm. For comparison, the extinction efficiency factors of the single sphere with radius equal to the cluster components are listed in Table 1. As shown, with enhancing the absorption of the host medium, the extinction efficiency factor of the single sphere decreases at a i = 0.2 and 1.0 μm, but increases at a i = 2.0 μm. A similar trend is found for the two-sphere cluster. In addition, for the case of a i = 0.2 μm and m h = 0.01 and 0.02, the extinction efficiency factor is negative for the single sphere, but may become positive for the two-sphere cluster. This phenomenon can be understood from the negative extinction and the dependent scattering effects. The concept of extinction in an absorbing host medium is the difference between two readings of a well-collimated detector in the absence and in the presence of particles [40,41]. Therefore, under the conditions where the particle is less attenuated than the host medium, the extinction efficiency factor of the particle is negative [45,66]. However, the dependent scattering effects may enhance the extinction of light by particles, so that the host medium can be less attenuated than the two-sphere cluster. As a consequence, the extinction efficiency factor is no longer negative. On the other hand, it is found in Figure 6f that when m h = 0.02, the amplitude of the interference oscillations increases exponentially with increasing c / a i . This phenomenon has been observed and discussed in previous studies of individual particles [40,41]. For a nonabsorbing particle embedded in an absorbing host medium, the exponentially increasing amplitude of the oscillations with increasing particle size is due to the fact that the directly transmitted light is no longer subject to the exponential attenuation over the path length given by the particle diameter [40,41]. With this explanation, the exponential increase in the amplitude of the oscillations of two-sphere clusters can also be understood.
Considering that the complex refractive index of a material typically varies with the wavelength, we also investigate the extinction efficiency factors Q ext of randomly oriented silica (SiO2) clusters embedded in polyethylene (PE), as presented in Figure 7a–d for different particle numbers N s , and in Figure 7d–f for different particle radii a i . The complex refractive indices of SiO2 and PE are taken from [67] and [68], respectively. The particle volume fraction is fixed at f v = 0.1. The extinction efficiency factors presented in the figures are averaged over 132 different orientations and 10 different configurations to obtain random-orientation averaging results. In practical applications, particles with radii in the order of a few nanometers to micrometers may be used [6,69]; however, we only focus on small particles with radii ranging from 40 to 160 nm. In addition, we consider the wavelength range of 3–12 µm because it is especially useful for designing functional coatings, and the resonance peak of SiO2 is also in this wavelength range [70,71]. As demonstrated in Figure 7, within the wavelength range of 3.0–8.0 μm, the host medium absorption can lead to negative extinction, which significantly deviates from the results for the nonabsorbing host medium. In contrast, within the wavelength range of 8.0–12.0 μm, the effect of host medium absorption on Q ext can be neglected due to the extremely small values of m h and the strong absorption of the particle. Moreover, the comparison between the clusters with different particle radii (see Figure 7d–f) reveals that the differences caused by neglecting the host medium absorption become more pronounced as the particle radius a i increases from 40 to 160 nm. Notably, the changes in the extinction will, in turn, affect the applicability conditions of the ISA, as we shall discuss in the following section.

4.2. Effect of Host Medium Absorption on the Dependent Scattering between Particles

In this subsection, we investigate the effect of host medium absorption on the applicability conditions of the independent scattering approximation (ISA). The ISA is widely used in many analyses involving particles [2]. Similar to previous studies [15,16,18,27], we use the model of a small volume filled with monodisperse spherical particles to investigate the applicability conditions of the ISA. Each cluster consists of N s = 27 identical particles. The complex refractive index of the host medium is m h = 1.5154 + i m h with m h ranging from 0.0 to 0.01, and the incident wavelength is λ = 4.0 μm, unless otherwise specified. Note that the value of 1.5154 is selected as the refractive index (real part of the complex refractive index) of polyethylene (PE) at λ = 4.0 μm. The particle radius a i ranges from 40 to 160 nm, which is small compared to the incident wavelength.
Figure 8 presents the color images of the nondimensional extinction cross section η ext as functions of relative refractive index m p , i / m h and the ratio of average interparticle clearance to wavelength c = c a v e / λ , with m p , i / m h varying from 0.8 to 2.0, and c varying from 0.055 to 4.91. The particle radius is a i = 40 nm, and the particle absorption index (imaginary part of the complex refractive index) is m p , i = 0.0. The white regions in the images represent the value of η ext lying between 0.95 and 1.05, indicating that the ISA is satisfied. It should be noted that Figure 8 does not encompass the instance where the complex refractive index of the host medium is equal to that of the particle ( m p , i / m h = 1.0 and m h = 0.0). As seen in Figure 8a, when the host medium is nonabsorbing ( m h = 0.0), the particle refractive index m p , i has a negligible impact on η ext and the applicability conditions of ISA. This observation agrees with previous results and can be attributed to the fact that for small particles, the phase shift is negligible and the scattering is isotropic [15]. Moreover, when particles are close to each other, the values of η ext are always greater than 1.0, suggesting that the interactions between particles enhance the particle extinction. In contrast, when the host medium is absorbing, as illustrated in Figure 8b–d, the values of η ext can be greater than 1.0, less than 1.0, or even negative (cause by the negative extinction), depending on m p , i / m h . This means that the interactions between particles may decrease particle extinction ( η ext < 1.0), which differs significantly from that observed in nonabsorbing host media. Meanwhile, it is found that as m p , i / m h decreases from 2.0 to 0.8, η ext may display a sharp transition from a value of greater than 1.0 to a value of less than 0.0 (negative value). For example, when m h = 0.001 and c = 0.055 (see Figure 8c), the values of η ext at m p , i / m h = 1.56 and 1.58 are −396.84 and 113.04, respectively. When m h = 0.01 and c = 0.055 (see Figure 8d), the value of η ext is 15.95 at m p , i / m h = 1.66, but decreases to −123.29 at m p , i / m h = 1.64.
The change in η ext will inevitably affect the critical c that makes the ISA applicable (0.95 η ext 1.05). To better illustrate the effect of host medium absorption on the ISA criterion, we also plot the nondimensional extinction cross section η ext as a function of the ratio of average interparticle clearance to wavelength c for a i = 40 nm, m p , i = 1.38, 2.1216 and 3.0308 ( m p , i / m h = 0.911, 1.4 and 2.0), and (a) m h = 0.0, (b) m h = 0.0001, (c) m h = 0.001 and (d) m h = 0.01, respectively, as shown in Figure 9. The blue region indicates that the ISA is satisfied for all the three m p , i / m h considered. Two dashed lines corresponding to η ext = 0.95 and 1.05, respectively, are also plotted for comparison. As expected, when the host medium is nonabsorbing (see Figure 9a), the particle refractive index m p , i has little impact on the critical c . In this case, the ISA is applicable regardless of the particle refractive index, provided that c 2.52. Nevertheless, it can be seen from Figure 9b–d that when the host medium is absorbing, the ISA criterion is correlated with particle refractive index m p , i , and even a minor change in the host medium absorption index m h may result in considerable changes in the critical c . For instance, when m p , i = 1.38 ( m p , i / m h = 0.911), the criteria ensuring that the ISA is accurate within 5 percent are approximately c = 2.52, 0.61, 0.22 and 0.055, for m h = 0.0, 0.0001, 0.001 and 0.01, respectively. When m p , i = 3.0308 ( m p , i / m h = 2.0), the criteria are approximately c = 2.52, 2.1, 1.64 and 0.66, for m h = 0.0, 0.0001, 0.001 and 0.01, respectively. We note that when m p , i = 1.38 and 3.0308, the critical c decreases with the increase in m h , which agrees with the traditional view that absorption suppresses dependent scattering [29]. However, when m p , i = 2.1216 ( m p , i / m h = 1.4), as m h increases from 0.0 to 0.0001, the critical c increases from 2.52 to 4.1. The larger critical c means that the intensity attenuation in the absorbing host is not the only factor to influence the dependent scattering.
The preceding discussions have shown that the host medium absorption can either decrease or increase the critical c . This is reasonable because host medium absorption can affect the applicability conditions of the ISA from two perspectives. On the one hand, the intensity of light is attenuated as it travels through an absorbing medium. In this regard, the host medium absorption suppresses the dependent scattering and thereby renders the ISA applicable at smaller c . On the other hand, the host medium absorption also influences the extinction efficiency factor of individual particles (or particle clusters), which, in turn, affects the critical c . These two different effects result in the dependence of the critical c on particle refractive index m p , i in an absorbing host medium. Although the variation in the critical c with m p , i is complex. It can be seen from Figure 8 and Figure 9 that in regimes where the value of η ext changes from larger than 1.0 to less than 0.0 (negative), a larger c may be required to ensure the accuracy of the ISA.
We also investigate particle clusters with different particle radii. The color images of the nondimensional extinction cross sections η ext as functions of m p , i / m h and the ratio of average interparticle clearance to wavelength c are illustrated in Figure 10 for a i = 80 nm and in Figure 11 for a i = 160 nm. It is obvious that both the nondimensional extinction cross section η ext and the critical c are almost independent of the particle refractive index m p , i when the host medium is nonabsorbing, but are correlated with the particle refractive index m p , i when the host medium is absorbing. A comparison among particles with different radii demonstrates that in absorbing host media, the particle radius has significant impacts on both η ext and the critical c . Furthermore, when the host medium is absorbing, it may need a larger c to ensure the accuracy of the ISA in regimes where the value of η ext changes from larger than 1.0 to less than 0.0, which is similar to that observed for the case of a i = 40 nm.
Overall, in absorbing host media, the applicability conditions of the ISA are affected by the particle radius a i , particle refractive index m p , i (or m p , i / m h ) and host medium absorption index m h , making it difficult to establish an accurate and comprehensive criterion of ISA based only on limited calculation results. However, for all the cases considered in Figure 8, Figure 9, Figure 10 and Figure 11, the criteria c 2.52 for nonabsorbing host media and c 4.8 for absorbing host media ensure that the accuracy of the ISA is kept within 5 percent.
Additionally, to illustrate the effect of host medium absorption on the ISA criterion in practical applications, we also investigate the case of silica (SiO2) particles embedded in polyethylene (PE) and its nonabsorbing counterpart. Three incident wavelengths of λ = 3.35, 5.0 and 10.0 μm are considered. The nondimensional extinction cross sections η ext of multiple-sphere clusters composed of N s = 27 identical particles as a function of the ratio of average interparticle clearance to wavelength c are presented in Figure 12. The radii of the particles are a i = 40, 80 and 160 nm, respectively. The results for particles embedded in PE are calculated using the actual complex refractive indices, whereas the results for particles embedded in nonabsorbing PE are calculated by assuming m h = 0.0.
When λ = 3.35 μm, the host medium exhibits strong absorption with m h = 0.014027, and the value of m p , i / m h is 0.9908. According to above discussions, under this condition, the strong absorption of the host medium should significantly suppress the interaction between particles, allowing the ISA to be applicable for very small values of c . As a result, the critical c for nonabsorbing host media should also be larger than the critical c for absorbing host media. These assertions are confirmed by the results depicted in Figure 12a–c. It can be seen from Figure 12a–c that the values of η ext always locate between 0.95 and 1.05, suggesting that the accuracy of the ISA is always kept within 5 percent.
When λ = 5.0 μm, the values of m h and m p , i / m h are 0.0063785 and 0.8814, respectively. From Figure 12d–f, it can be easily observed that the values of η ext can be less than 1.0 or less than 0.0, which will no doubt affect the applicability conditions of the ISA. Therefore, the applicability conditions of the ISA must be carefully considered in this case. These observations are consistent with the above finding that in this regime, the interactions between particles may result in a decreased or even negative particle extinction cross section.
When λ = 10.0 μm, the host medium exhibits extremely weak absorption with m h = 0.000023606 and a value of m p , i / m h = 1.6463. We can thus expect that the effects of host medium absorption on η ext and the applicability conditions of the ISA are negligible. These predictions are supported by the observations in Figure 12g–i, where the values of η ext obtained for nonabsorbing host media are almost identical to those obtained for absorbing host media and the ISA criteria for nonabsorbing host media can be regarded as the criteria for absorbing host media.
It is also found from Figure 12 that in nonabsorbing host media, the critical c changes with a i and m p , i / m h , which is inconsistent with our previous conclusion that the ISA criteria for nonabsorbing host media are independent of particle radius and particle refractive index. This discrepancy stems from the fact that the absorption index of SiO2 particle is nonzero at the wavelengths considered. The issue of the effect of particle absorption on ISA criteria is beyond the scope of the present work. Interested readers can refer to Ma et al. [16] for detailed discussions.
In summary, for small particles, the host medium absorption has a significant impact on the dependent scattering between particles. The generalized criteria for absorbing host media depend on the particle refractive index, size, and the host medium absorption index and may differ significantly from the conventional ones for nonabsorbing host media. Moreover, since it is impossible to encompass all factors, many relevant factors, such as polydispersity, spatial distributions and large particles, are not considered in this paper. These factors may also affect applicability conditions of the ISA. For example, the effects of host medium absorption on dependent scattering between particles are different for small particles and large particles (see Supplementary Materials, Figure S1). Thus, a lot of work remains to be done to establish an accurate and comprehensive criterion applicable to both absorbing and nonabsorbing host media.

5. Conclusions

The aim of this paper was to investigate the extinction properties and the applicability conditions of the independent scattering approximation (ISA) for multiple nanoparticles embedded in absorbing host media. The multi-sphere method (MSM) is employed to simultaneously consider the effects of host medium absorption and dependent scattering.
First, the effects of host medium absorption on the extinction properties of particles are investigated with consideration of the dependent scattering effect. The results for two-sphere clusters at a fixed orientation indicate that enhancing the absorption of the host medium has a more pronounced effect where the incident light propagates along the line connecting the centers of the two spheres. Moreover, the results for randomly oriented multiple-sphere silica clusters (with particle radius ranging from 40 to 160 nm and particle volume fraction of 0.1) embedded in polyethylene suggest that the host medium absorption can lead to negative extinction within the wavelength range of 3.0–8.0 μm, which differs significantly from the results for nonabsorbing host media.
Then, the effects of host medium absorption on the interactions between (nonabsorbing) particles and the applicability conditions of the ISA are investigated for a wide range of parameters. The particle radius ranges from 40 to 160 nm, which is small compared to the incident wavelength. The ratio of particle refractive index to medium refractive index ranges from 0.8 to 2.0, and the host medium absorption index ranges from 0 to 0.01. When the host medium is nonabsorbing, the nondimensional extinction cross section η ext and the ISA criteria are almost independent of the particle refractive index and particle radius, and the interactions between particles always enhance particle extinction. In contrast, when the host medium is absorbing, the nondimensional extinction cross section and the ISA criteria correlate with particle size, particle refractive index and host medium absorption index, and the interactions between particles may suppress particle extinction and even lead to negative extinction. Moreover, it is found that in regimes where the value of η ext changes from larger than 1.0 to less than 0.0 (negative), a larger c (the ratio of interparticle clearance to wavelength) may be required to ensure the accuracy of the ISA, which is inconsistent with the traditional view that absorption suppresses the dependent scattering. In this result, the criteria for absorbing host media differ significantly from the conventional ones for transparent host media. Specifically, for all the cases considered in this work, the criteria c 2.52 for nonabsorbing host media and c 4.8 for absorbing host media ensure that the accuracy of the ISA is kept within 5 percent.

Supplementary Materials

The following supporting information can be downloaded at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/photonics10070782/s1, Figure S1: The nondimensional extinction cross section η ext of 27-sphere cluster versus the ratio of average interparticle clearance to wavelength c with incident wavelength of λ = 4.0 μm, particle radius of a i = 500, 750 and 1000 nm, and particle refractive index of m p , i = 1.38, 2.1216 and 3.0308 ( m p , i / m h = 0.911, 1.4, and 2.0). The complex refractive index of the host medium is m h = 1.5154 + i m h with m h = 0.0, and 0.001. (a) a i = 500 nm, m h = 0.0, (b) a i = 500 nm, m h = 0.001, (c) a i = 750 nm, m h = 0.0, (d) a i = 750 nm, m h = 0.001, (e) a i = 1000 nm, m h = 0.0, and (f) a i = 1000 nm, m h = 0.001. The blue region represents that the ISA is satisfied for all the three m p , i / m h considered. The results are averaged over 132 different orientations and 10 different configurations.

Author Contributions

Conceptualization, J.Z., S.Z. and L.L.; methodology, J.Z. and S.Z.; software, J.Z.; investigation, J.Z., S.Z. and L.L.; data curation, J.Z., S.Z. and L.L.; writing—original draft preparation, J.Z.; writing—review and editing, J.Z., S.Z. and L.L.; supervision, S.Z. and L.L.; funding acquisition, S.Z. and L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant 52076123) and the China Postdoctoral Science Foundation (grant 2021M691906).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data underlying the results presented in this paper are not publicly available at this time, but may be obtained from the authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bohren, C.F.; Huffman, D.R. Absorption and Scattering of Light by Small Particles; Wiley: New York, NY, USA, 1983. [Google Scholar]
  2. Howell, J.R.; Mengüc, M.P.; Siegel, R. Thermal Radiation Heat Transfer; CRC Press: Boca Raton, FL, USA, 2015. [Google Scholar]
  3. Theobald, D.; Yu, S.; Gomard, G.; Lemmer, U. Design of Selective Reflectors Utilizing Multiple Scattering by Core–Shell Nanoparticles for Color Conversion Films. ACS Photonics 2020, 7, 1452–1460. [Google Scholar] [CrossRef]
  4. Gentle, A.R.; Smith, G.B. Radiative Heat Pumping from the Earth Using Surface Phonon Resonant Nanoparticles. Nano Lett. 2010, 10, 373–379. [Google Scholar] [CrossRef] [PubMed]
  5. Schertel, L.; Siedentop, L.; Meijer, J.-M.; Keim, P.; Aegerter, C.M.; Aubry, G.J.; Maret, G. The Structural Colors of Photonic Glasses. Adv. Opt. Mater. 2019, 7, 1900442. [Google Scholar] [CrossRef] [Green Version]
  6. Yalçın, R.A.; Blandre, E.; Joulain, K.; Drévillon, J. Colored Radiative Cooling Coatings with Nanoparticles. ACS Photonics 2020, 7, 1312–1322. [Google Scholar] [CrossRef]
  7. Hwang, V.; Stephenson, A.B.; Barkley, S.; Brandt, S.; Xiao, M.; Aizenberg, J.; Manoharan, V.N. Designing angle-independent structural colors using Monte Carlo simulations of multiple scattering. Proc. Natl. Acad. Sci. USA 2021, 118, e2015551118. [Google Scholar] [CrossRef]
  8. Boutghatin, M.; Assaf, S.; Pennec, Y.; Carette, M.; Thomy, V.; Akjouj, A.; Djafari Rouhani, B. Impact of SiO2 Particles in Polyethylene Textile Membrane for Indoor Personal Heating. Nanomaterials 2020, 10, 1968. [Google Scholar] [CrossRef]
  9. Zhang, C.; Zhang, J.; Wu, X.; Huang, M. Numerical analysis of light reflection and transmission in poly-disperse sea fog. Opt. Express 2020, 28, 25410–25430. [Google Scholar] [CrossRef]
  10. Ma, L.; Hu, L.; Jia, C.; Wang, C.; Liu, L. Quantitative Evaluation of the Phase Function Effects on Light Scattering and Radiative Transfer in Dispersed Systems. Photonics 2022, 9, 584. [Google Scholar] [CrossRef]
  11. Ma, L.; Hu, K.; Wang, C.; Yang, J.-Y.; Liu, L. Prediction and Inverse Design of Structural Colors of Nanoparticle Systems via Deep Neural Network. Nanomaterials 2021, 11, 3339. [Google Scholar] [CrossRef]
  12. Fei, T.; Lin, L.; Li, X.; Yang, J.-Y.; Zhao, J.; Liu, L. Modeling Effect of Bubbles on Time-Dependent Radiation Transfer of Microalgae in a Photobioreactor for Carbon Dioxide Fixation. Photonics 2022, 9, 864. [Google Scholar] [CrossRef]
  13. Mishchenko, M.I. Electromagnetic Scattering by Particles and Particle Groups: An Introduction; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar] [CrossRef] [Green Version]
  14. Penttilä, A.; Markkanen, J.; Väisänen, T.; Räbinä, J.; Yurkin, M.A.; Muinonen, K. How much is enough? The convergence of finite sample scattering properties to those of infinite media. J. Quant. Spectrosc. Radiat. Transf. 2021, 262, 107524. [Google Scholar] [CrossRef]
  15. Galy, T.; Huang, D.; Pilon, L. Revisiting independent versus dependent scattering regimes in suspensions or aggregates of spherical particles. J. Quant. Spectrosc. Radiat. Transf. 2020, 246, 106924. [Google Scholar] [CrossRef]
  16. Ma, L.; Zhai, J.; Wang, C. Investigation of the single scattering approximation through direct electromagnetic scattering simulation. OSA Contin. 2021, 4, 2496–2509. [Google Scholar] [CrossRef]
  17. Chen, J.; Zhao, C.Y.; Wang, B.X. Effect of nanoparticle aggregation on the thermal radiation properties of nanofluids: An experimental and theoretical study. Int. J. Heat Mass Transfer 2020, 154, 119690. [Google Scholar] [CrossRef]
  18. Aoyu, Z.; Fuqiang, W.; Ziming, C.; Huaxu, L.; Xuhang, S. Radiative property investigation of dispersed particulate medium with the consideration of non-uniform particle size distribution and dependent scattering effects. Int. J. Heat Mass Transfer 2022, 186, 122488. [Google Scholar] [CrossRef]
  19. Kokhanovsky, A. Springer Series in Light Scattering: Volume 1: Multiple Light Scattering, Radiative Transfer and Remote Sensing; Springer: Cham, Switzerland, 2018. [Google Scholar]
  20. Kokhanovsky, A.A. Optics of Light Scattering Media: Problems and Solutions; Springer: London, UK, 2001. [Google Scholar]
  21. Cartigny, J.D.; Yamada, Y.; Tien, C.L. Radiative Transfer With Dependent Scattering by Particles: Part 1—Theoretical Investigation. J. Heat Transfer 1986, 108, 608–613. [Google Scholar] [CrossRef]
  22. Yamada, Y.; Cartigny, J.D.; Tien, C.L. Radiative Transfer With Dependent Scattering by Particles: Part 2—Experimental Investigation. J. Heat Transfer 1986, 108, 614–618. [Google Scholar] [CrossRef]
  23. Drolen, B.L.; Tien, C.L. Independent and dependent scattering in packed-sphere systems. J. Thermophys. Heat Transfer 1987, 1, 63–68. [Google Scholar] [CrossRef]
  24. Quirantes, A.; Arroyo, F.; Quirantes-Ros, J. Multiple Light Scattering by Spherical Particle Systems and Its Dependence on Concentration: A T-Matrix Study. J. Colloid Interface Sci. 2001, 240, 78–82. [Google Scholar] [CrossRef] [Green Version]
  25. Hohmann, A.; Voit, F.; Schäfer, J.; Kienle, A. Multiple scattering of polarized light: Influence of absorption. Phys. Med. Biol. 2014, 59, 2583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Mishchenko, M.I.; Goldstein, D.H.; Chowdhary, J.; Lompado, A. Radiative transfer theory verified by controlled laboratory experiments. Opt. Lett. 2013, 38, 3522–3525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Mishchenko, M.I.; Liu, L.; Videen, G. Conditions of applicability of the single-scattering approximation. Opt. Express 2007, 15, 7522–7527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Yalcin, R.A.; Lee, T.; Kashanchi, G.N.; Markkanen, J.; Martinez, R.; Tolbert, S.H.; Pilon, L. Dependent Scattering in Thick and Concentrated Colloidal Suspensions. ACS Photonics 2022, 9, 3318–3332. [Google Scholar] [CrossRef]
  29. Mishchenko, M.I.; Dlugach, J.M. Multiple scattering of polarized light by particles in an absorbing medium. Appl. Opt. 2019, 58, 4871–4877. [Google Scholar] [CrossRef]
  30. Ma, L.X.; Xie, B.W.; Wang, C.C.; Liu, L.H. Radiative transfer in dispersed media: Considering the effect of host medium absorption on particle scattering. J. Quant. Spectrosc. Radiat. Transf. 2019, 230, 24–35. [Google Scholar] [CrossRef]
  31. Ivanenko, Y.; Gustafsson, M.; Nordebo, S. Optical theorems and physical bounds on absorption in lossy media. Opt. Express 2019, 27, 34323–34342. [Google Scholar] [CrossRef] [Green Version]
  32. Mundy, W.C.; Roux, J.A.; Smith, A.M. Mie scattering by spheres in an absorbing medium. J. Opt. Soc. Am. 1974, 64, 1593–1597. [Google Scholar] [CrossRef]
  33. Chýlekt, P. Light scattering by small particles in an absorbing medium. J. Opt. Soc. Am. 1977, 67, 561–563. [Google Scholar] [CrossRef]
  34. Quinten, M.; Rostalski, J. Lorenz-Mie Theory for Spheres Immersed in an absorbing host medium. Part. Part. Syst. Charact. 1996, 13, 89–96. [Google Scholar] [CrossRef]
  35. Sudiarta, I.W.; Chylek, P. Mie-scattering formalism for spherical particles embedded in an absorbing medium. J. Opt. Soc. Am. A 2001, 18, 1275–1278. [Google Scholar] [CrossRef]
  36. Fu, Q.; Sun, W. Mie theory for light scattering by a spherical particle in an absorbing medium. Appl. Opt. 2001, 40, 1354–1361. [Google Scholar] [CrossRef] [PubMed]
  37. Mishchenko, M.I. Electromagnetic scattering by a fixed finite object embedded in an absorbing medium. Opt. Express 2007, 15, 13188–13202. [Google Scholar] [CrossRef] [PubMed]
  38. Mishchenko, M.I. Multiple scattering by particles embedded in an absorbing medium. 1. Foldy—Lax equations, order-of-scattering expansion, and coherent field. Opt. Express 2008, 16, 2288–2301. [Google Scholar] [CrossRef] [Green Version]
  39. Mishchenko, M.I. Multiple scattering by particles embedded in an absorbing medium. 2. Radiative transfer equation. J. Quant. Spectrosc. Radiat. Transf. 2008, 109, 2386–2390. [Google Scholar] [CrossRef]
  40. Mishchenko, M.I.; Videen, G.; Yang, P. Extinction by a homogeneous spherical particle in an absorbing medium. Opt. Lett. 2017, 42, 4873–4876. [Google Scholar] [CrossRef] [Green Version]
  41. Mishchenko, M.I.; Dlugach, J.M. Scattering and extinction by spherical particles immersed in an absorbing host medium. J. Quant. Spectrosc. Radiat. Transf. 2018, 211, 179–187. [Google Scholar] [CrossRef] [Green Version]
  42. Mishchenko, M.I.; Dlugach, J.M. Plasmon resonances of metal nanoparticles in an absorbing medium. OSA Contin. 2019, 2, 3415–3421. [Google Scholar] [CrossRef]
  43. Peck, R.L.; Brolo, A.G.; Gordon, R. Absorption leads to narrower plasmonic resonances. J. Opt. Soc. Am. B 2019, 36, F117–F122. [Google Scholar] [CrossRef]
  44. Zhai, J.; Ma, L.; Xu, W.; Liu, L. Effect of host medium absorption on the radiative properties of dispersed media consisting of optically soft particles. J. Quant. Spectrosc. Radiat. Transf. 2020, 254, 107206. [Google Scholar] [CrossRef]
  45. Zhang, S.; Zhang, W.; Liu, L. Scattering by a charged sphere embedded in an absorbing medium. J. Quant. Spectrosc. Radiat. Transf. 2020, 246, 106908. [Google Scholar] [CrossRef]
  46. Khlebtsov, N.G. Extinction, absorption, and scattering of light by plasmonic spheres embedded in an absorbing host medium. Phys. Chem. Chem. Phys. 2021, 23, 23141–23157. [Google Scholar] [CrossRef] [PubMed]
  47. Peck, R.; Khademi, A.; Ren, J.; Hughes, S.; Brolo, A.G.; Gordon, R. Plasmonic linewidth narrowing by encapsulation in a dispersive absorbing material. Phys. Rev. Res. 2021, 3, 013014. [Google Scholar] [CrossRef]
  48. Zhang, S.; Dong, J.; Zhang, W.; Luo, M.; Liu, L. Extinction by plasmonic nanoparticles in dispersive and dissipative media. Opt. Lett. 2022, 47, 5577–5580. [Google Scholar] [CrossRef] [PubMed]
  49. Zhang, S.; Zhang, W.; Dong, J.; Liu, L. Optical theorem of an infinite circular cylinder in weakly absorbing media. Phys. Rev. A 2022, 105, 023516. [Google Scholar] [CrossRef]
  50. Dong, J.; Zhang, W.; Liu, L. Discrete dipole approximation method for electromagnetic scattering by particles in an absorbing host medium. Opt. Express 2021, 29, 7690–7705. [Google Scholar] [CrossRef] [PubMed]
  51. Yurkin, M.A.; Moskalensky, A.E. Open-source implementation of the discrete-dipole approximation for a scatterer in an absorbing host medium. J. Phys. Conf. Ser. 2021, 2015, 012167. [Google Scholar] [CrossRef]
  52. Khlebtsov, N.G. Extinction and scattering of light by nonspherical plasmonic particles in absorbing media. J. Quant. Spectrosc. Radiat. Transf. 2022, 280, 108069. [Google Scholar] [CrossRef]
  53. Wang, C.C.; Ma, L.X. Effect of host medium absorption on polarized radiative transfer in dispersed media. Appl. Opt. 2019, 58, 7157–7164. [Google Scholar] [CrossRef]
  54. Loiko, N.A.; Miskevich, A.A.; Loiko, V.A. Optical characteristics of a monolayer of identical spherical particles in an absorbing host medium. J. Opt. Soc. Am. A 2023, 40, 682–691. [Google Scholar] [CrossRef]
  55. Hergert, W.; Wriedt, T. The Mie Theory: Basics and Applications; Springer: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  56. Xu, Y. Electromagnetic scattering by an aggregate of spheres. Appl. Opt. 1995, 34, 4573–4588. [Google Scholar] [CrossRef]
  57. Mackowski, D.W.; Mishchenko, M.I. A multiple sphere T-matrix Fortran code for use on parallel computer clusters. J. Quant. Spectrosc. Radiat. Transf. 2011, 112, 2182–2192. [Google Scholar] [CrossRef]
  58. Markkanen, J.; Yuffa, A.J. Fast superposition T-matrix solution for clusters with arbitrarily-shaped constituent particles. J. Quant. Spectrosc. Radiat. Transf. 2017, 189, 181–188. [Google Scholar] [CrossRef] [Green Version]
  59. Egel, A.; Pattelli, L.; Mazzamuto, G.; Wiersma, D.S.; Lemmer, U. CELES: CUDA-accelerated simulation of electromagnetic scattering by large ensembles of spheres. J. Quant. Spectrosc. Radiat. Transf. 2017, 199, 103–110. [Google Scholar] [CrossRef] [Green Version]
  60. Mishchenko, M.I.; Yang, P. Far-field Lorenz–Mie scattering in an absorbing host medium: Theoretical formalism and FORTRAN program. J. Quant. Spectrosc. Radiat. Transf. 2018, 205, 241–252. [Google Scholar] [CrossRef]
  61. Mackowski, D.W. Calculation of total cross sections of multiple-sphere clusters. J. Opt. Soc. Am. A 1994, 11, 2851–2861. [Google Scholar] [CrossRef]
  62. Mackowski, D.W. Calculation of the Scattering Properties for a Cluster of Spheres. User Guide Accompanying the SCSMFO.FOR Code. 1999. Available online: ftp://ftp.eng.auburn.edu/pub/dmckwski/scatcodes/scsmfo.ps (accessed on 5 June 2023).
  63. Mackowski, D.W. Analysis of radiative scattering for multiple sphere configurations. Proc. Math. Phys. Eng. Sci. 1991, 433, 599–614. [Google Scholar]
  64. Tsang, L.; Kong, J.A.; Ding, K.H. Scattering of Electromagnetic Waves: Theories and Applications; John Wiley & Sons: New York, NY, USA, 2000. [Google Scholar]
  65. Arfken, G.B.; Weber, H.J.; Harris, F.E. Mathematical Methods for Physicists; Academic Press: Boston, MA, USA, 2013. [Google Scholar]
  66. Bohren, C.F.; Gilra, D.P. Extinction by a spherical particle in an absorbing medium. J. Colloid Interface Sci. 1979, 72, 215–221. [Google Scholar] [CrossRef]
  67. Palik, E.D. Handbook of Optical Constants of Solids; Academic Press: San Diego, CA, USA, 1998. [Google Scholar]
  68. David, M.; Disnan, D.; Lardschneider, A.; Wacht, D.; Hoang, H.T.; Ramer, G.; Detz, H.; Lendl, B.; Schmid, U.; Strasser, G.; et al. Structure and mid-infrared optical properties of spin-coated polyethylene films developed for integrated photonics applications. Opt. Mater. Express 2022, 12, 2168–2180. [Google Scholar] [CrossRef]
  69. Chen, M.; Pang, D.; Chen, X.; Yan, H. Investigating the effective radiative cooling performance of random dielectric microsphere coatings. Int. J. Heat Mass Transfer 2021, 173, 121263. [Google Scholar] [CrossRef]
  70. Baneshi, M.; Maruyama, S.; Komiya, A. Infrared Radiative Properties of Thin Polyethylene Coating Pigmented With Titanium Dioxide Particles. J. Heat Transfer 2009, 132, 023306. [Google Scholar] [CrossRef]
  71. Wang, G.-H.; Zhang, Y.; Zhang, D.-H.; Fan, J.-P. Design and calculation of low infrared transmittance and low emissivity coatings for heat radiative applications. Int. J. Miner. Metall. Mater. 2012, 19, 179–184. [Google Scholar] [CrossRef]
Figure 1. Illustration of the coordinate systems. The central point of the cluster O serves as the origin of the primary coordinate system, and the central point of the i th ( j th ) sphere O i ( O j ) serves as the origin of the i th ( j th ) coordinate system. R i j is the distance from O j to O i , and Θ i j and Φ i j are the directions from O j to O i .
Figure 1. Illustration of the coordinate systems. The central point of the cluster O serves as the origin of the primary coordinate system, and the central point of the i th ( j th ) sphere O i ( O j ) serves as the origin of the i th ( j th ) coordinate system. R i j is the distance from O j to O i , and Θ i j and Φ i j are the directions from O j to O i .
Photonics 10 00782 g001
Figure 2. Schematic of light scattering by particles when the incident direction is along the + z axis. The scattering plane contains the scattering direction n ^ s (or e ^ r ) and the incident direction n ^ i (or z ^ ). The scattered field is resolved into components parallel ( e ^ s ) and perpendicular ( e ^ s ) to the scattering plane, respectively.
Figure 2. Schematic of light scattering by particles when the incident direction is along the + z axis. The scattering plane contains the scattering direction n ^ s (or e ^ r ) and the incident direction n ^ i (or z ^ ). The scattered field is resolved into components parallel ( e ^ s ) and perpendicular ( e ^ s ) to the scattering plane, respectively.
Photonics 10 00782 g002
Figure 3. Spherical volume elements filled with N s = 27 and N s = 100 randomly distributed identical spheres.
Figure 3. Spherical volume elements filled with N s = 27 and N s = 100 randomly distributed identical spheres.
Photonics 10 00782 g003
Figure 4. Comparison of MSM and DDA method extinction efficiency factors Q ext for (ad) two-sphere clusters with different ratios of interparticle clearance to radius of sphere 1, c / a 1 , and (e) three-sphere and (f) eight-sphere clusters with different sphere radii a i . The complex refractive indices of the spheres and the host medium are m p = 1.5 and m h = 1.0 + 0.02 i , respectively. N d is the number of the discretized sub-volumes.
Figure 4. Comparison of MSM and DDA method extinction efficiency factors Q ext for (ad) two-sphere clusters with different ratios of interparticle clearance to radius of sphere 1, c / a 1 , and (e) three-sphere and (f) eight-sphere clusters with different sphere radii a i . The complex refractive indices of the spheres and the host medium are m p = 1.5 and m h = 1.0 + 0.02 i , respectively. N d is the number of the discretized sub-volumes.
Photonics 10 00782 g004
Figure 5. Comparison of MSM and DDA method scattering phase functions p θ for (ad) two-sphere clusters, (e) three-sphere cluster and (f) eight-sphere cluster. The complex refractive indices of the particles and the host medium are m p = 1.5 and m h = 1.0 + 0.02 i , respectively.
Figure 5. Comparison of MSM and DDA method scattering phase functions p θ for (ad) two-sphere clusters, (e) three-sphere cluster and (f) eight-sphere cluster. The complex refractive indices of the particles and the host medium are m p = 1.5 and m h = 1.0 + 0.02 i , respectively.
Photonics 10 00782 g005
Figure 6. The extinction efficiency factors Q ext of two-sphere clusters versus the ratio of interparticle gap to sphere radius c / a i with a i = 0.2, 1.0 and 2.0 μm, respectively. The incident light (λ = 2 π μm) propagates (ac) perpendicularly to and (df) along the line connecting the centers of the two spheres. The complex refractive indices of the particle and the host medium are m p = 1.59 and m h = 1.0 + i m h with m h = 0.0, 0.001, 0.01 and 0.02, respectively.
Figure 6. The extinction efficiency factors Q ext of two-sphere clusters versus the ratio of interparticle gap to sphere radius c / a i with a i = 0.2, 1.0 and 2.0 μm, respectively. The incident light (λ = 2 π μm) propagates (ac) perpendicularly to and (df) along the line connecting the centers of the two spheres. The complex refractive indices of the particle and the host medium are m p = 1.59 and m h = 1.0 + i m h with m h = 0.0, 0.001, 0.01 and 0.02, respectively.
Photonics 10 00782 g006
Figure 7. The extinction efficiency factors Q ext of SiO2 nanoparticle clusters embedded in PE and its nonabsorbing counterpart with different particle numbers N s and particle radii a i . Each cluster consists of N s identical particles with a particle volume fraction of f v = 0.1. (a) N s = 8 and a i = 40 nm, (b) N s = 16 and a i = 40 nm, (c) N s = 27 and a i = 40 nm, (d) N s = 100 a i = 40 nm, (e) N s = 100 and a i = 80 nm, and (f) N s = 100 and a i = 160 nm. The results are averaged over 132 different orientations and 10 different configurations.
Figure 7. The extinction efficiency factors Q ext of SiO2 nanoparticle clusters embedded in PE and its nonabsorbing counterpart with different particle numbers N s and particle radii a i . Each cluster consists of N s identical particles with a particle volume fraction of f v = 0.1. (a) N s = 8 and a i = 40 nm, (b) N s = 16 and a i = 40 nm, (c) N s = 27 and a i = 40 nm, (d) N s = 100 a i = 40 nm, (e) N s = 100 and a i = 80 nm, and (f) N s = 100 and a i = 160 nm. The results are averaged over 132 different orientations and 10 different configurations.
Photonics 10 00782 g007
Figure 8. Color images of the nondimensional extinction cross sections η ext of 27-sphere clusters versus m p , i / m h and the ratio of average interparticle clearance to wavelength c . The results are averaged over 132 different orientations and 10 different configurations. The radius of the particle is a i = 40 nm, and the incident wavelength is λ = 4.0 μm. The complex refractive index of the host medium is m h = 1.5154 + i m h with (a) m h = 0.0, (b) m h = 0.0001, (c) m h = 0.001 and (d) m h = 0.01. The white regions represent the ISA being satisfied ( 0.95 η ext 1.05 ).
Figure 8. Color images of the nondimensional extinction cross sections η ext of 27-sphere clusters versus m p , i / m h and the ratio of average interparticle clearance to wavelength c . The results are averaged over 132 different orientations and 10 different configurations. The radius of the particle is a i = 40 nm, and the incident wavelength is λ = 4.0 μm. The complex refractive index of the host medium is m h = 1.5154 + i m h with (a) m h = 0.0, (b) m h = 0.0001, (c) m h = 0.001 and (d) m h = 0.01. The white regions represent the ISA being satisfied ( 0.95 η ext 1.05 ).
Photonics 10 00782 g008
Figure 9. The nondimensional extinction cross section η ext of 27-sphere cluster versus the ratio of average interparticle clearance to wavelength c with incident wavelength of λ = 4.0 μm, particle radius of a i = 40 nm, and particle refractive index of m p , i = 1.38, 2.1216 and 3.0308 ( m p , i / m h = 0.911, 1.4, and 2.0), respectively. The complex refractive index of the host medium is m h = 1.5154 + i m h with (a) m h = 0.0, (b) m h = 0.0001, (c) m h = 0.001 and (d) m h = 0.01. The blue region indicates that the ISA is satisfied for all the three m p , i / m h considered. The results are averaged over 132 different orientations and 10 different configurations.
Figure 9. The nondimensional extinction cross section η ext of 27-sphere cluster versus the ratio of average interparticle clearance to wavelength c with incident wavelength of λ = 4.0 μm, particle radius of a i = 40 nm, and particle refractive index of m p , i = 1.38, 2.1216 and 3.0308 ( m p , i / m h = 0.911, 1.4, and 2.0), respectively. The complex refractive index of the host medium is m h = 1.5154 + i m h with (a) m h = 0.0, (b) m h = 0.0001, (c) m h = 0.001 and (d) m h = 0.01. The blue region indicates that the ISA is satisfied for all the three m p , i / m h considered. The results are averaged over 132 different orientations and 10 different configurations.
Photonics 10 00782 g009
Figure 10. Color images of the nondimensional extinction cross sections η ext of 27-sphere clusters versus m p , i / m h and the ratio of average interparticle clearance to wavelength c . The results are averaged over 132 different orientations and 10 different configurations. The radius of the particle is a i = 80 nm, and the incident wavelength is λ = 4.0 μm. The complex refractive index of the host medium is m h = 1.5154 + i m h with (a) m h = 0.0, (b) m h = 0.0001, (c) m h = 0.001 and (d) m h = 0.01. The white regions indicate that the ISA is satisfied (0.95   η ext 1.05).
Figure 10. Color images of the nondimensional extinction cross sections η ext of 27-sphere clusters versus m p , i / m h and the ratio of average interparticle clearance to wavelength c . The results are averaged over 132 different orientations and 10 different configurations. The radius of the particle is a i = 80 nm, and the incident wavelength is λ = 4.0 μm. The complex refractive index of the host medium is m h = 1.5154 + i m h with (a) m h = 0.0, (b) m h = 0.0001, (c) m h = 0.001 and (d) m h = 0.01. The white regions indicate that the ISA is satisfied (0.95   η ext 1.05).
Photonics 10 00782 g010
Figure 11. As in Figure 10, but for a i = 160 nm. c is the ratio of average interparticle clearance to wavelength.
Figure 11. As in Figure 10, but for a i = 160 nm. c is the ratio of average interparticle clearance to wavelength.
Photonics 10 00782 g011
Figure 12. The nondimensional extinction cross sections η ext of 27-sphere SiO2 clusters embedded in PE and its nonabsorbing counterpart versus the ratio of average interparticle clearance to wavelength c . The results are averaged over 132 different orientations and 10 different configurations. The incident wavelengths are (ac) λ = 3.35 μm, (df) λ = 5.0 μm, and (gi) λ = 10.0 μm. The incident wavelengths are λ = 3.35, 5.0 and 10.0 μm, respectively. The particle radii are a i = 40, 80 and 160 nm.
Figure 12. The nondimensional extinction cross sections η ext of 27-sphere SiO2 clusters embedded in PE and its nonabsorbing counterpart versus the ratio of average interparticle clearance to wavelength c . The results are averaged over 132 different orientations and 10 different configurations. The incident wavelengths are (ac) λ = 3.35 μm, (df) λ = 5.0 μm, and (gi) λ = 10.0 μm. The incident wavelengths are λ = 3.35, 5.0 and 10.0 μm, respectively. The particle radii are a i = 40, 80 and 160 nm.
Photonics 10 00782 g012
Table 1. The extinction efficiency factors Q ext of the single particle with a i = 0.2, 1.0 and 2.0 μm, respectively. The complex refractive indices of the particle and the host medium are m p = 1.59 and m h = 1.0 + i m h with m h = 0.0, 0.001, 0.01 and 0.02, respectively.
Table 1. The extinction efficiency factors Q ext of the single particle with a i = 0.2, 1.0 and 2.0 μm, respectively. The complex refractive indices of the particle and the host medium are m p = 1.59 and m h = 1.0 + i m h with m h = 0.0, 0.001, 0.01 and 0.02, respectively.
m h = 0.0 m h = 0.001 m h = 0.01 m h = 0.02
a i = 0.2 μm0.000488550.00043460−0.000050653−0.00058887
a i = 1.0 μm0.297010.296190.288940.28123
a i = 2.0 μm2.42202.42422.44442.4675
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhai, J.; Zhang, S.; Liu, L. Extinction and Independent Scattering Criterion for Clusters of Spherical Particles Embedded in Absorbing Host Media. Photonics 2023, 10, 782. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics10070782

AMA Style

Zhai J, Zhang S, Liu L. Extinction and Independent Scattering Criterion for Clusters of Spherical Particles Embedded in Absorbing Host Media. Photonics. 2023; 10(7):782. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics10070782

Chicago/Turabian Style

Zhai, Jinan, Shangyu Zhang, and Linhua Liu. 2023. "Extinction and Independent Scattering Criterion for Clusters of Spherical Particles Embedded in Absorbing Host Media" Photonics 10, no. 7: 782. https://0-doi-org.brum.beds.ac.uk/10.3390/photonics10070782

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop