AMie米散射理论基础

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

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

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

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

5、散射的严格解， 得出了任意直径、任意成分的均匀粒子的散射规律，这就是着名的米氏理论4-6。 根据米散射理论，当入射光强为I0,粒子周围介质中波长为人的自然光平行入射到 直径为D的各向同性真球形粒子上时，在散射角为8 ,距离粒子r处的散射光和散射 系数分别为：从上式中可以看到，因为是各向同性的粒子，散射光强的分布和小角无关。同时， 上式中：i1、i2为散射光的强度函数；s1、s琳为散射光白振幅函数？为粒子的尺寸参数(a= 九D/入);m=m1+im2为粒子相对周围介质的折射率，当虚部不为零时，表示粒子有 吸收。对于散射光的振幅函数，有：式中an bn为米散射系数，其表达式为：其中：J q-L (

6、 - i_L / - j是半奇阶的第一类贝塞尔函数 户二f”是第二类汉克尔函 数;Pn(cos8 )是第一类勒让彳惠函数；P(1)n(cos8 )是第一类缔合勒让德函数。Mie散射理论Mie散射理论是麦克斯韦方程对处在均匀介质中的均匀颗粒在平面单色波照射下的严格数学解。由Mie散射知道，距离散射体r处p点的散射光强为 式中：为光波波长；I0为入射光强；Isca为散射光强；为散射角;为偏振光的偏 振角 式中：S|()和$()是振幅函数；an和bn是与贝塞尔函数和汉克尔函数有关的函 数；n和n是连带勒让得函数的函数，仅与散射角 有关。其中式中：n()和n()分别是贝塞尔函数和第一类汉克尔函数；n(

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

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

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

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

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

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

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

14、究表明：对于确定的iR% ,想认工的值随着x的增加很快趋向于一个确定值。八加5二一8 + 26.22用30i4474m.n+ 0W04m# - GOOOI75m%如果在任意角度下Sp S2的实部和虚部的相对误差超过10 5时，认为对和S2的向后递推失败。(而此时，QscaQext并不受影响，因为当G, S2的相对 误差达到10 5时，QscaQext的相对误差总维持在10 10以下。)/恤如)=13.78m IOEjhr* + 3.9, 壬利b和S2/式呐)=1&35温153,对散射强度岛|，岛|2和偏正度(|S2p-|S1|2)/(|S2P+S1p)连分式算法总结：Mie散射计算的核心是计算

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

16、采用连分式的算法。具体方案如下 令。=41（ a ）/Jn（ a ）,根据贝塞耳差积公式:由以上二式整理得：_/J/二f夕)/ Cm (以)Y心fa)4J M - 1(上式中Cn的计算是采用类似于Dn的连分式的形式，计算中可调用同一函数计 算。若已知初值：Jo = sin/ 口； Ji = sin/ H co&t/ 谢；Yo = - cos/ a; Yi （a） = cosW/ 4ni/ aw；这样就可计算出各级JnffiYn。W川iamJ.Lentz关于连分式的文章：j= W + =十比（号）。和）2血苴中鹿+ % = % z = mx以4，】（Z）为基础，采用贝塞尔函数比值的连分式表示法：

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

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