
==== Front
Eur Phys J C Part Fields
Eur Phys J C Part Fields
The European Physical Journal. C, Particles and Fields
1434-6044
1434-6052
Springer Berlin Heidelberg Berlin/Heidelberg

12150
10.1140/epjc/s10052-023-12150-1
Regular Article - Theoretical Physics
Generalized teleparallel de Sitter geometries
http://orcid.org/0000-0001-5482-6703
Coley A. A. 1
http://orcid.org/0000-0002-1046-439X
Landry A. a.landry@dal.ca

1
http://orcid.org/0000-0003-2403-3176
van den Hoogen R. J. 2
http://orcid.org/0000-0003-2677-8649
McNutt D. D. 3
1 https://ror.org/01e6qks80 grid.55602.34 0000 0004 1936 8200 Department of Mathematics and Statistics, Dalhousie University, Halifax, NS B3H 3J5 Canada
2 https://ror.org/01wcaxs37 grid.264060.6 0000 0004 1936 7363 Department of Mathematics and Statistics, St. Francis Xavier University, Antigonish, NS B2G 2W5 Canada
3 https://ror.org/00wge5k78 grid.10919.30 0000 0001 2259 5234 Department of Mathematics and Statistics, UiT: The Arctic University of Norway, 9037 Tromso, Norway
30 10 2023
30 10 2023
2023
83 10 97724 7 2023
14 10 2023
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ 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://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3. SCOAP3 supports the goals of the International Year of Basic Sciences for Sustainable Development.
Theories of gravity based on teleparallel geometries are characterized by the torsion, which is a function of the coframe, derivatives of the coframe, and a zero curvature and metric compatible spin-connection. The appropriate notion of a symmetry in a teleparallel geometry is that of an affine symmetry. Due to the importance of the de Sitter geometry and Einstein spaces within General Relativity, we shall describe teleparallel de Sitter geometries and discuss their possible generalizations. In particular, we shall analyse a class of Einstein teleparallel geometries which have a 4-dimensional Lie algebra of affine symmetries, and display two one-parameter families of explicit exact solutions.

http://dx.doi.org/10.13039/501100000038 Natural Sciences and Engineering Research Council of Canada W.F. James Chair of Studies in the Pure and Applied Sciences at St.F.XAtlantic Association of Research in Maths Science (AARMS)Norwegian Financial Mechanismissue-copyright-statement© EDP Sciences, Societa Italiana di Fisica (SIF) and Springer-Verlag GmbH, DE, part of Springer Nature 2023
==== Body
pmcIntroduction

Teleparallel geometries, which are characterized by the torsion (which is a function of the coframe and its derivatives), and a zero curvature and metric compatible spin-connection, provide an alternative to Riemmannian geometries in which to formulate a theory of gravity. A particular subclass of teleparallel gravitational theories is dynamically equivalent to General Relativity (GR) and is often referred to as the Teleparallel Equivalent to General Relativity (TEGR) [1]. In the covariant approach to the more general class of F(T) teleparallel gravity theories [2], the teleparallel geometry is defined in a gauge invariant manner as a geometry having a spin connection with vanishing curvature. The spin connection can be zero in a very special class of “proper” frames and is non-zero in all other frames [1, 2].

The torsion scalar, T, is defined in terms of the torsion tensor1 T=12TμνaSaμν,

where the superpotential, Sμνa, is constructed from the torsion tensor by2 Saμν=12Taμν+Taνμ-Taμν-haνTϕϕμ+haμTϕϕν.

The variations of the Lagrangian h2κF(T)+LMatt, which include a non-trivial spin-connection [2, 3], lead to Lorentz covariant field equations (FEs) and the resulting theory is locally Lorentz invariant [4]. Treating the spin-connection as an independent field, the Lagrangian can then be written using Lagrange multipliers to impose the two constraints of zero curvature and metric compatibility. Assuming an orthonormal frame, the corresponding variations with respect to (w.r.t.) the Lagrange multipliers lead to the following equations (eqns.)3 ωbμa=Λca∂μΛbcandω(ab)μ=0,

where Λba∈SO(1,3) (and Λbc≡(Λ-1)bc). Variations of the Lagrangian with respect to the coframe can be decomposed into symmetric and antisymmetric parts: 4a κΘ(ab)=F′′(T)S(ab)ν∂νT+F′(T)G~ab+12gabF(T)-TF′(T),

4b 0=F′′(T)S[ab]ν∂νT,

where G~ab is the Einstein tensor calculated explicitly from the metric [2]. Note that any quantity with an overtilde is computed using the Levi-Civita connection in terms of the metric. In the above FEs, Θ(ab) is the matter energy momentum tensor (and Θ[ab]=0 due to the invariance of the FEs under SO(1, 3)).

The determination of new solutions in a given F(T) theory is vital for investigating the validity of the theory. Furthermore, distinguishing a new solution from a known solution can be difficult from the perspective of coordinates, as this necessitates finding the coordinate transformation that brings one solution to the other. As an alternative, it is possible to use torsion invariants to determine the equivalence of two solutions. The Cartan–Karlhede algorithm is an algorithmic approach to generate a geometrically preferred frame and connection which can characterize a given solution in terms of invariant quantities [5].

The algorithm exploits the canonical form of the q-th order torsion tensors (that is the torsion tensor and its q-th covariant derivatives) to fix the frame. At the zeroth iteration, the torsion tensor is computed and the Lorentz frame freedom is used to adapt the frame so that the torsion tensor takes on a canonical form. The number of functionally independent components, t0, are recorded and the remaining frame freedom, which does not affect the torsion tensor, is denoted as the 0-th order linear isotropy group, H0.

For the q-th iteration, the q-th order torsion tensor is computed and the q-1 linear isotropy group, Hq-1 is used to determine a canonical form for the q-th torsion tensor and the number of functionally independent components, tq, are recorded along with the q-th linear isotropy group, Hq. The algorithm stops at iteration p+1, when tp=tp+1 and Hp=Hp+1. The subgroup Hp of the Lorentz group is then called the linear isotropy group of the geometry.

The resulting invariant quantities that characterize the geometry are then {tq}q=0p+1, {Hq}q=0p+1 and the set of components T={Tbca,∇d1Tbca,…,∇dp+1…∇d1Tbca}. The set T are torsion invariants known as Cartan invariants. If two solutions share the same discrete sequences, with Cartan invariants T and T′, respectively, equivalence of the solutions can be proven in a non-algorithmic way by solving the equations arising from equating the components of T and T′.

In teleparallel geometries, the tetrad (or (co)frame) and the corresponding spin-connection replace the metric as the primary objects of study. Consequently, the appropriate notion of symmetry in a teleparallel geometry is that of an affine symmetry (and not a metric symmetry; i.e, a Killing vector). An affine frame or intrinsic symmetry on the frame bundle is a diffeomorphism from the manifold to itself which is characterized by the existence of a vector field, X, satisfying [5, 6]:5 LXha=λabhbandLXωbca=0,

where ωbca denotes the spin-connection relative to the geometrically preferred invariantly defined frame ha determined by the Cartan–Karlhede algorithm and λab is an element of the linear isotropy group determined by the algorithm.

There do not exist any teleparallel geometries admitting a maximal group of affine frame symmetries other than Minkowski space [5, 7]. If a four-dimensional teleparallel geometry has a non-zero torsion, then the maximum dimension of the group of affine symmetries is at most seven [5]. Therefore, in analogy with de Sitter geometries in GR, we define the teleparallel de Sitter (TdS) geometry as that nontrivial teleparallel geometry which has a seven dimensional group of symmetries that is also a subgroup of the group of the Killing symmetries of the de Sitter metric [8].

Due to the importance of the de Sitter geometry within GR, we shall study TdS geometries and their generalizations in this paper. After a brief introduction to the TdS geometry, we discuss generalizations of this geometry mathematically and describe some physical applications, with particular emphasis on possible analogues of Einstein spaces. In particular, we shall discuss a class of Einstein teleparallel geometries which have a 4-dimensional Lie algebra of affine symmetries (in the case where the Einstein space parameter λ is equal to zero). We display the governing symmetric FEs explicitly in the first of the two solutions of the antisymmetric FEs. We shall investigate power law solutions and display two one-parameter families of explicit solutions (where we can implicitly find the form of F(T) for each parameter value).

Teleparallel “de Sitter” (TdS)

We will work in a coframe in which the tangent space metric has the form gab=ηab=Diag[-1,1,1,1], which still allows a O(1, 3) subgroup of GL(4,R) of residual gauge transformations which leaves the metric gab=ηab invariant. We can then restrict attention to proper ortho-chronous Lorentz subgroups, SO(1, 3) or SO(1,3)+. Within this orthonormal gauge choice, the resulting field equations transform homogeneously under the remaining O(1, 3) (or SO(1, 3) or SO(1,3)+) Lorentz gauge transformations.

Using the Cartan–Karlhede algorithm [5], we can determine the G7 geometries obtained by requiring that the Cartan invariants are all constant. This follows from the formula for the dimension of the affine frame symmetry group, N=s+4-tp, where s is the dimension of the linear isotropy group and tp is the number of functionally independent invariants at the conclusion of the algorithm. If all of the Cartan invariants are constant, then tp=0 and the dimension of the linear isotropy group is three, yielding a seven-dimensional affine frame symmetry group.

Solving the differential equations arising from the requirement that the only non-trivial components of the Cartan invariants are Tabc yields the general case in which a(t)=A0eH0t is the scale factor, and H0 is a non-zero constant in the teleparallel k=0 Robertson–Walker geometry [9]. The special teleparallel geometry case where H0=0, a(t)=A0 a non-zero constant, k=±1 and the additional affine frame symmetry X7=∂t can be considered as the analogue of the Einstein static geometry in GR [8] . In the general case (H0≠0), the affine frame symmetry is of the form6 X7=-1H0∂t+r∂r.

The resulting Lie algebra of {XI}I=17={X1,X2,X3,X4,X5,X6,X7} is given by[X1,X5]=X3,[X1,X6]=X2,[X1,X7]=X1,[X2,X4]=-X3,[X2,X6]=-X1,[X2,X7]=X2,[X3,X4]=X2,[X3,X5]=-X1,[X3,X7]=X3,[X4,X5]=-X6,[X4,X6]=X5,[X5,X6]=-X4.

This is a subalgebra of the Lie algebra for the group of metric (Killing) symmetries of de Sitter metric. We therefore define the Teleparallel de Sitter (TdS) geometry as the teleparallel geometry with a G7 Lie group of affine symmetries which is the semi direct product of the one-dimensional subgroup of O(1, 4) and the six dimensional Euclidean group E(3). Note in this geometry the covariant derivative of the torsion tensor is zero.

Properties

The diagonal de Sitter co-frame is7 ha=dteH0tdreH0trdθeH0trsin(θ)dϕ,

which has the corresponding spin-connection one form8 ωba=000000-dθ-sin(θ)dϕ0dθ0-cos(θ)dϕ0sin(θ)dϕcos(θ)dϕ0,

or the non-trivial components of the spin-connection are (all indices down)9 ω23θ=-ω32θ=-1,ω24ϕ=-ω42ϕ=-sin(θ),ω34ϕ=-ω43ϕ=-cos(θ)

(alternatively, one could use the proper de Sitter coframe given in [8]). Using Eqs. (7) and (8) we have that10 Ta=H00h1∧h2h1∧h3h1∧h4.

The non-trivial components of the torsion tensor are11 Ttr2=-Trt2=H0eH0t,Ttθ3=-Tθt3=H0eH0tr,Ttϕ4=-Tϕt4=H0eH0trsin(θ).

The torsion tensor is decomposable into three irreducible parts [5]:12 Tabc=23tabc-tacb-13gabVc-gacVb+ϵabcdAd,

where the vector, axial and tensor components are, respectively: 13a Va=Tbab,

13b Aa=16ϵabcdTbcd,

13c t(ab)c=12Tabc+Tbac-16gcaVb+gcbVa+13gabVc.

The dual superpotential two form ∗Sa is14 ∗Sa=0-2H0e2H0tr2sin(θ)dθ∧dϕ2H0e2H0trsin(θ)dr∧dϕ-2H0e2H0trdr∧dθ.

The non-trivial components of ∗Sa are15 ∗Sθϕ2=-∗Sϕθ2=-2H0e2H0tr2sin(θ),∗Srϕ3=-∗Sϕr3=2H0e2H0trsin(θ),∗Srθ4=-∗Sθr4=-2H0e2H0tr,

or the non-dual non-trivial components are16 S212=-S221=-2H0,S313=-S331=-2H0,S414=-S441=-2H0.

The part of the superpotential that comes into the FEs in terms of coordinates is S(ab)μ and S[ab]μ. The non trivial components of S[ab]μ are17 S[12]r=-S[21]r=H0eH0t,S[13]θ=-S[31]θ=H0eH0tr,S[14]ϕ=-S[41]ϕ=H0eH0trsin(θ).

The non trivial components of S(ab)μ are18 S(12)r=S(21)r=-H0eH0t,S(13)θ=S(31)θ=-H0eH0tr,S(14)ϕ=S(41)ϕ=-H0eH0trsin(θ),S22t=S33t=S44t=-2H0.

The torsion scalar is T=6H02, the axial part of the torsion scalar is zero, and the magnitude of the vector part is -9H02 (k=0,T=6H02; for reference, the vector part is V2=-9H02 and the axial part A2=0).

Assuming a comoving perfect fluid as the matter source, the FEs in the TdS geometry reduce to19 κρ=-κP=-12F(T0)+6F′(T0)H02,

where necessarily ρ and P are constants. The equations formally reduce to their GR counterparts only when F(T)=T. We note that the effective equation of state20 ωeff=Pρ=-1,

is the same as its GR counterpart; however, the effective cosmological constant Λeff≡κρ depends in principle on the two parameters F(T0) and F′(T0). The energy-momentum conservation eqn.21 ρ˙+3a˙a(ρ+P)=0,

is identically satisfied.

Applications

If T=const., then the FEs for F(T) teleparallel gravity are equivalent to a rescaled version of TEGR (which looks like GR with a cosmological constant and a rescaled coupling constant) [2]. In the case of TEGR, where F(T)=T, Eq. (4b) vanishes. For F(T)≠T, the variation of the gravitational Lagrangian by the flat spin-connection is equivalent to the antisymmetric part of the FEs in Eq. (4b) [2]. In addition, since the canonical energy momentum is symmetric, Θ[ab]=0, then the antisymmetric part of the FEs (4b) limit the possible solutions of spacetimes.

Cosmological models in flat (k=0) TRW models have been studied in [10, 11] (also see references within). For example, simple power-law scale factor solutions in specific F(T) models (using a priori ansatz such as a polynomial functions) have been investigated. Futhermore, the late time behaviour of the usual de Sitter spacetime, and in particular its stability (as a fixed point), has been studied in [10] (see also [12]).

Stability

Let us simply consider the linear perturbations of the scalar quantities T,ρ,P in the TdS solution: 22a T=6H02+T1,

22b κρ=-F(T0)2+6F′(T0)H02+κρ1,

22c κP=--F(T0)2+6F′(T0)H02+βκρ1,

where we have used the zeroth order TdS expressions above with T0≡6H02, and we assume the parameter β≥-13 from the energy conditions.

To first order the antisymmetric FEs yield23 S0[ab]γT1,γ=0,

where the non trivial zeroth order components of S[ab]μ are given (by Eqs. (22a) to (22c)) above, so that trivially T1,γ=0 and hence T1=T1(t).

The FEs and the conservation law then imply that 24a -12F′(T0)+6F′′(T0)H02T1=κρ1,

24b -3H0(1+β)ρ1=ρ1˙.

We immediate have that ρ1=ρ1(t), so that from the conservation eqn. ρ1→0 to the future (since H0>0), whence (except in the degenerate case in which the term in square brackets above might be zero), T1→0, and hence T→T0, a constant, so that the resulting geometry is formally equivalent to that of GR and the corresponding stability results follow.

The function F(T)

A number of particular examples of F(T) theories studied in the literature include polynomial functions in T, and especially quadratic T2 theory [10, 11]. Recently, it has been shown that the theory [13]25 F(T)=-Λ+T+γT2,

can alleviate a variety of cosmological tensions. It is also of interest to study the theory with26 F(T)=T+γTβ~,

where β~>0 ( where β~ is potentially a positive integer greater than 2 for sufficient differentiability), or an exponential function27 F(T)=F(0)exp2γT-1∼F(0)T+γT2+...

for small T.

Analogues of Einstein spaces

In Riemannian geometries the subclass of Einstein spaces with R~ab=λgab, and which contain the de Sitter geometry as a special case, are of interest [14]. We note that for Einstein spaces, in which the algebraic structure of the Ricci tensor R~ab is special, the frame components of the Ricci tensor R~ab are constants and its covariant derivative is zero.

It may be of interest to study analogous particular teleparallel geometries with special properties (but non-constant T) which include the TdS geometry. Possible examples might include teleparallel geometries with special algebraic, topological or symmetry properties. Such examples often lead to R~ab=λgab and hence to a subcase of the usual Einstein spaces (which implies the energy momentum tensor Θab=λgab).

Algebraic

Examples include: The components of the torsion tensor are constant or its covariant derivative is zero.

The algebraic properties of some tensor constructed from the torsion tensor are special.

In addition, the torsion tensor can be invariantly decomposed algebraically into its vectorial part Va (Eq. (13a)), axial part Aa (Eq. (13b)) and tensorial part Tabc (Eq. (13c)). Special invariantly defined subcases occur when the axial and tensorial parts are zero (or tensorial and vector parts are zero). For example, the TdS geometry has vanishing axial part zero [6].

Topological

Special classes of spacetimes have topological properties such as product (or warped product) structures, and hence decomposable and reducible structures [15, 16]. Unlike the usual Riemannian case, it is the frame and spin-connection (and/or the torsion tensor) that should have these splitting properties (rather than the curvature and metric).

In the Riemannian case, the de Rham decomposition theorem states that if the holonomy group of a simply-connected Riemannian manifold preserves a proper subspace of the tangent space (i.e., is reducible), then the tangent space is decomposable (into holonomy invariant subspaces) and the manifold is (locally) isometric to a product manifold (i.e., a Riemannian manifold is locally a product of Riemannian manifolds with irreducible holonomy algebras). The Lorentzian case was studied in [17], utilizing Wu’s theorem [18, 19], that asserts that every simply-connected, complete semi Riemannian manifold is isometric to a product of simply-connected, complete semi-Riemannian manifolds, of which one can be flat and all others are indecomposable or “weakly-irreducible” (i.e., with no nondegenerate invariant subspace under holonomy representation).

Symmetries

The teleparallel de Sitter (TdS) geometry has a G7 Lie group of affine symmetries which is the semi direct product of a one-dimensional subgroup O(1, 4) and a six dimensional Euclidean group E(3). We could define an appropriate generalization of TdS by considering spaces with a subalgebra of affine symmetries. For example, if we consider the subalgebra generated by the 3 rotational affine symmetries X4-X6 (spherical symmetry) and the special affine symmetry X7=-1H0∂t+r∂r (Eq. (6)), we obtain a subclass of geometries with a 4-dimensional Lie subgroup as its affine symmetry group, which we shall define as Einstein Teleparallel (ET) geometries. The TdS geometry is then a special case in which there are an additional 3 affine symmetries (X1-X3) [6, 9].

Einstein teleparallel geometries

Relative to the representation for the isotropy group used in [6], we can determine the most general frame and spin-connection with the three spherically symmetric affine frame symmetry generators X4-X6 along with the fourth additional affine symmetry X7 defined above (see [6] and also the case λ=0 in [9]). We can also choose a new coordinate system to “diagonalize” the frame. The resulting veilbein is:28 hμa=A1(t,r)0000A2(t,r)0000A3(t,r)0000A3(t,r)sin(θ).

With this choice of invariant symmetry frame, we can now obtain the most general metric compatible connection [6]:29 ω341=X1(t,r),ω342=X2(t,r),ω233=ω244=X3(t,r),ω234=-ω243=X4(t,r),ω121=X5(t,r),ω122=X6(t,r),ω133=ω144=X7(t,r),ω134=-ω143=X8(t,r),ω344=-cos(θ)A3sin(θ).

Finally, to determine the most general connection for a teleparallel geometry we must impose the flatness condition. The resulting equations can be solved so that any spherically symmetric teleparallel geometry is defined by (the three arbitrary functions in the veilbein (28) along with) the following spin-connection components [9]:30 X1=-∂tχA1,X2=-∂rχA2,X3=cosh(ψ)cos(χ)A3,X4=cosh(ψ)sin(χ)A3,X5=-∂tψA1,X6=-∂rψA2,X7=sinh(ψ)cos(χ)A3,X8=sinh(ψ)sin(χ)A3,

where χ and ψ are arbitrary functions of the coordinates t and r.

Thus, to determine the conditions of the inclusion of the new symmetry generator, we need only examine the conditions LX7ha=0, from which we obtain:31 A1(t,r)=f1(z),A2(t,r)=f2(z)r,A3(t,r)=f3(z),

where z≡re-C0t and C0≡-H0 is a non-zero constant. The Lie derivatives of the metric g and the spin-connection, respectively, yield LX7g=0 and LX7ωbca=0, as desired. It therefore follows that the arbitrary functions in the spin-connection’s components must be functions of the form:32 χ(t,r)=χ(z),ψ(t,r)=ψ(z).

For general situations in which the parameter λ≠0, we have that the torsion scalar T=T(z,r) will also contain terms in rλ as studied in [9] for T(z,r)=T0(z)r2λ. From the FEs we obtain the following equation in term of r-2λF′(T):33 κP+ρ=2-g4h2+2h3T0′(z)+2m3T0(z)h4T0′(z)+m4T0(z)+g2+g3r-2λF′(T),≡2H(z)r-2λF′(T).

We note that in general with λ≠0, H(z)≠0 [9].

For λ=0, as considered here, we obtain that T(z,r)=T(z). In general, without substituting in the solutions of the antisymmetric FEs solutions (and the symmetric FEs), we have that H(z)≠0. However, we will see below that, using the solution to the antisymmetric FEs, we find that H(z)=0. This then implies that ρ+P=0. This is consistent with the fact that this is intended as an anologue of an Einstein space.

Antisymmetric FEs and solutions

For λ=0, the antisymmetric FEs become (when the torsion scalar T(z)≠const and F(T) is not linear in T): 34a 0=sinχ(z)C0f2(z)sinhψ(z)+f1(z)coshψ(z),

34b 0=cosχ(z)C0f2(z)coshψ(z)+f1(z)sinhψ(z).

There are two solutions: A: cosχ(z)=0: sinχ(z)=δ=±1, χ=k+12π and tanhψ(z)=-f1(z)C0f2(z).

B: sinχ(z)=0: cosχ(z)=δ=±1, χ=kπ and tanhψ(z)=-C0f2(z)f1(z).

General form of symmetric FEs

The general symmetric FEs can be expressed in the following form: 35a k1∂zln(F′(T))=g1,

35b κP-F(T)2=2F′(T)-k2∂zln(F′(T))+g2,

35c κρ+F(T)2=2F′(T)-k3∂zln(F′(T))+g3,

35d k4∂zln(F′(T))=g4,

where the functions gi and ki (i=1-4) are explicitly displayed below for solution A (the functions for solution B are displayed in the Appendix). By combining the pair of Eqs. (35b) and (35c) and then Eqs. (35a) and (35d), respectively, we obtain: 36a ∂zln(F′(T))=g1k1=g4k4,

36b κP+ρ=2F′(T)-k2+k3∂zln(F′(T))+g2+g3≡F′(T)H(z).

The general torsion scalar expression respecting Eqs. (36a) and (36b) is the following:37 T(z)=1f12f22f32[-4f1f2zcosχ(f3f1′+f3f2ψ′coshψ+sinhψψ′f1f3+C0f2f3′)+4zsinχcoshψχ′f12f2f3+4C0zsinχsinhψχ′f1f22f3-2f12z2f3′2+f22-4z2f1f3f1′f3′+4C02z2f2f3′f2′f3+f2f3′2].

Coordinates

There is a class of transformations that maintain the “diagonal” form of the frame in (28). However, in general, this will also lead to a change in the form of the connection and the FEs. Therefore, in order to maintain the comoving nature of the perfect fluid source, we restrict ourselve to transformations of the form t→f(t) and r→g(r) (i.e., redefining the r and t coordinates). This may simplify the eqns. further. Even then, the explicit forms for the symmetry vectors and the similarity variable z (for example) does change. We shall not change coordinates here.

However, since the static case corresponds formally to C0=0, we assume explicitly that C02≠0 here. In this case we can then effect a simple time change to set C02=1, which we shall do hereafter. This has the effect of changing the metric, -dt2→-1H02dt2. In the exact TdS solutions C02=1 and the metric changes accordingly.

Solution A

The symmetric FE components (see Eqs. (35a) to (35d)) for solution A (with C02=1) are explicitly given by (the symmetric FEs components for solution B are displayed in the Appendix): 38a g1=[-z2f1f2f3f3′′f12+f22-z2f12f2f32f1′′+z2f1f22f32f2′′+z2f3′2-f1f23+f13f2+f3f3′f22+f12z2f2f1′+f1z2f2′-zf2+f32f1′f12-f22z2f2′-zf12f2+f1f22zf32f2′-f12f2],

38b g2=1f13f22f32[-z2f1f22f3f3′′+z2f1f3′2f12-f22+f3f3′(z2f1′f22+2f12-f1f2z2f2′+zf2)],

38c g3=1f12f23f32[-z2f3′′f3f12f2+f2f3′2z2f22-f12-f3′(-z2f3f2′2f22+f12+f1f2z2f3f1′+zf1f3)],

38d g4=[lnf1f2zf3′]′,

39a k1=z2f1f2f3[f3′f12+f22+f3f1f1′-f2f2′],

39b k2=z2f3′f12f3,

39c k3=z2f3′f22f3,

39d k4=1.

We simplify the symmetric FEs for solution A into Eq. (36a), to obtain:40 [-z2f1f2f3f3′′f12+f22-z2f12f2f32f1′′+z2f1f22f32f2′′+z2f3′2-f1f23+f13f2+f3f3′f22+f12z2f2f1′+f1z2f2′-zf2+f32f1′f12-f22z2f2′-zf12f2+f1f22zf32f2′-f12f2][z2f1f2f3×[f3′f12+f22+f3f1f1′-f2f2′]]-1=-lnzf3′f1f2′.

Using the eqns. above and the symmetric FE, we find that the right-hand side of Eq. (36b) vanishes, thus giving P+ρ=0. So, for λ=0, we only have κP(t,r)=-κρ(t,r)≡-Λ0, regardless of the form of F(T(z)). This is consistent with the fact that this is intended as an anologue of an Einstein space. This means that we have the following relations:41 k2+k3g4=g2+g3,-lnzf3′f1f2′=g4=g2+g3k2+k3.

By putting the Eqs. (40) and (41) together, we arrive at the following superrelation:42 -lnzf3′f1f2′=g4=g2+g3k2+k3=g1k1.

The torsion scalar expression for solution A simplifies and is given explicitly by:43 T(z)=2z2f3′2(z)f22(z)-f12(z)+4z2f3(z)f3′(z)f2(z)f2′(z)-f1(z)f1′(z)-2f12(z)f22(z)f12(z)f22(z)f32(z).

Explicit equations: summary

First we satisfied the Antisymmetric FEs (Eqs. (34a) and (34b)) leading to solutions A and B. Then we set C02=1 in all of our physical quantities. For the symmetric FEs, we have Eqs. (35a) to (35d) with gi and ki (i=1-4) given by Eqs. (38a) to (39d), respectively. We restrict to solution A hereafter, whence Eq. (36a) becomes Eq. (40) (see also Eqs. (41) and (42)). After finding H(z)=0, Eq. (36b) implies that ρ+P=0. Then the torsion scalar T(z) is given explicitly by Eq. (43). The FE described by Eq. (35d), using k4=1 and the expression above for g4, is replaced by the first integral44 dFdT=B0f1f2zf3′,

which is a function of z, and where B0 is a constant. The remaining FE described by Eq. (35c) yields an eqn. for F(T(z)) in terms of z. We have that κρ=-κP=Λ0, where we have used solution A and set C02=1. In principle, the three remaining FEs, Eqs. (35c), (40) and (44) are then 3 ODEs for the 3 functions fi(z) for a given F(T(z)).

For example, for the quadratic function F(T) explicitly given by Eq. (25), we find that F′=1+2γT,F″=2γ which, when substituted into Eq. (44), yields45 T=-12γ+B02γf1f2zf3′,

whence using (43) yields a first order ODE for the fi. Equation (35c) then yields a quadratic eqn. for T in terms of functions of z (including a term dTdz) which upon using Eq. (43) yields a second order ODE for the fi. Equation (40) constitutes another second order ODE for the fi.

We note again that all potential solutions considered here (other than TdS explicitly) are not analogues of GR solutions and are consequently new.

Power law solutions

To find solutions to the symmetric part of the FEs describing an Einstein teleparallel geometry, we use the following ansatz:46 f1(z)=a0za,f2(z)=b0zb,f3(z)=c0zc,

where a0,b0,c0≠0. Using the solution A, from Eq. (35c) we obtain:47 Λ0=-F(T(z))2+2F′(T(z))[c(2b+c)a02z-2a-c(2a+c)b02z-2b].

We need the exact expressions for T(z) and F(T(z)) before simplifying further. Using g1=g4k1, we obtain from Eq. (40):48 0=za+b[a02c02c2-2a2+acz2a+2c-b02c02c2-2b2+bcz2b+2c-a02b02z2a+2b].

For the torsion scalar, we obtain:49 T(z)=2[c(2b+c)a02z-2a-c(2a+c)b02z-2b-1c02z-2c].

The derivative T′(z) is:50 T′(z)=-4c(2b+c)aa02z-2a-1-(2a+c)bb02z-2b-1-1c02z-2c-1.

From Eq. (44) we obtain:51 dF(T)dT=B1cza+b-c,

where B1 is a constant (the degenerate case c=0 is not included). For c=a+b we get exactly the TEGR solution. We will consider a+b≠c for a non-trivial F(T) solution.

We multiply by the derivative of the torsion scalar to obtain:52 dF(T(z))dz=dF(T)dTT′(z)=-4B1(c+2b)aa02zb-c-a-1-(c+2a)bb02za-c-b-1-1c02za+b-3c-1.

From this equation, we may integrate dF(T(z))dz w.r.t z (depending on the power of z – some situations may lead to z-1 and care should be taken). From Eq. (48), we can solve the symmetric FEs for a, b and c and then verify Eq. (47) by taking into account Eqs. (49) to (52) for each possible situation. In principle, there may be some free parameters remaining.

The possible solutions

If we consider Eq. (48) only there are, in principle, the following possibilities: a=b (with a02b02≠0), which is not possible because we need to satisfy: 53 0=c2-2a2+acc02[a02-b02]and0=a02b02.

This case includes the special case a=b=c.

a=c: which leads to: 54 0=-[c02c2-2b2+bc+a02]b02,

which gives the following relation between a and b: 55 (a+2b)(a-b)=-a02c02,

where, according to Eq. (55), -2b<a<+b for viable solutions.

b=c: which leads to 56 0=[c02b2-2a2+ab-b02].

Equation (56) leads to the following relation between a and b: 57 (b+2a)(b-a)=b02c02,

where it follows that -b2<a<+b for viable solutions.

Case I: a=c

For a=c and a02c02=-(a+2b)(a-b) (Eq. (55)), we then compute Eqs. (49), (51) and (52): 58a T(z)=2[(a+2b)(2a-b)z-2aa02-3a2b02z-2b],

58b dF(T)dT=F′(T(z))=B1azb,

58c dF(T(z))dz=-4B1(a+2b)(2a-b)a02zb-2a-1-3abb02z-b-1.

When we attempt to integrate Eq. (58c), we see that there may be special cases (b=0 or a=b2) which yield a z-1 power term, but in each case the corresponding coefficient is zero and therefore no power of z-1 terms ever occurs in Eq. (58c). From this, we can integrate Eq. (58c) giving the relation:59 F(T(z))=-4B1-(a+2b)a02zb-2a+3ab02z-b+C,

where C is a constant. By substituting Eq. (59) into Eq. (47), we find that:60 Λ0+C2=2B1(a+2b)a02-(a+2b)a02zb-2a=0.

This means that we have C=-2Λ0 in Eq. (59) (b≠0). This constitutes a one parameter (b say) family of solutions. It is possible to locally find the inverse z=z(T) from Eq. (58a) to explicitly find F(T).

As a simple albeit special example, if a=b2, Eqs. (58a) and (58b) become: 61a zb=32bb0-T,

61b dF(T)dT=24B12b0-T≡-F12-T,

where F1 is a constant. We integrate Eq. (61b) w.r.t. to T to obtain:62 F(T)=F1-T+C,

for T≤0. By substituting Eq. (62) into Eq. (47), we get that:63 C=-5bb0F12a0223-2Λ0≡F0-2Λ0,

where F0 is also a constant.

Case II: b=c

For b=c and b02c02=(b+2a)(b-a) (Eq. (57)), we must evaluate Eqs. (49), (51) and (52) (for a≠0): 64a T(z)=2[3b2a02z-2a-(2b-a)(2a+b)b02z-2b],

64b dF(T)dT=F′(T(z))=B1bza,

64c dF(T(z))dz=-4B13baa02z-a-1-(2b-a)(2a+b)b02za-2b-1.

When we attempt to integrate Eq. (64c), we see that there may be special cases when a=0 or a=2b which yield z-1 power terms, but in each case the corresponding coefficient is zero and therefore such a power of z-1 never occurs in Eq. (64c). We can consequently integrate Eq. (64c), giving the relation:65 F(T(z))=4B13ba02z-a-(2a+b)b02za-2b+C.

By substituting Eq. (65) into Eq. (47), we find that:66 Λ0+C2=-2B1(2a+b)b02-(2a+b)b02za-2b=0.

This again implies that C=-2Λ0 as an integrating constant in Eq. (65). In general, it is possible to locally find the inverse z=z(T) from Eq. (64a) then use it to compute F(T). Note that from Eqs. (64a) and (65), in the TEGR case we must have a=0 (which corresponds to the special case a+b=c), where67 z-2b=3b022a02-b024b2T.

Concluding remarks

Since the de Sitter geometry (and Einstein spaces in more generality) within GR have a number of important mathematical and physical applications, we have studied TdS geometries and their generalizations in theories of gravity based on teleparallel geometries. In particular, we have investigated a class of Einstein teleparallel geometries which have a 4-dimensional Lie algebra of affine symmetries. We solved the resulting antisymmetric FEs and we displayed all of the remaining governing eqns. explictly. We then investigated power law solutions and displayed two one-parameter families (within solution class A) of explicit Einstein teleparallel solutions (where we can implicitly find the form of F(T) for each parameter value). Very few explicit non-trivial exact solutions are known within teleparallel gravity, and the two one-parameter families of exact power law solutions obtained (being possible analogues of Einstein spaces) may have a number of important applications in cosmology and astrophysics.

Expressions for solution B

The functions in the symmetric FEs (Eqs. (35a) to (35d)) for solution B are explicitly given by:68 g1=[-z2f1f2f3f3′′f12+f22-z2f12f2f32f1′′+z2f1f22f32f2′′+z2f3′2-f1f23+f13f2+f3f3′f12+f22z2f1f2′-zf1f2+f32f1′z2f2′f12-f22-zf12f2+f1f22zf32f2′-f12f2],g2=1f13f22f32[-z2f1f22f3f3′′+z2f3′2-f1f22+f13+2f3′(z2f3f1′f222+f12-f1f22[z2f3f2′-δz(f12-f221-f22f12)+zf2f3])+δzf1f2f3f1f1′-f2f2′1-f22f12],g3=1f12f23f32[-z2f12f2f3f3′′+f3′2z2f23-f12f2-f3′[-z2f3f2′2f22+f12+f1f2(z2f1′f3+[δzf2(f12-f22f11-f22f12)+zf1f3])]-δzf1f22f3(f1f1′-f2f2′f11-f22f12)],g4=[lnf1f2zf3′]′,

69 k1=z2f1f2f3f3′f12+f22+f1f2f3z2f1f3f1′-f2z2f3f2′-δzf12-f221-f22f12,k2=zf12f3δf21-f22f12+zf3′,k3=zf22f3δf21-f22f12+zf3′,k4=1zf3′δf21-f22f12+zf3′.

From Eq. (36a) we explictly get:70 z2f1f2f3f3′f12+f22+f1f2f3z2f1f3f1′-f2×z2f3f2′-δzf12-f221-f22f12-1×[-z2f1f2f3f3′′f12+f22-z2f12f2f32f1′′+z2f1f22f32f2′′+z2f3′2-f1f23+f13f2+f3f3′f12+f22z2f1f2′-zf1f2+f32f1′z2f2′f12-f22-zf12f2+f1f22zf32f2′-f12f2]=1+δf2zf3′1-f22f12-1[lnf1f2zf3′]′.

Using Eqs. (68) and (69) and the symmetric FE, we find that the right-hand side of Eq. (36b) vanishes, thus giving P+ρ=0. So, for λ=0, we have P(t,r)=-ρ(t,r)≡-Λ0, regardless of the form of F(T(z)). This is consistent with the fact that this is intended as an anologue of an Einstein space. This means that we have the following relation:71 1+δf2zf3′1-f22f12-1lnf1f2zf3′′=g4k4=g2+g3k2+k3.

By putting the Eqs. (70) and (71) together, we arrive at the following superrelation:72 1+δf2zf3′1-f22f12-1lnf1f2zf3′′=g4k4=g2+g3k2+k3=g1k1.

The expression for the torsion scalar in the second solution then simplifies further.

Acknowledgements

AAC and RJvdH are supported by the Natural Sciences and Engineering Research Council of Canada. RJvdH is also supported by the W.F. James Chair of Studies in the Pure and Applied Sciences at St.F.X. AL is supported by an AARMS fellowship. DDMcN is supported by the Norwegian Financial Mechanism 2014-2021 (project registration number 2019/34/H/ST1/00636).

Data Availability Statement

This manuscript has no associated data or the data will not be deposited. [Authors’ comment: There is no data as such; everything is done analytically and displayed].
==== Refs
References

1. Aldrovandi R Pereira JG Teleparallel Gravity, Fundamental Theories of Physics 2013 Dordrecht Springer
2. Krssak M van den Hoogen RJ Pereira JG Boehmer CG Coley AA Teleparallel theories of gravity: illuminating a fully invariant approach Class. Quantum Gravity 2019 36 18 183001 10.1088/1361-6382/ab2e1f
3. Krssak M Saridakis EN The covariant formulation of f(T) gravity Class. Quantum Gravity 2016 33 115009 10.1088/0264-9381/33/11/115009
4. Krssak M Pereira JG Spin connection and renormalization of teleparallel action Eur. Phys. J. C 2015 75 519 10.1140/epjc/s10052-015-3749-2
5. Coley AA van den Hoogen RJ McNutt DD Symmetry and equivalence in teleparallel gravity J. Math. Phys. 2020 61 072503 10.1063/5.0003252
6. D.D. McNutt, A.A. Coley, R.J. van den Hoogen, A frame based approach to computing symmetries with non-trivial isotropy groups. J. Math. Phys. 64, 032503 (2023). arXiv:2302.11493 [gr-qc]
7. Hohmann M Järv L Krssak M Pfeifer C Modified teleparallel theories of gravity in symmetric spacetimes Phys. Rev. D 2019 100 084002 10.1103/PhysRevD.100.084002
8. Coley AA van den Hoogen RJ McNutt DD Symmetric teleparallel geometries Class. Quantum Gravity 2022 39 22LT01 10.1088/1361-6382/ac994a
9. A.A. Coley, R.J. van den Hoogen, A. Landry, D.D. McNutt, Preprint: Spherically symmetric teleparallel geometries (in preparation)
10. Bahamonde S Dialektopoulos KF Escamilla-Rivera C Farrugia G Gakis V Hendry M Hohmann M Said JL Mifsud J Di Valentino E Teleparallel gravity: from theory to cosmology Rep. Prog. Phys. 2023 86 026901 10.1088/1361-6633/ac9cef
11. Cai Y-F Capozziello S De Laurentis M Saridakis EN f(T) teleparallel gravity and cosmology Rep. Prog. Phys. 2016 79 106901 10.1088/0034-4885/79/10/106901 27599606
12. Bohmer CG Jensko E Lazkoz R Cosmological dynamical systems in modified gravity Eur. Phys. J. C 2022 82 6 500 10.1140/epjc/s10052-022-10412-y
13. E.N. Saridakis, Solving both H0 and σ8 tensions in f(T) gravity. (2023). arXiv:2301.06881 [gr-qc]
14. Besse AL Einstein Manifolds 1987 Berlin Springer Verlag
15. Kobayashi S Nomizu K Foundations of Differential Geometry and its Applications 1963 New York Interscience
16. Hall GS Symmetries and Curvature Structure in General Relativity 2004 Singapore World Science
17. Brannlund J Coley AA Hervik S Supersymmetry, holonomy and Kundt spacetimes Class. Quantum Gravity 2008 25 195007 10.1088/0264-9381/25/19/195007
18. Wu H Pac. J. Math. 1967 20 351 10.2140/pjm.1967.20.351
19. Wu H Ill. J. Math. 1964 8 291
