Symmetric teleparallel cosmology with boundary corrections
Abstract
We investigate the geometrodynamical effects of introducing the boundary term in symmetric teleparallel gravity. Specifically, we consider a homogeneous and isotropic universe in , where is the non-metricity scalar, and is the boundary term that relates the non-metricity and Ricci scalars. For the connection in the coincidence gauge, we find that the field equations are of fourth-order, and the fluid components introduced by the boundary are attributed to a scalar field. In the coincidence gauge, the cosmological field equations are equivalent to those of teleparallelism with a boundary term. Nevertheless, for the connection defined in the non-coincidence gauge, the geometrodynamical fluid consists of three scalar fields. We focus on the special case of theory, and we determine a new analytic cosmological solution that can explain the late-time acceleration of the universe and provide a geometric mechanism for the unification of dark energy with dark matter.
I Introduction
General Relativity is a well-tested theory gw1 ; however, it faces challenges on cosmological scales ten1 ; ten1a ; ten2 . Recent cosmological observations rr1 ; Teg ; Kowal ; Komatsu ; suzuki11 suggest that the universe is currently experiencing an acceleration phase driven by an exotic matter source known as dark energy jo1 ; jo . On the other hand, to address the homogeneity and flatness problems, it has been proposed that the universe underwent an inflationary epoch in its early stages driven by the inflaton field inf1 .
The nature and physical properties of the inflaton and dark energy remain subjects of debate among cosmologists. Cosmologists have proposed different models to explain cosmic acceleration, which can be broadly categorized into two groups. In the first category, a matter source with negative pressure is introduced into the field equations of General Relativity ratra ; de1 ; de2 ; de3 ; de4 ; de5 . In the second category of models, cosmologists focus on determining a new gravitational theory by modifying the gravitational Action Integral de01 ; de02 ; de03 .
The consideration of quantum-gravitational effects in the one-loop approximation in gravity df1 ; df2 leads to the introduction of the term in the gravitational Lagrangian. This quadratic theory of gravity qua1 ; qua2 ; qua3 has been used as a mechanism to explain inflation star ; barc ; stt1 . Indeed, when the term dominates, it leads to a de Sitter expansion bar1 . The success of the term in describing acceleration has given rise to a family of theories known as -gravity Buda . For a comprehensive review, we recommend referring to Sotiriou .
The symmetric teleparallel -theory of gravity has garnered the attention of cosmologists lav2 ; lav3 . -theory is defined within the framework of symmetric teleparallel theory nester ; nester1 , where the physical space is described by the metric tensor. However, the fundamental connection in this theory is not the Levi-Civita connection but a symmetric and flat connection, leading to the definition of a nonzero non-metricity scalar . When is linear, the gravitational theory becomes equivalent to General Relativity, known as Symmetric Teleparallel General Relativity (STGR) nester . In STGR, the Ricci scalar and the non-metricity scalar differ only by a boundary term , which means that the variation of the boundary term is neglected, resulting in equivalent field equations. There is a plethora of studies in the literature on -theory; for example, see fq1 ; fq2 ; fq3 ; fq4 ; fq5 ; fq6 ; fq7 ; fq8 ; fq9 ; fq10 ; fq11 , and references therein.
In the process of constructing a gravitational theory, the boundary term has been utilized to modify the gravitational Action Integral. Inspired by teleparallelism bah , the symmetric teleparallel theory with a boundary term, expressed as -gravity, was introduced in ftc1 ; ftc2 . It was discovered that the introduction of the boundary term further modifies the field equations of the gravitational model. In this study, we aim to investigate the effects of the boundary term on cosmological dynamics.
Because the definition of the symmetric and flat connection is not unique, there can be multiple different connections that describe the same gravitational model. The corresponding non-metricity constants differ by a boundary term, which means that they lead to the same equations in STGR. However, this is not true in the case of a nonlinear Lagrangian. As a result, the definition of the connection affects the boundary term, implying that -gravity depends on the definition of the connection.
Recently, in min , it was proven that -theory of gravity admits a minisuperspace description, and the theory is of fourth- or sixth-order, depending on the connection used. The higher-order derivatives can be attributed to a scalar field, and the gravitational Lagrangian can be written as that of a higher-dimensional second-order dynamical system. This approach is applied in this study in the context of -gravity for a homogeneous and isotropic geometry. The structure of the paper is as follows.
In Section II, we present the definition of symmetric teleparallel gravity and its generalization, the -theory, where represents the non-metricity scalar. The extension of symmetric teleparallel theory with a boundary term is discussed in Section III. For the -theory, we discuss when it is equivalent to teleparallel theory and when it approaches the limit of -gravity. In Section IV, we present the gravitational field equations for a homogeneous and isotropic universe. Specifically, we employ Lagrange multipliers to demonstrate that the field equations in -theory can be described within a minisuperspace framework. It’s important to note that the definition of the symmetric and flat connection in a FLRW universe is not unique. Therefore, we provide the field equations for all four different families of connections. One novel aspect of the Lagrange multiplier approach is the introduction of scalar fields to account for the higher-order derivatives in the gravitational model. Consequently, -cosmology is a fourth-order gravitational theory, equivalent to teleparallel -cosmology, when the symmetric and flat connection is defined in the coincidence gauge. However, when the connection is defined in the non-coincidence gauge, the field equations become eighth-order and are described by three scalar fields.
Furthermore, in Section V, we focus on the particular case of cosmology, where the gravitational Action Integral is modified by nonlinear terms of the boundary. For this specific gravitational model, the order of the gravitational field equations is reduced by two in the non-coincidence gauge. Consequently, the geometric fluid is described by two scalar fields. In Section VI, we present a new analytic solution for symmetric teleparallel cosmology with a boundary term in the non-coincidence gauge. This new solution is capable of describing the CDM universe at the present time and possesses a de Sitter attractor. Finally, in Section VII, we summarize our findings and draw our conclusions.
II Symmetric teleparallel gravity
We examine a gravitational model described by the four-dimensional metric tensor and the covariant derivative which is defined using the generic connection , such that the autoparallels are defined as auto
Connection determines the nature of the geometry.
For the general connection we can define the the Riemann tensor
(1) |
the torsion tensor
(2) |
and the non-metricity tensor Eisenhart
In General Relativity, is recognized as the Levi-Civita connection, denoted as . Consequently, in this framework, and . Therefore, the primary scalar in General Relativity is the Ricci scalar .
On the other hand, in the Teleparallel Equivalent of General Relativity (TEGR) de03 , the connection is replaced by the antisymmetric Weitzenb”ock connection, resulting in and . In this context, the torsion scalar takes on the role of the fundamental geometric object in teleparallel gravity.
In the theory under consideration, which is STGR, possesses the property of being both flat and torsionless. This implies that and . Additionally, it inherits the symmetries of the metric tensor . Thus, the nonmetricity scalar defined as nester
(3) |
is the fundamental geometric quantity of gravity.
Tensor is defined as nester
(4) |
which is written with the help of the traces111Parenthesis in the indices denote symmetrization, that is, ; and is the Kroncker delta. and .
The Ricci scalar correspond to the Levi-Civita connection of the metric tensor , and the nonmetricity scalar for a symmetric and flat connection differ by a boundary term defined as .
The gravitational Action Integral of STGR reads nester
(5) |
from where it follows that STGR is dynamically equivalent to GR.
However, when nonlinear components of the non-metricity scalar are introduced in the gravitational Action, as in -gravity, this equivalence is lost. Nevertheless, the corresponding gravitational theory does not possess any dynamical equivalence with General Relativity or its generalization, -gravity.
The Action Integral in symmetric teleparallel -gravity is defined lav2 ; lav3
(6) |
The resulting field equations are
(7) |
(8) |
Equations (7) represent the modified Einstein field equations in -gravity, while equation (8) corresponds to the equation of motion for the connection. . When equation (8) holds true for a specific connection at all times, that connection is referred to as the “coincidence gauge”. However, if the equation (8) is not always satisfied for a particular connection, then that connection is defined within the so-called “non-coincidence gauge”as discussed in lav3 .
III Symmetric teleparallel boundary gravity
Recently, a generalization of the theory was introduced ftc1 ; ftc2 by incorporating a boundary term into the gravitational Action Integral. More precisely, this extended framework includes the gravitational Action Integral as follows
(9) |
where .
In the notation of ftc2 , the boundary term is noted as . This gravitational theory is equivalent to General Relativity, when is a linear function, that is, ; the limit of is recovered when ; while the fourth-order theory of gravity is recovered when .
The gravitational field equations which correspond to the Action Integral (9) are
(10) |
has been inspired by the teleparallel boundary gravity bah , where now is the boundary term relates the Ricciscalar and the torsion scalar for the Weitzenböck connection Weitzenb23 . In a similar way, recovers -gravity when bah , while the limit of GR is recovered when and -theory follows when . There are many astrophysical and cosmological applications of -theory in the literature ftb1 ; ftb2 ; ftb3 ; ftb4 ; ftb5 ; ftb6 . From these results it is clear that the boundary plays an important role in geometric description of dark energy. Thus, the generalization if symmetric teleparallel theory seems natural.
theory of gravity is a fourth-order theory, similar to -theory, while when the resulting field equations are of second-order ftb7 . The order of symmetric teleparallel -theory depends on the connection which is used for the definition of the nonmetricity scalar . For the connection defined in the coincidence gauge is a second-order theory, while for a connection in the non-coincidence gauge -theory is a sixth-order theory of gravity.
To comprehend the degrees of freedom within -theory, we turn our attention to a cosmological model representing an isotropic and homogeneous universe. Employing the Lagrange multipliers method, we introduce scalar fields that account for the dynamical degrees of freedom within the -theory. Consequently, our analysis reveals that the theory introduces either one or three scalar fields.
IV Isotropic and Homogeneous Universe
The isotropic and homogeneous universe is described by the FLRW line element
(11) |
in which is the lapse function and is the scale factor denotes the radius of the universe. Hence, , where is the Hubble function. Parameter is the spatial curvature, for , the universe is spatially flat, corresponds to a closed FLRW geometry and describes an open universe.
FLRW spacetime admits a sixth-dimensional Killing algebra consisted by the vector fields
(12) |
(13) |
IV.1 Non-zero spatial curvature
For the FLRW with , there exist a unique connection defined in the non-coincidence gauge with non-zero components Heis2 ; Zhao
(14) |
and
For this connection, the calculation of the non-metricity scalar is produced
(15) |
From the Levi-Civita connection of spacetime (11) we calculate the Ricciscalar
(16) |
Consequently, the boundary termftc2 is
(17) |
In order to derive the gravitational field equations we apply the mathematical manipulation introduced in min and we introduce the scalar field such that .
We introduce in (9) the Lagrange multipliers and such that
(18) |
Variation with respect to the non-metricity scalar and the boundary term , gives and .
Thus, by replacing in (18) it follows
(19) |
Hence, integration by parts gives
Then
It follows
The minisuperspace Lagrangian is
(20) |
or equivalently
(21) |
in which , and .
In the four-dimensional space , we calculate . This implies that the field equations are of eighth-order and are described by the three scalar fields . The Lagrangian function (21) represents a singular dynamical system in which variation with respect to the lapse function yields the modified Friedmann equation.
Moreover, variation with respect to the dynamical variables gives four second-order differential equations.
For a constant lapse function, i.e. , the cosmological field equation are
(22) |
(23) |
(24) |
(25) |
(26) |
IV.2 Spatially flat case
For the spatially flat case, , it has been found that there exist three families of connections. The common non-zero coefficients of the three connections are
while the additional components for connections and are Heis2 ; Zhao
and
The non-metricity scalars for each connection are calculated
(32) |
(33) |
and
(34) |
Connection is the one defined in the coincidence gauge, while connections and are defined in the non-coincidence gauge. We observe that connection in the limit reduces to that of , that is, and .
We proceed with the derivation of the minisuperspace Lagrangian and the field equations for each family of connections in -theory of gravity.
IV.2.1 Connection
For connection , and the Ricciscalar (16) for the spatially flat FLRW geometry we derive the boundary term
(35) |
For the coincidence gauge, scalars and have the same functional form with the torsion scalar and the boundary . As a result gravity for the connection is equivalent with the teleparallel -gravity.
The Lagrangian of the field equations is
(36) |
or equivalently in scalar field description
(37) |
where similarly as before , and .
For the Lagrangian (37) in the three-dimensional space we derive . Hence, the field equations are of fourth-order.
We select the constant lapse function , and we derive the field equations
(38) |
(39) |
(40) |
(41) |
Modified Friedmann’s equations are written in the equivalent form
with energy density and pressure for the effective fluid
(42) | ||||
(43) |
and .
IV.2.2 Connection
For the non-coincidence connection it follows the boundary term
(44) |
Hence, by introducing Lagrange multipliers as we did for the generic connection , we determine the point-like Lagrangian
(45) |
or equivalently
(46) |
in which , , and .
The field equations are of eight-order described by three scalar fields. For , the equations of motions for the scale factor and the three scalar fields are
(47) |
(48) |
(49) |
(50) |
Hence. the effective geometric fluid has the following energy density and pressure components
(51) | ||||
(52) |
such that
and .
IV.2.3 Connection
We set in (21), thus, the minisuperspace Lagrangian function is
(53) |
in which , and , and the field equations are
(54) |
(55) |
(56) |
(57) |
(58) |
V -Cosmology
Based on the above analysis, we observe that -theory introduces a varying parameter , . Of special interest are the models, in which . and is always a constant. In this theory the nonlinear terms of the Action Integral follow correspond to boundary corrections. This approach has been previously explored in teleparallel -gravity, yielding numerous interesting results. Specifically, has the potential to explain both the late and early-time acceleration phases of the universe ftb5 .
In the theory, the number of field equations is reduced by one, as is not a dynamic parameter but a constant, i.e., . Below, we provide the sets of field equations in -theory for the four different families of connections.
V.1 Connection
For the connection defined in the coincidence gauge, the minisuperspace Lagrangian is
(59) |
where and .
Thus, for , the field equations are
(60) | ||||
(61) |
in which
(62) | ||||
(63) |
and the scalar field satisfies the equation of motion (41). Because the theory is equivalent to the teleparallel model, the results of the latter theory are valid for the the symmetric teleparallel theory with boundary term.
V.2 Connection
For connection defined in the non-coincidence gauge, the minisuperspace Lagrangian reads
(64) |
in which , and .
For , modified Friedmann’s equations are
(65) | ||||
(66) |
with
(67) | ||||
(68) |
where the scalar fields and satisfy the field equations
(69) | ||||
(70) |
V.3 Connection
For the connection the minisuperspace Lagrangian becomes
(71) |
with ,, and .
Furthermore, for , the gravitational field equations are
(72) | ||||
(73) |
where
(74) | ||||
(75) |
(76) |
where the scalar fields and satisfy the equations of motion
(77) |
(78) |
V.4 Connection
Finally, for the fourth connection where curvature is nonzero, the minisuperspace Lagrangian is written
(79) |
in which ,, and .
The gravitational field equations in the presence of curvature are
(80) | ||||
(81) |
The energy density and pressure components of the geometric fluid are defined as
(82) | ||||
(83) |
For the scalar fields we derive the equations of motion
(84) |
(85) |
From the above results we remark that in gravity the field equations are of fourth-order in the coincidence gauge and of sixth-order in the non-coincidence gauge.
VI New solution in the non-coincidence gauge
We focus on the field equations (65)-(70) for the connection in the case of theory. From equation (70) we construct the conservation law
(86) |
For the exponential potential , i.e. ; we are able to write the second conservation law
(87) |
In order to determine the latter conservation law we applied the method of variational symmetries which has been widely used in modified theories of gravity, for more details we refer the reader to nos1 ; nos2 and references therein.
With the use of the two conservation laws and of the constraint equation (65), the second Friedmann equation reads
(88) |
For , we are able to determine the closed form solution
(89) |
Consequently, the Hubble function is derived
(90) |
This analytic solution describes a universe with a cosmological constant, dark matter and a stiff fluid. Indeed when is neglected, i.e. , the limit of CDM universe is recovered. Furthermore, we calculate the deceleration parameter
from where it follows that for large values of , , that is, the de Sitter universe is recovered. Acceleration point is occurred when .
In order to solve equation (88) we apply the Lie symmetry analysis nos1 . We find that equation (88) is invariant under the action of the elements of a two-dimensional Lie algebra consisted by the vector fields and . The application of the Lie invariants of indicate the existence of the exact solution with constraint equation .
On the other hand, the application of the Lie symmetry vector provides the reduced equation
(91) |
This is an Abel type equation.
In Figs. 1 and 2 we present the qualitative evolution of the deceleration parameter and of the function , as it is given after the numerical simulation of equation (88). We observe that the de Sitter universe is a future attractor for the cosmological model, and in the limit of the de Sitter universe connection takes the form of .


VII Conclusions
We conducted a study on FLRW cosmology within the framework of symmetric teleparallel theory, considering boundary corrections in the gravitational Lagrangian. In this context, and for the four different families of connections, we determined the minisuperspace description of the field equations. To achieve this, we introduced Lagrange multipliers, enabling us to express the higher-order derivatives of the field equations in terms of scalar fields. As a result, we were able to recast the cosmological field equations into the equivalent form of multi-scalar field cosmology. In the case of connections defined in the non-coincidence gauge, -gravity is characterized by three scalar fields. However, in the limiting case of the model, the field equations are described by two scalar fields. Conversely, for the connection defined in the coincidence gauge, there exists only one scalar field, and the field equations are of fourth-order. It’s worth noting that -gravity introduces two scalar fields when the connection is defined in the non-coincidence gauge.
This scalar field description and the derivation of the minisuperspace representation are crucial for further investigations into the dynamic evolution of physical variables within the theory. Moreover, the minisuperspace Lagrangian can be employed to establish the Hamiltonian formalism of the model and derive the Wheeler-DeWitt equation of quantum cosmology.
To illustrate the practical application of the minisuperspace description, we employed the method of variational symmetries and successfully determined an integrable cosmological model. We were able to express the analytic solution in terms of the Abel equation. This particular cosmological model not only accounts for cosmic acceleration but also includes a dark matter component in the Hubble function.
These results suggest that -theory holds promise as a viable cosmological framework. However, one notable implication is the significant increase in degrees of freedom introduced by this theory. Therefore, the new scalar fields must be capable of describing a wide range of cosmological phenomena. In future work, we plan to investigate whether boundary correction terms in symmetric teleparallel theory can resolve cosmological tensions and whether the theory can provide explanations for eras in the cosmological history beyond late-time acceleration.
Data Availability Statements: Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.
Acknowledgements.
The author thanks the support of Vicerrectoría de Investigación y Desarrollo Tecnológico (Vridt) at Universidad Católica del Norte through Núcleo de Investigación Geometría Diferencial y Aplicaciones, Resolución Vridt No - 098/2022.References
- (1) B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration), Phys. Rev. Lett. 116, 061102 (2016)
- (2) K. Akiyama et al. (The Event Horizon Telescope Collaboration), ApJL 875, L1 (2019)
- (3) E. Abdalla et. al. JHEAp 34, 49 (2022)
- (4) E. Di Valentino etl. al. Astropart. Phys. 131, 102605 (2021)
- (5) E. Di Valentino etl. al. Astropart. Phys. 131, 102604 (2021)
- (6) A.G. Riess et al., Astron. J. 116, 1009 (1998)
- (7) M. Tegmark et al., Astrophys. J. 606, 702 (2004)
- (8) M. Kowalski et al., Astrophys. J. 686, 749 (2008)
- (9) E. Komatsu et al., Astrophys. J. Suppl. Ser. 180, 330 (2009)
- (10) N. Suzuki et al., Astrophys. J. 746, 85 (2012)
- (11) P. Brax, Rep. Prog. Phys. 81, 016902 (2018)
- (12) J. Yoo and Y. Watanabe, Int. J. Mod. Phys. D 21, 1230002 (2012)
- (13) A.H. Guth, Phys. Rev. D 23, 347 (1981)
- (14) P. Ratra and L. Peebles, Phys. Rev. D 37, 3406 (1988)
- (15) V. Gorini, A. Kamenshchik and U. Moschella, Phys. Rev. D 67, 063509 (2003)
- (16) J.S. Bagla, H.K. Jassal and T. Padmanabhan, Phys.Rev. D 67, 063504 (2003)
- (17) M.C. Bento, O. Bertolami and A.A. Sen, Phys. Rev. D 70, 083519 (2004)
- (18) J.D. Barrow and A. Paliathanasis, Phys. Rev. D 94, 083518 (2016)
- (19) C. Deffayet, X. Gao, D.A. Steer and G. Zahariade, Phys. Rev. D 84, 064039 (2011)
- (20) T. Clifton, P.G. Ferreira, A. Padilla and C. Skordis, Phys. Rept. 513, 1 (2012)
- (21) S. Nojiri, S.D. Odintsov and V.K. Oikonomou, Phys.Rept. 692, 1 (2017)
- (22) R. Ferraro and F. Fiorini, Phys. Rev. D 75, 084031 (2007)
- (23) H. Nariai, Prog. Theor. Phys. 46, 433 (1971)
- (24) F.C. Michel, Ann. Phys. 76, 281 (1973)
- (25) H. Nariai and K. Tomita, Prog. Theor. Phys. 46, 776 (1971)
- (26) G.V. Bicknell, J. Phys. A.: Math. Nucl. Gen. 7, 1061 (1974)
- (27) J.D. Barrow, Nucl. Phys. B 296, 679 (1988)
- (28) A.A. Starobinsky, Phys. Lett. B 91, 99 (1980)
- (29) J.D. Barrow and S. Cotsakis, Phys. Lett. B 214, 515 (1988)
- (30) V.R. Ivanov, S.V. Ketov, E.O. Pozdeeva and S.Yu. Vernov, JCAP 03, 058 (2022)
- (31) J.D. Barrow and A.C. Ottewill, J. Phys. A: Math. Gen. 16, 2757 (1983)
- (32) H.A. Buchdahl, Mon. Not. Roy. Astron. Soc. 150, 1 (1970)
- (33) T.P. Sotiriou and V. Faraoni Rev. Mod. Phys. 82 451 (2010)
- (34) L. Heisenberg, Phys. Reports 796, 1 (2019)
- (35) J.B. Jimenez, L. Heisenberg, T. Koivisto and S. Pekar, Phys. Rev. D 101, 103507 (2020)
- (36) J.M. Nester and H-J Yo, Chinese Journal of Physics 37, 113 (1999)
- (37) M. Hohmann, Phys. Rev. D 104, 124077 (2021)
- (38) R. Solanki, A. De, P.K. Sahoo, Phys. Dark Univ. 36, 100996 (2022)
- (39) A. Lymperis, JCAP 11, 018 (2022)
- (40) L. Atayde and N. Frusciante, Phys. Rev. D 104, 064052 (2021)
- (41) N. Dimakis, A. Paliathanasis, M. Roumeliotis and T. Christodoulakis, Phys. Rev. D 106, 043509 (2022)
- (42) B.J. Barros, T. Barreiro, T. Koivisto and N.J. Nunes, Phys. Dark Univ. 30, 100616 (2020)
- (43) O. Sokoliuk, S. Arora, S. Praharaj, A. Baransky and P.K. Sahoo, Mon. Not. Roy. Astron. Soc. 522, 252 (2023)
- (44) A. Paliathanasis, Phys. Dark Universe 41, 101255 (2023)
- (45) S.K. Maurya, K.N. Singh, M. Govender, G. Mustafa and S. Ray, The effect of gravitational decoupling on constraining the mass and radius for the secondary component of GW190814 and other self-bound strange stars in f(Q)-gravity theory, [arXiv:2309.10130]
- (46) W. Wang, H. Chen and T. Katsuragawa, Phys. Rev. D 105, 024060 (2022)
- (47) F.K. Anagnostopoulos, S. Basilakos and E.N. Saridakis, Phys. Lett. B 822, 136634 (2021)
- (48) F.K. Anagnostopoulos, V. Gakis, E.N. Saridakis, and S. Basilakos, EPJC 83, 58 (2023)
- (49) S. Bahamonde, C.G. Bohmer and M. Wright, Phys. Rev. D 10, 104042 (2015)
- (50) S. Capozziello, V. De Falco and C. Ferrara, The role of the boundary term in f(Q,B) symmetric teleparallel gravity, (2023) [arXiv:2307.13280]
- (51) A. De, T.-H. Loo and E.N. Saridakis, Non-metricity with bounday terms: f(Q,C) gravity and cosmology, (2023) [arXiv:2308.00652]
- (52) A. Paliathanais, N. Dimakis and T. Christodoulakis, Minisuperspace description of f(Q)-cosmology, (2023) [arXiv:2308.15207]
- (53) Y.N. Obuklov and D. Puetzeld, Phys. Rev. D 104, 044031 (2021)
- (54) L. P. Eisenhart, “Non-Riemannian Geometry”, American Mathematical Society, Colloquium Publications Vol. VIII, New York, (1927)
- (55) R. Weitzenböck, Invarianten Theorie, Nordhoff, Groningen (1923)
- (56) M. Wright, Phys. Rev. D 93, 103002 (2016)
- (57) M. Caruana, G. Farrugia and J.L. Said, EPJC 80, 640 (2020)
- (58) S. Bahamonde, A. Golovnev, M.-J. Guzmán, J.L. Said and C. Pfeifer, JCAP 01, 037 (2022)
- (59) A. R. P. Moreira, J. E. G. Silva and C.A.S. EPJC 81, 298 (2021)
- (60) A. Paliathanasis, JCAP 1708, 027 (2017)
- (61) C. Escamilla-Rivera, J.L. Said, EPJC 80, 677 (2020)
- (62) A. Paliathanasis and G. Leon, EPJC 81, 718 (2021)
- (63) F. D’ Ambrosio, L. Heisenberg and S. Kuhn, Revisiting cosmologies in teleparallelism, Class. Quantum Grav. 39 025013 (2022)
- (64) D. Zhao, Covariant formulation of f(Q) theory, Eur. Phys. J. C 82, 303 (2022)
- (65) M. Tsamparlis and A. Paliathanasis, Symmetry 10, 233 (2018)
- (66) P.A. Terzis, N. Dimakis, T. Christodoulakis, A. Paliathanasis and M. Tsamparlis, J. Geom. Phys. 101, 52 (2016)