Cosmology and thermodynamics of FLRW universe with bulk viscous stiff fluid

We consider a cosmological model dominated by a stiff fluid with a constant bulk viscosity. We classify all the possible cases of the universe predicted by the model and analyze the scale factor and the density as well as the curvature scalar. We find that when the dimensionless constant bulk viscous parameter is in the range 0<ζ¯<6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0 < \bar{\zeta }<6$$\end{document} the model begins with a big bang and makes a transition from the decelerating expansion epoch to an accelerating epoch and then tends to the de Sitter phase as t→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ t\rightarrow \infty $$\end{document}. The transition into the accelerating epoch would be in the recent past when 4<ζ¯<6.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$4<\bar{\zeta }<6.$$\end{document} For ζ¯>6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\zeta }>6$$\end{document} the model does not have a big bang and shows an increase in the fluid density and scalar curvature as the universe expands which eventually saturates as the scale factor a→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a \rightarrow \infty $$\end{document} in the future. We have analyzed the model with statefinder diagnostics and find that the model is different from the Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM model but approaches the Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM point as a→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a \rightarrow \infty $$\end{document}. We have also analyzed the status of the generalized second law of thermodynamics with an apparent horizon as the boundary of the universe and found that the law is generally satisfied when 0≤ζ¯<6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0 \le \bar{\zeta }<6$$\end{document}, and for ζ¯>6\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{\zeta }>6$$\end{document} the law is satisfied when the scale factor is larger than a minimum value.


Introduction
Recent studies on the current acceleration of the universe using type Ia supernovae data [1][2][3] have shown that about 72 % of the energy density of the universe is in the form of an exotic component, capable of producing a negative pressure, called dark energy. Cosmological data from a wide range of other sources-the cosmic microwave background radiation [4,5], baryon acoustic oscillations [6], cluster gas fractions [7], and gamma ray bursts [8,9] have all confirmed this conclusion. Concerning the remaining part of the energy density, it was concluded even before the discovery of dark energy a e-mail: titus@cusat.ac.in that 23 % of it consists of weakly interacting matter called dark matter. The evidence for this comes from a variety of observational tests including weak [10] and strong [11] lensing, large scale structure [12], as well as from supernovae and cosmic microwave background [13][14][15] studies. In spite of the fact that all these observational data establish the existence of the components like dark matter and dark energy, the existence of other exotic fluid components has not been ruled out. For example several models predict the existence of an exotic component called dark radiation in the universe [16]. Another exotic fluid which is predicted by several models is a stiff fluid, a fluid with an equation of state p s = ρ s , where p s and ρ s are the normal pressure and density of the stiff fluid, respectively. The equation of state parameter of this fluid assumes the largest value (equal to 1) consistent with causality because the speed of sound in this fluid is equal to the speed of light.
The model with a stiff fluid was first studied by Zeldovich [17]. In recent years a large number of models have been proposed for studying various cosmological properties of stiff fluids. In certain models with self-interacting dark matter components, the self-interaction between the dark matter particles is characterized by the exchange of vector mesons via minimal coupling. In such models the self-interaction energy is shown to behave like a stiff fluid [18]. A stiff fluid has also been considered in certain cosmological models based on Horava-Lifshitz gravity. In Horava-Lifshitz gravity theories a "detailed balancing" condition was imposed as a convenient simplification and the usefulness of this detailed balancing condition was discussed in Refs. [19][20][21]. The stiff fluid appears in such models when this detailed balancing condition is relaxed [22][23][24][25]. Cosmological models with a stiff fluid, based on Horava-Lifshitz gravity, have been studied in Refs. [26,27]. The existence of a stiff fluid have also been found as exact non-singular solutions in certain inhomogeneous cosmological models [28][29][30][31]. The relevance of the stiff fluid equation of state to the matter content of the universe in the early stage of the universe was investigated in Ref. [32]. In Ref. [33], the authors deal with some symmetry related properties of a bulk viscous stiff fluid where they argue that global conformal symmetry is incompatible with a non-vanishing entropy production, but conformal symmetries related to a specific internal interaction that macroscopically corresponds to negative pressure may be compatible with the production of entropy. In such cases one can separate the total pressure into a normal pressure and an effective negative pressure which satisfies the equation of state of the stiff fluid.
The decrease in the density of stiff fluid in the universe is found to be faster than that of radiation and matter; hence its effect on the expansion would be larger in the initial stage of the universe. Primordial nucleosynthesis is an event that took place in the early phase of the universe and a limit on the density of the stiff fluid can be obtained from big bang nucleosynthesis constraints as in Ref. [34].
In a homogeneous and isotropic universe, bulk viscosity is the unique viscous effect capable of modifying the background dynamics. From a theoretical point of view, bulk viscosity arises in a system due to its deviations from local thermodynamic equilibrium [35]. In cosmology, bulk viscosity arises as an effective pressure, restoring the system to its thermal equilibrium, which was broken when universe expanded too fast so that the system did not get enough time to restore the local thermal equilibrium [36][37][38][39]. Several years ago, before the discovery of the present acceleration of the universe, it has been proposed, in the context of the inflationary scenario of the early universe that a bulk viscous fluid can produce acceleration in the expansion of the universe [40][41][42][43][44][45][46]. Recently investigations were made of the possibility of bulk viscous matter causing the recent acceleration of the universe [47][48][49][50].
In the present work we study a stiff fluid dominated cosmological model with bulk viscosity. We consider a universe with effectively one component, the stiff fluid. Such a one component description may not comply with the most general realistic situation. However, this will help us in understanding the essential features of a universe which is dominated by a bulk viscous stiff fluid. We assume a stiff fluid with the equation of state p s = ρ s and the bulk viscosity characterized by a constant viscosity coefficient which is the simplest parametrization for the bulk viscosity. We derive the Hubble parameter, density, equation of state, and deceleration parameter, and we analyze their behavior with regard to further possibilities including the recent acceleration of the universe. The paper is organized as follows. In Sect. 2 we give the basic equations of the FLRW (Friedmann-Lemaitre-Robertson-Walker) universe and derive the general equation for the Hubble parameter in a bulk viscous stiff fluid dominated universe. We classify the different cases depending on the value of the bulk viscous parameter and analyze the evolu-tion of various cosmological parameters. Section 3 contains the statefinder diagnosis of the model. In Sect. 4 we present the status of the generalized second law of thermodynamics. This is followed by the presentation of our conclusions in Sect. 5.

Stiff fluid with bulk viscosity
Stiff fluid cosmological models have generated interest because in these fluids the speed of light is equal to the speed of sound and the governing equations have the same characteristics as that of the gravitational field [51]. The equation of state of the stiff fluid is given as [17] This equation of state corresponds to a special case of the models investigated by Masso and others [52].
In cosmological models the effect of bulk viscosity can be shown to be an added correction to the net pressure p s as where ζ is the constant coefficient of viscosity and H is the Hubble parameter. The form of the above equation was originally proposed by Eckart [53] in the context of relativistic dissipative processes occurring in thermodynamic systems which have gone out of local thermal equilibrium. Later Landau and Lifshitz [54] proposed an equivalent formulation. However, Eckart's theory has got the shortcoming that it describes all the equilibria as unstable [55] and signals can propagate through the fluid with superluminal velocities [56]. Later Israel and Stewart [57,58] proposed a more general theory which avoids these problems and from which Eckart theory appears as the first order limit. However, because of the simple form of Eckart theory, it has been used by several authors to characterize the bulk viscous fluid. For example the Eckart approach has been used in models explaining the recent acceleration of the universe with a bulk viscous fluid [59][60][61][62][63]. Moreover, Hiscock et al. [64] have shown that Eckart theory is to be favored over the Israel-Stewart model in explaining the inflationary acceleration of FLRW universe with bulk viscous fluid. These motivate the use of Eckart's results, especially when one tries to look at the phenomenon of the recent acceleration of the universe. At this point one may also take note of the more general formulation than the Israel-Stewart by Pavon et al. dealing with the thermodynamic equilibrium as in Ref. [65] We consider the flat FLRW universe favored by the recent WMAP observation [66] with the scale factor where a(t) is the scale factor, t is the cosmic time, and (r, θ, φ) are the comoving coordinates. The corresponding dynamical equations are and the conservation equation iṡ where we have adopted the standard units convention, 8π G = 1; an over-dot represents a derivative with respect to cosmic time. From the dynamical equations (4) and (5), we can formulate a first order differential equation for the Hubble parameter by using Eqs. (1), (2), and (6) aṡ The above equation can be expressed in terms of the variable x = log a, suitably integrated and the final result can be written in terms of the scale factor as whereζ = 3ζ /H 0 is the dimensionless bulk viscous coefficient, H 0 is the present value of the Hubble parameter, and we have made the assumption that the present value of the density parameter of the stiff fluid s0 = 1 for a stiff fluid dominated universe.

Classification and evolution of the bulk viscous stiff fluid dominated model
The equation for the Hubble parameter shows that for different values of the viscosity coefficientζ we get different models. In this section we classify different models of the universe which arise due to the different values of the dimensionless viscosity coefficient. We analyze the behavior of the scale factor, the density, and other parameters in these different cases.
This shows that the density of the non-viscous stiff fluid decays more rapidly than the non-relativistic matter or radiation in a FLRW universe, which implies that the effect of the stiff fluid on the expansion of the universe would be larger at early times. So the limit on the density of the stiff fluid can be obtained by considering its effect on the big bang nucleosynthesis. Dutta et al. [34] made an investigation in this regard and found that the change in the primordial abundance of helium-4 is proportional to the ratio ρ s /ρ R , where ρ R is the radiation density. Consequently they found a limit on the non-viscous stiff fluid density as ρ s /ρ R < 30 when the temperature of the universe was around 10 MeV. The evolution of the scale factor can be obtained by integrating the Hubble parameter as The second order derivative of the scale factor with time is This shows that the universe will undergo an eternal deceleration in this case. The behavior of the density following from Eq. (9) reveals that as the scale factor a(t) → 0, the density ρ s → ∞. This implies the existence of a singularity at the beginning of the universe. This fact can be further verified by calculating the curvature scalar for the flat FLRW universe using the equation [ Using the equation for the Hubble parameter and its time derivative it can easily be shown that R ∼ H 2 , which according to the equation H = H 0 a −3 implies that the curvature scalar R → ∞ as a → 0 at the origin, confirming the presence of the initial singularity. So it can be concluded that in this case the universe had a big bang. The time elapsed since the big bang, t B , is found to be Also it is evident from the behavior of the density that, as a(t) → ∞ the density ρ s → 0. In this respect apart from the difference in the dependence on the scale factor, the general behavior of the non-viscous stiff fluid is the same as that of non-relativistic matter or relativistic radiation.

Case
The Hubble parameter is given by Eq. (8). Following the dynamical equations the density of the bulk viscous stiff fluid in this case is given by The evolution of the density with scale factor is given in Fig. 1. This equation shows that as the scale factor a(t) → 0 the density ρ s → ∞, indicating that there is a singularity at the origin. The presence of the singularity is further confirmed by calculating the curvature scale using Eq. (12) and it is which shows that R → ∞ as a(t) → 0, confirming the presence of the initial singularity. So the model of the universe in this case does have a big bang. For finding the scale factor, Eq. (8) can be put in the form which can be integrated to obtain the scale factor as This equation for the scale factor reveals that the time elapsed since the big bang is and hence the age of the universe since the big bang is Taking H 0 = 100 h km/s/Mpc, with h = 0.74, the age of the universe evaluated as per the above equation is around 13.8 Gyr forη = 5.7, a value which is very close to that predicted by the CMB anisotropy data [69]. A plot of the evolution of the scale factor is given in Fig. 2. The scale factor Eq. (17) shows that as t → ∞ the scale factor approaches a form like that of the de Sitter universe, In the early stages of the evolution of the universe, when ζ H 0 [t − t 0 ]/2 < 1, the scale factor can be approximated as These equations of the scale factor at the respective limits show that the universe had an earlier deceleration phase followed by an acceleration phase in the later stage of the evolution. This means that the bulk viscous stiff fluid behaves almost like cold dark matter in the past and behaves like dark energy in the later phase of the accelerating universe, indicating that this model can naturally unify the dark matter and dark energy at least at the background level. The value of the scale factor or the red shift at which the transition from the decelerated to the accelerated expansion occurs depends on the viscosity coefficientζ as shown below. From the Hubble parameter we can calculate the derivative ofȧ with respect a(t) as Equating this to zero, we obtain the transition scale factor as [49] a T = 2(6 −ζ ) ζ and the corresponding transition redshift is From Eqs. (23) and (24) it is clear that forζ = 4 the transition from the decelerated phase to the accelerated phase occurs at z T = 0, a T = 1, which corresponds to the present stage of the universe. In the range 0 <ζ < 4 the transition between the decelerated and the accelerated phase takes place in the future corresponds to z T < 0, a T > 1. The transition takes place in the past of the universe (z T > 0, a T < 1) when 4 <ζ < 6. Whenζ = 0 the value of z T becomes −1 Fig. 3 Evolution of the deceleration parameter q with redshift z and the value of the scale factor a T becomes infinity in the future, this implies that there is no transition to accelerated expansion within a finite time, and the universe is always decelerating. Forζ = 6 the transition takes place at a time corresponding to a T → 0, closer to the big bang.
As a further clarification of the conclusions in the above paragraph we evaluate the deceleration parameter and the equation of state parameter of the bulk viscous stiff fluid in this case. A positive value of the deceleration parameter characterizes a decelerating universe, while a negative value characterizes an accelerating universe. The deceleration parameter q can be evaluated using the equation Using the Hubble parameter from Eq. (8), the deceleration parameter in terms of the redshift z is where we took a = (1 + z) −1 . The evolution of the deceleration parameter is shown in Fig. 3. It is clear from the figure that the deceleration parameter q → −1 in the far future of the evolution of the universe as z → −1 for any positive value of the dimensionless bulk viscous parameterζ . The transition redshift z T can be obtained by equating q to zero, and it leads to Eq. (24). Forζ = 0 the deceleration parameter will be 2, corresponding to a universe dominated by a non-viscous stiff fluid. Forζ = 6 the parameter q = −1 corresponds to the de Sitter phase. So for 0 <ζ < 6 the deceleration parameter is always decreasing from q(a = 0) = 2 to q(a = ∞) = −1, with a transition from positive to negative values corresponding to the transition from deceleration to acceleration in the expansion of the universe. The deceleration parameter today, i.e. for z = 0, is found to be This agrees with our earlier results in Eqs. (23) and (24) that forζ = 4 the universe would enter the accelerating phase from the decelerated expansion at the present time.
Forζ < 4, corresponding to q > 0, we have a decelerating universe today and forζ > 4 corresponding to q < 0, we have an accelerating universe today. From the current observational results [66,69], the present value of the deceleration parameter is around −0.64 ± 0.03, from which the bulk viscous coefficient is seen to beζ > 4 for a universe dominated by a bulk viscous stiff fluid. This analysis shows that if universe is dominated by a bulk viscous stiff fluid, it can take the role of the conventional dark energy, causing the recent acceleration of the universe for a bulk viscous coefficient in the range 4 <ζ < 6. The evolution of the equation of state ω s of the stiff fluid with bulk viscosity can be studied by calculating it using the relation [70] where h = H/H 0 , the weighted Hubble parameter. Evaluating ω s in terms of the redshift z gives The evolution of the equation of state is as shown in Fig. 4. As z → −1, (a → ∞) the equation of state parameter ω s → −1 in the future, which corresponds to the de Sitter universe, irrespective of the value of the viscosity coefficient. Forζ = 0 the equation of state parameter becomes ω s = 1, implying the equation of state for the non-viscous stiff fluid, p s = ρ s . Forζ = 6 the ω s becomes equal to −1. In the range 0 <ζ < 6 the equation of state varies from +1 to −1, making possible a transition from positive to negative values. Even though a negative value of ω s leads to a negative pressure, for the universe to be in the accelerating phase we must have ω s < −1/3.
The present value of ω s is found to be This equation reveals that ω s makes a transition from positive values to negative values at the present time ifζ = 3. While considering the evolution of the q(a = 1) parameter, we have shown that q makes a transition to negative values, giving a universe with an accelerated expansion for ζ = 4. The negativity of the parameter q implies that the universe is accelerating and at the same time the equation of state parameter must be less than −1/3 for the universe to be accelerated [68]. From Eq. (30) it is clear that ω s (a = 1) < −1/3 only forζ ≥ 4. The current observational value of the equation of state parameter of the fluid responsible for the recent acceleration is around −0.94 ± 0.1 [69] and from Eq. (30) we can infer that in a universe dominated by bulk viscous stiff fluid, the corresponding value of the bulk viscous coefficient isζ > 4 causing the recent acceleration. Thus the analysis of the evolution of ω s also shows that the bulk viscous stiff fluid can replace the conventional dark energy in causing the recent acceleration, for 4 <ζ < 6.

Case 3:ζ > 6
Equations (8) and (17) can be used in this case too for assessing the behaviors of the Hubble parameter and scale factor. Forζ > 6 these equations show that the resulting universe will always be accelerating. That is, there is no decelerating epoch at all. When t → ∞ the universe tends to the de Sitter phase. But when t − t 0 → −∞ the scale factor tends to a finite minimum value (see Fig. 2) instead of zero and it is given as The corresponding derivativesȧ andä are zero and hence in this limit the universe becomes an Einstein static universe. As the universe evolves, the scale factor increases monotonically. So there is no big bang in this case and the age of the universe is not properly defined. The curvature scalar can be obtained using Eq. (15). At a = a min , bothä and H are zero, hence the curvature scalar is also zero and it increases as the universe expands, attaining the maximum value R = 5 9 (H 0ζ ) 2 when a → ∞. The density of the bulk viscous stiff fluid follows the same behavior as the curvature scalar (see Fig. 1); the density is zero when a = a min and attains the maximum value (H 0ζ ) 2 /12 as a → ∞.

Statefinder analysis for 4 <ζ < 6
In the analysis in Sect. 2.1.2 we have concluded that there is a transition from decelerated to accelerated expansion in the recent past when 4 <ζ < 6. This gives us hope for the discovery of the recent acceleration of the universe in the context of a universe dominated by a bulk viscous stiff fluid. The behaviors of the scale factor, the q parameter, and the equation of state all show that the bulk viscous stiff fluid is playing the role of dark energy. So we analyze the model using statefinder parameters to compare it with the standard dark energy models. Statefinder parameters [71] are sensitive tools to discriminate various dark energy models, and they are defined as Using the equations for the Hubble parameter (8) and the deceleration parameter (26), the r -s parameter equations can be expressed as The equations show that in the limit a → ∞ the statefinder parameters (r, s) → (1, 0), a value similar to the CDM model of the universe, which shows that the model resembles the CDM model in the far future of the evolution of the universe. A plot of the present model in the r − s plane is shown in Fig. 5, for the bulk viscous coefficientη = 5, and we also found that the plot for other values ofζ show the same behavior. The plot reveals that the (r, s) trajectory is lying in the region corresponding to r > 1 s < 0, a feature similar to that of the generalized Chaplygin gas model of dark energy [72]. On the other hand in comparison with the holographic dark energy model with event horizon as the IR-cut-off [73,74] whose r -s evolution starts in the region s ∼ 2/3, r ∼ 1, and ends on the CDM point, the present model starts in the region r > 1, s < 0, and ends on the CDM point in the r -s plane. Equations (34) and (35) show that forζ = 0, (r, s) = (10, 2) and for higher values ofζ the (r, s) parameter values decrease. The values of the statefinder parameters for the present stage of the universe dominated by a bulk viscous stiff fluid, corresponding to a = 1 (z = 0), is This shows that asζ increases the present values of (r, s) decrease. In Fig. 5 the present position of the universe is marked and it corresponds to (r, s) = (1.25, −0.08), which is different from the CDM model. This means the model presented here is distinguishably different from the CDM model of the universe.

Entropy and generalized second law of thermodynamics
Bulk viscosity may be the only dissipative effect occurring in a homogeneous and isotropic universe. Any covariant description of dissipative fluids, where no matter creation takes place, must satisfy the conservation equation, where T μν is the energy momentum tensor of the fluid in the universe and the semicolon denotes the covariant derivative. The energy momentum tensor in covariant form is given as [75,76], where u μ is the velocity of the observer who measures the pressure p , the form of which is given in Eq. (2). The conservation equation with the above form of the energy momentum tensor will lead to Eq. (6). The bulk viscosity causes the generation of local entropy in the FLRW universe [75,76]. The viscous entropy generation in the early universe was studied in Ref. [77]. It has been established that during the evolution of the universe the sum of the entropies of the fluid within the universe and that of the horizon must always be greater than or equal to zero. This important result is known as the generalized second law (GSL) of thermodynamics. The status of the GSL for the flat FLRW universe with matter and a cosmological vacuum was discussed in Ref. [78]. The status of the GSL in a flat universe with viscous dark energy was discussed in Ref. [79], where the authors have shown that the GSL is valid in the FLRW universe with the apparent horizon as the boundary.
In this section we analyze the validity of GSL in the present model of the universe dominated by a bulk viscous stiff fluid with the apparent horizon taken as the boundary of the universe. The GSL can be formally stated aṡ where S s is the entropy of the stiff fluid and S h is that of the apparent horizon of the universe. The entropy of the stiff fluid within the horizon of the universe is related to its energy density and pressure through the Gibb's relation [80], where V = 4π/3H 3 is the volume of the universe within the apparent horizon with radius r = H −1 and T is the temperature of the fluid within the horizon. We take the temperature T = H/2π equal to the Hawking temperature of the horizon with the assumption that the fluid within the horizon is in equilibrium with the horizon, so there is no effective flow of the fluid toward the horizon. Using the dynamical equation and the net pressure in Eq. (2), the time evolution of the entropy of the bulk viscous stiff fluid within the horizon becomeṡ The entropy of the apparent horizon is given by the Bakenstein-Hawking formula [81][82][83], where A = 4π H 2 is the area of the apparent horizon. Hence the time rate of the horizon entropy iṡ From Eqs. (41) and (43) the GSL condition given by Eq. (39) is satisfied iḟ Using Eq. (7) the above condition becomes As far as H is positive in an expanding the universe, it is evident that the GSL is satisfied in a bulk viscous stiff fluid dominated universe with the apparent horizon as boundary. From Eq. (8), the required condition for the validity of GSL is For 0 <ζ ≤ 6 the above condition is fulfilled and consequently the GSL is well satisfied. But whenζ > 6, the above condition is satisfied only when a ≥ a min as given by Eq. (31) and it is clear from the plot given in Fig. 6 of the Hubble parameter with the scale factor.

Conclusions
In this paper we present a study of the bulk viscous stiff fluid dominated universe model with a constant bulk viscous coefficientζ . A stiff fluid is an exotic fluid with the equation of state parameter ω s = 1, first studied by Zeldovich [17]. We have analyzed the different possible phases of the model according to the value of the dimensionless bulk viscous parameterζ takingζ ≥ 0. Forζ ≥ 0 the model predicts an expanding universe in general. Forζ = 0 the model reduces to a non-viscous stiff fluid dominated universe beginning with a big bang, which is always decelerating with the density varying as ρ ∼ a −6 .
For 0 <ζ < 6 the model corresponds to a universe starting with a big bang and undergoing an initial decelerated expansion followed by a transition to an accelerated phase of expansion at later time. Forζ = 4 the transition from decelerated to accelerated expansion epoch takes place at the present time. For 0 <ζ < 4 the transition to the accelerated expansion phase will take place in the future, but for 4 <ζ < 6 this transition would have occurred in the past. This shows that the bulk viscous stiff fluid can cause the recent acceleration of the universe. From the behavior of the scale factor we have obtained the age of the universe as t 0 − t B = −(2/H 0ζ ) ln(1 −ζ /6).
We have also studied the evolution of the deceleration parameter q and the equation of state parameter ω s for 0 <ζ < 6. Forζ = 4 the deceleration parameter enters the negative region at the present time resulting in a currently accelerated universe. Forζ < 4 the value of q will enter the negative region in the future. For 4 <ζ < 6 it would have entered the negative region in past, implying that the universe made a transition from the decelerated to its accelerated phase in the past. In general q → −1 as a → ∞, corresponding to the de Sitter model of the universe. However, forζ > 6, q is always negative, implying eternal acceleration without any transition from decelerated to accelerated epoch.
The behavior of ω s shows that its value changes from positive to negative when 0 <ζ < 6, which implies a transition from decelerated to accelerated epoch. It is always negative whenζ > 6, which results in an eternally accelerated universe. But irrespective of the value of the viscous coefficient, ω s → −1 as z → −1 (a → ∞). The equation giving the present value of ω s indicates that it would be negative today ifζ > 3, which, however, does not correspond to an accelerating universe. For an accelerating universe ω s < − 1 3 , for whichζ > 4 according to the equation of the present value of ω s .
Statefinder analysis of the model for the range 4 <ζ < 6 predicts the recent acceleration of the universe. The current position of the model in the r -s plane is found to be (r 0 , s 0 ) = (1.25, −0.08) which is different from the CDM model. However, as a → ∞ the statefinder parameters (r, s) → (1, 0), which corresponds to the CDM point.
Whenζ > 6 we have found that as (t 0 − t) → −∞ the scale factor tends to a minimum as given by Eq. (31). In this case the model does not have a big bang. The density and the curvature scalar increase as the universe expands and attain their maxima as a → ∞.
We have analyzed the status the GSL in the present model and found that the GSL of thermodynamics is generally valid with the apparent horizon as the boundary when 0 <ζ < 6. But whenζ > 6 the GSL is satisfied only if the scale factor a > a min where a min is given by Eq. (31).
Finally, the question whether the models presented here can successfully replace the conventional dark energy models and thus could lead to a unification of the dark sectors can be answered only after a full length analysis of the growth of a linear perturbation and structure formation within these models has been carried out. It is clear that the bulk viscous stiff fluid will behave differently from the cold dark matter regarding linear perturbation growth. The bulk viscous pressure can effectively lead to decay of the density perturbation initially and in the later stages it will effectively damp out the density perturbations [47]. Thus, it may be doubtful whether this model can effectively explain structure formation. However, an exact opinion regarding this can only be obtained by checking whether this model can effectively predict the CMBR temperature fluctuations.