光学仪器  2021, Vol. 43 Issue (4): 71-77   PDF    
具有亚波长结构的突然自聚焦光束的频谱调控研究
施雨1,2, 耿滔1,2     
1. 上海理工大学 上海市现代光学系统重点实验室,上海 200093;
2. 上海理工大学 光电信息与计算计工程学院,上海 200093
摘要: 为了提高具有亚波长结构的圆对称艾里光束的自聚焦特性,将光束的工作波长扩大至太赫兹波段,并在频谱空间中分别加入高通和带通滤波进行调控。对两种滤波效果进行对比发现,通过滤波函数的参数变化可以控制空间频谱的分布,从而对光束的传播特性进行调制。当选取合适的参数时,焦斑尺寸减小,焦点强度增加,突然自聚焦特性会大幅增强,同时可实现焦点位置的控制。比较两种调制方法后发现,带通滤波后的光束拥有更好的自聚焦特性。
关键词: 太赫兹    突然自聚焦    高通滤波    带通滤波    傅里叶空间    
Spatial spectrum modulation for abruptly autofocusing beams with subwavelength structures
SHI Yu1,2, GENG Tao1,2     
1. Shanghai Key Laboratory of Modern Optical System, University of Shanghai for Science and Technology, Shanghai 200093, China;
2. School of Optical-Electrical and Computer Engineering, University of Shanghai for Science and Technology, Shanghai 200093, China
Abstract: In order to improve the autofocusing property of circular Airy beam with subwavelength structures, the working wavelength of the beam is expanded to terahertz band. This paper studies the propagation properties of the terahertz circular Airy beam by high-pass and band-pass filtering in Fourier space. It is found that the parameter changes of the high-pass and band-pass functions can control the distribution of the spatial frequency spectrum. Thus, the propagation characteristics of the beam are modulated. When the parameters are chosen appropriately, the size of focal spot is reduced, the maximum focal intensity and especially the abruptly autofocusing property could be greatly enhanced. Meanwhile, the focal position can be controlled. Compared with the high-pass filtering, the beam after the band-pass filtering has better abruptly autofocusing characteristics.
Key words: terahertz    abruptly autofocusing    high-pass filter    band-pass filter    Fourier space    
引 言

2007年,有限能量的艾里光束(Airy beam, AB)第一次被提出并实验实现[1-2],经过不断地研究,发现其拥有无衍射、自加速、自修复等特点[3-6]。随后,研究者将AB作了径向对称处理[7-8],产生了圆艾里光束(circular Airy beam, CAB)并发现CAB具有独特的自聚焦特性,具体表现是光束在焦点前能保持相当低的强度分布,但是在能量传播到焦点时光强会突然提升至数十甚至数百倍。CAB也因为这种独特的性质在生物医学、激光加工等领域有着巨大的潜在应用价值[9-15],而自聚焦焦斑的峰值强度、焦斑的大小和位置对这些应用有着至关重要的影响。

目前,已经有很多研究者在提升CAB自聚焦能力的方向做出了努力。Lin等[16]发现去除CAB初始面的主光环,能增强光束的自聚焦能力;Deng等[17]发现在CAB的频谱面采用高通滤波的方法可以更好地改善光束的自聚焦能力;最近,我们发现与高通滤波相比,在频谱面引入高斯包络调制,能够使得光束的自聚焦能力获得进一步提升[18]。可见光的光源虽然更易获得,但目前CAB的研究主要集中在傍轴近似区域,即自聚焦焦距远远大于光束截面半径。而自聚焦焦距与光波的波长成反比,即波长增大,焦距减小。因此,当波长增大到太赫兹频段时,可以有效地缩短自聚焦焦距,焦距与光束截面半径大小相当,光波的传播由傍轴区域迈向离轴区域,且此时CAB的初始光场分布具有与波长尺度相当甚至小于波长的精细结构,由于焦距的减小,光束的会聚角增大,可获得强度更大的焦点强度。对于可见光而言,波长远小于初始光斑的空间结构,因此其频谱都在传输波频率范围内。而对于太赫兹波而言,其波长与初始光斑的空间结构已处于同一数量级,特高频分量已进入倏逝波频段,无法传播到远场,因此与可见光只需要低频滤波不同,太赫兹波需要同时滤去低频和特高频分量才能有效地提升光束的自聚焦能力。He等[19]将工作波长从可见光移到了太赫兹波段,使得初始光束具有了接近波长尺寸甚至亚波长尺寸的光强分布空间结构,光束的自聚焦焦斑尺寸接近了波长数量级。

本文将上述两种方法相结合,将工作波长移到太赫兹波段,考虑到频谱分布的不对称性,在频谱面引入高通滤波器(high-pass filter, HPF),以减少光束的低频分量。同时,由于太赫兹CAB在初始面具有亚波长的空间结构,使得特高频分量进入到倏逝波范围,无法传输到远场,因此本文又引入低通滤波器(low-pass filter, LPF),滤去特高频分量。这样HPF和LPF共同组成了带通滤波器(band-pass filter, BPF)。以下将HPF调制的CAB简称为HCAB,BPF调制的CAB简称为BCAB。本文详细探讨了CAB经BPF滤波后的光束自聚焦特性,并与HPF的结果进行了对比,结果发现使用BPF确实可以增强CAB的聚焦特性,且比单独使用HPF的结果更好。在合适的参数下,BPF调制后的聚焦强度可提高至CAB的7倍,半高全宽也有所改善,且可实现焦距的控制,所以该方法是一种提高太赫兹CAB自聚焦性能的有效方法。本文的研究以期能对太赫兹成像和太赫兹光谱领域的研究提供有益的参考。

1 理 论

CAB在入射平面处的电场表达式为[17]

${u_0}\left( r \right) = {E_0}Ai\left( {\frac{{{r_0} - r}}{\omega }} \right)\exp \left[ {\alpha \left( {\frac{{{r_0} - r}}{\omega }} \right)} \right]$ (1)

式中: ${E_0}$ 为振幅常数; $Ai\left( * \right)$ 为艾里函数; ${r_0}$ 为主光环的半径; $\omega $ 为控制光斑宽度; $\alpha $ 为指数截断因子; $r$ 为圆坐标。

目前由于CAB在频谱面上还没有确切的解析式,因此研究通常使用Hankel变换方法获得数值解,表达式为

$g\left( k \right) = \int_0^\infty {{u_0}\left( r \right){J_0}\left( {kr} \right)r{\rm{d}}r} $ (2)

式中: $k$ 为径向空间频率; ${J_0}\left( * \right)$ 为零阶贝塞尔函数。再次通过Hankel变换,可以得到CAB的传输方程为

$u\left( {r,z} \right) = \int_0^\infty {g\left( k \right){J_0}\left( {kr} \right)\exp \left( {2i {\text{π}} z\sqrt {{\lambda ^{ - 2}} - {k^2}} } \right)k{\rm{d}}k} $ (3)

式中: $z$ 为沿光轴的传播距离; $\lambda $ 为光束的波长。为了便于比较聚焦能量,文中所有的计算都令光束的初始光强为1,即 $2 {\text{π}} {\int {\left| {{u_0}\left( r \right)} \right|} ^2}r{\rm{d}}r = $ $ 1$ 。CAB的参数为 ${r_0} = 5\;{\rm{mm}}$ $\omega = 0.3\;{\rm{mm}}$ $\lambda = $ $ 0.4\;{\rm{mm}}$ $\alpha = 0.1$ ,并保持不变。

为了减少CAB的低频分量需引入HPF,用 ${f_{\rm{H}}}$ 表示为

${f_{\rm{H}}} = \frac{1}{{1 + \exp \left[ { - {\beta _{\rm{H}}}\left( {k - {k_{\rm{H}}}} \right)} \right]}}$ (4)

式中: ${\;\beta _{\rm{H}}}$ 为控制增长速度; ${k_{\rm{H}}}$ 为截止位置[17]

另一方面,为了滤去进入倏逝波的特高频成分需引入LPF,用 ${f_{\rm{L}}}$ 表示为

${f_{\rm{L}}} = \frac{1}{{1 + \exp \left[ {{\beta _{\rm{L}}}\left( {k - {k_{\rm{L}}}} \right)} \right]}}$ (5)

式中: ${\;\beta _{\rm{L}}}$ 为控制衰减速度; ${k_{\rm{L}}}$ 为截止位置。

那么,由 ${f_{\rm{H}}}$ ${f_{\rm{L}}}$ 可组成BPF,用 ${f_{\rm{B}}}$ 表示为

${f_{\rm{B}}} = {f_{\rm{H}}}{f_{\rm{L}}}$ (6)

图1给出了 ${f_{\rm{H}}}$ ${f_{\rm{B}}}$ 的示意图,其中, ${\;\beta _{\rm{H}}} = 15$ ${k_{\rm{H}}} = 1.35$ ${\beta _{\rm{L}}} = 13$ ${k_{\rm{L}}} = 1.88$ 。经这两种滤波器滤波后,光束的频谱变为 ${g_{\rm{H}}}\left( k \right) = {C_{\rm{H}}}g\left( k \right){f_{\rm{H}}}\left( k \right)$ ${g_{\rm{B}}}\left( k \right) = {C_{\rm{B}}}g\left( k \right){f_{\rm{B}}}\left( k \right)$ ,其中 ${C_{\rm{H}}}$ ${C_{\rm{B}}}$ 为频谱的振幅常数,使滤波后的总能量保持不变,以满足能量守恒。图2给出了 $g\left( k \right)$ ${g_{\rm{H}}}\left( k \right)$ ${g_{\rm{B}}}\left( k \right)$ 的计算结果,反映CAB,HCAB和BCAB的能量归一化傅里叶频谱,模拟使用的参数同上,其中F表示复振幅的大小。从图2可看出,当CAB被 ${f_{\rm{H}}}$ 调制后,频谱的前峰逐渐消失(点线);当CAB的频谱被 ${f_{\rm{B}}}$ 修改后,频谱的前峰和后峰都在慢慢消失(短横线)。

图 1 HPF(实线)和BPH(横线)的函数示意图 Figure 1 The function of HPF(solid line)and BPH(dash line)

图 2 CAB(实线)、HCAB(点线)和BCAB(短横线)的能量归一化傅里叶频谱图 Figure 2 Energy normalized Fourier spectrum of CAB(solid line), HCAB(dotted line), and BCAB(short dash line)
2 模拟分析

根据式(3),令z=0可以得到初始平面处CAB、HCAB、BCAB的强度分布,其强度随初始平面上的径向距离 $r$ 变化而变化如图3所示。其中, ${k_{\rm{H}}} = 1.35$ ${k_{\rm{L}}} = 1.88$ 保持不变,图3(a) ${\;\beta _{\rm{H}}} = {\beta _{\rm{L}}} = 40$ 图3(b) ${\;\beta _{\rm{H}}} = 15$ ${\;\beta _{\rm{L}}} = 13$ 。显然,HCAB和BCAB的最大强度远小于CAB,HCAB强度最小。由于来自调制函数的衍射效应,有一段小的峰出现在HCAB和BCAB的第一个环之前。HCAB和BCAB最大强度的位置移动到后面的峰。HCAB和BCAB第8个峰成为最高峰。BCAB中有近5个峰,相较于HCAB和CAB有更少的峰值。(a)、(b)对比可以看出 ${\beta _{\rm{H}}}$ ${\beta _{\rm{B}}}$ 的不同会对初始平面有影响,当 ${\beta _{\rm{H}}}$ ${\beta _{\rm{L}}}$ 较大时,相比选取小的参数,随着半径越来越大,HCAB会出现不规则的锯齿段,这是由于 ${\beta _{\rm{H}}}$ ${\beta _{\rm{L}}}$ 的取值过大引起的。

图 3 初始平面处CAB(实线)、HCAB(点线)和BCAB(短横线)的强度分布 Figure 3 The intensity distributions of CAB(solid line), HCAB(dotted line), and BCAB(short dash line) on the initial plane

CAB、HCAB和BCAB在数值模拟中的传播强度曲线如图4所示,为了清楚地显示传播特性,对其做了对数处理。图中CAB最重要的自聚焦特性在HCAB和BCAB中仍然保持不变。但是通过比较发现,HCAB和BCAB的焦点位置都发生了变化,相较于CAB能在更近的位置聚焦。更为有趣的是,当在傅里叶空间中减小高频的部分空间频谱时(见图4(b)(c)),焦斑纵向长度也随之减小,焦斑长度与初始环数成正比。与(a)、(b)相比,(c)中光束汇聚到焦点的能量更多。

图 4 CAB,HCAB和BCAB在自由空间中的传播强度 Figure 4 Propagation intensity of CAB, HCAB and BCAB in free space

CAB焦点在z = 12 mm处,模拟计算可知 ${\beta _{\rm{H}}}$ ${\beta _{\rm{B}}}$ 的变化对焦点的位置没有太大影响,而随着 ${k_{\rm{H}}}$ ${k_{\rm{L}}}$ 的变化,HCAB和BCAB焦距的变化可以从图5表现。显然,随着 ${k_{\rm{H}}}$ ${k_{\rm{L}}}$ 的增大,HCAB和BCAB焦距不断减小。其中,BCAB的焦距始终小于CAB。

图 5 不同参数下调制得到的焦点位置 Figure 5 The focus position obtained by modulation under different parameters

在实际应用中,焦点处的峰值强度是一个极其重要的参数。通过 ${{{{\left( {{I_{{\rm{HCAB}}}}} \right)}_{\max }}} / {{{\left( {{I_{{\rm{CAB}}}}} \right)}_{\max }}}}$ ${{{{\left( {{I_{{\rm{BCAB}}}}} \right)}_{\max }}} / {{{\left( {{I_{{\rm{CAB}}}}} \right)}_{\max }}}}$ 反映调制后不同参数下光束的聚焦能力,其中 ${\left( {{I_{{\rm{CAB}}}}} \right)_{\max }}$ ${\left( {{I_{{\rm{HCAB}}}}} \right)_{\max }}$ ${\left( {{I_{{\rm{BCAB}}}}} \right)_{\max }}$ 分别表示CAB、HCAB和BCAB的焦斑最大光强。如图6所示,图6(a)中HCAB的焦点光强随着 ${k_{\rm{H}}}$ 的变大不断增强,同时 ${\;\beta _{\rm{H}}}$ 变大其强度也随之增大,但是不如 ${k_{\rm{H}}}$ 敏感。当 ${k_{\rm{H}}}$ 超过1.35时,强度开始减小,即HCAB想要增强聚焦强度,并不意味着一直增大 ${k_{\rm{H}}}$ 。(b)中BCAB的焦点强度最大值在 ${k_{\rm{L}}}$ 取1.8到2之间,并且强度随着 ${\;\beta _{\rm{L}}}$ 的增大而增大,但是并不明显,其中 ${k_{\rm{L}}} = 1.88$ 为极值点。所以,虽然更大的 ${\;\beta _{\rm{H}}}$ ${\;\beta _{\rm{L}}}$ 意味着更好的聚焦特性,但是基于实际应用考虑, ${\;\beta _{\rm{H}}}$ ${\;\beta _{\rm{L}}}$ 不能取值过大,故选取参数 ${k_{\rm{H}}} = 1.35$ ${\;\beta _{\rm{H}}} = $ $ 15$ ${k_{\rm{L}}} = 1.88$ ${\;\beta _{\rm{L}}} = 13$ 。(a),(b)对比可以发现,BCAB的聚焦特性比CAB增大了7.46倍,HCAB仅5.28倍,即为了获取更大的聚焦强度,BPF的调制效果比HPF要好。Deng等[17]发现HPF的调制效果要优越于BPF,与本文结果相反,这是由于波长的不同。与可见光相比,在THz波段大量能量聚集在倏逝波段,无法传播到远场,造成了能量的消耗。而BPF将能量从倏逝波向前转移,减少了能量的消耗,增大了聚焦特性,故在THz波段BPF更有效果。

图 6 不同参数下调制后模拟得到的突然自聚焦能力 Figure 6 Abruptly autofocusing capability obtained by simulation after modulation under different parameters

假设 ${I_{\max }}$ ${\left( {{I_0}} \right)_{\max }}$ 分别表示CAB,HCAB和BCAB任一横截面上的最大强度和初始平面的最大强度。突然自聚焦能力可通过 ${{{I_{\max }}} / {{{\left( {{I_0}} \right)}_{\max }}}}$ 分析,其随着z的变化而变化的 ${{{I_{\max }}}/ {{{\left( {{I_0}} \right)}_{\max }}}}$ 图7所示。CAB,HCAB和BCAB的焦距分别为12 mm,8.82 mm和9.18 mm。在聚焦之前,光轴处光强很小,但是在达到焦点附近光强突然增加,BCAB和HCAB在焦点处的光强可以达到初始面上CAB的100多倍,光束能量主要集中在中心光斑处,体现出特有的突然自聚焦特性,并且在聚焦了一段距离后,光强会很快衰减。其中BCAB的光束能量更为集中,拥有更好的自聚焦特性。

图 7 CAB(实线)、HCAB(点线)和BCAB(短横线)的自聚焦特性对比 Figure 7 Comparison of autofocusing characteristics between CAB(solid line), HCAB(dotted line), and BCAB(short dash line)

焦平面上的光场分布也很重要,因为它们可以直接显示焦点的尺寸和光强分布,图8绘制了焦平面上CAB,HCAB和BCAB的归一化强度分布。从图8可以看出,HCAB和BCAB的焦斑尺寸比CAB小。CAB的半高全宽为0.296 mm,HCAB的半高全宽为0.218 mm,BCAB的半高全宽为0.225 mm。通过调制,两种调制方法都减小了半高全宽,但HCAB的效果相较于BCAB要稍好,但不明显。

图 8 焦面上CAB(实线)、HCAB(点线)和BCAB(短横线)的归一化强度分布 Figure 8 Normalized strength distributions of CAB(solid line), HCAB(dotted line), and BCAB(short dash line)on focal plane

为了讨论波长分布对调制结果的影响,设一有限带宽光源,中心波长为0.4 mm,谱宽为0.2~0.6 mm。带通滤波的参数针对中心波长0.4 mm优化,图9给出了波长分别为0.2 mm、0.3 mm、0.4 mm、0.5 mm和0.6 mm的太赫兹波经带通调制后的结果,其中 ${\left( {{I_{{\rm{CAB}}}}} \right)_{\max }}$ ${\left( {{I_{{\rm{BCAB}}}}} \right)_{\max }}$ 分别表示CAB和BCAB的焦斑最大光强。从图中可以看出,虽然带通滤波器的参数只针对0.4 mm做了优化,但其他波长的波经频谱滤波后自聚焦能力也获得了较大幅度的提升,提升效果随着工作波长逐渐远离中心波长而下降。因此,这种频谱滤波方法对宽谱光源也能起到很好的提升效果。

图 9 宽谱中不同波长下BCAB的聚焦能力 Figure 9 Focusing ability of BCAB at different wavelengths in broad spectrum
3 结 论

本文研究了CAB在THz波段的传播特性,在傅里叶空间利用HPF和BPF对CAB进行调制,产生了HCAB和BCAB。并且通过数值仿真对结果进行对比,发现BPF确实可以增强CAB的聚焦特性并且比HPF的结果要好,证明了合适的参数下BPF调制后的聚焦强度可提高7倍,半高全宽也有改善,且可实现焦距的控制,并且对宽光谱的实际应用也具有现实意义,所以在THz波段带通滤波是一种提高自聚焦性能的有效方法。

参考文献
[1] SIVILOGLOU G A, BROKY J, DOGARIU A, et al. Observation of accelerating Airy beams[J]. Physical Review Letters, 2007, 99(21): 213901. DOI:10.1103/PhysRevLett.99.213901
[2] SIVILOGLOU G A, CHRISTODOULIDES D N. Accelerating finite energy Airy beams[J]. Optics Letters, 2007, 32(8): 979–981. DOI:10.1364/OL.32.000979
[3] MINOVICH A E, KLEIN A E, NESHEV D N, et al. Airy plasmons: non-diffracting optical surface waves[J]. Laser & Photonics Reviews, 2014, 8(2): 221–232.
[4] 钱义先, 毛红行. 可调谐广义艾里光束产生及传播特性[J]. 光子学报, 2018, 47(2): 0226001.
[5] 韩坤, 高爱华, 戢凯文, 等. 双复合离散型艾里涡旋加速光的传输特性[J]. 光子学报, 2018, 47(3): 0326001.
[6] 吴云龙, 聂劲松, 邵立, 等. 用傅里叶分析法研究艾里光束远场传播特性[J]. 光子学报, 2017, 46(3): 0314004.
[7] EFREMIDIS N K, CHRISTODOULIDES D N. Abruptly autofocusing waves[J]. Optics Letters, 2010, 35(23): 4045–4047. DOI:10.1364/OL.35.004045
[8] PAPAZOGLOU D G, EFREMIDIS N K, CHRISTODOULIDES D N, et al. Observation of abruptly autofocusing waves[J]. Optics Letters, 2011, 36(10): 1842–1844. DOI:10.1364/OL.36.001842
[9] KOULOUKLIDIS A D, PAPAZOGLOU D G, FEDOROV V Y, et al. Phase memory preserving harmonics from abruptly autofocusing beams[J]. Physical Review Letters, 2017, 119(22): 223901. DOI:10.1103/PhysRevLett.119.223901
[10] MANOUSIDAKI M, PAPAZOGLOU D G, FARSARI M, et al. Abruptly autofocusing beams enable advanced multiscale photo-polymerization[J]. Optica, 2016, 3(5): 525–530. DOI:10.1364/OPTICA.3.000525
[11] LIU K, KOULOUKLIDIS A D, PAPAZOGLOU D G, et al. Enhanced terahertz wave emission from air-plasma tailored by abruptly autofocusing laser beams[J]. Optica, 2016, 3(6): 605–608. DOI:10.1364/OPTICA.3.000605
[12] PANAGIOTOPOULOS P, PAPAZOGLOU D G, COUAIRON A, et al. Sharply autofocused ring-Airy beams transforming into non-linear intense light bullets[J]. Nature Communications, 2013, 4: 2622. DOI:10.1038/ncomms3622
[13] JIANG Y F, HUANG K K, LU X H. Radiation force of abruptly autofocusing Airy beams on a Rayleigh particle[J]. Optics Express, 2013, 21(20): 24413–24421. DOI:10.1364/OE.21.024413
[14] EFREMIDIS N K, PALTOGLOU V, VON KLITZING W. Accelerating and abruptly autofocusing matter waves[J]. Physical Review A, 2013, 87(4): 043637. DOI:10.1103/PhysRevA.87.043637
[15] HUANG L, GUO H L, LI J F, et al. Optical trapping of gold nanoparticles by cylindrical vector beam[J]. Optics Letters, 2012, 37(10): 1694–1696. DOI:10.1364/OL.37.001694
[16] LI N, JIANG Y F, HUANG K K, et al. Abruptly autofocusing property of blocked circular Airy beams[J]. Optics Express, 2014, 22(19): 22847–22853. DOI:10.1364/OE.22.022847
[17] DENG X Q, ZHENG G L, XU S X, et al. Fourier-space modification of a circular Airy beam and the creation of an optical needle[J]. Optics Communications, 2017, 402: 159–162. DOI:10.1016/j.optcom.2017.05.027
[18] GENG T, ZHANG X X. Propagation properties of the circular Airy beam with a Gaussian envelope in Fourier space[J]. Optics Express, 2020, 28(2): 2447–2455. DOI:10.1364/OE.384143
[19] HE J W, WANG S, XIE Z W, et al. Abruptly autofocusing terahertz waves with meta-hologram[J]. Optics Letters, 2016, 41(12): 2787–2790. DOI:10.1364/OL.41.002787