Viscosity-Induced Crossing of the Phantom Divide in the Dark Cosmic Fluid

Choosing various natural forms for the equation-of-state parameter w and the bulk viscosity \zeta, we discuss how it is possible for a dark energy fluid to slide from the quintessence region across the divide w=-1 into the phantom region, and thus into a Big Rip future singularity. Different analytic forms for \zeta, as powers of the scalar expansion, are suggested and compared with experiments.


Introduction
The discovery of the accelerated universe [1,2] has led to new concepts and ideas in cosmology, in particular, the concept of a dark energy (for a recent review, see [3]). A characteristic feature of dark energy is the pressure of a negative pressure (i.e., a positive tensile stress) in the cosmic fluid. About 73% of the total mass/energy in the universe consists of dark energy, whereas only 27% consists of a combination of dark matter and baryonic matter. Even modifications of gravity theory itself is a topic that has attracted considerable interest (a review of this kind of theory can be found in [4]).
Usually, one takes the equation of state for the cosmic fluid in the homogeneous form p = wρ, with p the pressure and ρ the mass/energy. Thus w = 0 corresponds to the pressure-less fluid. When w < 0, strange thermodynamic effects are encountered. Thus when the borderline w = −1/3 is crossed, marking the transition into the so-called quintessence region, the strong energy condition ρ + 3p ≥ 0 becomes violated [5]. Even more peculiar properties are encountered when the borderline w = −1 (called the phantom divide) is crossed. A characteristic feature here is that a singularity of the universe may occur, in a finite span of time. It is called a Big Rip. The possibility for such a fate of our universe was first pointed out by Caldwell et al. about 10 years ago [6,7], and has later been studied by a number of authors; cf., for instance, [8]. There are actually several variants of future singularity theories, implying a more 'soft' approach of the universe to this particular limit. Thus scenarios of Little Rip [9], Pseudo Rip [10], and Quasi Rip [11] have recently been discussed in the literature. Our references to the literature are here very limited; a much more extensive overview is given, for instance, in the recent paper [12]. Recent astronomical observations indicate that the value of w lies close to −1, cf. [13], and so a detailed analysis of the behavior of a dark energy fluid is of obvious physical interest. In view of the dominance of the dark energy fluid component in the universe we shall for simplicity consider a model containing one dark component only. Moreover, as an essential point we shall analyze the influence from a viscosity in the cosmic fluid. Most of the earlier cosmological theories have assumed the fluid to be nonviscous. From a hydrodynamicist's point of view this is actually somewhat surprising, since viscosity effects so often play a role in ordinary fluid mechanics. In accordance with common usage we shall take the universe to be spatially isotropic, meaning that it is the bulk viscosity, called ζ, and not the shear viscosity, that becomes relevant. One important property on which we shall focus attention in the following, is that when ζ is taken to be positive it becomes possible for the fluid to slide from the quintessence region (i.e. −1 < w < −1/3) through the phantom divide into the phantom region and thereafter into the future singularity. This was first pointed out in [14]. Obviously, the magnitude of ζ will be important for this transition process. Specifically, we shall consider the following two points: • What is the influence from the equation-of-state parameter w for this process?
• What is the influence from the use of different analytic forms for the bulk viscosity ζ?
These questions will be dealt with in sections 2 and 3 below. Some comparison with astronomical observations will also be made in section 4. Some more references to papers on viscous cosmology are [15] and [16] (these are early papers, the latter being an extensive review up to 1990). Later works can be found in [17,18,19,20,21]. We also mention two works where transition through the phantom barrier was considered in a more general context [22,23].

General formalism, and the case when w is constant
Let g µν be the general metric such that the diagonal components are (−, +, +, +) in the Minkowski case, and let h µν = g µν + U µ U ν be the projection tensor with U µ the fluid's four-velocity. Then, since the shear viscosity is assumed to be zero, the energy-momentum tensor can be written in the simple form where θ = U λ ;λ is the scalar expansion. The FRW metric in comoving coordinates is where a(t) is the scale factor. In this metric θ = 3H, where H =ȧ/a is the Hubble parameter. Defining κ = 8πG we can write the Friedmann equations as The energy conservation equation T 0ν ;ν = 0 implieṡ Let us now consider the equation of state for the dark energy fluid, and first assume that where α is a constant. Thus α = 0 corresponds to the presence of a cosmological constant Λ in conventional relativity, while α > 0 corresponds to the phantom region.
We can now derive the governing equations for the scalar expansion, or equivalently, for the energy density. The governing equation for ρ, taking (8) into account, becomesρ which has the solution Here the integration is taken from present time t = 0 when the density is ρ 0 , to an arbitrary time t in the future.
We consider now different assumptions for the form of the bulk viscosity.
(i) ζ equal to a constant. Let us assume with ζ 0 a constant. From the above equations we get where θ 0 is the present-time expansion and t c the 'viscosity time' The density will vary with time as We can now make the following important observation: If the universe starts from a state lying in the phantom region, α > 0, it will inevitably be developing into a future singularity of the Big Rip type, at a finite time By contrast, if it starts from the quintessence region, α < 0, the universe will never encounter a future singularity. Both θ(t) and ρ(t) tend to finite values as t → ∞. The scale factor a(t) → 0.
(ii) ζ proportional to θ. Let us now make the ansatz This is physically reasonable, as the viscosity may be expected to increase during the violent motions of the cosmic fluid towards the future singularity. Equation (10) yields now From this we see the following: If the universe starts from the phantom region α > 0 at t = 0, it will inevitably end up in a future Big Rip singularity at a finite time, irrespective of the value of the parameter τ 1 . If the the starting point lies in the quintessence region, however, the fate of the universe will depend on how viscous the universe is. The point is whether −|α| + 3κτ 1 is positive or negative. Thus, if τ 1 is larger than a critical value given by the Big Rip singularity (ρ = ∞) occurs. If this condition is not met, ρ(t) → 0 as t → ∞. This property of the universe was pointed out also earlier, in [14].
(iii) ζ proportional to θ 2 . A natural generalization of the above ansatz is to consider the case when ζ(ρ) is proportional to the square of the expansion, We assume that τ 2 , like τ 1 above, are positive quantities, as viscosities should be positive for dissipative processes. (In practice, one would expect that a linear combination of (16) and (19) occurs, but for simplicity we consider here the ansatz (19) alone.) From (10) we now get with If α > 0 (phantom region) the universe thus runs into a Big Rip singularity, ρ = ∞, at a finite time the logarithmic term in (20) fades away. If α < 0 (quintessence region), the situation becomes however complicated. If −|α| + A √ ρ 0 > 0 at the initial instant t = 0, the logarithmic term in (20) fades away when ρ → ∞, but the expression for t becomes negative because of the factor 1/α in the first term in (20). That is unacceptable, since we are looking at the development of the universe in the future only. If −|α| + A √ ρ 0 < 0, a logarithmic singularity (t → −∞) is encountered when √ ρ = |α|/A. We conclude that the case α < 0 is hardly of physical interest here.

A more general form for the equation-ofstate parameter w = w(ρ)
We now make some remarks on the case when w is still taken to be a function of ρ, but has a more general form. Let us start with the ansatz where α and β are unspecified constants to begin with. Here β is nondimensional, while the dimension of α is [α] = cm 4(β−1) in geometric units. The expression (23) means that The previous case (8) corresponds to the choice β = 1. Equation (10) becomes now replaced by Looking for a mathematically simple and at the same time a physical reasonable form for the viscosity, we see that the following form with τ a positive constant, is most natural. Then for β = 1 the case (ii) in the previous section is recovered with τ = τ 1 , and for β = 3/2 the case (iii) is recovered with τ = τ 2 . We see that with (26) the expression between square parentheses in (25) becomes independent of ρ, and so In order to obtain a convergent integral over ρ when the upper limit is chosen as ρ = ∞, one must have β > 1/2. For a Big Rip to occur in a finite time t one must in addition have the condition α + (3κ) β τ > 0 satisfied. We see that the universe possesses the same ability to slide through the phantom divide w = −1 as we saw before: The universe may start from a point in the quintessence region, α < 0, and yet run into a Big Rip singularity if the coefficient τ in (26) is large enough. The condition for Big Rip is seen from (27) to be

Conclusions
Choosing various forms for the equation-of-state parameter w = w(ρ) and the bulk viscosity ζ = ζ(ρ), our main objective has been to discuss the possibilities the dark energy universe has to slide from the quintessence region w > −1 into the phantom region w < −1 and thus into the future Big Rip singularity. The sliding process is thus viscosity-generated. If w is assumed constant, set equal to −1−α in (8), the universe possesses this property in a natural way if ζ(ρ) is taken to be proportional to the scalar expansion θ; cf. (16).
If w(ρ) has the more complicated form (24) the same property persists, if ζ(ρ) is taken to have the form (26), what is a natural generalization. For a Big Rip to occur, the coefficient β in (24) must be larger than 1/2. Finally, it is of physical interest to investigate possible relationships between the assumptions made above, and observations in cosmology. As we know, for an imperfect fluid the viscosity is generated by molecular interactions and can be represented as a functions of macroscopic thermodynamical variables such as temperature. Thus, we can assume the form ζ = ζ(T ). It is natural to make use of conventional kinetic theory. One possibility is to adopt the Chapman-Enskog formula for a dilute fluid (cf. the Chapman-Cowling volume [24]) according to which, for low temperatures (T < 300 K), we can approximate ζ ∝ T 1/2 . Another possibility is to choose the Sutherland formula, implying ζ ∝ T 3/2 . An analysis if this sort was recently given by Wang and Meng [25], comparing with astronomical observations. The temperature was identified with that of cosmic microwave radiation (CMB), with T 0 = 2.73 K the present CMB temperature and z the redshift. One could thus write with ζ 0 an effective constant and α = 1/2 or α = 3/2 in the Chapman-Enskog or Sutherland cases, respectively. Comparing with different observational data sets, Wang and Meng were able to give approximate values for the quantity 12πGζ 0 T α 0 for the two cases mentioned. We reproduce here the values inferred from the SNe Ia data: These numerical estimates are of obvious physical interest. For our purpose here the main issue is however to make a comparison between the exponents: Assume first that the universe is flat and matter dominated, ρ ∝ a −3 , so that a(t) = 2.3 × 10 −12 t 2/3 , T (t) = 10 12 t −2/3 K.
In order to get ζ ∝ θ, our preferred option above, we thus have to set α = 2, not so very far from the Sutherland value 3/2 after all. The case ζ ∝ θ 2 corresponds now to α = 4.