Mathematical Problems in Engineering
Volume 2015 (2015), Article ID 408586, 10 pages
http://dx.doi.org/10.1155/2015/408586

Research Article

Bifurcation Analysis and Solutions of a Higher-Order Nonlinear Schrödinger Equation

Yi Li,1,2 Wen-rui Shan,1,2 Tianping Shuai,1,2 and Ke Rao1,2

¹State Key Laboratory of Information Photonics and Optical Communications, Beijing University of Posts and Telecommunications, Beijing 100876, China
²School of Science, Beijing University of Posts and Telecommunications, Beijing 100876, China

Received 2 October 2014; Revised 8 December 2014; Accepted 8 December 2014

Academic Editor: Reza Jazar

Copyright © 2015 Yi Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The purpose of this paper is to investigate a higher-order nonlinear Schrödinger equation with non-Kerr term by using the bifurcation theory method of dynamical systems and to provide its bounded traveling wave solutions. Applying the theory, we discuss the bifurcation of phase portraits and investigate the relation between the bounded orbit of the traveling wave system and the energy level. Through the research, new traveling wave solutions are given, which include solitary wave solutions, kink wave solutions, and periodic wave solutions.

1. Introduction

In the past decades, communication systems have scored a great growth of the transmission capacity. Due to the undamped and unchanged characteristics in a far distance, optical solitons are the focus of many research groups during the past decades and stand a good chance to be the main information carriers in telecommunications in the future. Ignoring optical losses, the wave dynamics of nonlinear pulse propagation in a monomode fiber is described by the nonlinear Schrödinger equation [1, 2], which accounts for the group velocity dispersion and self-phase modulation. To increase the bit rate, it is often desirable to use shorter femtosecond pulses. However, when short pulses are considered, the equation can no longer represent the propagations of light pulses in fibers because higher-order dispersion terms and the non-Kerr nonlinearity effects cannot be neglected. This phenomenon can be expressed by a higher-order nonlinear Schrödinger equation [3]:

where , , , , and are real constants. is a slowly varying envelope amplitude, represents the normalized retarded time (in the group velocity frame), and represents the normalized distance along the direction of propagation. comes from the group velocity dispersion (GVD). is proportional to the nonlinear index which originates from the Kerr effect. is the coefficients with the relevant work of the third-order dispersion. is related to self-steepening due to stimulated Raman scattering. The coefficient of the last term that is proportional to has its origin in the delayed Raman response . Generally speaking, can be estimated from the slope of the Raman gain and is defined as the first moment of the nonlinear response function [4]. In fact, should be an imaginary number, but many analytical studies have been done when is real, such as Painlevé property [5], inverse scattering transform [6], Hirota direct method, and conservation laws [7]. These researches verify its integrable nature and have obtained many exact wave solutions.

Laser spectroscopic techniques have been widely used in all fields of science. It can help us observe the physical processes in materials and molecules which occur on a femtosecond time scale by using ultrashort lasers. The pulses can also be applied in telecommunication and ultrafast signal routing systems. Research indicates that non-Kerr nonlinear effects begin to have some effects when the pulse width becomes narrower and the intensity of the incident light field becomes stronger. The influence is described by the NLS family of equations with nonlinear terms [8]. The nonlinearity due to fifth-order susceptibility can be obtained in many optical materials such as semiconductors and some transparent organic materials. Actually, it is also important to include some additional higher-order perturbation effects into the HNLS equation to analyze the solitary wave solution in a non-Kerr nonlinear medium.

In this paper, with the aid of Mathematica, we study the new traveling wave solutions for a higher-order NLS equation that contains the non-Kerr nonlinear terms, which describes propagation of very short pulses in highly nonlinear optical fibers by using different elliptic functions. The bifurcation theory method is widely used to solve differential equations [9–12]. By using this method of dynamical systems, we obtain the explicit expressions of the bounded traveling wave solutions for the equation and investigate the relation between the bounded orbit of the traveling wave system and the energy level . The new solutions correspond to the orbits on phase portraits and they include solitary, kink, and periodic wave solutions. Note that the existence of solitary wave solutions depends essentially on the model coefficients and therefore on the specific nonlinear features of the medium.

2. Bifurcation and Phase Portraits

We consider the higher-order NLS equation with non-Kerr term [13]:

where , , , , , , , and are real constants. , , and represent the coefficients of quintic non-Kerr nonlinearities. The quintic nonlinearities arise from the expansion of the refractive index of the light pulse. The polarizations induced through the susceptibilities give the cubic and quintic (non-Kerr) terms in a nonlinear Schrödinger equation. When , (2) reduces to (1).

Assume that (2) has the form of the exact solution

where is the real-valued function and the parameters of , , , and are real constants to be determined later.

Substituting (3) into (2) and removing the exponential term, we change (2) into the form

The real and imaginary parts of (4), respectively, are

Integrating (6), we can get

where is a constant. Equations (5) and (7) can be reduced to an equation if and

By solving (8), we get

or

Then (5) and (7) reduce to the following planar dynamic system:

where , , and .

Obviously, the above system (11) has the first integral

We suppose that

is the coefficient matrix of the linearized system (11) at an equilibrium point and

is the Jacobian determinant. By the bifurcation theory of planar dynamical system, we know that if , then the equilibrium point is a saddle point; if and , then it is a center point; if and , then it is a node; if and Poincaré index of the equilibrium point is , then it is a cusp point. By using the above facts to do qualitative analysis, we have the following.

(1)If and or , , , and , then system (11) has only one equilibrium point . It is easy to find . So it is a saddle point (see Figure 1).(2)If , , , and , then system (11) has only one equilibrium point . It is easy to find that and . So it is a center point (see Figure 2).(3)If , , , and , then system (11) has three equilibrium points and . , , and Poincaré index of is equal to zero. So is a saddle point; are cusp points (see Figure 3).(4)If , , , and , then system (11) has three equilibrium points and . , , , and Poincaré index of is equal to zero. So is a center point; are cusp points (see Figure 4).(5)If and , then system (11) has three equilibrium points and . , , and . So is a center point; are saddle points (see Figure 5).(6)If and , then system (11) has three equilibrium points and . , , and . So is a saddle point; are center points (see Figure 6).(7)If , , , and , then system (11) has five equilibrium points , , and . , , , and . So and are saddle points; are center points. Furthermore, if , then (see Figure 7); if , then (see Figure 8); if , then (see Figure 9).(8)If , , , and , then system (11) has five equilibrium points , , and . , , , , and . So and are center points; are saddle points (see Figure 10).

Figure 1: The phase portraits of (11) when and or , , , and .

Figure 2: The phase portraits of (11) when , , , and .

Figure 3: The phase portraits of (11) when , , , and .

Figure 4: The phase portraits of (11) when , , , and .

Figure 5: The phase portraits of (11) when and .

Figure 6: The phase portraits of (11) when and .

Figure 7: The phase portraits of (11) when , , , and .

Figure 8: The phase portraits of (11) when , , , and .

Figure 9: The phase portraits of (11) when , , , and .

Figure 10: The phase portraits of (11) when , , , and .

For a fixed , the curve

is called a level curve with energy level [14]. Obviously, each orbit of (12) is a branch of certain energy curve. For convenience, we name the orbit as the orbit with energy level .

To facilitate further analysis, we investigate the relation between the bounded orbit of (12) and the energy level .

Put

If , , and , it is easy to obtain the five extreme points of as follows:

Let

Therefore we can easily draw the graphics of the function in Figures 11, 12, and 13.

Figure 11: The phase portraits of (11) when , , , and .

Figure 12: The phase portraits of (11) when , , , and .

Figure 13: The phase portraits of (11) when , , , and .

Observe that the energy curve is equivalent to the curve defined by . According to the above analysis, we have the following.

Case 1. Assume that , , , and ; system (11) has three saddle points , , and two center points . There exist two homoclinic orbits with energy level to the saddle point (corresponding to red closed lines in the left half-plane and the right half-plane in Figure 7) and two families of periodic orbits with energy level () around the centers (corresponding to orange closed lines in the left half-plane and the right half-plane in Figure 7) which lie inside of above two homoclinic orbits, respectively. Besides, system (11) has a singular closed orbit with energy level to the closed curve (corresponding to the green closed lines in Figure 7) and a family of periodic orbits with energy level () (corresponding to the family of periodic orbits enclosing the equilibrium points and ) which lie inside of above singular closed orbit. It implies that for (11) there exist three families of periodic wave solutions, two solitary wave solutions, and two kink wave solutions.

Case 2. Assume that , , , and ; system (11) has three saddle points , , and two center points . There exist a singular closed orbit with energy level to the closed curve (corresponding to the red closed lines in Figure 8) and two families of periodic orbits with energy level () around the centers (corresponding to green closed lines in the left half-plane and the right half-plane in Figure 8) which lie inside of above singular closed orbit. It means that there exist two families of periodic wave solutions and four kink wave solutions.

Case 3. Assume that , , , and ; system (11) has three saddle points , , and two center points . There exist two homoclinic orbits with energy level at the saddle points (corresponding to red closed lines in the left half-plane and in the right half-plane in Figure 9) and two families of periodic orbits with energy level () around the centers (corresponding to blue closed lines in the left half-plane and the right half-plane in Figure 9) which lie inside of above two homoclinic orbits, respectively.

3. Traveling Wave Solutions of the HNLS Equation

(1) Solitary, Kink, and Periodic Wave Solutions When , , , and . (i) If , as is seen in Figure 7, there are two symmetric homoclinic orbits connected at the saddle point . In -plane the expressions of the homoclinic orbits are given as

where

Substituting (19) into and integrating them along the homoclinic orbits, noticing that and , we get two solitary wave solutions:

which correspond to the family of homoclinic orbits in the left half-plane and in the right half-plane shown in Figure 7.

(ii) If , as is seen in Figure 7, there are two heteroclinic orbits connected at the saddle points and . In -plane the expressions of the heteroclinic orbits are given as

where

Substituting (22) into and integrating them along the heteroclinic orbits, noticing that and , we get two kink wave solutions:

which correspond to the family of heteroclinic orbits and shown in Figure 7.

(iii) If , as is seen in Figure 7, there are two families of periodic orbits inside the homoclinic orbits. In -plane the expressions of the periodic orbits are given as ():

where

Substituting (25) into and integrating them along the periodic orbits, noticing that and , we get two families of periodic wave solutions:

which correspond to two families of periodic orbits inside the homoclinic orbits in the right half-plane and in the left half-plane shown in Figure 7.

When the energy level , the above periodic solutions tend to the solitary wave solutions (21).

(iv) If , as is seen in Figure 7, there is a family of periodic orbits enclosing the equilibrium points , and . In -plane the expressions of the periodic orbits are given as ():

where

Substituting (28) into and integrating them along the periodic orbits, noticing that and , we get a family of periodic wave solutions:

which correspond to a family of periodic orbits enclosing three equilibrium points shown in Figure 7.

When the energy level , the above periodic solutions also tend to the solitary wave solutions (21).

(2) Kink and Periodic Wave Solutions When , , , and . (i) If , as is seen in Figure 8, there are four heteroclinic orbits. In -plane the expressions of the heteroclinic orbits are given as

Substituting (31) into and integrating them along the heteroclinic orbits, noticing that and , we get four kink wave solutions:

which correspond to , , , and in Figure 8.

(ii) If , as is seen in Figure 8, there are two families of periodic orbits inside the homoclinic orbits. Thus we have the same solutions as (27).

(3) Solitary and Periodic Wave Solutions When , , , and . (i) If , as is seen in Figure 9, there are two symmetric homoclinic orbits connected at the saddle points and . In -plane the expressions of the homoclinic orbits are given as

where and satisfy (23).

Substituting (33) into and integrating them along the homoclinic orbits, noticing that and , we get two solitary wave solutions:

(ii) If , as is seen in Figure 9, there are two families of periodic orbits inside the homoclinic orbits. Thus we have the same solutions as (27).

We will draw the figures of some solutions under the special conditions. Figure 14 is a bell solitary wave solution of the field given by (21) when the energy level . Figure 15 shows the kink wave solution of the field given by (21) when the energy level . Figure 16 is a plot for the periodic wave solution of the field given by (27) when the energy level . Figure 17 is a periodic wave solution of the field given by (30) when the energy level . Figure 18 shows periodic wave solutions of the field given by (27) when . We can see that when , or we say , these periodic wave solutions will be reduced back to zero. When , or we say , the periodic solutions tend to the solitary wave solutions (21). Figure 19 is a periodic wave solution of the field when changes. When the energy level , the above periodic solutions also tend to the solitary wave solutions (21). Figures 20, 21, 22, 23, 24, and 25 show the relationship near critical points between , and , , and , respectively.

Figure 14: when , , , , and .

Figure 15: when , , , , and .

Figure 16: when , , , , , and .

Figure 17: when , , , , , and .

Figure 18: when , , and .

Figure 19: when , , and .

Figure 20: when and .

Figure 21: when and .

Figure 22: when and .

Figure 23: when and .

Figure 24: when and .

Figure 25: when and .

4. Conclusion

By using the bifurcation theory method of dynamical systems, we successfully obtain 10 phase portraits for the corresponding dynamic system of (2). Through analysing three of these phase portraits, we get new traveling wave solutions, including solitary wave solutions, kink wave solutions, and periodic wave solutions. The solutions are new and have not been investigated. By the dependent variable transformations, four linear forms for (2) can also be obtained. We will continue to consider the problem by using Hirota method in the future.

Conflict of Interests

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

Acknowledgments

The authors express their sincere thanks to the editor and reviewers for their valuable comments. This work has been supported by the National Natural Science Foundation of China under Grant no. 11272023, by the Open Fund of State Key Laboratory of Information Photonics and Optical Communications (Beijing University of Posts and Telecommunications) under Grant no. IPOC2013B008, and by the Fundamental Research Funds for the Central Universities of China under Grant no. 2011BUPTYB02.

References

G. P. Agrawal, Nonlinear Fiber Optics, Academic Press, San Diego, Calif, USA, 1995.

A. Hasegawa and Y. Kodama, Solitons in Optical Communications, Oxford University Press, Oxford, UK, 1995.

Y. Kodama and A. Hasegawa, “Nonlinear pulse propagation in a monomode dielectric guide,” IEEE Journal of Quantum Electronics, vol. 23, no. 5, pp. 510–524, 1987.

G. P. Agrawal, Nonlinear Fiber Optics, Academic Press, New York, NY, USA, 2006.

K. Porsezian and K. Nakkeeran, “Optical solitons in presence of Kerr dispersion and self-frequency shift,” Physical Review Letters, vol. 76, no. 21, pp. 3955–3958, 1996.

S. Palacios, A. Guinea, J. Fernández-Díaz, and R. D. Crespo, “Dark solitary waves in the nonlinear Schrödinger equation with third order dispersion, self-steepening, and self-frequency shift,” Physical Review E, vol. 60, no. 1, pp. R45–R47, 1999.

W.-J. Liu, “Soliton interaction in the higher-order nonlinear Schrödinger equation investigated with Hirota’s bilinear method,” Physical Review E, vol. 77, Article ID 066605, 2008.

R. Radhakrishnan, A. Kundu, and M. Lakshmanan, “Coupled nonlinear Schrödinger equations with cubic-quintic nonlinearity: integrability and soliton interaction in non-Kerr media,” Physical Review E—Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics, vol. 60, no. 3, pp. 3314–3323, 1999.

Z. Zhang, F.-L. Xia, and X.-P. Li, “Bifurcation analysis and the travelling wave solutions of the Klein-Gordon-Zakharov equations,” Pramana, vol. 80, no. 1, pp. 41–59, 2013.

M. Song, B. S. Ahmed, and A. Biswas, “Topological soliton solution and bifurcation analysis of the Klein-Gordon-Zakharov equation in -dimensions with power law nonlinearity,” Journal of Applied Mathematics, vol. 2013, Article ID 972416, 7 pages, 2013.

M. Song, Z. Liu, E. Zerrad, and A. Biswas, “Singular solitons and bifurcation analysis of quadratic nonlinear Klein-Gordon equation,” Applied Mathematics & Information Sciences, vol. 7, no. 4, pp. 1333–1340, 2013.

Z. S. Wen, “Bifurcation of solitons, peakons, and periodic cusp waves for theta-equation,” Nonlinear Dynamics, vol. 77, no. 1-2, pp. 247–253, 2014.

A. Choudhuri and K. Porsezian, “Higher-order nonlinear Schrödinger equation with derivative non-Kerr nonlinear terms: a model for sub-10-fs-pulse propagation,” Physical Review A, vol. 88, Article ID 033808, 2013.

J. B. Li and H. H. Dai, On the Study of Singular Nonlinear Traveling Wave Equations: Dynamica System Approach, Science Press, Beijing, China, 2007.