The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In par- ticular, we consider the relative Poincare variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator Mean Exponential Growth factor of Nearby Orbits (MEGNO), a measure of the maximal Lyapunov exponent, that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best-fitting solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best-fitting solutions are found in dynamically active region of the phase space. The long-term stability of the system is determined by a net of low-order two-body and three-body mean mo- tion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myr of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behaviour.

[1]  S. Ferraz-Mello,et al.  Resonant Structure of the Outer Solar System in the Neighborhood of the Planets , 2001 .

[2]  Kevin P. Rauch,et al.  Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions , 1999 .

[3]  D. Nesvorný,et al.  Frequency modified fourier transform and its application to asteroids , 1996 .

[4]  P. Cincotta.  Arnold diffusion: an overview through dynamical astronomy , 2002 .

[5]  J. Lissauer.  Chaotic motion in the Solar System , 1999 .

[6]  M. Suzuki,et al.  General theory of fractal path integrals with applications to many‐body theories and statistical physics , 1991 .

[7]  Three-Body Mean Motion Resonances and the Chaotic Structure of the Asteroid Belt , 1998 .

[8]  S. Ferraz-Mello,et al.  Regular motions in extra-solar planetary systems , 2004, astro-ph/0402335.

[9]  Edmund Taylor Whittaker,et al.  A Treatise on the Analytical Dynamics of Particles and Rigid Bodies: INDEX OF TERMS EMPLOYED , 1988 .

[10]  Archie E. Roy,et al.  Predictability, stability, and chaos in N-body dynamical systems , 1991 .

[11]  R. Paul Butler,et al.  Two New Planets from the Anglo-Australian Planet Search , 2001 .

[12]  J. Laskar,et al.  High order symplectic integrators for perturbed Hamiltonian systems , 2000 .

[13]  Alessandro Morbidelli,et al.  Modern celestial mechanics : aspects of solar system dynamics , 2002 .

[14]  J. Chambers.  A hybrid symplectic integrator that permits close encounters between massive bodies , 1999 .

[15]  R. Devaney.  Celestial mechanics. , 1979, Science.

[16]  Gianfranco Capriz,et al.  On the basic laws of reacting mixtures of structured continua , 1974 .

[17]  Victor Brumberg,et al.  Essential Relativistic Celestial Mechanics , 1991 .

[18]  M. Konacki,et al.  Dynamical Properties of the Multiplanet System around HD 169830 , 2004 .

[19]  C. Froeschlé,et al.  The Lyapunov characteristic exponents-applications to celestial mechanics , 1984 .

[20]  G. Benettin,et al.  Lyapunov Characteristic Exponents for smooth dynamical systems and for hamiltonian systems; a method for computing all of them. Part 1: Theory , 1980 .

[21]  Victor Szebehely,et al.  Review of concepts of stability , 1984 .

[22]  The web of three-planet resonances in the outer Solar System: II. A source of orbital instability for Uranus and Neptune , 2006 .

[23]  P. M. Cincotta,et al.  Simple tools to study global dynamics in non-axisymmetric galactic potentials – I , 2000 .

[24]  K. Goździewski.  A dynamical analysis of the HD 37124 planetary system , 2003 .

[25]  F. Pirani.  MATHEMATICAL METHODS OF CLASSICAL MECHANICS (Graduate Texts in Mathematics, 60) , 1982 .

[26]  Maciej Konacki,et al.  Orbital configurations and dynamical stability of multiplanet systems around sun-like stars HD 202206, 14 Herculis, HD 37124, and HD 108874 , 2006 .

[27]  R. Duncombe.  Dynamics of the Solar System , 1979 .

[28]  A. Morbidelli,et al.  An Analytic Model of Three-Body Mean Motion Resonances , 1998 .

[29]  Harold F. Levison,et al.  A Multiple Time Step Symplectic Algorithm for Integrating Close Encounters , 1998 .

[30]  A. Deprit.  Elimination of the nodes in problems ofn bodies , 1983 .

[31]  R. Paul Butler,et al.  Five New Multicomponent Planetary Systems , 2005 .

[32]  Symplectic Tangent map for Planetary Motions , 1999 .

[33]  J. Laskar.  Analytical Framework in Poincare Variables for the Motion of the Solar System , 1991 .

[34]  Henri Poincaré,et al.  Leçons de mécanique céleste , 1905 .

[35]  Carles Simó,et al.  Phase space structure of multi-dimensional systems by means of the mean exponential growth factor of nearby orbits , 2003 .

[36]  Frequency Analysis of a Dynamical System , 1993 .

[37]  M. A. Hendry,et al.  Chaotic Worlds: from Order to Disorder in Gravitational N-Body Dynamical Systems , 2006 .

[38]  B. Melendo,et al.  Synchronous motion in the Kinoshita problem - Application to satellites and binary asteroids , 2005 .

[39]  S. Ferraz-Mello,et al.  Resonances in the motion of planets, satellites and asteroids , 1985 .

[40]  J. Wisdom,et al.  Symplectic maps for the N-body problem. , 1991 .

[41]  On the Origin of Chaos in the Asteroid Belt , 1998 .