Scattering of light by polydisperse, randomly oriented, finite circular cylinders

We use the T-matrix method, as described by Mishchenko [Appl. Opt. 32, 4652 (1993)], to compute rigorously light scattering by finite circular cylinders in random orientation. First we discuss numerical aspects of T -matrix computations specific for finite cylinders and present results of benchmark computations for a simple cylinder model. Then we report results of extensive computations for polydisperse, randomly oriented cylinders with a refractive index of 1.53 + 0.008i, diameter-to-length ratios of 1/2, 1/1.4, 1, 1.4, and 2, and effective size parameters ranging from 0 to 25. These computations parallel our recent study of light scattering by polydisperse, randomly oriented spheroids and are used to compare scattering properties of the two classes of simple convex particles. Despite the significant difference in shape between the two particle types (entirely smooth surface for spheroids and sharp rectangular edges for cylinders), the comparison shows rather small differences in the integral photometric characteristics (total optical cross sections, single-scattering albedo, and asymmetry parameter of the phase function) and the phase function. The general patterns of the other elements of the scattering matrix for cylinders and aspect-ratio-equivalent spheroids are also qualitatively similar, although noticeable quantitative differences can be found in some particular cases. In general, cylinders demonstrate much less shape dependence of the elements of the scattering matrix than do spheroids. Our computations show that, like spheroids and bispheres, cylinders with surface-equivalent radii smaller than a wavelength can strongly depolarize backscattered light, thus suggesting that backscattering depolarization for nonspherical particles cannot be universally explained by using only geometric-optics considerations.     

Scattering of electromagnetic waves by spheroidal particles: a novel approach exploiting the T matrix computed in spheroidal coordinates

A method other than the extended-boundary-condition method (EBCM) to compute the T matrix for electromagnetic scattering is presented. The separation-of-variables method (SVM) is used to solve the electromagnetic scattering problem for a spheroidal particle and to derive its T matrix in spheroidal coordinates. A transformation is developed for transforming the T matrix in spheroidal coordinates into the corresponding T matrix in spherical coordinates. The T matrix so obtained can be used for analytical calculation of the optical properties of ensembles of randomly oriented spheroids of arbitrary shape by use of an existing method to average over orientational angles. The optical properties obtained with the SVM and the EBCM are compared for different test cases. For mildly aspherical particles the two methods yield indistinguishable results. Small differences appear for highly aspherical particles. The new approach can be used to compute optical properties for arbitrary values of the aspect ratio. To test the accuracy of the expansion coefficients of the spheroidal functions for arbitrary arguments, a new testing method based on the completeness relation of the spheroidal functions is developed.     

Scattering matrix of infrared radiation by ice finite circular cylinders

Scattering matrix characteristics of polydisperse, randomly oriented, small ice crystals modeled by finite circular cylinders with various ratios of the length to diameter (L/D) ratio are calculated by use of the exact T-matrix approach, with emphasis on the thermal infrared spectral region that extends from the atmospheric short-wave IR window to the far-IR wavelengths to as large as 30 microm. The observed ice crystal size distribution and the well-known power-law distribution are considered. The results of the extensive calculations show that the characteristics of scattering matrix elements of small ice circular cylinders depend strongly on wavelengths and refractive indices, particle size distributions, and the L/D ratios. The applicability of the power-law distribution and particle shapes for light scattering calculations for small ice crystals is discussed. The effects of the effective variance of size distribution on light scattering characteristics are addressed. It seems from the behavior of scattering matrix elements of small ice crystals that the combination of 25 and 3.979 microm has some advantages and potential applications for remote sensing of cirrus and other ice clouds.     

Scattering properties of needlelike and platelike ice spheroids with moderate size parameters

We use the current advanced version of the T-matrix method to compute the optical cross sections, the asymmetry parameter of the phase function, and the scattering matrix elements of ice spheroids with aspect ratios up to 20 and surface-equivalent-sphere size parameters up to 12. We demonstrate that platelike and needlelike particles with moderate size parameters possess unique scattering properties: their asymmetry parameters and phase functions are similar to those of surface-equivalent spheres, whereas all other elements of the scattering matrix are typical of particles much smaller than the wavelength (Rayleigh scatterers). This result may have important implications for optical particle sizing and remote sensing of the terrestrial and planetary atmospheres.     

Laboratory measurements and T-matrix calculations of the scattering matrix of rutile particles in water

We present experimentally determined scattering matrix elements of birefringent rutile particles in water as a function of the scattering angle for a wavelength of 633 nm (in air). These elements are compared with the results of T-matrix calculations for prolate spheroids. For the diagonal matrix elements the results of the T-matrix calculations are in good agreement with those of the measurements. A good fit for the whole matrix, including the off-diagonal elements, is obtained when we compensate for the birefringence of the rutile particles by performing the computations for spheroids with a slightly larger length/width ratio than measured.     

Improved T-matrix computations for large, nonabsorbing and weakly absorbing nonspherical particles and comparison with geometrical-optics approximation

We show that the use of a matrix inversion scheme based on a special lower triangular-upper triangular factorization rather than on the standard Gaussian elimination significantly improves the numerical stability of T-matrix computations for nonabsorbing and weakly absorbing nonspherical particles. As a result, the maximum convergent size parameter for particles with small or zero absorption can increase by a factor of several and can exceed 100. We describe an improved scheme for evaluating Clebsch-Gordon coefficients with large quantum numbers, which allowed us to extend the analytical orientational averaging method developed by Mishchenko [J. Opt. Soc. Am. A 8, 871 (1991)] to larger size parameters. Comparisons of T-matrix and geometrical optics computations for large, randomly oriented spheroids and finite circular cylinders show that the applicability range of the ray-tracing approximation depends on the imaginary part of the refractive index and is different for different elements of the scattering matrix.     

The T-matrix revisited

We consider electromagnetic scattering from penetrable cylinders of general cross section. After summarizing the basic T-matrix equations the low-frequency case is examined, which leads for nonmagnetic materials to the exact result T=iR-R(2) in the Rayleigh limit, satisfying both reciprocity and energy constraints. Here elements of R are given by integrals of regular wave functions over the cylinder surface. A "Rayleigh expansion" is then found that is convergent throughout the Rayleigh region and the lower end of the resonance region and requires no matrix inversion. For bodies of high aspect ratio, there is a problem with significance loss during numerical integration, due to large oscillatory terms. A class of surfaces has now been found for which these terms can be removed, however, enabling us to treat aspect ratios up to 1,000:1. These methods are expected to apply also in three dimensions.     

Use of circular cylinders as surrogates for hexagonal pristine ice crystals in scattering calculations at infrared wavelengths

We investigate the errors associated with the use of circular cylinders as surrogates for hexagonal columns in computing the optical properties of pristine ice crystals at infrared (8-12-microm) wavelengths. The equivalent circular cylinders are specified in terms of volume (V), projected area (A), and volume-to-area ratio that are equal to those of the hexagonal columns. We use the T-matrix method to compute the optical properties of the equivalent circular cylinders. We apply the finite-difference time-domain method to compute the optical properties of hexagonal ice columns smaller than 40 microm. For hexagonal columns larger than 40 microm we employ an improved geometric optics method and a stretched scattering potential technique developed in previous studies to calculate the phase function and the extinction (or absorption) efficiency, respectively. The differences between the results for circular cylinders and hexagonal columns are of the order of a few percent. Thus it is quite reasonable to use a circular cylinder geometry as a surrogate for pristine hexagonal ice columns for scattering calculations at infrared (8-12-microm) wavelengths. Although the pristine ice crystals can be approximated as circular cylinders in scattering calculations at infrared wavelengths, it is shown that optical properties of individual aggregates cannot be well approximated by those of individual finite columns or cylinders.     

Light scattering at oblique incidence on two coaxial cylinders

A solution for the problem of a plane wave at oblique incidence on two coaxial cylinders is presented. The solution of the wave equation is determined for various geometric regions, and boundary conditions are applied at the material interfaces. The resulting solution consists of a system of eight equations in eight unknown coefficients. Expressions for two of the Mueller-scattering matrix elements (S(11) and S(12)) and the extinction, scattering, and backscattering cross sections are derived. A numerical algorithm for the solution is developed and implemented. The algorithm is tested for several limiting cases: homogeneous, hollow, and metal-core cylinders at various angles of incidence for TM and TE waves. Comparisons of the results of the algorithm with the results of studies reported in the literature are made. The comparisons are favorable, achieving good agreement with published work. For two coaxial cylinders, the numerical calculations show that if one is to use light scattering as a diagnostic tool, both of the Mueller-scattering matrix elements S(11) and S(12) must be measured simultaneously. In addition, the backscattering cross section is very sensitive for monitoring change in the radii of the cylinders.     

Application of the exact solution for scattering by an infinite cylinder to the estimation of scattering by a finite cylinder

A new algorithm for cylindrical Bessel functions that is similar to the one for spherical Bessel functions allows us to compute scattering functions for infinitely long cylinders covering sizes ka = 2πa/λ up to 8000 through the use of only an eight-digit single-precision machine computation. The scattering function and complex extinction coefficient of a finite cylinder that is seen near perpendicular incidence are derived from those of an infinitely long cylinder by the use of Huygens's principle. The result, which contains no arbitrary normalization factor, agrees quite well with analog microwave measurements of both extinction and scattering for such cylinders, even for an aspect ratio p = l/(2a) as low as 2. Rainbows produced by cylinders are similar to those for spherical drops but are brighter and have a lower contrast.     

Light extinction and absorption by arbitrarily oriented finite circular cylinders by use of geometrical path statistics of rays

From the geometrical path statistics of rays in an anomalous-diffraction theory (ADT) [Opt. Lett. 28, 179 (2003)] closed-form expressions for the geometrical path distribution of rays and analytical formulas for the optical efficiencies of finite circular cylinders oriented in an arbitrary direction with respect to the incident light are derived. The characteristics of the shapes of the cylinders produce unique features in the geometrical path distributions of the cylinders compared with spheroids. Gaussian ray approximations, which depend only on the mean and the mean-squared geometrical paths of rays, of the optical efficiencies of finite circular cylinders and spheroids are compared with the exact optical efficiencies in ADT. The influence of the difference in shape between cylinders and spheroids on the optical efficiencies in ADT is illustrated by their respective geometrical path distributions of rays.     

Scattering of light from arbitrarily oriented finite cylinders

An iterative approach to the scattering of light from a finite dielectric cylinder first developed by Shifrin and extended by Acquista is applied to cases where the phase shift is <2, and the cylinder is arbitrarily oriented. It is found that the first 2 orders of the iteration converge to within 1% when the aspect ratio (length/diameter) of the cylinder is as small as 20. The results are compared to the exact theory for infinite cylinders, and the effects of finite size are calculated and discussed.     

Modeling the tensile stress–strain response of carbon nanotube/polypropylene nanocomposites using nonlinear representative volume element

This paper presents a finite element model for predicting the mechanical behavior of polypropylene (PP) composites reinforced with carbon nanotubes (CNTs) at large deformation scale. Existing numerical models cannot predict composite behavior at large strains due to using simplified material properties and inefficient interfaces between CNT and polymer. In this work, nonlinear representative volume elements (RVE) of composite are prepared. These RVEs consist of CNT, PP matrix and non-bonded interface. The nonlinear material properties for CNT and polymer are adopted to solid elements. For the first time, the interface between CNT and matrix is simulated using contact elements. This interfacial model is capable enough to simulate wide range of interactions between CNT and polymer in large strains. The influence of adding CNT with different aspect ratio into PP is studied. The mechanical behavior of composites with different interfacial shear strength (ISS) is discussed. The success of this new model was verified by comparing the simulation results for RVEs with conducted experimental results. The results shows that the length of CNT and ISS values significantly affect the reinforcement phenomenon.     

Calculation of the single-scattering properties of randomly oriented hexagonal ice columns: a comparison of the T-matrix and the finite-difference time-domain methods

We calculated the scattering and absorption properties of randomly oriented hexagonal ice columns using T-matrix theory, employing analytic orientation averaging, and the finite-difference time-domain method, which uses a numerical procedure to simulate random orientation. The total optical properties calculated are the extinction efficiency, absorption efficiency, single-scattering albedo, and the asymmetry parameter. The optical properties are calculated at the wavelengths of 0.66, 8.5, and 12 mum, up to a size parameter of 20 at 0.66 mum and 15 at the two other wavelengths. The phase-matrix elements P11, P12, and P22 are also calculated and compared, up to a size parameter of 20 at 0.66 mum and 15 at 12.0 mum. The scattering and absorption solutions obtained from the two independent electromagnetic methods are compared and contrasted, as well as the central processing unit time and memory load for each size parameter. It is found that the total optical properties calculated by the two methods are well within 3% of each other for all three wavelengths and size parameters. In terms of the phase-matrix elements it is found that there are some differences between the T-matrix and the finite-difference time-domain methods appearing in all three elements. Differences between the two methods for the P11 element are seen particularly at scattering angles from approximately 120 degrees to 180 degrees ; and at the scattering angle of 180 degrees , relative differences are less than 16%. At scattering angles less than 100 degrees , agreement is generally within a few percent. Similar results are also found for the P12 and P22 elements of the phase matrix. The validity of approximating randomly oriented hexagonal ice columns by randomly oriented equal surface area circular cylinders is also investigated in terms of the linear depolarization ratio.     

Extension of the discrete sources method to light scattering by highly elongated finite cylinders

Abstract

Light scattering by long finite cylinders is applied for particle characterization, investigation of scattering and absorption properties of interstellar dust, ice crystals and many other fields. In recent years many methods have been used to solve this problem, but usually their applicability is restricted to aspect ratios of about 10. In this paper a renewed algorithm of the discrete sources method is described, which allows computer simulation of light scattering by highly elongated cylinders with aspect ratios up to 100 and length up to 40 μm.     

Light scattering simulation by oblate disc spheres using the null field method with discrete sources located in the complex plane

The T-matrix method, which is also known as the null field method (NFM) or extended boundary condition method (EBCM), has established itself as a well known and highly regarded method for calculating light scattering by non-spherical particles. Its biggest advantage is the possibility to obtain all information about the scattering characteristics of the particle and to store it into one matrix. This enables one to do additional investigations with low efforts. Unfortunately the standard NFM fails to converge for particles with extremely non-spherical particle shapes, like long cylinders or coin-like flat cylinders. In this paper we investigate light scattering by finite particles in the form of an oblate disc sphere, which can be described as flat cylinders with a rounded edge. We use an advanced form of the T-matrix method—the null field method with discrete sources (NFM-DS). By presenting light scattering results we would like to demonstrate the potential this advanced NFM-DS offers. It allows one to calculate particle shapes with aspect ratios (relation between radius and thickness of the particle) up to 100:1 and size parameters (relation between radius and wavelength) up to 30.     

Light scattering by a core-mantle spheroidal particle

A solution of the electromagnetic scattering problem for confocal coated spheroids has been obtained by the method of separation of variables in a spheroidal coordinate system. The main features of the solution are (i) the incident, scattered, and internal radiation fields are divided into two parts: an axisymmetric part independent of the azimuthal angle ? and a nonaxisymmetric part that with integration over ? gives zero; the diffraction problems for each part are solved separately; (ii) the scalar potentials of the solution are chosen in a special way: Abraham's potentials (for the axisymmetric part) and a superposition of the potentials used for spheres and infinitely long cylinders (for the nonaxisymmetric part). Such a procedure has been applied to homogeneous spheroids [Differential Equations 19, 1765 (1983); Astrophys. Space Sci. 204, 19, (1993)] and allows us to solve the light scattering problem for confocal spheroids with an arbitrary refractive index, size, and shape of the core or mantle. Numerical tests are described in detail. The efficiency factors have been calculated for prolate and oblate spheroids with refractive indices of 1.5 + 0.0 i, 1.5 + 0.05 i for the core and refractive indices of 1.3 + 0.0 i, 1.3 + 0.05i for the mantle. The effects of the core size and particle shape as well as those of absorption in the core or mantle are examined. It is found that the efficiency factors of the coated and homogeneous spheroids with the volume-averaged refractive index are similar to first maximum.     

Sensitivity of the backscattering Mueller matrix to particle shape and thermodynamic phase

The Mueller matrix (M) corresponding to the phase matrix in the backscattering region (scattering angles ranging from 175 degrees to 180 degrees) is investigated for light scattering at a 0.532-microm wavelength by hexagonal ice crystals, ice spheres, and water droplets. For hexagonal ice crystals we assume three aspect ratios (plates, compact columns, and columns). It is shown that the contour patterns of the backscattering Mueller matrix elements other than M11, M44, M14, and M41 depend on particle geometry; M22 and M33 are particularly sensitive to the aspect ratio of ice crystals. The Mueller matrix for spherical ice particles is different from those for nonspherical ice particles. In addition to discriminating between spherical and nonspherical particles, the Mueller matrix may offer some insight as to cloud thermodynamic phase. The contour patterns for large ice spheres with an effective size of 100 microm are substantially different from those associated with small water droplets with an effective size of 4 microm.     

Light scattering by polydispersions of randomly oriented spheroids with sizes comparable to wavelengths of observation

We report the results of an extensive study of the scattering of light by size and size-shape distributions of randomly oriented prolate and oblate spheroids with the index of refraction 1.5 + 0.02i typical of some mineral terrestrial aerosols. The scattering calculations have been carried out with Waterman's T-matrix approach, as developed recently by Mishchenko [J. Opt. Soc. Am. A 8, 871 (1991); Appl. Opt. 32, 4562 (1993)]. Our main interest is in light scattering by polydisperse models of nonspherical particles because averaging over sizes provides more realistic modeling of natural ensembles of scattering particles and washes out the interference structure and ripple typical of monodisperse scattering patterns, thus enabling us to derive meaningful conclusions about the effects of particle nonsphericity on light scattering. Following Hansen and Travis [Space Sci. Rev. 16, 527 (1974)], we show that scattering properties of most physically plausible size distributions of randomly orientednonspherical part cles depend primarily on the effective equivalent-sphere radius and effective variance of the distribution, the actual shape of the distribution having a minor influence. To minimize the computational burden, we have adopted a computationally convenient power law distribution of particle equivalent-sphere radii n(r) α r(-3),r(1) ≤ r≤r(2). The effective variance of the size distribution is fixed at 0.1, and the effective size parameter continuously varies from 0 to 15. We present results of computer calculations for 24 prolate and oblate spheroidal shapes with aspect ratios from 1.1 to 2.2. The elements of the scattering matrix for the whole range of size parameters and scattering angles are displayed in the form of contour plots. Computational results are compared with analogous calculations for surface-equivalent spheres, and the effects of particle shape on light scattering are discussed in detail.     

水下刚硬体声散射场计算及分析

针对水下任意形状刚硬体,建立了其声散射T矩阵法计算模型。利用物体的几何对称性,提出了各种简化计算Q矩阵的方法。以两端有半椭球顶的有限长圆柱体为例,分别计算了具有不同长宽比或者是不同端部几何形状的有限长圆柱体的反向散射和前向散射函数。深入讨论了有限长圆柱体的镜反射波和Franz波的干涉物理模型,以及其圆柱体中段长度对散射场的影响,并对得到的计算结果给出了合理的物理解释。计算结果同时也表明T矩阵法是一种很简便实用的计算水下物体远场散射的方法。