AMie米散射理论基础

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

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

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

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

5、据米散射理论，当入射光强为I0，粒子周围介质中波长为人的自然光平行入射到直径为D的各向同性真球形粒子上时，在散射角为9,距离粒子r处的散射光和散射系数分别为从上式中可以看到，因为是各向同性的粒子，散射光强的分布和小角无关。同时，上式中：H=s(m,0)a)Xj*(rtjr61a)J2=si(mt0a)X57(mt01i1、i2为散射光的强度函数；si、s2称为散射光的振幅函数；a为粒子的尺寸参数(a=兀D/入);m=mi+im2为粒子相对周围介质的折射率，当虚部不为零时，表示粒子有吸收。对于散射光的振幅函数，有：M2k+1$】1«/«-I-1)(+九七，n=1n(n+1/O

6、O-*7内1/工ts2/1匕f"”工打+«=in(n+17式中an、bn为米散射系数，其表达式为：jI#*1I/JIJJ中i2.f0/加aJ-mg式fa，.f加ai，打，n)，n(mn)-n(a)3rlfma)_，心uf3if酬)界mZn(a)St/rna)其中:riaQ=看尸吟d尸fcosQ17=浦dfco助S=肃产产feMJ几十之一,是半奇阶的第一类贝塞尔函数；'二,是第二类汉克尔函数Pn(cosO)是第一类勒让彳惠函数；P(1)n(cosO)是第一类缔合勒让德函数。Mie散射理论Mie散射理论是麦克斯韦方程对处在均匀介质中的均匀颗粒在平面单色波照射下的严格数学

7、解。由Mie散射知道，距离散射体r处p点的散射光强为/空=log产g(a中/(Qsi(9)|2sin3(P-|si(9)|2cos2(P式中：九为光波波长；I0为入射光强；Isca为散射光强；8为散射角;中为偏振光的偏振角8、12打+51(3=白汽5+D麻+412曾+irfn(+1)he兀跖1<30S2(0=式中：§。)和0(6)是振幅函数；an和bn是与贝塞尔函数和汉克尔函数有关的函数；nn和7n是连带勒让得函数的函数，仅与散射角日有关。其中_W(o()31加/)-?收51(o()9?64&)OL6(仪)中1(所0()加二(0()可。*0()_.一(0)5丁a)-T；

8、(x)加加a)”一mG(仪)；所仪).£n(0()(Z?(ma)式中：中n(u)和、(。)分别是贝塞尔函数和第一类汉克尔函数；甲；3)和名n(。)是外3)和际(”)的导数；口为无因次直径，口=nD/，-D为颗粒的实际直径；九是入射光的波长；m是散射颗粒相对于周围介质的折射率，它是一个复数，虚部是颗粒对光的吸收的量化。由以上公式可见,Mie散射计算的关键是振幅函数6(e)和S20),它们是一个无穷求和的过程，理论上无法计算。求解振幅函数的关键是计算anffibn,所以Mie散射的计算难点是求解an和bn。Mie散射理论的数值计算通过以上分析可知，Mie散射计算的核心是求解an和bn,我

9、们编制程序也是围绕它进行编写。在an和bn的表达式中n(«)Wn(a),sn(a)和8n(口)满足下列递推关系："1(X)-皿式MWi<P；(od=-"+Wi(odG(ft)=：G4】(X)-G+2(oO£/i(o0=-&G(ft)+G.i这些函数的初始值为；q?I(0t)=CO50(W(K)=smttci(ot)=cosOf-ismttG(ot)=51noe-icostt与散射角有关的平n(a)和8np)满足下列递推公式：r-=茂oosB山sin'8-JTt-icost>-TO-2n-1w-1TT斤二二*-1）后*1+TTn

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

11、入，其中a为球形粒子的半径，散射粒子相对周围介质的折射率为m=m1+i*m2。则散射光垂直于散射面和平行于散射面的两个分量的振幅函数为s：叫(幻yj加姆mW'Ja、g，"%=(a)wfn()na)-(a)gf加为Pnfa)卜;f理0)-里；(仪)里淑仪)&mn(a)Wn(ma)-5/幻甲/m©尸J'(8.伪Tn=sinek"n六?，&初以上式中：/7=F上FJn+1/2(z)和Yn+1/2(z)分别为半整数阶的第一类，第二类贝塞尔函数Pn(cos8)为一阶n次第一类缔合勒让彳惠函数；Pn(cos9)为第一类勒让德函数。在数值模拟过程

12、中选取初始下：fsinz+/cosz)-fcos-;sinz)弧=-sinz-cosz丐=8S8J=1微粒子对光的散射和吸收是电磁波与微粒子相互作用的重要特征，而微粒对电磁辐射的吸收与散射与粒子的线度有密切关系，对于不同线度的粒子必须应用不同的散射理论。Mie散射理论主要用于从亚微米至微米的尺寸段；在微米以下至纳米的光散射则近似为形式更明晰简单的瑞利散射定律，散射光强烈依赖于光波长入(I入-4);而对大于微米至毫米的大粒子则近似为意义明确的夫朗和费衍射规律了。Mie散射理论给出了球型粒子在远场条件下的散射场振幅an、bn以及粒子内部电磁场振幅cn、dn的计算表达式，通常称为Mie散射系数mFm

13、x)3K皿1工)f博xjjf)"injn(inx)xhJ(x)-f艾)mxj/mx>Jtt.短fxi/幻1-3短mx”d7再5(xhn/x)-方(x)mxjn(wxjxh?G)-"i工号#CvJfxjn(nix)xh(x)'小切力f工)1式中m表示微粒子外部介质的相对折射率,x=ka,a为球的半径，k=2冗/入称为波数，仙为相对磁导率，即球的磁导率与介质磁导率的比值,jn(x)和hn(x)分别为第一类虚宗量球Bessel函数和Hankell函数。散射系数，消光系数及偏振状态下散射相位函数：ksca=A冥2*+1)(aookg=彳X72n+1/Ruf。有+bjP

14、(6)=nV(an2+1bnV散射截面（Tsca（散射率Qsc、吸收截面（Tabs（吸收率Qab、消光截面（Text（消光率Qext）、后向散射截面（7b（后向散射率Qb）以及辐射压力（Tpr（辐射压力效率Qpr）o其表达式如下：其中i为sca、abs、ext、pr分别表示散射、吸收、消光、辐射压力。按照能量守恒定律有：0H=Qsca+Qahcr或笃k=+%为Qpr（辐射压力效率的计算公式）：Qpr-Qsxt-Qsca<COS8.>。如<8与白>=改式官田;+瓦b；2+元一；.工I“k打+1hM力f理+1，）Qb（后向散射系数）：Qb=2|彳2"+1J（-Vn

15、（a-bn）这些都是无穷级数求和，在实际计算过程中必须取有限项,Bohren和Huffman给出了级数项最大值取舍的标准："2=-V+4A+2对于单位振幅入射波经微粒散射后，具散射场振幅的大小与散射角有关，在球坐标系下，远场散射振幅的大小为：ikrEtf=ZCO*S2(8S6)-ikrikr£二-sinS'$i(ws0)ikr其中S1和S2为散射辐射电场在垂直及平行于散射面的两个偏振分量。微球内部场振幅计算公式颗粒内部电场强度为：8Ei=丁;j田叙d,N2/nfn+1，其中MdMn和Ne1n为矢量波球谐函数，在球坐标系中定义如下：M1=cos4Feos4人、-sin

16、巾hrn(cos0)jn(n+I)co§6*sin/rmx)(rmx)L力(QQsd)rmxNd«=COS小-sin4L吸收截面Qabs具有损耗介质颗粒的吸收截面为：仃ah:二k其中&是粒子相对介电常数的虚部80%左=k/X%式中mn、nn为：,门1rz*xi/rmx)7vnfCOSP>rjnxl/&frm、：mfr沅013Ffco5处rmx:/JkL，经整理可得：f*Cn2+dn2)r2dr0徵+I2实际上由Mie散射理论可知，上式中的积分项为电场强度的平方对角度8、小全空间积分的平均值，即：于是吸收效率为。心=eIJ）工”式中x'=rk=z

17、/m。当xn1时即瑞利散射情况，颗粒的内部平土§场强为常数，其值为:ImprovedMiescatteringalgorithmsW.J.WiscombeMie计算存在的问题就是如何最有效地构造Mie计算，同时保证准确性和避免数值的不稳定性和病态。Mie计算以耗时著称，首先无穷项级数N的求和，例如:100Rm的水滴在0.5Nm的可见光散射情况下，大约需1260项求和。其次，典型的计算都希望能对一系列半径（如对尺寸分布求积分）、一系列波长（如对太阳光谱求积分）及一系列折射率求和（如通过散射参量反推折射率）。Qut=FE(2m+I+bn).工上n"12N=WE(2以任+|垢,X

18、Sn=IWn(n+2)tad*Rd。”以界+Jl-lRe(an6；)»N2n+1S3石许/v2n+1rS2(M)=£-lanTn(u)+fenirn()fn(n+I)当折射率虚部mim很大时，用向后循环法求An很不稳定。而向前递推总是稳定的（但向后递推安全时，总是优先选择，因为其计算速度很快）。得出允许向后递推的经验标准：用正确的向前地推与相对应的向后地推做比较，当发现对Qsch和g的相对误差超过10-6时，认为计算失败。对于一对确定的（x,mRe）,我们采用向后递推寻找第一个循环失败的mfm研究表明：对于确定的41Rs,工的值随着x的增加很快趋向于一个确定值。xmfm&g

19、t;minCcmim)allx.8+26.22191)0.4474m也+0.00204喊日-0.000175mM.£4002K02000ieoo*400S4CT）4000BOO2coI£34S6TS如果在任意角度下6、s2的实部和虚部的相对误差超过105时，认为对S1和S2的向后递推失败。（而此时，QscaQxt并不受影响，因为当S1,$2的相对误差达到10V时，QscaQext的相对误差总维持在10-10以下。）scadxi又tS1和S2=13,78nI410,8mr*+39=16,35mL+15.04,对散射强度|Sip+|S2p和偏正度(展-|Si|2)/(|S/+S

20、i|Z)连分式算法总结：Mie散射计算的核心是计算an和bn必，为fma)-川"h，幻小用00刃用刃加二九(ma)打心川5(加勿九，幻山州二仇=,*mZn(a)九(mg-fM(a)ipn(根戊，其中6n(a)=aJn(a),En(a)=aJn(a)+iaYn(a),Jn和Yn分别是第一和二类贝塞耳函数，a称为当量直径，a=2兀r/入，r是球形颗粒的真实半径，入是入射光的波长，m为折射率中网=(H2)sj*L)，=/2)叫加+J1)以11AFDn(p)=ln0fp>=i.dpWJP)式中p为函数任一自变量。贝塞耳函数递推关系式：f>t/l也,九，幻=if幻-、*i岛fa)E

21、丹/勾=?打4a为Dn(m0)/mE/R;Dm/加色，/?力一咒/a九，a)-“h-ifml)望fn偿)+a小/a).1niDn(m()+n/町Wja)Mie散射计算中Jn、Yn、Dn的计算是关键和难点。对于Dn，我们采用的是Lentz的连分式的算法D“旅矽=.能Lentz证明有如下关系fma);f©G7faka、J鞭(ma)a1工程电/其中,1厂一»"+卜05.人我们注意到当上78时,"不。所以可以利用上式累积相乘直到满足精度要求。(可根据精度要求例如10-7来确定所要达到的k值)对于Jn、Yn的生成本文也采用连分式的算法。具体方案如下令Cn=Jn-1

22、(a)/Jn(a)，根据贝塞耳差积公式:JJRY人"勾J八"刃Y"幻=才？由以上二式整理得:J八1/C-fa)Y相ifa)餐)Y”®=JW上式中Cn的计算是采用类似于Dn的连分式的形式，计算中可调用同一函数计算。若已知初值：Jo=sdna/a；a)=sin/后cos/cm；,v2.Yq=-cos/a;Yi(=-cos/)-sin度/四；这样就可计算出各级Jn和YnW川iamJ.Lentz关于连分式的文章:九=（打2/2产心龙幺Q）=（值/2）啊心.姑+（-1）丐上心使W/）-2/其中n+%=%N=m工。以ArG）为基础，采用贝塞尔函数比值的连分式表示法：

23、式力,利用此法可产生所有的4，工匕），尽管耗时，但能减少存储需求。同时可通过计算高阶值，使用下面的递推公式，从后往前算出其他值。兀卜）_2"_j（z）1g=；Q）芸jrr-lQ）1+-2（+1U-1箝=2K不像一般的函数，贝塞尔函数的比值一旦超过可控制的边界，就不再增长，初始的高阶出工（胃）值决定了所有低阶值的准确性，因此，采用新方法计算准确的初始比值是必要的。+2（+2）源1-2B+3）父1+.处于分母位置的+号表示分母上加上一个特殊的连分式。类似于上式/中的表示形式。定义一种新的符号：n1,111/=,Gi+I+也Lentz给出了n阶部分收敛值为：例如：实变量=L0,虚数炉=9.

24、5:计算过程:）1-1,2,3,.事%+工=-20.94736842%工（T产以廿十言一1）/aaJ=19%=-21,"J=几1_(19)(-20.94736842)(223”226131)卜24.95g43131)(26.95993017?Jv=(-21)22?95238095)(-24.95643154)(26.95993017)=18.95228198,米散射学习目前所遇到的困难：到底怎样的计算结果才算正确，如何能找到一个米散射计算结果准确又有效的数据库，来验证自己算法及程序的正确性。倒退式算法的总结：由于Mie级数的收敛速度随特征值X的增大而减慢，因此不同的万值即使在相同的计

25、算精度要求下所需级数项数M也不一样.1979年美国学者WJWiscomt/s】在大量计算的基础上，参考前人的工作总结出一个X；的经验计算公式.利用该公式给出的项数可以使整个计算误差小于10.Wiscombe公式是工+4立3+Ix+4.05x3+10.02w工S88<x<42001工+4%5+24200三#W20000当X小于0.02时可以利用lluykigh公式计算,而当X比20000更大时.几何光学即可适用.Dn的计算采用Dave的倒推式：八/b-ikn/tncl4Dn（/椁）£>,1（如=ma由于Dn函数有很强的收敛性，对于Dn的倒推计算的初值的选取有很强的随意性。因为当n-o"Dn(ma)一0，所