# Nonlinear Dynamics and Exact Traveling Wave Solutions of the Higher-Order Nonlinear Schrodinger Equation with Derivative Non-Kerr Nonlinear Terms.

1. IntroductionNonlinear Schrodinger (NLS) equation is one of the most important nonlinear models in mathematical physics and has many applications in nonlinear optics, plasma physics, condensed matter physics, photonics, and Bose-Einstein condensates. In particular, in the studies on optical fibers, the NLS equation is very important. As is well known, the higher-order nonlinear Schrodinger (HNLS) equation describes the propagation of picosecond or femtosecond optical pulse in fibers. Therefore, the study of the higher-order nonlinear Schrodinger (HNLS) equation is the hotspot in the study of nonlinear scientific fields.

In this paper, we consider the following higher-order nonlinear Schrodinger equation with derivative non-Kerr nonlinear terms [1]:

[mathematical expression not reproducible], (1)

where E(z, t) is the slowly varying envelope of the electric field, the subscripts z, t are the spatial and temporal partial derivatives in related time coordinates, and [a.sub.1], [a.sub.2], [a.sub.3], [a.sub.4], and [a.sub.5] are the real parameters related to the group velocity, self-phase modulation, third-order dispersion, self-steeping, and self-frequency shift arising from stimulated Raman scattering, respectively, [2]. The terms related to coefficients [a.sub.6], [a.sub.7], [a.sub.8] in (1) represent the quintic non-Kerr nonlinearities. The investigation of this equation has raised great interest due to its wide range of applications. In [1], Choudhuri and Porsezian investigated the periodic wave solutions, the bright and dark solitary wave solutions of (1). In [3], Choudhuri and Porsezian investigated the Dark-in-the-Bright (DITB) solitary wave solution of (1). In [4], Choudhuri and Porsezian have studied the modulational instability (MI) of (1) with forth-order dispersion in context of optics and presented an analytical expression for MI gain to show the effects of non-Kerr nonlinearities and higher-order dispersions on MI gain spectra, and so on. However, we notice that the dynamics of the traveling wave solutions of (1) have not be studied. It is meaningful and necessary to consider the dynamical behavior of (1) and to find all possible exact solutions of (1). In the present paper, we will use the dynamical system method to investigate the travelling wave solutions of the higher-order nonlinear Schrodinger equation with derivative non-Kerr nonlinear terms.

Firstly, to investigate the existence of travelling wave solution of HNLS equation in presence of non-Kerr terms, we begin with scaling the variables of (1) in the form

E = [b.sub.1][psi], z = [b.sub.2]x, t = [b.sub.3][tau] (2)

and choosing [b.sub.1], [b.sub.2], and [b.sub.3] such that the coefficients corresponding to group velocity dispersion (GVD), self-phase modulation (SPM), and third-order dispersion (TOD) become unity. Thus (1) becomes

[mathematical expression not reproducible], (3)

where [mathematical expression not reproducible]. We have chosen [b.sub.1] = [square root of [absolute value of ([a.sup.3.sub.1]/[a.sup.2.sub.2][a.sup.3.sub.3])]], [b.sub.2] = [a.sup.2.sub.3]/[a.sup.3.sub.1] and [b.sub.3] = [a.sub.3]/[a.sub.1] in writing (3).

To obtain the exact travelling wave solutions of (3), we consider the travelling wave solutions of the following form:

[psi](x, [tau]) = P([xi])[e.sup.i[eta]], [xi] = vx + [tau], [eta] = kx - c[tau], (4)

where k, c, and v are travelling wave parameters. Substituting (4) into (3), canceling [e.sup.in], and separating the real and imaginary parts, we have

[mathematical expression not reproducible]. (5)

Equating the two equations, we get the following conditions:

c = 3[c.sub.1] + 2[c.sub.2] - 3/6([c.sub.1] + [c.sub.2]) = 5[c.sub.4] + 4[c.sub.5] - 5[c.sub.3]/10[c.sub.4] + 12[c.sub.5], k = (1 - 3c)(v - 2c + 3[c.sup.2]) - [c.sup.2] + [c.sup.3], (6)

with constraint relations

[c.sub.3] = 3/5, [c.sub.4] = 3[c.sub.1]/5, [c.sub.5] - [c.sub.2]/2. (7)

As a result of the freedom of these parameters which is consistency, under condition (6) and (7), (5) is simplified to the following equation:

[P.sup.n] - (v - 2c + 3[c.sup.2])P + 3[c.sub.1] + 2[c.sub.2]/3 [P.sup.3] + 5[c.sub.4] + 4[c.sub.5]/5 [p.sup.5] = 0. (8)

Denote that [d.sub.1] = v - 2c + 3[c.sup.2], [d.sub.2] = -(3[c.sub.1] + 2[c.sub.2])/6, and [d.sub.3] = -(5[c.sub.4] + 4[c.sub.5])/15. Thus, (8) has the following form:

[P.sup.n] - [d.sub.1]P - 2[d.sub.2][P.sup.3] - 3[d.sub.3][P.sup.5] = 0, (9)

which corrsponds to the two-dimensional Hamiltonian system:

dP/d[xi] = y, dy/d[xi] = [d.sub.1]P + 2[d.sub.2][P.sup.3] + 3[d.sub.3][P.sup.5] (10)

with the Hamiltonian

H(P, y) = 1/2 [y.sup.2] - [d.sub.1]/2 [P.sup.2] - [d.sub.2]/2 [P.sup.4] - [d.sub.3]/2 [P.sup.6] = h. (11)

According to the Hamiltonian, we can get all kinds of phase portraits in the parametric space. Because the phase orbits defined the vector fields of system (10) and determined all their travelling wave solutions of (3), by investigating the bifurcations of phase portraits of system (10), we can seek the travelling wave solutions of (3) [5-11]. The detailed calculation procedure can be found in the technical appendix. The rest of this paper is built up as follows. In Section 2, we give all phase portraits of system (10) and discuss the dynamics of phase portraits of system (10). In Section 3, according to the dynamics of the phase orbits of system (10) given by Section 2, we obtain all possible bounded travelling wave solutions of (3). Finally, a conclusion is given in Section 4.

2. Nonlinear Dynamics of Phase Portraits of System (10)

In this section, we consider the bifurcations of the phase orbits of system (10) on the phase plane (P, y) as the parameters [d.sub.1], [d.sub.2], [d.sub.3] are changed.

Firstly, we consider the distribution of the equilibrium points of system (10). Obviously, the zeros of the function f(P) = [d.sub.1] + 2[d.sub.2][P.sup.3] + 3[d.sub.3][P.sup.5] appeared and the second term of system (10) is the abscissas of equilibrium points of system (10) on the phase plane (P, y). We write that [mathematical expression not reproducible] are the same sign, according to Vieta theorem, that is, the relationship between the root and coefficient of quadratic equation with one unknown; we have the following proposition.

Proposition 1. (1) For the case of [d.sub.3] = 0, because [d.sub.2] = (-3[c.sub.1] - 2[c.sub.2])/6 and [d.sub.3] = (-3[c.sub.1] - 2[c.sub.2])/15 are the same sign, [d.sub.2] = [d.sub.3] =0. system (10) has only one equilibrium point U0(0,0). (2) Forthe case of [d.sub.3] [not equal to] 0, when [DELTA] < 0 or [d.sub.1] [d.sub.3] [greater than or equal to] 0, system (10) has only one equilibrium point [U.sub.0](0, 0). (3) For the case of [d.sub.3] [not equal to] 0, when [delta] > 0 and [d.sub.1][d.sub.3] < 0, system (10) has three equilibrium points, [U.sub.0](0, 0) and [U.sub.1] [+ or -] ([+ or -][square root of ([f.sub.1], 0)] or [U.sub.2] [+ or -] ([+ or -][square root of ([f.sub.1], 0)]).

Then, we consider the type of the equilibrium points of system (10). Let M([U.sub.i], 0) (i = 0, ..., 2) be the Jacobin matrix of system (10) at an equilibrium point ([U.sub.i], 0) and let J([U.sub.i], 0) (i = 0, ..., 2) be the Jacobin determinant. Thus, we have

Trace ([U.sub.i], 0) = 0, (i = 0, ..., 2), J(0, 0) = -[d.sub.1], J([+ or -][square root of ([f.sub.1], 0)] = -4[f.sub.1][square root of ([DELTA])] J([+ or -][square root of ([f.sub.2], 0)] = -4[f.sub.1][square root of ([DELTA])]. (12)

For H(P, y) = (1/2)[y.sup.2] - ([d.sub.1]/2)[P.sup.2] - ([d.sub.2]/2)[P.sup.4] - ([d.sub.3]/2)[P.sup.6] = h, we write that

[mathematical expression not reproducible]. (13)

By the bifurcation theory of dynamical system, for an equilibrium point of a planar integrable system, the equilibrium point is a saddle point when J < 0. The equilibrium point is a center point when J > 0 and Trace M([U.sub.i], 0) = 0. The equilibrium point is a cusp when J = 0 and the Poincare index of the equilibrium point is 0.

Finally, by using the above discussion, the bifurcations of phase portraits of system (10) for the case of [d.sub.3] = 0, [d.sub.3] > 0, and [d.sub.3] < 0 are shown in Figures 1-3 with the aid of Maple.

3. Exact Travelling Wave Solutions of (3)

In this section, corresponding to all phase orbits given by Section 2, through qualitative analysis and the Jacobian elliptic functions [12], we discuss the exact travelling wave solutions of (3). Since only bounded travelling waves are meaningful to a physical model, here we just pay our attention to the bounded solutions of (3).

3.1. The Case of [d.sub.3] = 0 (See Figure 1). Suppose that [d.sub.1] <0 (see Figure 1(b)). When h [member of] (0, +[infinity]), there exists a family of periodic travelling wave solutions which correspond to a family of periodic orbits of (10). It follows from (11) that [y.sup.2] = -[d.sub.1](-(2h/[d.sub.1]) - [P.sup.2]). By using the first equation of (10), we have

[P.sub.1][xi] = [square root of (-2h/[d.sub.1])] cos ([square root of (-[d.sub.1] [xi])] (14)

Thus, we obtain the periodic travelling wave solutions of (3) as follows:

[[psi].sub.1](x,[tau]) = [square root of (-2h/[d.sub.1])] cos ([square root of (-[d.sub.1])](vx + [tau]))[e.sup.i(kx-c[tau])]. (15)

3.2. The Case of [d.sub.3] > 0 (See Figure 2)

(1) Suppose That [d.sub.1] < 0 (See Figure 2(b)). When h [member of] (0, [h.sub.1]), there exists a family of periodic travelling wave solutions which correspond to a family of periodic orbits of (10). It follows from (11) that [y.sup.2] = [d.sub.3]([[gamma].sub.1] - [P.sup.2])([P.sup.2] - [[gamma].sub.2])([[gamma].sub.3] - [P.sup.2]), where [[gamma].sub.1] > [[gamma].sub.3] > 0 > [[gamma].sub.3] . Therefore, we have

[mathematical expression not reproducible]. (16)

Thus, we obtain the periodic travelling wave solutions of (3) as follows:

[mathematical expression not reproducible]. (17)

(2) Suppose That [d.sub.1] <0 (See Figure 2(b)). When h = [h.sub.1], there exists a kink wave solution and an anti-kink wave solution which correspond to two heteroclinic orbits of (10). It follows from (11) that [y.sup.2] = [d.sub.3][([f.sub.1] - [P.sup.2]).sup.2]([P.sup.2] - [[gamma].sub.4]), where [[gamma].sub.4] = -([d.sub.2] + 2[square root of ([DELTA])]/3[d.sub.3], [f.sub.1] > 0 > [[gamma].sub.4]. Therefore, we have

[P.sub.3]([xi]) = [+ or -]([f.sub.1] + 2[f.sub.1]([f.sub.1] - [[gamma].sub.4])/[[gamma].sub.4] - 2[f.sub.1] + [[gamma].sub.4] cosh(2[square root of ([d.sub.3][f.sub.1]([f.sub.1] - [[gamma].sub.4])[xi])]). (18)

Thus, we obtain a kink wave solution and an anti-kink wave solution of (3) as follows:

[[psi].sub.3](x, [tau]) = [+ or -]([f.sub.1] + 2[f.sub.1]([f.sub.1] - [[gamma].sub.4])/[[gamma].sub.4] - 2[f.sub.1] + [[gamma].sub.4] cosh(2[square root of ([d.sub.3][f.sub.1]([f.sub.1] - [[gamma].sub.4])(vx + [tau])]) x [e.bar]i(kx-c[tau]) (19)

3.3. The Case of [d.sub.3] < 0 (See Figure 3)

3.3.1. The Travelling Wave Solutions Corresponding to Figure 3(a)

(1) Suppose That [d.sub.2] < 0, [d.sub.1] < [d.sup.2.sub.2]/4[d.sub.3] or [d.sub.2] < 0, [d.sub.1] < [d.sup.2.sub.2]/3[d.sub.3]. When h [member of] (0, +[infinity]), there exists a family of periodic travelling wave solutions which correspond to a family of periodic orbits of (10). Therefore, we have

[mathematical expression not reproducible]. (20)

Thus, we obtain the periodic travelling wave solutions of (3) as follows:

[mathematical expression not reproducible], (21)

where [A.sup.2.sub.1] = [([[gamma].sub.1] - [[alpha].sub.1]).sup.2] + [[beta].sup.2.sub.1], [B.sup.2.sub.1] = [[alpha].sup.2.sub.1] + [[beta].sup.2.sub.1], [[alpha].sub.1] = ([[gamma].sub.2] + [[gamma].sub.3])/2, and [[beta].sub.1] = -[([[gamma].sub.2] - [[gamma].sub.3]).sup.2]/4.

(2) Suppose That [d.sub.2] < 0, [d.sup.2.sub.2]/4[d.sub.3] < [d.sub.1] [less than or equal to] 0. When h = [h.sub.3], there exists a family of periodic travelling wave solutions which correspond to a family of periodic orbits of (10). It follows from (11) that [y.sup.2] = -[d.sub.3][([P.sup.2] - [f.sub.1]).sup.2]([[gamma].sub.4] - [P.sup.2]), where [[gamma].sub.4] = -([d.sub.1] + 2[square root of ([DELTA])])/3[d.sub.3], [[gamma].sub.4] > 0 > [f.sub.1], and [h.sub.3] = ([square root of [[gamma].sub.4]], 0) = 2[d.sub.1][d.sub.2][d.sub.3] - 2[d.sup.3.sub.2] + 2[[DELTA].sup.3/2])/54[d.sup.2.sub.3] Therefore, we have

[P.sub.5]([xi]) = ([f.sub.i] + 2[f.sub.1]([[gamma].sub.4] - [f.sub.1])/2[f.sub.1] - [[gamma].sub.4] - [[gamma].sub.4] cos (2[square root of ([d.sub.3][f.sub.1]([[gamma].sub.4] - [f.sub.1])[xi])]). (22)

Thus, we obtain the periodic travelling wave solutions of (3) as follows:

[[psi].sub.5](x, [tau]) = [([f.sub.1] + 2[f.sub.1]([[gamma].sub.4] - [f.sub.1])/2[f.sub.1] - [[gamma].sub.4] - [[gamma].sub.4] cos(2[square root of [d.sub.3][f.sub.1]([[gamma].sub.4] - [f.sub.1])(vx + [tau])]).sup.1/2] [e.sup.i(kx-c[tau])]. (23)

(3) Suppose That [d.sub.2] < 0, [d.sup.2.sub.2]/4[d.sub.3] < [d.sub.1] [less than or equal to] 0. When h [member of] (0, [h.sub.3]), there exists a family of periodic travelling wave solutions which correspond to a family of periodic orbits of (10). It follows from (11) that [[gamma].sub.2] = -[d.sub.3]([P.sup.2] - [[gamma].sub.1])([[gamma].sub.2] - [P.sup.2])([P.sup.2] - [[gamma].sub.3]), where [[gamma].sub.2] > 0 > [[gamma].sub.3] > [[gamma].sub.1]. Therefore, we have

[mathematical expression not reproducible]. (24)

Thus, we obtain the periodic travelling wave solutions of (3) as follows:

[mathematical expression not reproducible]. (25)

(4) Suppose That [d.sub.2] < 0, [d.sup.2.sub.2]/4[d.sub.3] < [d.sub.1] [less than or equal to] 0. When h [member of] ([h.sub.3], +[infinity]), there exists a family of periodic travelling wave solutions which correspond to a family of periodic orbits of

(10). It follows from (11) that [y.sup.2] = -[d.sub.3]([P.sup.2] - [[gamma].sub.1])([[gamma].sub.2] - [P.sup.2])([P.sup.2] - [[gamma].sub.3]), where [[gamma].sub.2] > 0 and [[gamma].sub.1] and [[gamma].sub.3] are conjugate complex number. Therefore, we have

[mathematical expression not reproducible]. (26)

Thus, we obtain the periodic travelling wave solutions of (3) as follows:

[mathematical expression not reproducible], (27)

where [A.sup.2.sub.2], = [([[gamma].sub.2] - [[alpha].sub.2]).sup.2] + [[beta].sup.2.sub.2], [B.sup.2.sub.2] = + [[alpha].sup.2.sub.2] + [[beta].sup.2.sub.2], [[alpha].sub.2] = ([[gamma].sub.1] + [[gamma].sub.3])/2, and [[beta].sub.2] = -[([[gamma].sub.1] - [[gamma].sub.3]).sup.2]/4.

(5) Suppose That [d.sub.2] < 0, [d.sub.1] = [d.sup.2.sub.2]/4[d.sub.3]. Equation (3) has a family of periodic travelling wave solutions with the same solutions as (27).

3.3.2. The Travelling Wave Solutions Corresponding to Figure 3(b)

(1) Suppose That [d.sub.1] > 0. When h [member of] ([h.sub.2], 0), there exists two family of periodic travelling wave solutions which correspond to two family of periodic orbits of (10). It follows from (11) that [y.sup.2] = -[d.sub.3]([P.sup.2] - [[gamma].sub.1])([[gamma].sub.2] - [P.sup.2])([P.sup.2] - [[gamma].sub.3]), where [[gamma].sub.2] > [[gamma].sub.3] > 0 > [[gamma].sub.1].

Therefore, we have

[mathematical expression not reproducible]. (28)

Thus, we obtain the periodic travelling wave solutions of (3) as follows:

[mathematical expression not reproducible]. (29)

(2) Suppose That [d.sub.1] > 0. When h = 0, there exists a pair of solitary wave solutions which correspond to two homoclinic orbits of (10). It follows from (11) that [y.sup.2] = [d.sub.3][P.sup.2])([[gamma].sub.6] - [P.sup.2]), where [[gamma].sub.5] > [[gamma].sub.6] > 0, [[gamma].sub.5] = (-[d.sub.2] + [square root of ([d.sup.2.sub.2] - 4[d.sub.1][d.sub.3]))]/2[d.sub.3], and [[gamma].sub.6] = (-[d.sub.2] -[square root of ([d.sup.2.sub.2] - 4[d.sub.1][d.sub.3])])/2[d.sub.3]. Therefore, we have

[P.sub.9] ([xi])

= [+ or -][(2[[gamma].sub.5][[gamma].sub.6]/[[gamma].sub.5] + [[gamma].sub.6] + ([[gamma].sub.5] - [[gamma].sub.6]) cosh (2[square root of [d.sub.3][[gamma].sub.5] [[gamma].sub.6][xi]]))).sup.1/2]. (30)

Thus, we obtain the solitary wave solutions of (3) as follows:

[mathematical expression not reproducible] (31)

(3) Suppose That [d.sub.1] > 0. When h [member of] (0, +[infinity]), h [member of] (0, [h.sub.3]), h = [h.sub.3], and h [member of] ([h.sub.3], +[infinity]), (3) has a family of periodic travelling wave solutions with the same solutions as (21), (23), (25), and (27), respectively.

Through the approach of dynamical system, we have studied the exact travelling wave solutions of the higher-order nonlinear Schrodinger equation with derivative non-Kerr nonlinear terms. Firstly, through the travelling wave transformation, the higher-order nonlinear Schrodinger equation with derivative non-Kerr nonlinear terms is reduced a planer Hamiltonian system. Then with the aid of Maple, the bifurcations of phase portraits of the planer Hamiltonian system are given. By studying the bifurcations of phase portraits of the planer Hamiltonian system, we obtain exact travelling wave solutions of the higher-order nonlinear Schroodinger equation with derivative non-Kerr nonlinear terms, which contain solitary wave solutions, kink and anti-kink wave solutions, and periodic travelling wave solutions. Note that our solutions in this paper are different from the given ones in previous references [1-4]. We have the hyperbolic function solutions, trigonometric function solutions, and the Jacobian elliptic function solutions. From the above discussions, obviously the dynamical system method is very powerful method to seek exact travelling wave solutions for nonlinear travelling wave equations. The method is concise, direct, and effective which reduces the large amount of calculations. It is a good method which allows us to solve complicated nonlinear evolution equations in mathematical physics.

4. Conclusion

By the results of Sections 2 and 3 and considering (6) and (7), we obtain the following main conclusion of this paper.

Theorem 2. Suppose that [d.sub.1], [d.sub.3], and [d.sub.3] are given by Section 1 and [f.sub.1], [f.sub.3], [A.sub.1], [B.sub.1], and [[gamma].sub.1],(i = 1, ..., 6) are defined by Sections 2 and 3. The higher-order nonlinear Schrodinger equation with derivative non-Kerr nonlinear terms has the following 9 classes exact bounded travelling wave solutions.

(1) When [d.sub.1] < 0, [d.sub.2] = 0, and [d.sub.3] = 0,

[[psi].sub.1](x, [tau]) = [square root of (- 2h/[d.sub.1])] cos ([square root of (-[d.sub.1])](vx + [tau]))[e.sup.i(kx-c[tau])]. (32)

(2) When [d.sub.1] < 0, [d.sub.2] > 0, and [d.sub.3] > 0,

[mathematical expression not reproducible]. (33)

(3) When [d.sub.1] < 0, [d.sub.2] > 0, and [d.sub.3] > 0,

[mathematical expression not reproducible]. (34)

(4) When [d.sub.1] < [d.sup.2.sub.2]/4[d.sub.3], [d.sub.2] < 0, [d.sub.3] < 0 or [d.sub.1] < [d.sup.2.sub.2]/3[d.sub.3], [d.sub.2] < 0, [d.sub.3] < 0,

[mathematical expression not reproducible]. (35)

(5) When [d.sup.2.sub.2]/4[d.sub.3] < [d.sub.1] [less than or equal to] 0, [d.sub.2] < 0, and [d.sub.3] < 0,

[[psi].sub.5](x, [tau]) = [([f.sub.1] + 2[f.sub.1]([[gamma].sub.4] - [f.sub.1])/2[f.sub.1] - [[gamma].sub.4] - [[gamma].sub.4] cos (2[square root of ([d.sub.3][f.sub.1]([[gamma].sub.4] - [f.sub.1]))](vx + [tau])).sup.1/2] x [e.sup.i(kx-c[tau])]. (36)

(6) When [d.sup.2.sub.2]/4[d.sub.3] < [d.sub.1] [less than or equal to] 0, [d.sub.2] < 0, and [d.sub.3] < 0,

[mathematical expression not reproducible].(37)

(7) When [d.sup.2.sub.2]/4[d.sub.3] < [d.sub.1] [less than or equal to] 0, [d.sub.2] < 0, and [d.sub.3] < 0,

[mathematical expression not reproducible]. (38)

(8) When [d.sub.1] > 0, [d.sub.2] < 0, and [d.sub.3] < 0,

[mathematical expression not reproducible]. (39)

(9) When [d.sub.1] > 0, [d.sub.2] < 0, and [d.sub.3] < 0,

[[psi].sub.9](x, [tau]) = [+ or -] [(2[[gamma].sub.5][[gamma].sub.6]/[[gamma].sub.5] + [[gamma].sub.6] + ([[gamma].sub.5] - [[gamma].sub.6]) cosh (2[square root of ([d.sub.3][[gamma] .sub.5][[gamma].sub.6])](vx + [tau]))).sup.1/2]x [e.sup.i(kx-c[tau])]. (40)

Appendix

Here we describe the dynamical system method for finding traveling wave solutions of nonlinear wave equations. A (n + 1)-dimensional nonlinear partial differential equation is given as follows:

[mathematical expression not reproducible]. (A.1)

The main steps of the dynamical system method are as follows.

Step 1 (reduction of (A.1)). Making a transformation M(t, [x.sub.1], [x.sub.2], ..., [x.sub.n]) = [phi]([xi]), [xi] = [[summation].sup.n.sub.i=1] [k.sub.i][x.sub.i] - ct, (A.1) can be reduced to a nonlinear ordinary differential equation:

D([xi], [phi], [[phi].sub.[xi]], [[phi].sub.[xi][xi]], [[phi].sub. [xi][xi][xi]], ...) = 0 (A.2)

where [k.sub.i] are nonzero constant and c is the wave speed. Integrating several times for (A.2), if it can be reduced to the following second-order nonlinear ordinary differential equation,

E([phi], [[phi].sub.[xi]], [[phi].sub.[xi][xi]]) = 0, (A.3)

then let [[phi].sub.[xi]] = d[phi]/d[xi] = y, and (A.3) can be reduced to a two dimensional dynamical system:

d[phi]/d[xi] = y, dy/d[xi] = f([phi], y), (A.4)

where f([phi], y) is an integral expression or a fraction. If f([phi], y) is a fraction such as f([phi], y) = F([phi], y)/g([phi]) and g([[phi].sub.s]) = 0, dy/d[xi] does not exist when 0 = 0S. Then we will make a transformation d[xi] = g([phi])d[[zeta]; thus system (A.4) can be rewritten as

d[phi]/d[zeta] = g([phi])y, dy/d[zeta] = F([phi], y), (A.5)

where [zeta] is a parameter. If (A.1) can be reduced to the above system (A.4) or (A.5), then we can go on to the next step.

Step 2 (discussion of bifurcations of phase portraits of system (A.4)). If system (A.4) is an integral system, systems (A.4) and (A.5) can be reduced the differential equation:

dy/d[phi] = f([phi], y)/y, dy/d[phi] = F([phi], y)/g([phi])y = f([phi], y)/y, (A.6)

and then systems (A.4) and (A.5) have the same first integral (that is Hamiltonian) as follows:

H ([phi], y) = h, (A.7)

where [zeta] is an integral constant. According to the first integral, we can get all kinds of phase portraits in the parametric space. Because the phase orbits defined the vector fields of system (A.4) (or system (A.5)) and determined all their travelling wave solutions of (A.1), we can investigate the bifurcations of phase portraits of system (A.4) (or system (A.5)) to seek the travelling wave solutions of (A.1). Usually, a periodic orbit always corresponds to a periodic wave solution; a homoclinic orbit always corresponds to a solitary wave solution; a heteroclinic orbit (or so called connecting orbit) always corresponds to kink (or anti-kink) wave solution. When we find all phase orbits, we can get the value of h or its range.

Step 3 (calculation of the first equation of system (A.4)). After h is determined, we can get the following relationship from (A.7):

y = y([phi],h); (A.8)

that is, dtyldS, = y(Q, h). If the expression (A.8) is an integral expression, then substituting it into the first term of (A.4) and integrating it, we obtain

[mathematical expression not reproducible] (A.9)

where [phi](0) and 0 are initial constants. Usually, the initial constants can be taken by a root of (A.8) or inflection points of the travelling waves. Taking proper initial constants and integrating (A.9), through the Jacobian elliptic functions, we can obtain the exact travelling wave solutions of (A.1).

From the above description of "three-step method," we can see that solutions of (A.1) can be obtained by studying and solving the dynamical system simplified by (A.1). Therefore, this approach is called dynamical system method. The different nonlinear wave equations correspond to different dynamical systems. The different dynamical systems correspond to different travelling wave solutions. This is the whole process of the dynamical system method.

http://dx.doi.org/10.1155/2016/7405141

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] A. Choudhuri and K. Porsezian, "Higher-order nonlinear Schroodinger equation with derivative non-Kerr nonlinear terms: a model for sub-10-fs-pulse propagation," Physical Review A, vol. 88, no. 3, Article ID 033808, 2013.

[2] F. Yan andH. Liu, "Bifurcation of phase and exact traveling wave solutions of a higher-order nonlinear Schroodinger equation," International Journal of Bifurcation and Chaos, vol. 22, no. 5, Article ID 1250121, 11 pages, 2012.

[3] A. Choudhuri and K. Porsezian, "Dark-in-the-Bright solitary wave solution of higher-order nonlinear Schrodinger equation with non-Kerr terms," Optics Communications, vol. 285, no. 3, pp. 364-367, 2012.

[4] A. Choudhuri and K. Porsezian, "Impact of dispersion and non-Kerr nonlinearity on the modulational instability of the higher-order nonlinear Schrodinger equation," Physical Review A, vol. 85, no. 3, Article ID 033820, 2012.

[5] J. Li, Singular Nonlinear Travelling Wave Equations: Bifurcations and Exact Solutions, Science Press, Beijing, China, 2013.

[6] J. Li and H.-H. Dai, On the Study of Singular Nonlinear Travelling Wave Equations: Dynamical System Approach, Science Press, Beijing, China, 2007.

[7] H. Liu and J. Li, "Symmetry reductions, dynamical behavior and exact explicit solutions to the Gordon types of equations," Journal of Computational and Applied Mathematics, vol. 257, pp. 144-156, 2014.

[8] J.-B. Li and T.-L. He, "Exact traveling wave solutions and bifurcations in a nonlinear elastic rod equation," Acta Mathematicae Applicatae Sinica--English Series, vol. 26, no. 2, pp. 283-306, 2010.

[9] H. Li, K. Wang, and J. Li, "Exact traveling wave solutions for the Benjamin-Bona-Mahony equation by improved Fan sub-equation method," Applied Mathematical Modelling, vol. 37, no. 14-15, pp. 7644-7652, 2013.

[10] S. Xie, L. Wang, and Y. Zhang, "Explicit and implicit solutions of a generalized Camassa-Holm Kadomtsev-Petviashvili equation," Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 3, pp. 1130-1141, 2012.

[11] J. Li, "Exact explicit travelling wave solutions for (n +1)-dimensional Klein-Gordon-Zakharov equations," Chaos, Solitons & Fractals, vol. 34, no. 3, pp. 867-871, 2007.

[12] P. F. Byrd and M. D. Friedman, Handbook of Elliptic Integrals for Engineers and Scientists, Springer, Berlin, Germany, 1971.

Heng Wang, (1) Longwei Chen, (1) Hongjiang Liu, (2) and Shuhua Zheng (1)

(1) College of Statistics and Mathematics, Yunnan University of Finance and Economics Kunming, Yunnan 650221, China

(2) City and Environment College, Yunnan University of Finance and Economics Kunming, Yunnan 650221, China

Correspondence should be addressed to Longwei Chen; 1187411801@qq.com

Received 9 July 2015; Revised 8 November 2015; Accepted 10 December 2015

Academic Editor: Giacomo Innocenti

Caption: FIGURE 1: The bifurcations of phase portraits of (10) for [d.sub.3] = 0. (a) [d.sub.2] = 0, [d.sub.1] > 0. (b) [d.sub.2] = 0, [d.sub.1] < 0.

Caption: FIGURE 2: The bifurcations of phase portraits of (10) for [d.sub.3] > 0. (a) [d.sub.1] [greater than or equal to] 0, [d.sub.2] > 0 or [d.sub.1] > [d.sup.2.sub.2]/3[d.sub.3], [d.sub.2] > 0. (b) [d.sub.1] <0, [d.sub.2] > 0.

Caption: FIGURE 3: The bifurcations of phase portraits of (10) for [d.sub.3] < 0. (a) [d.sub.1] [less than or equal to] 0, [d.sub.2] < 0 or [d.sub.1] < [d.sup.2.sub.2]/3[d.sub.3], [d.sub.2] < 0. (b) [d.sub.1] > 0, [d.sub.2] < 0.

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Wang, Heng; Chen, Longwei; Liu, Hongjiang; Zheng, Shuhua |

Publication: | Mathematical Problems in Engineering |

Date: | Jan 1, 2016 |

Words: | 5035 |

Previous Article: | Synchronization for a Class of Uncertain Fractional Order Chaotic Systems with Unknown Parameters Using a Robust Adaptive Sliding Mode Controller. |

Next Article: | A Hybrid Algorithm of Particle Swarm Optimization and Tabu Search for Distribution Network Reconfiguration. |

Topics: |