Abstract and Applied Analysis

Volume 2014 (2014), Article ID 650781, 12 pages

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

Research Article

Melnikov’s Criteria and Chaos Analysis in the Nonlinear Schrödinger Equation with Kerr Law Nonlinearity

Jiuli Yin, Liuwei Zhao, and Lixin Tian

Nonlinear Scientific Research Center, Faculty of Science, Jiangsu University, Zhenjiang, Jiangsu 212013, China

Received 27 January 2014; Accepted 30 April 2014; Published 10 June 2014

Academic Editor: Guanglu Zhou

Copyright © 2014 Jiuli Yin 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 dynamics of the nonlinear Schrödinger equation with Kerr law nonlinearity with two perturbation terms are investigated. By using Melnikov method, the threshold values of chaotic motion under periodic perturbation are given. Moreover we also study the effects of the parameters of system on dynamical behaviors by using numerical simulation. The numerical simulations, including bifurcation diagram of fixed points, chaos threshold diagram of system in three-dimensional space, maximum Lyapunov exponent, and phase portraits, are also plotted to illustrate theoretical analysis and to expose the complex dynamical behaviors. In particular, we observe that the system can leave chaotic region to periodic motion by adjusting controller e, amplitude , and frequency of external forcing which can be considered a control strategy, and when the frequenciesy and approach the maximum frequency of disturbance, the system turmoil intensifies and control intensity increases.

1. Introduction

Fiber-optic signal transmission has a very wide application in the real life, which makes our lives more convenient and quicker. This is known to all. The fiber-optic signal plays an important role in real life. Optical solitary waves have been the subject of intense current research, which is motivated by their important applications in the areas of high-capacity fiber telecommunications and all-optical switches due to their capability of propagating over a long distance without attenuation and changing their shapes [1].

As is well known, the fiber-optic signal plays an important role in real life. It is the optical soliton propagation to maintain shape, amplitude and velocity constant light pulses for a long time. The use of optical solitons can achieve ultra-long-distance, large-capacity optical communications.

Fiber-optic signals make it seem that the signal propagation cannot exist in pure environment. It is always influenced by external environmental perturbations. It seems that chaos may be unavoidable with external perturbations and has been observed in many practical applications such as engineering, biology, industry, and production. Besides, many other systems with external perturbations have been widely investigated by using analytic methods and numerical simulations [2–5].

However, there is little to be subject to external disturbances as a model of optical signals. In fact, the fiber-optic signal propagation will be affected by lots of interferences in the process of real propagation. Moreover, the dynamic characteristics of nonlinear system in two-frequency interference will be richer and more complicated. We therefore think it is worthwhile to undertake a detailed discussion for system (1) in order to point out which range of parameter corresponds to certain behavior. Consequently, the two-frequency disturbance has more practical meaning.

One of the famous fiber-optic models is the external periodic perturbations of the nonlinear Schrödinger equation with Kerr law nonlinearity [6]

where is the third-order dispersion and and are the nonlinear dispersion. Here    () are real parameters and is a positive constant, , where the new variables , () denote the amplitude and the frequency of the parametric excitation, respectively. Here and () are real parameters. For this reason, ,   ().

Equation (1) describes the propagation of optical solitons in nonlinear optical fibers that exhibits a Kerr law nonlinearity. Equation (1) has an important application in various fields such as semiconductor materials, optical fiber communications, plasma physics, fluid and solid mechanics [6–12], and many other applications, can readily be described by the model, or analogous ones. Equation (1) has received more attention [6–10, 13, 14]. For example, by the first integral method, Taghizadeh et al. obtained the new exact solutions of the perturbed NLSE with Kerr law nonlinearity [6]. Wiggins [12] used the ()-expansion method to construct traveling wave solutions for NLSEs. Masemola et al. [9] obtained new solutions of nonlinear partial differential equation by a direct algebraic method. Miao and Zhang [11] constructed exact solutions of perturbed nonlinear Schrödinger’s equation with Kerr law nonlinearity by the simplest equation method.

The study of perturbed nonlinear Schrödinger’s equation with Kerr law nonlinearity is of fundamental and even practical interest. Eminent characteristics of perturbed nonlinear Schrödinger’s equation with Kerr law nonlinearity have a rich content of nonlinear properties which are suitable for a detailed investigating of various dynamical states.

In this paper, we will study the chaotic behavior and control of (1) under two-frequency perturbations. We study the bifurcation and dynamical behaviors depending on the parameters by using bifurcation and chaos theories in [15–21] and numerical simulation. By applying Melnikov method [15, 20, 21], we prove the criterion of existence of chaos under periodic perturbation. Our interests have the following two points.

(i)What if there appears a system chaos under the two-frequency interference? If chaos appears, how does one design a controller to suppress chaos owing to the complex nonlinear item of system? Compared with the Duffing system [4, 5, 12], it is easy to see that there is no damping in our system. Once perturbed with external forcing, the system may easily tend to the chaotic state. Therefore, we will select the controller which has the same function with the damping and we consider the following equation:

where denotes the controller’s strength.(ii)We also study the effects of the parameters of system on dynamical behaviors by using numerical simulation. The numerical simulations, including bifurcation diagram of fixed points, chaos threshold diagram of system in three-dimensional space, maximum Lyapunov exponent, and phase portraits, are also plotted to illustrate theoretical analysis and to expose the complex dynamical behaviors. In particular, we observe that the system can leave chaotic region to periodic motion by adjusting controller , amplitude , and frequency of external forcing which can be considered a control strategy.

This paper is organized as follows. In Section 2, we briefly describe the fixed points and phase portraits for the two types of unperturbed system of (1). In Section 3, the conditions of existence of chaos under periodic perturbation resulting from the homoclinic bifurcations are performed. The numerical investigations are given in Section 4. The bifurcation analysis is given in Section 5. Finally, we give remark to conclude this paper in Section 6.

2. Theoretic Analysis of System

2.1. The NLSE with Kerr Law Nonlinearity Equation

Assume that (1) has traveling wave solutions in the form [6]

where is the propagation speed of a wave ().

Substituting (3) into (1), we get

where    (), , and are positive constants and prime meaning differentiation with respect to .

By virtue of [6], we have

where denote the parametric of linear and nonlinear terms. Equation (5) is the fiber-optic signal transmission system in ideal environment. But it seems that the signal propagation cannot exist in pure environment. It is always influenced by external environmental perturbations.

Without loss of generality in (5), we set . By the transformations , , (6) can be transformed into first-order nonautonomous equation:

2.2. Fixed Points and Phase Portraits for Unperturbed System

If , (6) is considered an unperturbed system and can be written as

The system (7) is a Hamiltonian system with Hamiltonian function

And the function

is called the potential function. By the analysis of the fixed points () and their stabilities for system (7), we can easily obtain the following results.

Lemma 1. When , the unperturbed system has three equilibrium points: two centers and and one saddle . The saddle is connected to itself by two symmetric homoclinic orbits as shown in Figure 1(b).

(a)

(b)

Figure 1: (a) The two-well potential function of the system; (b) the corresponding phase space portraits.

3. Melnikov Theoretic Analysis

In this section, we discuss the chaotic behaviors of (6). Melnikov theory has proved to be a simple, elegant, and successful alternative to characterizing the complex dynamics of multistable oscillators. This section, thus, develops a global analysis technique, known as Melnikov’s method, to find the necessary conditions for homoclinic bifurcation to occur. For a detailed derivation of Melnikov’s method, there are several texts of varying rigor and sophistication to which the reader is referred [12, 22–24].

In this section, we discuss the chaotic behaviors of (6) in which and are assumed to be small parameters. Transformations of and are done in order to apply the first-order perturbation scheme of the Melnikov theory. Hence, the system of (6) may be written as

3.1. Melnikov Criterion for Chaos

The unperturbed system for system (7) has homoclinic orbits. When the perturbation is added, the closed homoclinic orbits break and may have transverse homoclinic or heteroclinic orbits. By the Smale-Birkhoff Theorem [12, 25], the existence of such orbits results in chaotic dynamics. We therefore apply the Melnikov method to system (7) for finding the criteria of the existence of homoclinic or heteroclinic bifurcation and chaos.

The Melnikov method derives a function to describe the first order distance between perturbed stable and unperturbed manifolds. Suppose that the unperturbed homoclinic or heteroclinic orbit is written as . Satisfying the conditions for a double-well potential gives rise to a homoclinic orbit in the phase space for . The homoclinic trajectory can be found by setting . Solving for the resulting displacement and differentiating to determine velocity, the homoclinic trajectory is given as follows:

and then the Melnikov function for system (7) can be given by

where is the cross-section time of the Poincaré map and can be interpreted as the initial time of the forcing term.

Because it is difficult to give analytical expression of , we will compute numerically in Section 5. We note that is a function of time from to . We therefore choose that the initial condition is and would be an odd function of time for the homoclinic orbit and an even function for heteroclinic orbit.

For the homoclinic orbits , the Melnikov function can be simplified as

where and are functions of the frequencies and , respectively.

Using the previous results and Melnikov’s theorem [26, 27], the following is stated: if and for some and some set of parameters, then horseshoes exist and chaos occurs [12, 25]. If has a simple zero and the corresponding critical parameter value is

where is a constant once is given, then in the system with fractional order displacement (11) the deterministic chaos may appear for certain parameter values which satisfy the relation

Remark 2. Using the Melnikov criterion for the appearance of the intersection between the perturbed and unperturbed separatrices, therefore, for fixed frequencies and , the system always produces Smale commutation of chaos.

3.2. Control of Chaos

Because the fiber-optic transmission system in the chaotic state is very sensitive to its initial condition and chaos often causes irregular behavior, chaos is undesirable. It is not hard to see that system (6) is similar to Duffing system, except the absence of damping in the former. Therefore, we will select the controller that has the same function with the damping.

Equation (2) is equivalent to the following system:

Equation (16) can be transformed into first-order nonautonomous equation:

where , , and .

Now, the Melnikov function for system (17) can be given by

and the Melnikov function can be simplified as

Thus, if

then there is a such that and , , and the following lemma can be obtained.

Lemma 3. The homoclinic bifurcation will occur at

Remark 4. For (20) the distance between the stable and unstable manifolds of the homoclinic point () is zero and the manifolds intersect transversely forming the transverse homoclinic orbits. The presence of such orbits implies that, for certain parameters (21) and the countable infinity of unstable periodic orbits, an uncountable set of bounded nonperiodic orbits and a dense orbit are the main characteristics of the chaotic motion.

4. Numerical Simulations

In this section, we give numerical simulations to support the theoretical results obtained in the previous sections and to find other new dynamics.

The interesting problem is to analyze the parameter regions for optical fiber signals’ stable propagation of the controlled system. The controlled fiber-optic transmission system has several parameters and each of them plays different and virtual roles in the system. We will analyze the influence on optic-fiber signals propagation of controlled system (21) when the parameter of system changes with the fixed controller.

Let us study the intersections of the invariant manifolds of the saddle point. It is known that these intersections are the necessary conditions for the existence of chaos. Since the Melnikov function theory measures the distance between the perturbed stable manifold and unstable manifolds, a homoclinic tangency will occur when a real solution can be found for some time such that the function has a simple zero. This means that only necessary conditions for the appearance of strange attractors are obtained from Poincaré-Melnikov-Arnold analysis, and therefore one always has the chance of finding sufficient conditions for the elimination of even transient chaos. Then the general necessary condition for which the invariant manifolds intersect is given by

We give numerical simulations to support the theoretical results as shown in Figures 2–6.

Figure 2: Chaotic threshold in space.

Figure 3: Chaotic threshold in space.

Figure 4: Chaotic threshold in space.

Figure 5: Chaotic threshold in space.

Figure 6: Chaotic threshold in space.

The chaotic threshold in space with , , , and is shown in Figure 2. The chaotic threshold in space with , , , and is shown in Figure 3. The chaotic threshold in space with , , , and is shown in Figure 4. The chaotic threshold in space with , , , and is shown in Figure 5. The chaotic threshold in space with , , , and is shown in Figure 6. When the value of is below the surface, the system may be chaotic state. In order to show what happens to the solutions and attractors as one crosses these bifurcation surfaces, we chose the parameter values from the regions , and then the system (17) can exhibit chaos.

Remark 5. This implies that if is sufficiently small, the reduced equation (17) has transverse homoclinic orbits resulting in possible chaotic dynamic. With and constant, we study chaotic threshold as a function of only the frequency parameter . A typical plot of is shown in Figure 2. The qualitative form of this function remains the same as and are the values for which the potential is two-well. Another point is that, when the ratio tends to zero, this means that the external amplitude tends to infinity and then the Melnikov theory is not valid for these values.

5. Bifurcation Analysis

In this section, we give numerical simulations to support the theoretical results of control of chaos.

5.1. Perturbation and Controlled System Analysis

We will discuss the behaviors of the fiber-optic signal transmission under perturbation and we draw the bifurcation diagram of (6) in () plane and the corresponding maximal Lyapunov exponents in Figure 7. System (6) is integrated using the Runge-Kutta technique of order four to conduct numerical simulation. Numerical calculations have been made for the selected parameter values , , , , and with the initial conditions and . Since the equation is nonlinear, its solution therefore admits the possibility of periodic and chaotic orbits.

(a)

(b)

Figure 7: (a) Bifurcation diagram of system (6) in () plane . (b) Maximum Lyapunov exponents corresponding to (a).

In Figure 7, we can observe that the value of Lyapunov exponents is positive (a positive top Lyapunov exponent for a bounded attractor is usually a sign of chaos), so the system easily converts to chaos even if there is small perturbation. It indicates that the system is outside perturbed, as the optical fiber transmission signal is very vulnerable to the phenomenon of chaos. So an appropriate controller is needed to satisfy the practical applications of fiber-optic propagation. Such maps can be used to suppress chaotic dynamics of the system.

Next the corresponding numerical simulations are performed for the case of periodic perturbation of system (16); we select . In Figure 8(a) after a large band of chaotic regime for one can find a sequence of backward period-doubling bifurcations as a route to periodic motion after a periodic regime; another bifurcation takes place at a critical value where another large band of chaotic regime with small periodic window occurs. At and , the system displays periodic behavior after backward period-doubling bifurcations. Such maps can be used to suppress chaotic dynamics of the system.

(a)

(b)

(c)

Figure 8: (a) Bifurcation diagram of system (17) in () plane . (b) Maximum Lyapunov exponents corresponding to (a). (c) Enlarged view of (b).

From Figure 8(b) one can see that for smaller values of controller , the top Lyapunov exponent is positive. As increases, the top Lyapunov exponent changes from positive value to negative value, signifying the suppressing of homoclinic chaos motion. Beyond the threshold for onset of chaotic motion, there are some “periodic windows;” this can be the feature of the transient chaos. However, the larger control parameters are responsive for the decreasing signal amplitude. It is easy to see that the signal cannot propagate normally and might leak from the media, which is called escape.

Remark 6. According to the above analysis, the increasing of the controller’s coefficient makes the system stable, but escape occurs when crosses a certain value.

5.2. Parameters Analysis

Now we wish to find other interesting bifurcation structures and dynamics of system (17). The bifurcation parameters are considered in the following five cases:

(i)varying in range and fixing , , , , and for several values of ,(ii)varying in range and fixing , , , , and for several values of ,(iii)varying in range and fixing , , , , and for several values of ,(iv)varying in range and fixing , , , , and for several values of ,(v)varying in range and fixing , , , , and for several values of .

For cases (i) and (ii), we give bifurcation diagrams in () and () plane of system (17) for and and several values of in Figures 9 and 10, respectively. We observe that a wide chaotic region in Figure 9(a). The coexistences of periodic orbits and chaotic motions in chaotic windows and intermittent mechanism in Figure 9(b). Figure 9(c) shows that the onset of chaos after period-n and chaotic regions with complex periodic windows, Figure 9(d) also shows the onset of chaos and chaotic regions with small periodic windows. Comparing Figures 10(a)–10(d), we find that, as is increased, the complexity of dynamical behaviors is decreased, and as there is only a period-one orbit and there is no bifurcation (see Figure 10(d)).

(a)

(b)

(c)

(d)

Figure 9: Bifurcation diagrams of system (17) in () plane for case (iii): (a) , (b) , (c) , and (d) .

(a)

(b)

(c)

(d)

Figure 10: Bifurcation diagrams of system (17) in () plane for case (iii): (a) , (b) , (c) , and (d) .

For case (iii), the bifurcation diagrams in () plane and the corresponding maximum Lyapunov exponents are given in Figure 11, where we show that the maximum Lyapunov exponents are all negative in Figure 11(d) for which are corresponding to quasiperiodic solutions in Figure 11(b). Figures 11(a) and 11(c) exhibit the process of the quasiperiodic route to chaos and chaos converges to quasiperiodic orbit.

(a)

(b)

(c)

(d)

Figure 11: Bifurcation diagrams of system (17) in () plane for case (iii): (a) ; (b) . The maximum Lyapunov exponents: (c) corresponding to (a), (d) corresponding to (b).

For case (iv), the bifurcation diagrams in Figures 12(a)–12(d) in () plane for various values of are plotted for showing the effect of parameter on dynamical behaviors. When is increased a little, the chaotic behavior occurrence has changed and the amplitude had decreased a little in Figure 12(b). But, as is further increased, there is only the period-five orbit and a period-doubling bifurcation of period-one.

(a)

(b)

(c)

(d)

Figure 12: Bifurcation diagrams of system (17) in () plane for case (iii): (a) , (b) , (c) , and (d) .

For case (v), the bifurcation diagrams in Figures 13(a)–13(d) in () plane for various values of are plotted for showing the effect of parameter on dynamical behaviors. We show that the chaotic behavior occurs alternately in smaller region of and the chaotic motions behavior occurs alternately in bigger region as a is increased. Therefore, we find that the parameter plays a very important role for dynamical behaviors from comparing. It can be considered as an control strategy of chaos by adjusting the parameter .

(a)

(b)

(c)

(d)

Figure 13: Bifurcation diagrams of system (17) in () plane for case (iii): (a) , (b) , (c) , and (d) .

6. Conclusions

In this paper, we have investigated the behaviors of the nonlinear Schrödinger equation with Kerr law nonlinearity with two perturbation terms and find many complex and interesting dynamical behaviors by using analytic and numerical methods. We conclude that chaos occurs easily due to the absence of damping in the system. This phenomenon will cause the distortion in the process of information transmission. One can add a controller to suppress the chaos. The efficiency of this controller was demonstrated. The complex fiber-optic transmission system of the perturbed NLSE with Kerr law nonlinearity was controlled. What is more is that we discussed the sensitivity to be controlled and found the practical parameters regions. Some of the results obtained for other NLSEs can be extended correspondingly.

Conflict of Interests

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

Acknowledgment

This work is supported by the National Nature Science Foundation of China (no. 11101191).

References

Y. JiuLi, “A duoble optical solitary wave in a nonlinear Schrödinger-type equation,” Chinese Physics B, vol. 22, no. 6, Article ID 060205, 2013.

S. Özgül, M. Turan, and A. Yildirim, “Exact traveling wave solutions of perturbed nonlinear Schrödinger's equation (NLSE) with Kerr law nonlinearity,” Optik, vol. 123, no. 24, pp. 2250–2253, 2012.

Z. Jing and R. Wang, “Complex dynamics in Duffing system with two external forcings,” Chaos, Solitons and Fractals, vol. 23, no. 2, pp. 399–411, 2005.

Z. Jing, J. Huang, and J. Deng, “Complex dynamics in three-well duffing system with two external forcings,” Chaos, Solitons and Fractals, vol. 33, no. 3, pp. 795–812, 2007.

R. Wang, J. Deng, and Z. Jing, “Chaos control in duffing system,” Chaos, Solitons and Fractals, vol. 27, no. 1, pp. 249–257, 2006.

N. Taghizadeh, M. Mirzazadeh, and A. Mahmoodirad, “Application of Kudryashov method for high-order nonlinear Schrödinger equation,” Indian Journal of Physics, vol. 87, no. 8, pp. 781–785, 2013.

Z.-Y. Zhang, Z.-H. Liu, X.-J. Miao, and Y.-Z. Chen, “New exact solutions to the perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity,” Applied Mathematics and Computation, vol. 216, no. 10, pp. 3064–3072, 2010.

H. Moosaei, M. Mirzazadeh, and A. Yildirim, “Exact solutions to the perturbed nonlinear Schrödinger's equation with kerr law nonlinearity by using the first integral method,” Nonlinear Analysis: Modelling and Control, vol. 16, no. 3, pp. 332–339, 2011.

P. Masemola, A. H. Kara, and A. Biswas, “Optical solitons and conservation laws for driven nonlinear Schrödingers equation with linear attenuation and detuning,” Optics and Laser Technology, vol. 45, no. 1, pp. 402–405, 2013.

A. R. Shehata, “The traveling wave solutions of the perturbed nonlinear Schrödinger equation and the cubic-quintic Ginzburg Landau equation using the modified (G′/G)-expansion method,” Applied Mathematics and Computation, vol. 217, no. 1, pp. 1–10, 2010.

X.-J. Miao and Z.-Y. Zhang, “The modified (G′/G)-expansion method and traveling wave solutions of nonlinear the perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 11, pp. 4259–4267, 2011.

S. Wiggins, Introduction to Applied Nonlinear Dynamical Systems and Chaos, vol. 2 of Texts in Applied Mathematics, Springer, New York, NY, USA, 1990.

W.-X. Ma and M. Chen, “Direct search for exact solutions to the nonlinear Schrödinger equation,” Applied Mathematics and Computation, vol. 215, no. 8, pp. 2835–2842, 2009.

J. Belmonte-Beitia and G. F. Calvo, “Exact solutions for the quintic nonlinear Schrödinger equation with time and space modulated nonlinearities and potentials,” Physics Letters A: General, Atomic and Solid State Physics, vol. 373, no. 4, pp. 448–453, 2009.

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

M. Lakshmanan and K. Murali, Chaos in Nonlinear Oscillators: Controlling and Synchronization, World Scientific, 1996.

E. Ott, Chaos in Dynamical Systems, Cambridge University Press, 2002.

S. N. Rasband, Chaotic Dynamics of Nonlinear Systems, John Wiley & Sons, New York, NY, USA, 1990.

J. M. T. Thompson and H. B. Stewart, Nonlinear Dynamics and Chaos, John Wiley & Sons, New York, NY, USA, 1986.

S. Wiggins, Global Bifurcation and Chaos: Analytical Methods, Springer, 1988.

Y. Uyaroğlu and S. Emiroğlu, “Passivity-based chaos control and synchronization of the four dimensional Lorenz-Stenflo system via one input,” Journal of Vibration and Control, 2013.

J. Guckenheimer and P. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, vol. 42 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1983.

D. W. Jordan and P. Smith, Nonlinear Ordinary Differential Equations: An Introduction for Scientists and Engineers, Oxford University Press, New York, NY, USA, 2007.

A. H. Nayfeh and B. Balachandran, Introduction to Nonlinear Dynamics: Analytical, Computational, and Experimental Methods, John Wiley & Sons, New York, NY, USA, 1995.

J. Guckenheimer and P. Homel, Nonlinear Oscillations, Dynamical System, and Bifurcation of Vector Fields, Springer, 1992.

S. C. Stanton, B. P. Mann, and B. A. M. Owens, “Melnikov theoretic methods for characterizing the dynamics of the bistable piezoelectric inertial generator in complex spectral environments,” Physica D: Nonlinear Phenomena, vol. 241, no. 6, pp. 711–720, 2012.

V. K. Melnikov, “On the stability of a center for time-periodic perturbations,” Trudy Moskovskogo Matematičeskogo Obščestva, vol. 12, pp. 3–52, 1963 (Russian).