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.

Using $N$-body simulations for the first time, we show that the anisotropic primordial non-Gaussianity (PNG) causes a scale-dependent modification, given by $1/k^2$ at small $k$ limit, in the three-dimensional power spectra of halo shapes (intrinsic alignments), whilst the conventional power spectrum of halo number density field remains unaffected. We discuss that wide-area imaging and spectroscopic surveys observing the same region of the sky allow us to constrain the quadrupole PNG coefficient $f_{\rm NL}^{s=2}$ at a precision comparable with that of the cosmic microwave background.

The spatial clustering of active galactic nuclei (AGNs) is considered to be one of the important diagnostics for the understanding of the underlying processes behind their activities complementary to measurements of the luminosity function (LF). We analyse the AGN clustering from a recent semi-analytic model performed on a large cosmological $N$-body simulation covering a cubic gigaparsec comoving volume. We have introduced a new time-scale of gas accretion on to the supermassive black holes to account for the loss of the angular momentum on small scales, which is required to match the faint end of the observed X-ray LF. The large simulation box allows us accurate determination of the auto-correlation function of the AGNs. The model prediction indicates that this time-scale plays a significant role in allowing massive haloes to host relatively faint population of AGNs, leading to a higher bias factor for those AGNs. The model predictions are in agreement with observations of X-ray selected AGNs in the luminosity range $10^{41.5}~\mathrm{erg} \ \mathrm{s}^{-1} \leq L_{2-10\mathrm{keV}} \leq 10^{44.5}~\mathrm{erg} \ \mathrm{s}^{-1}$, with the typical host halo mass of $10^{12.5-13.5} h^{-1}\,{\rm M}_{\odot}$ at $z \lesssim 1$. This result shows that the observational clustering measurements impose an independent constraint on the accretion time-scale complementary to the LF measurements. Moreover, we find that not only the effective halo mass corresponding to the overall bias factor, but the extended shape of the predicted AGN correlation function shows remarkable agreement with those from observations. Further observational efforts towards the low luminosity end at $z \sim 1$ would give us stronger constraints on the triggering mechanisms of AGN activities through their clustering.

An optical cluster finder inevitably suffers from projection effects, where it misidentifies a superposition of galaxies in multiple halos along the line-of-sight as a single cluster. Using mock cluster catalogs built from cosmological N-body simulations, we quantify the impact of these projection effects with a particular focus on the observables of interest for cluster cosmology, namely the cluster lensing and the cluster clustering signals. We find that "observed" clusters, i.e. clusters identified by our cluster finder algorithm, exhibit lensing and clustering signals that deviate from expectations based on a statistically isotropic halo model -- while both signals agree with halo model expectations on small scales, they show unexpected boosts on large scales, by up to a factor of 1.2 or 1.4 respectively. We identify the origin of these boosts as the inherent selection bias of optical cluster finders for clusters embedded within filaments aligned with the line-of-sight, and show that a minority ($\sim 30\%$) of such clusters within the entire sample is responsible for this observed boost. We discuss the implications of our results on previous studies of optical cluster, as well as prospects for identifying and mitigating projection effects in future cluster cosmology analyses.

The long-wavelength coherent overdensity and tidal force, which are not direct observables for a finite-volume survey, affect time evolution of cosmic structure formation and therefore clustering observables through the mode coupling. In this paper we develop an "anisotropic" separate universe (SU) simulation technique to simulate large-scale structure formation taking into account the effect of large-scale tidal force into the anisotropic expansion of local background. We modify the TreePM N-body simulation code to implement the anisotropic SU simulations, and then study the "response" function of matter power spectrum that describes how the matter power spectrum responds to the large-scale tidal effect as a function of wavenumber and redshift for a given global cosmology. We test and validate the SU simulation results from the comparison with the perturbation theory predictions and the results from high-resolution PM simulations. We find that the response function displays characteristic scale dependences over the range of scales down to nonlinear scales, up to k ~ 6 h/Mpc.

We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter.

The Consistency Relations for the Large Scale Structure provide a link between the amplitude of Baryonic Acoustic Oscillations in the squeezed bispectrum (BS) and in the power spectrum (PS). This relation depends on the large scale bias of the considered tracer, $b_\alpha$, and on the growth rate of structures, $f$. Remarkably, originating from basic symmetry principles, this relation is exact and independent on the underlying cosmological model. By analysing data from large volume simulations, both for dark matter and for haloes, we illustrate how BS and PS measurements can be used to extract $b_\alpha$ and $f$ without the need of any theoretical approximation scheme for the computation of the BS and the PS. We show that, combining measurements of the squeezed BS with the quadrupole to monopole ratios for the PS at large scales can successfully break the $b_\alpha -f$ degeneracy. We forecast that this method, applied to a Euclid-like survey, will be able to measure bias, and then the growth rate, at better than $10\%$ level, with no extra assumption.

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.

An accurate theoretical template of the redshift-space galaxy power spectrum, if applicable out to nonlinear scales, enables us to extract more stringent and robust constraints on cosmological parameters from the measured galaxy clustering. In this work we develop a simulation-based template, so-called emulator, for the redshift-space power spectrum of dark matter halos. Using the redshift-space halo power spectra measured from the Dark Quest $N$-body simulation suite that covers 101 flat-geometry $w$CDM cosmologies around the Planck $\Lambda$CDM model, we feed these data into a feed-forward neural network to build the fast and accurate emulation of the power spectrum from the linear to nonlinear scales up to $k \simeq 0.6 \, h \, {\rm Mpc}^{-1}$. Our emulator achieves about 1 and 5% fractional accuracies in predicting the monopole and quadrupole moments of the power spectrum, respectively, for halos of $\sim 10^{13}h^{-1}M_\odot$ that correspond to host halos of the SDSS LOWZ- and CMASS-like galaxies, where the achieved accuracies are sufficient compared to the statistical errors of SDSS volume. The validation and performance of the emulator are given by the comparison of the emulator predictions with the power spectra directly measured from the simulations for validation sets that are not used in the training. We demonstrate that the emulator outputs can be used to make model predictions for the redshift-space power spectrum of galaxies by employing user-fed models for the halo-galaxy connection, such as the halo occupation distribution. The emulator allows us to easily incorporate the Finger-of-God effect due to the virial motions of galaxies and the Alcock--Paczy\'{n}ski distortions. Our code can compute the redshift-space galaxy power spectrum in a CPU subseconds, and is ready to perform the emulator-based cosmological analysis for the exiting and upcoming galaxy redshift surveys.

Consistency Relations (CR) for the Large Scale Structure are exact equalities between correlation functions of different order. These relations descend from the equivalence principle and hold for primordial perturbations generated by single-field models of inflation. They are not affected by nonlinearities and hold also for biased tracers and in redshift space. We show that Baryonic Acoustic Oscillations (BAO) in the bispectrum (BS) in the squeezed limit are suppressed with respect to those in the power spectrum (PS) by a coefficient that depends on the BS configuration and on the bias parameter (and, in redshift space, also on the growth rate). We test these relations using large volume N-body simulations and show that they provide a novel way to measure large scale halo bias and, potentially, the growth rate. Since bias is obtained by comparing two directly observable quantities, the method is free from theoretical uncertainties both on the computational scheme and on the underlying cosmological model.

We study the predictions for the matter redshift-space power spectrum and correlation function of a Lagrangian-space Gaussian ansatz introduced in a previous work. This model is a natural extension of the Zeldovich approximation, where the displacement and velocity power spectra are determined by the equations of motion, instead of being set equal to the linear power spectrum. It does not contain any free parameter. As for the real-space statistics, we find that this Lagrangian-space approach is much more efficient for the correlation functions than for the power spectra. The damping of the BAO oscillations is well recovered but there is a large smooth drift from the simulations in the power spectra. The multipoles of the correlation functions are well recovered on BAO scales, with an accuracy of $2\%$ for $\xi^s_0$ down to $10 h^{-1}$ Mpc, and of $3\%$ for $\xi^s_2$ down to $26 h^{-1}$ Mpc, at $z \geq 0.35$.

We use a suite of $N$-body simulations to study intrinsic alignments (IA) of halo shapes with the surrounding large-scale structure in the $\Lambda$CDM model. For this purpose, we develop a novel method to measure multipole moments of the three-dimensional power spectrum of the $E$-mode field of halo shapes with the matter/halo distribution, $P_{\delta E}^{(\ell)}(k)$ (or $P^{(\ell)}_{{\rm h}E}$), and those of the auto-power spectrum of the $E$ mode, $P^{(\ell)}_{EE}(k)$, based on the $E$/$B$-mode decomposition. The IA power spectra have non-vanishing amplitudes over the linear to nonlinear scales, and the large-scale amplitudes at $k\lesssim 0.1~h~{\rm Mpc}^{-1}$ are related to the matter power spectrum via a constant coefficient ($A_{\rm IA}$), similar to the linear bias parameter. We find that the cross- and auto-power spectra $P_{\delta E}$ and $P_{EE}$ at nonlinear scales, $k\gtrsim 0.1~h~{\rm Mpc}^{-1}$, show different $k$-dependences relative to the matter power spectrum, suggesting a violation of the nonlinear alignment model commonly used to model contaminations of cosmic shear signals. The IA power spectra exhibit baryon acoustic oscillations, and vary with halo samples of different masses, redshifts and cosmological parameters ($\Omega_{\rm m}, S_8$). The cumulative signal-to-noise ratio for the IA power spectra is about 60\% of that for the halo density power spectrum, where the super-sample covariance is found to give a significant contribution to the total covariance. Our results demonstrate that the IA power spectra of galaxy shapes, measured from imaging and spectroscopic surveys for an overlapping area of the sky, can be powerful tools to probe the underlying matter power spectrum, the primordial curvature perturbations, and cosmological parameters, in addition to the standard galaxy density power spectrum.

Recent constraints on the splashback radius around optically selected galaxy clusters from the redMaPPer cluster-finding algorithm in the literature have shown that the observed splashback radius is $\sim 20\%$ smaller than that predicted by N-body simulations. We present analyses on the splashback features around $\sim 3000$ optically selected galaxy clusters detected by the independent cluster-finding algorithm CAMIRA over a wide redshift range of $0.1<z_{\rm cl}<1.0$ from the second public data release of the Hyper Suprime-Cam (HSC) Subaru Strategic Program covering $\sim 427~{\rm deg}^2$ for the cluster catalog. We detect the splashback feature from the projected cross-correlation measurements between the clusters and photometric galaxies over the wide redshift range, including for high redshift clusters at $0.7<z_{\rm cl}<1.0$, thanks to deep HSC images. We find that constraints from red galaxy populations only are more precise than those without any color cut, leading to $1\sigma$ precisions of $\sim 15\%$ at $0.4<z_{\rm cl}<0.7$ and $0.7<z_{\rm cl}<1.0$. These constraints are more consistent with the model predictions ($\lesssim 1\sigma$) than their $20\%$ smaller values as suggested by the previous studies with the redMaPPer ($\sim 2\sigma$). We also investigate selection effects of the optical cluster-finding algorithms on the observed splashback features by creating mock galaxy catalogs from a halo occupation distribution model, and find that such effects to be sub-dominant for the CAMIRA cluster-finding algorithm. We also find that the redMaPPer-like cluster-finding algorithm induces a smaller inferred splashback radius in our mock catalog, especially at lower richness, which can well explain the smaller splashback radii in the literature. In contrast, these biases are significantly reduced when increasing its aperture size.

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.

An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.

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.

To understand the nature of the accelerated expansion of the Universe, we need to combine constraints on the expansion rate and growth of structure. The growth rate is usually extracted from three dimensional galaxy maps by exploiting the effects of peculiar motions on galaxy clustering. However, theoretical models of the probability distribution function (PDF) of galaxy pairwise peculiar velocities are not accurate enough on small scales to reduce the error on theoretical predictions to the level required to match the precision expected for measurements from future surveys. Here, we improve the modelling of the pairwise velocity distribution by using the Skew-T PDF, which has nonzero skewness and kurtosis. Our model accurately reproduces the redshift-space multipoles (monopole, quadrupole and hexadecapole) predicted by N-body simulations, above scales of about $10\,h^{-1}{\rm Mpc}$. We illustrate how a Taylor expansion of the streaming model can reveal the contributions of the different moments to the clustering multipoles, which are independent of the shape of the velocity PDF. The Taylor expansion explains why the Gaussian streaming model works well in predicting the first two redshift-space multipoles, although the velocity PDF is non-Gaussian even on large scales. Indeed, any PDF with the correct first two moments would produce precise results for the monopole down to scales of about $10\,h^{-1}{\rm Mpc}$, and for the quadrupole down to about $30\,h^{-1}{\rm Mpc}$. An accurate model for the hexadecapole needs to include higher-order moments.

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.

Clustering properties and peculiar velocities of halos in large-scale structure carry a wealth of cosmological information over a wide range of scales from linear to nonlinear scales. We use halo catalogs in a suite of high-resolution $N$-body simulations to construct mock catalogs of galaxies that resemble the SDSS-like luminous early-type galaxies at three redshift bins in the range $0.15\le z\le 0.7$. To do this we include 10 nuisance parameters to model variations in halo-galaxy connections for each redshift bin; the halo occupation distribution, and the spatial and velocity distributions of galaxies in the host halos. We evaluate the Fisher information matrix for the redshift-space power spectrum of SDSS-like galaxies using different sets of the mock catalogs that are generated from changes in each of model parameters; cosmological parameters ($\sigma_8$ and $\Omega_{\rm m}$), the halo-galaxy connection parameters, and the cosmological distances ($D_{\rm A}$ and $H$ parameters at each redshift bin) for modeling an apparent geometrical distortion of the redshift-space power spectrum (the Alcock-Paczynski effect). We show that combining the monopole and quadrupole power spectra of galaxies allows for accurate estimations of the cosmological parameters and the cosmological distances, even after marginalization over the halo-galaxy parameters, by lifting the parameter degeneracies that are otherwise inevitable if either of the two spectra alone is used. When including the galaxy power spectrum information up to $k=0.3\,h\,$Mpc$^{-1}$, we find about factor of 6 gain in the cosmological information content of ($\sigma_8,\Omega_{\rm m},D_{\rm A}$'s and $H$'s) compared to $k=0.2\,h\,$Mpc$^{-1}$.

Detecting primordial non-Gaussianity on mildly non-linear scales requires precise modelling of late-time structure formation. Accurately predicting the impact of non-linear gravitational collapse, non-linear tracer bias and baryonic physics on the variance and higher order moments of the cosmic density field is challenging, as they strongly depend on the tails of the probability distribution function (PDF) of density fluctuations. A way around this problem is to directly analyse the bulk of the PDF instead. For this purpose we devise a new method to predict the impact of general non-Gaussian initial conditions on the late-time density PDF. With this formalism we show that - even when marginalizing over potential ignorance of the amplitude and slope of the non-linear power spectrum - an analysis of the PDF at mildly non-linear densities can measure the amplitude of different primordial bispectrum shapes to an accuracy of $\Delta f_{\mathrm{NL}}^{\mathrm{loc}} = \pm 3.1\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{equi}} = \pm 10.0\ ,\ \Delta f_{\mathrm{NL}}^{\mathrm{ortho}} = \pm 17.0\ $. This assumes a joint analysis of the PDF on smoothing scales of $15$Mpc/$h$ and $30$Mpc/$h$ in a survey volume of $V=100(\mathrm{Gpc}/h)^3$ at $z=1$, analysing only densities of $\delta(15\mathrm{Mpc}/h) \in [-0.4, 0.5]$ ($\approx 87\%$ of probability) and $\delta(30\mathrm{Mpc}/h) \in [-0.3, 0.4]$ ($\approx 95\%$ of probability). Note that a formalism closely related to ours was already successfully applied to observational data \citep{Gruen2018, Friedrich2018}, demonstrating that the methodology developed here can indeed be carried over to real data analysis.

We examine the impact of fiber assignment on clustering measurements from fiber-fed spectroscopic galaxy surveys. We identify new effects which were absent in previous, relatively shallow galaxy surveys such as Baryon Oscillation Spectroscopic Survey . Specifically, we consider deep surveys covering a wide redshift range from z=0.6 to z=2.4, as in the Subaru Prime Focus Spectrograph survey. Such surveys will have more target galaxies than we can place fibers on. This leads to two effects. First, it eliminates fluctuations with wavelengths longer than the size of the field of view, as the number of observed galaxies per field is nearly fixed to the number of available fibers. We find that we can recover the long-wavelength fluctuation by weighting galaxies in each field by the number of target galaxies. Second, it makes the preferential selection of galaxies in under-dense regions. We mitigate this effect by weighting galaxies using the so-called individual inverse probability. Correcting these two effects, we recover the underlying correlation function at better than 1 percent accuracy on scales greater than 10 Mpc/h.

Galaxy imaging surveys enable us to map the cosmic matter density field through weak gravitational lensing analysis. The density reconstruction is compromised by a variety of noise originating from observational conditions, galaxy number density fluctuations, and intrinsic galaxy properties. We propose a deep-learning approach based on generative adversarial networks (GANs) to reduce the noise in the weak lensing map under realistic conditions. We perform image-to-image translation using conditional GANs in order to produce noiseless lensing maps using the first-year data of the Subaru Hyper Suprime-Cam (HSC) survey. We train the conditional GANs by using 30000 sets of mock HSC catalogs that directly incorporate observational effects. We show that an ensemble learning method with GANs can reproduce the one-point probability distribution function (PDF) of the lensing convergence map within a $0.5-1\sigma$ level. We use the reconstructed PDFs to estimate a cosmological parameter $S_{8} = \sigma_{8}\sqrt{\Omega_{\rm m0}/0.3}$, where $\Omega_{\rm m0}$ and $\sigma_{8}$ represent the mean and the scatter in the cosmic matter density. The reconstructed PDFs place tighter constraint, with the statistical uncertainty in $S_8$ reduced by a factor of $2$ compared to the noisy PDF. This is equivalent to increasing the survey area by $4$ without denoising by GANs. Finally, we apply our denoising method to the first-year HSC data, to place $2\sigma$-level cosmological constraints of $S_{8} < 0.777 \, ({\rm stat}) + 0.105 \, ({\rm sys})$ and $S_{8} < 0.633 \, ({\rm stat}) + 0.114 \, ({\rm sys})$ for the noisy and denoised data, respectively.

It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between $0<z<5$. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, $z>1$. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within $\sim 30^{\circ}$ in the three dimensional space and $\sim 20^{\circ}$ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the {\it slowest collapsing} (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch.

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.

Constraining the relation between the richness $N$ and the halo mass $M$ over a wide redshift range for optically-selected clusters is a key ingredient for cluster-related science in optical surveys, including the Subaru Hyper Suprime-Cam (HSC) survey. We measure stacked weak lensing profiles around 1747 HSC CAMIRA clusters over a redshift range of $0.1\leq z_{\rm cl}\leq 1.0$ with $N\geq 15$ using the HSC first-year shear catalog covering $\sim$$140$ ${\rm deg^2}$. The exquisite depth and image quality of the HSC survey allow us to measure lensing signals around the high-redshift clusters at $0.7\leq z_{\rm cl}\leq 1.0$ with a signal-to-noise ratio of 19 in the comoving radius range $0.5\lesssim R\lesssim 15 h^{-1}{\rm Mpc}$. We constrain richness-mass relations $P(\ln N|M,z)$ of the HSC CAMIRA clusters assuming a log-normal distribution without informative priors on model parameters, by jointly fitting to the lensing profiles and abundance measurements under both Planck and WMAP cosmological models. We show that our model gives acceptable $p$-values when we add redshift dependent terms which are proportional to $\ln (1+z)$ and $[\ln (1+z)]^{2}$ into the mean and scatter relations of $P(\ln N|M,z)$. Such terms presumably originate from the variation of photometric redshift errors as a function of the redshift. We show that the constraints on the mean relation $\langle M|N \rangle$ are consistent between the Planck and WMAP models, whereas the scatter values $\sigma_{\ln M|N}$ for the Planck model are systematically larger than those for the WMAP model. We also show that the scatter values for the Planck model increase toward lower richness values, whereas those for the WMAP model are consistent with constant values as a function of richness. This result highlights the importance of the scatter in the mass-richness relation for cluster cosmology.

We perform an ensemble of $N$-body simulations with $2048^3$ particles for 101 flat $w$CDM cosmological models sampled based on a maximin-distance Sliced Latin Hypercube Design. By using the halo catalogs extracted at multiple redshifts in the range of $z=[0,1.48]$, we develop Dark Emulator, which enables fast and accurate computations of the halo mass function, halo-matter cross-correlation, and halo auto-correlation as a function of halo masses, redshift, separations and cosmological models, based on the Principal Component Analysis and the Gaussian Process Regression for the large-dimensional input and output data vector. We assess the performance of the emulator using a validation set of $N$-body simulations that are not used in training the emulator. We show that, for typical halos hosting CMASS galaxies in the Sloan Digital Sky Survey, the emulator predicts the halo-matter cross correlation, relevant for galaxy-galaxy weak lensing, with an accuracy better than $2\%$ and the halo auto-correlation, relevant for galaxy clustering correlation, with an accuracy better than $4\%$. We give several demonstrations of the emulator. It can be used to study properties of halo mass density profiles such as the mass-concentration relation and splashback radius for different cosmologies. The emulator outputs can be combined with an analytical prescription of halo-galaxy connection such as the halo occupation distribution at the equation level, instead of using the mock catalogs, to make accurate predictions of galaxy clustering statistics such as the galaxy-galaxy weak lensing and the projected correlation function for any model within the $w$CDM cosmologies, in a few CPU seconds.

The distributions of the pairwise line-of-sight velocity between galaxies and their host clusters are segregated according to the galaxy's colour and morphology. We investigate the velocity distribution of red-spiral galaxies, which represents a rare population within galaxy clusters. We find that the probability distribution function of the pairwise line-of-sight velocity $v_{\rm{los}}$ between red-spiral galaxies and galaxy clusters has a dip at $v_{\rm{los}} = 0$, which is a very odd feature, at 93\% confidence level. To understand its origin, we construct a model of the phase space distribution of galaxies surrounding galaxy clusters in three-dimensional space by using cosmological $N$-body simulations. We adopt a two component model that consists of the infall component, which corresponds to galaxies that are now falling into galaxy clusters, and the splashback component, which corresponds to galaxies that are on their first (or more) orbit after falling into galaxy clusters. We find that we can reproduce the distribution of the line-of-sight velocity of red-spiral galaxies with the dip with a very simple assumption that red-spiral galaxies reside predominantly in the infall component, regardless of the choice of the functional form of their spatial distribution. Our results constrain the quenching timescale of red-spiral galaxies to a few Gyrs, and the radius where the morphological transformation is effective as $r \sim 0.2 h^{-1} \rm{Mpc}$.

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.

The galaxy bispectrum, measuring excess clustering of galaxy triplets, offers a probe of dark energy via baryon acoustic oscillations (BAOs). However up to now it has been severely underused due to the combinatorically explosive number of triangles. Here we exploit interference in the bispectrum to identify triangles that amplify BAOs. This approach reduces the computational cost of estimating covariance matrices, offers an improvement in BAO constraints equivalent to lengthening BOSS by 30%, and simplifies adding bispectrum BAO information to future large-scale redshift survey analyses.

We use a simple cosmological model with two parameters $(A_s, n_s)$ to illustrate the impact of using Halofit on error forecast based on Fisher information matrix for a $h^{-3} Gpc^3$ volume survey. We show that Halofit fails to reproduce well the derivatives of the power spectrum with respect to the cosmological parameters despite the good fit produced for its amplitude. We argue that the poor performance on the derivatives prediction is a general feature of this model and we exhibit the response function for the Halofit to show how it compares with the same quantity measured on simulations. The analytic structure of the Halofit response function points towards the origin of its weak performance at reproducing the derivatives of the non-linear power spectrum, which translates into unreliable Fisher information matrices.

In this study, we measure the response of matter and halo projected power spectra $P^{\rm 2D}_{\rm XY}(k)$ (X, Y are matter and/or halos), to a large-scale density contrast, $\delta_{\rm b}$, using separate universe simulations. We show that the fractional response functions, i.e., $\mathrm{d}\ln P^{\rm 2D}_{\rm XY}(k)/\mathrm{d}\delta_{\rm b}$, are identical to their respective three-dimensional power spectra within simulation measurement errors. Then, using various $N$-body simulation combinations (small-box simulations with periodic boundary conditions and sub-volumes of large-box simulations) to construct {mock observations of projected fields}, we study how super-survey modes, in both parallel and perpendicular directions to the projection direction, affect the covariance matrix of $P^{\rm 2D}_{\rm XY}(k)$, known as super-sample covariance (SSC). Our results indicate that the SSC term provides dominant contributions to the covariances of matter-matter and matter-halo spectra at small scales but does not provide significant contributions in the halo-halo spectrum. We observe that the large-scale density contrast in each redshift shell causes most of the SSC effect, and we did not observe a SSC signature arising from large-scale tidal field within the levels of measurement accuracy. We also develop a response approach to calibrate the SSC term for cosmic shear correlation function and galaxy--galaxy weak lensing, and validate the method by comparison with the light-cone ray-tracing simulations. Our method provides a reasonably accurate, albeit computationally inexpensive, way to calibrate the covariance matrix for clustering observables available from wide-area galaxy surveys.

We develop a novel approach in exploring the joint dependence of halo bias on multiple halo properties using Gaussian process regression. Using a $\Lambda$CDM $N$-body simulation, we carry out a comprehensive study of the joint bias dependence on halo structure, formation history and environment. We show that the bias is a multivariate function of halo properties that falls into three regimes. For massive haloes, halo mass explains the majority of bias variation. For early-forming haloes, bias depends sensitively on the recent mass accretion history. For low-mass and late-forming haloes, bias depends more on the structure of a halo such as its shape and spin. Our framework enables us to convincingly prove that $V_\mathrm{max}/V_\mathrm{vir}$ is a lossy proxy of formation time for bias modelling, whereas the mass, spin, shape and formation time variables are non-redundant with respect to each other. Combining mass and formation time largely accounts for the mass accretion history dependence of bias. Combining all the internal halo properties fully accounts for the density profile dependence inside haloes, and predicts the clustering variation of individual haloes to a $20\%$ level at $\sim 10\mathrm{Mpc}h^{-1}$. When an environmental density is measured outside $1\mathrm{Mpc}h^{-1}$ from the halo centre, it outperforms and largely accounts for the bias dependence on the internal halo structure, explaining the bias variation above a level of $30\%$.

Velocity dispersions have been employed as a method to measure masses of clusters. To complement this conventional method, we explore the possibility of constraining cluster masses from the stacked phase space distribution of galaxies at larger radii, where infall velocities are expected to have a sensitivity to cluster masses. First, we construct a two component model of the three-dimensional phase space distribution of haloes surrounding clusters up to 50 $h^{-1}$Mpc from cluster centres based on $N$-body simulations. We find that the three-dimensional phase space distribution shows a clear cluster mass dependence up to the largest scale examined. We then calculate the probability distribution function of pairwise line-of-sight velocities between clusters and haloes by projecting the three-dimensional phase space distribution along the line-of-sight with the effect of the Hubble flow. We find that this projected phase space distribution, which can directly be compared with observations, shows a complex mass dependence due to the interplay between infall velocities and the Hubble flow. Using this model, we estimate the accuracy of dynamical mass measurements from the projected phase space distribution at the transverse distance from cluster centres larger than $2h^{-1}$Mpc. We estimate that, by using $1.5\times 10^5$ spectroscopic galaxies, we can constrain the mean cluster masses with an accuracy of 14.5\% if we fully take account of the systematic error coming from the inaccuracy of our model. This can be improved down to 5.7\% by improving the accuracy of the model.

We investigate the splashback features of dark-matter halos based on cosmic density and velocity fields. Besides the density correlation function binned by the halo orientation angle which was used in the literature, we introduce, for the first time, the corresponding velocity statistic, alignment momentum correlation function, to take into account the asphericity of halos. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity. On halo scales, $x\sim R_\mathrm{200m} \sim 1Mpc/h$, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. We also find that the splashback radius determined from the density correlation becomes $\sim 3.5\%$ smaller than that from the momentum correlation, with their correlation coefficient being 0.605. Moreover, the orientation-dependent splashback feature due to halo asphericity is measured when the density profile is determined by dark-matter particles, which can be used as a test of collisional CDM since the halo shape is predicted to be rounder in such a model.

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.

While various observations measured ellipticities of galaxy clusters and alignments between orientations of the brightest cluster galaxies and their host clusters, there are only a handful of numerical simulations that implement realistic baryon physics to allow direct comparisons with those observations. Here we investigate ellipticities of galaxy clusters and alignments between various components of them and the central galaxies in the state-of-the-art cosmological hydrodynamical simulation Horizon-AGN, which contains dark matter, stellar, and gas components in a large simulation box of $(100 h^{-1}$ Mpc$)^3$ with high spatial resolution ($\sim1$ kpc). We estimate ellipticities of total matter, dark matter, stellar, gas surface mass density distributions, X-ray surface brightness, and the Compton $y$-parameter of the Sunyaev-Zel'dovich effect, as well as alignments between these components and the central galaxies for 120 projected images of galaxy clusters with masses $M_{200}>5\times10^{13}M_{\odot}$. Our results indicate that the distributions of these components are well aligned with the major-axes of the central galaxies, with the root mean square value of differences of their position angles of $\sim 20^\circ$, which vary little from inner to the outer regions. We also estimate alignments of these various components with total matter distributions, and find tighter alignments than those for central galaxies with the root mean square value of $\sim 15^\circ$. We compare our results with previous observations of ellipticities and position angle alignments and find reasonable agreements. The comprehensive analysis presented in this paper provides useful prior information for analyzing stacked lensing signals as well as designing future observations to study ellipticities and alignments of galaxy clusters.

We present weak-lensing measurements using the first-year data from the Hyper Suprime-Cam Strategic Survey Program on the Subaru telescope for eight galaxy clusters selected through their thermal Sunyaev-Zel'dovich (SZ) signal measured at 148 GHz with the Atacama Cosmology Telescope Polarimeter experiment. The overlap between the two surveys in this work is 33.8 square degrees, before masking bright stars. The signal-to-noise ratio of individual cluster lensing measurements ranges from 2.2 to 8.7, with a total of 11.1 for the stacked cluster weak-lensing signal. We fit for an average weak-lensing mass distribution using three different profiles, a Navarro-Frenk-White profile, a dark-matter-only emulated profile, and a full cosmological hydrodynamic emulated profile. We interpret the differences among the masses inferred by these models as a systematic error of 10\%, which is currently smaller than the statistical error. We obtain the ratio of the SZ-estimated mass to the lensing-estimated mass (the so-called hydrostatic mass bias $1-b$) of $0.74^{+0.13}_{-0.12}$, which is comparable to previous SZ-selected clusters from the Atacama Cosmology Telescope and from the {\sl Planck} Satellite. We conclude with a discussion of the implications for cosmological parameters inferred from cluster abundances compared to cosmic microwave background primary anisotropy measurements.

We study the dependence of surface mass density profiles, which can be directly measured by weak gravitational lensing, on the orientation of haloes with respect to the line-of-sight direction, using a suite of $N$-body simulations. We find that, when major axes of haloes are aligned with the line-of-sight direction, surface mass density profiles have higher amplitudes than those averaged over all halo orientations, over all scales from $0.1$ to $100\,\mathrm{Mpc}/h$ we studied. While the orientation dependence at small scales is ascribed to the halo triaxiality, our results indicate even stronger orientation dependence in the so-called two-halo regime, up to $100\,\mathrm{Mpc}/h$. The orientation dependence for the two-halo term is well approximated by a multiplicative shift of the amplitude and therefore a shift in the halo bias parameter value. The halo bias from the two-halo term can be overestimated or underestimated by up to $\sim 30 \%$ depending on the viewing angle, which translates into the bias in estimated halo masses by up to a factor of two from halo bias measurements. The orientation dependence at large scales originates from the anisotropic halo-matter correlation function, which has an elliptical shape with the axis ratio of $\sim 0.55$ up to $100\, \mathrm{Mpc}/h$. We discuss potential impacts of halo orientation bias on other observables such as optically selected cluster samples and a clustering analysis of large-scale structure tracers such as quasars.

We constrain the scaling relation between optical richness ($\lambda$) and halo mass ($M$) for a sample of SDSS redMaPPer galaxy clusters within the context of the {\it Planck} cosmological model. We use a forward modeling approach where we model the probability distribution of optical richness for a given mass, $P(\ln \lambda| M)$. To model the abundance and the stacked lensing profiles, we use an emulator specifically built to interpolate the halo mass function and the stacked lensing profile for an arbitrary set of halo mass and redshift, which is calibrated based on a suite of high-resolution $N$-body simulations. We apply our method to 8,312 SDSS redMaPPer clusters with $20\le \lambda \le 100$ and $0.10\le z_{\lambda}\le0.33$, and show that the log-normal distribution model for $P(\lambda|M)$, with four free parameters, well reproduces the measured abundances and lensing profiles simultaneously. The constraints are characterized by the mean relation, $\left\langle \ln{\lambda}\right\rangle(M)=A+B\ln(M/M_{\rm pivot})$, with $A=3.207^{+0.044}_{-0.046}$ and $B=0.993^{+0.041}_{-0.055}$ (68\%~CL), where the pivot mass scale $M_{\rm pivot}=3\times 10^{14} h^{-1}M_\odot$, and the scatter $\sigma_{\mathrm{\ln\lambda}|M}=\sigma_0+q\ln(M/M_{\rm pivot})$ with $\sigma_0=0.456^{+0.047}_{-0.039}$ and $q=-0.169^{+0.035}_{-0.026}$. We find that a large scatter in halo masses is required at the lowest richness bins ($20\le \lambda \lesssim 30$) in order to reproduce the measurements. Without such a large scatter, the model prediction for the lensing profiles tends to overestimate the measured amplitudes. This might imply a possible contamination of intrinsically low-richness clusters due to the projection effects. Such a low-mass halo contribution is significantly reduced when applying our method to the sample of $30\le \lambda \le 100$.

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.

We study a new procedure to measure the sound horizon scale via Baryonic Acoustic Oscillations (BAO). Instead of fitting the measured power spectrum (PS) to a theoretical model containing the cosmological informations and all the nonlinear effects, we define a procedure to project out (or to "extract") the oscillating component from a given nonlinear PS. We show that the BAO scale extracted in this way is extremely robust and, moreover, can be reproduced by simple theoretical models at any redshift. By using N-body simulations, we discuss the effect of the nonlinear evolution of the matter field, of redshift space distortions and of scale-dependent halo bias, showing that all these effects can be reproduced with sub-percent accuracy. We give a one-parameter theoretical model based on a simple (IR) modification of 1-loop perturbation theory, which reproduces the BAO scale from measurements of halo clustering in redshift space at better than $0.1\%$ level and does not need any external UV input, such as coefficients measured from N-body simulations.

We present 108 full-sky gravitational lensing simulation data sets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological $N$-body simulations. The data sets include full-sky convergence and shear maps from redshifts $z=0.05$ to $5.3$ at intervals of $150 \, h^{-1}{\rm Mpc}$ comoving radial distance (corresponding to a redshift interval of $\Delta z \simeq 0.05$ at the nearby universe), enabling the construction of a mock shear catalog for an arbitrary source distribution up to $z=5.3$. The dark matter halos are identified from the same $N$-body simulations with enough mass resolution to resolve the host halos of the Sloan Digital Sky Survey (SDSS) CMASS and Luminous Red Galaxies (LRGs). Angular positions and redshifts of the halos are provided by a ray-tracing calculation, enabling the creation of a mock halo catalog to be used for galaxy-galaxy and cluster-galaxy lensing. The simulation also yields maps of gravitational lensing deflections for a source redshift at the last scattering surface, and we provide 108 realizations of lensed cosmic microwave background (CMB) maps in which the post-Born corrections caused by multiple light scattering are included. We present basic statistics of the simulation data, including the angular power spectra of cosmic shear, CMB temperature and polarization anisotropies, galaxy-galaxy lensing signals for halos, and their covariances. The angular power spectra of the cosmic shear and CMB anisotropies agree with theoretical predictions within $5\%$ up to $\ell = 3000$ (or at an angular scale $\theta > 0.5$ arcsin). The simulation data sets are generated primarily for the ongoing Subaru Hyper Suprime-Cam survey but are freely available for download at http://cosmo.phys.hirosaki-u.ac.jp/takahasi/allsky_raytracing.

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.

We present a statistical study of the redshift evolution of the cluster galaxy population over a wide redshift range from 0.1 to 1.1, using $\sim 1900$ optically-selected CAMIRA clusters from $\sim 232$~deg$^2$ of the Hyper Suprime-Cam (HSC) Wide S16A data. Our stacking technique with a statistical background subtraction reveals color-magnitude diagrams of red-sequence and blue cluster galaxies down to faint magnitudes of $m_z\sim 24$. We find that the linear relation of red-sequence galaxies in the color-magnitude diagram extends down to the faintest magnitudes we explore with a small intrinsic scatter $\sigma_{\rm int}(g-r)<0.1$. The scatter does not evolve significantly with redshift. The stacked color-magnitude diagrams are used to define red and blue galaxies in clusters for studying their radial number density profiles without resorting to photometric redshifts of individual galaxies. We find that red galaxies are significantly more concentrated toward cluster centers and blue galaxies dominate the outskirt of clusters. We explore the fraction of red galaxies in clusters as a function of redshift, and find that the red fraction decreases with increasing distances from cluster centers. The red fraction exhibits a moderate decrease with increasing redshift. The radial number density profiles of cluster member galaxies are also used to infer the location of the steepest slope in the three dimensional galaxy density profiles. For a fixed threshold in richness, we find little redshift evolution in this location.

Non-Gaussianities of dynamical origin are disentangled from primordial ones using the formalism of large deviation statistics with spherical collapse dynamics. This is achieved by relying on accurate analytical predictions for the one-point probability distribution function (PDF) and the two-point clustering of spherically-averaged cosmic densities (sphere bias). Sphere bias extends the idea of halo bias to intermediate density environments and voids as underdense regions. In the presence of primordial non-Gaussianity, sphere bias displays a strong scale dependence relevant for both high and low density regions, which is predicted analytically. The statistics of densities in spheres are built to model primordial non-Gaussianity via an initial skewness with a scale-dependence that depends on the bispectrum of the underlying model. The analytical formulas with the measured nonlinear dark matter variance as input are successfully tested against numerical simulations. For local non-Gaussianity with a range from $f_{\rm NL}=-100$ to $+100$ they are found to agree within 2\% or better for densities $\rho\in[0.5,3]$ in spheres of radius 15 Mpc$/h$ down to $z=0.35$. The validity of the large deviation statistics formalism is thereby established for all observationally relevant local-type departures from perfectly Gaussian initial conditions. The corresponding estimators for the amplitude of the nonlinear variance $\sigma_8$ and primordial skewness $f_{\rm NL}$ are validated using a fiducial joint maximum likelihood experiment. The influence of observational effects and the prospects for a future detection of primordial non-Gaussianity from joint one- and two-point densities-in-spheres statistics are discussed.

We investigate the effects of intrinsic alignments (IA) of dark-matter halo shapes on cosmic density and velocity fields from cluster to cosmic scales beyond 100 Mpc/h. Besides the density correlation function binned by the halo orientation angle which was used in the literature, we introduce, for the first time, the corresponding two velocity statistics, the angle-binned pairwise infall momentum and momentum correlation function. Using large-volume, high-resolution N-body simulations, we measure the alignment statistics of density and velocity, both in real and redshift space. We find that the alignment signal is not amplified by redshift-space distortions at linear scales. Behaviors of IA 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 a scale larger than those with the density field, x>100 Mpc/h. On halo scales, x~ R_{200m} ~ 1 Mpc/h, we detect a sharp steepening in the momentum correlation associated with the physical halo boundary, or the splashback feature, which is found more prominent than in the density correlation. Our results indicate that observations of IA with the velocity field can provide additional information on cosmological models from large scales and on physical sizes of halos from small scales.

We develop a method to simulate galaxy-galaxy weak lensing by utilizing all-sky, light-cone simulations and their inherent halo catalogs. Using the mock catalog to study the error covariance matrix of galaxy-galaxy weak lensing, we compare the full covariance with the "jackknife" (JK) covariance, the method often used in the literature that estimates the covariance from the resamples of the data itself. We show that there exists the variation of JK covariance over realizations of mock lensing measurements, while the average JK covariance over mocks can give a reasonably accurate estimation of the true covariance up to separations comparable with the size of JK subregion. The scatter in JK covariances is found to be $\sim10\%$ after we subtract the lensing measurement around random points. However, the JK method tends to underestimate the covariance at the larger separations, more increasingly for a survey with a higher number density of source galaxies. We apply our method to the the Sloan Digital Sky Survey (SDSS) data, and show that the 48 mock SDSS catalogs nicely reproduce the signals and the JK covariance measured from the real data. We then argue that the use of the accurate covariance, compared to the JK covariance, allows us to use the lensing signals at large scales beyond a size of the JK subregion, which contains cleaner cosmological information in the linear regime.

We investigate the information content of various cosmic shear statistics on the theory of gravity. Focusing on the Hu-Sawicki-type $f(R)$ model, we perform a set of ray-tracing simulations and measure the convergence bispectrum, peak counts and Minkowski functionals. We first show that while the convergence power spectrum does have sensitivity to the current value of extra scalar degree of freedom $|f_{\rm R0}|$, it is largely compensated by a change in the present density amplitude parameter $\sigma_{8}$ and the matter density parameter $\Omega_{\rm m0}$. With accurate covariance matrices obtained from 1000 lensing simulations, we then examine the constraining power of the three additional statistics. We find that these probes are indeed helpful to break the parameter degeneracy, which can not be resolved from the power spectrum alone. We show that especially the peak counts and Minkowski functionals have the potential to rigorously (marginally) detect the signature of modified gravity with the parameter $|f_{\rm R0}|$ as small as $10^{-5}$ ($10^{-6}$) if we can properly model them on small ($\sim 1\, \mathrm{arcmin}$) scale in a future survey with a sky coverage of 1,500 squared degrees. We also show that the signal level is similar among the additional three statistics and all of them provide complementary information to the power spectrum. These findings indicate the importance of combining multiple probes beyond the standard power spectrum analysis to detect possible modifications to General Relativity.

We characterize the non-sphericity of galaxy clusters by the projected axis ratio of spatial distribution of star, dark matter, and X-ray surface brightness (XSB). We select 40 simulated groups and clusters of galaxies with mass larger than 5E13 Msun from the Horizon simulation that fully incorporates the relevant baryon physics, in particular, the AGN feedback. We find that the baryonic physics around the central region of galaxy clusters significantly affects the non-sphericity of dark matter distribution even beyond the central region, approximately up to the half of the virial radius. Therefore it is very difficult to predict the the probability density function (PDF) of the projected axis ratio of XSB from dark-matter only N-body simulations as attempted in previous studies. Indeed we find that the PDF derived from our simulated clusters exhibits much better agreement with that from the observed X-ray clusters. This indicates that our present methodology to estimate the non-sphericity directly from the Horizon simulation is useful and promising. Further improvements in both numerical modeling and observational data will establish the non-sphericity of clusters as a cosmological test complementary to more conventional statistics based on spherically averaged quantities.

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.

We revisit the non-sphericity of cluster-mass scale halos from cosmological N-body simulation on the basis of triaxial modelling. In order to understand the difference between the simulation results and the conventional ellipsoidal collapse model (EC), we first consider the evolution of individual simulated halos. The major difference between EC and the simulation becomes appreciable after the turn-around epoch. Moreover, it is sensitive to the individual evolution history of each halo. Despite such strong dependence on individual halos, the resulting nonsphericity of halos exhibits weak but robust mass dependence in a statistical fashion; massive halos are more spherical up to the turn-around, but gradually become less spherical by z = 0. This is clearly inconsistent with the EC prediction; massive halos are usually more spherical. In addition, at z=0, inner regions of the halos are less spherical than outer regions, i.e., the density distribution inside the halos is highly inhomogeneous and therefore not self-similar. Since most of previous fitting formulae for the PDF of axis ratio of triaxial ellipsoids have been constructed under the self-similarity assumption, they are not accurate. Indeed, we compute the PDF of projected axis ratio a1/a2 directly from the simulation data without the self-similarity assumption, and find that it is very sensitive to the assumption. The latter needs to be carefully taken into account in direct comparison with observations, and therefore we provide an empirical fitting formula for the PDF of a1/a2. Our preliminary analysis suggests that the derived PDF of a1/a2 roughly agrees with the current weak-lensing observations. More importantly, the present results will be useful in future exploration of the non-sphericity of clusters in X-ray and optical observations.

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.

Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations may ease the assembly of massive BHs. In particular, they produce more low-mass halos at high redshift, thus altering the production of metals and ultra-violet flux, believed to be important factors in BH formation. Additionally, a higher number of progenitors and of nearly equal-mass halo mergers would boost the mass increase provided by BH-BH mergers and merger-driven accretion. We use a set of two cosmological simulations, with either Gaussian or scale-dependent non-Gaussian primordial fluctuations to perform a proof-of-concept experiment to estimate how BH formation and growth are altered. We estimate the BH number density and the fraction of halos where BHs form, for both simulations and for two popular scenarios of BH formation (remnants of the first generation of stars and direct collapse in the absence of metals and molecular hydrogen). We find that the fractions of halos where BHs form are almost identical, but that non-Gaussian primordial perturbations increase the total number density of BHs for the both BH formation scenarios. We also evolve BHs using merger trees extracted from the simulations and find that non-Gaussianities increase both the BH masses and the number of the most massive BHs.

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.

We present the redshift-space generalization of the equal-time angular-averaged consistency relations between $(\ell+n)$- and $n$-point polyspectra of the cosmological matter density field. Focusing on the case of $\ell=1$ large-scale mode and $n$ small-scale modes, we use an approximate symmetry of the gravitational dynamics to derive explicit expressions that hold beyond the perturbative regime, including both the large-scale Kaiser effect and the small-scale fingers-of-god effects. We explicitly check these relations, both perturbatively, for the lowest-order version that applies to the bispectrum, and nonperturbatively, for all orders but for the one-dimensional dynamics. Using a large ensemble of $N$-body simulations, we find that our squeezed bispectrum relation is valid to better than $20\%$ up to $1h$Mpc$^{-1}$, for both the monopole and quadrupole at $z=0.35$, in a $\Lambda$CDM cosmology. Additional simulations done for the Einstein-de Sitter background suggest that these discrepancies mainly come from the breakdown of the approximate symmetry of the gravitational dynamics. For practical applications, we introduce a simple ansatz to estimate the new derivative terms in the relation using only observables. Although the relation holds worse after using this ansatz, we can still recover it within $20\%$ up to $1h$Mpc$^{-1}$, at $z=0.35$ for the monopole. On larger scales, $k = 0.2 h\mathrm{Mpc}^{-1}$, it still holds within the statistical accuracy of idealized simulations of volume $\sim8h^{-3}\mathrm{Gpc}^3$ without shot-noise error.

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.

[Abridged] Understanding how the intergalactic medium (IGM) was reionized at z > 6 is one of the big challenges of current high redshift astronomy. It requires modelling the collapse of the first astrophysical objects (Pop III stars, first galaxies) and their interaction with the IGM, while at the same time pushing current observational facilities to their limits. The observational and theoretical progress of the last few years have led to the emergence of a coherent picture in which the budget of hydrogen-ionizing photons is dominated by low-mass star-forming galaxies, with little contribution from Pop III stars and quasars. The reionization history of the Universe therefore critically depends on the number density of low-mass galaxies at high redshift. In this work, we explore how changes in the statistical properties of initial density fluctuations affect the formation of early galaxies. Following Habouzit et al. (2014), we run 5 N-body simulations with Gaussian and (scale-dependent) non-Gaussian initial conditions, all consistent with Planck constraints. By appealing to a galaxy formation model and to a population synthesis code, we compute the far-UV galaxy luminosity function down to M_UV = -14 at redshift 7 < z < 15. We find that models with strong primordial non-Gaussianities on < Mpc scales show a far-UV luminosity function significantly enhanced in low-mass galaxies. We adopt a reionization model calibrated from state-of-the-art hydrodynamical simulations and show that such non-Gaussianities leave a clear imprint on the Universe reionization history and electron Thomson scattering optical depth tau_E. Although current uncertainties in the physics of reionization and on the determination of tau_E still dominate the signatures of non-Gaussianities, our results suggest that tau_E could ultimately be used to constrain the statistical properties of initial density fluctuations.

Primordial non-Gaussianities provide an important test of inflationary models. Although the Planck CMB experiment has produced strong limits on non-Gaussianity on scales of clusters, there is still room for considerable non-Gaussianity on galactic scales. We have tested the effect of local non-Gaussianity on the high redshift galaxy population by running five cosmological N-body simulations down to z=6.5. For these simulations, we adopt the same initial phases, and either Gaussian or scale-dependent non-Gaussian primordial fluctuations, all consistent with the constraints set by Planck on clusters scales. We then assign stellar masses to each halo using the halo - stellar mass empirical relation of Behroozi et al. (2013). Our simulations with non-Gaussian initial conditions produce halo mass functions that show clear departures from those obtained from the analogous simulations with Gaussian initial conditions at z>~10. We observe a >0.3 dex enhancement of the low-end of the halo mass function, which leads to a similar effect on the galaxy stellar mass function, which should be testable with future galaxy surveys at z>10. As cosmic reionization is thought to be driven by dwarf galaxies at high redshift, our findings may have implications for the reionization history of the Universe.

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 explicitly test the equal-time consistency relation between the angular-averaged bispectrum and the power spectrum of the matter density field, employing a large suite of cosmological $N$-body simulations. This is the lowest-order version of the relations between $(\ell+n)$-point and $n$-point polyspectra, where one averages over the angles of $\ell$ soft modes. This relation depends on two wave numbers, $k'$ in the soft domain and $k$ in the hard domain. We show that it holds up to a good accuracy, when $k'/k\ll 1$ and $k'$ is in the linear regime, while the hard mode $k$ goes from linear ($0.1\,h\mathrm{Mpc}^{-1}$) to nonlinear ($1.0\,h\mathrm{Mpc}^{-1}$) scales. On scales $k\lesssim 0.4\,h\mathrm{Mpc}^{-1}$, we confirm the relation within the statistical error of the simulations (typically a few percent depending on the wave number), even though the bispectrum can already deviate from leading-order perturbation theory by more than $30\%$. We further examine the relation on smaller scales with higher resolution simulations. We find that the relation holds within the statistical error of the simulations at $z=1$, whereas we find deviations as large as $\sim 7\%$ at $k \sim 1.0\,h\mathrm{Mpc}^{-1}$ at $z=0.35$. We show that this can be explained partly by the breakdown of the approximation $\Omega_\mathrm{m}/f^2\simeq1$ with supplemental simulations done in the Einstein-de Sitter background cosmology. We also estimate the impact of this approximation on the power spectrum and bispectrum.

We model the apparent clustering anisotropy of Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey using subhalos identified in cosmological $N$-body simulations. We first conduct a Markov-chain Monte Carlo analysis on the parameters characterizing subhalos hosting LRGs assuming a specific $\Lambda$CDM cosmology on which we run the simulations. We show that simple models with central and satellite subhalos can explain the observed multipole moments of the power spectrum up to hexadecapole on large scales ($k\lesssim0.3~h\mathrm{Mpc}^{-1}$). A satellite fraction of $20$ to $30$ per cent is favored weakly depending on the detail of the model. The fraction is shown to be robust when we adopt a more refined model based on the halo occupation number from the literature. We then vary cosmological parameters controlling the anisotropy in redshift-space effectively by deforming the simulation box (the Alcock-Paczynski effect) and changing the amplitude of the velocities (the redshift-space distortions). We demonstrate that we can constrain the geometry of the universe, the structure growth rate, and the parameters characterizing LRGs simultaneously. This is a step toward cosmological analysis with realistic bias description beyond empirical bias functions with nuisance parameters.

Redshift space distortion (RSD) observed in galaxy redshift surveys is a powerful tool to test gravity theories on cosmological scales, but the systematic uncertainties must carefully be examined for future surveys with large statistics. Here we employ various analytic models of RSD and estimate the systematic errors on measurements of the structure growth-rate parameter, $f\sigma_8$, induced by non-linear effects and the halo bias with respect to the dark matter distribution, by using halo catalogues from 40 realisations of $3.4 \times 10^8$ comoving $h^{-3}$Mpc$^3$ cosmological N-body simulations. We consider hypothetical redshift surveys at redshifts z=0.5, 1.35 and 2, and different minimum halo mass thresholds in the range of $5.0 \times 10^{11}$ -- $2.0 \times 10^{13} h^{-1} M_\odot$. We find that the systematic error of $f\sigma_8$ is greatly reduced to ~5 per cent level, when a recently proposed analytical formula of RSD that takes into account the higher-order coupling between the density and velocity fields is adopted, with a scale-dependent parametric bias model. Dependence of the systematic error on the halo mass, the redshift, and the maximum wavenumber used in the analysis is discussed. We also find that the Wilson-Hilferty transformation is useful to improve the accuracy of likelihood analysis when only a small number of modes are available in power spectrum measurements.

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.

We study how well the Gaussian approximation is valid for computing the covariance matrices of the convergence power and bispectrum in weak gravitational lensing analyses. We focus on its impact on the cosmological parameter estimations by comparing the results with and without non-Gaussian error contribution in the covariance matrix. We numerically derive the covariance matrix as well as the cosmology dependence of the spectra from a large set of N-body simulations performed for various cosmologies and carry out Fisher matrix forecasts for tomographic weak lensing surveys with three source redshifts. After showing the consistency of the power and bispectra measured from our simulations with the state-of-the-art fitting formulas, we investigate the covariance matrix assuming a typical ongoing survey across 1500 deg^2 with the mean source number density of 30 arcmin^{-2} at the mean redshift z_s=1.0. Although the shape noise contributes a significant fraction to the total error budget and it mitigates the impact of the non-Gaussian error for this source number density, we find that the non-Gaussian error degrades the cumulative signal-to-noise ratio up to the maximum multipole of 2000 by a factor of about 2 (3) in the power (bi-) spectrum analysis. Its impact on the final cosmological parameter forecast with 6 parameters can be as large as 15% in the size of the one-dimensional statistical error. This can be a problem in future wide and deep weak lensing surveys for precision cosmology. We also show how much the dark energy figure of merit is affected by the non-Gaussian error contribution and demonstrate an optimal survey design with a fixed observational time.

The comparison between dynamical mass and lensing mass provides a targeted test for a wide range of modified gravity models. In our previous paper we showed, through numerical simulations, that the measurement of the line-of-sight velocity dispersion around stacked massive clusters whose lensing masses are known allows for stringent constraints on modified gravity on scales of 2 - 15 Mpc/h. In this work we develop a semi-analytical approach based on the halo model to describe the phase-space distribution and the line-of-sight velocity dispersion for different tracers. The model distinguishes contributions from the halo pairwise velocity and the virial velocity within halos. We also discuss observational complications, in particular the contribution from Hubble flow, and show how our model can incorporate these complications. We then incorporate the effects of modified gravity (specifically, f(R) and braneworld models), and show that the model predictions are in excellent agreement with modified gravity simulations. More broadly, the phase-space distribution provides a sensitive test of our understanding of hierarchical structure formation when confronted with observations via this model.

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.

We revisit an analytical model to describe the halo-matter cross-power spectrum and the halo auto-power spectrum in the weakly nonlinear regime, by combining the perturbation theory (PT) for matter clustering, the local bias model, and the halo bias. Nonlinearities in the power spectra arise from the nonlinear clustering of matter as well as the nonlinear relation between the matter and halo density fields. By using the "renormalization" approach, we express the nonlinear power spectra by a sum of the two contributions: the nonlinear matter power spectrum with the effective linear bias parameter, and the higher-order PT spectra having the halo bias parameters as the coefficients. The halo auto-power spectrum includes the residual shot noise contamination that needs to be treated as additional free parameter. The term(s) of the higher-order PT spectra and the residual shot noise cause a scale-dependent bias function relative to the nonlinear matter power spectrum in the weakly nonlinear regime. We show that the model predictions are in good agreement with the spectra measured from a suit of high-resolution $N$-body simulations up to $k\simeq 0.2 h$/Mpc at $z=0.35$, for different halo mass bins.

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 cosmological information contained in the Minkowski Functionals (MFs) of weak gravitational lensing convergence maps. We show that the MFs provide strong constraints on the local type primordial non-Gaussianity parameter f_NL. We run a set of cosmological N-body simulations and perform ray-tracing simulations of weak lensing, to generate 100 independent convergence maps of 25 deg^2 field-of-view for f_NL = -100, 0 and 100. We perform a Fisher analysis to study the degeneracy among other cosmological parameters such as the dark energy equation of state parameter w and the fluctuation amplitude sigma_8. We use fully nonlinear covariance matrices evaluated from 1000 ray-tracing simulations. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam survey with the proposed survey area of 1500 deg^2, the primordial non-Gaussianity can be constrained with a level of f_NL ~ 80 and w ~ 0.036 by weak lensing MFs. If simply scaled by the effective survey area, a 20000 deg^2 lensing survey using Large Synoptic Survey Telescope will give constraints of f_NL ~ 25 and w ~ 0.013. We show that these constraints can be further improved by a tomographic method using source galaxies in multiple redshift bins.

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.

In General Relativity, the average velocity field of dark matter around galaxy clusters is uniquely determined by the mass profile. The latter can be measured through weak lensing. We propose a new method of measuring the velocity field (phase space density) by stacking redshifts of surrounding galaxies from a spectroscopic sample. In combination with lensing, this yields a direct test of gravity on scales of 1-30 Mpc. Using N-body simulations, we show that this method can improve upon current constraints on f(R) and DGP model parameters by several orders of magnitude when applied to upcoming imaging and redshift surveys.

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing statistics. We compare our predictions for the weak-lensing convergence power spectrum and bispectrum with numerical simulations and fitting formulas proposed in previous works. We find that this model provides better agreement with simulations than published fitting formulas. This shows that building on systematic and physically motivated models is a promising approach. Moreover, this makes explicit the link between the weak-lensing statistics and the underlying properties of the 3D matter distribution, as a function of scale $\ell$. Thus, we obtain the contributions to the lensing power spectrum and bispectrum that arise from perturbative terms (complete up to one-loop) and nonperturbative terms (e.g., "1-halo" term). Finally, we show that this approach recovers the dependence on cosmology (for realistic scenarios).

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing configuration-space statistics. We compared our predictions for the weak-lensing convergence two-point and three-point correlation functions with numerical simulations and fitting formulas proposed in previous works. We also considered the second- and third-order moments of the smoothed convergence and of the aperture-mass. As in our previous study of Fourier-space weak-lensing statistics, we find that our model agrees better with simulations than previously published fitting formulas. Moreover, we recover the dependence on cosmology of these weak-lensing statistics and we can describe multi-scale moments. This approach allows us to obtain the quantitative relationship between these integrated weak-lensing statistics and the various contributions to the underlying 3D density fluctuations, decomposed over perturbative, two-halo, or one-halo terms.

We investigate the clustering of halos in cosmological models starting with general local-type non-Gaussian primordial fluctuations. We employ multiple Gaussian fields and add local-type non-Gaussian corrections at arbitrary order to cover a class of models described by frequently-discussed f_nl, g_nl and \tau_nl parameterization. We derive a general formula for the halo power spectrum based on the peak-background split formalism. The resultant spectrum is characterized by only two parameters responsible for the scale-dependent bias at large scale arising from the primordial non-Gaussianities in addition to the Gaussian bias factor. We introduce a new inequality for testing non-Gaussianities originating from multi fields, which is directly accessible from the observed power spectrum. We show that this inequality is a generalization of the Suyama-Yamaguchi inequality between f_nl and \tau_nl to the primordial non-Gaussianities at arbitrary order. We also show that the amplitude of the scale-dependent bias is useful to distinguish the simplest quadratic non-Gaussianities (i.e., f_nl-type) from higher-order ones (g_nl and higher), if one measures it from multiple species of galaxies or clusters of galaxies. We discuss the validity and limitations of our analytic results by comparison with numerical simulations in an accompanying paper.

We investigated the spatial correlation between the arrival directions of the highest energy cosmic rays (HECRs) detected by the Akeno Giant Air Shower Array (AGASA) with energies above $4 \times 10^{19}$ eV and the positions of galaxies observed by the Sloan Digital Sky Survey (SDSS) within $z = 0.024$. We systematically tested the dependence of the correlation on the redshift ranges and properties of the galaxies, i.e., absolute luminosity, color, and morphology, to understand where HECR sources are and what objects are HECR sources. In the systematic survey, we found potential signals of the positive correlation at small angular scale ($<10^{\circ}$) with the (non-penalized) chance probability less than 5% in intermediate redshift ranges. Then, we estimated penalized probabilities to compensate the trial effects of angular scan, and found that the strongest correlation is produced by early-type galaxies in $0.012 \leq z < 0.018$ at 90% C.L. The possible origin of HECRs which correlating galaxies imply is also discussed.

We investigate how unified models should be built to be able to predict the matter-density bispectrum (and power spectrum) from very large to small scales and that are at the same time consistent with perturbation theory at low $k$ and with halo models at high $k$. We use a Lagrangian framework to decompose the bispectrum into "3-halo", "2-halo", and "1-halo" contributions, related to "perturbative" and "non-perturbative" terms. We describe a simple implementation of this approach and present a detailed comparison with numerical simulations. We show that the 1-halo and 2-halo contributions contain counterterms that ensure their decay at low $k$, as required by physical constraints, and allow a better match to simulations. Contrary to the power spectrum, the standard 1-loop perturbation theory can be used for the perturbative 3-halo contribution because it does not grow too fast at high $k$. Moreover, it is much simpler and more accurate than two resummation schemes investigated in this paper. We obtain a good agreement with numerical simulations on both large and small scales, but the transition scales are poorly described by the simplest implementation. This cannot be amended by simple modifications to the halo parameters, but we show how it can be corrected for the power spectrum and the bispectrum through a simple interpolation scheme that is restricted to this intermediate regime. Then, we reach an accuracy on the order of 10% on mildly and highly nonlinear scales, while an accuracy on the order of 1% is obtained on larger weakly nonlinear scales. This also holds for the real-space two-point correlation function.

We investigate the building of unified models that can predict the matter-density power spectrum and the two-point correlation function from very large to small scales, being consistent with perturbation theory at low $k$ and with halo models at high $k$. We use a Lagrangian framework to re-interpret the halo model and to decompose the power spectrum into "2-halo" and "1-halo" contributions, related to "perturbative" and "non-perturbative" terms. We describe a simple implementation of this model and present a detailed comparison with numerical simulations, from $k \sim 0.02$ up to $100 h$Mpc$^{-1}$, and from $x \sim 0.02$ up to $150 h^{-1}$Mpc. We show that the 1-halo contribution contains a counterterm that ensures a $k^2$ tail at low $k$ and is important not to spoil the predictions on the scales probed by baryon acoustic oscillations, $k \sim 0.02$ to $0.3 h$Mpc$^{-1}$. On the other hand, we show that standard perturbation theory is inadequate for the 2-halo contribution, because higher order terms grow too fast at high $k$, so that resummation schemes must be used. We describe a simple implementation, based on a 1-loop "direct steepest-descent" resummation for the 2-halo contribution that allows fast numerical computations, and we check that we obtain a good match to simulations at low and high $k$. Our simple implementation already fares better than standard 1-loop perturbation theory on large scales and simple fits to the power spectrum at high $k$, with a typical accuracy of 1% on large scales and 10% on small scales. We obtain similar results for the two-point correlation function. However, there remains room for improvement on the transition scale between the 2-halo and 1-halo contributions, which may be the most difficult regime to describe.

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.

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 how primordial non-Gaussianity affects the pairwise velocity probability density function (PDF) using an analytical model and cosmological N-body simulations. We adopt the local type non-Gaussian models characterized by f_{nl}, and examine both the linear velocity difference PDF and the linear pairwise velocity PDF. We show explicitly how f_{nl} induces correlations between originally independent velocities along the parallel and the perpendicular to the line of separation directions. We compare the model results with measurements from N-body simulations of the non-Gaussian models. Linear theory fails to predict the PDF in the f_{nl} models. Therefore we develop an analytic model based on the Zeldovich approximation to describe the evolution of the velocity PDF. Our analytical model and simulation results show remarkably good agreement in both the parallel and the perpendicular directions for the PDF profiles, as well as the change in the PDF due to primordial non-Gaussianity. The agreement is particularly good for relatively small separations (< 10 h^{-1} Mpc). The inclusion of the evolution of the velocity PDF is important to obtain a good description on the signature of primordial non-Gaussianity in the PDF. Our model provides the foundation to constrain f_{nl} using the peculiar velocity in future surveys.

We report on a new merging group of galaxies, Suzaku J1552+2739 at z ~ 0.08 revealed by a SUZAKU observation. The group was found by observing a junction of galaxy filaments optically identified in the Sloan Digital Sky Survey spectroscopic data. Suzaku J1552+2739 exhibits an irregular morphology and presents several peaks in its X-ray image. A bright elliptical galaxy, observable in the central peak, allows the localization of the group at z = 0.083. We found a significant hot spot visible in the X-ray hardness map, close to the second peak. The spectroscopic temperature is T = 1.6+0.4-0.1 keV within R500 = 0.6 Mpc and T = 3 - 5 keV in the hot spot. We interpret those results as Suzaku J1552+2739 being located in the center of a major merging process. The observation of a galaxy group showing multiple X-ray peaks and a hot spot at the same time is rare and we believe in particular that the study of Suzaku J1552+2739 potentially presents a significant interest to better understand the dynamical and thermal evolution of the intragroup and intracluster medium, as well as its relation with surrounding environment.

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

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.

We investigate correlation between the arrival directions of ultra-high-energy cosmic rays (UHECRs) and the large-scale structure (LSS) of the Universe by using statistical quantities which can find the angular scale of the correlation. The Infrared Astronomical Satellite Point Source Redshift Survey (IRAS PSCz) catalog of galaxies is adopted for LSS. We find a positive correlation of the highest energy events detected by the Pierre Auger Observatory (PAO) with the IRAS galaxies inside $z=0.018$ within the angular scale of $\sim 15^{\circ}$. This positive correlation observed in the southern sky implies that a significant fraction of the highest energy events comes from nearby extragalactic objects. We also analyze the data of the Akeno Giant Air Shower Array (AGASA) which observed the northern hemisphere, but the obvious signals of positive correlation with the galaxy distribution are not found. Since the exposure of the AGASA is smaller than the PAO, the cross-correlation in the northern sky should be tested using a larger number of events detected in the future. We also discuss the correlation using the all-sky combined data sets of both the PAO and AGASA, and find a significant correlation within $\sim 8^{\circ}$. These angular scales can constrain several models of intergalactic magnetic field. These cross-correlation signals can be well reproduced by a source model in which the distribution of UHECR sources is related to the IRAS galaxies.

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.

We investigate the damping of the baryon acoustic oscillations in the matter power spectrum due to the quasinonlinear clustering and redshift-space distortions by confronting the models with the observations of the Sloan Digital Sky Survey luminous red galaxy sample. The chi-squared test suggests that the observed power spectrum is better matched by models with the damping of the baryon acoustic oscillations rather than the ones without the damping.

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 investigate the damping of the baryon acoustic oscillations (BAO) signature in the matter power spectrum due to the quasi-nonlinear clustering of density perturbations. On the basis of the third order perturbation theory, we construct a fitting formula of the damping in an analytic way. This demonstrates that the damping is closely related with the growth factor and the amplitude of the matter power spectrum. Then, we investigate the feasibility of constraining the growth factor through a measurement of the damping of the BAO signature. An extension of our formula including higher order corrections of density perturbations is also discussed.

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.

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