Scattering and Radiative Properties of Morphologically Complex Carbonaceous Aerosols: A Systematic Modeling Study

: This paper provides a thorough modeling-based overview of the scattering and radiative properties of a wide variety of morphologically complex carbonaceous aerosols. Using the numerically-exact superposition T -matrix method, we examine the absorption enhancement, absorption Ångström exponent (AAE), backscattering linear depolarization ratio (LDR), and scattering matrix elements of black-carbon aerosols with 11 different model morphologies ranging from bare soot to completely embedded soot–sulfate and soot–brown carbon mixtures. Our size-averaged results show that ﬂuffy soot particles absorb more light than compact bare-soot clusters. For the same amount of absorbing material, the absorption cross section of internally mixed soot can be more than twice that of bare soot. Absorption increases as soot accumulates more coating material and can become saturated. The absorption enhancement is affected by particle size, morphology, wavelength, and the amount of coating. We refute the conventional belief that all carbonaceous aerosols have AAEs close to 1.0. Although LDRs caused by bare soot and certain carbonaceous particles are rather weak, LDRs generated by other soot-containing aerosols can reproduce strong depolarization measured by Burton et al. for aged smoke. We demonstrate that multi-wavelength LDR measurements can be used to identify the presence of morphologically complex carbonaceous particles, although additional observations can be needed for full characterization. Our results show that optical constants of the host/coating material can signiﬁcantly inﬂuence the scattering and absorption properties of soot-containing aerosols to the extent of changing the sign of linear polarization. We conclude that for an accurate estimate of black-carbon radiative forcing, one must take into account the complex morphologies of carbonaceous aerosols in remote sensing studies as well as in atmospheric radiation computations.

Mishchenko et al. [84] have demonstrated that specific complex morphologies of soot-containing aerosols can reproduce the observed spectral dependence of the LDRs and thereby confirmed the potential of using depolarization measurements to identify and characterize smoke particles.
It has recently been suggested that making passive multiangle photopolarimetric measurements and especially combining active and passive remote sensing can improve the detection and accurate quantitative characterization of BC aerosols [87,88]. In this case the knowledge of the entire scattering matrix as a function of BC morphology becomes essential.
Consistent with the above rationale, the primary objective of this paper is to perform a thorough and systematic theoretical overview of the scattering and radiative properties of a wide variety of morphologically complex carbonaceous aerosols, including their spectral dependence. Our goal is to obtain a much improved quantitative understanding of how the absorption enhancement, AAE, LDR, and scattering matrix elements depend on BC-particle microphysics and quantify the likely range of their variability. To this end, we use the highly efficient and numerically-exact superposition T-matrix method (STMM), which is a first-principles scattering methodology rendering direct computer solutions of the macroscopic Maxwell equations [73][74][75][76]89]. Given the enormous diversity and variability of morphologies of carbonaceous aerosols, any single study cannot be exhaustive. Yet we believe that our systematic analysis is uniquely representative by virtue of featuring 11 soot-particle morphological models ranging from bare (both fluffy and compact) to partially embedded to fully-engulfed soot. We specifically analyze the factors that are potentially critical but have often been left out or downplayed in previous studies, i.e., the dependence of light scattering and absorption on (i) the distribution of soot inside a host, (ii) the fraction of soot enclosed inside the host, (iii) the volume fraction of the host/coating material, (iv) the coating thickness of each individual spherule, and (v) the range of monomer radii. We assume that soot and soot-containing aerosols obey a power law size distribution with the effective volume-equivalent radii consistent with observations. In addition, two different types of coating substance, viz., sulfate and BrC, are considered. The calculations are carried out at multiple wavelengths and for spectrally dependent refractive indices for both the BC and the host/coating material.

Model Aerosol-Particle Morphologies
It has been demonstrated that the typical overall morphology of a relatively fresh dry BC particle is well represented by a fractal cluster described by the following statistical scaling law [81,90]: where a is the monomer mean radius; k 0 is the fractal prefactor; D f is the fractal dimension; N S is the number of monomers in the cluster; and R g , called the radius of gyration, is a measure of the overall cluster radius. Aged soot-containing aerosols can exhibit varying degrees of deviation from this primordial morphology. Figure 1 depicts the 11 morphological models of soot and soot-containing aerosols considered in this study, ranging from bare soot to completely embedded soot mixtures. Models 1-3 represent bare soot particles. Model 1 is a freshly emitted fluffy soot aggregate with k 0 = 1.2 and D f = 1.8, Model 2 represents a somewhat aged soot aggregate with k 0 = 1.2 and D f = 2.6, and Model 3 shows a compact soot cluster whose overall shape is constrained by an imaginary spheroid with an aspect ratio of 2.
Specifically, we assumed that the distribution of soot and soot-containing particles shown in Figure  1 over sizes could be parameterized by the standard power-law [112] where r is a "characteristic" aerosol-particle size. The effective size eff r and the effective variance  [112]. We assume throughout the paper that the effective variance is fixed at 0.1, which implies that =1 In Models 4 and 5, the BC material is semi-externally mixed with a spherical host aerosol. Model 4 shows a sulfate or BrC host particle and an aged soot aggregate with k 0 = 1.2 and D f = 2.6 in physical contact, while Model 5 represents a host particle randomly "dusted" with small soot spherules [91].
Model 6 represents the scenario wherein a compact soot fractal is partially embedded in a sulfate or BrC host. Models 7-11 represent soot particles completely embedded in a spherical host or uniformly coated by a sulfate or BrC material. Model 7 is an internal mixture wherein all soot material is compacted into a single large spherical inclusion. The simplest yet highly idealistic version of this model is the conventional concentric core-shell configuration discussed above. In Model 8, the individual primary soot monomers are quasi-randomly distributed inside the spherical host. Model 9 represents the so-called closed-cell morphology resulting from the accumulation of a refractory material around individual soot monomers forming a compact cluster (with k 0 = 1.2 and D f = 2.6). Model 10 represents a spherical sulfate or BrC particle hosting a fully engulfed compact soot cluster. Model 11 is composed of two spherical sulfate or BrC particles in contact, each encapsulating a compact soot aggregate.
Note that for Models 1-10, the amount of soot material is kept fixed, while Model 11 is a semi-external mixture of two Model 10 particles and therefore has twice the amount of soot.
We assume that the effective soot monomer radius is equal to 0.02 µm [20,31,92]. Two different values for the number of monomers are adopted in this study: N S = 125 and 422. Therefore, the corresponding cumulative volume-equivalent radii of soot are equal to 0.1 µm (N S = 125) and 0.15 µm (N S = 422), respectively, consistent with field measurements. Indeed, according to Ueda et al. [30], the peak refractory BC mass-equivalent-sphere diameters were around 0.2 µm for all the samples they collected. In another study [93], the mode for the externally mixed soot (without coating or with minimal coating) size distribution was 0.15-0.2 µm in terms of the equivalent-sphere diameter. The values of N S ranging from 132 to 498 and the area-equivalent radii within 264 to 590 nm were observed for individual freshly emitted wildfire carbonaceous aerosols by China et al. [31].
Aged soot has often been found to be internally mixed with sulfate aerosols [93,94]. Sulfates represent the host or coating material most commonly used in modeling studies [19,44,45,47,78,84,86,95]. Soot can also be mixed with or embedded in weakly absorbing carbonaceous materials derived from combustion of biomass and biofuels and termed brown carbon, or BrC [69]. Therefore, Remote Sens. 2018, 10, 1634 6 of 26 both sulfates and BrC are assumed to form the host or the coating in this study. We calculated light scattering and absorption by carbonaceous aerosols at three wavelengths λ = 355, 532, and 1064 nm representing the spectral channels of the HSRL-2 [14]. The refractive indices of soot at the three wavelengths were estimated according to Equations (19a) and (19b) of Ref. [96]. The corresponding sulfate refractive indices were interpolated from the tabulated values at a 50% relative humidity in Ref. [97]. The BrC refractive index at 532 nm was taken from Ref. [98]. The spectrally dependent refractive indices for all the chemical components used in this study are listed in Table 1. The soot volume fraction values considered in this study are f = 7% and 15%. It should be emphasized that since we kept the soot amount constant, the respective volume-equivalent radii of the host were R = 0.2426 and 0.1882 µm when the soot volume fractions were 7% and 15%, while N S = 125. When N S was equal to 422, the corresponding R values were 0.3640 and 0.2823 µm. The selection of these parameters was consistent with field measurements. Specifically, it was reported by Vester et al. [99] that the average soot volume fraction for particles of 0.15-0.8 µm in diameter at an urban background site in Mainz, Germany fell between 7-27%. According to Adachi and Buseck [93], the median value of the soot volume fraction for particle aerodynamic diameters 0.05-0.3 µm in Mexico City was 15%. Ueda et al. [29] found that the median soot volume fraction values were less than 20% for 0.2-0.4 µm and less than 10% for 0.4-0.7 µm particles at Cape Hedo. In a more recent study [30], the soot volume fractions were found to be 21-50% and ≈10%, respectively, for particles with diameters ≤0.6 and >0.6 µm.
We used the algorithm described by Mackowski [100,101] to generate fractal-parameterequivalent soot clusters composed of 125 or 422 equal-size monomers. The fractal prefactor k 0 was fixed at 1.2. For freshly emitted fluffy soot represented by Model 1, we assumed the fractal dimension D f = 1.8, while for other models, when applicable, the D f value was set to be 2.6, thereby corresponding to a compact aged soot cluster.

Optical Modeling
For accurate calculations of the scattering and radiative properties of morphologically complex soot particles, we needed a first-principles methodology based on direct computer solutions of the macroscopic Maxwell equations. The STMM and the DDA [53,102,103] appeared to be the most frequently used techniques. However, other approaches, such as the finite difference time domain method [85,104,105], the pseudo-spectral time domain method [106,107], and the finite element-boundary integral technique [108,109], have also been applied (see, e.g., Ref. [83] and numerous references therein).
Having a reasonable combination of accuracy, efficiency, and versatility in mind, we used the STMM in this study. Importantly, a recent implementation of this technique [76] has made it applicable to arbitrary configurations of spherical domains wherein any of the spheres can be located at points that are either internal or external to the other spheres (i.e., such as the configurations exemplified by . Furthermore, the corresponding program can efficiently be run on distributed-memory computer clusters as well as on serial platforms. Overall, the STMM programs used enable fast and accurate calculations of ensemble-averaged absorption and scattering characteristics of a representative subset of morphologically complex carbonaceous aerosols. Consistent with Ref. [110], we assume throughout this paper that all nonspherical and/or internally heterogeneous aerosols are randomly oriented and form a statistically isotropic and mirror-symmetric ensemble. The key single-scattering characteristics we were looking for were the (ensemble-averaged) scattering, C sca , and absorption, C abs , cross sections per particle, as well as the elements of the normalized 4 × 4 Stokes scattering matrix. The latter has the following well-known block-diagonal structure [73,111]: where Θ ∈ [0 • , 180 • ] is the scattering angle. The (1, 1) element of the scattering matrix, a 1 (Θ), is the conventional phase function; it is normalized according to and describes the angular distribution of the scattered intensity in the case of unpolarized incident light. The conventional LDR is defined as while the degree of linear polarization of the scattered light for the case of unpolarized incident light is given by −b 1 (Θ)/a 1 (Θ).
To make our numerical results representative of naturally polydisperse particle groups and further suppress interference and resonance effects typical of electromagnetic scattering by monodisperse particles [73], we supplemented orientation averaging with size averaging. Specifically, we assumed that the distribution of soot and soot-containing particles shown in Figure 1 over sizes could be parameterized by the standard power-law [112] where r is a "characteristic" aerosol-particle size. The effective size r eff and the effective variance v eff of the size distribution are then defined according to and respectively, where G = r 2 r 1 drn(r)πr 2 (9) [112]. We assume throughout the paper that the effective variance is fixed at 0.1, which implies that r 1 = 0.546765 × r eff and r 2 = 1.65323 × r eff . The computation of the corresponding size averages is based on the following formulas [73]: and similarly for the other elements of the normalized scattering matrix, where r i and w i (with 1 ≤ i ≤ N G ) were the division points and weights of a Gaussian quadrature formula on the interval [r 1 , r 2 ]. The number of quadrature points was kept at 100 for all models. According to this approach, the spherule radii for the bare soot aerosols (Models 1, 2, and 3) varied in the range [0.546765 × 0.02 µm, 1.65323 × 0.02 µm] which corresponded to an effective monomer radius of 0.02 µm. In the case of composite particles, the sizes of all the components were scaled simultaneously so that the volume fraction of soot remained the same at all quadrature division points.
In addition to size averaging, one could also think of averaging over other "random" parameters, for example over random locations of the soot cluster embedded in a large spherical host in Model 10. We found, however, that the "randomized" scattering and radiative characteristics of this model were virtually indistinguishable from those computed for the quasi-central position of the cluster inside the host. This result may imply that it was the combination of the quasi-random morphology of the soot cluster and size averaging that served as the dominant "randomization" factor.

Numerical Results for Monodisperse Carbonaceous Aerosols
Before analyzing realistic size-averaged results, it was useful to examine the range of sensitivity of the radiative and scattering characteristics introduced above to microphysical parameters of the individual members of an ensemble. Indeed, this helped us to understand how the ensemble-averaged results could differ from those for the individual members of the ensemble and to what extent the latter could be representative of the former. Figure 2 illustrates the variation of the absorption cross section and the AAE calculated for different monodisperse soot-containing particles (partially) coated by sulfates. For Models 6, 7, and 9, the soot volume fraction was fixed at f = 7%. All soot aggregates consisted of N S = 125 monomers having the monodisperse radius a = 0.02 µm.
The top row of panels displays the absorption cross section and AAE as functions of the fraction of soot monomers encapsulated inside the sulfate host for Model 6 particles. It is seen that C abs nearly doubles as more and more soot spherules get engulfed by the host.
The second row of panels shows the sensitivity of the absorption cross section and AAE to the position of the spherical soot inclusion within the host (Model 7) quantified as the distance from the center. R denotes the radius of the host particle and r is the radius of the soot inclusion; the value of Rr is thus the farthest distance from the center that the soot inclusion can move without crossing the outer surface. Given the specified model parameters, the corresponding R and r values were equal to 0.2426 µm and 0.1 µm, respectively. C abs was maximal when the soot inclusion was at the center of the host, i.e., in the case of the concentric core-shell model frequently used in climate models. As the inclusion moved away from the center toward the outer boundary, absorption gradually decreased. This is consistent with our previous findings reported in Ref. [113].
The third row of panels shows that C abs did not necessarily increase monotonously with increasing coating thickness for Model 9 aerosols. Finally, for Model 10 particles (the bottom row of panels), C abs grew with the size of the sulfate host and could become saturated as the host reached about 95% of the total particle volume.
The AAE computed for the {532 nm, 1064 nm} pair of wavelengths (the right-hand column in Figure 2) also displays a wide range of variability, from about 0.7 to about 1.5. For the specific range of model parameters used, the AAE values for the core-shell Model 7 particles were the closest to unity, i.e., to the "standard" value adopted in the majority of previous studies [2,[59][60][61][62][63][64]. The corresponding AAEs increased from ≈0.93 to ≈0.98 as the soot inclusion moved away from the center toward the boundary of the sulfate host. However, the majority of AAE values in Figure 2 were not close to unity but rather were quite variable depending on the particle microphysical parameters. For example, the AAE curve for Model 9 particles shows values of about 1.05 for an aged compact soot aggregate with minimal coating, but as the spherules accumulated sulfate mantles and grew, so did the AAE until it reached the maximum value about 1.5, when the mantle thickness was about five times the monomer radius; then it started to decline.
Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 9 of model parameters used, the AAE values for the core-shell Model 7 particles were the closest to unity, i.e., to the "standard" value adopted in the majority of previous studies [2,[59][60][61][62][63][64]. The corresponding AAEs increased from ≈0.93 to ≈0.98 as the soot inclusion moved away from the center toward the boundary of the sulfate host. However, the majority of AAE values in Figure 2 were not close to unity but rather were quite variable depending on the particle microphysical parameters. For example, the AAE curve for Model 9 particles shows values of about 1.05 for an aged compact soot aggregate with minimal coating, but as the spherules accumulated sulfate mantles and grew, so did the AAE until it reached the maximum value about 1.5, when the mantle thickness was about five times the monomer radius; then it started to decline.     Figure 2 and depicts the LDR at the 355 nm (left-hand column) and 532 nm (right-hand column) wavelengths for different model parameters. It is clear that linear depolarization was very sensitive to particle size, morphology, sulfate amount, and heterogeneity. Although the outer shape of Model 7 and 10 particles was perfectly spherical, the corresponding LDR values deviated from zero to some extent owing to the existence of the (eccentric) soot impurity. For Model 6 aerosols, the LDRs reached the maximum value when about 35% of the soot material was engulfed by the sulfate host. The LDRs were very close to zero for Model 9 particles when there was minimal or no coating, but could reach nearly 40% when the coating thickness was about seven times the spherule radius (i.e., ≈0.14 µm). The strong linear depolarization induced by Model 9 aerosols (as well as by other large particles of certain morphologies) can potentially explain the high LDR values identified in recent lidar observations [14,79,80].
Remote Sens. 2018, 10, x FOR PEER REVIEW 10 of 10 was very sensitive to particle size, morphology, sulfate amount, and heterogeneity. Although the outer shape of Model 7 and 10 particles was perfectly spherical, the corresponding LDR values deviated from zero to some extent owing to the existence of the (eccentric) soot impurity. For Model 6 aerosols, the LDRs reached the maximum value when about 35% of the soot material was engulfed by the sulfate host. The LDRs were very close to zero for Model 9 particles when there was minimal or no coating, but could reach nearly 40% when the coating thickness was about seven times the spherule radius (i.e., ≈ 0.14 µm). The strong linear depolarization induced by Model 9 aerosols (as well as by other large particles of certain morphologies) can potentially explain the high LDR values identified in recent lidar observations [14,79,80].   Figure 4 shows the LDRs at 355 nm and the {532 nm, 1064 nm} AAEs as functions of soot monomer radius a for different Model 2, 6, 9, and 11 soot and soot-sulfate particles and a ranging from 0.01 µm to 0.034 µm. The number of soot spherules was N S = 125 for Models 2, 6, and 9 and N S = 250 for Model 11; the soot volume fraction for Models 6, 9, and 11 was fixed at f = 7%. Since the soot volume fraction was fixed, the overall particle sizes shrunk or expanded accordingly with the size of the soot monomers. For example, for Model 6 aerosols, the sulfate host radius changed from 0.2426 µm to 0.1213 µm as a decreased from 0.02 µm down to 0.01 µm, while for Model 9 aerosols, the radius of the soot cores was always kept at 0.4121 times that of the shells.
Remote Sens. 2018, 10, x FOR PEER REVIEW 11 of 11 volume fraction was fixed, the overall particle sizes shrunk or expanded accordingly with the size of the soot monomers. For example, for Model 6 aerosols, the sulfate host radius changed from 0.2426 µm to 0.1213 µm as a decreased from 0.02 µm down to 0.01 µm, while for Model 9 aerosols, the radius of the soot cores was always kept at 0.4121 times that of the shells. Both the LDRs and the AAEs were highly sensitive to the change of particle morphology and size. For the same amount of soot, depending on the specific particle morphology, the difference in LDR values could easily exceed a factor of ten. Overall, the LDR values tended to increase with increasing particle size. When the soot spherule radius a was about 0.01 µm, there was virtually no depolarization for Model 6 aerosols with roughly 50% of soot material engulfed by the host. The LDR values rose quickly to over 40% with four local maxima and minima as a increased from about 0.01 to over 0.03 µm. Interestingly, a similar pattern was observed for Model 11 aerosols, but with stronger oscillations. Such non-uniqueness, wherein different a values can induce the same LDR value if all the other parameters are kept fixed, can potentially cause an ill-posed retrieval problem in a remotesensing study based on lidar measurements alone. Therefore, other measurements besides the LDR ones may be needed to accurately identify aerosol microphysics.  Both the LDRs and the AAEs were highly sensitive to the change of particle morphology and size. For the same amount of soot, depending on the specific particle morphology, the difference in LDR values could easily exceed a factor of ten. Overall, the LDR values tended to increase with increasing particle size. When the soot spherule radius a was about 0.01 µm, there was virtually no depolarization for Model 6 aerosols with roughly 50% of soot material engulfed by the host. The LDR values rose quickly to over 40% with four local maxima and minima as a increased from about 0.01 to over 0.03 µm. Interestingly, a similar pattern was observed for Model 11 aerosols, but with stronger oscillations. Such non-uniqueness, wherein different a values can induce the same LDR value if all the other parameters are kept fixed, can potentially cause an ill-posed retrieval problem in a remote-sensing study based on lidar measurements alone. Therefore, other measurements besides the LDR ones may be needed to accurately identify aerosol microphysics.
The AAEs can either increase or decrease with increasing soot monomer radius, as the right-hand column of Figure 4 illustrates. One cannot help but notice that the AAE value was about 1.16 for the bare soot aggregate (Model 2, k 0 = 1.2, D f = 2.6, and N S = 125) even when the soot spherules were as small as ≈0.01 µm in radius and the cumulative soot volume-equivalent radius was r = 0.05 µm. The AAEs gradually decreased to the "conventional" value 1.0 when a = 0.022 µm, and then to ≈0.85 when a exceeded 0.03 µm. The deviation of the AAEs from 1.0 can be explained by our use of spectrally varying refractive indices and, more importantly, by the inability of the Rayleigh approximation to describe light scattering and absorption by morphologically complex soot aggregates, even though the individual soot spherules were much smaller than the incident wavelength. It is also evident from Figure 4 that once the soot aerosols became coated with refractive materials, the AAE values were rather different from those of bare soot. Coating could either increase or decrease the AAE values depending on the specific particle morphology. Figure 5 shows the size-averaged absorption cross sections at λ = 355 nm (top row), 532 nm (middle row), and 1064 nm (bottom row) for Model 1-10 aerosols. The left-hand and right-hand columns depict the results for aerosols containing N S = 125 and N S = 422 soot spherules, respectively, with an effective soot-monomer radius of a = 0.02 µm. We assumed that 50% of the soot spherules were engulfed by the sulfate host for Model 6 particles. In the case of Model 7, all soot spherules were "compacted" into a single centrally positioned spherical volume-equivalent soot core, thereby representing the conventional concentric core-shell model [39]. The blue and red curves show C abs for the soot volume fractions f = 7% and 15%, respectively. The effective cumulative volume-equivalent soot radii corresponding to a = 0.02 µm were equal to 0.1 and 0.15 µm. For the bare soot aggregates (Models 1-3), the results shown by the blue and red curves were identical.

Ensemble-Averaged Numerical Results
It is clear from Figure 5 that freshly emitted fluffy soot particles (Model 1) tended to absorb more light than the other types of bare soot with a more compact structure (Models 2 and 3). As pointed out by Liu et al. [35], this may be the consequence of the shielding effect of the outer layer of monomers that blocks light from penetrating deeply inside a compact aggregate. The semi-external mixture represented by Model 5 was a stronger absorber than that represented by Model 4 owing to the more direct exposure of the soot spherules to the incident light. The role of the strongest light absorbers, however, belonged to soot-sulfate internal mixtures. Figure 5 shows that for the same amount of absorbing material, C abs of the internal mixtures could be more than twice that of bare soot depending on how the soot was distributed throughout the host particle. At λ = 355 nm, the size-averaged absorption cross section was about 0.043 µm 2 for Model 3 particles with N S = 125 and a = 0.02 µm, but jumped to about 0.091 µm 2 for Model 8 aerosols with a 7% soot volume fraction. For N S = 422, C abs increased from about 0.113 µm 2 to about 0.263 µm 2 when the particle morphology changed from Model 3 to Model 9.
The absorption enhancement was influenced by particle size (left-hand vs right-hand columns), morphology (different models), wavelength (different rows), and amount of coating (blue vs red curves). No single morphology served as the universally strongest absorber: depending on many factors accommodated in Figure 5, it could be Model 7, 8, 9, or 10. It was always true, however, that for the same amount of absorbing material, soot-containing aerosols with a thicker coating absorbed more light than their counterparts with less coating, unless the amount of coating was so large that absorption got saturated. Remote Sens. 2018, 10, x FOR PEER REVIEW 13 of 13 It is clear from Figure 5 that freshly emitted fluffy soot particles (Model 1) tended to absorb more light than the other types of bare soot with a more compact structure (Models 2 and 3). As pointed out by Liu et al. [35], this may be the consequence of the shielding effect of the outer layer of monomers that blocks light from penetrating deeply inside a compact aggregate. The semi-external mixture represented by Model 5 was a stronger absorber than that represented by Model 4 owing to the more direct exposure of the soot spherules to the incident light. The role of the strongest light absorbers, however, belonged to soot-sulfate internal mixtures. Figure 5 shows that for the same amount of absorbing material, Cabs of the internal mixtures could be more than twice that of bare soot depending on how the soot was distributed throughout the host particle. At λ = 355 nm, the sizeaveraged absorption cross section was about 0.043 µm 2 for Model 3 particles with NS = 125 and a = 0.02 µm, but jumped to about 0.091 µm 2 for Model 8 aerosols with a 7% soot volume fraction. For NS = 422, Cabs increased from about 0.113 µm 2 to about 0.263 µm 2 when the particle morphology changed from Model 3 to Model 9.
The absorption enhancement was influenced by particle size (left-hand vs right-hand columns), morphology (different models), wavelength (different rows), and amount of coating (blue vs red curves). No single morphology served as the universally strongest absorber: depending on many factors accommodated in Figure 5, it could be Model 7, 8, 9, or 10. It was always true, however, that for the same amount of absorbing material, soot-containing aerosols with a thicker coating absorbed   Figure 5 and illustrates the variability of the AAEs, calculated for three different pairs of wavelengths, with varying aerosol morphology. Again, the effective soot monomer radius a was fixed at 0.02 µm. Our numerically-exact results show that the AAE was quite sensitive to the choice of pair of wavelengths as well as to the particle morphology and size. It could also be sensitive to the amount of coating material. The fresh-soot AAEs (Model 1) were always larger than the AAEs of the other two types of bare soot (Models 2 and 3). By comparing the results on the left with those on the right, we concluded that the AAE could significantly decrease with increasing aerosol-particle size and soot volume fraction. This was consistent with the finding by Liu et al. [71]. When we compared the AAEs for "soot-aggregate-equivalent" Model 2 and 10 particles, we could not conclude definitively that coating increased (or decreased) the AAE, as both scenarios were possible. Among the morphological types considered in this study, the concentric core-shell model produced the smallest AAE values: it could become as low as 0.1 when the AAE was calculated for the {355 nm, 532 nm} pair of wavelengths. Although the size-averaged AAEs were close to 1.0 for some aerosol particles with N S = 125, the overall range of AAEs was quite wide: from as low as 0.1 to about 1.7. We thus cannot help but conclude that using a flat value of 1.0 is likely to introduce significant errors in soot-aerosol absorption analyses.  The concentric core-shell Model 7 was spherically symmetric and hence caused zero depolarization. Figure 7 depicts the size-averaged LDRs calculated for seven other models containing N S = 125 (left-hand panels) or 422 (right-hand panels) BC spherules with an effective radius a = 0.02 µm. Model 11 was a semi-external mixture of two Model 10 particles and therefore contained twice as much soot. In most cases, the solid lines lay above the corresponding dotted ones, thereby indicating that particles with more sulfate material generated larger LDR values than their BC-equivalent counterparts with less sulfate material. The only exception was Model 10, which appears to be inconsistent with the bottom panels of Figure 3. It should be kept in mind, however, that Figure 3 shows the results computed for monodisperse aerosols and may not be representative of size-averaged LDRs.
As expected, the LDRs caused by Model 4 and 5 particles are rather insignificant and may often be within the instrumental noise level. As such, they hardly contain much useful information. However, we did see larger LDRs for other morphologies such as Models 6 and 11. Figure 7 confirms the potential of exploiting polarization backscattering lidars for optical characterization of morphologically complex particulates, especially when depolarization data are available at multiple wavelengths. For example, a spherical sulfate host with a partially imbedded (in this case 50%) soot aggregate (Model 6) with N S = 422, f = 7% and a semi-external mixture of two heterogeneous sulfate particles each encapsulating a soot aggregate (Model 11) with N S = 125, f = 7% both produced LDR values at λ = 355, 532, and 1064 nm within the range observed for an aged smoke plume by Burton et al. [14] (δ L = 0.203 ± 0.053 at 355 nm, 0.093 ± 0.026 at 532 nm, and 0.018 ± 0.009 at 1064 nm). It is worth mentioning that Model 9 particles were also able to generate spectral LDRs consistent with the lidar measurements by Burton et al. provided that the sulfate coating of each individual soot spherule was of the right thickness [84]. As expected, the LDRs caused by Model 4 and 5 particles are rather insignificant and may often be within the instrumental noise level. As such, they hardly contain much useful information. However, we did see larger LDRs for other morphologies such as Models 6 and 11. Figure 7 confirms the potential of exploiting polarization backscattering lidars for optical characterization of morphologically complex particulates, especially when depolarization data are available at multiple wavelengths. For example, a spherical sulfate host with a partially imbedded (in this case 50%) soot aggregate (Model 6) with NS = 422, f = 7% and a semi-external mixture of two heterogeneous sulfate Active lidar measurements (especially multi-wavelength lidar systems [114]) coupled with advanced passive photopolarimetric observations may provide the capability to identify and Remote Sens. 2018, 10, 1634 16 of 26 characterize aerosols with nonspherical geometries such as soot and dust particles [87,88,115]. The fidelity of such retrievals obviously depends on the adequacy of modeling the Stokes scattering matrix. As an example, Figure 8 illustrates the differences in the size-averaged phase function a 1 (Θ) and linear polarization −b 1 (Θ)/a 1 (Θ) between different soot and soot-sulfate aerosol models. The number of soot spherules was N S = 125, and their effective radius a was fixed at 0.02 µm. Again, in the case of the concentric core-shell Model 7, all the soot spherules were consolidated into one single volume-equivalent spherical soot core. The solid and dotted curves depict the results calculated for soot volume fractions of f = 7% and 15%, respectively. As before, we assumed that 50% of soot spherules were engulfed by the host particles in the case of Model 6. measurements by Burton et al. provided that the sulfate coating of each individual soot spherule was of the right thickness [84].
Active lidar measurements (especially multi-wavelength lidar systems [114]) coupled with advanced passive photopolarimetric observations may provide the capability to identify and characterize aerosols with nonspherical geometries such as soot and dust particles [87,88,115]. The fidelity of such retrievals obviously depends on the adequacy of modeling the Stokes scattering matrix. As an example, Figure 8 illustrates the differences in the size-averaged phase function 1 ( ) a Θ and linear polarization -1 1 ( ) ( ) b Θ a Θ between different soot and soot-sulfate aerosol models. The number of soot spherules was NS = 125, and their effective radius a was fixed at 0.02 µm. Again, in the case of the concentric core-shell Model 7, all the soot spherules were consolidated into one single volume-equivalent spherical soot core. The solid and dotted curves depict the results calculated for soot volume fractions of f = 7% and 15%, respectively. As before, we assumed that 50% of soot spherules were engulfed by the host particles in the case of Model 6.   Figure 9, reveal a considerable degree of dependence of both a 1 (Θ) and −b 1 (Θ)/a 1 (Θ) on aerosol-particle morphology. For bare soot aggregates represented by Models 1-3, we could not help but notice that the degree of linear polarization for unpolarized incident light is essentially Rayleigh-like, with a nearly 100% maximum at Θ ≈ 90 • , thereby implying the dominant role of single scattering by the subwavelength soot monomers. At the same time, the corresponding phase functions had a pronounced forward-scattering peak typical of wavelength-sized particles. This remarkable scattering behavior of large aggregates composed of Rayleigh-sized monomers was first identified by West [116,117]. As the soot aggregates acquired more sulfate material or got attached to a sulfate host, the −b 1 (Θ)/a 1 (Θ) became significantly reduced, the position of maximal polarization shifted toward larger scattering angles around 130 • , and a negative-polarization branch could develop at backscattering angles. This behavior became even more prominent with increasing N S , as shown in Figure 9.
Remote Sens. 2018, 10, x FOR PEER REVIEW 18 of 18 Figure 8 and its companion, Figure 9, reveal a considerable degree of dependence of both 1 ( ) a Θ and -1 1 ( ) ( ) b Θ a Θ on aerosol-particle morphology. For bare soot aggregates represented by Models 1-3, we could not help but notice that the degree of linear polarization for unpolarized incident light is essentially Rayleigh-like, with a nearly 100% maximum at Θ ≈ 90°, thereby implying the dominant role of single scattering by the subwavelength soot monomers. At the same time, the corresponding phase functions had a pronounced forward-scattering peak typical of wavelength-sized particles. This remarkable scattering behavior of large aggregates composed of Rayleigh-sized monomers was first identified by West [116,117]. As the soot aggregates acquired more sulfate material or got attached to a sulfate host, the -1 1 ( ) ( ) b Θ a Θ became significantly reduced, the position of maximal polarization shifted toward larger scattering angles around 130°, and a negative-polarization branch could develop at backscattering angles. This behavior became even more prominent with increasing NS, as shown in Figure 9. Our computations of light scattering and absorption by soot particles hosted or coated by BrC material were less comprehensive and will be expanded once consistent and accurate BrC refractive indices at λ = 355 and 1064 nm become available. Based on our limited data shown in Figure 10, it is clear that even with the same size and morphology, the physical properties of the BrC material can fundamentally change the overall optical properties of soot-containing aerosols. Indeed, comparison of the dotted and the corresponding solid curves reveals significant differences in the phase function and the degree of linear polarization for the soot-equivalent aerosols, to the extent of the polarization switching sign. Our computations of light scattering and absorption by soot particles hosted or coated by BrC material were less comprehensive and will be expanded once consistent and accurate BrC refractive indices at λ = 355 and 1064 nm become available. Based on our limited data shown in Figure 10, it is clear that even with the same size and morphology, the physical properties of the BrC material can fundamentally change the overall optical properties of soot-containing aerosols. Indeed, comparison of the dotted and the corresponding solid curves reveals significant differences in the phase function and the degree of linear polarization for the soot-equivalent aerosols, to the extent of the polarization switching sign.
Remote Sens. 2018, 10, x FOR PEER REVIEW 19 of 19 Figure 10. As in Figure 8, except that now the solid and dotted curves represent the results calculated for soot-containing aerosols with brown-carbon (BrC) and sulfate (Sul) host/coating, respectively. Figure 11 shows the corresponding LDR differences at λ = 532 nm. We can see that the LDRs for Models 5, 8, and 10 were rather small and showed little dependence on the host/coating material. However, the LDR for Model 11 particles could be significantly altered by the nature of the host substance.  Figure 8, except that now the solid and dotted curves represent the results calculated for soot-containing aerosols with brown-carbon (BrC) and sulfate (Sul) host/coating, respectively. Figure 11 shows the corresponding LDR differences at λ = 532 nm. We can see that the LDRs for Models 5, 8, and 10 were rather small and showed little dependence on the host/coating material. However, the LDR for Model 11 particles could be significantly altered by the nature of the host substance. Remote Sens. 2018, 10, x FOR PEER REVIEW 20 of 20 Figure 11. LDR differences at λ = 532 nm between soot-containing aerosols hosted/coated by sulfate (blue lines) and brown-carbon (red lines) materials. The soot volume fraction is fixed at f = 7%. The number of soot spherules is NS = 125 (left-hand panel) and 422 (right-hand panel).

Conclusions and Future Work
We have used the numerically-exact and highly efficient superposition T-matrix method to investigate the absorption enhancement, the absorption Ångström exponent, the backscattering linear depolarization ratio, and essential scattering matrix elements of carbonaceous aerosols represented by 11 distinct morphologies ranging from bare soot to completely embedded soot-sulfate and soot-BrC internal mixtures.
Our sensitivity analysis for monodisperse aerosols shows that for the commonly used core-shell model (Model 7), absorption was maximal when the soot inclusion was centrally positioned and steadily weakened as the soot inclusion moved closer to the boundary of the host. We demonstrated numerically that for Model 6 (a soot aggregate partially imbedded in the spherical sulfate host), absorption nearly doubled as the percentage of the monomers engulfed by the host increased from zero to 100%. For the Model 10 morphology (a spherical sulfate particle hosting a completely imbedded compact soot cluster), absorption increased with the volume fraction of the sulfate material and could get saturated as this fraction reached ≈ 95%.
Our size-averaged results demonstrated once again that freshly emitted fluffy soot particles (Model 1) tended to absorb more light compared with the other types of bare soot aerosols having more compact morphologies (Models 2 and 3). Model 3 soot clusters, which mimic the cocoon shape observed by Ueda et al. [30], were the least efficient absorbers. For the same amount of absorbing material, the absorption cross section of internally mixed soot could be more than twice that of bare soot. In general, the absorption enhancement was influenced by particle size, morphology, wavelength, and amount of sulfate or BrC material. As a result, no single morphological model served as a universally strongest absorber. It is worth emphasizing that although many previous studies have indicated that the BC absorption enhancement predicted by the concentric core-shell morphology is biased high [19,[49][50][51][52], which may indeed be the case when the calculations and/or measurements are performed at a longer wavelength, such as 1064 nm, the absorption by the concentric core-shell geometry at λ = 355 nm was actually weaker than that by the other "equivalent" soot-sulfate internal mixtures.
We believe that the absorption enhancement is also strongly affected by the choice of constituent refractive indices. In this study, we used a set of pre-defined spectrally varying refractive indices. In the future, it may be instructive to consider refractive indices that represent the higher and lower bounds of BC absorption.
Our analysis of the AAE largely validated the conclusions drawn by Liu et al. [71] and extended their study to a great extent. In general, the AAE was quite sensitive to the choice of the wavelength pair, particle morphology, and size, and the constituent refractive indices. While the conventional assumption is that the AAE is ≈1.0 for isolated BC particles, our computations of the size-averaged Figure 11. LDR differences at λ = 532 nm between soot-containing aerosols hosted/coated by sulfate (blue lines) and brown-carbon (red lines) materials. The soot volume fraction is fixed at f = 7%. The number of soot spherules is N S = 125 (left-hand panel) and 422 (right-hand panel).

Conclusions and Future Work
We have used the numerically-exact and highly efficient superposition T-matrix method to investigate the absorption enhancement, the absorption Ångström exponent, the backscattering linear depolarization ratio, and essential scattering matrix elements of carbonaceous aerosols represented by 11 distinct morphologies ranging from bare soot to completely embedded soot-sulfate and soot-BrC internal mixtures.
Our sensitivity analysis for monodisperse aerosols shows that for the commonly used core-shell model (Model 7), absorption was maximal when the soot inclusion was centrally positioned and steadily weakened as the soot inclusion moved closer to the boundary of the host. We demonstrated numerically that for Model 6 (a soot aggregate partially imbedded in the spherical sulfate host), absorption nearly doubled as the percentage of the monomers engulfed by the host increased from zero to 100%. For the Model 10 morphology (a spherical sulfate particle hosting a completely imbedded compact soot cluster), absorption increased with the volume fraction of the sulfate material and could get saturated as this fraction reached ≈95%.
Our size-averaged results demonstrated once again that freshly emitted fluffy soot particles (Model 1) tended to absorb more light compared with the other types of bare soot aerosols having more compact morphologies (Models 2 and 3). Model 3 soot clusters, which mimic the cocoon shape observed by Ueda et al. [30], were the least efficient absorbers. For the same amount of absorbing material, the absorption cross section of internally mixed soot could be more than twice that of bare soot. In general, the absorption enhancement was influenced by particle size, morphology, wavelength, and amount of sulfate or BrC material. As a result, no single morphological model served as a universally strongest absorber. It is worth emphasizing that although many previous studies have indicated that the BC absorption enhancement predicted by the concentric core-shell morphology is biased high [19,[49][50][51][52], which may indeed be the case when the calculations and/or measurements are performed at a longer wavelength, such as 1064 nm, the absorption by the concentric core-shell geometry at λ = 355 nm was actually weaker than that by the other "equivalent" soot-sulfate internal mixtures.
We believe that the absorption enhancement is also strongly affected by the choice of constituent refractive indices. In this study, we used a set of pre-defined spectrally varying refractive indices. In the future, it may be instructive to consider refractive indices that represent the higher and lower bounds of BC absorption.
Our analysis of the AAE largely validated the conclusions drawn by Liu et al. [71] and extended their study to a great extent. In general, the AAE was quite sensitive to the choice of the wavelength pair, particle morphology, and size, and the constituent refractive indices. While the conventional assumption is that the AAE is ≈1.0 for isolated BC particles, our computations of the size-averaged AAE for bare-soot aerosols exhibited values ranging from ≈0.68 to ≈1.26. Furthermore, when soot was semi-externally or internally mixed with sulfates, the range of AAE values extended from ≈0.1 to ≈1.7. Thus, using the flat value ≈1.0 can be thoroughly inappropriate.
Our computations confirmed and further illustrated that size-averaged LDRs can be very sensitive to particle size, shape, and morphology. The LDRs produced by bare-soot aerosols and some types of internal soot mixtures were rather small and hence hardly useful. However, Model 6, 9, and 11 aerosols containing large amounts of refractory materials could cause much larger LDRs and could even reproduce exceptionally high, spectrally dependent values such as those measured by Burton et al. [14] for a plume of aged smoke. Our results demonstrate that spectral LDR measurements can be used to identify the presence of morphologically complex smoke particles, even though additional observations (e.g., with a passive polarimeter [87,88,[118][119][120][121] or a bistatic polarization lidar [122][123][124][125]) may be required to narrow down the plausible ranges of particle morphology (including size) and composition.
The scattering phase function and the degree of linear polarization also varied greatly among different carbonaceous aerosol models having the same amount of absorbing material. While bare-soot aerosols exhibited classical bell-shaped polarization curves typical of Rayleigh scattering, mixing them semi-externally or internally with sulfate or BrC materials could modify the scattering-angle dependence of polarization quite dramatically. This result emphasized the importance of using appropriate aerosol models in analyses of remote-sensing observations and in situ measurements.
In summary, the results of our extensive computations based on numerically-exact solutions of the Maxwell equations further demonstrate that the scattering and radiative properties of particulate soot can vary significantly from the moment it is injected into the atmosphere as it goes through different processes of aging all the way until it settles. Our study confirms the importance of explicit morphological and optical modeling of soot and soot-containing aerosols in atmospheric-radiation computations as well as in remote-sensing studies.
Being as extensive and systematic as it is, our analysis is still not exhaustive. For example, it needs to be supplemented by computations of light scattering and absorption at λ = 355 nm and 1064 nm for BC-BrC semi-external and internal mixtures. It is also important to study the effects of other hosting/coating materials with different refractive indices. Furthermore, significant effort is required in order to incorporate complex morphological models of BC aerosols in retrieval algorithms and climate-related computations of atmospheric radiation. Lastly, extensive comparisons of modeling results with experiments must be a priority.