Modeling single-scattering properties of small cirrus particles by use of a size-shape distribution of ice spheroids and cylinders

In this study, we model single-scattering properties of small cirrus crystals using mixtures of polydisperse, randomly oriented spheroids and cylinders with varying aspect ratios and with a refractive index representative of water ice at a wavelength of 1.88 μm. The Stokes scattering matrix elements averaged over wide shape distributions of spheroids and cylinders are compared with those computed for polydisperse surface-equivalent spheres. The shape-averaged phase function for a mixture of oblate and prolate spheroids is smooth, featureless, and nearly flat at side-scattering angles and closely resembles those typically measured for cirrus. Compared with the ensemble-averaged phase function for spheroids, that for a shape distribution of cylinders shows a relatively deeper minimum at side-scattering angles. This may indicate that light scattering from realistic cirrus crystals can be better represented by a shape mixture of ice spheroids. Interestingly, the single-scattering properties of shape-averaged oblate and prolate cylinders are very similar to those of compact cylinders with a diameter-to-length ratio of unity. The differences in the optical cross sections, single-scattering albedo, and asymmetry parameter between the spherical and the nonspherical particles studied appear to be relatively small. This may suggest that for a given optical thickness, the influence of particle shape on the radiative forcing caused by a cloud composed of small ice crystals can be negligible.     

Scattering phase functions of horizontally oriented hexagonal ice crystals

Finite-difference time domain (FDTD) solutions are first compared with the corresponding T-matrix results for light scattering by circular cylinders with specific orientations. The FDTD method is then utilized to study the scattering properties of horizontally oriented hexagonal ice plates at two wavelengths, 0.55 and 12 μm. The phase functions of horizontally oriented ice plates deviate substantially from their counterparts obtained for randomly oriented particles. Furthermore, we compute the phase functions of horizontally oriented ice crystal columns by using the FDTD method along with two schemes for averaging over the particle orientations. It is shown that the phase functions of hexagonal ice columns with horizontal orientations are not sensitive to the rotation about the principal axes of the particles. Moreover, hexagonal ice crystals and circular cylindrical ice particles have similar optical properties, particularly, at a strongly absorbing wavelength, if the two particle geometries have the same length and aspect ratio defined as the ratio of the radius or semi-width of the cross section of a particle to its length. The phase functions for the two particle geometries are slightly different in the case of weakly absorbing plates with large aspect ratios. However, the solutions for circular cylinders agree well with their counterparts for hexagonal columns.     

Can particle shape information be retrieved from light-scattering observations using spheroidal model particles?

We address the question if and how observations of scattered intensity and polarisation can be employed for retrieving particle shape information beyond a simple classification into spherical and nonspherical particles. To this end, we perform several numerical experiments, in which we attempt to retrieve shape information of complex particles with a simple nonspherical particle model based on homogeneous spheroids. The discrete dipole approximation is used to compute reference phase matrices for a cube, a Gaussian random sphere, and a porous oblate and prolate spheroid as a function of size parameter. Phase matrices for the model particles, homogeneous spheroids, are computed with the T-matrix method. By assuming that the refractive index and the size distribution is known, an optimal shape distribution of model particles is sought that best matches the reference phase matrix. Both the goodness of fit and the optimal shape distribution are analysed. It is found that the phase matrices of cubes and Gaussian random spheres are well reproduced by the spheroidal particle model, while the porous spheroids prove to be challenging. The “retrieved” shape distributions, however, do not correlate well with the shape of the target particle even when the phase matrix is closely reproduced. Rather, they tend to exaggerate the aspect ratio and always include multiple spheroids. A most likely explanation why spheroids succeed in mimicking phase matrices of more irregularly shaped particles, even if their shape distributions display little similarity to those of the target particles, is that by varying the spheroids’ aspect ratio one covers a large range of different phase matrices. This often makes it possible to find a shape distribution of spheroids that matches the phase matrix of more complex particles.     

Sizing of spheroidal and cylindrical particles in a binary mixture by measurement of scattered light intensity: application of neural networks

The problem on the retrieval of sizes of an individual optically soft particle taken from binary mixtures of either oblate and prolate spheroids or cylinders and oblate spheroids is considered. It is based on multiangle scattered light intensity data. The multilevel neural networks method with a linear activation function and the method of the discrimination functions are used. Neural networks to retrieve characteristics of cylinders, oblate and prolate spheroids are designed. The errors in retrieved particle characteristics are investigated for the radius of an equivolume sphere in the range of 0.3-, shape parameter of spheroidal and cylindrical particles from -0.5 to 0.5 and 0 to 0.5, respectively.     

An Ellipsoidal Model for Small Multilayer Particles

This paper presents an ellipsoidal model that is constructed for small layered nonspherical particles and methods for constructing “effective” multilayer ellipsoids, the light-scattering properties of which would be close to the corresponding properties of original particles. In the case of axisymmetric particles, prolate or oblate spheroids (ellipsoids of revolution) are implied. Numerical calculations of the polarizability and scattering cross sections of small layered nonspherical particles, including nonconfocal (similar) spheroids, Chebyshev particles, and pseudospheroids, are performed by different approximate and rigorous methods. Approximate approaches involve the use of an ellipsoidal model, in which the polarizability of a layered particle is determined in two ways. In the first case, the polarizability is calculated in the approximation of confocal spheroids, while, in the second case, it is sought as a linear combination of the polarizabilities of embedded spheroids proportionally to the volumes of layers. Among rigorous methods, the extended boundary conditions method and the generalized separation of variables method are applied. On the basis of a comparison of the results obtained with rigorous and approximate approaches, their drawbacks and advantages are discussed.     

Polarization properties of light multiply scattered by non-spherical Rayleigh particles

Multiple scattering of incoherent polarized light propagating through a random medium comprised of spheroidal Rayleigh particles is studied using Monte Carlo simulations. Two approaches are taken for the implementation of the simulation: the first uses individual realizations of particle orientation and the second, an accelerated method, averages over the particle orientation. These different methods produce results that are indistinguishable within statistical errors. The depolarization of light is examined in both transmission and backscatter for media comprised of spheroids of different polarizability ratios. In media containing spheroidal particles the depolarization is greater than that for spherical particles. Media containing prolate spheroids are more depolarizing than media comprising oblate particles of the same polarizability ratio. The extra depolarization due to asphericity is much less pronounced in the multiple scattering regime than for single scattering.     

Polarization properties of light multiply scattered by non-spherical Rayleigh particles

Abstract

Multiple scattering of incoherent polarized light propagating through a random medium comprised of spheroidal Rayleigh particles is studied using Monte Carlo simulations. Two approaches are taken for the implementation of the simulation: the first uses individual realizations of particle orientation and the second, an accelerated method, averages over the particle orientation. These different methods produce results that are indistinguishable within statistical errors. The depolarization of light is examined in both transmission and backscatter for media comprised of spheroids of different polarizability ratios. In media containing spheroidal particles the depolarization is greater than that for spherical particles. Media containing prolate spheroids are more depolarizing than media comprising oblate particles of the same polarizability ratio. The extra depolarization due to asphericity is much less pronounced in the multiple scattering regime than for single scattering.     

Radiative properties of cirrus clouds in the infrared (8–13 μm) spectral region

Atmospheric radiation in the infrared (IR) 8–13 μm spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time-domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 μm. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 and 10,000 μm over the wavelengths ranging from 8 to 13 μm. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 μm. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8–13 μm spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, the extinction and absorption efficiencies are computed for hexagonal ice crystals with sizes ranging from 1 to 10,000 μm at 12 wavelengths between 8 and 13 μm.

Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Ice crystals are assumed to be hexagonal columns randomly oriented in space. The bulk scattering properties are parameterized through the use of second-order polynomial functions for the extinction efficiency and the single-scattering albedo and a power-law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 μm whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 μm. For effective sizes larger than 100 μm, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing studies of ice clouds involving IR window bands.     

Modeling the radiative properties of nonspherical soil-derived mineral aerosols

Mineral dust aerosols have complex nonspherical shapes and varying composition. This study utilizes data on morphology (size and shape) and composition of dust particles to determine the extent to which the optical properties of real particles differ from those of spheres. A method for modeling the optical properties of complex particle mixtures is proposed. The method combines dust particle composition-shape-size (CSS) distributions reconstructed from the electron microscopy data, effective medium approximations and discrete dipole approximation. The method is used to compute optical characteristics of realistic dust mixtures representative of Saharan and Asian dust. We demonstrate that considered CSS distributions result in various differences in the extinction coefficient, single scattering albedo, asymmetry parameter and the scattering phase function relative to the volume-equivalent spheres and the mixtures of the randomly oriented oblate and prolate spheroids. Implications of these differences for radiation/climate modeling and remote sensing are discussed.     

随机取向非球形沙尘气溶胶粒子的光学特性

利用离散偶极子近似法分析了一种随机取向旋转椭球体沙尘气溶胶粒子模型在尺度参数变化范围为0.1~23时(波长0.55!m对应有效半径为0.01~2!m)的光学特性,研究了沙尘粒子非球形性程度对其光学特性的影响,并考察了非球形粒子的随机取向能否用等体积球体来代替。就随机取向单分散和多分散旋转椭球体沙尘气溶胶而言,粒子非球形特征越明显,消光效率因子、不对称因子和单次散射反照率基本上偏离其等体积球体越大;对于相同的非球形,不对称因子偏离其等体积球体的相对偏差要比消光效率因子和单次散射反照率要大。非球形粒子的随机取向并不能使其光学特性严格等效为其等体积球体的光学特性。如果粒子形状偏离球体较小,则非球形粒子的随机取向的平均效果能使其消光效率因子、不对称因子和单次散射反照率近似用等体积球体的对应光学参量来等效;而如果粒子形状偏离球形较大,仅有单次散射反照率可以近似用等体积球体的单次散射反照率来等效,例如,轴半径比为16的旋转椭球体沙尘粒子的单次散射反照率偏离其等体积球体仅在3%以内。     

Application of high gradient magnetic separation principles to magnetic drug targeting

A hypothetical magnetic drug targeting system, utilizing high gradient magnetic separation (HGMS) principles, was studied theoretically using FEMLAB simulations. This new approach uses a ferromagnetic wire placed at a bifurcation point inside a blood vessel and an externally applied magnetic field, to magnetically guide magnetic drug carrier particles (MDCP) through the circulatory system and then to magnetically retain them at a target site. Wire collection (CE) and diversion (DE) efficiencies were defined and used to evaluate the system performance. CE and DE both increase as the strength of the applied magnetic field (0.3–2.0 T), the amount of ferromagnetic material (iron) in the MDCP (20–100%) and the size of the MDCP (1–10 μm radius) increase, and as the average inlet velocity (0.1–0.8 m s⁻¹), the size of the wire (50–250 μm radius) and the ratio (4–10) of the parent vessel radius (0.25–1.25 mm radius) to wire radius decrease. The effect of the applied magnetic field direction (0° and 90°) on CE and DE was minimal. Under these plausible conditions, CEs as high as 70% were obtained, with DEs reaching only 30%; however, when the MDCPs were allowed to agglomerate (4–10 μm radius), CEs and DEs of 100% were indeed achieved. These results reveal that this new magnetic drug targeting approach for magnetically collecting MDCPs at a target site, even in arteries with very high velocities, is feasible and very promising; this new approach for magnetically guiding MDCPs through the circulatory system is also feasible but more limited. Overall, this study shows that magnetic drug targeting, based on HGMS principles, has considerable promise as an effective drug targeting tool with many potential applications.     

Light scattering modeling of small feldspar aerosol particles using polyhedral prisms and spheroids

The use of simplified particle shapes for modeling scattering by irregularly shaped mineral-dust particles is studied using polyhedral prisms and spheroids as model particles. Simulated phase matrices averaged over shape and size distributions at wavelength 633 nm are compared with a laboratory-measured phase matrix of feldspar particles with known size distribution with effective radius of . When an equi-probable shape distribution is assumed, prisms and oblate spheroids agree with measurements to a similar degree, whereas prolate spheroids perform markedly better. Both spheroids and prisms perform much better than spheres. When an automatic fitting method is applied for finding optimal shape distributions, it is found that the most elongated spheroids are most important for good fits, whereas nearly-spherical spheroids are generally of very little importance. The phase matrices for the different polyhedral prisms, on the other hand, are found to be similar, thus their shape-averaged phase matrices are insensitive to the shape distribution assumed. For spheroids, a simple parameterization for the shape distribution, where weights increase with increasing departure from spherical shape, is proposed and tested. This parameterization improves the fit of most phase matrix elements attained with an equi-probable shape distribution, and it performs particularly well for reproducing the measured phase function.     

Using simple particle shapes to model the Stokes scattering matrix of ensembles of wavelength-sized particles with complex shapes: possibilities and limitations

We investigate to what extent the full Stokes scattering matrix of an ensemble of wavelength-sized particles with complex shapes can be modeled by employing an ensemble of simple model shapes, such as spheres, spheroids, and circular cylinders. We also examine to what extent such a simple-shape particle model can be used to retrieve meaningful shape information about the complex-shaped particle ensemble. More specifically, we compute the Stokes scattering matrix for ensembles of randomly oriented particles having several polyhedral prism geometries of different sizes and shape parameters. These ensembles serve as proxies for size-shape mixtures of particles containing several different shapes of higher geometrical complexity than the simple-shaped model particles we employ. We find that the phase function of the complex-shaped particle ensemble can be accurately modeled with a size distribution of volume-equivalent spheres. The diagonal elements of the scattering matrix are accurately reproduced with a size-shape mixture of spheroids. A model based on circular cylinders accurately fits the full scattering matrix including the off-diagonal elements. However, the modeling results provide us with only a rough estimate of the effective shape parameter of the complex-shaped particle ensemble to be modeled. They do not allow us to infer detailed information about the shape distribution of the complex-shaped particle ensemble.     

Particle sizing by multiangle light-scattering data using the high-order neural networks

The problem of retrieval of size and refractive index of a spherical particle by angular dependence of scattered light in scanning flow cytometry is considered. For its solution, the high-order neural networks are used. We restricted the range of angles available for measurement from 10° to 60°. The retrieval errors of characteristics of nonabsorbing particles were investigated at the ranges of the radius and relative refractive index 0.6–10.6 μm, and 1.02–1.38, respectively.     

Retrieval of cirrus particle sizes using a split-window technique: a sensitivity study

This paper examines the sensitivity of retrieved ice particle sizes using split-window method to the light scattering program for the single scattering calculation. We find that for randomly oriented hexagonal ice particles the retrieval algorithm using the anomalous diffraction theory (ADT) significantly overestimates the mean effective ice particle sizes, D_ge. The retrieved D_ge based on the geometric optics method (GOM) and Mie theory agrees with reference results within 20% when D_ge<30 μm. Based on the speculation that there is no “tunneling” for complex particles, some recent studies suggest that the ADT is an appropriate method to simulate the absorption coefficient for irregularly shaped particles in the infrared. In this study, however, we find that the overestimation of D_ge due to the ADT is largely caused by the neglect of refraction and reflection processes, instead of by the neglect of “tunneling” in the absorption calculations. By considering complex particle shapes such as aggregates with surface roughness, we further find that the retrieved D_ge based on the GOM is not sensitive to the particle shapes. Note that both ADT and GOM do not consider the “tunneling”, but the retrieved D_ge based on the ADT is about two times larger than those based on the GOM. “Tunneling” plays a significant role in the retrieved D_ge only when the D_ge is larger than 30–35 μm. In this study, we also examine the sensitivities of retrieved D_ge to the ice particle size distributions assumed in the retrieval algorithm and to the errors in the emissivities. It is found that when the D_ge is larger than 30–40 μm, the retrieved D_ge becomes very sensitive to the uncertainties related to the ice particle size distributions and to the errors in cirrus emissivities derived from measurements.     

Severe loss of precision in calculations of T-matrix integrals

A severe loss of precision is unravelled in the numerical calculation of surface integrals that appear in the Extended Boundary Condition Method (EBCM), to calculate the T-matrix elements of axisymmetric particles. We systematically study the occurrence of numerical cancellations for three basic particle shapes, namely cylinders, spheroids, and offset spheres, with typical sizes, aspect ratios and materials often studied as benchmark examples in the literature. The cancellations are evidenced both for spheroids and offset spheres, and are particularly pronounced in the latter case. The resulting loss of precision is independent from the commonly asserted problems of matrix inversion. We show that the origin of these severe cancellations can be further studied and understood by numerical investigations of the scaling of the integrands and integrals with respect to the particle size parameter. This allows us to develop a detailed mathematical proof of these cancellations. The results suggest that the EBCM method, in its usual formulation, suffers important numerical instabilities which reduce the domain of convergence for specific particle shapes that are commonly used for testing and benchmarking the method.     

Interaction between ultra short pulses and a dense scattering medium by Monte Carlo simulation: consideration of particle size effect

With an enough short-pulse incident to an individual particle, elementary scattering modes can be observed: internal or external reflection, refraction and diffraction. Simulation of pulse propagation in dense scattering medium is usually computed for large observation time, so that time delays of pulse interaction with the particles are negligible compared to propagation times between particle. A Monte Carlo method is proposed to compute the propagation of an incident 100 fs laser pulse in dense medium taking into account time-dependent scattering characteristics of particle: observation time of scattered light is less than 5000 fs. Two extreme cases are exemplified: predominance of direct and single-scattered photons appears in a thin time window for small particles (1 μm). On the contrary multiple scattering is always predominant and scrambles the transmitted signal for large particles (100 μm).     

On the impact of non-sphericity and small-scale surface roughness on the optical properties of hematite aerosols

We perform a comparative modelling study to investigate how different morphological features influence the optical properties of hematite aerosols. We consider high-order Chebyshev particles as a proxy for aerosol with a small-scale surface roughness, and spheroids as a model for nonspherical aerosols with a smooth boundary surface. The modelling results are compared to those obtained for homogeneous spherical particles. It is found that for hematite particles with an absorption efficiency of order unity the difference in optical properties between spheres and spheroids disappears. For optically softer particles, such as ice particles at far-infrared wavelengths, this effect can be observed for absorption efficiencies lower than unity. The convergence of the optical properties of spheres and spheroids is caused by absorption and quenching of internal resonances inside the particles, which depend both on the imaginary part of the refractive index and on the size parameter, and to some extent on the real part of the refractive index. By contrast, small-scale surface roughness becomes the dominant morphological feature for large particles. This effect is likely to depend on the amplitude of the surface roughness, the relative significance of internal resonances, and possibly on the real part of the refractive index. The extinction cross section is rather insensitive to surface roughness, while the single-scattering albedo, asymmetry parameter, and the Mueller matrix are strongly influenced. Small-scale surface roughness reduces the backscattering cross section by up to a factor of 2-3 as compared to size-equivalent particles with a smooth boundary surface. This can have important implications for the interpretation of lidar backscattering observations.     

Feasibility of a Lidar Utilizing the Glory for Measuring Particle Size of Water Clouds

A new lidar method for measuring water cloud particle size is proposed, and the feasibility of the measurement is discussed. The method utilizes the phenomenon known as the glory which is observed in open air. The proposed lidar consists of a multicolor laser transmitter and two receiver systems looking at the scattering from the target cloud with different scattering angles. Results of the theoretical study show that a system with five laser wavelengths (355, 532, 750, 1064 and 1500 nm) and two receivers located at scattering angles of 180 and 177.5–179 deg is useful for measuring particle size (mode radius of the size distribution) in a range of 4 to 12μm.     

Near-field calculation based on the T-matrix method with discrete sources

We propose a novel approach to compute the field scattered by a particle in the near-zone, in the framework of the transition matrix formulation. This method is based on the expansion of the total near-field in terms of discrete sources vector spherical wave functions and turns out to be particularly effective to model strongly elongated or flattened particles with high permittivity. The performances are evaluated on gold spheroids with several aspect ratios. This method results very useful to understand the complex behavior of enhanced plasmon fields in resonant metallic nanostructures.