Can dark energy explain a high growth index?

Ícaro B. S. Cortês [email protected] Departamento de Engenharia de Computação e Automação, Universidade Federal do Rio Grande do Norte, Caixa Postal 1524, CEP 59078-970, Natal, Rio Grande do Norte, Brazil.    Ronaldo C. Batista Escola de Ciências e Tecnologia, Universidade Federal do Rio Grande do Norte, Caixa Postal 1524, CEP 59078-970, Natal, Rio Grande do Norte, Brazil.
(November 1, 2024)
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 γ𝛾\gammaitalic_γ, the so-called growth index. It is well-known that the ΛΛ\Lambdaroman_ΛCDM cosmology predicts γ=0.55𝛾0.55\gamma=0.55italic_γ = 0.55. However, using observational data, Ref. [1] measured a much higher γ=0.6330.024+0.025𝛾subscriptsuperscript0.6330.0250.024\gamma=0.633^{+0.025}_{-0.024}italic_γ = 0.633 start_POSTSUPERSCRIPT + 0.025 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.024 end_POSTSUBSCRIPT, excluding the ΛΛ\Lambdaroman_ΛCDM value within 3.7σ3.7𝜎3.7\sigma3.7 italic_σ. In this work, we analyze whether Dark Energy (DE) with the Equation of State (EoS) parameter described by the CPL parametrization can significantly modify γ𝛾\gammaitalic_γ with respect to the ΛΛ\Lambdaroman_Λ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 γ𝛾\gammaitalic_γ distribution, we use data from 32323232 Cosmic Chronometers, H(zH(zitalic_H ( italic_z), 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 γ𝛾\gammaitalic_γ samples compatible with the value determined in Ref. [1]. Therefore, explaining the measured value of γ𝛾\gammaitalic_γ is a challenge for DE models. Moreover, we present new fitting functions for γ𝛾\gammaitalic_γ, 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, H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, obtained locally by using the distance ladder methods [4] and globally, assuming the ΛΛ\Lambdaroman_ΛCDM model, from Cosmic Microwave Background (CMB) data [5]. There is also the S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT tension (where S8σ8Ωm0/0.3subscript𝑆8subscript𝜎8subscriptΩ𝑚00.3S_{8}\equiv\sigma_{8}\sqrt{\Omega_{m0}/0.3}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT ≡ italic_σ start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT square-root start_ARG roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT / 0.3 end_ARG, σ8subscript𝜎8\sigma_{8}italic_σ start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT is the amplitude of matter fluctuations at 8h1Mpc8superscript1Mpc8h^{-1}\text{Mpc}8 italic_h start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT Mpc and Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT 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 S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT 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]

f=dlnδmdlna,𝑓𝑑subscript𝛿𝑚𝑑𝑎f=\frac{d\ln\delta_{m}}{d\ln a}\,,italic_f = divide start_ARG italic_d roman_ln italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT end_ARG start_ARG italic_d roman_ln italic_a end_ARG , (1)

which can be parametrized as [2]

fΩmγ,similar-to-or-equals𝑓superscriptsubscriptΩ𝑚𝛾f\simeq\Omega_{m}^{\gamma}\,,italic_f ≃ roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT , (2)

where γ𝛾\gammaitalic_γ is a constant that depends weakly on the Equation of State (EoS) parameter of Smooth DE (SDE) models, w(t)=pde(t)/ρde(t)𝑤𝑡subscript𝑝𝑑𝑒𝑡subscript𝜌𝑑𝑒𝑡w\left(t\right)=p_{de}\left(t\right)/\rho_{de}\left(t\right)italic_w ( italic_t ) = italic_p start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_t ) / italic_ρ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_t ), and Ωm=Ωm(a)subscriptΩ𝑚subscriptΩ𝑚𝑎\Omega_{m}=\Omega_{m}\left(a\right)roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( italic_a ) is the matter energy density parameter. For the ΛΛ\Lambdaroman_ΛCDM cosmology, Refs. [8, 2] found γ=0.55𝛾0.55\gamma=0.55italic_γ = 0.55. However, using observational data, Ref. [1] found a much higher γ=0.633𝛾0.633\gamma=0.633italic_γ = 0.633 (implying a suppression of growth), excluding the ΛΛ\Lambdaroman_ΛCDM value by 3.7σ3.7𝜎3.7\sigma3.7 italic_σ. In the same work, it was also shown that a higher γ𝛾\gammaitalic_γ values reduces the S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT tension from 3.2σ3.2𝜎3.2\sigma3.2 italic_σ to 0.9σ0.9𝜎0.9\sigma0.9 italic_σ.

This high value of γ𝛾\gammaitalic_γ naturally raises the question of what theoretical mechanisms are capable of producing it. In the case of SDE, it was shown that γ𝛾\gammaitalic_γ can be accurately described by the fitting function [2],

γ(w1)={0.55+0.02(1+w1),w110.55+0.05(1+w1),w1>1,𝛾subscript𝑤1cases0.550.021subscript𝑤1subscript𝑤110.550.051subscript𝑤1subscript𝑤11\gamma(w_{1})=\begin{cases}0.55+0.02(1+w_{1}),&w_{1}\leq-1\\ 0.55+0.05(1+w_{1}),&w_{1}>-1\end{cases},italic_γ ( italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = { start_ROW start_CELL 0.55 + 0.02 ( 1 + italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) , end_CELL start_CELL italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≤ - 1 end_CELL end_ROW start_ROW start_CELL 0.55 + 0.05 ( 1 + italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) , end_CELL start_CELL italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > - 1 end_CELL end_ROW , (3)

where w1=w(z=1)subscript𝑤1𝑤𝑧1w_{1}=w\left(z=1\right)italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = italic_w ( italic_z = 1 ). This indicates that, for w11subscript𝑤11w_{1}\leq-1italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≤ - 1, we would need w1=3.15subscript𝑤13.15w_{1}=3.15italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 3.15 to get γ=0.633𝛾0.633\gamma=0.633italic_γ = 0.633, in contradiction with the parametrization condition. For w1>1subscript𝑤11w_{1}>-1italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT > - 1, we would need w1=0.66subscript𝑤10.66w_{1}=0.66italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0.66, which would generate a large DE density around z=1𝑧1z=1italic_z = 1. 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 γ>0.6𝛾0.6\gamma>0.6italic_γ > 0.6.

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 cs=1subscript𝑐𝑠1c_{s}=1italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 1, 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, cssubscript𝑐𝑠c_{s}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT 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 cs0subscript𝑐𝑠0c_{s}\rightarrow 0italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT → 0, 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 γ0.633similar-to-or-equals𝛾0.633\gamma\simeq 0.633italic_γ ≃ 0.633. We also find that CDE models can not raise the values of γ𝛾\gammaitalic_γ significantly, therefore concluding that DE models are quite challenged by observed value of γ𝛾\gammaitalic_γ.

As we will show, for the case of CDE, this result has a simple explanation. When w(t)>1𝑤𝑡1w\left(t\right)>-1italic_w ( italic_t ) > - 1, DE perturbations are have the tendency of being correlated with matter perturbations (δdeδmproportional-tosubscript𝛿𝑑𝑒subscript𝛿𝑚\delta_{de}\propto\delta_{m}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ∝ italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT) and anti-correlated for w(t)<1𝑤𝑡1w\left(t\right)<-1italic_w ( italic_t ) < - 1 (δdeδmproportional-tosubscript𝛿𝑑𝑒subscript𝛿𝑚\delta_{de}\propto-\delta_{m}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ∝ - italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT). In other to raise γ𝛾\gammaitalic_γ, DE perturbations must be anti-correlated with δmsubscript𝛿𝑚\delta_{m}italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT, causing a decrease in the gravitational potential and consequently slowing down the growth rate (higher γ𝛾\gammaitalic_γ). However, in the case w(t)<1𝑤𝑡1w\left(t\right)<-1italic_w ( italic_t ) < - 1, the DE energy density decreases rapidly at high-z𝑧zitalic_z, thus the overall impact is very limited. On the other hand, the case w(t)>1𝑤𝑡1w\left(t\right)>-1italic_w ( italic_t ) > - 1 can easily enhance the matter growth, thus providing a significantly lower γ𝛾\gammaitalic_γ values, as already shown in [16].

Based on the samples of γ𝛾\gammaitalic_γ that we have computed, we where able to construct a new parametrization γ=𝛾absent\gamma=italic_γ =γ(w1)𝛾subscript𝑤1\gamma\left(w_{1}\right)italic_γ ( italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ), which is more accurate than (3) and valid for larger parameter space. Moreover, for the first time, we present a γ𝛾\gammaitalic_γ 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

ds2=c2dt2+a(t)2[dr2+rd2Ω2],𝑑superscript𝑠2superscript𝑐2𝑑superscript𝑡2𝑎superscript𝑡2delimited-[]𝑑superscript𝑟2𝑟superscript𝑑2superscriptΩ2ds^{2}=-c^{2}dt^{2}+a(t)^{2}\left[dr^{2}+r{}^{2}d\Omega^{2}\right],italic_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = - italic_c start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_a ( italic_t ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT [ italic_d italic_r start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_r start_FLOATSUPERSCRIPT 2 end_FLOATSUPERSCRIPT italic_d roman_Ω start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] , (4)

where a(t)𝑎𝑡a(t)italic_a ( italic_t ) is the scale factor. As so, the square of the Hubble function can be written as

H2=(a˙a)2=H02(Ωm0a3+Ωde(a)),superscript𝐻2superscript˙𝑎𝑎2superscriptsubscript𝐻02subscriptΩ𝑚0superscript𝑎3subscriptΩ𝑑𝑒𝑎H^{2}=\left(\frac{\dot{a}}{a}\right)^{2}=H_{0}^{2}\left(\Omega_{m0}a^{-3}+% \Omega_{de}\left(a\right)\right)\,,italic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ( divide start_ARG over˙ start_ARG italic_a end_ARG end_ARG start_ARG italic_a end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT italic_a start_POSTSUPERSCRIPT - 3 end_POSTSUPERSCRIPT + roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_a ) ) , (5)

where Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT is the matter (dark matter plus baryons) density parameter now and Ωde(a)subscriptΩ𝑑𝑒𝑎\Omega_{de}\left(a\right)roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_a ) 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:

w(a)=w0+wa(1a).𝑤𝑎subscript𝑤0subscript𝑤𝑎1𝑎w(a)=w_{0}+w_{a}(1-a).italic_w ( italic_a ) = italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 1 - italic_a ) . (6)

Thus, the DE density parameter is given by

Ωde(a)=(1Ωm0)a3(1+w0+wa)exp(3wa(1a)).subscriptΩ𝑑𝑒𝑎1subscriptΩ𝑚0superscript𝑎31subscript𝑤0subscript𝑤𝑎3subscript𝑤𝑎1𝑎\Omega_{de}(a)=(1-\Omega_{m0})a^{-3(1+w_{0}+w_{a})}\exp\left(-3w_{a}(1-a)% \right).roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_a ) = ( 1 - roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT ) italic_a start_POSTSUPERSCRIPT - 3 ( 1 + italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ) end_POSTSUPERSCRIPT roman_exp ( - 3 italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ( 1 - italic_a ) ) . (7)

Besides the general case of free w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and wasubscript𝑤𝑎w_{a}italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT, we will also analyze the ΛΛ\Lambdaroman_Λ limit (w0=1subscript𝑤01w_{0}=-1italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 1 and wa=0subscript𝑤𝑎0w_{a}=0italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 0) and constant EoS case (free ww0𝑤subscript𝑤0w\equiv w_{0}italic_w ≡ italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and wa=0subscript𝑤𝑎0w_{a}=0italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 0). We refer to these models as CPL, ΛΛ\Lambdaroman_ΛCDM and w𝑤witalic_wCDM, respectively. In the general case, we have four free parameters: h=H0/(100km s1Mpc1)subscript𝐻0100superscriptkm s1superscriptMpc1h=H_{0}/\left(100\text{km s}^{-1}\text{Mpc}^{-1}\right)italic_h = italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / ( 100 km s start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT Mpc start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ), Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT, w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and wasubscript𝑤𝑎w_{a}italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT. Next, we discuss how to minimally constrain these parameters using H(z)𝐻𝑧H\left(z\right)italic_H ( italic_z ) 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 χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT given by

χ2=χnosys2+χsys2,superscript𝜒2superscriptsubscript𝜒nosys2superscriptsubscript𝜒sys2\chi^{2}=\chi_{{\rm nosys}}^{2}+\chi_{{\rm sys}}^{2}\,,italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = italic_χ start_POSTSUBSCRIPT roman_nosys end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_χ start_POSTSUBSCRIPT roman_sys end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (8)

where χnosys2superscriptsubscript𝜒nosys2\chi_{{\rm nosys}}^{2}italic_χ start_POSTSUBSCRIPT roman_nosys end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT is associated with the 17 data-points without systematics, provided by [28, 29, 30, 31, 32, 33], which reads

χnosys2=i(ΔHiσi)2,superscriptsubscript𝜒nosys2subscript𝑖superscriptΔsubscript𝐻𝑖subscript𝜎𝑖2\chi_{{\rm nosys}}^{2}=\sum_{i}\left(\frac{\Delta H_{i}}{\sigma_{i}}\right)^{2% }\,,italic_χ start_POSTSUBSCRIPT roman_nosys end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = ∑ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ( divide start_ARG roman_Δ italic_H start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG start_ARG italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (9)

where ΔHi=Hmodel(zi)Hdata(zi)Δsubscript𝐻𝑖subscript𝐻modelsubscript𝑧𝑖subscript𝐻datasubscript𝑧𝑖\Delta H_{i}=H_{\text{model}}\left(z_{i}\right)-H_{{\rm\text{data}}}\left(z_{i% }\right)roman_Δ italic_H start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = italic_H start_POSTSUBSCRIPT model end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) - italic_H start_POSTSUBSCRIPT data end_POSTSUBSCRIPT ( italic_z start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT ) and σisubscript𝜎𝑖\sigma_{i}italic_σ start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT is the corresponding uncertainty of data points. The other part of the χ2superscript𝜒2\chi^{2}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT associated with data points with systematics is given by

χsys2=ΔHT𝑪𝒐𝒗1ΔH,superscriptsubscript𝜒sys2Δsuperscript𝐻𝑇𝑪𝒐superscript𝒗1Δ𝐻\chi_{{\rm sys}}^{2}=\Delta\vec{H}^{T}\boldsymbol{Cov}^{-1}\Delta\vec{H}\,,italic_χ start_POSTSUBSCRIPT roman_sys end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = roman_Δ over→ start_ARG italic_H end_ARG start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT bold_italic_C bold_italic_o bold_italic_v start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT roman_Δ over→ start_ARG italic_H end_ARG , (10)

where ΔHΔ𝐻\Delta\vec{H}roman_Δ over→ start_ARG italic_H end_ARG is a vector with ΔHiΔsubscript𝐻𝑖\Delta H_{i}roman_Δ italic_H start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT components associated the 15 data points provided by Refs. [25, 26, 27], and 𝑪𝒐𝒗𝑪𝒐𝒗\boldsymbol{Cov}bold_italic_C bold_italic_o bold_italic_v 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 θ=(h,Ωm0,w0,wa)𝜃subscriptΩ𝑚0subscript𝑤0subscript𝑤𝑎\theta=\left(h,\Omega_{m0},w_{0},w_{a}\right)italic_θ = ( italic_h , roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT , italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT ), given by

(θ|D)P(θ)L(D|θ),proportional-toconditional𝜃𝐷𝑃𝜃𝐿conditional𝐷𝜃\mathcal{L}\left(\theta|D\right)\propto P\left(\theta\right)L\left(D|\theta% \right)\,,caligraphic_L ( italic_θ | italic_D ) ∝ italic_P ( italic_θ ) italic_L ( italic_D | italic_θ ) , (11)

where

L(D|θ)exp(χ22)proportional-to𝐿conditional𝐷𝜃superscript𝜒22L\left(D|\theta\right)\propto\exp\left(-\frac{\chi^{2}}{2}\right)italic_L ( italic_D | italic_θ ) ∝ roman_exp ( - divide start_ARG italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ) (12)

is the likelihood of the data D𝐷Ditalic_D given θ𝜃\thetaitalic_θ and P(θ)𝑃𝜃P\left(\theta\right)italic_P ( italic_θ ) indicates the priors assumed for the parameters, listed in Tab. 1. We also implement an additional flat prior to exclude cosmologies in which Ωde(z=1000)>0.01subscriptΩ𝑑𝑒𝑧10000.01\Omega_{de}\left(z=1000\right)>0.01roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z = 1000 ) > 0.01, where z=1/a1𝑧1𝑎1z=1/a-1italic_z = 1 / italic_a - 1 is the redshift. This condition is similar to the assumption w0+wa<0.1subscript𝑤0subscript𝑤𝑎0.1w_{0}+w_{a}<-0.1italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT < - 0.1 described by [34]. Allowing for DE models with non-negligible energy density at high-z𝑧zitalic_z, 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.

Table 1: Priors for all parameters.
Parameter Prior
hhitalic_h 𝒰[0.5,1]𝒰0.51\mathcal{U}[0.5,1]caligraphic_U [ 0.5 , 1 ]
Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT 𝒰[0.01,0.99]𝒰0.010.99\mathcal{U}[0.01,0.99]caligraphic_U [ 0.01 , 0.99 ]
w𝑤witalic_w 𝒰[3,1/3]𝒰313\mathcal{U}[-3,-1/3]caligraphic_U [ - 3 , - 1 / 3 ]
w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT 𝒰[3,1]𝒰31\mathcal{U}[-3,1]caligraphic_U [ - 3 , 1 ]
wasubscript𝑤𝑎w_{a}italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT 𝒰[3,3]𝒰33\mathcal{U}[-3,3]caligraphic_U [ - 3 , 3 ]

The results for ΛΛ\Lambdaroman_ΛCDM, w𝑤witalic_wCDM 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, Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT and hhitalic_h, in the ΛΛ\Lambdaroman_ΛCDM model. For w𝑤witalic_wCDM this is not case, mainly because w𝑤witalic_w is being limited by the prior assumed, but note that the posterior also shows some preference for w1similar-to-or-equals𝑤1w\simeq-1italic_w ≃ - 1. Fig. 3 also shows that CPL parameters, w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and wasubscript𝑤𝑎w_{a}italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT, can not be constrained by CC data. The posteriors of w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and wasubscript𝑤𝑎w_{a}italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT are dominated by the flat priors shown in table 1 together with the condition Ωde(z=1000)<0.01subscriptΩ𝑑𝑒𝑧10000.01\Omega_{de}\left(z=1000\right)<0.01roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z = 1000 ) < 0.01. The mean values for the marginalized 1D distributions are given in Tab. 2. The number of samples in our chains is N1000τsimilar-to𝑁1000𝜏N\sim 1000\tauitalic_N ∼ 1000 italic_τ, where τ𝜏\tauitalic_τ is the largest auto-correlation time of the parameters. This is well beyond the suggested N>50τ𝑁50𝜏N>50\tauitalic_N > 50 italic_τ described in [23] as convergence criteria.

Refer to caption
Figure 1: Marginalized posterior distributions for the ΛΛ\Lambdaroman_ΛCDM model parameters.
Refer to caption
Figure 2: Marginalized posterior distributions for the w𝑤witalic_wCDM models, with a shared background and distinct Smooth and Clustering Dark Energy growth index distributions. The dashed lines indicate the values of the ΛΛ\Lambdaroman_ΛCDM parameters.
Refer to caption
Figure 3: Marginalized posterior distributions for the CPL models, with a shared background and distinct Smooth and Clustering Dark Energy growth index distributions. The dashed lines indicate the values of the ΛΛ\Lambdaroman_ΛCDM parameters.

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 γ𝛾\gammaitalic_γ. As we will see, despite this large variation of the EoS parameters, γ𝛾\gammaitalic_γ 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, δmρm/ρ¯m1subscript𝛿𝑚subscript𝜌𝑚subscript¯𝜌𝑚1\delta_{m}\equiv\rho_{m}/\overline{\rho}_{m}-1italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ≡ italic_ρ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT / over¯ start_ARG italic_ρ end_ARG start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT - 1, is solely affected by the background expansion, as described by the equation

δm′′+(21+3w(a)Ωde(a)2)δma32Ωm(a)a2δm=0,superscriptsubscript𝛿𝑚′′213𝑤𝑎subscriptΩ𝑑𝑒𝑎2superscriptsubscript𝛿𝑚𝑎32subscriptΩ𝑚𝑎superscript𝑎2subscript𝛿𝑚0\delta_{m}^{\prime\prime}+\left(2-\frac{1+3w\left(a\right)\Omega_{de}\left(a% \right)}{2}\right)\frac{\delta_{m}^{\prime}}{a}-\frac{3}{2}\frac{\Omega_{m}% \left(a\right)}{a^{2}}\delta_{m}=0,italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT + ( 2 - divide start_ARG 1 + 3 italic_w ( italic_a ) roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_a ) end_ARG start_ARG 2 end_ARG ) divide start_ARG italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_ARG start_ARG italic_a end_ARG - divide start_ARG 3 end_ARG start_ARG 2 end_ARG divide start_ARG roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( italic_a ) end_ARG start_ARG italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0 , (13)

where the prime denotes a derivative with respect to the scale factor. Starting the integration at a high redshift, zi=1000subscript𝑧𝑖1000z_{i}=1000italic_z start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT = 1000, initial values can be computed using the analytical solution for a matter-dominated Einstein-de-Sitter (EdS) universe δmi=δmi/aisuperscriptsubscript𝛿𝑚𝑖subscript𝛿𝑚𝑖subscript𝑎𝑖\delta_{mi}^{\prime}=\delta_{mi}/a_{i}italic_δ start_POSTSUBSCRIPT italic_m italic_i end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT = italic_δ start_POSTSUBSCRIPT italic_m italic_i end_POSTSUBSCRIPT / italic_a start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT, where we assumed δmi>0subscript𝛿𝑚𝑖0\delta_{mi}>0italic_δ start_POSTSUBSCRIPT italic_m italic_i end_POSTSUBSCRIPT > 0.

In the case of CDE, δmsubscript𝛿𝑚\delta_{m}italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT is also affected by DE perturbations, δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT, whose effect is maximized for DE models with negligible sound speed. In this work, we consider the extreme case cs=0subscript𝑐𝑠0c_{s}=0italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0 based on the fluid description of DE and phenomenologically allowing for EoS which can be non-phantom (w(t)>1𝑤𝑡1w(t)>-1italic_w ( italic_t ) > - 1), phantom (w(t)<1𝑤𝑡1w(t)<-1italic_w ( italic_t ) < - 1) or transit between these regimes. In this context, matter and DE perturbations obey the following system of equations [15]:

{δm+qa2=0,δde3waδde+(1+w)qa2=0,q12(1+3wΩde)qa+32(Ωmδm+Ωdeδde)=0,casessuperscriptsubscript𝛿𝑚𝑞superscript𝑎2absent0superscriptsubscript𝛿𝑑𝑒3𝑤𝑎subscript𝛿𝑑𝑒1𝑤𝑞superscript𝑎2absent0superscript𝑞1213𝑤subscriptΩ𝑑𝑒𝑞𝑎32subscriptΩ𝑚subscript𝛿𝑚subscriptΩ𝑑𝑒subscript𝛿𝑑𝑒absent0\begin{cases}\delta_{m}^{\prime}+\frac{q}{a^{2}}&=0\,,\\ \delta_{de}^{\prime}-3\frac{w}{a}\delta_{de}+(1+w)\frac{q}{a^{2}}&=0\,,\\ q^{\prime}-\frac{1}{2}(1+3w\Omega_{de})\frac{q}{a}+\frac{3}{2}(\Omega_{m}% \delta_{m}+\Omega_{de}\delta_{de})&=0\,,\end{cases}{ start_ROW start_CELL italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + divide start_ARG italic_q end_ARG start_ARG italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_CELL start_CELL = 0 , end_CELL end_ROW start_ROW start_CELL italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - 3 divide start_ARG italic_w end_ARG start_ARG italic_a end_ARG italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT + ( 1 + italic_w ) divide start_ARG italic_q end_ARG start_ARG italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_CELL start_CELL = 0 , end_CELL end_ROW start_ROW start_CELL italic_q start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT - divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( 1 + 3 italic_w roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ) divide start_ARG italic_q end_ARG start_ARG italic_a end_ARG + divide start_ARG 3 end_ARG start_ARG 2 end_ARG ( roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT + roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ) end_CELL start_CELL = 0 , end_CELL end_ROW (14)

where qθ/H𝑞𝜃𝐻q\equiv\theta/Hitalic_q ≡ italic_θ / italic_H and θ=v𝜃𝑣\theta=\vec{\nabla}\cdot\vec{v}italic_θ = over→ start_ARG ∇ end_ARG ⋅ over→ start_ARG italic_v end_ARG 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 cs=0subscript𝑐𝑠0c_{s}=0italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0. To compute the initial conditions for DE perturbations, we use the solution valid for constant w𝑤witalic_w in matter-dominated era [36, 37, 38, 39]

δde=1+w13wδm.subscript𝛿𝑑𝑒1𝑤13𝑤subscript𝛿𝑚\delta_{de}=\frac{1+w}{1-3w}\delta_{m}.italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT = divide start_ARG 1 + italic_w end_ARG start_ARG 1 - 3 italic_w end_ARG italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT . (15)

Although solution (15) only gives a general qualitative behavior at low-z𝑧zitalic_z [39], it can help us to predict the influence of CDE on the growth of δmsubscript𝛿𝑚\delta_{m}italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT. When w(a)<1𝑤𝑎1w\left(a\right)<-1italic_w ( italic_a ) < - 1 we have the tendency δdeδmproportional-tosubscript𝛿𝑑𝑒subscript𝛿𝑚\delta_{de}\propto-\delta_{m}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ∝ - italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT, while for w(a)>1𝑤𝑎1w\left(a\right)>-1italic_w ( italic_a ) > - 1 we have δdeδmproportional-tosubscript𝛿𝑑𝑒subscript𝛿𝑚\delta_{de}\propto\delta_{m}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ∝ italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT. 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 δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT and δmsubscript𝛿𝑚\delta_{m}italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT might take some time to achieve the expected behavior. The general trend is clear: phantom EoS will induce negative δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT, lowering the clustering power encoded in the last term of Eq. (14), whereas non-phantom EoS enhances growth. Therefore, the desired higher values of γ𝛾\gammaitalic_γ could be obtained in phantom CDE models.

We solve Eqs. (13) and (14), then determine the growth rate f=dlnδm/dlna𝑓𝑑subscript𝛿𝑚𝑑𝑎f=d\ln\delta_{m}/d\ln aitalic_f = italic_d roman_ln italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT / italic_d roman_ln italic_a and fit a constant γ𝛾\gammaitalic_γ assuming f=Ωmγ(a)𝑓superscriptsubscriptΩ𝑚𝛾𝑎f=\Omega_{m}^{\gamma}(a)italic_f = roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT ( italic_a ) for all the posterior distribution obtained and the measured value. In this process, it is important to ask whether the fitted values of γ𝛾\gammaitalic_γ can accurately describe f𝑓fitalic_f. 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 gδ/a𝑔𝛿𝑎g\equiv\delta/aitalic_g ≡ italic_δ / italic_a. Accuracies better than 0.05%percent0.050.05\%0.05 % for ΛΛ\Lambdaroman_ΛCDM with Ωm0[0.22,1]subscriptΩ𝑚00.221\Omega_{m0}\in[0.22,1]roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT ∈ [ 0.22 , 1 ] were reported, and 1%percent11\%1 % when lowering the limit to Ωm0=0.01subscriptΩ𝑚00.01\Omega_{m0}=0.01roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT = 0.01. For w1𝑤1w\neq-1italic_w ≠ - 1, the paper reported an accuracy within 0.4%percent0.40.4\%0.4 % for w=0.5𝑤0.5w=-0.5italic_w = - 0.5 and for dynamical EoS 0.25%percent0.250.25\%0.25 % when w0=0.82subscript𝑤00.82w_{0}=-0.82italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = - 0.82 and wa=0.58subscript𝑤𝑎0.58w_{a}=0.58italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT = 0.58. Here, we first will check how accurately a constant γ𝛾\gammaitalic_γ can reproduce f𝑓fitalic_f. In the next section, we will propose and analyze a new fitting function for γ𝛾\gammaitalic_γ.

In this work, we report the distribution of percent residuals given by the Root Mean Square (RMS) of 1Ωmγnun(z)/fnun(z)1superscriptsubscriptΩ𝑚subscript𝛾nun𝑧subscript𝑓nun𝑧1-\Omega_{m}^{\gamma_{{\rm nun}}}\left(z\right)/f_{{\rm nun}}\left(z\right)1 - roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_γ start_POSTSUBSCRIPT roman_nun end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ( italic_z ) / italic_f start_POSTSUBSCRIPT roman_nun end_POSTSUBSCRIPT ( italic_z ), where fnunsubscript𝑓nunf_{{\rm nun}}italic_f start_POSTSUBSCRIPT roman_nun end_POSTSUBSCRIPT is the numerical solution and γnunsubscript𝛾nun\gamma_{{\rm nun}}italic_γ start_POSTSUBSCRIPT roman_nun end_POSTSUBSCRIPT the corresponding fitted gamma value. We compute the RMS percent residuals along 10101010 evenly spaced values of 0z20𝑧20\leq z\leq 20 ≤ italic_z ≤ 2, and show its distribution in Fig. 4. As can be seen, a very small fraction of models have residuals greater than 0.5%percent0.50.5\%0.5 %. The worst case occurs for CDE-CPL, reaching up to 4%percent44\%4 %, but still for a small fraction of the realizations. In general, all models present a concentration of residuals close to 0.15%percent0.150.15\%0.15 %. We also have checked that the largest residuals are mainly associated with low Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT values. This demonstrates that, considering the CPL EoS, the constant γ𝛾\gammaitalic_γ 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-z𝑧zitalic_z, the constant γ𝛾\gammaitalic_γ parametrization might not be adequate [16].

Refer to caption
Figure 4: Histogram of RMS percent residuals for the constant growth index fit, with respect to the numerical solution of f𝑓fitalic_f for each model, normalized by relative frequency.

IV Results for γ𝛾\gammaitalic_γ

The ΛΛ\Lambdaroman_ΛCDM posteriors shown in Fig. 1 display a great determination of the growth index at γΛ=0.5514±0.0010subscript𝛾Λplus-or-minus0.55140.0010\gamma_{\Lambda}=0.5514\pm 0.0010italic_γ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT = 0.5514 ± 0.0010, 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 γm=0.6330.024+0.025subscript𝛾𝑚superscriptsubscript0.6330.0240.025\gamma_{m}=0.633_{-0.024}^{+0.025}italic_γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.633 start_POSTSUBSCRIPT - 0.024 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.025 end_POSTSUPERSCRIPT, excluding the ΛΛ\Lambdaroman_ΛCDM value within 3.7σ3.7𝜎3.7\sigma3.7 italic_σ, also showing that this value effectively solves the S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT tension. Assuming normally distributed probabilities for the measured value and the ΛΛ\Lambdaroman_ΛCDM one, with mean and standard deviation given in Tab. 2, a simply quantification of the tension between these values is given by

#σ=|γmγΛ|σm2+σΛ2=3.26#𝜎subscript𝛾𝑚subscript𝛾Λsuperscriptsubscript𝜎𝑚2superscriptsubscript𝜎Λ23.26\#\sigma=\frac{|\gamma_{m}-\gamma_{\Lambda}|}{\sqrt{\sigma_{m}^{2}+\sigma_{% \Lambda}^{2}}}=3.26\,# italic_σ = divide start_ARG | italic_γ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT - italic_γ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT | end_ARG start_ARG square-root start_ARG italic_σ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_σ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG end_ARG = 3.26 (16)

where we used σm=0.025subscript𝜎𝑚0.025\sigma_{m}=0.025italic_σ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.025. If we assume the usual value γΛ=0.55subscript𝛾Λ0.55\gamma_{\Lambda}=0.55italic_γ start_POSTSUBSCRIPT roman_Λ end_POSTSUBSCRIPT = 0.55 and neglect its uncertainty, we get 3.32σ3.32𝜎3.32\sigma3.32 italic_σ.

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 ΛΛ\Lambdaroman_ΛCDM background to produce the γ𝛾\gammaitalic_γ constraints, assessing the possible values of γ𝛾\gammaitalic_γ in more general background and perturbative models will help us to identify scenarios in which γ𝛾\gammaitalic_γ can be increased and how likely this can happen.

We first analyze the correlations between the γ𝛾\gammaitalic_γ and the EoS parameters and how γ𝛾\gammaitalic_γ changes with respect to the ΛΛ\Lambdaroman_ΛCDM value. The clearest case is for w𝑤witalic_wCDM model, Fig. 2, but the following analysis also holds for CPL model. For w>1𝑤1w>-1italic_w > - 1 (non-phantom EoS) we have lower Ωm(z)subscriptΩ𝑚𝑧\Omega_{m}\left(z\right)roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( italic_z ) in the past with respect to the w=1𝑤1w=-1italic_w = - 1 case. In the SDE scenario, this behavior causes a suppression of growth, giving a higher γ𝛾\gammaitalic_γ. In the case of CDE, this correlation is inverted because w>1𝑤1w>-1italic_w > - 1 induces positive δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT, which will act as an extra source of the gravitational potential, enhancing the growth and lowering γ𝛾\gammaitalic_γ. For w<1𝑤1w<-1italic_w < - 1 (phantom EoS) we have higher Ωm(z)subscriptΩ𝑚𝑧\Omega_{m}\left(z\right)roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT ( italic_z ) in the past with respect to w=1𝑤1w=-1italic_w = - 1. In SDE case, this induces a lower γ𝛾\gammaitalic_γ. Again, CDE inverts this correlation because now negative δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT is induced, reducing the source of the gravitational potential, consequently increasing γ𝛾\gammaitalic_γ. 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 w1subscript𝑤1w_{1}italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT, these correlations become more evident in Fig. 7. Therefore we can summarize impact of DE model in the γ𝛾\gammaitalic_γ , with respect to the ΛΛ\Lambdaroman_ΛCDM value, as follows:

  • SDE, non-phantom EoS: γ>0.55𝛾0.55\gamma>0.55italic_γ > 0.55

  • SDE, phantom EoS: γ<0.55𝛾0.55\gamma<0.55italic_γ < 0.55

  • CDE, non-phantom EoS: γ<0.55𝛾0.55\gamma<0.55italic_γ < 0.55

  • CDE, phantom EoS: γ>0.55𝛾0.55\gamma>0.55italic_γ > 0.55

It is also interesting to check the frequency of positive or negative occurrences of δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT for the CDE scenario. In Fig. 5, we present the distribution of δde/δmsubscript𝛿𝑑𝑒subscript𝛿𝑚\delta_{de}/\delta_{m}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT / italic_δ start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT at z=0𝑧0z=0italic_z = 0 and z=0.5𝑧0.5z=0.5italic_z = 0.5. We can see a small preference for negative δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT, which is associated with the allowed values of the EoS parameters as follows. For non-phantom EoS, Ωde(z)subscriptΩ𝑑𝑒𝑧\Omega_{de}\left(z\right)roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z ) can be large at intermediate and high-z𝑧zitalic_z, a situation that is disfavored by data, e.g., [35]. In our analysis, this fact is mainly implemented with the prior Ωde(z=1000)<0.01subscriptΩ𝑑𝑒𝑧10000.01\Omega_{de}\left(z=1000\right)<0.01roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z = 1000 ) < 0.01. On the other hand, for phantom EoS, Ωde(z)subscriptΩ𝑑𝑒𝑧\Omega_{de}\left(z\right)roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z ) is very small at intermediate and high-z𝑧zitalic_z and can still can induce an adequate accelerated expansion at low-z𝑧zitalic_z. Therefore, the allowed parameter space for the EoS parameters has a larger fraction of phantom realizations. As a consequence, given the correlations between w𝑤witalic_w and δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT explained earlier, we see some preference for δde<0subscript𝛿𝑑𝑒0\delta_{de}<0italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT < 0.

Refer to caption
Figure 5: Distribution of perturbation ratios for w𝑤witalic_wCDM and CPL models at z=0𝑧0z=0italic_z = 0 (left panel) and z=0.5𝑧0.5z=0.5italic_z = 0.5 (right panel).

In Fig. 6 we show a direct comparison between the γ𝛾\gammaitalic_γ 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 γ<0.55𝛾0.55\gamma<0.55italic_γ < 0.55, whereas CDE prefers γ>0.55𝛾0.55\gamma>0.55italic_γ > 0.55. The γ𝛾\gammaitalic_γ distributions for non-ΛΛ\Lambdaroman_Λ models have greater variance and asymmetry. However, none of these models provide a significant fraction of realizations around γ0.633similar-to-or-equals𝛾0.633\gamma\simeq 0.633italic_γ ≃ 0.633.

Considering the 2σ2𝜎2\sigma2 italic_σ lower limit of the measured value as γmin=0.63320.025=0.583subscript𝛾min0.63320.0250.583\gamma_{\rm min}=0.633-2\cdot 0.025=0.583italic_γ start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT = 0.633 - 2 ⋅ 0.025 = 0.583, only the SDE-w𝑤witalic_wCDM model has values higher than this, what happens only above percentile 96 of its γ𝛾\gammaitalic_γ distribution. If we cut this γ𝛾\gammaitalic_γ distribution with Ωm0>0.2subscriptΩ𝑚00.2\Omega_{m0}>0.2roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT > 0.2, γminsubscript𝛾min\gamma_{\rm min}italic_γ start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT occurs only above percentile 99. As can be seen in Fig. 6, all the other models have much smaller probabilities of producing γminsubscript𝛾min\gamma_{\rm min}italic_γ start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT.

Refer to caption
Figure 6: Distribution for γ𝛾\gammaitalic_γ in different models. The vertical dashed line marks the measured value γ=0.633𝛾0.633\gamma=0.633italic_γ = 0.633, the grey bands indicate its 1σ1𝜎1\sigma1 italic_σ and 2σ2𝜎2\sigma2 italic_σ intervals, assuming a normal distribution.

We could expect that phantom CDE models are able to produce high γ𝛾\gammaitalic_γ values because, as explained earlier, the associated negative δdesubscript𝛿𝑑𝑒\delta_{de}italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT operates in this direction. However, given that the gravitational potential depends on ΩdeδdesubscriptΩ𝑑𝑒subscript𝛿𝑑𝑒\Omega_{de}\delta_{de}roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT and, in the phantom case, ΩdesubscriptΩ𝑑𝑒\Omega_{de}roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT is small at intermediate and high-z𝑧zitalic_z, the actual impact of these models on γ𝛾\gammaitalic_γ is very limited. Consequently, although some preference for γ>0.55𝛾0.55\gamma>0.55italic_γ > 0.55 can be seen in Fig. 6, the highest possible values virtually never reach γ0.58greater-than-or-equivalent-to𝛾0.58\gamma\gtrsim 0.58italic_γ ≳ 0.58.

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 γ=0.6330.024+0.025𝛾subscriptsuperscript0.6330.0250.024\gamma=0.633^{+0.025}_{-0.024}italic_γ = 0.633 start_POSTSUPERSCRIPT + 0.025 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.024 end_POSTSUBSCRIPT, found in Ref. [1]. Bear in mind that in order to produce the γ𝛾\gammaitalic_γ distributions we only used CC data, the priors in Tab. 1 and the condition Ωde(z=1000)<0.01subscriptΩ𝑑𝑒𝑧10000.01\Omega_{de}\left(z=1000\right)<0.01roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z = 1000 ) < 0.01. Therefore, adding more background data in the analysis can only constrain more the parameter space and consequently the γ𝛾\gammaitalic_γ distribution.

The mean values for the marginalized 1D distributions are listed in Tab. 2. Note that, as the γ𝛾\gammaitalic_γ distribution is heavy tailed for CDE models, and reported mean value can lie outside the 1σ1𝜎1\sigma1 italic_σ interval. We verified that the best fit model parameters lie within 1σ1𝜎1\sigma1 italic_σ of the 1D marginalized distributions.

Table 2: Marginalized 1D constraints and 68%percent6868\%68 % C.L. intervals for the parameters of each cosmology.
Parameter ΛΛ\Lambdaroman_ΛCDM w𝑤witalic_wCDM CPL
hhitalic_h 0.689±0.041plus-or-minus0.6890.0410.689\pm 0.0410.689 ± 0.041 0.7270.092+0.054superscriptsubscript0.7270.0920.0540.727_{-0.092}^{+0.054}0.727 start_POSTSUBSCRIPT - 0.092 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.054 end_POSTSUPERSCRIPT 0.7330.10+0.069superscriptsubscript0.7330.100.0690.733_{-0.10}^{+0.069}0.733 start_POSTSUBSCRIPT - 0.10 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.069 end_POSTSUPERSCRIPT
Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT 0.3240.072+0.049superscriptsubscript0.3240.0720.0490.324_{-0.072}^{+0.049}0.324 start_POSTSUBSCRIPT - 0.072 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.049 end_POSTSUPERSCRIPT 0.296±0.076plus-or-minus0.2960.0760.296\pm 0.0760.296 ± 0.076 0.3190.082+0.060superscriptsubscript0.3190.0820.0600.319_{-0.082}^{+0.060}0.319 start_POSTSUBSCRIPT - 0.082 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.060 end_POSTSUPERSCRIPT
w𝑤witalic_w or w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - - - 1.420.37+0.97superscriptsubscript1.420.370.97-1.42_{-0.37}^{+0.97}- 1.42 start_POSTSUBSCRIPT - 0.37 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.97 end_POSTSUPERSCRIPT 1.470.72+0.91superscriptsubscript1.470.720.91-1.47_{-0.72}^{+0.91}- 1.47 start_POSTSUBSCRIPT - 0.72 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.91 end_POSTSUPERSCRIPT
wasubscript𝑤𝑎w_{a}italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT - - - - - - 0.62.0+1.1superscriptsubscript0.62.01.1-0.6_{-2.0}^{+1.1}- 0.6 start_POSTSUBSCRIPT - 2.0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 1.1 end_POSTSUPERSCRIPT
γ𝛾\gammaitalic_γ (SDE) 0.5514±0.0010plus-or-minus0.55140.00100.5514\pm 0.00100.5514 ± 0.0010 0.54650.021+0.0054superscriptsubscript0.54650.0210.00540.5465_{-0.021}^{+0.0054}0.5465 start_POSTSUBSCRIPT - 0.021 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.0054 end_POSTSUPERSCRIPT 0.54120.015+0.0056superscriptsubscript0.54120.0150.00560.5412_{-0.015}^{+0.0056}0.5412 start_POSTSUBSCRIPT - 0.015 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.0056 end_POSTSUPERSCRIPT
γ𝛾\gammaitalic_γ (CDE) - - - 0.545+0.0061+0.019superscriptsubscript0.5450.00610.0190.545_{+0.0061}^{+0.019}0.545 start_POSTSUBSCRIPT + 0.0061 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.019 end_POSTSUPERSCRIPT 0.5560.00039+0.013superscriptsubscript0.5560.000390.0130.556_{-0.00039}^{+0.013}0.556 start_POSTSUBSCRIPT - 0.00039 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + 0.013 end_POSTSUPERSCRIPT

New fits for γ(w)𝛾𝑤\gamma(w)italic_γ ( italic_w )

With the posterior distributions obtained, we were able find a more accurate and general fitting function for γ𝛾\gammaitalic_γ, given by

γ(w1)=aw1+becw1+d,𝛾subscript𝑤1𝑎subscript𝑤1𝑏superscript𝑒𝑐subscript𝑤1𝑑\gamma(w_{1})=aw_{1}+be^{cw_{1}}+d,italic_γ ( italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = italic_a italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_b italic_e start_POSTSUPERSCRIPT italic_c italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT + italic_d , (17)

The coefficients for each model are listed in Tab. 3. The tested range of w1subscript𝑤1w_{1}italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is [3,1/3]313[-3,-1/3][ - 3 , - 1 / 3 ]. The main difference between SDE and CDE models is the sign of the coefficient in the exponential term.

Refer to caption
Refer to caption
Figure 7: Different fits of γ(w1)𝛾subscript𝑤1\gamma(w_{1})italic_γ ( italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) for w𝑤witalic_wCDM models (left panels) and for CPL models (right panels) with contours for γ(w1)𝛾subscript𝑤1\gamma(w_{1})italic_γ ( italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) for SDE (blue) and CDE (orange) in the the upper panels. The lower panels show histograms, normalized by relative frequency, of the RMS percent residual between the γ𝛾\gammaitalic_γ given by fitting functions and the corresponding numerical value, obtained from the parametrization f=Ωmγ𝑓superscriptsubscriptΩ𝑚𝛾f=\Omega_{m}^{\gamma}italic_f = roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_γ end_POSTSUPERSCRIPT.
Table 3: Coefficients for the fits of γ(w1)=aw1+becw1+d𝛾subscript𝑤1𝑎subscript𝑤1𝑏superscript𝑒𝑐subscript𝑤1𝑑\gamma(w_{1})=aw_{1}+be^{cw_{1}}+ditalic_γ ( italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) = italic_a italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + italic_b italic_e start_POSTSUPERSCRIPT italic_c italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT + italic_d for each collapse model and cosmology.
Model Cosmology a𝑎aitalic_a b𝑏bitalic_b c𝑐citalic_c d𝑑ditalic_d
SDE w𝑤witalic_wCDM 0.01010.01010.01010.0101 0.12660.12660.12660.1266 2.80132.80132.80132.8013 0.55320.55320.55320.5532
SDE CPL 0.00680.00680.00680.0068 0.13460.13460.13460.1346 2.84532.84532.84532.8453 0.54930.54930.54930.5493
CDE w𝑤witalic_wCDM 0.00470.00470.00470.0047 0.63200.6320-0.6320- 0.6320 3.85703.85703.85703.8570 0.57060.57060.57060.5706
CDE CPL 0.00450.00450.00450.0045 0.46730.4673-0.4673- 0.4673 3.22503.22503.22503.2250 0.57380.57380.57380.5738

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 w𝑤witalic_wCDM model (left panels), presenting a mode of residuals at 0.1%percent0.10.1\%0.1 %, while the linear fit of Eq. (3) presents a mode around 0.25%percent0.250.25\%0.25 %. In any case, very few realizations have with residuals >1%absentpercent1>1\%> 1 %. 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 1%percent11\%1 % while the exponential fit, Eq. (17), peaks at 0.2%percent0.20.2\%0.2 %. The main reason for this improvement is that the linear fit does not perform well for w1<2subscript𝑤12w_{1}<-2italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < - 2, as can be seen in the top panels of Fig. 4.

Our fitting function for γ𝛾\gammaitalic_γ 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 0.1%percent0.10.1\%0.1 % for both w=const.𝑤const.w=\text{const.}italic_w = const. and CPL parametrizations.

As an example of the use of our fit, in Tab. 4, we present the γ𝛾\gammaitalic_γ 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 γ𝛾\gammaitalic_γ value, but note that CDE can provide a slightly higher γ𝛾\gammaitalic_γ because the w0wasubscript𝑤0subscript𝑤𝑎w_{0}w_{a}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT values have a phantom behavior at high an intermediate z𝑧zitalic_z. It is interesting to note that all these three DESI data analysis provide w11.2similar-to-or-equalssubscript𝑤11.2w_{1}\simeq-1.2italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≃ - 1.2, thus the γ𝛾\gammaitalic_γ values for SDE or CDE are quite similar.

Table 4: Values of γ𝛾\gammaitalic_γ, given by Eq. (17), for DESI best-fit values for CPL parametrization using three different data combinations described in Ref. [40].
w0subscript𝑤0w_{0}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT wasubscript𝑤𝑎w_{a}italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT γ𝛾\gammaitalic_γ-SDE γ𝛾\gammaitalic_γ-CDE
DESI+CMB+Pahnth. -0.872 -0.75 0.546 0.559
DESI+CMB+Union3 -0.65 -1.27 0.544 0.561
DESI+CMB+DESY5 -0.727 -1.05 0.545 0.560

V Conclusions

In this paper, we have analyzed the possibilities of DE models described by the CPL parametrization of producing a high value of γ𝛾\gammaitalic_γ, compatible with γ=0.6330.024+0.025𝛾subscriptsuperscript0.6330.0250.024\gamma=0.633^{+0.025}_{-0.024}italic_γ = 0.633 start_POSTSUPERSCRIPT + 0.025 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.024 end_POSTSUBSCRIPT obtained in Ref. [1]. As summarized in Fig. 6, both smooth DE and clustering DE have γ𝛾\gammaitalic_γ distributions with small or negligible overlap with the 2σ2𝜎2\sigma2 italic_σ interval of the measurement. This result depends only on the following assumptions: (1) H(z)𝐻𝑧H\left(z\right)italic_H ( italic_z ) must be compatible with CC data, (2) Ωde(z=1000)<0.01subscriptΩ𝑑𝑒𝑧10000.01\Omega_{de}\left(z=1000\right)<0.01roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z = 1000 ) < 0.01 and (3) the priors described in Tab. 1. The combination of these assumptions imposes very loose constraints on w0wasubscript𝑤0subscript𝑤𝑎w_{0}w_{a}italic_w start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_w start_POSTSUBSCRIPT italic_a end_POSTSUBSCRIPT parameters, allowing for a vast, but meaningful, exploration of the possible values of γ𝛾\gammaitalic_γ.

We have also analyzed the correlations between γ𝛾\gammaitalic_γ, the EoS parameters and the clustering properties of DE. In particular, we can expect that phantom CDE models can raise the value of γ𝛾\gammaitalic_γ. However, since the impact of such models depend on Ωde(z)δdesubscriptΩ𝑑𝑒𝑧subscript𝛿𝑑𝑒\Omega_{de}\left(z\right)\delta_{de}roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z ) italic_δ start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT, the actual change is very limited because Ωde(z)subscriptΩ𝑑𝑒𝑧\Omega_{de}\left(z\right)roman_Ω start_POSTSUBSCRIPT italic_d italic_e end_POSTSUBSCRIPT ( italic_z ) decays rapidly with z𝑧zitalic_z. Non-phantom SDE can also give higher γ𝛾\gammaitalic_γ, but these realizations are correlated with very low Ωm00.1similar-to-or-equalssubscriptΩ𝑚00.1\Omega_{m0}\simeq 0.1roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT ≃ 0.1.

We have also proposed a new fitting function for γ=γ(w1)𝛾𝛾subscript𝑤1\gamma=\gamma(w_{1})italic_γ = italic_γ ( italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ), with overall accuracies better than the linear fit of Eq. (3) covering the interval of 3w11/33subscript𝑤113-3\leq w_{1}\leq-1/3- 3 ≤ italic_w start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ≤ - 1 / 3 . For the first time, we produced a γ𝛾\gammaitalic_γ fitting function for CDE models, Eq. (17) with coefficients given in Tab. 3. These fits can be useful for a fast estimation of γ𝛾\gammaitalic_γ for both SDE and CDE models described by CPL parametrization.

If the high γ𝛾\gammaitalic_γ 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 cssubscript𝑐𝑠c_{s}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT on small scales. Since we have considered the two limiting cases of smooth DE (cs=1subscript𝑐𝑠1c_{s}=1italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 1) and full clustering DE (cs=0subscript𝑐𝑠0c_{s}=0italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 0), it seems unlikely that intermediate or time-varying cssubscript𝑐𝑠c_{s}italic_c start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT values can produce a significantly higher γ𝛾\gammaitalic_γ. A possible alternative is to consider more general EoS parametrizations. However, as demonstrated in our study, large variations on the EoS parameters and Ωm0subscriptΩ𝑚0\Omega_{m0}roman_Ω start_POSTSUBSCRIPT italic_m 0 end_POSTSUBSCRIPT produce much narrower γ𝛾\gammaitalic_γ distributions. Therefore, in principle, even more general EoS parametrizations should have some difficulties in explaining such high γ𝛾\gammaitalic_γ values, with the disadvantage of introducing more parameters. If this ‘γ𝛾\gammaitalic_γ 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