Dark energy and accelerating cosmological evolution from osculating Barthel–Kropina geometry

Finsler geometry is an important extension of Riemann geometry, in which each point of the spacetime manifold is associated with an arbitrary internal variable. Two interesting Finsler geometries with many physical applications are the Randers and Kropina type geometries. A subclass of Finsler geometries is represented by the osculating Finsler spaces, in which the internal variable is a function of the base manifold coordinates only. In an osculating Finsler geometry, we introduce the Barthel connection, with the remarkable property that it is the Levi–Civita connection of a Riemannian metric. In the present work we consider the gravitational and cosmological implications of a Barthel–Kropina type geometry. We assume that in this geometry the Ricci type curvatures are related to the matter energy–momentum tensor by the standard Einstein equations. The generalized Friedmann equations in the Barthel–Kropina geometry are obtained by considering that the background Riemannian metric is of Friedmann–Lemaitre–Robertson–Walker type. The matter energy balance equation is also derived. The cosmological properties of the model are investigated in detail, and it is shown that the model admits a de Sitter type solution and that an effective dark energy component can also be generated. Several cosmological solutions are also obtained by numerically integrating the generalized Friedmann equations. A comparison of two specific classes of models with the observational data and with the standard Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM model is also performed, and it is found that the Barthel–Kropina type models give a satisfactory description of the observations.


Introduction
Finsler geometry [1][2][3][4][5][6] is a well-studied generalization of the Riemann geometry which is constructed from the general line element ds = F x 1 , x 2 , . . . , x n ; dx 1 , dx 2 , . . . , dx n = F(x, y), where F(x, y) > 0, ∀y = 0, is a function on the tangent bundle T (M), and homogeneous of degree 1 in y. In fact, however, Finsler geometry is not just an extension of Riemannian geometry, but it is Riemannian geometry without the quadratic line element restriction [7]. If in Riemannian geometry F 2 = g I J (x)dx I dx J [8], in a Finsler geometry the differential line element ds at x in general is a function of both x and y, and it is defined according to √ dx · dx = g AB (x, y)dx A dx B 1/2 , respectively. In the present paper, capital Latin indices take values in the range (0, 1, 2, 3), while small Latin letters take values in the range (1, 2, 3). Finsler spaces have a much richer mathematical structure than their Riemannian counterparts. In a Finsler space, one can define three kinds of curvature tensors R κ νλμ , S κ νλμ , P κ νλμ , and five torsion tensors [5], which shows that Finsler geometry indeed has a higher degree of complexity than Riemann geometry.
Finsler spaces have a large number of physical applications. Randers [9] unified theory of gravity and electromagnetism, initially constructed in the framework of fivedimensional general relativity, proved to be an example of a Finsler geometry, with ds = (α + β)du, where α = g I K (x)y I y K 1/2 , and β = b I (x)y I , where y I = dx I /du. Quantum mechanics in its hydrodynamical formulation can also be interpreted in terms of Finslerian geometry [10][11][12][13].
The Einstein gravitational field equations, G μν = R μν − (1/2)Rg μν = κ 2 T μν , where R μν is the contraction of the Riemann curvature tensor, R μν = R κ μκν , R is the Ricci scalar R = R μ μ , T μν is the matter energy-momentum tensor, and κ 2 = 8π G/c 4 is the gravitational coupling constant, are the mathematical representation of one of the most beautiful and successful physical theories ever proposed: general relativity. However, general relativity has recently begun to face a number of strong theoretical challenges. The data obtained from the observations of the type Ia supernovae in [14][15][16][17][18][19][20] suggest that the universe is at present in an accelerated, de Sitter type expansionary phase. These observations have led to many observational and theoretical works trying to under-stand the present-day cosmological dynamics (for a review of the cosmic acceleration problem see [21]). The Planck satellite studies of the cosmic microwave background [22], together with the investigations of baryon acoustic oscillations [23][24][25], also confirmed the accelerating expansion of the universe. But to explain these remarkable discoveries, we need a fundamental change in our present understanding of the gravitational interaction. The simplest explanation for the accelerated expansion is obtained by reintroducing in the gravitational field equations the cosmological constant proposed by Einstein in 1917 [26] to obtain a static model of the cosmos. The Einstein field equations together with the cosmological constant provide very good fits to all the observational data, but with the very high price of introducing a parameter whose physical (or geometric) nature is not (yet) known (for detailed discussions on the cosmological constant problem see [27][28][29]. Hence, in order to solve some of the theoretical problems of cosmology without resorting to the cosmological constant, the existence of a new essential component of the universe, called dark energy, was postulated. Many dark energy models have been proposed (for extensive reviews see [30][31][32][33][34][35]). Perhaps the simplest dark energy model can be constructed by using a single scalar field, φ, in the presence of a self-interaction potential V (φ), with the gravitational action given by where R is the Ricci scalar, and M p denotes the Planck mass. The dark energy models obtained in this way are called quintessence models [36][37][38][39][40].
The dark matter problem represents another fundamental problem in present-day astrophysics and cosmology (see [41][42][43] for detailed reviews of the recent results on the search for dark matter, and for its properties). The presence of dark matter at galactic and extragalactic scales is required to obtain an explanation of two basic astrophysical/astronomical observations, the behavior of the galactic rotation curves and the virial mass problem in clusters of galaxies, respectively. The observations of the galactic rotation curves [44][45][46][47] convincingly show that neither Newtonian gravity nor Einstein's general relativity can describe galactic dynamics. To explain the properties of the galactic rotation curves and to solve the missing mass problem in clusters of galaxies, it is necessary to postulate the existence of a dark (invisible) component of the universe that interacts only gravitationally with baryonic matter, and which resides in a spherically symmetric halo around the galaxies. Usually, dark matter is described as a cold, pressureless cosmic fluid. Many candidates have been proposed for the dark matter particle, including weakly interacting massive particles (WIMPs), axions, neutrinos, gravitinos, and neutralinos (for extensive reviews of the dark matter particle candidates see [48][49][50][51]). Dark matter particles interact with normal baryonic matter very weakly; however, their interaction cross sections are expected to be nonzero. Thus, their direct experimental detection in laboratory experiments may be possible.
Therefore, in the present-day view of the cosmological evolution, the local dynamics and the global expansion of the universe are dominated by two major components, dark energy and cold dark matter, respectively, with baryonic matter playing an insignificant role in late-time cosmological evolution. The simplest theoretical model which can fully explain the late-time accelerating de Sitter expansion is obtained by reintroducing in the Einstein gravitational field equations the cosmological constant . The extension of the Einstein gravitational field equations represents the theoretical foundation of the currently dominant standard cosmological paradigm, the cold dark matter ( CDM) model, with cold dark matter also playing a fundamental role. Despite its theoretical simplicity, the CDM model fits the cosmological observations very well [52][53][54][55].
However, apart from its theoretical and fundamental aspects, the CDM model is also confronted with some interesting (and as yet unsolved) observational problems. The most important of these problems is the "Hubble tension," which originated from the important differences obtained between the numerical values of the Hubble constant, H 0 , as determined by the Planck satellite from the measurements of the cosmic microwave background (CMB) [55], and the values measured directly by using cosmological and astrophysical observations in the local universe [56][57][58]. To illustrate the differences, the SH0ES determinations of H 0 give the value H 0 = 74.03 ± 1.42 km/s/Mpc [56], while the analysis of the CMB, originating in the early universe, gives, via the Planck satellite data, H 0 = 67.4 ± 0.5 km/s/Mpc [54], a value differing from the SH0ES result by ∼ 5σ .
Hence, the search for alternative explanations for the cosmological dynamics and the nature of the two mysterious major components of the universe is a critical task for presentday theoretical physics. One of the attractive possibilities for solving the gravitational puzzles is to go beyond the framework of the Riemannian geometry on which general relativity is constructed, and look for more general geometries that could describe the intricate nature of the gravitational phenomena. In this direction, one very promising candidate is Finsler geometry.
Horváth [59] and Horváth and Moór [60] were the first to attempt to formulate a relativistic theory of gravitation by using Finsler geometry, with their work later extended in [61] and [62], respectively. The physical and cosmological implications of the Finsler geometry have been investigated from different points of view in . In particular, in [100], the cosmological evolution in an osculating point Barthel-Randers type geometry was considered. In this type of geometry, each point of the spacetime manifold is associated with an arbitrary point vector field. For the Barthel-Randers geometry, the connection is given by the Levi-Civita connection of the associated osculating Riemann metric. This Finsler type geometry is assumed to describe the physical properties of the gravitational field via the standard Einstein equations, as well as the cosmological dynamics. The generalized Friedmann equations in the Barthel-Randers geometry were obtained by considering that the background Riemannian metric in the Randers line element is of Friedmann-Lemaître-Robertson-Walker type. The matter energy balance equation was obtained, and it was investigated from the perspective of the thermodynamics of irreversible processes in the presence of particle creation. The cosmological properties of the Barthel-Randers model were investigated in detail, and it was shown that the model permits a de Sitter type solution, and that an effective cosmological constant can also be generated. Several exact cosmological solutions were also obtained, and a comparison of three specific models with the observational data and with the standard CDM model was performed by fitting the observed values of the Hubble parameter. It turns out that the Barthel-Randers models give a satisfactory description of the observations. Alternative mathematical formulations within a Finsler geometric framework are also possible. For example, in [91] the cosmological implications of scalar-tensor theories that are effectively obtained from the Lorentz fiber bundle of a Finsler-like geometry were investigated. In this approach, it is assumed that in the Finsler space one can define a nonlinear connection with local components N (α) μ x ν , φ (β) , where x ν , φ (β) are the local coordinates, with x ν , ν = 0, . . . 3, the coordinates on the local manifold, and φ (β) , β = 1, 2 the coordinates on the fiber. The nonlinear connection uniquely splits the total space T E into a horizontal distribution H E and a vertical distribution V E, so that T E = H E ⊕ ⊕V E. The metric tensor can then be constructed as where for g μν a Lorentzian signature (−, +, +, +) was adopted, with the fiber variables φ (α) playing the role of internal variables. The gravitational field equations are obtained from a variational principle, with the action constructed as (2) . For the Lagrangian density, two forms were considered, the first being L G =R − V (φ)/φ, wherẽ R = R −2 φ/φ +∂ μ φ∂ μ φ/2φ 2 is the curvature for the case of a holonomic basis, and V (φ) is the scalar field potential. A second possible choice considered in [91] is L G =R, wherē R is the curvature of a non-holonomic basis. Matter can also be added to the gravitational action, and thus the total action becomes where L M is the matter Lagrangian. By using both choices of action, the cosmological implications of the model were investigated, and it was shown that an effective dark energy sector does appear from the geometric structure of the Finsler model. Moreover, an explicit interaction between the matter and dark energy sectors is generated. As applied to the cosmological evolution of the universe, it was shown that one can obtain a sequence of successive matter-and dark energydominated epochs, with the parameter of the equation of state of the dark energy being quintessence-like, phantom-like, or experiencing the phantom-divide crossing. An exponential de Sitter solution can also be obtained. A generalized scalar-tensor theory obtained from vector bundle constructions was considered in [98], where its kinematic, dynamical, and cosmological consequences were studied. The theory is characterized by a mathematical structure constructed over a pseudo-Riemannian spacetime base manifold, together with a fiber structure with two scalar fields. Hence, one obtains a six-dimensional vector bundle endowed with a nonlinear connection as the geometric arena for the description of the gravitational interaction. The geodesics and the Raychaudhuri and field equations are obtained using both the Palatini and the metric approach. For the cosmological metric, we adopt the expression [98] (1) ⊗ δφ (1) + δφ (2) ⊗ δφ (2) , (4) and the corresponding Friedmann equations become [98] where and respectively, where W + is the contribution coming from the non-diagonal components of the field equations. Hence, the Finsler geometric structure generates extra terms in the modified Friedmann equations corresponding to an effective dark energy sector, also leading to an interaction of the dark mat-ter sector with the metric. The parameter of the dark energy equation of state can take values in either the quintessence or the phantom regime, or that indicate the phantom-divide crossing.
It is the goal of the present paper to extend the investigations of the cosmological properties of the Finsler spaces, as initiated in [100], by considering a systematic investigation of another Finsler type geometry which is based on the Kropina metric [101,102], a particular class of the general (α, β) metrics with F = α 2 /β. The Randers metric also belongs to the class of (α, β) metrics. We assume that the true description of the gravitational interaction can be accomplished by using the mathematical properties of the Kropina spaces in a mathematical framework introduced by Barthel [103,104], in which one can consider a Finsler space as an n-dimensional point space which is locally Minkowskian (a space that is flat, homogeneous, and anisotropic) but, in general, is not locally Euclidean. Note that a general Finsler space is both inhomogeneous and anisotropic.
The starting point of the theory of the gravitational interaction based on the point Finsler spaces is the assumption that the gravitational field can be represented by a Riemannian metric g(x), satisfying the Einstein gravitational field equations. The next step in the construction of a Finsler type gravitational theory is the non-localization (anisotropization) of the background geometry, by attaching to each point x = x I , I = 0, 1, 2, 3, an internal variable y = y J , J = 0, 1, 2, 3. Under the assumption that y is a vector, the nonlocal gravitational field can be described by using a Finsler type geometry F 4 or, alternatively, by the geometry of a general vector bundle. In nonlocal Finslerian geometry, the metric tensor depends on both the local coordinates x and the vector y, so thatĝ =ĝ(x, y).
However, for many realistic physical processes, one can assume that the variable y is a function of the position y = Y (x). Thus, the Finslerian metric becomesĝ = g (x, Y (x))). Hence, in this particular type of physical model, the Finslerian metric tensorĝb is a function only of x. In this way, we define a manifold called the osculating Finsler manifold. The geometric properties of the osculating Finsler space are described by using the Barthel connection [103,104]. In the case of the Kropina geometry (as well as in all (α, β) geometries) we have the remarkable result that the Barthel connection is the Levi-Civita connection of the Riemannian metricĝ(x, y(x)). In the present work we systematically and rigorously investigate the cosmological implications of the Barthel-Kropina geometry by adopting for the background Riemann metric the Friedmann-Lemaître-Robertson-Walker (FLRW) form. In the Barthel-Kropina-FLRW geometry we obtain the generalized Friedmann equations of the model. The properties of the cosmological equations are investigated in detail, and several evolutionary scenarios are constructed, describing different types of dynamics. The model admits a de Sitter type solution, describing an accelerating universe but in the presence of nonzero matter density. Different cosmological solutions corresponding to particular choices of the coefficients of the one-form β are also considered, and their behavior is compared with the observational data as well as with the predictions of the standard CDM model. An effective dark energy model can also be generated from the generalized Friedmann equations, and its properties are analyzed in detail. Decelerating solutions can also be obtained, and they may provide, for example, alternative descriptions to the radiation-dominated phase of standard cosmology. The overall comparison of the considered models with the observational data indicates that they provide a satisfactory description of the cosmological dynamics, at least on a qualitative level.
The present paper is organized as follows. In Sect. 2, we review the necessary elements and definitions of the Finsler geometry and point and osculating Finsler spaces. The Barthel connection is also introduced in an intuitive way. General definitions of the curvature tensors are also presented, and the specific form of the Einstein gravitational field equations is postulated. The cosmological framework involving the definitions of the Riemann metric, the one-form β, and the thermodynamic parameters is constructed in Sect. 3. The Finslerian metric tensor coefficients, the Christoffel symbols, and the curvature tensors are obtained for the considered geometry. Based on this, geometric results are derived for the generalized Friedmann equations describing the cosmological evolution in the Barthel-Kropina-FLRW geometry. Particular cosmological models are considered in Sect. 4. It is shown that the model admits a de Sitter type solution, in the presence of a nonzero matter density. Cosmological models corresponding to two specific choices of the function η, giving the only nonzero coefficient of the oneform β, are considered in detail. The predictions of the model with a specific dependence on the Hubble function of β are compared with observational data and the standard CDM model. The possibility of obtaining an effective dark energy model equivalent to a time-varying cosmological constant is analyzed in Sect. 5. We discuss our results and presents conclusions in Sect. 6. The calculations of all the geometric quantities are presented in five appendices. In Appendix A, the details of the computation of the Finsler metric tensor are explicitly shown. The expressions of the Christoffel symbols in the Barthel-Kropina-FLRW geometry are derived in Appendix B. The components of the Ricci tensors are obtained in Appendix C. The Ricci scalar is computed in Appendix D. Finally, the explicit computations leading to the Einstein gravitational field equations and to the generalized Friedmann equations are presented in Appendix E.

A brief review of the Finsler, (α, β), and Barthel-Kropina geometries
In the present section we will briefly introduce the basics of the Finsler geometry, the (α, β)-metrics, the Barthel connection, and the Barthel-Kropina geometry.
It is important to emphasize that Finsler geometry already appears in classical Newtonian mechanics when dissipative effects are present. If the equations of motion of a dynamical system, defined on an n-dimensional differentiable manifold M, can be obtained from a regular Lagrangian L via the Euler-Lagrange equations, given by where F i are the external forces, the above Euler-Lagrange equations are equivalent to a system of second-order differential equations of the form where each function G i x j , y j , t is C ∞ in in the neighborhood of some initial conditions (x) 0 , (y) 0 , t 0 . Equation (10) can naturally be interpreted as describing geodesic motion in a Finsler space.

Finsler geometry and particular Finsler spaces
A basic assumption of the present-day theoretical physics is that space and time are unified in a single structure (the spacetime) that can be described mathematically as a fourdimensional differentiable manifold M, endowed with a pseudo-Riemannian metric tensor g I J , where I, J, K . . . = 0, 1, 2, 3. On the world line of a standard clock, the spacetime interval between two events x I and x I + dx I is given, according to the chronological hypothesis, by ds = g I J dx I dx J 1/2 [105,106]. But from a mathematical point of view, one can go beyond the Riemannian mathematical structures of the spacetime manifold. One of the most important metric extensions of the Riemann geometry is the geometry introduced by Finsler [1][2][3][4][5][6].
From a general mathematical point of view, Finsler spaces are metric spaces in which the interval ds between two nearby points x = (x I ) and x + dx = (x I + dx I ) is given by where F, the Finsler metric function, is positively homogeneous of degree 1 in dx, and thus has the property To permit the use of local coordinates in mathematical computations, the Finsler metric function F is usually written in terms of the canonical coordinates of the tangent bundle (x, y) = (x I , y I ), where y = y I ∂ ∂ x I is any tangent vector y at x. Then the Finsler metric tensorĝ I J is defined aŝ thus allowing us to write Eq. (11) as dŝ 2 =ĝ I J (x, y) y I y J . Riemann spaces are particular cases of Finsler spaces, corresponding toĝ I J (x, y) = g I J (x), y I = dx I , and ds 2 = g I J (x)dx I dx J , respectively. Another important geometric quantity, the Cartan tensor C(x, y), is defined according tô

Randers, Kropina and general (α, β) metrics
The Randers space, which has many applications in physics, is a special kind of Finsler space [9], with the Finsler metric function given by where g I J (x) is the fundamental tensor (metric) of a Riemannian space, and A I (x)dx I is a linear one-form on the tangent bundle T M. Kropina [101,102] considered Finsler spaces with metrics given by These early results were generalized by Matsumoto [107,108], who introduced the notion of the (α, β) metrics. A Finsler metric function F(x, y) is called an (α, β) metric when F is a positively homogeneous function F(α, β) of the first degree in two variables α (x, y) = g I J (x)dx I dx J 1/2 and β (x, y) = A I (x)y I , respectively.
In the following we will assume that α is a Riemannian metric; that is, it has the basic properties of being nondegenerate (regular) and positive-definite. The Randers and the Kropina metrics belong to the class of (α, β) metrics, with F = α + β for the case of the Randers metric and can introduce general (α, β) metrics having the form where s = β/α, Denoting L = F 2 /2, we obtain for the fundamental metric tensor of the (α, β) space the expression where and the indices α, β of L indicate partial differentiation with respect to α and β. Alternatively, one can compute the components of the Finslerian metric tensor by using the formula obtained in [109]. Let α = g I J (x)y I y J , β = A I (x)y I , where = ±1, and F is given by Eq. (17). Then the Hessian g can be obtained aŝ (20) where α I := α y I and

The Barthel connection and the osculating Finsler spaces
In the following we will briefly review the fundamental mathematical concepts on which the gravitational applications of the Kropina geometries are based, namely, the Barthel connection and the osculation of Finsler spaces.

The Barthel connection
Let us now assume that (M, F) is a Finsler space and Y (x) = 0 is a vector field defined on M. We can now introduce a specific structure An important class of Finsler spaces is the point Finsler spaces. Based on the definition introduced by Barthel [103,104], we consider in the following a Finsler space as an ndimensional point space, which is locally Minkowskian and, in general, not locally Euclidean. A general Finsler space is both inhomogeneous and anisotropic, while a Minkowski space is flat and homogeneous, but anisotropic. Accordingly, we shall call the Finsler n-space a Barthel-Finsler space or, for short, a point Finsler space.
Given a point vector field Y I (x) and a Finsler metric tensorĝ(x, y), one can construct the absolute differential of Y according to the definition [110] Barthel by writing the expressions in the second term of Eq. (24) as and defining the Barthel connection as [110] b The Barthel connection depends on the field to which it is applied, leading to a situation that is very different from the properties of the connection in Riemannian geometry. For anisotropic metrics, all properties can depend on the direction, and in the case of the Barthel connection, the dependence is only on the direction of the field, and not on its magnitude. The Barthel connection is the simplest one that preserves the metric function by the parallel transport. Moreover, for Finsler vector fields, depending on both x and y, it allows a natural passage to the standard Cartan geometry of Finsler spaces. Therefore, we can consider the Barthel connection as the connection of a point Finsler space.
From the general theory of the Finsler connections [5,6] it follows that these connections, unlike the usual Levi-Civita connection of a Riemannian metric, or more general affine connections, do not live on the base manifold M but on the total space of the tangent bundle. This essential difference sometimes induces major dissimilarities in the geometric theory of Riemannian and Finsler manifolds.
To fill in this gap, one of the possibilities is to evaluate the Finsler connections on a nowhere-vanishing tangent vector field (x, Y (x)), provided the topology of the base manifold M allows the existence of such vector fields. By evaluating the fundamental tensor g i j (x, y) of (M, F) at (x, Y (x)), one obtains a Riemannian metric g Y on M, with its own Levi-Civita connection. Likewise, by evaluating the Cartan connection of (M, F) at (x, Y (x)) one obtains an affine connection on M, in general different from the Levi-Civita connection of g Y , which is the Barthel connection already introduced above. Barthel himself arrived at this connection when studying the parallel transport on M.

The osculating Riemannian metric
In this subsection we will introduce the osculating Riemannian metric associated with a Finsler metric (M, F). The concept of osculating Riemannian spaces of Finsler spaces was introduced by Nazim [111], and it was later studied extensively in [112]. In simple terms, the osculation means that a rather complicated structure (Finsler geometry, Finsler connection) has an associated simpler structure (a Riemannian metric, an affine, or linear connection), with the simpler structure assumed to approximate the former in some sense. Based on this link, one can obtain results on the more complicated structure.
If we fix such a local section Y of π M : T M → M, all geometric objects defined on the manifold T M can be pulled back to M. Sinceĝ I J • Y is a function on U , we can definê The For the osculating Riemannian metric (28), the Christoffel symbols of the first kind are defined according tô By using the derivative law of composed functions we obtain In the case of a general (α, β) metric for the Cartan tensor, we obtain respectively, where we have denoted Hence, if Y is a nonvanishing global section of T M, so that Y (x) = 0, ∀x ∈ M, one can always define the osculating Riemannian manifold (M,ĝ i j ).
In the case of an (α, β) metric, let us consider the vector field Y = A, where A I = g I J A J . Since the vector field A is globally nonvanishing on M, it follows that β has no zero points. Hence, we can introduce the A-osculating Rieman- Letã be the length of A with respect to α. Hence, Therefore, the A-osculating Riemannian metric becomeŝ Moreover, we have and consequently from Eq. (30) we obtain the important result thatĈ I J K (x, A) = 0. On other hand, for Y = A, we obtain Hence, for a Finsler space with (α, β)-metric the linear A-connection with A I = (g I J A J ) (the Barthel connection), is the Levi-Civita connection of the A-Riemannian metric.

The curvature tensor and its contractions
The curvature tensor of an affine connection with local coefficients A BC (x) is given by The Barthel connection with local coefficients b A BC (x) is an affine connection, and hence its curvature tensor must be given by the above formula, with A BC (x) = b A BC (x) . As we have already mentioned, in the case of the Kropina metric F = α 2 /β, the Barthel connection coincides with the Levi-Civita connection of the osculating metricĝ AB (x) = g AB (x, A(x)), where A I (x) are the components of β, and g AB is the fundamental tensor of F. Hence, since b A BC = γ A BC , whereγ A BC are the Levi-Civita coefficients, we obtain for the curvature tensors the expressionŝ and respectively, where A, B, C, D, E ∈ {0, 1, 2, 3}. For further details on the definitions of the affine connections and the curvature tensors, see [100] and references therein. The contractions of the curvature tensor lead to the generalized Ricci tensor and Ricci scalar, respectively given bŷ andR respectively.

Cosmological evolution in Barthel-Kropina geometry
In the present section, we will consider the general framework of the cosmological evolution in a Barthel-Kropina type geometry. We will begin by formulating a theoretical framework for the geometric and physical quantities of our model, in which the geometric and physical quantities are defined. As a next step, we will consider the computation of the geometric quantities (metric, Christoffel symbols and curvatures) in the adopted geometry. Finally, the generalized Friedmann equations describing the cosmological evolution are obtained.

Metric and thermodynamic quantities
In order to systematically construct the cosmology of the Barthel-Kropina geometry, and to facilitate the physical interpretation of the results, we adopt the following assumptions.

The Riemannian metric g I J (x) in α is given by the Friedmann-Lemaître-Robertson-Walker (FLRW) metric,
where we have adopted a coordinate system with (x 0 = ct, The space-like components of A vanish. The cosmological principle, together with the diagonal nature of the FLRW metric, imposes another mathematical condition on the one-form β. More exactly, we require that the space-like components of A vanish so that A 1 = A 2 = A 3 = 0. If this condition were not satisfied, we could perform a spatial rotation, thus obtaining a preferred direction, for example, in the x coordinates. But such a cosmological model would contradict the observationally well-confirmed, large-scale spatial isotropy of the universe. Therefore, we assume that the vector A has only one time-like independent component, A 0 x 0 . Hence, we represent the one-form field β as where η x 0 is an arbitrary function of time to be determined from the gravitational field equations. A frame comoving with matter does exist. Similarly to the standard general relativistic case and its Riemannian geometry, we assume that in the Barthel-Kropina geometry one can introduce a comoving frame, in which the motion of all observers take place with the Hubble flow, defined by the Riemannian metric g I J (x). Thermodynamic properties. We assume that the thermodynamic properties of the cosmological matter can be described by two physical quantities only, the energy density ρc 2 and the thermodynamic pressure p, respectively, defined in the standard way. From assumptions c and d it follows that the only nonvanishing components of the matter energymomentum tensorT AB areT 0 0 = ρc 2 ,T 00 =ĝ 00T 0 0 and

The gravitational field equations We postulate that the Einstein gravitational field equations can be formulated in a Barthel-Kropina geometry aŝ
where κ 2 = 8π G/c 4 is a constant, and G and c are the Newtonian gravitational constant and the speed of light, respectively.T B D is the matter energy-momentum tensor, constructed with the help of the usual thermodynamic quantities as defined in the standard Riemann geometry, and of the Finslerian metric tensorĝ B D .
In the present investigation, we have adopted for the Einstein equations (43) the simplest possible mathematical form, guided by the analogy with standard general relativity. Moreover, the adopted field equations have a well-defined Riemannian limit, corresponding to the caseĝ(x, y) → g(x), when we directly recover standard general relativity, without the cosmological constant term. There are many extensions of general relativity, and other forms of the field equations, inspired by these extensions, are also possible. For example, in the f (R) modified gravity theory, introduced in [113], the gravitational Lagrangian is an arbitrary function of the Ricci scalar R. One could assume in the framework of a Finsler variational principle that the field equations follow from a Finsler type Lagrangian of the form f R . Alternatively, one could construct the field equations in a less rigorous approach directly by substituting in the standard general relativistic equations the Riemannian quantities with their Finslerian counterparts. Modified gravity theories with curvature-matter coupling [114][115][116] could also be extended to a Finslerian geometric framework. Riemannian geometries with torsion [117,118] or non-metricity [119][120][121] can also represent a source of inspiration for obtaining their Finslerian analogues.
Hence, based on the above assumptions, the metric properties of Barthel-Kropina-FLRW model are given by

Geometric quantities in Barthel-Kropina-FLRW geometry
As a first step in constructing a geometric theory in the Barthel-Kropina geometry, we need to obtain the components of the Finslerian metric tensor (the Hessian), which are given by (for the full computational details see Appendix A) Next, we need to obtain the Barthel-Kropina-FLRW Finsler metric, using the geometric properties introduced in the previous section. We immediately obtain For i, j ∈ {1, 2, 3} we find Hence, we obtain the nonvanishing components of the Barthel-Kropina-FLRW metric tensorĝ I J as given bŷ The inverse components of the Barthel-Kropina-FLRW metric tensorĝ I J can be obtained aŝ For the Christoffel symbols associated with the Barthel-Kropina-FLRW metric, we obtain (the detailed calculations are presented in Appendix B) where H = a a . The nonvanishing components of the Ricci tensor are given by (see Appendix C for the calculational details) and respectively. For the Ricci scalarR =ĝ I JR I J we find (see Appendix D) Finally, we obtain the Einstein tensor componentŝ and respectively. For the explicit calculations of these components see Appendix E.

The generalized Friedmann equations
Once the geometric quantities in the Barthel-Kropina-FLRW geometry are known, we can write the full set of Einstein gravitational field equations relating geometry and matter. By taking into account our assumptions about the physical nature of the cosmological quantities, the Einstein field equationsĜ 00 = 8π G/c 4 ĝ 00 ρ andĜ ii = − 8π G/c 4 ĝ ii p lead to the system of the generalized Friedmann equations, given by and respectively. By replacing the term −3 η 2 with the help of Eq. (52), Eq. (53) takes the alternative form

General relativistic limit and the energy balance equation
A first interesting property of the cosmological equations of the Barthel-Kropina-FLRW model is that in the limit η → ±1/a, β = (1, 0, 0, 0), the generalized Friedmann Eqs. (52) and (53) do reduce to the standard Friedmann equations of general relativity, a result which is easy to check. From the Friedmann equations of standard general relativistic cosmology, it follows that the matter energy density ρ satisfies the conservation equatioṅ where the dot denotes the derivative with respect to the time, and H = cH. For p = 0, the conservation equation can be immediately integrated to give where ρ 0 = ρ 0b + ρ 0DM is the present-day matter density of the universe, with ρ 0b and ρ 0DM denoting the present-day density of the baryonic and dark matter, respectively, and we have introduced the redshift variable, defined as 1/z = 1+a. We now introduce the generalized matter density parameter m , defined as By considering a three-component universe filled with ordinary matter, dark matter, and dark energy, respectively, in the standard CDM model, the Hubble function is given by = ρ b /ρ cr , DM = ρ DM /ρ cr and = /ρ cr denote the density parameters of the baryonic matter, dark matter, and dark energy, respectively. Here, ρ cr denotes the critical density of the universe, defined as ρ cr = 3H 2 0 /8π G. An important observational parameter, the deceleration parameter is defined as In the CDM model, the deceleration parameter is obtained as To compare the predictions of the Barthel-Kropina-FLRW cosmological model with observations and with the CDM model, we adopt for the matter and dark energy density parameters the values DM = 0.2589, b = 0.0486, and = 0.6911, respectively [22]. Then, for the total matter density parameter m = DM + b , we obtain m = 0.3089. The present-day value of the deceleration parameter has the value q(0) = −0.5381, indicating that at present the universe is in an accelerating phase.

The Riemannian limit of the Kropina metric
We now consider the limiting case of the Finsler metric. Since in the case of the considered Kropina metric, F = α 2 /β, with α 2 = g I J (x)dy I dy j and β = A I dy I , the standard general relativistic Riemann geometry is recovered in the limit β → α, which gives lim β→α F = α = g I J (x)y I y J . The condition α = β gives immediately g I J y I y J = A I A J y I y J , or lim α→βĝI J = A I A J . This means that the limiting Riemann metric g I J is actually an inner product, and hence it becomes a degenerate Riemann metric in the sense that det |g I J | = 0; that is, rankg I J = n − 1.
Equivalently, by taking the limit β → α in Eq. (A1), giving the metric tensor of the Barthel-Kropina geometry, one can show that lim β→αĝI J = g I J = A I A J . The fact that in the Riemannian limit the Kropina metric leads to a degenerate Riemann metric should not come as a surprise. Indeed, the original Kropina metric is degenerate at the origin of T x M, and by taking the limit β → α, we obtain a degenerate Riemann metric.
This also means that in this limit, the Kropina indicatrix degenerates into the Riemannian unit circle of the degenerate metric g I J . Since lim β→αĝ (x, y) = g I J (x), under the same limit the Y -osculating metricĝ (x, Y (x)) =ĝ(x) tends to g I J (x), which is degenerate, as already explained. Moreover, let us observe that the condition g I J → A I A J is incompatible with the choice g I J = diag 1, −a 2 , −a 2 , −a 2 and b I = (b 0 , 0, 0, 0), respectively, except for the case a 2 = 0, which is obviously meaningless in the present context.

The energy balance equation
We now multiply Eq. (52) by a 3 , and take the derivative of the resulting relation with respect to x 0 . Thus, we obtain After multiplying Eq. (53) by 3a 2 a , we obtain By adding Eqs. (63) and (64), we obtain the energy balance equation in the Barthel-Kropina-FLRW cosmological model as The energy balance equation can be reformulated, using the generalized Friedmann equations, as

Particular cosmological models
In the present section, we will investigate the cosmological implications of the generalized Friedmann equations (52) and (53) (53) is automatically satisfied. Hence, a de Sitter type exponential expansion can take place only in the presence of matter. Thus, we assume a nonzero matter energy density, but a vanishing thermodynamic pressure for the cosmological matter, with p = 0. In the case H = H 0 = constant, Eq. (53) becomes By introducing a new variable u, defined as η = u, we immediately obtain η = du/dx 0 = (du/dη) dη/dx 0 = u(du/dη). Hence, Eq, (67) becomes with the general solution given by where C 1 is an arbitrary constant of integration. Therefore, from we obtain the general solution of Eq. (67) as where C 2 is an integration constant that can be taken as zero without any loss of generality by a rescaling of the time coordinate x 0 . Moreover, in order to obtain a physically consistent solution, we will adopt the minus sign in the above equation. Hence, a de Sitter type exponential expansion of the Riemann metric corresponds to During the de Sitter phase, the matter energy density varies as At the initial moment x 0 = 0, the matter energy density is while for x 0 → ∞, we obtain 8π G c 2 lim In order to obtain a monotonically decreasing matter energy density during the de Sitter phase satisfying the condition ρ x 0 x 0 =0 > lim x 0 →∞ ρ x 0 , the integration constant C 1 must satisfy the condition C 1 > 1/2.

Conservative cosmological evolution
We will now consider conservative Barthel-Kropina-FLRW models, that is, models in which the matter conservation equation is identically satisfied. In this case, the system of generalized Friedmann equations reduces to the Friedmann equations (52) and (53), together with the constraint

Linear barotropic fluid cosmological models
As a first example of a conservative Barthel-Kropina-FLRW cosmological model, we consider the case of a barotropic fluid-dominated universe, with p = (γ − 1)ρc 2 , where γ = constant and 1 ≤ γ ≤ 2. In this case, Eq. (76) immediately gives where ρ 0 is an integration constant. Then, Eq. (52) gives where we have assumed that η is a monotonically decreasing function of time. By taking into account the linear barotropic equation of state, using Eq. (52), the pressure can be obtained as After substituting the above expression of the pressure in Eq. (53), we obtain the equation which can be integrated immediately to obtain the first integral where C is an arbitrary integration constant. Using the expression (79) for η , we obtain By taking the derivative of the above expression of η with respect to x 0 , and equating it with η obtained from the first Friedmann equation as given by Eq. (79), we obtain for a the equation with the general solution given by a x 0 ∝ x 0 2/3γ . Hence, conservative Barthel-Kropina-FLRW cosmological models exactly mimic their general relativistic counterparts, describing the decelerating evolution of a matter-dominated universe satisfying a linear barotropic equation of state. For a dust universe we obtain a x 0 = x 0 2/3 , while for a radiation-dominated universe, the scale factor varies as a x 0 = x 0 1/2 .

Nonconservative Barthel-Kropina-FLRW cosmological models
In the following we consider cosmological models in which the matter density and pressure do not satisfy the conservation Eq. (76), and thus for a linear barotropic equation of state, the expression of the energy density is different from the one given by Eq. (78). However, we will consider that matter still satisfies a linear barotropic equation of state. Hence, by eliminating the pressure from the second Friedmann equation (53) with the help of the equation of state and the first Friedmann equation (52), we obtain which can be immediately integrated to give the first integral where C is an arbitrary constant of integration, and we have used the standard mathematical relation For the matter density we find the expression 8π G c 2 ρ = a 2 η 2 = C 2 η 3γ .
In order to obtain a monotonically decreasing matter density, η must also be a decreasing function of time. In terms of η, the Hubble function is obtained as Eqs. (87) and (88) give the full solution of the Friedmann equations for the cosmological evolution in the Barthel-Kropina-FLRW geometry. Once the functional form of η is specified, the evolution of the scale factor and the other cosmological parameters can be immediately obtained.

η = η(a)
As a first example of a nonconservative Barthel-Kropina-FLRW cosmological model, we consider the case in which η is a function of the scale factor only, η = η(a). Then, Eq. (87) immediately gives the evolution equation which can be integrated to give where x 0 0 is an arbitrary integration constant. In the following we will consider the case of a dust universe with γ = 1. The simplest possible cosmological model can be obtained for η(a) = a 1 a −n , where a 1 and n are constants. Then we immediately obtain where we assume 5n < 2, n < 0.4. For the Hubble function, we find while for the deceleration parameter we obtain the expression q = −5n/2. The universe is in an accelerated expansionary state, and the scale factor has a power law dependence on the time coordinate.

η = η (H)
As another particular class of Barthel-Kropina-FLRW cosmological models, we consider the case in which η is a function of H, the Hubble function, so that η = η(H). Hence, under this assumption, Eq. (87) takes the form At this moment we introduce as the independent variable the redshift z, defined as 1 + z = 1/a. Hence, we obtain In the redshift variable, Eq. (94) becomes In the following we adopt for η(H) the simple form where a 1 , b 1 , and n are constants. We also rescale the Hubble function according to H = (1/c)H = (1/c)H 0 h, where H 0 is the present-day value of the Hubble function. Moreover, we rescale the coefficient a 1 as a 1 → a 1 /H n 0 . We consider the late-time evolution of the universe in its dust phase, with zero pressure, and therefore we take γ = 1. Hence, Eq. (96) takes the form For the matter density parameter, we obtain where we have denoted C 1 = C 2 c 2 ρ c /3H 2 0 ρ 0 . The variation in the normalized Hubble function, the deceleration parameter, and the matter density parameter are presented in Figs. 1, 2, and 3 for four values of n (n = 1, 2, 3, 4) and for different values of a 1 , b 1 , and c 1 . In the figures we have also included the observational data for the Hubble function, together with their error bars [124,125], as well as the predictions of the CDM model.
As one can see from Fig. 1, for the adopted functional form of η, the Barthel-Kropina-FLRW model gives an acceptable description of the observational data for h(z), up to a redshift of around z ≈ 3, and it can also reproduce, at least qualitatively, the predictions of the CDM model. However, with increasing z, at redshifts higher than 3, important deviations from the predictions of the CDM model do appear, with the dimensionless Hubble functions of the CDM model increasing much faster than in Barthel-Kropina-FLRW cosmology.
For the adopted set of the model parameters, the comparison of the deceleration parameter curves, presented in Fig. 2, as predicted by the Barthel-Kropina-FLRW model and by the CDM model indicate the existence of significant differences between the two models, especially at higher redshifts. Both models predict a late-time accelerating behavior, with  At a qualitative level, the Barthel-Kropina-FLRW cosmology can also reproduce the behavior of the total matter density parameter m in the CDM model. At low redshifts 0 ≤ z < 1, the predictions of the two models basically agree. However, the quantitative differences, already observed in the case of the deceleration parameter, do also exist for this physical parameter. In the range 1 ≤ z ≤ 3, the behavior of the matter density parameter strongly depends on the model parameters, CDM, and at higher redshifts, both lower and higher matter density parameters can be obtained, thus leading to matter densities that are different from the predictions of the CDM model.

Dark energy in the Barthel-Kropina-FLRW cosmology
We will now consider that the function η can be represented in a general form as where f is an arbitrary function to be determined from the field equations. In the limit f → 0, η → 1/a, and as we have already seen, we recover the standard Friedmann equations of general relativity. With this representation of η, the generalized Friedmann equations take the form and where and respectively. An effective dark energy term, satisfying the condition p DE = wρ DE , w = constant can be obtained if the function f satisfies the equation The dynamical cosmological dark energy is thus dependent on the scale factor and on the properties of η. In the following we will consider only the late-time evolution of the universe, and hence we take p = 0. To simplify the mathematical formalism, we introduce the dimensional time By introducing the critical density ρ c = 3H 2 0 /8π G, and by defining the matter density parameter as m = ρ/ρ c , from the first Friedmann equation we obtain In terms of the redshift variable 1 + z = 1/a, the system of equations (107)-(109) become For the redshift dependence of the matter density parameter, we find For the deceleration parameter, we obtain The variation in the Hubble function, obtained by numerically solving the system of differential equations Eqs. (111)- (113), is represented for different values of the parameter w of the dark energy equation of state in Fig. 4. As one can see from Fig. 4, the model gives a good description of the observational data, as well as of the CDM model up to a redshift of at least z = 5. For higher redshifts, some differences with respect to the standard cosmology do appear. The numerical results do not show a significant influence of the variation in the parameter of the dark energy equation of state w, with w = 1.85 predicting almost the same evolution as w = 3.05, but depend strongly on the initial conditions for f and u = d f /dτ . The variation in the deceleration parameter as a function of redshift as predicted by the Barthel-Kropina-FLRW dark energy model is presented in Fig. 5. Even though on a qualitative level the behavior of the deceleration parameter of the Barthel-Kropina-FLRW model can reproduce the predictions of the CDM model, significant quantitative differences still appear. In the Finslerian approach, the evolution of the deceleration parameter is strongly dependent on the parameter w of the dark energy equation of state. At higher redshifts, the numerical values of the Barthel-Kropina-FLRW are slightly lower than those obtained from CDM, and almost independent of w, while at low redshifts the behavior of q essentially depends on the numerical value of w. The critical redshift z cr indicating the transition from deceleration to acceleration, given by q (z cr ) = 0, also takes numerical values that are different from the CDM predictions. The critical redshift is also strongly dependent on the initial conditions at the present time for f and its derivative.
The density parameter of the total energy matter content of the universe is represented as a function of the redshift in Fig. 6. In this case there are significant differences, especially at high redshifts, between the predictions of the Barthel-Kropina-FLRW and the CDM model, with the Finsler type cosmological model predicting much higher matter densities at higher redshifts. In the present model, the matter energy density not only depends on the Hubble function, as in standard general relativity, but also includes the contribution from the one-form field β. Since there is no independent conservation of the matter energy-momentum tensor, the increase in the total matter density at high redshifts can also be interpreted as describing a particle creation process via transfer of energy from the η field to ordinary matter, or equivalently, to the increase in the contribution of the dark energy. On the other hand, in the present model we cannot distinguish between baryonic matter, dark matter, and dark energy, and our definition of the density contains all three components. Hence, m includes the contributions of all matter and energy forms of the universe.
The variation as a function of the redshift of the function η is represented in Fig. 7. The coefficient of the one-form β of the Finsler metric is a monotonically increasing function of z (a monotonically decreasing function of time), and it is a linear function of the redshift. At small redshifts, the evolution of η is almost independent of the numerical values of w, but some small differences do appear at higher redshifts.

Discussion and final remarks
In this paper we have investigated the cosmological implications of a particular Finsler type geometry in which the fundamental function is obtained as the ratio of a Riemannian metric α and a one-form β, so that F = α 2 /β. The corresponding geometry is called Kropina geometry, and its properties have been intensively investigated from a mathematical point of view. We have assumed that the physical cosmological metric of the universe is the FLRW metric, from which the metric α 2 is constructed. From a geometric point of view, we have adopted the mathematical formalism of the osculating Finsler spaces, in which the Finsler metricĝ(x, y) is localized via the functional relation y = y(x), and thus a Riemann type metric g(x, y(x)) is generated. In our study we have considered Aosculating Riemannian manifoldsĝ(x, A), where A I are the components of the one-form β. Moreover, we assume that the Finsler space is an n-dimensional point Finsler space (locally Minkowskian but generally non-Euclidean). The connection of a point Finsler space is the Barthel connection, a connection that depends on the field to which it is applied. For an (α, β) metric, the Barthel connection is the Levi-Civita connection of the A -Riemannian metricĝ(x, A). From a physical point of view, we have postulated that the gravitational phenomena can be described by the A-Riemannian metric via the standard Einstein gravitational field equations, in which the Riemannian metric g(x) is substituted byĝ(x, A). Hence, the present approach is based on an extended (but still Riemannian) metric, with the connection and curvature constructed in the standard way.
By adopting as the α metric the standard FLRW form, the generalized Friedmann equations can be obtained by a number of straightforward calculations. Despite the apparent complexity of the mathematical model, the generalized Friedmann equations have a very simple mathematical form, and they are obtained in terms of two geometric quantities, the scale factor a and the component A 0 = aη of the one-form β. The generalized Friedmann equations have the remarkable property of giving, in the limit η → 1/a, the standard Friedmann equations of general relativity. Thus, the Barthel-Kropina-FLRW model represents a nontrivial deformation of the standard general relativistic cosmology. In the present work we have performed a systematic study of the cosmological properties of the Finslerian type model, which, despite its close relation to standard general relativistic cosmology, also has some very different properties. The generalized Friedmann equations of the Barthel-Kropina-FLRW model do not admit a vacuum solution for ρ = p = 0, since the first Friedmann equation gives η = 0, η = constant, with the second field equation identically satisfied. This situation has some similarities with standard general relativity, in which the vacuum Friedmann equations 3H 2 = 8π G/c 2 ρ and 2Ḣ +3H 2 = − 8π G/c 4 p give for ρ = p = 0 the solution H = 0, a = constant. Thus, the cosmological background solution of the vacuum gives the Minkowski geometry. On the other hand, in the Brathel-Kropina-FLRW model, the vacuum solution requires η = 0, η = constant, a condition that does not fix any background spacetime geometry, which remains of the FLRW type, but with arbitrary scale factor a. However, the model admits a de Sitter type solution for the case of a pressureless cosmological fluid in the presence of a nonzero matter density. The behavior of ρ in this case depends, as one can see from Eq. (73), on an arbitrary integration constant C 1 . If this constant is taken as C 1 = 1, lim x 0 →0 ρ = ∞, and therefore the de Sitter evolution begins from a singular state. Finite initial density states are also possible, and in the large time limit, both the matter density and η tend towards some constant values.
In the Barthel-Kropina-FLRW model, the matter energy density is generally not conserved, and the energy density and pressure satisfy the balance equation (65). Conservative cosmological models, obtained by imposing the condition of the energy conservation ρ + 3H ρc 2 + p = 0, lead to the standard general relativistic evolution for the scale factor in the presence of matter satisfying a linear barotropic equation of state, a = x 0 2/3γ . Thus, the cosmology of the conservative Barthel-Kropina-FLRW model is identical to the standard general relativistic one. Hence, accelerated expansion is not possible in the presence of matter satisfying the ordinary matter conservation equations. On the other hand, a large number of accelerating cosmological models can be obtained using the general energy balance equation (65), once the functional form of η is fixed. For the sake of the illustration of the qualitative properties of the model, we have considered two particular cosmological scenarios, with η(a) ∝ a −n and η = η(H ). The first choice of η leads to an eternally accelerating universe, with q = constant < 0, while the second choice leads to a plethora of cosmological models whose properties depend on the choice of the function η(h). We have considered and analyzed in detail a particular class of models obtained by choosing η(H ) = a 1 H n + b 1 . The cosmological parameters of this model have been compared with both the observational data and the predictions of the CDM model. Even without the use of a proper fitting procedure for the analysis of the observational data, the Barthel-Kropina-FLRW model can give at least an acceptable qualitative description of the observations of the Hubble function and reproduce the standard CDM model predictions for H . However, significant differences between the Barthel-Kropina-FLRW and the CDM model do appear in the evolution of the deceleration and matter density parameters. In order to test the validity of the present model, a detailed statistical analysis of its predictions with several classes of observational data is necessary.
An effective dark energy model can also be constructed within the framework of the Barthel-Kropina-FLRW geometric theory. Since in the limit η → 1/a one fully recovers the standard Friedmann equations, one can consider small deviations from the standard Friedmann cosmology by adding a correction term f to 1/a in the expression of η. This procedure allows the reformulation of the generalized Friedmann equations as the standard ones plus some new terms, which are essentially geometric. They can be interpreted as a geometric dark energy, and they contribute effective energy density ρ DE and pressure p DE terms to the standard baryonic thermodynamic quantities. By imposing a linear equation of state for the dark energy terms, p DE = wρ DE , the system of the cosmological equations can be closed and reformulated as a first-order dynamical system in the redshift space. The solutions of the system can be obtained numerically for different numerical values of the parameter w of the dark energy equation of state and for different initial conditions of the function f and its derivative f . We have performed a detailed comparison of the solutions of the system with the observational data, and the predictions of the CDM model. The Barthel-Kropina-FLRW model gives a good description of the observational data, and of the CDM model up to at least a redshift of z = 5. From the numerical simulations it follows that the numerical values of w have a relatively small influence on the cosmological behavior, which is strongly influenced by the present-day values of f and f .
Bouncing cosmological models represent an attractive alternative view of the early universe, and many such types of models have been proposed to explain the origin of the universe (see [126] for a review). Bouncing cosmologies have the attractive feature of providing a solution of the singularity problem that generally appears in the standard cosmological models. In the framework of the Finsler and Finslerlike geometries, the problem of the bouncing cosmologies was investigated in [127]. As a first general result, it was shown that in general very special relativity and in Finslerlike gravity on the tangent bundle, the bounce cannot be easily obtained. But in the Finsler-Randers geometry, by adopting for the scale factor the expression a(t) = a b 1 + Bt 2 1/3 , the scalar anisotropy induced by the geometry can satisfy the bounce conditions. Thus, bouncing solutions can be obtained in this class of gravitational theory. In theories constructed by using a nonlinear connection in which a scalar field does appear with an induced scalar-tensor structure, bouncing solutions can also be obtained. Interestingly enough, in the case of the non-holonomic basis, if one imposes to the quan- μ , a specific expression in terms of the scalar field and the Hubble function, then this N 0 can generate the bouncing scale factor.
Bouncing solutions can also be obtained, at least in principle, in the Barthel-Kropina cosmological model. By fixing in advance the form of the scale factor so that it has the required bouncing properties, the generalized Friedmann equations (52) and (53) can be reduced to a (strongly nonlinear) system of differential equations for η, in which the matter energy density and pressure terms are also present. The investigation of this system can be done only by using numerical (or perturbative) methods. However, in order to obtain bouncing solutions, a nonzero matter energy density must be present in the very early universe.
There are at least three possibilities for a possible theoretical understanding of the vast amount of cosmological data that have radically changed our view of the universe. The first approach is called the dark components model, and it generalizes the Einstein gravitational field equations by adding two new terms in the total cosmological energy momentum tensor which correspond to dark energy and dark matter, respectively. In this pathway, the gravitational dynamics are described by the generalized Einstein equation [128] where T bar μν , T DM μν (φ, ψ μ , . . .), and T DE μν (φ, ψ μ , . . .) represent the energy-momentum tensors of baryonic matter, dark matter, and dark energy, respectively. Hence, in this extension of general relativity, dark energy is a form of matter, an interpretation suggested by the equivalence principle of the theory of relativity that states the equivalence of mass and energy. The energy-momentum tensors of dark energy and dark matter can be realized by using scalar or vector fields. The simplest dark constituent model can be obtained by assuming that dark energy is a scalar field φ with a selfinteraction potential V (φ). For reviews and extensive discussions of the dark component models see [30][31][32][33][34][35]).
A second approach to the geometric description of the gravitational interaction is the dark gravity approach. The dark gravity formalism assumes a purely geometric description of the gravitational phenomena, and it is based on a modification of the geometric structure of the Einstein field equations, still formulated in a Riemannian geometry. In the dark gravity approach, one can formulate the Einstein gravitational field equations generally as where T (geom) μν g μν , R, R, . . . is a term generating a geometrically effective energy-momentum tensor, constructed from the Riemannian metric. The geometric term T (geom) μν g μν , R, R, . . . can describe dark matter, dark energy, or even both. An interesting example of dark gravity is the f (R) theory [113], in which the standard Hilbert-Einstein action S = R/κ 2 + L m √ −gd 4 x is generalized to an action is an arbitrary analytical function of the Ricci scalar R. For reviews of dark gravity theories and their applications, see [117,[129][130][131][132].
There is a third avenue for the understanding of the present-day cosmological data, called the dark coupling theory. In this approach, the standard Hilbert-Einstein Lagrangian density, having a simple additive structure in the geometric and matter Lagrangian, is replaced by a Lagrangian with a more general algebraic system. In the dark coupling theory, one looks for the maximal extension of the standard Hilbert-Einstein gravitational Lagrangian, and this extension can be obtained by assuming that the gravitational Lagrangian density is an arbitrary analytical function of the curvature scalar R describing the geometric properties of the spacetime and the matter Lagrangian L m . Other thermodynamic or geometric parameters can also be naturally included in specific models. The dark coupling approach leads natu-rally to the existence of a non-minimal coupling between curvature and matter.
In the dark coupling approach, the Einstein gravitational field equations can be formulated as with the effective energy-momentum tensor of the theory T (coup) μν g μν , R, L m , T, R, T, . . . obtained by considering the maximal extension of the Hilbert-Einstein Lagrangian. Moreover, a non-additive curvature-matter algebraic structure is introduced, describing the couplings between matter and spacetime geometry. The dark coupling theories were considered in [133], where a gravitational action of the form S= [ f 1 (R)+(1+λ f 2 (R))L m ] √ −gd 4 x was proposed. This action was generalized in [134] and [134], leading finally to the f (R, L m ) gravity theory [135], in which the gravitational Lagrangian density is given by an arbitrary function of the Ricci scalar and of the matter One can also couple curvature and matter via the trace of the matter energymomentum tensor, as in the f (R, T ) theory, with action given . For a detailed presentation of the theories with curvature-matter coupling see [114].
However, a fourth possibility for the description of the gravitational dynamics and evolution, including the accelerating expansion, is also possibly in the framework of the dark geometry approach, in which one assumes that the true geometry of the nature is beyond the Riemann one and that the extra terms generated by the non-Riemannian mathematical structures may be responsible for the existence of dark matter and dark energy. One possible dark geometry candidate is Weyl geometry, with conformally invariant gravitational models explaining the present and early dynamics of the universe [137]. In the present work, we have presented another example of a dark geometry, the Finsler type Barthel-Kropina-FLRW geometry, in which an effective dark energy can be generated from the mathematical structures underlying the geometry. The gravitational field equations are postulated as having a similar form as the Einstein equations in Riemann geometry, but with the curvature tensors replaced by their Finslerian counterparts. There is a close relation between the Riemannian cosmological evolution equations and the Finslerian ones, and this relation allows the natural introduction of a geometric dark energy term in the gravitational formalism.
Hence, the Barthel-Kropina-FLRW model not only may represent an attractive alternative to the standard CDM model, but could also open new avenues for the understanding of the complex relation between mathematics and physical reality. a grant of the Romanian Ministry of Education and Research, CNCS-UEFISCDI, project number PN-III-P4-ID-PCE-2020-2255 (PNCDI  III).

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: There are no data associated with this article because of its theoretical and formal nature.] 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 .

Appendix A: Computation of the Barthel-Kropina metric
1 Method 1 We recall the formula for the fundamental tensor of this Kropina metric in [138] g I J (x, y) = 2α 2 β 2 g I J (x) + where g 0I := g I J y J (see [138]). Note that the index zero in the g terms means contraction by y I . Using the above formula, we will compute the components of the metric tensor under the conditions specified below.
Step 2. We evaluate this matrix for arbitrary x and the specific direction Then the Kropina fundamental tensor reads 2 Method 2 Let α = g I J (x)y I y J and β = A I (x)y I , and The Hessianĝ where α I := α y I and ρ = φ 2 − sφφ , ρ 0 = φφ + φ φ , ρ 1 = −s(φφ + φ φ ) + φφ . We now consider y i = A i , then Therefore, Hence, the nonvanishing components ofĝ I J arê

Appendix C: Computation of the Ricci tensor
We start by computing the derivatives of the nonvanishing Christoffel symbols, Next, we compute the Ricci tensor. Recall the formula for the Ricci tensor, that is, Thus we obtain and

Appendix E: The generalized Friedmann equations
Recall the formula of the Einstein tensor, It follows that The rest of the components of the Einstein tensor vanish.