## Cluster counts. II. Tensions, massive neutrinos, and modified gravity

Stéphane Ili´ c 1 2 ; , Ziad Sakr 2 3 ; , and Alain Blanchard 2

- 1 CEICO, Institute of Physics of the Czech Academy of Sciences, Na Slovance 2, Praha 8 Czech Republic
- 2 IRAP, Université de Toulouse, CNRS, CNES, UPS, Toulouse, France
- 3 Université St Joseph; UR EGFEM, Faculty of Sciences, Beirut, Lebanon

Received XXXX-XX-XX; accepted XXXX-XX-XX

## ABSTRACT

The Lambda cold dark matter ( GLYPH&lt;3&gt; CDM) concordance model is very successful at describing our Universe with high accuracy and only a few parameters. Despite its successes, a few tensions persist; most notably, the best-fit GLYPH&lt;3&gt; CDM model, as derived from the Planck cosmic microwave background (CMB) data, largely overpredicts the abundance of Sunyaev-Zel'dovich (SZ) clusters when using their standard mass calibration. Whether this is the sign of an incorrect calibration or the need for new physics remains a matter of debate. In this study, we examined two simple extensions of the standard model and their ability to release the aforementioned tension: massive neutrinos and a simple modified gravity model via a non-standard growth index GLYPH&lt;13&gt; . We used both the Planck CMB power spectra and SZ cluster counts as datasets, alone and in combination with local X-ray clusters. In the case of massive neutrinos, the cluster-mass calibration (1 GLYPH&lt;0&gt; b ) is constrained to 0 585 : + 0 031 : GLYPH&lt;0&gt; 0 037 : (68% limits), more than 5 GLYPH&lt;27&gt; away from its standard value (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8. We found little : correlation between neutrino masses and cluster calibration, corroborating previous conclusions derived from X-ray clusters; massive neutrinos do not alleviate the cluster-CMB tension. With our simple GLYPH&lt;13&gt; model, we found a large correlation between the calibration and the growth index GLYPH&lt;13&gt;; but contrary to local X-ray clusters, SZ clusters are able to break the degeneracy between the two parameters thanks to their extended redshift range. The calibration (1 GLYPH&lt;0&gt; b ) was then constrained to 0 602 : + 0 053 : GLYPH&lt;0&gt; 0 065 : , leading to an interesting constraint on GLYPH&lt;13&gt; = 0 60 : GLYPH&lt;6&gt; 0 13. When both massive neutrinos and modified gravity were allowed, preferred values remained centred on standard : GLYPH&lt;3&gt; CDM values, but a calibration (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8 was allowed (though only at the 2 : GLYPH&lt;27&gt; level) provided P m GLYPH&lt;23&gt; GLYPH&lt;24&gt; 0 34 eV and : GLYPH&lt;13&gt; GLYPH&lt;24&gt; 0 8. We : conclude that massive neutrinos do not relieve the cluster-CMB tension, and that a calibration close to the standard value (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8 : would call for new physics in the gravitational sector.

Key words. Galaxies: clusters: general - large-scale structure of Universe - cosmological parameters - cosmic background radiation

## 1. Introduction

The Lambda cold dark matter ( GLYPH&lt;3&gt; 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 GLYPH&lt;3&gt; 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 GLYPH&lt;3&gt; CDMmodel with unprecedented accuracy.

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 GLYPH&lt;0&gt; b ) = 0 8 in the fiducial case. :

However, some tension is also noted, particularly between the GLYPH&lt;3&gt; CDMcosmological 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 et al. 2014, 2016b). This particular tension appears only when the standard mass calibration (1 GLYPH&lt;0&gt; 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 GLYPH&lt;27&gt; 8 (the presentday amplitude of matter fluctuations) versus GLYPH&lt;10&gt; 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).

A critical aspect of this tension is that the mass of a cluster (including its dark matter halo) is not directly observable. The

Should the aforementioned tension be confirmed, we would be forced to consider extensions or alternatives to the standard GLYPH&lt;3&gt; CDMmodelofcosmology. 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 e GLYPH&lt;11&gt; ects, 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 &amp;Pastor 2012, for a review on the e GLYPH&lt;11&gt; ect of neutrinos in cosmol-

ogy). Alternatively, one could consider replacing the standard cosmological constant with a di GLYPH&lt;11&gt; erent 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 GLYPH&lt;13&gt; model, can reconcile the Planck cluster mass calibration with the Planck CMB cosmology at the expense of a relatively high GLYPH&lt;13&gt; GLYPH&lt;24&gt; 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 GLYPH&lt;0&gt; 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 GLYPH&lt;13&gt; 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. c (2018). We refer the interested reader to those two articles for more details, as we only provide the main steps here.

Wecomputed 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 selfconsistently (see Sakr et al. 2018). We did, however, use the 'critical M 500 ' 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 M 500 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 e GLYPH&lt;11&gt; ect of the CDM prescription in the analysis; as discussed in Sakr et al. (2018), we include its e GLYPH&lt;11&gt; ects for the sake of completeness, even though it led to only marginal di GLYPH&lt;11&gt; erences. As seen in Fig. 1, it led to slightly lower values of GLYPH&lt;27&gt; 8 ; and concurrently high values of GLYPH&lt;10&gt; 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

$$\underset { \text{ct} \text{ CMR} } { \text{rom use} } \ \ f ( a ) = \Omega _ { m } ( a ) ^ { \gamma },$$

where GLYPH&lt;13&gt; , the growth index, is here a free parameter generalising the approximation GLYPH&lt;13&gt; GLYPH&lt;24&gt; 0 545 in standard gravity (Peebles 1980). : The growth rate f is defined as f ( a ) GLYPH&lt;17&gt; d ln G a ( ) = d ln a , where G is the (scale-independent) linear growth factor of matter density perturbations GLYPH&lt;14&gt; m defined as GLYPH&lt;14&gt; m ( a ) = GLYPH&lt;14&gt; m ; 0 G a ( ) = G (0).

In all extensions to the standard model considered here, we assume the non-linear regime would follow the same threshold GLYPH&lt;14&gt; c ; and the same virial density, GLYPH&lt;1&gt; 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

$$3, \text{ while} \quad T = A _ { T - M } ( h \, M _ { \Delta } ) ^ { 2 / 3 } \left ( \frac { \Omega _ { m } \Delta ( z ) } { 1 7 8 } \right ) ^ { 1 / 3 } ( 1 + z ),$$

where AT GLYPH&lt;0&gt; M is the normalisation parameter, GLYPH&lt;1&gt; is the density contrast chosen for the definition of a cluster, expressed with respect to the total background matter density 1 of the Universe at redshift z , and M GLYPH&lt;1&gt; 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 ¯ ) and mass Y M 500

$$\stackrel { \dots } { \text{akr et al.} } \\ \text{icles for } \ E ^ { - \beta } ( z ) \left [ \frac { D _ { A } ^ { 2 } ( z ) \bar { Y } _ { 5 0 0 } } { 1 0 ^ { - 4 } \, \text{Mpc} ^ { 2 } } \right ] = Y _ { * } \left [ \frac { h } { 0. 7 } \right ] ^ { - 2 + \alpha } \left [ \frac { ( 1 - b ) \, M _ { 5 0 0 } } { 6 \times 1 0 ^ { 1 4 } \, M _ { \odot } } \right ] ^ { \alpha },$$

where the quantity DA refers to the angular diameter distance and E z ( ) GLYPH&lt;17&gt; H z ( ) = H 0 . We note here the aforementioned hydrostatic mass bias (1 GLYPH&lt;0&gt; b ), while GLYPH&lt;11&gt; , GLYPH&lt;12&gt; , and Y GLYPH&lt;3&gt; are additional free parameters in the SZ scaling law. A fourth free parameter exists, namely GLYPH&lt;27&gt; ln Y ; it corresponds to the scatter of the log-normal distribution that the measured Y 500 is assumed to follow (and of which the mean ¯ Y 500 is given by Eq. (3)) 2 .

1 Eq. (2) is di GLYPH&lt;11&gt; erent when the contrast density GLYPH&lt;1&gt; c is expressed with respect to the critical density at redshift z .

2 Planck Collaboration et al. (2016b) 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 GLYPH&lt;12&gt; parameter to its best-fit value ( = 2 3), letting = GLYPH&lt;27&gt; ln Y free with Gaussian priors derived from its fit ( = 0 075 : GLYPH&lt;6&gt; 0 01), and leaving : GLYPH&lt;11&gt; completely free. We do, however, fix Y GLYPH&lt;3&gt; to its best-fit value (log Y GLYPH&lt;3&gt; = GLYPH&lt;0&gt; 0 186) instead of leaving it free, : noting that it is degenerate with the bias parameter (1 GLYPH&lt;0&gt; b ).

## 2.2. Datasets and numerical tools

We followed a standard MCMC analysis using two publicly available codes: the CosmoMC package (Lewis &amp; 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 et al. 2016c), 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 c clusters sample (Planck Collaboration et al. 2016b), separated or combined.

The 2015 Planck CMB (Planck Collaboration et al. 2016a) 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 et al. 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 neutric nos are already implemented in both the CLASS and CAMB codes, we implemented there the other GLYPH&lt;3&gt; CDM extension we considered, namely, the GLYPH&lt;13&gt; model. Additionally, we modified our X-ray and SZ clusters likelihood modules accordingly to support the massive neutrinos case and GLYPH&lt;13&gt; 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 GLYPH&lt;3&gt; CDMmodel is illustrated in Fig. 1 in the GLYPH&lt;10&gt; m GLYPH&lt;0&gt; GLYPH&lt;27&gt; 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 GLYPH&lt;0&gt; 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 e GLYPH&lt;11&gt; ect 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 GLYPH&lt;10&gt; m and lower values of GLYPH&lt;27&gt; 8 , elongating the two in the direction of the GLYPH&lt;10&gt; m GLYPH&lt;0&gt; GLYPH&lt;27&gt; 8 degeneracy, with the CDM prescription shifting the contours slightly to higher GLYPH&lt;10&gt; m and lower GLYPH&lt;27&gt; 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 Xray 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 GLYPH&lt;10&gt; b and H 0 to realistic values for cluster contours. We checked that considering wider priors results in more elongated contours in the direction of the GLYPH&lt;10&gt; m GLYPH&lt;0&gt; GLYPH&lt;27&gt; 8 degeneracy, but without relieving the tension.

σ

m

Fig. 1. Confidence contours (68 and 95%) in the GLYPH&lt;10&gt; m GLYPH&lt;0&gt; GLYPH&lt;27&gt; 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 GLYPH&lt;0&gt; b ) is fixed to its fiducial value of 0.8, with priors added to GLYPH&lt;10&gt; b from BBN constraints, and to H 0 from BAO data. Two additional dashed and dash-dotted contours are shown to illustrate the e GLYPH&lt;11&gt; ects of changing the mass function used (from Despali et al. 2016 to Tinker et al. 2008) and removing the CDM prescription, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-31148087-2554-43ea-b6c7-6dec38ecdc17.jpg)

## 3.1. The GLYPH&lt;3&gt; 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 GLYPH&lt;3&gt; 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 GLYPH&lt;0&gt; b ) = 0 58 : GLYPH&lt;6&gt; 0 03 and : note a clear correlation between the (1 GLYPH&lt;0&gt; b ) calibration factor for SZ clusters and the AT GLYPH&lt;0&gt; M calibration for X-ray clusters (at virial mass), indicating the consistency of the constraints obtained from the two types of cluster samples.

## 3.2. Are massive neutrinos a solution?

The possible role of massive neutrinos in alleviating the clusterCMB 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 GLYPH&lt;0&gt; 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 GLYPH&lt;0&gt; b ) = 0 585 : + 0 031 : GLYPH&lt;0&gt; 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 GLYPH&lt;24&gt; 28%) compared

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-40fd611c-45a2-4b89-a2b3-3ca16074e1fe.jpg)

-

m

Fig. 2. Confidence contours (68 and 95%) for cosmological parameters GLYPH&lt;10&gt; m and GLYPH&lt;27&gt; 8 , as well as mass calibration parameters (1 GLYPH&lt;0&gt; b ) and AT GLYPH&lt;0&gt; M when combining CMB data with X-ray and or SZ clusters data for D16 / mass function.

to the CMB-only constraint (black curve). The standard Planck mass calibration ((1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8) is : formally rejected at more than 5 GLYPH&lt;27&gt; . However, for illustration purposes, we show the consequence of putting a strong (Gaussian) prior on (1 GLYPH&lt;0&gt; b ) consistent with the standard Planck calibration, namely (1 GLYPH&lt;0&gt; b ) = 0 8 : GLYPH&lt;6&gt; 0 01 : (green contours). In this case, the favoured neutrino masses are found to be non-zero: P m GLYPH&lt;23&gt; = 0 70 : GLYPH&lt;6&gt; 0 15 eV, but the correspond-: ing best-fit model has a GLYPH&lt;1&gt; GLYPH&lt;31&gt; 2 GLYPH&lt;25&gt; 20 compared to the case without a prior on (1 GLYPH&lt;0&gt; 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 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 6, significantly lower than the standard : calibration (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8. :

## 3.3. Introducing modified gravity

If the standard calibration (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8 is consolidated by ob-: servations, the conclusion of the previous sections leads us to examine alternative theories to GLYPH&lt;3&gt; 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 CMBseen in the standard cosmological picture. Such a lower GLYPH&lt;27&gt; 8 would naturally come from a late-time deviation of the growth of matter fluctuations with respect to GLYPH&lt;3&gt; 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 GLYPH&lt;13&gt; following Eq. (1). We follow the same approach as in

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-47cca5e5-dead-4123-843c-e507fcea0789.jpg)

-

Fig. 3. Confidence contours (68 and 95%) and posterior distributions for neutrino masses P m GLYPH&lt;23&gt; and mass calibration parameter (1 GLYPH&lt;0&gt; b ) when combining CMB and SZ cluster data, for D16 mass function (blue). The e GLYPH&lt;11&gt; ect of adding X-ray clusters data is shown in orange, and the e GLYPH&lt;11&gt; ect of a strong prior on (1 GLYPH&lt;0&gt; 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.

previous sections and perform an MCMC analysis to explore the space of the standard cosmological parameters, adding the calibration (1 GLYPH&lt;0&gt; b ) and the index GLYPH&lt;13&gt; as additional free parameters under the constraints of our data sets. As previously found in Sakr et al. (2018), we expect increasing GLYPH&lt;13&gt; to yield a lower GLYPH&lt;27&gt; 8 ; and thereby a higher calibration (1 GLYPH&lt;0&gt; 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 GLYPH&lt;0&gt; b ) shows a clear correlation with GLYPH&lt;13&gt; , as anticipated. However, a significant di GLYPH&lt;11&gt; erence 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 GLYPH&lt;13&gt; are rejected by the SZ cluster data, resulting in closed contours in the GLYPH&lt;13&gt; versus (1 GLYPH&lt;0&gt; 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 GLYPH&lt;0&gt; b ) = 0 602 : + 0 053 : GLYPH&lt;0&gt; 0 065 : , and the standard calibration (1 GLYPH&lt;0&gt; b ) = 0 8 is rejected at more than 99%, while : GLYPH&lt;13&gt; = 0 60 : GLYPH&lt;6&gt; 0 13 : is consistent with the growth rate in the standard GLYPH&lt;3&gt; CDM model ( GLYPH&lt;24&gt; 0 545). This illustrates the fact that even a simple modified : version of the growth rate cannot easily accommodate the calibration (1 GLYPH&lt;0&gt; b ) = 0 8. Again for illustration purposes, we show the : consequences of putting a strong prior on (1 GLYPH&lt;0&gt; b ) ; consistent with standard Planck calibration ( = 0 8 : GLYPH&lt;6&gt; 0 01, green contours). The : value of GLYPH&lt;13&gt; is then found to be constrained to GLYPH&lt;13&gt; = 0 917 : GLYPH&lt;6&gt; 0 088 : ; consistent with conclusions derived from the X-ray sample alone (Sakr et al. 2018).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-c48e07b2-b224-4edb-8509-69c544edbc2e.jpg)

-

Fig. 4. Confidence contours (68 and 95%) and posterior distributions for growth index GLYPH&lt;13&gt; and mass calibration parameter (1 GLYPH&lt;0&gt; b ) when combining CMB and SZ cluster data for D16 mass function (blue). The e GLYPH&lt;11&gt; ect of adding X-ray clusters data is shown in orange, and the e GLYPH&lt;11&gt; ect of a strong prior on (1 GLYPH&lt;0&gt; 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 GLYPH&lt;13&gt; . As can be seen from Fig. 5 and 6, when combining CMB and SZ clusters, a degeneracy is present between the e GLYPH&lt;11&gt; ects of massive neutrinos and modified gravity. The preferred model remains, however, close to the standard picture: (1 GLYPH&lt;0&gt; b ) = 0 651 : + 0 043 : GLYPH&lt;0&gt; 0 12 : , GLYPH&lt;13&gt; = 0 65 : + 0 13 : GLYPH&lt;0&gt; 0 16 : and P m GLYPH&lt;23&gt; GLYPH&lt;20&gt; 0 197, although an SZ calibration of : (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8 is allowed within the two : GLYPH&lt;27&gt; posterior contours. The addition of X-ray clusters reduces the contours, although in a marginal way. When constraining (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8 with a strong : prior, the preferred parameters are still quite far away from the standard model: GLYPH&lt;13&gt; = 0 834 : + 0 080 : GLYPH&lt;0&gt; 0 094 : and P m GLYPH&lt;23&gt; = 0 34 : + 0 13 : GLYPH&lt;0&gt; 0 15 : eV, but we note that a model with GLYPH&lt;13&gt; GLYPH&lt;24&gt; 0.65 and P m GLYPH&lt;23&gt; GLYPH&lt;24&gt; 0 1 is within the two : GLYPH&lt;27&gt; contours.

## 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 GLYPH&lt;3&gt; CDM picture, but allowed the present-day amplitude of matter fluctuations GLYPH&lt;27&gt; 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 GLYPH&lt;24&gt; 1 (the redshift range containing both the SZ and X-ray cluster samples) follows the standard GLYPH&lt;3&gt; CDM growth rate, but modified at earlier times through a process about which we do not assume anything. In practice, we rescale the GLYPH&lt;3&gt; CDM matter power spectrum for z 2 [0 ; 1] by a constant factor, so that we obtain the desired input GLYPH&lt;27&gt; 8 while keeping the same value of the primordial As . Such

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-bc5fa3a8-6803-44af-835e-e92dfd51acea.jpg)

-

Fig. 5. Confidence contours (68 and 95%) and posterior distributions for growth index GLYPH&lt;13&gt; and mass calibration parameter (1 GLYPH&lt;0&gt; b ) when combining CMB and SZ clusters data for D16 mass function (blue). The e GLYPH&lt;11&gt; ect of adding X-ray clusters data is shown in orange, and the e GLYPH&lt;11&gt; ect of a strong prior on (1 GLYPH&lt;0&gt; b ) centred around the Planck calibration value of 0.8 in green.

an approach allows for (comparatively) more freedom than the GLYPH&lt;13&gt; 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 GLYPH&lt;0&gt; b ) versus AT GLYPH&lt;0&gt; 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 ( GLYPH&lt;3&gt; CDM + free GLYPH&lt;27&gt; 8 at present epoch). We observe that all contours lay along the correspondence line 3 (dashed black) between (1 GLYPH&lt;0&gt; b ) and AT GLYPH&lt;0&gt; M ; and that in all cases the preferred SZ mass calibration (1 GLYPH&lt;0&gt; b ) remains around 0.6 (see 1D posteriors in Fig. 8) even in the case where GLYPH&lt;27&gt; 8 is treated as a free parameter. In the latter case, a value of (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 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 GLYPH&lt;27&gt; 8 .

## 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 GLYPH&lt;0&gt; b ) for SZ clusters, and AT GLYPH&lt;0&gt; M for X-ray clusters) as free parameters without any priors. Neither massive neutrinos nor a modified gravity model (parametrised by the GLYPH&lt;13&gt; growth index) allow us to reconcile the Planck data (CMB and SZ clusters) with the standard SZ mass calibration value of

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

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-4284358c-a366-4012-a0e2-8212bd154dd3.jpg)

-

Fig. 6. Confidence contours (68 and 95%) and posterior distributions for neutrino masses P m GLYPH&lt;23&gt; and mass calibration parameter (1 GLYPH&lt;0&gt; 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.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-9eecc19d-3529-4f57-a024-09e34359020a.jpg)

Fig. 7. Confidence contours (68 and 95%) in (1 GLYPH&lt;0&gt; b ) vs. AT GLYPH&lt;0&gt; 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 GLYPH&lt;13&gt; , and free GLYPH&lt;27&gt; 8 . The black dashed line corresponds to the (1 GLYPH&lt;0&gt; b ) -AT GLYPH&lt;0&gt; M relation leading to the same derived (SZ and X-ray) masses for a selected sample of clusters (see text for details).

(1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8. Consequently, this implies that there is no sim-: ple 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 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 8 value). Despite this observation, we ex-: plore 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 GLYPH&lt;27&gt; 8 ). Nevertheless, the SZ cluster counts com-

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00163v2-figure-af6596a1-ce8c-4497-ac2f-c89e8288afca.jpg)

-

Fig. 8. Marginalised posterior distributions for AT GLYPH&lt;0&gt; M and (1 GLYPH&lt;0&gt; 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 GLYPH&lt;13&gt; , and free GLYPH&lt;27&gt; 8 . For the GLYPH&lt;3&gt; CDMand GLYPH&lt;23&gt;GLYPH&lt;13&gt; CDM models, the dashed lines show the same posterior obtained when using the Planck 2018 CMB data instead of the 2015 one.

bined with CMB data still prefer a low calibration ((1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 0 6) : in all considered models. We therefore find it striking that no cosmological model appears to prefer the (comparatively) high calibration of (1 GLYPH&lt;0&gt; b ) GLYPH&lt;24&gt; 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 GLYPH&lt;3&gt; CDM case notably, using this new CMB data yields a slight di GLYPH&lt;11&gt; erence on the preferred values of (1 GLYPH&lt;0&gt; b ) = 0 622 : GLYPH&lt;6&gt; 0 028 : ; which still leaves the level of tension above 6 GLYPH&lt;27&gt; (see Fig. 8). Therefore, we do not expect significant changes to the other results of our study and to our conclusions.

Acknowledgements. 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-11LABX-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

Anderson, L., Aubourg, É., Bailey, S., et al. 2014, MNRAS, 441, 24 Arnaud, M., Pratt, G. W., Pi GLYPH&lt;11&gt; aretti, R., et al. 2010, A&amp;A, 517, A92 Audren, B., Lesgourgues, J., Benabed, K., &amp; Prunet, S. 2013, JCAP, 1302, 001 Bull, P., Akrami, Y., Adamek, J., et al. 2016, Physics of the Dark Universe, 12, 56

Castorina, E., Sefusatti, E., Sheth, R. K., Villaescusa-Navarro, F., &amp; Viel, M. 2014, J. Cosmology Astropart. Phys., 2, 049

Clifton, T., Ferreira, P. G., Padilla, A., &amp; Skordis, C. 2012, Phys. Rep., 513, 1 Costanzi, M., Villaescusa-Navarro, F., Viel, M., et al. 2013, J. Cosmology Astropart. Phys., 12, 012

Despali, G., Giocoli, C., Angulo, R. E., et al. 2016, MNRAS, 456, 2486 Dvorkin, C., Wyman, M., Rudd, D. H., &amp; Hu, W. 2014, Phys. Rev. D, 90, 083503 Font-Ribera, A., Kirkby, D., Busca, N., et al. 2014, J. Cosmology Astropart. Phys., 5, 027

Ili´, S., Blanchard, A., &amp; Douspis, M. 2015, A&amp;A, 582, A79

c

Ishak, M. 2019, Living Reviews in Relativity, 22, 1

Lesgourgues, J. &amp; Pastor, S. 2012, ArXiv e-prints

Lewis, A. 2013, Phys. Rev. D, 87, 103529

Lewis, A. &amp; Bridle, S. 2002, Phys. Rev. D, 66, 103511

L'Huillier, B., Shafieloo, A., &amp; Kim, H. 2018, MNRAS, 476, 3263

Linder, E. V. 2005, Phys. Rev. D, 72, 043529

Mantz, A. B., von der Linden, A., Allen, S. W., et al. 2015, MNRAS, 446, 2205 Peebles, P. J. E. 1980, The large-scale structure of the universe (Princeton University Press)

Planck Collaboration, Adam, R., Ade, P. A. R., et al. 2016a, A&amp;A, 594, A1 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2014, A&amp;A, 571, A20 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016b, A&amp;A, 594, A24 Planck Collaboration, Aghanim, N., Arnaud, M., et al. 2016c, A&amp;A, 594, A11 Roncarelli, M., Carbone, C., &amp; Moscardini, L. 2015, Monthly Notices of the Royal Astronomical Society, 447, 1761

Rozo, E., Ryko GLYPH&lt;11&gt; , E. S., Bartlett, J. G., &amp; Evrard, A. E. 2013, arXiv e-prints, arXiv:1302.5086

Sakr, Z., Ili´, S., Blanchard, A., Bittar, J., &amp; Farah, W. 2018, A&amp;A, 620, A78 c

Salvati, L., Douspis, M., &amp; Aghanim, N. 2018, A&amp;A, 614, A13 Tinker, J., Kravtsov, A. V., Klypin, A., et al. 2008, The Astrophysical Journal, 688, 709

Zarrouk, P., Burtin, E., Gil-Marín, H., et al. 2018, Monthly Notices of the Royal Astronomical Society, 477, 1639
Stéphane Ilić Ziad Sakr Alain Blanchard
Cluster counts. II. Tensions, massive neutrinos, and modified gravity
The $\Lambda$CDM concordance model is very successful at describing our
Universe with high accuracy and few parameters. Despite its successes, a few
tensions persist; most notably, the best-fit $\Lambda$CDM model, as derived
from the Planck CMB data, largely overpredicts the abundance of SZ clusters
when using their standard mass calibration. Whether this is a sign of an
incorrect calibration or the need for new physics remains a matter of debate.
Here we examined two simple extensions of the standard model and their ability
to release this tension: massive neutrinos and a simple modified gravity model
via a non-standard growth index $\gamma$. We used both the Planck CMB and SZ
cluster counts as datasets, with or without local X-ray clusters. In the case
of massive neutrinos, the SZ calibration $(1-b)$ is constrained to
$0.59^{+0.03}_{-0.04}$ (68\%), more than 5$\sigma$ away from its standard value
$\sim0.8$. We found little correlation between $\sum m_\nu$ and $(1-b)$,
corroborating previous conclusions derived from X-ray clusters; massive
neutrinos do not alleviate the cluster-CMB tension. With our simple $\gamma$
model, we found a large correlation between calibration and growth index but
contrary to local X-ray clusters, SZ clusters are able to break the degeneracy
between the two thanks to their extended $z$ range. The calibration $(1-b)$ was
then constrained to $0.60^{+0.05}_{-0.07}$, leading to an interesting
constraint on $\gamma=0.60\pm 0.13$. When both massive neutrinos and modified
gravity were allowed, preferred values remained centred on standard
$\Lambda$CDM values, but $(1-b)\sim0.8$ was allowed (though only at the
$2\sigma$ level) provided $\sum m_\nu\sim0.34 $ eV and $\gamma\sim0.8$. We
conclude that massive neutrinos do not relieve the cluster-CMB tension and that
a calibration close to the standard value $0.8$ would call for new physics in
the gravitational sector.