Finsler-Randers-Sasaki gravity and cosmology

We present for the first time a Friedmann-like construction in the framework of an osculating Finsler-Randers-Sasaki geometry. In particular, we consider a vector field in the metric on a Lorentz tangent bundle, and thus the curvatures of horizontal and vertical spaces, as well as the extra contributions of torsion and non-linear connection, provide an intrinsic richer geometrical structure, with additional degrees of freedom, that lead to extra terms in the field equations. Applying these modified field equations at a cosmological setup we extract the generalized Friedmann equations for the horizontal and vertical space, showing that we obtain an effective dark energy sector arising from the richer underlying structure of the tangent bundle. Additionally, as it is common in Finsler-like constructions, we obtain an effective interaction between matter and geometry. Finally, we consider a specific model and we show that it can describe the sequence of matter and dark-energy epochs, and that the dark-energy equation of state can lie in the quintessence or phantom regimes, or cross the phantom divide.


I. INTRODUCTION
General Relativity (GR) has been proved a successful theory of gravity, tested with high precision at Earth-based and Solar System experiments (perihelion preccession of Mercury, gravitational redshift, Shapiro time-delay effect, etc) [1].Nevertheless, at the theoretical level one faces the problem of non-renormalizability [2], since GR cannot be incorporated in a quantum description [3].Additionally, at the cosmological level we still have the open issues of dark matter and dark energy [4,5], as well as possible tensions between predictions and observations, such as the H 0 [6] and σ 8 tensions [7] (for a review see [8]).Hence, a large amount of research was devoted in the construction of various gravitational modifications, namely theories that possess general relativity at as a particular limit, but which in general exhibit richer behavior, theoretically and cosmologically improved [9][10][11][12].
However, one can proceed to more radical geometrical modifications, namely use Finsler and Finsler-like geometries, which have richer structure than Riemannian framework, and use them in order to construct gravitational theories .These modified theories of gravity use generalized metric structures, where a vector field is incorporated in the geometrical construction, and have contributed with different directions in the development of locally anisotropic models for the gravitational field theory and cosmology.
In Finsler and Finsler-like geometries more than one connection and curvature appear, which depend on the position and velocity, in contrast to GR in which there is only the Levi-Civita connection and the curvature of the Riemannian space.Therefore, gravity can be studied in a different way in the framework of an 8-dimensional Lorentz tangent bundle or a vector bundle which includes the observer (velocity/tangent vector) with extra internal/dynamical degrees of freedom [25][26][27][28][29][30]54], as well as in an oscullating Riemannian and Barthel framework [55][56][57].
Concerning this approach, all kinds of generalized metric theories belong to the larger class of the so called "anisotropic field theories", since Lorentz violations, velocity fields and torsions produce anisotropies in the space and the matter sector [40,42,58,59].Hence, these internal geometrical anisotropies, which should not be confused with spacetime anisotropies that may appear in Riemannian geometry too (e.g. in Bianchi cases) are induced by internal direction/velocity y-variables in addition to the position x-variables in the structure of the base manifold.In these lines, geometrical anisotropies can be considered as a "potential" or a tidal field in the matter sector [60].In cases where an anisotropy is included in the metric structure of spacetime, as it appears in Finsler and Finsler-like cosmologies, it is incorporated in the effective energy-momentum tensor of the anisotropic structure, which could potentially lead to energy exchange between geometry and matter [61].Finally, similarly to general relativity, geometrical effects are produced not only by the distribution of mass-energy but also by its motion [62].
In this work we propose a novel geometrical structure, namely that of Finsler-Randers-Sasaki (F-R-S) type, in order to extract generalized gravitational field equations.Then, applying them in a cosmological framework we construct F-R-S cosmology, which is characterized by modified Friedmann equations with new terms that depend on the underlying geometry and the tangent bundle features.As we will show, these terms can lead to interesting cosmological implications, and describe the thermal history of the Universe, as well as the effective dark energy sector.The plan of the work is the following: In Section II we provide the basic concepts of F-R-S geometry, and we extract the general gravitational field equations.Then, in Section III we proceed to the application at a cosmological framework, and we extract the modified Friedmann equations for the horizontal and vertical subspaces, investigating also specific examples.Finally, in Section IV we discuss our main results.

II. FINSLER-RANDERS-SASAKI GEOMETRY AND GRAVITY
In this section we present the basics of Finsler-Randers-Sasaki geometry and gravity.We will start by introducing some geometrical aspects from Finsler geometry and the oscullating Riemannian metric, and then we will use it to construct a gravitational theory.

A. Finsler-Randers-Sasaki geometry with oscullating Riemannian metric
We consider an n-dimensional bundle M , as well as its tangent bundle T M , with a fibered and differentiable (smooth) metric function F (x, y) with the following properties: 1. F is continuous on T M and smooth on T M ≡ T M \ {0}, namely the tangent bundle minus the null set 2. F is positively homogeneous of first degree on its second argument: 3. For each x ∈ M the fundamental metric tensor : is non-singular, with i, j = 0, 1, ..., n − 1.
A Lorentz tangent bundle T M over a spacetime 4-dimensional manifold M is a fibered 8-dimensional manifold with local coordinates {x µ , y a }, where the Greek indices of the spacetime variables x are κ, λ, µ, ν, . . .= 0, . . ., 3 and the Latin indices of the fiber variables y are a, b, . . ., f = 0, . . ., 3.An extended Lorentzian structure on T M can be provided if the background manifold is equipped with a Lorentz metric tensor of signature (−1, . . ., 1).As it is known, a metric following the above three properties is called a Finsler metric [22,23].Additionally, one can introduce the oscullating Riemannian metric on a differentiable manifold [63].In particular, this can be defined by a tangent vector field Y : U → T U , where U ⊂ M is an open neighborhood on M with the property Y (x) = 0 ∀x ∈ U .In such a case the metric can be defined by the relation: In the following we will consider that all non-vanishing global vector fields, defined on the spacetime manifold, satisfy M y(x) = Y (x) [64].The pair (U, g ij (x, y(x)) is called Y -oscullating Riemannian metric associated to (M, F ) manifold.
As it is known, the length of a curve c in a Finsler space is given by the integral with y = dx dτ and τ an affine parameter along the curve.In a Finsler-Randers (FR) spacetime [65] the metric is given by the relation [66,67]: where a µν (x) is a Riemannian pseudo-metric and f α a covector with ||f α || ≪ 1.Note that the 1-form f α dx α can be interpreted as the "energy" produced by the anisotropic force field f α , and hence due to (4) and ( 5) the integral b a F (x, dx) represents the "total work" that a particle requires in order to move along a path with proper time τ .We proceed by writing the corresponding Lagrangian function of an FR space with an oscullating Riemannian metric, which is with ||f α || ≪ 1.We mention here that in the above expression the second term f α (x)y α (x) can be interpreted as the "power" that is produced due to propagation of particles through the force-field f α (x).Now, from (2),( 3) and ( 6) we can extract the metric tensor v αβ (x, y(x)) as where with L = −g αβ y α (x)y β (x) and ||A γ || << 1. Due to relation ( 8), the metric ( 7) can be called "weak Finslerian metric" [67].Hence, as we can see from (7), the term h αβ (x, y) can be considered as a perturbation.For convenience, and in order to make notation lighter, in the following we will write y instead of y(x).
Let us now introduce the Sasaki-type metric on T M [68].Such a metric has the form: In our approach we consider that the Finslerian metric v αβ (x, y) is given by ( 7), and the unified adapted frame is defined in the form and and where E A is the adapted basis of the tangent space T x M .Furthermore, we define the dual basis E A = (dx µ , δy α ) with where E A is the adapted basis of the cotangent bundle T * M and N α λ are the components of the nonlinear connection with α, λ = (0, 1, 2, 3).The nonzero coefficients of a canonical and distinguished d−connection D on T M read as [22]: Finally, concerning the non-linear connection, we choose a Cartan-type of the form: which is known to have interesting applications [69].
We now have all the geometrical quantities in order to calculate the curvature tensors.In particular, in such a framework the Riemann and Ricci curvature tensors of the horizontal space are defined as [22,23]: where Ω α νκ represents the curvature of the nonlinear connection and is defined as Moreover, the curvature tensors of the vertical space are given by: Thus, the generalized Ricci scalar curvature in the adapted basis is: where In the same lines, one can define the torsion tensor as where L α βν is given in (14).

B. Finsler-Randers-Sasaki gravity
Having presented the Finsler-Randers-Sasaki geometrical framework, we can use it in order to construct a gravitational theory.An Einstein-Hilbert-like action on T M can be defined as for some closed subspace N ⊂ T M , where L M is the standard matter Lagrangian and κ is the gravitational constant.Note that while the absolute value of the metric determinant |G| is |G| = √ −g √ −v, with g, v the determinants of the metrics g µν , v αβ respectively, as it follows from the form of (9).
Performing variation of the above action in terms of g µν , v αβ and N α κ we extract the field equations as (the details are presented in Appendix A): where we have defined the "energy-momentum tensors" where δ µ ν and δ α β are the Kronecker symbols.Note that the second field equation ( 29) is simplified to: where in our case S αβ and S are zero and the mixed Cartan symbols C µ νγ are also zero as we can see from (15).Additionally, we mention that the third field equation ( 30) is identically zero since the coefficients L ρ µν are identical to the classical to Christoffel coefficients of the FRW model, thus they do not depend on the vertical coordinates y α .The mixed Cartan coefficients C µ νγ are zero as we mentioned and the mixed energy-momentum tensor Z κ α is also zero from relation (33) since we do not consider the matter field to depend on the non-linear connection N α κ .Let us make some remarks in order to give a physical interpretation of equations ( 28), ( 29) and (30).As one can see, they may contain a source of local-matter creation and contribute to the anisotropic energy-momentum tensors T µν and Y αβ of the horizontal and vertical spaces.Hence, the energy-momentum tensor T µν includes additional information of the action of the local anisotropy of matter fields.Y αβ , on the other hand, is an object with no equivalent in Riemannian gravity, and it incorporates more information of intrinsic anisotropy, which is produced from the vertical metric structure v αβ , and it includes additional gravitational field in the framework of the osculating tangent bundle.Finally, the energy-momentum tensor Z κ α reflects the dependence of matter fields on the nonlinear connection N α µ , a structure which induces an interaction between internal and external spaces.This tensor is different from T µν and Y αβ , which depend on just the external or internal structure respectively.
Lastly, we introduce a new tensor, quantifying the covariant derivative of the torsion tensor, namely Using this definition we can re-write (28) as: where we have omitted from the above relation the term T γ κγ T β λβ since it is of second order.Hence, taking the covariant derivative of (36) we extract the continuity equation for our theory, namely As expected, and as we discussed in the Introduction, F-R-S gravity, similarly to other Finsler-like models, gives rise to an effective interaction between geometry and matter, which can have interesting cosmological implications.

III. COSMOLOGY
In this section we apply the Finsler-Randers-Sasaki geometry and the oscullating Riemannian framework at a cosmological setup, and using the general field equations we derive the generalized Friedmann equations.Then we will provide specific examples.

A. General case
We consider the usual homogeneous and isotropic Friedmann-Robertson-Walker (FRW) metric with a(t) the usual scale factor and k = 0, ±1 the spatial curvature, and substituting it in (9) we find: For simplicity, in the following we focus on the spatially-flat case k = 0. We consider the energy momentum tensor for a perfect fluid in the horizontal and the vertical space: where ρ m and p m are respectively the energy density and pressure of the matter perfect fluid, while u µ and y α are the velocities of the fluid in the horizontal and vertical space respectively.We notice from relations ( 41) and ( 7), (8), that the energy momentum tensor Y αβ of the vertical space constitutes an anisotropic perturbation of the horizontal energy momentum tensor T µν .In order to proceed, we have to consider an ansatz for A γ and y γ .Firstly, it proves convenient to introduce the following scalars: and as we can see, W 0 represents the time-anisotropic contribution of our space while W 1 ,W 2 ,W 3 express the directional components of the anisotropic contribution.In this work we will focus on the case considering y 2 = y 3 = 0 in order to have only a dependance on the parameter t: with y 0 , y 1 constants, and A 0 (t), A 1 (t) time-dependent functions, in agreement with FRW symmetries.
As we observe, in F-R-S cosmology we obtain extra terms in the Friedmann equations, arising from the richer geometrical structure.In particular, the anisotropic torsion terms quantified by the covector field A µ introduce additional degrees of freedom on the tangent bundle of spacetime, which provide the extra contributions in the Friedmann equations.In the case where the internal geometrical structure disappears, namely when W 0 and W 1 become zero, the above equations recover the standard Friedmann equations.Additionally, as it was discussed above, in the scenario at hand we obtain an interaction between geometry and matter, which is now clear by the form of ( 52), and which in the case W 0 = W 1 = 0 recovers the standard conservation equation too.
We can re-write the Friedmann equations ( 48), (49) in their standard form with H(t) = ȧ(t)/a(t) the Hubble function, and where we have introduced an effective dark energy density and pressure of the form where we have emitted the explicit time-dependence of the various quantities in order to make the notation lighter.Thus, the dark-energy equation-of-state parameter is defined as Additionally, the conservation equation ( 52) can be written as where the interaction term Q is given by Hence, differentiating (48) and inserting into ( 49) using ( 59) we also acquire One can now clearly see that in the scenario of F-R-S gravity and cosmology we obtain an interaction between geometry and matter, and therefore an interaction between the effective dark energy and matter sectors.Such an interaction is common in Finsler-like cosmologies [27,51,53,59,70] and it is very interesting since interacting cosmologies [71][72][73][74][75] are known to have many advantages, including solving the coincidence problem [76,77] as well as alleviating cosmological tensions [78,79].

B. Specific model
For completeness, in this subsection we will examine a specific model.As one can see from the general Friedmann equations ( 48), (49), or equivalently from the effective dark-energy sector ( 55), (56), the appearance of the arbitrary functions W 0 (t) and W 1 (t), i.e. of A 0 (t), A 1 (t), makes the resulting cosmological phenomenology very capable.The only point that one should be careful is that the argument of the square root in L(t) should be positive, and thus y 1 should be suitably smaller than y 0 .
We present the evolution of Ω m (z) and Ω DE (z) in the upper graph of Fig. 1.As we can see, we can recover the universe thermal history, i.e. the succession of matter and dark energy epochs.Moreover, in the middle graph Fig. 1 we show the evolution of the corresponding effective dark-energy equation-of-state parameter w DE (z) according to (57).For this specific example w DE lies in the quintessence regime.Nevertheless, note that according to the form of ( 55), (56), one could have other scenarios, in which w DE can be phantom-like, or experience the phantom divide crossing during the evolution.Finally, for completeness, in the lower graph of Fig. 1 we depict the deceleration parameter q, defined as q = −1 − Ḣ H 2 .As we can see, the transition from acceleration to deceleration happens at z tr ≈ 0.7, in agreement with observations.Upper graph: The effective dark energy density parameter ΩDE (green-solid), and the matter density parameter Ωm (blue-dashed), as a function of the redshift z, for F-R-S cosmology under the ansätze ( 46), (47), with y 0 = 1, y 1 = 10 −6 in units where κ = 1.Middle graph: The corresponding dark-energy equation-of-state parameter wDE.Lower graph: The corresponding deceleration parameter q.We have set the initial conditions ΩDE(z = 0) ≡ ΩDE0 ≈ 0.69 [80].

IV. DISCUSSION AND CONCLUDING REMARKS
In this work we presented for the first time a Friedmann-like construction in the framework of an osculating Finsler-Randers-Sasaki geometry, building the corresponding gravitational theory and applying it at a cosmological setup.In particular, we considered a vector field in the metric on the total structure of a Lorentz tangent bundle, which depends on the position coordinates.In this approach, the curvatures of horizontal and vertical spaces, the extra contributions of torsion, non-linear connection and the vector field, that depend on x and y(x), provide an intrinsic richer geometrical structure, with additional degrees of freedom, that lead to extra terms in the field equations.
Applying these modified field equations at a cosmological setup, considering explicit ansätze for the F-R-S metric functions, we extracted the generalized Friedmann equations for the horizontal and vertical parts of R-S spacetime, in which we have the appearance of extra terms, which can be collected to build an effective dark energy sector.Hence, in the framework of F-R-S geometry and gravity, we obtain an effective dark energy density and pressure arising from the richer underlying structure of the tangent bundle.Additionally, as it is common in Finsler-like constructions, we acquire an effective interaction between the matter and the geometrical sectors, and in particular of the extra F-R-S degree of freedom and matter energy density.
We elaborated the generalized Friedmann equations numerically for a specific model, replicating the thermal evolution of the universe, encompassing distinct matter and dark energy epochs.Our analysis revealed that the dark-energy equation-of-state parameter could occupy the quintessence or phantom regime, or undergo a phantom-divide crossing during the evolution.
Several crucial steps remain to be comprehensively explored.Firstly, a thorough investigation into cosmological applications is imperative, involving data confrontation from Type Ia Supernovae (SNIa), Baryon Acoustic Oscillations (BAO), and Cosmic Microwave Background (CMB) observations.Additionally, one could consider different and more complicated ansätze for the F-R-S metric functions.Another interesting subject is the investigation of spherically symmetric and black hole solutions in the theory at hand.These essential and intriguing inquiries are reserved for future research projects.
where we have also used the Leibniz rule.Applying the Stokes theorem again and eliminating the new boundary terms, we find As a last step, the matter part of the action yields Finally, combining equations (A1)-(A6), (A8), (A10), (A11) and setting ∆K = 0, we result to the field equations, namely where we have defined the "energy-momentum tensors" where δ µ ν and δ α β are the Kronecker symbols.

Appendix B: Friedmann equations
In this Appendix we show how the general field equations ( 28)- (30), under the cosmological metric (38) and (39) give rise to the Friedmann equations on the horizontal and vertical space.
Let us start from the horizontal space.First we will calculate the trace of the torsion T β νβ that is required inside (35), (36).From the torsion definition (25) we have that: If we substitute the non-linear connection (17) and the components of L α ακ from ( 14) we find: where h is the trace of the metric h αβ from (8).From the above relation we can see that the torsion tensor is of first order in terms of the weak metric h αβ .
As a last step we take the time component of (36) and its trace, inserting all the above expressions, and substituting the FRW metric (38).After some straightforward manipulations we obtain two equations, namely ȧ a where ρ m and p m are respectively the energy density and pressure of the perfect fluid energy-momentum tensor T µν .Lastly, inserting (B11),(B12) we obtain In order to extract the cosmological equations arising from the vertical space we begin from ( 7) and (8).From those we can calculate : v αβ (x, y(x)) = g αβ (x) − h αβ (x, y(x)), (B17) where h αβ (x, y(x)) = 1 L A γ (y β (x)g αγ + y γ (x)g αβ + y α (x)g βγ ) + 1 L 3 A γ y α (x)y β (x)y γ (x).(B18) The vertical energy-momentum tensor is: Fig. 1.Upper graph: The effective dark energy density parameter ΩDE (green-solid), and the matter density parameter Ωm (blue-dashed), as a function of the redshift z, for F-R-S cosmology under the ansätze (46),(47), with y 0 = 1, y 1 = 10 −6 in units where κ = 1.Middle graph: The corresponding dark-energy equation-of-state parameter wDE.Lower graph: The corresponding deceleration parameter q.We have set the initial conditions ΩDE(z = 0) ≡ ΩDE0 ≈ 0.69[80].