Cosmological curvature acceleration

Observed accelerated behavior of the Hubble flow, usually related to the existence of a mysterious unknown fluid dubbed dark energy, can be framed in natural extensions of General Relativity that work at infrared scales. This approach gives rise to a geometric view of the phenomenon that does not require any new particle to be addressed. In this review paper, without claiming for completeness, we will give the main ingredients of such an approach considering models like f(R) gravity, with R the Ricci curvature scalar, and, in general, higher-order theories of gravity including other curvature invariants. In this context, we will review also Teleparallel Equivalent General Relativity and its extensions like f(T) gravity, where T is the torsion scalar. We will show that the curvature picture and the torsion picture can be related each other giving an equivalent paradigm where the observed cosmological acceleration can be addressed by extending the gravitational sector.


Introduction
Explaining the observed accelerated phase of the late Universe [1,2] is till now an argument of debate. According to the present view, approximately 70% of the energy content of the Universe is a mysterious form of dark energy, exotic, invisible, and unclustered. Many classes of models have been proposed to address the problem. They can be classified as follows: (i) models with cosmological constant Λ; (ii) dark energy related to some fundamental field; and (iii) possible modifications of gravity. The cosmological constant Λ gives rise to the cosmic acceleration and dominate today the Universe. It seems to be the most obvious explanation. However, the cosmological constant is affected by the well known cosmological constant problem and the 8coincidence problem [3].
The second class of models postulates the existence of a dark energy fluid (the equation of state is P ≈ −ρ (ρ and P are the energy density and pressure of the fluid, respectively), and it is assumed to dominate the late Universe. Such a fluid should be related to some fundamental field. Dark energy could even be phantom energy (the equation of state is P < −ρ). Many dark energy models have been proposed. They are not definitely convincing, as well as they are not free of finetuning problems. A third possibility consists in modia e-mail: capozziello@na.infn.it b e-mail: lambiase@sa.infn.it (corresponding author) fying gravity at the large scales, i.e. at infrared scales. Here we focus on this class of models.
The most popular and simplest theory proposed to modified General Relativity (GR) is the f (R) gravity. More than a modification, it is an extension of GR where the Einstein theory is a particular case. In fact, the usual Einstein-Hilbert action is modified to where f (R) is a non-linear function of its argument [4,5]. Here R denoted the Ricci curvature of the metric tensor g αβ with metric determinant g, κ ≡ 8πG, is the Newton coupling [6]. It is worth noticing that, in principle, the metric tensor contains several degrees of freedom: tensor, vector, and scalar, massless or massive.
In GR, only massless spin 2 graviton propagates, while its modifications give rise to the occurrence of other degrees of freedom. For example, the change R → f (R), in the action, induces, in addition to the massless graviton, the appearance of a massive scalar mode which can be responsible of the cosmic acceleration (this is analogous to the inflaton field driving the accelerated expansion of the early Universe). Including in the gravitational Lagrangian quadratic terms depending on Ricci and Riemann tensor, and/or other curvature invariants, can induce massive gravitons and vector degrees of freedom enriching a lot the dynamics [7,8].
A comment is in order at this point: in modifying gravity, it may occur that there could be violations of experimental limits on the parametrized-post-Newtonian (PPN) parameters at terrestrial and Solar System scales [9], instabilities, ghosts and, as in any newly proposed theory, the Cauchy problem could be ill-posed. All these problems are studied in literature, and we will briefly focus on them, being mainly interested on the cosmological consequences.
The prototype of f (R) gravity working at late epoch is where μ is a mass scale of the order of the present value of the Hubble parameter μ ∼ H 0 ∼ 10 −33 eV. Although ruled out by its weak-field limit [10] and by a violent instability [11], this model gives the idea underlying modified gravity: the 1/R correction is negligible in comparison with R at the high curvatures of the early universe, and kicks in only as R → 0, late in the history of the universe. Early and late cosmic acceleration can be successfully described by assuming this geometric approach, allowing at the same time to address shortcomings of the standard cosmological model [12][13][14][15][16][17]. The paradigm is that R 2 term dominates at early epoch, R −1 dominate at late cosmological scales while R is working at local scales.
The main advantage is that the left-hand side of the Einstein field equations get modified and the introduction of ad hoc matter contributions to the energymomentum is avoided [18][19][20][21][22][23][24]. This 'geometric picture', generally referred as Extended Theories of Gravity (ETG) [23], applies to all approaches in which curvature or torsion invariants enter into the effective gravitational action allowing the recovery of GR in some limit or as a particular case [23,25].
The general view is that terms depending on geometric invariant are significative at a given epoch [24], being directly related to important phenomena in the Universe evolution (the process of structure formation, early and late-time cosmic acceleration, inflation, and so on. Refs.  and references therein give an idea of possibilities.
From a fundamental physics point of view, such curvature invariants are derived as quantum corrections from renormalisation of gravity in curved spacetime [65]. Other curvature invariants, as R μν R μν , R μνσρ R μνσρ have also been taken into account in literature [66][67][68][69][70][71][72], in particular the Gauss-Bonnet topological invariant G which is related to the emergence of the trace anomaly in curved spacetime [65,73]. In general, if both R and G are present in the gravitational action, all the curvature budget is considered, if we do not take into account further derivative terms like R and others [8]. As shown in [74], the presence of nonlinear G terms in the action gives rise to a further inflationary episode that can be connected to the observed large-scale structure process. In these models, one has an early Gdominated phase followed by the usual R-dominated phase. In general, the presence of Gauss-Bonnet topological invariant can solve some problems of f (R) gravity, as discussed in [75][76][77][78][79][80][81][82][83]. The advantage of considering a Gauss-Bonnet non-minimal coupling is due to the fact that an improved phantom-quintessence phase of the late universe occurs thanks to such a term. The Gauss-Bonnet curvature becomes dominant and then the phantom phase is a transient: this means that the Big Rip singularity is avoided.
An important issue to be mentioned is related to the gravitational wave event GW170817 [84]. In fact, many alternative theories of gravity can be discarded owing the the upper bound on the wave propagation giving by |c g /c − 1| ≤ 5 · 10 −16 [85]. Interestingly f (R) gravity remains a viable theory, while GR modified with a scalar field f (φ)G [86], where f (φ) is a function of a phantom field and G is the Gauss-Bonnet invariant seems to be excluded by the observations. However, not only f (R) remains a viable theory, but remarkably, F (R, G) can be retained because it is a singularity free theory where the Gauss-Bonnet contribution enhances the reliable behavior of curvature quintessence [87] and inflation [74].
In this short review paper, without claiming for completeness, we restrict ourselves to curvature based theories of gravity in their various formulations, and we focus the possibility that they can address the accelerated behavior of the present observed Universe. In other context, these curvature based theories can also be used as an alternative to dark matter in galaxies and clusters of galaxies [89][90][91][92][93][94][95][96].
We adopt a conservative point of view considering these theories to test the principle that modifying gravity is a viable approach to explain the observed accelerated behavior of the cosmic Hubble fluid. However, none can claim that any modified gravity model, proposed thus far, is the correct final one, or has exceptional support from the observational data. While it is true that many models, in particular f (R), pass all the available experimental tests and fit the cosmological data, the same is true for many dark energy models, and it is currently impossible to use observational data to discriminate between most of them, and between dark energy and modified gravity models.
The layout of the paper is the following. Different formulations of f (R) gravity, like metric, Palatini and metric-affine ones, are discussed in Sect. 2. Viability criteria are reported in Sect. 3 and equivalent between f (R) gravity and scalar-tensor theories are sketched in Sect. 4. Cosmology in curvature based modified gravity is discussed in Sect. 5. In particular, we take into account f (R) gravity and F (R, G) gravity improved with Gauss-Bonnet invariant. This last case is important because the introduction of the topological invariant exhausts all the quadratic curvature invariants that can be considered into the gravitational action. We close the section considering the f (T, B) extension of teleparallel gravity, where T is the torsion scalar and B is the boundary term. As we will see, a linear combination of T and B is capable of restoring the Ricci curvature scalar giving rise to an equivalent picture of f (R) gravity. Conclusions and outlooks are reported in Sect. 6.

Different formulations of f (R) gravity
As discussed, f (R) gravity is the straightforward extension of GR. It has quite a long history (see [97] for an historical review). In his theory, Weyl (1919) added a term quadratic in the Weyl tensor to the Einstein-Hilbert Lagrangian [98]. Later, during the the 1960's and 1970's, it was realized that quadratic corrections to S EH were necessary to improve the renormalizability of GR [99], and, some years later, Starobinsky pointed out that quadratic corrections are needed to fuel inflation without scalar fields [18]. Non-linear corrections are also motivated by string theories [100,101].
As proposed in literature, f (R) gravity comes in three versions: (1) metric (or second order) formalism; (2) Palatini (or first order) formalism; and (3) metric-affine formalism. Here, we will sketch them and their relations with dark energy.

Metric formalism
In the metric formalism [4,5], the action is Variation with respect to the metric tensor g αβ yields the field equations where a prime denotes differentiation with respect to R. The first two terms on the right hand side introduce fourth order derivatives of the metric, hence the name "fourth order gravity" sometimes given to these theories. The trace of Eq. (5) yields where T ≡ T α α is the trace of the matter stress-energy tensor. This second order differential equation for f (R) differs deeply from the trace of the Einstein equations R = −κ T which, instead, relates algebraically the Ricci scalar to T . We already see that f (R) is indeed a dynamical variable, the scalar degree of freedom contained in the theory.
Formally, one can rewrite the field Eq. (5) in the form of effective Einstein equations considering as an effective stress-energy tensor containing geometric terms. Of course, as usual when adopting this procedure, T αβ satisfies generalized energy conditions [102] and the effective energy density is, in general, nonpositively defined. As is clear from these equations, in f (R) gravity one can define an effective gravitational coupling G eff ≡ G/f (R) in a way analogous to what is done in scalar-tensor theories. Hence, f (R) must be positive in order for the graviton to carry positive kinetic energy. It is possible to demonstrate that (7) is a perfect fluid whose origin is geometric [103].

The Palatini formalism
In the Palatini approach, both the metric g αβ and the connection Γ α βγ are independent variables, i.e., the connection is not the metric connection of g αβ . While in GR, the metric and Palatini variations produce the same (Einstein) equations, this is no longer true for nonlinear Lagrangians. 1 After metric f (R) theories were proposed as alternatives to dark energy, also the Palatini version has been considered for the same purpose, originally for the model f (R) = R − μ 4 /R [105]. The Palatini action reads In this approach, we have to take into account two Ricci tensors: R αβ , which is constructed using the Levi-Civita connection of the metric g αβ , andR αβ which is the Ricci tensor of the non-metric connection Γ α βγ . The latter gives rise to the scalarR ≡ g αβR αβ . The matter part of the action does not depend explicitly on the connection Γ, but only on metric and matter fields, collectively denoted by ψ (m) . The variation of the Palatini action (8) yields the field equations derivatives of f , in contrast with Eq. (5). Variation with respect to the independent connection produces the field equatioñ where∇ δ denotes the covariant derivative associated to the non-metric connection Γ. Traces of Eqs. (9) and (10) yield Equation (12) tells us that∇ γ is the covariant derivative of the metricg conformally related to g αβ . In this sense, the Palatini approach can be considered a bimetric formulation. Note that Eq. (11) is an algebraic equation for f (R): hence, this quantity is non-dynamical, contrary to f (R) gravity in metric formalism. This lack of dynamics has important consequences which we will sketch below. However, metric and Palatini formalism can be connected each other taking into account conserved currents [106]. An important remark is in order at this point. It is worth noticing that, in Palatini formalism, non-linear couplings between matter and gravity naturally emerge as it is easy to see from Eq. (11). Such an equation, sometimes called the structural equation or the master equation [107], points out that Palatini field equations contain nonlinear terms in the stress-tensor and its derivatives. If we solve algebraically Eq. (11), according to the form of f (R), the solution R is a function of the stress-energy tensor trace T . As a consequence, the conformal factor in Eq. (13) is a function of T . Using this result in the field Eq. (10), these equations contain nonlinear terms in T and its derivatives. This picture immediately tells us that the the theory reduces to GR in vacuum (T = 0) and the so called universality of the Einstein equations is restored [108]. In this perspective, the theory gives cosmic accelerated expansion only if some effective cosmological constant is recovered, at late times, or strong deviations from GR emerge into dynamics at early times [50].

The metric-affine formalism
In metric-affine formalism [109], also the matter component depends on the connection. It is Here Γ can be also non-symmetric. In this case, torsion associated with matter has to be considered. The approach was originally introduced to consider the spin of elementary particles which generate the torsion. A detailed discussion of f (R) gravity with torsion in metric-affine formalism is reported in [110]. In the following, we shall focus on metric and Palatini f (R) gravity in view of cosmological applications.

Viability criteria for f (R) gravity
Besides cosmological applications, f (R) gravity and any modification/extension of GR must also pass tests imposed by Solar System and terrestrial experiments. In other words, any modification of GR must satisfy some viability criteria to be a self-consistent formulation for a relativistic theory of gravitation. In summary they are: • having the correct Newtonian and post-Newtonian limit; • cosmological dynamics has to be reproduced; • not suffer for instabilities and ghosts; • giving rise to cosmological perturbations compatible with data from the cosmic microwave background and large scale structure; • having a well-posed and well-formulated Cauchy problem.
Failing also even a single one of these criteria is a statement that the theory does not work. These viability criteria are examined in the following.

Weak field limit and solar system tests
In order to assess whether the limits set on the PPN parameter γ by the available Solar System experiments hold or not, one needs to find the weak-field solution of the field equations and compute this parameter. It is possible to consider a static, spherically symmetric, non-compact body which constitutes a perturbation of a background de Sitter solution. The line element is written as in Schwarzschild coordinates, where dΩ 2 is the line element on the unit 2-sphere and Ψ and Φ are post-Newtonian potentials. These are of small amplitude, |Ψ(r)| , |Φ(r)| << 1, and one considers small (noncosmological) scales so that H 0 r << 1 while expanding the Ricci scalar around the constant curvature of the background de Sitter space, that is R(r) The PPN parameter γ is given by γ = −Φ(r)/Ψ(r) [9]. Three assumptions can be made [10]: where m is the effective mass of the scalar degree of freedom of the theory, i.e., it is assumed that this scalar field is light and has a range larger than the size of the Solar System (we remind the reader that there are no experimental constraints of scalars in the range m −1 < 0.2 mm). 3. For matter composing the spherical body, the pressure is negligible, P 0, so that T = T 0 + T 1 −ρ.
The first and the last assumption are not stringent, but the second one is, as it will be clear below.
The trace Eq. (6) yields the equation for the Ricci scalar perturbation where is the effective mass squared of the scalar. Equation (17) coincides with the expression obtained in the gaugeinvariant stability analysis of de Sitter space and in propagator calculations.
If mr << 1, the solution of the linearized field equations is Therefore, the PPN parameter is in gross violation of the experimental limits [150] This result would be the end of metric f (R) gravity if the assumptions made in the calculation were satisfied. However, this is not the case for assumption 2): mr is not always less than unity due to the Chameleon Effect . This consists in the effective mass m depending on the curvature or, alternatively, the matter density of the environment. The scalar degree of freedom can be short-ranged (say m > 10 −3 eV, corresponding to a range λ < 0.2 mm) at Solar System densities and evade the experimental constraints, while being longranged at cosmological densities and thus it is able to affect the cosmological dynamics [147,151]. Although at a first glance the Chameleon Effect could be seen as a contrived and fine-tuned mechanism, f (R) gravity is rather complicated and the effective range does indeed depend on the environment. The Chameleon Mechanism is well-known and accepted in quintessence models, in which it was discovered for the scalar field potential V (φ) ≈ 1/φ [152,153]. Many forms of the function f (R) are known to exhibit the chameleon mechanism and pass the observational tests [154]. For example, the model is compatible with the PPN limits if μ ∼ 10 −50 eV∼ 10 −17 H 0 [151]. It is obvious that a correction term ∼ R n with n < 1 to the Einstein-Hilbert Lagrangian R will come to dominate as R → 0 + (for example, for n = 1/2, √ R > R as R → 0). The model (22) is compatible with the experimental data but it could be essentially indistinguishable from a dark energy model. Hope of discriminating between dark energy and f (R) models, or between different modified gravities, relies on the study of the growth of cosmological perturbations.

Cosmological dynamics
In order to be acceptable and compatible with observations, a cosmological model has to exhibit some minimal requirements which are: (i) an early inflationary epoch or an alternative mechanism to solve horizon, flatness, and monopole problems together with a mechanism to generate density perturbations; (ii) a following radiation-dominated era and (iii) a matter-dominated era, and, finally, (iv) the model has to be compatible with the present acceleration.
f (R) gravity has mainly be proposed to address this sequence in a self-consistent way. The future eras are usually found to be either an eternal de Sitter attractor phase, or a Big Rip singularity truncating the history of the universe at a finite time.
In particular, smooth transitions between different eras are required. It has been pointed out that the exit from the radiation era, in particular, may have problems in many models [111][112][113][114][115], a warning that has to be considered in building up f (R) cosmologies. In general, exit from the radiation, or any era can be achieved for given f (R) models. Consider, for example, the so-called "designer f (R) gravity": one can prescribe a desired expansion history of the universe by a choice of the scale factor a(t) and then integrate a differential equation determining the function f (R) producing this a(t) [116][117][118][119][120]. In general, this function is not unique and can assume complicated forms and not only the usual R − μ/R n or other simple forms. In summary, the recipe is a sort of "inverse scattering" where the form of the model is achieved starting from the desired behavior needed for matching the observations.

Instabilities and ghosts
The prototype model for discussing instabilities is again f (R) = R − μ 4 /R with μ ∼ H 0 ∼ 10 −33 eV. Shortly after it was proposed, this model was found to a suffer catastrophic instability [11]. The stability analysis was later generalized to any metric f (R) theory [121] and the extension to even more general gravitational theories has been pursued in [122][123][124]. The approach consists in parameterizing the deviations from GR as where is a small positive constant with the dimensions of a squared mass and the function ϕ is dimensionless. The trace equation for the Ricci scalar R takes the form Next, one expands around a de Sitter background and writes the metric locally as while the scalar degree of freedom R is expanded as with R 1 a perturbation. To first order, the trace equation yields the dynamical equation for R 1 The last term on the left hand side is dominated by the term in −1 and gives the effective mass squared of R 1 from which one deduces that the theory is The case of GR is excluded by the assumption f = 0, but one can extend the stability criterion for metric f (R) gravity to be f ≥ 0, including GR.
As an example, the prototype model f (R) = R − μ 4 /R, which has f < 0 is unstable. The time scale for the onset of this instability is dictated by the smallness of the scale μ which corresponds to ∼ 10 −26 s [11], making this an explosive instability.
One can give a physical interpretation of this result as follows [125]: remembering that the effective gravitational coupling is G eff = G N /f (R), with G N the Newton constant, if dG eff /dR = −f G N /(f ) 2 > 0 (which corresponds to f < 0), then G eff increases with R and a large curvature causes gravity to become stronger, which, in turn, causes a larger R, in a positive feedback mechanism driving the system away. If instead dG eff /dR < 0, then a negative feedback damps the increase in the gravitational coupling strength.
Palatini f (R) gravity, by contrast, is described by second order field equations, the trace equation f (R)R− 2f (R) = κ T . It is not a differential equation but rather a non-dynamical algebraic one and, therefore, there is no instability [126].
The previous analysis for metric f (R) gravity, obtained with the local expansion (25), is necessarily limited to short wavelengths (compared to the curvature radius), but it can be extended to the longest wavelengths [127][128][129]. This is necessarily more complicated because these modes suffer for the well-known gauge-dependence problems of cosmological perturbations. In order to solve it, a covariant and gauge-invariant formalism is needed. Let us assume that the background space is de Sitter and consider the general action which contains a mixture of f (R) and scalar-tensor gravity. On a Friedmann-Lemaître-Robertson-Walker background, the field equations become hold. An analysis [127][128][129], using the covariant and gauge-invariant Bardeen-Ellis-Bruni-Hwang formalism [130][131][132][133] in the version given by Hwang [134][135][136][137][138][139][140][141] for alternative gravitational theories, yields the stability condition of de Sitter space in metric f (R) gravity with respect to inhomogeneous perturbations This condition is obtained in the zero momentum limit. It coincides with the stability condition with respect to homogeneous perturbations [125].
At this point it is worth checking that the equivalence between metric f (R) gravity and Brans-Dicke theory with ω = 0 holds also at the level of perturbations. Previous issues [142,143] on this equivalence have been addressed.
For a Brans-Dicke model, with ω = 0, the stability condition of de Sitter solution with respect to inhomogeneous perturbations is given again by Eq. (34). For stability with respect to homogeneous perturbations, Eq. (34) holds provided that stability against local perturbations, expressed by f 0 > 0, is assumed. Therefore, there is a complete equivalence between metric f (R) gravity and Brans-Dicke theory with ω = 0 also at the level of cosmological perturbations.
From a more conceptual point of view, instabilities can be related to ghosts that are massive states of negative norm which cause a lack of unitarity and are common for generalizations of Einstein's gravity. The good news is that f (R) gravity is ghost-free. More general theories of the form f R, R αβ R αβ , R αβγδ R αβγδ , . . . , in general, contain ghost fields. A possible exception (under certain conditions [144,145]) is the case in which the extra terms appear in the Gauss-Bonnet combination G = R 2 − 4R αβ R αβ + R αβγδ R αβγδ , as in f = f (R, G). In this case, the field equations are of second order and there are no ghosts [146,147].
Another important issue is related with the weakfield limit. Early works on the weak-field limit of both metric and Palatini f (R) gravity was subject to errors and incompleteness (see [23] for details). A satisfactory treatment for the prototype model f (R) = R − μ 4 /R in the metric formalism was given in [148] and then generalized to arbitrary forms of the function f (R) in [10,149].

Cosmological perturbations dynamics
The expansion history of the universe alone is not sufficient to discriminate between various models, but the growth of structures depends on the theory of gravity and has the potential to achieve this goal. For example, in [155][156][157][158][159][160][161][162][163], it is assumed an expansion history a(t) typical of the ΛCDM model and found that vector and tensor modes are not affected by f (R) at lowest order, while perturbation scalar modes depend on the specific f (R) model. It is also found the condition f (R) > 0 for the stability of scalar perturbations, in agreement with the arguments discussed above. The most interesting results are that f (R) corrections lower the large angle anisotropies of the cosmic microwave background and produce correlations between cosmic microwave background and galaxy surveys different from those of dark energy models.
Overall, the study of structure formation in modified gravity is still work in progress, and often it is performed within the context of specific models, some of which are already in trouble because they do not pass the weak-field limit or the stability constraints. A similar situation holds for all Palatini f (R) models, and for this reason, their weak-field limit and cosmological perturbations are not discussed here.

The Cauchy problem
A physical theory must have predictive power and, therefore, a well-posed and well-formulated initial value problem. GR satisfies this requirement for "reasonable" forms of matter [6]. The well-posedness of the Cauchy problem for vacuum f (R) gravity was briefly discussed for special metric models in earlier papers [165]. Thanks to the equivalence between f (R) gravity and scalartensor theory when f (R) = 0, the Cauchy problem can be reduced to the analogous one for Brans-Dicke gravity with ω = 0, −3/2. Specifically, the initial value problem is well-posed and demonstrated for some scalartensor theories in [164,165]. A general analysis is performed in [166]. A general discussion of the Cauchy problem in the case of Palatini f (R) gravity is reported in [167][168][169]. Here, the well-posedness and the wellformulation are discussed according to the fluids sourcing the field equations.

Equivalence of f (R) gravity with scalar-tensor theories
An important issue is related with the fact that f (R) gravity can be a particular scalar-tensor theory. Specifically, if f (R) = 0, the metric f (R) gravity is equivalent to an ω = 0 Brans-Dicke theory 2 [170]. On the other hand, in Palatini formalism, f (R) gravity is equivalent to Brans-Dicke with ω = −3/2. This equivalence is discussed in Refs. [171][172][173][174][175][176]. Let us give details of both formulations in the following.

Metric formalism
Assuming that f (R) = 0 and starting with the Lagrangian (2), one can introduce the auxiliary scalar field φ and considers the action where The action (35) trivially reduces to (2) for metric f (R) gravity if φ = R. Vice-versa, the variation of (35) with respect to g αβ gives 2 The general form of the Brans-Dicke action is SBD = 1 2κ while varying with respect to φ yields and φ = R under the assumption f = 0. Hence, the scalar field φ = R is dynamical and satisfies the trace equation This scalar is massive and the analysis of small perturbations of de Sitter space allows one to compute explicitly its mass squared where a zero subscript denotes quantities evaluated at the constant curvature of the de Sitter background. It turns out to be more convenient considering the scalar field ψ ≡ f (φ), which satisfies ). It is clear, therefore, that the theory contains a scalar degree of freedom, and the action is recognized as an ω = 0 Brans-Dicke theory. This theory, called "massive dilaton gravity" was originally introduced in the 1970's in order to generate a Yukawa term in the Newtonian limit [177]. The assumption f = 0 can be seen as the requirement that the change of variable R → ψ(R) be invertible.

The Palatini formalism
In the Palatini case, the discussion of the equivalence with a Brans-Dicke theory proceeds in an analogous way. One begins with the action (8) and introduces φ = R and ψ ≡ f (φ). Then, apart from a boundary term that can be neglected, the action is rewritten, in terms of the metric g αβ and of its Ricci tensor R αβ , as where we used the fact that, sinceg αβ = ψ g αβ , the Ricci curvatures of g αβ andg αβ are related bỹ The action (43) is recognized as a Brans-Dicke theory with ω = −3/2. As a final consideration, we can state that f (R) gravity can be reduced in any case to the more general class of scalar-tensor theories and, in general, higher-order theories in curvature invariants can be reduced to multiscalar tensor theories [178].

Cosmology in extended gravity
One of the main purpose for introducing Extended Gravity is to achieve the cosmic accelerated expansion observed at late time [4]. However, besides the first models related to f (R) gravity, the approach has been expanded in view to consider also other geometric invariants like the Gauss-Bonnet one G, the scalar torsion T and so on.
Here we present and comment the cosmological equations for different Extended Gravity models. We focus in particular on f (R), Gauss-Bonnet f (R, G), and f (T ) cosmologies. The aim is to show that the equation of states for the models provide a sort of geometric adiabatic index which negative which can naturally explaining the accelerated phase of the today universe. In what follows, we shall refer to the spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) metric, which is usually adopted as the kinematic description of the universe, ds 2 = −dt 2 + a 2 (t) dx 2 + dy 2 + dz 2 . (45)

f (R) cosmology and the accelerated expansion
From the field equations (5), the cosmological equations can be derived. Here an overdot denotes differentiation with respect to the comoving time t. It is assumed that f > 0 to keep the effective gravitational coupling positive, and that f > 0 to avoid local instabilities [121,179]. The effective energy density and pressure of the f (R) fluid are . (49) The effective density ρ ef f is non-negative in a spatially flat FLRW universe in the limit ρ → 0. f (R) gravity can produce accelerated expansion without any material dark energy counterpart. In vacuo, the above equations assume the form [4,180] but, if a material cosmological fluid is present, it couples to gravity with the effective strength κ/f . One can define the effective EoS parameter as (52) In this approach, a metric f (R) model mimics the de Sitter equation of state w eff = −1 when By introducing explicitly the scalar degree of freedom of metric f (R) gravity, φ ≡ f (R), the effective EoS parameter becomes and a de Sitter solution corresponds toφ = f (R)Ṙ = 0. The above mechanism can, in principle, not only address the dark energy behavior but also the cosmological constant problem related to the fact that Λ does not evolve. Considering the above equations, the evolution of the gravitational vacuum state is allowed by the Extended Gravity. It is worth stressing that equations describing spatially homogeneous and isotropic cosmologies are of fourth order in the scale factor a(t). When matter is absent (a situation of interest in early time inflation or in a late universe completely dominated by f (R) corrections), a(t) appears only through the Hubble parameter H ≡ȧ/a. In this situation, it is convenient to adopt H, instead of a, as the dynamical variable. First, H is a cosmological observable; second, the above field equations are of third order in H. This elimination of a is not possible when the spatial 4s are not flat or when a fluid with density ρ(a) is present.
The dynamical fields of the theory are the metric g μν and the massive scalar degree of freedom φ ≡ f (R). As previously noted by several authors, quadratic corrections to the Hilbert-Einstein action introduce a massive scalar field [181][182][183][184][185][186]. This result is relevant for any metric f (R) theory [187][188][189].
In principle, the metric tensor contains various degrees of freedom: spin two modes, vector and scalar modes, and all of these can be massless or massive. GR contains only a massless graviton but when nonlinear corrections depending on R, R μν R μν , R μνρσ R μνρσ are included in the action, other modes have to be considered [8]. In the case of f (R) gravity, these further degrees of freedom include only a massive scalar mode, which is dynamical in the metric formalism but not in the Palatini formalism.
f (R) gravity can achieve cosmic acceleration through the effective equation of state parameter w eff −1, as it is well known from R 2 -inflation. This is possible also in the late universe, and it has even been attempted to unify early inflation and late time acceleration in the context of modified gravity [190][191][192][193]. However, modelling the late-time cosmic acceleration should not spoil the successes of the standard cosmological model which requires early inflation, a radiation era allowing Big Bang nucleosynthesis, a matter era during which matter overdensities can grow and form structures, and the present accelerated epoch leading to an uncertain future era the prediction of which is model-dependent (a de Sitter attractor solution or a Big Rip singularity are common predictions as discussed above). The transitions between consecutive eras must be smooth.
Smoothness may not be guaranteed in all f (R) models and the radiation-matter transition, in particular, can originate problems in specific f (R) models, including f = R − μ 2(n+1) /R n , n > 0 [194][195][196][197]. However, the prototypical toy model f (R) = R − μ 4 /R, which reportedly could not terminate the radiation era, was analyzed in detail with singular perturbation methods [198] and a sufficiently long matter era was found. In general, although a caveat exists about terminating the radiation era and allowing a sufficiently long matter era, one can always find choices of the function f (R) which achieve the correct cosmological dynamics (or any prescribed evolutionary history) by first assigning the desired form of the scale factor a(t) and then by integrating a differential equation for f (R) that produces the desired scale factor as discussed above. The solution f (R) solution is not unique [199][200][201][202], which shows that observational data providing information on (a segment of) the cosmic expansion history a(t) cannot suffice for a reconstruction of the function f (R) specifying the theory of gravity. Additional information is necessary, and it may come from the growth history of cosmological density perturbations, which depends on the theory of gravity.
Analytical solutions of FLRW cosmological equations are rare in literature [203] and phase space analyses are often necessary. They originated with pre-1998 studies of R 2 -inflation and are not limited to spatially flat FLRW spaces [18, [204][205][206].

The Gauss-Bonnet cosmology
A straightforward generalization of f (R) gravity involves the inclusion in the action of the Gauss-Bonnet topological invariant Considering the action F (R, G) means to take into account all the second order curvature invariants connected by the relation (56). Following the lines of Ref. [78] and adopting physical units such that c = k B = = 1, we can write the action It is worth saying that the action (57) contains all the possible curvature invariants that can be derived starting from the Riemann tensor. The variation of action (57) with respect to the metric provides the following field equations where κ = 8πG, T (matter) μν is the stress energy tensor describing the ordinary matter. The extra term T (GB) μν , containing all extra curvature terms, is defined as where ∇ is the covariant derivative and is the d'Alembert operator. As for f (R) case, Einstein's gravity is immediately recovered for F (R, G) = R.
In the following, for the sake of simplicity, we denote by the partial derivatives with respect to R and G. From the above definition (59), it is clear that the form of F (R, G) determines the dynamical behavior of the theory. In particular, the term (61) can be considered as an effective double-field potential for the theory where R and G act as two different scalar fields whose regimes can lead different phases of the cosmological evolution. This can be easily seen considering the trace of field equations Clearly we have the combination of two Klein-Gordon equations where the two fields can lead different phases of dynamics. Below, we are going to derive the cosmological equations in order to discuss the quintessential behavior. Assuming as above a spatially flat FLRW, we can write the Friedmann equations as The system of cosmological equations becomes selfconsistent considering the definition of the Ricci curvature scalar and the Gauss-Bonnet invariant in terms of the scale factor a and then the Hubble parameter H. As derived in [78], they are related to the Lagrange multipliers that have to be introduced in the action (57). We have We can also rewrite the total energy density and pressure, ρ (tot) and P (tot) , due to R and Gauss-Bonnet contributions as [83]: For cosmic acceleration, ρ (tot) +3P (tot) < 0, and assuming that all matter components have non-negative pressure, we can write the equation of state w (GB) = P (GB) /ρ (GB) from the geometry terms as where the quintessence behavior is obtained for −1 ≤ w (GB) < 0 while the phantom behaviour is achieved for w (GB) < −1.
Specifically, the Gauss-Bonnet term plays the role of a geometrical dark energy as in the case of f (R) gravity and then it contributes to the effective cosmological density according to the formula (69) This means that the effective value of the Gauss-Bonnet contribution can be "measured" by evaluating the standard matter and radiation contributions at the various epochs, that is However, the form of F (R, G) determines the evolution of w (GB) . A viable F (R, G) form can be achieved by considering the so called Noether Symmetry Approach [78]. It can be shown that symmetries select the form of the function to be where n is any real number and F 0 is a constant. For power law solutions of the form a(t) = a 0 t s [80], we obtain the relations n 1 = 1 + s 2 and n 2 = 1 1 + 2s(s − 1) − 2s , (72) such that Eq. (68) can be written in terms of n and s as: where n can assume the values n 1 or n 2 . Asymptotically w (GB) behaves like the ΛCDM model with w (GB) → −1.
By a cosmographic analysis reported in Ref. [87], it is possible to trace the whole cosmic history in terms of F (R, G) gravity. It emerges that Gauss-Bonnet field is dominant in early epochs while Ricci scalar field dominates in the late universe. Finally, it should be noted that F (R, G) = R + G, obtained simply adding the Gauss-Bonnet term to the Ricci scalar R, does not differ from GR in 4-dimensional spacetime as it can be rewritten as a boundary term. A detailed discussion on this point is reported in [88].

Recovering f (R) gravity from teleparallel gravity
In GR, gravitational interaction is described in terms of spacetime curvature. However, this is not the only possibility because equivalent formulations of GR can be achieved in flat spacetimes where gravity is described either by torsion or by non-metricity. This fact points out the existence of further representations of the same underlying theory [207]. Particularly relevant are the Teleparallel Equivalent of General Relativity (TEGR) and its generalization, f (T ) gravity, which do not require the Equivalence Principle as its foundation [25]. However, also in this framework an analogue description of accelerated cosmological dynamics can be achieved as in the case of curvature.
Let us sketch the main ingredients of TEGR approach and how f (R) gravity can be recovered also in this framework.
In this theory, the dynamical variables are the vierbeins or tetrad fields e a μ . They form an orthonormal basis for the tangent space at each point x μ of the spacetime manifold. The tetrads e n μ and their inverses E μ m fulfill the orthogonality conditions E μ m e n μ = δ n m , E ν m e m μ = δ ν μ . In terms of the tetrad fields, the metric tensor can be written as g μν = e a μ e b ν η ab , where η ab denotes the Minkowski metric. The idea underlying TEGR is to build up a theory with a geometry endorsed with torsion and with a globally flat curvature. This can be realized by defining the torsion tensor in terms of the curvatureless Weitzenböck connection W μ a ν = ∂ μ e a ν , that is In addition, it is convenient to define the contorsion tensor and the tensor 2S σ μν = K σ μν −δ μ σ T ν +δ ν σ T μ . The combination T = S σ μν T σ μν is the torsion scalar, 3 which is a topological object. The action corresponding to the TEGR is given by where e = det(e a μ ) = √ −g is the volume element of the metric. The Ricci scalar R and the torsion scalar T differs by a boundary term where B = (2/e)∂ μ (eT μ ) = ∇ μ T μ . It can be shown that the TEGR action reproduces the same field equations as GR being (76) equivalent to the Hilbert-Einstein action. The straightforwardly generalization of (76) is where f (T ) is a smooth function of the torsion scalar. Some comment are in order: (1) In this generalized theory is not possible to find the teleparallel equivalent of f (R) gravity since the boundary term in (77) contributes to the field equations. (2) Since the torsion field T is not invariant under local Lorentz transformations, this theory is also not invariant under Lorentz transformations. (3) The theory is a second order theory and hence, mathematically, it is easier than f (R) gravity. An important remark is in order at this point. As we said, the dynamical variables in TEGR are the tetrad fields e a μ , while, in the metric formulation, this role is fulfilled by the metric tensor g μν . This means the breaking of Lorentz invariance. However, it is worth noticing that TEGR and its generalizations can be reduced to equivalent metric theories, where the further degrees of freedom of the geometrical background are encoded in the metric potentials (see ref. [25], for more details). It means that passing from the anholonomic system of variables given by e a μ , to the holonomic system of spacetime coordinates, assigning g μν , the Lorentz invariance is restored.
The above action (78) can be immediately generalized by assuming that the function f (T ) depends also on the boundary term B. The action can be then written as [208] where f (T, B) is a smooth function of the scalar torsion T and of the (scalar) boundary term B. The motivation of this action is due to the fact that from the only f (T ) gravity, it is not possible to find an equivalent theory of its metric counterpart, the f (R) gravity. From the above action, instead, one can easily show that f (R) gravity can be recovered by assuming f (T, B) = f (−T + B) = f (R). The variation of the action (79) with respect to the tetrad field yields the field equations where f T = ∂f /∂T , f B = ∂f /∂B, ∇ σ is the covariant derivative with respect to the Levi-Civita connection and Θ λ a is the energy-momentum tensor. Let us consider f (T, B) cosmology in a flat FLRW Universe. In this perspective, one can construct the tetrad field e a μ = diag − 1, a(t), a(t), a(t) . Since f (T, B) is not invariant under Lorentz transformations, the choice of the tetrad deserves some caution (for instance, the unwanted condition f T T = 0 appears when one considers a flat diagonal FLRW tetrad in spherical coordinates). According to Ref. [209], the above choice of vierbein fields is a good tetrad. The Friedmann equations read As above, dots represent derivation with respect to the cosmic time and f T and f B are the derivatives of f (T, B) function with respect to T and B. Immediately, the cosmological equations can be written as 3H 2 = κ(ρ (matter) + ρ eff ), 3H 2 +Ḣ = −κ(P (matter) + P eff ), (84) and then the adiabatic factor is again w eff = P eff /ρ eff which means that the above considerations for f (R) gravity can be easily recovered.

Conclusions
In this short review paper we have discussed the possibility to address the cosmic history by Extended Gravity, i.e. extensions of GR aimed to address the gravitational interaction at ultraviolet and infrared scales. In particular, a straightforward extension of Einstein theory, f (R) gravity, can potentially explain the cosmic accelerated expansion by curvature both at early and late epochs. The approach is based on the fact that the further degrees of freedom, related to curvature, behave as a scalar field capable of leading the expansions in the various cosmic epochs without considering further material ingredients like dark energy and dark matter.
Taking into account some viability criteria in order to avoid inconsistencies and dangerous instabilities of the models, large classes of theories, involving also other geometric invariants, can be considered. The underlying philosophy is that geometry can give a comprehensive picture of the universe ranging from inflation up today. Clearly, the final cosmological model will be derived once a consistent amount of data, tracing back the evolution at any epoch, will be matched with methods of the precision cosmology. This inverse scattering approach is a sort of Designer Extended Gravity which will be more and more reliable if it will steam out from some fundamental theory.
Funding Open access funding provided by Università degli Studi di Salerno within the CRUI-CARE Agreement.

Data availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: In this paper we have not used data. The paper treat theoretical topics.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .