No Arabic abstract
Mass measurements of astronomical objects are most wanted but still elusive. We need them to trace the formation and evolution of cosmic structure but we can get direct measurements only for a minority. This lack can be circumvented with a proxy and a scaling relation. The twofold goal of estimating the unbiased relation and finding the right proxy value to plug in can be hampered by systematics, selection effects, Eddington/Malmquist biases and time evolution. We present a Bayesian hierarchical method which deals with these issues. Masses to be predicted are treated as missing data in the regression and are estimated together with the scaling parameters. The calibration subsample with measured masses does not need to be representative of the full sample as far as it follows the same scaling relation. We apply the method to forecast weak lensing calibrated masses of the Planck, redMaPPer and MCXC clusters. Planck masses are biased low with respect to weak lensing calibrated masses, with a bias more pronounced for high redshift clusters. MCXC masses are under-estimated by ~20 per cent, which may be ascribed to hydrostatic bias. Packages and catalogs are made available with the paper.
We simultaneously present constraints on the stellar-to-halo mass relation for central and satellite galaxies through a weak lensing analysis of spectroscopically classified galaxies. Using overlapping data from the fourth data release of the Kilo-Degree Survey (KiDS), and the Galaxy And Mass Assembly survey (GAMA), we find that satellite galaxies are hosted by halo masses that are $0.53 pm 0.39$ dex (68% confidence, $3sigma$ detection) smaller than those of central galaxies of the same stellar mass (for a stellar mass of $log(M_{star}/M_{odot}) = 10.6$). This is consistent with galaxy formation models, whereby infalling satellite galaxies are preferentially stripped of their dark matter. We find consistent results with similar uncertainties when comparing constraints from a standard azimuthally averaged galaxy-galaxy lensing analysis and a two-dimensional likelihood analysis of the full shear field. As the latter approach is somewhat biased due to the lens incompleteness and as it does not provide any improvement to the precision when applied to actual data, we conclude that stacked tangential shear measurements are best-suited for studies of the galaxy-halo connection.
Mapping the underlying density field, including non-visible dark matter, using weak gravitational lensing measurements is now a standard tool in cosmology. Due to its importance to the science results of current and upcoming surveys, the quality of the convergence reconstruction methods should be well understood. We compare three methods: Kaiser-Squires (KS), Wiener filter, and GLIMPSE. KS is a direct inversion, not accounting for survey masks or noise. The Wiener filter is well-motivated for Gaussian density fields in a Bayesian framework. GLIMPSE uses sparsity, aiming to reconstruct non-linearities in the density field. We compare these methods with several tests using public Dark Energy Survey (DES) Science Verification (SV) data and realistic DES simulations. The Wiener filter and GLIMPSE offer substantial improvements over smoothed KS with a range of metrics. Both the Wiener filter and GLIMPSE convergence reconstructions show a 12 per cent improvement in Pearson correlation with the underlying truth from simulations. To compare the mapping methods abilities to find mass peaks, we measure the difference between peak counts from simulated {Lambda}CDM shear catalogues and catalogues with no mass fluctuations (a standard data vector when inferring cosmology from peak statistics); the maximum signal-to-noise of these peak statistics is increased by a factor of 3.5 for the Wiener filter and 9 for GLIMPSE. With simulations we measure the reconstruction of the harmonic phases; the phase residuals concentration is improved 17 per cent by GLIMPSE and 18 per cent by the Wiener filter. The correlation between reconstructions from data and foreground redMaPPer clusters is increased 18 per cent by the Wiener filter and 32 per cent by GLIMPSE.
The weak distortions produced by gravitational lensing in the images of background galaxies provide a method to measure directly the distribution of mass in the universe. However this technique requires high precision measurements of the lensing shear and cautious corrections for systematic effects. Kaiser, Squires, & Broadhurst (1995) proposed a method to calibrate the ellipticity-shear relation in the presence of Point Spread Function (PSF) anisotropies and camera distortions. We revisit the KSB method and show that both the PSF and the camera distortions can be corrected for using source moments, as opposed to ellipticities. We clarify the applicability of some of the approximations made in this method. We derive expressions for the corrections which only involve the galaxy moments. We derive an explicit relation between the shear and the average ellipticity. We discuss the shortcomings of the method, and test its validity using numerical simulations. As an application of the method, we repeat the analysis of the HST WFPC2 camera performed by Hoekstra et al. (1998). We confirm the presence of sizable (10%) PSF ellipticities at the edge of the WFPC2 chips. We also show that the PSF ellipticity varies by as much as 2% over time. We use these measurements to correct the shape of galaxies in the HST Survey Strip (``Groth Strip). By considering the dependence of the ellipticities on object size, we show that, after corrections, the residual systematic uncertainty for galaxies with radii greater than 0.15 arcsec, is about 0.4%, when averaged over each chip. We discuss how these results provide good prospects for measuring weak lensing by large-scale structure with deep HST surveys.
Unbiased and precise mass calibration of galaxy clusters is crucial to fully exploit galaxy clusters as cosmological probes. Stacking of weak lensing signal allows us to measure observable-mass relations down to less massive halos halos without extrapolation. We propose a Bayesian inference method to constrain the intrinsic scatter of the mass proxy in stacked analyses. The scatter of the stacked data is rescaled with respect to the individual scatter based on the number of binned clusters. We apply this method to the galaxy clusters detected with the AMICO (Adaptive Matched Identifier of Clustered Objects) algorithm in the third data release of the Kilo-Degree Survey. The results confirm the optical richness as a low scatter mass proxy. Based on the optical richness and the calibrated weak lensing mass-richness relation, mass of individual objects down to ~10^13 solar masses can be estimated with a precision of ~20 per cent.
We explore the effect of massive neutrinos on the weak lensing shear bispectrum using the Cosmological Massive Neutrino Simulations. We find that the primary effect of massive neutrinos is to suppress the amplitude of the bispectrum with limited effect on the bispectrum shape. The suppression of the bispectrum amplitude is a factor of two greater than the suppression of the small scale power-spectrum. For an LSST-like weak lensing survey that observes half of the sky with five tomographic redshift bins, we explore the constraining power of the bispectrum on three cosmological parameters: the sum of the neutrino mass $sum m_ u$, the matter density $Omega_m$ and the amplitude of primordial fluctuations $A_s$. Bispectrum measurements alone provide similar constraints to the power spectrum measurements and combining the two probes leads to significant improvements than using the latter alone. We find that the joint constraints tighten the power spectrum $95%$ constraints by $sim 32%$ for $sum m_ u$, $13%$ for $Omega_m$ and $57%$ for $A_s$ .