Phase Space description of Nonlocal Teleparallel Gravity

We study cosmological solutions in nonlocal teleparallel gravity or $f(T)$ theory, where $T$ is the torsion scalar in teleparallel gravity. This is a natural extenstion of the usual teleparallel gravity with nonlocal terms. In this work the phase space portrait proposed to describe the dynamics of an arbitrary flat, homogeneous cosmological background with a number of matter contents, both in early and late time epochs. The aim was to convert the system of the equations of the motion to a first order autonomous dynamical system and to find fixed points and attractors using numerical codes. For this purpose, firstly we derive effective forms of cosmological field equations describing the whole cosmic evolution history in a homogeneous and isotropic cosmological background and construct the autonomous system of the first order dynamical equations. In addition, we investigate the local stability in the dynamical systems called"the stable/unstable manifold"by introducing a specific form of the interaction between matter, dark energy, radiation and a scalar field. Furthermore, we explore the exact solutions of the cosmological equations in the case of de Sitter spacetime. In particular, we examine the role of an auxiliary function called"gauge"$\eta$ in the formation of such cosmological solutions and show whether the de Sitter solutions can exist or not. Moreover, we study the stability issue of the de Sitter solutions both in vacuum and non-vacuum spacetimes. It is demonstrated that for nonlocal $f(T)$ gravity, the stable de Sitter solutions can be produced even in vacuum spacetime.


I. INTRODUCTION
It has been supported that in addition to the inflationary stage [1] in the early universe, currently the expansion of the universe is also accelerating by various cosmological observations including Type Ia Supernovae [2], cosmic microwave background (CMB) radiation [3], large scale structure [4], baryon acoustic oscillations (BAO) [5] as well as weak lensing [6]. We have two representative explanations for such a late-time cosmic acceleration. One approach is to introduce "dark energy" (DE) in the context of general relativity. The other approach is to consider the modification of gravity on the large scale (for reviews on not only DE problem but also modified gravity theories, see, for example, [7]).
There is a possible candidate for a theory of gravitation alternative to general relativity, namely, teleparallel gravity, which is described by using the Weitzenböck connection [8]. In teleparallel gravity, there exists torsion. This is opposite to the case of general relativity, in which the Levi-Civita connection is used. The torsion scalar T represents the Lagrangian density of teleparallel gravity. It can be extended to a function of T , that is, f (T ) gravity (for a recent review, see, for instance, [9]). This idea is similar to that of f (R) gravity [10], where R is the scalar curvature. Inflation in the early universe [11] and the late-time cosmic acceleration [12] can be realized in f (T ) gravity. Various cosmological and astrophysical considerations in f (T ) gravity have widely been executed [13]. It is known that in f (T ) gravity, the local Lorentz invariance is broken [14], and the relevant investigations on this point have been discussed [15].
On the other hand, in Ref. [16], there has been considered a way of modifying gravitation, the so-called nonlocal gravity, which comes from quantum effects. Furthermore, in order to unify inflation in the early universe and the late-time accelerated expansion of the universe, non-local gravity has been modified by adding an f (R) term in Ref. [17]. In addition, a possible solution for the cosmological constant problem through the nonlocal property of gravitation [18] has been proposed. Moreover, a physical mechanism by which a cosmological constant is screened in the framework of nonlocal gravity has been investigated [19][20][21]. It has also been indicated that in nonlocal gravity, there is the issue of ghosts [19]. Various aspects of nonlocal gravity have widely been explored [22] (for a recent review on nonlocal gravity, see, e.g. [42]). It is worth noting that nonlocal terms ✷T was first used in the framework of modified teleparallel gravity in Ref. [43] Furthermore , the nonlocal deformations of teleparallel gravity have been analyzed in Refs. [44], [45]. This theory is called nonlocal f (T ) gravity, which can be considered as an extension of nonlocal general relativity to the Weitzenböck spacetime. It has been discussed that there is a possibility to distinguish teleparallel gravity from general relativity by future experiments detecting nonlocal effects. In this paper, we investigate exact cosmological solutions in nonlocal f (T ) gravity. We analyze the autonomous system of the first order dynamical equations by deriving effective forms of cosmological field equations in a homogeneous and isotropic cosmological background, describing the whole evolution history of the universe. Moreover, we propose a specific form of the interaction between matter, dark energy, radiation and a scalar field and examine the local stability in the dynamical systems, which is called "the stable/unstable manifold". As a result, it is demonstrated that the system has a stable attractor. Furthermore, we study exact solutions of the cosmological equations in the case of de Sitter spacetime. Particularly, we explore the role of an auxiliary function called "gauge" η in the formation of such cosmological solutions and show whether the de Sitter solutions can exist or not in this scenario. In addition, we consider the stability problem of the de Sitter solutions both in vacuum and non-vacuum spacetimes and find that even in vacuum spacetime, the stable de Sitter solutions can be produced in the framework of nonlocal f (T ) gravity.
The organization of the paper is the following. In Sec. 2, we explain the framework of nonlocal f (T ) gravity. In Sec. 3, we explore the cosmological background and effective field equations in nonlocal f (T ) gravity. In Sec. 4, the interaction term and phase portrait are analyzed. In Sec. 5, the de Sitter solution is derived and its stability is examined in the following Sec. 6. Finally, conclusions are provided in Sec. 7.

II. FORMAL FRAMEWORK OF NONLOCAL f (T ) GRAVITY
Let us develop the formalism of nonlocal modified gravity with torsion T in a same manner as the nonlocal f (R) gravity is developed [17]. We suppose that the possible action for gravity with matter contents is given in terms of classical gauge invariant action as follows: (1) where κ = 8πG, G is Newtonian gravitational constant, L m is matter Lagrangian. To describe the geometry of spacetime in teleparallel gravity, it is com-monly used the tetrads formalism where the metric can be written in an orthogonal frame e µ a , in a such manner that g µν = e µ a e ν b η ab , where Greek alphabets run from µ, ν = 0...3, the flat Minkowski metric is denoted by η ab . Note that e µ a e a ν = δ µ ν and ✷ −1 is considered as an integral over the entirely spacetime manifold. The local operator ✷ = ∇ µ ∇ µ is called d'Alembert operator defined as ✷ = e −1 ∂ α (e∂ α ) here e = det(e µ a ) = −det(g αβ ) and T is torsion scalar and it is defined in a same form as f (T ) gravity.
It is always possible to reduce nonlocal theories to scalar-tensor equivalent theories and it is easy to do that for our model given in (1) using two auxiliary (non ghost) The reason that those fields are considered as non ghost is that , the norm of them defined as ||φ|| = Σ |φ| 2 ed 4 x is always positive definite and never becomes complex as long as the metric and its torsion T remains real numbers. As long as we work in Riemanninan manifolds this condition will be hold and we can safely use them as an appropriate set of auxiliary fields.
The new form for the reduced action is written as follows: Note that in (3) the action function f (φ) is supposed to have any desired form. Formally if we take the case: ξ = 1 and f (φ) = 2 in action of theory given by Eq. (2), then the action is reduced to the one which is equivalently of Teleparallel Gravity for T = 0. Classical tests for GR prove a very good agreement with observations. As a result it is very important to know whether this nonolocal teleparallel gravity has GR limit or not. At the level of action we already demonstrate it. By an enough good choosing of the function f (φ) we can recover GR as a limiting case. The form of equations of motion is presented in [44] : ✷φ − T = 0.
here the tensor E β α is defined through the variation of e as follows δe = eE β α e a β , Θ β a is the energy-momentum tensor of matter contents defined by Θ β a = e −1 δ(eLm) δe a β and ✷ ≡ e −1 ∂ µ (e∂ µ ). In Refs.
[44] , the authors investigated cosmological data analysis on a suitable chosen function f (φ) = A exp(nφ) and later in Ref. [45] , using Noether symmetry approach. In our paper we will fix f (φ) in another simple/adequate form in next section.

III. COSMOLOGICAL BACKGROUND AND EFFECTIVE FIELD EQUATIONS
The aim of this section is to write equations of motion for a cosmological background in the presence of matter fields in an effective form. Let us suppose that the non singular, physical metric of spacetime is given in the form of a Friedman-Lemaitre-Robertson-Walker (FLRW) metric given by ds 2 = dt 2 −a(t) 2 (dx b dx b ), where b = 1, 2, 3 is spatial coordinate and a(t) is scale factor and measures expansion of the whole cosmological Universe as well as its acceleration/deceleration phase. The corresponding suitable, diagonal tetrads basis is given by e a µ = diag 1, a(t), a(t), a(t) . The set of FLRW equations and the equations for the scalar fields are written as follow: The matter energy-momentum tensor is given in terms of a diagonal tensors for matter , dark energy, radiation as follows: where e a µ e µ a = δ a b is unit matrix. The matter budget of our model is drak matter density ρ m , radiation field density ρ r and scalar field density ρ φ . In order to preserve the acceleration expansion and the existence of late time de Sitter cosmology we inserted a non zero cosmological constant Λ with energy density ρ Λ . As an attempt to keep simplicity we assume that all matter contents are given in barotropic forms, where we define the equation of state (EoS) parameter w a for each fluid component, namely matter, radiation and cosmological constant and as a result for any component of matter field we have a linear EoS , i.e, p a = w a ρ a . Here the roman index a refers to different matter contents. Namely we denote it by a = {m, Λ, r} where m is for matter, Λ is for DE and r is for radiation field. Note that neither φ nor ξ are considered as the DE. The reason is that both fields play the role of auxiliary fields. We can't make guarantee that whether the fields φ, ξ will be ghost or not. Actually the appearance of ghost scalar fields in the non local theories for gravity is an important issue and should be addressed adequately. For example in the non local extensions of the GR, when the action is corrected by nonlocal terms 1 ✷ R or higher order terms, one must count the number of the degrees of freedom of the localized form of the Lagrangian. Additionally, one needs to check the equivalence between local and nonlocal representations of theory both at the action level and equations of motion levels. It is possible to make a categorization based on the first form of the auxiliary fields. Based on this classification we can find the number of algebraic constraints which they will limit our ability to write the local or nonlocal representations of theory. Although in nonlocal extensions of the GR, as long as we have a linear term we can ascertain the equivalence between frames. However, with higher order terms this equivalence is broken. That means in a general nonlocal GR when we have only curvature terms our theory may suffer from ghosts. In nonlocal extensions of the teleparallel gravity (TEGR), we can deduce the same as long as we make the theory using a nonlocal action made by linear scalar torsion T , and, hence, no ghost will appear. That is because it was proved that the Einstein-Hilbert action GR is dynamically equivalent to the TEGR at level of action as well as the equations of motion [46]. In our study with the nonlocal term which we will opt in next paragraph the model will be ghost free. Although probably the scalar field φ will not be a ghost but still we do not have any strong reason to keep it as the only sector for the acceleration expansion in our model. For this reason we also keep the cosmological constant Λ and its energy density ρ Λ . The first challenge is to choose one suitable form for f (φ). Note that ✷ξ = −f ′ (φ)✷φ. It is illustrative to expand and write this equation in the following equivalent form : Here, we suppose that f ′′ (φ) = 0. Thus one suitable class of models is: Note that for A = 0 , B = 2 the results reduce to the GR as a limiting case. Note that now, ✷(ξ + f (φ)) = 0, and we have a freedom to take ξ + f (φ) = Ψ, where Ψ is a harmonic function over R 4 . A possible option is to consider Ψ = 2Aφ, consequently the set of eqs. (6,7) are simply written in the following forms: Here parameter B measures the difference between TEGR and nonlocal theory respectively. Because we will study time evolution of the energy densities, it is adequate to rewrite cosmological equation , presented in Eq. (12) in the following forms, in the above equation, we have defined Note that the other density functions can not be written explicitly in terms of the scale factor a or scalar fields, till the time which we will present continuity equations for all the matter components. In our scenario we assumed that different matter contents interact with each other through some interaction forms which will present in next section.

A. Hartman-Grobman linearizion theorem
To investigate the phase space analysis it is needed to reduce the system of equations to an autonomous system of first order differential equations in the form d X dN = f ( X), where N plays the role of time and X is a vector field with density functions as components.In many works, investigation of various aspects of dynamical systems in cosmology of modified gravity is discussed [41].
The Hartman-Grobman linearizion theorem provides a powerful technique to study the local stability and the portrait of the phase space, when we have a set of hyperbolic fixed points . Let X(t) ∈ R n be a non trivial solution to the following system of first order differential equations, called flow, here g( X) is a locally Lipschitz, one-to one continuous map g : R n → R n . Let X * denotes the location of the fixed points of the dynamical system (18), and the corresponding Jacobian matrix, which we denote as J (g), is equal to, In order to have stable fixed points for system (18) it is enough to set all eigenvalues of the Jacobian matrix so that λ i satisfies Re(λ i ) = 0. The Hartman theorem predicts the existence of a homeomorphism F : It is proved that (20) is a topologically conjugate flow to the one system given in Eq. (18).

B. Building the cosmological autonomous system of equations
Now we study a model of interacting matter contents, where the continuity equation for each energy density ρ a is given by the following form: where a = {m, Λ, r, φ} and Γ a is the interaction function given by the general form Γ a = Γ a (Ω m , Ω Λ , Ω r , Ω φ ) and it satisfies Σ 4 a=1 Γ a = 0. In f (T ) gravity, such interacting models are widely studied in the literatures, namely [48]- [51]. In Ref. [52] authors showed that the total gravitational energy is transferred from dark matter ρ m to dark energy ρ Λ , and the cosmological coincidence problem in the Lambda-Cold Dark Matter (ΛCDM) model is slightly assuaged.
In comparison to matter, DE and radiation energy densities, let us define an auxiliary scalar energy density as It is important to mention here that the auxiliary field φ is not a physical field. Consequently the kinetic term could be treated as tachyonic field as well as pressureless dust matter. In this paper we consider w φ as a free parameter to be adjusted using observational data. In this case we can write the following equation for the ratio between pressure and density , called effective EoS equation, as follows: where we prescribed the form of f (φ) as it is given in Eq. (11) and we supposed that w m = w φ = 0, w r = 1 3 , w Λ ∈ (−1, − 1 3 ). It is easy to rewrite (21) using the definition of in the following set of first order differential equations where we used (22) in it, Recall a = {m, Λ, r, φ} and we use slow-roll variable N = log( a a0 ) = − ln(1 + z) (the derivatives will be taken with respect to N ) and z is redshift. This is an autonomous system and should be analyzed in the vicinity of critical points where dΩa dN | c = 0 using techniques developed in Sec. (III A).
In terms of the variables (23) the Friedmann equation (12) becomes the restriction : Note that due to the interaction term in the model, the density parameters Ω m , Ω Λ , Ω r , Ω φ should be interpreted very strictly as effective density parameters. We mention here that the above constraint guaranteed the existence of possible cosmological attractors, because actually the shape of density functions remains typically the same and the full 4-dimensional configuration space constructed using density functions defines a shape invariant manifold and it defines the attractor solution in the dynamical system.
The effective EoS for system is defined

IV. INTERACTION TERM AND PHASE PORTRAIT
The general linear dependent model for interaction could be in the following form: There are some criticisms about interacting models of DE, however the thermal properties of this model in var-ious gravities have been discussed in the literature [53]. Furthermore, in Ref. [54], the authors proposed a systematic scheme to construct the interaction form Γ a in a self consistent manner both in the perturbed form and in the background. They proved that in the perturbation formalism , there are possible ways to break the degeneracy between the interaction, DE EoS and DM abundance.
With this choice, the system of equations (24) are written in the following form: Let us study a class of these models where the DE interacts with both matter Ω m and scalar field components. . Based on our former notation given in (27), our interaction model is parametrized as follows: The autonomous system of first order differential equations for density functions are written in the following forms, These equations are related to the dynamics and the interaction form, characterizing the main properties of our model.

A. The critical(fixed) points
We stress here that the high dimensionality of the phase space, where the system is described using dynamical systems presented in previous section, restricts us to have an effective graphical description of the phase space, and thus we will focus our investigations only on the analytical results.
To make the dynamical analysis we first need to find the critical(fixed) points of the system by setting the left hand side of equations (30)- (33) to zero. Then we use the Hartman theorem to find the type and stability of each point [55].
The location of the fixed points P = (Ω m , Ω Λ , Ω r , Ω φ ) and their corresponding eigenvalues of the dynamical system are in the following table, where the stability of the fixed points is determined by evaluating the eigenvalues of the Jacobian matrix associated with the system.
The corresponding Jacobian matrix, which we denote as J (g), is equal to, The Eigenvalues and their stability for each point are written as following: • Stability for point A: The enough and sufficient condition to have A 1 as a stable fixed point for system is that all eigenvalues of Jacobian matrix λ i must satisfy Re(λ i ) = 0, i.e., From these it is found that the stability occurs at: The corresponding effective EoS behaves like w ef f = w φ . Depending on the w φ , EoS evolves from larger than −1 to less than −1, that is, it crosses the phantom divide line of w ef f = −1.
• Stability for point B This is unstable critical point and the corresponding effective EoS, w ef f = w Λ is always larger than −1 and it crosses the phantom divide line of when w Λ = −1.
• Stability for point C: Stability condition is −1 < w φ , 2b − 1 < w φ . Consequently C can be stable conditionally. The corresponding effective EoS is w ef f = −1 is located at the crosses the phantom divide line.
• Stability for point D: We obviously conclude that is unstable. The corresponding effective EoS,is given w ef f = 0 and is always larger than −1 and it can not crosses the phantom divide line.
In Fig. 1, we plot the time evolution of the density functions for b = 0.5, 0.7, 0.9, w Λ = −1/3, w φ = 0, where the horizontal axis shows log(1 + z) and the vertical axis does the value of the density functions. We can observe that the density functions of matter Ω m and the cosmological constant Ω Λ increase in time, whereas the density functions of radiation Ω r and Ω φ decrease in time. For low redshift values, 0 < z < 0.1, the densities Ω φ , Ω r are monotonically increasing functions, but matter and cosmological constant density functions decrease. At the present redshift z ∼ 0, Ω φ ∼ Ω r are negligible in comparison to the matter and cosmological constant densities. This confirms our remarkable observation about the scalar field φ that it cannot play the role of DE. So, the density of the scalar field is almost negligible at the present time. At distinct values of the redshift shown as z * , densities of matte, radiation and scalar field become equal, i.e. Ω m ∼ Ω r ∼ Ω φ . This occurs at z * ≈ 0.1. Furthermore, there is an era when z † ∼ 0.4 in which Ω m ∼ Ω λ , shows another equilibrium among matter and cosmological constant. These behaviors of the density functions are compatible with the observations.

A. Observation of a type of deceleration to acceleration phase transition
The deceleration parameter q is defined as If the expansion of the Universe is decelerating, q > −1, while if the cosmic expansion is accelerating, q < −1. In our model using (22) we obtain We find numerical solutions by using h = 0.7127 +0.013 −0.015 km/s/Mpc, Ω Λ = 0.7018 +0.018 −0.02 , and Ω m0 = 0.2981 +0.02 −0.018 , with χ 2 min = 707.4, H 0 = 73.24 ± 1.74 km/s/Mpc. In Fig. 2, we depict the time evolution of the deceleration parameter q for b = 0.5, 0.7, 0.9. Here, the horizontal axis shows log(1 + z) and the vertical axis shows the value of q. All of the curves meet. From Fig. 2, it is found that in the past for lower values of redshift, the value of q evolved from larger than −1 to less than −1, namely, the expansion phase of the Universe changed from the deceleration to the acceleration. This is consistent with the observations.

B. Effective Equation of State (EoS) of the Universe
Effective EoS of the Universe was defined in Eq. (26). We note that in the DE dominated stage, the value of the EoS of DE can be regarded as the effective EoS of the Universe w ef f . In Fig. 3, we show the time evolution of the effective EoS w ef f for b = 0.5, 0.7, 0.9 , where the horizontal axis shows log(1 + z) and the vertical axis shows the value of w ef f . From Fig. 3, it is found that for low redshift values , the value of w ef f became less than − 1 3 and therefore the cosmic expansion phase of the Universe changed from the deceleration to the acceleration. It is also seen that in our model, the value of w ef f evolves from larger than −1 to less than −1; that is, it crosses the phantom divide line of w ef f = −1. The value of the Eos parameter at the present redshift is around w ef f ∼ −1.6 < −1 shows an acceleration expansion beyond the phantom line.

VI. DE SITTER SOLUTION
Cosmological models usually have de Sitter (dS) solution where the Hubble parameter is constant (or almost constant in inflationary scenarios) H = H 0 as trivial solution. In GR such solution (used in inflationary mechanism as well as late time cosmology) becomes accessible when the dominant energy density ρ ≈ ρ 0 , which means that to have dS we need matter fields with very slowly varying energy density. It is not possible to find dS solution as an empty space solution in GR . But in modified gravity because of the geometrical terms (curvature R or torsion T ) it will be possible to find dS as an (almost) exact solution for field equations. In forthcoming sections, we look for dS solution both in empty and matter contents cases in model defined by Eqs. (6,6).
Let us firstly perform a little investigation on the equations of motion. As a result of continuity eq, we have two additional Klein-Gordon like dissipative eqs. (21) in non interacting case, we can find equations of motion for φ, ξ. If we suppose that in (11), B = 0, A = −1, they are written as following: ξ + 3Hξ + 6H 2 = 0.
Note that always with f (φ) = −φ, we have ∇ µ ∇ µ (φ − ξ) = 0, and we have a "gauge freedom" to write fields φ, ξ as follows: An exact solution for η, in FLRW background is given by, Note thatη(t) ∼ ρ m (t). In this case we can interpretη(t) as cold dark matter density. It is possible to take η 0 = 0 or η 0 = 0. We will study both cases in next subsections.
From first equation we find φ(t) = 4H 0 t + φ 0 . If we substitute it in second equation we obtain H 0 = 0. In f (T ) gravity stability for Einstein Universe is well studied in [64]. This is just Einstein static Universe and it proves that no dS solution exist.
Let us check whether this solution is stable or not. We make perturbation around the solution given by (H, φ) = (0, φ 0 ). The equation is given by the following: substituting the zeroth order solution we find Exact solution for perturbation function is When t → ∞, we clearly observe that perturbation is growing up linearly and consequently the system becomes unstable under infinitesimal field and background perturbations.

Case η(t) = 0
: When η = η0 a(t) 3 dt, we haveφ −ξ = η0 a(t) 3 . Thus the system (38,39) and (22) is reduced to the following system: Note that in dS phase, a(t) = a 0 e H0t . An exact solution for fields pair φ, ξ is given as follows: Here C 1 , C 2 are arbitrary integration constants and ∆ = 16H 0 a 6 0 + a 3 0 H 2 0 e −6H0t . The solutions given in (50,51) are exact solutions for dS phase of our model under study. Let us study its stability under perturbations in time representation of fields and backgrounds. Later in Sec. (VII A) another equivalent analysis using slow-roll coordinate N will be introduced.

B. Case of matter contents
Now we study exact solutions in dS phase when ρ a = 0. As a general case we consider the model given in (6,7) for general densities ρ m = ρ 0 m a(t) −3 .
Exact solutions provide that H 0 = 0 is the only possible solution. Thus similar to the empty case, still we just have Einstein static Universe. If We perturb the system around the solutions given above, we have the following system of equations: where φ 0 = C is a constant. The system is clearly asymptotically unstable, consequently no stable static Einstein Universe exists. When η = η0 a(t) 3 dt, the system (38,39) and (22) is reduced to the following system: Exact solutions are given in terms of first integrals: Therefore we have dS solution. We also conclude here that the system behaves asymptotically, φ ∼ ξ ≈ 2H 0 t, consequently we can find the perurbated solutions when t → ∞ diverge and consequently system becomes unstable.

VII. STABILITY OF DE SITTER VIA HARTMAN-GROBMAN LINEARIZION THEOREM
Stability of dS solution plays an essential role in inflationary scenario to have thermalization phase. Because nonlocal theory supposed to be an alternative theories for inflation we will study the stability in this context.
In this case because of a(t) the system becomes non autonomous but still we can study the local stability in the vicinity of a critical point for t ≥ t s . Using the time coordinate N and by redefining ϕ ≡ φ ′ , α ≡ ξ ′ we have the following system of differential equations: ϕ ′ = −3(4 H 2 a 3 α ϕ + 2 H 2 a 3 ϕ 2 + κ ϕ ρ 0 m + 4 κ ρ 0 m ) H 2 a 3 α ϕ + 2 κ ρ 0 m (73) The corresponding linearized system near the unique physically accepted critical point {H = 0, α = −4, φ = −4} called a proper node, or a star point(Actually it defines Einstein static solution) has a triplet proper node degenerated eigen value λ = − 3 2 . The critical point is asymptotically stable , and it shows that system is asymptotically stable .

VIII. CONCLUSIONS
In this paper, we have considered the exact cosmological solutions in nonlocal f (T ) gravity, which can be regarded as an extension of nonlocal general relativity to the Weitzenböck spacetime. We have explored the autonomous system of the first order dynamical equations by deriving effective forms of cosmological field equations in a homogeneous and isotropic cosmological background to describe the whole evolution history of the universe. Furthermore, we have introduced a specific form of the interaction between matter, DE, radiation and a scalar field and analyzed the local stability in the dynamical systems, which is the so-called "the stable/unstable manifold". It has been found that the system has a stable(unstable) attractor solutions. In addition, we have investigated the exact solutions of the cosmological equations in the case of de Sitter spacetime. We have demonstrated whether the de Sitter solutions can exist or not in this scenario by examining the role of an auxiliary function called "gauge" η in the formation of such cosmological solutions. Moreover, we have studied the stability problem of the de Sitter solutions both in vacuum and non-vacuum spacetimes. It has been shown that for nonlocal f (T ) gravity, we can obtain the stable de Sitter solutions even in vacuum spacetime.