Systematic Biases in Weak Lensing Cosmology with the Dark Energy Survey

PhD thesis submitted to the University of Manchester, School of Physics and Astronomy, August 2017. Abstract: This thesis presents a practical guide to applying shear measurements as a cosmological tool. We first present one of two science-ready galaxy shape catalogues from Year 1 of the Dark Energy Survey (DES Y1), which covers 1500 square degrees in four bands griz, with a median redshift of 0.59. We describe the shape measurement process implemented by the DES Y1 im3shape catalogue, which contains 21.9M high-quality r-band bulge/disc fits. In Chapter 3 a new suite of image simulations, referred to as hoopoe, are presented. The hoopoe dataset is tailored to DES Y1 and includes realistic blending, spatial masks and variation in the point spread function. We derive shear corrections, which we show are robust to changes in calibration method, galaxy binning and variance within the simulated dataset. Sources of systematic uncertainty in the simulation-based shear calibration are discussed, leading to a final estimate of the 1 sigma uncertainties in the residual multiplicative bias after calibration of 0.025. Chapter 4 describes an extension of the analysis on the hoopoe simulations into a detailed investigation of the impact of galaxy neighbours on shape measurement and shear cosmology. Four mechanisms by which neighbours can have a non-negligible influence on shear measurement are identified. These effects, if ignored, would contribute a net multiplicative bias of m ~ 0.03 - 0.09 in DES Y1, though the precise impact will depend on both the measurement code and the selection cuts applied. We use the cosmological inference pipeline of DES Y1 to explore the cosmological implications of neighbour bias and show that omitting blending from the calibration simulation for DES Y1 would bias the inferred clustering amplitude S8 = sigma_8 (Omega_m /0.3)^0.5 by 1.5 sigma towards low values. Finally, we use the hoopoe simulations to test the effect of neighbour-induced spatial correlations in the multiplicative bias. We find the cosmological impact to be subdominant to statistical error at the current level of precision. Another major uncertainty in shear cosmology is the accuracy of our ensemble redshift distributions. Chapter 5 presents a numerical investigation into the combined constraining power of cosmic shear, galaxy clustering and their cross-correlation in DES Y1, and the potential for internal calibration of redshift errors. Introducing a moderate uniform bias into the redshift distributions used to model the weak lensing (WL) galaxies is shown to produce a > 2 sigma bias in S8. We demonstrate that this cosmological bias can be eliminated by marginalising over redshift error nuisance parameters. Strikingly, the cosmological constraint of the combined dataset is largely undiminished by the loss of prior information on the WL distributions. We demonstrate that this implicit self-calibration is the result of complementary degeneracy directions in the combined data. In Chapter 6 we present the preliminary results of an investigation into galaxy intrinsic alignments. Using the DES Y1 data, we show a clear dependence in alignment amplitude on galaxy type, in agreement with previous results. We subject these findings to a series of initial robustness tests. We conclude with a short overview of the work presented, and discuss prospects for the future

Modelling Galaxy Cluster Triaxiality in Stacked Cluster Weak Lensing Analyses

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

Systematic Bias in Cosmic Shear: Beyond the Fisher Matrix

We describe a method for computing the biases that systematic signals introduce in parameter estimation using a simple extension of the Fisher matrix formalism. This allows us to calculate the offset of the best fit parameters relative to the fiducial model, in addition to the usual statistical error ellipse. As an application, we study the impact that residual systematics in tomographic weak lensing measurements. In particular we explore three different types of shape measurement systematics: (i) additive systematic with no redshift evolution; (ii) additive systematic with redshift evolution; and (iii) multiplicative systematic. In each case, we consider a wide range of scale dependence and redshift evolution of the systematics signal. For a future DUNE-like full sky survey, we find that, for cases with mild redshift evolution, the variance of the additive systematic signal should be kept below 10^-7 to ensure biases on cosmological parameters that are sub-dominant to the statistical errors. For the multiplicative systematics, which depends on the lensing signal, we find the multiplicative calibration m0 needs to be controlled to an accuracy better than 10^-3. We find that the impact of systematics can be underestimated if their assumes redshift dependence is too simplistic. We provide simple scaling relations to extend these requirements to any survey geometry and discuss the impact of our results for current and future weak lensing surveys.Comment: Submitted to MNRAS. 11 pages, including 11 figures and 4 table