## Global properties of the growth index of matter inhomogeneities in the Universe

- R. Calderon , D. Felbacq , R. Gannouji , D. Polarski , A.A. Starobinsky 1 ∗ 1 † 2 ‡ 1 § 3 4 , ¶ 1 Laboratoire Charles Coulomb, Universit´ Montpellier &amp; CNRS e
- UMR 5221, F-34095 Montpellier, France
- 2 Instituto de F´ ısica, Pontificia Universidad Cat´lica de Valpara´ ıso, o Casilla 4950, Valpara´ ıso, Chile 3
- 4 Moscow 101000, Russia
- Landau Institute for Theoretical Physics RAS, Moscow, 119334, Russia National Research University Higher School of Economics,

September 6, 2019

## Abstract

We perform here a global analysis of the growth index γ behaviour from deep in the matter era till the far future. For a given cosmological model in GR or in modified gravity, the value of γ (Ω m ) is unique when the decaying mode of scalar perturbations is negligible. However, γ ∞ , the value of γ in the asymptotic future, is unique even in the presence of a nonnegligible decaying mode today. Moreover γ becomes arbitrarily large deep in the matter era. Only in the limit of a vanishing decaying mode do we get a finite γ , from the past to the future in this case. We find further a condition for γ (Ω m ) to be monotonically decreasing (or increasing). This condition can be violated inside general relativity (GR) for varying w DE though generically γ (Ω m ) will be monotonically decreasing (like ΛCDM), except in the far future and past. A bump or a dip in G eff can also lead to a significant and rapid change in the slope dγ d Ω m . On a ΛCDM background, a γ substantially lower (higher) than 0 55 . with a negative (positive) slope reflects the opposite evolution of G eff . In DGP models, γ (Ω m ) is monotonically increasing except in the far future. While DGP gravity becomes weaker than GR in the future and w DGP → -1, we still get γ DGP ∞ = γ Λ CDM ∞ = 2 3 . In contrast, despite G DGP eff → G in the past, γ does not tend to its value in GR because dG DGP eff = 0.

/negationslash

∗ email:rodrigo.calderon-bruni@umontpellier.fr

$$\stackrel { \nu, \infty, \infty } { \frac { d G _ { \text{eff} } ^ { D G P } } { d \Omega _ { m } } } \Big | _ { - \infty } \frac { \bar { 7 } } { 7 } \\ \text{elier.fr}$$

† email:didier.felbacq@umontpellier.fr

‡ email:radouane.gannouji@pucv.cl

§ email:david.polarski@umontpellier.fr

¶ email:alstar@landau.ac.ru

## 1 Introduction

Despite intensive activity in recent years, the late-time accelerated expansion rate of the Universe remains a theoretical challenge. A wealth of theoretical models and mechanisms were put forward for its solution, see the reviews [1]. Remarkably, the simplest model - GR with a non-relativistic matter component (mainly a non-baryonic one) and a cosmological constant Λ - is in fair agreement with observational data, especially on large cosmic scales. Notwithstanding the theoretical problems it raises, this model provides a benchmark for the assessment of other dark energy (DE) models phenomenology. One can make progress by exploring carefully the phenomenology of the proposed models and comparing it with observations [2]. Tools which can efficiently discriminate between models, or between classes of models (e.g. [3]) are then needed. The growth index γ which provides a representation of the growth of density perturbations in dust-like matter, is an example of such phenomenological tool. Its use was pioneered long time ago in order to discriminate spatially open from spatially flat universes [4] and then generalized to other cases [5]. It was later revived in the context of dark energy [6], with the additional incentive to single out DE models beyond GR. The growth index has a clear and important signature in the presence of Λ: it approaches 6/11 for z /greatermuch 1 and it changes little from this value even up to z = 0. This behaviour still holds for smooth noninteracting DE models inside GR with a constant equation of state w DE . A strictly constant γ is very peculiar [7, 8]. On the other hand, this behaviour is strongly violated for some models beyond GR, see e.g. [9, 10]. Surprisingly, important properties of the growth index behaviour can be understood by making a connection with a strictly constant γ . Further, while present observations probe low redshifts z /lessorsimilar 2, still as we will see, more insight is gained looking at the global evolution of γ . Before starting this analysis, we first review the basic formalism for the study of the growth index γ .

## 2 The growth index

In this section, we review the basic equations and definitions concerning the growth index γ . Let us consider the dynamics on cosmic scales at the linear level of density perturbations δ m = δρ m /ρ m in the dust-like matter component. Deep inside the Hubble radius their evolution obeys the following equation

$$\ddot { \delta } _ { m } + 2 H \dot { \delta } _ { m } - 4 \pi G \rho _ { m } \delta _ { m } = 0$$

where H t ( ) ≡ ˙ ( a t /a t ) ( ) stands for the Hubble parameter and a t ( ) is the scale factor of a Friedmann-Lemaˆ ıtre-Robertson-Walker (FLRW) Universe filled with standard dust-like matter and DE (we neglect radiation at the matter and DE dominated stages).

For vanishing spatial curvature, we have for z /lessmuch z eq

$$h ^ { 2 } ( z ) = \Omega _ { m, 0 } \ ( 1 + z ) ^ { 3 } + ( 1 - \Omega _ { m, 0 } ) \ \exp \left [ 3 \int _ { 0 } ^ { z } d z ^ { \prime } \ \frac { 1 + w _ { D E } ( z ^ { \prime } ) } { 1 + z ^ { \prime } } \right ] \, \quad \ \ ( 2 )$$

with h z ( ) ≡ H H 0 , w DE ( z ) ≡ p DE ( z /ρ ) DE ( z ), z = a 0 a -1, and finally Ω m, 0 ≡ ρ m, 0 ρ cr 0 , . Equality (2) will hold for all FLRW models inside GR. It is still valid for many models beyond GR with appropriate definition of the dark energy sector. We recall the definition

$$\Omega _ { m } = \Omega _ { m, 0 } \frac { a _ { 0 } ^ { 3 } } { a ^ { 3 } } h ^ { - 2 }$$

and the useful relation

$$\intertext { u o n } w _ { D E } = \frac { 1 } { 3 ( 1 - \Omega _ { m } ) } \ \frac { d \ln \Omega _ { m } } { d \ln a } \.$$

Introducing the growth function f ≡ d ln δ m d ln a and using (4), equation (1) leads to the equivalent nonlinear first order equation [11]

$$\frac { d f } { d N } + f ^ { 2 } + \frac { 1 } { 2 } \left ( 1 - \frac { d \ln \Omega _ { m } } { d N } \right ) f = \frac { 3 } { 2 } \ \Omega _ { m }$$

with N ≡ ln a . Clearly f = p if δ m ∝ a p (with p constant). In particular, when the growing mode dominates, f → 1 in ΛCDM for large z and f = 1 in the Einsteinde Sitter universe. In the peculiar situation where the decaying mode dominates, f → -3 2 in ΛCDM for large z and f = -3 2 in the Einstein-de Sitter universe. We will return later to the consequences of a nonnegligible decaying mode. Also, if the dependence f ( N ) is known from observations, the Hubble function H z ( ) and finally the scale factor a t ( ) dependence can be determined unambiguosly in an analytical form [12] (see [13] for implementation of this to recent observational data).

The following parametrization has been intensively used and investigated in the context of dark energy

$$f = \Omega _ { m } ( z ) ^ { \gamma ( z ) }$$

where γ is dubbed growth index, though in general γ z ( ) is a genuine function. It can even depend on scales for models where the growth of matter perturbations has a scale dependence. The representation (6) has attracted a lot of interest with the aim to discriminate between DE models based on modified gravity theories and the ΛCDM paradigm. It turns out that the growth index is quasi-constant for the standard ΛCDM from the past untill today. Such a behaviour holds for smooth non-interacting DE models when w DE is constant, too [7].

In many DE models outside GR the modified evolution of matter perturbations is recast into

$$\ddot { \delta } _ { m } + 2 H \dot { \delta } _ { m } - 4 \pi G _ { \text{eff} } \rho _ { m } \delta _ { m } = 0$$

where G eff is some effective gravitational coupling appearing in the model. For example, for effectively massless scalar-tensor models [14], G eff is varying with time but it has no scale dependence while its value today is equal to the usual Newton's constant G . Introducing for convenience the quantity

$$g \equiv \frac { G _ { \text{eff} } } { G }$$

equation (7) is easily recast into the modified version of Eq. (5), viz.

$$\frac { d f } { d N } + f ^ { 2 } + \frac { 1 } { 2 } \left ( 1 - \frac { d \ln \Omega _ { m } } { d N } \right ) f = \frac { 3 } { 2 } \ g \ \Omega _ { m }$$

where the same GR definition Ω m = 8 πGρ m 3 H 2 is used. Some subtleties can arise if the defined energy density of DE becomes negative. When the growth index γ is strictly constant, it is straightforward to deduce from (9) that w DE = w with

$$\overline { w } \ = \ - \frac { 1 } { 3 ( 2 \gamma - 1 ) } \ \frac { 1 + 2 \Omega _ { m } ^ { \gamma } - 3 g \Omega _ { m } ^ { 1 - \gamma } } { 1 - \Omega _ { m } }$$

$$\equiv \ - \frac { \theta ^ { \cdot } _ { 1 } \, ^ { \prime } } { 3 ( 2 \gamma - 1 ) } \ F ( \Omega _ { m } ; g, \gamma ).$$

The quantity w (Ω m , γ ) defines a function in the (Ω m , γ ) plane which will be useful even when γ is not constant. The case g = 1 reduces to GR and we will simply write

$$F ( \Omega _ { m } ; g = 1, \gamma ) \equiv F ( \Omega _ { m } ; \gamma ) \.$$

Below, for any quantity v , v ∞ , resp. v -∞ , will denote its (limiting) value for N →∞ in the DE dominated era (Ω m → 0), resp. N →-∞ (generically Ω m → 1). We have in particular from (10) for g = 1 (GR)

$$\gamma \ = \ \frac { 3 \overline { w } _ { \infty } - 1 } { 6 \overline { w } _ { \infty } } \quad \ \quad \ \ ( 1 3 )$$

$$\gamma \ = \ \frac { 3 ( 1 - \stackrel { \sim } { w } _ { - \infty } ) } { 5 - 6 \overline { w } _ { - \infty } } \.$$

Here we assume that w &lt; 0 to get matter dominated stage in the past and dark energy dominated stage in the future. In addition, Eq. (13) requires that w ∞ &lt; -1 3 to have 0 &lt; γ &lt; 1, otherwise w ∞ becomes infinite.

As it was found in [8], these relations between a constant γ and the corresponding asymptotic values w ∞ and w -∞ apply also for the dynamical γ obtained for an arbitrary but given w DE . In the latter case, we obtain for g = 1

$$\gamma _ { \infty } \ = \ \frac { 3 w _ { \infty } - 1 } { 6 w _ { \infty } } \quad _ { \theta } \quad,$$

$$\gamma _ { - \infty } \ = \ \frac { 3 ( 1 - w _ { - \infty } ) } { 5 - 6 w _ { - \infty } }$$

with w ∞ , resp. w -∞ , the asymptotic value of w DE in the future, resp. past. We will return later to this important property.

## 3 Global analysis in the (Ω m, γ ) plane

It is natural to consider Ω m as the fundamental integration variable. In this case the evolution equation for γ obtained from (9) using (6) yields

$$2 \alpha \Omega _ { m } \ln ( \Omega _ { m } ) \frac { d \gamma } { d \Omega _ { m } } + \alpha ( 2 \gamma - 1 ) + F ( \Omega _ { m } ; g, \gamma ) = 0$$

where we have defined

$$\alpha \equiv 3 w _ { D E } \.$$

We restrict ourselves in this section to GR i.e. g = 1. This equation is well defined provided w DE is a given function of Ω m only. In particular (17) can be readily used for constant w DE . For arbitrary w DE ( a ) we can either express a as a function of Ω m or write (17) using the integration variable a . Note that the function a (Ω m ) is one-to-one for w DE &lt; 0. We consider here all the solutions to equation (17) on the entire Ω m interval. Each integral curve of equation (17) is the envelope of its tangent vectors

$$\left ( \begin{array} { c } 1 \\ \frac { d \gamma } { d \Omega _ { m } } \end{array} \right )$$

The collection of all these tangent vectors defines a vector field that can be written

$$\left ( \begin{array} { c } 2 \alpha \Omega _ { m } ( 1 - \Omega _ { m } ) \ln ( \Omega _ { m } ) \\ - \alpha ( 2 \gamma - 1 ) ( 1 - \Omega _ { m } ) - \tilde { F } ( \Omega _ { m } ; \gamma ) \end{array} \right )$$

where we have defined

$$\tilde { F } ( \Omega _ { m } ; \gamma ) \equiv ( 1 - \Omega _ { m } ) F ( \Omega _ { m } ; \gamma ) = 1 + 2 \Omega _ { m } ^ { \gamma } - 3 \Omega _ { m } ^ { 1 - \gamma } \.$$

For convenience, we have written the vector field in this way in order to have explicitly regular functions, so this vector field is well-defined everywhere for (Ω m , γ ) ∈ [0 , 1] × R . The associated unit vector field is represented in fig.(1) for w DE = -1.

The integral curves of this vector field (i.e. the phase portrait), are obtained by solving the autonomous differential system

$$\frac { d \Omega _ { m } } { d s } \ & = \ 2 \alpha \Omega _ { m } ( 1 - \Omega _ { m } ) \ln ( \Omega _ { m } ) \\ \frac { d \gamma } { d s } \ & = \ - \alpha ( 2 \gamma - 1 ) ( 1 - \Omega _ { m } ) - \tilde { F } ( \Omega _ { m } ; \gamma ) \\ - \. \quad. \quad. \quad. \quad.$$

where s ∈ R + is a dummy variable parametrizing the curves.

The vector field corresponding to w DE = -1 is displayed on fig.(1) and will be considered now for concreteness. The growth index γ (Ω m ) which starts at γ (1) ≡ γ -∞ = 6 11 is the only curve / γ (Ω m ) which is finite everywhere. This curve corresponds to the presence solely of the growing mode of eq.(1) or equivalently to the limit of a vanishing decaying mode. Let us consider this point in more details. Generally, in the presence of two modes δ ≡ δ 1 + δ 2 , we have

$$f = \frac { \delta _ { 1 } } { \delta } f _ { 1 } + \frac { \delta _ { 2 } } { \delta } f _ { 2 }$$

where f 1 &gt; 0, resp. f 2 &lt; 0, corresponds to the growing mode δ 1 , resp. decaying mode δ 2 . Inspection of (23) shows that f &lt; f 1 if δ δ 1 2 &gt; 0 while f &lt; f 1 or f &gt; f 1 when δ δ 1 2 &lt; 0. For concreteness we consider the situation around today. If both modes are positive we have f &lt; f 1 and hence γ &gt; γ 1 . As we go back in time (Ω m → 1), f → f 2 and so γ → ∞ when f ≈ 0. Let us consider now the situation δ 1 × δ 2 &lt; 0

Figure 1: The (unit time-oriented) vector field tangent to the solutions γ of eq.(17) is displayed for w DE = -1. The red curve corresponds to the solution for γ , with boundary condition γ (1) ≡ γ -∞ = 6 11. It is the only solution which is finite on the / entire Ω m interval, physically it corresponds to the presence of the growing mode of (1) only. The difference can be important in the far past only. All the other solutions correspond to solutions containing also the decaying mode solution. As explained in the text, these solutions (if physical) tend to the same value of γ , namely γ ∞ = 2 3 , for Ω m → 0, and they diverge for Ω m → 1. Some integral curves γ (Ω m ) are shown on figure (2).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00117v3-figure-38040b27-fc43-4cfa-a8e9-7bc9763ab81f.jpg)

with | δ 1 | &gt; δ | 2 | today. Now we have today f &gt; f 1 and therefore γ &lt; γ 1 . When we go backwards, f → ∞ at some point where the absolute value of both modes become equal, hence in that case γ → -∞ . These different situations are shown on figures (1) and (2).

Another interesting point concerns the asymptotic future (Ω m → 0). Inside GR, the solution to eq. (17) gives γ → γ ∞ , eq. (15), with

$$f \in \frac { 1 } { a ^ { 2 } H } \in a ^ { - \frac { 1 } { 2 } ( 1 - 3 w _ { \infty } ) } \to 0 \.$$

The crucial point is that the same asymptotic behaviour is obtained for all cases where the decaying mode is present, up to a change of the prefactor in (24) which depends on the initial conditions and on the amplitude of the decaying mode with respect to the growing mode. This immediately follows from the fact that one can neglect the last term in eq.(1) in this regime, so that

$$\delta _ { m } = \delta _ { \infty } + \text{const} \int \, \frac { d t } { a ^ { 2 } } = \delta _ { \infty } + \text{const} \int \, \frac { d N } { a ^ { 2 } H } \.$$

The same occurs for the more general modified gravity eq.(7), apart from the case of an anomalously large growth of G eff in the future. On the other hand we have in the future ( w ∞ &lt; -1 3 )

$$\Omega _ { m } & \sim a ^ { 3 w _ { \infty } } \.$$

Using the definition of γ , eq. (6), it is straightforward to obtain

$$\gamma = \frac { \ln f } { \ln \Omega _ { m } } \rightarrow \gamma _ { \infty }$$

for all curves with γ today smaller or larger than the value obtained when we integrate eq.(17) with γ (1) = γ -∞ (= 6 11 for w DE = -1). This is somehow complementary to the results obtained in [15] where the limit of a vanishing decaying mode was considered, however allowing for models beyond GR. Of course, in standard cosmological scenarios, the decaying mode is negligible already deep in the matter era so this result is essentially of mathematical interest. This is nicely illustrated with figure (2) for w DE = -1, all curves tend to the same limit γ ∞ = 2 3 but only one curve, corresponding to the limit of a purely growing mode, tends to the finite value γ -∞ = 6 11 in the past.

## 4 Global evolution of the slope

In this section we want to consider more closely the slope of γ (Ω m ) when γ evolves with time and this will allow us to find an interesting connection with the constant γ case. It is seen from (17) that a solution can satisfy dγ d Ω m = 0 in the ( γ, Ω ) plane m only on the curve Γ(Ω m ) defined as

$$\overline { w } = w _ { D E }$$

where w DE is the true equation of state (EoS) of the DE component of the system under consideration, while w is given in (10). We can view w (Ω m , γ ) as a given function of both variables Ω m and γ , while w DE is the EoS of DE given in function of Ω m . Hence in each value of Ω m , the curve Γ(Ω m ) takes that value of γ such that w = w DE , namely

$$\overline { w } \left ( \Omega _ { m }, \Gamma ( \Omega _ { m } ) \right ) = w _ { D E } ( \Omega _ { m } ) \.$$

The fact that Γ(Ω m ) is not constant just means that the solution γ (Ω m ) of (17) is not constant.

The curve Γ(Ω m ) defined by (28) satisfies by construction in the asymptotic future

$$\overline { w } \left ( 0, \Gamma ( 0 ) \right ) = w _ { \infty } \.$$

Hence we have in view of equalities (15) and (13)

$$\Gamma ( 0 ) = \frac { 3 w _ { \infty } - 1 } { 6 w _ { \infty } } = \gamma _ { \infty } \.$$

Analogously, we get in the asymptotic past

$$\overline { w } \left ( 1, \Gamma ( 1 ) \right ) = w _ { - \infty }$$

with, using (16) and (14)

$$\Gamma ( 1 ) = \frac { 3 ( 1 - w _ { - \infty } ) } { 5 - 6 w _ { - \infty } } = \gamma _ { - \infty } \.$$

It is exactly here that we can use the crucial property reminded at the end of Section 2. Indeed, as it was shown in [8], the solution γ (Ω m ) of (17) starts at γ -∞ in the past and tends to γ ∞ in the future (whence our choice of notation). In other words the curves Γ(Ω m ) and γ (Ω m ) meet at their endpoints (see figure (2)).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00117v3-figure-4fa25fe9-33cc-43ce-8470-ca07b54db139.jpg)

/negationslash

Figure 2: The curves γ (Ω m ) (red), with boundary condition γ (1) ≡ γ -∞ = 6 11, and / Γ(Ω m ) (purple), are shown for w DE = -1. Above the curve Γ(Ω m ), dγ d Ω m &gt; 0 while dγ d Ω m &lt; 0 below it. We see that Γ(Ω m ) is monotonically decreasing while the endpoints γ ∞ and γ -∞ , with γ ∞ &gt; γ -∞ , are identical for both curves. Therefore γ (Ω m ) cannot cross Γ(Ω m ) at 0 &lt; Ω m &lt; 1. The dashed lines are integral curves γ (Ω m ) however with γ (1) = γ -∞ = 6 11. / Still, as explained in the text while these curves diverge in the past, they all converge to the same limit γ ∞ in the future. Physically, all the dashed curves shown here correspond to density perturbations with a non negligible decaying mode. Note that those dashed curves which cross the curve Γ(Ω m ) satisfy dγ d Ω m = 0 where they cross.

It is also seen from (17) that all points ( γ, Ω ) above the curve Γ(Ω m m ) satisfy dγ d Ω m &gt; 0; on the other hand in the region below Γ(Ω m ) we have dγ d Ω m &lt; 0. Let us assume that (28) represents a monotonically decreasing function of Ω , hence m

γ ∞ &gt; γ -∞ . This holds for example for all constant EoS of DE inside GR. Then it is clear that γ (Ω m ) cannot cross the curve Γ(Ω m ) at some value γ 1 with γ -∞ &lt; γ 1 &lt; γ ∞ . If it did, γ (Ω m ) would cross with a zero derivative and have its minimum at γ 1 . This is in contradiction with γ -∞ &lt; γ 1 .

This also shows that a change in the sign of the slope of γ (Ω m ) is possible only if Γ(Ω m ) is not a monotonically decreasing function of Ω m . Note that the slope of γ (Ω m ) does not vanish at Ω m = 1 and Ω m = 0 even though γ (Ω m ) and Γ(Ω m ) meet there. Actually it even diverges in Ω m = 0 [8]. Indeed, we have that the prefactor of the first term of (17) vanish there leaving the slope a priori unspecified.

For a given growth function f , a lower Ω m implies a lower γ . For a given Ω m on the other hand, a decrease in f induces an increase in γ . All this is clearly understood from the equality γ = ln f ln Ω m . When γ (Ω m ) is monotonically descreasing, and this is the case for generic w DE inside GR, it is the second effect which prevails.

## 5 Systems with a non monotonic γ (Ω m )

We can consider now several cases where the growth index is not a monotonically decreasing function of Ω m .

## 5.1 Varying equation of state w DE

The first case that comes to one's mind is an oscillating equation of state w DE . Indeed, it is easy to show that Γ(Ω m ) is no longer monotonically decreasing in this case and that γ may oscillate provided the oscillations in w DE , and hence in Γ(Ω m ), are sufficiently pronounced. Observations however do not seem to favour such oscillations. We rather want to investigate whether a smoothly varying EoS can lead to a change in slope of γ . For this purpose we use the parametrization in terms of x ≡ a a 0 [16]

$$w _ { D E } = w _ { 0 } + w _ { a } ( 1 - x ) \.$$

We consider the cases with

$$w _ { - \infty } = w _ { 0 } + w _ { a } < 0 \quad \ \ \ w _ { a } > 0 \.$$

The first inequality is to enforce the domination of matter (and the disppearance of DE) in the past. We also restrict our attention to w a &gt; 0 in order to have a varying EoS leading to the opposite case in the future, namely DE domination and disappearance of matter. Also we will not impose the restriction w DE ≥ -1 valid for quintessence (a self-interacting scalar field minimally coupled to gravity) and admit arbitrarily negative values of w DE . We consider conservative values 0 &lt; w a ≤ 0 5 . and w 0 ≈ -1. We obtain that γ x ( ) is quasi-constant deep in the matter dominated stage and starts decreasing (with the expansion) in the past typically at z /similarequal 10, this decrease being stronger as w a is larger. In the future, a bump is obtained typically at z /similarequal -0 9, which is higher for lower . w a , this bump would essentially disappear for w a = 0 8. . Finally in the asymptotic future γ → 1 2 as we expect for w ∞ →-∞ .

γ

Figure 3: The evolution of γ is shown for CPL varying w DE . We see that γ (Ω m ) is essentially decreasing except in the remote future and, for w a /greaterorsimilar 0 2, in the far past. . Note the value γ ∞ = 1 2 as expected for w ∞ = -∞ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00117v3-figure-ce8158f6-923f-4259-9513-f98bb18b8529.jpg)

In function of Ω m , γ will be increasing deep in the DE domination typically till Ω m /similarequal 10 -5 , then decreasing till Ω m /similarequal 0 7, and finally increasing again for . w a /greaterorsimilar 0 2. . To summarize, γ (Ω m ) is essentially slowly varying and decreasing from the past to the future in this case covering in particular the range probed by observations. The significant changes in slope are pushed around Ω m ≈ 0. While this is interesting in itself, from an observational point of view models with w a /lessorsimilar 0 2 follow essentially the . ΛCDM phenomenology regarding their growth index. Strong departure from ΛCDM takes place in the remote future only. These features are shown on figure 3.

## 5.2 Beyond GR: a bump or a dip in g

There is however another way to have a change of slope of γ , and possibly in the range probed by observations. This is in the framework of modified gravity (see e.g. [17]) when g is substantially varying and this can take place even when w DE is almost constant and close to -1. We have to use the modified version of (17) with F (Ω m ; γ ) replaced by F (Ω m ; g, γ ) where g can be some arbitrary function, no longer equal to one as in GR. As we have already said earlier, a decreasing Ω m tends to decrease γ while a decreasing f tends to increase it, see eq.(27). A substantially varying g can work in both ways depending on how it boosts or dampens the growth and on the time this process takes place. In other words it can strongly affect f even if the background evolution remains 'standard'.

- a) Bump in g : A boost in g was already found some time ago for f ( R ) modified

gravity. Let us consider more closely this case. When the bump in g has reached its maximum ( ≈ 4 3 ) in the past, g is already starting on low redshifts its decrease towards its asymptotic value g ∞ = 1 while still being higher today than one. This explains the low value γ 0 ≈ 0 42 obtained in [9, 10] but also the large . negative slope dγ dz ≈ -0 25. . It would also be possible in principle in f ( R ) models to have a bump whose maximum is reached in the future with g today already substantially larger than one. In that case γ 0 could be again substantially lower than 0 55 however now . the slope on low redshifts would be positive, increasing with z . This provides a good example for which a global analysis gives more insight: γ today lower than γ Λ CDM with dγ dz &lt; 0 occurs if we have presently passed a bump in g .

b) Dip in g : While in f ( R ) models we have g ≥ 1, one can have modified gravity models with gravity weaker than GR. Then a decrease of g on low redshifts yields an increasing γ on these redshifts. This is in agreement with results found in [18]. In this case too, a global analysis gives more insight: on a background with w DE = -1, γ today higher than γ Λ CDM with dγ dz &gt; 0 tells us that we have presently passed a dip in g .

Roughly speaking, we have γ increasing, resp. decreasing when g is decreasing, resp. increasing, so they have opposite behaviours. We note finally a shift in the location of the extrema of g and of γ and as expected a dependence on the location and on the width of the bump in g . These properties are shown on figures 4.

Figure 4: a) Left-hand panel: a schematic representation is given when there is a bump, here chosen to be slightly in the past (we take Ω m, 0 = 0 3) while the background . expansion is fixed ( w DE = -1). We see that γ decreases first (with expansion) when g goes up before increasing when g decreases back to one. It is this last behaviour which was found today in many f ( R ) models. Hence the sign of the derivative tells us that we (just) passed the maximum of the bump. Note that γ bump and γ Λ CDM shown here have the same asymptotic values both in past and future because g → 1 in a smooth way and they have identical w DE = -1. The behaviour on small redshifts and in the future captures essentially the f ( R ) phenomenology. b) Right-hand panel: the same as on the left panel but now with a dip in g starting today. At the present time, γ &gt; γ Λ CDM and γ is decreasing with expansion reflecting that g is increasing and we have passed the minimum of the dip.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00117v3-figure-f199a270-8958-485f-9970-063f64b6a485.jpg)

## 5.3 Beyond GR: DGP model

Another interesting example is provided by DGP brane models [19]. Let us note first that in this model both g and w DE are given explicitly in terms of Ω m as follows [20]

$$w ^ { D G P } & = - [ 1 + \Omega _ { m } ] ^ { - 1 } & g ^ { D G P } & = 1 - \frac { 1 } { 3 } \frac { 1 - \Omega _ { m } ^ { 2 } } { 1 + \Omega _ { m } ^ { 2 } } \. & \quad \ ( 3 6 )$$

Hence it is straightforward to integrate γ (Ω m ) using the integration variable Ω m as in (17). It is seen from (36) that w DGP tends to w DGP ∞ = -1 in the future while g DGP tends to g DGP ∞ = 2 3 in the future. So in the future the expansion looks like ΛCDM, however cosmic perturbations feel a weaker gravity. In the past w DGP tends to w DGP -∞ = -1 2 while g DGP tends to g DGP -∞ = 1, so the gravitational driving force for the perturbations growth tends to its GR value. There are many interesting features when we solve for γ (Ω m ). Let us start with the behaviour in the future. As

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00117v3-figure-1d4c723c-c5ec-4220-8891-85efbc061372.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00117v3-figure-8eadc15b-0122-4460-9e67-bb0ea5b00418.jpg)

/negationslash

Figure 5: a)Left-hand panel: The behaviour of γ in DGP models is shown and compared to that in ΛCDM. We see that γ DGP &gt; Γ DGP and hence γ DGP (Ω m ) is an increasing function, except for the tiny interval Ω m /lessorsimilar 10 -3 . While w DGP -∞ = -1 2 , γ DGP -∞ = 11 16 and not 9 16 as one would have in GR. This is because the DGP model does not tend to GR in a smooth enough way, namely ( dg DGP d Ω m ) -∞ = 1 3 = 0. In the future however, while DGP does not tend to GR with g DGP ∞ = 2 3 and w DGP ∞ = -1, the same asymptotic value γ DGP ∞ = 2 3 is obtained as in GR for w ∞ = -1 because gravity in DGP is weaker than in GR in the far future. Note that all four curves tend to 2 3 in the future. b) Right-hand panel; The same in function of z . We see that a value very close to 11 16 is quickly reached in the past. As 11 16 &gt; 2 3 , Γ DGP (Ω m ) could be monotonically increasing, which it is except for the tiny range Ω m /lessorsimilar 10 -3 where γ DGP (Ω m ) is a decreasing function.

/negationslash gravity becomes weaker than GR in the future, γ DGP tends to γ GR ∞ corresponding to w GR ∞ = w DGP ∞ = -1. Note that this can still be so even if g grows but not too quickly, see the discussion in [15]. This can also be understood by looking at the asymptotic form of equation (17) for Ω m → 0. Indeed, the same equation is obtained in GR when g → constant, also for g = 1, and even more generally if g Ω 1 -γ m → 0. Hence we obtain

$$\gamma _ { \infty } ^ { D G P } = \frac { 2 } { 3 } \.$$

In other words the same future asymptotic value is obtained as for ΛCDM.

We turn now our attention to the behaviour in the past. Curiously, while g DGP -∞ = 1 we would naively expect to have the same relation between γ DGP -∞ and w DGP -∞ = -1 2 as between γ GR -∞ and w GR -∞ = -1 2 . However this is not so due to the non-vanishing derivative

$$\left ( \frac { d g ^ { D G P } } { d \Omega _ { m } } \right ) _ { - \infty } = \frac { 1 } { 3 } \.$$

In order to find γ (Ω m ), we can use the same method as in [8] and write eq.(17) for Ω m → 1. After some calculation, solving for γ (Ω m ) we obtain in the asymptotic past for the finite constant solution γ -∞ for arbitrary modified gravity models

$$\gamma _ { - \infty } = \frac { 3 \left ( w _ { - \infty } - 1 - d \right ) } { 6 w _ { - \infty } - 5 }$$

where we have set

$$d \equiv \left ( \frac { d g } { d \Omega _ { m } } \right ) _ { - \infty } \.$$

The GR result is recovered setting d = 0. Note that a similar relation was found between a constant γ corresponding to w -∞ = w -∞ in an arbitrary modified gravity model [8]. We finally obtain for the DGP model

$$\gamma _ { - \infty } ^ { D G P } = \frac { 1 1 } { 1 6 } \.$$

So the (finite) γ -∞ corresponding to w -∞ = -1 2 has the value γ DGP -∞ = 11 16 for a DGP model instead of 9 16 for GR. The derivation given here relies solely on the properties of the DGP model, eq. (36), and exhibits the origin of this anomalous value in an explicit way (see [6], [21] for other approaches).

This value for γ DGP -∞ corresponds to the constant γ DGP (and g = g DGP ) yielding identical equations of state, w -∞ = w DGP -∞ , in the past. We can summarize this in the following way

$$\gamma _ { - \infty } ^ { D G P } \left ( w _ { - \infty } ^ { D G P } \right ) = \gamma \left ( \overline { w } _ { - \infty } = w _ { - \infty } ^ { D G P } \right ) \. \quad \text{(4)} \\ \text{and side corresponds to a varying $\gamma - the true $\gamma$ obtained for $w^{DGP}$ and}$$

The left hand side corresponds to a varying γ - the true γ obtained for w DGP and g DGP - while the right hand side corresponds to a constant γ in a modified gravity model with g = g DGP . A similar equality holds in the future and we have shown it here with the DGP model

$$\gamma _ { \infty } ^ { D G P } \left ( w _ { \infty } ^ { D G P } \right ) = \gamma \left ( \overline { w } _ { \infty } = w _ { \infty } ^ { D G P } \right ) \. & & ( 4 3 ) \\ \text{allities were found in [8] for CR and imnlv that the curve for } \mathbb { R } \ \mathbb { R } \ \mathbb { R } \ \mathbb { R }$$

Analogous equalities were found in [8] for GR and imply that the curves Γ(Ω m ) and γ (Ω m ) meet at Ω m = 0 and Ω m = 1, a property that we have used in the previous section. We have generalized this result to a class of modified gravity models including DGP models.

As a corollary, we have the following important result: any modified gravity model with g (1) = 1, i.e. tending to GR in the past, however in a smooth way satisfying d = 0, will have the same limit γ GR -∞ ( w -∞ ) as in GR.

As we can see on figure (5), the curve Γ DGP (Ω m ) and γ DGP (Ω m ) intersect indeed at Ω m = 0 and Ω m = 1. In contrast to generic models inside GR, here we have γ DGP ∞ = 2 3 &lt; γ DGP -∞ = 11 16 . Hence Γ(Ω m ) cannot be monotonically decreasing, at best it would be monotonically increasing. If that were the case, γ DGP (Ω m ) would be monotonically increasing as well, always lying above Γ DGP . Actually this is what happens between Ω m /similarequal 10 -3 and Ω m = 1. It is only in the far future, Ω m /lessorsimilar 10 -3 , that the slope of γ DGP (Ω m ) is negative (i.e. γ DGP increases with expansion) and that the behaviour is similar to ΛCDM.

## 6 Summary and conclusion

The growth index γ allows to distinguish efficiently the phenomenology of dark energy (DE) models and has been used extensively for this purpose (see e.g. [22]). In this work we have performed an analysis of the growth index evolution from deep in the matter era till the asymptotic future. In this way global properties are exhibited. While from an observational point of view the main focus lies in the low redshift behaviour of the growth index γ still, a global analysis yields some interesting insight and results. Some of the properties found here become transparent when we perform a global analysis of the growth index evolution. We have shown that when the growth index had a bump, resp. dip, in the recent past while the background evolution is similar to ΛCDM, today it is substantially lower, resp. larger, than 0 55 . with a negative, resp. positive, slope dγ d Ω m reflecting that the gravitational coupling G eff of the underlying modified gravity model is already decreasing, resp. increasing, with the expansion. The behaviour with a bump is a schematic representation of many f ( R ) models [9, 10], the second case was considered in e.g. [18].

Using results valid for a constant growth index, we suggest a condition giving the global sign of the slope dγ d Ω m : when the curve Γ(Ω m ) introduced in Section 4 is monotonically decreasing then we have globally dγ d Ω m &lt; 0. This is the case in particular for models inside GR with a constant equation of state w DE = constant. Another interesting point concerns the value of the growth index for a given cosmological model at a given time. Actually the growth index γ can take a range of values. What is really meant by the value of γ (Ω m ) which corresponds to a given model is its value when the decaying mode of the perturbations mode tends to zero. We show that while the presence of a substantial decaying mode does not change the value of γ in the asymptotic future, this leads (as expected) to a divergence in the asymptotic past. It is only in the limit of a vanishing decaying mode that γ takes a finite value from γ -∞ in the asymptotic past (Ω m → 1) 6 11 for ΛCDM - up to γ ∞ in the asymptotic future (Ω m → 0) 2 3 for ΛCDM.

/negationslash

We have studied further the global behaviour of γ for the DGP model. In contrast to generic models in GR, we have γ DGP -∞ = 11 16 &gt; γ DGP ∞ = 2 3 with w DGP -∞ = -1 2 and w DGP ∞ = -1. While g DGP → 2 3 in the future, we have γ DGP ∞ = γ Λ CDM ∞ , so we get the same relation as in GR between γ DGP ∞ and w DGP ∞ = -1. Interestingly, while g DGP → 1 in the past so this model tends to GR in the past, γ DGP -∞ = 9 16 , the value expected in GR for w -∞ = -1 2 . This is because the DGP model does not tend to

/negationslash

GR in a way which is 'smooth enough'. Indeed, it satisfies dg DGP d Ω m -∞ = 1 3 = 0. As a corrolary we find that any modified gravity model which tends to GR in the past yields a γ -∞ which is the same function of w -∞ as in GR provided dg DGP d Ω m -∞ = 0. Finally we find that γ DGP (Ω m ) is monotonically increasing from the past until the far future (Ω m ≈ 10 -3 ) where it crosses the curve Γ(Ω m ) in accordance with the condition mentioned above. The results presented in this work indicate that a measurement of γ on a significant part of the expansion could give interesting constraints and a deeper insight into the physics governing the Universe dynamics.

## 7 Acknowledgements

A.A.S. was partly supported by the project number 0033-2019-0005 of the Russian Ministry of Science and Higher Education.

## References

- [1] V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. D 9 , 373 (2000); T. Padmanabhan, Phys. Rep. 380 , 235 (2003); P. J. E. Peebles and B. Ratra, Rev. Mod. Phys. 75 , 559 (2003); E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15 , 1753 (2006); V. Sahni and A. A. Starobinsky, Int. J. Mod. Phys. 15 , 2105 (2006); M. Li, X.-D. Li, S. Wang and Y. Wang, Commun. Theor. Phys. 56 , 525 (2011).
- [2] D. H. Weinberg, M. J. Mortonson, D. J. Eisenstein, C. Hirata, A. G. Riess and E. Rozo, Phys. Rept. 530 , 87 (2013); L. Amendola et al. , Living Rev. Rel. 16 , 6 (2013); P. Bull et al. , Phys. Dark Univ. 12 , 56 (2016).
- [3] V. Sahni, A. Shafieloo and A. A. Starobinsky, Astrophys. J. 793 , L40 (2014).
- [4] P. J. E. Peebles, Astrophys. J. 284 , 439 (1984).
- [5] O. Lahav, P. B. Lilje, J. R. Primack and M. J. Rees, MNRAS 251 , 128 (1991).
- [6] E. V. Linder and R. N. Cahn, Astropart. Phys. 28 481 (2007).
- [7] D. Polarski and R. Gannouji, Phys. Lett. B 660 , 439 (2008).
- [8] D. Polarski, A. A. Starobinsky and H. Giacomini, JCAP 1612 , 037 (2016).
- [9] R. Gannouji, B. Moraes and D. Polarski, JCAP 0902 , 034 (2009).
- [10] H. Motohashi, A. A. Starobinsky and J. Yokoyama, Progr. Theor. Phys. 123 , 887 (2010).
- [11] L. Wang and P. J. Steinhardt, Astrophys. J. 508 , 483 (1998).
- [12] A. A. Starobinsky, JETP Lett. 68 , 757 (1998).

- [13] B. L'Huillier, A. Shafieloo, D. Polarski and A. A. Starobinsky, arXiv:1906.05991.
- [14] B. Boisseau, G. Esposito-Far`se, D. Polarski and A. A. Starobinsky, Phys. Rev. e Lett. 85 , 2236 (2000).
- [15] E. V. Linder, D. Polarski, Phys. Rev. D 99 , 2, 023503 (2019)
- [16] M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10 , 213 (2001); E. V. Linder, Phys. Rev. Lett. 90 , 091301 (2003).
- [17] A. De Felice and S. Tsujikawa, Living Rev. Rel. 13 , 3 (2010); T. Clifton, P. G. Ferreira, A. Padilla and C. Skordis, Phys. Rept. 513 , 1 (2012); A. Joyce, B. Jain, J. Khoury and M. Trodden, Phys. Rept. 568 , 1 (2015); Mustapha Ishak, Living Rev. Rel. 22 , 1 (2019).
- [18] R. Gannouji, L. Kazantzidis, L. Perivolaropoulos and D. Polarski, Phys. Rev. D 98 , 10, 104044 (2018).
- [19] G. Dvali, G. Gabadadze and M. Porrati, Phys.Lett. B 485 , 208 (2000).
- [20] A. Lue, R. Scoccimarro and G. D. Starkman, Phys.Rev. D 69 , 124015 (2004)
- [21] Q.-G. Huang, Eur. Phys. J. C 74 , 2964 (2014).
- [22] V. Acquaviva, A. Hajian, D. N. Spergel, S. Das, Phys. Rev. D 78 Hao Wei, Phys. Lett. B664 1 (2008);
- S. Nesseris, L. Perivolaropoulos, Phys. Rev. D 77 , Phys.Rev. D 78 , 123010 (2008);
- 043514 (2008); 023504 (2008); Yungui Gong,
- Puxun Wu, Hong Wei Yu, Xiangyun Fu, JCAP 0906, 019 (2009);
- R. Bean, M. Tangmatitham, Phys. Rev. D 81 , 083534 (2010);
- Seokcheon Lee, Kin-Wang Ng, Phys. Lett. B 688 , 1 (2010);
- A. Bueno belloso, J. Garcia-Bellido, D. Sapone, JCAP 1110, 010 (2011);
- S. Nesseris, S. Basilakos, E.N. Saridakis, L. Perivolaropoulos, Phys. Rev. D 88
- 103010 (2013);
- K. Bamba, Antonio Lopez-Revelles, R. Myrzakulov, S.D. Odintsov, L. Sebastiani, Class. Quant. Grav. 30 015008 (2013);
- S. Basilakos, J. Sol`, Phys. Rev. D a 92 , no.12, 123501 (2015);
- I. de Martino, M. De Laurentis, S. Capozziello, Universe 1, no.2, 123 (2015);
- J. N. Dossett, M. Ishak, D. Parkinson, T. M. Davis, Phys.Rev. D 92 , no.2, 023003 (2015);
- A. B. Mantz et al., Mon. Not. Roy. Astron. Soc. 446, 2205 (2015);
- Alberto Bailoni, Alessio Spurio Mancini, Luca Amendola, arXiv:1608.00458;
- Xiao-Wei Duan, Min Zhou, Tong-Jie Zhang, arXiv:1605.03947;
- B. Wang, E. Abdalla, F. Atrio-Barandela, D. Pavon, Rept. Prog. Phys. 79, no.9, 096901 (2016);
- N. Nazari-Pooya, M. Malekjani, F. Pace, D. Mohammad-Zadeh Jassur, Mon. Not. Roy. Astron. Soc. 458, no.4, 3795 (2016);
- M. Malekjani, S. Basilakos, Z. Davari, A. Mehrabi, M. Rezaei, Mon. Not. Roy.

Astron. Soc. 464, 1192 (2017);

R. Gannouji, D. Polarski, Phys.Rev. D 98 , no.8, 083533 (2018)
R. Calderon D. Felbacq R. Gannouji D. Polarski A. A. Starobinsky
Global properties of the growth index of matter inhomogeneities in the universe
We perform here a global analysis of the growth index $\gamma$ behaviour from
deep in the matter era till the far future. For a given cosmological model in
GR or in modified gravity, the value of $\gamma(\Omega_{m})$ is unique when the
decaying mode of scalar perturbations is negligible. However,
$\gamma_{\infty}$, the value of $\gamma$ in the asymptotic future, is unique
even in the presence of a nonnegligible decaying mode today. Moreover $\gamma$
becomes arbitrarily large deep in the matter era. Only in the limit of a
vanishing decaying mode do we get a finite $\gamma$, from the past to the
future in this case. We find further a condition for $\gamma(\Omega_{m})$ to be
monotonically decreasing (or increasing). This condition can be violated inside
general relativity (GR) for varying $w_{DE}$ though generically
$\gamma(\Omega_{m})$ will be monotonically decreasing (like $\Lambda$CDM),
except in the far future and past. A bump or a dip in $G_{\rm eff}$ can also
lead to a significant and rapid change in the slope
$\frac{d\gamma}{d\Omega_{m}}$. On a $\Lambda$CDM background, a $\gamma$
substantially lower (higher) than $0.55$ with a negative (positive) slope
reflects the opposite evolution of $G_{\rm eff}$. In DGP models,
$\gamma(\Omega_{m})$ is monotonically increasing except in the far future.
While DGP gravity becomes weaker than GR in the future and $w^{DGP}\to -1$, we
still get $\gamma_{\infty}^{DGP}= \gamma_{\infty}^{\Lambda CDM}=\frac{2}{3}$.
In contrast, despite $G^{DGP}_{\rm eff}\to G$ in the past, $\gamma$ does not
tend to its value in GR because $\frac{dG^{DGP}_{\rm
eff}}{d\Omega_{m}}\Big|_{-\infty}\ne 0$.