J0相关径向偏振涡旋光束紧聚焦的空间相干特性

饶连周

(1.三明学院 机电工程学院,三明 365004；2.三明学院 福建省装备智能控制重点实验室，三明 365004)

引 言

众所周知，由于涡旋光束的独有特性，人们对偏振涡旋光束通过大数值孔径的紧聚焦特性进行了广泛的研究[1-7]。特别是紧聚焦的涡旋光束会产生波长尺度的光斑或暗核，这些结果可以应用于许多领域。近年来，部分相干光通过大数值孔径的紧聚焦特性也引起了极大的研究兴趣[8-20]。例如，FISCHER和VISSER[10]研究了部分相干紧聚焦光场的光谱相干特性；GUO等人[11]研究了部分相干径向偏振光涡旋光束紧聚焦的空间相干特性；还有一些研究人员[15-17]分析了J0相关偏振光束紧聚焦的传输特性等等。然而，通过大数值孔径的J0相关径向偏振光束的空间相干特性鲜见相关研究报道。本文中，将基于矢量德拜衍射理论，研究J0相关径向偏振光束的紧聚焦场，分析焦点区域光场的光谱相干度。结果表明，焦点区域相干度不仅与数值孔径大小和相干参量有关，而且还与拓扑电荷数有关。特别是，相干度还具有相位异常点。

1 J0相关径向偏振光束的紧聚焦

一个完全相干且径向偏振的光束，通过大数值孔径聚焦后，在焦点区域的电场可以表示为[21]:

exp[ik(zcosθ+rsinθcos(φ-ψ)]dφdθ

(1)

式中，r,ψ和z是焦点区域中的径向、角向和纵向坐标,φ为入射光束的方位角,波数k=2π/λ,A是常数，α=arcsin(dNA)是由透镜的数值孔径(numerical aperture,NA)决定的最大角度，θ是入射光束的会聚角,并从0变化到α，B(θ,φ)是入射场的切趾函数。

部分相干涡旋光束的相干特性可以用交叉谱密度矩阵来描述，矩阵元素的表达式如下[22]:

Wij(r1,r2)=〈Ei*(r1,ψ1,z1)×

Ej(r2,ψ2,z2)〉,(i,j=x,y,z)

(2)

式中，Ei*(r1,ψ1,z1)和Ej(r2,ψ2,z1)表示两点处电场的笛卡尔分量，*表示复共轭，〈〉表示系综平均值。将(1)式代入(2)式，焦点区域交叉密度矩阵的元素就可以表示为:

exp[-ikr1sinθ1cos(φ1-ψ1)]exp[ikr2sinθ2×

cos(φ2-ψ2)]Pi(θ1,φ1)Pj(θ2,φ2)×

exp[ik(z2cosθ2-z1cosθ1)]dθ1dθ2dφ1dφ2,

(i,j=x,y,z)

(3)

(4)

式中,B0(θ1,φ1,θ2,φ2)=B*(θ1,φ1)B(θ2,φ2)称为入射光束的切趾相干函数。假设通过大数值孔径聚焦的光束是一个J0相关高斯-谢尔模型涡旋光束，这种光束的切趾相干函数可以表示为[17，23]:

B0(θ1,φ1,θ2,φ2)=B*(θ1,φ1)B(θ2,φ2)=

Jm(dNA-1δsinθ2)exp[i(n-m)(φ2-φ1)]

(5)

式中,σ=w0/d表示光束的相对光斑大小，w0为束腰半径,d为大数值孔径的半径，δ=βd为相干参量,n和m为序数，n>m。

光谱相干度是描述紧聚焦光束特性的重要指标，其定义式[10,17]可以表示为:

(6)

式中,Tr表示矩阵对角元素之和。将(4)式和(5)式代入(3)式，即可得到焦区焦点区域的交叉谱密度Wij(r1,ψ1,z1;r2,ψ2,z2)。通过设置r1=r2=r，ψ1=ψ2=ψ，z1=z2=z，就可得到焦点区域的光强分布:

I(r,ψ,z)=TrW(r,ψ,z)==Wxx(r,ψ,z)+

Wyy(r,ψ,z)+Wzz(r,ψ,z)

(7)

下面通过(3)式～(6)式，进一步分析J0相关径向偏振涡旋光束在光轴和焦平面上的光谱相干度。

2 光轴上一对点的光谱相干度

现在考虑一对点在光轴上，即r1=(0,0,z1),r2=(0,0,z2)。通过使用以下公式:

πin-m[exp[i(n-m+1)ψ]Jn-m+1(krsinθ)+

exp[i(n-m-1)ψ]Jn-m-1(krsinθ)]

(8)

πin-m+1[exp[i(n-m+1)ψ]Jn-m+1(krsinθ)-

exp[i(n-m-1)ψ]Jn-m-1(krsinθ)]

(9)

由(3)式可得到交叉谱密度的表达式:

(10)

其中,

Jn+1(dNA-1δsinθ)exp(ikzcosθ)dθ

(11)

Jn-1(dNA-1δsinθ)exp(ikzcosθ)dθ

(12)

将z1=z2=z代入(7)式和(10)式，可得z轴上的总光强表达式为:

(13)

从上式中可以注意到，对于空间完全相干的涡旋光束(即δ=0)，只有拓扑荷数n=1时在光轴上有光强，也就是说，拓扑荷数n≠1时,在光轴上是零光强。然而，对于部分相干涡旋光束，则不存在这种零光强。通过分析(11)式和(12)式中的贝塞尔函数Jn+1和Jn-1，就很容易理解这些现象。而从(10)式、(13)式也很容易领悟一些对称性质，即:

Wij(0,0,-z1;0,0,-z2)=

Wij(0,0,z1;0,0,z2)*,(i,j=x,y,z)

(14)

I(0,0,-z)=I(0,0,z)

(15)

Fig.1 Real parts and modulus of μ(0,0,0;0,0,z)for different values of topological charge n and coherence parameter δ

应用(14)式 、(15)式和光谱相干度定义(6)式，可以获得以下对称关系:

μ(0,0,z1;0,0,z2)=μ(0,0,-z1;0,0,-z2)*

(16)

因此，应用(10)式、(13)式和(6)式，就可以确定z轴上的光强分布和纵向相干度。从(10)式可以看出，如果z1保持不变，μ(0,0,z1;0,0,z2)是z2的振荡函数，而振荡周期会随拓扑荷数n和相干参量δ变化。例如由图1所示，给出纵向相干度μ(0,0,0;0,0,z)实数部分和模随拓扑荷数n和相干参量δ的变化曲线。可以看到，在图1a～图1b中的μ(0,0,0;0,0,z)实数部分，其曲线的振荡周期略大于波长，约等于1.2λ～1.35λ。随拓扑荷数n的增加振荡周期会相应增大，但随相干参量δ的增加振荡周期只略有变化。从图1c～图1d中可知，在焦点区域的纵向相干度μ(0,0,0;0,0,z)的模随拓扑荷数n的增加或随相干参量δ的减小而增大。图1中其余参量σ=2,dNA=0.8。

图2为不同拓扑荷数n的纵向相干度|μ(0,0,z1;0,0,z2)|分布图,以及z1(λ)与z2(λ)的关系图。从图中可以清楚地看出，纵向相干度呈对称分布，其最小值出现在一对轴对称点z1，z2上，并且这个最小值点将随着拓扑荷数的增加而逐渐离开焦点。其余参量同图1。

Fig.2 Modulus of μ(0,0,z1;0,0,z2)for different values of topological charge n and coherence parameter δ

图3中给出了|μ(0,0,-z;0,0,z)|随拓扑荷数n、相干参量δ和数值孔径dNA变化的曲线图。从图中可见，在焦点区域的纵向相干度|μ(0,0,-z;0,0,z)|会随拓扑荷数n的增加而增大，但是也会随δ或dNA的减小而增加。特别是，当n≥5 (见图3a)或dNA≤0.4(见图3c)时，纵向相干度在焦点附近约等于1。其余参量同图1。

Fig.3 Modulus of μ(0,0,-z;0,0,z) for different values of topological charge n and coherence parameter δ

总之，在焦点区域纵向相干度随拓扑荷数n的增大而增大，也会随相干参量δ或数值孔径dNA的减小而增大,也就是说拓扑荷数越大，或入射光束的相干性越好(即δ越小)，或数值孔径越小(即dNA值越小)，都会使焦点附近光束的相干性越好。类似的研究见参考文献[12]和参考文献[17]。

3 焦平面上一对点的光谱相干度

接下来研究位于焦平面上的一对点的光谱相干度。其中，一个点位于几何焦点上,而另一个点位于焦平面上并与z轴有一定的距离。由于聚焦场的旋转对称性，可以假设这另一个点位于x轴上，研究相干度也不会失去一般性,即r1=(0,0,0),r2=(x,0,0)。将(4)式和(5)式代入(3)式，可得在焦平面上聚焦场的交叉谱密度:

(17)

其中,

Ji(dNA-1δsinθ1)Ji(dNA-1δsinθ2)×

J0(kxsinθ2)dθ1dθ2

(18)

Ji(dNA-1δsinθ1)Ji(dNA-1δsinθ2)×

J2(kxsinθ2)dθ1dθ2

(19)

同时，可获得焦平面上的总光强表达式为:

Qm-1,x(x)]2+[Qm+1,y(x)+Qm-1,y(x)]2}

(20)

其中,

Jm(dNA-1δsinθ)Jn-m+1(kxsinθ)dθ

(21)

Jm(dNA-1δsinθ)Jn-m-1(kxsinθ)dθ

(22)

因此，应用(17)式、(20)式和(6)式，就可以确定焦平面的光强分布和横向相干度。从(17)式中也可看出，交叉谱密度和横向相干度都是实数；而从(10)式中也可看出，纵向相干度是复数。

在图3中，给出了聚焦场的光强和横向相干度μ(0,0,0;x,0,0)随拓扑荷数n、相干参量δ和数值孔径dNA变化的曲线图。从图4a可以清楚地看出，当拓扑荷数n=1时，焦平面上有一小光斑；而当拓扑荷数n>1时，焦平面上有一小暗核，可作为光镊囚禁粒子。从图4b～图4f可以观察到，横向相干度会随着n,δ或者dNA的增大而减小。也就是说拓扑荷数越大，或入射光束的相干性越差(即δ越大)，或数值孔径越大(即dNA值越大)，都会使焦平面上光束的相干性越差，横向相干度衰减越快。此外，当n≥2或者δ≥3时(如图4b、图4d～图4f所示)，对于图中描绘的曲线与线μ(0,0,0;x,0,0)=0相交情况，横向相干度至少有两个零点，这种点也被称为光谱相干度的相位奇点或相干奇点。在这些点上，J0相关径向偏振涡旋聚焦场与点r1=(0,0,0)是完全不相干的。这说明部分相干矢量场也拥有相干奇点[12-17]。图4中其余参量同图1。

Fig.4 a—intensity distribution b～f—the degree of coherence μ(0,0,0;x,0,0) in the focal plane

4 结 论

基于矢量德拜理论和相干理论，研究了J0相关径向偏振涡旋光束通过大数值孔径透镜聚焦后，在光轴和焦平面上的空间相干特性。推导出在光轴和焦平面上的相干度表达式。数值模拟结果表明，在焦点区域纵向相干度随拓扑荷数n的增大而增大，随相干参量δ或数值孔径dNA的增大而减小。特别是，当n≥5或者dNA≤0.4时，纵向相干度在焦点附近约等于1。研究还发现，焦平面的横向相干度也会随拓扑荷数n、相干参量δ或者数值孔径dNA的增大而减小。此外，焦平面的横向相干度还具有相位奇点，这意味着在这种点上J0相关矢量涡旋聚焦场与几何焦点是完全不相干的。这些研究结果对J0相关矢量涡旋光束在许多领域的应用具有重要意义。