The phase portrait of the Hamiltonian system associated to a Pinchuk map

In this paper we describe the global phase portrait of the Hamiltonian system associated to a Pinchuk map in the Poincaré disc. In particular, we prove that this phase portrait has 15 separatrices, five of them singular points, and 7 canonical regions, six of them of type strip and one annular.


INTRODUCTION
As far as we know, the simplest class of non-injective polynomial local diffeomorphisms of R 2 are the Pinchuk maps, constructed by Pinchuk (1994).The existence of these maps disproves the real Jacobian conjecture, that a polynomial local diffeomorphism of R 2 is globally injective.One open problem is to know what exactly fails in this conjecture.
One of the most known conditions for a local diffeomorphism to be a global one is that it is proper.The asymptotic variety of a map of R 2 is the set of points where the map is not proper (i.e., points that are limits of the map under sequences tending to infinity).In particular, a local diffeomorphism is a global diffeomorphism if and only if this set is empty.Gwoździewicz (2000) and Campbell (arXiv:math/9812032 in 1998, 2011) calculated the asymptotic variety of two Pinchuk maps in details.Our aim in this paper is to do a similar work, i.e., to describe a Pinchuk map, but now from a different point of view.
Let U ⊂ R 2 be an open connected set.Let F = (p, q) : U ⊂ R 2 → R 2 be a C 2 local diffeomorphism.Let H F (x, y) = (p(x, y) 2 + q(x, y) 2 )/2 and consider the Hamiltonian system ẋ = −(H F ) y (x, y), ẏ = (H F ) x (x, y), (1) Lemma 1.The singular points of system (1) are the zeros of F, each of them is a center of system (1).
The following is a generalization of the characterization of global invertibility of polynomial maps given by Sabatini (1998).This version is due to Braun andLlibre (arXiv:1706.02643 in 2017).
Theorem 2. Let z 0 ∈ U such that F(z 0 ) = (0, 0).The center z 0 of system (1) is global if and only if (i) F is globally injective and (ii) F(U) = R 2 or F(U) is an open disc centered at the origin.
In the special case that F is a polynomial map and U = R 2 , it follows that F(R 2 ) = R 2 provided F is injective (Białynicki-Birula and Rosenlicht 1962).Hence, in this case, z 0 is a global center of (1) if and only if F is globally injective.An application of this result was given by Braun et al. (2016).
Since the phase portrait on the Poincaré sphere of a Hamiltonian polynomial vector field having a global center is simple, i.e., at the infinite either it does not have singular points, or the infinite singular points are formed by two degenerate hyperbolic sectors (for Hamiltonian vector fields, the infinity contains only isolated singular points), it is interesting to know how complex can be the phase portrait of a non-global center of a Hamiltonian system (1).
In this paper we provide the qualitative global phase portrait of the Hamiltonian system (1) when F is given by the Pinchuk map considered by Campbell (1998Campbell ( , 2011)), after a translation in the target in order to have only a point z 0 such that F(z 0 ) = (0, 0).More precisely, we prove the following result.
Theorem 3. Let F = (p, q) : R 2 → R 2 , where (p, q + 208) : R 2 → R 2 is the Pinchuk map considered by Campbell (1998Campbell ( , 2011) ) (see the definition below).Then the phase portrait of the Hamiltonian system (1) in the Poincaré disc is topologically equivalent to the phase portrait given in Fig. 1.To prove Theorem 3, we first study the infinite singular points of system (1).These infinite singular points are very degenerate, and we apply homogeneous and quasi-homogeneous blow ups to study them.Then, we complete the proof of Theorem 3 by proving that the separatrix configuration of system ( 1) is qualitatively the one presented in Fig. 1.
We think that a good understanding of what fails in the real Jacobian conjecture could be interesting to investigate a related problem, the Jacobian conjecture in R 2 , that a polynomial local diffeomorphism whose Jacobian determinant is constant is globally injective.This conjecture remains unsolved until now.For the Jacobian conjecture we address the reader to the works of Bass et al. (1982) and Van Den Essen (2000).

INJECTIVITY, CENTERS AND A PINCHUK MAP
We begin with the proof of Lemma 1.
Proof of Lemma 1.Let z 0 be a singular point of the Hamiltonian system (1).We have , which is true if and only in F(z 0 ) = (p, q)(z 0 ) = (0, 0) because the Jacobian determinant of F is nowhere zero.
The point z 0 is a center of the Hamiltonian system (1) because it is an isolated minimum of H F .Now, we select the map F that we are going to work in this paper.Let t = xy − 1, h = t(xt + 1) and f = (xt + 1) 2 (t 2 + y).A Pinchuk map is a non-injective polynomial map with nowhere zero Jacobian determinant of the form (P, Q) : R 2 → R 2 such that P = h + f and Q = −t 2 − 6th(h + 1) − u(h, f ), where u is chosen so that det D(P, Q)(x, y) = t 2 + (t + f (13 + 15h)) 2 + f 2 .The following is the Pinchuk map studied by Campbell (1998Campbell ( , 2011)): According to Campbell (2011), the points (−1, −163/4) and (0, 0) ∈ R 2 have no inverse image under (p, q), all the other points of the curve which is a parametrization of the asymptotic variety of (p, q), have exactly one inverse image under this map, and the points of R 2 \ γ(R) have two inverse images.Hence, in particular, the point (0, 208) has precisely one inverse image under (p, q).We consider the map F = (p, q) : R 2 → R 2 given by the translation p(x, y) = p(x, y), q(x, y) = q(x, y) − 208. (2) Observe that F is a Pinchuk map according to our above-definition.Moreover, now there exists exactly one point z 0 ∈ R 2 such that F(z 0 ) = (0, 0).From Lemma 1 the point z 0 is the only finite singular point of system (1), corresponding to a non-global center of this system according to Theorem 2. Further, the curve s ∈ R, is the asymptotic variety of F, whose points have exactly one inverse image over F, but the points (−1, −995/4) and (0, 208), which have none.
From now on, we restrict our attention to the specific Pinchuk map (2).
We first calculate the coordinates of the point z 0 .Observe that xt + 1 = x 2 y − x + 1 is a factor of p.If this factor annihilates, then h = 0 and q = −t 2 − 208 < 0. The other factor of p is We observe that g(0, y) = y and q(0, y) = 50y − 799/4 do not annihilate at the same time, thus the first coordinate of the point z 0 is not 0.Moreover, since the leading coefficient of q(x, y) as a polynomial in y is −75x 15 , it follows that the first coordinate of z 0 will be a point where the resultant in y between g(x, y) and q(x, y) is zero.This resultant is the cubic c(x) = 31008391 − 11757152x − 155580672x 2 + 2239078400x 3 multiplied by −x 36 /64.The discriminant of c(x) is negative, so it has only one real root, which will be the first coordinate of the point z 0 .
Repeating a similar reasoning now looking g and q as polynomials in x, we calculate their resultant and obtain that its zero is the only real root of the cubic c(y) = 1789023641600 + 100675956992y + 26252413280y 2 + 1506138481y 3 , which will be the second coordinate of the point z 0 .
Hence z 0 = (−0, 22568337..., −17, 491214...) approximately.Since z 0 is a center, the only finite singular point of system (1), near z 0 the phase portrait of this system is simple.Indeed, since z 0 is the minimum point of H F , it follows that the gradient of H F points outward of each closed orbit of the center, and so each closed orbit of the center rotates in counterclockwise around z 0 .
In the following section we shall investigate the infinite of system (1).

THE INFINITE OF SYSTEM
In this section, we will use results and notations on the Poincaré compactification of polynomial vector fields of R 2 .In particular U i , V i , i = 1, 2, 3, are the canonical local charts of the Poincaré sphere S 2 .
For details on this technique we refer the reader to Chapter 5 of (Dumortier et al. 2006) or to (González Velasco 1969).
We call a singular point of a vector field linearly zero when the linear part of the vector field at this point is identically zero.
We begin by proving a general fact about the infinite singular points of Hamiltonian systems of the form (1). Writing , where H i is the homogeneous part of degree i of the polynomial H, it is simple to conclude that the infinite singular points (u, 0) of system (1) in the local charts U 1 and U 2 are An Acad Bras Cienc (2018) 90 (3) the points satisfying H d+1 (1, u) = 0 and H d+1 (u, 1) = 0, respectively.Let (u, 0) be an infinite singular point of system (1) and assume it is in the chart U 1 .The linear part of the vector field at (u, 0) is Therefore, (u, 0) is a linearly zero singular point.This proves the following result.
Lemma 4. The infinite singular points of the Hamiltonian system (1) are linearly zero.Now, we return to the Pinchuk map F defined by (2).Observe that the highest homogeneous part of H F (x, y) is 5625x 30 y 20 /2.Thus, the origins of the charts U 1 ,V 1 and U 2 ,V 2 are the infinite singular points of the Hamiltonian system (1), each of them linearly zero from Lemma 4.
We will use the quasi-homogeneous directional blow up technique to desingularize each of these infinite singular points.An exposition about blow-ups can be found in (Álvarez et al. 2011), see also Chapter 3 of (Dumortier et al. 2006).We now recall the directional blow up transformations.
By the quasi-homogeneous blow up in the positive (resp.negative) x-direction with weights α and β , or simply (α, β )-blow up in the positive (resp.negative) x-direction, we mean the transformation which carries the variables (x 1 , y 1 ) to the variables (x 2 , y 2 ) according to the formulas respectively.Similarly, by the quasi-homogeneous blow up in the positive (resp.negative) y-direction with weights α and β , or simply (α, β )-blow up in the positive (resp.negative) y-direction, we mean the transformations Clearly if α (resp.β ) is odd, then, the blow up in the positive x-direction (respec.y-direction) provides the information of the respectively negative blow ups.Also, if β is odd, the x-directional blow ups swap the second and third quadrants, while the y-directional blow ups swap the third and the fourth quadrants if α is odd.After the (α, β )-blow up in the x-direction, a system ẋ1 , in the positive and negative directions according to ±.Similarly, the (α, β )-blow up in the y-direction transforms ẋ1 2 ), in the positive and in the negative directions according to ±.
After the blow up in the x-direction (resp.y-direction) we cancel a common appearing factor x k 2 (y k 2 ) for a suitable k.So, if k is odd, the direction of the orbits are reversed in x 2 < 0 (y 2 < 0).
The weights α and β are chosen analyzing the Newton polygon of (P, Q), see the construction in (Álvarez et al. 2011).
The application of (α, β )-blow ups with αβ = 1 usually reduces the number of blow ups necessary for studying the local phase portrait of a linearly zero singular point.
To make the exposition clearer, we shall apply the most part of the blow ups in the x-direction.So, sometimes, we will first apply a xy-change, (x 1 , y 1 ) → (y 1 , x 1 ) = (x 2 , y 2 ), before making the blow-up.
In the next two subsections we will desingularize the origin of the charts U 1 and U 2 , respectively.We will denote the coordinates of the system in the step i of the algorithm as the variables (w i , z i ), so that after either a wz-change, a translation or a blow up, the new obtained system will be written in the variables (w i+1 , z i+1 ).In each step, we will denote the system ẇi Since the Hamiltonian system (1) with the polynomials p and q given by ( 2) has degree 49, it follows that for the calculations in each step of the algorithm we have to deal with polynomials of very high degree.So, we persuade these calculations with the algebraic manipulator Mathematica.We do not show in each step the whole expressions of the systems (P i , Q i ) because this would be impractical.

THE ORIGIN OF THE CHART U 1
We write the compactification of system (1) in the chart U 1 in the variables (w 0 , z 0 ), as (P 0 , Q 0 ).From Lemma 4, the singular point (0, 0) is linearly zero.
We first apply a wz-change and write the new system in the variables (w 1 , z 1 ) as (P 1 , Q 1 ).
The Newton polygon of system (P 1 , Q 1 ) has only one compact edge contained in the straight line x + 2y = 38.We apply (1, 2)-blow ups in the positive w-direction and in the positive and negative z-directions obtaining systems (P 2 , Q 2 ) and (P ± 2 , Q ± 2 ), in the variables (w 2 , z 2 ) and (w ± 2 , z ± 2 ), after canceling the common factors w 38 2 and (w ± 2 ) 38 , respectively.The first terms of these systems have the following expressions: The only singular point of (P 2 , Q 2 ) over the line w 2 = 0 is the linearly zero singular point (0, 1).The origin of the systems (P ± 2 , Q ± 2 ) are saddles as depicted in the planes w + 2 z + 2 and w − 2 z − 2 of Fig. 2. The reader can follow a schema of each step of the calculations in Fig. 2. We just need to analyze the origin of the systems (P ± 2 , Q ± 2 ), because the other singularities over the lines z ± 2 = 0 will correspond to the singularity (0, 1) of (P 2 , Q 2 ).Cienc (2018) 90 (3) Figura 2: We now analyze this linearly zero singularity.We first do a translation bringing this point to the origin, obtaining the new system (P 3 , Q 3 ) in the variables (w 3 , z 3 ).We also apply a wz-change obtaining the system (P 4 , Q 4 ) in the variables (w 4 , z 4 ).The Newton polygon of this system has two compact edges.We choose the one contained in the straight line x + y = 11.This compact edge has the point of negative abscissa (−1, 12), thus, concerning (1, 1)-blow ups, it follows from Proposition 3.2 of (Álvarez et al. 2011) that w 4 is not a characteristic direction, and so we only need to apply a w-directional (1, 1)-blow up, obtaining the system (P 5 , Q 5 ) in the variables (w 5 , z 5 ), after canceling the common factor w 11 5 .The first terms of (P 5 , Q 5 ) are:

An Acad Bras
Over the line w 5 = 0, the singular points of (P 5 , Q 5 ) are (0, 0) and two points of the form (0, z 5 ), with z 5 the two real solutions of The discriminant of this quartic equation is negative, thus, it has two real solutions.Those are approximately z 5 = −70722.424... and z 5 = −1.6611121....The singular point (0, 0) is linearly zero and the other two singular points are saddles, as represented in the w 5 z 5 -plane of Fig. 2. Now, we study the linearly zero point (0, 0) of (P 5 , Q 5 ).It is clear from (4) that the characteristic equation of (P 5 , Q 5 ) is identically zero, so (0, 0) is a dicritical singular point.We apply (1, 1)-blow ups in both the wand z-directions obtaining systems (P 6 , Q 6 ) and (P y 6 , Q y 6 ) in the variables (w 6 , z 6 ) and (w y 6 , z y 6 ), after canceling the factors w 9 6 and (z y 6 ) 9 , respectively.System (P y 6 , Q y 6 ) does not have (0, 0) as a singular point, so we just need to consider system (P 6 , Q 6 ) over the line w 6 = 0. We have P 6 (0, z 6 ) = 1 4 (z 6 + 1)(2z 6 + 1) 217160z  (Isaacson and Keller 1994); in the software Mathematica, the Sturm theorem is programed by the instruction CountRoots) we see that the polynomial of degree 6 multiplying (z 6 + 1)(2z 6 + 1)/4 in P 6 (0, z 6 ) has no real roots, so, the only singular points are (0, −1/2) and (0, −1).The first one is a weak focus and the second one is a saddle, as depicted in the plane w 6 z 6 of Fig. 2. Since the origin of (P 5 , Q 5 ) is dicritical, it follows that each orbit crossing the line w 6 = 0 will correspond to two orbits tending to (0, 0) in positive or negative directions.
We now begin the process of blowing down.
It is simple to conclude that the phase portrait of the system (P 5 , Q 5 ) close to the origin is qualitatively the one depicted in (a) of Fig. 3. Consequently, by considering also the information close to the other two singular points in the line w 5 = 0 (see the plane w 5 z 5 of Fig. 2), we can understand the behavior near the origin of system (P 4 , Q 4 ).We then apply a wz-change and conclude that the behavior of system (P 3 , Q 3 ) near the origin is the one presented qualitatively in (b) of Fig. 3.
By translating (0, 0) to (0, 1) and by using the information provided by the saddles of planes w ± 2 z ± 2 , we make the blow downs with α = 1 and β = 2, obtaining the origin of system (P 1 , Q 1 ).We then finally apply a wz-change and conclude that the origin of system (P 0 , Q 0 ) is qualitatively as drawn in Fig. 4.  As in the calculations made above, we write the compactified vector field in the chart U 2 as ( ẇ0 , ż0 ) = (P 0 , Q 0 ).The Newton polygon of (P 0 , Q 0 ) has two compact edges: one of them contained in the straight line 3x + 2y = 87.We apply a (3, 2)-blow up in the w-direction, obtaining the system ( ẇ1 , ż1 ) = (P 1 , Q 1 ) after canceling the factor w 87 1 .The first terms of P 1 and Q 1 are: The polynomials P 1 and Q 1 have degree 61.
It is clear that at (0, 0) we have a saddle.The other singular point of (P 1 , Q 1 ) in the line w 1 = 0 is (0, −1), and it is a linearly zero point.See the w 1 z 1 -plane of Fig. 5.The reader can follow the steps of the calculations in the schema shown in this figure.We just warn that, differently of Fig. 2, we already draw the final phase portrait of each step, including the behavior close to the linearly zero points (information that we will know only after persuading all the blow ups).We also apply (3, 2)-blow ups in the positive and negative z-directions, obtaining the systems ẇ±

wz-change
, respectively, with linearly zero singular points at (w ± 1 , z ± 1 ) = (0, 0).The polynomials P ± 1 and Q ± 1 have degree 30 and Q ± 1 is a factor of z ± 1 .We do not need to analyze the other singular points over the lines z ± 1 = 0, as the information provided by them is already contained in the w-directional blow up.We desingularize these points applying (1, 1)-blow ups in the w-direction.Here, we do not need to apply blow ups in the z-directions because the characteristic equations of the systems are and so w ± 1 = 0 are not characteristic directions.We obtain the systems (P ± 2 , Q ± 2 ) after canceling a factor (w ± 2 ) 5 .The polynomials P ± 2 and Q ± 2 have degree 45, and up to order 2 they have the same expressions: Thus, at (0, 0) the systems have a saddle, as depicted in the planes w ± 2 z ± 2 of Fig. 5.Moreover, any other singularity of the form (0, z ± 2 ) must satisfy By using Sturm's theorem, we conclude that this equation has no real solution.Now, we desingularize the point (0, −1) of system (P 1 , Q 1 ).First, we apply a translation to bring this point to the origin, obtaining the system (P 2 , Q 2 ) in the variables (w 2 , z 2 ).We also apply a wz-change obtaining the system (P 3 , Q 3 ) in the variables (w 3 , z 3 ).The Newton polygon of this system has only one compact edge contained in the line x + 2y = 11, and this edge has points of negative abscissa, so, concerning (1, 2)blow ups we just need to apply them in the w-direction, according to Proposition 3.2 of (Álvarez et al. 2011).Hence, we apply a (1, 2)-blow up in the positive w-direction, obtaining the system (P 4 , Q 4 ) in the variables (w 4 , z 4 ), after canceling a factor of w 11 4 .These polynomials have degree 90, and their first terms are: Clearly, (0, 0) is a singularity corresponding to a saddle.The other singular point in the line w 4 = 0 is (0, ξ ), where ξ is the only real root of the cubic which is approximately ξ = 18.8848.....This cubic has only one real root because its discriminant is negative.A calculation shows that (0, ξ ) is linearly zero.See the plane w 4 z 4 in Fig. 5. Now, we apply a translation to bring the point (0, ξ ) to the origin, obtaining system (P 5 , Q 5 ) written in the variables (w 5 , z 5 ).Since ξ is not a rational number, we do this translation with a parameter x, and As before the Newton polygon of this system has only one compact edge contained in the straight line x + y = 1.Moreover, the characteristic equation does not have z 9 = 0 as a solution.Here, we just apply a (1, 1)-blow up in the positive z-direction, obtaining system (P 10 , Q 10 ) in the variables (w 10 , z 10 ), after canceling the factor z 10 (here we do not use the superscript y as this is the only system in this step).The degree of this new system is 234, but as we are going to see, just the origin is a singular point in the line z 10 = 0.The first terms of P 10 and Q 10 are: with x = ξ .Now, over the line z 10 = 0, the singular points of (P 10 , Q 10 ) are (0, 0) and the points (w 10 , 0), with w 10 the real roots of with x = ξ .The discriminant of this equation after applying the polynomial remainder is It is simple to conclude that ∆(ξ ) < 0, thus, only (0, 0) is a singular point of (P 10 , Q 10 ) in z 10 = 0.This singular point is the saddle depicted in the plane w 10 z 10 of Fig. 5. Since the behavior near each appearing singular points in each step above is very simple, the blow down of each step is also very simple: following the arrays in Fig. 5, it is easy to conclude that the origin of U 2 has a degenerate hyperbolic sector as shown in the w 0 z 0 -plane of Fig. 5.

THE GLOBAL PHASE PORTRAIT
We begin with a background on separatrices and canonical regions of the Poincaré compactification p(X ) in the Poincaré disc D of a polynomial system ẋ = X (x).Let ϕ be the flow of p(X ) defined in D. As usual we denote by (U, ϕ) the flow of p(X ) on an invariant subset U ⊂ D. Two flows (U, ϕ) and (V, ψ) are said to be topologically equivalent if there exists a homeomorphism h : U → V sending orbits of (U, ϕ) onto orbits of (V, ψ) preserving or reversing the orientation of all the orbits.
Following Markus (1954), we say that the flow (U, ϕ) is parallel if it is topologically equivalent to one of the following flows: (i) the flow defined in R 2 by the system ẋ = 1, ẏ = 0; (ii) the flow defined in R 2 \ {(0, 0)} by the system in polar coordinates ṙ = 0, θ = 1; and (iii) the flow defined in R 2 \ {(0, 0)} by the system in polar coordinates ṙ = r, θ = 0. Parallel flows topologically equivalent to (i), (ii) and (iii) are called strip, annular and spiral (or radial), respectively.
We denote by γ x the orbit of p(X ) passing through x when t = 0 with maximal interval I x , and the positive (resp.negative) orbit of γ x by γ + x = {γ x (t) | t ∈ I x and t ≥ 0} (resp.γ − x = {γ x (t) | t ∈ I x and t ≤ 0}).Then we set a ± (x) = γ ± x \ γ ± x , here as usual γ ± x denotes the closure of γ ± x .Observe that a − (x) differs from α(x) in the case of periodic orbits and singular points: indeed, a − (x) = / 0 and α(x) = γ x in this case (similarly for a + (x) and ω(x)).An orbit γ x of p(X ) is called a separatrix of p(X ) if it is not contained in an open neighborhood U such that (U, ϕ) is parallel and such that both a ± (x) = a ± (y) for all y ∈ U and U \U consists of a + (x), a − (x) and exactly two orbits γ y and γ z such that a ± (x) = a ± (y) = a ± (z).
If X is a polynomial vector field it is known that the separatrices of p(X ) are (i) the finite and infinite singular points of p(X ); (ii) the orbits of p(X ) contained in the boundary S 1 of D; (iii) the limit cycles of p(X ); and (iv) the separatrices of the hyperbolic sectors of the finite and infinite singular points of p(X ).Moreover, if p(X ) has finitely many finite and infinite singular points and finitely many limit cycles, then p(X ) has finitely many separatrices.We call each connected component of the complement of the union of separatrices a canonical region of p(X ).Neumann (1975) proved that each canonical region of a vector field p(X ) is parallel.
To the union of the separatrices of p(X ) together with an orbit belonging to each canonical region of p(X ) we call a separatrix configuration of p(X ).We say that the separatrix configurations S 1 and S 2 of p(X 1 ) and p(X 2 ) are topologically equivalent if there exists an orientation preserving homeomorphism from D to D which transforms orbits of S 1 onto orbits of S 2 .The following is the Markus-Neumann-Peixoto classification theorem (Markus 1954, Neumann 1975, Neumann and O'Brien 1976, Peixoto 1973, Dumortier et al. 2006) for the Poincaré compactification in the Poincaré disc of polynomial systems.
Theorem 5 (Markus-Neumann-Peixoto).Let p(X 1 ) and p(X 2 ) be the Poincaré compactification of two polynomial systems ẋ = X 1 (x) and ẋ = X 2 (x), respectively.The flows of p(X 1 ) and p(X 2 ) on the Poincaré disc are topological equivalent if and only if the separatrix configurations of p(X 1 ) and p(X 2 ) are topological equivalent.
Hence, in order to qualitatively describe the phase portrait on the Poincaré disc of system (1) it is enough to qualitatively describe its separatrix configuration.This was done in Fig. 1, where we have drawn the separatrices other than singular points with bold lines.The other lines are orbits contained in its respective canonical regions.We observe from Fig. 1 that system (1) has 15 separatrices, five of them singular points, and 7 canonical regions, six of them of type strip and the one formed by the closed orbits surrounding z 0 , annular.
Below, we prove Theorem 3 by proving that Fig. 1 is a separatrix configuration of system (1).
From the previous sections we conclude that close enough to the singular points, the phase portrait of system (1) is qualitatively the one presented in Fig. 6.For further references we label the hyperbolic, parabolic and elliptic sectors presenting in the origins of the charts U 1 and V 1 in Fig. 6 as h 1 , h 2 , h 3 , h 4 , p 1 , p 2 and e 1 , e 2 , respectively.
From the definition of system (1), each of its orbits is a connected component of a level set of H F = (p 2 + q 2 )/2 (because the only singular point of this system is the center z 0 ), which in turn is the inverse image under F = (p, q) of circles surrounding the point (0, 0).Since F preserves orientation (because the Jacobian determinant of F is positive), each orbit of ( 1) is carried onto a curve contained in a circle with counterclockwise orientation.As we have seen, the curve β (s) defined in (3) is the asymptotic variety of F.Moreover, the points β (0) = (−1, −995/4) and β (1) = (0, 208) of this curve have no inverse image under F, all the other points of this curve have exactly one inverse image and the other points of R 2 have precisely two inverse images.Acting as Campbell (arXiv:math/9812032 1998), we delete from the curve β (s) the points β (0) and β (1), obtaining three curves: Campbell (1998), the inverse image under F of each C i is a curve that divides the plane into two connected components.We call D i the inverse image of C i , i = 1, 2, 3.The set D 1 ∪D 2 ∪D 3 is called by Campbell (1998) the asymptotic flower of F. It follows that R 2 \ (D 1 ∪ D 2 ∪ D 3 ) is formed by 4 connected components, each of them mapped twice onto each of the two connected components of R 2 \ {β (s)}.Each curve C i has a natural orientation, given by its parametrization (it is the opposite orientation used by Campbell (1998)).So, each curve D i also has a natural orientation (recall that F preserves orientation).The graphics of C i and D i , i = 1, 2, 3, are given in (a) and (b) of Fig. 7, respectively.As in (Campbell 1998(Campbell , 2011) ) the axes in (a) have different scales.Following Campbell (1998), we label the regions as R (right) and L (left) of the curves C i and D i .+ 432312s 2 − 86565s + 116423 , and this polynomial of degree 6 multiplying (1 − s)s(s + 1) has no real zeros by Sturm's theorem, it follows that the curves C 1 , C 2 and C 3 are transversal to the circles centered at (0, 0) = β (−1).As a consequence the curves D 1 , D 2 and D 3 are transversal to the non-singular orbits of system (1).In particular, the image of a non-periodic orbit of system (1) has α-and ω-limits contained in the curve β (s).Below, we will say that the image of an orbit starts or finishes at β (s 0 ) meaning that its α-or ω-limit is β (s 0 ), respectively.Moreover, through each point in the intersection of C 1 ∪C 2 ∪C 3 with a circle, it crosses exactly one image of an orbit of system (1).We call S 1 and S 2 the circles centered at (0, 0) and containing the points β (1) and β (0), respectively.The point z 0 , being the inverse image under F of (0, 0) = β (−1), is contained in the curve D 1 .The images under F of the closed orbits surrounding z 0 are circles surrounding (0, 0) contained in the bounded region defined by S 1 .Thus the boundary of the period annulus of the center z 0 corresponds to the arc of circle contained in S 1 , starting and finishing at the point β (1).This means that the boundary of the period annulus is an orbit that goes to infinity through the region labeled by L in (b) of Fig. 7.In particular, in the Poincaré disc, this orbit tends to the origin of the chart V 1 .Then, analyzing the possibilities in Fig. 6, we see that this orbit contains the two separatrices of the hyperbolic sector h 2 .This period annulus is an annular canonical region.
Now we analyze the parabolic sectors p 1 and p 2 .Close to the two points of D 1 cut by the orbit giving the boundary of the period annulus of the center (i.e., the orbit connecting the two separatrices of the hyperbolic sector h 2 ), and outside the period annulus, there must exist orbits cutting D 1 .Analyzing the images of these orbits, they are contained in circles surrounding the circle S 1 .So, there are two possibilities for the images of these orbits: either they are arcs starting and finishing at a point of the curve C 2 , or they are arcs starting at the curve C 3 and finishing at the curve C 2 or C 3 .At a first glance both of these possibilities are compatible with the parabolic sectors p 1 and p 2 in Fig. 6.We claim that the correct possibility is the first one.Indeed, we can increase the radii of these circles containing the images of the orbits of p 1 and p 2 until we achieve the circle S 2 .If we are in the second possibility, the orbit whose image is contained in S 2 and starts at a point of C 3 will contain the separatrice of the end of the parabolic sector p 2 .But, this orbit will not contain the separatrice of the end of the parabolic sector p 1 , because we can continue drawing arcs starting at C 3 with radii bigger than the radius of S 2 .Thus, the parabolic sector p 1 will not finish, a contradiction with the nature of the vector field at the origin of the chart V 1 , as shown in Fig. 6.This proves the claim.
So, the image of the orbits of the parabolic sectors p 1 and p 2 are arcs starting and finishing at a point of the curve C 2 .And since we can continue drawing these arcs until we arrive at circle S 2 , this means that the parabolic sectors p 1 and p 2 are connected, and the image of the orbit containing the separatrices that separate p 1 from h 1 and p 2 from h 3 is contained in the arc of S 2 starting and finishing at the point β (0).The region connecting p 1 to p 2 is a strip canonical region, see Fig. 1.Now, since the image of this last orbit cuts the curves C 3 and C 1 , there must exist orbits near it whose images cross C 3 and C 1 .The only possibility is that those images are arcs of circles starting at the curve C 1 , An Acad Bras Cienc (2018) 90 (3) PHASE PORTRAIT OF A PINCHUK MAP 2615 rotating a complete turn crossing C 3 and C 1 and continue up to finishing in the curve C 3 .We call these orbits the big orbits.A big orbit whose image is contained in a circle close enough to S 2 enters both the hyperbolic sectors h 1 and h 3 .We have to see where the big orbits start and finish.
The orbits whose images are arcs of the circles with radii smaller than the radius of S 2 , starting and finishing at C 1 and contained in the region R correspond to an elliptic sector with boundary formed by an orbit having image contained in the arc of S 2 starting at C 1 an finishing at β (0).Close to this boundary and out of the elliptic sector there must exist orbits whose images start at C 1 .These orbits are the big orbits.Hence, it follows that this elliptic sector is e 1 and that the big orbits start at the origin of V 1 , see Fig. 1.
The orbits whose images are arcs starting and finishing at C 1 and contained in the region L form the elliptic sector e 2 .Clearly its boundary is formed by the two orbits containing the separatrices of the hyperbolic sector h 4 .The image of these orbits are the arcs starting at C 1 and finishing at β (1) and starting at β (1) and finishing at C 1 , respectively.
In particular, this means that the big orbits must finish at the origin of the chart U 1 , below the hyperbolic sector h 4 .Since their images are contained in the circles bigger than S 2 , there exist orbits whose images are arcs contained in the circles between S 1 and S 2 , starting at C 1 , crossing C 2 and finishing at C 3 .These orbits produce a parabolic sector between h 3 and e 2 , and give rise to a strip canonical region as presented in Fig. 1.
The big orbits also produce a strip canonical region.
The orbits of the strip canonical region placed above the hyperbolic sector h 4 have their images contained in arcs of circles with radii bigger than the radius of S 1 , starting at C 3 and finishing at C 1 .
The elliptic sectors e 1 and e 2 form another two strip canonical regions.Hence we have 7 canonical regions, six of them are strip and one is annular.Analyzing Fig. 1, we see there are 6 finite orbits that are separatrices.The infinite has another 4 orbits.Hence, since there are 5 singular points, we have 15 separatrices in the separatrix configuration of system (1) in the Poincaré disc.

Figure 1 -
Figure 1 -The qualitative global phase portrait of system (1) in the Poincaré disc.

Figure 3 -
Figure3-The origin of system (P 5 , Q 5 ) in (a) and the origin of system (P 3 , Q 3 ) in (b).

Figure 5 -
Figure 5 -The sequence of blow downs in the study of the origin of the chart U 2 .

Figure 6 -
Figure6-The phase portrait of system (1) near the singular points.

Figure 7 -
Figure 7 -The asymptotic variety of F in (a) and the asymptotic flower of F in (b).