In modern cosmology, the precision of the theoretical prediction is increasingly required. In cosmological $N$-body simulations, the effect of higher-order Lagrangian perturbation on the initial conditions appears in terms of statistical quantities of matter density field. We have considered the effect of third-order Lagrangian perturbation (3LPT) on the initial conditions, which can be applied to Gadget-2 code. Then, as statistical quantities, non-Gaussianity of matter density field has been compared between cases of different order perturbations for the initial conditions. Then, we conclude that we should apply the initial conditions with 3LPT for predictions of precision cosmology,

In this study, an $N$-body simulation code was developed for self-gravitating systems with a limited first-order post-Newtonian approximation. The code was applied to a special case in which the system consists of one massive object and many low-mass objects. Therefore, the behavior of stars around the massive black hole could be analyzed. A graphics processing unit (GPU) was used to accelerate the code execution, and it could be accelerated by several tens of times compared to a single-core CPU for $N \simeq 10^4$ objects.

We study the initial conditions for cosmological $N$-body simulations for precision cosmology. In general, Zel'dovich approximation has been applied for the initial conditions of $N$-body simulations for a long time. These initial conditions provide incorrect higher-order growth. These error caused by setting up the initial conditions by perturbation theory is called transients. We investigated the impact of transient on non-Gaussianity of density field by performing cosmological $N$-body simulations with initial conditions based on first-, second-, and third-order Lagrangian perturbation theory in previous paper. In this paper, we evaluates the effect of the transverse mode in the third-order Lagrangian perturbation theory for several statistical quantities such as power spectrum, non-Gaussianty, and multifractal dimensions. Then we clarified that the effect of the transverse mode in the third-order Lagrangian perturbation theory is quite small.

We have derived fourth-order perturbative equations in Lagrangian perturbation theory for a cosmological dust fluid. These equations are derived under the supposition of Newtonian cosmology in the Friedmann-Lema\^{i}tre-Robertson-Walker Universe model. Even if we consider the longitudinal mode in the first-order perturbation, the transverse mode appears in the third-order perturbation. Furthermore, in this case, six longitudinal-mode equations and four transverse-mode equations appear in the fourth-order perturbation. The application of the fourth-order perturbation leads to a precise prediction of the large-scale structure.

We will construct a theory which can explain the dynamics toward the steady state self-gravitating systems (SGSs) where many particles interact via the gravitational force. Real examples of SGS in the universe are globular clusters and galaxies. The idea is to represent an interaction by which a particle of the system is affected from the others by a special random force. That is, we will use a special Langevin equation, just as the normal Langevin equation can unveil the dynamics toward the steady state described by the Maxwell-Boltzmann distribution. However, we cannot introduce the randomness into the system without any evidence. Then, we must confirm that each orbit is random indeed. Of course, it is impossible to understand orbits of stars in globular clusters from observations. Thus we use numerical simulations. From the numerical simulations of SGS, grounds that we use the random noise become clear. The special Langevin equation includes the additive and the multiplicative noise. By using the random process, we derive the non-Maxwellian distribution of SGS especially around the core. The number density can be obtained through the steady state solution of the Fokker-Planck equation corresponding to the random process. We exhibit that the number density becomes equal to the density profiles around the core by adjusting the friction coefficient and the intensity of the multiplicative noise. Moreover, we also show that our model can be applied in the system which has a heavier particle, corresponding to the black hole in a globular cluster.

We derive the non-Maxwellian distribution of self-gravitating $N$-body systems around the core by a model based on the random process with the additive and the multiplicative noise. The number density can be obtained through the steady state solution of the Fokker-Planck equation corresponding to the random process. We exhibit that the number density becomes equal to that of the King model around the core by adjusting the friction coefficient and the intensity of the multiplicative noise. We also show that our model can be applied in the system which has a heavier particle. Moreover, we confirm the validity of our model by comparing with our numerical simulation.

We study second-order cosmological perturbations in scalar-tensor models of dark energy that satisfy local gravity constraints, including f(R) gravity. We derive equations for matter fluctuations under a sub-horizon approximation and clarify conditions under which first-order perturbations in the scalar field can be neglected relative to second-order matter and velocity perturbations. We also compute the skewness of the matter density distribution and find that the difference from the LCDM model is only less than a few percent even if the growth rate of first-order perturbations is significantly different from that in the LCDM model. This shows that the skewness provides a model-independent test for the picture of gravitational instability from Gaussian initial perturbations including scalar-tensor modified gravity models.

We study the effect of modified gravity on weak lensing in a class of scalar-tensor theory that includes $f(R)$ gravity as a special case. These models are designed to satisfy local gravity constraints by having a large scalar-field mass in a region of high curvature. Matter density perturbations in these models are enhanced at small redshifts because of the presence of a coupling $Q$ that characterizes the strength between dark energy and non-relativistic matter. We compute a convergence power spectrum of weak lensing numerically and show that the spectral index and the amplitude of the spectrum in the linear regime can be significantly modified compared to the $\Lambda$CDM model for large values of $|Q|$ of the order of unity. Thus weak lensing provides a powerful tool to constrain such large coupling scalar-tensor models including $f(R)$ gravity.

We explore the initial conditions for cosmological N-body simulations suitable for calculating the skewness and kurtosis of the density field. In general, the initial conditions based on the perturbation theory (PT) provide incorrect second-order and higher-order growth. These errors implied by the use of the perturbation theory to set up the initial conditions in N-body simulations are called transients. Unless these transients are completely suppressed compared with the dominant growing mode, we can not reproduce the correct evolution of cumulants with orders higher than two, even though there is no problem with the numerical scheme. We investigate the impact of transients on the observable statistical quantities by performing $N$-body simulations with initial conditions based on Lagrangian perturbation theory (LPT). We show that the effects of transients on the kurtosis from the initial conditions, based on second-order Lagrangian perturbation theory (2LPT) have almost disappeared by $z\sim5$, as long as the initial conditions are set at $z > 30$. This means that for practical purposes, the initial conditions based on 2LPT are accurate enough for numerical calculations of skewness and kurtosis.

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.

We explore the possibility of putting constraints on dark energy models with statistical property of large scale structure in the non-linear region. In particular, we investigate the $w$ dependence of non-Gaussianity of the smoothed density distribution generated by the nonlinear dynamics. In order to follow the non-linear evolution of the density fluctuations, we apply N-body simulations based on $P^3 M$ scheme. We show that the relative difference between non-Gaussianity of $w=-0.8$ model and that of $w=-1.0$ model is $0.67 %$ (skewness) and $1.2 %$ (kurtosis) for $R=8h^{-1}$ Mpc. We also calculate the correspondent quantities for $R=2h^{-1}$ Mpc, $3.0 %$ (skewness) and $4.5 %$ (kurtosis), and the difference turn out to be greater, even though non-linearity in this scale is so strong that the complex physical processes about galaxy formation affect the galaxy distribution. From this, we can expect that the difference can be tested by all sky galaxy surveys with the help of mock catalogs created by selection functions, which suggests that non-Gaussianity of the density distribution potentially plays an important role for extracting information on dark energy.

We study the behavior of Lagrangian perturbative solutions. For a spherical void model, the higher order the Lagrangian perturbation we consider, the worse the approximation becomes in late-time evolution. In particular, if we stop to improve until an even order is reached, the perturbative solution describes the contraction of the void. To solve this problem, we consider improving the perturbative solution using Shanks transformation, which accelerates the convergence of the sequence. After the transformation, we find that the accuracy of higher-order perturbation is recovered and the perturbative solution is refined well. Then we show that this improvement method can apply for a $\Lambda$CDM model and improved the power spectrum of the density field.

According to recent observations, the existence of the dark energy has been considered. Even though we have obtained the constraint of the equation of the state for dark energy ($p = w \rho$) as $-1 \le w \le -0.78$ by combining WMAP data with other astronomical data, in order to pin down $w$, it is necessary to use other independent observational tools. For this purpose, we consider the $w$ dependence of the non-Gaussianity of the density distribution generated by nonlinear dynamics. To extract the non-Gaussianity, we follow a semi-analytic approach based on Lagrangian linear perturbation theory, which provides an accurate value for the quasi-nonlinear region. From our results, the difference of the non-Gaussianity between $w = -1$ and $w= -0.5$ is about 4% while that between $w = -1$ and $w= -0.8$ is about $0.9 %$. For the highly non-linear region, we estimate the difference by combining this perturbative approach with N-body simulation executed for our previous paper. From this, we can expect the difference to be more enhanced in the low-$z$ region, which suggests that the non-Gaussianity of the density distribution potentially plays an important role for extracting the information of dark energy.

For the evolution of density fluctuation in nonlinear cosmological dynamics, adhesion approximation (AA) is proposed as a phenomenological model, which is especially useful for describing nonlinear evolution. However, the origin of the artificial viscosity in AA is not clarified. Recently, Buchert and Dom\'{\i}nguez report if the velocity dispersion of the dust fluid is regarded as isotropic, it works on a principle similar to viscosity or effective pressure, and they consider isotropic velocity dispersion as the origin of the artificial viscosity in AA. They name their model the Euler-Jeans-Newton (EJN) model. In this paper, we focus on the velocity distribution in AA and the EJN model and examine the time evolution in both models. We find the behavior of AA differs from that of the EJN model, i.e., although the peculiar velocity in the EJN model oscillates, that in AA is monotonically decelerated due to viscosity without oscillation. Therefore it is hard to regard viscosity in AA as effective pressure in the EJN model.

Lagrangian perturbation theory for cosmological fluid describes structure formation in the quasi-nonlinear stage well. In a previous paper, we presented a third-order perturbative equation for Lagrangian perturbation with pressure. There we considered only the longitudinal modes for the first-order perturbation. In this paper, we generalize the perturbation, i.e., we consider both the longitudinal and the transverse modes for the first-order perturbation. Then we derive third-order perturbative equations and solutions.

We discuss various analytical approximation methods for the evolution of the density fluctuation in the Universe. From primordial density fluctuation, the large-scale structure is formed via its own self-gravitational instability. For this dynamical evolution, several approaches have been considered. In Newtonian cosmology, in which matter motion is described by Newtonian dynamics with cosmic expansion, Lagrangian description for the cosmic fluid is known to work rather well for quasi-nonlinear clustering regime. In this paper, we briefly review Lagrangian perturbation theory in Newtonian cosmology.

We study the one-point probability distribution functions (PDFs) of the peculiar velocity and the density fluctuation in a cosmological fluid. Within the perturbative approach to the structure formation scenario, the effect of ``pressure'' has recently been an area of research interest. The velocity dispersion of the cosmological fluid creates effective ``pressure'' or viscosity terms. From this viewpoint, because the pressure reflects a nonlinear effect of the motion of the fluid, the pressure model would include nonlinear effects. Here we analyze the Lagrangian linear perturbation PDFs for both the Zel'dovich approximation and the pressure model. We find that the PDFs of the peculiar velocity remain Gaussian, even if we consider the pressure. For the PDFs of the density fluctuation, the occurrence of non-Gaussianity depends on the ``equation of state'' for the fluid. Therefore we distinguish the ``equation of state'' using the PDFs.

Lagrangian perturbation theory for cosmological fluid describes structure formation in the quasi-nonlinear stage well. We present a third-order perturbative equation for Lagrangian perturbation with pressure in both the longitudinal and transverse modes. Then we derive the perturbative solution for simplest case.

We present the phase diagram, in both the microcanonical and the canonical ensemble, of the Self-Gravitating-Ring (SGR) model, which describes the motion of equal point masses constrained on a ring and subject to 3D gravitational attraction. If the interaction is regularized at short distances by the introduction of a softening parameter, a global entropy maximum always exists, and thermodynamics is well defined in the mean-field limit. However, ensembles are not equivalent and a phase of negative specific heat in the microcanonical ensemble appears in a wide intermediate energy region, if the softening parameter is small enough. The phase transition changes from second to first order at a tricritical point, whose location is not the same in the two ensembles. All these features make of the SGR model the best prototype of a self-gravitating system in one dimension. In order to obtain the stable stationary mass distribution, we apply a new iterative method, inspired by a previous one used in 2D turbulence, which ensures entropy increase and, hence, convergence towards an equilibrium state.

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.

Basing our discussion on the Lagrangian description of hydrodynamics, we studied the evolution of density fluctuation for nonlinear cosmological dynamics. Adhesion approximation (AA) is known as a phenomenological model that describes the nonlinear evolution of density fluctuation rather well and that does not form a caustic. In addition to this model, we have benefited from discussion of the relation between artificial viscosity in AA and velocity dispersion. Moreover, we found it useful to regard whether the velocity dispersion is isotropic produces effective `pressure' or viscosity terms. In this paper, we analyze plane- and spherical-symmetric cases and compare AA with Lagrangian models where pressure is given by a polytropic equation of state. From our analyses, the pressure model undergoes evolution similar to that of AA until reaching a quasi-nonlinear regime. Compared with the results of a numerical calculation, the linear approximation of the pressure model seems rather good until a quasi-nonlinear regime develops. However, because of oscillation arising from the Jeans instability, we could not produce a stable nonlinear structure.

We analyzed the performance of a perturbation theory for nonlinear cosmological dynamics, based on the Lagrangian description of hydrodynamics. In our previous paper, we solved hydrodynamic equations for a self-gravitating fluid with pressure, given by a polytropic equation of state, using a perturbation method. Then we obtained the first-order solutions in generic background universes and the second-order solutions for a wider range of polytrope exponents. Using these results, we describe density fields with scale-free spectrum, SCDM, and LCDM models. Then we analyze cross-correlation coefficient of the density field between N-body simulation and Lagrangian linear perturbation theory, and the probability distribution of density fluctuation. From our analyses, for scale-free spectrum models, the case of the polytrope exponent 5/3 shows better performance than the Zel'dovich approximation and the truncated Zel'dovich approximation in a quasi-nonlinear regime. On the other hand, for SCDM and LCDM models, improvement by the effect of the velocity dispersion was small.

We extensively develop a perturbation theory for nonlinear cosmological dynamics, based on the Lagrangian description of hydrodynamics. We solve hydrodynamic equations for a self-gravitating fluid with pressure, given by a polytropic equation of state, using a perturbation method up to second order. This perturbative approach is an extension of the usual Lagrangian perturbation theory for a pressureless fluid, in view of inclusion of the pressure effect, which should be taken into account on the occurrence of velocity dispersion. We obtain the first-order solutions in generic background universes and the second-order solutions in wider range of a polytropic index, whereas our previous work gives the first-order solutions only in the Einstein-de Sitter background and the second-order solutions for the polytropic index 4/3. Using the perturbation solutions, we present illustrative examples of our formulation in one- and two-dimensional systems, and discuss how the evolution of inhomogeneities changes for the variation of the polytropic index.

We study an effect of inhomogeneity of density distribution of the Universe. We propose a new Lagrangian perturbation theory with a backreaction effect by inhomogeneity. The inhomogeneity affects the expansion rate in a local domain and its own growing rate. We numerically analyze a one-dimensional plane-symmetric model, and calculate the probability distribution functions (PDFs) of several observed variables to discuss those statistical properties. We find that the PDF of pairwise peculiar velocity shows an effective difference from the conventional Lagrangian approach, i.e. even in one-dimensional plane symmetric case, the PDF approaches an exponential form in a small relative-velocity region, which agree with the N-body simulation.

We formulate a perturbative approximation to gravitational instability, based on Lagrangian hydrodynamics in Newtonian cosmology. We take account of `pressure' effect of fluid, which is kinematically caused by velocity dispersion, to aim hydrodynamical description beyond shell crossing. Master equations in the Lagrangian description are derived and solved perturbatively up to second order. Then, as an illustration, power spectra of density fluctuations are computed in a one-dimensional model from the Lagrangian approximations and Eulerian linear perturbation theory for comparison. We find that the results by the Lagrangian approximations are different from those by the Eulerian one in weakly non-linear regime at the scales smaller than the Jeans length. We also show the validity of the perturbative Lagrangian approximations by consulting difference between the first-order and the second-order approximations.

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.

Two-point correlation function of galaxy distribution shows that the structure in the present Universe is scale-free up to a certain scale (at least several tens Mpc), which suggests that a fractal structure may exist. If small primordial density fluctuations have a fractal structure, the present fractal-like nonlinear structure below the horizon scale could be naturally explained. We analyze the time evolution of fractal density perturbations in Einstein-de Sitter universe, and study how the perturbation evolves and what kind of nonlinear structure will come out. We assume a one-dimensional collisionless sheet model with initial Cantor-type fractal perturbations. The nonlinear structure seems to approach some attractor with a unique fractal dimension, which is independent of the fractal dimensions of initial perturbations. A discrete self-similarity in the phase space is also found when the universal nonlinear fractal structure is reached.

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.