• ### A large sample of shear selected clusters from the Hyper Suprime-Cam Subaru Strategic Program S16A wide field mass maps(1802.10290)

Feb. 28, 2018 astro-ph.CO
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.
• ### X-ray properties of high-richness CAMIRA clusters in the Hyper Suprime-Cam Subaru Strategic Program field(1802.08692)

Feb. 23, 2018 astro-ph.CO
We present the first results of a pilot X-ray study of 17 rich galaxy clusters at $0.14<z<0.75$ in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) field. Diffuse X-ray emissions from these clusters were serendipitously detected in the XMM-Newton fields of view. We systematically analyze X-ray images and emission spectra of the hot intracluster gas by using the XMM-Newton archive data. The frequency distribution of the offset between the X-ray centroid or peak and the position of the brightest cluster galaxy was derived for the optically-selected cluster sample. The fraction of relaxed clusters estimated from the X-ray peak offsets is $30\pm13$%, which is smaller than that of the X-ray cluster samples such as HIFLUGCS. Since the optical cluster search is immune to the physical state of X-ray-emitting gas, it is likely to cover a larger range of the cluster morphology. We also derived the luminosity-temperature relation and found that the slope is marginally shallower than those of X-ray-selected samples and consistent with the self-similar model prediction of 2. Accordingly, our results show that the X-ray properties of the optically-selected clusters are marginally different from those observed in the X-ray samples. The implication of the results and future prospects are briefly discussed.
• ### The first-year shear catalog of the Subaru Hyper Suprime-Cam SSP Survey(1705.06745)

Dec. 26, 2017 astro-ph.CO, astro-ph.IM
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.
• ### Multiwavelength study of X-ray Luminous Clusters in the Hyper Suprime-Cam Subaru Strategic Program S16A field(1705.06852)

Oct. 18, 2017 astro-ph.CO
We present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.
• ### Full-sky Gravitational Lensing Simulation for Large-area Galaxy Surveys and Cosmic Microwave Background Experiments(1706.01472)

Sept. 26, 2017 gr-qc, astro-ph.CO
We present 108 full-sky gravitational lensing simulation datasets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological $N$-body simulations. The datasets 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 datasets 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.
• ### Robust covariance estimation of galaxy-galaxy weak lensing: validation and limitation of jackknife covariance(1607.08679)

June 28, 2017 astro-ph.CO
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.
• ### The Galaxy--Halo Connection in High-Redshift Universe: Details and Evolution of Stellar-to-Halo Mass Ratios of Lyman Break Galaxies on CFHTLS Deep Fields(1612.06869)

April 5, 2017 astro-ph.CO, astro-ph.GA
We present the results of clustering analyses of Lyman break galaxies (LBGs) at $z\sim3$, $4$, and $5$ using the final data release of the Canada--France--Hawaii Telescope Legacy Survey (CFHTLS). Deep- and wide-field images of the CFHTLS Deep Survey enable us to obtain sufficiently accurate two-point angular correlation functions to apply a halo occupation distribution analysis. Mean halo masses, calculated as $\langle M_{h} \rangle=10^{11.7}-10^{12.8}h^{-1}M_{\odot}$, increase with stellar-mass limit of LBGs. The threshold halo mass to have a central galaxy follows the same increasing trend with the low-$z$ results, whereas the threshold halo mass to have a satellite galaxy shows higher values at $z=3-5$ than $z=0.5-1.5$ over the entire stellar mass range. Satellite fractions of dropout galaxies, even at less massive haloes, are found to drop sharply from $z=2$ down to less than $0.04$ at $z=3-5$. These results suggest that satellite galaxies form inefficiently within dark haloes at $z=3-5$ even for less massive satellites with $M_{\star}<10^{10}M_{\odot}$. We compute stellar-to-halo mass ratios (SHMRs) assuming a main sequence of galaxies, which is found to provide consistent SHMRs with those derived from a spectral energy distribution fitting method. The observed SHMRs are in good agreement with the model predictions based on the abundance-matching method within $1\sigma$ confidence intervals. We derive observationally, for the first time, $M_{{\rm h}}^{{\rm pivot}}$, which is the halo mass at a peak in the star-formation efficiency, at $3<z<5$, and it shows a little increasing trend with cosmic time at $z>3$. In addition, $M_{{\rm h}}^{{\rm pivot}}$ and its normalization are found to be almost unchanged during $0<z<5$. Our study shows an observational evidence that galaxy formation is ubiquitously most efficient near a halo mass of $M_{{\rm h}}\sim10^{12}M_{\odot}$ over cosmic time.
• ### The very wide-field $gzK$ galaxy survey -- II. The relationship between star-forming galaxies at $z \sim 2$ and their host haloes based upon HOD modelling(1602.01845)

Feb. 4, 2016 astro-ph.GA
We present the results of an halo occupation distribution (HOD) analysis of star-forming galaxies at $z \sim 2$. We obtained high-quality angular correlation functions based on a large sgzK sample, which enabled us to carry out the HOD analysis. The mean halo mass and the HOD mass parameters are found to increase monotonically with increasing $K$-band magnitude, suggesting that more luminous galaxies reside in more massive dark haloes. The luminosity dependence of the HOD mass parameters was found to be the same as in the local Universe; however, the masses were larger than in the local Universe over all ranges of magnitude. This implies that galaxies at $z \sim 2$ tend to form in more massive dark haloes than in the local Universe, a process known as downsizing. By analysing the dark halo mass evolution using the extended Press--Schechter formalism and the number evolution of satellite galaxies in a dark halo, we find that faint Lyman break galaxies at $z \sim 4$ could evolve into the faintest sgzKs $(22.0 < K \leq 23.0)$ at $z \sim 2$ and into the Milky-Way-like galaxies or elliptical galaxies in the local Universe, whereas the most luminous sgzKs $(18.0 \leq K \leq 21.0)$ could evolve into the most massive systems in the local Universe. The stellar-to-halo mass ratio (SHMR) of the sgzKs was found to be consistent with the prediction of the model, except that the SHMR of the faintest sgzKs was smaller than the prediction at $z \sim 2$. This discrepancy may be attributed that our samples are confined to star-forming galaxies.
• ### Probing cosmology with weak lensing selected clusters II: Dark energy and f(R) gravity models(1508.02104)

Nov. 4, 2015 astro-ph.CO
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.
• ### Probing cosmology with weak lensing selected clusters I: Halo approach and all-sky simulations(1504.05672)

Sept. 4, 2015 astro-ph.CO
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.
• ### Properties of Weak Lensing Clusters Detected on Hyper Suprime-Cam 2.3 Square Degree Field(1504.06974)

April 27, 2015 astro-ph.CO
• ### Detailed cluster lensing profiles at large radii and the impact on cluster weak lensing studies(1101.0650)

May 18, 2011 astro-ph.CO
Using a large set of ray-tracing in N-body simulations, we examine lensing profiles around massive dark haloes in detail, with a particular emphasis on the profile at around the virial radii. We compare radial convergence profiles, which are measured accurately in the ray-tracing simulations by stacking many dark haloes, with our simple analytic model predictions. Our analytic models consist of a main halo, which is modelled by the Navarro-Frenk-White (NFW) density profile with three different forms of the truncation, plus the correlated matter (2-halo term) around the main halo. We find that the smoothly truncated NFW profile best reproduces the simulated lensing profiles, out to more than 10 times the virial radius. We then use this analytic model to investigate potential biases in cluster weak lensing studies in which a single, untruncated NFW component is usually assumed in interpreting observed signals. We find that cluster masses, inferred by fitting reduced tangential shear profiles with the NFW profile, tend to be underestimated by ~5-10% if fitting is performed out to ~10'-30'. In contrast, the concentration parameter is overestimated typically by ~20% for the same fitting range. We also investigate biases in computing the signal-to-noise ratio of weak lensing mass peaks, finding them to be <4% for significant mass peaks. In the Appendices, we provide useful formulae for the smoothly truncated NFW profile.
• ### Simulations of Wide-Field Weak Lensing Surveys II: Covariance Matrix of Real Space Correlation Functions(1009.2558)

April 11, 2011 astro-ph.CO
Using 1000 ray-tracing simulations for a {\Lambda}-dominated cold dark model in Sato et al. (2009), we study the covariance matrix of cosmic shear correlation functions, which is the standard statistics used in the previous measurements. The shear correlation function of a particular separation angle is affected by Fourier modes over a wide range of multipoles, even beyond a survey area, which complicates the analysis of the covariance matrix. To overcome such obstacles we first construct Gaussian shear simulations from the 1000 realizations, and then use the Gaussian simulations to disentangle the Gaussian covariance contribution to the covariance matrix we measured from the original simulations. We found that an analytical formula of Gaussian covariance overestimates the covariance amplitudes due to an effect of finite survey area. Furthermore, the clean separation of the Gaussian covariance allows to examine the non-Gaussian covariance contributions as a function of separation angles and source redshifts. For upcoming surveys with typical source redshifts of z_s=0.6 and 1.0, the non-Gaussian contribution to the diagonal covariance components at 1 arcminute scales is greater than the Gaussian contribution by a factor of 20 and 10, respectively. Predictions based on the halo model qualitatively well reproduce the simulation results, however show a sizable disagreement in the covariance amplitudes. By combining these simulation results we develop a fitting formula to the covariance matrix for a survey with arbitrary area coverage, taking into account effects of the finiteness of survey area on the Gaussian covariance.
• ### Optimal filtering of optical and weak lensing data to search for galaxy clusters: application to the COSMOS field(1006.0610)

Dec. 16, 2010 astro-ph.CO
Galaxy clusters are usually detected in blind optical surveys via suitable filtering methods. We present an optimal matched filter which maximizes their signal-to-noise ratio by taking advantage of the knowledge we have of their intrinsic physical properties and of the data noise properties. In this paper we restrict our application to galaxy magnitudes, positions and photometric redshifts if available, and we also apply the filter separately to weak lensing data. The method is suitable to be naturally extended to a multi-band approach which could include not only additional optical bands but also observables with different nature such as X-rays. For each detection, the filter provides its significance, an estimate for the richness and for the redshift even if photo-z are not given. The provided analytical error estimate is tested against numerical simulations. We finally apply our method to the COSMOS field and compare the results with previous cluster detections obtained with different methods. Our catalogue contains 27 galaxy clusters with minimal threshold at 3-sigma level including both optical and weak-lensing information.
• ### MOIRCS Deep Survey. IX. Deep Near-Infrared Imaging Data and Source Catalog(1012.2115)

Dec. 13, 2010 astro-ph.CO
We present deep J-, H-, and Ks-band imaging data of the MOIRCS Deep Survey (MODS), which was carried out with Multi-Object Infrared Camera and Spectrograph (MOIRCS) mounted on the Subaru telescope in the GOODS-North region. The data reach 5sigma total limiting magnitudes for point sources of J=23.9, H=22.8, and Ks=22.8 (Vega magnitude) over 103 arcmin^2 (wide field). In 28 arcmin^2 of the survey area, which is ultra deep field of the MODS (deep field), the data reach the 5sigma depths of J=24.8, H=23.4, and Ks=23.8. The spatial resolutions of the combined images are FWHM ~ 0.6 arcsec and ~ 0.5 arcsec for the wide and deep fields in all bands, respectively. Combining the MODS data with the multi-wavelength public data taken with the HST, Spitzer, and other ground-based telescopes in the GOODS field, we construct a multi-wavelength photometric catalog of Ks-selected sources. Using the catalog, we present Ks-band number counts and near-infrared color distribution of the detected objects, and demonstrate some selection techniques with the NIR colors for high redshift galaxies. These data and catalog are publicly available via internet.
• ### Shear Power Spectrum Reconstruction using Pseudo-Spectrum Method(1004.3542)

Nov. 17, 2010 physics.data-an, astro-ph.CO
We develop a pseudo power spectrum technique for measuring the lensing power spectrum from weak lensing surveys in both the full sky and flat sky limits. The power spectrum approaches have a number of advantages over the traditional correlation function approach. We test the pseudo spectrum method by using numerical simulations with square-shape boundary that include masked regions with complex configuration due to bright stars and saturated spikes. Even when 25% of total area of the survey is masked, the method recovers the E-mode power spectrum at a sub-percent precision over a wide range of multipoles 100<l<10000. The systematic error is smaller than the statistical errors expected for a 2000 square degree survey. The residual B-mode spectrum is well suppressed in the amplitudes at less than a percent level relative to the E-mode. We also find that the correlated errors of binned power spectra caused by the survey geometry effects are not significant. Our method is applicable to the current and upcoming wide-field lensing surveys.