Including massive neutrinos in thermal Sunyaev Zeldovich power spectrum and cluster counts analyses

We consistently include the effect of massive neutrinos in the thermal Sunyaev Zeldovich (SZ) power spectrum and cluster counts analyses, highlighting subtle dependencies on the total neutrino mass and data combination. In particular, we find that using the transfer functions for cold dark matter (CDM) + baryons in the computation of the halo mass function, instead of the transfer functions including neutrino perturbations, as prescribed in recent work, yields an ≈0.25 per cent downward shift of the σ_8 constraint from tSZ power spectrum data, with a fiducial neutrino mass Σm_ν = 0.06 eV. In ΛCDM, with an X-ray mass bias corresponding to the expected hydrostatic mass bias, i.e. (1 − b) ≃ 0.8, our constraints from Planck SZ data are consistent with the latest results from SPT, DES-Y1, and KiDS+VIKING-450. In νΛCDM, our joint analyses of Planck SZ with Planck 2015 primary CMB yield a small improvement on the total neutrino mass bound compared to the Planck 2015 primary CMB constraint, as well as (1 − b) = 0.64 ± 0.04 (68 per cent CL). For forecasts, we find that competitive neutrino mass measurements using cosmic variance limited SZ power spectrum require masking the heaviest clusters and probing the small-scale SZ power spectrum up to ℓ_max ≈ 10^4. Although this is challenging, we find that SZ power spectrum can realistically be used to tightly constrain intracluster medium properties: we forecast a 2 per cent determination of the X-ray mass bias by combining CMB-S4 and our mock SZ power spectrum with ℓ_max = 10^3.