A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it is not clearly seen in lower metallicity region ([Fe/H] < -0.25). By using a large number of chemo-dynamical 2D test-particle simulations with a rotating bar and/or spiral arms, we find that the observed [Fe/H] dependence of the Hercules stream is a natural consequence of the inside-out formation of the stellar disc and the existence of highly non-closed orbits in the rotating frame of the bar or spiral arms. Our successful models that reproduce the observed properties of the Hercules stream include not only fast-bar-only and fast-bar+spiral models, but also slow-bar+spiral models. This indicates that it is very difficult to estimate the pattern speed of the bar or spiral arms based only on the observations of the Hercules stream in the Solar neighbourhood. As a by-product of our simulations, we make some predictions about the locations across the Galactic plane where we can observe velocity bimodality that is not associated with the Hercules stream. These predictions can be tested by the Gaia Data Release 2, and such a test will improve our understanding of the evolution of the Milky Way stellar disc.

We develop a method for identifying a compact object in binary systems with astrometric measurements and apply it to some binaries. Compact objects in some high-mass X-ray binaries and gamma-ray binaries are unknown, which is responsible for the fact that emission mechanisms in such systems have not yet confirmed. The accurate estimate of the mass of the compact object allows us to identify the compact object in such systems. Astrometric measurements are expected to enable us to estimate the masses of the compact objects in the binary systems via a determination of a binary orbit. We aim to evaluate the possibility of the identification of the compact objects for some binary systems. We then calculate probabilities that the compact object is correctly identified with astrometric observation (= confidence level) by taking into account a dependence of the orbital shape on orbital parameters and distributions of masses of white dwarfs, neutron stars, and black holes. We find that the astrometric measurements with the precision of 70 $\mu$as for $\gamma$~Cas allows us to identify the compact object at 99% confidence level if the compact object is a white dwarf with 0.6 M$_{\odot}$. In addition, we can identify the compact object with the precision of 10 $\mu$as at 97% or larger confidence level for LS~I~+61$^{\circ}$~303 and 99% or larger for HESS~J0632+057. These results imply that the astrometric measurements with the 10-$\mu$as precision level can realize the identification of compact objects for $\gamma$~Cas, LS~I~+61$^{\circ}$~303, and HESS~J0632+057.

We investigate how accurately phase space distribution functions (DFs) in galactic models can be reconstructed by a made-to-measure (M2M) method, which constructs $N$-particle models of stellar systems from photometric and various kinematic data. The advantage of the M2M method is that this method can be applied to various galactic models without assumption of the spatial symmetries of gravitational potentials adopted in galactic models, and furthermore, numerical calculations of the orbits of the stars cannot be severely constrained by the capacities of computer memories. The M2M method has been applied to various galactic models. However, the degree of accuracy for the recovery of DFs derived by the M2M method in galactic models has never been investigated carefully. Therefore, we show the degree of accuracy for the recovery of the DFs for the anisotropic Plummer model and the axisymmetric St\"{a}ckel model, which have analytic solutions of the DFs. Furthermore, this study provides the dependence of the degree of accuracy for the recovery of the DFs on various parameters and a procedure adopted in this paper. As a result, we find that the degree of accuracy for the recovery of the DFs derived by the M2M method for the spherical target model is a few percent, and more than ten percent for the axisymmetric target model.

We present post-Newtonian $N$-body simulations on mergers of accreting stellar-mass black holes (BHs), where such general relativistic effects as the pericenter shift and gravitational wave (GW) emission are taken into consideration. The attention is concentrated on the effects of the dynamical friction and the Hoyle-Lyttleton mass accretion by ambient gas. We consider a system composed of ten BHs with initial mass of $30~M_\odot$. As a result, we show that mergers of accreting stellar-mass BHs are classified into four types: a gas drag-driven, an interplay-driven, a three body-driven, or an accretion-driven merger. We find that BH mergers proceed before significant mass accretion, even if the accretion rate is $\sim10$ Eddington accretion rate, and then all BHs can merge into one heavy BH. Using the simulation results for a wide range of parameters, we derive a critical accretion rate ($\dot{m}_{\rm c}$), below which the BH growth is promoted faster by mergers. Also, it is found that the effect of the recoil by the GW emission can reduce $\dot{m}_{\rm c}$ especially in gas number density higher than $10^8~{\rm cm}^{-3}$, and enhance the escape probability of merged BHs. Very recently, a gravitational wave event, GW150914, as a result of the merger of a $\sim 30~M_\odot$ BH binary has been detected (Abbott et al. 2016). Based on the present simulations, the BH merger in GW150914 is likely to be driven by three-body encounters accompanied by a few $M_\odot$ of gas accretion, in high-density environments like dense interstellar clouds or galactic nuclei.

We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitational wave emission. As a result, we find that multiple BHs are able to merge into one BH within 100 Myr in a wide range of BH density. The merger mechanism is revealed to be categorized into three types: gas-drag-driven merger (type A), interplay-driven merger (type B), and three-body-driven merger (type C). We find the relation between the merger mechanism and the ratio of the gas mass within the initial BH orbit ($M_{gas}$) to the total BH mass (${\Sigma}M_{BH}$). Type A merger occurs if $M_{gas} \gtrsim 10^5 {\Sigma}M_{BH}$, type B if $M_{gas} \lesssim 10^5 {\Sigma}M_{BH}$, and type C if $M_{gas} \ll 10^5 {\Sigma}M_{BH}$. Supposing the gas and BH density based on the recent numerical simulations on first stars, all the BH remnants from first stars are likely to merge into one BH through the type B or C mechanism. Also, we find that multiple massive BHs ($M_{BH} = 10^4 M_{\odot}$) distributed over several parsec can merge into one BH through the type B mechanism, if the gas density is higher than $5\times 10^6$ cm$^{-3}$. The present results imply that the BH merger may contribute significantly to the formation of supermassive BHs at high redshift epochs.

A moment approach for orbit determinations of astrometric binaries from astrometric observations alone has been recently studied for a low signal-to-noise ratio (Iwama et al. 2013, PASJ, 65, 2). With avoiding a direct use of the time-consuming Kepler equation, temporal information is taken into account to increase the accuracy of statistical moments. As numerical tests, 100 realizations are done and the mean and the standard deviation are also evaluated. For a semi-major axis, the difference between the mean of the recovered values and the true value decreases to less than a tenth in the case of $10000$ observed points. Therefore, the present moment approach works better than the previous one for the orbit determinations when one has a number of the observed points. The present approach is thus applicable to Cyg X-1.

A phase-space distribution function of the steady state in galaxy models that admits regular orbits overall in the phase-space can be represented by a function of three action variables. This type of distribution function in Galactic models is often constructed theoretically for comparison of the Galactic models with observational data as a test of the models. On the other hand, observations give Cartesian phase-space coordinates of stars. Therefore it is necessary to relate action variables and Cartesian coordinates in investigating whether the distribution function constructed in galaxy models can explain observational data. Generating functions are very useful in practice for this purpose, because calculations of relations between action variables and Cartesian coordinates by generating functions do not require a lot of computational time or computer memory in comparison with direct numerical integration calculations of stellar orbits. Here, we propose a new method called a torus-fitting method, by which a generating function is derived numerically for models of the Galactic potential in which almost all orbits are regular. We confirmed the torus-fitting method can be applied to major orbit families (box and loop orbits) in some two-dimensional potentials. Furthermore, the torus-fitting method is still applicable to resonant orbit families, besides major orbit families. Hence the torus-fitting method is useful for analyzing real Galactic systems in which a lot of resonant orbit families might exist.

To determine the local dark matter density (LDMD) of the solar system is a classical problem in astronomy. Recently, a novel method of determining the LDMD from stellar distribution and vertical velocity dispersion profiles perpendicular to the Galactic plane was devised. This method has the advantage of abolishing conventional approximations and using only a few assumptions. Our aims are to carefully scrutinize this method and to examine influences by uncertainties of astrometric observations. We discuss how the determinations of the LDMD vary with observational precisions on parallax, proper motion, and line-of-sight velocity measurements. To examine the influences by the observational imprecision, we created mock observation data for stars that are dynamical tracers based on an analytical galaxy model and applied parametrized observational errors to the mock data. We evaluated the accuracy of determining the LDMD by applying the method to the mock data. In addition, we estimated a sample size and observational precision required to determine the LDMD with accuracy. We find that the method is capable of determining the LDMD with accuracy if the sample size and observational precisions are satisfactory. The random errors of parallaxes and proper motions can cause systematic overestimation of the LDMD. We estimate the required precisions of the parallax measurements to be approximately 0.1-0.3 milliarcseconds at 1 kpc away from the Sun; the proper motion precisions do not seem to be as important as the parallaxes. From these results, we expect that using the Hipparcos catalog would overestimate the LDMD because of the imprecise parallax measurements if this method is applied; however, we emphasize the capability of the method. We expect that Gaia will provide data precise enough to determine the LDMD.

We investigate the properties of HI-rich galaxies detected in blind radio surveys within the hierarchical structure formation scenario using a semi-analytic model of galaxy formation. By drawing a detailed comparison between the properties of HI-selected galaxies and HI absorption systems, we argue a link between the local galaxy population and quasar absorption systems, particularly for Damped Ly-alpha absorption (DLA) systems and sub-DLA systems. First, we evaluate how many HI-selected galaxies exhibit HI column densities as high as those of DLA systems. We find that HI-selected galaxies with HI masses M(HI) > 10^8 solar masses have gaseous disks that produce HI column densities comparable to those of DLA systems. We conclude that DLA galaxies where the HI column densities are as high as those of DLA systems, contribute significantly to the population of HI-selected galaxies at M(HI) > 10^8 solar masses. Second, we find that star formation rates (SFRs) correlate tightly with HI masses rather than B- (and J-) band luminosities. In the low-mass range M(HI) < 10^8 solar masses, sub-DLA galaxies replace DLA galaxies as the dominant population. The number fraction of sub-DLA galaxies relative to galaxies reaches 40%-60% at HI masses 10^8 solar masses and 30%-80% at 10^7 solar masses. The HI-selected galaxies at 10^7 solar masses are a strong probe of sub-DLA systems that place stringent constraints on galaxy formation and evolution.

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 construct the Numerical Galaxy Catalog ($\nu$GC), based on a semi-analytic model of galaxy formation combined with high-resolution N-body simulations in a $\Lambda$-dominated flat cold dark matter ($\Lambda$CDM) cosmological model. The model includes several essential ingredients for galaxy formation, such as merging histories of dark halos directly taken from N-body simulations, radiative gas cooling, star formation, heating by supernova explosions (supernova feedback), mergers of galaxies, population synthesis, and extinction by internal dust and intervening HI clouds. As the first paper in a series using this model, we focus on basic photometric, structural and kinematical properties of galaxies at present and high redshifts. Two sets of model parameters are examined, strong and weak supernova feedback models, which are in good agreement with observational luminosity functions of local galaxies in a range of observational uncertainty. Both models agree well with many observations such as cold gas mass-to-stellar luminosity ratios of spiral galaxies, HI mass functions, galaxy sizes, faint galaxy number counts and photometric redshift distributions in optical pass-bands, isophotal angular sizes, and cosmic star formation rates. In particular, the strong supernova feedback model is in much better agreement with near-infrared (K'-band) faint galaxy number counts and redshift distribution than the weak feedback model and our previous semi-analytic models based on the extended Press-Schechter formalism. (Abridged)

We investigate chemical enrichment in Damped Lyman alpha (DLA) systems in the hierarchical structure formation scenario using a semi-analytic model of galaxy formation. The model developed by Nagashima, Totani, Gouda and Yoshii takes into account various selection effects on high-redshift galaxies and can show fundamental observational properties of galaxies, such as luminosity functions and number-magnitude/redshift relations. DLA systems offer the possibilities of measuring metal abundance more accurately than faint galaxies. For example, recent measurements of zinc abundance can provide good evidence for understanding the processes of metal pollution and star formation in DLA systems because zinc is virtually unaffected by dust depletion. Here we focus on this advantage for observation in order to explore the metallicity evolution in DLA systems at high redshifts. We can consistently show the metallicity evolution for reasonable models which also reproduce fundamental properties of local galaxy population. This result suggests that the chemical evolution of DLA systems can be consistently reconciled with the observational features of typical galaxies. We also investigate other properties of DLA systems (column density distribution and mass density of cold gas), and find that star formation in massive galaxies should be more active than that in low-mass ones. This is consistent with the results by Nagashima et al. and Cole et al. in which the star formation timescale is set by reproducing cold gas mass fraction in local spiral galaxies. Finally we discuss host galaxies associated with DLA systems. We conclude that they primarily consist of sub-L* and/or dwarf galaxies from the observations.

We numerically investigate the long--time evolution of density perturbations after the first appearance of caustics in an expanding cosmological model with one--dimensional `single--wave' initial conditions. Focussing on the time--intervals of caustic appearances and the spatial distribution of caustics at subsequent times, we find that the time--intervals of caustic appearances approach a constant, i.e., their time--subsequent ratio converges to 1; it is also found that the spatial distribution of caustics at a given time features some universality rules, e.g., the ratio between the position of the nearest caustic from the center and that of the second nearest caustic from the center approaches a constant. Furthermore we find some rules for the mass distribution for each caustic. Using these universality constants we are in the position to predict the spatial distribution of caustics at an arbitrary time in order to give an estimate for the power spectral index in the fully--developed non--dissipative turbulent (`virialized') regime.

We investigate the environments of quasars such as number distribution of galaxies using a semi-analytic model which includes both galaxy and quasar formations based on the hierarchical clustering scenario. We assume that a supermassive black hole is fueled by accretion of cold gas and that it is a source of quasar activity during a major merger of the quasar host galaxy with another galaxy. This major merger causes spheroid formation of the host galaxy. Our model can reproduce not only general form of the galaxy luminosity functions in the local Universe but also the observed relation between a supermassive black hole mass and a spheroid luminosity, the present black hole mass function and the quasar luminosity functions at different redshifts. Using this model, we predict the mean number of quasars per halo, bias parameter of quasars and the probability distribution of the number of galaxies around quasars. In our model, analysis of the mean number of quasars per halo shows that the spatial distribution of galaxies is different from that of quasars. Furthermore, we found from calculation of the probability distribution of galaxy numbers that at $0.2 \lesssim z \lesssim 0.5$, most quasars are likely to reside in galaxy groups. On the other hand, at $1 \lesssim z \lesssim 2$ most quasars seem to reside in more varied environments than at a lower redshift; quasars reside in environments ranging from small groups of galaxies to clusters of galaxies. Comparing these predictions with observations in future will enable us to constrain our quasar formation model.

Number counts of galaxies are re-analyzed using a semi-analytic model (SAM) of galaxy formation based on the hierarchical clustering scenario. Faint galaxies in the Subaru Deep Field (SDF) and the Hubble Deep Field (HDF) are compared with our model galaxies. We have determined the astrophysical parameters in the SAM that reproduce observations of nearby galaxies, and used them to predict the number counts and redshifts of faint galaxies for three cosmological models, the standard cold dark matter (CDM) universe, a flat lambda-CDM, and an open CDM. The novelty of our SAM analysis is the inclusion of selection effects arising from the cosmological dimming of surface brightness of high-z galaxies, and from the absorption of visible light by internal dust and intergalactic HI clouds. As was found in our previous work, in which the UV/optical HDF galaxies were compared with our model galaxies, we find that our SAM reproduces counts of near-IR SDF galaxies in low-density models, and that the standard CDM universe is not preferred, as suggested by other recent studies. Moreover, we find that simple prescriptions for (1) the timescale of star formation being proportional to the dynamical time scale of the formation of galactic disks, (2) the size of galactic disks being rotationally supported with the same specific angular momentum as that of surrounding dark halo, and (3) the dust optical depth being proportional to the metallicity of cold gas, cannot completely explain all of observed data. Improved prescriptions incorporating mild z-dependence for those are suggested from our SAM analysis.

We examined effects of the ultraviolet background radiation (UVB) on the colour--magnitude relation (CMR) of elliptical galaxies in clusters of galaxies in the hierarchical clustering scenario by using a semi-analytic model of galaxy formation. In our model the UVB photoionizes gas in dark haloes and suppresses the cooling of the diffuse hot gas onto galaxy discs. By using a semi-analytic model without the effect of the UVB, Kauffmann & Charlot found that the CMR can be reproduced by strong supernova heating because such supernova feedback suppresses the chemical enrichment in galaxies especially for small galaxies. We find that the CMR also becomes bluer because of the UVB, in a different way from the effect of supernova feedback. While the supernova feedback suppresses the chemical enrichment by a similar mechanism to galactic wind, the UVB suppresses the cooling of the hot gas. This fact induces the suppression of the metallicity of the intracluster medium (ICM). In our model we find that the existence of the UVB can plausibly account for an observed ICM metallicity that is equal to nearly 0.3 times the solar value, and that in this case we can reproduce the CMR and the metallicity of the ICM simultaneously.

Number counts of galaxies are re-analyzed using a semi-analytic model (SAM) of galaxy formation based on the hierarchical clustering scenario. We have determined the astrophysical parameters in the SAM that reproduce observations of nearby galaxies, and used them to predict the number counts and redshifts of faint galaxies for three cosmological models for (1) the standard cold dark matter (CDM) universe, (2) a low-density flat universe with nonzero cosmological constant, and (3) a low-density open universe with zero cosmological constant. The novelty of our SAM analysis is the inclusion of selection effects arising from the cosmological dimming of surface brightness of high-redshift galaxies, and also from the absorption of visible light by internal dust and intergalactic \ion{H}{1} clouds. Contrary to previous SAM analyses which do not take into account such selection effects, we find, from comparison with observed counts and redshifts of faint galaxies in the Hubble Deep Field (HDF), that the standard CDM universe is {\it not} preferred, and a low-density universe either with or without cosmological constant is favorable, as suggested by other recent studies. Moreover, we find that a simple prescription for the time scale of star formation (SF), being proportional to the dynamical time scale of the formation of the galactic disk, is unable to reproduce the observed number- redshift relation for HDF galaxies, and that the SF time scale should be nearly independent of redshift, as suggested by other SAM analyses for the formation of quasars and the evolution of damped Ly-$\alpha$ systems.

We present a method, based on the kinematic analysis of the Galactic disk stars, to clarify whether the internal motions of the stellar system in spiral arms follow those expected in the density wave theory. The method relies on the comparison with the linear relation between the phases of spatial positions and epicyclic motions of stars, as drawn from the theory. The application of the method to the 78 Galactic Cepheids near the Sun, for which accurate proper motions are available from the {\it Hipparcos} Catalogue, has revealed that these Cepheids hold no correlation between both phases, thereby implying that their motions are in contradiction with the theoretical predictions. Possible reasons for this discrepancy are discussed and future prospects are outlined.

In order to constrain parameters in galaxy formation theories, especially those for a star formation process, we investigate cold gas in elliptical galaxies. We calculate the detection rate of cold gas in them using a semi-analytic model of galaxy formation and compare it with observations. We show that the model with a long star formation time-scale (~20 Gyr) is inconsistent with observations. Thus, some mechanisms of reducing the mass of interstellar medium, such as the consumption of molecular gas by star formation and/or reheating from supernovae, are certainly effective in galaxies. Our model predicts that star formation induced when galaxies in a halo collide each other reduces the cold gas left until the present. However, we find that the reduction through random collisions of satellite (non-central) galaxies in mean free time-scale in a halo is not required to explain the observations. This may imply that the collisions and mergers between satellite galaxies do not occur so often in clusters or that they do not stimulate the star formation activity as much as the simple collision model we adopted. For cD galaxies, the predicted detection rate of cold gas is consistent with observations as long as the transformation of hot gas into cold gas is prevented in halos whose circular velocities are larger than 500 km s^-1. Moreover, we find that the cold gas brought into cDs through captures of gas-rich galaxies is little. We also show that the fraction of galaxies with observable cold gas should be small for cluster ellipticals in comparison with that for field ellipticals.

We have investigated the relation between the two-point spatial correlation function and the density profile of the dark matter halo in the strongly non-linear regime. It is well known that when the density fluctuation grows into dark matter halo whose density profile is $\rho \propto r^{-\epsilon}$(${3/2}<\epsilon<3$) on almost all mass scales, the two-point spatial correlation function obeys a power law with the power index $\gamma = 2\epsilon -3$ in the strongly non-linear regime. We find the form of the two-point spatial correlation function, which does not obey the power law when the power index $\epsilon$ is smaller than 3/2, such as the density profile $\rho \propto r^{-1}$ around the center of the halo which is proposed by Navarro, Frenk & White (1996,1997). By using the BBGKY equation in the strongly non-linear regime, it is also found that velocity parameter $h \equiv - < v > / \dot{a}x$ is not a constant even in the strongly non-linear regime ($\tilde{x} \equiv x/x_{nl} \to 0$) although it is a constant when $\epsilon > 3/2$ and then the two-point spatial correlation function can be regarded as the power law. The velocity parameter $h$ becomes 0 at the non-linear limit of $\tilde{x} \to 0$, that is, the stable clustering hypothesis cannot be satisfied when $\epsilon < 3/2$.

We investigate the effects of the UV background radiation on galaxy formation by using the semi-analytic model including the photoionization process. The semi-analytic model is based on Cole et al. (1994) and we use almost the same parameters of their `fiducial' model. We find that the UV background mainly affects the formation of dwarf galaxies. Because of the suppression of star formation, the number density of small objects corresponding to dwarf galaxies decreases compared to the case of no UV radiation when the UV background exists until the present epoch. When the UV background vanishes at a low redshift, the number density of small objects is hardly changed but the colour becomes bluer, compared to the case of no UV radiation, because stars are newly formed after the UV background vanishes. On the other hand, the UV radiation hardly affects massive galaxies. This is because the massive galaxies are formed by mergers of small galaxies.

We examined the physical origin of the colour-magnitude relation of elliptical galaxies in clusters of galaxies in the hierarchical clustering scenario by using a semi-analytic model of galaxy formation. We particularly investigate the effects of the star formation process, the supernova feedback and the UV background on the formation and evolution of elliptical galaxies. It is suggested by observation that the colour-magnitude relation reflects the metallicity-luminosity relation of galaxies, and that for reproducing the colour-magnitude relation it is important to suppress the chemical enrichment of stars in faint galaxies. For suppressing the chemical enrichment, we propose some solutions. One is assuming a strong supernova feedback. Another is introducing the effect of the UV background radiation. These solutions degenerate, but they may be distinguished by studying other properties of galaxies, like the luminosity function, colour distribution, etc. Moreover, we find that the supernova feedback model in the starburst which make the bulge component is important when the star formation time-scale in disc is longer than the age of the Universe. On the other hand, when the star formation time-scale is short, the resulting colour-magnitude relation does not depend on such a starburst model because most of stars are formed not by the burst but by the ordinary star formation process in 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.

The effects of merging histories of proto-objects on the angular momentum distributions of the present-time dark matter haloes are analysed. An analytical approach to the analysis of the angular momentum distributions assumes that the haloes are initially homogeneous ellipsoids and that the growth of the angular momentum of the haloes halts at their maximum expansion time. However, the maximum expansion time cannot be determined uniquely, because in the hierarchical clustering scenario each progenitor, or subunit, of the halo has its own maximum expansion time. Therefore the merging history of the halo may be important in estimating its angular momentum. Using the merger tree model by Rodrigues & Thomas, which takes into account the spatial correlations of the density fluctuations, we have investigated the effects of the merging histories on the angular momentum distributions of dark matter haloes. It was found that the merger effects, that is, the effects of the inhomogeneity of the maximum expansion times of the progenitors which finally merge together into a halo, do not strongly affect the final angular momentum distributions, so that the homogeneous ellipsoid approximation happens to be good for the estimation of the angular momentum distribution of dark matter haloes. This is because the effect of the different directions of the angular momenta of the progenitors cancels out the effect of the inhomogeneity of the maximum expansion times of the progenitors.

We have investigated time evolutions of power spectra of density fluctuations for long time after the first appearance of caustics in the expanding one-dimensional universe. It is found that when an initial power spectrum is sale-free with a power index $n$, a self-similarity of the time evolution of the power spectrum is achieved. We find that the power spectrum can be separated roughly into three regimes according to the shape of the power spectrum: the linear regime ($k < k_{nl}$ : the regime {\cal 1}),the single-caustic regime($k_{nl} < k < k_{snl}$ : the regime 2), and the multi-caustics regime($k > k_{snl}$ : the regime 3). The power index of the power spectrum in each regime has the values of $n,-1$, and $\mu$ which depends on $n$, respectively. Even in the case of an initial power-law spectrum with a cutoff scale, there might be the possibility of the self-similar evolution of the power spectrum after the appearance of the caustics. It is found, however, the self-similarity is not achieved in this case. The shape of the power spectrum on scales smaller than the cutoff scale can be separated roughly in two regimes: the virialized regime ($k_{cut}< k < k_{cs}$ : the regime 4), and the smallest-single-caustic regime ($ k > k_{cs}$ : the regime 5). The power index of the power spectrum is $\nu$ which may be determined by the distribution of singular points in the regime 4. In the regime 5, the value of the power index is -1. Moreover we show the general property about the shape of a power spectrum with a general initial condition.

We investigate stability of the scale-invariant solutions of the BBGKY equations for two-point spatial correlation functions of the density fluctuations in the strongly non-linear regime. In the case that the background skewness of the velocity field is equal to 0, we found that there is no local instability in the strongly non-linear regime. The perturbation does not grow nor does it decay. It has an only marginal stable mode. This result means that no special value of the power index of the two-point spatial correlation function are favored in terms of the stability of the solutions. In other words, the argument about the stability does not determine the power index of the two-point spatial correlation functions in the strongly non-linear regime.

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.

The effects of spatial correlations of density fluctuations on merger histories of dark matter haloes (so-called `{\it merger trees}') are analysed. We compare the mass functions of dark haloes derived by a new method for calculating merger trees, that proposed by Rodrigues \& Thomas (RT), with those given by other methods such as the Block model, the Press-Schechter formula and our own formula in which the mass functions are analytically expressed in a way that takes into consideration the spatial correlations. It is found that the mass functions given by the new method are well fit by those given by our formula. We believe that new method (RT) {\it naturally} and correctly takes into account the spatial correlations of the density fluctuations due to a calculated, grid-based realisation of the density fluctuations and so is very useful for estimating the merger tree accurately in a way that takes into consideration spatial correlations. Moreover, by applying our formula, we present an analytic expression which reproduces the mass function derived by the Block model. We therefore show clearly why and how the mass functions given by the new method and the Block model are different from each other. Furthermore, we note that the construction of merger trees is sensitive to the criterion of collapse and merging of overlapped haloes in cases in which two or more haloes happen to overlap. In fact, it is shown that the mass function is very much affected when the criterion of overlapping is changed.

We have investigated the scale-invariant solutions of the BBGKY equations for spatial correlation functions of cosmological density fluctuations and the mean relative peculiar velocity in the strongly nonlinear regime. It is found that the solutions for the mean relative physical velocity depend on the three-point spatial correlation function and the skewness of the velocity fields. We find that the stable condition in which the mean relative physical velocity vanishes on the virialized regions is satisfied only under the assumptions which Davis \& Peebles took in his paper. It is found, however, that their assumptions may not be general in real. The power index of the two-point correlation function in the strongly nonlinear regime depends on the mean relative peculiar velocity, the three-point correlation function and the skewness. If the self-similar solutions exist, the power index in the strongly nonlinear regime is related to the power index of the initial power spectrum and its relation depends on the three-point correlation function and the skewness through the mean relative peculiar velocity. Furthermore it is found that the mean relative physical velocity should have the values between 0 and the Hubble expanding velocity. When the mean relative physical velocity equals to the Hubble expanding velocity, there might exist self-similar solutions in which the power index of the two-point correlation function in the strongly nonlinear regime is independent of the initial power index $n$.

The pseudo-spectral method is proposed for following the evolution of density and velocity fluctuations at the weakly non-linear stage in the expanding universe with a good accuracy. In this method, the evolution of density and velocity fluctuations is integrated in the Fourier Space with using FFT. This method is very useful to investigate accurately the non-linear dynamics in the weakly non-linear regime. Because the pseudo-spectral method works directly in the Fourier space, it should be especially useful for examining behavior in the Fourier domain, for an example, the effects of the non-linear coupling of different wave modes on the evolution of the power spectrum. I show the results of this analysis both in one and three dimensional systems.

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 Press-Schechter(PS) formalism for the mass function of the collapsed objects are reanalyzed. The factor of two in the Press-Schechter formalism is argued to be correct in the sharp $k$-space filter even when we use the another approach proposed by Jedamzik(1994) in the cloud-in-cloud problem, which is different from the previous approach by Peacock & Heavens(1990) and Bond et al.(1991). The spatial correlation of the density fluctuations, however, had been neglected in the cloud-in-cloud problem. The effects of this spatial correlation is analyzed by using the Jedamzik formalism and it is found that this effect alter the PS mass function especially on larger mass scales. Furthermore the exact formula of deriving mass function is shown. We also find that the probability of the overlap of the collapsed objects can be neglected on very small mass scales while it might not be neglected on other mass scales.

We investigate the distance-redshift relation in the simple void model. As discussed by Moffat and Tatarski, if the observer stays at the center of the void, the observed Hubble parameter is not so different from the background Hubble parameter. However, if the position of observer is off center of the void, we must consider the peculiar velocity correction which is measured by the observed dipole anisotropy of cosmic microwave background. This peculiar velocity correction for the redshift is crucial to determine the Hubble parameter and we shall discuss this effect. Further the results of Turner et al by the N-body simulation will be also considered.

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.