基于单一元件产生太赫兹无衍射波束技术研究.pdf

文章编号：1002-2082(2024)01-0033-06基于单一元件产生太赫兹无衍射波束技术研究韩美瑶1，郝丛静2，黄倩1，涂思语1，刘劲松1，杨振刚1，王可嘉1（1.华中科技大学武汉光电国家研究中心光学与电子信息学院，湖北武汉430074；2.中国航空工业集团公司北京长城航空测控技术研究所，北京100111）摘摘 要：要：提出了一种新型的光学元件，可以产生与元件具有一定距离的无衍射波束。普通轴棱锥的横截面为等腰三角形，通过在两个腰上限定抛物线进行元件设计，采用 3D 打印技术制作了不同结构参数的光学元件。角谱理论的仿真计算结果以及在 0.1THz 频率下进行的实验结果，都清晰地表明了产生的无衍射波束与元件之间具有一定距离。利用空间频谱对元件的波束行为进行分析，结果表明，该元件产生了一种新型的无衍射波束，根据元件的空频形状，将元件产生的无衍射波束命名为“梳状空频无衍射波束”。关键词：关键词：无衍射波束；抛物轴棱锥；角谱理论；空间频谱中图分类号：TN015文献标志码：ADOI：10.5768/JAO202445.0101005Terahertz diffraction-free beam generation technology based on single elementHANMeiyao1，HAOCongjing2，HUANGQian1，TUSiyu1，LIUJinsong1，YANGZhengang1，WANGKejia1（1.SchoolofOpticsandElectronicInformation,WuhanNationalResearchCenterofOpto-Electronics,HuazhongUniversityofScienceandTechnology,Wuhan430074,China；2.BeijingGreatWallAviationMeasurementandControlTechnologyResearchInstitute,AviationIndustryCorporationofChina,Beijing100111,China）Abstract：Anewopticalelementwasproposed,whichcouldgenerateadiffraction-freebeamwithacertaindistancefromtheelement.Thecrosssectionofanordinaryaxialpyramidwasanisoscelestriangleandtheelementwasdesignedbydefiningaparabolaonthetwowaists.The3Dprintingtechnologywasusedtofabricate optical elements with different structures.The results of simulation calculation using angularspectrum theory and the experiment conducted at 0.1 THz frequency clearly show that there is a certaindistancebetweenthegenerateddiffraction-freebeamandtheelement.Thebeambehavioroftheelementwasanalyzedaccordingtothespatialfrequencyspectrum.Theresultsshowthattheelementgeneratesanewtypeofdiffraction-freebeam,andaccordingtothespace-frequencyshapeoftheelement,thediffraction-freebeamgeneratedbytheelementisnamedcomb-likespace-frequencydiffraction-freebeam.Key words：diffraction-freebeam；parabolicpyramid；angularspectrumtheory；spatialfrequencyspectrum引言1987 年，DURNIN 首次提出了无衍射光束的概念1-2，自此无衍射光束成为光学领域研究的热点之一，各种不同的无衍射光束被产生出来，例如贝塞尔光束1、马丢光束3、艾里光束4以及抛物光束5等。理论上，无衍射光束是指，在传输过程中横截面光场分布保持不变的光束，即不会发生扩束现象。理想的无衍射光束是不可能实现的，因为它需要占据整个自由空间，且需要无穷的能量。因此在现实中，已产生出的无衍射光束均为收稿日期：2023-02-14；修回日期：2023-05-24基金项目：装备预研重点实验室基金（61422160107）作者简介：韩美瑶（1998），女，硕士研究生，主要从事太赫兹无衍射波束研究。E-mail：通信作者：王可嘉（1978），男，博士，副教授，主要从事太赫兹科学与技术研究。E-mail：第45卷第1期应用光学Vol.45No.12024年1月JournalofAppliedOpticsJan.2024准无衍射光束。由于具有无衍射1、自修复6和提供光学拉力7等诸多特性，无衍射光束已经广泛应用于光学操纵8-9、光学成像10-11、光通信12-13、精密测量14和激光加工15等相关领域。作为第一个被产生出的无衍射光束，贝塞尔光束是现在应用最多、最为著名的一类无衍射光束。迄今为止，人们已经开发出许多种产生贝塞尔光束的方法，例如环缝法1、法布里-珀罗腔16、轴棱锥法17、超表面法18-19、计算机全息技术20-22和空间光调制器23等。其中轴棱锥法具有制备方法简单、能量转换效率高等优点，被广泛使用。例如在红外、可见光波段，利用抛光工艺人们可以很轻易地获得玻璃轴棱锥24-25。对于微波、太赫兹波而言，使用增材制造技术可以快速制备出轴棱锥26-28，进而获得长波贝塞尔波束。实验结果表明，利用单一轴棱锥产生的贝塞尔光束总是紧贴轴棱锥。形像地说，对于轴棱锥而言是没有所谓的“焦距”概念，这一缺点限制了贝塞尔光束的一些实际应用，尤其是一些远程成像和检测应用等。为解决这一问题，2020 年 XIANG等人29提出了一种透镜轴棱锥级联结构，产生的米级长度太赫兹无衍射波束不再紧贴产生单元，而是具有百毫米量级的投送距离。需要指出的是，多级联透镜-轴棱锥结构器件数较多，器件间间距较大，因此损耗较大。2022 年，ZHANG 等人30研发了一种椭圆轴棱锥，产生出了具有“焦距”的太赫兹“贝塞尔波束”。借鉴 ZHANG 等人的设计理念，本文设计并制备出了抛物轴棱锥，同样也产生出了具有“焦距”的太赫兹无衍射波束。利用角谱理论和空间频谱对抛物轴棱锥进行分析，发现抛物轴棱锥与普通轴棱锥的空频形状不同，即抛物轴棱锥产生出的无衍射波束不是普通轴棱锥产生的贝塞尔波束，而是一种全新的无衍射波束梳状空频无衍射波束。1 理论设计与仿真计算2a=R/cos普通轴棱锥的横截面是一个等腰三角形，如图 1(a)所示。通过在等腰三角形的两个等边上限定抛物线设计轴棱锥结构，称其为抛物轴棱锥，如图 1(b)所示。等腰三角形的底边和底角分别为2R 和，两个腰的长度为。当底边 2R和底角 被确定时，两个腰的长度 a 也会被确定，此时改变抛物轴棱锥结构的唯一变量为 b。在此k=a/b引入一个无量纲量代替 b。R2Rb2a(a)普通轴棱锥(b)抛物轴棱锥(c)抛物轴棱锥坐标旋转(d)抛物轴棱锥相位延迟RLOOOzABC图 1 棱锥结构示意图Fig.1 Structure diagram of pyramid0=26 mmEin=exp(2/02)=x2+y2将光学元件放置于频率 0.1THz、束腰半径的准直高斯波束的束腰处，高斯波束可表示成：，其中为柱坐标系下的极轴。将光学元件看作薄器件，不考虑吸收损耗，则出射面的光场可表示为Eout=Einexp(j)（1）式中：j 为虚数单位；，为波束通过光学元件的相位延迟。如图 1(a)所示，普通轴棱锥的相位延迟为A=K(n1)(R)tan（2）K=2/式中：表示波矢；n 为材料折射率。接下来计算抛物轴棱锥的相位延迟。如图 1(c)所示，通过坐标系旋转和平移，将 坐标系转化为 z 坐标系。由于抛物轴棱锥是轴对称的，在此仅以 z 轴上半部分为例。坐标系逆时针旋转 后得到 坐标系：=cossinsincos=cos+sinsin+cos（3）坐标系通过平移转换为 z 坐标系，可表示为z=+acos+asin（4）34应用光学第45卷第1期将公式(3)带入公式(4)，整理可得 坐标系转化为 z 坐标系的转换公式：z=cos+sin+acossin+cos+asin（5）=(b/a2)2+b(A,zA)(A,A)AB(B,zB)(B,B)抛物线在 坐标系中的表达式为，设 A 点为抛物线上一动点，在 z 和 坐标系中的坐标分别为和。与 轴垂直，与 轴相交于 B 点，在 z 和 坐标系中的坐标分别为和。所以抛物轴棱锥的厚度函数为?AB?=zAzB=zA（6）将公式(5)带入公式（6），可得：?AB?=Asin+Acos+asin（7）=(a2/b)tan=(a2/b)tan2+b根 据 抛 物 线 参 数 方 程和，抛物轴棱锥的厚度函数可表示为?AB?=a2btansina2btan2cos+bcos+asin（8）0,arctan(b/a)(arctan(b/a),式中：参数。?AB?对 求导，即：?AB?=a2bsinsec22a2bcostansec2（9）?AB?/=0令，可得：tan=12tan（10）?AB?/=0?AB?当 满足公式(10)时，取最大值，即：?AB?max=a24btansin+bcos+asin（11）L=?AB?maxL=?AB?+?AC?如图 1(d)所示，设且，抛物轴棱锥右侧的竖直虚线为出射面，则抛物轴棱锥的相位延迟为=nK?AB?+K?AC?=K(n1)?AB?+KL RKL R（12）K=2/式中：表示波矢；n 为材料折射率。在此仅考虑相对的相位延迟，将公式(12)的上下两式同时减去 KL，再将公式(8)代入可得：=K(n1)(a2btansina2btan2cos+bcos+asin)（13）a=R/2cosk=a/b根据抛物轴棱锥的几何结构可知，代入（13）式可得：=K(n1)(tan+tan)tanRk2+R2k+Rtan2（14）0,arctan(1/k)arctan(1/k),式中。因此，当抛物轴棱锥的底边 2R 和底角 确定后，k 值是改变抛物轴棱锥相位延迟的唯一变量。2R=76.2 mm=10将公式(2)和公式(14)代入公式(1)，即可计算出普通轴棱锥和抛物轴棱锥的出射波束。两种光学元件的底边（3in），底角。光学元件利用 3D 打印技术制作，当频率为 0.1THz时 3D 打印材料的折射率 n 为 1.645。z=0 x=y=0利用角谱理论计算的仿真结果如图 2 所示。图 2(a)为利用公式(2)和公式(14)计算出的普通轴棱锥和抛物轴棱锥的相位分布；图 2(b)为普通轴棱锥和抛物轴棱锥的出射波束在 x-z 面的光强分布。其中对应光学元件的出射面，红色曲线为在时在 z 方向的光强变化，称其为中心光强曲线。由图 2(b1)可知，普通轴棱锥产生的无衍射波束紧贴元件，由图 2(b2)图 2(b4)可知，抛物轴棱锥的出射波束从圆环状波束逐渐转变为由一个主瓣和多个旁瓣组成的无衍射波束，即无衍射波束与元件具有一定的距离。随着 k 值的增大，5000普通轴棱锥普通轴棱锥211010归一化强度1010(b1)(a1)(a2)(a3)(a4)(b2)(b3)(b4)050500505005050050142k=2k=3k=5k=2k=3x/mmk=55888100129200z/mm3004000(a)相位分布(b)在x-z面的光强分布图 2 普通轴棱锥和抛物轴棱锥的仿真结果Fig.2 Simulation results of normal axial pyramid andparabolic pyramid应用光学2024，45（1）韩美瑶，等：基于单一元件产生太赫兹无衍射波束技术研究35k=2 k=3k=5抛物轴棱锥的无衍射波束中的旁瓣数量逐渐减少，能量逐渐被聚集在中心主瓣。在此，类比于透镜中焦距的定义，将抛物轴棱锥的出射面与中心光强曲线最大值处之间的距离定义为“焦距”，因此，和时抛物轴棱锥的焦距分别为58mm，88mm 和 129mm。2 频谱分析和实验验证f利用空间频谱进一步分析抛物轴棱锥产生的无衍射波束。普通轴棱锥和抛物轴棱锥都是中心对称的，所以其空间频谱也具有中心对称性，仅需观察空间频谱沿方向的分布即可，如图 3 所示。由图 3 可知，相比于普通轴棱锥的圆环状空频，抛物轴棱锥的空间频谱是由一个幅值最大的主频和以原点为圆心、幅值逐渐降低的侧边频构成。因此，抛物轴棱锥产生的波束并不是普通轴棱锥产生的贝塞尔波束，而是一种新型的无衍射波束，根据抛物轴棱锥的空频形状，将其所产生的波束命名为“梳状空频无衍射波束”。1.0轴棱锥k=2k=3k=50.80.6归一化强度0.40.200.30.20.1f/mm100.10.20.3k=2图 3 普通轴棱锥和、3、5 时抛物轴棱锥的空间频谱k=2Fig.3 Spatial frequency spectrum of normal axis pyramidand parabolic pyramid when,3,5=102R=76.2 mm k=2k=3以上仅为仿真计算和频谱分析结果，下面用实验验证抛物轴棱锥的波束行为。图 4(a)为实验装置示意图，太赫兹源为单频耿氏二极管（SPACEKLABS,GW-900），频率为 0.1THz，功率为 25mW。使用一个直径 76.2mm、焦距 100mm 的透镜（高密度聚乙烯材料）对太赫兹源发射的发散波束进行准直，准直后高斯波束的束腰半径为 26mm。抛物轴棱锥利用 3D 打印技术制作，打印精度为 0.01mm，打印材料为 LY1101，该材料在 0.1THz 频率下的折射率和吸收系数分别为 1.645 和 1.2cm1。选择底角，底边，、的抛物轴棱锥进行制作，图 4(b)为抛物轴棱锥的三维模型图和实物图。太赫兹探测器为肖特基二极管，将其固定在三维步进电机上进行单点逐行扫描。38.1 mm38.1 mmk=2太赫兹源 准直透镜(a)实验装置示意图(b)抛物轴棱锥的模型图和实物图抛物轴棱锥探测器k=2k=3k=3图 4 实验验证Fig.4 Experimental verificationk=2k=3z=0k=2z=42 mmk=2k=3抛物轴棱锥的实验结果如图 5 所示。由于仿真空间的有限性，实验结果与仿真计算结果存在一定的差别，但出射波束的光强分布仍吻合较好。图 5（a）和图 5（b）分别为和时抛物轴棱锥的光强分布。图 5 中对应抛物轴棱锥的出射面，红色曲线为中心光强曲线。由图 5 可知，抛物轴棱锥产生的梳状空频无衍射波束与出射面具有一定的距离。以的抛物轴棱锥为例，出射波束的起始阶段呈空心圆环状，中心光强在处取极小值，之后圆环状波束逐渐转变为由一个中心主瓣和多个旁瓣组成的无衍射波束，中心主瓣沿传播方向逐渐扩展，如图 5(a)所示。由图 5 可知，和时抛物轴棱锥的焦距分别为 154mm和 172mm。在利用角谱理论进行计算时，将抛物轴棱锥看做薄器件，不考虑其吸收损耗，但实验中抛物轴棱锥具有一定厚度且存在20005050050y/mmx/mm154 mmz=42 mm250 mm350 mm450 mm500400z/mm(b)k=3的抛物轴棱锥在x-y面和x-z面的光强分布600800归一化强度110020005050050y/mmx/mm172 mmz=58 mm250 mm350 mm450 mm500400z/mm600800归一化强度1100(a)k=2的抛物轴棱锥在x-y面和x-z面的光强分布图 5 抛物轴棱锥实验结果Fig.5 Experimental results of parabolic pyramid36应用光学第45卷第1期能量吸收，因此，对于抛物轴棱锥的焦距，仿真结果与实验结果存在一定误差。3 结论通过改变普通轴棱锥的几何结构，设计了一种新型的光学元件抛物轴棱锥。该元件可以产生具有“焦距”的无衍射波束，在 0.1THz 频率下利用 3D 打印技术制备了不同结构参数的抛物轴棱锥。仿真计算与实验结果均表明，抛物轴棱锥产生的无衍射波束与元件具有一定的距离。根据空间频谱，发现抛物轴棱锥产生了一种新型的无衍射波束，称其为“梳状空频无衍射波束”。这种利用抛物轴棱锥产生的“梳状空频无衍射波束”可应用于太赫兹成像领域，用于扩大成像系统的景深。参考文献：DURNINJ.Exactsolutionsfornondiffractingbeams:I.The scalar theoryJ.Journal of the Optical Society ofAmericaA，1987，4（4）：651-654.1DURNINJ,MICELIJrJJ,EBERLYJH.Diffraction-free beamsJ.Physical Review Letters，1987，58（15）：1499-1501.2GUTIRREZ-VEGA J C,ITURBE-CASTILLO M D,CHVEZ-CERDAS.Alternativeformulationforinvari-ant optical fields:Mathieu beamsJ.Optics Letters，2000，25（20）：1493-1495.3刘会龙,夏菁,陈懂,等.可调谐艾里光束及其传输特性J.应用光学，2016，37（6）：847-853.LIUHuilong,XIAJing,CHENDong,etal.TunableAirybeam and its propagation propertieJ.Journal of Ap-pliedOptics，2016，37（6）：847-853.4BANDRESMA,GUTIRREZ-VEGAJC,CHVEZ-CERDA S.Parabolic nondiffracting optical wavefieldsJ.OpticsLetters，2004，29（1）：44-46.5BOUCHALZ,WAGNERJ,CHLUPM.Self-reconstruc-tionofadistortednondiffractingbeamJ.OpticsCom-munications，1998，151（4/5/6）：207-211.6CHENJ,NGJ,LINZ,etal.OpticalpullingforceJ.NaturePhotonics，2011，5（9）：531-534.7ARLTJ,GARCS-CHVEZV,SIBBETTW,etal.Op-ticalmicromanipulationusingaBessellightbeamJ.Op-ticsCommunications，2001，197（4/5/6）：239-245.8ASHKINA,DZIEDZICJM,BJORKHOLMJE,etal.9Observationofasingle-beamgradientforceopticaltrapfor dielectric particlesJ.Optics Letters，1986，11（5）：288-290.HUY,CHENZ,XIANGL,etal.Extendeddepth-of-fieldall-opticalphotoacousticmicroscopywithadualnon-dif-fracting Bessel beamJ.Optics Letters，2019，44（7）：1634-1637.10XIANGF,YUM,ZHANGM,etal.Terahertzptycho-graphyusingalong-distancediffraction-freebeamastheprobeJ.OpticsLetters，2022，47（2）：238-241.11MPHUTHIN,GAILELEL,LITVINI,etal.Free-spaceoptical communication link with shape-invariant orbitalangular momentum Bessel beamsJ.Applied Optics，2019，58（16）：4258-4264.12王旭,王中宇,毛方儒,等.自由空间无衍射光通信模拟系统及其传输性能J.应用光学，2012，33（3）：644-648.WANGXu,WANGZhongyu,MAOFangru,etal.Free-spacecommunicationsimulationsystemanditstransmis-sion performance based on non-diffracting beamJ.JournalofAppliedOptics，2012，33（3）：644-648.13BAUERT,ORLOVS,PESCHELU,etal.Nanointerfero-metricamplitudeandphasereconstructionoftightlyfo-cusedvectorbeamsJ.NaturePhotonics，2014，8（1）：23-27.14BHUYANMK,VELPULAPK,COLOMBIERJP,etal.Single-shothighaspectratiobulknanostructuringoffusedsilica using chirp-controlled ultrafast laser BesselbeamsJ.Applied Physics Letters，2014，104（2）：021107.15COXAJ,DIBBLEDC.Nondiffractingbeamfromaspa-tially filtered Fabry Perot resonatorJ.Journal of theOpticalSocietyofAmericaA，1992，9（2）：282-286.16SCOTTG,MCARDLEN.Efficientgenerationofnearlydiffraction-freebeamsusinganaxiconJ.OpticalEngin-eering，1992，31（12）：2640-2643.17CHENWT,KHORASANINEJADM,ZHUAY,etal.Generation of wavelength-independent subwavelengthBessel beams using metasurfacesJ.Light:Science&Applications，2016，5（1）：39-44.18PFEIFFERC,GRBICA.ControllingvectorBesselbeamswith metasurfacesJ.Physical Review Applied，2014，2（4）：044012-1-7.19VIEIRATA,YEPESISV,SUAREZRAB,etal.Op-ticalreconstructionofnon-diffractingbeamsviaphotore-fractive holographyJ.Applied Physics B，2017，20应用光学2024，45（1）韩美瑶，等：基于单一元件产生太赫兹无衍射波束技术研究37123（4）：1-10.VASARAA,TURUNENJ,FRIBERGAT.Realizationofgeneralnondiffractingbeamswithcomputer-generatedhologramsJ.JournaloftheOpticalSocietyofAmericaA，1989，6（11）：1748-1754.21VIEIRATA,ZAMBONI-RACHEDM,GESUALDIMRR.Modelingthespatialshapeofnondiffractingbeams:ExperimentalgenerationofFrozenwavesviaholographicmethodJ.OpticsCommunications，2014，315：374-380.22ARRIZNV,MNDEZG,SNCHEZ-DE-LA-LLAVED.Accurate encoding of arbitrary complex fields withamplitude-onlyliquidcrystalspatiallightmodulatorsJ.OpticsExpress，2005，13（20）：7913-7927.23程治明,吴逢铁,范丹丹,等.绿光LED产生无衍射Bessel光束J.中国科学:物理学力学天文学，2012，42（8）：805-811.CHENGZhiming,WUFengtie,FANDandan,etal.Non-diffraction bessel-like beam generated by green LEDsourceJ.Science China:Physica,Mechanica&Astro-nomica，2012，42（8）：805-811.24卢文和,吴逢铁,郑维涛.透镜轴棱锥产生近似无衍射贝塞尔光束J.光学学报，2010，30（6）：1618-1621.25LUWenhe,WUFengtie,ZHENGWeitao.Generationofnon-diffractionBessel-likedbeamusingalensaxiconJ.ActaOpticaSinica，2010，30（6）：1618-1621.BITMAN A,MOSHE I,ZALEVSKY Z.Improvingdepth-offieldinbroadbandTHzbeamsusingnondiffract-iveBesselbeamsJ.OpticsLetters，2012，37（19）：4164-4166.26WUZ,WANGX,SUNW,etal.Vectorcharacterizationofzero-orderterahertzBesselbeamswithlinearandcircu-larpolarizationsJ.ScientificReports，2017，7（1）：1-11.27WEIX,LIUC,NIUL,etal.Generationofarbitraryor-derBesselbeamsvia3Dprintedaxiconsattheterahertzfrequency rangeJ.Applied Optics，2015，54（36）：10641-10649.28XIANGF,LIUD,XIAOL,etal.Generationofameter-scaleTHz diffraction-free beam based on multiple cas-cadedlens-axicondoublets:detailedanalysisandexperi-mentaldemonstrationJ.OpticsExpress，2020，28（24）：36873-36883.29ZHANGM,LIUW,HUANGQ,etal.Opticalelementtogeneratezero-orderquasi-Besselbeamwith“focallen-gth”J.OpticsLetters，2022，47（3）：553-556.3038应用光学第45卷第1期