The structure of geomagnetism is very complex, and there are still some problems left in magnetohydrodynamics (MHD) simulations. Recently, the macro spin model has been suggested. This is the idea that geomagnetism is described by interaction with many local dynamo elements (macro-spins). This model can reproduce some features of geomagnetism by solving equations of motion with only one variable. In this paper, we study this model to make more general by adding one variable. In this result, our model becomes possible for several things that are not treated in the previous study, for example, migration of the North (or South) Magnetic Pole, comparison with observed data of magnetic field distributions expressed in two directions, etc. Moreover, as a result of an application to the sun, we could reproduce periodic polarity reversals and the power index of the power spectrum, etc. In addition, we investigate the statistical properties of the pole migration.

We solve the equations of motion of a one-dimensional planar Heisenberg (or Vaks-Larkin) model consisting of a system of interacting macro-spins aligned along a ring. Each spin has unit length and is described by its angle with respect to the rotational axis. The orientation of the spins can vary in time due to random forcing and spin-spin interaction. We statistically describe the behaviour of the sum of all spins for different parameters. The term "domino model" in the title refers to the interaction among the spins. We compare the model results with geomagnetic field reversals and find strikingly similar behaviour. The aggregate of all spins keeps the same direction for a long time and, once in a while, begins flipping to change the orientation by almost 180 degrees (mimicking a geomagnetic reversal) or to move back to the original direction (mimicking an excursion). Most of the time the spins are aligned or anti-aligned and deviate only slightly with respect to the rotational axis (mimicking the secular variation of the geomagnetic pole with respect to the geographic pole). Reversals are fast compared to the times in between and they occur at random times, both in the model and in the case of the Earth's magnetic field.

Geomagnetism is characterized by intermittent polarity reversals and rapid fluctuations. We have recently proposed a coupled macro-spin model to describe these dynamics based on the idea that the whole dynamo mechanism is described by the coherent interactions of many small dynamo elements. In this paper, we further develop this idea and construct a minimal model for magnetic variations. This simple model naturally yields many of the observed features of geomagnetism: its time evolution, the power spectrum, the frequency distribution of stable polarity periods, etc. This model has coexistent two phases; i.e. the cluster phase which determines the global dipole magnetic moment and the expanded phase which gives random perpetual perturbations that yield intermittent polarity flip of the dipole moment. This model can also describe the synchronization of the spin oscillation. This corresponds to the case of sun and the model well describes the quasi-regular cycles of the solar magnetism. Furthermore, by analyzing the relevant terms of MHD equation based on our model, we have obtained a scaling relation for the magnetism for planets, satellites, sun, and stars. Comparing it with various observations, we can estimate the scale of the macro-spins.

Many scaling relations are observed for self-gravitating systems in the universe. We explore the consistent understanding of them from a simple principle based on the proposal that the collision-less dark matter fluid terns into a turbulent state, i.e. dark turbulence, after crossing the caustic surface in the non-linear stage. The dark turbulence will not eddy dominant reflecting the collision-less property. After deriving Kolmogorov scaling laws from Navier-Stokes equation by the method similar to the one for Smoluchowski coagulation equation, we apply this to several observations such as the scale-dependent velocity dispersion, mass-luminosity ratio, magnetic fields, and mass-angular momentum relation, power spectrum of density fluctuations. They all point the concordant value for the constant energy flow per mass: $0.3 cm^2/sec^3$, which may be understood as the speed of the hierarchical coalescence process in the cosmic structure formation.

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.

Self gravitating systems (SGS) in the Universe are generally thought to be non-extensive, and often show long-tails in various distribution functions. In principle, these non-Boltzmann properties are naturally expected from the peculiar property of gravity, long-range and unshielded. Therefore the ordinary Boltzmann statistical mechanics would not be applicable for these self gravitating systems (SGS) in its naive form. In order to step further, we quantitatively investigate the above two properties, non-extensivity and long-tails, by explicitly introducing various models of statistical mechanics. We use the data of CfA II South redshift survey and apply the count-in-cell method. We study four statistical mechanics, (1) Boltzmann, (2) Fractal, (3) R\'enyi, and (4) Tsallis, and use Akaike information criteria (AIC) for the fair comparison.

Large-scale astrophysical systems are non-extensive due to their long-range force of gravity. Here we show an approach toward the statistical mechanics of such self-gravitating systems (SGS). This is a generalization of the standard statistical mechanics based on the new definition of entropy; Tsallis statistical mechanics. Developing the composition of entropy and the generalized Euler relation, we investigate the galaxy distributions in count-in-cell method. This is applied to the data of CfA II South redshift survey.

We have searched for molecular gas in the cD galaxy of a poor cluster of galaxies AWM7 using Nobeyama 45 m telescope. We do not detect CO emission in the galaxy. Our limit of molecular gas in the inner 7.5 kpc is M_H_2< 4 x 10^8 M_sun. We estimate the total mass of molecular gas left in the cD galaxy when the gas deposited by a cooling flow once becomes molecular gas and the molecular gas is continuously evaporated by the ambient hot gas. The observational limit of molecular gas requires f >~ 10^{-3}, where f is the ratio of the heat conduction rate to that of Spitzer. However, this contradicts recent X-ray observations showing f<10^{-5}. Thus, the non-detection of CO cannot be explained by the evaporation, and most of the cooled gas predicted by a cooling flow model may not change into molecular gas in the cD galaxy. Moreover, we estimate the evaporation time of molecular clouds brought to a cD galaxy through the capture of gas-rich galaxies and find that these clouds should not be evaporated if f <~ 10^{-3}-10^{-4}. Therefore, the non-detection of CO in a cD galaxy could constrain the total mass of the molecular clouds brought into it.

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 explore the cosmological model in which a late-time mild inflation is realized after the star formation epoch. Non-vanishing curvature coupling of a classical boson field yields this mild inflation without a cosmological constant. Accordingly the lifetime of the present Universe is remarkably increased in our model. Thus we show that the present observed high value of the Hubble parameter $H_0 \approx 70-80{\rm km/sec/Mpc}$ is compatible with the age of the oldest stars $14{\rm Gyr}$ without introducing the cosmological constant or the open Universe model. Moreover in our model, the local Hubble parameter becomes larger than the global one. Thus we show that the present observed local Hubble parameter measured by using the Cepheid variables is compatible with the global Hubble parameter measured by using the Sunyaev-Zeldovich effect. We also examine the energy conditions, evolution of the effective gravitational constant and the nucleosynthesis process.