Coupled dark energy: a dynamical analysis with complex scalar field

The dynamical analysis for coupled dark energy with dark matter is presented, where a complex scalar field is taken into account and it is considered in the presence of a barothropic fluid. We consider three dark energy candidates: quintessence, phantom and tachyon. The critical points are found and their stabilities analyzed, leading to the three cosmological eras (radiation, matter and dark energy), for a generic potential. The results presented here enlarge the previous analyses found in the literature.


Introduction
Observations of Type IA Supernova indicates that the universe undergoes an accelerated expansion [1,2], which is dominant at present times (∼ 68%) [3]. Despite of ordinary matter, the remaining 27% is an unknown form of matter that interacts in principle only gravitationally, known as dark matter. The nature of the dark sector is still mysterious and it is one of the biggest challenges in the modern cosmology. The simplest dark energy candidate is the cosmological constant, whose equation of state w Λ = p Λ /ρ Λ = −1 is in agreement with the Planck results [3]. This attempt, however, suffers from the so-called cosmological constant problem, a huge discrepancy of 120 orders of magnitude between the theoretical prediction and the observed data.
Among a wide range of alternatives, a scalar field is a viable candidate to be used. Its usage includes the canonical scalar field, called "quintessence" [4][5][6][7][8], and the scalar field with the opposite-sign in the kinetic term, known as "phantom" [9,10]. Beyond the real scalar field case, a complex quintessence was also used in ref. [11] to account the acceleration of the universe. The U(1) symmetry associated with this complex scalar leads to a more sofisticated structure for the dark sector, and unless the stardard model of particle a rlandim@if.usp.br physics is a very special case in the nature, there is no reason (apart from simplicity) not to consider a richer physics of the dark sector 1 .
Another possibility of noncanonical scalar field is the tachyon and it comes from string theory. In the bosonic string, its ground state is the tachyon field, whereas in supersymmetric string theory a real tachyon is present in non-BPS Dp-branes, while the complex tachyon appears in a braneanti-brane system [13]. The tachyon potential has a minimum [14,15] and at this minimum the tachyon field behaves like a pressureless gas [16]. As soon as tachyon condensation in string theory had been proposed, tachyon was also regarded as a dark energy candidate [17][18][19].
Still regarding the dynamical dark energy, there exists the possibility of interaction between dark energy and dark matter [20,21], since their densities are comparable and, depending on the coupling used, the interaction can alleviate the coincidence problem [22,23]. This approach was applied to phantom and tachyon as well, in refs. [24][25][26][27][28][29].
When a scalar field is in the presence of a barothropic fluid (with equation of state w m = p m /ρ m ) the relevant evolution equations can be converted into an autonomous system and the asymptotic states of the cosmological models can be analyzed. Such approach was done for uncoupled dark energy (quintessence, tachyon field and phantom field for instance [30][31][32][33][34]) and coupled dark energy [21,28,[35][36][37][38][39][40][41]. Since the complex scalar field and the coupled dark energy are generalizations of the real field and the uncloupled case, respectively, we aim to study both possibilities together, in the light of the linear dynamical systems theory. Thus, we investigate in this paper the critical points that come from the evolution equations for the complex scalar field (quintessence, phantom and tachyon), considering the possibility of interaction between the two components of the dark sector. The dynamical equations are derived and the critical point are analyzed, to find out what kind of universe can come up. As we have said, this is a natural extension of the previous works [21,28,37,39]. As a result, we have found no new fixed points, however there are some crucial differences if compared with the case of real field.
The rest of the paper is organized as follows. In section 2 we present the basics of the interacting dark energy and the dynamical analysis theory, where we show the couplings we have used. In sections 3 and 4 we use the dynamical system theory for the canonical (and phantom) and the tachyon field, respectively. Within the respective section we show the critical points with their stabilities, and also the viable sequence of cosmological eras (radiation-matter-dark energy). Section 5 is reserved for conclusions. We use Planck units (h = c = 1 = M pl = 1) throughout the text.

Interacting dark energy and dynamical analysis
We consider that dark energy is described by a scalar field with energy density ρ φ and pressure p φ , and with an equation of state given by w φ = p φ /ρ φ . We assume that the scalar field is coupled with dark matter, in such a way that total energy-momentum is still conserved. In the flat Friedmann-Robertson-Walker background with a scale factor a, the continuity equations for both components and for radiation arė respectively, where H =ȧ/a is the Hubble rate, Q is the coupling and the dot is a derivative with respect to the cosmic time t. The indices m and r stand for matter and radiation, respectively. The case of Q > 0 corresponds to dark energy transformation into dark matter, while Q < 0 is the transformation in the opposite direction. In principle, the coupling can depend on several variables Q = Q(ρ m , ρ φ ,φ , H,t, . . .), so we assume for the canonical scalar (quintessence) and for the phantom field the coupling is Q = Qρ mφ [20,21], while for the tachyon field the coupling is Q = Qρ m ρ φφ /H [39], where Q is a positive constant. The case with negative Q is similar and we will not consider it here because the minus sign of the case Q < 0 can be absorved into theφ , instead of considering Q < 0.
To deal with the dynamics of the system, we will define dimensionless variables. The new variables are going to characterize a system of differential equations in the form where X is a column vector of dimensionless variables and the prime is the derivative with respect to log a, where we set the present scale factor a 0 to be one. The critical points X c are those ones that satify X = 0. In order to study stability of the fixed points, we consider linear perturbations U around them, thus X = X c +U. At the critical point the perturbations U satisfy the following equation where J is the Jacobian matrix. The eigenvalues of J determine if the critical points are stable (if all eigenvalues are negative), unstable (if all eigenvalues are positive) or saddle points (if at least one eigenvalue is positive and the others are negative, or vice-versa).

Quintessence and phantom dynamics
The complex scalar field Φ can be written as Φ = φ e iθ , where φ is the absolute value of the scalar field and θ is a phase. Both canonical and phantom fields are described by the Lagrangian where V (|Φ|) is the potential for the complex scalar and we consider it depends only on the absolute value of the scalar field φ ≡ |Φ|. We have ε = +1 for the canonical field (quintessence) and ε = −1 for the phantom field. For a homogeneous field φ ≡ φ (t) and θ ≡ θ (t), in an expanding universe with Friedmann-Robertson-Walker metric with scale factor a ≡ a(t), the equations of motion are εθ + 3H + 2φ where the prime denotes derivative with respect to φ . In the uncloupled case Eq. (8) gives rise to the effective potential in (7), which is d dφ ω 2 2a 6 1 φ 2 +V (φ ) [11], where ω is an integration constant interpreted as angular velocity. The first term in the brackets drives φ away from zero and the factor a −6 may make the term decreases very fast, provided that φ does not decrease faster than a −3/2 .
We assume the interaction between the scalar field with dark matter through the coupling Qρ mφ and it enters in the right-hand side of Eq. (7).
In the presence of matter and radiation, the Friedmann equations for the canonical (phantom) field are and the equation of state becomes We are now ready to proceed the dynamical analysis of the system.

Autonomous system
The dimensionless variables are defined as The dark energy density parameter is written in terms of these new variables as so that Eq. (9) can be written as where the matter and radiation density parameter are defined by (13) and (14) we have that x 1 , x 2 and y are restricted in the phase plane by the relation due to 0 ≤ Ω φ ≤ 1. Notice that if y = 0 the restriction (15) forbids the possibility of phantom field (ε = −1) because for this case Ω φ < 0. The equation of state w φ becomes which is a trivial extension of the real scalar field case. The total effective equation of state is with an accelerated expansion for w e f f < −1/3. The dynamical system for the variables x 1 , x 2 , x 3 , y, z and λ are where saddle or stable Table 2 Eigenvalues and stability of the fixed points for the quintessence (phantom) field.

Critical points
The fixed points of the system are obtained by setting The fixed points are shown in Table 1. Notice that x 3 and y cannot be negative and recall that Ω r = z 2 . Some of the fixed points do not exist for the phantom field because for those cases Ω φ is negative. The eingenvalues of the Jacobian matrix were found for each fixed point in Table 1. The results are shown in Table  2.
The eigenvalues µ 4e and µ 5e are At the first sight, one might think that the linear analysis would not give a complete description of the stability, because all fixed points, but (f), have at least one eigenvalue equals zero. However, as pointed out in [42], fixed points that have at least one positive and one negative eigenvalue are always unstable, and methods such as center manifold [42] should be used to analyze the stability of the critical points that can be stable [((g) and (h)]. Even so, for almost all fixed points, but (a) and (c), x 3 = 0, which means φ → ∞.  (16), the effective equation of state (17) and the density parameter for dark energy (13).
However, this limit implies that x 2 ∝ φθ /H → ∞ as well, provided that H is finite. This issue occured for the points (b) and (d) to (h), as can be seen in the Table 1 whose mathematical inconsistency indicates that the these critical points are not physically acceptable. Since all critical points are similar to those ones found in the literature [21,28,30] we reproduce the main results in the appendix 6 for the sake of completeness, which are valid for the case of real scalar field.

Summary
From the eight fixed points presented in the quintessence (phantom) case, only (a) and (c) are physically viable and they describe the sequence: radiation → matter. Both of them are unstable, however there does not exist a point that describe the dark-energy-dominated universe. Thus, the extra degree of freedom due to the phase θ spoils the physically acceptable fixed points that exist for the case of real scalar field, indicating that the dynamical system theory is not a good tool when one tries to analyze the complex quintessence (phantom).

Tachyon dynamics
The complex tachyon field Φ = φ e iθ , where φ is the absolute value of the tachyon field and θ is a phase, is described by the Born-Infeld Lagrangian where V (φ ) is the tachyon potential, which depends only on the absolute value of the scalar field φ ≡ |Φ|. For a homogeneous field φ ≡ φ (t) and θ ≡ θ (t), in an expanding universe with Friedmann-Robertson-Walker metric, the Lagrangian becomes where a ≡ a(t) is the scale factor. The equations of motion for φ and θ are respectivelÿ where the prime denotes time derivative with respect to φ . When the phase θ is zero, we recover the well-known equation of motion for the tachyon field. The Friedmann equations for the complex tachyon, in the presence of barothorpic fluids, are The equation of state for the tachyon field yields thus, the tachyon behavior is between the cosmological constant one (w φ = −1) and matter one (w φ = 0). The interaction between the tachyon and the dark matter is driven by the coupling Q = Qρ m ρ φφ /H, which in turn we consider that it modifies the right-hand side of Eq. (28). With this form of coupling, the time dependence of the coupling is implicit in the Hubble parameter H. We are now ready to proceed the dynamical analysis of the system. 5

Autonomous system
The dimensionless variables for the case of tachyon field are Sinceφ and θ are dimensionless variables, φ has dimension of time. The dark energy density parameter is written in terms of these new variables as so that Eq. (30) can be written as where the matter and radiation density parameter are defined by Ω i = ρ i /(3H 2 ), with i = m, r. From Eqs. (34) and (35) we have that x 1 , x 2 and y are restricted in the phase plane by the relation due to 0 ≤ Ω φ ≤ 1. In terms of these new variables the equation of state w φ is which is clearly a trivial extension for the complex scalar field. The total effective equation of state is with an accelerated expansion for w e f f < −1/3. The dynamical system for the variables x 1 , x 2 , x 3 , y, z and λ are
The fixed points y c , x f and y f are shown below The eingenvalues of the Jacobian matrix were found for each fixed point in Table 3. The results are shown in Table  4.
The points (a1)-(a4) correspond to a matter-dominated solution, since Ω m = 1 and w e f f = 0. They are saddle points because at least one eigenvalue has an opposite sign. The point (e) is in principle a dark-energy-dominated solution with Ω φ = 1 and w e f f = w φ = −1, whose existence is restrict to λ = 0. However a careful analysis shows that the Jacobian matrix for this critical point has zero eigenvector, thus it cannot be considered. Points (f1) and (f2) are also a dark-energy-dominated solution (Ω φ = 1) whose equation of state depends on λ , which in turn can be either constant or zero. The case with constant λ are shown in the Table 3 and an accelerated expansion occurs for λ 2 < 2/ √ 3. For λ = 0 we recover the point (e). The eigenvalues µ 2 , µ 4 and µ 5 of the fixed point (f1) and (f2) are always negative. For these points λ 2 y 2 c 3 ≤ 1, then the first eigenvalue is also negative if Q = 0, λ < 0 or √ 3Qλ y c < 3. Therefore, the point (f2) describes a dark-energy-dominated universe and can lead to a late-time accelerated universe if the requirement µ 1 < 0 is  (37), the effective equation of state (38) and the density parameter for dark energy (34). Table 4 Eigenvalues and stability of the fixed points for the tachyon field.
statisfied. On the other hand, (f1) is a saddle point. The effective equation of state depends only on λ , so the coupling Q only changes the property of the fixed point. The point (g) is also a saddle point with a matter-dominated solution, however, different from (a1)-(a4), the equation of state for the dark energy w φ is no longer zero, but depends on Q, leading to an universe with accelerated expansion for Q 2 < 2/3. For this point the coupling is restrict to values 0 ≤ Q 2 ≤ 1.
The last fixed points (h1) and (h2) are valid for x f = 0 4 , for Q = 0 and for constant λ , and its behavior depends on Q. In order to have x 2 f ≤ 1, we must have Q > 0 for the case with plus sign in x f (46), while we have Q < 0 for the minus sign case. We restrict our attention for the plus sign case. When Q → ∞, x f → 0 and y f → 1, in agreement with the 4 The case for x f = 0 is the fixed point (e). restriction (36). In addtition, as pointed out in ref. [39], the fixed points exists for some values of λ > 0 and Q, due to Eq. (36). Both fixed points have similar behaviour, however, (h1) is a saddle point, while (h2) is stable. Such difference is due to the eigenvalue µ 3 ( Table 4). The eigenvalues µ 4 and µ 5 are always negative because w e f f is between zero and minus one. The first eigenvalue is also negative because , since x f is always positive. Therefore, the point (h2) can lead to a late-time accelerated universe, depending on the value of λ and Q.
As in the case of quintessence and phantom, the fixed points that have x 3 = 0 [(a1), (a2), (a4), (b), (f1), (g) and (h1)] indicate that φ → ∞ and therefore x 2 ≡ φθ → ∞ as well. However this limit is in contradiction to what is pre-sented in Table 3 for x 2 showing that these seven critical points are not physically acceptable.
All fixed points reproduce the previous results in the literature [28,32,39,43] and they are generalizations of those analyses, with same stability behaviour for the critical points. This indicates that the degree of freedom due to the complex scalar has no effect on the stability and on the evolution of the system of equations, when compared with the case of real scalar field.

Summary
The critical points showed in the tachyonic case describe the three phases of the universe: the radiation-dominated era, the matter-dominated era, and the present dark-energydominated universe. The matter-dominated universe can be described by the saddle point (a3). There are two points that can represent the radiation-dominated era: (c) and (d). The two points are saddle, with the additional difference that the point (d) has an equation of state for dark energy equals to minus one. A tachyonic-dominated universe is described by the point point (f2) and (h2). The point (f2) can be stable only if the coupling is zero or λ < 0. The last fixed point (h2) is stable and can describe an accelerated universe depending on the value of λ and Q.
From all the critical points, the cosmological transition radiation → matter → dark energy is achieved considering the following sequence of fixed points: (c) or (d) → (a3) → (f2) [λ dependent] or (h2) [Q and λ dependent]. Although the sequence is viable, the form of the potential dictates whether the fixed points are allowed or not. Among several possibilities in the literature, the potential V (φ ) ∝ φ −n , for instance, leads to a dynamically changing λ (either if λ → 0 for 0 < n < 2, or λ → ∞ for n > 2) [19]. A dynamically changing λ is allowed for the fixed points (a3), (c) and (d).

Conclusions
In this paper we studied coupled dark energy using a complex scalar field, in the light of the dynamical system theory. There were analyzed three possibilities: quintessence, phantom and tachyon field. All three possibilities are known in the literature for the real field [21,28,37,39], and for uncoupled and complex quintessence field [33]. Thus, a natural question that arises is how a complex scalar field changes the previous results and if there are new fixed points due to the complex field. Although some equations for the dimensionless variables are trivial extensions of the real field case (e.g. the equation of state for the scalar field), the differential equations were generalized. All fixed points found here are in agreement with the previous results, with no new fixed points, however there are some crucial differences. For the quintessence and the phantom there is a contradiction between the fixed points x 2 and x 3 when the latter is zero. This situation occurs for almost all fixed points and the only two exceptions are unstable points that represent respectively the radiation and matter era, so the dark-energy-dominated universe is absent. Therefore the extra degree of freedom spoils the results known in the case of real scalar field. For the tachyon field all the critical points are also similar to the real field case, with same stabilities. Therefore, the extra degree of freedom due to the complex tachyon field plays no role on the stability of the critical points. Although the results presented here enlarge the previous results found in the literature, with the generalization of the equations of motion, the dynamical system theory does not provide further information in what is already known for the case of real scalar field, letting open the possibility of studing complex scalar fields by other ways of analysis.

Appendix
In this appendix we reproduce the results in the literature [21,28,30] regarding the fixed points presented in Table 1, for the sake of completeness.
The fixed point (a) is a saddle point which describes a matter-dominated universe, however it is valid only for Q = 0. The other possibility of matter-dominated universe with Q = 0 arises from the fixed point (b). This point is called "φmatter-dominated epoch" (φ MDE) [21] and it can be either unstable or a saddle point. However, due to Ω φ = 2Q 2 /3 1, the condition Q 2 1 should hold in order to the point be responsible for the matter era. Thus, µ 1 and µ 2 are negative, while µ 5 is always positive and µ 4 is positive for Q(λ + Q) > −3/2. Therefore (b) is a saddle point. The radiation-dominated universe is described by the critical points (c), (d) and (e), only for the quintessence field. The first two points are saddle, as it is easily seen in table 2, and the last one had its stability described numerically in [21]. However, both (d) and (e) are not suitable to describe the universe we live in, due to nucleosynthesis constraints [44,45]. The nucleosynthesis bound Ω BBN φ < 0.045 [46] implies Q 2 > 3.7 for the point (d) and λ 2 > 88.9 for the point (e). Thus, the requirement for the point (d) is not consistent with the condition of point (b) and the constraint on λ 2 does not allow a scalar field attractor, as we will se soon. Therefore, the only viable cosmological critical point for the radiation era is (c).
The point (f) is an unstable or saddle point and it does not describe an accelerated universe. The last possibility for the matter era is the point (g), with eingenvalues showed in table 2 and Eq. (25).
Since w e f f 0 for |λ | |Q|, the fixed point is either stable or stable spiral, hence the universe would not exit from the matter dominance.
On the other hand, the point (g) can lead to an accelerated universe, for the quintessence field case (ε = +1), provided that 3 < λ (λ + Q), because Ω φ ≤ 1, and Q > λ /2, from w e f f < −1/3. Regarding λ > 0, the two eigenvalues µ 1 and µ 2 are always negative and since Q > 3/λ − λ , the behaviour of µ 4d,5d depends on the second term in the square root of (25) From the condition 3 < λ (λ + Q) we have A < 0, and if A < 1 the fixed point is stable. Otherwise, i.e. A > 1, the critical point is a stable spiral. Thus, the value of the coupling dictates which behaviour the fixed point will have: stable for 3/λ − λ < Q < Q * or stable spiral for Q > Q * , where Q * is the solution of A = 1. However, even in the case where one can get Ω φ 0.7 [45,47], there are no allowed region in the (Q, λ ) plane corresponding to the transition from φ MDE to scaling attractor [21]. Thus, it is hard to gather the conditions for the point φ MDE and the point (g). For the case of the phantom field (ε = −1), the condition y 2 > 0 implies 2Q(Q + λ ) > 3. Hence, µ 4 < 0 and µ 5 > 0, and (g) is a saddle point. The last fixed point (h) leads to an accelerated universe provided that λ 2 < 2. With this condition, the eingenvalues µ 2 , µ 4 and µ 5 are always negative. The first eigenvalue µ 1 is also always negative for the phantom field, and it is for the quintessence field with the condition λ (λ + Q) < 3. Therefore, the point is stable if the previous conditions are satisfied.