Schwarzschild-like solutions in Finsler–Randers gravity

In this work, we extend for the first time the spherically symmetric Schwarzschild and Schwarzschild–De Sitter solutions with a Finsler–Randers-type perturbation which is generated by a covector Aγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_\gamma $$\end{document}. This gives a locally anisotropic character to the metric and induces a deviation from the Riemannian models of gravity. A natural framework for this study is the Lorentz tangent bundle of a spacetime manifold. We apply the generalized field equations to the perturbed metric and derive the dynamics for the covector Aγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_\gamma $$\end{document}. Finally, we find the timelike, spacelike and null paths on the Schwarzschild–Randers spacetime, we solve the timelike ones numerically and we compare them with the classic geodesics of general relativity. The obtained solutions are new and they enrich the corresponding literature.


Introduction
In the context of research on generalized metric spaces, Finsler, Lagrange and Finsler-like geometries have played an important role in modified theories of gravity, general relativity and cosmology. In the last two decades, a development of these promising topics of research has extended the limits of general relativity and cosmology by including locally-anisotropic approaches. Over the last two decades field equations have been thoroughly studied in the context of Finsler, Lagrange, generalized Finsler and Finslerlike geometries as well as in scalar-tensor theories. The basic feature of these theories is the presence of extra terms in the equations of motion due to the intrinsic geometrical spacetime anisotropy. In this framework the term "spacetime a e-mail: alktrian@phys.uoa.gr (corresponding author) b e-mail: svasil@academyofathens.gr c e-mail: manoliskapsabelis@yahoo.gr d e-mail: pstavrin@math.uoa.gr anisotropy" is related to the Lorentz violation feature of the geometry. This approach may provide the necessary platform in understanding one of the most crucial problems in cosmology which is related with the underlying mechanism of cosmic acceleration and thus of dark energy. There is a lot of work in the literature on "spacetime anisotropic" geometries and below we briefly present some relevant works.
Generalized Einstein field equations have been studied in the Finsler, Lagrange, generalized Finsler and Finsler-like spaces, for an osculating gravitational approach in which the second variable y(x) is a tangent/vector field [1][2][3] and in Finsler cosmology [4][5][6][7][8]. Different sets of generalized Einstein field equations were derived for the aforementioned spaces in the framework of a tangent bundle [9][10][11][12][13][14][15][16][17] and for the momentum space on the cotangent bundle [18][19][20][21][22][23]. Additionally, Lorentz invariance violation in Finsler/Finslerlike spacetime and in Finsler cosmology in very special relativity has also been studied in a large series of papers [5,[24][25][26][27]. Investigations on generalized scalar-tensor theories with Finsler-like structure modeled on a vector bundle with two-internal fibers have also been done [28][29][30]. Also, the causality problem and light cones for different types of Finsler spaces have been investigated in [31][32][33][34]. In this context the Raychaudhuri equations with locally anisotropic internal variables which are providing extra effective terms have been derived in [16,18,[35][36][37]. Investigations of the extended Friedmann equations in Finsler spaces with extra internal degrees of freedom [4] and dynamical analysis (critical points) [38,39] provide a better understanding of the dynamical properties of the Finsler-Randers cosmological models. To this end, articles in the framework of the weak field and pp-waves in Finsler spacetime can be found in [16,35,40] and potentially they can be used in order to test the performance of the Finslerian gravitational theory against current observations of gravitational waves.
It is well known in gravitation and cosmology that the Schwarzschild metric constitutes a fundamental ingredient of general relativity. This metric describes the most general spherical symmetric solution of the Einstein field equations in a region of spacetime where the energy-momentum tensor vanishes.
In the context of a Finsler/generalized Finsler space, an extension of a locally anisotropic perturbation of a Finsler type Schwarzschild metric has been studied by different authors. In some of these works, possible observational predictions are given based on the direction-dependent structure of spacetime. We would like to point out that our study is realized in a Schwarzschild-Randers model which is different from the works of other researchers on the Finslerian extensions of a classic Schwarzschild metric [41][42][43][44][45][46][47][48][49][50][51]. In our approach, we use sufficiently generalized Einstein field equations on a Lorentz tangent bundle of a spacetime manifold.
In Sect. 2 we introduce the basic framework and geometric structures of our model.
In Sect. 3 we study a Schwarzschild metric in a special Finsler-like spacetime of Randers type. This study provides a locally anisotropic perturbation of the classical Schwarzschild metric of the Riemannian structure in a natural way, induced by a covector field of the base manifold. In addition, the geometrical setting that we use, namely the framework of a Lorentz tangent bundle of a Riemannian spacetime, contains additional degrees of freedom compared to classic gravity. The generalized field equations which have been derived in [17] are applied on the perturbed metric of our Schwarzschild-Randers spacetime and are solved for the covector field.
In Sect. 4, we study particle paths for our generalized spacetime. We follow the approach in [17] which takes into account the effect of internal degrees of freedom on the point particle dynamics. We apply the solution of the covector derived in the previous sections and we obtain an explicit form for the path equations which is an extension of classical geodesics of general relativity. Finally, as an application, we solve the timelike paths numerically and compare them with the geodesics of general relativity.

Preliminaries
The natural background space for a locally anisotropic gravity is the tangent bundle of a differentiable Lorentzian spacetime manifold called a Lorentz Tangent Bundle (we will refer to it as T M hereafter) [15,22]. T M is itself an 8-dimensional differentiable manifold, so we can define coordinate charts and tensors on it in the usual way. We briefly present the basics for this structure, for more details see Sect. 2.
The Lorentz tangent bundle T M is locally covered by a coordinate map {x μ , y α } where the range of values for the indices of the x variables is κ, λ, μ, ν, . . . = 0, . . . , 3 and the range of values for the indices of the y variables is α, β, . . . , θ = 4, . . . , 7. An adapted basis on T M is and its dual basis is dx μ , δy α = dy α + N α ν dx ν . The bundle T M is equipped with a Sasaki-type metric G: where the metric of the horizontal space (h-space) g μν and the metric of the vertical space (v-space) v αβ are defined to be of Lorentzian signature (−, +, +, +). In the rest of this work, the following homogeneity conditions will be assumed: g μν (x, ky) = g μν (x, y), v αβ (x, ky) = v αβ (x, y), k > 0. These conditions are met when the following relations hold: where the functions F g , F v satisfy the following conditions: 3. The form defines a non-degenerate matrix: with g αβ =δ μ αδ ν β g μν and the sign in the rhs of (2), (3) is chosen so that the resulting metric has the corrrect signature.
When the above conditions are met, our metric is called a (pseudo-)Finsler metric. Details about the connection, curvature and torsion structures can be found in Sect. 2.

Field equations
In this chapter we consider a general metric (1) on the tangent bundle of a spacetime manifold constituted from a horizontal and a vertical part. In this consideration, we set the horizontal part to be the classical Schwarzschild metric and the vertical part to be a Randers-type perturbation of the Schwarzschild metric, so the form of our metric will be a Schwarzschild-Randers metric. We also consider a similar case for a Schwarzschild-De Sitter-Randers metric. This is a different type of metric than the Finsler-Randers one which has been considered in the framework of cosmological study by different authors [4,8,37,39,[52][53][54][55][56]. In these cases, the horizontal part is a Friedmann-Robertson-Walker model. In both models the velocity and the vertical part play the role of intrinsic anisotropic perturbation of the traditional metrics. We already know the form of the Schwarzschild metric, so we need to find an explicit form for the Randers-type perturbation which contains the additional information of local anisotropy.

Field equations on the Lorentz tangent bundle
In this paragraph, we present a set of field equations for the dynamic variables of our generalized framework. These equations are: where L M is the Lagrangian of the matter fields, δ μ ν and δ α β are the Kronecker symbols, |G| is the absolute value of the determinant of the total metric (1), and (13) are torsion components, where L α βν is defined in (A.14). From the form of (1) it follows that √ |G| = √ −g √ −v, with g, v the determinants of the metrics g μν , v αβ respectively. This relation was used in (10)- (12). Equations (7)-(9) are derived from an extension of the Hilbert-Einstein action on the eight-dimensional Lorentz tangent bundle and constitute a generalization of the Einstein field equations of general relativity. They are appropriate for the study of locally anisotropic models of gravity with internal degrees of freedom. For details on their derivation from a Hilbert-like action on the Lorentz tangent bundle see [17].
We will make some comments in order to give a physical interpretation in relation to the Eqs. (10)- (12). Lorentz violations produce anisotropies in the space and the matter sector. These act as a source of local anisotropy and can contribute to the torsion, connection, curvature components and to the energy-momentum tensors of the horizontal and vertical space T μν and Y αβ . As a result, the curvatures R μν and R as well as the energy-momentum tensor T μν contain the additional information of local anisotropy of the metric and the matter fields. S αβ , S and Y αβ , on the other hand, are objects with no equivalent in Riemannian gravity. They contain more information about local anisotropy which is produced from the metric v αβ which includes additional internal structure of spacetime. Finally, the nonlinear connection N α μ , a structure which induces an interaction between internal and external spaces [29], can also contribute to all the above-mentioned objects, while the energy-momentum tensor Z κ α shows the variation of L M with respect to N α μ and it reflects the dependence of matter fields on the nonlinear connection [17]. This is different from T μν and Y αβ which depend on just the external or internal structure respectively.
Notice that the field equations (7)-(9) reduce to the usual Einstein field equations of general relativity (GR) in the limit: We observe that when the metric g μν (x, y) is reduced to a Riemannian one according to (14) then the connection coefficients L μ κλ (x, y) (eq. (A.13)) reduce to the Christoffel symbols: and the Cartan torsion tensor (A.15) vanishes: From (15) it becomes clear that the internal structure (y variable) of spacetime doesn't provide additional information about the geometry, since the metric v αβ of the internal structure reduces to the metric g μν (x) of base spacetime (x variable). Moreover, we see that the Cartan-type [22] nonlinear connection (16) ensures that T α να = 0. Indeed, from Eq. (13) and from the conditions (14-16) we get Now, as long as conditions (14) and (15) Based on (20), it is straightforward to calculate the curvatures R μν and R from (A.29) and (A.30) respectively in the GR limit of the Lorentz tangent bundle gravity and we find that they reduce to the Ricci tensor R μν and Ricci scalar R of general relativity for the metric g μν . As expected, the curvatures S αβ and S (Eqs. (A.26) and (A.28)) both vanish in this limit. From the above relations (14)(15)(16)(17)(18)(19)(20) it follows that, in the GR limit, the corresponding field equations (7-9) boil down to where R μν and R are the Ricci tensor and scalar of general relativity for the metric g μν , as we mentioned above. Of course in Eq. (21) G denotes the Newton's constant and κ = 8π G c 4 , while the energy momentum tensor T μν is given by (10). From the the trace of Eq. (21) we obtain R = −κ T , hence Eq. (22) gives This is the GR limit for the energy-momentum tensor Y αβ . Finally, from (12) and (23) we conclude that in the GR limit the matter fields have no direct dependence on the nonlinear connection.

Schwarzschild-Randers spacetime
As we mentioned above, the horizontal part g μν of the metric (1) will be taken to be the Schwarzschild metric so that where R s = 2G M is the Schwarzschild radius (we have set the speed of light constant c = 1).
In the following, we assume a function F v of Randers type from which we will derive v αβ by using (3): (26) where g αβ = g μνδ μ αδ ν β is the Schwarzschild metric and A γ (x) is a covector which will be determined by the equations. We focus on the timelike subspace of the internal yspace with respect to the Schwarzschild metric (g αβ (x)y α y β < 0) hence the minus sign under the square root. We take A γ (x) to be a weak term (|A γ (x)| 1), hence we neglect high order terms from the calculations. In addition, we consider as the appropriate form for the nonlinear connection the one given in (16) i.e. its GR limit value on the tangent bundle: This choice will give us a locally anisotropic gravitational model which deviates minimally from general relativity due to the extra Randers term A γ y γ in (26).
We calculate the metric tensor v αβ of (26) from (3): and we find where we have set a = −g αβ y α y β . From (29) we see that the metric v αβ takes the form where we have set Its inverse is v βγ = g βγ − w βγ so that v αβ v βγ = g αβ g βγ = δ γ α to first order in w αβ . The total metric over the tangent bundle is then written as We remark that, as we can see from (30), the metric v αβ (x, y) is a Finslerian perturbation of the Riemannian metric g αβ (x).
We observe that if we let w αβ → 0 then the field equations (7)-(9) reduce to the Einstein field equations of general relativity.
Next we will calculate the terms for (7) and (8). From the definitions (A.29) and (A.30) we see that when g μν has no explicit dependence on y then R μν and R reduce to the classical Ricci tensor and scalar of general relativity. Additionally, since g μν (x) is the Schwarzschild metric, both R μν and R are zero. Here we have assumed vacuum solutions, so the energy momentum tensors are zero and the Eqs. (7) and (8) become Field equation (9) gives us no additional information since all three terms vanish identically in our case. We can simplify Eq. (34) by calculating C μ μδ from (A.15) and we find that it is zero since the metric g μν depends only on x. Then by taking the trace of the remaining terms in (34) we can show that S αβ and S are also zero, so the field equation (33) becomes We substitute (A.14) and (27) in (13) and after some calculations we get with w = g αβ w αβ . The above relation (36) shows us that the torsion is of first order on w αβ so the terms T γ κγ T β λβ from (35) are omitted. Then by taking the trace of the remaining terms in (35) we have the equation that follows: Substituting the latter equation to (35) we find By the definition of the covariant derivative in (A.10), Eq. becomes Using Eqs. (31) and (36) we get T α να in terms of A β : It is straightforward to show that with Using (40) and (41) we get from relation (43) we can calculate (39): Relation (45) must hold for every y. Since the expression in parentheses does not depend on y, we conclude that it must identically vanish: Remark. By comparing (A.14) and (42) we get a relation of the form L α βμ = E α βμ + O(A). If we use this in (44) and (46) we get the equation D (μ D ν) A γ = 0, which is equivalent to the system of Eqs. (44) and (46). Contracting this with g μν gives In order to fully determine the metric (32) for the respective subspace of the internal space, we need g μν (x) and A γ (x). The first is already defined in (25), so we need to solve (46) for A γ (x) to get a full expression for the metric in our space. If we use the definitions (42) and (44) on relation (46), we get the equation that we need to solve for A(x): where g μν (x) is the Schwarzschild metric (25). Once we get A(x) from (48), we can calculate w αβ (x, y) from (31) and then use the result to calculate the full metric (32). A similar analysis holds for the spatial subspace of the internal space. In that case, instead of (26) and (28) we have and v αβ = 1 2 where a = g μν y μ y ν for the spacelike sector of g μν taking into consideration the signature of the metric. Following the same steps as above, we reach the same equation for A γ , i.e. Eq. (48). Therefore, solving this equation will give us the metric for both the timelike and spacelike (with respect to g αβ ) sub-spaces of the internal space. We will solve (48) analytically with separation of variables, see Appendix B for more details. After calculations we find the solution withÃ 4 a constant. This is a timelike covector since g αβ A α A β = −(Ã 4 ) 2 < 0. It is interesting to mention that the horizon of the Schwarzschild-Randers metric is correlated with that of Schwarzschild. Practically, the quantity A γ can be seen as a distortion factor which quantifies the deviation from the pure Schwarzschild solution. Obviously, the solution (51) on small spherical scales (r ∼ R s ) tends to zero. On the other hand, for r R s the Schwarzschild -Randers metric tends asymptotically to Minkowski. Finally, we see that this metric has a singularity in r = 0 similarly with the classic Schwarzschild one. The Schwarzschild-Randers model will be further studied for intrinsic singularities and horizons in a future research.

Schwarzschild-De Sitter-Randers spacetime
We will follow the same procedure as in the previous paragraph but for a Schwarzschild-Randers spacetime with a cosmological horizon, namely a Schwarzschild-De Sitter-Randers spacetime. In this scenario, we take the horizontal part of the metric (1) to be: while, as before, the metric tensor v αβ will be derived from the Lagrangian (26) and relation (28), where g αβ =δ μ αδ ν β is now given by (52). The latter is a static spherically symmetric vacuum solution for the classical Einstein field equations with a cosmological constant : with R μν and R the Ricci tensor and scalar of general relativity. In accordance, we introduce a cosmological constant term to the field equations (7) in vacuum: The tensors R μν and R reduce to the standard Ricci tensor and scalar of general relativity for the metric (52) since the latter has no direct dependence on y. Additionally, from (8) in vacuum we get S = 0 the same way as in the previous paragraph. Therefore, using (53) in Eq. (54) we get relation (35) again. It is obvious that the procedure is the same as before and only the explicit form of g μν (x) changes. As such, we reach the same equation for A γ , namely relation (48) with g μν given by (52). Again, by separation of variables, one finds (see Appendix B):

Paths in the Schwarzschild-Randers spacetime
Now that we have A γ and hence the full metric, we can study particle trajectories in T M. A Lagrangian for point particles in the total space of the Lorentz tangent bundle has been proposed in [17]: L(x,ẋ, y) = (ag μνẋ μẋ ν + bδ α μ v αβẋ μ y β + cv αβ y α y β ) 1/2 (56) with a, b, c constants. The associated equations of motion are (g κν + zδ α κ δ β ν v αβ )ẍ κ + (γ νκλ + zσ νκλ )ẋ κẋ λ = 0 and with and z = −b 2 /4ac is a constant. The Christoffel symbols of the first kind for the metric g κν (x) are The termδ α κδ β ν v αβ is the metric of the v-space lowered down to the h-space via the generalized Kronecker symbols which are defined asδ μ α =δ α μ = 1 for a = μ+4 and equal to zero otherwise. 1 We will write for convenienceδ α κδ β ν v αβ = v κν and similarlyδ α κδ β ν w αβ = w κν . We define g κν = g κν + zv κν and we observe that its inverse is g μν = (1 + z) −2 (g μν + zv μν ) in the sense that g κν g μν = δ μ κ to first order in w μν . Contracting (57) with g μν gives g μν g κνẍ κ +(g μν γ νκλ + zg μν σ νκλ )ẋ κẋ λ = 0 where γ μ κλ = g μν γ νκλ and σ μ κλ = g μν σ νκλ . After some straightforward calculations, eq. (61) gives The horizontal part of the tangent vector on the paths iṡ x μ = dx μ /ds with s an affine parameter along the path defined as [17]: with s 0 , λ 0 and λ 1 constants and λ is an arbitrary parameter of the path. The sign of g μν (x, y) is determined by the tangent vector of the path, specifically if dx ν /dλ is timelike with respect to g μν (x, y) (g μν (x, y) dx μ dλ dx ν dλ < 0) then we get "−", likewise for a spacelike tangent vector with respect to g μν (x, y) (g μν (x, y) dx μ dλ dx ν dλ > 0) we get "+". The paths (62) will play the role for our model that the geodesics play for general relativity. As is the case for the latter, we need a classification of path segments with respect to their character i.e. timelike, null and spacelike. We define: -Timelike segment: g μν (x)ẋ μẋ ν < 0 at every point -Null segment: g μν (x)ẋ μẋ ν = 0 at every point -Spacelike segment: g μν (x)ẋ μẋ ν > 0 at every point Therefore, the character of the path is determined by the metric tensor g μν (x) of the horizontal subspace. We define the proper time τ as 1 In general,δ α μ andδ ν β can be used to lift an object from the horizontal to the vertical subspace of T T M or lower down one from the vertical to the horizontal subspace. This allows us to perform algebraic operations between components of tensors belonging to different sub-spaces of T T M.
where τ 0 is constant. By comparing relations (64) and (65) we see that the parameter s on the paths (62) cannot be written as an affine transformation of the proper time in general.
We remark that Eq. (62) reduces to the classic geodesics equation of general relativity when the perturbation w αβ goes to zero, as it should.

Timelike paths
To begin, we rewrite the perturbation (31) as where a = −g μν y μ y ν and we have lowered down A γ and y γ using the generalized Kronecker deltas. It is straightforward to show Using (63), (66) and (68) we calculatẽ Now, if we take into account that g μν u μ u ν = −1, the above relation gives Substituting the relations (58), (60), (66) and (70) into (62) we geẗ This is the generalized path equation for the timelike sector of the metric g μν (x).
Remark. If we set a = 1 at some fixed point then (71) can be written as with F κν = ∂ ν A κ − ∂ κ A ν the field strength tensor of A ν and we have set z 1+z := − e m where e the electric charge and m the mass of the particle. If we ignore the r.h.s of the above equation then (72) will have the same form as the equation of a charged particle subject to the Lorentz force with an electromagnetic vector potential A γ in the Riemannian setting. A similar equation which is derived from a Finsler-Randers Lagrangian and contains a Lorentz force term has been studied in [37]. However, in our more generalized setting we also get the r.h.s. perturbation term which depends on A ν and its first derivatives. Therefore, a possible relation between our Schwarzschild-Randers metric and the Lorentz force requires further investigation and goes beyond the scope of this work. Now, it is known that we can always approach a timelike path (geodesic) with a proper time parameter broken null path with the same endpoints [58]. In this approximation it is considered that the number of null path segments with infinitesimal distance between two neighboring points increases following the timelike path. Therefore, the final null path of zero length (with respect to g μν ) approaches the timelike path (71), however the parameter along them is replaced by an appropriate affine one.
Substituting to (71) the solution (51) we get the explicit form of the timelike paths components for r > R s : with f = 1 − R S r , R S the Schwarzschild radius.

Application
As an application for our model we present a numerical solution (using the differential equation solver of Mathematica) of the timelike path equations (73-76) for an appropriate choice of parameters and initial values. For this application we consider that θ = π 2 while t, r and φ are the variables. Below we present two figures which clarify the difference between Schwarzschild-Finsler-Randers (S-F-R) paths (red line) and general relativity (GR) geodesics (blue line). Notice, that the full analysis of timelike paths as well as the applications to Astrophysics will be studied in a forthcoming paper.
From the figures we observe that in the case of S-F-R model the timelike path reaches a higher maximum distance which is somewhat larger than the path provided by General Relativity. We notice that the time it takes for the S-F-R path to reach the Schwarzschild radius is more than in GR (Fig. 1). From these observations we see that in our model Fig. 1 This is an r, t graph of the timelike paths that we find using our theoretical model Schwarzschild-Finsler-Randers (S-F-R) in comparison to the geodesics of general relativity (GR) the maximum radial distance of the orbit is greater than that of GR and also the rate at which the particle falls is slower (Fig. 2). As expected, this deviation from the GR geodesic is produced from the right hand side of (73-76), where the extra terms in the S-F-R model act as a force that opposes gravity. Since the extra terms of the rhs of (73-76) are taken to be small the corresponding deviation from the GR geodesic is relatively small.

Spacelike paths
For completeness, we will find the spacelike paths from (62) following the same procedure as above. From (50) we get where, as before, we have lowered down w αβ to the horizontal space using the generalized Kronecker symbols and we have set u ν = y ν /a where a = g μν y μ y ν . Following the same steps as for the timelike section of g μν and taking into account that g μν u μ u ν = 1, (62) gives in this work constitute a first qualitative astrophysical application of our model. We plan to make a more detailed study on the S-F-R model's predictions in a future work in which the full data will be provided.] 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: Basic structures on the Lorentz tangent bundle
The The adapted basis on the total space T T M is defined as where N α μ are the components of a nonlinear connection. The curvature of the nonlinear connection is defined as The nonlinear connection induces a split of the total space T T M into a horizontal distribution T H T M and a vertical distribution T V T M. The above-mentioned split is expressed with the Whitney sum: The horizontal distribution or h-space is spanned by δ μ , while the vertical distribution or v-space is spanned by∂ α . Under a local coordinate transformation on the base manifold, the adapted basis vectors transform as: 2 The adapted dual basis of the adjoint total space T * T M is {E A } = {dx μ , δy α } with the definition The transformation rule for {dx μ , δy α } is: The bundle T M In this work, we consider a distinguished connection (dconnection) D on T M. This is a linear connection with coefficients {Γ A BC } = {L μ νκ , L α βκ , C μ νγ , C α βγ } which preserves by parallelism the horizontal and vertical distributions: From these, the definitions for partial covariant differentiation follow as usual, e.g. for X ∈ T T M we have the definitions for covariant h-derivative and covariant v-derivative A d-connection can be uniquely defined given that the following conditions are satisfied: We use the symbol D instead of D for a connection satisfying the above conditions, and call it a canonical and distinguished d-connection. Metric compatibility translates into the conditions: The coefficients of canonical and distinguished d-connection are (A.14) 2 The generalized Kronecker symbols are defined as:δ μ α =δ α μ = 1 for a = μ + 4 and equal to zero otherwise. The h-curvature tensor of the d-connection in the adapted basis and the corresponding h-Ricci tensor have, respectively, the components The v-curvature tensor of the d-connection in the adapted basis and the corresponding v-Ricci tensor have, respectively, the components The generalized Ricci scalar curvature in the adapted basis is defined as where R = g μν R μν , S = v αβ S αβ (A. 28) In the main text, we use the more convenient definitions After calculations we end up with c (5) = 0 and by substitution to (B.38) and (B.39) we find with k 5 being a constant. Also, like before, if we calculate the μ = 2, ν = 2, γ = 5 and μ = 3, ν = 3, γ = 5 we get no dependence on θ and φ. Therefore, we find withÃ 5 a constant of integration. If we put this solution in the μ = 1, ν = 1 equation, we getÃ 5 = 0. For μ = 0, ν = 0, γ = 6 we separate variables like before A 6 = R 6 (r )T 6 (t) and we get two equations: with c 2 (6) the separation constant. For μ = 0, ν = 1, γ = 6 like before we find that c (6) = 0 and by substitution to (B.43) we find that R 6 (r ) =R 6 r withR 6 a constant of integration. We setR 6 to zero to keep our solution finite at infinity, so we end up with A 6 = 0. For A 7 we set μ = 0, ν = 0, γ = 7 and we find the same equations as for A 6 . That leads to A 7 = 0 as well.
To sum up, we have found the following solution for A γ from Eq. Solution for the Schwarzschild-De Sitter-Randers spacetime We will modify the solution (51) and see if it satisfies (48) for the metric (52). An obvious ansatz is to replace the term 1 − R S r in (51) with 1 − R S r − 3 r 2 . Doing this we get which is verified to be a solution of (48) for the metric (52).