Bulk viscosity, interaction and the viability of phantom solutions

We study the dynamics of a bulk viscosity model in the Eckart approach for a spatially flat Friedmann-Robertson-Walker (FRW) universe. We have included radiation and dark energy, assumed as perfect fluids, and dark matter treated as an imperfect fluid having bulk viscosity. We also introduce an interaction term between the dark matter and dark energy components. Considering that the bulk viscosity is proportional to the dark matter energy density and imposing a complete cosmological dynamics, we find bounds on the bulk viscosity in order to reproduce a matter-dominated era (MDE). This constraint is independent of the interaction term. Some late time phantom solutions are mathematically possible. However, the constraint imposed by a MDE restricts the interaction parameter, in the phantom solutions, to a region consistent with a null value, eliminating the possibility of late time stable solutions with $w<-1$. From the different cases that we study, the only possible scenario, with bulk viscosity and interaction term, belongs to the quintessence region. In this latter case, we find bounds on the interaction parameter compatible with latest observational data.


Introduction
Since the discovered of the present stage of acceleration of the Universe [1,2] many candidates have been proposed to explain such observational result [3,4,5,6,7]. Among them, the cosmological constant, w Λ = −1, remains not only as the simplest alternative but also as consistent with the latest observational data [8]. Despite this, the ΛCDM model is not able to explain those results that still point to a phantom Universe [8], w < −1.
Following the dissipative approach, in [19] was shown that phantom solutions can be obtained by accepting the existence of bulk viscosity whithin the Eckart theory in the ΛCDM model 1 . This result was obtained by using multiple observational tests and considering that the bulk viscosity of some fluid depends on its own energy density, namely ζ j = ζ j (ρ j ). This ansatz avoid the degeneracy problem associate when the bulk viscosity is taken as ζ j = ζ j (H) [19].
Send offprint requests to: 1 Either the bulk viscosity was acting on the radiation, or on the pressureless matter, the crossing of the phantom divide is possible.
However, in [21], the same scenario was studied, from the dynamical system point of view, finding that viscous phantom solutions with stable behavior are not allowed in the framework of complete cosmological dynamics [20,29]. In the present paper we work along these lines by including an interaction between the dark matter and the dark energy. This kind of interaction mechanism has shown to be compatible with the current data [30]. In the context of viscous fluids, the interaction between dark matter and dark energy was studied in [20]. It has been shown that, under the ansatz ζ j = ζ j (H), low-redshift data favors a positive definite value of the bulk viscosity whereas, that high-redshift data prefers negative value of the bulk viscosity. This latter result is in tension with the local second law of thermodynamics (LSLT) [31,32], which it states that for an expanding universe ζ ≥ 0 [33].
In the present work we are interested in extend the results obtained in [19,21] by taking into account an interaction term between dark energy and dark matter and, at the same time, extend the results in [20] by exploring a different functional form for the bulk viscosity 2 .
The paper is organized as follow: in Section 2 we present the field equation of the model. We take into account the contribution of pressureless matter, radiation and dark energy. The first matter fluid is considered as an imperfect fluid, having bulk viscosity in the framework of the Eckart theory [25], whereas the remaining fluid obeying barotropic equation of state (EOS). The bulk viscosity coefficient is taken to be proportional to the dark matter energy density 3 . In Section 3, we study the evolution of the field equations from the perspective of the equivalent autonomous system. We focus our attention in a particular form for the interaction term between the dark matter and dark energy components. A detailed discussion about the viability of a complete cosmological dynamics [20,29] is provided. Important constraints on the bulk viscosity and interaction parameter are obtained. Finally, Section 4 is devoted to conclusions.

The model
We study a cosmological model in a spatially flat FRW background metric, in which the matter components are radiation, dark matter and dark energy. We assumed that the dark matter fluid presents bulk viscosity in the framework of the Eckart theory, whereas the radiation and dark energy are assumed as perfect fluids. Following this set up, the Friedmann constraint, the conservation equations for the matter fluids and the Raychadury equation can be written as: where G is the Newton gravitational constant, H the Hubble parameter, (ρ dm , ρ r , ρ de ) are the energy densities of dark matter, radiation and DE fluid components respectively. Whereas, γ de is the barotropic index of the EOS of DE, which is defined from the relationship p de = (γ de − 1)ρ de , where p de is the pressure of DE. The term Q in (3)(4) is the interaction term between the dark matter and the dark energy components, while 9H 2 ζ in Eq. (3) corresponds to the bulk viscous pressure of the dark matter fluid, with ζ the bulk viscous coefficient. We assume the bulk viscous coefficient ζ to be proportional to the energy density of the dark matter component in the form: where ρ dm0 is the present day value of the dark matter energy density.

The autonomous system
In order to study the dynamical properties of the system (2-4,5), we introduce the following dimensionless phase 3 More specifically ζ ∝ ρ space variables to build an autonomous dynamical system: using the Friedmann constraint (1) is possible to reduce one degree of freedom, namely Ω r = 1 − x − y. Then the equation of motion can be written as where the derivatives are with respect to the e-folding number N ≡ ln a and we have introduced the dimensionaless parameter where, in order to guarantees nonviolation of the LSLT [33,31,32], ξ > 0 4 In addition, in orden to achieve an autonomous system from (8-9) we must define the interaction function Q. If the interaction term is taken as Q = 3Hf (ρ m , ρ de ) [34,35,36,37,38,20,30], then we can introducce a new function hence, the system (8-9) can be written as a two-dimensional autonomous system Imposing the conditions that radiation, dark matter and DE components be positive, definite, and bounded at all times, we can define the phase space of Eqs. (12)(13) as Moreover, we can introduce other cosmological parameter of interest, like the deceleration parameter (q = −(1 + H/H 2 )) and the total effective EOS (w ef f ) in terms of the dimensionless phase space variables (7):  (14), and stability of the critical points of the autonomous system (12)-(13) for γ de = 0 and z = 3βx. The eigenvalues of the linear perturbation matrix associated to each of the following critical points are displayed in Table 2. We have introduced the definitions Table 2. Eigenvalues and some basic physical parameters for the critical points listed in Table 1, see also Eqs. (15) and (16).

Dynamics of the autonomous system
Despite the autonomous system (12)(13) allows to study the dynamics of (2-4,5) for general interaction functions of the form z = z(x, y), we will focus our attention in the particular case of z = 3βHρ de = 3βx [34,35,39,30]. The choice of this particular form is motivated by the requirement of a critical point associated with a MDE in order to explain the structure formation. A simple inspection of (12)(13) shows that the latter requirement implies that (x, y) = (0, 1) thus z(0, 1) = 0. Hence, only those interaction function that fulfill this condition are able to allow a MDE 5 . We also will restrict our analysis to the case γ de = 0. The full set of critical points of (12-13) are summarized in Table 1, whereas the corresponding eigenvalues of the linear perturbation matrix are given in Table 2 3.1.1 Critical points and stability P 1 represents a decelerating solution(q = 1, w ef f = 1/3) dominated by the radiation component, Ω r = 1, and exists unrestrictedly of the sing/value of the interaction and bulk viscosity parameters. However, its stability behavior depends on the value of the interaction parameter β, namely, (i) unstable if β < 4/3 or (ii) saddle if β > 4/3. Critical point P 2 corresponds to a pure dark matterdomination period (Ω m = 1) and always exists. 6 If the condition ξ 2 0 ≪ 0 is satisfied, then this point corresponds 5 See a similar analysis in the case of the ansatz ζi = ζi(H) in [20] 6 Recall that P3, like P1, exists independently of the values/sign of the bulk viscosity and the existence of interaction to the standar matter-domination period, namely w ef f ≈ 0 and q ≈ 1/2. Otherwise w ef f is negative and can behave as an acelerated solution if ξ 0 > 1 3 or even as a phantom solution if ξ 0 > 1. As Table 1 and 2 show, these acelerated solutions are possible in absence of dark energy (x = Ω de = 0). From the stability point of view, P 3 displays two different behaviors, that is: P 3 represents an scaling solutions between dark matter and dark energy components and exists when A background level, P 3 is able to mimic accelerated solutions 7 in the phantom and de Sitter regions, namely: i) phantom region (w ef f < −1) Fig. 2 for more details.
ii) de Sitter region (w ef f = −1) between the dark components. This point has a similar behavior that points 2a in [20] and P2 in [21] 7 Unlike the previous critical points (P1-P2), is not possible to reproduce, in the region of existence, decelerated solutions such as pressureless matter (w ef f = 0) or radiation (w ef f = 1/3).  1. Only if β = 0. We are not considered this case here because that means null interaction between dark matter and dark energy 8 . 8 The case with β = 0 was studied in [21]  Critical point P 4 corresponds to a scaling solution between dark matter and dark energy. This point exists in the region In the existence regions, P 4 is able to mimic only accelerated solutions, namely iii) phantom region (w ef f < −1) Fig. 3 shows, otherwise is stable. 3. If ξ 0 = 2∧β = −1 then w ef f = −2, being an unrealistic value for the effective EOS parameter. iv) de Sitter region (w ef f = −1) 1. As in P 3 , β = 0 leads to a de Sitter solution. As we mention before, this is discarded because it requires a null interaction between dark matter and dark energy 9 v) quintessence region (−1 < w ef f < −1/3). Fig. 4 shows, otherwise is stable.

Cosmological evolution
According to current observational data, any model that aims to make a complete description of the evolution of the Universe must have to follow the complete cosmological paradigm ( [20,21,29]). This paradigm impose transitions between three different evolutions eras from early times to late times, namely: i) radiation-dominated era (RDE), ii) matter-dominated era (MDE) at intermediate stage of evolution, and iii) accelerated expansion. Every one of these statement can be translate into critical point connected by heteroclinic orbits [40,41,42,43].
The condition for a purely RDE (Ω r = 1) is always satisfied by P 1 , independently of the value of the bulk viscosity parameter ξ 0 . Its unstable behavior, given that β < 4/3, guarantees that it can be the source of any solution in the phase space.
For intermediate stages of cosmic evolution, the presence of MDE is needed in order to describe the formation of structures. This matter-dominated period can be recovered by P 2 . This critical point exists independently of the value of the bulk viscosity parameter but a background level, for a no null value of ξ 0 , it behaves as a decelerating solution 10 if 0 < ξ 0 < 1/3 but, if the bulk viscosity takes a sufficiently small value, ξ 0 ≪ 0, is possible to recover 9 Recall that the null interaction case was developed in [21] 10 As the existence of this critical point is also independently of the interaction between dark matter and dark energy (β = 0), this results recover the behavior of P2 in [21] Fig. 5. Vector field in the plane (x, y) for the autonomous system (12)-(13) with γ de = 0. The free parameters have been chosen as (ξ0, β)=(0.0008, 0.038). In this case, the quintessence solution, P4, is the late time attractor of the system, representing an accelerated solution (w ef f = −0.96). The transition from the RDE (P1) to P4 allows for the selection of appropriate initial conditions to recover a true MDE (P3) with w ef f ≃ 0 according with condition (17) [44].
w ef f ≈ 0 and q ≈ 1/2. However this latter statement may be in tension with the recents constraints to the dark matter EOS, which state that −0.000896 < w dm < 0.00238 at the 3σ level [44] using the lastest Planck data realease [8] 11 . Thus, only tiny contribution of bulk viscosity is allowed in order to recover a true MDE with P 2 : this constraints on ξ 0 are also consistent with those obtained in [17,18] in the absence of interaction between dark matter and dark energy. As we mentioned in subsection 3.1.1, P 2 is able to reproduce an accelerated solution given that ξ 0 > 1/3. However, as Table 1 and 1 show, this possible behavior has to be rule out because of the impossibility of finding another critical point Concerning the late time evolution of the Universe, the model has two more critical point capable of providing accelerated solutions, namely P 3 and P 4 . Both represent scaling solutions between dark matter and dark energy. As was discussed in the previous subsection 3.1.1, from the mathematical point of view, is possible to obtain phantom, de-Sitter and quintessence solutions with saddle or stable behaviors depending of the values of the free parameters (ξ 0 , β). If the interaction parameter is negative (β < 0), meaning an energy transfer from dark matter to dark energy, is possible to obtain a late time transition between two phantom solutions: P 3 case i)1 (saddle) → P 4 case iii)1. (stable). This transition requires −ξ 2 0 /4 ≤ β < 0 but, if we also demand a previos stages of RDE and MDE we must impose condition (17), leading to an almost null value for the interaction parameter − 2.00704 * 10 −7 < β < 0, thus the phantom solutions P 3 and P 4 tends to de Sitter solutions w ef f = −1 (β = 0). The rest of the late time phantom solutions demand very large values of the bulk viscosity parameter, ξ 0 > 1, compared to those allowed by (17) in order to recover a true MDE, hence they are ruled out. The only possible late time scenario with a non null value of the interaction parameter corresponds to a stable quintessence solution (P 4 ). This solution requires If we impose the condition (17) to ensure a true MDE and, take into account the latest constraint on the value of the dark energy EOS [8], the following tiny region is obtained for the interaction parameter Fig. 5 shows some example orbits in the plane (x = Ω de , y = Ω m ) to illustrate the above scenario.

Concluding remarks
In this work we studied the dynamics of model of the universe filled with radiation, dark matter and dark energy. The dark matter component was treated as an imperfect fluid having bulk viscosity whereas the remaining fluids were considered as a perfect fluids. The bulk viscosity was taken as proportional to the dark matter density ζ ∝ ρ 1 2 m [19] and, we introduce an interaction term between the dark matter and the dark energy components with the objetive of extend previous results developed in [21]. This new term was taken as Q = 3Hρ de [34,35,39,30].
Recall that the ansatz on the bulk viscosity used in [39,20] (ζ ∝ H) is different from the used in this work. Thus the results obtained now are new compared with those obtained in [39,20] and extend those obtained in [21] by the introduction of the interaction term.
We performed a dynamical system analysis of the model in order to investigate its asymptotic evolution and behavior. The imposition of a transition from a RDE to an accelerated dominated solution, passing through a true MDE reduce the possible values of the bulk viscosity parameter to a tiny region 0 < ξ 0 < 0.000896. This finding is independent of the value of β and support those obtained in [17,18,21] with no interaction between dark matter and dark energy. The presence of an interaction between dark matter and dark energy allows, from the mathematical point of view, to obtain stable(saddle) late time accelerated solutions in the phantom, de Sitter and quintessence regions. However, the requirement of a true MDE imposes strong constraints on the interaction parameter β in the case of late time phantom solutions. In both cases, regardless of the direction of energy transfer between dark matter and dark energy, the interaction parameter is consistent with a null value, hence the de Sitter solution will be the late time attractor. Moreover, the impossibility of having late time accelerated solutions, caused solely by the viscous matter (P 2 ), found in [21] with β = 0, was extended to this new scenario with interaction between dark matter and dark energy.
The only favorable scenario with a no null value of the interaction parameter, 0 < β ≤ 0.039, is described by the late time stable quintessence solution P 4 . This solution is able to fulfill the complete cosmological paradigm, that is a transition between P 1 (RDE) → P 2 (MDE) →P 4 . Recall that this quintessence solution is compatible with the latest constraint on the values of the dark energy EOS [8].