-
基于RIEZNIK[15]和CERQUEIRA[16]等人的报道所提出的四波混频非线性效应的分析方法,将光纤激光的传输看作是多个纵模模式演变后的统计结果。利用四波混频非线性效应过程对此进行模拟,通过采用分步傅里叶方法对非线性薛定谔方程求解,结合掺镱光纤的实际参数条件,模拟掺镱光纤中非单模模式激光传输过程,得到相应的输出光谱,对掺镱光纤激光输出光谱的光谱特性进行分析,分析不同参数下四波混频非线性效应对掺镱光纤激光光谱展宽特性的影响。
-
四波混频效应是指光纤中不同波长的光波间相位匹配,即能量和波矢量都守恒的状态下,发生的非线性相互作用,并且在其它波长处产生边带,该效应的发生还与色散参数密切相关。四波混频效应源于光波光子湮灭新光子产生的过程。比如在3束光波入射到光纤的情况下,其光子分别以ω1,ω2和ω3表示,其中两束光波光子ω1和ω2湮灭,同时产生两个新光子ω3和ω4,新光子ω4可产生新的光波,而ω3则用于产生参量增益。当入射光波中ω1=ω2时,这种特殊情况称作简并四波混频。
四波混频过程可以用光场中非线性极化来描述[16]:
$ \begin{array}{c} P\left(\omega_4\right) \propto E_1 E_2 E_3^* \exp \left\{-\mathrm{i}\left[\beta\left(\omega_1\right)+\right.\right. \\ \left.\left.\beta\left(\omega_2\right)-\beta\left(\omega_3\right)\right] z\right\} \end{array} $
(1) 式中,β(ωj)为光纤在频率ωj(j=1, 2, 3, 4)处的传输常数,z为沿着光纤的位置坐标,E1、E2、E3表示频率ω1、ω2、ω3处的光场振幅,E3*是E3的共轭。当非线性极化的相位与光纤中传输的新光子光场相位一致时,将发生四波混频现象,这要求波数失配为0,即κ=0[17]:
$ \kappa=\Delta \kappa_{\mathrm{m}}+\Delta \kappa_{\mathrm{w}}+\Delta \kappa_{\mathrm{NL}}=0 $
(2) 式中,Δκm、Δκw和ΔκNL分别表示材料色散、波导色散和非线性效应引起的波数失配,(2)式满足能量守恒基本定律。
-
考虑ω1+ω2=ω3+ω4[16]的情况,角频率分别为ω1和ω2的抽运光入射到光纤中,当实现相位匹配时,主要发生两个现象: 一是角频率分别为ω3和ω4的斯托克斯光和反斯托克斯光将从噪声产生;二是如果弱信号光ω3与抽运光同时入射到光纤中,ω3会被放大,产生新光子ω4。将以上非线性过程用耦合非线性薛定谔方程进行描述[17]:
$ \begin{array}{c} \frac{\partial A_1}{\partial z}=\mathrm{i} \gamma\left[\left(\left|A_1\right|^2+2 \sum\limits_{m \neq 1}\left|A_m\right|^2\right) A_1+\right. \\ \left.2 A_2^* A_3 A_4 \exp (\mathrm{i} \Delta \beta z)\right] \end{array} $
(3) $ \begin{array}{c} \frac{\partial A_2}{\partial z}=\mathrm{i} \gamma\left[\left(\left|A_2\right|^2+2 \sum\limits_{m \neq 2}\left|A_m\right|^2\right) A_2+\right. \\ \left.2 A_1^* A_3 A_4 \exp (\mathrm{i} \Delta \beta z)\right] \end{array} $
(4) $ \begin{array}{c} \frac{\partial A_3}{\partial z}=\mathrm{i} \gamma\left[\left(\left|A_3\right|^2+2 \sum\limits_{m \neq 3}\left|A_m\right|^2\right) A_3+\right. \\ \left.2 A_4^* A_1 A_2 \exp (-\mathrm{i} \Delta \beta z)\right] \end{array} $
(5) $ \begin{array}{c} \frac{\partial A_4}{\partial z}=\mathrm{i} \gamma\left[\left(\left|A_4\right|^2+2 \sum\limits_{m \neq 4}\left|A_m\right|^2\right) A_4+\right. \\ \left.2 A_3^* A_1 A_2 \exp (-\mathrm{i} \Delta \beta z)\right] \end{array} $
(6) 式中,γ表示光纤的非线性系数,Am(m=1, 2, 3, 4)表示光波场的缓变振幅函数。以上方程是在暂时忽略走离效应、光纤色散以及损耗的影响下,将耦合非线性薛定谔方程组简化所得。而当光波在无损光纤中传播时,通常可以采用广义非线性薛定谔方程对其进行描述[17]:
$ \mathrm{i} \frac{\partial A}{\partial z}=-\frac{1}{2} \beta_2 \frac{\partial^2 A}{\partial t^2}+\frac{\mathrm{i}}{6} \beta_3 \frac{\partial^3 A}{\partial t^3}+\gamma|A|^2 A $
(7) 式中,A=A(z, t)表示光波场的缓变包络,β2和β3表示2阶和3阶色散系数,t表示脉冲的平均群速度,此处不考虑高阶色散系数如β4、β5等。
在分步傅里叶方法(split-step Fourier method, SSFM)中,将光纤按间隔长度h(h可沿光纤变化)分段,脉冲在光纤中传输,每一空间步长内对该非先线性薛定谔方程按照无色散情况进行求解,此处假设无色散(β3=β2=0)情况,直接在时域求解,可得非线性薛定谔方程的解:A(z, t)=A(0, t)exp[iγ|A(0, t)|2z][15],使用SSFM对其进行求解计算时,算法程序的精度和执行时间分别取决于空间步长h和时间步长dt的选择,RIEZNIK等人[15]对空间步长如何选择已经提出了有效可行的方法,本文的数值模拟中将沿用这种优化SSFM计算的方法,使数值结果更加准确。
-
掺镱光纤激光器分析的基本结构如图 1所示。图中,Ain处表示初始输入光谱,HR(high reflectivity)表示高反射率光纤光栅,OC(output couple)表示输出耦合光纤光栅,输入光谱设定为由两台抽运光源(表示为ω1,ω2)和光信噪比(optical signal-to-noise ratio, OSNR)为80 dB@12.5 GHz的高斯白噪声组成[16],输入到掺镱光纤中,进行光放大,模拟模型中采用纤芯和内包层直径分别为20 μm和400 μm的LMA-YDF-20/400-HP-M光纤作为增益光纤,L表示增益光纤长度。
在数值模拟中,除了增益光纤类型和抽运源的选择,光纤长度的选择对四波混频非线性效应的影响也很关键。大量的实验和理论研究结果表明,非线性过程的效率可以通过较大的非线性系数或者增加光纤长度实现,但在使用几百米乃至上千米长度的光纤时,非线性过程效率提高的同时,光纤中的随即双折射、群速度色散变化和相位失配现象也会随之增大。此外,在第1.1节中已经提到相位匹配是四波混频发生的前提,这也意味着相位失配ΔβL要尽量的小。虽然理论上增加光纤长度可以提高四波混频效应的效率,但随着光纤长度的增加,宽带噪声也会随之放大,对四波混频会产生更加严重的影响。而采用较短的光纤,此时几乎可以认为光纤构成了无损耗系统,可以显示多个频率间的周期性传输。另外,通过使用较短长度的光纤,也能够允许一定的相位失配。
由于LMA-YDF-20/400-HP-M掺镱光纤的非线性系数γ和2阶色散系数β2较小,结合实验中光纤长度选取情况,增益光纤长度设为28 m。参考光纤厂商数据,数值模拟中选用的参数如表 1所示。
表 1 输出光谱模拟所用的参数
Table 1. Parameters used for output spectrum simulation
symbol physical meaning of parameter value/unit L fiber length 28 m γ nonlinear Kerr coefficients 1.04 W-1·km-1 ROSNR optical signal-to-noise ratio 80 dB@12.5 GHz β2 second-order dispersion coefficient 0.01 ps2·m-1 S0 dispersion slope 0.2 ps·nm-2·km-1 α fiber loss 0.003 m-1 λ0 zero dispersion wavelength of fiber 1494.5 nm 模拟中依次改变抽运功率、抽运波长、激光中心波长、初始光场的注入模式数以及光纤长度等模拟参数,观察在不同参数变化时,四波混频非线性过程对掺镱光纤激光输出光谱展宽的影响。
-
根据上述四波混频非线性过程模型对该掺镱光纤激光的输出光谱进行了仿真,研究不同实验参数下四波混频对掺镱光纤激光光谱展宽的影响。在仿真过程中通过对抽运源的功率进行设置,得到不同抽运功率下掺镱光纤激光输出光谱,抽运功率从最低功率100 W,逐渐变化增加至2000 W。为更好地观察激光输出光谱展宽随抽运功率的变化,对光谱数据进行了高斯拟合数据处理,图 2a、图 2b和图 2c即880 nm、976 nm和1018 nm波长的抽运源抽运时得到1080 nm激光的输出光谱。从模拟结果可以观察到,随着抽运功率的增大,输出光谱线宽逐渐展宽,输出光谱的半峰全宽随抽运功率的变化曲线如图 2c所示。线宽随抽运功率增加而展宽,曲线尽管有一定起伏,但随着抽运功率的增大输出光谱线宽展宽趋势显著。
图 2 不同抽运功率下1080 nm掺镱光纤激光输出光谱
Figure 2. Output spectra of 1080 nm ytterbium-doped fiber lasers at different pump powers
根据BABIN等人[18]曾提出的将广义非线性薛定谔方程与光纤中传输模式振幅表达形式相结合,得到相应的激光传输的简化波动方程:
$ \begin{aligned} & {\left[\delta(\omega)+2 \alpha L+\delta_{\mathrm{NL}}\right] I(\omega)=2 g_{\mathrm{R}} L \bar{P} I(\omega)+} \\ & \frac{\delta_{\mathrm{NL}}}{I^2} \int I\left(\omega_1\right) I\left(\omega_2\right) \times I\left(\omega_1+\omega_2-\omega\right) \mathrm{d} \omega_1 \mathrm{~d} \omega_2 \end{aligned} $
(8) 式中,ω、ω1和ω2分别表示不同的光子频率;δ(ω)表示光纤布喇格光栅的有效损耗函数,该损耗曲线接近抛物线形式;α表示光纤损耗系数;P表示光纤中传输的平均抽运光功率;δNL表示四波混频带来的非线性相移损耗;I表示光纤中传输的斯托克斯光功率。右边第1项为喇曼增益项,gR表示喇曼增益系数,本文中仅考虑腔内增益g,所以(8)式可表示为:
$ \begin{aligned} & {\left[\delta(\omega)+2 \alpha L+\delta_{\mathrm{NL}}\right] I(\omega)=2 g L \bar{P} I(\omega)+} \\ & \frac{\delta_{\mathrm{NL}}}{I^2} \int I\left(\omega_1\right) I\left(\omega_2\right) \times I\left(\omega_1+\omega_2-\omega\right) \mathrm{d} \omega_1 \mathrm{~d} \omega_2 \end{aligned} $
(9) 式中,光纤中斯托克斯光传输功率谱I(ω)可写作如下形式[18]:
$ I(\omega ) = \frac{{2I}}{{\pi \mathit{\Gamma }\cosh (2\omega /\mathit{\Gamma })}} $
(10) 光谱宽度Γ表示为[18]:
$ \mathit{\Gamma } = \frac{2}{\pi }\sqrt {\frac{{2{\delta _{{\rm{NL}}}}}}{{{\delta _2}}}} $
(11) 式中,δ2表示光纤布喇格光栅的有效损耗系数,结合LAPOINTE等人[10]的报道,非线性相移损耗表示为:
$ \delta_{\mathrm{NL}}=\sqrt{\frac{2}{3}} \frac{\gamma P L}{1+\left(4 \beta L / 3 \delta_2\right)^2} $
(12) 在低色散短距离光纤激光系统中[4βL/(3δ2)]2可以忽略,则$\delta_{\mathrm{NL}}=\sqrt{2 / 3} \gamma P L$,P表示光纤中传输的斯托克斯平均光功率,$P=\int_0^L P(x) \mathrm{d} x / L$,P(x)表示光纤中的输入光功率,β表示光纤的正常色散系数,结合(11)式和(12)式可知,输出光谱宽度Γ与输入光功率的变化成正比,随着输入光功率的增加,掺镱光纤激光的输出光谱线宽发生展宽。所以在模拟中改变激光系统的输入抽运功率,会导致输出光谱线宽展宽,如图 2所示,模拟结果也与参考文献[19]中1150 nm掺镱光纤激光的实验研究的结果一致。
-
从图 2中的结果发现,不同抽运波长下掺镱光纤激光输出光谱展宽并不完全一致。短波长抽运如880 nm和976 nm,与1018 nm波长抽运进行比较,1018 nm和976 nm波长抽运下,光纤激光输出光谱线宽展宽明显快于880 nm波长抽运时的光谱线宽展宽,因此,通过使用不同波长的抽运源来研究抽运波长参数对掺镱光纤激光光谱展宽特性的影响。
抽运波长选用范围为808 nm~1018 nm,依然输出1080 nm激光,得到不同抽运波长下输出激光光谱展宽情况, 如图 3所示。从图 3a可以观察到,抽运功率的增加,不同抽运波长抽运时得到的输出光谱逐渐展宽,不同抽运波长相应的激光输出光谱展宽情况不同。当抽运功率小于600 W时,808 nm,880 nm,915 nm,938 nm,976 nm和1018 nm波长抽运所得激光输出光谱线宽随抽运功率变化的展宽较小,但展宽趋势一致。这是因为当增益光纤长度一定、抽运功率较小时,色散和非线性效应对输出激光线宽展宽贡献较少,需要更多的能量促进吸收转换光能量过程来实现光放大,色散和非线性效应随之累积放大,使得输出激光有明显的线宽展宽,所以600 W及以下的功率抽运时,各波长抽运下的输出光谱展宽现象并不显著。抽运功率增加至2000 W时,808 nm,915 nm,938 nm,976 nm和1018 nm波长抽运下的激光输出光谱线宽随之继续展宽,高功率时976 nm和1018 nm波长抽运下激光线宽有明显的展宽趋势,而808 nm波长抽运下,激光线宽展宽速度较其它抽运波长下激光线宽展宽速度则趋于平缓,这一现象在图 3b中可以更为直观地观察到。这是因为镱离子在808 nm处的吸收系数是6种抽运波长中最低的,如图 4所示,镱离子在其它波长处都有较高的吸收系数,甚至在915 nm和976 nm两处存在吸收峰,随着抽运功率的增加,输出激光线宽固然展宽,增益光纤掺杂粒子对抽运光的吸收利用效率也会使得不同波长光源抽运时输出激光线宽展宽现象不完全一致。
-
四波混频非线性过程中,抽运功率和抽运波长一定的情况下,不同输出激光波长的输出光谱展宽情况也会有所不同,为了探究四波混频非线性过程中不同中心波长对输出光谱展宽特性的影响,在仿真模拟中采用976 nm波长抽运产生不同中心波长的输出激光,分别实现1045 nm、1064 nm、1080 nm和1090 nm激光输出。图 5a和图 5b分别是100 W和1700 W功率抽运时的4种激光输出光谱。抽运功率增加时输出光谱强度有明显增强,这与之前的研究结果[19],即激光输出功率会随着抽运功率的提高而增大是一致的。得到的不同中心波长激光的输出光谱,如图 5c和图 5d所示。抽运功率增加,4种中心波长激光的输出光谱线宽都在展宽,对其进行数据处理,得到1045 nm,1064 nm,1080 nm和1090 nm激光的光谱线宽展宽曲线斜率分别为0.00099,0.00097,0.00090和0.00105,光谱展宽速度都较为接近,但是1080 nm和1090 nm的激光输出光谱展宽始终比1045 nm和1064 nm激光的光谱展宽现象明显,这一点在图 5d中可以更加清晰地看到,相同抽运功率下,实现长波长激光输出时,光谱线宽相较于较短波长更加展宽。
图 5 采用976 nm波长抽运源,不同中心波长激光的输出光谱线宽变化
Figure 5. Output spectra of lasers with different central wavelengths when pumped with a 976 nm pump source
光纤中平均抽运光功率P可由以下形式表达[18]:
$ \bar{P}=P_0 \frac{1-\exp \left(-2 \alpha_{\mathrm{p}} L-2 \frac{\lambda}{\lambda_{\mathrm{p}}} g L I\right)}{\left(\alpha_{\mathrm{p}}+\frac{\lambda}{\lambda_{\mathrm{p}}} g I\right) L} $
(13) 式中,P0表示输入抽运功率,αp表示光纤中抽运功率损耗系数,g为腔内增益系数,λ和λp分别表示中心波长和抽运波长,结合斯托克斯光谱表达式(10)式,非线性相移损耗还可以表示为[18]:
$ \delta_{\mathrm{NL}}=\sqrt{\frac{2}{3}} \frac{\gamma I L}{1+\left[4 \beta L /\left(3 \delta_2\right)\right]^2} $
(14) 根据(13)式,可以得到模型中固定抽运光功率下,光纤中斯托克斯光功率与激光中心波长成正比关系,即I∝λ,由(11)式可知,光谱宽度与非线性相移损耗的平方根成正比Γ∝δNL1/2,在(14)式中有δNL∝I,则光谱宽度与激光中心波长的平方根成比关系Γ∝λ1/2。所以长波长激光光谱线宽展宽始终较大,但由于光谱宽度是与波长的平方根成正比,中心波长变化对光谱宽度的影响被削弱,所以本文中模拟的4种中心波长的输出光谱线宽展宽趋势相近。
-
实际应用中,掺镱光纤激光功率提升一般采用主振荡功率放大结构,光谱展宽现象对放大器的性能影响不容忽视,不同的初始输入模式对光放大和功率放大过程的光谱展宽特性的影响也有差异。根据HUANG等人[20]的报道,在光谱中每一个纵模模式可以由一个复振幅表示,所有模式的演化可以用耦合复振幅方程来描述。本文中将模拟计算输入的初始光场假设是多个振幅相同的模式构成,共同注入掺镱光纤。在模拟中保持模间距不变,改变模式数目N,得到不同模式数下的掺镱光纤激光输出光谱,研究光纤中注入光场模式数变化对光纤激光光谱展宽造成的影响,结果如图 6所示。从图 6a中可知,相同抽运功率下,随着注入模式数的增加,输出光谱线宽逐渐展宽,在高功率抽运下,输出光谱展宽的现象非常明显,在图 6b中可以更为清晰地观察到。图 6b中展示了在相同模式数目时掺镱光纤激光输出光谱线宽的变化,光谱线宽随着抽运功率的增加快速展宽,这与图 2c中的结果也是一致的。这是因为在实际激光系统中,由于众多纵模模式的相位分布和强度不同,模式演化过程往往比较复杂。为了对非线性过程下光谱展宽进行更加直观的估计,文中模型假设注入的所有纵模具有相同的相位φ0和相等的振幅$\sqrt{P_{\mathrm{in}} / N}$,Pin表示光纤中输入的光功率,将非线性强度SNL表示为耦合复振幅方程,即前面(3)式~(6)式中的非线性项Ai*AkAl(i=2, 1, 4, 3;k=3, 3, 1, 1;l=4, 4, 2, 2)之和,可以得到非线性强度SNL与模式振幅成正比关系,进而与模式数成正比关系,即SNL∝N3/2。所以在图 6结果中,掺镱光纤激光光谱线宽随着模式数N的增加展宽。
四波混频对掺镱光纤激光光谱展宽特性的影响
Effect of four-wave mixing on spectrum broadening characteristics of ytterbium-doped fiber laser
-
摘要: 为了探究四波混频非线性过程对高功率掺镱光纤激光光谱展宽特性的影响,采用纤芯和内包层直径分别为20 μm和400 μm的掺镱光纤作为增益介质,基于四波混频非线性过程建立掺镱光纤激光光谱分析模型,分析了抽运功率、抽运波长、中心波长和注入模式数参数变化对掺镱光纤激光光谱展宽特性的影响,得到了1080 nm激光的输出光功率谱以及输出激光线宽随各参数变化的曲线。结果表明,四波混频非线性过程中,抽运功率的增大会导致掺镱光纤激光输出光谱线宽展宽;不同抽运波长抽运的输出光谱线宽展宽速率不同;不同中心波长的输出激光光谱展宽程度相近;掺镱光纤激光光谱线宽随着注入模式数增加逐渐展宽。该分析结果为掺镱光纤激光光谱展宽特性研究提供了理论支持,对高功率光纤激光的工程应用有一定的参考意义。Abstract: In order to investigate the influence of the nonlinear process of four-wave mixing on the spectral broadening characteristics of high-power ytterbium-doped fiber laser, ytterbium-doped fiber with core and inner cladding diameters of 20 μm and 400 μm, resperctively, was used as the gain medium. Based on the nonlinear process of four wave mixing, an ytterbium-doped fiber laser spectral analysis model was established. The effects of pump power, pump wavelength, center wavelength, and injection mode number on the spectral broadening characteristics of an ytterbium-doped fiber laser were analyzed. The output power spectrum of 1080 nm laser and the curve of output laser linewidth with the above parameters were obtained. The simulation results show that in the nonlinear process of four-wave mixing, the increase of pump power leads to the broadening of the output spectral linewidth of the ytterbium-doped fiber laser; the output spectral linewidth broadening rate of different pump wavelengths is different; the broadening degree of the output laser spectrum of the long band central wavelength is similar; the spectral linewidth of ytterbium-doped fiber laser gradually broadens with the increase of the number of injection modes. The simulation results provide theoretical support for the study of ytterbium-doped fiber laser spectral broadening characteristics, and have certain reference significance for the engineering application of high-power fiber laser.
-
表 1 输出光谱模拟所用的参数
Table 1. Parameters used for output spectrum simulation
symbol physical meaning of parameter value/unit L fiber length 28 m γ nonlinear Kerr coefficients 1.04 W-1·km-1 ROSNR optical signal-to-noise ratio 80 dB@12.5 GHz β2 second-order dispersion coefficient 0.01 ps2·m-1 S0 dispersion slope 0.2 ps·nm-2·km-1 α fiber loss 0.003 m-1 λ0 zero dispersion wavelength of fiber 1494.5 nm -
[1] 党文佳, 李哲, 卢娜, 等. 0.9~1.0 μm近红外连续光纤激光器的研究进展[J]. 中国光学, 2021, 14(2): 264-274. DANG W J, LI Zh, LU N, et al. Research progress of 0.9~1.0 μm near-infrared continuous-wave fiber lasers[J]. Chinese Optics, 2021, 14(2): 264-274(in Chinese). [2] 吴春婷, 常奥磊, 温雅, 等. 单掺Nd3+双波长全固态激光器研究进展[J]. 发光学报, 2020, 41(4): 414-428. WU Ch T, CHANG A L, WEN Y, et al. Research progress of Nd3+-doped dual-wavelength all-solid-state laser[J]. Chinese Journal of Luminescence, 2020, 41(4): 414-428(in Chinese). [3] 李虎, 郭子龙, 杨文婷, 等. 空芯光纤多模干涉型光纤液位传感技术研究[J]. 激光技术, 2022, 46(1): 120-124. LI H, GUO Z L, YANG W T, et al. Research of the liquid level sensing technology based on a hollow fiber multimode interference optical fiber[J]. Laser Technology, 2022, 46(1): 120-124(in Chinese). [4] 刘伟, 肖虎, 王小林, 等. 掺Yb光纤激光器输出光谱特性研究[J]. 中国激光, 2013, 40(9): 0902006. LIU W, XIAO H, WANG X L, et al. Study on output spectral characteristic of Yb-doped fiber lasers[J]. Chinese Journal of Lasers, 2013, 40(9): 0902006(in Chinese). [5] BABIN S A, CHURKIN D V, ISMAGULOV A E, et al. Turbulence-induced square-root broadening of the Raman fiber laser output spectrum[J]. Optics Letters, 2008, 33(6): 633-635. doi: 10.1364/OL.33.000633 [6] BEDNYAKOVA A E, GORBUNOV O A, POLITKO M O, et al. Generation dynamics of the narrowband Yb-doped fiber laser[J]. Optics Express, 2013, 21(7): 8177-8182. doi: 10.1364/OE.21.008177 [7] SMIRNOV S V, CHURKIN D V. Modeling of spectral and statistical properties of a random distributed feedback fiber laser[J]. Optics Express, 2013, 21(18): 21236-21241. doi: 10.1364/OE.21.021236 [8] GORBUNOV O A, SUGAVANAM S, CHURKIN D V. Revealing statistical properties of quasi-CW fibre lasers in bandwidth-limited measurements[J]. Optics Express, 2014, 22(23): 28071-28076. doi: 10.1364/OE.22.028071 [9] CHURKIN D V, KOLOKOLOV I V, PODIVILOV E V, et al. Wave kinetics of random fibre lasers[J]. Nature Communications, 2015, 6(1): 1-6. [10] LAPOINTE M A, PICHE M. Linewidth of high-power fiber lasers[C]// Photonics North 2009. Quebec, Canada: International Society for Optics and Photonics, 2009: 232-239. [11] TURITSYN S K, BEDNYAKOVA A E, FEDORUK M P, et al. Modeling of CW Yb-doped fiber lasers with highly nonlinear cavity dynamics[J]. Optics Express, 2011, 19(9): 8394-8405. doi: 10.1364/OE.19.008394 [12] KABLUKOV S I, ZLOBINA E A, PODIVILOV E V, et al. Output spectrum of Yb-doped fiber lasers[J]. Optics Letters, 2012, 37(13): 2508-2510. doi: 10.1364/OL.37.002508 [13] KABLUKOV S I, ZLOBINA E A, PODIVILOV E V, et al. Modeling and measurement of ytterbium fiber laser generation spectrum[C]// Laser Sources and Applications. Brussels, Belgium: International Society for Optics and Photonics, 2012: 36-44. [14] 肖虎. 掺镱光纤激光级联泵浦技术研究[D]. 长沙: 国防科学技术大学, 2012: 10. XIAO H. Study on tandem pumping technology of ytterbium-doped fiber lasers[D]. Changsha: National University of Defense Technology, 2012: 10(in Chinese). [15] RIEZNIK A A, TOLISANO T, CALLEGARI F A, et al. Uncertainty relation for the optimization of optical-fiber transmission systems simulations[J]. Optics Express, 2005, 13(10): 3822-3834. doi: 10.1364/OPEX.13.003822 [16] CERQUEIRA S A, BOGGIO J M C, RIEZNIK A A, et al. Highly efficient generation of broadband cascaded four-wave mixing products[J]. Optics Express, 2008, 16(4): 2816-2828. doi: 10.1364/OE.16.002816 [17] AGRAWAL G P. Nonlinear fiber optics, third[J]. Lecture Notes in Physics, 2001, 18(1): 195-211. [18] BABIN S A, CHURKIN DV, KABLUKOV S I, et al. Four-wave-mixing-induced turbulent spectral broadening in a long Raman fiber laser[J]. Journal of the Optical Society of America, 2007, B24(8): 1729-1738. [19] 谌鸿伟, 沈炎龙, 陶蒙蒙, 等. 1150 nm掺镱光纤激光器输出特性实验研究[J]. 光子学报, 2016, 45(10): 1014001. CHEN H W, SHEN Y L, TAO M M, et al. Experimental investigation of an 1150 nm Yb-doped fiber laser[J]. Acta Phtonica Sinica, 2016, 45(10): 1014001(in Chinese). [20] HUANG Zh H, LIANG X B, LI Ch Y, et al. Spectral broadening in high-power Yb-doped fiber lasers employing narrow-linewidth multilongitudinal-mode oscillators[J]. Applied Optics, 2016, 55(2): 297-302. doi: 10.1364/AO.55.000297