Compact focusing properties of radial vector beam with vortex phase encoding

矢量光束是偏振态在空间呈非均匀分布的光束。矢量光束在高数值孔径下的紧聚焦特性吸引了很多人的研究^[1-6]，例如:径向偏振光在紧聚焦的情况下会在沿轴方向产生纵向分量；而角向偏振光在紧聚焦的情况下会在焦平面产生中空的光场分布。由于矢量光束聚焦场的这些新颖特性，使其在粒子加速^[7-8]与捕获^[9-10]、扫描光学显微镜^[11-12]、激光切割与加工^[13-14]等领域应用广泛。

当光场相位围绕着中心点螺旋增加或减小时，称这种相位分布为涡旋相位。LI等人^[15]利用矢量衍射理论研究了涡旋相位的聚焦特性。研究表明，通过调整拓扑荷m^[16]以及聚焦系统，能够在聚焦场实现从单个焦点到两个重叠焦点的改变。这些特性适用于诸如光学捕获^[17-18]和粒子加速^[19]方面的相关应用。而最近的研究表明，当在矢量光束上加载涡旋相位后，聚焦场的强度分布将发生变化，展现出更加与众不同的特性。例如通过在角向偏振光场加载涡旋相位，HAO等人^[20]得到了比同等条件下的径向偏振光更小的焦点。

于是作者对径向矢量高斯光束加载涡旋相位的聚焦特性进行研究。基于矢量衍射理论，推导了径向矢量高斯光束编码单个涡旋的聚焦场表达式，并研究了聚焦场分布随拓扑荷变化的规律。在此基础上，进一步推导了径向矢量高斯光束编码对称的两个涡旋相位的聚焦场表达式，并给出了数值模拟结果。研究结果所展示的聚焦场特性对于矢量光束的调控具有重要的理论意义和实际价值。

矢量光束在高数值孔径物镜下的聚焦场由RICHARDS和WOLF提出的矢量Debye积分理论^[21]所描述。根据这个理论，对于径向偏振光，焦点附近的电场分布可以表示为：

式中，θ是物镜对焦点的最大孔径角, 也即物镜数值孔径; ϕ是极坐标下的角度; A(θ, ϕ)是出瞳切趾函数; λ是入射光波长; k是波数; α=arcsinθ是最大收敛角，由θ所决定; p, φ和z表示焦点处的柱矢量坐标; ${\mathit{\boldsymbol{\widehat e}}_\rho } $, ${{{\mathit{\boldsymbol{\widehat e}}}_\varphi }} $和${{{\mathit{\boldsymbol{\widehat e}}}_z}} $是柱坐标中的单位向量。二重积分的上下限一个为0~2π，它表示沿着角向积分一周；另一个为0~最大收敛角α，它表示沿着径向的积分。这两个积分范围共同确定对焦场有贡献的入射场的范围。

当径向矢量高斯光束编码单个涡旋相位时，其电场可以表示为：

式中，r表示极坐标下角向距离, m表示涡旋相位的拓扑荷, w为光束尺寸的常数。

由于物镜的设计通常满足正弦条件，于是有r=fsinθ^[22], f是高数值孔径物镜的焦距。将(3)式代入(1)式并化简，得到聚焦场${\mathit{\boldsymbol{\widehat e}}_\rho } $, ${{{\mathit{\boldsymbol{\widehat e}}}_\varphi }} $和${{{\mathit{\boldsymbol{\widehat e}}}_z}} $的电场分量为:

式中，$P\left( \theta \right) = \left\{ {\exp \left[ { - \left( {{\rm{i}}\mathit{k}z\cos \theta + {f^2}{{\left( {\sin \theta } \right)}^2}/{w^2}} \right)} \right]} \right\} \times \sqrt {\cos \theta } \sin \theta , $ J_m(x)是m阶第1类贝塞尔函数。该一重积分的范围就是径向需要积分的范围，即从0~α。此一重积分表示入射场的一个环状区域对聚焦场的贡献。

根据(4)式~(6)式, 取参量λ=632.8nm, w=2mm, θ=0.9，模拟计算了矢量光束编码单个涡旋相位的聚焦场。图 1中给出当拓扑荷m取0, 1, 2和3时所对应的聚焦场强度分布, 它展示了拓扑荷m对于焦场光斑图样的影响。由于入射光是径向偏振光，因此当拓扑荷m=0时，焦场只有一个焦点，如图 1a所示。随着拓扑荷数增加，在聚焦场的中心出现光学暗点，即2维中空聚焦场，并且中空聚焦场的半径也逐渐增大，如图 1c和图 1d所示。为了便于对比分析，将焦平面上x=0、沿y方向的强度分布显示图 2中。它更加直观地显示了聚焦场能量分布随拓扑荷m变化的数量关系。图 2说明，随着拓扑荷m的增加，聚焦场中心的光强逐渐变弱；当m=2和m=3时，能够产生2维中空聚焦场，并且中空场的半径也随着m的增加而增加。与角向矢量光的紧聚焦场相比，径向矢量光束编码单个涡旋时的紧聚焦场同时具有径向分量和纵向分量，这进一步提升了光学微操控^[23-24]的灵活性。

Figure 1.  Compact focusing field of radial vector Gaussian beam with single vortex phase encoding

Figure 2.  Intensity distribution of radial polarization component in x=0 plane along y direction

在径向矢量光束上编码两个对称的涡旋相位，并且两个涡旋相位的拓扑荷m相同。令拓扑荷m均为1的两个涡旋相位分别位于入射平面x=a和x=-a处。这样，入射光场的电场表达为：

两个涡旋相位之间的距离由a确定。当a=0时，两个涡旋相重叠，这与编码单涡旋相位拓扑荷m=2的情况相同。随着a逐渐增大，两个涡旋相位逐渐分离，聚焦场的分布也随着a的改变而出现新的变化。

将(7)式代入(1)式中并化简，聚焦场${\mathit{\boldsymbol{\widehat e}}_\rho } $, ${{{\mathit{\boldsymbol{\widehat e}}}_\varphi }} $和${{{\mathit{\boldsymbol{\widehat e}}}_z}} $的电场分量分别表示为：

式中, J₀(x), J₁(x), J₂(x)和J₃(x)分别为0阶、1阶、2阶以及3阶第1类贝塞尔函数。

根据(8)式~(10)式的推导结果，取参量λ=632.8nm, θ=0.9, 计算在两个涡旋之间的间距，即不同a值情况下焦场的强度分布，结果如图 3所示。图 3a为a=0w的情况，这等效于拓扑荷m为2的径向矢量光束编码单涡旋相位的聚焦场，它产生了一个2维中空场。随着两个涡旋相位之间距离的增加，聚焦场中开始出现两个沿y方向对称分布的光学暗点，如图 3b所示。随着涡旋之间距离的继续增加，这两个光学暗点之间的距离也随着增加，但其能量轮廓则不断下降，如图 3c所示。当a增加到1.0w时，两个光学暗点几乎消失，焦平面出现了一个亮的条形光班。总之，当两个涡旋重叠时，只有一个光学暗点出现；当两个涡旋相位逐渐分离开时，光学暗点开始分成两个并且随着涡旋相位之间距离的持续增加而消失。此外，对于特定的a值，两个光学暗点的能量分布和形状具有相同的特征。

Figure 3.  Compact focusing field of radial vector Gaussian beam with double vortex phase encoding

研究并给出了聚焦场的特性与涡旋相位的拓扑荷和间距之间的关系。发现，当编码单个涡旋相位并且拓扑荷为2和3时，聚焦场中存在一个2维中空能量分布，这种聚焦场适用于对单个粒子的微操纵，并且这种2维中空场能量分布的半径与拓扑荷m的大小相关，即通过调节拓扑荷的大小可调控中空场的大小，这增加了操控微粒大小的自由度。而当编码两个涡旋相位(拓扑荷m=1)时，聚焦场的分布依赖于两个涡旋相位之间的间距。当两个涡旋相位同时位于入射光束的中心时，聚焦场只会产生一个光学暗点，而随着两个涡旋相位中心的逐渐远离，此前的一个光学暗点逐渐分裂成两个暗点，这意味着这种聚焦场可以同时操控两个折射率不同于周围煤质的粒子。同时也发现，两个暗点的间距与涡旋相位的间距有关，因此可以通过调节涡旋相位间距来调整暗点的间距，从而增加对不同间距的双微粒的操控的自由度。