Issue 
A&A
Volume 530, June 2011



Article Number  A68  
Number of page(s)  6  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201016273  
Published online  11 May 2011 
COSMOS weaklensing constraints on modified gravity^{⋆}
^{1}
Centro de Astronomia e Astrofísica da Universidade de Lisboa, Tapada da Ajuda, 1349018 Lisboa, Portugal
email: tereno@fc.ul.pt
^{2}
Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, The Netherlands
Received: 6 December 2010
Accepted: 15 March 2011
The observed acceleration of the universe, explained through dark energy, could alternatively be explained through a modification of gravity that would also induce modifications in the evolution of cosmological perturbations.
We use new weak lensing data from the COSMOS survey to test for deviations from general relativity (GR). The departure from GR is parametrized in a modelindependent way that consistently parametrizes the twopoint cosmic shear amplitude and growth.
Using CMB priors, we perform a likelihood analysis. We find constraints on the amplitude of the signal that do not indicate a deviation from GR.
Key words: gravitational lensing: weak / largescale structure of Universe / cosmological parameters
Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data archives at the Space Telescope European Coordination Facility and the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 526555.
© ESO, 2011
1. Introduction
The ΛCDM model is currently the best fit to the available cosmological data. In this framework, the period of accelerated expansion of the Universe is due to dark energy, a dominant component of the energy density with a negative equation of state. A possible alternative, worth to investigate observationally, is to replace the standard gravitational theory of general relativity (GR) by a modification that admits selfaccelerating solutions. Such gravity theory must converge to GR both at high redshift, when the expansion is not accelerated, and on small scales, since GR has consistently passed local tests of gravity (Will 2006).
The various classes of modified gravity and dark energy models (Uzan 2007) cannot be distinguished based only on observations of the cosmological background. They are only distinguishable due to their different predictions for cosmological structure formation (Zhang et al. 2007). The evolution of cosmological background and perturbations have been studied in various classes of physically motivated modified gravity theories (Nojiri & Odintsov 2010). Most prominent today are theories where the modifications arise from extra dimensions, with matter fields confined to a 4dimensional brane embedded in a higher dimensional bulk, such as the DvaliGabadadzePorrati (DGP) model (Dvali et al. 2000; Deffayet et al. 2002; Lue et al. 2004), and also the socalled f(R) theories (Nojiri & Odintsov 2003; Carroll et al. 2004; Song et al. 2007). In the latter, the field equations are modified because the corresponding Einstein action no longer depends linearly on the Ricci scalar (R) but there are extra terms on R.
The standard way of studying cosmological structure formation is to introduce scalar perturbations in the FriedmannRobertsonWalker (FRW) spacetime metric. In the Newtonian or longitudinal gauge (Mukhanov et al. 1992), where by construction the metric remains diagonal, there are only two scalar perturbations: the scalar fields Ψ (the Newtonian potential) and Φ (the spatial curvature potential). They are sourced by the cosmological fluid perturbations contained in the stressenergy tensor. Conservation of the stressenergy tensor imposes evolution equations for its components. In addition, the equations of the theory of gravitation, relating the metric components with the energy components, bring extra evolution and constraint equations.
In modified gravity, structure formation has been studied with various parametric and nonparametric approaches (Bertschinger & Zukin 2008; Silvestri & Trodden 2009). Parametric approaches include direct parametrization of the linear growth (Linder 2005; Linder & Cahn 2007); the socalled parametrized postFriedmann framework (PPF) (Hu & Sawicki 2007) and several related phenomenological postGR parametrizations, which give a modelindependent description of deviations to GR in a way that is convenient for the testing of cosmological predictions (Caldwell et al. 2007; Amendola et al. 2008; Daniel et al. 2008, 2010; Pogosian et al. 2010); and phenomenological applications to parametrizations of particular theories of modified gravity (Acquaviva et al. 2008; Zhao et al. 2009; Song et al. 2010).
The various metric and energy perturbation variables in modified gravity scenarios may be measured by different probes (Huterer & Linder 2007; Jain & Zhang 2008; Schmidt 2008; Song & Doré 2009; Guzik et al. 2010; Beynon et al. 2010). In particular, weak lensing by largescale structure (cosmic shear), which probes the sum Φ + Ψ, is well suited for such studies (Uzan & Bernardeau 2001; Knox et al. 2006). In cosmic shear measurements, the gravitational information is contained in the equations relating the measured matter density perturbations with the two potentials, or equivalently in one of those equations (a generalized Poisson equation) and an equation relating the two potentials (the anisotropy equation). A parametrization of those two relations is thus well suited for cosmic shear measurements, which will provide a degenerate estimation of those two quantities. This is the approach we follow in this work.
We compute cosmic shear twopoint correlation models in tomographic bins, using standard cosmological parameters plus two modified gravity parameters, and analyse them against the new COSMOS tomographic cosmic shear measurements of Schrabback et al. (2010). Other cosmic shear data have already been used in recent studies of modified gravity, using similar or different parametrizations (Doré et al. 2007; Thomas et al. 2009; Daniel et al. 2009; Bean & Tangmatitham 2010; Zhao et al. 2010; Daniel & Linder 2010). The standard parameters are constrained by using WMAP7year CMB data (Komatsu et al. 2011). Our goal is to study the impact of this phenomenological model in the degeneracy between the amplitude of the cosmic shear signal and the linear growth of structure. In other words, to investigate if, given the strong cosmological priors imposed by the WMAP7 results, there is room to compensate a possible alteration of the predicted redshiftdependent amplitude of the cosmic shear effect with a modification of the growth of matter perturbations, which, phenomenologically, would keep open the possibility of an alternate gravity theory as an origin for an effective dark energy behaviour.
In the next section, we present the parametrization of modified gravity used. The new COSMOS dataset is described in Sect. 3. The statistical analysis and its results are presented and discussed in Sects. 4 and 5. We conclude in Sect. 6 with a reminder of the main assumptions made.
2. Parametrization
Images of background galaxies are distorted (by convergence and shear) in the perturbed FRW spacetime. In the Newtonian gauge, considering scalar perturbations and a flat background, this spacetime is described by (1)The local deflection of null geodesics, with respect to unperturbed ones, depends on the traveling time of a light ray, and thus on the sum of the metric perturbations, (2)Integrating the multiple deflections over cosmological distances, and introducing the lens equation, the cosmic convergence is written, to first order in the amplification matrix and for a flat Universe, as (Bartelmann & Schneider 2001) (3)where (θ,w) are comoving spatial coordinates and the Laplacian of the lenses potentials at comoving distances w′ is integrated up to the source at distance w.
The equations of the theory of gravitation at perturbation level relate the metric perturbations (the potentials Ψ, Φ) to the perturbations of the cosmological fluid (matter density contrast δ, pressure perturbation δp, fluid velocity v, and stress or anisotropic pressure σ). In GR these equations are the linearized Einstein equations. The phenomenological approach parametrizes deviations from GR by introducing extra functions in the linearized Einstein equations.
Following Ma & Bertschinger (1995), a combination of the 0 − 0 and 0 − i linearized Einstein equations relates Φ with the comoving matter density perturbation: it is the generalized Poisson equation. In Fourier space, with k^{2} being the spatial Laplacian, it writes, (4)where Δ is the comoving density perturbation Δ = δ + 3H(1 + w)v/k^{2}, w is the equationofstate parameter w = p/ρ, and the second equality is written for dark matter inserting Ω_{m}. The deviation from GR is parametrized by introducing the mass screening phenomenological postGR function Q(a,k), replacing the gravitational constant by an effective function, (5)Notice that a hypothetical clustering dark energy component would also add a contribution to the Poisson equation. Factorizing out the dark matter density contrast, a multiplicative term, function of δ_{DE}/δ, would appear, and play the same role as Q.
The i − j Einstein equation gives a second constraint, relating the difference Φ − Ψ with the stress σ: it is the anisotropy equation, (6)In GR the two potentials are identical. The deviation from GR is parametrized by introducing the gravitational slip phenomenological postGR function η(a,k), (7)Again, a hypothetical dark energy component with anisotropic stress σ would have the same effect as Qη (Koivisto & Mota 2006).
The two remaining independent Einstein equations, for example 0 − i and i − i, are evolution equations for Ψ and Φ, and will not be needed when considering the quasistatic approximation, valid for subhorizon scales.
The linearized Einstein equations allow thus for 2 independent postGR functions. From Eqs. (5) and (7), we can also write modified equations for the Laplacians of Φ + Ψ and Ψ. They are (8)and (9)These equations introduce two derived postGR functions that are useful for cosmic shear studies: Σ(a) = Q(a)(1 + η(a)/2) and Γ(a) = Q(a)(1 + η(a)).
One approach, taken e.g. by Daniel & Linder (2010) and Zhao et al. (2010), to study the postGR modifications as function of scale and redshift is to parametrize them with various constant parameters in redshift and scale bins. Here, we consider them to be scaleindependent in the limited range of scales provided by the COSMOS data, although it is clear that in a broad range of scales there is scale dependence at least to differentiate the GR behaviour on solar system scales from the possible modifications on cosmological scales.
Regarding the redshift dependence, we adopt a functional form using the scaling argument of Caldwell et al. (2007), according to which the departure from GR originates an effective smooth dark energy Λ with constant w = −1. The deviations in the growth history scale with the deviations in the expansion history, such that they do not turn on too early, deviating from GR at early times, or too late. In this approach, the background evolution is identical to ΛCDM, no further modification is assumed in the Friedmann equations, and the postGR functions take the GR values at a ≪ 1 (Q = 1, η = 0 or Σ = 1, Γ = 1) and deviate from GR proportionally to ρ_{DE}(a)/ρ_{m}(a). This implies that (10)This dependence for the fundamental functions Q(a) and η(a) implies a different redshift dependence for the derived functions Σ(a) and Γ(a). Only in the case of small parameter values, (Q_{a} < 1, η_{a} < 1), do the derived functions have a similar scaling, (11)where Σ_{a} = Q_{a} + η_{a}/2 and Γ_{a} = Q_{a} + η_{a}. We will work in this ansatz, which reduces the range of values allowed for the postGR parameters. In particular, we will consider models such that the difference between the two potentials is less than 50% and the effective gravitational constant does not differ more than 50% from the GR value. Large deviations are better studied with direct parametrizations in redshift bins, in the absence of a wellmotivated functional form. Also, the approach of modifying the linearized Einstein equations should break down for too large deviations.
Going back to the convergence, Eq. (3), we compute its power spectrum (Bartelmann & Schneider 2001) and, inserting Eqs. (5) and (7), obtain the cosmic shear power spectrum in the following form, (12)The amplitude of the power spectrum depends on the derived function Σ(a). This equation considers the correlation of source galaxies from two redshift bins, i,j, with redshift distributions p (w), which are lensed with geometrical efficiency factors given by, (13)Deviations to GR thus modify the cosmic shear amplitude via Σ. They further impact the cosmic shear power spectrum via a modification of the matter power spectrum, P_{δ}, present in Eq. (12). The conservation of stressenergy provides the evolution equations for the matter perturbations. The evolution equations for δ and v depend on the stressenergy components, the time evolution of the potentials and on the Laplacian of the Newtonian potential (Ma & Bertschinger 1995). Inserting the constraint equations Eqs. (5) and (7), and writing for the comoving density Δ, a great number of terms arise. The solution is a scaledependent growth in the linear regime, even for scaleindependent postGR parameters (Daniel et al. 2010).
For the case of subhorizon scales, assuming time derivatives of the perturbations are negligible compared to the spatial derivatives, and Δ ≈ δ, the evolution equation of the cold dark matter perturbation greatly simplifies. It only contains a Hubble drag term and a term on k^{2}Ψ which means, from Eq. (9), that the growth depends directly on the function Γ(a).
In this regime, the growth solution can be cast in a form that explicitly separates the expansion history and gravitational dependencies of the growth (Linder 2005), (14)This introduces the growth parameter, which contains the gravitational dependence, and has the value γ = 6/11 in GR (Wang & Steinhardt). In the case of ΛCDM background evolution and scaling evolution of Γ(a), it takes the value (Linder & Cahn 2007; Daniel et al. 2010), (15)
3. COSMOS cosmic shear data
The new COSMOS cosmic shear catalogue (Schrabback et al. 2010, hereafter S10) contains around 450 000 galaxies with magnitude i < 26.7 in an area of 1.64 deg^{2}, corresponding to a density of 76 galaxies/arcmin^{2}.
The cosmic shear analysis of S10 used the COSMOS30 photometric redshift catalogue (Ilbert et al. 2009), which provides individual photometric redshifts for galaxies of magnitude i < 25. Those were grouped in 5 redshift bins. The highest magnitude galaxies (i > 24) of the lowest redshift bin (z < 0.6) suffer from catastrophic errors due to the dificulty to locate the 4000 Å break in their spectra. They are thus likely to contain highredshift galaxies and were removed from the analysis. In addition, the redshift distribution of galaxies with i > 25, for which individual redshifts were not available, was estimated in S10 and those galaxies were placed in a single wide redshift bin. Table 1 shows the 6 redshift bins used, containing the remaining roughly 420 000 galaxies.
Description of the redshift bins used in the analysis, with their amount of contributing galaxies and mean redshifts.
The data used in this work consists on a vector of shear correlations ξ_{+} and ξ_{ − } computed with the estimator given in Eq. (13) of S10. We use all the crosscorrelations between redshift bins (12, 13, 14, 15, 23, 24, 25, 34, 35, 45, 16, 26, 36, 46, 56) and the autocorrelation of the wide sixth bin. The autocorrelations of the narrow bins (11, 22, 33, 44, 55) were not used to minimize the effect of intrinsic alignments. In all cases, the shear correlations were computed on 5 angular separation bins ranging from 0.7 to 21 arcmin. The data vector has thus (15 + 1) × 5 × 2 = 160 components. S10 found the Bmodes to be consistent with zero at all scales. In this case, we work directly with the ξ_{ ± } correlations and there is no need to make an E/B mode decomposition of the shear field, which has the disadvantage of being nonlocal.
The 160 × 160 covariance matrix of the shear correlations was estimated in S10 from 288 realisations of COSMOSlike fields obtained from raytracing through the Millenium simulation (Hilbert et al. 2009). The inverse covariance matrix, needed for the likelihood analysis, was debiased by applying the correction of Hartlap et al. (2007).
4. Analysis
We build models on a threedimensional grid consisting of (Q_{a},η_{a}) and a redshiftcalibration parameter f_{z}.
The parameter f_{z} is a nuisance parameter that allows to include an uncertainty on the redshift calibration of bin 6. Following S10, we marginalize over this uncertainty by including shifted redshift distributions on the grid of models, such that each model is evaluated assuming various redshift distributions p (f_{z} z) in the range 0.9 < f_{z} < 1.1.
We choose to work with the two original phenomenological parameters, to be able to have an uniform coverage of the (Q_{a},η_{a}) plane. We probe a relatively narrow range of values: −0.5 < Q_{a} < 0.5 and −0.5 < η_{a} < 0.5. As discussed in Sect. 2, this choice allows the redshiftdependence of the postGR derived parameters to scale with the background expansion. The grid mapped on the plane of the derived parameters, the ones that have a direct impact on the cosmic shear signal, i.e., the amplitude Σ_{a} and the growth Γ_{a} or γ, introduces a correlation between these 2 parameters. Its range is such that Σ(a = 1) can deviate from the GR value, Σ(a = 1) = 1, by up to 70%, while γ extends to a maximum deviation of 40% of the GR value of γ = 0.55. (Hereafter we denote the values of the postGR parameters at a = 1 by Q, η, Σ and Γ.)
The grid is computed for each point of a sample of standard cosmological models built by thinningout the ΛCDM WMAP7 year Monte Carlo Markov Chains available in the WMAP archive^{1}. In this way, the analysis has strong priors on the standard cosmological parameters. The average σ_{8} and Ω_{m} values of the models used are σ_{8} = 0.807 and Ω_{m} = 0.271.
Fig. 1 Correlation function . Top panel shows the measurement (points with error bars), the bestfit GR model for a nontomographic analysis (Ω_{m} = 0.27, σ_{8} = 0.76) (solid), modified gravity models with the same standard parameters but (Σ = 1.6, γ = 0.55) (dashed), or (Σ = 1, γ = 0.65) (dotted). Bottom panel shows the deviation between the same dashed and dotted models with respect to the GR model. 

Open with DEXTER 
For each model, we evaluate the various shear power spectra using Eq. (12). The matter power spectrum is computed using the transfer function of Eisenstein & Hu (1998) and the linear growth of Eq. (14) with γ given in Eq. (15). Since the small scales of the COSMOS data are in the quasilinear and nonlinear regime, we need a nonlinear prescription to correct the linear power spectrum. We compute the nonlinear power spectrum using the formula of Smith et al. (2003), which was built from Nbody simulations of ΛCDM. Finally, the shear correlations are computed from (16)where J_{n} is the nth order Bessel function of the first kind.
As an example, Fig. 1 shows the ξ_{ + } autocorrelation function for the wide redshift bin. In the models shown, Σ and γ take the extreme values used in the grid, when the other parameter is kept fixed at GR. Consistently with the approximations described before, the effect is scaleindependent. We notice the impact of γ is very small on the COSMOS angular scales whereas the impact of Σ clearly dominates. This illustrates why cosmic shear analyses were not yet able to constrain γ, while tighter constraints were obtained using the mass function of galaxy clusters (Rapetti et al. 2010). For cosmic shear studies, the effect of the growth parameter would be more important on larger scales, where the linear growth is scaledependent due to the k^{2} factor in the Poisson equation (see also Daniel et al. 2010).
The likelihood of each model of the grid is evaluated with respect to the previously described COSMOS data vector of dimension 160. The likelihood values are then used to update the original weights of the chain, in order to obtain the joint CMBcosmic shear weights according to the method of importance sampling (Lewis & Bridle 2002).
The resulting marginalised Probability Distribution Functions are shown in Fig. 2, where it is clear that the phenomenological parameters Q and η are not constrained by the data in the range of values probed.
Fig. 2 Marginalised onedimensional PDFs obtained for the phenomenological postGR parameters. 

Open with DEXTER 
Fig. 3 The dotted lines show the limits of the probed region. The 68% C.L. likelihood contour essentially fills up that region, showing that for these data the range of models chosen constitutes a significant prior. 

Open with DEXTER 
5. Results and discussion
We consider now the derived parameters, Σ and γ, the ones that directly affect the cosmic shear signal. Figure 3 shows the result in the (Σ,γ) plane. In this analysis, given the scaleindependent approximations made (both in the choice of scaleindependent parameters and in the use of small scales where the modifications to growth are scaleindependent), the postGR models differ only in their amplitudes and redshiftevolution of the amplitudes. Σ and γ are thus degenerate and we constrain the following combination of parameters: (17)This result combines information from the data, the cosmological prior, and the prior given by the range of the grid. The latter introduces a correlation between Σ and γ that arises from the fact that the uncorrelated grid on the (Q,η) plane of the phenomenological parameters maps onto a correlated grid on the (Σ,γ) plane, as shown in Fig. 3.
The corresponding onedimensional PDFs of both Σ and γ are centered on the GR values. This happens both for the marginalised cases and when fixing one of the parameters at the GR value. This shows that the amplitudes of the S10 COSMOS cosmic shear measurements agree with those of the WMAP 7year CMB measurements. They both show a preference for the same value of σ_{8}. In other words, given the strong priors of WMAP the amplitude of the cosmic shear models is set, with the prefered ones lying close to GR. This can also be seen as an indirect test on the amplitude of systematics. Indeed, if there were systematics not accounted for in the data, or in the estimation of the redshift distribution, they would make the signal less compatible with the WMAP priors and could show up in the results as evidence for modified gravity. Figure 4 shows the PDF of Σ, marginalised over γ, which corresponds to a 1σ constraint of Σ = 1.02 ± 0.38.
Fig. 4 Marginalised onedimensional PDF of Σ. 

Open with DEXTER 
The finding that these postGR parameters are centered on the GR values is quite robust, making use of the information on the redshift evolution of the amplitude from the tomographic correlations. By repeating the analysis using only the wide bin autocorrelations (which contains 60% of the galaxies) we found a slight change in the (σ_{8},Σ) degeneracy, shown in Fig. 5, and a shift on the mean of the posterior distribution of Σ to .
Fig. 5 Marginalised likelihood contours on the (Σ,σ_{8}) plane at 68% and 95% confidence levels, obtained from the full analysis (filled contours), and using only the “66” correlation (open, dashed lines). 

Open with DEXTER 
6. Summary and conclusions
We used the new COSMOS cosmic shear data to test possible deviations from GR. The departure from GR is parametrized in a phenomenological modelindependent way extensively used in the literature. The two phenomenological parameters allow to derive two other parameters that consistently parametrize the modifications to twopoint cosmic shear amplitude and growth, and are scaleindependent and have a scaling redshift dependence with the ΛCDM background.
A more realistic parametrization should take into account three regimes of modified gravity (Hu & Sawicki 2007), allowing for compatibility with the background expansion on large scales and agreement with GR on subgalactic scales. Hu & Sawicki (2007) give a fitting formula for the nonlinear matter power spectrum across the regimes. That formula was confirmed for particular classes of theories by perturbation theory (Koyama et al. 2009) and Nbody simulations (Oyaizu et al. 2008). In particular for f(R), according to Koyama et al. (2009) it shows systematic deviations of at least 10% from the formula of Smith et al. (2003), for scales smaller than k ≈ 0.8 h M_{pc}^{1}. COSMOS has an average source redshift of z ≈ 1.5 and is mainly sensitive to lensing structures at around half comoving distance to the source. Therefore, this limiting scale corresponds to θ ≈ 20 arcmin and the COSMOS range is affected. We do not apply this result to our modelindependent study, but in general a systematic increase in the amplitudes of the models implies that cosmic shear becomes more sensitive to Σ, leading to a more precise determination of this parameter for the same data.
Our results use WMAP priors obtained for ΛCDM models. This is consistent with using a ΛCDM background evolution. The CMB power spectrum is affected by the Integrated SachsWolfe effect essentially at the recent epoch when the modified gravity effects become important and the potentials change significantly over time. This affects only the scales large enough compared with the photons crossing time. Since the effect is at low redshift those large physical scales correspond to very large angular scales. In this way, the peak structure is not changed and the parameters of modified gravity are not degenerate with the standard parameters at perturbation level, allowing to use ΛCDM priors.
We did not find evidence for inconsistency with GR, in agreement with other cosmological tests in the linear regime (see Jain & Khoury 2010, for a review). In particular, we do not find the evidence for deviation reported by Daniel & Linder (2010) when using smallscale CFHTLS data. Our result using the new tomographic correlations of COSMOS corroborates the results obtained by Daniel & Linder (2010) using a different COSMOS catalogue (Massey et al. 2007). We cannot however make a direct comparison of the constraints, since we assume a different redshiftdependence and different priors.
Our analysis relied on the tomographic amplitude of 2pt cosmic shear, and did not allow to break the Σ − γ degeneracy. The scaledependence of the modifications, which is more distinctive on scales larger than those probed by COSMOS, should be explored with future wide highprecision weak lensing surveys. Such an analysis will have a better chance of reaching our original goal, i.e., to investigate if given strong cosmological priors there would be room to compensate a modification of the predicted amplitude of the cosmic shear signal with a modification of the growth of matter perturbations.
Acknowledgments
It is a pleasure to thank the COSMOS team for producing the data used in this work, in particular Jan Hartlap, Benjamin Joachimi, Martin Kilbinger and Patrick Simon. We thank Tommaso Giannantonio for fruitful discussions and the anonymous referee for sharing his/her worries. I.T. is grateful to AIfA – University of Bonn and the Leiden Observatory for hospitality. I.T. is funded by an FCT (Portugal) fellowship and acknowledges support of the European Programmes “DUEL” RTN and “CALICE” ERG. T.S. acknowledges support from the Netherlands Organization for Scientific Research (NWO).
References
 Acquaviva, V., Hajian, A., Spergel, D. N., & Das, S. 2008, Phys. Rev. D, 78, 043514 [NASA ADS] [CrossRef] [Google Scholar]
 Amendola, L., Kunz, M., & Sapone, D. 2008, J. Cosmology Astropart. Phys., 4, 13 [NASA ADS] [CrossRef] [Google Scholar]
 Bartelmann, M., & Schneider, P. 2001, Phys. Rep., 340, 291 [NASA ADS] [CrossRef] [Google Scholar]
 Bean, R., & Tangmatitham, M. 2010, Phys. Rev. D, 81, 083534 [NASA ADS] [CrossRef] [Google Scholar]
 Beynon, E., Bacon, D. J., & Koyama, K. 2010, MNRAS, 403, 353 [NASA ADS] [CrossRef] [Google Scholar]
 Bertschinger, E., & Zukin, P. 2008, Phys. Rev. D, 78, 024015 [NASA ADS] [CrossRef] [Google Scholar]
 Caldwell, R., Cooray, A., & Melchiorri, A. 2007, Phys. Rev. D, 76, 023507 [NASA ADS] [CrossRef] [Google Scholar]
 Carroll, S., Duvvuri, V., Trodden, M., & Turner, M. S. 2004, Phys. Rev. D, 70, 043528 [NASA ADS] [CrossRef] [Google Scholar]
 Daniel, S. F., & Linder, E. V. 2010, Phys. Rev. D, 82, 103523 [NASA ADS] [CrossRef] [Google Scholar]
 Daniel, S. F., Caldwell, R. R., Cooray, A., & Melchiorri, A. 2008, Phys. Rev. D, 77, 103513 [NASA ADS] [CrossRef] [Google Scholar]
 Daniel, S. F., Caldwell, R. R., Cooray, A., Serra, P., & Melchiorri, A. 2009, Phys. Rev. D, 80, 023532 [NASA ADS] [CrossRef] [Google Scholar]
 Daniel, S. F., Linder, E. V., Smith, T. L., et al. 2010, Phys. Rev. D, 81, 123508 [NASA ADS] [CrossRef] [Google Scholar]
 Deffayet, C., Dvali, G., & Gabadadze, G. 2002, Phys. Rev. D, 65, 044023 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Doré, O., Martig, M., Mellier, Y., et al. 2007 [arXiv:0712.1599] [Google Scholar]
 Dvali, G., Gabadadze, G., & Porrati, M. 2000, Phys. Lett. B, 485, 208 [NASA ADS] [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Eisenstein, D. J., & Hu, W. 1998, ApJ, 496, 605 [NASA ADS] [CrossRef] [Google Scholar]
 Guzik, J., Jain, B., & Takada, M. 2010, Phys. Rev. D, 81, 023503 [NASA ADS] [CrossRef] [Google Scholar]
 Hartlap, J., Simon, P., & Schneider, P. 2007, A&A, 464, 399 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hilbert, S., Hartlap, J., White, S. D. M., & Schneider, P. 2009, A&A, 499, 31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hu, W., & Sawicki, I. 2007, Phys. Rev. D, 76, 104043 [NASA ADS] [CrossRef] [Google Scholar]
 Huterer, D., & Linder, E. V. 2007, Phys. Rev. D, 75, 023519 [NASA ADS] [CrossRef] [Google Scholar]
 Ilbert, O., Capak, P., Salvato, M., et al. 2009, ApJ, 690, 1236 [NASA ADS] [CrossRef] [Google Scholar]
 Jain, B., & Khoury, J. 2010, Ann. Phys., 325, 1479 [NASA ADS] [CrossRef] [Google Scholar]
 Jain, B., & Zhang, P. 2008, Phys. Rev. D, 78, 063503 [NASA ADS] [CrossRef] [Google Scholar]
 Knox, L., Song, Y.S., & Tyson, J. A. 2006, Phys. Rev. D, 74, 023512 [NASA ADS] [CrossRef] [Google Scholar]
 Koivisto, T., & Mota, D. F. 2006, Phys. Rev. D, 73, 083502 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Koyama, K., Taruya, A., & Hiramatsu, T. 2009, Phys. Rev. D, 79, 123512 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, A., & Bridle, S. 2002, Phys. Rev. D, 66, 103511 [NASA ADS] [CrossRef] [Google Scholar]
 Linder, E. V. 2005, Phys. Rev. D, 72, 043529 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Linder, E. V., & Cahn, R. N. 2007, Astropart. Phys., 28, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Lue, A., Scoccimarro, R., & Starkman, G. D. 2004, Phys. Rev. D, 69, 124015 [NASA ADS] [CrossRef] [Google Scholar]
 Ma, C.P., & Bertschinger, E. 1995, ApJ, 455, 7 [NASA ADS] [CrossRef] [Google Scholar]
 Massey, R., Rhodes, J., Leauthaud, A., et al. 2007, ApJS, 172, 239 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Mukhanov, V. F., Feldman H. A., & Brandenberger, R. H. 1992, Phys. Rep., 215, 206 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Nojiri, S., & Odintsov, S. D. 2003, Phys. Rev. D, 68, 123512 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Nojiri, S., & Odintsov, S. D. 2010 [arXiv:1011.0544] [Google Scholar]
 Oyaizu, H., Lima, M., & Hu, W. 2008, Phys. Rev. D, 78, 123524 [NASA ADS] [CrossRef] [Google Scholar]
 Pogosian, L., Silvestri, A., Koyama, K., & Zhao, G.B. 2010, Phys. Rev. D, 81, 104023 [NASA ADS] [CrossRef] [Google Scholar]
 Rapetti, D., Allen, S. W., Mantz, A., & Ebeling, H. 2010, MNRAS, 406, 1796 [NASA ADS] [Google Scholar]
 Schmidt, F. 2008, Phys. Rev. D, 78, 043002 [NASA ADS] [CrossRef] [Google Scholar]
 Schrabback, T., Hartlap, J., Joachimi, B., et al. 2010, A&A, 516, 63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Silvestri, A., & Trodden, M. 2009, Rep. Prog. Phys., 72, 096901 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Smith, R. E., Peacock, J. A., Jenkins, A., et al. 2003, MNRAS, 341, 1311 [NASA ADS] [CrossRef] [Google Scholar]
 Song, Y.S., & Doré, O. 2009, J. Cosmology Astropart. Phys., 3, 25 [NASA ADS] [CrossRef] [Google Scholar]
 Song, Y.S., Hu, W., & Sawicki, I. 2007, Phys. Rev. D, 75, 044004 [NASA ADS] [CrossRef] [Google Scholar]
 Song, Y.S., Hollenstein, L., CalderaCabral, G., & Koyama, K. 2010, J. Cosmology Astropart. Phys., 4, 18 [NASA ADS] [CrossRef] [Google Scholar]
 Thomas, S. A., Abdalla, F. B., & Weller, J. 2009, MNRAS, 395, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Uzan, J.P. 2007, Gen. Rel. Grav., 39, 307 [CrossRef] [Google Scholar]
 Uzan, J.P., & Bernardeau, F. 2001, Phys. Rev. D, 64, 083004 [NASA ADS] [CrossRef] [Google Scholar]
 Wang, L., & Steinhardt, P. J. 1998, ApJ, 508, 483 [NASA ADS] [CrossRef] [Google Scholar]
 Will, C. 2006, Liv. Rev. Rel., 9, 3 [Google Scholar]
 Zhang, P., Liguori, M., Bean, R., & Dodelson, S. 2007, Phys. Rev. Lett., 99, 141302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Zhao, G.B., Pogosian, L., Silvestri, A., & Zylberberg, J. 2009, Phys. Rev. D, 79, 083513 [NASA ADS] [CrossRef] [Google Scholar]
 Zhao, G.B., Giannantonio, T., Pogosian, L., et al. 2010, Phys. Rev. D, 81, 103510 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Description of the redshift bins used in the analysis, with their amount of contributing galaxies and mean redshifts.
All Figures
Fig. 1 Correlation function . Top panel shows the measurement (points with error bars), the bestfit GR model for a nontomographic analysis (Ω_{m} = 0.27, σ_{8} = 0.76) (solid), modified gravity models with the same standard parameters but (Σ = 1.6, γ = 0.55) (dashed), or (Σ = 1, γ = 0.65) (dotted). Bottom panel shows the deviation between the same dashed and dotted models with respect to the GR model. 

Open with DEXTER  
In the text 
Fig. 2 Marginalised onedimensional PDFs obtained for the phenomenological postGR parameters. 

Open with DEXTER  
In the text 
Fig. 3 The dotted lines show the limits of the probed region. The 68% C.L. likelihood contour essentially fills up that region, showing that for these data the range of models chosen constitutes a significant prior. 

Open with DEXTER  
In the text 
Fig. 4 Marginalised onedimensional PDF of Σ. 

Open with DEXTER  
In the text 
Fig. 5 Marginalised likelihood contours on the (Σ,σ_{8}) plane at 68% and 95% confidence levels, obtained from the full analysis (filled contours), and using only the “66” correlation (open, dashed lines). 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.