1、文章编号：1002-2082(2023)05-1102-07被动锁模光纤激光器中高阶色散对输出脉冲特性影响陈炳焜，高博（吉林大学通信工程学院，吉林长春130012）摘摘 要：要：构建了被动锁模光纤激光器的理论模型，通过数值仿真从时域和频域两个方面研究了高阶色散对输出脉冲特性的影响，在时域方面，三阶色散导致脉冲形状发生畸变，脉冲时域轮廓不对称，时域偏移等现象；在频域方面，三阶色散导致边带位置及强度发生变化。同时，研究发现四阶色散可以展宽脉冲，改变脉冲的边带强度。本文的研究结果对被动锁模掺镱光纤激光器的实际应用具有指导意义。关键词：关键词：激光器；被动锁模激光器；三阶色散；四阶色散中图分类号：TN
2、248文献标志码：ADOI：10.5768/JAO202344.0507001Influence of high-order dispersion on output pulse characteristics in passivelymode-locked fiber laserCHENBingkun，GAOBo（CollegeofCommunicationEngineering,JilinUniversity,Changchun130012,China）Abstract：Atheoreticalmodelofpassivelymode-lockedfiberlaserwasconstruc
3、ted,andtheinfluencesofhigh-orderdispersiononoutputpulsecharacteristicswerestudiedbynumericalsimulationinbothtimedomainandfrequencydomain.Intermsoftimedomain,thethird-orderdispersionresultedindistortionofpulseshape,asymmetryofpulsetime-domainprofileandtemporalshift.Intermsoffrequencydomain,thethird-o
4、rderdispersioncausedchangesinpositionandintensityofsidebands.Atthesametime,itwasfoundthatthefourth-orderdispersioncouldbroadenthepulseandchangesidebandsintensityofthepulse.Theresultsareinstructiveforthepracticalapplicationsofthepassivelymode-lockedytterbium-dopedfiberlasers.Key words：lasers；passivel
5、ymode-lockedlaser；third-orderdispersion；fourth-orderdispersion引言超短脉冲目前已广泛应用于材料加工1-3、光学频率梳4-6、传感器7-8、光谱学9-10和其他领域11-12。随着超快光学领域的快速发展，被动锁模光纤激光器因其结构紧凑、脉冲稳定和成本低而得到广泛研究13-15。先前的理论和实验结果表明，高阶色散是影响激光器输出脉冲的重要因素。三阶色散(third-orderdispersion，TOD)会导致脉冲展宽或压缩、脉冲失真、脉冲混沌等现象。三阶色散可以增加锁模脉冲的持续时间16，其能够压缩脉冲是因为相速度和群速度匹配的色散波
6、与锁模脉冲的共振作用17。由于三阶色散的作用，脉冲不对称并且中心波长发生偏移18。同时，研究人员也发现三阶色散可以抑制脉冲不稳定性，它可以改变主脉冲与群速度和相速度匹配色散波之间的相互作用19。在近零色散点附近，通过增大三阶色散的参数值，发现了脉冲混沌状态20。一些有关四阶色散(fourth-orderdispersion，FOD)的研究也大量展开。研究人员发现四阶色散和自相位调制作用可以产生孤子，同时四阶色散可以导致脉冲展宽21。RUNGEAFG 等人发现纯四次孤子脉冲，脉冲的能量大幅收稿日期：2023-02-27；修回日期：2023-04-11基金项目：吉林省教育厅科技项目（JJKH202
7、31170KJ）作者简介：陈炳焜（1998），女，硕士研究生，主要从事光纤激光器研究。E-mail：通信作者：高博（1980），男，博士，副教授，博导，主要从事光纤激光器研究。E-mail：第44卷第5期应用光学Vol.44No.52023年9月JournalofAppliedOpticsSep.2023度提升22。本文在被动锁模掺镱光纤激光器中对高阶色散参数进行了理论和仿真研究。首先，介绍了被动锁模光纤激光器的模型和相应的方程。其次，分析不同三阶色散值下孤子特性，研究三阶色散对孤子特性的影响。时域方面，三阶色散造成脉冲轮廓不对称及时域偏移等现象；光谱方面，三阶色散导致边带位置及强度发生变化。
8、最后，研究了四阶色散对被动锁模光纤激光器输出脉冲特性的影响。通过分析脉冲时域演化图及光谱图，发现脉冲宽度随四阶色散发生改变，造成脉宽展宽，也对光谱边带产生影响。纯四次孤子具有较高的能量，而本文关于高阶色散的研究对未来研究纯四次孤子有一定意义。1 数值仿真模型被动锁模掺镱激光器的腔设计如图 1 所示。掺镱光纤(Yb-dopedgainfiber，YDF)由泵浦二极管泵浦，在可饱和吸收体(saturableabsorber，SA)之后，使用耦合器输出一部分脉冲。另一部分光脉冲继续传输，通过色散补偿光纤(dispersioncomp-ensationfiber，DCF)、单 模 光 纤(single
9、-modefiber，SMF)，再进入光隔离器，从而形成闭合环路。YDFWDMSMFDCF输出耦合器隔离器泵浦SA图 1 被动锁模掺镱光纤激光器示意图Fig.1 Schematic of passively mode-locked ytterbium-dopedfiber laser被动锁模掺镱光纤激光器可以由非线性薛定谔方程 NLSE 建模，如下所示：AZ=g2Ai22U2T2+36U3T3+i|A|2A（1）A(z,t)23gg(w)式中：是 NLSE 描述的脉冲时间复包络；是非线性系数；是群速度色散；是三阶色散参数；是掺杂光纤的增益系数，可以用抛物线形增益谱表示：g(w)=g01+(w/
10、g)2+Epulse/Esat（2）g0gEsatEpulse=TR/2rTR/2|A(z,T)|2dTTR式中：是小信号增益；是增益带宽；是增益饱和能量；阐述了脉冲能量。是仿真时域窗口，使用的可饱和吸收体模型如下：R=Runsat+Rsat111+P()Psat（3）RsatRunsatP()Psatg0gEsatRunsatRsatPsat222式中：是饱和反射系数；是非饱和反射系数；是光脉冲时域中每点对应的瞬时光功率；是可饱和吸收体的饱和功率。该模型采用分布傅里叶算法求解。下面的数值用于本次仿真：为8m1，为 0.5W1km1，为 15Thz，为 20pJ；为 0.6，为 0.3，为 2
11、.6625W；耦合器分光比为 2080；掺镱光纤长度为 0.25m，为 25ps2/km；色散补偿光纤长度为 3.5m，为25ps2/km；单模光纤长度为 6m，为 25ps2/km。仿真中心波长为1030nm，这与掺镱光纤的实际波长一致。2 仿真分析与讨论2.1 三阶色散对孤子特性影响的数值仿真iUZ=22U2T2+i63UT3U(z,T)=12+rU(0,w)exp(i22w2z+i63w3ziwT)dw当腔内群速度色散接近零，三阶色散是一个重要因素。为了研究三阶色散对输出脉冲的影响，保持其他参数不变，线性减小三阶色散从 0到15ps3/km，间隔为0.2ps3/km。通过线性降低三阶色散
12、数值，观察到脉冲的时域轮廓变得不对称，形状出现畸变。当三阶色散数值降到6ps3/km 时，脉冲出现不对称的震荡拖尾结构。这是因为忽略非线性效应，只考虑色散效应，可以解 NLSE 方程：，并 得 到：，所以只考虑2这一参数，脉冲时域轮廓不变。当考虑到 3时，脉冲时域轮廓发生畸变，并在脉冲前沿或者后沿产生不对称震荡结构。脉冲的时域轮廓如图 2 所示。1.000.750.500.250505时间/ps强度/(a.u.)1015201510时间/ps强度/(a.u.)55时间/ps强度/(a.u.)5时间/ps强度/(a.u.)(a)3=1 ps3/km(b)3=3.6 ps3/km(c)3=6 ps
13、3/km(d)3=15 ps3/km应用光学2023，44（5）陈炳焜，等：被动锁模光纤激光器中高阶色散对输出脉冲特性影响11035时间/ps强度/(a.u.)1.000.750.500.2501510时间/ps强度/(a.u.)5051.000.750.500.25050时间/ps强度/(a.u.)5101.000.750.500.25050时间/ps强度/(a.u.)51015(a)3=1 ps3/km(b)3=3.6 ps3/km(c)3=6 ps3/km(d)3=15 ps3/km图 2 不同三阶色散下，脉冲时域轮廓Fig.2 Pulse temporal profile under
14、different third-orderdispersion此外，当三阶色散为正时，震荡发生在脉冲的后沿；当三阶色散为负时，震荡发生在脉冲的前沿，如图 3 所示。数值仿真结果表明，当三阶色散从6ps3/km 变化到15ps3/km 时，脉冲持续时间从 0.78ps 到 1.1ps 不断增加，而脉冲峰值功率从 31.1307W减小到 24.3268W，这影响了脉冲的传输质量。继续观察在三阶色散影响下的脉冲频谱，将三阶色散从1ps3/km 降低到15ps3/km，间隔为0.2ps3/km。最初脉冲顶部开始出现倾斜，如图 4(a)所示。当三阶色散降低到3.4ps3/km，脉冲出现凯利边带。随着三阶
15、色散继续减小，脉冲频谱变得更加不对称，频谱形状近似三角形，如图 4(b)所示。1.000.750.500.2505时间/ps强度/(a.u.)05101.000.750.500.25010时间/ps强度/(a.u.)505(a)3=15 ps3/km(b)3=15 ps3/km图 3 三阶色散绝对值相同时，脉冲时域轮廓Fig.3 Pulse temporal profile with the same absolute valueof third-order dispersion1.00.60.80.40.201 0151 0201 0251 0301 0351 0401 045波长/nm强度
16、/(a.u.)1.00.60.80.40.201 0201 0251 0301 0351 040波长/nm强度/(a.u.)(a)3=1 ps3/km(b)3=6 ps3/km图 4 不同三阶色散下，脉冲频谱Fig.4 Pulse frequency spectrum under different third-order dispersion1104应用光学第44卷第5期图 5 记录了边带的位置分布和强度。当三阶色散从3.4ps3/km 降低到15ps3/km，即从右向左读取横坐标。左边带位置和主峰位置从 1030nm左侧接近 1030nm，右边带位置从 1030nm 右侧接近1030nm。
17、图 5 在三阶色散为6ps3/km 这一点时，变化很大，推测这与该点的脉冲形状变化有关。在3.4ps3/km6ps3/km 内，左边带强度持续增加，然后在6ps3/km15ps3/km 内，左边带强度持续降低。在整个区域内，右边带强度持续降低。左右边带强度的差异可以解释如下：在三阶色散的影响下，边带频率相对频率中心不对称，再结合有限增益带宽的影响，导致一侧的边带强度更强23。此外，当三阶绝对值相同时，脉冲频谱关于1030nm 镜像对称。1 0281 0271 0261 0251 0241 0231 022141512108643.4三阶色散3/(ps3km1)(a)左边带位置随三阶色散的变化(
18、b)左边带强度随三阶色散的变化(c)右边带位置随三阶色散的变化(d)右边带强度随三阶色散的变化(e)主峰位置随三阶色散的变化左边带位置/nm0.50.40.30.20.10141512108643.4三阶色散3/(ps3km1)左边带强度/(a.u.)1 033.51 033.01 032.51 032.0141512108643.4三阶色散3/(ps3km1)右边带位置/nm1.00.80.60.40.2141512108643.4三阶色散/(ps km1)右边带强度/(a.u.)1 029.51 029.01 028.51 028.01 027.51 027.0141512108643.4
19、三阶色散/(ps km1)主峰位置/nm141512108643.4三阶色散/(ps km1)(a)左边带位置随三阶色散的变化(b)左边带强度随三阶色散的变化(c)右边带位置随三阶色散的变化(d)右边带强度随三阶色散的变化(e)主峰位置随三阶色散的变化左边带位置/nm0.50.40.30.20.1141512108643.4三阶色散/(ps km1)左边带强度/(a.u.)1 033.51 033.01 032.51 032.0141512108643.4三阶色散/(ps km1)右边带位置/nm1.00.80.60.40.20141512108643.4三阶色散3/(ps3km1)右边带强度
20、/(a.u.)1 029.51 029.01 028.51 028.01 027.51 027.0141512108643.4三阶色散3/(ps3km1)主峰位置/nm图 5 脉冲边带的位置和强度随三阶色散的变化Fig.5 Position and intensity of pulse sidebands vary withthird-order dispersion观察三阶色散下的脉冲演化，可以发现三阶色散的一个作用是提前或者延迟脉冲，使得脉冲在时域发生偏移24。在这里定义：脉冲时域偏移=偏移时间/循环数，其符号为正时，脉冲在时域向右移动，符号为负时，脉冲在时域向左移动。当三阶色散从1ps3
21、/km 变化到15ps3/km 时，观察到脉冲在时域向左偏移，并且偏移程度越来越大，如图 6 所示。在3.4ps3/km 和6ps3/km 的两点变化可能归因于脉冲形状的变化。0.100.10.20.30.40.5151050三阶色散3/(ps3km1)时域偏移/(ps/循环)图 6 三阶色散对脉冲偏移的影响Fig.6 Influence of third-order dispersion on pulse shift应用光学2023，44（5）陈炳焜，等：被动锁模光纤激光器中高阶色散对输出脉冲特性影响11052.2 四阶色散对孤子特性影响的数值仿真进一步研究四阶色散对脉冲的影响，4是四阶色散
22、参数。考虑到三阶色散比四阶色散对脉冲的影响更大，让 3参数为零，去消除三阶色散的影响，2以及其他参数仍然保持最初仿真数值。线性降低 4数值，从10ps4/km 降低到100ps4/km，间隔为10ps4/km。四阶色散参数和光脉冲传输方程参照相关文献22,25。当考虑到四阶色散时，光脉冲传输方程如下：AZ=g2Ai22U2T2+i424U4T4+i|A|2A（4）四阶色散可以使脉冲展宽21。这是因为负四阶色散使脉冲的蓝移快于红移，与群速度色散使脉冲展宽的原理相似。当四阶色散从10ps4/km降低到100ps4/km 时，脉宽从 0.62ps 展宽到 1.42ps，如图 7 所示。1.61.41
23、.21.00.80.6100806040200四阶色散/(ps4km1)脉冲宽度/ps图 7 脉宽时间随四阶色散的变化Fig.7 Variation of pulse width time with fourth-orderdispersion图 8 记录了四阶色散为20ps4/km 情况下，脉冲的演化过程以及频谱。在小饱和能量 Esat为 20pJ 的条件下，受到四阶色散影响的脉冲频谱依旧呈对称结构，这与三阶色散导致脉冲频谱的不对称有所不同。同时脉冲也未产生长拖尾结构，也未发生时域偏移。1 000(a)60080040020001050510时间/ps(a)脉冲演化图(b)脉冲频谱循环圈数(
24、b)1 0251 035波长/nm强度/(a.u.)(a)1055时间/ps(a)脉冲演化图(b)脉冲频谱循环圈数1.0(b)0.60.80.40.201 0201 0251 0301 0351 040波长/nm强度/(a.u.)4=20 ps4km图 8 当时的脉冲演化图与脉冲频谱4=20 ps4kmFig.8 Pulse evolution diagram and pulse spectrum while在四阶色散影响下，脉冲的频谱呈现对称结构的凯利边带，如图8(b)所示。当四阶色散从10ps4/km变化到100ps4/km 时，记录边带强度的变化。随着四阶色散的减小，边带强度逐渐减小，变
25、化趋势如图 9 所示。0.80.60.40.2010080200四阶色散4/(ps4km1)强度/(a.u.)6040图 9 边带强度随四阶色散的变化Fig.9 Variation of sideband intensity with fourth-orderdispersion3 结论本文研究了被动锁模掺镱光纤激光器中高阶色散对输出脉冲的影响。首先，探索了三阶色散对输出孤子脉冲特性的影响。通过数值仿真，发现三阶色散在时域上改变了脉冲时域轮廓，并可以提前或者延迟脉冲。在频域上，三阶色散造成了频谱的不对称，影响了脉冲频谱边带的位置和强度。最后，探索了四阶色散对锁模脉冲特性的影响，给出孤子脉冲的时
26、域演化图及频谱图。仿真过程中，脉冲未产生时域偏移和长拖尾结构，这与三阶色散有所不同。通过仿真发现，脉冲宽度随四阶色散的变化发生改变，四阶色散类似于1106应用光学第44卷第5期群速度色散，也可以展宽脉冲，同时引起频谱边带的强度发生改变。四阶色散可以产生高能的纯四次孤子，关于四阶色散的研究对理解纯四次孤子有一定的指导意义。这项工作有利于更好地理解高阶色散，对光纤激光器的实际应用具有重要意义。参考文献：SUGIOKA K,CHENG Y.Ultrafast lasersreliabletools for advanced materials processingJ.Light:Sci-ence&A
27、pplications，2014，3（4）：e149.1GATTASSRR,MAZURE.Femtosecondlasermicroma-chining in transparent materialsJ.Nature Photonics，2008，2（4）：219-225.2刘友强,曹银花,潘飞,等.激光加工用半导体激光器的光束变换J.光学精密工程，2012，20（3）：455-461.LIUYouqiang,CAOYinhua,PANFei,etal.Beamtrans-formationofdiodelasersusedinlaserprocessingJ.Op-ticsandPreci
28、sionEngineering，2012，20（3）：455-461.3赵力杰,周艳宗,夏海云,等.飞秒激光频率梳测距综述J.红外与激光工程，2018，47（10）：1006008.ZHAOLijie,ZHOUYanzong,XIAHaiyun,etal.Over-viewofdistancemeasurementwithfemtosecondopticalfrequency combJ.Infrared and Laser Engineering，2018，47（10）：1006008.4谢建东,严利平,陈本永,等.可调谐激光器激光波长宽范围自动偏频锁定J.光学精密工程，2021，29（2）
29、：211-219.XIEJiandong,YANLiping,CHENBenyong,etal.Auto-matic offset-frequency locking of external cavity diodelaserinwidewavelengthrangeJ.OpticsandPrecisionEngineering，2021，29（2）：211-219.5高宇炜,方守龙,武腾飞,等.双飞秒激光频率梳光谱测量技术研究进展J.应用光学，2021，42（1）：157-175.GAOYuwei,FANGShoulong,WUTengfei,etal.Re-search progress
30、of double femtosecond laser frequencycombspectroscopymeasurementtechnologyJ.JournalofAppliedOptics，2021，42（1）：157-175.6KAMATAM,OBARAM,GATTASSRR,etal.Opticalvibrationsensor fabricated by femtosecond laser mi-cromachiningJ.AppliedPhysicsLetters，2005，87（5）：1442.7ZHANGC,JIANGZS,FUSN,etal.Femtosecondlase
31、renabled selective micro-holes drilling on the multicore-8fiberfacetfordisplacementsensorapplicationJ.OpticsExpress，2019，27（8）：10777.XUC,WISEFW.RecentadvancesinfibrelasersfornonlinearmicroscopyJ.NaturePhotonics，2013，7（11）：875-882.9OZEKIY,ASAIT,SHOUJW,etal.Multicolorstimu-latedRamanscatteringmicrosco
32、pywithfastwavelength-tunableYbfiberlaserJ.IEEEJournalofSelectedTop-icsinQuantumElectronics，2018，25（1）：1-11.10SCHWEITZER C,BREZIN A,COCHENER B,etal.Femtosecond laser-assisted versus phacoemulsificationcataract surgery(FEMCAT):a multicentre participant-maskedrandomisedsuperiorityandcost-effectivenesst
33、ri-alJ.TheLancet，2020，395（10219）：212-224.11POPMINTCHEVT,CHENMCH,POPMINTCHEVD,etal.BrightcoherentultrahighharmonicsinthekeVX-rayregimefrommid-infraredfemtosecondlasersJ.Sci-ence，2012，336（6086）：1287-1291.12况庆强,桑明煌,聂义友,等.被动锁模掺铒光纤激光器中的有理数谐波锁模J.光学精密工程，2009，17（11）：2719-2723.KUANGQingqiang,SANGMinghuang,NI
34、EYiyou,etal.Rational harmonic mode-locking of passively mode-lockederbium-dopedfiberlaserJ.OpticsandPrecisionEngineering，2009，17（11）：2719-2723.13TANGDY,ZHAOLM,ZHAOB,etal.Mechanismofmultisolitonformationandsolitonenergyquantizationinpassively mode-locked fiber lasersJ.Physical ReviewA，2005，72（4）：0438
35、16.14赵德双,刘永智,黄绣江,等.飞秒被动锁模光纤激光器的稳定性研究J.应用光学，2005，26（4）：26-29.ZHAODeshuang,LIUYongzhi,HUANGXiujiang,etal.Stability of passively mode-locked femtosecond fiberlaserJ.JournalofAppliedOptics，2005，26（4）：26-29.15HAUS H A,MOORES J D,NELSON L E.Effect ofthird-orderdispersiononpassivemodelockingJ.Op-ticsLetter
36、s，1993，18（1）：51-53.16LOGVINY,KALOSHAVP,ANISH.Third-orderdis-persion impact on mode-locking regimes of Yb-dopedfiber laser with photonic bandgap fiber for dispersioncompensationJ.OpticsExpress，2007，15（3）：985.17DENNISML,DULINGIN.Third-orderdispersioninfemtosecond fiber lasersJ.Optics Letters，1994，19（2
37、1）：1750-1752.18LOGVINY,ANIS H.Suppression of multi-period in-19应用光学2023，44（5）陈炳焜，等：被动锁模光纤激光器中高阶色散对输出脉冲特性影响1107stabilitiesbythird-orderdispersioninmode-lockedYb-dopedfiberlasersJ.JournaloftheOpticalSocietyofAmericaB，2008，25（4）：622-632.SAKAGUCHIH,SKRYABINDV,MALOMEDBA.Stationary and oscillatory bound s
38、tates of dissipativesolitonscreatedbythird-orderdispersionJ.OpticsLet-ters，2018，43（11）：2688-2691.20BLANCO-REDONDOA,DESTERKECM,SIPEJE,etal.Pure-quartic solitonsJ.Nature Communications，2016，7：10427.21RUNGEAFJ,HUDSONDD,TAMKKK,etal.The22pure-quartic soliton laserJ.Nature Photonics，2020，14（8）：492-497.OLI
39、VIERM,ROYV,PICHM.Third-orderdispersionandboundstatesofpulsesinafiberlaserJ.OpticsLet-ters，2006，31（5）：580-582.23PICH M,CORMIER J F,ZHU X N.Bright opticalsolitoninthepresenceoffourth-orderdispersionJ.Op-ticsLetters，1996，21（12）：845-847.24ZHANGZX,LUOM,CHENJX,etal.Pulsatingdy-namicsinapure-quarticsolitonfiberlaserJ.OpticsLet-ters，2022，47（7）：1750-1753.251108应用光学第44卷第5期
