Open Access
Add to BibSonomy
Abstract
Pure-quartic soliton (PQS) fiber lasers provide a promising avenue for exploring novel soliton interaction dynamics and generating high-energy pulses. Here, we present the numerical observation of vector PQSs generation and the evolution dynamics in a mode-locked fiber laser, using the coupled Ginzburg-Landau equations. We investigate the buildup dynamics of vector PQSs in a mode-locked laser with birefringent fibers, passing through three stages: energy amplification, energy pulsation owing to the cross-phase modulation (XPM) effect, and finally stabilization. Depending on the strength of the cavity-birefringence, the evolution of PQSs in non-polarization-maintaining fibers reveals that both the elliptical-polarization vector PQSs and near-linear-polarization vector PQSs can be formed by the energy conservation and balance between the two orthogonal directions. Additionally, we observe the transition process from vector PQSs to scalar PQSs with higher cavity-birefringence, resulting from the failure compensation of the walk-off via the soliton trapping effect between the two orthogonal components. These results provide valuable insights into the ultrafast transient process of vector solitons and enhance the understanding of PQS generation in fiber lasers.
© 2023 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
Pure-quartic solitons (PQSs), as a specific type of solitons in the nonlinear system, are generated through a balance between the fourth-order dispersion and nonlinear effects [1–4]. PQSs in fiber lasers possess outstanding advantages in energy scaling compared to traditional quadratic solitons. Furthermore, they are relatively insensitive to noise and perturbations, which potentially makes them useful for long-distance communication applications. However, the lack of a suitable experimental platform, such as a piece of fiber with both quadratic and cubic dispersion close to zero, makes limited research in this field. Magnus et al. numerically presented an exact localized analytical solution to the pulse propagation in a piece of fiber that includes the negative second-order and fourth-order dispersion [5]. Michel et al. considered the effect of third-order dispersion on the optical soliton propagation in the presence of fourth-order dispersion [6]. Until 2020, the first PQS fiber laser with both net second-order and third-order dispersion being zero was experimentally demonstrated, using a spatial light modulator to adjust the net cavity dispersion. It implies that the pulse energy of the PQS scales with the third power of the inverse pulse duration [7]. Since then, a burgeoning corpus of research about PQSs in fiber lasers has been conducted such as the pulsating PQS soliton [8,9], the dissipative PQS [10], the creeping and erupting PQS soliton [11], the dissipative soliton resonance [12], etc.
Due to the material imperfection and fabrication imperfection, intrinsic birefringence always exists in optical fibers. The walk-off effect between the orthogonal polarization solitons typically occurs in non-polarization-maintaining fibers with weak birefringence, while the soliton trapping effect can shift their frequency in opposite directions and exchange their energy through four-wave mixing and cross-phase modulation (XPM), resulting in the formation of a vector soliton in this case [13–16]. Numerous types of vector solitons have been numerically predicted and experimentally demonstrated, such as the group-velocity-locked vector solitons [17,18], polarization-locked vector solitons [19,20], polarization rotation vector solitons [21,22], domain wall vector solitons [23,24]. In principle, a vector soliton is a self-sustaining entity with both polarization and spatial degrees of freedom, giving rise to many interesting phenomena that can be observed by adjusting the cavity-birefringence intensity [25–27]. As for the PQS, its long oscillating tail in the time domain can easily promote the XPM between the two orthogonal components in a weak-birefringence fiber, resulting in fascinating and unique nonlinear dynamics. Consequently, there is a great potential for conducting research on the generation and internal motion of the vector PQS in fiber lasers.
In this study, a numerical model of the vector PQS fiber laser is established based on the coupled Ginzburg-Landau equations. Starting from a pair of weak pulses with orthogonal polarization states, the stable vector PQS can be obtained with weak cavity-birefringence. The simulation results show that the gain competition caused by the XPM effect between the two orthogonal components plays a crucial role in the vector PQS buildup process, which is demonstrated by the obvious complementary energy oscillation processes. Enhancing the strength of the cavity-birefringence, both elliptical-polarization vector PQSs and near-linear-polarization PQSs can be generated. Furthermore, the transition from vector PQSs to scalar PQSs has also been observed due to the imbalance between walk-off and the trapping effect between the two orthogonal components. The results indicate that weakly birefringent fibers also support the generation and propagation of the vector PQSs. The gain competition between the two orthogonal components caused by the XPM effect remains the primary reason for the formation of vector solitons, and the strength of cavity-birefringence dominates the buildup dynamics of vector PQSs.
2. Numerical methods
Our simulation model, as shown in Fig. 1, is composed of six parts in a fiber laser, which includes a pump laser source, a wavelength division multiplexer (WDM), a 2-m Erbium-doped fiber (EDF), a saturable absorber (SA), a 2-m passive single-mode fiber (SMF), and an optical coupler (OC). The pump light is injected into the cavity through the WDM for signal amplification in the EDF. The SA is adopted to initialize the mode-locking operation. 10% of the intra-cavity power is emitted from the OC.
Fig. 1. Configuration of the vector PQS fiber laser. WDM, wavelength division multiplexer; EDF, Erbium-doped fiber; SA, saturable absorber; SMF, single-mode fiber; OC, optical coupler.
The laser pulse propagates in the cavity is governed by the coupled Ginzburg-Landau equations [28–30]:
Here Ax and Ay are the slowly varying amplitudes of two orthogonal polarized pulse envelopes along slow axis and fast axis of fibers. z and T are the propagation distance and the transformation time for the retarded frame. Our simulation is based on the fiber with linear birefringence and includes the second-order dispersion β2, the third-order dispersion β3, and the fourth-order dispersion β4. Δβ and Δβ1 are related to the linear birefringence which can be written as:
where the w0 is the central optical frequency, c is the light speed and Bm is the strength of modal birefringence defined as the difference between the modal refractive indices along the slow axis x and the fast axis y:The coupled Ginzburg-Landau equations are numerically solved using the fourth-order Runge–Kutta method in the interaction picture [31]. Apart from the birefringence related parameters, the two axes share the same values for other parameters, including the fiber loss α, the gain bandwidth Ωg, the nonlinear coefficient γ, and the saturable gain g of the fiber, which can be written as [32,33]:
Table 1. Parameters used in the simulations.
3. Results and discussions
3.1 Elliptical-polarization vector PQSs
Adjustment of the cavity-birefringence has a significant impact on the generation and evolution dynamics of vector solitons. To observe the buildup dynamics of the vector PQSs, different values of the modal birefringence Bm have been adopted. In this study, we first investigate the evolution dynamics of the vector soliton when Bm is 0.1 × 10−6, as shown in Fig. 2. The simulation starts with the initial weak pulse: A = sqrt(P0) × sech(t/tos), where P0 is 0.0001 W and tos is 2 ps.
Fig. 2. The pulse evolution dynamics in (a) the slow axis and (b) the fast axis. (c) The pulse shapes output from the 2000th round trip. The spectra evolution dynamics in (d) the slow axis and the fast axis. (f) The spectra output from the 2000th round trip. Bm is 0.1 × 10−6.
We observe that the pulse intensity output from the two orthogonal axes increases rapidly at the beginning. However, the two orthogonal pulses exhibit pulsating dynamics with energy complementarity in the time domain and gradually become stable after the 1500th round trip. The evolution of the spectra corresponding to Figs. 2(a) and 2(b) is shown in Figs. 2(d) and 2(e), which have a similar dynamic as observed in the time domain. Figures 2(c) and 2(f) depict the pulse and spectrum shapes at the 2000th round trip. It is evident that the laser in the two orthogonal directions has nearly the same shape and intensity in the time and frequency domains, indicating that a coherent coupling process of stable PQSs has been established [36].
In order to have a deep insight into the vector PQS generation, the energy and the phase evolution of the vector PQS versus the round trip number is plotted in Fig. 3. The energy of the two solitons increases quickly to 0.3 nJ around the 20th round trip. And then, the total energy of these two solitons nearly remains the same. However, the energy of the pulse in the fast axis starts to decrease, while that in the slow axis increases in a symmetrical trend. When the energy of the pulse in the slow axis increases to the first vertex around 0.58 nJ at the 186th round trip, the pulse in the fast axis also reaches a minimum value (0.04 nJ). Subsequently, the energy of the pulse in the slow axis features multiple minimum-maximum values and becomes stable after the 1500th round trip. Similarly, the pulse energy in the fast axis exhibits multiple maximum-minimum values and tends to be stable after the 1500th round trip. Notably, the energy dynamics of the two axes are complementary, as the total energy of the two solitons remains constant after the 20th round trip. The evolution process reveals that the buildup of the vector PQS is dominated by gain competition induced by the XPM effect, and that a balance within the pulse packet is formed during energy exchange [37]. After the 1500th round trip, the gain competition forces the two PQSs to have the same energy and constant phase difference. Furthermore, the two PQSs keep a synchronous evolution in the time domain because both positions are near 0 ps all the time. This is because the walk-off between the two pulses has been fully compensated by the soliton trapping effect. The evolution of phase difference between pulses in two axes closely corresponds to the changes in energy relationship. As we can see, the turning points of phase evolution correspond to the intersection points of pulse energy, which is because the nonlinear phase shift of the two orthogonal pulses is directly proportional to the pulse energy. Considering the final output intensity of the two axes and the fixed phase difference, we refer to the PQSs as “elliptical-polarization vector PQSs”.
Fig. 3. (a) The energy of PQSs output from the fast/slow axis and their total energy, and (b) the phase difference evolution dynamics with the increase of the round trips when Bm is 0.1 × 10−6.
Enhanced Bm of 1.0 × 10−6 leads to the generation of another type of “elliptical-polarization vector PQS”, as depicted in Fig. 4. It is obvious that the two PQSs experience the similar amplification and pulsating dynamics before the stable status is achieved. However, their spectrum shapes are different, indicating the incoherently coupled elliptical-polarization vector PQSs have been obtained [38]. The underlying physics is primarily associated with walk-off and soliton trapping between the two pulses. The cavity-birefringence directly results in the walk-off between the two PQSs in the time domain owing to the different group velocity delay. At the beginning of the evolution (before the 23rd round trip), the weak soliton trapping effect resulting from the weak energy of the two pulses cannot bind them together. Consequently, the two pulses evolve with the mismatched time delay and show opposite directions in the time domain, i.e., the pulses in the slow and fast axes evolve in separate negative and positive time directions, as shown in Fig. 5. However, the increasing pulse energy leads to the accumulation of nonlinear interactions, and gradually compensating for the walk-off with the enhanced soliton trapping effect. Hence, the two pulses have nearly identical evolution traces in the time domain with respect to the round trip number. Around the 180th round trip, the energy of the pulse in the slow axis exceeds that in the fast axis, which results in the pulse in the fast axis gradually following the evolving direction with the pulse in the slow axis and moving towards the negative time direction. With the round trip number increased to ∼500, the pulse in the fast axis regains control of the evolving direction of the pulse packet as its energy is larger than that in the slow axis. Due to the combined effects of walk-off and soliton trapping, the temporal positions and the spectrum shapes of the two orthogonal PQSs exhibit significant differences.
Fig. 4. (a) The pulse and (b) spectrum evolution dynamics in the slow axis. (c) The pulse shapes output from the 7000th round trip. (d) The pulse and (e) spectrum evolution dynamics in the fast axis. (f) The spectra output from the 7000th round trip. Bm is 1 × 10−6.
Fig. 5. The detailed evolution dynamics of the two pulses in (a) the slow axis and (b) the fast axis in the time domain within 1200 round trips.
3.2 Near-linear-polarization vector PQSs
The situation changes significantly when the cavity-birefringence continues to increase. For example, Fig. 6 illustrates another type of evolution dynamics for vector PQSs in the time domain when the Bm is 1.5 × 10−6. Both PQSs show a creeping dynamic toward the negative time direction in Figs. 6(a) and 6(b). However, it is interesting that the two PQSs have the completely different evolution dynamics. The PQS in the slow axis initially evolves towards the left with a slow speed, but then turns quicker after the 200th round trip. On the other hand, the PQS in the fast axis undergoes energy amplification before the 200th round trip and then gradually decreases in energy, but the time direction evolves first towards the right and then towards the left. As mentioned in the previous section, the evolution direction in the time domain is determined by the pulse with the stronger energy. Initially, the two pulses evolve in opposite directions in the time domain because of birefringence. With the increase of the round trips, the increased energy of the pulse in the slow axis and the energy exchange between the two pulses can enhance the soliton trapping effect between them, which can gradually compensate for the walk-off. However, the energy exchange between the two axes seems irreversible in this situation. Although the creeping speed of the pulses on the two axes is mutually restrained before the 300th round trip, afterwards, the pulse in the fast axis is completely captured by the pulse in the slow axis, which has taken away almost all the energy. Considering the fact that the energy of the pulse in the slow axis is far greater than that of the fast axis, the pulse in the fast axis will follow the evolution direction with the pulse in the slow axis, i.e., both pulses evolve together towards the negative time direction.
Fig. 6. The pulse evolution versus the round trip number in (a) the slow axis and (b) the fast axis. (c) The pulse shapes and (d) spectra output from the 1000th round trip. Bm is 1.5 × 10−6.
Figure 7 shows the pulse energy evolution in the two axes. As mentioned above, the pulse energy in the two axes increases quickly in the first 20 round trips. Then the pulse energy evolves into a complementary pulsating dynamic and tends to stabilize. Thus, in this case, the gain competition also helps to redistribute the energy in the vector PQSs packet. However, the increased cavity-birefringence causes a different gain-competition process with the case of Bm =0.1 × 10−6. The PQS in the slow axis has an advantage in obtaining energy over the PQS in the fast axis, resulting in higher pulse energy. Since the pulse energy in the slow axis is much higher than that in another axis at the end, we can name it as the “near-linear polarization vector PQS”.
Fig. 7. The energy of the PQSs output from the fast/slow axis and their total energy evolution dynamics with the increase of the round trips when the Bm is 1.5 × 10−6.
3.3 Transition from vector PQSs into scalar PQSs
The cavity-birefringence can bring a group-velocity difference between the two solitons, while for the vector soliton, the soliton trapping effect can bind them together. The above results show that the walk-off between the two PQSs induced by the group-velocity difference has been compensated by the trapping effect between them. However, the higher the cavity-birefringence, the larger the walk-off. Therefore, it is interesting to observe what will happen with larger cavity-birefringence.
Figure 8 presents the evolution of the two PQSs in the slow/fast axis when Bm =3 × 10−6. It can be seen that the intensity of the PQSs increases quickly within 22 round trips. And then, the energy difference increases gradually until the pulse intensity in the slow axis is nearly 1.6 × 109 times larger than that in the fast axis, because the energy in the slow axis continues to increase, while the pulse energy in the fast axis starts to decrease. However, it should be noted that the walk-off between the two solitons starts quickly since the 22nd round trip, indicating that energy exchange decreases with the increase of round trips despite the existence of oscillating tails in the time domain. The energy evolution is similar to that in the case of Bm =1.5 × 10−6, but the two PQSs consistently maintained their initial creeping directions in the time domain. The more round trips, the less the overlap between the two PQSs. As the two pulses completely move away from each other, the XPM effect will disappear and the pulse in one axis will fail to capture the pulse in the other axis. Eventually, in this case, the uncompensated walk-off and the huge energy difference between the two solitons demonstrate that the vector PQS have evolved into the scalar PQS in the end.
To gain deeper insights into PQS evolution in birefringent fibers and investigate what the difference is between the conventional nonlinear Schrodinger equation soliton (NLSES) and the PQS, we have made another simulation of the NLSES for comparison, as shown in Fig. 9. In the NLSES simulation model, β2-EDF = β2-SMF = -10 ps2/km and β4-EDF = β4-SMF = 0 ps4/km. The other parameters are the same as those in the case for elliptical-polarization vector PQSs with Bm of 0.1 × 10−6. It can be observed that even after undergoing the same roundtrip (i.e., 2000 roundtrips), a stable vector soliton state cannot be achieved. The temporal positions along the two axes nearly remain the same, which indicates the two pulses have trapped together. From the perspective of pulse energy on both axes, the energy exchange between the two axes has been successfully established. The evolution dynamic is similar to the PQS evolution with Bm of 1 × 10−6 described in Fig. 5 while with a slower speed.
To construct the vector soliton, it is crucial to not only concentrate on the formation of the soliton itself, but also consider the interaction between two orthogonal polarization solitons and their ability to develop into a stable bound state. So, vector solitons primarily rely on the trapping effect via XPM between the orthogonal polarization solitons, while the presence of birefringence can induce the walk-off effect and weaken the trapping. The key point here is the degree of the pulse interaction induced by XPM during the soliton formation. Through our simulations, we have discovered that a significant characteristic distinguishing the PQS from the conventional NLSES is their intensity distribution and the long oscillating tail feature. This distinction becomes more evident when observed on a log scale, as shown in Figs. 10(a) and 10(b). It is important to note that when two orthogonal PQSs interact through XPM, the energy distribution of these oscillations in the tail can mutually influence each other. In other words, through XPM, the two orthogonal PQSs actively participate in each other's formation process. So, we observe similar characteristics in terms of oscillation period and amplitude (inserted picture in Fig. 10(b)). With the enhanced birefringence for the near-linear polarization PQS, the degree of their interaction still ensures the trapping of the pulses on the two axes with an increased energy difference. The major energy is concentrated in the slow axis, and the pulse in the slow axis exhibits characteristics of a PQS with the oscillating tail. The pulse in the fast axis has lost its symmetry due to XPM, but a portion of the tail can still be observed on its left side in Fig. 10(c). When the soliton trapping fails with higher birefringence, the transition from vector PQSs into scalar PQSs will occur, as shown in Fig. 10(d).
Fig. 10. The pulse shapes of the last round trip in log scale for (a) NLSES, (b) the elliptical-polarization vector PQS, (c) the near-linear-polarization vector PQS, and (d) the scalar PQS.
Additionally, according to the formation mechanism [3], the maximum energy of the PQS is usually higher than that of the conventional NLSES under similar cavity structure, which will lead to stronger XPM. And the oscillating tail of PQSs can indeed facilitate XPM. Thus, we can conclude that it is much easier to obtain the stable vector PQS state than the NLSES case.
4. Conclusions
In summary, our numerical simulations reveal the formation and evolution dynamics of stable vector PQSs in the mode-locked fiber laser with weak cavity-birefringence. With appropriate birefringence strengths, both elliptical-polarization vector PQSs and near-linear-polarization vector PQSs can be obtained, which is determined by the results of the gain competition between the two orthogonal components. Moreover, the further enhanced cavity-birefringence can result in the vector PQSs evolving into scalar PQSs because the walk-off between the two solitons cannot be fully compensated by the soliton trapping effect. Our findings open new perspectives for vector soliton pulse evolution dynamics and will bring new insights into the design and application of vector PQS fiber lasers.
Funding
National Natural Science Foundation of China (62175230, 62225507, U2033211); the CAS Project for Young Scientists in Basic Research (YSBR-065); Innovation Fund Denmark (2105-00039B); Villum Fonden (00037822).
Disclosures
The authors declare no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
References
1. A. B. Redondo, C. M. Sterke, J. E. Sipe, T. F. Krauss, B. J. Eggleton, and C. Husko, “Pure-quartic solitons,” Nat. Commun. 7(1), 10427 (2016). [CrossRef]
2. D. Grassani, E. Tagkoudi, H. Guo, C. Herkommer, F. Yang, T. J. Kippenberg, and C. Brès, “Mid infrared gas spectroscopy using efficient fiber laser driven photonic chip-based supercontinuum,” Nat. Commun. 10(1), 1553 (2019). [CrossRef]
3. C. M. Sterke, A. F. J. Runge, D. D. Hudson, and A. B. Redondo, “Pure-quartic solitons and their generalizations-Theory and experiments,” APL Photonics 6(9), 091101 (2021). [CrossRef]
4. K. K. K. Tam, T. J. Alexander, A. B. Redondo, and C. M. Sterke, “Stationary and dynamical properties of pure-quartic solitons,” Opt. Lett. 44(13), 3306–3309 (2019). [CrossRef]
5. M. Karlsson and A. Höök, “Soliton-like pulses governed by fourth order dispersion in optical fibers,” Opt. Commun. 104(4-6), 303–307 (1994). [CrossRef]
6. M. Piché, J. F. Cormier, and X. Zhu, “Bright optical soliton in the presence of fourth-order dispersion,” Opt. Lett. 21(12), 845–847 (1996). [CrossRef]
7. A. F. J. Runge, D. D. Hudson, K. K. K. Tam, C. M. Sterke, and A. B. Redondo, “The pure-quartic soliton laser,” Nat. Photonics 14(8), 492–497 (2020). [CrossRef]
8. Z. X. Zhang, M. Luo, J. X. Chen, L. H. Chen, M. Liu, A. P. Luo, W. C. Xu, and Z. C. Luo, “Pulsating dynamics in a pure-quartic soliton fiber laser,” Opt. Lett. 47(7), 1750–1753 (2022). [CrossRef]
9. S. Yang, Z. Zhu, Y. Qi, L. Jin, L. Li, and X. Lin, “Internal motion within pulsating pure-quartic soliton molecules in a fiber laser,” Chaos, Solitons Fractals 172, 113544 (2023). [CrossRef]
10. Z. C. Qian, M. Liu, A. P. Luo, Z. C. Luo, and W. C. Xu, “Dissipative pure-quartic soliton fiber laser,” Opt. Express 30(12), 22066–22073 (2022). [CrossRef]
11. Y. Han, B. Gao, G. Wu, J. Huo, H. Wen, Y. Li, L. Liu, and C. Ma, “Creeping and erupting dynamics in a pure-quartic soliton fiber laser,” Opt. Express 31(2), 1787–1798 (2023). [CrossRef]
12. Y. Zhang, C. Jin, C. Tao, S. Luo, Q. Ling, Z. Guan, D. Chen, and Y. Cui, “Dissipative pure-quartic soliton resonance in an Er-doped fiber laser,” Opt. Commun. 538, 129479 (2023). [CrossRef]
13. Y. Qi, S. Yang, J. Wang, L. Li, Z. Bai, Y. Wang, and Z. Lv, “Recent advance of emerging low-dimensional materials for vector soliton generation in fiber lasers,” Mater. Today Phys. 23, 100622 (2022). [CrossRef]
14. M. N. Islam, C. D. Poole, and J. P. Gordon, “Soliton trapping in birefringent optical fibers,” Opt. Lett. 14(18), 1011–1013 (1989). [CrossRef]
15. Y. Q. Du, X. W. Shu, and P. Y. Cheng, “Numerical simulations of fast-axis instability of vector solitons in mode-locked fiber lasers,” Opt. Express 25(2), 1131–1141 (2017). [CrossRef]
16. X. Wang, Z. Wang, Y. Song, D. Fan, and J. Liu, “Real-time dynamics of noise-like vector pulses in a figure-eight fiber laser,” Opt. Express 30(6), 9137–9146 (2022). [CrossRef]
17. J. H. Lin, Z. Y. Zhang, Z. Y. Li, P. C. Peng, Y. F. Song, and H. Zhang, “Abundant dynamics of group velocity locked vector solitons from Er-doped fiber laser based on GO/PVA film,” Nanophotonics 11(21), 4831–4841 (2022). [CrossRef]
18. Y. Luo, L. Li, D. Liu, Q. Sun, Z. Wu, Z. Xu, D. Tang, S. Fu, and L. Zhao, “Group velocity locked vector dissipative solitons in a high repetition rate fiber laser,” Opt. Express 24(16), 18718–18726 (2016). [CrossRef]
19. S. T. Cundiff, B. C. Collings, N. N. Akhmediev, J. M. S. Crespo, K. Bergman, and W. H. Knox, “Observation of Polarization-Locked Vector Solitons in an Optical Fiber,” Phys. Rev. Lett. 82(20), 3988–3991 (1999). [CrossRef]
20. L. M. Zhao, D. Y. Tang, X. Wu, H. Zhang, and H. Y. Tam, “Coexistence of polarization-locked and polarization-rotating vector solitons in a fiber laser with SESAM,” Opt. Lett. 34(20), 3059–3061 (2009). [CrossRef]
21. Y. F. Song, H. Zhang, D. Y. Tang, and D. Y. Shen, “Polarization rotation vector solitons in a graphene mode-locked fiber laser,” Opt. Express 20(24), 27283–27289 (2012). [CrossRef]
22. S. V. Sergeyev, M. Eliwa, and H. Kbashi, “Polarization attractors driven by vector soliton rain,” Opt. Express 30(20), 35663–35670 (2022). [CrossRef]
23. S. Yang, Q. Zhang, Z. Zhu, P. Yin, Y. Ge, L. Li, L. Jin, L. Zhang, and H. Zhang, “Recent advances and challenges on dark solitons in fiber lasers,” Opt. Laser Technol. 152, 108116 (2022). [CrossRef]
24. M. Liu, Y. Qi, S. Yang, Z. Bai, B. Yan, J. Ding, Y. Wang, and Z. Lu, “Switchable L-band dual-wavelength dark–bright pulse pair generation from an Er-doped mode-locked fiber laser with SMF–GIMF–SMF as the saturable absorber,” Appl. Phys. B 128(10), 190 (2022). [CrossRef]
25. K. Zhao, C. Gao, X. Xiao, and C. Yang, “Vector quartic solitons in birefringent fibers,” Opt. Lett. 46(4), 761–764 (2021). [CrossRef]
26. C. B. Tabi, H. Tagwo, C. G. L. Tiofack, and T. C. Kofané, “Pure quartic modulational instability in weakly nonlocal birefringent fibers,” Opt. Lett. 47(21), 5557–5560 (2022). [CrossRef]
27. C. M. Sterke and A. B. Redondo, “Even-order dispersion solitons: A pedagogical note,” Opt. Commun. 541, 129560 (2023). [CrossRef]
28. G. P. Agrawal, Nonlinear Fiber Optics, 5th ed. (Academic, 2013).
29. X. Hu, J. Guo, J. Wang, J. Ma, L. Zhao, S. Yoo, and D. Tang, “Novel optical soliton molecules formed in a fiber laser with near-zero net cavity dispersion,” Light: Sci. Appl. 12(1), 38 (2023). [CrossRef]
30. L. Jin, Q. Zhang, B. Zhang, Z. Gao, S. Yang, and L. Li, “Numerical analysis of hybrid mode-locking stability in a Ho-doped fiber laser,” Opt. Express 32(2), 1141–1153 (2023). [CrossRef]
31. J. Hult, “A fourth-order Runge-Kutta in the interaction picture method for simulating supercontinuum generation in optical fibers,” J. Lightwave Technol. 25(12), 3770–3775 (2007). [CrossRef]
32. J. Zheng, S. Yang, Z. Zhu, K. Lau, and L. Li, “72-fs Er-doped Mamyshev oscillator,” J. Lightwave Technol. 40(7), 2123–2127 (2022). [CrossRef]
33. Y. Wang, L. Su, S. Wang, L. Hua, L. Li, D. Shen, D. Tang, A. Komarov, M. Klimczak, S. Fu, M. Tang, X. Tang, and L. Zhao, “Breach and recurrence of dissipative soliton resonance during period-doubling evolution in a fiber laser,” Phys. Rev. A 102(1), 013501 (2020). [CrossRef]
34. S. Yang, Q. Zhang, Z. Zhu, Y. Qi, L. Li, and X. Lin, “Self-parametric amplification of a bound-state dissipative soliton assisted by spectral filtering effect in the negative dispersion regime,” Chaos, Solitons Fractals 164, 112733 (2022). [CrossRef]
35. S. Yang, F. Li, M. Gong, L. Zhang, Z. Zhu, H. Shen, and S. Chen, “Generation of Q-switched and mode-locked pulses based on PbS/CdS saturable absorbers in an Er-doped fiber laser,” J. Mater. Chem. C 10(15), 5956–5961 (2022). [CrossRef]
36. H. Zhang, D. Y. Tang, L. M. Zhao, and N. Xiang, “Coherent energy exchange between components of a vector soliton in fiber lasers,” Opt. Express 16(17), 12618–12623 (2008). [CrossRef]
37. D. Y. Tang, L. M. Zhao, B. Zhao, and A. Q. Liu, “Mechanism of multisoliton formation and soliton energy quantization in passively mode-locked fiber lasers,” Phys. Rev. A 72(4), 043816 (2005). [CrossRef]
38. H. Zhang, D. Tang, L. Zhao, Q. Bao, and K. P. Loh, “Vector dissipative solitons in graphene mode locked fiber lasers,” Opt. Commun. 283(17), 3334–3338 (2010). [CrossRef]
