Modulation instability,rogue waves and conservation laws in higher-order nonlinear Schr?dinger equation

Min-Jie Dong and Li-Xin Tian,2,*

1 School of Mathematical Sciences,Nanjing Normal University,Nanjing,Jiangsu 210023,China

2Nonlinear Scientific Research Center,Jiangsu University,Zhenjiang,Jiangsu 212013,China

Abstract In this paper,the modulation instability (MI),rogue waves(RWs) and conservation laws of the coupled higher-order nonlinear Schr?dinger equation are investigated.According to MI and the 2 × 2 Lax pair,Darboux-dressing transformation with an asymptotic expansion method,the existence and properties of the one-,second-,and third-order RWs for the higher-order nonlinear Schr?dinger equation are constructed.In addition,the main characteristics of these solutions are discussed through some graphics,which are draw widespread attention in a variety of complex systems such as optics,Bose-Einstein condensates,capillary flow,superfluidity,fluid dynamics,and finance.In addition,infinitely-many conservation laws are established.

Keywords: higher-order nonlinear Schr?dinger equation,modulation instability,rogue waves,conservation laws

1.Introduction

Rogue waves(RWs)are regarded as one of the highly unsafe phenomena in the ocean [1],and are usually called monster waves or extreme waves in the ocean.In addition to the ocean,RWs can exist in optical fibers,Bose-Einstein condensates [2],and finance [3,4].The peregrine soliton as the classic NLS equation is seen as a possible mathematical explanation in 2010 [2].

Modulation instability(MI)is one of the main reasons for the formation of RWs in nonlinear dispersion systems,and its formation is intimately related to baseband MI.MI is found in the environment of water waves and is also known as the instability of Benjamin Fair [5].MI can be interpreted as the exponential growth of the initial sine wave perturbation of the plane wave solution.In addition to water waves,MI also exists in plasma,nonlinear optics [6] and Bose-Einstein condensate.However,according to [7-9],we know that not every type of MI leads to the formation of the RW.The baseband MI is defined as a condition where the cw background is unstable relative to a disturbance with an infinitesimal frequency by Baronio et al in 2014 [10].The passband MI is defined as the interference in which disturbances are experienced in a spectral region where the limit frequency does not include the zero frequency.This shows that in different nonlinear wave models,the existence conditions of RW solutions are consistent with baseband MI conditions.Experiments in optical fibers are designed to provide evidence of the passband and baseband polarization MI in the defocused Manakov system by Frisquet et al in 2015 [11].Recently,Yang research group studied MI and related localized wave excitations in general high-order NLSE,see [12,13].On the other hand,the MI developed from localized perturbations can also studied by the breathers exactly,namely the super-regular breathers,see [14-16].

Conservation laws play an important role in discussing the integrability of soliton equations.Since Dogan discovered the infinite number of conservation law of the KdV equation[17],many methods have been developed.For example,the infinite number of conservation laws or conserved quantities of a continuous system can be obtained through the B?klund transformation,the couple of Ricatti equations [18],and the scattering problem [19].Recently,Tsuchida and Wadati proposed a graceful trace identity as a component extension to describe the protection of multi-component cases [20].In addition,all the above methods have been extended to discrete soliton lattices.

In this paper,we shall use the Darboux-dressing transformation (DDT) to concentrate the rational solutions of the higher-order nonlinear Schr?dinger (HONLS) equation [21].

where ψ is the wave function in optics,and t is the time variable,x is the propagation variable.The superscripts denote the complex conjugate and the subscripts represent the partial derivatives.Within our knowledge,the HONLS equation contains three completely integrable equations,which are α=0,δ=0;α ≠ 0,δ=0;α=0,δ ≠ 0.These equations have been studied by many scholars from different perspectives [22-25].

To facilitate calculation,Liu and his colleagues assumed r=iψ andq=,the equation (1) can be transformed into following coupled higher-order nonlinear Schr?dinger(CHONLS) equation [26]

To our knowledge,the MI,RWs and conservation laws of CHONLS equations have not been reported in the existing literatures.In this paper,the breather solution and the RWs of this equation through DDT are mainly studied.The change of the parameters α and δ will affect the propagation direction of the wave,which will shown in detail with the figure.We also study the relationship between the RWs and the MI of the equation,and the conservation law.

The remainder of our article is constructed as follows.In section 2,according to the modulating instability,the linear stability of nonlinear plane waves with variable coefficients in the presence of small perturbations will be analyzed.In section 3,the DDT of the Lax pair system will be presented.In section 4,the new breather wave and RWs of equation(2)will be systematically derived.The first-order,second-order,and third-order accurate RWs are given and their dynamic characteristics will be analyzed.In section 5,the conservation laws will be constructed.Our conclusions will be drawn in section 6.

2.Modulation instability

In this section,we pay attention to MI on the plane wave state of the equation (2),which is believed to be the cause of the formation of RWs.According to [12-16],the plane wave solutions of equation (2) have the following form

where a denotes the amplitude,b denotes the frequency,both of them are real constants,and c is the plane wave number.Substituting equation (3) into (2),it can be obtained that

In order to perform the linear stability analysis,we add a small perturbation term in the the plane wave solution

where ε is a small parameter,Q and R are the function of x and t.Substituting the above solution (5) into the (2),then collecting the terms in the first order of ε,we can get the following perturbation equation for Q and R

Noting the linearity of the above equation (6) with respect to R and Q,we assume R and Q as

where Λ is the wave number,Ω is the modulation frequency,while U and V are small parameters.Putting equation(8)into(6),the following dispersion relation for the perturbations are obtained

Figure 1.Gain spectra of modulation instability for the parameter values:(a)α=1,δ=1; (b)α=0,δ=0;(c)α=1,δ=0;(d)α=0,δ=1 ,and for the different power levels as legend.

where

the existence conditions of the solutions for U and V (i.e.UV ≠ 0) can be obtained

Solving equation (11),the following dispersion relation are obtained

The power gain at any frequency Ω is obtained from equation (12) and is given by

in which G stands for the gain with (Λ2-4a2)Figure 1 shows the gain spectra at three power levels.

3.Lax pair and the DDT

Due to complete integrability [26],equation (2) can be cast into the following 2 × 2 linear eigenvalue problem

where Φ= (φ1,φ2)′(′mean a matrix transpose)is the vector eigenfunction,φ1and φ2are the complex functions of(x,t),M and N are the 2 × 2 square matrices,and

with i2=-1,J=diag(i,-i) and

the compatibility condition Mt- Nx+ MN - NM=0.

Theorem 1.The following unified DDT yields

where

and I2=diag(1,1),Φ is a special vector for the lax pair with λ=λ1.

It is necessary to show that the linear system can be rewritten in next form

and transformation between potential functions reads as

where commutator [A,B]=AB - BA.

4.Breather wave solutions and RWs

In what follows,we will construct breathing waves and RWs of the equation (2) based on the DDT.

4.1.Breather wave solutions

Based on the [27-34],the corresponding solution of the Lax pair can be sought in a new form

where Z is an arbitrary complex vector.It is easy to obtain that R and S meet

where

Through complex calculations,the exponential matrices F can be written as

where

Similarly,the G is obtained

Figure 2.Breather solutions in equation (2) for parameters:α =,d = (a),(b),(c)λ = i; (d),(e),(f)λ =

where

Whenλ=-±ia,the exponential exp (iRx+ iSt)can be transformed into a combination of polynomial and exponential functions.

where

4.2.Rogue waves

According the Theorem,the first-order RWs arrive at

Figure 3.One-order RWs in equation (2) for parameters (a),(d):α =,d =; (b),(e):α =,d =; (c),(f):α =,d=

where

with

It is easy to observe that the expression of R1contains two arbitrary parameters α and δ,which are the third-order dispersion coefficient and the four-order dispersion coefficient respectively.With the increase of α and δ,the crest of the RWs deflect clockwise and its width decreases.Figure 3 illustrates the above dynamic characteristics.

Similarly,the second-order RWs arrive at

where

The F1and G1expressed as

where

Figure 4.Second-order RWs in equation (2) for parameters (a),(d):α =; (b),(e)α =; (c),(f):α=

Figure 5.Second-order RWs in equation(2)for parameters α = (a):s0=0 s1=1,(b):s0=0 s1=10,(c):s0=0 s1=100.

Figure 6.The third-order RWs in equation (2) for parameters (a),(d):α =; (b),(e):α=; (c),(f):α=

Figure 7.The third-order RWs in equation(2)for parameters α= .(a)s0=0 s1=1 s2=1,(b)s0=0 s1=10 s2=10,(c)s0=0 s1=100 s2=100.

As shown in figure 4,a three-dimensional diagram of the second-order RWs are plotted,according to which we observe that the main peak is surrounded by four lower peaks when s1=0.However,when amplify s1,the single peak can split into three peaks,which are symmetric about the straight line t=0 in the(x,t)-plane of figure 5.The solution is called the‘three sisters’ or a ‘RW triplet’.

Due to the complex expressions showing of the higher-order solutions,the third-order RWs are plotted on here.Figures 7 and 8 illustrate two different spatial and temporal distribution patterns of the third-order solution under the same parameter values and different sj(j=1,2,3),from which,we can see that the wave presents a ring pattern with the internal second-order fundamental mode.The distance between a ring pattern and inner secondorder fundamental pattern increases with the increase of s2.Figure 7 is the ring-triangle distribution,figure 8 shows a fivepass distribution including six first-order basic patterns.

5.Conservation laws

In this section,according to [35-37] we present infinitely many independent conservation laws as a further support of the integrability for equation (2),from its Lax pairs.By introducing the functionthe Ricatti equation obtained

Supposing Γ=Γ-1λ + Σn=0Γnλ-n,where Γnare the functions of x and t to be determined,substituting it into expression (39) and equating the coefficients of the same power of Γ to zero,the recurrence relations are obtained,

Figure 8.The third-order RWs in equation (2) for parameters α = s0 =0,s1=0,(a) s2=10,(b) s2=100,(c) s2=1000.

From the compatibility conditionone can get

Substituting expressions(5)into(40),the infinitely-many conservation laws for System (2) derived as

with

where Rnand Snrepresent the conserved fluxes and conserved densities,respectively.

6.Conclusions

In this paper,we studied the MI of equation (2),which described MI for the possible generation mechanism of RWs.Based on the DDT method,we studied the first-,second-and third-order RW solutions under different forms of equation(2).In order to help the readers better understand the solutions,figures 2-7 give the breather wave and RWs by looking for the appropriate parameters,respectively.Based on Lax pair (15),an infinite number of conservation laws (41)was constructed to prove the integrability of equation (2).

Acknowledgments

We express our sincere thanks to the editor and the referees for their valuable comments.This work is supported by the National Natural Science Foundation of China (Grant No.71690242,No.11731014,No.12001241) and the Basic Research Program of Jiangsu Province (Grant No.BK20200885).