# Bifurcation of Limit Cycles and Center Conditions for Two Families of Kukles-Like Systems with Nilpotent Singularities

Authors: Peiluan Li; Yusen Wu; Xiaoquan Ding
Journal: Journal of Function Spaces and Applications
 (2013)

---

## Abstract

We solve theoretically the center problem and the cyclicity of the Hopf bifurcation for two families of Kukles-like systems with their origins being nilpotent and monodromic isolated singular points.


---

## Body

## 1. Introduction

One of the main open problems in the qualitative theory of planar analytic differential systems is the case when the singular point is monodromic (the orbits move around the singular point). A monodromic point of an analytic system is either a center (i.e., a singular point with a punctures neighborhood filled with periodic orbits) or a focus (i.e., a singular point with a neighborhood where all the orbits are spirals which arrive at the equilibrium point in forward or backward time). The problem of distinguishing when a monodromic singular point is either a center or a focus is called thecenter problem.The so-called center problem for planar vector fields has been intensively and extensively studied over the last century and is also closely related to Hilbert’s 16th problem. The singular point can be elementary or not in the sense whether the corresponding Jacobian matrix has zero eigenvalues. If the eigenvalues of the quoted matrix are imaginary with real part null, the singular point may be a focus or a center, which is known as the celebratedPoincaré-Lyapunov center problem and has been theoretically solved by Poincaré [1] and Liapunov [2]. If the matrix of the linear part at the singular point has its two eigenvalues equal to zero, but it is not identically null, by Andreev [3] we know what is the behavior of the solutions in a neighborhood of the singular point, except if it is a center of a focus (nilpotent center problem).Yet for the bifurcation of limit cycles and center problem of nilpotent singular points in a planar vector field, its intrinsic dynamics is still far away from understanding due to the complexity and technical difficulties. Therefore, it is natural to restrict our study to nilpotent singularities. An analytic system of differential equations in the plane having an isolated nilpotent singularity, in some suitable coordinates, can be written as(1)x˙=y+X(x,y),y˙=Y(x,y),
with X(x,y) and Y(x,y) real analytic functions without constant nor linear terms defined in a certain neighborhood of the singularity.The study of nilpotent singularities and their unfolding for vector fields is important not only for mathematical interest but also for practical reasons. To solve the finite cyclicity problem in the second part of Hilbert’s 16th problem, Dumortier et al. [4, 5] presented an impressive list of 121 graphics that occur in quadratic systems, among which many are nilpotent cases, studied in [6–11]. On the other hand, the center problem for nilpotent singularities is another interesting issue and has been partially solved; see [3, 12–19], and so forth.Liu and Li [20] originated a method in studying the center problem of third-order nilpotent singularities by integral factor method, where a new definition of the focal values, quasi-Lyapunov constants, is provided. At the same time, the equivalence of quasi-Lyapunov constants and focal values is proved. A linear recursive formula for computing quasi-Lyapunov constants is also derived. Afterward, they proved that if the third-order nilpotent origin is a m-order weakened focus, then, by a small perturbation for the unperturbed system, there exist m limit cycles in a neighborhood of the origin. Meanwhile, the origin becomes an elementary critical point and two complex singular points.In the present paper, we are interested in the study of two families of Kukles-like systems with third-order nilpotent part of the form(y,-2x3) and, respectively, nonlinear part given by quartic and quintic homogeneous polynomial with an isolated singular point at the origin; that is,
(2)dxdt=y,dydt=-2x3+b40x4+b31x3y+b22x2y2+b13xy3+b04y4,(3)dxdt=y,dydt=-2x3+b50x5+b41x4y+b32x3y2+b23x2y3+b14xy4+b05y5.This paper will be organized as follows. In Section2 we summarize some general definitions and results about third-order nilpotent singularities. In Section 3 we present the conditions for third-order nilpotent singularities to be a center and give the proof of the cyclicity for system (2). In the last section, we present the parallel results for system (3).


## 2. The Formula to Compute the Quasi-Lyapunov Constants

The origin is a third-order monodromic critical point if and only if the system is of the following form:(4)dxdt=y+μx2+∑i+2j=3∞‍aijxiyj=X(x,y),dydt=-2x3+2μxy+∑i+2j=4∞‍bijxiyj=Y(x,y).The quasi-Lyapunov constant is defined and a way of computing them is provided in [20].Theorem 1.
For any positive integers and a given number sequence
(5){c0β},β≥3,
one can construct successively the terms with the coefficients cαβ satisfying α≠0 of the formal series
(6)M(x,y)=y2+∑α+β=3∞‍cαβxαyβ=∑k=2∞‍Mk(x,y),
such that
(7)(∂X∂x+∂Y∂y)M-(s+1)(∂M∂xX+∂M∂yY)=∑m=3∞‍ωm(s,μ)xm,
where for all k, Mk(x,y) is a k-homogeneous polynomial of x,y and sμ=0.It is easy to see that (7) is linear with respect to the function M, so that we can easily find in [20] the following recursive formulae for the calculation of cαβ and ωm(s,μ).Theorem 2.
Forα≥1, α+β≥3 in (6) and (7), cαβ can be uniquely determined by the recursive formula
(8)cαβ=1(s+1)α(Aα-1,β+1+Bα-1,β+1).
For m≥1, ωm(s,μ) can be uniquely determined by the recursive formula
(9)ωm(s,μ)=Am,0+Bm,0,
where
(10)Aαβ=∑k+j=2α+β-1‍[k-(s+1)(α-k+1)]×akjcα-k+1,β-j,Bαβ=∑k+j=2α+β-1‍[j-(s+1)(β-j+1)]×bkjcα-k,β-j+1.
Notice that in (10), one sets
(11)c00=c10=c01=0,c20=c11=0,c02=1,cαβ=0,ifα<0orβ<0.We see from Theorem1 that if the origin of system (4) is s-class or ∞-class, then, by choosing {cαβ}, such that
(12)ω2k+1(s,μ)=0,k=1,2,…,
we can obtain a solution group of {cαβ} of (12); thus, we have
(13)λm=ω2m+4(s,μ)2m-4s-1.Clearly, the recursive formulas presented by Theorem2 are linear with respect to all cαβ. Therefore, it is convenient to realize the computations of quasi-Lyapunov constants by using computer algebraic system like MATHEMATICA.


## 3. Center Conditions and Limit Cycles of System (2)

By Theorems1 and 2, we directly have the following.Lemma 3.
Assume thats is a natural number. One can derive a power series (6) for system (2) under which (7) is satisfied, where
(14)c[0,0]=0,c[1,0]=0,c[0,1]=0,c[2,0]=0,c[1,1]=0,c[0,2]=1.
In addition, for any natural numbers α,β, c[α,β] is given by the following recursive formula:
(15)c[α,β]=(-b40(1+s)(2+β)c[-5+α,2+β]+b31(1-(1+s)(1+β))c[-4+α,1+β]+2(1+s)(2+β)c[-4+α,2+β]+b22(2-(1+s)β)c[-3+α,β]+b13(3-(1+s)(-1+β))×c[-2+α,-1+β]+b04(4-(1+s)(-2+β))×c[-1+α,-2+β]-b40(1+s)(2+β)c[-5+α,2+β])×(α(s+1))-1,
and, for any natural number m, ωm is given by the following recursive formula:
(16)ωm=-b40(1+s)c[-4+m,1]+b31c[-3+m,0]+2(1+s)c[-3+m,1]+b22(3+s)c[-2+m,-1]+b13(3+2(1+s))c[-1+m,-2]+b04(4+3(1+s))c[m,-3].The straightforward computing according to Lemma3 with MATHEMATICA gives
(17)ω6=0,ω7=12(b31-2b31s+6c03+6sc03),ω8~-110b40b31(-3+4s),ω9~-2b13(-1+s).To makeω9=0, we will consider two different cases: b13=0 and s=1.Case 1 (b13=0).
Taking expression (13) into account, the first four quasi-Lyapunov constants of system (2) are given by
(18)λ1=0,λ2~110b40b31,λ3~121b31b22,λ4~-29b31b04.Now we can prove the following result.Theorem 4.
In this case, the origin of system (2) is a center if and only if the first four quasi-Lyapunov constants vanish; that is, one of the following two series is satisfied:
(19)b31=b13=0,(20)b40=b22=b13=b04=0.Proof.
It is direct to obtain conditions (19) and (20) by imposing λi=0,  i=1,2,3,4.
On the other hand, when condition (19) holds, system (2) goes over to
(21)dxdt=y,dydt=-2x3+b40x4+b22x2y2+b04y4,
whose vector field is symmetric with respect to the x-axis. When condition (20) holds, system (2) goes over to
(22)dxdt=y,dydt=x3(-2+b31y),
whose vector field is symmetric with respect to the y-axis. This ends the proof.Case 2 (s=1).
Taking expression (13) into account, the first five quasi-Lyapunov constants of system (2) are given by
(23)λ1=0,λ2~110b40b31,λ3~1105(5b31b22+9b40b13),λ4~16300(680b22b13-1400b31b04-117b403b13),λ5~18085000b13×(476000b04-286495b402b22+18018b405).For the same reason as before, we obtain the following.Theorem 5.
In this case, the origin of system (2) is a center if and only if the first five quasi-Lyapunov constants vanish; that is, the following series is satisfied:
(24)b40=b22=b04=0.Proof.
It is direct to obtain conditions (24) by imposing λi=0,  i=1,2,3,4,5.
On the other hand, when condition (24) holds, system (2) goes over to
(25)dxdt=y,dydt=x(-2x2+b31x2y+b13y3),
whose vector field is symmetric with respect to the y-axis. Then, we have proved the theorem.It is not hard to see that condition (20) is a particular case of condition (24); in summary, the following result characterizes the center of system (2).Theorem 6.
System (2) has a center at the origin if and only if one of the following two sets of conditions holds:
(
I
)
b
31
=
b
13
=
0,
(
I
I
)
b
40
=
b
22
=
b
04
=
0.Remark 7.
Indeed, in the light of the theory of [21], the origin is not only a center but also an analytic center under condition (I).Referring to the expressions given in (18) and (23), we have the following.Lemma 8.
The origin of system (2) is a fine focus of maximum order five. Further, it is of order five if and only if
(26)b40=b31=b22=0,b13b04≠0.Consider the perturbed system of (2)
(27)dxdt=δ(ε)x+y,dydt=2δ(ε)y-2x3+b21(ε)x2y+b40(ε)x4+b31(ε)x3y+b22(ε)x2y2+b13xy3+b04y4.Theorem 9.
Under a small perturbation of system (2), for a small parameter δ, in a neighborhood of the origin of system (27), there exist exactly 5 small amplitude limit cycles enclosing the elementary node O(0,0).The next example shows this fact.Example 10.
Take(28)δ(ε)=-ε30,b21(ε)=ε20,b40(ε)=ε6,b31(ε)=-ε6,b22(ε)=-ε2,b13=b04=1.
Straightforward computations by using expression (23) give the first 5 quasi-Lyapunov constants of the origin of system (27):
(29)λ1=13ε20+o(ε20)≈0.333333ε20+o(ε20),λ2~-110ε12+o(ε12)=-0.1ε12+o(ε12),λ3~335ε6+o(ε6)≈0.0857143ε6+o(ε6),λ4~-34315ε2+o(ε2)≈-0.107937ε2+o(ε2),λ5~681155+o(1)≈0.0588745+o(1).
Thereby, for 0<ε≪1, system (27) has 5 limit cycles Γk:r=r~(θ,hk(ε)) in a small neighborhood of the origin, where hk(ε)=o(εk),  k=1,2,3,4,5.


## 4. Center Conditions and Limit Cycles of System (3)

Similar to Lemma3, we have the following.Lemma 11.
Assume thats is a natural number. One can derive a power series (6) for system (3) under which (7) is satisfied, where
(30)c[0,0]=0,c[1,0]=0,c[0,1]=0,c[2,0]=0,c[1,1]=0,c[0,2]=1.
In addition, for any natural numbers α,β, c[α,β] is given by the following recursive formula:
(31)c[α,β]=(-b50(1+s)(2+β)c[-6+α,2+β]+b41(1-(1+s)(1+β))c[-5+α,1+β]+b32(2-(1+s)β)c[-4+α,β]+2(1+s)(2+β)c[-4+α,2+β]+b23(3-(1+s)(-1+β))×c[-3+α,-1+β]+b14(4-(1+s)(-2+β))×c[-2+α,-2+β]+b05(5-(1+s)(-3+β))×c[-1+α,-3+β]-b50(1+s)(2+β)c[-6+α,2+β])×(α(s+1))-1,
and, for any natural number m, ωm is given by the following recursive formula:
(32)ωm=-b50(1+s)c[-5+m,1]+b41c[-4+m,0]+b32(3+s)c[-3+m,-1]+2(1+s)c[-3+m,1]+b23(3+2(1+s))c[-2+m,-2]+b14(4+3(1+s))c[-1+m,-3]+b05(5+4(1+s))c[m,-4].Taking expression (13) into account, the first four quasi-Lyapunov constants of system (3) are given by
(33)λ1=0,λ2~15b41,λ3~27b23,λ4~43b05.Hence, we arrive at the following.Theorem 12.
System (3) has a center at the origin if and only if the following set of conditions holds:
(
I
I
I
)
b
41
=
b
23
=
b
05
=
0.Proof.
It is direct to obtain condition (III) by imposingλi=0,  i=1,2,3,4.
On the other hand, when condition (III) is satisfied, system (3) reduces to
(34)dxdt=y,dydt=x(-2x2+b50x4+b32x2y2+b14y4),
whose vector field is symmetric with respect to the origin. The claim follows.Remark 13.
Indeed, in the light of the theory of [21], the origin is not only a center but also an analytic center under condition (III).We obtain the following by expression (33).Lemma 14.
The origin of system (3) is a fine focus of maximum order four. Further, it is of order four if and only if
(35)b41=b23=0,b05≠0.Consider the perturbed system of (3):
(36)dxdt=δ(ε)x+y,dydt=2δ(ε)y-2x3+b21(ε)x2y+b50x5+b41(ε)x4y+b32x3y2+b23(ε)x2y3+b14xy4+b05y5.Theorem 15.
Under a small perturbation of system (3), for a small parameter δ, in a neighborhood of the origin of system (36), there exist exactly 4 small amplitude limit cycles enclosing the elementary node O(0,0).Example 16.
Take(37)δ(ε)=ε20,b21(ε)=-ε12,b41(ε)=ε6,b23(ε)=-ε2,b05=1.
Straightforward computations by using expression (33) give the first 4 quasi-Lyapunov constants of the origin of system (36):
(38)λ1=-13ε12+o(ε12)≈-0.333333ε12+o(ε12),λ2~15ε6+o(ε6)=0.2ε6+o(ε6),λ3~-27ε2+o(ε2)≈-0.285714ε2+o(ε2),λ4~43+o(1)≈1.33333+o(1).
Thereby, for 0<ε≪1, system (36) has 4 limit cycles Γk:r=r~(θ,hk(ε)) in a small neighborhood of the origin, where hk(ε)=o(εk),  k=1,2,3,4.



---

*Source: 403836-2013-11-06.xml*