arxiv: v1 [nlin.cd] 25 Jul 2013

Transcription

1 Evoltion of the tangent vectors and localization of the stable and nstable manifolds of hyperbolic orbits by Fast Lyapnov Indicators Massimiliano Gzzo Dipartimento di Matematica Via Trieste, Padova, Italy arxiv: v [nlin.cd] 25 Jl 203 Elena Lega Université de Nice Sophia Antipolis, CNRS UMR 7293 Observatoire de la Côte d Azr Bv. de l Observatoire, B.P. 4229, Nice cedex 4, France elena.lega@oca.e Abstract The Fast Lyapnov Indicators are fnctions defined on the tangent fiber of the phase space of a discrete or continos dynamical system, by sing a finite nmber of iterations of the dynamics. In the last decade, they have been largely sed in nmerical comptations to localize the resonances in the phase space and, more recently, also the stable and nstable manifolds of normally hyperbolic invariant manifolds. In this paper, we provide an analytic description of the growth of tangent vectors for orbits with initial conditions which are close to the stable-nstable manifolds of a hyperbolic saddle point of an area preserving map. The representation explains why the Fast Lyapnov Indicator detects the stable-nstable manifolds of all fixed points which satisfy a certain condition. If the condition is not satisfied, a sitably modified Fast Lyapnov Indicator can be still sed to detect the stable-nstable manifolds. The new method allows for a detection of the manifolds with a nmber of precision digits which increases linearly with respect to the integration time. We illstrate the method on the critical problem of detection of the so called tbe manifolds of the Lyapnov orbits of L, L 2 in the circlar restricted three body problem. Introdction Since the first detection of chaotic motions in 964 Henon Heiles [7], several indicators have been largely sed to characterize the different dynamics of dynamical systems. Many dynamical indicators, sch as the Lyapnov characteristic exponents and the more recently introdced finite time chaos indicators sch as the Finite Time Lyapnov Exponent FTLE [3], Fast Lyapnov Indicator FLI [7], Mean Exponential Growth of Nearby Orbits MEGNO [4], are defined by the local divergence of nearby initial conditions, that is by the variational dynamics. For example, for a discrete dynamical system defined by the map Φ : M M z Φz, with M R n open invariant, by denoting with DΦ z the tangent map of Φ at z: DΦ z : R n R n

2 v DΦ z v, 2 the characteristic Lyapnov exponent of a point z M and a vector v R n \0 is defined by the limit λz, v = lim T + T log DΦ T z v, 3 v and the largest Lyapnov exponent of z is the maximm of λz, v for v 0. As a matter of fact, the nmerical estimation of the characteristic Lyapnov exponents see [2] relies on extrapolation of finite time comptations, since compters cannot integrate on infinite time intervals. The so called finite time chaos indicators sch as the FTLE, the FLI and the MEGNO have been afterwards introdced as srrogate indicators of the largest Lyapnov exponent, with the aim to discriminate between reglar orbits and chaotic orbits sing time intervals which are significantly smaller than the time interval reqired for a reliable estimation of the largest characteristic Lyapnov exponent [7], [4]. For example, the fnction Fast Lyapnov Indicator of z and v is simply defined by DΦ T l T z, v = log z v, 4 v and depends parametrically on the integer T > 0, as well as on the choice of a norm on R n. The definition of finite time chaos indicators was jstified by the possibility of their systematic nmerical comptation over large grids of initial conditions in the phase space in a reasonable comptational time. We remar that, specifically in Celestial Mechanics, the nmerical detection of the resonances of a system sing dynamical indicators, both formlated sing the Lyapnov exponent theory or alternatively the Forier analysis sch as the freqency analysis [9, 2, 20], is one of the major tools for stdying its long term instability for recent examples, see [27, 28, 26, 25, 8, 9, 33]. The papers [5],[], focsed and proved properties of the finite time chaos indicators, specifically the FLI, which are lost by taing the limit of l T z, v/t, ths differentiating the se of these indicators from the parent largest Lyapnov characteristic exponent. Specifically, since [5],[], the FLI has been sed to discriminate reglar motions of different natre: for example the motions which are reglar becase are spported by a KAM tors from the reglar motions in the resonances of a system. This property of the FLI improved a lot the precision in the nmerical localization of different types of resonant motions, the so called Arnold web, and provided the technical tool for the first nmerical comptations of diffsion along the resonances of qasi integrable systems in exponentially long times [22, 2, 6, 4, 6], as depicted in the celebrate Arnold s paper []. More recently, the FLI has been sccessflly sed to compte the stable and nstable manifolds of normally hyperbolic invariant manifolds of the standard map and its generalizations [0, 3], and of the three body problem [32, 23, 5]. In these cases it happens that, depending on the choice of the parameter T, finite pieces of the stable and nstable manifolds appear as sharp local maxima of the FLI. As a matter of fact, the possibility of sharp detection of the stable and nstable manifolds of a fixed point, or periodic orbit, with a FLI comptation is not general and trns ot to be a property of the manifolds. A model example is represented by the stable and nstable manifold of the fixed point 0, 0 of the symplectic map sinϕ + I Φϕ, I = ϕ + I, I + σ cosϕ + I + 2 2, 5 where ϕ, I M = 2πS R are the phase space variables, σ = ± is a parameter: for σ = the FLI may be sed for excellent detection of the manifolds; for σ = the FLI does not provide any detection. To explain this fact, in this paper we provide a representation for the growth of tangent vectors for orbits with initial conditions close to the stable manifold of a saddle fixed point. To better illstrate the theory, we consider a two dimensional area preserving map with a saddle fixed point z, bt the techniqes which we se the local stable manifold theorem and Lipschitz estimates can be sed also in the higher dimensional cases. The two dimensional case allows s to treat also Poincaré sections of the circlar restricted three body problem. 2

3 Let s denote by z the saddle point of the map, and by W s, W its stable and nstable manifold. We consider a point z s W s, a tangent vector v R 2, and we provide estimates abot the norm of the tangent vector DΦ T z v, for points z / W s which are close to z s. As it is sal, the same argments applied to the inverse map Φ, allow to reformlate the reslt by exchanging the role of the stable manifold with that of the nstable manifold. For the points z which are the sitably close to z s W s, the orbit Φ z follows closely the orbit Φ z s for any T, and DΦ z v remains close to DΦ zs v as well. The most interesting sitation happens for the points z which are little more distant from the stable manifold: their orbit i follows closely the orbit Φ z s only for smaller than some K 0 < T ; ii then remains close to the hyperbolic fixed point for a nmber of iterations which increases logarithmically with respect to some distance between z and z s, see Section 2, iii then follows closely the orbit of a point on the nstable manifold W in the remaining iterations. It is dring the process iii that the growth of the tangent vector DΦ z v can be significantly different from the growth of DΦ zs v, and the difference may be possibly sed to characterize the distance of z from the stable manifold. As a matter of fact, with evidence any difference may exist only de to the non linearity of the map Φ. In Section 2 we provide a representation for sch a difference, and we discss a condition which garantees the desired scaling of the FLI with respect to the distance of z from the stable manifold. If this condition is satisfied, the comptation of the FLI on a grid of initial conditions provides a sharp detection of the stable and nstable manifolds see Section 3: typically, the time T sed for the FLI comptation, which is the time needed by the orbits with initial condition z to approach the fixed point z, trns ot to be proportional to the nmber of precision digits of the detection. At the light of the representation provided in Section 2, we propose a generalization of the FLI which weaens a lot the condition for the detection of the stable and nstable manifold. For any smooth and positive fnction : R 2 R + we define the modified FLI indicator of z M, v R 2 at time T > 0, as the T th element of the seqence DΦzj v j l = ln v, l j+ = l j + z j ln, 6 v j where z j := Φ j z and v j := DΦ j z j v. The traditional FLI is obtained with the choice z = for any z M. We consider the alternative case of fnctions z which are test fnctions of some neighborhood B M of the fixed point, and precisely with z = for z B, and z = 0 for z otside a given open set V B. When the diameter of the set B is small, bt not necessarily extremely small, the comptation of the modified FLI indicator allows to refine the localization of the fixed point by many orders of magnitde. Therefore, at variance with the traditional FLI indicator, the modified indicators are proposed as a general tool for the nmerical detection of the stable and nstable manifolds. An illstration of the potentialities of these indicators is given in Section 3, where we provide comptations of the stable and nstable manifolds and their heteroclinic intersections, of the Lyapnov orbits arond L, L 2 of the circlar restricted three body problem. The application is particlarly critical, since these manifolds are located in a region of the phase space close to the singlarity de to the secondary mass. The paper is strctred as follows. In Section 2 we provide the representation for the evoltion DΦ T z v of the norm of tangent vector v for points z / W s which are sitably close to the stable manifold, and we also discss a sfficient condition for the FLI to detect sharply the stable and nstable manifolds of the map. In Section 3 we provide an illstration of the method for the comptation of the stable and nstable manifolds of the Lyapnov orbits arond L, L 2 of the circlar restricted three body problem; in Section 4 we provide the proof of Proposition. In Section 5 we formlate and prove two technical lemmas. Finally, Conclsions are provided in Section 6. 3

4 2 Evoltion of the tangent vectors close to the stable manifolds of the saddle points of two dimensional area preserving maps We consider a smooth two dimensional area preserving map: Φz = Az + fz, 7 where A is a 2 2 diagonal matrix with A = λ >, A 22 = /λ and f is at least qadratic in z, z 2, that is f i 0, 0 = 0 and fi z j 0, 0 = 0, for any i, j =, 2. Therefore, the origin is a saddle fixed point. We need to introdce some constants which characterize the analytic properties of Φ. We denote by λ Φ, λ Φ, λ DΦ the Lipschitz constants of Φ, Φ, DΦ respectively defined with respect to the norm := max{, 2 }, in the set BR = {z : z R}. Also, we set η sch that, for any z BR, we have fz η z 2, Df z η z, D 2 f z η fz fz η max{ z, z } z z, where D 2 f z denotes the Hessian matrix of f at the point z and, by denoting with Φ z = A z + fz the inverse map, we also have fz η z 2, D f D z η z, 2 fz η Moreover, since Φ is a diffeomorphism, we have fz fz η max{ z, z } z z. σ = min min DΦ zv > 0. 8 z BR v = By the local stable manifold theorem, we consider e neighborhood Br of the origin where the local stable and nstable manifolds W l s, W l are Cartesian graphs over the z 2 and z axes respectively, that is W l s = {z : z 2 r, z = w s z 2 } W l = {z : z r, z 2 = w z } with w s 0 = w 0 = 0, w s0 = w 0 = 0 and, by possibly increasing η, and w s z 2 η z 2 2, w z η z 2 w s ξ w s ξ λ w max{ ξ, ξ } ξ ξ w ξ w ξ λ w max{ ξ, ξ } ξ ξ. We denote by W s, W the stable and nstable manifolds of the origin. We consider a point z s W s, a tangent vector v R 2, and we provide estimates abot the norm of the tangent vector DΦ T z v, for points z / W s which are sitably close to z s, precisely in a crve, with z 0 = z s and z = ε. Let s consider a small δ := δ0 T, with δ 0 satisfying δ 0 min 6 max, η 2 e 3 λ 2, r. λ 2 Then, we consider the minimm T s := T s δ sch that Φ Ts z s Bδ 2δ 2. Typically, one has T s ln/δ. For all ε, we have see Lemma 5.2: Φ T s Φ Ts z s λ T s Φ ε 9 4

5 W z s s.. z ε π ε.. Φ Ts zε Φ. Ts+T ε ζ ε W Figre : Illstration of z s, ; of Φ Ts and its parallel projection π ε on the local stable manifold; of Φ Ts+Tε and its parallel projection ζ ε on the local nstable manifold. DΦ T s v DΦ Ts z s v DΦ T s z s v λ T s ε, 0 where λ = maxλ Φ, DΦ + λ DΦ /σ. We consider only the small ε satisfying λ Ts ε < δ 2, so that Φ Ts Bδ δ 2, are close to Φ Ts z s and DΦ T s v are close to DΦ T s z s v. We rename the vector DΦ Ts z s v as follows: w = w s + w = DΦ Ts z s v, where w s, w are the orthogonal projections of w over the stable and nstable spaces of the matrix A, i.e. the z 2 and z axes, respectively. We need a condition which ensres that v is not close to some special contracting direction. Precisely, we assme that the initial vector v is sch that w s w = w. In particlar, for any 0, we have A w = λ w. Let s denote by π ε = w s Φ Ts 2, Φ Ts 2 Ws l the parallel projection of Φ Ts on the local stable manifold see figre, that is the point on Ws l with z 2 = Φ Ts 2, and by ε = Φ Ts w s Φ Ts 2 the distance between Φ Ts and the point π ε. Since ε depends continosly on ε, 0 = 0, and the local stable manifold is invariant, there exists ε sch that ε is strictly monotone increasing fnction of ε [0, ε ]. We have also see Section 4: ε + λ w λ Ts Φ ε, so that if + λ w λ Ts ε < δ 2 we have π ε Bδ. We se ε to parameterize the distance of from the stable manifold W s, and we introdce the time [ T ε = ln eδ ] 2 ln λ ε 5

6 which, as we will prove see Lemma 4., is reqired by the orbit with initial condition Φ Ts to exit from Bδ. We also denote by ζ ε = Φ Ts+Tε, w Φ Ts+Tε W l the parallel projection of Φ Ts+Tε over the local nstable manifold. Proposition Let s consider any large T satisfying with αt Ts eδλ ε 3 σ Ts δ 0 max, η + λ w λ Ts T 2 4 αt Ts eλ By denoting with ε 0 the constant sch that T > T s + α α = ε0 then, for any ε ε 0, if T ε T T s we have ln λ ln λ + ln λ. 5 αt Ts = eδλ, 6 if αt T s T ε < T T s we have DΦ T zε v A T Ts w λ T T s w T, w = DΦ Ts z s v, 7 DΦ T zε v DΦ T z s v + j DΦ ζ ε T λ j, j = T T s T ε. 8 The proof is reported in Section 4. Remar. Conditions 3, 4 and 5 may be all satisfied by times T which are sitably large, bt not necessarily extremely large, becase of the presence of the exponentials in 3 and 4, and becase of the typical dependence T s δ ln/δ ln T. Therefore, the proposition is meaningfl also for ε 0 which are small, bt not necessarily extremely small. Moreover, from the definition of ε 0, apart from a small difference de to the se of the integer part in the definition of T ε, we have T ε0 αt T s, and T T s T ε T := T T s α. For z s W s, and for all the points which are so close to the stable manifold that T ε T T s, the FLI is approximated by ln A T T s w = T Ts ln λ + ln w. Therefore, the only possibility for the FLI to strongly decrease by increasing ε is that, for αt T s T ε < T T s, we have an exponential decrement of DΦ j ζ ε /λ j with respect to j. The assmption which garantees a desired scaling of the FLI with respect to ε is sp ε:αt T s T ε T T s DΦ T Ts Tε ζ ε 9 Cλ T Ts Tε 6

7 with some C <, so that we have ln DΦ T zε v ln DΦ T zs v T Ts T ε ln C + ln +. T From the definition of T ε, we have therefore a linear decrement of the FLI with respect to ln ε, p to the maximm vale of T T s T ε αt T s. Therefore, at the exponentially small distance from the manifold 6 the FLI has decreased of a qantity which is proportional to integration time T, and conversely, the differences of nits in the FLI vale typically determines a proportional nmber of precision digits in the localization of the stable manifold. With evidence, condition 9 may be satisfied if DΦ z has an absolte maximm for z T Φ W l. For example, the condition may be satisfied for the map 5 with σ =, since the origin is a local strict maximm for DΦ z, z W, while it is not satisfied for σ =, since in this case the origin is a local strict minimm for DΦ z, z W. In any case, it is not practical to verify if condition 9 is satisfied by a certain choices of the parameters. Therefore, at the light of the above analysis, we consider a generalization of the FLI indicators which depend on a fnction : R 2 R + as follows: let s consider z M, v T z M, and T > 0. Then, we consider l T z, v defined as the T th element of the seqence DΦzj v j l = ln v, l j+ = l j + z j ln, 20 v j where z j := Φ j z and v j := DΦ j z j v. The sal FLI is obtained by z = for any z M. We consider the alternative case of fnctions z which are test fnctions of some neighborhood B M of the fixed point, and precisely with z = for z B, and z = 0 for z otside a given open set V B. We remar that the set B needs to be small, bt not necessarily extremely small. For example, if B Bδ, we only need, in V \B, DΦ z z Cλ for some C <. The fnction described above depends on a specific hyperbolic fixed point. If one is interested in the stable or nstable manifolds of more fixed points or hyperbolic periodic orbits, with the same nmerical integration of the variational eqations, forward and bacward in time, one may compte the FLI indicators related to the different fixed points withot increasing significantly the comptational time, and se the reslts to find, for example, homoclinic and heteroclinic intersections between the different manifolds. If instead, one is interested in determining with a single nmerical integration the largest nmber of manifolds in some finite domain B, one can divide the domain B in many small sets B j, j N, and compte the N indicators FLI j adapted to the sets B j. This procedre increases the comptational time only logarithmically with N, since the time reqired for the nmerical localization of a point in one of the sets B j increases logarithmically with N. Then, the portrait of all the manifolds is obtained by representing, for any initial condition, the maximm between all the FLI j. Therefore, at variance with the traditional FLI indicator, the modified indicators are proposed as a general tool for the nmerical detection of the stable and nstable manifolds. 3 A nmerical example: the tbe manifolds of L and L 2 in the planar circlar restricted three body problem The circlar restricted three-body problem describes the motion of a massless body P in the gravitation field of two massive bodies P and P 2, called primary and secondary body respectively, which rotate niformly arond their common center of mass. In a rotating frame xoy, the eqations of motion of P are: { ẍ = 2ẏ + x µ x+µ ÿ = r 3 2ẋ + y µ y r 3 µ x +µ r2 3 µ y r

8 Figre 2: Projection on the plane x-y of the Lyapnov orbits related to the points L and L 2, for the vale C = of the Jacobi constant. The shaded area represents a region of the orbit plane which is forbidden for this vale of the Jacobi constant. where the nits of masses, lengths and time have been chosen so that the masses of P and P 2 are µ and µ µ /2 respectively, their coordinates are µ, 0 and µ, 0 and their revoltion period is 2π. We denoted by r 2 = x + µ 2 + y 2 and by r 2 2 = x + µ 2 + y 2. As it is well nown, eqations 2 have an integral of motion, the so called Jacobi constant, defined by: Cx, y, ẋ, ẏ = x 2 + y µ r + 2 µ r 2 ẋ 2 ẏ 2, 22 and five eqilibria sally denoted by L,..., L 5. Here we consider µ = , which corresponds to the Jpiter Sn mass ratio vale, and a vale of the Jacobi constant slightly smaller than Cx L2, 0, 0, 0 := C 2. As it is extensively explained in [8], in these conditions, one may find particlarly interesting dynamics, which we briefly smmarize. The eqilibrim points L, L 2 are partially hyperbolic, and their center manifolds WL c, WL c 2 are two dimensional, and foliated near L, L 2 respectively by periodic orbits called Lyapnov orbits. For vales C of the Jacobi constant slightly smaller than C 2, there exist one Lyapnov orbit related to L and one Lyapnov orbit related to L 2 respectively with Jacobi constant eqal to C see figre 2. The Lyapnov orbits are hyperbolic, and transverse intersections of their stable and nstable manifolds sally called tbe manifolds prodce the complicate dynamics related to the heteroclinic chaos. The nmerical comptation of the tbe manifolds has been afforded in several papers, and has important implications also for modern space mission design see [29], [8]. In this Section we analyze the FLI method for the detection of the tbe manifolds introdced in [24, 5] at the light of the theoretical analysis performed in Section 2, and we show that the method allows for a detection of the manifolds with a nmber of precision digits which increase linearly with respect to the integration time. Moreover, the modified FLI allows s to compte the manifolds with a precision limited only by the rond off of the nmerical comptations. We report here three nmerical experiments. In the first one we illstrate the nmerical precision of the FLI method in the determination of the stable tbe manifold of a Lyapnov periodic orbit arond L ; in the second one, we provide some snapshots of the stable tbe manifold of the Lyapnov periodic orbit arond L 2 and the nstable tbe manifold of the Lyapnov periodic orbit arond L, obtained by extending the integration time; in the third one we illstrate the nmerical precision of the FLI method for the localization of a heteroclinic intersection between these two manifolds. We remar that these comptations are particlarly critical since the tbe manifolds are located in a region of the phase space 8

9 close to the singlarity at x, y = µ, 0. In these circmstances, the nmerical comptation of both eqations of motions 2 and their variational eqations becomes critical, and several approaches have been introdced see [32, 23, 3, 5]. For the comptation of the tbe manifolds, we find particlarly sefl to define the variational eqation in the space of the variables obtained by reglarizing eqations 2 with respect to the secondary mass, as in [3, 5]. Precisely, we consider the Levi Civita reglarization defined by the space transformation { x µ = y = and by the fictitios time s related to t by dt = r 2 ds. The eqations of motion in the variables, 2, and fictitios time s are see for example [30]: { = 4 [a + b + c 2 ] 2 = 4 [a b c ] with: a = 2 µ r C + x 2 + y 2 b = 4y + 2r 2 x 2 µr2x +µ r 3 c = 2r 2 y 4x 2 µr2y r 3 where C denotes the vale of the Jacobi constant, and the primed derivatives denote derivatives with respect the fictitios time s. To define the FLI, we first write 24 as a system of first order differential eqations: = v 2 = v 2 v = 4 [a + b 26 + c 2 ] v 2 = 4 [a b 2 + c ] and we introdce its compact form: with ξ =, 2, v, v 2. The variational eqations of 27 are therefore: { 25 ξ = F ξ 27 ξ = F ξ w = F ξ ξw, 28 where w R 4 represents a tangent vector. Following [5], we here consider the reglarized FLI indicator defined by F LIξ0, w0, T = log wst 29 where ξs, ws denotes the soltion of the variational eqations 28 with initial condition ξ0, w0 and st is the fictitios time which corresponds to the physical time T for that orbit. The indicator 29 will be compted also for negative times T < 0. FLI detection of the tbe manifolds. In order to test the precision of the FLI method in the localization of the tbe manifolds, we consider a point z s = x s, y s, ẋ s, ẏ s WL s in the stable tbe manifold of the Lyapnov orbit arond L see Figre 3, and we compte the traditional and modified FLIs for a set of many initial conditions. with x0, y0 = x s, ẋ s see Fig.3, log y0 y s in the interval [ 25, ] and y0 obtained from the vale of the Jacobi constant C = The integration times are respectively T = 5 and T = 25. We appreciate a localization of the manifold determined by a linear decrement of the FLI with respect to log y0 y s. The time T = 5 allows s to localize the manifold with a precision of order 0 5, which is greatly improved by sing T = 25. We obtain a good localization of the manifold already with the traditional FLI, see Figre 4, althogh the 9

10 Figre 3: Projection on the plane x, y of an orbit with initial condition z s = x s, y s, ẋ s, ẏ s W s L, with x s = , y s = , ẋ s = , and ẏ s obtained from the Jacobi constant C = The shaded area represents a region of the orbit plane which is forbidden for the vale C of the Jacobi constant. irreglarities in the FLI crve limit the precision of the localization to 0 22, higher than the nmerical rond off precision. Then, we considered a modified FLI defined by eqations 6 with fnction z which is a test fnction of a neighborhood of the Lyapnov orbit γ arond L. Precisely, we se a test fnction defined by: z = 2 if z γ r 2 z γ [cos r 2 π + ] if r 2 < z γ 3r 2 0 if z γ > 3r 2 where z γ denotes the distance between z and the Lyapnov orbit γ we set r = 0 3 in the following comptations. Also in this case the time T = 5 allows s to localize the manifold with a precision of order 0 5, while the time T = 25 allows s to localize the manifold more precisely than The se of the modified FLI has eliminated the irreglarities in the crves of Figre 4, and improved the precision of the localization. As a matter of fact, the precision of the localization is redced to the rond off sed for the nmerical comptation. Snapshots of tbe manifolds of WL and WL s 2. Motivated by these reslts, we obtained sharp representations of the intersections WL s 2 Σ, WL Σ of the stable tbe manifold WL s 2 of the Lyapnov orbit γ 2 arond L 2 and of the nstable tbe manifold WL of the Lyapnov orbit γ arond L with the two dimensional section of the phase space defined by Σ = {x, y, ẋ, ẏ : y = 0, ẏ 0 : Cx, 0, ẋ, ẏ = C}. 3 Any point z Σ is parameterized and identified by its two components x, ẋ. The representation of the manifolds are obtained by compting the modified FLIs on refined grids of initial conditions x, ẋ on Σ for different integration times T. The stable manifold WL s 2 is obtained by compting the modified FLI 30 0

11 30 25 t=5 t=25 20 FLI logy0-y s 0 Figre 4: Vales of the traditional FLI compted on a set of 960 initial conditions with x0, y0 = x s, ẋ s see Fig.3, log y0 y s in the interval [ 25, ] and ẏ0 obtained from the Jacobi constant C = The integration times are respectively T = 5 and T = 25, the negative vales correspond to initial conditions with y0 < y s. We appreciate a localization of the manifold determined by a linear decrement of the FLI with respect to log y0 y s. The time T = 5 allows s to localize the manifold with a precision of order 0 5, while the time T = 25 allows s to localize the manifold more precisely than t=5 t=25 8 FLI logy0-y s 0 Figre 5: Vales of the modified FLI defined by eqations 6 with fnction z which is a test fnction of a neighborhood of the Lyapnov orbit arond L. The initial conditions are the same 960 initial conditions considered in Figre 4, that is x0, y0 = x s, ẋ s see Fig.3, log y0 y s in the interval [ 25, ] and y0 obtained from the Jacobi constant C = The integration times are respectively T = 5 and T = 25, the negative vales correspond to initial conditions with y0 < y s. We appreciate a localization of the manifold determined by a linear decrement of the FLI with respect to log y0 y s. The time T = 5 allows s to localize the manifold with a precision of order 0 5, while the time T = 25 allows s to localize the manifold more precisely than The se of the modified FLI has eliminated the irreglarities in the crves of Figre 4, and improved he precision of the localization. As a matter of fact, the precision of the localization is redced to the rond off sed for the nmerical comptation.

12 Figre 6: Representation of the modified FLIs compted on a grid of initial conditions reglarly spaced on x, ẋ the axes on the pictre the other initial conditions are y = 0 and ẏ is compted from the Jacobi constant C = , compted with integration time T = 5. In order to represent both manifolds on the same pictre, we represent with a color scale the weighted average 34 of the two indicators FLI, FLI 2 with weight w = 00. The yellow crves on the pictre correspond to different lobes of the manifolds. on a time T 2, sing a test fnction defined by if z γ 2 r2 2 z γ2 z = 2 [cos r 2 2 π + ] if r 22 < z γ 2 3r2 2 0 if z γ 2 > 3r2 2 where z γ 2 denotes the distance between z and the Lyapnov orbit γ 2 and r 2 = The nstable manifold WL is obtained by compting the modified FLI on a negative time T, sing a test fnction defined by if z γ r 2 z γ z = 2 [cos r 2 π + ] if r 2 < z γ 3r if z γ > 3r 2 where z γ denotes the distance between z and the Lyapnov orbit γ and r = 0 3. In sch a way, for any x, ẋ we compte the modified FLIs: FLI, FLI 2. The representation of both manifolds on the same pictre is obtained by representing with a color scale a weighted average of the two indicators: 32 wfli + FLI w + The reslts are represented in Figres 6 and 7 for T = 5 and T = 00 respectively. We clearly appreciate different lobes of both manifolds already for the shorter integration time T = 5. The longer time T = 00 allows s to appreciate additional lobes, which contain initial condition approaching the manifolds only after several revoltion periods of Jpiter. Localization of heteroclinic intersections. The detection of both manifolds W L and W s L 2 on the same pictre see Figre 6 and Figre 7 allows s to obtain a precise localization of the heteroclinic 2

13 Figre 7: Representation of the modified FLIs compted on a grid of initial conditions reglarly spaced on x, x the axes on the pictre the other initial conditions are y = 0 and y is compted from the Jacobi constant C = , compted with an integration time T = 00. In order to represent both manifolds on the same pictre, we represent with a color scale the weighted average 34 of the two indicators FLI, FLI2 with weight w = 500. The yellow crves on the pictre correspond to different lobes of the manifolds. De to the integration time which is mch longer than the time sed in Figre 6, many additional lobes of the tbe manifolds of both γ and γ2 appear on this figre. Their corresponding initial conditions approach the manifolds only after several revoltion periods of Jpiter. 3

14 Figre 8: Comptation of the averaged FLI 34 on a grid of initial conditions centered on the point z he,5 of coordinates: x he,5 = , y he,5 = 0, ẋ he,5 = The velocity ẏ he,5 is obtained from the Jacobi constant C = The integration time is T = 8. The vales of the FLI are provided as the average between FLI and FLI 2. A sharp detection of both manifolds appears thans to the differentiation of the FLI vales on this refined grid. The maximm vale of the FLI in this pictre provides a new refined initial condition for the orbit plotted in Fig.9. intersections points, which we denote by z he. Precisely, the intersection between the two yellow crves in the box of Fig.6 corresponds to an intersection point z he. Of corse, accordingly to the resoltion of the comptation, at first we are only able to determine a point z he, in the box which is close z he. To improve the localization of z he we compte again the modified FLIs on a refined grid of points in the box of Fig.6, and we obtain a new point z he,2 the point with the maximm vale of the averaged FLI 34 closer to the intersection point. The procedre is iterated by compting again the FLIs on zoomed ot grids of initial conditions centered on z he,j with j = 2,...5, with increasing integration times to increase the nmber of precision digits in the localization of the heteroclinic point. In Fig.8 we plot the FLI vales compted on a grid of initial conditions centered on the point z he,5, sing the integration time T = 8. The maximm vale of the FLI in this pictre provides a new refined initial condition that we sed to compte the heteroclinic orbit shown in Fig.9. The convergence of the forward bacward integration towards the Lyapnov orbit related to L 2 L clearly shows the validity of the method for the precise localization of heteroclinic orbits. 4 Proofs Proof of Proposition. We first remar that 3 implies ε 0 ε, and condition 4 implies ε 0 δ 2 0 max, η + λ w λ Ts T 2. The proof of Proposition is a conseqence of the following: Lemma 4. For any ε, δ satisfying δ min max, η + λ w λ Ts ε δ , r 2, η, 4e 2 λ η λ, 2e 3 λ 2 η T ε 36 4

15 Figre 9: Projection on the plane x, y of the heteroclinic orbit fond throgh the maximm of the FLI see text. The conditions are : x he 0 = , y he 0 = 0, ẋ he 0 = and ẏ he 0 obtained from the Jacobi constant C = Ble points: forward integration, the orbit converges to the Lyapnov orbit related to L 2. Red points: bacward integration, the orbit converges to the Lyapnov orbit related to L. we have: δ η ε and the tangent vector DΦ Ts+Tε v satisfies DΦ T s+t ε Φ Ts+Tε e 2 λ δ + η ε 37 Φ Ts+Tε 2 e 2 λ δ + 2η ε 38 Φ Ts+Tε 2 ζ ε 2 4e 3 λ δ, 39 v A Tε w + w s λ T ε δ λ Tε η maxe, η 4e2 λ 2 λ + w. 40 First, as we anticipated in Section 2, the time T ε can be identified as the time reqired by the orbit with initial condition Φ Ts to exit from Bδ and to arrive at the small distance 4e 3 λ δ/λ Tε from the local nstable manifold. If T ε T T s, we can repeat the proof of Lemma 4. by limiting all the estimates to the time interval [0, T ], and obtaining DΦ T zε v A T Ts w λ T T s δ η maxe, η 4e2 λ 2 λ + w, 4 so that 7 is proved. If T ε < T T s, we need an estimate of the growth of the tangent vectors in the remaining time interval [T s + T ε, T ], and we obtain it by comparison with the growth of the tangent vectors of the orbits with initial condition in the point ζ ε on the nstable manifold. We will provide estimates of the FLI for T ε in the interval: ln λ T T s δ T ε < T T s δ. 42 ln λ + ln λ 5

16 ln λ Let s consider j = T T s T ε {, T T s δ ln λ+ln λ }. First, we have In fact, since DΦ T v = DΦ j Φ DΦ T v DΦ j ζ ε DΦ Ts+Tε v 4e 3 λ j λ δ DΦ j ζ ε DΦ Ts+Tε v, 43 Ts+Tε zεdφts+tε sing Lemmas 4. and 5.2 we obtain DΦ T v DΦ j ζ ε DΦ Ts+Tε v λ j Φ T s+t ε ζ ε DΦ j Therefore, we have λ Tε v = DΦ j ζ ε DΦz Ts+Tε ε v + DΦ jφ Ts+Tε zε DΦjζε DΦ Ts+Tε v, ζ ε DΦ Ts+Tε DΦ T v DΦ j ζ ε DΦ Ts+Tε v + 4e 3 λ j λ δ. λ Tε v 4e 3 λ j λ δ DΦ j ζ ε DΦ Ts+Tε v. We now analyze and compare the FLI for initial conditions at different distances from the stable manifold. We have DΦ T zε v DΦ T z s v DΦ j ζ ε DΦ Ts+Tε v DΦ T z s v + 4e 3 λ j λ δ. λ Tε Using ineqalities 4 and 40, we obtain DΦ T zε v DΦ T zs v DΦ j ζ ε + δ η maxe, η 4e2 λ 2 λ + λ j + 4e 3 λ j λ δ. 44 δ η maxe, η 4e2 λ 2 λ + λ Tε In fact, from 4 and A T Ts w + w s = λ T Ts w, we obtain that for all ε with T ε T T s, inclding z 0 = z s, we have DΦ T v λ T Ts λ Tε δ η maxe, η 4e2 λ 2 λ + w. 45 From 40, for all ε with T ε T T s, we have: DΦ j ζ ε DΦ Ts+Tε v DΦ j DΦ T s+t ε ζ ε v Since for all ε with we have sing also DΦ j A T ε ζ ε w + λ T ε δ η maxe, η 4e2 λ 2 λ + w DΦ j ζ ε λ T ε w + δ η maxe, η 4e2 λ 2 λ ln λ T ε T T s, ln λ + ln λ λ j λ Tε = λt Ts Tε λ Tε, δ η maxe, η 4e2 λ 2 λ + 2eT, 4e3 λ δ < 4T 6

17 we have + δ η maxe, η 4e2 λ 2 λ + δ η maxe, η 4e2 λ 2 λ + so that, from 44, we immediately obtain e 3 λ λ j Proof of. We have: ε = Φ Ts w s Φ Ts 2 z ε Φ Ts Φ Ts z s + δ λ Tε λ Ts Φ ε + λ wλ Ts Φ ε + λ wλ Ts Φ ε. +, T w s Φ Ts 2 z s w s Φ Ts 2 Proof of Lemma 4.. We consider the segments which join Φ π ε and Φ Φ Ts, and define = Φ Φ Ts Φ π ε 2 = Φ 2 Φ Ts Φ 2π ε. We prove that, for all the sch that Φ π ε, Φ Φ Ts BAδ, for A >, we have 2 <. In fact, we have 0 = ε, 0 2 = 0; then, if 2 <, we have 2 = Φ 2 Φ Ts Φ 2π ε = Φ2 Φ Φ Ts Φ 2 Φ π ε as soon as Therefore, we have λ 2 + f2 Φ Φ Ts f 2 Φ π ε λ 2 + Aηδ + 2 λ + 2Aηδ <. + 2Aηδ < λ and as soon as T and We obtained = Φ Φ Ts Φ π ε λ + Aηδ + 2 λ + 2ηAδ λ + 2ηAδ 0 = λ + 2ηAδ ε eλ λ ε λ Aη + 2 λ 2Aηδ λ 2ηAδ ε λ e λ ε 2ηAδ λ 2. e λ ε eλ ε. We now provide an estimate of Φ π ε and Φ 2π ε. We consider the segment which joins the origin 0, 0 and Φ π ε and define δ = Φ π ε, δ 2 = Φ 2 π ε. 7

18 We have δ < δ 2 for any. In fact, since π ε Bδ and π ε W l s, then Φ π ε Bδ for any and we have δ = Φ π ε = ws Φ 2π ε η Φ 2 π ε 2 ηδδ 2 < δ 2 as soon as For = 0 we have δ 0 = π ε = w s π ε 2 η Then, we have as soon as and Therefore, from we have Φ Ts ηδ <. 2 z s 2 ηδ 2 z s ηδ 2, δ2 0 = Φ Ts 2 z s δ δ 2 = Φ 2 π ε = Φ2 Φ π ε λ Φ 2 π ε + η Φ π ε 2 λ δ 2 + ηδ 2 2 Finally, from 2 < we have From the definition of T ε, we have and therefore we have Φ Ts δ + ηδ δ2 0 + ηδ λ λ 2 δ λ + ηλ δ λ eδ Therefore, since ε + λ w λ Ts Φ ε, as soon as ηλ δ e, δ = w s Φ 2π ε η Φ 2π ε 2 = ηδ2 2 ηe λ δ. e λ ε Φ Φ Ts Φ π ε eλ ε e λ ε ηe λ δ Φ Φ Ts eλ ε + ηe λ δ. Φ 2 Φ Ts + Φ 2 π ε eλ ε + λ eδ. eδ λ Tε < λ eδ, ε ε δ η ε δ ηe δ Φ λ Tε Φ Ts e 2 λ δ + η ε Φ 2 Φ Ts e 2 λ δ + 2η ε. max, η + λ w λ Ts ε < δ 2 we have Φ Φ Ts e 2 λ δ + 2η ε e 2 λ δ + 2η + λ w λ Ts Φ ε < e2 λ δ + δ 2 < 2e 2 λ δ = Aδ for A = 2e 2 λ. The thresholds conditions on δ become δ 4e 2 λ η λ, δ 8e 2 η 8 T ε, δ 2e 3 λ 2 η T ε.

19 We now consider the point ζ ε = Φ Ts+Tε, w Φ Ts+Tε and the segments which join Φ ζ ε and Φ Φ Ts+Tε, for T ε. We already now that Φ Φ Ts+Tε BAδ, and ζ ε = Φ Ts+Tε, w Φ Ts+Tε maxaδ, ηaδ 2 Aδ as soon as ηδ. By definition of local nstable manifold, we have Φ ζ ε BAδ. We define = Φ +Ts+Tε Φ ζ ε in particlar we have 2 = Φ +Ts+Tε 2 Φ 2 ζ ε, 0 = 0, 0 2 := ε. By repeating the above argments sing the inverse map Φ x, we have < 2 for any and: ε = e λ ε 2 eλ ε Φ Ts+Tε 2 ζ ε 2 e It remains to prove 40. For any T ε, we have: and Φ Φ Ts eλ ε + max, ηe λ T ε =0 Φ Φ Ts T ε =0 T ε 2eλ maxe, ηδ λ =0 λ Tε δ e 2 e 2 λ λ Tε Tε 2 e λ λ Tε λ Tε 2Aδ. δ + max, ηe λ δ + max, ηe λ δ λ 2eλ maxe, η λ δ. so that, by sing also lemma 5., we have: DΦ Tε ΦTs zεdφts v A Tε DΦ Ts v DΦ Tε Φ Ts zε DΦ T s ATε v ηλ Tε + 2ηe 2 δ T Tε ε =0 Φ Φ Ts DΦ T s v We have η maxe, η 2e2 λ 2 λ λtε δ DΦ T s v. DΦ Ts+Tε v A Tε w DΦ Tε ΦTs zεdφts v A Tε DΦ Ts v + A Tε DΦ Ts v w η maxe, η 2e2 λ 2 λ λtε δ DΦ Ts v + A Tε DΦ Ts v w η maxe, η 2e2 λ 2 λ λtε δ w + DΦ T s v w + λ T ε DΦ T s v w 9

20 and sing 0 and w = w we obtain DΦ T s+t ε v A Tε w λ T ε and, since λ Ts ε δ 2 δ : DΦ Ts+Tε v A Tε w + w s λ Tε δ 5 Two Technical Lemmas η maxe, η 2e2 λ 2 λ δ + λts ε + λ Ts ε w η maxe, η 4e2 λ 2 λ + w. 47 In this Section we prove two technical Lemmas which we obtain by sing Lipschitz ineqalities for Φ and DΦ. Lemma 5. Let U R n be a neighborhood of 0, and Φ : U R n be a smooth map: Φz = Az + fz with f i 0,..., 0 = 0, fi z j 0,..., 0 = 0 for any i, j and, for any z BR, satisfying fz η z 2, Df z η z, DΦ z l. Then, for any z, K sch that Φ z BR for any = 0,..., K, we have and DΦ z A η j=0 λ j Φ j z λ + η Φ j+ z... λ + η Φ z 48 DΦ z A η λ + η max j Φ j z Φ j z. 49 j=0 Proof of Lemma 5.. For = we have DΦ z A = Df z η z. For generic K, since A = λ, we have DΦ z A DΦ Φ zdφ z A = DΦ Φ zdφ z A + DΦ Φ z AA DΦΦ z DΦ z A + DΦΦ z A A = A + Df Φ z DΦ z A + DΦ Φ z A A λ + η Φ z DΦ z Assming that 48 is valid for, we have DΦ z A λ + η Φ z 2 η +η Φ z λ 2 = η j=0 = η j=0 λ j λ j From 48 we immediately obtain 49. j=0 λ j A + η Φ z λ. Φ j z λ + η Φ j+ z... λ + η Φ 2 z Φ j z λ + η Φ j+ z... λ + η Φ z + η Φ z λ Φ j z λ + η Φ j+ z... λ + η Φ z. 20

21 Lemma 5.2 Let U R n be a neighborhood of 0, and Φ : U R n be a smooth map with finite Lipschitz constants λ Φ, λ DΦ for Φ and DΦ respectively. For any initial conditions z 0, z 0, their time evoltions z = Φ z 0, z = Φ z 0 satisfy z T z T λ T Φ z 0 z 0 50 and for any v 0, the time evoltion of the tangent vectors satisfies with where σ = min z U min v = DΦ z v. v T = DΦ T z 0 v, v T = DΦ T z 0 v v T v T v T λ T z 0 z 0. 5 λ = max λ Φ, DΦ + λ DΦ σ Proof of Lemma 5.2. We prove 50 by indction on T. If T = we have Le s assme z T z Then, we have v = z z = Φz 0 Φz 0 λ Φ z 0 z 0. T λ T Φ z 0 z 0. z T z T = Φz T Φz T λφ z T z Then, let s prove 5 by indction on T. If T =, we have v v DΦz 0 DΦ z v DΦz 0 0 DΦ z v 0 v = v T v v = λ T Φ z 0 z 0. DΦz 0 DΦ z By Lipschitz estimate and ineqality 8 we have: DΦz 0 DΦ z v 0 DΦz v 0 DΦ z 0 λdφ z 0 z DΦz v 0 min DΦz v v 0 = σ > 0, v = and therefore we obtain v v v λ DΦ σ z 0 z 0 λ z 0 z 0. v 0 v v 0 DΦz v. 0 We now assme that 5 is satisfied for T, that is: v T v T v λ T z 0 z T Then, let s consider v T v T v T = DΦ z T v T DΦ z T v v T T 2

22 DΦ z T v T v T DΦ z T DΦ z v T T v T + v T v sp DΦ z T v T v T z v T + DΦ z T DΦ z T v T v T = sp DΦ z v T v T z v v T + λ DΦ z T z T v T v T. T Using 8: v T v T = v T DΦ z v T and 50, 52, we obtain v T v T v T = sp z DΦ z T σ sp z DΦ z λ T z 0 z 0 + λ DΦ σ σ λ T + λ DΦ λ T Φ σ λt Φ z 0 z z 0 z 0 λ T z 0 z Conclsions In this paper we have explained why the FLI indicators, sitably modified by the introdction of test fnctions, may be sed for high precision comptations of the stable and nstable manifolds of dynamical systems, inclding the critical comptations of the so called tbe manifolds of the restricted three body problem. An advantage of the FLI method is that it does not reqires a preliminary high precision localization of the hyperbolic fixed points or periodic orbits to provide high precision comptations of their stable and nstable manifolds. This is particlarly sefl for practical applications, since additional pertrbations can be easily inclded in the nmerical comptations. Acnowledgments Part of the comptations have been done on the Mesocentre SIGAMM machine, hosted by the Observatoire de la Cote d Azr. References [] V.I. Arnold, Instability of dynamical systems with several degrees of freedom. Sov. Math. Dol., 6: , 964. [2] Benettin G. Galgani L. and Strelcyn J.M. Kolmogorov entropy and nmerical experiments. Physical Review A, Vol. 4, n. 6, , 976. [3] Celletti A., Lega E., Stefanelli L. and Froeschlé C. Some reslts on the global dynamics of the reglarized restricted three body problem with dissipation. Cel. Mech. and Dyn. Astr., 09, , 20. [4] P. Cincotta, C. Simó, Simple tools to stdy global dynamics in non-axisymmetric galactic potentials - I. Astron. Astrophys. Sp. 47,

23 [5] C. Froeschlé, M. Gzzo and E. Lega, Graphical Evoltion of the Arnold Web: From Order to Chaos. Science, 289, n. 5487: [6] C. Froeschlé, M. Gzzo and E. Lega, Local and global diffsion along resonant lines in discrete qasi integrable dynamical systems. Cel. Mech. and Dyn. Astron., 92, -3: , [7] C. Froeschlé, E. Lega, and R. Gonczi. Fast Lyapnov indicators. Application to asteroidal motion. Celest. Mech. and Dynam. Astron., 67: 4 62, 997. [8] M. Gzzo, The web of three planets resonances in the oter Solar System. Icars, vol. 74, n.., , [9] M. Gzzo M., The web of three-planet resonances in the oter solar system II: a sorce of orbital instability for Urans and Neptne. Icars, 8, , [0] Gzzo M., Chaos and diffsion in dynamical systems throgh stable nstable manifolds, in Space Manifolds Dynamics: Novel Spaceways for Science and Exploration, proceedings of the conference: Novel spaceways for scientific and exploration missions, a dynamical systems approach to affordable and sstainable space applications held in Fcino Space Centre Avezzano 5 7 October Editors: Perozzi and Ferraz Mello. Springer [] M. Gzzo, E. Lega and C. Froeschlé, On the nmerical detection of the effective stability of chaotic motions in qasi-integrable systems. Physica D, 63, -2: [2] M. Gzzo, E. Lega and C. Froeschlé, First Nmerical Evidence of Arnold diffsion in qasi integrable systems. DCDS B, 5, 3: [3] Gzzo M., Lega E. and Froeschlé C., A nmerical stdy of the topology of normally hyperbolic invariant manifolds spporting Arnold diffsion in qasi-integrable systems. Physica D, 82, , [4] M. Gzzo M., E. Lega and C. Froeschlé, First nmerical investigation of a conjectre by N.N. Nehoroshev abot stability in qasi-integrable systems. Chaos, 2, Isse [5] M. Gzzo M., E. Lega, On the identification of mltiple close-enconters in the planar circlar restricted three body problem. Monthly Notices of the Royal Astronomical Society, 428, , 203. [6] M. Gzzo M., E. Lega, The nmerical detection of the Arnold web and its se for long-term diffsion stdies in conservative and wealy dissipative systems, Chaos, vol. 23, 02324, 203. [7] Hénon M. and Heiles C.: The Applicability of the Third Integral of Motion: Some Nmerical Experiments. The Astronomical Jornal, 69, p , 964. [8] Koon W.S., Lo M.W., Marsden J.E. and Ross S.D. Dynamical Systems, the three body problem and space mission design. Marsden Boos. ISBN , [9] J. Lasar. The chaotic motion of the Solar System. A nmerical estimate of the size of the chaotic zones. Icars, 88:266 29, 990. [20] J. Lasar, C. Froeschlé, and A. Celletti. The measre of chaos by the nmerical analysis of the fndamental freqencies. Application to the standard mapping. Physica D, 56:253, 992. [2] J. Lasar. Freqency analysis for mlti-dimensional systems. Global dynamics and diffsion. Physica D, 67:257 28, 993. [22] E. Lega, M. Gzzo and C. Froeschlé, Detection of Arnold diffsion in Hamiltonian systems. Physica D, 82:

24 [23] E. Lega, M. Gzzo and C. Froeschlé, A nmerical stdy of the hyperbolic manifolds in a priori nstable systems. A comparison with Melniov approximations., Cel. Mech. and Dyn. Astron., 07, 5-27, 200. [24] E. Lega, M. Gzzo and C. Froeschlé, Detection of close enconters and resonances in three-body problems throgh Levi-Civita reglarization, Monthly Notices of the Royal Astronomical Society, 48, 07-3, 20. [25] T.A Mitcheno and S. Ferraz Mello, Resonant strctre of the oter solar system in the neighborhood of the planets. A.J. 22, , 200. [26] P. Robtel, Freqency map analysis and qasiperiodic decompositions, in Hamiltonian systems and Forier analysis, Editor: Benest et al., in Hamiltonian systems and Forier analysis, 7998, Adv. Astron. Astrophys., Camb. Sci. Pbl., Cambridge [27] P. Robtel and J. Lasar, Freqency map and global dynamics in the Solar System I. Icars, [28] P. Robtel and F. Gabern, The resonant strctre of Jpiter s Trojan asteroids I. Long term stability and diffsion. Monthly Notices of the Royal Astronomical Society., [29] C. Simó, Dynamical systems methods for space missions on a vicinity of collinear libration points, in Simó, C., editor, Hamiltonian Systems with Three or More Degrees of Freedom S Agaró, 995, volme 533 of NATO Adv. Sci. Inst. Ser. C Math. Phys. Sci., pages , Dordrecht. Klwer Acad. Pbl., 999. [30] Szebehely V. Theory of orbits. Academic Press, New Yor, 967. [3] X.Z. Tang and A.H. Boozer, Finite time Lyapnov exponent and advection-diffsion eqation. Phys. D, 95, 3-4, [32] Villac B.F., Using FLI maps for preliminary spacecraft trajectory design in mlti-body environments. Cel. Mech. and Dyn. Astron., 02, 29-48, [33] B.H. Wayne, A.V. Malyh and C.M. Danforth, The interplay of chaos between the terrestrial and giant planets. Monthly Notices of the Royal Astronomical Society Volme 407, Isse 3, September 200, Pages: