Cosmography and the redshift drift in Palatini $f({\cal R})$ theories

We present an application to cosmological models in $f({\cal R})$ theories within the Palatini formalism of a method that combines cosmography and the explicit form of the field equations in the calculation of the redshift drift. The method yields a sequence of constraint equations which lead to limits on the parameter space of a given $f({\cal R})$-model. Two particular families of $f({\cal R})$-cosmologies capable of describing the current dynamics of the universe are explored here: (i) power law theories of the type $f({\cal R})={\cal R}-\beta /{\cal R}^n$, and (ii) theories of the form $f({\cal R})={\cal R}+\alpha \ln{{\cal R}} -\beta$. The constraints on $(n,\beta)$ and $(\alpha,\beta)$, respectively, limit the values to intervals that are narrower than the ones previously obtained. As a byproduct, we show that when applied to General Relativity, the method yields values of the kinematic parameters with much smaller errors that those obtained directly from observations.

the usual Einstein-Hilbert Lagrangian density by an arbitrary function of the Ricci curvature scalar R [1,2,3,4,5]. Cosmological models built with them seek to describe current astronomical data without the use of the so-called dark energy, which in the scope of the standard cosmological model amounts up to ∼ 70% of the total matter content of the observable universe [6].
Within f (R)-theories the dynamics of the gravitational degrees of freedom is governed by the field equations derived from the minimisation of the corresponding gravitational action. Different variational principles (usually referred as formalisms) can be considered in order to get such equations, according to the role attributed to the connection Γ [1]. Among them, the Palatini formulation is based on the assumption that the metric and the connection are independent fields. In such a case, the corresponding Riemann and Ricci tensors are constructed with a connection a priori independent of the metric [7]. The Einstein-Hilbert Lagrangian is then replaced by a function f (R), where R is defined as R ≡ g µν R µν (Γ ), and R µν is the Ricci tensor defined in terms of the independent connection.
The dependence of the function f with the scalar curvature R introduces a set of constants, characterising a particular family of f (R)-cosmologies. The particular form of f as well as a range of suitable values for the constants must be chosen taking into account several criteria, such as the appropriate sequence of cosmological eras [8], the correct dynamics of cosmological density perturbations [9], as well as the correct weakfield limit at both the Newtonian and post-Newtonian levels [10], the well-posedness of the Cauchy problem [11], and the correct fit of cosmological observables [12,13,14]. 1 Among the observable quantities with the potential of discriminating between f (R)-models and those based on GR, particular attention was recently given to the time variation of the cosmological redshift z due to the variation of the expansion rate of the universe, namely the redshift drift (RD). The possibility of using this observable as a test of cosmological models was first proposed by Sandage [16], and later developed by other authors [17,18,19,20,21]. The redshift drift (a.k.a. Sandage-Loeb effect) was considered for many years of little use in the task of distinguishing cosmological models because of the difficulties associated to its measurement. However, this observable may be important for Cosmology, since it allows the test of the Copernican Principle [19], and has the potential to distinguish different cosmological models [20]. As first discussed in [22], the redshift drift can be also used to limit the values of the otherwise arbitrary constants of a given theory of gravity by resorting to its series expansion in powers of z. Such expansion can be computed using two different approaches. The first one is of a cosmographical type, i.e. independent of the dynamics of the subjacent theory and only based on the assumed symmetries of the space-time. The coefficients of the series expansion depend, in this case, on the so-called cosmographical kinematic parameters defined in terms of the time derivatives of the scale factor, the values of which follow from different observations. The second approach yields a series expansion which explicitly depends on the dynamics of the scale factor through the gravitational field equations of a given theory. The subsequent term-by-term comparison of the series expansions leads to a sequence of constraint equations, which explicitly relate different orders of derivatives of the scale factor (through the kinematic parameters) and the function f (R) and its derivatives (and, consequently, the constants of an specific family of f (R)-models) evaluated at z = 0. It is worth pointing out that such a comparison does not actually depend on the actual measurements of the redshift drift.
The constraint equations mentioned above can be used in two directions: (i) to get theoretical estimations of the kinematic parameters if General Relativity is assumed, and (ii) to constrain the space-parameter of particular f (R)-models if the kinematic parameters are derived from independent observational data. In this regard, cosmography has been widely used to distinguish between cosmological models (see for instance [23,24,25,26,27,28,29,30]). In particular, the application 1 The Dolgov-Kawasaki instability does not take place in in f (R) theories, see [15]. of the comparative method coming from the redshift drift series was used in [22] within the metric formalism to constrain the space parameter of f (R)-models. The main goal of this work is to extend the analysis of [22] to models built with f (R) theories in the Palatini formalism. We shall apply the results to the particular cases of power-law gravity (f (R) = R − β/R n ) and logarithmic gravity (f (R) = R + α ln (R) − β).
The paper is organised as follows. General considerations about the RD and the cosmographical approach from which constraints are obtained are presented in Section 2. The application of these ideas to the standard cosmological model are given in Section 3. In Section 4, we present the features of f (R)-theories within the Palatini formalism and apply the method to set limits on the parameter space of the two above-mentioned f (R) functions, and we compare our results with previously reported limits. Our final remarks are presented in Section 5.

The redshift drift
The redshift of a photon emitted by a source at time t and observed at time t 0 , is defined as follows: Due to the variation of the expansion rate of the universe, the redshift of a source is indeed a function of time. Then, a second photon emitted at t ′ = t + ∆t will have a redshift z(t ′ ). This time variation of the redshift is the so-called redshift drift, and can be expressed up to first order as [18] ∆z where ∆t 0 is the time delay between the two observed photons. Note that this observable depends neither on specific features of the source (such as its absolute luminosity) nor on the definition of a standard ruler. We shall present next the two approaches to the redshift drift, each of which involve a series expansion of this observable in terms of z. The first one is based on geometric and kinematic properties of the metric, namely a cosmographic treatment, whereas the other takes into account the dynamics imposed by a chosen theory of gravity through the cosmological field equations.
A cosmographic approach.The geometric and kinematic properties of the metric (14) are characterised by the socalled kinematic parameters, defined as the coefficients of the series expansion of the scale factor around t 0 . In particular, where a dot indicates derivatives w.r.t. time, and the sub-index 0 indicates that all quantities are evaluated at z = 0. The aim of the cosmographic approach in this setting is then to compute the redshift drift in terms of these quantities. 2 By performing a series expansion of H(z) in terms of the above defined kinematic parameters, the redshift drift can be written as [22] The coefficients of this equation, which follows from the cosmographic approach, depend purely on the properties of the FLRW metric, in the sense that no dynamical evolution for the scale factor was assumed.
A dynamical approach.We shall compute the dynamical counterpart of Eq. (3), now considering the dynamics obeyed by the scale factor a(t). Let us consider again the general expression for the redshift drift given by Eq.
(2). The dynamics enters through the Hubble parameter function, H(z), the evolution of which is completely determined by the gravitational field equations once a theory of gravity is chosen. Using the expansion of H(z) in a Taylor series of the redshift, together with the chain rule dH/dz = dH/dt · dt/dz, we obtain Note that Eq. (4) for the redshift drift involves the gravitational field equations through H and its derivatives.

Cosmography for the ΛCDM-model
Cosmography is a useful standpoint to study assumptions of cosmological models which are based entirely on the Cosmological Principle [33], and provides valuable model-independent information about the evolution of the scale factor and its derivatives (see also [34,35,36] for an extended discussion). This mathematical framework is inherently kinematic in the sense that it relies only on geometrical assumptions for the metric and is independent of the dynamics obeyed by the scale factor. We will focus in this section on the cosmographic information provided by the redshift drift in the case of the ΛCDM cosmology. The term-by-term comparison of both series for different powers of z leads to a sequence of equations (actually, an infinite number of them) for the kinematic parameters in terms of the cosmological parameters. The first members of the sequence are These relations allow to estimate values for the kinematic parameters as a function of the value of Ω m,0 . This estimation has the advantage of being independent of the restrictions associated to the convergence and the truncation of the Taylor series usually implemented in cosmography (see discussion in [37,25]). Taking for the dimensionless energy densities the values Ω m,0 = 0.315 ± 0.007 and Ω Λ,0 = 0.6847 ± 0.0073 [6], we obtain The above values are in agreement with other estimations of the kinematic parameters coming from different data sets [38,39,40], but present significantly smaller errors. This is particularly convenient for the case of higher-order parameters, such as j 0 , which have large errors when estimated by other methods.

The Palatini formalism
Let us consider a general action with the form where f (R) is an arbitrary function of the curvature The matter Lagrangian density L m depends on the matter fields ψ m , the metric g and its first derivatives, but does not depend on the affine connection Γ , which appears only in the gravitational action. The energy-momentum tensor is conserved since the total (gravitational plus matter) action is diffeomorphism-invariant, and gravity and matter are minimally coupled by assumption [41,42]. The Palatini f (R) gravity is then a metric theory (in the sense that the matter is minimally coupled to the metric and not coupled to any other fields), and hence the energy-momentum tensor T µν and its conservation laws will remain the ones of GR [1].
The variation of the action with respect to the metric and the connection yields, respectively, [1] where f ′ ≡ df /dR, and∇ ρ represents the derivative operator associated to the independent connection Γ ρ µν , which is assumed symmetric (torsion-less). We use units such c = 1 and κ 2 = 8πG.
The trace of Eq. (10) yields an algebraic equation for the Palatini Ricci scalar R in terms of the trace of the energy-momentum tensor T , The trace of Eq. (11), written as is used to define the connection of the metric h µν ≡ f ′ (R)g µν , conformal to g µν . Then,∇ ρ h µν = 0 and the connection is actually the Christoffel symbols of the conformal metric h µν . The General Relativity case is recovered when f ′ (R) = 1 (see Eq.13). Let us consider cosmological solutions by assuming the homogeneous and isotropic flat FLRW metric, described by the line element The modified Friedmann equation for a dust-dominated universe then becomes [8]

Constraints on particular f (R)-models
Two different approaches to compute a series expansion of the redshift drift were presented in the Section II, namely the cosmographic and dynamical treatments. As exemplified by the ΛCDM model, the term-by-term comparison of both series expansion leads to a sequence of constraint equations to be satisfied for a given cosmological model. In the case of f (R)-models these equations relate kinematic quantities with the parameters of the f (R) function. Therefore, in this section we will consider observational values for the kinematic parameters (already estimated in the literature from different data sets) to constrain the parameter space of the f (R) functions.
From the comparison of the linear terms in z in both expansion series (3) and (4), we obtain a relation of the form with G a lengthy algebraic function of its arguments 3 . This expression, together with the trace Eq. (12), directly implies a constraint on the space-parameter of a given f (R). Note that the absence of H 0 in the constraint equation 16 reduces the sources of error. We shall work with the values q 0 = −0.57 +0.10 −0.08 , estimated in [39] using the supernova type Ia JLA compilation [44], and Ω m,0 = 0.315 ± 0.007 from the last results of Planck Collaboration [6].
We present next the analysis of two f (R)-cosmological models: (i) power-law type models, characterized by the function f (R) = R − β/R n , and (ii) logarithmic type models of the form f (R) = R + α ln R − β. These particular choices of the function f (R) allow the occurrence of three cosmological phases (radiation-, matter-, and de Sitter-dominated eras) [8], and describe the accelerated expansion without the introduction of nonstandard sources of matter. Moreover, these theories satisfy the general criteria for cosmological viability discussed in Sect. 1. 4 It is important to remark that both models have been shown to accommodate a phase of late-time accelerated expansion also by the integration of the corresponding equations of motion, and for a range of parameters compatible with our findings. Model (i) presents a transition from non-accelerated to accelerated expansion at around z = 1, as shown in [46,47] by means of the evolution of the effective equation of state parameter w ef f in terms of z, and also in [48,47] by studying the behaviour of q(z). The same can be said about model (ii), for which the function w ef f (z) was analyzed in [46,47], and the deceleration factor q(z) was presented in [47]. 3 Additional constraints can be obtained by taking into account higher-order terms of both series, which involve higher order derivatives of the scale factor (that is, additional kinematic parameters). See for instance [43]. 4 In addition to these models, the cosmology of the socalled exponential gravity, described by the function f (R) = R − αnH 2 0 1 − e −R/(αH 2 0 ) , has been widely studied within the Palatini formalism [45]. However we have verified that this type of models does not satisfy the necessary condition imposed by the correct sequence of cosmological eras.

Power law gravity: f (R) = R − β/R n
Theoretical and observational studies determining the viability of power-law type gravity within the Palatini formalism have been presented by many authors [49,50,51,52,53]. We restrict our analysis here to the function f (R) = R − β/R n , where n is a dimensionless parameter and β is reported in units of H 2(n−1) 0 . Using dynamical analysis, it was shown that this type of theories can reproduce the sequence of radiation-dominated, matterdominated and de-Sitter eras for n > −1 and β > 0 [8]. Concerning the cosmology described by these models, a study of the evolution of the density perturbations has been presented in [9]. Numerical studies of the dynamic of flat models were also presented in [54].
In this work we present the constraint relation for the (n, β) space, which must be satisfied by the parameters to be consistent with the cosmographic and dynamical approaches of the redshift drift. The restrictive condition coming from Eq. (16) is shown in Fig. 1, together with the best-fit values previously reported and their associated error bars. The shadowed regions indicate the error propagation associated to the Ω m,0 and q 0 values. Since the reported error for the q 0 estimation is large, we also include a forecast corresponding to a measurement of q 0 with an error of 10% of its value. The plot shows that mildly improved estimations of the deceleration parameter would narrow even more the allowed region on the parameter space. Our constraint relation is consistent with the ΛCDM model, which is recovered for (n, β) = (0, 4.1082 ± 0.0438) according to the latest CMB data [6], and is in good agreement with other limits for this model. Furthermore, note that the error associated to our constraint is less than that of all but one of the best-fit values, and one of these is discarded by the constraint.   Fig. 1 Constraint on the parameter space (n, β) for theories of the type f (R) = R − β/R n coming form the cosmographic approach to the redshift drift, together with the best-fit values previously reported and the ΛCDM case. The shadowed regions indicate the error propagation associated to the Ω m,0 and q 0 values, and a forecast corresponding to a measurement of q 0 with an error of 10% of its value.

Logarithmic gravity
This type of theories within the Palatini formalism was firstly studied in [60], where it was shown that logarithmic terms in the action may be of use in the description of the accelerated expansion and reduce to the standard Friedmann evolution for high redshifts. Furthermore, it was shown that logarithmic gravity is one of the simplest forms capable of reproducing a suitable sequence of cosmological eras [8]. Limits on the parameter space (α, β), where α and β are computed in units of H 2 0 , have been less studied than that in the previous case. Estimations using different sets of data coming from SN, BAO and CMB observations were obtained in [8], with best-fit values (α, β) = (0.11 +1.75 −1.11 , 4.62 +3.54 −5.58 ). A more recent work has reported (α, β) = (−0.48 +1.26 −2.67 , 3.58) [46] (see also [47]).
The restriction given by Eq. (16) for the case of logarithmic gravity is shown in Fig. 2. As in the previous analysis, the shadowed regions indicate the error propagation associated to the Ω m,0 and q 0 parameters, together with a forecast corresponding to a measurement of q 0 with an error of 10% of its value. In this case our constraint significantly improves the limits in the parameter space and is also consistent with the ΛCDM model.

Conclusion
Alternatives to General Relativity may be the key to describe the accelerated expansion of the universe without the use of the so-called dark energy. In particular,   Fig. 2 Constraint on the parameter space (α, β) for theories of the type f (R) = R + α ln R − β coming from the cosmographic approach to the redshift drift, together with the best-fit values previously reported and the ΛCDM case. The shadowed region indicates the error propagation associated to the Ω m,0 and q 0 values, and a forecast corresponding to a measurement of q 0 with an error of 10% of its value.
cosmological models based on f (R)-theories of gravity within the Palatini formalism have been constructed to deal with observational data sets and solar-system constraints, yielding a satisfactory description of the late time dynamics of the universe. The investigation of cosmological observables capable to distinguish between different models can also be used as a tool to set constraints on the parameters of such theories. In this direction, a novel method for f (R)-theories within the metric formalism was recently proposed in [22]. Its essence is the comparison of two Taylor expansions of a given observable. While the first one is of the cosmographic type, namely based on the assumed symmetries of the space-time and independent of the dynamics obeyed by the scale factor, the second takes into account the specific dependence of H with z, determined by the field equations of a given gravity theory. The order-by-order comparison of these two series leads to a sequence of constraints relating the parameters of the theory. Furthermore, the method does not rely on the actual measurement of the observable, but on the condition that both series coincide term-by-term. This last feature can be taken as an advantage for exploring limits imposed by yet to-be-measured quantities as is the case of the redshift drift. Using values for the cosmological parameters Ω m,0 and Ω Λ,0 obtained from observations, we have shown that a direct application of this method to GR yields theoretical estimations of the kinematic parameters with an error much smaller than the one that follows straight from the cosmographic approach.
The main goal of this work was to apply the method to f (R)-cosmological models within the Palatini for-malism. Using two series expansions of the redshift drift we obtained a constraint relation, given by Eq.(16), on the parameter space of such models in terms of the so-called kinematic parameters. Two particular f (R)models were studied: (i) those based on power-law type theories, characterised by a Lagrangian function of the form f (R) = R − β/R n , and (ii) logarithmic-gravity type models of the form f (R) = R + α ln R − β. These two particular models allow the occurrence of the right sequence of cosmological eras, and have been previously studied in the literature. We have used here the redshift drift expansions to impose an independent constraint relation on their parameters in terms of q 0 and Ω m,0 . Our results, presented in Figs. 1 and 2, are in good agreement with previously reported values, but with considerably smaller errors. These could be reduced even more taking into account future improved estimations of the kinematic parameters, as shown by the forecast obtained using an error of ten percent of the value of q 0 .
To close, we would like to emphasise that the bounds obtained by the method presented here do not depend on the actual measurement of the redshift drift. Such bounds can also be considered together with those coming from theoretical considerations (such as the energy conditions [43]), in order to decide whether a given f (R) theory furnishes an appropriate description of the current state of the universe.