Issue 
A&A
Volume 571, November 2014
Planck 2013 results



Article Number  A20  
Number of page(s)  20  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201321521  
Published online  29 October 2014 
Planck 2013 results. XX. Cosmology from Sunyaev–Zeldovich cluster counts
^{1}
APC, AstroParticule et Cosmologie, Université Paris Diderot,
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris Cité, 10 rue Alice Domon et Léonie
Duquet, 75205
Paris Cedex 13,
France
^{2}
Aalto University Metsähovi Radio Observatory,
Metsähovintie 114, 02540
Kylmälä,
Finland
^{3}
Academy of Sciences of Tatarstan, Bauman Str., 20, Kazan, 420111
Republic of Tatarstan,
Russia
^{4}
African Institute for Mathematical Sciences,
68 Melrose Road, Muizenberg,
7701
Rondebosch Cape Town, South
Africa
^{5}
Agenzia Spaziale Italiana Science Data Center, c/o ESRIN, via
Galileo Galilei, 00044
Frascati,
Italy
^{6}
Agenzia Spaziale Italiana, Viale Liegi 26,
00198
Roma,
Italy
^{7}
Astrophysics Group, Cavendish Laboratory, University of
Cambridge, J J Thomson
Avenue, Cambridge
CB3 0HE,
UK
^{8}
Astrophysics & Cosmology Research Unit, School of
Mathematics, Statistics & Computer Science, University of
KwaZuluNatal, Westville Campus,
Private Bag X54001, 4000
Durban, South
Africa
^{9}
Atacama Large Millimeter/submillimeter Array, ALMA Santiago Central
Offices, Alonso de Cordova 3107,
Vitacura, Casilla
763 0355
Santiago,
Chile
^{10}
CITA, University of Toronto, 60 St. George St., Toronto, ON
M5S 3H8,
Canada
^{11}
CNRS, IRAP, 9 Av.
colonel Roche, BP
44346, 31028
Toulouse Cedex 4,
France
^{12}
California Institute of Technology, Pasadena, California, USA
^{13}
Centre for Theoretical Cosmology, DAMTP, University of Cambridge,
Wilberforce Road, Cambridge
CB3 0WA UK
^{14}
Centro de Astrofísica, Universidade do Porto,
Rua das Estrelas, 4150762
Porto,
Portugal
^{15}
Centro de Estudios de Física del Cosmos de Aragón (CEFCA), Plaza San
Juan, 1, planta 2, 44001
Teruel,
Spain
^{16}
Computational Cosmology Center, Lawrence Berkeley National
Laboratory, Berkeley,
California,
USA
^{17}
Consejo Superior de Investigaciones Científicas
(CSIC), 28006
Madrid,
Spain
^{18}
DSM/Irfu/SPP, CEASaclay, 91191
GifsurYvette Cedex,
France
^{19}
DTU Space, National Space Institute, Technical University of
Denmark, Elektrovej
327, 2800
Kgs. Lyngby,
Denmark
^{20}
Département de Physique Théorique, Université de
Genève, 24 quai E.
Ansermet, 1211
Genève 4,
Switzerland
^{21}
Departamento de Física Fundamental, Facultad de Ciencias,
Universidad de Salamanca, 37008
Salamanca,
Spain
^{22}
Departamento de Física, Universidad de Oviedo,
Avda. Calvo Sotelo s/n,
33007
Oviedo,
Spain
^{23}
Department of Astronomy and Astrophysics, University of
Toronto, 50 Saint George Street,
Toronto, Ontario,
Canada
^{24}
Department of Astronomy and Geodesy, Kazan Federal
University, Kremlevskaya Str.,
18, 420008
Kazan,
Russia
^{25}
Department of Astrophysics/IMAPP, Radboud University
Nijmegen, P.O. Box
9010, 6500 GL
Nijmegen, The
Netherlands
^{26}
Department of Electrical Engineering and Computer Sciences,
University of California, Berkeley, California, USA
^{27}
Department of Physics & Astronomy, University of British
Columbia, 6224 Agricultural Road,
Vancouver, British
Columbia, Canada
^{28}
Department of Physics and Astronomy, Dana and David Dornsife College
of Letter, Arts and Sciences, University of Southern California,
Los Angeles, CA
90089,
USA
^{29}
Department of Physics and Astronomy, University College
London, London
WC1E 6BT,
UK
^{30}
Department of Physics and Astronomy, University of
Sussex, Brighton
BN1 9QH,
UK
^{31}
Department of Physics, Gustaf Hällströmin katu 2a, University of
Helsinki, 00014
Helsinki,
Finland
^{32}
Department of Physics, Princeton University,
Princeton, New Jersey, USA
^{33}
Department of Physics, University of California,
Berkeley, California, USA
^{34}
Department of Physics, University of California,
One Shields Avenue, Davis, California, USA
^{35}
Department of Physics, University of California,
Santa Barbara, California, USA
^{36}
Department of Physics, University of Illinois at Urbana
Champaign, 1110 West Green
Street, Urbana,
Illinois,
USA
^{37}
Dipartimento di Fisica e Astronomia G. Galilei, Università degli
Studi di Padova, via Marzolo
8, 35131
Padova,
Italy
^{38}
Dipartimento di Fisica e Scienze della Terra, Università di
Ferrara, Via Saragat
1, 44122
Ferrara,
Italy
^{39}
Dipartimento di Fisica, Università La Sapienza,
P. le A. Moro 2, 00185
Roma,
Italy
^{40}
Dipartimento di Fisica, Università degli Studi di
Milano, Via Celoria,
16, 20133
Milano,
Italy
^{41}
Dipartimento di Fisica, Università degli Studi di
Trieste, via A. Valerio
2, 34127
Trieste,
Italy
^{42}
Dipartimento di Fisica, Università di Roma Tor
Vergata, via della Ricerca Scientifica,
1, 00133
Roma,
Italy
^{43}
Discovery Center, Niels Bohr Institute, Blegdamsvej 17, 2100
Copenhagen,
Denmark
^{44}
Dpto. Astrofísica, Universidad de La Laguna (ULL),
38206, La Laguna, Tenerife, Spain
^{45}
European Southern Observatory, ESO Vitacura, Alonso de Cordova 3107, Vitacura, Casilla
19001
Santiago,
Chile
^{46}
European Space Agency, ESAC, Planck Science Office, Camino bajo del
Castillo, s/n, Urbanización Villafranca del Castillo, 28691 Villanueva de la
Cañada, Madrid,
Spain
^{47}
European Space Agency, ESTEC, Keplerlaan 1,
2201 AZ
Noordwijk, The
Netherlands
^{48}
Finnish Centre for Astronomy with ESO (FINCA), University of
Turku, Väisäläntie
20, 21500
Piikkiö,
Finland
^{49}
Haverford College Astronomy Department, 370 Lancaster Avenue, Haverford, Pennsylvania, USA
^{50}
Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University
of Helsinki, 00014
Helsinki,
Finland
^{51}
INAF – Osservatorio Astrofisico di Catania, Via S. Sofia
78, 95123
Catania,
Italy
^{52}
INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio
5, 35122
Padova,
Italy
^{53}
INAF – Osservatorio Astronomico di Roma, via di Frascati
33, 00040
Monte Porzio Catone,
Italy
^{54}
INAF – Osservatorio Astronomico di Trieste, via G.B. Tiepolo
11, Trieste,
Italy
^{55}
INAF/IASF Bologna, via Gobetti 101, 40129
Bologna,
Italy
^{56}
INAF/IASF Milano, via E. Bassini 15, 20133
Milano,
Italy
^{57}
INFN, Sezione di Bologna, via Irnerio 46,
40126
Bologna,
Italy
^{58}
INFN, Sezione di Roma 1, Università di Roma Sapienza,
Piazzale Aldo Moro 2,
00185
Roma,
Italy
^{59}
IPAG: Institut de Planétologie et d’Astrophysique de Grenoble,
Université Joseph Fourier, Grenoble 1/CNRSINSU, UMR 5274, 38041
Grenoble,
France
^{60}
ISDC Data Centre for Astrophysics, University of
Geneva, Ch. d’Ecogia
16, 1290
Versoix,
Switzerland
^{61}
IUCAA, Post Bag 4, Ganeshkhind, Pune University
Campus, 411 007
Pune,
India
^{62}
Imperial College London, Astrophysics group, Blackett
Laboratory, Prince Consort
Road, London,
SW7 2AZ,
UK
^{63}
Infrared Processing and Analysis Center, California Institute of
Technology, Pasadena,
CA
91125,
USA
^{64}
Institut Néel, CNRS, Université Joseph Fourier Grenoble
I, 25 rue des
Martyrs, 38042
Grenoble,
France
^{65}
Institut Universitaire de France, 103 bd SaintMichel, 75005
Paris,
France
^{66}
Institut d’Astrophysique Spatiale, CNRS (UMR 8617) Université
ParisSud 11, Bâtiment
121, 91405
Orsay,
France
^{67}
Institut d’Astrophysique de Paris, CNRS (UMR 7095),
98bis boulevard Arago,
75014
Paris,
France
^{68}
Institute for Space Sciences, 07745
BucharestMagurale,
Romania
^{69}
Institute of Astronomy and Astrophysics, Academia
Sinica, 10617
Taipei,
Taiwan
^{70}
Institute of Astronomy, University of Cambridge,
Madingley Road, Cambridge
CB3 0HA,
UK
^{71}
Institute of Theoretical Astrophysics, University of
Oslo, Blindern,
0315
Oslo,
Norway
^{72}
Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, 38200, La Laguna, Tenerife, Spain
^{73}
Instituto de Física de Cantabria (CSICUniversidad de
Cantabria), Avda. de los Castros
s/n, 39005
Santander,
Spain
^{74}
Jet Propulsion Laboratory, California Institute of
Technology, 4800 Oak Grove
Drive, Pasadena,
California,
USA
^{75}
Jodrell Bank Centre for Astrophysics, Alan Turing Building, School
of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13
9PL, UK
^{76}
Kavli Institute for Cosmology Cambridge,
Madingley Road, Cambridge, CB3 0HA, UK
^{77}
LAL, Université ParisSud, CNRS/IN2P3, 91898
Orsay,
France
^{78}
LERMA, CNRS, Observatoire de Paris, 61 avenue de
l’Observatoire, 75014
Paris,
France
^{79}
Laboratoire AIM, IRFU/Service d’Astrophysique – CEA/DSM – CNRS –
Université Paris Diderot, Bât. 709, CEASaclay, 91191
GifsurYvette Cedex,
France
^{80}
Laboratoire Traitement et Communication de l’Information, CNRS (UMR
5141) and Télécom ParisTech, 46 rue
Barrault, 75634
Paris Cedex 13,
France
^{81}
Laboratoire de PhysiqueSubatomique et de Cosmologie, Université
Joseph Fourier Grenoble I, CNRS/IN2P3, Institut National Polytechnique de
Grenoble, 53 rue des
Martyrs, 38026
Grenoble Cedex,
France
^{82}
Laboratoire de Physique Théorique, Université ParisSud 11 &
CNRS, Bâtiment 210,
91405
Orsay,
France
^{83}
Lawrence Berkeley National Laboratory, Berkeley, California, USA
^{84}
MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStr. 1, 85741
Garching,
Germany
^{85}
MaxPlanckInstitut für Extraterrestrische Physik,
Giessenbachstraße, 85748
Garching,
Germany
^{86}
McGill Physics, Ernest Rutherford Physics Building, McGill
University, 3600 rue University, Montréal, QC,
H3A 2T8,
Canada
^{87}
MilliLab, VTT Technical Research Centre of Finland, Tietotie
3, Espoo,
Finland
^{88}
Moscow Institute of Physics and Technology, Institutsky
per., 9,
141700
Dolgoprudny,
Russia
^{89}
Niels Bohr Institute, Blegdamsvej 17, 2100
Copenhagen,
Denmark
^{90}
Observational Cosmology, Mail Stop 36717, California Institute of
Technology, Pasadena,
CA, 91125, USA
^{91}
Optical Science Laboratory, University College London,
Gower Street, London, UK
^{92}
SBITPLPPC, EPFL, CH1015
Lausanne,
Switzerland
^{93}
SISSA, Astrophysics Sector, via Bonomea 265,
34136
Trieste,
Italy
^{94}
SUPA, Institute for Astronomy, University of Edinburgh, Royal
Observatory, Blackford
Hill, Edinburgh
EH9 3HJ,
UK
^{95}
School of Physics and Astronomy, Cardiff University,
Queens Buildings, The Parade,
Cardiff, CF24 3AA, UK
^{96}
School of Physics and Astronomy, University of
Nottingham, Nottingham
NG7 2RD,
UK
^{97}
Space Research Institute (IKI), Russian Academy of
Sciences, Profsoyuznaya Str,
84/32, 117997
Moscow,
Russia
^{98}
Space Sciences Laboratory, University of California,
Berkeley, California, USA
^{99}
Special Astrophysical Observatory, Russian Academy of
Sciences, Nizhnij Arkhyz,
Zelenchukskiy region, 369167
KarachaiCherkessian Republic,
Russia
^{100}
Stanford University, Dept of Physics, Varian Physics Bldg, 382 Via Pueblo
Mall, Stanford,
California,
USA
^{101}
SubDepartment of Astrophysics, University of Oxford,
Keble Road, Oxford
OX1 3RH,
UK
^{102}
Tübitak National Observatory, Akdeniz University
Campus, 07058
Antalya,
Turkey
^{103}
Theory Division, PHTH, CERN, 1211
Geneva 23,
Switzerland
^{104}
UPMC Univ Paris 06, UMR7095, 98bis boulevard Arago, 75014
Paris,
France
^{105}
Université de Toulouse, UPSOMP, IRAP, 31028
Toulouse Cedex 4,
France
^{106}
University Observatory, Ludwig Maximilian University of
Munich, Scheinerstrasse
1, 81679
Munich,
Germany
^{107}
University of Granada, Departamento de Física Teórica y del Cosmos,
Facultad de Ciencias, 18071
Granada,
Spain
^{108}
University of Miami, Knight Physics Building, 1320 Campo Sano Dr.,
Coral Gables, Florida, USA
^{109}
Warsaw University Observatory, Aleje Ujazdowskie 4, 00478
Warszawa,
Poland
Received:
20
March
2013
Accepted:
4
February
2014
We present constraints on cosmological parameters using number counts as a function of redshift for a subsample of 189 galaxy clusters from the Planck SZ (PSZ) catalogue. The PSZ is selected through the signature of the SunyaevZeldovich (SZ) effect, and the subsample used here has a signaltonoise threshold of seven, with each object confirmed as a cluster and all but one with a redshift estimate. We discuss the completeness of the sample and our construction of a likelihood analysis. Using a relation between mass M and SZ signal Y calibrated to Xray measurements, we derive constraints on the power spectrum amplitude σ_{8} and matter density parameter Ω_{m} in a flat ΛCDM model. We test the robustness of our estimates and find that possible biases in the Y–M relation and the halo mass function are larger than the statistical uncertainties from the cluster sample. Assuming the Xray determined mass to be biased low relative to the true mass by between zero and 30%, motivated by comparison of the observed mass scaling relations to those from a set of numerical simulations, we find that σ_{8} = 0.75 ± 0.03, Ω_{m} = 0.29 ± 0.02, and σ_{8}(Ω_{m}/ 0.27)^{0.3} = 0.764 ± 0.025. The value of σ_{8} is degenerate with the mass bias; if the latter is fixed to a value of 20% (the central value from numerical simulations) we find σ_{8}(Ω_{m}/0.27)^{0.3} = 0.78 ± 0.01 and a tighter onedimensional range σ_{8} = 0.77 ± 0.02. We find that the larger values of σ_{8} and Ω_{m} preferred by Planck’s measurements of the primary CMB anisotropies can be accommodated by a mass bias of about 40%. Alternatively, consistency with the primary CMB constraints can be achieved by inclusion of processes that suppress power on small scales relative to the ΛCDM model, such as a component of massive neutrinos. We place our results in the context of other determinations of cosmologicalparameters, and discuss issues that need to be resolved in order to make further progress in this field.
Key words: cosmological parameters / largescale structure of Universe / galaxies: clusters: general
© ESO, 2014
1. Introduction
This paper, one of a set associated with the 2013 release of data from the Planck^{1} mission (Planck Collaboration I 2014), describes the constraints on cosmological parameters using number counts as a function of redshift for a sample of 189 galaxy clusters.
Within the standard picture of structure formation, galaxies aggregate into clusters of galaxies at late times, forming bound structures at locations where the initial fluctuations create the deepest potential wells. The study of these galaxy clusters has played a significant role in the development of cosmology over many years (see, for example, Perrenod 1980; Oukbir & Blanchard 1992; White et al. 1993; Carlberg et al. 1996; Voit 2005; Henry et al. 2009; Vikhlinin et al. 2009b; Allen et al. 2011a). More recently, as samples of clusters have increased in size and variety, number counts inferred from tightlyselected surveys have been used to obtain detailed constraints on the cosmological parameters.
The early galaxy cluster catalogues were constructed by eye from photographic plates with a “richness” (or number of galaxies) attributed to each cluster (Abell 1958; Abell et al. 1989). As time has passed, new approaches for selecting clusters have been developed, most notably using Xray emission due to thermal Bremsstrahlung radiation from the hot gas that makes up most of the baryonic matter in the cluster. Xray cluster surveys include both the NORAS (Böhringer et al. 2000) and REFLEX (Böhringer et al. 2004) surveys, based on ROSAT satellite observations, which have been used as source catalogues for higherprecision observations by the Chandra and XMMNewton satellites, as well as surveys with XMMNewton, including the XMM Cluster Survey (XCS, Mehrtens et al. 2012) and the XMM Large Scale Structure survey (XMMLSS, Willis et al. 2013).
To exploit clusters for cosmology, a key issue is how the properties used to select and characterize the cluster are related to the total mass of the cluster, since this is the quantity most readily predicted using theoretical models. Galaxies account for a small fraction of the cluster mass, and the scatter between richness and mass appears to be large. However, there are a number of other possibilities. In particular, there are strong correlations between the total mass and both the integrated Xray surface brightness and Xray temperature, making them excellent mass proxies.
The SunyaevZeldovich (SZ) effect (Sunyaev & Zeldovich 1970; Zeldovich & Sunyaev 1969) is the inverse Compton scattering of cosmic microwave background (CMB) photons by the hot gas along the line of sight, and this is most significant when the line of sight passes through a galaxy cluster. It leads to a decrease in the overall brightness temperature in the RayleighJeans portion of the spectrum and an increase in the Wien tail, with a null around 217 GHz (see Birkinshaw 1999 for a review). The amplitude of the SZ effect is given by the integrated pressure of the gas within the cluster along the line of sight. Evidence both from observation (Marrone et al. 2012; Planck Collaboration Int. III 2013) and from numerical simulations (Springel et al. 2001; da Silva et al. 2004; Motl et al. 2005; Nagai 2006; Kay et al. 2012) suggests that the SZ effect is an excellent mass proxy. A number of articles have discussed the possibility of using SZselected cluster samples to constrain cosmological parameters (Barbosa et al. 1996; Aghanim et al. 1997; Haiman et al. 2001; Holder et al. 2001; Weller et al. 2002; Diego et al. 2002; Battye & Weller 2003).
Fig. 1 Distribution on the sky of the Planck SZ cluster subsample used in this paper, with the 35% mask overlaid. 
This paper describes the constraints on cosmological parameters imposed by a high signaltonoise (S/N) subsample of the Planck SZ Catalogue (PSZ, see Planck Collaboration XXIX 2014, henceforth Paper I, for details of the entire catalogue) containing nearly 200 clusters (shown in Fig. 1). This subsample has been selected to be pure, in the sense that all the objects within it have been confirmed as clusters via additional observations, either from the literature or undertaken by the Planck collaboration. In addition all objects but one have a measured redshift, either photometric or spectroscopic. This is the largest SZselected sample of clusters used to date for this purpose. We will show that it is the systematic uncertainties from our imperfect knowledge of cluster properties that dominate the overall uncertainty on cosmological constraints.
The Planck cluster sample is complementary to those from observations using the South Pole Telescope (SPT, Carlstrom et al. 2011) and the Atacama Cosmology Telescope (ACT, Swetz et al. 2011), whose teams recently published the first large samples of SZselected clusters (Reichardt et al. 2013; Hasselfield et al. 2013). The resolution of Planck at the relevant frequencies is between 5 and 10 arcmin, whereas that for ACT and SPT is about 1 arcmin, but the Planck sky coverage is much greater. This means that Planck typically finds larger, more massive, and lowerredshift clusters than those found by SPT and ACT.
Our strategy is to focus on number counts of clusters, as a function of redshift, above a high S/N threshold of seven and to explore the robustness of the results. We do not use the observed SZ brightness of the clusters, due to the significant uncertainty caused by the sizeflux degeneracy as discussed in Paper I. Accordingly, our theoretical modelling of the cluster population is directed only at determining the expected number of clusters in each redshift bin exceeding the S/N threshold. The predicted and observed numbers of clusters are then compared in order to obtain the likelihood. In the future, we will make use of the SZestimated mass and a larger cluster sample to extend the analysis to broader cosmological scenarios.
This paper is laid out as follows. We describe the theoretical modelling of the redshift number counts in Sect. 2, while Sect. 3 presents the PSZ cosmological sample and selection function used in this work. The likelihood we adopt for putting constraints on cosmological parameters is given in Sect. 4. Section 5 presents our results on cosmological parameter estimation and assesses their robustness. We discuss how they fit in with other cluster and cosmological constraints in Sect. 6, before providing a final summary. A detailed discussion of our calibration of the SZ flux versus mass relation and its uncertainties is given in Appendix A.
2. Modelling cluster number counts
2.1. Model definitions
We parameterize the standard cosmological model as follows. The densities of various components are specified relative to the presentday critical density, with Ω_{X} = ρ_{X}/ρ_{crit} denoting that for component X. These components always include matter, Ω_{m}, and a cosmological constant Ω_{Λ}. For this work we assume that the Universe is flat, that is, Ω_{m} + Ω_{Λ} = 1, and the optical depth to reionization is fixed at τ = 0.085 except in the CMB and SZ analyses. The presentday expansion rate of the Universe is quantified by the Hubble constant H_{0} = 100 h km s^{1} Mpc^{1}.
The cluster number counts are very sensitive to the amplitude of the matter power spectrum. When studying cluster counts it is usual to parametrize this in terms of the density variance in spheres of radius 8 h^{1} Mpc, denoted σ_{8}, rather than overall power spectrum amplitude, A_{s}. In cases where we include primary CMB data we use A_{s} and compute σ_{8} as a derived parameter. In addition to the parameters above, we allow the other standard cosmological parameters to vary: n_{s} representing the spectral index of density fluctuations; and Ω_{b}h^{2} quantifying the baryon density.
The number of clusters predicted to be observed by a survey in a given redshift interval [ z_{i},z_{i + 1} ] can be written (1)with (2)where dΩ is the solid angle element and M_{500} is the mass within the radius where the mean enclosed density is 500 times the critical density. The quantity is the survey completeness at a given location (l,b) on the sky, given by (3)Here P(z,M_{500}  Y_{500},θ_{500}) is the distribution of (z,M_{500}) for a given (Y_{500},θ_{500}), where Y_{500} and θ_{500} are the SZ flux and size of a cluster of redshift and mass (z,M_{500}).
This distribution is obtained from the scaling relations between Y_{500}, θ_{500}, and M_{500}, discussed later in this section. Note that depends on cosmological parameters through P(z,M_{500}  Y_{500},θ_{500}), while the completeness in terms of the observables, χ(Y_{500},θ_{500},l,b), does not depend on the cosmology as it refers directly to the observed quantities.
Fig. 2 Blind S/N versus S/N reextracted at the Xray position using the Xray size, for the MMF3 detections of Planck clusters that are associated with known Xray clusters in the reference cosmological sample. In contrast to the blind SZ flux, the blind S/N is in good agreement with S/N measured using Xray priors. 
For the present work, we restrict our analysis to the quantity dN/dz that measures the total counts in redshift bins. In particular, we do not use the blind SZ flux estimated by the cluster candidate extraction methods that, as detailed in Planck Collaboration VIII (2011), is found to be significantly higher than the flux predicted from Xray measurements. In contrast to the blind SZ flux, the blind S/N is in good agreement with the S/N measured using Xray priors. Figure 2 shows the blind S/N (S/N_{blind}) versus the S/N reextracted at the Xray position and using the Xray size (S/N_{X}). The clusters follow the equality line. In Sect. 3, we use the (S/N_{blind}) values to define our cosmological sample, while for the predicted counts (defined in Sect. 2) we use the completeness based on S/N_{X}. Our analysis relies on the good match between these two quantities^{2}.
To carry out a prediction of the counts expected in a survey, given cosmological assumptions, we therefore need the following inputs:

a mass function that tells us the number distribution of clusterswith mass and redshift;

scaling relations that can predict observable quantities from the mass and redshift;

the completeness of the survey in terms of those observables, which tells us the probability that a model cluster would make it into the survey catalogue.
These are described in the remainder of this section and in the next.
2.2. Mass function
Our main results use the mass function from Tinker et al. (2008), giving the number of haloes per unit volume: (4)where (5)and ρ_{m}(z = 0) is the mean matter density at z = 0. The coefficients A, a, b and c are tabulated in Tinker et al. (2008) for different overdensities, Δ_{mean}, with respect to the mean cosmic density, and depend on z. Here we use Δ_{critical} = 500 relative to the critical density, so we compute the relevant mass function coefficients by interpolating the Tinker et al. (2008) tables for haloes with Δ_{mean} ≡ Δ_{critical}/ Ω_{m}(z) = 500/Ω_{m}(z), where Ω_{m}(z) is the matter density parameter at redshift z.
The quantity σ is the standard deviation, computed in linear perturbation theory, of the density perturbations in a sphere of radius R, which is related to the mass by M = 4πρ_{m}(z = 0)R^{3}/ 3. It is given by (6)where P(k,z) is the matter power spectrum at redshift z, which we compute for any given set of cosmological parameters using CAMB (Lewis et al. 2000), and W(x) = 3(sinx − xcosx) /x^{3} is the filter function of a spherical top hat of radius R.
The quantity dN/(dz dM_{500} dΩ) in Eq. (2) is computed by multiplying the mass function dN(M_{500},z)/dM_{500} by the volume element dV/(dz dΩ).
As a baseline we use, except where stated otherwise, the Tinker et al. (2008) mass function, but we consider an alternative mass function as a crosscheck. In a recent publication by Watson et al. (2013), a new mass function is extracted from the combination of large cosmological simulations (typical particle numbers of 3000^{3} to 6000^{3}) with a very large dynamic range (size from 11 h^{1} to 6000 h^{1} Mpc), which extends the maximum volume probed by Tinker et al. by two orders of magnitude. The two mass functions agree fairly well, except in the case of the most massive objects, where Tinker et al.’s mass function predicts more clusters than Watson et al.’s. The Tinker et al. mass function might be derived from volumes that are not large enough to properly sample the rarer clusters. These rare clusters are more relevant for Planck than for groundbased SZ experiments, which probe smaller areas of the sky. The Watson et al. mass function is used only in Sect. 5.3, which deals with mass function uncertainties.
2.3. Scaling relations
A key issue is to relate the observed SZ flux, Y_{500}, to the mass M_{500} of the cluster. As we show in Sect. 5, cosmological constraints are sensitive to the normalization, slope and scatter of the assumed Y_{500}–M_{500} relation. We thus paid considerable attention to deriving the most accurate scaling relations possible, with careful handling of statistical and systematic uncertainties, and to testing their impact on the derived cosmological parameters.
The baseline relation is obtained from an observational calibration of the Y_{500}–M_{500} relation on onethird of the cosmological sample. The calibration uses , the mass derived from the Xray Y_{X}–M_{500} relation, as a mass proxy. Here Y_{X} is the Xray analogue of the SZ signal introduced by Kravtsov et al. (2006), as defined in Appendix A. Y_{500} is then measured interior to , the radius corresponding to . The mean bias between and the true mass, (1 − b), is assumed to account for all possible observational biases (departure from hydrostatic equilibrium (HE), absolute instrument calibration, temperature inhomogeneities, residual selection bias, etc.) as discussed in full in Appendix A. In practice, the plausible range for this mean bias (1 − b) was estimated by comparing the observed relation with predictions from several sets of numerical simulations, as detailed in Appendix A.
The large uncertainties on (1 − b) are due to the dispersion in predictions from the various simulation sets. This is a major factor limiting the accuracy of our analysis. A value (1 − b) = 0.8 could be considered as a best guess given available simulations, with no clear dependence on mass or redshift. From one cluster to the next the ratio of to the true mass is expected to be stochastic, contributing to the scatter in the Y_{500}–M_{500} relation given below. A conspiracy of all possible sources of bias (departure from HE, absolute instrument calibration, temperature inhomogeneities, residual selection bias) would seem necessary to lead to a significantly lower value of (1 − b). This apparently implausible possibility needs to be excluded through tests using other probes such as baryon and gas fractions, gas pressure, etc. As a baseline we take (1 − b) to vary within the range [ 0.7,1.0 ] with a flat prior. We also consider, when analysing systematic uncertainties on the derived cosmological parameters, a case where the bias is fixed to the value (1 − b) = 0.8.
As detailed in Appendix A, we derive a baseline relation for the mean SZ signal from a cluster of given mass and redshift in the form (7)where D_{A}(z) is the angulardiameter distance to redshift z and E^{2}(z) = Ω_{m}(1 + z)^{3} + Ω_{Λ}. The coefficients Y_{∗}, α, and β are given in Table 1.
Summary of scalinglaw parameters and error budget.
Equation (7) has an estimated intrinsic scatter^{3}σ_{log Y} = 0.075, which we take to be independent of redshift (see Appendix A). This is incorporated by drawing the cluster’s Y_{500} from a lognormal distribution (8)where is given by Eq. (7). Inclusion of this scatter increases the number of clusters expected at a given S/N; since the cluster counts are a steep function of M_{500} in the range of mass in question, there are more clusters that scatter upwards from below the limit given by the zeroscatter scaling relation than those that scatter downwards.
In addition to Eq. (7) we need a relation between θ_{500} (in fact , the angular size corresponding to the physical size ), the aperture used to extract Y_{500}, and M_{500}. Since and θ_{500} = R_{500}/D_{A}, this can be expressed as (9)where θ_{∗} = 6.997 arcmin.
2.4. Limiting mass
One can use Eqs. (7) and (9) to compute the limiting mass at a point on the sky where the noise level, σ_{Y}, has been computed as described in Sect. 3. As the latter is not homogeneous on the sky, we show in Fig. 3 the limiting mass, defined at 50% completeness, as a function of redshift for three different zones, deep, medium, and shallow, covering respectively, 3.5%, 47.8%, and 48.7% of the unmasked sky. For each line a S/N cut of 7 has been adopted.
Fig. 3 Limiting mass as a function of z for the selection function and noise level computed for three zones (deep, blue; medium, cyan; shallow, green), and on average for the unmasked sky (dashed black). 
2.5. Implementation
We have implemented three independent versions of the computation of counts and constraints. The differences in predicted counts are of the order of a few percent, which translates to less than a tenth of 1σ on the cosmological parameters of interest.
3. The Planck cosmological samples
3.1. Sample definition
The reference cosmological sample is constructed from the PSZ Catalogue published in Planck Collaboration XXIX (2014) and made public with the first release of Planck cosmological products. It is based on the SZ detections performed with the matched multifilter (MMF) method MMF3 (Melin et al. 2006), which relies on the use of a filter of adjustable width θ_{500} chosen to maximize the S/N of the detection. In order to ensure a high purity and to maximize the number of redshifts, the cosmological sample was constructed by selecting the SZ detections above a S/N threshold of 7 outside Galactic and point source masks covering 35% of the sky, as discussed in Paper I. From the original PSZ, only the information on S/N (for the selection) and redshift are used.
This sample contains 189 candidates. All but one are confirmed bona fide clusters with measured redshifts, including 184 spectroscopic redshifts. Among these confirmed clusters 12 were confirmed with followup programmes conducted by the Planck collaboration (see Paper I for details). The remaining nonconfirmed cluster candidate is a highreliability class i candidate, meaning that its characterization as a cluster is supported by data in other wavebands (see Paper I for details). It is thus considered as a bona fide cluster. The distribution on the sky of this baseline cosmological sample is shown in Fig. 1.
In addition to our reference sample, we consider two other samples drawn from the PSZ for consistency checks. One is based on the detections from the second implementation of the MMF algorithm, MMF1, described in Paper I. It contains 188 clusters with S/N> 7 and no missing redshifts, with almost complete overlap with the baseline sample (187 clusters in common). The third sample considered in the present study is also based on MMF3 detections but with a higher S/N cut of S/N> 8. It allows us to test selection effects and to probe the consistency of the results as a function of the S/N cut. It contains 136 clusters, all with measured redshifts.
The selection function for each of these samples is constructed as described in the next section.
3.2. Completeness
The completeness of the reference cosmological sample is computed with two distinct and complementary approaches: a semianalytic approach based on the assumption of Gaussian uncertainties, and a computational approach based on Monte Carlo cluster injection into real sky maps.
The completeness χ can be evaluated analytically by setting the probability of the measured SZ flux, Y_{500}, to be Gaussian distributed with a standard deviation equal to the noise, σ_{Y500}(θ_{500},l,b), computed for each size θ_{500} of the MMF filter and at each position (l,b) on the sky: (10)where X = 7 is the S/N threshold and the error function is defined as usual by (11)χ_{erf}(Y_{500},θ_{500},l,b) thus lies in the range [ 0,1 ] and gives the probability for a cluster of flux Y_{500} and size θ_{500} at position (l,b) to be detected at S/N ≥ X.
Fig. 4 Noise map σ_{Y500}(θ_{500}) for θ_{500} = 6arcmin. The PSZ is limited by instrumental noise at high ( b  > 20deg) Galactic latitude (deeper at ecliptic poles) and foreground noise at low Galactic latitude. The scale of the map ranges from 0.5 to 2 times the mean noise of the map, which is ⟨ σ_{Y500}(6arcmin) ⟩ = 2.2 × 10^{4} arcmin^{2}. 
The noise estimate σ_{Y500}(θ_{500},l,b) is a byproduct of the detection algorithm and can be written in the form (see e.g., Melin et al. 2006) (12)with F_{θ500}(k) being a vector of dimension N_{freq} (the six highest Planck frequencies here) containing the beamconvolved cluster template scaled to the known SZ frequency dependence. The cluster template assumed is the nonstandard universal pressure profile from Arnaud et al. (2010). P(k,l,b) is the noise power spectrum (dimension N_{freq}×N_{freq}) directly estimated from the data at position (l,b). Figure 4 shows σ_{Y500}(θ_{500},l,b) for θ_{500} = 6 arcmin in a Mollweide projection with the Galactic mask used in the analysis applied. As expected, the noise at high Galactic latitude is lower than in the areas contaminated by diffuse Galactic emission. The ecliptic pole regions have the lowest noise level, reflecting the longer Planck integration time in these highredundancy areas.
The Monte Carlo (MC) completeness is calculated by injecting simulated clusters into real sky maps following the method presented in Paper I, with the modifications that the 65% Galactic dust mask and a S/N> 7 threshold are applied to match the cosmological sample definition. The Monte Carlo completeness encodes effects not probed by the erf approximation, including the variation of cluster pressure profiles around the fiducial pressure profile used in the MMF, spatiallyvarying and asymmetric effective beams, and the effects of correlated nonGaussian uncertainties in the estimation of (Y_{500},θ_{500}). As shown in Fig. 5, the erfbased formula for the completeness is a good approximation to the Monte Carlo completeness. The agreement is best for the typical sizes probed by Planck (5 to 10 arcmin), though the two determinations of the completeness start to deviate for small and large sizes, due to beam and profile effects, respectively. For simplicity, we chose the erf formulation as the baseline. The effect of using the Monte Carlo completeness instead is discussed in Sect. 5.2.
Fig. 5 Completeness averaged over the unmasked sky as a function of Y_{500} for two different filter sizes, θ_{500} = 6 and 15.3 arcmin. The dashed lines show the semianalytic approximation of Eq. (10), while the diamonds with errors show the completeness estimated by the MC injection technique. 
Bestfit cosmological parameters for various combinations of data and analysis methods.
4. Likelihood and Markov chain Monte Carlo
4.1. The likelihood
We now have all the information needed to predict the counts in redshift bins for our theoretical models. To obtain cosmological constraints with the PSZ sample presented in Sect. 3, we construct a likelihood function based on Poisson statistics (Cash 1979): (13)where is the probability of finding N_{i} clusters in each of N_{b} bins given an expected number of n_{i} in each bin in redshift. The likelihood includes bins that contain no observed clusters. As a baseline, we assume bins in redshift of Δz = 0.1 and we checked that our results are robust when changing the bin size between 0.05 and 0.2. The modelled expected number n_{i} depends on the bin range in redshift and on the cosmological parameters, as described in Sect. 2. It also depends on the scaling relations and the selection function of the observed sample. The parameters of the scaling relations between flux (or size) and mass and redshift are taken to be Gaussian distributed with central values and uncertainties stated in Table 1, and with the scatter in Y_{500} incorporated into the method via the lognormal distribution with width σ_{log Y}.
In the PSZ, the redshifts have been collected from different observations and from the literature. Individual uncertainties in redshift are thus spread between 0.001 and 0.1. Most of the clusters in the cosmological sample have spectroscopic redshifts (184 out of 189) and we checked that the uncertainties in redshift are not at all dominant in our error budget; they are thus neglected. The cluster without known redshift is incorporated by scaling the counts by a factor 189/188, i.e., by assuming its redshift is drawn from the distribution defined by the other 188 objects.
4.2. Markov chain Monte Carlo
In order to impose constraints on cosmological parameters from our sample(s) given our modelled expected number counts, we modified CosmoMC (Lewis & Bridle 2002) to include the likelihood described above. We mainly study constraints on the spatiallyflat ΛCDM model, varying Ω_{m}, σ_{8}, Ω_{b}, H_{0}, and n_{s}, but also adding in the total neutrino mass, ∑ m_{ν}, in Sect. 6. In each of the runs, the nuisance parameters (Y_{∗}, α, σ_{log Y}) follow Gaussian priors, with the characteristics detailed in Table 1, and are marginalized over. The bias (1 − b) follows a flat prior in the range [0.7,1]. The redshift evolution of the scaling, β, is fixed to its reference value unless stated otherwise.
Fig. 6 Planck SZ constraints (+BAO+BBN) on ΛCDM cosmological parameters in red. The black lines show the constraints upon substituting HST constraints on H_{0} for the BAO constrainsts. Contours are 68 and 95% confidence levels. 
4.3. External datasets
When probing the six parameters of the ΛCDM model, we combine the Planck clusters with the Big Bang nucleosynthesis (BBN) constraints from Steigman (2008), Ω_{b}h^{2} = 0.022 ± 0.002. We also use either the H_{0} determination from HST by Riess et al. (2011), H_{0} = (73.8 ± 2.4) km s^{1} Mpc^{1}, or baryon acoustic oscillation (BAO) data. In the latter case we adopt the combined likelihood of Hinshaw et al. (2013) and Planck Collaboration XVI (2014), which uses the radial BAO scales observed by 6dFGRS (Beutler et al. 2011), SDSSDR7rec and SDSSDR9rec (Padmanabhan et al. 2012; Anderson et al. 2012), and WiggleZ (Blake et al. 2012).
5. Constraints from Planck clusters: ΛCDM
5.1. Results for Ω_{m} and σ_{8}
Cluster counts in redshift for our Planck cosmological sample are not sensitive to all parameters of the ΛCDM model. We focus first on (Ω_{m},σ_{8}), assuming that n_{s} follows a Gaussian prior centred on the bestfit Planck CMB value^{4} (n_{s} = 0.9603 ± 0.0073). We combine our SZ counts likelihood with the BAO and BBN likelihoods discussed earlier. We also incorporate the uncertainties on scaling parameters in Table 1. Allowing the bias to range uniformly over the interval [ 0.7,1.0 ], we find the expected degeneracy between the two parameters, σ_{8}(Ω_{m}/ 0.27)^{0.3} = 0.764 ± 0.025,^{5} with central values and uncertainties of Ω_{m} = 0.29 ± 0.02 and σ_{8} = 0.75 ± 0.03 (Table 2 and Fig. 6, red contours). The cluster counts as a function of redshift for the bestfit model are plotted in Fig. 7. When fixing the bias to (1 − b) = 0.8, the constraint on Ω_{m} remains unchanged while the constraint on σ_{8} becomes stronger: σ_{8}(Ω_{m}/ 0.27)^{0.3} = 0.78 ± 0.01 and σ_{8} = 0.77 ± 0.02 (Fig. 8).
To investigate how robust our results are when changing our priors, we repeat the analysis substituting the HST constraints on H_{0} for the BAO results. Figure 6 (black contours) shows that the main effect is to change the bestfit value of H_{0}, leaving the (Ω_{m},σ_{8}) degeneracy almost unchanged. In the following robustness tests, we assume a fixed mass bias, (1 − b) = 0.8, to better identify the effect of each of our assumptions.
5.2. Robustness to the observational sample
To test the robustness of our results, we performed the same analysis with different subsamples drawn from our cosmological sample or from the PSZ, as described in Sect. 3, following that section’s discussion of completeness. Figure 9 shows the likelihood contours of the three samples (blue, MMF3 S/N> 8; red, MMF3 S/N> 7; black, MMF1 S/N> 7) in the (Ω_{m},σ_{8}) plane. There is good agreement between the three samples. Obviously the three samples are not independent, as many clusters are common, but the noise estimates for MMF3 and MMF1 are different leading to different selection functions. Table 2 summarizes the bestfit values.
Fig. 7 Distribution in redshift for the Planck cosmological cluster sample. The observed number counts (red), are compared to our bestfit model prediction (blue). The dashed and dotdashed histograms are the bestfit models from the Planck SZ power spectrum and Planck CMB power spectrum fits, respectively. The cyan long dashed histogram is the best fit CMB and SZ when the bias is free (see Sect. 6.3). The uncertainties on the observed counts, shown for illustration only, are the standard deviation based on the observed counts, except for empty bins where we show the inferred 84% upper limit on the predicted counts assuming a Poissonian distribution. See Sect. 6 for more discussion. 
Fig. 8 Comparison of the constraints using the mass functions of Watson et al. (black) and Tinker et al. (red), with b fixed to 0.8. When relaxing the constraints on the evolution of the scaling law with redshift (blue), the contours move along the degeneracy line. Allowing the bias to vary uniformly in the range [ 0.7,1.0 ] enlarges the constraints perpendicular to the σ_{8}–Ω_{m} degeneracy line due to the degeneracy of the number of clusters with the mass bias (purple). Contours are 95% confidence levels here. 
Fig. 9 95% contours for different robustness tests: MMF3 with S/N cut >7 in red; MMF3 with S/N cut >8 in blue; and MMF1 with S/N cut >7 in black; and MMF3 with S/N cut at 7 but adopting the MC completeness in purple. 
We perform the same analysis as on the baseline cosmological sample (SZ and BAO and BBN), but employing a different computation of the completeness function using the Monte Carlo method described in Sect. 3. Figure 9 shows the change in the 2D likelihoods when this alternative approach is adopted. The Monte Carlo estimation (in purple), being close to the analytic one, gives constraints that are similar, but enlarge the contour along the (Ω_{m},σ_{8}) degeneracy.
5.3. Robustness to cluster modelling
A key ingredient in the modelling of the number counts is the mass function. Our main results adopt the Tinker et al. mass function as the reference model. We compare these to results from the Watson et al. mass function to evaluate the impact of uncertainty in predictions for the abundance of the most massive and extreme clusters. Figure 8 shows the 95% contours when adopting the different mass functions. The main effect is to change the orientation of the degeneracy between Ω_{m} and σ_{8}, moving the bestfit values by less than 1σ.
Constraints from clusters on σ_{8}(Ω_{m}/ 0.27)^{0.3}.
We also relax the assumption of standard evolution of the scalings with redshift by allowing β to vary with a Gaussian prior taken from Planck Collaboration X (2011), β = 0.66 ± 0.5. Once again, the contours move along the σ_{8}–Ω_{m} degeneracy direction (shown in blue in Fig. 8).
As shown in Appendix A, the estimation of the mass bias is not trivial and there is a large scatter amongst simulations. The red and purple contours compare the different constraints when fixing the mass bias to 0.8 and when allowing it to vary uniformly in the range [ 0.7,1.0 ] respectively. Modelling of the cluster observablemass relation is clearly the limiting factor in our analysis.
6. Discussion
Our main result is the constraint in the (Ω_{m},σ_{8}) plane for the standard ΛCDM model imposed by the SZ counts, which we have shown is robust to the details of our modelling. We now compare this result first to constraints from other cluster samples, and then to the constraints from the Planck analysis of the skymap of the Compton yparameter (Planck Collaboration XXI 2014) and of the primary CMB temperature anisotropies (Planck Collaboration XVI 2014).
6.1. Comparison with other cluster constraints
We restrict our comparison to some recent analyses exploiting a range of observational techniques to obtain cluster samples and mass calibrations.
Benson et al. (2013) used 18 galaxy clusters in the first 178 deg^{2} of the SPT survey to find σ_{8}(Ω_{m}/ 0.25)^{0.3} = 0.785 ± 0.037 for a spatiallyflat model. They break the degeneracy between σ_{8} and Ω_{m} by incorporating primary CMB constraints, deducing that σ_{8} = 0.795 ± 0.016 and Ω_{m} = 0.255 ± 0.016. In addition, they find that the dark energy equation of state is constrained to w = − 1.09 ± 0.36, using just their cluster sample along with the same HST and BBN constraints used here. Subsequently, Reichardt et al. (2013) reported a much larger cluster sample and used this to improve on the statistical uncertainties on the cosmological parameters (see Table 3). Hasselfield et al. (2013) use a sample of 15 high S/N clusters from ACT, in combination with primary CMB data, to find σ_{8} = 0.786 ± 0.013 and Ω_{m} = 0.250 ± 0.012 when assuming a scaling law derived from the universal pressure profile.
Strong constraints on cosmological parameters have been inferred from Xray and optical richness selected samples. Vikhlinin et al. (2009b) used a sample of 86 wellstudied Xray clusters, split into low and highredshift bins, to conclude that Ω_{Λ}> 0 with a significance about 5σ and that w = −1.14 ± 0.21. Rozo et al. (2010) used the approximately 10^{4} clusters in the Sloan Digital Sky Survey (SDSS) MaxBCG cluster sample, which are detected using a colour–magnitude technique and characterized by optical richness. They found that σ_{8}(Ω_{m}/ 0.25)^{0.41} = 0.832 ± 0.033. The fact that this uncertainty is similar to those quoted above for much smaller cluster samples signifies, once again, that cluster cosmology constraints are now limited by modelling, rather than statistical, uncertainties.
Fig. 10 Comparison of constraints (68% confidence interval) on σ_{8}(Ω_{m}/ 0.27)^{0.3} from different experiments of large–scale structure (LSS), clusters, and CMB. The solid line ACT point assumes the same universal pressure profile as this work. Probes marked with an asterisk have an original power of Ω_{m} different from 0.3. See text and Table 3 for more details. 
Table 3 and Fig. 10 show some current constraints on the combination σ_{8}(Ω_{m}/ 0.27)^{0.3}, which is the main degeneracy line in cluster constraints. This comparison is only meant to be indicative, as a more quantitative comparison would require full consideration of modelling details which is beyond the scope of this work. Cosmic shear (Kilbinger et al. 2013), Xrays (Vikhlinin et al. 2009b), and MaxBCG (Rozo et al. 2010) each have a different slope in Ω_{m}, being 0.6, 0.47, and 0.41, respectively (instead of 0.3), as they are probing different redshift ranges. We have rescaled when necessary the best value and errors to quote numbers with a pivot Ω_{m} = 0.27. Hasselfield et al. (2013) have derived “clusteronly” constraints from ACT by adopting several different scaling laws, shown in blue and dashed blue in Fig. 10. The constraint assuming the universal pressure profile is highlighted as the solid symbol and error bar. For SPT we show the “clusteronly” constraints from Reichardt et al. (2013). For our own analysis we show our baseline result for SZ and BAO and BBN with a prior on (1 − b) distributed uniformly in [ 0.7,1 ]. The figure thus demonstrates good agreement amongst all cluster observations, whether in optical, Xrays, or SZ. Table 3 compares the different data and assumptions of the different clusterrelated publications.
6.2. Consistency with the Planck ymap
In a companion paper (Planck Collaboration XXI 2014), we performed an analysis of the SZ angular power spectrum derived from the Planckymap obtained with a dedicated componentseparation technique. For the first time, the power spectrum has been measured at intermediate scales (50 ≤ ℓ ≤ 1000). The same modelling as in Sect. 2 and Taburet et al. (2009, 2010) has been used to derive bestfit values of Ω_{m} and σ_{8}, assuming the universal pressure profile (Arnaud et al. 2010), a bias 1 − b = 0.8, and the bestfit values for other cosmological parameters from Planck Collaboration XVI (2014)^{6}. The best model obtained, shown in Fig. 7 as the dashed line, demonstrates the consistency between the PSZ number counts and the signal observed in the ymap.
6.3. Comparison with Planck primary CMB constraints
We now compare the PSZ cluster constraints to those from the analysis of the primary CMB temperature anisotropies given in Planck Collaboration XVI (2014) (see Footnote 6). In that analysis σ_{8} is derived from the standard six ΛCDM parameters.
The Planck primary CMB constraints, in the (Ω_{m},σ_{8}) plane, differ significantly from our own, in particular through favouring a higher value of σ_{8}, (see Fig. 11). For (1 − b) = 0.8, this leads to a factor of 2 larger number of predicted clusters than is actually observed (see Fig. 7). There is therefore some tension between the results from the Planck CMB analysis and the current cluster analysis. Figure 10 illustrates this with a comparison of three analyses of primary CMB data alone (Planck Collaboration XVI 2014; Story et al. 2013; Hinshaw et al. 2013) and cluster constraints in terms of σ_{8}(Ω_{m}/ 0.27)^{0.3}.
Fig. 11 2D Ω_{m}–σ_{8} likelihood contours for the analysis with Planck CMB only (red); Planck SZ and BAO and BBN (blue) with (1 − b) in [ 0.7,1 ]. 
It is possible that the tension results from a combination of some residual systematics with a substantial statistical fluctuation. Enough tests and comparisons have been made on the Planck data sets that it is plausible that at least one discrepancy at the two or three sigma level will arise by chance. Nevertheless, it is worth considering the implications if the discrepancy is real.
As we have noted, the modelling of the cluster gas physics is the most important uncertainty in our analysis, in particular through its influence on the mass bias (1 − b). While we have argued for a preferred value of (1 − b) ≃ 0.8 based on comparison of our Y_{500}–M_{500} relation to those derived from a number of different numerical simulations, and we suggest a plausible range of (1 − b) from 0.7 to 1, a significantly lower value would substantially alleviate the tension between CMB and SZ constraints. We have undertaken a joint analysis using the CMB likelihood presented in Planck Collaboration XV (2014) and the cluster likelihood presented in the present paper, sampling (1 − b) in the range [ 0.1,1.5 ]. This results in a “measurement” of (1 − b) = 0.59 ± 0.05. We show in Fig 7 the SZ cluster counts predicted by the Planck’s bestfit primary CMB model for (1 − b) = 0.59. Clearly, this substantial reduction in (1 − b) is enough to reconcile our observed SZ cluster counts with Planck’s bestfit primary CMB parameters.
Such a large bias is difficult to reconcile with numerical simulations, and cluster masses estimated from Xrays and from weak lensing do not typically show such large offsets (see Appendix A). Systematic discrepancies in the relevant scaling relations have, however, been identified and studied in stacking analyses of Xray, SZ, and lensing data for the very large MaxBCG cluster sample, e.g., Planck Collaboration XII (2011), Biesiadzinski et al. (2012), Draper et al. (2012), Rozo et al. (2012), and Sehgal et al. (2013), suggesting that the issue is not yet fully settled from an observational point of view. The uncertainty reflects the inherent biases of the different mass estimates. Systematic effects arising from instrument calibration constitute a further source of uncertainty – in Xray mass determinations, temperature estimates represent the main source of systematic uncertainty in mass, as the mass scales roughly with T^{3/2}. Other biases in the determination of massobservable scaling relations come from the object selection process itself (e.g., Mantz et al. 2010; Allen et al. 2011b; Angulo et al. 2012). This may be less of a concern for SZ selected samples because of the expected small scatter between the measured quantity and the mass. An improbable conspiracy of all sources of bias seems required to lead to a sufficientlylow effective value of (1 − b) to reconcile the SZ and CMB constraints. This possibility needs to be carefully examined with probes based on a variety of physical quantities and derived from a wide range of types of observation, including masses, baryon and gas fractions, etc.
A different mass function may also help reduce the tension. Mass functions are calibrated against numerical simulations that may still suffer from volume effects for the largest haloes, as shown in the difference between the Tinker et al. (2008) and Watson et al. (2013) mass functions. This does not seem sufficient, however, given the results presented in Fig. 8.
One might instead ask whether the Planck data analysis could somehow have missed a nonnegligible fraction of the total number of clusters currently predicted to have S/N> 7, resulting in a lower observed number count distribution. This is linked to a possible underestimate of the true dispersion about the Y_{500}–Y_{X} relation at a given M_{500}. It would be necessary for Planck to have missed ~40 percent of the clusters with predicted SZ S/N> 7 in order for the SZ and CMB number count curves in Fig. 7 to be in agreement. Increasing the dispersion about the Y_{500}–M_{500} relation and allowing it to correlate strongly with the scatter in Xray properties (in particular, Y_{X}) would raise the possibility that our calibration procedures (which are based on Xray and SZ selected clusters assuming the scatter in Y and Y_{x} at fixed M_{500} to be small and to be uncorrelated with cluster dynamical state) might produce a relation which is biased high. A sufficientlylarge effect seems, however, to require a level of scatter and a degree of correlation with cluster structure which are inconsistent with the predictions of current hydrodynamical simulations (see the discussion in Appendix A).
Alternatively, the discrepancy may reflect a need to extend the minimal ΛCDM model in which the σ_{8} constraints are derived from the primary CMB analysis. Any extension would need to modify the power spectrum on the scales probed by clusters, while leaving the scales probed by primary CMB observations unaffected. The inclusion of neutrino masses, quantified by their sum over all families, ∑ m_{ν}, can achieve this (see Marulli et al. 2011 and Burenin 2013 for reviews of how cosmological observations can be affected by the inclusion of neutrino masses). The SPT collaboration (Hou et al. 2014) recently considered such a possibility to mitigate their tension with WMAP7 primary CMB data. There is an upper limit of ∑ m_{ν} < 0.93 eV from the Planck primary CMB data alone^{7}. If we combine the Planck CMB (Planck+WP) likelihood and the cluster count data using a fixed value (1 − b) = 0.8, then we find a 2.8σ preference for the inclusion of neutrino masses with ∑ m_{ν} = (0.53 ± 0.19) eV, as shown in Fig. 12. If, on the other hand, we adopt a more conservative point of view and allow (1 − b) to vary between 0.7 and 1.0, this preference drops to 1.9σ with ∑ m_{ν} = (0.40 ± 0.21) eV. Adding BAO data to the compilation lowers the value of the required mass but increases the significance, yielding ∑ m_{ν} = (0.20 ± 0.09) eV, due to a breaking of the degeneracy between H_{0} and ∑ m_{ν}.
Fig. 12 Marginalized posterior distribution for ∑ m_{ν} from: Planck CMB data alone (black dotted line); Planck CMB and SZ with 1 − b in [ 0.7,1 ] (red); Planck CMB and SZ and BAO with 1 − b in [ 0.7,1 ] (blue); and Planck CMB and SZ with 1 − b = 0.8 (green). 
As these results depend on the value and allowed range of (1 − b), better understanding of the scaling relation is the key to further investigation. This provides strong motivation for further study of the relationship between Y and M. Over the past few years we have moved into an era where systematic uncertainties dominate to an increasing extent over statistical uncertainties. Observed mass estimates using different methods are of improving quality; for instance, Xray mass proxies can be measured to better than 10% precision (e.g., Vikhlinin et al. 2006; Arnaud et al. 2007). In this context, systematic calibration uncertainties are playing an increasingly important role when using the cluster population to constrain cosmology.
7. Summary
We have used a sample of nearly 200 clusters from the PSZ, along with the corresponding selection function, to place strong constraints in the (Ω_{m},σ_{8}) plane. We have carried out a series of tests to verify the robustness of our constraints, varying the observed sample choice, the estimation method for the selection function, and the theoretical methodology, and have found that our results are not altered significantly by those changes.
The relation between the mass and the integrated SZ signal plays a major role in the computation of the expected number counts. Uncertainties in cosmological constraints from clusters are no longer dominated by small number statistics, but by the gas physics and sample selection biases. Uncertainties in the Y–M relation include contributions from Xray instrument calibration, Xray temperature measurements, inhomogeneities in cluster density and temperature profiles, and selection effects. Considering several ingredients of the gas physics of clusters, numerical simulations predict scaling relations with 30% scatter in amplitude (at a fiducial mass of 6 × 10^{14}M_{⊙}), and suggest a mass bias between the true mass and the estimated mass of . Adopting the central value we found constraints on Ω_{m} and σ_{8} that are in good agreement with previous measurements using clusters of galaxies.
Comparing our results with Planck primary CMB constraints within the ΛCDM cosmology reveals some tension. This can be alleviated by permitting a large mass bias (1 − b ≃ 0.60), which is however significantly larger than expected. Alternatively, the tension may indicate a need for an extension of the base ΛCDM model that modifies its power spectrum shape. For example the inclusion of nonzero neutrino masses helps in reconciling the primary CMB and cluster constraints, a fit to Planck CMB and SZ and BAO yielding ∑ m_{ν} = (0.20 ± 0.09) eV.
Cosmological parameter determination using clusters is currently limited by the knowledge of the observable–mass relations. In the future our goal is to increase the number of dedicated followup programmes to obtain better estimates of the mass proxy and redshift for most of the S/N> 5Planck clusters. This will allow for improved determination of the scaling laws and the mass bias, increase the number of clusters that can be used, and allow us to investigate an extended cosmological parameter space.
Planck (http://www.esa.int/Planck) is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states (in particular the lead countries France and Italy), with contributions from NASA (USA) and telescope reflectors provided by a collaboration between ESA and a scientific consortium led and funded by Denmark.
Table 2 of Planck Collaboration XVI (2014).
For Planck CMB we took the constraints from the Planck+WP case, Col. 6 of Table 2 of Planck Collaboration XVI (2014). The baseline model includes massive neutrinos with ∑ m_{ν} = 0.06 eV.
Planck Collaboration XVI (2014), Table 10, Col. 3.
Acknowledgments
The development of Planck has been supported by: ESA; CNES and CNRS/INSUIN2P3INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN, JA and RES (Spain); Tekes, AoF and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); and PRACE (EU). A description of the Planck Collaboration and a list of its members, including the technical or scientific activities in which they have been involved, can be found at http://www.sciops.esa.int/index.php?project=planck&page=Planck_Collaboration
References
 Abell, G. O. 1958, ApJS, 3, 211 [NASA ADS] [CrossRef] [Google Scholar]
 Abell, G. O., Corwin, Jr., H. G., & Olowin, R. P. 1989, ApJS, 70, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Aghanim, N., de Luca, A., Bouchet, F. R., Gispert, R., & Puget, J. L. 1997, A&A, 325, 9 [NASA ADS] [Google Scholar]
 Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 [NASA ADS] [CrossRef] [Google Scholar]
 Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011a, Ann. Rev. Astron. Astrophys., 49, 409 [Google Scholar]
 Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011b, ARA&A, 49, 409 [NASA ADS] [CrossRef] [Google Scholar]
 Anderson, L., Aubourg, E., Bailey, S., et al. 2012, MNRAS, 427, 3435 [Google Scholar]
 Angulo, R. E., Springel, V., White, S. D. M., et al. 2012, MNRAS, 426, 2046 [CrossRef] [Google Scholar]
 Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 474, L37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Barbosa, D., Bartlett, J. G., Blanchard, A., & Oukbir, J. 1996, A&A, 314, 13 [NASA ADS] [Google Scholar]
 Battaglia, N., Bond, J. R., Pfrommer, C., & Sievers, J. L. 2012, ApJ, 758, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Battye, R. A., & Weller, J. 2003, Phys. Rev. D, 68, 083506 [NASA ADS] [CrossRef] [Google Scholar]
 Becker, M. R., & Kravtsov, A. V. 2011, ApJ, 740, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Benson, B. A., de Haan, T., Dudley, J. P., et al. 2013, ApJ, 763, 147 [NASA ADS] [CrossRef] [Google Scholar]
 Beutler, F., Blake, C., Colless, M., et al. 2011, MNRAS, 416, 3017 [NASA ADS] [CrossRef] [Google Scholar]
 Biesiadzinski, T., McMahon, J., Miller, C. J., Nord, B., & Shaw, L. 2012, ApJ, 757, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Birkinshaw, M. 1999, Phys. Rep., 310, 97 [NASA ADS] [CrossRef] [Google Scholar]
 Blake, C., Brough, S., Colless, M., et al. 2012, MNRAS, 425, 405 [NASA ADS] [CrossRef] [Google Scholar]
 Böhringer, H., Voges, W., Huchra, J. P., et al. 2000, ApJS, 129, 435 [NASA ADS] [CrossRef] [Google Scholar]
 Böhringer, H., Schuecker, P., Guzzo, L., et al. 2004, A&A, 425, 367 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Burenin, R. A. 2013, Astron. Lett., 39, 357 [Google Scholar]
 Carlberg, R. G., Yee, H. K. C., Ellingson, E., et al. 1996, ApJ, 462, 32 [Google Scholar]
 Carlstrom, J. E., Ade, P. A. R., Aird, K. A., et al. 2011, PASP, 123, 568 [NASA ADS] [CrossRef] [Google Scholar]
 Cash, W. 1979, ApJ, 228, 939 [NASA ADS] [CrossRef] [Google Scholar]
 da Silva, A. C., Kay, S. T., Liddle, A. R., & Thomas, P. A. 2004, MNRAS, 348, 1401 [NASA ADS] [CrossRef] [Google Scholar]
 Diego, J. M., MartínezGonzález, E., Sanz, J. L., Benitez, N., & Silk, J. 2002, MNRAS, 331, 556 [NASA ADS] [CrossRef] [Google Scholar]
 Draper, P., Dodelson, S., Hao, J., & Rozo, E. 2012, Phys. Rev. D, 85, 023005 [NASA ADS] [CrossRef] [Google Scholar]
 Fabjan, D., Borgani, S., Rasia, E., et al. 2011, MNRAS, 416, 801 [NASA ADS] [CrossRef] [Google Scholar]
 Haiman, Z., Mohr, J. J., & Holder, G. P. 2001, ApJ, 553, 545 [NASA ADS] [CrossRef] [Google Scholar]
 Hasselfield, M., Hilton, M., Marriage, T. A., et al. 2013, J. Cosmol. Astropart. Phys., 7, 8 [Google Scholar]
 Henry, J. P., Evrard, A. E., Hoekstra, H., Babul, A., & Mahdavi, A. 2009, ApJ, 691, 1307 [NASA ADS] [CrossRef] [Google Scholar]
 Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Holder, G., Haiman, Z., & Mohr, J. J. 2001, ApJ, 560, L111 [NASA ADS] [CrossRef] [Google Scholar]
 Hou, Z., Reichardt, C. L., Story, K. T., et al. 2014, ApJ, 782, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Kay, S. T., Peel, M. W., Short, C. J., et al. 2012, MNRAS, 422, 1999 [NASA ADS] [CrossRef] [Google Scholar]
 Kilbinger, M., Fu, L., Heymans, C., et al. 2013, MNRAS, 430, 2200 [NASA ADS] [CrossRef] [Google Scholar]
 Krause, E., Pierpaoli, E., Dolag, K., & Borgani, S. 2012, MNRAS, 419, 1766 [NASA ADS] [CrossRef] [Google Scholar]
 Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 [NASA ADS] [CrossRef] [Google Scholar]
 Lau, E. T., Kravtsov, A. V., & Nagai, D. 2009, ApJ, 705, 1129 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, A., & Bridle, S. 2002, Phys. Rev., D66, 103511 [Google Scholar]
 Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 473 [NASA ADS] [CrossRef] [Google Scholar]
 Mahdavi, A., Hoekstra, H., Babul, A., & Henry, J. P. 2008, MNRAS, 384, 1567 [NASA ADS] [CrossRef] [Google Scholar]
 Mahdavi, A., Hoekstra, H., Babul, A., et al. 2013, ApJ, 767, 116 [NASA ADS] [CrossRef] [Google Scholar]
 Mantz, A., Allen, S. W., Ebeling, H., Rapetti, D., & DrlicaWagner, A. 2010, MNRAS, 406, 1773 [NASA ADS] [Google Scholar]
 Marrone, D. P., Smith, G. P., Okabe, N., et al. 2012, ApJ, 754, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Marulli, F., Carbone, C., Viel, M., Moscardini, L., & Cimatti, A. 2011, MNRAS, 418, 346 [NASA ADS] [CrossRef] [Google Scholar]
 Mazzotta, P., Rasia, E., Moscardini, L., & Tormen, G. 2004, MNRAS, 354, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Mehrtens, N., Romer, A. K., Hilton, M., et al. 2012, MNRAS, 423, 1024 [NASA ADS] [CrossRef] [Google Scholar]
 Melin, J., Bartlett, J. G., & Delabrouille, J. 2006, A&A, 459, 341 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Meneghetti, M., Rasia, E., Merten, J., et al. 2010, A&A, 514, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Motl, P. M., Hallman, E. J., Burns, J. O., & Norman, M. L. 2005, ApJ, 623, L63 [NASA ADS] [CrossRef] [Google Scholar]
 Nagai, D. 2006, ApJ, 650, 538 [NASA ADS] [CrossRef] [Google Scholar]
 Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007, ApJ, 668, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Oukbir, J., & Blanchard, A. 1992, A&A, 262, L21 [NASA ADS] [Google Scholar]
 Padmanabhan, N., Xu, X., Eisenstein, D. J., et al. 2012, MNRAS, 427, 2132 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Perrenod, S. C. 1980, ApJ, 236, 373 [NASA ADS] [CrossRef] [Google Scholar]
 Piffaretti, R. & Valdarnini, R. 2008, A&A, 491, 71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VIII. 2011, A&A, 536, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IX. 2011, A&A, 536, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration X. 2011, A&A, 536, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XI. 2011, A&A, 536, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XII. 2011, A&A, 536, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. I. 2012, A&A, 543, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. III. 2013, A&A, 550, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. IV. 2013, A&A, 550, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration I. 2014, A&A, 571, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration II. 2014, A&A, 571, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration III. 2014, A&A, 571, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IV. 2014, A&A, 571, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration V. 2014, A&A, 571, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VI. 2014, A&A, 571, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VII. 2014, A&A, 571, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VIII. 2014, A&A, 571, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IX. 2014, A&A, 571, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration X. 2014, A&A, 571, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XII. 2014, A&A, 571, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2014, A&A, 571, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIV. 2014, A&A, 571, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XV. 2014, A&A, 571, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVII. 2014, A&A, 571, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVIII. 2014, A&A, 571, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIX. 2014, A&A, 571, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XX. 2014, A&A, 571, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXI. 2014, A&A, 571, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXII. 2014, A&A, 571, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIII. 2014, A&A, 571, A23 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Planck Collaboration XXIV. 2014, A&A, 571, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXV. 2014, A&A, 571, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVI. 2014, A&A, 571, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVII. 2014, A&A, 571, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVIII. 2014, A&A, 571, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIX. 2014, A&A, 571, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXX. 2014, A&A, 571, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXXI. 2014, A&A, 571, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Poole, G. B., Babul, A., McCarthy, I. G., et al. 2007, MNRAS, 380, 437 [NASA ADS] [CrossRef] [Google Scholar]
 Pratt, G. W., Croston, J. H., Arnaud, M., & Böhringer, H. 2009, A&A, 498, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rasia, E., Meneghetti, M., Martino, R., et al. 2012, New J. Phys., 14, 055018 [NASA ADS] [CrossRef] [Google Scholar]
 Reichardt, C. L., Stalder, B., Bleem, L. E., et al. 2013, ApJ, 763, 127 [NASA ADS] [CrossRef] [Google Scholar]
 Ricker, P. M., & Sarazin, C. L. 2001, ApJ, 561, 621 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Macri, L., Casertano, S., et al. 2011, ApJ, 730, 119 [NASA ADS] [CrossRef] [Google Scholar]
 Ritchie, B. W., & Thomas, P. A. 2002, MNRAS, 329, 675 [NASA ADS] [CrossRef] [Google Scholar]
 Rozo, E., Wechsler, R. H., Rykoff, E. S., et al. 2010, ApJ, 708, 645 [NASA ADS] [CrossRef] [Google Scholar]
 Rozo, E., Vikhlinin, A., & More, S. 2012, ApJ, 760, 67 [NASA ADS] [CrossRef] [Google Scholar]
 Sehgal, N., Bode, P., Das, S., et al. 2010, ApJ, 709, 920 [NASA ADS] [CrossRef] [Google Scholar]
 Sehgal, N., Addison, G., Battaglia, N., et al. 2013, ApJ, 767, 38 [NASA ADS] [CrossRef] [Google Scholar]
 Sembolini, F., Yepes, G., De Petris, M., et al. 2013, MNRAS, 434, 2718 [NASA ADS] [CrossRef] [Google Scholar]
 Springel, V., White, M., & Hernquist, L. 2001, ApJ, 562, 1086 [NASA ADS] [CrossRef] [Google Scholar]
 Steigman, G. 2008 [arXiv:0807.3004] [Google Scholar]
 Story, K. T., Reichardt, C. L., Hou, Z., et al. 2013, ApJ, 779, 86 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1970, Ap&SS, 7, 20 [Google Scholar]
 Swetz, D. S., Ade, P. A. R., Amiri, M., et al. 2011, ApJS, 194, 41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Taburet, N., Aghanim, N., Douspis, M., & Langer, M. 2009, MNRAS, 392, 1153 [NASA ADS] [CrossRef] [Google Scholar]
 Taburet, N., Douspis, M., & Aghanim, N. 2010, MNRAS, 404, 1197 [NASA ADS] [Google Scholar]
 Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709 [NASA ADS] [CrossRef] [Google Scholar]
 Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 691 [NASA ADS] [CrossRef] [Google Scholar]
 Vikhlinin, A., Burenin, R. A., Ebeling, H., et al. 2009a, ApJ, 692, 1033 [NASA ADS] [CrossRef] [Google Scholar]
 Vikhlinin, A., Kravtsov, A. V., Burenin, R. A., et al. 2009b, ApJ, 692, 1060 [NASA ADS] [CrossRef] [Google Scholar]
 Voit, G. M. 2005, Rev. Mod. Phys., 77, 207 [NASA ADS] [CrossRef] [Google Scholar]
 von der Linden, A., Allen, M. T., Applegate, D. E., et al. 2014, MNRAS, 439, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Watson, W. A., Iliev, I. T., D’Aloisio, A., et al. 2013, MNRAS, 433, 1230 [Google Scholar]
 Weller, J., Battye, R. A., & Kneissl, R. 2002, Phys. Rev. Lett., 88, 231301 [NASA ADS] [CrossRef] [Google Scholar]
 White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 1993, Nature, 366, 429 [NASA ADS] [CrossRef] [Google Scholar]
 Willis, J. P., Clerc, N., Bremer, M. N., et al. 2013, MNRAS, 430, 134 [NASA ADS] [CrossRef] [Google Scholar]
 Yang, H.Y. K., Bhattacharya, S., & Ricker, P. M. 2010, ApJ, 725, 1124 [NASA ADS] [CrossRef] [Google Scholar]
 Zeldovich, Y. B., & Sunyaev, R. A. 1969, Ap&SS, 4, 301 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, Y.Y., Okabe, N., Finoguenov, A., et al. 2010, ApJ, 711, 1033 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Calibration of the Y_{500}–M_{500} relation
A cluster catalogue is a list of positions and measurements of observable physical quantities. Its scientific utility depends largely on our ability to link the observed quantities to the underlying mass, in other words, to define an observable proxy for the mass. Planck detects clusters through the SZ effect. This effect is currently the subject of much study in the cluster community, chiefly because numerical simulations indicate that the sphericallyintegrated SZ measurement is correlated particularly tightly with the underlying mass. In other words, this measurement potentially represents a particularly valuable mass proxy.
To establish a mass proxy, one obviously needs an accurate measurement both of the total mass and of the observable quantity in question. However, even with highly accurate measurements, the correlation between the observable quantity and the mass is susceptible to bias and dispersion, and both of these effects need to be taken into account when using cluster catalogues for scientific applications.
The aim of this Appendix is to define a baseline relation between the measured SZ flux, Y_{500}, and the total mass M_{500}. The latter quantity is not directly measurable. On an individual cluster basis, it can be inferred from dynamical analysis of galaxies, from Xray analysis assuming hydrostatic equilibrium (HE), or from gravitational lensing observations. However, it is important to note that all observed mass estimates include inherent biases. For instance, numerical simulations suggest that HE mass measurements are likely to underestimate the true mass by 10–15 percent due to neglect of bulk motions and turbulence in the intracluster medium (ICM, e.g., Nagai et al. 2007; Piffaretti & Valdarnini 2008; Meneghetti et al. 2010), an effect that is commonly referred to in the literature as the “hydrostatic mass bias”. Similarly, simulations indicate that weak lensing mass measurements may underestimate the mass by 5 to 10 percent, owing to projection effects or the use of inappropriate mass models (e.g., Becker & Kravtsov 2011). Instrument calibration systematic effects constitute a further source of error. For Xray mass determinations, temperature estimates represent the main source of systematic uncertainty, as the mass at a given density contrast scales roughly with T^{3/2}. Other biases in the determination of massobservable scaling relations come from the object selection process itself (e.g., Allen et al. 2011b; Angulo et al. 2012). A classic example is the Malmquist bias, where bright objects near the flux limit are preferentially detected. This effect is amplified by Eddington bias, the mass function dictating that many more lowmass objects are detected compared to highmass objects. Both of these biases depend critically on the distribution of objects in mass and redshift, and on the dispersion in the relation between the mass and the observable used for sample selection. This is less of a concern for SZ selected samples than for Xray selected samples, the SZ signal having much less scatter at a given mass than the Xray luminosity. However for precise studies it should still be taken into account.
On the theoretical side, numerous Y_{500}–M_{500} relations have been derived from simulated data, as discussed below. The obvious advantage of using simulated data is that the relation between the SZ signal and the true mass can be obtained, because the “real” value of all physical quantities can be measured. The disadvantage is that the “real” values of measurable physical quantities depend strongly on the phenomenological models used to describe the different nongravitational processes at work in the ICM.
Nevertheless, the magnitude of the bias between observed and true quantities can only be assessed by comparing multiwavelength observations of a wellcontrolled cluster sample to numerical simulations. Thus, ideally, we would have full followup of a complete Planck cluster sample. For large samples, however, full followup is costly and time consuming. This has led to the widespread use of mass estimates obtained from massproxy relations. These relations are generally calibrated from individual deep observations of a subset of the sample in question (e.g., Vikhlinin et al. 2009a), or from deep observations of objects from an external dataset (e.g., the use of the REXCESS relations in Planck Collaboration XI 2011).
For the present paper, we will rely on mass estimates from a mass–proxy relation. In this context, the M_{500}–Y_{X} relation is clearly the best to use. Y_{X}, proposed by Kravtsov et al. (2006), is defined as the product of M_{g,500}, the gas mass within R_{500}, and T_{X}, the spectroscopic temperature measured in the [ 0.15–0.75 ] R_{500} aperture. In the simulations performed by Kravtsov et al. (2006), Y_{X} was extremely tightly correlated with the true cluster mass, with a logarithmic dispersion of only 8 percent. Observations using masses derived from Xray hydrostatic analysis indicate that Y_{X} does indeed appear to have a low dispersion (Arnaud et al. 2007; Vikhlinin et al. 2009a). Furthermore, the local M_{500}–Y_{X} relation for Xray selected relaxed clusters has been calibrated to high statistical precision (Arnaud et al. 2010; Vikhlinin et al. 2009a), with excellent agreement achieved between various observations (see e.g., Arnaud et al. 2007). Since simulations suggest that the Y_{500}–M_{500} relation is independent of dynamical state, calibrating the Y_{500}–M_{500} relation via a lowscatter mass proxy, itself calibrated on clusters for which the HE bias is expected to be minimal, is a better approach than using HE mass estimates for the full sample, since the latter can be highly biased for very unrelaxed objects.
We approach the determination of the Y_{500}–M_{500} relation in two steps. We first calibrate the Y_{500}–proxy relation. This is combined with the Xray calibrated relation, between the proxy and M_{500}, to define an observationbased Y_{500}–M_{500} relation. In the second step, we assess possible biases in the relation by directly comparing the observationbased relation with that from simulations. This approach, rather than directly assessing the HE mass bias, allows us to avoid complications linked to the strong dependence of the HE bias on cluster dynamical state, and thus on the cluster sample (real or simulated). The final output from this procedure is a relation between Y_{500} and M_{500}, with a full accounting of the different statistical and systematic uncertainties that go into its derivation, including bias.
In the following, all relations are fit with a power law in logspace using the orthogonal BCES method (Akritas & Bershady 1996), which takes into account the uncertainties in both variables and the intrinsic scatter. All dispersions are given in log_{10}.
Appendix A.1: Baseline massproxy relation
As a baseline, we use the relation between Y_{X} and the Xray hydrostatic mass established for 20 local relaxed clusters by Arnaud et al. (2010): (A.1)assuming standard evolution, and where the uncertainties are statistical only. For easier comparison with the Y_{500}–M_{500} relation given below, the normalization for Y_{X} expressed in 10^{4} Mpc^{2} is 10^{− 0.171 ± 0.018}. The HE mass is expected to be a biased estimator of the true mass, (A.2)where all of the possible observational biases discussed above (departure from HE, absolute instrument calibration, temperature inhomogeneities, residual selection bias) have been subsumed into the bias factor (1 − b). The form of the Y_{X}–M_{500} relation is thus (A.3)where σ_{A} and σ_{α} are the statistical uncertainties on the normalization and slope and b is the bias between the true mass and the observed mass used to calibrate the relation. The bias is a poorlyknown stochastic variable with substantial variation expected between clusters. In our case, b represents the mean bias between the observed mass and the true mass.
Parameters for the Y_{500}–M_{500} relation, expressed as .
The mass proxy is defined from the bestfit Y_{X}– relation (A.4)For any cluster, , together with the corresponding Y_{X} and , can be estimated iteratively about this relation from the observed temperature and gas mass profile, as described in Kravtsov et al. (2006). The calibration of the Y_{X}–M_{500} relation is equivalent to a calibration of the –M_{500} relation, which relates the mass proxy, , to the mass via (A.5)In addition to the bias factor, there are statistical uncertainties on the slope and normalization of the relation, as well as intrinsic scatter around the relation, linked to the corresponding statistical uncertainties and scatter of the Y_{X}– relation.
Appendix A.2: Relation between Y_{500} and M
Appendix A.2.1: Bestfit relation
We first investigate the relationship between Y_{500} and , the mass estimated iteratively from Eq. A.4, with parameters given by the bestfit Arnaud et al. (2010) relation (Eq. (A.1)). Full Xray followup of the Planck SZ cosmological cluster sample is not yet available. Our baseline sample is thus a subset of 71 detections from the Planck cosmological cluster sample, detected at S/N> 7, for which good quality XMMNewton observations are available. The sample consists of data from our previous archival study of the Planck Early SZ (ESZ) clusters (Planck Collaboration XI 2011), of Planckdetected LoCuSS clusters presented by Planck Collaboration Int. III (2013), and from the XMMNewton validation programme (Planck Collaboration IX 2011; Planck Collaboration Int. I 2012; Planck Collaboration Int. IV 2013). The corresponding subsamples include 58, 4, and 9 clusters, respectively. The Xray data were reanalysed in order to have a homogeneous data set; measurement differences are negligible with respect to previouslypublished values. Uncertainties on Y_{X}, , and include those due to statistical errors on the Xray temperature and the gas mass profile.
The SZ signal is estimated within a sphere of radius centred on the position of the Xray peak, as detailed in e.g., Planck Collaboration XI (2011). The reextraction procedure uses matched multifilters (MMF) and assumes that the ICM pressure follows the universal profile shape derived by Arnaud et al. (2010) from the combination of the REXCESS sample with simulations. The extraction is undertaken on the 15.5 month Planck survey data set, and so statistical precision on the SZ signal is improved with respect to previouslypublished values. The uncertainty on Y_{500} includes statistical uncertainties on the SZ signal derived from the MMF, plus the statistical uncertainty on the aperture . The latter uncertainty is negligible compared to the statistical error on the SZ signal. The resulting relation for these 71 clusters from the cosmological sample is (A.6)This agrees within 1σ with the results from the sample of 62 clusters from the ESZ sample with archival XMMNewton data published in Planck Collaboration XI (2011). The slope and normalization are determined at slightly higher precision, due to the better quality SZ data. The derived intrinsic scatter (Table A.1) is significantly smaller. This is a consequence of: a more robust treatment of statistical uncertainties; propagation of gas mass profile uncertainties in the Y_{X} error budget; and, to a lesser extent, the propagation of uncertainties to Y_{500} estimates.
Appendix A.2.2: Effects of Malmquist bias
The fitted parameters are potentially subject to selection effects such as Malmquist bias, owing to part of the sample lying close to the selection cut. For the present sample, we use an approach adapted from that described in Vikhlinin et al. (2009a) and Pratt et al. (2009), where each data point is rescaled by the mean bias for its flux, and the relation refitted using the rescaled points. The method is described in more detail in Paper I. For the baseline cosmological sample of 71 systems, the biascorrected Y_{500}– relation is (A.7)The bestfit relation, together with Malmquist bias corrected data points, is plotted in Fig. A.1.
Fig. A.1 Best scaling relation between Y_{500} and M_{500}, and the data points utilized after correction of the Malmquist bias. 
The correction decreases the effective Y_{500} values at a given mass, an effect larger for clusters closer to the S/N threshold. The net effect is small, a roughly 1σ decrease of the normalization and a slight steepening of the powerlaw slope (Table A.1).
Appendix A.2.3: Stability of slope and normalization
The slope and normalization of this relation are robust to the choice of sample (Table A.1). We compared our results to those obtained from:

An extended sample of 78 clusters withS/N> 7 (71 in common with the baseline sample). This is built from all objects falling in the 84% sky mask used to define the SZ catalogue Planck Collaboration XXIX (2014), and for which XMMNewton data have been published by the Planck Collaboration (Planck Collaboration IX 2011; Planck Collaboration XI 2011; Planck Collaboration Int. I 2012; Planck Collaboration Int. III 2013; Planck Collaboration Int. IV 2013).

The original 62 clusters from the ESZ sample published in Planck Collaboration XI (2011), with updated SZ signal measurements obtained from 15.5 month Planck data (62 in common with the baseline sample). These objects are all known from Xray surveys and all lie at z < 0.5. We use them to test fit robustness to the inclusion of nonXray selected, higherredshift systems.
As indicated in Table A.1, there is agreement within 1 σ between the various samples. The results are also in agreement with the relation obtained from a simple combination of the Y_{500}–Y_{X} relation (discussed in Paper I) and the Y_{X}– relation (Eq. (A.1) above).
Appendix A.3: The observationbased Y_{500}–M_{500} relation
Appendix A.3.1: Combination of the Y_{500}–M and the M–M_{500} relations
We now combine Eq. (A.7) with the –M_{500} relation. This will not change the bestfit parameters, but will increase their uncertainties. As the determinations of the two relations are independent, we added quadratically the uncertainties in the bestfit parameters of the Y_{500}– (Eq. (A.6)) and –M_{500} (Eq. (A.5), with values from Eq. (A.1)) relations. Our bestfit Y_{500}–M_{500} relation is then (A.8)Thus inclusion of the statistical uncertainty in the – relation doubles the uncertainty on the normalization and increases the uncertainty on the slope by 40%. Note that we have implicitly assumed here that the scatter around the two relations is uncorrelated.
Appendix A.3.2: Effect of use of an external dataset
The above results assume a mass estimated from the baseline Y_{X}–M_{500} relation, derived by Arnaud et al. (2010) from an external dataset of 20 relaxed clusters (Eq. (A.1)). How does this relation compare to the individual hydrostatic Xray masses of the Planck cosmological cluster sample? Of the 71 clusters in the baseline sample:

58 objects have temperature profile information extending tovarious fractions of R_{500}, of which

10 coolcore objects have temperature profiles measurements at least out to R_{500}.
Thus, while spatiallyresolved temperature profiles are available for 58 of the 71 clusters with XMMNewton observations, we must be careful in interpretation of these data. The Arnaud et al. relation was derived from a carefully chosen data set consisting of relaxed, coolcore objects having wellconstrained temperature profiles out to around R_{500}, i.e., the type of object for which it makes sense to undertake a hydrostatic mass analysis. Many clusters of the Planck sample are merging systems for which such an analysis would give results that are difficult to interpret. In addition, few of the Planck sample have spatiallyresolved temperature profiles out to R_{500}. However, as given in Table A.1, the bestfit Y_{X}– relation for the 10 coolcore clusters that are detected to R_{500} agrees with Eq. (A.8) within 1σ. Moreover, the relation for the 58 Planck clusters with HE mass estimates, derived regardless of dynamical state and radial detection extent, also agrees within 1σ (albeit with greatly increased scatter). We are thus confident that the masses estimated from an externallycalibrated Y_{X}– relation are applicable to the present data set.
Fig. A.2 Left: comparison of Y_{500}–M_{500} relations from 12 simulations undertaken by six different groups with the updated observational Y_{500}– result from Planck, Eq. (A.8). Right: ratio of each simulated Y_{500}–M_{500} relation relative to Eq. (A.8). The different scaling laws are taken from Kay et al. (2012), Battaglia et al. (2012), Yang et al. (2010), Sehgal et al. (2010), Krause et al. (2012), Nagai (2006), Sembolini et al. (2013) and Planck Collaboration XI (2011). 
Appendix A.3.3: Dispersion about the observed relations
A key issue is the dispersion around the mean relation. We first estimate the intrinsic scatter of the Y_{500}– relation by combining the intrinsic scatter of the Y_{500}– relation and that of the – relation. This estimate is applicable to relaxed objects only, since the Y_{500}– relation has been measured using a sample of such systems. If the scatter about the input relations is independent, this gives (A.9)where α is the slope of the Y_{500}– relation. As the HE mass estimate introduces extra scatter as compared to the true mass (Kay et al. 2012), the dispersion about the Y_{500}–M_{500} relation is expected to be smaller than that of the Y_{500}– relation (although this will depend on correlations between the scatter in the –M_{500} and –Y_{500} relations). The above expression thus also provides an estimate of the scatter of the Y_{500}–M_{500} relation, again for relaxed objects. While merging events are expected to induce shocks in the ICM, leading to higher temperatures and thus an increase in Y_{500}, current simulations suggest that this is a weak effect. This may be due to the relatively short duration of the shocking phase during a merger (e.g., Ricker & Sarazin 2001; Ritchie & Thomas 2002; Poole et al. 2007). Thus, further assuming that the intrinsic scatter of the Y_{500}–M_{500} relation is the same for the total relaxed and unrelaxed population, as indicated by numerical simulations (Kravtsov et al. 2006; Kay et al. 2012), Eq. (A.9) gives a conservative estimate of the intrinsic scatter of the Y_{500}–M_{500} relation.
The intrinsic dispersion about our baseline Y_{X}– relation (Eq. (A.1)), taken from Arnaud et al. (2010), is not measurable; neither is it measurable for the bestfit ChandraY_{X}– relation published in Vikhlinin et al. (2009a). Using a smaller sample of 10 systems, Arnaud et al. (2007) measured an intrinsic scatter of (9 percent), in excellent agreement with the results of the simulations of Nagai et al. (2007) for the scatter of the –Y_{X} relation for relaxed clusters (8.7 percent, their Table 4). It is somewhat larger than the intrinsic scatter of the relation between the true mass and Y_{X} derived by Kravtsov et al. (σ_{log M500  YX} = 5 − 7 percent) but close to the results of Fabjan et al. (2011), who find σ_{log M500  YX} = 0.036–0.046. We thus take as a conservative estimate . The intrinsic dispersion about the Y_{500}– relation for our data is . This value is three times larger than the results of Kay et al. (2012). Partly this is due to the presence of outliers in our dataset (as discussed in Paper I), and it may also be due to projection effects in observed data sets (Kay et al. 2012).
Our final observational estimate of the intrinsic scatter is then σ_{log Y500  M500} < 0.074 or 18 percent, similar to the predictions from Kay et al. (2012) and Sehgal et al. (2010). These predictions depend both on the numerical scheme and specific physics assumptions, with values varying by a factor of two in the typical range 0.04 to 0.08 (references in Sect. A.4.1 below).
Appendix A.4: Assessing the bias from comparison with numerical simulations
The final piece of the jigsaw consists of assessing the bias b in Eq. (A.2). Since the relation has been calibrated using the HE mass for a sample of relaxed clusters, b represents the bias between and the true mass for this category of clusters. In principle, this can be assessed through comparison with numerical simulations. However, this approach is hampered by two difficulties. The first is the exact definition of “relaxed”, since it is almost impossible to select such clusters from observations and simulations according to the same criteria. The second is the specific implementation of the HE equation, which can differ substantially between observations (e.g., the use of forward fitting using parametric models, etc.) and simulations (e.g., the use of mock observations, etc.). Thus the amplitude of the bias that is found will depend not only on physical departures from HE, but also on technical details in the approach to data analysis.
Here we use a different approach that avoids these pitfalls, assessing the bias b by comparing directly the estimated Y_{500}–M_{500} relations with those found from numerical simulations. We then discuss the consistency of the resulting bias estimate with the HE bias expected from simulations and from absolute calibration uncertainties.
Appendix A.4.1: Comparison of simulated Y_{500}–M_{500} relations and data
We first compared the Y_{500}–M_{500} relations from 14 different analyses done by seven groups (Nagai 2006; Yang et al. 2010; Sehgal et al. 2010; Krause et al. 2012; Battaglia et al. 2012; Kay et al. 2012; Sembolini et al. 2013). We translated these simulations results into a common cosmology and, where necessary, converted cylindrical relations into spherical measurements assuming a ratio of Y_{500,cyl}/Y_{500,sph} = 0.74/0.61 ≃ 1.2, as given by the Arnaud et al. (2010) pressure profile.
The lefthand panel of Fig. A.2 shows the different Y_{500}–M_{500} relations rescaled to our chosen cosmology. The simulations use various different types of input physics, and the resulting Y_{500}–M_{500} relations depend strongly on this factor. The only obvious trend is a mild tendency for adiabatic simulations to find nearly selfsimilar slopes (1.66). Runs with nongravitational processes tend to find slightly steeper slopes, but this is not always the case (e.g., the Krause et al. 2012 simulations). The righthand panel of Fig. A.2 shows the ratio of each simulation Y_{500}–M_{500} relation to the PlanckY_{500}– relation given in Eq. (A.8). All results have been rescaled to account for the differences in baryon fraction between simulations. At our reference pivot point of M_{500} = 6 × 10^{14}M_{⊙}, all simulations are offset from the measured relation. There is also a clear dependence on mass arising from the difference in slope between the majority of the simulated relations and that of the Planck relation. The Planck slope is steeper, possibly indicating the stronger effect of nongravitational processes in the real data.
Appendix A.4.2: Quantification of the mass bias
We define the mass bias b between the “true” and observed M_{500} values, following Eq. (A.2) and explicitly allowing for possible mass dependence of the bias, i.e, . Both masses are defined at a fixed density contrast of 500, so that the relations between observed and true mass and radius read
where “true” denotes simulated quantities, and “obs” denotes quantities estimated at the apertures derived from observations. The corresponding Y_{500}–M_{500} relations are In our case, Y_{500} is measured interior to as opposed to . The ratio depends on the radial variation of Y_{500} for scaled radii, , which is close to 1. For a GNFW universal profile (Arnaud et al. 2010), we find that it can be well fit by a power law of the form (1 − b)^{−1/4}. Combining Eq. (A.10), Eq. (A.12) and Eq. (A.13) we then arrive at (A.14)The relation makes it clear that a mass dependence of the bias naturally translates into a different slope of the observed and true Y_{500}–M_{500} relations.
Fig. A.3 Dependence of (1 − b) on mass. Note that this value is strongly dependent on the baryon fraction f_{b} (see text). 
The bias b can then be estimated from a comparison of observed and simulated relations, with the caveat that differences can also arise from imperfect modelling of cluster physics within the simulations. For the ensemble of simulations shown in Fig. A.2, the right panel shows the ratio of the observed and simulated relations as a function of mass. Figure A.3 shows the corresponding variation of (1 − b) as a function of mass M from Eq. (A.14) for the observed slope α = 1.79. This is mass dependent due to the difference in slopes between the simulated and observed relations. At a pivot point of M_{500} = 6 × 10^{14}M_{⊙}, the median value of A_{true}/A_{obs} is 0.74, implying (1 − b) = 0.81. However, there is a large amount of scatter in the predictions from simulations. As a consequence, (1 − b) can vary from 0.74 to 0.97 at M_{500} = 6 × 10^{14}M_{⊙}. Note that the above results depend significantly on the baryon fraction f_{b}. For example, assuming the WMAP7 value f_{b} = 0.167, the median value of (1 − b) is 0.86 at the pivot point of M_{500} = 6 × 10^{14}M_{⊙}.
Appendix A.4.3: Consistency with HE bias predictions and absolute calibration uncertainties
Taken at face value, the bias we derive above of (1 − b) ≃ 0.8 implies that the HE mass used to calibrate the Y_{500}– relation is offset from the true mass by around 20 percent. Is this reasonable?
We can first compare HE Xray and weak lensing (WL) masses. Although as mentioned above both measurements are expected to be biased, such comparisons are useful because the mass measurements involved are essentially independent. In addition measurements for moderately large sample sizes (tens of systems) are now starting to appear in the literature. However, at present there is little consensus, with some studies finding good agreement between HE and WL masses (e.g., Vikhlinin et al. 2009a; Zhang et al. 2010), some finding that HE masses are lower than WL masses, (e.g., Mahdavi et al. 2008), and some even finding that HE masses are higher than WL masses (Planck Collaboration Int. III 2013). The key point in such analyses is rigorous data quality on both the Xray and optical sides. Most recent work points to relatively good agreement between Xray and WL masses, with M^{HE}/M^{WL} ≃ 0.9 on average, and M^{HE}/M^{WL} ≃ 1 for relaxed systems (Mahdavi et al. 2013; von der Linden et al. 2014).
According to cosmological numerical simulations, the measurement bias induced by Xray measurements relative to the “true” values can be caused by two main effects. The first is the classic “hydrostatic mass bias” due to nonthermal pressure support from turbulence/random motions, etc. However, the exact details are very modeldependent. The HE bias expected from simulations varies substantially, depending on the details of the numerical scheme, the input physics, and the approach used to calculate the HE masses (e.g., Rasia et al. 2012). In addition, the amount of bias is different depending on the dynamical state of the object, relaxed systems having less bias than unrelaxed systems. The majority of numerical simulations predict HE biases of 10 to 20 percent (Nagai et al. 2007; Piffaretti & Valdarnini 2008; Lau et al. 2009; Kay et al. 2012; Rasia et al. 2012).
Temperature inhomogeneities constitute the second contributor to Xray measurement bias. In the presence of large amounts of cool gas, a singletemperature fit to a multitemperature plasma will yield a result that is biased towards lower temperatures (e.g., Mazzotta et al. 2004). The presence of temperature inhomogeneities will depend on the dynamical state. While this effect can be investigated with simulations, estimates of its impact vary widely, owing to differences in numerical schemes and the different implementations of the input physics. For instance, simulations with heat conduction consistently predict smoother temperature distributions, thus Xray spectroscopic biases are minimal in this case. On the other hand, “adiabatic” simulations predict longlasting highdensity coolcore type phenomena, which will lead to significant biases in singletemperature fits. Estimates of biasing due to temperature inhomogeneities can range up to 10 or 15 percent (e.g., Rasia et al. 2012).
Finally, for HE mass estimates obtained from Xray data, instrument calibration uncertainties also play a significant role in introducing uncertainties in mass estimates. For instance, the difference in calibration between XMMNewton and Chandra can induce differences in Y_{X}. This is typically 5 percent, from a comparison of XMMNewton based values published by Planck Collaboration XI (2011) to Chandra values for 28 ESZ clusters by Rozo et al. (2012). This can lead to differences of up to 10 percent in the mass derived from Y_{X}, owing to the dependence of the mass on Y_{X}.
Thus our adopted baseline value of (1 − b) ≃ 0.8, ranging from 0.7–1, appears to encompass our current ignorance of the exact bias.
Appendix A.5: Conclusions
In summary the baseline is (A.15)with an intrinsic scatter of σ_{log Y} = 0.075 and a mean bias . The statistical uncertainty on the normalization is about 5% and the error budget is dominated by the systematic uncertainties.
All Tables
Bestfit cosmological parameters for various combinations of data and analysis methods.
All Figures
Fig. 1 Distribution on the sky of the Planck SZ cluster subsample used in this paper, with the 35% mask overlaid. 

In the text 
Fig. 2 Blind S/N versus S/N reextracted at the Xray position using the Xray size, for the MMF3 detections of Planck clusters that are associated with known Xray clusters in the reference cosmological sample. In contrast to the blind SZ flux, the blind S/N is in good agreement with S/N measured using Xray priors. 

In the text 
Fig. 3 Limiting mass as a function of z for the selection function and noise level computed for three zones (deep, blue; medium, cyan; shallow, green), and on average for the unmasked sky (dashed black). 

In the text 
Fig. 4 Noise map σ_{Y500}(θ_{500}) for θ_{500} = 6arcmin. The PSZ is limited by instrumental noise at high ( b  > 20deg) Galactic latitude (deeper at ecliptic poles) and foreground noise at low Galactic latitude. The scale of the map ranges from 0.5 to 2 times the mean noise of the map, which is ⟨ σ_{Y500}(6arcmin) ⟩ = 2.2 × 10^{4} arcmin^{2}. 

In the text 
Fig. 5 Completeness averaged over the unmasked sky as a function of Y_{500} for two different filter sizes, θ_{500} = 6 and 15.3 arcmin. The dashed lines show the semianalytic approximation of Eq. (10), while the diamonds with errors show the completeness estimated by the MC injection technique. 

In the text 
Fig. 6 Planck SZ constraints (+BAO+BBN) on ΛCDM cosmological parameters in red. The black lines show the constraints upon substituting HST constraints on H_{0} for the BAO constrainsts. Contours are 68 and 95% confidence levels. 

In the text 
Fig. 7 Distribution in redshift for the Planck cosmological cluster sample. The observed number counts (red), are compared to our bestfit model prediction (blue). The dashed and dotdashed histograms are the bestfit models from the Planck SZ power spectrum and Planck CMB power spectrum fits, respectively. The cyan long dashed histogram is the best fit CMB and SZ when the bias is free (see Sect. 6.3). The uncertainties on the observed counts, shown for illustration only, are the standard deviation based on the observed counts, except for empty bins where we show the inferred 84% upper limit on the predicted counts assuming a Poissonian distribution. See Sect. 6 for more discussion. 

In the text 
Fig. 8 Comparison of the constraints using the mass functions of Watson et al. (black) and Tinker et al. (red), with b fixed to 0.8. When relaxing the constraints on the evolution of the scaling law with redshift (blue), the contours move along the degeneracy line. Allowing the bias to vary uniformly in the range [ 0.7,1.0 ] enlarges the constraints perpendicular to the σ_{8}–Ω_{m} degeneracy line due to the degeneracy of the number of clusters with the mass bias (purple). Contours are 95% confidence levels here. 

In the text 
Fig. 9 95% contours for different robustness tests: MMF3 with S/N cut >7 in red; MMF3 with S/N cut >8 in blue; and MMF1 with S/N cut >7 in black; and MMF3 with S/N cut at 7 but adopting the MC completeness in purple. 

In the text 
Fig. 10 Comparison of constraints (68% confidence interval) on σ_{8}(Ω_{m}/ 0.27)^{0.3} from different experiments of large–scale structure (LSS), clusters, and CMB. The solid line ACT point assumes the same universal pressure profile as this work. Probes marked with an asterisk have an original power of Ω_{m} different from 0.3. See text and Table 3 for more details. 

In the text 
Fig. 11 2D Ω_{m}–σ_{8} likelihood contours for the analysis with Planck CMB only (red); Planck SZ and BAO and BBN (blue) with (1 − b) in [ 0.7,1 ]. 

In the text 
Fig. 12 Marginalized posterior distribution for ∑ m_{ν} from: Planck CMB data alone (black dotted line); Planck CMB and SZ with 1 − b in [ 0.7,1 ] (red); Planck CMB and SZ and BAO with 1 − b in [ 0.7,1 ] (blue); and Planck CMB and SZ with 1 − b = 0.8 (green). 

In the text 
Fig. A.1 Best scaling relation between Y_{500} and M_{500}, and the data points utilized after correction of the Malmquist bias. 

In the text 
Fig. A.2 Left: comparison of Y_{500}–M_{500} relations from 12 simulations undertaken by six different groups with the updated observational Y_{500}– result from Planck, Eq. (A.8). Right: ratio of each simulated Y_{500}–M_{500} relation relative to Eq. (A.8). The different scaling laws are taken from Kay et al. (2012), Battaglia et al. (2012), Yang et al. (2010), Sehgal et al. (2010), Krause et al. (2012), Nagai (2006), Sembolini et al. (2013) and Planck Collaboration XI (2011). 

In the text 
Fig. A.3 Dependence of (1 − b) on mass. Note that this value is strongly dependent on the baryon fraction f_{b} (see text). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.