MIE理论小结汇编.docxVIP

§1 Rayleigh 散射与Mie 散射 天空呈蓝色和线偏振这两个特性在上个世纪曾经是很重要的科学之谜，最后由Rayleigh做出解释，Rayleigh注意到了产生这种散射的散射元并不是当时一般人们所认为的水或冰，而是气体分子本身造成这样的散射。 Rayleigh散射具有（1）散射光的强度和波长的四次方成反比，因此晴朗天空基本上是蓝色而不是太阳光的颜色；（2）散射光的空间分布与观测方向有一个简单的关系；（3）在散射的前半球和后半球具有相同的散射强度；（4）方向的散射光几乎是全偏振的。 空气分子的Rayleigh 散射代表大气光路散射的最小值。在低层大气，粒子的Mie散射处于主导地位，但平均来说，随着高度的增加，Mie散射的减小比Rayleigh散射快。这是因为霾气溶胶的标高一般近似1Km，而恒常气体层的标高一般为6—9Km甚至更高。因此除了几个稀薄的粒子层之外，气溶胶只在一个有限高度范围内分布，而Rayleigh散射却一直到很高的高度上仍有影响。但到这样的高度，大气透明度已经很高了，故限制视程的主要原因是几何因子而不是光学因子。[1] 相对于云雾粒子，我们主要考虑这种由直径大于波长的0.03倍的粒子造成的散射，即Mie散射。从很小的粒子开始，当其半径相对于波长而言逐渐加大时，就逐渐发生从Rayleigh散射向Mie散射的过度。Mie散射具有（1）散射光强度随角度分布变得十分复杂，粒子相对于波长的尺度越大，分布越复杂。（2）当粒子的尺度加大时，前向散射与后向散射之比随之增加，结果使前向散射的波瓣增大。（3）当粒子尺度比波长大时，散射过程和波长的依赖关系就不密切了，这一点可以从云一般是发白的现象推测到。白色的云和蓝色天空反映了两种不同类型的散射。当r 0.03时， Rayleigh近似式和Mie散射公式相比，误差在1％以内。[1]§2单球的Mie散射Mie理论自1908年被提出，它给出了均匀介质球引起平面电磁波散射的精确解。考虑半径为a,相对背景的折射率m,平面电磁波沿z方向入射，电矢量沿轴方向极化。考虑到在球坐标系中，入射平面波可写为：。球内外场满足波动方程：式中m 为球粒子的复折射率。将入射电场与磁场，球形粒子内部场与，以及散射场与分别用矢量球谐函数展开，其中包含四个待定系数。再利用边界条件： ， 将入射场，内场和散射场带入边界条件，并利用勒让德函数和三角函数的正交性，获得四个待定系数满足的线性方程组。通过求解这个线性方程组可以得到：式中， 为球贝塞尔函数，为第一类汉克尔函数。则在远离球粒子的散射场：其中和 称为振幅函数（amplitude function), 根据Mie理论, 可由下面的对称级数给出：其中称为尺寸参数（size parameter）,定义为,它代表粒子周长和入射波长的比值，此外，m 为折射指数（index of refraction）为散射角。和称为角度系数（angular coefficient），由下式给出： 其中，是连带Legendre函数。 可以发现当m趋于1时，，即将变为零，这就是当粒子消失时，散射场自然也会消失。同样得到散射系数与消光系数以及入射散射的斯托克斯参量（stokes parameters）：[2]散射系数：消光系数：单次反照率 斯托克斯参量矩阵：其中，这些矩阵元中，只有三个是完全独立的，即：.§3散射系数和散射截面的计算[2]据此Aden在1951年首先在计算单粒子散射问题时，引进“对数导数”，即。我们利用这一结果，改写成下列形式：而对数导数满足递推关系：。对于此式，Kattawar 和Plass在1967年发表的论文中指出在向下生成时是“数值稳定的”，而另两个函数，的计算采用向上生成。至于截断次数Nstop的取值，根据Wiscombe在1979年和1980年的研究和实验成果，我们可以取为最接近的整数。的初值和NMX的选取，需要，这是向下生成稳定性的要求。我们将取值，并取，这样一般说来计算结果是合理的。和的生成公式是相同的，即取双精度初值为，以及单精度初值最后根据复数得到。考虑计算的效率问题，我们只需要计算出，然后利用关系式和，求的的。最后，根据前面的公式计算出散射矩阵元并计算了散射问题的吸收系数，散射系数和消光系数，以及极化程度。其中，计算出的散射矩阵元对所有的散射角应该满足关系。§4单球散射结果分析。（1）令折射指数m=1.5,考虑消光系数的计算结果随尺寸参数的变换。 由图可以看出在x=4.5处，即波长为粒子半径的4/3时，可以看出强烈的散射发生，此外，这个效率因子随着x的增加一面振动一面减小，最后接近于2，也就是说，对于大粒子来说，入射于其截面的能量的两倍被削弱了，在这种情况下只有散射，因为m为实数，这种现像是几何光学作用，即反射，折射。（2）参考书中的初始数据集