Breather and Rogue Wave solutions of a Generalized Nonlinear Schrödinger Equation

Breather and Rogue Wave solutions of a Generalized Nonlinear Schrödinger Equation


In this paper, using the Darboux transformation, we demonstrate the generation of first-order breather and higher-order rogue waves from a generalized nonlinear Schrödinger equation with several higher-order nonlinear effects representing femtosecond pulse propagation through nonlinear silica fiber. The same nonlinear evolution equation can also describes the soliton-type nonlinear excitations in classical Heisenberg spin chain. Such solutions have a parameter , denoting the strength of the higher-order effects. From the numerical plots of the rational solutions, the compression effects of the breather and rogue waves produced by are discussed in detail.

05.45.Yv, 42.65.Tg, 42.65.Sf, 02.30.Ik

I Introduction

It is well known that one of the most challenging aspects of modern science and technology is the nonlinear nature of the system, which is considered to be fundamental to the understanding of many natural phenomena. In recent years, nonlinear science has emerged as a powerful subject for explaining the mysteries of the challenging nature. Nonlinearity is a fascinating occurrence of nature whose importance has been well appreciated for many years, in the context of large amplitude waves or high-intensity laser pulses observed in various fields ranging from fluids to solid state, chemical, biological, nonlinear optical, and geological systems (1); (2); (3); (4); (5); (6); (7); (8); (9). This fascinating subject has branched out in almost all areas of science, and its applications are percolating through the whole of science. In general, nonlinear phenomena are often modeled by nonlinear evolution equations exhibiting a wide range of high complexities in terms of different linear and nonlinear effects. In recent years, the advent of high-speed computers, many advanced mathematical software, and development of many sophisticated and systematic analytical methods in the study of the nonlinear phenomena and also supported by many experiments have encouraged both theoretical and experimental research. In the past few decades, nonlinear science has experienced an explosive growth by the invention of several exciting and fascinating new concepts, such as solitons, dispersion-managed solitons, dromions, rogue waves, similaritons, supercontinuum generation, complete integrability, fractals, chaos, etc. (1); (2); (3); (4); (5); (6); (7); (8); (9); (10); (11). Many of the completely integrable nonlinear partial differential equations (NPDEs) admit one of the most striking aspects of nonlinear phenomena called soliton, which describe soliton as a universal character, and they are of great mathematical as well as physical interest, too. The study of the solitons and other related issues of the construction of the solutions to a wide class of NPDEs have become one of the most exciting and extremely active areas of research in science and technology for many years.

In addition to several developments in soliton theory, recent developments in modulational instability (MI) have also been widely used to explain why experiments involving white coherent light supercontinuum generation (SCG), admit a triangular spectrum. Such universal triangular spectra can be well described by the analytical expressions for the spectra of Akhmediev breather solutions at the point of extreme compression. In the context of the NLS equation, Peregrine already in Ref. (12) had identified the role of MI in the formation of patterns resembling freak waves or rogue wave (RW); these theoretical results were later supported by several experiments. Rogue waves in the ocean are localized large amplitude waves on a rough background, which have two remarkable characteristics: (1) “appear from nowhere and disappear without a trace” (13), (2) exhibit one dominant peak. RWs have recently appeared in several areas of science. Particularly in photonic crystal fibers, RWs have been well established in connection with SCG (14). This actually has stimulated research for RWs in other physical systems and has paved the way for many important applications, including the control of RWs by means of SCG (15); (16), as well as studies in superfluid Helium (17), Bose-Einstein condensates (18), plasmas (19); (20), microwave (21), capillary phenomena (22), in telecommunication data streams (23), inhomogeneous media (24), water experiments (25); (26), and so on. More recently, Kibler et al. (27), using their elegant experimental apparatus in optical fibers, were able to generate femtosecond pulses with strong temporal and spatial localization and near-ideal temporal Peregrine soliton characteristics.

In the recent past, several equations have been shown to admit the rogue wave solutions. For example, in addition to the NLS equation, the Hirota equation (28); (29); (30), the first-type derivative NLS equation (31), the third-type DNLS equation (32), the Fokas-Leneels equation (33), the NLS-MB equations (34), the Hirota Maxwell-Bloch(MB) equations (35), the Sasa-Satsuma equation (36), the discrete Ablowitz-Ladik and Hirota equation (37), the two-component NLS equations (38); (39); (40), the three-components NLS equations (41), the variable coefficient NLS (42); (43); (44); (45); (46), the variable coefficient derivative NLS (47), and the variable coefficient higher-order NLS (VCHNLS) (48), the rogue waves in dissipative systems (49) are a few of the nonlinear evolution equations that admit RWs. From the above studies, it is clear that one of the possible generating mechanisms (50) for the higher-order RW is the interaction of the multiple breathers possessing the same and the very particular frequency of the underlying equation.

In recent years, there has been a considerable interest in the study of the nonlinear excitations of the spin chains with competing bilinear and biquadratic interactions. In particular, the complete integrability and nonlinear excitations of spin chains with spin magnitude has been established if suitable polynomials in () are added to the original bilinear Heisenberg spin Hamiltonian. In this connection and also from the mathematical point of view, it is of interest to study the influence of the biquadratic interactions on the integrability of the Heisenberg bilinear spin chain in the classical limit as well. Considering the above points, one of the authors of this paper, has investigated the integrability aspects of a classical one-dimensional isotropic biquadratic Heisenberg spin chain in its continuum limit up to order in the lattice parameter through a classical differential geometric approach (51); (52); (53) and investigated the soliton and integrability aspects of the corresponding generalized nonlinear Schrödinger Equation (GNLSE). This equation is given by (51)


Here, is the complex envelope and denotes the strength of higher-order linear and nonlinear effects. When we consider the propagation of ultra short pulse propagation through optical fiber, i.e, less than 100-fs pulses, it has been shown that higher-order dispersion, self-steepening, self-frequency, and quintic effects should be included in the model. The above equation has been shown to be integrable and admits exact soliton solutions and also gauge equivalent to Heisenberg spin chain equation. Thus, it is an interesting problem to find how these higher-order terms will affect the breather and rogue wave in an associated optical system and spin system by means of changing the value of . It is our prime aim to answer this problem in this paper.

The paper is organized as follows. In Sec. II, the Lax pair and the Darboux transformation (DT) are introduced. In Sec. III, we shall give the first-order breather and it’s limit of infinitely large period. In Sec. IV, higher-order rogue waves are given. The compression effects on the breather and rogue waves produced by the higher-order terms of GNLSE are given in Secs. III and IV, respectively. Section V is devoted to conclusions.

Ii Lax Pair and Darboux transformation

As discussed above, in this section, we would like to recall the Lax pair of GNLSE (51); (52); (53)and to show its Darboux transformation. According to the AKNS formalism, Lax pair for Eq. (1) is written as


with the following matrices:

and . Here,


denotes the eigenfunction of Lax pair Eq. (2) associated with .

The Lax pair of GNLSE provides a basis for the solvability of this equation by means of the Darboux transformation. To construct the -fold DT, it is necessary to introduce eigenfunctions , associated with eigenfunction , and satisfy corresponding reduction condition as we have done for the NLS equation (54); (55). Furthermore, a similar -fold DT determinant representation derived for the NLS equation as given in (54). For example, from the one-fold DT, we get



Iii The first-order breather and its limit

In this section, we first solve the eigenfunctions associated with a periodic seed , and then use it to get a first-order breather by using the determinant representation of one-fold DT in Eq. (3). Further, this breather implies a first-order rogue wave in the limit of infinitely large period. These two solutions have explicitly such that we can use it to study the effects of breather and rogue waves affected by the higher-order terms.

Considering a periodic solution in the following form:


with . By the method of separation of variables and the superposition principle, we have the following eigenfunction associated with :


with, ,, ,

Using the one-fold DT, a first-order breather is constructed in the form


with , . This is a periodic traveling wave. The coefficient can affect the period of the breather through .

It is trivial to find

which is the height of peaks of this breather. Obviously, the height is independent of and . This does not mean that cannot affect the properties of the breather. In fact, we can see from Eq. (7) that actually controls the period of the breather. This observation can be clearly seen in Fig. 1: the number of peaks on same time interval is increasing when goes up from to with a constant gap (we skip the figure for ) . These pictures clearly show that the resulting breather is compressed by the higher-order effects due to the presence of than case. In addition to the above, when the value of increases, the number of peaks also increases. We use a short interval in Fig. 1(a) and 1(b) to avoid too many peaks in it.

Now we can consider what will happen in a breather when its period goes to infinity. According to the explicit expression in Eq. (7) of the first-order breather, this limit can be realized by setting ; i.e., . For simplicity, on setting , the limit of the breather solution is obtained as


with . This is nothing but a first-order rogue wave possessing asymptotic height when and go to infinity. Further, we find , which denotes the height of a first-order RW. We can see from that, like in the case of breather compression, is also responsible for compression effect of RW in the time direction, which is clearly seen in Fig. 2 with , respectively.

Iv Higher-order Rogue Waves

The limit method in Eq. (8) is not applicable for the higher-order breather when . We can overcome this problem by using the coefficient of the Taylor expansion in the determinant representation of a higher-order breather (50); (55). Similar to the case of NLS (50), the first-order rogue wave of GNLSE is given by



Note that reduces to the when .

What follows is the second-order rogue wave given by the Taylor expansion when as the case of NLS (50). There are two patterns for the second-order RW. The first one is called the fundamental pattern possessing a highest peak surrounded by four small equal peaks in two sides. Setting , then the Taylor expansion in the determinant of provides



and are defined in Eq. (8). By comparing two cases with and in Figure 3, the compression effect in direction is shown clearly. The second is a triangular pattern, which consists of three equal peaks. Setting ,, then an explicit formula of this pattern is



and X,T are defined in Eq. (8). Figure 4 is plotted for the to show its compression effect. Because of the explicit appearance of in , , , and one is not be able to derive from the corresponding RWs of the NLSE by a mere scalar transformation of .

Next, we construct the third-order RW of GNLSE by substituting present in (50). There exists three patterns: a fundamental pattern with , a triangular pattern with , and a circular pattern with . The explicit form of the fundamental pattern with the third-order RW is given in the Appendix, and, for brevity, the very lengthy forms of the other two cases are deleted. Note that includes the -dependence explicitly. This fact shows that the second-order RW solution cannot be obtained from the corresponding solution of the NLS equation by a scalar transformation of . Furthermore, to show the compression effect on the RWs, Figs. (5-7) are plotted for different parametric choices.

V Conclusions

In this paper, we considered the integrable version of the generalized nonlinear Schrödinger equation with several higher-order nonlinear terms, which describes ultra short pulse propagation through nonlinear silica fiber and soliton-type nonlinear excitations in classical Heisenberg spin chain. Using Daurboux transformation and periodic seed solutions, we have constructed the first-order breather solution and also discussed the behavior of these solutions with an infinitely large period. Finally, we have also constructed the first-order, second-order and third-order rogue wave solutions by the Taylor expansion. All of these solutions have parameter denoting the contribution of higher-order nonlinear terms. The compressed effects of these solutions are discussed through numerical plots by increasing the value of . This new phenomenon of the rogue wave is useful for us to observe or analyze its evolution in some complicated physical system. Another advantage of our results of this paper is that, as GNLSE is equivalent to spin chain, the rogue wave nature of spin systems can also be explained through suitable geometrical and gauge equivalence methods. In addition, as higher-order linear and nonlinear effects in optical fibers are playing key roles in explaining the generation and propagation of ultra short pulse through silica wave guides, we hope that our results with all these higher-order effects can be observed in real experiments in the near future.


This work is supported by the NSF of China under Grants No. 10971109 and No. 11271210, K. C. Wong Magna Fund in Ningbo University. J.H. is also supported by the Natural Science Foundation of Ningbo under Grant No. 2011A610179. J.H. thanks A. S. Fokas for his support during his visit at Cambridge. K.P. thanks the DST, DAE-BRN, and CSIR, Government of India, for the financial support through major projects. L.W. is also supported by the Natural Science Foundation of China, under Grants No. 11074136 and No. 11101230, and the Natural Science Foundation of Zhejiang province under Grant No. 2011R09025-06.


Figure 1: (Color online) The dynamical evolution of the first-order breather on the () plane. When the value of increases, the number of peaks on same interval of time also increases.
Figure 2: (Color online) The dynamical evolution of the first-order rogue wave on the () plane. For larger values of , it is clear that the compression effects in direction are quite high.
Figure 3: (Color online) The dynamical evolution of the second-order rogue wave on the () plane. Comparing (a) and (b) with (c) and (d) indicates effective high compression in the direction.
Figure 4: (Color online) The dynamical evolution of the second-order rogue wave on the () plane. It is shown from (b),(c), and (d) that rogue wave compression increases as the value of increases.
Figure 5: (Color online) The dynamical evolution of the third-order rogue wave on the () plane. By comparison with (a) and (b), (c) and (d) are highly compressed.
Figure 6: (Color online)The dynamical evolution of the third-order rogue wave on the () plane. It is shown from (b), (c), and (d) that rogue wave is compressed more while increasing the value of .
Figure 7: (Color online) The dynamical evolution of the third-order rogue wave on the () plane. It is shown from (b), (c), and (d) that rogue wave is changing its shape and also compression increases by increasing the value of .


  1. G. P. Agrawal, Nonlinear Fiber Optics, (Academic Press, San Diego 2006).
  2. A. Hasegawa, M. Matsumoto, Optical solitons in fibers, (Springer-Verlag, Berlin, 2003).
  3. L.F. Mollenauer, J.P. Gordon, Solitons in Optical Fibers: Fundamentals And Applications, (Academic Press, New York, 2006).
  4. K. Porsezian and V. C. Kuriakose, Optical Solitons: Theory and Experiment, Lecture Notes in Physics, Vol. 613, (Springer-Verlag, Berlin, February 2003).
  5. K. Porsezian and V. C. Kuriakose, Solitons in Nonlinear Optics: Advances and Applications, European Journal of Physics-Special Topics, (Springer-Verlag, Berlin 2009).
  6. J. R. Taylor, Optical Solitons: Theory and Experiment, (Cambridge University Press, Cambridge, 1992).
  7. F. Abdullaev, S.Darmanyan, P.Khabibullaev, Optical Solitons, (Springer Series in Nonlinear Dynamics, Berlin, 1993).
  8. F. Abdullaev, Theory of Solitons in Inhomogeneous Media, (John Wiley, New York, 1994).
  9. Y. S. Kivshar, G.P. Agrawal, Optical Solitons: From Fibres to Photonic Crystals, (Academic Press, San Diego, 2003).
  10. G. P. Agrawal, Applications of Nonlinear Fibre Optics, 2nd ed., (Academic Press, New York, 2006).
  11. N. N. Akhmediev and A. Ankiewicz, Solitons, (Chapman Hall, London, 1997).
  12. D. H. Peregrine, J. Austral. Math. Soc. B 25, 16 (1983).
  13. N. Akhmediev, A. Ankiewicz, and M. Taki, Phys. Lett. A 373, 675 (2009).
  14. D. R. Solli, C. Ropers, P.Koonath, and B. Jalali, Nature (London) 450, 1054 (2007).
  15. D. R. Solli, C. Ropers, and B. Jalali,Phys. Rev. Lett. 101, 233902 (2008).
  16. J. M. Dudley, G. Genty, and S. Coen, Rev. Mod. Phys. 78, 1135 (2006); J. M. Dudley, G. Genty, F. Dias, B. Kibler, and N. Akhmediev, Opt. Express 17, 21497 (2009).
  17. A. N. Ganshin, V. B. Efimov, G. V. Kolmakov, L. P. Mezhov-Deglin, and P. V. E. McClintock, Phys. Rev. Lett. 101, 065303 (2008).
  18. Yu. V. Bludov, V. V. Konotop, and N. Akhmediev, Phys. Rev. A. 80, 033610 (2009).
  19. M. S. Ruderman, Eur. Phys. J. Special Topics 185, 57 (2010).
  20. W. M. Moslem, P. K. Shukla, and B. Eliasson, Euro. Phys. Lett. 96, 25002 (2011).
  21. R. Höhmann, U. Kuh, H.-J. Stockmann, L. Kaplan, and E. J. Heller, Phys. Rev. Lett. 104, 093901 (2010).
  22. M. Shats, H. Punzmann, andH.Xia,Phys.Rev. Lett. 104, 104503 (2010).
  23. A. I. Chervanyov, Phys. Rev. E 83, 061801(R) (2011).
  24. F. T. Arecchi, U. Bortolozzo, A.Montina, and S. Residori, Phys. Rev. Lett. 106, 153901 (2011).
  25. A. Chabchoub, N. P. Hoffmann, and N. Akhmediev, Phys. Rev. Lett. 106, 204502 (2011).
  26. A. Chabchoub, N. P. Hoffmann, and N. Akhmediev, J. Geophys. Res. 117, C00J02 (2012).
  27. B. Kibler, J. Fatome, C. Finot, G. Millot, F. Dias, G. Genty, N. Akhmediev, and J. M. Dudley, Nat. Phys. 6, 790 (2010).
  28. A. Ankiewicz, J. M. Soto-Crespo, and N. Akhmediev, Phys. Rev. E. 81, 046602 (2010).
  29. Y. S. Tao and J. S. He, Phys. Rev. E. 85, 026601 (2012).
  30. G. G. Yang, L. Li, and S. T. Jia, Phys. Rev. E 85, 046608 (2012).
  31. S. W. Xu, J. S. He, and L. H. Wang, J. Phys. A 44, 305203 (2011).
  32. S. W. Xu and J. S. He, J. Math. Phys. 53, 063507 (2012).
  33. J. S. He, S. W. Xu, and K. Porsezian, J. Phys. Soc. Jpn. 81, 124007 (2012).
  34. J. S. He, S. W. Xu, and K. Porsezian, J. Phys. Soc. Jpn. 81, 033002 (2012).
  35. C. Z. Li, J. S. He, and K. Porsezian, Phys. Rev. E 87, 012913 (2013).
  36. U. Bandelow and N. Akhmediev, Phys. Rev. E 86, 026606 (2012).
  37. A. Ankiewicz, N. Akhmediev, and J. M. Soto-Crespo, Phys. Rev. E 82, 026602 (2010).
  38. Yu. V. Bludov, V. V. Konotop, and N. Akhmediev, Eur. J. Phys. Special Topics 185, 169 (2010).
  39. B. L. Guo and L. M. Lin, Chin. Phys. Lett. 28, 110202 (2011).
  40. D. Meschede, F. Steglich, W. Felsch, H. Maletta, and W. Zinn, Phys. Rev. Lett. 109, 044102 (2012).
  41. Z. Y. Qin and M. Gu, Phys. Rev. E 86, 036601 (2012).
  42. A. Ankiewicz, N. Devine, and N. Akhmediev, Phys. Lett. A 373, 3997 (2009).
  43. M. Taki, A. Mussot, A. Kudlinski, E. Louvergneaux, M. Kolobov, and M. Douay, Phys. Lett. A 374, 691 (2010).
  44. Z. Y. Yan, Phys. Lett. A. 374, 672 (2010).
  45. Y. Y. Wang, J. S. He, and Y. S. Li, Commun. Theor. Phys. 56, 995 (2011).
  46. L.Wen, L. Li, Z. D. Li, S.W. Song, X. F. Zhang, andW.M. Liu, Eur. Phys. J. D 64, 473 (2011).
  47. S. W. Xu, J. S. He, and L. H. Wang, Europhys. Lett. 97, 30007(2012).
  48. C. Q. Dai, G. Q. Zhou, and J. F. Zhang, Phys. Rev. E 85, 016603 (2012).
  49. J. M. Soto-Crespo, Ph. Grelu, and N. Akhmediev, Phys. Rev. E 84, 016604 (2011).
  50. J. S. He, H. R. Zhang, L. H. Wang, K. Porsezian, and A. S. Fokas, Phys. Rev. E 87, 052914 (2013).
  51. M. Lakshmanan, K. Porsezian, and M. Daniel, Phys. Lett. A 133, 483 (1988).
  52. K. Porsezian, M. Daniel, and M. Lakshmanan, J. Math. Phys. 33, 1807 (1992).
  53. K. Porsezian, Phys. Rev. E 55, 3785 (1997).
  54. V. B. Matveev and M. A. Salle, Darboux Transformations and Solitons (Springer, Berlin, 1991).
  55. J. S. He, L. Zhang, Y. Cheng, and Y. S. Li, Sci. China A 12, 1867 (2006).