We present a Next-to-next-to-leading (fifth or NNLO) order calculation for the covariance matrix of the matter power spectrum, taking into account the effect of survey window functions. Using the grid-based calculation scheme for the standard perturbation theory, GridSPT, we quickly generate multiple realizations of the nonlinear density fields to fifth order in perturbation theory, then estimate the power spectrum and the covariance matrix from the sample. To the end, we have obtained the non-Gaussian covariance originated from the one-loop trispectrum without explicitly computing the trispectrum. By comparing the GridSPT calculations with the N-body results, we show that NNLO GridSPT result reproduces the N-body results on quasi-linear scales where SPT accurately models nonlinear matter power spectrum. Incorporating the survey window function effect to GridSPT is rather straightforward, and the resulting NNLO covariance matrix also matches well with the N-body results.

We provide a new fitting formula of the matter bispectrum in the nonlinear regime calibrated by high-resolution cosmological $N$-body simulations of $41$ cold dark matter ($w$CDM, $w=$ constant) models around the Planck 2015 best-fit parameters. As the parameterization in our fitting function is similar to that in Halofit, our fitting is named BiHalofit. The simulation volume is sufficiently large ($> 10 \, {\rm Gpc}^3$) to cover almost all measurable triangle bispectrum configurations in the universe. The function is also calibrated using one-loop perturbation theory at large scales ($k<0.3 \, h \, {\rm Mpc}^{-1}$). Our formula reproduced the matter bispectrum to within $10 \, (15) \, \%$ accuracy in the Planck 2015 model at wavenumber $k< 3 \, (10) \, h \, {\rm Mpc}^{-1}$ and redshifts $z=0-3$. The other $40$ $w$CDM models obtained poorer fits, with accuracy approximating $20 \, \%$ at $k<3 \, h \, {\rm Mpc}^{-1}$ and $z=0-1.5$ (the deviation includes the $10 \, \%$-level sample variance of the simulations). We also provide a fitting formula that corrects the baryonic effects such as radiative cooling and active galactic nucleus feedback, using the latest hydrodynamical simulation IllustrisTNG. We demonstrate that our new formula more accurately predicts the weak-lensing bispectrum than the existing fitting formulas. This formula will assist current and future weak-lensing surveys and cosmic microwave background lensing experiments. Numerical codes of the formula are available, written in Python, C and Fortran.

The observed galaxy distribution via galaxy redshift surveys appears distorted due to redshift-space distortions (RSD). While one dominant contribution to RSD comes from the Doppler effect induced by the peculiar velocity of galaxies, the relativistic effects, including the gravitational redshift effect, are recently recognized to give small but important contributions. Such contributions lead to an asymmetric galaxy clustering along the line of sight, and produce non-vanishing odd multipoles when cross-correlating between different biased objects. However, non-zero odd multipoles are also generated by the Doppler effect beyond the distant-observer approximation, known as the wide-angle effect, and at quasi-linear scales, the interplay between wide-angle and relativistic effects becomes significant. In this paper, based on the formalism developed by Taruya et al., we present a quasi-linear model of the cross-correlation function taking a proper account of both the wide-angle and gravitational redshift effects, as one of the major relativistic effects. Our quasi-linear predictions of the dipole agree well with simulations even at the scales below $20\,h^{-1}\,$Mpc, where non-perturbative contributions from the halo potential play an important role, flipping the sign of the dipole amplitude. When increasing the bias difference and redshift, the scale where the sign flip happens is shifted to a larger scale. We derive a simple approximate formula to quantitatively account for the behaviors of the sign flip.

We show that it is possible to build effective matter density power spectra in tomographic cosmic shear observations that exhibit the Baryonic Acoustic Oscillations (BAO) features once a nulling transformation has been applied to the data. The precision with which the amplitude and position of these features can be reconstructed is quantified in terms of sky coverage, intrinsic shape noise, median source redshift and number density of sources. BAO detection in Euclid or LSST like wide surveys will be possible with a modest signal-to-noise ratio. It would improve dramatically for slightly deeper surveys.

We show that the spatial correlation of the intrinsic alignments (IAs) of galaxies, measured in galaxy redshift surveys, offers a precision route to improve the geometrical and dynamical constraints on cosmology. The IA has been treated as a contaminant against cosmological probes such as weak gravitational lensing experiments. However, the large-scale correlation of IAs is expected to follow the coherent large-scale matter inhomogeneities. Here, making use of its anisotropic nature, we show that the large-scale IA correlations help to improve the measurements of the geometric distances and growth of structure. In combination with the conventional galaxy clustering statistics, we find that constraints on equation-of-state parameter for dark energy and Hubble parameter can be tighter than those from the clustering statistics alone by a factor of more than 1.5.

Investigating intrinsic alignments (IAs) of galaxy shapes is important not only to constrain cosmological parameters unbiasedly from gravitational lensing but also to extract cosmological information complimentary to galaxy clustering analysis. We derive simple and useful formulas for the various IA statistics, including the intrinsic ellipticity--ellipticity correlation, the gravitational shear--intrinsic ellipticity correlation, and the velocity-intrinsic ellipticity correlation functions. The angular dependence of each statistic is explicitly given, namely the angle between the line-of-sight direction and the separation vector of two points. It thus allows us to analyze anisotropies of baryon acoustic oscillations encoded in the IA statistics, and we can extract the maximum cosmological information using the Alcock-Pacyznski and redshift-space distortion effects. We also provide these formulas for the intrinsic ellipticities decomposed into E and B modes.

Using the motion of accreting particles onto halos in cosmological $N$-body simulations, we study the radial phase-space structures of cold dark matter (CDM) halos. In CDM cosmology, formation of virialized halos generically produces radial caustics, followed by multi-stream flows of accreted dark matter inside the halos. In particular, the radius of the outermost caustic called the splashback radius exhibits a sharp drop in the slope of the density profile. Here, we focus on the multi-stream structure of CDM halos inside the splashback radius. To analyze this, we use and extend the SPARTA algorithm developed by Diemer. By tracking the particle trajectories accreting onto the halos, we count their number of apocenter passages, which is then used to reveal the multi-stream flows of the dark matter particles. The resultant multi-stream structure in radial phase space is compared with the prediction of the self-similar solution by Fillmore & Goldreich for each halo. We find that $\sim30\%$ of the simulated halos satisfy our criteria to be regarded as being well fitted to the self-similar solution. The fitting parameters in the self-similar solution characterizes physical properties of the halos, including the mass accretion rate and the size of the outermost caustic (i.e., the splashback radius). We discuss in detail the correlation of these fitting parameters and other measures directly extracted from the $N$-body simulation.

The power spectrum response function of the large-scale structure of the Universe describes how the evolved power spectrum is modified by a small change in initial power through non-linear mode coupling of gravitational evolution. It was previously found that the response function for the coupling from small to large scales is strongly suppressed in amplitude, especially at late times, compared to predictions from perturbation theory (PT) based on the single-stream approximation. One obvious explanation for this is that PT fails to describe the dynamics beyond shell-crossing. We test this idea by comparing measurements in $N$-body simulations to prescriptions based on PT but augmented with adaptive smoothing to account for the formation of non-linear structures of various sizes in the multi-stream regime. We first start with one-dimensional (1D) cosmology, where the Zel'dovich approximation provides the exact solution in the single stream regime. Similarly to the 3D case, the response function of the large-scale modes exhibits a strong suppression in amplitude at small scales which cannot be explained by the Zel'dovich solution alone. However, by performing adaptive smoothing of initial conditions to identify haloes of different sizes and solving approximately post-collapse dynamics in the 3-streams regime, agreement between theory and simulations drastically improves. We extend our analyses to the 3D case using PINOCCHIO algorithm, in which similar adaptive smoothing is implemented on the Lagrangian PT fields to identify haloes and is combined with a spherical halo prescription to account for post-collapse dynamics. Again, a suppression is found in the coupling between small- and large-scale modes and the agreement with simulations is improved.

There is a growing interest of using the intrinsic alignment (IA) of galaxy images as a tool to extract cosmological information complimentary to galaxy clustering analysis. Recently, Okumura & Taruya derived useful formulas for the intrinsic ellipticity--ellipticity correlation, the gravitational shear--intrinsic ellipticity correlation, and the velocity--intrinsic ellipticity correlation functions based on the linear alignment (LA) model. In this paper, using large-volume $N$-body simulations, we measure these alignment statistics in real and redshift space and compare them to the LA and nonlinear alignment model predictions. We find that anisotropic features of baryon acoustic oscillations in the IA statistics can be accurately predicted by our models. The anisotropy due to redshift-space distortions (RSDs) is also well described in the large-scale limit. Our results indicate that one can extract the cosmological information encoded in the IA through the Alcock-Paczynski and RSD effects.

Redshift-space distortions (RSD) in galaxy redshift surveys generally break both the isotropy and homogeneity of galaxy distribution. While the former aspect is particularly highlighted as a probe of growth of structure induced by gravity, the latter aspect, often quoted as wide-angle RSD but ignored in most of the cases, will become important and critical to account for as increasing the statistical precision in next-generation surveys. However, the impact of wide-angle RSD has been mostly studied using linear perturbation theory. In this paper, employing the Zel'dovich approximation, i.e., first-order Lagrangian perturbation theory for gravitational evolution of matter fluctuations, we present a quasi-linear treatment of wide-angle RSD, and compute the cross-correlation function. The present formalism consistently reproduces linear theory results, and can be easily extended to incorporate relativistic corrections (e.g., gravitational redshift).

The kinematic Sunyaev-Zel'dovich effect enables us to directly probe the density-weighted velocity field up to very large cosmic scales. We investigate the effects of intrinsic alignments (IAs) of dark-matter halo shapes on cosmic density and velocity fields on such large scales. In the literature IAs have been detected up to $\sim 100~h^{-1}{\rm Mpc}$ using the gravitational shear-intrinsic ellipticity correlation and the alignment correlation function. In this paper we introduce the corresponding various velocity statistics: the (density-weighted) velocity-intrinsic ellipticity correlation as well as the alignment pairwise infall momentum, momentum correlation function, and density-weighted pairwise velocity dispersion. We derive theoretical expressions for these velocity alignment statistics for the first time based on the assumption that the density fluctuation is a Gaussian random field. Using large-volume, high-resolution $N$-body simulations, we measure the alignment statistics of density and velocity fields. The behaviors of IAs in the velocity statistics are similar to those in the density statistics, except that the halo orientations are aligned with the velocity field up to scales larger than those with the density field, $r\gg 100~h^{-1}{\rm Mpc}$, because of a factor of the wave number in the linear relation between the density and velocity fields in Fourier space, $v\propto \delta/k$. We show that the detected IAs of the velocity field can be well predicted by the linear alignment model. We also demonstrate that the baryon acoustic oscillation features can be detected in both the conventional and alignment velocity statistics. Our results indicate that observations of IAs with the velocity field on large scales can provide additional information on cosmological models, complementary to those with the density field.

Direct detection of gravitational waves (GWs) from compact binary systems suggests that the merger rate of such events is large, and the sum of their GWs can be viewed as stochastic signals. Because of its random nature, cross-correlating the signals from multiple detectors is essential to disentangle the GWs from instrumental noise. However, the global magnetic fields in the Earth-ionosphere cavity produce the environmental disturbances at low-frequency bands, known as Schumann resonances, and coupled with GW detectors, they potentially contaminate the stochastic GW signal as a correlated noise. Previously, we have presented a simple analytical model to estimate its impact on the detection of stochastic GWs. Here, extending the analysis to further take account of the effects of anisotropic lightning source distributions, we present a comprehensive study of the impact of correlated magnetic noise at low-frequency bands, including non-tensor-type GWs, as well as circularly polarized tensor-type GWs. We find that as opposed to a naive expectation, the impact of correlated magnetic noise does not always increase with anisotropies in the lighting source distribution. Even in the presence of large anisotropies, there is a robust detector pair for which the amplitude of correlated magnetic noise becomes comparable to or well below detectable amplitude of stochastic GWs. The results indicate that the properties of the correlated magnetic noise depend crucially on both the geometrical and geographical setup of the detector's pair, and Virgo and KAGRA would be potentially the most insensitive detector pair against the correlated magnetic for both tensor- and non-tensor-type stochastic GWs.

Recent studies have demonstrated that {\em secondary} non-Gaussianity induced by gravity will be detected with a high signal-to-noise (S/N) by future and even by on-going weak lensing surveys. One way to characterise such non-Gaussianity is through the detection of a non-zero three-point correlation function of the lensing convergence field, or of its harmonic transform, the bispectrum. A recent study analysed the properties of the squeezed configuration of the bispectrum, when two wavenumbers are much larger than the third one. We extend this work by estimating the amplitude of the (reduced) bispectrum in four generic configurations, i.e., {\em squeezed, equilateral, isosceles} and {\em folded}, and for four different source redshifts $z_s=0.5,1.0,1.5,2.0$, by using an ensemble of all-sky high-resolution simulations. We compare these results against theoretical predictions. We find that, while the theoretical expectations based on widely used fitting functions can predict the general trends of the reduced bispectra, a more accurate theoretical modelling will be required to analyse the next generation of all-sky weak lensing surveys. The disagreement is particularly pronounced in the squeezed limit.

We present robust constraints on the stochastic gravitational waves (GWs) at Mpc scales from the cosmic microwave background (CMB) data. CMB constraints on GWs are usually characterized as the tensor-to-scalar ratio, assuming specifically a power-law form of the primordial spectrum, and are obtained from the angular spectra of CMB. Here, we relax the assumption of the power-law form, and consider to what extent one can constrain a monochromatic GW at shorter wavelengths. Previously, such a constraint has been derived at the wavelengths larger than the resolution scale of the CMB measurements, typically above $100$Mpc (below $10^{-16}$Hz in frequency). However, GWs whose wavelength is much shorter than $100$Mpc can imprint a small but non-negligible signal on CMB anisotropies at observed angular scales, $\ell<1000$. Here, using the CMB temperature, polarization, and lensing data set, we obtain the best constraints to date at $10^{-16}-10^{-14}$Hz of the GWs produced before the time of decoupling, which are tighter than those derived from the astrometric measurements and upper bounds on extra radiations. In the future, the constraints on GWs at Mpc scales will be further improved by several orders of magnitude with the precision $B$-mode measurement on large scales, $\ell<100$.

We study the accuracy with which cosmological parameters can be determined from real space power spectrum of matter density contrast at weakly nonlinear scales using analytical approaches. From power spectra measured in $N$-body simulations and using Markov chain Monte-Carlo technique, the best-fitting cosmological input parameters are determined with several analytical methods as a theoretical template, such as the standard perturbation theory, the regularized perturbation theory, and the effective field theory. We show that at redshift 1, all two-loop level calculations can fit the measured power spectrum down to scales $k \sim 0.2 \, h \, \mathrm{Mpc}^{-1}$ and cosmological parameters are successfully estimated in an unbiased way. Introducing the Figure of bias (FoB) and Figure of merit (FoM) parameter, we determine the validity range of those models and then evaluate their relative performances. With one free parameter, namely the damping scale, the regularized perturbation theory is found to be able to provide the largest FoM parameter while keeping the FoB in the acceptance range.

Cosmic microwave background (CMB) lensing is an integrated effect whose kernel is greater than half the peak value in the range $1<z<5$. Measuring this effect offers a powerful tool to probe the large-scale structure of the Universe at high redshifts. With the increasing precision of ongoing CMB surveys, other statistics than the lensing power spectrum, in particular the lensing bi-spectrum, will be measured at high statistical significance. This will provide ways to improve the constraints on cosmological models and lift degeneracies. Following on an earlier paper, we test analytical predictions of the CMB lensing bi-spectrum against full-sky lensing simulations, and discuss their validity and limitation in detail. The tree-level prediction of perturbation theory agrees with the simulation only up to $\ell\sim 200$, but the one-loop order allows capturing the simulation results up to $\ell\sim 600$. We also show that analytical predictions based on fitting formulas for the matter bi-spectrum agree reasonably well with simulation results, although the precision of the agreement depends on the configurations and scales considered. For instance, the agreement is at the $10\%$-level for the equilateral configuration at multipoles up to $\ell\sim2000$, but the difference in the squeezed limit raises to more than a factor of two at $\ell\sim2000$. This discrepancy appears to come from limitations in the fitting formula of the matter bi-spectrum. We also find that the analytical prediction for the post-Born correction to the bi-spectrum is in good agreement with the simulation. We conclude by discussing the bi-spectrum prediction in some theories of modified gravity.

We consider the growth of primordial dark matter halos seeded by three crossed initial sine waves of various amplitudes. Using a Lagrangian treatment of cosmological gravitational dynamics, we examine the convergence properties of a high-order perturbative expansion in the vicinity of shell-crossing, by comparing the analytical results with state-of-the-art high resolution Vlasov-Poisson simulations. Based on a quantitative exploration of parameter space, we study explicitly for the first time the convergence speed of the perturbative series, and find, in agreement with intuition, that it slows down when going from quasi one-dimensional initial conditions (one sine wave dominating) to quasi triaxial symmetry (three sine waves with same amplitude). In most cases, the system structure at collapse time is, as expected, very similar to what is obtained with simple one-dimensional dynamics, except in the quasi-triaxial regime, where the phase-space sheet presents a velocity spike. In all cases, the perturbative series exhibits a generic convergence behavior as fast as an exponential of a power-law of the order of the expansion, allowing one to numerically extrapolate it to infinite order. The results of such an extrapolation agree remarkably well with the simulations, even at shell-crossing.

The apparent distribution of large-scale structures in the universe is sensitive to the velocity/potential of the sources as well as the potential along the line-of-sight through the mapping from real space to redshift space (redshift-space distortions, RSD). Since odd multipoles of the halo cross-correlation function vanish when considering standard Doppler RSD, the dipole is a sensitive probe of relativistic and wide-angle effects. We build a catalogue of ten million haloes (Milky-Way size to galaxy-cluster size) from the full-sky light-cone of a new "RayGalGroupSims" N-body simulation which covers a volume of ($2.625~h^{-1}$Gpc)$^3$ with $4096^3$ particles. Using ray-tracing techniques, we find the null geodesics connecting all the sources to the observer. We then self-consistently derive all the relativistic contributions (in the weak-field approximation) to RSD: Doppler, transverse Doppler, gravitational, lensing and integrated Sachs-Wolfe. It allows us, for the first time, to disentangle all contributions to the dipole from linear to non-linear scales. At large scale, we recover the linear predictions dominated by a contribution from the divergence of neighbouring line-of-sights. While the linear theory remains a reasonable approximation of the velocity contribution to the dipole at non-linear scales it fails to reproduce the potential contribution below $30-60~h^{-1}$Mpc (depending on the halo mass). At scales smaller than $\sim 10~h^{-1}$Mpc, the dipole is dominated by the asymmetry caused by the gravitational redshift. The transition between the two regimes is mass dependent as well. We also identify a new non-trivial contribution from the non-linear coupling between potential and velocity terms.

Gravity-induced non-Gaussianity in the large-scale structure of the Universe, characterised by higher-order statistics such as the bispectrum (three-point cumulant), is expected to contain rich cosmological information. A measurement of the bispectrum will not only improve the cosmological constraints, but also give us the possibility to probe gravity on cosmological scales. In this paper, we present a framework to numerically calculate the one-loop matter bispectrum based on standard perturbation theory (SPT). This approach allows general modifications to the standard $\Lambda$CDM model to be easily implemented. We demonstrate the performance of the bispectrum calculation in three representative cases, namely the Vainshtein-screened Dvali-Gabadadze-Porrati (DGP) model, the chameleon-screened Hu-Sawicki $f(R)$ model and the phenomenological dark scattering (DS) momentum-exchange model. The predicted bispectra are then compared with measured results from a set of cosmological $N$-body simulations, and the impact of possible systematics arising from simplified or approximate treatments in the perturbative calculation is studied in detail. We find that the one-loop bispectrum calculation offers significantly more information on general screening and momentum exchange effects than the leading-order bispectrum calculation. Further, the accuracy of the one-loop prediction is shown to be comparable to non-linear fitting formulas over a wide range of wavenumbers ($k\lesssim0.3\,h$$\mbox{Mpc}^{-1}$) even at lower redshifts, $z\lesssim 1$.

Perturbation theory (PT) calculation of large-scale structure has been used to interpret the observed non-linear statistics of large-scale structure at the quasi-linear regime. In particular, the so-called standard perturbation theory (SPT) provides a basis for the analytical computation of the higher-order quantities of large-scale structure. Here, we present a novel, grid-based algorithm for the SPT calculation, hence named GridSPT, to generate the higher-order density and velocity fields from a given linear power spectrum. Taking advantage of the Fast Fourier Transform, the GridSPT quickly generates the nonlinear density fields at each order, from which we calculate the statistical quantities such as non-linear power spectrum and bispectrum. Comparing the density fields (to fifth order) from GridSPT with those from the full N-body simulations in the configuration space, we find that GridSPT accurately reproduces the N-body result on large scales. The agreement worsens with smaller smoothing radius, particularly for the under-dense regions where we find that 2LPT (second-order Lagrangian perturbation theory) algorithm performs well.

Cosmological structures grow differently in theories of gravity which are modified as compared to Einstein's General relativity (GR). Cosmic microwave background (CMB) fluctuation patterns at the last scattering surface are lensed by these structures along the photon path to the observer. The observed CMB pattern therefore keeps trace of the growth history of structures. We show that observations of the CMB lensing bi-spectrum offer an interesting way to constrain deviations from GR in a broad class of scalar-tensor theories of gravity called "beyond Horndeski". We quantify how the constraints on generic parameters describing the deviations from GR depend on the effective multipole range of the analysis. Our results further indicate that an accurate nonlinear correction of the matter bi-spectrum in the modified gravity considered is necessary when the bi-spectrum is used to probe scales beyond a multipole $\ell_{\rm max} \gtrsim 1500$. We also found that the results are insensitive to details of the implementation of the screening mechanism, at very small scales. We finally demonstrate the potential of the lensing bi-spectrum to provide a blind reconstruction of the redshift evolution of our modified gravity parameters by combining the analysis of CMB and low-z source lensing data.

We investigate how the consistency relations of large-scale structures are modified when the initial density field is not Gaussian. We consider both scenarios where the primordial density field can be written as a nonlinear functional of a Gaussian field and more general scenarios where the probability distribution of the primordial density field can be expanded around the Gaussian distribution, up to all orders over $\delta_{L0}$. Working at linear order over the non-Gaussianity parameters $f_{\rm NL}^{(n)}$ or $S_n$, we find that the consistency relations for the matter density fields are modified as they include additional contributions that involve all-order mixed linear-nonlinear correlations $\langle \prod \delta_L \prod \delta \rangle$. We derive the conditions needed to recover the simple Gaussian form of the consistency relations. This corresponds to scenarios that become Gaussian in the squeezed limit. Our results also apply to biased tracers, and velocity or momentum cross-correlations.

The observed power spectrum in redshift space appears distorted due to the peculiar motion of galaxies, known as redshift-space distortions (RSD). While all the effects in RSD are accounted for by the simple mapping formula from real to redshift spaces, accurately modeling redshift-space power spectrum is rather difficult due to the non-perturbative properties of the mapping. Still, however, a perturbative treatment may be applied to the power spectrum at large-scales, and on top of a careful modeling of the Finger-of-God effect caused by the small-scale random motion, the redshift-space power spectrum can be expressed as a series of expansion which contains the higher-order correlations of density and velocity fields. In our previous work [JCAP 8 (Aug., 2016) 050], we provide a perturbation-theory inspired model for power spectrum in which the higher-order correlations are evaluated directly from the cosmological $N$-body simulations. Adopting a simple Gaussian ansatz for Finger-of-God effect, the model is shown to quantitatively describe the simulation results. Here, we further push this approach, and present an accurate power spectrum template which can be used to estimate the growth of structure as a key to probe gravity on cosmological scales. Based on the simulations, we first calibrate the uncertainties and systematics in the pertrubation theory calculation in a fiducial cosmological model. Then, using the scaling relations, the calibrated power spectrum template is applied to a different cosmological model. We demonstrate that with our new template, the best-fitted growth functions are shown to reproduce the fiducial values in a good accuracy of 1 \% at $k<0.18 \hompc$ for cosmologies with different Hubble parameters.

We present numerical measurements of the power spectrum response function of the gravitational growth of cosmic structures, defined as the functional derivative of the nonlinear spectrum with respect to the linear counterpart, based on $1,400$ cosmological simulations. We develop a simple analytical model based on a regularization of the standard perturbative calculation. Using the model prediction, we show that this function gives a natural way to interpolate the nonlinear power spectrum over cosmological parameter space from single or multi-step interpolations. We demonstrate that once an accurate numerical spectrum template is available for one (or a small number of) cosmological model(s), it doubles the range in $k$ for which percent level accuracy can be obtained even for large change in the cosmological parameters. The python package RESPRESSO we developed to make those predictions is publicly available.

After the first direct detection of gravitational waves (GW), detection of stochastic background of GWs is an important next step, and the first GW event suggests that it is within the reach of the second-generation ground-based GW detectors. Such a GW signal is typically tiny, and can be detected by cross-correlating the data from two spatially separated detectors if the detector noise is uncorrelated. It has been advocated, however, that the global magnetic fields in the Earth-ionosphere cavity produce the environmental disturbances at low-frequency bands, known as Schumann resonances, which potentially couple with GW detectors. In this paper, we present a simple analytical model to estimate its impact on the detection of stochastic GWs. The model crucially depends on the geometry of the detector pair through the directional coupling, and we investigate the basic properties of the correlated magnetic noise based on the analytic expressions. The model reproduces the major trend of the recently measured global correlation between the GW detectors via magnetometer. The estimated values of the impact of correlated noise also match those obtained from the measurement. Finally, we give an implication to the detection of stochastic GWs including upcoming detectors, KAGRA and LIGO India. The model suggests that LIGO Hanford-Virgo and Virgo-KAGRA pairs are possibly less sensitive to the correlated noise, and can achieve a better sensitivity to the stochastic GW signal in the most pessimistic case.

The large-scale matter distribution in the late-time Universe exhibits gravity-induced non-Gaussianity, and the bispectrum, three-point cumulant is expected to contain significant cosmological information. In particular, the measurement of the bispectrum helps to tighten the constraints on dark energy and modified gravity through the redshift-space distortions (RSD). In this paper, extending the work by Taruya, Nishimichi & Saito (2010, Phys.Rev.D 82, 063522), we present a perturbation theory (PT) based model of redshift-space matter bispectrum that can keep the non-perturbative damping effect under control. Characterizing this non-perturbative damping by a univariate function with single free parameter, the PT model of the redshift-space bispectrum is tested against a large set of cosmological $N$-body simulations, finding that the predicted monopole and quadrupole moments are in a good agreement with simulations at the scales of baryon acoustic oscillations (well beyond the range of agreement of standard PT). The validity of the univariate ansatz of the damping effect is also examined, and with the PT calculation at next-to-leading order, the fitted values of the free parameter is shown to consistently match those obtained from the PT model of power spectrum by Taruya, Nishimichi & Saito (2010).

We develop a new perturbation theory (PT) treatment that can describe gravitational dynamics of large-scale structure after shell-crossing in the one-dimensional cosmological case. Starting with cold initial conditions, the motion of matter distribution follows at early stages the single-stream regime, which can, in one dimension, be described exactly by the first-order Lagrangian perturbation, i.e. the Zel'dovich solution. However, the single-stream flow no longer holds after shell-crossing and a proper account of the multi-stream flow is essential for post-collapse dynamics. In this paper, extending previous work by Colombi (2015, MNRAS 446, 2902), we present a perturbative description for the multi-stream flow after shell-crossing in a cosmological setup. In addition, we introduce an adaptive smoothing scheme to deal with the bulk properties of phase-space structures. The filtering scales in this scheme are linked to the next-crossing time in the post-collapse region, estimated from our PT calculations. Our PT treatment combined with adaptive smoothing is illustrated in several cases. Predictions are compared to simulations and we find that post-collapse PT with adaptive smoothing reproduces the power spectrum and phase-space structures remarkably well even at small scales, where Zel'dovich solution substantially deviates from simulations.

In order to infer the impact of the small-scale physics to the large-scale properties of the universe, we use a series of cosmological $N$-body simulations of self-gravitating matter inhomogeneities to measure, for the first time, the response function of such a system defined as a functional derivative of the nonlinear power spectrum with respect to its linear counterpart. Its measured shape and amplitude are found to be in good agreement with perturbation theory predictions except for the coupling from small to large-scale perturbations. The latter is found to be significantly damped, following a Lorentzian form. These results shed light on validity regime of perturbation theory calculations giving a useful guideline for regularization of small scale effects in analytical modeling. Most importantly our result indicates that the statistical properties of the large-scale structure of the universe are remarkably insensitive to the details of the small-scale physics, astrophysical or gravitational, paving the way for the derivation of robust estimates of theoretical uncertainties on the determination of cosmological parameters from large-scale survey observations.

The first discovery of the gravitational wave (GW) event, GW150914, suggests a higher merger rate of black-hole (BH) binaries. If this is true, a number of BH binaries will be observed via the second-generation GW detectors, and the statistical properties of the observed BH binaries can be scrutinized. A naive but important question to ask is whether the spatial distribution of BH binaries faithfully traces the matter inhomogeneities in the Universe or not. Although the BH binaries are thought to be formed inside the galaxies in most of the scenarios, there is no observational evidence to confirm such a hypothesis. Here, we estimate how well the second-generation GW detectors can statistically confirm the BH binaries to be a tracer of the large-scale structure by looking at the auto- and cross-correlation of BH binaries with photometric galaxies and weak lensing measurements, finding that, with a 3 year observation, the $>3\sigma$ detection of non-zero signal is possible if the BH merger rate today is $\dot{n}_0\gtrsim100$ Gpc$^{-3}$yr$^{-1}$ and the clustering bias of BH binaries is $b_{\rm BH,0}\gtrsim1.5$.

We present a simple numerical scheme for perturbation theory (PT) calculations of large-scale structure. Solving the evolution equations for perturbations numerically, we construct the PT kernels as building blocks of statistical calculations, from which the power spectrum and/or correlation function can be systematically computed. The scheme is especially applicable to the generalized structure formation including modified gravity, in which the analytic construction of PT kernels is intractable. As an illustration, we show several examples for power spectrum calculations in $f(R)$ gravity and $\Lambda$CDM models.

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.

Gravitational waves (GWs) from compact binary stars at cosmological distances are promising and powerful cosmological probes, referred to as the GW standard sirens. With future GW detectors, we will be able to precisely measure source luminosity distances out to a redshift $z\sim5$. To extract cosmological information, previously proposed cosmological studies using the GW standard sirens rely on source redshift information obtained through an extensive electromagnetic follow-up campaign. However, the redshift identification is typically time consuming and rather challenging. Here, we propose a novel method for cosmology with the GW standard sirens free from the redshift measurements. Utilizing the anisotropies of the number density and luminosity distances of compact binaries originated from the large-scale structure, we show that, once GW observations will be well established in the future, (i) these anisotropies can be measured even at very high redshifts ($z\geq 2$), where the identification of the electromagnetic counterpart is difficult, (ii) the expected constraints on the primordial non-Gaussianity with the Einstein Telescope would be comparable to or even better than the other large-scale structure probes at the same epoch, and (iii) the cross-correlation with other cosmological observations is found to have high-statistical significance, providing additional cosmological information at very high redshifts.

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

The apparent anisotropies of the galaxy clustering in observable redshift space provide a unique opportunity to simultaneously probe cosmic expansion and gravity on cosmological scales via the Alcock--Paczynski effect and redshift-space distortions. While the improved theoretical models have been proposed and developed to describe the apparent anisotropic clustering at weakly non-linear scales, the applicability of these models is still limited in the presence of the non--perturbative smearing effect caused by the randomness of the relative velocities. Although the cosmological constraint from the anisotropic clustering will be certainly improved with a more elaborate theoretical model, we here consider an alternative way by using the statistical power of both the power spectrum and bispectrum at large scales. Based on the Fisher matrix analysis, we estimate the benefit of combining the power spectra and bispectra, finding that the constraints on the cosmic expansion and growth of structure will be improved by a factor of two. This compensates for the loss of constraining power using the power spectrum alone due to the randomness of the relative velocities.

We present a method of lensing reconstruction on the full sky, by extending the optimal quadratic estimator proposed by Okamoto & Hu (2003) to the case including the curl mode of deflection angle. The curl mode is induced by the vector and tensor metric perturbations, and the reconstruction of the curl mode would be a powerful tool to not only check systematics in the estimated gradient mode but also probe any vector and tensor sources. We find that the gradient and curl modes can be reconstructed separately, thanks to the distinctive feature in the parity symmetry between the gradient and curl modes. We compare our estimator with the flat-sky estimator proposed by Cooray et al (2005). Based on the new formalism, the expected signal-to-noise ratio of the curl mode produced by the primordial gravitational-waves and a specific model of cosmic strings are estimated, and prospects for future observations are discussed.

We discuss the future detectability of gravitational-wave induced lensing from high-sensitivity cosmic microwave background (CMB) experiments. Gravitational waves can induce a rotational component of the weak-lensing deflection angle, usually referred to as the curl mode, which would be imprinted on the CMB maps. Using the technique of reconstructing lensing signals involved in CMB maps, this curl mode can be measured in an unbiased manner, offering an independent confirmation of the gravitational waves complementary to the B-mode polarization experiments. Based on the Fisher matrix analysis, we first show that with the noise levels necessary to confirm the consistency relation for the primordial gravitational waves, the future CMB experiments will be able to detect the gravitational-wave induced lensing signals. For a tensor-to-scalar ratio of $r < 0.1$, even if the consistency relation is difficult to confirm with a high significance, the gravitational-wave induced lensing would be detected at more than $3\,\sigma$ significance level. Further, we point out that high-sensitivity experiments will be also powerful to constrain the gravitational waves generated after the recombination epoch. Compared to the B-mode polarization, the curl mode is particularly sensitive to gravitational waves generated at low redshifts ($z < 10$) with a low frequency ($k < 10^{-3}$ Mpc$^{-1}$), and it could give a much tighter constraint on their energy density $\Omega_{\rm GW}$ by more than three orders of magnitude.

The most widely used Galactic extinction map (Schlegel, Finkbeiner, & Davis 1998, SFD) is constructed assuming that the observed FIR fluxes entirely come from the Galactic dust. According to the earlier suggestion by Yahata et al. (2007), we consider how far-infrared (FIR) emission of galaxies affects the SFD map. We first compute the surface number density of SDSS DR7 galaxies as a function of the $r$-band extinction, $A_{r,\rm SFD}$. We confirm that the surface densities of those galaxies positively correlate with $A_{r,\rm SFD}$ for $A_{r,\rm SFD}<0.1$, as first discovered by Yahata et al. (2007) for SDSS DR4 galaxies. Next we construct an analytic model to compute the surface density of galaxies taking account of the contamination of their FIR emission. We adopt a log-normal probability distribution for the ratio of $100\mu {\rm m}$ and $r$-band luminosities of each galaxy, $y \equiv (\nu L)_{100\mu {\rm m}}/(\nu L)_r$. Then we search for the mean and r.m.s values of $y$ that fit the observed anomaly using the analytic model. The required values to reproduce the anomaly are roughly consistent with those measured from the stacking analysis of SDSS galaxies (Kashiwagi, Yahata, & Suto 2013). Due to the limitation of our statistical modeling, we are not yet able to remove the FIR contamination of galaxies from the extinction map. Nevertheless the agreement with the model prediction suggests that the FIR emission of galaxies is mainly responsible for the observed anomaly. While the corresponding systematic error in the Galactic extinction map is 0.1 to 1mmag, it is directly correlated with galaxy clustering, and thus needs to be carefully examined in precision cosmology.

Redshift-space distortions (RSD) offer an attractive method to measure the growth of cosmic structure on large scales, and combining with the measurement of the cosmic expansion history, it can be used as cosmological tests of gravity. With the advent of future galaxy redshift surveys aiming at precisely measuring the RSD, an accurate modeling of RSD going beyond linear theory is a critical issue in order to detect or disprove small deviations from general relativity (GR). While several improved models of RSD have been recently proposed based on the perturbation theory (PT), the framework of these models heavily relies on GR. Here, we put forward a new PT prescription for RSD in general modified gravity models. As a specific application, we present theoretical predictions of the redshift-space power spectra in f(R) gravity model, and compare them with N-body simulations. Using the PT template that takes into account the effects of both modifications of gravity and RSD properly, we successfully recover the fiducial model parameter in N-body simulations in an unbiased way. On the other hand, we found it difficult to detect the scale dependence of the growth rate in a model-independent way based on GR templates.

An explicit full nulling scheme for cosmic shear observations is presented. It makes possible the construction of shear maps from extended source distributions for which the lens distance distribution is restricted to a definite interval. Such a construction allows to build totally independent shear maps, at all scales, that can be taken advantage of to constrain background cosmological parameters and systematics using the full statistical power of cosmic shear observations. Another advantage of such construction is that, as the lens redshift distribution can be made arbitrarily narrow, scale mixing due to projection effects can be limited allowing controlled predictions on the large scale shear power spectrum from perturbation theory calculations.

Based on the multi-point propagator expansion, we present resummed perturbative calculations for cosmological power spectra and correlation functions in the context of modified gravity. In a wide class of modified gravity models that have a screening mechanism to recover general relativity (GR) on small scales, we apply the eikonal approximation to derive the governing equation for resummed propagator that partly includes the non-perturbative effect in the high-k limit. The resultant propagator in the high-k limit contains the new corrections arising from the screening mechanism as well as the standard exponential damping. We explicitly derive the expression for new high-k contributions in specific modified gravity models, and find that in the case of f(R) gravity for a currently constrained model parameter, the corrections are basically of the sub-leading order and can be neglected. Thus, in f(R) gravity, similarly to the GR case, we can analytically construct the regularized propagator that reproduces both the resummed high-k behavior and the low-k results computed with standard perturbation theory, consistently taking account of the nonlinear modification of gravity valid at large scales. With the regularized multi-point propagators, we give predictions for power spectrum and correlation function at one-loop order, and compare those with N-body simulations in f(R) gravity model. As an important application, we also discuss the redshift-space distortions and compute the anisotropic power spectra and correlation functions.

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.

Baryon acoustic oscillations (BAO), known as one of the largest cosmological objects, is now recognized as standard cosmological tool to measure geometric distances via the Alcock-Paczynski effect, by which the observed BAO exhibits characteristic anisotropies in addition to the redshift distortions. This implies that once we know the correct distances to the observed BAO, the tip points of baryon acoustic peaks in the anisotropic correlation function of galaxies, $\xi(\sigma,\pi)$, can form a great circle (hereafter 2D BAO circle) in the $\sigma$ and $\pi$ plane, where $\sigma$ and $\pi$ are the separation of galaxy pair parallel and perpendicular to the line-of-sight, respectively. This 2D BAO circle remains unchanged under the variations of the unknown galaxy bias and/or coherent motion, while it varies transversely and radially with respect to the variations of $D_A$ and $H^{-1}$, respectively. Hereby the ratio between transverse distance $D_A$ and the radial distance $H^{-1}$ reproduces the intrinsic shape of 2D BAO circle, which is {\it a priori} given by the known broadband shape of spectra. All BAO peaks of $\xi(\sigma,\pi)$ are precisely calculated with the improved theoretical model of redshift distortion. We test this broadband Alcock--Paczynski method using BOSS--like mock catalogues. The transverse and radial distances are probed in precision of several percentage fractional errors, and the coherent motion is observed to match with the fiducial values accurately.

The anisotropic galaxy clustering on large scales provides us with a unique opportunity to probe into the gravity theory through the redshift-space distortions (RSDs) and the Alcock-Paczynski effect. Using the multipole power spectra up to hexadecapole (ell=4), of the Luminous Red Galaxy (LRG) sample in the data release 7 (DR7) of the Sloan Digital Sky Survey II (SDSS-II), we obtain simultaneous constraints on the linear growth rate f, angular diameter distance D_A, and Hubble parameter H at redshift z = 0.3. For this purpose, we first extensively examine the validity of a theoretical model for the non-linear RSDs using mock subhalo catalogues from N-body simulations, which are constructed to match with the observed multipole power spectra. We show that the input cosmological parameters of the simulations can be recovered well within the error bars by comparing the multipole power spectra of our theoretical model and those of the mock subhalo catalogues. We also carefully examine systematic uncertainties in our analysis by testing the dependence on prior assumption of the theoretical model and the range of wavenumbers to be used in the fitting. These investigations validate that the theoretical model can be safely applied to the real data. Thus, our results from the SDSS DR7 LRG sample are robust including systematics of theoretical modeling; f(z = 0.3) sigma_8(z = 0.3) =0.49+-0.08(stat.)+-0.04(sys.), D_A (z = 0.3) =968+-42(stat.)+-17(sys.)[Mpc], H (z = 0.3) =81.7+-5.0(stat.)+-3.7(sys.)[km/s/Mpc]. We believe that our method to constrain the cosmological parameters using subhaloes catalogues will be useful for more refined samples like CMASS and LOWZ catalogues in the Baryon Oscillation Spectroscopic Survey in SDSS-III.

The angle between the stellar spin and the planetary orbit axes (spin-orbit angle) is supposed to carry valuable information on the initial condition of the planet formation and the subsequent migration history. Indeed current observations of the Rossiter- McLaughlin effect have revealed a wide range of spin-orbit misalignments for transiting exoplanets. We examine in detail the tidal evolution of a simple system comprising a Sun-like star and a hot Jupiter adopting the equilibrium tide and the inertial wave dissipation effects simultaneously. We find that the combined tidal model works as a very efficient realignment mechanism; it predicts three distinct states of the spin-orbit angle (i.e., parallel, polar, and anti-parallel orbits) for a while, but the latter two states eventually approach the parallel spin-orbit configuration. The intermediate spin-orbit angles as measured in recent observations are difficult to be achieved. Therefore the current model cannot reproduce the observed broad distribution of the spin-orbit angles, at least in its simple form. This indicates that the observed diversity of the spin-orbit angles may emerge from more complicated interactions with outer planets and/or may be the consequence of the primordial misalignment between the proto-planetary disk and the stellar spin, which requires future detailed studies.

The KOI-94 system is a closely-packed, multi-transiting planetary system discovered by the Kepler space telescope. It is known as the first system that exhibited a rare event called a "planet-planet eclipse (PPE)," in which two planets partially overlap with each other in their double-transit phase. In this paper, we constrain the parameters of the KOI-94 system with an analysis of the transit timing variations (TTVs). Such constraints are independent of the radial velocity (RV) analysis recently performed by Weiss and coworkers, and valuable in examining the reliability of the parameter estimate using TTVs. We numerically fit the observed TTVs of KOI-94c, KOI-94d, and KOI-94e for their masses, eccentricities, and longitudes of periastrons, and obtain the best-fit parameters including $m_{\rm c} = 9.4_{-2.1}^{+2.4} M_{\oplus}$, $m_{\rm d} = 52.1_{-7.1}^{+6.9} M_{\oplus}$, $m_{\rm e} = 13.0_{-2.1}^{+2.5} M_{\oplus}$, and $e \lesssim 0.1$ for all the three planets. While these values are mostly in agreement with the RV result, the mass of KOI-94d estimated from the TTV is significantly smaller than the RV value $m_{\rm d} = 106 \pm 11 M_{\oplus}$. In addition, we find that the TTV of the outermost planet KOI-94e is not well reproduced in the current modeling. We also present analytic modeling of the PPE and derive a simple formula to reconstruct the mutual inclination of the two planets from the observed height, central time, and duration of the brightening caused by the PPE. Based on this model, the implication of the results of TTV analysis for the time of the next PPE is discussed.

We present constraints on a cosmic string network with a measurement of weak gravitational lensing from CMB temperature map. The cosmic string network between observer and last scattering surface of CMB photons generates vector and/or tensor metric perturbations, and the deflection of CMB photons by these gravitational fields has curl mode which is not produced by the scalar metric perturbations. In this paper, we use the power spectrum of curl mode obtained from Planck to constrain the string tension, G\mu, and the reconnection probability, P. In demonstrating the parameter constraints with Planck curl mode, we also measure the lensing power spectrum from the Atacama Cosmology Telescope (ACT) 2008 season data, which have better angular resolution with lower instrumental noise on a much smaller chunk of the sky. Assuming P=1, the upper bound on tension is G\mu = 6.6\times 10^{-5} with 2\sigma (95% C.L.), using curl mode from Planck, which is weaker than that from the small-scale temperature power spectrum. For small values of P, however, the constraint from curl mode becomes tighter compared to that from temperature power spectrum. For P\lsim 10^{-2}, we obtain the constraint on the combination of the string parameters as G\mu P^{-1}\leq 3.4\times 10^{-5} at more than 2\sigma (95% C.L.).

We systematically derive full-sky formulae for the weak lensing power spectra generated by scalar, vector and tensor perturbations from the total angular momentum (TAM) method. Based on both the geodesic and geodesic deviation equations, we first give the gauge-invariant expressions for the deflection angle and Jacobi map as observables of the CMB lensing and cosmic shear experiments. We then apply the TAM method, originally developed in the theoretical studies of CMB, to a systematic derivation of the angular power spectra. The TAM representation, which characterizes the total angular dependence of the spatial modes projected along a line-of-sight, can carry all the information of the lensing modes generated by scalar, vector, and tensor metric perturbations. This greatly simplifies the calculation, and we present a complete set of the full-sky formulae for angular power spectra in both the E-/B-mode cosmic shear and gradient-/curl-mode lensing potential of deflection angle. Based on the formulae, we give illustrative examples of non-vanishing B-mode cosmic shear and curl-mode of deflection angle in the presence of the vector and tensor perturbations, and explicitly compute the power spectra.

We present a new approach to computing the matter density power spectrum, from large linear scales to small highly nonlinear scales. Instead of explicitly computing a partial series of high-order diagrams, as in perturbative resummation schemes, we embed the standard perturbation theory within a realistic nonlinear Lagrangian-space ansatz. We also point out that an "adhesion-like" regularization of the shell-crossing regime is more realistic than a "Zel'dovich-like" behavior, where particles freely escape to infinity. This provides a "cosmic web" power spectrum with good small-scale properties that provide a good matching with a halo model on mildly nonlinear scales. We obtain a good agreement with numerical simulations on large scales, better than 3% for $k\leq 1 h$Mpc$^{-1}$, and on small scales, better than 10% for $k \leq 10 h$Mpc$^{-1}$, at $z \geq 0.35$, which improves over previous methods.

Using a full implementation of resummed perturbation theory (PT) from a multi-point propagator expansion, we put forward new theoretical predictions for the two-point statistics of matter fluctuations in redshift space. The predictions consistently include PT corrections up to the two-loop order and are based on an improved prescription of the redshift-space distortions (RSD) that properly takes into account their non-Gaussian impact from a systematic low-k expansion. In contrast to the previous studies that partly used standard PT calculations, the present treatment is able to provide a consistent prediction for both power spectra and correlation functions. These results are compared with $N$-body simulations with which a very good agreement is found up to the quadrupole moment. The theoretical predictions for the hexadecapole moment of the power spectra are however found to significantly departs from the numerical results at low redshift. We examine this issue and found it to be likely related to an improper modeling of the RSD damping effects on which this moment shows large dependence.

The cosmological power spectrum of the coherent matter flow is measured exploiting an improved prescription for the apparent anisotropic clustering pattern in redshift space. New statistical analysis is presented to provide an optimal observational platform to link the improved redshift distortion theoretical model to future real datasets. The statistical power as well as robustness of our method are tested against 60 realizations of 8 Gpc/h^3 dark matter simulation maps mocking the precision level of upcoming wide--deep surveys. We showed that we can accurately extract the velocity power spectrum up to quasi linear scales of k~0.1 h/Mpc at z = 0.35 and up to k~0.15 h/Mpc at higher redshifts within a couple of percentage precision level. Our understanding of redshift space distortion is proved to be appropriate for precision cosmology, and our statistical method will guide us to righteous path to meet the real world.

We explore the properties of two-point cosmic propagators when Perturbation Theory (PT) loop corrections are consistently taken into account. We show in particular how the interpolation scheme proposed in arXiv:1112.3895 can be explicitly used up to two-loop order introducing the notion of regular parts for the contributing terms. Extending the one-loop results, we then derive and give semi analytical forms of the two-loop contributions for both the cosmic density and velocity propagators. These results are tested against numerical simulations and shown to significantly improve upon one-loop results for redshifts above 0.5. We found however that at lower redshift two-loop order corrections are too large partly due to a strong sensitivity of those terms to the small scale modes. We show that this dependence is expected to be even larger for higher order loop corrections both from theoretical investigations and numerical tests, the latter obtained with Monte Carlo evaluations of the three-loop contributions. This makes small-scale regularization schemes necessary for the use of such higher order corrections.

Based on a suite of state-of-the-art high-resolution $N$-body simulations, we revisit the so-called halofit model (Smith et al. 2003) as an accurate fitting formula for the nonlinear matter power spectrum. While the halofit model has been frequently used as a standard cosmological tool to predict the nonlinear matter power spectrum in a universe dominated by cold dark matter, its precision has been limited by the low-resolution of $N$-body simulations used to determine the fitting parameters, suggesting the necessity of improved fitting formula at small scales for future cosmological studies. We run high-resolution $N$-body simulations for 16 cosmological models around the Wilkinson Microwave Anisotropy Probe (WMAP) best-fit cosmological parameters (1, 3, 5, and 7 year results), including dark energy models with a constant equation of state. The simulation results are used to re-calibrate the fitting parameters of the halofit model so as to reproduce small-scale power spectra of the $N$-body simulations, while keeping the precision at large scales. The revised fitting formula provides an accurate prediction of the nonlinear matter power spectrum in a wide range of wavenumber ($k \leq 30h$\,Mpc$^{-1}$) at redshifts $0 \leq z \leq 10$, with 5% precision for $k\leq1 h$ Mpc$^{-1}$ at $0 \leq z \leq 10$ and 10% for $1 \leq k\leq 10 h$ Mpc$^{-1} $ at $0 \leq z \leq 3$. We discuss the impact of the improved halofit model on weak lensing power spectra and correlation functions, and show that the improved model better reproduces ray-tracing simulation results.

We study the observational signature of vector metric perturbations through the effect of weak gravitational lensing. In the presence of vector perturbations, the non-vanishing signals for B-mode cosmic shear and curl-mode deflection angle, which have never appeared in the case of scalar metric perturbations, naturally arise. Solving the geodesic and geodesic deviation equations, we drive the full-sky formulas for angular power spectra of weak lensing signals, and give the explicit expressions for E-/B-mode cosmic shear and gradient-/curl-mode deflection angle. As a possible source for seeding vector perturbations, we then consider a cosmic string network, and discuss its detectability from upcoming weak lensing and CMB measurements. Based on the formulas and a simple model for cosmic string network, we calculate the angular power spectra and expected signal-to-noise ratios for the B-mode cosmic shear and curl-mode deflection angle. We find that the weak lensing signals are enhanced for a smaller intercommuting probability of the string network, $P$, and they are potentially detectable from the upcoming cosmic shear and CMB lensing observations. For $P\sim 10^{-1}$, the minimum detectable tension of the cosmic string will be down to $G\mu\sim 5\times 10^{-8}$. With a theoretically inferred smallest value $P\sim 10^{-3}$, we could even detect the string with $G\mu\sim 5\times 10^{-10}$.

We present a specific prescription for the calculation of cosmological power spectra, exploited here at two-loop order in perturbation theory (PT), based on the multi-point propagator expansion. In this approach power spectra are constructed from the regularized expressions of the propagators that reproduce both the resummed behavior in the high-k limit and the standard PT results at low-k. With the help of N-body simulations, we show that such a construction gives robust and accurate predictions for both the density power spectrum and the correlation function at percent-level in the weakly non-linear regime. We then present an algorithm that allows accelerated evaluations of all the required diagrams by reducing the computational tasks to one-dimensional integrals. This is achieved by means of pre-computed kernel sets defined for appropriately chosen fiducial models. The computational time for two-loop results is then reduced from a few minutes, with the direct method, to a few seconds with the fast one. The robustness and applicability of this method are tested against the power spectrum cosmic emulator from which a wide variety of cosmological models can be explored. The fortran program with which direct and fast calculations of power spectra can be done, RegPT, is publicly released as part of this paper.

We present an investigation of spin-orbit angles for planetary system candidates reported by Kepler. By combining the rotational period $P_s$ inferred from the flux variation due to starspots and the projected rotational velocity $V\sin I_s$ and stellar radius obtained by a high resolution spectroscopy, we attempt to estimate the inclination $I_s$ of the stellar spin axis with respect to the line-of-sight. For transiting planetary systems, in which planetary orbits are edge-on seen from us, the stellar inclination $I_s$ can be a useful indicator of a spin-orbit alignment/misalignment. We newly conducted spectroscopic observations with Subaru/HDS for 15 KOI systems, whose lightcurves show periodic flux variations. After detailed analyses of their lightcurves and spectra, it turned out that some of them are binaries, or the flux variations are too coherent to be caused by starspots, probably representing ellipsoidal variations, and consequently we could constrain stellar inclinations $I_s$ for eight systems. Among them, KOI-262 and 280 are in good agreement with $I_s=90^\circ$ suggesting a spin-orbit alignment, while at least one system, KOI-261, shows a possible spin-orbit misalignment. We also obtain a small $I_s$ for KOI-1463, but the transiting companion seems to be a star rather than a planet. The results for KOI-257, 269, 367, and 974 are ambiguous, and can be explained with either misalignments or moderate differential rotation. Since our method can be applied to any system having starspots regardless of the planet size, future observations will allow for the expansion of the parameter space in which the spin-orbit relations are investigated.

Proposed space-based gravitational-wave (GW) detectors such as DECIGO and BBO will detect $\sim10^6$ neutron-star (NS) binaries and determine the luminosity distances to the binaries with high precision. Combining the luminosity distances with cosmologically-induced phase corrections on the GWs, cosmological expansion out to high redshift can be measured without the redshift determinations of host galaxies by electromagnetic observation and can be a unique probe for dark energy. This article is based on the results obtained in [1] where we investigated constraining power of the GW standard siren without redshift information on the equation of state of dark energy with future space-based GW detectors. We also compare the results with those obtained with other instruments and methods.

Proposed space-based gravitational-wave (GW) detectors such as DECIGO and BBO will detect ~10^6 neutron-star (NS) binaries and determine the luminosity distances to the binaries with high precision. Combining the luminosity distances with cosmologically-induced phase corrections on the GWs, cosmological expansion out to high redshift can be measured without the redshift determinations of host galaxies by electromagnetic observation and be a unique probe for dark energy. On the other hand, such a NS-binary foreground should be subtracted to detect primordial GWs produced during inflation. Thus, the constraining power on dark energy and the detectability of the primordial gravitational waves strongly depend on the detector sensitivity and are in close relation with one another. In this paper, we investigate the constraints on the equation of state of dark energy with future space-based GW detectors with/without identifying the redshifts of host galaxies. We also study the sensitivity to the primordial GWs, properly dealing with the residual of the NS binary foreground. Based on the results, we discuss the detector sensitivity required to achieve the forementioned targeted study of cosmology.

We report on observations of the Rossiter-McLaughlin (RM) effect for the XO-3 exoplanetary system. The RM effect for the system was previously measured by two different groups, but their results were statistically inconsistent. To obtain a decisive result we observed two full transits of XO-3b with the Subaru 8.2-m telescope. By modeling these data with a new and more accurate analytic formula for the RM effect, we find the projected spin-orbit angle to be \lambda=37.3 deg \pm 3.0 deg, in good agreement with the previous finding by Winn et al. (2009). In addition, an offset of ~22 m/s was observed between the two transit datasets. This offset could be a signal of a third body in the XO-3 system, a possibility that should be checked with future observations. We also attempt to search for a possible signature of the stellar differential rotation in the RM data for the first time, and put weak upper limits on the differential rotation parameters.

We present an improved formula for the anomalous radial velocity of the star during planetary transits due to the Rossiter-McLaughlin (RM) effect. The improvement comes from a more realistic description of the stellar absorption line profiles, taking into account stellar rotation, macroturbulence, thermal broadening, pressure broadening, and instrumental broadening. Although the formula is derived for the case in which radial velocities are measured by cross-correlation, we show through numerical simulations that the formula accurately describes the cases where the radial velocities are measured with the iodine absorption-cell technique. The formula relies on prior knowledge of the parameters describing macroturbulence, instrumental broadening and other broadening mechanisms, but even 30% errors in those parameters do not significantly change the results in typical circumstances. We show that the new analytic formula agrees with previous ones that had been computed on a case-by-case basis via numerical simulations. Finally, as one application of the new formula, we reassess the impact of the differential rotation on the RM velocity anomaly. We show that differential rotation of a rapidly rotating star may have a significant impact on future RM observations.

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.

We measure the halo power spectrum in redshift space from cosmological N-body simulations, and test the analytical models of redshift distortions particularly focusing on the scales of baryon acoustic oscillations (BAOs). Remarkably, the measured halo power spectrum in redshift space exhibits a large-scale enhancement in amplitude relative to the real-space clustering, and the effect becomes significant for the massive or highly biased halo samples. These findings cannot be simply explained by the so-called streaming model frequently used in the literature. By contrast, a physically-motivated perturbation theory model developed in the previous paper reproduces the halo power spectrum very well, and the model combining a simple linear scale-dependent bias can accurately characterize the clustering anisotropies of halos in two dimensions, i.e., line-of-sight and its perpendicular directions. The results highlight the significance of non-linear coupling between density and velocity fields associated with two competing effects of redshift distortions, i.e., Kaiser and Finger-of-God effects, and a proper account of this effect would be important in accurately characterizing the BAOs in two dimensions.

Proposed space-based gravitational-wave detectors such as BBO and DECIGO can detect ~10^6 neutron-star binaries and determine luminosity distance to the binaries with a high precision. Combining the luminosity distance and electromagnetically-derived redshift, one would be able to probe cosmological expansion out to high redshift. In this paper, we show that the Hubble parameter as a function of redshift can be directly measured with monopole and dipole components of the luminosity distance on the sky. As a result, the measurement accuracies of the Hubble parameter in each redshift bin up to z=1 are 3-14 %, 1.5-8 %, and 0.8-4% for the observation time 1 yr, 3 yr, and 10 yr, respectively.

We present forecast results for constraining the primordial non-Gaussianity from photometric surveys through a large-scale enhancement of the galaxy clustering amplitude. In photometric surveys, the distribution of observed galaxies at high redshifts suffers from the gravitational-lensing magnification, which systematically alters the number density for magnitude-limited galaxy samples. We estimate size of the systematic bias in the best-fit cosmological parameters caused by the magnification effect, particularly focusing on the primordial non-Gaussianity. For upcoming deep and/or wide photometric surveys like HSC, DES and LSST, the best-fit value of the non-Gaussian parameter, fNL, obtained from the galaxy count data is highly biased, and the true values of fNL would typically go outside the 3-sigma error of the biased confidence region, if we ignore the magnification effect in the theoretical template of angular power spectrum. The additional information from cosmic shear data helps not only to improve the constraint, but also to reduce the systematic bias. As a result, the size of systematic bias on fNL would become small enough compared to the expected 1-sigma error for HSC and DES, but it would be still serious for deep surveys with z_m > 1.5, like LSST. Tomographic technique improves the constraint on fNL by a factor of 2-3 compared to the one without tomography, but the systematic bias would increase.

We compare the model power spectrum, computed based on perturbation theory (PT) with the power spectrum of luminous red galaxies (LRG) measured from the SDSSDR7 catalog, assuming a flat, CDM-dominated cosmology. The model includes the effects of massive neutrinos, nonlinear matter clustering and nonlinear, scale-dependent galaxy bias in a self-consistent manner. We first test the accuracy of PT-model by comparing the model predictions with the halo power spectrum in real- and redshift-space measured from simulations without massive neutrinos. We show that the PT-model with bias parameters being properly adjusted can fairly well reproduce the simulation results. As a result the best-fit parameters obtained from the hypothetical parameter fitting recover, within statistical uncertainties, the input cosmological parameters in simulations, including an upper bound on neutrino mass, if the power spectrum information up to k~0.15h/Mpc is used. However, for the redshift-space power spectrum, the best-fit cosmological parameters show a sizable bias from the input values if using the information up to k~0.2h/Mpc, probably due to nonlinear redshift distortion effect. Given these tests, we decided, as a conservative choice, to use the LRG power spectrum up to k=0.1h/Mpc in order to minimize possible unknown nonlinearity effects. In combination with the recent results from Wilkinson Microwave Background Anisotropy Probe (WMAP), we derive a robust upper-bound on the sum of neutrino masses, given as m_nu,tot < 0.81eV (95% C.L.), marginalized over other parameters including nonlinear bias parameters and dark energy equation of state parameter. The neutrino mass limit is improved by a factor of 1.85 compared to the limit from the WMAP5 alone, m_nu,tot < 1.5eV.

Baryon acoustic oscillations (BAOs) imprinted in the galaxy power spectrum can be used as a standard ruler to determine angular diameter distance and Hubble parameter at high redshift galaxies. Combining redshift distortion effect which apparently distorts the galaxy clustering pattern, we can also constrain the growth rate of large-scale structure formation. Usually, future forecast for constraining these parameters from galaxy redshift surveys has been made with a full 2D power spectrum characterized as function of wavenumber $k$ and directional cosine $\mu$ between line-of-sight direction and wave vector, i.e., $P(k,\mu)$. Here, we apply the multipole expansion to the full 2D power spectrum, and discuss how much cosmological information can be extracted from the lower-multipole spectra, taking a proper account of the non-linear effects on gravitational clustering and redshift distortion. The Fisher matrix analysis reveals that compared to the analysis with full 2D spectrum, a partial information from the monopole and quadrupole spectra generally degrades the constraints by a factor of $\sim1.3$ for each parameter. The additional information from the hexadecapole spectrum helps to improve the constraints, which lead to an almost comparable result expected from the full 2D spectrum.

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.

We discuss the synergy of the cosmic shear and CMB lensing experiments to simultaneously constrain the neutrino mass and dark energy properties. Taking fully account of the CMB lensing, cosmic shear, CMB anisotropies, and their cross correlation signals, we clarify a role of each signal, and investigate the extent to which the upcoming observations by a high-angular resolution experiment of CMB and deep galaxy imaging survey can tightly constrain the neutrino mass and dark energy equation-of-state parameters. Including the primary CMB information as a prior cosmological information, the Fisher analysis reveals that the time varying equation-of-state parameters, given by w(a)=w_0+w_a(1-a), can be tightly constrained with the accuracies of 5% for w_0 and 15% for w_a, which are comparable to or even better than those of the stage-III type surveys neglecting the effect of massive neutrinos. In other words, including the neutrino mass in the parameter estimation would not drastically alter the Figure-of-Merit estimates of dark energy parameters from the weak lensing measurements. For the neutrino mass, a clear signal for total neutrino mass with 0.1 eV can be detected with 2-sigma significance. The robustness and sensitivity of these results are checked in detail by allowing the setup of cosmic shear experiment to vary as a function of observation time or exposure time, showing that the improvement of the constraints very weakly depends on the survey parameters, and the results mentioned above are nearly optimal for the dark energy parameters and the neutrino mass.

We study the halo bispectrum from non-Gaussian initial conditions. Based on a set of large $N$-body simulations starting from initial density fields with local type non-Gaussianity, we find that the halo bispectrum exhibits a strong dependence on the shape and scale of Fourier space triangles near squeezed configurations at large scales. The amplitude of the halo bispectrum roughly scales as $f_nl^2$. The resultant scaling on the triangular shape is consistent with that predicted by Jeong & Komatsu based on perturbation theory. We systematically investigate this dependence with varying redshifts and halo mass thresholds. It is shown that the $f_nl$ dependence of the halo bispectrum is stronger for more massive haloes at higher redshifts. This feature can be a useful discriminator of inflation scenarios in future deep and wide galaxy redshift surveys.

We present an improved prescription for matter power spectrum in redshift space taking a proper account of both the non-linear gravitational clustering and redshift distortion, which are of particular importance for accurately modeling baryon acoustic oscillations (BAOs). Contrary to the models of redshift distortion phenomenologically introduced but frequently used in the literature, the new model includes the corrections arising from the non-linear coupling between the density and velocity fields associated with two competitive effects of redshift distortion, i.e., Kaiser and Finger-of-God effects. Based on the improved treatment of perturbation theory for gravitational clustering, we compare our model predictions with monopole and quadrupole power spectra of N-body simulations, and an excellent agreement is achieved over the scales of BAOs. Potential impacts on constraining dark energy and modified gravity from the redshift-space power spectrum are also investigated based on the Fisher-matrix formalism. We find that the existing phenomenological models of redshift distortion produce a systematic error on measurements of the angular diameter distance and Hubble parameter by 1~2%, and the growth rate parameter by ~5%, which would become non-negligible for future galaxy surveys. Correctly modeling redshift distortion is thus essential, and the new prescription of redshift-space power spectrum including the non-linear corrections can be used as an accurate theoretical template for anisotropic BAOs.

Characterizing the surfaces of rocky exoplanets via the scattered light will be an essential challenge to investigate the existence of life on habitable exoplanets. We present a simple reconstruction method for fractional areas of different surface types from photometric variations, or colors, of a second Earth. We create mock light curves for Earth without clouds using empirical data. Then these light curves are fitted to the isotropic scattering model consisting of 4 surface types: ocean, soil, snow and vegetation. In an idealized situation where the photometric errors are only photon shot noise, we are able to reproduce the fractional areas of those components fairly well. We may be even able to detect a signature of vegetation from the distinct feature of photosynthesis on the Earth, known as the red edge. In our reconstruction method, Rayleigh scattering due to the atmosphere has an important effect, and for terrestrial exoplanets with atmosphere similar to our Earth, it is possible to estimate the presence of oceans and an atmosphere simultaneously.

We obtain analytical expressions for the velocity anomaly due to the Rossiter- McLaughlin effect, for the case when the anomalous radial velocity is obtained by cross-correlation with a stellar template spectrum. In the limit of vanishing width of the stellar absorption lines, our result reduces to the formula derived by Ohta et al. (2005), which is based on the first moment of distorted stellar lines. Our new formula contains a term dependent on the stellar linewidth, which becomes important when rotational line broadening is appreciable. We generate mock transit spectra for four existing exoplanetary systems (HD17156, TrES-2, TrES- 4, and HD209458) following the procedure of Winn et al. (2005), and find that the new formula is in better agreement with the velocity anomaly extracted from the mock data. Thus, our result provides a more reliable analytical description of the velocity anomaly due to the Rossiter-McLaughlin effect, and explains the previously observed dependence of the velocity anomaly on the stellar rotation velocity.

We study the non-linear evolution of baryon acoustic oscillations in the matter power spectrum and correlation function from the improved perturbation theory (PT). Based on the framework of renormalized PT, we apply the {\it closure approximation} that truncates the infinite series of loop contributions at one-loop order, and obtain a closed set of integral equations for power spectrum and non-linear propagator. The resultant integral expressions keep important non-perturbative properties which can dramatically improve the prediction of non-linear power spectrum. Employing the Born approximation, we then derive the analytic expressions for non-linear power spectrum and the predictions are made for non-linear evolution of baryon acoustic oscillations in power spectrum and correlation function. A detailed comparison between improved PT results and N-body simulations shows that a percent-level agreement is achieved in a certain range in power spectrum and in a rather wider range in correlation function. Combining a model of non-linear redshift-space distortion, we also evaluate the power spectrum and correlation function in correlation function. In contrast to the results in real space, the agreement between N-body simulations and improved PT predictions tends to be worse, and a more elaborate modeling for redshift-space distortion needs to be developed. Nevertheless, with currently existing model, we find that the prediction of correlation function has a sufficient accuracy compared with the cosmic-variance errors for future galaxy surveys with volume of a few (Gpc/h)^3 at z>=0.5.

Future or ongoing galaxy redshift surveys can put stringent constraints on neutrinos masses via the high-precision measurements of galaxy power spectrum, when combined with cosmic microwave background (CMB) information. In this paper we develop a method to model galaxy power spectrum in the weakly nonlinear regime for a mixed dark matter (CDM plus finite-mass neutrinos) model, based on perturbation theory (PT) whose validity is well tested by simulations for a CDM model. In doing this we carefully study various aspects of the nonlinear clustering and then arrive at a useful approximation allowing for a quick computation of the nonlinear power spectrum as in the CDM case. The nonlinear galaxy bias is also included in a self-consistent manner within the PT framework. Thus the use of our PT model can give a more robust understanding of the measured galaxy power spectrum as well as allow for higher sensitivity to neutrino masses due to the gain of Fourier modes beyond the linear regime. Based on the Fisher matrix formalism, we find that BOSS or Stage-III type survey, when combined with Planck CMB information, gives a precision of total neutrino mass constraint, sigma(m_nu,tot) 0.1eV, while Stage-IV type survey may achieve sigma(m_nu,tot) 0.05eV, i.e. more than a 1-sigma detection of neutrino masses. We also discuss possible systematic errors on dark energy parameters caused by the neutrino mass uncertainty. The significant correlation between neutrino mass and dark energy parameters is found, if the information on power spectrum amplitude is included. More importantly, for Stage-IV type survey, a best-fit dark energy model may be biased and falsely away from the underlying true model by more than the 1-sigma statistical errors, if neutrino mass is ignored in the model fitting.

We present a formalism to calculate the non-linear matter power spectrum in modified gravity models that explain the late-time acceleration of the Universe without dark energy. Any successful modified gravity models should contain a mechanism to recover General Relativity (GR) on small scales in order to avoid the stringent constrains on deviations from GR at solar system scales. Based on our formalism, the quasi non-linear power spectrum in the Dvali-Gabadadze-Porratti (DGP) braneworld models and $f(R)$ gravity models are derived by taking into account the mechanism to recover GR properly. We also extrapolate our predictions to fully non-linear scales using the Parametrized Post Friedmann (PPF) framework. In $f(R)$ gravity models, the predicted non-linear power spectrum is shown to reproduce N-body results. We find that the mechanism to recover GR suppresses the difference between the modified gravity models and dark energy models with the same expansion history, but the difference remains large at weakly non-linear regime in these models. Our formalism is applicable to a wide variety of modified gravity models and it is ready to use once consistent models for modified gravity are developed.

We present a new numerical scheme to treat the non-linear evolution of cosmological power spectra. Governing equations for matter power spectra have been previously derived by a non-perturbative technique with closure approximation. Solutions of the resultant closure equations just correspond to the resummation of an infinite class of perturbation corrections, and they consistently reproduce the one-loop results of standard perturbation theory. We develop a numerical algorithm to solve closure evolutions in both perturbative and non-perturbative regimes. The present numerical scheme is particularly suited for examining non-linear matter power spectrum in general cosmological models, including modified theory of gravity. As a demonstration, we study weakly non-linear evolution of power spectrum in a class of modified gravity models, as well as various dark energy models.

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.

In a general metric theory of gravitation in four dimensions, six polarizations of a gravitational wave are allowed: two scalar and two vector modes, in addition to two tensor modes in general relativity. Such additional polarization modes appear due to additional degrees of freedom in modified theories of gravitation or theories with extra dimensions. Thus, observations of gravitational waves can be utilized to constrain the extended models of gravitation. In this paper, we investigate detectability of additional polarization modes of gravitational waves, particularly focusing on a stochastic gravitational-wave background, with laser-interferometric detectors on the Earth. We found that multiple detectors can separate the mixture of polarization modes in detector outputs, and that they have almost the same sensitivity to each polarization mode of stochastic gravitational-wave background.

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.

We present theoretical predictions for photometric and spectroscopic signatures of rings around transiting extrasolar planets. On the basis of a general formulation for the transiting signature in the stellar light curve and the velocity anomaly due to the Rossiter effect, we compute the expected signals analytically for a face-on ring system, and numerically for more general configurations. We study the detectability of a ring around a transiting planet located at a=3AU for a variety of obliquity and azimuthal angles, and find that it is possible to detect the ring signature both photometrically and spectroscopically unless the ring is almost edge-on (i.e., the obliquity angle of the ring \theta is much less than unity). We also consider the detectability of planetary rings around a close-in planet, HD 209458b ($\theta\approx 90^\circ-i_{\rm orbit} \approx 3^\circ.32$), and Saturn ($\theta \approx 26^\circ.7$) as illustrative examples. While the former is difficult to detect with the current precision (photometric precision of 10^-4 and radial velocity precision of 1 m/s), a marginal detection of the latter is possible photometrically. If the future precision of the radial velocity measurement reaches even below 0.1 m/s, they will be even detectable from the ground-based spectroscopic observations.

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 present new spectroscopic and photometric observations of the HAT-P-1 planetary system. Spectra obtained during three transits exhibit the Rossiter-McLaughlin effect, allowing us to measure the angle between the sky projections of the stellar spin axis and orbit normal, \lambda = 3.7 +/- 2.1 degrees. The small value of \lambda for this and other systems suggests that the dominant planet migration mechanism preserves spin-orbit alignment. Using two new transit light curves, we refine the transit ephemeris and reduce the uncertainty in the orbital period by an order of magnitude. We find a upper limit on the orbital eccentricity of 0.067, with 99% confidence, by combining our new radial-velocity measurements with those obtained previously.

We present the first attempt to analytically study the nonlinear matter power spectrum for a mixed dark matter (cold dark matter plus neutrinos of total mass ~0.1eV) model based on cosmological perturbation theory. The suppression in the power spectrum amplitudes due to massive neutrinos is, compared to the linear regime, enhanced in the weakly nonlinear regime where standard linear theory ceases to be accurate. We demonstrate that, thanks to this enhanced effect and the gain in the range of wavenumbers to which the PT prediction is applicable, the use of such a nonlinear model may enable a precision of sigma(m_nu,tot) ~ 0.07eV in constraining the total neutrino mass for the planned galaxy redshift survey, a factor of 2 improvement compared to the linear regime.

The Stokes V parameter characterizes asymmetry of amplitudes between right- and left-handed waves, and non-vanishing value of the V parameter yields a circularly polarized signal. Cosmologically, V parameter may be a direct probe for parity violation in the universe. In this paper, we theoretically investigate a measurement of this parameter, particularly focusing on the gravitational-wave backgrounds observed via ground-based interferometers. In contrast to the traditional analysis that only considers the total amplitude (or equivalently $\Omega_{GW}$), the signal analysis including a circular-polarized mode has a rich structure due to the multi-dimensionality of target parameters. We show that, by using the network of next-generation detectors, separation between polarized and unpolarized modes can be performed with small statistical loss induced by their correlation.

We apply a non-linear statistical method in turbulence to the cosmological perturbation theory and derive a closed set of evolution equations for matter power spectra. The resultant closure equations consistently recover the one-loop results of standard perturbation theory and beyond that, it is still capable of treating the non-linear evolution of matter power spectra. We find the exact integral expressions for the solutions of closure equations. These analytic expressions coincide with the renormalized one-loop results presented by Crocce & Scoccimarro (2006,2007). By constructing the non-linear propagator, we analytically evaluate the non-linear matter power spectra based on the first-order Born approximation of the integral expressions and compare it with those of the renormalized perturbation theory.

We report a measurement of the Rossiter--McLaughlin effect in the transiting extrasolar planetary system TrES-1, via simultaneous spectroscopic and photometric observations with the Subaru and MAGNUM telescopes. By modeling the radial velocity anomaly that was observed during a transit, we determine the sky-projected angle between the stellar spin axis and the planetary orbital axis to be $\lambda = 30 \pm 21$ [deg]. This is the third case for which $\lambda$ has been measured in a transiting exoplanetary system, and the first demonstration that such measurements are possible for relatively faint host stars ($V \sim 12$, as compared to $V \sim 8$ for the other systems). We also derive a time of mid-transit, constraints on the eccentricity of the TrES-1b orbit ($e = 0.048 \pm 0.025$), and upper limits on the mass of the Trojan companions ($\lesssim$14 $M_{\oplus}$) at the 3$\sigma$ level.

An acoustic oscillation of the primeval photon-baryon fluid around the decoupling time imprints a characteristic scale in the galaxy distribution today, known as the baryon acoustic oscillation (BAO) scale. Several on-going and/or future galaxy surveys aim at detecting and precisely determining the BAO scale so as to trace the expansion history of the universe. We consider nonlinear and redshift-space distortion effects on the shifts of the BAO scale in $k$-space using perturbation theory. The resulting shifts are indeed sensitive to different choices of the definition of the BAO scale, which needs to be kept in mind in the data analysis. We present a toy model to explain the physical behavior of the shifts. We find that the BAO scale defined as in Percival et al. (2007) indeed shows very small shifts ($\lesssim$ 1%) relative to the prediction in {\it linear theory} in real space. The shifts can be predicted accurately for scales where the perturbation theory is reliable.

We discuss the polarization signature of primordial gravitational waves imprinted in cosmic microwave background (CMB) anisotropies. The high-energy physics motivated by superstring theory or M-theory generically yield parity violating terms, which may produce a circularly polarized gravitational wave background (GWB) during inflation. In contrast to the standard prediction of inflation with un-polarized GWB, circularly polarized GWB generates non-vanishing TB and EB-mode power spectra of CMB anisotropies. We evaluate the TB and EB-mode power spectra taking into account the secondary effects and investigate the dependence of cosmological parameters. We then discuss current constraints on the circularly polarized GWB from large angular scales (l < 16) of the three year WMAP data. Prospects for future CMB experiments are also investigated based on a Monte Carlo analysis of parameter estimation, showing that the circular polarization degree, varepsilon, which is the asymmetry of the tensor power spectra between right- and left-handed modes normalized by the total amplitude, can be measured down to |varepsilon| 0.35(r/0.05)^{-0.6}.

We show that pairs of widely separated interferometers are advantageous for measuring the Stokes parameter V of a stochastic background of gravitational waves. This parameter characterizes asymmetry of amplitudes of right- and left-handed waves and generation of the asymmetry is closely related to parity violation in the early universe. The advantageous pairs include LIGO(Livingston)-LCGT and AIGO-Virgo that are relatively insensitive to Omega_GW (the simple intensity of the background). Using at least three detectors, information of the intensity Omega_GW and the degree of asymmetry V can be separately measured.

We discuss a robust data analysis method to detect a stochastic background of gravitational waves in the presence of non-Gaussian noise. In contrast to the standard cross-correlation (SCC) statistic frequently used in the stochastic background searches, we consider a {\it generalized cross-correlation} (GCC) statistic, which is nearly optimal even in the presence of non-Gaussian noise. The detection efficiency of the GCC statistic is investigated analytically, particularly focusing on the statistical relation between the false-alarm and the false-dismissal probabilities, and the minimum detectable amplitude of gravitational-wave signals. We derive simple analytic formulae for these statistical quantities. The robustness of the GCC statistic is clarified based on these formulae, and one finds that the detection efficiency of the GCC statistic roughly corresponds to the one of the SCC statistic neglecting the contribution of non-Gaussian tails. This remarkable property is checked by performing the Monte Carlo simulations and successful agreement between analytic and simulation results was found.

We consider nonlinear biasing models of galaxies with particular attention to a correlation between linear and quadratic biasing coefficients, b_1 and b_2. We first derive perturbative expressions for b_1 and b_2 in halo and peak biasing models. Then we compute power spectra and bispectra of dark matter particles and halos using N-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their b_1 and b_2. We find that the values of those coefficients at linear regimes (k<0.2h/Mpc) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, Q, for equilateral triangles, are insensitive to the values of b_1 implying that b_2 indeed correlates with b_1. The present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between b_1 and b_2 are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

We develop a numerical scheme to make a high-frequency skymap of gravitational-wave backgrounds (GWBs) observed via space-based interferometer. Based on the cross-correlation technique, the intensity distribution of anisotropic GWB can be directly reconstructed from the time-ordered data of cross-correlation signals, with full knowledge of detector's antenna pattern functions. We demonstrate how the planned space interferometer, LISA, can make a skymap of GWB for a specific example of anisotropic signals. At the frequency higher than the characteristic frequency $f_*=1/(2\pi L)$, where $L$ is the arm-length of the detector, the reconstructed skymap free from the instrumental noise potentially reaches the angular resolution up to the multipoles $\ell\sim10$. The presence of instrumental noises degrades the angular resolution. The resultant skymap has angular resolution with multipoles $\ell\leq 6\sim7$ for the anisotropic signals with signal-to-noise ratio S/N$>5$.

Future missions of gravitational-wave astronomy will be operated by space-based interferometers, covering very wide range of frequency. Search for stochastic gravitational-wave backgrounds (GWBs) is one of the main targets for such missions, and we here discuss the prospects for direct measurement of isotropic and anisotropic components of (primordial) GWBs around the frequency 0.1-10 Hz. After extending the theoretical basis for correlation analysis, we evaluate the sensitivity and the signal-to-noise ratio for the proposed future space interferometer missions, like Big-Bang Observer (BBO), Deci-Hertz Interferometer Gravitational-wave Observer (DECIGO) and recently proposed Fabry-Perot type DECIGO. The astrophysical foregrounds which are expected at low frequency may be a big obstacle and significantly reduce the signal-to-noise ratio of GWBs. As a result, minimum detectable amplitude may reach h^2 \ogw = 10^{-15} \sim 10^{-16}, as long as foreground point sources are properly subtracted. Based on correlation analysis, we also discuss measurement of anisotropies of GWBs. As an example, the sensitivity level required for detecting the dipole moment of GWB induced by the proper motion of our local system is closely examined.

We present a perturbative reconstruction method to make a skymap of gravitational-wave backgrounds (GWBs) observed via space-based interferometer. In the presence of anisotropies in GWBs, the cross-correlated signals of observed GWBs are inherently time-dependent due to the non-stationarity of the gravitational-wave detector. Since the cross-correlated signal is obtained through an all-sky integral of primary signals convolving with the antenna pattern function of gravitational-wave detectors, the non-stationarity of cross-correlated signals, together with full knowledge of antenna pattern functions, can be used to reconstruct an intensity map of the GWBs. Here, we give two simple methods to reconstruct a skymap of GWBs based on the perturbative expansion in low-frequency regime. The first one is based on harmonic-Fourier representation of data streams and the second is based on "direct" time-series data. The latter method enables us to create a skymap in a direct manner. The reconstruction technique is demonstrated in the case of the Galactic gravitational wave background observed via planned space interferometer, LISA. Although the angular resolution of low-frequency skymap is rather restricted, the methodology presented here would be helpful in discriminating the GWBs of galactic origins by those of the extragalactic and/or cosmological origins.

We discuss the gravitational wave background (GWB) from a cosmological population of gamma-ray bursts (GRBs). Among various emission mechanisms for the gravitational waves (GWs), we pay a particular attention to the vast anisotropic neutrino emissions from the accretion disk around the black hole formed after the so-called failed supernova explosions. The produced GWs by such mechanism are known as burst with memory, which could dominate over the low-frequency regime below \sim 10Hz. To estimate their amplitudes, we derive general analytic formulae for gravitational waveform from the axisymmetric jets. Based on the formulae, we first quantify the spectrum of GWs from a single GRB. Then, summing up its cosmological population, we find that the resultant value of the density parameter becomes roughly \Omega_{GW} \approx 10^{-20} over the wide-band of the low-frequency region, f\sim 10^{-4}-10^1Hz. The amplitude of GWB is sufficiently smaller than the primordial GWBs originated from an inflationary epoch and far below the detection limit.

In this paper, a quantitative characterization for the evolutionary sequence of stellar self-gravitating system is investigated, focusing on the pre-collapse stage of the long-term dynamical evolution. In particular, we consider the quasi-equilibrium behaviors of the N-body systems in the setup of the so-called Antonov problem, i.e., self-gravitating N-body system confined in an adiabatic wall and try to seek a possible connection with thermostatistics of self-gravitating systems. For this purpose, a series of long-term N-body simulations with various initial conditions are performed. We found that a quasi-equilibrium sequence away from the thermal equilibrium can be characterized by the one-parameter family of the stellar models. Especially, the stellar polytropic distribution satisfying the effective equation of state $P\propto\rho^{1+1/n}$ provides an excellent approximation to the evolutionary sequence of the N-body system. Based on the numerical results, we discuss a link between the quasi-equilibrium state and the generalized thermostatistics by means of the non-extensive entropy.

We determine the stellar, planetary, and orbital properties of the transiting planetary system HD 209458, through a joint analysis of high-precision radial velocities, photometry, and timing of the secondary eclipse. Of primary interest is the strong detection of the Rossiter-McLaughlin effect, the alteration of photospheric line profiles that occurs because the planet occults part of the rotating surface of the star. We develop a new technique for modeling this effect, and use it to determine the inclination of the planetary orbit relative to the apparent stellar equator [lambda = (-4.4 +/- 1.4) degrees], and the line-of-sight rotation speed of the star [v*sin(I) = (4.70 +/- 0.16) km/s]. The uncertainty in these quantities has been reduced by an order of magnitude relative to the pioneering measurements by Queloz and collaborators. The small but nonzero misalignment is probably a relic of the planet formation epoch, because the expected timescale for tidal coplanarization is larger than the age of the star. Our determination of v*sin(I) is a rare case in which rotational line broadening has been isolated from other broadening mechanisms.

We present a first morphometric investigation of a preliminary sample from the SDSS of 154287 galaxies with apparent magnitude 14.5<m_r<17.5 and redshift 0.001<z<0.4. We measure the Minkowski Functionals, which are a complete set of morphological descriptors. To account for the complicated wedge--like geometry of the present survey data, we construct isodensity contour surfaces from the galaxy positions in redshift space and employ two complementary methods of computing the Minkowski Functionals. We find that the observed Minkowski Functionals for SDSS galaxies are consistent with the prediction of a Lambda--dominated spatially--flat Cold Dark Matter model with random--Gaussian initial conditions, within the cosmic variance estimated from the corresponding mock catalogue. We expect that future releases of the SDSS survey will allow us to distinguish morphological differences in the galaxy distribution with regard to different morphological type and luminosity ranges.

A transiting extrasolar planet sequentially blocks off the light coming from the different parts of the disk of the host star in a time dependent manner. Due to the spin of the star, this produces an asymmetric distortion in the line profiles of the stellar spectrum, leading to an apparent anomaly of the radial velocity curves, known as the Rossiter - McLaughlin effect. Here, we derive approximate but accurate analytic formulae for the anomaly of radial velocity curves taking account of the stellar limb darkening. The formulae are particularly useful in extracting information of the projected angle between the planetary orbit axis and the stellar spin axis, \lambda, and the projected stellar spin velocity, V sin I_s. We create mock samples for the radial curves for the transiting extrasolar system HD209458, and demonstrate that constraints on the spin parameters (V sin I_s, \lambda) may be significantly improved by combining our analytic template formulae and the precision velocity curves from high-resolution spectroscopic observations with 8-10 m class telescopes. Thus future observational exploration of transiting systems using the Rossiter - McLaughlin effect is one of the most important probes to better understanding of the origin of extrasolar planetary systems, especially the origin of their angular momentum.

We discuss the sensitivity to anisotropies of stochastic gravitational-wave backgrounds (GWBs) observed via space-based interferometer. In addition to the unresolved galactic binaries as the most promising GWB source of the planned Laser Interferometer Space Antenna (LISA), the extragalactic sources for GWBs might be detected in the future space missions. The anisotropies of the GWBs thus play a crucial role to discriminate various components of the GWBs. We study general features of antenna pattern sensitivity to the anisotropies of GWBs beyond the low-frequency approximation. We show that the sensitivity of space-based interferometer to GWBs is severely restricted by the data combinations and the symmetries of the detector configuration. The spherical harmonic analysis of the antenna pattern functions reveals that the angular power of the detector response increases with frequency and the detectable multipole moments with effective sensitivity h_{eff} \sim 10^{-20} Hz^{-1/2} may reach $\ell \sim$ 8-10 at $f \sim f_*=10$ mHz in the case of the single LISA detector. However, the cross correlation of optimal interferometric variables is blind to the monopole (\ell=0) intensity anisotropy, and also to the dipole (\ell=1) in some case, irrespective of the frequency band. Besides, all the self-correlated signals are shown to be blind to the odd multipole moments (\ell=odd), independently of the frequency band.

We discuss the cosmological evolution of gravitational waves (GWs) after inflation in a brane-world cosmology embedded in five-dimensional anti-de Sitter (AdS_5) bulk spacetime. In a brane-world scenario, the evolution of GWs is affected by the non-standard cosmological expansion and the excitation of the Kaluza-Klein modes (KK-modes), which are significant in the high-energy regime of the universe. We numerically solve the wave equation of GWs in the Poincare coordinates of the AdS_5 spacetime. Using a plausible initial condition from inflation, we find that, while the behavior of GWs in the bulk is sensitive to the transition time from inflation to the radiation dominated epoch, the amplitude of GWs on the brane is insensitive to this time if the transition occurs early enough before horizon re-entry. As a result, the amplitude of GWs is suppressed by the excitation of KK-modes which escape from the brane into the bulk, and the effect may compensate the enhancement of the GWs by the non-standard cosmological expansion. Based on this, the influence of the high-energy effects on the stochastic background of GWs is discussed.

After introducing the fundamental properties of self-gravitating systems, we present an application of Tsallis' generalized entropy to the analysis of their thermodynamic nature. By extremizing the Tsallis entropy, we obtain an equation of state known as the stellar polytrope. For a self-gravitating stellar system confined within a perfectly reflecting wall, we discuss the thermodynamic instability caused by its negative specific heat. The role of the extremum as a quasi-equilibrium is also demonstrated from the results of N-body simulations.

We calculate the one-point probability distribution function (PDF) for cosmic density in non-linear regime of the gravitational evolution. Under the local approximation that the evolution of cosmic fluid fields can be characterized by the Lagrangian local dynamics with finite degrees of freedom, the analytic expressions of PDF are derived taking account of the smoothing effect. The validity and the usefulness of the local approximation are then discussed comparing those results with N-body simulations in a Gaussian initial condition. Adopting the ellipsoidal collapse model (ECM) and the spherical collapse model (SCM) as Lagrangian local dynamics, we found that the PDFs from the local approximation excellently match the simulation results in the case of the cold dark matter initial spectrum. As for the scale-free initial spectra given by $P(k)\propto k^n$, N-body result suffers from spurious numerical effects, which prevent us to give a detailed comparison. Nevertheless, at the quality of N-body data, the model predictions based on the ECM and the SCM quantitatively agree with N-body results in cases with spectral index n<0. For the index n>=0, choice of the Lagrangian local dynamics becomes crucial for an accurate prediction and a more delicate modeling is required, however, we find that the model prediction based on the ECM provides a better approximation to the N-body results of cumulants and PDFs.

We study the evolution of gravitational waves(GWs) after inflation in a brane-world cosmology embedded in five-dimensional anti-de Sitter spacetime. Contrary to the standard four-dimensional results, the GWs at the high-energy regime in brane-world model suffer from the effects of the non-standard cosmological expansion and the excitation of the Kaluza-Klein modes(KK-modes), which can affect the amplitude of stochastic gravitational wave background significantly. To investigate these two high-energy effects quantitatively, we numerically solve the wave equation of the GWs in the radiation dominated epoch at relatively low-energy scales. We show that the resultant GWs is suppressed by the excitation of the KK modes. The created KK modes are rather soft and escape away from the brane to the bulk gravitational field. The results are also compared to the semi-analytic prediction from the low-energy approximation and the evolved amplitude of GWs on the brane reasonably matches the numerical simulations.

With particular attention to the recently postulated introduction of a non-extensive generalization of Boltzmann-Gibbs statistics, we study the long-term stellar dynamical evolution of self-gravitating systems on timescales much longer than the two-body relaxation time. In a self-gravitating N-body system confined in an adiabatic wall, we show that the quasi-equilibrium sequence arising from the Tsallis entropy, so-called stellar polytropes, plays an important role in characterizing the transient states away from the Boltzmann-Gibbs equilibrium state.

We present a numerical analysis of genus statistics for dark matter halo catalogs from the Hubble volume simulation. The huge box-size of the Hubble volume simulation enables us to carry out a reliable statistical analysis of the biasing properties of halos at a Gaussian smoothing scale of R_G>30Mpc/h with a cluster-mass scale of between 7*10^{13}Msolar/h and 6*10^{15}Msolar/h. A detailed comparison of the genus for dark matter halos with that for the mass distribution shows that the non-Gaussianity induced by the halo biasing is comparable to that by nonlinear gravitational evolution, and both the shape and the amplitude of the genus are almost insensitive to the halo mass at R_G>30Mpc/h. In order to characterize the biasing effect on the genus, we apply a perturbative formula developed by Matsubara (1994). We find that the perturbative formula well describes the simulated halo genus at R_G>50Mpc/h. The result indicates that the biasing effect on the halo genus is well approximated by nonlinear deterministic biasing up to the second-order term in the mass density fluctuation. The two parameters describing the linear and quadratic terms in the nonlinear biasing accurately specify the genus for galaxy clusters.

We present a general framework to treat the evolution of one-point probability distribution function (PDF) for cosmic density $\delta$ and velocity-divergence fields $\theta$. In particular, we derive an evolution equation for the one-point PDFs and consider the stochastic nature associated with these quantities. Under the local approximation that the evolution of cosmic fluid fields can be characterized by the Lagrangian local dynamics with finite degrees of freedom, evolution equation for PDFs becomes a closed form and consistent formal solutions are constructed. Adopting this local approximation, we explicitly evaluate the one-point PDFs $P(\delta)$ and $P(\theta)$ from the spherical and the ellipsoidal collapse models as the representative Lagrangian local dynamics. In a Gaussian initial condition, while the local density PDF from the ellipsoidal model almost coincides with the that of the spherical model, differences between spherical and ellipsoidal collapse model are found in the velocity-divergence PDF. Importantly, the joint PDF of local density, $P(\delta,t;\delta',t')$, evaluated at the same Lagrangian position but at the different times $t$ and $t'$ from the ellipsoidal collapse model exhibits a large amount of scatter. The mean relation between $\delta$ and $\delta'$ does fail to match the one-to-one mapping obtained from spherical collapse model. Moreover, the joint PDF $P(\delta;\theta)$ from the ellipsoidal collapse model shows a similar stochastic feature, both of which are indeed consistent with the recent result from N-body simulations.

We revisit the issues on the thermodynamic property of stellar self-gravitating system arising from Tsallis' non-extensive entropy. Previous papers (Taruya & Sakagami, Physica A 307 (2002) 185 (cond-mat/0107494); ibid. (2002) in press (cond-mat/0204315)) have revealed that the extremum-state of Tsallis entropy characterized by the so-called stellar polytrope has consistent thermodynamic structure, which predicts the thermodynamic instability due to the negative specific heat. However, their analyses heavily relies on the old Tsallis formalism using standard linear mean values. In this paper, extending our previous study, we focus on the equilibrium structure based on the standard framework by means of the normalized q-expectation values. It then turns out that the new extremum-state of Tsallis entropy essentially remains unchanged from the previous result, i.e., the stellar quasi-equilibrium distribution can be described by the stellar polytrope. While the thermodynamic stability for a system confined in an adiabatic wall completely agrees with the previous study and thereby the stability/instability criterion remains unchanged, the stability analysis reveals a new equilibrium property for the system surrounded by a thermal bath. In any case, the stability/instability criteria are consistently explained from the presence of negative specific heat and within the formalism, the stellar polytrope is characterized as a plausible non-extensive meta-equilibrium state.

We present a simple model based on the dark halo approach which provides a useful way to understand key points determining the shape of the non-Gaussian tails of the dark matter one-point probability distribution function(PDF). In particular, using the scale-free models with power-law profile of dark halos, we derive a simple analytic expression for the one-point PDF. It is found that the shape of the PDF changes at the characteristic value of $\delta_*$ which is defined by the smoothed density of a halo with the characteristic mass $M_*$ at the epoch. In cold dark matter models with top-hat smoothing filters, the characteristic smoothed density at present time typically takes the value $\delta_*\gg 1$ for a small smoothing scale $\rth\sim 1$Mpc$/h$ and conversely $\delta_*\ll 1$ for a large smoothing scale $\rth > 10$Mpc$/h$. On the range $\delta/\delta_*<1$, the shape of the PDF is almost solely determined by the outer slope of halos and scales as a power-law. The resultant non-Gaussian tails of PDF then resemble the log-normal PDFs in that range and show a good agreement with N-body simulations, which can be ascribed to the universality of the outer slope of the halo profile. In contrast, tails of one-point PDF in the range $\delta/\delta_*>1$ basically follow the steep exponential tails of the halo mass function, which exhibit a strong sensitivity to both the outer slope of the halo profile and the initial power spectrum. Based on these results, the discussion on the PDF of galaxy distribution and the application to the weak lensing statistics are also presented.

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 investigate the thermodynamic properties of stellar self-gravitating system arising from the Tsallis generalized entropy. In particular, physical interpretation of the thermodynamic instability, as has been revealed by previous paper(Taruya & Sakagami, cond-mat/0107494, Physica A 307, 185 (2002)), is discussed in detail based on the non-extensive thermostatistics. Examining the Clausius relation in a quasi-static experiment, we obtain the standard result of thermodynamic relation that the physical temperature of the equilibrium non-extensive system is identified with the inverse of the Lagrange multiplier, $T_{phys}=1/\beta$. Using this relation, the specific heat of total system is computed, and confirm the common feature of self-gravitating system that the presence of negative specific heat leads to the thermodynamic instability. In addition to the gravothermal instability discovered previously, the specific heat shows the curious divergent behavior at the polytrope index $n>3$, suggesting another type of thermodynamic instability. Evaluating the second variation of free energy, we find that the marginal stability condition indicated from the specific heat can be exactly recovered from the second variation of free energy. Thus, the stellar polytropic system is consistently characterized by the non-extensive thermostatistics as a plausible thermal equilibrium state. We also clarify the non-trivial scaling behavior appeared in specific heat and address the origin of non-extensive nature in stellar polytrope.

We derive the pairwise peculiar velocity distribution function of dark matter particles applying the dark matter halo approach. Unlike the previous work, we do not assume a Gaussian velocity distribution function of dark matter in a single halo, but compute it self-consistently with the assumed density profile for dark matter halo. The resulting distribution function is well approximated by an exponential distribution which is consistent with the previous observational, numerical and theoretical results. We also compute the pairwise peculiar velocity dispersion for different density profiles, and provide a practical fitting formula. We apply an empirical biasing scheme into our model and present prediction for pairwise peculiar velocity dispersion of {\it galaxies}, and reproduce the previous results of simulations using our semi-analytical method.

The statistical property of the weak lensing fields is studied quantitatively using the ray-tracing simulations. Motivated by the empirical lognormal model that characterizes the probability distribution function(PDF) of the three-dimensional mass distribution excellently, we critically investigate the validity of lognormal model in the weak lensing statistics. Assuming that the convergence field, $\kappa$, is approximately described by the lognormal distribution, we present analytic formulae of convergence for the one-point PDF and the Minkowski functionals. Comparing those predictions with ray-tracing simulations in various cold dark matter models, we find that the one-point lognormal PDF can describe the non-Gaussian tails of convergence fields accurately up to $\nu\sim10$, where $\nu$ is the level threshold given by $\nu\equiv\kappa/\var^{1/2}$, although the systematic deviation from lognormal prediction becomes manifest at higher source redshift and larger smoothing scales. The lognormal formulae for Minkowski functionals also fit to the simulation results when the source redshift is low. Accuracy of the lognormal-fit remains good even at the small angular scales, where the perturbation formulae by Edgeworth expansion break down. On the other hand, lognormal models does not provide an accurate prediction for the statistics sensitive to the rare events such as the skewness and the kurtosis of convergence. We therefore conclude that the empirical lognormal model of the convergence field is safely applicable as a useful cosmological tool, as long as we are concerned with the non-Gaussianity of $\nu\simlt5$ for low source redshift samples.

The distribution of differential time delays \Delta t between images produced by strong gravitational lensing contains information on the mass distributions in the lensing objects as well as on cosmological parameters such as H_0. We derive an explicit expression for the conditional probability distribution function of time delays P(\Delta t|\theta), given an image separation between multiple images \theta, and related statistics. We consider lensing halos described by the singular isothermal sphere (SIS) approximation and by its generalization as proposed by Navarro, Frenk, & White (NFW) which has a density profile \rho \propto r^{-\alpha} in the innermost region. The time delay distribution is very sensitive to these profiles; steeper inner slopes tend to produce larger time delays. For example, if H_0=70 km/s/Mpc, a \Lambda-dominated cosmology and a source redshift z_S=1.27 are assumed, lenses with \theta=5'' produce a time delay of \Delta t[yr]=1.5^{+1.7}_{-0.9}, 0.39^{+0.37}_{-0.22}, 0.15^{+0.11}_{-0.09}, and 0.071^{+0.054}_{-0.038} (50% confidence interval), for SIS, generalized NFW with \alpha=1.5, \alpha=1.0, and \alpha=0.5, respectively. At a fixed image separation, the time delay is determined by the difference in the lensing potential between the position of the two images, which typically occur at different impact parameters. Although the values of \Delta t are proportional to the inverse of H_0, P(\Delta t|\theta) is rather insensitive to all other cosmological model parameters, source redshifts, magnification biases and so on. A knowledge of P(\Delta t|\theta) will also be useful in designing the observing program of future large scale synoptic variability surveys and for evaluating possible selection biases operating against large splitting lens systems.

We present a first physical application of Tsallis' generalized entropy to the thermodynamics of self-gravitating systems. The stellar system confined in a spherical cavity of radius $r_e$ exhibits an instability, so-called gravothermal catastrophe, which has been originally investigated by Antonov (1962) and Lynden-Bell & Wood (1968) on the basis of the maximum entropy principle for the phase-space distribution function. In contrast to previous analyses using the Boltzmann-Gibbs entropy, we apply the Tsallis-type generalized entropy to seek the equilibrium criteria. Then the distribution function of Vlassov-Poisson system can be reduced to the stellar polytrope system. Evaluating the second variation of Tsallis entropy and solving the zero eigenvalue problem explicitly, we find that the gravothermal instability appears in cases with polytrope index $n>5$. The critical point characterizing the onset of instability are obtained, which exactly matches with the results derived from the standard turning-point analysis. The results give an important suggestion that the Tsallis entropy is indeed applicable and viable to the long-range nature of the self-gravitating system.

The well-known self-similar solution for the two-point correlation function of the density field is valid only in an Einstein-de Sitter universe. We attempt to extend the solution for non-Einstein-de Sitter universes. For this purpose we introduce an idea of quasi-self-similar evolution; this approach is based on the assumption that the evolution of the two-point correlation is a succession of stages of evolution, each of which spans a short enough period to be considered approximately self-similar. In addition we assume that clustering is stable on scales where a physically motivated `virialization condition' is satisfied. These assumptions lead to a definite prediction for the behavior of the two-point correlation function in the strongly nonlinear regime. We show that the prediction agrees well with N-body simulations in non-Einstein-de Sitter cases, and discuss some remaining problems.

We quantitatively study the probability distribution function (PDF) of cosmological nonlinear density fluctuations from N-body simulations with Gaussian initial condition. In particular, we examine the validity and limitations of one-point and two-point log-normal PDF models against those directly estimated from the simulations. We find that the one-point log-normal PDF describes very accurately the cosmological density distribution even in the nonlinear regime (the rms variance \sigma_{nl} \simlt 4 and the over-density \delta \simlt 100). Furthermore the two-point log-normal PDFs are also in good agreement with the simulation data from linear to fairly nonlinear regime, while slightly deviate from them for \delta \simlt -0.5. Thus the log-normal PDF can be used as a useful empirical model for the cosmological density fluctuations. While this conclusion is fairly insensitive to the shape of the underlying power spectrum of density fluctuations P(k), models with substantial power on large scales, i.e., n\equiv d\ln P(k)/d \ln k \simlt -1, are better described by the log-normal PDF. On the other hand, we note that the one-to-one mapping of the initial and the evolved density fields consistent with the log-normal model does not approximate the broad distribution of their mutual correlation even on average. Thus the origin of the phenomenological log-normal PDF approximation still remains to be understood.

We study the arc statistics of gravitational lensing generated by dark matter halos in order to probe their density profile. We characterize the halo profile by two parameters, the inner slope of the central cusp $\alpha$, and the median amplitude of the concentration parameter, $c_{\rm norm}$, for a halo of mass $10^{14}h^{-1}M_\odot$ at $z=0$, and compute the numbers of tangential and radial arcs produced by gravitational lensing of galaxy clusters. We find that the number of arcs divided by the number of halos is a good statistic which is sensitive to both $c_{\rm norm}$ and $\alpha$ with very weak dependence on the cosmological parameters. If the arc samples with well-defined selection criteria for the clusters become available, one can strongly constrain both $c_{\rm norm}$ and $\alpha$. While our tentative comparison with the existing observational data indicates that the inner density profile of dark halos is indeed as steep as predicted by recent simulations ($\alpha \sim 1.5$), the homogeneous samples of tangential and radial arcs are required for more quantitative discussions.

We perform an extensive analysis of nonlinear and stochastic biasing of galaxies and dark halos in spatially flat low-density CDM universe using cosmological hydrodynamic simulations. We compare their biasing properties with the predictions of an analytic halo biasing model. Dark halos in our simulations exhibit reasonable agreement with the predictions only on scales larger than 10h^{-1}Mpc, and on smaller scales the volume exclusion effect of halos due to their finite size becomes substantial. Interestingly the biasing properties of galaxies are better described by extrapolating the halo biasing model predictions. We also find the clear dependence of galaxy biasing on their formation epoch; the distribution of old populations of galaxies tightly correlates with the underlying mass density field, while that of young populations is slightly more stochastic and anti-biased relative to dark matter. The amplitude of two-point correlation function of old populations is about 3 times larger than that of the young populations. Furthermore, the old population of galaxies reside within massive dark halos while the young galaxies are preferentially formed in smaller dark halos. Assuming that the observed early and late-type galaxies correspond to the simulated old and young populations of galaxies, respectively, all of these segregations of galaxies are consistent with observational ones for the early and late-type of galaxies such as the morphology--density relation of galaxies.

We derive an analytic formula of genus for dark matter halos assuming that the primordial mass density field obeys the random-Gaussian statistics. In particular, we for the first time take account of the nonlinear nature of the halo biasing, in addition to the nonlinear gravitational evolution of the underlying mass distribution. We examine in detail the extent to which the predicted genus for halos depends on the redshift, smoothing scale of the density field, and the mass of those halos in representative cold dark matter models. In addition to the full model predictions, we derive an explicit perturbation formula for the halo genus which can be applied to a wider class of biasing model for dark matter halos. With the empirical relation between the halo mass and the luminosity/temperature of galaxy clusters, our prediction can be compared with the upcoming cluster redshift catalogues in optical and X-ray bands. Furthermore, we discuss the detectability of the biasing and cosmological model dependence from the future data.

We quantify the degree of nonlinearity and stochasticity of the clustering of biased objects, using cosmological N-body simulations. Adopting the peaks and the halos as representative biasing models, we focus on the two-point correlation of the biased objects, dark matter and their cross-correlation. Especially, we take account of the effect of redshift-space distortion and attempt to clarify the scale-dependence and the time-dependence by analyzing the biasing factor and the cross-correlation factor. On small scales, stochasticity and nonlinearity become appreciable and strongly object-dependent, especially in redshift space due to the pair-wise velocity dispersion of the biased objects. Nevertheless, an approximation of deterministic linear biasing $\delta_{\rm obj}\simeq b_{\rm obj} \delta_{\rm mass}$ works reasonably well even in the quasi-linear regime $r > 10 h^{-1}$ Mpc, and linear redshift-space distortion explains the clustering amplitudes in redshift space in this regime.

We revisit the feasibility of a cosmological test with the geometric distortion focusing on an ambiguous factor of the evolution of bias. Starting from defining estimators for the spatial two-point correlation function and the power spectrum, in a rigorous manner, we derive useful formulas for the two-point clustering statistics which take the light-cone effect and the redshift-space distortions into account. Then we investigate how the predicted correlation functions are affected by the redshift-space distortions and the bias assuming quasar samples which roughly match the 2dF survey. Effect of the bias is crucial, in particular, on anisotropic component of the clustering statistics in redshift-space. By surveying behavior of the predicted correlation functions to parameters of a phenomenological model for the evolution of bias, it is shown that the correlation functions on the light-cone is sensitive to a mean amplitude of the bias and rather insensitive to the speed of its redshift-evolution. This feature is useful for an analysis of the geometric test.

We discuss the systematic effects arising from the cosmological redshift-space (geometric) distortion on the statistical analysis of isodensity contour using high-redshift catalogs. Especially, we present a simple theoretical model for isodensity statistics in cosmological redshift-space, as a generalization of the work by Matsubara (1996). The statistical quantities considered here are the two- and three-dimensional genus of isodensity contour, the area of surface, the length of contour intersecting with a plane and the number of the crossing points of isodensity contour on a line. We give useful analytic formulae for the isodensity statistics, which take into account the corrections from the geometric distortion, the nonlinear clustering and the nonlinear velocity distortion phenomenologically. We then demonstrate how the geometric distortion and the nonlinear corrections alter shapes of the statistical quantities on the basis of plausible cosmological models. Our results show that the nonlinear correction can be sensitive to a choice of the redshift-space coordinate as increasing the redshift. The low-dimensional quantities such as two-dimensional genus systematically yield anisotropy due to the geometric and velocity distortions and their angle-dependence shows the $10\sim20%$ difference of amplitude. Sensitivity for typical high-redshift samples are estimated in an analytic manner, and the influence of the light-cone effect for the isodensity statistics is also discussed.

We discuss the influence of the cosmological background density field on the spherical infall model. The spherical infall model has been used in the Press-Schechter formalism to evaluate the number abundance of clusters of galaxies, as well as to determine the density parameter of the universe from the infalling flow. Therefore, the understanding of collapse dynamics play a key role for extracting the cosmological information. Here, we consider the modified version of the spherical infall model. We derive the mean field equations from the Newtonian fluid equations, in which the influence of cosmological background inhomogeneity is incorporated into the averaged quantities as the {\it backreaction}. By calculating the averaged quantities explicitly, we obtain the simple expressions and find that in case of the scale-free power spectrum, the density fluctuations with the negative spectral index make the infalling velocities slow. This suggests that we underestimate the density parameter $\Omega$ when using the simple spherical infall model. In cases with the index $n>0$, the effect of background inhomogeneity could be negligible and the spherical infall model becomes the good approximation for the infalling flows. We also present a realistic example with the cold dark matter power spectrum. There, the anisotropic random velocity leads to slowing down the mean infalling velocities.

We propose a physical model for nonlinear stochastic biasing of one-point statistics resulting from the formation epoch distribution of dark halos. In contrast to previous works on the basis of extensive numerical simulations, our model provides for the first time an analytic expression for the joint probability function. Specifically we derive the joint probability function of halo and mass density contrasts from the extended Press-Schechter theory. Since this function is derived in the framework of the standard gravitational instability theory assuming the random-Gaussianity of the primordial density field alone, we expect that the basic features of the nonlinear and stochastic biasing predicted from our model are fairly generic. As representative examples, we compute the various biasing parameters in cold dark matter models as a function of a redshift and a smoothing length. Our major findings are (1) the biasing of the variance evolves strongly as redshift while its scale-dependence is generally weak and a simple linear biasing model provides a reasonable approximation roughly at $R\simgt 2(1+z)\himpc$, and (2) the stochasticity exhibits moderate scale-dependence especially on $R\simlt 20\himpc$, but is almost independent of $z$. Comparison with the previous numerical simulations shows good agreement with the above behavior, indicating that the nonlinear and stochastic nature of the halo biasing is essentially understood by taking account of the distribution of the halo mass and the formation epoch.

We consider constraints on the structure formation model based on non-Gaussian fluctuations generated during inflation, which have $\chi_m^2$ distributions. Using three data sets, the abundance of the clusters at $z=0$, moderate $z$ and the correlation length, we show that constraints on the non-Gaussianity and the amplitude of fluctuations and the density parameter can be obtained. We obtain an upper bound for $\Omega_m$ and a lower bound for the non-Gaussianity and the amplitude of the fluctuations. Using the abundance of clusters at $z \sim 0.6$, for the spectrum parameterized by cold dark matter (CDM) shape parameter $\Gamma=0.23$, we obtain an upper bound for the density parameter $\Omega_m \sim 0.5$ and lower bounds for the amplitude $\sigma_8 \sim 0.7$ and for the non-Gaussianity of fluctuations $G \sim 2$ $(m \sim 200)$, where G=1 for Gaussian.

Distribution of galaxies may be a biased tracer of the dark matter distribution and the relation between the galaxies and the total mass may be stochastic, non-linear and time-dependent. Since many observations of galaxy clustering will be done at high redshift, the time evolution of non-linear stochastic biasing would play a crucial role for the data analysis of the future sky surveys. In this paper, we develop the weakly non-linear analysis and attempt to clarify the non-linear feature of the stochastic biasing. We compute the one-loop correction of the power spectrum for the total mass, the galaxies and their cross correlation. Assuming the local functional form for the initial galaxy distribution, we investigate the time evolution of the biasing parameter and the correlation coefficient. On large scales, we first find that the time evolution of the biasing parameter could deviate from the linear prediction in presence of the initial skewness. However, the deviation can be reduced when the initial stochasticity exists. Next, we focus on the quasi-linear scales, where the non-linear growth of the total mass becomes important. It is recognized that the scale-dependence of the biasing dynamically appears and the initial stochasticity could affect the time evolution of the scale-dependence. The result is compared with the recent N-body simulation that the scale-dependence of the halo biasing can appear on relatively large scales and the biasing parameter takes the lower value on smaller scales. Qualitatively, our weakly non-linear results can explain this trend if the halo-mass biasing relation has the large scatter at high redshift.

It is believed that the biasing of the galaxies plays an important role for understanding the large-scale structure of the universe. In general, the biasing of galaxy formation could be stochastic. Furthermore, the future galaxy survey might allow us to explore the time evolution of the galaxy distribution. In this paper, the analytic study of the galaxy-mass density relation and its time evolution is presented within the framework of the stochastic biasing. In the weakly non-linear regime, we derive a general formula for the galaxy-mass density relation as a conditional mean using the Edgeworth expansion. The resulting expression contains the joint moments of the total mass and galaxy distributions. Using the perturbation theory, we investigate the time evolution of the joint moments and examine the influence of the initial stochasticity on the galaxy-mass density relation. The analysis shows that the galaxy-mass density relation could be well-approximated by the linear relation. Compared with the skewness of the galaxy distribution, we find that the estimation of the higher order moments using the conditional mean could be affected by the stochasticity. Therefore, the galaxy-mass density relation as a conditional mean should be used with a caution as a tool for estimating the skewness and the kurtosis.

We study the adiabatic density perturbation in the {\it oscillating inflation}, proposed by Damour and Mukhanov. The recent study of the cosmological perturbation during reheating shows that the adiabatic fluctuation behaves like as the perfect fluid and no significant amplification occurs on super-horizon scales. In the oscillating inflation, however, the accelerated expansion takes place during the oscillating stage and there might be a possibility that the parametric amplification on small scales affects the adiabatic long-wavelength perturbation. We analytically show that the density perturbation neglecting the metric perturbation can be amplified by the parametric resonance and the instability band becomes very broad during the oscillating inflation. We examined this issue by solving the evolution equation for perturbation numerically. We found that the parametric resonance is strongly suppressed for the long wave modes comparable to the Hubble horizon. The result indicates that the metric perturbation plays a crucial role for the evolution of scalar field perturbation. Therefore, in the single field case, there would be no significant imprint of the oscillating inflation on the primordial spectrum of the adiabatic perturbation. However, it could be expected that the oscillating inflation in the multi-field system gives the enormous amplification on large scales, which may lead to the production of the primordial black holes.

It is generally believed that the spatial distribution of galaxies does not trace that of the total mass. The understanding of the bias effect is therefore necessary to determine the cosmological parameters and the primordial density fluctuation spectrum from the galaxy survey. The deterministic description of bias may not be appropriate because of the various stochasticity of galaxy formation process. In nature, the biasing is epoch dependent and recent deep survey of the galaxy shows the large biasing at high redshift. Hence, we investigate quasi non-linear evolution of the stochastic bias by using the tree level perturbation method. Especially, the influence of the initial cross correlation on the evolution of the skewness and the bi-spectrum is examined in detail. We find that the non-linear bias can be generated dynamically. The small value of the initial cross correlation can bend the $\dg$-$\dm$ relation effectively and easily lead to the negative curvature ($b_2<0$). We also propose a method to predict the bias, cross correlation and skewness at high redshift. As an illustration, the possibility of the large biasing at high redshift is discussed. Provided the present bias parameter as $b=1.5$ and $\Omega=1.0$, we predict the large scale bias as $b=4.63$ at $z=3$ by fitting the bi-spectrum to the Lick catalog data. Our results will be important for the future deep sky survey.

We investigate gauge invariant cosmological perturbations in a spatially flat Friedman-Robertson-Walker universe with scalar fields. It is well known that the evolution equation for the gauge invariant quantities has exact solutions in the long-wavelength limit. We find that these gauge invariant solutions can be obtained by differentiating the background solution with respect to parameters contained in the background system. This method is very useful when we analyze the long-wavelength behavior of cosmological perturbation with multiple scalar fields.

We investigate the cosmological perturbation of two-scalar field model during the reheating phase after inflation. Using the exact solution of the perturbation in long-wavelength limit, which is expressed in terms of the background quantities, we analyze the behavior of the metric perturbation. The oscillating inflaton field gives rise to the parametric resonance of the massless scalar field and this leads to the amplification of the iso-curvature mode of the metric perturbations.

The evolution of the cosmological perturbation during the oscillatory stage of the scalar field is investigated. For the power law potential of the inflaton field, the evolution equation of the Mukhanov's gauge invariant variable is reduced to the Mathieu equation and the density perturbation grows by the parametric resonance.

As an application of Gradient Expansion (long-wavelength) approximation, we studied the inhomogeneous universe including the gravitational wave(GW). For a plane-symmetric cosmological model, we could implement the 2nd order expansion and analyze the purely non-linear perturbation of GW. The result indicates that the non-linear effect plays an important role in anisotropic stage of the universe. We also confirm that the perturbation from the linear source term in gradient expansion is consistent with the result of linear perturbation.

Using the long wave perturbation scheme(gradient expansion), the effect of inhomogeneity on the inflationary phase is investigated. We solved the perturbation equation of which source term comes from inhomogeneity of a scalar field and a seed metric. The result indicates that sub-horizon scale inhomogeneity strongly affects the onset of inflation.