Cite this Article
Bai Cheng-Lin, Cai Yue-Jin, Luo Qing-Long. Breathers and Rogue Waves Derived from an Extended Multi-dimensional N-Coupled Higher-Order Nonlinear Schrödinger Equation in Optical Communication Systems. Communications in Theoretical Physics, 2018, 70(3): 255  
Permissions
Breathers and Rogue Waves Derived from an Extended Multi-dimensional N-Coupled Higher-Order Nonlinear Schrödinger Equation in Optical Communication Systems
Bai Cheng-Lin1, 3, Cai Yue-Jin1, 2, †, Luo Qing-Long1, 3
School of Physics Science and Information Engineering, Liaocheng University, Liaocheng 252059, China
School of Computer and Information, Anqing Normal University, Anqing 246133, China
Shandong Province Key Laboratory of Optical Communications Science & Technology, Liaocheng 252059, China

 

† Corresponding author. E-mail: caiyj456@163.com

Supported by the National Natural Science Foundation of China under Grant No. 61671227 and the Natural Science Foundation of Shandong Province in China under Grant No. ZR2014AM018

Abstract
Abstract

In this paper, an extended multi-dimensional N-coupled higher-order nonlinear Schrödinger equation (N-CHNLSE), which can describe the propagation of the ultrashort pulses in wavelength division multiplexing (WDM) systems, is investigated. By the bilinear method, we construct the breather solutions for the extended (1+1), (2+1) and (3+1)-dimensional N-CHNLSE. The rogue waves are derived as a limiting form of breathers with the aid of symbolic computation. The effect of group velocity dispersion (GVD), third-order dispersion (TOD) and nonlinearity on breathers and rogue waves solutions are discussed in the optical communication systems.

Keyword:N-coupled higher-order nonlinear Schrödinger equations;multi-dimensional;rogue waves;breathers;bilinear transformation
1 Introduction

Rogue waves, first emerged in the studies of oceans, are also called freak waves or extreme waves. It has been ten years since rogue waves were introduced into optics by Solli et al.[1] Then in recent years, the study of breathers and rogue waves have attracted much interests in the researches of nonlinear optics,[2–4] fluid mechanics,[5] Bose-Einstein condensate (BEC)[6] and so on. Various mathematical methods such as the Darboux transformation[7–9] and the bilinear method[10–13] were applied to the construction of the explicit expressions of rogue wave solutions.

The propagation of nonlinear pulse can be well described by the nonlinear Schrödinger equation (NLSE):

where q represents the complex envelope amplitudes, ς and κ are the normalized distance along the direction of the propagation and retarded time. Breathers and rogue waves for NLSE were reported in many literatures.[14–17] However in modern optical systems, in order to increase bit rates, the width optical pulses have already come to be femtosecond.[18–19] When we consider the higher-order effects in the high speed communication systems, the ultrashort pulses can be governed by the higher-order nonlinear Schrödinger equation (HNLSE).[20–23] Also the technology of multiplexing has become a significant topic, which draws much attention of researchers.[24–25]

In wavelength division multiplexing (WDM) systems, there are more fields need to be considered simultaneously. Therefore we intend to investigate an extended multi-dimensional N-coupled higher-order nonlinear Schrödinger equations (N-CHNLSE), which can describe the propagation of optical breathers and rogues waves in optical fibers.

where qj (j = 1, 2, …, N) are the complex envelope amplitudes of the channel j in all N transmission channels, N is the number of all multiplexing channels in WDM system, k″ and k′″ describe the effect of group velocity dispersion (GVD) and third-order dispersion (TOD) respectively, β is the effective nonlinear coefficient, γ and γs describe self-steepening (SS) and the low frequency component of stimulated Raman scattering (SRS). The parameter c is the coefficient of the group velocity, which we have extended into the normal CHNLSE in order to make a balance between TOD and nonlinearity (see expression (9)) when we construct breathers solutions. When ∇ is defined as (∂/∂x), (∂/∂x, ∂/∂y), or (∂/x, ∂/∂y,∂/∂z), Eq. (2) respectively means the (1+1), (2+1) and (3+1)-dimensional extended N-CHNLSE. Here, x, y, z, t are normalized distances along the each propagation direction in the space and retarded time.

Recently, breathers and rogue waves have been discovered for coupled nonlinear Schrödinger equations (CNLSE),[26–28] even in multi-dimensions.[29] Nonautonomous multi-peak solitons, breathers and rogue waves for the inhomogeneous nonlinear Schrödinger equations have been studied in Refs. [30–31] including the dynamics of Peregrine combs and Peregrine walls.[32] Breathers and rogue waves have been obtained for CHNLSE.[33–34] To our knowledge, breathers and rogue waves for the extended multi-dimensional N-coupled higher-order nonlinear Schrödinger Eq. (2) have not been derived yet. The structure of this paper is organized as follows. In Sec. 2, comprehensive bilinear forms and breathers solutions for Eq. (2) are derived with the help of symbolic computation. Discussion on features of breathers and rogue waves solutions is arranged in Sec. 3. Section 4 is allotted for a conclusion.

2 Bilinear Forms and Solutions

The bilinear method used to obtain the soliton solutions can also be applied to get the breathers and rogue waves for nonlinear Schrödinger equations as shown in this paper. Here we set the bilinear transformation as:

where σj (j = 1,2, …,N) are real amplitude of background waves, ω is real angular frequency, g is complex, while f is real. We construct g and f as:

where ξ1 and ξ2 are complex parameters, M is a real parameter, and the star denotes complex conjugation. For the convenience of description, we choose m + 1 (m = 1, 2, 3) to represent (1+1), (2+1) and (3+1) dimensions of equation. Here, we define m as the number of space coordinate axes in this study. In expression (4), θ has the different definitions, in (m+1)-dimensional equations

where φ represents the phase factors. Without loss of generality while creating the comprehensive bilinear forms for the extended multi-dimensional N-CHNLSE, we intend to set ξ1 = ξ2 and p = p1 = p2 = p3. In addition, p and δ are real parameters. Combining expressions (3), (4), (5) with Eq. (2), the bilinear forms can then be written as

with γ + γs = 0, Pm = m. λ is a real parameter to be determined. Also the following provision is established in our study:

While Dz and Dt are Hirota bilinear operators in Eq. (6) defined

According to the conditions we employed above plus substituting expressions (4), (5), (7) into bilinear forms (6), we will obtain

Now the solutions of the extended multi-dimensional N-CHNLSE are obtained, when we combine all the conditions in this section. The solutions for (m+1)-dimensional Eq. (2) can be written as follows:

The breathers and rogue waves solutions for the (1+1), (2+1) and (3+1)-dimensional extended N-CHNLSE can be described as combinations of expressions (10a) with (10b), (10c), or (10d) respectively.

3 Discussions on Breathers and Rogue Waves for Multi-dimensional Eq. (2)

In this section, we will discuss the characters of the solutions (10) with the effects of GVD, TOD, and nonlinearity. Due to space limitations in this paper, we just take (1+1) and (2+1)-dimensional Eq. (2) as examples to discuss.

3.1 Breathers and Rogue Waves for (1 + 1)-Dimensional Eq. (2)

According to solutions we got in expression (10), the breathers solutions for (1+1)-dimensional 2-CHNLSE are shown in Fig. 1, where breathers of |q1| and |q2| are displayed. It is noticed that the breathers are localized in time and periodic in the transmission direction, which were termed as “Akhmediev” breather in Ref. [35]. Furthermore, the breathers we got in this paper are all considered to be Akhmediev breathers. Figure 1(c) is the profiles of breathers in two channels, which shows the differences in amplitudes with the same period between each transport channel. It shows that each signal propagates with same speed but different amplitudes. Besides, we found there will be rogue waves solutions, when the parameters come to limit in our study.

Fig. 1 (Color online) Breathers for (1+1)-dimensional 2-CHNLSE via solution (10) with σ1 = 0.6, σ2 = 0.8, k′″ = − 1, kk″ = 1, p = 1, γ = 1, φ = 0. (a) |q1|, (b) |q2|, (c) The profiles of breathers q1 and q2 with t = 0.

To generate rogue waves, we restrict the parameter p tending to zero with the following provision:

Then the rogue waves solution for (1+1)-dimensional 2-CHNLSE can be defined as:

In Figs. 2(a) and 2(b), we will see rogue waves of |q1| and |q2|, which shows that they are completely symmetric about x = 0 and t = 0. The rogue waves appear as localized large amplitude waves on a rough background, which have two obvious characteristics. They not only appear from nowhere and disappear without a trace, but also exhibit a dominant peak.[36] The amplitudes of rogue waves are obtained as follows on the axis of symmetry, with the condition of t = 0:

The effects of TOD and nonlinearity on the widths of rogue waves are explicated in Fig. 3(a). Here, we just take k′″γ < 0, to keep the parameters δ we have defined before real. Furthermore, the width of rogue waves for (1+1)-dimensional N-CHNLSE are generated as

where σ has been defined in Sec. 2. In Fig. 3(b), the width of rogue waves decrease rapidly with increasing number of channels when we set k′″ = − 1, γ = 1.

Fig. 3 (Color online) The evolution of width of rogue waves (a) via expression (14) with via σ1 = 0.6, σ2 = 0.8, (b) via expression (15) with k′″ = −1, γ = 1. σn = 0.6, n = 1,2, …,N.

The localized characters including length and width of the rogue wave have been studied for other equations.[37] Here we define the length of rogue waves at height 2σn. So the length of the rogue waves can be written as

By changing the value of GVD k″ in Fig. 4(a) we shall discover the changes in the length of rogue wave comparing with Fig. 2(a). When there is no group velocity dispersion in optical communication system, the rogue wave turns to be a W-shaped soliton wave in Fig. 4(b) with an infinite length. The amplitude of W-shaped soliton wave have been depicted in Fig. 4(c) at t = 0. Comparing Fig. 4(c) with Fig. 2(c), there is no difference in amplitude of q1. It proves that the value of GVD does not affect the amplitudes of signals in optical system.

Fig. 4 (Color online) Evolution of length via solution (12a) with σ1 = 0.6, σ2 = 0.8, k′″ = − 1, γ = 1. (a) Rogue wave with k″ = 0.2, (b) W-shaped soliton wave with k″ = 0, (c) The profiles of W-shaped soliton wave with k″ = 0, t = 0.

In Figs. 5(a), 5(b), and 5(c), the solution status of breathers for (1+1)-dimensional N-CHNLSE are changing with increasing number of channels. Figure 5(d) describes the profiles of a period of breather illustrates that the peaks, valleys and width of breathers are different depending on the number of channels in multiplexing systems. Exactly speaking, the number increased with the rising peaks, descending valleys and narrowing width.

Fig. 5 (Color online) Evolution of breathers for (1+1)-dimensional N-CHNLSE via expression (10) with k′″ = − 1, k″ = 1, p = 1, γ = 1. σn = 0.6, n = 1,2, …, N. (a) N = 4, (b) N = 8, (c) N = 16, (d) The profiles of a period of breather with t = 0.
3.2 Breathers and Rogue Waves for (2+1)-Dimensional Eq. (2)

The breather solutions for (2+1)-dimensional N-CHNLSE are described by Eqs. (10a) and (10c). Figures 6(a) and 6(b) display the breather solutions of q1 in the x−t and y−t planes. In both x, y directions, signal transmission form are basically the same. In Fig. 6(c), it is noticed that there are periodic linear wave solutions in x−y plane. In addition, the profiles of q1 will be shown in Fig. 7(a) along y = x. In Fig. 7(b), the amplitude of q1 is formed as breathers with the continuous-time. To generate the rogue waves, by the means mentioned above, the rogue wave solution for (2+1)-dimensional 2-CHNLSE can be written as follows.

Fig. 6 (Color online) The solutions of q1 for (2+1)-dimensional 2-CHNLSE via expression (10) with σ1 = 0.6, σ2 = 0.8, k′″ = − 1, k″ = 1, p = 0.8, γ = 1. φ = 0. (a) Breather in the x−t plane, (b) Breather in the y−t plane, (c) The periodic linear waves in the x−y plane.
Fig. 7 (Color online) The periodic linear waves solutions with σ1 = 0.6, σ2 = 0.8, k′″ = − 1, k″ = − 1, p = 1, γ = 1. φ = 0. y = x. (a) The profiles of periodic linear waves solutions in the discrete-time, (b) The breathers formed in the continuous-time.

In Figs. 8(a) and 8(b), the rogue wave solutions are shown in x−t and y−t planes. The signals travel as rogue waves in different derections. However, it is noticed that there is W-shaped soliton wave in x−y plane depicted in Fig. 8(a). The amplitudes of rogue waves in x−t and y−t planes are obtained as follows with the condition t = 0:

Fig. 8 (Color online) The solutions of q1 for (2+1)-dimensional 2-CHNLSE via solutions (16) with σ1 = 0.6, σ2 = 0.8, k′″ = − 1, k″ = − 1, γ = 1. φ = 0. (a) Rogue wave in the x−t plane, (b) Rogue wave in the y−t plane, (c) The W-shapedsoliton waves in the x−y plane.

According to Eq. (18), we get different amplitudes in Fig. 9(a), when we change the values of y. In x−y plane, when we choose y = x, the profiles of the W-shaped soliton wave in different time are displayed in Fig. 9(b). Similarly, the amplitude is formed as rogue wave with respect to t and x in Fig. 9(c) with the continuous-time. The widths of rogue waves are expressed as

which are all the same in x−t, y−t, and x−y planes. In the same way, the width of rogue waves for (2+1)-dimensional N-CHNLSE is described as

From expressions (14), (15), (19), and (20), it is indicated the additional dimension does not affect widths of rogue wave solutions.

Fig. 9 (Color online) (a) Amplitudes ofrogue waves in the x−t plane, (b) The profiles of W-shaped soliton waves in the x−y plane, (c) The rogue wave formed in the continuous-time.
4 Conclusions

In this study, our main objective has focused on an extended multi-dimensional N-CHNLSE, which governs the propagation properties of ultrashort pulses in N optical fields. Based on a series of settings, we have constructed the comprehensive bilinear forms (6) and the solutions (10) for (1+1), (2+1), and (3+1)-dimensional N-CHNLSE. According to the solutions, this paper puts emphasis on discussing the features of breather and rogue waves for (1+1) and (2+1)-dimensional N-CHNLSE. The features of the solutions are principally amplitude, peak, valley, width and length which are revealed the relations to GVD, TOD, nonlinearity and the number of space dimensions and transmission channels. In addition, as for improving the comprehensive degree of bilinear forms and solutions, there are more interesting special cases we have not shown in this paper and will be researched later.

Reference
[1] Solli D. R. Ropers C. Koonath P. Jalali B. Nature (London) 450 2007 1054
[2] Akhmediev N. Dudley J. M. Solli D. R. Turitsyn S. K. J. Opt. 15 2013 0201
[3] Kibler B. Fatome J. Finot C. et al. Nat. Phys. 6 2010 790
[4] Flach S. Gorbach A. V. Phys. Rep. 467 2008 1
[5] Chai J. Tian B. Sun W. R. Xie X. Y. Comput. Math. Appl. 71 2016 2060
[6] Wu X. Y. Tian B. Zhen H. L. Mod. Opt. 63 2016 1
[7] Wang L. H. Porsezian K. He J. S. Phys. Rev. E 87 2013 053202
[8] Guo B. Ling L. Liu Q. P. Phys. Rev. E 85 2012 026607
[9] Guo R. Zhao H. H. Wang Y. Nonlinear Dynam. 83 2016 2475
[10] Hirota R. J. Math. Phys. 14 1973 805
[11] Cai Y. J. Bai C. L. Luo Q. L. Commun. Theor. Phys. 67 2017 273
[12] Xu Z. X. Chow K. W. Appl. Math. Lett. 56 2016 72
[13] Wazwaz A. M. Appl. Math. Lett. 58 2016 1
[14] Mirzazadeh M. Eslami M. Zerrad E. et al. Nonlinear Dynam. 81 2015 1
[15] Yan Z. Phys. Lett. A 374 2010 672
[16] He J. S. Zhang H. R. Wang L. H. et al. Phys. Rev. E 87 2013 052914
[17] Kedziora D. J. Ankiewicz A. Akhmediev N. Phys. Rev. E 85 2012 066601
[18] Agrawal G. P. Nonlinear Fiber Optics 5 Academic Press New York 2013
[19] Liu C. Yang Z. Y. Zhao L. C. et al. Phys. Rev. E 94 2016 042221
[20] Ankiewicz A. Soto-Crespo J. M. Akhmediev N. Phys. Rev. E 81 2010 046602
[21] Choudhuri A. Porsezian K. Phys. Rev. A 85 2012 1431
[22] Triki H. Biswas A. Milović D. Belić M. Opt. Commun. 366 2016 362
[23] Soto-Crespo J. M. Devine N. Hoffmann N. P. Akhmediev N. Phys. Rev. E 90 2014 032902
[24] Adhikary A. Nam J. Ahn J. Y. Caire G. IEEE T Inform. Theory 59 2012 6441
[25] Richardson D. J. Fini J. M. Nelson L. E. Nat. Photon 7 2013 354
[26] Liu C. Yang Z. Y. Zhao L. C. Yang W. L. Phys. Rev. A 89 2014 055803
[27] Cheng X. Wang J. Li J. Nonlinear Dynam. 77 2014 545
[28] Guo B. Ling L. Chin. Phys. Lett. 28 2011 110202
[29] Manikandan K. Senthilvelan M. Kraenkel R. A. Eur. Phys. J. B 89 2016 218
[30] Wang L. Li X. Qi F. H. Zhang L. L. Ann. Phys. 359 2015 97
[31] Cai L. Y. Wang X. Wang L. Nonlinear Dynam. 90 2017 2221
[32] Wang L. Wang Z. Q. Sun W. R. Commun. Nonlinear Sci. 47 2017 190
[33] Qi F. H. Xu X. G. Wang P. Appl. Math. Lett. 54 2016 60
[34] Zhang Y. Nie X. J. Zha Q. Phys. Lett. A 378 2014 191
[35] Akhmediev N. Korneev V. I. Theor. Math. Phys. 69 1986 1089
[36] Akhmediev N. Ankiewicz A. Taki M. Phys. Lett. A 373 2009 675
[37] He J. S. Wang L. H. Li L. J. et al. Phys. Rev. E 89 2014 062917