-
自相似脉冲可由光纤激光器、色散渐减光纤或非线性渐增光纤产生,本文中采用掺铒锁模激光器来实现自相似脉冲的演化。分别构建激光器和MZI的数值分析模型,便于分析脉冲在演化和压缩时的传输特性。
-
自相似锁模激光器能够将任何类型的输入脉冲演化成自相似锁模脉冲。这种激光器由抽运光源、掺铒增益光纤(erbium-doped gain fiber, EDF)、滤波器、波分复用耦合器(wave division multiplexing coupler, WDM)、可饱和吸收体(saturable absorber, SA)、耦合器(output coupler, OC)、色散补偿光纤(dispersion compensation fiber, DCF)及单模光纤(single-mode fiber, SMF)组成,如图 1所示。
增益光纤采用掺铒光纤,脉冲在增益光纤中的传输方程表示为[22]:
$ \begin{gathered} \frac{\partial u(z, T)}{\partial z}=-\mathrm{i} \frac{\beta_{\mathrm{EDF}, 2}}{2} \frac{\partial^2 u(z, T)}{\partial T^2}+ \\ \mathrm{i} \gamma_{\mathrm{EDF}}|u(z, T)|^2 u(z, T)+ \\ \frac{g_{\text {pulse }}}{2}\left[u(z, T)+T_2{ }^2 \frac{\partial^2 u(z, T)}{\partial T^2}\right] \end{gathered} $
(1) 式中,z为光脉冲在光纤中的传播距离,u(z, T)为脉冲包络的慢变幅值,时间量度T在一参照系中随脉冲以群速度vg(T=t-z/vg)运动而测量, t为时间, βEDF, 2为增益光纤的2阶色散,βEDF, 2=23 ps2/km,γEDF为增益光纤的非线性系数,γEDF=4.7 W-1·km-1,gpulse为增益总量,gpulse=g0/(1+Epulse/Esat), 小信号增益系数g0=30.4 /m,Epulse为光脉冲的能量,Esat为增益饱和能量,Esat=60 pJ,增益带宽时间T2=λ2/(2πcΔλ),λ为输入脉冲的波长,c为光速,Δλ为增益带宽。
耦合器将一部分的脉冲输出,输出的比例为10%。
可饱和吸收体透射率表示为:
$ T_{\mathrm{sa}}=1-L_0 \times\left(1+\frac{P}{P_{\mathrm{sat}}}\right) $
(2) 式中,P为脉冲的瞬时功率,Psat是可饱和吸收体的饱和功率,Psat=1500 W,L0为调制深度,L0=0.73。
色散补偿光纤起改变腔内净色散的作用,可用非线性薛定谔方程表示[23]:
$ \begin{gathered} \frac{\partial u(z, T)}{\partial z}=-\frac{\mathrm{i} \beta_{\mathrm{DCF}, 2}}{2} \frac{\partial^2 u(z, T)}{\partial T^2}+ \\ \mathrm{i} \gamma_{\mathrm{DCF}}|u(z, T)|^2 u(z, T) \end{gathered} $
(3) 式中,βDCF, 2、γDCF分别为色散补偿光纤对应的2阶色散系数和非线性系数,βDCF, 2=-130 ps2/km,γDCF=1.3 W-1·km-1。
脉冲在SMF中的演化不考虑光纤的损耗和高阶色散,采用非线性薛定谔方程来表示:
$ \begin{gathered} \frac{\partial u(z, T)}{\partial z}+\frac{\mathrm{i} \beta_2}{2} \frac{\partial^2 u(z, T)}{\partial T^2}= \\ \mathrm{i} \gamma|u(z, T)|^2 u(z, T) \end{gathered} $
(4) 式中,β2为单模光纤的2阶色散系数, γ为单模光纤的非线性系数。对公式进行分步傅里叶法来研究激光器输出脉冲在SMF中传输特性。
将(4)式改写成如下形式:
$ \frac{\partial u(z, T)}{\partial z}=(\widehat{D}+\widehat{N}) u(z, T) $
(5) 式中,$\widehat{D}$是微分算符,表示线性介质的色散,$widehat{N}$是非线性算符,决定脉冲传输过程中光纤非线性效应的影响。
$ \widehat{D}=-\frac{\mathrm{i} \beta_2}{2} \frac{\partial^2}{\partial T^2} $
(6) $ \widehat{N}=\mathrm{i} \gamma|u(z, T)|^2 $
(7) -
MZI由两种单模光纤组成(如图 2所示),自相似脉冲经过第1个耦合器后被分成两个脉冲,与非线性环形镜不同,两个分脉冲分别进入不同的物理路径,通过改变双臂的参数可以破坏MZI的平衡,从而得到自相位调制感应相移,实现光开关功能。
本文中设计了一个非对称的MZI,输入脉冲经过第1个耦合器后进入上下臂的两路脉冲的光信号振幅为u1和u2如下式所示[24]:
$ u_1=\sqrt{p_1} u_{\text {in }} $
(8) $ u_2=\mathrm{i} \sqrt{1-p_1} u_{\text {in }} $
(9) 式中,p1为第1个耦合器的功分比,uin为自相似脉冲包络的慢变振幅。两脉冲在MZI中传输并于第2个耦合器处发生干涉,两个输出端口的振幅u3和u4分别为:
$ u_3=\sqrt{p_2} u_1{ }^{\prime}+\mathrm{i} \sqrt{1-p_2} u_2{ }^{\prime} $
(10) $ u_4=\mathrm{i} \sqrt{1-p_2} u_1{ }^{\prime}+\sqrt{p_2} u_2{ }^{\prime} $
(11) 式中,u1′为上臂输出脉冲的振幅,u2′为下臂输出脉冲的振幅,u3为输出的压缩脉冲的振幅,u4则为脉座的振幅,p2为第2个耦合器的功分比。
MZI的上臂采取SMF级联的方法,能够使脉冲在保证峰值功率增大的前提下,减小脉宽。关于SMF1和SMF2的最佳长度选取规则如下:上臂的脉冲在到达第2个耦合器时的脉宽最窄,峰值功率最高;下臂的光纤长度要使第2个耦合器干涉后输出的脉座最小。
-
采用低功率的高斯脉冲作为抽运光源,为了能够更好地反馈脉冲在自相似激光器中的演化过程,抽运光源的中心波长采用光纤通信的常用波长。将一个半峰全宽(full width at half maximum, FWHM)为2.8 ps、峰值功率为10-10 W以及中心波长为1550 nm的高斯脉冲输入到自相似锁模激光器,循环30圈后获得自相似脉冲输出,如图 3所示。
从图 3a可以看出, 高斯脉冲在激光器的前几圈演化过程中,由于初始的峰值功率较小,波形图变化不明显,随后脉冲峰值功率突然急剧增加后降低,这是因为脉冲在循环的过程中受到大量自相位调制感应频率啁啾的作用,较弱的色散效应会引起脉冲整形,整个激光器腔的净色散呈正色散,脉冲展宽,脉冲的前后沿变陡。随着循环次数的增加,脉冲的前后沿逐渐变缓,脉冲趋于稳定,达到自相似锁模脉冲的输出条件。图 3b所示的啁啾线性部分也是突然减小后增大的过程。从图 3c脉冲的频谱图可以看出,相比于前5圈,第6和第7圈频谱的强度突然增大,由于入射的高斯脉冲为无啁啾的,随着圈数的增加,自相位调制产生的频率分量展宽了频谱,频谱的振荡峰的个数增加,且最外层的峰的强度最大。在第21圈后,从时域图可以看出脉冲波形呈抛物线型,脉冲频谱与脉冲的时域形状相关,振荡的结构变得不明显, 最后几圈的频谱趋于一个强度较低且稳定的抛物线型。最终从激光器中获得了峰值功率和半峰全宽分别为26.9941 W和4.044 ps的自相似脉冲。
进一步引入自相似因子S来衡量激光器输出的脉冲是否符合自相似脉冲的演化要求,如下式所示[25]:
$ S=\frac{\int_{-\infty}^{\infty}|u(z, T)|^2-\left|u^{\prime}(z, T)\right|^2 \mathrm{~d} T}{\int_{-\infty}^{\infty}\left|u^{\prime}(z, T)\right|^2 \mathrm{~d} T} $
(12) 式中,u(z, T)为数值仿真的激光器输出脉冲的包络振幅,u′(z, T)为具有与u(z, T)相同的峰值功率和半峰全宽的抛物线脉冲。根据多次仿真模拟,当S≤0.1时, 输出脉冲符合自相似脉冲的演化要求。经过计算,输出脉冲的自相似因子为0.04,符合自相似脉冲的要求。最终输出的自相似脉冲和抛物线脉冲的波形及啁啾如图 3d所示。
-
采用两种单模光纤来构建MZI,两种单模光纤的区别在于2阶色散的值,SMF1和SMF2的β2分别为-20 ps2/km和-1.25 ps2/km,两个耦合器的功分比均为0.55。
-
分脉冲u1在上臂上传输是为了获得峰值功率较高和脉宽较窄的压缩脉冲,确定了耦合器的功分比,将u1在SMF1中传输,可以获得传输过程的时域图。如图 4a所示,脉冲的脉宽随着传输距离的增大而逐渐变窄,峰值功率在z=8.16 m达到最高后急剧下降。
引入基座能量比例Pe来评价脉冲基座的大小,基座能量比例Pe是传输脉冲的总能量与双曲正割脉冲总能量的相对差值,双曲正割脉冲是与输出脉冲具有
相同峰值功率和的无基座脉冲。基座能量比例的表达式如下式所示[26]:
$ P_{\mathrm{e}}(\%)=\frac{\left|E_{\text {total }}-E_{\text {sech }}\right|}{E_{\text {total }}} \times 100 \% $
(13) 式中, Esech=2PpeakTFWHM/1.763代表双曲正割脉冲能量, Ppeak为双曲正割脉冲的峰值功率,TFWHM表示双曲正割脉冲的半峰全宽,Etotal是一个常量,代表传输脉冲的总能量。经计算,锁模激光器输出的自相似脉冲的基座能量比例为15.0354%, 当脉冲的总能量越接近双曲正割脉冲时,Pe越小,说明输出脉冲得到基座能量越小。
从图 4b可以看出, 脉冲的峰值功率随着传输距离的增加而增大,基座能量比例曲线存在两个极小值,分别是0.0431%(在6.324 m处)和0.0543%(在8.160 m处),针对这两个点(A和B)来讨论脉冲后续在SMF2的传输特性。
将这两个脉冲传输到足够长的SMF2中,脉冲的时域图如图 5所示。红色虚线和黑色实线区域分别为脉冲在SMF1和SMF2中的演化过程,两组脉冲在SMF2的峰值功率呈现周期性地增大后减小的趋势,关于两组脉冲在SMF2的峰值功率最高处的数据如表 1所示。相比于极小值点A, 极小值点B达到其最大峰值功率所需的SMF2更短,脉宽更窄,基座能量比例更小,脉冲的整体质量更高,更有利于后期输出压缩脉冲。
表 1 极小值点A和B对应的光纤长度和输出脉冲的参数
Table 1. Optical fiber length and output pulse parameters corresponding to minimum points A and B
minimum point A B length of SMF1/m 6.324 8.160 peak power P1/W 75.2949 413.1303 FWHM/fs 346.9 126.3 Pe/% 0.0431 0.0543 length of SMF2/m 12.000 2.160 peak power P2/W 981.63 1969.74 FWHM/fs 30.9 22.5 Pe/% 41.76 15.0896 -
为了得到质量更高的压缩脉冲,在保持耦合器功分比不变以及上下臂的SMF1长度保持一致的前提下,调整SMF2的长度,观察在特定光纤长度内脉冲参数的变化。已知脉冲在SMF2的传输过程如图 5b所示,脉冲演化过程中最大峰值功率出现在z=2.160 m处, 选取1.296 m~3.024 m区间的SMF2长度来分析SMF2对后期输出压缩脉冲的影响。表 2中为在SMF2的区间均等分所选取的数据点,便于分析各个点对应脉冲的差异。
表 2 数据点所对应的SMF2长度
Table 2. SMF2 length corresponding to the data point
data point 1 2 3 4 5 6 7 8 9 length of SMF2/m 1.296 1.512 1.728 1.944 2.160 2.376 2.592 2.808 3.024 图 6a显示了对应SMF2长度u1′和u3的峰值功率大小。通过对比SMF2的输出脉冲u1′和最终获得的压缩脉冲u3,两个脉冲峰值功率曲线都是先增大后减小的走势,可以发现, 在2.160 m(即第6个数据点)前,压缩脉冲的峰值功率要高于上臂的输出脉冲,2.160 m之后压缩脉冲的峰值功率降低至200 W以下。如图 6b所示,从u1′的演化趋势可以看出,在随着光纤长度的增加,脉冲先被压缩后展宽,与峰值功率曲线的演化趋势相反,在2.592 m前,压缩脉冲的脉宽要大于u1′,这是因为下臂的输出脉冲脉宽较大,两个脉冲相互干涉后,脉冲的基座部分通过干涉相消的方法去除,脉冲的中央部分干涉相长。图 6c为基座能量比例曲线。表明了压缩脉冲的基座明显比未干涉前要小,且两条曲线的差值从较大变为较小,在2.160 m前,上臂的脉冲峰值功率增大的同时基座能量比例减小,前6个数据点的压缩脉冲能够保持在10.5%以下。
用压缩品质因子Q来综合评价压缩脉冲质量[27]:
$ Q=\frac{R}{F}=\frac{P_{\text {out }} T_0}{P_{\text {in }} T_{\text {out }}} $
(14) 式中,品质因子R是输出脉冲与输入脉冲的峰值功率比,Pin为上臂输入脉冲的峰值功率,Pout为压缩脉冲的峰值功率,压缩因子F是输出脉冲与输入脉冲半峰全宽的比值,Tout为压缩脉冲的半峰全宽,T0为上臂输入脉冲的半峰全宽,Q越大,压缩效果越好。
从图 7中可以看出, 脉冲的压缩品质因子曲线存在极大值,位于光纤长度2.16 m处,因此,该长度的SMF2经过耦合器干涉后输出的脉冲压缩效果最好。
-
输入脉冲在单模光纤中传输存在脉宽小于1 ps的情况,即使β2≠0,也必须包含β3项,这是因为参量Δω/ω0(Δω为谱宽,ω0为中心频率)没有足够小到可把β3项简而化之。脉冲在SMF中传输特性用包含3阶色散的非线性薛定谔方程表示:
$ \begin{gathered} \frac{\partial u(z, T)}{\partial z}=-\frac{\mathrm{i} \beta_2}{2} \frac{\partial^2 u(z, T)}{\partial T^2}+ \\ \frac{\beta_3}{6} \frac{\partial^3 u(z, T)}{\partial T^3}+\mathrm{i} \gamma|u(z, T)|^2 u(z, T) \end{gathered} $
(15) 式中,β3为3阶色散系数,等号的右侧分别代表着2阶色散、3阶色散和非线性效应。取3阶色散系数为0.01 ps3/km,图 8为考虑3阶色散的上臂脉冲演化的时域波形图和上臂输出脉冲与压缩脉冲的对比波形图。
图 8 a—上臂的输出时域图 b—输出脉冲与压缩脉冲的波形图
Figure 8. a—output time-domain diagram of the upper arm b—waveform diagram of output pulse and compression pulse
从图 8a可以看出,脉冲的峰值功率相比于只考虑2阶色散有所提升,但是3阶色散效应会使脉冲的波形不再对称,且产生带有振荡结构的拖尾。从基座上分析,SMF1段输出脉冲的基座能量比例仅为0.1104%,在SMF2中压缩,基座能量比例增长至2.9114%,通过耦合器的干涉作用,上下臂的脉冲相互干涉,最后输出的压缩脉冲的基座能量比例高达至17.653%,干涉仪没有达到减小基座的要求。这是因为脉冲的3阶色散给脉冲带来位移,上下臂的输出脉冲中心不在同一直线上,从压缩脉冲的波形图可以得知脉冲的前沿变缓,后沿的基座得到有效的消除。
以0.001 ps3/km为步长, 将3阶色散系数从0取到0.01 ps3/km, 对应的压缩脉冲的峰值功率、半峰全宽以及基座能量比例曲线如图 9所示。3阶色散系数的增大对压缩脉冲的峰值功率增长有一定程度的抑制。半峰全宽曲线在3阶色散系数较小的时候呈现较缓的增长,当3阶色散系数大于0.004 ps3/km时,脉宽不断减小。基座能量比例曲线的最低点出现在0.001 ps3/km处,随着3阶色散系数的增大,曲线出现小范围的减小,但基座能量比例均大于11%, 为了减少3阶色散对压缩脉冲的影响,3阶色散系数应小于0.001 ps3/km。
基于Mach-Zehnder干涉仪的自相似激光器压缩系统设计
Design of self-similar laser compression system based on Mach-Zehnder interferometer
-
摘要: 为了研究自相似脉冲在Mach-Zehnder干涉仪的压缩特性, 采用非线性薛定谔方程对自相似脉冲的演化和压缩进行了模拟, 分析了基于级联单模光纤的Mach-Zehnder干涉仪的光纤参数对脉冲压缩的影响。结果表明, 在不考虑高阶色散的情况下, 当上臂的两种单模光纤长度分别为8.16 m和2.16 m、下臂的单模光纤长度为8.16 m时, 获得半峰全宽为27.85 fs、峰值功率为1860.59 W、基座能量比例为10.241%的最佳压缩脉冲; 考虑高阶色散时, 脉冲在单模光纤中传输呈现出峰值功率增大、基座增大的现象, 且脉冲右移不利于输出基座较小的压缩脉冲; 当3阶色散系数小于0.001 ps3/km时, 利用Mach-Zehnder干涉仪来压缩能获得质量较好的飞秒脉冲。该研究结果对于自相似脉冲的压缩研究具有一定的参考价值。Abstract: In order to study the compression characteristics of the self-similar pulse in the Mach-Zehnder interferometer(MZI), the evolution and compression of the self-similar pulse were simulated by using the nonlinear Schrödinger equation. The influence of the fiber parameters on pulse compression of the MZI based on the cascade single-mode fiber was analyzed. The results show that without considering the high-order dispersion, when the length of the two single-mode fibers in the upper arm is 8.16 m and 2.16 m, respectively, and the length of the single-mode fiber in the lower arm is 8.16 m, the optimal compression pulse is obtained with the full width at half maximum of 27.85 fs, peak power of 1860.59 W and the pedestal energy ratio is 10.241%. When high-order dispersion is considered, it is found that pulse transmission in single-mode fiber presents the phenomenon of the peak power increase and the pedestal energy ratio increase, and the right shift of pulse is not favorable to output compression pulse with a small pedestal. When the third-order dispersion is less than 0.001 ps3/km, femtosecond pulses of good quality can be obtained by using MZI compression. The result of this study has certain reference values for self-similar pulse compression research.
-
表 1 极小值点A和B对应的光纤长度和输出脉冲的参数
Table 1. Optical fiber length and output pulse parameters corresponding to minimum points A and B
minimum point A B length of SMF1/m 6.324 8.160 peak power P1/W 75.2949 413.1303 FWHM/fs 346.9 126.3 Pe/% 0.0431 0.0543 length of SMF2/m 12.000 2.160 peak power P2/W 981.63 1969.74 FWHM/fs 30.9 22.5 Pe/% 41.76 15.0896 表 2 数据点所对应的SMF2长度
Table 2. SMF2 length corresponding to the data point
data point 1 2 3 4 5 6 7 8 9 length of SMF2/m 1.296 1.512 1.728 1.944 2.160 2.376 2.592 2.808 3.024 -
[1] SHAH L, FERMANN M. High-power ultrashort-pulse fiber amplifiers[J]. IEEE Journal of Selected Topics in Quantum Electronics, 2007, 13(3): 552-558. doi: 10.1109/JSTQE.2007.896096 [2] 王竞, 李建中, 温伟峰, 等. 利用自相关方法实现光脉冲时间延迟精确测量[J]. 中国光学, 2015, 8(2): 270-276. WANG J, LI J Zh, WEN W F, et al. Precisely measuring for optical pulse time delay using autocorrelation[J]. Chinese Optics, 2015, 8(2): 270-276(in Chinese). [3] 夏照远, 钱静, 王关德, 等. 超短脉冲激光加工金属玻璃研究进展[J]. 激光与光电子学进展, 2021, 58(15): 1516027. XIA Zh Y, QIAN J, WANG G D, et al. Research progress on ultrashort pulsed laser processing of metallic glasses[J]. Laser & Optoelectronics Progress, 2021, 58(15): 1516027(in Chinese). [4] 邹贵生, 林路禅, 肖宇, 等. 超快激光纳米连接及其在微纳器件制造中的应用[J]. 中国激光, 2021, 48(15): 1502001. ZOU G Sh, LIN L Ch, XIAO Y, et al. Ultrafast laser nanojoining and its applicantions in the manufacturing of micro-nano devices[J]. Chinese Journal of Lasers, 2021, 48(15): 1502001(in Chinese). [5] 刘宇翔, 张瑞康, 王欢, 等. 1.5-μm波段25-GHz重频亚皮秒脉冲输出半导体锁模激光器(特邀)[J]. 光子学报, 2022, 51(2): 0251211. LIU Y X, ZHANG R K, WANG H, et al. 25 GHz semiconductor mode-locked laser with subpicosecond pulse output in the 1.5 μm band(Invited)[J]. Acta Photonica Sinica, 2022, 51(2): 0251211(in Chin-ese). [6] SALTARELLI F, GRAUMANN I J, LANG L, et al. Power scaling of ultrafast oscillators: 350-W average-power sub-picosecond thin-disk laser[J]. Optics Express, 2019, 27(22): 31465. doi: 10.1364/OE.27.031465 [7] MULLER M, ALESHIRE C, KLENKE A, et al. 10.4 kW coherently combined ultrafast fiber laser[J]. Optics Letters, 2020, 45(11): 3083-3086. doi: 10.1364/OL.392843 [8] WANG W, WU H, LIU C, et al. Multigigawatt 50 fs Yb∶CALGO regenerative amplifier system with 11 W average power and mid-infrared generation[J]. Photonics Research, 2021, 9(8): 1439-1445. doi: 10.1364/PRJ.425149 [9] MENG Y, ZHANG S, JIN C, et al. Enhanced compression of femtosecond pulse in hollow-core photonic bandgap fibers[J]. Optics Communications, 2010, 283(11): 2411-2415. doi: 10.1016/j.optcom.2010.02.016 [10] LI Q, HUANG H. Effective pulse compression in dispersion decreasing and nonlinearity increasing fibers[J]. Optics Communications, 2015, 342: 36-43. doi: 10.1016/j.optcom.2014.11.086 [11] YE F, HUANG J, GANDHI M S A, et al. Nearly self-similar pulse compression of high-repetition-rate pulse trains in tapered silicon waveguides[J]. Journal of Lightwave Technology, 2021, 39(14): 4717-4724. doi: 10.1109/JLT.2021.3077607 [12] YU T, LIU X, PRYAMIKOV A, et al. Femtosecond pulse compression with pedestal suppression in a sagnac interferometer constructed of anti-resonant hollow core fiber[J]. IEEE Photonics Journal, 2021, 13(2): 1-9. [13] 李沐霖, 张巧芬, 史圣达. 基于啁啾补偿技术的自相似脉冲压缩光纤设计[J]. 激光技术, 2021, 45(5): 566-570. LI M L, ZHANG Q F, SHI Sh D. Design of self-similar pulse compression fiber based on chirp compensation technology[J]. Laser Technology, 2021, 45(5): 566-570(in Chinese). [14] 赵羽, 刘永智, 赵德双, 等. 光纤激光器锁模技术研究进展[J]. 激光技术, 2009, 33(2): 162-165. ZHAO Y, LIU Y Zh, ZHAO D Sh, et al. Evolution of mode-locked technology of fiber lasers[J]. Laser Technology, 2009, 33(2): 162-165(in Chinese). [15] WANG F, LI Q. Cascaded single mode fibers for higher-order soliton compression at 2 μm[J]. Applied Optics, 2020, 59(17): E17-E22. doi: 10.1364/AO.385682 [16] 徐永钊, 张耿, 叶海, 等. 级联单模光纤中初始啁啾对高阶孤子脉冲压缩的影响[J]. 发光学报, 2016, 37(11): 1360-1366. XU Y Zh, ZHANG G, YE H, et al. Effect of initial frequency chirp on pulse compression of higher-order solitons in cascaded sing-mode fibers[J]. Chinese Journal of Luminescence, 2016, 37(11): 1360-1366(in Chinese). [17] 曹文华, 徐平, 刘颂豪, 等. 马赫-曾德尔型色散渐减光纤干涉仪的孤子效应脉冲压缩研究[J]. 光学学报, 2011, 31(4): 0419001. CAO W H, XU P, LIU S H, et al. Soliton-effect pulse compression in a dispersion-decreasing fiber-based Mach-Zehnder interferometer[J]. Acta Optica Sinica, 2011, 31(4): 0419001 (in Chinese). [18] FERMANN M E, KRUGLOV V I, THOMSEN B C, et al. Self-similar propagation and amplification of parabolic pulses in optical fibers[J]. Physical Review Letters, 2000, 84(26): 6010-6013. [19] ZHANG Q F, DENG Y H. Influence of gain coefficient on the self-similar pulses propagation in a dispersion-decreasing fiber[J]. Optik, 2016, 127(12): 5110-5114. [20] ZHANG Q F. Investigation on generalized analytical solution of similariton chirp in different tapered DDF and NIF[J]. Optical and Quantum Electronics, 2021, 53(1): 54. [21] 汪徐德, 周正, 李素文, 等. 掺镱光纤放大器中脉冲自相似演化特性分析[J]. 激光技术, 2012, 36(1): 8-12. WANG X D, ZHOU Zh, LI S W, et al. Self-similar pulse evolution in ytterbium doped fiber amplifiers[J]. Laser Technology, 2012, 36(1): 8-12(in Chinese). [22] DU Y, SHU X. Transformation from conventional dissipative solitons to amplifier similaritons in all-normal dispersion mode-locked fiber lasers[J]. IEEE Photonics Journal, 2018, 10(1): 1-11. [23] HIROOKA T, NAKAZAWA M. Parabolic pulse generation by use of a dispersion-decreasing fiber with normal group-velocity dispersion[J]. Optics Letters, 2004, 29(5): 498-500. [24] DORAN N J, WOOD D. Nonlinear-optical loop mirror[J]. Optics Letters, 1988, 13(1): 53-58. [25] ZHANG Q F, LI H, WU L M, et al. Research on evolution region of self-similar pulses in a dispersion-decreasing fiber[J]. Optical and Quantum Electronics, 2019, 51(6): 190. [26] CAO W H, WAI P K. Picosecond soliton transmission by use of concatenated gain-distributed nonlinear amplifying fiber loop mirrors[J]. Applied Optics, 2005, 44(35): 7611-7620. [27] 李沐霖, 张巧芬, 史圣达, 等. 基于啁啾补偿的自相似脉冲压缩光栅对的设计[J]. 光通信研究, 2021(4): 56-60. LI M L, ZHANG Q F, SHI Sh D, et al. Design of self-similar pulse compression grating pair based on chirp compensation[J]. Study on Optical Communications, 2021(4): 56-60(in Chinese).