A Way Forward for Cosmic Shear: Monte-Carlo Control Loops

Monthly Notices of The Royal Astronomical Society, 2010

We present the results of the GREAT08 Challenge, a blind analysis challenge to infer weak gravitational lensing shear distortions from images. The primary goal was to stimulate new ideas by presenting the problem to researchers outside the shear measurement community. Six GREAT08 Team methods were presented at the launch of the Challenge and five additional groups submitted results during the 6 month competition. Participants analyzed 30 million simulated galaxies with a range in signal to noise ratio, point-spread function ellipticity, galaxy size, and galaxy type. The large quantity of simulations allowed shear measurement methods to be assessed at a level of accuracy suitable for currently planned future cosmic shear observations for the first time. Different methods perform well in different parts of simulation parameter space and come close to the target level of accuracy in several of these. A number of fresh ideas have emerged as a result of the Challenge including a re-examination of the process of combining information from different galaxies, which reduces the dependence on realistic galaxy modelling. The image simulations will become increasingly sophisticated in future GREAT challenges, meanwhile the GREAT08 simulations remain as a benchmark for additional developments in shear measurement algorithms.

The Astrophysical Journal, 2017

We infer gravitational lensing shear and convergence fields from galaxy ellipticity catalogs under a spatial process prior for the lensing potential. We demonstrate the performance of our algorithm with simulated Gaussian-distributed cosmological lensing shear maps and a reconstruction of the mass distribution of the merging galaxy cluster Abell 781 using galaxy ellipticities measured with the Deep Lens Survey. Given interim posterior samples of lensing shear or convergence fields on the sky, we describe an algorithm to infer cosmological parameters via lens field marginalization. In the most general formulation of our algorithm we make no assumptions about weak shear or Gaussian distributed shape noise or shears. Because we require solutions and matrix determinants of a linear system of dimension that scales with the number of galaxies, we expect our algorithm to require parallel high-performance computing resources for application to ongoing wide field lensing surveys.

Astronomy & Astrophysics, 2021

Upcoming surveys will map the growth of large-scale structure with unprecented precision, improving our understanding of the dark sector of the Universe. Unfortunately, much of the cosmological information is encoded on small scales, where the clustering of dark matter and the effects of astrophysical feedback processes are not fully understood. This can bias the estimates of cosmological parameters, which we study here for a joint analysis of mock Euclid cosmic shear and Planck cosmic microwave background data. We use different implementations for the modelling of the signal on small scales and find that they result in significantly different predictions. Moreover, the different non-linear corrections lead to biased parameter estimates, especially when the analysis is extended into the highly non-linear regime, with the Hubble constant, H0, and the clustering amplitude, σ8, affected the most. Improvements in the modelling of non-linear scales will therefore be needed if we are to r...

The Astrophysical Journal, 2013

We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 sq. degree fields with two National Optical Astronomy Observatory (NOAO) 4-meter telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the pointspread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmologydependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the Ω M − σ 8 likelihood contour tightly, providing a joint constraint of Ω M = 0.262 ± 0.051 and σ 8 = 0.868 ± 0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the Ω M − σ 8 degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7-year (WMAP7) likelihood data, we obtain Ω M = 0.278 ± 0.018 and σ 8 = 0.815 ± 0.020.

Monthly Notices of the Royal Astronomical Society, 2012

Weak gravitational lensing has the potential to constrain cosmological parameters to high precision. However, as shown by the Shear TEsting Programmes (STEP) and GRavitational lEnsing Accuracy Testing (GREAT) Challenges, measuring galaxy shears is a nontrivial 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 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 order 1-10% 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 this method can effectively reduce noise bias so that shear measurement reaches the level of accuracy required for estimating cosmic shear in upcoming lensing surveys.

The Astrophysical Journal, 2003

We present a measurement of cosmic shear on scales ranging from 10 ′′ to 2 ′ in 347 WFPC2 images of random fields. Our result is based on shapes measured via image fitting and on a simple statistical technique; careful calibration of each step allows us to quantify our systematic uncertainties and to measure the cosmic shear down to very small angular scales. The WFPC2 images provide a robust measurement of the cosmic shear signal decreasing from 5.2% at 10 ′′ to 2.2% at 130 ′′ .

Astronomy & Astrophysics, 2012

Context. Forthcoming experiments will enable us to determine tomographic shear spectra at a high precision level. Most predictions about them have until now been biased on algorithms yielding the expected linear and non-linear spectrum of density fluctuations. Even when simulations have been used, so-called Halofit (Smith et al 2003) predictions on fairly large scales have been needed. Aims. We wish to go beyond this limitation. Methods. We perform N-body and hydrodynamical simulations within a sufficiently large cosmological volume to allow a direct connection between simulations and linear spectra. While covering large length-scales, the simulation resolution is good enough to allow us to explore the high-ℓ harmonics of the cosmic shear (up to ℓ ∼ 50000), well into the domain where baryon physics becomes important. We then compare shear spectra in the absence and in presence of various kinds of baryon physics, such as radiative cooling, star formation, and supernova feedback in the form of galactic winds. Results. We distinguish several typical properties of matter fluctuation spectra in the different simulations and test their impact on shear spectra. Conclusions. We compare our outputs with those obtainable using approximate expressions for non-linear spectra, and identify substantial discrepancies even between our results and those of purely N-body results. Our simulations and the treatment of their outputs however enable us, for the first time, to obtain shear results taht are fully independent of any approximate expression, also in the high-ℓ range, where we need to incorporate a non-linear power spectrum of density perturbations, and the effects of baryon physics. This will allow us to fully exploit the cosmological information contained in future high-sensitivity cosmic shear surveys, exploring the physics of cosmic shears via weak lensing measurements.

2016

Physical Review D

Monthly Notices of the Royal Astronomical Society, 2017

Weak lensing has become an increasingly important tool in cosmology and the use of galaxy shapes to measure cosmic shear has become routine. The weak-lensing distortion tensor contains two other effects in addition to the two components of shear: the convergence and rotation. The rotation mode is not measurable using the standard cosmic shear estimators based on galaxy shapes, as there is no information on the original shapes of the images before they were lensed. Due to this, no estimator has been proposed for the rotation mode in cosmological weak-lensing surveys, and the rotation mode has never been constrained. Here, we derive an estimator for this quantity, which is based on the use of radio polarization measurements of the intrinsic position angles of galaxies. The rotation mode can be sourced by physics beyond cold dark matter (CDM), and also offers the chance to perform consistency checks of CDM and of weak-lensing surveys themselves. We present simulations of this estimator and show that, for the pedagogical example of cosmic string spectra, this estimator could detect a signal that is consistent with the constraints from Planck. We examine the connection between the rotation mode and the shear B modes and thus how this estimator could help control systematics in future radio weak-lensing surveys.