It has passed almost twenty years since the accelerated expansion of the universe was first observed by distance measurements using type Ia supernovae (SNeIa)   . 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)  , Baryon Acoustic Oscillations (BAO) and Large Scale Structure (LSS) surveys    . 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 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  , LSST  , and Euclid  will provide conclusive evidence on the subject.
In the standard cosmological model (LCDM), the energy density of the universe at present day is made up of 69% of DE described by a cosmological constant L, 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 LCDM model has proved to agree very well with the observations  , there is no understanding of the physical mechanism that determines the origin and the magnitude of the cosmological constant L and hence of why and when the universe accelerates  . 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 10120 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 L, scalar fields have been extensively explored as a source of DE  -  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) potential giving a good fit to cosmological observations. The IPL potentials may arise from the non-perturbative dynamics of a non-abelian gauge group with colors and massless fields  -  . Unlike ordinary quintessence where the scalar field representing the DE is subdominant 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.
2. 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  -   . Particularly, IPL potentials
proposed by RP    have been widely investigated giving an equivalent fit as LCDM    . The evolution of the energy density , pressure and equation of state (EoS) 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 at present time and must be adjusted by the cosmological observations      . 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, , the present day energy densities of matter and radiation and the redshift z is given by with . 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 is a decreasing function.
We define the scale factor and mass parameter M as the scale when the scalar field takes the value of M (i.e. ) and at this time we have
The parameters , M and the DE energy density parameter at are related by
where the initial conditions at are , , , and .
For particles in thermal equilibrium as long as they are relativistic we have with T the temperature with an energy density given by , with the degrees of freedom of , and an EoS . Notice that at we have , i.e. . Therefore, for scales we have and the evolution of has an EoS while for the scalar fields are no longer relativistic and its 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  :
where 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 is the conformal expansion rate and . The scalar field enters the perturbation equations of the other fluids  via and through extra source terms proportional to and . For example, the CDM overdensities evolve according to:
where the sum runs over all the fluids with sound speed .
We explore the parameter space using the CosmoMC  and CAMB  codes properly modified to account for the field background and linear perturbation
dynamics described in section 2. In our analysis, we take and find
the best fit point for each case varying the physical densities of baryons ( ) and cold dark matter ( ), the optical depth ( ), the spectral index ( ), the amplitude of scalar perturbations ( ), and the density parameter of the scalar field at that time ( ) and M, where and stand for the densities of radiation (photons and three massless neutrino species) and matter (baryons and CDM), respectively; is given by the solution of the system of Equations ((2) and (3)) at present time. Table 1 displays the best-fit values with the corresponding goodness of fit ( ) of some selected parameters from measurements on the CMB temperature anisotropy spectrum  , BAO  , and SNeIa  data. As a reference, we include the results for the LCDM model. Notice that the IPL models improve the BAO measurements by 20% compared to LCDM and have an equivalent fit for CMB and SNIa. This shows that a dynamical DE has a better fit than LCDM, 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 ( ) and the scale of energy M increase as we go to larger n. The same trend is observed in and the DE density parameter at present. Nevertheless, none of these two is larger than in LCDM.
Table 1. Best fit point for each dark energy model. is expressed in km×s−1Mpc−1.
The evolution of the matter, radiation and DE densities is shown in Figure 1. Since our model contains at early times ( ) relativistic particles, its energy density evolves as . However, the scalar field quickly dilutes just after the phase transition at , since the EOS leaps to 1 so that (c.f. Figure 2). Finally, the scalar field remains subdominant until late times, where slightly but noticeably different from a cosmological constant since . Notice that at early times (small a), exposing the naturalness and coincidence problems of the LCDM model. The evolution of the scalar field EoS is depicted in Figure 2. Initially we have (dotted line) 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), L (dashed red) and IPL scalar field (solid blue)) and the vertical dotted line marks the transition epoch . Notice that at early times but not in IPL models.
Figure 2. We show the evolution of the EoS as afunction of a. Notice the steep transition from to for and a leap to for to finally grow to at present time for all three cases .
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 ( for the three models) and finally grows to at recent times. The inner panel zooms into the late-time behaviour.
It is interesting to note that in  an IPL potential with was also studied but with different initial conditions for . In  the evolution of the EoS close to present time is a decreasing function from to . They assumed that the potential V has reached the tracking regime  . As we can see, different initial conditions have a distinctive evolution of the EoS since in our case we have a growing EoS from to at present time (see Table 1). Clearly the choice of initial conditions of is also important.
4. Cosmological Implications
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 (c.f. Table 1), where is the coomoving sound horizon at the drag epoch, and 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 is larger in IPL in the range sensitive to BAO and SNIa measurements. Moreover, since the accuracy of the angular size of the sound horizon at recombination obtained from CMB measurements forces IPL and LCDM to have the same angular distance ―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 in IPL than in LCDM.
4.2. 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 produced by this extra radiation (c.f. Table 1) leads to a more damped tail of the CMB spectrum   at high multipoles . 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 enhance the power    . When fitting to the data, some of these effects are compensated or even mimicked by a proper combination of the cosmological
Figure 3. Evolution of (in km×s−1×Mpc−1) for IPL models and LCDM. Lower panel shows the relative differences.
Figure 4. CMB Power Spectrum for IPL models and LCDM. Lower panel shows the residuals with respect to LCDM. Grey dots correspond to Planck 2015 measurements  .
parameters. Figure 4 shows the fit of the IPL and LCDM models to the Planck 2015 CMB temperature anisotropy spectrum. Although the models with and have a slightly better fit to the data, we note that the amount of early extra radiation introduced by the model with is in unavoidable tension with Big Bang Nucleosynthesis constraints on primordial helium   and deuterium   .
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 LCDM lie within the 6% in the linear regime . Taking the best fit point as input parameters, we can make predictions on quantities that characterize structure formation such as the redshift space distortion (RSD) parameter shown in Figure 6 as a function of z. Like the matter power spectrum, the IPL models suppress by nearly a constant fraction in the range . 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.
Here we consider three different IPL potentials models with a power and . We have seen that they have a better cosmological fit than the standard LCDM model, reducing BAO 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. Matter Power Spectrum for IPL models and LCDM. Lower panel shows the relative differences.
Figure 6. as a function of z for IPL models and LCDM. Lower panel shows the relative differences.
We acknowledge financial support from UNAM DGAPA IN103518 and Conacyt Fronteras 281.