Issue 
A&A
Volume 550, February 2013



Article Number  A131  
Number of page(s)  24  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201220040  
Published online  07 February 2013 
Planck intermediate results
V. Pressure profiles of galaxy clusters from the SunyaevZeldovich effect^{⋆}
^{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} Agenzia Spaziale Italiana Science Data Center, c/o ESRIN, via Galileo Galilei, Frascati, Italy
^{5} Agenzia Spaziale Italiana, Viale Liegi 26, Roma, Italy
^{6} Astrophysics Group, Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK
^{7} Atacama Large Millimeter/submillimeter Array, ALMA Santiago Central Offices, Alonso de Cordova 3107, Vitacura, Casilla 763 0355, Santiago, Chile
^{8} CITA, University of Toronto, 60 St. George St., Toronto, ON M5S 3H8, Canada
^{9} CNRS, IRAP, 9 Av. colonel Roche, BP 44346, 31028 Toulouse Cedex 4, France
^{10} California Institute of Technology, Pasadena, California, USA
^{11} Centre of Mathematics for Applications, University of Oslo, Blindern, Oslo, Norway
^{12} Centro de Astrofísica, Universidade doPorto, Rua das Estrelas, 4150762 Porto, Portugal
^{13} Centro de Estudios de Física del Cosmos de Aragón (CEFCA), Plaza San Juan, 1 planta 2, 44001 Teruel, Spain
^{14} Computational Cosmology Center, Lawrence Berkeley National Laboratory, Berkeley, California, USA
^{15} Consejo Superior de Investigaciones Científicas (CSIC), Madrid, Spain
^{16} DSM/Irfu/SPP, CEASaclay, 91191 GifsurYvette Cedex, France
^{17} DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 327, 2800 Kgs. Lyngby, Denmark
^{18} Département de Physique Théorique, Université de Genève, 24 Quai E. Ansermet, 1211 Genève 4, Switzerland
^{19} Departamento de Física Fundamental, Facultad de Ciencias, Universidad de Salamanca, 37008 Salamanca, Spain
^{20} Departamento de Física, Universidad de Oviedo, Avda. Calvo Sotelo s/n, Oviedo, Spain
^{21} Department of Astronomy and Geodesy, Kazan Federal University, Kremlevskaya Str. 18, 420008 Kazan, Russia
^{22} Department of Astrophysics, IMAPP, Radboud University, PO Box 9010, 6500 GL Nijmegen, The Netherlands
^{23} Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, British Columbia, Canada
^{24} Department of Physics and Astronomy, Dana and David Dornsife College of Letter, Arts and Sciences, University of Southern California, Los Angeles, CA 90089, USA
^{25} Department of Physics and Astronomy, University of Sussex, Brighton BN1 9QH, UK
^{26} Department of Physics, Gustaf Hällströmin katu 2a, University of Helsinki, Helsinki, Finland
^{27} Department of Physics, Princeton University, Princeton, New Jersey, USA
^{28} Department of Physics, University of California, Berkeley, California, USA
^{29} Department of Physics, University of California, One Shields Avenue, Davis, California, USA
^{30} Department of Physics, University of California, Santa Barbara, California, USA
^{31} Department of Physics, University of Illinois at UrbanaChampaign, 1110 West Green Street, Urbana, Illinois, USA
^{32} Department of Statistics, Purdue University, 250 N. University Street, West Lafayette, Indiana, USA
^{33} Dipartimento di Fisica e Astronomia G. Galilei, Università degli Studi di Padova, via Marzolo 8, 35131 Padova, Italy
^{34} Dipartimento di Fisica, Università La Sapienza, P. le A. Moro 2, Roma, Italy
^{35} Dipartimento di Fisica, Università degli Studi di Milano, via Celoria 16, Milano, Italy
^{36} Dipartimento di Fisica, Università degli Studi di Trieste, via A. Valerio 2, Trieste, Italy
^{37} Dipartimento di Fisica, Università di Ferrara, via Saragat 1, 44122 Ferrara, Italy
^{38} Dipartimento di Fisica, Università di Roma Tor Vergata, via della Ricerca Scientifica 1, Roma, Italy
^{39} Dipartimento di Matematica, Università di Roma Tor Vergata, via della Ricerca Scientifica 1, Roma, Italy
^{40} Discovery Center, Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark
^{41} Dpto. Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain
^{42} European Southern Observatory, ESO Vitacura, Alonso de Cordova 3107, Vitacura, Casilla 19001, Santiago, Chile
^{43} European Space Agency, ESAC, Planck Science Office, Camino bajo del Castillo, s/n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, Madrid, Spain
^{44} European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands
^{45} GEPI, Observatoire de Paris, Section de Meudon, 5 Place J. Janssen, 92195 Meudon Cedex, France
^{46} Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University of Helsinki, Helsinki, Finland
^{47} INAF – Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, Padova, Italy
^{48} INAF – Osservatorio Astronomico di Roma, via di Frascati 33, Monte Porzio Catone, Italy
^{49} INAF – Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, Trieste, Italy
^{50} INAF Istituto di Radioastronomia, via P. Gobetti 101, 40129 Bologna, Italy
^{51} INAF/IASF Bologna, via Gobetti 101, Bologna, Italy
^{52} INAF/IASF Milano, via E. Bassini 15, Milano, Italy
^{53} INFN, Sezione di Roma 1, Universit‘a di Roma Sapienza, Piazzale Aldo Moro 2, 00185, Roma, Italy
^{54} IPAG: Institut de Planétologie et d’Astrophysique de Grenoble, Université Joseph Fourier, Grenoble 1/CNRSINSU, UMR 5274, 38041 Grenoble, France
^{55} IUCAA, Post Bag 4, Ganeshkhind, Pune University Campus, Pune 411 007, India
^{56} Imperial College London, Astrophysics group, Blackett Laboratory, Prince Consort Road, London, SW7 2AZ, UK
^{57} Infrared Processing and Analysis Center,California Institute of Technology, Pasadena, CA 91125, USA
^{58} Institut Néel, CNRS, Université Joseph Fourier Grenoble I, 25 rue des Martyrs, Grenoble, France
^{59} Institut Universitaire de France, 103 bd SaintMichel, 75005 Paris, France
^{60} Institut d’Astrophysique Spatiale, CNRS (UMR8617) Université ParisSud 11, Bâtiment 121, Orsay, France
^{61} Institut d’Astrophysique de Paris, CNRS (UMR7095), 98 bis Boulevard Arago, 75014 Paris, France
^{62} Institut de Ciències de l’Espai, CSIC/IEEC, Facultat de Ciències, Campus UAB, Torre C5 par2, Bellaterra 08193, Spain
^{63} Institute for Space Sciences, BucharestMagurale, Romania
^{64} Institute of Astro and Particle Physics, Technikerstrasse 25/8, University of Innsbruck, 6020 Innsbruck, Austria
^{65} Institute of Astronomy and Astrophysics, Academia Sinica, Taipei, Taiwan
^{66} Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
^{67} Institute of Theoretical Astrophysics, University of Oslo, Blindern, Oslo, Norway
^{68} Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, La Laguna, Tenerife, Spain
^{69} Instituto de Física de Cantabria (CSICUniversidad de Cantabria), Avda. de los Castros s/n, Santander, Spain
^{70} Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, California, USA
^{71} Jodrell Bank Centre for Astrophysics, Alan Turing Building, School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
^{72} Kavli Institute for Cosmology Cambridge, Madingley Road, Cambridge, CB3 0HA, UK
^{73} LAL, Université ParisSud, CNRS/IN2P3, Orsay, France
^{74} LERMA, CNRS, Observatoire de Paris, 61 avenue de l’Observatoire, Paris, France
^{75} Laboratoire AIM, IRFU/Service d’Astrophysique – CEA/DSM – CNRS – Université Paris Diderot, Bât. 709, CEASaclay, 91191 GifsurYvette Cedex, France
^{76} Laboratoire Traitement et Communication de l’Information, CNRS (UMR 5141) and Télécom ParisTech, 46 rue Barrault, 75634 Paris Cedex 13, France
^{77} Laboratoire de Physique Subatomique et de Cosmologie, Université Joseph Fourier Grenoble I, CNRS/IN2P3, Institut National Polytechnique de Grenoble, 53 rue des Martyrs, 38026 Grenoble Cedex, France
^{78} Laboratoire de Physique Théorique, Université ParisSud 11 & CNRS, Bâtiment 210, 91405 Orsay, France
^{79} Lawrence Berkeley National Laboratory, Berkeley, California, USA
^{80} MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStr. 1, 85741 Garching, Germany
^{81} MaxPlanckInstitut für Extraterrestrische Physik, Giessenbachstraße, 85748 Garching, Germany
^{82} MilliLab, VTT Technical Research Centre of Finland, Tietotie 3, Espoo, Finland
^{83} National University of Ireland, Department of Experimental Physics, Maynooth, Co. Kildare, Ireland
^{84} Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark
^{85} Observational Cosmology, Mail Stop 36717, California Institute of Technology, Pasadena, CA 91125, USA
^{86} Optical Science Laboratory, University College London, Gower Street, London, UK
^{87} SISSA, Astrophysics Sector, via Bonomea 265, 34136 Trieste, Italy
^{88} School of Physics and Astronomy, Cardiff University, Queens Buildings, The Parade, Cardiff, CF24 3AA, UK
^{89} Space Research Institute (IKI), Profsoyuznaya 84/32, Moscow, Russia
^{90} Space Research Institute (IKI), Russian Academy of Sciences, Profsoyuznaya Str 84/32, Moscow 117997, Russia
^{91} Space Sciences Laboratory, University of California, Berkeley, California, USA
^{92} Stanford University, Dept of Physics, Varian Physics Bldg, 382 via Pueblo Mall, Stanford, California, USA
^{93} TÜBİTAK National Observatory, Akdeniz University Campus, 07058, Antalya, Turkey
^{94} UPMC Univ Paris 06, UMR7095, 98 bis Boulevard Arago, 75014 Paris, France
^{95} Université Denis Diderot (Paris 7), 75205 Paris Cedex 13, France
^{96} Université de Toulouse, UPSOMP, IRAP, 31028 Toulouse Cedex 4, France
^{97} Universities Space Research Association, Stratospheric Observatory for Infrared Astronomy, MS 2113, Moffett Field, CA 94035, USA
^{98} University Observatory, Ludwig Maximilian University of Munich, Scheinerstrasse 1, 81679 Munich, Germany
^{99} University of Granada, Departamento de Física Teórica y del Cosmos, Facultad de Ciencias, Granada, Spain
^{100} Warsaw University Observatory, Aleje Ujazdowskie 4, 00478 Warszawa, Poland
Received: 17 July 2012
Accepted: 29 October 2012
Taking advantage of the allsky coverage and broadfrequency range of the Planck satellite, we study the SunyaevZeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R_{500}. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3 × R_{500}, i.e., at a density contrast of about 50–100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and Xray measurements inside R_{500}. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the Xray constraints from XMMNewton in the region in which the profiles overlap (i.e., [0.1–1] R_{500}), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is slightly flatter than most predictions from numerical simulations. Combining the SZ and Xray observed profiles into a joint fit to a generalised pressure profile gives bestfit parameters [P_{0},c_{500},γ,α,β ] = [6.41,1.81,0.31,1.33,4.13 ] . Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3 R_{500}. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.
Key words: cosmology: observations / galaxies: clusters: general / galaxies: clusters: intracluster medium / submillimeter: general / Xrays: general
Appendices are available in electronic form at http://www.aanda.org
© ESO, 2013
1. Introduction
In a pure hierarchical gravitational collapse scenario in the concordance ΛCDM cosmology, concentrations of matter (“halos”) are fully characterised by their redshift and mean matter density, which in turn are related to the power spectrum of initial density fluctuations (Peebles 1980). The scalefree dark matter (DM) collapse drives the evolution of halo concentration across cosmic time (see e.g., recent work by Bhattacharya et al. 2011; Gao et al. 2012), and the ensuing similarity yields both a universal dark matter distribution and simple global scaling relations that should describe the entire halo population (Bertschinger 1985; Kaiser et al. 1995; Navarro et al. 1995; Evrard et al. 1996; Navarro et al. 1997; Voit 2005; Arnaud et al. 2005). However, the observable properties of clusters are determined by the visible baryonic component, which is subject to more complex physical processes related to galaxy formation and feedback. As the main baryonic reservoir in massive halos, the hot gas in the intracluster medium (ICM) is the natural target for studying the physical processes at play and their link to the underlying cluster DM content. Modelling and understanding the baryon physics and disentangling the effect of various feedback processes is a very challenging task (see e.g., Borgani & Kravtsov 2011, for a review), underlining the need for appropriate observational constraints.
The ICM attains Xray emitting temperatures due to gravitational heating within the halo potential well. Xray emission is proportional to the square of the gas density, thus it probes the denser regions of the hot gas. The thermal SunyaevZeldovich (SZ) effect (Sunyaev & Zeldovich 1972), due to inverse Compton scattering of CMB photons by the ICM, is proportional to the thermal gas pressure integrated along the line of sight. It is sensitive to density and/or temperature variations, such as shocks and compression. These two independent observational probes are thus complementary and allow us to further constrain the physics of the ICM.
Recent Xray observations based primarily on representative samples have returned a consistent picture of the scaling and structural properties of halos, from high mass clusters down to the low mass group regime (see, e.g., Böhringer et al. 2007; Vikhlinin et al. 2006; Croston et al. 2008; Pratt et al. 2009; Arnaud et al. 2010; Sun et al. 2011; Sun 2012). In parallel, SZ observations with instruments such as the SunyaevZeldovich Array (SZA, Muchovej et al. 2007), the South Pole Telescope (SPT, Carlstrom et al. 2011), the Atacama Cosmology Telescope (ACT, Swetz et al. 2011), and Planck (Tauber et al. 2010)^{1} have recently started to deliver on the promise of SZ observations for cluster studies.
Building on earlier works on smaller cluster samples (Benson et al. 2004; Andersson et al. 2011), recent results from Planck have underlined the consistency between the Xray and SZ view of the ICM within R_{500}^{2}. These constraints were achieved using three different approaches. The first, detailed in Planck Collaboration (2011e), involved binaveraging of an Xray selected sample from the MetaCatalogue of Xray detected Clusters of galaxies (MCXC, Piffaretti et al. 2011). The second, described in Planck Collaboration (2011f), concerned the comparison of SZ measurements of 62 local clusters from the Early Release Compact Source Catalogue (ERCSC, Planck Collaboration 2011b,c) with good quality XMMNewton archive data. Finally Planck Collaboration (2011d) examined the scaling properties of 21 newlydetected Planck clusters confirmed by XMMNewton. These studies have provided wellconstrained scaling relations (e.g., Y_{500} − M_{500}, Y_{500} − L_{X,500}) in the local Universe, including the first measurement of their intrinsic scatter, to be used as a reference for future evolution and cosmology studies. Other dedicated investigations of smaller samples (e.g., Andersson et al. 2011; Sifon et al. 2012) have reached conclusions similar to Planck early results regarding the agreement between the Xray and SZ view of the ICM.
Beyond global properties, similarity of shape in cluster radial quantities has long been shown in Xray observations (e.g., Arnaud et al. 2001; Pointecouteau et al. 2005; Vikhlinin et al. 2006; Croston et al. 2008; Arnaud et al. 2010), in the optical (e.g., Rines & Diaferio 2006; Wojtak & Lokas 2010), and with weaklensing (e.g., Postman et al. 2012). For SZ studies, the radial distribution of the thermal ICM pressure in massive halos is of particular interest. Recent Xray observations of REXCESS, a representative sample of nearby objects, have shown that the ICM pressure distribution, when scaled appropriately, follows an approximately “universal” shape (Arnaud et al. 2010, A10 hereafter) up to a radius of R_{500}. The small differences in shape for the pressure profiles between relaxed and unrelaxed clusters, especially in the central regions, does not seem to have an impact on the integrated thermal content within R_{500}, or on other global properties, as noted by Planck Collaboration (2011f).
A10 also compared their Xray pressure profiles to those predicted from numerical simulations. Although limited to R_{500}, good agreement was seen between radial pressure profile observations and predictions from various numerical simulations within this radius. However, there are presently few observational constraints beyond R_{500}, and consequently in this region the shape of the “universal pressure profile” was extrapolated according to the predictions from numerical simulations (Borgani et al. 2004; Nagai et al. 2007b; Piffaretti & Valdarnini 2008).
The thermodynamical state of the gas beyond R_{500} bears the signature of the complex physics taking place in the outer parts of the clusters. Characterisation of the gas in cluster outskirts is also necessary to unveil the level of thermal pressure in the cluster periphery, thus constraining their dynamical state. These constraints are crucial for our understanding of the formation and evolution of massive halos. In the cluster surroundings, accretion along filaments is threedimensional and nonspherical (Tozzi & Norman 2001) and the initial conditions of the accretion shock are driven by the thermodynamical state of the (preshocked) infalling material (Voit et al. 2002, 2003). The physical origin of any possible preheating is still unclear, but it may smooth the continuous accretion of substructures. Thus the entropy production at the accretion shock may be boosted (Voit & Ponman 2003; Borgani et al. 2005). Feedback may also improve the degree of thermalisation of the ICM after the shock (with a residual kinetic energy of about 10% beyond R_{500}, Kay et al. 2004).
Observations with Xray telescopes have only recently started to provide insight into the physical properties of the gas in the cluster outskirts beyond R_{500} (George et al. 2009; Reiprich et al. 2009; Urban et al. 2011). The steep decline of the Xray emission with clustercentric radius makes such observations extremely challenging with current instruments. A deep Xray study of the Perseus cluster with Suzaku by Simionescu et al. (2011) indicated that the gas mass fraction in the cluster outskirts was overestimated, possibly due to the effect of gas clumping and the nonvirialised state of the ICM in these regions, as suggested in some numerical simulations (Nagai et al. 2007b).
The different sensitivity of the SZ effect to the radial ICM distribution means that SZ observations have the potential to contribute greatly to the discussion on cluster outskirts. The radial pressure distribution of the first SZ cluster samples have recently been presented based on observations by SPT (Plagge et al. 2010), ACT (Sehgal et al. 2011) and SZA/CARMA (Bonamente et al. 2012). These studies confirmed that the ICM properties as seen by SZ and Xray observations are consistent at least out to R_{500}. Beyond this radius the SZ effect offers the interesting possibility to further constrain the thermal pressure support. Planck is highly competitive in this regard. It is the only SZ experiment with a full sky coverage, able to map even nearby clusters to their outermost radii and offering the possibility of an indepth statistical study through the combination of many observations. In this paper we present constraints on the thermal pressure support derived using SZ observations from the Planck survey. Following our previous methodology for scaling relations, we investigate these issues from a statistical point of view, working with data from 62 local clusters selected from the Planck ESZ sample (Planck Collaboration 2011c), for which there are good quality XMMNewton archival data.
The paper is organised as follows. In the next section we describe the Planck mission and data, together with the sample of galaxy clusters used in this study. In Sect. 3 we recall the basic formalism on the SZ effect and the parametrisaton of the cluster pressure profile we use throughout the paper. Section 4 is devoted to a detailed description of the processing involved in the reconstruction of the SZ and pressure profiles from Planck data. In Sects. 5 and 6 we present the stacked profile of our sample and its best analytical representation. We discuss the comparison with other observational and theoretical constraints in Sect. 7, before presenting our conclusions in Sect. 8.
Throughout the paper we adopt a ΛCDM cosmology with H_{0} = 70 km s^{1} Mpc^{1}, Ω_{M} = 0.3 and Ω_{Λ} = 0.7 . The quantity E(z) is the ratio of the Hubble parameter at redshift z to its present value, H_{0}, i.e., E(z)^{2} = Ω_{M}(1 + z)^{3} + Ω_{Λ}.
2. Data
2.1. Planck data
Planck (Tauber et al. 2010; Planck Collaboration 2011a) is the third generation space mission to measure the anisotropy of the CMB. Planck observes the sky in nine frequency bands covering 30–857 GHz, with high sensitivity and angular resolution from 31′ to 5′. The Low Frequency Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; Mennella et al. 2011) covers the 30, 44, and 70 GHz bands with amplifiers cooled to 20 K. The High Frequency Instrument (HFI; Lamarre et al. 2010; Planck HFI Core Team 2011a) covers the 100, 143, 217, 353, 545, and 857 GHz bands with bolometers cooled to 0.1 K. Early astrophysics results, based on data taken between 13 August, 2009 and 7 June, 2010 (Planck HFI Core Team 2011b; Zacchei et al. 2011), are given in Planck Collaboration 2011c,d,e,f,g. Intermediate astrophysics results are now based on data taken between 13 August, 2009 and 27 November, 2010.
We used the full sky maps in the nine Planck frequency bands provided in healpix (Górski et al. 2005) N_{side} = 2048 full resolution. An error map associated to each frequency band is obtained from the difference of the first half and second half of the Planck rings for a given position of the satellite. The resulting jackknife maps are basically free from astrophysical emission, whilst being a good representation of the statistical instrumental noise and systematic error. We adopted a circular Gaussian as the beam pattern for each frequency, as described in Planck HFI Core Team (2011b) and Zacchei et al. (2011). Uncertainties in flux measurements due to beam corrections, map calibrations and uncertainties in bandpasses are expected to be small, as discussed extensively in Planck Collaboration (2011c), Planck Collaboration (2011d) and Planck Collaboration (2011e).
2.2. Cluster sample
The first allsky coverage by the Planck satellite (Planck Collaboration 2011a) has allowed the detection of dozens of clusters via their SZ signature on the CMB. The Early Release SZ sample (ESZ, Planck Collaboration 2011c) comprises 189 clusters of galaxies characterised over the submillimetre to centimetre wavelength range. The sample consists of SZ clusters and candidates detected with signaltonoise ratios (S/N) spanning from 6 to 29 in the first allsky survey. The sample was thoroughly validated by Planck internal quality assessment, external Xray and optical data crosscorrelations, and a multifrequency followup programme for confirmation. The ESZcluster sample spans over a decade in mass, from 0.9 to 15 × 10^{14} M_{⊙}, which is essentially the full galaxy cluster mass range.
In the following, we focus on a subsample of 62 clusters from the ESZ catalogue. This subsample is defined and extensively described in Planck Collaboration (2011f), where it was used to calibrate scaling relations between the SZ and Xray cluster properties. All 62 clusters were already known in Xrays (found in the MCXC, Piffaretti et al. 2011), and have good XMMNewton archive data, allowing for a high quality Xray data analysis. Masses and radii for the sample were estimated using the M_{500}–Y_{X} relation of A10 (see also Pratt et al. 2010), assuming standard evolution (see Eq. (9)).
The radius R_{500} was calculated iteratively as described in Kravtsov et al. (2006). While this sample is neither representative nor complete, it is the largest, highestquality SZXray data set currently available. The majority of objects lie at a redshift lower than 0.3 (and all have z < 0.5) and cover approximately a decade in mass. In angular size, i.e., θ_{500}, they range between 3.7 and 22.8 arcmin, with a median value of 7.6 arcmin. These extended, SZbright clusters are thus ideal targets to investigate the spatial distribution of ICM thermal pressure support in clusters of galaxies by means of the spatial distribution of their SZ signal.
3. Basic formalism
3.1. The SunyaevZeldovich effect
The inverse Compton effect for thermalised electrons under a blackbody radiation field is treated in the Kompaneets (1957) equation. It was studied and characterised for the case of thermal electrons in clusters of galaxies by Sunyaev & Zeldovich (1970, 1972). The ensuing spectral distortion of the CMB spectrum in the direction of clusters is named the SunyaevZeldovich (SZ) effect.
In the following, we leave aside the kinetic SZ effect, a doppler effect, resulting from the cluster peculiar motion within the comoving reference frame of the Hubble flow. Also, the relativistic corrections on the SZ spectrum, of the order of kT_{e}/m_{e}c^{2}, are not relevant for our study and are thus neglected (see, e.g., Pointecouteau et al. 1998; Challinor & Lasenby 1998; Sazonov & Sunyaev 1998).
The intensity of the SZ effect is characterised by the dimensionless Comptonisation parameter y, the product of the average fractional energy transferred per collision and the average number of collisions: (1)where σ_{T} is the Thomson crosssection, m_{e} the mass of the electron and c the speed of light. P is the pressure produced by the plasma of thermal electrons along the line of sight.
Assuming the clusters are spherically symmetric, we can express the profile in the Comptonisation parameter as a geometrical projection of the spherical pressure profile along the axis of a cylinder: (2)As y is dimensionless we have y(θ) ≡ y(r). Actual SZ measurements derive from the convolution of the yprofile on the sky with the instrument spatial response, f_{PSF}: (3)If the SZ brightness, which is proportional to y, is indeed independent of redshift, the SZ flux is proportional to the integrated Comptonisation parameter and thus depends on the source distance via (4)where D_{A} is the angular diameter distance.
In the following, the observed value of the integrated Comptonisation parameter is given: in units of Mpc^{2} when expressed in the source intrinsic reference frame, i.e., Y(R); and in arcmin^{2} when expressed in the Planck satellite reference frame, i.e. Y(θ). Analogously we can define the total integrated pressure within a sphere of radius r, and express it in Y units. We define this “pseudo” or “spherical” integrated Comptonisation parameter as (5)
3.2. The scaled pressure profile
As mentioned in Sect. 1, in a hierarchical scenario of structure formation the halo population is selfsimilar in scale and structure. Profiles of physical quantities are universal once scaled according to their radius and reference quantities defined at a given density contrast, δ. Within this selfsimilar framework, we adopt throughout the paper the value of δ = 500. The scaled pressure profile thus reads (6)We adopt for the pressure profile the analytical formulation given by Nagai et al. (2007a) for the generalised NavarroFrenk White (Navarro et al. 1997; Nagai et al. 2007a, GNFW) profile (7)where x = r/R_{500}, and the model is defined by the following parameters: P_{0}, normalisation; c_{500}, concentration parameter defined at the characteristic radius R_{500}; and the slopes in the central (x ≪ 1/c_{500}), intermediate (x ~ 1/c_{500}) and outer regions (x ≫ 1/c_{500}), given by γ, α and β, respectively.
A10 have fitted this analytical profile to a combination of observed pressure profiles derived for the REXCESS sample (Böhringer et al. 2007) from high quality XMMNewton data, together with three sets of predicted profiles from numerical simulations of structure formation implementing DM and baryon physics (i.e., radiative cooling and recipes for feedback). Measurements of the gas mass (integrating the density profile) and the Xray spectroscopic temperature lead to the quantity Y_{X,500} = M_{gas,500}T_{X,500}. This parameter links to the actual integrated parameter as seen from the SZ, i.e., Y_{500}, as (8)where μ_{e} = 1.148, is the mean molecular weight of electrons for a 0.3 solar abundance plasma and m_{p} is the proton mass. Here A_{XSZ} = 0.924 ± 0.004 from Eq. (19) of A10 derived from the REXCESS sample, and A_{XSZ} = 0.95 ± 0.04 from the best fit SZ scaling relation between Y_{500} and Y_{X,500} from Planck Collaboration (2011f) using the same sample as in this work. For consistency purposes we use this latest value in the following. The predicted and measured values for Y_{500} are compatible within the ± 1σ limits.
Y_{X,500} is a good proxy for the cluster total mass (Kravtsov et al. 2006; Nagai et al. 2007a; Arnaud et al. 2007). The nonstandard scaling relation fitted against the REXCESS data is provided by their Eq. (2): (9)In turn the characteristic pressure P_{500} scales with the cluster total mass, reflecting the mass variation expected in the standard selfsimilar model, purely based on gravitation (Eq. (5) in A10): (10)The M − Y_{X} relation given in Eq. (9) deviates from the standard selfsimilar case (which has a slope of 3/5). From the definition of P_{500}, any deviation from the standard selfsimilar scaling will appear as a variation of the scaled pressure profiles. As shown by A10 in their Eq. (9), this variation can be expressed as a function of the total mass. At δ = 500, this mass dependence is almost constant with radius, and it can be approximated by (11)
3.3. Conversions and normalisations
In the following, we measure SZ profiles from the Planck data (Sect. 4.2). All through this work, we compare these observed profiles with the predictions that arise from the Xray constraints. For all clusters we assume the universal pressure profile shape as derived from A10. We parametrise it according to the quantities derived from the Xray analysis R_{500} and Y_{X,500} (see previous section). We can then derive predicted pressure profiles (see Eq. (7)), and Comptonisation parameter profiles (Eqs. (2) and (4)). The later convolved by the instrument beam (Eq. (3)) is directly comparable to our observed Planck profiles.
We apply a statistical approach to the observed SZ profile and to the pressure profile, averaging individual profiles in our sample once scaled. Profiles are scaled in radius according to R_{500}. The observed SZ profiles (or predicted SZ profiles from Xray constraints) are dimensionless, and we therefore normalise them by the quantity (12)In the following we used the value of Y_{500} as defined in Eq. (5).
To translate our measured SZ profiles, we deconvolve and deproject them for each cluster following Eqs. (2) and (3). In practice this operation is performed as described in Sect. 6.1. The resulting profiles are converted to pressure such that (13)
4. Reconstruction of the SZ profile
4.1. SZ signal extraction methods
The thermal SZ maps were recovered from a combination of Planck channels, making use of either real or Fourier space methods. We applied three methods based on internal linear combination (ILC) algorithms to the Planck data: (1) the Modified Internal Linear Combination Algorithm (MILCA, Hurier et al. 2010); (2) Needlet Internal Linear Combination (NILC, Delabrouille et al. 2009); and (3) the Generalized Morphological Component Analysis (GMCA Bobin et al. 2008). A description of each method is given in Appendix A.
To optimise the reconstructed SZ map we have to take into account that: (i) the brightness of the SZ effect (as an increment or decrement) is maximum in the submillimetre to millimetre range; (ii) the typical angular size of the cluster is a few arcmin; and (iii) the final resolution of the reconstructed SZ map is determined by the lowest resolution of the combined frequency maps. We therefore restricted ourselves to the use of the six PlanckHFI channels in the SZ map reconstruction process. The final resolution of the reconstructed SZ maps is that of the lowestresolution 100 GHz channel, or 10 arcmin. SZ maps are in units of the (dimensionless) Comptonisation parameter. Integrated Comptonisation parameters are expressed in arcmin^{2} for the observed values and in Mpc^{2} for the intrinsic value in the source reference frame. The associated SZ noise maps are built in a fully consistent way from the individual frequency error maps (see Sect. 2).
4.2. SZ profile computation
4.2.1. Profile extraction
Individual cluster profiles were computed from the reconstructed allsky SZ maps. Uncertainties were obtained from the allsky reconstructed SZ error maps, which were derived by applying the same reconstruction methods to the frequency error maps (see Sect. 2). For instance, in the case of a linear method (e.g., ILC based), the coefficients used for the linear combination of the frequency maps were propagated to the combination of the frequency error maps.
We extracted a square patch of side 20 × θ_{500} around each cluster position from the allsky SZ map. Patches were projected from healpix (Górski et al. 2005) to a tangential projection. For each patch, the pixel size was adapted so that it was constant in scaled units of θ_{500} over the full cluster sample. This unavoidably leads to an oversampling of the pixels or, more precisely, to a redundancy of the original allsky map pixels in the reprojected patch. We produced an associated patch that tracks this redundancy. We also extracted equivalent patches from the associated allsky SZ error map, and from the allsky variance map.
We then computed a profile from each SZ map patch. These were calculated on a regular radial grid with bins of width Δr/R_{500} = 0.25, allowing us to sample each cluster profile with four points within R_{500}. The y value of a bin was defined as the mean of the values of the pixels falling in each annulus. We subtracted a background offset from the maps prior to the profile computation. The offset value was estimated from the surrounding region of each cluster where r > 7 R_{500}. The uncertainty associated with this baseline offset subtraction was propagated into the uncertainty of each bin of the radial SZ profile.
4.3. Expected SZ profile
We used the GNFW pressure profile shape characterised by A10 from the REXCESS sample and numerical simulations to generate a pressure model for each cluster of our sample. We kept the A10 best fit values of c_{500} = 1.18 and the three slopes, α = 1.05, β = 5.49, γ = 0.31. We derived the normalisation P_{0} from the observed value of Y_{X,500} (see Eq. (5)), using a conversion of 0.95 between Y_{X} and Y_{SZ} as described in Sect. 3.2. Finally R_{500} was fixed to the values derived from the XMMNewton analysis. Both Y_{X} and R_{500} are reported in Planck Collaboration (2011f). We computed the projection matrix and the PSF redistribution matrix (for a Gaussian beam of 10 arcmin FWHM) as expressed in Eqs. (2) and (3), respectively, to multiply each pressure profile model. The derived y profile model is directly comparable to the Planck measurements. The average profile model across the sample was derived similarly to the observed stacked profile (see above), and is used to compare to the observed stacked profile.
4.3.1. Accounting for correlation between points
A certain level of correlation is introduced between the points of a radial profile derived in this manner. We account for this correlation in the covariance matrix of each profile, which is computed as follows.
For each map patch, we masked radii θ < 7 θ_{500} centred on the cluster and computed the power spectrum of the noise on the remaining area. The use of the surrounding regions of the cluster to characterise the noise properties allows us to account for the effect of astrophysical contamination in the cluster vicinity, as well as the instrument noise and systematics. These sources of contamination were, by construction, excluded from the allsky frequency error maps (see Sect. 2). The power spectrum drawn from the SZ error map of each patch is therefore, systematically, slightly lower than the one computed over the cluster surroundings.
We then simulated m = 500 realisations of the noise patch based on this power spectrum, accounting for the variance map and sky pixel redundancy and assuming an inhomogeneous correlated Gaussian noise. We then extracted a noise profile from each realisation, reproducing the baseline background subtraction used for the observed profile. The covariance matrix was built from all the simulated noise profiles, i.e., , where P_{n} is an npoints × m matrix of simulated noise profiles. Considering two points as correlated when their correlation coefficient is larger than 0.3, the typical level of pointtopoint correlation in our profiles is about 20% (16%, 21% and 28% for MILCA, NILC and GMCA, respectively).
Fig. 1 Average scaled SZ profiles for each of the nine Planck frequency bands (decreasing frequency from left to right and from top to bottom). The points within each individual profile are correlated at about the 20% level (see Sect. 4.3.1). At each frequency, the stacked radial profile is obtained from the average of the 62 individual cluster profiles scaled by R_{500} and Φ_{i}/ ⟨ Φ_{i} ⟩ , in the x and yaxis directions, respectively (see Sect. 5.1). The lightred shaded area in each panel indicates the statistical uncertainty about the average observed profile. The solid black line corresponds to the prediction of the universal pressure profile (A10) for our cluster sample in each band, and the two dashed lines indicate its associated dispersion. 
4.3.2. Stacking procedure
Given that all profiles are computed on the same grid in scaled radius, they can easily be stacked together. Each profile y_{i} was rescaled by R_{500,i} and Φ_{i} (as given in Eq. (12)), respectively in the X and Y axis directions. The stacked profile and associated covariance matrix are given by (14)For the computation of Φ, we used Y_{sph,500} in Mpc^{2} as given in Eq. (5). For simulated data (Appendix B), Φ was derived from the Planck M_{500}–Y_{500} relation (Planck Collaboration 2011d); for observed data (Sect. 5.3), Φ was measured directly from the Planck data.
When for a given bin, values for all clusters are dominated by the signal, we assumed a lognormal distribution of their scatter and stacked them in logarithmic space. The logarithm of the stacked value and associated error were then translated back into linear space. In other terms, for bin j, if all the clusters satisfied the condition , we stacked all the measurements in the logarithmic space; otherwise, we did it in the linear space.
We computed the statistical and weighted (i.e., 1/σ^{2}) average profiles, as well as the median profile of the sample, and checked their consistency. Over a radial range of 3 × R_{500}, the average of the relative error of the weighted mean and median profiles with respect to the statistical average profile is smaller than 10% and 5%, respectively. Based on this agreement, we used the statistical average to compute the stacked profile throughout this study.
Of the 62 clusters in our ESZXMM sample, only two are spatially close and thus are potentially physically connected: A3528 and A3532. These two clusters are members of the Shapley supercluster. Each of those clusters was masked over an area of radius 3 R_{500} when processing the other. Thus, we consider that their profiles and covariance matrices are independent.
At the end of our data processing, the stacked covariance matrix encompasses the statistical errors due to instrumental noise, astrophysical fluctuations at the cluster location and systematic effects (e.g., instrumental and arising from data processing). For the purposes of our study we have also computed the dispersion across our sample for each position within our profiles. Both uncertainties are propagated throughout the analysis, and are crosscompared in the following sections. In Figs. 3, 4 and 6 the error bars shown on the Planck data points are purely statistical and correspond to the square root of the diagonal elements of the covariance matrix (i.e., Eq. (14)). The rms scatter of individual profiles around the mean is indicated by a coloured band.
5. Planck galaxy cluster SZ profile
5.1. Frequency stacked profiles
We first looked for the raw cluster signature in each of the nine Planck frequencies. We extracted a brightness profile at each frequency from the raw Planck maps for each of the 62 clusters in the ESZXMM sample. These were rescaled in units of R_{500} in the xaxis direction and in units of Φ_{i}/ ⟨ Φ_{i} ⟩ in the yaxis direction (see Sect. 3.3), and then stacked. The resulting average brightness profiles are shown in Fig. 1, where the shaded area in each panel depicts the statistical error about the mean flux value in each bin.
The cluster signal appears clearly in most of the nine frequencies, following the SZ thermal spectral signature. The signal is positive at 545 and 353 GHz, compatible with zero at 217 GHz, and then negative down to 30 GHz. The 62 clusters of our sample are strong SZ sources, as expected for objects detected at S/N > 6 (Planck Collaboration 2011c). The stacking procedure mostly averages out the effect of foreground and background contaminants; however, they contribute to the final dispersion in the stacked profiles. Owing to their Gaussian nature, CMB fluctuations are more easily washed out in the stacking procedure than the Galactic and extragalactic dust emission in the high frequency range, and Galactic and point source radio emission (i.e., freefree and synchrotron) in the lower frequency range.
We also plotted for each frequency the expected profile computed as described in Sect. 4.3. This comparison shows a first order agreement between the Planck measurements and the expected SZ profiles derived from Xray constraints assuming a GNFW pressure profile shape.
5.2. Stacked SZ image
Figure 2 presents the stacked average image of the 62 clusters of our ESZXMM sample, plotted on a logarithmic scale. Before averaging, each individual SZ map was normalised identically to the profiles (see Eq. (12)) and randomly rotated by 0°, 90°, 180° or 270°. The stacked image was then renormalised by ⟨ Φ_{i} ⟩ in order to be expressed in y units. The image shows that the SZ signal is detected out to ~ 3 × R_{500}, marked by the white circle in the figure. The extent of existing observational constraints on the gas distribution in clusters in SZ and Xrays, i.e., R_{500}, is marked by the black circle. We also generated a jackknifetype map by averaging the scaled individual SZ maps, m_{i}, as ∑ ( − 1)^{i}m_{i}. The rms values of both the stacked SZ and onoff maps beyond 5 × R_{500} are compatible: 4.6 × 10^{7} and 4.8 × 10^{7}, and, in turn, compatible with the rms value below this radius on the jackknife stacked map: 3.5 × 10^{7}. This check demonstrates that our average map is not strongly affected by residuals or biases. The radial profile computed from the stacked map is fully compatible with the stacked profile discussed below.
Fig. 2 Stacked scaled image (size is 10 R_{500} on a side) using a logarithmic stretch in y. Individual maps are rescaled by Φ_{i} before averaging, and then multiplied by ⟨ Φ_{i} ⟩ . The black and white circles mark the loci of 1 and 3 × R_{500}, respectively. 
5.3. Observed stacked SZ profile
We derived the observed stacked Planck SZ profiles for the ESZXMM sample as described in Sect. 4.2. Y_{500} is obtained from the algorithms used for blind detection of SZ clusters in the Planck survey (Planck Collaboration 2011c), namely the Powell Snakes (PwS, Carvalho et al. 2009, 2011) and multifrequency matched filter (MMF3, Melin et al. 2006) algorithms. In both cases the algorithms implement a universal pressure profile shape (A10) with the position of the cluster fixed to the XMMNewton coordinates and the size of the universal pressure profile fixed to R_{500}. The fluxes from both methods (i.e., PwS and MMF3) are consistent over the whole sample. The median value for the ratio PwS/MMF3 is 0.96 ± 0.05 (see Appendix B.2). To further validate the above, and as a consistency check with previous Planck results, we fitted each individual SZ profile with a projected, PSFconvolved universal pressure profile. We fixed R_{500} to the best fitting Xray value from Planck Collaboration (2011f) and only fitted the normalisation, Y_{500}. Given the agreement, we used the MMF3 values to compute Φ for each cluster (Eq. (12)).
Appendix B presents a detailed investigation of the convergence between stacked profiles derived from each of the three different methods of SZ signal reconstruction. We find that all three methods agree remarkably well over the entire radial range, both for simulated and observed SZ profiles. The three SZ reconstruction methods lead to profiles fully compatible with each other. Across the range of radii over which the profiles have been computed (i.e., [0–10] × R_{500}), accounting for the correlated errors of each profile, the reduced χ^{2} of NILC and GMCA with respect to MILCA are 0.48 and 0.62, respectively^{3}.
The left panel of Fig. 3 shows the stacked MILCA profile for the 62 clusters in the ESZXMM sample compared to the individual profiles (see Sect. 4.2). The SZ signal is statistically detected over more than two decades in intensity, and out to a remarkably large radius of ~3 R_{500}, reaching far into the cluster outskirts. Assuming that the virial theorem can be applied (which at such large radius might be breaking down), the outer radius of our statistical detection corresponds to a density contrast of δ ~ 50. More conservatively, we can assume that we are statistically probing the average SZ and pressure distribution down to regions of density contrast of δ ~ 50–100. The dispersion about the mean profile dominates the statistical uncertainties. It is minimal (by construction) at ~20% around R_{500}, but increases towards the centre and the outskirts to ~50% and ~65% at 0.3 and 2 × R_{500}, respectively. At the most external radius of our statistical detection, i.e., ~3 × R_{500}, the dispersion of the individual profiles about the mean is more than 100%, as at these large radii, the individual SZ profiles are fully dominated by noise.
The right panel of Fig. 3 shows the correlation coefficient matrix for the stacked profiles presented on the left panel, and illustrates the degree of correlation between the points in the profiles.
5.4. Comparison with expectations from universal profile
The stacked model used for comparison was computed as described in Sect. 4.3. When considering the statistical errors only, the measured stacked profile is significantly above the model. Taking into account the error on the model and the correlated errors between points of the reconstructed SZ profile, we obtained a reduced χ^{2} value of 3.53 within 3 R_{500}. If we omit the error on the model, this value becomes 7.48 As our tests on simulations show a very good agreement between the input model and the output reconstructed profiles (see Appendix B), the difference observed here between the measured SZ profiles and the predicted model from Xray constraints is not an artefact of the method, but an intrinsic difference. This difference is significant at a 2–3σ level from around R_{500} out to 3 R_{500}.
6. The galaxy cluster pressure profile from Planck
Fig. 3 Left: individual profiles (grey lines) and scaled stacked profile (red points) for our sample of 62 clusters. The lightred shaded area marks the dispersion about the average stacked profile, with its upper and lower limits highlighted by the two dashed red lines. The points within each individual profile are correlated at about the 20% level (see Sect. 4.3.1). The error bars on the Planck data points are purely statistical and correspond to the square root of the diagonal elements of the covariance matrix. The solid black line (labeled in the legend as “Expected from P_{A10}”) is the average stacked profile obtained from the expected SZ individual profiles drawn from the universal pressure profile of A10, parameterised according R_{500} and Y_{500} derived from the XMMNewton data analysis presented in Planck Collaboration (2011f). The two dotted black lines indicate the associated dispersion about this model profile. The bottom panel shows the value of χ at each point with respect to the expectations from the universal profile taking into account the measured errors. Right: matrix of the correlation coefficient for the observed radial stacked profiles. Both panels are illustrated here for the MILCA results. 
Fig. 4 Left: Planck pressure profile obtained from the average of the individual pressure profiles across our sample (red points) shown together with the stacked pressure profile derived from the XMMNewton data for the same sample (purple points). The dispersions about the SZ and Xray profiles are depicted, respectively, by the red and purple shaded areas. Our best fit GNFW profile is shown as a solid black line and that of A10 as a solid green line. The lower panel shows the χ profile of these two best models taking into account the statistical errors and dispersion about the observed profile. Right: stacked profile from Planck and XMMNewton for the subsamples of coolcore (open circles) and non coolcore (open stars) clusters within the ESZXMM sample. The shaded areas are identical to the one shown in the left panel. Our best fit models for each subsample are shown as black solid and dashed lines for the coolcore and non coolcore clusters, respectively (see Table 1). The best A10 fit for coolcore and non coolcore clusters are shown as green solid and dashed lines, respectively. The lower panel gives the χ profiles of our best fit models with respect to the measured profiles and associated errors (including dispersion). On both panels the Planck data points are correlated at about the 20% level (see Sect. 4.3.1). As for the stacked SZ profile, the error bars on the Planck points are purely statistical and correspond to the square root of the diagonal elements of the covariance matrix. For the XMMNewton points they correspond to the statistical error on mean. 
6.1. PSF deconvolution and deprojection
The pressure profile is derived for each cluster in our sample by applying a deconvolution and deprojection to the observed SZ profiles. As already mentioned, we assume spherical symmetry. The convolution by the instrumental beam and the geometrical projection of a spherical pressure profile into a Comptonisation parameter profile are expressed in Eqs. (3) and (2), respectively. We have applied a deconvolution plus deprojection algorithm adapted from the method described by Croston et al. (2006). This method allows us to perform a straightforward real space deconvolution and deprojection using a regularisation procedure originally applied to Xray surface brightness profiles to derive the underlying density profile. The method is adapted to lower statistics (i.e., lower number of bins per radial profile). In the present case, the PSF response matrix was computed for a circular Gaussian beam with FWHM = 10 arcmin, i.e., the angular resolution of our SZ maps. SZ cylindrical profiles deconvolved from the PSF were then deprojected into 3D profiles (assuming spherical symmetry for the clusters), and converted to pressure following Eq. (13).
We used a Monte Carlo (MC) method to propagate the measured error over the SZ profiles onto the PSFcorrected and deprojected profiles, accounting for the covariance between points expressed in the form of the matrix C_{i} for the cluster i. For each cluster i, its covariance matrix was Cholesky decomposed (i.e., ). We assumed correlated Gaussian noise, and generated m = 10 000 realisations of the SZ profile, P_{i}, with respect to this decomposition (i.e., , where x is an array of random numbers following a normal distribution). With this method we make sure we properly sample the noise properties of each profile. Each realisation was then rebinned logarithmically in radius and run through our deconvolution and deprojection algorithm. We conservatively chose to feed to the algorithm values of the weight per point equal to the inverse of the square of the flux dispersion in this bin (i.e., ). The weights remained the same for all realisations. These m realisations of the pressure profile were used to compute the covariance matrix of the cluster pressure profile similarly to what is described in Sect. 4.3.1. With this MC process, we ensure a proper propagation of the correlated noise through the PSFcorrection and deprojection of the SZ profile.
6.2. The stacked pressure profile
As we did for the SZ profiles, the 62 pressure profiles were rescaled in the x and yaxis directions according to R_{500} and P_{500}, and then stacked together.
To further compare with A10’s results, we accounted for the slight difference in mass range between the REXCESS sample and ours via the factor (Eq. (11), and see also Sun et al. 2011). We divided our stacked profile by the average value across our ESZXMM sample, i.e., ⟨ f(M) ⟩ = 1.09.
Our three reconstruction methods lead to compatible stacked pressure profiles. However, we have accounted for their fluctuations by adding to the diagonal terms of the covariance matrix of the stacked profile the maximum pointtopoint difference between the MILCA profile and the two other methods. In parallel to our Planck SZ analysis, we have also derived the pressure profiles of our 62 clusters from the analysis of XMMNewton data. We followed the method applied to the REXCESS sample and presented in A10. Details of the XMMNewton analysis and results are provided in Démocles et al. (in prep.).
The Planck and XMMNewton stacked pressure profiles derived from our sample of 62 ESZXMM clusters, are displayed in Fig. 4. The two sets of data agree remarkably well. They are fully compatible within their respective dispersions, and they overlap in the radial range ~(0.3−1) R_{500}. We recall that the points from the Planck data suffer a degree of correlation at about the 20% level, whereas all the points in the XMMNewton profile can be considered as uncorrelated. These two independently derived pressure profiles are fully complementary. Indeed, in the central region the pressure profiles are very well constrained from the Xray data, whilst the Planck measurements hardly reach down to 0.1 R_{500} due to the moderate spatial resolution of our SZ maps. At larger radii, the extent of the Xray observations is limited to radii smaller than ~(0.7−1.0) R_{500} because of the XMMNewton field of view and of the quick drop in Xray surface brightness with increasing radius. The Planck profile extends far beyond this radius.
With this joint constraint, we bring for the first time a comprehensive observational view of the distribution of the average thermal pressure distribution in clusters of galaxies out to a density contrast of δ ~ 50–100.
Following A10, we have investigated two subsamples of our ESZXMM clusters. Keeping the coolcore (CC) versus non coolcore (nonCC) classification as provided in Planck Collaboration (2011f), we computed the stacked profiles for the 22 CC versus 40 nonCC clusters. Both stacked profiles are displayed in the right panel of Fig. 4. They are, as expected, different in the central parts, with a more peaked profile for CC systems and a shallower one for nonCC clusters. However, in the outer parts, our observed profiles for the two subsamples have very similar slopes.
Fig. 5 Marginalised posterior likelihood distribution for our best MCMC fit to a GNFW pressure profile with free parameters [P_{0},c_{500},α,β] and γ = 0.31. The best fit values of the parameters are marked as white crosses and solid red lines, respectively, on the 2D and 1D posterior distributions. On the 2D distributions, the white triangles mark the value of the A10 parameters, while the red, purple and green solid lines define, respectively, the 68, 95 and 99% confidence levels. 
6.3. Bestfit to the generalised NFW profile
We have combined the Planck and XMMNewton pressure profiles into a joint fit to the generalised NFW model (Eq. (7)). The pressure profiles derived over our sample from the Planck and XMMNewton data cover a radial range from 0.02 × R_{500} to 3 × R_{500}, with an overlap in the radial interval of (0.1−1) × R_{500}. In order to find the best parametric representation of our observed stacked profile, we ran a Monte Carlo Markov chain analysis to find the maximum likelihood solution parametrising the GNFW profile. The Planck and XMMNewton data are two different and fully independent observational data sets. We therefore computed the likelihood as the product of the two independent likelihoods: ℒ = ℒ_{X} × ℒ_{SZ}. For the Planck data, we have accounted for the correlation between points through the use of the covariance matrix (see Sect. 6.1). We also accounted for the dispersion about the average Xray and SZ profiles across the sample by quadratically adding it to the errors bars of the XMMNewton points and to the diagonal elements of the covariance matrix for the Planck profile.
We performed our MCMC analysis in log–log space for various combinations of free parameters. Fixed parameters were assigned the values from the best fit provided by A10. Our best fit parameters are reported in Table 1. In order to find the best possible analytical representation of our data, we assessed the quality of the fit via the value of the reduced χ^{2} using only the statistical errors, . The configuration with four and five free parameters are of equivalent quality, the former being slightly better (i.e., and 1.0, respectively). Unsurprisingly the degeneracy between free parameters increases with their number, translating into a drastic increase in the parameter covariance. We therefore adopted the four free parameter approach as our best parametrisation of a GNFW profile of our data, i.e., [P_{0},c_{500},α,β] free, with γ fixed to 0.31. The associated values are marked in bold face in Table 1, and are given purposefully with high precision with respect to the errors on parameters. These uncertainties are illustrated by the corresponding marginalised posterior likelihood distributions that are shown for individual (1D) and pairs of parameters (2D) in Fig. 5.
Best fit parameters for a generalised NFW pressure profile, with different numbers of fixed parameters.
In parallel, following exactly the same procedure, we also fitted the stacked profiles for the subsamples of coolcore and non coolcore clusters. The best fit parameters are also reported in Table 1. Both the stacked profiles and best fit model are shown on the right panel of Fig. 4.
We further adopted this four free parameter configuration to fit each of our measured individual pressure profiles, combining XMMNewton and Planck pressure data. As describe in the above section for the stacked profile, we accounted for the pointtopoint differences between our three reconstruction methods (see Appendix A) in the error budget. Results are reported in Appendix C, Fig. C.1 and Table C.1.
7. Discussion
7.1. The observed pressure profile
7.1.1. The core of clusters
Our present results are in very good agreement with the universal profile derived by A10 outside the core and down to R_{500}. Within the core, i.e., R < 0.15 R_{500}, our observed profile lies significantly below the A10 profile. In comparison to the REXCESS sample which is Xray selected, our ESZXMM sample is SZ selected (except for the intersection with the XMMNewton archives) and is thus, a priori, closer to being a massselected sample. It therefore contains more dynamically disturbed clusters (e.g., A2163, RXJ2228+2036, etc.). As has also been indicated by results in the validation followup of Planck clusters with XMMNewton (Planck Collaboration 2011d, 2012, 2013b), this seems to indicate that Xray selection undersamples the population of morphologically disturbed clusters. When splitting our sample into CC and nonCC clusters, the average profile for the CC clusters (and its associated best fit) is in better agreement with the universal profile given in A10. However, it is still slightly lower in the central parts than the A10 profile for CC clusters (see right panel of Fig. 4). There might thus still be differences between the population sampled by our CC clusters and those from A10. On the contrary, our nonCC cluster profile agrees well with the A10 nonCC profile, except for the very central parts (R < 0.04 R_{500}).
7.1.2. The inner profile
Observationally, it is hard to accurately determine the distribution of the SZ signal and, thus, of the underlying thermal pressure. Early single target SZ studies were limited to the clusters’ inner regions (i.e., within R_{500}), and were not really competitive with the Xray measurements in terms of resolution and sensitivity (see e.g., Pointecouteau et al. 2001; Halverson et al. 2009; Korngut et al. 2011). However, the combination of SZ and Xray tracers has already shown its potential in terms of structural studies of clusters (Pointecouteau et al. 2002; Kitayama et al. 2004; Jia et al. 2008; Basu et al. 2010). The first studies working on SZ cluster samples were presented by the SPT and ACT collaborations (Plagge et al. 2010; Sehgal et al. 2011), with 15 and 9 high significance clusters, respectively. The SPT team found good agreement between the Xray predicted signal and their SZ measurement within R_{500}. Similarly, in a recent study, Bonamente et al. (2012) used a sample of 25 clusters observed with SZA/CARMA (Muchovej et al. 2007) and found an excellent agreement when modelling the SZ emission over their sample, either with the universal pressure profile from A10 or with the model from Bulbul et al. (2010). Our SZ and Xray data are in excellent agreement over the radial range (0.1–1) R_{500}, providing tighter observational constraints.
Our results, together with the SPT, and ACT results, have to be compared with other works that reported a difference between SZ and Xray measurements. Two papers based on WMAP data reported a measured SZ signal lower than the expected signal from their Xray properties. Based on a WMAP (one year) analysis for 31 randomly selected nearby clusters, Lieu et al. (2006) found a ratio of one fourth between the measured SZ signal and the expected signal from Xray constraints. An analysis based on WMAP7 data for 49 clusters (Komatsu et al. 2011) found this ratio to be 0.5–0.7. However, these claims did not agree with other independent WMAP5 statistical analysis for 893 clusters (Melin et al. 2011). Neither do we, obtaining instead an excellent agreement between Xray and SZ properties within R_{500}, in particular, when using the Planck frequencies overlapping those of WMAP (see profiles for 30, 44 and 70 GHz in Fig. 1). The WMAP fluxes derived by Lieu et al. (2006) and Komatsu et al. (2011) are also discrepant with fluxes from OVRO/BIMA (Bonamente et al. 2006). These differences could be due, for instance, to complex effects involving largescale beams and/or nonlinear gains (e.g. Whitbourn et al. 2011, and references therein).
7.1.3. The outskirts
Beyond R_{500}, the observational situation is even more wanting. Only a few Xray observations with XMMNewton, Chandra and/or Suzaku constrain the density, temperature or gas fraction profiles out to ~R_{200} ≃ 1.4 R_{500} (George et al. 2009; Reiprich et al. 2009; Urban et al. 2011; Simionescu et al. 2011; Walker et al. 2012). This type of Xray detection remains very challenging and requires very long exposure times, as at larger radii the Xray emission is extremely faint. We recall that beyond R_{500} the universal pressure profile from A10 was constrained by predictions from numerical simulations, not by observations. Probing the gas with the SZ effect is therefore a powerful alternative, as shown, for instance, by the SPT average emission over 15 clusters (Plagge et al. 2010) or the WMAP statistical analysis over ~700 clusters by AtrioBarandela et al. (2008), though the latter is affected by the limited resolution and sensitivity of the survey. Conversely, from the tentative analysis of WMAP3 data for 193 clusters with Xray temperature above 3 keV, Afshordi et al. (2007) provided constraints for the cluster pressure profile (significant out to ~1.3 R_{500}). In agreement with these earlier works, the Planck measurements of the present study are the first to allow a precise description of the thermal pressure distribution out to the cluster outskirts.
As shown in Appendix C, Planck also resolves some individual pressure profiles. The work presented by the Planck Collaboration (2013c) on the Coma cluster is even more striking, as the SZ profile detection reaches beyond 3 × R_{500}, i.e., as far as the statistical measurement presented in this work. The derived best fit model for the Planck observations of Coma is in full agreement with our stacked result.
While in the inner parts the profiles of CC and nonCC clusters are significantly different, and although the CC profiles lie just below the nonCC beyond ~1 R_{500}, the two are compatible within our statistical errors. The same is true when comparing our bestfit models and A10’s for CC and nonCC clusters. This suggests that across our sample the average differences between the outskirts of the two types of clusters are smaller than the scatter between clusters.
We have also investigated the change in integrated Comptonisation parameter derived by comparing our best fit of a GNFW model to that of A10. Within a fixed aperture R_{500}, known from ancillary data, we have computed Y_{500} for each cluster in our sample and for both sets of parameters (see Eq. (B.2)). The average ratio across the sample between the value of Y_{500} from our best profile and that of A10’s is 1.02 ± 0.03. The difference is thus marginal as the two profiles are alike, and consistent with the ratio of 1.01 between the two GNFW parametrisations at R_{500}. This consistency demonstrates the robustness of the previouslypublished Planck SZ analysis and scaling relations where the A10 profile was adopted as a fiducial model for known R_{500} values (Planck Collaboration 2011d,e,f, 2012, 2013a,b). The ratio for the values of Y for the two profiles is 1.13 and 1.19 within fixed apertures of 3 R_{500} and 5 R_{500} respectively, as expected given the slightly flatter outer slope derived in this work with respect to that of A10.
7.2. Comparison with theoretical predictions
7.2.1. Sets of numerical data
We have investigated three sets of simulated clusters in order to compare to our combined Planck and XMMNewton pressure profile. All are taken from numerical simulations of structure formation in a ΛCDM cosmology:

1.
We first compared to the set of combined simulations used byA10 to derive their universal pressure profile together with theREXCESS data. They comprise 93, 14, and 88simulated clusters withM > 10^{14} M_{⊙}, from Borgani et al. (2004), Nagai et al. (2007a) and Piffaretti & Valdarnini (2008), respectively. We refer to this set of simulations as B04+N07+P08 hereafter.

2.
The second set is built from the pressure profiles of 64 massive clusters with M > 3 × 10^{14} M_{⊙} from a simulation by Dolag (in prep.).

3.
The third comprises the 40 most massive clusters (i.e., M > 3 × 10^{14} M_{⊙}) provided from the numerical simulations described in Battaglia et al. (2010, 2012).
All simulations include treatment of radiative cooling, star formation and energy feedback from supernova explosions. Simulations by Dolag (in prep.) and Battaglia et al. (2010, 2012) have prescriptions for AGN feedback. The different simulation sets use different techniques and different implementations of the physical processes. This ensures a fair description of uptodate theoretical predictions, hopefully bracketing the plausible range on the thermal pressure profiles distribution in clusters of galaxy.
We kept the B04+N07+P08 sample as used by A10, even though its mass limit is slightly smaller than our ESZXMM sample. The higher mass threshold for the two other samples enforce the best possible match to our data. As for the observed profiles (see Sect. 6.2), the simulated cluster profiles were renormalised by the factor ⟨ f(M) ⟩ . The average value for this factor is 1.0, 1.08 and 1.03 for the B04+N07+P08, Dolag (in prep.) and Battaglia et al. (2010) simulations, respectively. Finally, we accounted for the differences in the definition of R_{500} and M_{500} between observations and simulations. We corrected from the differences between the true mass and the hydrostatic masses which impact the estimation of both R_{500} and P_{500} (see Eq. (10)). We refer to Sect. 4.1 of A10 for a complete discussion. Assuming a hydrostatic bias r_{cor} = M_{true}/M_{HE} = 1.15 (Kay et al. 2004; Piffaretti & Valdarnini 2008), we divided the values of R_{500} and P_{500} derived from the simulations by and , respectively.
All three sets of data are reported in Fig. 6. It is well beyond the scope of this paper to extensively discuss the comparison between theoretical predictions. Nevertheless, we note that they agree within their respective dispersions across the whole radial range. The Dolag (in prep.) and Battaglia et al. (2010, 2012) profiles best agree within the central part, and are flatter than the B04+N07+P08 profile. This is likely due to the implementation of AGN feedback, which triggers energy injection at the clusters centre, balancing radiative cooling and thus stopping the gas cooling. In the outer parts where cooling is negligible, the B04+N07+P08 and Dolag (in prep.) profiles are in perfect agreement. The Battaglia et al. (2010) profile is slightly higher, but still compatible within its dispersion with the two other sets. Here again it might be due to the specific implementation of the simulatons.
Fig. 6 Measured pressure profiles from Planck and XMMNewton, displayed as in Fig. 4. The three shaded areas mark the dispersions about the average profiles for three samples of simulated clusters: in blue, the simulations from Borgani et al. (2004), Nagai et al. (2007a) and Piffaretti & Valdarnini (2008), which were used in A10 to derive the universal pressure profile together with the REXCESS data; in green, the simulated sample of clusters from Dolag (in prep.); and in brown the simulated clusters from Battaglia et al. (2012). The corresponding average profiles for each set of simulations are plotted as solid lines using the same colour scheme. In the bottom panel, we present the χ profiles between the observed profile and the simulated average profiles, taking into account their associated dispersion. See Sect. 6.2 for more details. 
Fig. 7 Gas mass fraction profile derived from the combined Planck and XMMNewton pressure profile, assuming for the temperature profile: (H1) the average best fit model across the sample from Xray spectroscopy (red line and striped area); or (H2) the same but extrapolating beyond R_{500} to a constant value equal to the average temperature measured in the last radial bin across the sample (blue line and striped area). The green dashed curve marks the expected gas fraction profile assuming the A10 pressure profile and kT(r) as in hypothesis (1). The star gives f_{gas}(r = R_{500}) for REXCESS clusters with M_{500} > 5 × 10^{14} M_{⊙} (Pratt et al. 2009). In maroon we reproduce the gas mass fraction profile derived from Suzaku measurements for the Perseus cluster (Simionescu et al. 2011). The solid and dashed black lines mark the cosmic baryon fraction expected from CMB measurements (Komatsu et al. 2011) and the expected gas fraction, assuming that 12% of baryons are in stars, respectively. The shaded blue and red areas translate the dispersion in the pressure profile across the ESZXMM sample as shown in Fig. 4, assuming hypothesis (H1) for kT(r). See Sect. 7.3. 
7.2.2. Comparison with the Planck pressure profile
Figure 6 compares our SZ+Xray observational constraints to the theoretical predictions described above. Overall, our observed pressure profile lies within the scatter across the profiles from these various sets of simulated clusters.
In the central parts, the data points lie on the lower envelope of the scatter of the B04+N07+P08 simulations, similarly to the REXCESS sample profiles (left panel of Figs. 4 and 7 in A10). This flatter inner slope is more compatible with that of the clusters from simulations which include the effect of AGN feedback (Dolag, in prep.; Battaglia et al. 2012), although our measured points lie above them.
With increasing radius (i.e., R > R_{500}) both SZ and Xray profiles are marginally compatible with the stacked profiles for simulated clusters. Interestingly our profile lies above the B04+N07+P08 and Dolag (in prep.) simulated clusters, and below the Battaglia et al. (2012) clusters. In the outer parts of halos, the predicted pressure in numerical simulations is essentially sustained by gravitational heating. The general agreement (within the dispersion) with our measured pressure profile favours the equilibrium between ions and electrons in the outskirts, i.e., both populations have an equivalent temperature. Indeed, a difference would result in a drop in thermal pressure (Rudd & Nagai 2009). This may suggest that we have the correct global understanding of the intracluster thermal pressure support. However, the discrepancies between simulations and the average observed profile in the clusters’ outer parts calls for a deeper investigation of the gas distribution and structure in clusters. Thus we need more detailed modelling of the baryon physics in cluster outskirts, where incomplete virialisation of the intracluster medium affects the thermal pressure. Even if a universal behaviour of cluster quantities is anticipated from the simplest model of gravitational collapse, dispersion is also expected, even without considering the effects of nongravitational physics, as shown by numerical simulations of halo formation (e.g., Navarro et al. 1997; Reed et al. 2011; Bhattacharya et al. 2011; Gao et al. 2012). This is plainly reflected in the large level of dispersion that we observe across our sample of simulated clusters.
On the observational side, our investigation of the subsample of CC and nonCC clusters did not show differences in the outer shape of the pressure profile (CC and nonCC cluster profiles are compatible within the measured statistical error bars). This might indicate that the aforementioned dispersion is not mainly driven by the dynamical state of clusters, which likely affects primarily the inner parts of halos. If our Planck constraints over the whole ESZXMM sample give a flatter outer shape with respect to the average simulated pressure profiles by B04+N07+P08 and Dolag (in prep.), we need nonetheless to keep in mind that our current sample may be affected by selection biases, as discussed in Planck Collaboration (2011f).
7.3. Constraint on the gas mass fraction
As already mentioned, Xray measurements hardly reach density contrasts of ~ 200 (e.g., George et al. 2009; Reiprich et al. 2009; Urban et al. 2011; Simionescu et al. 2011), so modelling beyond R_{200} only relies on predictions from numerical simulations. We used our Planck and XMMNewton derived pressure profile to investigate the average gas mass fraction distribution across our sample. Assuming the ICM to be a perfect gas, the thermal pressure is the genuine product of the temperature and the density, i.e., P(r) ∝ n_{e}(r) × kT(r). Following previous Xray works (Pratt et al. 2007; Arnaud et al. 2010; Démoclès et al. 2010), we used the best fit to the spectroscopic temperature profile (Démocles et al., in prep.) with the analytical function proposed by Vikhlinin et al. (2006) for each of our clusters. Beyond the reach of the Xray observations (i.e., [0.5–1.2] R_{500}), we have worked with two hypotheses for the temperature: (H1) the extrapolation of the best fit model (note tha ); or (H2) a constant temperature fixed to the last radial bin value. We have then derived the density profiles, and integrated them over the cluster volume to obtain the gas mass radial distribution. Meanwhile, we modelled each total mass profile with an NFW profile with the mass, M_{500}, derived from Eq. (9), and the concentration estimated from the c − M relation by Bhattacharya et al. (2011). We checked that the average of the individual NFW mass models agrees with the NFW model derived from the average values of M_{500} and c_{500} across the sample. The average gas fraction profile is computed as the ratio of the average gas mass model and the total mass model profiles.
Figure 7 shows the resulting gas mass fraction distribution for the two hypotheses on the temperature at large radii (red and blue solid lines). The red and blue striped areas picture the statistical errors measured from the SZ+Xray pressure profile. The red and blue shaded areas overplot the XMMNewton and Planck dispersion onto the pressure profile (see Fig. 4) using hypothesis (H1) for the temperature, and provide an illustration of the clustertocluster variation of f_{gas} within our sample.
From the stacked extrapolated temperature models, the average temperature is 1.3 keV at 3 R_{500}, which leads to f_{gas} ~ 0.2. An underestimation of this temperature by a third (i.e., kT ~ 1 keV) would increase the gas fraction likewise (i.e., f_{gas} ~ 0.3). For the lower bound, hypothesis (H2) fixes the temperature of the averaged profile beyond R_{500} to ~ 4 keV, which likely overestimates the true average value of the temperature across our sample in the (2–3) R_{500} radial range. Therefore we are confident that both hypotheses fairly bracket the range for the gas fraction distribution in massive clusters out to 3 R_{500}. In other words, our constraints on the pressure profile imply that the temperature required to flatten the gas fraction profile at 3 R_{500} to the expected cosmic value lies between 1 and 4 keV and is likely closer to the lower value.With hypothesis (H1) for kT(r), the green dashed line gives the gas fraction for the A10 pressure profile. As for our SZ+Xray pressure profile, it leads to f_{gas} values above the expected gas fraction out to 3 R_{500}.
At R_{500}, the measurement of f_{gas} for massive clusters presented by Pratt et al. (2009) (see also Giodini et al. 2009, for massive clusters with (4–10) × 10^{14} M_{⊙}) is lower by a factor 1.3 with respect to our estimate (although compatible within ± 2σ). We can also compare our result to the gas fraction profile reported by Simionescu et al. (2011) using Suzaku measurements in the direction of the Perseus cluster. The derived gas fraction profile is in very good agreement with ours within R_{500}. Beyond this radius, it rises within our plausible dispersion envelope (red and blue shaded areas in Fig. 7). Note however, that such a time consuming Xray observation sampled less than 5% of the sky area covered by the cluster and as a specific measurement is neither representative nor out of the allowed range for the gas fraction in cluster outskirts.
The constrained interval for the gas fraction from our observed pressure profile encompasses the cosmic gas fraction derived from CMB measurements (solid black line in Fig. 7, from WMAP7 results by Komatsu et al. 2011) and the expected hot gas fraction (dashed black line) assuming that 12% of cluster baryons are in stars (Gonzalez et al. 2007; Giodini et al. 2009).
8. Conclusion
From the Planck nominal mission (i.e., 14 months of survey), making use of its full spectral coverage of the SZ spectrum, we have extracted and reconstructed the SZ signal distribution in clusters of galaxies for a sample of 62 massive nearby clusters. All are individually detected in the survey with high significance (Planck Collaboration 2011c) and were previously used to investigate the total integrated SZ flux and the SZ scaling relations (Planck Collaboration 2011f). We have scaled and averaged together all the SZ profiles in the sample into a stacked profile. We have statistically detected the SZ signal out to 3 × R_{500}, providing, for the first time, stringent observational constraints on the ICM gas beyond R_{500} and out to a density contrast of δ ~ 50–100.
From the Planck SZ raw profiles, we have derived the underlying thermal pressure profiles of the ICM. Together with the pressure profile derived from the XMMNewton data, we have provided for the first time a comprehensive observational view of the distribution of thermal pressure support in clusters from 0.01 out 3 × R_{500}. We have fitted these unique measurements to a generalised NFW profile. Our best analytical representation over this wide radial range is given by the parameters [P_{0},c_{500},γ,α,β ] = [6.41,1.81,0.31,1.33,4.13 ] .
Our observational measurements further confirm the agreement of the SZ and Xray constraints on the intracluster gas properties within the inner part of the clusters commonly considered to be virialised, i.e., R_{500}, as found in Planck Collaboration (2011e) and Planck Collaboration (2011f). Overall it also agrees with a wide range of simulated clusters through the whole radial range, although in the central regions it matches best the numerical predictions that implement prescriptions for AGN feedback. The statistical nature of our stacked detection at large radii provides the average trend for the thermal pressure distribution in cluster outskirts out to 3 × R_{500}, which is slightly flatter than most theoretical predictions.
In conjunction with Xray constraints on the temperature profile of our clusters, we have derived the profile of the gas mass fraction out to the cluster outskirts. From reasonable hypotheses on the gas temperature to account for the lack of constraints beyond R_{500}, we have bracketed a range for the gas fraction in the cluster outer regions, which is compatible with the cosmic baryon fraction and the expected gas fraction in halos.
The processes governing the thermodynamical state of the outer regions in clusters still need to be understood from the theoretical and observational points of view (see Kravtsov & Borgani 2012, for a review and references). Issues such as gas clumping will affect the pressure estimation (i.e., effect on the Xray surface brightness, see Roncarelli et al. 2006), departures from hydrostatic equilibrium and contribution from nonthermal pressure (due to magnetic fields and/or cosmic rays) will modify the gas fraction. In this regard SZ observations provide a straightforward description of the thermal pressure distribution in massive halos, and a clear path to the gas fraction determination. Future SZ instruments with increased spatial resolution and sensitivity, and retaining the ability to map clusters of galaxies out to large radii will certainly provide us with further details and insight. In the meantime, with Planck we bring unique observational constraints that are extremely valuable to further test and understand the physics at play in the outskirts of clusters.
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.
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 and JA (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 DEISA (EU). We acknowledge the use of the Healpix software.
References
 Afshordi, N., Lin, Y.T., Nagai, D., & Sanderson, A. J. R. 2007, MNRAS, 378, 293 [NASA ADS] [CrossRef] [Google Scholar]
 Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 [NASA ADS] [CrossRef] [Google Scholar]
 Andersson, K., Benson, B. A., Ade, P. A. R., et al. 2011, ApJ, 738, 48 [NASA ADS] [CrossRef] [Google Scholar]
 Arnaud, M., Neumann, D. M., Aghanim, N., et al. 2001, A&A, 365, L80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 441, 893 [NASA ADS] [CrossRef] [EDP Sciences] [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]
 AtrioBarandela, F., Kashlinsky, A., Kocevski, D., & Ebeling, H. 2008, ApJ, 675, L57 [NASA ADS] [CrossRef] [Google Scholar]
 Basu, K., Zhang, Y.Y., Sommer, M. W., et al. 2010, A&A, 519, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Battaglia, N., Bond, J. R., Pfrommer, C., Sievers, J. L., & Sijacki, D. 2010, ApJ, 725, 91 [NASA ADS] [CrossRef] [Google Scholar]
 Battaglia, N., Bond, J. R., Pfrommer, C., & Sievers, J. L. 2012, ApJ, 758, 75 [NASA ADS] [CrossRef] [Google Scholar]
 Benson, B. A., Church, S. E., Ade, P. A. R., et al. 2004, ApJ, 617, 829 [NASA ADS] [CrossRef] [Google Scholar]
 Bersanelli, M., Mandolesi, N., Butler, R. C., et al. 2010, A&A, 520, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bertschinger, E. 1985, ApJS, 58, 39 [NASA ADS] [CrossRef] [Google Scholar]
 Bhattacharya, S., Habib, S., & Heitmann, K. 2011 [arXiv:1112.5479] [Google Scholar]
 Bobin, J., Moudden, Y., Starck, J.L., Fadili, J., & Aghanim, N. 2008, Stat. Meth., 5, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Böhringer, H., Schuecker, P., Pratt, G. W., et al. 2007, A&A, 469, 363 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bonamente, M., Joy, M. K., LaRoque, S. J., et al. 2006, ApJ, 647, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Bonamente, M., Hasler, N., Bulbul, E., et al. 2012, New J. Phys., 14, 025010 [NASA ADS] [CrossRef] [Google Scholar]
 Borgani, S., & Kravtsov, A. 2011, Adv. Sci. Lett., 4, 204 [CrossRef] [Google Scholar]
 Borgani, S., Murante, G., Springel, V., et al. 2004, MNRAS, 348, 1078 [NASA ADS] [CrossRef] [Google Scholar]
 Borgani, S., Finoguenov, A., Kay, S. T., et al. 2005, MNRAS, 361, 233 [NASA ADS] [CrossRef] [Google Scholar]
 Bulbul, G. E., Hasler, N., Bonamente, M., & Joy, M. 2010, ApJ, 720, 1038 [NASA ADS] [CrossRef] [Google Scholar]
 Carlstrom, J. E., Ade, P. A. R., Aird, K. A., et al. 2011, PASP, 123, 568 [NASA ADS] [CrossRef] [Google Scholar]
 Carvalho, P., Rocha, G., & Hobson, M. P. 2009, MNRAS, 393, 681 [NASA ADS] [CrossRef] [Google Scholar]
 Carvalho, P., Rocha, G., Hobson, M. P., & Lasenby, A. 2011 [arXiv:1112.4886] [Google Scholar]
 Challinor, A., & Lasenby, A. 1998, ApJ, 499, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Croston, J. H., Arnaud, M., Pointecouteau, E., & Pratt, G. W. 2006, A&A, 459, 1007 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Croston, J. H., Pratt, G. W., Böhringer, H., et al. 2008, A&A, 487, 431 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Delabrouille, J., Cardoso, J.F., Le Jeune, M., et al. 2009, A&A, 493, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Delabrouille, J., Betoule, M., Melin, J.B., et al. 2012, A&A, submitted [arXiv:1207.3675] [Google Scholar]
 Démoclès, J., Pratt, G. W., Pierini, D., et al. 2010, A&A, 517, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Eriksen, H. K., Banday, A. J., Górski, K. M., & Lilje, P. B. 2004, ApJ, 612, 633 [NASA ADS] [CrossRef] [Google Scholar]
 Evrard, A. E., Metzler, C. A., & Navarro, J. F. 1996, ApJ, 469, 494 [NASA ADS] [CrossRef] [Google Scholar]
 Gao, L., Frenk, C. S., Jenkins, A., Springel, V., & White, S. D. M. 2012, MNRAS, 419, 1721 [NASA ADS] [CrossRef] [Google Scholar]
 George, M. R., Fabian, A. C., Sanders, J. S., Young, A. J., & Russell, H. R. 2009, MNRAS, 395, 657 [NASA ADS] [CrossRef] [Google Scholar]
 Giodini, S., Pierini, D., Finoguenov, A., et al. 2009, ApJ, 703, 982 [NASA ADS] [CrossRef] [Google Scholar]
 Gonzalez, A. H., Zaritsky, D., & Zabludoff, A. I. 2007, ApJ, 666, 147 [NASA ADS] [CrossRef] [Google Scholar]
 Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Halverson, N. W., Lanting, T., Ade, P. A. R., et al. 2009, ApJ, 701, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Hurier, G., Hildebrandt, S. R., & MaciasPerez, J. F. 2010 [arXiv:1007.1149] [Google Scholar]
 Jia, S. M., Böhringer, H., Pointecouteau, E., Chen, Y., & Zhang, Y. Y. 2008, A&A, 489, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kaiser, N., Squires, G., & Broadhurst, T. 1995, ApJ, 449, 460 [NASA ADS] [CrossRef] [Google Scholar]
 Kay, S. T., Thomas, P. A., Jenkins, A., & Pearce, F. R. 2004, MNRAS, 355, 1091 [NASA ADS] [CrossRef] [Google Scholar]
 Kitayama, T., Komatsu, E., Ota, N., et al. 2004, PASJ, 56, 17 [NASA ADS] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Kompaneets, A. S. 1957, Sov. Phys. – JETP Lett., 4, 730 [Google Scholar]
 Korngut, P. M., Dicker, S. R., Reese, E. D., et al. 2011, ApJ, 734, 10 [NASA ADS] [CrossRef] [Google Scholar]
 Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 [NASA ADS] [CrossRef] [Google Scholar]
 Kravtsov, A. V., & Borgani, S. 2012, ARA&A, 50, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Lamarre, J., Puget, J., Ade, P. A. R., et al. 2010, A&A, 520, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Leach, S. M., Cardoso, J.F., Baccigalupi, C., et al. 2008, A&A, 491, 597 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lieu, R., Mittaz, J. P. D., & Zhang, S.N. 2006, ApJ, 648, 176 [NASA ADS] [CrossRef] [Google Scholar]
 Mandolesi, N., Bersanelli, M., Butler, R. C., et al. 2010, A&A, 520, A3 [Google Scholar]
 Melin, J., Bartlett, J. G., & Delabrouille, J. 2006, A&A, 459, 341 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Melin, J.B., Bartlett, J. G., Delabrouille, J., et al. 2011, A&A, 525, A139 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mennella, A., Butler, R. C., Curto, A., et al. 2011, A&A, 536, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Muchovej, S., Mroczkowski, T., Carlstrom, J. E., et al. 2007, ApJ, 663, 708 [NASA ADS] [CrossRef] [Google Scholar]
 Nagai, D., Kravtsov, A. V., & Vikhlinin, A. 2007a, ApJ, 668, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Nagai, D., Vikhlinin, A., & Kravtsov, A. V. 2007b, ApJ, 655, 98 [NASA ADS] [CrossRef] [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1995, MNRAS, 275, 56 [Google Scholar]
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 [NASA ADS] [CrossRef] [Google Scholar]
 Peebles, P. J. E. 1980, The largescale structure of the universe (Princeton University Press) [Google Scholar]
 Piffaretti, R., & Valdarnini, R. 2008, A&A, 491, 71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Piffaretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., & Melin, J.B. 2011, A&A, 534, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Plagge, T., Benson, B. A., Ade, P. A. R., et al. 2010, ApJ, 716, 1118 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration 2011a, A&A, 536, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2011b, A&A, 536, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2011c, A&A, 536, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2011d, A&A, 536, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2011e, A&A, 536, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2011f, A&A, 536, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2011g, A&A, 536, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2012, A&A, 543, A102 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2013a, A&A, 550, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2013b, A&A, 550, A130 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2013c, A&A, in press, DOI: 10.1051/00046361/201220247 [Google Scholar]
 Planck HFI Core Team 2011a, A&A, 536, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck HFI Core Team 2011b, A&A, 536, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pointecouteau, E., Giard, M., & Barret, D. 1998, A&A, 336, 44 [NASA ADS] [Google Scholar]
 Pointecouteau, E., Giard, M., Benoit, A., et al. 2001, ApJ, 552, 42 [NASA ADS] [CrossRef] [Google Scholar]
 Pointecouteau, E., Hattori, M., Neumann, D., et al. 2002, A&A, 387, 56 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pointecouteau, E., Arnaud, M., & Pratt, G. W. 2005, A&A, 435, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Postman, M., Coe, D., Benítez, N., et al. 2012, ApJS, 199, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Pratt, G. W., Böhringer, H., Croston, J. H., et al. 2007, A&A, 461, 71 [NASA ADS] [CrossRef] [EDP Sciences] [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]
 Pratt, G. W., Arnaud, M., Piffaretti, R., et al. 2010, A&A, 511, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Reed, D. S., Koushiappas, S. M., & Gao, L. 2011, MNRAS, 415, 3177 [NASA ADS] [CrossRef] [Google Scholar]
 Reiprich, T. H., Hudson, D. S., Zhang, Y.Y., et al. 2009, A&A, 501, 899 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Remazeilles, M., Delabrouille, J., & Cardoso, J.F. 2011, MNRAS, 418, 467 [NASA ADS] [CrossRef] [Google Scholar]
 Rines, K., & Diaferio, A. 2006, AJ, 132, 1275 [NASA ADS] [CrossRef] [Google Scholar]
 Roncarelli, M., Ettori, S., Dolag, K., et al. 2006, MNRAS, 373, 1339 [NASA ADS] [CrossRef] [Google Scholar]
 Rudd, D. H., & Nagai, D. 2009, ApJ, 701, L16 [NASA ADS] [CrossRef] [Google Scholar]
 Sazonov, S. Y., & Sunyaev, R. A. 1998, ApJ, 508, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Sehgal, N., Trac, H., Acquaviva, V., et al. 2011, ApJ, 732, 44 [NASA ADS] [CrossRef] [Google Scholar]
 Sifon, C., Menanteau, F., Hasselfield, M., et al. 2012, ApJ, submitted [arXiv:1201.0991] [Google Scholar]
 Simionescu, A., Allen, S. W., Mantz, A., et al. 2011, Science, 331, 1576 [NASA ADS] [CrossRef] [Google Scholar]
 Sun, M. 2012, New J. Phys., 14, 045004 [NASA ADS] [CrossRef] [Google Scholar]
 Sun, M., Sehgal, N., Voit, G. M., et al. 2011, ApJ, 727, L49 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1970, Comm. Astrophys. Space Phys., 2, 66 [Google Scholar]
 Sunyaev, R. A., & Zeldovich, Y. B. 1972, Comm. Astrophys. Space Phys., 4, 173 [Google Scholar]
 Swetz, D. S., Ade, P. A. R., Amiri, M., et al. 2011, ApJS, 194, 41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tauber, J. A., Mandolesi, N., Puget, J., et al. 2010, A&A, 520, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tozzi, P., & Norman, C. 2001, ApJ, 546, 63 [NASA ADS] [CrossRef] [Google Scholar]
 Urban, O., Werner, N., Simionescu, A., Allen, S. W., & Böhringer, H. 2011, MNRAS, 414, 2101 [NASA ADS] [CrossRef] [Google Scholar]
 Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 691 [NASA ADS] [CrossRef] [Google Scholar]
 Voit, G. M. 2005, Rev. Mod. Phys., 77, 207 [NASA ADS] [CrossRef] [Google Scholar]
 Voit, G. M., & Ponman, T. J. 2003, ApJ, 594, L75 [NASA ADS] [CrossRef] [Google Scholar]
 Voit, G. M., Bryan, G. L., Balogh, M. L., & Bower, R. G. 2002, ApJ, 576, 601 [NASA ADS] [CrossRef] [Google Scholar]
 Voit, G. M., Balogh, M. L., Bower, R. G., Lacey, C. G., & Bryan, G. L. 2003, ApJ, 593, 272 [NASA ADS] [CrossRef] [Google Scholar]
 Walker, S. A., Fabian, A. C., Sanders, J. S., George, M. R., & Tawara, Y. 2012, MNRAS, 2785 [Google Scholar]
 Whitbourn, J. R., Shanks, T., & Sawangwit, U. 2011, MNRAS, submitted [arXiv:1107.2654] [Google Scholar]
 Wojtak, R., & Lokas, E. L. 2010, MNRAS, 408, 2442 [NASA ADS] [CrossRef] [Google Scholar]
 Zacchei, A., Maino, D., Baccigalupi, C., et al. 2011, A&A, 536, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
Online material
Appendix A: SZ map reconstruction methods
In this Appendix we discuss the different SZ map reconstruction methods applied to the Planck data in Sect. 4.1.

1.
MILCA (Hurier et al. 2010):The thermal SZ signal reconstruction is performed on thesix Planck allsky maps from 100 GHzto 857 GHz. MILCA (Modified InternalLinear Combination Algorithm) is a component separation approach aiming at extracting a chosen component(here the thermal SZ signal) from a multichannel set ofinput maps. It is mainly based on the well known InternalLinear Combination approach (e.g., Eriksenet al. 2004), that searches for thelinear combination of the input maps that minimises the varianceof the final reconstructed map by imposing spectral constraints.In this paper, we applied MILCA using two constraints: preservation of the thermal SZ signal (knowing the SZ spectral signature);and removal of the CMB contamination in the final SZ map (alsomaking use of the well known spectrum of the CMB). In addition,to compute the weights of the linear combination, we have usedthe extra degrees of freedom to minimise residuals from othercomponents (2 degrees) and from the noise (2 degrees). The noisecovariance matrix was estimated from the frequency error maps(see Sect. 2).

2.
NILC: Needlet ILC performs a linear combination of the observed maps which has minimum variance under the constraint of offering unit response to the component of interest (here the thermal SZ, whose frequency scaling is known). The weights of the ILC depend on the covariance between the various observations, and can be computed, for example, on domains of the observed pixels (pixel space) or angular scales (spherical harmonics). In the case of NILC, covariances (and hence weights for component separation) are computed independently in domains of a needlet decomposition (spherical wavelet frame). The needlet decomposition provides localisation of the ILC filters both in pixel and in multipole space, allowing us to deal with local contamination conditions varying both in position and in scale. NILC was developed to extract a CMB map from WMAP data (Delabrouille et al. 2009) and was also tested for SZ effect extraction in Leach et al. (2008). Multicomponent extensions have been investigated by Remazeilles et al. (2011).

3.
GMCA (Bobin et al. 2008) is a blind source separation method developed for separating sources from instantaneous linear mixtures. The components are assumed to be sparsely represented (i.e., have a few significant samples in a specific basis) in a socalled sparse representation Ψ (typically wavelets). The assumption that the components have a sparse representation in the wavelet domain is equivalent to assuming that most components have a certain spatial regularity. These components and their spectral signatures are then recovered by minimising the number of significant coefficients in Ψ. Recently, LGMCA has been further introduced to analyse the CMB data in a local and multiscale manner (Bobin et al. 2008). More precisely, the multichannel data are analysed in four frequency bands in spherical harmonics (i.e., wavelet bands). In each wavelet band, GMCA is applied locally on small patches with banddependent sizes. In each band, the observations are analysed at the same banddependent resolution. In Bobin et al. (2008), it was shown this local analysis approach enhances the separation quality. The spectral signatures of CMB, freefree and SZ are assumed to be known.
Fig. A.1 Comparison between the SZ profile reconstruction from MILCA, NILC and GMCA The comparison is presented in the form of the averaged stacked radial profile from our 62 clusters of the ESZXMM sample, individually scaled respectively in the x and yaxis directions according to R_{500} and Φ (Eq. (12)). The lightred shaded area marks the dispersion about the average stacked profile for the MILCA method. The points within each individual profile are correlated at about the ~ 20% level (see Sect. 4.2) and the plotted errors correspond to the square root of the diagonal elements of the covariance matrix of the profile. The solid black line (labelled in the legend as “Expected from P_{A10}”) is the average stacked profile obtained from the expected SZ individual profiles drawn from the universal pressure profile by A10, parameterised according R_{500} and Y_{500} derived from the XMMNewton data analysis presented in Planck Collaboration (2011f). The two dotted lines indicate the associated dispersion across the sample. The bottom panel shows the value of χ at each point of the observed profiles with respect to the expectations from the universal profile of A10. 
Appendix B: Validation on simulations
B.1. Profiles
We have validated our reconstruction of the SZ signal and profile computation methods with simulations. We used the Planck sky model (PSM, Delabrouille 2012) to simulate the whole sky as seen by Planck. The PSM includes most astrophysical components acting as foreground or background contamination with respect to the SZ signal (i.e., Galactic dust, freefree and synchrotron emission, extragalactic sources in the infrared and radio parts of the spectrum and CMB). For the instrumental noise we added to the simulation the allsky error map drawn from the Planck jackknife maps at each frequency (see Sect. 2), thus mimicking the noise properties of the actual survey (noise and systematic effects). For the SZ component, we built a full sky template for all the 1743 clusters from the MCXC metacatalogue (Piffaretti et al. 2011), which includes all 62 clusters in our sample. Each cluster was modelled assuming spherical symmetry and a thermal pressure distribution following the universal profile from A10. We recall that the universal pressure profile provided by these authors is the best fit of a generalised NFW profile (Nagai et al. 2007a) to the median profile derived from the REXCESS pressure profile and the predictions from three different sets of numerical simulations (i.e., from Borgani et al. 2004; Nagai et al. 2007a,b; Piffaretti & Valdarnini 2008). In our case, the universal pressure profile is parametrised using the values of R_{500} and M_{500} from the MCXC, together with the M–Y_{X} scaling relation from Arnaud et al. (2007, 2010) to estimate Y_{500}.
We processed the simulated Planck sky maps through each of our SZ map reconstruction methods. For our sample of 62 clusters, we extracted the patches and computed the individual profiles as described in Sect. 4.2. We stacked the simulated profiles for the three methods after renormalisation, making use of the aforementioned R_{500} and Y_{500} values from the MCXC prescription.
The resulting stacked profiles are displayed in Fig. A.1. All three tested methods agree remarkably well over the whole radial range. Taking into account the correlated errors of each profile, the χ^{2} over the radial range (0–10) × R_{500} for NILC and GMCA with respect to MILCA are 0.78 and 1.54, respectively. The reconstructed profiles are also in good agreement with the reconstructed stacked model. Accounting for the errors on the model and the correlated errors between points on the reconstructed profiles, the reduced χ^{2} for MILCA, NILC and GMCA with respect to the stacked model profile are 0.80, 0.66 and 1.03 within 3 × R_{500}. This excellent agreement also translates into a small relative error below 15% in the same range of radii. Beyond R ~ 2 R_{500}, the three reconstructed profiles drop slightly below the input model, even though they still agree within their dispersions and errors. This drop is caused by an intrinsic overprescription of point sources in the PSM simulations with respect to real data. Despite our careful masking, the unidentified or unresolved extra sources in the PSM affect the power spectrum of the noise, increasing the background level and making it more difficult to detect the SZ signal further out in the simulations than in real data. We stress that this is an effect limited to our test simulations, and not relevant for the real data.
B.2. Consistency between Y_{500} measurements
Fig. B.1 Comparison of Y_{500} values from the Planck cluster detection method (MMF3) against our best fit for each of our 62 clusters (from NILC in this case). In both cases the distribution of thermal pressure is assumed to follow a universal pressure profile (A10) for which the scaled radius R_{500} is fixed to the XMMNewton best fit value. See Sect. 6.3 for details. 
To further validate the above, and as a consistency check with previous Planck results, we fitted each individual SZ profile with a projected, PSFconvolved universal pressure profile. We fixed R_{500} to the best fitting Xray value from Planck Collaboration (2011f) and only fitted the normalisation, Y_{500}. For a given profile, P, with an associated covariance matrix, C, and a model M, the χ^{2} statistic can be expressed as (B.1)The solution that minimises χ^{2}, Y_{500}, and its associated uncertainty, σ_{Y500}, are analytic: We compared the resulting value of Y_{500} to those obtained from the detection algorithms PwS and MMF3. Both detection algorithms were run under the same assumptions, i.e., universal pressure profile with the “nonstandard” parameterisation and a fixed size (R_{500}). The fluxes were derived from the Planck nominal mission survey maps. These fluxes are in agreement with those extracted from the Planck first year survey (Planck Collaboration 2011f). For instance, for MMF3, the average ratio is 1.04 ± 0.09, compatible within the average relative errors with the fluxes extracted from the nominal mission.
The relation between the integrated Comptonisation parameters derived from the SZ profiles and from the detection methods is displayed for MMF3 vs NILC fuxes in Fig. B.1. A BCES linear fit accounting for the errors in both x and yaxis directions (Akritas & Bershady 1996) yields a best fitting slope and normalisation of 1.00 ± 0.02 and 0.01 ± 0.01 (), respectively. There is 4% dispersion about this relation. Furthermore, the median and associated deviation for the onetoone ratios of PwS and MMF3 fluxes compared to the fitted NILC profile values are 1.03 ± 0.11 and 0.98 ± 0.08, respectively. The median ratios for the errors are 1.10 ± 0.36 are 1.64 ± 0.34. The errors associated with our best fit values are in very good agreement with those from PwS; however, they are smaller than those derived from the MMF3 method. This difference with MMF3 might be explained by the different treatment of the noise. Whilst we estimated the noise from the surroundings of the clusters, defined according to the scaled radius of each target, and PwS also used a local computation of the noise, the MMF3 algorithm used a larger area to characterise the noise properties. Our tests against the PwS and MMF3 fluxes emphasise our grasp on the characterisation of the local noise for each target, and further assess the reliability of our estimates of the statistical errors in our stacked profiles.
Such agreement between the fluxes derived from our fit over the cluster profiles and those from the detection algorithm, strongly confirms the reliability of our SZ flux estimate and further strengthen the overall result stated in the Planck collaboration early papers that the SZ and the Xray fluxes within R_{500} are fully consistent (Planck Collaboration 2011c,d,e,f).
Appendix C: Pressure profiles of individual clusters
In this appendix we present the individual pressure profiles obtained from Planck and XMMNewton data for each of the 62 clusters in our sample (see Fig. C.1). We show the agreement with the universal pressure profile (A10) and provide the best fitting model for each cluster to the GNFW profile in Table C.1. The following hard limit priors were adopted during the fitting procedure: 0 < P_{0} < 100; 0 < c_{500} < 10; 0 < α < 10; and 0 < β < 15.
Fig. C.1 Measured individual pressure profiles for the ESZXMM sample using the Planck and XMMNewton data. As for the stacked pressure profile, the maximum difference pointtopoint between the three SZ reconstruction methods have been taken into account in the Planck measurement error budget (Sect. 6.3). The best fit model on each individual profile is shown as a blue solid line. The black and green solid curves mark the best fit model to the stacked SZ+Xray pressure profile (Sect. 6.2) and the A10 best fit model, respectively. 
Fig. C.1 continued. 
Fig. C.1 continued. 
Fig. C.1 continued. 
Best fit parameters for a generalised NFW pressure profile for individual clusters.
All Tables
Best fit parameters for a generalised NFW pressure profile, with different numbers of fixed parameters.
Best fit parameters for a generalised NFW pressure profile for individual clusters.
All Figures
Fig. 1 Average scaled SZ profiles for each of the nine Planck frequency bands (decreasing frequency from left to right and from top to bottom). The points within each individual profile are correlated at about the 20% level (see Sect. 4.3.1). At each frequency, the stacked radial profile is obtained from the average of the 62 individual cluster profiles scaled by R_{500} and Φ_{i}/ ⟨ Φ_{i} ⟩ , in the x and yaxis directions, respectively (see Sect. 5.1). The lightred shaded area in each panel indicates the statistical uncertainty about the average observed profile. The solid black line corresponds to the prediction of the universal pressure profile (A10) for our cluster sample in each band, and the two dashed lines indicate its associated dispersion. 

In the text 
Fig. 2 Stacked scaled image (size is 10 R_{500} on a side) using a logarithmic stretch in y. Individual maps are rescaled by Φ_{i} before averaging, and then multiplied by ⟨ Φ_{i} ⟩ . The black and white circles mark the loci of 1 and 3 × R_{500}, respectively. 

In the text 
Fig. 3 Left: individual profiles (grey lines) and scaled stacked profile (red points) for our sample of 62 clusters. The lightred shaded area marks the dispersion about the average stacked profile, with its upper and lower limits highlighted by the two dashed red lines. The points within each individual profile are correlated at about the 20% level (see Sect. 4.3.1). The error bars on the Planck data points are purely statistical and correspond to the square root of the diagonal elements of the covariance matrix. The solid black line (labeled in the legend as “Expected from P_{A10}”) is the average stacked profile obtained from the expected SZ individual profiles drawn from the universal pressure profile of A10, parameterised according R_{500} and Y_{500} derived from the XMMNewton data analysis presented in Planck Collaboration (2011f). The two dotted black lines indicate the associated dispersion about this model profile. The bottom panel shows the value of χ at each point with respect to the expectations from the universal profile taking into account the measured errors. Right: matrix of the correlation coefficient for the observed radial stacked profiles. Both panels are illustrated here for the MILCA results. 

In the text 
Fig. 4 Left: Planck pressure profile obtained from the average of the individual pressure profiles across our sample (red points) shown together with the stacked pressure profile derived from the XMMNewton data for the same sample (purple points). The dispersions about the SZ and Xray profiles are depicted, respectively, by the red and purple shaded areas. Our best fit GNFW profile is shown as a solid black line and that of A10 as a solid green line. The lower panel shows the χ profile of these two best models taking into account the statistical errors and dispersion about the observed profile. Right: stacked profile from Planck and XMMNewton for the subsamples of coolcore (open circles) and non coolcore (open stars) clusters within the ESZXMM sample. The shaded areas are identical to the one shown in the left panel. Our best fit models for each subsample are shown as black solid and dashed lines for the coolcore and non coolcore clusters, respectively (see Table 1). The best A10 fit for coolcore and non coolcore clusters are shown as green solid and dashed lines, respectively. The lower panel gives the χ profiles of our best fit models with respect to the measured profiles and associated errors (including dispersion). On both panels the Planck data points are correlated at about the 20% level (see Sect. 4.3.1). As for the stacked SZ profile, the error bars on the Planck points are purely statistical and correspond to the square root of the diagonal elements of the covariance matrix. For the XMMNewton points they correspond to the statistical error on mean. 

In the text 
Fig. 5 Marginalised posterior likelihood distribution for our best MCMC fit to a GNFW pressure profile with free parameters [P_{0},c_{500},α,β] and γ = 0.31. The best fit values of the parameters are marked as white crosses and solid red lines, respectively, on the 2D and 1D posterior distributions. On the 2D distributions, the white triangles mark the value of the A10 parameters, while the red, purple and green solid lines define, respectively, the 68, 95 and 99% confidence levels. 

In the text 
Fig. 6 Measured pressure profiles from Planck and XMMNewton, displayed as in Fig. 4. The three shaded areas mark the dispersions about the average profiles for three samples of simulated clusters: in blue, the simulations from Borgani et al. (2004), Nagai et al. (2007a) and Piffaretti & Valdarnini (2008), which were used in A10 to derive the universal pressure profile together with the REXCESS data; in green, the simulated sample of clusters from Dolag (in prep.); and in brown the simulated clusters from Battaglia et al. (2012). The corresponding average profiles for each set of simulations are plotted as solid lines using the same colour scheme. In the bottom panel, we present the χ profiles between the observed profile and the simulated average profiles, taking into account their associated dispersion. See Sect. 6.2 for more details. 

In the text 
Fig. 7 Gas mass fraction profile derived from the combined Planck and XMMNewton pressure profile, assuming for the temperature profile: (H1) the average best fit model across the sample from Xray spectroscopy (red line and striped area); or (H2) the same but extrapolating beyond R_{500} to a constant value equal to the average temperature measured in the last radial bin across the sample (blue line and striped area). The green dashed curve marks the expected gas fraction profile assuming the A10 pressure profile and kT(r) as in hypothesis (1). The star gives f_{gas}(r = R_{500}) for REXCESS clusters with M_{500} > 5 × 10^{14} M_{⊙} (Pratt et al. 2009). In maroon we reproduce the gas mass fraction profile derived from Suzaku measurements for the Perseus cluster (Simionescu et al. 2011). The solid and dashed black lines mark the cosmic baryon fraction expected from CMB measurements (Komatsu et al. 2011) and the expected gas fraction, assuming that 12% of baryons are in stars, respectively. The shaded blue and red areas translate the dispersion in the pressure profile across the ESZXMM sample as shown in Fig. 4, assuming hypothesis (H1) for kT(r). See Sect. 7.3. 

In the text 
Fig. A.1 Comparison between the SZ profile reconstruction from MILCA, NILC and GMCA The comparison is presented in the form of the averaged stacked radial profile from our 62 clusters of the ESZXMM sample, individually scaled respectively in the x and yaxis directions according to R_{500} and Φ (Eq. (12)). The lightred shaded area marks the dispersion about the average stacked profile for the MILCA method. The points within each individual profile are correlated at about the ~ 20% level (see Sect. 4.2) and the plotted errors correspond to the square root of the diagonal elements of the covariance matrix of the profile. The solid black line (labelled in the legend as “Expected from P_{A10}”) is the average stacked profile obtained from the expected SZ individual profiles drawn from the universal pressure profile by A10, parameterised according R_{500} and Y_{500} derived from the XMMNewton data analysis presented in Planck Collaboration (2011f). The two dotted lines indicate the associated dispersion across the sample. The bottom panel shows the value of χ at each point of the observed profiles with respect to the expectations from the universal profile of A10. 

In the text 
Fig. B.1 Comparison of Y_{500} values from the Planck cluster detection method (MMF3) against our best fit for each of our 62 clusters (from NILC in this case). In both cases the distribution of thermal pressure is assumed to follow a universal pressure profile (A10) for which the scaled radius R_{500} is fixed to the XMMNewton best fit value. See Sect. 6.3 for details. 

In the text 
Fig. C.1 Measured individual pressure profiles for the ESZXMM sample using the Planck and XMMNewton data. As for the stacked pressure profile, the maximum difference pointtopoint between the three SZ reconstruction methods have been taken into account in the Planck measurement error budget (Sect. 6.3). The best fit model on each individual profile is shown as a blue solid line. The black and green solid curves mark the best fit model to the stacked SZ+Xray pressure profile (Sect. 6.2) and the A10 best fit model, respectively. 

In the text 
Fig. C.1 continued. 

In the text 
Fig. C.1 continued. 

In the text 
Fig. C.1 continued. 

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.