Weak Lensing Analysis Research Articles

The First Weak-lensing Analysis with the James Webb Space Telescope: SMACS J0723.3–7327

Utilizing the James Webb Space Telescope Early Release NIRCam Observations, we perform a weak-lensing analysis of the massive galaxy cluster SMACS J0723.3–7327 (z = 0.39). We investigate the spatial variation of the point-spread function (PSF) from the stars in the mosaic image. Our measurements show that the PSF for both modules has very small spatial and temporal variation with average complex ellipticity components of e 1 = 0.007 ± 0.001 and e 2 = 0.029 ± 0.001 in the observed north-up reference frame. We create PSF models through a principal component analysis of the stars and show that they properly account for the ellipticity of the PSF with residual shapes of e 1 = (0.3 ± 3.5) × 10−4 and e 2 = (1.8 ± 4.0) × 10−4. We select background galaxies by their photometric redshift and measure galaxy shapes by model fitting. Our weak-lensing source catalog achieves ∼215 galaxies arcmin−2. We map the convergence field of SMACSJ0723 and detect the cluster with a peak significance of 12.2σ. The mass distribution is found to elongate in the east–west direction with an extension to the northeast edge of the field of view where a candidate substructure is found in the Chandra X-ray imaging. We fit the tangential shear with a Navarro–Frenk–White model and estimate the mass of the cluster to be M 500 = 7.9 ± 1.1 × 1014 M ⊙ (M 200 = 11.4 ± 1.5 × 1014 M ⊙), which agrees with existing mass estimates. Combining the multiwavelength evidence from literature with our weak-lensing analysis, we hypothesize that SMACSJ0723 is observed near first pericenter passage and we identify candidate radio relics.

A simulation-based inference pipeline for cosmic shear with the Kilo-Degree Survey

ABSTRACT The standard approach to inference from cosmic large-scale structure data employs summary statistics that are compared to analytic models in a Gaussian likelihood with pre-computed covariance. To overcome the idealizing assumptions about the form of the likelihood and the complexity of the data inherent to the standard approach, we investigate simulation-based inference (SBI), which learns the likelihood as a probability density parameterized by a neural network. We construct suites of simulated summary statistics, exactly Gaussian distributed for validation purposes, for the most recent Kilo-Degree Survey (KiDS) weak gravitational lensing analysis and demonstrate that SBI recovers the full 12-dimensional KiDS posterior distribution with just under 104 simulations. We optimize the simulation strategy by initially covering the parameter space by a hypercube, followed by batches of actively learnt additional points. The data compression in our SBI implementation is robust to suboptimal choices of fiducial parameter values and of data covariance. Together with a fast simulator, SBI is therefore a competitive and more versatile alternative to standard inference.

Weak lensing tomographic redshift distribution inference for the Hyper Suprime-Cam Subaru Strategic Program three-year shape catalogue

ABSTRACTWe 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 cover a spatial subset of the shear catalogue. 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.

Beyond the ultradeep frontier fields and legacy observations (BUFFALO): a high-resolution strong+weak-lensing view of Abell 370

ABSTRACT The HST treasury programme BUFFALO provides extended wide-field imaging of the six Hubble Frontier Fields galaxy clusters. Here we present the combined strong and weak-lensing analysis of Abell 370, a massive cluster at z = 0.375. From the reconstructed total projected mass distribution in the 6 arcmin × 6 arcmin BUFFALO field-of-view, we obtain the distribution of massive substructures outside the cluster core and report the presence of a total of seven candidates, each with mass ∼5 × 1013 M⊙. Combining the total mass distribution derived from lensing with multiwavelength data, we evaluate the physical significance of each candidate substructure, and conclude that five out of the seven substructure candidates seem reliable, and that the mass distribution in Abell 370 is extended along the north-west and south-east directions. While this finding is in general agreement with previous studies, our detailed spatial reconstruction provides new insights into the complex mass distribution at large cluster-centric radius. We explore the impact of the extended mass reconstruction on the model of the cluster core and in particular, we attempt to physically explain the presence of an important external shear component, necessary to obtain a low root-mean-square separation between the model-predicted and observed positions of the multiple images in the cluster core. The substructures can only account for up to half the amplitude of the external shear, suggesting that more effort is needed to fully replace it by more physically motivated mass components. We provide public access to all the lensing data used as well as the different lens models.

The Chromatic Point-spread Function of Weak Lensing Measurement in the Chinese Space Station Survey Telescope

Weak gravitational lensing is a powerful tool in modern cosmology. To accurately measure the weak lensing signal, one has to control the systematic bias on a small level. One of the most difficult problems is how to correct the smearing effect of the Point-Spread Function (PSF) on the shape of the galaxies. The chromaticity of PSF for a broad-band observation can lead to new subtle effects. Since the PSF is wavelength-dependent and the spectrum energy distributions between stars and galaxies are different, the effective PSF measured from the star images will be different from those that smear the galaxies. Such a bias is called color bias. We estimate it in the optical bands of the Chinese Space Station Survey Telescope from simulated PSFs, and show the dependence on the color and redshift of the galaxies. Moreover, due to the spatial variation of spectra over the galaxy image, another higher-order bias exists: color gradient bias. Our results show that both color bias and color gradient bias are generally below 0.1% in CSST. Only for small-size galaxies, one needs to be careful about the color gradient bias in the weak lensing analysis using CSST data.

A general framework for removing point-spread function additive systematics in cosmological weak lensing analysis

ABSTRACT 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 modelling the impact of PSF systematics on the cosmic shear 2PCF and mitigating its impact on cosmological analysis. Our framework can detect PSF leakage and modelling 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 Ωm and S8 by ∼0.3σ. The traditional second moment-based model can only correct for a 0.1σ bias, leaving the contamination largely uncorrected. We conclude it is necessary to model both PSF second and fourth moment contaminations for HSC Y3 cosmic shear analysis. We also reanalyse the HSC Y1 cosmic shear analysis with our updated systematics model and identify a 0.07σ bias on Ω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.

The MillenniumTNG project: intrinsic alignments of galaxies and haloes

ABSTRACT The intrinsic alignment (IA) of observed galaxy shapes with the underlying cosmic web is a source of contamination in weak lensing surveys. Sensitive methods to identify the IA signal will therefore need to be included in the upcoming weak lensing analysis pipelines. Hydrodynamical cosmological simulations allow us to directly measure the intrinsic ellipticities of galaxies, and thus provide a powerful approach to predict and understand the IA signal. Here we employ the novel, large-volume hydrodynamical simulation MTNG740, a product of the MillenniumTNG (MTNG) project, to study the IA of galaxies. We measure the projected correlation functions between the intrinsic shape/shear of galaxies and various tracers of large-scale structure, w+g, w+m, w++ over the radial range $r_{\rm p} \in [0.02 , 200]\, h^{-1}{\rm Mpc}$ and at redshifts z = 0.0, 0.5, and 1.0. We detect significant signal-to-noise IA signals with the density field for both elliptical and spiral galaxies. We also find significant intrinsic shear–shear correlations for ellipticals. We further examine correlations of the intrinsic shape of galaxies with the local tidal field. Here we find a significant IA signal for elliptical galaxies assuming a linear model. We also detect a weak IA signal for spiral galaxies under a quadratic tidal torquing model. Lastly, we measure the alignment between central galaxies and their host dark-matter haloes, finding small to moderate misalignments between their principal axes that decline with halo mass.

UNCOVERing the extended strong lensing structures of Abell 2744 with the deepest JWST imaging

ABSTRACT We present a new parametric lens model for the massive galaxy cluster Abell 2744 based on new ultra-deep JWST imaging taken in the framework of the UNCOVER program. These observations constitute the deepest JWST images of a lensing cluster to date, adding to existing deep Hubble Space Telescope (HST) images and the recent JWST Early Release Science and Director’s Discretionary Time data taken for this field. The wide field of view of UNCOVER (∼45 arcmin2) extends beyond the cluster’s well-studied central core and reveals a spectacular wealth of prominent lensed features around two massive cluster sub-structures in the north and north-west, where no multiple images were previously known. We identify 75 new multiple images and candidates of 17 sources, 43 of which allow us, for the first time, to constrain the lensing properties and total mass distribution around these extended cluster structures using strong lensing (SL). Our model yields an effective Einstein radius of θE, main = 23.2 ± 2.3 arcsec for the main cluster core (for zs = 2), enclosing a mass of M(< θE, main) = (7.7 ± 1.1) × 1013 M⊙, and θE, NW = 13.1 ± 1.3 arcsec for the newly discovered north-western SL structure enclosing M(< θE, NW) = (2.2 ± 0.3) × 1013 M⊙. The northern clump is somewhat less massive with θE, N = 7.4 ± 0.7 arcsec enclosing M(< θE, N) = (0.8 ± 0.1) × 1013 M⊙. We find the northern sub-structures of Abell 2744 to broadly agree with the findings from weak lensing analyses and align with the filamentary structure found by these previous studies. Our model in particular reveals a large area of high magnification values between the various cluster structures, which will be paramount for lensed galaxy studies in the UNCOVER field. The model is made publicly available to accompany the first UNCOVER data release.

Modelling galaxy cluster triaxiality in stacked cluster weak lensing analyses

ABSTRACT Counts of galaxy clusters offer a high-precision probe of cosmology, but control of systematic errors will determine the accuracy of this measurement. Using Buzzard simulations, we quantify one such systematic, the triaxiality distribution of clusters identified with the redMaPPer optical cluster finding algorithm, which was used in the Dark Energy Survey Year-1 (DES Y1) cluster cosmology analysis. We test whether redMaPPer selection biases the clusters’ shape and orientation and find that it only biases orientation, preferentially selecting clusters with their major axes oriented along the line of sight. Modelling the richness–mass relation as log-linear, we find that the log-richness amplitude ln (A) is boosted from the lowest to highest orientation bin with a significance of 14σ, while the orientation dependence of the richness-mass slope and intrinsic scatter is minimal. We also find that the weak lensing shear-profile ratios of cluster-associated dark haloes in different orientation bins resemble a ‘bottleneck’ shape that can be quantified with a Cauchy function. We test the correlation of orientation with two other leading systematics in cluster cosmology – miscentering and projection – and find a null correlation. The resulting mass bias predicted from our templates confirms the DES Y1 finding that triaxiality is a leading source of bias in cluster cosmology. However, the richness-dependence of the bias confirms that triaxiality does not fully resolve the tension at low-richness between DES Y1 cluster cosmology and other probes. Our model can be used for quantifying the impact of triaxiality bias on cosmological constraints for upcoming weak lensing surveys of galaxy clusters.

Extremely expensive likelihoods: a variational-Bayes solution for precision cosmology

ABSTRACTWe present a variational-Bayes solution to compute non-Gaussian posteriors from extremely expensive likelihoods. Our approach is an alternative for parameter inference when Markov chain Monte Carlo (MCMC) sampling is numerically prohibitive or conceptually unfeasible. For example, when either the likelihood or the theoretical model cannot be evaluated at arbitrary parameter values, but only previously selected values, then traditional MCMC sampling is impossible, whereas our variational-Bayes solution still succeeds in estimating the full posterior. In cosmology, this occurs e.g. when the parametric model is based on costly simulations that were run for previously selected input parameters. We demonstrate the applicability of our posterior construction on the KiDS-450 weak lensing analysis, where we reconstruct the original KiDS MCMC posterior at 0.6 per cent of its former numerical posterior evaluations. The reduction in numerical cost implies that systematic effects which formerly exhausted the numerical budget could now be included.

Efficient computation of the super-sample covariance for stage IV galaxy surveys

Super-sample covariance (SSC) is an important effect for cosmological analyses that use the deep structure of the cosmic web; it may, however, be nontrivial to include it practically in a pipeline. We solve this difficulty by presenting a formula for the precision (inverse covariance) matrix and show applications to update likelihood or Fisher forecast pipelines. The formula has several advantages in terms of speed, reliability, stability, and ease of implementation. We present an analytical application to show the formal equivalence between three approaches to SSC: (i) at the usual covariance level, (ii) at the likelihood level, and (iii) with a quadratic estimator. We then present an application of this computationally efficient framework for studying the impact of inaccurate modelling of SSC responses for cosmological constraints from stage IV surveys. We find that a weak-lensing-only analysis is very sensitive to inaccurate modelling of the scale dependence of the response, which needs to be calibrated at the ∼15% level. The sensitivity to this scale dependence is less severe for the joint weak-lensing and galaxy clustering analysis (also known as 3×2pt). Nevertheless, we find that both the amplitude and scale-dependence of the responses have to be calibrated at better than 30%.

KiDS-Legacy calibration: Unifying shear and redshift calibration with the SKiLLS multi-band image simulations

We present SKiLLS, a suite of multi-band image simulations for the weak lensing analysis of the complete Kilo-Degree Survey (KiDS), dubbed KiDS-Legacy analysis. The resulting catalogues enable joint shear and redshift calibration, enhancing the realism and hence accuracy over previous efforts. To create a large volume of simulated galaxies with faithful properties and to a sufficient depth, we integrated cosmological simulations with high-quality imaging observations. We also improved the realism of simulated images by allowing the point spread function (PSF) to differ between CCD images, including stellar density variations and varying noise levels between pointings. Using realistic variable shear fields, we accounted for the impact of blended systems at different redshifts. Although the overall correction is minor, we found a clear redshift-bias correlation in the blending-only variable shear simulations, indicating the non-trivial impact of this higher-order blending effect. We also explored the impact of the PSF modelling errors and found a small yet noticeable effect on the shear bias. Finally, we conducted a series of sensitivity tests, including changing the input galaxy properties. We conclude that our fiducial shape measurement algorithm, lensfit, is robust within the requirements of lensing analyses with KiDS. As for future weak lensing surveys with tighter requirements, we suggest further investments in understanding the impact of blends at different redshifts, improving the PSF modelling algorithm and developing the shape measurement method to be less sensitive to the galaxy properties.

Fast and realistic large-scale structure from machine-learning-augmented random field simulations

ABSTRACT Producing thousands of simulations of the dark matter distribution in the Universe with increasing precision is a challenging but critical task to facilitate the exploitation of current and forthcoming cosmological surveys. Many inexpensive substitutes to full N-body simulations have been proposed, even though they often fail to reproduce the statistics of the smaller non-linear scales. Among these alternatives, a common approximation is represented by the lognormal distribution, which comes with its own limitations as well, while being extremely fast to compute even for high-resolution density fields. In this work, we train a generative deep learning model, mainly made of convolutional layers, to transform projected lognormal dark matter density fields to more realistic dark matter maps, as obtained from full N-body simulations. We detail the procedure that we follow to generate highly correlated pairs of lognormal and simulated maps, which we use as our training data, exploiting the information of the Fourier phases. We demonstrate the performance of our model comparing various statistical tests with different field resolutions, redshifts, and cosmological parameters, proving its robustness and explaining its current limitations. When evaluated on 100 test maps, the augmented lognormal random fields reproduce the power spectrum up to wavenumbers of $1 \, h \, \rm {Mpc}^{-1}$, and the bispectrum within 10 per cent, and always within the error bars, of the fiducial target simulations. Finally, we describe how we plan to integrate our proposed model with existing tools to yield more accurate spherical random fields for weak lensing analysis.

The eROSITA Final Equatorial-Depth Survey (eFEDS)

Context.We present the results of a systematic X-ray analysis of optically rich galaxy clusters detected by the Subaru Hyper Suprime-Cam (HSC) survey in the eROSITA Final Equatorial-Depth Survey (eFEDS) field.Aims.Through a joint analysis of the SRG (Spectrum Roentgen Gamma)/eROSITA and Subaru/HSC surveys, we aim to investigate the dynamical status of the optically selected clusters and to derive the cluster scaling relations.Methods.The sample consists of 43 optically selected galaxy clusters with a richness >40 in the redshift range of 0.16–0.89. We systematically analyzed the X-ray images and emission spectra using the eROSITA data. We identified the brightest cluster galaxy (BCG) using the optical and far-infrared databases. We evaluated the cluster’s dynamical status by measuring three quantities: offset between the X-ray peak and BCG position, the gas concentration parameter, and the number of galaxy-density peaks. We investigated the luminosity–temperature and mass–luminosity relations based on eROSITA X-ray spectra and HSC weak-lensing data analyses.Results.Based on these three measurements, we estimated the fraction of relaxed clusters to be 2(< 39)%, which is smaller than that of the X-ray-selected cluster samples. After correcting for a selection bias due to the richness cut, we obtained a shallowL−Tslope of 2.1 ± 0.5, which is consistent with the predictions of the self-similar model and the baseline model incorporating a mass–concentration relation. TheL−Mslope of 1.5 ± 0.3 is in agreement with the above-cited theoretical models as well as the data on the shear-selected clusters in the eFEDs field.Conclusions.Our analysis of high-richness optical clusters yields a small fraction of relaxed clusters and a shallow slope for the luminosity–temperature relation. This suggests that the average X-ray properties of the optical clusters are likely to be different from those observed in the X-ray samples. Thus, the joint eROSITA and HSC observations are a powerful tool in extending the analysis to a larger sample and understanding the selection effect with a view to establish cluster scaling relations.

Hubble Space Telescope and Hyper-Suprime-Cam Weak-lensing Study of the Equal-mass Dissociative Merger CIZA J0107.7+5408

A dissociative merger is formed by the interplay of ram pressure and gravitational forces, which can lead to a spatial displacement of the dark matter and baryonic components of the recently collided subclusters. CIZA J0107.7+5408 is a nearby (z = 0.105) dissociative merger that hosts two X-ray brightness peaks and a bimodal galaxy distribution. Analyzing MMT/Hectospec observations, we investigate the line-of-sight and spatial distribution of cluster galaxies. Utilizing deep, high-resolution Hubble Space Telescope Advanced Camera for Surveys imaging and large field-of-view Subaru Hyper-Suprime-Cam observations, we perform a weak-lensing analysis of CIZA J0107.7+5408. Our weak-lensing analysis detects a bimodal mass distribution that is spatially consistent with the cluster galaxies but significantly offset from the X-ray brightness peaks. Fitting two Navarro–Frenk–White halos to the lensing signal, we find an equal-mass merger with subcluster masses of M ⊙ and M ⊙. Moreover, the mass-to-light ratios of the subclusters, and , are found to be consistent with each other and within the range of mass-to-light ratios found for galaxy clusters.

Robust sampling for weak lensing and clustering analyses with the Dark Energy Survey

ABSTRACTRecent cosmological analyses rely on the ability to accurately sample from high-dimensional posterior distributions. A variety of algorithms have been applied in the field, but justification of the particular sampler choice and settings is often lacking. Here, we investigate three such samplers to motivate and validate the algorithm and settings used for the Dark Energy Survey (DES) analyses of the first 3 yr (Y3) of data from combined measurements of weak lensing and galaxy clustering. We employ the full DES Year 1 likelihood alongside a much faster approximate likelihood, which enables us to assess the outcomes from each sampler choice and demonstrate the robustness of our full results. We find that the ellipsoidal nested sampling algorithm multinest reports inconsistent estimates of the Bayesian evidence and somewhat narrower parameter credible intervals than the sliced nested sampling implemented in polychord. We compare the findings from multinest and polychord with parameter inference from the Metropolis–Hastings algorithm, finding good agreement. We determine that polychord provides a good balance of speed and robustness for posterior and evidence estimation, and recommend different settings for testing purposes and final chains for analyses with DES Y3 data. Our methodology can readily be reproduced to obtain suitable sampler settings for future surveys.

Impact of point spread function higher moments error on weak gravitational lensing – II. A comprehensive study

ABSTRACT Weak lensing is one of the most powerful probes for dark matter and dark energy science, although it faces increasing challenges in controlling systematic uncertainties as the statistical errors become smaller. The point spread function (PSF) needs to be precisely modeled to avoid systematic error on the weak lensing measurements. The weak lensing biases induced by errors in the PSF model second moments, i.e. its size and shape, are well-studied. However, Zhang et al. showed that errors in the higher moments of the PSF may also be a significant source of systematics for upcoming weak lensing surveys. Therefore, this work comprehensively investigate the modelling quality of PSF moments from the 3rd to 6th order, and propagate the PSFEx higher moments modelling error in the HSC survey data set to the weak lensing shear–shear correlation functions and their cosmological analyses. The overall multiplicative shear bias associated with errors in PSF higher moments can cause an ∼0.1σ shift on the cosmological parameters for LSST Y10, while the associated additive biases can induce 1σ uncertainties in cosmology parameter inference for LSST Y10, if not accounted. We compare the PSFEx model with PSF in Full FOV, and find similar performance in modelling the PSF higher moments. We conclude that PSF higher moment errors of the future PSF models should be reduced from those in current methods, otherwise needed to be explicitly modeled in the weak lensing analysis.

Accurate dark-matter halo elongation from weak-lensing stacking analysis

ABSTRACT Shape estimates that quantifying the halo anisotropic mass distribution are valuable parameters that provide information on their assembly process and evolution. Measurements of the mean shape estimates for a sample of cluster-sized haloes can be used to test halo formation scenarios, as well as improving the modelling of potential biases in constraining cosmological parameters using these systems. In this work, we test the recovery of halo cluster shapes and masses applying weak-lensing stacking techniques. To this end, we use lensing shear and a new dark-matter halo catalogue, derived from the light-cone output of the cosmological simulation MICE-GC. We perform this study by combining the lensing signals obtained for several samples of haloes, selected according to their mass and redshift, taking into account the main directions of the dark-matter distributions. In the analysis, we test the impact of several potential introduced systematics, such as the adopted modelling, the contribution of the neighbouring mass distribution, miscentring, and misalignment effects. Our results show that when some considerations regarding the halo relaxation state are taken into account, the lensing semi-axial ratio estimates are in agreement within a 5 per cent with the mean shapes of the projected dark-matter particle distribution of the stacked haloes. The presented methodology provides a useful tool to derive reliable shapes of galaxy clusters and to contrast them with those expected from numerical simulations. Furthermore, our proposed modelling, that takes into account the contribution of neighbouring haloes, allows to constraint the elongation of the surrounding mass distribution.

Cosmological inference from an emulator based halo model. II. Joint analysis of galaxy-galaxy weak lensing and galaxy clustering from HSC-Y1 and SDSS

We present high-fidelity cosmology results from a blinded joint analysis of galaxy-galaxy weak lensing (ΔΣ) and projected galaxy clustering (wp) 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<z<0.7. We define luminosity-limited samples of SDSS galaxies to serve as the tracers of wp in three spectroscopic redshift bins, and as the lens samples for ΔΣ. For the ΔΣ measurements, we select a single sample of 4×106 source galaxies over 140 deg2 from HSC-Y1 with photometric redshifts (photo z) greater than 0.75, enabling a better handle of photo-z errors by comparing the ΔΣ amplitudes for the three lens redshift bins. The deep, high-quality HSC-Y1 data enable significant detections of the ΔΣ signals, with integrated signal-to-noise ratio S/N∼15 in the range 3≤R/[h−1 Mpc]≤30 for the three lens samples, despite the small area coverage. For cosmological parameter inference, we use an input galaxy-halo connection model built on the dark emulator package (which uses an ensemble set of high-resolution N-body simulations and enables fast, accurate computation of the clustering observables) with a halo occupation distribution that includes nuisance parameters to marginalize over modeling uncertainties. We model the ΔΣ and wp measurements on scales from R≃3 and 2 h−1 Mpc, respectively, up to 30 h−1 Mpc (therefore excluding the baryon acoustic oscillations information) assuming a flat ΛCDM cosmology, marginalizing over about 20 nuisance parameters and demonstrating the robustness of our results to them. With various tests using mock catalogs described in Miyatake et al. [preceding paper, Phys. Rev. D 106, 083519 (2022)], we show that any bias in the clustering amplitude S8≡σ8(Ωm/0.3)0.5 due to uncertainties in the galaxy-halo connection is less than ∼50% of the statistical uncertainty of S8, unless the assembly biaseffect is unexpectedly large. Our best-fit models have S8=0.795−0.042+0.049 (mode and 68% credible interval) for the flat ΛCDM model; we find tighter constraints on the quantity S8(α=0.17)≡σ8(Ωm/0.3)0.17=0.745−0.031+0.039.23 MoreReceived 29 November 2021Accepted 22 August 2022DOI:https://doi.org/10.1103/PhysRevD.106.083520© 2022 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasCosmologyDark energyGravitational lensesGravitation, Cosmology & Astrophysics

ShapePipe: A new shape measurement pipeline and weak-lensing application to UNIONS/CFIS data

Context. The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) is an ongoing collaboration that will provide the largest deep photometric survey of the northern sky in four optical bands to date. As part of this collaboration, the Canada-France Imaging Survey (CFIS) is observing r-band data with an average seeing of 0.65 arcsec, which is complete to magnitude 24.5 and thus ideal for weak-lensing studies. Aims. We perform the first weak-lensing analysis of CFIS r-band data over an area spanning 1700 deg2 of the sky. We create a catalogue with measured shapes for 40 million galaxies, corresponding to an effective density of 6.8 galaxies per square arcminute, and demonstrate a low level of systematic biases. This work serves as the basis for further cosmological studies that will use the full UNIONS survey of 4800 deg2 when completed. Methods. Here we present SHAPEPIPE, a newly developed weak-lensing pipeline. This pipeline makes use of state-of-the-art methods such as NGMIX for accurate galaxy shape measurement. Shear calibration is performed with metacalibration. We carry out extensive validation tests on the point spread function (PSF) and on the galaxy shapes. In addition, we create realistic image simulations to validate the estimated shear. Results. We quantify the PSF model accuracy and show that the level of systematics is low as measured by the PSF residuals. Their effect on the shear two-point correlation function is sub-dominant compared to the cosmological contribution on angular scales &lt; 100′. The additive shear bias is below 5 × 10−4, and the residual multiplicative shear bias is at most 10−3 as measured on image simulations. Using complete orthogonal sets of E-/B-mode integrals (COSEBIs), we show that there are no significant B-modes present in second-order shear statistics. We present convergence maps and see clear correlations of the E-mode with known cluster positions. We measure the stacked tangential shear profile around Planck clusters at a significance higher than 4σ.