# Cluster gas fraction as a test of gravity

###### Abstract

We propose a new cosmological test of gravity, by using the observed mass fraction of X-ray emitting gas in massive galaxy clusters. The cluster gas fraction, believed to be a fair sample of the average baryon fraction in the Universe, is a well-understood observable, which has previously mainly been used to constrain background cosmology. In some modified gravity models, such as gravity, gas temperature in a massive cluster is determined by the effective mass of that cluster, which can be larger than its true mass. On the other hand, X-ray luminosity is determined by the true gas density, which in both modified gravity and CDM models depends mainly on and hence the true total cluster mass. As a result, the standard practice of combining gas temperatures and X-ray surface brightnesses of clusters to infer their gas fractions can, in modified gravity models, lead to a larger – in gravity this can be larger – value of than that inferred from other observations such as the CMB. A quick calculation shows that the Hu-Sawicki model with is in tension with the gas fraction data of the 42 clusters analysed by Allen et al. (2008). We also discuss the implications for other modified gravity models.

###### keywords:

^{†}

^{†}pagerange: Cluster gas fraction as a test of gravity–Cluster gas fraction as a test of gravity

^{†}

^{†}pubyear: 2015

## 1 Introduction

In recent years, attempts to understand the origin of the accelerated cosmic expansion (Riess et al., 1998; Perlmutter et al., 1999) have led to a large number of theoretical models (Copeland et al., 2006). Apart from the current standard model, which assumes that the acceleration is caused by a cosmological constant (hence the name -cold-dark-matter, or CDM), these model can be roughly put in two categories: dark energy – which replaces by some dynamical field, and modified gravity – which assumes that there is no exotic matter species beyond the standard CDM model but gravity is not described by General Relativity (GR) on cosmological scales (see, e.g., Clifton et al., 2012; Joyce et al., 2015, for some recent reviews). Some of the models in the latter class, such as the chameleon theory (Khoury & Weltman, 2004), of which the well-known gravity (Carroll et al., 2005) is an example, have been active research topics in recent years.

Ultimately, any new cosmological model or theory of gravity should be put to rigorous tests against observational data. For this reason, a number of tests have been proposed or applied in the past to examine the viability of the models (see, e.g., Lombriser, 2014; De Martino et al., 2015, for some recent reviews in gravity). The present paper shall follow the same line to propose a test using observations of galaxy clusters.

Galaxy clusters are the largest gravitationally bound and virialised objects in our Universe. The most massive clusters observed today typically have masses in the range of , of which the dominant component is dark matter. These are homes to galaxies, stars and eventually lives, which together hold the vast majority of the information that can be extracted from cosmological and astrophysical observations. In dark energy or modified gravity theories, the different cosmic expansion histories and gravitational laws between particles can have sizeable effects on how the clusters form and evolve. Schmidt, Vikhlinin & Hu (2009); Cataneo et al. (2014), based on this observation, have placed constraints on gravity using cluster abundance data. In theories such as gravity, massive and massless particles feel different strengths of gravity, thus allowing these theories to be constrained by comparing the so-called dynamical and lensing masses of clusters (Schmidt, 2010; Zhao, Li & Koyama, 2011). Combining with lensing observations, Terukina et al. (2014); Wilcox et al. (2015) obtained even stronger constraints on gravity.

Given the abundant information associated with these rich objects, one expects that they will provide a wealth of other potential tests of new cosmological models. The arrival of the era of precision cosmology lends this perspective both more interest and more support. In this paper, we propose to utilise the observationally inferred mass fraction of hot X-ray emitting gas in galaxy clusters as a new test of gravity.

The gas fraction in clusters, , is a well established and understood observable in the standard cosmological model, which can be used to place strong constraints on background cosmology (see, e.g., Allen et al., 2004, 2008, and following-up works, for some examples). The basic assumptions (or approximations) are (i) clusters are in hydrostatic equilibrium between thermal pressure and gravity and (ii) as the largest objects in the Universe, the cluster baryon fraction, dominantly contributed by gas, is a faithful representation of the cosmological average baryon fraction , in which and are respectively the fractional mass density of baryons and all matter (White et al., 1993; Eke, Navarro & Frenk, 1998). Using (i), one can find, from the observed X-ray temperature and surface brightness profiles, the mass profiles of baryons and all matter inside a cluster, and consequently – the profile of gas fraction. Combined with (ii), this can have a say about provided that is measured elsewhere (e.g., from big bang nucleosynthesis or the cosmic microwave background (CMB)).

In modified gravity theories, the hydrostatic equilibrium inside clusters is changed as a result of the different law of gravity. Hence, a cluster can have a higher dynamical mass, with the same baryonic mass inside, leading to a that is lower than the true which is related to . If a cosmologist wishes to infer from in a modified gravity universe, a correction has to be done to the end result, which can lead to inconsistency with other observational determinations of , and therefore a constraint on the gravity theory. We shall demonstrate the potential constraints from this test using gravity as an example.

This paper is organised as following: in Sect. 2 we briefly describe the gravity theory and its equations which will be used in the discussion below. In Sect. 3 we give a more detailed account of the physics related to the gas fraction test described above. Then we present a numerical example in Sect. 4, which shows how current data of cluster gas fraction alone can give powerful constraints on gravity. We discuss the results and their implications in Sect. 5.

## 2 The gravity theory

This section is devoted to a quick overview of gravity. It will be kept brief and only include essential equations, given that there are already many papers in the literature covering this topic.

gravity (Carroll et al., 2005) is a simple generalisation of GR, by replacing the Ricci scalar, , in the Einstein-Hilbert action with an algebraic function

(1) |

where is the reduced Planck mass, with being Newton’s constant, is the determinant of the metric and is the Lagrangian density for (normal plus dark) matter fields. The model is defined by specifying the functional form of .

The action in Eq. (1) leads to a modified Einstein equation

(2) |

in which , are respectively the Einstein and Ricci tensors, , the covariant derivative compatible to the metric , and is the matter energy momentum tensor. Eq. (2) can be considered as the standard Einstein equation of GR with an extra scalar field, , whose dynamics is governed by

(3) |

where is the mass density of baryons and dark matter. As we are interested in late times, photons and neutrinos will be neglected.

On scales well inside the Hubble radius, and for the models to be considered, it is safe to work with the quasi-static approximation (Bose, Hellwing & Li., 2015), in which the scalar equation becomes

(4) |

where denotes the three dimensional gradient, is the scale factor, and an overbar takes the background value of a quantity. Notice that can be expressed as a function of by inverting .

Similarly, the modified Poisson equation in this limit reads as

(5) |

where is the Newtonian potential.

Eq. (4) implies two limits of the behaviour of gravity:

(i) When is small, or more accurately, when , it recovers the well-known GR solution , and so Eq. (5) reduces to GR as well. This is the chameleon (Khoury & Weltman, 2004) regime which any viable model must be in to pass the stringent solar system and terrestrial tests of gravity.

(ii) When , the second term on the rhs of Eq. (5) is negligible compared with the first term, so that we have a gravity that is stronger than in GR. This is ususally known as the non-chameleon, or unscreened, regime.

It is evident that the unscreened regime mostly happens where is shallow, or in extensive regions of low density. On large scales, matter density is close to the cosmological average, and so the total gravity in gravity is enhanced within scales comparable to the Compton wavelength of the scalar field (which in most models of interest is in the range Mpc). This naturally leads to an enhanced large-scale structure formation, and features such as over-abundant and more massive galaxy clusters - a topic which has been extensively studied previously. This will also be the topic that we focus on in this paper.

## 3 Cluster gas fraction

Galaxy clusters are the largest bound objects in the Universe, whose masses are dominated by the dark matter component, with the baryonic masses dominated by X-ray emitting intracluster gas, which is heated to temperatures of the order of keV during virialisation. It is the mass fraction of this gas component that we will employ to test the theory of gravity here.

In this section, we shall first give a brief overview of how the baryon fraction can be estimated observationally, and how it can be used to constrain cosmological models and their parameters. Then, we will discuss how this process might be affected if the underlying theory of gravity is modified. For simplicity, we shall neglect other baryonic components than the intracluster gas in our analysis unless otherwise stated.

### 3.1 The standard CDM model

In the standard cosmological scenario, halo density profiles can be universally described by the Navarro, Frenk & White (1997, NFW) fitting formula, which is often expressed as

(6) |

in which is the halo mass density as a function of the distance, , from the halo’s centre, is a characteristic density and is the scale radius. We shall assume that the halo is spherically symmetric and well relaxed throughout the analysis, unless otherwise stated.

The mass of the halo can be obtained by integrating the NFW profile from to , in which is the edge of the halo and is defined as the radius within which the average mass density is , with the critical density at the redshift when the halo is identified. This leads to

(7) |

Observationally, the total and baryonic masses of a cluster can be obtained by measuring its X-ray surface brightness profile, and the temperature profile of its X-ray gas. For a dynamically relaxed system that consists of dark matter and baryonic gas, a hydrostatic equilibrium can be achieved, which satisfies the following equation

(8) |

in which is the total mass of dark matter and gas within radius , and , are respectively the density and pressure of the gas at . For simplicity, we neglect non-thermal pressure in our discussion (the effects of non-thermal pressure, however, are taken into account in the error budget when modelling the relation between and , cf. Eq. (19) below).

For an ideal thermal gas, its pressure and density are related to its temperature, , as

(9) |

where is the Boltzmann constant, and the mass and number densities of the gas particles are connected by , where is the proton mass and the mean molecular weight. Applying Eqs. (8, 9), we obtain

(10) |

which, if evaluated at the halo edge (), gives

(11) |

In the mean time, the X-ray emission of galaxy clusters is produced mainly by thermal bremsstrahlung radiation, leading to a X-ray surface brightness profile that also depends on the gas density and temperature profiles, and . In terms of gas density, the gas mass can be expressed as

(12) |

where is the gas density at the cluster centre, and we have assumed that with some function describing the profile and a characteristic scale. In the isothermal model, for example, we have

(13) |

where is a dimensionless constant.

The bolometric luminosity is given by (e.g., Sasaki, 1996)

(14) |

where we have neglected the proportionality coefficients which are irrelevant here.

Combining Eqs. (12, 14), we have

(15) |

where and are respectively the luminosity and angular diameter distances and are related by . Here we have used where is the X-ray flux, and , in which denotes the angle spanned by at redshift . Note that (or equivalently the surface brightness) and are the observed quantities in this description. In real observations, one has the surface brightness and temperature profiles, or equivalently and , using which Eqs. (10, 14) can be solved simultaneously to find and (note that the innermost regions of clusters are often excluded due to complicated processes such as cooling flow. For example, the study of Allen et al. (2008) uses only the data within to measure ). These provide the necessary information to find

(16) |

Often in observations, people quote the value of at as the cluster gas fraction. We note in passing that here enters the picture only through its gravitational effect on gas particles, as this fact is important for the discussion of modified gravity below.

To see how this can be used to constrain background cosmology, we note that, once the right-hand side of Eq. (11) is known by observations, we also have a fixed numerical value of the left-hand side and hence have the relation

(17) |

Eqs. (15, 17) imply that the measured cluster gas fraction depends on redshift and angular diameter distance in the following specific way:

(18) |

As mentioned earlier, we expect the cluster gas fraction to be a reasonably fair sample of the mean cosmological baryon fraction and is therefore roughly independent of for massive clusters at low . This is a reasonable assumption which should hold regardless of the cosmological model/parameters. The measurement of the apparent gas fraction, as described above, involves the angular diameter distance and is indeed dependent on both the cosmological model and its parameters. Therefore, with incorrect cosmological models or parameters, the constancy of the true gas fraction (hereafter ) is not guaranteed to be reflected in the observed value . This provides a powerful test (Allen et al., 2004) of background cosmology and can be used to constrain cosmological parameters.

Were galaxy clusters perfectly fair samples of the average matter components in the Universe, their baryon fraction would just be . The true situation, however, is more complicated. To take the complexities into account, Allen et al. (2008), improving on the earlier work of Allen et al. (2004), propose the following model of the relation between cluster gas fractions and ,

(19) |

in which:

is a constant accounting for systematic effects such as the calibration of instrument and X-ray modelling – it is assumed to be in Allen et al. (2008);

models the non-thermal pressure support in galaxy clusters which can cause a bias in the estimate of of about %;

is the so-called ’depletion’ factor which is inspired by the observation that the baryon fraction at in non-radiative simulations (Eke, Navarro & Frenk, 1998) is actually smaller than , with and small indicating a weak redshift evolution below ;

accounts for the fact that a small fraction of baryons can be in the form of stars, with and describing its redshift evolution.

### 3.2 Modified gravity scenarios

In many modified gravity theories, including gravity, the way in which the trajectories of massive test bodies – e.g., galaxies, stars and gas particles – respond to the underlying matter distribution is different. This change of the dynamics of test bodies is sometimes described as the change of the dynamical mass of matter. Massless particles, such as photons, could behave differently: in some theories, such as the Galileon model (Nicolis et al., 2009; Deffayet et al., 2009) and the K-mouflage model (Brax & Valageas., 2014a, b), photons can also feel a different mass of matter, but in other models, for example gravity and the Dvali, Gabadadze & Porrati (2000, DGP) model, photon trajectories depend on the matter distribution in essentially the same way as in GR, since the conformal coupling does not affect geodesics of massless particles. For distinction, the mass felt by photons is usually called the lensing mass. The differences in the dynamical and lensing masses of galaxy clusters have been used to constrain gravity in, e.g., Terukina et al. (2014); Wilcox et al. (2015)

As mentioned above, gas particles, like dark matter particles and galaxies, do feel the dynamical mass of a cluster. In gravity, the same cluster can have a dynamical mass times its value in GR. This maximum enhance factor of 1/3, however, is not necessarily realised in all clusters, because of the chameleon screening (Khoury & Weltman, 2004). The screening helps to reduce the difference between the dynamical and lensing masses, especially for more massive clusters. Consequently, constraints relying on the dynamical masses of clusters are in general weaker than those coming from astrophysical considerations. Nevertheless, they have cleaner physics than that of astrophysical observables – which can often depend on whether the considered astrophysical system lives inside a screened cluster – and are amongst the tightest constraints obtained using cosmological data (Terukina et al., 2014).

Evidently, without a reliable measurement of the lensing mass, it is difficult to tell whether one observes a cluster of true mass in GR, or one with a smaller mass in gravity, since both have the same dynamical masses. This has motivated He et al. (2015a) to propose the concept of effective haloes. Briefly speaking, the idea is to redefine the right-hand side of the modified Poisson equation in gravity, Eq. (5), so that it can be rewritten as

(20) |

In the above, is the Newtonian potential that determines the dynamics of massive bodies and is the density perturbation of (dark plus baryonic) matter. is the effective density field, with which the Poisson equation takes exactly the same form as in GR [cf. Eq. (20)]. In this way, the complicated new physics in gravity is absorbed into , and with that solved (e.g., in numerical simulations) one can in principle proceed assuming GR as the true theory of gravity. He & Li (2015b), for example, show with hydrodynamical simulations that cluster gas temperatures depend only on the masses of the corresponding effective haloes, and that with certain rescaling depending on the effective halo mass the scaling relations – such as the - relation with the X-ray luminosity – in gravity can be derived reliably using existing knowledge of GR.

In the left panel of Fig. 1, we present the gas temperature profiles for standard haloes in GR and effective haloes in gravity, both in the mass bin . Though there are differences in the inner regions – which could be due to different halo density profiles or screening – we notice that beyond kpc the two agree very well. He & Li (2015b) find that the average gas temperatures in the two also show very good agreement, and indeed the temperature-mass relation is barely distinguishable in the two models, provided that effective haloes are used in gravity.

The fact that the cluster gas temperatures depend on the mass of the effective haloes is as expected, since for relaxed systems the virial temperature depends on the Newtonian potential, which does not distinguish between standard (GR) and effective ( gravity) density fields. For a polytropic gas with an equation of state in which the constant , the hydrostatic equation implies that (see, e.g., Mo, van den Bosch & White, 2011) the temperature can be analytically expressed as a function of the potential .

Let us now consider two haloes, one identified in the standard dark matter field in a model with GR as the gravity theory, another from the effective density field in gravity. The profiles of the two haloes are the same so that Eq. (20) sees no difference in them. Since gravity only enters the picture through Eq. (10), we make the following two observations/predictions:

the gas temperature profiles are the same in these two haloes;

the two sides of the spherical hydrostatic equation, Eq. (10), are the same for the two haloes.

Since and in Eq. (10) are the same, we conclude that is also the same in the haloes. This, however, does not necessarily mean that the haloes have identical gas density profiles, because we can rescale by a constant factor without changing . To confirm this, in the middle panel of Fig. 1 we compare the gas density profile of effective haloes with that of CDM haloes of the same mass, and find that the two show a constant shift by 1/4 beyond kpc. Note that the simulated haloes in the plot do not have perfectly identical total – standard or effective – mass profiles, which is why in the middle panel of Fig. 1 the red crosses and red dashed line do not agree on scales below kpc. However, as mentioned above Eq. (22), in real observations, such innermost regions are not used in the determination of anyway.

As a result, to obtain the gas density profile, we need further, independent, information to fix its normalisation – as opposed to its shape – which brings us back to the measurements of cluster X-ray surface brightness. Inspecting the equation for the cluster X-ray luminosity, Eq. (14), we notice that the luminosity density (i.e., the integrand) depends on (i) the physical gas density , and (ii) the gas temperature which, as we have seen above, depends on the total mass of the effective halo. Consequently, should the physical gas densities be the same for effective and CDM haloes of the same mass, there would be no difference in their X-ray surface brightness profiles.

However, despite the standard (GR) and effective ( gravity) haloes above having identical gas temperature and halo mass, their actual (physical, or lensing) masses are different, and it is important to remember that characterises the amount of the actual mass in the Universe. If, as we have assumed so far, the gas fraction in clusters is a fair sample of the cosmological value, it would be the ratio of the gas mass and actual halo’s mass that satisfies Eq. (19). Gas fractions inferred observationally, in the way described in the previous subsection, are in fact the ratio of the gas mass and that of the effective halo. If we denote the ratio of the effective and actual masses of a halo by , then , depending on the actual mass and environment of the halo, its redshift, as well as the model parameters^{1}^{1}1More accurately speaking, effective haloes are identified from the effective density field, in Eq. (20), and standard haloes are identified from the physical density field . They do not necessarily share the same physical particles. Here, for simplicity, when talking about the effective and actual masses of some halo, we mean the masses of the effective and standard haloes that would be considered as matched haloes in the two catalogues.. Here, as we are interested in the most massive clusters, with , we can for simplicity neglect the impacts of the halo’s environment, so that mainly depends on , i.e., , for a given halo mass.

Consider the extreme case in which as example. The apparent inferred from X-ray cluster observations would be

(21) |

This is also confirmed by hydro simulations, as shown in the right panel of Fig. 1. In that plot, we see that if haloes are defined using their actual masses in gravity, then they share the same as CDM haloes of the same mass (black solid vs. blue dotted lines). For effective haloes, on the other hand, their profiles are a constant downward shift by from the results of CDM haloes of the same masses (red dashed line vs. red crosses), which is just what Eq. (21) predicts.

It is worthwhile to pause for a moment and try to understand the physics behind the behaviour of Fig. 1. It may seem surprising that, although the gas density profiles in CDM and standard haloes are significantly different within kpc, their gas fraction profiles are very close to each other. This suggests that the two also have different dark matter (or total) mass profiles to cancel the differences in . In other words, follows , which is the physical dark matter mass density, in the same way in CDM and standard haloes, even though they have different mass and even more different potential profiles. To understand this, we note that in Eq. (10) the on the left-hand side and on the right-hand side are both modified in standard haloes. Indeed, if one assumes hydrostatic equilibrium and , then the gas density profile can be written as (Komatsu & Seljak, 2001):

(22) |

in which . Eq. (22) is derived under the assumption of self-similarity of gas density profiles, but the key point therein, that the modified gravity effects on and can be cancelled out, is not affected by this assumption. If this cancellation happens, then the gas density profile is determined by the total mass profile under the assumption of hydrostatic equilibrium, regardless of the theory of gravity^{2}^{2}2Note in Eq. (22) it is that is determined by . The normalisation of will then be fixed by the total gas fraction inside the halo..

In clusters, gas is heated by accretion shocks during the assembly of the halo, a process which involves the conversions of energy from gravitational to kinetic (that of the cold accreted gas) and then to thermal (via shocks). Assuming a complete thermalisation, the post-shock gas temperature is proportional to , with the infall speed of the accreted gas (e.g., Mo, van den Bosch & White, 2011). Consequently, energy conservation implies that the final gas temperatures in the central regions will be affected in the same way as of the cold gas and hence in modified gravity. Of course, this is only an approximation, and the cancellation of the effects of modified gravity on and depend on various factors including the screening and formation history of a cluster, which is not expected to be complete. However, Fig. 1 suggests that it works pretty well for the haloes we use here. We checked explicitly that it works slightly less well for more massive haloes, for which the agreement between the in CDM and standard haloes is slightly less perfect – this may be because those haloes became unscreened only very recently.

The argument above in theory also applies to effective haloes, for which is the same as in GR, but the effects of modified gravity are incorporated in . However, in the effective halo case the normalisation is different because of the different total gas fraction (see footnote 2) – although the shape is the same – hence the nearly constant rescaling of the red dashed curve compared with the black solid one in Fig. 1.

Coming back to the discussion prior to the previous three paragraphs, our result suggests two possible tests of gravity:

(i) If an observer actually lives in a universe shaped by gravity, then the true cluster gas fraction is given by . Assuming that Eq. (19) still holds for , the observer will need to do the following transformation to get the true :

(23) |

As depends only on the actual observational data, the observer will obtain the same value as an observer in a standard GR universe would do. The resulting might then be too large to be compatible with other constraints, such as the one from the CMB.

(ii) Alternatively, if one takes the measured by other probes as the true value and starts from there, then Eq. (21) implies that the observed cluster gas fraction will be smaller than what the CDM model and simulations predict. Because of the time dependence of (see above), if the model parameters happen to take the values for to evolve from 1 to between and the present for the clusters of interest, there may also be an apparent decrease of as decreases, by a maximum of .

## 4 Numerical Examples

In this section, we use a simplified example to illustrate the power of the cluster gas fraction test proposed above. For this, we will use the gas fraction data of the 42 clusters studied by Allen et al. (2008, Table 3). As described above, these data are obtained by fitting the gas temperature and X-ray surface brightness profiles of these clusters simultaneously, assuming NFW profiles for the total mass in clusters. We have also found that, in the context of gravity, as long as we use effective haloes, the dynamics of gas particles can be calculated using standard gravity theory. Therefore, in this work we can directly take the data of Allen et al. (2008) as , bearing in mind that the cluster mass inferred therein would be the effective mass and therefore can be different from for unscreened clusters, cf. Eq. (21).

To obtain an estimation of the mean and standard deviation of from each cluster, random samples of size are drawn for each parameter or data: and . Of these, is taken, for a given cluster, from Table 3 of Allen et al. (2008), and is assumed to satisfy a Gaussian distribution with mean and standard deviation given by Allen et al. (2008). The other parameters and their distributions are shown in Table 1. We therefore obtain realisations of , from which its mean and standard deviation can be calculated. This procedure is repeated for all 42 clusters.

param | physical effect described | mean stddev | prior |
---|---|---|---|

overall calibration | Gaussian | ||

non-thermal pressure | Uniform | ||

gas bias: normalisation | Uniform | ||

gas bias: evolution | Uniform | ||

stellar fraction: normalisation | Gaussian | ||

stellar fraction: evolution | Uniform |

The estimation of the effects of modified gravity, i.e., the factor , is more complicated, since it depends on the cluster mass, density profile, environment, redshift, as well as the parameters. Because the main purpose of this paper is to illustrate the basic idea, we shall leave a full analysis using real cluster data for future work, and instead adopt a simplified modelling. The cluster masses are assumed to be the same, with , for all 42 clusters, since this is a typical value for massive X-ray clusters. The cluster’s radius ( or ) is taken to be Mpc and its concentration parameter, (or ) is . The cluster is assumed to live on the cosmological background, so that the ratio of its effective and actual masses can be approximated as (see, e.g., Li, Zhao & Koyama, 2012)

(24) |

where is the background value of at redshift , and

(25) |

is the Newtonian potential at the edge of an NFW halo. We adopt the model by Hu & Sawicki (2007) with , for which

(26) |

The values and are taken from the latest results of Planck Collaboration, Ade et al. (2015). As a result, the physics of modified gravity is completely governed by , which is the present-day value of . Once this is specified, we can obtain , and therefore infer given and of a cluster.

In Fig. 2 we show the result obtained from for 4 different cases: standard CDM (upper left), and gravity models with (upper right), (lower left) and (lower right). For comparison, we have also, in each panel, plotted the mean value (solid line) and 1 confidence level (dotted) of from Planck Collaboration, Ade et al. (2015)^{3}^{3}3Note that the models studied here have practically identical CMB power spectra as the CDM model with the same and . As a result, using CMB data only, the constraints on cosmological parameters such as would be the same in all these models. Because of this, the CMB constraints are less model-dependent.. The results from the 42 clusters, with 1 errors, are shown as symbols.

A quick naked-eye inspection shows that the method and the CMB observation give compatible if one assumes the CDM paradigm (upper left). The model with (lower right), on the other hand, leads to a significantly higher value of than what CMB says, and is therefore inconsistent. The other two cases are more interesting: for (lower left), increases to at , while for (upper right) only increases to at . In both cases, however, the inferred values of are still substantially larger than the Planck result, especially for the low- clusters. This shows that cluster gas fraction can be a potentially powerful test of gravity, using X-ray observations only. In such tests, lensing data can be a useful addition, but is not necessary.

The test can be done in an alternative way. For this, we assume the value of obtained from CMB observations, and check what value of an observer would have found if living in a universe. The idea is that, if this value differs too much from what our observers have told us (e.g., in Allen et al., 2008), then it would place a constraint on the extent to which the assumed model can deviate from standard CDM. As in the previous case, we draw random samples of size for the parameters from which we find realisations of . Then, by modelling modified gravity effects using Eq. (24), we compute the mean for the 4 models shown in Fig. 2, and these are shown as curves in Fig. 3 together with the observed values of from Allen et al. (2008). Again, we note that current data favour CDM over all three variants of gravity. As clusters are less screened at late times, we find that low- data is more useful in constraining the model than high- data.

## 5 Discussion and Conclusions

In this paper, we proposed a new cosmological test of gravity, by inferring the cosmic baryon fraction from the apparent gas fractions of massive clusters, and comparing with the results from other, less model-dependent, measurements such as the CMB. In theories with a stronger gravity, the apparent gas fraction is smaller than that in CDM for a fixed . Reversely, if the observed value is fixed, we would find a higher value of than in GR, that can be inconsistent with the model-independent measurements. Taking the Hu-Sawicki model as an example: our quick calculation shows that model parameters are in tension with the gas fraction data of the 42 clusters from Allen et al. (2008), though a more rigorous constraint will be left for future work.

has been a rather widly-used observable (e.g., White et al., 1993), and its power in constraining cosmology – in particular dark energy models – is convincingly demonstrated in various previous works (e.g., Sasaki, 1996; Allen et al., 2004, 2008). The inclusion of baryons opens a new dimension for tests of gravity, since ultimately most cosmological observables can be tracked back to lights emitted by interactions involving baryons. In the mean time, the physics of the X-ray-emitting hot gas in massive clusters is relatively clean, making it easier both for the modelling and to use the observational data. As an example, the assumption that gas temperature depends on the gravitational potential and our main conclusions that (i) – the true gas fraction – is unchanged with modified gravity while (ii) is changed are supported by hydrodynamical simulations in gravity (e.g., He & Li, 2015b). Some uncertainties remain in relating to , but these have been included in the error budget estimate above. Furthermore, within our current state of understanding, slightly changing its modelling (e.g., from Allen et al. (2004) to Allen et al. (2008)) does not change results drastically.

Here, we would like to emphasise the use of effective haloes (He et al., 2015a) in our analysis. Though the idea has a similar origin as that of the dynamical mass of halo (e.g., Schmidt, 2010), there are fundamental differences. Dynamical mass is a certain attribute of a given halo which is defined in the standard way, while effective halo is a completely new way to define and identify haloes. Given an effective halo, all gravitational effect can be calculated from GR, and in particular this means that the way in which is currently extracted from observational data – and the resulting results – can be directly used for our purpose. Thus, with a little extra effort from the people who generate a halo catalogue, the analyses of end users can be made much more straightforward, and this provides an efficient bridge between simulators, theorists and observers.

One may naturally wonder about the generality of this method. As a cosmological test, it relies on galaxy clusters being totally or partially unscreened. Because we are talking about massive clusters which tend to be better screened, this test, like most other cosmological ones, will probably not be able to constrain to substantially smaller than the quoted values here. However, it does provide a fairly clean test – with good observational data available – that has the potential to place one of the strongest constraints from cosmology on gravity. Also, one can always combine and other observables, such as lensing (Terukina et al., 2014), cluster scaling relations (e.g. Arnold, Puchwein & Springel, 2014), and cluster gas pressure profiles (De Martino et al., 2014), to place joint, and likely stronger, constraints. In principle, the test would be more powerful if observational data for smaller galaxy clusters (e.g., those in the mass range ) and galaxy groups are included, because these objects are less screened and so gravity deviates more from GR in general. This, however, requires a better understanding of the feedbacks in different models, which are not well studied so far.

The test can be applied not only to gravity and the more general chameleon theory, but also to similar models such as dilatons (Brax et al., 2011) and symmetrons (Hinterbichler & Khoury, 2010). These models are all featured by a universal coupling of all matter species to a scalar field, that effectively enhances the gravity for all particles (at least in unscreened regimes). There are models in which only certain matter species, e.g.. dark matter, experiences the scalar coupling: therein, baryons can still feel a different gravity depending on how the dark matter particle mass evolves with time, in which case the proposed test does apply. In addition to these theories, we have mentioned above the DGP, Galileon and k-mouflage models. In the first two classes, the deviations from GR are strongly suppressed inside dark matter haloes (Barreira et al., 2013, 2014b), and so we do not expect the new test to work. For K-mouflage, as is for the so-called non-local gravity (Maggiore & Mancarella, 2014; Dirian et al., 2014; Barreira et al., 2014a), there can be a time evolution of Newton’s constant inside clusters, making it possible to use this test. However, one needs to bear in mind that in many of these theories the background evolution history is also modified, and that can affect the test (whether it leads to degeneracies or stronger constraints can only be told by a case-by-case study in the future).

As mentioned earlier, the aim of this paper is to illustrate the main idea of using as a test of gravity theories, and therefore we have made a simplified estimate and have not quoted any numerical results on the confidence levels of the constrained . A more complete and rigorous analysis will require one to relax the simplification that all observed clusters share the same mass, radius and concentration, and use the real observational results of these for all clusters. If the cluster mass is obtained from its dynamical effects, we also need to account for the fact that different clusters may have experienced different degrees of screening, and so a more accurate modelling of the screening is needed to compute the cluster mass profile. These will be left for future work. We note that hydrodynamical simulations for modified gravity theories started to appear recently (e.g., Arnold, Puchwein & Springel, 2014; Hammami et al., 2015; He & Li, 2015b), and such works will be useful for improving the constraining power of this test in the future.

## Acknowledgments

We thank Sownak Bose, Vince Eke, Claudio Llinares and Gongbo Zhao for discussions and comments. The work has used the DiRAC Data Centric system at Durham University, operated by the Institute for Computational Cosmology on behalf of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment was funded by BIS National E-infrastructure capital grant ST/K00042X/1, STFC capital grant ST/H008519/1, and STFC DiRAC Operations grant ST/K003267/1 and Durham University. DiRAC is part of the National E-Infrastructure. BL acknowledges support by the UK STFC Consolidated Grant No. ST/L00075X/1 and No. RF040335. JHH is supported by the Italian Space Agency (ASI) via contract agreement I/023/12/0. LG acknowledges support from NSFC grants Nos. 11133003 and 11425312, the Strategic Priority Research Program The Emergence of Cosmological Structure of the Chinese Academy of Sciences (No. XDB09000000), MPG partner Group family, and an STFC and Newton Advanced Fellowship.

## References

- Allen et al. (2004) Allen S. W., Schmidt R. W., Ebeling H., Fabian A. C., van Speybroeck L., 2004, MNRAS, 353, 457
- Allen et al. (2008) Allen S. W., Rapetti D. A., Schmidt R. W., Ebeling H., Morris R. G., Fabian A. C., 2008, MNRAS, 383, 879
- Arnold, Puchwein & Springel (2014) Arnold C., Puchwein E., Springel V., 2014, MNRAS, 440, 833
- Barreira et al. (2013) Barreira A., Li B., Hellwing W. A., Baugh C. M., Pascoli S., 2013, JCAP, 10, 027
- Barreira et al. (2014a) Barreira A., Li B., Hellwing W. A., Baugh C. M., Pascoli S., 2014, JCAP, 09, 031
- Barreira et al. (2014b) Barreira A., Li B., Hellwing W. A., Lombriser L., Baugh C. M., Pascoli S., 2014b, JCAP, 04, 029
- Brax & Valageas. (2014a) Brax P., Valageas P., 2014a, PRD, 90, 023507
- Brax & Valageas. (2014b) Brax P., Valageas P., 2014b, PRD, 90, 023508
- Brax et al. (2011) Brax P., van de Bruck C., Davis A. C., Shaw D. J., 2010, PRD, 82, 063519
- Bose, Hellwing & Li. (2015) Bose S., Hellwing W. A., Li B., 2015, JCAP, 02, 034
- Carroll et al. (2005) Carroll S. M., de Felice A., Duvvuri V., Easson D. A., Trodden M., Turner M. S., 2005, PRD, 71, 063513
- Cataneo et al. (2014) Cataneo M., Rapetti D., Schmidt F., Mantz A. B., Allen S. W.., Applegate D. E., Kelly P. L., von der Linden A., Morris R. G., 2014, arXiv:1412.0133
- Clifton et al. (2012) Clifton T., Ferreira P. G., Padilla A., Skordis C., 2012, Phys. Rept., 513, 1
- Copeland et al. (2006) Copeland E. J., Sami M., Tsujikawa S., 2006, IJMPD, 15, 1753
- De Martino et al. (2014) De Martino I., De Laurentis M., Atrio-Barandela F., Capozzilello S., 2014, MNRAS, 442, 921
- De Martino et al. (2015) De Martino I., De Laurentis M., Capozzilello S., 2015, in Special Issue Modified Gravity Cosmology: From Inflation to Dark Energy; arXiv:1507.06123
- Deffayet et al. (2009) Deffayet C., Esposito-Farese G., Vikman A., 2009, PRD, 79, 084003
- Dirian et al. (2014) Dirian Y., Foffa S., Khosravi N., Kunz M., Maggiore M., 2014, JCAP, 06, 033
- Dvali, Gabadadze & Porrati (2000) Dvali G., Gabadadze G., Porrati M., 2000, PLB, 485, 208
- Eke, Navarro & Frenk (1998) Eke V., Navarro J. F., Frenk C. S., 1998, ApJ, 503, 569
- Hammami et al. (2015) Hammami A., Llinares C., Mota D. F., Winther, H. A., 2015, MNRAS, 449, 3635
- He et al. (2015a) He J.-h., Hawken A. J., Li B., Guzzo L., 2015a, PRL, in press
- He & Li (2015b) He J.-h., Li B., 2015b, in preparation
- Hinterbichler & Khoury (2010) Hinterbichler K., Khoury J., 2010, PRL, 104, 231301
- Hu & Sawicki (2007) Hu W., Sawicki I., 2007, PRD, 76, 064004
- Joyce et al. (2015) Joyce A., Jain B., Khoury J., Trodden M., 2015, Phys. Rept., 568, 1
- Khoury & Weltman (2004) Khoury J., Weltman A., 2004, PRD, 69, 044026
- Komatsu & Seljak (2001) Komatsu E., Seljak U., 2001, MNRAS, 327, 1353
- Li, Zhao & Koyama (2012) Li B., Zhao G.-B., Koyama K., 2012, MNRAS, 421, 3481
- Lombriser (2014) Lombriser L., 2014, Annalen der Physik, 526, 259
- Maggiore & Mancarella (2014) Maggiore M., Mancarella M., 2014, PRD, 90, 023005
- Mo, van den Bosch & White (2011) Mo H. J., van den Bosch F., White S. D. M., 2011, Galaxy Formation and Evolution, Cambridge University Press
- Navarro, Frenk & White (1997) Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 493
- Nicolis et al. (2009) Nicolis A., Rattazzi R., Trincherini E., 2009, PRD, 79, 064036
- Perlmutter et al. (1999) Perlmutter S. et. al., 1999, ApJ, 517, 565
- Planck Collaboration, Ade et al. (2015) Planck Collaboration. Ade P. A. R., Aghanim N., Arnaud M., et al., 2015, arXiv:1502.01589
- Riess et al. (1998) Riess A. G. et. al., 1998, Astron. J., 116, 1009
- Sasaki (1996) Sasaki, S., 1996, PASJ, 48, L119
- Schmidt (2010) Schmidt F., 2010, PRD, 81, 103002
- Schmidt, Vikhlinin & Hu (2009) Schmidt F., Vikhlinin A., Hu W., 2009, PRD, 80, 083505
- Terukina et al. (2014) Terukina A., Lombriser L., Yamamoto K., Bacon D., Koyama K., Nichol R. C. 2014, JCAP, 04, 013
- White et al. (1993) White S. D. M., Navarro J. F., Evrard A. E., Frenk C. S., 1993, Nature, 366, 429
- Wilcox et al. (2015) Wilcox H., Bacon D., Nichol R. C., et al., 2015, arXiv:1504.03937
- Zhao, Li & Koyama (2011) Zhao G.-B., Li B., Koyama K., 2011, PRL, 107, 071303