Astron. Astrophys. Suppl. Ser. 147, 205-228

Simple tools to study global dynamics in non-axisymmetric galactic potentials - I

P.M. Cincotta¹ - C. Simó²

Send offprint request: P.M. Cincotta

1 - Facultad de Ciencias Astronómicas y Geofísicas, Universidad Nacional de La Plata, Paseo del Bosque, 1900 La Plata, Argentina
2 - Departament de Matemàtica Aplicada i Anàlisi, Universitat de Barcelona, Gran Via 585, 08007 Barcelona, Spain

Received July 30, 1999; accepted September 4, 2000

Abstract:

In a first part we discuss the well-known problem of the motion of a star in a general non-axisymmetric 2D galactic potential by means of a very simple but almost universal system: the pendulum model. It is shown that both loop and box families of orbits arise as a natural consequence of the dynamics of the pendulum. An approximate invariant of motion is derived. A critical value of the latter sharply separates the domains of loops and boxes and a very simple computation allows to get a clear picture of the distribution of orbits on a given energy surface. Besides, a geometrical representation of the global phase space using the natural surface of section for the problem, the 2D sphere, is presented. This provides a better visualization of the dynamics.

In a second part we introduce a new indicator of the basic dynamics, the Mean Exponential Growth factor of Nearby Orbits (MEGNO), that is suitable to investigate the phase space structure associated to a general Hamiltonian. When applied to the 2D logarithmic potential it is shown to be effective to obtain a picture of the global dynamics and, also, to derive good estimates of the largest Lyapunov characteristic number in realistic physical times. Comparisons with other techniques reveal that the MEGNO provides more information about the dynamics in the phase space than other wide used tools.

Finally, we discuss the structure of the phase space associated to the 2D logarithmic potential for several values of the semiaxis ratio and energy. We focus our attention on the stability analysis of the principal periodic orbits and on the chaotic component. We obtain critical energy values for which connections between the main stochastic zones take place. In any case, the whole chaotic domain appears to be always confined to narrow filaments, with a Lyapunov time about three characteristic periods.

Key words: galaxies: dynamics -- stellar dynamics -- methods: analytical-numerical -- chaos