Forecast and analysis of the cosmological redshift drift.

The cosmological redshift drift could lead to the next step in high-precision cosmic geometric observations, becoming a direct and irrefutable test for cosmic acceleration. In order to test the viability and possible properties of this effect, also called Sandage-Loeb (SL) test, we generate a model-independent mock data set in order to compare its constraining power with that of the future mock data sets of Type Ia Supernovae (SNe) and Baryon Acoustic Oscillations (BAO). The performance of those data sets is analyzed by testing several cosmological models with the Markov chain Monte Carlo (MCMC) method, both independently as well as combining all data sets. Final results show that, in general, SL data sets allow for remarkable constraints on the matter density parameter today [Formula: see text] on every tested model, showing also a great complementarity with SNe and BAO data regarding dark energy parameters.

Introduction

Within the general relativity (GR) framework no reliable explanation to the current acceleration of the universe exists which is simpler than a -term or cosmological constant [1]. It behaves as a fluid with negative pressure [2], thus driving gravitational repulsion. This is of course also the kind of behavior displayed by the plethora of other possible fluids so far proposed to try to accommodate data better than a cosmological constant. In broad terms, these settings, causing the universe to accelerate, are usually included in the so-called dark energy theories (for reviews, see [3-7]). There are other theoretical routes with different levels of complexity (not necessarily unrelated [8]) which venture to modify GR.

The background expansion of the universe can be measured with a lot of different probes: luminosity distances from Type Ia Supernovae (SNe) [9-12]; the acoustic peaks in the Cosmic Microwave Background (CMB) [13, 14]; their counterpart imprinted in clustered matter, i.e. Baryon Acoustic Oscillations (BAO) [15-19]; through the matter power spectrum obtained from weak lensing [20, 21]. Usually a time integral along redshift connects those data with the expansion rate/history of the universe, the Hubble parameter H(z). It is enough to constrain the geometry and energy content of the universe quite satisfactorily.

On the other hand it is expected that the expansion of the universe will make the redshift of a given astrophysical object exhibit a drift over time, which should in principle be amenable to giving an accurate description of that very same expansion once an underlying model is chosen. While looking for a possible temporal variation of the redshift of extra-galactic sources, Sandage came in 1962 [22] to the conclusion that it should indeed occur. However, the limited technological resources on deck at that epoch lead to the inference that a measurement time interval of the order of years would be required for a signal detection. When new spectroscopic techniques became available to astrophysicists, Loeb considered the concept [23] anew in 1998, he concluded that the new technology would allow a reduction in the observation time interval of a few decades. This cosmological redshift drift measurement, also called the Sandage–Loeb (SL) test, would then provide a direct proof of the accelerated expansion of the universe. In fact, this temporal variation is directly related to the expansion rate at the source redshift, being thus a direct measurement for the Hubble function.

The last results of the Planck survey [13] have made us enter an ultra high-precision cosmology era. Other future surveys are scheduled which should further improve the accuracy of cosmological measurements, for example Euclid [24], Wide-Field Infrared Survey Telescope (W-First) [25] or Square Kilometer Array (SKA) [26].

Thus, in the near future available resources will allow us to start thinking about the next level of cosmological observational data to which the cosmological redshift drift will contribute, complementing the previously cited surveys. However, even with future precision radio telescopes, the measurement of the SL effect represents a difficult enterprise [27] as it demands several years of observation (usually some decades) to register enough signal-to-noise ratio so as to yield a possible reliable detection of the cosmological redshift drift signal. Best candidate objects for a feasible detection of this faint signal are good Hubble flow tracers as far as possible [28]. As put forward by Loeb [23], an auspicious target would be the Lyman- forest measurements of distant quasars (QSO). With the spectroscopic techniques that are operational in the near future, like the CODEX (COsmic Dynamics and EXo-earth experiment) experiment [29] proposed for the European-Extremely Large Telescope (E-ELT) or radio telescopes as SKA [30], these observations will grant access to direct measurements of the Hubble parameter up to redshifts , a so far not yet observed redshift range. Thus, the SL test will open a new “cosmological window”.

Due to near future possibilities to measure the cosmological redshift drift, this type of observations recently has drawn some attention. The reconstruction of the theoretical SL signal that different cosmological models would produce has been explored quite extensively [31, 32]. It turns out that the range and variety of the different cosmological redshift drift signals created by various models is remarkable: from those created by different proposals for dark energy’s equation of state or modified gravity [33, 34] to the ones created by backreaction in an inhomogeneous universe without the presence of dark energy [35]; from the peculiar signals for Lemaître–Tolman–Bondi models [36, 37] to even a null signal [38] for the Universe, or other several exotic scenarios [39-44]. SL signals have been used as a hypothetical geometric cosmic discriminant [45-47] to show the corresponding improvement in the constraints that can be achieved due to the degeneracy breaking (around of improvements for dark energy parameters and even for matter density). SL mock data sets have been applied with similar results as cosmic observational discriminators to test other various models, like interactive dark energy models [48, 49], modified gravity [50, 51] and other exotic cosmologies [52, 53]. Their power to differentiate models has been exploited also in the context of the model-independent approach of cosmography [54, 55]. Besides, some new approaches [56] can lead to ambitious ideas, such as real-time cosmology [57].

We stress again the fact that the measurement of the cosmological redshift drift is not an easy pursuit. It requires a lot of planning due to the large observation time interval of the survey. Thus, foreseeing the contribution and behavior of this type of measurements is important we precisely carry out a quite thorough forecast analysis of cosmological redshift drift constraints on various cosmological models. The analysis includes a comparison between the proposed SL data with other future planned surveys, generating mock data based on the given specifications. Furthermore, unlike previous work, all mock data sets are generated in a fully model-independent way, without choosing a fiducial cosmological model to generate the points. In Sect. 2 the mathematical formalism of the cosmological redshift drift is introduced and the details of the mock data sets we use for our predictions are given. We find it convenient to produce a SL data set but also use auxiliary SNe and BAO data. In Sect. 3, we explain our MCMC procedure which will eventually constrain the cosmological models we have chosen as reference. Finally, in Sect. 4, the outcomes of the statistical analysis are presented and discussed. Then the main conclusions are summarized.

Cosmological redshift drift

A preliminary straightforward calculation introduces the main observable quantity we are going to focus on, i.e. the cosmological redshift drift, (see for example [29] or [33]). In a homogeneous and isotropic universe with a Friedmann–Robertson–Walker metric a source at rest emitting electromagnetic waves isotropically, without any (significant) peculiar velocity, is considered. Thus, the comoving distance between the source and an observer can be considered fixed. If the source emits electromagnetic waves during time and these are detected by the observer in the interval , where is the emission time and is the time they reach the observer, the following relation is satisfiedprovided the universe through which the waves travel is a spatially flat Friedmann–Robertson–Walker spacetime. If the time intervals are small the above expression leads to the well-known redshift relation between the emitted and the observed radiationwhere is the redshift of the source as measured at a certain observation time . Other waves can be emitted by the source time later, specifically at time . They will be observed at . Concerning these waves Eq. (2) has to be modified regarding the new time periods and redshift. Thus, the observer can measure the difference between the redshifts observed at and :Within the approximation, the first ratio can be expanded to linear orderInserting Eq. 2 into the first order expansion in Eq. 4 an approximated expression for the redshift variation can be obtained,Under the assumption that the observation time is today, we normalize by letting the corresponding scale factor satisfy . Then, using both the Friedmann equation and the known redshift equation Eq. (2), the above expression can be rewritten in terms of the Hubble parameter with being the Hubble constant today. This redshift variation can be expressed as a spectroscopic velocity shift . Using the dimensionless Hubble parameter , the final expression can be found:

Sandage–Loeb mock data set

In order to generate our SL observational mock data set in a fully model-independent manner, we try to derive a Hubble function from a phenomenological distance modulus, in a fashion similar to [58]. We propose this observable because it is well measured by Type Ia Supernovae (SNe) and can be extended to high redshifts, even if with lower precision, by Gamma Ray Bursts (GRBs, Mayflower sample) [59]. We model this phenomenological distance modulus thus:where is an ad hoc proposed function (among many) mimicking the luminosity distance. This phenomenological function is then fitted using the SNe data set Union 2.1 [11] for the low-redshift regime and the GRBs sample calibrated by the Padé method [59] for the high-redshift one. Once is fitted, other observational quantities relevant to our work can easily be obtained. For instance, the Hubble function can be derived recalling the relationwhere, in the spatially flat universe we are considering, the dimensionless luminosity distance is defined asand stores all the information related to the constants involved, such as the speed of light c, the Hubble constant and the SNe absolute magnitude. By comparing both distance moduli, from Eq. (9) and from Eq. (8), one realizes that the dimensionless luminosity distance is equivalent to the function . Thus, the dimensionless Hubble function isOnce such a phenomenological dimensionless Hubble parameter is obtained, we can “mimick” all the cosmological probes we need for our analysis as they are all related to it. In this way cosmological-model-independent mock data sets can be created where the only intrinsic information that is used for is that it has to be able to fit present data (in this case, SNe and GRBs). Of course, some arbitrariness lies behind the choice of the phenomenological function ; we have tried to use the most general type of functions possible and we have selected the best one based on a simple best-fitting (minimum ) criterion. The best performing function we have found iswhere the values for the parameters are shown in Table 1. It can be seen in Fig. 1, in the top left panel, that this function fits the distance modulus points of the Union 2.1 [11] and Mayflower [59] data sets as much satisfactorily as a CDM with Planck values, (sixth column of Table 4 in [13]). In the top right panel, we also compare the expansion rate function H which can be derived from Eq. (12) with the same Planck CDM and with data from cosmic chronometers [60]. In the bottom left panel, the comparison between angular diameter distance derived from Eq. (12) and the same Planck CDM is done, with the data coming as comoving angular diameter distance from galaxy clustering (BAO+FS column of Table 7 in [61]) and physical angular diameter distance coming from quasar cross-correlation (Eq. (21) in [18]). Finally, in the bottom right panel, it can be seen that the difference between our model and the Planck CDM is minimal for the case of the distance modulus and small for both the Hubble function and the angular diameter distance , all over the redshift range we cover with our mock data in our analysis.

Table 1

Parameter values of

	Estimate	Standard error
a	43.2025	0.146659
b	2.29876	1.60875
c	0.92048	0.969826
d	1.05317	0.62311
e	0.814751	0.922533

Fig. 1

Top left panel: Comparison between the selected phenomenological function given in Eq. (12) (solid red) with the Planck CDM (dashed blue) described in the text. Gray dots and bars are distance modulus values and related errors for SNe Union and Mayflower GRBs samples and black ones are for our generated SNe mock data. Top right panel: comparison between the H(z) function derived from Eq. (12) (solid red) with that corresponding to the Planck CDM (dashed blue) described in the text. Gray dots and bars are expansion rate values and related errors from cosmic chronometers and black ones our generated mock data. Bottom left panel: comparison of the function derived from Eq. (12) (solid red) with that corresponding to the Planck CDM (dashed blue) case described in the text. Gray dots and bars are angular diameter distances values and related errors from BOSS and SDSS and black ones our generated mock data. Bottom right panel: relative residuals between our model and the Planck CDM for the Hubble function (dashed blue), the distance modulus (solid black) and the angular diameter distance (dotted red)

Table 4

Parameter results of the slow-roll model

Data set	h	Ω_m	Ω_b	δw	χmin2	χred2
BAO	0.678-0.006+0.006	0.341-0.008+0.009	0.0485-0.0009+0.0010	0.074-0.035+0.036	5.20	0.520
SNe	0.675-0.006+0.006	0.330-0.012+0.010		-0.260-0.138+0.128	1383.10	0.509
SNe + BAO	0.688-0.006+0.005	0.324-0.006+0.006	0.0468-0.0008+0.0008	0.005-0.030+0.030	1402.05	0.513
SL (24 years)	0.674-0.006+0.006	0.324-0.007+0.007	0.0471-0.0008+0.0008	0.231-0.539+0.252	13.12	0.504
SNe + BAO + SL (24 years)	0.687-0.005+0.005	0.325-0.002+0.002		0.014-0.029+0.028	1417.62	0.513
SL (28 years)	0.674-0.006+0.006	0.323-0.006+0.006	0.0472-0.0008+0.0008	0.282-0.330+0.199	17.85	0.686
SNe + BAO + SL (28 years)	0.686-0.005+0.005	0.325-0.002+0.002		0.017-0.027+0.028	1423.21	0.515
SL (32 years)	0.674-0.006+0.006	0.323-0.005+0.005	0.0473-0.0007+0.0008	0.306-0.254+0.171	23.31	0.897
SNe + BAO + SL (32 years)	0.686-0.005+0.005	0.325-0.002+0.002		0.021-0.029+0.029	1429.59	0.518

Once we have our , we only need to specify a fiducial value for the Hubble constant to insert in Eq. (7), whose effect is only the rescaling of the velocity shift value. We fix the value of from the TT, TE, EE + lowP + lensing baseline model of Planck [13]. Then, concerning the SL data, the points lie in the redshift range , randomly distributed within the following bins: (13 points), (7 points), (4 points), (3 points) and (3 points). In this way we try to mimic the reduction of the number of data points while increasing the redshift as in [62].

According to Monte Carlo simulations carried out to eventually mimic results from CODEX [29, 63], the standard deviation on the measured spectroscopic velocity shift can be estimatedwhere x is 1.7 for , and 0.9 beyond that redshift, S / N is the spectral signal-to-noise ratio of Ly-, is the number of observed quasars and their redshift. The error for the mock data is given by assuming a fixed number of integration time hours which yields a value of for the signal-to-noise ratio and for the number of quasars observed [28]. We also introduce some noise to disperse the data points around the fiducial value derived from , using a Gaussian distribution centered on such values, and with a standard deviation corresponding to the expected error on the SL observation, , obtained by error propagation from the fitted parameters of the selected function Eq. (12).

Note that the magnitude of the observed cosmological redshift drift is proportional to the observation period although the error does not depend on it. Thus, once a data set for some given observational time is created, any new mock data set with different observation period can easily be calculated byWe use the three observation periods of 24, 28 and 32 years, which are the most illustrative among the data sets tested. The resulting data sets for SL test can be seen on Fig. 2.

Fig. 2

Datasets for SL test based on for different observation periods: blue (circle and dashed line) for 24 years, red (triangle and dotted line) for 28 years and green (square and solid line) for 32 years

Auxiliary mock datasets

Additional future mock data sets are included alongside the cosmological redshift drift data set in order to better constrain models. Basically, the reason why we introduce in the picture these other probes is our interest in studying and quantifying the relative performance of SL with respect to more standard and used probes and our aim of finding out whether the cosmological redshift drift data have some degree of complementarity with them, thus providing eventual tighter constraints. These auxiliary mock data sets are created from the same model-independent function of Eq. (8).

W-First SNe

The first mock data set we produce is a SNe catalog based on the W-First forecast [25] which includes 2725 SNe randomly picked in redshift bins of spread through a redshift range of according to the distribution given by [25].

Given that in the SNe case one measures the distance modulus, direct use can be made of the fitted function Eq. (8) to generate the mock data points. As in the SL case some Gaussian noise is introduced to disperse the data points around the mean value.

To create the error bars for this catalog and the dispersion for the Gaussian noise the information given in [25] is used. The statistical errors they account for are the following: the photometric measurement error, ; the intrinsic luminosity dispersion ; and the gravitational lensing magnification . Besides, they assume a systematic error . Thus the total error per SNe iswhere and is the number of SNe in the bin. The data set generated for the W-First SNe survey is shown on Fig. 1.

Euclid BAO

The second data set considered is BAO. We choose the future Euclid survey [24] as the experiment to reproduce. The two quantities considered are the angular diameter distancenormalized by the sound horizon , and the Hubble parameter times the sound horizon, , where the value of Mpc, consistent with the previous , is used [13].

Both the angular diameter distance and the Hubble parameter are reconstructed, again using Eq. (8). It is already discussed that the Hubble parameter can be inferred as in Eq. (11), once a value for is decided upon. Instead, in order to derive the angular diameter distance from Eq. (8), its definition and its relation with the luminosity distance are used, thus leading toThe redshift values of the data set are taken from [64]. They specifically are the central redshifts of 15 bins with width, spread from to . The error in each redshift value for both and is build from the percentage error also given in [64]. Finally, some Gaussian noise is introduced using the error from each bin as dispersion when generating the points and . The resulting data sets can be seen in Fig. 1 before normalizing the observables by the comoving sound horizon .

Testing models

Within the Bayesian framework, we aim to find how SL constrains the probability distribution function of some cosmological parameters. For that purpose the posterior distribution is needed, or equivalently the likelihood, which can be straightforwardly computed with MCMC sampling while minimizing the function. The knowledge of the posterior probability gives a better and more complete information as regards the parameters, including the full correlation among them.

Thus, once we have the mock data sets, we build the function for each observable and, once all contributions are summed, we minimize the total in order to perform our statistical analysis. The contribution for the spectroscopic velocity shift is simplywhere follows from Eq. (7). The errors are given by Eq. (13). The errors are arranged into a diagonal covariance matrix. Depending on whether the SL surveys will use overlapping redshift bins or not, the error could be more realistically given by a non-diagonal covariance matrix. As we lack such information, we adopt the optimistic diagonal covariance matrix assumption, always keeping in mind that it could lead to a general underestimation of the global errors on the cosmological parameters. The period of observation as specified above changes depending on the mock SL survey tested. In the case of the contribution of SNe we havewhere the error is given by Eq. (15). We can marginalize over the parameter by expanding the in Eq. (19) with respect to aswhereThen, integrating out of the likelihood we can retrievewhere has now no dependence on the parameter. We have to point out that also in this case we are using a diagonal covariance matrix because it is not possible to forecast out-of-diagonal terms. This may lead to underestimated errors on the cosmological parameters. With BAO we have two correlated measurements to contribute to the total ; these are and . With the Hubble parameter from our phenomenological fit and the angular diameter defined in the previous section, the comoving sound horizon readswhere the sound speed is , with and [65]. The comoving sound horizon is evaluated at photon-decoupling epoch redshift given by the fitting formula [66]with where and are the baryon and matter content of the universe and . The BAO contribution is independently calculated for each redshift, . However, taking into account the correlation of the magnitudes, each term at each redshift has the formwhere and are the differences between the model predicted and the mock generated measurements, The correlation between the two magnitudes and in each redshift is fixed as [67]. Since CMB data are not used, SNe data are marginalized over the parameter and BAO data do not give information about it (because and do not basically depend on it), the parameters and the combination cannot be well constrained. Thus, we also include a Gaussian prior for and , with and for both derived from Planck [13].

The minimization of the function was performed using the MCMC method [68-70], with a Wolfram Mathematica self-developed code based on the Metropolis–Hastings algorithm. In order to see the contribution of each mock data set to the total we also have run chains for each data set separately. In this way, the cosmological redshift drift data are compared with those from the other future surveys. Thus it is found whether it will be useful and up to what extent. Moreover, for a round analysis regarding the viability of the Sandage–Loeb test and the performance of the future (mock) surveys several dark energy scenarios are put to the test.

CDM

The first model we test is the extremely well-known CDM model [71, 72], which has no degree of freedom in the dark energy equation of state and whose dimensionless Hubble parameter is given bytaking with [73]and using [65]. We enforce , and as physical priors. The same is done for all the other models analyzed in this paper. The results of the Bayesian analysis for the CDM model can be seen in Table 2 and Fig. 4.

Table 2

Parameter results of the CDM model

Data set	h	Ω_m	Ω_b	χmin2	χred2
BAO	0.689-0.002+0.002	0.335-0.008+0.008	0.0467-0.0004+0.0004	9.47	0.861
SNe	0.675-0.007+0.006	0.301-0.008+0.008		1387.41	0.510
SNe + BAO	0.689-0.002+0.002	0.324-0.006+0.006	0.0472-0.0004+0.0004	1402.07	0.513
SL (24 years)	0.674-0.006+0.006	0.328-0.003+0.003	0.0467-0.0004+0.0004	14.51	0.538
SNe + BAO + SL (24 years)	0.689-0.002+0.002	0.324-0.002+0.003		1417.82	0.513
SL (28 years)	0.673-0.006+0.006	0.328-0.003+0.003	0.0467-0.0004+0.0004	19.73	0.731
SNe + BAO + SL (28 years)	0.689-0.002+0.002	0.325-0.002+0.002		1423.49	0.515
SL (32 years)	0.673-0.006+0.006	0.328-0.002+0.002	0.0468-0.0004+0.0004	25.75	0.954
SNe + BAO + SL (32 years)	0.689-0.002+0.002	0.325-0.002+0.002		1430.03	0.518

Fig. 4

CDM model; solid contours limit regions and clear contours region. Purple for the BAO, green for SNe and red for SL data set, gray for SNe + BAO and blue SL + BAO + SNe. First set (left) for 24 years, second (middle) for 28 and third (right) for 32 years

Quiessence

The second model tested is quiessence [74, 75] with a single degree of freedom in the dark energy equation of state parameter (i.e. no redshift dependence). Its dimensionless Hubble parameter is given bywhere all the parameters except w are built like in the CDM model and have the same priors. The parameter w has the prior . This range was chosen after having verified that expanding it further has no influence on results. Table 3 and Fig. 5 show the results for quintessence model.

Table 3

Parameter results of the quiessence model

Data set	h	Ω_m	Ω_b	w	χmin2	χred2
BAO	0.677-0.006+0.006	0.336-0.008+0.008	0.0485-0.0009+0.0010	-0.948-0.025+0.024	5.10	0.510
SNe	0.675-0.007+0.006	0.341-0.015+0.013		-1.244-0.122+0.123	1383.13	0.509
SNe + BAO	0.686-0.006+0.006	0.323-0.006+0.006	0.0472-0.0008+0.0008	-0.987-0.023+0.022	1401.74	0.513
SL (24 years)	0.674-0.007+0.006	0.323-0.011+0.011	0.0474-0.0008+0.0008	-0.888-0.660+0.141	12.73	0.490
SNe + BAO + SL (24 years)	0.685-0.005+0.005	0.324-0.002+0.003		-0.982-0.021+0.020	1417.07	0.513
SL (28 years)	0.674-0.007+0.006	0.321-0.009+0.009	0.0475-0.0008+0.0008	-0.845-0.176+0.102	17.32	0.666
SNe + BAO + SL (28 years)	0.684-0.005+0.005	0.325-0.002+0.002		-0.979-0.021+0.020	1422.49	0.515
SL (32 years)	0.674-0.006+0.007	0.320-0.008+0.007	0.0476-0.0008+0.0008	-0.830-0.119+0.084	22.62	0.870
SNe + BAO + SL (32 years)	0.684-0.005+0.005	0.325-0.002+0.002		-0.977-0.020+0.020	1428.75	0.517

Fig. 5

Quiessence model; solid contours limit regions and clear contours region. Purple for the BAO, green for SNe and red for SL data set, gray for SNe + BAO and blue SL + BAO + SNe. First set (left) for 24 years, second (middle) for 28 and third (right) for 32 years

Slow-roll dark energy

We consider another one-parameter dark energy model, coming from the slow-roll dark energy scenario described in [76]. Its dimensionless Hubble parameter, taking into account a radiation component [77, 78], is given byFor we impose a prior of the same width as that of the parameter w of quiessence. However, as is supposed to have its mean value at , its prior is designed accordingly. Thus, we take . The results for the slow-roll dark energy model can be found in Table 4 and Fig. 6.

Fig. 6

Slow-roll model; solid contours limit regions and clear contours region. Purple for the BAO, green for SNe and red for SL data set, gray for SNe + BAO and blue SL + BAO + SNe. First set (left) for 24 years, second (middle) for 28 and third (right) for 32 years

CPL

We are also interested in testing models of dark energy whose equation of state parameter w has more than one degree of freedom. As our first two-parameter dark energy model we take the CPL model [79, 80], its dimensionless Hubble parameter beingwhere all the terms except and are built like in previous models and with the same priors. The parameter has the same prior as w does in quiessence; and we take for the second parameter. We demand in this case in order to have an equation of state for the DE component which is negative in the asymptotic past. Table 5 and Fig. 7 give the results of our Bayesian analysis for the CPL model.

Table 5

Parameter results of the CPL model

Data set	h	Ω_m	Ω_b	w_0	w_a	χmin2	χred2
BAO	0.677-0.007+0.006	0.395-0.075+0.035	0.0486-0.0010+0.0010	-0.558-0.467+0.311	-1.722-1.457+2.035	3.37	0.374
SNe	0.675-0.006+0.006	0.356-0.038+0.030		-1.165-0.182+0.203	-0.555-1.708+1.232	1383.07	0.509
SNe + BAO	0.678-0.006+0.006	0.281-0.011+0.011	0.0484-0.0009+0.0010	-1.187-0.032+0.040	1.022-0.195+0.119	1386.45	0.508
SL (24 years)	0.674-0.006+0.007	0.328-0.017+0.008	0.0476-0.0008+0.0009	-1.026-1.518+0.498	-0.001-1.584+1.044	12.19	0.488
SNe + BAO+SL (24 years)	0.684-0.006+0.006	0.314-0.008+0.006		-1.117-0.066+0.066	0.602-0.286+0.286	1412.27	0.512
SL (28 years)	0.674-0.006+0.006	0.324-0.017+0.010	0.0476-0.0008+0.0009	-0.937-0.627+0.423	-0.021-1.541+0.869	16.56	0.662
SNe + BAO + SL (28 years)	0.683-0.006+0.006	0.315-0.008+0.006		-1.113-0.068+0.066	0.590-0.276+0.293	1417.52	0.513
SL (32 years)	0.674-0.006+0.006	0.319-0.019+0.010	0.0477-0.0008+0.0009	-0.895-0.256+0.372	0.240-1.147+0.637	21.65	0.866
SNe + BAO + SL (32 years)	0.683-0.006+0.006	0.315-0.008+0.006		-1.120-0.062+0.066	0.628-0.278+0.297	1423.36	0.516

Fig. 7

CPL model; solid contours limit regions and clear contours region. Purple for the BAO, green for SNe and red for SL data set, gray for SNe + BAO and blue SL + BAO + SNe. First set (left) for 24 years, second (middle) for 28 and third (right) for 32 years

Lazkoz–Sendra pivotal dark energy

Another model with two parameters for the equation of state for DE [81] is considered which can be understood easily as a perturbative departure from CDM up to second order in redshift. Even though it is a different parametrization as compared to CPL it can be also expressed in terms of the parameters and with the same interpretation: is the value of equation of state of the dark energy at present, whereas is its value in the asymptotic past. Specifically, the Lazkoz–Sendra pivotal dark energy parametrization has the following dimensionless Hubble parameter:where all the relative densities are built like in the CPL case, all parameters also having the same priors as in CPL, including and . In case of the Lazkoz–Sendra model the results of the Bayesian analysis are shown in Table 6 and Fig. 8.

Table 6

Parameter results of the Lazkoz–Sendra pivotal model

Data set	h	Ω_m	Ω_b	w_0	w_a	χmin2	χred2
BAO	0.677-0.006+0.006	0.391-0.065+0.034	0.0485-0.0009+0.0009	-0.661-0.333+0.244	-2.114-1.820+2.401	3.33	0.370
SNe	0.675-0.007+0.006	0.361-0.031+0.023		-1.170-0.143+0.150	-1.063-2.004+1.565	1383.13	0.509
SNe + BAO	0.680-0.006+0.006	0.295-0.008+0.009	0.0481-0.0009+0.0009	-1.093-0.025+0.028	0.934-0.196+0.112	1388.46	0.508
SL (24 years)	0.673-0.006+0.007	0.331-0.014+0.005	0.0477-0.0008+0.0008	-1.056-2.390+0.510	-0.616-2.186+1.416	12.53	0.501
SNe + BAO + SL (24 years)	0.683-0.006+0.006	0.311-0.006+0.006		-1.088-0.029+0.037	0.826-0.278+0.185	1408.54	0.510
SL (28 years)	0.674-0.007+0.007	0.325-0.015+0.009	0.0477-0.0008+0.0009	-0.884-0.449+0.349	-0.173-1.775+0.917	16.77	0.671
SNe + BAO + SL (28 years)	0.682-0.006+0.005	0.312-0.006+0.006		-1.089-0.029+0.037	0.839-0.279+0.177	1413.75	0.512
SL (32 years)	0.674-0.006+0.006	0.322-0.015+0.009	0.0479-0.0008+0.0008	-0.848-0.227+0.335	-0.080-1.498+0.815	21.93	0.877
SNe + BAO + SL (32 years)	0.682-0.005+0.005	0.312-0.006+0.005		-1.087-0.028+0.038	0.835-0.274+0.177	1419.53	0.514

Fig. 8

Lazkoz–Sendra pivotal model; solid contours limit regions and clear contours region. Purple for the BAO, green for SNe and red for SL data set, gray for SNe + BAO and blue SL + BAO + SNe. First set (left) for 24 years, second (middle) for 28 and third (right) for 32 years

Results and conclusions

In the summary tables for each model are found, the minimum value of is presented and the constraints for all the free parameters and the reduced . As explained in previous section the -minimization is done using different combinations of data sets. In the tables we first show the results from using BAO and SNe separately and those from joining both. We then move on to present the results from SL only. Finally, the results for the total SNe + BAO + SL combination are presented. When using SL data, each data set with different observation years is treated separately. In this way the performance of the cosmological redshift drift data sets can be clearly analyzed. For each model we also show the confidence contours for the most interesting cosmological parameters. Each MCMC round is tested for statistical convergence using the method described in [82].

In the CDM scenario we find that the cosmological redshift drift data provide remarkably good constraints on : when those data are used alone we get standard deviations on which are 2–3 times smaller than those from the SNe+BAO combination. Considering the broad priors taken for in all cases and the negligible correlation between the Hubble constant h and ,1 we conclude that the result for the matter density is not influenced by any prior and is solely given by the data.

Indeed, the SL data sets always do better in constraining than the SNe data and, depending on the model and on the years of observation, even better than the BAO data set. Once we combine the SL data set with the other two, the cosmological redshift drift is still helpful, even though the BAO + SNe data set already greatly improves the constraints in the parameter space. In general, it is clear that the cosmological redshift drift data considerably helps to constrain the parameter in all the models.

Regarding the dark energy parameters, it can be observed that for most cases the 24 years of observation for SL is not enough to properly constrain them. This can for example clearly be seen from the contours of the parameters and in Figs. 7 and 8. With 28 years SL data, the regions improve noticeable and with 32 years of observation both and regions are well constrained for all the DE parameters. The best example is seen in Fig. 8. However, similar behavior can be appreciated in the rest of the models. Besides, it is clear that increasing the observation years improves the overall constraining ability of the cosmological redshift drift data sets. It is worth to note that in all these cases the contours of the SL data set are almost perpendicular to the contours of the SNe and BAO data sets, thus showing a great complementarity between SL and the rest of the data sets [83], as for example in the –w plane for the quiessence model Fig. 5 or for the dynamical dark energy models Figs. 6, 7 and 8. This is very important because it means that even the cosmological redshift drift data set with the lowest observation period noticeably contributes to improved dark energy insights when used as cosmological probe together with other kind of observations.

However, if one focuses on the and parameters two things can be noted: first, that the best fit for the SNe + BAO case is completely different from the values derived from the SNe and the BAO separate analysis (this is more evident for than ); second, the errors on the and parameters slightly increase when adding cosmological redshift drift data to the SNe+BAO data. Both trends might have an explanation. Concerning the first one, if we look at the top panel of Fig. 3 (this is true for the CPL case but also for the LS model), we can see how unsatisfactorily the SNe and BAO contours overlap: the borders of the confidence levels show a small overlap in a region which is far from the best fit expected for each of them when considered separately. This reduces the constraints on the parameters in a considerable way and shifts the best fit estimations (not only in but also on ). Note also that this behavior is somehow expected and might be counter-productive in the future, as explained in [84]. Anyway, we must also remember we are working with mock data and not real data. Thus the potential future goodness of the joint use of SNe and BAO at present and maybe in the near future is not put at stake. Moreover, we have to remember that in order to gain more insights into a dynamical dark energy model we need to improve the number and the quality of data at high redshift; that is the reason behind pushing SNe observations to higher redshifts [85] for example or employing BAO data at . But the strongest hints about the dynamical nature of the dark energy might come from data like SL which are able to cover a larger and deeper redshift range. The second issue discussed above should be exactly connected to this: if we check again the bottom panel of Fig. 3 we can see how the SL data set alone, which should be more sensitive to a dynamical dark energy, determines a consistent shift in the parameter with respect to the SNe + BAO case but with smaller uncertainty with respect to both SNe and BAO data separately, which eventually ends in a slightly large error for this parameter for the total SNe + BAO + SL sample.

Fig. 3

Contours in the – plane for CPL; solid contours are for regions and clear contours are for regions. Top panel: purple is for the BAO data; green for SNe and gray for SNe + BAO. Bottom panel: red is for 32 years SL data; gray for SNe + BAO; blue for SNe + BAO + SL

In the case of models with a single DE parameter whose equation of state is fixed during time high-redshift SL data are also helpful. In the extreme case when SL data are added to the SNe+BAO data set, even the SL data with lowest observational period help constrain the single parameter of DE. However, it is also remarkable how every data set separately constrains the single DE parameter to a different value. Taking into account that the redshift range of each data set is quite different, the fact that they separately have a different value for the parameter could be evidence for a time evolution in the equation of state of DE. This is a clear example of another application for the SL observation, where its high-redshift data could easily test the time evolution of the equation of state of DE once compared to the results of other data sets coming from different sources.

Much of what has been stated above can be easily inferred upon closer examination of the various contours plots. However, these plots are more useful for analyzing the correlation between different parameters. As stated previously, in most of the contour plots a different correlation angle can be seen for the cosmological redshift drift data compared to the other data sets. Thus, it clearly emerges that SL data sets will be of utmost importance in breaking degeneracies among cosmological parameters. Besides, considering the high-redshift data that will be available thanks to cosmological redshift drift we conclude that it can be a cosmic observable worthy to consider.