Growth index of matter perturbations in the light of Dark Energy Survey

We study how the cosmological constraints from growth data are improved by including the measurements of bias from Dark Energy Survey (DES). In particular, we utilize the biasing properties of the DES Luminous Red Galaxies (LRGs) and the growth data provided by the various galaxy surveys in order to constrain the growth index (γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}) of the linear matter perturbations. Considering a constant growth index we can put tight constraints, up to ∼10%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 10\%$$\end{document} accuracy, on γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}. Specifically, using the priors of the Dark Energy Survey and implementing a joint likelihood procedure between theoretical expectations and data we find that the best fit value is in between γ=0.64±0.075\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma =0.64\pm 0.075$$\end{document} and 0.65±0.063\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.65\pm 0.063$$\end{document}. On the other hand utilizing the Planck priors we obtain γ=0.680±0.089\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma =0.680\pm 0.089$$\end{document} and 0.690±0.071\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.690\pm 0.071$$\end{document}. This shows a small but non-zero deviation from General Relativity (γGR≈6/11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma _{\mathrm{GR}}\approx 6/11$$\end{document}), nevertheless the confidence level is in the range ∼1.3-2σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 1.3-2\sigma $$\end{document}. Moreover, we find that the estimated mass of the dark-matter halo in which LRGs survive lies in the interval ∼6.2×1012h-1M⊙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 6.2 \times 10^{12} h^{-1} M_{\odot }$$\end{document} and 1.2×1013h-1M⊙\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.2 \times 10^{13} h^{-1} M_{\odot }$$\end{document}, for the different bias models. Finally, allowing γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} to evolve with redshift [Taylor expansion: γ(z)=γ0+γ1z/(1+z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma (z)=\gamma _{0}+\gamma _{1}z/(1+z)$$\end{document}] we find that the (γ0,γ1)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\gamma _{0},\gamma _{1})$$\end{document} parameter solution space accommodates the GR prediction at ∼1.7-2.9σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim 1.7-2.9\sigma $$\end{document} levels.


Introduction
The past and present analysis of various cosmological data (SNIa, Cosmic Microwave Background-CMB, Baryonic Acoustic Oscillations-BAOs, Hubble parameter measurements etc) converge to the following cosmological paradigm, the observed Universe is spatially flat and the cosmic fluid consists of ∼ 4% luminous (baryonic) matter, ∼ 26% dark a e-mail: svasil@academyofathens.gr b e-mail: fotis-anagnostopoulos@hotmail.com (corresponding author) matter and ∼ 70% some sort of dark energy (hereafter DE) which plays a key role in explaining the accelerated expansion of the universe (cf. [2,13,15,24,31] and references therein). Despite the fact that there is an agreement among the majority of cosmologists concerning the ingredients of the cosmic fluid however, there are different explanations regarding the physical mechanism which causes the accelerated expansion of the universe. In brief, the general avenue that one can design in order to study cosmic acceleration is to treat DE either as a new field in nature or as a modification of General Relativity (see for a review [4,18,20]).
It has been proposed [7,38,57] that in order to discriminate between scalar field DE and modified gravity one may utilize the evolution of the linear growth of matter fluctuations δ m (z) = δρ m /ρ m . In particular, we introduce the growth rate of clustering, which is given by f (a) = dlnD dlna γ m (a), where D(a) = δ m (a)/δ m (a = 1) is the linear growth factor (normalized to unity at the present epoch), a(z) = (1 + z) −1 is the scale factor of the universe, m (a) is the dimensionless matter density parameter and γ is the so called growth index [45,65].
In fact the determination of the growth index is considered one of the main targets in these kind of studies because it can be used in order to test general relativity (GR) on extragalactic scales [28], even in a model independent fashion [43]. Indeed, in the literature one may find a large family of studies in which the functional form of the growth index is given analytically for several cosmological models namely, scalar field [11,26,30,39,42,52,66], f (R) [27,61], f (T ) [6] Finsler-Randers [12], running vacuum models [10], clustered and Holographic dark energy [41].
From the view point of large scale structure, the study of the distribution of matter on extragalactic using different mass tracers (galaxies, AGNs, clusters of galaxies etc) provides important constraints on theories of structure formation. Specifically, owing to the fact that gravity reflects, via gravitational instability, on the physics of clustering [45] it is natural to utilize the clustering/biasing properties of the extragalactic mass tracers in constraining cosmological models (see [8,36,40]) as well as to test the validity of GR on cosmological scales (see [11,32]). Following the above lines, in the current article we combine the linear bias data of Luminous Red Galaxies (hereafter LRGs; [23]), recently released by the DES group, with the growth rate data as provided by Sagredo et al. [48], in order to place constraints on (γ , M h ). Notice that M h is the dark matter halo in which LRG live.
The structure of the paper is as follows. In Sect. 2 we present the DES bias data and the growth data. In Sect. 3 we provide the family of basic bias models, while in Sect. 4 we discuss the evolution of linear matter fluctuations. The outcome of our analysis is presented in Sect. 5, while our main conclusions can be found in Sect. 6.

DESY1 red galaxies bias data and growth data
In a sequence of previous theoretical articles we have proposed to use the biasing properties of extragalactic sources in order to constrain the growth index of matter fluctuations [11]. Therefore, in the light of recent Dark Energy Survey (DES) bias data, we attempt to compare the predictions of the most popular linear bias models (see below) with the data. Specifically, the DES bias data [23] were extracted in the context of the angular correlation function (ACF) using the 1-year DES sample of ∼ 6.6 × 10 5 LRGs as tracers of the LSS. In Table 1 we list the numerical values of the DES bias data with the corresponding errors.
The aim of our work is the following: if we accept that the background expansion is given by the concordance CDM model then we are interested to check the growth index of matter fluctuations. Specifically, we restrict the present analysis to the most popular expansion models. First we utilize the DES/Planck/JLA/BAO CDM cosmology, namely m0 = 1 − 0 = 0.301, h = 0.682, b0 = 0.048, n = 0.973 with σ 8 = 0.801 [1] and second we use the Planck TT+TE+EE+low+lensing CDM cosmology, hence m0 = 1 − 0 = 0.3153, h = 0.6736, b0 = 0.0493, n = 0.9649, and σ 8 = 0.811 [2]. In this context, the normalized Hubble parameter of the CDM model is written as In addition to DES bias data, we use in our analysis the growth data and the corresponding covariances as collected by Sagredo et al. [48] (see their Table 1 and references therein). This sample contains 22 entries for which the product f (z)σ 8 (z) is available as a function of redshift, where f (z) is the growth rate of clustering. 1 It is well known that the product f σ 8 is almost a model-independent parametrization of expressing the observed growth history of the universe [55].

Bias models
Let us first briefly present the main bias models. In particular, from the so called merging bias family we include here the models of Sheth et al. [51], Jing [34] and De Simone et al. [21]. For these models the bias factor is written as a function of the peak-height parameter, ν = δ c (z)/σ (M h , z) where δ c is the linearly extrapolated density threshold above which structures collapse. Here we use the accurate fitting formula of Weinberg and Kamionkowski [67] to estimate δ c (z). Moreover, the mass variance is written as 3 is the top-hat smoothing kernel with R = [3M h /(4πρ m )] 1/3 , M h is the halo mass and ρ m is the present value of the mean matter density, namely ρ m 2.78 × 10 11 m M /Mpc 3 . The quantity P(k) is the CDM linear power spectrum given by P(k) = P 0 k n T 2 (k) where n is the spectral index of the primordial power spectrum and T (k) is the CDM transfer function provided by Eisenstein and Hu [22]: 1 By definition the estimator f (z)σ 8 (z) is independent from linear bias and thus it is not affected by the dark matter halo. Indeed, the observed growth rate of structure ( f obs (z) = βb) is derived from the redshift space distortion parameter β(z) and the linear bias. Observationally, using the anisotropy of the spatial correlation function one can estimate the β(z) parameter. The linear bias factor can be defined as the ratio of the variances of the tracer (galaxies, Luminous Red Galaxies etc) and underlying mass density fields, smoothed at 8h where σ 8,tr (z) is measured directly from the sample. Combining the above definitions we arrive at f σ 8 = βσ 8,tr , hence the growth rate data are not not affected from bias, implying that the two data sets used in the present analysis are not correlated.
(2) the normalization of the power spectrum becomes Below we provide some details concerning the bias models.

SMT and JING
Sheth et al. ( [51], hereafter SMT) based on the ellipsoidal collapse model they found the following bias formula with Using N-body simulations they evaluated the free parameters of the model, α = 0.707, b = 0.5, c = 0.6 Also, Jing [34] proposed the following bias form

DMR
De Simone et al. [21] (hereafter DMR) generalized the original Press-Schether formalism incorporating a non-Markovian extension with a stochastic barrier. In this model, the critical value for spherical collapse was assumed to be a stochastic variable, whose scatter reflects a number of complicated aspects of the underlying dynamics. Therefore, the bias factor is

BPR
In addition to merging bias models we shall use the generalized model of Basilakos et al. [11] (hereafter BPR). This form of bias is valid for any dark energy model including those of modified gravity. In this case, using the hydrodynamic equations of motion, linear perturbation theory and the Friedmann-Lemaitre solutions a second differential equation of bias is derived Basilakos et al. [11]. The solution of the differential equation is given by: with J (z) = z 0 1+y E(y) dy where b 0 is the bias factor at the present time. The integration constants b 0 and C 2 can be found in Basilakos et al. [11], namely (11) and Finally, we would like to mention that the performance of the above bias models in reproducing the DES bias data has been studied in Papageorgiou et al. [44]. At this point some comments are in order. The choice of linear bias, used in the current work, is dictated by the fact that the DES bias data were extracted in the linear bias regime. According to Krause et al. [35] the scale of ∼ 8h −1 Mpc utilized by the DES team implies that the impact of non-linear effects on clustering and thus on biasing is almost negligible (see also [23]). In this context, we use the common assumption that each extragalactic mass tracers is hosted by a dark matter halo of a given mass. Possible double tracer occupancy of the dark matter halo is ignored. It is interesting to mention here that the linear bias approach (also used by DES team) relates a mass tracer, being a galaxy, an AGN, a LRG or a cluster of galaxies, with a host dark matter halo within which the mass tracer forms and evolves. The bias models themselves follow the linear evolution of the host halo and not the internal evolution of the astrophysical processes of the tracer. Thus the assumption is that the effects of nonlinear gravity and hydrodynamics (merging, feedback mechanisms, etc.) can be ignored in the linear-regime (see [19,25]).
However, in our work following the notations of [44] we include in our bias analysis the nonlinear power spectrum of [60] (and references therein). Notice that the halo-fit model of [60] is applicable in the wave-number range of k < 30hMpc −1 and the redshift range of 0 ≤ z ≤ 10. In brief, the nonlinear power spectrum consists of two terms: is the dimensionless nonlinear power spectrum 2 (k) = k 3 P(k)/(2π 2 ). The first term in the above formula is the so called one-halo term which dominates the power spectrum at small scales and it describes the matter correlation within the same dark matter halo. Also, the second is the two-halo term, which dominates the power spectrum at large scales and it describes the matter correlation between two distinct halos (for more details see appendix).

Evolution of liner growth
In this section we discuss the main points of the linear growth of matter fluctuations via which the growth index, γ , enters in the current analysis. Focusing on sub-horizon scales the differential equation that governs the linear matter perturbations [39,56,61,62] and references therein) is given bÿ where ρ m ∝ a −3 is the matter density, G eff = G N Q(t) with G N being the Newton's gravitational constant, while the effects of modified gravity are encapsulated in the quantity Q(t). Of course for those DE models which adhere to General Relativity G eff reduces to G N , hence Q(a) = 1.
The solution of the aforementioned equation (13) is δ m ∝ D(a), where D(a) is the growth factor. For any type of gravity the growth rate of clustering is given by the following useful parametrization [39,45,65] and thus we have where m (a) = m0 a −3 /E 2 (a) and γ is the growth index. Notice that the growth factor is normalized to unity at the present epoch. Now, inserting the operator d/dt = H d/d ln a and Eq. (14) into Eq. (13) we arrive at Considering the concordance CDM model, namely Q(a) = 1 it is easy to show thaṫ where w(a) = −1. In this case the Hubble parameter H (a) = H 0 E(a), where E(a) is given by Eq. (1) and H 0 is the Hubble constant. 2 Generally speaking the growth index may not be a constant but rather evolve with redshift; γ ≡ γ (z). In this framework, 2 For the comoving distance and for the dark matter halo mass we use the traditional parametrization H 0 = 100hkm/s/Mpc. Of course, when we treat the power spectrum shape parameter we utilize h ≡h = 0.68 [1].
substituting Eq. (14) into Eq. (16) we find where the prime denotes derivative with respect to redshift.

The likelihood analysis
In this section we provide the statistical method that we adopt in order to constrain the growth index, presented in the previous section. We implement a standard χ 2 minimization analysis in order to constrains either the (γ , M h ) parameter space. Specifically, in our case the situation is as follows: (1) For the DES biasing cosmological probe we use where the various forms of the theoretical b th are given in Sect. 3. Notice that σ bi = σ 2 i + σ 2 z , where σ i and σ z = 0.075 are the uncertainties of the observed bias and redshift respectively (see Table 1). 2. Regarding the analysis of the growth-rate data we use where M = { f σ 8,obs (z 1 ) − f σ 8 (z 1 , p 2 ), . . . , f σ 8,obs (z n ) − f σ 8 (z n , p 2 )} and C −1 cov is the inverse covariance matrix [48]. The theoretical growth-rate is given by: The vectors p 1 and p 2 provide the free parameters that enter in deriving the theoretical expectations. The first vector includes the free parameters which are related to the expansion and the environment of the parent dark matter halo in which the LRGs DES galaxies live. Since we wish to perform a joint likelihood analysis of the two cosmological probes and owing to the fact that likelihoods are defined as L i ∝ exp(−χ 2 i /2), the overall likelihood function becomes which is equivalent to: Based on the above we will provide our results for each free parameter that enters in the p 1 vector. The uncertainty of each fitted parameter will be estimated after marginalizing one parameter over the other, providing as its uncertainty the range for which χ 2 (≤ 1σ ).
As a further quality measure over the fits, we have used the AIC [3] criterion, in a modified form that is appropriate for small data sets [37]. Considering Gaussian errors AIC is given by where N is the total number of data and k is the number of fitted parameters (see also [37]). Of course, a smaller value of AIC implies a better model-data fit. In order to test the performance of the different bias models in fitting the data we need to utilize the model pair difference, namely AIC = AIC model − AIC min . From one hand, the restriction AIC ≤ 2 indicates consistency between the two comparison models. On the other hand, the inequalities 4 < AIC < 7 indicate a positive evidence against the model with higher value of AIC model Burnham and Anderson [17], while the condition AIC ≥ 10 suggests a strong such evidence.

Observational constraints
Below, we provide a qualitative discussion of our constraints, giving the reader the opportunity to appreciate the new results of the current study.

Constant growth index
Here we focus on the 1 parametrization, which means that the parameter space contains the following free parameters (γ , M h ). The presentation of our constraints is provided in Table 2 for the case of DES/Planck/JLA/BAO reference cosmology (see Sect. 2). The Table includes the goodness of fit  statistics (χ 2  min , AIC), for the specific bias models. Also, in Fig. 1 we present the 1σ , 2σ and 3σ confidence contours in the (γ , M h ) plane.
In particular, we find: We observe that the aforementioned bias models provide very similar results (within 1σ errors) as far as the growth index is concerned. The corresponding best fit values show a small but non-zero deviation from the theoretically predicted value of GR γ GR ≈ 6/11 (see solid lines of Fig. 1), where the range of the confidence level is ∼ 1.3σ − 1.7σ . Such a small discrepancy between the predicted and observationally fitted value of γ has also been discussed by other authors. For example recently, Zhao et al. [70] found γ = 0.656 +0.042 −0.046 , while [69] obtained γ = 0.628 +0.036 −0.039 . Also, similar results can be found in previous papers [14,29,49,68] in which the tension can reach to ∼ 2.5σ .
Furthermore, we find that the best bias model is the SMT, however the inequality AIC ≤ 2 indicates that the SMT bias model is statistically equivalent with rest of the models. The second result is that the differences of the bias models are absorbed in the fitted value of the DM halo mass in which LRGs live, and which ranges from ∼ 6.2 × 10 12 h −1 M − 1.2 × 10 13 h −1 M , for the different bias models and in the case of DESY1COSMO bias. As it can also be seen from Table 2, our derived mass of the host DM halo mass is consistent with that of Papageorgiou et al. [44], while Sawangwit et al. [50] and Pouri et al. [47] found M h (1.9 − 2) × 10 13 h −1 M .
In order to complete the present investigation we repeat the likelihood procedure in the case of Planck TT+TE+EE+ low+lensing CDM cosmology, hence m0 = 1 − 0 = 0.3153, h = 0.6736, b0 = 0.0493, n = 0.9649, and σ 8 = 0.811 [2]. Specifically, for the explored bias models we obtain (see also Table 2 Table 2 Observational constraints for the joint analysis of bias (see Table 1) and growth rate data: The 1st column shows the expansion model (see Sect. 2; Abbott et al. [1]), the 2nd column indicates the bias models (see Sect. 3), the 3rd column corresponds to γ and the 4 rth provides the fitted DM halo mass. The remaining columns present the goodness-of-fit statistics χ 2 min , AIC and AIC= AIC ,i − AIC min . Obviously, our statistical results remain quite robust (within 1σ ) against the choice of the undelying expansion [1,2]. Moreover, as it can be seen from Fig. 2 the growth index of the Planck TT + TE + EE + low + lensing CDM cosmology deviates with respect to that of GR (γ GR ≈ 6/11) at ∼ 1.5 − 2σ levels. Fig. 2 The iso-likelihood contours for 1σ -2σ -3σ levels in the (γ , M h ) parameter space for different bias models. Upper row: From left to right, JING and SMT models. Lower row: DMR and BPR models. For further details regarding the models, please see the relevant section. Notice that we use the Planck TT+TE+EE+low+lensing CDM cosmology [2]. The best fit values are given in Table 2. The vertical dashed line corresponds to γ GR ≈ 6/11

Constraints on γ (z)
In this section we implement the overall likelihood procedure in the (γ 0 , γ 1 ) parameter space. Based on the considerations discussed in the previous section the statistical vector takes the form p 1 In Fig. 3 we plot the results of our statistical analysis in the (γ 0 , γ 1 ) plane for the SMT bias model, since we have verified that using the other bias models we get similar contours. The predicted (γ (th) 0 , γ (th) 1 ) CDM values are indicated by the solid point, while the star corresponds to our best fit values. In brief for the DES/Planck/JLA/BAO cosmology we find γ 0 = 0.630 ± 0.072, γ 1 = 0.040 ± 0.403 with Fig. 3 Iso-likelihood contours for χ 2 = −2lnL/L max equal to 2.30, 6.18 and 11.83, corresponding to 1σ , 2σ and 3σ confidence levels in the (γ 0 , γ 1 ) plane in the case of 1 parametrization. The bias model is that of SMT, while the star corresponds to the best-fit point and the dot to the theoretical CDM point (γ 0 , γ 1 ) = (0.556, −0.011). In the right panel we present the contours that obtained using the Planck TT + TE + EE + low + lensing CDM cosmology [2], while in the left the contours obtained using DES/Planck/JLA/BAO CDM cosmology χ 2 min = 15.033 (AIC=19.533), while in the case of Planck TT + TE + EE + low + lensing CDM cosmology we get γ 0 = 0.670 ± 0.073, γ 1 = 0.100 ± 0.422 with χ 2 min = 14.988 (AIC=19.488). Notice that for the sake of simplicity we have marginalized the likelihood analysis over the LRG dark matter halo, namely log(M h / h −1 M ) = 13.00 and 13.07 respectively (see SMT model in Table 2).
We conclude that the joint statistical analysis put tight constraints γ 0 , however for γ 1 the corresponding error bars remain quite large. Also the range of deviation from GR is 1.7 − 2.9σ . We argue that with the next generation of data (mainly from Euclid) we will be able to test whether the growth index of matter fluctuations depends on time.

Conclusions
In the context of the concordance CDM model, testing the validity of general relativity (GR) on extragalactic scales is considered one of the most important tasks in cosmological studies, hence it is crucial to minimize the amount of priors needed to successfully complete such an effort. One such prior is the growth index (γ ) of matter perturbations. It is well known that a necessary step toward testing GR is to measure γ at the ∼ 1% accuracy level. Obviously, in order to control the systematic effects that possibly affect individual methods and tracers of the growth of matter perturbations we need to have independent estimations of γ .
Despite the fact that the CDM+GR model is found to be in a very good agreement with the majority of cosmological data [2], nonetheless the model seems to be currently in tension with some recent measurements, 3 related with the Hubble constant H 0 and the present value of the mass variance at 8h −1 Mpc σ 8 . Whether these tensions are the result of yet unknown systematic errors or hint some underlying 3 For recent reviews see [54,64]. new Physics is still unclear. In the light of the latter results, an intense debate is taking place in the literature and the aim of the present article is to contribute to this debate.
In this work we used the biasing properties of the Luminous Red Galaxies, recently released by the group of Dark Energy Survey (DES), together with growth rate data in order to constrain the growth index of matter perturbations. Specifically, in the framework of concordance +GR cosmology, we study the ability of four bias models to fit the DES bias data. Then we combined bias in a joint analysis with the growth rate of matter fluctuations to place constraints on the parameters.
Considering a constant growth index we placed constraints, up to ∼ 10% accuracy, on the growth index. Specifically, using the priors of the Dark Energy Survey we found that the constraints remain mostly unaffected by using different forms of bias. In particular, we obtained 0.640 ± 0.071, 0.650 ± 0.063, 0.640 ± 0.066 and γ = 0.640 ± 0.075 for SMT [51], JING [34], DMR [21] and BPR [11] bias models. Also utilizing the Planck priors we got γ = 0.680 ± 0.076, γ = 0.690 ± 0.071, γ = 0.680 ± 0.075 and γ = 0.680 ± 0.089 for the aforementioned bias factors. Obviously, we found a small but non-zero deviation from GR (γ GR ≈ 6/11), where the confidence level lies in the interval ∼ 1.3σ − 2σ . Such a small discrepancy between the predicted and observationally fitted value of γ has also been reported in several studies [14,29,49,[68][69][70]. Moreover, the intrinsic differences of the bias models are absorbed in the fitted value of the dark-matter halo mass in which LRGs survive, and which belongs in the range ∼ 6.2 × 10 12 h −1 M − 1.2 × 10 13 h −1 M .
Under the assumption that the growth index varies with time, namely γ (z) = γ 0 + γ 1 z/(1 + z), we showed that the (γ 0 , γ 1 ) parameter solution space accommodates the GR (γ 0 , γ 1 ) values at ∼ 1.7σ (∼ 2.9σ ) level utilizing the DES/Planck/JLA/BAO (Planck) priors. Similar to previous studies, we placed tight constraints on γ 0 , however the corresponding uncertainties of γ 1 remain large. The next generation (mainly from Euclid) of dynamical data are expected to improve the constraints on γ 1 , hence the validity of general relativity on extragalactic scales will be effectively checked. vide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .