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.

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.

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.

Cross-correlation analysis of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) provides a powerful probe of cosmology and astrophysics of the intra-cluster medium. We present the measurement of the cross-correlation of tSZ and WL from Planck and Subaru Hyper-Suprime Cam. The combination enables us to study cluster astrophysics at high redshift. We use the tSZ-WL cross-correlation and the tSZ auto-power spectrum measurements to place a tight constraint on the hydrostatic mass bias, which is a measure of the degree of non-thermal pressure support in galaxy clusters. With the prior on cosmological parameters derived from the analysis of the cosmic microwave background anisotropies by Planck and taking into account foreground contributions both in the tSZ auto-power spectrum and the tSZ-WL cross-correlation, the hydrostatic mass bias is estimated to be $26.9^{+8.9}_{-4.4} \%$ ($68\%$ C.L.), which is consistent with recent measurements by mass calibration techniques.

We present a novel statistic to extract cosmological information in weak lensing data: the lensing minima. We also investigate the effect of baryons on the cosmological constraints from peak and minimum counts. Using the \texttt{MassiveNuS} simulations, we find that lensing minima are sensitive to non-Gaussian cosmological information and are complementary to the lensing power spectrum and peak counts. For an LSST-like survey, we obtain $95\%$ credible intervals from a combination of lensing minima and peaks that are significantly stronger than from the power spectrum alone, by $44\%$, $11\%$, and $63\%$ for the neutrino mass sum $\sum m_\nu$, matter density $\Omega_m$, and amplitude of fluctuation $A_s$, respectively. We explore the effect of baryonic processes on lensing minima and peaks using the hydrodynamical simulations \texttt{BAHAMAS} and \texttt{Osato15}. We find that ignoring baryonic effects would lead to strong ($\approx 4 \sigma$) biases in inferences from peak counts, but negligible ($\approx 0.5 \sigma$) for minimum counts, suggesting lensing minima are a potentially more robust tool against baryonic effects. Finally, we demonstrate that the biases can in principle be mitigated without significantly degrading cosmological constraints when we model and marginalize the baryonic effects.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

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.

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.

Gravitational waves from inspiraling compact objects provide us with information of the distance scale since we can infer the absolute luminosity of the source from analysis of the wave form, which is known as standard sirens. The first detection of the gravitational wave signal of the binary black hole merger event by Advanced LIGO has opened up the possibility of utilizing standard sirens as cosmological probe. In order to extract information of the distance-redshift relation, we cross-correlate weak lensing, which is an unbiased tracer of matter distribution in the Universe, with the projected number density of gravitational wave sources. For weak lensing, we employ tomography technique to efficiently obtain information of large-scale structures at wide ranges of redshifts. Making use of the cross-correlations along with the auto-correlations, we present forecast of constraints on four cosmological parameters, i.e., Hubble parameter, matter density, the equation of state parameter of dark energy, and the amplitude of matter fluctuation. To fully explore the ability of cross-correlations, which require large overlapping sky coverage, we consider the specific case with the upcoming surveys by \textit{Euclid} for weak lensing and Einstein Telescope for standard sirens. We show that cosmological parameters can be tightly constrained solely by these auto- and cross-correlations of standard sirens and weak lensing. For example, the $1\text{-}\sigma$ error of Hubble parameter is expected to be $\sigma (H_0) = 0.33 \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$. Thus, the proposed statistics will be a promising probe into the distance scale.

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.

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.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

Recent detections of the cross-correlation of the thermal Sunyaev-Zel'dovich (tSZ) effect and weak gravitational lensing (WL) enable unique studies of cluster astrophysics and cosmology. In this work, we present constraints on the amplitude of the non-thermal pressure fraction in galaxy clusters, $\alpha_0$, and the amplitude of the matter power spectrum, $\sigma_8$, using measurements of the tSZ power spectrum from Planck, and the tSZ-WL cross-correlation from Planck and the Red Cluster Sequence Lensing Survey. We fit the data to a semi-analytic model with the covariance matrix using $N$-body simulations. We find that the tSZ power spectrum alone prefers $\sigma_8 \sim 0.85$ and a large fraction of non-thermal pressure ($\alpha_0 \sim 0.2$-$0.3$). The tSZ-WL cross-correlation on the other hand prefers a significantly lower $\sigma_8 \sim 0.6$, and low $\alpha_0 \sim 0.05$. We show that this tension can be mitigated by allowing for a steep slope in the stellar-mass-halo-mass relation, which would cause a reduction of the gas in low-mass halos. In such a model, the combined data prefer $\sigma_8 \sim 0.7$ and $\alpha_0 \sim 0.2$, consistent with predictions from hydrodynamical simulations.

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

Light gravitinos of mass $\lesssim \mathcal{O} (10)$ eV are of particular interest in cosmology, offering various baryogenesis scenarios without suffering from the cosmological gravitino problem. The gravitino may contribute considerably to the total matter content of the Universe and affect structure formation from early to present epochs. After the gravitinos decouple from other particles in the early Universe, they free-stream and consequently suppress density fluctuations of (sub-)galactic length scales. Observations of structure at the relevant length-scales can be used to infer or constrain the mass and the abundance of light gravitinos. We derive constraints on the light gravitino mass using the data of cosmic microwave background (CMB) lensing from Planck and of cosmic shear from the Canada France Hawaii Lensing Survey, combined with analyses of the primary CMB anisotropies and the signature of baryon acoustic oscillations in galaxy distributions. The obtained constraint on the gravitino mass is $m_{3/2} < 4.7$ eV (95% C.L.), which is substantially tighter than the previous constraint from clustering analysis of Ly-$\alpha$ forests.

The top-hat spherical collapse model (TSC) is one of the most fundamental analytical frameworks to describe the non-linear growth of cosmic structure. TSC has motivated, and been widely applied in, various researches even in the current era of precision cosmology. While numerous studies exist to examine its validity against numerical simulations in a statistical fashion, there are few analyses to compare the TSC dynamics in an individual object-wise basis, which is what we attempt in the present paper. We extract 100 halos at z = 0 from a cosmological N-body simulation according to the conventional TSC criterion for the spherical over-density. Then we trace back their spherical counter-parts at earlier epochs. Just prior to the turn-around epoch of the halos, their dynamics is well approximated by TSC, but their turn-around epochs are systematically delayed and the virial radii are larger by ~ 20 percent on average relative to the TSC predictions. We find that this systematic deviation is mainly ascribed to the non-uniformity/inhomogeneity of dark matter density profiles and the non-zero velocity dispersions, both of which are neglected in TSC. In particular, the inside-out-collapse and shell-crossing of dark matter halos play an important role in generating the significant velocity dispersion. The implications of the present result are briefly discussed.

We study the impact of baryonic physics on cosmological parameter estimation with weak lensing surveys. We run a set of cosmological hydrodynamics simulations with different galaxy formation models. We then perform ray-tracing simulations through the total matter density field to generate 100 independent convergence maps of 25 ${\rm deg}^2$ field-of-view, and use them to examine the ability of the following three lensing statistics as cosmological probes; power spectrum, peak counts, and Minkowski functionals. For the upcoming wide-field observations such as Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of 1400 ${\rm deg}^2$, these three statistics provide tight constraints on the matter density, density fluctuation amplitude, and dark energy equation of state, but parameter bias is induced by the baryonic processes such as gas cooling and stellar feedback. The bias can be significant when the statistical errors become small in future observations with a much larger survey area. We find the bias is induced in different directions in the parameter space depending on the statistics employed. While the two-point statistic, i.e. power spectrum, yields robust results against baryonic effects, the overall constraining power is weak compared with peak counts and Minkowski functionals. On the other hand, using one of peak counts or Minkowski functionals, or combined analysis with multiple statistics, results in biased parameter estimate. The bias can be as large as $1\sigma$ for HSC survey, and will be more significant for upcoming wider area surveys. We suggest to use an optimized combination so that the baryonic effects on parameter estimation are mitigated. Such `calibrated' combination can place stringent and robust constraints on cosmological parameters.