Search

Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the 'renoising' noise bias correction approach introduced by Metacalibration. We implement this analytical noise bias correction within the AnaCal shear estimation framework and propose several enhancements to the noise bias correction algorithm. We evaluate the improved AnaCal using simulations designed to replicate Rubin LSST imaging data. These simulations feature semi-realistic galaxies and stars, complete with representative distributions of magnitudes and Galactic spatial density. We conduct tests under various observational challenges, including cosmic rays, defective CCD columns, bright star saturation, bleed trails, and spatially variable point spread functions. Our results indicate a multiplicative bias in weak lensing shear recovery of less than a few tenths of a percent, meeting LSST DESC requirements without requiring calibration from external image simulations. Additionally, our algorithm achieves rapid processing, handling one galaxy in less than a millisecond., Comment: 15 pages, 11 figures, submitted to MNRAS

The analyses of the next generation cosmological surveys demand an accurate, efficient, and differentiable method for simulating the universe and its observables across cosmological volumes. We present Hamiltonian ray tracing (HRT) -- the first post-Born (accounting for lens-lens coupling and without relying on the Born approximation), three-dimensional (without assuming the thin-lens approximation), and on-the-fly (applicable to any structure formation simulations) ray tracing algorithm based on the Hamiltonian formalism. HRT performs symplectic integration of the photon geodesics in a weak gravitational field, and can integrate tightly with any gravity solver, enabling co-evolution of matter particles and light rays with minimal additional computations. We implement HRT in the particle-mesh library $\texttt{pmwd}$, leveraging hardware accelerators such as GPUs and automatic differentiation capabilities based on $\texttt{JAX}$. When tested on a point-mass lens, HRT achieves sub-percent accuracy in deflection angles above the resolution limit across both weak and moderately strong lensing regimes. We also test HRT in cosmological simulations on the convergence maps and their power spectra., Comment: 26 pages, 8 figures; as published in JCAP

Recent cosmic shear analyses have exhibited inconsistencies of up to $1\sigma$ between the inferred cosmological parameters when analyzing summary statistics in real space versus harmonic space. In this paper, we demonstrate the consistent measurement and analysis of cosmic shear two-point functions in harmonic and real space using the $i${\sc Master} algorithm. This algorithm provides a consistent prescription to model the survey window effects and scale cuts in both real space (due to observational systematics) and harmonic space (due to model limitations), resulting in a consistent estimation of the cosmic shear power spectrum from both harmonic and real space estimators. We show that the $i$\textsc{Master} algorithm gives consistent results using measurements from the HSC Y1 mock shape catalogs in both real and harmonic space, resulting in consistent inferences of $S_8=\sigma_8(\Omega_m/0.3)^{0.5}$. This method provides an unbiased estimate of the cosmic shear power spectrum, and $S_8$ inference that has a correlation coefficient of 0.997 between analyses using measurements in real space and harmonic space. We observe the mean difference between the two inferred $S_8$ values to be 0.0004, far below the observed difference of 0.042 for the published HSC Y1 analyses and well below the statistical uncertainties. While the notation employed in this paper is specific to photometric galaxy surveys, the methods are equally applicable and can be extended to spectroscopic galaxy surveys, intensity mapping, and CMB surveys., Comment: 17 pages, 6 figures. For submission to MNRAS

The baryonic feedback effect is considered as a possible solution to the so-called $S_8$ tension indicated in cosmic shear cosmology. The baryonic effect is more significant on smaller scales, and affects the cosmic shear two-point correlation functions (2PCFs) with different scale- and redshift-dependencies from those of the cosmological parameters. In this paper, we use the Hyper Suprime-Cam Year 3 (HSC-Y3) data to measure the cosmic shear 2PCFs ($\xi_{\pm}$) down to 0.28 arcminutes, taking full advantage of the high number density of source galaxies in the deep HSC data, to explore a possible signature of the baryonic effect. While the published HSC analysis used the cosmic shear 2PCFs on angular scales, which are sensitive to the matter power spectrum at $k\lesssim 1~h{\rm Mpc}^{-1}$, the smaller scale HSC cosmic shear signal allows us to probe the signature of matter power spectrum up to $k\simeq 20~h{\rm Mpc}^{-1}$. Using the accurate emulator of the nonlinear matter power spectrum, DarkEmulator2, we show that the dark matter-only model can provide an acceptable fit to the HSC-Y3 2PCFs down to the smallest scales. In other words, we do not find any clear signature of the baryonic effects or do not find a systematic shift in the $S_8$ value with the inclusion of the smaller-scale information as would be expected if the baryonic effect is significant. Alternatively, we use a flexible 6-parameter model of the baryonic effects, which can lead to both enhancement and suppression in the matter power spectrum compared to the dark matter-only model, to perform the parameter inference of the HSC-Y3 2PCFs. We find that the small-scale HSC data allow only a fractional suppression of up to 5 percent in the matter power spectrum at $k\sim 1~h{\rm Mpc}^{-1}$, which is not sufficient to reconcile the $S_8$ tension., Comment: 30 pages, 16 figures

We aim to participate in the calibration of the X-ray photon count rate to halo mass scaling relation of galaxy clusters selected in the first eROSITA All-Sky Survey on the Western Galactic Hemisphere (eRASS1) using KiDS-1000 weak-lensing (WL) data. We measure the radial shear profiles around eRASS1 galaxy clusters using background galaxies in KiDS-1000, as well as the cluster member contamination. Furthermore we provide consistency checks with the other stage-III WL surveys who take part in the eRASS1 mass calibration, DES Y3 and HSC-Y3. We determine the cluster member contamination of eRASS1 clusters present in KiDS-1000 based on source number density profiles, where we account for the obscuration caused by cluster galaxies. The extracted shear profiles, together with the contamination model and the lens sample selection, are then analysed through a Bayesian population model. We calibrate the WL mass bias parameter by analysing realistic synthetic shear profiles from mock cluster catalogues. Our consistency checks between KiDS-1000 and DES Y3 & HSC-Y3 include the comparison of contamination-corrected density contrast profiles employing the union of background sources around common clusters, as well as the individual scaling relation results. We present a global contamination model for eRASS1 clusters in KiDS-1000 and the calibration results of the X-ray photon count rate to halo mass relation. The results of the WL mass bias parameter show that the uncertainty of the multiplicative shear bias dominates the systematic error budget at low clusters redshifts while the uncertainty of our contamination model does at high ones. The cross-checks between the three WL surveys show that they are statistically consistent with each other. This enables for the first time cosmological constraints from clusters calibrated by three state-of-the-art WL surveys. (abridged), Comment: 22 pages, 23 figures

We present $\texttt{Miko}$, a catalog-to-cosmology pipeline for general flat-sky field-level inference, which provides access to cosmological information beyond the two-point statistics. In the context of weak lensing, we identify several new field-level analysis systematics (such as aliasing, Fourier mode-coupling, and density-induced shape noise), quantify their impact on cosmological constraints, and correct the biases to a percent level. Next, we find that model misspecification can lead to both absolute bias and incorrect uncertainty quantification for the inferred cosmological parameters in realistic simulations. The Gaussian map prior infers unbiased cosmological parameters, regardless of the true data distribution, but it yields overconfident uncertainties. The log-normal map prior quantifies the uncertainties accurately, although it requires careful calibration of the shift parameters for unbiased cosmological parameters. We demonstrate systematics control down to the $2\%$ level for both models, making them suitable for ongoing weak lensing surveys., Comment: 18 pages, 19 figures; see figures 12 and 17 for the key results; published in PRD

Continuing work presented in Li et al. (2021), we performed a series of tests to our high-resolution three-dimensional mass map reconstruction algorithm \splinv{}. We test the mass reconstruction accuracy against realistic mock catalogs generated using shear field produced by triaxial halos with the inner density profile of $\rho \propto r^{-1}$ and of $\rho \propto r^{-1.5}$. The galaxy shape noise is modeled based on the Year-1 Subaru Hyper Suprime-Cam (HSC) Survey. After reviewing mathematical details of our algorithm and dark matter halo models, we determine an optimal value of the coefficient of the adaptive LASSO regression penalty term for single halo reconstruction. We successfully measure halo masses for massive triaxial halos; the mass determination accuracy is 5 percent for halos with $M = 10^{14.6}~M_\odot$ at $0.0625\leq z \leq 0.2425$, and 5 percent for those with $10^{14.8}~M_\odot$ at $0.0625\leq z \leq 0.4675$, and 20 percent for $M= 10^{15.0} ~M_\odot$ and $M=10^{15.2}~M_\odot$ in the redshift range $0.0625\leq z \leq 0.4675$. The redshift estimate accuracy is consistently below $\Delta z /z \leq 0.05$ for the above halo masses in the range $0.1525\leq z \leq 0.4675$. We also demonstrate that the orientation of triaxial halos and systematic error in our halo model do not affect reconstruction result significantly. Finally, we present results from reconstruction of mass distribution using shear catalogs produced by multiple halos, to show \splinv{}'s capability using realistic shear maps from ongoing and future galaxy surveys., Comment: 20 pages, 20 figures

We present cosmology results obtained from a blind joint analysis of the abundance, projected clustering, and weak lensing of galaxy clusters measured from the Sloan Digital Sky Survey (SDSS) redMaPPer cluster catalog and the Hyper-Suprime Cam (HSC) Year3 shape catalog. We present a full-forward model for the cluster observables, which includes empirical modeling for the anisotropic boosts on the lensing and clustering signals of optical clusters. We validate our analysis via mock cluster catalogs which include observational systematics, such as the projection effect and the effect of baryonic feedback, and find that our analysis can robustly constrain cosmological parameters in an unbiased manner without any informative priors on our model parameters. The joint analysis of our observables in the context of the flat $\Lambda$CDM model results in cosmological constraints for $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m} / 0.3}=0.816^{+0.041}_{-0.039}$. Our result is consistent with the $S_8$ inference from other cosmic microwave background- and large scale structure-based cosmology analyses, including the result from the \emph{Planck} 2018 primary CMB analysis., Comment: v1: 22 pages, 15 figures, Comments welcome

Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-percent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (FPFS), using JAX, which is publicly available as part of a new suite of analytic shear algorithms called AnaCal. This code can analytically calibrate the shear response of any nonlinear observable constructed with the FPFS shapelets and detection modes utilizing auto-differentiation (AD), generalizing the formalism to include a family of shear estimators with corrections for detection and selection biases. Using the AD capability of JAX, it calculates the full Hessian matrix of the non-linear observables, which improves the previously presented second-order noise bias correction in the shear estimation. As an illustration of the power of the new AnaCal framework, we optimize the effective galaxy number density in the space of the generalized shear estimators using an LSST-like galaxy image simulation for the ten-year LSST. For the generic shear estimator, the magnitude of the multiplicative bias $|m|$ is below $3\times 10^{-3}$ (99.7% confidence interval), and the effective galaxy number density is improved by 5%. We also discuss some planned future additions to the AnaCal software suite to extend its applicability beyond the FPFS measurements., Comment: 9 pages, 7 figures, Accepted for publication in MNRAS

We present cosmological parameter constraints from a blind joint analysis of three two-point correlation functions measured from the Year 3 Hyper Suprime-Cam (HSC-Y3) imaging data, covering 416 deg$^2$, and the SDSS DR11 spectroscopic galaxies spanning the redshift range $[0.15, 0.70]$. We subdivide the SDSS galaxies into three volume-limited samples separated in redshift, each of which acts as a large-scale structure tracer characterized by the measurement of the projected correlation function, $w_{\rm p}(R)$. We also use the measurements of the galaxy-galaxy weak lensing signal $\Delta \Sigma(R)$ for each of these SDSS samples which act as lenses for a secure sample of source galaxies selected from the HSC-Y3 shape catalog based on their photometric redshifts. We combine these measurements with the cosmic shear correlation functions, $\xi_{\pm}(\vartheta)$, measured for our HSC source sample. We model these observables with the minimal bias model of the galaxy clustering observables in the context of a flat $\Lambda$CDM cosmology. We use conservative scale cuts, $R>12$ and $8~h^{-1}$Mpc, for $\Delta\Sigma$ and $w_{\rm p}$, respectively, where the minimal bias model is valid, in addition to conservative prior on the residual bias in the mean redshift of the HSC photometric source galaxies. Our baseline analysis yields $S_8=0.775^{+0.043}_{-0.038}$ (68% C.I.) for the $\Lambda$CDM model, after marginalizing over uncertainties in other parameters. Our value of $S_8$ is consistent with that from the Planck 2018 data, but the credible interval of our result is still relatively large. Our results are statistically consistent with those of a companion paper, which extends this analysis to smaller scales with an emulator-based halo model., Comment: 29 pages, 17 figures, 4 tables, version accepted for publication in PRD together with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

We present cosmology results from a blinded joint analysis of cosmic shear, $\xi_{\pm}(\vartheta)$, galaxy-galaxy weak lensing, $\Delta\!\Sigma(R)$, and projected galaxy clustering, $w_{\rm p}(R)$, measured from the Hyper Suprime-Cam three-year (HSC-Y3) shape catalog and the Sloan Digital Sky Survey (SDSS) DR11 spectroscopic galaxy catalog - a 3$\times$2pt cosmology analysis. We define luminosity-cut samples of SDSS galaxies to serve as the tracers of $w_{\rm p}$ and as the lens samples for $\Delta\!\Sigma$ in three spectroscopic redshift bins spanning the range $0.15

We use the Sloan Digital Sky Survey (SDSS) BOSS galaxies and their overlap with approximately 416 sq. degree of deep $grizy$-band imaging from the Subaru Hyper Suprime-Cam Survey (HSC). We measure three two-point correlations that form the basis of the cosmological inference presented in our companion papers, Miyatake et al. and Sugiyama et al. We use three approximately volume limited subsamples of spectroscopic galaxies by their $i$-band magnitude from the SDSS-BOSS: LOWZ (0.1

We measure weak lensing cosmic shear power spectra from the three-year galaxy shear catalog of the Hyper Suprime-Cam (HSC) Subaru Strategic Program imaging survey. The shear catalog covers $416 \ \mathrm{deg}^2$ of the northern sky, with a mean $i$-band seeing of 0.59 arcsec and an effective galaxy number density of 15 $\mathrm{arcmin}^{-2}$ within our adopted redshift range. With an $i$-band magnitude limit of 24.5 mag, and four tomographic redshift bins spanning $0.3 \leq z_{\mathrm{ph}} \leq 1.5$ based on photometric redshifts, we obtain a high-significance measurement of the cosmic shear power spectra, with a signal-to-noise ratio of approximately 26.4 in the multipole range $300<\ell<1800$. The accuracy of our power spectrum measurement is tested against realistic mock shear catalogs, and we use these catalogs to get a reliable measurement of the covariance of the power spectrum measurements. We use a robust blinding procedure to avoid confirmation bias, and model various uncertainties and sources of bias in our analysis, including point spread function systematics, redshift distribution uncertainties, the intrinsic alignment of galaxies and the modeling of the matter power spectrum. For a flat $\Lambda$CDM model, we find $S_8 \equiv \sigma_8 (\Omega_m/0.3)^{0.5} =0.776^{+0.032}_{-0.033}$, which is in excellent agreement with the constraints from the other HSC Year 3 cosmology analyses, as well as those from a number of other cosmic shear experiments. This result implies a $\sim$$2\sigma$-level tension with the Planck 2018 cosmology. We study the effect that various systematic errors and modeling choices could have on this value, and find that they can shift the best-fit value of $S_8$ by no more than $\sim$$0.5\sigma$, indicating that our result is robust to such systematics., Comment: 35 pages, 18 figures, 6 tables, for coordinated submission to PRD with other HSC Y3 weak lensing cosmology papers - see https://hsc-release.mtk.nao.ac.jp/doc/index.php/wly3/

We present cosmological constraints from the Subaru Hyper Suprime-Cam (HSC) first-year weak lensing shear catalogue using convolutional neural networks (CNNs) and conventional summary statistics. We crop 19 $3\times3\,\mathrm{{deg}^2}$ sub-fields from the first-year area, divide the galaxies with redshift $0.3\le z\le1.5$ into four equally-spaced redshift bins, and perform tomographic analyses. We develop a pipeline to generate simulated convergence maps from cosmological $N$-body simulations, where we account for effects such as intrinsic alignments (IAs), baryons, photometric redshift errors, and point spread function errors, to match characteristics of the real catalogue. We train CNNs that can predict the underlying parameters from the simulated maps, and we use them to construct likelihood functions for Bayesian analyses. In the $\Lambda$ cold dark matter model with two free cosmological parameters $\Omega_\mathrm{m}$ and $\sigma_8$, we find $\Omega_\mathrm{m}=0.278_{-0.035}^{+0.037}$, $S_8\equiv(\Omega_\mathrm{m}/0.3)^{0.5}\sigma_8=0.793_{-0.018}^{+0.017}$, and the IA amplitude $A_\mathrm{IA}=0.20_{-0.58}^{+0.55}$. In a model with four additional free baryonic parameters, we find $\Omega_\mathrm{m}=0.268_{-0.036}^{+0.040}$, $S_8=0.819_{-0.024}^{+0.034}$, and $A_\mathrm{IA}=-0.16_{-0.58}^{+0.59}$, with the baryonic parameters not being well-constrained. We also find that statistical uncertainties of the parameters by the CNNs are smaller than those from the power spectrum (5--24 percent smaller for $S_8$ and a factor of 2.5--3.0 smaller for $\Omega_\mathrm{m}$), showing the effectiveness of CNNs for uncovering additional cosmological information from the HSC data. With baryons, the $S_8$ discrepancy between HSC first-year data and Planck 2018 is reduced from $\sim2.2\,\sigma$ to $0.3\text{--}0.5\,\sigma$., Comment: 22 pages, 14 figures

Cosmological weak lensing measurements rely on a precise measurement of the shear two-point correlation function (2PCF) along with a deep understanding of systematics that affect it. In this work, we demonstrate a general framework for detecting and modeling the impact of PSF systematics on the cosmic shear 2PCF, and mitigating its impact on cosmological analysis. Our framework can describe leakage and modeling error from all spin-2 quantities contributed by the PSF second and higher moments, rather than just the second moments, using the cross-correlations between galaxy shapes and PSF moments. We interpret null tests using the HSC Year 3 (Y3) catalogs with this formalism, and find that leakage from the spin-2 combination of PSF fourth moments is the leading contributor to additive shear systematics, with total contamination that is an order of magnitude higher than that contributed by PSF second moments alone. We conducted a mock cosmic shear analysis for HSC Y3, and find that, if uncorrected, PSF systematics can bias the cosmological parameters $\Omega_m$ and $S_8$ by $\sim$0.3$\sigma$. The traditional second moment-based model can only correct for a 0.1$\sigma$ bias, leaving the contamination largely uncorrected. We conclude it is necessary to model both PSF second and fourth moment contamination for HSC Y3 cosmic shear analysis. We also reanalyze the HSC Y1 cosmic shear analysis with our updated systematics model, and identify a 0.07$\sigma$ bias on $\Omega_m$ when using the more restricted second moment model from the original analysis. We demonstrate how to self-consistently use the method in both real space and Fourier space, assess shear systematics in tomographic bins, and test for PSF model overfitting., Comment: 33 pages, 28 figures, match the refereed version in MNRAS

We present posterior sample redshift distributions for the Hyper Suprime-Cam Subaru Strategic Program Weak Lensing three-year (HSC Y3) analysis. Using the galaxies' photometry and spatial cross-correlations, we conduct a combined Bayesian Hierarchical Inference of the sample redshift distributions. The spatial cross-correlations are derived using a subsample of Luminous Red Galaxies (LRGs) with accurate redshift information available up to a photometric redshift of $z < 1.2$. We derive the photometry-based constraints using a combination of two empirical techniques calibrated on spectroscopic- and multiband photometric data that covers a spatial subset of the shear catalog. The limited spatial coverage induces a cosmic variance error budget that we include in the inference. Our cross-correlation analysis models the photometric redshift error of the LRGs to correct for systematic biases and statistical uncertainties. We demonstrate consistency between the sample redshift distributions derived using the spatial cross-correlations, the photometry, and the posterior of the combined analysis. Based on this assessment, we recommend conservative priors for sample redshift distributions of tomographic bins used in the three-year cosmological Weak Lensing analyses., Comment: 23 pages, 11 figures, 1 table, updated to match the accepted version in the MNRAS

Shear estimation bias from galaxy detection and blending identification is now recognized as an issue for ongoing and future weak lensing surveys. Currently, the empirical approach to correcting for this bias involves numerically shearing every observed galaxy and rerunning the detection and selection process. In this work, we provide an analytical correction for this bias that is accurate to subpercent level and far simpler to use. With the interpretation that smoothed image pixel values and galaxy properties are projections of the image signal onto a set of basis functions, we analytically derive the linear shear responses of both the pixel values and the galaxy properties (i.e., magnitude, size and shape) using the shear responses of the basis functions. With these derived shear responses, we correct for biases from shear-dependent galaxy detection and galaxy sample selection. With the analytical covariance matrix of measurement errors caused by image noise on pixel values and galaxy properties, we correct for the noise biases in galaxy shape measurement and the detection/selection process to the second-order in noise. The code used for this paper can carry out the detection, selection, and shear measurement for ~1000 galaxies per CPU second. The algorithm is tested with realistic image simulations, and we find, after the analytical correction (without relying on external image calibration) for the detection/selection bias of about $-4\%$, the multiplicative shear bias is $-0.12 \pm 0.10\%$ for isolated galaxies; and $-0.3 \pm 0.1\%$ for blended galaxies with Hyper Suprime-Cam observational condition., Comment: 23 pages, 22 figures, Accepted in MNRAS

Recovering credible cosmological parameter constraints in a weak lensing shear analysis requires an accurate model that can be used to marginalize over nuisance parameters describing potential sources of systematic uncertainty, such as the uncertainties on the sample redshift distribution $n(z)$. Due to the challenge of running Markov Chain Monte-Carlo (MCMC) in the high dimensional parameter spaces in which the $n(z)$ uncertainties may be parameterized, it is common practice to simplify the $n(z)$ parameterization or combine MCMC chains that each have a fixed $n(z)$ resampled from the $n(z)$ uncertainties. In this work, we propose a statistically-principled Bayesian resampling approach for marginalizing over the $n(z)$ uncertainty using multiple MCMC chains. We self-consistently compare the new method to existing ones from the literature in the context of a forecasted cosmic shear analysis for the HSC three-year shape catalog, and find that these methods recover similar cosmological parameter constraints, implying that using the most computationally efficient of the approaches is appropriate. However, we find that for datasets with the constraining power of the full HSC survey dataset (and, by implication, those upcoming surveys with even tighter constraints), the choice of method for marginalizing over $n(z)$ uncertainty among the several methods from the literature may significantly impact the statistical uncertainties on cosmological parameters, and a careful model selection is needed to ensure credible parameter intervals., Comment: 16 pages, 9 figures, Accepted by MNRAS

We assess the dark matter halo masses of luminous AGNs over the redshift range 0.2 to 1.2 using galaxy-galaxy lensing based on imaging data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We measure the weak lensing signal of a sample of 8882 AGNs constructed using HSC and WISE photometry. The lensing detection around AGNs has a signal-to-noise ratio of 15. As expected, we find that the lensing mass profile is consistent with that of massive galaxies ($M_{*}\sim 10.8~M_\odot$). Surprisingly, the lensing signal remains unchanged when the AGN sample is split into low and high stellar mass hosts. Specifically, we find that the excess surface density (ESD) of AGNs, residing in galaxies with high stellar masses, significantly differs from that of the control sample. We further fit a halo occupation distribution model to the data to infer the posterior distribution of parameters including the average halo mass. We find that the characteristic halo mass of the full AGN population lies near the knee ($\rm log(M_h/h^{-1}M_{\odot})=12.0$) of the stellar-to-halo mass relation (SHMR). Illustrative of the results given above, the halo masses of AGNs residing in host galaxies with high stellar masses (i.e., above the knee of the SHMR) falls below the calibrated SHMR while the halo mass of the low stellar mass sample is more consistent with the established SHMR. These results indicate that massive halos with higher clustering bias tend to suppress AGN activity, probably due to the lack of available gas., Comment: 17 pages,10 figures, accepted for publication in the Astrophysical Journal

We present cosmological parameter constraints from a blinded joint analysis of galaxy-galaxy weak lensing, $\Delta\!\Sigma(R)$, and projected correlation function, $w_\mathrm{p}(R)$, measured from the first-year HSC (HSC-Y1) data and SDSS spectroscopic galaxies over $0.150.75$ for the $\Delta\!\Sigma$ measurements, selected based on their photometric redshifts. For theoretical template, we use the "minimal bias" model for the cosmological clustering observables for the flat $\Lambda$CDM cosmological model. We compare the model predictions with the measurements in each redshift bin on large scales, $R>12$ and $8~h^{-1}\mathrm{Mpc}$ for $\Delta\!\Sigma(R)$ and $w_\mathrm{p}(R)$, respectively, where the perturbation theory-inspired model is valid. When we employ weak priors on cosmological parameters, without CMB information, we find $S_8=0.936^{+0.092}_{-0.086}$, $\sigma_8=0.85^{+0.16}_{-0.11}$, and $\Omega_\mathrm{m}=0.283^{+0.12}_{-0.035}$ for the flat $\Lambda$CDM model. Although the central value of $S_8$ appears to be larger than those inferred from other cosmological experiments, we find that the difference is consistent with expected differences due to sample variance, and our results are consistent with the other results to within the statistical uncertainties. (abriged), Comment: 24 pages, 19 figures, 4 tables, to be submitted to Phys. Rev. D

We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing ($\Delta\!\Sigma$) and projected galaxy clustering ($w_{\rm p}$) measured from the Hyper Suprime-Cam Year-1 (HSC-Y1) data and spectroscopic Sloan Digital Sky Survey (SDSS) galaxy catalogs in the redshift range $0.15

Dedicated 'Stage IV' observatories will soon observe the entire extragalactic sky, to measure the 'cosmic shear' distortion of galaxy shapes by weak gravitational lensing. To measure the apparent shapes of those galaxies, we present an improved version of the Fourier Power Function Shapelets (FPFS) shear measurement method. This now includes analytic corrections for sources of bias that plague all shape measurement algorithms: including noise bias (due to noise in nonlinear combinations of observable quantities) and selection bias (due to sheared galaxies being more or less likely to be detected). Crucially, these analytic solutions do not rely on calibration from external image simulations. For isolated galaxies, the small residual $\sim$$10^{-3}$ multiplicative bias and $\lesssim$$10^{-4}$ additive bias now meet science requirements for Stage IV experiments. FPFS also works accurately for faint galaxies and robustly against stellar contamination. Future work will focus on deblending overlapping galaxies. The code used for this paper can process $>$$1000$ galaxy images per CPU second and is available from https://github.com/mr-superonion/FPFS., Comment: 12 pages, 8 figures, 3 tables; MNRAS (in press)

The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5\sigma$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/., Comment: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

We present the first weak-lensing mass calibration and X-ray scaling relations of galaxy clusters and groups selected in the $eROSITA$ Final Equatorial Depth Survey (eFEDS) observed by Spectrum Roentgen Gamma/$eROSITA$ over a contiguous footprint with an area of $\approx140$ deg$^2$, using the three-year (S19A) weak-lensing data from the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. In this work, a sample of $434$ optically confirmed galaxy clusters (and groups) at redshift $0.01\lesssim z \lesssim1.3$ with a median of $0.35$ is studied, of which $313$ systems are uniformly covered by the HSC survey to enable the extraction of the weak-lensing shear observable. In a Bayesian population modelling, we perform a blind analysis for the weak-lensing mass calibration by simultaneously modelling the observed count rate $\eta$ and the shear profile $g$ of individual clusters through the count rate-to-mass-and-redshift ($\eta$--$M_{500}$--$z$) and weak-lensing mass-to-mass-and-redshift ($M_{\mathrm{WL}}$--$M_{500}$--$z$) relations, respectively, while accounting for the bias in these observables using simulation-based calibrations. As a result, the count rate-inferred and lensing-calibrated cluster mass is obtained from the joint modelling of the scaling relations, as the ensemble mass spanning a range of $10^{13}h^{-1}M_{\odot}\lesssim M_{500}\lesssim10^{15} h^{-1}M_{\odot}$ with a median of $\approx10^{14} h^{-1}M_{\odot}$ for the eFEDS sample. With the mass calibration, we further model the X-ray observable-to-mass-and-redshift relations, including the rest-frame soft-band and bolometric luminosity ($L_{\mathrm{X}}$ and $L_{\mathrm{b}}$), the emission-weighted temperature $T_{\mathrm{X}}$, the mass of intra-cluster medium $M_{\mathrm{g}}$, and the mass proxy $Y_{\mathrm{X}}$, which is the product of $T_{\mathrm{X}}$ and $M_{\mathrm{g}}$. (abridged), Comment: Accepted for publication in A&A. Submitted to A&A for the Special Issue: The Early Data Release of eROSITA and Mikhail Pavlinsky ART-XC on the SRG Mission

We present the galaxy shear catalog that will be used for the three-year cosmological weak gravitational lensing analyses using data from the Wide layer of the Hyper Suprime-Cam (HSC) Subaru Strategic Program (SSP) Survey. The galaxy shapes are measured from the $i$-band imaging data acquired from 2014 to 2019 and calibrated with image simulations that resemble the observing conditions of the survey based on training galaxy images from the Hubble Space Telescope in the COSMOS region. The catalog covers an area of 433.48 deg$^2$ of the northern sky, split into six fields. The mean $i$-band seeing is 0.59 arcsec. With conservative galaxy selection criteria (e.g., $i$-band magnitude brighter than 24.5), the observed raw galaxy number density is 22.9 arcmin$^{-2}$, and the effective galaxy number density is 19.9 arcmin$^{-2}$. The calibration removes the galaxy property-dependent shear estimation bias to a level: $|\delta m|<9\times 10^{-3}$. The bias residual $\delta m$ shows no dependence on redshift in the range $01$ degree scales and will require mitigation during the inference of cosmological parameters using cosmic shear measurements., Comment: 39 pages, 31 figures, 5 tables; PASJ (in press)

Using photometric galaxies from the HSC survey, we measure the stellar mass density profiles for satellite galaxies as a function of the projected distance, $r_p$, to isolated central galaxies (ICGs) selected from SDSS/DR7 spectroscopic galaxies at $z\sim0.1$. By stacking HSC images, we also measure the projected stellar mass density profiles for ICGs and their stellar halos. The total mass distributions are further measured from HSC weak lensing signals. ICGs dominate within $\sim$0.15 times the halo virial radius ($0.15R_{200}$). The stellar mass versus total mass fractions drop with the increase in $r_p$ up to $\sim0.15R_{200}$, beyond which they are less than 1\% while stay almost constant, indicating the radial distribution of satellites trace dark matter. The total stellar mass in satellites is proportional to the virial mass of the host halo, $M_{200}$, for ICGs more massive than $10^{10.5}M_\odot$, i.e., $M_{\ast,\mathrm{sat}} \propto M_{200}$, whereas the relation between the stellar mass of ICGs $+$ stellar halos and $M_{200}$ is close to $M_{\ast,\mathrm{ICG+diffuse}}\propto M_{200}^{1/2}$. Below $10^{10.5}M_\odot$, the change in $M_{200}$ is much slower with the decrease in $M_{\ast,\mathrm{ICG+diffuse}}$. At fixed stellar mass, red ICGs are hosted by more massive dark matter halos and have more satellites. At $M_{200}\sim10^{12.7}M_\odot$, both $M_{\ast,\mathrm{sat}}$ and the fraction of stellar mass in satellites versus total stellar mass, $f_\mathrm{sat}$, tend to be slightly higher around blue ICGs, perhaps implying the late formation of blue galaxies. $f_\mathrm{sat}$ increases with the increase in both $M_{\ast,\mathrm{ICG+diffuse}}$ and $M_{200}$, and scales more linearly with $M_{200}$. We provide best-fitting formulas for these scaling relations and for red and blue ICGs separately., Comment: accepted by ApJ - comments welcome - data available upon request

We use the Hyper Suprime-Cam Subaru Strategic Program S19A shape catalog to construct weak lensing shear-selected cluster samples. From aperture mass maps covering $\sim 510$~deg$^2$ created using a truncated Gaussian filter, we construct a catalog of 187 shear-selected clusters that correspond to mass map peaks with the signal-to-noise ratio larger than 4.7. Most of the shear-selected clusters have counterparts in optically-selected clusters, from which we estimate the purity of the catalog to be higher than 95\%. The sample can be expanded to 418 shear-selected clusters with the same signal-to-noise ratio cut by optimizing the shape of the filter function and by combining weak lensing mass maps created with several different background galaxy selections. We argue that dilution and obscuration effects of cluster member galaxies can be mitigated by using background source galaxy samples and adopting the filter function with its inner boundary larger than about $2'$. The large samples of shear-selected clusters that are selected without relying on any baryonic tracer are useful for detailed studies of cluster astrophysics and cosmology., Comment: 13 pages, 9 figures, accepted for publication in PASJ, cluster catalogs are available at https://github.com/oguri/cluster_catalogs

We propose a novel method to reconstruct high-resolution three-dimensional mass maps using data from photometric weak-lensing surveys. We apply an adaptive LASSO algorithm to perform a sparsity-based reconstruction on the assumption that the underlying cosmic density field is represented by a sum of Navarro-Frenk-White halos. We generate realistic mock galaxy shape catalogues by considering the shear distortions from isolated halos for the configurations matched to Subaru Hyper Suprime-Cam Survey with its photometric redshift estimates. We show that the adaptive method significantly reduces line-of-sight smearing that is caused by the correlation between the lensing kernels at different redshifts. Lensing clusters with lower mass limits of $10^{14.0} h^{-1}M_{\odot}$, $10^{14.7} h^{-1}M_{\odot}$, $10^{15.0} h^{-1}M_{\odot}$ can be detected with 1.5-$\sigma$ confidence at the low ($z<0.3$), median ($0.3\leq z< 0.6$) and high ($0.6\leq z< 0.85$) redshifts, respectively, with an average false detection rate of 0.022 deg$^{-2}$. The estimated redshifts of the detected clusters are systematically lower than the true values by $\Delta z \sim 0.03$ for halos at $z\leq 0.4$, but the relative redshift bias is below $0.5\%$ for clusters at $0.4

We conduct a comprehensive and statistical study of the luminosity functions (LFs) for satellite galaxies, by counting photometric galaxies from HSC, DECaLS and SDSS around isolated central galaxies (ICGs) and paired galaxies from the SDSS/DR7 spectroscopic sample. Results of different surveys show very good agreement. The satellite LFs can be measured down to $M_V\sim-10$, and for central primary galaxies as small as $8.5<\log_{10}M_\ast/M_\odot<9.2$ and $9.2<\log_{10}M_\ast/M_\odot<9.9$, implying there are on average 3--8 satellites with $M_V<-10$ around LMC-mass ICGs. The bright end cutoff of satellite LFs and the satellite abundance are both sensitive to the magnitude gap between the primary and its companions, indicating galaxy systems with larger magnitude gaps are on average hosted by less massive dark matter haloes. By selecting primaries with stellar mass similar to our MW, we discovered that i) the averaged satellite LFs of ICGs with different magnitude gaps to their companions and of galaxy pairs with different colour or colour combinations all show steeper slopes than the MW satellite LF; ii) there are on average more satellites with $-15-12$. Hence the MW and its satellite system are statistically atypical of our sample of MW-mass systems. In consequence, our MW is not a good representative of other MW-mass galaxies. Strong cosmological implications based on only MW satellites await additional discoveries of fainter satellites in extra-galactic systems. Interestingly, the MW satellite LF is typical among other MW-mass systems within 40~Mpc in the local Universe, perhaps implying the Local Volume is an under-dense region., Comment: accepted by MNRAS, with more detailed discussions about the scatter

We apply the Fourier Power Function Shapelets (FPFS) shear estimator to the first year data of the Hyper Suprime-Cam survey to construct a shape catalog. The FPFS shear estimator has been demonstrated to have multiplicative bias less than $1\%$ in the absence of blending, regardless of complexities of galaxy shapes, smears of point spread functions (PSFs) and contamination from noise. The blending bias is calibrated with realistic image simulations, which include the impact of neighboring objects, using the COSMOS Hubble Space Telescope images. Here we carefully test the influence of PSF model residual on the FPFS shear estimation and the uncertainties in the shear calibration. Internal null tests are conducted to characterize potential systematics in the FPFS shape catalog and the results are compared with those measured using a catalog where the shapes were estimated using the re-Gaussianization algorithms. Furthermore, we compare various weak lensing measurements between the FPFS shape catalog and the re-Gaussianization shape catalog and conclude that the weak lensing measurements between these two shape catalogs are consistent with each other within the statistical uncertainty., Comment: 24 pages, 19 figures, accepted by ApJS

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/., Comment: 26 pages, 23 figures, accepted for publication in PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

We study the faint stellar halo of isolated central galaxies, by stacking galaxy images in the HSC survey and accounting for the residual sky background sampled with random points. The surface brightness profiles in HSC $r$-band are measured for a wide range of galaxy stellar masses ($9.2<\log_{10}M_\ast/M_\odot<11.4$) and out to 120 kpc. Failing to account for the stellar halo below the noise level of individual images will lead to underestimates of the total luminosity by $\leq 15\%$. Splitting galaxies according to the concentration parameter of their light distributions, we find that the surface brightness profiles of low concentration galaxies drop faster between 20 and 100 kpc than those of high concentration galaxies. Albeit the large galaxy-to-galaxy scatter, we find a strong self-similarity of the stellar halo profiles. They show unified forms once the projected distance is scaled by the halo virial radius. The colour of galaxies is redder in the centre and bluer outside, with high concentration galaxies having redder and more flattened colour profiles. There are indications of a colour minimum, beyond which the colour of the outer stellar halo turns red again. This colour minimum, however, is very sensitive to the completeness in masking satellite galaxies. We also examine the effect of the extended PSF in the measurement of the stellar halo, which is particularly important for low mass or low concentration galaxies. The PSF-corrected surface brightness profile can be measured down to $\sim$31 $\mathrm{mag}/\mathrm{arcsec}^2$ at 3-$\sigma$ significance. PSF also slightly flattens the measured colour profiles., Comment: accepted by MNRAS - Significant changes have been made compared with the first version, including discussions on the extended PSF wings, robustness of our results to source detection and masking thresholds and more detailed investigations on the indications of positive colour gradients

The tilt, rotation, or offset of each CCD with respect to the focal plane, as well as the distortion of the focal plane itself, cause shape distortions to the observed objects, an effect typically known as field distortion (FD). We point out that FD provides a unique way of quantifying the accuracy of cosmic shear measurement. The idea is to stack the shear estimators from galaxies that share similar FD-induced shape distortions. Given that the latter can be calculated with parameters from astrometric calibrations, the accuracy of the shear estimator can be directly tested on real images. It provides a way to calibrate the multiplicative and additive shear recovery biases within the scientific data itself, without requiring simulations or any external data sets. We use the CFHTLenS images to demonstrate the accuracy of the Fourier_Quad shear recovery method. We highlight some details in our image processing pipeline, including background removal, source identification and deblending, astrometric calibration, star selection for PSF reconstruction, noise reduction, etc.. We show that in the shear ranges of -0.005 < g_1 < 0.005 and -0.008 < g_2 < 0.008, the multiplicative biases are at the level of < 0.04. Slight additive biases on the order of 5E-4 (6 sigma) are identified for sources provided by the official CFHTLenS catalog (not using its shear catalog), but are minor (4 sigma) for source catalog generated by our Fourier_Quad pipeline., Comment: 11 pages, 12 figures, ApJ in press

We reinterpret the shear estimator developed by Zhang & Komatsu (2011) within the framework of Shapelets and propose the Fourier Power Function Shapelets (FPFS) shear estimator. Four shapelet modes are calculated from the power function of every galaxy's Fourier transform after deconvolving the Point Spread Function (PSF) in Fourier space. We propose a novel normalization scheme to construct dimensionless ellipticity and its corresponding shear responsivity using these shapelet modes. Shear is measured in a conventional way by averaging the ellipticities and responsivities over a large ensemble of galaxies. With the introduction and tuning of a weighting parameter, noise bias is reduced below one percent of the shear signal. We also provide an iterative method to reduce selection bias. The FPFS estimator is developed without any assumption on galaxy morphology, nor any approximation for PSF correction. Moreover, our method does not rely on heavy image manipulations nor complicated statistical procedures. We test the FPFS shear estimator using several HSC-like image simulations and the main results are listed as follows. (i) For simulations which only contain isolated galaxies, the amplitude of multiplicative bias is below 1 % . (ii) For more realistic simulations which also contain blended galaxies, the blended galaxies are deblended by the first generation HSC deblender before shear measurement. Multiplicative bias of (-5.71 +- 0.31) % is found. The blending bias is calibrated by image simulations. Finally, we test the consistency and stability of this calibration., Comment: 18 pages, 3 figures, accepted by MNRAS

We study the possibility of using quadrupole moments of auto-convolved galaxy images to measure cosmic shear. The autoconvolution of an image corresponds to the inverse Fourier transformation of its power spectrum. The new method has the following advantages: the smearing effect due to the Point Spread Function (PSF) can be corrected by subtracting the quadrupole moments of the auto-convolved PSF; the centroid of the auto-convolved image is trivially identified; the systematic error due to noise can be directly removed in Fourier space; the PSF image can also contain noise, the effect of which can be similarly removed. With a large ensemble of simulated galaxy images, we show that the new method can reach a sub-percent level accuracy in general conditions, albeit with increasingly large stamp size for galaxies of less compact profiles., Comment: 8 pages, 6 figures, ApJ accepted