A Redshift estimates of confirmed candidates

Planck Intermediate Results. IV. The XMM-Newton validation programme for new Planck galaxy clusters

Key Words.:
Cosmology: observations Galaxies: clusters: general Galaxies: clusters: intracluster medium Cosmic background radiation, X-rays: galaxies: clusters

We present the final results from the XMM-Newton validation follow-up of new Planck galaxy cluster candidates. We observed 15 new candidates, detected with signal-to-noise ratios between 4.0 and 6.1 in the 15.5-month nominal Planck survey. The candidates were selected using ancillary data flags derived from the ROSAT All Sky Survey (RASS) and Digitized Sky Survey all-sky maps, with the aim of pushing into the low SZ flux, high- regime and testing RASS flags as indicators of candidate reliability. Fourteen new clusters were detected by XMM-Newton, ten single clusters and two double systems. Redshifts from X-ray spectroscopy lie in the range 0.2 to 0.9, with six clusters at . Estimated masses () range from to . We discuss our results in the context of the full XMM-Newton validation programme, in which 51 new clusters have been detected. This includes four double and two triple systems, some of which are chance projections on the sky of clusters at different redshifts. We find that association with a source from the RASS-Bright Source Catalogue is a robust indicator of the reliability of a candidate, whereas association with a source from the RASS-Faint Source Catalogue does not guarantee that the SZ candidate is a bona fide cluster. Nevertheless, most Planck clusters appear in RASS maps, with a significance greater than being a good indication that the candidate is a real cluster. Candidate validation from association with SDSS galaxy overdensity at is also discussed. The full sample gives a Planck sensitivity threshold of arcmin, with indication for Malmquist bias in the relation below this threshold. The corresponding mass threshold depends on redshift. Systems with at are easily detectable with Planck. The newly-detected clusters follow the relation derived from X-ray selected samples. Compared to X-ray selected clusters, the new SZ clusters have a lower X-ray luminosity on average for their mass. There is no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. In particular, there is no significant evolution of the ratio.

1 Introduction

The Planck1 satellite has been surveying the millimetre sky since 2009. Its two instruments together cover nine frequency bands: the Low Frequency Instrument (LFI; Mandolesi et al. 2010; Bersanelli et al. 2010; Mennella et al. 2011) at 30, 44, and 70 GHz, and the High Frequency Instrument (HFI; Lamarre et al. 2010; Planck HFI Core Team 2011) at 100, 143, 217, 353, 545, and 857 GHz. Before the HFI coolant ran out in January 2012, Planck had successfully performed nearly 5 surveys of the entire sky.

Planck allows the detection of galaxy clusters by their imprint on the cosmic microwave background (CMB) via the Sunyaev-Zeldovich (SZ) effect, a characteristic spectral distortion of the CMB due to inverse Compton scattering of photons by hot electrons in the intra-cluster medium (Sunyaev & Zeldovich 1972). The SZ signal of galaxy clusters is expected to correlate tightly with cluster mass (e.g., da Silva et al. 2004) and its surface brightness is independent of redshift. SZ selected cluster samples are thus particularly well-suited for statistical studies of the galaxy cluster population, either as a probe of the physics of structure formation, or for cosmological studies based on cluster abundance as a function of mass and redshift. Compared to other SZ surveys, such as those with the Atacama Cosmology Telescope (ACT, Marriage et al. 2011) or the South Pole Telescope (SPT, Carlstrom et al. 2009), the Planck survey covers an exceptionally large volume; indeed, it is the first all-sky survey since the ROSAT All-Sky Survey (RASS) in the X-ray domain. Planck allows the detection of clusters below the flux limit of RASS based catalogues at redshifts typically greater than (Planck Collaboration I 2012, Fig. 9). The first Planck SZ catalogue, the Early SZ (ESZ) sample, was published in Planck Collaboration VIII (2011). It contains 189 clusters and candidates detected at high signal-to-noise ratio () in the all-sky maps from the first ten months of observations, 20 of which were previously unknown. At the release of the ESZ sample, 12 of those 20 had been confirmed as new clusters, 11 using XMM-Newton validation observations undertaken in Director’s Discretionary Time (DDT) via an agreement between the XMM-Newton and Planck Project Scientists.

All cluster surveys include false detections. For Planck, these are mainly due to inhomogeneous, non-isotropic, and highly non-Gaussian fluctuations (galactic dust emission, confusion noise as result of the unsubtracted point sources, etc.) in the complex microwave astrophysical sky. After identification of known clusters, a follow-up programme is required for cluster confirmation and redshift estimation. It is essential to build as pure as possible an initial candidate sample in order for such a programme to be efficient and manageable. For this we rely both on internal Planck candidate selection and assessment of the SZ signal quality, and on cross-correlation with ancillary data, as described in Planck Collaboration VIII (2011). Beyond simple confirmation of new clusters, the XMM-Newton validation programme aims to refine this validation process and to yield a better understanding of the new objects that Planck is detecting. It consists of snapshot exposures ( ks), sufficient for unambiguous discrimination between clusters and false candidates (Planck Collaboration IX 2011), for a total allocated time of 500 ks for 50 candidates.

In the first two follow-up programmes, described by Planck Collaboration IX (2011), we observed 25 candidates in total and helped to define the selection criteria for the ESZ sample. They yielded the confirmation of 17 single clusters, two double systems, and two triple systems2. The observations showed that the new clusters are on average less X-ray-luminous and more morphologically disturbed than their X-ray-selected counterparts of similar mass, suggesting that Planck may be revealing a non-negligible population of massive, dynamically perturbed objects that are under-represented in X-ray surveys. However, despite their particular properties, the new clusters appear to follow the relation established for X-ray selected objects, where , introduced by Kravtsov et al. (2006), is the product of the gas mass and temperature.

In the third follow-up programme, described in Planck Collaboration I (2012), we observed 11 candidates with lower SZ detection levels () than the previous programmes () in order to investigate the internal SZ quality flags. Probing lower SZ flux than previous campaigns, the third programme also demonstrated the capability of Planck to find new clusters below the RASS limit and up to high , including the blind detection at of PLCK G266.627.3, confirmed by XMM-Newton to be an cluster at (Planck Collaboration XXVI 2011). We also detected tentative evidence for Malmquist bias in the relation, with a turnover at arcmin.

In the fourth and last XMM-Newton validation programme, presented here, we further probe the low SZ flux, high redshift regime. The sample includes 15 candidates, detected at signal-to-noise ratios between 4 and 6.1 in the 15.5-month nominal survey data. We use the results from all XMM-Newton validation observations to address the use of ancillary RASS information as an indicator of candidate reliability (Sec. 5). The evolution of cluster SZ/X-ray properties is discussed in Sec. 6. This paper, together with Planck Collaboration IX (2011); Planck Collaboration XXVI (2011); Planck Collaboration I (2012), presents our complete analysis of the DDT XMM-Newton validation programme.

We adopt a CDM cosmology with   km s Mpc, , and . The factor is the ratio of the Hubble constant at redshift to its present-day value. The quantities and are the total mass and radius corresponding to a total density contrast , as compared to , the critical density of the Universe at the cluster redshift; . The SZ flux is characterised by , where is the spherically integrated Compton parameter within , and is the angular-diameter distance to the cluster. Thus, as defined here, has units of solid angle and is given in arcmin in Table 2.

Name RA Dec S/N OBSID Filter Clean fraction Category Confirmed
(deg) (deg) (ks EPN) (EMOS/EPN)
PLCK G348.425.5 291.257 49.426 6.12 3 A 0679180101 t t t 10.6 1.0 0.9 PHZ Y
PLCK G329.522.7 278.270 65.570 5.84 3 B 0679181501 m m t 8.5 1.0 1.0 Y
PLCK G219.934.4 73.680 20.269 5.74 2 A 0679180501 t t t 9.5 1.0 0.9 PHZ Y
PLCK G352.124.0 290.233 45.842 5.63 2 C 0679180201 m m t 8.5 1.0 1.0 PHZ Y
PLCK G305.944.6 5.946 72.393 5.40 3 B 0679180301 t t t 10.1 0.4 0.2 Y
PLCK G196.745.5 55.759 8.704 5.21 3 B 0679180401 m m m 9.0 1.0 0.8 Y
PLCK G208.674.4 30.044 24.897 5.01 3 B 0679180601 t t t 2.9 0.8 0.8 Y
PLCK G130.117.0 22.678 45.288 4.93 3 C 0679180801 t t t 7.5 1.0 1.0 Y
PLCK G239.940.0 71.683 37.029 4.76 3 B 0679181001 t t t 9.9 1.0 1.0 Y
PLCK G310.527.1 201.148 35.245 4.77 2 B 0679180901 t t t 12.9 0.9 0.8 PHZ
PLCK G196.468.3 34.921 19.263 4.73 2 B 0679181101 t t t 11.8 0.7 0.4 PHZ
PLCK G204.715.9 113.614 14.295 4.57 3 A 0679180701 t t t 9.4 1.0 1.0 Y
PLCK G011.240.4 315.233 33.107 4.47 1 C 0679181201 t t t 8.5 1.0 1.0 PHZ Y
PLCK G147.316.6 44.099 40.291 4.41 3 B 0679181301 t t t 10.7 0.9 0.6 PHZ Y
PLCK G210.620.4 120.218 11.093 4.01 1 C 0679181401 t t t 8.5 1.0 1.0 SDSS
Table 1: Summary of ancillary information used in selecting candidates for XMM observations, and log of the XMM-Newton observations. Column (1): Planck source name. Columns (2) & (3): Right ascension and declination of the Planck source (J2000). Columns (4) & (5): Signal-to-noise ratio of the Planck cluster candidate detection with the MMF3 algorithm in the Planck-maps, and number of methods blindly detecting the candidate. Column (6): quality grade of the SZ detection (A is best). Column (7)–(10): XMM-Newton observation identification number, filter used, on-source exposure time with the EPN camera, and fraction of useful time after cleaning for periods of high background due to soft proton flares (EMOS and EPN camera, respectively). Column (11): category resulting from the pre-selection of the candidates. (12) Confirmed clusters are flagged. indicates double projected systems.

2 Sample selection

2.1 Planck catalogue

In this paper, candidates were chosen from the catalogue derived from the first 15.5 months of data (the “nominal” mission). The processing status, calibration, and map versions were those of March 2011. The detection and quality assessment of the cluster candidates followed the general procedure described in Planck Collaboration VIII (2011). Briefly, a blind cluster search was performed with three methods: the matched multi-frequency filter “MMF3” developed by Melin et al. (2006); an independent matched multi-frequency filter “MMF1”; and the PowellSnakes algorithm (PWS; Carvalho et al. 2009, 2011). Candidates then underwent internal SZ quality checks, removing spurious detections (e.g., association with artefacts or galactic sources), and assessment of the SZ signal detection. The signal assessment included quantitative criteria such as the signal-to-noise ratio and the number of methods blindly detecting the candidate, , as well as a qualitative assessment based on visual inspection of the frequency maps, reconstructed SZ images, and the frequency spectra for each cluster. The latter procedure is summarised in an SZ quality grade, , as described in Planck Collaboration I (2012).

Previously known clusters were identified via cross-correlation with catalogues and NED/Simbad queries. Possible counterparts were searched for within a 5 radius of the Planck position, allowing us to assign two further external reliability flags:

  • association of a FSC (Faint Source Catalogue) or a BSC (Bright Source Catalogue) RASS source (Voges et al. 1999, 2000) or an excess of counts (with corresponding signal-to-noise ratio) in the RASS [0.5–2] keV image.

  • galaxy over-density in the Digitized Sky Survey (DSS) red plates3, from a visual check. In the Sloan Digital Sky Survey (SDSS) area4, two independent galaxy detection algorithms were applied to the DR7 galaxy catalogues (Fromenteau et al., in preparation; Li & White, in preparation). Both algorithms use photometric redshift information. Quality match criteria were assigned based on cluster richness or the over-density signal-to-noise ratio.

2.2 XMM-Newton target selection

The resulting targets are listed together with their SZ quality flags in Table 1. The range of signal-to-noise ratios, , is wide, with nearly uniform coverage, so that the validation results can be useful for defining the final signal-to-noise ratio for the Planck Cluster Catalogue. We considered lower signal-to-noise ratios than the previous validation programme, with 9 targets at and a median S/N of , as compared to 5.1 previously (for 10.5 months of survey data). A priori, this allows us to reach lower flux or higher redshift. To further push the sample towards high redshift, we discarded candidates with estimated size greater than . Although the large positional uncertainty of Planck candidates makes the search for a DSS counterpart non-trivial, the brightest galaxies of clusters at are generally visible in DSS (e.g., Fassbender et al. 2011). We thus also used DSS images to select high- clusters. Half the targets, labelled PHZ (potentially at high z) in Table 1, have no visible counterpart in DSS red plates. These are obviously riskier candidates, particularly those with low or .

As previous validation observations have shown, the association of a SZ candidate with a RASS FSC or BSC source is not in itself sufficient to confirm the candidate, as chance association with a point source is always a possibility. Conversely, a candidate with no counterpart in the RASS catalogue may well be a bona fide cluster. With this campaign, in combination with the previous observations, we also aim to address the use of RASS data as an indicator of candidate reliability. In the sample of 36 candidates observed previously, thirteen candidates were associated with a BSC source and seventeen candidates with an FSC source. Only six SZ candidates had no FSC/BSC counterpart, of which the three confirmed candidates were detected in RASS at a signal to noise ratio of . To better span the range of external RASS flags, we chose ten candidates with no FSC or BSC counterpart, six of which correspond to a RASS . Of the remaining five candidates, only one is associated with a BSC source and four are associated with an FSC source. The RASS association for all XMM-Newton validation targets is summarised in Table. 3.

Figure 1: XMM-Newton [0.3–2] keV energy band images of the three unconfirmed cluster candidates centred on the SZ position (yellow cross). The red circles indicate the presence of an extended source. Green squares in the right panel are positions of galaxies in the SDSS over-density.

Finally one candidate, PLCK G210.620.4, was specifically chosen to further test our SDSS-based confirmation of very poor SZ candidates. PLCK G210.620.4 is the lowest SZ signal-to-noise candidate, detected at by one method only, with a grade and no significant signal in RASS data. However, the galaxy-detection algorithms (Sec. 2.1) that we used indicated that the candidate is associated with an SDSS galaxy over-density at .

3 XMM-Newton observations and data analysis






Figure 2: XMM-Newton [0.3-2] keV energy band images of confirmed cluster candidates. North is up and East is to the left. Image sizes are on a side, where is estimated from the relation of Arnaud et al. (2010) assuming standard evolution. Images are corrected for surface brightness dimming with , divided by the emissivity in the energy band, taking into account galactic absorption and instrument response, and scaled according to the self-similar model. The colour table is the same for all clusters, so that the images would be identical if clusters obeyed strict self-similarity. A yellow cross indicates thePlanck position and a red/green plus sign the position of a RASS-BSC/FSC source. The clusters are sorted according their estimated redshift. For the double systems (last two rows) the middle and right panels show the two components and the left panel the wavelet-filtered overall image.

Candidates were observed between 31 July 2011 and 13 October 2011. The observation identification number and observation setup are given in Table 1. Due to a slew failure in the satellite revolution 2132, the PLCK G208.674.4 observation was incomplete, with an EPN exposure time of  ks. The target was observed initially at the end of its summer visibility window, and could only be reobserved five months later. It was replaced with an additional visible candidate, PLCK G329.522.7.

Calibrated event lists were produced with v11.0 of XMM-SAS. Data that were affected by periods of high background due to soft proton flares were omitted from the analysis (Pratt et al. 2007); clean observing time after flare removal is given in Table 1. The status of each SZ candidate is also given in Table 1: 12 of the 15 candidates are confirmed to be real clusters, among which two are double systems. XMM-Newton images of unconfirmed candidates are shown in Fig. 1; confirmed candidates are shown in Fig. 2.

We derived redshifts and physical parameters of the confirmed candidates as described in Planck Collaboration VIII (2011); Planck Collaboration I (2012). Cleaned XMM-Newton data were pattern-selected. Each photon was then assigned a weight equivalent to the ratio of the effective area at the photon energy and position to the central effective area, computed with SAS task evigweight. Images and spectra were extracted using this weight, assuring full vignetting correction (see Arnaud et al. 2001). Bright point sources were excised from the data and the background was handled as described in Pratt et al. (2010). The particle-induced background (PB) was estimated using a stacked event list built from observations obtained with the filter wheel in closed position. The cosmic X-ray background was modeled using a PB-subtracted spectrum of an annular region external to the cluster emission.

In the spectroscopic analysis, the hydrogen column density was fixed at the 21-cm value of Kalberla et al. (2005). The redshift was estimated by fitting an absorbed redshifted thermal model to the spectrum extracted within a circular region corresponding to the maximum X-ray detection significance. The quality of the estimate was characterised by the quality flag as introduced in Planck Collaboration IX (2011). was set to when the redshift could not be constrained due to the lack of line detection. corresponds to ambiguous estimate, when the spectral fit as a function of exhibited several minima that could not be distinguished at the confidence level. corresponds to a well constrained redshift (i.e., a single minimum).

Surface brightness profiles centred on the X-ray peak were extracted from bins in the [0.3–2] keV band for each instrument independently, background subtracted, co-added and rebinned to per bin. 3D gas density profile were obtained using the regularised non-parametric method of direct deprojection and PSF deconvolution of the surface brightness profile developed by Croston et al. (2006). Global cluster parameters are estimated self-consistently within via iteration about the relation of Arnaud et al. (2010), assuming standard evolution,

The quantity , is defined as the product of , the gas mass within , and , the spectroscopic temperature measured in the [0.15–0.75]  aperture. In addition, , the X-ray luminosity inside , is calculated as described in Pratt et al. (2009). The SZ flux was then re-extracted, being calculated with the X-ray position and size fixed to the refined values derived from the high-quality XMM-Newton observation. The X-ray properties of the clusters and resulting refined values are listed in Table 2.

4 XMM-Newton validation outcome

4.1 False cluster candidates

For the three candidates shown in Fig. 1, no obvious extended X-ray sources were found within of the Planck position. We followed the maximum likelihood procedure described by Planck Collaboration VIII (2011) to find all extended sources in the field detected at the level. We then assessed whether they could be the counterpart of the Planck candidate from their position and X-ray flux, using the relation between the X-ray flux in the [0.1–2.4] keV band, , and the SZ flux established by Planck Collaboration I (2012):

(1)

Here is the correction, neglecting its temperature dependance.

Name RA Dec
[h:m:s] [d:m:s] [kpc] [keV]
PLCK G219.934.4 04:54:45.4 20:17:06.6 0.66 2 2.7 1048
PLCK G348.425.5 19:24:56.1 49:27:02.1 0.25 2 6.4 1020
PLCK G352.124.0 19:20:59.7 45:51:02.2 0.77 1 2.0 925
PLCK G305.944.6 00:23:38.9 72:24:06.1 0.30 2 6.7 1178
PLCK G208.674.4 02:00:16.4 24:54:54.4 0.90 1 3.2 1012
PLCK G130.117.0 01:30:51.3 45:17:54.9 0.20 2 6.3 963
PLCK G239.940.0 04:46:47.2 37:03:49.7 0.74 1 4.4 1033
PLCK G204.715.9 07:34:27.3 14:16:50.2 0.34 2 4.6 1121
PLCK G011.240.4 21:00:37.6 33:08:05.7 0.46 2 2.9 833
PLCK G147.316.6 02:56:25.3 40:17:18.7 0.62 1 2.4 1042
PLCK G329.522.7 A 18:33:00.3 65:33:20.0 0.24 2 4.9 917
PLCK G329.522.7 B 18:33:33.7 65:26:39.2 0.46 2 2.6 872
PLCK G196.745.5 A 03:42:54.2 08:40:58.2 0.57 1 5.6 820
PLCK G196.745.5 B 03:43:02.4 08:46:09.8 0.42 1 2.4 826
5
Table 2: X-ray and SZ properties of the confirmed Planck sources. Columns (2) & (3): Right ascension and declination of the peak of the X-ray emission (J2000). Column (4): redshift from X-ray spectral fitting. Column (5): Quality flag for the X-ray redshift measurement (see Sec. 3). Column (6): Total EPIC count rates in the [0.3-2] keV band, within the maximum radius of detection given in Column (7). Columns (8)–(14): is the radius corresponding to a density contrast of 500, estimated iteratively from the relation, is the product of the gas mass within and the spectroscopic temperature , and is the total mass within . is the luminosity within in the [0.1-2.4] keV band. is the spherically integrated Compton parameter measured with Planck, centred on the X-ray peak, interior to the estimated with the X-ray observations.

Plck g196.468.3 and PLCK G310.527.1

PLCK G196.468.3 was classified as PHZ (potentially at high ). Analysis of the XMM-Newton data on PLCK G196.468.3 revealed two extended sources at 9.9 and 11.8 from the SZ position. The former corresponds to a RASS-FSC source. Both sources are too far away to be the X-ray counterpart of the Planck candidate. A RASS-FSC source is located at 5.2 from the SZ position and likely contributes to the signal derived from RASS data at the Planck source location. However, the comparison of its surface brightness profile with the XMM-Newton PSF shows that it is consistent with a point source. We thus conclude that PLCK G196.468.3 is a false detection.

PLCK G310.527.1 was also classified as PHZ. Two extended X-ray sources were detected at 10.5 and 2.5 from the SZ position, respectively. The former is too far away to be the X-ray counterpart, while the latter is very weak. Analysis of the surface brightness profile confirmed that it is extended. The detection radius is small, and the spectrum extracted from this region is too poor to put robust constraints on the redshift or the temperature. However, using the relation (Eq. 1) and the measured X-ray flux, we can put an upper limit on assuming a redshift as high as and taking into account a factor of two dispersion around the relation. For a temperature of and , we derive a flux within the detection radius of  erg s cm. Assuming that this flux is close to the total, this gives an upper limit on the SZ flux of , nearly an order of magnitude smaller than the Planck value . Moreover, the SZ significance drops under 2 when the flux is re-extracted at the X-ray position. We conclude that this candidate is also a false detection.

Both of these false candidates were detected by two methods, with a medium quality grade of and at and , respectively. A quality grade is thus not sufficient to ensure candidate validity at these signal-to-noise ratios. On the other hand, all candidates down to that have been followed up by XMM-Newton have been confirmed.

Plck g210.620.4

PLCK G210.620.4 is associated with an SDSS cluster. The SDSS search algorithm identified a galaxy over-density of 77 members at a photometric redshift of , consistent with the spectroscopic redshift of the brightest cluster galaxy (BCG) at . The barycentre of the concentration and the BCG are located 1.5 and 5 from the Planck position (see Fig. 1), respectively. The X-ray analysis revealed the presence of an extended source, centred on the BCG, detected at 3.3 in the [0.3–2] keV image. However, the source is very faint and more reminiscent of a group of galaxies than of a rich cluster. This is confirmed by the X-ray spectroscopic analysis. Extracting and fitting the spectrum with an absorbed thermal model at , we measured a temperature within the detection radius of and a flux of . Using Eq. 1 as above, the upper limit on the corresponding SZ flux is , more than 10 times lower than the Planck value of . The X-ray source is too weak to be the Planck counterpart and we conclude that the candidate is not a cluster.

In the previous XMM-Newton validation run, the two candidates potentially associated with SDSS clusters were confirmed, including PLCK G193.346.1 at . This showed that SDSS can robustly confirm candidates up to such high . It is instructive to compare PLCK G210.620.4 with PLCK G193.346.1. In both cases the search algorithm found a rich concentration of galaxies, as expected for Planck clusters. The masses, reconstructed from the luminosity function, are and , respectively, i.e., the false candidate has a larger mass. In both cases, the galaxy distribution appears rather loose (compare Fig. 1 right panel and Planck Collaboration I 2012, Fig. 5). The XMM-Newton observation revealed that PLCK G193.346.1 is a double peaked cluster, i.e., a dynamically perturbed cluster with an ICM distribution consistent with the galaxy morphology. In view of the XMM-Newton image, the galaxy concentration at the location of PLCK G210.620.4 is likely a filamentary structure where only the part around the BCG is virialised and contains gas that is hot enough to emit in X-rays. This would also explain the large offset between the BCG position and the galaxy concentration barycentre, which is much larger than in the case of PLCK G193.346.1. These two cases illustrate the difficulty of distinguishing between massive clusters and pre-virialised structures with rather shallow SDSS data at high . Beyond luminosity and mass estimates, important diagnostics include the offset between the SZ, BCG, and barycentre, as well as the galaxy distribution morphology, if available, and other ancillary data, such as significant RASS emission. These factors must all be considered for firm confirmation of low signal-to-noise-ratio SZ detections. On the other hand, we cannot be sure that the apparent SZ signal is purely due to noise, and cannot exclude a contribution from the pre-virialised structure itself, especially if it corresponds to a warm filament along the line of sight.

4.2 Confirmed candidates

Twelve of the 15 candidates are confirmed as real clusters, of which two are double systems as shown in Fig. 2. Physical parameters are given in Table 2. For the two double systems, the cluster closest to the Planck position is labelled A and the other is labelled B in Table 2.

Single clusters

The redshifts of eight clusters are well constrained by the XMM-Newton spectrum (quality flag of ). Three of these clusters, PLCK G219.934.4, PLCK G011.240.4 and PLCK G348.425.5 were classified as PHZ. The first two are indeed at and , respectively, but PLCK G348.425.5 is at . Knowing the precise cluster location with XMM-Newton, we re-examined the DSS image. A bright galaxy is indeed located exactly at the position of the X-ray peak; however, the field is crowded and there is no obvious galaxy concentration around that BCG. This explains our initial mis-classification.

Figure 3: A composite image of the central of PLCK G147.316.6, based on imaging data from NOT/MOSCA ( and ) and TNG/DOLORES ( and ). Boxes: cluster galaxies spectroscopically confirmed with Gemini (excluding the two galaxies at ). North is up and East to the left. The green contours are isocontours of the wavelet filtered XMM-Newton image. The white contours show the luminosity distribution of the red sequence galaxies indicated by red symbols in Fig. 15, smoothed with a Gaussian filter. The plotted contour levels are at (10,20,30) times the rms variation in the luminosity distribution.

The redshift determination for three single clusters is more uncertain. There are several minima that cannot be distinguished at the confidence level (). As proposed by Planck Collaboration I (2012), we used the X-ray versus SZ properties to eliminate unphysical solutions, as well as DSS data. This is detailed in Appendix A.1. The XMM-Newton analysis gives three possible redshifts for PLCK G147.316.6: 0.4, 0.62, and 1.1, the last being the best-fitting value. The cluster has an interesting double-peaked morphology. It is likely an on-going merger of two nearly equal mass systems (Fig. 3). The analysis of imaging data obtained with the Telescopio Nazionale Galileo La Palma (TNG) telescope and the Nordic Optical Telescope, as well as spectroscopic data obtained at Gemini, are detailed in Appendix A.2. We confirm a redshift of .

The spectral analysis of PLCK G208.674.4 gives a single minimum at , in very good agreement with SZ versus X-ray properties. However we assign a quality flag of since the statistical quality of the spectrum is poor due to the short exposure time. Furthermore the DSS image is ambiguous: although there is no visible galaxy at the X-ray maximum, the centroid of the large scale X-ray emission is close to a bright DSS galaxy.

In summary, of the seven candidates we classified as PHZ, two are false, four are indeed at , and one is at a low redshift of . In addition to those clusters which were classified as PHZ, two further clusters, PLCK G239.940.0 and PLCK G208.674.4, are most likely at high .

Multiple systems

In PLCK G196.745.5, two clusters, separated by arcmin, lie within the Planck position error box: PLCK G196.745.5A at 2.34 and PLCK G196.745.5 B at 3.9 from the SZ position. In view of the Planck resolution, to depending on frequency (Mennella et al. 2011; Planck HFI Core Team 2011), both clusters certainly contribute to the SZ signal. It is likely a chance association, although given the uncertainty in the redshifts, a binary system cannot be ruled out (see Appendix A.1).

Figure 4: Distance of blind SZ position to X-ray position, , as a function of , normalised to the cluster size for single confirmed systems. The clusters are colour-coded according to redshift.
Figure 5: Histogram of the distance between the X-ray peak determined from the XMM-Newton validation observations and the Planck SZ position for all clusters (orange filled) and those associated with a source from the RASS Faint Source Catalogue or Bright Source Catalogue (red hatched). The histogram of the distance between the X-ray peak and the RASS source position is plotted for comparison (blue hatched).

In PLCK G329.522.7, the cluster PLCK G329.522.7 A lies about 1 from the Planck position, while the second object is about 8 away. From the values and redshift estimates, cluster B is expected to have a flux 1.8 times smaller than that of cluster A, thus contributing to the total flux. Its contribution to the blind signal may differ, as the blind signal is extracted using a single component model found roughly peaked at cluster A. Indeed comparison of such a single component extraction with that using a double component model (with flux ratio fixed to the X-ray constraint) suggests a contamination from cluster B of about . In summary, PLCK G329.522.7 A is the main contributor to the SZ detection, although PLCK G329.522.7 B certainly contributes. The redshifts of the two clusters are well determined, and , respectively, showing that they are not physically related. This double system is thus a chance association on the sky.

Overall, we have found four double systems and two triple systems among the 43 Planck candidates confirmed by XMM-Newton, i.e., 14 % multiple systems. Since the XMM-Newton validation follow-up observations are neither representative nor complete, this fraction of multiple systems cannot be extrapolated to the population at large; however, it is more than five times larger than the fraction of cluster pairs separated by less than 10 ( objects) in the whole MCXC X-ray catalogue compilation (Piffaretti et al. 2011). This is clearly a selection effect due to confusion in the large Planck beam, which it might be necessary to take into account for a precise estimate of the selection function.

Figure 6: The new SZ-discovered Planck single objects compared to clusters from the ROSAT All-Sky Survey catalogues in the  plane. Green points represent Planck clusters previously confirmed with XMM-Newton (Planck Collaboration IX 2011; Planck Collaboration I 2012) and red points are the newly confirmed single clusters. The X-ray luminosity is calculated in the [0.1–2.4] keV band. Catalogues shown are REFLEX (Böhringer et al. 2004), NORAS (Böhringer et al. 2000), BCS (Ebeling et al. 1998), eBCS (Ebeling et al. 2000) and MACS (Ebeling et al. 2007). The solid line is the REFLEX flux limit, the dotted line is the HIFLUCGS flux limit of and the dashed line is from the MACS flux limits.

4.3 Planck position reconstruction uncertainty

The Planck position reconstruction uncertainty is driven by the spatial resolution of the instruments. The positions determined by the Planck detection algorithm are compared to the precise XMM-Newton positions in Fig. 4 and Fig. 5, where we put together all validation observations of single systems. The mean offset between the Planck and the XMM-Newton position is , with a median value of , as expected from Planck sky simulations (Planck Collaboration VIII 2011, Fig. 7 left). For 70 and 86 % of the clusters, this offset is less than and , respectively. The assumed positional uncertainty of up to is certainly conservative and an offset of is actually very unlikely. This needs to be taken into account when searching for possible counterparts in ancillary data or follow-up observations.

The offsets of five sources are greater than . Three of those objects are very diffuse, likely dynamically unrelaxed systems, at relatively low , including the prominent outlier PLCK G18.7+23.6 at (Fig. 4, purple point). As noted by Planck Collaboration IX (2011) a real, physical offset between the X-ray and SZ peak may contribute to the overall offset for this type of cluster. In all cases but one, the offset remains smaller than the cluster size . The notable exception is PLCK G11.240.4 (Fig. 4). The XMM-Newton position of this cluster is or from the Planck position. The peak in the SZ reconstructed map is also away from from the Planck position. This cluster is detected by only one method and has a low quality grade , being located in a particularly noisy region of the Planck map. This is likely to complicate the estimate of the cluster position.

Finally, we note that the position reconstruction uncertainty is on average smaller than for the ESZ sample that peaks at (Planck Collaboration VIII 2011, Fig. 7 right). This is likely the result of the higher redshift range considered here. Indeed, at this redshift the sources are more compact and their position is easier to reconstruct. Furthermore, possible physical offsets are expected to become negligible as they become unresolved.

Figure 7: Relation between apparent SZ signal () and the corresponding normalised parameter for single systems confirmed with XMM-Newton (green and red points). Black points show clusters in the Planck-ESZ sample with XMM-Newton archival data as presented in Planck Collaboration XI (2011); The blue lines denote the scaling relations predicted from the REXCESS X-ray observations (Arnaud et al. 2010). The grey area corresponds to median values in bins with standard deviation.

4.4 New clusters in the and plane and Planck sensitivity

The present validation sample covers a wide range of redshift, , and SZ flux, . As expected from the lower signal-to-noise ratio considered and the deeper sky coverage (Sec. 2), the range is lower than that of the previous validation sample, . Although not perfect, the strategy to preferentially select high- clusters was successful, with five clusters found at , including three PHZ candidates. The full XMM-Newton validation sample (single objects only) is shown in the - plane in Fig. 6. We continue to populate the higher part of the - plane and confirm Planck can detect clusters well below the X-ray flux limit of RASS-based catalogues, ten times lower than REFLEX at high , and below the limit of the most sensitive RASS survey (MACS). The figure makes obvious the gain in redshift coverage as compared to the RASS-based catalogues.

Figure 8: The new SZ-discovered Planck single objects (blue, red and green symbols) in the plane. For comparison, black points show known clusters from the ESZ Planck catalogue with archival XMM-Newton data (Planck Collaboration XI 2011). are estimated from and the  relation of Arnaud et al. (2010).
Name S/N Ass. Run Confirmed
RASS
PLCK G271.231.0 18.4 4.82 0.03 B 2 Y
PLCK G286.631.3 10.7 3.23 0.07 B 2 Y
PLCK G018.7+23.6 8.8 6.35 0.09 B 2 Y
PLCK G305.944.6 7.9 2.38 0.05 B 4 Y
PLCK G234.220.5 7.2 2.45 0.02 B 3 Y
PLCK G285.023.7 7.1 3.80 0.04 B 2 Y
PLCK G060.1+15.6 7.0 2.86 0.06 B 3 Y
PLCK G268.528.1 6.6 0.48 0.02 F 3 Y
PLCK G171.940.7 6.1 5.78 0.06 B 2 Y
PLCK G266.627.3 5.6 0.84 0.02 F 3 Y
PLCK G241.228.7 5.1 1.28 0.02 B 2 Y
PLCK G019.1+31.2 5.1 2.93 0.04 B 3 Y
PLCK G277.851.7 4.9 1.63 0.03 F 1 Y
PLCK G208.674.4 4.3 0.45 0.02 F 4 Y
PLCK G250.0+24.1 4.2 0.73 0.04 F 1 Y
PLCK G286.338.4 4.1 1.51 0.04 B 1 Y
PLCK G285.617.2 3.8 1.24 0.02 F 2 Y
PLCK G130.117.0 3.7 1.92 0.03 4 Y
PLCK G200.928.2 3.7 0.77 0.03 F 3 Y
PLCK G235.6+23.3 3.6 0.86 0.02 F 3 Y
PLCK G262.2+34.5 3.5 1.15 0.02 F 3 Y
PLCK G004.519.5 3.3 2.00 0.03 B 1 Y
PLCK G272.9+48.8 3.2 2.60 0.10 F 2 Y
PLCK G205.063.0 3.0 1.44 0.02 F 2 Y
PLCK G348.425.5 2.9 1.72 0.02 F 4 Y
PLCK G292.5+22.0 2.8 2.22 0.04 2 Y
PLCK G100.230.4 2.8 1.27 0.03 F 2 Y
PLCK G226.116.9 2.3 F 1
PLCK G193.346.1 2.2 0.45 0.01 F 3 Y
PLCK G204.7+15.9 2.0 1.32 0.02 4 Y
PLCK G287.0+32.9 1.9 4.01 0.05 2 Y
PLCK G147.316.6 1.8 0.59 0.06 F 4 Y
PLCK G011.240.4 1.8 0.27 0.01 4 Y
PLCK G210.6+17.1 1.7 0.86 0.01 3 Y
PLCK G196.468.3 1.7 4
PLCK G070.821.5 1.6 1
PLCK G262.740.9 1.3 2.26 0.02 F 2 Y
PLCK G113.174.4 1.2 F 3
PLCK G343.443.4 1.2 F 1
PLCK G239.940.0 1.0 0.66 0.01 F 4 Y
PLCK G352.124.0 0.7 0.52 0.01 4 Y
PLCK G219.934.4 0.7 0.53 0.02 4 Y
PLCK G317.454.1 0.4 1
PLCK G310.5+27.1 0.9 4
PLCK G210.6+20.4 1.0 4
6
Table 3: RASS information for single confirmed clusters and false candidates. (1) Name of the candidate. (2) signal-to-noise ratio of the RASS count rate in the [0.5–2] keV band, measured within a region of radius centred on the SZ candidate position. (3) Flux in the [0.1–2.4] keV band as measured with XMM-Newton within . (4) association with a source from the RASS Faint Source Catalogue (F) or Bright Source catalogue (B) published by Voges et al. (1999, 2000). (5) number of the XMM-Newton validation run. (6) Confirmed clusters are flagged.

We confirm our previous results on the  relation. Most clusters are consistent with the REXCESS prediction:

(2)

with . However, all clusters below a normalised lie above the predicted relation and the bin average deviation increases with decreasing (Fig. 7). As noted by Planck Collaboration I (2012), this is an indication of Malmquist bias.

Figure 8 shows the new Planck clusters confirmed with XMM-Newton in the plane (single objects only). The minimum mass increases with redshift, an indication of an increase of the mass detection threshold with . Such an increase is expected from the fact that clusters are not resolved by Planck at high ; however, we clearly confirm that Planck can detect clusters above . Two clear outliers in the plane are evident in Fig. 8. They correspond to the lowest flux clusters PLCK G11.240.4 and PLCK G268.528.1 at and , respectively (Fig. 7), lying in the region most affected by the Malmquist bias. PLCK G11.240.4 is the cluster mentioned in Sec. 4.3, which is detected with a large offset between the Planck position and the X-ray peak, due to its lying in a region with a noisy background. The blind signal is two times higher than the signal extracted at the X-ray position. This is a clear case of a detection boosted by specific local noise conditions.

5 Using RASS data in the construction of the Planck cluster catalogue

5.1 Position refinement

The positions of the associated FSC and BSC source are indicated in the individual XMM-Newton image of each candidate in Fig. 2, and for previous observations, in Fig. 3 and Fig. 2 published in Planck Collaboration IX (2011) and Planck Collaboration I (2012). Comparing the positions of the SZ candidates and their FSC/BSC counterparts with the X-ray peaks determined from the XMM-Newton validation observations, we notice that the FSC/BSC position is a better estimate of the position of the cluster than the position returned by Planck alone. Most of the FSC/BSC sources are located within 1′ of the XMM-Newton position versus 2′ for the Planck-SZ position (see Fig. 5). Thus, the association with a faint or bright RASS source can be used to refine the SZ position estimate.

5.2 X-ray flux estimate

Figure 9 summarises the comparison between RASS and XMM-Newton unabsorbed fluxes computed in the [0.1–2.4] keV band. The XMM-Newton flux is given in Table 3. Fluxes measured in an aperture of 5′ centred on the Planck candidate position from RASS images are referred to as “blind.” Here the RASS count rate is converted to flux assuming a typical redshift of , temperature of , and the 21–cm value. All other fluxes are recomputed in an aperture corresponding to , centred on the X-ray peak as determined from the XMM-Newton validation observations, and using the measured temperature and redshift to convert XMM-Newton or RASS count rates to flux.

These figures indicate that the RASS blind fluxes and the RASS fluxes measured within are in relatively good agreement, with a slight underestimate at high fluxes (left panel). RASS and XMM-Newton fluxes measured within are also in relatively good agreement, although with a slight underestimate together with increased dispersion at low fluxes (middle panel). As a result, RASS blind fluxes slightly underestimate the “true” XMM-Newton flux measured within , by at . The underestimate increases with decreasing S/N (right panel).

In view of this agreement, we conclude that the RASS blind flux can be used to estimate the exposure time required for X-ray follow-up of a Planck candidate, once confirmed at other wavelengths. The main limitation is the statistical precision on the RASS estimate.

Figure 9: Relations between unabsorbed X-ray fluxes measured in the  keV band. Blind fluxes are measured in a 5′ aperture centred on the Planck position; all other fluxes are measured in an aperture corresponding to centred on the XMM-Newton X-ray peak. Left panel: Blind RASS flux vs RASS flux. Middle panel: RASS flux vs XMM-Newton flux. Right panel: Blind RASS flux vs XMM-Newton flux.
Figure 10: Density map of the RASS-Faint Source Catalogue (FSC) with XMM-Newton validation results overplotted. The source density map has been normalised by the median of the pixel density distribution. The source density directly reflects the RASS scanning strategy, with the largest exposure and source density at the Ecliptic poles. Cyan pluses (): confirmed candidates associated with a BSC source. Other confirmed candidates are plotted in green, and false candidates are plotted in red. Pluses (): good association with a FSC source. Crosses (): mis-association with an FSC source. Circles (): no association with a FSC/BSC source. Confirmed candidates with no association are mostly located in low density regions corresponding to the shallower part of the RASS survey.

5.3 Candidate reliability

The association of an SZ candidate with a RASS-B/FSC source is neither a necessary nor a sufficient condition for an SZ candidate to be a bona fide cluster. Putting together the results from all XMM-Newton validation observations for a total of 51 Planck cluster candidates, we find that three of the eight false candidates are associated with an FSC source, while eleven candidates are confirmed without association with a RASS-FSC/BSC source. On the other hand, it is striking that PLCK G266.627.3, the most distant cluster of the sample, with a , is detected at a in RASS, and is in fact found in the RASS Faint Source Catalogue.

RASS source density

It is important to underline that the RASS is not homogeneous, and that neither the BSC nor the FSC are flux-limited or complete in any way. Using the RASS-BSC and FSC, we computed the source density map of each catalogue and the associated probability that a Planck candidate will be associated with a B/FSC source within a radius of 5′. The method is described in Appendix B, and the resulting probabilities are given in Table 4.

Figure 10 and Fig. 16 show the RASS-FSC and BSC source density maps with all XMM-Newton validation observations overplotted. The faint source distribution directly reflects the RASS scanning strategy, as evident in Figure 10. In this context, the probability of chance association is also an indication of how well covered the region is and thus on the depth of the X-ray observation at this position. We found a mean probability of association with an FSC source of  % over the whole sky, where is the area corresponding to a circle of 5′ and is the mean density at the position of the candidate, respectively. The corresponding mean probability of association with a BSC source is  %. However, in the best-covered regions of the RASS the probability can reach 95 % for the FSC and 9 % for the BSC, while in the least-covered regions these probabilities drop to 0.4 % and 0.2 %, respectively.

Figure 11: Relation between RASS blind fluxes and SZ fluxes, , for single systems confirmed with XMM-Newton (all validation observations). The RASS flux is the unabsorbed flux computed in the  keV band and measured in a 5′ aperture centred on the Planck position. The points are colour-coded as a function of redshift. Squares are candidates associated with a FSC source while diamonds are candidates associated with a BSC source

BSC source association

All 12 candidates associated with a BSC source are confirmed. This is not surprising. For the BSC, the probability of chance association is relatively low, varying from less than 1 % to 9 %, depending on the sky region. For one cluster, PLCKG305.944.6, the XMM-Newton validation observation reveals that a point source is located at the position of the BSC source. However, the source is labelled as extended in the BSC, and in fact the X-ray emission likely corresponds to a blend of the point source and extended cluster emission that was not resolved with the large ROSAT PSF. This is supported by a comparison of the XMM-Newton and RASS images.

Thus we conclude that the correspondence of a Planck SZ candidate with a RASS-BSC source is a very good indication of there being a real cluster at this position.

FSC source association

For the FSC catalogue, on the contrary, the conclusion is more uncertain because of the larger probability of chance association. Most (18 of 21, i.e., more than 85 %) of the candidates associated with a faint source are indeed confirmed. For the triple system PLCK G214.6+36.9, the FSC source is classified as extended. Its position as given in the RASS catalogue lies between the three clusters and is close to that of a bright XMM-Newton point source. The FSC source is thus in fact a blend of the cluster(s) and of the point source. In only one case, PLCK G193.346.1, does the FSC not correspond to the cluster emission. The XMM-Newton and RASS data shows that it is a point source located away from the cluster centre.

Taking into account PLCK G193.346.1 and the three false candidates associated with an FSC source, we found four cases of mis-associations out of 51 candidates, i.e., 8 %. This is consistent with the mean probability of chance association of 6 % computed above; however, the association with an FSC source is still an indicator of reliability even in the regions of high probability of chance association. For instance, the two highest-redshift clusters () are correctly associated with a faint source, despite both being in the ecliptic pole region where the probability of false association is high. The scanning strategies of Planck and RASS are very similar in that both surveys are deeper in the same regions. In well-covered regions, the association with the faint source catalogue allows us to probe less massive or higher redshift potential clusters. A possible indicator of false association might be the distance between the FSC source and the SZ position, although no strict criterion can be applied. Seventy-five per cent of the false associations correspond to a distance greater than , compared to 2 out of 16 (13 %) for true associations.

No association

Sixteen candidates are not associated with a B/FSC source. Five of these candidates are false and eleven candidates are true sources with no B/FSC source association. As mentioned above, the association with a B/FSC is not necessary for an SZ candidate to be a bona fide cluster. However, we note that the median probability of FSC chance association, a measure of survey depth as discussed Sec. 5.3.1, is 2.1 % for clusters without association, to be compared to 6.7 % for associated clusters (see also Fig. 10). These true clusters with no B/FSC counterpart are located in the shallower part of the RASS survey, which likely explains why they are not associated.

Figure 12: Scaling properties of Planck clusters, colour-coded as a function of redshift. In all figures, and are estimated from the  relation of Arnaud et al. (2010). Top left panel: The scaled density profiles of the new clusters confirmed with XMM-Newton observations. The radii are scaled to . The density is scaled to the mean density within . The thick lines denote the mean scaled profile for each sub-sample. The black line is the mean profile of the REXCESS sample (Arnaud et al. 2010). Other panels: Scaling relations. Squares show the new clusters confirmed with XMM-Newton observations. Points show clusters in the Planck-ESZ sample with XMM-Newton archival data as presented in Planck Collaboration XI (2011). Relations are plotted between the intrinsic Compton parameter, , and the mass (top right panel), between the X-ray luminosity and (bottom left panel) and between mass and luminosity (bottom right panel). Each quantity is scaled with redshift, as expected from standard self-similar evolution. The lines in the left and middle panel denotes the predicted scaling relations from the REXCESS X-ray observations (Arnaud et al. 2010). The line in the right panel is the Malmquist bias corrected relation from the REXCESS sample (Pratt et al. 2009; Arnaud et al. 2010). The new clusters are on average less luminous at a given , or more massive at a given luminosity, than X-ray selected clusters. There is no evidence of non-standard evolution.
Figure 13: Ratio of the Compton parameter to the normalised parameter. Left panel: Variation as a function of redshift. The dotted line is the REXCESS prediction (Arnaud et al. 2010). The full line is the best fit power law and the grey shaded area indicates the uncertainty. Clusters with normalised (green points) were excluded from the fit, to minimise Malmquist bias. Right panel: Histogram of the ratio without and with low flux clusters.

RASS flux and signal-to-noise limit for candidate validation

Unassociated and associated candidates follow the same general correlation between the RASS blind flux, , and the SZ flux, (Fig. 11). This correlation presents some dispersion, with deviations from the mean as large as a factor of three. This is expected from the large statistical errors, as well as from the intrinsic dispersion and dependence of the ratio (Planck Collaboration I 2012) and the difference between the blind and true X-ray fluxes (Sec. 5.2).

Because of this large dispersion, it is not possible to determine a strict RASS flux (or signal-to-noise ratio) limit below which a candidate should be discarded. However, we note that all new clusters have an X-ray flux greater than (grey area on Fig. 11). This flux is consistent with the threshold , as defined from the region affected by the Malmquist bias (see Fig. 7). This RASS flux limit is more than 10 times lower than the REFLEX flux limit of , but still detectable with RASS7. For the confirmed candidates, the minimum signal-to-noise ratio computed from RASS data is . Below that limit, all the candidates were false. All candidates with RASS  are confirmed, and only one false candidate is found for RASS . The latter is an SZ candidate detected at low Planck .

RASS reliability flag

In view of the above results, we conclude the following regarding the most relevant RASS reliability flags:

  • Positional association of a Planck SZ candidate with a RASS-BSC source is a very strong indication that the candidate is a cluster;

  • Positional association of a Planck SZ candidate with a RASS-FSC source at is a good indication of a real cluster;

  • An SZ candidate with no signal at all in RASS is false at very high confidence. Obviously, candidates with low signal-to-noise ratio in a well-covered region are particularly likely to be false.

6 A preview of cluster evolution

With this new XMM-Newton validation campaign, we have now assembled a sample of 37 new single Planck clusters covering a redshift range . With only snapshot XMM-Newton observations, the global properties and density profile of each object are measured accurately enough to allow a first assessment of evolution with redshift. The structural and scaling properties of the sample are illustrated in Fig. 12. We considered three redshift bins, (10 clusters), (19 clusters) and (8 clusters). We confirm our previous finding regarding the scaling properties of these new Planck selected clusters, and do not find any evidence of departure from standard self-similar evolution.

The average scaled density profile (top left panel of Fig. 12) is similar for each bin and is flatter than that of REXCESS, a representative sample of X-ray selected clusters (Arnaud et al. 2010). Once scaled as expected from standard evolution, the new clusters in each redshift bin follow the same trends in scaling properties (Fig. 12): they are on average less luminous at a given , or more massive at a given luminosity, than X-ray selected clusters. On the other hand, they follow the  relation predicted from REXCESS data (Eq. 2).

To study possible evolution with , we plot in Fig. 13 the ratio as function of , including the 62 clusters of the Planck-ESZ sample with XMM-Newton archival data (Planck Collaboration XI 2011). We exclude clusters at low flux, , to minimise possible Malmquist bias (see Sec. 4.4). The best fitting power law gives a slope , with a normalisation of at . The relation is thus consistent with a constant ratio at the REXCESS value of . A histogram of the ratio shows a peak exactly at the REXCESS position. The distribution is skewed towards high ratios, the skewness decreasing if low flux clusters are excluded. This skewness might be intrinsic to the cluster population. It might also reflect a residual effect of the Malmquist bias, clusters with intrinsic high ratio being preferentially detected in SZ surveys.

7 Conclusions

We have presented results on the final 15 Planck galaxy cluster candidates observed as part of a 500 ks validation programme undertaken in XMM-Newton Director’s Discretionary Time. The sample was derived from blind detections in the full 15.5-month nominal Planck survey, and includes candidates detected at . External flags including RASS and DSS detection were used to push the sampling strategy into the low-flux, high-redshift regime and to better assess the use of RASS data for candidate validation. This last phase of the follow-up programme yielded 14 clusters from 12 Planck candidate detections (two candidates are double systems) with redshifts between 0.2 and 0.9, with six clusters at . Their masses, estimated using the relation, range from to . We found an interesting double peaked cluster, PLCK G147.316.6, that is likely an ongoing major merger of two systems of equal mass. Optical observations with NOT, TNG, and Gemini confirmed a redshift of 0.65.

The full XMM-Newton validation follow-up programme detailed in this paper and in Planck Collaboration IX (2011); Planck Collaboration I (2012) comprises 51 observations of Planck cluster candidates. The efficiency of validation with XMM-Newton stems both from its high sensitivity, allowing easy detection of clusters in the Planck mass and redshift range, and from the tight relation between X-ray and SZ properties, which probe the same medium. The search for extended XMM-Newton emission and a consistency check between the X-ray and SZ flux is then sufficient for unambiguous discrimination between clusters and false candidates. We have confirmed the relation between the X-ray flux and the SZ flux, as a function of redshift, and estimated its typical scatter. This relation is used in the validation procedure. By contrast, optical validation is hampered by the relatively large Planck source position uncertainty and the large scatter between the optical observables (such as galaxy number) and the mass (or SZ signal), both of which increase the chance of false associations.

The programme yielded 51 bona fide newly-discovered clusters, including four double systems and two triple systems. There are eight false candidates. Thirty-two of the 51 individual clusters have high quality redshift measurements from the Fe K line. For other cases, the spectral fitting yields several minima as a function of , that cannot be distinguished at the confidence level. We showed that the relation between the X-ray and SZ properties can be used to further constrain the redshift. The new clusters span the redshift range 0.09 to 0.97 and cover more than one decade in , from to .  of single systems is in the range (. These observations provided a first characterisation of the new objects that Planck is detecting:

  • The newly-detected clusters follow the relation derived from X-ray selected samples. This is consistent with the prediction that both quantities are tightly related to the cluster mass.

  • New SZ selected clusters are X-ray underluminous on average compared to X-ray selected clusters, and more morphologically disturbed. The dispersion around the relation may be larger than previously thought and dynamically perturbed (merging) clusters might be under-represented in X-ray surveys. This has implications for statistical studies of X-ray selected samples, either to constrain cosmological models from cluster number counts or to probe the physics of structure formation from the cluster scaling properties. As discussed in detail by Angulo et al. (2012), precise knowledge of the actual scatter between the mass and the observable used in the detection is critical in both applications.

  • We found no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. In particular, there is no significant evolution of the ratio.

Beyond new cluster confirmation and characterisation, we checked the pertinence of the validation process based on Planck internal quality assessments and cross-correlation with ancilliary data. There are eight false candidates in total, all of which were found at . These failures underline the importance of the number of methods detecting the clusters and were used to refine our internal quality flag definitions. All candidates with are confirmed. Galaxy overdensity in SDSS data can confirm candidates up to , although it remains difficult to distinguish between massive clusters and pre-virialised structures at high . The quality of the SZ detection, ancillary data such as significant RASS emission, and the offsets between SZ, BCG, and other positions, must all be considered for firm confirmation. Using the full sample of 51 observations, we investigated the use of RASS-based catalogues and maps for Planck catalogue construction, finding that:

  • Planck clusters appear almost always to be detectable in RASS maps, although there is not a one-to-one correspondence between a RASS-BSC or FSC source and the presence of a cluster.

  • Association of a cluster candidate with a RASS-BSC source is a very strong indication that it is a real cluster.

  • Whether or not there is a RASS-BSC or FSC source, in the RASS maps is a good indication of a true candidate, while is a good indication of a false candidate.

  • The association with a faint or bright RASS source can be used to refine the SZ position estimate. The RASS blind flux can be used to estimate the exposure time required for X-ray follow-up of a Planck candidate, once confirmed at other wavelengths. The main limitation is the statistical precision on the RASS estimate.

The XMM-Newton validation observations could also be used for the verification of Planck performances, showing that:

  • The mean offset between the Planck position and the cluster position is , as expected from Planck sky simulations, and this offset is less than for  % of the clusters.

  • Planck can detect clusters well below the X-ray flux limit of RASS based catalogues, ten times lower than REFLEX at high , and below the limit of the most sensitive RASS survey (MACS).

  • The Planck sensitivity threshold for the nominal survey is arcmin, with an indication of Malmquist bias in the relation below this threshold. The corresponding mass threshold depends on redshift, but Planck can detect systems with at .

  • Overall, there is a high fraction of double/triple systems in the XMM-Newton validation follow-up sample, illustrating the problems of confusion in the Planck beam.

These results illustrate the potential of the all-sky Planck survey to detect the most massive clusters in the Universe. Their characterisation, and the determination of their detailed physical properties, depends on a vigorous follow-up programme, which we are currently undertaking.

Acknowledgements.
The Planck Collaboration thanks Norbert Schartel for his support of the validation process and for granting discretionary time for the observation of Planck cluster candidates. The present work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and the USA (NASA). This research has made use of the following databases: SIMBAD, operated at the CDS, Strasbourg, France; the NED database, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration; BAX, which is operated by IRAP, under contract with the Centre National d’Etudes Spatiales (CNES); and the SZ repository operated by IAS Data and Operation centre (IDOC) under contract with CNES. Based on photographic data obtained using The UK Schmidt Telescope. We further used observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundación Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias (Science Program ID AOT24/11-A24DDT3), on observations made with the Nordic Optical Telescope, operated on the island of La Palma jointly by Denmark, Finland, Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias (Science Program ID 43-016), observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the Science and Technology Facilities Council (United Kingdom), the National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), Ministério da Ciência e Tecnologia (Brazil) and Ministerio de Ciencia, Tecnologìa e Innovaciòn Productiva (Argentina), Gemini Science Program ID: GN-2011B-Q-41. The authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of Mauna Kea has always had within the indigenous Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain. A description of the Planck Collaboration and a list of its members, indicating which technical or scientific activities they have been involved in, can be found at http://www.rssd.esa.int/Planck_Collaboration. The Planck Collaboration acknowledges the support of: ESA; CNES and CNRS/INSU-IN2P3-INP (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).

Appendix A Redshift estimates of confirmed candidates

a.1 Refinement of the XMM-Newton redshift estimate for cases

The redshift determination from XMM-Newton spectral analysis is uncertain for five clusters. There are several minima that cannot be distinguished at the confidence level (). As proposed by Planck Collaboration I (2012), we estimated the and ratios as a function of and compared them to expected values, to eliminate unphysical solutions.

Three possible redshifts were found for PLCK G352.124.0, 0.12, 0.4, and 0.77. The ratio method enables us to exclude the low redshift solution. The solution yields a ratio twice higher than expected, at the limit of the observed dispersion. Furthermore, we confirmed that there is no evidence of galaxy concentrations in the DSS red image at the precise XMM-Newton cluster location. We thus adopt the highest value, , confirming the cluster to be at high .

The best fitting redshift for PLCK G239.940.0, , yields the ratio closest to expectation and is adopted in the further analysis. The lowest solution is very unlikely, yielding a ratio twice as high as expected. The other possible solution is : there are some very faint objects in the DSS images at the XMM-Newton position, although whether those are galaxies is unclear.

Figure 14: EPIC spectra (data points with errors) of PLCK G147.316.6. Only data points above 2 keV are shown for clarity, but data down to 0.3 keV are used in the spectral fitting. The redshift estimate is ambiguous, with the distribution showing three minima. Left panel: the best-fitting thermal model (solid lines) at with the position of the redshifted Fe K line marked. Right panel: Same for the second best solution at , consistent with the optical redshift.

In the case of PLCK G147.316.6, all three redshift solutions, 0.4, 0.62, and 1.03, yield a ratio within the observed dispersion. The best fitting value, , and the second best solution, , are consistent at the confidence level, with values of 125.9 and 128.7 for 132 degree of freedom, respectively. The two models are shown in Fig. 14. The optical measurement is described below (Sec. A.2).

The redshifts of the two components in PLCK G196.745.5 are uncertain. The and ratio methods cannot be used for such double systems, since the individual SZ components are unresolved by Planck. Of the two solutions, and for PLCK G196.745.5A, the latter can be excluded: a clear concentration galaxies at the XMM-Newton location is visible in the DSS images, which thus cannot be at such high (see Sec. 2.2). For PLCK G196.745.5B we adopted the best fitting value, .

Figure 15: A vs.  colour-magnitude diagram of non-stellar objects in the field of PLCK G147.316.6, observed with NOT/MOSCA. Galaxies plotted as red squares, in the region defined by and , form the red sequence constituted by early-type galaxies in the cluster.

a.2 Optical redshift estimate of PLCK G147.316.6

The optical data for PLCK G147.316.6 were taken using Director’s Discretionary Time with DOLORES (Device Optimized for the LOw RESolution), a low resolution spectrograph and imager permanently installed at the TNG telescope (Telescopio Nazionale Galileo La Palma). The camera is equipped with a pixel CCD covering a field of view of (pixel scale of 0.252 per pixel). Exposure times of 3000 s in the and bands were split into 10 single exposures of 300 s each. Exposure times of 4000 s in the band were split into eight separate exposures. Taking advantage of the dither-offsets between single exposures, no separate sky images were required. The images were bias and flat field corrected using IRAF8. For astrometric calibration we used astrometry.net. The average seeing derived from the final images is 0.84, 0.85, and 0.84 in the , , and -bands, respectively. In the final images, we reach signal-to-noise ratios (over the PSF area) of 11, 23, and 8 for unresolved sources of 24th magnitude. The colour composite image allows us to pre-identify the cluster members.

The cluster was also observed using the 2.56-m Nordic Optical Telescope with the MOSCA camera, a mosaic of pixel CCDs. This camera covers a total field of 7.7 7.7, and was used in binned mode. This gives a pixel scale of 0217 per binned pixel. Total exposure times of 900 s were split into 3 dithered exposures of 300 s in each of the SDSS - and -bands in photometric conditions. The telescope was pointed such that the two peaks of the X-ray emission from the cluster would fall in the centreof the mosaic CCD chip that has the best cosmetic quality (named “CCD7”). After standard basic reduction and image registration, the combined images had FWHM of and in the and bands, respectively. Photometric calibration was based on an ensemble of stars in a field located inside the SDSS footprint, observed at similar airmass immediately following the observations of PLCK G147.3