We propose the self-organized relaxation process which drives a collisionless self-gravitating system (SGS) to the equilibrium state satisfying local virial (LV) relation. During the violent relaxation process, particles can move widely within the time interval as short as a few free fall times, because of the effective potential oscillations. Since such particle movement causes further potential oscillations, it is expected that the system approaches the critical state where such particle activities, which we call gravitational fugacity, is independent of the local position as much as possible. Here we demonstrate that gravitational fugacity can be described as the functional of the LV ratio, which means that the LV ratio is a key ingredient estimating the particle activities against gravitational potential. We also demonstrate that LV relation is attained if the LV ratio exceeds the crticaial value $b=1$ everywhere in the bound region during the violent relaxation process. The local region which does not meet this criterion can be trapped into the pre-saturated state. However, small phase-space perturbation can bring the inactive part into the LV critical state.

We investigate the merging process in N-body self-gravitating system from the viewpoints of the local virial relation which is the relation between the local kinetic energy and the local potential. We compare both the density profile and the phase space density profile in cosmological simulations with the critical solutions of collisionless static state satisfying the local virial (LV) relation. We got the results that the critical solution can explain the characteristic density profile with the appropriate value of anisotropy parameter $\beta \sim 0.5$. It can also explain the power law of the phase space density profile in the outer part of a bound state. However, it fails in explaining the central low temperature part which is connected to the scale invariant phase space density. It can be well fitted to the critical solution with the higher value of$\beta \sim 0.75$. These results indicate that the LV relation is not compatible with the scale invariant phase space density in cosmological simulation.

We demonstrate that the quasi-equilibrium state in self-gravitating $N$-body system after cold collapse are uniquely characterized by the local virial relation using numerical simulations. Conversely assuming the constant local virial ratio and Jeans equation for spherically steady state system, we investigate the full solution space of the problem under the constant anisotropy parameter and obtain some relevant solutions. Especially, the local virial relation always provides a solution which has a power law density profile in both the asymptotic regions $r\to 0$ and $\infty$. This type of solutions observed commonly in many numerical simulations. Only the anisotropic velocity dispersion controls this asymptotic behavior of density profile.

We propose two hypotheses which characterize the quasi-equilibrium state that realizes after a cold collapse of self-gravitating system. The first hypothesis is the linear temperature-mass (TM) relation, which yields a characteristic non-Gaussian velocity distribution. The second is the local virial (LV) relation, which, combining the linear TM relation, determines a unique mass density profile as $\rho(r)={\rho}%_{0}r^{-4}e^{-r_0/r}$. Although this density profile is unphysical in the central region, the region is just inner a few percent around the center of a bound region in cumulative mass, which is beyond the resolution of our numerical simulations. Hence posing two hypotheses is compatible to the numerical simulations for almost the whole region of the virialized bound state. Actually, except for this inner part, this density profile fits well to the data of cold collapse simulations. Two families of spherical and isotropic models, polytropes and King models, are examined from a view point of these two hypotheses. We found that the LV relation imposes a strong constraint on these models: only polytropes with index $n \sim 5$ such as Plummer's model are compatible with the numerical results characterized by the two hypotheses among these models. King models with the concentration parameter $% c \sim 2$ violate the LV relation while they are consistent with the $R^{1/4} $ law for the surface brightness. Hence the above characteristics can serve as a guideline to build up the models for the bound state after a cold collapse, besides the conventional criteria concerning the asymptotic behavior.

The collisionless quasi-equilibrium state realized after the cold collapse of self-gravitating systems has two remarkable characters. One of them is the linear temperature-mass (TM) relation, which yields a characteristic non-Gaussian velocity distribution. Another is the local virial (LV) relation, the virial relation which holds even locally in collisionless systems through phase mixing such as cold-collapse. A family of polytropes are examined from a view point of these two characters. The LV relation imposes a strong constraint on these models: only polytropes with index $n \sim 5$ with a flat boundary condition at the center are compatible with the numerical results, except for the outer region. Using the analytic solutions based on the static and spherical Jeans equation, we show that this incompatibility in the outer region implies the important effect of anisotropy of velocity dispersion. Furthermore, the velocity anisotropy is essential in explaining various numerical results under the condition of the local virial relation.

We study the velocity distribution in spherical collapses and cluster-pair collisions by use of N-body simulations. Reflecting the violent gravitational processes, the velocity distribution of the resultant quasi-stationary state generally becomes non-Gaussian. Through the strong mixing of the violent process, there appears a universal non-Gaussian velocity distribution, which is a democratic (equal-weighted) superposition of many Gaussian distributions (DT distribution). This is deeply related with the local virial equilibrium and the linear mass-temperature relation which characterize the system. We show the robustness of this distribution function against various initial conditions which leads to the violent gravitational process. The DT distribution has a positive correlation with the energy fluctuation of the system. On the other hand, the coherent motion such as the radial motion in the spherical collapse and the rotation with the angular momentum suppress the appearance of the DT distribution.

Fractal structures and non-Gaussian velocity distributions are characteristic properties commonly observed in virialized self-gravitating systems such as galaxies or interstellar molecular clouds. We study the origin of these properties using the one-dimensional ring model which we newly propose in this paper. In this simple model, $N$ particles are moving, on a circular ring fixed in the three dimensional space, with mutual interaction of gravity. This model is suitable for accurate symplectic integration method by which we find the phase transition in this system from extended-phase to collapsed-phase through an interesting phase (\halo-phase) which has negative specific heat. In this intermediate energy scale, there appear scaling properties, non-thermal and non-Gaussian velocity distributions. In contrast, these peculiar properties are never observed in other \gas and \core phases. Particles in each phase have typical time scales of motion determined by the cutoff length $\xi$, the ring radius $R$ and the total energy $E$. Thus all relaxation patterns of the system are determined by these three time scales.

We study statistical mechanics of the self--gravitating system applying the cluster expansion method developed in solid state physics. By summing infinite series of diagrams, we derive a complex free energy whose imaginary part is related to the relaxation time of the system. Summation of another series yields two--point correlation function whose correlation length is essentially given by the Jeans wavelength of the system.

We apply the renormalization group (RG) method to examine the observable scaling properties in Newtonian cosmology. The original scaling properties of the equations of motion in our model are modified for averaged observables on constant time slices. In the RG flow diagram, we find three robust fixed points: Einstein-de Sitter, Milne and Quiescent fixed points. Their stability (or instability) property does not change under the effect of fluctuations. Inspired by the inflationary scenario in the early Universe, we set the Einstein-de Sitter fixed point with small fluctuations, as the boundary condition at the horizon scale. Solving the RG equations under this boundary condition toward the smaller scales, we find generic behavior of observables such that the density parameter $\Omega$ decreases, while the Hubble parameter $H$ increases for smaller averaging volume. The quantitative scaling properties are analyzed by calculating the characteristic exponents around each fixed point. Finally we argue the possible fractal structure of the Universe beyond the horizon scale.

We study the motion of test particle in static axisymmetric vacuum spacetimes and discuss two criteria for strong chaos to occur: (1) a local instability measured by the Weyl curvature, and (2) a tangle of a homoclinic orbit, which is closely related to an unstable periodic orbit in general relativity. We analyze several static axisymmetric spacetimes and find that the first criterion is a sufficient condition for chaos, at least qualitatively. Although some test particles which do not satisfy the first criterion show chaotic behavior in some spacetimes, these can be accounted for the second criterion.