We present a weak lensing cluster search using Hyper Suprime-Cam Subaru Strategic Program (HSC survey) first-year data. We pay special attention to the dilution effect of cluster member and foreground galaxies on weak lensing signals from clusters of galaxies; we adopt the globally normalized weak lensing estimator which is least affected by cluster member galaxies, and we select source galaxies by using photometric redshift information to mitigate the effect of foreground galaxies. We produce six samples of source galaxies with different low-$z$ galaxy cuts, construct weak lensing mass maps for each of the source sample, and search for high peaks in the mass maps that cover an effective survey area of $\sim$120 deg$^2$. We combine six catalogs of high peaks into a sample of cluster candidates which contains 124 high peaks with signal-to-noise ratios greater than five. We cross-match the peak sample with the public optical cluster catalog constructed from the same HSC survey data to identify cluster counterparts of the peaks. We find that 107 out of 124 peaks have matched clusters within 5 arcmin from peak positions. Among them, we define a sub-sample of 64 secure clusters that we use to examine dilution effects on our weak lensing cluster search. We find that source samples with the low-$z$ galaxy cuts mitigate the dilution effect on weak lensing signals of high-$z$ clusters ($z > 0.3$), and thus combining multiple peak catalogs from different source samples improves the efficiency of weak lensing cluster searches.

Line intensity mapping (LIM) is an emerging observational method to study the large-scale structure of the Universe and its evolution. LIM does not resolve individual sources but probes the fluctuations of integrated line emissions. A serious limitation with LIM is that contributions of different emission lines from sources at different redshifts are all confused at an observed wavelength. We propose a deep learning application to solve this problem. We use conditional generative adversarial networks to extract designated information from LIM. We consider a simple case with two populations of emission line galaxies; H$\rm\alpha$ emitting galaxies at $z = 1.3$ are confused with [OIII] emitters at $z = 2.0$ in a single observed waveband at 1.5 $\rm\mu$m. Our networks trained with 30,000 mock observation maps are able to extract the total intensity and the spatial distribution of H$\rm\alpha$ emitting galaxies at $z = 1.3$. The intensity peaks are successfully located with 74% precision. The precision increases to 91% when we combine the results of 5 networks. The mean intensity and the power spectrum are reconstructed with an accuracy of $\sim$10%. The extracted galaxy distributions at a wider range of redshift can be used for studies on cosmology and on galaxy formation and evolution.

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.

Cosmic rays in galaxy clusters are unique probes of energetic processes operating with large-scale structures in the Universe. Precise measurements of cosmic rays in galaxy clusters are essential for improving our understanding of non-thermal components in the intracluster-medium (ICM) as well as the accuracy of cluster mass estimates in cosmological analyses. In this paper, we perform a cross-correlation analysis with the extragalactic gamma-ray background and the thermal Sunyaev-Zel'dovich (tSZ) effect in the cosmic microwave background. The expected cross-correlation signal would contain rich information about the cosmic-ray-induced gamma-ray emission in the most massive galaxy clusters at $z\sim0.1-0.2$. We analyze the gamma-ray background map with 8 years of data taken by the Large Area Telescope onboard Fermi satellite and the publicly available tSZ map by Planck. We confirm that the measured cross-correlation is consistent with a null detection, and thus it enables us to put the tightest constraint on the acceleration efficiency of cosmic ray protons at shocks in and around galaxy clusters. We find the acceleration efficiency must be below 5\% with a $2\sigma$ confidence level when the hydrostatic mass bias of clusters is assumed to be 30\%, and our result is not significantly affected by the assumed value of the hydrostatic mass bias. Our constraint implies that the non-thermal cosmic-ray pressure in the ICM can introduce only a $\le 3\%$ level of the hydrostatic mass bias, highlighting that cosmic rays alone do not account for the mass bias inferred by the Planck analyses. Finally, we discuss future detectability prospects of cosmic-ray-induced gamma rays from the Perseus cluster for the Cherenkov Telescope Array.

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.

Searches for dark matter annihilation signals have been carried out in a number of target regions such as the Galactic Center and Milky Way dwarf spheroidal galaxies (dSphs), among a few others. Here we propose low surface brightness galaxies (LSBGs) asnovel targets for the indirect detection of dark matter emission. In particular, LSBGs are known to have very large dark matter contents and be less contaminated by extragalactic gamma-ray sources (e.g., blazars) compared to star forming galaxies. We report on an analysis that uses eight LSBGs (detected by Subaru Hyper Suprime-Cam survey data) with known redshifts to conduct a search for gamma-ray emission at the positions of these new objects in Fermi Large Area Telescope data. We found no excesses of gamma-ray emission and set constraints on the dark matter annihilation cross-section. We exclude (at the 95% C.L.) dark matter scenarios predicting a cross-section higher than 10^-23[cm^3/s] for dark matter particles of mass 10 GeV self-annihilating in the b_b channel. Although this constraint is weaker than the ones reported in recent studies using other targets, we note that in the near future, the number of detections of new LSBGs will increase by a few orders of magnitude. We forecast that with the use of the full catalog of soon-to-be-detected LSBGs the constraint will reach cross-section sensitivities of ~ 3*10^-25 [cm^3/s] for dark matter particles with masses less than 10 GeV.

We present measurements of cosmic shear two-point correlation functions (TPCFs) from Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) first-year data, and derived cosmological constraints based on a blind analysis. The HSC first-year shape catalog is divided into four tomographic redshift bins ranging from $z=0.3$ to 1.5 with equal widths of $\Delta z =0.3$. The unweighted galaxy number densities in each tomographic bin are 5.9, 5.9, 4.3, and 2.4 arcmin$^{-2}$ from lower to higher redshifts, respectively. We adopt the standard TPCF estimators, $\xi_\pm$, for our cosmological analysis, given that we find no evidence of the significant B-mode shear. The TPCFs are detected at high significance for all ten combinations of auto- and cross-tomographic bins over a wide angular range, yielding a total signal-to-noise ratio of 19 in the angular ranges adopted in the cosmological analysis, $7'<\theta<56'$ for $\xi_+$ and $28'<\theta<178'$ for $\xi_-$. We perform the standard Bayesian likelihood analysis for cosmological inference from the measured cosmic shear TPCFs, including contributions from intrinsic alignment of galaxies as well as systematic effects from PSF model errors, shear calibration uncertainty, and source redshift distribution errors. We adopt a covariance matrix derived from realistic mock catalogs constructed from full-sky gravitational lensing simulations that fully account for survey geometry and measurement noise. For a flat $\Lambda$ cold dark matter model, we find $S_8 \equiv \sigma_8\sqrt{\Omega_m/0.3}=0.804_{-0.029}^{+0.032}$, and $\Omega_m=0.346_{-0.100}^{+0.052}$. We carefully check the robustness of the cosmological results against astrophysical modeling uncertainties and systematic uncertainties in measurements, and find that none of them has a significant impact on the cosmological constraints.

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.

Upcoming multi-wavelength astronomical surveys will soon discover all massive galaxy clusters and provide unprecedented constraints on cosmology and cluster astrophysics. In this paper, we investigate the constraining power of the multi-band cluster surveys, through a joint analysis of three observables associated with clusters of galaxies, including thermal Sunyaev-Zel'dovich (tSZ) effect in cosmic microwave background (CMB), X-ray emission of ionized gas, and gravitational weak lensing effect of background galaxies by the cluster's gravitational potential. We develop a theoretical framework to predict and interpret two-point correlation statistics among the three observables using a semi-analytic model of intracluster medium (ICM) and halo-based approach. In this work, we show that the auto- and cross-angular power spectra in tSZ, X-ray and lensing statistics from upcoming missions (eROSITA, CMB-S4, and LSST) can help break the degeneracy between cosmology and ICM physics. These correlation statistics are less sensitive to selection biases, and are able to probe ICM physics in distant, faint and small clusters that are otherwise difficult to be detected individually. We show that the correlation statistics are able to provide cosmological constraints comparable to the conventional cluster abundance measurements, while constraining cluster astrophysics at the same time. Our results indicate that the correlation statistics can significantly enhance the scientific returns of upcoming multi-wavelength cluster surveys.

A mass of dark matter halo is commonly defined as the spherical over-density (SO) mass with respect to a reference density, whereas the time evolution of an SO mass can be affected by the redshift evolution of the reference density as well as the physical mass accretion around halos. In this study, we directly measure the amount of pseudo evolution of the SO masses of cluster-sized halos by the changes in the reference density from a time series of $N$-body simulations for the first time. We find that the $52\pm19\%$ difference in the virial SO masses between $z=0$ and $1$ can be accounted for by the pseudo evolution of clusters with a virial mass of $10^{14}\, h^{-1}M_{\odot}$ at $z=0$. The amount of pseudo evolution is found to be correlated with the age and density environment of a galaxy cluster. The stacked mass density profiles of cluster-sized halos with a greater amount of pseudo evolution in the SO mass shows the higher concentration ${\it and}$ greater linear bias parameter that is a counter-example of the known secondary halo bias due to concentration on the scale of clusters. We discuss how more concentrated clusters can show larger clustering amplitudes than their less concentrated counterparts and argue that the presence of rich filamentary structures plays a critical role in determining the linear halo bias of galaxy clusters.

Weak gravitational lensing is a powerful probe of the large-scale cosmic matter distribution. Wide-field galaxy surveys allow us to generate the so-called weak lensing maps, but actual observations suffer from noise due to imperfect measurement of galaxy shape distortions and to the limited number density of the source galaxies. In this paper, we explore a deep-learning approach to reduce the noise. We develop an image-to-image translation method with conditional adversarial networks (CANs), which learn efficient mapping from an input noisy weak lensing map to the underlying noise field. We train the CANs using $30000$ image pairs obtained from $1000$ ray-tracing simulations of weak gravitational lensing. We show that the trained CANs reproduce the true one-point probability distribution function (PDF) of the noiseless lensing map with a bias less than $1\sigma$ on average, where $\sigma$ is the statistical error. We perform a Fisher analysis to make forecast for cosmological parameter inference with the one-point lensing PDF. By our denoising method using CANs, the first derivative of the PDF with respect to the cosmic mean matter density and the amplitude of the primordial curvature perturbations becomes larger by $\sim50\%$. This allows us to improve the cosmological constraints by $\sim30-40\%$ with using observational data from ongoing and upcoming galaxy imaging surveys.

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 report the weak lensing measurement of filaments between Sloan Digital Sky Survey (SDSS) III/Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy pairs at $z\sim0.55$, using the Subaru Hyper Suprime-Cam (HSC) first-year galaxy shape catalogue. Despite of the small overlap of $140$ deg$^2$ between these surveys we detect the filament lensing signal at 3.9$\sigma$ significance, which is the highest signal-to-noise lensing measurement of filaments between galaxy-scale halos at this redshift range. We derive a theoretical prediction and covariance using mock catalogues based on full-sky ray-tracing simulations. We find that the intrinsic scatter of filament properties and the fluctuations in large scale structure along the line-of-sight are the primary component of the covariance and the intrinsic shape noise from source galaxies no longer limits our lensing measurement. This fact demonstrates the statistical power of the HSC survey due to its deep observations and high number density of source galaxies. Our result is consistent with the theoretical prediction and supports the "thick" filament model. As the HSC survey area increases, we will be able to study detailed filament properties such as the dark matter distributions and redshift evolution of filaments.

Two decades after its discovery, cosmic acceleration remains the most profound mystery in cosmology and arguably in all of physics. Either the Universe is dominated by a form of dark energy with exotic physical properties not predicted by standard model physics, or General Relativity is not an adequate description of gravity over cosmic distances. WFIRST emerged as a top priority of Astro2010 in part because of its ability to address the mystery of cosmic acceleration through both high precision measurements of the cosmic expansion history and the growth of cosmic structures with multiple and redundant probes. We illustrate in this white paper how mission design changes since Astro2010 have made WFIRST an even more powerful dark energy facility and have improved the ability of WFIRST to respond to changes in the experimental landscape. WFIRST is the space-based probe of DE the community needs in the mid-2020s.

We use the full-sky ray-tracing weak lensing simulations to generate 2268 mock catalogues for the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalogue. Our mock catalogues take into account various effects as in the real data: the survey footprints, inhomogeneous angular distribution of source galaxies, statistical uncertainties in photometric redshift (photo-$z$) estimate, variations in the lensing weight, and the statistical noise in galaxy shape measurements including both intrinsic shapes and the measurement errors. We then utilize our mock catalogues to evaluate statistical uncertainties expected in measurements of cosmic shear two-point correlations $\xi_{\pm}$ with tomographic redshift information for the HSC survey. We develop a quasi-analytical formula for the Gaussian sample variance properly taking into account the number of source pairs in the survey footprints. The standard Gaussian formula significantly overestimates or underestimates the mock results by $50\%$ level. We also show that different photo-$z$ catalogues or the six disconnected fields, rather than a consecutive geometry, cause variations in the covariance by $\sim 5\%$. The mock catalogues enable us to study the chi-square distribution for $\xi_{\pm}$. We find the wider distribution than that naively expected for the distribution with the degrees-of-freedom of data vector used. Finally, we propose a method to include non-zero multiplicative bias in mock shape catalogue and show the non-zero multiplicative bias can change the effective shape noise in cosmic shear analyses. Our results suggest an importance of estimating an accurate form of the likelihood function (and therefore the covariance) for robust cosmological parameter inference from the precise measurements.

The cross-correlation study of the unresolved $\gamma$-ray background (UGRB) with galaxy clusters has a potential to reveal the nature of the UGRB. In this paper, we perform a cross-correlation analysis between $\gamma$-ray data by the Fermi Large Area Telescope (Fermi-LAT) and a galaxy cluster catalogue from the Subaru Hyper Suprime-Cam (HSC) survey. The Subaru HSC cluster catalogue provides a wide and homogeneous large-scale structure distribution out to the high redshift at $z=1.1$, which has not been accessible in previous cross-correlation studies. We conduct the cross-correlation analysis not only for clusters in the all redshift range ($0.1 < z < 1.1$) of the survey, but also for subsamples of clusters divided into redshift bins, the low redshift bin ($0.1 < z < 0.6$) and the high redshift bin ($0.6 < z < 1.1$), to utilize the wide redshift coverage of the cluster catalogue. We find the evidence of the cross-correlation signals with the significance of 2.0-2.3$\sigma$ for all redshift and low-redshift cluster samples. On the other hand, for high-redshift clusters, we find the signal with weaker significance level (1.6-1.9$\sigma$). We also compare the observed cross-correlation functions with predictions of a theoretical model in which the UGRB originates from $\gamma$-ray emitters such as blazars, star-forming galaxies and radio galaxies. We find that the detected signal is consistent with the model prediction.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

Cross correlation with thermal Sunyaev-Zel'dovich (tSZ) effect in cosmic microwave background observation and weak gravitational lensing effect in galaxy imaging survey opens a new window on constraining matter contents in the Universe at redshift less than 1. In this paper, we study the impact of radio sources and cosmic infrared background (CIB) on observed tSZ-lensing correlation. Assuming the best-fit model of CIB by the Planck satellite, we estimate that the residual correlation of CIB with large-scale structures will be of an order of $2\%$ of expected tSZ-lensing correlation from intracluster medium in current lensing surveys. On the other hand, despite large model uncertainties, we find that correlation of lensing and radio sources can induce a negative correction for the observed tSZ-lensing correlation with a $\sim10-30\%$ level. This is originated from positive cross correlation with radio sources and lensing at $\sim100\, {\rm GHz}$ frequency, whereas tSZ-lensing correlation should show a negative value in temperature fluctuations at that frequency. We also show that the negative correction by radio-lensing correlation can solve the tension between recent measurements of tSZ-lensing correlation and an expected signal from "universal" gas pressure profile of nearby galaxy clusters, when the radio sources with flat-spectral index are assumed to populate massive cluster-sized dark matte halos. Our results indicate that minor source population in radio bands can play an important role in determining observed tSZ-lensing correlation at $<10$ arcmin.

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.

The extragalactic $\gamma$-ray background (EGB) arises from the accumulation of $\gamma$-ray emissions from resolved and unresolved extragalactic sources as well as diffuse processes. It is important to study the statistical properties of the EGB in the context of cosmological structure formation. Known astrophysical $\gamma$-ray sources such as blazars, star-forming galaxies, and radio galaxies are expected to trace the underlying cosmic matter density distribution. We explore the correlation of the EGB from Fermi-LAT data with the large-scale matter density distribution from the Subaru Hyper Suprime-Cam (HSC) SSP survey. We reconstruct an unbiased surface matter density distribution $\kappa$ at $z<1$ by applying weak-gravitational lensing analysis to the first-year HSC data. We then calculate the $\gamma - \kappa$ cross-correlation. Our measurements are consistent with a null detection, but a weak correlation is found at angular scales of 30-60 arcmin, especially when distant source galaxies at $z > 1$ are used for the lensing $\kappa$ reconstruction. The large-scale correlation suggests strong clustering of high-redshift $\gamma$-ray sources such as blazars. However, the inferred bias factor of $4-5$ is larger by about a factor of two than results from other clustering analyses. The final HSC data covering 1,400 squared degrees will play an essential role to determine accurately the blazar bias at $z > 0.5$.

Stacked lensing is a powerful means of measuring the average mass distribution around large-scale structure tracers. There are two stacked lensing estimators used in the literature, denoted as $\Delta\Sigma$ and $\gamma_+$, which are related as $\Delta\Sigma=\Sigma_{\rm cr}\gamma_+$, where $\Sigma_{\rm cr}(z_l,z_s)$ is the critical surface mass density for each lens-source pair ($z_l$ and $z_s$ are lens and source redshifts, respectively). In this paper we derive a formula for the covariance matrix of $\Delta\Sigma$-estimator focusing on `weight' function to improve the signal-to-noise ($S/N$). We assume that the lensing fields and the distribution of lensing objects obey the Gaussian statistics. With this formula, we show that, if background galaxy shapes are weighted by an amount of $\Sigma_{\rm cr}^{-2}(z_l,z_s)$, the $\Delta\Sigma$-estimator maximizes the $S/N$ in the shot noise limited regime. We also show that the $\Delta\Sigma$-estimator with the weight $\Sigma_{\rm cr}^{-2}$ gives a greater $(S/N)^2$ than that of the $\gamma_+$-estimator by about 5--25\% for lensing objects at redshifts comparable with or higher than the median of source galaxy redshifts for hypothetical Subaru HSC and DES surveys. However, for low-redshift lenses such as $z_l<0.3$, the $\gamma_+$-estimator has higher $(S/N)^2$ than $\Delta\Sigma$. We also discuss that the $(S/N)^2$ for $\Delta\Sigma$ at large separations in the sample variance limited regime can be boosted, by up to a factor of 1.5, if one adopts a weight of $\Sigma_{\rm cr}^{-\alpha}$ with $\alpha>2$. Our formula allows one to explore how the combination of the different estimators can approach an optimal estimator in all regimes of redshifts and separation scales.

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.

Gravitational lensing is a powerful probe of the mass distribution of galaxy clusters and cosmology. However, accurate measurements of the cluster mass profiles are limited by uncertainties in cluster astrophysics. In this work, we present a physically motivated model of baryonic effects on the cluster mass profiles, which self-consistently takes into account the impact of baryons on the concentration as well as mass accretion histories of galaxy clusters. We calibrate this model using the Omega500 hydrodynamical cosmological simulations of galaxy clusters with varying baryonic physics. Our model will enable us to simultaneously constrain cluster mass, concentration, and cosmological parameters using stacked weak lensing measurements from upcoming optical cluster surveys.

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.

We present the result of searching for clusters of galaxies based on weak gravitational lensing analysis of the $\sim 160$~deg$^2$ area surveyed by Hyper Suprime-Cam (HSC) as a Subaru Strategic Program. HSC is a new prime focus optical imager with a 1.5 diameter field of view on the 8.2-meter Subaru telescope. The superb median seeing on the HSC $i$-band images of $0.56$ arcsec allows the reconstruction of high angular resolution mass maps via weak lensing, which is crucial for the weak lensing cluster search. We identify 65 mass map peaks with signal-to-noise (SN) ratio larger than 4.7, and carefully examine their properties by cross-matching the clusters with optical and X-ray cluster catalogs. We find that all the 39 peaks with SN$>5.1$ have counterparts in the optical cluster catalogs, and only 2 out of the 65 peaks are probably false positives. The upper limits of X-ray luminosities from ROSAT All Sky Survey (RASS) imply the existence of an X-ray under-luminous cluster population. We show that the X-rays from the shear selected clusters can be statistically detected by stacking the RASS images. The inferred average X-ray luminosity is about half that of the X-ray selected clusters of the same mass. The radial profile of the dark matter distribution derived from the stacking analysis is well modeled by the Navarro-Frenk-White profile with a small concentration parameter value of $c_{500}\sim 2.5$, which suggests that the selection bias on the orientation or the internal structure for our shear selected cluster sample is not strong.

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 and characterize the catalog of galaxy shape measurements that will be used for cosmological weak lensing measurements in the Wide layer of the first year of the Hyper Suprime-Cam (HSC) survey. The catalog covers an area of 136.9 deg$^2$ split into six fields, with a mean $i$-band seeing of $0.58$ arcsec and $5\sigma$ point-source depth of $i\sim 26$. Given conservative galaxy selection criteria for first year science, the depth and excellent image quality results in unweighted and weighted source number densities of 24.6 and 21.8 arcmin$^{-2}$, respectively. We define the requirements for cosmological weak lensing science with this catalog, then focus on characterizing potential systematics in the catalog using a series of internal null tests for problems with point-spread function (PSF) modeling, shear estimation, and other aspects of the image processing. We find that the PSF models narrowly meet requirements for weak lensing science with this catalog, with fractional PSF model size residuals of approximately $0.003$ (requirement: 0.004) and the PSF model shape correlation function $\rho_1<3\times 10^{-7}$ (requirement: $4\times 10^{-7}$) at 0.5$^\circ$ scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star-galaxy shape correlations reveal additive systematics on $>1^\circ$ scales that are sufficiently large as to require mitigation in cosmic shear measurements. Finally, we discuss the dominant systematics and the planned algorithmic changes to reduce them in future data reductions.

Weak lensing three-point statistics are powerful probes of the structure of dark matter halos. We propose to use the correlation of the positions of galaxies with the shapes of background galaxy pairs, known as the halo-shear-shear correlation (HSSC), to measure the mean halo ellipticity and the abundance of subhalos in a statistical manner. We run high-resolution cosmological $N$-body simulations and use the outputs to measure the HSSC for galaxy halos and cluster halos. Non-spherical halos cause a characteristic azimuthal variation of the HSSC, and massive subhalos in the outer region near the virial radius contribute to $\sim10\%$ of the HSSC amplitude. Using the HSSC and its covariance estimated from our $N$-body simulations, we make forecast for constraining the internal structure of dark matter halos with future galaxy surveys. With 1000 galaxy groups with mass greater than $10^{13.5}\, h^{-1}M_{\odot}$, the average halo ellipticity can be measured with an accuracy of ten percent. A spherical, smooth mass distribution can be ruled out at a $\sim5\sigma$ significance level. The existence of subhalos whose masses are in 1-10 percent of the main halo mass can be detected with $\sim10^4$ galaxies/clusters. We conclude that the HSSC provides valuable information on the structure of dark halos and hence on the nature of dark matter.

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 an analysis of galaxy-galaxy weak gravitational lensing (GGL) in chameleon $f(R)$ gravity - a leading candidate of non-standard gravity models. For the analysis we have created mock galaxy catalogues based on dark matter haloes from two sets of numerical simulations, using a halo occupation distribution (HOD) prescription which allows a redshift dependence of galaxy number density. To make a fairer comparison between the $f(R)$ and $\Lambda$CDM models, their HOD parameters are tuned so that the galaxy two-point correlation functions in real space (and therefore the projected two-point correlation functions) match. While the $f(R)$ model predicts an enhancement of the convergence power spectrum by up to $\sim30\%$ compared to the standard $\Lambda$CDM model with the same parameters, the maximum enhancement of GGL is only half as large and less than 5\% on separations above $\sim1$-$2h^{-1}$Mpc, because the latter is a cross correlation of shear (or matter, which is more strongly affected by modified gravity) and galaxy (which is weakly affected given the good match between galaxy auto correlations in the two models) fields. We also study the possibility of reconstructing the matter power spectrum by combination of GGL and galaxy clustering in $f(R)$ gravity. We find that the galaxy-matter cross correlation coefficient remains at unity down to $\sim2$-$3h^{-1}$Mpc at relevant redshifts even in $f(R)$ gravity, indicating joint analysis of GGL and galaxy clustering can be a powerful probe of matter density fluctuations in chameleon gravity. The scale dependence of the model differences in their predictions of GGL can potentially allow to break the degeneracy between $f(R)$ gravity and other cosmological parameters such as $\Omega_m$ and $\sigma_8$.

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 propose to use degree-scale angular clustering of fast radio bursts (FRBs) to identify their origin and the host galaxy population. We study the information content in autocorrelation of the angular positions and dispersion measures (DM) and in cross-correlation with galaxies. We show that the cross-correlation with Sloan Digital Sky Survey (SDSS) galaxies will place stringent constraints on the mean physical quantities associated with FRBs. If $\sim$10,000 FRBs are detected with $\lesssim \rm deg$ resolution in the SDSS field, the clustering analysis with the intrinsic DM scatter of $100\, {\rm pc}/{\rm cm}^3$ can constrain the global abundance of free electrons at $z\lt1$ and the large-scale bias of FRB host galaxies (the statistical relation between the distribution of host galaxies and cosmic matter density field) with fractional errors (with a $68\%$ confidence level) of $\sim10\%$ and $\sim20\%$, respectively. The mean near-source dispersion measure and the delay time distribution of FRB rates relative to the global star forming rate can be also determined by combining the clustering and the probability distribution function of DM. Our approach will be complementary to high-resolution ($\ll {\rm deg}$) event localization using e.g., VLA and VLBI for identifying the origin of FRBs and the source environment. We strongly encourage future observational programs such as CHIME, UTMOST, and HIRAX to survey FRBs in the SDSS field.

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.

The number density of local maxima of weak lensing field, referred to as weak-lensing peak counts, can be used as a cosmological probe. However, its relevant cosmological information is still unclear. We study the relationship between the peak counts and other statistics in weak lensing field by using 1000 ray-tracing simulations. We construct a local transformation of lensing field $\cal K$ to a new Gaussian field $y$, named local-Gaussianized transformation. We calibrate the transformation with numerical simulations so that the one-point distribution and the power spectrum of $\cal K$ can be reproduced from a single Gaussian field $y$ and monotonic relation between $y$ and $\cal K$. Therefore, the correct information of two-point clustering and any order of moments in weak lensing field should be preserved under local-Gaussianized transformation. We then examine if local-Gaussianized transformation can predict weak-lensing peak counts in simulations. The local-Gaussianized transformation is insufficient to explain weak-lensing peak counts in the absence of shape noise. The prediction by local-Gaussianized transformation underestimates the simulated peak counts with a level of $\sim20-30\%$ over a wide range of peak heights. Local-Gaussianized transformation can predict the weak-lensing peak counts with a $\sim10\%$ accuracy in the presence of shape noise. Our analyses suggest that the cosmological information beyond power spectrum and its moments would be necessary to predict the weak-lensing peak counts with a percent-level accuracy, which is an expected statistical uncertainty in upcoming wide-field galaxy surveys.

We derive constraints on dark matter (DM) annihilation cross section and decay lifetime from cross-correlation analyses of the data from Fermi-LAT and weak lensing surveys that cover a wide area of $\sim660$ squared degrees in total. We improve upon our previous analyses by using an updated extragalactic $\gamma$-ray background data reprocessed with the Fermi Pass 8 pipeline, and by using well-calibrated shape measurements of about twelve million galaxies in the Canada-France-Hawaii Lensing Survey (CFHTLenS) and Red-Cluster-Sequence Lensing Survey (RCSLenS). We generate a large set of full-sky mock catalogs from cosmological $N$-body simulations and use them to estimate statistical errors accurately. The measured cross correlation is consistent with null detection, which is then used to place strong cosmological constraints on annihilating and decaying DM. For leptophilic DM, the constraints are improved by a factor of $\sim100$ in the mass range of O(1) TeV when including contributions from secondary $\gamma$ rays due to the inverse-Compton upscattering of background photons. Annihilation cross-sections of $\langle \sigma v \rangle \sim 10^{-23}\, {\rm cm}^3/{\rm s}$ are excluded for TeV-scale DM depending on channel. Lifetimes of $\sim 10^{25}$ sec are also excluded for the decaying TeV-scale DM. Finally, we apply this analysis to wino DM and exclude the wino mass around 200 GeV. These constraints will be further tightened, and all the interesting wino DM parameter region can be tested, by using data from future wide-field cosmology surveys.

The thermal Sunyaev-Zel'dovich (tSZ) effect signal is widely recognized as a robust mass proxy of galaxy clusters with small intrinsic scatter. However, recent observational calibration of the tSZ scaling relation using weak lensing (WL) mass exhibits considerably larger scatter than the intrinsic scatter predicted from numerical simulations. This raises a question as to whether we can realize the full statistical power of ongoing and upcoming tSZ-WL observations of galaxy clusters. In this work, we investigate the origin of observed scatter in the tSZ-WL scaling relation, using mock maps of galaxy clusters extracted from cosmological hydrodynamic simulations. We show that the inferred intrinsic scatter from mock tSZ-WL analyses is considerably larger than the intrinsic scatter measured in simulations, and comparable to the scatter in the observed tSZ-WL relation. We show that this enhanced scatter originates from the combination of the projection of correlated structures along the line of sight and the uncertainty in the cluster radius associated with WL mass estimates, causing the amplitude of the scatter to depend on the covariance between tSZ and WL signals. We present a statistical model to recover the unbiased cluster scaling relation and cosmological parameter by taking into account the covariance in the tSZ-WL mass relation from multi-wavelength cluster surveys.

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.

We study the effect of $f(R)$ gravity on the statistical properties of various large-scale structures which can be probed in weak gravitational lensing measurements. A set of ray-tracing simulations of gravitational lensing in $f(R)$ gravity enables us to explore cosmological information on (i) stacking analyses of weak lensing observables and (ii) peak statistics in reconstructed lensing mass maps. For the $f(R)$ model proposed by Hu \& Sawicki, the measured lensing signals of dark matter haloes in the stacking analysis would show a $\simlt10\%$ difference between the standard $\Lambda$CDM and the $f(R)$ model when the additional degree of freedom in $f(R)$ model would be $|f_{\rm R0}|\sim10^{-5}$. Among various large-scale structures to be studied in stacking analysis, troughs, i.e, underdensity regions in projected plane of foreground massive haloes, could be promising to constrain the model with $|f_{\rm R0}|\sim10^{-5}$, while stacking analysis around voids is found to be difficult to improve the constraint of $|f_{\rm R0}|$ even in future lensing surveys with a sky coverage of $\sim1000$ square degrees. On the peak statistics, we confirm the correspondence between local maxima and dark matter haloes along the line of sight, regardless of the modification of gravity in our simulation. Thus, the number count of high significance local maxima would be useful to probe the mass function of dark matter haloes even in the $f(R)$ model with $|f_{\rm R0}|\simlt10^{-5}$. We also find that including local minima in lensing mass maps would be helpful to improve the constant on $f(R)$ gravity down to $|f_{\rm R0}|=10^{-5}$ in ongoing weak lensing surveys.

Measurements of the cross-correlation between the extragalactic gamma-ray background (EGB) and large-scale structure provide a novel probe of dark matter on extragalactic scales. We focus on luminous red galaxies (LRGs) as optimal targets to search for the signal of dark matter annihilation. We measure the cross-correlation function of the EGB taken from the Fermi Large Area Telescope with the LRGs from the Sloan Digital Sky Survey. Statistical errors are calculated using a large set of realistic mock LRG catalogs. The amplitude of the measured cross-correlation is consistent with null detection. Based on an accurate theoretical model of the distribution of dark matter associated with LRGs, we exclude dark matter annihilation cross-sections over $\langle \sigma v\rangle =3\times10^{-25}-10^{-26}\, {\rm cm}^3 \,{\rm s}^{-1}$ for a 10 GeV dark matter. We further investigate systematic effects due to uncertainties in the Galactic gamma-ray foreground emission, which we find to be an order of magnitude smaller than the current statistical uncertainty. We also estimate the contamination from astrophysical sources in the LRGs by using known scaling relations between gamma-ray luminosity and star-formation rate, finding them to be negligibly small. Based on these results, we suggest that LRGs remain ideal targets for probing dark matter annihilation with future EGB measurement and galaxy surveys. Increasing the number of LRGs in upcoming galaxy surveys such as LSST would lead to big improvements of factors of several in sensitivity.

Ongoing and future wide-field galaxy surveys can be used to locate a number of clusters of galaxies with cosmic shear measurement alone. We study constraints on cosmological models using statistics of weak lensing selected galaxy clusters. We extend our previous theoretical framework to model the statistical properties of clusters in variants of cosmological models as well as in the standard $\Lambda$CDM model. Weak lensing selection of clusters does not rely on the conventional assumption such as the relation between luminosity and mass and/or hydrostatic equilibrium, but a number of observational effects compromise robust identification. We use a large set of realistic mock weak-lensing catalogs as well as analytic models to perform a Fisher analysis and make forecast for constraining two competing cosmological models, $w$CDM model and $f(R)$ model proposed by Hu & Sawicki, with our lensing statistics. We show that weak lensing selected clusters are excellent probe of cosmology when combined with cosmic shear power spectrum even in presence of galaxy shape noise and masked regions. With the information of weak lensing selected clusters, the precision of cosmological parameter estimate can be improved by a factor of $\sim1.6$ and $\sim8$ for $w$CDM model and $f(R)$ model, respectively. Hyper Suprime-Cam survey with sky coverage of $1250$ squared degrees can constrain the equation of state of dark energy $w_{0}$ with a level of $\Delta w_0 \sim0.1$. It can also constrain the additional scalar degree of freedom in $f(R)$ model with a level of $|f_{R0}| \sim5\times10^{-6}$, when constraints from cosmic microwave background measurements are incorporated. Future weak lensing surveys with sky coverage of $20,000$ squared degrees will place tighter constraints on $w_{0}$ and $|f_{R0}|$ even without cosmic microwave background measurements.

We explore a variety of statistics of clusters selected with cosmic shear measurement by utilizing both analytic models and large numerical simulations. We first develop a halo model to predict the abundance and the clustering of weak lensing selected clusters. Observational effects such as galaxy shape noise are included in our model. We then generate realistic mock weak lensing catalogs to test the accuracy of our analytic model. To this end, we perform full-sky ray-tracing simulations that allow us to have multiple realizations of a large continuous area. We model the masked regions on the sky using the actual positions of bright stars, and generate 200 mock weak lensing catalogs with sky coverage of ~1000 squared degrees. We show that our theoretical model agrees well with the ensemble average of statistics and their covariances calculated directly from the mock catalogues. With a typical selection threshold, ignoring shape noise correction causes overestimation of the clustering of weak lensing selected clusters with a level of about 10%, and shape noise correction boosts the cluster abundance by a factor of a few. We calculate the cross-covariances using the halo model with accounting for the effective reduction of the survey area due to masks. The covariance of the cosmic shear auto power spectrum is affected by the mode-coupling effect that originates from sky masking. Our model and the results can be readily used for cosmological analysis with ongoing and future weak lensing surveys.

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.

We propose a novel method to select satellite galaxies in outer regions of galaxy groups or clusters using weak gravitational lensing. The method is based on the theoretical expectation that the tangential shear pattern around satellite galaxies would appear with negative values at an offset distance from the center of the main halo. We can thus locate the satellite galaxies statistically with an offset distance of several lensing smoothing scales by using the standard reconstruction of surface mass density maps from weak lensing observation. We test the idea using high-resolution cosmological simulations. We show that subhalos separated from the center of the host halo are successfully located even without assuming the position of the center. For a number of such subhalos, the characteristic mass and offset length can be also estimated on a statistical basis. We perform a Fisher analysis to show how well upcoming weak lensing surveys can constrain the mass density profile of satellite galaxies. In the case of the Large Synoptic Survey Telescope with a sky coverage of 20,000 square degrees, the mass of the member galaxies in the outer region of galaxy clusters can be constrained with an accuracy of ~0.1 dex for galaxy clusters with mass 10^14 Msun/h at z=0.15. Finally we explore the detectability of tidal stripping features for subhalos having a wide range of masses of 10^11-10^13 Msun/h.

We present the first measurement of the cross-correlation of weak gravitational lensing and the extragalactic gamma-ray background emission using data from the Canada-France-Hawaii Lensing Survey and the Fermi Large Area Telescope. The cross-correlation is a powerful probe of signatures of dark matter annihilation, because both cosmic shear and gamma-ray emission originate directly from the same DM distribution in the universe, and it can be used to derive constraints on dark matter annihilation cross-section. We show that the measured lensing-gamma correlation is consistent with a null signal. Comparing the result to theoretical predictions, we exclude dark matter annihilation cross sections of <sigma v> =10^{-24}-10^{-25} cm^3 s^-1 for a 100 GeV dark matter. If dark matter halos exist down to the mass scale of 10^-6 M_sun, we are able to place constraints on the thermal cross sections <sigma v> ~ 3 x 10^{-26} cm^3 s^-1 for a 10 GeV dark matter annihilation into tau^{+} tau^{-}. Future gravitational lensing surveys will increase sensitivity to probe annihilation cross sections of <sigma v> ~ 5 x 10^{-26} cm^3 s^-1 even for a 100 GeV dark matter. Detailed modeling of the contributions from astrophysical sources to the cross correlation signal could further improve the constraints by ~ 40-70 %.

The measurement of cosmic shear using weak gravitational lensing is a challenging task that involves a number of complicated procedures. We study in detail the systematic errors in the measurement of weak lensing Minkowski Functionals (MFs). Specifically, we focus on systematics associated with galaxy shape measurements, photometric redshift errors, and shear calibration correction. We first generate mock weak-lensing catalogs that directly incorporate the actual observational characteristics of the Canada-France-Hawaii Lensing Survey (CFHTLenS). We then perform a Fisher analysis using the large set of mock catalogs for various cosmological models. We find that the statistical error associated with the observational effects degrades the cosmological parameter constraints by a factor of a few. The Subaru Hyper Suprime-Cam (HSC) survey with a sky coverage of ~1400 deg2 will constrain the dark energy equation of the state parameter with an error of Delta w_0 ~ 0.25 by the lensing MFs alone, but biases induced by the systematics can be comparable to the 1-sigma error. We conclude that the lensing MFs are powerful statistics beyond the two-point statistics, only if well-calibrated measurement of both the redshifts and the shapes of source galaxies is performed. Finally, we analyze the CFHTLenS data to explore the ability of the MFs to break degeneracies between a few cosmological parameters. Using a combined analysis of the MFs and the shear correlation function, we derive the matter density Omega_m0= 0.256+0.054-0.046.

We study weak lensing properties of filaments that connect clusters of galaxies through large cosmological $N$-body simulations. We select 4639 halo pairs with masses higher than $10^{14}h^{-1}\mathrm{M}_\odot$ from the simulations and investigate dark matter distributions between two haloes with ray-tracing simulations. In order to classify filament candidates, we estimate convergence profiles and perform profile fitting. We find that matter distributions between haloes can be classified in a plane of fitting parameters, which allow us to select straight filaments from the ray-tracing simulations. We also investigate statistical properties of these filaments, finding them to be consistent with previous studies. We find that $35\%$ of halo pairs possess straight filaments, $4\%$ of which can directly be detected at $S/N\geq2$ with weak lensing. Furthermore, we study statistical properties of haloes at the edges of filaments. We find that haloes are preferentially elongated along filamentary structures and are less massive with increasing filament masses. However, the dependence of these halo properties on masses of straight filaments is very weak.

We explore the discovery potential of light gravitino mass m_{3/2} by combining future cosmology surveys and collider experiments. The former probe the imprint of light gravitinos in the cosmic matter density field, whereas the latter search signatures of a supersymmetry breaking mechanism. Free-streaming of light gravitinos suppresses the density fluctuations at galactic and sub-galactic length scales, where weak gravitational lensing can be used as a powerful probe. We perform numerical simulations of structure formation to quantify the effect. We then run realistic ray-tracing simulations of gravitational lensing to measure the cosmic shear in models with light gravitino. We forecast the possible reach of future wide-field surveys by Fisher analysis; the light gravitino mass can be determined with an accuracy of m_{3/2}=4\pm 1 eV by a combination of the Hyper Suprime Cam survey and cosmic microwave background anisotropy data obtained by Planck satellite. The corresponding accuracy to be obtained by the future Large Synoptic Survey Telescope is \delta m_{3/2}=0.6 eV. Data from experiments at Large Hadron Collider at 14 TeV will provide constraint at m_{3/2} \simeq 5 eV in the minimal framework of gauge-mediated supersymmetry breaking (GMSB) model. We conclude that a large class of the GMSB model can be tested by combining the cosmological observations and the collider experiments.

Sky masking is unavoidable in wide-field weak lensing observations. We study how masks affect the measurement of statistics of matter distribution probed by weak gravitational lensing. We first use 1000 cosmological ray-tracing simulations to examine in detail the impact of masked regions on the weak lensing Minkowski Functionals (MFs). We consider actual sky masks used for a Subaru Suprime-Cam imaging survey. The masks increase the variance of the convergence field and the expected values of the MFs are biased. The bias then affects the non-Gaussian signals induced by the gravitational growth of structure. We then explore how masks affect cosmological parameter estimation. We calculate the cumulative signal-to-noise ratio S/N for masked maps to study the information content of lensing MFs. We show that the degradation of S/N for masked maps is mainly determined by the effective survey area. We also perform simple chi^2 analysis to show the impact of lensing MF bias due to masked regions. Finally, we compare ray-tracing simulations with data from a Subaru 2 deg^2 survey in order to address if the observed lensing MFs are consistent with those of the standard cosmology. The resulting chi^2/n_dof = 29.6/30 for three combined MFs, obtained with the mask effects taken into account, suggests that the observational data are indeed consistent with the standard LambdaCDM model. We conclude that the lensing MFs are powerful probe of cosmology only if mask effects are correctly taken into account.

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 examine scatter and bias in weak lensing selected clusters, employing both an analytic model of dark matter haloes and numerical mock data of weak lensing cluster surveys. We pay special attention to effects of the diversity of dark matter distributions within clusters. We find that peak heights of the lensing convergence map correlates rather poorly with the virial mass of haloes. The correlation is tighter for the spherical overdensity mass with a higher mean interior density. We examine the dependence of the halo shape on the peak heights, and find that the rms scatter caused by the halo diversity scales linearly with the peak heights with the proportionality factor of 0.1-0.2. The noise originated from the halo shape is found to be comparable to the source galaxy shape noise and the cosmic shear noise. We find the significant halo orientation bias, i.e., weak lensing selected clusters on average have their major axes aligned with the line-of-sight direction. We compute the orientation bias using an analytic triaxial halo model and obtain results quite consistent with the ray-tracing results. We develop a prescription to analytically compute the number count of weak lensing peaks taking into account all the main sources of scatters in peak heights. We find that the improved analytic predictions agree well with the simulation results for high S/N peaks. We also compare the expected number count with our weak lensing analysis results for 4 sq deg of Subaru/Suprime-Cam observations and find a good agreement.