We investigated the dynamical evolution of clusters of galaxies in virial equilibrium using Fokker-Planck models and self-consistent N-body models. In particular, we focused on the growth of a common halo, which is a cluster-wide halo formed by matter stripped from galaxies, and the development of a central density cusp. The Fokker-Planck models include the effects of two-body gravitational encounters both between galaxies and between galaxies and common halo particles. The effects of tidal mass stripping from the galaxies due to close galaxy-galaxy encounters and accompanying dissipation of the orbital kinetic energies of the galaxies were also taken into account in the Fokker-Planck models. We find that the results of the Fokker-Planck models are in excellent agreement with those of the N-body models regarding the growth of the common halo mass and the evolution of the cluster density profiles. In the central region of the cluster, a shallow density cusp, approximated by $\rho (r) \propto r^{-\alpha}$ ($\alpha \sim$ 1), develops. This shallow cusp results from the combined effects of two-body relaxation and tidal stripping. The cusp steepness, $\alpha$, weakly depends on the relative importance of the tidal stripping. When the effect of stripping is important, the central velocity dispersion decreases as the central density increases and, consequently, a shallow ($\alpha <2$) cusp is formed. In the limit of no stripping, usual gravothermal core collapse occurs, i.e. the central velocity dispersion increases as the central density increases with a steep ($\alpha >2$) cusp left. We conclude from our consideration of the origin of the cusp demonstrated here that shallow cusps should develop in real galaxy clusters.

We investigated the evolution of clusters of galaxies using self-consistent N-body simulations. We varied the initial model of galaxies and clusters, and studied the dependence of evolution on initial conditions. We found that the growth rate of the common halo depends only weakly on galaxy models. On the other hand, the growth rate depends strongly on cluster models. Initially the growth rate is higher for cluster models with higher central density. However, this high growth rate drops in a few crossing times of the cluster, and after several crossing times, roughly half of the total mass is in the common halo in all models we considered. In the central region of clusters density cusps with the profile \rho ~ r^{-1.2} develops regardless of the models of clusters and galaxies. We also found that the galaxies evolved so as to satisfy the relation between the masses of galaxies m_{gx} and their velocity dispersion \sigma_{gx} expressed as m_{gx} \propto \sigma_{gx}^{3~4} for all galaxy models as a consequence of their dynamical evolution through galaxy-galaxy interactions. We discuss the relation between our result and the observed Faber-Jackson relation.

We investigated the evolution of clusters of galaxies using self-consistent $N$-body simulations in which each galaxy was modeled by many particles. We carried out simulations for about 20 cases using different initial conditions. In all simulations, clusters were initially in virial equilibrium. We found that more than half of the total mass escaped from individual galaxies within a few crossing times of the cluster, and that a diffuse halo was formed. The growth rate of the common halo depended on the size of individual galaxies only weakly. The stripping of the mass from galaxies was mainly due to the interaction of galaxies, not due to the effect of the tidal field of the cluster potential. The amount of stripped mass was larger for galaxies in the central region than for those in the outer region, since the interactions were more frequent in the central region. As a result, a positive correlation between the distance from the center and the mass of the galaxy developed. The volume-density profile of the common halo is expressed as $\rho\propto r^{-1}$ in the central region. This mass distribution is consistent with the mass distribution in clusters estimated using X-ray observations.