Asymptotic cosmological regimes in scalar-torsion gravity with a perfect fluid

We consider cosmological dynamics of nonminimally coupled scalar field in the scalar-torsion gravity in the presence of a hydrodynamical matter. Potential of the scalar field have been chosen as power-law with negative index, this type of potentials is usually used in quintessence scenarios. We identify several asymptotic regimes, including de Sitter, kinetic dominance, kinetic tracker and tracker solution and study conditions for their existence and stability. We show that for each combination of coupling constant and potential power index one of regimes studied in the present paper is stable to the future.


Introduction
Recently a new class of theories modifying General Relativity (GR) has started to attract great attention mainly in order to explain observation data indicating the accelerated expansion of the late Universe [1]. It is based on the rather old theory attributed to Einstein himself [2], however, forgotten for decades since it is in fact not a separate theory from the viewpoint of equation of motion, but rather a reformulation of GR. Instead of torsion-free Levi-Civita connections it uses curvature free Weitzenböck connections [3]. It is well known that zerocurvature connections allow for the existence of a path-independent definition of vector parallel transport, so this theory has got the name of the Teleparallel Equivalent of General Relativity (TEGR) [4,5,6,7]. What, however, has become clear rather recently is that the well-known modifications of GR, such as a theory of a scalar field nonminimally coupled to gravity, when constructed on the base of TEGR, lead to different equations of motion (for the reason of this see [8]). Such a theory has been intensively investigated in many papers during several recent years.
Generalizations of TEGR are usually constructed in the same way as modifications of GR. It is possible to generalize the Lagrangian replacing in it the torsion scalar T with a function f (T ) [9,10]. Different types of cosmological scenarios appear in scalar-torsion gravity, namely, in the class of models with nonminimal coupling between the torsion scalar and the scalar field φ of the form ξT F (φ), where ξ is a coupling constant, F (φ) is some function of a scalar field [11,12,13,14]. Recently some other modifications of TEGR including a nonminimal derivative coupling to torsion [15] and analogs of the Gauss-Bonnet invariant have appeared [16]. There are also modifications of TEGR with no direct analogs of GR modification [8].
In our recent paper [17] we applied dynamical system methods to scaler-torsion theory in order to find some cosmological asymptotic regimes and describe the corresponding phase portraits, focusing mainly on growing scalar field potentials. The main qualitative result of that paper is that scalar-torsion coupling leads to much less variety of possible dynamical regimes than scalar-curvature coupling [18,19]. In the present paper we give an heuristic explanation of this feature, as well as show that more possibilities appears for decreasing scalar fields potentials, usually used in quintessence models. As such models have been created for a description of the late Universe, when the usual matter is important as well, we add hydrodynamical matter to the scalar field. We will use units with = c = 1.
The paper is organized as follows: in Sect. 2 we present equation of motion in initial and expansion normalized variables, in Sect. 3 the analysis of stationary points and corresponding asymptotic solutions is given, and Sect. 4 gives a brief summary of results obtained.

Main equations
We start with describing basic objects of teleparallel gravity. In this theory the dynamical variables are the tetrad fields (also called the vierbein fields) e A (x µ ); here Greek indices relate to space-time and capital Latin indices belong to the tangent space-time. The metric tensor is expressed in terms of the tetrad as where η AB = diag(1, −1, −1, −1). For the definition of parallel transport of a vector the Weitzenböck connection [3] is used, Then the torsion tensor and the torsion scalar are given by On the other hand, the curvature of the Weitzenböck connection vanishes identically.
We consider the cosmological model with the following action: where e = √ −g is the determinant of the tetrad, K = 8πG, φ is a canonical scalar field, V (φ) its potential, and B(φ) its arbitrary nonminimal coupling with the torsion scalar T , S m is the matter action. In the classical scalar-curvature theory the action has the same form except for the curvature scalar R replacing torsion scalar T . For the spatially flat Friedmann-Robertson-Walker tetrad e A µ = diag(1, a(t), a(t), a(t)) (the corresponding metric is ds 2 = dt 2 − a 2 (t)dl 2 ) the system of field equations is [17] Here a(t) is the scale factor, H(t) ≡˙a a is the Hubble parameter, the prime denotes the derivative with respect to φ. We have used that in the chosen tetrad T = −6H 2 . The matter equation of state is p = ωρ, where ω ∈ [−1; 1].

The effective potential
In the standard scalar-curvature theory of a nonminimal coupling the conformal transformation to Einstein frame is often used. In the Einstein frame the theory is equivalent to GR with a scalar field as a source, so the evolution of the scalar field in an expanding Universe is, as usual, directed to the minimum of the potential, which is not true in the initial nonminimal formulation, called the Jordan frame. This conformal transformation is usually divided into two steps: a redifinition of the potential via V eff (φ) = V (φ)/U 2 (φ), where U(φ) = 1 + KξB(φ) in our notations, and a redifinition of the scalar field in order to get the canonical kinetic term. The combination denoted here as V eff (φ) is rather interesting by itself. It is a conformal invariant, so it is not changed by any conformal transformation [20]. Moreover, it can give important information as regards the dynamics of the system without the second step to the Einstein frame (which is usually much more technically complicated than very easily calculated first step). Namely, de Sitter solutions corresponds to minima of the effective potential, their stability is determined by sign of the second derivative of the effective potential in a standard way [19,20].
In the scalar-torsion theory the Einstein frame does not exist [21]. However, it is possible to introduce some analog of effective potential with the same properties as in the scalar-curvature theory. Indeed, taking into account (6) we rewrite (8) We see that the effective potential, of the form has a derivative with respect to φ which coincides with numerator of the last term in (9). The de Sitter solution H = H 0 , φ = φ 0 , in this model exists for From this system it is follows that Now we add small perturbations to the de Sitter solution: φ(t) = φ 0 + δφ,φ(t) = δφ, H(t) = H 0 + δH. Substituting these perturbations to (9) we get in the first order New variables s = δφ, r = δφ are introduced and the first-order system of differential equations is written in the formṡ = r, We find eigenvalues for the system (15) So, the de Sitter solution is stable in the minima of the effective potential, as expected. It is worth to note that the second equation in Eq. (13) is just the "balanced solution" studied in [22], so from computational point of view the effective potential gives nothing new. However, from heuristic point of view often it is much easier to visualize the locations of the minima of some function instead of doing calculations. For example, for positive power-law potential and coupling functions B(φ) ∼ φ N and V (φ) ∼ φ n de Sitter solution exists only for negative n with 0 < −n < N. This fact has been established computationally in [17] and becomes now a trivial consequence of the form of V eff (φ) = U(φ)V (φ).

Dimensionless variables
We introduce new dimensionless variables and also dimensionless parameters Choosing the power-law function B(φ) = φ N and the potential Note the useful relation between A, x, m, namely, It is useful to introduce auxiliary variables which are expressed through dimensionless variables and parameters from the system (6)-(8), From (Eq. 6) using (17) we get Taking the derivative of the variables y, m, A with respect to ln(a) ( ′ = d d ln a ) we obtain the following first-order system of differential equations: and finally substituting (19), (21), and (22) we get If N = 2 (b = 1), then instead of (19), and (21) we have and the first-order system of differential equations We consider only the case of N = 2 in the present paper.
3 Stationary points and corresponding regimes

Stationary points analysis
Solving the system (27) with vanishing left-hand sides, we find the following stationary points: 1. x = 0, y = 1, z = 0, m = 0, A = − c Kξ . We calculate the eigenvalues for the Jacobian matrix associated with the system (27) to find so this point is stable for non-phantom matter. For this point the quantity Y =Ḣ H 2 = 0 is found using (26) and, therefore,Ḣ = 0. Then we find the time dependence of the scale factor, As the coordinate of this fixed point is A = N φ N 1+Kξφ N → − n Kξ = 0, then the scalar field approaches a constant φ = φ 0 .
This is a de Sitter solution. The constants φ 0 and H 0 are found with the substitution of this solution to the initial system of Eqs. (6)- (8): Clearly, this regime exists only for 0 > n > −2. Since it is stable it can be used for describing late-time acceleration of our Universe. Definitely, the unnaturally low value of H 0 (in natural units) should transform to very low value of some of the parameters of the theory. Note, however, that in the theory under investigation we have more possibilities to get a very low H 0this may happen due to either very low V 0 , a very low ξ or a potential power index n being very close to −2. For the last possibility, if we denote n + 2 = ǫ << 1, we get H 0 ∼ ǫ √ V 0 ξ, so small corrections to the φ −2 potential would do the job.
Kξ . The corresponding eigenvalues are Using (26) the quantity Y at this stationary point is obtained to be Y stat = −3 − 2 √ 6ξ and now we can restore the time dependence of the Hubble parameter, H(t) = − 1 Ystat(t−t 0 ) , and the scale factor The corresponding behavior of the scalar field is obtained applying the definition of the variable m =φ Hφ , thenφ φ = m statȧ a , where m stat -the coordinate of a stationary point, 3. x = 1, y = 0, z = 0, m = − √ 6ξ, A = 2 Kξ . The eigenvalues for this fixed point are We get asymptotic behavior a(t) and φ(t) calculating the quantity Y stat = −3 + 2 √ 6ξ and using m = m stat = − √ 6ξ at this point, For ξ = 3 8 instead of this solution (35) we have the exponential one (because in this case Y = 0, The points 2 and 3 represent a situation when kinetic term of the scalar field dominates.
Kξ . The eigenvalues are calculated, Similar to the previous points, using Y stat = 1 2 ξ(4 − c 2 ) and m stat = −ξ(c + 2) we find a(t), φ(t) (c = n), (42) Kξ . We obtain the eigenvalues for this point, c)) . (43) give us the behavior a(t), φ(t), and ρ(t) (c = n) This is a tracker when the potential and kinetic terms of the scalar field as well as the density of the matter are constant with respect to each other.

Conditions of existence and stability of obtained solutions
With the substitution of the power-law solution a(t) = a 0 (t − t 0 ) α , φ(t) = φ 0 (t − t 0 ) αβ , and ρ(t) = ρ 0 (t − t 0 ) −3(1+ω)α to the initial system of equations (6), (7), (8), where B(φ) = φ 2 , V (φ) = V 0 φ n , we find the conditions of the existence and stability of the six solutions obtained from the previous subsection. We consider only the cases with n < 0, ξ > 0, and only positive values of φ. It is worth to recall that in general a fixed point does not necessary correspond to some solution of the initial system. First of all, all regimes found above (except the de Sitter solution, which is an exact, and not only an asymptotic solution) exist in the φ → ∞ limit where we can neglect the Einstein term in comparison with the term originating from a nonminimal coupling in the denominator of the expansion in normalized variables we use. When the field is small (which is realized near a cosmological singularity for decreasing power-law potentials) the Universe expands according to GR, since the correction terms are less important. Such regimes for quintessence potentials are well known and are not included in our analysis here. As for the large φ regimes, studied in the present paper, we require that the omitted terms in the equations of motion be negligible for large φ. Otherwise the regime cannot be realized as an asymptotic solution. For example, suppose we have a vacuum asymptotic solution. If the influence of matter grows with growing φ, this means that the asymptotic solution corresponding to this particular expansion dynamics is absent if we add any amount of matter. This does not mean that such a regime has no physical meaning at all -the Universe may follow it as transient one, if the amount of matter is small enough. However, such a situation does not belong to the asymptotic regimes, and we will not list it in the following.
The analysis described above (with stability results got from the corresponding eigenvalues) leads to the following results for vacuum solutions summarized in Table 1 (points 4 and 6, being non-vacuum, are not included).
We can see from this table that in the vacuum case there are three possible future stable regimes: • For n > −2 a de Sitter solution exists and is stable.
When the matter is taken into account, the situation is described in the Table 2 (note that the de Sitter solution in the presence of matter is an asymptotic solution).

Conclusions
In the present paper we considered the evolution of the Universe in a teleparallel version of a nonminimally coupled scalar field theory, focusing mostly on quintessence scalar field potentials. In our previous paper [17] we already pointed out the radical differences in the dynamics of a "classical" and teleparallel cases of this theory. Here we provided a qualitative explanation of this difference using the notion of an effective scalar field potential. It appears as a first step from a Jordan to an Einstein frame (to be followed by redefinition of the scalar field in order to transform the kinetic term to the canonical form) in the standard theory, though it already gives an intuitively clear picture of possible dynamical regimes without this second step. It is interesting that in the teleparallel version, where the Einstein frame does not exist at all [21], the effective potential nevertheless can be introduced in the same manner. The difference in the forms of this effective potential (being V (φ)/U 2 (φ) for the classical case and V (φ)U(φ) for the teleparallel case, where V (φ) is the scalar field potential and U(φ) is the coupling function which includes the Einstein-Hilbert term) results in very different dynamics for the same scalar field potentials.
For the standard scalar-curvature coupling, the presence of U(φ) [which is usually taken in the form U(φ) = 1 + 8πGξφ 2 ] in the denominator leads to run-away solution for growing V (φ), corresponding to V eff (φ) decreasing for large φ. This may happen if U 2 (φ) grows more rapidly than V with φ [19]. In the boundary case of U 2 (φ) ∼ V (φ) for large φ the effective potential is asymptotically flat, which is the reason for the existence of viable Higgs inflation models [23].
None of these features can exist for the effective potential in the form of U(φ)V (φ) as in scalar-torsion theory. For the usual quadratic coupling (we study only this form of the coupling in the present paper) V eff (φ) is growing for an arbitrarly growing V (φ), so the dynamics inevitably pushes the scalar field toward zero, as we have seen in our previous paper. On the contrary, for the decreasing potentials studied in the present paper the effective potential can have a minimum, resulting in de Sitter solution. For a quadratic coupling this happens if V (φ) scales less steep than φ −2 . This condition for a de Sitter solution to exist has been obtained earlier in [17] and in a more general context in [22], and we now have a very easy way to explain it.
Steeper potentials result in unlimited growing of the scalar field, similar to the quintessence scenario in GR. We have identified a tracker solution which is always an attractor in its range of existence. It is shown that when the coupling constant ξ is low enough for the tracker to exist (see the corresponding condition in Table 2), a vacuum regime is stable, so the scalar field dominates at late time. We also list several unstable asymptotic regimes. A full description of the scalar field dynamics requires matching solutions studied in the present paper with GR solutions valid for small φ. We leave this for future work.