Análisis comparativo de la estructura orbital de modelos de galaxias triaxiales no cuspidales y cuspidales utilizando distintas expansiones para el potencial gravitatorio

Abstract

Un método para generar modelos de galaxias elípticas consiste en partir de un sistema de partículas con una configuración adecuada y seguir su evolución con un código de N-cuerpos hasta que alcanza el equilibrio. En trabajos anteriores hemos utilizado el código de L. A. Aguilar para generar modelos tanto no cuspidales como cuspidales, ofreciendo los segundos una mejor representación de las galaxias elípticas. Sin embargo, el método propuesto por Hernquist & Ostriker es más adecuado para los casos con cúspides centrales, puesto que utiliza una expansión del potencial en funciones radiales que ya tiene carácter cuspidal. Aquí consideramos dos modelos, uno no cuspidal y otro cuspidal, obtenidos con el método de Aguilar y cuyos potenciales ajustamos ahora con la expansión de Hernquist & Ostriker. Obtuvimos las ecuaciones variacionales correspondientes y las empleamos para determinar los exponentes de Lyapunov de órbitas ya estudiadas con el método de Aguilar y separar las caóticas. Las órbitas regulares las clasificamos mediante el análisis de sus frecuencias orbitales. Presentamos una comparación de los resultados obtenidos con los dos métodos para ambos tipos de órbitas.

One way to create models of elliptical galaxies is to start from a system of particles with a suitable distribution and to follow its evolution with an N-body code until it reaches equilibrium. In previous works we used the code of L.A. Aguilar to generate noncuspy and cuspy models, the latter giving a better representation of elliptical galaxies. However, the method proposed by Hernquist and Ostriker is more suitable for the cases with central cusps, because it uses an expansion of the potential in radial functions that already reflects the cuspy character. Here we consider two models, one non-cuspy and another cuspy, obtained with the method of Aguilar and whose potentials we now fit with the Hernquist and Ostriker expansion. We obtained the corresponding variational equations and we used them to determine the Lyapunov exponents of orbits already investigated with the method of Aguilar and to separate the chaotic ones. We classified the regular orbits through the analysis of their orbital frequencies. Here we present a comparison of the results obtained with the two methods for both types of orbits.