Growth rate in the dynamical dark energy models

Dark energy models with a slowly rolling cosmological scalar field provide a popular alternative to the standard, time-independent cosmological constant model. We study the simultaneous evolution of background expansion and growth in the scalar field model with the Ratra–Peebles self-interaction potential. We use recent measurements of the linear growth rate and the baryon acoustic oscillation peak positions to constrain the model parameter \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}α that describes the steepness of the scalar field potential.


Introduction
Cosmological observations now convincingly show that the expansion of the Universe is accelerating [1][2][3][4]. One of the possible explanations of this empirical fact is that the energy density of the Universe is dominated by the so-called dark energy (DE) [5,6], a component with effective negative pressure.
The simplest DE candidate is a time-independent cosmological constant , and the corresponding cosmological model, the so-called CDM model, is considered to be a concordance model. This simple model, however, suffers from fine tuning and coincidence problems [7,8]. In the attempt of constructing a more natural model of DE many alternative scenarios have been proposed [9][10][11][12][13][14][15].
One of the alternatives to a cosmological constant are the models of a dynamical scalar field. In these models a spatially uniform cosmological scalar field, slowly rolling down its almost flat self-interaction potential, plays the role of a time-dependent cosmological constant. This family of models avoid the fine tuning problem, having a more natural explanation for the observed low energy scale of DE [17,18,[39][40][41]. For the scalar field models (the so-called a e-mail: olga.avsajanishvili@iliauni.edu.ge φCDM model) the equation of state P φ = wρ φ (with P φ and ρ φ the pressure and energy density of the scalar field) is time dependent, w = w(t), and unlike the cosmological constant, w(t) = −1, although at late-times it approaches −1. When the scalar field energy density starts to dominate the energy budget of the Universe, the Universe expansion starts to accelerate [19,20]. Even though at low redshifts the predictions of the model are very close to the ones of the cosmological constant, the two models ( CDM and the dynamical DE model) predict different observables over a wide range of redshifts.
The scalar field models can be classified via their effective equation of state parameter. The models with −1 < w < −1/3 are referred to as quintessence models, while the models with w < −1 are referred to as phantom models. The quintessence models can be divided in two broad classes: tracking quintessence, in which the evolution of the scalar field is slow, and thawing quintessence, in which the evolution is fast compared to the Hubble expansion [21][22][23][24].
In tracking models the scalar field exhibits tracking solutions in which the energy density of the scalar field scales as the dominant component at the time; therefore the DE is subdominant but closely tracks first the radiation and then matter for most of the cosmic evolution. At some point in the matter domination epoch the scalar field becomes dominant, which results in its effective negative pressure and accelerated expansion [25,26]. The simplest example of such a model is provided by a scalar field with an inverse-powerlaw potential energy density V φ ∝ φ −α , α > 0 [27], the so-called Ratra-Peebles model.
The scalar field models predict a different background expansion history and a growth rate compared to the cosmological constant model ones. Thus the scalar field model can be distinguished from the CDM model through high precision measurements of distances and growth rates over a wide redshift range [28][29][30][31][32][33][34][35][36][37][38].
In this paper we study generic predictions of slowly rolling scalar field models by taking the Ratra-Peebles model as a representative example. We present a self-consistent and effective way of solving the joint equations for the background expansion and the growth rate. We use a compilation of recent growth rate and baryon acoustic oscillation (BAO) peak measurements to put constraints on the parameter α describing the steepness of the scalar field's potential.
This paper is organized as follows. In Sect. 2 we investigate in detail the dynamics and the energy of the φCDM models. In Sect. 3 we study the influence of the φCDM models on the growth factor of matter density perturbations. In Sect. 4 a comparison is presented of the obtained theoretical results with observational data. We discuss our results and conclude in Sect. 5. We use the natural units with c =h = 1 throughout this paper.

Background equations
We assume the presence of a self-interacting scalar field φ minimally coupled to gravity on cosmological scales. The action of this scalar field is given by where M pl = G −1/2 is the Planck mass, with G being the Newtonian gravitational constant; V (φ) is the field's potential. Note that in this presentation the scalar field φ is dimensionless, and the potential V (φ) has the M 2 pl dimension. Following [27] we will assume that the self-interacting potential has a power-law functional form: where α > 0 is a model parameter that determines the steepness of the scalar field potential. Compliance with current observational data requires α ≤ 0.7 [39][40][41]. The larger value of α induces the stronger time dependence of the equation of state parameter w φ , while α=0 corresponds to the CDM case. Another model parameter κ > 0 is a positive dimensionless constant which is related to α (see the appendix and Ref. [42] for its dependence on α). We assume the flat and isotropic Universe that is described by the standard Friedmann-Lemaître-Robertson-Walker homogeneous cosmological spacetime model (FLRW) ds 2 = −dt 2 + a(t) 2 dx 2 , and we normalize the scale factor to be equal to 1 at present time, a today = a 0 = 1, i.e. a = 1/(1+z), where z is the redshift.
Using the action for the scalar field, Eq. (1), we obtain the Klein-Gordon equation (equation of motion) for the scalar fieldφ where an over-dot represents the derivative with the respect of physical time t; H (a) = H 0 E(a) =ȧ/a is the Hubble parameter and H 0 is its value today. The flatness of the Universe requires that the total energy density of the Universe is equal to the critical energy density, i.e. ρ tot = ρ cr = 3H 2 0 M 2 pl /(8π). We also introduce the energy density parameters for each component as i = ρ i /ρ cr (where the index i denotes the individual components, such as radiation, matter or the scalar field).
The energy density and pressure of the scalar field are given by The corresponding equation of state is given by where r 0 and m0 are the radiation and matter (including all non-relativistic components, except neutrinos, which were relativistic at the early stages) density parameters today, while ν is the total neutrino energy density which scales as ∝ a −4 before neutrinos becoming non-relativistic, and thereafter evolves as a −3 . The scalar field energy density parameter is given by To ensure the flatness of the Universe, we require that where ν0 and φ0 are the current energy density parameters for neutrinos and the scalar field, respectively. Since in the standard cosmological scenario the neutrino density is believed to be negligible compared to the matter and DE densities at low redshifts, we will ignore this component in our computations from now on (as well we neglect the radiation contribution to today's energy density).

Initial conditions
We integrate the set of equations Eqs. (3) and (6) numerically, starting from a very early moment a in = 5 × 10 −5 to the present time a 0 = 1. For the scalar field we assume the following initial conditions: where a prime denotes differentiation with respect to the scale factor a. We also used a(t) ∝ t 1/2 as consistent with a radiation dominated epoch. These initial conditions were derived from Eq. (3) (for details see Appendix A). We fix the values of the parameters, m0 = 0.315, φ0 = 0.685, h = 0.673, to the best-fit values obtained by Planck collaboration [43].   The evolution of E(a) for different values the α parameters is shown on Fig. 3. As we can expect the expansion of the Universe occurs more rapidly with increasing value of the α parameter, the CDM limit corresponding to the slowest rate of the expansion. The value of the α parameters affects also the redshift of the equality between matter and scalar field energy densities (see Fig. 4); with larger values of α the scalar field domination begins earlier and vice versa.

Growth factor of matter density perturbations in dark energy models
We use the linear perturbation equations for matter overdensities [47,48] to describe the evolution of small overdensities in a homogeneous expanding Universe, where δ ≡ δρ m /ρ m , with ρ m and δρ m the density and overdensity of the matter component, respectively.
We define by D(a) = δ(a) δ(a i ) the linear growth rate, which shows how much the perturbations have grown since initial moment a in . We normalize the growth rate so that D(a in ) = 1. The fractional matter density f 1 (a) ≡ m (a) as a function The γ parameter is also referred to as the growth index [50], and it depends on both the model of DE and the theory of gravity. In general relativity (GR) the time dependence of the γ index can be fit by [50] For the CDM model (with w = −1), the growth index is γ = 0.55 [50,51]. The φCDM model has been tested through the growth rate in Ref. [52]. In more complex coupled dark energy models, the growth rate has been studied in Refs.  [53][54][55]. The measured value of γ in conjunction with tight constraints on the other cosmological parameter can be used to test the validity of GR; see Refs. [56,57] for recent studies to use the linear growth rate data to determine the deviation of the theory of gravity on extragalactic scales from the standard GR.

The results of computations of the growth factor of matter density perturbations in φCDM dark energy model
We present the solutions of the growth Eq. (10) in RP models on Fig. 5.
We have checked that the power-law approximation Eq. (11) works well for the scalar field DE. The effective value of the growth index γ depends on α and is slightly higher than the CDM limit of 0.55 (Fig. 6).  [59]. On the bottom panel are shown the constraints, obtained after adding BAO measurements and CMB distance prior as in [60] for the BAO/CMB distance prior

Comparison with observations
The φCDM models generically predict a faster expansion rate and a slower rate of growth at low redshifts. Tight measurements of the expansion rate, distance-redshift relationship and the growth rate at multiple redshift ranges can be used to simultaneously constrain the background dynamics and the growth of structure and discriminate between φCDM and CDM models. For the rest of this section we will concentrate specifically on the discriminative power of the growth rate and BAO measurements from galaxy surveys. For simplicity we will assume that the spatial curvature is known precisely and measurements from [59] and obtain a posterior likelihood function of the parameters α and m . To do this we apply the same method as [59]; we numerically solve Eq. (10) for series of α and m values and compute a χ 2 value, where f m is the measured value of the growth rate, f th a theoretically computed value, and σ 2 f one standard deviation error of the measurement. Assuming that the likelihood is Gaussian we have The 1σ and 2σ confidence contours resulting from this likelihood are presented on the top panel of Fig. 5. The likelihood contours in the αm plane obtained from the growth rate data alone are highly degenerate. If we fix α = 0 we get m = 0.278 ± 0.03, which is within 1σ of the best-fit value obtained by the Planck collaboration [43]. Values of m < 0.2 are ruled out at more than 2σ confidence level, but large values of m are still allowed as long as α is large.
To break the degeneracy between the m and α parameters we now add a compilation of low-redshift BAO measurements from [60]. We follow the same approach as [60]; we compute the angular distance and a distance scale at a series of redshifts and construct a combination η(z) ≡ d A (z bao )/D V (z bao ) where H (z) is the Hubble parameter and H 0 is a Hubble constant. Assuming Gaussianity of the error bars we again compute the χ 2 , and a likelihood function where X = η th − η m and C is the covariance matrix of the measurements. To marginalize over the parameter H 0 in L bao we take a Gaussian prior of H 0 = 74.3 ± 2.1 from [61]. We assume that L f and L bao are independent and the combined likelihood is simply a product of the two. The results are presented on the bottom panel of Fig. 5. The addition of BAO measurements breaks the degeneracy in the growth rate data.
m is now constrained to be within 0.26 < m < 0.34 at 1σ confidence level. For the α parameter we get 0 ≤ α ≤ 1.3 at 1σ confidence level.

Discussion and conclusions
We explored observable predictions of the scalar field DE model. We showed that the model differs from CDM in a number of ways that are generic and do not depend on the specific values of model parameters. For example, in scalar field models the expansion rate of the Universe is always faster and the DE dominated epoch sets in earlier than in CDM model when other cosmological parameters are kept fixed. The two models also differ in their predictions for the growth rate, where the scalar field model generically predicts a slower growth rate than CDM.
We used a compilation of BAO, growth rate, and the distance prior from the CMB to constrain the model parameters of the scalar field model. We find that if only the growth rate data is used there is a strong degeneracy between m and α, where higher values of α are allowed as long as the m parameter is large as well. When combining these constraints with the constraints coming from a distance-redshift relationship (BAO data and the distance prior from CMB) the degeneracy is broken and we get m = 0.30 ± 0.04 and α < 1.30 with a best-fit value of α = 0.00.
where ρ ≡ ρ φ and we assume that it is the energy density of a single component that was dominant at t < t in the Universe. Assuming ρ(t) = ρ (t/t ) β we get β = −2α/(α + 2). On the other hand, assuming that the component energy density is ρ at a = a , and accounting for the dominance of the component at this epoch, we have where n = 1/2 is for radiation and n = 2/3 for the matter dominant epochs. Expressing 1/t 2 through Eq.
Since Eq. (24) must be valid for any t , we imply the freedom of our choice and use t = M −1 pl . Finally we have for n = 1/2 and n = 2/3, respectively: