Can dark energy explain a high growth index?
Abstract
A promising way to test the physics of the accelerated expansion of the Universe is by studying the growth rate of matter fluctuations, which can be parameterized by the matter energy density parameter to the power , the so-called growth index. It is well-known that the CDM cosmology predicts . However, using observational data, Ref. [1] measured a much higher , excluding the CDM value within . In this work, we analyze whether Dark Energy (DE) with the Equation of State (EoS) parameter described by the CPL parametrization can significantly modify with respect to the CDM one. Besides the usual Smooth DE (SDE) scenario, where DE perturbations are neglected on small scales, we also consider the case of Clustering Dark Energy (CDE), which has more potential to impact the growth of matter perturbations. In order to minimally constrain the background evolution and assess the largest meaningful distribution, we use data from Cosmic Chronometers, ), data points. In this context, we found that both SDE and CDE models described by the CPL parametrization can not provide a significant number of samples compatible with the value determined in Ref. [1]. Therefore, explaining the measured value of is a challenge for DE models. Moreover, we present new fitting functions for , which are more accurate and general than the one proposed in Ref. [2] for SDE, and, for the first time, fitting functions for CDE models.
I Introduction
The accelerated expansion of the universe is still a big question in Cosmology. It can be explained either by a modified theory of gravity, or a unknown component of the universe, labeled as Dark Energy (DE) [3]. In the pursue of more accurate data to answer this question, some tensions arose. The most significant one is known as the Hubble tension, which reflects the difference between measurements of the present expansion rate, , obtained locally by using the distance ladder methods [4] and globally, assuming the CDM model, from Cosmic Microwave Background (CMB) data [5]. There is also the tension (where , is the amplitude of matter fluctuations at and is the matter density parameter now), which is related to the difference between the values of these parameters inferred from CMB [5] and measurements of galaxy clustering and weak gravitational lensing, as discussed in [6].
The tension is directly related to the growth of cosmic perturbations. A particular simple and promising way to test how the physics of the accelerated expansion impacts the evolution of matter perturbations is by analyzing the growth rate of matter perturbations [7]
(1) |
which can be parametrized as [2]
(2) |
where is a constant that depends weakly on the Equation of State (EoS) parameter of Smooth DE (SDE) models, , and is the matter energy density parameter. For the CDM cosmology, Refs. [8, 2] found . However, using observational data, Ref. [1] found a much higher (implying a suppression of growth), excluding the CDM value by . In the same work, it was also shown that a higher values reduces the tension from to .
This high value of naturally raises the question of what theoretical mechanisms are capable of producing it. In the case of SDE, it was shown that can be accurately described by the fitting function [2],
(3) |
where . This indicates that, for , we would need to get , in contradiction with the parametrization condition. For , we would need , which would generate a large DE density around . This simple extrapolation exercise indicates that it should be very difficult for SDE models to produce a high growth index. As we will show, with very loose constraints on the background evolution based on Cosmic Chronometers (CC) data, SDE with CPL parametrization has an almost negligible probability of producing .
The usual assumption of SDE, in which one usually neglects DE perturbations on small scales, is based on Quintessence models [9, 10, 11]. In this case, DE is described by a minimally coupled canonical scalar field. The linear perturbations of this field propagate with sound speed , not allowing its perturbations to grow significantly on small scales. This approximation is well justified even in the nonlinear regime [12]. The simplest generalization of this scenario can be done by describing DE as a non-canonical minimally coupled scalar field, called k-essence [13, 14]. In this case, can be chosen and DE pertubations can grow at small scales, see [15] for a more detailed discussion. In this work, we will consider the limiting case of , in which DE pertubations have the maximal potential to grow and impact the evolution of matter perturbations. We refer to this scenario as Clustering DE (CDE), which growth index has already been studied in [16, 17, 18, 19].
In this paper, we explore how SDE and CDE can change the growth index. We confirm the expectation based on the fitting formula (3) that SDE model can not provide . We also find that CDE models can not raise the values of significantly, therefore concluding that DE models are quite challenged by observed value of .
As we will show, for the case of CDE, this result has a simple explanation. When , DE perturbations are have the tendency of being correlated with matter perturbations () and anti-correlated for (). In other to raise , DE perturbations must be anti-correlated with , causing a decrease in the gravitational potential and consequently slowing down the growth rate (higher ). However, in the case , the DE energy density decreases rapidly at high-, thus the overall impact is very limited. On the other hand, the case can easily enhance the matter growth, thus providing a significantly lower values, as already shown in [16].
Based on the samples of that we have computed, we where able to construct a new parametrization , which is more accurate than (3) and valid for larger parameter space. Moreover, for the first time, we present a parametrization for the case of CDE.
The plan of this paper is as follows. In Sect. II, we define the background cosmology, the data, the parameters priors and statistics used to constrain the background evolution. In Sect. III, we present the equations for the evolution of matter and DE pertubations. The results of the sampling, analysis and suitable fittings are shown in Sect. IV and the conclusions are organized in Sect. V.
II Background cosmology
In this work, we assume General Relativity and a flat universe described Friedmann-Lemaître-Robertson-Walker metric, in which the line element is represented by
(4) |
where is the scale factor. As so, the square of the Hubble function can be written as
(5) |
where is the matter (dark matter plus baryons) density parameter now and is the DE energy parameter, which depends on the EoS assumed.
We consider that DE EoS is given by Chevallier-Polarski-Linder (CPL) [20, 21] parametrization:
(6) |
Thus, the DE density parameter is given by
(7) |
Besides the general case of free and , we will also analyze the limit ( and ) and constant EoS case (free and ). We refer to these models as CPL, CDM and CDM, respectively. In the general case, we have four free parameters: , , and . Next, we discuss how to minimally constrain these parameters using data.
Cosmic Chronometers Data
In order to minimally constraint the background evolution and the model parameters, we make use 32 of the most recent CC data available, compiled by [22]. We use the Python library Emcee [23] as a Monte Carlo Markov Chain (MCMC) sampler and Python library GetDist [24] to analyse and plot the posteriors distributions.
There is a well-know discussion about the systematic errors in CC data, [25, 26, 27]. In these works, systematics of 15 of the 32 data points have been analyzed. Here, we split the data set in two groups (with and without systematics), and given by
(8) |
where is associated with the 17 data-points without systematics, provided by [28, 29, 30, 31, 32, 33], which reads
(9) |
where and is the corresponding uncertainty of data points. The other part of the associated with data points with systematics is given by
(10) |
where is a vector with components associated the 15 data points provided by Refs. [25, 26, 27], and is the corresponding covariance matrix for BC3 model, available at https://gitlab.com/mmoresco/CCcovariance.
We sample the posterior distribution of the vector of parameters , given by
(11) |
where
(12) |
is the likelihood of the data given and indicates the priors assumed for the parameters, listed in Tab. 1. We also implement an additional flat prior to exclude cosmologies in which , where is the redshift. This condition is similar to the assumption described by [34]. Allowing for DE models with non-negligible energy density at high-, besides being highly disfavored by data [35], invalidates the Einstein-de-Sitter initial conditions used to solve the evolution of matter and DE perturbations, which will be discussed in detail in the next section.
Parameter | Prior |
---|---|
The results for CDM, CDM and CPL models are presented in figures 1, 2 and 3, respectively. Regarding the background parameters, we can see in Fig. 1 that CC data can constraint the two free parameters, and , in the CDM model. For CDM this is not case, mainly because is being limited by the prior assumed, but note that the posterior also shows some preference for . Fig. 3 also shows that CPL parameters, and , can not be constrained by CC data. The posteriors of and are dominated by the flat priors shown in table 1 together with the condition . The mean values for the marginalized 1D distributions are given in Tab. 2. The number of samples in our chains is , where is the largest auto-correlation time of the parameters. This is well beyond the suggested described in [23] as convergence criteria.
We recall that the main purpose of this analysis is to provide a wide, but still meaningful, distribution of background parameters that can be used to compute the linear growth of matter perturbations and . As we will see, despite this large variation of the EoS parameters, can be determined with much smaller variation. In the next section we describe the perturbative equations for SDE and CDE models.
III Growth of Cosmological Perturbations
In a universe with SDE, the linear evolution of the matter density contrast, , is solely affected by the background expansion, as described by the equation
(13) |
where the prime denotes a derivative with respect to the scale factor. Starting the integration at a high redshift, , initial values can be computed using the analytical solution for a matter-dominated Einstein-de-Sitter (EdS) universe , where we assumed .
In the case of CDE, is also affected by DE perturbations, , whose effect is maximized for DE models with negligible sound speed. In this work, we consider the extreme case based on the fluid description of DE and phenomenologically allowing for EoS which can be non-phantom (), phantom () or transit between these regimes. In this context, matter and DE perturbations obey the following system of equations [15]:
(14) |
where and is the divergence of the peculiar velocity of the DE-matter fluid. Note that the equation for the peculiar velocity is the same for matter and CDE because we are considering . To compute the initial conditions for DE perturbations, we use the solution valid for constant in matter-dominated era [36, 37, 38, 39]
(15) |
Although solution (15) only gives a general qualitative behavior at low- [39], it can help us to predict the influence of CDE on the growth of . When we have the tendency , while for we have . If the EoS is always phantom or non-phanton, these relations are always valid. In the case that the EoS transits between phantom and non-phanton or vice-versa, the actual relation between and might take some time to achieve the expected behavior. The general trend is clear: phantom EoS will induce negative , lowering the clustering power encoded in the last term of Eq. (14), whereas non-phantom EoS enhances growth. Therefore, the desired higher values of could be obtained in phantom CDE models.
We solve Eqs. (13) and (14), then determine the growth rate and fit a constant assuming for all the posterior distribution obtained and the measured value. In this process, it is important to ask whether the fitted values of can accurately describe . In the paper [2], where the fit (3) was introduced for the SDE model described by CPL parametrization, the accuracy of the fit was reported with respect to the growth variable . Accuracies better than for CDM with were reported, and when lowering the limit to . For , the paper reported an accuracy within for and for dynamical EoS when and . Here, we first will check how accurately a constant can reproduce . In the next section, we will propose and analyze a new fitting function for .
In this work, we report the distribution of percent residuals given by the Root Mean Square (RMS) of , where is the numerical solution and the corresponding fitted gamma value. We compute the RMS percent residuals along evenly spaced values of , and show its distribution in Fig. 4. As can be seen, a very small fraction of models have residuals greater than . The worst case occurs for CDE-CPL, reaching up to , but still for a small fraction of the realizations. In general, all models present a concentration of residuals close to . We also have checked that the largest residuals are mainly associated with low values. This demonstrates that, considering the CPL EoS, the constant parametrization for the growth rate is very reliable, even in the case of CDE. However, for more complex EoS parametrizations and for realizations including a non-negligible DE density at high-, the constant parametrization might not be adequate [16].
IV Results for
The CDM posteriors shown in Fig. 1 display a great determination of the growth index at , agreeing with the previous results from Refs. [8, 2]. However, using data from Cosmological Microwave Background, galaxy surveys and Baryon Acoustic Oscillation data, Ref. [1] found a value of , excluding the CDM value within , also showing that this value effectively solves the tension. Assuming normally distributed probabilities for the measured value and the CDM one, with mean and standard deviation given in Tab. 2, a simply quantification of the tension between these values is given by
(16) |
where we used . If we assume the usual value and neglect its uncertainty, we get .
Now let us check whether more general SDE models and their CDE counterparts can explain the observed value. Although the analysis in Ref. [1] assumes the CDM background to produce the constraints, assessing the possible values of in more general background and perturbative models will help us to identify scenarios in which can be increased and how likely this can happen.
We first analyze the correlations between the and the EoS parameters and how changes with respect to the CDM value. The clearest case is for CDM model, Fig. 2, but the following analysis also holds for CPL model. For (non-phantom EoS) we have lower in the past with respect to the case. In the SDE scenario, this behavior causes a suppression of growth, giving a higher . In the case of CDE, this correlation is inverted because induces positive , which will act as an extra source of the gravitational potential, enhancing the growth and lowering . For (phantom EoS) we have higher in the past with respect to . In SDE case, this induces a lower . Again, CDE inverts this correlation because now negative is induced, reducing the source of the gravitational potential, consequently increasing . These correlations also hold for CPL, but since many realizations shown in Fig. 3 include transitions from phantom to non-phantom EoS and vice-versa, they are not very clearly visualized. However, in terms of , these correlations become more evident in Fig. 7. Therefore we can summarize impact of DE model in the , with respect to the CDM value, as follows:
-
•
SDE, non-phantom EoS:
-
•
SDE, phantom EoS:
-
•
CDE, non-phantom EoS:
-
•
CDE, phantom EoS:
It is also interesting to check the frequency of positive or negative occurrences of for the CDE scenario. In Fig. 5, we present the distribution of at and . We can see a small preference for negative , which is associated with the allowed values of the EoS parameters as follows. For non-phantom EoS, can be large at intermediate and high-, a situation that is disfavored by data, e.g., [35]. In our analysis, this fact is mainly implemented with the prior . On the other hand, for phantom EoS, is very small at intermediate and high- and can still can induce an adequate accelerated expansion at low-. Therefore, the allowed parameter space for the EoS parameters has a larger fraction of phantom realizations. As a consequence, given the correlations between and explained earlier, we see some preference for .
In Fig. 6 we show a direct comparison between the distributions obtained. As can be seen, the SDE models are similar, regardless of the EoS parametrization. The same happens for the CDE models. Given the larger fraction of phantom realizations, SDE models have a slight preference for , whereas CDE prefers . The distributions for non- models have greater variance and asymmetry. However, none of these models provide a significant fraction of realizations around .
Considering the lower limit of the measured value as , only the SDE-CDM model has values higher than this, what happens only above percentile 96 of its distribution. If we cut this distribution with , occurs only above percentile 99. As can be seen in Fig. 6, all the other models have much smaller probabilities of producing .
We could expect that phantom CDE models are able to produce high values because, as explained earlier, the associated negative operates in this direction. However, given that the gravitational potential depends on and, in the phantom case, is small at intermediate and high-, the actual impact of these models on is very limited. Consequently, although some preference for can be seen in Fig. 6, the highest possible values virtually never reach .
Therefore, our main result is that both SDE and CDE models described by the CPL EoS parametrization have very small probabilities of providing a growth index compatible with close to , found in Ref. [1]. Bear in mind that in order to produce the distributions we only used CC data, the priors in Tab. 1 and the condition . Therefore, adding more background data in the analysis can only constrain more the parameter space and consequently the distribution.
The mean values for the marginalized 1D distributions are listed in Tab. 2. Note that, as the distribution is heavy tailed for CDE models, and reported mean value can lie outside the interval. We verified that the best fit model parameters lie within of the 1D marginalized distributions.
Parameter | CDM | CDM | CPL |
or | - - - | ||
- - - | - - - | ||
(SDE) | |||
(CDE) | - - - |
New fits for
With the posterior distributions obtained, we were able find a more accurate and general fitting function for , given by
(17) |
The coefficients for each model are listed in Tab. 3. The tested range of is . The main difference between SDE and CDE models is the sign of the coefficient in the exponential term.
Model | Cosmology | ||||
---|---|---|---|---|---|
SDE | CDM | ||||
SDE | CPL | ||||
CDE | CDM | ||||
CDE | CPL |
The quality of the fits are demonstrated in the upper panels of Fig. 7, where our fits are well contained within the posteriors contours. The lower panels show the residuals of the parametrization of Eq. (3) and ours. We computed these residuals with the same RMS definition used for Fig. 4.
Let us first analyze the results for SDE in Fig. 4. Our fit is most accurate for the CDM model (left panels), presenting a mode of residuals at , while the linear fit of Eq. (3) presents a mode around . In any case, very few realizations have with residuals . For the CPL models (right panels), our fit is noticeable better than the one of Eq. (3). For instance, the linear fit residuals, Eq. (3), peaks around while the exponential fit, Eq. (17), peaks at . The main reason for this improvement is that the linear fit does not perform well for , as can be seen in the top panels of Fig. 4.
Our fitting function for can also be used to describe the values for the CDE case, which are shown with orange lines and contours in Fig. 4. As can be seen, it also performs very well in models with clustering DE. The modes of the RMS residuals is around for both and CPL parametrizations.
As an example of the use of our fit, in Tab. 4, we present the values associated with the best-fit parameters obtained by DESI-BAO analysis [40]. As expected from the analysis already presented here, neither SDE nor CDE models are able to produce a high value, but note that CDE can provide a slightly higher because the values have a phantom behavior at high an intermediate . It is interesting to note that all these three DESI data analysis provide , thus the values for SDE or CDE are quite similar.
V Conclusions
In this paper, we have analyzed the possibilities of DE models described by the CPL parametrization of producing a high value of , compatible with obtained in Ref. [1]. As summarized in Fig. 6, both smooth DE and clustering DE have distributions with small or negligible overlap with the interval of the measurement. This result depends only on the following assumptions: (1) must be compatible with CC data, (2) and (3) the priors described in Tab. 1. The combination of these assumptions imposes very loose constraints on parameters, allowing for a vast, but meaningful, exploration of the possible values of .
We have also analyzed the correlations between , the EoS parameters and the clustering properties of DE. In particular, we can expect that phantom CDE models can raise the value of . However, since the impact of such models depend on , the actual change is very limited because decays rapidly with . Non-phantom SDE can also give higher , but these realizations are correlated with very low .
We have also proposed a new fitting function for , with overall accuracies better than the linear fit of Eq. (3) covering the interval of . For the first time, we produced a fitting function for CDE models, Eq. (17) with coefficients given in Tab. 3. These fits can be useful for a fast estimation of for both SDE and CDE models described by CPL parametrization.
If the high value found in Ref. [1] is validated by other observations and analysis, it poses a significant challenge to DE models based on minimally coupled scalar fields that can be described by the CPL EoS and a constant on small scales. Since we have considered the two limiting cases of smooth DE () and full clustering DE (), it seems unlikely that intermediate or time-varying values can produce a significantly higher . A possible alternative is to consider more general EoS parametrizations. However, as demonstrated in our study, large variations on the EoS parameters and produce much narrower distributions. Therefore, in principle, even more general EoS parametrizations should have some difficulties in explaining such high values, with the disadvantage of introducing more parameters. If this ‘ tension’ persists, we might be seeing an early evidence of modified gravity or non-standard DM.
Acknowledgements.
We thank Valerio Marra for useful discussions. IBSC thanks Universidade Federal do Rio Grande do Norte for the scientific initiation fellowship, N° 01/2023 (PIBIC-PROPESQ) project PIJ20915-2023.References
- Nguyen et al. [2023] N.-M. Nguyen, D. Huterer, and Y. Wen, Evidence for Suppression of Structure Growth in the Concordance Cosmological Model, Phys. Rev. Lett. 131, 111001 (2023), arXiv:2302.01331 [astro-ph.CO] .
- Linder [2005] E. V. Linder, Cosmic growth history and expansion history, Phys. Rev. D72, 043529 (2005), astro-ph/0507263 .
- Abdalla et al. [2022] E. Abdalla et al., Cosmology intertwined: A review of the particle physics, astrophysics, and cosmology associated with the cosmological tensions and anomalies, JHEAp 34, 49 (2022), arXiv:2203.06142 [astro-ph.CO] .
- Riess et al. [2021] A. G. Riess, W. Yuan, L. M. Macri, D. Scolnic, D. Brout, S. Casertano, D. O. Jones, Y. Murakami, L. Breuval, T. G. Brink, A. V. Filippenko, S. Hoffmann, S. W. Jha, W. D. Kenworthy, G. Anand, J. Mackenty, B. E. Stahl, and W. Zheng, A comprehensive measurement of the local value of the hubble constant with 1 km/s/mpc uncertainty from the hubble space telescope and the sh0es team, The Astrophysical Journal Letters 934, L7 (2021), arXiv:2112.04510 [astro-ph.CO] .
- Aghanim et al. [2020] N. Aghanim et al. (Planck), Planck 2018 results. vi. cosmological parameters, Astron. Astrophys. 641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO] .
- Di Valentino et al. [2021] E. Di Valentino et al., Cosmology intertwined iii: and , Astropart. Phys. 131, 102604 (2021), arXiv:2008.11285 [astro-ph.CO] .
- Amendola et al. [2013] L. Amendola et al. (Euclid Theory Working Group), Cosmology and fundamental physics with the Euclid satellite, Living Rev.Rel. 16, 6 (2013), arXiv:1206.1225 [astro-ph.CO] .
- Wang and Steinhardt [1998] L.-M. Wang and P. J. Steinhardt, Cluster abundance constraints on quintessence models, Astrophys. J. 508, 483 (1998), astro-ph/9804015 .
- Peebles and Ratra [1988] P. J. E. Peebles and B. Ratra, Cosmology with a time variable cosmological ’constant’, Astrophys. J. 325, L17 (1988).
- Ratra and Peebles [1988] B. Ratra and P. J. E. Peebles, Cosmological consequences of a rolling homogeneous scalar field, Phys. Rev. D37, 3406 (1988).
- Wetterich [1988] C. Wetterich, Cosmology and the fate of dilatation symmetry, Nucl. Phys. B302, 668 (1988).
- Batista et al. [2023] R. C. Batista, H. P. de Oliveira, and L. R. W. Abramo, Spherical collapse of non-top-hat profiles in the presence of dark energy with arbitrary sound speed, JCAP 02, 037, arXiv:2210.14769 [astro-ph.CO] .
- Armendariz-Picon et al. [1999] C. Armendariz-Picon, T. Damour, and V. F. Mukhanov, k - inflation, Phys. Lett. B 458, 209 (1999), arXiv:hep-th/9904075 .
- Armendariz-Picon et al. [2001] C. Armendariz-Picon, V. F. Mukhanov, and P. J. Steinhardt, Essentials of k-essence, Phys. Rev. D63, 103510 (2001), arXiv:astro-ph/0006373 .
- Batista [2021] R. C. Batista, A Short Review on Clustering Dark Energy, Universe 8, 22 (2021), arXiv:2204.12341 [astro-ph.CO] .
- Batista [2014] R. C. Batista, Impact of dark energy perturbations on the growth index, Phys.Rev. D89, 123508 (2014), arXiv:1403.2985 [astro-ph.CO] .
- Mehrabi et al. [2015a] A. Mehrabi, S. Basilakos, and F. Pace, How clustering dark energy affects matter perturbations, Mon. Not. Roy. Astron. Soc. 452, 2930 (2015a), arXiv:1504.01262 [astro-ph.CO] .
- Mehrabi et al. [2015b] A. Mehrabi, S. Basilakos, M. Malekjani, and Z. Davari, Growth of matter perturbations in clustered holographic dark energy cosmologies, Phys. Rev. D 92, 123513 (2015b), arXiv:1510.03996 [astro-ph.CO] .
- Mehrabi et al. [2015c] A. Mehrabi, M. Malekjani, and F. Pace, Can observational growth rate data favor the clustering dark energy models?, Astrophys. Space Sci. 356, 129 (2015c), arXiv:1411.0780 [astro-ph.CO] .
- Chevallier and Polarski [2001] M. Chevallier and D. Polarski, Accelerating universes with scaling dark matter, Int. J. Mod. Phys. D10, 213 (2001), gr-qc/0009008 .
- Linder [2003] E. V. Linder, Exploring the expansion history of the universe, Phys. Rev. Lett. 90, 091301 (2003), astro-ph/0208512 .
- Favale et al. [2023] A. Favale, A. Gomez-Valent, and M. Migliaccio, Cosmic chronometers to calibrate the ladders and measure the curvature of the universe. a model-independent study, Monthly Notices of the Royal Astronomical Society 523, 3406 (2023), arXiv:2301.09591 [astro-ph.CO] .
- Foreman-Mackey et al. [2012] D. Foreman-Mackey, D. W. Hogg, D. Lang, and J. Goodman, emcee: The mcmc hammer, Publications of the Astronomical Society of the Pacific 125, 306 (2012), arXiv:1202.3665 [astro-ph.IM] .
- Lewis [2019] A. Lewis, GetDist: a Python package for analysing Monte Carlo samples, (2019), arXiv:1910.13970 [astro-ph.IM] .
- Moresco et al. [2012] M. Moresco, A. Cimatti, R. Jimenez, L. Pozzetti, G. Zamorani, M. Bolzonella, J. Dunlop, F. Lamareille, M. Mignoli, H. Pearce, P. Rosati, D. Stern, L. Verde, E. Zucca, C. Carollo, T. Contini, J.-P. Kneib, O. L. Fevre, S. Lilly, V. Mainieri, A. Renzini, M. Scodeggio, I. Balestra, R. Gobat, R. McLure, S. Bardelli, A. Bongiorno, K. Caputi, O. Cucciati, S. de la Torre, L. de Ravel, P. Franzetti, B. Garilli, A. Iovino, P. Kampczyk, C. Knobel, K. Kovac, J.-F. L. Borgne, V. L. Brun, C. Maier, R. Pello, Y. Peng, E. Perez-Montero, V. Presotto, J. Silverman, M. Tanaka, L. Tasca, L. Tresse, D. Vergani, O. Almaini, L. Barnes, R. Bordoloi, E. Bradshaw, A. Cappi, R. Chuter, M. Cirasuolo, G. Coppa, C. Diener, S. Foucaud, W. Hartley, M. Kamionkowski, A. Koekemoer, C. Lopez-Sanjuan, H. McCracken, P. Nair, P. Oesch, A. Stanford, and N. Welikala, Improved constraints on the expansion rate of the universe up to z 1.1 from the spectroscopic evolution of cosmic chronometers, Journal of Cosmology and Astroparticle Physics 2012 (08), 006.
- Moresco [2015] M. Moresco, Raising the bar: new constraints on the hubble parameter with cosmic chronometers at z 2, Monthly Notices of the Royal Astronomical Society: Letters 450, L16 (2015).
- Moresco et al. [2016] M. Moresco, L. Pozzetti, A. Cimatti, R. Jimenez, C. Maraston, L. Verde, D. Thomas, A. Citro, R. Tojeiro, and D. Wilkinson, A 6measurement of the hubble parameter at z 0.45: direct evidence of the epoch of cosmic re-acceleration, Journal of Cosmology and Astroparticle Physics 2016 (05), 014.
- Zhang et al. [2014] C. Zhang, H. Zhang, S. Yuan, S. Liu, T.-J. Zhang, and Y.-C. Sun, Four new observationalh(z) data from luminous red galaxies in the sloan digital sky survey data release seven, Research in Astronomy and Astrophysics 14, 1221 (2014).
- Jimenez et al. [2003] R. Jimenez, L. Verde, T. Treu, and D. Stern, Constraints on the equation of state of dark energy and the hubble constant from stellar ages and the cosmic microwave background, The Astrophysical Journal 593, 622 (2003).
- Simon et al. [2005] J. Simon, L. Verde, and R. Jimenez, Constraints on the redshift dependence of the dark energy potential, Physical Review D 71, 123001 (2005).
- Ratsimbazafy et al. [2017] A. L. Ratsimbazafy, S. I. Loubser, S. M. Crawford, C. M. Cress, B. A. Bassett, R. C. Nichol, and P. Vaisanen, Age-dating luminous red galaxies observed with the southern african large telescope, Monthly Notices of the Royal Astronomical Society 467, 3239 (2017).
- Stern et al. [2010] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski, and S. A. Stanford, Cosmic chronometers: constraining the equation of state of dark energy. i: H(z) measurements, Journal of Cosmology and Astroparticle Physics 2010 (02), 008.
- Borghi et al. [2022] N. Borghi, M. Moresco, and A. Cimatti, Toward a better understanding of cosmic chronometers: A new measurement of h(z) at z 0.7, The Astrophysical Journal Letters 928, L4 (2022).
- Huterer and Linder [2006] D. Huterer and E. V. Linder, Separating dark physics from physical darkness: Minimalist modified gravity vs. dark energy, Phys.Rev.D75:023519,2007 75, 023519 (2006), arXiv:astro-ph/0608681 [astro-ph] .
- Gómez-Valent et al. [2021] A. Gómez-Valent, Z. Zheng, L. Amendola, V. Pettorino, and C. Wetterich, Early dark energy in the pre- and postrecombination epochs, Phys. Rev. D 104, 083536 (2021), arXiv:2107.11065 [astro-ph.CO] .
- Abramo et al. [2009] L. Abramo, R. Batista, L. Liberato, and R. Rosenfeld, Physical approximations for the nonlinear evolution of perturbations in inhomogeneous dark energy scenarios, Phys.Rev. D79, 023516 (2009), arXiv:0806.3461 [astro-ph] .
- Sapone et al. [2009] D. Sapone, M. Kunz, and M. Kunz, Fingerprinting dark energy, Phys. Rev. D80, 083519 (2009), arXiv:0909.0007 [astro-ph.CO] .
- Creminelli et al. [2009] P. Creminelli, G. D’Amico, J. Norena, and F. Vernizzi, The effective theory of quintessence: the w<-1 side unveiled, JCAP 0902, 018, arXiv:0811.0827 [astro-ph] .
- Batista and Pace [2013] R. Batista and F. Pace, Structure formation in inhomogeneous Early Dark Energy models, JCAP 1306, 044, arXiv:1303.0414 [astro-ph.CO] .
- Adame et al. [2024] A. G. Adame et al. (DESI), Desi 2024 vi: Cosmological constraints from the measurements of baryon acoustic oscillations, (2024), arXiv:2404.03002 [astro-ph.CO] .