We present a statistical weak-lensing magnification analysis on an optically selected sample of 3029 \texttt{CAMIRA} galaxy clusters with richness $N>15$ at redshift $0.2\leq z <1.1$ in the Subaru Hyper Suprime-Cam (HSC) survey. We use two distinct populations of color-selected, flux-limited background galaxies, namely the low-$z$ and high-$z$ samples at mean redshifts of $\approx1.1$ and $\approx1.4$, respectively, from which to measure the weak-lensing magnification signal by accounting for cluster contamination as well as masking effects. Our magnification bias measurements are found to be uncontaminated according to validation tests against the "null-test" samples for which the net magnification bias is expected to vanish. The magnification bias for the full \texttt{CAMIRA} sample is detected at a significance level of $9.51\sigma$, which is dominated by the high-$z$ background. We forward-model the observed magnification data to constrain the normalization of the richness-to-mass ($N$--$M$) relation for the \texttt{CAMIRA} sample with informative priors on other parameters. The resulting scaling relation is $N\propto {M_{500}}^{0.92\pm0.13} (1 + z)^{-0.48\pm0.69}$, with a characteristic richness of $N=\left(17.72\pm2.60\right)$ and intrinsic log-normal scatter of $0.15\pm0.07$ at $M_{500} = 10^{14}h^{-1}M_{\odot}$. With the derived $N$--$M$ relation, we provide magnification-calibrated mass estimates of individual \texttt{CAMIRA} clusters, with the typical uncertainty of $\approx39\%$ and $\approx32\%$ at richness$\approx20$ and $\approx40$, respectively. We further compare our magnification-inferred $N$--$M$ relation with those from the shear-based results in the literature, finding good agreement.

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.

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

We present results from a set of simulations designed to constrain the weak lensing shear calibration for the Hyper Suprime-Cam (HSC) survey. These simulations include HSC observing conditions and galaxy images from the Hubble Space Telescope (HST), with fully realistic galaxy morphologies and the impact of nearby galaxies included. We find that the inclusion of nearby galaxies in the images is critical to reproducing the observed distributions of galaxy sizes and magnitudes, due to the non-negligible fraction of unrecognized blends in ground-based data, even with the excellent typical seeing of the HSC survey (0.58" in the $i$-band). Using these simulations, we detect and remove the impact of selection biases due to the correlation of weights and the quantities used to define the sample (S/N and apparent size) with the lensing shear. We quantify and remove galaxy property-dependent multiplicative and additive shear biases that are intrinsic to our shear estimation method, including a $\sim 10$ per cent-level multiplicative bias due to the impact of nearby galaxies and unrecognized blends. Finally, we check the sensitivity of our shear calibration estimates to other cuts made on the simulated samples, and find that the changes in shear calibration are well within the requirements for HSC weak lensing analysis. Overall, the simulations suggest that the weak lensing multiplicative biases in the first-year HSC shear catalog are controlled at the 1 per cent level.

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.

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 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 optimized source galaxy selection schemes for measuring cluster weak lensing (WL) mass profiles unaffected by cluster member dilution from the Subaru Hyper Suprime-Cam Strategic Survey Program (HSC-SSP). The ongoing HSC-SSP survey will uncover thousands of galaxy clusters to $z\lesssim1.5$. In deriving cluster masses via WL, a critical source of systematics is contamination and dilution of the lensing signal by cluster {members, and by foreground galaxies whose photometric redshifts are biased}. Using the first-year CAMIRA catalog of $\sim$900 clusters with richness larger than 20 found in $\sim$140 deg$^2$ of HSC-SSP data, we devise and compare several source selection methods, including selection in color-color space (CC-cut), and selection of robust photometric redshifts by applying constraints on their cumulative probability distribution function (PDF; P-cut). We examine the dependence of the contamination on the chosen limits adopted for each method. Using the proper limits, these methods give mass profiles with minimal dilution in agreement with one another. We find that not adopting either the CC-cut or P-cut methods results in an underestimation of the total cluster mass ($13\pm4\%$) and the concentration of the profile ($24\pm11\%$). The level of cluster contamination can reach as high as $\sim10\%$ at $R\approx 0.24$ Mpc/$h$ for low-z clusters without cuts, while employing either the P-cut or CC-cut results in cluster contamination consistent with zero to within the 0.5% uncertainties. Our robust methods yield a $\sim60\sigma$ detection of the stacked CAMIRA surface mass density profile, with a mean mass of $M_\mathrm{200c} = (1.67\pm0.05({\rm {stat}}))\times 10^{14}\,M_\odot/h$.

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.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

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.

In this paper, we describe the optical imaging data processing pipeline developed for the Subaru Telescope's Hyper Suprime-Cam (HSC) instrument. The HSC Pipeline builds on the prototype pipeline being developed by the Large Synoptic Survey Telescope's Data Management system, adding customizations for HSC, large-scale processing capabilities, and novel algorithms that have since been reincorporated into the LSST codebase. While designed primarily to reduce HSC Subaru Strategic Program (SSP) data, it is also the recommended pipeline for reducing general-observer HSC data. The HSC pipeline includes high level processing steps that generate coadded images and science-ready catalogs as well as low-level detrending and image characterizations.

The Subaru Strategic Program (SSP) is an ambitious multi-band survey using the Hyper Suprime-Cam (HSC) on the Subaru telescope. The Wide layer of the SSP is both wide and deep, reaching a detection limit of i~26.0 mag. At these depths, it is challenging to achieve accurate, unbiased, and consistent photometry across all five bands. The HSC data are reduced using a pipeline that builds on the prototype pipeline for the Large Synoptic Survey Telescope. We have developed a Python-based, flexible framework to inject synthetic galaxies into real HSC images called SynPipe. Here we explain the design and implementation of SynPipe and generate a sample of synthetic galaxies to examine the photometric performance of the HSC pipeline. For stars, we achieve 1% photometric precision at i~19.0 mag and 6% precision at i~25.0 in the i-band. For synthetic galaxies with single-Sersic profiles, forced CModel photometry achieves 13% photometric precision at i~20.0 mag and 18% precision at i~25.0 in the i-band. We show that both forced PSF and CModel photometry yield unbiased color estimates that are robust to seeing conditions. We identify several caveats that apply to the version of HSC pipeline used for the first public HSC data release (DR1) that need to be taking into consideration. First, the degree to which an object is blended with other objects impacts the overall photometric performance. This is especially true for point sources. Highly blended objects tend to have larger photometric uncertainties, systematically underestimated fluxes and slightly biased colors. Second, >20% of stars at 22.5< i < 25.0 mag can be misclassified as extended objects. Third, the current CModel algorithm tends to strongly underestimate the half-light radius and ellipticity of galaxy with i>21.5 mag.

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.