Most of the galaxies harbor supermassive Black Holes (SMBH) in their center. Some of them are observed in very high redshifts. We explore the possibility that SMBH form from the coherent waves of Bose-Einstein condensate (BEC) which are supposed to form the dark matter. We first study the isotropic and anisotropic collapses of BEC. We find the BEC wave can easily collapse to form SMBH but the realistic amount of angular momentum completely prevents the collapse. We further explore the Axion case with attractive interaction and find the moderate mass ratio between the SMBH and the dark halo around it. We also obtain the mass distribution function of BH within a single galaxy.

Super-rotations of the planetary atmosphere are reconsidered from the dynamical point of view. In particular, we emphasize that the super-rotation appears spontaneously without any explicit force. Although the super-rotation violates the bilateral symmetry (east-west reflection symmetry) of the system, this violation is spontaneous. Constructing a minimal model that derives the super-rotation, we clarify the condition for the super-rotation to appear. We find that the flow is always determined autonomously so that the flow speed becomes maximum or the temperature difference smallest. After constraining the parameters of the model from observations, we compare our model with the others most of which demand the explicit symmetry violation due to the planetary rotation.

Inflation in the early Universe is a grand phase transition which have produced the seeds of all the structures we now observe. We focus on the non-equilibrium aspect of this phase transition especially the inevitable infrared (IR) divergence associated to the the quantum and classical fields during the inflation. There is a long history of research for removing this IR divergence for healthy perturbation calculations. On the other hand, the same IR divergence is quite relevant and have developed the primordial density fluctuations in the early Universe. We develop a unified formalism in which the IR divergence is clearly separated from the microscopic quantum field theory but only appear in the statistical classical structure. We derive the classical Langevin equation for the order parameter within the quantum field theory through the instability of the de Sitter vacuum during the inflation. This separation process is relevant in general to develop macroscopic structures and to derive the basic properties of statistical mechanics in the quantum field theory.

Supermassive black holes (SMBH) of size $10^{6-10}M_{\odot}$ are common in the Universe and they define the center of the galaxies. A galaxy and the SMBH are generally thought to have co-evolved. However, the SMBH cannot evolve so fast as commonly observed even at redshift $z>6$. Therefore SMBH must form first before galaxy. Our goal is to clarify how this mature SMBH forms galaxy. Furthermore we clarify the mechanism how the SMBH designs variety of structures of galaxies. We explore a natural hypothesis that the SMBH has been formed mature at $z\approx10$ before stars and galaxies. The SMBH forms energetic jets and outflows which trigger massive star formation in the ambient gas. They eventually construct globular clusters and classical bulge as well as the body of elliptical galaxies. We propose simple models which implement these processes along with the standard $\Lambda$CDM-model. We point out that the globular clusters and classical bulges have a common origin but are in different phases. The same is true for the elliptical and spiral galaxies. Physics behind these phase division is the runaway star formation process with strong feedback to SMBH. This is similar to the forest-fire model that displays self-organized criticality. Finally we speculate several observational predictions that may help to test the present arguments.

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 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 develop our novel model of cosmology based on the Bose-Einstein condensation. This model unifies the Dark Energy and the Dark Matter, and predicts multiple collapse of condensation, followed by the final acceleration regime of cosmic expansion. We first explore the generality of this model, especially the constraints on the boson mass and condensation conditions. We further argue the robustness of this model over the wide range of parameters of mass, self coupling constant and the condensation rate. Then the dynamics of BEC collapse and the preferred scale of the collapse are studied. Finally, we describe possible observational tests of our model, especially, the periodicity of the collapses and the gravitational wave associated with them.

Flat density cores have been obtained for a limited number of clusters of galaxies by strong gravitational lensing. This paper explores the possibility that the degeneracy pressure of fermionic dark matter accounts for the flat top density profiles. This is a case study of A1689 for which the density profile has been obtained from the inner region out to 1Mpc by the combination of strong and weak lensing. In the case that the dark matter consists of the mixture of degenerate relic neutrinos and collisionless cold dark matter particles, the acceptable mass range for relic neutrinos is between 1 and 2 eV, if the ratio of the two kinds of dark matter particles is fixed to its cosmic value.

We do not know 96% of the total matter in the universe at present. In this paper, a cosmological model is proposed in which Dark Energy (DE) is identified as Bose-Einstein Condensation (BEC) of some boson field. Global cosmic acceleration caused by this BEC and multiple rapid collapses of BEC into black holes etc. (=Dark Matter (DM)) are examined based on the relativistic version of the Gross-Pitaevskii equation. We propose (a) a novel mechanism of inflation free from the slow-rolling condition, (b) a natural solution for the cosmic coincidence ('Why Now?') problem through the transition from DE into DM, (c) very early formation of highly non-linear objects such as black holes, which might trigger the first light as a form of quasars, and (d) log-z periodicity in the subsequent BEC collapsing time. All of these are based on the steady slow BEC process.

Why does the Zel'dovich approximation (ZA) work well to describe gravitational collapse in the universe? This problem is examined by focusing on its dependence on the dimensionality of the collapse. The ZA is known to be exact for a one-dimensional collapse. We show that the ZA becomes progressively more accurate in the order of three-, two-, and one-dimensional collapses. Furthermore, using models for spheroidal collapse, we show that the ZA remains accurate in all collapses, which become progressively lower dimensional with the passage of time. That is, the ZA is accurate because the essence of the gravitational collapse is incorporated in the ZA.

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.

Equilibrium configurations of weakly interacting fully degenerate fermionic dark matter are considered at various scales in the Universe. We treat the general situations for the gravity from Newtonian to general relativity and the degeneracy from nonrelativistic to relativistic. A dimensionless equilibrium configuration is specified by a single parameter regardless of particle properties, the Fermi velocity at the center, and the scalings of mass and length are specified by the rest mass and statistical weight of the dark matter particle. We focus our attention to the flat-top nature of the mass column density profile of the cluster of galaxies, A1689, recently reported by Broadhurst et al. using gravitational lensing. We convert the column density profile to a volume density profile assuming spherical symmetry and derive a 3D encircled mass profile of A1689, which is compared with the model profiles of degenerate fermion structures. The flat-top profile is reproduced. The corresponding fermion mass ranges from 2 eV to 30 eV depending on the actual scale of the degenerate structure. If massive neutrinos are the dominant dark matter, the rest mass will be about 4.7 or 2.3 eV respectively for Majorana or Dirac neutrinos. The mass and size of the degenerate structure are $10^{14}M_\odot$ and 100 kpc for Majorana neutrinos, and 5$\times10^{14}M_\odot$ and 300 kpc for Dirac neutrinos. If we identify the fermions as heavier sterile neutrinos, they yield the characteristic mass hierarchy of black holes; giant black hole at the center of a galaxy and the intermediate mass black holes. Thus we propose the possibility that the mass hierarchy of fermions determines that of black holes in the Universe.

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.

We study a cosmological model in which the boson dark matter gradually condensates into dark energy. Negative pressure associated with the condensate yields the accelerated expansion of the Universe and the rapid collapse of the smallest scale fluctuations into many black holes, which become the seeds of the first galaxies. The cycle of gradual sedimentation and rapid collapse of condensate repeats many times and self-regularizes the ratio of dark energy and dark matter to be order one.

According to N-body numerical simulations, the radial pairwise velocities of galaxies have negative skewness in the quasi-nonlinear regime. To understand its origin, we calculate the probability distribution function of the radial pairwise velocity using the Zel'dovich approximation, i.e., an analytical approximation for gravitational clustering. The calculated probability distribution function is in good agreement with the result of N-body simulations. Thus the negative skewness originates in relative motions of galaxies in the clustering process that the infall dominates over the expansion.

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.

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.

Gamma ray bursts (GRBs) are known to have short-time variability and power-law behavior with the index -1.67 in the power spectrum density. Reanalyzing the expanded data, we have found a) the power-law comes from the global profile of the burst and not from the self-similar shots nor rapid fluctuations in the luminosity profile. b) The power indices vary from burst to burst and the value -1.67 is given simply as the mean value of the distribution; there is no systematic correlation among GRBs to yield the power law.

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 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.

Among several analytic approximations for the growth of density fluctuations in the expanding Universe, Zel'dovich approximation in Lagrangian coordinate scheme is known to be unusually accurate even in mildly non-linear regime. This approximation is very similar to the Pad\'e approximation in appearance. We first establish, however, that these two are actually different and independent approximations with each other by using a model of spheroidal mass collapse. Then we propose Pad\'e-prescribed Zel'dovich-type approximations and demonstrate, within this model, that they are much accurate than any other known nonlinear approximations.

Among various analytic approximations for the growth of density fluctuations in the expanding Universe, Zel'dovich approximation and its extensions in Lagrangian scheme are known to be accurate even in mildly non-linear regime. The aim of this paper is to investigate the reason why these Zel'dovich-type approximations work accurately beyond the linear regime from the following two points of view: (1) Dimensionality of the system and (2) the Lagrangian scheme on which the Zel'dovich approximation is grounded. In order to examine the dimensionality, we introduce a model with spheroidal mass distribution. In order to examine the Lagrangian scheme, we introduce the Pad\'e approximation in Eulerian scheme. We clarify which of these aspects supports the unusual accuracy of the Zel'dovich-type approximations. We also give an implication for more accurate approximation method beyond the Zel'dovich-type approximations.

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.

Weak first-order phase transitions proceed with percolation of new phase. The kinematics of this process is clarified from the point of view of subcritical bubbles. We examine the effect of small subcritical bubbles around a large domain of asymmetric phase by introducing an effective geometry. The percolation process can be understood as a perpetual growth of the large domain aided by the small subcritical bubbles.

We derivate the Langevin and the Fokker-Planck equations for the radius of $O(3)$-symmetric subcritical bubbles as a phenomenological model to treat thermal fluctuation. The effect of thermal noise on subcritical bubbles is examined. We find that the fluctuation-dissipation relation holds and that in the high temperature phase the system settles down rapidly to the thermal equilibrium state even if it was in a nonequilibrium state initially. We then estimate the typical size of subcritical bubbles as well as the amplitude of fluctuations on that scale. We also discuss their implication to the electroweak phase transition.

We estimate the amplitude of thermal fluctuations by calculating the typical size of subcritical bubbles in cosmological electroweak phase transition and show that this thermal fluctuation effect drastically changes dynamics of the phase transition from the ordinary first order type with supercooling. From this fact, we conclude that the standard electroweak baryogenesis scenario associated with such a first order transition does not work in the minimal standard model in certain conditions.