A new approach on the stability analysis in ELKO cosmology

In this work it has been developed a new approach to study the stability of a system composed by an ELKO field interacting with dark matter, which could give some contribution in order to alleviate the cosmic coincidence problem. It is assumed that the potential which characterizes the ELKO field is not specified, but it is related to a constant parameter δ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta $$\end{document}. The strength of the interaction between matter and ELKO field is characterized by a constant parameter β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} and it is also assumed that both ELKO field and matter energy density are related to their pressures by equations of state parameters ωϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _\phi $$\end{document} and ωm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _m$$\end{document}, respectively. The system of equations is analyzed by a dynamical system approach. It has been found the conditions of stability between the parameters δ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta $$\end{document} and β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} in order to have stable fixed points for the system for different values of the equation of state parameters ωϕ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _\phi $$\end{document} and ωm\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega _m$$\end{document}, and the results are presented in form of tables. The possibility of decay of the ELKO field into dark matter or vice versa can be read directly from the tables, since the parameters δ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta $$\end{document} and β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document} satisfy some inequalities. It allows us to constrain the potential assuming that we have a stable system for different interactions terms between the ELKO field and dark matter. The cosmic coincidence problem can be alleviated for some specific relations between the parameters of the model.

Besides the usual scalar fields, some recent works have shown that there are some classes of Non-Standard Spinors with interesting properties which could be useful in order to describe both DM and DE. One of such spinors is called ELKO, from German Eigenspinoren des Ladungsk onjugationsoperators, which has the property to be an eigenspinor of charge conjugation and parity, possessing nonlocality [45][46][47][48][49][50][51][52][53][54][55][56][57][58][59][60][61][62]. It satisfies (CPT) 2 = −I and also has some other properties, like having spin 1/2 and a mass dimension 1, which makes it a good candidate to a particle with small probability of interacting with Standard Model particles, exactly as desired for the DM particles and maybe also to DE. The searching for ELKO dark matter at the CERN LHC has also been addressed recently [63,64].
The ELKO spinor field, although being a fermionic field, can be factored out in a time dependent scalar field multiplied by a four components spinor field. Its scalar part drives the evolution of the pressure and energy density associated to the ELKO field, and it is through these physical quantities that ELKO contributes to cosmological equations. Here we will consider the ELKO field as a candidate to DE interacting with a DM field. The recent concern in this kind of exotic non-standard field as a good candidate to DE or DM is due to the fact that the scalar part of ELKO spinors has a much richer structure than the standard scalar field when applied to cosmology [17][18][19][20][21][22][23][24][25][26][27]. The coupled system of equations involving the interaction of DE and DM plus the Friedmann equations are much more involved, even for a standard scalar field [26,27]. Such a system is much more involved when dealing with ELKO fields. Due to its complexity, dynamical system analysis has been applied in order to deal with ELKO field as a possible candidate to DM or DE. Although exact solutions are not found even within this method for most cases, obtaining stability points that turn the system well behaved provides useful information as regards some undetermined parameters of the system. In our specific case we are interested in some conditions to be satisfied by the potential of the ELKO spinor field in order to keep the system stable.
Maintaining system stability around some fixed points indicates that the involved system of equations has stable solutions for that specific region where the variables are evolving. In our case it can be interpreted as possible solutions where the ELKO field and DM can coexist, indicating a solution to the cosmic coincidence problem around these fixed points.
An interesting aspect concerning the ELKO field dynamics is that the choice of variables is an important question when one is looking for stable points of the dynamical system. Such a search for adequate variables depends in general on the structure of the physical quantities involved, like the energy density, pressure and Friedmann constraints. In some recent works [65][66][67], different choices of variables for interacting systems concerning the ELKO field have shown that there are no stable points in order to explain the cosmic coincidence problem. In [68] a new method of analysis was proposed, based on a constant parameter that leads to stable points under some conditions. In the present work the same method as [68] has been applied for a new set of dynamic variables and stable fixed points have been found for the system, including new interaction terms. This opens the possibility to alleviate the cosmological coincidence problem considering the ELKO field interacting with DM if the conditions found here are satisfied.
The new set of dynamical variables for the ELKO field used in this work is independent of the choice of potential. The potential dependence is set to a constant parameter δ. Besides that it is assumed that the pressure and energy density of the ELKO field satisfy an equation of state of the form p φ = ω φ ρ φ , where p φ and ρ φ are the ELKO field pressure and energy density, respectively. The dark matter content is described by an energy density ρ m that satisfies an equation of state of the type p m = ω m ρ m . There is also assumed an interaction between ELKO field and DM, characterized by a constant parameter β. The conditions for stability between δ and β have been studied, depending on the type (radiation, dust, vacuum or ultrarelativistic matter) of the thermodynamic equation of state parameters ω m and ω φ .
Let us finish mentioning that the most general ELKO theory applied to cosmology should include torsion [49,51,69,70]. This generalization is possible constructing an action where connections are no more symmetric at all. Such an antisymmetric part of the connection is defined as the torsion tensor, which is a rotation of a vector when it is parallel transported. Nevertheless problems with torsion have not been properly studied yet, although in some works we can see that torsion could play a cosmological constant role, accelerating the universe. It is also not well known if torsional effects appear to affect the dynamics only at the Planck scale.
The paper is organized as follows. In Sect. 2 we start presenting a general scenario including torsion and then we restrict ourselves to the torsionless case with the basic equations for a two fluid model, namely DM and ELKO field equations related to our cosmological applications in a spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) background. We present the pressure and energy density expressions, as well as the Friedmann equations and the conservation equations including the interaction between dark matter and ELKO field. We also define the variables concerning the dynamical system equations to be analyzed. Section 3 contains the main results. We study the stability of the dynamical system by imposing the restriction that the potential of the ELKO field is related to a constant parameter δ. This analysis allows us to study the conditions to be satisfied by δ and β for each equation of state parameters of DM and the ELKO field, in order to maintain stability. For each kind of interaction we present the results in the form of tables, specifying the equation of state parameter and conditions of stability to be satisfied. For DM and the ELKO field we restrict the equation of state parameters to vacuum, dust, radiation and ultrarelativistic matter. In Sect. 4 we conclude with some remarks.

ELKO cosmology as a dynamical system
We begin introducing the general Einstein-Cartan action: where the tilde denotes the presence of torsion terms in the covariant derivatives. We have also κ 2 ≡ 8π G with the normalization c = 1. The matter source (DM and ELKO) is present in the matter Lagrangian,L mat . In this work we have divided such a Lagrangian into two parts. One part concerns the DM (including baryonic matter),L DM , with an attractive behavior. The other part is playing the role of a DE fluid, here represented by the ELKO spinor LagrangianL λ : The covariant derivative acting on a spinor in the presence of torsion is where a = 1 4 abc [γ b γ c , γ b γ c ], and K abc is known as the contortion tensor, which represents the antisymmetric part of the Christoffel symbol. Then it is possible to find an analogous of the Einstein equation by variation of the action (1): and also an equation for the spin angular momentum tensor τ i j k , obtained by variation of the Einstein-Cartan action with respect to the contortion tensor K i j k , which is related to the torsion by T μν As well as in the ordinary Einstein equation, the right hand side of (4) represents the matter source, however, with a new form due to torsion. It is given by where σ i j corresponds to the energy momentum tensor related toL mat . In a flat FLRW metric, ds 2 = dt 2 −a(t) 2 (dx 2 +dy 2 +dz 2 ), where a(t) is the scale factor, the Christoffel symbols are i ti =ȧ/a and t ii =ȧa, where i = x, y, z, and whereρ andp are the energy density and the pressure of the matter sources, in our case composed by DM and ELKO field, which can be approximated by a perfect fluid on large scales. If the ELKO field density depends only on cosmological time, as it must obey cosmological principle, the allowed components of the torsion tensor are reduced to where h and f are functions carrying the torsion contributions. Finally, the Friedmann equations with torsion are given by with H =ȧ/a.
From the above general treatment, we see that even assuming homogeneity and isotropy, the torsion functions h(t) and f (t) yet carry non trivial contributions to the time component of the evolution equations. Thus, even for an ELKO field that behaves dynamically as a scalar, its contributions should be taken into account. When it is used in cosmology, it is important to emphasize that they are spinors, nevertheless with important consequences: the most important of them is the fact that a totally coupled ELKO field must have torsion, and with torsion additional terms are present in the Friedmann equations, as seen above.
However, in order to study the possibility of ELKO field interaction with DM, intending to alleviate the cosmic coincidence problem, we will restrict our study to the torsionless case. It is evident that the presence of torsion should lead to more general results.
When the torsion terms are dropped out, it is easy to see that we obtain the ordinary Friedmann equation. Before doing this, let us just make a simplification on the ELKO field. As has been done in recent works [52][53][54]57,65], we will restrict the ELKO spinor field to the form λ ≡ φ(t)ξ and where the subscripts m and φ indicate, respectively, the DM and ELKO field contributions to the energy density and pressure.
Thus, the FLRW equations in a spatially flat background without torsion are given by We will assume that the DM and ELKO field satisfy equations of state of the form p m = ω m ρ m and p φ = ω φ ρ φ , respectively. Due to the homogeneity of the field (∂ i φ = 0), the equation of motion that follows from (2) is substantially simplified tö where V ,φ ≡ dV /dφ. The pressure and energy density of the ELKO field are given by [54]: The continuity equations for DM and ELKO field are: where Q stands for a possible interaction term between the DM and ELKO field. If Q = 0, there is no interaction and both components evolve separately. If Q > 0, the ELKO field decays into DM, an interesting scenario as regards inflation, and if Q < 0, DM decays into the ELKO field (or DE). Exact solutions to the system of coupled equations (14)- (18) are not available. For this reason we make a stability analysis based on a dynamical system approach.
We define the new variables where we have redefined the potential asṼ = V + 3 8 H 2 φ 2 . The Friedmann equation (12) can be written as a constraint equation or in terms of the densities parameters, Equations (13), (17), and (18) can be written as a dynamical system of the form where λ =˙ṼṼ , Q 1 = κ 2 Q 6H 3 , and stands for the derivative with respect to N ≡ ln a, such that f =ḟ /H for any function f .
It is easy to see that the simple change of variables (19) should not be sufficient to rewrite the Friedmann equations (10) and (11) in a simple form with the torsion terms. At least two new variables would be necessary concerning f and h, and also the system of Eqs. (22)-(24) would be much more complicated in the torsion case.

Stability analysis
The stability analysis of the above set of dynamical equations consists of finding fixed pointsx,v, andȳ that make x , v , and y equal to zero. In the last section we have defined a three dimensional system according to our variable choice.
However, due to the Friedmann constraint (20) the system can be reduced to a two dimensional one.
Before we proceed, let us examine carefully the dynamical system (22)- (24). We see that, in addition to the dynamical variables x, v, and y we also have the factor λ 2H =Ṽ

2Ṽ
. The presence of the derivative shows that such term is also a dynamical variable, which should also be taken into account. However,Ṽ depends on the potential, but the potential is not specified in our analysis, thus we cannot deal with this new variable. In order to avoid this problem, we set the additional assumption related to the potential, − λ 2H ≡ δ, where δ is a constant parameter. Such an imposition just reflects our ignorance on the potential V (φ).
Now it is easy to see that the resulting dynamical system, for a given interaction Q 1 (x, v, y), is written in terms of the dynamical variables x, v, y, and the constants ω m , ω φ , and δ. In order to analyze the stability of this system around the fixed pointsx,v, andȳ we must study the system satisfying x = 0, v = 0, and y = 0. Notice that Eq. (24) is independent of the interaction Q 1 , thus the condition y = 0 can be achieved The first condition can be satisfied only ifṼ = 0 (see (19)), but from the definition of λ this leads to a divergent λ. So, we restrict ourselves to the second condition, namely where y = 0 even forȳ = 0. By imposing the above condition on the dynamical system (22)- (24) we are left with a 2 × 2 system, since y = 0 is always satisfied: where we have also used the Friedmann constraint (20). In order to study such a dynamical system it is worth to define its linearized matrix, with which one can determine the stability of a fixed point by just analyzing its determinant and trace. Such a mechanism is ensured by the Hartmann-Grobman theorem [71]. Thus, in the neighborhood of the fixed points we take infinitesimal displacements of the variables from the fixed points, x →x + δx and y →v + δv, so that where M is given by and All variables carry a bar over them to show that the matrix M is taken at the stable points which solve the system. There is a simple way to know if the system described by the matrix M is stable or not. It depends on the values of the determinant ( ) and also on the trace (τ ) of such matrix. When > 0, both eigenvalues have the same sign and if they are positive, the solution increases with time, indicating that the solutions diverge from the fixed point and consequently this point is classified as unstable. On the other hand, when both eigenvalues are negative, the solution goes to zero and the fixed point is stable. In order to know what kind of fixed point we are dealing with, it is necessary to check the value of the trace of matrix M. When τ > 0, it means that both eigenvalues are greater than zero, describing unstable points. However, when τ < 0 we see that they are negative and the point is stable. For the case where < 0 we see that both eigenvalues have opposite signs and then the fixed point is in fact a saddle point. Finally, when = 0, at least one of the eigenvalues is zero and consequently nothing can be said about the stability of system.
Let us return to the dynamical system. Together with Eqs. (26) and (27), the assumption − λ 2H = δ leads to the new constraint for the fixed pointv, according to (25): But it is easy to see that such a constraint already determines the value of the fixed pointv, since it depends only on the fixed parameters ω m , ω φ and δ. Thus, in order to also satisfy the dynamical equation (27), we have verified that this restriction tells us that Q 1 , which represents the interaction, could not assume an arbitrary value, since Eq. (27) would not be always solved for an arbitrary Q 1 . In other words, we have found that when Q 1 does not depend on x, not all fixed v that make Eq. (27) vanish also satisfy (34), except for some very specific relations among the parameters ω m , ω φ , and δ. Such a restriction on the interaction term Q 1 is not so strong, since the variable x is proportional toφ, which characterizes the time variation of the field φ, which is reasonable for an interacting theory.
In the following, the stability conditions will be analyzed for different interaction terms between DM and ELKO field. The interaction terms are characterized by a dimensionless coupling constant β. We will search for the stability conditions between the parameters δ and β for different equation of state parameters ω m and ω φ . Besides stability conditions characterized by negative eigenvalues of the matrix of perturbation M, we impose the additional reality condition on the parameters (19), namely we will imposex 2 > 0, v 2 > 0, andȳ 2 > 0. As particular cases, we will discuss the physical content concerning the present time, characterized byv 2 = m 0.315,x 2 +ȳ 2 = φ 0.685, and ω m = 0 according to recent observations based on the CDM model [72]. We will also analyze the inflationary phase, where we believe there is no matter contribution, corresponding tov 2 → 0.

Q 1 = 0 and Q 1 = β
For the case Q 1 = 0 there is no interaction between the DM and the ELKO field, thus they evolve independently. The fixed points that follow from the analysis of the system (22) to (24) and satisfy x = 0, v = 0, and y = 0 are given just byx = 1,v = 0, andȳ = 0, which does not represent a scaling solution, in the sense that does not admit a mixture of fluids. Besides that our model is valid just forȳ = 0 and v given by (34).
For the case Q 1 = β, a constant interaction term, we have scaling solutions of the formx = 0,v = 0, andȳ = 0, which could admit a mixture of the fluids, but the condition y = 0 shows that the potential part of the ELKO field is null according to (19), leading toṼ = 0, and as discussed earlier, this leads to a divergence in the λ term, but we have defined it as proportional to the constant δ, so we will discard such a kind of fixed point in our analysis. We are interested only in fixed points that satisfyx = 0,v = 0, andȳ = 0, which are scaling solutions and do not have null potential contributions.

Q 1 = βx 2
Such an interaction between DM and ELKO field corresponds to Q = β Hφ 2 , where β is a dimensionless parameter.
From the analysis of the system of Eqs. (26) and (27), the fixed points are given by [x,ȳ,v], with andv is given by (34), which is valid for all interactions.
The determinant and the trace τ of the matrix of the linearized system (28) to (33) are given by The stability of the fixed points, namely > 0 and τ < 0, is related with the values of δ, β, ω m , and ω φ . Table 1 presents the above conditions plus the reality conditions for the fixed points, namelyx 2 > 0,v 2 > 0, andȳ 2 > 0 for some specific values of the equation of state parameter for both DM (ω m ) and the ELKO field (ω φ ). Now let us analyze Table 1. We are interested in two different epochs, namely the inflation and the late time acceleration epochs. The first one corresponds to an universe without DM and totally filled with the ELKO field. This means v 2 = m = 0 and ω m = 0, which leads to δ = 3 2 (1 + ω φ ). By replacing intox andȳ we obtainx 2 = 0 andȳ 2 = 1, which shows that all the contributions should come from the potential partṼ and also we should haveφ = 0 from (19), but our interaction term Q ∼φ, thus we conclude that such a condition cannot be applied to inflation.
On the other hand, for the present time such variables are given byv 2 = m = 0.315 andx 2 +ȳ 2 = φ = 0.685 according to the CDM model, and besides that we must have ω m = 0. From (34) we conclude that: (i) δ 0.47 if ω φ = −1-from the corresponding cell in Table 1 (ω φ = −1, ω m = 0), we see that if β 0.47 the system is stable around the fixed points, and such a positive value of β corresponds to decay of ELKO field into DM; (ii) δ 1.84 if ω φ = 1/3, but there is no solution for β from Table 1 in this case; and (iii) δ 2.53 if ω φ = 1, which has no solution for β too. We conclude that, for the present time, the only possibility in order to have stable fixed points is β 0.47, leading to decay of ELKO field into DM for an equation of state parameter ω φ = −1, that is, the ELKO field behavior must be of vacuum type. Notice that other stability conditions are possible if the equation of state parameter of dark matter is of radiation or ultrarelativistic type.
The case ω m = ω φ has no physical meaning since both fluids have the same equation of state parameter; thus they are thermodynamically identical.
For this interaction we have Q = 1 3 κ 2 Hβρ mφ 2 . The fixed points are [x,ȳ,v], with andv is given by (34). From the linearized matrix one finds Table 2 shows the stability conditions for some specific values of ω m and ω φ .
For the inflationary epoch (v 2 = 0 and ω m = 0), we have the same condition for δ, namely δ = 3 2 (1 + ω φ ). From the corresponding cell in Table II it is easy to see that ω φ = −1 is the only possible condition of stability, which leads to δ = 0, but such a value of δ is not possible from Table II if ω m = 0.

Q
andv given by (34). From the linearized matrix one finds Table 3 shows the stability conditions for some types of ω m and ω φ .
For the inflationary epoch (v 2 = 0 and ω m = 0) we havē x 2 = 0 andȳ 2 = 1, which shows that the contribution comes only from the potential part. The kinetic part is null, but the interaction is proportional toφ 2 ; thus this interaction does not apply to inflation.
andv given by (34). The determinant and trace are given by Table 4 shows the stability conditions for some specific values of ω m and ω φ . For the inflationary epoch (v 2 = 0 and ω m = 0), we have the same condition for δ, namely δ = 3 2 (1 + ω φ ). Contrary to the previous cases, we see that ω φ = −1 is not a stable solution. For ω φ = 1/3 we have δ = 2, which is a stable solution corresponding to β = −1/2, representing the decay of DM into the ELKO field. But we have inferredv 2 = 0, thus there is no matter to decay at the inflation epoch. For ω φ = 1 we have δ = 3, but the stability condition requires δ < 3 from the corresponding cell in Table IV. We conclude that such an interaction does not present stable points for the inflation.
For the present time we havev 2 = 0.315 and ω m = 0. As in the previous case, we have: (i) δ 0.47 if ω φ = −1, which has no stable solution; (ii) δ 1.84 if ω φ = 1/3, and it is easy to see from Table IV that such a value of δ is possible for a negative value of β -for instance, if β = −1/2 the condition for δ is 1.69 δ 2; and (iii) δ 2.53 if ω φ = 1, and for this condition we also have stable solution for a negative β, as can be seen from Table IV. If β = −3/2 for instance, the condition on δ is 2.07 δ 3.0, which includes δ 2.53. Thus, contrary to the previous cases, if the ELKO equation of state parameter is of radiation or ultrarelativistic type, the system presents stable solutions if β is negative, which corresponds to the decay of DM into ELKO field. The present acceleration of the universe can be understood in this model as the decay of dark matter into ELKO particles. This is a very interest-ing scenario which also alleviates the cosmic coincidence problem.

Concluding remarks
In this work we have developed a new approach to study the stability of a system composed by an ELKO field interacting with DM, which could give some contribution in order to alleviate the cosmic coincidence problem. Since recent works [65][66][67][68] have not found stable points for such system for different dynamic variables and interactions terms, we are led to believe (without demonstration) that the system ELKO-DM does not allow stable points. Based on these results, we have supposed there to be an additional constraint to the dynamical system, namely that the potential of the ELKO field is related to a constant parameter δ; then we have analyzed the stability conditions for such a new system. We have also assumed that both the ELKO field and the dark matter energy density are related to the pressure by equations of state parameters ω φ and ω m , respectively. We have found different stability conditions relating the parameter δ and the interaction parameter β, which states if the decay is from DM to ELKO (β < 0) or from ELKO to DM (β > 0). Different values of ω φ and ω m corresponding to vacuum, dust, radiation, and ultrarelativistic equation of state parameter are presented in Tables 1, 2 For the first three tables, corresponding to the interactions Q 1 = βx 2 , Q 1 = βv 2 x 2 , and Q 1 = β(x 2 + y 2 )x 2 , the conditions for stable fixed points in order to satisfy the present stage of acceleration (with ω m = 0) are given by positive β and δ parameters and also require an equation of state parameter for the ELKO field of vacuum type (ω φ = −1). Positive values of β means the decay of ELKO field into DM particles. Such conditions could alleviate the cosmological coincidence problem. The inflationary phase cannot be driven for these interactions if we set ω m = 0 andv 2 = 0. Other pos- Table 4 Stability conditions for some equation of state parameters of DM and ELKO field, corresponding to the interaction Q 1 = β(v 2 − y 2 ). We have defined the following parameters: δ 1 = 3 4 − β − 1 4 9 + 16β 2 ;

Radiation
(ω m = 1/3) For the last interaction, namely Q 1 = β(v 2 − y 2 ), we have the opposite. The conditions for stable fixed points in order to satisfy the present stage of acceleration are given by negative β and positive δ parameters. The equation of state parameter for the ELKO field must be of radiation (ω φ = 1/3) or ultrarelativistic (ω φ = 1) type. The case ω φ = −1 does not present stable fixed points. Negative values of β mean the decay of DM particles into ELKO field. Such conditions also could alleviate the cosmological coincidence problem, and the equation of state parameter of ELKO field is not of exotic type. The inflationary phase cannot be driven for this interaction too. Other possibilities are allowed if the equation of state parameters of DM and ELKO field are of radiation or ultrarelativistic type. Another interesting aspect that follows from Table IV is that, contrary to the previous cases, there are no stable fixed points if ω φ < ω m .
In such an analysis the interaction Q must be proportional to the variable x, otherwise the relations among the parameters must be very restrictive. But such condition is not so strong, since the variable x is proportional toφ, which characterizes a time dependence of the field φ, and it is reasonable for an interacting theory. Notice that all the interactions studied are proportional toφ 2 . For all the interactions analyzed here there are conditions of stability in order to alleviate the cosmic coincidence problem. Such kinds of interactions and conditions on the parameters β and δ open possibilities for future searches concerning the interaction between DM and ELKO field for specific potentials satisfying the conditions presented in the tables.
Finally, a general theory that includes torsion in ELKO cosmology is of great interest, and in fact it may be the one responsible for the late time acceleration of the universe [69,70]. From a dynamical system approach the system of dynamical equation should be rewritten for new variables, certainly with some interesting consequences, but this is left for future.