ABSTRACT
The next generation of weak lensing surveys will measure the matter distribution of the local universe with unprecedented precision, allowing the resolution of non-Gaussian features of the ...convergence field. This encourages the use of higher-order mass-map statistics for cosmological parameter inference. We extend the forward-modelling based methodology introduced in a previous forecast paper to match these new requirements. We provide multiple forecasts for the $w$CDM parameter constraints that can be expected from stage 3 and 4 weak lensing surveys. We consider different survey setups, summary statistics and mass map filters including wavelets. We take into account the shear bias, photometric redshift uncertainties, and intrinsic alignment. The impact of baryons is investigated and the necessary scale cuts are applied. We compare the angular power spectrum analysis to peak and minima counts as well as Minkowski functionals of the mass maps. We find a preference for Starlet over Gaussian filters. Our results suggest that using a survey setup with 10 instead of 5 tomographic redshift bins is beneficial. Adding cross-tomographic information improves the constraints on cosmology and especially on galaxy intrinsic alignment for all statistics. In terms of constraining power, we find the angular power spectrum and the peak counts to be equally matched for stage 4 surveys, followed by minima counts and the Minkowski functionals. Combining different summary statistics significantly improves the constraints and compensates the stringent scale cuts. We identify the most ‘cost-effective’ combination to be the angular power spectrum, peak counts and Minkowski functionals following Starlet filtering.
Shear peak statistics has gained a lot of attention recently as a practical alternative to the two-point statistics for constraining cosmological parameters. We perform a shear peak statistics ...analysis of the Dark Energy Survey (DES) Science Verification (SV) data, using weak gravitational lensing measurements from a 139 deg super( 2) field. We measure the abundance of peaks identified in aperture mass maps, as a function of their signal-to-noise ratio, in the signal-to-noise range ... To predict the peak counts as a function of cosmological parameters, we use a suite of N-body simulations spanning 158 models with varying ... and ..., fixing ..., to which we have applied the DES SV mask and redshift distribution. In our fiducial analysis we measure ..., after marginalizing over the shear multiplicative bias and the error on the mean redshift of the galaxy sample. We introduce models of intrinsic alignments, blending and source contamination by cluster members. These models indicate that peaks with ... would require significant corrections, which is why we do not include them in our analysis. We compare our results to the cosmological constraints from the two-point analysis on the SV field and find them to be in good agreement in both the central value and its uncertainty. We discuss prospects for future peak statistics analysis with upcoming DES data. (ProQuest: ... denotes formulae/symbols omitted.)
In this paper, we present results from the weak-lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the ...level of scrutiny employed in weak-lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches.
GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the point spread function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations, including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear, we present a general pseudo-C
ℓ formalism that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations.
The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by other shape measurement methods. The best methods achieve sub-per cent average biases. We find a strong dependence on accuracy as a function of signal-to-noise ratio, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can be brought. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.
GalSim is a collaborative, open-source project aimed at providing an image simulation tool of enduring benefit to the astronomical community. It provides a software library for generating images of ...astronomical objects such as stars and galaxies in a variety of ways, efficiently handling image transformations and operations such as convolution and rendering at high precision. We describe the GalSim software and its capabilities, including necessary theoretical background. We demonstrate that the performance of GalSim meets the stringent requirements of high precision image analysis applications such as weak gravitational lensing, for current datasets and for the Stage IV dark energy surveys of the Large Synoptic Survey Telescope, ESA’s Euclid mission, and NASA’s WFIRST-AFTA mission. The GalSim project repository is public and includes the full code history, all open and closed issues, installation instructions, documentation, and wiki pages (including a Frequently Asked Questions section). The GalSim repository can be found at https://github.com/GalSim-developers/GalSim.
Abstract
We use weak-lensing shear measurements to determine the mean mass of optically selected galaxy clusters in Dark Energy Survey Science Verification data. In a blinded analysis, we split the ...sample of more than 8000 redMaPPer clusters into 15 subsets, spanning ranges in the richness parameter 5 ≤ λ ≤ 180 and redshift 0.2 ≤ z ≤ 0.8, and fit the averaged mass density contrast profiles with a model that accounts for seven distinct sources of systematic uncertainty: shear measurement and photometric redshift errors; cluster-member contamination; miscentring; deviations from the NFW halo profile; halo triaxiality and line-of-sight projections. We combine the inferred cluster masses to estimate the joint scaling relation between mass, richness and redshift,
${\cal M}(\lambda ,z) \propto M_0 \lambda ^{F} (1+z)^{G}$
. We find
$M_0 \equiv \langle M_{200\mathrm{m}}\,|\,\lambda =30,z=0.5 \rangle = 2.35 \pm 0.22\ \rm {(stat)} \pm 0.12\ \rm {(sys)} \times \ 10^{14}\ \mathrm{M}_{{\odot }}$
, with
$F = 1.12\,\pm \,0.20\ \rm {(stat)}\, \pm \, 0.06\ \rm {(sys)}$
and
$G = 0.18\,\pm \, 0.75\ \rm {(stat)}\, \pm \, 0.24\ \rm {(sys)}$
. The amplitude of the mass–richness relation is in excellent agreement with the weak-lensing calibration of redMaPPer clusters in SDSS by Simet et al. and with the Saro et al. calibration based on abundance matching of SPT-detected clusters. Our results extend the redshift range over which the mass–richness relation of redMaPPer clusters has been calibrated with weak lensing from z ≤ 0.3 to z ≤ 0.8. Calibration uncertainties of shear measurements and photometric redshift estimates dominate our systematic error budget and require substantial improvements for forthcoming studies.
Abstract
Weak gravitational lensing has the potential to constrain cosmological parameters to high precision. However, as shown by the Shear Testing Programmes and Gravitational lensing Accuracy ...Testing challenges, measuring galaxy shears is a non-trivial task: various methods introduce different systematic biases which have to be accounted for. We investigate how pixel noise on the image affects the bias on shear estimates from a maximum likelihood forward model-fitting approach using a sum of co-elliptical Sérsic profiles, in complement to the theoretical approach of an associated paper. We evaluate the bias using a simple but realistic galaxy model and find that the effects of noise alone can cause biases of the order of 1-10 per cent on measured shears, which is significant for current and future lensing surveys. We evaluate a simulation-based calibration method to create a bias model as a function of galaxy properties and observing conditions. This model is then used to correct the simulated measurements. We demonstrate that, for the simple case in which the correct range of galaxy models is used in the fit, the calibration method can reduce noise bias to the level required for estimating cosmic shear in upcoming lensing surveys.
Weak lensing by large-scale structure is a powerful probe of cosmology and of the dark universe. This cosmic shear technique relies on the accurate measurement of the shapes and redshifts of ...background galaxies and requires precise control of systematic errors. Monte Carlo control loops (MCCL) is a forward modeling method designed to tackle this problem. It relies on the ultra fast image generator (UFig) to produce simulated images tuned to match the target data statistically, followed by calibrations and tolerance loops. We present the first end-to-end application of this method, on the Dark Energy Survey (DES) Year 1 wide field imaging data. We simultaneously measure the shear power spectrum Cℓ and the redshift distribution n(z) of the background galaxy sample. The method includes maps of the systematic sources, point spread function (PSF), an approximate Bayesian computation (ABC) inference of the simulation model parameters, a shear calibration scheme, and a fast method to estimate the covariance matrix. We find a close statistical agreement between the simulations and the DES Y1 data using an array of diagnostics. In a nontomographic setting, we derive a set of Cℓ and n(z) curves that encode the cosmic shear measurement, as well as the systematic uncertainty. Following a blinding scheme, we measure the combination of Ωm, σ8, and intrinsic alignment amplitude AIA, defined as S8DIA=σ8(Ωm/0.3)0.5DIA, where DIA=1−0.11(AIA−1). We find S8DIA=0.895−0.039+0.054, where systematics are at the level of roughly 60% of the statistical errors. We discuss these results in the context of earlier cosmic shear analyses of the DES Y1 data. Our findings indicate that this method and its fast runtime offer good prospects for cosmic shear measurements with future wide-field surveys.
We measure the weak lensing masses and galaxy distributions of four massive galaxy clusters observed during the Science Verification phase of the Dark Energy Survey (DES). This pathfinder study is ...meant to (1) validate the Dark Energy Camera (DECam) imager for the task of measuring weak lensing shapes, and (2) utilize DECam's large field of view to map out the clusters and their environments over 90 arcmin. We conduct a series of rigorous tests on astrometry, photometry, image quality, point spread function (PSF) modelling, and shear measurement accuracy to single out flaws in the data and also to identify the optimal data processing steps and parameters. We find Science Verification data from DECam to be suitable for the lensing analysis described in this paper. The PSF is generally well behaved, but the modelling is rendered difficult by a flux-dependent PSF width and ellipticity. We employ photometric redshifts to distinguish between foreground and background galaxies, and a red-sequence cluster finder to provide cluster richness estimates and cluster–galaxy distributions. By fitting Navarro–Frenk–White profiles to the clusters in this study, we determine weak lensing masses that are in agreement with previous work. For Abell 3261, we provide the first estimates of redshift, weak lensing mass, and richness. In addition, the cluster–galaxy distributions indicate the presence of filamentary structures attached to 1E 0657−56 and RXC J2248.7−4431, stretching out as far as 1°(approximately 20 Mpc), showcasing the potential of DECam and DES for detailed studies of degree-scale features on the sky.
We measure the redshift evolution of galaxy bias for a magnitude-limited galaxy sample by combining the galaxy density maps and weak lensing shear maps for a ∼116 deg2 area of the Dark Energy Survey ...(DES) Science Verification (SV) data. This method was first developed in Amara et al. and later re-examined in a companion paper with rigorous simulation tests and analytical treatment of tomographic measurements. In this work we apply this method to the DES SV data and measure the galaxy bias for a i < 22.5 galaxy sample. We find the galaxy bias and 1σ error bars in four photometric redshift bins to be 1.12 ± 0.19 (z = 0.2–0.4), 0.97 ± 0.15 (z = 0.4–0.6), 1.38 ± 0.39 (z = 0.6–0.8), and 1.45 ± 0.56 (z = 0.8–1.0). These measurements are consistent at the 2σ level with measurements on the same data set using galaxy clustering and cross-correlation of galaxies with cosmic microwave background lensing, with most of the redshift bins consistent within the 1σ error bars. In addition, our method provides the only σ8 independent constraint among the three. We forward model the main observational effects using mock galaxy catalogues by including shape noise, photo-z errors, and masking effects. We show that our bias measurement from the data is consistent with that expected from simulations. With the forthcoming full DES data set, we expect this method to provide additional constraints on the galaxy bias measurement from more traditional methods. Furthermore, in the process of our measurement, we build up a 3D mass map that allows further exploration of the dark matter distribution and its relation to galaxy evolution.
Convolutional Neural Networks (CNNs) are a cornerstone of the Deep Learning toolbox and have led to many breakthroughs in Artificial Intelligence. So far, these neural networks (NNs) have mostly been ...developed for regular Euclidean domains such as those supporting images, audio, or video. Because of their success, CNN-based methods are becoming increasingly popular in Cosmology. Cosmological data often comes as spherical maps, which make the use of the traditional CNNs more complicated. The commonly used pixelization scheme for spherical maps is the Hierarchical Equal Area isoLatitude Pixelisation (HEALPix). We present a spherical CNN for analysis of full and partial HEALPix maps, which we call DeepSphere. The spherical CNN is constructed by representing the sphere as a graph. Graphs are versatile data structures that can represent pairwise relationships between objects or act as a discrete representation of a continuous manifold. Using the graph-based representation, we define many of the standard CNN operations, such as convolution and pooling. With filters restricted to being radial, our convolutions are equivariant to rotation on the sphere, and DeepSphere can be made invariant or equivariant to rotation. This way, DeepSphere is a special case of a graph CNN, tailored to the HEALPix sampling of the sphere. This approach is computationally more efficient than using spherical harmonics to perform convolutions. We demonstrate the method on a classification problem of weak lensing mass maps from two cosmological models and compare its performance with that of three baseline classifiers, two based on the power spectrum and pixel density histogram, and a classical 2D CNN. Our experimental results show that the performance of DeepSphere is always superior or equal to the baselines. For high noise levels and for data covering only a smaller fraction of the sphere, DeepSphere achieves typically 10% better classification accuracy than the baselines.Finally, we show how learned filters can be visualized to introspect the NN. Code and examples are available at https://github.com/SwissDataScienceCenter/DeepSphere.