A COMPARISON BETWEEN METHODS TO COMPUTE LYAPUNOV EXPONENTS G. TANCREDI, A.SA NCHEZ
|
|
- Solomon Armstrong
- 5 years ago
- Views:
Transcription
1 THE ASTRONOMICAL JOURNAL, 121:1171È1179, 21 February ( 21. The American Astronomical Society. All rights reserved. Printed in U.S.A. A COMPARISON BETWEEN METHODS TO COMPUTE LYAPUNOV EXPONENTS G. TANCREDI, A.SA NCHEZ Departamento de Astronom a, Facultad de Ciencias, Universidad de la Repu blica, Igua 4225, 114 Montevideo, Uruguay; gonzalo=ðsica.edu.uy, andrea=ðsica.edu.uy AND F. ROIG Instituto Astronoü mico e Geof sico, Universidade Sa8 o Paulo, Avenida Miguel Este fano 42, Sa8 o Paulo, Brasil; froig=iagusp.usp.br Received 2 January 25; accepted 2 October 31 ABSTRACT Lyapunov characteristic exponents measure the rate of exponential divergence between neighboring trajectories in the phase space. For a given autonomous dynamical system, the maximum Lyapunov characteristic exponent (hereafter LCE) is computed from the solution of the variational equations of the system. There are many dynamical systems in which the formulation and solution of the variational equations is a cumbersome task. In those cases an alternative procedure, Ðrst introduced by Benettin et al., is to replace the variational solution by computing two neighbor trajectories (the test particle and its shadow) and calculating the mutual distance. In this paper, we deal with a comparison between these two di erent techniques for the calculation of LCE: the variational method and the two-particle method. We point out a problem that can appear when the two-particle method is used, which can lead to a false estimation of a positive LCE. The explanation of this phenomenon can be analyzed in two di erent situations: (1) for relatively large initial separations the two-particle method is not a good approximation to the solution of the variational equations, and (2) for small initial separations the two-particle method have problems related to the machine precision, even when the separation can be many order of magnitudes larger than the machine precision. We show some examples of false estimates of the LCE that have already appeared in the literature using the two-particle method, and Ðnally we present some suggestions to be taken into account when this method has to be used. Key words: methods: numerical È solar system: general 1. INTRODUCTION where d(t) is the distance in the phase space between two trajectories (a reference trajectory and a shadow one), initially separated by a distance d. For the result to be mean- Lyapunov characteristic exponents are a widely used tool for the estimation of chaoticity in dynamical systems. They ingful, the exponential divergence should be measured in basically measure the rate of exponential divergence the immediate neighborhood of the reference trajectory; between neighboring trajectories in the phase space. For a thus the LCE should now be deðned as given autonomous dynamical system x5 \ f (x), the maximum Lyapunov characteristic exponent (hereafter c \ lim s(t). (3) LCE) is deðned as c \ lim s(t), with t?= t?= d? s(t) \ 1 o m(t) o Then, for practical calculations, a periodic renormalization ln, (1) of the distance d(t) is necessary in order to avoid the shadow t o m o to remain far away from the reference trajectory. The renormalization method introduced by Benettin et al. (1976) con- where m(t) is the solution of the variational equations of the system, m5 \ + fém and m \ m(). Then the system is said sists in approaching the shadow trajectory to the reference to be regular x (integrable) if c \ and chaotic if c[. The one along the separation vector down to the initial distance numbers s(t) are called Lyapunov characteristic indicators d. If this renormalization is made at a Ðxed time interval q, (LCI; Froeschle, Froeschle, & Lohinger 1993) and for practhen we can write tical calculations it is not possible to take t ] O. Thus, the evolution of s(t) is followed up to some time and we plot s \ 1 k kq ; k ln d(iq), (4) d log s(t) versus log t. If the curve has a negative constant i/ slope, the system is regular; if it presents an inñection of the and thus slope which comes close to and the plot converges to a c \ lim s. (5) certain value of s(t), the system is chaotic. k There are many dynamical systems in which the formulation and solution of the variational equations is a cumber- The renormalization is also necessary to avoid k?= d? some task. In those cases, an alternative procedure to saturation,ïï which is caused because a chaotic region is calculate c, Ðrst introduced by Benettin, Galgani, & normally bounded in the phase space, and thus the distance Strelcyn (1976), is to replace the deðnition in equation (1) by between the two trajectories cannot grow by more than the s(t) \ 1 average size of the chaotic region. It is worth noting that all d(t) ln, (2) these problems are avoided with the solution of the variational equations, since it is a linear solution, and so its t d 1171
2 1172 TANCREDI, SAŠ NCHEZ, & ROIG Vol. 121 behavior is independent of the initial condition and of the size of the chaotic region (for a discussion see, for example, Wisdom 1983). Both the solution of variational equations (hereafter variational method) and the simultaneous integration of two trajectories (hereafter two-particle method) have been extensively used in the literature to calculate the LCE. The aim of this paper is to point out a problem that can appear when the two-particle method is used, which can lead to a false estimation of a positive LCE ( 2). In 3we discuss the dependence of the false estimates of the LCEs with di erent initial parameters, and we present a conjecture to explain this phenomenon. We will show some examples from the literature where false estimates have already appeared ( 4). In the conclusion ( 5) we give some recommendations in case one is forced to use the two-particle method. 2. FALSE ESTIMATE OF LCEs We recall that the Lyapunov characteristic indicators (the LCIs) are Ðnite estimates of the Lyapunov characteristic exponents (the LCEs). We then use the term LCI when we refer to the result of a numerical computation and to LCE when we refer to the expected theoretical value. The main di erence between the variational and twoparticle methods is clear: the solution m(t) is simply the linearized distance d(t). Roughly speaking, both methods can be compared with the approximation of a function obtained by using the tangent at a point or the secant by two points, respectively. There is no a priori reason to think that the results obtained with each method could be di erent. At least they should be qualitatively similar, provided that we are far away from a singularity of the system, in which case the linear approach of the distance is not valid and the deðnition of LCE (eq. [1]) becomes meaningless. However, we can Ðnd many cases of quasi-regular or even regular (integrable) systems in which the results of two methods are not the same. As an example, let us consider the two body problem (a massless particle in a central gravitational Ðeld) and calculate the LCI using the two methods. For this calculation we used the integration algorithm Radau by Everhart (1985). We adapted the integrator to solve either the di erential equations for the particle and its shadow (two-particle method) or the di erential equations for the particle and the variational equations (variational method). For the two-particle method, if the particle has a position and a velocity vector p we create the shadow by, an homothesis of factor (1 ] b), i.e., p 1 \ (1 ] b)p. (6) FIG. 1.ÈLCI for the solution of the two-body problem calculated using the two-particle method and the variational method using the RADAU integrator. The initial separation was b \ d /p \ 1~8 and the renormali- zation time step was q \ 1 yr. techniques, with always the same qualitative result, though there are important di erences in the quantitative results among the set of variables and numerical algorithms used (see 3). Figure 2 presents curves log m versus log t for di erent numerical algorithms; the calculations were done for the same initial conditions, renormalization time step, and The initial distance d is then simply d \ b o p o. (7) The initial orbits for the following calculations in the two body problem correspond to the orbit of Ceres around the Sun (semimajor axis a \ 2.77 AU, eccentricity e \.76). As we can see in Figure 1, the variational method gives the actual result, but the two-particle method gives a positive LCI, which is clearly wrong. We made several experiments using di erent values of the distance, d, and of the renor- malization time step, q. We used also di erent sets of variables and coordinate systems and even di erent integration FIG. 2.ÈLCI for the solution of the two-body problem calculated using the variational equations with renormalization (nearly straight downward line), and the two-particle method with di erent numerical integration algorithms: Radau (top curve), the Keplerian algorithm of this paper (middle curve), SWIFT LYAP (bottom curve). For all cases d /p \ 1~8 and q \ 1 yr.
3 No. 2, 21 METHODS TO COMPUTE LYAPUNOV EXPONENTS 1173 initial separation. We compare results using the Radau integrator presented above, the subroutine to compute Lyapunov exponents with the two-body method that comes with the Swift package (SWIFT LYAP; Levison & Duncan 1994),1 and a very simple numerical algorithm based on the solution of KeplerÏs equation, hereinafter called Keplerian algorithm. This algorithm consists in the following steps: 1. Given the position and velocity of the particle, we create the shadow by an homothesis of factor (1 ] b), as explained above. 2. The initial mutual distance d between the particle and shadow is computed. 3. The position and velocity of particle and shadow are transformed into orbital elements a semimajor axis, e eccentricity, u argument of perihelion, M mean anomaly. 4. During a time step q, a, e and u remain constant, and the mean anomaly evolves as: M i`1 \ M i ] nq, (8) where n is the mean motion. 5. After solving KeplerÏs equation, the orbital elements are transformed back to position and velocity in Cartesian coordinates. The mutual distance at the end of the time step is then computed. 6. Afterward the renormalization procedure is applied, the shadow particle is brought back to a distance d along the separation vector in position and velocity. Return to point 3. This method is obviously much faster than the numerical solution of the gravitational di erential equation, letting us make many tests in short CPU time. Note that all numerical algorithms that use the twoparticle method give positive values of the LCI (see Fig. 2). The other alternative is to apply the renormalization procedure whenever the separation distance reaches a Ðxed value, rather than at Ðxed time steps. The results of applying this method are shown in Figure 3 for two di erent initial separations. The renormalization was done when the distance reached 1 times its original value. For similar initial separation we can compare the results with those of Figure 1 to conclude that, though we get a smaller positive LCI, we still obtain a false estimate Some T entative Explanations Holman & Murray (1996) presented a possible explanation for the overestimation of LCEs observed in quasiregular systems when the two-particle method is used. They argued that the false positive LCE is generated in the twoparticle method because after the i-th renormalization the mutual distance d between the particles is given approx- imately by i d i ^ d (1 ] aqn) exp (jq), (9) where d is the initial distance, q is the interval between successive renormalizations, a and n are constants associated with the initial power-law transient separation, and j is ÈÈÈÈÈÈÈÈÈÈÈÈÈÈÈ 1 Available at FIG. 3.ÈLCI for the solution of the two-body problem calculated using the two-particle method with renormalization at Ðxed separation distance. The initial separation was b \ d /p \ 1~6 for the top curve and b \ d /p \ 1~8 for the lower one. The renormalization was done whenever the distance reached hundred times its original value. the maximum LCE. Combining this expression with the renormalization formula (eq. 4) they obtained ln 1 ] (aqn) s ^ j ]. (1) k q Then, they claimed that for those trajectories with j D or j > ln (1 ] aqn)/q (i.e., quasi-regular trajectories) there is an overestimation of the LCE that is larger for smaller values of q. This does not a ect the strong chaotic trajectories, since in that case j? ln (1 ] aqn)/q. Although this is in good agreement with the behavior observed in their experiments, it is possible to prove that their explanation is not valid. This is because a in equation (9) is not actually a constant as they assumed, but rather it depends on k. We will show this in the following. Consider that the modulus of the distance between the two trajectories evolves with time following the law d(t) \ d (1 ] j r tn) exp (j c t). (11) In this expression, j and j are constants associated with the rate of divergence r of the c regular and chaotic motion, respectively. If the trajectory is regular, then j \ and the distance will grow according to (1 ] j tn). In c other words, an initial error d in the reference trajectory r will propagate following a polynomial of order n (it is common to say that the error propagates linearly with time). On the other hand, if the trajectory is chaotic, j [. Thus, for very short times the polynomial divergence dominates, c but for longer times the exponential divergence becomes more important. It is clear that j is the maximum LCE. Now, suppose c that we make a renormalization at a Ðxed interval q and, in between the steps, the distance evolves
4 1174 TANCREDI, SAŠ NCHEZ, & ROIG Vol. 121 following the function d(t) deðned in equation (11). Then at t \ q the distance before renormalization is d \ d(q); at t \ 2q is d \ d d(2q)/d(q); consequently at 1 t \ (i ] 1)q would be 2 d[(i ] 1)q] d \ d i`1 d(iq) or in terms of equation (11), (12) [1 ] j qn(i ] 1)n] exp [j (i ] 1)q] d \ d r c. (13) i`1 (1 ] j qnin) exp (j iq) r c Using the binomial theorem, we can write this expression as d \ d (1 ] j qna ) exp (j q), (14) i`1 r i`1 c where a \ 1 i`1 1 ] j qnin ; n An r jb in~j. (15) j/1 It is simple to see that lim a \, since a D 1/i for large i?= i i values of i. Applying equation (14) to the renormalization formula (eq. [4]), we obtain s \ 1 k kq ; k ln [(1 ] jr qna ) exp (j q)] \ j ] S, (16) i c c k i/1 where S \ 1 k kq ; k ln (1 ] j qna ). (17) r i i/1 LetÏs show that lim S \. First, we have that for k?= k large values of i we can write ln (1 ] j r qna i ) D ln A 1 ] j r qn 1 ib ^ jr qn 1 i and thus for large values of k we obtain ] ÉÉÉ, (18) S D 1 k kq ; k j qn r ] ÉÉÉ \ ln (j r qnk), (19) i kq i/1 which clearly goes to zero for k ] O. Then, c \ lim s \ j. If k?= a is k assumed c to be constant ( \ a), equation (16) becomes i identical to equation (1) (the one assumed by Holman & Murray 1996); but we have just demonstrated that the assumption is not valid, a goes to zero for large k. It is worth noting that an expression i similar to equation (16) can be found when applying the variational method. Note that the modulus of the solution of the variational equations can be also represented by an expression similar to equation (11), i.e., m(t) \ m (1 ] j r tn) exp (j c t), (2) where j is associated with the maximum eigenvalue of the system and c j and n are related to the degeneracy of some eigenvalues. Then, r a periodic renormalization of m(t) combined with equation (4) (now substituting d by m), leads again to equation (16). Since in practise the renormalization procedure is carried out up to a Ðnite k, then a will always have a Ðnite nonzero value. In this case, the overestimation k of the LCE observed in quasi-regular systems could be justiðed in the framework of the explanation given by Holman & Murray (1996). However, since lim a \, this eventual overestimation should decrease (or k?= even k disappear) as soon as we extend our calculation to larger values of k, i.e., advancing the trajectories for a larger time span keeping the same q between renormalizations. Of course this is not the behavior observed in quasi-regular trajectories when using the two-particle method. In fact the false positive LCE is independent of the total time span of the calculation. On the other hand, if we assume for a moment that in practice ÏÏ a is almost constant for any k, as Holman & Murray (1996) k assumed, then we should also observe an overestimation of the LCE when applying the renormalization formula (eq. [4]) in the variational method, since the asymptotic behavior of s is the same for both methods (eq. [16]). Again, this is not what k is observed in the experiments. For example, in Figure 2 the almost downward straight line corresponds to a calculation where we apply the renormalization procedure to the variational method exactly in the same way as for the two-particle method. It is possible to verify that the variational method never produces false positive LCEs. Then, we conclude that: (1) the explanation given by Holman & Murray (1996) is not necessary correct since the assumption a \ const. is wrong, and (2) even assuming that in practise a could be e ectively constant, their explanation does not correctly account for the discrepancy observed between the two-particle and variational method, respectively. We stress the fact that, based on the explanation given by Holman & Murray (1996), the rescaling technique should also lead to a false LCE when the variational solution is used. In the next section, we will show that, in fact, the origin of the false estimates of LCEs in the two-particle method relies in the accumulation of round-o errors during the computation of the mutual distances between the two particles in the course of the successive renormalizations. These errors are not present in the variational solution since this does not involve the calculation of mutual distances. 3. DEPENDENCE OF LCI WITH d AND q Before formulating a conjecture to explain the false estimation of LCE in regular systems when using the twoparticle method, let us study the dependence of the LCI with the initial distance between the two trajectories d. We use the Keplerian algorithm for the following calculations. Figure 4 shows the LCI of a two-body problem calculated using the two-particle method with di erent values of d. Although all the LCIs give wrong estimates (they should be zero), we note that for very large and very small values of d the estimation is worse, and there seems to be an optimal ÏÏ value of d /p \ 1~8 for which the estimation is better. A more detailed analysis of the dependence of the LCIs with d is presented in Figure 5. The two-particle method is again applied to the two-body problem for a wide set of initial separation. The Ðnal value of s at t \ 17 yr is plotted against log (d /p ). The curves from top to bottom correspond to calculations in: (1) double precision with renormalization time step, q \ 1 yr; (2) double precision with q \ 1 yr; and (3) quadruple precision with q \ 1
5 No. 2, 21 METHODS TO COMPUTE LYAPUNOV EXPONENTS 1175 FIG. 4.ÈLCI for the solution of the two-body problem calculated using the two-particle Keplerian algorithm of this paper. Di erent curves correspond to di erent initial separations d. The renormalization time step was 1 yr in all cases. yr. [Note that for d /p \ 1~1,1~15, and 1~2 we just give upper limits to the LCEs, since the curves log s(t) versus log t do not level o.] Two distinct regimes can be distinguished: one for large values of d (the right-hand branch of the V-shaped curve in Fig. 5) and another one for small values of d (left branch); the change from one regime to the other depends on the precision of the calculation and slightly on the renormalization time step. For the right branch, larger values of s correspond to larger d. This is because for large values of d we are far away from the condition of neighboring trajectories. The distance, d, is no more equivalent to its linear approach m and the deðnition equation (2) is no longer applicable (recall that the deðnition of LCEs uses m not d!). In the deðnition of the LCE, the distances are measured on the surface tangent to the manifold, while in the two-particle method the tangent vector is approximated by a secant. This secant vector belongs to another space, which is itself tangent to a manifold di erent than the original one. This manifold deðnes a di erent dynamical problem, and then its associated LCE does not need to be equal to that of the original problem. Under this interpretation, it is not expected that the spurious LCIs have any dependence on the precision of the calculations nor on the renormalization time step, because it is just a geometrical problem. In fact, the convergence of the right-hand branches of the middle and bottom curves proves that the machine precision has no e ect on the spurious values. Also, we see that for very large values of d, (e.g., for d /p Z 1~3), a good agreement in the Ðnal value of the LCE is found for di erent renormalization steps (see Fig. 5 and the three upper curves of Fig. 6). These numerical results are in disagreement with the interpretation by Holman 7 Murray (1996), since they would expect a dependence with the renormalization step. As it is expected from the deðnition of the LCEs, the spurious values become closer to the true value (i.e., zero) as d ]. This tendency is at least true up to a certain d, where the curve s versus d has a turning point (see Fig. 5). We found that in our particular problem the turning point occurs for values of d many orders of magnitudes higher FIG. 5.ÈLCI for the solution of the two-body problem as a function of the initial separation d. The calculation was made with the Keplerian algorithm of this paper, using di erent renormalization time steps and di erent machine precision. For some calculations we give only an upper limit of the estimated LCI (dotted line). See text for explanation. FIG. 6.ÈLCI for the two-body problem calculated with the Keplerian algorithm for di erent values of d and q. The top three curves correspond to d /p \ 1~3 and the three bottom ones correspond to d /p \ 1~1.
6 1176 TANCREDI, SAŠ NCHEZ, & ROIG Vol. 121 than the machine precision. When using a renormalization step of 1 yr, the turning point occurs at d /p D 1~8 in double precision calculations (where the machine precision is D 1~16)2; while it occurs at d /p D 1~15 in quadruple precision calculations (where the precision is D 1~34). To understand the cause of this problem we will make the following experiments with the results of the Keplerian algorithm: let us project the separation vector connecting the particle and its shadow on the phase space of orbital elements a, e, u, M, where a is the semimajor axis, e the eccentricity, u the argument of perihelion, and M the mean anomaly. After a certain time interval, the separation vector is generally in the direction of the mean anomaly axis, and hence its projection dm on this axis after renormalization has in modulus almost the size of the initial separation. On ÈÈÈÈÈÈÈÈÈÈÈÈÈÈÈ 2 We deðne machine precision as the quantity p for which (1 ] p) [ 1is consider to be equal to. the other hand, the projections in the other axis are many orders of magnitudes smaller. Between two renormalizations, a, e, and u su er no changes, while M evolves linearly with time (see the description of the algorithm in 2). Then, after a renormalization time q we can write da \ da, (21) de \ de, (22) du \ du, (23) while the separation dm is the sum of two components dm \ dm ] qdn \ dm [ 3 q da. (24) 2 a5@2 (For the particle we use the following standard units: semimajor axis a \ 1, mean motion n \ 1, period T \ 2n). In Figure 7 we present the evolution of da and dm with time for an initial distance d /p \ 1~1 and two renorma- FIG. 7.ÈSeparation vector d is projected onto the phase space of orbital elements a, e, u, M and the components da and dm are shown here as a function of time. (a) For a renormalization time step of 1 yr; (b) for a renormalization time step of 1 yr. da is expressed in units of 1~13 and dm in 1~9 radian. Calculations were made with the Keplerian algorithm for an initial separation d /p \ 1~1.
7 No. 2, 21 METHODS TO COMPUTE LYAPUNOV EXPONENTS 1177 lization steps q \ 1 and 1 yr. As mentioned above dm is of the order of d. Considering that a for the particle and the shadow are on the order of 1, the di erence da could not be smaller than the machine precision. Again we can see in Figure 7 that da is generally of that order, but it presents random Ñuctuations on signiðcant digits that, following equation (24), produce spurious contributions to dm. The random Ñuctuations are generated in the renormalization procedure because the shadow is approached to the particle along the separation vector (in Cartesian coordinates), and this separation tends to be aligned with the orbit so da ] (numerically means, to the machine precision). But, depending on the position on the orbit, the values of da could be much larger than the machine precision. The accumulation of these spurious values cause a false estimation of the LCE. For a smaller d, the contribution of da on the values of dm is larger, and hence larger LCIs are obtained. Since these spurious values are generated in the renormalization procedure, it is expected that the Ðnal values of the LCIs depend also on the renormalization time step. This is clear if we look at the left-hand branches of the top and middle curves in Figure 5, and also the three lower curves in Figure 6, where we can see a strong dependence of the spurious LCIs with the renormalization time step. Comparing the evolution of da and dm with time for q \ 1 and 1 yr (Fig. 7), we conclude that the shorter the renormalization time step the larger the number of spurious values of da and hence of dm. This is obvious since, for a given time span, there are more renormalizations for a shorter renormalization time step, and there are more chances for large spurious values to occur. For the numerical integrations using the gravitational di erential equations instead of the Keplerian algorithm, the problem is slightly di erent. In Figure 1 we present results for a single value of q and d. Similar plots as the ones presented in Figure 5 can be obtained in the integrations of the gravitational di erential equations, with a V - shaped dependence of the LCI with d. The explanation for the right-hand branch is the same as described above. In the left-hand branch the role played by da is now played by the separation of the velocity vector dv. The velocity vectors for the particle and the shadow tend to be aligned, i.e., o dv Š o ]. Since dr is on the order of d, the renormaliza- tion procedure introduces large Ñuctuations on dv. Consider that VŠ 2\k A2 r [ 1 ab ; (25) hence A 2V dv \ k [ 2 r2 dr [ da a2b. (26) Large spurious values of dv would imply large values of da and hence of dm, leading to values of dr? dr. (Note that da and dm are not calculated in the numerical integrations, but obviously they are implicitly a ecting dr and dv.) The accumulation of these spurious values cause the false estimate of the LCE. 4. EXAMPLES IN THE LITERATURE The two-particle method has been extensively used in many research Ðelds for the computation of the Lyapunov exponents. As we mentioned earlier, it is much easier to implement and in some cases it is the only alternative available. As an example we present a few cases where this method has been used in the study of the orbital dynamics of solar system bodies; e.g., Sussman & Wisdom (1988), Franklin, Lecar, & Soper (1989), Lecar, Franklin, & Murison (1992), Whipple & Shelus (1993), Franklin, Lecar, & Murison (1993), Saha & Tremaine (1993), Murison, Lecar, & Franklin (1994), Franklin (1994), Mu ller & Dvorak (1995), Tancredi (1995, 1998). These papers deal with many di erent problems, like the chaotic behavior of Pluto, main belt asteroids, NEAs, Jupiter family comets, and outer Jovian satellites. Holman & Murray (1996) already showed that the computations of the Lyapunov exponents of some asteroids in the outer belt made by Murison et al. (1994) with the twoparticle method were mostly wrong. Murison et al. (1994) results were correct only in the cases where the Lyapunov times were less than a few thousand years (the Lyapunov time is the inverse of the Lyapunov exponent). Their shadow orbits were started at a distance of 1~6 in phase space. In Figure 5 we observe that for this initial separation, spurious values of the Lyapunov time on the order of a few thousand years are obtained depending on the renormalization time. The two-particle method is not suitable to analyze the slow chaos observed in the asteroid main belt. Sussman & Wisdom (1988) used the two-particle method to compute the Lyapunov time of Pluto but they were very careful in the implementation. They used a high-precision computer (which allows them to use initial separations of the order of 1~2); but they also tested their results with the variational method. We repeat their calculations with the variational method and we obtain a Lyapunov time of the same order. Tancredi (1995) and Tancredi (1998) studied the dynamical evolution of Jupiter family comets and NEAs, respectively. He obtained Lyapunov times clustered in the range 5È15 yr and generally lower than 1 yr. As observed in Figure 5, for such high values of the LCI (1~2È1~3 yr~1), the two-particle method do not introduce errors for initial relative separation smaller than 1~4, as he had used. Whipple & Shelus (1993) and Saha & Tremaine (1993) studied the chaotic evolution of the outer Jovian satellites. The Ðrst authors centered on J-VtIII Pasiphae, while the later authors considered the four retrograde satellites (J- VIII Pasiphae, J-IX Sinope, J-XI Carme, J-XII Ananke). The integrations of Whipple & Shelus (1993) were not long enough to estimate a positive Lyapunov exponent. Figure 2c of Saha & Tremaine (1993) shows the variation of c(t) with time in a logarithmic scale. As stated by the authors, All four curves show some sign of levelling o, but only Sinope is clearly chaotic.ïï The integrations by Saha & Tremaine (1993) were performed with a symplectic algorithm, they used a initial separation of neighboring orbits of da/ a \ 1~4 and a renormalization time of 512 yr. They also mentioned, A better procedure would have been to integrate the tangent equations, but the integrator we use is not suited to this task.ïï We reproduce their calculations using the same initial conditions and two di erent algorithms: (1) the Radau integrator and the variational equations; (2) the Radau integrator and the two-particle method with the distances calculated in the elements phase space. With the variational method we obtain a positive LCI for Pasiphae and Sinope, but no levelling o for Carme and
8 1178 TANCREDI, SAŠ NCHEZ, & ROIG Vol. 121 Ananke (Fig. 8a). With the two-particle method, a similar LCE is obtained for Sinope; for Pasiphae we obtain a slightly higher LCE, and the curves for Carme and Ananke show signs of levelling o similar to the results obtained by Saha & Tremaine (1993) (compare their Fig. 2c and our Fig. 8b). One should also compare the slopes of the corresponding curves in Figures 8a and b and note that the slopes in the two-particle method are higher than in the variational method. The two-particle method gives false estimates of the Lyapunov exponent higher than the actual values. As stated by the authors the implementation of the variational method is not a straightforward task in a symplectic integrator, whereas it is much easier in a general integrator of ordinary di erential equations such as Radau or Bulirsh- Stoer. In recent years the simplectic integrator has been widely used within the solar system dynamicist community, in particular after the public release of the package Swift developed by Levison & Duncan (1994). This package includes two codes to compute the LCI: one with the twoparticle method that it is not recommended and the other implements the variational method but it is still under test.3 Recently, Murray & Holman (1999) measured the Lyapunov exponent of the system composed by the Sun ] the giant planets by comparing pairs of integrations with very small di erences in the initial conditions. They introduced a shift in the coordinate of Uranus of the order of 1~14 AU. By using this tiny separation they hoped to avoid the requirement of renormalization. Nevertheless, since they obtained a Lyapunov time of D 7 ] 16 yr in a D 2 ] 18 yr integration, following equation 11 the Ðnal separation should be of the order of.25 AU. Two criticism can be pointed out to this procedure to calculate LCIs: (1) with such tiny initial separations that are less than two order of magnitudes greater than the machine precision, the errors ÈÈÈÈÈÈÈÈÈÈÈÈÈÈÈ 3 Available at in the integration algorithm becomes very relevant; in particular if they claim a relative error in energy of the order of 1~9; (2) though, one would expect that the separation at the end of the integration would mainly be along the angular variables, it may be large enough to question the validity of the assumption of sampling the same chaotic region, in particular on view of the quasi-resonance structure of the system. A detailed analysis of the results of Murray & Holman (1999) is out of the scope of this paper, we just point out the inconvenience to avoid renormalization. 5. CONCLUSIONS We have shown in this paper that in certain cases it is possible to obtain a spurious estimation of the LCEs when using the two-particle renormalization method of Benettin et al. (1976). This behavior can be explained by two major causes: (1) When the initial distance d /p is larger than a certain value (of the order of 1~8 in double precision in our integrations of heliocentric orbits) the deðnition of distance between two neighboring orbits over the manifold does not match the deðnition of linearized distance over the tangent space to the manifold. Then the deðnition of LCE as given by equation (2) is no longer valid. (2) On the other hand, if d /p is much smaller than this value, the renormalization procedure for the two particles is able to introduce large numeric errors in the calculation of the mutual distance; e.g., the value of some parameters can fall near to the machine precision and they can indirectly a ect the value the principal components of the mutual distance, leading to spurious large values. Since the errors are introduced in the renormalization procedure, it is natural to expect a dependence of the estimated LCEs with the number of renormalizations performed (i.e., with the renormalization time step) in the sense that the smaller the step the worse the estimation. However, the renormalization does not provoke the same e ect in the variational method; because, as we FIG. 8.ÈLCI for the outer Jovian satellites Pasiphae ( P ÏÏ), Sinope ( S ÏÏ), Carme ( C ÏÏ) and Ananke ( A ÏÏ) computed with (left) the Radau integrator and the variational equations; (right) the Radau integrator and the two-particle method with the distances calculated in the elements phase space.
9 No. 2, 21 METHODS TO COMPUTE LYAPUNOV EXPONENTS 1179 explain in 2.1, the exact solution given by the variational equations lead to a correct result even if a renormalization procedure is applied. We can conclude that the two-particle method is not recommended to calculate LCIs in those cases where the solution can fall in a region of regular or quasi-regular motion of the phase space. In these cases the variational method should be preferable. For a region of strong chaoticity, the LCIs calculated with the two-particle method give acceptable values. Nevertheless, it is recommended to make several test with di erent initial separations and renormalization time steps. To be conðdent on the estimate of the LCE, no variation on the Ðnal value of the LCI should be observed among the di erent runs. One of the authors (F. R.) whish to thanks the Sa8 o Paulo State Science Foundation (FAPESP) for partially supporting this work. Benettin, G., Galgani, L., & Strelcyn, J.-M. 1976, Phys. Rev. A, 14, 2338 Everhart, E. 1985, Dynamics of Comets: Their Origin and Evolution, ed. A. Carusi & G. B. Valsecchi (Dordrecht: Reidel), 185 Franklin, F. 1994, AJ, 17, 189 Franklin, F., Lecar, M., & Murison, M. 1993, AJ, 15, 2336 Franklin, F., Lecar, M., & Soper, P. 1989, Icarus, 79, 223 Froeschle, C., Froeschle, Ch., & Lohinger, E. 1993, Celest. Mech. Dyn. Astron., 56, 37 Holman, M., & Murray, N. 1996, AJ, 112, 1278 Lecar, M., Franklin, F., & Murison, M. 1992, AJ, 14, 123 REFERENCES Levison, H. F., & Duncan, M. J. 1994, Icarus, 18, 18 Mu ller, P., & Dvorak, R. 1995, A&A, 3, 289 Murison, M. A., Lecar, M., & Franklin, F. A. 1994, AJ, 18, 2323 Murray, N., & Holman, M. 1999, Science, 283, 1877 Saha, P., & Tremaine, S. 1993, Icarus, 16, 549 Sussman, J., & Wisdom, J. 1988, Science, 241, 433 Tancredi, G. 1995, A&A, 299, 288 ÈÈÈ. 1998, Celest. Mech. Dyn. Astron., 7, 181 Whipple, A. L., & Shelus, P. J. 1993, Icarus, 11, 265 Wisdom, J. 1983, Icarus, 56, 51
On the reliability of computation of maximum Lyapunov Characteristic Exponents for asteroids
Dynamics of Populations of Planetary Systems Proceedings IAU Colloquium No. 197, 25 Z. Knežević and A. Milani, eds. c 25 International Astronomical Union DOI: 1.117/S1743921348646 On the reliability of
More informationLong simulations of the Solar System: Brouwer s Law and chaos
ANZIAM J. 46 (E) pp.c1086 C1103, 2005 C1086 Long simulations of the Solar System: Brouwer s Law and chaos K. R. Grazier W. I. Newman James M. Hyman Philip. W. Sharp (received 28 October 2004, revised 18
More informationHuge Quantum Gravity Effects in the Solar System
Huge Quantum Gravity Effects in the Solar System arxiv:1005.3024v1 [hep-th] 17 May 2010 Don N. Page Theoretical Physics Institute Department of Physics, University of Alberta Room 238 CEB, 11322 89 Avenue
More informationEFFICIENT INTEGRATION OF HIGHLY ECCENTRIC ORBITS BY QUADRUPLE SCALING FOR KUSTAANHEIMO-STIEFEL REGULARIZATION
The Astronomical Journal, 128:3108 3113, 2004 December # 2004. The American Astronomical Society. All rights reserved. Printed in U.S.A. EFFICIENT INTEGRATION OF HIGHLY ECCENTRIC ORBITS BY QUADRUPLE SCALING
More informationDynamical behaviour of the primitive asteroid belt
Mon. Not. R. Astron. Soc. 293, 405 410 (1998) Dynamical behaviour of the primitive asteroid belt Adrián Brunini Observatorio Astronómico de La Plata, Profoeg, Paseo del Bosque, (1900) La Plata, Argentina
More informationDIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES
DIFFUSION OF ASTEROIDS IN MEAN MOTION RESONANCES KLEOMENIS TSIGANIS and HARRY VARVOGLIS Section of Astrophysics, Astronomy and Mechanics, Department of Physics, University of Thessaloniki, 540 06 Thessaloniki,
More informationA SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION
3 A SYMPLECTIC MAPPING MODEL FOR THE STUDY OF 2:3 RESONANT TRANS-NEPTUNIAN MOTION K.G. HADJIFOTINOU and JOHN D. HADJIDEMETRIOU Department of Physics, University of Thessaloniki, 540 06 Thessaloniki, Greece
More informationChaos and stability in the vicinity of a Jovian planet
BEREA COLLEGE Chaos and stability in the vicinity of a Jovian planet by Shiblee Ratan Barua Berea College /22/28 It has been widely known that the influence of large bodies (the Sun, the terrestrial and
More informationTowards stability results for planetary problems with more than three bodies
Towards stability results for planetary problems with more than three bodies Ugo Locatelli [a] and Marco Sansottera [b] [a] Math. Dep. of Università degli Studi di Roma Tor Vergata [b] Math. Dep. of Università
More informationFROM THE SCATTERED DISK TO THE OORT CLOUD The Extended Scattered Disk
FROM THE SCATTERED DISK TO THE OORT CLOUD The Extended Scattered Disk Julio A. Fernández Departamento de Astronomía, Facultad de Ciencias, Montevideo, URUGUAY Adrián Brunini, Tabaré Gallardo, Rodney Gomes
More informationTisserand parameter (J) is more suitable for long-term orbital evolution. Following Carusi et al. (1987), Levison &
THE ASTRONOMICAL JOURNAL, 115:848È854, 1998 February ( 1998. The American Astronomical Society. All rights reserved. Printed in U.S.A. COLLISIONAL PROBABILITY OF PERIODIC COMETS WITH THE TERRESTRIAL PLANETS:
More informationPLANETARY MIGRATION AND PLUTINO ORBITAL INCLINATIONS R. S. GOMES
THE ASTRONOMICAL JOURNAL, 120:2695È2707, 2000 November ( 2000. The American Astronomical Society. All rights reserved. Printed in U.S.A. PLANETARY MIGRATION AND PLUTINO ORBITAL INCLINATIONS R. S. GOMES
More informationarxiv:astro-ph/ v3 2 Mar 2006
Mon. Not. R. Astron. Soc. 000, 1?? (2006) Printed 5 February 2008 (MN LATEX style file v2.2) Pluto s moon system: survey of the phase space I. I. Nagy, Á. Süli and B. Érdi Department of Astronomy, Eötvös
More informationEVOLUTIONS OF SMALL BODIES IN OUR SOLAR SYSTEM
EVOLUTIONS OF SMALL BODIES IN OUR SOLAR SYSTEM Dynamics and collisional processes Dr. Patrick MICHEL Côte d'azur Observatory Cassiopée Laboratory, CNRS Nice, France 1 Plan Chapter I: A few concepts on
More informationarxiv:astro-ph/ v1 7 Feb 2007
Is the Outer Solar System Chaotic? Wayne B. Hayes arxiv:astro-ph/0702179v1 7 Feb 2007 February 4, 2008 Wayne Hayes Computer Science Department University of California, Irvine Irvine, California 92697-3435
More informationSymplectic Correctors for Canonical Heliocentric N-Body Maps
Symplectic Correctors for Canonical Heliocentric N-Body Maps J. Wisdom Massachusetts Institute of Technology, Cambridge, MA 02139 wisdom@poincare.mit.edu Received ; accepted 2 ABSTRACT Symplectic correctors
More informationEXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS FOR EACH OF MULTIPLE BODIES
The Astrophysical Journal, 687:1294Y1302, 2008 November 10 # 2008. The American Astronomical Society. All rights reserved. Printed in U.S.A. EXTENDING NACOZY S APPROACH TO CORRECT ALL ORBITAL ELEMENTS
More informationThe stability of planets in the Alpha Centauri system
The stability of planets in the Alpha Centauri system Paul A. Wiegert 1, Matt Holman 2 1 Department of Astronomy, University of Toronto, Toronto, Canada 2 Canadian Institute for Theoretical Astrophysics,
More informationThe Mars 1:2 Resonant Population. Tabaré Gallardo.
1 The Mars 1:2 Resonant Population Tabaré Gallardo Departamento de Astronomía, Facultad de Ciencias, Iguá 4225, 11400 Montevideo, Uruguay gallardo@fisica.edu.uy Icarus, in press May 2007 2 Abstract An
More informationarxiv:astro-ph/ v1 16 Sep 1996
The stability of planets in the Alpha Centauri system Paul A. Wiegert 1, Matt Holman 2 1 Department of Astronomy, University of Toronto, Toronto, Canada 2 Canadian Institute for Theoretical Astrophysics,
More informationMOTIONS 1. INTRODUCTION 2. KOLMOGOROV-SMIRNOV REGULARIZED ORBITAL
THE ASTRONOMICAL JOURNAL, 120:3333È3339, 2000 December ( 2000. The American Astronomical Society. All rights reserved. Printed in U.S.A. LONG-TERM INTEGRATION ERROR OF KUSTAANHEIMO-STIEFEL REGULARIZED
More informationA survey of near-mean-motion resonances between Venus and Earth
Celest Mech Dyn Astr (2) 7:63 76 DOI 7/s569--9266-6 ORIGINAL ARTICLE A survey of near-mean-motion resonances between Venus and Earth Á. Bazsó R. Dvorak E. Pilat-Lohinger V. Eybl Ch. Lhotka Received: 3
More informationLong period comet encounters with the planets: an analytical approach. G. B. Valsecchi IAPS-INAF, Roma, Italy IFAC-CNR, Sesto Fiorentino, Italy
Long period comet encounters with the planets: an analytical approach G. B. Valsecchi IAPS-INAF, Roma, Italy IFAC-CNR, Sesto Fiorentino, Italy The engine of cometary orbital evolution The orbital evolution
More informationPlanetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune
PASJ: Publ. Astron. Soc. Japan 54, 493 499, 2002 June 25 c 2002. Astronomical Society of Japan. Planetary Perturbations on the 2 : 3 Mean Motion Resonance with Neptune Tetsuharu FUSE Subaru Telescope,
More informationAlignment indices: a new, simple method for determining the ordered or chaotic nature of orbits
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSA: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 34 (2001) 10029 10043 PII: S0305-4470(01)25097-5 Alignment indices: a new, simple method for determining
More informationarxiv: v2 [astro-ph.ep] 1 Feb 2010
Noname manuscript No. (will be inserted by the editor) Where are the Uranus Trojans? R. Dvorak Á. Bazsó L.-Y. Zhou arxiv:911.2366v2 [astro-ph.ep] 1 Feb 21 Received: date / Accepted: date Abstract The area
More informationResonance and chaos ASTRONOMY AND ASTROPHYSICS. II. Exterior resonances and asymmetric libration
Astron. Astrophys. 328, 399 48 (1997) ASTRONOMY AND ASTROPHYSICS Resonance and chaos II. Exterior resonances and asymmetric libration O.C. Winter 1,2 and C.D. Murray 2 1 Grupo de Dinâmica Orbital e Planetologia,
More informationChaos Indicators. C. Froeschlé, U. Parlitz, E. Lega, M. Guzzo, R. Barrio, P.M. Cincotta, C.M. Giordano, C. Skokos, T. Manos, Z. Sándor, N.
C. Froeschlé, U. Parlitz, E. Lega, M. Guzzo, R. Barrio, P.M. Cincotta, C.M. Giordano, C. Skokos, T. Manos, Z. Sándor, N. Maffione November 17 th 2016 Wolfgang Sakuler Introduction Major question in celestial
More informationResonance In the Solar System
Resonance In the Solar System Steve Bache UNC Wilmington Dept. of Physics and Physical Oceanography Advisor : Dr. Russ Herman Spring 2012 Goal numerically investigate the dynamics of the asteroid belt
More information1. INTRODUCTION. THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( The American Astronomical Society. All rights reserved. Printed in U.S.A.
THE ASTRONOMICAL JOURNAL, 121:1730È1735, 2001 March ( 2001. The American Astronomical Society. All rights reserved. Printed in U.S.A. ON THE SIZE DEPENDENCE OF THE INCLINATION DISTRIBUTION OF THE MAIN
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More informationHAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS. 1. Introduction
HAMILTONIAN STABILITY OF SPIN ORBIT RESONANCES IN CELESTIAL MECHANICS ALESSANDRA CELLETTI 1 and LUIGI CHIERCHIA 2 1 Dipartimento di Matematica, Università di Roma Tor Vergata, Via della Ricerca Scientifica,
More informationThe interplay of chaos between the terrestrial and giant planets
Mon. Not. R. Astron. Soc. (2010) doi:10.1111/j.1365-2966.2010.17027.x The interplay of chaos between the terrestrial and giant planets Wayne B. Hayes, 1,2 Anton V. Malykh 1 and Christopher M. Danforth
More informationThe precession of the perihelion of Mercury explained by Celestial Mechanics of Laplace Valdir Monteiro dos Santos Godoi
The precession of the perihelion of Mercury explained by Celestial Mechanics of Laplace Valdir Monteiro dos Santos Godoi valdir.msgodoi@gmail.com ABSTRACT We calculate in this article an exact theoretical
More informationChaos and Liapunov exponents
PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic
More informationTheory of mean motion resonances.
Theory of mean motion resonances. Mean motion resonances are ubiquitous in space. They can be found between planets and asteroids, planets and rings in gaseous disks or satellites and planetary rings.
More informationA REGION VOID OF IRREGULAR SATELLITES AROUND JUPITER
The Astronomical Journal, 136:909 918, 2008 September c 2008. The American Astronomical Society. All rights reserved. Printed in the U.S.A. doi:10.1088/0004-6256/136/3/909 A REGION VOID OF IRREGULAR SATELLITES
More informationCharacterization of Attractors in Gumowski-Mira Map Using Fast Lyapunov Indicators
Letter Forma, 21, 151 158, 2006 Characterization of Attractors in Gumowski-Mira Map Using Fast Lyapunov Indicators L. M. SAHA 1 *, M. K. DAS 2 and Mridula BUDHRAJA 3 1 Department of Mathematics, Zakir
More informationKepler, Newton, and laws of motion
Kepler, Newton, and laws of motion First: A Little History Geocentric vs. heliocentric model for solar system (sec. 2.2-2.4)! The only history in this course is this progression: Aristotle (~350 BC) Ptolemy
More informationChaotic diffusion of small bodies in the Solar System
Chaotic diffusion of small bodies in the Solar System Kleomenis Tsiganis (tsiganis@obs-nicefr) and Alessandro Morbidelli Observatoire de la Côte d Azur - Nice, CNRS, BP 4229, 06304 Nice Cedex 4, France
More informationarxiv: v1 [astro-ph.ep] 3 Apr 2018
Astronomy& Astrophysics manuscript no. zanardi_ c ESO 28 September 2, 28 The role of the general relativity on icy body reservoirs under the effects of an inner eccentric Jupiter M. Zanardi, 2, G. C. de
More informationLong-term integrations and stability of planetary orbits in our Solar system
Mon. Not. R. Astron. Soc. 336, 483 500 (2002) Long-term integrations and stability of planetary orbits in our Solar system Takashi Ito and Kiyotaka Tanikawa Data Analysis Center and Theoretical Astrophysics
More informationMean-Motion Resonance and Formation of Kirkwood Gaps
Yan Wang Project 1 PHYS 527 October 13, 2008 Mean-Motion Resonance and Formation of Kirkwood Gaps Introduction A histogram of the number of asteroids versus their distance from the Sun shows some distinct
More informationRESONANCE OVERLAP IS RESPONSIBLE FOR EJECTING PLANETS IN BINARY SYSTEMS
ACCEPTED FOR PUBLICATION IN APJ Preprint typeset using LATEX style emulateapj v. 11/26/04 RESONANCE OVERLAP IS RESPONSIBLE FOR EJECTING PLANETS IN BINARY SYSTEMS LAWRENCE R. MUDRYK AND YANQIN WU Department
More informationSurfing on the edge: chaos versus near-integrability in the system of Jovian planets
Mon. Not. R. Astron. Soc. 386, 295 306 (2008) doi:0./j.365-2966.2008.3024.x Surfing on the edge: chaos versus near-integrability in the system of Jovian planets Wayne B. Hayes Computer Science Department,
More informationResonancias orbitales
www.fisica.edu.uy/ gallardo Facultad de Ciencias Universidad de la República Uruguay Curso Dinámica Orbital 2018 Orbital motion Oscillating planetary orbits Orbital Resonances Commensurability between
More informationarxiv:astro-ph/ v1 29 Jul 2006
On a possible regularity connecting physical characteristics of a Solar system planet and elements of its orbit. arxiv:astro-ph/0607658v1 29 Jul 2006 Z.G.Murzahanov, V.P.Merezhin e-mail: dulkyn@mail.ru
More informationThe behavior of the AZ and QZ solutions in Figure 3 shows only dierences of
years; the latter is magnied by about a factor of 50 for better visualization. The reduced solution has a high and a low eccentricity state with regular and fast transitions between them. Its phase is
More informationUnit 1: The Earth in the Universe
Unit 1: The Earth in the Universe 1. The Universe 1.1. First ideas about the Universe 1.2. Components and origin 1.3. Sizes and distances 2. The Solar System 3. The planet Earth 3.1. Movements of the Earth
More informationInvestigation of evolution of orbits similar to that of (4179) Toutatis
Investigation of evolution of orbits similar to that of (4179) Toutatis Wm. Robert Johnston May 1994 Abstract This project compared the orbital evolution of asteroid (4179) Toutatis to that of several
More informationA study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1
Astronomy & Astrophysics manuscript no. Eris c ESO 2013 March 27, 2013 A study upon Eris I. Describing and characterizing the orbit of Eris around the Sun I. Breda 1 Faculty of Sciences (FCUP), University
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationLecture 13. Gravity in the Solar System
Lecture 13 Gravity in the Solar System Guiding Questions 1. How was the heliocentric model established? What are monumental steps in the history of the heliocentric model? 2. How do Kepler s three laws
More informationPaweł Kankiewicz, Ireneusz Włodarczyk
Meeting on Asteroids and Comets in Europe Vienna, Austria, 12-14.05.2006 The stability of orbits of effective Mars Crossers Paweł Kankiewicz, Ireneusz Włodarczyk Astrophysics Division, Institute of Physics,
More informationAstr 2320 Tues. Jan. 24, 2017 Today s Topics Review of Celestial Mechanics (Ch. 3)
Astr 2320 Tues. Jan. 24, 2017 Today s Topics Review of Celestial Mechanics (Ch. 3) Copernicus (empirical observations) Kepler (mathematical concepts) Galileo (application to Jupiter s moons) Newton (Gravity
More informationA STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3:1 RESONANCE. Érica C. Nogueira, Othon C. Winter
A STUDY OF CLOSE ENCOUNTERS BETWEEN MARS AND ASTEROIDS FROM THE 3: RESONANCE Érica C. Nogueira, Othon C. Winter Grupo de Dinâmica Orbital e Planetologia UNESP -- Guaratinguetá -- Brazil Antonio F.B. de
More informationDynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly
First Kepler s Law The secondary body moves in an elliptical orbit, with the primary body at the focus Valid for bound orbits with E < 0 The conservation of the total energy E yields a constant semi-major
More informationBROCK UNIVERSITY. 1. The observation that the intervals of time between two successive quarter phases of the Moon are very nearly equal implies that
BROCK UNIVERSITY Page 1 of 10 Test 1: November 2014 Number of pages: 10 Course: ASTR 1P01, Section 2 Number of students: 30 Examination date: 10 November 2014 Time limit: 50 min Time of Examination: 9:00
More information1. Solar System Overview
Astronomy 241: Foundations of Astrophysics I 1. Solar System Overview 0. Units and Precision 1. Constituents of the Solar System 2. Motions: Rotation and Revolution 3. Formation Scenario Units Text uses
More informationON THE DYNAMICAL STABILITY OF THE SOLAR SYSTEM
The Astrophysical Journal, 683:1207 1216, 2008 August 20 # 2008. The American Astronomical Society. All rights reserved. Printed in U.S.A. A ON THE DYNAMICAL STABILITY OF THE SOLAR SYSTEM Konstantin Batygin
More informationThe orbit evolution of 32 plutinos over 100 million year
A&A 368, 700 705 (2001) DOI: 10.1051/0004-6361:20010056 c ESO 2001 Astronomy & Astrophysics The orbit evolution of 32 plutinos over 100 million year X.-S. Wan 1 and T.-Y. Huang 1,2 1 Astronomy Department,
More informationDynamics of the conservative and dissipative. spin orbit problem
Dynamics of the conservative and dissipative spin orbit problem Alessandra Celletti Dipartimento di Matematica, Università di Roma Tor Vergata, Via della Ricerca Scientifica, I-0033 Roma (Italy), celletti@mat.uniroma2.it
More informationAccretion of Planets. Bill Hartmann. Star & Planet Formation Minicourse, U of T Astronomy Dept. Lecture 5 - Ed Thommes
Accretion of Planets Bill Hartmann Star & Planet Formation Minicourse, U of T Astronomy Dept. Lecture 5 - Ed Thommes Overview Start with planetesimals: km-size bodies, interactions are gravitational (formation
More informationGravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the
Gravity: Motivation An initial theory describing the nature of the gravitational force by Newton is a product of the resolution of the Geocentric-Heliocentric debate (Brahe s data and Kepler s analysis)
More informationDynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth.
Dynamical evolution of asteroid 2009 HE60, a quasi-satellite of Earth. P. Caspari ª ª Swinburne Astronomy Online, Student 6593658 Not yet received Not yet available online Abstract This paper announces
More informationChapter 29. The Solar System. The Solar System. Section 29.1 Models of the Solar System notes Models of the Solar System
The Solar System Chapter 29 The Solar System Section 29.1 Models of the Solar System 29.1 notes Models of the Solar System Geocentric: : Earth-centered model of the solar system. (Everything revolves around
More informationAST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements
AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends
More informationBROCK UNIVERSITY. 1. The observation that the intervals of time between two successive quarter phases of the Moon are very nearly equal implies that
BROCK UNIVERSITY Page 1 of 10 Test 1: November 2014 Number of pages: 10 Course: ASTR 1P01, Section 2 Number of students: 961 Examination date: 7 November 2014 Time limit: 50 min Time of Examination: 17:00
More informationarxiv:astro-ph/ v1 27 Sep 1995
THE PHASE SPACE STRUCTURE NEAR NEPTUNE RESONANCES IN THE KUIPER BELT arxiv:astro-ph/9509141v1 27 Sep 1995 RENU MALHOTRA Lunar and Planetary Institute, 3600 Bay Area Blvd., Houston, TX 77058. Electronic
More informationAstronomy 241: Review Questions #2 Distributed: November 7, 2013
Astronomy 241: Review Questions #2 Distributed: November 7, 2013 Review the questions below, and be prepared to discuss them in class. For each question, list (a) the general topic, and (b) the key laws
More informationAstronomy Section 2 Solar System Test
is really cool! 1. The diagram below shows one model of a portion of the universe. Astronomy Section 2 Solar System Test 4. Which arrangement of the Sun, the Moon, and Earth results in the highest high
More information+ (38 yr) 1 yr. = 742 mo. 1 yr The number of days in 742 months is
ASTR 101 Homework 2 Solutions 3-44 Chinese Calendar The traditional Chinese lunar calendar has 12 months in most years but adds a thirteenth month to 22 of every 60 years. How many days does this give
More informationThe effects of YORP on the spin rate distribution of the Near Earth Objects
The effects of YORP on the spin rate distribution of the Near Earth Objects A. Rossi, F. Marzari, D.J. Scheeres ISTI CNR, Pisa Università di Padova University of Colorado Binary Asteroid Dynamics Workshop
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationAstronomy 1001/1005 Midterm (200 points) Name:
Astronomy 1001/1005 Midterm (00 points) Name: Instructions: Mark your answers on this test AND your bubble sheet You will NOT get your bubble sheet back One page of notes and calculators are allowed Use
More informationas the orbits of distant planetoids are expected to be randomized over billions of year by the gravity of the four giant planets.
Dynamics of Extreme Outer Solar System Objects Principle investigators/supervisors: Denis Gaidashev (Matematiska Institution); Erik Zackrisson (Institutionen fo r fysik och astronomi) 1. Purpose and aims
More informationLong-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators
Long-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators by L. Anselmo and C. Pardini (Luciano.Anselmo@isti.cnr.it & Carmen.Pardini@isti.cnr.it)
More informationStability Dynamics Habitability of Planets in Binary Systems
Stability Dynamics Habitability of Planets in Binary Systems Facts and Questions: Binary and Multi-Star Systems in the Solar neighbourhood: >60% Do the gravitational perturbations allow planetary formation?
More informationPADEU. Pulsating zero velocity surfaces and capture in the elliptic restricted three-body problem. 1 Introduction
PADEU PADEU 15, 221 (2005) ISBN 963 463 557 c Published by the Astron. Dept. of the Eötvös Univ. Pulsating zero velocity surfaces and capture in the elliptic restricted three-body problem F. Szenkovits
More informationObservational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws
Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial
More informationGravitation and the Motion of the Planets
Gravitation and the Motion of the Planets 1 Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets go around
More informationDynamical Evolution of Ecliptic Comets
Duncan et al.: Dynamical Evolution of Ecliptic Comets 193 Dynamical Evolution of Ecliptic Comets Martin Duncan Queen s University Harold Levison and Luke Dones Southwest Research Institute Ecliptic comets
More informationSTABILITY OF HYPOTHETICAL TROJAN PLANETS IN EXOPLANETARY SYSTEMS
STABILITY OF HYPOTHETICAL TROJAN PLANETS IN EXOPLANETARY SYSTEMS Bálint Érdi, Georgina Fröhlich, Imre Nagy and Zsolt Sándor Department of Astronomy Loránd Eötvös University Pázmány Péter sétány 1/A H-1117
More informationWILLIAM I. NEWMAN, UCLA PHILIP W. SHARP, U. OF AUCKLAND BRUCE G. BILLS, JPL
CATASTROPHIC METEORITE 1 IMPACTS ON EARTH WILLIAM I. NEWMAN, UCLA PHILIP W. SHARP, U. OF AUCKLAND BRUCE G. BILLS, JPL Solar system with different scales for distance and for size: presents role of giant
More informationOn the numerical integration of variational equations
On the numerical integration of variational equations Haris Skokos Max Planck Institute for the Physics of Complex Systems Dresden, Germany E-mail: hskokos@pks.mpg.de, URL: http://www.pks.mpg.de/~hskokos/
More informationSolar System revised.notebook October 12, 2016 Solar Nebula Theory
Solar System revised.notebook The Solar System Solar Nebula Theory Solar Nebula was a rotating disk of dust and gas w/ a dense center dense center eventually becomes the sun start to condense b/c of gravity
More informationLecture Outlines. Chapter 15. Astronomy Today 7th Edition Chaisson/McMillan Pearson Education, Inc.
Lecture Outlines Chapter 15 Astronomy Today 7th Edition Chaisson/McMillan Chapter 15 The Formation of Planetary Systems Units of Chapter 15 15.1 Modeling Planet Formation 15.2 Terrestrial and Jovian Planets
More informationLinearized methods for ordinary di erential equations
Applied Mathematics and Computation 104 (1999) 109±19 www.elsevier.nl/locate/amc Linearized methods for ordinary di erential equations J.I. Ramos 1 Departamento de Lenguajes y Ciencias de la Computacion,
More informationMoon Obs #1 Due! Moon visible: early morning through afternoon. 6 more due June 13 th. 15 total due June 25 th. Final Report Due June 28th
Moon Obs #1 Due! Moon visible: early morning through afternoon 6 more due June 13 th 15 total due June 25 th Final Report Due June 28th Our Solar System Objectives Overview of what is in our solar system
More informationPubl. Astron. Obs. Belgrade No. 96 (2017), COMPUTATION OF TRANSIT ORBITS IN THE THREE-BODY-PROBLEM WITH FAST LYAPUNOV INDICATORS
Publ. Astron. Obs. Belgrade No. 96 (017), 71-78 Invited Lecture COMPUTATION OF TRANSIT ORBITS IN THE THREE-BODY-PROBLEM ITH FAST LYAPUNOV INDICATORS M. GUZZO 1 and E. LEGA 1 Università degli Studi di Padova,
More informationDynamically Unstable Planetary Systems Emerging Out of Gas Disks
EXTREME SOLAR SYSTEMS ASP Conference Series, Vol. 398, 2008 D. Fischer, F. A. Rasio, S. E. Thorsett, and A. Wolszczan Dynamically Unstable Planetary Systems Emerging Out of Gas Disks Soko Matsumura, Edward
More informationEarth Science Unit 6: Astronomy Period: Date: Elliptical Orbits
Earth Science Name: Unit 6: Astronomy Period: Date: Lab # 5 Elliptical Orbits Objective: To compare the shape of the earth s orbit (eccentricity) with the orbits of and with a circle. other planets Focus
More informationAnalysis of frozen orbits for solar sails
Trabalho apresentado no XXXV CNMAC, Natal-RN, 2014. Analysis of frozen orbits for solar sails J. P. S. Carvalho, R. Vilhena de Moraes, Instituto de Ciência e Tecnologia, UNIFESP, São José dos Campos -
More informationKozai-Lidov oscillations
Kozai-Lidov oscillations Kozai (1962 - asteroids); Lidov (1962 - artificial satellites) arise most simply in restricted three-body problem (two massive bodies on a Kepler orbit + a test particle) e.g.,
More informationLocated Between The Orbits Of
The Asteroids Shown In The Data Table Are Located Between The Orbits Of The asteroid belt lies between the orbits of In science fiction movies, spaceships are often shown dodging through large numbers
More informationTides and Lagrange Points
Ast111, Lecture 3a Tides and Lagrange Points Arial view of Tidal surge in Alaska (S. Sharnoff) Tides Tidal disruption Lagrange Points Tadpole Orbits and Trojans Tidal Bulges Tides Tidal Force The force
More informationSOLAR SYSTEM 2019 SAMPLE EXAM
SOLAR SYSTEM 2019 SAMPLE EXAM Team Name: Team #: No calculators are allowed. All questions are of equal weight unless otherwise noted. Turn in all materials when you have completed the test! Make sure
More informationTransneptunian objects. Minor bodies in the outer Solar System. Transneptunian objects
Transneptunian objects Minor bodies in the outer Solar System Planets and Astrobiology (2016-2017) G. Vladilo Around 1980 it was proposed that the hypothetical disk of small bodies beyond Neptune (called
More informationJoseph Castro Mentor: Nader Haghighipour
ON THE POSSIBILITY OF ADDITIONAL PLANETS IN THE γ CEPHEI BINARY-PLANETARY SYSTEM Joseph Castro Mentor: Nader Haghighipour ABSTRACT Results of the simulations of the dynamical stability of additional hypothetical
More informationAstronomy 241: Foundations of Astrophysics I. The Solar System
Astronomy 241: Foundations of Astrophysics I. The Solar System Astronomy 241 is the first part of a year-long introduction to astrophysics. It uses basic classical mechanics and thermodynamics to analyze
More information