Higher-order gravity in higher dimensions: Geometrical origins of four-dimensional cosmology?

Determining cosmological field equations represents a still very debated matter and implies a wide discussion around different theoretical proposals. A suitable conceptual scheme could be represented by gravity models that naturally generalize Einstein Theory like higher order gravity theories and higher dimensional ones. Both of these two different approaches allow to define, at the effective level, Einstein field equations equipped with source-like energy momentum tensors of geometrical origin. In this paper, it is discussed the possibility to develop a five dimensional fourth order gravity model whose lower dimensional reduction could provide an interpretation of cosmological four dimensional matter-energy components. We describe the basic concepts of the model, the complete field equations formalism and the 5-Dim to 4-Dim reduction procedure. Five dimensional $f(R)$ field equations turn out to be equivalent, on the four dimensional hypersurfaces orthogonal to the extra-coordinate, to an Einstein like cosmological model with three matter-energy tensors related with higher derivative and higher dimensional counter-terms. By considering a gravity model $f(R)=f_0R^n$ it is investigated the possibility to obtain five dimensional power law solutions. The effective four dimensional picture and the behaviour of the geometrically induced sources are finally outlined in correspondence to simple cases of such higher dimensional solutions.


Introduction
Type Ia supernovae (SNeIa) observations depicted a latetime speeding up universe [1,2,3], driven by an unknown component, whose properties can be ascribed to some sort of exotic fluid. This elusive component has been addressed to a form of dark energy that exhibits an anti-gravitational negative equation of state (EoS) [4]. The combination of standard matter and dark energy provides an Einsteinlike Universe where gravitational attraction is counterbalanced by such an unusual gravitational source therefore reproducing cosmological observations. Soon after this discover a plethora of theoretical proposals has been suggested in order to explain dark energy origin. However, at today, a well endowed and self consistent physical interpretation is so far unknown [5,6,7]. The puzzling quest of a satisfactory explanation about these phenomenological results led cosmologists to explore several research lines. In particular, people followed two main directions. From one side, standard Einstein gravity has been reviewed by introducing a new cosmological component: the cosmological constant or a whatsoever well behaved fluid with a negative EoS. Differently, adopting an unconventional point of view, they have been considered modified gravity models that generalize Einstein theory: i.e. scalar-tensor gravity [8,9,10], f (R) theories [11,12,13], DGP gravity ⋆ e-mail:atroisi@unisa.it [14], braneworld scenarios [15,16,17], induced matter theory [18,19,20] and so on.
Among others f (R) theories of gravity received a considerable attention. Such theories, which represent a natural generalization of Einstein gravity, are obtained by relaxing the hypothesis of linearity contained in the Hilbert-Einstein Lagrangian. Several results, sometimes controversial, have been obtained in this framework [21]. These models have been satisfactory checked with cosmological observations [22,23,24] and suggest intriguing peculiarities in the low energy and small velocity limit [25] since the gravitational potential displays a Yukawa like correction [26]. In order to be viable f (R) gravity theories have to satisfy some minimal prescriptions. In particular, they have to match cosmological observations avoiding instabilities and ghost-like solutions and they have to evade solar system tests [27]. Constraints from energy conditions represent a further argument to settle a suitable gravity Lagrangian [28]. However, at this stage, a fully satisfactory fourth order gravity model is far to be achieved. As a matter of fact, such kind of models have been considered also from other points of view. For example, recently, models with a non-minimal coupling between gravity and the matter sector have been taken into account [29,30].
On the other side, historically, attempts for a unification of gravity with other interactions stimulated the search for theoretical schemes based on higher dimensions, i.e. beyond our conventional four dimensional (4-D) spacetime. Nordstrøm [31], who was the first to formulate a unified theory based on extra dimensions, Kaluza [32] and Klein [33] developed a five-dimensional (5-D) version of General Relativity (GR) in which electrodynamics derivates as a counter-effect of the extra-dimension. Successively, a lot of work has been dedicated to such theoretical proposal both considering fifth-coordinate compactification and large extra-dimensions, allowing for non compact mechanisms [19]. This last approach led to the so called Space-Time-Matter (STM) or Induced-Matter Theories (IMT) [18,20] and to braneworld theories [15]. The basic characteristic of IMT is that 5-D vacuum field equations can be recast after the reduction procedure as 4-D Einstenian field equations with a source of geometrical origin. In such an approach, the matter-energy source of 4-D spacetime represents a manifestation of extra dimensions. Generalizations of Kaluza-Klein theory represent themselves a suitable scheme to frame modern cosmological observations. In fact, after SNIe observations, several attempts have been pursued to match Induced Matter Theories with the dark energy puzzle [34,35,36,37].
In this paper we try to merge the two approaches. Since vacuum fourth order gravity theories, as well as IMT, can be cast as an Einstein model with a matter-energy source is of geometrical origin [38,39], it seems a significant proposal to confront the two theoretical schemes. The genuine idea underlying this work is to investigate 5-D f (R)gravity from the point of view of a completely geometric self-consistent approach. In this scheme, all matterenergy sources will represent a byproduct of the dimensional reduction related to higher order terms and higher dimensional quantities. Such an approach determines a completely different conceptual framework with respect to standard five dimensional fourth order gravity models. The main purpose is to recover along this scheme both dark matter and dark energy dynamical effects. In the standard realm, in fact, extra dimensions and curvature counter-terms will only play the role of mending dark energy phenomenology in presence of ordinary matter. Along this orthodox paradigm, 5-D models of fourth order gravity have been studied in presence of perfect fluid sources [42,43] under peculiar assumptions on the metric potentials. In a similar fashion, accelerating 4-D cosmologies, induced by generalized 5-D f (R) gravity, have been also studied considering a curvature-matter coupling [44]. Furthermore, a vacuum fourth order Kaluza-Klein theory, defined in term of the Gauss-Bonnet invariant, has been studied from the point of view field equations predictions. In particular, it has been investigated, in the cylinder approximation, the propagation of its electromagnetic degrees of freedom [40] and the particle spectrum in the linear regime [41].
Attempts to meet f (R)-gravity and IMT have been also developed in time. For example it has been investigated the possibility to draw information on the Space-Time-Matter tensor starting from the energy-momentum tensor induced by higher order curvature counter-terms [45]. On the other side a new effective coupled F ( (4) R, ϕ) gravity theory has been proposed as a consequence of a five di-mensional f (R) model [46]. Within such a paper it has been also demonstrated that the Dolgov-Kawasaki stability criterion of the 5-D f (R) theory remains the same as in usual f (R) theories: f ′′ (R) > 0, with prime that indicates the derivative with respect to the Ricci scalar. In our work we try to develop a more general framework, at first 5-D f (R)-gravity is discussed in a complete analytical form. We provide general field equations in the f (R)-IMT approach and, by exploiting the reduction procedure from 5-D to the 4-D ordinary spacetime (considering a suitable extension of Friedmann-Robertson-Walker metric), we obtain a complete set of new Einstein-like field equations with matter-energy sources induced by higher order derivative terms and higher dimensions (top-down reduction). In other words, we develop a fully geometrical cosmology where ordinary matter and dark components could be addressed, in principle, as the outgrowth of the GR formalism adopted within the top-down reduction mechanism. The paper is organized as follows. In section II we provide a brief review of 4-D f (R)-gravity. Section III is devoted to summarize Induced Matter Theory. The 5-D f (R)-gravity formalism is outlined in section IV, and, thereafter, we discuss the 5-D to 4-D reduction procedure of field equations. In section V it is described a routine to find 5-D power law solutions. Section VI is dedicated to the analysis of the results and, in particular, to their interpretation from the point of view of the 4-D induced cosmologies. Finally, Section VII is left to conclusions.

The Cosmological Equations for f (R) Gravity
Dark energy models are based on the underlying assumption that Einstein's General Relativity is indeed the correct theory of gravity. However, adopting a different point of view, both cosmic speed up and dark matter can be addressed to a breakdown of GR. As a matter of fact, one should consider the possibility to generalize the Hilbert-Einstein (H-E) Lagrangian. With these premises in mind, the choice of the gravity Lagrangian can be settled by means of data with the only prescription of adopting an "economic" strategy, i.e. only minimal generalization of H-E action are taken into account. The so called f (R)gravity [11,21,47], which consider an analytic function in term of Ricci scalar and provides fourth order field equations, is based on this conceptual scheme. It has to be reminded that higher order gravity theories represent the natural effective result of several theoretical schemes [48], i.e. quantum field theories on curved spacetimes and Mtheory. In addition, they have been widely studied as inflationary models in the early universe [49,50].
Fourth order gravity is favored by Ostrogradski theorem. It has been in fact demonstrated [51] that f (R)-Lagrangians are the only metric-based, local and potentially stable modifications of gravity among the several that can be constructed by means of the curvature tensor and possibly by means of its covariant derivatives.
Let us consider the generic f (R) action, we have: where we have used natural units 8πG = c = = 1, g is the determinant of the metric and R is the Ricci scalar and ψ characterizes the matter fields. Varying with respect to the metric we get field equations where, the prime, as said in the introduction, denotes derivative with respect to R and, as usual, Tracing (2) we obtain where T M = g µν T M µν , so that the relation between the Ricci scalar and T M is obtained by means a differential equation, differently than GR where R = − T M . This result suggest thatf (R)-gravity field equations admit a larger variety of solutions than Einstein's theory. In particular, T M = 0 solutions will no longer imply Ricci flat cosmologies and therefore that R = 0. Actually, a suitable property of this model is that field equations (2) can be recast in the Einstein form [11,38]: with that represents a curvature stress -energy tensor induced by higher order derivative terms (the terms f ′ (R) ;µν render the equations of fourth order). The limit f (R) → R reduces Eqs. (5) to the standard second-order Einstein field equations. The Einstein-like form of fourth order field equations (5) suggests that higher order counter-terms can play the role of a source-like component within gravity field equations. In practice, it is possible to postulate that geometry can play the role of a mass-energy component when higher than second order quantities are taken into account. Indeed, the trace equation (4) propagates a scalar-like degree of freedom. As a matter of fact, in principle, one can imagine to address Universe dark phenomenology to such a kind of effective fluid. From the cosmological point of view a relevant role is played by the barotropic factor of the curvature matter-energy fluid. In the vacuum case this quantity will discriminate accelerating solutions from standard matter ones. In the case of a Friedman spacetime, for example, one has [38,39]: where H =ȧ/a stands for the Hubble parameter and dot means a time derivative. In the following we will resort to a similar quantity in order to check the behaviour of the cosmological fluids induced by 5-D f (R)-gravity on the 4-D hypersurfaces that slice the higher dimensional spacetime along the extracoordinate. In our conventions we will consider latin indices like A, B, C, etc. running from 0 to 4, latin indices like i, j, etc., running from 1 to 3 and greek indices taking values from 0 to 3.

Induced Matter Theory
Modeling a unification theory that contains gravity and particle physics forces, typically, implies the resort to extradimensional models [19,52]. Among these 5-D Kaluza-Klein theory [32,33] and its modern revisitations Induced-Matter and membrane theory represent significant approaches; in addition these kind of models represent the low energy limit of more sophisticate theories (i.e. supergravity) [53,54]. The main difference between the seminal approach of Kaluza and Klein and Induced Matter Theory is related with the role of the extra-dimension. In its first conception the fifth dimension was "rolled up" to a very small size, answering the question of why we do not "see" the fifth dimension. Modern theories like IMT postulate that we are constrained to live in a smaller 4-D hypersurface embedded in a higher dimensional spacetime. The key-point for embedding the 4-D Einstein theory into the 5-D Kaluza-Klein Induced-Matter-Theory is represented by Campbell-Magaard (CM) theorem. The problem of embedding a Riemannian manifold in a Ricci-flat space was studied by Campbell soon after GR discovery [55], and finally demonstrated by Magaard in 1963 [56]. Further on, Tavakol and coworkers [57] used these studies to establish mathematical well endowed bases to the 4-D reinterpretation of 5-D Kaluza-Klein theory that is dubbed Induced-Matter Theory [58,60]. In such a framework, the 5-D to 4-D reduction procedure determines on the 4-D hypersurface an Einstein gravity theory plus induced matter components of geometrical origin. The CM theorem can be formalized as follows: Any analytic Riemannian space V n (x µ , t) can be locally embedded in a Ricci-flat Riemannian space V n+1 x A , t .
Here the "smaller" space has dimensionality n with µ = 0, ..., n − 1, while the "host" space has dimensionality N = n + 1 with A running from 0 to n, the extracoordinate can be both space-like and time-like. We are interested to the case N = 5. It is important to remark that CM theorem is a local embedding theorem. Therefore, more general issues related with global embedding, i.e. initial-value problems, stability or general induced solutions [59], cannot be resort to this achievement. However, for our purposes, the theorem guarantees the right analytic framework in order to frame 4-D matter phenomenology in relation to 5-D field equations [60].
Because of CM theorem it is possible to write down the 5-D metric as follows: where x µ = (x 0 , x 1 , x 2 , x 3 ) are 4-D coordinates, g µν turns out to be the spacetime metric and y is the fifth coordinate, ǫ = ±1 allows for space-like or time-like extra dimension. Along the paper we adopt a spacetime signature (+, −, −, −). With these premises in mind ordinary 4D spacetime result as a hypersurface Σ y : y = y 0 = constant, orthogonal to the 5D extra-coordinate basis vector If one considers vacuum 5-D field equations (Ricci flat): R AB = 0, the CM theorem, suggests a natural reduction process to a 4-D pseudo-Riemannian spacetime. In practice, one can obtain an Einstein-like 4-D model 1 G µν = T µν [58,60], once the right hand side of these equations fulfills the definition Here, a comma is the ordinary partial derivative, a semicolon denotes the ordinary 4D covariant and starred quantities describe terms derived with respect to the fifth coordinate. This quantity defines the so called Induced Matter Tensor; the only hypothesis underlying this achievement has been relaxing the cylinder condition within the Kaluza-Klein scheme (independence on the fifth coordinate).
As it is possible to observe from the previous result, fourth order gravity theories and Kaluza-Klein IMT models provide the same conceptual scheme. Both of the approaches allow to obtain an Einstein-like gravity model where matter-energy sources are of geometrical origin. In particular, these effective matter-energy tensors descend, respectively, from the higher order derivative contributions and from the higher dimensional counter-terms. Therefore, it seems that deviations from GR can be naturally recast as sources of standard Einstenian models. In that respect let us notice that IMT and 5-D f (R) gravity are built in a different philosophy with respect to conventional higher dimensional approaches like braneworld models. In fact, despite the same working scenario: bulk universe with non trivial dependence on the extra-coordinate, 4-D metric obtained evaluating the background metric at specific 4-D hypersurfaces, matter fields confined in the 4-D spacetime, there are intrinsic conceptual differences. 5-D f (R) gravity and Induced Matter Theory are based on the hypothesis that standard matter is nothing else than a 4-D manifestation of geometrical deviations with respect to GR. Braneworlds model our universe as a four dimensional singular hypersurface, the brane, embedded in a five dimensional anti De Sitter spacetime. In such a case the motivation for a non compactified extra-dimension is to solve the hierarchy problem. The differences in term of physical motivations for large extra-dimensions imply also different technical approaches. Within IMT and 5-D f (R) gravity one considers a Ricci flat vacuum bulk and develops the 4-D physics as a byproduct. In braneworlds it is assumed the opposite point of view. In this case, one deals with suitable 4-D solutions on the brane for some matter distributions and these solutions are matched with an appropriate 5-D bulk considering Israel junction conditions. Nevertheless, despite such conceptual differences it has been showed in time that brane theory, IMT and 5-D f (R) gravity can have a suitable matching. In particular, IMT and braneworlds have been put in strict analogy [61]. The key point is to consider the induced matter approach in term of the spacetime extrinsic curvature. In this case it is possible to write down braneworld-like field equations with a 4-D matter-energy source and a brane tension that are defined in term of the extrinsic curvature. On the same line very recently it has been shown that Z 2 braneworld-like solutions can be obtained as particular maximally symmetric solutions of 5-D f (R) gravity with matter [46]. In future works one can imagine to deepen the interconnections that seem to arise among these different higher dimensional approaches. Such considerations suggest, in principle, the intriguing possibility that GR experimental shortcomings could be actually addressed to the effective property of modified gravity models. In the following we will exploit such approach to study a general fourth order 5-D model, where higher order derivative gravity is "merged" with higher dimensions.

A 5-D f (R)-gravity model and its 4-D reduction
A f (R) theory of gravity in five dimensions can be described by the action where R (5) is the 5-D Ricci scalar, L m (g AB , ψ) is a Lagrangian density for matter fields denoted, as in the 4-D case, by ψ and g (5) represents the determinant of the 5-D metric tensor g AB . 5-D field equations can be obtained by action (12) varying with respect to the metric: is the energy-momentum tensor for matter sources, while ∇ A is the 5D covariant derivative, (5) = g AB ∇ A ∇ B is the 5D D'Alambertian operator. In order to simplify the notation we have defined f ,R (R (5) Together the field equations (13) one can obtain the trace where AB . Actually, 5-D field equations (13) can be naturally recast as a generalization of Einstein 5-D equations in the same manner of the ordinary 4-D formalism given in section 2. In fact, with some algebra and isolating 5-D Einstein tensor one gets : , by fact determining a 5-D new source of geometrical origin in the right member of field equations This fact, which can resemble only the byproduct of mathematical trick, determines significant physical consequences at 5-D and, as a consequence, also in the ordinary spacetime. In fact, as we will see in the following, the 4-D dimensional reduction of Eqs. (15) and (16) implies a "new" set of Einstein-like equations with three matter-energy components all of geometrical origin. One of these quantities descends from the 4-D reduction of the Einstein tensor and the two terms derivate from the relative 4-D reduction applied to T Since we want to explore how geometric counter-terms can effectively mimic cosmological sources we neglect ordinary matter, i.e. here on we will assume T Let us now develop the 5-D to 4-D reduction procedure of our higher dimensional framework. At first, by assuming the metric (31), it is possible to draw the reduction rules for the differential operators [36] : ,R , (18) ,R = f here, again, the asterisk denotes partial derivative with respect to the extra coordinate (i.e., ∂/∂y = * ); D α is the four-dimensional covariant derivative defined on the hypersurface Σ y=y0 , calculated with g µν , and the usual D'Alambertian is referred to 4-D quantities. In the same shape, all the quantities that are not labelled in term of the 5th coordinate will be intended to refer to ordinary spacetime. On this basis, it is possible to rewrite the 5-D field equations (15) by separating the spacetime part (µ = 0, .., 3; ν = 0, .., 3) from the extra-coordinate (y) one : Our purpose is now to disentangle the extra-coordinate dependence from the Einstein tensor. In this way it will be possible to isolate the 4-D part of Einstein tensor on the left member and move on the r.h.s all quantities that depend on higher derivative terms and on the extra-coordinate. All the higher order and higher dimensional counter-terms will play the role of effective source terms on the 4-D hypersurface. It is evident that since the fourth order Lagrangian f (R (5) ) depends on the 5-D Ricci scalar, this dependence cannot be completely untwined until the Lagrangian dependence is not specified. The reduction rules for the Einstein tensor [60] give . (23) The next step is to calculate the 5-D Ricci scalar 44 . After the substitution of the expression (23) within Eq.(21) and some tedious algebra one gets the relation ,R ) Φf that reminds an analogous expression given in [36] for the 4-D reduction of 5-D Brans-Dicke theory. Further manipulations of (24) allow to obtain a relatively simpler expression : that can be, finally, used in the reduction of Ricci scalar to collect some terms. Therefore, by using also the trace equation, one has : It is easy to check that this result reproduces the analogous relation of 5-D GR [19] when f (R (5) ) → R (5) . We remember that in such a case vacuum Einstein field equations reduce to the system R Here: is the usual IMT tensor, deriving from the 5-D to 4-D reduction of Einstein tensor, is the curvature tensor written preserving the 4-D form. In such a case extra-dimension contributions are still hidden in the scalar curvature nested within the definition of the gravity Lagrangian and of its derivatives. Finally, ,R ) Φf represents a mixed tensor containing terms that depend explicitly on the fifth coordinate and on the derivative of f ,R with respect to y. Therefore, we have obtained a framework where 4-D gravity is fully geometrized. In fact, within such a scheme matter-energy sources are related only with geometrical counter-effects deriving either from higher dimensional metric quantities or from higher order derivative ones.
In order to complete our discussion one should take into account also the off-diagonal equation. By considering that R [60], one obtains that coincides with similar expressions achieved elsewhere [36,46] and in the limit f (R (5) ) → R (5) gives back, in absence of ordinary matter, the conservation law P β α ;β = 0 [19]. Equation (30) resembles a more general conservation equation which relates the spacetime derivatives of P β α and * f (5) ,R . It has been conjectured [65] that the spacetime components of the field equations relate geometry with the macroscopic properties of matter, while the extracoordinate part ( α 4 ) and ( 4 4 ) might describe their microscopic ones. With these hypotheses, within fourth order gravity Kaluza Klein models, one would have that microscopic properties of matter are influenced by spacetime derivatives of the scalar degree of freedom intrinsic in f (R) gravity. Furthermore, in our case, P β α dynamics depends also on Φ, that, in turn (see Eq. (25)), shows an evolution driven by the higher order gravity counter-terms. It seems that there is a strict interconnection between the extra-dimension properties and the intrinsic f (R) scalar degree of freedom. For example, looking to equation (25), it seems that f ′ (R) guarantees a sort of Machian effect for this kind of models. However, all these considerations represent, at this stage, nothing more than speculations and, therefore, we do not discuss this issue further here.

5-D f(R) gravity cosmological solutions
Let us now verify what kind of solutions can be derived for our 5-D f (R)-gravity model. Once general cosmological solutions have been achieved, it is be possible to draw 5-D effects on Σ y=y0 hypersurfaces, deriving the effective 4-D picture. In particular, one can determine the spacetime effectual matter-energy behaviour of the different geometrical components described in the previous section. In order to search for field equations solutions we assume the 5-D metric: where k = 0, +1, −1 is referred to the 3-D spacetime curvature and (t, r, θ, φ) are the usual coordinates for spherically symmetric spatial sections. It is important to notice that, since the metric choice, our f (R ,R n 2 + 3ȧṅ an 3 +ḟ The last equation (35) recalls, again, the result obtained in [36] for 5-D Brans-Dicke vacuum solutions. Therefore it seems a suitable choice to pursue the same approach performed in this work to get some particular solutions of our equations system (32)- (35). If we consider the metric coefficients as separable functions of their arguments : n(t, y) = N (y), a(t, y) = P (y)Q(t), from Eq. 0 4 one gets : This equation suggests the possibility to look for power law cosmological solutions. In the following section we will discuss some examples in this sense.

Preliminary power law solutions
In order to look for 5-D f (R) gravity power law solutions we consider a set of metric potentials defined as follows [36]: with N 0 , A 0 , Φ 0 some constants with appropriate units and α, β, γ, δ, σ that represent the unknown parameters asked to satisfy 5-D field equations. Up to now we have developed a completely general scheme, without any assumption on the gravity Lagrangian. However, to completely define the solving algorithm one has to make a choice about the Ricci scalar function entering the gravity action. To remain conservative with the solution procedure, we consider a power law function of the Ricci scalar f (R) = f 0 R n . Such a model has been extensively studied in literature at four dimensions both in cosmology (curvature quintessence) [11,38,39] and in the low energy limit [25,64]. In addition, the 5-D phenomenology of a power law fourth order gravity has been investigated in the standard approach in presence of a perfect fluid matter source [43,44]. In particular, in such a case, it is assumed a compact fifth dimension and it is supposed that the homogeneously distributed fluid does not travel along the fifth dimension. As a matter of fact the matter energy density and the pressure experienced by a four-dimensional observer have to be an integrated throughout the compact extra-dimensional ring. In our approach we have overcome this point of view, we relax the hypotheses on the extracoordinate and, above all, we discard ordinary matter in favor of a completely geometric self-consistent approach. According with cosmological observations [62,63], we will study field equations considering a flat spacetime geometry (k = 0).
A suitable recipe to study Eqs. (32)- (35), is to plug the metric functions (38) within the last of these relations and to search for its solutions. This approach allows to obtain, with some efforts, a number of constraints on the model parameters. Then, one can try to satisfy also the other more complicated field equations (32)- (34). By inserting metric potentials within Eq. 0 4 , one obtains a rather cumbersome expression : Ξ(α, β, γ, δ, σ, n, t, y) Υ (α, β, γ, δ, σ, n, t, y) = 0 where Ξ(α, β, γ, δ, σ, n, t, y) = y −1−2δ 3(βγ + α(−β + δ)) N 2 This result, evidently, suggests the possibility to obtain, in principle, much more solutions than the Brans-Dicke case studied in [36]. To simplify our search, as a preliminary step, we make some trivial hypotheses about model parameters, leaving a complete analysis of field equations solutions to a future dedicated work.

Kaluza-Klein GR limit
As a first step we want to verify that the Kaluza-Klein GR limit is recovered. To perform standard KK limit we have to settle n = 1 and, in addition, following customary approaches to the model, we assume σ = 0.

Standard 4-D f (R) gravity limit
What about standard 4-D fourth order gravity? One should expect that this framework has to represent a natural subcase of 5-D f (R) gravity. In order to get this limit one has to assign a solutions set with δ = 0, β = 0, σ = 0 and γ = 0, that means, no dependence on y and no dynamics on the fifth coordinate. Starting with this assumptions and by considering equations (32)- (33) and (35), one obtains that exactly matches the well known power law f (R) gravity solution in the case of vacuum field equations [11]. However, in the 5-D case, to fully satisfy the theory we have to fulfill an equation more, Eq. (34). This request, by fact, settles the power law index n. Therefore the only allowed combination is δ = 0, β = 0, σ = 0, γ = 0, n = 5/4 and α = 1/2.

Cylinder-type solutions
We can now relax some hypotheses in order to look for generalizations of 5-D GR. To search for such a kind of solutions we start, as a first simple case, from the assumption of no dependence on the extra-coordinate. This means to study our theory, looking for solutions of Eqs.(32)-(35), in the limit of cylinder condition. Differently than the standard case, in the f (R) generalization of Induced Matter Theory this condition, as we will see, does not imply radiation as the only possible kind of induced matter. At first, let us gibe a look to Eq. (25). In the case of 5-D f (R) gravity, independence on the fifth coordinate does not mean to have a massless Klein-Gordon equation for the extra-coordinate potential Φ. Higher order gravity counter-terms will play the role of a mass term and Φ will play a completely new role into the dynamics : In addition, as a consequence, T IMT µν will be no more traced to zero [19], therefore the Induced Matter Tensor can span more general kinds of matter in relation to the underlying cosmological solution.

4-D induced f (R) gravity: a top-down geometrization of matter
The effective 4-D picture induced by higher dimensions can be achieved once cosmological solutions obtained in Sec.5 are plugged within the 5-D to 4-D reduction framework previously outlined. In particular the 4-D setting is achieved, as already said, by considering hypersurfaces Σ y=y0 that slice the 5-D universe along the fifth coordinate. In this scheme, the extra-coordinate effects on 4-D physical quantities are evaluated taking y = y 0 = const. The resulting cosmological model is a 4-D induced matter f (R) theory, where field equations are given in the form of Einstein equations equipped with three different matterenergy sources of geometrical origin. Deviations from standard GR provide, in the effective 4-D spacetime description, matter-energy sources and one could wonder if these contributions are related with the elusive nature of Dark Energy and Dark Matter. Actually, we study the 4-D induced cosmologies in correspondence of the solutions given in Sec.5. In particular, we are interested to investigate the behaviour of the three different induced cosmological fluids (11), (28), (29) on the y = y 0 = const hypersurfaces. One can notice that for all of these quantities it holds

4-D GR limit
If one considers solutions obtained along the standard KK GR limit, the related 4-D effective dynamics are described by well known cosmological models. In particular, the solution (41) represents the standard radiation universe with p = 1/3ρ. In such a case the three matter-energy tensors collapse into one matter-energy source, in particular the curvature quintessence tensor and the mixed tensor respectively vanish. On the other side, the solution (42) represents a spatially flat FRW metric once the extracoordinate is fixed : The induced matter tensor plays, again, the role of the only one effective cosmological fluid. In fact, the other two quantities combine to give a cancelling result. The cosmological energy density scales as ρ = ρ 0 t −2 , with ρ 0 = 3α 2 N 2 0 y 2 0 , while the barotropic factor is ω = 2 − 3α 3α . This result is in complete agreement with previous achievements on non compactified KK gravity [36].

4-D cosmology from standard f (R) gravity limit
We have seen that the attempt to get the standard 4-D f (R) gravity limit is frustrated by the higher equations number of our model. In fact, it is possible to get the right 4-D power law solution only if the extra-coordinate equation is neglected. When this equations enters the game the power law index n gets fixed. The behaviour of the three matter-energy tensors describing the top-down effect of higher dimensions and higher derivative terms on the 4-D gravity hypersurface confirmate such a result. In fact, if we evaluate ω curv on the 4-D f (R) solution (43) we obtain that, again, matches the already known result for this kind of model [11]. However, since the fifth equation select the allowed values for n, we know that the admissible solution is indeed δ = 0, β = 0, σ = 0, γ = 0, n = 5/4 and α = 1/2. The effective, 4-D, consequent cosmology shows ρ IMT = ρ Mix = 0, ρ curv = 3 4N 2 0 t 2 and ω curv = 1/3. As a matter of fact, we find a new radiation like universe with n = 5/4. The radiation component derivates from the higher derivatives counter-terms while the Induced Matter Tensor vanishes.

4-D Cylinder models
Cosmologies induced by cylinder solutions are quite interesting and allow different dynamics. In particular, the matter-energy sector, characterized by means of the three matter-energy sources of geometrical origin, assumes a variety of significant behaviours. A relevant aspect of these solutions is that the mixed tensor T Mix µν provides a cosmological constant-like source. In fact, as it is possible to observe from the definition (29), in such a case (no extracoordinate dependence) it is T 0 Mix depending on the solution we consider, we obtain an inverse square law cosmological constant. This behaviour is favored in string cosmologies [67] and in time-varying Λ theories that accommodate a large Λ for early times and a negligible Λ for late times [68]. On the other side, the effective barotropic factors of the other two matter-energy sources depend on the parameter n in relation to the solutions (45a)-(45c). In order to have a schematic portrait of these result we have plotted the behaviour of each ω K with respect to n in Fig.1. It is possible to observe that the value of the gravity Lagrangian power index n determines the possibility to have different cosmological components. Both the standard matter sector and the dark energy one can be mimicked by the geometrically induced fluids according with previous results in this sense obtained in the framework of IMT [58,69,70] and fourth order gravity [22]. Very interesting is the coexistence of different regimes with the same value of n.
Together the behaviour of ω K we display, for all the components (see Fig.2), the energy density value ρ K and the overall sum of these quantities 5 . It is evident that energy conditions are violated for certain values of n. For example it is easy to observe that the Weak Energy Condition, ρ i ≥ 0, for the first solution, is not satisfied by ρ IMT and alternatively by the other two energy densities ρ Curv and ρ Mix . For the second solution we have again no n intervals where all components satisfy this energy prescription, since ρ IMT is quite always well behaved but ρ Curv and ρ Mix are alternatively negative. Finally, for the third solution we have a small region around n ≈ 1.5 where all energy densities are positive, while the WEC is violated elsewhere. Although single matter-energy sources can violate the energy conditions, the total geometrized matter satisfies the WEC (dotted lines in each panel of Fig.2) for all the ranges of the parameter space in which cylinder solutions are defined. In addition, it has been recently showed that cosmological models with negative energy components can be suitably matched with observations if these "disturbing" quantities are coupled with greater magnitude of positive energy forms [71]. Therefore the total top down geometrized matter is physical and preserves conceptions about energy conditions in 4-D. Such a result is quite customary when dealing with Induced Matter Theory models [19]. 5 In order to achieve these plots we settle the quantity N 2 0 t 2 = 1. Finally, let us discuss the cosmological evolution of the four-dimensional sections determined by each of solution (45a)-(45c) and the behaviour of the extra dimension. Figure 3 allow to draw significant informations about this behaviour describing the dependence of a(t, y) and Φ(t, y) with respect to n in the different cases. In order to display an accelerating dynamics, a power law scale factor requires that the power index has to be negative or greater than one. In particular, universes expanding at an accelerating rate need α > 1 if we consider the definition given in Eqs. (38).
Actually, since our approach in the 4-D limit determines a multi-fluid model, the acceleration condition for the cosmological dynamics turns out to be characterized by the general relation for the scale factor Therefore, barotropic fluids like ones outlined by our model, with ρ tot = K ρ K and p tot = K p K , imply that such a relation can be written as and the condition for accelerating cosmologies becomes Υ = K (1 + 3ω K ) ρ K < 0. As one can see in Fig. 4 for our model the regions of the parameters space that fulfill this requirement exactly determine accelerating expansions (the deceleration parameter q = −ä ȧ a 2 is negative). Let us remember that the observations suggest q ≈ −0.5 [72,73], albeit models considering power law solutions can also allow for higher limits i.e. q ≈ −0.3 [74]. A peculiar aspect of our model is represented by the possibility that singles energy densities can achieve negative defined values. Of course, as said above, the overall evolution is still physical since the effective energy density budget is always positive definite. This behaviour, however, can affect the effective cosmological evolution. For example, it can prevent the possibility to obtain models with speeding up evolution that accord with cosmological observations. In fact, it is possible to observe (cfr. Fig.1 and Fig.4) that there are regions of the parameter space where the induced cosmological fluids all assume negative barotropic factors while cosmological dynamics is still decelerating (q ≥ 0). On the other side, this peculiarity can also determine suitable conditions to frame models that experience a cosmic speeding up.
It is evident that, at this stage, our results can be only considered as a toy model achievement. Power law solutions can only describe limited stages of universe evolution. More general analyses are required in order to understand the real predictions of the model. In particular, the possibility to have a coexistence of negative defined equations of state with both decelerating dynamics and accelerating ones can represent an intriguing perspective to deepen in the framework of more general cosmological solutions. Let us discuss in some detail the results contained in Fig.1 This solution is allowed in the interval 1 2 < n < 13 10 . Values of n < 1 give three negative EoS fluids within the cosmic pie. The GR limit (n = 1) implies a trivial behaviour with vanishing ρ K . For n ∼ 1.3 we have two standard matter components and the cosmological constant like source driven by T Mix µν . In such a case ρ IMT , ρ Curv > 0 while ρ Mix < 0, the whole cosmic mass-energy budget is positively defined as already said above. Looking to the cosmological dynamics one can notice that this solution implies a negative power law index for the scale factor within the interval 1 2 < n < 1 while 0 ≤ α < 1 for 1 < n < 13 10 . According with Fig.3-4 this behaviour indicates accelerated contracting universes in the first interval and standard matter-like dominated cosmologies in the second one. The extradimension is always decreasing since is 0 < γ < 1. It is interesting to notice that the ordinary matter-like regime is achieved: in the interval 1 < n < 1.25 with two of the three fluids that behave as sources with a negative defined EoS (the other one is a stiffed fluid with ω IMT > 1) and in the range 1.25 < n < 13 10 with two standard matter-like sources (0 < ω IMT < 1 and ω Curv > 1) together the cosmological constant contribute. The divergence of the barotropic factor ω Curv for n = 1.25 is regulated by the correspondent vanishing of ρ Curv .
-Solution (α 2 , γ 2 )-(45b); We have again 1 2 < n < 13 10 . For this solution, the scale factor shows a standard matter rate since is always 0 < α < 1 2 , there is no acceleration (q(α 2 , γ 2 ) > 0). The induced matter component ρ IMT is always a well behaved thermodynamical fluid with 0 ≤ ω IMT ≤ 2/3. The curvature induced sources can have a negative EoS or a positive one with −5/6 ≤ ω Curv ≤ 3/2. For values n 1.5 we have two standard matter components and a cosmological constant like source. Energy densities can be positive either negative defined. Around n ∼ 1.2 we have two standard matter-like sources with positive ρ i , while ρ Mix is negative, ρ tot is always ≥ 0. For n = 1, the two curvature induced sources collapse to a zerovalued energy density, while the induced matter tensor plays the role of a radiation fluid with ω IMT = 1/3 as in usual KK-GR cylinder models. The fifth dimension can experience a reduction when 1 2 < n < 5 4 and an increasing evolution for 5 4 < n < 13 10 , is static when Admissible values for n ∈ R − {5/2}. This solution explicitly resembles the standard 4-D f (R)-gravity solution given in [38]. The scale factor and the extradimension show the same dynamics since α 3 = γ 3 . There are two interesting regions around n ∼ 0.2 and n ∼ 1.6. In the first case we have two standard matter components, dust-like and radiation-like, with ω Curv ≈ 0, ω IMT ≈ 0.26 together a cosmological constant-like one with ω Mix = −1. The correspondent scale factor is however non accelerating since q is positive. The sec-  ond region shows a dust-like source with ω Curv ≈ 0 and two dark energy-like components, one of them phantom-like with ω IMT ≈ −2.6. In such a case the 4-D section shows an accelerating expansion since q < 0 (q 0 = −0.5 for n ≈ 1.69) as well as an increasing extradimension. The zeros of ω IMT (dust-like behaviour of this source) are trivial since the correspondent value of the energy density is vanishing. Actually, this solution allows a plethora of different behaviours. In particular, expanding solutions with no acceleration are obtained for − 1 2 ≤ n < 1 2 and 1 < n ≤ 3 2 with two standard matter-like sources plus the effective cosmological constant contribute given by T Mix µν ; on the other side accelerated expansions are given by models with n < − 1 2 and 3 2 < n < 5 2 . In these in-tervals T Curv µν plays the role of a baryonic component and the other two matter-energy sources act as dark energy-like sources. The energy densities have, again, both negative and positive values with ρ tot ≥ 0. For n > 1.5 we have all positive ρ K with ω K ∼ −1, in such a case the three induced components behave all as cosmological constant-like sources.
As a matter of fact, fourth order gravity provide, in the cylindricity case, a completely different framework with respect to the standard Induced-Matter Theory. The standard scheme implies, in the 4-D limit, that induced matter behaves as a radiation-like fluid (p = ρ/3) [19] according with the result obtained with an initial condition of cylindricity by Kaluza [32]. In the case of higher order Kaluza-Klein models the energy-momentum tensors en- tering the 4-D induced field equations can describe more general matter sources. Therefore one can obtain dynamically significant cosmologies also in presence of a Killing symmetry for the extra-coordinate. Of course, the significance of such solutions have to be corroborated with cosmological and astrophysical data.

Conclusions
We discussed the possibility to develop a higher order and higher dimensional gravity model. Within such an approach, after the 5-D to 4-D reduction procedure, it has been possible to suggest a new interpretive scheme for the 4-D cosmological phenomenology. Specifically, we have considered a vacuum 5-D f (R)-gravity model. The 4-D effective reduction provides, in this case, a GR-like cos-mology characterized by means of three induced matterenergy tensors of geometrical origin. Within this matter geometrization paradigm the cosmological fluid sources can be achieved considering: the standard Induced Matter Tensor of non-compactified KK theory T IMT µν , the Curvature Quintessence tensor T Curv µν coming from higher order derivative terms of f (R)-gravity and T Mix µν a mixed quantity that arises from both higher dimensional and higher order curvature counter-terms. This last quantity embodies the 5-D to 4-D reduction of the two combined approaches. We have obtained the complete field equations formalism and in order to check model prediction we have worked out a routine to search for power law solutions in the case f (R) = f 0 R n . The whole scheme provides the right KK-GR limit once f (R) → R. In such a case, on shell, radiation dominated cosmologies (p = 1/3ρ) are obtained if no dependence on the extra-coordinate is taken into account. More generally, the KK-GR limit furnishes solutions that agree with similar results already present in literature. In particular, the curvature induced matter-energy tensors vanish while the Induced Matter Tensor, that collects all metric terms depending on extracoordinate, turns out to be the only one source of 4-D field equation. Relaxing the hypotheses on n, i.e. by assuming a generic power law fourth order gravity model, it is possible to search for more general solutions. As a first preliminary step, cylindricity has been taken into account. A more general study has been left to a dedicated forthcoming publication. Differently from the standard KK-GR theory, the combination of fourth order gravity with higher dimensions implies, in such a case, non trivial solutions that can give rise to interesting cosmologies. We have obtained a set of solutions, parameterized by the f (R) power index n, that determine a variety of possible barotropic fluids in the 4-D mass-energy sector. In relation to the value of n one can have both standard matter components and dark energy-like ones. The cylinder condition implies that T Mix µν , according to its definition, behaves as an effective cosmological constant-like term with ω Mix = −1. Actually, in such a case, the energy density ρ Mix depends on time with an inverse square law. This fact suggests the intriguing possibility to have a time varying cosmological constant. The cosmological behaviour of the 4-D sections shows that it is possible to have both accelerating and non accelerating expansions. In particular, non accelerated expansion can be obtained also in presence of subdominant dark-energy like components. On the other side, it can be easily framed a dark energy dominated speeding up universe generalizing results obtained in the standard 4-D f (R)-gravity. The interesting aspect of this last solutions is that geometry can provide both standard matter-like sources and dark-energy-like ones, one of this can mimic a phantom fluid. The extra-dimension follows an independent evolution, that, in the third solution, is equivalent to one of the scale factor. Therefore, in such a case, expanding accelerated evolutions are characterized also by a growing fifth coordinate. The price to pay for geometrizing matter-energy sources is that cosmological energy densities can assume negative values violating the WEC. We have shown that although this is true for some values of n, also in interesting regions of the parameters space, the whole mass-energy budget ρ tot is always positive definite and, therefore, total topdown geometrized matter is physical.
The fact that geometrically induced matter-energy tensors could, alternatively, play the role of standard matter components and of dark energy-like sources can represent an intriguing perspective in order to interpret cosmological observations. In particular, this achievement seems to deserve some more studies in the direction of dark energy and dark matter interpretation. Of course, more insights are necessary in order to evaluate such a theoretical scheme. In particular, the matter sector (i.e. matter-like components), which has been satisfactory checked within standard IMT [69,75], require a careful investigation since matter properties should be completely induced by gravity. In addition, model predictions have also to be checked also against cosmological and astrophysical data.
Furthermore, if f (R)-Kaluza-Klein gravity determine, in the 5-D to 4-D reduction of cylinder solutions, cosmological models with a suitable top-down matter geometrization, in principle, more general solutions could allow to develop cosmological models with a wider phenomenology. On the other side, the fifth dimension role requires a careful analysis when the extra-coordinate Killing symmetry is discarded. In fact, if the cylinder condition is relaxed one naturally gets into the game of a length scale related with the extra-coordinate. In this sense, it has been shown that 5-D f (R) gravity suggests that, eventually, there is a strict interconnection between the fifth dimension dynamics and the higher derivative counter-terms.