the Creative Commons Attribution 4.0 License.

Special issue: SKYNET – the international network for aerosol, clouds,...

**Research article**
07 Jan 2019

**Research article** | 07 Jan 2019

# A shape model of internally mixed soot particles derived from artificial surface tension

Hiroshi Ishimoto, Rei Kudo, and Kouji Adachi

**Hiroshi Ishimoto et al.**Hiroshi Ishimoto, Rei Kudo, and Kouji Adachi

- Meteorological Research Institute, Tsukuba, 305-0052, Japan

- Meteorological Research Institute, Tsukuba, 305-0052, Japan

**Correspondence**: Hiroshi Ishimoto (hiroishi@mri-jma.go.jp)

**Correspondence**: Hiroshi Ishimoto (hiroishi@mri-jma.go.jp)

Received: 25 Jul 2018 – Discussion started: 20 Aug 2018 – Revised: 12 Nov 2018 – Accepted: 11 Dec 2018 – Published: 07 Jan 2019

To retrieve the physical properties of aerosols from multi-channel ground-based and satellite measurements, we developed a shape model of coated soot particles and created a dataset of their optical properties. Bare soot particles were assumed to have an aggregate shape, and two types of aggregates with different size–shape dependences were modeled using a polyhedral Voronoi structure. To simulate the detailed shape properties of mixtures of soot aggregates and adhered water-soluble substances, we propose a simple model of surface tension derived from the artificial surface potential. The light-scattering properties of the modeled particles with different volume fractions of water-soluble material were calculated using the finite-difference time-domain method and discrete-dipole approximation. The results of the single-scattering albedo and asymmetry factors were compared to those of conventional internally mixed spheres (i.e., effective medium spheres based on the Maxwell-Garnett approximation and simple core-shell spheres). In addition, the lidar backscattering properties (i.e., lidar ratios and linear depolarization ratios) of the modeled soot particles were investigated. For internally mixed soot particles, the lidar backscattering properties were sensitive to the shape of the soot particles and the volume mixing ratio of the assumed water-soluble components. However, the average optical properties of biomass smoke, which have been reported from in situ field and laboratory measurements, were difficult to explain based on the individually modeled particle. Nonetheless, our shape model and its calculated optical properties are expected to be useful as an alternative model for biomass smoke particles in advanced remote sensing via multi-channel radiometer and lidar measurements.

- Article
(7714 KB) -
Supplement
(54 KB) - BibTeX
- EndNote

During the atmosphere aging process of emitted combustion products, soot particles tend to become hydrophilic and form mixtures with weakly light-absorbing materials (Mikhailov et al., 2006; Adachi et al., 2007; Moteki and Kondo, 2007; Shiraiwa et al., 2007, 2010; Adachi and Buseck, 2008). Because of significant enhancements in light absorption and scattering, it has been suggested that soot particles in a mixing state are the second most important contributor to global warming after carbon dioxide (Jacobson, 2001; Ramanathan et al., 2008). For climate monitoring and numerical prediction via atmospheric data assimilation, precise estimation of the amount of mixed soot particles and the fraction of soot (black carbon) from satellite- and ground-based remote-sensing measurements is important (Kahnert et al., 2013). Thus, understanding the optical properties of internally mixed soot particles is essential to improve the retrieval accuracy of atmospheric soot particles (Hara et al., 2018). The light-scattering properties of internally mixed particles depend strongly on the complex refractive index of each mixing component. Furthermore, the shape of the incorporated soot particles and overall particle shape in the mixing state significantly alter some of the scattering properties. In particular, particle shape is important for the interpretation of lidar backscattering measurements. Many shape models have been proposed for internally mixed soot particles and their light-scattering properties using the discrete-dipole approximation (DDA) and T-matrix methods (Adachi et al., 2010; Scarnato et al., 2013; Cheng et al., 2014; Dong et al., 2015; Liu et al., 2016; Mishchenko et al., 2016; Moteki, 2016; Wu et al., 2016, 2017; Kahnert, 2017; Zhang et al., 2017; Luo et al., 2018). However, the relationship between mixing state or morphology and light-scattering properties is not well defined.

As an alternative model for soot particles in a mixing state, we developed a new shape model of internally mixed particles. We assumed that bare soot particles with fractal-like shapes were mixed with water-soluble (WS) components. Moreover, we considered hydrophilic soot particles with high wettability due to atmospheric aging. The shapes of the mixed particles were determined by applying artificial potential calculations of the surface tension of the WS components. The numerical results of the light-scattering properties of the modeled particles at visible and near-infrared wavelengths are discussed.

## 2.1 Soot model

Bare soot particles are commonly described as fractal aggregates formed from
primary particles (or monomers) with a degree of overlapping and necking
between neighboring primary particles (Yon et al., 2015; Okyay et al., 2016).
Primary particles have a diameter of 20–50 nm (Bond and Bergstrom, 2006),
and the fractal dimension of the aggregates depends on emissions' conditions
and atmospheric aging. For example, newly generated soot aggregates often
form lace-like structures with relatively small fractal dimensions, whereas
aged soot aggregates tend to be compact and to be characterized by large
fractal dimensions (Mikhailov et al., 2001, 2006; Zhang et al., 2016). In
this study, we modeled soot aggregates using spatial Poisson–Voronoi
tessellation. The basic methods to make the aggregate model are the same as
those of dust particles and ice aggregates described in our previous works
(Ishimoto et al., 2010, 2012a, b; Baran et al., 2018). A spatial
Poisson–Voronoi tessellation was produced from randomly distributed
nucleation points in the numerical field, and polyhedral cells overlapping
with the assumed fractal frame were selected from the tessellation (Ishimoto
et al., 2012a). We regarded the cells (i.e., Voronoi cells) as the primary
soot aggregate particles to mimic overlapping and necking between neighboring
primary particles. To ensure that the size of the Voronoi cells was
relatively uniform, we applied a Matérn hard-core point field (Ohser and
Mücklich, 2000) to the spatial distribution of the nucleation points
(Ishimoto et al., 2012b). Aggregate particles of different sizes were
produced in the same manner but by changing the relative size of the fractal
frame within the same tessellation. Figure 1 shows two sets of aggregate
particles (Types A and B) created from fractal frames of different shapes. To
calculate the light-scattering properties, Type A aggregates of 10 sizes and
Type B aggregates of 13 sizes (where nos. 1, 2 and 3 were used for both
sets) were prepared. The mean radius of the primary particle *a*=20 nm was
used as a typical value (Wu et al., 2015, 2017; Mishchenko et al., 2016; Luo
et al., 2018), and the total size of each aggregate was corrected by
adjusting the average cell volume of the aggregate as ${V}_{\mathrm{c},\mathrm{agg}}=\mathrm{4}/\mathrm{3}\mathit{\pi}{a}^{\mathrm{3}}$. Although the aggregates shown in Fig. 1 were created using a
box-counting approach to maintain a fractal relationship between size and
shape, this fractal relationship differed from that of the commonly used
fractal dimension *D*_{f} used to describe individual soot particles.
The fractal dimension *D*_{f} is defined from the number of monomers
*N*, fractal prefactor *k*_{0} and gyration radius *R*_{g} as follows
(Adachi et al., 2007):

For numerical simulation of the light-scattering properties of aged soot
particles, various *D*_{f} values for a typical prefactor of
*k*_{0}=1.2 have been proposed in the literature, such as
*D*_{f}= 2.0–2.5 (Nyeki and Colbeck, 1995), 1.9–2.6 (Adachi et al.,
2007), 2.5 (He et al., 2015), 2.6 (Mishchenko et al., 2016) and 2.5–3.0
(Zhang et al., 2017). The corresponding *D*_{f} values for our modeled
soot particles derived from the calculation of *R*_{g} are plotted in
Fig. 2. Applying the fractal prefactor *k*_{0}=1.2 resulted in aggregate
particles in the range of *D*_{f}= 1.9–2.5 (Type A) and 2.5–3.0
(Type B) given a normalized gyration radius of $\mathrm{ln}\left({R}_{\mathrm{g}}/a\right)\le \mathrm{3.6}$.

## 2.2 Artificial surface tension of mixed soot and water-soluble components

According to microscopic images of internally mixed soot particles, soot particles are often entirely encapsulated in a spherical shell and completely covered by WS components (Reid and Hobbs, 1998; Reid et al., 2005a). Thus, WS components behave like a liquid on the particle surface, and the surface tension of WS components is important to describing the overall shapes of mixed particles. Surface tension is the result of intermolecular forces at the microscopic scale, described as the Lennard-Jones potential (Becker et al., 2014). Although some approaches to simulate surface tension in a discrete particle system, such as smoothed particle hydrodynamics, have been proposed (Tartakovsky and Meakin, 2005; Leinonen and von Lerber, 2018), a general approach to reproduce the variety of fluid effects has not yet been developed (Akinci et al., 2013). In this study, we examined a virtual and simple potential field on the surface of the modeled particle to simulate the morphological effects of surface tension. We assumed that the dynamic behavior of liquid could be simulated by the movement of liquid elements from locations of high potential to locations of low potential. From a simple molecular-scale explanation of surface tension, the potential at a surface point becomes smaller as the number of surrounding molecules increases. The modeled particles were projected in a three-dimensional Cartesian grid space to define surface points for the potential calculations. The artificial surface potential was defined as an analog of the microscopic surface potential.

where *U*_{j} is the artificial potential at surface point *j*, and *U*_{j}
is defined as the negative value of the total number of grid points of
material with distance *r*≤*d*. We used this simple model for the surface
potential calculations to focus on the shape of the mixture in discrete grid
space. An equilibrium mixing state was simulated assuming that the applied
WS components preferentially accumulate at grid points of lower potential.
The shape of internally mixed particles under an arbitrary volume mixing
ratio was determined based on iterative calculations of the surface
potential and by adding WS elements. In the artificial potential
calculations in Eq. (2), the setting of the length *d* is important. We applied two
steps for the potential calculations and adhesion of WS components for each
iterative calculation.

where *l* is the grid length of the space, *d*_{1} and
*d*_{2} are the lengths *d* for first and second steps and
*N*_{add} is the number of grid points for WS adhesion, which are
chosen from the total surface points *N*_{s} (note that a surface
point is defined as an empty grid and neighbor of an occupied grid). Small
*d*_{1} values resulted in a coated particle with a thin layer at a
local scale. By contrast, WS components tended to accumulate in the same
region on particles when a large *d*_{1} was applied. We assumed the
parameter *d*_{1} to be the minimum scale to derive isotropic potential in
the discrete Cartesian grid space and determined the value
*d*_{1}=3 *l* from our results for shapes of coated particles in
preliminary calculations. The second step in Eq. (3) is important,
particularly when WS components cover the entire soot aggregate, and the
value *d*_{2} is determined to ensure that the overall shape of the
mixed particle is spherical. The correction length *d*_{cor} is the
minimum length at which the curvature of the sphere with a radius *R* can be
discriminated from the calculated artificial potential in the grid field, and
*d*_{cor} is estimated from the following relationships.

where *N*_{V} is the total number of grid points occupied by the material,
and *γ* is the effective skin depth (*γ*=1.51*l*). Figure 3 presents a schematic diagram describing *d*_{cor} in a two-dimensional
case.

The artificial potential was calculated for each iterative step, and
*N*_{add} of WS elements was added, starting with the grid points of
lowest potential. The internally mixed soot model (i.e., the particle shape
modeled based on artificial surface tension, AST hereafter) was created for
different values of the volume ratio
${V}_{\mathrm{r}}={V}_{\mathrm{ws}}/{V}_{\mathrm{soot}}$, where *V*_{ws} and
*V*_{soot} are the volumes of the WS and soot materials,
respectively. Figure 4 shows the results of several mixed soot models for
${V}_{\mathrm{r}}\sim \mathrm{0},\mathrm{2},\mathrm{5},\mathrm{10},\mathrm{20}$. The total number of iterations was
approximately 1000 (2000) to create *V*_{r}∼10 (*V*_{r}∼20)
particles. For simplicity, we neglected the difference in materials (i.e.,
soot or WS) in the potential calculations of Eq. (2), which implicitly
assumes that soot material is sufficiently hydrophilic, with high WS
wettability. Although we made assumptions in the modeling of these
particles, their shapes are similar to those observed by electron microscopy
(Adachi et al., 2010; Fig. 5). The shape model of the particle was
ultimately defined using a three-dimensional rendering technique.

We calculated the light-scattering properties using the finite-difference time-domain (FDTD) method (Cole, 2005; Taflove and Hagness, 2005; Ishimoto et al., 2012a) and DDA (DDSCAT version 7.3; https://code.google.com/archive/p/ddscat/ (last access: 25 December 2018); Draine and Flatau, 1994). Particles defined by discrete points as the input for the light-scattering calculations were reconstructed from the shape model described in Sect. 2, and the grid length was set within the range $\mathrm{2}\le \mathrm{2}a/l\le \mathrm{10}$ for correct reproduction of the defined aggregates within the numerical convergence criterion. Although soot particles are generally small aerosols, large computational resources for FDTD and DDA calculations are necessary to estimate the exact scattering properties of the mixed particles due to the fine structure of soot and additional volume contributed by WS material. Using the same discretized particle under an appropriate convergence condition, the scattering properties calculated by FDTD would have approximately the same accuracy as those derived from DDA (Yurkin et al., 2007). However, the numerical cost of the two methods differs depending on the particle shape and the refractive index of the particle material (Yurkin et al., 2007). The memory and CPU time of the DDA calculations mainly depend on the number of discretized points (i.e., dipoles) and depend on particle volume. Meanwhile, the numerical costs of the FDTD calculations are determined by the size of the discretized numerical field that encloses the particle. For fractal-shaped particles of relatively small fractal dimensions, DDA calculations are faster than FDTD ones because of the small relative volume with respect to size. By contrast, the numerical cost of DDA calculations drastically increases as the volume of attached WS elements increases, whereas the FDTD method can output results within similar CPU times given similar total sizes. Therefore, we used both numerical methods for the light-scattering calculations in accordance with the size and shape of the particles. In the numerical environment of our nonparallel computation, light-scattering calculations performed using FDTD were faster than those using DDA for AST-B particles (nos. 10–13).

As indicated in Sect. 2 and shown in Figs. 1 and 2, 10 sizes of bare soot aggregates with a volume-equivalent sphere radius ${r}_{\mathrm{eq}}=\mathrm{0.02}-\mathrm{0.20}\phantom{\rule{0.125em}{0ex}}\mathrm{\mu}\mathrm{m}$ for Type A and 13 sizes with ${r}_{\mathrm{eq}}=\mathrm{0.02}-\mathrm{0.42}\phantom{\rule{0.125em}{0ex}}\mathrm{\mu}\mathrm{m}$ for Type B were prepared, and internally mixed particles of ${V}_{\mathrm{r}}\sim \mathrm{0},\mathrm{2},\mathrm{5},\mathrm{10},\mathrm{20}$ were numerically created for each soot particle. The corresponding size ranges of the mixed soot particles were ${r}_{\mathrm{eq}}=\mathrm{0.02}-\mathrm{0.55}\phantom{\rule{0.125em}{0ex}}\mathrm{\mu}\mathrm{m}$ for Type A (AST-A) and ${r}_{\mathrm{eq}}=\mathrm{0.02}-\mathrm{1.15}\phantom{\rule{0.125em}{0ex}}\mathrm{\mu}\mathrm{m}$ for Type B (AST-B). Assuming a synthetic analysis using satellite- and ground-based multi-channel radiometer and lidar measurements, 10 wavelengths from near-ultraviolet to near-infrared (340, 355, 380, 400, 500, 532, 675, 870, 1020 and 1064 nm) were selected. We used a spectral refractive index dataset by Chang and Charalampopoulos (1990) for the bare soot material. For the WS components, the dependence of the refractive index on relative humidity was considered, and relative humidity values of 0 %, 50 %, 90 % and 98 % in the software package Optical Properties of Aerosols and Clouds (OPAC) were applied (Hess et al., 1998). The outputs included the results of the light-scattering properties with those averaged over 88 orientations for the FDTD method and 100 orientations for the DDA method.

As examples of the numerical results, Figs. 6 and 7 present the size
(*r*_{eq}) dependence of single-scattering albedo (*ω*) and the
asymmetry factor (*g*) at a wavelength of *λ*=532 nm for
AST-A. The complex refractive indices for soot and WS are listed in Table 1,
for which the refractive index at a relative humidity of 50 % was applied for
WS. For comparison, the results of *ω* and *g* with the same *V*_{r}
but derived via Mie calculations for spheres of the effective refractive
index calculated using the Maxwell-Garnett mixing rule (MG) (Bohren and
Huffman, 1983) and for spheres with a soot-core/WS-shell structure (CS) were
also plotted.

For *ω* and *g* at *λ*=532 nm, the results of the AST
model were approximately consistent with previous modeled results (Dong et
al., 2015; Liu et al., 2016). The results of *ω* for AST-A at *V*_{r}=0
markedly differed from those of MG and CS due to the volume-equivalent sphere
approximation adopted in the MG and CS treatments. The results of *ω* for
AST-A with *V*_{r}≥2 showed a similar trend to the MG and CS results such that
the MG results were closer to AST-A than CS results were. Regarding the asymmetry
factor, *g* depends mainly on the particle size *r*_{eq} and is less
sensitive to the mixing ratio *V*_{r}. The derived *g* for AST-A fell between
the results of MG and of CS, with the MG results closer to those of AST-A.
Because the primary particle of the assumed soot (*a*=0.02 µm)
was smaller than the wavelength and the Type A aggregates were fractal shapes
with relatively small *D*_{f} (Figs. 1 and 2), the effective medium
theory based on the MG mixing rule offered a better approximation than
the CS approximation did for the AST-A model.

Compared to *ω* and *g*, the backscattering properties of particles are
sensitive to particle shape and mixing state. Therefore, lidar measurements
could potentially offer information on the validity of the particle model.
The calculated lidar ratio *L* and linear depolarization ratio
*δ*_{L} for the AST-A and AST-B particles at wavelengths of 355,
532 and 1064 nm are plotted in Figs. 8–10. For a single particle, *L* is
calculated from *ω* and the normalized phase function *P*_{11} in the
backscattering direction, and *δ*_{L} is derived from the
*P*_{11} and *P*_{22} components of the scattering matrix. Here, we omitted
the backscattering of *P*_{12} for *δ*_{L} due to the assumption
of random particle orientation.

For lidar ratios, the MG and CS sphere results for the same volume of bare
soot aggregates and the same *r*_{eq} range of *V*_{r}≤20 as that for
AST-A and AST-B are also plotted in Figs. 8 and 9. Because backscattering
*P*_{11} is sensitive to the size of spherical particles, calculations were
performed at a step size of 0.01 µm. The lidar ratios of CS particles
tended to be smaller than those of MG particles, and the difference between
MG and CS was significant at wavelengths of 355 and 532 nm. Furthermore,
the lidar ratios of AST-A and AST-B particles were approximately between the
MG and CS results for particles of the same *r*_{eq} and mixing ratios,
particularly for entirely encapsulated *V*_{r}≥1 particles. As denoted in
the asymmetry factor results, MG and CS corresponded to two extreme cases of
mixture soot materials within the WS shell; our lidar ratio results for
fractal-like soot shapes among AST-A and AST-B are reasonable.

Due to strong absorption properties, bare soot particles showed small
depolarization ratios *δ*_{L} (Fig. 10). As the *V*_{r} value increased, *δ*_{L} of the AST particles increased through the mixing
of weak absorbing material (i.e., WS), and *δ*_{L} began to decrease
for larger values of *V*_{r} due to their spherical shape. This indicates
that the smooth surface created by a thin WS coating was rather ineffective
for reducing *δ*_{L} when the overall shape was highly nonspherical.
As a result, AST particles of the same bare soot size as those shown in
Fig. 10 (data of same color) had a peak depolarization ratio. To some
extent, peak values of *δ*_{L} were larger for larger bare soot
particles at shorter wavelengths. This implies that the particle size
parameter *x*_{eq} ($=\mathrm{2}\mathit{\pi}{r}_{\mathrm{eq}}/\mathit{\lambda}$) is an important factor for the
evaluation of *δ*_{L} as well as shape nonsphericity.

According to the biomass smoke measurement results, the size of dry particles
was 0.10–0.16 µm at count median diameter (0.25–0.30 µm
at volume median diameter) (Reid et al., 2005a). Laboratory and in situ
measurements for aged biomass smoke yielded estimates of *ω*≥0.75
and *g*∼0.6 at *λ*=532 nm (Reid et al., 2005b; Pokhrel et al.,
2016). For lidar measurements, typical values of *L*∼70 sr and
*δ*_{L}∼7 % at *λ*=532 nm have been reported
(Groß et al., 2015a, b). Among our optical property results for AST-A,
particle no. 4 (*V*_{r}=20) had values of $\left(\mathit{\omega},g,L,{\mathit{\delta}}_{\mathrm{L}}\right)=\left(\mathrm{0.87},\mathrm{0.67},\mathrm{69}\phantom{\rule{0.125em}{0ex}}\mathrm{sr},\mathrm{0}\phantom{\rule{0.125em}{0ex}}\mathit{\%}\right)$, and particle no. 5 (*V*_{r}=2) had values of $\left(\mathit{\omega},g,L,{\mathit{\delta}}_{\mathrm{L}}\right)=\left(\mathrm{0.54},\mathrm{0.63},\mathrm{188}\phantom{\rule{0.125em}{0ex}}\mathrm{sr},\mathrm{6}\phantom{\rule{0.125em}{0ex}}\mathit{\%}\right)$. Using the refractive index of a typical soot particle, relatively
high values of *V*_{r} would be necessary to simulate *ω*≥0.75. In the AST model, the soot aggregate was entirely encapsulated by WS
components, and the overall shape became spherical under large *V*_{r}
values. Although such spherical particles have consistent lidar ratios of
*L*∼70 sr, the depolarization ratio becomes *δ*_{L}∼0
for spherical shapes. By contrast, large depolarization ratios can occur if
internally mixed particles are highly nonspherical. Moreover, the reported
spectral dependence of depolarization ratios $\left({\mathit{\delta}}_{L,\mathrm{355}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}}\sim \mathrm{20}\phantom{\rule{0.125em}{0ex}}\mathit{\%},{\mathit{\delta}}_{L,\mathrm{532}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}}\sim \mathrm{9}\phantom{\rule{0.125em}{0ex}}\mathit{\%},{\mathit{\delta}}_{L,\mathrm{1064}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}}\sim \mathrm{2}\phantom{\rule{0.125em}{0ex}}\mathit{\%}\right)$ from airborne measurements
for smoke plumes (Burton et al., 2015; Mishchenko et al., 2016) can be
explained, for example, by the AST-A particle no. 6 (${V}_{\mathrm{r}}=\mathrm{5}-\mathrm{10}$;
Fig. 10a–c). However, such nonspherical mixed particles tended to have
relatively large lidar ratios of *L*≥100 sr at $\mathit{\lambda}=\mathrm{355},\mathrm{532}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}$ in the AST-A model. In the comparison between AST-A and
AST-B particles, variation in *L* and *δ*_{L} for AST-B was
greater than for AST-A. At *λ*=355 nm, the measurement-derived
lidar ratio was ${L}_{\mathrm{355}\phantom{\rule{0.125em}{0ex}}\mathrm{nm}}=\mathrm{76}\pm \mathrm{12}\phantom{\rule{0.125em}{0ex}}\mathrm{sr}$ (Groß et
al., 2015b). The value was difficult to explain for MG particles, but was
relatively consistent with the large *V*_{r} observed for AST-A and AST-B
particles and some CS particles. By contrast, the calculated depolarization
ratios of AST-A and AST-B particles at both *λ*=355 and 532 nm
exceeded 30 % for larger aggregates with a high amount of coating. These
ratios were higher than those measured for biomass burning
(7 %–16 %; Groß et
al., 2015b), suggesting that the depolarization ratios of nonspherical soot
particles were easily enhanced when coated with weakly absorbing WS material.
A similar effect was previously reported (Kahnert, 2017). According to
Kanngießer and Kahnert (2018), the speed of transition from film coating
to spherical growth (Pei et al., 2018) is a morphological parameter that
strongly affects the depolarization ratio. A rapid transition to spherical
growth may cause a smaller depolarization ratio of internally mixed soot
particles.

Overall, the results shown in Figs. 8–10 suggest that the effects of internal mixing on lidar backscattering are strongly related to changes in absorption and shape properties due to mixing and bare soot particle size. However, it was difficult to simulate average smoke optical properties $\left(\mathit{\omega},g,L,{\mathit{\delta}}_{\mathrm{L}}\right)$ and their spectral dependence solely using AST-A or AST-B particles.

Among the field measurements, the observed optical properties were the
average of the particle size distribution. Our results indicate that the
presence of large mixed soot particles may enhance the bulk *δ*_{L} of
smoke.

Typical values of *V*_{r} for mixed biomass smoke particles likely
depend on the relative humidity, concentrations of WS components and
the size and shape of soot aggregates. The mixing of different aerosol types, such as
dust in biomass burning (Groß et al., 2011, 2013), may be important for
the interpretation of measured optical properties. Improved retrieval
calculations based on realistic aerosol simulations that consider particle
size and shape distribution and other types of aerosol contamination are expected
using our AST model for internally mixed soot particles.

For satellite- and ground-based remote-sensing analysis applications, we developed a shape model of internally mixed soot particles and calculated the optical properties of particles at visible and near-infrared wavelengths. Fractal-like structures extracted using spatial Voronoi tessellation were considered to mimic necking and overlapping between neighboring primary particles. We created two types of polyhedral aggregates with different particle shape–size dependences (Type A and Type B) to account for the effects of compaction on soot aggregate shape during the atmospheric aging process. Then the artificial surface potential for the particle in a Cartesian grid space was defined, and the surface tension of the WS components on the soot aggregate was simulated, assuming that the soot was hydrophilic with high wettability. Based on a simple assumption of the behavior of WS elements, the shapes of internally mixed soot particles dependent on the amount of attached dissolved material were determined from iterative calculations. Overall, an aggregate coated with a thin film was simulated given a relatively small volume of added WS components, whereas the soot aggregate was covered with a spherical shell given a large number of WS components. The optical properties of the developed internally mixed particles were calculated using the FDTD and DDA methods while considering the spectral dependence on the refractive indices.

For single-scattering albedo and asymmetry, AST model results were similar
to those for MG and CS, except when the WS ratio was low (*V*_{r}∼0).
Due to shape irregularities, the lidar ratios of AST particles in random
orientation were less sensitive to particle size than those of MS and CS
particles. The *L* values of AST particles were approximately between those
of MG and CS particles. For AST particles, *δ*_{L} increased as the
amount of weakly absorbing material (i.e., WS) increased, and *δ*_{L}
decreased as particle shape became spherical. As a result, irregularly
shaped soot particles had peak depolarization ratios during the internal
mixing process. Maximum *δ*_{L} values tended to be larger when the
size parameters of nonspherical mixed particles were large. Following
comparisons with reported optical features, we determined that average
optical properties $\left(\mathit{\omega},g,L,{\mathit{\delta}}_{\mathrm{L}}\right)$ and their
spectral dependences for measured biomass burning aerosols cannot be
simulated by a single particle in our internally mixed soot model. In
particular, AST particles tended to have larger lidar ratios and
depolarization ratios than those obtained by field measurements. The
inconsistencies in lidar backscattering properties between the model results
and measurements may remain unresolved even after size-averaging the optical
properties. One possible explanation for this phenomenon is that soot
aggregates change their shapes to become more compact (i.e., become
spherical) during the WS adhesion process in the atmosphere. The
contribution of the core-shell type of internally mixed soot particles and
the use of the nonspherical model may be necessary to simulate observed
results for burned biomass aerosols.

In addition to its use in direct analyses for the field measurement results, the dataset of the AST particle optical properties can be used for the parameterization of conventional spherical particle models, such as the MG model, CS model, and their combined model. Furthermore, the dataset will be useful for determining the shape property conditions of smoke particles observed from multi-sensor measurements, including lidar backscattering.

The data in the figures are available upon request from the corresponding author (Hiroshi Ishimoto).

The supplement related to this article is available online at: https://doi.org/10.5194/amt-12-107-2019-supplement.

HI, RK and KA created particle shape model and finalized the paper. Light scattering properties were calculated by HI. The optical database was designed by RK.

The authors declare that they have no conflict of interest.

This article is part of the special issue “SKYNET – the international network for aerosol, clouds, and solar radiation studies and their applications (AMT/ACP inter-journal SI)”. It is not associated with a conference.

This research was supported by the Global Change Observation Mission-Climate
(GCOM-C) research project of the Japan Aerospace Exploration Agency (JAXA),
and by a grant for the Global Environmental Research Coordination System from
the Ministry of the Environment, Japan.

Edited
by: Teruyuki Nakajima

Reviewed by: three anonymous referees

Adachi, K. and Buseck, P. R.: Internally mixed soot, sulfates, and organic matter in aerosol particles from Mexico City, Atmos. Chem. Phys., 8, 6469–6481, https://doi.org/10.5194/acp-8-6469-2008, 2008.

Adachi, K., Chung, S. H., Friedrich, H., and Buseck, P. R.: Fractal parameters of individual soot particles determined using electron tomography: Implications for optical properties, J. Geophys. Res.-Atmos., 112, 1–10, https://doi.org/10.1029/2006JD008296, 2007.

Adachi, K., Chung, S. H., and Buseck, P. R.: Shapes of soot aerosol particles and implications for their effects on climate, J. Geophys. Res.-Atmos., 115, 1–9, https://doi.org/10.1029/2009JD012868, 2010.

Adachi, K., Sedlacek, A. J., Kleinman, L., Chand, D., Hubbe, J. M., and Buseck, P. R.: Volume changes upon heating of aerosol particles from biomass burning using transmission electron microscopy, Aerosol Sci. Technol., 52, 46–56, https://doi.org/10.1080/02786826.2017.1373181, 2018.

Akinci, N., Akinci, G., and Teschner, M.: Versatile surface tension and adhesion for SPH fluids, ACM Trans. Graph., 32, 182, https://doi.org/10.1145/2508363.2508395, 2013.

Baran, A. J., Ishimoto, H., Sourdeval, O., Hesse, E., and Harlow, C.: The applicability of physical optics in the millimetre and sub-millimetre spectral region, Part II: Application to a three-component model of ice cloud and its evaluation against the bulk single-scattering properties of various other aggregate models, J. Quant. Spectrosc. Ra., 206, 83–100, https://doi.org/10.1016/j.jqsrt.2017.10.027, 2018.

Becker, S., Urbassek, H. M., Horsch, M., and Hasse, H.: Contact Angle of Sessile Drops in Lennard-Jones Systems, Langmuir, 30, 13606–13614, https://doi.org/10.1021/la503974z, 2014.

Bohren, C. F. and Huffman, D. R.: Absorption and scattering of light by small particles, John Wiley & Sons, New York, 1983.

Bond, T. C. and Bergstrom, R. W.: Light absorption by carbonaceous particles: An investigative review, Aerosol Sci. Technol., 40, 27–67, https://doi.org/10.1080/02786820500421521, 2006.

Burton, S. P., Hair, J. W., Kahnert, M., Ferrare, R. A., Hostetler, C. A., Cook, A. L., Harper, D. B., Berkoff, T. A., Seaman, S. T., Collins, J. E., Fenn, M. A., and Rogers, R. R.: Observations of the spectral dependence of linear particle depolarization ratio of aerosols using NASA Langley airborne High Spectral Resolution Lidar, Atmos. Chem. Phys., 15, 13453–13473, https://doi.org/10.5194/acp-15-13453-2015, 2015.

Chang, H. and Charalampopoulos, T. T.: Determination of the Wavelength Dependence of Refractive Indices of Flame Soot, Proc. R. Soc. A Math. Phys. Eng. Sci., 430, 577–591, https://doi.org/10.1098/rspa.1990.0107, 1990.

Cheng, T., Gu, X., Wu, Y., and Chen, H.: Effects of atmospheric water on the optical properties of soot aerosols with different mixing states, J. Quant. Spectrosc. Ra., 147, 196–206, https://doi.org/10.1016/j.jqsrt.2014.06.002, 2014.

Cole, J. B.: High accuracy nonstandard finite-difference time-domain algorithms for computational electromagnetics: Applications to optics and photonics, in: Advances in the Applications of Nonstandard Finite Difference Schemes, edited by: Mickens, R. E., 89–189, World Scientific, Singapore, 2005.

Dong, J., Zhao, J. M., and Liu, L. H.: Morphological effects on the radiative properties of soot aerosols in different internally mixing states with sulfate, J. Quant. Spectrosc. Ra., 165, 43–55, https://doi.org/10.1016/j.jqsrt.2015.06.025, 2015.

Draine, B. T. and Flatau, P. J.: Discrete-Dipole Approximation For Scattering Calculations, J. Opt. Soc. Am. A, 11, 1491, https://doi.org/10.1364/JOSAA.11.001491, 1994.

Groß, S., Tesche, M., Freudenthaler, V., Toledano, C., Wiegner, M., Ansmann, A., Althausen, D., and Seefeldner, M.: Characterization of Saharan dust, marine aerosols and mixtures of biomass-burning aerosols and dust by means of multi-wavelength depolarization and Raman lidar measurements during SAMUM 2, Tellus B, 63, 706–724, https://doi.org/10.1111/j.1600-0889.2011.00556.x, 2011.

Groß, S., Esselborn, M., Weinzierl, B., Wirth, M., Fix, A., and Petzold, A.: Aerosol classification by airborne high spectral resolution lidar observations, Atmos. Chem. Phys., 13, 2487–2505, https://doi.org/10.5194/acp-13-2487-2013, 2013.

Groß, S., Freudenthaler, V., Schepanski, K., Toledano, C., Schäfler, A., Ansmann, A., and Weinzierl, B.: Optical properties of long-range transported Saharan dust over Barbados as measured by dual-wavelength depolarization Raman lidar measurements, Atmos. Chem. Phys., 15, 11067–11080, https://doi.org/10.5194/acp-15-11067-2015, 2015a.

Groß, S., Freudenthaler, V., Wirth, M., and Weinzierl, B.: Towards an aerosol classification scheme for future EarthCARE lidar observations and implications for research needs, Atmos. Sci. Lett., 16, 77–82, https://doi.org/10.1002/asl2.524, 2015b.

Hara, Y., Nishizawa, T., Sugimoto, N., Osada, K., Yumimoto, K., Uno, I., Kudo, R., and Ishimoto, H.: Retrieval of Aerosol Components Using Multi-Wavelength Mie-Raman Lidar and Comparison with Ground Aerosol Sampling, Remote Sens., 10, 937, https://doi.org/10.3390/rs10060937, 2018.

He, C., Liou, K.-N., Takano, Y., Zhang, R., Levy Zamora, M., Yang, P., Li, Q., and Leung, L. R.: Variation of the radiative properties during black carbon aging: theoretical and experimental intercomparison, Atmos. Chem. Phys., 15, 11967–11980, https://doi.org/10.5194/acp-15-11967-2015, 2015.

Hess, M., Koepke, P., and Schult, I.: Optical Properties of Aerosols and Clouds: The Software Package OPAC, B. Am. Meteor. Soc., 79, 831–844, https://doi.org/10.1175/1520-0477(1998)079<0831:OPOAAC>2.0.CO;2, 1998.

Ishimoto, H., Zaizen, Y., Uchiyama, A., Masuda, K., and Mano, Y.: Shape modeling of mineral dust particles for light-scattering calculations using the spatial Poisson-Voronoi tessellation, J. Quant. Spectrosc. Ra., 111, 2434–2443, https://doi.org/10.1016/j.jqsrt.2010.06.018, 2010.

Ishimoto, H., Masuda, K., Mano, Y., Orikasa, N., and Uchiyama, A.: Irregularly shaped ice aggregates in optical modeling of convectively generated ice clouds, J. Quant. Spectrosc. Ra., 113, 632–643, https://doi.org/10.1016/j.jqsrt.2012.01.017, 2012a.

Ishimoto, H., Zaizen, Y., Masuda, K., Mano, Y., and Uchiyama, A.: Shape modeling of dust and soot particles for remote sensing applications considering the geometrical features of sampled aerosols, in: Technical Reports of the Meteorological Research Institute, 40–43, available at: http://www.mri-jma.go.jp/Publish/Technical/index_en.html (last access: 26 December 2018), 2012b.

Jacobson, M. Z.: Strong radiative heating due to the mixing state of black carbon in atmospheric aerosols, Nature, 409, 695–697, https://doi.org/10.1038/35055518, 2001.

Kahnert, M.: Optical properties of black carbon aerosols encapsulated in a shell of sulfate: comparison of the closed cell model with a coated aggregate model, Opt. Express, 25, 24579, https://doi.org/10.1364/OE.25.024579, 2017.

Kahnert, M., Nousiainen, T., and Lindqvist, H.: Models for integrated and differential scattering optical properties of encapsulated light absorbing carbon aggregates, Opt. Express, 21, 7974–7993, 2013.

Kanngießer, F. and Kahnert, M.: Calculation of optical properties of light-absorbing carbon with weakly absorbing coating: A model with tunable transition from film-coating to spherical-shell coating, J. Quant. Spectrosc. Ra., 216, 17–36, https://doi.org/10.1016/j.jqsrt.2018.05.014, 2018.

Leinonen, J. and von Lerber, A.: Snowflake Melting Simulation Using Smoothed Particle Hydrodynamics, J. Geophys. Res.-Atmos., 123, 1811–1825, https://doi.org/10.1002/2017JD027909, 2018.

Liu, F., Yon, J., and Bescond, A.: On the radiative properties of soot aggregates – Part 2: Effects of coating, J. Quant. Spectrosc. Ra., 172, 134–145, https://doi.org/10.1016/j.jqsrt.2015.08.005, 2016.

Luo, J., Zhang, Y., and Zhang, Q.: A model study of aggregates composed of spherical soot monomers with an acentric carbon shell, J. Quant. Spectrosc. Ra., 205, 184–195, https://doi.org/10.1016/j.jqsrt.2017.10.024, 2018.

Mikhailov, E. F., Vlasenko, S. S., Krämer, L., and Niessner, R.: Interaction of soot aerosol particles with water droplets: Influence of surface hydrophilicity, J. Aerosol Sci., 32, 697–711, https://doi.org/10.1016/S0021-8502(00)00101-4, 2001.

Mikhailov, E. F., Vlasenko, S. S., Podgorny, I. A., Ramanathan, V., and Corrigan, C. E.: Optical properties of soot-water drop agglomerates: An experimental study, J. Geophys. Res.-Atmos., 111, D07209, https://doi.org/10.1029/2005JD006389, 2006.

Mishchenko, M. I., Dlugach, J. M., and Liu, L.: Linear depolarization of lidar returns by aged smoke particles, Appl. Opt., 55, 9968–9973, https://doi.org/10.1364/AO.55.009968, 2016.

Moteki, N.: Discrete dipole approximation for black carbon-containing aerosols in arbitrary mixing state: A hybrid discretization scheme, J. Quant. Spectrosc. Ra., 178, 306–314, https://doi.org/10.1016/j.jqsrt.2016.01.025, 2016.

Moteki, N. and Kondo, Y.: Effects of mixing state on black carbon measurements by laser-induced incandescence, Aerosol Sci. Technol., 41, 398–417, https://doi.org/10.1080/02786820701199728, 2007.

Nyeki, S. and Colbeck, I.: Fractal Dimension Analysis of Single, In-Situ, Restructured Carbonaceous Aggregates, Aerosol Sci. Technol., 23, 109–120, https://doi.org/10.1080/02786829508965298, 1995.

Ohser, J. and Mücklich, F.: Statistical Analysis of Microstructures in Materials Science, John Wiley & Sons, Chichester, 2000.

Okyay, G., Héripré, E., Reiss, T., Haghi-Ashtiani, P., Auger, T., and Enguehard, F.: Soot aggregate complex morphology: 3D geometry reconstruction by SEM tomography applied on soot issued from propane combustion, J. Aerosol Sci., 93, 63–79, https://doi.org/10.1016/j.jaerosci.2015.11.009, 2016.

Pei, X., Hallquist, M., Eriksson, A. C., Pagels, J., Donahue, N. M., Mentel, T., Svenningsson, B., Brune, W., and Pathak, R. K.: Morphological transformation of soot: investigation of microphysical processes during the condensation of sulfuric acid and limonene ozonolysis product vapors, Atmos. Chem. Phys., 18, 9845–9860, https://doi.org/10.5194/acp-18-9845-2018, 2018.

Pokhrel, R. P., Wagner, N. L., Langridge, J. M., Lack, D. A., Jayarathne, T., Stone, E. A., Stockwell, C. E., Yokelson, R. J., and Murphy, S. M.: Parameterization of single-scattering albedo (SSA) and absorption Ångström exponent (AAE) with EC∕OC for aerosol emissions from biomass burning, Atmos. Chem. Phys., 16, 9549–9561, https://doi.org/10.5194/acp-16-9549-2016, 2016.

Ramanathan, V. and Carmichael, G.: Global and regional climate changes due to black carbon, Nat. Geosci., 1, 221–227, https://doi.org/10.1038/ngeo156, 2008.

Reid, J. S., Koppmann, R., Eck, T. F., and Eleuterio, D. P.: A review of biomass burning emissions part II: intensive physical properties of biomass burning particles, Atmos. Chem. Phys., 5, 799–825, https://doi.org/10.5194/acp-5-799-2005, 2005a.

Reid, J. S., Eck, T. F., Christopher, S. A., Koppmann, R., Dubovik, O., Eleuterio, D. P., Holben, B. N., Reid, E. A., and Zhang, J.: A review of biomass burning emissions part III: intensive optical properties of biomass burning particles, Atmos. Chem. Phys., 5, 827–849, https://doi.org/10.5194/acp-5-827-2005, 2005b.

Reid, S. and Hobbs, P. V: Physical and optical properties of young smoke from individual biomass fires in Brazil, J. Geophys. Res., 103, 32013–32030, 1998.

Scarnato, B. V., Vahidinia, S., Richard, D. T., and Kirchstetter, T. W.: Effects of internal mixing and aggregate morphology on optical properties of black carbon using a discrete dipole approximation model, Atmos. Chem. Phys., 13, 5089–5101, https://doi.org/10.5194/acp-13-5089-2013, 2013.

Shiraiwa, M., Kondo, Y., Moteki, N., Takegawa, N., Miyazaki, Y., and Blake, D. R.: Evolution of mixing state of black carbon in polluted air from Tokyo, Geophys. Res. Lett., 34, 2–6, https://doi.org/10.1029/2007GL029819, 2007.

Shiraiwa, M., Kondo, Y., Iwamoto, T., and Kita, K.: Amplification of light absorption of black carbon by organic coating, Aerosol Sci. Technol., 44, 46–54, https://doi.org/10.1080/02786820903357686, 2010.

Taflove, A. and Hagness, S. C.: Computational electrodynamics The finite-difference time-domain method, 3rd edn., Artech House, Boston, 2005.

Tartakovsky, A. and Meakin, P.: Modeling of surface tension and contact angles with smoothed particle hydrodynamics, Phys. Rev. E, 72, 026301, https://doi.org/10.1103/PhysRevE.72.026301, 2005.

Wu, Y., Cheng, T., Zheng, L., Chen, H., and Xu, H.: Single scattering properties of semi-embedded soot morphologies with intersecting and non-intersecting surfaces of absorbing spheres and non-absorbing host, J. Quant. Spectrosc. Ra., 157, 1–13, https://doi.org/10.1016/j.jqsrt.2015.02.006, 2015.

Wu, Y., Cheng, T., Zheng, L., and Chen, H.: Optical properties of the semi-external mixture composed of sulfate particle and different quantities of soot aggregates, J. Quant. Spectrosc. Ra., 179, 139–148, https://doi.org/10.1016/j.jqsrt.2016.03.012, 2016.

Wu, Y., Cheng, T., Zheng, L., and Chen, H.: Sensitivity of mixing states on optical properties of fresh secondary organic carbon aerosols, J. Quant. Spectrosc. Ra., 195, 147–155, https://doi.org/10.1016/j.jqsrt.2017.01.013, 2017.

Yon, J., Bescond, A., and Liu, F.: On the radiative properties of soot aggregates part 1: Necking and overlapping, J. Quant. Spectrosc. Ra., 162, 197–206, https://doi.org/10.1016/j.jqsrt.2015.03.027, 2015.

Yurkin, M. A., Hoekstra, A. G., Brock, R. S., and Lu, J. Q.: Systematic comparison of the discrete dipole approximation and the finite difference time domain method for large dielectric scatterers, Opt. Express, 15, 17902–17911, 2007.

Zhang, X., Mao, M., Yin, Y., and Wang, B.: Absorption enhancement of aged black carbon aerosols affected by their microphysics: A numerical investigation, J. Quant. Spectrosc. Ra., 202, 90–97, https://doi.org/10.1016/j.jqsrt.2017.07.025, 2017.

Zhang, Y., Zhang, Q., Cheng, Y., Su, H., Kecorius, S., Wang, Z., Wu, Z., Hu, M., Zhu, T., Wiedensohler, A., and He, K.: Measuring the morphology and density of internally mixed black carbon with SP2 and VTDMA: new insight into the absorption enhancement of black carbon in the atmosphere, Atmos. Meas. Tech., 9, 1833–1843, https://doi.org/10.5194/amt-9-1833-2016, 2016.