On the Chebyshev Property of Certain Abelian Integrals Near a Polycycle

Dumortier and Roussarie formulated in (Discrete Contin Dyn Syst 2:723–781, 2009) a conjecture concerning the Chebyshev property of a collection I0,I1,…,In\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$I_0,I_1,\ldots ,I_n$$\end{document} of Abelian integrals arising from singular perturbation problems occurring in planar slow-fast systems. The aim of this note is to show the validity of this conjecture near the polycycle at the boundary of the family of ovals defining the Abelian integrals. As a corollary of this local result we get that the linear span ⟨I0,I1,…,In⟩\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\langle I_0,I_1,\ldots ,I_n \rangle $$\end{document} is Chebyshev with accuracy k=k(n)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k=k(n)$$\end{document}.


Introduction and Statement of the Main Result
This paper is concerned with the problem of studying when a collection of Abelian integrals form an extended complete Chebyshev system (see Definition 2.1). This type of problem arises in the context of the so-called infinitesimal Hilbert's 16th problem proposed by Arnold [1]. In the present paper we are interested in a conjecture formulated by Dumortier and Roussarie [4], where the authors consider singular perturbation problems occurring in planar slow-fast systems depending on parameters. B J. Villadelprat Fig. 1 Phase portrait in the Poincaré disc of the differential system −y∂ x + (x − 1 2 y 2 )∂ y , which has H (x, y) = e −x (1 + x − 1 2 y 2 ) as a first integral They investigate the number of limit cycles that appear near a slow-fast Hopf point, i.e., its cyclicity. Their main results show that under very general conditions this cyclicity is finite and, modulo the aforementioned conjecture, provide its sharp upper bound.
In order to give a precise statement of their conjecture let us consider the function H : R 2 −→ R given by H (x, y) = e −x (1 + x − 1 2 y 2 ). It can be checked that the level sets {H (x, y) = h} for h ∈ (0, 1) are ovals γ h surrounding the origin. The family {γ h } h∈(0,1) form a period annulus and its boundary has two connected components, the parabola y 2 = 2(x + 1) and the origin (0, 0), which are the level sets h = 0 and h = 1, respectively (see Fig. 1). Let us define the family of Abelian integrals With this notation, the conjecture posed by Dumortier and Roussarie [4] is the following: Conjecture For each n 0, (I 0 , I 1 , . . . , I n ) is an ECT-system on [h 0 , 1] for any h 0 ∈ (0, 1).
The validity of this conjecture for n 2 is proved in [5]. Moreover, taking advantage of the analyticity of I k at h = 1, in that paper is also proved that for each n 0 there exists ε > 0 such that (I 0 , I 1 , . . . , I n ) is an ECT-system on (1−ε, 1], see [5,Corollary 3.5]. In other words, that the conjecture is true near the center. The aim of the present paper is to show the same property near the polycycle. More concretely our main result is the following: Theorem A For each n 0 there exists ε > 0 such that (I 0 , I 1 , . . . , I n ) is an ECTsystem on (0, ε).
Besides giving another evidence in support of the validity of the conjecture of Roussarie and Dumortier, Theorem A also implies that for each n 0 there exists a non-negative integer k = k(n) such that the linear span I 0 , I 1 , . . . , I n is Chebyshev with accuracy k in (0, 1], see Corollary 2.7. In this regard note that the validity of the conjecture is equivalent to k(n) = 0 for all n 0. For other results about the analytic properties of the Abelian integrals I k (h) the reader is referred to the paper of Françoise and Xiao [6].
The proof of Theorem A is based on an accurate study of the asymptotic properties of the Abelian integrals I k (h) at h = 0, which in the Poincaré disc is a polycycle. In the literature (see [8] and references there in) there are results that provide the expression of the coefficients in the asymptotic development of Abelian integrals near certain polycycles. However we cannot apply them because in our case the singular point at the polycycle is very degenerated. We bypass this problem by composing the Abelian integral with an infinitely flat function. Although this is certainly an ad hoc solution we hope that it can be adapted to tackle similar problems. That being said, let us stress that by no means the cyclicity of the polycycle can be deduced from Theorem A. An upper bound for the cyclicity is proved in [4], see Theorem 5.11, and the authors mention that it is only one unit higher than the sharp upper bound that they expect to hold.

Proof of the Main Result
has at most k − 1 isolated zeros on I counted with multiplicities.
To prove Theorem A we first need to show several technical lemmas dealing with some particular integrals. For reader's convenience we advance the idea of the proof in order to understand the role of these integrals. To this end note that, for each h ∈ (0, 1), the oval γ h intersects the x-axis at two points ( Thus, since Note that in principle x ± (h) are well-defined analytic functions for h ∈ (0, 1). However, due to f (−1) = 0, it follows that x − (h) can be extended analytically to h = 0 by setting x − (0) = −1. This will be a key point in our proof because it shows that, contrary to I + k , the function I − k is analytic at h = 0. To take advantage of this we will compose I k with a diffeomorphism which is infinitely flat at h = 0. More concretely, we perform the change of variable given by h = f (1/s) = e −1/s (1 + 1/s) and definê Since With this aim in view we also definê The technical lemmas that we mentioned before are concerned with the behaviour of I + k (s) as s −→ 0 + . In the statement of the first one we use the notation Taking this into account and performing the change and this proves the result.
Next, setting we get the following identity.
Proof Performing the change of variable x = 1 − e − t s in (2) we get Next we make the change u = (1 + s)x + s log(1 − x), which on account of s From this equality, an easy computation shows the validity of the result.
Indeed, the fact that M( · ; s) is an upper bound for g( · ; s) follows from the convexity of the latter and that, by construction, ∂ Since φ(s) is a flat function at s = 0, it is clear that so it is the first summand in ∂ s J α (s). Taking this into account, it is clear that the result will follow once we prove that With this aim in view we first note that, for all 0, where we set c • = c + if α − > 0 and c • = c − otherwise. Here we use that 1 − s < 1−ϕ(s), due to lim s→0 + ϕ(s) s = 0, and that consequently c − (s)x g(x; s) c + (s)x for all x ∈ (0, 1 − s). From the above inequality, using in the first summand the fact that x −→ − log(1−x) x is increasing and applying the Cauchy-Schwartz inequality in the second one, we get with One can verify that a primitive of where Q 2 is a polynomial of degree 2 with integer coefficients. On account of this we can assert that L 2 (s) 1/2 tends to zero as s −→ 0 + . Let us study next the behaviour of s L 1 (s) 1/2 . If β := 2(α − ) 0 then g β (x; s) 1, so that L 1 (s) tends to zero as and Thus, by means of an easy integration, if β = −1 then we get Cs min(β+2, 1) for some constant C > 0. (Here we use that lim s→0 + ω(s) = 1.) Consequently s L 1 (s) 1/2 Cs min(α+1, + 1 2 ) , which tends to zero as s −→ 0 + . Likewise, if β = −1 then which also tends to zero as s −→ 0 + . So far we have proved that Taking this into account, from (3) we can assert that if > 0 then lim s→0 + s J α (s) = 0, as desired.
It only remains to be proved the assertion for = 0, i.e. that lim s→0 + J α (s) = 1 1+α . To this end we claim that N − (s) J α (s) N + (s) for some functions with lim s→0 + N ± (s) = 1 1+α . The existence of N + follows from taking = 0 in (3) and (4) and the fact that lim s→0 + c • (s) = 1. Let us show the existence of N − . Setting c = c + if α ∈ (−1, 0) and c = c − if α 0, we get that Once again in the second inequality above we use the piecewise linear bounds and take 1 − s < 1 − ϕ(s) into account. Since it is clear that lim s→0 + N − (s) = 1 α+1 , the claim follows and this completes the proof of the result. f 0 , f 1 , . . . , f k−1 be analytic functions on an open interval I of R.

Definition 2.5 Let
Then The following well-known result (see for instance [9,10]) will constitute the last ingredient in the proof of Theorem A. Proof of Theorem A As we already mentioned, we shall prove that for each n 0 there exists ε > 0 such that (Î 0 ,Î 1 , . . . ,Î n ) is an ECT-system on (0, ε), where recall thatÎ k (s) := I k f (1/s) with f (x) = e −x (1 + x). This is equivalent to the assertion we want to prove because s −→ f (1/s) is a diffeomorphism from (0, +∞) to (0, 1) with lim s→0 + f (1/s) = 0. Following the notation we have introduced so far, see (1), where a k is a constant and N k is a flat function. Thus Lemma 2.2 shows thatÎ Then, by applying Lemma 2.3 and an easy computation, we can assert that where againâ k is a constant andN k is a flat function. The combination of this with Lemma 2.4 shows that we can writeÎ k (s) = (2/s) Observe that Î k (s) n k=0 is an ECT-system on (0, ε) if, and only if, s n−k L k (s) n k=0 is an ECT-system on (0, ε). (This follows by multiplying each function in the first system by s n− 1 2 .) Thanks to Lemma 2.6, the latter assertion is equivalent to show that W i (s) = 0 for all s ∈ (0, ε) and i = 0, 1, . . . , n, where .
For all k ∈ {1, 2, . . . , i + 1}, we multiply the kth row by s k and the kth column by s −k . It is clear that the value of the determinant remains unchanged after these transformations. In doing so the ( , k)-entry of the corresponding matrix is equal to which tends to k ! i−k+ 1 2 as s −→ 0 + thanks to (5). Since k = 0 in case that > k, the limit as s −→ 0 + of the corresponding matrix is upper triangular and lim s→0 + W i (s) = i k=0 k! i−k+ 1 2 = 0. Hence there exists ε > 0 such that W i (s) = 0 for s ∈ (0, ε) and i = 0, 1, . . . , n. This proves the validity of the result.
This completes the proof of the result.