A hydrodynamical halo model for weak-lensing cross correlations

On the scale of galactic haloes, the distribution of matter in the cosmos is affected by energetic, non-gravitational processes; so-called baryonic feedback. A lack of knowledge about the details of how feedback processes redistribute matter is a source of uncertainty for weak-lensing surveys, which accurately probe the clustering of matter in the Universe over a wide range of scales. We develop a cosmology-dependent model for the matter distribution that simultaneously accounts for the clustering of dark matter, gas and stars. We inform our model by comparing it to power spectra measured from the BAHAMAS suite of hydrodynamical simulations. As well as considering matter power spectra, we also consider spectra involving the electron-pressure field, which directly relates to the thermal Sunyaev-Zel'dovich (tSZ) effect. We fit parameters in our model so that it can simultaneously model both matter and pressure data and such that the distribution of gas as inferred from tSZ has influence on the matter spectrum predicted by our model. We present two variants; one that matches the feedback-induced suppression seen in the matter-matter power spectrum at the per-cent level and a second that matches the matter-matter data slightly less well (~2 per cent), but that is able to simultaneously model the matter-electron pressure spectrum at the ~15 per-cent level. We envisage our models being used to simultaneously learn about cosmological parameters and the strength of baryonic feedback using a combination of tSZ and lensing auto- and cross-correlation data

Galaxy cluster mass estimation with deep learning and hydrodynamical simulations

We evaluate the ability of Convolutional Neural Networks (CNNs) to predict galaxy cluster masses in the BAHAMAS hydrodynamical simulations. We train four separate single-channel networks using: stellar mass, soft X-ray flux, bolometric X-ray flux, and the Compton $y$ parameter as observational tracers, respectively. Our training set consists of $\sim$4800 synthetic cluster images generated from the simulation, while an additional $\sim$3200 images form a validation set and a test set, each with 1600 images. In order to mimic real observation, these images also contain uncorrelated structures located within 50 Mpc in front and behind clusters and seen in projection, as well as instrumental systematics including noise and smoothing. In addition to CNNs for all the four observables, we also train a `multi-channel' CNN by combining the four observational tracers. The learning curves of all the five CNNs converge within 1000 epochs. The resulting predictions are especially precise for halo masses in the range $10^{13.25}M_{\odot}<M<10^{14.5}M_{\odot}$, where all five networks produce mean mass biases of order $\approx$1\% with a scatter of $\lesssim$20\%. The network trained with Compton $y$ parameter maps yields the most precise predictions. We interpret the network's behaviour using two diagnostic tests to determine which features are used to predict cluster mass. The CNN trained with stellar mass images detect galaxies (not surprisingly), while CNNs trained with gas-based tracers utilise the shape of the signal to estimate cluster mass

The BAHAMAS project: the CMB--large-scale structure tension and the roles of massive neutrinos and galaxy formation

Recent studies have presented evidence for tension between the constraints on Omega_m and sigma_8 from the cosmic microwave background (CMB) and measurements of large-scale structure (LSS). This tension can potentially be resolved by appealing to extensions of the standard model of cosmology and/or untreated systematic errors in the modelling of LSS, of which baryonic physics has been frequently suggested. We revisit this tension using, for the first time, carefully-calibrated cosmological hydrodynamical simulations, which thus capture the back reaction of the baryons on the total matter distribution. We have extended the BAHAMAS simulations to include a treatment of massive neutrinos, which currently represents the best motivated extension to the standard model. We make synthetic thermal Sunyaev-Zel'dovich effect, weak galaxy lensing, and CMB lensing maps and compare to observed auto- and cross-power spectra from a wide range of recent observational surveys. We conclude that: i) in general there is tension between the primary CMB and LSS when adopting the standard model with minimal neutrino mass; ii) after calibrating feedback processes to match the gas fractions of clusters, the remaining uncertainties in the baryonic physics modelling are insufficient to reconcile this tension; and iii) invoking a non-minimal neutrino mass, typically of 0.2-0.4 eV (depending on the priors on the other relevant cosmological parameters and the datasets being modelled), can resolve the tension. This solution is fully consistent with separate constraints on the summed neutrino mass from the primary CMB and baryon acoustic oscillations, given the internal tensions in the Planck primary CMB dataset

Probing hot gas around luminous red galaxies through the Sunyaev-Zel'dovich effect

We construct the mean thermal Sunyaev-Zel'dovich (tSZ) Comptonization y profile around Luminous Red Galaxies (LRGs) in the redshift range 0.16 < z < 0.47 from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) using the Planck y map. The mean central tSZ signal for the full sample is y ~ 1.8 * 10^(-7) and we detect tSZ emission out to ~30 arcmin, which is well beyond the 10 arcmin angular resolution of the y map and well beyond the virial radii of the LRGs. We compare the measured profile with predictions from the cosmo-OWLS suite of cosmological hydrodynamical simulations. This comparison agrees well for models that include feedback from active galactic nuclei (AGN), but not with hydrodynamic models without this energetic feedback mechanism. This suggests that an additional heating mechanism is required over SNe feedback and star formation to explain the y data profile. We also compare our results with predictions based on the halo model with a universal pressure profile (UPP) giving the y signal. The predicted profile is consistent with the data, but only if we account for the clustering of haloes via a two-halo term and if halo masses are estimated using the mean stellar-to-halo mass (SHM) relation of Coupon et al. (2015) or Wang et al.(2016) estimated from gravitational lensing measurements. We also discuss the importance of scatter in the SHM relation on the model predictions

The Shear TEsting Programme 1: Weak lensing analysis of simulated ground-based observations

The Shear TEsting Programme, STEP, is a collaborative project to improve the accuracy and reliability of all weak lensing measurements in preparation for the next generation of wide-field surveys. In this first STEP paper we present the results of a blind analysis of simulated ground-based observations of relatively simple galaxy morphologies. The most successful methods are shown to achieve percent level accuracy. From the cosmic shear pipelines that have been used to constrain cosmology, we find weak lensing shear measured to an accuracy that is within the statistical errors of current weak lensing analyses, with shear measurements accurate to better than 7%. The dominant source of measurement error is shown to arise from calibration uncertainties where the measured shear is over or under-estimated by a constant multiplicative factor. This is of concern as calibration errors cannot be detected through standard diagnostic tests. The measured calibration errors appear to result from stellar contamination, false object detection, the shear measurement method itself, selection bias and/or the use of biased weights. Additive systematics (false detections of shear) resulting from residual point-spread function anisotropy are, in most cases, reduced to below an equivalent shear of 0.001, an order of magnitude below cosmic shear distortions on the scales probed by current surveys. Our results provide a snapshot view of the accuracy of current ground-based weak lensing methods and a benchmark upon which we can improve. To this end we provide descriptions of each method tested and include details of the eight different implementations of the commonly used Kaiser, Squires and Broadhurst (1995) method (KSB+) to aid the improvement of future KSB+ analyses

Dissecting the thermal Sunyaev-Zeldovich-gravitational lensing cross-correlation with hydrodynamical simulations

We use the cosmo-OWLS suite of cosmological hydrodynamical simulations, which includes different galactic feedback models, to predict the cross-correlation signal between weak gravitational lensing and the thermal Sunyaev-Zeldovich (tSZ) y-parameter. The predictions are compared to the recent detection reported by van Waerbeke and collaborators. The simulations reproduce the weak lensing-tSZ cross-correlation, ξyk (θ), well. The uncertainty arising from different possible feedback models appears to be important on small scales only ($\theta \lesssim 10$ arcmin), while the amplitude of the correlation on all scales is sensitive to cosmological parameters that control the growth rate of structure (such as $\sigma_8$, $\Omega_m$ and $\Omega_b$). This study confirms our previous claim (in Ma et al.) that a significant proportion of the signal originates from the diffuse gas component in low-mass ($M_{\rm{halo}} \lesssim 10^{14} M_{\odot}$) clusters as well as from the region beyond the virial radius. We estimate that approximately 20$\%$ of the detected signal comes from low-mass clusters, which corresponds to about 30$\%$ of the baryon density of the Universe. The simulations also suggest that more than half of the baryons in the Universe are in the form of diffuse gas outside halos ($\gtrsim 5$ times the virial radius) which is not hot or dense enough to produce a significant tSZ signal or be observed by X-ray experiments. Finally, we show that future high-resolution tSZ-lensing cross-correlation observations will serve as a powerful tool for discriminating between different galactic feedback models

A Search for Warm/Hot Gas Filaments Between Pairs of SDSS Luminous Red Galaxies

We search the Planck data for a thermal Sunyaev-Zel'dovich (tSZ) signal due to gas filaments between pairs of Luminous Red Galaxies (LRG's) taken from the Sloan Digital Sky Survey Data Release 12 (SDSS/DR12). We identify $\sim$260,000 LRG pairs in the DR12 catalog that lie within 6-10 $h^{-1} \mathrm{Mpc}$ of each other in tangential direction and within 6 $h^{-1} \mathrm{Mpc}$ in radial direction. We stack pairs by rotating and scaling the angular positions of each LRG so they lie on a common reference frame, then we subtract a circularly symmetric halo from each member of the pair to search for a residual signal between the pair members. We find a statistically significant (5.3$\sigma$) signal between LRG pairs in the stacked data with a magnitude $\Delta y = (1.31 \pm 0.25) \times 10^{-8}$. The uncertainty is estimated from two Monte Carlo null tests which also establish the reliability of our analysis. Assuming a simple, isothermal, cylindrical filament model of electron over-density with a radial density profile proportional to $r_c/r$ (as determined from simulations), where $r$ is the perpendicular distance from the cylinder axis and $r_c$ is the core radius of the density profile, we constrain the product of over-density and filament temperature to be $\delta_c \times (T_{\rm e}/10^7 \, {\rm K}) \times (r_c/0.5h^{-1} \, {\rm Mpc}) = 2.7 \pm 0.5$. To our knowledge, this is the first detection of filamentary gas at over-densities typical of cosmological large-scale structure. We compare our result to the BAHAMAS suite of cosmological hydrodynamic simulations (McCarthy et al. 2017) and find a slightly lower, but marginally consistent Comptonization excess, $\Delta y = (0.84 \pm 0.24) \times 10^{-8}$

Next Generation Virgo Cluster Survey. XXI. The weak lensing masses of the CFHTLS and NGVS RedGOLD galaxy clusters and calibration of the optical richness

We measured stacked weak lensing cluster masses for a sample of 1323 galaxy clusters detected by the RedGOLD algorithm in the Canada–France–Hawaii Telescope Legacy Survey W1 and the Next Generation Virgo Cluster Survey at $0.2\lt z\lt 0.5$, in the optical richness range $10\lt \lambda \lt 70$. This is the most comprehensive lensing study of a $\sim 100 \%$ complete and $\sim 80 \%$ pure optical cluster catalog in this redshift range. We test different mass models, and our final model includes a basic halo model with a Navarro Frenk and White profile, as well as correction terms that take into account cluster miscentering, non-weak shear, the two-halo term, the contribution of the Brightest Cluster Galaxy, and an a posteriori correction for the intrinsic scatter in the mass–richness relation. With this model, we obtain a mass–richness relation of $\mathrm{log}{M}_{200}/{M}_{\odot }\,=(14.46\pm 0.02)+(1.04\pm 0.09)\mathrm{log}(\lambda /40)$ (statistical uncertainties). This result is consistent with other published lensing mass–richness relations. We give the coefficients of the scaling relations between the lensing mass and X-ray mass proxies, L X and T X, and compare them with previous results. When compared to X-ray masses and mass proxies, our results are in agreement with most previous results and simulations, and consistent with the expected deviations from self-similarity

CFHTLenS: The Environmental Dependence of Galaxy Halo Masses from Weak Lensing

We use weak gravitational lensing to analyse the dark matter halos around satellite galaxies in galaxy groups in the CFHTLenS dataset. This dataset is derived from the CFHTLS-Wide survey, and encompasses 154 sq. deg of high-quality shape data. Using the photometric redshifts, we divide the sample of lens galaxies with stellar masses in the range 10^9 Msun to 10^10.5 Msun into those likely to lie in high-density environments (HDE) and those likely to lie in low-density environments (LDE). Through comparison with galaxy catalogues extracted from the Millennium Simulation, we show that the sample of HDE galaxies should primarily (~61%) consist of satellite galaxies in groups, while the sample of LDE galaxies should consist of mostly (~87%) non-satellite (field and central) galaxies. Comparing the lensing signals around samples of HDE and LDE galaxies matched in stellar mass, the lensing signal around HDE galaxies clearly shows a positive contribution from their host groups on their lensing signals at radii of ~500--1000 kpc, the typical separation between satellites and group centres. More importantly, the subhalos of HDE galaxies are less massive than those around LDE galaxies by a factor 0.65 +/- 0.12, significant at the 2.9 sigma level. A natural explanation is that the halos of satellite galaxies are stripped through tidal effects in the group environment. Our results are consistent with a typical tidal truncation radius of ~40 kpc

CFHTLenS: Cosmological constraints from a combination of cosmic shear two-point and three-point correlations

Higher-order, non-Gaussian aspects of the large-scale structure carry valuable information on structure formation and cosmology, which is complementary to second-order statistics. In this work we measure second- and third-order weak-lensing aperture-mass moments from CFHTLenS and combine those with CMB anisotropy probes. The third moment is measured with a significance of $2\sigma$. The combined constraint on $\Sigma_8 = \sigma_8 (\Omega_{\rm m}/0.27)^\alpha$ is improved by 10%, in comparison to the second-order only, and the allowed ranges for $\Omega_{\rm m}$ and $\sigma_8$ are substantially reduced. Including general triangles of the lensing bispectrum yields tighter constraints compared to probing mainly equilateral triangles. Second- and third-order CFHTLenS lensing measurements improve Planck CMB constraints on $\Omega_{\rm m}$ and $\sigma_8$ by 26% for flat $\Lambda$CDM. For a model with free curvature, the joint CFHTLenS-Planck result is $\Omega_{\rm m} = 0.28 \pm 0.02$ (68% confidence), which is an improvement of 43% compared to Planck alone. We test how our results are potentially subject to three astrophysical sources of contamination: source-lens clustering, the intrinsic alignment of galaxy shapes, and baryonic effects. We explore future limitations of the cosmological use of third-order weak lensing, such as the nonlinear model and the Gaussianity of the likelihood function