Spatial structure of regular and chaotic orbits in a self-consistent triaxial stellar system (Q1776102)

There are two base approaches for numerical investigation of self- consistent models of stellar systems. The first one was developed by \textit{M. Schwarzschild} [A numerical model for a triaxial stellar system in dynamical equilibriu. Astrophys. J. 232, 236--247 (1979), doi:10.1086/157282] and consists of: 1. choosing a mass distribution, 2. finding the potential produced by that distribution , 3. determining the orbits sustained by that potential, 4. finding in which proportion those orbits should be added to yield the initial mass distribution. The second approach was introduced by \textit{L. S. Sparke} and \textit{J. A. Sellwood} [Dissection of an \(N\)-body bar. Mon. Not. R. Acad. Sci. 225, 653--675 (1987)]. They employed \(N\)-body codes to create stationary stable systems, then froze and smoothed the potential and, finally, they investigated the orbits of samples of the bodies whose distributions produce those potentials. The authors follow the second approach and in their previous papers investigated the spatial distribution of regular and chaotic stellar orbits inside self-consistent model of elliptical galaxies and galactic satellites. Here, the authors extend that model to investigate a self-consistent triaxial model and, especially, to recognize the different kinds of regular orbits.