## Dynamics of Generalized Tachyon Field in Teleparallel Gravity

Behnaz Fazlpour a 1 and Ali Banijamali b 2

- a Department of Physics, Babol Branch, Islamic Azad University, Babol, Iran
- b Department of Basic Sciences, Babol University of Technology, Babol, Iran

## Abstract

We study dynamics of generalized tachyon scalar field in the framework of teleparallel gravity. This model is an extension of tachyonic teleparallel dark energy model which has been proposed in [26]. In contrast with tachyonic teleparallel dark energy model that has no scaling attractors, here we find some scaling attractors which means that the cosmological coincidence problem can be alleviated. Scaling attractors present for both interacting and non-interacting dark energy, dark matter cases.

PACS numbers: 95.36.+x, 98.80.-k, 04.50.kd

Keywords: Generalized Tachyon Field; Teleparallel gravity; Phase-space analysis.

1 b.fazlpour@umz.ac.ir

2 a.banijamali@nit.ac.ir

## 1 Introduction

The usual proposal to explain the late-time accelerated expansion of our universe is an unknown energy component, dubbed as dark energy. The natural choice and most attractive candidate for dark energy is the cosmological constant but it is not well accepted because of the cosmological constant problem [1] as well as the age problem [2]. Thus, many dynamical dark energy models as alternative possibilities have been proposed. Quintessence, phantom, k-essence, quintom and tachyon field are the most familiar dark energy models in the literature (for reviews on dark energy models, see [3]). The tachyon field arising in the context of string theory [4] and its application in cosmology both as a source of early inflation and late-time cosmic acceleration has been extensively studied [5-8].

The so-called 'Teleparallel Equivalent of General Relativity' or Teleparallel Gravity was first constructed by Einstein [9-12]. In this formulation one uses the curvature-less weitzenbock connection instead of the torsion-less Levi-Civita connection. The relevant lagrangian in teleparallel gravity is the torsion scalar T which is constructed by contraction of the torsion tensor. We recall that the Einstein-Hilbert Lagrangian R is constructed by contraction of the curvature tensor. Since teleparallel gravity with torsion scalar as lagrangian density is completely equivalent to a matter-dominated universe in the framework of general relativity, it can not be accelerated. Thus one should generalize teleparallel gravity either by replacing T with an arbitrary function -the so-called f ( T ) gravity [13-15] or by adding dark energy into teleparallel gravity allowing also a non-minimal coupling between dark energy and gravity. Note that both approaches are inspired by the similar modifications of general relativity i.e. f ( R ) gravity [16, 17] and non-minimally coupled dark energy models in the framework of general relativity [18-20].

Recently Geng et al. [21, 22] have included a non-minimal coupling between quintessence and gravity in the context of teleparallel gravity. This theory has been called 'teleparallel dark energy' and its dynamics was studied in [23-25]. Tachyonic teleparallel dark energy is a generalization of teleparallel dark energy by inserting a non-canonical scalar field instead of quintessence in the action [26]. Phase-space analysis of this model has been investigated in [27]. On the other hand, there is no physical argument to exclude the interaction between dark energy and dark matter. The interaction between these completely different component of our universe has same important consequences such as addressing the coincidence problem [28]. In this paper we consider generalized tachyon field as responsible for dark energy in the framework of teleparallel gravity. We will be interested in performing a dynamical analysis of such a model in FRW space time. In such a study we investigate our model for both interacting and non-interacting cases. The basic equations are presented in section 2. In section 3 the evolution equations are translated in the language of the autonomous dynamical system by suitable transformation of the basic variables. Subsection 3.1 deals with phase-space analysis as well as the cosmological implications of the equilibrium points of the model in non-interacting dark energy dark matter case. In subsection 3.2 an interaction between dark energy and dark matter has been considered an critical points and their behavior extracted. Section 5 is devoted to a short summary of our results.

## 2 Basic Equations

Our model is described by the following action as a generalization of tachyon teleparallel dark energy model [26],

$$S = \int d ^ { 4 } x e \left [ \mathcal { L } _ { T } + \mathcal { L } _ { \varphi } + \mathcal { L } _ { m } \right ], \\ \mathcal { L } _ { T } = \frac { T } { 2 \kappa ^ { 2 } }, \\ \mathcal { L } _ { \varphi } = \xi f ( \varphi ) T - V ( \varphi ) ( 1 - 2 X ) ^ { \beta }, \\ \mathrm e \, \text{the orthonormal components of the tet} \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ }$$

$$\mathcal { L } _ { \varphi } = \xi f ( \varphi ) T - V ( \varphi ) ( 1 - 2 X ) ^ { \beta }, \\ = \mathcal { \L i } \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \cdots \quad T \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \collections \quad.$$

where e = det e ( i µ ) = √ -g ( e i µ are the orthonormal components of the tetrad) while T 2 κ 2 is the Lagrangian of teleparallelism with T as the torsion scalar (for an introductory review of teleparallelism see [11]). L ϕ

shows a non-minimal coupling of generalized tachyon field ϕ with gravity in the framework of teleparallel gravity and X = 1 2 ∂ µ ϕ∂ µ ϕ . The second part in L ϕ is the Lagrangian density of the generalized tachyon field which has been studied in Ref [29]. f ( ϕ ) is the non-minimal coupling function, ξ is a dimensionless constant measuring the non-minimal coupling and L m is the matter Lagrangian. For β = 1 2 our model reduced to tachyonic teleparallel dark energy discussed in [26]. Here we consider the case β = 2 for two reasons. The first is that for arbitrary β our equations will be very complicated and one can not solve them analytically and the second is that for β = 2 we will obtain interesting physical results as we will see below.

Furthermore, due to complexity of tachyon dynamics Ref [30] has proposed an approach based on a redefinition of the tachyon field as follows,

$$\varphi \to \phi = \int d \varphi \sqrt { V ( \varphi ) } \Leftrightarrow \partial \varphi = \frac { \partial \phi } { \sqrt { V ( \phi ) } }. \\ \text{closed autonomous system and perform the phase-space analysis of the model we apply}$$

In order to obtain a closed autonomous system and perform the phase-space analysis of the model we apply (2) in (1) for β = 2 that leads to the following action:

$$S = \int d ^ { 4 } x e \left [ \frac { T } { 2 \kappa ^ { 2 } } + \xi f ( \phi ) T - V ( \phi ) ( 1 - \frac { 2 X } { V ( \phi ) } ) ^ { 2 } + \mathcal { L } _ { m } \right ].$$

In a spatially-flat FRW space-time,

$$d s ^ { 2 } = d t ^ { 2 } - a ^ { 2 } ( t ) ( d r ^ { 2 } + r ^ { 2 } d \Omega ^ { 2 } ), \\. \quad \.$$

and a vierbein choice of the form e i µ = diag (1 , a, a, a ), the corresponding Friedmann equations are given by,

$$H ^ { 2 } & = \frac { 1 } { 3 } ( \rho _ { \phi } + \rho _ { m } ), & & ( 5 )$$

$$\dot { H } = - \frac { 1 } { 2 } \left ( \rho _ { \phi } + P _ { \phi } + \rho _ { m } + P _ { m } \right ), \\ \intertext { \text{lo normal and a dot stand for the derivative with respect to the cosmic time} }$$

where H = ˙ a a is the Hubble parameter and a dot stands for the derivative with respect to the cosmic time t . In these equations, ρ m and P m are the matter energy density and pressure respectively.

The effective energy density and pressure of generalized tachyon dark energy read,

$$\rho _ { \phi } = V ( \phi ) + 2 \dot { \phi } ^ { 2 } - 3 \frac { \dot { \phi } ^ { 4 } } { V ( \phi ) } - 6 \xi H ^ { 2 } f ( \phi ),$$

and

$$P _ { \phi } = - V ( \phi ) + 2 \xi ( 3 H ^ { 2 } + 2 \dot { H } ) f ( \phi ) + 1 0 \xi H f _ {, \phi } \dot { \phi } + \dot { \phi } ^ { 2 } ( 2 - \frac { \dot { \phi } ^ { 2 } } { V ( \phi ) } ),$$

where f ,φ = df dφ .

The equation of motion of the scalar field can be obtained by variation of the action (3) with respect to φ ,

$$\ddot { \phi } + 3 \mu ^ { - 2 } \nu ^ { 2 } H \dot { \phi } + \frac { 1 } { 4 } \nu ^ { 2 } ( 1 + \frac { 3 \dot { \phi } ^ { 4 } } { V ^ { 2 } ( \phi ) } ) V _ {, \phi } + 6 \xi \nu ^ { 2 } H ^ { 2 } f _ {, \phi } = - \frac { Q } { \dot { \phi } },$$

with Q a general interaction coupling term between dark energy and dark matter, µ = 1 √ 1 -2 X V and ν = 1 √ 1 -6 X V . In (7), (8) and (9) we have used the useful relation,

$$T & = - 6 H ^ { 2 }, & ( 1 0 )$$

which simply arises from the calculation of torsion scalar for the FRW metric (4). The scalar field evolution (9) expresses the continuity equation for the field and matter as follows

$$\dot { \rho } _ { \phi } + 3 H ( 1 + \omega _ { \phi } ) \rho _ { \phi } = - Q,$$

$$\dot { \rho } _ { m } + 3 H ( 1 + \omega _ { m } ) \rho _ { m } = Q,$$

where ω φ = P φ ρ φ is the equation of state parameter of dark energy which is attributed to the scalar field φ . The barotropic index is defined by γ ≡ 1 + ω m with 0 &lt; γ &lt; 2.

Although, dynamics of tachyonic teleparallel dark energy has been studied in [27], no scaling attractors found. Here we are going to perform a phase-space analysis of generalized tachyonic teleparallel dark energy and as we will see below some interesting scaling attractors appear in such theory.

## 3 Cosmological Dynamics

In order to perform phase-space and stability analysis of the model, we introduce the following auxiliary variables:

$$x \equiv \frac { \dot { \phi } } { \sqrt { V } }, \ y \equiv \frac { \sqrt { V } } { \sqrt { 3 } H }, \ u \equiv \sqrt { f }. \quad \quad ( 1 3 )$$

The auxiliary variables allow us to straightforwardly obtain the density parameter of dark energy and dark matter

$$\Delta _ { m } \equiv \frac { \rho _ { m } } { 3 H ^ { 2 } } = 1 - \Omega _ { \phi },$$

$$\Omega _ { \phi } \equiv \frac { \rho _ { \phi } } { 3 H ^ { 2 } } & = \mu ^ { - 2 } y ^ { 2 } ( 1 + 3 x ^ { 2 } ) - 2 \xi u ^ { 2 }, \\ & \wedge \rho _ { m } \quad. \quad \wedge \quad \ w.$$

while the equation of state of the field reads

ω

φ

φ

≡

P

φ

ρ

where α ≡ f ,φ √ f and

$$r _ { \varphi } ^ { r _ { \varphi } } = \frac { - \mu ^ { - 4 } y ^ { 2 } + 2 \xi u [ \frac { 5 \sqrt { 3 } } { 3 } \alpha x y + u ( 1 - \frac { 2 } { 3 } s ) ] } { \mu ^ { - 2 } y ^ { 2 } ( 1 + 3 x ^ { 2 } ) - 2 \xi u ^ { 2 } },$$

$$s = - \frac { \dot { H } } { H ^ { 2 } } = ( 2 \xi u ^ { 2 } + 1 ) ^ { - 1 } \left [ 5 \sqrt { 3 } \alpha \xi u x y + 6 \mu ^ { - 2 } x ^ { 2 } y ^ { 2 } - \frac { 3 } { 2 } \gamma \mu ^ { - 2 } y ^ { 2 } ( 1 + 3 x ^ { 2 } ) \right ] + \frac { 3 \gamma } { 2 }. \quad \quad ( 1 7 ) \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

Another quantities with great physical significance namely the total equation of state parameter and the deceleration parameter are given by

$$\check { \omega } _ { t o t } & \equiv \frac { P _ { \phi } + P _ { m } } { \rho _ { \phi } + \rho _ { m } } = \mu ^ { - 2 } y ^ { 2 } ( 4 x ^ { 2 } - \gamma ( 1 + 3 x ^ { 2 } ) ) + 2 \xi u [ \frac { 5 \sqrt { 3 } } { 3 } \alpha x y + u ( \gamma - \frac { 2 } { 3 } s ) ] + \gamma - 1, \quad \quad ( 1 8 )$$

and

˙

H

1

H

q

≡ -

1

-

2

=

2

+ 3

2

ω

tot

$$\cdots \quad - \quad - \quad - \\ = \frac { 3 } { 2 } \mu ^ { - 2 } y ^ { 2 } ( 4 x ^ { 2 } - \gamma ( 1 + 3 x ^ { 2 } ) ) + \xi u [ 5 \sqrt { 3 } \alpha x y + u ( 3 \gamma - 2 s ) ] + \frac { 3 \gamma } { 2 } - 1. \quad \ ( 1 9 ) \\ \alpha \text{ auxiliary variables} ( 1 \cdots \, \text{$1\cup$ the solution amounts} ( 5 \cdots \, \text{$6\cup$ and} \, ( 0 \cup \, \text{can be exact as mathematical event} }$$

Using auxiliary variables (13) the evolution equations (5), (6) and (9) can be recast as a dynamical system of ordinary differential equations

$$\mu _ { \nu } ^ { - } \mu _ { \nu } ^ { - } \mu _ { \nu } ^ { - } \mu _ { \nu } ^ { 2 } ( 1 + 3 x ^ { 4 } ) y - 4 \alpha \xi \nu ^ { 2 } u y ^ { - 1 } - 2 \sqrt { 3 } \mu ^ { - 2 } \nu ^ { 2 } x ] - \hat { Q },$$

$$y ^ { \prime } & = \left ( \, - \, \frac { \sqrt { 3 } } { 2 } \lambda x y + s \right ) y,$$

$$u ^ { \prime } & = \frac { \sqrt { 3 } \alpha x y } { 2 },$$

where ˆ = Q Q ˙ φH √ V ( φ ) , λ ≡ -V ,φ κV and prime in equations (20)-(22) denotes differentiation with respect to the so-called e-folding time N = ln a .

From now we concentrate on exponential scalar field potential of the form V = V e 0 -kλφ and the non-minimal coupling function of the form f ( φ ) ∝ φ 2 . These choices lead to constant λ and α respectively.

The next step is the introduction of interaction term Q to obtain an autonomous system out of equations (20)-(22). The fixed points ( x , y c c , u c ) for which x ′ = y ′ = u ′ = 0 depend on the choice of the interaction term Q and two general possibilities will be treated in the sequel. The stability of the system at a fixed point can be obtained from the analysis of the determinant and trace of the perturbation matrix M . Such a matrix can be constructed by substituting linear perturbations x → x c + δx, y → y c + δy and u → u c + δu about the critical point ( x , y c c , u c ) into the autonomous system (20)-(22). The 3 × 3 matrix M of the linearized perturbation equations of the autonomous system is shown in appendix A. Therefor, for each critical point we examine the sign of the real part of the eigenvalues of M . According to the usual dynamical system analysis, if the eigenvalues are real and have opposite signs, the corresponding critical point is a saddle point. A fixed point is unstable if the eigenvalues are positive and it is stable for negative real part of the eigenvalues.

In the following subsections we will study the dynamics of generalized tachyon field with different interaction term Q . Without lose of generality we assume γ = 1 for simplicity.

## 3.1 The case for Q = 0

The first case Q = 0 clearly means there is no interaction between dark energy and background matter. In this case, there are two critical points presented in Table 1. From equations (14) and (16) one can obtain the corresponding values of density parameter Ω φ and equation of state of dark energy ω φ at each point. Also, using equation (19) we can find the condition required for acceleration ( q &lt; 0) at each point. These parameters and conditions have been shown in Table 1. The stability and existence conditions of critical points A 10 and A 2 are presented in Table 2. We mention that the corresponding eigenvalues of perturbation matrix M at critical points A 10 and A 2 are considerably involved and here we do not present their explicit expressions but we can find sign of them numerically.

Critical point A 1 : This critical point is a scaling attractor if ξ &gt; 0, α &lt; 0 and λ &lt; 0. Thus, it can give the hope alleviate the cosmological coincidence problem. A 1 is a saddle point for α &gt; 0 and ξ λ &gt; 0.

Critical point A 2 : A 2 can also be a scaling attractor of the model or a saddle point under the same conditions as for A 1 .

In figure 1 we have chosen the values of the parameters ξ , λ and α , such that A 1 become a stable attractor of the model. Plots in figure 1 show the phase-space trajectories on x -y , x -u and u -y planes from left to right respectively. The same plots are shown in figure 2 for critical point A 2 . Note that the values of the parameter have chosen in the way that A 2 become a stable point of the model. In figure 3, the corresponding 3 dimensional phase-space trajectories of the model have been presented. One can see that A 1 and A 2 are stable attractor of the model in the left and right plots respectively.

| label   | ( x c , y c , u c )     | Ω φ                          | ω φ                                                   | acceleration                                                                                                                         |
|---------|-------------------------|------------------------------|-------------------------------------------------------|--------------------------------------------------------------------------------------------------------------------------------------|
| A 1     | 0 , 2 √ λ 1 , λλ 1 2 αξ | 4 λ 1 ( 1 - λ 2 8 ξα 2 λ 1 ) | 16 ξ 2 α 4 ( λ 2 λ 2 1 +2 ξα 2 ) ( λ 2 λ 1 - 8 ξα 2 ) | λ 2 > ξα 2 (4 λ 1 - 1) λ 2 1 ( 1 + √ 1 - 32 λ 1 (1 - 4 λ 1 ) 2 ) or λ 2 < ξα 2 (4 λ 1 - 1) λ 2 1 ( 1 - √ 1 - 32 λ 1 (1 - 4 λ 1 ) 2 ) |
| A 2     | 0 , 2 √ λ 2 , λλ 2 2 αξ | 4 λ 2 ( 1 - λ 2 8 ξα 2 λ 2 ) | 16 ξ 2 α 4 ( λ 2 λ 2 2 +2 ξα 2 ) ( λ 2 λ 2 - 8 ξα 2 ) | λ 2 > ξα 2 (4 λ 2 - 1) λ 2 2 ( 1 + √ 1 - 32 λ 2 (1 - 4 λ 2 ) 2 ) or λ 2 < ξα 2 (4 λ 2 - 1) λ 2 2 1 - √ 1 - 32 λ 2 (1 - 4 λ 2 ) 2     |

(

)

Table 1: Location of the critical points and the corresponding values of the dark energy density parameter Ω φ and equation of state ω φ and the condition required for an accelerating universe for Q = 0. Here λ 1 = α ( 4 αξ + √ 16 α ξ 2 2 -2 λ ξ 2 ) λ 2 and λ 2 = α ( 4 αξ -√ 16 α ξ 2 2 -2 λ ξ 2 ) λ 2 .

Table 2: Stability and existence conditions of the critical points of the model for Q = 0.

| label   | stability                                                                 | existence                      |
|---------|---------------------------------------------------------------------------|--------------------------------|
| A 1     | saddle point if α > 0 and ξ λ > 0 stable point if ξ > 0 , α < 0 and λ < 0 | for all ξ < 0 or ξ ≥ λ 2 8 α 2 |
| A 2     | saddle point if α < 0 and ξ λ < 0 stable point if ξ > 0 , α > 0 and λ > 0 | for all ξ < 0 or ξ ≥ λ 2 8 α 2 |

| label   | ( x c , y c , u c )        | Ω φ                    | ω φ                                                | acceleration                            |
|---------|----------------------------|------------------------|----------------------------------------------------|-----------------------------------------|
| B 1     | 0 , 2 √ 2 λαξθ 1 λ , θ 1   | 2 ξθ 1 ( 4 α λ - θ 1 ) | 2 ( ξθ 3 1 + 2 α λ ) ( θ 1 - 4 α λ ) (1+2 ξθ 2 1 ) | λ > 8 α ( ξθ 2 1 - 1) θ 1 (8 ξθ 2 1 +1) |
| B 2     | 0 , 2 √ 2 λαξθ 2 λ , θ 2   | 2 ξθ 2 ( 4 α λ - θ 2 ) | 2 ( ξθ 3 2 + 2 α λ ) ( θ 2 - 4 α λ ) (1+2 ξθ 2 2 ) | λ > 8 α ( ξθ 2 2 - 1) θ 2 (8 ξθ 2 2 +1) |
| B 3     | 0 , - 2 √ 2 λαξθ 1 λ , θ 1 | 2 ξθ 1 ( 4 α λ - θ 1 ) | 2 ( ξθ 3 1 + 2 α λ ) ( θ 1 - 4 α λ ) (1+2 ξθ 2 1 ) | λ > 8 α ( ξθ 2 1 - 1) θ 1 (8 ξθ 2 1 +1) |
| B 4     | 0 , - 2 √ 2 λαξθ 2 λ , θ 2 | 2 ξθ 2 ( 4 α λ - θ 2 ) | 2 ( ξθ 3 2 + 2 α λ ) θ 2 - 4 α λ (1+2 ξθ 2 2 )     | λ > 8 α ( ξθ 2 2 - 1) θ 2 (8 ξθ 2 2 +1) |

(

)

Table 3: Location of the critical points and the corresponding values of the dark energy density parameter Ω φ and equation of state ω φ and the condition required for an accelerating universe for Q = βκρ m ˙ . φ Here θ 1 = ( 4 αξ + √ 16 α ξ 2 2 -2 λ ξ 2 ) 2 λξ and θ 2 = ( 4 αξ -√ 16 α ξ 2 2 -2 λ ξ 2 ) 2 λξ .

## 3.2 The case for Q = βκρ m ˙ φ

This deals with the most familiar interaction term extensively considered in the literature (see e.g. [25, 31-34]). Here ˆ Q in terms of auxiliary variables is ˆ Q = √ 3 βy -1 Ω m . Inserting such an interaction term in equations (20)-(22) and setting the left hand sides of the equations to zero lead to the critical points B 1 , B 2 , B 3 and B 4 presented in Table 3. In the same table we have provided the corresponding values of Ω φ and ω φ as well as the condition needed for accelerating universe at each fixed points.

The stability and existence conditions for each point presented in Table 4. Since the corresponding eigenvalues of the fixed points are complicated we do not give them here but one can obtain their signs numerically and so concludes about the stability properties of the critical points.

Critical point B 1 : This point exists for λ &gt; 0 and ξ ≥ λ 2 8 α 2 . However it is an unstable saddle point.

Critical point B 2 : The critical point B 2 exists for λ &gt; 0 and ξ &lt; 0 or ξ ≥ λ 2 8 α 2 . This point is a scaling attractor of the model if α &gt; 0 and ξ &gt; 0. Figure 4 shows clearly such a behavior of the model for suitable choices of ξ , λ and α .

Critical point B 3 : This point exists for negative values of λ and ξ or when ξ ≥ λ 2 8 α 2 . Also, it is a stable point if α &lt; 0 and ξ &gt; 0 and a saddle point if α &gt; 0 and ξ &lt; 0. The values of parameters have been chosen in figure 5 such that B 3 become a attractor of the model as it is clear from phase-space trajectories.

Critical point B 4 : The point B 4 exists for λ &lt; 0 and ξ ≥ λ 2 8 α 2 . It is a stable point if α &lt; 0 and ξ &gt; 0. In figure 6 values of the parameter ξ and α are those satisfy these constraints and so B 4 becomes a attractor point for phase-plane trajectories. The corresponding 3-dimensional phase-space trajectories of the model for attractor points B 2 (left), B 3 (middle) and B 4 (right) are plotted in figure 7.

## 4 Conclusion

A model of dark energy with non-minimal coupling of quintessence scalar field with gravity in the framework of teleparallel gravity was called teleparallel dark energy [21]. If one replaces quintessence by tachyon field in such a model then tachyonic teleparallel dark energy will be constructed [26].

Moreover, although dark energy and dark matter scale differently with the expansion of our universe, ac-

Table 4: Stability and existence conditions of the critical points of the model for Q = βκρ m ˙ . φ

| label   | stability                                                       | existence                                    |
|---------|-----------------------------------------------------------------|----------------------------------------------|
| B 1     | saddle point if α > 0 and ξ > 0                                 | λ > 0 and ξ ≥ λ 2 8 α 2                      |
| B 2     | saddle point if α < 0 and ξ < 0 stable point if α > 0 and ξ > 0 | λ > 0 and for all ξ < 0 or ξ λ 2 2           |
| B 3     | saddle point if α > 0 and ξ < 0 stable point if α < 0 and ξ > 0 | ≥ 8 α λ < 0 and for all ξ < 0 or ξ λ 2 8 α 2 |
| B 4     | stable point if α < 0 and ξ > 0                                 | ≥ λ < 0 and ξ ≥ λ 2 8 α 2                    |

cording to the observations [35] we are living in an epoch in which dark energy and dark matter densities are comparable and this is the well-known cosmological coincidence problem [24]. This problem can be alleviated in most dark energy models via the method of scaling solutions in which the density parameters of dark energy and dark matter are both non-vanishing over there.

In this paper we investigated the phase-space analysis of generalized tachyon cosmology in the framework of teleparallel gravity. Our model described by action (1) which generalizes tachyonic teleparallel dark energy model proposed in [26]. We found some scaling attractors in our model for the case β = 2. These scaling attractors are A 1 and A 2 when there is no interaction between dark energy and dark matter. B 2 , B 3 and B 4 are scaling attractors in the case that dark energy interacts with dark matter through the interacting term Q = βκρ m ˙ . φ Our results show that generalized tachyon field represents interesting cosmological behavior in compare with ordinary tachyon fields in the framework of teleparallel gravity because there is no scaling attractor in the latter model. So, generalized tachyon field gives us the hope that cosmological coincidence problem can be alleviated without fine-tunings. One can study our model for different kinds of potential and other famous interaction term between dark energy and dark matter.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-d234e64f-ef31-45d6-9418-6b044f6a661c.jpg)

Figure 1: From left to right, the projections of the phase-space trajectories on the x -y , x -u and u -y planes with ξ = 0 5, . λ = -0 6 and . α = -2 for Q = 0. For these values of the parameters, point A 1 is a stable attractor of the model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-15ccf6b5-6756-4b3f-87b6-c16542a0c8f5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-0a7b842f-f0fe-4829-98c9-60e8f5c0b01a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-c617fe38-0771-4312-aff8-c0ddc498e016.jpg)

Figure 2: From left to right, the projections of the phase-space trajectories on the x -y , x -u and u -y planes with ξ = 0 5, . λ = 0 6 and . α = 2 for Q = 0. For these values of the parameters, point A 2 is a stable attractor of the model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-d59c8723-07f2-48fe-b6a1-f744ad92c9dd.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-b989f47b-9a54-49f2-b021-2c71447a378b.jpg)

Figure 3: 3-dimensional phase-space trajectories of the model for Q = 0 with stable attractors A 1 (left) and A 2 (right) . The values of the parameters are those mentioned in figure 1 and 2 respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-1c62754f-938b-473d-9239-79b3d23ada20.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-fd227e59-6cf6-4cf1-8169-5b047732f72a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-ec6c02cb-ed7f-4fb4-9c1e-95926c9abe21.jpg)

Figure 4: From left to right, the projections of the phase-space trajectories on the x -y , x -u and u -y planes with ξ = 0 5, . λ = 0 6, . α = 2 and β = 1 5 for . Q = βκρ m ˙ . φ For these values of the parameters, point B 2 is a stable attractor of the model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-6c5b9172-2d57-4865-9a77-92cdf8baeb6c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-fc875f31-b538-413b-9dd0-0e889962c1e5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-12c8c6a9-a192-40ff-9663-fab516cfc7a7.jpg)

Figure 5: From left to right, the projections of the phase space trajectories on the x -y , x -u and u -y planes with ξ = 0 5, . λ = -0 6, . α = -2 and β = 1 5 for . Q = βκρ m ˙ . φ For these values of the parameters, point B 3 is a stable attractor of the model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-9c08d76b-bf3c-486d-b74b-c36d53c45e4f.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-7b3f22cc-735b-451a-b309-ae832fe3527b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-9ad49114-5bf2-4502-993a-57e9de9ae4e8.jpg)

Figure 6: From left to right, the projections of the phase space trajectories on the x -y , x -u and u -y planes with ξ = 0 5, . λ = -0 6, . α = -2 and β = 1 5 for . Q = βκρ m ˙ . φ For these values of the parameters, point B 4 is a stable attractor of the model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-405ef7e1-f946-4305-913f-cb50c8e82f13.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-a69c76ad-2377-419f-a6eb-c52cd9abecf7.jpg)

Figure 7: 3-dimensional phase-space trajectories of the model for Q = βκρ m ˙ φ with stable attractors B 2 (left), B 3 (middle) and B 4 (right). The values of the parameters are those mentioned in figure 4, 5 and 6 respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-e874b83a-026d-4b87-b633-be71a935828c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-9b0ac53a-fcce-4644-a4cf-1c7e4c9b14ac.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0203v1-figure-3d341446-66a4-4799-986f-cb43821c52df.jpg)

## 5 Appendix: Perturbation Matrix Elements

The elements of 3 × 3 matrix M of the linearized perturbation equations for the real and physically meaningful critical points ( x , y c c , u c ) of the autonomous system (20)- (22) read,

$$& \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ M _ { 1 1 } = 3 \nu _ { c } ^ { 2 } \left ( \frac { \sqrt { 3 } } { 2 } \lambda x _ { c } y _ { c } ( 2 x _ { c } ^ { 2 } + \nu _ { c } ^ { 2 } ( 1 + 3 x _ { c } ^ { 4 } ) ) - 6 \mu _ { c } ^ { - 2 } \nu _ { c } ^ { 2 } x _ { c } ^ { 2 } - 4 \sqrt { 3 } \alpha \xi u _ { c } x _ { c } \nu _ { c } ^ { 2 } y _ { c } ^ { - 1 } \right ) + \sqrt { 3 } \lambda x _ { c } y _ { c } - 3 + M _ { 1 1 }, \ \ ( 2 3 )$$

$$M _ { 1 2 } = \frac { \sqrt { 3 } } { 4 } \left ( \lambda ( 2 x _ { c } ^ { 2 } + \nu _ { c } ^ { 2 } ( 1 + 3 x _ { c } ^ { 4 } ) ) + 8 \alpha \xi u _ { c } \nu _ { c } ^ { 2 } y _ { c } ^ { - 2 } \right ) + \mathcal { M } _ { 1 2 },$$

$$M _ { 1 3 } = - 2 \sqrt { 3 } \alpha \xi \nu _ { c } ^ { 2 } y _ { c } ^ { - 1 } + \mathcal { M } _ { 1 3 },$$

$$M _ { 2 1 } = \frac { 2 y _ { c } ^ { 2 } \left ( \sqrt { 3 } \alpha \xi u _ { c } + 3 x _ { c } y _ { c } \nu _ { c } ^ { - 2 } \right ) } { \left ( 2 \xi u _ {$$

$$M _ { 2 2 } = \frac { 2 y _ { c } ( - \frac { 9 } { 4 } \mu _ { c } ^ { - 4 } y _ { c } + 5 \sqrt { 3 } \alpha \xi x _ { c } u _ { c } ) } { ( 2 \xi u _ { c$$

$$M _ { 2 3 } = \frac { 6 \xi u _ { c } y _ { c } ^ { 2 } ( - \frac { 1 0 \sqrt { 3 } } { 3 } \alpha \xi u _ { c } x _ { c } + \mu _ { c } ^ { - 4 } y _ { c }$$

$$M _ { 3 1 } = \frac { \sqrt { 3 } \alpha y _ { c } } { 2 }, \ M _ { 3 2 } = \frac { \sqrt { 3 } \alpha x _ { c } } { 2 }, \ M _ { 3 3 } = 0, \quad \quad \quad$$

where in the case of Q = 0 we have M 11 = M 12 = M 13 = 0 and in the case of Q = βκρ m ˙ φ we have

$$\mathcal { M } _ { 1 1 } = 4 \sqrt { 3 } \beta \nu _ { c } ^ { - 2 } x _ { c } y _ { c },$$

$$\mathcal { M } _ { 1 2 } = 2 \sqrt { 3 } \beta \mu _ { c } ^ { - 2 } ( 1 + 3 x _ { c } ^ { 2 } ), \\ \mathcal { M } _ { 1 3 } = - 4 \sqrt { 3 } \beta \xi u _ { c } y _ { c } ^ { - 1 }. \\ \text{es of the matrix } M \text{ for each critical noint. one determines its stability conditions.}$$

Examining the eigenvalues of the matrix M for each critical point, one determines its stability conditions.

## References

- [1] S. Weinberg, Rev. Mod. Phys. 61 , 1 (1989); S. M. Carroll, Living Rev. Rel. 4 , 1 (2001).
- [2] R.-J. Yang and S. N. Zhang, Mon. Not. R. Astron. Soc. 407 , 1835 (2010).
- [3] Y. -F. Cai, E. N. Saridakis, M. R. Setare and J. -Q. Xia, Phys. Rept. 493 , 1 (2010); K. Bamba, S. Capozziello, S. Nojiri and S. D. Odintsov, Astrophys. Space Sci. 342 , 155 (2012).
- [4] Y. -F. Cai, E. N. Saridakis, M. R. Setare and J. -Q. Xia, Phys. Rept. 493 , 1 (2010); K. Bamba, S. Capozziello, S. Nojiri and S. D. Odintsov, Astrophys. Space Sci. 342 , 155 (2012).
- [5] A. Sen, Phys. Scripta 117 , 70-75 (2005).
- [6] G. W. Gibbons, Phys. Lett. B 537 , 1 (2002).
- [7] M. R. Garousi, Nucl. Phys. B 584 , 284 (2000).
- [8] A. Sen, Mod. Phys. Lett. A 17 , 1797 (2002).
- [9] A. Unzicker and T. Case, [arXiv:physics/0503046].
- [10] K. Hayashi and T. Shirafuji, Phys. Rev. D 19 , 3524 (1979) [ Addendum-ibid. D 24 , 3312 (1982)].
- [11] R. Aldrovandi and J. G. Pereira, Teleparallel Gravity: An Introduction , Springer, Dordrecht (2013).
- [12] J. W. Maluf, Annalen Phys. 525 , 339 (2013).
- [13] R. Ferraro and F. Fiorini, Phys. Rev. D 75 , 084031 (2007).
- [14] G. R. Bengochea and R. Ferraro, Phys. Rev. D 79 , 124019 (2009).
- [15] E. V. Linder, Phys. Rev. D 81 , 127301 (2010).
- [16] A. De Felice and S. Tsujikawa, Living Rev. Rel. 13 , 3 (2010).
- [17] S. i. Nojiri and S. D. Odintsov, Phys. Rept. 505 , 59 (2011).
- [18] B. L. Spokoiny, Phys. Lett. B 147 , 39 (1984); F. Perrotta, C. Baccigalupi and S. Matarrese, Phys. Rev. D 61 , 023507 (1999).
- [19] V. Faraoni, Phys. Rev. D 62 , 023504 (2000); E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70 , 043539 (2004).
- [20] O. Hrycyna and M. Szydlowski, JCAP 0904 , 026 (2009); O. Hrycyna and M. Szydlowski, Phys. Rev. D 76 , 123510 (2007); R. C. de Souza and G. M. Kremer, Class. Quant. Grav. 26 , 135008 (2009); A. A. Sen and N. Chandrachani Devi, Gen. Rel. Grav. 42 , 821 (2010).
- [21] C. -Q. Geng, C. -C. Lee, E. N. Saridakis, Y. -P. Wu, Phys. Lett. B 704 , 384-387 (2011).
- [22] C. -Q. Geng, C. -C. Lee and E. N. Saridakis, JCAP 1201 , 002 (2012).
- [23] C. Xu, E. N. Saridakis and G. Leon, [arXiv:1202.3781 [gr-qc]].

- [24] H. Wei, Phys. Lett. B 712 , 430 (2012).
- [25] G. Otalora, JCAP 1307 , 044 (2013).
- [26] A. Banijamali and B. Fazlpour, Astrophys. Space Scie. 342 , 229 (2012).
- [27] G. Otalora, Phys. Rev. D 88 , 063505 (2013).
- [28] S. H. Pereira, A. Pinho S. S. and J. M. Hoff da Silva, [arXiv:1402.6723 [gr-qc]].
- [29] S. Unnikrishnan, Phys. Rev. D 78 , 063007 (2008); R. Yang and J. Qi, Eur. Phys. J. C 72 , 2095 (2012).
- [30] I. Quiros, T. Gonzalez, D. Gonzalez, Y. Napoles, R. Garcia-Salcedo and C. Moreno, Class. Quant. Grav. 27 , 215021 (2010).
- [31] E. J. Copeland, A. R. Liddle and D. Wands, Phys. Rev. D 57 , 4686 (1998).
- [32] L. Amendola, Phys. Rev. D 60 , 043501 (1999); Z. K. Guo, R. G. Cai and Y. Z. Zhang, JCAP 0505 , 002 (2005).
- [33] H. Wei, Nucl. Phys. B 845 , 381 (2011).
- [34] L. P. Chimento, A. S. Jakubi, D. Pavon and W. Zimdahl, Phys. Rev. D 67 , 083513 (2003).
- [35] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15 , 1753 (2006); J. Frieman, M. Turner and D. Huterer, Ann. Rev. Astron. Astrophys. 46 , 385 (2008); S. Tsujikawa, arXiv:1004.1493 [astro-ph.CO].
Behnaz Fazlpour Ali Banijamali
Dynamics of Generalized Tachyon Field in Teleparallel Gravity
We study dynamics of generalized tachyon scalar field in the framework of
teleparallel gravity. This model is an extension of tachyonic teleparallel dark
energy model which has been proposed in [26]. In contrast with tachyonic
teleparallel dark energy model that has no scaling attractors, here we find
some scaling attractors which means that the cosmological coincidence problem
can be alleviated. Scaling attractors present for both interacting and
non-interacting dark energy, dark matter cases.