李顺 李正军† 屈檀 李海英 吴振森

1) (西安电子科技大学物理学院,西安 710071)

2) (西安电子科技大学,信息感知技术协同创新中心,西安 710071)

3) (西安电子科技大学电子工程学院,西安 710071)

基于广义洛伦兹Mie 理论,研究了单轴各向异性球形粒子对两束具有任意传播和极化方向的零阶贝塞尔波束的传播和散射特性,并与单零阶贝赛尔波束入射单轴各向异性球形粒子时的传播和散射特性进行了对比研究.利用球矢量波函数的正交关系及坐标旋转定理,导出了任意传播和极化方向零阶贝塞尔波束的球矢量波函数的展开形式,通过矢量叠加得到了总入射场的展开系数.基于傅里叶变换方法和切向连续的边界条件,得到了单轴各向异性球内部电磁场的球矢量波函数展开式,并导出了散射系数解析表达式.将零阶贝塞尔波束退化成平面波,通过将其入射到单轴各向异性球形粒子的雷达散射截面角分布与文献结果进行对比,验证了本文理论及程序的正确性.数值分析了入射角、锥角及极化角等参数对雷达散射截面角分布的影响.本文理论和数值结果希望能应用于多波束入射下各向异性粒子、生物细胞等复杂粒子体系的散射、粒径分析以及光学俘获等特性的研究中.

1 引言

自从Durnin[1]提出贝塞尔波束以来,由于其非衍射和自重构特性,受到越来越多的关注,并广泛应用于不同的领域,如光捕获和操作、粒径分析、光导和对准等[2-5].作为激光束,许多学者对贝塞尔波束在球面和柱面坐标系中的轮廓和描述进行了深入研究[6-11].基于该描述,许多学者进一步研究了球形粒子对贝塞尔波束的散射问题.Marson[12]使用平面波分解方法研究了球体对零阶贝塞尔的散射特性,Ma 和Li[13]研究了球形粒子对非偏振贝塞尔波束的散射.入射贝塞尔波束被认为是Cimár等[8]提出的平面波的叠加,并在球面矢量波函数方面进行了展开,实际上它是轴棱镜生成的贝塞尔波束.基于广义洛伦兹Mie 理论[14],Li 等[15]研究了由棱镜生成的贝塞尔波束照射的球形粒子的散射.用积分局域近似法(integral local approximation,ILA),Ambrosio 和Herández-Figueroa[10]给出了零阶贝塞尔波束的球矢量波函数展开系数.然而,贝塞尔波束由标量波理论描述,只有当圆锥角比较小时才能提供令人满意的结果.Mishra[11]则利用矢量波理论来描述任意圆锥角的贝塞尔波束,得到了比较满意无衍射特性的贝塞尔波束.使用这种矢量描述,Mitri[16]研究了各向同性球形粒子对零阶贝塞尔波束的散射特性.利用数值方法,Cui 等[17]研究了任意形状的均匀介电粒子对零阶贝塞尔波束的散射特性.Klimov[18]导出了零阶贝塞尔波束在球体上散射的解析解,Wang 等[19]研究了零阶贝塞尔波束激发的受激拉曼散射显微成像用于散射组织中深层目标的成像.但是上述研究主要针对简单的各向同向粒子的散射问题,关于各向异性粒子对零阶贝塞尔波束散射特性的研究还比较少.

近年来,各向异性介质的散射特性因其在光信号处理、雷达散射截面(radar cross section,RCS)控制和微波器件制造等领域的应用越来越广泛而备受关注.Stout 等[20]使用微分理论给出了由任意各向异性介质构成的任意形状物体的散射解.Wong等[21]和Qiu 等[22]研究了单轴各向异性(uniaxial anisotropy,UA)球形粒子对平面波的散射特性.通过引入傅里叶变换的方法,Geng 等[23]利用解析方法研究了UA 球形粒子对平面波的散射特性,Wang 等[24]研究了旋转各向异性介质球对离轴高斯波束的散射特性.近十年,Yuan 等[25,26]深入研究了在轴、离轴和任意入射高斯波束照射下UA 球体的散射特性.Wang 等[27]基于T矩阵方法研究了UA 介质球对高斯波束的散射特性,并与解析方法进行了比较验证.Qu 等[28]研究了UA 球体对离轴零阶贝塞尔波束的散射特性.

然而这些研究中,入射波仅限于平面波或单个贝塞尔波束,其中各向异性球体对双零阶贝塞尔波束散射特性的研究还是非常少的.特别地,由于主光轴的存在,零阶贝塞尔波束离轴斜入射到UA 粒子上时,其散射特性与在轴和离轴入射情况下UA粒子的散射特性有很大不同.并且,与单波束光势阱相比,双波束光势阱或驻波光势阱用于非常小的同位素球形粒子及非各向同性粒子的光学捕获更加具有优势[29].这就要求对于非各向同性球形粒子对双聚焦波束的散射理论及特性问题要研究透彻.基于此,本文主要研究了UA 球形粒子对于双任意方向入射零阶贝塞尔波束的散射问题,并详细讨论了双零阶贝塞尔波束的电场强度分布及作用于UA 球形粒子RCS 角分布的情况.值得注意的是,本文的研究与广义洛伦兹Mie 理论不同的地方主要有两点: 一是入射波束从高斯波束入射扩展到任意方向传播双零阶贝塞尔波束入射,其球矢量波函数的展开利用了坐标旋转定理,展开系数是通过积分方法及矢量加法定理得到的;二是内场的展开要比各向同性球形粒子要复杂得多,基于傅里叶变换方法根据球矢量波函数L,M和N的正交完备性进行展开.因而,UA 球体对于双任意方向入射零阶贝塞尔波束散射问题的研究是广义洛伦兹Mie 理论的深入的扩展.

2 散射理论

2.1 零阶贝塞尔波束的描述

作为一种典型的有形波束,零阶贝塞尔波束的标量形式可以看作是标量波动方程的一类精确解.假设时谐因子为 e-iωt(ω为角频率),零阶贝塞尔波束电场在波束坐标系Oxyz中可表示为[1]

其中E0为归一化振幅,J0(·) 为零阶柱阶贝塞尔函数,参数R=以及ϕ=tan-1(y/x) 分别为横向面平面 (x,y)上的半径和方位角,kR=k0sinC0和kz=k0cosC0分别为波矢k的横向分量和纵向分量,C0为波束的半锥角.

一般来说,当半锥角小于10°时,上述标量表示(1)式足以精确地描述零阶贝塞尔波束的无衍射特性.但是,当半锥角比较大时: 例如kR≈k0,(1)式所表示的贝塞尔波束就会与实际情况出现比较大的误差.因此,Mishra[11]通过两个矢势来构建一个电磁场,得到具有对称性的电场强度和磁场强度的表达式.矢量描述是通过引入矢量势A来定义零阶贝塞尔波束的电磁场[30]:

假设矢势A沿x轴方向偏振,通过求解亥姆霍兹方程,可以得到

对于矢量波方法,为得到圆对称的零阶贝塞尔波束,设沿y方向偏振的矢势A′,则

E′和H′的表达式为

(5)式与(2)式结果相加再取平均,即可得到沿z轴传播、x轴偏振的零阶贝塞尔波束在其波束坐标系Oxyz中的电磁场表示式[11,16]:

2.2 双零阶贝塞尔波束的展开

如图1(a)所示,两束零阶贝塞尔波束离轴斜入射到UA 介质球上,两波束中心分别为O1和O2.O1x1y1z1和O2x2y2z2分别是以各自波束中心为原点,传播方向为z1和z2轴建立的直角坐标系.Oxyz是以UA 介质球心O为原点建立的直角坐标系,设定各向异性介质球的主光轴沿z轴方向,半径为a.设定在粒子坐标系Oxyz中,两波束中心O1和O2的坐标分别为 (x1,y1,z1) 和 (x2,y2,z2),两束零阶贝塞尔波束的半锥角分别用C1和C2表示,传播方向分别用k1和k2表示.如图1(b)所示,设定α1为第一束零阶贝塞尔波束的入射角,是波束传播方向k1与粒子坐标系Oxyz中z轴的夹角;β1为第一束零阶贝塞尔波束的极化角,是波束传播方向k1在粒子坐标系Oxyz下xOy面上的投影与x轴的夹角,这样可以表示出零阶贝塞尔波束的传播和极化方向的任意性.类似地,可以设定α2和β2为第二束零阶贝塞尔波束的入射角和极化角.这里与图1(b)的情况相同,只是角度的表示不一样,不再重复画出示意图.值得注意的是,由于UA 介质存在主光轴,当波束的入射方向与主光轴不同时会产生不同的散射特性.入射角定义为0°—180°的范围,相对各向异性球形粒子的主光轴z轴而言,零阶贝塞尔波束是以任意方向入射UA 介质球的.

图1 双零阶贝塞尔波束离轴斜入射UA 介质球的示意图Fig.1.Schematic diagram of a uniaxial anisotropic sphere illuminated by off-axis obliquely incident double zero-order Bessel beams.

将(6)式和(7)式中相应的坐标表示 (x,y,z) 分别替代成坐标表示 (x1,y1,z1) 和 (x2,y2,z2),即可以得到两束零阶贝塞尔波束在各自坐标系O1x1y1z1和O2x2y2z2的电磁场矢量表示式.由于零阶贝塞尔波束的矢量表示式本身就是从麦克斯韦方程组推导得来的,因此可以将这两束零阶贝塞尔波束根据球矢量波函数进行展开.以第一束零阶贝塞尔波束为例,可根据球矢量波函数在粒子坐标系Oxyz下展开:

其中,ε0和µ0分别为在真空传播时的介电常数和磁导率,ψn(k0a)=k0ajn(k0a) 为黎卡提贝塞尔函数.如果已知原始零阶贝塞尔波束在粒子坐标系Oxyz下的径向电磁场,则可以通过求解方程(11)式中的二重积分来得到展开系数.需要注意的是,零阶贝塞尔波束的入射方向不平行于z轴,因此无法通过简单的坐标平移得到波束的径向电磁场.可以通过坐标系的旋转理论来解决这个问题.如图2 所示,以第一束零阶贝塞尔波束中心O1为原点建立与粒子坐标系Oxyz相互平行的中间直角坐标系O1x10y10z10.一般地,任意两个非平行的直角坐标系都可以通过欧勒角α,β,γ角相互旋转得到[34,35].为了简化,这里只用两个角即α1,β1.首先将直角坐标系O1x1y1z1绕y1轴逆时针旋转α1角,得 到z1与z10轴重合;再将O1x1y1z1以z1为轴逆时针旋转β1角,即可以使直角坐标系O1x1y1z1与O1x10y10z10相互重合.

通过图2 的旋转关系,可以得到两个坐标系下坐标的关系:

图2 两直角坐标的旋转关系示意图Fig.2.Schematic diagram of rotation relationship between two rectangular coordinate systems.

由于中间坐标系O1x10y10z10与粒子坐标系Oxyz的坐标之间有如下关系:

因此可以得到第一束零阶贝塞尔波束在粒子坐标系下的电磁场径向分量,再代入(11)式即可以获得第一束零阶贝塞尔波束在粒子坐标系Oxyz下的展开系数.

类似地,建立一个平行于粒子坐标系Oxyz的中间坐标系O2x20y20z20,第二束零阶贝塞尔波束也可以根据球矢量波函数在粒子坐标系Oxyz下展开:

同样可以得到第二束零阶贝塞尔波束在粒子坐标系下的展开系数:

对于两束入射的零阶贝塞尔波束,可通过两个电磁场的矢量叠加来计算总入射场,则总入射场展开系数为

2.3 双零阶贝塞尔波束的传播特性

图3 为不同半锥角下沿z轴正、负方向相同极化角时,反向传播的双零阶贝塞尔波束在xOy面的电场强度分布情况.零阶贝塞尔波束电场振幅E0设定为1,计算参数为:α1=β1=β2=0°,α2=180°,z1=1 μm,z2=-1 μm (后文若无特殊说明,E0=1,z1与z2保持不变).与单束零阶贝塞尔波束相比,反向传播双零阶贝塞尔波束在xOy面的电场强度分布形成了4 个明显的孤岛.同时,半锥角越大,反向传播双零阶贝塞尔波束在xOy的场值整体越强,4 个孤岛更加向中心区域集中,周边区域的电场强度分布也会形成更多的孤岛.

如图4 所示,计算了不同半锥角下沿z轴正、负方向不同极化角时,反向传播的双零阶贝塞尔波束在xOy面的电场强度分布情况.计算参数为:α1=β1=0°,α2=β2=180°.与极化角相同时反向传播双零阶贝塞尔波束(图3)相比,形成的孤岛有所减少,而孤岛中心强度却明显增加,形成孤岛的位置和方位也发生了较大的改变.在与粒子相互作用中,两束零阶贝塞尔波束极化角不同时,散射特性应该会不同.

图5 给出了不同半锥角下斜入射反向传播的双零阶贝塞尔波束在xOy面的电场强度分布情况.计算参数为:α1=30°,α2=210°,β1=β2=0°.虽然双零阶贝塞尔波束的传播方向是反向的,但是电场强度分布并没有形成孤岛.因为斜入射情况下,双波束相向传播到波束中心时,应该在与其传播方向相互垂直的面上形成孤岛而不是xOy面上.随着半锥角的增大,双零阶贝塞尔波束在中心区域的电场渐渐集中,强度增强,与图3 和图4 所示的电场强度分布非常不同.

图3 相同极化角下反向传播双零阶贝塞尔波束的电场强度分布 (a) C1=C2=5°;(b) C1=C2=15°;(c)C1=C2=30°Fig.3.Electric intensity distribution of back propagating double zero-order Bessel beams with identical polarization angles:(a) C1=C2=5°;(b) C1=C2=15°;(c) C1=C2=30°.

图4 不同极化角下反向传播双零阶贝塞尔波束的电场强度分布 (a) C1=C2=5°;(b) C1=C2=15°;(c)C1=C2=30°Fig.4.Electric intensity distribution of back propagating double zero-order Bessel beams with different polarization angles:(a) C1=C2=5°;(b) C1=C2=15°;(c) C1=C2=30°.

图6 展示了不同半锥角下斜入射非反向传播的双零阶贝塞尔波束在xOy面的电场强度分布情况.计算参数为:α1=30°,α2=120°,β1=β2=0°.由于双零阶贝塞尔波束的传播方向不是反向,电场强度在xOy面的分布没有形成较为明显的孤岛,这与非反向传播的波束并不会形成驻波的事实相符合.随着半锥角的增大,双零阶贝塞尔波束中心区域的电场渐渐分散,与图5 所示的电场强度分布刚好相反.与图3—图5 相比,图6 所展示的非反向传播的双零阶贝塞尔波束场强分布与反向传播的双零阶贝塞尔波束场强分布有着明显的区别,这会导致入射粒子的散射特性的不同.

图5 斜入射时反向传播双零阶贝塞尔波束的电场强度分布 (a) C1=C2=5°;(b) C1=C2=15°;(c)C1=C2=30°Fig.5.Electric intensity distribution of back propagating double zero-order Bessel beams with oblique incidence: (a) C1=C2=5° ;(b) C1=C2=15°;(c) C1=C2=30°.

图6 斜入射时非反向传播双零阶贝塞尔波束的电场强度分布 (a) C1=C2=5°;(b) C1=C2=15°;(c)C1=C2=30°Fig.6.Electric intensity distribution of non-back propagating double zero-order Bessel beams with oblique incidence: (a)C1=C2=5°;(b) C1=C2=15°;(c) C1=C2=30°.

2.4 UA 介质球的内场展开

对于半径为a的UA 介质球,其介电常数和磁导率分别表示为

在无源情况下,由UA 介质的麦克斯韦方程组可以得到电场矢量波动方程:

(19)式即是无源UA 介质中电场矢量所满足的微分波动方程.该波动方程由于是张量,电场的各个分量耦合在一起,无法用分离变量法来求解,可以引入傅里叶变换来求解[23]:

通过引入傅里叶变换解本征方程,UA 介质球的内部电磁场可以用球矢量波函数表示为[36]

2.5 通过边界条件求解散系数

在r=a处,由电磁场在边界上切向分量连续有

其中t表示电场磁场的切向分量.将入射场、内场及散射场的展开详细表达式代入边界条件,于是在r=a处,可得

求解上述方程组(27)—(30),消去散射系数可得

其中,

化简后可以得到关于内场展开系数Gmr′q的方程组:

联立方程组求解内场展开系数Gmr′q将其代入(27)式和(30)式,可求解出散射系数的表达式为[37]

将散射系数代入(23)式和(24)式中得到散射场的电场强度.由RCS 的定义可得

3 数值结果和讨论

当半锥角为0 时,零阶贝塞尔波束将退化为平面波形式.为了验证程序的正确性,我们将单零阶贝塞尔波束退化成平面波,计算了入射UA 介质球RCS 的角分布,并且与文献[38]给出的结果进行了比较.如图7 所示,发现二者结果符合较好,验证了程序与理论公式推导的正确性.图7 使用的参数为:εt=2.5ε0,εz=1.0ε0,µt=µz=1.0µ0,λ=1.064 μm,α=β=0°,a=1 μm,C0=0°.Eplane 对应xOz面,称为E面;H-plane 对应yOz面,称为H面.

图7 UA 介质球对单零阶贝塞尔波束的RCS 角分布Fig.7.Normalized RCS of a UA dielectric sphere illuminated by single zero-order Bessel beam versus the scattering angle.

图8 计算了反向传输的双零阶贝塞尔波束退化成平面波入射到UA 球形粒子的RCS 的角分布,并与文献[39]结果进行了对比,发现两者也符合得比较好,进一步验证了本文理论和程序的正确性.计算过程中所使用的参数为εt=2.5ε0,εz=1.0ε0,µt=µz=1.0µ0,λ=1.064 μm,a=1 μm,α1=β1=β2=0°,α2=180°,x1=y1=x2=y2=0,z1=1 μm,z2=-1 μm.比较图7 和图8 可以看出,当双零阶贝塞尔波束照射UA 球形粒子时,产生的散射效应是增加的.双零阶贝塞尔波束入射情况下,由于波束在粒子中存在干涉效应,所以E面和H面的峰值也不是特别明显.

图8 双零阶贝塞尔波束退化成平面波入射UA 介质球的RCS 角分布 (a) E-plane;(b) H-planeFig.8.Angular distribution of the RCS of a UA dielectric sphere by double zero-order Bessel beams when degenerating into plane waves: (a) E-plane;(b) H-plane.

如图9 所示,计算了UA 介质球对不同离轴情形单零阶贝塞尔波束的RCS 角分布,计算参数为:λ=1.064 μm,εt=5.3495,εz=4.9284,µt=µz=1,n0=1.0,α1=β1=0°,C1=C1=10°.当入射波束沿x轴方向偏移时,E面和H面上的RCS 值都在减小,E面出现了不对称的情况,但是H面上的RCS 的角分布对称性和形状都没有太大改变.当单波束沿y轴方向偏移时,E面和H面上的RCS值也都在减小,与沿x轴方向偏移时的情形不同,此时E面上的RCS 角分布的对称性和形状都没有太大的改变,但是在H面出现不对称的情况.之所以出现RCS 值减小的情况,是因为波束离轴入射粒子,波束中心距离球心横向位置越远时,波束照射球体的强度会越小.而沿x轴和y轴偏移会出现不一样的情况,是因为方位角为ϕ=arctan(y/x),x和y的改变会影响到散射系数.

图9 UA 介质球对不同离轴情形入射的单零阶贝塞尔波束的RCS 角分布 (a) E-plane;(b) H-planeFig.9.Angular distribution of the RCS of a UA dielectric sphere illuminated by different off-axis incident zero-order Bessel beams:(a) E-plane;(b) H-plane.

如图10 所示,计算了UA 介质球对不同离轴双零阶贝塞尔波束的RCS 角分布,计算参数为:λ=1.064 μm,εt=5.3495,εz=4.9284,µt=µz=1,n0=1.0,α1=β1=β2=0°,α2=180°,C1=C2=10°.可以明显地看到,无论波束是沿x轴方向还是沿y轴方向偏离,RCS 值都会减小.这是因为离轴入射时波束中心偏离球心的横向位置越远,波束照射到球体上的强度就会越小,这与单波束的情况类似.与图9 所示单离轴零阶贝塞尔波束入射情况相比,图10 中RCS 角分布发生了明显地改变,尤其是双零阶贝塞尔波束沿正y轴和负y轴离轴入射UA 介质球时,RCS 出现了两个极小值点.

图10 UA 介质球对离轴入射的双零阶贝塞尔波束的RCS 角分布 (a) E-plane;(b) H-planeFig.10.Angular distribution of the RCS of a UA dielectric sphere illuminated by different off-axis incident double zero-order Bessel beams: (a) E-plane;(b) H-plane.

图11 计算了第二束零阶贝塞尔波束入射角与第一束零阶贝塞尔波束的不同时,双波束非反向入射UA 球形粒子的RCS 角分布,所用参数为:λ=1.064 μm,εt=5.3495,εz=4.9284,µt=µz=1.0µ0,n0=1.0,x1=y1=z1=0,x2=y2=z2=0,α1=β1=β2=0°,C1=C2=20°.可以看出,随着第二束零阶贝塞尔波束入射角不断的变化而产生的RCS 角分布是不同的.当入射角都是0°时,两束零阶贝塞尔波束的传播方向相同,此时RCS 值最大.但是两束波束由于在粒子之间有相互作用的关系,所得的RCS 值并不是单零阶贝塞尔波束入射UA球形粒子时RCS 的2 倍.随着第二束零阶贝塞尔波束的入射角不断增大,即第二束零阶贝塞尔波束是斜入射时,在E面上峰值没有正入射时大,并且峰值出现的次数会随着入射角的不断增大而增加.这说明在双波束的入射角度变大时,其相互干涉的效应不断增大.而在H面上,由于双零阶贝塞尔波束的极化方向垂直于H面,所以H面上的RCS 角分布总是对称的.

图11 UA 介质球对不同入射角下双零阶贝塞尔波束的RCS 角分布 (a) E-plane;(b) H-planeFig.11.Angular distribution of the RCS of a UA dielectric sphere illuminated by double zero-order Bessel beams with different incident angles: (a) E-plane;(b) H-plane.

图12 计算了双反向传输的零阶贝塞尔波束在不同锥角入射的情况下UA 介质球的RCS 角分布,参数如下:λ=1.064 μm,a=0.5 μm,εt=5.3495,εz=4.9284,µt=µz=1.0µ0,n0=1.0,x1=y1=x2=y2=0,z1=z2=1 μm,α1=β1=β2=0°,α2=180°.可以看出,随着零阶贝塞尔波束锥角的逐渐增大,在E面和H面的RCS 值都是逐渐减小.这是因为随着锥角的增大,波束之间的干涉效应在逐渐增加.当两束波束都沿z轴传播时,所得到的RCS 值是对称的.

图12 双反向传输的零阶贝塞尔波束在不同锥角入射的情况下UA 介质球的RCS 角分布 (a) E-plane;(b) H-planeFig.12.Angular distribution of the RCS of a UA dielectric sphere illuminated by back propagating double zero-order Bessel beams with different conic angles: (a) E-plane;(b) H-plane.

4 结论

基于广义洛伦兹Mie 理论,研究了UA 球形粒子对两束具有任意传播和极化方向的零阶贝塞尔波束的传播和散射特性.利用坐标旋转关系及特殊函数的正交关系,导出了任意传播和极化方向的双零阶贝塞尔波束的球矢量波函数的展开形式及展开系数的解析表达式,详细讨论了单、双零阶贝塞尔波束的电场强度分布特性.通过将单、双零阶贝塞尔波束退化成平面波的特殊情况下的RCS 角分布与文献结果进行对比,验证了本文理论及程序的正确性.数值分析了各参数对双零阶贝塞尔波束入射UA 球形粒子的RCS 角分布的影响,并与单零阶贝塞尔波束入射UA 球形粒子的散射特性进行对比.由于干涉效应,当粒子被双零阶贝塞尔波束照射时,RCS 值并不是单零阶贝塞尔波束入射时的两倍.且第二束零阶贝塞尔波束在入射角增大的时候干涉效应逐渐增强,说明使用双波束与各向异性介质产生的相互作用会更强.本文理论和数值结果希望能为多波束入射下非各向同性、生物细胞等复杂粒子的散射特性、粒径分析以及光学微操控等提供理论依据.