New family of cubic Hamiltonian centers

We characterize the 11 non-topological equivalent classes of phase portraits in the Poincaré disc of the new family of cubic polynomial Hamiltonian differential systems with a center at the origin and Hamiltonian H=12((x+ax2+bxy+cy2)2+y2),\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} H= \frac{1}{2} ( (x + a x^2 + b x y + c y^2)^2+y^2 ), \end{aligned}$$\end{document}with a2+b2+c2≠0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^2+b^2+c^2\ne 0$$\end{document}.


Introduction
For a given family of real planar polynomial differential systems depending on parameters one of the main problems is the characterization of their centers and their phase portraits. The notion of center goes back to Poincaré in [16]. He defined it for differ-B Martín-Eduardo Frías-Armenta martineduardofrias@gmail.com; eduardo@gauss.mat.uson.mx Jaume Llibre jllibre@mat.uab.cat ential systems on the real plane; i.e., a singular point surrounded by a neighborhood fulfilled of closed orbits with the unique exception of the singular point.
The classification of the centers of the real polynomial differential systems started with the quadratic ones with the works of Kapteyn [10,11], Bautin [3], Vulpe [21], Schlomiuk [18,19],Żoładek in [23], etc. Schlomiuk et al. in [20] described a brief history of the problem of the center in general and it includes a list of 30 papers covering the history of the center for the quadratic polynomial differential systems (see pages 3, 4 and 13).
While the centers and all their phase portraits have been characterized for all the quadratic polynomial differential systems, this is not the case for the polynomial differential systems of degree larger than 2, but for such systems there are many partial results. Thus, the centers for cubic polynomial differential systems of the form linear with homogeneous nonlinearities of degree 3 were classified by Malkin [12], and Vulpe and Sibirskii [22], and their phase portraits when they are Hamiltonian have been classified by Colak et al. in [6,7]. Moreover, for polynomial differential systems which are linear with homogeneous nonlinearities of degree k > 3, the centers are not classified, but there are partial results for k = 4, 5; see Chavarriga and Giné [4,5], respectively.
Nowadays, we are very far from obtaining a complete classification of the centers for the class of all polynomial differential systems of degree 3. In any case, some interesting results on some subclasses of cubic polynomial differential systems are the ones of Rousseau and Schlomiuk [17] and the ones ofŻoładek [24,25].
Note that such cubic polynomial differential systems have terms of degree 1, 2 and 3 and consequently are not studied in the previously mentioned papers dedicated to the centers of cubic polynomial differential systems. It is clear that the origin (0, 0) is an isolated minimum of the Hamiltonian function H given in (1), and so the origin is a center of the Hamiltonian system (2) because near (0, 0) the curves H (x, y) = constant are closed.
In what follows, we shall talk interchangeably of the cubic polynomial Hamiltonian differential system (2), or of its associated cubic polynomial Hamiltonian vector field X H = (−y − (bx + 2cy)(x + ax 2 + bx y + cy 2 )) ∂ ∂ x + ((1 + 2ax + by)(x + ax 2 + bx y + cy 2 )) ∂ ∂ y . Let X be a polynomial vector field. See information on the notations and definitions used in this paper on the Poincaré disc and on the compactified vector field p(X ) associated with X in Sect. 2.1.
We say that two polynomial vector fields X and Y on R 2 are topologically equivalent if there is a homeomorphism on the Poincaré disc preserving the infinity and carrying orbits of the vector field p(X ) into orbits of the vector field p(Y), preserving or reversing simultaneously the sense of all orbits.
Our main result can be stated as follows.
Theorem 1 is proved in Sect. 5. In Sects. 3 and 4, we study the finite and infinite singular points of the system (2).
We must mention that our study has only used the program P4 for doing the pictures of the phase portraits and that we studied the local phase portraits of finite and infinite singular points in the Poincaré disc and the separatrices of the system analytically. We also must say that we prefer to provide the phase portraits in the Poincaré disc, which we have studied analytically, drawn with the program P4, because then we have the analytical qualitative study of these phase portraits and also its quantitative study for This phase portrait corresponds to the values a = 1, b = 1 and c = 3. This phase portrait has three canonical regions the considered values of the parameters. Usually, this is not possible because some separatrices become too close and we cannot distinguish between them, but for the Hamiltonian systems studied here, only two phase portraits (see Figs. 3 and 10) have one center which becomes too close to the separatrices, but since we know analytically that they are centers, these two phase portraits can be well understood.

Notations and basic results
In this section, we recall the basic definitions, notations and results that we will need for the analysis of the local phase portraits of the finite and infinite singular points of the polynomial Hamiltonian systems (2), and for doing their phase portraits in the Poincaré disc.
We denote by P 3 (R 2 ) the set of polynomial vector fields in R 2 of the form X (x, y) = (P(x, y), Q(x, y)), where P and Q are real polynomials in the variables x and y, such that the maximal degree of P and Q is 3.

Poincaré compactification
In this subsection, we provide a brief summary of the Poincaré compactification of a polynomial vector field X of degree 3, presenting the formulas that we need for This phase portrait corresponds to the values a = −1/2, b = 8 and c = −7. It has two centers, one saddle and one cusp. One of the centers is difficult to see, but it is clear where it is. We do not draw a periodic orbit in each canonical region having as inner boundary a center, because there is not sufficient space. This phase portrait has four canonical regions studying the phase portraits of system (2). For the proof of these formulas and all the details on the Poincaré compactification, see for instance Chapter 5 of [8].
The Poincaré disc is the image of the northern hemisphere of the Poincaré sphere S 2 = {y = (y 1 , y 2 , y 3 ) ∈ R 3 : y 2 1 + y 2 2 + y 2 3 = 1} into the plane y 3 = 0 under the projection (y 1 , y 2 , y 3 ) → (y 1 , y 2 ). The interior of the Poincaré disc is diffeomorphic to the plane R 2 , where we have a copy X of our polynomial vector field X through the mentioned diffeomorphism, and the boundary of this disc S 1 corresponds to the infinity of R 2 . There is a unique analytic extension of the vector field X to the Poincaré disc called the Poincaré compactification p(X ) of X . The infinity S 1 is invariant under the flow of p(X ).
On the sphere S 2 , we take for i = 1, 2, 3 the six local charts U i = {y ∈ S 2 : y i > 0} and V i = {y ∈ S 2 : y i < 0}. Let F i (y) = −G i (y) = (y j /y i , y k /y i ) with j < k and j, k = i. We denote by (u, v) the value of F i (y) or G i (y) for any i = 1, 2, 3 (thus, (u, v) means different coordinates in the distinct local charts). Easy computations give for p(X ) the next expressions: The expression for V i is the same as that for U i except for a change of sign. In the local charts with subindices i = 1, 2, v = 0 always denotes the points of S 1 , i.e., the points of the infinity. In what follows, we shall omit the factor (z) by a convenient scaling of the vector field p(X ). Thus, we have a polynomial vector field for p(X ) in every local chart.
The singular points of p(X ) which are in the interior of the Poincaré disc are called the finite singular points, which correspond with the singular points of X , and the singular points of p(X ) which are in S 1 are called the infinite singular points of X .
We note that studying the infinite singular points of the local chart U 1 , we obtain also the ones of the local chart V 1 . It only remains to be seen if the origin of the local chart U 2 and consequently the origin of the local chart V 2 are infinite singular points.
This phase portrait corresponds to the values a = −3/5, b = 5 and c = −8, and it has three centers and two saddles as finite singular points. One of the centers is difficult to see in the figure, but it is clear where it is. We do not draw a periodic orbit in the canonical region having as inner boundary the left center, because there is not sufficient space. This phase portrait has five canonical regions

Singular points
Let X H be the vector field given in (3) with a 2 + b 2 + c 2 = 0, and let p(X H ) be its Poincaré compactification. Let q be a singular point of p(X H ).
If some of the two eigenvalues λ 1 and λ 2 of the linear part of the vector field, p(X H ) at the singular point q is not zero, then this singular point is elementary.
The local phase portrait of a finite elementary singular point q with both eigenvalues non-zero for a Hamiltonian system is well known, it is a saddle if λ 1 λ 2 = 0 (which in fact can only be λ 1 λ 2 < 0) and a center if Reλ 1 = Reλ 2 = 0; for more details, see for instance, [2]. We recall that the flow of a Hamiltonian system in the plane preserves the area, so their finite singular points cannot have parabolic and elliptic sectors, only hyperbolic sectors, or they must be centers.
The local phase portrait of an infinite elementary singular point q with both eigenvalues non-zero can be studied using, for instance, Theorem 2.15 of [8].
Let q be an elementary singular point with an eigenvalue zero. Then, q is called a semi-hyperbolic singular point. The local phase portrait of a semi-hyperbolic singular point of a Hamiltonian system only can be a saddle; see again for example [2]. The local phase portrait of an infinite semi-hyperbolic singular point q can be described using, for instance, Theorem 2.19 of [8].
When both eigenvalues are zero, but the linear part of p(X H ) at the singular point q is not the zero matrix, we say that q is a nilpotent singular point. The local phase portrait of a nilpotent singular point can be studied using Theorem 3.5 of [8]. If q is nilpotent and finite, then it only can be a saddle, a cusp, or a center.
Finally, if the Jacobian matrix of p(X H ) at the singular point q is identically zero and q is isolated inside the set of all singular points, then we say that q is a linearly zero singular point. The study of the local phase portraits of such singular points needs special changes of variables called blow-ups; see for more details Chapter 3 of [8], or [1].

Local phase portraits of the singular points
We want to classify in the Poincaré disc, modulo topological equivalence, the phase portraits of p(X H ) being X − H the polynomial Hamiltonian vector fields given in (3). For doing that, we must start by classifying the local phase portraits at all finite and infinite singular points of p(X H ) in the Poincaré disc. This will be made by using the techniques described in Sect. 2.2.

Phase portraits in the Poincaré disc
In this subsection, we shall see how to characterize the phase portraits of p(X H ) in the Poincaré disc for the polynomial Hamiltonian vector fields X H of (3).
A separatrix of p(X H ) is an orbit which is either a singular point, or a trajectory which lies in the boundary of a hyperbolic sector of a finite or infinite singular point, or any orbit contained in the infinity S 1 . Neumann [14] proved that the set formed by all separatrices of p( The open connected components of D 2 \ S( p(X H )) are called canonical regions of X H or of p(X H ). A separatrix configuration is the union of S( p(X H )) plus one solution chosen from each canonical region. Two separatrix configurations S( p(X H 1 )) and S( p(X H 2 )) are topologically equivalent if there is an orientation preserving or reversing homeomorphism, which maps the trajectories of S( p(X H 1 )) into the trajectories of S( p(X H 2 )). The following result is due to Markus [13], Neumann [14] and Peixoto [15], who proved it independently.

Theorem 2
The phase portraits in the Poincaré disc of two compactified polynomial differential systems p(X H 1 ) and p(X H 2 ) are topologically equivalent if and only if their separatrix configurations S( p(X H 1 )) and S( p(X H 2 )) are topologically equivalent. In summary, for drawing the complete phase portrait of a Hamiltonian vector field X H in the Poincaré disc, we must draw all the separatrices of p(X H ) plus one orbit in each canonical region.

The finite singular points
In the next lemma, we characterize the finite singular points of our cubic Hamiltonian system (2).

Lemma 3
The following statements hold. (2) has three finite singular points: and the unique real solution q 1 of the system Moreover, p 1 and p 2 are centers, and q 1 is a saddle.
Moreover, p 1 and p 2 are centers, q 1 is a saddle and q 2 is a cusp. (2) has five finite singular points: p 1 , p 2 and the three real solutions q 1 , q 2 and q 3 of system (6). Moreover, p 1 , p 2 and a q i are centers and the q j with j = i is a saddle.
with p 1 being a center and q 1 a saddle.
if they exist.
where p 1 and p 2 are centers and p 3 is a saddle.
Proof The singular points of system (2) are the solutions of the system The singular points p 1 and p 2 are obtained solving the system y = 0, x + ax 2 + bx y + cy 2 = 0.
For solving system (6), we first eliminate the variable x between the previous two equations and get while eliminating the variable y we obtain The discriminant of these two cubic equations are is positive, zero, or negative, then it follows that system (6) has a unique solution, two solutions, or three solutions, respectively. See for more details [9]. Adding to these solutions the two solutions p 1 and p 2 , we obtain the singular points which appear in statements (a), (b) and (c) of the lemma, respectively. The singular points of the remaining statements of the lemma follow easily by direct computations.
The local phase portraits of the singular points which appear in all the statements of the lemma have been studied easily using the results of Sect. 2.2.

The infinite singular points
Let q be an infinite singular point of some polynomial vector field and let h be a hyperbolic sector of q. We say that h is degenerated if its two separatrices are contained in infinity (i.e., in S 1 using the notation of Sect. 2.1), otherwise h is called nondegenerated.
In the following lemma, we characterize the infinite singular points of our cubic Hamiltonian system (2), using the notations of Sect. 2.1.

Lemma 4
The following statements hold. Proof From (4) the Hamiltonian system (2) in the local chart U 1 becomeṡ Therefore, if c = 0 and b 2 − 4ac > 0, this system has two infinite singular points in this chart, namely If c = 0 and b 2 − 4ac = 0, this system has only one infinite singular point in the chart U 1 . If c = 0 and b = 0, again the system has only one infinite singular point in U 1 . Finally if b 2 − 4ac < 0, then there are no infinite singular points in U 1 .
The charts U 1 and its symmetric V 1 cover all the infinity except the origins of the local charts U 2 and V 2 .
From (5) the Hamiltonian system (2) in the local chart U 2 iṡ So the origin of the chart U 2 is an infinite singular point if and only if c = 0. All the infinite singular points have linear part identically zero. So, to know its local phase portrait, we must do blowups; see Sect. 2.2.
From the different statements of the lemma, we see that there are only two different topological kinds of local phase portraits at the infinite singular points: (i) the local phase portrait on the Poincaré sphere at an infinite singular point is formed by two degenerate hyperbolic sectors; (ii) or the local phase portrait on the Poincaré sphere at an infinite singular point is formed by one hyperbolic and one elliptic sector separated by two parabolic ones, and these two parabolic sectors contain the infinity; (iii) or the local phase portrait on the Poincaré sphere at an infinite singular point is formed by two non-degenerate hyperbolic sectors separated by two parabolic sectors, and these two parabolic sectors contain the infinity.
We only prove using blowups one local phase portrait of type (i), and the remainder local phase portraits of this type which appear in the statements of the lemma are proved in a similar way.
For proving one local phase portrait of type (i), we have chosen the origin of the local chart U 2 when a = c = 0 and b = 0 (see Fig. 9); thus, we must study the local phase portrait at the origin of the systeṁ Since the linear part at the origin is identically zero, for studying its local phase portrait, we do a directional blowup in the direction of the v-axis, i.e., (u, v) → (u, w) with w = v/u, then system (9) in the new variables (u, w), after eliminating a common factor u between the components (u,ẇ), on performing a scaling of the independent variable, becomesu We note that this change of variables has blow up the origin of system (9) to the whole w-axis. The unique singular point of system (10) on the w-axis is the origin which is a saddle. Going back through the changes of coordinates, the scaling and the blowup, and taking into account that for system (9) we haveu| u=0 = −v 2 , the local phase portrait at the origin is either the one that we want to prove, i.e., the one described in (i), or it has some hyperbolic sector with both separatrices tangent to the v-axis. To eliminate this second possibility, we do another directional blowup, but now in the direction of the v-axis, i.e., (u, v) → (w, v) with w = u/v, then system (9) in the new variables (w, v), after eliminating a common factor w between the components (ẇ,v) doing a scaling of the independent variable, may be written aṡ This system has no singular points on the w-axis. This implies that when we go back through the changes of variables, the scaling and the blowup, the mentioned second possibility does not occur. This completes the proof of the local phase portrait (i).
The local phase portrait of type (ii) only appears in statement (g) when a = 0 and bc = 0, and at the origin of the local chart U 1 (see Fig. 7); thus, we must study the local phase portrait at the origin of the systeṁ Since the origin is a linearly zero singular point, we do the directional blowup (u, v) → (u, w) with w = v/u, then system (11) in the new variables (u, w), after eliminating a common factor u between the components (u,ẇ), on performing a scaling of the independent variable, becomeṡ We note that this change of variables has blown up the origin of system (11) to the whole w-axis. There are two singular points of system (12) on the w-axis, the origin which is a saddle and the point (−b, 0) which is linearly zero. For studying the local phase portrait of the singular point (−b, 0), first we translate it to the origin of coordinates doing the change (u, w) → (u, W = w + b). So, system (12) in the variables (u, W ) iṡ We do the directional blowup (u, W ) → (u, z) with z = W/u, then system (13) in the new variables (u, z), after eliminating a common factor u 2 between the components (u,Ẇ ), on performing a scaling of the independent variable, becomeṡ This system has no singular points on the z-axis, but it has a contact point on this axis at the point (0, −c); more precisely in that point, one orbit is tangent to the z-axis and does not cross it. Going back through the changes of coordinates, the scaling and the blowup, we get that the local phase portrait at the origin of system (13) has one elliptic sector and one hyperbolic sector, in opposite quadrants of the plane (u, W ) separated by two parabolic sectors. Consequently, this is the local phase portrait at the singular point (−b, 0) of system (12). Now, going back through the change of coordinates from system (12) to system (11), we obtain the local phase portrait (ii). The local phase portrait of type (iii) only appears in statement (e) when c = a = 0 and b = 0, and at the origin of the local chart U 1 (see Fig. 9); thus, we must study the local phase portrait at the origin of the systeṁ Since the origin is a linearly zero singular point, we do the directional blowup (u, v) → (u, w) with w = v/u; then the system (14) in the new variables (u, w), after eliminating a common factor u between the components (u,ẇ), performing a scaling of the independent variable, becomeṡ We note that this change of variables has blown up the origin of system (14) to the whole w-axis. There are two singular points of system (15) on the w-axis, the origin of which is a saddle and the point (−b, 0) which is linearly zero. For studying the local phase portrait of the singular point (−b, 0), first we translate it to the origin of coordinates by doing the change (u, w) → (u, W = w + b). So system (12) in the variables (u, W ) iṡ We do the directional blowup (u, W ) → (u, z) with z = W/u, then system (16) in the new variables (u, z), after eliminating a common factor u 2 between the components (u,Ẇ ), performing a scaling of the independent variable, becomeṡ The unique singular point of this system on the z-axis is the origin, which is a saddle. Now, we assume that b > 0, and the proof in the case b < 0 is analogous. This saddle has two unstable separatrices on the u-axis, and the two stable ones tangent to a straight line through the origin which has a negative slope. Going back through the changes of coordinates, the scaling and the blowup, we get that the local phase portrait at the origin of system (16), which in W ≥ 0 has un unstable parabolic sector with the orbits tangent to the W -axis with the exception of the two ones on the u-axis, and in W ≤ 0 it has two hyperbolic sectors with both separatrices tangent to the u-axis, separated by one parabolic sector. Therefore, this is the local phase portrait at the singular point (−b, 0) of system (15). Now, going back through the change of coordinates from system (15) to system (14), we obtain the local phase portrait (iii).

Global phase portraits
Taking into account the results on the finite and infinite singular points given in Sects. 3 and 4, respectively, and drawing the curves H (x, y) = h passing through the finite singular points, we shall obtain the different phase portraits of the Hamiltonian system (2) that we describe in what follows. We recall that the figures of the phase portraits in the Poincaré disc are drawn with the program P4 described in Chapters 9 and 10 of [8]. In this way as we have mention in Sect. 1, we have the analytical qualitative study of the phase portraits in the Poincaré disc that we have done and also its quantitative study for the considered values of the parameters. We separate the proof of Theorem 1 in nine cases. Case 1: Assume that b(b 2 − 4ac) = 0 and a(27ab 2 + 4(2a − c) 3 ) > 0. Then, by Lemma 3(a), system (2) has three finite singular points, two centers and one saddle.
We need to distinguish two subcases according to the infinite singular points. First, the subcase b 2 − 4ac > 0. If c = 0, then by Lemma 4(a) the system has two pairs of infinite singular points in the local charts U 1 ∪ V 1 , and each of these infinite singular points is formed by two degenerate hyperbolic sectors. If c = 0, then, by Lemma 4(d), the system has two pairs of infinite singular points: one pair in U 1 ∪ V 1 and the other at the origins of the charts U 2 and V 2 . Again, all the local phase portraits at these infinite singular points are formed by two degenerate hyperbolic sectors. So near infinity, the orbits rotate, they cannot spiral because they must be contained in some level curve H (x, y) = h, and H (x, y) is a polynomial; therefore near infinity, the orbits are closed, i.e., are periodic orbits.
Additionally, we have two centers with their periodic orbits surrounding them. In short, we have three continuous periodic orbits, one reaching infinity and the other two reaching the two centers. The common boundaries of these continuums are separatrices which only can belong to the separatrices of the unique saddle. This is checked drawing the level curve of the Hamiltonian H on the finite saddle. Thus, we get that the phase portrait of system (2) in this subcase is topologically equivalent to the one of Fig. 1.
Second, the subcase b 2 − 4ac < 0. By Lemma 4(c), the system has no infinite singular points. Using the same arguments of the previous case, we obtain that the phase portrait of the system now is topologically equivalent to the one of Fig. 2. Case 2: Assume that ab(b 2 − 4ac) = 0 and 27ab 2 + 4(2a − c) 3 = 0. Then, by Lemma 3(b), system (2) has four finite singular points, two centers, one saddle and one cusp. The singular points at infinity are the same as in Case 1. Therefore in this case, the phase portrait of the system is topologically equivalent to the one of Fig. 3 if b 2 − 4ac > 0, and to the one of Fig. 4 if b 2 − 4ac < 0. Case 3: Assume that b(b 2 − 4ac) = 0 and a(27ab 2 + 4(2a − c) 3 ) < 0. Then by Lemma 3(c), system (2) has five finite singular points, three centers and two saddles. Again, we need to distinguish two subcases according the infinite singular points.
So the behavior at the infinite singular points is as in the first subcase of Case 1, i.e., all the infinite singular points are formed by two degenerate hyperbolic sectors. So the orbits near infinity rotate, but since they cannot be spirals they are periodic orbits. So near infinity, we have a continuum of periodic orbits and also three more continuums of periodic orbits surrounding the three centers. We shall see that there is another continuum of periodic orbits surrounding two centers. We note that the values of the Hamiltonian H are different in both saddles, so separatrices of different saddles cannot connect and can only connect between the same saddle. We study the two level curves H (x, y) = h passing through the two saddles and they are qualitatively the ones which appear in Fig. 5. Removing all the separatrices of these two saddles, we obtain five canonical regions, four of them are filled with a continuum of periodic orbits, the continuum of the infinity and the three continuum of the centers; in the remaining canonical region, we know that the orbits near its two boundaries rotated, but since they cannot spiral, they are periodic orbits. By continuity, all the orbits of this canonical region are periodic (in fact, from Neumann's result [14] if a canonical region has a periodic orbit, all the orbits of the canonical region are periodic). Hence, the phase portrait of the system is topologically equivalent to the one of Fig. 5. Subcase 3.2: b 2 −4ac < 0 with no infinite singular points, then using the same kind of arguments as in the previous subcase the phase portrait of the system is topologically equivalent to the one of Fig. 6. Case 4: Assume that a = 0 and bc = 0. Then, by Lemma 3(d), system (2) has two finite singular points, one center and one saddle. By Lemma 4(f), the system has two pairs of infinite singular points in the local charts U 1 ∪V 1 , the local phase portrait of the singular points of one of these pairs is formed by two degenerate hyperbolic sectors, and the local phase portrait at the singular points of the other pair is formed by one non-degenerate hyperbolic sector and one elliptic sector separated by two parabolic sectors; these two parabolic sectors contain the infinity. So, first evaluating the level curve γ of the Hamiltonian H on the saddle, and after the level curve near the saddle but inside the region limited by γ encircling the center, we obtain that the phase portrait of the system is topologically equivalent to the one in Fig. 7. Case 5: Assume that a = 0, bc = 0 and b 2 + c 2 = 0. Then, by Lemma 3(e), system (2) has a unique finite singular point, a center at the origin of coordinates. If b = 0, by Lemma 4(b), the system has a unique pair of infinite singular points in U 1 ∪ V 1 and their local phase portraits are formed by two degenerate hyperbolic sectors. So the phase portrait of the system is topologically equivalent to the one of Fig. 8. If c = 0, by Lemma 4(e), we have two pairs of infinite singular points, one pair at the origins of the local charts U 1 and V 1 , and the local phase portrait of the singular points of this pair is formed by two non-degenerate hyperbolic sectors separated by two parabolic sectors; these two parabolic sectors contain the infinity. The other pair is at the origins of the local charts U 2 and V 2 , and the local phase portrait of the singular points of this pair is formed by two degenerate hyperbolic sectors. Consequently, the phase portrait of the system is topologically equivalent to the one of Fig. 9. Case 6: Assume that ac = 0, b = 0 and a(c − 2a) > 0. Then, by Lemma3(f), system (2) has five finite singular points, three centers and two saddles. We note that now the two saddles belong to the same level curve of the Hamiltonian H . Since a(c−2a) > 0, it follows that ac > 0. Then, by Lemma 4(c), there are no infinite singular points, and the phase portrait of the system is topologically equivalent to the one of Fig. 10. Case 7: Assume that ac = 0, b = 0 and a(c − 2a) ≤ 0. Then, by Lemma3(f), system (2) has three finite singular points, two centers and one saddle. If ac > 0, by Lemma 4(c), there are no infinite singular points, and the phase portrait is topologically equivalent to the one of Fig. 2). If ac < 0, then by Lemma 4(a), there is a pair of infinite singular points in U 1 ∪ V 1 , and the local phase portrait at each of these infinite singular points is formed by two degenerate hyperbolic sectors. Hence, the phase portrait of system (2) is topologically equivalent to the one in Fig. 1. Case 8: Assume that a = 0 and b = c = 0. Then, by Lemma 3(g), system (2) has three finite singular points, two centers and one saddle. By Lemma 4f, the system has only one pair of infinite singular points, the origins of the local charts U 2 and V 2 . The local phase portrait at these infinite singular points is formed by two degenerate hyperbolic sectors. Consequently, the phase portrait of the system is topologically equivalent to the one in Fig. 11. Case 9: Assume that ab = 0 and b 2 − 4ac = 0. Then, by Lemma 3(h), system (2) has three finite singular points, two centers and one saddle. By Lemma 4(c), the system has a unique pair of infinite singular points in U 1 ∪ V 1 , and the local phase portraits at these infinite singular points are formed by two degenerate hyperbolic sectors. So the phase portrait of this system is topologically equivalent to the one in Fig. 11. This completes the proof of Theorem 1.