A.-Mie米散射理论基础.docVIP

A.-Mie米散射理论基础

PAGE

PAGE17

面波散射的精确解。一般把粒子直径与入射光波长相当的微粒子所造成的散射称为米散射。米散射适合于任何粒子尺度,只是当粒子直径相对于波长而言很小时利用瑞利散射、很大时利用夫琅和费衍射理论就可以很方便的近似解决问题。米散射理论最早是由G1Mie在研究胶体金属粒子的散射时建立的。

1908年,米氏通过电磁波的麦克斯韦方程,解出了一个关于光散射的严格解,得出了任意直径、任意成分的均匀粒子的散射规律,这就是著名的米氏理论[4-6]。根据米散射理论,当入射光强为I0,粒子周围介质中波长为λ的自然光平行入射到直径为D的各向同性真球形粒子上时,在散射角为θ,距离粒子r处的散射光和散射系数分别为:

从上式中可以看到,因为是各向同性的粒子,散射光强的分布和φ角无关。同时,上式中:

i1、i2为散射光的强度函数;s1、s2称为散射光的振幅函数;a为粒子的尺寸参数(a=πD/λ);m=m1+im2为粒子相对周围介质的折射率,当虚部不为零时,表示粒子有吸收。对于散射光的振幅函数,有:

式中an、bn为米散射系数,其表达式为:

其中:

是半奇阶的第一类贝塞尔函数;是第二类汉克尔函数;Pn(cosθ)是第一类勒让德函数;P(1)n(cosθ)是第一类缔合勒让德函数。

Mie散射理论

Mie散射理论是麦克斯韦方程对处在均匀介质中的均匀颗粒在平面单色波照射下的严格数学解。由Mie散射知道,距离散射体r处p点的散射光强为

式中:为光波波长;I0为入射光强;Isca为散射光强;为散射角;为偏振光的偏振角。

式中:和是振幅函数;an和bn是与贝塞尔函数和汉克尔函数有关的函数;和是连带勒让得函数的函数,仅与散射角有关。其中

式中:和分别是贝塞尔函数和第一类汉克尔函数;和是和的导数;为无因次直径,,D为颗粒的实际直径;是入射光的波长;m是散射颗粒相对于周围介质的折射率,它是一个复数,虚部是颗粒对光的吸收的量化。由以上公式可见,Mie散射计算的关键是振幅函数和,它们是一个无穷求和的过程,理论上无法计算。求解振幅函数的关键是计算an和bn,所以Mie散射的计算难点是求解an和bn。

Mie散射理论的数值计算

通过以上分析可知,Mie散射计算的核心是求解an和bn,我们编制程序也是围绕它进行编写。在an和bn的表达式中,,和

满足下列递推关系:

这些函数的初始值为;

与散射角有关的和满足下列递推公式:

有了这些递推公式可以很方便地通过计算机程序求解。但是对于n的大小,因为计算机不可能计算无穷个数据,所以n在计算之前就要被确定。

散射理论基础与Matlab实现

若散射体为均匀球体,如图1所示,照射光为线偏振平面波,振幅为E,光强I0,沿z轴传播,其电场矢量沿x轴振动。散射体位于坐标原点O,P为观测点。散射光方向(OP方向)与照射光方向(z轴)所组成的平面称为散射面,照射光方向至散射光方向之间的夹角θ称为散射角,而x轴至OP在xy平面上投影线(OP′)之间的夹角φ称为极化角。观测点与散射体相距r。根据经典的Mie散射理论,散射粒子的尺度参数为α=2πa/λ,其中a为球形粒子的半径,散射粒子相对周围介质的折射率为m=m1+i*m2。则散射光垂直于散射面和平行于散射面的两个分量的振幅函数为:

以上式中：

Jn+1/2(z)和Yn+1/2(z)分别为半整数阶的第一类,第二类贝塞尔函数。

P(1)n(cosθ)为一阶n次第一类缔合勒让德函数;Pn(cosθ)为第一类勒

让德函数。在数值模拟过程中选取初始下:

微粒子对光的散射和吸收是电磁波与微粒子相互作用的重要特征,而微粒对电磁辐射的吸收与散射与粒子的线度有密切关系,对于不同线度的粒子必须应用不同的散射理论。Mie散射理论主要用于从亚微米至微米的尺寸段;在微米以下至纳米的光散射则近似为形式更明晰简单的瑞利散射定律,散射光强烈依赖于光波长λ(I～λ-4);而对大于微米至毫米的大粒子则近似为意义明确的夫朗和费衍射规律了。

Mie散射理论给出了球型粒子在远场条件下的散射场振幅an、bn以及粒子内部电磁场振幅cn、dn的计算表达式,通常称为Mie散射系数

式中m表示微粒子外部介质的相对折射率,x=κa,a为球的半径,κ=2π/λ称为波数,μ为相对磁导率,即球的磁导率与