## Abstract

Ultrafast Tm-doped fibre lasers have been actively studied for the last decade due to their potential applications in precise mid-IR spectroscopy, LIDARs, material processing and more. The majority of research papers is devoted to the comparison between a numerical modelling and experimental results; however, little attention is being paid to the comprehensive description of the mathematical models and parameters of the active and passive components forming cavities of Tm-doped all-fibre lasers. Thus, here we report a numerical model of a stretched-pulsed Tm-doped fibre laser with hybrid mode-locking and compare it with experimental results. The key feature of the developed numerical model is employment of the experimentally measured dispersion coefficients and optimisation of some model parameters, such as the bandwidth of the spectral filter spectral filtering and the saturation power of the active fibre, for a conformity with the experiment. The developed laser emits 331.7 fs pulses with a 23.8 MHz repetition rate, 6 mW of average power, 0.25 nJ of pulse energy, and a 21.66 nm spectral bandwidth at a peak wavelength of 1899.5 nm. The numerical model characteristics coincide with experimentally achieved spectral width, pulse duration, and average power with inaccuracy of 4.7%, 5.4%, and 22.9%, respectively. Moreover, in the discussion of the work the main possible reasons influencing this inaccuracy are highlighted. Elimination of those factors might allow to increase accuracy even more. We show that numerical model has a good agreement with the experiment and can be used for development of ultrafast Tm-doped fibre laser systems.

## Introduction

Mode-locked lasers based on Tm-doped fibres have been intensively studied in the last decade due to a high demand on sources emitting ultrashort pulses at the wavelengths of around \(2\,\upmu {\text {m}}\) for a wide variety of applications, such as vibrational spectroscopy for medicine, environmental sensing and supercontinuum generation^{1,2}. The short pulse duration (less than 150 fs), the sufficient peak power, and the stable generation regime are the essential parameters of a laser for the broadband and coherent supercontinuum generation in highly nonlinear fibres^{3}.

The diversity of mode-locking mechanisms with high performance in the desired spectral region around \(2\,\upmu {\text {m}}\)—such as nonlinear polarisation evolution (NPE)^{4,5,6}, nonlinear loop mirror^{7,8}, single wall carbon nanotubes (SWCNT)^{9,10,11}, \({\text {Bi}}_2 {\text {Te}}_3\)^{12}, SESAM^{13,14,15}, \({\text {MoS}}_2\)^{16}, graphene^{17}—allows to achieve ultrashort pulses with required duration, power level along with high stability and regime repeatability. A huge variety of generation regimes were achieved in Tm-doped fibre systems, such as: soliton^{7,10,11,15,17,18,19,20,21}, stretched pulse^{5,6,9,10,15,22}, dissipative soliton^{4,8,9,10,13,15,23,24}, noise-like pulse^{14,23,25} and bound solitons^{12,16,26}. Soliton Tm-doped fibre lasers have relatively narrow spectral width and, thus, long pulse duration compared with stretched-pulse lasers. Solitons with relatively short pulse duration (190 fs) were achieved in the Tm-doped fibre laser with small absolute value of group delay dispersion (GDD)^{27}; however, the achieved pulse energy is only 20 pJ. Use of fibres with normal group velocity dispersion (GVD) (e.g. germanosilicate fibres with small core diameter (less than \(5\,\upmu {\text {m}}\) ) and high concentration of germanium (more than 20%)) or bulk diffraction grating based compressors allow to improve the performance of Tm-doped fibre systems at \(2\,\upmu {\text {m}}\) via dispersion management^{28}. The dispersion management allows to achieve pulse energies of more than 2 nJ^{22,29}.

The numerical modelling has always been a great tool for the design of mode-locked lasers and the pulse evolution analysis. Various numerical models based on different modifications of the nonlinear Schrödinger equation for simulation of intracavity pulse propagation were proposed^{30,31,32,33} and implemented^{34,35,36,37} to study ultrafast lasers during last three decades. The majority of these papers were focused on studying of Er- and Yb-doped fibre lasers because of the current relative maturity of these laser sources. However, the numerical modelling of Tm-doped ultrafast laser sources is of the current interest and was presented in^{20,21,29,38,39,40,41}. Despite the comprehensive analysis of generation regimes and the comparison between the numerical modelling and the experimental results, little attention is being paid to the detailed description of the used mathematical model and the parameters of active and passive components of cavity (GVDs of fibres, nonlinear coefficients of fibres, spectral characteristics of the laser cavity components, parameters of the active fibres, parameters of SAs, etc.).

Thus, in this work we present the numerical model of the stretched-pulse Tm-doped all-fibre laser with hybrid mode-locking that is verified by the experimental data. The numerical simulation based on the nonlinear Schrödinger equation solved by split-step Fourier method is used. Moreover, observation of the intracavity pulse evolution is conducted for the analysis of the regime performance. The stretched-pulse generation is achieved by obtaining slightly anomalous GDD using normal GVD fibres with high the germanium oxide concentration (30%) and the high numerical aperture for both dispersion management and total cavity nonlinearity magnification. The SWCNT saturable absorber and the NPE technique are chosen as the mode-locking mechanisms. Based on our experience, the use of a polariser in a laser cavity providing the NPE mechanism allows to prevent the generation of a continuous wave (CW) component along with ultrashort pulses, leading to a more efficient subsequent amplification of the pulses. The parameters of the spectral filter and active fibre are varied to correspond to the experiment results. The discrepancy of the pulse duration and the width of the spectrum obtained in the model and the experiment is 5.4% and 4.7%, respectively. The error of the modelled value of the output power is 22.9%. A potential reason for these discrepancies are discussed in the paper.

The remainder of this paper is organised as follows: Section 2 describes the laser system and gives a description of the theoretical model along with the parameters of the cavity components used in the numerical simulation. Section 3 describes the results achieved both numerically and experimentally and shows the dynamics of intracavity pulse evolution. In Section 4 the reasons for discrepancies between experimental and numerical results are discussed, and possible solutions to avoid errors are suggested.

## Laser setup and numerical model

Firstly, we describe the setup of the developed ultrafast Tm-doped all-fibre laser. Figure 1 shows a scheme of the thulium-doped all-fibre ring laser with hybrid mode-locking. The radiation of a CW erbium/ytterbium-co-doped fibre laser operating at a wavelength of 1550 nm is coupled with wavelength-division multiplexer (WDM) and comes through a segment of the step-index (\(\Delta {\text {n}}=0.012\), core diameter \({{\text {d}}}_{\mathrm{c}}= 10\,\upmu {\text {m}}\)) thulium-doped aluminosilicate (0.8 wt% thulium, 3.6 wt% aluminum) glass fibre with anomalous GVD of \(-\,70.56\,{{\text {ps}}}^{2}/{\text {km}}\) at 1900 nm, where it is dramatically absorbed (more than 98%). Small signal gain is shown in Fig. 2a. The form of the small signal gain was taken from the consideration of its similarity with emission cross section^{42}.

An isolator-polariser is used for unidirectional generation and in combination with two polarisation controllers (PCs) enables the NPE technique which works as the mode-locking mechanism together with the saturable absorber based on the polymer film with SWCNT located between two FC/APC connectors. An output coupler removes 81% of the total power from the cavity and has 5% of total insertion loss. An additional isolator on the output is used to prevent the beam to reflect back to the cavity. An FC/APC connector is used for convenient connection with measurement instruments. The additional fibre with the high concentration of germanium (30% germanium oxide, core diameter \({{\text {d}}}_{\mathrm{c}}=2.2\,\upmu {\text {m}}\)) (Hi–Ge) with a GVD value of \(108\,{{\text {ps}}}^{2}/{\text {km}}\) and nonlinear coefficient of 3.9 1/W/km is used for dispersion management and makes GDD slightly anomalous. It further helps to achieve the non-solitonic behaviour of the pulse. The dependencies of the cavity fibres dispersion parameters on the wavelength are depicted in Fig. 2b. The dispersion of the fibres was measured using the setup described in^{43}. In this section we present a general laser scheme, all exact fibre lengths are presented in the next section. The cavity GDD value is \({-}1.74 \times 10^{-3}\,{\hbox {ps}}^{2}\). The numerical model is based on the modified nonlinear Schrödinger equation and can be represented as a consecutive propagation through the cavity elements. The split-step Fourier method is used to solve the modified nonlinear Schrödinger equation for passive and active fibres^{31}:

Here *A*(*z*, *T*) is a slowly varying complex envelope of the pulse propagating through fibre with length z and time \(T = t-\beta _1z\), which represents the moving frame in real time *t* with an envelope group velocity \(\beta _1^{-1}\). Moreover, the amplitude is normalised such that \(|A(z,T)|^2\) gives the instantaneous power in watts. \(\beta _k\) are the parameters of \(k^{th}\) order of dispersion, \(\gamma\) is the nonlinear coefficient, \(\alpha\) is the attenuation constant, and *g* is the gain which spectral dependence is included in the frequency domain through the following parameter^{34}:

where \(g_0\) is the small signal gain which is depicted in Fig. 2a, \(E_{sat}\) is the saturation energy, which is defined through the saturation power and the repetition frequency as \(E_{sat} = P_{sat}/f_{rep}\), and the total pulse energy \(E_0\), defined as an integral of a field envelope over the time \(E_0 = \int _{ }{ }|A|^2dT\). All the parameters of the passive and active fibres used in the numerical model, including dispersion coefficients at the wavelength of 1900 nm, the loss or gain characteristics, are shown in Table 1. The nonlinear coefficients were calculated using the formula \(\gamma = \omega _0 n_2/(c A_{eff})\)^{31}, where \(\omega _0\) is the central frequency, \(n_2\) is the nonlinear refractive index and is used from the paper^{44}, *c* is the speed of light, and \(A_{eff}\) is the effective mode area and is calculated from the fibre V parameter^{45,46}.

To estimate the transmission functions of the saturable absorbers (SA) in the time domain the following equation was used:

where \({q_{sat}}\) is a saturated losses of a SA, and \({q(T,|A(T)|^2)}\) is a function of saturation of a SA. For the SWCNT SA the following differential equation for a slow SA^{47} is solved with the \(4^{th}\)-order Runge-Kutta method^{48}:

Here \({q_{uns}}\) and \({q_{sat}}\) are the levels of unsaturated and saturated losses taken from the insert of Fig. 9 in^{49}, \(\tau\) is the relaxation time of the SA taken from Fig. 7 in^{49} and \(P_{abs}\) is the saturation power of the SA calculated from the insert in Fig. 9 in^{49}. The similar characteristics of carbon nanotubes with polymer film base are reported in^{50,51}. The unsaturated loss here are referred to the maximum total losses of SA in the case of the low incident signal power. In the case of the NPE SA with relatively small relaxation time the equation (5) can be simplified to the equation^{52}:

Moreover, to introduce the spectral selectivity of the cavity components the spectral filter with Gaussian shape with the spectral bandwidth \(\Omega\) and spectral shift of the filter \(\omega _0\) is included in the model after NPE SA:

The coupler is described by the simple transmission function \(Q=(1-R_{out}-\alpha )\), where \(R_{out}\) is the output coupling efficiency and \(\alpha\) is the level of the internal loss of the coupler equal to 5%. The parameters of the passive components used in the model are depicted in Table 2. The achieved splice loss between fibres with different mode-fields (e.g. the splice between the Hi–Ge fibre and both Tm-doped and SMF fibres) are experimentally measured and are approximately equal to 0.63 dB.

## Experimental and numerical results

The stable mode-locking regime of the experimental setup is achieved by adjusting PCs at the incident pump power of 350 mW with an average output power of 6 mW. The measured characteristics are shown in Fig. 3. The achieved spectrum is shown in Fig. 3a (orange curve), the centre wavelength is 1899.5 nm, the spectral full width at half maximum (FWHM) is 21.66 nm. There are small side wings that are present in both time and spectral domains, and the possible reasons for their appearance are discussed in the last section. The experimentally achieved autocorrelation trace is depicted in Fig. 3b (orange curve) and has Gaussian shape with FWHM of 469 fs, that correspond to 331.7 fs FWHM of pulse. The pulse energy is equal to 0.25 nJ. However, as the time-bandwidth product (TBP) for bandwidth-limited Gaussian pulse is \(\approx \,0.44\)^{53,54} and the TBP for achieved pulse is equal to 0.595, the pulse has a chirp. Figure 3c shows the radio frequency spectrum, repetition rate of the pulses is 23.835 MHz and the signal-to-noise ratio is 60 dB. The insert in Fig. 3c shows the captured oscilloscope pulse trace with the amplitude modulation with a value of ± 2.85% relative to the average value.

Afterwords, the parameters of the spectral filter and the saturation power of the active fibre in the numerical model are manually varied that the spectrum, the average power, and the autocorrelation trace at the output of the FC/APC connector coincide most accurately with the measured values. The input field of the model is represented as the one-photon-per-mode noise which is amplified in the active fibre. The filter spectral bandwidth (\(\Omega\)) and the filter centre wavelength appears to be 5.63 rad × THz (10.8 nm), and 1901 nm, respectively, and the active fibre saturation power is 0.93 mW. In this particular model the total gain varies through the saturation power and, thus, defines the average power^{34}. A stable generation regime with a time window size of 20 ps and the number of the grid points of \({2^{11}}\) is achieved. It should be noted that the temporal and spectral resolutions of this numerical model is 10 fs and 50 GHz, respectively, however, during the analysis of the output power distributions in both spectral and temporal domain, they were interpolated with an increased resolution (in 100 times) leading to the decreased temporal and frequency grid step to 0.1 fs and 0.05 GHz. The increase of grid points to \({2^{13}}\) in time domain with the same bandwidth of the time window led to a relative change of pulse and spectral width less than 0.04%, while the computation time was increased by 3.7 times. The results of pulse generation in the designed laser scheme is presented in Fig. 4. This figure shows the spectral (a) and temporal (b) distributions of the intensity and the average power (c) at the laser output.

The convergence criterion is described in^{37}, and in our case the relative change \(\epsilon\) should be less \(1.2\cdot {10^{-6}}\) for at least 300 roundtrips. Thus, the stable pulse generation starts from the 430th roundtrip. The temporal shift of the pulse in the time window corresponds to the displacement of the central wavelength from 1900 nm (centre of the frequency window), and misalignment between the speed of the pulse and speed of the time frame due to nonzero dispersion parameter^{31}. The fibre lengths before measurement devices are also included in the model to directly compare the autocorrelation traces and spectra on the output of the experimental setup and theoretical model. The total length of SMF-28 at the laser output is 2.3 m, and contains an optical isolator with the total internal loss of 0.45 dB. The modeled average laser output power in the stable generation regime is 4.63 mW after the output isolator. The numerical autocorrelation trace and spectrum match with a high precision (Fig. 3a, b), the achieved pulse duration is 349.61 fs, the achieved spectrum FWHM is 22.67 nm.

Figure 5 shows the evolution of the three main pulse parameters as functions of the position in the cavity—the pulse duration, the spectral width, and the average power—and illustrates the principle characteristics of the laser operation. At the top of the figure there is a schematic of the laser which consists of the Tm-doped fibre (\({\text {Tm}}^{3+}\)), the nonlinear fibre with the high concentration of germanium (Hi–Ge), the standard telecommunication SMF-28 fibre, the optical isolator-polariser which enables the nonlinear polarisation evolution technique as the saturable absorber (NPE), the filter (F), the output coupler (OC), and the SWCNT as the saturable absorber (SA). Here we assume that the filter depicts spectral dependencies of all cavity components except fibre gain whose spectral dependency is included in the gain function; however, it is placed just after the NPE because this technique introduces the most of spectral selectivity in the cavity^{55,56,57}. The length of the finite elements (optical fibres) are shown in the horizontal axis of Fig. 5. The background colour depends on the pulse chirp: the red one refers to the negative chirp sign, and the green one refers to the positive chirp sign.

In the first 30 cm of fibre, the pulse undergoes the amplification, the spectral and temporal narrowing along with the negative chirp obtainment in the active fibre. After that the chirp of the pulse changes the sign. The temporal and spectral distributions at this point are shown in Fig. 6a, e, respectively. The temporal profiles of the pulse at all points are approximated using a Gaussian function. After that the pulse spectrum is narrowed along with temporal broadening until the end of the fibre. Temporal and spectral distributions at the end of the active fibre are shown in Fig. 6b, f, respectively. Afterwords, the pulse goes through Hi–Ge fibre. The length of this fibre is large enough to change the sign of the pulse chirp back to normal and enable the mechanism of spectral and temporal broadening together. The temporal and spectral distributions at the second zero-chirp point are shown in Fig. 6c, g. Then the pulse passes through a cascade of two SAs, filter and OC connected with each other by SMF-28 fibre. The additional cavity point is shown in Fig. 6: pulse at the output of the laser (d,h). On this basis, pulse dynamics proves the assumption about the stretched-pulse generation regime^{35}, as there are the pulse duration has 2 local maximums and local minimums, and the chirp changes the sign 2 times in the cavity.

The maximum spectral width of the generation regime is mostly determined by the interplay between the nonlinear processes and the normal dispersion that occurs in the Hi–Ge fibre. However, the optical filter partially limits the pulse spectrum. Even though the SAs do not determine pulse duration, the minimal pulse duration is defined by evolution in both active and passive fibres and the interplay between the pulse and spectral broadening in the Hi–Ge fibre. Moreover, the NPE technique contributes little in the pulse shaping, and only includes high spectral selectivity into the cavity, as the influence of the filter changes both spectral and temporal widths dramatically. The pulse compression factor, which is the ratio between the maximum and minimum pulse duration in the cavity, is approximately 10, thus the main pulse dynamics mechanism is the pulse breathing^{58} that arises due to the dispersion map and the nonlinear pulse evolution. The nonlinear phase (B-integral) in the cavity is 1.69 \(\pi\), which is consistent with the value for the stretched-pulse generation regime^{59}.

The achieved experimental results and the direct comparison with the numerical modelling show that the developed numerical model is accurate enough to predict output characteristics of the generation regime of ultrafast Tm-doped fibre lasers. Moreover, the presented numerical model has proved to be a good tool for the analysis of generation regimes in Tm-doped laser, as it shows intracavity evolution of pulse characteristics e.g. both spectral and temporal width and shape, as well as chirp and power distributions.

## Discussion

In this section we want to briefly discuss discrepancies between the model and the experiment and the possible reasons for that behavior. Compared to the theoretical data the measured pulse autocorrelation has small side wings, and there are side peaks in the measured spectrum. In addition, the measured output power (6 mW) is more than the theoretical one (4.7 mW). Moreover, the pulse train has the small amplitude modulations that are not predicted in the model. Since the energy of the pulse turned out to be larger in the experiment than in the model it can be assumed that the appearance of side wings and peaks in the autocorrelation and in the spectrum is associated with the nonlinear deformation of the pulse in both time and spectral domains in the measurement line. Figure 7 shows the propagation of the laser pulse achieved in the numerical modelling in the extended measurement fibre in spectral (a) and temporal (b) domain.

The comparison of the spectrum at the laser output and at the compression point (2.6 m) is shown in Fig. 7c. The autocorrelation trace of the pulse at the compression point with a Gaussian approximation is shown in Fig. 7d. The spectrum undergoes the broadening at the compression point and the small side wings appear in the time domain. The inset of Fig. 7d shows that in comparison to the Gaussian approximation there are lateral wings in the autocorrelation trace of the pulse in the compression point. We assume that all the above disparity between the model and the experiment are related to the following assumptions accepted in the model. Firstly, we use an optimised Gaussian filter in the cavity, despite the fact that each component of the cavity introduces its own spectral filtering to the pulse. Secondly, we presented a simplified model of the active medium; however, the laser rate equations are more accurate. Thirdly, errors in parameters of the cavity components can affect discrepancies, including the dispersion coefficients and the nonlinear coefficients of all fibres, parameters of SWCNT, etc. Finally, we neglect the birefringence and the effect of stimulated Raman scattering in media. Despite these assumptions the pulse behavior in time and spectral domains in the model and in the experiment has very good agreement. Analysis of the influence of the above mentioned effects on the pulse dynamics inside and beyond the cavity is might be a possible solution for accuracy increase for the problems of modelling of more complex Tm-doped fibre systems.

## Conclusion

In this paper the numerical model of an all-fibre Tm-doped ultrafast laser is presented and experimentally verified. The generated pulses have the 331.7 fs duration and the 21.66 nm spectral bandwidth. The evolution of the pulse within the cavity was studied with the adapted numerical model and showed vast spectral and temporal breathing during the pulse propagation which corresponds to the stretched-pulse generation regime. The developed numerical model has good agreement with the experimental results and can be used as the robust and versatile instrument for the Tm-doped fibre laser design and the pulse generation regime analysis. Discrepancies between the model and the experiment were discussed. The developed laser setup is a first step towards the system of supercontinuum generation in the mid-IR region.

## References

- 1.
Rudy, C. W., Digonnet, M. J. & Byer, R. L. Advances in 2-\(\mu\)m Tm-doped mode-locked fiber lasers.

*Opt. Fiber Technol.***20**, 642–649 (2014). - 2.
Donodin, A.

*et al.*Supercontinuum Generation in a \({\text{As}}_2 {\text{S}}_3\)-Silica Nanospike Waveguide Pumped by Tm-Doped Fiber Laser. In*2019 Conference on Lasers and Electro-Optics Europe & European Quantum Electronics Conference (CLEO/Europe-EQEC)*, 1–1 (IEEE, 2019). - 3.
Genty, G., Coen, S. & Dudley, J. M. Fiber supercontinuum sources (invited).

*J. Opt. Soc. Am. B***24**, 1771–1785. https://doi.org/10.1364/JOSAB.24.001771 (2007). - 4.
Sobon, G., Sotor, J., Martynkien, T. & Abramski, K. M. Ultra-broadband dissipative soliton and noise-like pulse generation from a normal dispersion mode-locked Tm-doped all-fiber laser.

*Opt. Express***24**, 6156–6161 (2016). - 5.
Haxsen, F.

*et al.*Stretched-pulse operation of a thulium-doped fiber laser.*Opt. Express***16**, 20471–20476 (2008). - 6.
Zhao, D.

*et al.*Generation of stretched pulse in dispersion-managed Tm-doped fiber laser. In*14th National Conference on Laser Technology and Optoelectronics (LTO 2019)*, vol. 11170, 111701U (International Society for Optics and Photonics, 2019). - 7.
Michalska, M. & Swiderski, J. All-fiber thulium-doped mode-locked fiber laser and amplifier based on nonlinear fiber loop mirror.

*Opt. Laser Technol.***118**, 121–125 (2019). - 8.
Zhao, J.

*et al.*100 W dissipative soliton resonances from a thulium-doped double-clad all-fiber-format MOPA system.*Opt. Express***24**, 12072–12081 (2016). - 9.
Wang, Y.

*et al.*Generation of stretched pulses and dissipative solitons at 2 \(\mu\)m from an all-fiber mode-locked laser using carbon nanotube saturable absorbers.*Opt. Lett.***41**, 3864–3867 (2016). - 10.
Wang, J.

*et al.*152 fs nanotube-mode-locked thulium-doped all-fiber laser.*Sci. Rep.***6**, 28885 (2016). - 11.
Filatova, S. A.

*et al.*Spectral and temporal dynamics of ultrashort pulses in a holmium-doped fibre amplifier.*Quantum Electron.***49**, 1108–1111. https://doi.org/10.1070/qel17162 (2019). - 12.
Yin, K.

*et al.*Soliton mode-locked fiber laser based on topological insulator \({\text{Bi}}_2{\text{Te}}_3\) nanosheets at 2 \(\mu\)m.*Photonics Res.***3**, 72–76 (2015). - 13.
Huang, C.

*et al.*Developing high energy dissipative soliton fiber lasers at 2 micron.*Sci. Rep.***5**, 13680 (2015). - 14.
Mashiko, Y., Fujita, E. & Tokurakawa, M. Tunable noise-like pulse generation in mode-locked Tm fiber laser with a SESAM.

*Opt. Express***24**, 26515–26520 (2016). - 15.
Gumenyuk, R., Vartiainen, I., Tuovinen, H. & Okhotnikov, O. G. Dissipative dispersion-managed soliton \(2~\upmu \text{m }\) thulium/holmium fiber laser.

*Opt. Lett.***36**, 609–611. https://doi.org/10.1364/OL.36.000609 (2011). - 16.
Wang, P., Zhao, K., Gui, L., Xiao, X. & Yang, C. Self-organized structures of soliton molecules in \(2{\text{-}}\upmu {\text{m }}\) fiber laser based on \({\text{MoS}}_2\) Saturable Absorber.

*IEEE Photonics Technol. Lett.***30**, 1210–1213 (2018). - 17.
Sobon, G.

*et al.*All-polarization maintaining, graphene-based femtosecond Tm-doped all-fiber laser.*Opt. Express***23**, 9339–9346 (2015). - 18.
Nelson, L., Ippen, E. & Haus, H. Broadly tunable sub-500 fs pulses from an additive-pulse mode-locked thulium-doped fiber ring laser.

*Appl. Phys. Lett.***67**, 19–21 (1995). - 19.
Kieu, K. & Wise, F. Soliton thulium-doped fiber laser with carbon nanotube saturable absorber.

*IEEE Photonics Technol. Lett.***21**, 128–130 (2008). - 20.
Sotor, J.

*et al.*All-polarization-maintaining, stretched-pulse Tm-doped fiber laser, mode-locked by a graphene saturable absorber.*Opt. Lett.***42**, 1592–1595. https://doi.org/10.1364/OL.42.001592 (2017). - 21.
Wienke, A.

*et al.*Ultrafast, stretched-pulse thulium-doped fiber laser with a fiber-based dispersion management.*Opt. Lett.***37**, 2466–2468. https://doi.org/10.1364/OL.37.002466 (2012). - 22.
Chernysheva, M. A., Krylov, A. A., Kryukov, P. G. & Dianov, E. M. Nonlinear amplifying loop-mirror-based mode-locked thulium-doped fiber laser.

*IEEE Photonics Technol. Lett.***24**, 1254–1256 (2012). - 23.
Wang, Q.

*et al.*All-fiber ultrafast thulium-doped fiber ring laser with dissipative soliton and noise-like output in normal dispersion by single-wall carbon nanotubes.*Appl. Phys. Lett.***103**, 011103 (2013). - 24.
Gumenyuk, R., Vartiainen, I., Tuovinen, H. & Okhotnikov, O. G. Dissipative dispersion-managed soliton 2 \(\mu\)m thulium/holmium fiber laser.

*Opt. Lett.***36**, 609–611 (2011). - 25.
Voropaev, V.

*et al.*Generation of multi-solitons and noise-like pulses in a high-powered thulium-doped all-fiber ring oscillator.*Sci. Rep.***9**, 1–11 (2019). - 26.
Chernysheva, M.

*et al.*Double-wall carbon nanotube hybrid mode-locker in Tm-doped fibre laser: a novel mechanism for robust bound-state solitons generation.*Sci. Rep.***7**, 44314 (2017). - 27.
Sharp, R. C., Spock, D. E., Pan, N. & Elliot, J. 190-fs passively mode-locked thulium fiber laser with a low threshold.

*Opt. Lett.***21**, 881–883. https://doi.org/10.1364/OL.21.000881 (1996). - 28.
Tamura, K., Ippen, E. & Haus, H. Pulse dynamics in stretched-pulse fiber lasers.

*Appl. Phys. Lett.***67**, 158–160 (1995). - 29.
Haxsen, F., Wandt, D., Morgner, U., Neumann, J. & Kracht, D. Pulse characteristics of a passively mode-locked thulium fiber laser with positive and negative cavity dispersion.

*Opt. Express***18**, 18981–18988. https://doi.org/10.1364/OE.18.018981 (2010). - 30.
Dudley, J. M. & Taylor, J. R.

*Supercontinuum Generation in Optical Fibers*(Cambridge University Press, Cambridge, 2010). - 31.
Agrawal, G. P. Chapter 2—pulse propagation in fibers. In Agrawal, G. P. (ed.)

*Nonlinear Fiber Optics (Sixth Edition)*, 27–55, https://doi.org/10.1016/B978-0-12-817042-7.00009-9 (Academic Press, 2019). - 32.
Laegsgaard, J. Mode profile dispersion in the generalized nonlinear Schrödinger equation.

*Opt. Express***15**, 16110–16123 (2007). - 33.
Mamyshev, P. & Chernikov, S. V. Ultrashort-pulse propagation in optical fibers.

*Opt. Lett.***15**, 1076–1078 (1990). - 34.
Yarutkina, I., Shtyrina, O., Fedoruk, M. & Turitsyn, S. Numerical modeling of fiber lasers with long and ultra-long ring cavity.

*Opt. Express***21**, 12942–12950 (2013). - 35.
Wise, F. W., Chong, A. & Renninger, W. H. High-energy femtosecond fiber lasers based on pulse propagation at normal dispersion.

*Laser Photonics Rev.***2**, 58–73 (2008). - 36.
Oktem, B., Ülgüdür, C. & Ilday, F. Ö. Soliton-similariton fibre laser.

*Nat. Photonics***4**, 307–311 (2010). - 37.
Schreiber, T., Ortaç, B., Limpert, J. & Tünnermann, A. On the study of pulse evolution in ultra-short pulse mode-locked fiber lasers by numerical simulations.

*Opt. Express***15**, 8252–8262 (2007). - 38.
Xia, H.

*et al.*Characteristics of dissipative solitons in an all-fiber thulium-doped fiber ring laser.*Opt. Eng.***52**, 054201 (2013). - 39.
Li, H.

*et al.*Pulse-shaping mechanisms in passively mode-locked thulium-doped fiber lasers.*Opt. Express***23**, 6292–6303 (2015). - 40.
Lyu, Y.

*et al.*Theoretical comparison of NPR and hybrid mode-locked soliton thulium-doped fiber lasers.*IEEE Photonics J.***9**, 1–11. https://doi.org/10.1109/JPHOT.2017.2653860 (2017). - 41.
Wang, P.

*et al.*Dissipative rogue waves among noise-like pulses in a Tm fiber laser mode locked by a monolayer \({\text{MoS}}_2\) saturable absorber.*IEEE J. Sel. Top. Quantum Electron.***24**, 1–7. https://doi.org/10.1109/JSTQE.2017.2749972 (2018). - 42.
Agger, S. D. & Povlsen, J. H. Emission and absorption cross section of thulium doped silica fibers.

*Opt. Express***14**, 50–57 (2006). - 43.
Levchenko, A. E., Kurkov, A. S. & Semenov, S. L. Measurement of dispersion in optical fibres with a microstructure cladding.

*Quantum Electron.***35**, 835 (2005). - 44.
Yatsenko, Y. & Mavritsky, A. D-scan measurement of nonlinear refractive index in fibers heavily doped with \({\text{GeO}}_2\).

*Opt. Lett.***32**, 3257–3259 (2007). - 45.
Hussey, C. & Martinez, F. Approximate analytic forms for the propagation characteristics of single-mode optical fibres.

*Electron. Lett.***21**, 1103–1104 (1985). - 46.
Mortensen, N. A., Folkenberg, J. R., Nielsen, M. D. & Hansen, K. P. Modal cutoff and the V parameter in photonic crystal fibers.

*Opt. Lett.***28**, 1879–1881 (2003). - 47.
Paschotta, R. & Keller, U. Passive mode locking with slow saturable absorbers.

*Appl. Phys. B***73**, 653–662 (2001). - 48.
Hult, J. A fourth-order Runge–Kutta in the interaction picture method for simulating supercontinuum generation in optical fibers.

*J. Lightw. Technol.***25**, 3770–3775 (2007). - 49.
Chernysheva, M.

*et al.*Carbon nanotubes for ultrafast fibre lasers.*Nanophotonics***6**, 1–30 (2017). - 50.
Hasan, T.

*et al.*Nanotube–polymer composites for ultrafast photonics.*Adv. Mater.***21**, 3874–3899 (2009). - 51.
Chernysheva, M.

*et al.*High power Q-switched thulium doped fibre laser using carbon nanotube polymer composite saturable absorber.*Sci. Rep.***6**, 24220 (2016). - 52.
Kurtner, F. X., Der Au, J. A. & Keller, U. Mode-locking with slow and fast saturable absorbers—What’s the difference?.

*IEEE J. Sel. Top. Quantum Electron.***4**, 159–168 (1998). - 53.
Diels, J.-C.M., Fontaine, J. J., McMichael, I. C. & Simoni, F. Control and measurement of ultrashort pulse shapes (in amplitude and phase) with femtosecond accuracy.

*Appl. Opt.***24**, 1270–1282. https://doi.org/10.1364/AO.24.001270 (1985). - 54.
Lim, H., Ilday, F. & Wise, F. Generation of 2-nJ pulses from a femtosecond ytterbium fiber laser.

*Opt. Lett.***28**, 660–662 (2003). - 55.
Li, P.

*et al.*Tunable picosecond SESAM mode-locking ytterbium-doped double-clad LMA PCF oscillator.*Opt. Commun.***317**, 62–66 (2014). - 56.
Li, P.

*et al.*Subpicosecond SESAM and nonlinear polarization evolution hybrid mode-locking ytterbium-doped fiber oscillator.*Appl. Phys. B***118**, 561–566 (2015). - 57.
Chong, A., Renninger, W. H. & Wise, F. W. All-normal-dispersion femtosecond fiber laser with pulse energy above 20 nJ.

*Opt. Lett.***32**, 2408–2410 (2007). - 58.
Chong, A., Renninger, W. H. & Wise, F. W. Properties of normal-dispersion femtosecond fiber lasers.

*JOSA B***25**, 140–148 (2008). - 59.
Fu, W., Wright, L. G., Sidorenko, P., Backus, S. & Wise, F. W. Several new directions for ultrafast fiber lasers.

*Opt. Express***26**, 9432–9463 (2018).

## Acknowledgements

A.D. is grateful for being funded by the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement No 814276. The authors acknowledge the financial support of Russian Science Foundation (Project No. 17-79-20431). The authors are grateful to Krylov A.A. for assistance in developing a mode-locked laser, and to Obraztsova E.D. for providing single-walled carbon nanotubes for mode-locking in a laser.

## Author information

### Affiliations

### Contributions

V.L., A.D., V.V., M.T. designed the experiment. V.V., A.D., D.B., D.V. conducted the experiments. A.D., V.V. conducted the modelling. V.V., A.D., M.T., and V.L. analysed the results. A.D., V.V., D.B., D.V. wrote the manuscript. V.L. and A.D. supervised the project. M.L. designed the fibre with high concentration of germanium oxide. A.K. and D.V. designed the pump source for the laser. All authors reviewed the manuscript.

### Corresponding authors

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Donodin, A., Voropaev, V., Batov, D. *et al.* Numerical model of hybrid mode-locked Tm-doped all-fibre laser.
*Sci Rep* **10, **17396 (2020). https://doi.org/10.1038/s41598-020-74194-7

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.