A.Mie米散射理论基础

1、精心整理米散射(Miescattering);又称 粗粒散射”。粒子尺度接近或大于入射光波长的粒子散射现象。德国 物理学家米(GustavMie,18681957)指出，具散射光强在各方向是不对称的，顺入射方向上的前向散射最强。粒子愈大，前向散射愈强。 米散射当球形粒子的尺度与波长可比拟时，必须考虑散射粒子体内电荷的三维分布。此散射情况下， 散射粒子应考虑为由许多聚集在一起的复杂分子构成，它们在入射电磁场的作用下，形成振荡的多极子，多极子辐射的电磁波相叠加，就构成散射波。又因为粒子尺度可与波长相比拟，所以入射波的相位在粒子上是不均匀的，造成了各子波在空间和时间上的相位差。 在子波组合产生散射波

2、的地 方，将出现相位差造成的干涉。这些干涉取决于入射光的波长、粒子的大小、折射率及散射角。当 粒子增大时，造成散射强度变化的干涉也增大。因此，散射光强与这些参数的关系，不象瑞利散射那样简单,而用复杂的级数表达，该级数的收敛相当缓慢。这个关系首先由德国科学家G.米得出，故称这类散射为米散射。它具有如下特点：散射强度比瑞利散射大得多，散射强度随波长的变化 不如瑞利散射那样剧烈。随着尺度参数增大，散射的总能量很快增加，并最后以振动的形式趋于一 定值。散射光强随角度变化出现许多极大值和极小值，当尺度参数增大时，极值的个数也增加。 当尺度参数增大时，前向散射与后向散射之比增大，使粒子前半球散射增大。当尺

3、度参数很小时， 米散射结果可以简化为瑞利散射；当尺度参数很大时，它的结果又与几何光学结果一致；而在尺度 参数比较适中的范围内，只有用米散射才能得到唯一正确的结果。所以米散射计算模式能广泛地描述任何尺度参数均匀球状粒子的散射特点。19世纪末，英国科学家瑞利首先解释了天空的蓝色：在清洁大气中，起主要散射作用的是大 气气体分子的密度涨落。分子散射的光强度和入射波长四次方成反比，因此在发生大气分子散射的 ：尸 C J I日光中，紫、蓝和青色彩光比绿、黄、橙和红色彩光为强，最后综合效果使天穹呈现蓝色。从而建 立了瑞利散射理论。20世纪初，德国科学家米从电磁理论出发，又称粗进一步解决了均匀球形粒子的散射问

4、题， 建立了米散射理论，粒散射理论。质点半径与波长接近时的散射，特点：粗粒散射与波长无关， 对各波长的散射能力相同，大气较混浊时，大气中悬浮较多的的尘粒与水滴时，天空呈灰白色。米散射理论是由麦克斯韦方程组推导出来的均质球形粒子在电磁场中对平面波散射的精确解。一般把粒子直径与入射光波长相当的微粒子所造成的散射称为米散射。 米散射适合于任何粒子尺度 只是当粒子直径相对于波长而言很小时利用瑞利散射、 很大时利用夫琅和费衍射理论就可以很方便 的近似解决问题。米散射理论最早是由 G1 Mie在研究胶体金属粒子的散射时建立的。190孙，米氏通过电磁波的麦克斯韦方程，解出了一个关于光散射的严格解，得出了任意

5、直径、任 精心整理意成分的均匀粒子的散射规律，这就是着名的米氏理论4-6。根据米散射理论，当入射光强为。粒子 周围介质中波长为人的自然光平行入射到直径为D的各向同性真球形粒子上时，在散射角为8,距离 粒子r处的散射光和散射系数分别为：从上式中可以看到，因为是各向同性的粒子，散射光强的分布和小角无关。同时，上式中：i1、i2为散射光的强度函数；s1、s2称为散射光的振幅函数？为粒子的尺寸参数（a=兀D/入）;m=m1+im2 为粒子相对周围介质的折射率，当虚部不为零时，表示粒子有吸收。对于散射光的振幅函数，有： 式中an、bn为米散射系数，其表达式为：其中：方，是半奇阶的第一类贝塞尔函数；耳照4

6、二）是第二类汉克尔函数；Pn（cos8）是第一类勒 让彳惠函数;P（1）n（cos9）是第一类缔合勒让德函数。Mie散射理论Mie散射理论是麦克斯韦方程对处在均匀介质中的均匀颗粒在平面单色波照射下的 严格数学解。由Mie散射知道，距离散射体r处p点的散射光强为式中：儿为光波波长；I0为入射光强；Isca为散射光强；6为散射角；中为偏振光的偏振角。式中：S1（e ）和S2（e ）是振幅函数；an?口 bn是与贝塞尔函数和汉克尔函数有关的函数 产n和7 n是连 带勒让得函数的函数，仅与散射角9有关。其中FF式中：邛n（a）和、（立）分别是贝塞尔函数和第一类汉克尔函数 9 n但）和名n （。）是平口

7、依）和 、（口）的导数；。为无因次直径F=nD/九，D为颗粒的实际直径 法是入射光的波长；m是散射颗粒 相对于周围介质的折射率，它是一个复数 虚部是颗粒对光的吸收的量化。由以上公式可见，Mie散射计算的关键是振幅函数S1（e）和S2（e）,它们是一个无穷求和的过程，理论上无法计算。求解振幅函 数的关键是计算anftbn,所以Mie散射的计算难点是求解an和bn。Mie散射理论的数值计算通过以上分析可知，Mie散射计算的核心是求解anffibn,我们编制程序也是围绕它进行编写。在an fr和bn的表达式中、（"），% （u） 3n （u）和、（二）满足下列递推关系：这些函数的初始值为；

8、精心整理精心整理与散射角有关的3)和,(口)满足下列递推公式：有了这些递推公式可以很方便地通过计算机程序求解。但是对于n的大小，因为计算机不可能计算无穷个数据，所以n在计算之前就要被确定。散射理论基础与Matlab实现若散射体为均匀球体，如图1所示，照射光为线偏振平面波，振幅为E,光强。沿z轴传播，其电场矢 量沿x轴振动。散射体位于坐标原点O,P为观测点。散射光方向(OP方向)与照射光方向(z轴)所组成 的平面称为散射面，照射光方向至散射光方向之间的夹角 8称为散射角，而x轴至OP在xy平面上投影 线(OP')之间的夹角小称为极化角。观测点与散射体相距r。根据经典的Mie散射理论，散射

9、粒子的 尺度参数为a =2冗a/入洪中a为球形粒子的半径，散射粒子相对周围介质的折射率为 m=m1+i*m2。则 散射光垂直于散射面和平行于散射面的两个分量的振幅函数.: IX + 支尸 / /为：以上式中：Jn+1/2(z)和Yn+1/2(z)分别为半整数阶的第一类，第二类贝塞尔函数。Pn(cos8 )为一阶n次第一类缔合勒让德函数；Pn(cos8)为第一类勒 让德函数。在数值模拟过程中选取初始下：微粒子对光的散射和吸收是电磁波与微粒子相互作用的重要特征，而微粒对电磁辐射的吸收与 散射与粒子的线度有密切关系，对于不同线度的粒子必须应用不同的散射理论。Mie散射理论主要用于从亚微米至微米的尺寸

10、段；在微米以下至纳米的光散射则近似为形式更明晰简单的瑞利散射定律 散射光强烈依赖于光波长 入(I入-4);而对大于微米至毫米的大粒子则近似为意义明确的夫朗和费 衍射规律了。Mie散射理论给出了球型粒子在远场条件下的散射场振幅an、bn以及粒子内部电磁场振幅cn、dn的计算表达式，通常称为Mie散射系数式中m表示微粒子外部介质的相对折射率，x= k a,a为球的半径i二2冗/人称为波数,为相对磁 导率，即球的磁导率与介质磁导率的比值，jn（x）和hn（x）分别为第一类虚宗量球Besse函数和Hankell 函数。散射系数，消光系数及偏振状态下散射相位函数：散射截面（Tsca（散射率Qsc,、吸收

11、截面crabs（吸收率QabS）、消光截面ext（消光率Qext）、后向散射 截面bb（后向散射率Qb）以及辐射压力（Tpr（辐射压力效率Qpr）。其表达式如下：其中i为sca abs、ext、pr分别表示散射、吸收、消光、辐射压力。按照能量守恒定律有：犬。frtfgz + bub+:Qb （后向散射系Qpr （辐射压力效率的计算公式）：这些都是无穷级数求和，在实际计算过程中必须取有限项，BohrenHuffman给出了级数项最大 值取舍的标准：对于单位振幅入射波经微粒散射后，具散射场振幅的大小与散射角有关，在球坐标系下，远场散 射振幅的大小为：其中S1和S2为散射辐射电场在垂直及平行于散射面

12、的两个偏振分量。微球内部场振幅计算公式颗粒内部电场强度为：其中M（1）o1n和N（1）e1n为矢量波球谐函数，在球坐标系中定义如下：吸收截面Qabs具有损耗介质颗粒的吸收截面为：其中e 是粒子相对介电常数的虚部，经整理可得：式中mn、nn为：实际上由Mie散射理论可知，上式中的积分项为电场强度的平方对角度8、小全空间积分的平均值即：于是吸收效率为：式中x'=rk二z/m。当xn1时即瑞利散射情况，颗粒的内部平均场强为常数，其值为:精心整理Mie计算存在的问题就是如何最有效地构造 Mie计算，同时保证准确性和避免数值的不稳定性和 病态。Mie计算以耗时着称，首先无穷项级数 N的求和，例如

13、：100Nm的水滴在0.5 Nm的可见光散 射情况下，大约需1260项求和。其次，典型的计算都希望能对一系列半径（如对尺寸分布求积分）、 一系列波长（如对太阳光谱求积分）及一系列折射率求和（如通过散射参量反推折射率）。当折射率虚部mim很大时，用向后循环法求 An很不稳定。而向前递推总是稳定的 （但向后递推安全时， 总是优先选择， 因为其计算速度很快）。得出允许向后递推的经验标准：用正确的向前地推与相对应的向后地推做比较，当发现对Q飙曰和g的相对误差超过10-6时，认为计算失败。对于一对确定的（x,mRe）,我们采用向后递推寻找第一个循环失败的mim研究表明：对于确定的 41R。，的值随着x的

14、增加很快趋向于一个确定值。Zi(mRcf 二一8 +- 0.4474m%+ 60020 W 0.0001751 一 对如果在任意角度下 蚪、S2的实部和虚部的相对误差超过10-5时，认为对§和S2的向后递推失败。（而此时，QsQe并不受影响，因为当G，S2的相对误差达到10一5时，QscaQext的相对误差 总维持在10-10以下。）3. + 3 9 对 S1 和 &/式- 16.35m+ 8.42/«%-15必对散射强度岛2 +岛/和偏正度（|S/|Si|2）/（|ST+Si|2）连分式算法总结：Mie散射计算的核心是计算an和bn其中3 n( a )= a Jn

15、(民)，七 n( a )=a Jn（ a ）+i a Yn（ a ）,Jn和Yn分别是第一和二类贝塞耳函数，a称为当量直径，a =2冗r/人,是球形颗粒的真实半径，入是入射光的波长，1为折射率式中p为函数任一自变量。贝塞耳函数递推关系式：Mie散射计算中Jn、Yn、Dn的计算是关键和难点。对于Dn,我们采用的是Lentz的连分式的算法：Lentz证明有如下关系：其中，社=r 1， L ko. 5人我们注意到当上T8时，精心整理精心整理/ "去。足80所以可以利用上式累积相乘直到满足精度要求。（可根据精度要求 例如10-7来确定所要达到的k值）对于Jn、Yn的生成本文也采用连分式的算法

16、。具体方案如下 ：令Cn=Jn-1（ OC ）/Jn（ 0C，根据贝塞耳差积公式：/I M （ 夕）/ C.，7 Y置 . 1，夕2上 /=T /"i由以上二式整理得：- -上式中Cn的计算是采用类似于Dn的连分式的形式，计算中可调用同一函数计算。Jo = siii/ 优；Ji f 戊）=sin/ "2 cos/ s； :若已知初值：，一.一:这样就可计算出各级JnffiYn。WilliamJ.Lentz关于连分式的文章： （ 1 1 X,J1,"%-%4*诙1力苴中n + % = V, w = m工只十'.一。以出工（幻为基础，采用贝塞尔函数比值的连分式

17、表示法：,利用此法可产生所有的 出工匕），尽管耗时，但能减少存储需求。同时可通过计算高阶J、N）/Jv-1（N）值，使用下面的递推公式，从后往前算出其他值。不像一般的函数，贝塞尔函数的比值一旦超过可控制的边界，就不再增长，初始的高阶A（n）值决定了所有低阶值的准确性，因此，采用新方法计算准确的初始比值是必要的。处于分母位置的+号表示分母上加上一个特殊的连分式。类似于上式 /中的表示形式。定义一种新的符号：Lentz给出了 n阶部分收敛值为：例如：实变量工=L°,虚数"=9.5:计算过程：米散射学习目前所遇到的困难：到底怎样的计算结果才算正确，如何能找到一个米散射计算结果准确

18、又有效的数据库，来验证自己算法及程序的正确性。精心整理精心整理倒退式算法的总结：当*小于0刃2时,可以利用RayE处公式计算，而当D比2m00更大时，几何光学即可适用 Dn的计算采用Dave的倒推式:由于Dn函数有很强的收敛性，对于Dn的倒推计算的初值的选取有很强的随意性。因为当n-oo时Dn(ma)-0,所以可以取0乍为初值。倒推起点选取大一些，可以保证Dn函数的收敛完全，但是同时却 增加了计算时间。所以必须选取一个最佳的选择标准。通过试算，作者认为最佳的上限为这里ml是复折射率的实部.同样，对于贝塞耳函数Jn的计算也可以用倒才t的方法计算产生：上式是一个普通的Jn的递推式，知道了 JnffiJn-1可以顺利地计算出所有的Jn序列值。为了避免计算Jn