Free Access
Issue
A&A
Volume 616, August 2018
Article Number A32
Number of page(s) 9
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201833313
Published online 08 August 2018

© ESO 2018

1. Introduction

Evidence for the recent acceleration of the universe has been steadily accumulating over the past two decades (Riess et al. 1998; Perlmutter et al. 1999). Nevertheless, the properties of the dark energy hypothesised to be responsible for this acceleration (not to mention its ultimate origin) remain mostly unknown, and a plethora of possible theoretical explanations and phenomenological models has been proposed (Copeland et al. 2006). In particular, one would like to know whether this is due to a cosmological constant (as first introduced, in a different context, by Einstein) or to a dynamical degree of freedom such as a scalar field.

One pragmatic observational approach to the problem consists of mapping the behaviour of dark energy as a function of redshift and trying to identify any dynamics, and this is commonly done in terms of the dark energy equation of state w(z) = p(z)/ρ(z), with a constant wΛ = −1 corresponding to a cosmological constant – see, for example, the recent review by Huterer & Shafer (2018). Indeed, it is known that the present-day value of the dark energy equation of state, w0, must be very close to −1 (Planck Collaboration XIV 2016), but the constraints at higher redshifts are weaker.

It is well known that the early universe went through several phase transitions. However, in some cosmological models these phase transitions can also occur in the more recent universe. Examples include models involving topological defects (an early suggestion being Hill et al. 1989), the vacuum metamorphosis scenario of Parker & Raval (1999, 2000), in which non-perturbative quantum effects can have observationally significant consequences at late times, and scalar field models with a non-canonical kinetic term (Mortonson et al. 2009). Such phase transitions would affect the properties of the universe’s contents, and in particular they could change the behaviour of the dark energy equation of state.

Several phenomenological parametrisations of the dark energy equation of state have been proposed and studied, which allow for a fast transition in the dark energy equation of state, including Bassett et al. (2002), Corasaniti & Copeland (2003), and Linder & Huterer (2005). Our goal in the present work is to revisit these and other more recent works (Lazkoz et al. 2016; Jaber & de la Macorra 2016, 2018; Marcondes & Pan 2017), focusing on low redshift constraints on these models, by using more recent data but also by using some simplifying physical arguments and hypotheses, to be further described in the next section.

Most of the previous analyses use a combination of low redshift data (such as Type Ia supernovas) and cosmic microwave background (CMB) measurements; it is well known that the latter are particularly constraining. In the present work we have restricted ourselves to using low redshift data: specifically we used data from Type Ia supernovas (Suzuki et al. 2012) and measurements of the Hubble parameter (Farooq et al. 2017). While this leads to somewhat weaker constraints, our goal is to complement the existing literature by quantifying the constraining power of the current low redshift data. That said, we also discuss the impact of the CMB and compare our results to those of previous works.

We note that recently Di Valentino et al. (2018) have studied a different vacuum phase transition scenario, where the model parameters are chosen such that the universe contains matter plus a cosmological constant at early times and the phase transition leads to phantom dark energy today (that is, to w0 < −1), which asymptotes to a cosmological constant in the far future. Relying mainly on CMB data, these authors suggest that such a model may alleviate the aforementioned Hubble constant tension. By contrast, the models we study always have canonical equations of state (w(z) ≥ −1), so in this sense the present work complements theirs although the motivations of both are somewhat different.

2. Two-parameter dark energy equation of state parametrisations

In order to parametrise a transition in the dark energy equation of state, one generically needs four free parameters: the values of the equation of state itself at early and late times (respectively denoted wi and wf), plus the characteristic redshift at which the transition occurs, zt, and a transition width Δ (describing whether the transition is fast or slow). We note that the present-day value of the dark energy equation of state, w0, can be expressed in terms of these four parameters, or alternatively replace one of them.

Three such phenomenological parametrisations have been described and studied in previous works:

  • Bassett et al. (2002) have parametrised the transition as a function of redshift,

    w ( z ) = w i + w b w i 1 + exp [ z z t Δ ] . $$ w{(z)}=w_\text{i}+\frac{w_\text{b}-w_\text{i}}{1+exp{\lbrack\frac{z-z_\text{t}}{\mathrm\Delta}\rbrack}}. $$(1)

    We note that Bassett et al. (2002) state that their parameter wb corresponds to the late time behaviour of their equation of state, a parameter which we denote wf (and in their paper they indeed denote it wf in the above equation). This is only approximately correct; the exact relation between the two is

    w b = w f ( 1 + e ( 1 + z t ) / Δ ) w i e ( 1 + z t ) / Δ . $$ w_\text{b}=w_\text{f}{(1+e^{-{(1+z_\text{t})/\mathrm\Delta}})}-w_\text{i}e^{-{(1+z_\text{t})/\mathrm\Delta}}. $$(2)

  • Corasaniti & Copeland (2003) suggest an analogous parametrisation, but expressed as a function of the scale factor,

    w ( a ) = w 0 + ( w i w 0 ) 1 + e a t / Δ 1 e 1 / Δ 1 exp [ 1 a Δ ] 1 + exp [ a t a Δ ] $$ w{(a)}=w_0+{(w_\text{i}-w_0)\frac{1+e^{a_t/\mathrm\Delta}}{1-e^{1/\mathrm\Delta}}}\frac{1-exp{\lbrack\frac{1-a}{\mathrm\Delta}\rbrack}}{1+exp{\lbrack\frac{a_\text{t}-a}{\mathrm\Delta}\rbrack}} $$(3)

  • Linder & Huterer (2005; see also Linder 2008) later suggested an alternative parametrisation, which has the practical advantage of allowing an analytic expression for the Hubble parameter,

    w ( a ) = w f + w i w f 1 + [ a a t ] 1 Δ $$ w{(a)}=w_\text{f}+\frac{w_\text{i}-w_\text{f}}{1+{\lbrack\frac a{a_\text{t}}\rbrack}^{1-\mathrm\Delta}} $$(4)

In what follows we have used the three parametrisations, as a means to ascertain how much the constraints derived from current data depend on this choice. On the other hand, we introduce two physically reasonable simplifying assumptions. Firstly, and considering the fact that current data can only tightly constrain two dark energy parameters (see for example the discussion in Linder & Huterer 2005) in the next section we restrict ourselves to two-parameter models whose dark energy equation of state behaves as matter in the asymptotic past, and as a cosmological constant in the asymptotic future: in other words, we choose two of the free parameters to have the values

w i = 0 $$ w_\text{i}=0 $$(5)

and

w f = 1. $$ w_\text{f}=-1. $$(6)

The two remaining parameters can then be chosen to be w0 and Δ—in other words, the present-day value of the dark energy equation of state and a measure of how fast this value changes. The transition redshift can then be obtained from both of these. For the last model the wi = 0 assumption is relaxed later on, leading to a three-parameter model.

Secondly, we restrict ourselves to canonical models, having a dark energy equation of state with w(z) ≥ −1; as previously mentioned, the case of phantom models, having w(z) ≥ −1, was recently studied by Di Valentino et al. (2018). Re-writing the above expressions in terms of our chosen free parameters, we obtain

  • For the Bassett et al. (2002) parametrisation,

    w ( z ) = w 0 ( 1 e 1 / Δ ) ( 1 + w 0 ) e z / Δ + ( w 0 + e 1 / Δ ) , $$ w{(z)}=\frac{w_0{(1-e^{-1/\mathrm\Delta})}}{(1+w_0)e^{z/\mathrm\Delta}+{(w_0+e^{-1/\mathrm\Delta})}}, $$(7)

    with the transition redshift given by

    z t = Δ log [ w 0 + e 1 / Δ ( 1 + w 0 ) ] ; $$ z_\text{t}=\mathrm\Delta log{\lbrack\frac{w_0+e^{-1/\mathrm\Delta}}{-{(1+w_0)}}\rbrack}; $$(8)

    we henceforth denote this specific parametrisation as Model B.

  • For the Corasaniti & Copeland (2003) parametrisation, henceforth denoted Model C, we have

    w ( a ) = w 0 ( 1 e a / Δ ) ( 1 + w 0 ) e ( 1 a ) / Δ + ( w 0 + e a / Δ ) , $$ w{(a)}=\frac{w_0{(1-e^{-a/\mathrm\Delta})}}{(1+w_0)e^{{(1-a)}/\mathrm\Delta}+{(w_0+e^{-a/\mathrm\Delta})}}, $$(9)

    with the transition redshift

    a t = 1 1 + z t = Δ log [ 1 + w 0 w 0 ( e 1 / Δ 1 ) 1 ] ; $$ a_\text{t}=\frac1{1+z_\text{t}}=\mathrm\Delta log{\left[\frac{1+w_0}{-w_0}{(e^{1/\mathrm\Delta}-1)-1}\right];} $$(10)

  • Finally, for the Linder & Huterer (2005) model, henceforth denoted Model L, we have

    w ( z ) = w 0 ( 1 + w 0 ) ( 1 + z ) 1 / Δ w 0 , $$ w{(z)}=\frac{w_0}{(1+w_0){(1+z)}^{1/\mathrm\Delta}-w_0}, $$(11)

    with a transition redshift

    z t = ( w 0 1 + w 0 ) Δ 1. $$ z_\text{t}={(\frac{-w_0}{1+w_0})}^\mathrm\Delta-1. $$(12)

We note that there are in principle choices of parameters for which one would obtain a negative transition redshift, implying that in such cases the effective epoch of transition only occurs in the future. However, as we see in the following section, such choices are observationally excluded. Figure 1 illustrates the redshift dependence of w(z) in our three models, for various choices of the free parameters w0 and Δ, while Fig. 2 shows how the transition redshift depends on the two free parameters; we note that different ranges of the parameters w0 and Δ are depicted in the various panels of this figure, since different parameter choices will lead to reasonable values of transition redshifts in each of the models.

thumbnail Fig. 1.

Examples of the redshift dependence of the dark energy equation of state in the three models considered in this work, for various choices of the free parameters in the models. Blue dashed, red dotted, and black solid lines correspond, respectively, to the Models B, C and L described in the text – cf. respectively, Eqs. (7), (9), and (11). The top left panel has w0 = −0.9 and Δ = 0.2, the top right one has w0 = −0.99 and Δ = 0.2, the bottom left one has w0 = −0.9 and Δ = 0.8, and the bottom right one has w0 = −0.99 and Δ = 0.8.

thumbnail Fig. 2.

Dependence of the transition redshift zt on the present-day value of the dark energy equation of state w0 and the transition width Δ in the three models considered in this work. The top, middle and bottom panels correspond, respectively, to the Models B, C, and L described in the text – cf. respectively, Eqs. (7)(12). In all panels the colourmap denotes the value of zt corresponding to each pair of values of w0 and Δ (we note that the colourmap and the ranges of w0 and Δ are different in each panel) and some specific contours are also identified by black lines with white labels.

We emphasise that the models have in common the fact that they are two-parameter models for dynamical dark energy, described by its current value, w0, and a measure of how fast this changes, Δ. Naturally, w0 has the same physical meaning in all three models. On the other hand, Δ plays the same qualitative role in all three models (describing how fast the dark energy equation of state evolves, near the present day), but it is clear that its quantitative meaning is different in each model, and therefore this will also be the case for the characteristic redshift of the transition, zt, which can be obtained from it and w0. This should also be clear from Figs. 1 and 2. Thus a direct comparison of the constraints on these parameters for the different models is instructive, but should not be taken literally.

We further assume flat Friedmann–Lemaître–Robertson–Walker models, so that the present-day values of the matter and dark energy densities (expressed as a function of the critical density) satisfy Ωm + Ωdark = 1. We neglect the contribution of the radiation density, since we will only be concerned with low-redshift data. Thus the Friedmann equation will have the general form

H 2 H 0 2 = Ω m ( 1 + z ) 3 + ( 1 Ω m ) exp [ 3 0 z 1 + w ( z ) 1 + z d z ] ; $$ \frac{H^2}{H_0^2}={\mathrm\Omega}_\text{m}{(1+z)}^3+{(1-{\mathrm\Omega}_\text{m})}exp{\left[3\int_0^z\frac{1+w{(z^')}}{1+z^'}dz^'\right]}; $$(13)

for the particular case of Model L, the substitution of Eq. (11) leads to the analytic expression

H 2 H 0 2 = Ω m ( 1 + z ) 3 + ( 1 Ω m ) [ ( 1 + w 0 ) ( 1 + z ) 1 / Δ w 0 ] 3 Δ . $$ \frac{H^2}{H_0^2}={\mathrm\Omega}_\text{m}{(1+z)}^3+{(1-{\mathrm\Omega}_\text{m})}{\lbrack{(1+w_0){(1+z)}^{1/\mathrm\Delta}}-w_0\rbrack}^{3\mathrm\Delta}. $$(14)

3. Low-redshift constraints on two-parameter models from current data

Most previous works studying these scenarios constrain them using a combination of low-redshift data (such as Type Ia supernovas) and CMB measurements. It is well known that in such cases the latter will carry most of the statistical weight in the analysis, given its high accuracy and the fact that the combination of the two observables encompasses a very wide redshift range. In the present work we will restrict ourselves to using low-redshift data, quantifying how constraining this data currently is. In the discussion section we will compare our results to those of other recent works which do use the CMB.

Specifically, we use the recent Union2.1 catalogue of Type Ia supernovas (Suzuki et al. 2012) together with the compilation of measurements of the Hubble parameter by Farooq et al. (2017). The latter has the advantage of extending our redshift lever arm, since it includes measurements up to red-shift z ~ 2.36, while the supernova data is all at z < 1.5 (and indeed mostly at z < 1). On the other hand, we emphasise that this compilation of Hubble parameter measurements is heterogeneous: some of the measurements come from galaxy clustering and baryon acoustic oscillations (BAO) observations, while others come from the so-called cosmic chronometers or differential age method (Jimenez & Loeb 2002). We note that it is not currently clear that possible systematics issues of the cosmic chronometers method are well understood and under control (Liu et al. 2016; Lopez-Corredoira et al. 2017; Concas et al. 2017; Lopez-Corredoira & Vazdekis 2018).

We have carried out a standard statistical likelihood analysis, assuming a flat logarithmic prior on (1 + w0), specifically log10(1 + w0) = [−3, −0.5], and a uniform prior on the transition width, Δ = ]0, 1]. Our choice of prior for w0 is motivated by the fact that (as will be seen in what follows) the data strongly prefers a value of very close to that of a cosmological constant, and that for values smaller than log10(1 + w0) = −3 the phase space volume has almost the same likelihood. It should be noted that this choice of a logarithmic prior for a parameter whose value is near zero may affect the resulting posterior, due to the effect of prior volume weighting.

We assume a flat universe, and further fix the present-day value of the matter density to be Ωm = 0.3, in agreement with both CMB and low-redshift data (Planck Collaboration XIII 2016; Abbott et al. 2017; Scolnic et al. 2018; Jones et al. 2018). We initially fix the Hubble constant to be H0 = 70 km s−1 Mpc−1, but we will subsequently relax this assumption.

Figure 3 shows the constraints on the three models in the two-dimensional (Δ, w0) plane, including the reduced chi-square for the comparison of the models (with each choice of parameters) with our data sets. One sees that the data strongly prefers a present equation of state close to a cosmological constant, as well as a comparatively large transition width. In other words, any large deviations from wΛ = −1 are only allowed at large red-shifts. Thus, despite the fact that we have only used relatively low redshift data, we effectively constrain such deviations to happen only well into the matter era.

thumbnail Fig. 3.

Constraints in the (Δ, w0) plane for models with a late-time phase transition in the dark energy equation of state. The top, middle and bottom panels correspond, respectively, to Models B, C, and L described in the text. In all cases the Hubble constant has been kept fixed at H0 = 70 km s−1 Mpc−1, the black contours denote the one, two and three-sigma confidence levels, and the colourmap depicts the reduced chi-square of the fit for each set of model parameters (a dark red colour corresponds to a reduced chi-square of three or larger).

The one-dimensional posterior likelihoods for Δ and w0 are shown in Fig. 4, from which we can infer the dependence of the constraints on the choice of parametrisation for the transition. Broadly speaking, Model B generally leads to steeper transitions for a given choice of model parameters than Model C (cf. the examples shown in Fig. 1), and therefore the former model leads to more stringent constraints than the latter. The behaviour of Model L to some extent interpolates between the other models (for faster transitions it is closer to Model B, while for slower ones it is more similar to Model C) and therefore leads to intermediate constraints. The two-sigma (95.4% confidence level) constraints on the two parameters are summarised in Table 1.

thumbnail Fig. 4.

Posterior likelihoods for the transition width Δ and the present-day dark energy equation of state w0, marginalising the other, for models with a late-time phase transition in the dark energy equation of state. The blue dashed, red dotted, and black solid lines correspond, respectively, to Models B, C, and L described in the text. In all cases the Hubble constant has been kept fixed at H0 = 70 km s−1 Mpc−1.

Table 1.

Two-sigma (95.4% confidence level) constraints on the dark energy equation of state parameters w0 and Δ (marginalised over the other), for the various two-parameter dark energy models described in the text.

Study of the impact of the assumption of a fixed Hubble constant is relevant to the present work. We have done this for Model L, assuming a uniform prior between H0 = 65 km s−1 Mpc−1 and H0 = 75 km s−1 Mpc−1, which span the range of reasonable values given the currently available data (Planck Collaboration XIII 2016; Riess et al. 2016). The results are shown in Fig. 5, and also shown for comparison in the last line of Table 1. We see that this has a minimal effect on the constraints (making them weaker by a very small amount), the reason being that both data sets being used are fully compatible with the previously chosen value of H0 = 70 km s−1 Mpc−1.

thumbnail Fig. 5.

Constraints in the (Δ, w0), (H0, Δ) and (H0, w0) planes (top, middle and bottom panels, respectively) for Model L, with the Hubble constant being a free parameter and the third parameter marginalised. In all cases the black contours denote the one, two and three-sigma confidence levels, while the colourmap depicts the reduced chi-square of the fit for each set of model parameters (a dark red colour corresponds to a reduced chi-square of three or larger).

4. Three-parameter extension

Our analysis in the previous section does contain one caveat. Our choice of wi = 0 in the above analysis implies that at early times we have an additional matter-like component1. In other words, such models have a non-zero amount of early dark energy (Wetterich 2004), which is constrained by CMB experiments such as Planck to sub-percent level (Planck Collaboration XIV 2016), while the constraints in the previous section would allow values as high as about ten percent. In the current section we therefore relax this assumption, and repeat the analysis for the three-parameter extension of the L model, with wi as an additional free parameter.

In this case Eq. (11) becomes

w ( z ) = w i ( 1 + w 0 ) ( 1 + z ) 1 / Δ ( w i w 0 ) ( 1 + w 0 ) ( 1 + z ) 1 / Δ + ( w i w 0 ) ; $$ w{(z)}=\frac{w_\text{i}{(1+w_0){(1+z)}^{1/\mathrm\Delta}-{(w_\text{i}-w_0)}}}{(1+w_0){(1+z)}^{1/\mathrm\Delta}+{(w_\text{i}-w_0)}}; $$(15)

we note that if either w0 = −1 or wi = −1 this parametrisation reduces to w(z) = −1 throughout. The effective redshift is now

z t = ( w i w 0 1 + w 0 ) Δ 1 , $$ z_\text{t}={\left(\frac{w_\text{i}-w_0}{1+w_0}\right)}^\mathrm\Delta-1, $$(16)

while the Friedmann equation becomes

H 2 H 0 2 = Ω m ( 1 + z ) 3 + ( 1 Ω m ) [ 1 + w 0 1 + w i ( 1 + z ) 1 / Δ + w i w 0 1 + w i ] 3 Δ ( 1 + w i ) . $$ \frac{H^2}{H_0^2}={\mathrm\Omega}_\text{m}{(1+z)}^3+{(1-{\mathrm\Omega}_\text{m})}{\left[\frac{1+w_0}{1+w_\text{i}}{(1+z)}^{1/\mathrm\Delta}+\frac{w_\text{i}-w_0}{1+w_\text{i}}\right]}^{3\mathrm\Delta{(1+w_\text{i})}}. $$(17)

We can now repeat the analysis of the previous section for this extended model. We still fix the cosmological parameters Ωm = 0.3 and H0 = 70 km s−1 Mpc−1, and maintain our assumptions of a flat logarithmic prior on (1 + w0), log10(1 + w0) = [−3, −0.5], and a uniform prior on the transition width, Δ = ]0, 1]. As for the early-time value of the dark energy equation of state, wi we separately consider the cases of a flat prior, wi = ]−1, 0], or a logarithmic prior, log10(1 + wi) = [−3, 0], again as a means of checking how sensitive the results are to these choices.

The results of these analyses are presented in Figs. 6 and 7 and Table 2. As expected, the extended parameter space weakens the constraints on individual parameters, though this occurs in an interesting way. For the case of a flat (uniform) prior one can still obtain two-sigma bounds on w0 and Δ, though in Table 2 we report one-sigma bounds since wi is unconstrained at two-sigma (but still constrained at one-sigma). On the other hand, for a log prior wi and w0 are still constrained at one sigma (but not two), while Δ is totally unconstrained.

thumbnail Fig. 6.

Constraints in the (w0, wi) and (Δ, wi) planes (left and right panels, respectively) for the three-parameter Model L. The top panels correspond to a uniform prior on wi, while the bottom ones correspond to a uniform prior on log10(1 + wi). In all cases the black contours denote the one, two and three-sigma confidence levels, while the colourmap depicts the reduced chi-square of the fit for each set of model parameters (dark red corresponds to a reduced chi-square of 1.5 or larger).

thumbnail Fig. 7.

Posterior likelihoods for the transition width Δ and the values of the dark energy equation of state in the asymptotic past, wi and the present day, w0, with the other parameters marginalised, in the three-parameter version of the L model. The black solid and green dashed lines correspond respectively to the cases of logarithmic and uniform priors on wi, as detailed in the text; for comparison the likelihood for the case of a constant wi = 0 is also shown by the red dotted lines. In all cases the Hubble constant has been kept fixed at H0 = 70 km s−1 Mpc−1.

Table 2.

One-sigma (68.3% confidence level) constraints on the dark energy equation of state parameters wi, w0 and Δ (marginalised over the others), for the three-parameter dark energy L model described in the text, with a uniform prior on wi or on log10(1 + wi).

The physical reason for this behaviour is clear. In the two-parameter case, with a fixed value wi = 0 the data strongly prefers a dark energy equation of state at low redshifts (and in particular today, w0) close to that of a cosmological constant, and therefore the putative transition is pushed to high redshifts (into the matter era). In the three parameter case, and despite the correlations between the three parameters, both values of the equation of state (wi and w0) are constrained to be close to wΛ = −1 (though at a lower level of statistical significance), and in such a case there is effectively no transition, so the transition width Δ is unconstrained. As is clear from Fig. 7, the choice of a linear or logarithmic prior for wi does affect the posteriors, due to the effects of volume weighting.

The fact that in the case of a logarithmic prior the constraint on wi is slightly stronger than that of w0 is worthy of comment. A simple but physically meaningful way to classify dynamical dark energy models is to divide them into freezing and thawing ones, depending of the behaviour of the dark energy equation of state w(z) (Caldwell & Linder 2005): qualitatively, in the former ones w(z) is approaching −1 today, while in the latter it is moving away from it. Recent observational constraints such as those from the Planck mission in combination with other probes (Planck Collaboration XIV 2016) may be interpreted as slightly preferring thawing models, as pointed out in Linder (2015). In our case, models with wi < w0 are thawing models so – with the caveats of the choice of prior and the comparatively low level of statistical significance – our results are also consistent with this interpretation.

5. Discussion

It is enlightening to compare our results with those of previous authors, as a means to ascertain how the various data sets and modelling assumptions influence the resulting constraints. We briefly discuss other relevant works in the context of our results.

The work of Bassett et al. (2002) provided the initial motivation for our own work. Their analysis has three free parameters, the transition redshift zt, a future value of the dark energy equation of state which is approximately wf (and indeed, given their other choices of parameters it effectively coincides with w0), and the dark energy density ΩQ. They fixed wi = 0 (as did we in the first part of our article), but further restricted the transition width to obey Δ = zt/30. Moreover, they assume a flat universe (so Ωm + ΩQ = 1) and a Hubble constant H0 = 65 km s−1 Mpc−1.

Using a combination of CMB, Type Ia supernova and large scale structure data then available, Bassett et al. (2002) find (at the 68.3% confidence level) z t = 2.0 0.8 + 2.2 $ z_\text{t}=2.0_{-0.8}^{+2.2}$, wf < −0.8 and Ω Q = 0.73 0.04 + 0.02 $ {\mathrm\Omega}_Q=0.73_{-0.04}^{+0.02} $. It is clear (and discussed therein) that these constraints mainly come from the CMB: the acoustic peaks and integrated Sachs-Wolfe effect are both affected in these models. The relatively small Type Ia supernova data set then available carries a fairly small statistical weight. In any case, their data sets have a relatively small sensitivity to high redshifts, so to some extent their constraint on zt can be interpreted as a lower limit, which given their imposed relation between zt and Δ would be Δ > 0.04, which is comparable to but weaker than our two-parameter result.

A different kind of transition was studied by Lazkoz et al. (2016), who considered a class of so-called unified dark energy models, in which dark matter and dark energy are assumed to be different manifestations of the same underlying component which behaves as the former at high redshifts and as the latter at low redshifts. In Lazkoz et al. (2016) the transition is phenomenologically described by a Heaviside-like function, with the behaviour of the dark energy equation of state before the transition assumed to be wi = 0, and the low-redshift behaviour being w0 = −1. By using a combination of CMB, galaxy clustering and Type Ia supernova data, they find a best-fit value for the effective redshift of the transition zt = 4.7 ± 0.6 (at the 68.3% confidence level), though from the point of view of Bayesian evidence there is no clear preference for this model over standard ΛCDM.

Another four-parameter model for dynamical dark energy was recently proposed by Jaber & de la Macorra (2016). This is an extension of the two-parameter CPL model (Chevallier & Polarski 2001; Linder 2003) with two additional parameters describing the characteristic redshift of the transition between the early and late time behaviours, and the steepness of this transition. Constraints on this model, from a combination of CMB and BAO data together with the local value of the Hubble parameter (Riess et al. 2016) have been recently discussed in Jaber & de la Macorra (2018). The four parameters are allowed to vary with stated priors, for example on zt = [0, 3], w0 = [−1, 0] and wi = [−1, 0]. The authors then find w 0 = 0.96 0.17 + .022 $ w_0=-0.96_{-0.17}^{+.022}$, and w i = 0.00 0.02 + 0.04 $ w_\text{i}=0.00_{-0.02}^{+0.04} $ (despite their statement on the priors on w0 and wi), with a preferred transition redshift z t = 1.3 0.4 + 1.4 $ z_\text{t}=1.3_{-0.4}^{+1.4} $. Again there is no statistical preference of this model with respect to CDM. The fact that in their Table 2 they find a very tight constraint w i = 0.00 0.02 + 0.04 $ w_\text{i}=0.00_{-0.02}^{+0.04} $ from combining CMB, BAO and H0 is intriguing, especially because looking at their Fig. 3 the one-sigma confidence level includes wi = −1. The latter is also consistent with their statement that if they just combine BAO with CMB or BAO with H0 they obtain upper bounds consistent with wi = −1.

The full four-parameter model of Linder & Huterer (2005) has also been studied recently by Marcondes & Pan (2017; who in addition study two other parametrisations). Here the four parameters are allowed to vary with ample priors (including both canonical and phantom equations of state), and the data used includes CMB, BAO, type Ia supernovas, cosmic chronometers, and the local value of the Hubble parameter. They find tight constraints on the current dark energy equation of state (specifically w f = 1.04 0.36 + 0.10 $ w_\text{f}=-1.04_{-0.36}^{+0.10} $ at the 68.3% confidence level for the model of Linder & Huterer 2005), much weaker constraints on wi (though consistent with wi = −1) and an effectively unconstrained transition redshift. While their results are consistent with ours, they also confirm the expectation that current data (even including the CMB) is not powerful enough to strongly constrain a four-parameter dark energy model.

More recently Durrive et al. (2018) have used CMB, Type Ia supernova and BAO data to constrain several classes of canonical quintessence dark energy models. One of the models studied coincides with our Model L since, as first discussed in Linder (2008), this phenomenological parametrisation is a good approximation to the evolution of the dark energy equation of state in the case of a particular sub-class of quintessence models known as scaling freezing models (for which there is no exact analytic expression for the evolution of w(z)). However, they fix the width of the transition to be Δ = 1/3, on the grounds that this is the value for which this parametrisation best approximates the dynamics of the said scaling freezing models.

Under these assumptions, Durrive et al. (2018) obtain a lower bound zt > 8.1, at the 95.4% confidence level, and marginalising over all other parameters. No constraints are explicitly given for the other parameters, but it is clear from the text that a cosmological constant is fully consistent with the data. As in the previous case, this constraint is again dominated by the CMB data. While this constraint is stronger than the one obtained in the previous section, our work does show that current low redshift data is already stringent enough to effectively push any hypothetical transition into the matter era. We allow more freedom to the behaviour of the dark energy equation of state (that is, we do not fix Δ), although we are more restrictive in fixing the other cosmological parameters.

Finally, let us again compare our results with those of Di Valentino et al. (2018), who have explored and constrained a particular class of vacuum metamorphosis models – first discussed in Caldwell et al. (2006) – where both asymptotic values of the dark energy equation of state correspond to a cosmological constant (wi = wf = −1, in the notation we have defined above) while during the low-redshift transition the equation of state is in the phantom branch, w(z) < −1. For this reason their results are not directly comparable with ours, other than in the general sense that opposite classes of models are being explored. Nevertheless, we note that by using CMB data together with local measurements of the Hubble constant, the authors claim that this class of models can alleviate the H0 tension that seemingly exists for the ΛCDM model.

In their model the effective transition redshift is not an independent parameter but depends on the matter density and the ratio of a mass-scale to the Hubble constant. Because of this, Di Valentino et al. (2018) do not explicitly provide constraints on the transition redshift, but it is curious to note that for their best-fit model the transition redshift would be zt ~ 1.2, which is again in the matter era and comparable to our constraints.

6. Conclusions

We have studied classes of phenomenological but physically motivated two and three parameter models where the dark energy equation of state undergoes a (possibly) rapid transition at low redshifts. Such transitions might be associated with the onset of the acceleration phase itself. By using Type Ia supernova and Hubble parameter measurements, we have provided updated constraints on the redshift and steepness of these possible transitions as well as on the values of the dark energy equation of state at the present day and in the asymptotic past in these models.

Our results confirm and update those of previous works, showing that such transitions are tightly constrained. The present-day value of the dark energy equation of state, w0 is always constrained to be close to that of a cosmological constant, wΛ = −1. In the two-parameter case where wi is fixed at zero, the low-redshift data pushes the transition to high redhsifts. Effectively, this transition must happen within the matter era. Nevertheless, a value of wi = 0 is also constrained by other data, and especially by the CMB. Relaxing this assumption and allowing wi to be a third free model parameter, then it is also constrained to be close to wΛ = −1 so there is effectively no transition. Our constraints have a mild dependence of the specific parametrisation being used and on the choice of priors for the model parameters, but should be relatively insensitive to the choice of the Hubble parameter H0.

We conclude that the dark energy equation of state near the present day must be very similar to that of a cosmological constant, and any significant deviations from this behaviour can only occur in the deep matter era. In particular, at the level of the behaviour of the dark energy equation of state there is effectively no room for a phase transition associated with the onset of the acceleration phase. Therefore, if one wants to identify possible deviations from a cosmological constant, it is important to develop observational tools that are capable of probing this behaviour in the deep matter era. Two such promising examples are tests of the stability of fundamental couplings (Martins 2017) and direct measurements of the expansion of the universe (Liske et al. 2008; Martins et al. 2016). We leave discussion of the potential of these observables to further constrain these models to subsequent work.

Acknowledgments

We are grateful to Ana Catarina Leite, Ferran Delgà, Sara Raposeiro, and Sofia Cardoso for helpful discussions on the subject of this work. This work was financed by FEDER—Fundo Europeu de Desenvolvimento Regional funds through the COMPETE 2020—Operacional Programme for Competitiveness and Internationalisation (POCI), and by Portuguese funds through FCT—Fundação para a Ciência e a Tecnologia in the framework of the project POCI-01-0145-FEDER-028987. C.J.M. is supported by an FCT Research Professorship, contract reference IF/00064/2012, funded by FCT/MCTES (Portugal) and POPH/FSE (EC). M.P.C. acknowledges financial support from Programa Joves i Ciència, funded by Fundació Catalunya-La Pedrera (Spain).


1

We thank the anonymous referee for emphasising this point.

References

  1. Abbott, T. M. C., et al. (DES Collaboration) 2017, ArXiv e-prints [arXiv:1708.01530] [Google Scholar]
  2. Bassett, B. A., Kunz, M., Silk, J., & Ungarelli, C. 2002, MNRAS, 336, 1217 [NASA ADS] [CrossRef] [Google Scholar]
  3. Caldwell, R., & Linder, E. V. 2005, Phys. Rev. Lett., 95, 141301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  4. Caldwell, R. R., Komp, W., Parker, L., & Vanzella, D. A. T. 2006, Phys. Rev. D, 73, 023513 [NASA ADS] [CrossRef] [Google Scholar]
  5. Chevallier, M., & Polarski, D. 2001, Int. J. Mod. Phys. D, 10, 213 [NASA ADS] [CrossRef] [Google Scholar]
  6. Concas, A., Pozzetti, L., Moresco, M., & Cimatti, A. 2017, MNRAS, 468, 1747 [NASA ADS] [CrossRef] [Google Scholar]
  7. Copeland, E. J., Sami, M., & Tsujikawa, S. 2006, Int. J. Mod. Phys. D, 15, 1753 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  8. Corasaniti, P. S., & Copeland, E. J. 2003, Phys. Rev. D, 67, 063521 [NASA ADS] [CrossRef] [Google Scholar]
  9. Di Valentino, E., Linder, E. V., & Melchiorri, A. 2018, Phys. Rev. D, 97, 043528 [NASA ADS] [CrossRef] [Google Scholar]
  10. Durrive, J.-B., Ooba, J., Ichiki, K., & Sugiyama, N. 2018, Phys. Rev. D, 97, 043503 [NASA ADS] [CrossRef] [Google Scholar]
  11. Farooq, O., Madiyar, F. R., Crandall, S., & Ratra, B. 2017, ApJ, 835, 26 [NASA ADS] [CrossRef] [Google Scholar]
  12. Hill, C. T., Schramm, D. N., & Fry, J. N. 1989, Comm. Nucl. Part. Phys., 19, 25, [407(1988)] [Google Scholar]
  13. Huterer, D., & Shafer, D. L. 2018, Rept. Prog. Phys., 81, 016901 [NASA ADS] [CrossRef] [Google Scholar]
  14. Jaber, M., & de la Macorra, A. 2016, ArXiv e-prints [arXiv:1604.01442] [Google Scholar]
  15. Jaber, M., & de la Macorra, A. 2018, Astropart. Phys., 97, 130 [NASA ADS] [CrossRef] [Google Scholar]
  16. Jimenez, R., & Loeb, A. 2002, ApJ, 573, 37 [NASA ADS] [CrossRef] [Google Scholar]
  17. Jones, D. O., Scolnic, D. M., Riess, A. G., et al. 2018, ApJ, 857, 51 [NASA ADS] [CrossRef] [Google Scholar]
  18. Lazkoz, R., Leanizbarrutia, I., & Salzano, V. 2016, Phys. Rev. D, 93, 043537 [NASA ADS] [CrossRef] [Google Scholar]
  19. Linder, E. V. 2003, Phys. Rev. Lett., 90, 091301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  20. Linder, E. V. 2008, Gen. Rel. Grav., 40, 329 [NASA ADS] [CrossRef] [Google Scholar]
  21. Linder, E. V. 2015, Phys. Rev. D, 91, 063006 [NASA ADS] [CrossRef] [Google Scholar]
  22. Linder, E. V., & Huterer, D. 2005, Phys. Rev. D, 72, 043509 [NASA ADS] [CrossRef] [Google Scholar]
  23. Liske, J., Grazian, A., Vanzella, E., et al. 2008, MNRAS, 386, 1192 [NASA ADS] [CrossRef] [Google Scholar]
  24. Liu, G. C., Lu, Y. J., Xie, L. Z., Chen, X. L., & Zhao, Y. H. 2016, A&A, 585, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Lopez-Corredoira, M., & Vazdekis, A. 2018, A&A, 614, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Lopez-Corredoira, M., Vazdekis, A., Gutierrez, C. M., & Castro-Rodriguez, N. 2017, A&A, 600, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Marcondes, R. J. F., & Pan, S. 2017, ArXiv e-prints [arXiv:1711.06157] [Google Scholar]
  28. Martins, C. J. A. P. 2017, Rep. Prog. Phys., 80, 126902 [Google Scholar]
  29. Martins, C. J. A. P., Martinelli, M., Calabrese, E., & Ramos, M. P. L. P. 2016, Phys. Rev. D, 94, 043001 [NASA ADS] [CrossRef] [Google Scholar]
  30. Mortonson, M. J., Hu, W., & Huterer, D. 2009, Phys. Rev. D, 80, 067301 [NASA ADS] [CrossRef] [Google Scholar]
  31. Parker, L., & Raval, A. 1999, Phys. Rev. D, 60, 123502; Erratum: 2003, 67, 029902 [NASA ADS] [CrossRef] [Google Scholar]
  32. Parker, L., & Raval, A. 2000, Phys. Rev. D, 62, 083503; Erratum: 2003, 67, 029903 [NASA ADS] [CrossRef] [Google Scholar]
  33. Perlmutter, S., et al. 1999, ApJ, 517, 565 [NASA ADS] [CrossRef] [Google Scholar]
  34. Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Planck Collaboration XIV. 2016, A&A, 594, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Riess, A. G., Filippenko, A. V., Challis, P., et al. 1998, AJ, 116, 1009 [NASA ADS] [CrossRef] [Google Scholar]
  37. Riess, A. G., Macri, L. M., Hoffmann, S. L., et al. 2016, ApJ, 826, 56 [NASA ADS] [CrossRef] [Google Scholar]
  38. Scolnic, D. M., Jones, D. O., Rest, A., et al. 2018, ApJ, 859, 101 [NASA ADS] [CrossRef] [Google Scholar]
  39. Suzuki, N., Rubin, D., Lidman, C., et al. 2012, ApJ, 746, 85 [NASA ADS] [CrossRef] [Google Scholar]
  40. Wetterich, C. 2004, Phys. Lett. B, 594, 17 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Two-sigma (95.4% confidence level) constraints on the dark energy equation of state parameters w0 and Δ (marginalised over the other), for the various two-parameter dark energy models described in the text.

Table 2.

One-sigma (68.3% confidence level) constraints on the dark energy equation of state parameters wi, w0 and Δ (marginalised over the others), for the three-parameter dark energy L model described in the text, with a uniform prior on wi or on log10(1 + wi).

All Figures

thumbnail Fig. 1.

Examples of the redshift dependence of the dark energy equation of state in the three models considered in this work, for various choices of the free parameters in the models. Blue dashed, red dotted, and black solid lines correspond, respectively, to the Models B, C and L described in the text – cf. respectively, Eqs. (7), (9), and (11). The top left panel has w0 = −0.9 and Δ = 0.2, the top right one has w0 = −0.99 and Δ = 0.2, the bottom left one has w0 = −0.9 and Δ = 0.8, and the bottom right one has w0 = −0.99 and Δ = 0.8.

In the text
thumbnail Fig. 2.

Dependence of the transition redshift zt on the present-day value of the dark energy equation of state w0 and the transition width Δ in the three models considered in this work. The top, middle and bottom panels correspond, respectively, to the Models B, C, and L described in the text – cf. respectively, Eqs. (7)(12). In all panels the colourmap denotes the value of zt corresponding to each pair of values of w0 and Δ (we note that the colourmap and the ranges of w0 and Δ are different in each panel) and some specific contours are also identified by black lines with white labels.

In the text
thumbnail Fig. 3.

Constraints in the (Δ, w0) plane for models with a late-time phase transition in the dark energy equation of state. The top, middle and bottom panels correspond, respectively, to Models B, C, and L described in the text. In all cases the Hubble constant has been kept fixed at H0 = 70 km s−1 Mpc−1, the black contours denote the one, two and three-sigma confidence levels, and the colourmap depicts the reduced chi-square of the fit for each set of model parameters (a dark red colour corresponds to a reduced chi-square of three or larger).

In the text
thumbnail Fig. 4.

Posterior likelihoods for the transition width Δ and the present-day dark energy equation of state w0, marginalising the other, for models with a late-time phase transition in the dark energy equation of state. The blue dashed, red dotted, and black solid lines correspond, respectively, to Models B, C, and L described in the text. In all cases the Hubble constant has been kept fixed at H0 = 70 km s−1 Mpc−1.

In the text
thumbnail Fig. 5.

Constraints in the (Δ, w0), (H0, Δ) and (H0, w0) planes (top, middle and bottom panels, respectively) for Model L, with the Hubble constant being a free parameter and the third parameter marginalised. In all cases the black contours denote the one, two and three-sigma confidence levels, while the colourmap depicts the reduced chi-square of the fit for each set of model parameters (a dark red colour corresponds to a reduced chi-square of three or larger).

In the text
thumbnail Fig. 6.

Constraints in the (w0, wi) and (Δ, wi) planes (left and right panels, respectively) for the three-parameter Model L. The top panels correspond to a uniform prior on wi, while the bottom ones correspond to a uniform prior on log10(1 + wi). In all cases the black contours denote the one, two and three-sigma confidence levels, while the colourmap depicts the reduced chi-square of the fit for each set of model parameters (dark red corresponds to a reduced chi-square of 1.5 or larger).

In the text
thumbnail Fig. 7.

Posterior likelihoods for the transition width Δ and the values of the dark energy equation of state in the asymptotic past, wi and the present day, w0, with the other parameters marginalised, in the three-parameter version of the L model. The black solid and green dashed lines correspond respectively to the cases of logarithmic and uniform priors on wi, as detailed in the text; for comparison the likelihood for the case of a constant wi = 0 is also shown by the red dotted lines. In all cases the Hubble constant has been kept fixed at H0 = 70 km s−1 Mpc−1.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.