A Counterpart to the Radial Orbit Instability in Triaxial Stellar Systems

Self-consistent solutions for triaxial mass models are highly non-unique. In general, some of these solutions might be dynamically unstable, making them inappropriate as descriptions of steady-state galaxies. Here we demonstrate for the first time the existence in triaxial galaxy models of an instability similar to the radial-orbit instability of spherical models. The instability manifests itself when the number of box orbits, with predominantly radially motions, is sufficiently large. N-body simulations verify that the evolution is due neither to chaotic orbits nor to departures of the model from self-consistency, but rather to a collective mode. The instability transforms the triaxial model into a more prolate, but still triaxial, configuration. Stable triaxial models are obtained when the mass contribution of radial orbits is reduced. The implications of our results for the shapes of dark-matter halos are discussed.

💡 Research Summary

The paper addresses the long‑standing problem of dynamical stability in self‑consistent triaxial galaxy models. While spherical systems are known to suffer from the radial‑orbit instability (ROI) when a large fraction of stars move on nearly radial orbits, it has been unclear whether an analogous collective instability can arise in fully three‑dimensional, triaxial configurations. Using the Schwarzschild orbit‑superposition technique, the authors construct a suite of triaxial mass models that differ primarily in the relative weight of box orbits – the orbits that support the central regions and are predominantly radial – versus tube (loop) orbits.

High‑resolution N‑body simulations are then employed to evolve each model forward in time. When the mass contribution of box orbits exceeds a critical threshold (roughly 20‑30 % of the total mass in the examples studied), the system undergoes a rapid, coherent deformation: the intermediate axis shortens while the major axis lengthens, producing a more prolate but still triaxial shape. This transformation is not driven by an increase in chaotic orbits, nor by a violation of the initial self‑consistency; Lyapunov analyses show that the majority of orbits remain regular, and the potential continues to satisfy the Schwarzschild constraints throughout the evolution. Instead, the authors identify a genuine collective mode – a global, low‑frequency oscillation of the mass distribution – that grows exponentially and saturates in the new configuration.

Conversely, when the box‑orbit fraction is reduced below the identified threshold, the models remain essentially unchanged over many dynamical times. Their axis ratios, velocity dispersion profiles, and orbital families stay close to the initial conditions, confirming that the instability is suppressed. The authors therefore propose a practical stability criterion for triaxial systems: the radial (box) orbit mass fraction must be limited to roughly one‑fifth of the total mass to avoid the ROI‑like mode.

The astrophysical implications are significant. Dark‑matter halos in cosmological ΛCDM simulations are typically triaxial, yet observational inferences (e.g., from weak lensing or stellar streams) often suggest rounder shapes. The present work suggests that if a halo’s inner region contains too many radial orbits, a collective instability could reshape it toward a more prolate geometry, potentially reconciling simulation predictions with observations. Moreover, the findings provide a diagnostic for constructing stable initial conditions in N‑body experiments that include baryonic components, gas inflows, or external tidal fields.

In summary, the study demonstrates for the first time that triaxial stellar systems possess an instability directly analogous to the spherical ROI, driven by an excess of box (radial) orbits. By quantifying the critical orbital composition and showing that the resulting mode is a genuine collective phenomenon, the authors offer both a deeper theoretical understanding of galaxy dynamics and practical guidance for modeling realistic, long‑lived triaxial structures.