Anti-dark and Mexican-hat solitons in the Sasa-Satsuma equation on the continuous wave background

Anti-dark and Mexican-hat solitons in the Sasa-Satsuma equation on the continuous wave background

Tao Xu, Min Li and Lu Li
1. College of Science, China University of Petroleum, Beijing 102249, China,
2. Department of Mathematics and Physics, North China Electric
Power University, Beijing 102206, China.
3. Institute of Theoretical Physics, Shanxi University, Taiyuan, Shanxi 030006, China.

In this letter, via the Darboux transformation method we construct new analytic soliton solutions for the Sasa-Satsuma equation which describes the femtosecond pulses propagation in a monomode fiber. We reveal that two different types of femtosecond solitons, i.e., the anti-dark (AD) and Mexican-hat (MH) solitons, can form on a continuous wave (CW) background, and numerically study their stability under small initial perturbations. Different from the common bright and dark solitons, the AD and MH solitons can exhibit both the resonant and elastic interactions, as well as various partially/completely inelastic interactions which are composed of such two fundamental interactions. In addition, we find that the energy exchange between some interacting soliton and the CW background may lead to one AD soliton changing into an MH one, or one MH soliton into an AD one.

PACS numbers: 05.45.Yv; 42.65.Tg; 42.81.Dp

Introduction. — As localized wave packets formed by the balance between the group-velocity dispersion and self-phase modulation [1], solitons in optical fibers have drawn considerable attention because of their robust nature and potential application in all-optical, long-distance communications [2]. In the picosecond regime, the model governing the propagation of optical solitons in a single-mode fibre is the celebrated nonlinear Schrödinger equation (NLSE) [1]. However, one has to take into account some higher-order linear and nonlinear effects for the ultrashort pulses propagating in high-bit-rate transmission systems [2]. The governing model for the femtosecond pulse propagation is the following higher-order NLSE [3]:


where , is a real small parameter, , and represent the third-order dispersion (TOD), self-steepening (SS, also known as Kerr dispersion) and stimulated Raman scattering (SRS) effects, respectively [2, 3].

With , Eq. (1) has two important integrable versions: (i) the Hirota equation (HE) [4], ; (ii) the Sasa-Satuma equation (SSE) [5], . The SSE usually takes the form [5]


Although there is a fixed relation among the higher-order terms, the SSE is thought to be more fundamental than the HE for applications in optical fibers because the former contains the SRS term [2, 3]. Up to now, many integrable properties of Eq. (2) have been detailed, like the inverse scattering transform scheme [5, 6], bilinear representation [7], Painlevé property [8], conservation laws [9], nonlocal symmetries [10], squared eigenfunctions [11], Bäcklund transformation [12] and Darboux transformation (DT) [13, 14].

The presence of the SRS term enriches the solitonic behavior in Eq. (2[5, 6, 7, 14, 18, 15, 16, 17]. Under the vanishing boundary condition (VBC), Eq. (2) with possesses the common single-hump soliton [5, 15], the double-hump soliton behaving like two in-phase solitons with a fixed separation [5, 6], and the multi-hump breather with the periodically-oscillating structure [6, 7]. Under the non-vanishing boundary condition, Eq. (2) with admits the common dark soliton, and the double-hole dark soliton displaying two symmetric dips with a fixed separation [16, 17]; while Eq. (2) with has the bright-like soliton which is linearly combined of a dark one and a bright one [18]. On the other hand, the soliton interaction behavior underlying in the SSE is far more abundant and complicated than that in the NLSE. Even with the VBC, the shape-changing interactions between soliton and breather have recently been found in Eq. (2) with  [14].

In this letter, we are trying to reveal some novel solitonic phenomena on a continuous wave (CW) background for Eq. (2) with . Via the DT technique [14], we obtain three families of single-soliton solutions which can display two completely-different profiles. The first type is the anti-dark (AD) soliton having the form of a bright soliton on a CW background, i.e., it looks like a dark soliton with reverse sign amplitude [19]. The second type takes the Mexican-hat (MH) shape, that is, one high hump carries two small dips which have a symmetrical distribution with respect to the hump, hence such new type of soliton is called the MH soliton. More importantly, we find that the femtosecond AD and MH solitons admit the resonant interaction, elastic interaction, as well as various partially/completely inelastic interactions which consist of the fundamental resonant and elastic interaction structures. To our knowledge, it is the first time that the coexistence of elastic and resonant soliton interactions has been found in the NLSE-type models. Physically, the resonant interaction of optical waves excited from a CW background can be used to realize a second-harmonic generation in the centro-symmetry optical fiber [20].

N-th iterated soliton solutions via the Darboux transformation. — With the simple CW solution ( and are both real constants) as a seed, we employ the DT-iterated algorithm presented in Ref. [14] to obtain the N-th iterated solution in the form




where , the block matrices , , , and the functions , , are given as


with , , , , , and being nonzero complex constants.

In order to obtain the solitonic structure from solution (3), we require all ’s () be real numbers, that is, and . For convenience of our analysis, we introduce the notations , and (, and use “” to represent the soliton interaction with asymptotic solitons as and ones as .

Figure 1: Evolution of the (a) AD soliton and (b) MH soliton plotted via solution (6), where , , , , , for (a) and for (b). (c) Transverse plots of AD (blue dotted line) and MH (red solid line) solitons at .

Anti-dark and Mexican-hat solitons. — For solution (3) with , we can obtain three families of single-soliton solutions under the reducible cases (). If (i.e., ), the solution can be written as


where , , . In this solution, the first part is a CW solution of Eq. (2), while the second part describes a soliton embedded in the CW background (Note that the denominator has no singularity because ).

The parameter implies that the embedded solution has the same phase as that of the CW solution. In this case, represents an AD soliton which displays the bright soliton profile on the CW pedestal [see Figs. 1 and 1]. The soliton velocity and width are, respectively, given by and , and reaches the maximum when . If , the embedded solution and CW solution have the same phases in the inner region , but their phases are opposite in the outer region or . Hence, the modulus of exhibits that one high hump is symmetrically accompanied with two small dips beneath the CW background, which looks like the MH shape [see Figs. 1 and 1]. The velocity and width of the MH soliton are the same as those of the AD one, but its maximum amplitude drastically increases to at the center of the hump, and drops to zero at the centers of two dips. The generation of the MH soliton could be explained as that the phase oppositeness makes some energy be transferred from the CW background to the embedded solution, and further leads to rising of one hump and sinking of two dips.

For the reducible cases and , we can obtain the other two families of single-soliton solutions as follows:


where , , , , and . Because , either or displays only the AD soliton profile, which is similar to the case in solution (6). Solutions (7) and (8) are also called the combined solitary wave solutions [18]. The combined dark and bright solitons have a constant phase difference or . Such phase difference causes a nonlinear phase shift, for example, the nonlinear phase shift in solution (7) can be given as

Figure 2: Numerical evolution of the AD soliton under the perturbation of a white noise with the maximal value . The initial pulse corresponds to solution (6) at with the parameters as , , , (a) , (b) .

The stability of solitons is a crucial issue for their applications in optical communication lines [2]. It has been shown in Ref. [18] that the solitons described by solutions (7) and (8) enjoy a good stability against finite-amplitude initial perturbations. Here, the numerical simulation is performed to examine the stability of solution (6) by the split-step Fourier method [2]. Fig. 2 shows that with the presence of a white noise, the AD soliton can propagate stably for dispersion lengths along the fiber. Also, we numerically simulate the evolution of the AD soliton when the TOD, SS and SRS terms do not obey the fixed relation in Eq. (2). Fig. 2 illustrates that the AD soliton still keeps its stable shape after propagating dispersion lengths when and the initial pulse is perturbed by a white noise. However, there is some radiation on the background of the MH soliton during the propagation if a white noise is added in the initial pulse. We note that the practical optics telecommunication system is usually dissipative because of the fiber loss/gain [21]. With the inclusion of a linear loss/gain term into Eq. (2), our numerical experiments show that both the AD and MH solitons together with the CW background decay/grow exponentially with the evolution of . Thus, the balance between the energy input and output also plays an important role in maintaining a long-lived optical soliton [22].

Figure 3: (a) Resonant -interaction of three AD solitons with and . (b) Resonant -interaction of two AD solitons and one MH soliton with and . The other parameters are chosen as , , , and .

Resonant and elastic interactions. — If () in solution (3) with , the phase difference (which is neither nor ) between the embedded solution and the CW solution results in that there are three asymptotic solitons appearing on top of the same CW background. The asymptotic expressions of the three solitons as have the same form in Eqs. (6)–(8) except that


Their wave numbers and frequencies can be respectively given as follows: , , , which exactly satisfy the three-soliton resonant conditions and . Associated with and , the solution can, respectively, exhibit the -resonant structure of two solitons merging into one soliton and the -resonant structure of one soliton diverging into two solitons, as shown in Figs. 3 and 3. When , the three resonant solitons all belongs to the AD case; while for , two are still the AD solitons but the other one is of the MH shape. That means that the CW background exchanges its energy with one interacting soliton, and causes such soliton changes its shape after resonant interaction, as shown in Fig. 3.

Figure 4: (a) Shape-preserving elastic -soliton interaction with . (b) Shape-changing elastic -soliton interaction with . The other parameters are chosen as , , , , , , , , and .

On the other hand, one can also obtain the elastic -soliton interactions by implementing the DT two times. Our analysis shows that there are four asymptotic solitons as in solution (3) with under the condition (). The expressions of asymptotic solitons are still of the form in Eqs. (6)–(8) except that there are some difference for the parameters and (; ) (details are omitted for saving the space). Each interacting soliton could be either the AD or MH one, depending on the concrete parametric choice. For example, Fig. 4 illustrates that the AD solitons display the standard elastic interaction, that is, they can completely recover their individual intensities and velocities after an interaction except for the phase shift in their envelops. Note that the phase shift, which corresponds to the instantaneous frequency at pulse peak being nonzero, will result in the relative motion of interacting solitons [23]. Also, the energy exchange may take place between some interacting soliton and the CW background, and result in the shape change of such soliton after interaction, as seen in Fig. 4. However, this kind of soliton interactions are still considered to be elastic in the sense that there is no energy exchange between two different solitons.

Partially and completely inelastic interactions. — For other cases in solution (3) with , one can obtain five different types of inelastic soliton interactions. If there is only one equal to , the solution can exhibit the - and -soliton interactions which are, respectively, associated with and [see Figs. 5 and 5]. In both the two cases, the numbers of interacting solitons as are not equal, but one soliton and one soliton [which are marked by the red arrows in Figs. 5 and 5] have the same velocities and intensities and differ only by the phases of their envelops. Accordingly, the - and -soliton interactions belong to the partially inelastic type. If none of ’s (, ) is equal to , the solution can display the -, - and -soliton interactions which are, respectively, associated with , and [see Figs. 66]. Such three interactions are of the completely inelastic type in the sense that the asymptotic solitons as totally differ from those as in the velocities and intensities.





Figure 5: (a) Inelastic -soliton interaction with , , , and . (b) Inelastic -soliton interaction with , , , and . The other parameters are chosen as , and .
Figure 6: (a) Inelastic -soliton interaction with , , , , , and . (b) Inelastic -soliton interaction with , , , , , and . (c) Inelastic -soliton interaction with , , , , , and . The other parameters are chosen as , , and .

As for the above five inelastic interactions, we find that in the near-field region the asymptotic solitons connect with one another via the “X”- and “Y”-type junctions, which correspond to the elastic and resonant interactions, respectively. For example, there are three resonant interactions and one elastic interaction in Fig. 5, and there appear three resonant interactions in Fig. 5. Therefore, the “X”-type elastic interaction and “Y”-type resonant interaction are two fundamental structures to form various complicated soliton interactions in Eq. (2), where the numbers, velocities and intensities of interacting solitons as are in general not the same.

Conclusion. — In this letter, via the DT method we have constructed new analytic soliton solutions for Eq. (2) which governs the propagation of femtosecond pulses in a monomode fiber with the TOD, SS and SRS effects. We have revealed that two new types of femtosecond solitons (i.e., the AD and MH solitons) occur in Eq. (2) with on a CW background. The numerical experiments have indicated that the AD soliton can propagate stably for a long distance with presence of a small initial perturbation or slight violation of the fixed ratio of parameters in Eq. (2). More importantly, we have obtained that the AD and MH solitons can exhibit both the resonant and elastic interactions. Such two fundamental interactions can generate various complicated structures, in which the numbers, velocities and intensities of interacting solitons are usually not the same before and after interaction. In addition, we have found that some interacting soliton may exchange its energy with the background in the interaction, which results in one AD soliton changing into an MH one, or one MH soliton into an AD one. It should be noted that changing the propagation direction of optical solitons is an important concept for realizing optical switching [24]. Therefore, as a self-induced Y-junction waveguide, the soliton resonant interaction might bring about some applications in all-optical information processing and routing of optical signals [2, 24]. In mathematics, our results will enrich the knowledge of soliton interactions in a (1+1)-dimensional integrable equation with the single field. It is also worthy of being studied to make a finer classification of soliton interactions in Eq. (2) with .

Acknowledgements. — This work has been supported by the Science Foundations of China University of Petroleum, Beijing (Grant No. BJ-2011-04), by the National Natural Science Foundations of China under Grant Nos. 11247267, 11371371, 11426105, 61475198, and by the Fundamental Research Funds of the Central Universities (Project No. 2014QN30).


  • [1] A. Hasegawa and F. Tappert, Appl. Phys. Lett. 23, 142 (1973); 23, 171 (1973).
  • [2] G. P. Agrawal, Nonlinear Fiber Optics (5th edition, Academic, Oxford, 2012).
  • [3] Y. Kodama, J. Stat. Phys. 39, 597 (1985); Y. Kodama and A. Hasegawa, IEEE J. Quantum Electron. 23, 510 (1987).
  • [4] R. Hirota, J. Math. Phys. 14, 805 (1973).
  • [5] N. Sasa and J. Satsuma, J. Phys. Soc. Jpn. 60, 409 (1991).
  • [6] D. Mihalache, L. Torner, F. Moldoveanu, N.-C. Panoiu, and N. Truta, Phys. Rev. E 48, 4699 (1993).
  • [7] C. Gilson, J. Hietarinta, J. Nimmo, and Y. Ohta, Phys. Rev. E 68, 016614 (2003).
  • [8] D. Mihalache, N. Truta, and L. C. Crasovan, Phys. Rev. E 56, 1064 (1997).
  • [9] J. Kim, Q. Han Park, and H. J. Shin, Phys. Rev. E 58, 6746 (1998).
  • [10] A. Sergyeyev and D. Demskoi, J. Math. Phys. 48, 042702 (2007).
  • [11] J. K. Yang and D. J. Kaup, J. Math. Phys. 50, 023504 (2009).
  • [12] O. C. Wright III, Chaos Solitons Fractals 33, 374 (2007).
  • [13] Y. S. Li and W. T. Han, Chin. Ann. Math. Ser. B 22, 171 (2001).
  • [14] T. Xu and X. M. Xu, Phys. Rev. E 87, 032913 (2013); T. Xu, D. H. Wang, M. Li, and H. Liang, Phys. Scr. 89, 075207 (2014).
  • [15] K. Porsezian and K. Nakkeeran, Phys. Rev. Lett. 76, 3955 (1996); M. Gedalin, T. C. Scott, and Y. B. Band, Phys. Rev. Lett. 78, 448 (1997).
  • [16] Y. Jiang and B. Tian, EPL 102, 10010 (2013).
  • [17] Y. Ohta, AIP Conference Proceedings 1212, 114 (2010).
  • [18] Z. H. Li, L. Li, H. P. Tian, and G. S. Zhou, Phys. Rev. Lett. 84, 4096 (2000).
  • [19] Yu. S. Kivshar, Phys. Rev. A 43, 1677 (1991); Yu. S. Kivshar, V. V. Afansjev, and A. W. Snyder, Opt. Commun. 126, 348 (1996).
  • [20] W. N. Cui, G. X. Huang, and B. Hu, Phys. Rev. E 70, 057602 (2004); W. N. Cui and G. X. Huang, Chin. Phys. Lett. 21, 2437 (2004).
  • [21] X. Liu, D. Han, Z. Sun, et al, Sci. Rep. 3, 2718 (2013); Y. Cui and X. Liu, Opt. Express 21, 18969 (2013).
  • [22] X. Liu, Phys. Rev. A 81, 023811 (2010); P. Grelu and N. Akhmediev, Nat. Photon. 6, 84 (2012).
  • [23] X. Liu, Phys. Rev. A 84, 053828 (2011).
  • [24] Yu. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals (Academic, San Diego, 2003).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description