Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Spatiotemporal dissipative solitons and vortices in a multi-transverse-mode fiber laser

Open Access Open Access

Abstract

We introduce a model for spatiotemporal modelocking in multimode fiber lasers, which is based on the (3+1)-dimensional cubic-quintic complex Ginzburg-Landau equation (cGLE) with conservative and dissipative nonlinearities and a 2-dimensional transverse trapping potential. Systematic numerical analysis reveals a variety of stable nonlinear modes, including stable fundamental solitons and breathers, as well as solitary vortices with winding number n = 1, while vortices with n = 2 are unstable, splitting into persistently rotating bound states of two unitary vortices. A characteristic feature of the system is bistability between the fundamental and vortex spatiotemporal solitons.

© 2019 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction and the model

New physical principles for generating and controlling ultrashort pulses are under constant scrutiny of the scientific community. A recent prominent topic of the research in this area is complex multi-mode spatiotemporal dynamics in optical fibers used for both generation and delivery of short pulses, see, e.g., [14] for review. As intensity of light is increased through tighter focusing in both time and space, it starts interacting with matter nonlinearly, which is crucially important for the generation of ultrashort pulses and their spectral and temporal shaping. While the structure of the short pulses is intrinsically multimode (in a cavity) and multi-wavelength in the propagation direction (in free space), their transverse multimode structure has been traditionally considered as a thing to avoid and a special care is commonly taken to ensure that a pulse is carried by a single transverse mode. Nevertheless, transversely multimode ultrashort pulses have recently become a topic of intense research. Their spatial complexity is a feature to be controlled and utilized [2], fiber-based devices being the most prominent setting in this context. In particular, multi-transverse-mode passive fibers have been used for boosting transmission capacity of communication lines by means of spatial-division multiplexing [5]. Strongly nonlinear propagation of pulses in graded-index (GRIN) fibers with a transversely parabolic refractive index has been investigated in the context of the beam-profile self-cleaning, as well as the generation of continuum radiation and solitons [614]. Spatiotemporal mode-locking (STML) of longitudinal and transverse modes in a fiber-laser system was recently reported in [41], with follow-up results contributing to this area. Those include observations of solitons and soliton bound states [15,16], a proposal for the multi-transverse-mode gain-saturation mechanism [17], and direct generation of vortex beams [18]. The latter result is directly relevant and motivational for the present work, though the above-mentioned recent STML results predominantly dealt with the normal-GVD (group-velocity-dispersion) case, while we are address the case of anomalous GVD in this work. Previously, this case has been chiefly considered in studies with suppressed transverse effects [1]. Multimode fibers are closely related to multi-core ones, since coupled cores introduce multiple transverse modes, whose number is equal the numbers of cores. Pulse compression, light bullets and spatiotemporal discrete vortices all have been studied in the multi-core fibers [1921]. STML in various multi-core waveguide systems has also been addressed by several groups [2224]. It is relevant to mention that studies of STML and soliton effects in (non-fiber) laser cavities had started decades ago—see, e.g., [2527].

Mode-locking effects and formation of dissipative solitons have been traditionally (and successfully) modeled using various realizations of the complex Ginzburg-Landau equation (cGLE), its ubiquitous version being one with the cubic-quintic nonlinearity [1,28]. The quintic term supports the subcritical transition to the lasing state, thereby providing conditions for mode-locking and the creation of stable solitons at subcritical pumping levels [28]. Locking of the longitudinal modes in single-transverse-mode fiber lasers is described by (1+1)-dimensional cGLE, while the consideration of multi-mode fibers opens the way to the realm of the (3+1)D cGLE. The lateral confinement of the fiber modes immediately requires inclusion of the trapping potential acting in two transverse dimensions, which is a basic parabolic potential for GRIN fibers. This potential is a feature creating a new playground relative to the quite extensive studies of "light bullets" [2932] and vortex solitons [3337] in (3+1)D cGLE. Below we use the (3+1)D cGLE including conservative and dissipative cubic-quintic terms and a 2D parabolic potential.

Measuring the propagation distance along the fiber, beam’s width, and time in units of suitably defined scales $L$, $w$ and $\tau$, the cGLE is written as

$$\begin{array}{r}i\frac{\partial \psi }{\partial z}+\frac{L}{2\tau ^{2}}k^{\prime \prime }\left( 1-i\beta \right) \frac{\partial ^{2}\psi }{\partial {t}^{2}}+\frac{L }{2k_{0}w^{2}}\left( \frac{\partial ^{2}}{\partial {x}^{2}}+\frac{\partial ^{2}}{\partial {y}^{2}}\right) \psi -\frac{k_{0}gw^{2}L}{2}\left( {x}^{2}+{y} ^{2}\right) \psi +iL\varepsilon ^{\prime }\psi\\ +L\left( \eta ^{\prime }-i\alpha ^{\prime }\right) |\psi |^{2}\psi +L\left( \nu ^{\prime }+i\mu ^{\prime }\right) |\psi |^{4}\psi =0.\end{array}$$
Here $z$, $t$ and $x,y$ are the dimensionless distance, the temporal coordinate, and transverse coordinate, and time is defined in the reference frame moving with the pulse’s group velocity. Further, $k^{\prime \prime }$ is the GVD coefficient, $\beta$ the spectral-filtering coefficient, $k_{0}$ the propagation constant of the fundamental mode, $g$ the GRIN parameter, and $\varepsilon \equiv L\varepsilon ^{\prime }>0$ the dimensionless linear loss over length $L$. If $\psi$ is the dimensionless field envelope, then the scaled (dimensionless) parameters $\eta =L\eta ^{\prime }$, $\alpha =L\alpha ^{\prime }>0$, $\nu =L\nu ^{\prime }$ and $\mu =L\mu ^{\prime }>0$ determine ratios between $L$ and the nonlinear lengths associated with the Kerr effect, cubic gain, quintic nonlinearity, and nonlinear saturation, respectively. We now fix $L=1/(k_{0}gw^{2})$, $w^{2}=L/k_{0}$ and $\tau ^{2}=Lk^{\prime \prime }$ to write the spatiotemporal cGLE (2) with the cubic-quintic terms in the scaled form:
$$i\frac{\partial \psi }{\partial z}+\frac{1}{2}\left( 1-i\beta \right) \frac{ \partial ^{2}\psi }{\partial t^{2}}+\frac{1}{2}\left( \frac{\partial ^{2}}{ \partial x^{2}}+\frac{\partial ^{2}}{\partial y^{2}}\right) \psi -\left( x^{2}+y^{2}-i\varepsilon \right) \psi +\left( \eta -i\alpha \right) |\psi |^{2}\psi +\left( \nu +i\mu \right) |\psi |^{4}\psi =0.$$
Nontrivial solutions can be supported by the cubic gain in the competition with the linear and quintic losses, if the gain coefficient exceeds a threshold value [38,39]
$$\alpha =A\alpha _{\mathrm{thr}},~\alpha _{\mathrm{thr}}\equiv \left( 2\sqrt{ \varepsilon \mu }\right)$$
with $A>1$. We here fix values $A=1.5$ and $\beta =\nu =0$, which adequately represent the generic case, while the linear loss and nonlinear gain, $\varepsilon$ and $\mu$, are varied.

Thus, in this work we focus on predicting STML scenarios in the setup based on the interplay of the cubic-quintic dissipative nonlinearity with the anomalous GVD, transverse diffraction, parabolic GRIN potential, and Kerr self-focusing, in the absence of temporal and spatial filtering. In fact, the prediction of patterns which are stable without the help of the filtering is an essential novelty of the analysis.

2. Numerical results

We look for stationary localized solutions, in the form of $\psi (x,y,t,z)=\phi (x,y,t)\exp \left ( iqz\right )$, where $q$ is a real propagation constant, and function $\phi (x,y,t)$ satisfies the equation following from the substitution of this ansatz in Eq. (2) with $\beta =\nu =0$:

$$-q\phi +\frac{1}{2}\left( \frac{\partial ^{2}}{\partial t^{2}}+\frac{ \partial ^{2}}{\partial x^{2}}+\frac{\partial ^{2}}{\partial y^{2}}\right) \phi -\left( x^{2}+y^{2}-i\varepsilon \right) \phi +\left( \eta -i\alpha \right) |\phi |^{2}\phi +i\mu |\phi |^{4}\phi =0.$$
Splitting complex Eq. (4) into real and imaginary parts, solutions for localized states have been produced by means of the modified squared-operator method [40]. The solutions are characterized by the value of their total energy, i.e., integral norm,
$$P=\int \int \int \left\vert \psi (x,y,t)\right\vert ^{2}dxdydt.$$

2.1 Fundamental dissipative solitons

To determine numerical solutions and the propagation constant simultaneously, the Gaussian input was used, $\phi _{0}=1.5~\mathrm {sech} \left ( x^{2}+y^{2}\right ) \mathrm {\exp }(-1/2t^{2})$. The existence and stability of numerically found fundamental (zero-vorticity) dissipative solitons, produced by Eq. (4), is summarized in Fig. 1(a), which displays the numerically found $q$ and $P$ (integral power) vs. the linear-loss parameter, $\varepsilon$, for a fixed value of the quintic-loss coefficient, $\mu =1$, while the cubic gain is fixed as per Eq. (3). The stability was identified by means of systematic simulations of perturbed evolution of the solitons in the framework of Eq. (2). Figure 1(b) summarizes the results by plotting stability boundaries for the dissipative solitons in the plane of $\left ( \mu ,\varepsilon \right )$. The shape and evolution of the stable solitons are illustrated by Fig. 2.

 figure: Fig. 1.

Fig. 1. (a) Numerically computed propagation constant $q$ and integral power $P$ [see Eq. (5)] of fundamental three-dimensional solitons versus $ \varepsilon$, for fixed $ \mu =1$ in Eq. (3). Red, blue and black segments correspond to stable fundamental solitons, moving breathers, and eventually destructed modes, respectively. (b) The stability chart for the stationary fundamental solitons and modes into which unstable solitons are spontaneously transformed. Indicated inside the stability area of solitons is the bistability region, in which solitary vortices with winding numbers are stable too, cf. Figure 5(b).

Download Full Size | PDF

 figure: Fig. 2.

Fig. 2. (a) The evolution of a stable fundamental soliton, for $ \mu =1, \varepsilon =0.5$, with $q=-0.8635,P=6.9898$, is displayed by means of isosurfaces of $\left \vert \psi \left ( x,y,z,t\right ) \right \vert$. (b) Left and right panels display spatial and temporal cross sections of the soliton’s shape, $| \psi (x,0,z,0)|$ and $| \psi (0,0,z,t)|$, respectively.

Download Full Size | PDF

Weakly unstable fundamental solitons can be easily found too. They spontaneously transform into quiescent or moving breathers, see examples in Fig. 3 and 4, respectively. The breathers are robust dynamical modes, which keep integrity in the course of indefinitely long subsequent evolution. Their stability regions in the parameter plane of $\left ( \mu ,\varepsilon \right )$ are also drawn in Fig. 1(b). Note, in particular, that the region of destruction of all modes (i.e., onset of spatial “ turbulence” as a result of the simulated evolution) occupies, quite naturally, an area in which the loss and gain coefficients [recall that the latter one is taken as per Eq. (3)] are too large to maintain the stability of regular modes, such as solitons and breathers. The same peculiarity is observed in stability charts displayed below for vortex states in Figs. 5(b) and 10(b).

 figure: Fig. 3.

Fig. 3. (a) The evolution of a robust quiescent breather, into which an unstable fundamental soliton transforms for $ \mu =0.3, \varepsilon =2.8$, which corresponding to $q=-2.0446,P=5.3081$, shown by means of the isosurface of $\left \vert \psi \left ( x,y,z,t\right ) \right \vert$. (b) Left and right panels show the evolution in the spatial and temporal of cross sections, $| \psi (x,0,z,0)|$ and $| \psi (0,0,z,t)|$, respectively.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. A robust moving breather for $ \mu =1.0, \varepsilon =1.8$, into which the respective unstable soliton, with $q=-1.7927,P=4.5920$, spontaneously transforms. (a) The evolution in the spatial and temporal cross sections, $| \psi (x,0,z,0)|$ (left) and $| \psi (0,0,z,t)|$ (right). (b) The spatiotemporal dynamics with temporal motion.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. (a) The numerically found values of $q$ (top panel) and $P$ (bottom panel) for stationary vortices with winding number $n=1$ versus $ \varepsilon$ for fixed $ \mu =1.0$. Red, blue, magenta, yellow and black segments correspond to stable vortices, splitting into fundamental solitons, chaotically oscillating modes, moving vortex breathers, and destruction of unstable vortices, respectively. (b) The stability chart for vortices with $n=1$ in the plane of $\left ( \mu , \varepsilon \right )$. The central shaded area in (b) is populated by chaotically oscillating localized modes, see an example in Fig. 8. Note that the stability area of vortices is completely covered by the region in which the fundamental solitons are stable [cf. Fig. 1(b)], i.e., it is a bistability area.

Download Full Size | PDF

2.2 Vortices with winding number $n=1$

Localized vortex states with winding number (topological charge) $n=1$ were constructed as solutions of Eq. (4), starting with input, written in polar coordinates $\left ( r,\theta \right )$ (in the plane of $\left ( x,y\right )$) as $\phi _{0}=2.5r~\mathrm {sech}\left ( r^{2}\right ) ~ \mathrm {sech}\left ( t\right ) \exp \left ( i\theta \right )$. As a result, two coexisting vortex families, one unstable and one stable, have been found. Figure 5(a) displays the respective numerically found values of $q$ and integral power $P$ versus $\varepsilon$. In particular, the stable (red) branch was found in the interval of $0<\varepsilon <0.4$ for $\mu =1.0$. Unstable vortex states split into fundamental solitons (in the blue segment), or transform into chaotically oscillating modes (magenta) or moving breathers, which keep the initial vorticity (yellow), or suffer destruction (black). The results are summarized in Fig. 5(b) in terms of the $\left ( \mu ,\varepsilon \right )$ plane, where the shaded area pertains to chaotically oscillating modes. Typical examples of the stable vortices, splitting, chaotically oscillating modes, and moving vortex breathers are shown in Fig. 6, 7, 8 and 9 , respectively.

 figure: Fig. 6.

Fig. 6. (a) The isosurface evolution of a stable vortex found at $ \mu =1.0, \varepsilon =0.2$, which corresponding to $q=-2.3130,P=15.9142$. (b) Its amplitude and phase patterns (c) Left and right panels display the evolution of the spatial and temporal cross sections $| \psi (x,0,z,0)|$ and $| \psi (0,0,z,t)|$, respectively.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. The spontaneous splitting of an unstable vortex with $n=1$ into two fundamental solitons, at $ \mu =1.0, \varepsilon =0.6$. (a) The evolution of $| \psi |$ is shown in the $t$- and $x$-cross sections (top and bottom panels, respectively). (b) The temporal shape of the moving fundamental solitons, produced by the splitting, at $z=400$ (the blue line), fitted to the standard soliton’s shape, $\mathrm {sech}(t)$.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. (a) The evolution of $| \psi \left ( 0,0,z,t\right ) |$ and $| \psi \left ( x,0,z,0\right ) |$ in the $t$- and $x$- cross-sections, illustrating the spontaneous transformation of an unstable stationary vortex with $n=1$ into a chaotically oscillating mode at $ \mu =1.0, \varepsilon =0.9$. (b) The respective amplitude and phase patterns at $z=1000$ (eventual values of the amplitude are small, as a result of the evolution).

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. A moving breather, into which an unstable stationary vortex, with $q=-2.9276,P=9.9546$, is spontaneously transformed (keeping its vorticity) at $ \mu =1.0, \varepsilon =1.2$. The evolution of cross sections $| \psi \left ( x,0,z,0\right ) |$ and $| \psi \left ( 0,0,z,t\right ) |$ is displayed in the left and right panels, respectively.

Download Full Size | PDF

An obviously important conclusion, demonstrated by the comparison of Figs. 5(b) and 1(b), is the bistability between the fundamental solitons and solitary vortices, as the solitons’ stability area completely covers that of the vortices. In fact, this feature may be interpreted as tri-stability, as stable vortices may have winding numbers $+1$ and $-1$. This fact may find applications to optical switching in the multimode fibers, and is of considerable interest in terms of general properties of vortex solitons in dissipative media. The switching may be controlled by adjusting the power: it is seen in the bottom panels of Figs. 1(a) and 5(a) that the power of a stable vortex is $\simeq 3$ times larger than the power of a coexisting stable fundamental soliton.

2.3 Vortices with $n=2$

Double vortices, with winding number $n=2$, were constructed as solutions of Eq. (5), starting with input $\phi _{0}=2.5r~\mathrm {sech} \left ( r^{2}\right ) ~\mathrm {sech}\left ( t\right ) \exp \left ( 2i\theta \right )$. Systematic simulations have demonstrated that no double vortices are stable. They spontaneously develop into various modes, which are robust in different areas in the plane of $\left ( \mu ,\varepsilon \right )$, as shown in Fig. 10. These are rotating bound states of two unitary vortices, into which the double vortex splits, which is shown in Fig. 11; “ wobbling breathers” , which lose the original vorticity and perform spontaneous shuttle motion in the $x$ direction, coupled to intrinsic vibrations, see Fig. 12; and moving breathers (not shown here in detail), which lose the initial vorticity. In addition, there is a parameter region in which the double vortex splits into $2$, $3$, or $4$ breathers (each without intrinsic vorticity), see details in Fig. 13.

 figure: Fig. 10.

Fig. 10. The stability chart in the plane of $\left ( \mu , \varepsilon \right )$ for vortices with $n=2$. The shaded area corresponds to the “ wobbling breathers” .

Download Full Size | PDF

 figure: Fig. 11.

Fig. 11. Spontaneous splitting of an unstable vortex with $n=2$ into a stably rotating bound state of two vortices with $n=1$ at $ \mu =0.4$, $ \varepsilon =0.1$ (a) The evolution of the field in the temporal and spatial cross sections, $| \psi (0,0,z,t)|$ and $| \psi (x,0,z,0)|$. (b) Amplitude and phase patterns of $ \psi (x,y)$ in the plane of $t=0$ and (bottom) the temporal profile at the central point, $x=y=0$, plotted at $z=500$ (c) The top, left and right panels zoom in on the evolution of the absolute value of the field, and spatial and temporal cross sections, respectively, in the interval of $300<z<320$. The established evolution (rotation) period is $T_{z}=2.2$.

Download Full Size | PDF

 figure: Fig. 12.

Fig. 12. The “ wobbling breather” at $ \mu =0.4$, $ \varepsilon =0.6$ (a) The evolution of the field in the temporal and spatial cross sections, $| \psi (0,0,z,t)|$ and $| \psi (x,0,z,0)|$. (b) The top, left and right panels zoom in on the evolution of the absolute value of the field, and spatial and temporal cross sections in the interval of $z=200-220$, respectively. The established evolution period is $T_{z}=4.5$.

Download Full Size | PDF

 figure: Fig. 13.

Fig. 13. (a) Splitting of an unstable breather with $n=2$ into $2$, $3$ and $4$ secondary breathers with zero intrinsic vorticity at $ \mu =0.4, \varepsilon =0.9$ (top), $ \mu =0.4, \varepsilon =1.4$ (middle), and $ \mu =0.2, \varepsilon =3.0$ (bottom), displayed by means of temporal cross sections, $| \psi \left ( 0,0,z,t\right ) |$. Panels in the right column display the respective evolution in the spatial cross section, $| \psi \left ( x,0,z,0\right ) |$. (b) The same as (a) but illustrating the isosurface evolution in 3D.

Download Full Size | PDF

3. Conclusions

In this work, we have proposed and investigated a model for the passively mode-locked laser, based on the graded-index nonlinear multimode fiber. The model is based on the $(3+1)$D cGLE (complex Ginzburg-Landau equation) with the complex cubic-quintic nonlinearity and the transverse two-dimensional trapping potential, which represents the graded-index structure. A novelty in comparison with previously studied models is that stable stationary and dynamical modes are found in the absence of the spatial and temporal filtering in the model.

The systematic numerical analysis reveals a variety of lasing regimes, including fundamental $(3+1)$D solitons and breathers, as well as stable solitary vortices with winding number $n=1$, that can also exhibit a splitting instability leading to formation of two fundamental solitons. The vortices with $n=2$ are found to be unstable against splitting into a rotating bound state of two unitary vortices. A remarkable feature of the system is the coexistence of between the stable fundamental solitons and vortices with $n=1$. A variety of operation regimes reported here, and a possibility of their implementation in practical devices is expected to assist the development of methods to optimize and tune temporal, spatial and spectral characteristics of short spatiotemporal pulses generated in multimode fiber lasers.

Funding

Thailand Research Fund (BRG6080017); Russian Foundation for Basic Research (17-02-00081); Israel Science Foundation (1287/17).

References

1. P. Grelu and N. Akhmediev, “Dissipative solitons for mode-locked lasers,” Nat. Photonics 6(2), 84–92 (2012). [CrossRef]  

2. W. Fu, L. G. Wright, P. Sidorenko, S. Backus, and F. W. Wise, “Several new directions for ultrafast fiber lasers [Invited],” Opt. Express 26(8), 9432–9463 (2018). [CrossRef]  

3. L. G. Wright, Z. M. Ziegler, P. M. Lushnikov, and Z. Zhu, “Multimode Nonlinear Fiber Optics: Massively Parallel Numerical Solver, Tutorial, and Outlook,” IEEE J. Sel. Top. Quantum Electron. 24(3), 1–16 (2018). [CrossRef]  

4. C. Jauregui, J. Limpert, and A. Tünnermann, “High-power fibre lasers,” Nat. Photonics 7(11), 861–867 (2013). [CrossRef]  

5. D. Richardson, J. Fini, and L. Nelson, “Space-division multiplexing in optical fibres,” Nat. Photonics 7(5), 354–362 (2013). [CrossRef]  

6. Z. Liu, L. G. Wright, D. N. Christodoulides, and F. W. Wise, “Kerr self-cleaning of femtosecond-pulsed beams in graded-index multimode fiber,” Opt. Lett. 41(16), 3675–3678 (2016). [CrossRef]  

7. K. Krupa, A. Tonello, B. M. Shalaby, M. Fabert, A. Barthé lémy, G. Millot, S. Wabnitz, and V. Couderc, “Spatial beam self-cleaning in multimode fibres,” Nat. Photonics 11(4), 237–241 (2017). [CrossRef]  

8. A. S. Ahsan and G. P. Agrawal, “Graded-index solitons in multimode fibers,” Opt. Lett. 43(14), 3345–3348 (2018). [CrossRef]  

9. L. G. Wright, W. H. Renninger, D. N. Christodoulides, and F. W. Wise, “Spatiotemporal dynamics of multimode optical solitons,” Opt. Express 23(3), 3492–3506 (2015). [CrossRef]  

10. O. V. Shtyrina, M. P. Fedoruk, Y. S. Kivshar, and S. K. Turitsyn, “Coexistence of collapse and stable spatiotemporal solitons in multimode fibers,” Phys. Rev. A 97(1), 013841 (2018). [CrossRef]  

11. L. G. Wright, D. N. Christodoulides, and F. W. Wise, “Controllable spatiotemporal nonlinear effects in multimode fibres,” Nat. Photonics 9(5), 306–310 (2015). [CrossRef]  

12. Z. Zhu, L. G. Wright, D. N. Christodoulides, and F. W. Wise, “Observation of multimode solitons in few-mode fiber,” Opt. Lett. 41(20), 4819–4822 (2016). [CrossRef]  

13. W. H. Renninger and F. W. Wise, “Optical solitons in graded-index multimode fibres,” Nat. Commun. 4(1), 1719 (2013). [CrossRef]  

14. K. Krupa, A. Tonello, A. Barthélémy, V. Couderc, B. M. Shalaby, A. Bendahmane, G. Millot, and S. Wabnitz, “Observation of geometric parametric instability induced by the periodic spatial self-imaging of multimode waves,” Phys. Rev. Lett. 116(18), 183901 (2016). [CrossRef]  

15. H. Qin, X. Xiao, P. Wang, and C. Yang, “Observation of soliton molecules in a spatiotemporal mode-locked multimode fiber laser,” Opt. Lett. 43(9), 1982–1985 (2018). [CrossRef]  

16. Y. Ding, X. Xiao, P. Wang, and C. Yang, “Multiple-soliton in spatiotemporal mode-locked multimode fiber lasers,” Opt. Express 27(8), 11435–11446 (2019). [CrossRef]  

17. T. Chen, Q. Zhang, Y. Zhang, and X. Li, “All-fiber passively mode-locked laser using nonlinear multimode interference of step-index multimode fiber,” Photonics Res. 6(11), 1033–1039 (2018). [CrossRef]  

18. D. Mao, M. Li, Z. He, X. Cui, H. Lu, W. Zhang, H. Zhang, and J. Zhao, “Optical vortex fiber laser based on modulation of transverse modes in two mode fiber,” APL Photonics 4(6), 060801 (2019). [CrossRef]  

19. I. S. Chekhovskoy, A. M. Rubenchik, O. V. Shtyrina, and M. P. Fedoruk, “Nonlinear combining and compression in multicore fibers,” Phys. Rev. A 94(4), 043848 (2016). [CrossRef]  

20. S. Minardi, F. Eilenberger, and Y. V. Kartashov, “Three-Dimensional Light Bullets in Arrays of Waveguides,” Phys. Rev. Lett. 105(26), 263901 (2010). [CrossRef]  

21. F. Eilenberger, K. Prater, S. Minardi, R. Geiss, U. Röpke, J. Kobelke, K. Schuster, H. Bartelt, S. Nolte, A. Tünnermann, and T. Pertsch, “Observation of discrete, vortex light bullets,” Phys. Rev. X 3(4), 041031 (2013). [CrossRef]  

22. H. G. Winful and D. T. Walton, “Passive mode locking through nonlinear coupling in a dual-core fiber laser,” Opt. Lett. 17(23), 1688–1690 (1992). [CrossRef]  

23. J. Proctor and J. N. Kutz, “Nonlinear mode-coupling for passive mode-locking: application of waveguide arrays, dual-core fibers, and/or fiber arrays,” Opt. Express 13(22), 8933–8950 (2005). [CrossRef]  

24. T. F. S. Büttner, D. D. Hudson, E. C. Mägi, A. C. Bedoya, T. Taunay, and B. J. Eggleton, “Multicore, tapered optical fiber for nonlinear pulse reshaping and saturable absorption,” Opt. Lett. 37(13), 2469–2471 (2012). [CrossRef]  

25. D. Auston, “Transverse mode locking,” IEEE J. Quantum Electron. 4(6), 420–422 (1968). [CrossRef]  

26. P. W. Smith, “Simultaneous phase-locking of longitudinal and transverse laser modes,” Appl. Phys. Lett. 13(7), 235–237 (1968). [CrossRef]  

27. A. M. Dunlop, W. J. Firth, and E. M. Wright, “Master equation for spatiotemporal beam propagation and Kerr lens mode-locking,” Opt. Commun. 138(1-3), 211–226 (1997). [CrossRef]  

28. N. Akhmediev and A. Ankiewicz, eds. Dissipative solitons, Vol. 661 of Lecture Notes in Physics (Springer, 2005).

29. S. Raghavan and G. P. Agrawal, “Spatiotemporal solitons in inhomogeneous nonlinear media,” Opt. Commun. 180(4-6), 377–382 (2000). [CrossRef]  

30. N. Akhmediev, J. M. Soto-Crespo, and P. Grelu, “Spatiotemporal optical solitons in nonlinear dissipative media: From stationary light bullets to pulsating complexes,” Chaos 17(3), 037112 (2007). [CrossRef]  

31. N. B. Aleksić, V. Skarka, D. V. Timotijević, and D. Gauthier, “Self-stabilized spatiotemporal dynamics of dissipative light bullets generated from inputs without spherical symmetry in three-dimensional Ginzburg-Landau systems,” Phys. Rev. A 75(6), 061802 (2007). [CrossRef]  

32. S. Chen, “Analytical spinless light-bullet solutions as attractive fixed points in the three-dimensional cubic-quintic complex Ginzburg-Landau equation,” Phys. Rev. A 86(3), 033829 (2012). [CrossRef]  

33. D. Mihalache, D. Mazilu, F. Lederer, and H. Leblond, “Stability limits for three-dimensional vortex solitons in the Ginzburg-Landau equation with the cubic-quintic nonlinearity,” Phys. Rev. A 76(4), 045803 (2007). [CrossRef]  

34. D. Mihalache, D. Mazilu, F. Lederer, and H. Leblond, “Collisions between coaxial vortex solitons in the three-dimensional cubic-quintic complex Ginzburg-Landau equation,” Phys. Rev. A 77(3), 033817 (2008). [CrossRef]  

35. B. Liu, Y.-F. Liu, and X.-D. He, “Impact of phase on collision between vortex solitons in three-dimensional cubic-quintic complex Ginzburg-Landau equation,” Opt. Express 22(21), 26203–26211 (2014). [CrossRef]  

36. N. A. Veretenov, S. V. Fedorov, and N. N. Rosanov, “Topological Vortex and Knotted Dissipative Optical 3D Solitons Generated by 2D Vortex Solitons,” Phys. Rev. Lett. 119(26), 263901 (2017). [CrossRef]  

37. B. A. Malomed, “(INVITED)Vortex solitons: Old results and new perspectives,” Phys. D (Amsterdam, Neth.) 399, 108–137 (2019). [CrossRef]  

38. T. Mayteevarunyoo, B. A. Malomed, and D. V. Skryabin, “One- and two-dimensional modes in the complex Ginzburg-Landau equation with a trapping potential,” Opt. Express 26(7), 8849–8865 (2018). [CrossRef]  

39. T. Mayteevarunyoo, B. A. Malomed, and D. Skryabin, “Vortex modes supported by spin-orbit coupling in a laser with saturable absorption,” New J. Phys. 20(11), 113019 (2018). [CrossRef]  

40. J. Yang, Nonlinear Waves in Integrable and Nonintegrable Systems (SIAM: Philadelphia, 2010).

41. L. G. Wright, D. N. Christodoulides, and F. W. Wise, “Spatiotemporal mode-locking in multimode fiber lasers,” Science 358(6359), 94–97 (2017). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (13)

Fig. 1.
Fig. 1. (a) Numerically computed propagation constant $q$ and integral power $P$ [see Eq. (5)] of fundamental three-dimensional solitons versus $ \varepsilon$, for fixed $ \mu =1$ in Eq. (3). Red, blue and black segments correspond to stable fundamental solitons, moving breathers, and eventually destructed modes, respectively. (b) The stability chart for the stationary fundamental solitons and modes into which unstable solitons are spontaneously transformed. Indicated inside the stability area of solitons is the bistability region, in which solitary vortices with winding numbers are stable too, cf. Figure 5(b).
Fig. 2.
Fig. 2. (a) The evolution of a stable fundamental soliton, for $ \mu =1, \varepsilon =0.5$, with $q=-0.8635,P=6.9898$, is displayed by means of isosurfaces of $\left \vert \psi \left ( x,y,z,t\right ) \right \vert$. (b) Left and right panels display spatial and temporal cross sections of the soliton’s shape, $| \psi (x,0,z,0)|$ and $| \psi (0,0,z,t)|$, respectively.
Fig. 3.
Fig. 3. (a) The evolution of a robust quiescent breather, into which an unstable fundamental soliton transforms for $ \mu =0.3, \varepsilon =2.8$, which corresponding to $q=-2.0446,P=5.3081$, shown by means of the isosurface of $\left \vert \psi \left ( x,y,z,t\right ) \right \vert$. (b) Left and right panels show the evolution in the spatial and temporal of cross sections, $| \psi (x,0,z,0)|$ and $| \psi (0,0,z,t)|$, respectively.
Fig. 4.
Fig. 4. A robust moving breather for $ \mu =1.0, \varepsilon =1.8$, into which the respective unstable soliton, with $q=-1.7927,P=4.5920$, spontaneously transforms. (a) The evolution in the spatial and temporal cross sections, $| \psi (x,0,z,0)|$ (left) and $| \psi (0,0,z,t)|$ (right). (b) The spatiotemporal dynamics with temporal motion.
Fig. 5.
Fig. 5. (a) The numerically found values of $q$ (top panel) and $P$ (bottom panel) for stationary vortices with winding number $n=1$ versus $ \varepsilon$ for fixed $ \mu =1.0$. Red, blue, magenta, yellow and black segments correspond to stable vortices, splitting into fundamental solitons, chaotically oscillating modes, moving vortex breathers, and destruction of unstable vortices, respectively. (b) The stability chart for vortices with $n=1$ in the plane of $\left ( \mu , \varepsilon \right )$. The central shaded area in (b) is populated by chaotically oscillating localized modes, see an example in Fig. 8. Note that the stability area of vortices is completely covered by the region in which the fundamental solitons are stable [cf. Fig. 1(b)], i.e., it is a bistability area.
Fig. 6.
Fig. 6. (a) The isosurface evolution of a stable vortex found at $ \mu =1.0, \varepsilon =0.2$, which corresponding to $q=-2.3130,P=15.9142$. (b) Its amplitude and phase patterns (c) Left and right panels display the evolution of the spatial and temporal cross sections $| \psi (x,0,z,0)|$ and $| \psi (0,0,z,t)|$, respectively.
Fig. 7.
Fig. 7. The spontaneous splitting of an unstable vortex with $n=1$ into two fundamental solitons, at $ \mu =1.0, \varepsilon =0.6$. (a) The evolution of $| \psi |$ is shown in the $t$- and $x$-cross sections (top and bottom panels, respectively). (b) The temporal shape of the moving fundamental solitons, produced by the splitting, at $z=400$ (the blue line), fitted to the standard soliton’s shape, $\mathrm {sech}(t)$.
Fig. 8.
Fig. 8. (a) The evolution of $| \psi \left ( 0,0,z,t\right ) |$ and $| \psi \left ( x,0,z,0\right ) |$ in the $t$- and $x$- cross-sections, illustrating the spontaneous transformation of an unstable stationary vortex with $n=1$ into a chaotically oscillating mode at $ \mu =1.0, \varepsilon =0.9$. (b) The respective amplitude and phase patterns at $z=1000$ (eventual values of the amplitude are small, as a result of the evolution).
Fig. 9.
Fig. 9. A moving breather, into which an unstable stationary vortex, with $q=-2.9276,P=9.9546$, is spontaneously transformed (keeping its vorticity) at $ \mu =1.0, \varepsilon =1.2$. The evolution of cross sections $| \psi \left ( x,0,z,0\right ) |$ and $| \psi \left ( 0,0,z,t\right ) |$ is displayed in the left and right panels, respectively.
Fig. 10.
Fig. 10. The stability chart in the plane of $\left ( \mu , \varepsilon \right )$ for vortices with $n=2$. The shaded area corresponds to the “ wobbling breathers” .
Fig. 11.
Fig. 11. Spontaneous splitting of an unstable vortex with $n=2$ into a stably rotating bound state of two vortices with $n=1$ at $ \mu =0.4$, $ \varepsilon =0.1$ (a) The evolution of the field in the temporal and spatial cross sections, $| \psi (0,0,z,t)|$ and $| \psi (x,0,z,0)|$. (b) Amplitude and phase patterns of $ \psi (x,y)$ in the plane of $t=0$ and (bottom) the temporal profile at the central point, $x=y=0$, plotted at $z=500$ (c) The top, left and right panels zoom in on the evolution of the absolute value of the field, and spatial and temporal cross sections, respectively, in the interval of $300<z<320$. The established evolution (rotation) period is $T_{z}=2.2$.
Fig. 12.
Fig. 12. The “ wobbling breather” at $ \mu =0.4$, $ \varepsilon =0.6$ (a) The evolution of the field in the temporal and spatial cross sections, $| \psi (0,0,z,t)|$ and $| \psi (x,0,z,0)|$. (b) The top, left and right panels zoom in on the evolution of the absolute value of the field, and spatial and temporal cross sections in the interval of $z=200-220$, respectively. The established evolution period is $T_{z}=4.5$.
Fig. 13.
Fig. 13. (a) Splitting of an unstable breather with $n=2$ into $2$, $3$ and $4$ secondary breathers with zero intrinsic vorticity at $ \mu =0.4, \varepsilon =0.9$ (top), $ \mu =0.4, \varepsilon =1.4$ (middle), and $ \mu =0.2, \varepsilon =3.0$ (bottom), displayed by means of temporal cross sections, $| \psi \left ( 0,0,z,t\right ) |$. Panels in the right column display the respective evolution in the spatial cross section, $| \psi \left ( x,0,z,0\right ) |$. (b) The same as (a) but illustrating the isosurface evolution in 3D.

Equations (5)

Equations on this page are rendered with MathJax. Learn more.

i ψ z + L 2 τ 2 k ( 1 i β ) 2 ψ t 2 + L 2 k 0 w 2 ( 2 x 2 + 2 y 2 ) ψ k 0 g w 2 L 2 ( x 2 + y 2 ) ψ + i L ε ψ + L ( η i α ) | ψ | 2 ψ + L ( ν + i μ ) | ψ | 4 ψ = 0.
i ψ z + 1 2 ( 1 i β ) 2 ψ t 2 + 1 2 ( 2 x 2 + 2 y 2 ) ψ ( x 2 + y 2 i ε ) ψ + ( η i α ) | ψ | 2 ψ + ( ν + i μ ) | ψ | 4 ψ = 0.
α = A α t h r ,   α t h r ( 2 ε μ )
q ϕ + 1 2 ( 2 t 2 + 2 x 2 + 2 y 2 ) ϕ ( x 2 + y 2 i ε ) ϕ + ( η i α ) | ϕ | 2 ϕ + i μ | ϕ | 4 ϕ = 0.
P = | ψ ( x , y , t ) | 2 d x d y d t .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.