We used a multipolar code to create, through the dissipationless collapses of systems of 1,000,000 particles, three self-consistent triaxial stellar systems with axial ratios corresponding to those of E4, E5 and E6 galaxies. The E5 and E6 models have small, but significant, rotational velocities although their total angular momenta are zero, that is, they exhibit figure rotation; the rotational velocity decreases with decreasing flattening of the models and for the E4 model it is essentially zero. Except for minor changes, probably caused by unavoidable relaxation effects, the systems are highly stable. The potential of each system was subsequently approximated with interpolating formulae yielding smooth potentials, stationary for the non-rotating model and stationary in the rotating frame for the rotating ones. The Lyapunov exponents could then be computed for randomly selected samples of the bodies that make up the different systems, allowing the recognition of regular and partially and fully chaotic orbits. Finally, the regular orbits were Fourier analyzed and classified using their locations on the frequency map. As it could be expected, the percentages of chaotic orbits increase with the flattening of the system. As one goes from E6 through E4, the fraction of partially chaotic orbits relative to that of fully chaotic ones increases, with the former surpassing the latter in model E4; the likely cause of this behavior is that 123 308 R. O. Aquilano et al. triaxiality diminishes from E6 through E4, the latter system being almost axially symmetric. We especulate that some of the partially chaotic orbits may obey a global integral akin to the long axis component of angular momentum. Our results show that is perfectly possible to have highly stable triaxial models with large fractions of chaotic orbits, but such systems cannot have constant axial ratios from center to border: a slightly flattened reservoir of highly chaotic orbits seems to be mandatory for those systems.

We used a multipolar code to create, through dissipationless collapses of systems of 10 6 particles, two cuspy self-consistent triaxial stellar systems with γ ≈ 1. One of the systems has an axial ratio similar to that of an E4 galaxy and it is only mildly triaxial (T = 0.914), while the other one is strongly triaxial (T = 0.593) and its axial ratio lies in between those of Hubble types E5 and E6. Both models rotate although their total angular momenta are zero, i.e., they exhibit figure rotation. The angular velocity is very small for the less triaxial model and, while it is larger for the more triaxial one, it is still comparable to that found by Muzzio (Celest Mech Dynam Astron 96 (2): [85][86][87][88][89][90][91][92][93][94][95][96][97] 2006) to affect only slightly the dynamics of a similar model. Except for minor evolution, probably caused by unavoidable relaxation effects of the N-body code, the systems are highly stable. The potential of each system was subsequently approximated with interpolating formulae yielding smooth potentials, stationary in frames that rotate with the models. The Lyapunov exponents could then be computed for randomly selected samples of the bodies that make up the two systems, allowing the recognition of regular and of partially and fully chaotic orbits. Finally, the regular orbits were Fourier analyzed and classified using their locations on the frequency map. Most of the orbits are chaotic, and by a wide margin: less than 30% of the orbits are regular in our most triaxial model. Regular orbits are dominated by tubes, long axis ones in the less 123 380 J. C. Muzzio et al. triaxial model and short axis tubes in the more triaxial one. Most of the boxes are resonant (i.e., they are boxlets), as could be expected from cuspy systems.

We used the N‐body code of Hernquist & Ostriker to build a dozen cuspy (γ≃ 1) triaxial models of stellar systems through dissipationless collapses of initially spherical distributions of 106 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 5000 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 in building such models with the method of Schwarzschild should be attributed to the method itself and not to physical causes.