# New exact solutions for high dispersive cubic-quintic nonlinear Schrodinger equation.

1. IntroductionPropagation of short pulses in optical fibers is governed by the well-known nonlinear Schrodinger equation (NLS) [1]. In recent years, There have been extensive study and application of NLS. The main purpose of this paper is to discuss the traveling wave solutions for a class of high dispersive cubic-Quintic nonlinear Schrodinger equations describing the ultrashort light pulse propagation as in the following:

[E.sub.z] = i[[[beta].sub.2]/2][E.sub.tt] + i[[gamma].sub.1] [[absolute value of E].sup.2]E + [[[beta].sub.3]/6] [E.sub.ttt] + [[[beta].sub.4]/24][E.sub.tttt] - i[[gamma].sub.2] [[absolute value of E].sup.4]E, (1)

where E(z, t) is the slowly varying envelope of the electric field, [[beta].sub.2] is the parameter of the group velocity dispersion, [[beta].sub.3] and [[beta].sub.4] are, respectively, the third-order and fourth-order dispersions, and [[gamma].sub.1] and [[gamma].sub.2] are the nonlinearity coefficients. When the higher order terms are ignored, we obtain the NLS. However, for femtosecond light pulses, whose duration is shorter than 10 fs, the last three terms are not ignored. Equation (1) was derived by Palacios and Fernandez-Diaz [2]. Azzouzi et al. [3] by using the extended hyperbolic auxiliary equation method in getting the exact explicit solutions to (1). He et al. [4] find the exact bright, dark, and gray analytical nonautonomous soliton solutions of the generalized CQNLSE with spatially inhomogeneous group velocity dispersion (GVD) and amplification or attenuation by the similarity transformation method under certain parametric conditions.

We will study (1) by using the improved Fan subequation method. As a result, more types of exact solutions to (1) are obtained, which include solitons, kink solutions, and Jacobian elliptic function solutions with double periods. The rest of this paper is organized as follows. In Section 2, we give the mathematical framework of the improved method. In Section 3, we apply it to the generalized equation (1) for finding more exact solutions. Finally, some conclusions are given.

2. The Ansatz Solution and Fan Subequation Method

The integrability of a nonlinear equation can be studied by applying the Painleve analysis. It is widely believed that possession of the Painleve property is a sufficient criterion for integrability. Moreover, there exists another technique which basically consists of expressing the solution in terms of an amplitude and a phase function as an approach to find exact solutions of nonlinear evolution equations. We will make use of this formalism looking for exact solution of (1) such as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)

where [phi]([xi]) is a real fuction and [xi] = [v.sub.0]z - vt. By inserting the expressions (2) into (1), and separating real and imaginary parts, we obtain

[l.sub.1][phi]'([xi]) + [l.sub.3][phi]'"([xi]) = 0, (3a)

[l.sub.1][phi]([xi]) + [l.sub.2][phi]"([xi]) + [l.sub.4][phi]""([xi]) + [[gamma].sub.1][[phi].sup.3]([xi]) - [[gamma].sub.2][[phi].sup.5]([xi]) = 0, (3b)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

Let [l.sub.1] = [l.sub.3] = 0, and we get

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)

Then, (3a) and (3b) become

[l.sub.0][phi]([xi]) + [l.sub.2][phi]"([xi]) + [l.sub.4][phi]""([xi]) + [[gamma].sub.1][[phi].sup.3]([xi]) - [[gamma].sub.2][[phi].sup.5]([xi]) = 0, (6)

We introduce auxiliary equation:

[phi]'[xi] [square root of ([c.sub.0] + [c.sub.1][phi]([xi]) + [c.sub.2][[phi].sup.2]([xi]) + [c.sub.3][[phi].sup.3([xi])] + [c.sub.4][[phi].sup.4]([xi]))], (7)

where [epsilon] = [+ or -] 1, which is known as Fan subequation method and proposed by Fan in [5]. This method is proposed to seek more types of exact solutions of nonlinear partial differential equations. Obviously, (7) is equivalent to the two-dimensional systems as follows:

d[phi]/d[xi] = y, dy/d[xi] = 1/2([c.sub.1] + 2[c.sub.1][phi] + 3[c.sub.3][[phi].sup.2] + 4[c.sub.4][[phi].sup.3]), (8)

which has the Hamiltonian function:

H([phi], y) = [y.sup.2] - ([c.sub.1][phi] + [c.sub.2][[phi].sup.2] + [c.sub.3][[phi].sup.3] + [c.sub.4][[phi].sup.4]) = [c.sub.0] (9)

One can easily find that [c.sub.0] corresponds to the Hamiltonian constant and (7) is equivalent to the Hamiltonian system (8). Thus, in order to search the exact solutions of (7) we need only to discuss (8). For a fixed [c.sub.0], (9) determines a set of orbits of (8). As [c.sub.0] varies, (9) defines different families of orbits of (8) which have different dynamical behavior. Below we will first study the bifurcation of phase portraits of (8) by making use of bifurcation method of dynamical systems and with the aid of the computer symbolic system Mathematica. Then according to the obtained bifurcation and the Hamiltonian function (9), we will gain many new exact solutions of (7) for all possible parameters [c.sub.j] [6].

Substituting (7) into (6), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10)

Setting all the coefficients of [phi]' (i = 0, 1, ..., 5) to zero, and solving the obtained algebraic equations, we find the following sets of solutions (I)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

and (II)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

and (III), [[gamma].sub.2] = 0:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

where [[omega].sub.0], v are any real number and v [not equal to] 0.

3. Exact Solutions of High Dispersive Cubic-Quintic Nonlinear Schrodinger Equation

3.1. Case (I) and Case (II). In this case, the Hamiltonian system (8) becomes

d[phi]/d[xi] = y, dy/d[xi] = [c.sub.2][phi] + 2[c.sub.4][[phi].sup.3], (14)

and the Hamiltonian function (9) reduces to

[H.sub.1] ([phi], y) = [y.sup.2] - [c.sub.2][[phi].sup.2] - [c.sub.4][[phi].sub.4] = [c.sub.0]. (15)

Now we discuss the bifurcations of phase portraits of (14). Obviously all the equilibrium points of (14) lie in the [phi]-axis and their abscissas are the real zeros of f([phi]) = [c.sub.2][phi] + 2[c.sub.4][[phi].sup.3]. Thus, we have the following proposition on the distribution of the equilibrium points of (14).

Proposition 1.

(1) For [c.sub.2][c.sub.4] < 0, (14) has three equilibria points at [E.sub.10](0,0), [E.sub.11]([[psi].sub.11], 0), and [E.sub.12]([[psi].sub.12], 0), where [[psi].sub.11] = -[square root of (-[c.sub.2]/(2[c.sub.4]))], [[psi].sub.12] = [square root of (-[c.sub.2]/(2[c.sub.4]))].

(2) For [c.sub.2][c.sub.4] [greater than or equal to] 0, (14) has a unique equilibrium at [E.sub.10](0, 0).

Using the bifurcation theory of dynamical systems [7-9], the ([c.sub.2], [c.sub.4])-plane was divided into four subregions:

D1: {[c.sub.2] > 0, [c.sub.4] > 0}; D2: {[c.sub.2] < 0, [c.sub.4] > 0}; D3: {[c.sub.2] < 0, [c.sub.4] < 0}; D4: {[c.sub.2] > 0, [c.sub.4] < 0}. (16)

The phase portraits of (14) are shown in Figure 1.

For the function defined by (15), we have

[h.sub.10] = [H.sub.1]([[psi].sub.10], 0) = 0, [h.sub.11] = [h.sub.12] = [H.sub.1]([[psi].sub.11], 0) = [H.sub.1]([[psi].sub.12], 0) = [c.sup.2.sub.2]/4[c.sub.4]. (17)

Let M([w.sub.e], [z.sub.e]) be the coefficient matrix of the system (14) at an equilibrium point ([w.sub.e], [z.sub.e]). Then, we have

J ([w.sub.e], 0) = det (M ([w.sub.e], 0)) = -([c.sub.2] + 6[c.sub.4][w.sup.2.sub.e]). (18)

By the theory of planar dynamical systems, we know that for an equilibrium point of a planar integrable system, if J < 0, then the equilibrium point is a saddle point; if J > 0 and Trace (M([w.sub.e])) = 0, then it is a center point; if J > 0 and [(Trace(M([w.sub.e],[z.sub.e]))).sup.2] - 4J([w.sub.e],[z.sub.e]) > 0, then it is a node; if J = 0 and the index of the equilibrium point is 0, then it is a cusp; otherwise, it is a high order equilibrium point.

Below, we will give explicit and exact solutions of (7) (also (6)). We always let [DELTA] = [c.sup.2.sub.2] - 4[c.sub.0][c.sub.4] and only pay attention to the bounded solutions of (6).

(1) Suppose that ([c.sub.2], [c.sub.4]) [member of] D4.

(a) For [c.sub.0] [member of] ([c.sup.2.sub.2]/(4[c.sub.4]), 0), (7) has periodic solutions:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (19)

where sn(x, k) and below cn(x, k) are Jacobian elliptic functions with modulus k [10]. The profiles of periodic solutions are shown in Figure 2.

Thus, we obtain the following solutions of (1):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (20)

(b) For [c.sub.0] = 0, we have solitary wave solutions of (7):

[[phi].sub.[+ or -] 2] = [+ or -] [square root of (-[c.sub.2]/[c.sub.4])] sech ([square root of [c.sub.2][xi]]. (21)

Thus, we obtain the following solutions of (1):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (22)

The profiles of solutions of (22) are shown in Figure 3.

Remark 2. By the expression of [c.sub.0], there always exists a [w.sub.0] such that [c.sub.0] = 0 if ([c.sub.2], [c.sub.4]) [member of] D4.

Because of the limitation of length, we omit the expression of E(z, i), beginning from here.

(c) For [c.sub.0] > 0, (6) has periodic solution:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (23)

(2) Suppose that ([c.sub.2], [c.sub.4]) [member of] D2.

(a) For [c.sub.0] [member of] (0, [c.sup.2.sub.2]/(4[c.sub.4])), system (14) has periodic solution:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (24)

(b) For [c.sub.0] = [c.sup.2.sub.2]/(4[c.sub.4]), there exists smooth kink wave solution of (14) as follows:

[[phi].sub.[+ or -] 5] = [+ or -] [square root of (-[c.sub.2]/2[c.sub.4])] tanh ([square root of (-[c.sub.2]/2)][xi]). (25)

(3) If ([c.sub.2], [c.sub.4]) [member of] D3, for [c.sub.0] > 0, we have the following periodic solution of (14):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (26)

(4) If ([c.sub.2], [c.sub.4]) [member of] D1, it is observed from Figure 1(a) that there is no bounded solution of system (14).

3.2. The Case (III). For this case, the Hamiltonian system (8) becomes

d[phi]/d[xi] = y, dy/d[xi] = f([phi]) = 1/2([c.sub.1] + 2[c.sub.2][phi] + 3[c.sub.3][[phi].sup.3], (27)

and the Hamiltonian function (9) reduces to

[H.sub.1] ([phi], y) = [y.sup.2] - [c.sub.1][phi] - [c.sub.2][[phi].sup.2] - [c.sub.3][[phi].sup.3] = [c.sub.0]. (28)

Similar to the previous discussion, we have the following proposition.

Proposition 3.

(1) For [c.sup.2.sub.2] - 3[c.sub.1][c.sub.3] > 0, (27) has two equilibria at [E.sub.21] ([[psi].sub.21], 0) and [E.sub.22] ([[psi].sub.22], 0), where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (29)

(2) For [c.sup.2.sub.2] - 3[c.sub.1][c.sub.3] = 0, (27) has a unique equilibrium at [E.sub.20]([[psi].sub.20], 0), where = [[psi].sub.20] = -[c.sub.2]/(3[c.sub.3]).

(3) For [c.sup.2.sub.2] - 3[c.sub.1][c.sub.3] < 0, (27) has no equilibrium.

Let [h.sub.2i] = [H.sub.1]([[psi].sub.2i], 0), i = 0, 1, 2, and notice that we need only to consider the case [c.sub.3] [greater than or equal to] 0 because of the invariance of (27) under the transformations [phi] [right arrow] -[psi], y [right arrow] -y, and [c.sub.3] [right arrow] -[c.sub.3].

(1) [c.sup.2.sub.2] - 3[c.sub.1][c.sub.3] > 0 and [c.sub.0] [member of] ([h.sub.21], [h.sub.22]). In this case,

[c.sub.3][[phi].sup.3] + [c.sub.2][[phi].sub.2] + [c.sub.1][phi] + [c.sub.0] = 0 (30)

has three mutually different real roots [[phi].sub.m] < [[phi].sub.l] < [[phi].sub.M]; thus,

[c.sub.3][[phi].sup.3] + [c.sub.2][[phi].sup.2] + [c.sub.1][phi] + [c.sub.0] = [c.sub.3] ([phi] - [[phi].sub.m]) ([phi] - [[phi].sub.l]) ([phi] - [[phi].sub.M]). (31)

Equation (7) has periodic wave solutions as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (32)

(2) [c.sup.2.sub.2] - 3[c.sub.1][c.sub.3] > 0 and [c.sub.0] = [h.sub.22]. In this case, [[phi].sub.22] is double root of (30); suppose that [[phi].sub.t] is other root of the equation, obviously [[phi].sub.t] < [[phi].sub.22], and we have a solitary wave solution of peak type of (7) as follows:

[[phi].sub.10] = [[phi].sub.t] + ([[phi].sub.22] - [[phi].sub.t])[tanh.sup.2] (1/2 [square root of ([c.sub.2] (([[phi].sub.22] - [[phi].sub.t])[xi])]). (33)

4. Conclusions

In this study, we apply bifurcation theory of dynamical systems and the Fan subequation method to investigate (1), and many new exact solutions have been obtained; most importantly, under more general conditions than [3], to the best of our knowledge, these solutions have not been reported in the literature. This method can help us find exact solutions of other types of nonlinear dispersion partial differential equations.

http://dx.doi.org/10.1155/2014/826746

Conflict of Interests

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

Acknowledgments

This work was supported by the National Natural Science Foundation of China (11261013 and 11361017), Natural Science Foundation of Guangxi (2012GXNSFAA053003 and 2013GXNSFAA019010), and the Foundation of Guangxi Key Lab of Trusted Software.

References

[1] G. P. Agrawal and C. Headley III, "Kink solitons and optical shocks in dispersive nonlinear media," Physical Review A, vol. 46, no. 3, pp. 1573-1577, 1992.

[2] S. L. Palacios and J. M. Fernandez-Diaz, "Black optical solitons for media with parabolic nonlinearity law in the presence of fourth order dispersion," Optics Communications, vol. 178, no. 4, pp. 457-460, 2000.

[3] F. Azzouzi, H. Triki, K. Mezghiche, and A. El Akrmi, "Solitary wave solutions for high dispersive cubic-quintic nonlinear Schrodinger equation," Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1304-1307, 2009.

[4] J. He, J. Zhang, M. Zhang, and C. Dai, "Analytical nonautonomous soliton solutions for the cubic-quintic nonlinear Schrodinger equation with distributed coefficients," Optics Communications, vol. 285, no. 5, pp. 755-760, 2012.

[5] E. Fan, "Uniformly constructing a series of explicit exact solutions to nonlinear equations in mathematical physics," Chaos, Solitons & Fractals, vol. 16, no. 5, pp. 819-839, 2003.

[6] Y. Xie, S. Tang, and D. Feng, "New exact solutions to the generalized KdV equation with generalized evolution," Pramana, vol. 78, no. 4, pp. 499-511, 2012.

[7] J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, Springer, New York, NY, USA, 1983.

[8] D. Luo, Bifurcation Theory and Methods of Dynamical Systems, World Scientific, Singapore, 1997.

[9] P. Lawrence, Differential Equations and Dynamical Systems, Springer, New york, NY, USA, 1991.

[10] P. F. Byrd and M. D. Friedman, Handbook of Elliptic Integrals for Engineers and Physicists, Springer, Berlin, Germany, 1954.

Yongan Xie and Shengqiang Tang

School of Mathematics and Computing Science, Guilin University of Electronic Technology, Guilin, Guangxi 541004, China

Correspondence should be addressed to Shengqiang Tang; tangsq@guet.edu.cn

Received 6 May 2014; Revised 8 July 2014; Accepted 18 July 2014; Published 5 August 2014

Academic Editor: Sheng Zhang

Printer friendly Cite/link Email Feedback | |

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

Author: | Xie, Yongan; Tang, Shengqiang |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 2703 |

Previous Article: | Several guaranteed descent conjugate gradient methods for unconstrained optimization. |

Next Article: | Spatial and temporal variation of annual precipitation in a river of the Loess Plateau in China. |

Topics: |