# Viscous Bianchi Type I Universe with Stiff Matter and Decaying Vacuum Energy Density

Authors: Raj Bali; Pratibha Singh; J. P. Singh
Journal: ISRN Mathematical Physics
 (2012)

---

## Abstract

We examine homogeneous and anisotropic Bianchi type I cosmological model with viscous stiff matter and time-varying cosmological termΛ which scales with Hubble parameter H. The resulting model approaches isotropy. The cosmological term Λ relaxes to a genuine cosmological constant, and the model in the absence of bulk viscosity tends to a deSitter universe asymptotically. Our scenario presents an initial epoch with decelerating expansion followed by late-time acceleration consistent with observations. Bulk viscosity advances the accelerating phase in the model and prevents the matter density to vanish for large times.


---

## Body

## 1. Introduction

Stiff fluid cosmological models create more interest in the study because, for these models, the speed of light equals the speed of sound, and its governing equations have the same characteristics as those of gravitational field [1]. Barrow [2] in his investigation has also pointed out that the entropy level of the universe makes it likely that its initial state was isotropic and quiescent (p=ωρ, ω∈(-1,0)) rather than chaotic only if the equation of state for high density of matter tends to stiff ρ=p (ρ being matter density and p, the isotropic pressure). Keeping in view the importance of stiff fluid models, Bali and Sharma [3], Bali et al. [4], and Mak and Harko [5, 6] have investigated cosmological models for stiff fluid distribution in different contexts.It has been argued for a long time that the dissipative process in early stages of cosmic expansion may well account for the high degree of isotropy we observe today. Dissipative effects involving both the bulk and shear viscosity play a significant role in the early evolution of the universe. The inclusion of dissipative terms in the energy-momentum tensor of cosmic fluid seems to be the best motivated generalization of the matter term of the gravitational field equations. Eckart [7] developed the first relativistic theory of nonequilibrium thermodynamics to the effect of bulk viscosity. Padmanabhan and Chitre [8] pointed out that the presence of bulk viscosity leads to inflationary-like situations in general relativity. Johri and Sudarshan [9]  have shown that bulk viscosity acts like a negative energy field in an expanding universe. Romano and Pavón [10] have investigated the evolution of Bianchi type I universe with viscous fluid. The effect of bulk viscosity on the cosmological evolution has been investigated by a number of authors namely, Saha [11–13], Singh et al. [14], Sahni and Starobinsky [15], Peebles and Ratra [16], Bali and Pradhan [17], Bali and Kumawat [18], and Bali [19].A wide range of observations suggest that universe possesses a non-zero cosmological constant. The cosmological constantΛ is the most favoured candidate of dark energy representing energy density of vacuum. Zel’dovich [20], Dreitlein [21], and Krauss and Turner [22] have studied its significance. Recently Barrow and Shaw [23] suggested that cosmological term corresponds to a very small value of the order 10-122 when applied to Friedmann universe. Linde [24] has investigated that Λ is a function of temperature and is related to the spontaneous symmetry-breaking process. A number of cosmological models in which Λ decays with time have been investigated by Sahni and Starobinsky [15], Bertolami [25], Beesham [26], Berman [27], Singh and Desikan [28], Abdussattar and Vishwakarma [29], Bronnikov et al. [30], Bali and Singh [31], and Ram and Verma [32].In this paper we investigate viscous Bianchi type I cosmological model containing stiff matter and a cosmological termΛ scaling with Hubble parameter H. We find that the model evolves with decelerating expansion in the initial epoch followed by late-time accelerating phase. The model approaches isotropy in the limit of large times, and shear viscosity accelerates the isotropization process. The presence of bulk viscosity prevents the matter density to vanish asymptotically. In the absence of viscosity, the model tends to a deSitter universe for large time [33].


## 2. Basic Equations

The line element for homogeneous, anisotropic Bianchi type I space-time with flat space sections in synchronous coordinate frame is given by(2.1)ds2=-dt2+A2(t)dx2+B2(t)dy2+C2(t)dz2.
We consider the matter component of source field to be viscous fluid described by the energy-momentum tensor
(2.2)TijM=(ρ+p-)vivj+p-gij-2ησij,
where p- is the effective pressure related to the equilibrium pressure p by
(2.3)p-=p-ζvi;i.
Here ρ is the energy density of matter, η≥0 and ζ≥0 are the coefficients of shear and bulk viscosity respectively, vi, the four-velocity vector of the fluid satisfying, vivi=-1, θ is the expansion scalar and σij is the shear tensor defined by
(2.4)θ=vi;i,σij=12(vi;khjk+vj;khik)-13θhij,
where hij=gij+vivj is the projection tensor.The shear viscosity characterizes a change in shape of a fixed volume of the fluid, whereas the bulk viscosity characterizes a change in volume of the fluid of a fixed shape.We choose gravitational units such that8πG=c=1. Since the vacuum has the symmetry of the background, its energy-momentum tensor has the form TijΛ=Λgij, where Λ is a function of time in a homogeneous space. In comoving coordinates system (vi=δ4i), it corresponds to a perfect fluid with energy density ρΛ=Λ and pressure pΛ=-Λ. The Einstein’s field equations with viscous matter and vacuum energy are given by
(2.5)Rij-12Rkkgij=-Tijtotal,
where Tijtotal=TijM+TijΛ=(ρt+pt)vivj+ptgij with the understanding that ρt=ρ+Λ and pt=p-Λ are total energy density and total pressure, respectively. The Bianchi identities require that Tijtotal has a vanishing divergence. The surviving components of the field equations (2.5) for the Bianchi type I metric are
(2.6)p-(ζ-23η)θ-Λ=-B¨B-C¨C-B˙C˙BC+2ηA˙A,(2.7)p-(ζ-23η)θ-Λ=-C¨C-A¨A-C˙A˙CA+2ηB˙B,(2.8)p-(ζ-23η)θ-Λ=-A¨A-B¨B-A˙B˙AB+2ηC˙C,(2.9)ρ+Λ=A˙B˙AB+B˙C˙BC+A˙C˙AC.
Here and henceforth, an overhead dot (·) denotes ordinary derivative with respect to cosmic time t, and semicolon (;) stands for covariant derivative. We define average scale factor R for Bianchi type I space-time as
(2.10)R3=ABC.
Generalized Hubble parameter H and generalized deceleration parameter q are defined as
(2.11)H=R˙R=13(H1+H2+H3),q=-RR¨R˙2=-H˙H2-1,
where H1=A˙/A, H2=B˙/B, and H3=C˙/C are directional Hubble's factors along x, y, and z directions, respectively. For Bianchi type I metric expressions for volume expansion scalar θ and shear tensor σij come out to be
(2.12)θ=3H,σ11=H1-H,σ22=H2-H,σ33=H3-H,σ44=0.
Magnitude σ of the shear tensor σij is given by
(2.13)σ2=12σijσij=16{(H1-H2)2+(H2-H3)2+(H3-H1)2}.
From (2.6) to (2.8), we obtain the following two equations after integration:
(2.14)A˙A-B˙B=k1ABCe-2∫ηdt,B˙B-C˙C=k2ABCe-2∫ηdt.
From (2.10) and (2.14), we obtain that
(2.15)A˙A=R˙R+2k1+k23R3e-2∫ηdt,B˙B=R˙R+k2-k13R3e-2∫ηdt,C˙C=R˙R-k1+2k23R3e-2∫ηdt.
We consider the shear viscosity η scaling with the expansion scalar θ [11], that is,
(2.16)η=3η0R˙R,η0 being constant and bulk viscosity ζ of the form given by Meng and Ren [34]
(2.17)ζ=ζ0+ζ1R˙R,
where ζ0 and ζ1 are constants. Using (2.16) in (2.15), we obtain that
(2.18)A˙A=R˙R+2k1+k23R3+6η0,B˙B=R˙R+k2-k13R3+6η0,C˙C=R˙R-k1+2k23R3+6η0.
Shear scalar σ in this case turns out to be
(2.19)σ=k3R3+6η0;k2=k12+k22+k1k2
implying that the shear viscosity accelerates the isotropization. Field equations (2.6)–(2.9) can be expressed in terms of H, and q as
(2.20)p--Λ=(2q-1)H2-σ2,(2.21)ρ+Λ=3H2-σ2.
From (2.20) and (2.21), we get the following:
(2.22)R¨R=-16(ρ+3p)-2σ23+Λ3+ζθ2,dθdt=Λ-2σ2-12(ρ+3p)-θ23+3ζθ2.
We observe that bulk viscosity and positive Λ drive the acceleration of the universe whereas active gravitational mass density and anisotropy are responsible for deceleration. Also, presence of bulk viscosity slows down the rate of decrease of volume expansion θ.Again from (2.20) and (2.21), we obtain that
(2.23)q=2-ρ-p+2Λ+ζθ2H2
implying that the presence of bulk viscosity lowers the value of deceleration parameter.


## 3. Solutions of the Field Equations and Discussion

We observe that (2.6)–(2.9), (2.16), and (2.17) are six equations in eight unknowns A, B, C, ρ, p, η, ζ, and Λ. We require two more equations to obtain a determinate solution. Following Schützhold [35, 36], we consider the decaying vacuum energy density
(3.1)Λ=aH,
where a is a positive constant and non-vacuum component of matter to be stiff fluid
(3.2)ρ=p.
Let
(3.3)Ω=Λρ
be the ratio between the vacuum and matter densities. From (2.21) and (3.1), we get the following:
(3.4)a=3HΩΩ+1(1-3σ2θ2),
which determines the precise value of a from the observed current values of H, Ω, and σ/θ.Equations (2.20) and (2.21) together with (3.1) and (3.2) give rise to
(3.5)H˙+3H2-aH-32ζH=0
determining the time evolution of Hubble parameter. Equation (3.5) with the use of (2.17) assumes the form
(3.6)H˙+(3-32ζ1)H2-(a+32ζ0)H=0.
We integrate (3.6) to obtain  that
(3.7)H=bd(1-e-bt),
where the integration constants are related to the choice of origin and b=a+(3/2)ζ0, d=3-(3/2)ζ1.On integration, (3.7) yields
(3.8)Rd=m(ebt-1),
where m is constant of integration. For the model, matter density ρ, vacuum density Λ, and cosmological parameters θ, σ, and q have the following expressions:
(3.9)ρ=b(3b-ad)+abde-btd2(1-e-bt)2-k23{m(ebt-1)}(6+12η0)/d,(3.10)Λ=abd(1-e-bt),(3.11)θ=3bd(1-e-bt),(3.12)σ=k3{m(ebt-1)}(3+6η0)/d,(3.13)q=de-bt-1.
We observe that the model has singularity at t=0. It evolves with a big bang at t=0, where ρ, Λ, θ, and σ all diverge. Matter density ρ varies as R-2d whereas vacuum energy density Λ as R-d. Thus, the matter density decays faster than the vacuum energy density. Bulk viscosity component ζ1 arrests the rate of decay. Energy density associated with anisotropy decays as R-6-12η0 contributing to matter production. Therefore, the presence of shear viscosity accelerates the decay of anisotropy energy. At early times Ω≈0 and q≈2-(3/2)ζ1 concluding that matter density is dominant in the beginning of the universe and deceleration of the expanding universe is reduced by the presence of bulk viscous component ζ1. In the limit of large times, we obtain that
(3.14)ρ≈b(3b-ad)d2,θ≈3bd,Λ≈abd,q≈-1,Ω≈ad3b-ad,σθ≈0.We observe that the presence of bulk viscosity does not allow the matter density to become zero in infinitely far future. Also, at late times cosmological termΛ tends to a genuine cosmological constant. Anisotropy σ/θ in the model vanishes asymptotically. The model evolves with decelerating expansion followed by accelerating one at late times. Initial anisotropy of the model dies out asymptotically. We also observe that our solution does not exactly tend to a deSitter universe due to bulk viscosity.As particular cases of the above solution, when only one componentζ0 or ζ1 of the bulk viscosity is present, we can obtain the physical quantities of the resulting models by replacing b by a when ζ0=0, ζ1≠0 and d by 3 when ζ0≠0, ζ1=0 in (3.8) and (3.13). In case only ζ0 is nonzero, we see that matter density varies as R-6 in accordance with the standard model whereas vacuum energy density decays as R-3. In this case q≈2 at the initial epoch. When only ζ1 is nonvanishing component of bulk viscosity ζ, matter density ρ decays as R-2d and vacuum density Λ as R-d. Since d=3-(3/2)ζ1<3, the presence of bulk viscosity component ζ1 slows down the rate of decay. In this case q≈2-(3/2)ζ1 at early times implying that deceleration of the expanding universe goes down by ζ1.We observe that initial epoch of the universe is dominated by matter density, and it decays faster than vacuum energy density. TimetΛ for which the cosmological term Λ turns out to be of the order of matter energy density is given by the condition Λ=ρ. For this condition, we obtain from (3.4) and (3.7)
(3.15)tΛ=1(3/2)ζ0+(3H0/2)(1-σ02/3H02)ln⁡{2d2d-3+σ02/H02-3ζ0/H}.
In the presence of bulk viscosity, the ratio Ω=Λ/ρ tends to a constant in the limit of large times. We also observe that the model starts with decelerating expansion, and expansion in the model changes from decelerating phase to accelerating one. Equating q=0, we can calculate the time tq when expansion changes from decelerating to accelerating phase. We use (3.4) and (3.13) to get  that
(3.16)tq=ln⁡d(3/2)ζ0+(3H0Ω0/(Ω0+1))(1-σ02/3H02).
Therefore, due to bulk viscosity, accelerating phase occurs earlier than in the case of perfect fluid. Time evolution of some cosmological parameters is shown graphically in Figures 1, 2, 3, and 4.Figure 1
Variation of vacuum energyΛ with cosmic time t.Figure 2
Variation of matter energy densityρ with cosmic time t.Figure 3
Variation of deceleration parameterq with cosmic time t.Figure 4
Variation of expansion anisotropyσ/θ with cosmic time t.


## 4. Conclusion

We have examined the possibility of viscous stiff matter distribution in the background of homogeneous anisotropic Bianchi type I space-time with a cosmological termΛ which scales with Hubble parameter H. Coefficient of shear viscosity η is assumed to vary as expansion scalar θ whereas bulk viscosity ζ is taken of the form ζ=ζ0+ζ1(R˙/R), ζ0, ζ1 being constants. Exact solutions of Einstein's field equations have been obtained. The resulting model evolves with decelerating expansion in the initial epoch followed by a late time accelerated expansion presenting an appropriate expansion history of the universe consistent with observations [37–39]. We observe that the model approaches isotropy for large values of t, and presence of shear viscosity accelerates the isotropization process. Cosmological term Λ being very large at initial times relaxes to a genuine cosmological constant asymptotically. In the absence of bulk viscosity, the model results in a deSitter universe for large values of t with H=Λ/3=a/3. We find that the initial phase is dominated by matter density which tends asymptotically to a constant. The ratio Ω between the vacuum and matter densities tends asymptotically to a constant. Presence of bulk viscosity prevents matter density to vanish for large t. We also observe that bulk viscous component ζ1 slows down the rate of decrease of volume expansion. Due to bulk viscosity, accelerating phase occurs earlier than the perfect fluid case.



---

*Source: 704612-2012-06-20.xml*