# Constraining Modified Gravity Models Using the Thermal Sunyaev- Zeldovich Effect

^{*}

**Corresponding Author:**Ivan De Martino, Department of Theoretical Physics and History of Science, University of the Basque Country UPV/EHU, Faculty of Science and Technology, Barrio Sarriena S/n, 48940 Leioa, Spain, Tel: +34 946 01 20 00, Email: [email protected]

*
Received Date: Sep 01, 2017 /
Accepted Date: Sep 25, 2017 /
Published Date: Sep 29, 2017 *

### Abstract

Some modified gravity theories add, in their weak field limit, a Yukawa-correction to the Newtonian gravitational potential. Such a correction usually depends on the two parameters, one that accounts for the modification of the gravitational constant, and another one representing the scale length on which the scalar field propagates. The thermal Sunyaev-Zeldovich temperature anisotropies can be used to test the modified gravitational potential well demonstrating that the Yukawa-like gravitational potential is able describe the distribution of the hot Intra Cluster Medium without accounting for a Dark Matter halo.

**Keywords:**
Alternative theories of gravity; MOG; Galaxy cluster; Dark matter; Dark energy; Cosmic microwave background; Sunyaev- Zeldovich effect

#### Introduction

Our knowledge about the Universe was increased through the centuries, and observations have quickly improved and became more accurate in the last three decades. We have entered the era of precision cosmology where the observables have been determined within a few percent accuracy [1-6]. The observations indicate that baryons only contribute a few percents to the total amount of the matter and energy in the Universe, and the latter is ongoing a period of accelerated expansion. These data are well accommodated within the framework of the concordance cosmological model (also named ΛCDM model). The model explains the evolution of the Universe from the first fractions of a second to the present day. To account for these observations, the concordance cosmological model that is entirely based on General Relativity (GR), assumes the existence of two extra energy-density components: i) the Dark Matter (DM), and ii) a cosmological constant Λ, equivalent to a perfect fluid with negative pressure, or its generalizations usually named Dark Energy (DE). The DM is characterized by a small temperature, it interacts only gravitationally with the other components, and it constitutes about a 26% of the total amount of energy-density of the Universe. It allow us to explain the emergence of the Large Scale Structure and the dynamics of self-gravitating systems in the framework of GR. The DE accounts for a ∼68% of the total energy-density, and it is required to explain the current period of accelerated expansion [6]. From one hand, the dynamical effects of both DM and DE on large scales are very well constrained; on the other hand, the lack of evidence of counterparts at the particle level can be interpreted as a breakdown of GR at scale beyond the Solar System. Thus, alternative models to GR have been proposed to explain both the dynamics of self-gravitating systems and the cosmological expansion history without resorting to extra components.

Broadly speaking, these alternative models, usually named Extended Theories of Gravity (ETGs), generalize the Hilbert-Einstein Lagrangian by including higher-order curvature invariants and minimally or nonminimally coupled terms between scalar fields and geometry. Such higher-order theories contain extra degrees of freedom that, in the weak field limit, can be recast as new gravitational scale lengths. The general paradigm for (2k+2)-order theories of gravity demonstrates that, a new characteristic scale length arises in Newtonian limit increasing the theory of two derivation orders [7]. Thus, gravity is not longer scale invariant, and DM and DE could be interpreted as the effect of high order theories on scale larger than the Solar system one. This review is focused on two ETGs: (a)the well known f (R)-gravity that replaces the Ricci scalar, R, in the Hilbert-Einstein action with a more general function of the curvature f (R) for comprehensive reviews see [8-12]; and (b) the more recent Scalar-Tensor-Vector Gravity theory (STVG), also known as MOdified Gravity (MOG), that adds scalar, tensor and massive vector fields to the standard Hilbert-Einstein action (for more details see [13,14]. Both ETGs introduce a Yukawa-like correction to the Newtonian gravitational potential in their weak field limit [13-16]. In f (R)-gravity, such correction term is characterized by two parameters the strength δ and the scale length L of the Yukawa-term that are related to the additional degree of freedom/scalar fields arising from the theory. In MOG theory, the mass of the vector field and its strength are governed by two running constants, α and μ, that are promoted to scalar fields and can be constrained by data.

Both theories have been tested from the astrophysical to the cosmological scales. Specifically, f (R)-gravity is able to describe the star formation and evolution [17,18], the emission of gravitational waves for binary systems [19,20], the galactic rotation curves in spiral galaxy [21], the dispersion velocity in elliptical galaxy [22], the dynamics of gas in galaxy clusters [23-25], and the emergence and evolution of the Large Scale Structure [26]. On the other hand, MOG is able to describe the dynamics of self-gravitating system from galactic to extragalactic scales[27-30] and, at the same time, it can explain the evolution of the Universe as a whole [31-33]. Despite their successes, there is not definitive proof about the need of modifying gravity; nowadays, ETGs just represent a valid alternative to the Ë CDM model to overcame some of its shortcomings. Neverthless, having alternatives demand to test them in all possible physical scenarios.

#### Literature Review

Here we briefly review the chance offered by the thermal Sunyaev- Zeldovich (TSZ, [34]) effect in galaxy clusters to constrain both f (R)and MOG theories. We will show that the predicted effect agrees with the observed one when the intra cluster gas is in hydrostatic equilibrium within the modified Newtonian potential. There is no need, in these models, for introducing a dominant DM halo component. This review is organized as follows: in sect. II, we briefly describe the main features of galaxy clusters and the TSZ effect; in sect. III, we describe the methodology used to test ETGs [24,25,30]. Specifically, we illustrate the data, the model and the statistical analysis; in sect. IV, we highlight the results and, finally, in sect. V we give the conclusion.

#### Cluster of Galaxies

Galaxy clusters are the largest virialized objects in the Universe, with a virial mass in the range from 10^{13} to (few) 10^{15} M_{e}. The mass of baryons in cluster is composed at least by two components, diffuse Intra Cluster Medium (ICM), and stars. Neverthless, most of baryons are not in galaxies but they are in the diffuse ICM [35-37].

Clusters contain from hundreds to up to one thousand galaxies within 2 Mpc from their center. They account for a a mere 3% of the total mass of the cluster, while hot ICM gas contains most of the baryons up to 12% of the total mass. The remaining 85% is usually associated to a DM halo. The ICM is highly rarefied: electron number densities are typically n_{e}~ 10^{−4} −10^{−2}cm^{−3} , butit has a temperature in the range from 10^{7} to 10^{8}K, thus clusters are strong X-rays sources with a luminosity L_{X}~ 10^{43} −10^{45}e / rg s . The size of the cluster are typically given in term of the radius r_{500}, defined as the radius at which the mean overdensity of the cluster is 500 times the critical density of the Universe at the same redshift and the corresponding enclosed mass is These definitions are rather convenient since scaling relations based on numerical simulations and X-ray observations allow to determine r_{500} for individual clusters [38,39]. Scaling relations can be also used to test the physics of clusters of galaxies[40-47].

The interest for obtaining a good description of the distribution of the ICM within the gravitational potential well generated by the DM distribution, led to propose phenomenological parameterizations based on the generalized Navarro-Frank-White profile (GNFW) [43,48]. The functional form of the proposed profile is

(1)

where for ΛCDM and

(2)

The parameters were firstly derived using 33 X-Ray selected massive clusters (M_{500}=[10^{14} −10^{15}]) up to redshift Z≤0.2. Their best fit values were the following ones: Recently, combining TSZ and X-Ray data, the Planck Collaboration had fit the pressure profiles out to 3r_{500} for a sample of 62 nearby massive clusters obtaining: [C500,α,β,γ,P0]=(1.81,1.33,4.13,0.31,6.41) [49-51]. Similarly, a study carried on 45 massive clusters with median mass and median redshift Z=0.42 found[C_{500},α,β,γ ,P_{0}]=(1.18,0.86,3.67,0.67,4.29) [52].

TSZ surveys, such as the ones carried out by the Atacama Cosmology Telescope (ACT, [53]) the South Pole Telescope (SPT, [54]) and the Planck satellite [55], will provide new observational dataset, and they will offer new possibilities of testing modified gravity models on Mpc scales by fitting the distribution of the ICM within the modified gravitational potential well.

**Thermal Sunyaev-Zeldovich effect**

Cosmic Microwave Background (CMB) photons cross galaxy clusters, and they are scattered off by the free electrons of the ICM via an inverse Compton scattering. Thus, galaxy clusters distort the CMB blackbody spectrum, and produce secondary temperature anisotropies that are imprinted on the CMB radiation. Two components can distinguished: the thermal component (TSZ, [34]) due to the thermal motion of the electrons in the potential well of the cluster, and the kinematic (KSZ, [56]) component due to the proper motion of the galaxy cluster as a whole with respect to the CMB rest frame. Neglecting relativistic corrections in the electron temperature (Te), the TSZ and KSZ effects are given by

(3)

where, is the direction of the cluster, K_{B} is the Boltzmann constant, m_{e} c^{2} is the electron annihilation temperature, c is the speed of light, is the reduced frequency of observation, and is the peculiar velocity of the cluster. We denoted the cluster optical depth to the SZ effect, with σT Thomson cross section, the electron density evaluated along the line of sight l; T_{0} is the current CMB blackbody temperature andG is the frequency dependence of TSZ effect that, in the non relativistic limit (T_{e} ≈few keV), is given by: For very hot clusters (T_{e} >10 keV), relativistic corrections must be included [57-59].

#### Methodology

The TSZ temperature anisotropies are usually expressed as the integral of the pressure profile along the line of sight:

(4)

To test modified theories of gravity using TSZ effect, the pressure profile (p_{e}) must be specified. It can be computed integrating the following system of differential equations

(5)

(6)

(7)

(8)

where in the case of f(R)-gravity, the effective gravitational potential takes the following functional form [15,16]

(9)

while, for MOG gravity model, it becomes

(10)

In both cases is the classical Newtonian potential. Let us remark that the model does not include any DM component, but it assumes that baryons follow the modified gravitational potential well described in Equations (9) and (10) for f(R) -gravity and MOG models, respectively. When integrating the system of equations (5)-(8), one is assuming that: (i) the gas is in hydrostatic equilibrium within the modified potential well; (ii) the gas distribution is spherically symmetric; and (iii) the state of gas can be described with a polytropic equation of state (equation (6)). Under these assumptions, the system of equations (5)-(8), plus the equation for the modified gravitational potential, constitutes a closed system that can be integrated numerically to compute the pressure profile in ETGs, and to use TSZ anisotropies to constraint the theoretical parameters.

**Foreground cleaned Planck 2013 nominal data: Coma cluster**

Planck 2013 Nominal maps were used to measure the TSZ cluster profile and constrain the parameters of the modified gravitational potential of both theories. The publicly available Planck Nominal maps contain the cosmological CMB signal, instrumental noise, TSZ and KSZ emissions, point and extended infrared sources, thermal dust and CO emissions, and other foreground contributions. Although TSZ has a unique frequency dependence, it is not possible to detect the TSZ anisotropies at cluster location in the Planck Nominal maps. The technique used to reduce such foregrounds emissions provides a cleaned patch of sky centered at the cluster position x, P(v,x) , of size 2° × 2°. Before the cleaning procedure is applied, publicly available Planck CO maps must be subtracted to each frequency channels [60], then the highest frequency channel could be used as a template for the thermal dust emission [61]. Then, the cleaning procedure uses the highest frequency channel as a thermal dust template; at the frequency channel v the cleaned patch is obtained as follows

(11)

where the weight w(v) are given by

(12)

with denoting the area where the weights are computed. Since the instruments operating at the frequencies (30-70 GHz) have angular resolution higher than 10 arcminutes, only High Frequency Instrument (HFI) data at frequencies from 100 to 857 GHz were considered.

As first attempt to build a new test, the study was particularized to the Coma cluster since is located close to the galactic pole where the foreground emission is comparatively low. The TSZ cluster profile was measured at 100, 143, and 353 GHz channels, since the data at 545 GHz were still dominated by the residuals of the dust emission. In **Figure 1**, it is shown a patch around the Coma cluster before (first row) and after (second row) the cleaning procedure was applied to remove foregrounds. The patches clearly show that the Nominal data are dominated by the the intrinsic CMB temperature fluctuations and dust emission, while TSZ signal dominates in the cleaned ones.

Patches of sky centered at the position of A1656 (Coma cluster) at 100-545 GHz. Patches have size of 2° × 2° First row: Planck Nominal maps. Second row: foreground cleaned data.

**MCMC priors**

The model parameters were constrained using a Monte Carlo Markov chain (MCMC) technique. The pipeline employed a Metropolis-Hastings sampling algorithm and used the Gelman-Rubin criteria to test the mixing and convergence of the runs. For each model, four independent chains have been run and the starting points of the parameter space have been randomly chosen. Each run contains at least 40,000 steps with an adapted step size to reach an optimal acceptance rate and to avoid poor mixing and slow convergence [62-66]. The parameter space explored for both f (R)-gravity and MOG models is given in **Table 1**.

Theory | Parameter | Units | Priors |
---|---|---|---|

4*90f(R) | P_{c} |
10^{-2}cm^{-3}keV |
[0.0,3.0] |

γ | - | ||

L | Mpc | [0.1,20.0] | |

δ | - | [-0.99,1.0] | |

*90MOG | P_{c} |
10^{-2}cm^{-3}keV |
[0.0,3.0] |

γ | - | ||

Mpc | [0.01,20.0] | ||

α | - |
[0.1,20.0] |

**Table 1:** MCMC priors of the explored parameter space.

The TSZ data have been measured averaging the emission over disc/rings out to ∼100 arcminutes from center of the Coma cluster. The galaxy cluster was assumed to be spherically symmetric and the ICM to be in hydrostatic equilibrium. The model predictions have been computed at the same apertures to compute the likelihood

(13)

In the previous equation, , where are the data and are the theoretical counterpart. Then, N=22 is the number of data points, p are the parameters of the model, C_{ij} denotes each Planck channel and C_{ij} is the correlation matrix [24,25,30].

#### Results And Discussion

MCMC were run to fit separately the data at each frequency channel. Then, the joint likelihood has been computed. The best fit parameters are summarized in **Table 2**, and the best fit models are shown in **Figure 2** [24,25,30].

**Figure 2:** Best fitting model of f(R)-gravity and MOG in panels left and right, respectively. The black lines represent the best-fit models, the blue lines show the GNFW profile with best fit parameter from [43], and the red line in the right panel illustrates the MOG model with parameter fixed to their “universal” values [75].

Theory | Parameter | Units | Results |
---|---|---|---|

*90f(R) | P_{c} |
10^{-2}cm^{-3}keV |
0.90 ± 0.04 |

γ | - | ||

L | Mpc | 2.19 ± 1.02 | |

δ | - | -0.48 ± 0.22 | |

*90MOG | P_{c} |
10^{-2}cm^{-3}keV |
0.77± 0.03 |

γ | - | ||

Mpc | |||

α | - |

**Table 2:** Best fit parameters for both f(R)-gravity and MOG.

f(R)-gravity MOG

Best fitting model of f(R) -gravity and MOG in panels left and right, respectively. The black lines represent the best-fit models, the blue lines show the GNFW profile with best fit parameter from [43], and the red line in the right panel illustrates the MOG model with parameter fixed to their “universal" values [67].

The analysis gave rise to some interesting results that can be summarized as follows:

• f(R)-theory: the strength of the modified potential in eq. (9) is δ ≠ 0 at the 95% of confidence level (CL). Therefore, the data are compatible with f(R)-gravity plus baryons. Next, the scale length of the potential L is not equal to zero at the 95% CL. This limit corresponds to a Newtonian gravitational potential generated by an effective mass Since the data favors models with would be analogue to the field generated by a cluster containing a large fraction of DM distributed like the baryonic gas [24,25].

• **MOG theory:** the strength of the Yukawa potential in eq. (10), α , is compatible at 68% CL with its universal value α ≅ 8.89 [14,67]. Whereas the universal value of scale length is ruled out at more than3.5σ. Therefore, the assumption that the parameters of the Yukawapotential can be assumed scale independent is also ruled out [30]. In fact, with this assumption, MOG is capable to fit only the central region (θ15 arcminutes) of the galaxy cluster, while it overestimates the TSZ emission at larger apertures: at θ~1degree the departure from the data is almost one order of magnitude (red dashed line in **Figure 2**).

For both theories, f (R) and MOG, the polytropic index is consistent at ∼1.5σ level with the value γ∼1.2 preferred by observations and numerical simulations within the ΛCDM concordance model. Since the physical state of the gas in a galaxy cluster is determined by its formation and evolution [68], the results could be interpreted as an indication that both f (R) and MOG could be able to explain the emergence of the large scale structure as well as the concordance model.

#### Conclusion

Several theories of gravity have been constrained using cluster of galaxies such as chameleon f (R) models [69-71], Galileon model [72], and K-mouflage gravity [73]. Here, we have reviewed the analysis that has been done to constrain those theories that modify the Newtonian potential in their weak filed limit by adding a Yukawa-like term. Specifically, we have reviewed the analysis on analytical f (R)-gravity and MOG models presented in [24,25,30]. Such a particular class of ETGs have been constrained using the TSZ temperature anisotropies due to the Coma cluster. With respect to other works, the model assumes that: (a) the hot ICM is in hydrostatic equilibrium within the modified gravitational potential wells described by the Equations (9) and (10) for f (R)-gravity and MOG, respectively; (b) the physical state of the gas is well described by a polytropic equation of state. These models have been tested using a foreground cleaned version of the 2013 data release of Planck Nominal maps. The measured TSZ profile at the location of Coma cluster has been used to constrain the model parameters employing a MCMC algorithm. Both analyses, for f (R) and MOG, show the capability of the two theories to describe galaxy cluster without accounting for a DM component. However, the analyses are based on the assumptions of hydrostatic equilibrium and spherical symmetry of the cluster. Although such assumptions hold in the intermediate region of the Coma cluster where models are tested, they do not hold in general due to the presence of substructures, turbulences, and physical processes that can heat up or cool down the gas in the cluster core [74-81]. The departure from the spherical symmetry could affect both the innermost and outermost regions of Coma cluster. This is figured out in the the degeneracy between the strength of the gravitational potential,δ and α , and the polytropic index of the gas [25,30]. Neverthless, a deeper study of this degeneracy should use N-body hydrodynamical simulations carried out for each set of parameter. Another limitation of the analysis was the angular resolution of the foreground cleaned data (FWHM =10 arcminutes). This could be overcame with the next-generation of full sky CMB missions such as COrE/PRISM [82] that will have a much higher angular resolution and frequency coverage providing a powerful tool to properly investigate the relation between the underlying theory of gravity and the baryonic processes.

#### References

- Riess AG, Strolger LG, Tonry J, Casertano S, Ferguson HC, et al.(2004) Type Ia supernova discoveries at
- Astier P, Guy J, Regnault N, Pain R, Aubourg E, et al. (2006) The supernova legacy survey: measurement of Ω
_{M}, Ω_{Λ}and*w*from the first year data set. Astron Astrophys 447: 31. - Clocchiatti A, Schmidt BP, Filippenko AV, Challis P, Coil AL, et al. (2006) Hubble space telescope and ground-based observations of type ia supernovae at redshift 0.5: Cosmological implications. Astrophys J 642: 1.
- Blake C, Kazin EA, Beutler F, Davis TM, Parkinson D et al. (2011) The wigglez dark energy survey: mapping the distance-redshift relation with baryon acoustic oscillations. Mon Not R Astron Soc 418: 1707.
- Hinshaw G, Larson D, Komatsu E, Spergel DN, Bennett CL, et al. (2013) Nine-year Wilkinson Microwave Anisotropy Probe (WMAP) observations: Cosmological parameter results. Astrophys J Suppl 208: 19H.
- Ade PAR, Aghanim N, Arnaud M, Ashdown M, Aumont J, et al. Planck collaboration (2015) Planck 2015 results. XIII. Cosmological parameters. Astron Astrophys 594: A13
- Quandt I, Schmidt HJ (1991) The Newtonian limit of fourth and higher order gravity. Astron Nachr 312: 97.
- Carroll SM, Duvvuri V, Trodden M, Turner MS (2004) Is cosmic speed-up due to new gravitational physics? Phys Rev D 70: 043528.
- Sotiriou TP, Faroni V (2010) f(R) theories of gravity. Rev Mod Phys 82: 451.
- Nojiri S, Odintsov SD (2006) Modified f(R) gravity consistent with realistic cosmology: From a matter dominated epoch to a dark energy universe. Phys Rev D 74: 086005.
- Nojiri S, Odintsov SD (2011) Unified cosmic history in modified gravity: From F(R) theory to Lorentz non-invariant models. Phys Rept 505: 59.
- Capozziello S, De Laurentis M (2011) Extended theories of gravity. Phys Rept 509: 167.
- Moffat JW (2006) Scalar tensor vector gravity theory. J Cosmol Astropart Phys 0603: 004.
- Moffat JW, Rahvar S (2013) The MOG weak field approximation and observational test of galaxy rotation curves. Mon Not R Astron Soc 436: 1439-1451.
- Capozziello S, Stabile A, Troisi A (2012) Fourth-order gravity and experimental constraints on eddington parameters. Modern Physics Lett A 24: 09.
- Capozziello S, De Laurentis M (2012) The dark matter problem from f(R) gravity viewpoint. Ann Phys 524:1.
- Capozziello S, De Laurentis M, Odintsov SD, Stabile A (2011) Hydrostatic equilibrium and stellar structure in f(R) gravity. Phys Rev D 83: 064004.
- Capozziello S, De Laurentis M, De Martino I, Formisano M, Odintsov SD (2012) Jeans analysis of self-gravitating systems in f(R) gravity. Phys Rev D 85: 044022.
- De Laurentis M, De Martino I (2013) Testing f (R) theories using the first time derivative of the orbital period of the binary pulsars. Mon Not R Astron Soc 43: 741.
- De Laurentis M, De Martino I (2015) Probing the physical and mathematical structure of f(R)-gravity by PSR J0348 + 0432. Int J Geom Methods Mod Phys 12: 1550040.
- Cardone VF, Capozziello S (2011) Systematic biases on galaxy haloes parameters from Yukawa-like gravitational potentials. Mon Not R Astron Soc 414: 1301.
- Napolitano NR, Capozziello S, Romanowsky AJ, Capaccioli M, Tortora C (2012) Testing Yukawa-like potentials from f(R)-gravity in elliptical galaxies. Astrophys J 748: 87.
- Capozziello S, De Filippis E, Salzano V (2009) Modelling clusters of galaxies by f(R) gravity. Mon Not R Astron Soc 394: 947.
- De Martino I, De Laurentis M, Atrio-Barandela F, Capozziello S (2014) Constraining f(R) gravity with Planck data on galaxy cluster profiles. Mon Not R Astron Soc 442: 921.
- De Martino I (2016) f(R)-gravity model of the Sunyaev-Zeldovich profile of the Coma cluster compatible with Planck data. Phys Rev D 93: 4043.
- De Martino I, De Laurentis M, Capozziello S (2015) Constraining f(R)gravity by the large scale structure. Universe 1: 123.
- Brownstein JR, Moffat JW (2006) Galaxy cluster masses without non-baryonic dark matter. Mon Not R Astron Soc 367: 527.
- Brownstein JR, Moffat JW (2006) Galaxy rotation curves without nonbaryonic dark matter. Astrophys J 636: 721.
- Brownstein JR, Moffat JW (2007) The bullet cluster 1E0657-558 evidence shows modified gravity in the absence of dark matter. Mon Not R Astron Soc 382: 29.
- De Martino I, De Laurentis M (2017) On the universality of MOG weak field approximation at galaxy cluster scale. Phys Lett B 770: 440.
- Moffat JW (2006) Late-time inhomogeneity and acceleration without dark energy. J Cosmol Astropart Phys 05: 001.
- Moffat JW (2007) A modified gravity and its consequences for the solar system, astrophysics and cosmology. Int J Mod Phys 16: 2075.
- Moffat JW, Toth VT (2009) The bending of light and lensing in modified gravity. Mon Not R Astron Soc 397: 1885.
- Sunyaev RA, Zeldovich YB (1972) The observations of relic radiation as a test of the nature of X-ray radiation from the clusters of galaxies. Comments on Astrophys Space Phys 4: 173.
- White SDM, Navarro JF, Evrard AE, Frenk CS (1993) Nature 366: 429-433.
- White DA, Fabian AC (1995) Einstein observatory evidence for the widespread baryon overdensity in clusters of galaxies. Mon Not R Astron Soc 273: 72-84.
- Lubin LM, Cen R, Bahcall NA, Ostriker JP (1996) The baryon fraction and velocity--temperature relation in galaxy clusters: Models versus observations Astrophys J 460: 10.
- Kaiser N (1986) Evolution and clustering of rich clusters. Mon Not R Astron Soc 222: 323
- Kravtsov AV, Borgani S (2012) Formation of galaxy clusters. Ann Rev AstronAstrophy 50: 353-409.
- Böhringer H, Schuecker P, Pratt G, Arnaud M, Ponman T, et al. (2007) The representative XMM-Newton cluster structure survey (REXCESS) of an X-ray luminosity selected galaxy cluster sample. A&A 469: 363-377.
- Bonamente M, Joy M, La Roque SJ, Carlstrom JE, Nagai D, et al. (2008) Scaling relations from Sunyaev-Zel’dovich effect and chandra X-Ray measurements of high-redshift galaxy clusters. Astrophys J 675: 106
- Marrone DP, Smith GP, Richard J, Marshall J, Bonamente M, et al. (2009) LoCuSS: A comparison of Sunyaev-Zel’dovich effect and gravitational-lensing measurements of galaxy clusters. Astrophys J 701: L114
- Arnaud M, Pratt GW, Piffaretti R, Böhringer H, Croston JH, et al. (2010) The universal galaxy cluster pressure profile from a representative sample of nearby systems (REXCESS) and the
- Melin J, Bartlett JG, Delabrouille J, Arnaud M, Piffaretti R, et al. (2011) The galaxy cluster YSZ - LX and YSZ-M relations from the WMAP 5-yr data Astron Astrophys 525: A139.
- Andersson K, Benson BA, Ade PAR, Aird KA, Armstrong B et al. (2011) X-ray properties of the first Sunyaev-Zel’dovich effect selected galaxy cluster sample from the south pole telescope. Astrophys J 738: 48
- De Martino I, Atrio-Barandela F (2016) SZ/X-ray scaling relations using X-ray data and Planck nominal maps. Mon Not R Astron Soc 461 3: 3222-3232
- Jimeno P, Diego JM, Broadhurst T, De Martino I, Lazkoz R (2017) Planck/SDSS cluster mass and gas scaling relations for a volume-complete redMaPPer sample. arXiv:1706.00395.
- Navarro JF, Frenk CS, White SDM (1997) A universal density profile from hierarchical clustering. Astrophys J 490: 493
- Planck Collaboration (2011a) Planck early results X: Statistical analysis of Sunyaev-Zeldovich scaling relations for X-ray galaxy clusters. A&A 536 A10.
- Ade PAR, Aghanim N, Arnaud M, Ashdown M, Aumont J, et al. (2011) Planck early results. XI. Calibration of the local galaxy cluster Sunyaev-Zeldovich scaling relations. Astron Astrophys 536: A11
- Ade PAR, Aghanim N, Arnaud M, Ashdown M, Atrio-Barandela F, et al. (2013) Planck intermediate results V: Pressure profiles of galaxy clusters from the Planck survey. Astron Astrophys 550: 131
- Sayers J, Czakon NG, Mantz A, Golwala SR, Ameglio S, et al. (2013). Sunyaev-Zel’dovich-measured pressure profiles from the Bolocam X-Ray/SZ galaxy cluster sample. Astrophys J 768: 177
- Kosowsky A (2003) The atacama cosmology telescope. New Astron Rev 47: 939
- Carlstrom JE, Ade PAR, Aird KA, Benson BA, Bleem LE, et al. (2011) The 10 meter south pole telescope. Publ Astron Soc Pac 123: 568
- Ade PAR, Aghanim N, Alves MIR, Armitage-Caplan C, Arnaud M, et al. (2014) Planck 2013 results. Overview of products and scientific results. Astron Astrophys 571: A1
- Sunyaev R, Zeldovich Y (1980) The velocity of clusters of galaxies relative to the microwave background - The possibility of its measurement. Mon Not R Astron Soc 190:413.
- Itoh N, Kohyama Y, Nozawa S (1998) Relativistic corrections to the Sunyaev-Zeldovich effect for clusters of galaxies. Astrophys J 5027
- Nozawa S, Itoh N, Kohyama Y (1998) Relativistic corrections to the Sunyaev-Zeldovich effect for clusters of galaxies. II. Inclusion of peculiar velocities. Astrophys J 508: 17
- Nozawa S, Itoh N, Suda Y, Ohhata Y (2006) An improved formula for the relativistic corrections to the kinematical Sunyaev-Zeldovich effect for clusters of galaxies. Nuovo Cimento 121: 487
- Ade PAR, Aghanim N, Alves MIR, Armitage-Caplan C, Arnaud M, et al. (2014) Planck 2013 results. XIII. Galactic CO emission. Astron Astrophys 571: 13
- De Martino I, Génova-Santos R, Atrio-Barandela F, Ebeling H, Kashlinsky A, et al. (2015) Constraining the redshift evolution of the cosmic microwave background blackbody temperature with PLANCK data. Astrophys J 808: 128
- Hastings WK (1970) Monte Carlo sampling methods using Markov chains and their applications. Biometrika 57: 97-109.
- Metropolis N, Rosenbluth AW, Rosenbluth MN, Teller AH (1953) Equation of state calculations by fast computing machines. J Chem Phys 21: 1087
- Gelman A, Rubin DB (1992) Inference from iterative simulation using multiple sequences. Statist Sci 7: 457
- Gelman A, Roberts GO, Gilks WR (1996) Efficient metropolis jumping rules. Bayesian statistics 5: 599-608
- Roberts GO, Gelman A, Gilks WR (1997) Weak convergence and optimal scaling of random walk Metropolis algorithms. Ann Appl Probab 7 1: 110
- Moffat JW, Rahvar S (2014) The MOG weak field approximation - II. Observational test of Chandra X-ray clusters. Mon Not R Astron Soc 441: 3724
- Bertschinger E (1985) Self-similar secondary infall and accretion in an Einstein-de Sitter universe. Astrophys J Suppl 58: 39
- Terukina A, Lombriser L, Yamamoto K, Bacon D, Koyama K, et al. (2014) Testing chameleon gravity with the Coma cluster. J Cosmol Astropart Phys 04: 013
- Wilcox H, Bacon D, Nichol RC, Rooney PJ, Terukina A, et al. (2015) The XMM cluster survey: testing chameleon gravity using the profiles of clusters. Mon Not R Astron Soc 452: 1171
- Li B, He JH, Gao L (2015) Cluster gas fraction as a test of gravity. Mon Not R Astron Soc 456: 146
- Terukina A, Yamamoto K, Okabe N, Matsushita K, Sasaki T (2015) Testing a generalized cubic Galileon gravity model with the Coma Cluster. J Cosmol Astropart Phys 10: 64
- Brax P, Rizzo LA, Valageas P (2015) K-mouflage effects on clusters of galaxies. Phys Rev D 92: 043519
- Arlen T, Aune T, Beilicke M, Benbow W, Bouvier A, et al. (2012) Constraints on cosmic rays, magnetic fields, and dark matter from gamma-ray observations of the coma cluster of galaxies with VERITAS and fermi. Astrophys J 757: 123
- Churazov E, Vikhlinin A, Zhuravleva I, Schekochihin A, Parrish I, et al. (2012) X-ray surface brightness and gas density fluctuations in the Coma cluster. Mon Not R Astron Soc 421: 1123
- Simionescu A, Werner N, Urban O, Allen SW, Fabian AC, et al. (2013) Thermodynamics of the Coma Cluster Outskirts. Astrophys J 775, 4
- Sato T, Matsushita K, Ota N, Sato K, Nakazawa K, et al. (2011) Suzaku observations of Iron K-lines from the intracluster medium of the coma cluster. Publ Astron Soc Jpn 63: 991
- Neumann DM, Lumb DH, Pratt GW, Briel UG (2003) The dynamical state of the coma cluster with XMM-Newton 400: 811-821
- Sanders JS, Fabian AC, Churazov E, Schekochihin AA, Simionescu A, et al. (2013) Linear structures in the core of the coma cluster of galaxies. Science 341: 6152
- Snowden SL, Mushotzky RF, Kuntz KD, Davis DS (2008) A catalog of galaxy clusters observed by XMM-Newton. Astron Astrophys 478: 615
- Wik DR, Sarazin CL, Finoguenov A, Matsushita K, Nakazawa K, et al. (2009) A Suzaku search for Nonthermal emission at hard X-Ray energies in the coma cluster. Astrophys J 696: 1700
- André P, Baccigalupi C, Banday A, Barbosa D, Barreiro B, et al. (2014) PRISM (Polarized Radiation Imaging and Spectroscopy Mission): an extended white paper. J Cosmol Astropart Phys 1402: 006

Citation: De Martino I (2017) Constraining Modified Gravity Models Using the Thermal Sunyaev-Zeldovich Effect. J Astrophys Aerospace Technol 5: 150. Doi: 10.4172/2329-6542.1000150

Copyright: © 2017 De Martino I. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.