We constrain the scaling relation between optical richness (λ) and halo mass (M) for a sample of Sloan Digital Sky Survey (SDSS) red-sequence Matched-filter Probabilistic Percolation (redMaPPer) ...galaxy clusters within the context of the Planck cosmological model. We use a forward modeling approach where we model the probability distribution of optical richness for a given mass, . 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 8312 SDSS redMaPPer clusters with 20 ≤ λ ≤ 100 and 0.10 ≤ zλ ≤ 0.33 and show that the lognormal distribution model for , with four free parameters, well reproduces the measured abundances and lensing profiles simultaneously. The constraints are characterized by the mean relation, , with and (68% CL), where the pivot mass scale Mpivot = 3 × 1014 h−1 M , and the scatter with and . We find that a large scatter in halo masses is required at the lowest-richness bins (20 ≤ λ 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 ≤ λ ≤ 100.
ABSTRACT
We present a statistical weak-lensing magnification analysis on an optically selected sample of 3029 CAMIRA (Cluster finding Algorithm based on Multiband Identification of Red-sequence ...gAlaxies) galaxy clusters with richness N > 15 at redshift 0.2 ≤ z < 1.1 in the Subaru Hyper Suprime-Cam survey. We use two distinct populations of colour-selected, flux-limited background galaxies, namely the low-z and high-z samples at mean redshifts of ≈1.1 and ≈1.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 CAMIRA sample is detected at a significance level of 9.51σ, 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 CAMIRA sample with informative priors on other parameters. The resulting scaling relation is N∝ (M500)0.92 ± 0.13(1 + z)−0.48 ± 0.69, with a characteristic richness of N = 17.72 ± 2.60 and intrinsic lognormal scatter of 0.15 ± 0.07 at M500 = 1014 h−1 M⊙. With the derived N–M relation, we provide magnification-calibrated mass estimates of individual CAMIRA clusters, with the typical uncertainty of ≈39 and ≈32 per cent at richness of ≈20 and ≈40, respectively. We further compare our magnification-inferred N–M relation with those from the shear-based results in the literature, finding good agreement.
ABSTRACT We show that the projected number density profiles of Sloan Digital Sky Survey photometric galaxies around galaxy clusters display 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, , at fixed richness and redshift. The sample with smaller has a smaller ratio of the splashback radius to the traditional halo boundary than the subsample with larger , indicative of different mass accretion rates for these subsamples. The same subsamples 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 difference in the clustering amplitudes of these samples using cluster-photometric galaxy cross-correlations. This measurement is a clear indication that halo clustering depends on parameters other than halo mass. If is related to the mass assembly history of halos, the measurement is a manifestation of the halo assembly bias. However, our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than the predictions of collisionless Λ cold dark matter simulations. We show that dynamical friction, cluster mis-centering, or projection effects are not likely to be the sole source of these discrepancies. However, further investigations regarding unknown catastrophic weak lensing or cluster identification systematics are warranted.
The Hyper Suprime-Cam software pipeline Bosch, James; Armstrong, Robert; Bickerton, Steven ...
Publications of the Astronomical Society of Japan,
01/2018, Letnik:
70, Številka:
SP1
Journal Article
Recenzirano
Odprti dostop
Abstract
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.
Abstract
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 (corresponding to statistical scatters of ∼0.01 and ∼0.06 mag respectively). For synthetic galaxies with single-Sérsic profiles, forced CModel photometry achieves 13% photometric precision at i ∼ 20.0 mag and 18% precision at i ∼ 25.0 in the i band (corresponding to statistical scatters of ∼0.15 and ∼0.22 mag respectively). We show that both forced point spread function 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. Secondly, >20% of stars at 22.5 < i < 25.0 mag can be misclassified as extended objects. Thirdly, the current CModel algorithm tends to strongly underestimate the half-light radius and ellipticity of galaxy with i > 21.5 mag.
Abstract
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 deg2 split into six fields, with a mean i-band seeing of 0${^{\prime\prime}_{.}}$58 and 5σ point-source depth of i ∼ 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 ρ1 < 3 × 10−7 (requirement: 4 × 10−7) at 0${^{\circ}_{.}}$5 scales. A variety of galaxy shape-related null tests are statistically consistent with zero, but star–galaxy shape correlations reveal additive systematics on >1° 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.
Abstract
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 ≲ 1.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 ∼900 clusters with richness larger than 20 found in ∼140 deg2 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 (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% ± 4%) and the concentration of the profile (24% ± 11%). The level of cluster contamination can reach as high as ∼10% at R ≈ 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 ∼60 σ detection of the stacked CAMIRA surface mass density profile, with a mean mass of M200c = 1.67 ± 0.05(stat) × 1014 M⊙/h.
We perform an ensemble of N-body simulations with 20483 particles for 101 flat wCDM 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 autocorrelation as a function of halo masses, redshift, separations, and cosmological models based on principal component analysis and 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 autocorrelation, 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 concentration-mass 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 galaxy-galaxy weak lensing and the projected correlation function for any model within the wCDM cosmologies, in a few CPU seconds.
Abstract
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 ≤ zcl ≤ 1.0 with N ≥ 15 using the HSC first-year shear catalog covering ∼140 deg2. The exquisite depth and image quality of the HSC survey allow us to measure lensing signals around high-redshift clusters at 0.7 ≤ zcl ≤ 1.0 with a signal-to-noise ratio of 19 within the comoving radius range $0.5\lesssim R\lesssim 15\, h^{-1}\:{\rm Mpc}$. We constrain the richness–mass relations P(ln N ∣ M, z) of 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 proportional to ln (1 + z) and ln (1 + z)2 in 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 redshift. We show that constraints on the mean relation 〈M ∣ N〉 are consistent between the Planck and WMAP models, whereas the scatter values σ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.