We present a mean-field approach to simulating merging processes of two spherical collisionless stellar systems. This approach is realized with a self-consistent field (SCF) method in which the full spatial dependence of the density and potential of a system is expanded in a set of basis functions for solving Poisson's equation. In order to apply this SCF method to a merging situation where two systems are moving in space, we assign the expansion center to the center of mass of each system, the position of which is followed by a mass-less particle placed at that position initially. Merging simulations over a wide range of impact parameters are performed using both an SCF code developed here and a tree code. The results of each simulation produced by the two codes show excellent agreement in the evolving morphology of the merging systems and in the density and velocity dispersion profiles of the merged systems. However, comparing the results generated by the tree code to those obtained with the softening-free SCF code, we have found that in large impact parameter cases, a softening length of the Plummer type introduced in the tree code has an effect of advancing the orbital phase of the two systems in the merging process at late times. We demonstrate that the faster orbital phase originates from the larger convergence length to the pure Newtonian force. Other application problems suitable to the current SCF code are also discussed.

We present hydrodynamic simulations of gas clouds inflowing from the disk to a few hundred parsec region of the Milky Way. A gravitational potential is generated to include realistic Galactic structures by using thousands of multipole expansions that describe 6.4 million stellar particles of a self-consistent Galaxy simulation. We find that a hybrid multipole expansion model, with two different basis sets and a thick disk correction, accurately reproduces the overall structures of the Milky Way. Through non-axisymmetric Galactic structures of an elongated bar and spiral arms, gas clouds in the disk inflow to the nuclear region and form a central molecular zone (CMZ)-like nuclear ring. We find that the size of the nuclear ring evolves into ~240 pc at T~1500 Myr, regardless of the initial size. For most simulation runs, the rate of gas inflow to the nuclear region is equilibrated to ~0.02 M_sun/yr. The nuclear ring is off-centered, relative to the Galactic center, by the lopsided central mass distribution of the Galaxy model, and thus an asymmetric mass distribution of the nuclear ring arises accordingly. The vertical asymmetry of the the Galaxy model also causes the nuclear ring to be tilted along the Galactic plane. During the first ~100 Myr, the vertical frequency of the gas motion is twice that of the orbital frequency, thus the projected nuclear ring shows a twisted, infinity-like shape.

The formation and evolution of tidal tails like those observed around some globular clusters and dwarf satellite galaxies is examined with an N-body simulation. In particular, we analyse in detail the evolving tidal features of a one-component satellite that is moving on a highly eccentric orbit in the external field of a host galaxy potential like our own. The results show that every time the satellite approaches apogalacticon, a fresh pair of tidal tails becomes notably prominent, and that eventually, the satellite possesses multiple tidal tails via repeating apocentre passages. Accordingly, the number of observed tidal arms can be used as a tracer of the number of orbital periods that such a system has completed around the centre of its host galaxy. By identifying the arm particles included in each of the first three consecutively formed pairs of tidal tails, we find that each pair of tidal tails is practically identical to one another regarding the energy and angular-momentum distributions. In addition, we demonstrate that the density profiles of these three pairs of tidal tails at their first apogalacticons after formation agree well with one another. It therefore follows that the multiplicity in tidal features originates from the repeated episode of tidal-arm formation in the course of the precessing motion of a satellite.

The characteristics of the galactic bars that are prone to suffer a damaging impact from a massive central black hole are examined using flat stellar disks. We construct three disk model groups that consist of exponential disks with one type of velocity distribution and Kuzmin-Toomre disks with two different types of exact equilibrium distribution function. For each disk model group, three disks that have different typical Toomre's Q values are evolved to form bars through dynamical instability. Once a bar is fully developed, a black hole (BH), whose mass is 1% of the disk mass, is adiabatically added at the center of the disk. Our results show that lower-amplitude bars, that is, weaker bars are dissolved more easily by that BH. We have found that this destructibility is rooted in the characteristic feature that the bar formed spontaneously becomes shorter in length and rounder in shape with decreasing bar amplitude. Since such weaker bars are found to originate from colder disks in each disk model group, it follows that for a given form of velocity structure, the coldness of an initial disk determines whether the bar produced in that disk is favorable to dissolution induced by a massive central BH. In addition, the existence of bar-dissolved galaxies of the kind studied here is also discussed.

The time evolution of barred structures is examined under the influence of the external forces exerted by a spherical halo and by prolate halos. In particular, galaxy disks are placed in the plane including the major axis of prolate halos, whose configuration is often found in cosmological simulations. N-body disks in fixed external halo fields are simulated, so that bars are formed via dynamical instability. In the subsequent evolution, the bars in prolate halos dissolve gradually with time, while the bar pattern in a spherical halo remains almost unchanged to the end of the simulation. The decay times of the bars suggest that they can be destroyed in a time smaller than a Hubble time. Our results indicate that this dissolution process could occur in real barred galaxies, if they are surrounded by massive dark prolate halos, and the configuration adopted here is not unusual from the viewpoint of galaxy formation. For a prolate halo model, an additional simulation that is restricted to two-dimensional in-plane motions has also ended up with similar bar dissolution. This means that the vertical motions of disk stars do not play an essential role in the bar dissolution demonstrated here.

A recent observation with the Hipparcos satellite and some numerical simulations imply that the interaction between an oblate halo and a disc is inappropriate for the persistence of galactic warps. Then, we have compared the time evolution of galactic warps in a prolate halo with that in an oblate halo. The haloes were approximated as fixed potentials, while the discs were represented by N-body particles. We have found that the warping in the oblate halo continues to wind up, and finally disappears. On the other hand, for the prolate halo model, the precession rate of the outer disc increases when the precession of the outer disc recedes from that of the inner disc, and vice versa. Consequently, the warping in the prolate halo persisted to the end of the simulation by retaining the alignment of the line of nodes of the warped disc. Therefore, our results suggest that prolate haloes could sustain galactic warps. The physical mechanism of the persistence of warp is discussed on the basis of the torque between a halo and a disc and that between the inner and outer regions of the disc.

It is shown that the cuspy density distributions observed in the cores of elliptical galaxies can be realized by dissipationless gravitational collapse. The initial models consist of power-law density spheres such as $\rho\propto r^{-1}$ with anisotropic velocity dispersions. Collapse simulations are carried out by integrating the collisionless Boltzmann equation directly, on the assumption of spherical symmetry. From the results obtained, the extent of constant density cores, formed through violent relaxation, decreases as the velocity anisotropy increases radially, and practically disappears for extremely radially anisotropic models. As a result, the relaxed density distributions become more cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity anisotropy could be a key ingredient for the formation of density cusps in a dissipationless collapse picture. The velocity dispersions increase with radius in the cores according to the nearly power-law density distributions. The power-law index, $n$, of the density profiles, defined as $\rho\propto r^{-n}$, changes from $n\approx 2.1$ at intermediate radii, to a shallower power than $n\approx 2.1$ toward the centre. This density bend can be explained from our postulated local phase-space constraint that the phase-space density accessible to the relaxed state is determined at each radius by the maximum phase-space density of the initial state.

The influence of central black holes on the dynamical evolution of bars in disk galaxies is examined. In particular, we use numerical simulation to estimate the minimum mass black hole (BH) needed to destroy a bar. Initially, bars form in the disks via dynamical instability. Thereafter, once a bar is fully developed, a BH is adiabatically added at the center of the disk. To mitigate the global effects of gravitational softening, Poisson's equation for the disk is solved by expanding the density and potential of the galaxy in a set of basis functions. Our results indicate that a bar can be completely destroyed in a time much smaller than a Hubble time if the central mass exceeds about 0.5% of the disk mass. Since the implied minimum BH mass for bar destruction is of order 10^{8.5} solar masses for a typical disk galaxy, this process should not be a rare phenomenon. The bar amplitude decreases gradually with time after the BH is added, and the rate at which the bar is destroyed increases with increasing BH mass. This suggests that bar destruction arises from scattering of stars that support the bar as they pass close to the center.

A particle method for reproducing the phase space of collisionless stellar systems is described. The key idea originates in Liouville's theorem which states that the distribution function (DF) at time t can be derived from tracing necessary orbits back to t=0. To make this procedure feasible, a self-consistent field (SCF) method for solving Poisson's equation is adopted to compute the orbits of arbitrary stars. As an example, for the violent relaxation of a uniform-density sphere, the phase-space evolution which the current method generates is compared to that obtained with a phase-space method for integrating the collisionless Boltzmann equation, on the assumption of spherical symmetry. Then, excellent agreement is found between the two methods if an optimal basis set for the SCF technique is chosen. Since this reproduction method requires only the functional form of initial DFs but needs no assumptions about symmetry of the system, the success in reproducing the phase-space evolution implies that there would be no need of directly solving the collisionless Boltzmann equation in order to access phase space even for systems without any special symmetries. The effects of basis sets used in SCF simulations on the reproduced phase space are also discussed.

First, we have ensured that spherical nonrotating collisionless systems collapse with almost retaining spherical configurations during initial contraction phases even if they are allowed to collapse three-dimensionally. Next, on the assumption of spherical symmetry, we examine the evolution of velocity dispersions with collapse for the systems which have uniform or power-law density profiles with Maxwellian velocity distributions by integrating the collisionless Boltzmann equation directly. The results show that as far as the initial contraction phases are concerned, the radial velocity dispersion never grows faster than the tangential velocity dispersion except at small radii where the nearly isothermal nature remains, irrespective of the density profiles and virial ratios. This implies that velocity anisotropy as an initial condition should be a poor indicator for the radial orbit instability. The growing behavior of the velocity dispersions is briefly discussed from the viewpoint that phase space density is conserved in collisionless systems.