Radial Orbit Instability Research Articles

Stirring N-body systems – II. Necessary conditions for the dark matter attractor

We study the evolution of the phase-space of collisionless N-body systems under repeated stirrings or perturbations, which has been shown to lead to a convergence towards a limited group of end states. This so-called attractor was previously shown to be independent of the initial system and environmental conditions. However the fundamental reason for its appearance is still unclear. It has been suggested that the origin of the attractor may be either radial infall (RI), the radial orbit instability (ROI), or energy exchange which, for instance, happens during violent relaxation. Here we examine the effects of a set of controlled perturbations, referred to as `kicks', which act in addition to the standard collisionless dynamics by allowing pre-specified instantaneous perturbations in phase-space. We first demonstrate that the attractor persists in the absence of RI and ROI by forcing the system to expand. We then consider radial velocity kicks in a rigid potential and isotropic velocity kicks, since there are no energy exchanges in these two recipes of kicks. We find that these kicks do not lead to the attractor, indicating that the energy exchange in a dynamic potential is indeed the physical mechanism responsible for the attractor.

Radially anisotropic models of collisionless spherical stellar systems without central singularity

A new class of the simplest equilibrium two-parameter distribution functions for spherical stellar systems with a radially anisotropic stellar velocity distribution is investigated. The models under consideration are a less singular counterpart of the so-called generalized polytropes, which in the past were among the most popular models in works on the equilibrium and stability of gravitating systems. In contrast to the well-known generalized polytropes, the proposed models have finite density and potential at the center. The absence of a singularity is necessary for a proper consideration of the radial orbit instability, which is the most important instability of spherical stellar systems. The main observed parameters of the proposed models (potential, density, anisotropy) are compared with those in well-known equilibrium models.

Equilibrium models of radially anisotropic spherical stellar systems with softened central potentials

We study a new class of equilibrium two-parametric distribution functions of spherical stellar systems with radially anisotropic velocity distribution of stars. The models are less singular counterparts of the so called generalized polytropes, widely used in works on equilibrium and stability of gravitating systems in the past. The offered models, unlike the generalized polytropes, have finite density and potential in the center. The absence of the singularity is necessary for proper consideration of the radial orbit instability, which is the most important instability in spherical stellar systems. Comparison of the main observed parameters (potential, density, anisotropy) predicted by the present models and other popular equilibrium models is provided.

StirringN-body systems: universality of end states

We study the evolution of the phase-space of collisionless N-body systems under repeated stirrings or perturbations. We find convergence towards a limited solution group, in accordance with Hansen 2010, that is independent of the initial system and environmental conditions, paying particular attention to the assumed gravitational paradigm (Newtonian and MOND). We examine the effects of changes to the perturbation scheme and in doing so identify a large group of perturbations featuring radial orbit instability (ROI) which always lead to convergence. The attractor is thus found to be a robust and reproducible effect under a variety of circumstances.

On the threshold anisotropy for unstable spherical stellar systems

The radial orbit instability generally arises in anisotropic collisionless stellar systems with the dominance of radial motions over transverse ones. Using the simplest anisotropic generalization of polytrope models for spherical clusters as an example, we show that the instability growth rates become exponentially small as the isotropic limit is approached. Given the finite lifetime of real astronomical objects, these systems should be assumed to become stable at some finite radial anisotropy.

Radial Orbit Instability: Review and Perspectives

This article presents elements about the radial orbit instability, which occurs in spherical self-gravitating systems with a strong anisotropy in the radial velocity direction. It contains an overview on the history of radial orbit instability. We also present the symplectic method we use to explore stability of equilibrium states directly related to the dissipation induced instability mechanism well known in theoretical mechanics and plasma physics.

Chaotic mixing and the secular evolution of triaxial cuspy galaxy models built with Schwarzschild’s method

We use both N-body simulations and integration in fixed potentials to explore the stability and the long-term secular evolution of self-consistent, equilibrium, non-rotating, triaxial spheroidal galactic models. More specifically, we consider Dehnen models built with the Schwarzschild method. We show that short-term stability depends on the degree of velocity anisotropy (radially anisotropic models are subject to rapid development of radial-orbit instability). Long-term stability, on the other hand, depends mainly on the properties of the potential, and in particular, on whether it admits a substantial fraction of strongly chaotic orbits. We show that in the case of a weak density cusp (gamma=1 Dehnen model) the N-body model is remarkably stable, while the strong-cusp (gamma=2) model exhibits substantial evolution of shape away from triaxiality, which we attribute to the effect of chaotic diffusion of orbits. The different behaviour of these two cases originates from the different phase space structure of the potential; in the weak-cusp case there exist numerous resonant orbit families that impede chaotic diffusion. We also find that it is hardly possible to affect the rate of this evolution by altering the fraction of chaotic orbits in the Schwarzschild model, which is explained by the fact that the chaotic properties of an orbit are not preserved by the N-body evolution. There are, however, parameters in Schwarzschild modelling that do affect the stability of an N-body model, so we discuss the recipes how to build a `good' Schwarzschild model.

Radial-orbit instability in modified Newtonian dynamics

The stability of radially anisotropic spherical stellar systems in modified Newtonian dynamics (MOND) is explored by means of numerical simulations performed with the N-body code N-MODY. We find that Osipkov-Merritt MOND models require for stability larger minimum anisotropy radius than equivalent Newtonian systems (ENSs) with dark matter, and also than purely baryonic Newtonian models with the same density profile. The maximum value for stability of the Fridman-Polyachenko-Shukhman parameter in MOND models is lower than in ENSs, but higher than in Newtonian models with no dark matter. We conclude that MOND systems are substantially more prone to radial-orbit instability than ENSs with dark matter, while they are able to support a larger amount of kinetic energy stored in radial orbits than purely baryonic Newtonian systems. An explanation of these results is attempted, and their relevance to the MOND interpretation of the observed kinematics of globular clusters, dwarf spheroidal and elliptical galaxies is briefly discussed.

Non-spherical similarity solutions for dark halo formation

We carry out fully 3-dimensional simulations of evolution from self-similar, spherically symmetric linear perturbations of a Cold Dark Matter dominated Einstein-de Sitter universe. As a result of the radial orbit instability, the haloes which grow from such initial conditions are triaxial with major-to-minor axis ratios of order 3:1. They nevertheless grow approximately self-similarly in time. In all cases they have power-law density profiles and near-constant velocity anisotropy in their inner regions. Both the power-law index and the value of the velocity anisotropy depend on the similarity index of the initial conditions, the former as expected from simple scaling arguments. Halo structure is thus not "universal" but remembers the initial conditions. On larger scales the density and anisotropy profiles show two characteristic scales, corresponding to particles at first pericentre and at first apocentre after infall. They are well approximated by the NFW model only for one value of the similarity index. In contrast, at all radii within the outer caustic the pseudo phase-space density can be fit by a single power law with an index which depends only very weakly on the similarity index of the initial conditions. This behaviour is very similar to that found for haloes formed from LCDM initial conditions and so can be considered approximately universal.

Radial orbit instability as a dissipation-induced phenomenon

This paper is devoted to radial orbit instability in the context of self-gravitating dynamical systems. We present this instability in the new frame of dissipation-induced instability theory. This allows us to obtain a rather simple proof based on energetics arguments and to clarify the associated physical mechanism.

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.

Caustics in growing cold dark matter haloes

We simulate the growth of isolated dark matter haloes from self-similar and spherically symmetric initial conditions. Our N-body code integrates the geodesic deviation equation in order to track the streams and caustics associated with individual simulation particles. The radial orbit instability causes our haloes to develop major-to-minor axis ratios approaching 10 to 1 in their inner regions. They grow similarly in time and have similar density profiles to the spherical similarity solution, but their detailed structure is very different. The higher dimensionality of the orbits causes their stream and caustic densities to drop much more rapidly than in the similarity solution. This results in a corresponding increase in the number of streams at each point. At 1% of the turnaround radius (corresponding roughly to the Sun's position in the Milky Way) we find of order 10^6 streams in our simulations, as compared to 10^2 in the similarity solution. The number of caustics in the inner halo increases by a factor of several, because a typical orbit has six turning points rather than one, but caustic densities drop by a much larger factor. This reduces the caustic contribution to the annihilation radiation. For the region between 1% and 50% of the turnaround radius, this is 4% of the total in our simulated haloes, as compared to 6.5% in the similarity solution. Caustics contribute much less at smaller radii. These numbers assume a 100 GeV c^-2 neutralino with present-day velocity dispersion 0.03 cm s^-1, but reducing the dispersion by ten orders of magnitude only doubles the caustic luminosity. We conclude that caustics will be unobservable in the inner parts of haloes. Only the outermost caustic might potentially be detectable.

THE RADIAL ORBIT INSTABILITY IN COLLISIONLESSN-BODY SIMULATIONS

Using a suite of self-gravitating, collisionless N-body models, we systematically explore a parameter space relevant to the onset and behavior of the radial orbit instability (ROI), whose strength is measured by the systemic axis ratios of the models. We show that a combination of two initial conditions, namely the velocity anisotropy and the virial ratio, determines whether a system will undergo ROI and exactly how triaxial the system will become. A third initial condition, the radial shape of the density profile, plays a smaller, but noticeable role. Regarding the dynamical development of the ROI, the instability a) begins after systems collapse to their most compact configuration and b) evolves fastest when a majority of the particles have radially anisotropic orbits while there is a lack of centrally-concentrated isotropic orbits. We argue that this is further evidence that self-reinforcing torques are the key to the onset of the ROI. Our findings support the idea that a separate orbit instability plays a role in halting the ROI.

NATURE AND NURTURE IN DARK MATTER HALOS

Cosmological simulations consistently predict specific properties of dark matter halos, but these have not yet led to a physical understanding that is generally accepted. This is especially true for the central regions of these structures. Recently, two major themes have emerged. In one, the dark matter halo is primarily a result of the sequential accretion of the primordial structure (i.e., nature), while in the other, dynamical relaxation (i.e., nurture) dominates at least in the central regions. Some relaxation is however required in either mechanism. In this paper, we accept the recently established scale-free substructure of halos as an essential part of both mechanisms. Consequently, a simple model for the central relaxation, based on a self-similar cascade of tidal interactions, is contrasted with a model based on the accretion of an adiabatically self-similar, primordial structure. We conclude that a weak form of this relaxation is present in the simulations, but that it is normally described as the radial orbit instability.

An instability in triaxial stellar systems

AbstractThe radial‐orbit instability is a collective phenomenon that has heretofore only been observed in spherical systems. We find that this instability occurs also in triaxial systems, as we checked by performing extensive N ‐body simulations whose initial conditions were obtained by sampling a self‐consistent triaxial model of a cuspy galaxy composed of luminous and dark matter (Capuzzo‐Dolcetta et al. 2007). N ‐body simulations show a time evolution of the galaxy that is not due to the development of chaotic motions but, rather, to the collective instability induced by an excess of radially biased (box) orbits. The instability quickly transforms such model into a more prolate configuration, with 0.64 < b/a < 0.77 and 0.6 < c/a < 0.7 for the dark halo and 0.64 < b/a < 0.77 and 0.59 < c/a < 0.67 for the luminous matter. Stable triaxial, cuspy galaxies with dark matter halo are obtained when the mass contribution of radially biased orbits to the solution is reduced. These results constitute the first evidence of radial orbit instability in triaxial stellar systems. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

The Role of the Radial Orbit Instability in Dark Matter Halo Formation and Structure

For a decade, N-body simulations have revealed a nearly universal dark matter density profile, which appears to be robust to changes in the overall density of the universe and the underlying power spectrum. Despite its universality, the physical origin of this profile has not yet been well understood. Semianalytic models by Barnes and coworkers have suggested that the density structure of dark matter halos is determined by the onset of the radial orbit instability (ROI). We have tested this hypothesis using N-body simulations of collapsing dark matter halos with a variety of initial conditions. For dynamically cold initial conditions, the resulting halo structures are triaxial in shape, due to the mild aspect of the instability. We examine how variations in initial velocity dispersion affect the onset of the instability, and find that an isotropic velocity dispersion can suppress the ROI entirely, while a purely radial dispersion does not. The quantity -->σ2/v2c is a criterion for instability, where regions with -->σ2/v2c 1 become triaxial due to the ROI or other perturbations. We also find that the radial orbit instability sets a scale length at which the velocity dispersion changes rapidly from isotropic to radially anisotropic. This scale length is proportional to the radius at which the density profile changes shape, as is the case in the semianalytic models; however, the coefficient of proportionality is different by a factor of ~2.5. We conclude that the radial orbit instability is likely to be a key physical mechanism responsible for the nearly universal profiles of simulated dark matter halos.

Retrograde precession of stellar orbits and gravitational loss-cone instability in galactic centers

We consider disk and spherical subsystems of stars with nearly radial orbits under conditions when the well-known radial orbit instability is not possible. This requires that the precession of stellar orbits be retrograde, i.e., in the direction opposite to the orbital rotation of stars. We show that an instability that is an analogue of the loss-cone instability known in plasma physics can then develop in the presence of a “loss cone” in the angular momentum distribution of stars, which ensures a deficit or even absence of stars with low angular momenta. Examples of systems with a loss cone are the centers of galaxies or star clusters with massive black holes. The instability can produce a flux of stars onto the galactic center, i.e., it can serve as a mechanism of fueling the nuclear activity of galaxies. Mathematically, the problem is reduced to analyzing simple characteristic equations that describe small perturbations in a disk and a sphere of radially highly elongated stellar orbits. In turn, these characteristics equations are derived through a number of successive simplifications of the general linearized Vlasov equations (i.e., the system that includes the collisionless Boltzmann kinetic equation and the Poisson equation) in action—angle variables.

Gravitational loss-cone instability

We study physical corollaries of the existing analogy between the simplest plasma traps (mirror traps) and star clusters surrounding massive black holes or dense galactic nuclei. There is a loss cone in the system through which plasma particles with low velocities transverse to the trap axis or, similarly, stars with low angular momenta (destroyed or absorbed by the central body) escape. The consequences of the “beam-like” deformation of the plasma distribution function in a trap are well known: a peculiar loss-cone instability producing a plasma flow into the loss cone develops as a result. We show that a similar gravitational loss-cone instability can also arise under certain conditions in the galactic case of interest to us. This instability is related to the slow precessional motions of highly elongated (nearly radial) stellar orbits and the main condition for its growth is that the precession of such orbits be retrograde (in the direction opposite to the orbital rotation of stars). Only under this condition do oscillations that can become unstable in the presence of a loss cone arise instead of the radial orbit instability (a variety of the Jeans instability in systems with highly elongated orbits) that takes place in the case of prograde precession. The instability produces a stream of stars onto the galactic center, i.e., serves as a mechanism of fueling the nuclear activity of galaxies. For a mathematical analysis, we have obtained relatively simple characteristic equations that describe small perturbations in a sphere of radially highly elongated stellar orbits. These characteristic equations are derived through a number of successive simplifications from a general linearized system of equations, including the collisionless Boltzmann kinetic equation and the Poisson equation (in action-angle variables). The central point of our analysis of the characteristic equations is preliminary detection of neutral modes (or proof of their absence in the case of stability).

Radial-orbit instability of a family of anisotropic Hernquist models with and without a supermassive black hole

We present a method to investigate the non-radial stability of a spherical anisotropic system that hosts a central supermassive black hole (SBH). Such systems have never been tested before for stability, although high anisotropies have been considered in the dynamical models that were used to estimate the masses of the central putative SBHs. A family of analytical anisotropic spherical Hernquist models with and without a black hole were investigated by means of N-body simulations. A clear trend emerges that the SBH has a significant effect on the overall stability of the system, that is, an SBH with a mass of a few per cent of the total mass of the galaxy can prevent or reduce the bar instabilities in anisotropic systems. Its mass determines not only the strength of the instability reduction, but also the time in which this occurs. These effects are most significant for models with strong radial anisotropies. Furthermore, our analysis shows that unstable systems with similar SBHs but with different anisotropy radii evolve differently: highly radial systems become oblate, while more isotropic models tend to form into prolate structures. In addition to this study, we also present a Monte Carlo algorithm to generate particles in spherical anisotropic systems.

On Universal Halos and the Radial Orbit Instability

The radial orbit instability drives dark matter halos toward a universal structure. This conclusion, first noted by Huss, Jain, and Steinmetz, is explored in detail through a series of numerical experiments involving the collapse of an isolated halo into the non-linear regime. The role played by the radial orbit instability in generating the density profile, shape, and orbit structure is carefully analyzed and, in all cases, the instability leads to universality independent of initial conditions. New insights into the underlying physics of the radial orbit instability are presented.