Stability of power-law discs — I. The Fredholm integral equation

The power-law discs are a family of infinitesimally thin, axisymmetric stellar discs of infinite extent. The rotation curve can be rising, falling or flat. The self-consistent power-law discs are scale-free, so that all physical quantities vary as a power of radius. They possess simple equilibrium distribution functions depending on the two classical integrals, energy and angular momentum. While maintaining the scale-free equilibrium force law, the power-law discs can be transformed into cut-out discs by preventing stars close to the origin (and sometimes also at large radii) from participating in any disturbance. This paper derives the homogeneous Fredholm integral equation for the in-plane normal modes in the self-consistent and the cut-out power-law discs. This is done by linearizing the collisionless Boltzmann equation to find the response density corresponding to any imposed density and potential. The normal modes — that is, the self-consistent modes of oscillation — are found by requiring the imposed density to equal the response density. In practice, this scheme is implemented in Fourier space, by decomposing both imposed and response densities in logarithmic spirals. The Fredholm integral equation then relates the transform of the imposed density to the transform of the response density. Numerical strategies to solve the integral equation and to isolate the growth rates and the pattern speeds of the normal modes are discussed.     

Finding how many isolating integrals of motion an orbit obeys

The correlation dimension, that is the dimension obtained by computing the correlation function of pairs of points of a trajectory in phase space, is a numerical technique introduced in the field of non-linear dynamics in order to compute the dimension of the manifold in which an orbit moves, without the need of knowing the actual equations of motion that give rise to the trajectory. This technique has been proposed in the past as a method to measure the dimension of stellar orbits in astronomical potentials, that is the number of isolating integrals of motion the orbits obey. Although the algorithm can in principle yield that number, some care has to be taken in order to obtain good results. We studied the relevant parameters of the technique, found their optimal values, and tested the validity of the method on a number of potentials previously studied in the literature, using the Smaller Alignment Index (SALI), Lyapunov exponents and spectral dynamics as gauges.     

Analytical solutions for the Newtonian gravitational field induced by matter within axisymmetric boundaries

An analytical method originally applied to the problem of the actuator disc in fluid mechanics has been applied to the closely analogous problem of constructing the classical Newtonian potential and attractions. The method can treat axisymmetric problems and also non-axisymmetric cases where matter is confined within axisymmetric boundaries. The potential and attractions for the generalized thin finite disc can be given in closed form in terms of elliptic integrals and elementary functions. For the general case of matter within an axisymmetric boundary, the potentials and attractions can be evaluated as one-dimensional integrals of albeit complex analytical expressions. These expressions represent the fields induced by matter in an extended region as a distribution of gravitating discs. For certain special cases, such as matter bounded by a circular cylinder and also for matter distributed in a spherical region, closed-form solutions can be given that appear to be new. Some non-axisymmetric results are also given for the thin disc of infinite radial extent.     

Self-consistent response of a galactic disc to vertical perturbations

We study the self-consistent, linear response of a galactic disc to vertical perturbations, as induced, say, by a tidal interaction. We calculate the self-gravitational potential corresponding to a non-axisymmetric, self-consistent density response of the disc using the Green's function method. The response potential is shown to oppose the perturbation potential because the self-gravity of the disc resists the imposed potential, and this resistance is stronger in the inner parts of a galactic disc. For the   m = 1  azimuthal wavenumber, the disc response opposes the imposed perturbation up to a radius that spans a range of 4–6 disc scalelengths, so that the disc shows a net warp only beyond this region. This physically explains the well known but so far unexplained observation that warps typically set in beyond this range of radii. We show that the inclusion of a dark matter halo in the calculation only marginally changes (by ∼10 per cent) the radius for the onset of warps. For perturbations with higher azimuthal wavenumbers, the net signature of the vertical perturbations can only be seen at larger radii – for example, beyond 7 exponential disc scalelengths for   m = 10  . Also, for the high- m cases, the magnitude of the negative disc response due to the disc self-gravity is much smaller. This is shown to result in corrugations of the mid-plane density, which explains the puzzling scalloping with   m = 10  detected in H  i in the outermost regions ∼30 kpc in the Galaxy.     

Chaos in cuspy triaxial galaxies with a supermassive black hole: a simple toy model

This paper summarises an investigation of chaos in a toy potential which mimics much of the behaviour observed for the more realistic triaxial generalisations of the Dehnen potentials, which have been used to model cuspy triaxial galaxies both with and without a supermassive black hole. The potential is the sum of an anisotropic harmonic oscillator potential, ${\text{V}}_{\text{0}} = \frac{1}{2}\left( {a^2 x^2 + b^2 y^2 + c^2 z^2 } \right)$ , and aspherical Plummer potential, ${\text{V}}_{\text{P}} = M_{BH} /\sqrt {r^2 + \varepsilon ^2 } $ , with $r^2 = x^2 + y^2 + z^2$ . Attention focuses on three issues related tothe properties of ensembles of chaotic orbits which impact on chaotic mixing and the possibility of constructing self-consistent equilibria:(1) What fraction of the orbits are chaotic? (2) How sensitive are the chaotic orbits, that is, how large are their largest (short time) Lyapunov exponents? (3) To what extent is the motion of chaotic orbits impeded by Arnold webs, that is, how 'sticky' are the chaotic orbits? These questions are explored as functions of the axis ratio a: b: c, black hole mass M _BH, softening length ε, and energy E with the aims of understanding how the manifestations of chaos depend onthe shape of the system and why the black hole generates chaos. The simplicity of the model makes it amenable to a perturbative analysis. That it mimics the behaviour of more complicated potentials suggests that much of this behaviour should be generic.     

Gas-driven evolution of stellar orbits in barred galaxies

We carry out a detailed orbit analysis of gravitational potentials selected at different times from an evolving self-consistent model galaxy consisting of a two-component disc (stars+gas) and a live halo. The results are compared with a pure stellar model, subject to nearly identical initial conditions, which are chosen so as to make the models develop a large-scale stellar bar. The bars are also subject to hose-pipe (buckling) instability which modifies the vertical structure of the disc. The diverging morphological evolution of both models is explained in terms of gas radial inflow, the resulting change in the gravitational potential at smaller radii, and the subsequent modification of the main families of orbits, both in and out of the disc plane.   We find that dynamical instabilities become milder in the presence of the gas component, and that the stability of planar and 3D stellar orbits is strongly affected by the related changes in the potential — both are destabilized, with the gas accumulation at the centre. This is reflected in the overall lower amplitude of the bar mode and in the substantial weakening of the bar, which appears to be a gradual process. The vertical buckling of the bar is much less pronounced and the characteristic peanut shape of the galactic bulge almost disappears when there is a substantial gas inflow towards the centre. Milder instability results in a smaller bulge, the basic parameters of which are in agreement with observations. We also find that the overall evolution in the model with a gas component is accelerated because of the larger central mass concentration and the resulting decrease in the characteristic dynamical time.     

Simulations of spheroidal systems with substructure: trees in fields

We present a hybrid technique of N -body simulation to deal with collisionless stellar systems having an inhomogeneous global structure. We combine a treecode and a self-consistent field code such that each of the codes models a different component of the system being investigated. The treecode is suited to treatment of dynamically cold or clumpy components, which may undergo significant evolution within a dynamically hot system. The hot system is appropriately evolved by the self-consistent field code. This combined code is particularly suited to a number of problems in galactic dynamics. Applications of the code to these problems are briefly discussed.     

Orbits supporting bars within bars

High-resolution observations of the inner regions of barred disc galaxies have revealed many asymmetrical, small-scale central features, some of which are best described as secondary bars. Because orbital time-scales in the galaxy centre are short, secondary bars are likely to be dynamically decoupled from the main kiloparsec-scale bars. Here we show that regular orbits exist in such doubly barred potentials, and that they can support the bars in their motion. We find orbits in which particles remain on loops : closed curves which return to their original positions after two bars have come back to the same relative orientation. Stars trapped around stable loops could form the building blocks for a long-lived, doubly barred galaxy. Using the loop representation, we can find which orbits support the bars in their motion, and the constraints on the sizes and shapes of self-consistent double bars. In particular, it appears that a long-lived secondary bar may exist only when an inner Lindblad resonance is present in the primary bar, and that it would not extend beyond this resonance.     

Hydrostatic gas distributions: global estimates of temperature and abundance

Estimating the temperature and metal abundance of the intracluster and the intragroup media is crucial to determine their global metal content and to determine fundamental cosmological parameters. When a spatially resolved temperature or abundance profile cannot be recovered from observations (e.g. for distant objects), or deprojection is difficult (e.g. due to a significant non-spherical shape), only global average temperature and abundance are derived. After introducing a general technique to build hydrostatic gaseous distributions of prescribed density profile in potential wells of any shape, we compute the global mass-weighted and emission-weighted temperature and abundance for a large set of barotropic equilibria and an observationally motivated abundance gradient. We also compute the spectroscopic-like temperature that is recovered from a single temperature fit of observed spectra. The derived emission-weighted abundance and temperatures are higher by 50 to 100 per cent than the corresponding mass-weighted quantities, with overestimates that increase with the gas mean temperature. Spectroscopic temperatures are intermediate between mass and luminosity-weighted temperatures. Dark matter flattening does not lead to significant differences in the values of the average temperatures or abundances with respect to the corresponding spherical case (except for extreme cases).     

Dynamical evolution of galactic disks driven by interaction with a satellite

Dynamical evolution of galactic disks driven by interaction with satellite galaxies, particularly the problem of the disk warping and thickening is studied numerically. One of the main purpose of the study is to resolve the long standing problem of the origin of the disk warping. A possible cause of the warp is interaction with a satellite galaxy. In the case of the Milky Way, the LMC has been considered as the candidate. Some linear analysis have already given a positive result, but one had to wait for a fully self-consistent simulation as a proof. I have accomplished the numerical simulations with a million particles, by introducing a hybrid algorithm, SCF-TREE. Those simulations give us quantitative estimates for the Milky Way system. We have found an example in which large warp amplitudes are developed. We also found that the warp amplitudes depend on the halo distribution. Among our three models, the most massive and spherical halo is preferable for the observable warp excitation. This revised version was published online in August 2006 with corrections to the Cover Date.     

Stability of power-law discs — II. The global spiral modes

This paper reports on the in-plane normal modes in the self-consistent and the cut-out power-law discs. Although the cut-out discs are remarkably stable to bisymmetric perturbations, they are very susceptible to one-armed modes. For this harmonic, there is no inner Lindblad resonance, thus removing a powerful stabilizing influence. A physical mechanism for the generation of the one-armed instabilities is put forward. Incoming trailing waves are reflected as leading waves at the inner cut-out, thus completing the feedback for the swing-amplifier. Growing three-armed and four-armed modes occur only at very low temperatures. However, neutral m  = 3 and m  = 4 modes are possible at higher temperatures for some discs. The rotation curve index β has a marked effect on stability. For all azimuthal wavenumbers, any unstable modes persist to higher temperatures and grow more vigorously if the rotation curve is rising (β < 0) than if the rotation curve is falling (β > 0). If the central regions or outer parts of the disc are carved out more abruptly, any instabilities become more virulent. The self-consistent power-law discs possess a number of unusual stability properties. There is no natural time-scale in the self-consistent disc. If a mode is admitted at some pattern speed and growth rate, then it must be present at all pattern speeds and growth rates. Our analysis — although falling short of a complete proof — suggests that such a two-dimensional continuum of non-axisymmetric modes does not occur and that the self-consistent power-law discs admit no global non-axisymmetric normal modes whatsoever. Without reflecting boundaries or cut-outs, there is no resonant cavity and no possibility of unstable growing modes. The self-consistent power-law discs certainly admit equi-angular spirals as neutral modes, together with a one-dimensional continuum of growing axisymmetric modes.     

Three-integral models for axisymmetric galactic discs

We present new equilibrium component distribution functions that depend on three analytic integrals in a Stäckel potential, and that can be used to model stellar discs of galaxies. These components are generalizations of two-integral ones and can thus provide thin discs in the two-integral approximation. Their most important properties are the partly analytical expression for their moments, the disc-like features of their configuration space densities (exponential decline in the galactic plane and finite extent in the vertical direction) and the anisotropy of their velocity dispersions. We further show that a linear combination of such components can fit a van der Kruit disc.     

Stellar dynamics in razor-thin discs with massive nuclear black holes

The bifurcations of orbit-averaged dynamics are studied in a class of razor-thin discs with central black holes. The model used here consists of a perturbed harmonic oscillator Hamiltonian augmented with a GM r potential. Through a sequence of conformal and canonical transformations, we reduce the phase-space flows of the system to a set of non-linear differential equations on a sphere. Based on the critical points of the averaged system, we classify orbit families and reveal the existence of six types of periodic motions: circular , long - and short-axis elliptical , long - and short-axis radial and inclined radial orbits. Long-axis elliptical orbits and their surrounding tubes have significant features: whilst they keep stars away from the centre, they elongate in the same direction as the density profile. These properties are helpful in the construction of self-consistent equilibria.     

Global m= 1 instabilities and lopsidedness in disc galaxies

Lopsidedness is common in spiral galaxies. Often, there is no obvious external cause, such as an interaction with a nearby galaxy, for such features. Alternatively, the lopsidedness may have an internal cause, such as a dynamical instability. In order to explore this idea, we have developed a computer code that searches for self-consistent perturbations in razor-thin disc galaxies and performed a thorough mode-analysis of a suite of dynamical models for disc galaxies embedded in an inert dark matter halo with varying amounts of rotation and radial anisotropy.
Models with two equal-mass counter-rotating discs and fully rotating models both show growing lopsided modes. For the counter-rotating models, this is the well-known counter-rotating instability, becoming weaker as the net rotation increases. The m = 1 mode of the maximally rotating models, on the other hand, becomes stronger with increasing net rotation. This rotating m = 1 mode is reminiscent of the eccentricity instability in near-Keplerian discs.
To unravel the physical origin of these two different m = 1 instabilities, we studied the individual stellar orbits in the perturbed potential and found that the presence of the perturbation gives rise to a very rich orbital behaviour. In the linear regime, both instabilities are supported by aligned loop orbits. In the non-linear regime, other orbit families exist that can help support the modes. In terms of density waves, the counter-rotating m = 1 mode is due to a purely growing Jeans-type instability. The rotating m = 1 mode, on the other hand, grows as a result of the swing amplifier working inside the resonance cavity that extends from the disc centre out to the radius where non-rotating waves are stabilized by the model's outwardly rising Q profile.     

Approximate third integrals for axisymmetric potentials using local Stäckel fits

We use a set of Stäckel potentials to obtain a local approximation for an effective third integral in axisymmetric systems. We present a study on the feasibility and effectiveness of this approach. We apply it to three trial potentials of various flattenings, corresponding to nearly ellipsoidal, discy and boxy density isophotes. In all three cases, a good fit to the potential requires only a small set of Stäckel potentials, and the associated Stäckel third integral provides a very satisfactory, yet analytically simple, approximation to the trial potentials effective third integral.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

Constraining black hole masses from stellar kinematics by summing over all possible distribution functions

When faced with the task of constraining a galaxy's potential given limited stellar kinematical information, what is the best way of treating the galaxy's unknown distribution function (DF)? Using the example of estimating black hole (BH) masses, I argue that the correct approach is to consider all possible DFs for each trial potential, marginalizing the DF using an infinitely divisible prior. Alternative approaches, such as the widely used maximum-penalized likelihood method, neglect the huge degeneracies inherent in the problem and simply identify a single, special DF for each trial potential.
Using simulated observations of toy galaxies with realistic amounts of noise, I find that this marginalization procedure yields significantly tighter constraints on BH masses than the conventional maximum-likelihood method, although it does pose a computational challenge which might be solved with the development of a suitable algorithm for massively parallel machines. I show that in practice the conventional maximum-likelihood method yields reliable BH masses with well-defined minima in their χ² distributions, contrary to claims made by Valluri, Merritt & Emsellem.     

Application of the global modal approach to the spiral galaxies

We have tested the applicability of the global modal approach in the density wave theory of spiral structure for a sample of spiral galaxies with measured axisymmetric background properties. We report here the results of the simulations for four galaxies: NGC 488, NGC 628, NGC 1566, and NGC 3938. Using the observed radial distributions for the stellar velocity dispersions and the rotation velocities we have constructed the equilibrium models for the galactic disks in each galaxy and implemented two kinds of stability analyses - the linear global analysis and 2D-nonlinear simulations. In general, the global modal approach is able to reproduce the observed properties of the spiral arms in the galactic disks. The growth of spirals in the galactic disks can be physically understood in terms of amplification by over-reflection at the corotation resonance. Our results support the global modal approach as a theoretical explanation of spiral structure in galaxies. This revised version was published online in August 2006 with corrections to the Cover Date.