The clustering amplitude of 7143 clusters from the Sloan Digital Sky Survey (SDSS) is found to increase with cluster mass, closely agreeing with the Gaussian random field hypothesis for structure formation. The amplitude of the observed cluster correlation exceeds the predictions from pure cold dark matter (CDM) simulation by $\simeq 6\%$ for the standard Planck-based values of the cosmological parameters. We show that this excess can be naturally accounted for by free streaming of light neutrinos, which opposes gravitational growth, so clusters formed at fixed mass are fewer and hence more biased than for a pure CDM density field. An enhancement of the cluster bias by 7\% matches the observations, corresponding to a total neutrino mass, $m_{\nu} = 0.119 \pm 0.034$ eV at 67\% confidence level, for the standard relic neutrino density. If ongoing laboratory experiments favor a normal neutrino mass hierarchy then we may infer a somewhat larger total mass than the minimum oscillation based value, $\sum m_{\nu} \simeq 0.056eV$, with 90\% confidence. Much higher precision can be achieved by applying our method to a larger sample of more distant clusters with weak lensing derived masses.

We construct a large, redshift complete sample of distant galaxy clusters by correlating Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12) redshifts with clusters identified with the red-sequence Matched-filter Probabilistic Percolation (redMaPPer) algorithm. Our spectroscopic completeness is 97% for ~ 7,000 clusters within the redMaPPer selection limit, $z \leqslant$ 0.325, so that our cluster correlation functions are much more precise than earlier work and not suppressed by photometric redshifts. We derive an accurate power-law mass-richness relation from the observed abundance with respect to the mass function from Millennium XXL (MXXL) simulation, adjusted to the Planck weighted cosmology. The number density of clusters is found to decline by 20% over the range 0.1 $< z <$ 0.3, in good agreement with the evolution predicted by MXXL. Our projected three-dimensional correlation function scales with richness, $\lambda$, rising from $r_0=$ 14 h$^{-1}$ Mpc at $\lambda\simeq$ 25, to $r_0=$ 22 h$^{-1}$ Mpc at $\lambda\simeq$ 60, with a gradient that matches MXXL when applying our mass-richness relation, whereas the observed amplitude of the correlation function at $\left<z\right>=$ 0.24 exceeds the MXXL prediction by 20% at the $\simeq$ 2.5$\sigma$ level. This tension cannot be blamed on spurious, randomly located clusters as this would reduce the correlation amplitude. Full consistency between the correlation function and the abundances is achievable for the pre-Planck values of $\sigma_8=$ 0.9, $\Omega_m=$ 0.25, and h = 0.73, matching the improved distance ladder estimate of the Hubble constant.

Using Planck satellite data, we construct SZ gas pressure profiles for a large, volume-complete sample of optically selected clusters. We have defined a sample of over 8,000 redMaPPer clusters from the Sloan Digital Sky Survey (SDSS), within the volume-complete redshift region 0.100 < z < 0.325, for which we construct Sunyaev-Zel'dovich (SZ) effect maps by stacking Planck data over the full range of richness. Dividing the sample into richness bins we simultaneously solve for the mean cluster mass in each bin together with the corresponding radial pressure profile parameters, employing an MCMC analysis. These profiles are well detected over a much wider range of cluster mass and radius than previous work, showing a clear trend towards larger break radius with increasing cluster mass. Our SZ-based masses fall ~24% below the mass-richness relations from weak lensing, in a similar fashion as the "hydrostatic bias" related with X-ray derived masses. We correct for this bias to derive an optimal mass-richness relation finding a slope 1.22 +/- 0.04 and a pivot mass log(M_500/M_0)= 14.432 +/- 0.041, evaluated at a richness lambda=60. Finally, we derive a tight Y_500-M_500 relation over a wide range of cluster mass, with a power law slope equal to 1.72 +/- 0.07, that agrees well with the independent slope obtained by the Planck team with an SZ-selected cluster sample, but extends to lower masses with higher precision.

A clean measurement of the evolution of the galaxy cluster mass function can significantly improve our understanding of cosmology from the rapid growth of cluster masses below z < 0.5. Here we examine the consistency of cluster catalogues selected from the SDSS by applying two independent gravity-based methods using all available spectroscopic redshifts from the DR10 release. First, we detect a gravitational redshift related signal for 20,119 and 13,128 clusters with spectroscopic redshifts contained in the GMBCG and redMaPPer catalogues, respectively, at a level of $\sim - 10$ km s$^{-1}$. This we show is consistent with the magnitude expected using the richness-mass relations provided by the literature and after applying recently clarified relativistic and flux bias corrections. This signal is also consistent with the richest clusters in the larger catalogue of Wen et al. (2012), corresponding to $M_{200m} \gtrsim 2 \times 10^{14}\,\mathrm{M}_\odot\,h^{-1}$, however we find no significant detection of gravitational redshift signal for less riched clusters, which may be related to bulk motions from substructure and spurious cluster detections. Second, we find all three catalogues generate mass-dependent levels of lensing magnification bias, which enhances the mean redshift of flux-selected background galaxies from the BOSS survey. The magnitude of this lensing effect is generally consistent with the corresponding richness-mass relations advocated for the surveys. We conclude that all catalogues comprise a high proportion of reliable clusters, and that the GMBCG and redMaPPer cluster finder algorithms favor more relaxed clusters with a meaningful gravitational redshift signal, as anticipated by the red-sequence colour selection of the GMBCG and redMaPPer samples.