Dynamical stability and evolution of the discs of Sc galaxies

We examine the local stability of galactic discs against axisymmetric density perturbations with special attention to the different dynamics of the stellar and gaseous components. In particular, the discs of the Milky Way and of NGC 6946 are studied. The Milky Way is shown to be stable, whereas the inner parts of NGC 6946, a typical Sc galaxy from the Kennicutt sample, are dynamically unstable. The ensuing dynamical evolution of the composite disc is studied by numerical simulations. The evolution is so fierce that the stellar disc heats up dynamically on a short time-scale to such a high degree, which seems to contradict the morphological appearance of the galaxy. The star formation rate required to cool the disc dynamically is estimated. Even if the star formation rate in NGC 6946 is at present high enough to meet this requirement, it is argued that the discs of Sc galaxies cannot sustain such a high star formation rate for extended periods.     

Dynamical behaviour of the primitive asteroid belt

In this paper we consider the dynamical evolution and orbital stability of objects in the asteroid belt. A simple physical model, including full gravitational perturbations from both giant planets, is used to compute the dynamical evolution of 1000 test particles simulating the primitive asteroids. The criterion of planet crossing (or close approach in the case of resonant objects) is used to reject particles from the simulation. 44 per cent of the particles survived for the whole time-span covered by the numerical integration (∼10⁹ yr).
The 4:1, 3:1 and to a lesser extent the 2:1 Kirkwood gaps are formed in ∼10⁷ yr of evolution, representing direct numerical evidence about their gravitational origin.
We found that the rms eccentricity and inclination of the sample experience a fast increase during the first 10⁶ yr. The final rms eccentricity is 0.11, ∼60 per cent smaller than the present rms eccentricity (0.17). Nevertheless, the gravitational action of the giant planets suffices to prevent the formation of large objects, allowing catastrophic collisions and the subsequent depletion of material from this zone of the Solar system. The excited eccentricity by Jupiter and Saturn may favour mutual encounters and the further increase of the relative velocities up to their present values.     

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.     

Dynamical studies of cluster pairs in the Magellanic Clouds

We performed N -body simulations of star cluster encounters with Hernquist's TREECODE in a CRAY YMP-2E computer under different initial conditions (relative positions and velocities, cluster sizes, masses and concentration degrees). The total number of particles per simulation ranged from 1024 to 20480. These models are compared with a series of isodensity maps of cluster pairs in the Magellanic Clouds. Evidence is found that during the interactions, transient morphological effects such as an expanded halo, isophotal deformation and isophotal twisting can occur as a result of tidal effects and dynamical friction. The simulations also show that different outcomes are possible depending on the initial parameters: (i) long-standing changes of concentration degree can occur after the collision; (ii) one member can disaggregate; or (iii) the pair can coalesce into a single cluster with a distinct structure compared with the original ones. These simulations can reproduce a wide range of morphological structures in observed cluster pairs.     

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.     

On a time-symmetric Hermite integrator for planetary N-body simulation

We describe a P(EC) ⁿ Hermite scheme for planetary N -body simulation. The fourth-order implicit Hermite scheme is a time-symmetric integrator that has no secular energy error for the integration of periodic orbits with time-symmetric time-steps. In general N -body problems, however, this advantage is of little practical significance, since it is difficult to achieve time-symmetry with individual variable time-steps. However, we can easily enjoy the benefit of the time-symmetric Hermite integrator in planetary N -body systems, where all bodies spend most of the time on nearly circular orbits. These orbits are integrated with almost constant time-steps even if we adopt the individual time-step scheme. The P(EC) ⁿ Hermite scheme and almost constant time-steps reduce the integration error greatly. For example, the energy error of the P(EC)² Hermite scheme is two orders of magnitude smaller than that of the standard PEC Hermite scheme in the case of an N  = 100,  m  = 10²⁵ g planetesimal system with the rms eccentricity 〈 e ²〉^1/2 ≲0.03.     

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.     

Stellar dynamics in a galactic centre surrounded by a massive accretion disc — I. Newtonian description

The long-term evolution of stellar orbits bound to a massive centre is studied in order to understand the cores of star clusters in central regions of galaxies. Stellar trajectories undergo tiny perturbations, the origins of which are twofold: (i) the gravitational field of a thin gaseous disc surrounding the galactic centre, and (ii) cumulative drag arising from successive interactions of the stars with the material of the disc. Both effects are closely related because they depend on the total mass of the disc, assumed to be a small fraction of the central mass. It is shown that, in contrast to previous works, most of the retrograde (with respect to the disc) orbits are captured by the central object, presumably a massive black hole. Initially prograde orbits are also affected, so that statistical properties of the central star cluster in quasi-equilibrium may differ significantly from those deduced in previous analyses.