Cosmological Finsler Spacetimes

Applying the cosmological principle to Finsler spacetimes, we identify the Lie Algebra of symmetry generators of spatially homogeneous and isotropic Finsler geometries, thus generalising Friedmann-Lema\^{i}tre-Robertson-Walker geometry. In particular, we find the most general spatially homogeneous and isotropic Berwald spacetimes, which are Finsler spacetimes that can be regarded as closest to pseudo-Riemannian geometry. They are defined by a Finsler Lagrangian built from a zero-homogeneous function on the tangent bundle, which encodes the velocity dependence of the Finsler Lagrangian in a very specific way. The obtained cosmological Berwald geometries are candidates for the description of the geometry of the universe, when they are obtained as solutions from a Finsler gravity equation.


Introduction
To describe the evolution of the whole universe in cosmology, one applies the cosmological principle (CP), which states that there exists no preferred spatial position and no preferred spatial direction on large scales. Applying this principle to general relativity leads to the spatially homogeneous and isotropic Friedmann-Lemaître-Robertson-Walker (FLRW) metric as the unique ansatz for the geometry of spacetime. It contains two free functions which depend only on time, the lapse function and the scale factor. The lapse function can be normalized to unity, by a suitable choice of the time coordinate. The scale factor remains the only free function to be determined as solution of the Einstein equations, sourced by a perfect fluid energy-momentum tensor. On the basis of this mathematical model for the universe, one has to conclude that only ∼5% of the Universe consists of standard model baryonic matter, while the rest is composed of what is nowadays called dark energy [1,2] and dark matter [3]. The standard approach to cosmology is excellently summarized for example in reference [4].
A promising approach for a geometric explanation of the dark matter and dark energy phenomenology is to use Finsler spacetime geometry for the description of the gravitational interaction, instead of pseudo-Riemannian geometry [5][6][7][8][9]. In particular, it has recently been suggested that Finsler geometry provides the correct mathematical framework [10] and extension of the Einstein equations [11] for the accurate determination of the gravitational field distribution of a kinetic gas.
In this article we apply the cosmological principle to Finsler spacetime geometry. Starting from a symmetry group which acts transitively on spatial equal time surfaces and which contains a local isotropy group acting transitively on spatial directions at each point, we find that a cosmological homogeneous and isotropic Finsler geometry is defined by a Finsler Lagrangian with a very specific dependence on arXiv:2003.02299v2 [gr-qc] 5 May 2020 the tangent bundle coordinates. Yet, since the Finsler Lagrangian is a 2-homogeneous function in its dependence on the directional variable of the tangent bundle, the demand of cosmological symmetry leaves large classes of allowed Finsler Lagrangians; the symmetry demand does not provide a strong limitation in this regard. Specific choices of Finsler geometries have been investigated in their capability to explain aspects of the cosmological dark matter and dark energy phenomenology [12][13][14][15][16][17][18].
A class of Finsler spacetime geometries, which can be regarded as closest to pseudo-Riemannian geometry, are the so-called Berwald spacetimes [19][20][21][22]. They are characterized by the fact that, in any local chart, the Chern-Rund connection coefficients on the tangent bundle only depend on the points of the base manifold, or, equivalently, the geodesic spray coefficients are quadratic in their dependence on the directional variables of the tangent bundle. In other words, on Berwald spacetimes, the Chern-Rund connection gives rise to an affine connection on the spacetime manifold. In general, this connection is not the Levi-Civita connection of any pseudo-Riemannian metric, but instead, it can be regarded as a non metric-compatible metric-affine connection without torsion [23].
We derive the most general cosmologically (spatially homogeneous and isotropic) Berwald spacetime geometry. It serves as simplest Finslerian candidate for the description of the geometry of the universe. The obtained Finsler Lagrangian contains one free function, which encodes the velocity dependence of the Finsler Lagrangian in a very specific way. The Berwald geometry we obtain is the minimal Finsler geometric extension of pseudo-Riemannian FLRW geometry.
We present our results in the following way. In Section 2, we introduce the definition of Finsler spacetimes and the mathematical language needed to discuss the cosmological principle and cosmologically symmetric Finsler spacetimes. We apply the cosmological principle to Finsler spacetime geometry and identify the symmetry generating vector fields in Section 3. In Section 4 we recall how to identify Berwald spacetimes and derive the most general cosmologically symmetric Berwald spacetime Finsler Lagrangian. Finally, we conclude in Section 5.

Berwald Finsler Spacetime Geometry
Throughout this article, we consider the tangent bundle TM of a 4-dimensional connected and oriented manifold M, equipped with manifold induced local coordinates, as follows. A point (x,ẋ) ∈ TM will have local coordinates of the form (x a ,ẋ a ), where x a are the local coordinates of the point x ∈ M anḋ x =ẋ a ∂ a ∈ T x M is the decomposition of the vectorẋ ∈ T x M in the natural basis. If there is no risk of confusion, we will sometimes suppress the indices of the coordinates. The symbol π denotes the canonical projection of the tangent bundle. The local coordinate bases of the tangent and cotangent spaces, T (x,ẋ) TM and T * (x,ẋ) TM, of the tangent bundle are {∂ a = ∂ ∂x a ,∂ a = ∂ ∂ẋ a } and {dx a , dẋ a }. A conic subbundle of TM is a non-empty open submanifold Q ⊂ TM\{0}, with the following properties: conic property: if (x,ẋ) ∈ Q, then, for any λ > 0 : (x, λẋ) ∈ Q.
By a Finsler spacetime we will understand in the following a pair (M, L), where M is a smooth n-dimensional manifold and the Finsler Lagrangian L : A → R is a smooth function on a conic subbundle A ⊂ TM, such that: • L is positively homogeneous of degree two with respect toẋ: L(x, λẋ) = λ 2 L(x,ẋ) for all λ ∈ R + ; • on A, the vertical Hessian of L, called L-metric, is non-degenerate, • there exists a conic subset T ⊂ A such that on T , L > 0, g L has Lorentzian signature (+, −, −, −) and on the boundary ∂T , L can be continuously extended as L| ∂T = 0. 1 This is a refined version of the definition of Finsler spacetimes in Reference [11] and basically covers, if one chooses A = T , the improper Finsler spacetimes defined in Reference [24].
The 1-homogeneous function F, which defines the point particle action for curves γ on M , is obtained from the Finsler-Lagrange function L as F = |L| and interpreted as proper time integral of observers. For clarity, we list the different sets which appear in the above definition and comment on their meaning: We like to point again the relation ∂T ⊂ N is demanded in our definition. An important building block of the geometry of Finsler spacetimes is the geodesic spray, locally given by the coefficients It defines the Finsler geodesic equation in arclength parametrizationẍ a + 2G a (x,ẋ) = 0, the canonical (Cartan) nonlinear connection coefficients G a b =∂ b G a and the Berwald linear connection coefficients G a bc =∂ c∂b G a . A Finsler spacetime is called of Berwald type [19,25], or simply Berwald spacetime, if and only if, in any local chart, the geodesic spray is quadratic in its dependence on the tangent space coordinatesẋ: This is equivalent to demanding that the canonical nonlinear connection coefficients are actually linear in theirẋ dependence, or that the Berwald linear connection coefficients are independent ofẋ. The latter means that the G a bc (x) define an affine connection on M. Next, we will determine the most general cosmologically symmetric Finsler spacetimes from the CP, before we derive the most general homogeneous and isotropic Berwald spacetime in Section 4.

The Cosmological Principle on Finsler Spacetimes
We will prove in the following that, applying the CP to Finsler spacetime geometry singles out 6 symmetry generating vector fields, which characterize spatially homogeneous and isotropic Finsler spacetimes. More precisely, these symmetry generating vector fields are the same as the ones defining cosmological symmetry in general relativity, that is, the same as in the case of pseudo-Riemannian geometry. 1 It is possible to equivalently formulate this property with opposite sign of L and metric g L of signature (−, +, +, +). We fixed the signature and sign of L here to simplify the discussion.
The key argument is that spatial homogeneity and isotropy ensure the existence of a maximally symmetric Riemannian metric on the 3-dimensional time slice hypersurfaces of a Finsler spacetime. The Finsler Lagrangians we will determine will involve this maximally symmetric Riemannian metric, albeit, in a non-trivial way.
We begin by assuming that the spacetime in consideration possesses a smooth global time function t : M → R, which assigns to each p ∈ M a time stamp and whose differential dt satisfies dt(X) > 0, ∀X ∈ T ∪ ∂T . The level sets of the time function Σ T = {p ∈ M|t(p) = T = constant} are interpreted as equal-time spatial hypersurfaces. We assume that these level sets are connected.
The demand of the existence of a symmetry group acting on the hypersurfaces Σ T is the core of the cosmological principle.

The Symmetry and the Isotropy Group
A Finsler spacetime geometry satisfies the CP if it is (see for example Wald [26] or Weinberg [4] for pseudo-Riemannian spacetimes): • spatially homogeneous, that is, for each fixed value T of the time function t and for any two points q 1 and q 2 in Σ T , there exists a diffeomorphism of M which maps q 1 to q 2 and preserves the Finsler Lagrangian L. In other words, there exists a Lie Group G of isometries of (M, L) acting transitively on each slice Σ T ; • spatially isotropic, that is, around each p ∈ M there exists a congruence of observer curves γ(t, s) with tangent vector field ∂γ ∂t ∈ T , such that for each two lines [ In other words, the stabilizer G p = {ψ ∈ G|ψ(p) = p} at p ∈ M acts transitively on the projective tangent spaces PT p Σ T of Σ T .
Fix T ∈ R. We start by a first remark: if a Lie group acts transitively on a connected manifold, then the connected component of the identity element is a Lie subgroup that still acts transitively on the respective manifold [29] (p. 398).
Since the slice Σ T is connected, there is no loss of generality if we assume that G is connected (in the contrary case, we can restrict our attention to the connected component G 0 of the identity in G, which still acts transitively on Σ T . The corresponding isotropy group will then be G p ∩ G 0 ). The connected component of the identity in G only contains orientation-preserving diffeomorphisms ψ : M → M.
We note that the homogeneity demand makes (Σ T , G) a homogeneous space; thus, Σ T is diffeomorphic to the quotient G/G p and the following relation between the dimensions of the involved sets holds We will now prove two lemmas to identify the dimension of the groups G and G p .
On a Finsler spacetime (M, L) satisfying the CP, the dimension of the space of Killing vectors at each point p ∈ M is at most 6 and the dimension of the isotropy group G p is at most 3.

Proof.
Locally, in a coordinate chart around each p ∈ M, the action of G on M is determined by the generating vector fields ξ = ξ a ∂ a ; the 1-parameter subgroup of G generated by such a vector field ξ acts on M by some diffeomorphisms ψ locally expressed as: ψ a (x) = x a + ξ a (x) + O( 2 ). On a Finsler spacetime (M, L), these diffeomorphsims are isometries if and only if ξ C (L) = 0, where ξ C = ξ a ∂ a +ẋ b ∂ b ξ a∂ a is the complete lift of ξ to the tangent bundle.
Using the Killing equation ξ C (L) = 0 on a Finsler spacetime (M, L), it was proven in Reference [28] (Section 5.4.4) that any Killing vector field ξ = ξ a ∂ a in local coordinates around any point p ∈ M, is locally uniquely determined by the values ξ a (p) and by the derivatives ∂ b ξ a (p), with a < b. For an n-dimensional Finsler spacetime this means, we can freely choose at most n values ξ a (p) and at most n(n−1) 2 values ∂ b ξ a (p); all in all, this gives at most n + n(n−1) 2 = n(n+1) 2 independent Killing vectors. Thus, for n = 4 the dimension of the isometry group can be at most 10.
Choose an arbitrary local chart such that the first coordinate is the time function x 0 = t. Since the diffeomorphisms ψ ∈ G preserve the time slices Σ T , we find that ξ(T) = 0, that is, ξ 0 and its derivatives identically vanish and that the remaining components ξ i , i = 1, 2, 3 can only depend on x i , i = 1, 2, 3. That is, we can only pick freely, at a given point, 3 values ξ a (p) and 3 derivatives ∂ b ξ a (p), with a < b, which gives us at most 6 degrees of freedom for the space of Killing vectors at each point. In particular, since the isotropy subgroup G p fixes the point p, we must have ξ a (p) = 0 for all a = 0, . . . , 3, hence we can only freely choose the 3 derivatives ∂ b ξ a (p) with a < b. Consequently, the dimension of G p cannot exceed 3, which proves our first remark.

Lemma 2.
The dimension of the isotropy group G p is at least 3.
Proof. Let us properly understand the statement "the isotropy group G p = {ψ ∈ G|ψ(p) = p} at p ∈ M acts transitively on the projective tangent spaces PT p Σ T of Σ T ". We will identify a transitive and effective action of a quotient group of G p on the projective tangent space PT p Σ T . This group is obtained by identifying elements of G p that have the same linear tangent map at p.
Any Lie group action on a manifold gives rise to an effective Lie group action on the respective manifold, by factorizing away the elements that provide trivial actions. In our case, assume the group G p does not act effectively on the (for the moment, non-projectivized) tangent space T p Σ T and denote by Id p the subset of G p whose elements provide trivial actions. Then, Id p is a normal subgroup of G p and the factor group acts effectively on T p Σ T , by the rule (dψ • Id p )v := dψ (v). The group G p can be identified with a subgroup of the general linear group GL 3 (R) (more precisely, of the connected component of GL 3 (R) consisting with matrices with positive determinants). This is justified as follows. Fix an arbitrary coordinate chart on TM. The action the subgroup G p fixing p on the tangent space T p Σ T is then expressed as a matrix multiplication, namely, by the Jacobian matrices of the diffeomorphisms ψ ∈ G p . Factorizing G p by the subgroup Id p actually means identifying as a single element those diffeomorphisms ψ, ψ of M which have the same values at p and the same Jacobian matrices (but whose higher order derivatives at p might differ). This way, the mapping from G p to GL 3 (R) associating to a class [ψ] ∈ G p the Jacobian matrix at p of ψ, is an injective homomorphism, which gives us the right of identifying G p as a subgroup of GL 3 (R).
Further, passing to the projectivised tangent spaces (that can be identified with the projective plane PR 3 once a choice of the basis of T p Σ T is made), the group that naturally acts is then the subgroup PG p of the projective group PGL 3 (R) obtained by factorizing GL 3 (R) by the group of rescalings A → λA, with λ ∈ R. Hence, on one hand, the isotropy request that G p acts transitively on the projective tangent space PT p Σ T implies that the group PG p (which is isomorphic to a subgroup of the projective linear group PGL 3 (R)), acts transitively and effectively on PT p Σ T . The latter is true since PGL 3 (R) acts effectively on PT p Σ T and so do all of its subgroups, hence in particular PG p . On the other hand, we notice that the only rescalings that can belong to G p are those with positive factors λ. But these cannot be L-isometries unless λ = 1, since L(λV) = λ 2 v for all v ∈ T p Σ T . That is, the subgroup of G p we have used for factorization is actually trivial, that is, There is actually a smaller subgroup of PGL 3 (R) which still acts transitively on PT p Σ T . To identify this subgroup, we mention the following result ( [27], [29], p. 398): If a connected Lie group acts on a compact manifold with finite fundamental group, then any maximal compact subgroup also acts transitively on the respective manifold.
The projectivized space PT p Σ T is connected, compact and its fundamental group is Z 2 . Since PT p Σ T is connected, the connected component of the identity in G p still acts transitively and effectively on it.
Denoting by H p a maximal compact subgroup of the connected component of the identity in PG p , we find that H p also acts transitively and effectively on PT p Σ T .
Finally, we use a result in Reference [29] (pp. 398-401), stating that any connected, compact Lie group acting transitively and effectively on the projective plane PR 3 is isomorphic to SO(3); that is, our subgroup H p of the isotropy group PG p must be isomorphic to SO(3) and thus: From Lemma 1 and Lemma 2 we find Theorem 1. On a Finsler spacetime satisfying the cosmological principle, the dimension dim G p of the isotropy group G p is 3 and the dimension dim G of the full symmetry group G is 6.

Remark 1.
The above result intuitively tells us that, on one hand, by taking into account only the linear approximation at p of the diffeomorphisms ψ and, on the other hand, by taking into account just a bounded region around the identity of the symmetry group, we do not lose any generators of the Lie algebra.
We can actually state a much stronger result.

Remark 2.
On Finsler spacetimes with cosmological symmetry, the isotropy group G p is compact.
Proof. First we show that G p is connected, and then we apply Cartan's Theorem on connected Lie groups.
• Under the above assumption that G is connected, the isotropy group G p is also connected. To see this, let us assume that G p is the disjoint union of at least two connected components, say, G p = A B, where A is the component of the identity. Since PT p Σ T is connected and G p acts transitively on it, it follows that A also acts transitively on PT p Σ T and see Reference [29] (p. 395), where H v is the stabilizer of a point , ∀v ∈ T p Σ T } consists of rescalings with positive factors; recalling that rescalings with positive factors cannot be L-isometries, it follows that H v is trivial, hence, G p coincides with its identity connected component A.

•
Further, we apply Cartan's classification theorem ( [29], p. 389 ), to G p . The Theorem states that any connected Lie group is the direct product between one of its maximal compact subgroups, say H, and a Euclidean space. Since G p acts transitively on the compact manifold PT p Σ T (which, as we have seen above, has finite fundamental group), its maximally compact subgroup H also acts transitively on PT p Σ T . But, the smallest compact group that can act transitively on PT p Σ T is 3-dimensional (more precisely, SO(3)), that is, dim H ≥ 3. Taking into account that the dimension of G p itself is 3, it means that G p = H, that is G p itself is compact.

The Symmetry Generators
To explicitly determine the generators of the groups G and G p , we use the above remark, which states that G p is compact. Then, we take into account that Σ T is a homogeneous manifold having a compact isotropy group and thus it must admit a G-invariant Riemannian metric h ([30], Example 1.3, p. 154). In other words, the generators of our group G are also Killing vector fields of h. But, for a 3 dimensional Riemannian manifold there can exist at most 6 Killing vector fields, hence Σ T is maximally symmetric and h is the corresponding maximally symmetric metric. In particular, h has constant scalar curvature k. The Riemannian metric h also gives us the possibility of having local spherical coordinates (t, r, θ, φ) around p. In these coordinates, see Reference [4,31,32], the isotropy group generators (which are the elements of the Lie algebra so (3)) are written as: and the generators of quasi translations, as: Finally, solving the Finsler Killing equation X C I (L) =, I = 1, ..., 6, where X C = X a ∂ a +ẋ b ∂ b X a∂ a yields that L(t, r, θ, φ,ṫ,ṙ,θ,φ) = L(t,ṫ, w), w 2 =ṙ is the most general spatially homogeneous and isotropic Finsler Lagrangian. The explicit calculation, as well as the expressions for the complete lifts X C I can be found for example in Reference [31].

Homogeneous and Isotropic Berwald Spacetimes
In order to find the desired Berwald Finsler spacetimes, we rewrite a generic Finsler Lagrangian in a specific way, which allows us to reduce the condition that a Finsler spacetime shall be Berwald, to a first order partial differential equation.

The Berwald Condition
Every Finsler spacetime Lagragian L can be written as L(x,ẋ) = g(ẋ,ẋ)Ω(x,ẋ), where g is an arbitrary pseudo-Riemannian metric, g(ẋ,ẋ) = g ab (x)ẋ aẋb and Ω = Ω(x,ẋ) is a 0-homogeneous function inẋ. In Reference [33], it was proven that L defines a Berwald Finsler geometry if and only if there exist a g and a (1, 2)-tensor field D on M, symmetric in its vector arguments, such that Ω satisfies the equation which we call the Berwald condition. Here, the indices were raised and lowered with the pseudo-Riemannian metric, that is,ẋ b =ẋ a g ab (x). The connection coefficients G a bc in (4) are then determined by the Christoffel symbols Γ of g and the tensor D as In case the expansion of L = g(ẋ,ẋ)Ω(x,ẋ) is done with a proper pseudo-Riemannian metric, where g(ẋ,ẋ) has a non trivial null structure, one might worry what happens on this null structure.
To avoid this problem observe that, if one has found one pseudo-Riemannian metric g, a tensor D and a factor Ω which solve (13), the Berwald condition is satisfied for any alternative expansion of L of this type, that is, for any other pseudo-Riemannian metricg, we can find a corresponding (1,2)-tensor field D and a factorΩ.
To see this, let us expand L = g(ẋ,ẋ)Ω(x,ẋ) =g(ẋ,ẋ)Ω(x,ẋ), where g and Ω satisfy (13) for a tensor field D. Then δ a L = ∂ a − G b a∂b L = 0 implies Using that G a b = (Γ a bc + D a bc )ẋ c and introducing the tensor field componentsD a bc = D a bc + Γ a bc − Γ a bc we find Hence, alsog,Ω andD satisfy the Berwald condition. Thus the metric factor in the expansion of the Finsler Lagrangian can be chosen arbitrarily; in particular, we can choose a positive definite metric g in order to avoid complications on the null structure of the metric.
In the following, we will insert the most general g and D which are compatible with a specific spacetime symmetry and solve the Berwald condition for Ω. This determines the most general Berwald Finsler spacetime for the desired spacetime symmetry.

The Cosmological Berwald Condition
To evaluate the Berwald condition, we consider Finsler Lagrangians L which are the form L(t,ṫ, w) = (ṫ 2 + σa(t) 2 w 2 )Ω(t,ṫ, w), where the metric factor g(ẋ,ẋ) = (ṫ 2 + σa(t) 2 w 2 ) is the most general homogeneous and isotropic (pseudo)-Riemannian metric with positive definite (σ = 1) resp. Lorentzian (σ = −1) signature. As we discussed, the Berwald condition (13) does not depend on this choice, which will also explicitly follow from the upcoming calculation. The possibility of choosing σ = 1 means that the set of null vectors of the chosen metric does not interfere with our result. On the conic bundle T , wherė t = 0, this expression can be nicely rewritten in terms of the 0-homogeneous variable s = w/ṫ L(t,ṫ, w) =ṫ 2 (1 + σa(t) 2 s 2 )Ω(t, 1, s) , which will be very convenient to evaluate the Berwald condition.

of 14
The second ingredient in this condition is the (1, 2)-tensor field D. The most general spatially homogeneous and isotropic such tensor field that is symmetric in its vector arguments, has the following nonzero components, see for example Reference [34], where b(t), c(t), d(t) are arbitrary functions of t.
Using the above expressions in the Berwald condition (13) yields two independent equations, which need to be solved to determine Ω(t, s) ≡ Ω(t, 1, s). Let denote derivatives with respect to the single argument of the functions a, b, c and d, then the spatial equations, if the index a assumes the values r, θ or ∂s Ω(t, s) = 0 (20) and the temporal equation, for which the index a assumes the value t is

Solving the Cosmological Berwald Condition
Solving the Berwald condition for the most general combination of Ω(t, s), b(t), c(t) and d(t) is not trivial. We will now classify all the solutions of the system (20) and (21).
Introducing the functions the Berwald condition becomes We now analyze the first equation and find several cases in which we obtain trivial solutions, in the sense that the Finsler Lagrangian is pseudo-Riemannian or zero. The only case that provides proper Finslerian solutions is M = N = 0, as we will see below.

Trivial Solutions
Trivial solutions arise in the following situations:

•
If N is different from zero, then we can divide the first equation in (24) by N, Now it is helpful to introduce the function (26) and to realize that it satisfies We can thus rewrite Equation (25) as After multiplication with A(t, s) we find where f (t) is an arbitrary function of t. Using the explicit form of the function A, we thus obtain the solution Constructing the Finsler Lagrangian L =ṫ 2 (1 − a(t) 2 s 2 )Ω(t, s) from this solution immediately yields with I(t) = f (t)(a (t) + a(t)d(t)) and J(t) = a(t) f (t)(σa(t)a (t) − c(t)). The just constructed Finsler Lagrangian is again quadratic in its dependence on the velocities and hence defines a pseudo-Riemannian spacetime geometry. In the particular case when M(t, s) is zero, we see from (25) that Ω is just a function of t.

•
If N = 0 and M = 0, then the first equation in (24) implies immediately that Ω(t, s) = 0 and thus the Finsler Lagrangian is L = 0.
From this analysis, we find that nontrivial cosmologically symmetric Berwald Finsler Lagrangians can only be obtained if M = N = 0.

Finslerian Solutions
Demanding that M = N = 0 leads to the equations for c(t) and d(t). Since s and t are independent variables, the s dependence in the latter must vanish, which immediately implies c(t) = σa(t)a (t). Plugging this into the remaining equations yields d(t) = − a (t) a(t) . Having solved (20) for general Ω(t, s), the remaining Equation (21) becomes This equation intertwines the t and s dependence of Ω and can be solved by a change of variables. For this purpose we substitute s by a new variable u, which is defined such that In other words, the function B(t) is given explicitly as the integral up to an undetermined constant of integration related to the choice of the lower bound t 0 . Replacing we find that Using our original Equation (33) becomes and so it contains derivatives with respect to t only. We can explicitly integrate this equation by introducing another function hence simplifying Equation (39) to After multiplication by C(t, u) we thus conclude with an arbitrary free function f (u) which depends only on u. Substituting back we obtain the general Finslerian solution of the Berwald condition Recall from the definition (35) that B(t) is determined only up to a multiplicative constant; however, this constant can simply be absorbed into the function f . With this result, we found the most general nontrivial Berwald Finsler spacetimes with cosmological symmetry: The only remnant of the two different possible choices for the metric in the decomposition of the Finsler Lagrangian appears in the overall factor σ, which can, without loss of generality, be absorbed into the free function f .
One further step of simplification can be done by introducing the new coordinatet(t) = t 0 dλB(λ), which implies dt dt = B(t) and thus, for the tangent bundle coordinates,ṫ =ṫB(t) ands = w/ṫ. In these new coordinates (45) becomes eventually there is only one free function f , which needs to be determined by Finsler gravitational field equations. The null structure of the above Finsler Lagrangian is determined by the zeros of f , and surely f has to be choosen in such a way that L defines a Finsler spacetime according its definition in Section 2.
To summarize our findings we have proven the following theorem: Theorem 2. If a Finsler spacetime Lagrangian L is of Berwald type and admits cosmological, spatial homogeneous and isotropic symmetry, then it falls into one of the following classes: 1. pseudo-Riemannian (quadratic inẋ), in which case it is, up to t coordinate redefinition, given by the FLRW metric, or 2. nontrivially Finslerian, in which case it is of the form (45).
As an explicit example, one may consider f = (c −s 2 )h(s) with c being a constant and h being a smooth and non-vanishing function. The resulting Finsler Lagrangian then is which satisfies all Finsler spacetime criteria. Still, the function h can be chosen freely and must be determined from the gravitational field equation. The latter is work in progress.

Discussion and Conclusions
The cosmological principle assumes the existence of a symmetry group G that acts transitively on spatial hypersurfaces of spacetime, and which contains a local isotropy group G p ⊂ G. In the context of Finsler geometry, so far, it has been assumed that the generators of this symmetry group have the form (10) and (11). So far this had not been derived from first principles. We closed this gap by showing that the dimension of G must be 6 and that the dimension of G p must be 3. Thus the dimension of these groups is the same in the pseudo-Finsler and in the pseudo-Riemannian setting. Moreover, with the help of an auxiliary metric on the spatial slices sharing the same symmetry group, which is guaranteed to exist, we could conclude that the symmetry generators indeed must have the assumed form.
Among the variety of possible Finsler geometric extensions of pseudo-Riemannian geometry as geometry of spacetime, Berwald spacetimes represent a most conservative generalization. Our discovery of the most general non-trivial cosmological, that is, spatially homogeneous and isotropic Berwald spacetimes reveals the class of geometries which extend the famous FLRW class of metrics into this realm. Most importantly, we found that cosmological Berwald geometries are parametrized by a free 0-homogeneous function on the tangent bundle, which intertwines the position and direction dependence of the Finsler Lagrangian in a very specific way. The resulting Finsler Lagrangian is L =ṫ 2 f (s) .
As the scale factor is determined by the Einstein equations on general relativity, the free function must be determined by suitable Finsler generalisations of the Einstein equations. Most of the suggested generalizations in the literature simplify significantly for Berwald geometries.
In particular, the ansatz (47) is an important step in the program of the description of the evolution of the universe in terms of a gravitational field distribution sourced by a kinetic gas. We argued in Reference [10] that the back reaction of a kinetic gas on the geometry of spacetime can be obtained directly from the 1-particle distribution function (1PDF) of the gas, when one employs Finsler geometry instead of pseudo-Riemannian geometry. The explicit form of the 1PDF will then determine the free function f , a derivation which is currently work in progress.

Author Contributions:
The authors have all contributed substantially to the derivation of the presented results as well as analysis, drafting, review, and finalization of the manuscript. All authors have read and agreed to the published version of the manuscript.