Constraints on Dynamical Dark Energy with Precision Cosmological Data

We present a Dark Energy (DE) model based on a scalar field with an inverse power law potential (IPL) ( ) 4 n n V M φ φ + − = . We consider three different models 1 2, 3 4 n n = = and 1 n = and we vary the value of M and the initial amount of energy density φ Ω at the scale factor c a . We obtain a time dependent equation of state (EoS) w p φ φ φ ρ = , with 1 3 wφ = at early times for a scale factor c a a < with a steep transition to 1 wφ = at 5 1.79 10 , c a − = × 7 8 9.3 10 ,5.48 10 − − × × lasting a long period of time and a subsequent descent to 1 wφ = − for 3 ~ 10 a − to finally grow to 0.906 wφ = − , 0.932 wφ = − , 0.924 wφ = − for 1 2, 3 4 n n = = , and 1 n = respectively. The values of M and ( ) c a φ Ω are ( ) eV 4.63,127.31,2465.46 M = and ( ) 0.038,0.148, c a φ Ω = 0.227 for 1 2, 3 4 n n = = and 1 n = respectively. We show the differences in the evolution of H, the CMB and Matter power spectra, and the redshift space distortion (RSD) 8 fσ parameter. Precision cosmological data allow us to test the dynamics of Dark Energy and we obtain in all three cases a reduction of 2 20% BAO χ  compared to ΛCDM with and an equivalent fit for CMB and SNIa data.


Introduction
It has passed almost twenty years since the accelerated expansion of the universe was first observed by distance measurements using type Ia supernovae (SNeIa) [1] [2].Since that time, the acceleration of the universe has been firmly established by a series of telescopes and satellites which gather information from other cosmological probes such as the Cosmic Microwave Background Radiation (CMB) [3], Baryon Acoustic Oscillations (BAO) and Large Scale Structure (LSS) surveys [4] [5] [6].In the standard paradigm of general relativity, the physical agent causing the late-time acceleration of the universe is referred to as the Dark Energy (DE).To elucidate the nature of the DE is one of the main goals of In the standard cosmological model (ΛCDM), the energy density of the universe at present day is made up of 69% of DE described by a cosmological constant Λ, 26% of dark matter (DM) while only 5% corresponds to the Standard Model (SM) particles consisting mainly of photons, neutrinos, and ordinary matter (protons, neutrons, and electrons).However, although the ΛCDM model has proved to agree very well with the observations [3], there is no understanding of the physical mechanism that determines the origin and the magnitude of the cosmological constant Λ and hence of why and when the universe accelerates [10].This lack of understanding is commonly expressed in terms of the "fine-tuning" and the "coincidence" problems; the former asks why the observations set the value of the DE density ρ Λ to almost 10 120 orders of magnitude below conservative estimations based on quantum field theory, while the "coincidence" problem inquires why the energy densities of the DE and the matter are of the same order of magnitude precisely at present time.This scenario has enforced the quest of other mechanisms to explain the nature of DE.
Alternative to Λ, scalar fields φ have been extensively explored as a source of DE [11]- [16] because of their intrinsic link with the Standard Model of particle physics and its extensions, opening up the possibility to explain the nature of DE from first principles.In this paper we study the constraints set by recent cosmological observations on a scalar field model with inverse power law (IPL) until recent times, our model considers the possibility that this component may play also an important role in the early universe, leading to cosmological imprints that can be tested against current observational information.However, here we are more interested in constraining IPL potentials with recent precision cosmological data than in a theoretical derivation of the IPL potential.
The outline of the paper is the following.In Section 2 we discuss the basic picture of the model and present the dynamical equations.The constraints set by the data are presented in Section 3. We study the cosmological consequences of the model in Section 4 and give our conclusions in Section 5.

Dynamical Dark Energy
Scalar field φ theories have been proposed as possible sources to describe DE and a wide range of models have been studied in recent years [11]- [16] [24].
Particularly, IPL potentials proposed by RP [25] [26] [27] have been widely investigated giving an equivalent fit as ΛCDM [28] [29] [30].The evolution of the energy density ( ) and equation of state (EoS) w p φ φ φ ρ = depends on the parameters n, M (with mass dimensions) and the initial conditions of φ .All these three parameters have an important impact in determining the size and slope of w φ at present time and must be adjusted by the cosmological observations [15] [16] [28] [29] [30].These quantities are free parameters to be adjusted by the cosmological observations or the choice of model.
The evolution of φ in a homogeneous flat universe described by the Friedmann-Lemaître-Robertson-Walker metric is completely determined by the Klein-Gordon equation where the dots stand for cosmic time derivatives, and the Hubble expansion rate: ( ) with ( ) the total energy density, mo ρ , ro ρ the present day energy densities of matter and radiation and the redshift z is given by ( ) The energy density and pressure for the scalar field φ are ( ) ( ) with an equation of state (EOS) and a mass given by ( ) Since the potential V in Equation ( 1) is an inverse power of φ , the evolution of φ is an increasing function of time, i.e. of the expansion of the universe, while the mass ( ) m φ is a decreasing function.
We define the scale factor c a and mass parameter M as the scale when the scalar field takes the value of M (i.e.
The parameters c a , M and the DE energy density parameter a are related by evolution is determined by Equations (( 2) and ( 3)) with the IPL potential in Equation ( 1).
The homogeneous background approximation must be refined by considering also the perturbations of the different fluids.We stay in the linear regime where the energy density and other quantities can be decomposed into a homogeneous part (commonly denoted with a bar) and a small position-dependent perturbation.
We solve the perturbed equations in the syncrhonous gauge defined by the line element [31]: where d dt a τ ≡ denotes the conformal time.The perturbations in the energy density and pressure of the scalar field are given by: where the primes stand for conformal time derivatives.The evolution of δφ in Fourier space is determined by: Here a a ′ ≡  is the conformal expansion rate and The scalar field enters the perturbation equations of the other fluids [31] via  and through extra source terms proportional to δφ and δφ′ .For example, the CDM overdensities c δ evolve according to: ( ) where the sum runs over all the fluids with sound speed 2 , s i i i c P δ δρ = .

Constraints
We explore the parameter space using the CosmoMC [32] and CAMB [33] 2χ ) of some selected parameters from measurements on the CMB temperature anisotropy spectrum [3], BAO [4], and SNeIa [34] data.As a reference, we include the results for the ΛCDM model.
Notice that the IPL models improve the BAO measurements by 20% compared to ΛCDM and have an equivalent fit for CMB and SNIa.This shows that a dynamical DE has a better fit than ΛCDM, where the BAO measurements allow us to constrain the dynamics of the DE.
We see that the transition from the radiation-like state to the scalar-field state takes place well deep in the radiation era for all the cases.However, the transition occurs later as we go to smaller values of n, while the scalar field density at c a ( c φ Ω ) and the scale of energy M increase as we go to larger n.The evolution of the matter, radiation and DE densities is shown in Figure 1.
Since our model contains at early times ( c a a < ) relativistic particles, its energy density evolves as  since the massless fields of the dark group evolve as radiation.Once the transition Figure 1.We show the evolution of energy densities ρ as a function of the scale factor a for radiation (dotted orange), dark matter (dashed black), Λ (dashed red) and IPL scalar field φ (solid blue)) and the vertical dotted line marks the transition epoch  to the scalar field occurs, the EOS leaps abruptly to 1 and stays at this value for some time, then drops to −1 around the decoupling epoch ( * 1090 z ≈ for the three models) and finally grows to 1 DEo w > − at recent times.The inner panel zooms into the late-time behaviour.
It is interesting to note that in [30] an IPL potential with 1 2 n = was also studied but with different initial conditions for φ .In [30] the evolution of the EoS close to present time is a decreasing function from 0.8 w = − to 0.87 o w = − .
They assumed that the potential V has reached the tracking regime [14].As we can see, different initial conditions have a distinctive evolution of the EoS since in our case we have a growing EoS from time (see Table 1).Clearly the choice of initial conditions of φ is also important.

Distances
The different late-time dynamics followed by the EOS leaves distinctive imprints on some quantities probed by the cosmological observables.The immediate consequence of such difference is the change in the amount of DE which modifies the size and evolution of the expansion rate H affecting the cosmological distances measured by the SNeIa, BAO and CMB observations.In this regard, we recall that SNeIa flux measurements provide an estimation of the luminosity distance given by ( ) ( ) ( ) , while BAO measurements are sensitive to the ratio ( ) ( ) 1), where drag r is the coomoving sound horizon at the drag epoch, ( ) ( ) ( ) ( ) ( ) ( ) ( ) is the angular diameter distance of the reconstructed cluster of galaxies.Figure 3 shows the comoving expansion rate ( ) ( ) for the three IPL models and ΛCDM.We see that ( ) ( )

( )
* A D z -differences less than 0.08%-and since the amount of matter is roughly the same-differences less than 0.2%-, the amount of φ ρ at present time must be reduced from ρ Λ , giving a smaller

CMB and Matter Power Spectrum
At perturbation level, the IPL model also leaves important imprints on the CMB power spectrum and the evolution of matter perturbations.When we run the models with the same parameters, a change in cosmological distances is reflected in a shift of the position of the peaks of the CMB power spectrum.Moreover, although the amount of extra early radiation vanishes once the transition occurs and therefore it does not have any direct influence on the physical processes from that time onwards, a larger value of primordial helium fraction P Y produced by this extra radiation (c.f.Table 1) leads to a more damped tail of the CMB spectrum [35] [36] at high multipoles 1000 l ≥ . On the other hand, the low-l region is sensitive to perturbation dynamics via the late-time ISW effect.In general, models with DE perturbations with a sound speed 2 0 s c > enhance the power [37] [38] [39].When fitting to the data, some of these effects are compensated or even mimicked by a proper combination of the cosmological  1 n = have a slightly better fit to the data, we note that the amount of early extra radiation introduced by the model with 1 n = is in unavoidable tension with Big Bang Nucleosynthesis constraints on primordial helium [40] [41] and deuterium [42] [43].
As far as the evolution of matter perturbations is concerned, Figure 5 shows the resulting matter power spectra.In general, the spectrum is suppressed in the IPL models, where the differences with respect to ΛCDM lie within the 6% in the linear regime .Future precision cosmological measurements will allow us to further constraint the dynamical DE models where RSD will play an important role as well as better distance measurements.

Conclusion
Here we consider three different IPL potentials

1 n
= and we vary the value of M and the initial amount of energy density φ Ω at the scale factor c a .We obtain a time dependent equation of state (EoS) period of time and a subsequent descent to 1 w φ = − for

1 n
= respectively.We show the differences modern cosmology.The amount and precision of observational data collected in the last decade have grown up impressively.It is now possible to test and discriminate different DE models and projects scheduled to start operating in the near future such as DESI [7], LSST [8], and Euclid [9] will provide conclusive evidence on the subject.

N
fit to cosmological observations.The IPL potentials may arise from the non-perturbative dynamics of a non-abelian gauge group massless fields[17]-[23].Unlike ordinary quintessence where the scalar field representing the DE is subdominant

=
) and at this time we have

(
particles in thermal equilibrium as long as they are relativistic we have m T  with T the temperature with an energy density given by ( ) the degrees of freedom of φ , and an EoS 1 3 w φ = .Notice that at c a we have ( ) e. ~m M T .Therefore, for scales c a a < we have m T < and the evolution of φ has an EoS 1 3 w φ = while for c a a > the scalar fields are no longer relativistic m T > and its the scalar field quickly dilutes just after the phase transition at c a , since the EOS leaps to 1 so that


at early times but not in IPL models.

Figure 2 .
Figure 2. We show the evolution of the EoS w φ as afunction of a. Notice the steep transition from 1 3 w φ = to 1 w φ = for

Figure 3
Figure 3. Evolution of

Figure 4 .
Figure 4. CMB Power Spectrum for IPL models and ΛCDM.Lower panel shows the residuals with respect to ΛCDM.Grey dots correspond to Planck 2015 measurements [3].
best fit point as input parameters, we can make predictions on quantities that characterize structure formation such as the redshift space distortion (RSD) parameter 8 f σ shown in Figure6as a function of z.Like the matter power spectrum, the IPL models suppress 8 f σ by nearly a constant fraction in the range 0.6 1.2 z < < We have seen that they have a better cosmological fit than the standard ΛCDM model, reducing BAO 2 χ by 20% and leaving CMB an SNIa with an equivalent fit.With more precise measurements coming in the near future we could be at the stage to determine the dynamics of Dark Energy.In particular BAO and RSD data will be vital to achieve this goal.

Figure 5 .
Figure 5. Matter Power Spectrum for IPL models and ΛCDM.Lower panel shows the relative differences.

8f
σ as a function of z for IPL models and ΛCDM.Lower panel shows the relative differences.

Table 1 .
Best fit point for each dark energy model.
The same trend is observed in o H and the DE density parameter at present.Nevertheless, none of these two is larger than in ΛCDM.o H is expressed in km⋅s −1 Mpc −1 .