Bulk viscous matter and recent acceleration of the universe based on causal viscous theory

The evolution of the bulk viscous matter dominated universe has been analysed using the full causal theory for the evolution of the viscous pressure in the context of the recent acceleration of the universe. The form of the viscosity is taken as ξ=αρ1/2.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi =\alpha \rho ^{1/2}. $$\end{document} We obtained analytical solutions for the Hubble parameter and scale factor of the universe. The model parameters have been computed using the observational data. The evolution of the prominent cosmological parameters was obtained. The age of the universe for the best estimated model parameters is found to be less than observational value. The viscous matter behaves like a stiff fluid in the early phase and evolves to a negative pressure fluid in the later phase. The equation of state is found to be stabilised with value ω>-1.\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega >- 1.$$\end{document} The local as well as generalised second law of thermodynamics is satisfied. The statefinder diagnostic shows that this model is distinct from the standard Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varLambda $$\end{document}CDM. One of the marked deviations seen in this model to be compared with the corresponding model using the Eckart approach is that in this model the bulk viscosity decreases with the expansion of the universe, while in the Eckart formalism it increases from negative values in the early universe towards positive values.


Introduction
The observational data indicates that the present universe is expanding as well as accelerating [1,2]. Many theoretical models have been proposed to interpret this recent acceleration either by modifying the right hand side of the Einstein's gravity equation with specific forms of the energy momentum tensor T μν , which can cause a negative pressure, or by modifying the left hand side, i.e. the geometry of the space time. In the first approach, one need an exotic cosmic a e-mail: jerinmohandk@cusat.ac.in b e-mail: athirasnair91@cusat.ac.in c e-mail: titus@cusat.ac.in component, dubbed "dark energy", with equation of state satisfying, ω < − 1 3 . The most successful model of the universe, which explains the recent acceleration of the universe, is the standard ΛCDM model. This model incorporates the cosmological constant Λ, characterised by the equation of state ω Λ = −1, as the dark energy. Even though this model fits substantially well with the observational data, it is faced with some drawbacks, mainly the coincidence problem and the cosmological constant problem. The model is unable to explain the observed coincidence between the densities of non-relativistic matter and cosmological constant of the current universe, known as the coincidence problem. The cosmological constant problem is regarding the large discrepancy between the theoretically predicted value of the cosmological constant and its observed value. The value of the cosmological constant predicted from field theoretical estimation is about 10 121 times larger than the observed value. To alleviate these problems, time varying dark energy models have been considered. For the various models of dynamical dark energy, one may refer the review [3] and the references therein. There are two main classes of dark energy models, the quintessence models [4,5], with equation of state ω > −1, and phantom dark energy models with ω < −1. Compared to the quintessence form, phantom dark energy leads to unusual cosmological scenarios, like big-rip [6] where the universe may undergo super-exponential expansion, which effectively rip away the structures in the long run of the expansion of the universe.
There are attempts to explain the recent acceleration without invoking the exotic dark energy component. It was shown by several authors that a bulk viscous dark matter can cause an accelerated expansion of the universe. The effect of bulk viscosity was primarily analysed in the context of acceleration in the early universe, the inflationary epoch [7]. In recent times, the effect of bulk viscous matter in causing the late acceleration was analysed by many [8][9][10][11][12][13][14]. In a recent review, Brevik et al. [15] have shown important implications and capabili-ties of viscosity in describing inflationary era and late time acceleration of the universe; see also Ref. [16].
The detailed mechanism for the origin of bulk viscosity in the universe is still not correctly understood. From the theoretical point of view, the bulk viscosity can originate due to the deviation from the local thermodynamic equilibrium. It manifest as an effective pressure to bring back the system to its thermal equilibrium, which was broken when the cosmological fluid expands (or contracts) too fast. The bulk viscosity pressure thus generated, ceases as soon as the fluid reaches the equilibrium condition.
There exist two main formalism to account for the bulk viscosity in cosmological theories, the non-causal theory where the dissipative perturbations propagate with infinite speed, and the causal theory, where the perturbations are propagating with finite speed. The non-causal formalism was developed by Eckart [17] and many used it in cosmology due to its easiness in analysing the evolutionary behaviour of the cosmological parameters. Later Landau and Lifshitz [18] gives an equivalent formalism. The causal formalism was developed mainly by Israel, Stewart and Hiscock [19][20][21][22][23].
In Eckart theory only the first order deviation from the equilibrium is considered, which effectively leads to the superluminal velocities of the dissipative signals, hence the theory is non-causal [19]. Moreover, the resulting equilibrium states are unstable. However, it illustrates a linear relationship between the bulk viscous pressure and the rate of expansion [24] of the universe. This facilitates the easy analytical method for the parameters in the context of expanding universe.
Based on the Eckart formalism, Brevik and Gorbunova [25] have shown that the viscosity associated with matter, proportional to the expansion rate, can drive the universe into a phantom epoch. Fabris et al. [8] have considered a model with viscous coefficient proportional to ρ ν (where ρ is the density and ν is a constant) and have shown that, for ν = − α + 1 2 , (α is defined by the Chaplygin gas equation for pressre p = −A/ρ α ) the model predicts a late acceleration similar to the generalised Chaplygin gas model of dark energy. They have also concluded that, even though the model is similar to the Chaplygin gas model at the background level, it does not show any oscillations in the power spectrum that plaugues the generalised Chaplygin gas model. This can be considered as a positive indication of bulk viscous models. Later Avelino et al. [26] have studied the bulk viscous matter model using the Eckart formalism, where the bulk viscosity were taken to be proportional to both the velocity and the acceleration of the universe. They have shown that the model can in general predict the late acceleration of the universe. These authors also addressed the asymptotic behaviour of this model and argued that the models are not stable asymptotically. Later in a more general analysis on bulk viscous matter dominated model of the uni-verse based on the Eckart formalism by Athira and Mathew [13], they have proved that the model has considerably good background evolution and is asymptotically stable if the bulk viscous coefficient is a constant. All these analyses were based on the non-causal theory of viscosity. But for physically sound conclusions, one must use the causal theory of bulk viscosity.
As mentioned earlier the causal theory of viscosity was proposed by Israel and Stewart [21,22], taking into account the higher order deviations from the equilibrium, especially the second order deviations, which results in the proper causal connection in the theory. Unlike Eckart theory, the equilibria arising are stable [27]. Moreover, the Eckart theory can be obtained from it as a first order approximation. In certain studies, a truncated version of this theory has been used, where they omit some divergence terms in the expression for the evolution of the bulk viscous pressure, which contains terms corresponding to second order deviations [23,28] from equilibrium.
Initial work where the full causal theory has been used in the context of the inflation occurring during the early period of the evolution of the universe. Maartens and Mendez [29] studied the early inflation caused by bulk viscous cosmic fluid using the full causal theory and found that the resulting solutions are thermodynamically consistent. Another interesting study of the cosmology of flat FLRW bulk viscous universe is in Ref. [30], where the authors find exact solutions corresponds to the early inflationary phase of the universe with a bulk viscous coefficient proportional to the Hubble parameter. In a later work by the same authors using full causal theory, a new class of exact solutions were found by reducing the evolutionary equations of the universe to Abel-type first order differential equations [31]. Zimdahl [32] combined the equivalence between cosmological particle creation and effective viscous fluid pressure using the Israel-Stewart model, finding that there exists an inherent self-limitation to the effective bulk viscous pressure due to the adiabatic particle production. One obtained solutions which indicates a transition from the inflationary phase to later non-inflationary epoch. Zakari and Jou [33] also studied the viscous driven inflationary epoch using the causal theory. Cooley et al. [34] analysed the entropy production in a viscous universe using Israel-Stewart theory. Our concern here is the use of causal viscous formalism in analysing the late accelerating epoch of the universe. We would say that such studies are comparatively less in number in the literature.
After the discovery of the late acceleration of the universe, the full causal theory of viscosity have been used to analyse the late stage of the universe having viscous cosmic components. Cataldo et al. [35] have analysed the possibility of late acceleration, using Israel-Stewart formalism of bulk viscosity. In this work, the authors have used an ansatz for the Hubble parameter (inspired from the non-causal theory), and they have shown that the universe might have undergone a transition to the phantom behaviour leading to big-rip singularity. Piattella et al. [36] have considered a bulk viscous universe using a full causal theory with the aim of unifying dark matter and dark energy. They have found numerical solutions to the gravitational potential using an ansatz for the viscous pressure depends on the density of the viscous fluid and compared it with the standard ΛCDM model. Their conclusion is that, in gross comparison with the standard ΛCDM model, the viscous model with the full causal theory leads to some disfavoured features compared to the truncated version of the model. So by and large in solving the viscous model using the causal theory, it seems that many have used some ansatz either for the Hubble parameter or for the viscous pressure. In the present study, we investigate the evolution of a bulk viscous matter dominated universe using the Israel-Stewart theory of bulk viscosity. We are trying to get the cosmic history by obtaining analytical solutions of the Friedmann equations.
The paper is organised as follows: in Sect. 2 the Hubble parameter for the bulk viscous matter dominated universe is obtained, the behaviour of the scale factor is analysed and the age of the universe is calculated. The evolution of cosmological parameters such as the deceleration parameter, the equation of state parameter, the matter density evolution and the curvature scalar have been discussed in Sect. 3. In Sect. 4, the validity of the local and generalised second law of thermodynamics has been investigated. Sections 5 and 6 deal with the statefinder analysis and estimation of the model parameters, respectively. The conclusions of this study are given in Sect. 7.

FLRW universe dominated with bulk viscous matter
A spatially flat, homogeneous and isotropic universe is described by the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric, where (r, θ, φ) are the co-moving coordinates, t is the cosmic time and a(t) is a scale factor of the universe. The Friedmann equations describing the evolution of the flat universe, dominated by bulk viscous matter, are where H =ȧ a is the Hubble parameter, ρ m is the matter density, P eff is the effective pressure. An overdot represents the derivative with respect to cosmic time t and we have taken c = 8π G = 1. The conservation equation for the viscous fluid iṡ In these equations P eff is given by where p is the normal pressure, given by p = (γ − 1)ρ, γ is the barotropic index and Π is the bulk viscous pressure. The radiation component is avoided and it is a rational simplification as long as we are concerned with the late time evolution of the universe. The bulk viscous pressure in Eckart's theory is of the form Π = −3H ξ , where ξ is the term representing the bulk viscosity of the fluid and as a transport coefficient it can be a function of the Hubble parameter of the universe. For sufficiently large ξ , it is possible that the negative pressure term can dominate and an accelerating cosmology can arise. According to Israel-Stewart causal theory, the effective pressure satisfies the condition where τ , ξ and T are the relaxation time, bulk viscosity and temperature, respectively, and they are functions of the density of the fluid in general, defined by the following equations [37]: where α, β and s are constant parameters satisfying the conditions, α ≥ 0 and β ≥ 0 and r = γ −1 γ . For relaxation time τ = 0, the differential equation for Π reduces to the simple Eckart equation for the viscous pressure. Avoiding the second term on the right hand side of the equation will result in the so-called truncated equation.
The Friedmann equation (2) can be combined with Eqs. (4) and (5) to express the bulk viscous pressure Π as and the time derivative iṡ Then the evolution of the bulk viscosity as given by Eq. (6) can be expressed as where we have used the density dependence of τ, ξ and T as given previously. We are considering non-relativistic matter for which γ = 1 and we also took s = 1 2 [38], implying that the bulk viscosity is directly proportional to the Hubble parameter. The above equation then takes the form where b 1 and b 2 are taken as For calculational purposes we change the variable from cosmic time t to x = ln a; then the above differential equation becomes On solving this we obtained the evolution of the Hubble parameter: where H 0 is the present Hubble parameter, and In Eqs. (17) and (18), is the dimensionless bulk viscous pressure parameter and constants C 1 and C 2 together satisfies, C 1 + C 2 = 1. However, we could obtain the behaviour of Hubble parameter by numerical methods. For this purpose, we evaluated the parameters, α,Π 0 and the present Hubble parameter H 0 in the present model using the cosmological data on supernovae (see Sect. 6). In parameters is shown. The limit of zero viscosity in the nonrelativistic matter implies α → 0, equivalently the viscous pressure Π → 0, and the Hubble parameter will reduces to, H ∼ a −3/2 , which corresponds to the ordinary (non-viscous) matter dominated phase. Equation (16) also shows that the Hubble parameter will become infinite as a → 0. Hence the density will also become infinite at the origin, which suggests the presence of the big bang at the origin.

Behaviour of the scale factor
The Hubble parameter given in Eq. (16) can be integrated, resulting in is the hyper-geometric function with respective arguments. It is to be noted that the hyper-geometric function on the left hand side itself depends on the scale factor. This equation can be used to assess the evolution of the scale factor. The nature of the evolution of a(t) is not explicitly evident from the above equation due to the appearance of hyper-geometric functions in the equation. The behaviour of the scale factor with H 0 (t − t 0 ) for the best estimated values of parameters is shown in Fig. 2. The scale factor is approximately linear at sufficiently early time, which corresponds to the decelerated epoch, and it evolves exponentially with time in the extreme future epoch, corresponding to the de Sitter phase. This asymptotic behaviour indicates the transition from the early decelerated phase to a later accelerated expansion of the universe. The figure also shows that as t → 0, the scale factor a → 0, indicating the presence of the big bang at the origin, hence the age of the universe is properly defined. On comparing Fig. 2 The evolution of the scale factor a with H 0 (t − t 0 ) for the best estimated values of the model parameters with the corresponding model using the Eckart formalism [13] (in which the viscosity is taken to be proportional to both the velocity and the acceleration of the universe) the evolution of the scale factor in the present model is almost similar.
The transition redshift z T , corresponding to the switchover from deceleration to acceleration, can be obtained as follows. From the Hubble parameter in Eq. (16), the derivative ofȧ with respect to a can be written as Equating this to zero, we can get the transition scale factor a T : then the transition redshift z T can be written For the best estimated values of the parameters, the transition redshift is z T ∼ 0.52 +0.010 −0.016 . This is within the WMAP range z T = (0.45 − 0.73) [39]. In Ref. [13], the authors have estimated the transition redshift for a bulk viscous universe using the Eckart formalism to be around 0.49. In the present study using the causal formalism, we have considered only the velocity dependence for the bulk viscous coefficient and the transition is found to have occurred slightly earlier.  2). On equating the scale factor to zero for t = t B , the big-bang time, the age t 0 − t B of the universe can be obtained. The age for different values of H 0 is shown in Fig. 3, where we have used the best estimate of the model parameters. An estimation using the scale factor equation will lead to a simple equation for the age: The age of the universe corresponds to the best estimates of α,Π 0 and H 0 and is found to be around 9.72 Gyr. This is considerably less than the standard value of the age, 13.74 Gyr, deduced from CMB anisotropy data [40], and 12.9±2.9 Gyr, from the oldest globular clusters [41]. Also the age from the present model is less than the age in the corresponding model using the Eckart non-causal formalism, around 10.9 Gyr [13].

The behaviour of deceleration parameter
The deceleration parameter gives a measure of the rate at which the expansion of the universe is taking place. If the deceleration parameter is positive, then the universe is in decelerating phase and vice versa. The deceleration parameter, q, can be expressed as On substituting the Hubble parameter and its derivative, the deceleration parameter becomes Using the best estimates of the model parameter, the coefficients take the values m 1 = 0.31 and m 2 = 5.29. Hence as a → ∞ the term a −m 2 decreases faster than the term a −m 1 , hence the a −m 2 term can be neglected, consequently the deceleration parameter in this limit becomes q → −1 + m 1 .
On the other hand, as a → 0, the term a −m 1 becomes negligibly small, and as a result the deceleration parameter becomes The deceleration parameter of the current epoch, corresponding to z = 0, is For the best estimated values of α andΠ 0 , the present value of the deceleration parameter is found to be q 0 ∼ −0.59 + 0.015 −0.016 , which is quite close to the WMAP value, q 0 = −0.60 [39]. By using Eckart theory and taking the velocity and acceleration dependence for the bulk viscous coefficient, we have q 0 ∼ −0.64 [13]. The evolution of q is shown in Fig. 4. From the figure it is seen that the deceleration parameter will be stabilised around −0.7 in the far future of the evolution of the universe and this is in confirmation with the previously obtained limit q = −1 + m 1 ∼ −0.7. So even though the present model is predicting a never-ending accelerating phase, the universe is not reaching the exact de Sitter phase and this is in marked deviation from the corresponding models using the Eckart formalism [13], in which the model evolves asymptotically to the de Sitter phase. The evolution of the parameter q shows that the model is within the quintessence class. In the earlier epoch, the value of q is positive and considerably large, around q ∼ 4, and it is in confirmation with the previously obtained asymptotic limit q = −1 + m 2 . In the corresponding model with the Eckart formalism, the deceleration parameter was found to be around q ∼ 2 in the remote past of the universe [13]. So the viscous matter behaves as a hard stiff fluid in the earlier epoch and evolves to a negative pressure fluid in the later phase. In this sense, the present model seems to be similar, except for the fact that in the Eckart formalism the model will ultimately evolve to de Sitter phase, while in the causal formalism the fluid is comparatively stiffer in the earlier phase but not ending with an exact de Sitter phase ultimately.

Evolution of equation of state parameter
The equation of state parameter ω has a significant effect on the future expansion profile of the universe. The universe enters the accelerating epoch when ω < − 1 3 . The equation of state can be obtained from the Hubble parameter using the relation where h = H H 0 and x = ln a. Substituting the expression for the Hubble parameter, the equation of state parameter is found to be Since as per the best estimates of the model parameters, m 2 > m 1 , in the limiting condition a → ∞, the equation of state parameter is ω → −1 + 2 3 m 1 ∼ −0.79, which corresponds to its quintessence nature. At a → 0, it becomes ω → −1 + 2 3 m 2 ∼ 2.5 and this corresponds to stiff fluid characteristics. So the viscous matter, being of a stiff fluid nature in the past, evolves to the characteristic of a fluid capable of providing negative pressure as the universe expands. In the Eckart formalism approach, the bulk viscous matter behaves as a stiff fluid with an equation of state equal to + 1 in the past and it becomes − 1 corresponding to the de Sitter phase. The present value of the equation of state parameter is and with the best estimated values of the model parameters, it comes to be around ω 0 ∼ −0.73 +0.01 −0.01 and slightly higher than the value obtained by the combined analysis of WMAP + BAO + H 0 + SN data, around − 0.93 [42,43]. The evolution of the equation of state parameter with redshift for best estimated values of the model parameters is shown in Fig. 5. At this juncture one should note the work by Brevik et al. [44] analysing the possibility of a little rip in which the equation of state approaches − 1 asymptotically from below. In their work they have used a very special equation is a chosen function of the density of bulk viscous matter and ξ(H ) is the bulk viscous pressure, depending on the Hubble parameter. They have shown the possibility of a little rip with some assumed form for f (ρ). On the other hand we have not used any such pre-assumed form for the equation of state. There are also some results telling us that the bulk viscous matter can lead to phantom nature in the early period of the universe [29][30][31]. In the present analysis using the causal viscous formalism due to Israel and Stewart, we get the result that the bulk viscous matter will behave like a strong stiff fluid in the early period and it shows the behaviour of the quintessence dark energy in the later universe, such that the equation of state is stabilised at around ω ∼ −0.79 in the far future of the evolution of the universe.

Evolution of the matter density
The matter density parameter is defined by where ρ crit = 3H 2 0 is the critical density. Using Eqs. (2) and From the above equation the present matter density parameter, Ω m 0 , can be obtained by taking the scale factor a = 1. We get since in the present model only matter is the major component. For zero bulk viscosity and bulk viscous pressure, the parameter takes the value C 1 = 1, C 2 = 0 and m 1 ∼ 3 2 , respectively, Ω m (a) ∼ a −3 , for the usual matter dominated  33), it is seen that as a → ∞ the density will go to Ω m → a −2m 1 , while as a → 0 it behaves as Ω m → a −2m 2 . The evolution of matter density parameter with scale factor for the best estimated values of α andΠ 0 is shown in Fig.  6. The figure shows that, as the scale factor a → 0, the matter density increases rapidly and approaches an infinitely large value. This implies the big bang at the beginning of the universe. The decreasing nature of density in the future depicts the absence of the big rip. In the overall way the evolution of density parameter in the present model is similar to that using the Eckart formalism [13].

Evolution of curvature scalar
The evolution of the curvature scalar of the universe enables one to confirm the occurrence of an initial singularity in the model. The curvature scalar R, for a flat universe, is defined as [45] R = −6(Ḣ + 2H 2 ).
Using the Hubble parameter, the evolution equation of the curvature scalar can be obtained: The above equation shows that the curvature scalar R → ∞ when a → 0, implying an initial singularity corresponding to the big bang. The behaviour of the curvature scalar with scale factor for the best estimated values of model parameters is shown in Fig. 7. This evolution of the curvature scalar suggests the presence of a big bang at the origin of the universe.

Entropy and second law of thermodynamics
In the FLRW universe the bulk viscosity causes the production of local entropy. The law of production of local entropy on the FLRW space-time is expressed as [46], where T is the temperature and ∇ ν S ν is the rate of generation of entropy in a volume unit. The condition for the validity of the second law of thermodynamics then becomes This in turn implies that the bulk viscosity must satisfy ξ ≥ 0. From Eq. (8), the bulk viscosity for s = 1/2 is Substituting Eqs. (2) and (16) in the above equation, we get For the best estimates of the model parameters m 2 > m 1 , the bulk viscosity has the following behaviour. As a → 0, the bulk viscosity satisfies ξ ∼ √ 3α H 0 C 2 a −m 2 . Since the magnitude of m 2 is large, ξ has a large positive value in the early phase. At a → ∞, it evolves to ξ ∼ √ 3α H 0 C 1 a −m 1 . Here the value ξ will be small in the future as the value m 1 is small and will be positive. Corresponding to the present epoch with a 0 = 1, the bulk viscosity becomes ξ = √ 3α H 0 . All these together implies that the viscosity ξ > 0 always. This in turn implies that the condition given in Eq. (38) always is satisfied, hence the second law of thermodynamics is satisfied throughout the evolution of the universe. The evolution of ξ with respect to the redshift z, for the best estimated values of parameters is shown in Fig. 8. Therefore the rate of entropy production is always positive. Here also the present  Fig. 8 The evolution of the bulk viscosity coefficient ξ with redshift z for the best estimated values of parameters model shows remarkable difference from the Eckart formalism model. In analysing the model in the Eckart formalism [13], it was found that the local second law of thermodynamics is violated during an early phase of the universe. Contrary to this, there is no violation of the local second law in the present causal model. In this sense the causal model is to be favoured over the one based on the Eckart formalism. Another distinct behaviour observed in the present model is that the bulk viscosity decreases with time. However, in the Eckart formalism, the bulk viscosity coefficient is increasing from a negative value region to a positive region [13]. The entropy production from the horizon can also be accounted, which leads to the generalised second law (GSL), which states that the sum of the total entropy of the fluid components of the universe and that of the horizon must always increase with time [47,48]. This can be expressed as where S m and S h represent the entropy of matter and that of the horizon, respectively. The apparent horizon radius r A , for a spatially flat FLRW universe, is given as [49] Using Eqs. (2), (4) and (5) the time derivative of r A is obtained: The entropy of the apparent horizon is proportional to the area of the Hubble horizon and is defined as [50] where A = 4πr 2 A is the area of the Hubble horizon. The time evolution of the entropy of the horizon iṡ For the temperature of the apparent horizon we use the relation [51] Using Eqs. (42), (43), (45) and (46), we can write To determine the change in entropy of the matter component, we can apply the Gibbs relation, where T m is the temperature of the bulk viscous matter, E = ρ m V, the total energy of the bulk viscous matter and V = In thermal equilibrium, the temperatures of the viscous matter and that of the horizon are equal, T m = T h . The Gibbs equation (49) can be re-written as The time variation of the total entropy can be obtained by adding Eqs. (47) and (50), Since the radius and the area of the apparent horizon are always positive, from Eq. (51) it is evident thatṠ h +Ṡ m ≥ 0, for a given temperature. Hence the GSL is satisfied. Therefore, in the model in the Eckart formalism [13] and in the present model, the generalised second law of thermodynamics is satisfied.

Statefinder diagnostic
In [ For the present bulk viscous matter dominated model, these parameters take the form here , The evolution of the statefinder parameters in the r -s plane is shown in Fig. 9. The plot exhibits that the trajectory begins from the second quadrant of the r -s plane, r > 0 and s < 0, and enters into the first quadrant, {r, s} > 0. For the present universe, the statefinder parameter takes the form Hence the present values of the statefinder pair is {r 0 , s 0 } = {0.582, 0.128} for the best estimated model parameters. This indicates that the present bulk viscous model is distinguishably different from the ΛCDM model. From Fig. 9 we can also infer that the bulk viscous model resembles the ΛCDM model during an early evolutionary phase of the universe and in the present time the model is moving away from ΛCDM model. In the first quadrant the trajectory is lying in the region r < 1 and s > 0; this represents the quintessence nature. This is a marked deviation from the corresponding model  [13] in the r -s plane.

The model parameter estimation using supernovae data
In this section we describe the evaluation of the model parameters by constraining it with the observational data on type Ia supernovae. We have used the "Union" SNe Ia data set [53], consisting of 307 type Ia supernovae from 13 independent data sets. Our aim here is to extract the best fit for the parameters α,Π 0 and the present value Hubble parameter H 0 . We obtained the parameter values by applying the χ 2 minimisation method. The luminosity distance d L in a flat universe is The difference between apparent and absolute magnitudes of supernovae depends on the distance. The equation that relates the theoretical distance moduli μ th , the apparent magnitude m, the absolute magnitude M and d L is given by The observational distance modulus μ i , obtained from SNe Ia data set is compared with μ th calculated using Eq. (59) corresponding to different values of the redshifts. The χ 2 function can be written as where n is the total number of data points and σ 2 i is the variance of the ith measurement. The best estimated values of the parameters α,Π 0 and H 0 have been obtained by χ 2 minimization. For comparison, the parameter values for the ΛCDM model have also been extracted using the same data. The best estimated parameter values is shown in Table 1 Table  1. It is these values which have been used to generate the evolutionary status of various cosmological parameters present in the previous sections.

Conclusions
In the present work, the evolution of flat FLRW bulk viscous non-relativistic matter dominated universe has been investigated. We have used relativistic second order full causal theory for the evolution of the bulk viscous pressure as given by Eq. (6). The bulk viscous coefficient is taken as ξ = αρ 1/2 , this means it is effectively depending on the expansion velocity of the universe. We solved the Friedmann equations analytically to obtain the Hubble parameter as given by Eq. (16). In the limit of zero viscosity, the present model reduces to the non-viscous matter dominated universe satisfying H ∼ a − 3 2 . We have also obtained the scale factor of the expansion. The asymptotic behaviour of the scale factor indicates the transition from an early decelerated to a late accelerated epoch, and its evolution as shown in Fig. 2 indicates the presence of a big bang at the origin of the universe. Hence the age of the universe is defined and is determined using the best estimated values of the parameters, and it is around 9.72 Gyr. This is considerably less than the age obtained from the CMB anisotropic data and from the data of the oldest globular clusters.
The evolution of the deceleration parameter q is obtained as shown in Fig. 4. From this the transition redshift was obtained, z T ∼ 0.52 +0.010 −0.016 . The present value of the deceleration parameter is obtained: q 0 ∼ −0.59 +0.015 −0.016 , and it is in the range obtained by WMAP data analysis. As a → ∞, q stabilises around the value − 0.7. So, unlike in the Eckart formalism approach, where the bulk viscous universe ultimately goes over to a de Sitter phase, the present model is lying well within the range of quintessence behaviour asymptotically. In its overall evolution, the deceleration parameter begins with q ∼ 4 in the remote past and stabilises around −0.7 in the far future of the universe.
The evolution of the equation of state parameter, ω, is obtained as in Eq. (30) and the variation of it with redshift is shown in Fig. 5. The present value of the equation of state parameter is found to be ω 0 ∼ −0.73 + 0.01 −0.01 . Hence the present universe is accelerating and acceleration has begun in the recent past. The value of ω 0 obtained in this model is slightly higher than that obtained from WMAP + BAO + H 0 + SN data. The future evolution of ω indicates a neverending acceleration phase but not approaching the de Sitter phase. This is a marked deviation from the corresponding model using the Eckart formalism, in which the expansion ultimately ends up with a de Sitter epoch. From the evolution of ω it was found that the equation of state starts with +2.5 in the remote past and evolves to −0.79 in the future stages. This indicates that the bulk viscous matter shows a stiff fluid characteristic in the earlier epoch and then evolves to the quintessence nature in the later stages.
The matter density parameter Ω m , obtained in the present model, is given in Eq. (33) and its evolution with respect to the scale factor is shown in Fig. 6. The evolution of the matter density as a → 0 indicates the presence of the big bang at the origin of the universe. The decrease in density along with the expansion of the universe suggests the absence of a big rip in the future. For zero bulk viscosity the density reduces to Ω ∼ a −3 , corresponding to the ordinary matter dominated era. In the overall way the behaviour of Ω m in the present model is similar to that obtained using the Eckart formalism [13].
The evolution of the curvature scalar R is given in Eq. (36) and is plotted with scale factor in Fig. 7. When a → 0, R → ∞, implying the presence of the big bang at the origin of the universe.
The evolution of the bulk viscosity in the present model shows that it starts with a large positive value in the early phase of the evolution and evolves to smaller values during the later evolutionary phase of the universe. Throughout the evolution it satisfies the condition ξ ≥ 0. Therefore, the entropy production is always positive and, hence, the local second law of thermodynamics is satisfied in the present model. Here also the present model differs from the model using Eckart theory [13]. In a model using the Eckart formalism, the coefficient of the bulk viscosity increases from negative to positive values as the universe expands, hence the local second law is violated in the early epoch. The generalised second law is valid in the present model like in the non-causal Eckart model of the bulk viscous universe.
In order to contrast the present model from the standard dark energy models the statefinder geometric diagnostic has been carried out. The evolution of the model in the r -s plane is shown in Fig. 5. The present value of the statefinder parameter pair is obtained: {r 0 , s 0 } = {0.582, 0.128}. This suggests that the present model is distinct from the standard ΛCDM model. The present bulk viscous model resembles the ΛCDM model in an early phase of the evolution. The evolution of the trajectory on the r -s plane is lying in the region r < 1 and s > 0, this indicates the quintessence nature of bulk viscous matter in the later stages of the expanding universe. In contrast, in the model using the Eckart formalism, the model will approach the ΛCDM model in the future and this is a marked difference from the causal approach.