Free Access
Issue
A&A
Volume 570, October 2014
Article Number L10
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201424387
Published online 21 October 2014

© ESO, 2014

1. Introduction

The Sunyaev-Zel’dovich (SZ) effect, which is the distortion caused by the high energy electrons of the intracluster medium on the cosmic microwave background (CMB) photons, has been used by Planck Collaboration XX (2014) to put constraints on cosmological parameters using Planck data. The cosmological constraint comes from matching the galaxy cluster abundance (per unit observable) to the cluster mass function (per unit mass), via a mass-observable relation. At present, the Planck mass-observable calibration is a two-step process (Planck Collaboration XX 2014; Planck Collaboration XXIX 2014): first, the YX mass proxy (YX = TXMgas, Kravtsov et al. 2006) is calibrated against hydrostatic masses (Arnaud et al. 2010), and, second, this calibration is transferred to the YSZ proxy using the measured YSZ values for 71 clusters (with YX values) observed by Planck (Planck Collaboration XX 2014; Planck Collaboration XXIX 2014). This calibration uses clusters spread over a sizeable redshift range (0 <z< 0.5, see Fig. 1) and assumes a given (self-similar) evolution for the mass proxies. The assumption of a self-similar evolution relies on the self-similar model (Kaiser 1991) and numerical simulations (Kravtsov et al. 2006) and lacks an observational determination, i.e. a calibration done using directly measured masses (i.e. weak lensing) at various redshifts.

The cosmological constraints derived by Planck using galaxy clusters differ from those derived (mostly) from the CMB (Planck Collaboration XX 2014). In particular, the best fit CMB cosmology predicts many more clusters than observed, suggesting that the current Planck-estimated masses may be underestimated.

Rozo et al. (2012) use Chandra data from Vikhlinin et al. (2009) and a subsample of an early Planck sample (Planck Collaboration VIII 2011) to argue that the calibration may have an intercept and, possibly, a slope, lower than derived in Planck Collaboration VIII (2011), which uses the same YSZ values but YX values derived from XMM data. The lower intercept agrees with X-ray expectations (Rozo et al. 2012). Rozo et al. (2012) also find no evidence of any evolution in the intercept of the YSZYX relation, giving indirect support to the self-similar evolution assumed by the Planck team but under a number of restrictive hypothesis. For example, they hold fixed the slope of the relation and they suppose to perfectly know the evolution of YX.

Recently, von der Linden et al. (2014) have compared Planck SZ mass estimates to weak-lensing masses and found that Planck masses need to be underestimated at a ~1.6σ level and, with modest evidence, need it even more so with increasing masses. Similar clues have also been found by Sereno et al. (2014) for the von der Linden et al. (2014) sample and for other samples as well. Analyses of these two works assume a self-similar evolution and therefore attribute the tilted scaling to a mass bias. However, Andreon & Congdon (2014) and Gruen et al. (2014) note that the mass slope and the evolution of the scaling may be collinear (degenerate)1. Therefore, the claimed mass bias can be a manifestation of an evolution in the YSZ-mass scaling.

Because the evolution of YSZ is unlikely to be perfectly known and because there are hints in the literature of possible mass-dependent biases that could instead be manifestations of redshift effects, in this paper we independent assess the Planck mass calibration. Our approach is free of evolutionary assumptions and breaks the mass-redshift degeneracy. More specifically, we fit the YSZ-mass allowing (solving for) evolutionary effects that have been neglected in the previous mass calibration. If ignored, these effects lead to a mass-dependent bias. In particular, we use the very same input data as used by the Planck team, in order to emphasize that differences in the results come from differences in the assumptions.

thumbnail Fig. 1

Mass vs. redshift plot of the Planck calibrating sample. Points are M500YX\hbox{$M^{YX}_{500}$} masses. The solid line marks the adopted threshold mass Mthr of our secondary analysis (Sect. 2.1).

We assume ΩM = 0.3, ΩΛ = 0.7, and H0 = 70 km s-1 Mpc-1. Results of stochastic computations are given in the form x ± y where x and y are the posterior mean and standard deviation. The latter also corresponds to 68% intervals, because we only summarized posteriors close to a Gaussian in that way. All log  are in base 10.

2. Revisiting the YSZ calibration

Our starting point is the catalogue of the YX-derived masses, M500YX\hbox{$M^{\rm YX}_{500}$}, and Malmquist-bias-corrected YSZ values delivered by the Planck team for the 71 calibrating clusters2. Their calibration sample has an increasing limiting mass with redshift (Fig. 1).

thumbnail Fig. 2

Mass-YSZ scaling (left-hand panel) and residuals (observed minus expected) as a function of redshift (right-hand panel). The solid line marks the mean fitted regression line. The shaded region marks the 68% uncertainty (highest posterior density interval) for the regression. In the left-hand panel, measurements are corrected for evolution. In the right-hand panel the dotted line shows the E2 /3 dependency assumed by the Planck team.

Following Planck Collaboration XX (2014) and Planck Collaboration XXIX (2014), we fit the data with the function logYSZ,z=logYSZ,z=0+γlogE(z)+s(log(M500YX/Mref))\begin{equation} \log Y_{{\rm SZ},z}=\log Y_{{\rm SZ},z=0} + \gamma \log E(z) + s \left(\log \left(M^{\rm YX}_{500}/M_{\rm ref}\right)\right) \end{equation}(1)allowing a log-normal scatter in YSZ|M500YX\hbox{$Y_{\rm SZ}|M^{\rm YX}_{500}$} around the mean relation, where Mref = 6 × 1014M. In contrast to Planck Collaboration XX (2014) and Planck Collaboration XXIX (2014), we leave the γ parameter free, rather than freezing it to the self-similar value (2 /3). This is an important difference, because we no longer assume the self-similar evolution, and we are trying to solve for the YSZ-mass scaling and its evolution. The fit is performed using the standard errors-in-variables Bayesian regression (Dellaportas & Stephens 1995; see Andreon & Hurn 2013, for an introduction)3.

The mass-YSZ-redshift fit results are shown in Fig. 2. We found that log YSZ scales with logM500YX\hbox{$\log M^{\rm YX}_{500}$} with slope s = 1.51 ± 0.07, with a negligible intrinsic scatter (logYSZ|M500YX=0.06±0.01\hbox{$\log Y_{\rm SZ}|M^{\rm YX}_{500} =0.06\,\pm\,0.01$} dex), logYSZ,z=(1.51±0.07)(log(M500YX/Mref))4.26±0.02+(2.5±0.4)logE(z)\begin{eqnarray} \log Y_{{\rm SZ},z} &=& (1.51\pm0.07) \left(\log \left(M^{\rm YX}_{500}/M_{\rm ref}\right)\right) -4.26\pm0.02 \nonumber \\ && +\,(2.5\pm0.4) \log E(z) \end{eqnarray}(2)with some covariance between the mass slope s and the evolutionary term γ, as shown in Fig. 3. The evolution of the YSZ-mass scaling is well constrained by the data, γ = 2.5 ± 0.4, mainly because a sizeable range of masses are available at a given redshift (because of the vertical spread at a fixed z in Fig. 1). Indeed, if the relation in Fig. 1 were scatterless, then a complete collinearity (degeneracy) between the YSZ-mass slope and the redshift evolution would be present, meaning the slope evolution could not be measured without assuming a slope for the relation. The scatter breaks the collinearity. This is not the case for the Gruen et al. (2014) sample, formed by just seven clusters, as also noted by the authors.

The evolutionary term γ derived from the data, γ = 2.5 ± 0.4, is much larger than the self-similar value assumed in the Planck Collaboration XX (2014) and Planck Collaboration XXIX (2014) fitting, γ = 2 /3. Indeed, γ = 2/3 is rejected by the data because the posterior probability of γ = 2/3 is extremely low (i.e. more than >500 times smaller than the modal value). Furthermore, γ = 2 /3 is just outside the boundary of the 99.9% interval of γ (Fig. 3), and only 0.01% of the Monte Carlo samples have lower γ. Expressed in another way, the γ value derived from the data differs from the one assumed in the Planck Collaboration XX (2014) and Planck Collaboration XXIX (2014) fit in a way that makes the new and old mass slopes different by 4.8 times the slope error quoted in Planck Collaboration XX (2014) and Planck Collaboration XXIX (2014). The disagreement between the evolution pointed out by the data and the one assumed in previous works (e.g. Rozo et al. 2012; Sereno et al. 2014; von der Linden et al. 2014; Planck Collaboration XX 2014; Planck Collaboration XXIX 2014) is also illustrated in the right-hand panel of Fig. 2: YSZ residuals (i.e. observed minus fitted) have a slope of ~2.5 (solid line), not 2/3 (dotted line).

Equation (3) has an immediate consequence: the mass slope of quantities like E−2/3YSZ or any quantity built from it like Planck masses, depends on how the studied clusters are distributed in the mass-redshift plane. If clusters are all at the same redshift, then the correct slope is recovered (because in such a case E(z) is a constant), but if the mass distribution changes with redshift (e.g. because of selection effects), a biased slope is found because of the mismatch between the true data dependency, E2.5, and the assumed evolution, E2 /3. The sample used in the YSZ-mass analysis of Rozo et al. (2012) has a very narrow redshift distribution: 70% of their sample is at 0.04 <z< 0.10. Their slope is, therefore, unaffected by holding the evolutionary term fixed to E2 /3. Instead, the samples used in the YSZ-mass analyses in von der Linden et al. (2014) and Sereno et al. (2014) are distributed over a sizeable redshift range. For these studies it still needs to be evaluated whether the effect interpreted as a mass bias is instead a consequence of the assumption of a self-similar evolution.

thumbnail Fig. 3

Bounds on the mass slope s and evolutionary parameter γ of the YZS-mass relation (68, 95, 99.7%, from inner to outer contours). The point with error bar marks the slope derived in Planck Collaboration XX (2014), Planck Collaboration XXIX (2014) for the assumed γ = 2 / 3.

thumbnail Fig. 4

Ratio between the recalibrated and original mass from Planck YSZ for the cosmological Planck sample (gray background) and for two redshift subsamples (black closed points and green open points). The blue line has a slope equal to 0.28.

In Fig. 4, we show the ratio between the recalibrated and original masses for the whole Planck cosmological sample and for two redshift subsamples. At a fixed redshift, points are aligned on a line of slope 0.28 (the difference in slope between the two calibrations), introducing a mass-dependent bias. The intercept decreases with increasing redshift or, equivalently, the abscissa (mass) at which the two calibrations equal each other increases with redshift. At a fixed redshift, the most massive clusters tend to have recalibrated masses five to ten percent higher than originally quoted (Fig. 4), in the same direction as suggested by von der Linden et al. (2014) and Sereno et al. (2014), by analyses that assume, however, a self-similar evolution. Instead, the originally derived Planck masses are biased high at the lowest masses entering in the cosmological sample (Fig. 4). Therefore, at face value, our new calibration implies more clusters than the original calibration, which already produces more clusters than observed for the Planck CMB cosmology. In other terms, the recalibration does not seem to solve the tension between the CMB- and cluster-based cosmologies. Nevertheless, we want to emphasize that YSZ proxy values of the cosmological sample are derived by the Planck team when assuming their YSZ-mass scaling, not ours. Therefore our conclusion about the role of evolution in solving the tension between the CMB- and cluster-based cosmologies should be re-evaluated after consistently deriving YSZ proxy values with the updated calibration (i.e. with Eq. (2)). A re-extraction of the YSZ values of the Planck cluster catalogue is beyond the scope of this work.

Finally, we would like to emphasize that the current Planck calibration is indirect because it is based on the mass surrogate YX, assumed to be self-similar evolving. Therefore, the found non-self-similar evolution of the YSZ-mass calibration can be genuine, i.e. related to a non-self-similar evolution of the SZ mass proxy itself, or induced by a possible non-self-similar evolution of the mass surrogate YX. In fact, the evolution of the YX-mass scaling lacks an observational determination, since our current knowledge is currently limited to a consistency check (Israel et al. 2014), at the differences of the richness, that has a directly measured (non-) evolution (Andreon & Congdon 2014).

2.1. Alternative treatment of the mass and selection function

We now show that the precise schema used to correct for the sample selection function and for the Malmquist bias has little practical consequences. To correct for the Malmquist bias, we now take the non-Malmquist bias-corrected data and compute the mass function and its evolution from the Multidark simulation (Prada et al. 2012; the data are available in CosmoSim, Riebe et al. 2013). We assume that the selection function is a step function with a threshold mass, Mthr, linearly increasing with E(z): logMthr=2(E(z)1.1)+14.55\begin{equation} \log M_{\rm thr}=2 \ (E(z)-1.1)+14.55 \end{equation}(3)as shown in Fig. 1. This approach properly deals with the Malmquist bias as a stochastic (variable for clusters of the same observed mass) and uncertain correction instead of a deterministic and perfectly known constant according to the Planck and our baseline analyses4. We found s = 1.50 ± 0.07 (vs. s = 1.51 ± 0.07 in the baseline analysis) and 2.1 ± 0.4 (vs. γ = 2.4 ± 0.5 ). The value γ = 2/3 is rejected at 99.9%, showing that our results are robust to the precise treatment of the mass+selection function.

3. Summary

We refitted the values of the 71 Planck clusters in the Planck calibration sample by allowing evolution to be determined from the data, i.e. letting the power of the E(z) term be free, instead of fixed to 2/3. We used the very same input numbers as the Planck team in order to emphasize that differences in the results come from differences in the assumptions. We also repeated the analysis twice, once adopting the Planck team corrections for Malmquist bias and once directly dealing with the mass and selection function. The two analyses consistently found (Eq. (2)) a shallower YSZ-mass relation with a stronger evolution than assumed by the Planck team, i.e. a mass-dependent bias. The assumed (by the Planck team) self-similar evolution, i.e. γ = 2/3, is rejected at 99.9% confidence by their own fitted data, which instead favour an E2.5 ± 0.4(z) evolution. This suggests caution in

interpreting results of works that assume a self-similar evolution and determine, or revisit, the Planck proxy-mass calibration under such an assumption because what is said to be a mass bias could be a neglected evolutionary term.

Our revised calibration is, at first sight, not useful for reducing the tension between the CMB- and cluster-based cosmologies because at face value it increases, instead of reducing, the already too large number of expected clusters for the Planck CMB cosmology. However, this conclusion should be re-evaluated after rederiving YSZ proxy values of the cosmological sample with a proxy-mass relation with an evolution consistent with the data (i.e. with Eq. (2)).

The bottom line is that the non-self-similar evolution of the SZ proxy cannot be ignored: γ = 2.5 (supported by the data) has to be preferred to γ = 2/3 (self-similar). Waiting for a re-derivation of the YSZ values, we suggest estimating the masses for Planck clusters using our Eq. (2) and the YSZ values listed in the Planck catalogue.


1

Sereno et al. (2014) also mention this possibility, but their analysis does not account for it.

2

The table is available at http://szcluster-db.ias.u-psud.fr

3

The fitting code is available at http://www.brera.mi.astro.it/~andreon/fitSZplanck.bug

4

The code used for the fitting is, after minor editing, the same Bayesian code as described in and distributed by Andreon & Bergé (2012), and used in Andreon & Condon (2014) for analysing the determination of the evolution of the richness-mass scaling. We assume weak priors on slope, intercept, and intrinsic scatter and solve for all variables at once. We account, as mentioned, for the scatter, evolution, mass, and selection function. We used a computationally inexpensive 10 million long Markov chain Monte Carlo, discarding the initial 10 thousand elements used for burn-in. By running multiple chains, we checked that convergence is already achieved with short chains.

Acknowledgments

I thank the referee for suggesting the swap of the baseline and secondary analysis to improve clarity and Daniel Gruen for pointing out his paper to me. I’m most grateful to Barbara Sartoris for comments on an early version of this draft and to the Planck team for making the data table used for the YSZ-mass calibration publicly available. I acknowledge Mauro Sereno for enlightening conversations on the subject of this paper.

References

  1. Planck Collaboration VIII. 2011, A&A, 536, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Planck Collaboration XX. 2014, A&A, 571, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Planck Collaboration XXIX. 2014, A&A, 571, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Andersson, K., Benson, B. A., Ade, P. A. R., et al. 2011, ApJ, 738, 48 [NASA ADS] [CrossRef] [Google Scholar]
  5. Andreon, S., & Bergé, J. 2012, A&A, 547, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Andreon, S., & Congdon, P. 2014, A&A, 568, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Andreon, S., & Hurn, M. A. 2013, Stat. Anal. Data Min., 6, 15 [Google Scholar]
  8. Andreon, S., Trinchieri, G., & Pizzolato, F. 2011, MNRAS, 412, 2391 [NASA ADS] [CrossRef] [Google Scholar]
  9. Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Dellaportas, P., & Stephens, D. 1995, Biometrics, 51, 1085 [CrossRef] [Google Scholar]
  11. Gruen, D., Seitz, S., Brimioulle, F., et al. 2014, MNRAS, 442, 1507 [NASA ADS] [CrossRef] [Google Scholar]
  12. Israel, H., Reiprich, T. H., Erben, T., et al. 2014, A&A, 564, A129 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Kaiser, N. 1991, ApJ, 383, 104 [NASA ADS] [CrossRef] [Google Scholar]
  14. Kravtsov, A. V., Vikhlinin, A., & Nagai, D. 2006, ApJ, 650, 128 [NASA ADS] [CrossRef] [Google Scholar]
  15. Prada, F., Klypin, A. A., Cuesta, A. J., Betancort-Rijo, J. E., & Primack, J. 2012, MNRAS, 423, 3018 [NASA ADS] [CrossRef] [Google Scholar]
  16. Riebe, K., Partl, A. M., Enke, H., et al. 2013, Astron. Nachr., 334, 691 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Rozo, E., Vikhlinin, A., & More, S. 2012, ApJ, 760, 67 [NASA ADS] [CrossRef] [Google Scholar]
  18. Sereno, M., Ettori, S., & Moscardini, L. 2014 [arXiv:1407.7869] [Google Scholar]
  19. von der Linden, A., Mantz, A., Allen, S. W., et al. 2014, MNRAS, 443, 1973 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1

Mass vs. redshift plot of the Planck calibrating sample. Points are M500YX\hbox{$M^{YX}_{500}$} masses. The solid line marks the adopted threshold mass Mthr of our secondary analysis (Sect. 2.1).

In the text
thumbnail Fig. 2

Mass-YSZ scaling (left-hand panel) and residuals (observed minus expected) as a function of redshift (right-hand panel). The solid line marks the mean fitted regression line. The shaded region marks the 68% uncertainty (highest posterior density interval) for the regression. In the left-hand panel, measurements are corrected for evolution. In the right-hand panel the dotted line shows the E2 /3 dependency assumed by the Planck team.

In the text
thumbnail Fig. 3

Bounds on the mass slope s and evolutionary parameter γ of the YZS-mass relation (68, 95, 99.7%, from inner to outer contours). The point with error bar marks the slope derived in Planck Collaboration XX (2014), Planck Collaboration XXIX (2014) for the assumed γ = 2 / 3.

In the text
thumbnail Fig. 4

Ratio between the recalibrated and original mass from Planck YSZ for the cosmological Planck sample (gray background) and for two redshift subsamples (black closed points and green open points). The blue line has a slope equal to 0.28.

In the text

Current usage metrics show cumulative count of Article Views (full-text 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 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.