Research Article | Open Access
The Center Conditions and Bifurcation of Limit Cycles at the Degenerate Singularity of a Three-Dimensional System
We investigate multiple limit cycles bifurcation and center-focus problem of the degenerate equilibrium for a three-dimensional system. By applying the method of symbolic computation, we obtain the first four quasi-Lyapunov constants. It is proved that the system can generate 3 small limit cycles from nilpotent critical point on center manifold. Furthermore, the center conditions are found and as weak foci the highest order is proved to be the fourth; thus we obtain at most 3 small limit cycles from the origin via local bifurcation. To our knowledge, it is the first example of multiple limit cycles bifurcating from a nilpotent singularity for the flow of a high-dimensional system restricted to the center manifold.
About dynamical behavior of the trajectories of three-dimensional system, bifurcation of limit cycles is one of major concerns; particularly, for Hopf bifurcation of a nondegenerate equilibrium with a pair of pure imaginary roots and a negative one, many investigations have been carried out in the past decades, for example, [1–4] for the three-dimensional chaotic systems, [5–8] for the three-dimensional Lotka-Volterra systems, and  for general three-dimensional systems.
However, up till now, study on bifurcation of limit cycles from the degenerate singularity for high-dimensional nonlinear dynamical systems is hardly seen in published references. In this paper, we investigate the following three-dimensional systems: where . Obviously system (1) has the Jacobian matrix at the origin as follows: One notices that the origin has two zero eigenvalues and one negative eigenvalue and the block matrix is nilpotent. From the center manifold theorem, for the system (1), there exists the center manifold: with , , and more the flow on the center manifold is governed by namely, where . It is usually called the reduced system and system (1) is topological equivalent to system (4) in the vicinity of the origin . Thus by investigating system (4), we can solve Hopf bifurcation of the origin on center manifold of the three-dimensional system (1). At the same time, we also discuss the center-focus problem for the flow restricted to the center manifold, which closely relates to the maximum number of limit cycles bifurcating from the origin.
As far as investigation about limit cycles bifurcation of a nilpotent critical point in planar system is concerned, it is similar to nondegenerate case that detecting nilpotent center and calculating the focal value are needed . There exist some available and classical ways, for instance, Poincare return map , Lyapunov function , and the normal form theory . At the same time, some good results on the cyclicity were obtained [11, 12, 15, 16]. Recently in [17, 18], the authors gave an inverse integral factor method of calculating the quasi-Lyapunov constants of the three-order nilpotent critical point; it is convenient to compute the higher order focal values and solve the center-focus problem of the equilibrium. Here we extend this method's application to the three-dimensional system (1) and consider its specific example as follows: where .
The rest of this paper is organized as follows. In Section 2, the corresponding quasi-Lyapunov constants are computed and the center conditions on the center manifold are determined. In Section 3, the multiple local bifurcations at the origin for system (5) are investigated, three limit cycles from the origin are obtained, and it is proved at most three small limit cycles from the origin via local bifurcation. In this work, the system and problem are all considered for the first time.
2. Quasi-Lyapunov Constants and Center Conditions
In this part, we firstly investigate the singular point quantities of the origin. For the center manifold of system (5), one can determine the formal expression and obtain the same form as system (4): Furthermore, we give the definition of quasi-Lyapunov constants for system (6) and the way of computing them; more details can be found in [17–19].
Lemma 1 (Theorems and 8.7.2 in ). For system (6), any positive integer and a given number sequence , we can derive successively and uniquely the terms of the following formal series with the coefficients satisfying ,
And, if , is determined by the following recursive formula:
and, for any positive integer , is determined by the following recursive formula:
and, when or , we have let .
Particularly, by choosing appropriate and number sequence , we can make and then let
Definition 2. The in (11) is called the th quasi-Lyapunov constant of the origin of system (6), , and more if all the quasi-Lyapunov constants vanish, that is, , , then the origin of system (5) is a center on the local center manifold at the origin.
Lemma 3. For system (6), the th focal value at the origin of system (6) is algebraic equivalent to ; that is, for any positive integer , if and hold, then , where is a defined constant (which is given in [17, 18]).
Remark 4. According to Lemma 1, one sees that each in (10) is related to only the coefficients of the first degree terms of system (6). Here we determine the above only up to a tenth degree polynomial with respect to , as follows: where is a homogeneous polynomial in , of degree , respectively, (which can be seen in the Appendix); thus and in system (6) are two polynomials with degree 11. And more all , of in (9) are given definitely by the coefficients , , (the specific is available in Email address of the corresponding author); the in (10) is given by the following specific form:
Applying the powerful symbolic computation function of Mathematica system and the recursive formulas in Lemma 1, and from Remark 4, we obtain the first quantities as follows: In the above expression of each , , we have already let .
Particularly, in order to make , , we let and choose
Theorem 5. For the flow on center manifold of system (5), we get the first 4 quasi-Lyapunov constants of the origin as follows:
Theorem 6. For system (5), the origin is a center on the local center manifold if and only if the following condition is satisfied:
Proof. The proof of the necessity is easy; then we omit it. Now we prove the sufficient condition; this technique derives from the method of Darboux (also see [20–23]). Obviously, if in the conditions (17) holds, then system (5) has the corresponding form as follows:
And more we can figure out easily one algebraic invariant surface for system (18): ; in fact, there exists a polynomial , the cofactor of , such that . One can observe that is tangent to the center eigenspace, the -plane, at the origin. Thus it forms a local center manifold in a neighborhood of the origin. From , we substitute it into the first and second equations of the system defined by system (18) and we have the differential equations
which is a Hamiltonian system with Hamiltonian function: . Therefore the origin is a center for systems (18) or (19) as shown in Figure 1(a).
And if in the conditions (17) holds, then system (5) has the corresponding form as follows: Clearly, the above is also a first integral of system (20); thus the origin is a center for the flow of system (5) restricted to a center manifold.
Theorem 7. For the origin of system (5) as a weak focus on center manifold, the highest is the fourth order and its first 4 focal values are as follows: where, for the expression of , we have let , and for , for .
3. Multiple Limit Cycle Bifurcation of the System
In this section, we apply the focal values obtained in last section to discuss multiple bifurcation of the equilibrium and demonstrate there exist 3 and at most 3 limit cycles.
Theorem 8. For the flow on center manifold of (5), the origin is a 4th-order fine focus; that is, , if and only if
Proof. Firstly, we prove the necessity, from the expressions (16) in Theorem 5; we let and consider and then only and are obtained at the same time. Then for the proof of the sufficiency, by substituting the above conditions into , one can get the conclusion easily. Furthermore we can transform the above conditions into the following form: , . Thus the proof is completed.
From Theorem 5, one calculates easily the Jacobian determinant with respect to the functions , and variables , , Considering the conditions (22) and when holds, we obtain . At the same time, one can get that holds.
Thus we take some appropriate perturbations for the coefficients of system (6) and get the perturbed system with the following form: where . According to Theorem 2.2 in , system (6) has 3 limit cycles in the neighborhood of the origin; that is, system (5) has 3 limit cycles on the center manifold necessarily. One can refer to [9, 17, 18] for more details about constructing of limit cycles. And using numerical method, we get an example of one stable limit cycle for system (24) as shown in Figure 1(b).
Therefore, from the above discussion and Theorem 7, we have the following.
Theorem 10. For system (5), 3 and at most 3 small amplitude limit cycles can be bifurcated from the origin on the center manifold.
In summary, based on precise symbolic computation, we have investigated deeply multiple limit cycles bifurcation in the vicinity of the degenerate equilibrium for a three-dimensional system. Firstly by computing its quasi-Lyapunov constants, we solved its center-focus problem on center manifold; thus the center conditions are found and as weak foci the highest order is proved to be the fourth. Then we obtained that the system can generate 3 small limit cycles from nilpotent critical point and at most 3 small limit cycles from the origin via local bifurcation. For the flow of a high-dimensional system restricted to the center manifold, only one interesting example of multiple limit cycles bifurcating from its nilpotent singularity was given here, but we believe that more outcomes on this problem will be shown in the near future.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This work was supported by Scientific Research Foundation of Guangxi Education Department (no. ZD2014131), Research Foundation of Hezhou University (no. HZUBS201302), and Natural Science Foundation of China grants nos. 11261013 and 11371373.
- C. Sparrow, The Lorenz Equations: Bifcations, Chaos, and Strange Attractors, Springer, New York, NY, USA, 1982.
- T. Li, G. Chen, and Y. Tang, “On stability and bifurcation of Chen's system,” Chaos, Solitons and Fractals, vol. 19, no. 5, pp. 1269–1282, 2004.
- L. F. Mello and S. F. Coelho, “Degenerate Hopf bifurcations in the Lü system,” Physics Letters A, vol. 373, no. 12-13, pp. 1116–1120, 2009.
- M. Messias, D. D. C. Braga, and L. F. Mello, “Degenerate Hopf bifurcations in Chua's system,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 19, no. 2, pp. 497–515, 2009.
- J. Hofbauer and J. W.-H. So, “Multiple limit cycles for three-dimensional Lotka-Volterra equations,” Applied Mathematics Letters, vol. 7, no. 6, pp. 65–70, 1994.
- Z. Lu and Y. Luo, “Two limit cycles in three-dimensional Lotka-Volterra systems,” Computers & Mathematics with Applications, vol. 44, no. 1-2, pp. 51–66, 2002.
- M. Gyllenberg and M. Yan, “Four limit cycles for a three-dimensional competitive Lotka–Volterra system with a heteroclinic cycle,” Computers & Mathematics with Applications, vol. 58, no. 4, pp. 649–669, 2009.
- Q. Wang, W. Huang, and B. Li, “Limit cycles and singular point quantities for a 3D Lotka-Volterra system,” Applied Mathematics and Computation, vol. 217, no. 21, pp. 8856–8859, 2011.
- Q. Wang, Y. Liu, and H. Chen, “Hopf bifurcation for a class of three-dimensional nonlinear dynamic systems,” Bulletin des Sciences Mathématiques, vol. 134, no. 7, pp. 786–798, 2010.
- J. Carr, Applications of Centre Manifold Theory, vol. 35 of Applied Mathematical Sciences, Springer, New York, NY, USA, 1981.
- M. J. Alvarez and A. Gasull, “Generating limit cycles from a nilpotent critical point via normal forms,” Journal of Mathematical Analysis and Applications, vol. 318, no. 1, pp. 271–287, 2006.
- M. J. Álvarez and A. Gasull, “Monodromy and stability for nilpotent critical points,” International Journal of Bifurcation and Chaos, vol. 15, no. 4, pp. 1253–1265, 2005.
- J. Chavarriga, H. Giacomin, J. Gine, and J. Llibre, “Local analytic integrability for nilpotent centers,” Ergodic Theory and Dynamical Systems, vol. 23, no. 2, pp. 417–428, 2003.
- R. Moussu, “Symetrie et forme normale des centres et foyers degeneres,” Ergodic Theory and Dynamical Systems, vol. 2, no. 2, pp. 241–251, 1982.
- A. F. Andreev, A. P. Sadovskii, and V. A. Tsikalyuk, “The center-focus problem for a system with homogeneous nonlinearities in the case of zero eigenvalues of the linear pact,” Differential Equations, vol. 39, no. 2, pp. 155–164, 2003.
- H. Giacomini, J. Giné, and J. Llibre, “The problem of distinguishing between a center and a focus for nilpotent and degenerate analytic systems,” Journal of Differential Equations, vol. 227, no. 2, pp. 406–426, 2006.
- Y. Liu and J. Li, “Bifurcations of limit cycles and center problem for a class of cubic nilpotent system,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 20, no. 8, pp. 2579–2584, 2010.
- Y. Liu and J. Li, “On third-order nilpotent critical points: integral factor method,” International Journal of Bifurcation and Chaos, vol. 21, no. 5, pp. 1293–1309, 2011.
- F. Li, Y. Liu, and Y. Wu, “Center conditions and bifurcation of limit cycles at three-order nilpotent critical point in a seventh degree Lyapunov system,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 6, pp. 2598–2608, 2011.
- A. Buică, I. A. García, and S. Maza, “Existence of inverse Jacobi multipliers around Hopf points in : emphasis on the center problem,” Journal of Differential Equations, vol. 252, no. 12, pp. 6324–6336, 2012.
- V. F. Edneral, A. Mahdi, V. G. Romanovski, and D. S. Shafer, “The center problem on a center manifold in ,” Nonlinear Analysis, vol. 75, no. 4, pp. 2614–2622, 2012.
- F. S. Dias and L. F. Mello, “Analysis of a quadratic system obtained from a scalar third order differential equation,” Electronic Journal of Differential Equations, vol. 2010, 6 pages, 2010.
- J. Llibre, “On the integrability of the differential systems in dimension two and of the polynomial differential systems in arbitrary dimension,” The Journal of Applied Analysis and Computation, vol. 1, no. 1, pp. 33–52, 2011.
Copyright © 2014 Shugang Song 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.