# Effect of Third-Order Dispersion on the Solitonic Solutions of the Schrödinger Equations with Cubic Nonlinearity

Authors: H. Chachou Samet; M. Benarous; M. Asad-uz-zaman; U. Al Khawaja
Journal: Advances in Mathematical Physics
 (2014)

---

## Abstract

We derive the solitonic solution of the nonlinear Schrödinger equation with cubic nonlinearity, complex potentials, and time-dependent coefficients using the Darboux transformation. We establish the integrability condition for the most general nonlinear Schrödinger equation with cubic nonlinearity and discuss the effect of the coefficients of the higher-order terms in the solitonic solution. We find that the third-order dispersion term can be used to control the soliton motion without the need for an external potential. We discuss the integrability conditions and find the solitonic solution of some of the well-known nonlinear Schrödinger equations with cubic nonlinearity and time-dependent coefficients. We also investigate the higher-order nonlinear Schrödinger equation with cubic and quintic nonlinearities.


---

## Body

## 1. Introduction

The nonlinear Schrödinger equation (NLSE) appears in different branches of physics [1–3]. In optical fiber communication systems, this type of equation describes the dynamics of an optical pulse through the fiber. Due to the advent of high power optical soliton, higher-order terms become very important in describing the soliton propagation accurately through the fiber. The role of the higher-order terms in the stability and dynamics of the solitons thus needs to be understood. To that end, we consider here the most general form of the NLSE with cubic nonlinearity, complex potential, and derivative terms. There are many different approaches to solve the nonlinear partial differential equations, some of which are developed very recently. A few of these approaches are similarity transformation method [4–6], Hirota bilinear method [7], the generalized and improved (G′/G)-expansion method [8, 9], the mapping method [10, 11], and Darboux Transformation (DT) method. We solve this general equation with the DT method. It has the beauty of generating nontrivial solutions from a trivial seed solution of a given integrable NLSE. We discuss the integrability conditions and solve some of the well-known NLSEs with cubic nonlinearity and time-dependent coefficients. The effect of the higher-order terms will then be evident. In some cases, we include a linear potential and a term describing the damping or gain in the system.The key ingredient to solve a NLSE using the DT is a set of matrices, commonly known as the Lax pair. We derive the Lax pair for a given NLSE using a systematic search approach [12]. The aim of this paper is to use the DT to get an exact solitonic solution for the NLSE and study the effect of the higher-order terms in the solution. Specifically, we consider NLSE up to cubic nonlinearity only; some of the well-known NLSEs with cubic nonlinearity are NLSE or Gross-Pitaevskii equation (GPE), Sasa-Satsuma equation (SSE), and Hirota equation (HE). One of the key factors here is that the coefficients are time dependent. The coefficients of different terms in the NLSE can be made time-dependent simply by varying the radius or the nonlinearity along the fiber. Alternatively, the coupling in wave guide arrays may be suitably changed [13]. The time-dependent coefficients introduce a time-dependent amplitude and a position dependent phase, respectively.We investigate the integrability conditions of the NLSE with cubic and quintic nonlinearity and find that it is not integrable in the Lax pair sense. Although we manage to find two matricesU0 and V0 that generate the given NLSE as their compatibility condition, we were unable to find the full set of matrices for the spectral parameters required to solve the given NLSE using the DT.The rest of the paper is organized as follows. In Section2 basic equations for deriving the Lax pair and DT are presented. In Section 3 solutions are found for some well-known NLSEs with cubic nonlinearity. Discussion and conclusions are in Section 4.


## 2. Lax Pair and Darboux Transformation

The Lax pair is the most important ingredient in solving a NLSE using the DT. Here we summarize the basic equations only. In order to derive the Lax pair we start by writing the following two equations:(1)Ψx=U0·Ψ+U1·Ψ·Λ,Ψt=V0·Ψ+V1·Ψ·Λ+V2·Ψ·Λ2+V3·Ψ·Λ3
for an auxiliary field Ψ. Here subscripts t and x represent the derivative with respect to time and position, respectively. The matrices U0, U1, V0, V1, V2, and V3 are functions of the solution of the given NLSE and its space derivatives and Λ denotes the spectral parameter matrix. The order of the matrices depends on the NLSE to be solved as will be seen next. The compatibility condition Ψxt=Ψtx leads to the following set of equations:
(2)U0t-V0x+[U0,V0]=0,(3)U1t-V1x+[U0,V1]+[U1,V0]=0,(4)V2x+[V2,U0]+[V1,U1]=0,(5)V3x+[V3,U0]+[V2,U1]=0,(6)[V3,U1]=0,
where [X,Y] denotes the commutator of X and Y. Equation (2) gives back the original NLSE, which is the condition that the Lax pair must satisfy. The most important and difficult part is to find U0 and V0 matrices for a given NLSE. The remaining matrices are determined from U0 and V0 matrices. We use a systematic search approach to find U0 and V0 matrices [12]. It is worth mentioning that for some higher-order NLSE we find U0 and V0 matrices but do not find the remaining matrices. Therefore, analytical solution of the given higher-order NLSE is not possible using the DT.After deriving the set of matrices for the given NLSE we use the DT to find the analytic solution. We solve (1) for the auxiliary field Ψ=Ψ0. Using this auxiliary field we find the solution of the original NLSE, which is given by
(7)U0[1]=U0+[U1,σ],
where U0[1] is the new U0 matrix and σ=Ψ0·Λ·Ψ0-1. This U0[1] contains the solution of the given NLSE. The beauty of the DT lies in the fact that a trivial seed solution, for example, zero or continuous wave solution of the given NLSE, leads to a nontrivial solution.


## 3. NLSE with Cubic Nonlinearity

In this section we consider the most general NLSE with qubic nonlinearity [1, 14–17]. In the presence of damping or gain and external potential, the general NLSE can be written as
(8)iψt+a1ψxx+a2|ψ|2ψ+ia3ψxxx+ia4|ψ|2ψx+ia5ψ2ψx*+iΓψ+Vextψ=0.
In [18, 19], the Lax pair of a special case of this equation was found. Here all the coefficients ai≡ai(t) are initially assumed to be complex functions of time, but it turns out that integrability restricts them to be real. The functions Vext≡Vext(x,t) and Γ≡Γ(t) are real and account for the external potential and the damping or gain in the system, respectively. Equation (8) describes the evolution of the electric field envelope of an optical pulse through an optical fiber. The second term denotes the chromatic or group velocity dispersion, the third term denotes the Kerr nonlinearity or self-phase modulation (arises due to intensity dependent refractive index), and the fourth term denotes the third-order dispersion. The fifth and the sixth terms are related to the self-steepening and self-frequency shift due to the Raman scattering.Although we have found two2×2 matrices for U0 and V0, other integrability conditions do come in while deriving the remaining matrices; in addition to that we need to use 3×3 matrices. Employing the systematic search approach we find the complete set of matrices for (8), which is rendered to the Appendix for convenience.When we useU0 and V0 in (2) we get back (8) with
(9)Γ=a˙22a2-a˙12a1,Vext=a3a˙1-a1a˙33a32x+2a1327a32.
The remaining two integrability conditions are
(10)a4=9a2a32a1,a5=3a2a32a1.
Therefore (8) is integrable only with three independent parameters a1, a2, and a3. One point worth mentioning here is that the damping or gain, Γ, and the linear part of the potential, Vext, arise due to the time variation of the coefficients a1, a2, and a3.Applying the DT on the Lax pair of (8) the following solitonic solution can be derived:
(11)ψ(x,t)=Acosh⁡(2x-b)ei(a1/3a3)(x+π/2),
where A=2a1/a2 and
(12)b=∫‍(2a123a3+8a3)dt-12log⁡2.
The amplitude A depends on the ratio of the coefficients of the chromatic dispersion and the nonlinearity. The phase is position dependent and depends on the ratio of the coefficients of the chromatic and third-order dispersions. For repulsive interaction, that is, a2<0, the equation supports diverging solutions, one of which is given by
(13)ψ(x,t)=Asinh⁡(2x-b)ei(a1/3a3)(x+π/2),
where A=2-a1/a2 and b is given by (12).In both the attractive and repulsive cases the third-order dispersion coefficient controls the phase of the envelope. One important factor to note that the self-steepening and the Raman scattering terms arise due to the third-order dispersion through the relations given by (10). The group velocity of the envelope depends on the coefficients a1 and a3. Therefore, we can control the group velocity, that is, the center of mass motion of the soliton, by properly choosing the chromatic dispersion and the third-order dispersion. This interesting fact can be used to manipulate soliton motion without the need of an external potential.Figure1 shows the solitonic solutions given by (11). For comparison, Figure 1(a) shows the solitonic solution of NLSE with constant coefficients, that is, without any damping and external potential. In Figure 1(b), we choose all the three coefficients a1–a3 to be time independent and take a3 to be very small compared to a1 and a2. In this particular case Γ=0 and Vext=2a12/27a32 and the center of mass velocity of the soliton is a constant. We assume trigonometric time dependence for a1 and a2 [20] and constant a3 in Figure 1(c). The soliton moves with a constant velocity with oscillation around it. In Figure 1(d), we assume all the coefficients be trigonometrically time dependent and the center of mass velocity oscillates. From Figures 1(b)-1(c), it is obvious that the center of mass velocity can be controlled by appropriate choice of the third-order dispersion coefficient a3.(a) Solitonic solution of NSLE with constant coefficients,a1=1 and a2=1 (color online). (b) The solution of (8) with constant coefficients, namely, a1=1, a2=1, and a3=0.05. (c) Trigonometric time dependence for a1=cos⁡(t), a2=cos⁡(t), and constant third-order dispersion a3=0.05. (d) All the three coefficients are trigonometrically time dependent a1=cos⁡(t), a2=cos⁡(t), and a3=0.05cos⁡(t).
(a)
(b)
(c)
(d)


### 3.1. Special Cases

#### 3.1.1. Nonlinear Schrödinger Equation, NLSE

If we seta3=a4=a5=0 in (8) we get the standard NLSE or GPE with damping or gaining and potential, which can be written as
(14)iψt+a1ψxx+a2|ψ|2ψ+iΓψ+Vextψ=0.
This is one of the most well-known NLSEs in the optical [1], Bose-Einstein [2], and ocean wave communities [3]. The Lax pair for (14) with time-independent coefficient and without damping and potential was first found by Zakharov and Shabat [21] and was then widely used in the Darboux transformation to find exact solutions. The damping is given by Γ=a˙2/(2a2)-a˙1/(2a1) and the potential is given by Vext=c12a1-c˙1x. Equation (14) supports solitonic solutions and one such solution is given by
(15)ψ(x,t)=Acosh⁡[c3x+b]ei[c32∫a1dt-c1x],
where A=-4c32a1/a2 and b=2c3∫c1a1dt. Similar results are also found by Kumar et al. [22]. Al Khawaja found solution of (14) with quadratic potential, periodic potential, and different types of solution with linear potential.It is well-known that the NLSE or GPE with time-dependent coefficients is integrable with potential up to quadratic power inx [23, 24]. Here, we find the Lax pair for NLSE with linear potential Vext=c12a1-c˙1x; this is due to the fact that we derive the Lax pair for the NLSE from that of the higher-order NLSE with cubic-quintic nonlinearity by taking appropriate limit of a3=a4=a5=a6=0 simultaneously [see (21)]. The higher-order terms in higher-order NLSE restrict integrability to linear potentials. One of the original motivations of our work is to solve the most general NLSE with cubic-quintic nonlinearity and find the solutions for other well-known NLSEs as special cases. However, we find that general NLSE with cubic-quintic nonlinearity is not solvable in the Lax pair sense.


#### 3.1.2. Hirota Equation, HE

If we seta5=0 in (8) we get the HE [25]; namely,
(16)iψt+a1ψxx+a2|ψ|2ψ+ia3ψxxx+ia4|ψ|2ψx+iΓψ+Vextψ=0.
In this case we find that the equation is integrable only if a4=3a2a3/a1, Γ=a˙4/(2a4)-a˙3/(2a3), and Vext=0. The one-soliton solution of the HE is given by
(17)ψ(x,t)=Acosh⁡[6x-b]e6i∫a1dt,
where A=-12a3/a4 and b=66∫a3dt. Similar result was also found by Dai and Zhang [26].


#### 3.1.3. Sasa-Satsuma Equation, SSE

If we seta1=a2=0 in (8) we get the SSE [27]; namely,
(18)iψt+ia3ψxxx+ia4|ψ|2ψx+ia5ψ2ψx*+iΓψ+Vextψ=0.
We find that the SSE is integrable for arbitrary values of a3, a4, and a5 with a4≠2a5 and Γ=Vext=0. The one-soliton solution of the SSE is given by
(19)ψ(x,t)=iAcosh⁡(2c2λx-b),
where A=2c1c2λ3/c3 and b=8c23(t-1)λ3a3. Similar result was found by Gilson et al. [28].We further investigate the integrability of the general higher-order NLSE equation with cubic and quintic nonlinearities [29–31]. In nonlinear optics, the nonlinear Schrödinger equation with quintic nonlinearity becomes important as the intensity and the inverse width of the optical pulse increase. In the presence of damping or gaining and an external potential, in nonlinear optics the general NLSE with cubic and quintic nonlinearities can be written as
(20)iψt+a1ψxx+a2|ψ|2ψ+ia3ψxxx+ia4|ψ|2ψx+ia5ψ2ψx*+a6|ψ|4ψ+ia7|ψ|4ψx+ia8|ψ|2ψ2ψx*+iΓψ+Vextψ=0.
Here, as for (8), all the coefficients ai(t) turn out to be real. The terms with coefficients a6, a7, and a8 account for the quintic nonlinearity and are generally known as non-Kerr nonlinearity. All the other terms are similar to those for (8). Using the systematic search approach we find the Lax pair for (20), which is given by
(21)U0=(c1ic2eGψ-ic2e-Gψ*c1),V0=(v11v12v21v22),
where the matrix elements vij are given by
(22)v11=-ic5|ψ|2+c22a3(ψ*ψx-ψψx*),v12=eG[ic4c5c2ψ-c5c2ψx-ic2a3ψxxlllllll-i3c2(a4+a5)|ψ|2ψ-i3c2a7|ψ|4ψc4c5c2],v21=e-G[-ic4c5c2ψ*-c5c2ψx*+ic2a3ψxx*lllllll+i3c2(a4+a5)|ψ|2ψ*+i3c2a7|ψ|4ψ*c4c5c2],v22=ic5|ψ|2-c22a3(ψ*ψx-ψψx*).
Here c2=(2a5-a4)/6a3, c4=-3a6/a7, c5=c22[a1(a4+a5)-3a2a3]/(3a5), and G=c3+ic4x, and c1 and c3 are two arbitrary time-independent constants. Using (21) in (2) we get (20) with the following integrability conditions: a7=9a3a5a6/(3a2a3-a1a4+2a1a5), a8=2a7/3, Γ=c˙2/c2, and Vext=c42c5/c22-c˙4x. This Lax pair is a fake one as we could not find the other matrices to solve the spectral problem given by (20). Therefore (20) is not integrable in the Lax pair sense. We also tried 3×3 matrices for U0 and V0, as we did for the NLSE with cubic nonlinearity, but we were unable to find the other matrices necessary to solve the spectral problem.


## 3.1. Special Cases

### 3.1.1. Nonlinear Schrödinger Equation, NLSE

If we seta3=a4=a5=0 in (8) we get the standard NLSE or GPE with damping or gaining and potential, which can be written as
(14)iψt+a1ψxx+a2|ψ|2ψ+iΓψ+Vextψ=0.
This is one of the most well-known NLSEs in the optical [1], Bose-Einstein [2], and ocean wave communities [3]. The Lax pair for (14) with time-independent coefficient and without damping and potential was first found by Zakharov and Shabat [21] and was then widely used in the Darboux transformation to find exact solutions. The damping is given by Γ=a˙2/(2a2)-a˙1/(2a1) and the potential is given by Vext=c12a1-c˙1x. Equation (14) supports solitonic solutions and one such solution is given by
(15)ψ(x,t)=Acosh⁡[c3x+b]ei[c32∫a1dt-c1x],
where A=-4c32a1/a2 and b=2c3∫c1a1dt. Similar results are also found by Kumar et al. [22]. Al Khawaja found solution of (14) with quadratic potential, periodic potential, and different types of solution with linear potential.It is well-known that the NLSE or GPE with time-dependent coefficients is integrable with potential up to quadratic power inx [23, 24]. Here, we find the Lax pair for NLSE with linear potential Vext=c12a1-c˙1x; this is due to the fact that we derive the Lax pair for the NLSE from that of the higher-order NLSE with cubic-quintic nonlinearity by taking appropriate limit of a3=a4=a5=a6=0 simultaneously [see (21)]. The higher-order terms in higher-order NLSE restrict integrability to linear potentials. One of the original motivations of our work is to solve the most general NLSE with cubic-quintic nonlinearity and find the solutions for other well-known NLSEs as special cases. However, we find that general NLSE with cubic-quintic nonlinearity is not solvable in the Lax pair sense.


### 3.1.2. Hirota Equation, HE

If we seta5=0 in (8) we get the HE [25]; namely,
(16)iψt+a1ψxx+a2|ψ|2ψ+ia3ψxxx+ia4|ψ|2ψx+iΓψ+Vextψ=0.
In this case we find that the equation is integrable only if a4=3a2a3/a1, Γ=a˙4/(2a4)-a˙3/(2a3), and Vext=0. The one-soliton solution of the HE is given by
(17)ψ(x,t)=Acosh⁡[6x-b]e6i∫a1dt,
where A=-12a3/a4 and b=66∫a3dt. Similar result was also found by Dai and Zhang [26].


### 3.1.3. Sasa-Satsuma Equation, SSE

If we seta1=a2=0 in (8) we get the SSE [27]; namely,
(18)iψt+ia3ψxxx+ia4|ψ|2ψx+ia5ψ2ψx*+iΓψ+Vextψ=0.
We find that the SSE is integrable for arbitrary values of a3, a4, and a5 with a4≠2a5 and Γ=Vext=0. The one-soliton solution of the SSE is given by
(19)ψ(x,t)=iAcosh⁡(2c2λx-b),
where A=2c1c2λ3/c3 and b=8c23(t-1)λ3a3. Similar result was found by Gilson et al. [28].We further investigate the integrability of the general higher-order NLSE equation with cubic and quintic nonlinearities [29–31]. In nonlinear optics, the nonlinear Schrödinger equation with quintic nonlinearity becomes important as the intensity and the inverse width of the optical pulse increase. In the presence of damping or gaining and an external potential, in nonlinear optics the general NLSE with cubic and quintic nonlinearities can be written as
(20)iψt+a1ψxx+a2|ψ|2ψ+ia3ψxxx+ia4|ψ|2ψx+ia5ψ2ψx*+a6|ψ|4ψ+ia7|ψ|4ψx+ia8|ψ|2ψ2ψx*+iΓψ+Vextψ=0.
Here, as for (8), all the coefficients ai(t) turn out to be real. The terms with coefficients a6, a7, and a8 account for the quintic nonlinearity and are generally known as non-Kerr nonlinearity. All the other terms are similar to those for (8). Using the systematic search approach we find the Lax pair for (20), which is given by
(21)U0=(c1ic2eGψ-ic2e-Gψ*c1),V0=(v11v12v21v22),
where the matrix elements vij are given by
(22)v11=-ic5|ψ|2+c22a3(ψ*ψx-ψψx*),v12=eG[ic4c5c2ψ-c5c2ψx-ic2a3ψxxlllllll-i3c2(a4+a5)|ψ|2ψ-i3c2a7|ψ|4ψc4c5c2],v21=e-G[-ic4c5c2ψ*-c5c2ψx*+ic2a3ψxx*lllllll+i3c2(a4+a5)|ψ|2ψ*+i3c2a7|ψ|4ψ*c4c5c2],v22=ic5|ψ|2-c22a3(ψ*ψx-ψψx*).
Here c2=(2a5-a4)/6a3, c4=-3a6/a7, c5=c22[a1(a4+a5)-3a2a3]/(3a5), and G=c3+ic4x, and c1 and c3 are two arbitrary time-independent constants. Using (21) in (2) we get (20) with the following integrability conditions: a7=9a3a5a6/(3a2a3-a1a4+2a1a5), a8=2a7/3, Γ=c˙2/c2, and Vext=c42c5/c22-c˙4x. This Lax pair is a fake one as we could not find the other matrices to solve the spectral problem given by (20). Therefore (20) is not integrable in the Lax pair sense. We also tried 3×3 matrices for U0 and V0, as we did for the NLSE with cubic nonlinearity, but we were unable to find the other matrices necessary to solve the spectral problem.


## 3.1.1. Nonlinear Schrödinger Equation, NLSE

If we seta3=a4=a5=0 in (8) we get the standard NLSE or GPE with damping or gaining and potential, which can be written as
(14)iψt+a1ψxx+a2|ψ|2ψ+iΓψ+Vextψ=0.
This is one of the most well-known NLSEs in the optical [1], Bose-Einstein [2], and ocean wave communities [3]. The Lax pair for (14) with time-independent coefficient and without damping and potential was first found by Zakharov and Shabat [21] and was then widely used in the Darboux transformation to find exact solutions. The damping is given by Γ=a˙2/(2a2)-a˙1/(2a1) and the potential is given by Vext=c12a1-c˙1x. Equation (14) supports solitonic solutions and one such solution is given by
(15)ψ(x,t)=Acosh⁡[c3x+b]ei[c32∫a1dt-c1x],
where A=-4c32a1/a2 and b=2c3∫c1a1dt. Similar results are also found by Kumar et al. [22]. Al Khawaja found solution of (14) with quadratic potential, periodic potential, and different types of solution with linear potential.It is well-known that the NLSE or GPE with time-dependent coefficients is integrable with potential up to quadratic power inx [23, 24]. Here, we find the Lax pair for NLSE with linear potential Vext=c12a1-c˙1x; this is due to the fact that we derive the Lax pair for the NLSE from that of the higher-order NLSE with cubic-quintic nonlinearity by taking appropriate limit of a3=a4=a5=a6=0 simultaneously [see (21)]. The higher-order terms in higher-order NLSE restrict integrability to linear potentials. One of the original motivations of our work is to solve the most general NLSE with cubic-quintic nonlinearity and find the solutions for other well-known NLSEs as special cases. However, we find that general NLSE with cubic-quintic nonlinearity is not solvable in the Lax pair sense.


## 3.1.2. Hirota Equation, HE

If we seta5=0 in (8) we get the HE [25]; namely,
(16)iψt+a1ψxx+a2|ψ|2ψ+ia3ψxxx+ia4|ψ|2ψx+iΓψ+Vextψ=0.
In this case we find that the equation is integrable only if a4=3a2a3/a1, Γ=a˙4/(2a4)-a˙3/(2a3), and Vext=0. The one-soliton solution of the HE is given by
(17)ψ(x,t)=Acosh⁡[6x-b]e6i∫a1dt,
where A=-12a3/a4 and b=66∫a3dt. Similar result was also found by Dai and Zhang [26].


## 3.1.3. Sasa-Satsuma Equation, SSE

If we seta1=a2=0 in (8) we get the SSE [27]; namely,
(18)iψt+ia3ψxxx+ia4|ψ|2ψx+ia5ψ2ψx*+iΓψ+Vextψ=0.
We find that the SSE is integrable for arbitrary values of a3, a4, and a5 with a4≠2a5 and Γ=Vext=0. The one-soliton solution of the SSE is given by
(19)ψ(x,t)=iAcosh⁡(2c2λx-b),
where A=2c1c2λ3/c3 and b=8c23(t-1)λ3a3. Similar result was found by Gilson et al. [28].We further investigate the integrability of the general higher-order NLSE equation with cubic and quintic nonlinearities [29–31]. In nonlinear optics, the nonlinear Schrödinger equation with quintic nonlinearity becomes important as the intensity and the inverse width of the optical pulse increase. In the presence of damping or gaining and an external potential, in nonlinear optics the general NLSE with cubic and quintic nonlinearities can be written as
(20)iψt+a1ψxx+a2|ψ|2ψ+ia3ψxxx+ia4|ψ|2ψx+ia5ψ2ψx*+a6|ψ|4ψ+ia7|ψ|4ψx+ia8|ψ|2ψ2ψx*+iΓψ+Vextψ=0.
Here, as for (8), all the coefficients ai(t) turn out to be real. The terms with coefficients a6, a7, and a8 account for the quintic nonlinearity and are generally known as non-Kerr nonlinearity. All the other terms are similar to those for (8). Using the systematic search approach we find the Lax pair for (20), which is given by
(21)U0=(c1ic2eGψ-ic2e-Gψ*c1),V0=(v11v12v21v22),
where the matrix elements vij are given by
(22)v11=-ic5|ψ|2+c22a3(ψ*ψx-ψψx*),v12=eG[ic4c5c2ψ-c5c2ψx-ic2a3ψxxlllllll-i3c2(a4+a5)|ψ|2ψ-i3c2a7|ψ|4ψc4c5c2],v21=e-G[-ic4c5c2ψ*-c5c2ψx*+ic2a3ψxx*lllllll+i3c2(a4+a5)|ψ|2ψ*+i3c2a7|ψ|4ψ*c4c5c2],v22=ic5|ψ|2-c22a3(ψ*ψx-ψψx*).
Here c2=(2a5-a4)/6a3, c4=-3a6/a7, c5=c22[a1(a4+a5)-3a2a3]/(3a5), and G=c3+ic4x, and c1 and c3 are two arbitrary time-independent constants. Using (21) in (2) we get (20) with the following integrability conditions: a7=9a3a5a6/(3a2a3-a1a4+2a1a5), a8=2a7/3, Γ=c˙2/c2, and Vext=c42c5/c22-c˙4x. This Lax pair is a fake one as we could not find the other matrices to solve the spectral problem given by (20). Therefore (20) is not integrable in the Lax pair sense. We also tried 3×3 matrices for U0 and V0, as we did for the NLSE with cubic nonlinearity, but we were unable to find the other matrices necessary to solve the spectral problem.


## 4. Conclusion

We employ the Darboux transformation to solve some of the NLSEs with cubic nonlinearity and time-dependent coefficients. We investigate the integrability of NLSE with damping and linear potential and find that it is integrable only with three independent coefficients, namely, dispersion, Kerr nonlinearity, and third-order dispersion coefficients. The coefficients for the Raman and self-steepening terms depend on the other three coefficients. The damping and linear part of the potential arise due to the time variation of the coefficients. This equation supports solitonic solution for focusing nonlinearity. For defocusing nonlinearity diverging nonsolitonic solution exists. In both the cases the coefficient for the third-order dispersion can control the center of mass motion of the solutions. This finding can be used to control the motion of the soliton in the optical fiber without an external potential.Using the same technique we are able to find the solitonic solution for some of the well-known NLSEs, namely, Hirota, Sasa-Satsuma, and Nonlinear Schrödinger equations. The Hirota and Sasa-Satsuma equations are not integrable with a linear potential, though the former is integrable with time-dependent coefficients. We check that the solutions of these equations are similar to those already found in the literature.We also investigate the integrability conditions for the NLSE with cubic and quintic nonlinearities, derivative terms, and complex potential. In this case, only fake Lax pairs were found indicating the nonintegrability of the equation in the Lax pair sense.



---

*Source: 323591-2014-09-15.xml*