Analytic expressions for the statistics of peaks of random fields with weak non-Gaussianity are provided. Specifically, the abundance and spatial correlation of peaks are represented by formulas which can be evaluated only by virtually one-dimensional integrals. We assume the non-Gaussianity is weak enough such that it is represented by linear terms of the bispectrum. The formulas are formally given in $N$-dimensional space, and explicitly given in the case of $N=1,2,3$. Some examples of peak statistics in cosmological fields are calculated for the cosmic density field and weak lensing field, assuming the weak non-Gaussianity is induced by gravity. The formulas of this paper would find a fit in many applications to statistical analyses of cosmological fields.

We investigate infrared dynamics of four-dimensional Einstein gravity in de Sitter space. We set up a general framework to investigate dynamical scaling relations in quantum/classical gravitational theories. The conformal mode dependence of Einstein gravity is renormalized to the extent that general covariance is not manifest. We point out that the introduction of an inflaton is necessary as a counterterm. We observe and postulate a duality between quantum effects in Einstein gravity and classical evolutions in an inflation (or quintessence) model. The effective action of Einstein gravity can be constructed as an inflation model with manifest general covariance. We show that $g=G_N H^2/\pi$: the only dimensionless coupling of the Hubble parameter $H^2$ and the Newton's coupling $G_N$ in Einstein gravity is screened by the infrared fluctuations of the conformal mode. We evaluate the one-loop $\beta$ function of $g$ with respect to the cosmic time $\log Ht$ as $\beta(g)=-(1/2)g^2$, i.e., $g$ is asymptotically free toward the future. The exact $\beta$ function with the backreaction of $g$ reveals the existence of the ultraviolet fixed point. It indicates that the de Sitter expansion stared at the Planck scale with a minimal entropy $S=2$. We have identified the de Sitter entropy $1/g$ with the von Neumann entropy of the conformal zero mode. The former evolves according to the screening of $g$ and the Gibbons-Hawking formula. The latter is found to increase by diffusion in the stochastic process at the horizon in a consistent way. Our Universe is located very close to the fixed point $g=0$ with a large entropy. We discuss possible physical implications of our results such as logarithmic decay of dark energy.

In the presence of the local-type primordial non-Gaussianity, it is known that the clustering of primordial black holes (PBHs) emerges even on super-horizon scales at the formation time. This effect has been investigated in the high-peak limit of the PBH formation in the radiation-dominated epoch in the literature. There is another possibility that the PBH formation takes place in the early matter-dominated epoch. In this scenario, the high-peak limit is not applicable because even initially small perturbations grow and can become a PBH. We first derive a general formula to estimate the clustering of PBHs with primordial non-Gaussianity without assuming the high-peak limit, and then apply this formula to a model of PBH formation in a matter-dominated epoch. Clustering is less significant in the case of the PBH formation in the matter-dominated epoch than that in the radiation-dominated epoch. Nevertheless, it is much larger than the Poisson shot noise in many cases. Relations to the constraints of the isocurvature perturbations by the cosmic microwave background radiation are quantitatively discussed.

In the peaks approach, the formation sites of observable structures in the Universe are identified as peaks in the matter density field. The statistical properties of the clustering of peaks are particularly important in this respect. In this paper, we investigate the large-separation expansion of the correlation function of peaks in Gaussian random fields. The analytic formula up to third order is derived, and the resultant expression can be evaluated by a combination of one-dimensional fast Fourier transforms, which are evaluated very fast. The analytic formula obtained perturbatively in the large-separation limit is compared with a method of Monte-Carlo integrations, and a complementarity between the two methods is demonstrated.

The biasing in the large-scale structure of the universe is a crucial problem in cosmological applications. The peaks model of biasing predicts a linear velocity bias of halos, which is not present in a simple model of local bias. We investigate the origin of the velocity bias in the peaks model from the viewpoint of the integrated perturbation theory, which is a nonlinear perturbation theory in the presence of general Lagrangian bias. The presence of the velocity bias in the peaks model is a consequence of the "flat constraint," ${\nabla}\delta = 0$; i.e., all the first spatial derivatives should vanish at the locations of peaks. We show that the velocity bias in the peaks model is systematically derived in the framework of the integrated perturbation theory, and then develop a formal theory to perturbatively trace the nonlinear evolution of biased objects with the flat constraint. A formula for the nonlinear velocity dispersion of peaks with the one-loop approximation is also derived.

In this paper, we explore a detectable imprint of massive fields with integer spins $s \geq 2$, which may be predicted from string theory. It was shown that such a massive non-zero spin field can generate the squeezed primordial bispectrum which depends on the angle between the two wavenumbers. We show that considering the contribution from the massive spin-2 field, the angular dependent primordial non-Gaussianity (PNG) yields a strong scale dependence in the bias parameter for the galaxy alignment, which becomes prominent at small scales. As another example of an angular dependent PNG, we also consider the primordial bispectrum where the angular dependence was introduced by a vector field, while breaking the global rotational symmetry. As a consequence, we find that the B-mode cosmic shear and non-diagonal components do not vanish. These aspects provide qualitative differences from the PNG sourced by massive non-zero spin fields.

The brightness temperature of the redshifted 21cm line brings rich information on the Inter Galactic Medium (IGM) from the Cosmic Dawn and Epoch of Reionization (EoR). While the power spectrum is a useful tool to statistically investigate the 21cm signal, the 21cm brightness temperature field is highly non-Gaussian, and the power spectrum is inadequate to characterize the non-Gaussianity. The Minkowski Functionals (MFs) are promising tools to extract non-gaussian features of the 21cm signal and give topological information such as morphology of ionized bubbles. In this work, we study the 21cm line signal in detail with MFs. To promote understanding of basic features of the 21cm signal, we calculate the MFs of not only the hydrogen neutral fraction but the matter density and spin temperature, which contribute to the brightness temperature fluctuations. We find that the structure of the brightness temperature depends mainly on the ionized fraction and the spin temperature at late and early stages of the EoR, respectively. Further, we investigate the redshift evolution of the MFs at $7 < z < 20$. Then, we find that, after the onset of reionization, the MFs reflect mainly the ionized bubble property. In addition, the MFs are sensitive to model parameters which are related to the topology of ionized bubbles and we consider the possibility of constraining the parameters by the future 21cm signal observations.

We explore the possible cosmological consequence of the gravitational coupling between the inflaton and axion-like fields. In view of the forthcoming cosmic microwave background (CMB) polarization and lensing data, we study the sensitivity of the CMB data on the cross-correlation between the curvature and axion isocurvature perturbations. Through a concrete example, we illustrate the explicit dependence of the scale dependent cross-correlation power spectrum on the axion parameters.

The impact of biasing schemes on the clustering of tracers of the large-scale structure is analytically studied in the weakly nonlinear regime. For this purpose, we use the one-loop approximation of the integrated perturbation theory together with the renormalized bias functions of various, physically motivated Lagrangian bias schemes. These include the halo, peaks and excursion set peaks model, for which we derive useful formulae for the evaluation of their renormalized bias functions. The shapes of the power spectra and correlation functions are affected by the different bias models at the level of a few percent on weakly nonlinear scales. These effects are studied quantitatively both in real and redshift space. The amplitude of the scale-dependent bias in the presence of primordial non-Gaussianity also depends on the details of the bias models. If left unaccounted for, these theoretical uncertainties could affect the robustness of the cosmological constraints extracted from galaxy clustering data.

We investigate the statistical power of higher-order statistics and cross-correlation statistics to constrain the primordial non-Gaussianity from the imaging surveys. In particular, we consider the local-type primordial non- Gaussianity and discuss how well one can tightly constrain the higher-order non-Gaussian parameters ($g_{\rm NL}$ and $\tau_{\rm NL}$) as well as the leading order parameter $f_{\rm NL}$ from the halo/galaxy clustering and weak gravitational lensing measurements. Making use of a strong scale-dependent behavior in the galaxy/halo clustering, Fisher matrix analysis reveals that the bispectra can break the degeneracy between non-Gaussian parameters ($f_{\rm NL}$, $g_{\rm NL}$ and $\tau_{\rm NL}$) and this will give simultaneous constraints on those three parameters. The combination of cross-correlation statistics further improves the constraints by factor of 2. As a result, upcoming imaging surveys like the Large Synoptic Survey Telescope have the potential to improve the constraints on the primordial non-Gaussianity much tighter than those obtained from the CMB measurement by Planck, giving us an opportunity to test the single-sourced consistency relation, $\tau_{\rm NL} \ge (36/25) f_{\rm NL}^2$.

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$ detection of RSD. Our constraint is consistent with the prediction of general relativity $f\sigma_8\sim 0.392$ within the $1-\sigma$ confidence level. When we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $f\sigma_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $f\sigma_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

Lensing of the Cosmic Microwave Background (CMB) changes the morphology of pattern of temperature fluctuations, so topological descriptors such as Minkowski Functionals can probe the gravity model responsible for the lensing. We show how the recently introduced two-to-two and three-to-one kurt-spectra (and their associated correlation functions), which depend on the power spectrum of the lensing potential, can be used to probe modified gravity theories such as $f({R})$ theories of gravity and quintessence models. We also investigate models based on effective field theory, which include the constant-$\Omega$ model, and low-energy Ho\vrava theories. Estimates of the cumulative signal-to-noise for detection of lensing-induced morphology changes, reaches ${\cal O}(10^3)$ for the future planned CMB polarization mission COrE$^{+}$. Assuming foreground removal is possible to $\ell_{max}=3000$, we show that many modified gravity theories can be rejected with a high level of significance, making this technique comparable in power to galaxy weak lensing or redshift surveys. These topological estimators are also useful in distinguishing {\em lensing} from other scattering secondaries at the level of the four-point function or trispectrum. Examples include the kinetic Sunyaev-Zel'dovich (kSZ) effect which shares, with lensing, a lack of spectral distortion. We also discuss the complication of foreground contamination from unsubtracted point sources.

In the standard perturbation theory (SPT) of self-gravitating Newtonian fluid in an expanding universe, recurrence relations for higher-order solutions are well known and play an important role both in practical applications and in theoretical investigations. The recurrence relations in Lagrangian perturbation theory (LPT), however, had not been known for a long time. Recently, two different kinds of recurrence relations in LPT have been proposed in limited cases. In this paper, we generalize those methods, and most generally derive the recurrence relations, which are capable to include any initial condition in general models of cosmology. The fastest-growing modes in the general relations are identified, and simplified recurrence relations with accurate approximation for the time-dependence are obtained.

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

General and explicit predictions from the integrated perturbation theory (iPT) for power spectra and correlation functions of biased tracers are derived and presented in the one-loop approximation. The iPT is a general framework of the nonlinear perturbation theory of cosmological density fields in presence of nonlocal bias, redshift-space distortions, and primordial non-Gaussianity. Analytic formulas of auto and cross power spectra of nonlocally biased tracers in both real and redshift spaces are derived and the results are comprehensively summarized. The main difference from previous formulas derived by the present author is to include effects of generally nonlocal Lagrangian bias and primordial non-Gaussianity, and the derivation method of the new formula is fundamentally different from the previous one. Relations to recent work on improved methods of nonlinear perturbation theory in literature are clarified and discussed.

We derive a formula for the halo/galaxy bispectrum on the basis of the integrated Perturbation Theory (iPT). In addition to the gravity-induced non-Gaussianity, we consider the non-Gaussianity of the primordial curvature perturbations, and investigate in detail the effect of such primordial non-Gaussianity on the large-scale halo/galaxy bispectrum. In iPT, the effects of primordial non-Gaussianity are wholly encapsulated in the linear (primordial) polyspectra, and we systematically calculate the contributions to the large-scale behaviors arising from the three types of primordial bispectrum (local-, equilateral-, and orthogonal-types), and primordial trispectrum of the local-type non-Gaussianity. We find that the equilateral- and orthogonal-type non-Gaussianities show distinct scale-dependent behaviors which can dominate the gravity-induced non-Gaussianity at very large scales. For the local-type non-Gaussianity, higher-order loop corrections are found to give a significantly large contribution to the halo/galaxy bispectrum of the squeezed shape, and eventually dominate over the other contributions on large scales. A diagrammatic approach based on the iPT helps us to systematically investigate an impact of such higher-order contributions to the large-scale halo/galaxy bispectrum.

In the context of upcoming large-scale structure surveys such as Euclid, it is of prime importance to quantify the effect of peculiar velocities on geometric probes. Hence the formalism to compute in redshift space the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmic fields is developed. Leveraging the partial isotropy of the target statistics, the Gram-Charlier expansion of the joint probability distribution of the field and its derivatives is reformulated in terms of the corresponding anisotropic variables. In particular, the cosmic non-linear evolution of the Minkowski functionals, together with the statistics of extrema are investigated in turn for 3D catalogues and 2D slabs. The amplitude of the non-Gaussian redshift distortion correction is estimated for these geometric probes. In 3D, gravitational perturbation theory is implemented in redshift space to predict the cosmic evolution of all relevant Gram-Charlier coefficients. Applications to the estimation of the cosmic parameters sigma(z) and beta=f/b1 from upcoming surveys is discussed. Such statistics are of interest for anisotropic fields beyond cosmology.

It is crucial to understand and model a behavior of galaxy biasing for future ambitious galaxy redshift surveys. Using 40 large cosmological N-body simulations for a standard LambdaCDM cosmology, we study the cross-correlation coefficient between matter and the halo density field, which is an indicator of the stochasticity of bias, over a wide redshift range 0\le z \le 3. The cross-correlation coefficient is important to extract information on the matter density field, e.g., by combining galaxy clustering and galaxy-galaxy lensing measurements. We compare the simulation results with integrated perturbation theory (iPT) proposed by one of the present authors and standard perturbation theory (SPT) combined with a phenomenological model of local bias. The cross-correlation coefficient derived from the iPT agrees with N-body simulation results down to r~15 (10) h^{-1}Mpc within 0.5 (1.0) % for all redshifts and halo masses we consider. The SPT with local bias does not explain complicated behaviors on quasilinear scales at low redshifts, while roughly reproduces the general behavior of the cross-correlation coefficient on fully nonlinear scales. The iPT is powerful to predict the cross-correlation coefficient down to quasilinear regimes with a high precision.

Anisotropic stress perturbations induced by primordial Gaussian vector fields create non-Gaussianity in curvature perturbations. We found that such non-Gaussianity closely resembles the local-type non-Gaussianity parametrized by $f_{\rm NL}$, and generates scale-dependent bias in large-scale structures. We also found a simple relationship between the scale-dependent bias and the power spectrum of the vector fields. When the vector fields are interpreted as primordial magnetic fields, the effective $f_{\rm NL}$ is shown to be always negative. The scale-dependent bias provides a new approach to probing primordial vector fields.

We analytically derive a more accurate formula for the power spectrum of the biased objects with the primordial non-Gaussianity parameterized not only by the non-linearity parameter fNL, but also by gNL and tauNL which characterize the trispectrum of the primordial curvature perturbations. We adopt the integrated perturbation theory which was constructed in Matsubara (2011). We discuss an inequality between fNL and tauNL in the context of the scale-dependent bias, by introducing a stochasticity parameter. We also mention higher order loop corrections into the scale-dependency of the bias parameter.

We analyze non-Gaussianity (NG) due to the primordial bispectrum and trispectrum using CMB temperature maps of WMAP 7-year data. We first apply the perturbative formulae of Minkowski functionals up to second-order NG derived by Matsubara (2010), which enable us to give limits on cubic NG parametrized with tau_NL and g_NL as well as various types of quadratic NG parametrized with f_NL. We find no signature of primordial NG in WMAP 7-year data, but give constraints on the local-type, equilateral-type, orthogonal-type f_NL: f_NL(loc)=20+-42, f_NL(eq)=-121+-208, f_NL(ort)=-129+-171, respectively, and tau_NL/10^4=-7.6+-8.7, and g_NL/10^5=-1.9+-6.4. We also find that these constraints are consistent with the limits from skewness and kurtosis parameters which characterize the perturbative corrections of MFs.

We apply the integrated perturbation theory (Matsubara 2011, PRD 83, 083518) to evaluate the scale-dependent bias in the presence of primordial non-Gaussianity. The integrated perturbation theory is a general framework of nonlinear perturbation theory, in which a broad class of bias models can be incorporated into perturbative evaluations of biased power spectrum and higher-order polyspectra. Approximations such as the high-peak limit or the peak-background split are not necessary to derive the scale-dependent bias in this framework. Applying the halo approach, previously known formulas are re-derived as limiting cases of a general formula in this work, and it is implied that modifications should be made in general situations. Effects of redshift-space distortions are straightforwardly incorporated. It is found that the slope of the scale-dependent bias on large scales is determined only by the behavior of primordial bispectrum in the squeezed limit, and is not sensitive to bias models in general. It is the amplitude of scale-dependent bias that is sensitive to the bias models. The effects of redshift-space distortions turn out to be quite small for the monopole component of the power spectrum, while the quadrupole component is proportional to the monopole component on large scales, and thus also sensitive to the primordial non-Gaussianity.

We study the lensing convergence power spectrum and its covariance for a standard LCDM cosmology. We run 400 cosmological N-body simulations and use the outputs to perform a total of 1000 independent ray-tracing simulations. We compare the simulation results with analytic model predictions. The semi-analytic model based on Smith et al.(2003) fitting formula underestimates the convergence power by ~30% at arc-minute angular scales. For the convergence power spectrum covariance, the halo model reproduces the simulation results remarkably well over a wide range of angular scales and source redshifts. The dominant contribution at small angular scales comes from the sample variance due to the number fluctuations of halos in a finite survey volume. The signal-to-noise ratio for the convergence power spectrum is degraded by the non-Gaussian covariances by up to a factor 5 for a weak lensing survey to z_s ~1. The probability distribution of the convergence power spectrum estimators, among the realizations, is well approximated by a chi-square distribution with broadened variance given by the non-Gaussian covariance, but has a larger positive tail. The skewness and kurtosis have non-negligible values especially for a shallow survey. We argue that a prior knowledge on the full distribution may be needed to obtain an unbiased estimate on the ensemble averaged band power at each angular scale from a finite volume survey.

The primordial non-Gaussianity of local type affects the clustering of dark matter halos, and the planned deep and wide photometric surveys are suitable for examining this class of non-Gaussianity. In our previous paper, we investigated the constraint from the cross correlation between CMB lensing potential and galaxy angular distribution on the primordial non-Gaussianity, without taking into account redshift slicing. To improve our previous analysis, in this paper, we add the galaxy lensing shear into our analysis and take into account redshift slicing to follow the redshift evolution of the clustering. By calculating 81 power spectra and using the Fisher matrix method, we find that the constraint on the primordial non-Gaussianity can be improved from {\Delta}fNL \sim 5.4 to 5.1 by including the galaxy-galaxy lensing shear cross correlations expected from the Hyper Suprime-Cam survey (HSC), in comparison with the constraint without any cross correlations. Moreover, the constraint can go down to {\Delta}fNL \sim 4.8 by including the galaxy-CMB lensing cross correlations from the ACTPol and Planck experiments.

Understanding a behavior of galaxy biasing is crucial for future galaxy redshift surveys. One aim is to measure the baryon acoustic oscillations (BAOs) within the precision of a few percent level. Using 30 large cosmological N-body simulations for a standard LCDM cosmology, we study the halo biasing over a wide redshift range. We compare the simulation results with theoretical predictions proposed by Matsubara (2008) which naturally incorporate the halo bias and redshift-space distortions into their formalism of perturbation theory with a resummation technique via the Lagrangian picture. The power spectrum and correlation function of halos obtained from Lagrangian resummation theory (LRT) well agree with N-body simulation results on scales of BAOs. Especially nonlinear effects on the baryon acoustic peak of the halo correlation function are accurately explained both in real and redshift space. We find that nonlinearity and scale dependence of bias are fairly well reproduced by 1-loop LRT up to k=0.35hMpc^{-1} (z=2 and 3) within a few percent level in real space and up to k=0.1hMpc^{-1} (z=2) and 0.15hMpc^{-1} (z=3) in redshift space. Thus, the LRT is very powerful for accurately extracting cosmological information in upcoming high redshift BAO surveys.

We present an improved prediction of the nonlinear perturbation theory (PT) via the Lagrangian picture, which was originally proposed by Matsubara (2008). Based on the relations between the power spectrum in standard PT and that in Lagrangian PT, we derive analytic expressions for the power spectrum in Lagrangian PT up to 2-loop order in both real and redshift spaces. Comparing the improved prediction of Lagrangian PT with $N$-body simulations in real space, we find that the 2-loop corrections can extend the valid range of wave numbers where we can predict the power spectrum within 1% accuracy by a factor of 1.0 ($z=0.5$), 1.3 (1), 1.6 (2) and 1.8 (3) vied with 1-loop Lagrangian PT results. On the other hand, in all redshift ranges, the higher-order corrections are shown to be less significant on the two-point correlation functions around the baryon acoustic peak, because the 1-loop Lagrangian PT is already accurate enough to explain the nonlinearity on those scales in $N$-body simulations.

Using 1000 ray-tracing simulations for a {\Lambda}-dominated cold dark model in Sato et al. (2009), we study the covariance matrix of cosmic shear correlation functions, which is the standard statistics used in the previous measurements. The shear correlation function of a particular separation angle is affected by Fourier modes over a wide range of multipoles, even beyond a survey area, which complicates the analysis of the covariance matrix. To overcome such obstacles we first construct Gaussian shear simulations from the 1000 realizations, and then use the Gaussian simulations to disentangle the Gaussian covariance contribution to the covariance matrix we measured from the original simulations. We found that an analytical formula of Gaussian covariance overestimates the covariance amplitudes due to an effect of finite survey area. Furthermore, the clean separation of the Gaussian covariance allows to examine the non-Gaussian covariance contributions as a function of separation angles and source redshifts. For upcoming surveys with typical source redshifts of z_s=0.6 and 1.0, the non-Gaussian contribution to the diagonal covariance components at 1 arcminute scales is greater than the Gaussian contribution by a factor of 20 and 10, respectively. Predictions based on the halo model qualitatively well reproduce the simulation results, however show a sizable disagreement in the covariance amplitudes. By combining these simulation results we develop a fitting formula to the covariance matrix for a survey with arbitrary area coverage, taking into account effects of the finiteness of survey area on the Gaussian covariance.

The standard nonlinear perturbation theory of the gravitational instability is extended to incorporate the nonlocal bias, redshift-space distortions, and primordial non-Gaussianity. We show that local Eulerian bias is not generally compatible to local Lagrangian bias in nonlinear regime. The Eulerian and Lagrangian biases are nonlocally related order by order in the general perturbation theory. The relation between Eulerian and Lagrangian kernels of density perturbations with biasing are derived. The effects of primordial non-Gaussianity and redshift-space distortions are also incorporated in our general formalism, and diagrammatic methods are introduced. Vertex resummations of higher-order perturbations in the presence of bias are considered. Resummations of Lagrangian bias are shown to be essential to handle biasing schemes in a general framework.

We study the sample variance of the matter power spectrum for the standard Lambda Cold Dark Matter universe. We use a total of 5000 cosmological N-body simulations to study in detail the distribution of best-fit cosmological parameters and the baryon acoustic peak positions. The obtained distribution is compared with the results from the Fisher matrix analysis with and without including non-Gaussian errors. For the Fisher matrix analysis, we compute the derivatives of the matter power spectrum with respect to cosmological parameters using directly full nonlinear simulations. We show that the non-Gaussian errors increase the unmarginalized errors by up to a factor 5 for k_{max}=0.4h/Mpc if there is only one free parameter provided other parameters are well determined by external information. On the other hand, for multi-parameter fitting, the impact of the non-Gaussian errors is significantly mitigated due to severe parameter degeneracies in the power spectrum. The distribution of the acoustic peak positions is well described by a Gaussian distribution, with its width being consistent with the statistical interval predicted from the Fisher matrix. We also examine systematic bias in the best-fit parameter due to the non-Gaussian errors. The bias is found to be smaller than the 1 sigma statistical error for both the cosmological parameters and the acoustic scale positions.

Recent studies have shown that the primordial non-Gaussianity affects clustering of dark matter halos through a scale-dependent bias and various constraints on the non-Gaussianity through this scale-dependent bias have been placed. Here we introduce the cross-correlation between the CMB lensing potential and the galaxy angular distribution to effectively extract information about the bias from the galaxy distribution. Then, we estimate the error of non-linear parameter, f_NL, for the on-going CMB experiments and galaxy surveys, such as Planck and Hyper Suprime-Cam (HSC). We found that for the constraint on f_NL with Planck and HSC, the wide field galaxy survey is preferable to the deep one, and the expected error on f_NL can be as small as: \Delta f_NL ~ 20 for b_0 = 2 and \Delta f_NL ~ 10 for b_0 = 4, where b_0 is the linear bias parameter. It is also found that future wide field galaxy survey could achieve \Delta f_NL ~ 5 with CMB prior from Planck if one could observe highly biased objects at higher redshift (z ~ 2).

Analytic formulas of Minkowski functionals in two-dimensional random fields are derived, including effects of second-order non-Gaussianity in the presence of both the bispectrum and trispectrum. The set of formulas provides a promising method to constrain the primordial non-Gaussianity of the universe by temperature fluctuations in the cosmic microwave background radiation. In a case of local-type non-Gaussianity, the Minkowski functionals are analytically given by powers of quadratic and cubic parameters, $f_{\rm NL}$ and $g_{\rm NL}$. Our formulas are not restricted to this particular model, and applicable to a wide class of non-Gaussian models. The analytic formulas are compared to numerical evaluations from non-Gaussian realizations of temperature maps, showing very good agreements.

We study the effect of primordial isocurvature perturbations on non-Gaussian properties of CMB temperature anisotropies. We consider generic forms of the non-linearity of isocurvature perturbations which can be applied to a wide range of theoretical models. We derive analytical expressions for the bispectrum and the Minkowski Functionals for CMB temperature fluctuations to describe the non-Gaussianity from isocurvature perturbations. We find that the isocurvature non-Gaussianity in the quadratic isocurvature model, where the isocurvature perturbation S is written as a quadratic function of the Gaussian variable sigma, S=sigma^2-<sigma^2>, can give the same signal-to-noise as f_NL=30 even if we impose the current observational limit on the fraction of isocurvature perturbations contained in the primordial power spectrum alpha. We give constraints on isocurvature non-Gaussianity from Minkowski Functionals using WMAP 5-year data. We do not find a significant signal of the isocurvature non-Gaussianity. For the quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature fraction alpha<0.070 (95% CL) for a scale invariant spectrum which is comparable to the limit obtained from the power spectrum.

We use 5000 cosmological N-body simulations of 1(Gpc/h)^3 box for the concordance LCDM model in order to study the sampling variances of nonlinear matter power spectrum. We show that the non-Gaussian errors can be important even on large length scales relevant for baryon acoustic oscillations (BAO). Our findings are (1) the non-Gaussian errors degrade the cumulative signal-to-noise ratios (S/N) for the power spectrum amplitude by up to a factor of 2 and 4 for redshifts z=1 and 0, respectively. (2) There is little information on the power spectrum amplitudes in the quasi-nonlinear regime, confirming the previous results. (3) The distribution of power spectrum estimators at BAO scales, among the realizations, is well approximated by a Gaussian distribution with variance that is given by the diagonal covariance component. (4) For the redshift-space power spectrum, the degradation in S/N by non-Gaussian errors is mitigated due to nonlinear redshift distortions. (5) For an actual galaxy survey, the additional shot noise contamination compromises the cosmological information inherent in the galaxy power spectrum, but also mitigates the impact of non-Gaussian errors. The S/N is degraded by up to 30% for a WFMOS-type survey. (6) The finite survey volume causes additional non-Gaussian errors via the correlations of long-wavelength fluctuations with the fluctuations we want to measure, further degrading the S/N values by about 30% even at high redshift z=3.

Employing the perturbative treatment of gravitational clustering, we discuss possible effects of primordial non-Gaussianity on the matter power spectrum. As gravitational clustering develops, the coupling between different Fourier modes of density perturbations becomes important and the primordial non-Gaussianity which intrinsically possesses a non-trivial mode-correlation can affect the late-time evolution of the power spectrum. We quantitatively estimate the non-Gaussian effect on power spectrum from the perturbation theory. The potential impact on the cosmological parameter estimation using the power spectrum are investigated based on the Fisher-matrix formalism. In addition, on the basis of the local biasing prescription, non-Gaussian effects on the galaxy power spectrum are considered, showing that the scale-dependent biasing arises from a local-type primordial non-Gaussianity. On the other hand, an equilateral-type non-Gaussianity does not induce such scale-dependence because of weaker mode-correlations between small and large Fourier modes.

We use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a \Lambda Cold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum.

The nonlinear perturbation theory of gravitational instability is extended to include effects of both biasing and redshift-space distortions, which are inevitable in predicting observable quantities in galaxy surveys. The precise determination of scales of baryon acoustic oscillations is crucial to investigate the nature of dark energy by galaxy surveys. We find that a local Lagrangian bias and redshift-space distortions are naturally incorporated in our formalism of perturbation theory with a resummation technique via the Lagrangian picture. Our formalism is applicable to any biasing scheme which is local in Lagrangian space, including the halo bias as a special case. Weakly nonlinear effects on halo clustering in redshift space are analytically given.We assume only a fundamental idea of the halo model: haloes form according to the extended Press-Schechter theory, and the spatial distributions are locally biased in Lagrangian space. There is no need for assuming the spherical collapse model to follow the dynamical evolution, which is additionally assumed in standard halo prescriptions. One-loop corrections to the power spectrum and correlation function of haloes in redshift space are explicitly derived and presented. Instead of relying on expensive numerical simulations, our approach provides an analytic way of investigating the weakly nonlinear effects, simultaneously including the nonlinear biasing and nonlinear redshift-space distortions. Nonlinearity introduces a weak scale dependence in the halo bias. The scale dependence is a smooth function in Fourier space, and the bias does not critically change the feature of baryon acoustic oscillations in the power spectrum. The same feature in the correlation function is less affected by nonlinear effects of biasing.

We present an analysis of the Minkowski Functionals (MFs) describing the WMAP three-year temperature maps to place limits on possible levels of primordial non-Gaussianity. In particular, we apply perturbative formulae for the MFs to give constraints on the usual non-linear coupling constant fNL. The theoretical predictions are found to agree with the MFs of simulated CMB maps including the full effects of radiative transfer. The agreement is also very good even when the simulation maps include various observational artifacts, including the pixel window function, beam smearing, inhomogeneous noise and the survey mask. We find accordingly that these analytical formulae can be applied directly to observational measurements of fNL without relying on non-Gaussian simulations. Considering the bin-to-bin covariance of the MFs in WMAP in a chi-square analysis, we find that the primordial non-Gaussianity parameter is constrained to lie in the range -70<fNL<91 at 95% C.L. using the Q+V+W co-added maps.

We critically examine how well the evolution of large-scale density perturbations is followed in cosmological $N$-body simulations. We first run a large volume simulation and perform a mode-by-mode analysis in three-dimensional Fourier space. We show that the growth of large-scale fluctuations significantly deviates from linear theory predictions. The deviations are caused by {\it nonlinear} coupling with a small number of modes at largest scales owing to finiteness of the simulation volume. We then develop an analytic model based on second-order perturbation theory to quantify the effect. Our model accurately reproduces the simulation results. For a single realization, the second-order effect appears typically as ``zig-zag'' patterns around the linear-theory prediction, which imprints artificial ``oscillations'' that lie on the real baryon-acoustic oscillations. Although an ensemble average of a number of realizations approaches the linear theory prediction, the dispersions of the realizations remain large even for a large simulation volume of several hundred megaparsecs on a side. For the standard $\Lambda$CDM model, the deviations from linear growth rate are as large as 10 percent for a simulation volume with $L = 500h^{-1}$Mpc and for a bin width in wavenumber of $\Delta k = 0.005h$Mpc$^{-1}$, which are comparable to the intrinsic variance of Gaussian random realizations. We find that the dispersions scales as $\propto L^{-3/2} \Delta k^{-1/2}$ and that the mean dispersion amplitude can be made smaller than a percent only if we use a very large volume of $L > 2h^{-1}$Gpc. The finite box size effect needs to be appropriately taken into account when interpreting results from large-scale structure simulations for future dark energy surveys using baryon acoustic oscillations.

We study the large-scale anisotropic two-point correlation function using 46,760 luminous red galaxies at redshifts 0.16 -- 0.47 from the Sloan Digital Sky Survey. We measure the correlation function as a function of separations parallel and perpendicular to the line-of-sight in order to take account of anisotropy of the large-scale structure in redshift space. We find a slight signal of baryonic features in the anisotropic correlation function, i.e., a ``baryon ridge'' which corresponds to a baryon acoustic peak in the spherically averaged correlation function which has already been reported using the same sample. The baryon ridge has primarily a spherical structure with a known radius in comoving coordinates. It enables us to divide the redshift distortion effects into dynamical and geometrical components and provides further constraints on cosmological parameters, including the dark energy equation-of-state. With an assumption of a flat $\Lambda$ cosmology, we find the best-fit values of $\Omega_{\rm m} = 0.218^{+0.047}_{-0.037}$ and $\Omega_{\rm b} = 0.047^{+0.016}_{-0.016}$ (68% C.L.) when we use the overall shape of the anisotropic correlation function of $40<s<200\himpc$ including a scale of baryon acoustic oscillations. When an additional assumption $\Omega_{\rm b}h^2=0.024$ is adopted, we obtain $\Omega_{\rm DE}=0.770^{+0.051}_{-0.040}$ and $w=-0.93^{+0.45}_{-0.35}$. These constraints are estimated only from our data of the anisotropic correlation function, and they agree quite well with values both from the cosmic microwave background (CMB) anisotropies and from other complementary statistics using the LRG sample. With the CMB prior from the 3 year WMAP results, we give stronger constraints on those parameters.

We develop a new approach to study the nonlinear evolution in the large-scale structure of the Universe both in real space and in redshift space, extending the standard perturbation theory of gravitational instability. Infinite series of terms in standard Eulerian perturbation theory are resummed as a result of our starting from a Lagrangian description of perturbations. Delicate nonlinear effects on scales of the baryon acoustic oscillations are more accurately described by our method than the standard one. Our approach differs from other resummation techniques recently proposed, such as the renormalized perturbation theory, etc., in that we use simple techniques and thus resulting equations are undemanding to evaluate, and in that our approach is capable of quantifying the nonlinear effects in redshift space. The power spectrum and correlation function of our approach are in good agreement with numerical simulations in literature on scales of baryon acoustic oscillations. Especially, nonlinear effects on the baryon acoustic peak of the correlation function are accurately described both in real space and in redshift space. Our approach provides a unique opportunity to analytically investigate the nonlinear effects on baryon acoustic scales in observable redshift space, which is requisite in constraining the nature of dark energy, the curvature of the Universe, etc., by redshift surveys.

A method of counts-in-cells analysis of galaxy distribution is investigated with arbitrary smoothing functions in obtaining the galaxy counts. We explore the possiblity of optimizing the smoothing function, considering a series of $m$-weight Epanechnikov kernels. The popular top-hat and Gaussian smoothing functions are two special cases in this series. In this paper, we mainly consider the second moments of counts-in-cells as a first step. We analytically derive the covariance matrix among different smoothing scales of cells, taking into account possible overlaps between cells. We find that the Epanechnikov kernel of $m=1$ is better than top-hat and Gaussian smoothing functions in estimating cosmological parameters. As an example, we estimate expected parameter bounds which comes only from the analysis of second moments of galaxy distributions in a survey which is similar to the Sloan Digital Sky Survey.

We derive analytical formulae for the Minkowski Functions of the cosmic microwave background (CMB) and large-scale structure (LSS) from primordial non-Gaussianity. These formulae enable us to estimate a non-linear coupling parameter, f_NL, directly from the CMB and LSS data without relying on numerical simulations of non-Gaussian primordial fluctuations. One can use these formulae to estimate statistical errors on f_NL from Gaussian realizations, which are much faster to generate than non-Gaussian ones, fully taking into account the cosmic/sampling variance, beam smearing, survey mask, etc. We show that the CMB data from the Wilkinson Microwave Anisotropy Probe should be sensitive to |f_NL|\simeq 40 at the 68% confidence level. The Planck data should be sensitive to |f_NL|\simeq 20. As for the LSS data, the late-time non-Gaussianity arising from gravitational instability and galaxy biasing makes it more challenging to detect primordial non-Gaussianity at low redshifts. The late-time effects obscure the primordial signals at small spatial scales. High-redshift galaxy surveys at z>2 covering \sim 10Gpc^3 volume would be required for the LSS data to detect |f_NL|\simeq 100. Minkowski Functionals are nicely complementary to the bispectrum because the Minkowski Functionals are defined in real space and the bispectrum is defined in Fourier space. This property makes the Minksowski Functionals a useful tool in the presence of real-world issues such as anisotropic noise, foreground and survey masks. Our formalism can be extended to scale-dependent f_NL easily.

Fourier methods are fundamental tools to analyze random fields. Statistical structures of homogeneous Gaussian random fields are completely characterized by the power spectrum. In non-Gaussian random fields, polyspectra, higher-order counterparts of the power spectrum, are usually considered to characterize statistical information contained in the fields. However, it is not trivial how the Fourier modes are distributed in general non-Gaussian fields. In this paper, distribution functions of Fourier modes are directly considered and their explicit relations to the polyspectra are given. Under the condition that any of the polyspectra does not diverge, the distribution function is expanded by dimensionless combinations of polyspectra and a total volume in which the Fourier transforms are performed. The expression up to second order is generally given, and higher-order results are also derived in some cases. A concept of N-point distribution function of Fourier modes are introduced and explicitly calculated. Among them, the one-point distribution function is completely given in a closed form up to arbitrary order. As an application, statistics of Fourier phases are explored in detail. A number of aspects regarding statistical properties of phases are found. It is clarified, for the first time, how phase correlations arise in general non-Gaussian fields. Some of our analytic results are tested against numerical realizations of non-Gaussian fields, showing good agreements.

We present a first analysis of the clustering of SDSS galaxies using the distribution function of the sum of Fourier phases. This statistic was recently proposed by one of authors as a new method to probe phase correlations of cosmological density fields. Since the Fourier phases are statistically independent of Fourier amplitudes, the phase statistic plays a complementary role to the conventional two-point statistics of galaxy clustering. In particular, we focus on the distribution functions of phase sum over three closed wavevectors as a function of the configurations of triangle wavevectors. We find that the observed distribution functions of phase sum are in good agreement with the lowest-order approximation from perturbation theory. For direct comparison with observations, we construct mock catalogs from N-body simulations taking account of several observational effects such as the survey geometry, the redshift distortion, and the discreteness due to the limited number of data. Indeed the observed phase correlations for the galaxies in the range of absolute magnitude -22<Mr<-18 agree well with those for Lambda CDM predictions with sigma8=0.9 evolved from Gaussian initial condition. This agreement implies that the galaxy biasing is approximately linear in redshift space. Instead, assuming that the galaxy biasing is described by a quadratic deterministic function at scale larger than 30Mpc/h, we can constrain the ratio of the quadratic biasing parameter b2 to the linear biasing parameter b1. We find that resulting b2/b1 is mainly dependent on sigma8 and that it is well fitted by b2/b1=0.54(+-0.06)-0.62(+-0.08)sigma8 in CDM models. Indeed, b2/b1 is nearly zero when sigma8=0.9.

We measure the topology of volume-limited galaxy samples selected from a parent sample of 314,050 galaxies in the Sloan Digital Sky Survey(SDSS), which is now complete enough to describe the fully three-dimensional topology and its dependence on galaxy properties. We compare the observed genus statistic G\nu_f) to predictions for a Gaussian random field and to the genus measured for mock surveys constructed from new large-volume simulations of the LambdaCDM cosmology. In this analysis we carefully examine the dependence of the observed genus statistic on the Gaussian smoothing scale R_G from 3.5 to 11 h^-1Mpc and on the luminosity of galaxies over the range -22.50<M_r<-18.5. The void multiplicity A_V is less than unity at all smoothing scales. Because A_V cannot become less than 1 through gravitational evolution, this result provides strong evidence for biased galaxy formation in low density environments. We also find clear evidence of luminosity bias of topology within the volume-limited sub-samples. The shift parameter Delta nu indicates that the genus of brighter galaxies shows a negative shift toward a ``meatball'' (i.e. cluster-dominated) topology, while faint galaxies show a positive shift toward a ``bubble'' (i.e. void-dominated) topology. The transition from negative to positive shift occurs approximately at the characteristic absolute magnitude M_{r*}=-20.4. Even in this analysis of the largest galaxy sample to date, we detect the influence of individual large-scale structures, as the shift parameter Delta nu and cluster multiplicity A_C reflect (at ~3 sigma) the presence of the Sloan Great Wall and a X-shaped structure which runs for several hundred Mpc across the survey volume.

We explore a signature of phase correlations in Fourier modes of dark matter density fields induced by nonlinear gravitational clustering. We compute the distribution function of the phase sum of the Fourier modes, theta_k1+theta_k2+theta_k3, for triangle wavevectors satisfying k1+k2+k3=0, and compare with the analytic prediction in perturbation theory recently derived by one of us. Using a series of cosmological N-body simulations, we extensively examine the time evolution and the dependence on the configuration of triangles and the sampling volume. Overall we find that the numerical results are remarkably consistent with the analytic formula from the perturbation theory. Interestingly the validity of the perturbation theory at a scale k is determined by P(k)/V_samp, the ratio of the power spectrum P(k) and the sampling volume V_samp, not by k^3P(k) as in the case of the conventional cosmological perturbation theory. Consequently this statistics of phase correlations is sensitive to the size of the sampling volume itself. This feature does not show up in more conventional cosmological statistics including the one-point density distribution function and the two-point correlation functions except as a sample-to-sample variation. Similarly if the sampling volume size V_samp is fixed, the stronger phase correlation emerges first at the wavevector where P(k) becomes largest, i.e., in linear regimes according to the standard cosmological perturbation theory, while the distribution of the phase sum stays fairly uniform in nonlinear regimes. The above feature can be naturally understood from the corresponding density structures in real space as we discuss in detail.

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at $8h^{-1}{\rm Mpc} < s < 500h^{-1}{\rm Mpc}$, with particular attention to its baryonic signature. Our sample consists of 19986 quasars extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is $0.72 \le z \le 2.24$ (the mean redshift is $\bar z = 1.46$) and the reddening-corrected $i$-band apparent magnitude range is $15.0 \le m_{i,{\rm rc}} \le 19.1$. Due to the relatively low number density of the quasar sample, the bump in the power spectrum due to the baryon density, $\Omega_{\rm b}$, is not clearly visible. The effect of the baryon density is, however, to distort the overall shape of the 2PCF.The degree of distortion makes it an interesting alternate measure of the baryonic signature. Assuming a scale-independent linear bias and the spatially flat universe, i.e., $\Omega_{\rm b} + \Omega_{\rm d} + \Omega_\Lambda =1$, where $\Omega_{\rm d}$ and $\Omega_\Lambda$ denote the density parameters of dark matter and the cosmological constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on $\Omega_{\rm b}$ and $\Omega_\Lambda$. Our result is fitted as $0.80- 2.8\Omega_{\rm b} < \Omega_\Lambda < 0.90 - 1.4\Omega_{\rm b}$ at the 2$\sigma$ confidence level, which is consistent with results from other cosmological observations such as WMAP. (abridged)

We analyze the Minkowski functionals with a large $N$-body simulation of a standard $\Lambda$CDM model, focusing on transition scales between linear and non-linear gravitational evolution. We numerically calculate the Minkowski functionals with sufficient accuracies to investigate the transition scales, 10--$50\himpc$. The results are compared with analytic formulae of linear and second-order perturbation theories. We first show that the skewness parameters of the density fluctuations, which are important in second-order analytic formulae of the Minkowski functionals, are in good agreement with the perturbation theory. Considering relative differences between the Minkowski functionals of the analytic formulae and that of the simulation data, we evaluate accuracy levels of the predictions of the perturbation theory. When the straightforward threshold $\nu$ by density value is used in Minkowski functionals, the accuracy of the second-order perturbation theory is within 10% for smoothing length $R > 15\himpc$, and within a several % for $R > 20\himpc$. The accuracies of the linear theory are 2--5 times worse than that. When the rescaled threshold by volume fraction, $\nu_{\rm f}$ is used, accuracies of both linear and second-order theories are within a few % on all scales of $10\himpc < R < 50\himpc$.

Recent development of galaxy surveys enables us to investigate the deep universe of high redshift. We quantitatively present the physical information extractable from the observable correlation function in deep redshift space in a framework of the linear theory. The correlation function depends on the underlying power spectrum, velocity distortions, and the Alcock-Paczy\'nski (AP) effect. The underlying power spectrum is sensitive to the constituents of matters in the universe, the velocity distortions are sensitive to the galaxy bias as well as the amount of total matter, and the Alcock-Paczy\'nski effect is sensitive to the dark energy components. Measuring the dark energy by means of the baryonic feature in the correlation function is one of the most interesting applications. We show that the ``baryon ridge'' in the correlation function serves as a statistically circular object in the AP effect. In order to sufficiently constrain the dark energy components, the redshift range of the galaxy survey should be as broad as possible. The survey area on the sky should be smaller at deep redshifts than at shallow redshifts to keep the number density as dense as possible. We illustrate an optimal survey design that are useful in cosmology. Assuming future redshift surveys of $z\simlt 3$ which are within reach of the present-day technology, achievable error bounds on cosmological parameters are estimated by calculating the Fisher matrix. According to an illustrated design, the equation of state of dark energy can be constrained within $\pm 5%$ error assuming that the bias is unknown and marginalized over. Even when all the other cosmological parameters should be simultaneously determined, the error bound for the equation of state is up to $\pm 10%$.

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

In modeling strong gravitational lens systems, one often adopts simple models, such as singular isothermal elliptical density plus lowest-order external perturbation. However, such simple models may mislead us if the real mass distribution is more complicated than that in the assumed models. In particular, assumptions on mass models are crucial in studying flux ratio anomalies that have been suggested as evidence for a cold dark-matter substructure. We reinvestigated four quadruple lens systems using power-law Fourier models, which have advantages of clear physical meanings and the applicability of a linear method, as well as a simple singular isothermal elliptical density model. We also investigated the effect of external perturbations, including a singular isothermal sphere, lowest-order expansion, and next-order expansion. We have found that the $\cos3\theta$ terms of the primary galaxy and/or of external perturbation significantly reduce the \chi^2 in PG 1115+080 and B1422+231. In particular, we could reproduce the flux ratios of B1422+231 with next-order external perturbation assuming 5% flux uncertainties, suggesting that external perturbation cannot be described by a simple singular isothermal sphere approximation. On the other hand, we could not fit B0712+472 and B2045+265 very well even with our models, although the \chi^2 were reduced compared with the case of using simple models. Our results clearly demonstrate that both the primary lens galaxy and the external perturbation are often more complicated than we usually assume.

We present estimates of cosmological parameters from the application of the Karhunen-Loeve transform to the analysis of the 3D power spectrum of density fluctuations using Sloan Digital Sky Survey galaxy redshifts. We use Omega_m*h and f_b = Omega_b/Omega_m to describe the shape of the power spectrum, sigma8 for the (linearly extrapolated) normalization, and beta to parametrize linear theory redshift space distortions. On scales k < 0.16 h/Mpc, our maximum likelihood values are Omega_m*h = 0.264 +/-0.043, f_b = 0.286 +/- 0.065, sigma8 = 0.966 +/- 0.048, and beta = 0.45 +/- 0.12. When we take a prior on Omega_b from WMAP, we find Omega_m*h = 0.207 +/- 0.030, which is in excellent agreement with WMAP and 2dF. This indicates that we have reasonably measured the gross shape of the power spectrum but we have difficulty breaking the degeneracy between Omega_m*h and f_b because the baryon oscillations are not resolved in the current spectroscopic survey window function.

Eigenmode analysis is one of the most promising methods of analyzing large data sets in ongoing and near-future galaxy surveys. In such analyses, a fast evaluation of the correlation matrix in arbitrary cosmological models is crucial. The observational effects, including peculiar velocity distortions in redshift space, light-cone effects, selection effects, and effects of the complex shape of the survey geometry, should be taken into account in the analysis. In the framework of the linear theory of gravitational instability, we provide the methodology to quickly compute the correlation matrix. Our methods are not restricted to shallow redshift surveys, arbitrarily deep samples can be dealt with as well. Therefore, our methods are useful in constraining the geometry of the universe and the dark energy component, as well as the power spectrum of galaxies, since ongoing and near-future galaxy surveys probe the universe at intermediate to deep redshifts, z ~ 0.2--5. In addition to the detailed methods to compute the correlation matrix in 3-dimensional redshift surveys, methods to calculate the matrix in 2-dimensional projected samples are also provided. Prospects of applying our methods to likelihood estimation of the cosmological parameters are discussed.

We present the general relationship between phase correlations and the hierarchy of polyspectra in the Fourier space, and the new theoretical understanding of the phase information is provided. Phase correlations are related to the polyspectra only through the non-uniform distributions of the phase sum $\theta_{\sbfm{k}_1} + ... + \theta_{\sbfm{k}_N}$ with closed wave vectors, $\bfm{k}_1 + ... + \bfm{k}_N = 0$. The exact relationship is given by the infinite series, which one can truncate in a consistent manner. The method to calculate the series to arbitrary order is explained, and the explicit expression of the first-order approximation is given. A numerical demonstration proves that the distribution of the phase sum is a robust estimator and provides an alternative statistic to search for the non-Gaussianity.

We show the apparent redshift-space clustering of galaxies in redshift range of 0.2--0.4 provides surprisingly useful constraints on dark energy component in the universe, because of the right balance between the density of objects and the survey depth. We apply Fisher matrix analysis to the the Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS), as a concrete example. Possible degeneracies in the evolution of the equation of state (EOS) and the other cosmological parameters are clarified.

We formulate a general method for perturbative evaluations of statistics of smoothed cosmic fields, and provide useful formulas in application of the perturbation theory to various statistics. This formalism is an extensive generalization of the method used by Matsubara (1994) who derived a weakly nonlinear formula of the genus statistic in a 3D density field. After describing the general method, we apply the formalism to a series of statistics, including genus statistics, level-crossing statistics, Minkowski functionals, and a density extrema statistic, regardless of the dimensions in which each statistic is defined. The relation between the Minkowski functionals and other geometrical statistics is clarified. These statistics can be applied to several cosmic fields, including 3D density field, 3D velocity field, 2D projected density field, and so forth. The results are detailed for second order theory of the formalism. The effect of the bias is discussed. The statistics of smoothed cosmic fields as functions of rescaled threshold by volume-fraction are discussed in the framework of second-order perturbation theory. In CDM-like models, their functional deviations from linear predictions plotted against the rescaled threshold are generally much smaller than that plotted against the direct threshold. There is still slight meat-ball shift against rescaled threshold, which is characterized by asymmetry in depths of troughs in the genus curve. A theory-motivated asymmetry factor in genus curve is proposed.

We present the first analysis of three-dimensional genus statistics for the SDSS EDR galaxy sample. Due to the complicated survey volume and the selection function, analytic predictions of the genus statistics for this sample are not feasible, therefore we construct extensive mock catalogs from N-body simulations in order to compare the observed data with model predictions. This comparison allows us to evaluate the effects of a variety of observational systematics on the estimated genus for the SDSS sample, including the shape of the survey volume, the redshift distortion effect, and the radial selection function due to the magnitude limit. The observed genus for the SDSS EDR galaxy sample is consistent with that predicted by simulations of a $\Lambda$-dominated spatially-flat cold dark matter model. Standard ($\Omega_0=1$) cold dark matter model predictions do not match the observations. We discuss how future SDSS galaxy samples will yield improved estimates of the genus.

We estimate how clustering in large-scale redshift surveys can constrain various cosmological parameters. Depth and sky coverage of modern redshift surveys are greater than ever, opening new possibilities for statistical analysis. We have constructed a novel maximum likelihood technique applicable to deep redshift surveys of wide sky coverage by taking into account the effects of both curvature and linear velocity distortions. The Fisher information matrix is evaluated numerically to show the bounds derived from a given redshift sample. We find that intermediate-redshift galaxies, such as the Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey, can constrain cosmological parameters, including the cosmological constant, unexpectedly well. The importance of the dense as well as deep sampling in designing redshift surveys is emphasized.

We present measurements of parameters of the 3-dimensional power spectrum of galaxy clustering from 222 square degrees of early imaging data in the Sloan Digital Sky Survey. The projected galaxy distribution on the sky is expanded over a set of Karhunen-Loeve eigenfunctions, which optimize the signal-to-noise ratio in our analysis. A maximum likelihood analysis is used to estimate parameters that set the shape and amplitude of the 3-dimensional power spectrum. Our best estimates are Gamma=0.188 +/- 0.04 and sigma_8L = 0.915 +/- 0.06 (statistical errors only), for a flat Universe with a cosmological constant. We demonstrate that our measurements contain signal from scales at or beyond the peak of the 3D power spectrum. We discuss how the results scale with systematic uncertainties, like the radial selection function. We find that the central values satisfy the analytically estimated scaling relation. We have also explored the effects of evolutionary corrections, various truncations of the KL basis, seeing, sample size and limiting magnitude. We find that the impact of most of these uncertainties stay within the 2-sigma uncertainties of our fiducial result.

We show how the cosmological constant can be estimated from redshift surveys at different redshifts, using maximum-likelihood techniques. The apparent redshift-space clustering on large scales (\simgt 20 \himpc) are affected in the radial direction by infall, and curvature influences the apparent correlations in the transverse direction. The relative strengths of the two effects will strongly vary with redshift. Using a simple idealized survey geometry, we compute the smoothed correlation matrix of the redshift-space correlation function, and the Fisher matrix for $\omm$ and $\oml$. These represent the best possible measurement of these parameters given the geometry. We find that the likelihood contours are turning, according to the behavior of the angular-diameter distance relation. The clustering measures from redshift surveys at intermediate-to-high redshifts can provide a surprisingly tight constraint on $\oml$. We also estimate confidence contours for real survey geometries, using the SDSS LRG and QSO surveys as specific examples. We believe that this method will become a practical tool to constrain the nature of the dark energy.

The topology of weak lensing fields is studied using the 2-dimensional genus statistic. Simulated fields of the weak lensing convergence are used to focus on the effect of nonlinear gravitional evolution and to model the statistical errors expected in observational surveys. For large smoothing angles, the topology is in agreement with the predictions from linear theory. On smoothing angles smaller than 10', the genus curve shows the non-Gaussian signatures of gravitational clustering and differs for open and flat cold dark matter models. Forthcoming surveys with areas larger than 10 square degrees should have adequate signal-to-noise to measure the non-Gaussian shape and the $\Omega$-dependence of the genus statistic.

The gravitational lensing, as well as the velocity field and the cosmological light-cone warp, changes the observed correlation function of high-redshift objects. We present an analytical expression of 3D correlation function, simultaneously including those three effects. When two objects are separated over several hundreds Mpc along a line of sight, the observed correlation function is dominated by the effect of gravitational lensing rather than the intrinsic clustering. For a canonical lambda-CDM model, the lensing signals in the galaxy-galaxy and galaxy-QSO correlations are beyond noise levels in large-scale redshift surveys like the Sloan Digital Sky Survey.

We present the first results of the Karhunen-Lo\`eve (KL) eigenmodes applied to real data of the Las Campanas Redshift Survey (LCRS) to simultaneously measure the values of the redshift-distortion parameter, $\beta = \Omega_0^{0.6}/b$, the linearly extrapolated normalization, $\sigma_8^{\rm L}$, and the CDM shape parameter, $\Gamma = \Omega_0 h$. The results of our numerical likelihood analysis indicate a low value of $\beta = 0.30 \pm 0.39$, a shape parameter $\Gamma = 0.16 \pm 0.10$, and a linearly extrapolated normalization $\sigma_8^{\rm L} = 0.79 \pm 0.08$, which are consistent with a low density universe, $\Omega_0 \simlt 0.5$.

A general formula for the correlation function in redshift space is derived in linear theory. The formula simultaneously includes wide-angle effects and cosmological distortions. The formula is applicable to any pair with arbitrary angle $\theta$ between lines of sight, and arbitrary redshifts, $z_1$, $z_2$, which are not necessarily small. The effects of the spatial curvature both on geometry and on fluctuation spectrum are properly taken into account, and thus our formula holds in a Friedman-Lema\^{\i}tre universe with arbitrary cosmological parameters $\Omega_0$ and $\lambda_0$. We illustrate the pattern of the resulting correlation function with several models, and also show that validity region of the conventional distant observer approximation is $\theta \le 10^\circ$.

If one wants to represent the galaxy number density at some point in terms of only the mass density at the same point, there appears the stochasticity in such a relation, which is referred to as ``stochastic bias''. The stochasticity is there because the galaxy number density is not merely a local function of a mass density field, but it is a nonlocal functional, instead. Thus, the phenomenological stochasticity of the bias should be accounted for by nonlocal features of galaxy formation processes. Based on mathematical arguments, we show that there are simple relations between biasing and nonlocality on linear scales of density fluctuations, and that the stochasticity in Fourier space does not exist on linear scales under a certain condition, even if the galaxy formation itself is a complex nonlinear and nonlocal precess. The stochasticity in real space, however, arise from the scale-dependence of bias parameter, $b$. As examples, we derive the stochastic bias parameters of simple nonlocal models of galaxy formation, i.e., the local Lagrangian bias models, the cooperative model, and the peak model. We show that the stochasticity in real space is also weak, except on the scales of nonlocality of the galaxy formation. Therefore, we do not have to worry too much about the stochasticity on linear scales, especially in Fourier space, even if we do not know the details of galaxy formation process.

The three-dimensional distribution of astronomical objects observed in redshift space significantly differs from the true distribution since the distance to each object cannot be determined by its redshift $z$ only; for $z \ll 1$ the peculiar velocity field contaminates the true recession velocity of the Hubble flow, while the true distance for objects at $z < 1$ sensitively depends on the (unknown and thus assumed) cosmological parameters. This hampers the effort to understand the true distribution of large-scale structure of the universe. In addition, all cosmological observations are carried out on a light-cone, the null hypersurface of an observer at $z=0$. This implies that their intrinsic properties and clustering statistics should change even within the survey volume. Therefore a proper comparison taking account of the light-cone effect is important to extract any cosmological information from redshift catalogues, especially for $z < 1$. We present recent theoretical development on the two effects -- the cosmological light-cone effect and the cosmological redshift-space distortion -- which should play key roles in observational cosmology in the 21st century.

Using a novel two-dimensional coordinate system, we have derived a particularly simple way to express the redshift distortions in galaxy redshift surveys with arbitrary geometry in closed form. This method provides an almost ideal way to measure the value of $\beta=\Omega_0^{0.6}/b$ in wide area surveys, since all pairs in the survey can be used for the analysis. In the limit of small angles, this result straightforwardly reduces to the plane-parallel approximation. This expansion can also be used together with more sophisticated methods such as for the calculation of Karhunen-Loeve eigenvectors in redshift space for an arbitrary survey geometry. Therefore, these results should provide for more precise methods in which to measure the large scale power spectrum and the value of $\beta$ simultaneously.

We have evaluated a systematic effect on counts-in-cells analysis of deep, wide-field galaxy catalogues induced by the evolution of clustering within the survey volume. A multiplicative correction factor is explicitly presented, which can be applied after the higher order correlation functions have been extracted in the usual way, without taking into account the evolution. The general theory of this effect combined with the ansatz describing the non-linear evolution of clustering in simulations enables us to estimate the magnitude of the correction factor in different cosmologies. In a series of numerical calculations assuming an array of cold dark matter models, it is found that, as long as galaxies are unbiased tracers of underlying density field, the effect is relatively small ($ \simeq 10%$) for the shallow surveys ($z <0.2$), while it becomes significant (order of unity) in deep surveys ($z \sim1$). Depending on the scales of interest, the required correction is comparable to or smaller than the expected errors of on-going wide-field galaxy surveys such as the SDSS and 2dF. Therefore at present, the effect has to be taken into account for high precision measurements at very small scales only, while in future deep surveys it amounts to a significant correction.

The observed two-point correlation functions of galaxies in redshift space become anisotropic due to the geometry of the universe as well as due to the presence of the peculiar velocity field. On the basis of linear perturbation theory, we expand the induced anisotropies of the correlation functions with respect to the redshift $z$, and obtain analytic formulae to infer the deceleration parameter $q_0$, the density parameter $\Omega_0$ and the derivative of the bias parameter $d\ln b/dz$ at $z=0$ in terms of the observable statistical quantities. The present method does not require any assumption of the shape and amplitude of the underlying fluctuation spectrum, and thus can be applied to future redshift surveys of galaxies including the Sloan Digital Sky Survey. We also evaluate quantitatively the systematic error in estimating the value of $\beta_0 \equiv \Omega_0^{0.6}/b$ from a galaxy redshift survey on the basis of a conventional estimator for $\beta_0$ which neglects both the geometrical distortion effect and the time evolution of the parameter $\beta(z)$. If the magnitude limit of the survey is as faint as 18.5 (in B-band) as in the case of the Sloan Digital Sky Survey, the systematic error ranges between -20% and 10% depending on the cosmological parameters. Although such systematic errors are smaller than the statistical errors in the current surveys, they will dominate the expected statistical error for future surveys.

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.

Evolution of genus density is calculated from Gaussian initial conditions using Zel'dovich approximation. A new approach is introduced which formulates the desired quantity in a rotationally invariant manner. It is shown that normalized genus density does not depend on the initial spectral shape but is a function of the fluctuation amplitude only.

We discuss the two potentially important effects which should be taken into account in the analysis of the cosmological redshift-space distortion especially at high redshifts; the effect of inhomogeneities in the light propagation and the evolution of bias. Although such inhomogeneities affect significantly the degree of the cosmological redshift distortion, one can determine the density parameter $\Omega_0$ from the distortion pattern at $z=0$, and the cosmological constant $\lambda_0$ from that at $z = (1\sim 2)$. In particular we show that in low-density ($\Omega_0\ll 1$) universes one can determine the value of $\lambda_0$ and even the evolution of the bias parameter almost insensitively to the degree of the clumpiness of the matter distribution in the universe.

We propose cosmological redshift-space distortion of correlation functions of galaxies and quasars as a probe of both the density parameter $\Omega_0$ and the cosmological constant $\lambda_0$. In particular, we show that redshift-space distortion of quasar correlation functions at $z\sim2$ can in principle set a constraint on the value of $\lambda_0$. This is in contrast to the popular analysis of galaxy correlation functions in redshift space which basically determines $\Omega_0^{0.6}/b$, where $b$ is the bias parameter, but is insensitive to $\lambda_0$. For specific applications, we present redshift-space distortion of correlation functions both in cold dark matter models and in power-law correlation function models, and discuss the extent to which one can discriminate between the different $\lambda_0$ models.

The Eulerian and Lagrangian second-order perturbation theories are solved explicitly in closed forms in $\Omega \neq 1$ and $\Lambda \neq 0$ {}Friedmann-Lema\^{\i}tre models. I explicitly write the second-order theories in terms of closed one-dimensional integrals. In cosmologically interested cases ($\Lambda = 0$ or $\Omega + \lambda = 1$), they reduce to elementary functions or hypergeometric functions. For arbitrary $\Omega$ and $\Lambda$, I present accurate fitting formula which are sufficient in practice for the observational cosmology. It is reconfirmed for generic $\Omega$ and $\Lambda$ of interest that second-order effect only weakly depends on these parameters.

The genus statistics is studied using large $N$-body simulations for several cosmological models. We consider the effects of nonlinear gravitational evolution, smoothing the particle data in fully nonlinear regime, and the redshift-space distortion on the genus curve. Detailed comparison between the theoretical prediction in weakly nonlinear theory and the appropriate simulation results shows that the analytic formula describes the behavior of genus in weakly nonlinear regime fairly accurately. We also find that the redshift-space distortion on genus statistics is small in linear and weakly nonlinear regimes. We conclude that if weakly nonlinear theory and direct numerical simulations are combined, the normalized genus curve $G(\nu)/G(0)$ is a powerful tool to directly check the random-Gaussian paradigm of the origin of the large-scale structure in the universe.

As a statistical measure to quantify the topological structure of the large-scale structure in the universe, the genus number is calculated for a number of non-Gaussian distributions in which the density field is characterized by a nontrivial function of some Gaussian-distributed random numbers. As a specific example, the formulae for the lognormal and the chi-square distributions are derived and compared with the results of $N$-body simulations together with the previously known formulae for the Gaussian distribution and second-order perturbation theory. It is shown that the lognormal formula fits most of the simulation data the best.

The gravitational evolution of the genus and other statistics of isodensity contours of the density field is derived analytically in a weakly nonlinear regime using second-order perturbation theory. The effect of final smoothing in perturbation theory on the statistics of isodensity contours is also evaluated. The resulting analytic expression for the genus is compared with $N$-body numerical simulations and exhibits a good agreement.

We develop technique concerning statistical analyses of the large-scale structure of the universe in the presence of biasing in the structure formation. We formulate the diagrammatic method to calculate the correlation functions of the nonlocally biased field from a generally non-Gaussian density field. Our method is based on generalized Wiener-Hermite expansion of the density field. The present formalism has not only general applicability, but also practical significance, too. To show the effectiveness of this method, we revisit the problems on biasing that have been considered previously, i.e., various approximations for peak statistics and hierarchical structure of correlation functions of locally biased field. Further analyses which have not been possible so far on these problems can be performed by our formalism. The gravitational evolution of primordial fluctuation or other structure-forming processes can be described by nonlocal biasing, so can be treated by our new formalism in principle.

The peculiar velocities of galaxies distort the clustering pattern anisotropically in redshift space. This effect on the statistics of isodensity contours is examined by linear theory. The statistics considered in this paper are the three- and two-dimensional genus of isodensity contours, the area of isodensity contours, the length of isodensity contours in the 2-dimensional slice and the level crossing statistics on the line. We find that all these statistics in redshift space as functions of density threshold of contours have the same shape as in real space. The redshift space distortion affects only amplitudes of these statistics. The three-dimensional genus and the area hardly suffer from the redshift space distortion for $0 \leq \Omega b^{-5/3} \leq 1$, where $b$ is a linear bias parameter. The other statistics are defined in one- or two-dimensional slices of the sample volume and depend on the direction of these slices relative to the line of sight. This dependence on direction of these statistics provides a way to determine the density parameter of the universe.

The gravitational evolution of the genus of the density field in large-scale structure is analytically studied in a weakly nonlinear regime using second-order perturbation theory. Weakly nonlinear evolution produces asymmetry in the symmetric genus curve for Gaussian initial density field. The effect of smoothing the density field in perturbation theory on the genus curve is also evaluated and gives the dependence of the asymmetry of the genus curve on spectra of initial fluctuations.