Free Access
Issue
A&A
Volume 631, November 2019
Article Number A96
Number of page(s) 6
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201936423
Published online 29 October 2019

© ESO 2019

1. Introduction

The Lambda cold dark matter (ΛCDM) model is remarkably successful at describing the majority of observations relevant to cosmology. These include probes of the background evolution of our Universe, its early times, and the dynamics and evolution of matter perturbations. Perhaps the most striking success of the ΛCDM model is its ability to explain the observed fluctuations of the cosmic microwave background (CMB), as measured most notably by the Planck satellite, which led to a determination of the cosmological parameters in the ΛCDM model with unprecedented accuracy.

However, some tension is also noted, particularly between the ΛCDM cosmological parameters favoured by Sunyaev–Zel’dovich (SZ) cluster counts, and those derived from the angular power spectra of CMB temperature and polarisation fluctuations (Planck Collaboration XX 2014; Planck Collaboration XXIV 2016). This particular tension appears only when the standard mass calibration (1 − b) = 0.8 of SZ clusters (defined as the ratio of the hydrostatic to true mass) is used, leading to a significant discrepancy in the σ8 (the present-day amplitude of matter fluctuations) versus Ωm (the dark matter density) plane, or equivalently a discrepancy on the calibration parameter when left free, compared to its standard value (see Fig. 3 in Sakr et al. 2018, and Fig. 1 in the present article).

thumbnail Fig. 1.

Confidence contours (68 and 95%) in the Ωm − σ8 plane for various datasets: Planck 2015 and 2018 CMB data, and Planck 2015 SZ cluster counts (using the D16 mass function). Contours are shown for the minimal neutrino case and three neutrinos with free masses (using the CDM prescription for cluster counts). For SZ clusters, the calibration parameter (1 − b) is fixed to its fiducial value of 0.8, with priors added to Ωb from BBN constraints, and to H0 from BAO data. Two additional dashed and dash-dotted contours are shown to illustrate the effects of changing the mass function used (from Despali et al. 2016 to Tinker et al. 2008) and removing the CDM prescription, respectively.

A critical aspect of this tension is that the mass of a cluster (including its dark matter halo) is not directly observable. The above tension is therefore present only when additional priors on cluster masses are used, and it more generally relies on the use of scaling relations between halo mass, redshift, and the observable cluster properties (known as mass proxies). Based on the assumption of hydrostatic equilibrium, the Planck Collaboration used X-ray observations from XMM-Newton in their analysis to derive masses of the intra-cluster medium (Arnaud et al. 2010) and corrected the value by 20% to account for the incomplete thermalisation of the gas. This correction is expressed as a “bias” b in mass, defined as (1 − b) = 0.8 in the fiducial case.

Should the aforementioned tension be confirmed, we would be forced to consider extensions or alternatives to the standard ΛCDM model of cosmology. Since the CMB itself is mostly sensitive to the physics of the early Universe, one can reconcile it with cluster observations by introducing a modification to the cosmological model that has a significant impact preferentially at late times. More specifically, a new theory leading to a lower growth rate of structures would be required in order to predict a lower abundance of clusters, which would better reflect the data.

Such growth can be achieved in several ways; one possibility is to add mass to neutrinos in the standard cosmology (beyond their minimal experimental mass, 0.06 eV). Among other effects, massive neutrinos slow down the growth of matter perturbations during the matter and dark energy-dominated eras on scales smaller than their free-streaming length (see Lesgourgues & Pastor 2012, for a review on the effect of neutrinos in cosmology). Alternatively, one could consider replacing the standard cosmological constant with a different form of dark energy. A fluid with a varying equation of state (the so-called quintessence models) or the addition of a scalar field can achieve the desired result, but a more radical approach consists of modifying the current theory of general relativity itself (see Clifton et al. 2012; Ishak 2019, for extensive reviews).

In Sakr et al. (2018), we pointed out a similar discrepancy between the abundance of local X-ray clusters and CMB cosmology. We found that massive neutrinos do not solve the issue, while a simple modified gravity model parametrised by a modified growth rate, the γ model, can reconcile the Planck cluster mass calibration with the Planck CMB cosmology at the expense of a relatively high γ ∼ 0.9 (a value potentially already excluded by current data; see L’Huillier et al. 2018; Zarrouk et al. 2018). In our study, we combined constraints from the redshift-mass distribution of SZ clusters with the latest CMB measurements from the Planck satellite, performing a Bayesian analysis through a Markov chain Monte Carlo (MCMC) analysis. We varied the parameters of the standard model, the cluster mass calibration (1 − b), and the additional parameters introduced in the context of the extended models we considered. They included the total mass of three massive degenerate neutrinos and the phenomenological γ parameter of our simple modified gravity model.

In Sect. 2, we describe the formalism used for predicting cluster abundances, as well as the extensions to the standard model we examined. In the same section, we also detail the datasets used in this work, and the implementation of the MCMC analysis, which samples the posterior probability distribution function. We present and discuss our results in Sect. 3, while we summarise and discuss our conclusions in Sect. 4.

2. Methodology and data

2.1. Cluster abundances as a cosmological probe

Clusters of galaxies provide tight constraints on cosmological models, as they (mostly) probe the growth of matter fluctuations directly throughout the history of the Universe. Their use as cosmological probes has been extensively studied and exploited in the past. The methodology that we adopted here is closely related to the approach followed by Ilić et al. (2015) and Sakr et al. (2018). We refer the interested reader to those two articles for more details, as we only provide the main steps here.

We computed the expected numbers of clusters (as a function of mass and redshift) in a given cosmological model via integrals over the mass function, with the latter being taken from Despali et al. (2016; thereafter D16). This more recent mass function was chosen over available alternatives (notably the popular Tinker et al. 2008, mass function), although this choice has only a minimal impact on our results and conclusions, as can be seen in Fig. 1 when comparing the filled (D16) and dash-dotted (Tinker) purple contours. The D16 mass function was optimised to work with the virial radius as a definition for clusters, which is the convention we used for X-ray clusters, although this choice had little impact on observables when modelling was done self-consistently (see Sakr et al. 2018). We did, however, use the “critical M500” convention (i.e. the density contrast of clusters is defined as 500 times the critical density of the Universe) for SZ clusters, as the available Planck SZ dataset adopts this particular convention. We used the analytical formulae provided by the D16 authors in order to convert their mass function between their fiducial (virial) cluster mass definition and any other definition (critical M500 here).

In the case of massive neutrinos, we used the so-called CDM prescription to amend the mass function (Costanzi et al. 2013; Castorina et al. 2014) unless otherwise stated. The filled and dashed orange contours of Fig. 1 illustrate the effect of the CDM prescription in the analysis; as discussed in Sakr et al. (2018), we include its effects for the sake of completeness, even though it led to only marginal differences. As seen in Fig. 1, it led to slightly lower values of σ8, and concurrently high values of Ωm, but with slightly increased error bars on both.

When considering a simple modified gravity model, as originally proposed by Linder (2005), we modelled the linear growth rate of structures like so

f ( a ) = Ω m ( a ) γ , $$ \begin{aligned} f(a)=\Omega _m(a)^\gamma , \end{aligned} $$(1)

where γ, the growth index, is here a free parameter generalising the approximation γ ∼ 0.545 in standard gravity (Peebles 1980). The growth rate f is defined as f(a)≡dlnG(a)/dlna, where G is the (scale-independent) linear growth factor of matter density perturbations δm defined as δm(a)=δm, 0G(a)/G(0).

In all extensions to the standard model considered here, we assume the non-linear regime would follow the same threshold δc, and the same virial density, Δv, as in the standard picture. Scaling laws are then used to relate theoretical cluster quantities (i.e. mass) to observable quantities. For the X-ray temperature T versus cluster mass M, we use the standard scaling relation

T = A T M ( h M Δ ) 2 / 3 ( Ω m Δ ( z ) 178 ) 1 / 3 ( 1 + z ) , $$ \begin{aligned} T=A_{T-M}(h\, M_\Delta )^{2/3}\left(\frac{\Omega _{m} \Delta (z)}{178}\right)^{1/3}(1+z) ,\end{aligned} $$(2)

where AT − M is the normalisation parameter, Δ is the density contrast chosen for the definition of a cluster, expressed with respect to the total background matter density1 of the Universe at redshift z, and MΔ is the mass of the cluster according to the same definition (as mentioned before, in our case for X-ray clusters, we used the virial radius and mass definition). Similarly, we used the following scaling relation between the expected mean SZ signal (denoted by the letter Y ¯ $ \bar{Y} $) and mass M500

E β ( z ) [ D A 2 ( z ) Y ¯ 500 10 4 Mpc 2 ] = Y [ h 0.7 ] 2 + α [ ( 1 b ) M 500 6 × 10 14 M ] α , $$ \begin{aligned} E^{-\beta }(z)\left[\frac{D_A^2(z) \bar{Y}_{500}}{10^{-4}\,\mathrm{Mpc} ^2}\right] = Y_*\left[ \frac{h}{0.7} \right]^{-2+\alpha } \left[\frac{(1{-}b)\, M_{500}}{6\times 10^{14}\,M_\odot }\right]^{\alpha }, \end{aligned} $$(3)

where the quantity DA refers to the angular diameter distance and E(z)≡H(z)/H0. We note here the aforementioned hydrostatic mass bias (1 − b), while α, β, and Y* are additional free parameters in the SZ scaling law. A fourth free parameter exists, namely σlnY; it corresponds to the scatter of the log-normal distribution that the measured Y500 is assumed to follow (and of which the mean Y ¯ 500 $ \bar{Y}_{500} $ is given by Eq. (3))2.

2.2. Datasets and numerical tools

We followed a standard MCMC analysis using two publicly available codes: the CosmoMC package (Lewis & Bridle 2002; Lewis 2013) and the Monte Python code (Audren et al. 2013) that respectively use the CAMB and CLASS Boltzmann codes for the computation of relevant cosmological quantities: the matter power spectrum and the power spectra of CMB anisotropies. Those two codes were used to explore our full-parameter space under the constraint of our data sets extracted from CMB measurements (Planck Collaboration XI 2016), baryon acoustic oscillations (BAO) data (Anderson et al. 2014; Font-Ribera et al. 2014), X-ray-selected clusters (Ilić et al. 2015) and an SZ clusters sample (Planck Collaboration XXIV 2016), separated or combined.

The 2015 Planck CMB (Planck Collaboration I 2016) and BAO (Anderson et al. 2014; Font-Ribera et al. 2014) likelihoods were already interfaced with the two aforementioned MCMC codes. The SZ likelihood (Planck Collaboration XX 2014) is present in the CosmoMC package, while we implemented it ourselves in the Monte Python code. We also used our dedicated module for computing the likelihood associated with the X-ray cluster sample (Ilić et al. 2015). While massive neutrinos are already implemented in both the CLASS and CAMB codes, we implemented there the other ΛCDM extension we considered, namely, the γ model. Additionally, we modified our X-ray and SZ clusters likelihood modules accordingly to support the massive neutrinos case and γ parametrisation, all according to the recipes described earlier.

2D confidence regions and 1D constraints on parameters were then derived (as is usual in a Bayesian formalism) from the posterior distribution resulting from the MCMC chains. All parameters limits mentioned in the text are to be read as 68% confidence limits.

3. Results

The long-standing, so-called discrepancy between SZ cluster counts and the Planck CMB measurements in the ΛCDM model is illustrated in Fig. 1 in the Ωm − σ8 plane, with contours from the two probes: in blue for the 2015 Planck CMB data, and purple for SZ cluster counts (the calibration parameter (1 − b) being fixed to the Planck Collaboration fiducial value of 0.8). One can see that the separation between the two contours is significantly large and remains essentially identical when considering the new Planck 2018 data release (green contours). Allowing massive neutrinos with varying mass has the effect of pushing both the CMB contours (light blue and light green contours respectively for the 2015 and 2018 Planck CMB data) and the SZ cluster contours (orange) to higher values of Ωm and lower values of σ8, elongating the two in the direction of the Ωm − σ8 degeneracy, with the CDM prescription shifting the contours slightly to higher Ωm and lower σ8.

The SZ clusters-CMB tension is therefore not relieved by the possible presence of massive neutrinos, leading to essentially identical conclusions in the case where only a sample of local X-ray clusters are used (Sakr et al. 2018). It should be noted that in Fig. 1 we added constraints from Big Bang nucleosynthesis (BBN) and BAO data to restrict Ωb and H0 to realistic values for cluster contours. We checked that considering wider priors results in more elongated contours in the direction of the Ωm − σ8 degeneracy, but without relieving the tension.

3.1. The ΛCDM case

In order to examine this tension in a quantitative way, we combine CMB data and the SZ cluster sample, adding or not adding the constraints obtained with the X-ray cluster sample, and freeing the SZ and X-ray calibration factors, as well as the slope of the mass-SZ relation and its dispersion. In the standard ΛCDM model (as defined by the Planck Collaboration, which includes one massive neutrino with a minimal mass of 0.06 eV) one can see from Fig. 2 that, in the cases of CMB data with SZ and/or X-ray clusters, the constraints on all cosmological parameters are essentially identical. We obtain (1 − b) = 0.58 ± 0.03 and note a clear correlation between the (1 − b) calibration factor for SZ clusters and the AT − M calibration for X-ray clusters (at virial mass), indicating the consistency of the constraints obtained from the two types of cluster samples.

thumbnail Fig. 2.

Confidence contours (68 and 95%) for cosmological parameters Ωm and σ8, as well as mass calibration parameters (1 − b) and AT − M when combining CMB data with X-ray and/or SZ clusters data for D16 mass function.

3.2. Are massive neutrinos a solution?

The possible role of massive neutrinos in alleviating the cluster-CMB tension is discussed in past literature, with somewhat different conclusions, (Rozo et al. 2013; Dvorkin et al. 2014; Mantz et al. 2015; Roncarelli et al. 2015; Bull et al. 2016; Salvati et al. 2018) often without clear indication as to whether conclusions depend on a calibration choice. In Fig. 3, we present the constraints on the SZ cluster calibration parameter (1 − b) from the combination of CMB data sets with the SZ cluster sample in the presence of massive neutrinos (blue contours). We obtain a constraint on the calibration ( 1 b ) = 0 . 585 0.037 + 0.031 $ (1{-}b) = 0.585^{+0.031}_{-0.037} $ essentially identical to the minimal neutrino case. Furthermore, no strong correlation appears in the calibration-neutrino masses plane, similarly to what is obtained from the X-ray sample alone (Sakr et al. 2018). We also notice that the constraint on neutrino masses is quite noticeably improved by combining clusters with the CMB (with 68% limits being reduced by about ∼28%) compared to the CMB-only constraint (black curve). The standard Planck mass calibration ((1 − b)∼0.8) is formally rejected at more than 5σ. However, for illustration purposes, we show the consequence of putting a strong (Gaussian) prior on (1 − b) consistent with the standard Planck calibration, namely (1 − b) = 0.8 ± 0.01 (green contours). In this case, the favoured neutrino masses are found to be non-zero: ∑mν = 0.70 ± 0.15 eV, but the corresponding best-fit model has a Δχ2 ≈ 20 compared to the case without a prior on (1 − b). This result is consistent with previous investigations obtaining non-zero neutrino masses using such a prior. It is therefore important to realise that Planck CMB and Planck SZ cluster counts are perfectly consistent with one another, but as mentioned before, this combination of data requires an SZ mass calibration of (1 − b)∼0.6, significantly lower than the standard calibration (1 − b)∼0.8.

thumbnail Fig. 3.

Confidence contours (68 and 95%) and posterior distributions for neutrino masses ∑mν and mass calibration parameter (1 − b) when combining CMB and SZ cluster data, for D16 mass function (blue). The effect of adding X-ray clusters data is shown in orange, and the effect of a strong prior on (1 − b) centred around the Planck calibration value of 0.8 in green. For comparison, the 1D posterior on neutrino masses from Planck CMB data alone is also shown in black.

3.3. Introducing modified gravity

If the standard calibration (1 − b)∼0.8 is consolidated by observations, the conclusion of the previous sections leads us to examine alternative theories to ΛCDM. One option is to consider a variable dark energy equation of state w instead of a simple cosmological constant. However, Salvati et al. (2018) found that a constant w does not lead to a significant change in the required calibration. Any proposed alternative explanation should lead somehow to a lower amplitude of matter fluctuations at low redshift without altering the remarkable agreement with the CMB seen in the standard cosmological picture. Such a lower σ8 would naturally come from a late-time deviation of the growth of matter fluctuations with respect to ΛCDM, a possibility that could naturally arise from a modification of the theory of general relativity for gravity. In the following, we will consider a phenomenological approach through the introduction of the growth index γ following Eq. (1). We follow the same approach as in previous sections and perform an MCMC analysis to explore the space of the standard cosmological parameters, adding the calibration (1 − b) and the index γ as additional free parameters under the constraints of our data sets. As previously found in Sakr et al. (2018), we expect increasing γ to yield a lower σ8, and thereby a higher calibration (1 − b) – corresponding to clusters being less massive. Similarly to the previous section, we use the combination of CMB and SZ cluster data, with or without adding our X-ray cluster sample to constrain the parameters of the model.

Our results are summarised in Fig. 4; the SZ calibration parameter (1 − b) shows a clear correlation with γ, as anticipated. However, a significant difference appears compared to the case where only the X-ray sample is used: as the SZ sample spans a wide range of (comparatively) higher redshifts, it provides additional constraints on the growth rate. As a consequence, more extreme values of γ are rejected by the SZ cluster data, resulting in closed contours in the γ versus (1 − b) plane (blue contours). The situation remains essentially unchanged when we add the X-ray sample data (orange contours). The constraint on the calibration reads ( 1 b ) = 0 . 602 0.065 + 0.053 $ (1{-}b)= 0.602^{+0.053}_{-0.065} $, and the standard calibration (1 − b)=0.8 is rejected at more than 99%, while γ = 0.60 ± 0.13 is consistent with the growth rate in the standard ΛCDM model (∼0.545). This illustrates the fact that even a simple modified version of the growth rate cannot easily accommodate the calibration (1 − b) = 0.8. Again for illustration purposes, we show the consequences of putting a strong prior on (1 − b), consistent with standard Planck calibration (=0.8 ± 0.01, green contours). The value of γ is then found to be constrained to γ = 0.917 ± 0.088, consistent with conclusions derived from the X-ray sample alone (Sakr et al. 2018).

thumbnail Fig. 4.

Confidence contours (68 and 95%) and posterior distributions for growth index γ and mass calibration parameter (1 − b) when combining CMB and SZ cluster data for D16 mass function (blue). The effect of adding X-ray clusters data is shown in orange, and the effect of a strong prior on (1 − b) centred around the Planck calibration value of 0.8 in green.

3.4. Allowing massive neutrinos and modified gravity

Here, we repeat the previous analysis, allowing both massive neutrinos and a modified gravity through the growth index γ. As can be seen from Figs. 5 and 6, when combining CMB and SZ clusters, a degeneracy is present between the effects of massive neutrinos and modified gravity. The preferred model remains, however, close to the standard picture: ( 1 b ) = 0 . 651 0.12 + 0.043 $ (1{-}b) = 0.651^{+0.043}_{-0.12} $, γ = 0 . 65 0.16 + 0.13 $ \gamma = 0.65^{+0.13}_{-0.16} $ and ∑mν ≤ 0.197, although an SZ calibration of (1 − b)∼0.8 is allowed within the two σ posterior contours. The addition of X-ray clusters reduces the contours, although in a marginal way. When constraining (1 − b)∼0.8 with a strong prior, the preferred parameters are still quite far away from the standard model: γ = 0 . 834 0.094 + 0.080 $ \gamma = 0.834^{+0.080}_{-0.094} $ and m ν = 0 . 34 0.15 + 0.13 $ \sum\nolimits m_\nu = 0.34^{+0.13}_{-0.15} $ eV, but we note that a model with γ ∼ 0.65 and ∑mν ∼ 0.1 is within the two σ contours.

thumbnail Fig. 5.

Confidence contours (68 and 95%) and posterior distributions for growth index γ and mass calibration parameter (1 − b) when combining CMB and SZ clusters data for D16 mass function (blue). The effect of adding X-ray clusters data is shown in orange, and the effect of a strong prior on (1 − b) centred around the Planck calibration value of 0.8 in green.

thumbnail Fig. 6.

Confidence contours (68 and 95%) and posterior distributions for neutrino masses ∑mν and mass calibration parameter (1 − b) when combining CMB and SZ clusters data. The legend and colours are identical to Fig. 5. For comparison, the 1D posterior on neutrino masses from Planck CMB data alone is also shown in black.

3.5. Freeing the amplitude of the late-time matter power spectrum

As a final illustration of what our current data sets actually tell us, we performed the same analysis once again using a combination of CMB and cluster data in the ΛCDM picture, but allowed the present-day amplitude of matter fluctuations σ8 to be an additional free parameter, no longer derived from an extrapolation at late times of the primordial amplitude As. More precisely, we assume that the growth rate between redshift zero and ∼1 (the redshift range containing both the SZ and X-ray cluster samples) follows the standard ΛCDM growth rate, but modified at earlier times through a process about which we do not assume anything. In practice, we rescale the ΛCDM matter power spectrum for z ∈ [0, 1] by a constant factor, so that we obtain the desired input σ8 while keeping the same value of the primordial As. Such an approach allows for (comparatively) more freedom than the γ parametrisation, as the latter only produces a very specific redshift evolution for the growth rate.

In Fig. 7, we show confidence contours in calibration parameters plane, for example, (1 − b) versus AT − M, constrained by the combination of CMB data and our two cluster samples: X-ray and SZ. The various contours shown refer to the cases we examined in the previous sections, with the addition of the new model described in this section (ΛCDM + free σ8 at present epoch). We observe that all contours lay along the correspondence line3 (dashed black) between (1 − b) and AT − M, and that in all cases the preferred SZ mass calibration (1 − b) remains around 0.6 (see 1D posteriors in Fig. 8) even in the case where σ8 is treated as a free parameter. In the latter case, a value of (1 − b)∼0.8 is still acceptable (meaning within reasonable confidence limits), but would thus require a significant modification of the growth rate at earlier times in order to reach the required current value of σ8.

thumbnail Fig. 7.

Confidence contours (68 and 95%) in (1 − b) vs. AT − M plane, for combination of CMB, SZ and X-ray cluster data in various cosmological models, with or without massive neutrinos, with or without free growth index γ, and free σ8. The black dashed line corresponds to the (1 − b)−AT − M relation leading to the same derived (SZ and X-ray) masses for a selected sample of clusters (see text for details).

4. Discussion and conclusions

In the present work, we have examined the constraints on the calibration constants in the scaling relations of the Planck SZ and a local X-ray cluster sample in the context of various cosmological models and assumptions. In our approach, we treat these calibrations (namely (1 − b) for SZ clusters, and AT − M for X-ray clusters) as free parameters without any priors. Neither massive neutrinos nor a modified gravity model (parametrised by the γ growth index) allow us to reconcile the Planck data (CMB and SZ clusters) with the standard SZ mass calibration value of (1 − b)∼0.8. Consequently, this implies that there is no simple solution to the so-called clusters-CMB tension, which may be more accurately described as a tension between Planck data and the empirical calibration of the mass-SZ observable (which yields the (1 − b)∼0.8 value). Despite this observation, we explore the possibility of more extended models allowing massive neutrinos, and a modification of gravity, as well as the approach of freeing the amplitude of matter fluctuations at low redshift (as measured by σ8). Nevertheless, the SZ cluster counts combined with CMB data still prefer a low calibration ((1 − b)∼0.6) in all considered models. We therefore find it striking that no cosmological model appears to prefer the (comparatively) high calibration of (1 − b)∼0.8.

On a final note, this work was finished within a week of the Planck 2018 CMB likelihood being published, hence the use of the older 2015 data instead. We did check, however, that for the ΛCDM case notably, using this new CMB data yields a slight difference on the preferred values of (1 − b) = 0.622 ± 0.028, which still leaves the level of tension above 6σ (see Fig. 8). Therefore, we do not expect significant changes to the other results of our study and to our conclusions.

thumbnail Fig. 8.

Marginalised posterior distributions for AT − M and (1 − b) calibration parameters, for combination of CMB, SZ and X-ray cluster data in various cosmological models: with or without massive neutrinos, with or without free growth index γ, and free σ8. For the ΛCDM and νγCDM models, the dashed lines show the same posterior obtained when using the Planck 2018 CMB data instead of the 2015 one.


1

Equation (2) is different when the contrast density Δc is expressed with respect to the critical density at redshift z.

2

Planck Collaboration XXIV (2016) provides best-fit values and uncertainties for all four of those parameters, derived from a fit based on X-ray observations. In our analysis, we follow the conventions used by the Planck Collaboration by fixing the β parameter to its best-fit value (=2/3), letting σlnY free with Gaussian priors derived from its fit (=0.075 ± 0.01), and leaving α completely free. We do, however, fix Y* to its best-fit value (log Y* = −0.186) instead of leaving it free, noting that it is degenerate with the bias parameter (1 − b).

3

We derived this line using a subset of the Planck SZ clusters sample, for which both temperature and SZ measurements are available.

Acknowledgments

SI was supported by the European Structural and Investment Fund and the Czech Ministry of Education, Youth and Sports (Project CoGraDS – CZ.02.1.01/0.0/0.0/15_003/0000437). ZS was supported by a grant of excellence from the Agence des Universités Francophones (AFU). This work has been carried out thanks to the support of the OCEVU Labex (ANR-11-LABX-0060) and the A*MIDEX project (ANR-11-IDEX-0001-02) funded by the “Investissements d’Avenir” French government program managed by the ANR. It was also conducted using resources from the IN2P3 Lyon computing center. We thank the referee for their relevant comments and recommendations.

References

  1. Anderson, L., Aubourg, É., Bailey, S., et al. 2014, MNRAS, 441, 24 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arnaud, M., Pratt, G. W., Piffaretti, R., et al. 2010, A&A, 517, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Audren, B., Lesgourgues, J., Benabed, K., & Prunet, S. 2013, JCAP, 1302, 001 [Google Scholar]
  4. Bull, P., Akrami, Y., Adamek, J., et al. 2016, Phys. Dark Univ., 12, 56 [Google Scholar]
  5. Castorina, E., Sefusatti, E., Sheth, R. K., Villaescusa-Navarro, F., & Viel, M. 2014, J. Cosmol. Astropart. Phys., 2, 049 [CrossRef] [Google Scholar]
  6. Clifton, T., Ferreira, P. G., Padilla, A., & Skordis, C. 2012, Phys. Rep., 513, 1 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  7. Costanzi, M., Villaescusa-Navarro, F., Viel, M., et al. 2013, J. Cosmol. Astropart. Phys., 12, 012 [Google Scholar]
  8. Despali, G., Giocoli, C., Angulo, R. E., et al. 2016, MNRAS, 456, 2486 [NASA ADS] [CrossRef] [Google Scholar]
  9. Dvorkin, C., Wyman, M., Rudd, D. H., & Hu, W. 2014, Phys. Rev. D, 90, 083503 [CrossRef] [Google Scholar]
  10. Font-Ribera, A., Kirkby, D., Busca, N., et al. 2014, J. Cosmol. Astropart. Phys., 5, 027 [Google Scholar]
  11. Ilić, S., Blanchard, A., & Douspis, M. 2015, A&A, 582, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Ishak, M. 2019, Liv. Rev. Relativ., 22, 1 [CrossRef] [Google Scholar]
  13. Lesgourgues, J., & Pastor, S. 2012, Adv. High Energy Phys., 2012, 608515 [CrossRef] [Google Scholar]
  14. Lewis, A. 2013, Phys. Rev. D, 87, 103529 [NASA ADS] [CrossRef] [Google Scholar]
  15. Lewis, A., & Bridle, S. 2002, Phys. Rev. D, 66, 103511 [NASA ADS] [CrossRef] [Google Scholar]
  16. L’Huillier, B., Shafieloo, A., & Kim, H. 2018, MNRAS, 476, 3263 [NASA ADS] [CrossRef] [Google Scholar]
  17. Linder, E. V. 2005, Phys. Rev. D, 72, 043529 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  18. Mantz, A. B., von der Linden, A., Allen, S. W., et al. 2015, MNRAS, 446, 2205 [NASA ADS] [CrossRef] [Google Scholar]
  19. Peebles, P. J. E. 1980, The Large-scale Structure of the Universe (Princeton, NJ: Princeton University Press) [Google Scholar]
  20. Planck Collaboration XX. 2014, A&A, 571, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Planck Collaboration I. 2016, A&A, 594, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Planck Collaboration XI. 2016, A&A, 594, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Planck Collaboration XXIV. 2016, A&A, 594, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Roncarelli, M., Carbone, C., & Moscardini, L. 2015, MNRAS, 447, 1761 [NASA ADS] [CrossRef] [Google Scholar]
  25. Rozo, E., Rykoff, E. S., Bartlett, J. G., & Evrard, A. E. 2013, ArXiv e-prints [arXiv:1302.5086] [Google Scholar]
  26. Sakr, Z., Ilić, S., Blanchard, A., Bittar, J., & Farah, W. 2018, A&A, 620, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Salvati, L., Douspis, M., & Aghanim, N. 2018, A&A, 614, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, ApJ, 688, 709 [NASA ADS] [CrossRef] [Google Scholar]
  29. Zarrouk, P., Burtin, E., Gil-Marín, H., et al. 2018, MNRAS, 477, 1639 [NASA ADS] [CrossRef] [Google Scholar]

All Figures

thumbnail Fig. 1.

Confidence contours (68 and 95%) in the Ωm − σ8 plane for various datasets: Planck 2015 and 2018 CMB data, and Planck 2015 SZ cluster counts (using the D16 mass function). Contours are shown for the minimal neutrino case and three neutrinos with free masses (using the CDM prescription for cluster counts). For SZ clusters, the calibration parameter (1 − b) is fixed to its fiducial value of 0.8, with priors added to Ωb from BBN constraints, and to H0 from BAO data. Two additional dashed and dash-dotted contours are shown to illustrate the effects of changing the mass function used (from Despali et al. 2016 to Tinker et al. 2008) and removing the CDM prescription, respectively.

In the text
thumbnail Fig. 2.

Confidence contours (68 and 95%) for cosmological parameters Ωm and σ8, as well as mass calibration parameters (1 − b) and AT − M when combining CMB data with X-ray and/or SZ clusters data for D16 mass function.

In the text
thumbnail Fig. 3.

Confidence contours (68 and 95%) and posterior distributions for neutrino masses ∑mν and mass calibration parameter (1 − b) when combining CMB and SZ cluster data, for D16 mass function (blue). The effect of adding X-ray clusters data is shown in orange, and the effect of a strong prior on (1 − b) centred around the Planck calibration value of 0.8 in green. For comparison, the 1D posterior on neutrino masses from Planck CMB data alone is also shown in black.

In the text
thumbnail Fig. 4.

Confidence contours (68 and 95%) and posterior distributions for growth index γ and mass calibration parameter (1 − b) when combining CMB and SZ cluster data for D16 mass function (blue). The effect of adding X-ray clusters data is shown in orange, and the effect of a strong prior on (1 − b) centred around the Planck calibration value of 0.8 in green.

In the text
thumbnail Fig. 5.

Confidence contours (68 and 95%) and posterior distributions for growth index γ and mass calibration parameter (1 − b) when combining CMB and SZ clusters data for D16 mass function (blue). The effect of adding X-ray clusters data is shown in orange, and the effect of a strong prior on (1 − b) centred around the Planck calibration value of 0.8 in green.

In the text
thumbnail Fig. 6.

Confidence contours (68 and 95%) and posterior distributions for neutrino masses ∑mν and mass calibration parameter (1 − b) when combining CMB and SZ clusters data. The legend and colours are identical to Fig. 5. For comparison, the 1D posterior on neutrino masses from Planck CMB data alone is also shown in black.

In the text
thumbnail Fig. 7.

Confidence contours (68 and 95%) in (1 − b) vs. AT − M plane, for combination of CMB, SZ and X-ray cluster data in various cosmological models, with or without massive neutrinos, with or without free growth index γ, and free σ8. The black dashed line corresponds to the (1 − b)−AT − M relation leading to the same derived (SZ and X-ray) masses for a selected sample of clusters (see text for details).

In the text
thumbnail Fig. 8.

Marginalised posterior distributions for AT − M and (1 − b) calibration parameters, for combination of CMB, SZ and X-ray cluster data in various cosmological models: with or without massive neutrinos, with or without free growth index γ, and free σ8. For the ΛCDM and νγCDM models, the dashed lines show the same posterior obtained when using the Planck 2018 CMB data instead of the 2015 one.

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.