Light scattering by a multilayered spheroidal particle

The light scattering problem for a confocal multilayered spheroid has been solved by the extended boundary condition method with a corresponding spheroidal basis. The solution preserves the advantages of the approach applied previously to homogeneous and core-mantle spheroids, i.e., the separation of the radiation fields into two parts and a special choice of scalar potentials for each of the parts. The method is known to be useful in a wide range of the particle parameters. It is particularly efficient for strongly prolate and oblate spheroids. Numerical tests are described. Illustrative calculations have shown that the extinction factors converge to average values with a growing number of layers and how the extinction varies with a growth of particle porosity.     

Effective medium theories for irregular fluffy structures: aggregation of small particles

The extinction efficiencies as well as the scattering properties of particles of different porosity are studied. Calculations are performed for porous pseudospheres with small size (Rayleigh) inclusions using the discrete dipole approximation. Five refractive indices of materials covering the range from 1.20+0.00i to 1.75+0.58i were selected. They correspond to biological particles, dirty ice, silicate, and amorphous carbon and soot in the visual part of the spectrum. We attempt to describe the optical properties of such particles using Lorenz-Mie theory and a refractive index found from some effective medium theory (EMT) assuming the particle is homogeneous. We refer to this as the effective model. It is found that the deviations are minimal when utilizing the EMT based on the Bruggeman mixing rule. Usually the deviations in the extinction factor do not exceed approximately 5% for particle porosity P = 0 - 0.9 and size parameters x(porous) = 2 pi r(s,porous)/lambda < or approximately = 25. The deviations are larger for scattering and absorption efficiencies and smaller for particle albedo and the asymmetry parameter. Our calculations made for spheroids confirm these conclusions. Preliminary consideration shows that the effective model represents the intensity and polarization of radiation scattered by fluffy aggregates quite well. Thus the effective models of spherical and nonspherical particles can be used to significantly simplify the computations of the optical properties of aggregates containing only Rayleigh inclusions.     

Light scattering by a coated sphere illuminated with a Gaussian beam

The intensity of light scattered by a coated sphere illuminated with an off-axis Gaussian beam is calculated. Results are shown for different beam positions with respect to the sphere. As the beam is shifted further away from the surface of the sphere, the higher-Q morphology-dependent resonances become increasingly important in the backscatter spectra, and the angular scattering intensity becomes smoother.

The scattered intensity depends on the beam position, the refractive indices of the core and coat, the radius of the core, and the thickness of the coat. As the beam is moved further away from the sphere, the effect of the core on the scattering intensity decreases. When the incident Gaussian beam is focused outside of a particle with a relatively small core, the scattering spectra and angular scattering patterns become similar to those of a homogeneous sphere having the refractive index of the coat. These calculated results suggest that measurements of spectral scattering and angular scattering patterns for several Gaussian beam positions could be useful for the characterization of coated spheres.

     

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.     

Finite-difference time-domain solution of light scattering by dielectric particles with large complex refractive indices

The finite-difference time-domain (FDTD) technique is examined for its suitability for studying light scattering by highly refractive dielectric particles. It is found that, for particles with large complex refractive indices, the FDTD solution of light scattering is sensitive to the numerical treatments associated with the particle boundaries. Herein, appropriate treatments of the particle boundaries and related electric fields in the frequency domain are introduced and examined to improve the accuracy of the FDTD solutions. As a result, it is shown that, for a large complex refractive index of 7.1499 + 2.914i for particles with size parameters smaller than 6, the errors in extinction and absorption efficiencies from the FDTD method are generally less than ~4%. The errors in the scattering phase function are less than ~5%. We conclude that the present FDTD scheme with appropriate boundary treatments can provide a reliable solution for light scattering by nonspherical particles with large complex refractive indices.     

Light scattering by a coated infinite cylinder in an absorbing medium

The scattering formulation for a coated infinite cylinder in an absorbing medium is presented in this paper. The cylinder is subjected to an arbitrarily polarized plane wave propagating in a general direction at the cylinder. The refractive index and magnetic permeability of the host medium, as well as those for the core and coating of the cylinder, can be real or complex. The scattering and extinction efficiencies and the scattering amplitudes are derived for both the near field and the far field. As the medium is absorbing, the "true" extinction and scattering efficiencies are derived based on the radiative energy outflow at the surface of the cylinder. The radiative efficiencies in the far field are denoted as "apparent" properties because they include absorption by the intervening medium. The influence of the refractive index and permeability of the host medium on the scattering properties of a coated cylinder is illustrated by numerical examples.     

Scattering of light by bispheres with touching and separated components

We use the T-matrix method as described by Mishchenko and Mackowski [Opt. Lett. 19, 1604 (1994)] to compute light scattering by bispheres in fixed and random orientations extensively. For all our computations the index of refraction is fixed at a value 1.5 + 0.005i, which is close to the refractive index of mineral tropospheric aerosols and was used in previous extensive studies of light scattering by spheroids and Chebyshev particles. For monodisperse bispheres with touching components in a fixed orientation, electromagnetic interactions between the constituent spheres result in a considerably more complicated interference structure in the scattering patterns than that for single monodisperse spheres. However, this increased structure is largely washed out by orientational averaging and results in scattering patterns for randomly oriented bispheres that are close to those for single spheres with size equal to the size of the bisphere components. Unlike other nonspherical particles such as cubes and spheroids, randomly oriented bispheres do not exhibit pronounced enhancement of side scattering and reduction of backscattering and positive polarization at side-scattering angles. Thus the dominant feature of light scattering by randomly oriented bispheres is the single scattering from the component spheres, whereas the effects of cooperative scattering and concavity of the bisphere shape play a minor role. The only distinct manifestations of nonsphericity and cooperative scattering effects for randomly oriented bispheres are the departure of the ratio F(22)/F(11) of the elements of the scattering matrix from unity, the inequality of the ratios F(33)/F(11) and F(44)/F(11), and nonzero linear and circular backscattering depolarization ratios. Our computations for randomly oriented bispheres with separated wavelengthsized components show that the component spheres become essentially independent scatterers at as small a distance between their centers as 4 times their radii.     

Light scattering by ellipsoids in a physical optics approximation

A physical optics approximation based on Presnel's laws is developed to calculate the intensity of light scattered by a three-axis ellipsoid of any orientation and any refractive index. Some results concerning totally reflecting spheres and dielectric spheroids are presented. An approach suitable for large scatterers is particularly good for small scattering angles. The angular intensities, i(1) and i(2), are then plotted versus θ for large axially oriented ellipsoids of various thicknesses. Theoretical small-angle light-scattering patterns are also presented and discussed. The data from one of them correspond to red cells in a shear flow.     

Feasibility of calculating morphology-dependent resonance peaks by numerical methods

We investigate the feasibility of numerically calculating morphology-dependent resonance (MDR) peaks. To do so, one has to calculate the scattering intensities numerically and determine how difficult it is to numerically predict the position and the magnitude of the MDR peaks. However, at present, in practice it is impossible to calculate MDR peaks with a personal computer because so much computing time is required. Therefore the surface values of the Debye potential and its derivative for a homogeneous sphere are obtained from Mie's analytical solution and then used in integral equations to give the scattering intensities at a specific position of infinity by numerical integrations. It is shown that if a sufficient number of surface elements are used, the MDR peaks are exactly calculated for a homogeneous sphere with refractive index of 1.5, 1.4, and 1.3 up to a size parameter of 20. One can conjecture the number of finite and boundary elements necessary to numerically compute accurate scattering intensities. It should be also noted that the number of surface elements necessary for exact integration shows peaks similar to MDR peaks with respect to the size parameter. Therefore one will need many more elements at the size parameter at which the MDR occurs.     

Electromagnetic Resonance at the Interface of Two Concentric Spheres

Abstract

Some results of a theoretical study on the interaction of an electromagnetic plane wave with a coated sphere are presented. The finely stratified sphere model (Kai, L., and Massoli, P. 1994, Appl. Optics, 33 501) is used for numerical calculations. It is found that when the refractive index of the core is larger than the refractive index of the coating, the Q _ext-x spectrum for a coated sphere differs from the spectrum of a homogeneous sphere by having both sharp minima and sharp maxima, all corresponding to strong electromagnetic resonances at the interface of the core and the coating.     

The coating of cadmium compounds with nickel compounds

Micrometer-sized particles of cadmium carbonate (cores) were prepared by homogeneous precipitation from aqueous solution containing urea in the presence of cadmium sulfate. These particles were then homogeneously coated with a layer of nickel hydroxy carbonate by heating their dispersion in aqueous solution, containing nickel sulfate and urea, at 85 °C for 70 min with constant agitation. The same solution mixture produced spheroids of nickel hydroxy carbonate (coating precursors), when heated under similar conditions in the absence of the cadmium carbonate particles. The existence of the coated layer on the cores and its composition was confirmed by various physical methods. The as-prepared carbonated solids (cores, coating precursors, and coated particles) were converted into their oxide forms by calcination at 700 °C. The cores became porous, whereas the coating precursors and coating layer disintegrated into smaller particles during the calcination process.     

Simultaneous retrieval of the complex refractive indices of the core and shell of coated aerosol particles from extinction measurements using simulated annealing

Simultaneously retrieving the complex refractive indices of the core and shell of coated aerosol particles given the measured extinction efficiency as a function of particle dimensions (core diameter and coated diameter) is much more difficult than retrieving the complex refractive index of homogeneous aerosol particles. Not only must the minimization be performed over a four-parameter space, making it less efficient, but in addition the absolute value of the difference between the measured extinction and the calculated extinction does not have an easily distinguished global minimum. Rather, there are a number of local minima to which almost all conventional retrieval algorithms converge. In this work, we develop a new (to our knowledge) retrieval algorithm that employs the numerical method known as simulated annealing with an innovative "temperature" schedule. This study is limited only to spherical particles with a concentric shell and to cases in which the diameter of both the core and the coated particle are known. We find that when the top ranking particle sizes according to their information content are combined from separate experiments to make up the particle size distribution, the simulated annealing retrieval algorithm is quite robust and by far superior to a greedy random perturbation approach often used.     

Single scattering of light by circular cylinders

We consider two topics pertaining to light scattering by circular cylinders. (A) Scattering properties of cylinders with increasing aspect ratio are studied. It is shown that the solution for finite cylinders does not converge to the solution for infinitely long cylinders if the aspect ratio increases. This is due to differences in the treatment of diffraction for finite and infinite cylinders. (B) Finite cylinders have sharp edges, so their scattering properties differ from those of spheroids having the same aspect ratio. To illustrate these differences we present scattering matrix elements of cylinders and spheroids for a large set of aspect ratios. To handle the large amount of data, the scattering matrix elements as functions of aspect ratio and scattering angle are presented in so-called three-dimensional figures.     

Revisit of an elliptic hole problem by using complex variable method

In this paper, the problem of an elliptic hole embedded in an infinite plate interacting with an arbitrary point load is revisited by using the complex variable method. Based on analytical continuation theorem, the continuity conditions across the interface are automatically satisfied in a straightforward manner. It is shown that the solution for an infinite domain with an elliptic hole can be obtained from the solution of the corresponding homogeneous problem merely by a simple algebraic expression. This relation is universal in the sense of being independent of the loading considered. The solution of the corresponding homogeneous problem is considered as the principal part of the complex potentials while the complementary part of the complex potentials can be obtained by using analytical continuation theorem. Different expressions of the complementary part of the complex potentials are presented in this paper which are all proved to be the same result.     

Finite element analysis of nonlinear orthotropic hyperelastic membranes

A finite element method is presented for geometrically and materially nonlinear orthotropic hyperelastic membranes. The constitutive relations are formulated in terms of the invariants of the 2D right Cauchy-Green strain tensor and the resulting system of nonlinear equations solved using a Newton-Raphson approach. Both axisymmetric and nonaxisymmetric versions of the method are developed and validated. Example problems are solved for isotropic and orthotropic membranes, and the effect of various parameters investigated. Finally, convergence studies are performed for various degrees of anisotropy.Also supported in part by the Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant F49620-95-1-0230.     

Retrieval of the scattering and microphysical properties of aerosols from ground-based optical measurements including polarization. I. Method

A method has been developed for retrieving the scattering and microphysical properties of atmospheric aerosol from measurements of solar transmission, aureole, and angular distribution of the scattered and polarized sky light in the solar principal plane. Numerical simulations of measurements have been used to investigate the feasibility of the method and to test the algorithm's performance. It is shown that the absorption and scattering properties of an aerosol, i.e., the single-scattering albedo, the phase function, and the polarization for single scattering of incident unpolarized light, can be obtained by use of radiative transfer calculations to correct the values of scattered radiance and polarized radiance for multiple scattering, Rayleigh scattering, and the influence of ground. The method requires only measurement of the aerosol's optical thickness and an estimate of the ground's reflectance and does not need any specific assumption about properties of the aerosol. The accuracy of the retrieved phase function and polarization of the aerosols is examined at near-infrared wavelengths (e.g., 0.870 mum). The aerosol's microphysical properties (size distribution and complex refractive index) are derived in a second step. The real part of the refractive index is a strong function of the polarization, whereas the imaginary part is strongly dependent on the sky's radiance and the retrieved single-scattering albedo. It is demonstrated that inclusion of polarization data yields the real part of the refractive index.     

Scattering by particles with radially variable refractive indices

The problem of scattering by a spherical particle whose refractive index arbitrarily depends on the distance from its center has been solved. A computational scheme for determining the scattering coefficients for the refractive-index profiles given by some piecewise smooth function is constructed, and the models to use it are presented. The simple algorithm for evaluating the scattered and internal field vectors is elaborated. The exact expression for the scattering cross section in terms of the generated Debye partial potentials that are dependent on the given refractive index has been obtained.     

Light scattering by aggregated red blood cells

In low flow rates, red blood cells (RBCs) fasten together along their axis of symmetry and form a so-called rouleaux. The scattering of He-Ne laser light by a rouleau consisting of n (2 < or = n < or = 8) average-sized RBCs is investigated. The interaction problem is treated numerically by means of an advanced axisymmetric boundary element--fast Fourier transform methodology. The scattering problem of one RBC was solved first, and the results showed that the influence of the RBC's membrane on the scattering patterns is negligible. Thus the rouleau is modeled as an axisymmetric, homogeneous, low-contrast dielectric cylinder, on the surface of which appears, owing to aggregated RBCs, a periodic roughness along the direction of symmetry. The direction of the incident laser light is considered to be perpendicular to the scatterer's axis of symmetry. The differential scattering cross sections in both perpendicular and parallel scattering planes and for all the scattering angles are calculated and presented in detail.     

Scattering of he-ne laser light by an average-sized red blood cell

The scattering of He-Ne laser light by an average-sized human red blood cell (RBC) is investigated numerically. The RBC is modeled as an axisymmetric, low-contrast dielectric, biconcave disk. The interaction problem is treated numerically by means of a boundary-element methodology. The differential scattering cross sections (DSCS's) corresponding to various cell orientations are calculated. The numerical results obtained for the exact RBC geometry are compared with those corresponding to a scattering problem in which the cell is assumed to be either a volume-equivalent sphere or an oblate spheroid. A parametric study demonstrating the dependence of the DSCS on the wavelength of the incident wave and the cell's refractive index is presented.     

确定轴向冲击下薄壁圆柱壳第二临界速度的新方法

研究薄壁圆柱壳的动态屈曲模式，有助于构造具有高吸能率的抗冲击结构。根据轴向冲击下的薄壁圆柱壳存在使其屈曲模式由轴对称转换为非轴对称的第二临界速度，且当冲击速度大于第二临界速度时薄壁圆柱壳的屈曲模式先呈现轴对称形式，然后随着冲击响应时间逐渐由轴对称形式转化为非轴对称形式这一理论，基于有限元仿真，比较撞击系统动能的时间历程和屈曲变形的时间历程，提出了用以确定第二临界速度的能量迭代法。应用此方法设计薄壁圆柱壳的动力屈曲结构可有效地减少试验次数，降低实验成本。该方法的可行性和正确性利用落锤实验得到了验证。