Dynamics in interacting scalar-torsion theory

Andronikos Paliathanasis1, 2, ∗

1Institute of Systems Science, Durban University of Technology, Durban 4000, South Africa 2Instituto de Ciencias F´ısicas y Matem´aticas, Universidad Austral de Chile, Valdivia 5090000, Chile (Dated: July 19, 2021)

In a spatially flat Friedmann–Lemaˆıtre–Robertson–Walker background space we consider a scalar-torsion gravitational model which has similar properties with the dilaton theory. This teleparallel model is invariant under a discrete transformation similar to the Gasperini-Veneziano duality transformation. Moreover, in the gravi- tational Action integral we introduce the Lagrangian function of a pressureless fluid source which is coupled to the teleparallel dilaton field. This specific gravitational theory with interaction in the dark sector of the universe is investigated by using methods of the dynamical system analysis. We calculate that the theory provides various areas of special interest for the evolution of the cosmological history. Infla- tionary scaling solutions and the de Sitter universe is recovered. Furthermore, we calculate that there exist an attractor which provides a stable solution where the two fluid components, the scalar field and the pressureless matter, contribute in the cosmological fluid. This solution is of special interest because it can describe the present epoch. Finally, the qualitative evolution of the cosmographic parameters is discussed.

arXiv:2107.05880v2 [gr-qc] 16 Jul 2021 PACS numbers: 98.80.-k, 95.35.+d, 95.36.+x

Keywords: Teleparallel; scalar field; dilaton field; scalar-torsion; interaction

∗Electronic address: [email protected] 2

1. INTRODUCTION

The detailed analysis of the recent cosmological data indicates that General Relativity may need to be modified in order to describe the observations, for a recent review we refer the reader to [1]. The late-time cosmic acceleration has been attributed to a fluid, the so- called dark energy, with has negative pressure and anti-gravity effects [2]. In order to explain the anti-gravitational effects, cosmologists have proposed the modification of the Einstein- Hiblert Action by using geometric invariants [3]. In this direction new geometrodynamical terms are introduced in the field equations which can drive the dynamics in order to explain, with a geometric approach, the observational phenomena [4, 5]. Teleparallelism [6, 7] includes a class of modify theories of gravity which have been widely studied the last years [8–15]. The fundamental invariant of teleparallelism is the torsion scalar of the antisymmetric connection which play the role important role, instead of the Levi-Civita connection in General Relativity. In previous studies it has been discussed that teleparallel gravity my violates the Lorentz symmetry [16], while Lorentz violation has not been observed yet, it is common in various subjects of quantum gravity [17]. However, re- cent studies has shown that Lorentz symmetry can be preserved in teleparallelism, see for instance [18, 19]. Specifically, in [18] it was found that the introduction of a scalar field in the gravitational Action Integral of teleparallelism preserves the Lorentz symmetry. For a recent review on teleparallelism we refer the reader in [20]. In the literature, teleparallel cosmology has been widely studied. For instance, in f (T ) teleparallel theory the cosmolog- ical perturbations were investigated in [21, 22], while in [22] it was found that f (T ) theory can mimicking dynamical dark energy models. The mechanism which describe the Higgs inflation era in scalar torsion theory was studied in [23]. An extension of the scalar torsion theories with the introduction of the boundary term was introduced in [24]. For other recent studies on teleparallelism we refer the reader in [25–29] and references therein. In this study we focus on the scalar-torsion or teleparallel dark energy models [30–32] which can be seen as the analogue of the scalar-tensor theories in teleparallelism. In scalar- torsion theory an scalar field is introduced in the gravitational Action Integral which is nonminimally coupled to the fundamental scalar of teleparallelism, the torsion scalar. There are various studies in the literature which indicate that the scalar-torsion theories can explain the recent cosmological observations [34, 35]. In the following we consider the existence of 3 a matter source with zero pressure in the gravitation Action which interacts with the scalar field [36–39]. In our analysis the interaction between the scalar field and the pressureless fluid are inspired by the interaction provided in the Weyl integrable theory [40, 41]. The plan of the paper is as follows. In Section 2 we introduce the model of our consideration which is the teleparallel dilaton model coupled to a pressureless fluid source. This model belongs to the family of scalar- torsion theories, the field equations are derived. Section 3 includes the new material of this study. We performed a detailed study of the asymptotic dynamics for the gravitational field equations for the model of our consideration. We determine the stationary points and we investigate their stability as we discuss the physical properties of the exact solutions described by the stationary points. This analysis provide important information about the cosmological viability of the proposed theory. It is clear that our model can explain the mayor eras of the cosmological evolution and it can be used as a dark energy candidate. Furthermore, in Section 4 we discuss the evolution of the cosmographic parameters as they are provided by our model. Finally, in Section 5 we summarize our results and we draw our conclusions.

2. TELEPARALLEL DILATON MODEL

The gravitational model of our consideration is an extension of the teleparallel dilaton theory known as scalar-torsion theory where the gravitational Action Integral is defined Z 1 4 h − φ  ω ;µ i S = d xe e 2 T + φ φ + V (φ) + L , (1) 16πG 2 ;µ m in which T is the torsion scalar of the antisymmetric curvatureless connection, φ xk
is a scalar field with potential function V (φ) , ω is a nonzero constant, Lm is the Lagrangian function for the additional matter source and e is the determinant of the vierbein fields. Action Integral (1) belongs to the family of gravitational models known as teleparallel dark energy models, or scalar torsion models [30, 31, 42, 43]. Scalar-torsion models can be seen as the analogue of scalar-tensor models in teleparallelism, in which instead of the Ricci scalar the torsion scalar T is used for the definition of the Action Integral. The gravitational field equations are of second-order, however, under conformal transformations the scalar- torsion theories are now equivalent with the quintessence model [42]. Under a conformal 4 transformation the scalar-torsion Action Integral is equivalent with that of a modified higher- order theory known as f (T,B) where B is the boundary term differs the torsion T and the Ricci scalar [42]. That makes the scalar-torsion and scalar-tensor theories to be totally different which means that there is not any conformal transformation which may connect the solutions of the two theories. In this study, we assume that the matter source and the − φ scalar field are interacting, that is, from (1) the mixed term e 2 Lm exists. In the case of a spatially flat Friedmann–Lemaˆıtre–Robertson–Walker metric (FLRW) background space ds2 = −dt2 + a2(t)(dx2 + dy2 + dz2). (2) the gravitational field equations for the (1) in the case of vacuum are invariant under a discrete transformation with origin the O (d, d) symmetry [44]. The resulting discrete trans- formation is a generalization of the Gasperini-Veneziano scale-factor duality transformation for the dilaton field in scalar tensor theory [45]. The existence of this discrete transformation it is important for the study of the pre-big bang period of the universe as it is described by string cosmology. However, when parameter ω is small, the Gasperini-Veneziano transfor- mation is recovered. However, in general the pre-big bang period for the teleparallel dilaton theory differs from that of sting cosmology by a term provided by the nonzero constant ω. In this study we assume the scalar field φ interact with the matter source. For the latter, we assume that of a pressureless fluid source, dust fluid, with energy density ρm, known as dark matter. Models with interaction in the dark sector of the universe has drawn the attention of cosmologists the last decade [46–50]. Indeed, there are various theoretical studies which shows that such models are viable, while from the analysis of the cosmological data it seems that the interacting models are supported by the observations [51–54]. The Action Integral (1) can be seen as the teleparallel extension of some scalar tensor models coupled to dark matter [40, 41]. 2 a˙ For the background space (2) it follows T = 6H , H = a , hence, from the Action integral (1) and for a dust fluid for the matter source we derive the modified gravitational field equations − φ  2 ω 2   e 2 6H − φ˙ + V (φ) − ρ = 0 , (3) 2 m   1 ω  2H˙ + 3H2 + φ˙2 − V (φ) − Hφ˙ = 0 , (4) 2 2 5

while for the matter source and the scalar field it follows

 ˙ ρ˙m + 3H − φ ρm = 0 , (5)

 ˙2 ¨ ˙ ω φ − 2φ − 6Hφ + 2 (V − V,φ) = 0. (6)

In the case where ρm = 0 and V (φ) = V0, the discrete transformation which keeps invariant the field equations is

p1 p2φ a → a e , φ → p3 ln a + p4φ (7)

1+3ω ω 12 1+3ω 1 in which p1 = − 1−3ω , p2 = 1−3ω , p3 = − 1−3ω and p4 = 1−3ω when ω 6= 3 . However, in the presence of the matter source the field equations do not remain invariant under the action of the discrete transformation. Moreover, when ω is near to zero the discrete transformation (7) becomes the Gasperini-Veneziano duality transformation [45]. The discrete transformation (7) follows from the presence of the O (d, d) symmetry for the Action Integral (1). In addition, the existence of the O (d, d) symmetry indicates the presence of a conservation law for the classical field equations which can be used in order to integrate and write the analytic solution in closed-form expression. Finally, the O (d, d) symmetry is preserved and in the quantization process of the theory. Hence, a quantum operator related to the O (d, d) symmetry is determined which helps us to solve the Wheeler-DeWitt equation [44]. Additionally, from equation (3) we define the energy density for the scalar field   ω ˙2 ρφ = 2 φ + V (φ) , thus in order ρφ ≥ 0 and do not have ghost terms we assume

ω > 0. The pressure component pφ of the scalar field from equation (4) it is defined as   1 ω ˙2 ˙ pφ = 2 2 φ − V (φ) − Hφ. In the following Section we study the general evolution for the cosmological history as it is provided by the dynamical system (3)-(6).

3. ASYMPTOTIC DYNAMICS

We define the new dimensionless variables r ω φ˙ V ρ V x = , y2 = , Ω = m , λ = ,φ (8) 12 H 6H2 m 6H2 V where the field equations are written as the following algebraic-differential system

2 2 Ωm = 1 − x − y (9) 6

dx 3 r 3 = x x2 − y2 − 1
+ y2 (1 − λ) , (10) dτ 2 ω ! dy 3 r 3 = y 1 + x2 − y2
− x (1 − λ) , (11) dτ 2 ω √ dλ 2 V,φφV = 2 3xλ (Γ (λ) − 1) , Γ(λ) = 2 , (12) dτ (V,φ) in which the new independent variable is dt = Hdτ. In addition, in the new variables the equation of state parameter for the effective fluid 2 H˙ weff = −1 − 3 H2 reads 2 2 2 weff (x, y, η, λ) = x − y − √ x. (13) 3ω As far as the scalar field potential is concerned, we consider the exponential potential

λ0φ V (φ) = V0e . (14)

For the exponential potential we infer that λ = λ0 is a constant and Γ (λ) = 1 such that the rhs of equation (12) is always zero. Thus, for the exponential potential the dimension of the dynamical system is two. Moreover, in order the matter source to be physically accepted, it follows that Ωm is constraint as, 0 ≤ Ωm ≤ 1, from (9) it follows that the variables (x, y) are constrained similarly as x2 + y2 ≤ 1. Furthermore, we observe that the equations are invariant under the change of variables y → −y, which means that without loss of generality we can work on the branch y ≥ 0.

3.1. Stationary points for the exponential potential

We proceed our analysis by considering the exponential scalar field potential. For ar- bitrary value of the free parameter λ the dynamical system (10)-(11) admits the following stationary points P = P (x, y):

Point A1 with coordinates A1 = (0, 0), in which Ωm (A1) = 1, and weff (A1) = 0. At the stationary point the scalar field is constant, that is φ = φ0. Hence, the point describes a 2 universe dominated by the dust fluid source, where the scale factor is a (t) = a0t 3 . In order to infer on the stability of the stationary points we linearize equations (10)-(11) around the stationary point, and we derive the eigenvalues of the linearized matrix, they are calculated 3 3 e1 (A1) = − 2 , e2 (A1) = 2 . Hence, because only one of the eigenvalues is negative, the stationary point is always a saddle point. 7

TABLE I: Stationary points of the field equations and their physical properties for the exponential potential.

1 Point (x, y) Existence Ωm weff weff < − 3

A1 (0, 0) Always 1 0 No, dust solution

A± (1, 0) Always 0 1 ∓ √2 A+ for ω > 4 2 3ω 2 3  q  (1−λ)2 2λ(λ−1) A ± √1−λ , 1 − ω > 1 (1 − λ)2 0 −1 + Yes under conditions 3 3ω 3ω 3 3ω  √ q  A 3ω , 3ω λ 6= 1 1 − 3ω 1 Yes under conditions 4 2(1−λ) 4(1−λ)2 2(1−λ)2 λ−1

2 2 ω < 3 (1 − λ)

TABLE II: Stationary points of the field equations and their stability properties for the exponential potential.

Point Eigenvalues Stability

3 3 A1 − 2 , 2 Saddle ± q 3 + 1 2 A2 3 ± ω (λ − 1) , 3 A2 saddle λ < 1 and ω < 3 (1 − λ) − 1 2 A2 saddle λ > 1 and ω < 3 (1 − λ) 2 A −3 + 2 (1 − λ)2 , − 3 + (1−λ) Stable for ω > 2 (1 − λ)2 3 ω √ ω 3 2 3 12ω−7(λ−1) A4 − 4 ± 4(λ−1) Attractor

± The stationary points A2 = (±1, 0) describe universes where only the scalar field con- ±
±
tributes in the cosmological solution, that is, Ωm A2 = 0. Because y A2 = 0, it means that the scalar field potential does not contribute in the total cosmological evolution, while only the kinematic part contributes, that is, V (φ) << φ˙2. The effective equation of state parameter is derived w A±
= 1 ∓ √2 . For point A+ it follows that w A+
< 1. eff 2 3ω 2 eff 2 + 4 Furthermore, the asymptotic solution at A2 describes an accelerated universe when ω > 3 . − −
On the other hand, the asymptotic solution at the stationary point A2 gives weff A2 > 1. As far as the stability properties of the points are concerned, we derive the eigenvalues ±
q 3 ±
+ e1 A2 = 3 ± ω (λ − 1) , e2 A2 = 3. Thus, point A2 is a saddle point when λ < 1 1 2 + − and ω < 3 (1 − λ) , otherwise A2 is a source. Point A2 a saddle point when λ > 1 and ω < 1 (1 − λ2). 3   q (1−λ)2 The stationary point A = √1−λ , 1 − is physically accepted when ω > 3 3ω 3ω 8

1 2 3 (1 − λ) . Point A3 describes a scaling solution for an ideal gas with parameter for the 2λ equation of state weff (A3) = −1 + 3ω (λ − 1), and Ωm (A3) = 0. The point describes an accelerated universe when the free parameters are constraint as {λ ≤ 0, ω > λ (λ − 1)}, {0 < λ ≤ 1} and {λ > 1, ω > λ (λ − 1)}. For λ = 0 or λ = 1, the stationary points describe the de Sitter universe i.e. weff (A3) = −1. The eigenvalues of the linearized system are 2 2 (1−λ)2 derived to be e1 (A3) = −3 + ω (1 − λ) , e2 (A3) = −3 + ω . Therefore, point A3 is an 2 2 attractor when ω > 3 (1 − λ) and for arbitrary value for the parameter λ.  √ q  Finally, the stationary points A = 3 , 3ω exist when λ 6= 1. The latter 4 2(1−λ) 4(1−λ)2 1 stationary points describe scaling solutions with weff (A4) = λ−1 , for λ 6= 0, or de Sitter

universes weff (A4) = −1 when λ = 0. In addition, the contribution of the dust fluid source is nonzero, i.e. Ω (A ) = 1− 3ω . The point is physically accepted, when 0 < Ω (A ) < 1, m 4 2(1−λ)2 m 4 that is, λ 6= 1, 0 < ω < 2 (1 − λ)2. The eigenvalues of the linearized system are derived √ 3 √ 2 2 3 12ω−7(λ−1) 3 12ω−7(λ−1) e1 (A4) = − 4 + 4(λ−1) , e2 (A4) = − 4 − 4(λ−1) . We infer that point A4 is always an attractor. The results for the stationary points and their physical properties are summarized in Table I, while in Table II we summarize the stability conditions for the stationary points. In Fig. 1 we present the phase-space portrait for the field equations in the plane (x, y) for

different values of the free parameters in which point A3 or A4 are attractors. In Figs. 2

we present the parametric plot with the evolution of the physical variables Ωm and weff for numerical solutions of the field equations.

The stationary point A4 is of special interest, because it can describe an accelerated universe where the dark matter and the dark energy contributes in the cosmological fluid. This point is always an attractor when it exists. Moreover, from the recent cosmological 1 observations we know that the deceleration parameter q = 2 (1 + 3weff ) for the Λ-cosmology

is qΛ ' −0.6 where Ωm ' 0.27 [55]. Thus, the effective equation of state parameter is Λ weff ' −0.73. Therefore, for (ω, λ) ' (0.91, −0.37) the solution of point A4 provide physical variables equal with that of the Λ-cosmology. 9

4. COSMOGRAPHIC PARAMETERS

The cosmographic approach is a model independent construction way of the cosmological physical variables [56]. Specifically, the scale factor it is written in the expansion form

a(t) 1 2 2 1 3 1 4 5 = 1 + H0(t − t0) − q0(t − t0) + j0(t − t0) + s0(t − t0) + O[(t − t0) ] , (15) a0 2 3! 4!

where H0 is the value of the Hubble function of today, q0 is the deceleration parameter of

today, j0 and s0 are the present value for the jerk and snap parameters. The H, q, j, s are kinematical quantities which are extracted directly from the spacetime. The kinematic quantities are defined as [57, 58] 1 da H(t) = , a dt 1 d2a 1 da−2 q(t) = − , a dt2 a dt 1 d3a 1 da−3 j(t) = , a dt3 a dt 1 d4a 1 da−4 s(t) = , a dt4 a dt or in terms of the Hubble function they can be written in the equivalent form H˙ q = −1 − (16) H2 H¨ j = 3 − 3q − 2 (17) H... H s = + 4j + 3q (q + 4) + 6, (18) H4 Thus, from the evolution of the cosmographic q, j, s, one can understand the expansion of the universe, the rate of acceleration and its derivative. In Figs. 3 and 4 we present the qualitative evolution of the cosmographic parameters for the model of our consideration for initial conditions near to the matter dominated era and for the values of the free parameters (λ, ω) same with that of the numerical solutions of Fig. 2. From the evolution of the cosmographic parameters presented in Figs. 3 and 4 we observe that while the q0 value for the ΛCDM is recovered that it is not true for the jerk and snap parameters which in general have different evolution. From the qualitative evolution we observe that our model can predict values for the cosmographic parameters as they given by the cosmological constraints [59]. 10

5. CONCLUSIONS

In this study we considered a scalar-torsion model known as teleparallel dilaton theory coupled to a pressureless fluid source, which we assume that it describes the dark matter. For this cosmological model the field equations are of second order and we investigated the evolution of the cosmological parameters in a spatially flat FLRW background space by determine the stationary points and study their stability. This kind of dynamical analysis is essential for the study of the general evolution of the dynamical system because it provide us with important information in order to infer about the cosmological viability of theory. For the model of our consideration we wrote the field equations into an equivalent algebraic-differential by using the H−normalization approach. Every stationary point of the dynamical system describes an exact asymptotic solution for the scale factor which corresponds to a specific epoch of the cosmological history. The stability properties of the stationary points are important to investigated because they tell us about the general evolution of the dynamical system.

For our model, and for the exponential scalar field potential, we found that the field equations admit four stationary points which describe four different cosmological epochs. ± Point A1 provides the exact solution of the unstable matter dominated era; points A2 describe universes dominated by the kinetic part of the scalar field where in contrary to the quintessence at this theory we calculate w A±
= 1 ∓ √2 . Hence for ω > 4 point eff 1 3ω 3 + A2 describes the cosmic acceleration. Point A3 is a scaling solution in general weff (A3) = 2λ(λ−1) −1 + 3ω , where for λ (λ − 1) = 0, the asymptotic solution at point A3 is the de Sitter solution. Finally, point A4 describes an asymptotic solution in which the dark matter and the scalar field contributes in the cosmological fluid. This solution is of special interest because it describes the present cosmological era. Moreover, in order to understand the evolution of the physical variables we presented the qualitative evolution of the cosmographic parameters from where we found that the values for some of the cosmographic parameters at the present era can be recovered. From the above results it is clear that while there are similarities of this model with the dilaton model of scalar tensor theory [60], the two theories are different in the general evolution of the cosmological history. In the same conclusion we end by comparing the 11 results of this work that previous studies on the Weyl Integrable theory, where interaction of the scalar field with the dust fluid is introduced in the gravitational Action Integral [61]. Consequently, the present model of study has interesting properties which can explain the cosmic history and deserves further study. In addition, in the absence of the dust fluid we recall that the theory admits a discrete symmetry which can be used to study the cosmic evolution in the pre-big bang era as in string cosmology. However, in contrary to string cosmology and the dilaton field for this teleparallel model the pre-big bang era is not recovered as reflection of the present epoch. In a future study, we plan to use the cosmological observations in order to constraint this specific theory as a dark energy candidate.

[1] E. Di Valentino, O. Mean, S. Pan, L. Visinelli, W. Yang, A. Melchiorri, D.F. Mota, A.G. Riess and J. Silk, In the Realm of the Hubble tension - a Review of Solutions, to appear in Class. Quantum Grav. [doi: 10.1088/1361-6382/ac086d] (2021) [arXiv:2103.01183] [2] J. Yoo and Y. Watanabe, Theoretical Models of Dark Energy, Int. J. Mod. Phys. D 21, 1230002 (2012) [3] T. Clifton, P.G. Ferreira, A. Padilla and C. Skordis, Modified Gravity and Cosmology, Physics Reports 513, 1 (2012) [4] S. Capozziello, Curvature quintessence, Int. J. Mod. Phys. D 11, 483 (2002) [5] S. Nojiri, S.D. Odintsov and V.K. Oikonomou, Physics Reports 693, 1 (2017) [6] A. Einstein 1928, Sitz. Preuss. Akad. Wiss. p. 217; ibid p. 224; A. Unzicker and T. Case, Translation of Einstein’s attempt of a unified field theory with teleparallelism, (2005) [physics/0503046] [7] K. Hayashi and T. Shirafuji, New general relativity, Phys. Rev. D 19, 3524 (1979)

[8] R.C. Nunes, Structure formation in f(T ) gravity and a solution for H0 tension, JCAP 05, 052 (2018) [9] R. Ferraro and F. Fiorini, Modified teleparallel gravity: Inflation without an inflaton, Phys. Rev. D 75, 084031 (2007) [10] A. Paliathanasis, J.D Barrow and P.G.L. Leach, Cosmological solutions of f(T ) gravity, Phys. Rev. D 94, 023525 (2016) 12

[11] L. Iorio, N. Radicella and M.L. Ruggiero, Constraining f(T ) gravity in the Solar System, JCAP 08, 021, (2015) [12] M. Khurshudyan, R. Myrzakulov and As. Khurshudyan, Interacting dark energy models in f(T ) gravity, Mod. Phys. Lett. A 32, 1750097 (2017) [13] A. Paliathanasis, de Sitter and Scaling solutions in a higher-order modified teleparallel theory, JCAP 08, 027 (2017) [14] C. Escamilla-Rivera, J.L. Said, Cosmological viable models in f(T,B) theory as solutions to

the H0 tension, EPJC 80, 677 (2020) [15] M. Hohmann, M. Krˇsˇs´ak,C. Pfeifer and U. Ualikhanova, Propagation of gravitational waves in teleparallel gravity theories, Phys. Rev. D 98, 124004 (2018) [16] G. Otalora and M.J. Reboucas, Violation of causality in f(T ) gravity, EPJC 77, 799 (2017) [17] D. Mattingly, Modern tests of Lorentz invariance, Living Reviews in Relativity 8, 5 (2005) [18] M. Hohmann, L Jarv and U. Ualikhanova, Covariant formulation of scalar-torsion gravity, Phys. Rev. D 97, 104011 (2018) [19] M. Krssak, R.J. van den Hoogen, J.G. Pereira, C.G. B¨ohmerand A.A. Coley, Teleparallel theories of gravity: illuminating a fully invariant approach, Class. Quantum Grav. 36, 183001 (2019) [20] S. Bahamonde, K.F. Dialektopoulos, C. Escamilla-Rivera, G. Farrugia, V. Gakis, M. Hendry, M.Hohmann, J.L. Said, J.Mifsud and E. Di Valentino, Teleparallel Gravity: From Theory to Cosmology, (2021) [arXiv:2106.13793] [21] S.H. Chen, J.B. Dent, S. Dutta and E.N. Saridakis, Cosmological perturbations in f (T) gravity, Phys. Rev. D 83, 023508 (2011) [22] J.B. Dent, S. Dutta and E.N. Saridakis, f (T) gravity mimicking dynamical dark energy. Background and perturbation analysis, JCAP 01, 009 (2011) [23] S. Raatikainen and S. Rasanen, Higgs inflation and teleparallel gravity, JCAP 12, 021 (2019) [24] S. Bahamonte, M. Marciu and P. Rudra, Generalised teleparallel quintom dark energy non- minimally coupled with the scalar torsion and a boundary term, JCAP 04, 056 (2018) [25] P.J. Pompeia, Scalar-multitensor approach to teleparallel modified theories of gravity, Phys. Rev. D 103, 124036 (2021) [26] J. Beltran Jimenez, T.S. Koivisto, Accidental gauge symmetries of Minkowski spacetime in Teleparallel theories, Universe 7, 142 (2021) 13

[27] L. Jarv and J. Lember, Global Portraits of Nonminimal Teleparallel Inflation, Universe 7, 178 (2021) [28] T. Harko, F.S.N. Lobo and E.N. Saridakis, Universe 7, 227 (2021) [29] H. Boumaza, Cosmology of cubic galileon in modified teleparallel gravity, EPJC 81, 114 (2021) [30] H. Wei, Phys. Dynamics of Teleparallel Dark Energy, Lett. B 712, 430 (2012) [31] C.-Q. Geng, C.-C. Lee, E.N. Saridakis and Y.-P. Wu, ”Teleparallel” Dark Energy, Phys. Lett. B 704, 384 (2011) [32] C. Xu, E.N. Saridakis and G. Leon, Phase-Space analysis of Teleparallel Dark Energy, JCAP 07, 005 (2012) [33] V. Faraoni, Cosmology in Scalar-Tensor Gravity, Springer Netherlands, (2004) [doi:10.1007/978-1-4020-1989-0] [34] C.-Q. Geng, C.-C Lei and E.N. Saridakis, JCAP 01, 002 (2012) [35] R. D’Agostino and O. Luongo, Phys. Rev. D 98, 124013 (2018) [36] L. Amendola, Coupled quintessence, Phys. Rev. D 62, 043511 (2000) [37] S. Pan and S. Chakraborty, A cosmographic analysis of holographic dark energy models, Int. J. Mod. Phys. D 23, 1450092 (2014) [38] S. Pan, G. S. Sharov and W. Yang, Field theoretic interpretations of interacting dark energy scenarios and recent observations, Phys. Rev. D 101, 103533 (2020) [39] W. Yang, S. Pan, E. Di Valentino, R.C. Nunes and S. Vagnozzi, JCAP 09, 019 (2018) [40] J.M. Salim and S.L. Sautu, Gravitational theory in Weyl integrable spacetime, Class. Quantum Grav. 13, 353 (1996) [41] J.M. Salim and S. Sautu, Spherically symmetric static solutions in Weyl integrable spacetime, Class. Quantum Grav. 15, 203 (1998) [42] M. Hohmann, Scalar-torsion theories of gravity I: general formalism and conformal transfor- mations, Phys. Rev. D 98, 064002 (2018) [43] M. Hohmann and C. Pfeifer, Scalar-torsion theories of gravity II: L (T, X, Y, φ) theory, Phys. Rev. D 98, 064003 (2018) [44] A. Paliathanasis, O (d, d) symmetry in teleparallel dark energy, Eur. Phys. J. Plus 136, 674 (2021) [arXiv:2105.08527] [45] M. Gasperini and G. Veneziano, Pre-big-bang in string cosmology, Astropart. Phys. 1, 317 (1993) 14

[46] M. Quartin, M. O. Calvao, S. E. Joras, R. R. R. Reis and I. Waga, Dark Interactions and Cosmological Fine-Tuning, JCAP 0805, 007 (2008) [47] L. Amendola, Perturbations in a coupled scalar field cosmology, Mon. Not. Roy. Astron. Soc. 312, 521 (2000) [48] F. Arevalo, A. P. R. Bacalhau and W. Zimdahl, Cosmological dynamics with nonlinear inter- actions, Class. Quant. Grav. 29, 235001 (2012) [49] A. Paliathanasis, S. Pan and W. Yang, Dynamics of nonlinear interacting dark energy models, IJMPD 28, 1950161 (2019) [50] Z. Razaei, Dark Matter-Dark Energy Interaction and the Shape of Cosmic Voids, The Astro- physical Journal. 902, 102 (2020) [51] L. Amendola, G. Camargo Campos and R. Rosenfeld, Consequences of dark matter-dark energy interaction on cosmological parameters derived from SNIa data, Phys. Rev. D 75, 083506 (2007) [52] G.S. Sharov, S. Bhattacharya, S. Pan, R.C. Nunesand S. Chakraborty, A new interacting two fluid model and its consequences, Mon. Not. Roy. Astron. Soc. 466, 3497 (2017) [53] W. Yang, S. Pan and A. Paliathanasis, Cosmological constraints on an exponential interaction in the dark sector, Mon. Not. Roy. Astron. Soc. 482, 1007 (2019) [54] S. Pan, W. Yang and A. Paliathanasis, Non-linear interacting cosmological models after Planck

2018 legacy release and the H0 tension, Mon. Not. Roy. Astron. Soc. 493, 3114 (2020) [55] Planck Collaboration: N. Aghanim et al., Planck 2018 results. VI. Cosmological parameters, A&A 641, A6 (2020) [56] S. Weinberg, Gravitation and cosmology: Principles and applications of the general theory of relativity, (Wiley, New York, 1972) [57] M. Visser, Jerk and the cosmological equation of state, Class. Quant. Grav. 21, 2603-2616 (2004) [58] Y. L. Bolotin, V. A. Cherkaskiy, O. Y. Ivashtenko, M. I. Konchatnyi and L. G. Zazunov, Applied Cosmography: A Pedagogical Review, arXiv:1812.02394 (2018) [59] F.Y. Wang, Z.G. Dai and S. Qi, Probing the cosmographic parameters to distinguish between dark energy and modified gravity models, A&A 507, 53 (2009) [60] O. Hrycyna and M. Szydlowski, Dynamical complexity of the Brans-Dicke cosmology, JCAP 12, 016 (2013) 15

[61] A. Paliathanasis, G. Leon and J.D. Barrow, Eur. Phys. J. C 80, 731 (2020) 16

λ=1/2,ω=1 λ=-1/2,ω=1

1.0 1.0

0.8 0.8

0.6 0.6 y y 0.4 0.4

0.2 0.2

0.0 0.0

-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 x x

λ=-1,ω=1 λ=1,ω=1

1.0 1.0

0.8 0.8

0.6 0.6 y y 0.4 0.4

0.2 0.2

0.0 0.0

-1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 x x

FIG. 1: Two-dimensional phase-space diagrams for the dynamical system (10)-(11) in the

1
plane (x, y) Left figure in the first row is for (λ, ω) = 2 , 1 with attractor point A3. Right 1
figure in the first row is for (λ, ω) = − 2 , 1 , with attractor point A4. For the figures of the second row, left figure is for (λ, ω) = (−1, 1) with attractor point A4, while the right figure is for

(λ, ω) = (1, 1) and attractor the de Sitter point A3. 17

0.0 λ=-1/5 0.0 ω=0.2 λ=-1/2 ω=0.5 -0.2 λ=0 -0.2 ω=0.7 λ=1/5 0.9 -0.4 ω= -0.4 weff -0.6 weff

-0.8 -0.6

-1.0 -0.8

0.0 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 Ωm Ωm

FIG. 2: Parametric plot for the qualitative evolution of the Ωm and of the effective equation of state parameter weff for numerical solution of the dynamical system (9)-(11) with initial conditions near 1 to the matter dominated era Ωm ' 1. In Left Fig., lines are for ω = 1, solid line is for λ = − 5 , 1 1 dashed line is for λ = − 2 , dotted line is for λ = 0 and dashed dot line is for λ = 5 . In the Right Fig. lines are for λ = 0, solid line is for ω = 0.2, dashed line is for ω = 0.5, dotted line is for

ω = 0.7 and dashed dot line is for ω = 0.9. The attractors is the de Sitter point A4. 18

λ=-1/5 ,ω=1 λ=-1/2 ,ω=1 1.0 1.0

0.5 0.5

0.0 0.0

-0.5 -0.5 q q Cosmographic Parameters Cosmographic Parameters j -1.0 j -1.0 s s

ln(a) ln(a)

λ=0 ,ω=1 λ=1/5 ,ω=1 1.0 1

0.5 0

0.0 -1 -0.5

-2 -1.0 q q Cosmographic Parameters Cosmographic Parameters j -1.5 -3 j s s

ln(a) ln(a)

FIG. 3: Qualitative evolution for the cosmographic parameters q, j, and s as provided by the field equations (10)-(11) for different values of the free parameters. The initial condition is a point near to the matter dominated solution A1. Solid line is for decellaration parameter q, dashed line is for the jerk parameter j and dotted line is for the snap parameter s. 19

λ=-0.37 ,ω=0.2 λ=-0.37 ,ω=0.5 1.0 1.0

0.5 0.5

0.0 0.0

-0.5 -0.5

-1.0 -1.0

-1.5 q -1.5 q Cosmographic Parameters Cosmographic Parameters j j -2.0 -2.0 s s

ln(a) ln(a)

λ=-0.37 ,ω=0.7 λ=-0.37 ,ω=0.9 1.0 1.0

0.5 0.5

0.0 0.0

-0.5 -0.5

-1.0 -1.0

-1.5 q -1.5 q Cosmographic Parameters Cosmographic Parameters j j -2.0 -2.0 s s

ln(a) ln(a)

FIG. 4: Qualitative evolution for the cosmographic parameters q, j, and s as provided by the field equations (10)-(11) for different values of the free parameters. The initial condition is a point near to the matter dominated solution A1. Solid line is for decellaration parameter q, dashed line is for the jerk parameter j and dotted line is for the snap parameter s.