Models of cuspy triaxial stellar systems. I. Stability and chaoticity
arXiv:1204.5428 · doi:10.1111/j.1365-2966.2012.21023.x
Abstract
We used the N-body code of Hernquist and Ostriker (1992) to build a dozen cuspy (γ\approx 1) triaxial models of stellar systems through dissipationless collapses of initially spherical distributions of 10^6 particles. We chose four sets of initial conditions that resulted in models morphologically resembling E2, E3, E4 and E5 galaxies, respectively. Within each set, three different seed numbers were selected for the random number generator used to create the initial conditions, so that the three models of each set are statistically equivalent. We checked the stability of our models using the values of their central densities and of their moments of inertia, which turned out to be very constant indeed. The changes of those values were all less than 3 per cent over one Hubble time and, moreover, we show that the most likely cause of those changes are relaxation effects in the numerical code. We computed the six Lyapunov exponents of nearly 5,000 orbits in each model in order to recognize regular, partially and fully chaotic orbits. All the models turned out to be highly chaotic, with less than 25 per cent of their orbits being regular. We conclude that it is quite possible to obtain cuspy triaxial stellar models that contain large fractions of chaotic orbits and are highly stable. The difficulty to build such models with the method of Schwarzschild (1979) should be attributed to the method itself and not to physical causes.
10 pages, 5 figures, published in MNRAS (The definitive version is available at www.blackwell-synergy.com)
References in corpus (2)
Cited by in corpus (11)
- A new code for orbit analysis and Schwarzschild modelling of triaxial stellar systems
- Milking the spherical cow: on aspherical dynamics in spherical coordinates
- On the relevance of chaos for halo stars in the Solar Neighbourhood
- Self-consistent triaxial models
- Discreteness effects, body chaos and the onset of radial-orbit instability
- Models of cuspy triaxial stellar systems. III: The effect of velocity anisotropy on chaoticity
- Partially chaotic orbits in a perturbed cubic force model
- Chaotic orbits obeying one isolating integral in a four dimensional map
- Models of cuspy triaxial stellar systems. IV: Rotating systems
- The Lyapunov exponents and the neighbourhood of periodic orbits
- Models of cuspy triaxial stellar systems. II. Regular orbits