arxiv: v1 [math.ds] 6 Jan 2015

Transcription

1 Global instability in the elliptic restricted three body problem Amadeu Delshams 1, Vadim Kaloshin, Abraham de la Rosa and Tere M. Seara 1 1 Universitat Politècnica de Catalunya, Amadeu.Delshams@upc.edu, Tere.M-Seara@upc.edu University of Maryland at College Park, vadim.kaloshin@gmail.com Universidad Autónoma de México, abraham.delarosa@gmail.com arxiv: v1 [math.ds] 6 Jan 15 January 7, 15 Abstract The planar) ERTBP describes the motion of a massless particle a comet) under the gravitational field of two massive bodies the primaries, say the Sun and Jupiter) revolving around their center of mass on elliptic orbits with some positive eccentricity. The aim of this paper is to show that there exist trajectories of motion such that their angular momentum performs arbitrary excursions in a large region. In particular, there exist diffusive trajectories, that is, with a large variation of angular momentum. The framework for proving this result consists on considering the motion close to the parabolic orbits of the Kepler problem between the comet and the Sun that takes place when the mass of Jupiter is zero. In other words, studying the so-called infinity manifold. Close to this manifold, it is possible to define a scattering map, which contains the map structure of the homoclinic trajectories to it. Since the inner dynamics inside the infinity manifold is trivial, two different scattering maps are used. The combination of these two scattering maps permits the design of the desired diffusive pseudo-orbits, which eventually give rise to true trajectories of the system with the help of shadowing techniques. Keywords: Elliptic Restricted Three Body problem, Arnold diffusion, splitting of separatrices, Melnikov integral. 1 Main result and methodology The planar) ERTBP describes the motion q of a massless particle a comet) under the gravitational field of two massive bodies the primaries, say the Sun and Jupiter) with mass ratio µ revolving around their center of mass on elliptic orbits with eccentricity e. In this paper we search for trajectories of motion which show a large variation of the angular momentum G = q q. In other words, we search for global instability diffusion is the term usually used) in the angular momentum of this problem. If the eccentricity vanishes, the primaries revolve along circular orbits, and such diffusion is not possible, since the planar and circular) RTBP is governed by an autonomous Hamiltonian with two degrees of freedom. This is not the case for the ERTBP, which is a +1/ degree-of-freedom Hamiltonian system with time-periodic Hamiltonian. Our main result is the following Theorem 1. There exist two constants C >, c > and µ = µ C, c) > such that for any < e < c/c and < µ < µ, and for any two values of the angular momentum in the region C G 1 < G c/e, there exists a trajectory of the ERTBP such that G) < G 1, GT ) > G for some T >. AD, AR, and TMS were partially supported by the Spanish MINECO-FEDER Grant MTM and the Catalan Grant 14SGR54. 1

2 This result will be a consequence of Theorem 1, the large C and the small constant c are explicitly computed C =, c = 1/8), and it is also shown the existence of trajectories of motion such that their angular momentum performs arbitrary excursions along the region C G 1 < G c/e. Let us recall related results about oscillatory motions and diffusion for the RTBP or the ERTBP. They hold close to a region when there is some kind of hyperbolicity in the Three Body Problem, like the Euler libration points [LMS85, CZ11, DGR1], collisions [Bol6], the infinity [LS8, Xia9, Xia9, Moe7, Rob84, MP94, MS14] or near mean motion resonances [FGKR14]. Among these papers, two were very influential for our computations, namely [LS8], the Laplace s method was used along special complex paths to compute several integrals, and [MP94], which contains asymptotic formulas for a scattering map on the infinity manifold for large values of eg. Another one, [GMS1], is very important for futute related work, since the proof of transversal manifolds of the infinity manifold is established for the RTBP for any µ, 1/]. Concerning the proof of our main result, let us first notice [LS8, GMS1] that, for a non-zero mass parameter small enough and zero eccentricity, the RTBP is not integrable, although for large G its chaotic zones have a size which is exponentially small in G. This phenomenon adds a first difficulty in proving the global instability of the angular momentum G in the ERTBP for large values of G. The framework for proving our result consists on considering the motion close to the parabolic orbits of the Kepler problem that takes place when the mass parameter is zero. To this end we study the infinity manifold, which turns out to be an invariant object topologically equivalent to a normally hyperbolic invariant manifold TNHIM). On this TNHIM, it is possible to define a scattering map, which contains the map structure of the homoclinic trajectories to the TNHIM. Unfortunately, the inner dynamics within the TNHIM is trivial, so it cannot be used combined with the scattering map to produce pseudo-orbits adequate for diffusion, and adds a second difficulty. Because of this, in this paper we introduce the use of two different scattering maps whose combination produces the desired diffusive pseudo-orbits, which eventually give rise to true trajectories of the system with the help of the shadowing results given in [GMS15]. The main issue to compute the two scattering maps consists on computing the Melnikov potential 51) associated to the TNHIM. The main difficulty for its computation comes from the fact that its size is exponentially small in the momentum G, so it is necessary to perform very accurate estimates for its Fourier coefficients. Such computations are performed in Theorem 7, and they involve a careful treatment of several Fourier expansions, as well as the computation of several integrals using Laplace s method along adequate complex paths, playing both with the eccentric and the true anomaly. To guarantee the convergence of the Fourier series, we have to assume that G is large enough G C), and e small enough Ge c). Under these two assumptions, the dominant part of the Melnikov potential consists on four harmonics, from which it is possible to compute the existence of two functionally independent scattering maps. The combination of these two scattering maps permits the design of the desired diffusive pseudoorbits, under the assumption of a mass parameter very small compared to eccentricity < µ < µ ), see 8)), which eventually give rise to true trajectories of the system with the help of shadowing techniques. It is worth noticing that since all the diffusive trajectories found in this paper shadow ellipses close to parabolas of the Kepler problem, that is, with a very large semi-major axis, their energy is close to zero, and the orientation of their semi-major axis only changes slightly at each revolution. The case of arbitrary eccentricity < e < 1 and arbitrary parameter mass parameter < µ < 1 remains open in this paper. Indeed, the case eg 1 involves the analysis of an infinite number of dominant Fourier coefficients of the Melnikov potential, as for the case eg > 1, the qualitative properties of the Melnikov function should be known without using its Fourier expansion. Larger values of the mass paratemer µ than those considered in this paper involve improving the estimates of the error terms of the splitting of separatrices in complex domains, as is usual when the splitting of separatrices is exponentially small. The computation of the explicit trajectories from the pseudoorbits found in this paper needs a suitable shadowing result given in [GMS15], which involves the translation to TNHIM of the usual shadowing techniques for NHIM. The plan of this paper is as follows. In Section we introduce the equations of the ERTBP, as

3 well as the McGehee coordinates to be used to study the motion close to infinity. In Section we recall the geometry of the Kepler problem, when the mass parameter vanishes, close to the infinity manifold and its associated separatrix. Next, in Section 4, we study the transversal intersection of the invariant manifolds for the ERTBP, as well as the scattering map associated, which depend on the Melnikov potential of the problem, whose concrete computation is deferred to Section 6. The global instability is proven in Section 5, using the computation of the Melnikov potential, and is based on the computation of two different scattering maps, whose combination gives rise to the diffusive trajectories in the angular momentum. Setting of the problem If we fix a coordinate reference system with the origin at the center of mass and call q S and q J the position of the primaries, then under the classical assumptions regarding time units, distance and masses normalization, the motion q of a massless particle under Newton s law of universal gravitation is given by d q dt = 1 µ) q S q q S q + µ q J q q J q 1) 1 µ is the mass of the particle at q S and µ the mass of the particle at q J. Introducing the conjugate momentum p = dq/dt and the self-potential function U µ q, t; e) = 1 µ q q S + µ q q J, ) equation 1) can be rewritten as a +1/ degree-of-freedom Hamiltonian system with time-periodic Hamiltonian H µ q, p, t; e) = p U µq, t; e). ) In the planar) ERTBP, the two primaries are assumed to be revolving around their center of mass on elliptic orbits with eccentricity e, unaffected by the motion q of the comet. In polar coordinates q = ρcos α, sin α), the equations of motion of the primaries are q S = µrcos f, sin f) q J = 1 µ)rcos f, sin f). 4) By the first Kepler s law the distance r between the primaries [Win41, p. 195] can be written as a function r = rf, e) r = 1 e 5) 1 + e cos f f = ft, e) is the so called true anomaly, which satisfies [Win41, p. ] df dt 1 + e cos f) =. 6) 1 e ) / Taking into account the expression 4) for the motion of the primaries, we can write explicitly the denominators of the self-potential function ) q q S = ρ µrρ cosα f) + µ r, 7) q q J = ρ + 1 µ)rρ cosα f) + 1 µ) r. 8) We now perform a standard polar-canonical change of variables q, p) ρ, α, P ρ, P α ) q = ρ cos α, ρ sin α), p = P ρ cos α P α r sin α, P ρ sin α P ) α r cos α 9) to Hamiltonian ). The equations of motion in the new coordinates are the associated to the Hamiltonian H µρ, α, P ρ, P α, t; e) = P ρ + P α ρ U µρ, α, t; e) 1)

4 with a self-potential U µ From now on we will write so that Hamiltonian 1) becomes U µρ, α, t; e) = U µ ρ cos α, ρ sin α, t; e). 11) G = P α, y = P ρ, H µρ, α, y, G, t; e) = y + G ρ U µρ, α, t; e). 1) Remark. In the planar) circular case e = RTBP), it is clear from equations 5) and 6) that r = 1 and f = t, and that the expressions for the distances 7) between the primaries depend on the time s and the angle α just through their difference α t. As a consequence, U µρ, α, s; ) as well as H µρ, α, y, G, s; ) depend also on s and α just through the same difference α t. This implies that the Jacobi constant H + G is a first integral of system..1 McGehee coordinates To study the behavior of orbits near infinity, we make the McGehee [McG7] non-canonical change of variables: ρ = x 1) for x >, which brings the infinity ρ = to the origin x =. In these McGehee coordinates, the equations associated to Hamiltonian 1) become dx dt = 1 4 x y dα dt = 1 4 x4 G the self-potential U µ is given by U µ x, α, t; e) = U µ/x, α, t; e) = x dy dt = 1 8 G x 6 x 4 dg dt = U µ α U µ x 1 µ + µ ) σ S σ J 14a) 14b) 15) with q q S = σ S = 1 µrx cosα f) µ r x 4, q q J = σ J = µ)rx cosα f) µ) r x 4. It is important to notice that the true anomaly f is present in these equations, so that the equation for f given in 6) should be added to have the complete description of the dynamics..1.1 Hamiltonian structure Under McGehee change of variables 1), the canonical form dρ dy + dα dg is transformed to ω = 4 dx dy + dα dg 16) x which, on x >, is a non-canonical) symplectic form. Therefore, expressing Hamiltonian 1) in McGehee coordinates H µ x, α, y, G, t; e) = y + x4 G U µ x, α, t; e), 17) 8 4

5 equations 14) can be written as ) dx dt = Hµ x 4 y dα dt = H µ G dy dt = x H ) µ 4 x dg dt = H µ α. 18a) 18b) Equivalently, we can write equations 18) as dz/dt = {z, H µ } in terms of the Poisson bracket ) {f, g} = x f g 4 x y f g + f g y x α G f G G α. 19) Geometry of the Kepler problem µ = ).1 The infinity manifold y h > h = h < Figure 1: Level curves of H in the x, y) plane, for fixed G > x For µ = and G >, Hamiltonian 17) becomes Duffing Hamiltonian H x, y, G) = y + x4 G U x) = y 8 + x4 G x 8 ) and is a first integral, since the system is autonomous. Moreover, H is also independent of e and α. Its associated equations are dx dt = 1 4 x y dα dt = 1 4 x4 G dy dt = 1 8 G x x4 dg dt = 1a) 1b) it is clear that G is a conserved quantity, which will be restricted to the case G > from now on, that is, G R +. The phase space, including the invariant locus x = is given by x, α, y, G) R T R R +. From equations 1) it is clear that E = {z = x =, α, y, G) R T R R + } ) is the set of equilibrium points of system 1). Moreover, for any fixed α T, G R, Λ α,g = {, α,, G)} is a parabolic equilibrium point, which is topologically equivalent to a saddle point, since it possesses stable and unstable 1-dimensional invariant manifolds. The union of such points is the -dimensional manifold of equilibrium points Λ = α,g Λ α,g. As we will deal with a time-periodic Hamiltonian, it is natural to work in the extended phase space z = z, s) = x, α, y, G, s) R T R R + T just by writing s instead of t in the Hamiltonian and adding the equation ds dt = 1 5

6 to systems 18) and 1). We write now the extended version of the invariant sets we have defined so far. For any α T, G R, the set Λ α,g = { z =, α,, G, s), s T} ) is a π-periodic orbit with motion determined by ds/dt = 1. The union of such periodic orbits is the -dimensional invariant manifold the infinity manifold) Λ = α,g Λ α,g = {, α,, G, s), α, G, s) T R + T} T R + T, 4) which is topologically equivalent to a normally hyperbolic invariant manifold TNHIM). Parameterizing the points in Λ by x = x α, G, s) = x α, G), s) =, α,, G, s) Λ T R + T the inner dynamics on Λ is trivial, since it is given by the dynamics on each periodic orbit Λ α,g : φ t, x ) =, α,, G, s + t) = x α, G), s + t) = x α, G, s + t). 5). The scattering map In the region of the phase space with positive angular momentum G, let us now look at the homoclinic orbits to the previously introduced invariant objects. The equilibrium points Λ α,g have stable and unstable 1-dimensional invariant manifolds γ α,g = W u Λ α,g ) = W s Λ α,g ) { = z = x, ˆα, y, G), H x, y, G) =, ˆα = α G H = } x y dx, as the -dimensional manifold of equilibrium points Λ has stable and unstable -dimensional invariant manifolds which coincide and are given by The surface γ = W u Λ ) = W s Λ ) = {z = x, α, y, G), H x, y, G) = }. γ α,g = W u Λ α,g ) = W s Λ α,g ) 6) { } x = z = x, ˆα, y, G, s), s T, H x, y, G) =, ˆα = α G y dx 7) is a -dimensional homoclinic manifold to the periodic orbit Λ α,g in the extended phase space. The 4-dimensional stable and unstable manifolds of the infinity manifold Λ coincide along the 4-dimensional homoclinic invariant manifold the separatrix), which is just the union of the homoclinic surfaces γ α,g : H = γ = W u Λ ) = W s Λ ) = α,g γ α,g = { z = x, α, y, G, s), α, G, s) T R + T, H x, α, y, G) = } 8) Due to the presence of the factor x /4 in front of equations 1), it is more convenient to parameterize the separatrix γ α,g given in 6) by the solutions of the Hamiltonian flow contained in H = in some time τ satisfying see [MP94]) dt dτ = G x. 9) 6

7 In this way, the homoclinic solution to the periodic orbit Λ α,g of system 1) can be written as x t; G) = G1 + τ ) 1/ a) α t; α, G) = α + π + arctan τ τ y t; G) = G1 + τ ) G t; G) = G s t; s) = s + t α and G are free parameters and the relation between t and τ is t = G τ + τ b) c) d) e) ), 1) which is equivalent to 9) on H. From the expressions above, we see that the convergence along the separatrix to the infinity manifold is power-like in τ and t: x, y, α α π G Gτ 6t, τ, t ±. ) We now introduce the notation z = z σ, α, G, s) = z σ, α, G), s) = x σ; G), α σ; α, G), y σ; G), G, s) γ ) so that we can parameterize any surface γ α,g as γ α,g = { z = z σ, α, G, s) = z σ, α, G), s), σ R, s T}. and we can parameterize the 4-dimensional separatrix as γ = W Λ ) = { z = z σ, α, G, s) = z σ, α, G), s), σ R, G R +, α, s) T }. 4) The motion on γ is given by φ t, z ) = z σ + t, α, G, s) = z σ + t, α, G), s + t) 5) and by equations ), 1) the following asymptotic formula follows: φ t, z ) φ t, x ) = z σ + t, α, G), s + t) x α, G), s + t) t ± 6) The scattering map S describes the homoclinic orbits to the infinity manifold Λ defined in 4)) to itself. Given x, x + Λ, we define S µ x ) := x + if there exists z W u µ Λ ) W s µ Λ ) such that φ t,µ z ) φ t,µ x ± ) for t ±. 7) In the case µ = the asymptotic relation 6) implies S x ) = x so that that the scattering map S : Λ Λ is the identity. 7

8 4 Invariant manifolds for the ERTBP µ > ) 4.1 The infinity manifold In order to analyze the structure of system 18), we will write H µ given in 17) as we have written U µ in 15) as H µ x, α, y, G, s; e) = H x, y, G) µ U µ x, α, s; e) 8) U µ x, α, s; e) = U x) + µ U µ x, α, s; e) = x + µ U µx, α, s; e), 9) and we proceed to study the dynamics as a perturbation of the limit case µ =. From 15), U x, α, s; e) = lim µ U µ x, α, s; e) x x ) r = [ 4 + x4 r + 4x r cosα f) ] + x 1/ cosα f) 4) r = rf, e) and f = fs, e) are given, respectively, in 5-6). For µ >, it is clear from equations 18) that the set E remains invariant and, therefore, so does the infinity manifold Λ, being again a TNHIM, and all the periodic orbits Λ α,g also persist. The inner dynamics on Λ is the same that in the case µ =, so that the parametrization x as well as its trivial dynamics 5) remain the same. 4. The scattering map From [McG7] we know that W s µ Λ ) and W u µ Λ ) exist for µ small enough and are 4-dimensional in the extended phase space. The existence of a scattering map will depend on the transversal intersection between these two manifolds. Let us take an arbitrary z = z, s) = z σ, α, G), s) γ as in ). Now, we have to construct points in W s µ Λ ) and W u µ Λ ) to measure the distance between them. It is clear from the definition of γ that ṽ = H z ), ) is orthogonal to γ = W u Λ ) = W s Λ ) at z and then if the normal bundle to γ is denoted by N z ) = { z + σ ṽ, σ R} we have that there exist unique points z s,u µ = z s,u µ, s) such that { z s,u µ } = W s,u µ Λ ) N z ). 41) The distance we want to compute between W s µ Λ ) and W u µ Λ ) is the signed magnitude given by d z, µ) = H z u µ) H z s µ). 4) We now introduce the Melnikov potential see [DG, DLS6]) Lα, G, s; e) = U x t; G), α t; α, G), s + t; e) dt 4) U is defined in 4). Thanks to the asymptotic behavior ) of the solutions along the separatrix and of the self potential close to the infinity manifold U x, α, s; e) = Ox 4 ) as x this integral is absolutely convergent, and will be computed in detail in Section 6. 8

9 Proposition. Given α, G, s) T R + T, assume that the function σ R Lα, G, s σ; e) R 44) has a non-degenerate critical point σ = σ α, G, s; e). Then for < µ small enough, close to the point z = z σ, α, G), s) γ see the parameterization in )), there exists a locally unique point z = z σ, α, G, s; e, µ) W s µ Λ ) W u µ Λ ) of the form z = z + Oµ). Also, there exist unique points x ± =, α ±,, G ±, s) =, α,, G, s) + Oµ) Λ such that Moreover, we have φ t,µ z ) φ t,µ x ± ) for t ±. 45) G + G = µ L α α, G, s σ α, G, s; e)) + Oµ ). 46) Proof. From equation ) we know that any point z γ have the form As in 41), we consider z = z σ, α, G, s). z s,u µ = z s,u µ, s) W s,u µ Λ ) N z ), and we are looking for z such that z s µ = z u µ. There must exist points x ± = z ±, s) Λ such that φ t,µ z s,u µ ) φ t,µ x ± ), 47) t ± moreover φ t,µ z s,u µ ) φ t, z ) = Oµ) for ±t see [McG7]). Since H does not depend on time, by 8) and the chain rule we have that d dt H φ t,µ z s,u µ )) = {H, H µ }φ t,µ z s,u µ )) = µ{h, U µ }φ t,µ z s,u µ ); e). Since H = in Λ, using 47) and the trivial dynamics on Λ we obtain H z s,u µ ) = µ ± Taylor expanding in µ and using the notation 5) H z u µ) H z s µ) = µ On the other hand, from 4) and then so that = µ Lα, G, s; e) = L α, G, s; e) = s L α, G, s σ; e) = s = {H, U µ }φ t,µ z s,u µ ); e) dt. {H, U }φ t, z ); e) dt + Oµ ) {H, U }z σ + t, α, G), s + t; e) dt + Oµ ). 48) U x ν s; G), α ν s; α, G), s; e) dν { U, H }z ν s, α, G), s; e) dν {H, U }z ν s + σ, α, G), s; e) dν {H, U }z t, α, G), s + t; e) dt 49) 9

10 and therefore, from 48) and 49) d z, µ) = H z u µ) H z s µ) = µ L s α, G, s σ; e) + Oµ ). For µ small enough, it is clear by the implicit function theorem that a non degenerate critical value σ of the function 44) gives rise to a homoclinic point z to Λ the manifolds Wµ Λ s ) and Wµ u Λ ) intersect transversally and has the desired form z = z + Oµ). Consider now the solution of system 18) represented by φ t,µ z ). By the fundamental theorem of calculus and 8) we have G + G = = µ = µ H µ α φ t,µ z )) dt = U α φ t, z ); e) dt + Oµ ) U µ α φ t,µ z ); e) dt U α z σ + t, α, G), s + t; e) dt + Oµ ) = µ L α α, G, s σ ; e) + Oµ ). Once we have found a critical point σ = σ α, G, s; e) of 44) on a domain of α, G, s), we can define the reduced Poincaré function see [DLS6]) L α, G; e) := Lα, G, s σ ; e) = Lα, G, s ; e) 5) with s = s σ. Note that the reduced Poincaré function does not depend on the s chosen, since by Proposition s L α, G, s σ α, G, s; e); e)). Note also that if the function 44) in Proposition has different non degenerate critical points there will exist different scattering maps. The next proposition gives an approximation of the scattering map in the general case µ > Proposition 4. The associated scattering map α +, G +, s + ) = S µ α, G, s) for any non degenerate critical point σ of the function defined in 44) is given by α, G, s) ) α µ L G α, G; e) + Oµ ), G + µ L α α, G; e) + Oµ ), s L is the Poincaré reduced function introduced in 5). Proof. By hypothesis we have a non degenerate critical point σ of 44). By definition 5), Proposition gives G + G = µ L α α, G) + Oµ ). as well as G = G + Oµ) to get the correspondence between G + and G that were looking for. The companion equation to 46) α + α = µ L G α, G) + Oµ ) is a direct consequence of the fact that the scattering map S µ is symplectic. Indeed, this is a standard result for a scattering map associated to a NHIM, and is proven in [DLS8, Theorem 8]. For what concerns our scattering map defined on a TNHIM, the only difference is that the stable contraction expansion) along Wµ s,u Λ ) is power-like ) instead of exponential with respect to time. Therefore we only have to check that Proposition 1 in[dls8] still 1

11 holds, namely that Area φ t,µ R)) when t for every -cell R in W s µ Λ ) parameterized by R : [, 1] [, 1] W s µ Λ ) in such a way that Rt 1, t ) W s µ Λ ), R, t ) Λ. But this is a direct consequence of the fact that the stable coordinates contract at least by C/ t see )) and the coordinates along Λ do not expand at all. Remark 5. In the planar) circular case e = RTBP), U µ x, α, s; e) depends on the time s and the angle α just through their difference α s, see Remark. From one readily obtains and therefore U µ α x, α, s; ) = U µ x, α, s; ) s L s α, G, s; ) = Lα, G, s; ) α L α α, G, s σ ; e) = L s α, G, s σ ; ) = and consequently the reduced Poincaré function L does not depend on α, and G + = G + Oµ ). But indeed G + G in the circular case, since there exists the first integral provided by the Jacobi constant C J = H µ + G and as H µ = on Λ, G + = G. Therefore in the circular case there is no possibility to find diffusive orbits studying the intersection of W s µ Λ ) and W u µ Λ ) since any scattering map preserves the angular momentum. 5 Global diffusion in the ERTBP We have already the tools to derive the scattering maps to the infinity manifold Λ, namely Proposition to find transversal homoclinic orbits to Λ and Proposition 4 to give their expressions. Both of them rely on computations on the Melnikov potential L. Inserting in the Melnikov potential introduced in 4) the expression for U in 4) we get Lα, G, s; e) = [ x [ 4 + x 4 r + 4x r cosα f) ] 1/ + x ) r cosα f) x ] dt 51) x and α, coordinates of the homoclinic orbit defined in ), are evaluated at t, as r and f, defined in 5) and 6), are evaluated at s + t. To evaluate the above Melnikov potential, we will compute its Fourier coefficients with respect to the angular variables α, s. Since x and r are even functions of t and f and α are odd, L is an even function of the angular variables α, s: L α, G, s; e) = Lα, G, s; e), and therefore L has a Fourier Cosine series with real coefficients L q,k : L = L q,k e iqs+kα) = L, + L,k cos kα + L q,k cosqs + kα). 5) q Z k Z k 1 q 1 k Z The concrete computation of the Fourier coefficients of the Melnikov potential 51) will be carried out in section 6. First, some accurate bounds will be obtained: 11

12 Lemma 6. Let G, q 1, k and l. Then L q,l B q,l and L,l B,l, B q, = 9 q e q e q G / e qg / B q,1 = 7 e q 1 + e) 4 G 7/ e qg / B q, 1 = 9 q e q e 1 q G 1/ e qg / B q,k = 5 k e q 1 + e) k G k 1/ e qg / B q, k = 5 q+k e q e k q G k 1/ e qg / B,l = 8 l e l G l 5) Directly from this lemma, we first see that the harmonics L q,l are exponentially small in G for q 1, so it will be convenient to split the Fourier expansion 5) as L = L + L 1 + L + = L + L 1 + L 54) L α, G; e) = L, + k 1 L,k cos kα, L q α, G, s; e) = k Z L q,k cosqs + kα), q 1. 55) The function L does not depend on the angle s and it contains the harmonics of L of order in s, which are of finite order in terms of G, L 1 the harmonics of first order, which are of order e qg /, and all the harmonics of L q for q are much exponentially smaller in G than those of L 1, so we will estimate L and L 1 and bound L. To this end, it will be necessary to sum the series in 55). From the bounds B q,k in 5) for the harmonics L q,k we get the quotients B q,k e) = B q,k G for k, B q, k+1) B q, k = 4eG for k q, B,l+1 B q,l = 4e for l, 56) G which indicate that, for fixed q, we will need at least the conditions G > 1 + e) and eg < 1/4 to ensure the convergence of the Fourier series. This is the main reason why we are going to restrict ourselves to the region G C large enough and eg c small enough along this paper to get the diffusive orbits. Among the harmonics L,k of order in s, by 56), the harmonic L, appears to be the dominant one, but we will also estimate L,1 to get information about the variable α, and bound the rest of harmonics L,k for k. Among the harmonics of first order L 1,k, again by 56), the five harmonics L 1,k for k are the only candidates to be the dominant ones, but the quotients from 5) B 1, 1 + e) = B 1, 1 8eG 4, B 1,1 1 + e)4 = B 1, 1 8eG, B 1, = e B 1, 1 G = eg G, 57) indicate that L 1, 1 and L 1, appear to be the two dominant harmonics of order 1. Summarizing, to compute the series 5) we estimate only the four harmonics L,, L,1, L 1, 1 and L 1,, and bound all the rest, providing the following result, whose proof will also be carried out in section 6. Theorem 7. For G, eg 1/8, the Melnikov potential 51) is given by with Lα, G, s; e) = L α, G; e) + L 1 α, G, s; e) + L α, G, s; e) 58) L α, G; e) = L, + L,1 cos α + E α, G; e) 59) L 1 α, G, s; e) = L 1, 1 coss α) + L 1, coss α) + E 1 α, G, s; e) 6) 1

13 the four harmonics above are given by and the error functions satisfy L, = π G 1 + E,) L,1 = 15πe 8G E,1) π L 1, 1 = 8G e G / 1 + E 1, 1 ) 61) L 1, = πeg / e G / 1 + E 1, ) 6) E, 1 G e E,1 1 G 4 + e E G e E 1, 17 G e E 14 e G 7 E 1 19 e G / [ G 7/ + e G 5/ + eg /] 6) L 8 G / e G / 64) The function L 1 contains only harmonics of first order in s, so we can write it as a cosine function in s. Introducing in the definition 6) of L 1 we can write p := L 1, L 1, 1 = 1eG 1 + E 1, E 1, =: 1eG 1 + E p ) 65) L 1,k L 1 = L 1, 1 cosqs + kα) L 1, 1 k Z = L 1, 1 coss α) p coss α) + = L 1, 1 R e is α) 1 pe iα + ) k 1, L 1,k L 1,k cosqs + kα) L 1, 1 e ik+1)α L 1, 1 k 1, = L 1, 1 R e is α) Be iθ) = L 1, 1 B coss α θ) 66) B = Bα, G; e) and θ = θα, G; e) [ π, π) are the modulus and the argument of the complex expression 1 pe iα + k 1, L 1,k L 1, 1 e ik+1)α =: Be iθ. 67) Writing also in polar form the quotient of the sum in 67) by the parameter p introduced in 65) Ee iφ := k 1, L 1,k pl 1, 1 e ik+1)α = k 1, L 1,k L 1, e ik+1)α, with E = Eα, G; e) and φ = φα, G; e) [ π, π), equation 67) for B and θ reads now as Be iθ = 1 pe iα + pee iφ 68) 1

14 or, equivalently, as the couple of real equations B cos θ = 1 p cos α + pe cos φ 69) B sin θ = p sin α pe sin φ. 7) The function E = Eα, G; e) is small, since, by 6), 6) and 65), E E 1 L 1, = E 1 pl 1, e)4 G 7/ + e G 5/ + eg / 8G/π πeg / = π π1 + e)4 G ) 1 G + eg ) = O G, eg ), 71) with an analogous bound for its derivative with respect to α. Writing equation 68) as Be iθ = Be i θ + pee iφ one gets the explicit formulae for B and θ B = 1 p cos α + p = 1 p) + 4p sin α, 7) ) p sin α θ = arctan π, π] B + 1 p cos α from which we see that B behaves like a distance to the point p = 1 and α =. The angle θ is not well defined when B =, but this happens only for α = and p = 1, that is, for G 1e) 1/. A totally analogous property holds for B: Lemma 8. Bα, G; e) > except for α = and p = 1 + k 1, L 1,k/L 1, 1. Proof. For B =, equation 7) reads as sin α = fα). 7) fα) = fα, G; e) := E sin φ = E sin φα, G; e). Since f + f/ α) < 1 due to 71), there are exactly two simple solutions of equation 7) in the interval [ π/, π/]; one is α,+ π/, π/) obtained as a fixed point of the contraction α = arcsin fα, G; e)), and a second α, π/, π/) fixed point of the contraction α = π arcsin fα, G; e)). Taking a closer look at equation 68), we see that if α changes to α, then φ, θ, B are solutions of 68) or, in other words, φ, θ are odd functions of α and B even. Therefore α =, π are the unique solutions of equation 7) for B =. Substituting α =, π in 69) for B =, only α = provides a positive p, which is then given by p = 1 + pe = 1 + k 1, L 1,k/L 1, 1. We are now in position to find critical points of the function s Lα, G, s; e). To this end we will check that s Lα, G, s; e) is indeed a cosine-like function, that is, with a non-degenerate maximum minimum) and no other critical points. By Theorem 7, The dominant part of the Melnikov potential L is given by L + L 1. By equation 58) and the bounds for the error term, for G large enough, the critical points in the variable s are well approximated by the critical points of the function L + L 1 and therefore will be close to s α θ =, π mod π) thanks to expression 66). For this purpose, we introduce L 1 = L 1α, G; e) = L 1, 1 B 74) B = Bα, G; e) is given in 67) and L 1, 1 is the harmonic computed in 61). With this notation the function L 1 can be written as a cosine function in s L 1 α, G, s; e) = L 1α, G; e) coss α θ), 75) 14

15 and differentiating the Melnikov potential 58) with respecto to s we get L s = L 1 sins α θ) + L = sins α θ) = 1 L s L 1 s 76) which is a equation of the form 7) for s α θ instead of α and f = L / s) /L 1. Therefore, as long as B > L / s) /L 1, 1 ), which by the estimate 61) ) for L 1, 1 and the bound 64) for L happens outside of a neighborhood of size O G / e G / of the point α =, G = G ) G 1e) 1/ is such that p = 1 + k 1, L 1,k /L 1, 1, 77) there exist exactly two no-degenerate critical points s,± of the function s Lα, G, s; e). Let us recall now that the Melnikov function L, as well as its terms L q are all expressed as Fourier Cosine series in the angles α and s, or equivalently, they are even functions of α, s). Consequently, L q / s is an odd function of α, s), and it is easy to check that each critical point s is an odd function of α. Moreover, using the Fourier Sine expansion of L q / s, one sees that it s is a critical point of s Lα, G, s; e), s + π too, so s, = s,+ + π. We state all this in the following proposition. Proposition 9. Let L be the Melnikov potential ) given in 58), G and eg 1/8. Then, except for a neighborhood of size O G / e G / of the point α =, G = G ) given in 77), s Lα, G, s; e) is a cosine-like function, and its critical points are given by s,+ = s,+α, G; e) = α + θ + ϕ +, s, = s,+ + π = α + θ + π + ϕ ± ) θ = θα, G; e) is given in 67) and ϕ + = O G / e G /. From the proposition above we know that there exist s, and s,+, non degenerate critical points of s Lα, G, s; e). Therefore, we can define two different reduced Poincaré functions 5) L ±α, G; e) = Lα, G, s,±; e) = L α, G; e) ± L 1α, G; e) + E ± α, G; e). By the symmetry properties of L q, it turns out that each L qα, G; e) = Lα, G, t,; e) is an even function of α. Moreover, since s, = s,+ + π, one has that L q = 1) q L q, so we can write the reduced Poincaré map as L ± = L ± L 1 + L ± L + L 4 ± 78) From the expression for the scattering map given in Proposition 4 we can define two different scattering maps, namely S ± α, G, s) = α + µ L ± G α, G; e) + Oµ ), G µ L ) ± α α, G; e) + Oµ ), s. 79) These two scattering maps are different since they depend on the two reduced Poincaré-Melnikov potentials L ±. From their expression 79), the scattering maps S ± follow closely the level curves of the Hamiltonians L ±. More precisely, up to Oµ ) terms, S ± is given by the time µ map of the Hamiltonian flow of Hamiltonian L ±. The Oµ ) remainder will be negligible as long as µ < L ± G, L ± α. Nevertheless, since we want to switch scattering maps, we will need to impose ) µ < L 1 = L 1, 1 B = O G 1/ e G /, 15

16 that is, µ exponentially small with respect to G in the region C G c/e which a fortiori is satisfied for < µ < µ = e c/e) /. 8) This is the relation between the eccentricity and the mass parameter that we need to guarantee that our main result holds. This kind of relation is typical in problems with exponential splitting, when the bound of the remainder, here Oµ ), is obtained through a direct application of the Melnikov method for the real system. To get better estimates for this remainder, one needs to bound this remainder for complex values of the parameter t or τ of the parameterization ) of the unperturbed separatrix. Such approach has recently been used for in the RTBP in [GMS1] and it is likely to work in the ERTBP, allowing us to consider any µ, 1/], that is, imposing no restrictions on the mass parameter. We want to show now that the foliations of L ± = constant are different, since this will imply that the scattering maps S ± are different. Even more, we will design a mechanism in which we will determine the places in the plane α, G) we will change from one scattering map to the other, obtaining trajectories with increasing angular momentum G. To check that the level curves of L + and L are different, and indeed transversal, we only need to check that their Poisson bracket is zero. Since L + and L are even functions of α, their Poisson bracket {L +, L +} will be an odd function of α, so we already know that it will have a factor sin α. Using equation 78) we can write {L +, L +} = {L + L 1 + L +, L L 1 + L } = {L, L 1} + E E contains only Poisson brackets of odd order E = {L, L } + {L 1, L }) [q/] q,q odd 5 q= {L q, L q q }. Therefore, by the bounds 5) for the harmonics L q,k, the error term E = O e G) is much ) exponentially smaller for large G than {L, L 1}, which is O e G / and we now compute. By splitting L, using 59), and L 1 = L 1, 1 B, using 61), 7) and 71), in their dominant and non-dominant parts L = L + E, L 1, 1 = L 1, E,1 ), B = B + E B, after a straightforward computation, we arrive at {L, L 1} = { L, L 1} + E J with d = [ {L, L 1} = L 1 B eg G cos α 5 48 B G πp sin α G 4 d [ G cos α + p B 4eG ] ] G. and a small error term E J = O G 5 + eg + e G + pe G peg + G 6 ) )) G 1/ e G / 5.1 Strategy for diffusion The previous lemma tells us that the level curves of L + and L are transversal in the region G and eg 1/8, except for the three curves α =, α = π and d = which, by the way, are also transversal to any of these level curves, see figure ). 16

17 Figure : Level Sets of L + L ) in Blue Red) and d= in Green Thus, apart from these curves, at any point in the plane α, G) the slopes dg/dα of the level curves of L + and L are different, and we are able to choose which level curve increases more the value of G, when both slopes are positive, or alternatively, to choose the level curve which decreases less the value of G, when both slopes are negative see Figure ). In the same way, we can find trajectories along which the angular momentum performs arbitrary excursions. More precisely, given an arbitrary finite sequence of values G i, i = 1,..., n we can find trajectories which satisfy GT i ) = G i, i = 1,..., n. Strictly speaking, This mechanism given by the application of scattering maps produce indeed pseudo-orbits, that is, heteroclinic connections between different periodic orbits in the infinity manifold which are commonly known as transition chains after Arnold s pioneering work [Arn64]. The existence of true orbits of the system which follow closely these transition chains relies on shadowing methods, which are standard for partially hyperbolic periodic orbits the so-called whiskered tori in the literature) lying on a normally hyperbolic invariant manifold NHIM). Such shadowing methods are equally applicable in our case, we have an infinity manifold Λ which is only topologically equivalent to a NHIM see [Rob88, Rob84, Moe, Moe7, GL6, GLS14]) and [GMS15]. With all these elements, we can finally state our main result Theorem 1. Let G 1 < G large enough and e >, µ > small enough. More precisely G 1 < G 1/8e) and < µ < µ = e 8e) /. Then, for any finite sequence of values G i G 1, G ), i = 1,..., n, there exists a trajectory of the ERTBP such that GT i ) = G i, i = 1,..., n for some < T i < T i+1. In particular, for any two values G 1 < G G 1, G ), there exists a trajectory such that G) < G 1, and GT ) > G for some time T >. 17

18 Figure : Zone of diffusion: Level curves of L + L ) in blue red) and diffusion trajectories in green. 6 Computation of the Melnikov potential: Proof of Theorem 7 The main difficulty to compute the Melnikov potential is that it is given by an integral 51) the coordinates of the separatrix x and α are given implicitly ) in terms of the time t through the variable τ 1), as r and f are given in terms of s+t through the differential equation 6) defining the true anomaly f. To evaluate the above Melnikov potential, we will compute its Fourier Cosine series 5) in the angles s, α. The next proposition gives formulas for its Fourier coefficients. To this we will consider the Fourier expansion of the functions: rft)) n e imft) = q Z c n,m q e iqt 81) which can be found in [MP94] and [Win41, p. 4]. Using that r is an even function and that f is and odd function, one readily sees that the above coefficients are real and indeed they satisfy c n, m q = c n,m q = c n,m q. 8) Once these coefficients c n,m q are introduced we can give explicit formulas for the Fourier coefficients of the Melnikov potential L. Proposition 11. The Melnikov potential given in 51) or in 5) can be written as L = q Z L q e iqs with L q = k Z L q,k e ikα. 8) 18

19 with L q, = l 1 L q,1 = l L q, 1 = l L q,k = l k L q, k = l k c l, q Nq, l, l) 84a) c l 1, 1 q Nq, l 1, l) 84b) c l 1,1 q Nq, l, l 1) 84c) c l k, k q Nq, l k, l) for k 84d) c l k,k q Nq, l, l k) for k 84e) and Nq, m, n) = m+n G m+n 1 ) ) 1/ 1/ m n e iqτ+τ /) G / τ i) m dτ 85) τ + i) n Proof. We write Melnikov potential 51) as: [ ] x L = L ) 1 + r cosα f) x dt, 86) L 1 = x [ 4 + x 4 r + 4x r cosα f) ] 1/ dt can be written as L 1 = Using we get L l k = S l k = x 1 l k+1 1 l+k x rft + s))eiα ft+s)) ) ) 1/ 1/ l k l ) ) 1/ 1/ k l l ) 1/ ) 1 + z) 1 1/ = z l l l= L 1 = L l k + k l k k< l k x 4l k+ ) 1/ 1 + x rft + s))e iα ft+s)) dt. 87) S l k [rft + s))] l k e ikα e ikft+s) dt; k l x 4l+k+ [rft + s))] l+k e ikα e ikft+s) dt; l k 1. With these expressions is easy to see that L cancels out the last term in the integral 86) and that L cancels the cosine term, and so S 1 1 L = l 1 L l + l Now we perform the change of variable t = G L l 1 + τ + τ l S l 1 + k>1 L l k + S k. l 88) l k ), dt = G 1 + τ )dτ 19 k< 1 l k

20 introduced in 1), and we use the formulas for x and α given in a) and b). In particular we will use that: x 4 = G 1 + τ ) = Gdτ, eiα = τ i τ + i eiα, and the expansion in Fourier series given in 81), obtaining L l k = e ikα = e ikα q Z S l k = e ikα = e ikα q Z l k G 4l k 1 ) ) 1/ 1/ l k l q Z e iq s cq l k, k e iqτ+τ /) G / dτ; k l τ i) l k) τ + i) l e iq s c l k, k q Nq, l k, l) 89a) l+k G 4l+k 1 1/ l ) 1/ k l ) q Z e iq s cq l+k, k e iqτ+τ /) G / dτ; l k 1 τ i) l τ + i) k l) e iq s c l+k, k q Nq, l, k l) 89b) substituting now equations 89a) and 89b) in the expansion 88) we get L = q Z e iqs l 1 c l, q Nq, l, l) + e iqs+α) q Z l + e iqs+kα) q Z k l k c l 1, 1 q Nq, l 1, l) + e iqs α) q Z l c l k, k q Nq, l k, l) + e iqs+kα) q Z k l k c l 1,1 q Nq, l, l 1) c l+k, k q Nq, l, k l) Now changing the indexes l l and k k in the third and fifth terms one obtain the desired formulas for the Fourier coefficients L q,k. In view of proposition 11) and formulas 84), to compute the dominant part of the Melnikov potential and obtain effective bounds of the errors we will need to estimate the constants c n,m q defined in 81) and the integrals Nq, m, n) defined in 85) for q and only for indices m, n satisfying n, m n + 1. Alternatively to 5), it will be very convenient to express the distance r between the primaries as r = 1 e cos E 91) in terms of the eccentric anomaly E, given by the Kepler equation [Win41, p. 194] This is done in the next three propositions. 9) t = E e sin E. 9) Proposition 1. Let n, m, q Z, n, q, m n + 1. Then the Fourier coefficients c n,m q in 81) satisfy { c n,m q+n+1 e q 1 e e m q m q 1 + e) n+1 m 1 Proof. In the integral formula for the Fourier coefficients c n,m q = 1 π π defined r n e imf e iqt dt 9) we change the variable of integration to the eccentric anomaly 9) dt = rde) to get c n,m q = 1 π π re if ) m r n+1 m e iqt de. 94)

21 To compute c n,m q from 94) we will use the identity see [Win41, p. ]) re if ) n = ae ie/ e 1 + e + 1 e a e ie/, a = which readily implies a + e 4a = 1, re if = a e ie e + e 4a e ie = ae ie/ e a e ie/ ) a e 4a = 1 e, a 4 + e 16a 4 = 1 e, a 4 e 16a 4 = 1 e. 95a) 95b) To bound the integral 9) for m we will consider two different cases: q m and m < q. Let us first consider the case q m. By the analyticity and periodicity of the integral we change the path of integration from IE) = to IE = lna /e): ) a E = u + i ln u [, π] e so that and then, by 95a), 95b), and 9) therefore e ie = e iu lna /e) = e a eiu re if = e eiu e + e e iu = ecos u 1) r = 1 e ) e a eiu + a e e iu = 1 e 4a eiu a e iu ) ) e e = 1 4a + a cos u + i 4a a sin u = 1 cos u + i 1 e sin u e it = a e iu ) e exp e 4a eiu a e iu = a e iu exp ) 1 e e cos u + i sin u re if = e1 cos u) e r = 1 cos u) + 1 e ) sin u = 1 cos u) e sin u e it = a e exp ) 1 e cos u a e e 1 e. Since a, substituting these bounds in 94) we find directly the desired result for q m. For the the case m < q we now perform the change of the integration variable through ) a E = v i ln, v [, π] e so that e ie = e iv+lna /e) = a e eiv and then, by 95a), 95b), and 9) ) re if = a a4 4 + e4 cos v e = e e eiv e + e 8a 4 e iv = e 16a 4 cos v e 1 + cos v) + i 1 e sin v ) a + i 4 e4 16a 4 r = 1 e 4a e iv a e iv = 1 cos v i 1 e sin v ) ) e it = a ee iv a exp e iv e 4a e iv = ee iv a exp 1 e cos u + i sin u ) ) sin v 1

22 therefore as before re if, r e it e a e 1 e. Since a 1, substituting these bounds in 94) we find the desired result for m < q. For m 1 we bound directly the integral 94) for E [, π]. Since e if = e it = 1 we have c n,m q 1 π π r n+1 de by noticing that r 1 + e) we conclude the proof of the bounds for the c n,m q. As we can see from equations 84) the Fourier coefficients of the Melnikov potential L depend on the function Nq, m, n) defined in 85), so to bound or to compute) these Fourier coefficients we need to bound or to compute) Nq, m, n). Introducing Nq, m, n) can be written as Iq, m, n) = Nq, m, n) = m+n G m+n 1 e iq G τ+ τ ) τ i) m τ + i) ) 1/ 1/ m n n dτ ) Iq, m, n). We will call τ ) hτ) = i + τ the variable term in the exponencial of the integral, so that 96) Iq, m, n) = e q G hτ) τ i) m dτ. 97) τ + i) n Since the integral Iq, m, n) involves an exponential, it will be useful the Laplace s method [Erd56] of integration. In particular on a complex path with Ihτ)) =. So, let us define the path Γ = Γ 1 Γ Γ Γ 4 Γ 5 98) < ε < 1, τ is a point such that Ihτ )) = that will be fixed later, and Γ 1 = {τ C Ihτ)) = } {τ C Rτ) R τ )} Γ 5 = {τ C Ihτ)) = } {τ C Rτ) R τ )} Γ = {τ C Ihτ)) = } {τ C R τ ) Rτ) } {τ C τ i cε} Γ 4 = {τ C Ihτ)) = } {τ C Rτ) R τ )} {τ C τ i cε} Γ = {τ C Ihτ)) } {τ C τ i = cε}. 99) τ Plane Γ 1 Γ 5 τ Γ Γ 4 C τ I [hτ)] i C Γ i c ε By means of the Cauchy-Goursat theorem and a limit argument, the integral Iq, m, n), defined in 97) over the real axis, is equal to the one taken over the path Γ thinking of τ as a complex number see [LS8]). In fact, by the same argument, its value does not depend on ε. The positive branch of the hyperbola defined by Ihτ)) = intersects the circumference of radius ε in two points that can be expressed as C ε and C ε defined by: C ε = Γ Γ 4 C ε = Γ Γ 1) Figure 4: The complex path Γ

23 Since the integral over Γ does not depend on ε, we will choose a particular value of ε to bound Iq, m, n) and consequently Nq, m, n) defined in 85). Later on, in proposition 17 we will just compute the ε-independent terms of this integral. It is not difficult to see that, if we define the function then uτ) = hi) hτ) = i τ + τ ) = τ i) i τ i), 11) uγ 1 Γ ), uγ 4 Γ 5 ) R +. Moreover, if τ Γ 1 Γ then τ + = τ Γ 4 Γ 5 and uτ ) = uτ + ). On the other hand one can see that u is an increasing function while moving along Γ 1 Γ or Γ 4 Γ 5 in the direction of increasing imaginary part. Therefore u has two inverses in R + : τ + and τ. Before writing them down let us notice that the point C ε defined in 1) can be written as C ε = i + ε e iθε with θ ε, π/) 1) and has the following expression in the coordinates u defined in 11) with uc ε ) = uc ε ) = ε 1 ε ieiθε = ε k ε 1) k ε = ε 1 = 1 + ε ) + ε ) ieiθε sin θ ε cos θ ε 1, since by construction, θ ε, π/) and then < sin θ ε. Now, we can write the inverses of the function u τ + :[uc ε ), + ) Γ 4 Γ 5 τ :[uc ε ), + ) Γ 1 Γ 14) u ξu) + iηu), u ξu) + iηu). The change 11) is useful over Γ 1 Γ and Γ 4 Γ 5, thus performing this change in 85), we have that for any ε > Nq, m, n) = d m,ne q G [ ] G m+n 1 [F m,nu) + Fm,nu)]e G q u du + i)e q G fm,nτ)dτ q 15) uc ε) Γ ) ) 1/ 1/ d m,n = i m+n n m 16) F m,nu) ± 1 = τ ± u) i) m+1 τ ± u) + i) n+1 17) f q m,nτ) = e q G hτ) τ i) m, 18) τ + i) n hτ) is given in 96). Now several helpful lemmas follow. Lemma 1. Let m, n Z, m, n and d m,n be defined by equation 16). Then d m,n e 1/ m+n if m + n >

24 Proof. Let s N, then ) 1/ = 1) s s s! 1 s 1 s 1 1 ) 1 ) ) s 1 = 1 [ s 1 ) s= 1 s) 1 s lim 1 1 s s s 1 s ) s= e 1/ ] ] Next lemma gives information about the functions F ± m,nu). Lemma 14. The function F ± m,nu) defined in 17) has the expansion the coefficients d m,n j F ± m,nu) = ± u) m 1 satisfy d m,n j= = 1/i) n+1, d m,n j d m,n j ± u) j 19) 4 Consequently, the series 19) is convergent for u < /. Proof. Let us introduce the function T ± m,nx) := ±)x m+1 F ± m,nx ) = ) m ) j+. 11) j= d m,n j ±x) j, which is well defined since u = x is a good change of variables in R + and has the two inverses x = ± u. To bound the coefficients d m,n j we use Cauchy formula: ±1) j d m,n j = 1 πi x =ε Applying the change of variables we obtain x = ± ±1) j d m,n j = 1 πi = 1 T ± m,nx) x j+1 τ i) i τ i) = ±τ i) τ i =ρ i j+1 m j+1 πm τ i) m j m j τ i =ρ dx = 1 πi x =ε iτ) m j F m,nx ± ) dx. x j m 1 i i τ i)) = ±τ iτ), 111) 1 1 iτ) τ i) m+1 τ + i) n+1 dτ iτ) 1 dτ τ i) j+1 τ + i) n τ + i) j+1 m Now, taking ρ = 1 and using that τ + i 1 and that τ + i 4 we have d m,n j 4 ) m ) j+, which is the desired bound. From this bound it is clear that the series defining T ± m,nx) is convergent for x < / and therefore the one for F ± m,nu) is convergent for u < /.. 4

25 From equation 19) we have F ± m,nu) = ± u) m 1 m j= the regular part of the function F ± m,nu) is given by g ± m,n± u) = ± u) m 1 d m,n j ± u) j + g ± m,n± u), 11) j=m+1 and d m,n j are defined by equation 19) and satisfy bounds 11). d m,n j ± u) j 11) Lemma 15. Let g ± m,n± u) as in equation 11), < β < 1 and < u < β /. Then Proof. It is clear from equation 11) that g ± m,n ± u) < 9 1 β m. g ± m,n± u) = s= d m,n s+m+1 ± u) s. Since by hypothesis < u < β / with β < 1, we can apply lemma 14 to get g ± m,n± u) which proves the lemma. ) m ) m+4 4 ) m ) m+4 4 s= s= Next proposition gives a bound for Nq, m, n). ) s u) s ) s β ) s 9 / = 1 β m Proposition 16. Let Nq, m, n) as defined in 15) for q >, m, n, m + n >, G > 1. Then Nq, m, n) n+m+ e q G m n 1/ e qg /. Proof. We will bound the integrals of Nq, m, n) choosing e q G G ε = G /, G > 1. We write down then, using 1) and that k ε > 1, F ± G q m,nu)e u du uc ε) F m,nu) e ± q G u du G k ɛ F m,nuc ± ε )) u G m+ [ du G )) n+1 qg e q F m,nu) e ± q G u du G ] G m. 114) It remains only the last integral of 15) the integrand is given in 18) and the domain Γ in 99). The path Γ can be parameterized by τ = i + G e iθ with θ [θ 1, θ ] = [π θ ε, θ ε ], 115) with θ ε given in 1). If we define τθ) hθ) = hτθ)) = i 5 ) + τθ),

26 a straightforward computation using 11) shows that and then, as G > 1: hθ) = G e iθ + 1 i G e iθ ) e q G hθ) = e q G e q cos θ+ 1 G sin θ) e q G e q 1+ 1 G ) e q G e q. 116) Note that, by 115), over Γ we have that τ i = G < 1 and therefore τ + i > 1, and we can bound the last integral of 15) using 116): e Γ q G hτ) θ τ i) m dτ τ + i) n = e q G hθ) τθ) i) m τθ) + i) n i G e iθ dθ θ θ 1 e q G hθ) θ G / m G dθ θ 1 θ 1 e q G e q G m G dθ πg m / e q G e q. 117) From lemma 1 and the bounds 114) and 117), we can bound Nq, m, n) by equation 15) as follows Nq, m, n) e 1/ m+n G q e G m n πe q) m+n+ e q G q e G m n 1. Next proposition provides an asymptotic expression for Nq, m, n). Proposition 17. Let n + m > and the constants d m,n j be defined by equation 19) and d n,m by equation 16). For q >, m, n and G > 1 we have Nq, m, n) = d m,ne q G [ m ] G m+n 1 1) s π q s 1 s 1)!! dm,n m s Gs + T q m,n + Rm,n q s= T q m,n 45 m+ G R q m,n 18 q m 1 G m. When s = the factor 1/s 1)!! in the formula above should be replaced by 1. Proof. We proceed as in the proof of proposition 16 changing the path of integration to the path Γ defined in 98) leading to the integral 15). The important fact is that the integral 15) does not depend on ε. So, we will compute only the ε-independent terms of that integral. We will follow a series of lemmas leading to the proof of the statement. Lemma 18. Let < ε < 1 and uc ε ) be as in equation 1), F ± m,n defined by 17). For any ε > small enough we have, if G > 1: uc ε) the constants d m,n j F m,nu)e ± G q u du = m j= uc ε) G q e u d m,n j are defined by equation 19) and Ê satisfies ± u) m 1+j du + Ê Ê 45 m+ G. Proof. Let us take u = β / with β = /.79. A simple calculation using 11) shows that τ ± u ) i = 1. By definition, for ε > small enough we have that < uc ε ) < u < u < /, so uc ε) F m,nu)e ± G q u du = u uc ε) F m,nu)e ± G q u du + Ê1 6