Linear response of galactic haloes to adiabatic gravitational perturbations

We determine the response of a self-similar isothermal stellar system to small adiabatic gravitational perturbations. For odd spherical harmonics the response is identical to the response of the analogous isothermal fluid system. For even spherical harmonics the response can be regarded as an infinite series of wavetrains in log r , implying alternating compression and rarefaction in equal logarithmic radius intervals. Partly because of the oscillatory nature of the solutions, tidal fields from external sources are not strongly amplified by an intervening isothermal stellar system, except at radii ≲10^−3.5 times the satellite radius; at some radii the stellar system can even screen the external tidal field in a manner analogous to Debye screening. As Weinberg has pointed out, individual resonances in a stellar system can strongly amplify external tidal fields over a limited radial range, but we cannot address this possibility because we examine only adiabatic perturbations. We also discuss the application of our method to the halo response caused by the slow growth of an embedded thin disc.     

The kinematics of the ionized gas in the Circinus galaxy

The TAURUS-2 Fabry–Perot interferometer, mounted on the 3.9-m Anglo-Australian Telescope, has been used to observe the Circinus galaxy. We have mapped the intensity and velocity distribution of the ionized hydrogen in the galaxy using the Balmer series Hα spectral line.
The semiresolved core (observed with a seeing disc of 30 pc) appears amorphous in shape, which is commonly observed in Seyfert 2 galaxies. Its peak coincides with the core position measured in the radio continuum, suggesting that ionized gas surrounds a non-thermal source.
A circumnuclear ring or spiral of radius 220 pc and a rotational velocity of 350 km s⁻¹ (assuming circular motions) surrounds the core. The inclination angle of this feature, i =40°±10°, is less than that of the previously observed radio continuum disc. The velocity channel maps obtained for the Hα ring show that the kinematics resemble those of a rotating ring and the intensity displays a complex structure indicative of several, unresolved, H II regions. We believe the ring to represent a circumnuclear starburst.
Our Hα data also show the presence of the previously detected [O III ] ionization cone to the north-west of the core, measuring more than 400 pc in length. We suggest that the ionization cone lies in a different plane from that of the starburst ring and is directed away from us. Several kinematic components of the core are derived and we calculate an outflow velocity in excess of 150–200 km s⁻¹ for gas above the core of Circinus. We also present evidence for inflowing ionized gas at the centre of Circinus.
The correlation of the Hα and radio continuum features is discussed, as well as the possible presence of a starburst-driven superwind in the Circinus galaxy.     

Orbit classification in arbitrary 2D and 3D potentials

A method of classifying generic orbits in arbitrary 2D and 3D potentials is presented. It is based on the concept of spectral dynamics introduced by Binney &38; Spergel that uses the Fourier transform of the time series of each coordinate. The method is tested using a number of potentials previously studied in the literature and is shown to distinguish correctly between regular and irregular orbits, to identify the various families of regular orbits (boxes, loops, tubes, boxlets, etc.), and to recognize the second-rank resonances that bifurcate from them. The method returns the position of the potential centre and, for 2D potentials, the orientation of the principal axes as well, should this be unknown. A further advantage of the method is that it has been encoded in a FORTRAN program that does not require user intervention, except for 'fine tuning' of search parameters that define the numerical limits of the code. The automatic character makes the program suitable for classifying large numbers of orbits.     

Frequency map analysis of the orbital structure in elliptical galaxies

We present an application of the frequency map analysis to an elliptical galaxy which is represented by a generalization of a double-power-law spherical mass model. The density distribution of this model varies as r ^−γ close to the centre and as r ⁻⁴ at large radii. We study the case with γ = 1, which is known as the 'weak-cusp' model and which represents well the density profile of the 'core' galaxies observed by the Hubble Space Telescope . The final objective of our work is to improve our understanding of the dynamics of elliptical galaxies in a similar way to Merritt &38; Fridman, finding the regions of stochasticity, looking for resonances that might play an important role in sustaining the triaxial morphology, and analysing the diffusion of orbits. To this end, we use the frequency map analysis of Laskar, which has been applied widely in the field of celestial mechanics but which is a relatively new technique in the area of galactic dynamics. Finally, we show some useful features of this method in understanding the global dynamical structure of the system.     

Orbital motion of a massive object in a dense stellar system

In this paper we show the positional oscillation of a massive object in a dense stellar system by numerical N -body simulations. We found that the central massive object, which at first is placed at rest at the centre of the surrounding spherical stellar system, promptly departs from the centre and rotates in accordance with the rotation of the stellar system, if the stellar system has an appreciable rotation. This oscillatory motion continues for a long time because of the absence of dynamical friction. Such a long-lasting oscillation may explain the asymmetric structure observed in the centres of M31 and NGC 4486B, may cause the secular flow of gaseous elements distributed in the central regions of galaxies on to the massive object, and may ignite activity in the centres of galaxies.     

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.     

Dynamical friction force exerted on spherical bodies

We present a rigorous calculation of the dynamical friction force exerted on a spherical massive perturber moving through an infinite homogenous system of field stars. By calculating the shape and mass of the polarization cloud induced by the perturber in the background system, which decelerates the motion of the perturber, we recover Chandrasekhar's drag force law with a modified Coulomb logarithm. As concrete examples, we calculate the drag force exerted on a Plummer sphere or a sphere with the density distribution of a Hernquist profile. It is shown that the shape of the perturber only affects the exact form of the Coulomb logarithm. The latter converges on small scales because encounters of the test and field stars with impact parameters less than the size of the massive perturber become inefficient. We confirm in this way the earlier results based on the impulse approximation of small angle scatterings.