No Arabic abstract
(Abridged) Weak gravitational lensing induces distortions on the images of background galaxies, and thus provides a direct measure of mass fluctuations in the universe. Since the distortions induced by lensing on the images of background galaxies are only of a few percent, a reliable measurement demands very accurate galaxy shape estimation and a careful treatment of systematic effects. Here, we present a study of a shear measurement method using detailed simulations of artificial images. The images are produced using realisations of a galaxy ensemble drawn from the HST Groth strip. We consider realistic observational effects including atmospheric seeing, PSF anisotropy and pixelisation, incorporated in a manner to reproduce actual observations with the William Herschel Telescope. By applying an artificial shear to the simulated images, we test the shear measurement method proposed by Kaiser, Squires & Broadhurst (1995, KSB). Overall, we find the KSB method to be reliable with several provisos. To guide future weak lensing surveys, we study how seeing size, exposure time and pixelisation affect the sensitivity to shear. In addition, we study the impact of overlapping isophotes of neighboring galaxies, and find that this effect can produce spurious lensing signals on small scales. We discuss the prospects of using the KSB method for future, more sensitive, surveys. Numerical simulations of this kind are a required component of present and future analyses of weak lensing surveys.
We present a study of weak lensing shear measurements for simulated galaxy images at radio wavelengths. We construct a simulation pipeline into which we can input galaxy images of known ellipticity, and with which we then simulate observations with eMERLIN and the international LOFAR array. The simulations include the effects of the CLEAN algorithm, uv sampling, observing angle, and visibility noise, and produce realistic restored images of the galaxies. We apply a shapelet-based shear measurement method to these images and test our ability to recover the true source ellipticities. We model and deconvolve the effective PSF, and find suitable parameters for CLEAN and shapelet decomposition of galaxies. We demonstrate that ellipticities can be measured faithfully in these radio simulations, with no evidence of an additive bias and a modest (10%) multiplicative bias on the ellipticity measurements. Our simulation pipeline can be used to test shear measurement procedures and systematics for the next generation of radio telescopes.
We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $lesssim 2%$ out to scales $ell approx 4000$. For lensing peak counts, the agreement is better than $5%$ for peaks with signal-to-noise $lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.
Weak gravitational lensing measurements are traditionally made at optical wavelengths where many highly resolved galaxy images are readily available. However, the Square Kilometre Array (SKA) holds great promise for this type of measurement at radio wavelengths owing to its greatly increased sensitivity and resolution over typical radio surveys. The key to successful weak lensing experiments is in measuring the shapes of detected sources to high accuracy. In this document we describe a simulation pipeline designed to simulate radio images of the quality required for weak lensing, and will be typical of SKA observations. We provide as input, images with realistic galaxy shapes which are then simulated to produce images as they would have been observed with a given radio interferometer. We exploit this pipeline to investigate various stages of a weak lensing experiment in order to better understand the effects that may impact shape measurement. We first show how the proposed SKA1-Mid array configurations perform when we compare the (known) input and output ellipticities. We then investigate how making small changes to these array configurations impact on this input-outut ellipticity comparison. We also demonstrate how alternative configurations for SKA1-Mid that are smaller in extent, and with a faster survey speeds produce similar performance to those originally proposed. We then show how a notional SKA configuration performs in the same shape measurement challenge. Finally, we describe ongoing efforts to utilise our simulation pipeline to address questions relating to how applicable current (mostly originating from optical data) shape measurement techniques are to future radio surveys. As an alternative to such image plane techniques, we lastly discuss a shape measurement technique based on the shapelets formalism that reconstructs the source shapes directly from the visibility data.
Given a foreground galaxy-density field or shear field, its cross-correlation with the shear field from a background population of source galaxies scales with the source redshift in a way that is specific to lensing. Such a source-scaling can be exploited to effectively measure geometrical distances as a function of redshift and thereby constrain dark energy properties, free of any assumptions about the galaxy-mass/mass power spectrum (its shape, amplitude or growth). Such a geometrical method can yield a ~ 0.03 - 0.07 f_{sky}^{-1/2} measurement on the dark energy abundance and equation of state, for a photometric redshift accuracy of dz ~ 0.01 - 0.05 and a survey with median redshift of ~ 1. While these constraints are weaker than conventional weak lensing methods, they provide an important consistency check because the geometrical method carries less theoretical baggage: there is no need to assume any structure formation model (e.g. CDM). The geometrical method is at the most conservative end of a whole spectrum of methods which obtain smaller errorbars by making more restrictive assumptions -- we discuss some examples. Our geometrical approach differs from previous investigations along similar lines in three respects. First, the source-scaling we propose to use is less demanding on the photometric redshift accuracy. Second, the scaling works for both galaxy-shear and shear-shear correlations. Third, we find that previous studies underestimate the statistical errors associated with similar geometrical methods, the origin of which is discussed.
Full ray-tracing maps of gravitational lensing, constructed from N-Body simulations, represent a fundamental tool to interpret present and future weak lensing data. However the limitation of computational resources and storage capabilities severely restrict the number of realizations that can be performed in order to accurately sample both the cosmic shear models and covariance matrices. In this paper we present a halo model formalism for weak gravitational lensing that alleviates these issues by producing weak-lensing mocks at a reduced computational cost. Our model takes as input the halo population within a desired light-cone and the linear power spectrum of the underlined cosmological model. We examine the contribution given by the presence of substructures within haloes to the cosmic shear power spectrum and quantify it to the percent level. Our method allows us to reconstruct high-resolution convergence maps, for any desired source redshifts, of light-cones that realistically trace the matter density distribution in the universe, account for masked area and sample selections. We compare our analysis on the same large scale structures constructed using ray-tracing techniques and find very good agreements both in the linear and non-linear regimes up to few percent levels. The accuracy and speed of our method demonstrate the potential of our halo model for weak lensing statistics and the possibility to generate a large sample of convergence maps for different cosmological models as needed for the analysis of large galaxy redshift surveys.