一层随机非球形粒子高阶散射矢量辐射传输迭代解的参数求逆

根据一层随机非球形粒子矢量辐射传输(VRT)方程的一阶与高阶Mueller矩阵解的差异,提出了在已知一些参数的条件下,用水平与垂直同极化和交叉极化后向散射测量反演非球形粒子复介电常数和单位面积粒子数的迭代方法,并用极化后向散射强度的模拟测量数据及机载合成孔径雷达(Air SAR)实测数据进行了反演试验,并对反演值与各参数的关系以及测量噪声对反演的影响进行了讨论.     

粗糙面上粒子层矢量辐射传输方程的二阶迭代解法

为了探究典型粗糙面上随机粒子层中能量传输的多次散射机制,提出了一种基于矢量辐射传输方程的建模二阶计算方法.该建模方法将建模场景(粗糙面上粒子层)在高度维(Z轴)划分为多个传输散射层,基于矢量辐射传输理论中的一阶迭代散射解,利用典型粗糙面的半经验半解析方法,求解出整个场景的二阶迭代散射解.同时,研究粒子层能量在粒子与粒子间的多次散射机制,以及粒子与地表粗糙面间的多次散射机制.数值结果表明,该二阶迭代解法相较于矢量辐射传输方程的一阶迭代散射解,能够更完整地探究互作用的散射特性,且可从能量传输角度解译建模场景中物体间的相干作用,从而可用于植被地物环境下的多次散射机制的解析以及散射系数变化趋势的预估.     

非均匀随机散射介质层的多次散射和热辐射的分层迭代解

为纳入矢量辐射传输方程的多次散射，将非均匀平行分层分布的随机散射介质划分成多个薄 层，利用各薄层的低阶Mueller矩阵解推导出整层的高阶散射迭代解.用这一方法计算了具有 占空比、粒子大小和温度廓线非均匀分布的平行分层随机密集球形粒子积雪层的极化热辐射 .数值结果与一层均匀随机散射层的离散坐标-特征值方法作了比较，讨论了多次散射和热辐 射与散射层各特征参数的关系，与积雪层微波遥感数据进行了对比. 关键词： VRT方程 平行分层 高阶散射解 迭代解     

非均匀散射层矢量辐射传输(VRT)方程高阶散射解的迭代法

将散射介质层在z轴方向划分成薄层，用薄层的一阶散射强度、Fourier变换和迭代方法求解散射介质整层的矢量辐射传输(VRT)方程的高阶散射解.该方法将一阶散射与高阶散射迭代结合起来，计算公式简明，可计算高阶迭代解，计算时间少.计算结果与一层均匀散射介质的VRT方程一阶Mueller矩阵解、半空间均匀散射介质二阶Mueller矩阵解、以及离散坐标-特征值特征矢量法的VRT热辐射的数值解作了全面的比较.提出并讨论了非均匀散射层主动与被动VRT方程的高阶解.本计算程序可以通用于非球形粒子多层结构及非均匀介质的散射和热辐射计算. 关键词： VRT方程 分层 迭代解     

非球形气溶胶粒子散射相函数经验公式

散射相函数是研究电磁波传输特性的重要参数，直接影响电磁波传输方程的简化程度和解的精度。基于电磁散射与辐射传输中的基本理论，对非球形粒子散射相函数的经验公式进行了研究。为了很好的模拟非球形粒子的后向散射峰值，提高辐射传输方程的简化程度和解的精度，提出了一种新的相函数经验公式。分析新的相函数对非球形粒子的适用性，以单个沙尘性气溶胶为例，计算了不同形状粒子的Henyey-Greenstein*相函数和新的相函数随角度的变化，并与T矩阵法的计算结果进行了对比，发现椭球形粒子的长短轴比和有限长圆柱形粒子的径长比大于0.5时，新的相函数在大角度后向散射部分与T矩阵法的吻合程度较高。考虑波长变化，对比了尺寸谱满足对数正态分布的四种气溶胶粒子的Henyey-Greenstein*相函数和新的相函数与T矩阵法的计算结果。研究表明，对于椭球形粒子和有限长圆柱形粒子，在大角度(大于90°)后向散射部分，除了0.694时的椭球形海洋性气溶胶，新的相函数均方根差较小的占100%，证明了新的相函数可以较好的模拟非球形粒子的后向散射特征。新的相函数对准确模拟辐射传输过程具有重要意义。     

紫外波段多分散系气溶胶散射相函数随机抽样方法研究

基于电磁散射与辐射传输中的基本理论,对紫外波段霾尺度范围内满足特定分布的多种气溶胶粒子的散射相函数进行了研究.提出了一种直接随机抽样拟合散射相函数的方法.比较了H-G相函数、改进的H-G相函数及随机抽样拟合的相函数与多分散系Mie相函数的偏离程度.数值计算了不同相函数拟合方法对应气溶胶的传输特性.计算结果表明,相函数的准确模拟计算对于蒙特卡罗方法等辐射传输问题的解决具有十分重要的意义.     

利用T-matrix计算非球形粒子散射特性的研究

本文用T-matrix方法计算了非球形气溶胶粒子的光学特性,得到了气溶胶粒子的消光截面、散射截面、吸收截面与气溶胶粒子形状的关系,不同形状气溶胶粒子的有着相同的散射相函数和不同的偏振度,非球形气溶胶粒子的散射相函数对其复折射指数的实部和虚部都不太敏感,而偏振相函数对其实部和虚部都比较敏感.此结论为研究大气辐射传输提供了较好的方法,尤其是偏振度与偏振相函数的提出为用偏振的方法进一步的反演气溶胶的光学参数提供了理论基础.     

改进的单次散射相函数解析表达式

单次散射相函数对电磁辐射传输模拟过程的准确性和计算效率有重要的影响.基于电磁散射与辐射传输中的基本理论,对单次散射相函数的解析表达式进行了研究,提出了一种新的单次散射相函数解析表达式.比较了单个粒子的Henyey-Greenstein相函数、Henyey-Greenstein*相函数与新的相函数随角度的分布,发现新的散射相函数提高了后向散射峰值,可以更合理地描述单个粒子的散射特性.按三种气溶胶粒子谱分布模式计算了Henyey-Greenstein*相函数和新的相函数对应的数值结果,并与多分散系Mie散射相函数进行对比,发现新的相函数提高了与多分散系Mie散射相函数的符合程度.研究表明,对于大角度(大于90°)后向散射,新的相函数与Mie散射相函数均方根差较小的占73.3%,高于Henyey-Greenstein*相函数的26.7%,证明了新的相函数可以显著提高后向散射峰值.新的相函数对准确模拟辐射传输过程具有重要意义.     

随机介质中的矢量辐射传输理论

本文讨论了极化电磁波在随机介质中多次散射,传输和热辐射的斯托克斯矢量的辐射传输理论。其中包括随机分布离散的球形和非球形粒子的矢量辐射传输方程,离散坐标-特征分析法,付利叶变换,迭代法等数值解。讨论了非球形粒子的穆勒矩阵。并研究了密集分布的散射粒子介质的辐射传输理论,考虑了密集粒子散射的相干性,计算了有效传播常数。理论及数值结果与实验作了很好的比较。     

随机介质中的矢量辐射传输理论

本文讨论了极化电磁波在随机介质中多次散射，传输和热辐射的斯托克斯矢量的辐射传输理论。其中包括随机分布离散的球形和非球形粒子的矢量辐射传输方程，离散坐标——特征分析法，付利叶变换，迭代法等数值解。讨论了非球形粒子的穆勒矩阵。并研究了密集分布的散射粒子介质的辐射传输理论，考虑了密集粒子散射的相干性，计算了有效传播常数。理论及数值结果与实验用了很好的比较。     

On scattering of light by small axially symmetric particles

Light scattering by small dielectric particles of an arbitrary axially symmetric shape is analyzed. A simple approximate expression that governs the polarizability of the particle is found under the assumption of field homogeneity inside of these particles. The expression includes four relatively simple one-dimensional integrals that can be calculated analytically for some types of particles (except for spheroids). A comparison with the numerical data obtained for various Chebyshev particles and finite cylinders showed that the obtained approximation yields acceptable results, even when the shape of scatterers is significantly different from spheroidal. For spheroids, our approximation coincides with the Rayleigh one.     

具有各向异性的梯度颗粒复合材料的有效介电响应

本文研究了具有梯度壳层椭球颗粒复合体系的有效介电响应，在稀释条件下用准静态近似方法, 推导了颗粒壳层具有任意介电梯度形式的椭球颗粒复合体系的有效介电常数和部分共振条件的一般表达式。并以壳层的介电常数为主轴的幂函数形式为例，得出了通过调节壳层的介电梯度形式、颗粒的结构和形状,可以提高该体系的有效介电常数和实现部分共振的结论。     

Implementation and investigation of iterative solvers in the Discrete Sources Method

The implementation of iterative methods as solvers for the Discrete Sources Method (DSM) is presented. In this method, light scattering computation linear systems with dense and relative small matrices are generated. The linear systems are traditionally solved using the QR-decomposition method. For large particles or particles with extreme geometries even this commonly stable solver can fail. In these cases, we expect that iterative methods can provide a satisfying solution nevertheless.We will present our investigation in two consecutive papers. Here, we study four different iterative solvers (RGMRES, BiCGStab, BiCGStab(l), and MinRes) considering the performance and the accuracy for typical light scattering problems. Using these iterative methods we increased the quality of a solution, especially for oblate spheroids with a higher aspect ratio. Preconditioning technique is considered in the following paper.     

Giant piezoelectric and dielectric enhancement in disordered heterogeneous systems

Effective complex piezoelectric and dielectric constants of disordered heterogeneous systems, such as statistical mixtures consisting of spheroidal particles of the same orientation but with random distribution in space, are studied. It is found for the first time that, in such systems, there exists giant piezoelectric enhancement accompanied by giant relaxation of piezoelectric coefficients and permittivity. The piezoelectric and dielectric spectra differ considerably from the Debye spectra and have a Cole-Cole character. The dependence of the effects considered on the aspect ratios of the spheroids is investigated. The physical mechanisms responsible for the anomalous behavior of the piezoelectric coefficients and permittivity are considered.     

The effect of the size,shape, and structure of metal nanoparticles on the dependence of their optical properties on the refractive index of a disperse medium

The effect of the size, shape, and structure of gold and silver nanoparticles on the dependence of their extinction and integral scattering spectra on the dielectric environment has been investigated. Calculations were performed using the Mie theory for spheres and nanoshells and the T-matrix method for chaotically oriented bispheres, spheroids, and s cylinders with hemispherical ends. The sensitivity of plasmon resonances to variations in the refractive index of the environment in the range 1.3–1.7 for particles of different equivolume size, as well as to variations in the thickness of the metal layer of nanoshells, was studied. For nanoparticles with an equivolume diameter of 15 nm, the maximal shifts of plasmon resonances due to variation in the refractive index of the environment are observed for bispheres and the shifts decrease in the series nanoshells, s cylinders or spheroids, and spheres. For particles 60 nm in diameter, the largest shifts of plasmon resonances occur for nanoshells and the shifts decrease in the series bispheres, s cylinders or spheroids, and spheres. All other conditions being the same, silver nanoparticles are more sensitive to the resonance tuning due to a change in the dielectric environment.     

On the validity of modeling rayleigh scatterers by spheroids

For single homogeneous lossy dielectric Rayleigh particles, the manner in which the particle shape affects the absorption and scattering is examined by computing the polarizability tensor elements using a general purpose program valid for any rotationally symmetric body. The results for a number of generic shapes are compared with those for the corresponding spheroids. When the dielectric constant is real and negative, the resonances attributable to bulk and to shape-dependent or surface polariton modes are tracked as functions of the length-to-width ratio of the particle and the bulk dielectric constant. The adequacy of using spheroid approximations for individual particles in a dispersion is also discussed.     

Surface plasmons in small layer metal particles

Positions of two surface plasma resonances in a layer metal particle containing a dielectric foreign nucleus are calculated within the random phase approximation. The results obtained give a possible explanation of a recent experimental observation by Genzel et al.¹ of two peak structure in the photoabsorption cross-section of small silver particles.     

Interaction dynamics of bright solitons in linearly coupled asymmetric systems

In this research work, composite media based on metamaterials including random distribution of spherical nanoparticles in a polymeric foam host are suggested to achieve negative effective refractive index in the visible spectrum. For this purpose structures including single, two and three layer spherical particles are investigated. Based on simulation results, media including single layer spheres (metallic and dielectric particles) and two layer nanospheres (core–shell particles consist of metallic core and dielectric shell) based on superposition of nanoparticles with different sizes and fill fractions are proposed for desired result. In this work, to obtain optimized band with negative RI media, superposition of three layer nanoparticles and doped semiconductor are designed.

     

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.