激光与光电子学进展, 2019, 56 (7): 070501, 网络出版: 2019-07-30  

用于产生空心焦场的二元环带相位板 下载: 927次

Binary Annular Phase Plate Used to Generate Hollow Focal Field
作者单位
安徽工业大学数理科学与工程学院, 安徽 马鞍山 243032
摘要
基于标量衍射理论,利用相干叠加的方法,研究了可以产生空心焦场的二元环带相位板。推导了各环带半径应满足的解析表达式,通过数值计算研究了环带半径对空心焦场的影响规律。研究结果表明,在焦平面上,相位板环带半径对旁瓣比和最大光强比的影响较大,对空心焦场半径的影响较小,增加环带数能够减小旁瓣比,增大最大光强比。在轴向上,环带数增加并不能增大最大光强比,环带数和环带半径对旁瓣比和最大光强比的影响较大。对于二元环带相位结构,旁瓣比和最大光强比仅取决于环带半径,和数值孔径无关;在焦面上,空心焦场半径和数值孔径成反比,在轴向上,空心焦场半径和数值孔径的平方成反比。
Abstract
Based on the scalar diffraction theory, a binary annular phase plate is studied by means of the coherent superposition method, which is used for the generation of a hollow focal field. The analytical expression of the radius for each ring is derived. The influence law of the ring radius on the hollow focal field is also investigated by the numerical calculation method. The research results show that, in the focal plane, the ring radius of this phase plate has a great influence on the side lobe ratio and the maximum light intensity ratio, but has a little influence on the radius of the hollow focal field. The increase in the number of rings reduces the side lobe radio, but increases the maximum light intensity ratio. In contrast, in the axial direction, the increase in the number of rings cannot increase the maximum light intensity ratio. Both the number of rings and the ring radius have relatively great influence on the side lobe ratio and the maximum light intensity ratio. For a binary annular phase structure, the side lobe ratio and the maximum light intensity ratio only depend on the ring radius and are independent of numerical aperture. In the focal plane, the radius of the hollow focal field is inversely proportional to the numerical aperture. In the axial direction, the radius of the hollow focal field is inversely proportional to the square of numerical aperture.

1 引言

空心焦场是指焦点处光强度为零,四周围绕着高强度光场的光学围笼,当轴向一段范围内中心都为暗场时可称为空心光束。空心焦场在光学捕获、光学成像和光学伪装等领域具有重要应用价值,如在光学操纵中可用于捕获低折射率的微粒[1-2],在超分辨荧光显微镜中可用作擦除光[3-4]

近年来,随着激光技术的发展,空心焦场的产生和应用受到广泛关注,刘双龙等[5]分析了相干拉曼散射显微镜中空心焦场的产生和应用,Peng等[6]利用方位偏振光产生的空心焦场捕获了低折射率的微粒,Xue等[7]研究了方位偏振光聚焦系统产生的空心焦场在荧光显微镜中的应用。目前,产生空心焦场的方法主要有轴棱镜法[8]、光纤法[9]、特殊偏振光合成法[10]和衍射光学元件法等[11-12]。其中,衍射光学元件法具有光路简单,体积小和重量轻的特点,日益受到学者们的重视[13-14]

本文利用标量衍射理论,研究了通过二元环带相位板产生空心焦场的方法,利用相干叠加的方法得到环带相位板的各环带半径的解析关系式,然后利用梯度迭代算法搜索得到不同环带数目的最优值,分别分析了轴向和焦平面上空心焦场的环带半径、最大光强比和旁瓣比的变化规律,利用数值拟合得到了空心焦场半径与物镜数值孔径之间的关系。

2 理论分析

根据标量衍射理论,物镜焦点附近光场的计算式为[15]

U(v,u)=201PρJ0()exp(iuρ2/2)ρdρ,(1)

式中P(ρ)为光瞳函数的复振幅透过率,ρ为环带半径,J0(·)代表零阶贝塞尔函数,vu分别为径向和轴向光学坐标,表达式分别为

v=kNAr,(2)u=k(NA)2z,(3)

式中k=2π/λ,NA为数值孔径,λ为光波波长,rz分别为径向和轴向坐标,如图1所示,图中F点焦点,P点为焦点附近的考察点。

图 1. 焦场计算示意图

Fig. 1. Schematic for calculation of focal field

下载图片 查看所有图片

二元环带相位板广泛应用于光学超分辨和光束整形等领域[16-17],二元环带相位板的光瞳函数可以表示为

P(ρ)=1,0<ρ<ρ1exp(),ρ1ρ<ρ21,ρ2ρ<ρ3    exp(),ρj-1ρ<ρj,(4)

式中ρj为第j个环带半径,j为正整数,φ为调制相位,通常为了不产生焦移,取φ=π。为了产生空心光束,将焦点处光强置为零,由(1)式得:

U(0,0)=201P(ρ)ρdρ(5)

假定单位振幅平面波入射,则P(ρ)=exp(iφ),(φ=0或π),那么第j个环带孔径在焦点处复振幅可以表示为

Uj(0,0)=2exp()ρj-1ρjρ=exp()(ρj2-ρj-12)(6)

对于多个环带结构的光瞳滤波器,焦点处的复振幅可看作多个环带在焦点处复振幅的叠加,对于二区相位结构,设φ=π,焦点处复振幅为

U2(0,0)=ρ12+(ρ22-ρ12)exp()=ρ12-(ρ22-ρ12),(7)

对于空心焦场,U(0,0)=0,将最大半径归一化,即ρ2=1,可求出

ρ1=2/2,(8)

这与文献[ 18]报道的二环结果一致。对于三环及以上的π相位结构,同样可求出产生空心焦场时各环带半径应满足的关系,如表1所示。

表 1. 产生空心焦场时各环带半径

Table 1. Radius of each ring when generating hollow focal field

Number of ringsRelation among normalized ring radii
22ρ12-1=0
32(ρ12-ρ22)+1=0
42(ρ12-ρ22+ρ32)-1=0
52(ρ12-ρ22+ρ32-ρ42)+1=0
62(ρ12-ρ22+ρ32-ρ42+ρ52)-1=0

查看所有表

表1结果可以归纳为

U(0,0)=2(ρ12-ρ22+ρ32-ρ42++ρj-22-ρj-12)+1,jisodd2(ρ12-ρ22+ρ32--ρj-22+ρj-12)-1,jiseven(9)

因此,要设计产生空心焦场的任意环带π相位板,可从(9)式中求出相应的环带半径,(9)式只满足焦点处光强为零,对于空心焦场的半径、光强分布和旁瓣大小等光学特性没有给出具体关系,还需求解(1)式,通过优化搜索算法寻找合适的结构。通过数值解法分析环带半径变化时空心焦场的变化规律可以利用梯度迭代算法[19],将最大半径归一化,先设定第环带半径变化范围,相邻两个环带半径间隔设置为0.01,最靠近外环的半径利用(9)式求出,利用(1)式求解光强分布,然后求出每个结构对应的最大光强、最大旁瓣强度和半径等特征值,通过设定特征值的阈值条件,得到与之相对应的环带结构。

3 结果分析与讨论

为了表征空心焦场的光强分布特性,利用参考文献[ 17]中的超分辨位相板的定义方法。引入空心焦场半径R、最大光强比S和旁瓣比M作为评价参数。R定义为焦点到最大光强点的距离;S定义为空心焦场最大光强和爱里斑中心最大光强的比值;M定义为旁瓣最大光强和空心焦场最大光强的比值。计算时,数值孔径取0.2,利用梯度迭代算法分别计算焦平面与轴向的光强分布,得到不同环带结构下RSM,分析环带半径对三者影响。

3.1 焦面光强分析

对于二环结构,由于其环带半径ρ1=0.707唯一,由(1)式可求出R=2.864λ,M=0.286,S=0.108,其在焦平面上的一维光强分布如图2所示。

图 2. 二环焦面光强分布

Fig. 2. Light intensity distribution in focal plane for two-ring structure

下载图片 查看所有图片

采用三环及以上的结构,由于各个环带之间的干涉作用,能量会重新分配,RMS随着环带半径的改变而变化。实际应用时,对于单环衍射场,希望M较小,S较大,对三环结构进行迭代计算,设定M<0.3可以得到满足条件的三环结构(对RS不做限制)。分析数据发现,三环结构可以产生两种半径的空心光束,一种R为5.4λ~5.7λ,另一种R为2.8λ~3.0λ。当R为5.4λ~5.7λ时,S变化较小,其值为0.61~0.65,M存在一个极小值;同样,R为2.8λ~3.0λ,S变化较小,但数值较大,M也存在一个极小值,相应的结果如表2所示,其在焦平面上的一维光强分布如图3所示。

表 2. 优化后的三环结构

Table 2. Optimized three-ring structure

ρ1ρ2R /λSM
0.350.7895.6490.0650.208
0.690.9882.920.0980.288
0.010.7072.8450.1080.287

查看所有表

图 3. 三环焦面光强分布。(a) R=5.65λ;(b) R=2.85λ

Fig. 3. Light intensity distributions in focal plane for three-ring structure. (a) R=5.65λ; (b) R=2.85λ

下载图片 查看所有图片

和二环结构进行比较可以发现,图3(b)和图2非常相似,而且RMS也差别不大,分析原因可发现,此时三环对应的结构和二环结构非常接近,如表2优化后的三环结构中ρ1=0.01,ρ2=0.707,由于第一环带非常小,对光场的调制作用很弱,因此结果和二环结构相近。图3(a)中R较大,此时最大值出现在干涉第二主极大处,而在图3(b)和图2中,最大值出现在干涉第一主极大处。

分析4环结构可以发现,增加环带数可以减小M。搜索时设定M<0.2,此时,环带半径对R影响较小,R为2.8λ~3.1λ,与2环和3环在干涉第一主极大处形成的空心光束半径差别不大,但对MS影响较大。当M<0.2时,四环结构S为0.04~0.09,S增大会引起M相应增大,存在几个最优结构,S较大时,对应M的极小值,具体数值如表3所示。和2环和3环的优化结构进行比较,可以发现,在半径R差别不大的情况下,4环结构能够显著的降低M,但同时S也有所减小。

表 3. 优化后的四环结构

Table 3. Optimized four-ring structure

ρ1ρ2ρ3R /λSM
0.60.730.822.880.070.12
0.580.690.802.840.070.12
0.570.700.8162.880.0650.12
0.550.650.7872.840.070.11

查看所有表

计算5环结构时,设定M<0.2和S>0.05。通过分析发现,5环结构和3环结构类似,能够产生两种半径的空心焦场,一种R为5.4λ~5.6λ,另一种为2.8λ~3.1λ,其最大值分别位于干涉第二主极大和第一主极大处。当R为5.4λ~5.6λ时,S为0.05~0.06,M为0.15~0.20;而R为2.8λ~3.1λ时,S为0.05~0.08,M为0.11~0.20。同样,增大S会使得M增大,存在最优的结构,可在S较大时,使得M极小,如表4所示。可以看出,5环和3环结构都可以产生两种R的空心焦场,5环结构可以得到更小的M,同时S略有降低。5环结构和4环结构比较,在相同的R下,M相差不大,而4环结构的S更大。

表 4. 优化后的5环结构

Table 4. Optimized five-ring structure

ρ1ρ2ρ3ρ4R /λSM
0.360.680.720.8285.520.050.15
0.540.640.760.983.010.0600.11
0.550.660.770.983.010.0590.10
0.560.670.780.982.910.0640.11

查看所有表

对于6环结构,在与5环结构相同的限制条件下进行优化搜索。分析结果发现,6环结构和2环、4环结构类似,只能产生一种半径的空心焦场,其R为2.7λ~3.1λ,但是6环结构对应的解更多。和4环结构相比较,6环结构能够略微增大S,但M变化不大,R也变化不大。同样,S增大会使得M增大,存在最优的结构,可在S较大时,使得M极小,如表5所示。

表 5. 优化后的6环结构

Table 5. Optimized six-ring structure

ρ1ρ2ρ3ρ4ρ5R /λSM
0.400.410.600.730.832.910.0660.10
0.390.400.580.690.802.810.0660.10
0.490.500.600.710.812.910.0710.11
0.500.510.610.730.822.910.0710.11

查看所有表

分析2环到6环的计算结果可以发现,环带数增加能够显著减小M,但并没有使得S明显增强。在M<0.3时,空心焦场环带半径和相位板环带数目的奇偶性有关。偶数环带只能产生一种半径的空心焦场,R为2.8λ~3.1λ,最大光强在干涉第一主极大处;奇数环带能够产生两种半径的空心焦场,R分别为5.4λ~5.6λ和2.8λ~3.1λ,其最大光强分别位于干涉第二主极大和第一主极大处。不同环带结构焦平面上的一维光强分布如图4所示。从图4可以看出,随着环带的增加,旁瓣的数目增加,M降低是由于第一旁瓣的能量转移到其他旁瓣,并没有转移到主瓣,进而导致S没有增大。

图 4. 不同环带结构的焦面光强分布。(a) 2环;(b) 3环;(c) 4环;(d) 5环;(e) 6环

Fig. 4. Light intensity distributions in focal plane for different ring structures. (a) 2-ring; (b) 3-ring; (c) 4-ring; (d) 5-ring; (e) 6-ring

下载图片 查看所有图片

为了确定物镜数值孔径对空心焦场特性的影响,分别选取2环到6环的优化结构,研究NA从0.1增大到0.6时RSM的变化规律,以4环和5环结构为例,三者随NA的变化曲线如图5所示。

图 5. 焦面上R、M和S随NA的变化。(a) 4环R值;(b) 5环R值;(c) 4环S、M值;(d) 5环S、M值

Fig. 5. R, M, S in focal plane versus NA. (a) R value for 4-ring; (b) R value for 5-ring; (c) S and M values for 4-ring; (d) S and M values for 5-ring

下载图片 查看所有图片

图5中对应的4环结构参数为ρ1=0.60,ρ2=0.73,ρ3=0.82,5环结构参数为ρ1=0.57,ρ2=0.65,ρ3=0.76,ρ4=0.99。可以看出,SM值不随NA的改变而变化,R值随NA的增加而减小。通过拟合R曲线,可以得到RNA的关系:

R=/(NA)(10)

图5(a)和图5(b)中m分别为0.5768和0.5856。分析其他环带结构也能得出相同结论,即RNA成反比,m和具体的结构参数有关,而SM值不随NA变化,均为常数。

3.2 轴向光强分析

入射波前的调制不仅可以改变焦面光强分布,也会影响子午面内的光强分布,可以形成多焦点、长焦深等光学现象。(1)式中,令ν=0,则可利用数值积分求出轴向上的光强分布。计算时直接利用优化焦面结构得到的环带相位板结构数据,计算轴向的光强分布。由于希望得到和焦面相似的光强分布,即中心光场为零,四周围绕高强度的光场,因此,和焦面类似,同样采用RSM值表征轴向上空心焦场的光学特性。通过设定RSM的阈值,可以搜索得到满足要求的结构数据。图6给出了2~6环不同结构调制下轴向的光强分布,其对应的相位板结构参数如表6所示,对应轴向上M取最小值的情况。从图6可以看出,改变各个环带的面积,由于环带之间的干涉作用,在焦点处光强为0,传播一段距离后,轴向光强达到峰值,具有较小的M[18]

图 6. 不同环带对应的轴向光强分布。(a) 2环;(b) 3环;(c) 4环;(d) 5环;(e) 6环

Fig. 6. Light intensity distributions in axial direction for different ring structures. (a) 2-ring; (b) 3-ring; (c) 4-ring; (d) 5-ring; (e) 6-ring

下载图片 查看所有图片

表 6. 图6中对应的环带半径和R、S、M值

Table 6. Ring radii and R, S, M values corresponding to Fig. 6

Number of ringsρ1ρ2ρ3ρ4ρ5R /λSM
20.70737.0460.5250.088
30.010.70737.0360.5250.088
40.590.650.75835.0340.3920.153
50.560.620.730.98138.0370.3320.174
60.50.540.690.730.77435.2840.3760.079

查看所有表

分析计算结果可知,环带数目对R值影响较小,偶数环带相位板只能产生一种R值,约为35λ,奇数环带相位板能产生两种R值,约为37λ~51λ

相位板的环带数目和环带半径对SM的影响较大,M值增大伴随S值减小,环带数增加并不能增大S值,但轴向S值要远大于焦面S值。焦面上M取极值时对应的环带结构和轴向取极值时对应的环带结构并不一致,即当焦面上M取极小值时,轴向上M并非极小值。如表5所示的6环结构中,当ρ1=0.49,ρ2=0.50,ρ3=0.60,ρ4=0.71,ρ5=0.81时,轴向的R=32.53λ,S=0.2714,M=0.6131,可见焦面M取极小值时,轴向M并非极小值。研究发现,当焦面M限定一个极大值时,在轴向上M存在一个极小值,对应S取得最大值,如对于4环结构,当焦面M<0.2时,在轴向上,S最大值为0.392,M最小值为0.153。轴向上,二环结构的S取极大值时,M取得较小值,增加环带数目,S不会增大,反而会减小,同时M会增大,虽然6环结构的M比2环结构略小,但其S要比2环结构小很多。

为了确定数值孔径对轴向光束的影响,以4环和5环结构为例,结合表6研究NA从0.1增大到0.6时RSM的变化规律,其变化规律如图7所示。

图 7. 轴向上R、M和S随NA的变化。(a) 4环R值;(b) 5环R值;(c) 4环S、M值;(d) 5环S、M值

Fig. 7. R, M, S in axial direction versus NA. (a) R value for 4-ring; (b) R value for 5-ring; (c) S and M values for 4-ring; (d) S and M values for 5-ring

下载图片 查看所有图片

可以看出,RMS值的变化规律和焦面上的变化规律相似,R值随NA的增加而减小,而SM值不随NA变化,均为常量。通过拟合R曲线,可以得到RNA的关系:

R=/(NA)2,(11)

对应图7(a)和(b)中m分别为1.393和1.515。

4 结论

利用环带二元相位板对入射光场相位进行调制,可在焦点区域形成多种新颖光束,经过理论推导得到了焦点光场与相位板环带半径的解析表达式,分析了相位板环带半径对空心焦场RSM值的影响规律。研究表明,增加相位板环带数目,并不能提高轴向上和焦平面上的S值,但可以降低焦面上的M值,空心焦场的R和相位板环带数目的奇偶性有关,并受相位板环带半径的影响。对于环带二元相位板结构,焦面上的R和数值孔径成反比,轴向上的R和数值孔径的平方成反比,而SM和数值孔径无关。

参考文献

[1] Xu H F, Zhang W J, Qu J, et al. Optical trapping Rayleigh dielectric particles with focused partially coherent dark hollow beams[J]. Journal of Modern Optics, 2015, 62(21): 1839-1848.

[2] Kuga T, Torii Y, Shiokawa N, et al. Novel optical trap of atoms with a doughnut beam[J]. Physical Review Letters, 1997, 78(25): 4713-4716.

[3] Liu H L, Lü Y F, Xia J, et al. Flat-topped vortex hollow beam and its propagation properties[J]. Journal of Optics, 2015, 17(7): 075606.

[4] Khonina S N, Golub I. Engineering the smallest 3D symmetrical bright and dark focal spots[J]. Journal of the Optical Society of America A, 2013, 30(10): 2029-2033.

[5] 刘双龙, 刘伟, 陈丹妮, 等. 超衍射极限相干反斯托克斯拉曼散射显微成像技术中空心光束的形成[J]. 物理学报, 2014, 63(21): 214601.

    Liu S L, Liu W, Chen D N, et al. Generation of dark hollow beams used in sub-diffraction-limit imaging in coherent anti-Stokes Raman scattering microscopy[J]. Acta Physica Sinica, 2014, 63(21): 214601.

[6] Peng F, Yao B L, Yan S H, et al. Trapping of low-refractive-index particles with azimuthally polarized beam[J]. Journal of the Optical Society of America B, 2009, 26(12): 2242-2247.

[7] Xue Y, Kuang C F, Li S, et al. Sharper fluorescent super-resolution spot generated by azimuthally polarized beam in STED microscopy[J]. Optics Express, 2012, 20(16): 17653-17666.

[8] Lin J H, Wei M D, Liang H H, et al. Generation of supercontinuum bottle beam using an axicon[J]. Optics Express, 2007, 15(6): 2940-2946.

[9] Chen Y H, Yan L, Rishøj L, et al. Dynamically tunable optical bottles from an optical fiber[J]. Optics Letters, 2012, 37(16): 3327-3329.

[10] Zeng T T, Ding J P. Three-dimensional multiple optical cages formed by focusing double-ring shaped radially and azimuthally polarized beams[J]. Chinese Optics Letters, 2018, 16(3): 031405.

[11] Wang H F, Gan F X. Phase-shifting apodizers for increasing focal depth[J]. Applied Optics, 2002, 41(25): 5263-5266.

[12] 陈国钧, 周巧巧, 纪宪明, 等. 用线偏振光产生可调矢量椭圆空心光束[J]. 物理学报, 2014, 63(8): 083701.

    Chen G J, Zhou Q Q, Ji X M, et al. Generation of the tunable vector ellipse hollow beam by using linearly polarized light beams[J]. Acta Physica Sinica, 2014, 63(8): 083701.

[13] Lin J, Chen R, Yu H C, et al. Generation of hollow beam with radially polarized vortex beam and complex amplitude filter[J]. Journal of the Optical Society of America A, 2014, 31(7): 1395-1400.

[14] 陈国钧, 周巧巧, 纪宪明, 等. π相位板产生矢量光束的高数值孔径聚焦特性研究[J]. 光学学报, 2014, 34(12): 1226001.

    Chen G J, Zhou Q Q, Ji X M, et al. Study on high-numerical-aperture-focused characteristics of vector beam produced by π phase plate[J]. Acta Optica Sinica, 2014, 34(12): 1226001.

[15] de Juana D M, Oti J E, Canales V F, et al. . Design of superresolving continuous phase filters[J]. Optics Letters, 2003, 28(8): 607-609.

[16] 蔡勋明, 赵晶云, 童红, 等. π相位光阑调制的高数值孔径聚焦特性研究[J]. 光学学报, 2017, 37(5): 0526001.

    Cai X M, Zhao J Y, Tong H, et al. Study on high-numerical-aperture-focused characteristics by π phase diaphragm modulation[J]. Acta Optica Sinica, 2017, 37(5): 0526001.

[17] Sales T R M, Morris G M. Diffractive superresolution elements[J]. Journal of the Optical Society of America A, 1997, 14(7): 1637-1646.

[18] Bokor N, Davidson N. Tight parabolic dark spot with high numerical aperture focusing with a circular π phase plate[J]. Optics Communications, 2007, 270(2): 145-150.

[19] Zhou C H, Jia J, Liu L R. Circular Dammann grating[J]. Optics Letters, 2003, 28(22): 2174-2176.

王伟, 吴建光, 莫绪涛, 黄仙山. 用于产生空心焦场的二元环带相位板[J]. 激光与光电子学进展, 2019, 56(7): 070501. Wei Wang, Jianguang Wu, Xutao Mo, Xianshan Huang. Binary Annular Phase Plate Used to Generate Hollow Focal Field[J]. Laser & Optoelectronics Progress, 2019, 56(7): 070501.

引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!