AMie米散射理論基礎(chǔ)

1、米散射(Miescattering);又稱粗粒散射”。粒子尺度接近或大于入射光波長的粒子散射現(xiàn)象。德國物理學(xué)家米(GustavMie,18681957)指出，具散射光強(qiáng)在各方向是不對稱的，順入射方向上的前向散射最強(qiáng)。粒子愈大，前向散射愈強(qiáng)。米散射當(dāng)球形粒子的尺度與波長可比擬時，必須考慮散射粒子體內(nèi)電荷的三維分布。此散射情況下，散射粒子應(yīng)考慮為由許多聚集在一起的復(fù)雜分子構(gòu)成，它們在入射電磁場的作用下，形成振蕩的多極子，多極子輻射的電磁波相疊加，就構(gòu)成散射波。又因?yàn)榱Ｗ映叨瓤膳c波長相比擬，所以入射波的相位在粒子上是不均勻的，造成了各子波在空間和時間上的相位差。在子波組合產(chǎn)生散射波的地方，將出現(xiàn)相位

2、差造成的干涉。這些干涉取決于入射光的波長、粒子的大小、折射率及散射角。當(dāng)粒子增大時，造成散射強(qiáng)度變化的干涉也增大。因此，散射光強(qiáng)與這些參數(shù)的關(guān)系，不象瑞利散射那樣簡單，而用復(fù)雜的級數(shù)表達(dá)，該級數(shù)的收斂相當(dāng)緩慢。這個關(guān)系首先由德國科學(xué)家G.米得出，故稱這類散射為米散射。它具有如下特點(diǎn)：散射強(qiáng)度比瑞利散射大得多，散射強(qiáng)度隨波長的變化不如瑞利散射那樣劇烈。隨著尺度參數(shù)增大，散射的總能量很快增加，并最后以振動的形式趨于一定值。散射光強(qiáng)隨角度變化出現(xiàn)許多極大值和極小值，當(dāng)尺度參數(shù)增大時，極值的個數(shù)也增加。當(dāng)尺度參數(shù)增大時，前向散射與后向散射之比增大，使粒子前半球散射增大。當(dāng)尺度參數(shù)很小時，米散射結(jié)果可以

3、簡化為瑞利散射；當(dāng)尺度參數(shù)很大時，它的結(jié)果又與幾何光學(xué)結(jié)果一致；而在尺度參數(shù)比較適中的范圍內(nèi)，只有用米散射才能得到唯一正確的結(jié)果。所以米散射計(jì)算模式能廣泛地描述任何尺度參數(shù)均勻球狀粒子的散射特點(diǎn)。19世紀(jì)末，英國科學(xué)家瑞利首先解釋了天空的藍(lán)色：在清潔大氣中，起主要散射作用的是大氣氣體分子的密度漲落。分子散射的光強(qiáng)度和入射波長四次方成反比，因此在發(fā)生大氣分子散射的日光中，紫、藍(lán)和青色彩光比綠、黃、橙和紅色彩光為強(qiáng)，最后綜合效果使天穹呈現(xiàn)藍(lán)色。從而建立了瑞利散射理論。20世紀(jì)初，德國科學(xué)家米從電磁理論出發(fā)，又稱粗進(jìn)一步解決了均勻球形粒子的散射問題，建立了米散射理論，粒散射理論。質(zhì)點(diǎn)半徑與波長九接近

4、時的散射，特點(diǎn)：粗粒散射與波長無關(guān)，對各波長的散射能力相同，大氣較混濁時，大氣中懸浮較多的的塵粒與水滴時，天空呈灰白色。米散射理論是由麥克斯韋方程組推導(dǎo)出來的均質(zhì)球形粒子在電磁場中對平面波散射的精確解。一般把粒子直徑與入射光波長相當(dāng)?shù)奈⒘Ｗ铀斐傻纳⑸浞Q為米散射。米散射適合于任何粒子尺度，只是當(dāng)粒子直徑相對于波長而言很小時利用瑞利散射、很大時利用夫瑯和費(fèi)衍射理論就可以很方便的近似解決問題。米散射理論最早是由G1Mie在研究膠體金屬粒子的散射時建立的。1908年，米氏通過電磁波白麥克斯韋方程,解出了一個關(guān)于光散射的嚴(yán)格解得出了任意直徑、任意成分的均勻粒子的散射規(guī)律，這就是著名的米氏理論4-6。根

5、據(jù)米散射理論，當(dāng)入射光強(qiáng)為I0，粒子周圍介質(zhì)中波長為人的自然光平行入射到直徑為D的各向同性真球形粒子上時，在散射角為9,距離粒子r處的散射光和散射系數(shù)分別為從上式中可以看到，因?yàn)槭歉飨蛲缘牧Ｗ?，散射光?qiáng)的分布和小角無關(guān)。同時，上式中：H=s(m,0)a)Xj*(rtjr61a)J2=si(mt0a)X57(mt01i1、i2為散射光的強(qiáng)度函數(shù)；si、s2稱為散射光的振幅函數(shù)；a為粒子的尺寸參數(shù)(a=兀D/入);m=mi+im2為粒子相對周圍介質(zhì)的折射率，當(dāng)虛部不為零時，表示粒子有吸收。對于散射光的振幅函數(shù)，有：M2k+1$】1«/«-I-1)(+九七，n=1n(n+1/O

6、O-*7內(nèi)1/工ts2/1匕f"”工打+«=in(n+17式中an、bn為米散射系數(shù)，其表達(dá)式為：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=肅產(chǎn)產(chǎn)feMJ幾十之一,是半奇階的第一類貝塞爾函數(shù)；'二,是第二類漢克爾函數(shù)Pn(cosO)是第一類勒讓彳惠函數(shù)；P(1)n(cosO)是第一類締合勒讓德函數(shù)。Mie散射理論Mie散射理論是麥克斯韋方程對處在均勻介質(zhì)中的均勻顆粒在平面單色波照射下的嚴(yán)格數(shù)學(xué)

7、解。由Mie散射知道，距離散射體r處p點(diǎn)的散射光強(qiáng)為/空=log產(chǎn)g(a中/(Qsi(9)|2sin3(P-|si(9)|2cos2(P式中：九為光波波長；I0為入射光強(qiáng)；Isca為散射光強(qiáng)；8為散射角;中為偏振光的偏振角8、12打+51(3=白汽5+D麻+412曾+irfn(+1)he兀跖1<30S2(0=式中：§。)和0(6)是振幅函數(shù)；an和bn是與貝塞爾函數(shù)和漢克爾函數(shù)有關(guān)的函數(shù)；nn和7n是連帶勒讓得函數(shù)的函數(shù)，僅與散射角日有關(guān)。其中_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)和、(。)分別是貝塞爾函數(shù)和第一類漢克爾函數(shù)；甲；3)和名n(。)是外3)和際(”)的導(dǎo)數(shù)；口為無因次直徑，口=nD/，-D為顆粒的實(shí)際直徑；九是入射光的波長；m是散射顆粒相對于周圍介質(zhì)的折射率，它是一個復(fù)數(shù)，虛部是顆粒對光的吸收的量化。由以上公式可見,Mie散射計(jì)算的關(guān)鍵是振幅函數(shù)6(e)和S20),它們是一個無窮求和的過程，理論上無法計(jì)算。求解振幅函數(shù)的關(guān)鍵是計(jì)算anffibn,所以Mie散射的計(jì)算難點(diǎn)是求解an和bn。Mie散射理論的數(shù)值計(jì)算通過以上分析可知，Mie散射計(jì)算的核心是求解an和bn,我

9、們編制程序也是圍繞它進(jìn)行編寫。在an和bn的表達(dá)式中n(«)Wn(a),sn(a)和8n(口)滿足下列遞推關(guān)系："1(X)-皿式MWi<P；(od=-"+Wi(odG(ft)=：G4】(X)-G+2(oO£/i(o0=-&G(ft)+G.i這些函數(shù)的初始值為；q?I(0t)=CO50(W(K)=smttci(ot)=cosOf-ismttG(ot)=51noe-icostt與散射角有關(guān)的平n(a)和8np)滿足下列遞推公式：r-=茂oosB山sin'8-JTt-icost>-TO-2n-1w-1TT斤二二*-1）后*1+TTn

10、-2H=017o=0TTo=7T1=0有了這些遞推公式可以很方便地通過計(jì)算機(jī)程序求解。但是對于n的大小,因?yàn)橛?jì)算機(jī)不可能計(jì)算無窮個數(shù)據(jù)，所以n在計(jì)算之前就要被確定。散射理論基礎(chǔ)與Matlab實(shí)現(xiàn)若散射體為均勻球體，如圖1所示，照射光為線偏振平面波，振幅為E,光強(qiáng)I0,沿z軸傳播，其電場矢量沿x軸振動。散射體位于坐標(biāo)原點(diǎn)O,P為觀測點(diǎn)。散射光方向（OP方向）與照射光方向（z軸）所組成的平面稱為散射面，照射光方向至散射光方向之間的夾角8稱為散射角，而x軸至OP在xy平面上投影線（OP'）之間的夾角小稱為極化角。觀測點(diǎn)與散射體相距r。根據(jù)經(jīng)典的Mie散射理論,散射粒子的尺度參數(shù)為a=2兀a/

11、入，其中a為球形粒子的半徑，散射粒子相對周圍介質(zhì)的折射率為m=m1+i*m2。則散射光垂直于散射面和平行于散射面的兩個分量的振幅函數(shù)為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)分別為半整數(shù)階的第一類，第二類貝塞爾函數(shù)Pn(cos8)為一階n次第一類締合勒讓彳惠函數(shù)；Pn(cos9)為第一類勒讓德函數(shù)。在數(shù)值模擬過程

12、中選取初始下：fsinz+/cosz)-fcos-;sinz)弧=-sinz-cosz丐=8S8J=1微粒子對光的散射和吸收是電磁波與微粒子相互作用的重要特征，而微粒對電磁輻射的吸收與散射與粒子的線度有密切關(guān)系，對于不同線度的粒子必須應(yīng)用不同的散射理論。Mie散射理論主要用于從亞微米至微米的尺寸段；在微米以下至納米的光散射則近似為形式更明晰簡單的瑞利散射定律，散射光強(qiáng)烈依賴于光波長入(I入-4);而對大于微米至毫米的大粒子則近似為意義明確的夫朗和費(fèi)衍射規(guī)律了。Mie散射理論給出了球型粒子在遠(yuǎn)場條件下的散射場振幅an、bn以及粒子內(nèi)部電磁場振幅cn、dn的計(jì)算表達(dá)式，通常稱為Mie散射系數(shù)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表示微粒子外部介質(zhì)的相對折射率,x=ka,a為球的半徑，k=2冗/入稱為波數(shù)，仙為相對磁導(dǎo)率，即球的磁導(dǎo)率與介質(zhì)磁導(dǎo)率的比值,jn(x)和hn(x)分別為第一類虛宗量球Bessel函數(shù)和Hankell函數(shù)。散射系數(shù)，消光系數(shù)及偏振狀態(tài)下散射相位函數(shù)：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其表達(dá)式如下：其中i為sca、abs、ext、pr分別表示散射、吸收、消光、輻射壓力。按照能量守恒定律有：0H=Qsca+Qahcr或篤k=+%為Qpr（輻射壓力效率的計(jì)算公式）：Qpr-Qsxt-Qsca<COS8.>。如<8與白>=改式官田;+瓦b；2+元一；.工I“k打+1hM力f理+1，）Qb（后向散射系數(shù)）：Qb=2|彳2"+1J（-Vn

15、（a-bn）這些都是無窮級數(shù)求和，在實(shí)際計(jì)算過程中必須取有限項(xiàng),Bohren和Huffman給出了級數(shù)項(xiàng)最大值取舍的標(biāo)準(zhǔn)："2=-V+4A+2對于單位振幅入射波經(jīng)微粒散射后，具散射場振幅的大小與散射角有關(guān)，在球坐標(biāo)系下，遠(yuǎn)場散射振幅的大小為：ikrEtf=ZCO*S2(8S6)-ikrikr£二-sinS'$i(ws0)ikr其中S1和S2為散射輻射電場在垂直及平行于散射面的兩個偏振分量。微球內(nèi)部場振幅計(jì)算公式顆粒內(nèi)部電場強(qiáng)度為：8Ei=丁;j田敘d,N2/nfn+1，其中MdMn和Ne1n為矢量波球諧函數(shù)，在球坐標(biāo)系中定義如下：M1=cos4Feos4人、-sin

16、巾hrn(cos0)jn(n+I)co§6*sin/rmx)(rmx)L力(QQsd)rmxNd«=COS小-sin4L吸收截面Qabs具有損耗介質(zhì)顆粒的吸收截面為：仃ah:二k其中&是粒子相對介電常數(shù)的虛部80%左=k/X%式中mn、nn為：,門1rz*xi/rmx)7vnfCOSP>rjnxl/&frm、：mfr沅013Ffco5處rmx:/JkL，經(jīng)整理可得：f*Cn2+dn2)r2dr0徵+I2實(shí)際上由Mie散射理論可知，上式中的積分項(xiàng)為電場強(qiáng)度的平方對角度8、小全空間積分的平均值，即：于是吸收效率為。心=eIJ）工”式中x'=rk=z

17、/m。當(dāng)xn1時即瑞利散射情況，顆粒的內(nèi)部平土§場強(qiáng)為常數(shù)，其值為:ImprovedMiescatteringalgorithmsW.J.WiscombeMie計(jì)算存在的問題就是如何最有效地構(gòu)造Mie計(jì)算，同時保證準(zhǔn)確性和避免數(shù)值的不穩(wěn)定性和病態(tài)。Mie計(jì)算以耗時著稱，首先無窮項(xiàng)級數(shù)N的求和，例如:100Rm的水滴在0.5Nm的可見光散射情況下，大約需1260項(xiàng)求和。其次，典型的計(jì)算都希望能對一系列半徑（如對尺寸分布求積分）、一系列波長（如對太陽光譜求積分）及一系列折射率求和（如通過散射參量反推折射率）。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)當(dāng)折射率虛部mim很大時，用向后循環(huán)法求An很不穩(wěn)定。而向前遞推總是穩(wěn)定的（但向后遞推安全時，總是優(yōu)先選擇，因?yàn)槠溆?jì)算速度很快）。得出允許向后遞推的經(jīng)驗(yàn)標(biāo)準(zhǔn)：用正確的向前地推與相對應(yīng)的向后地推做比較，當(dāng)發(fā)現(xiàn)對Qsch和g的相對誤差超過10-6時，認(rèn)為計(jì)算失敗。對于一對確定的（x,mRe）,我們采用向后遞推尋找第一個循環(huán)失敗的mfm研究表明：對于確定的41Rs,工的值隨著x的增加很快趨向于一個確定值。xmfm&g

19、t;minCcmim)allx.8+26.22191)0.4474m也+0.00204喊日-0.000175mM.£4002K02000ieoo*400S4CT）4000BOO2coI£34S6TS如果在任意角度下6、s2的實(shí)部和虛部的相對誤差超過105時，認(rèn)為對S1和S2的向后遞推失敗。（而此時，QscaQxt并不受影響，因?yàn)楫?dāng)S1,$2的相對誤差達(dá)到10V時，QscaQext的相對誤差總維持在10-10以下。）scadxi又tS1和S2=13,78nI410,8mr*+39=16,35mL+15.04,對散射強(qiáng)度|Sip+|S2p和偏正度(展-|Si|2)/(|S/+S

20、i|Z)連分式算法總結(jié)：Mie散射計(jì)算的核心是計(jì)算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分別是第一和二類貝塞耳函數(shù)，a稱為當(dāng)量直徑，a=2兀r/入，r是球形顆粒的真實(shí)半徑，入是入射光的波長，m為折射率中網(wǎng)=(H2)sj*L)，=/2)叫加+J1)以11AFDn(p)=ln0fp>=i.dpWJP)式中p為函數(shù)任一自變量。貝塞耳函數(shù)遞推關(guān)系式：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散射計(jì)算中Jn、Yn、Dn的計(jì)算是關(guān)鍵和難點(diǎn)。對于Dn，我們采用的是Lentz的連分式的算法D“旅矽=.能Lentz證明有如下關(guān)系fma);f©G7faka、J鞭(ma)a1工程電/其中,1廠一»"+卜05.人我們注意到當(dāng)上78時,"不。所以可以利用上式累積相乘直到滿足精度要求。(可根據(jù)精度要求例如10-7來確定所要達(dá)到的k值)對于Jn、Yn的生成本文也采用連分式的算法。具體方案如下令Cn=Jn-1

22、(a)/Jn(a)，根據(jù)貝塞耳差積公式:JJRY人"勾J八"刃Y"幻=才？由以上二式整理得:J八1/C-fa)Y相ifa)餐)Y”®=JW上式中Cn的計(jì)算是采用類似于Dn的連分式的形式，計(jì)算中可調(diào)用同一函數(shù)計(jì)算。若已知初值：Jo=sdna/a；a)=sin/后cos/cm；,v2.Yq=-cos/a;Yi(=-cos/)-sin度/四；這樣就可計(jì)算出各級Jn和YnW川iamJ.Lentz關(guān)于連分式的文章:九=（打2/2產(chǎn)心龍幺Q）=（值/2）啊心.姑+（-1）丐上心使W/）-2/其中n+%=%N=m工。以ArG）為基礎(chǔ)，采用貝塞爾函數(shù)比值的連分式表示法：

23、式力,利用此法可產(chǎn)生所有的4，工匕），盡管耗時，但能減少存儲需求。同時可通過計(jì)算高階值，使用下面的遞推公式，從后往前算出其他值。兀卜）_2"_j（z）1g=；Q）蕓jrr-lQ）1+-2（+1U-1箝=2K不像一般的函數(shù)，貝塞爾函數(shù)的比值一旦超過可控制的邊界，就不再增長，初始的高階出工（胃）值決定了所有低階值的準(zhǔn)確性，因此，采用新方法計(jì)算準(zhǔn)確的初始比值是必要的。+2（+2）源1-2B+3）父1+.處于分母位置的+號表示分母上加上一個特殊的連分式。類似于上式/中的表示形式。定義一種新的符號：n1,111/=,Gi+I+也Lentz給出了n階部分收斂值為：例如：實(shí)變量=L0,虛數(shù)爐=9.

24、5:計(jì)算過程:）1-1,2,3,.事%+工=-20.94736842%工（T產(chǎn)以廿十言一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,米散射學(xué)習(xí)目前所遇到的困難：到底怎樣的計(jì)算結(jié)果才算正確，如何能找到一個米散射計(jì)算結(jié)果準(zhǔn)確又有效的數(shù)據(jù)庫，來驗(yàn)證自己算法及程序的正確性。倒退式算法的總結(jié)：由于Mie級數(shù)的收斂速度隨特征值X的增大而減慢，因此不同的萬值即使在相同的計(jì)

25、算精度要求下所需級數(shù)項(xiàng)數(shù)M也不一樣.1979年美國學(xué)者WJWiscomt/s】在大量計(jì)算的基礎(chǔ)上，參考前人的工作總結(jié)出一個X；的經(jīng)驗(yàn)計(jì)算公式.利用該公式給出的項(xiàng)數(shù)可以使整個計(jì)算誤差小于10.Wiscombe公式是工+4立3+Ix+4.05x3+10.02w工S88<x<42001工+4%5+24200三#W20000當(dāng)X小于0.02時可以利用lluykigh公式計(jì)算,而當(dāng)X比20000更大時.幾何光學(xué)即可適用.Dn的計(jì)算采用Dave的倒推式：八/b-ikn/tncl4Dn（/槨）£>,1（如=ma由于Dn函數(shù)有很強(qiáng)的收斂性，對于Dn的倒推計(jì)算的初值的選取有很強(qiáng)的隨意性。因?yàn)楫?dāng)n-o"Dn(ma)一0，所