Continuing work initiated in earlier publications [Ichita, Yamada and Asada, Phys. Rev. D {\bf 83}, 084026 (2011); Yamada and Asada, Phys. Rev. D {\bf 86}, 124029 (2012)], we examine the post-Newtonian (PN) effects on the stability of the triangular solution in the relativistic three-body problem for general masses. For three finite masses, a condition for stability of the triangular solution is obtained at the first post-Newtonian (1PN) order, and it recovers previous results for the PN restricted three-body problem when one mass goes to zero. The stability regions still exist even at the 1PN order, though the PN triangular configuration for general masses is less stable than the PN restricted three-body case as well as the Newtonian one.

The physical state and properties of silicates at conditions encountered in the cores of gas giants, ice giants and of Earth like exoplanets now discovered with masses up to several times the mass of the Earth remains mostly unknown. Here, we report on theoretical predictions of the properties of silica, SiO$_2$, up to 4 TPa and about 20,000K using first principle molecular dynamics simulations based on density functional theory. For conditions found in the Super-Earths and in ice giants, we show that silica remains a poor electrical conductor up to 10 Mbar due to an increase in the Si-O coordination with pressure. For Jupiter and Saturn cores, we find that MgSiO$_3$ silicate has not only dissociated into MgO and SiO$_2$, as shown in previous studies, but that these two phases have likely differentiated to lead to a core made of liquid SiO$_2$ and solid (Mg,Fe)O.

We study the survival of substructures (clumps) within larger self-gravitating dark matter halos. Building on scaling relations obtained from N-body calculations of violent relaxation, we argue that the tidal field of galaxies and halos can only destroy substructures if spherical symmetry is imposed at formation. We explore other mechanisms that may tailor the number of halo substructures during the course of virialization. Unless the larger halo is built up from a few large clumps, we find that clump-clump encounters are unlikely to homogenize the halo on a dynamical timescale. Phase mixing would proceed faster in the inner parts and allow for the secular evolution of a stellar disk.

Non-linear effects in the dynamical evolution of a shearing sheet made of stars are studied. First the implications of hitherto neglected non-linearities of the Boltzmann equation for the dynamical evolution of the shearing sheet are investigated. Using a formalism developed previously on the basis of the linearized Boltzmann equation it is demonstrated that the inclusion of the non-linear term leads to a feedback cycle for swing amplified density waves in the unbounded shearing sheet. Such a feedback is unique to star disks and is not known for gas disks. In order to present concrete examples of the non-linear feedback cycle a SCF code was developed and numerical simulations of the dynamical evolution of the shearing sheet are performed. The numerical results seem to confirm the theoretical predictions. The evolution of the shearing sheet resembles closely and might actually explain the recurrent spiral instabilities found in large-scale numerical simulations of the dynamical evolution of galactic disks.

We perform N-body simulations of the dynamical evolution of a dwarf galaxy falling into the Milky Way galaxy in order to understand the formation scenario of the peculiar globular cluster $\omega$ Centauri. We use self-consistent models of the bulge and the disc of the Milky Way, as well as of the dwarf galaxy, and explore a range of dwarf models with different density distributions. Namely, we use King (1966) and Hernquist (1990) density profiles to model the density distribution in the dwarf. The central region of our King model has a density profile approximately $\propto r^{-2}$, while that of the Hernquist model is $\propto r^{-1}$. The difference in the dwarf's density distributions leads to distinct evolutionary scenarios. The King model dwarf loses its mass exponentially as a function of apocentric distance, with the mass-loss rate depending on the initial mass and size of the dwarf. Regardless of the initial mass and size, the King model dwarf remains more massive than $10^8$ \msun after a few Gyr of evolution. The Hernquist model dwarf experiences an accelerated mass loss, and the mass of the remnant falls below $10^8$ \msun within a few Gyr. By exploring an appropriate set of parameters, we find a Hernquist model that can attain the mass and orbital characteristics of $\omega$ Cen after a few Gyr.

We explore an accretion origin for Omega Cen by N-body modeling of the orbital decay and disruption of a Milky-Way dwarf satellite. This work is focused on studying a particular satellite model that aims to reproduce the present orbit of Omega Cen, as recently determined from absolute proper motions. The model satellite is launched from 58 kpc from the Galactic Center, on a radial, low-inclination orbit. We find that a capture scenario can produce an Omega Cen-like object with the current low-energy orbit of the cluster. Our best model is a nucleated galaxy with a Hernquist density profile that has a mass of 8 10**8 Msun, and a half-mass radius of 1.4 kpc.

Multi-scale interaction between the LMC, Galactic halo, and the disk is examined with N-body simulations, and precise amplitudes of the Galactic warp excitation are obtained. The Galactic models are constructed most realistically to satisfy available observational constraints on the local circular velocity, the mass, surface density and thickness of the disk, the mass and size of the bulge, the local density of the halo matter at the solar radius, and the mass and orbit of the LMC. By using hybrid algorithm (SCF-TREE) I have succeeded to follow the evolution with millions of particles. The orbiting satellite excites density enhancement as a wake, and the wake exerts a tidal force on the disk. Because of the additional torque from the wakes in the halo, the amplitudes of the induced warps are much larger than the classical estimate by [hunt1969], who considered only the direct torque from the LMC. The obtained amplitudes of m=0, 1, 2 warps in the larger halo model show very good agreement with the observed amplitude in the Milky Way. This result revives the LMC as a possible candidate of the origin of the Galactic warp. Our smaller halo model, however, yield only weak warps in all the harmonic modes. Therefore the halo still has significant influence on excitation of warp even in the interaction scenario for excitation of warps.

Motivated by the presence of numerous dark matter clumps in the Milky Way's halo as expected from the cold dark matter cosmological model, we conduct numerical simulations to examine the heating of the disk. We construct an initial galaxy model in equilibrium, with a stable thin disk. The disk interacts with dark matter clumps for about 5 Gyr. Three physical effects are examined : first the mass spectrum of the dark matter clumps, second the initial thickness of the galactic disk, and third the spatial distribution of the clumps. We find that the massive end of the mass spectrum determines the amount of disk heating. Thicker disks suffer less heating. There is a certain thickness at which the heating owing to the interaction with the clumps becomes saturates. We also find that the heating produced by the model which mimics the distribution found in Standard CDM cosmology is significant and too high to explain the observational constraints. On the other hand, our model that corresponds to the clump distribution in a $\Lambda$CDM cosmology produces no significant heating. This result suggests that the $\Lambda$CDM cosmology is preferable with respect to the Standard CDM cosmology in explaining the thickness of the Milky Way.

Using N-body numerical simulations, we examine possible mechanisms for the origin of the intracluster starlight features recently discovered in the Coma cluster of galaxies. We show that tidal interactions of a ``normal'' elliptical galaxy with a symmetrically distributed dark matter potential do not produce the observed intracluster starlight features. A head-on collision of two normal ellipticals does, however, explain the origin of the intracluster features with surface brightness about 26 mag/arcsec$^2$). Another possible explanation for the intracluster starlight features is galactic tidal interactions with massive intracluster substructures. The presence of substructures in Coma with a density of $\sim$ 0.1 times the density in the center of the galaxy, and with a total mass about three orders of magnitude larger than the mass of a normal elliptical galaxy, would account for the observations.

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.

The relation between relaxation, the time scale of Lyapunov instabilities, and the Kolmogorov-Sinai time in a one-dimensional gravitating sheet system is studied. Both the maximum Lyapunov exponent and the Kolmogorov-Sinai entropy decrease as proportional to $N^{-1/5}$. The time scales determined by these quantities evidently differ from any type of relaxation time found in the previous investigations. The relaxation time to quasiequilibria (microscopic relaxation) is found to coincide with the inverse of the minimum positive Lyapunov exponent. The relaxation time to the final thermal equilibrium differs to the inverse of the Lyapunov exponents and the Kolmogorov-Sinai time.

Cold collapse of a cluster composed of small identical clumps, each of which is in virial equilibrium, is considered. Since the clumps have no relative motion with respect to each other initially, the cluster collapses by its gravity. At the first collapse of the cluster, most of the clumps are destroyed, but some survive. In order to find the condition for the clumps to survive, we made systematic study in two-parameter space: the number of the clumps $N_c$ and the size of the clump $r_v$. We obtained the condition, $N_c \gg 1$ and $n_k \geq 1$, where $n_k$ is related to $r_v$ and the initial radius of the cluster $R_{ini}$ through the relation $R_{ini}/r_v = 2 N_c^{(n_k+5)/6}$. A simple analytic argument supports the numerical result. This $n_k$ corresponds to the index of the power spectrum of the density fluctuation in the cosmological hierarchical clustering, and thus our result may suggest that in the systems smaller than $2/(\Omega h^2)$Mpc, the first violent collapse is strong enough to sweep away all substructures which exist before the collapse.

Using N-body and SPH modeling we perform 3D numerical simulations of head-on collisions between gas rich disk galaxies, including collisions between counter-rotating disks and off-center collisions. Pure stellar intruders do not produce gaseous plumes similar to those seen in the Cartwheel and VII Zw466 complexes of interacting galaxies; the presence of gas in an intruder galaxy and radiative cooling are important for the formation of a gaseous plume extending from the disk of a target galaxy. A noticeable plume structure can be formed if the mass of an intruder is a few percent of the mass of the primary. The halo of the intruder is stripped in the collision, and dispersed particles form a broad stellar bridge connecting the two galaxies. The fraction of the intruder's halo dispersed in the collision depends on the total mass of the intruder, and low-mass intruders lose most of their mass.

This is the third paper of the series of our studies of the one-dimensional self-gravitating many-body systems. In this paper, we thus study the transition phenomena after the first transition from a quasiequilibrium. We found that irrespective of the initial conditions, the system wanders between many states which yield the same properties of the quasiequilibrium as the water-bag. This itinerancy is most prominent in the time scale $t\sim 10^5 \sim 10^6 t_c$. In the midway between two succeeding quasiequilibria, the system experiences a transient state, where one particle keeps exclusive high energy and its motion decouples with the others. Though the transient state is not the thermally relaxed equilibrium, its distribution quite resembles the isothermal distribution. Thus the macroscopic relaxation we discussed in the previous papers corresponds to the transition from one of quasiequilibria to a transient state. Averaging the behavior over a time scale much longer than the macroscopic relaxation time gives the isothermal distribution. Distribution of lifetime $\tau$ of the transient states yields a power-law distribution of $\tau^{-2}$. The transient state gives a clear example of chaotic itinerancy in conserved dynamical systems. The mechanism of the onset of itinerancy is examined.

Though one dimensional self-gravitating $N$-body systems have been studied for three decade, the nature of relaxation was still unclear. There were inconsistent results about relaxation time; some initial state relaxed in the time scale $T\sim N\,t_c$, but another state did not relax even after $T\sim N^2\,t_c$, where $t_c$ is the crossing time. The water-bag distribution was believed not to relax after $T\sim N^2\,t_c$. In our previous paper, however, we found there are two different relaxation times in the water-bag distribution;in the faster relaxation ( microscopic relaxation ) the equipartition of energy distribution is attains but the macroscopic distribution turns into the isothermal distribution in the later relaxation (macroscopic relaxation). In this paper, we investigated the properties of the two relaxation. We found that the microscopic relaxation time is $T\sim N\,t_c$, and the macroscopic relaxation time is proportional to $N\,t_c$, thus the water-bag does relax. We can see the inconsistency about the relaxation times is resolved as that we see the two different aspect of relaxations. Further, the physical mechanisms of the relaxations are presented.

The microscopic dynamics of one-dimensional self-gravitating many-body systems is studied. We examine two courses of the evolution which has the isothermal and stationary water-bag distribution as initial conditions. We investigate the evolution of the systems toward thermal equilibrium. It is found that when the number of degrees of freedom of the system is increased, the water-bag distribution becomes a quasi-equilibrium, and also the stochasticity of the system reduces. This results suggest that the phase space of the system is effectively not ergodic and the system with large degreees of freedom approaches to the near-integrable one.