A Viable Dark Fluid Model

TL;DR

This paper proposes a dark fluid cosmological model that interpolates between dust and dark energy states, constrained by observational data and theoretical conditions, providing exact solutions and insights into universe expansion.

Contribution

Introduces a generalized dark fluid model based on an extended equation of state, with analytical solutions and observational constraints to describe cosmic evolution.

Findings

Model can interpolate between dust and dark energy

Parameters constrained by supernova data and causality conditions

Exact solutions for luminosity distance and energy density evolution

Abstract

We Consider a cosmological model based on a generalization of the equation of state proposed by Nojiri and Odintsov [46] and \v{S}tefan\v{c}i'c [47], [48]. We argue that this model works as a dark fluid model which can interpolate between dust equation of state and the dark energy equation of state. We show how the asymptotic behavior of the equation of state constrained the parameters of the model. The causality condition for the model is also studied to constrain the parameters and the fixed points are tested to determine different solution classes. Observations of Hubble diagram of SNe Ia supernovae are used to further constrain the model. We present an exact solution of the model and calculate the luminosity distance and the energy density evolution. We also calculate the deceleration parameter to test the state of the universe expansion.

Tables

Table 1. Table 1: Asymptotic behavior of the equation of state

Case	ranges of $r$ and $s$	large $H$	small $H$
$I$	$r\succ 0$ $s\succ 0$	$\dot{H}\to \frac{-\alpha }{2\beta }{H}^{r-s}$	$\dot{H}\to \frac{-\alpha }{2}{H}^r$
$II$	$r\succ 0$ $s\prec 0$	$\dot{H}\to \frac{-\alpha }{2}{H}^r$	$\dot{H}\to \frac{-\alpha }{2\beta }{H}^{r-s}$
$III$	$r\prec 0$ $s\succ 0$	$\dot{H}\to \frac{-\alpha }{2\beta }{H}^{r-s}$	$\dot{H}\to \frac{-\alpha }{2}{H}^r$
$IV$	$r\prec 0$ $s\prec 0$	$\dot{H}\to \frac{-\alpha }{2}{H}^r$	$\dot{H}\to \frac{-\alpha }{2\beta }{H}^{r-s}$

Equations

$$p=-A/{\rho }^{\alpha }$$

$$p=A\rho -\frac{B}{{\rho }^{\alpha }}$$

$$d{s}^2=d{t}^2-a^2(t)(d{r}^2+r^{2}d{\Omega }^2)$$

$$T_{\mu \nu }=(\rho +P)U_{\mu }{U}_{\nu }-P{g}_{\mu \nu }$$

$$\frac{{\dot{a}}^2}{a^2}=\frac{8\pi G}{3}\rho$$

$$\frac{\ddot{a}}{a}=-\frac{4\pi G}{3}(\rho +3P)$$

$$\dot{\rho }=-3\frac{\dot{a}}{a}(\rho +P)$$

$$p=-\rho -f(\rho )$$

$$P=-\rho +\frac{\gamma {\rho }^n}{1+\delta {\rho }^m}$$

$$\omega =-1+\frac{\gamma {\rho }^{n-1}}{1+\delta {\rho }^m}$$

$$H^2=\frac{1}{3}\rho$$

$$\frac{\ddot{a}}{a}=H^2+\dot{H}=-\frac{1}{6}(\rho +3P)$$

$$\dot{H}=-\frac{1}{2}(\rho +P)$$

$$\dot{H}=-\frac{\alpha H^r}{2(1+\beta H^s)}$$

$$r=2n;s=2m;\alpha =3^n\gamma ;and\beta =3^m\delta$$

$$\frac{\alpha }{2}(t_0-t)=\frac{1}{1-r}\left({H_0}^{1-r}-H^{1-r}\right)+\frac{\beta }{s-r+1}\left({H_0}^{s-r+1}-H^{s-r+1}\right)$$

$$\ddot{a}=\dot{H}+H^2=0$$

$$\dot{H}=-H^2$$

$${c_s}^2=\frac{dp}{d\rho }$$

$$\frac{dp}{d\rho }\ge 0$$

$$p=-\frac{A}{{\rho }^{\alpha }}$$

$${c_s}^2=\frac{\alpha A}{{\rho }^{\alpha +1}}$$

$${c_s}^2=-1+\frac{\left[\left(n-m\right)\delta {\rho }^m+n\right]\gamma {\rho }^{n-1}}{{\left(1+\delta {\rho }^m\right)}^2}\le 1$$

$$\left(n-m\right)\frac{\gamma }{2\delta }{\rho }^{n-m-1}\le 1$$

$$n\le m+1$$

$$\frac{\left(n-m\right)}{2}\gamma \delta {\rho }^{n+m-1}+\frac{n}{2}\gamma {\rho }^{n-1}\le 1$$

$$n+m\ge 1andn\ge 1$$

$$1-m\le n\le 1+m$$

$$2(1-s)\le r-s\le 2$$

$$\dot{\rho }=-3H\frac{\alpha H^{s+2}}{1+\beta H^s}$$

$$\frac{dH}{da}=-\frac{1}{2a}\frac{\alpha H^{s+1}}{1+\beta H^s}$$

$$H(a)=a^{-\alpha /2\beta }\exp \left[\frac{1}{s}W\left(\frac{1}{\beta }{a}^{\frac{\alpha s}{2\beta }}{e}^{C\frac{\alpha s}{2\beta }}\right)e^{\frac{-\alpha }{2\beta }C}\right]$$

Full text

1 Introduction

Recent cosmological observations favoring the scenario of spatially flat and accelerated universe. Observations come from Hubble diagram of SNe Ia supernovae are best fitted by spatially flat accelerated cosmological models [1] - [4]. Besides, other observations from Cosmic Microwave Background Radiation (CMBR) and Baryon Acoustic Oscillation (BAO) are also strong evidences for this scenario [5] - [11]. This currently observed speeding up expansion of the universe leads to the fact that there is a repulsive force acting in the universe space-time.

According to our current understanding of the universe history, the radiation era was followed by a matter dominated era during which most of the universe structures like stars, galaxies and galaxy clusters we now observe were formed by gravitational instabilities. Before these ordinary matter structures were formed, a particular form of gravitating but non-baryonic sort of matter structures was formed. This sort of matter which is non-relativistic does not interact with radiation, and so is called Cold Dark Matter (CDM). The perturbations of this dark matter (DM) collapse first because its particles don’t couple to photons as they don’t interact with them. Their gravitational potentials then guide for baryon collapse [12] - [14]. Late in this matter dominated phase, a new but very peculiar form of invisible and non-gravitating "matter" started to dominate. This sort of "matter" does not interact with baryons, radiations, or any other sort of visible matter. It is now known as Dark Energy (DE).

Observational evidences show that, see for ex. [1] and [15] - [17], our universe at present contains approximately $4\%$ only of radiation and baryonic matter that is well known by the standard model of particle physics and can be detected in laboratory and so is considered visible. While about $26\%$ is the non-baryonic dark matter. The rest of our universe content, which is about $70\%$ is the exotic component known as dark energy. As a result, this dominated dark energy component antigravitates leading to the present observed universe expansion.

Although dark energy is an acceptable and effective description for such large scale homogeneous and isotropic accelerated universe, identifying its origin and nature is one of the most debated questions in physics and cosmology. Large variety of models are proposed to describe the nature of such dark energy [18], [19]. As the energy density due to DE remains unchanged throughout the universe space-time, the simplest choice was to revive the idea of Einstein about the cosmological constant $\Lambda$ [20], [21]. This together with Cold Dark Matter (CDM) provides the standard model for cosmology $\Lambda$CDM [21] - [24], which appears to be in a very good agreement with observational data. In Einstein’s field equation, this term has the concept of intrinsic energy density of vacuum [25]. However, this arises an important question, why the observed value of vacuum energy density is very far below that is predicted from particle physics?. This is the famous fine tuning problem. Another important question is why the energy densities of the vacuum and the matter are of the same order of magnitude today?, the famous coincidence problem.

Dynamical models are other alternative to the cosmological constant. Some of these are the minimal coupled scalar fields or quintessence [26] - [29], k-essence [30] - [33], Chaplygin gas [34] - [36], all have an equation of state parameter $\omega\geq-1$. These models adopt the idea of the possibility of unification of the two dark sectors, dark energy and dark matter. Models with $\omega\prec-1$ knowing as phantom energy models [37], [39] consider a growing vacuum energy with expansion. A more general class of models considers an equation of state parameter varies with the scale factor, or the red shift [20], [23], [24], [40], [41]. There are also some observational analysis which consider the variation of $\omega$ with the red shift and with time [42], [43]. Modification of the governing equation, or the so called modified gravity theory [44], [45], is also one way of thinking (For a general review of different DE models see [46]).

The idea of unification of DE and DM is a promising one and Chaplygin gas model is one of the most acclaimed models representing it. This model was proposed by Kamenshchik et. al. as early as 2001 [47] where they assumed that the universe, within the framework of standard cosmology, is filled with an exotic fluid obeys the Chaplygin gas equation of state (EoS)

[TABLE]

where $A$ is a positive constant and $\alpha=1$. This model has the advantage of smooth transition between different phases of the universe. Many authors then extended this model using the idea of generalizing or modifying the EoS of the background fluid in order to improve the behavior of the fluid throughout the universe evolution. In 2002, Bento et. al. [48] extended the model to a more general form known as the generalized Chaplygin gas (GCG) model where the $\alpha$ constant is considered in the range $0\prec\alpha\leq 1$. Other various models based on Chaplygin gas were also proposed such as the modified Chaplygin gas model introduced in 2002 by Benaoum [49] where he could interpolate between standard fluids at high energy densities and Chaplygin gas fluids at low energy densities using the EoS

[TABLE]

where $\alpha\geq 1$ and $A$ and $B$ are positive constants, and the hybrid Chaplygin gas model introduced by Bili$\acute{\texttt{c}}$ et. al. in 2005 [50] and leads to transient acceleration.

In this work, we have adopted a class of dynamical models in which the EoS has the advantage of interpolation between two different Chaplygin Gas models. Our fluid is a barotropic one with an EoS that can be considered as a correction to the vacuum EoS by a one asymptotes between two power laws, so that it has the advantage of interpolation between the two equations of state of the DE and DM, which might describe some sort of smooth phase transition. Another advantage of this model is that it is not only interpolating between dust and DE in early and late times, but also it has a more general EoS for DE that enables the cosmological constant as a special case. This has the more advantage of a general asymptotic solution than Chaplygin gas which goes to the cosmological constant. We present an exact solution of the equations which enables more accurate results for the calculations of the luminosity distance, while observations of Hubble diagram of SNe Ia supernovae are used to further constrain the model parameters. We study the energy density and the equation of state parameter evolutions. The analysis of the EoS shows that the dynamical properties of both dark sectors can be described through this model which lets our choice seems natural because of the smooth transition. We also calculate the deceleration parameter to examine the expansion of the universe due to the model.

The rest of the article is organized as follows; in the following section we have formulated the model. In section 3 the mathematical treatment of the model is made. In section 4 we study the phase diagram of the model and its classes of solutions. Section 5 treats the first case of study of our model where the causality constraints are examined and the dynamics of the universe is studied where the equation of the luminosity distance as a function of red shift is solved and parameters are constrained through Hubble diagram observations, the energy density and the deceleration parameter evolutions due to the model are examined. Section 6 treats the second case of study of our model. Section 7 concludes our work.

2 Building the Model

We are working through the standard FRW cosmology where the metric in the spatially flat geometry is

[TABLE]

where we consider units with $c=1$. Using the energy momentum tensor given by

[TABLE]

where in comoving coordinates $U_{\mu}={\delta^{0}}_{\mu}$, the Einstein’s equation will lead to the Friedman equations

[TABLE]

while the conservation equation ${{T_{\nu}}^{\mu}};\mu=0$ gives

[TABLE]

Now we have to consider an EoS to close the system. It is a fact that dark energy dominates the universe today and an adequate EoS for this component is $p=\omega\rho$. Observations indicate that the EoS parameter $\omega$ is close to $-1$ and people are now targeting to confirm whether $\omega$ is constant or varies as the universe expands. Owing to this fact, Nojiri and Odintsov [51] considered a general EoS of the form

[TABLE]

where $f(\rho)$ is an arbitrary function may be considered as a correction to the standard DE equation of state. $\check{\texttt{S}}$tefan$\check{\texttt{c}}$i$\acute{\texttt{c}}$ [52], [53] studied in details the case of $f(\rho)\propto\rho^{\alpha}$ with a constant $\alpha\neq 1$, and Nojiri and Odintsov [54], [55] studied the future singularity associated with this model.

In this work we adopted a barotropic fluid with an EoS which is also a correction to the standard dark energy EoS. The fluid pressure has a general form for the density dependence

[TABLE]

Where $\gamma$, $\delta$, $n$, and $m$ are constants which are considered as free parameters. This enables interpolation between different powers for the density, so that the phase transitions during the universe’s evolution are described smoothly. The EoS parameter is given by

[TABLE]

3 The model

Let’s now study the evolution of the physical quantities of the universe due to our model. Considering units with $8\pi G=1$, eqns (5)-(7) reduce to

[TABLE]

where $H=\frac{\dot{a}}{a}$ is the Hubble parameter. Combining eq(9) to eq(13) we get

[TABLE]

where again $\alpha$, $\beta$, $r$, and $s$ are constants related to those of (9) through the relations

[TABLE]

The solution of equation(14) indicates the evolution of the Hubble parameter with time. This equation is integrated to give

[TABLE]

However, our model is not describing a pure dark energy, in which case the values of the parameters are totally free and are not constrained. Instead, our model able to follow the smooth phase transition which takes place through the universe evolution. This restricts the parameter values to obey the limits satisfied by observations so that for early time we must have the equation of state of the perfect fluid.

Accordingly, we first have to study the asymptotic behavior of $\dot{H}$. From one side, this indicates to what extent our model agrees the real evolution of the universe, and from the other side may also fixes the limits of some of our model parameters. Table 1 shows this behavior for the different ranges of of the parameters $r$ and $s$.

One can see from the table that the values of negative $r$ must be precluded in our calculations as a result of the asymptotic behavior of the $\dot{H}$ function. This is because for asymptotically perfect fluid behavior of our universe, the $\dot{H}$ function at large $H$ forces $r$ to be positive if $s$ is negative, and forces $(r-s)$ to be positive if $s$ is positive, which does not satisfy constraints for the pre-assumptions.

On the other hand, the asymptotic behavior of $\dot{H}$ for the two cases of positive $r$ can attain the perfect fluid behavior at large $H$ if $r-s=2$ and $\alpha/\beta=3$ for the case of $+ve\;s$ and if $r=2$ and $\alpha=3$ for the case of $-ve\;s$. These two cases will be studied in details.

4 Analyzing the Model

The phase space analysis method enables information about the behavior of the cosmological model and its classes of solutions. In Fig.1 we show a phase diagram, $H-\dot{H}$ diagram, for our model agaist $\Lambda CDM$ model, where $h=\frac{H}{H_{0}}$ and $\tau=\frac{t}{t_{0}}$. In the Fig. we also plot the zero acceleration curve for which

[TABLE]

so that

[TABLE]

Due to [56], since the function $F(h)=dh/d\tau$ is continuous and differentiable and since we have future fixed point, our model has a unique solution for times $\tau\succ 0$ and has no future finite time singularities. instead, it evolves to a de Sitter space.

We can see from the Fig. that deceleration parameter changes sign from $+ve$ to $-ve$ as the universe evolved from matter dominated era to dark energy dominated era. We can also see that our model expects matter-DE equality time to be earlier than expected by $\Lambda CDM$ model.

5 Case I: $s$ is positive

This is the first case that match the dark fluid dynamics. The asymptotic behavior of the EoS at large $H$, early times, has the form $\dot{H}{\rightarrow}\frac{-\alpha}{2\beta}H^{r-s}$, so that for $r-s=2$ and $\alpha/\beta=3$ we have the perfect fluid EoS at early times. Let’s now examine the constraints that can be obtained by applying the causality condition to this model.

5.1 Causality Constraints

Spherical sound waves are lunched as a result of the overpressure due to the initial overdensity in the DM and gas. These oscillations have a speed which is an indication to the velocity by which such perturbations are transmitted. Accordingly, it must not exceed the speed of light, which is a condition of causality.

For a barotropic cosmic fluid, the speed of these sound waves is defined by the relation

[TABLE]

As a result, it is an imprint for the barotropic models to have a value of $dp/d\rho$ constrained by limits on the sound speed [57]. This is because $dp/d\rho$ must satisfy the condition

[TABLE]

to insure stability in perturbation growth, so that we must have ${c_{s}}^{2}\geq 0$. On the other hand, causality constrains ${c_{s}}^{2}$ not exceed the speed of light [58], [59], so that ${c_{s}}^{2}\leq 1$. This fact distinguishes the dynamical behavior of barotropic fluids from quintessence models which describe an imperfect fluid and generically have $\frac{dp}{d\rho}\prec 0$ [60]. As a prototype example for barotropic fluids is the Chaplygin gas [47] and the generalized Chaplygin gas models [48] which have an EoS

[TABLE]

where $A$ and $\alpha$ are $+ve$ constants. For these models we have

[TABLE]

which ensures ${c_{s}}^{2}\geq 0$ as $A$ and $\alpha$ are $+ve$.

Accordingly, applying constraints for barotropic fluids, one may get an indication for the available ranges for the values of the free parameters in the EoS. Now using (9), we get

[TABLE]

Now for our case where $n$ and $m$ are $+ve$, relation (23) for large $\rho$ gives the condition

[TABLE]

This is fully satisfied for

[TABLE]

While for small $\rho$, the condition will be

[TABLE]

This is satisfied if

[TABLE]

Of course, this last condition is automatically satisfied if the the condition in (25) is satisfied. Now Combining our conditions, (25) and (27), we get to the conclusion that for this case we must consider

[TABLE]

Which means

[TABLE]

In agreement with the conditions required by asymptotic behavior of the EoS at early time for this case. In fact, for the general values of $\rho$ we also have a large dependence on $\gamma$ and $\delta$. Accordingly, we have to make a check using the obtained optimized values for the parameters.

5.2 The dynamics of the universe

In this section we study the ability of our model to produce the physical quantities that can be really measured. This, of course, depends on the compatibility of the results of our model with the astrophysical data that are indeed observed.

5.2.1 The Luminosity Distance

One of the most important observed data is that of the Hubble diagram of type SNe Ia supernovae, which is a plot of the luminosity distance $D_{L}$, that determines the flux of the source, as a function of the red shift. We will use these data to constrain the parameters in our model. Due to the asymptotic and causality constraints, we have now only two out of four parameters that are free and needed to be optimized due to observations.

To satisfy the causality and asymptotic behavior constraints, we choose $(r-s)$ as $2$ while $(\alpha/\beta)$ has to be $3$, so that we can attain the perfect fluid behavior of $\dot{H}$ at large $H$. Let’s now start with our model by using the relation (14) for $(r-s)=2$ in (13), so that we have

[TABLE]

which in turn gives

[TABLE]

This has the solution

[TABLE]

The function $W(x)$ is the Lambert W-function. To fix the constant C we consider that at some time where $t=t_{0}$ we have $a=a_{0}$ and $H=H_{0}$, this gives

[TABLE]

So that

[TABLE]

Accordingly

[TABLE]

the luminosity distance $D_{L}$ is then given by

[TABLE]

where $D_{p}(z)$ is the proper distance given by

[TABLE]

Accordingly

[TABLE]

The integration of the above equation can not be solved analytically. Accordingly, we used numerical methods to solve the integral. This enabled exact solution for the relation controlling the dependence of the luminosity distance on the redshift. The comparison of our model against supernovae data relies on the $\chi^{2}-$statistics where

[TABLE]

where $\mu$ is the distance modulus. In Fig.2 we plot our model against data set from Conley et al. [61]. Data set are available at Supernova Cosmology Project. Results from $\Lambda CDM$ model [63] are also shown in the Fig. for comparison.

In terms of distance modulus, data points are

[TABLE]

where $M$ is the absolute magnitude which is assumed to be constant for a standard candle like SNe Ia, while $m$ is the apparent magnitude given by

[TABLE]

Using $H_{0}=67.8Km/s/MPc$ due to Plank $+$ WMAP $+$ BAO measurements [64], we obtained a best fit of our model at $\alpha=0.004$ and $s=1.369$ for which ${\chi^{2}}_{red}=1.020$. This result has a $Q-$value $=0.34$ which shows that data are consistent with the model. Fig.3 shows the confidence regions in parameters space.

Having the full set of parameters, let’s test the causality condition using these parameters. A plot of causality relation (23) using our parameters for this case is shown in Fig(4). We can see from the figure that the parameters satisfy the conditions for ${c_{s}}^{2}$ throughout the whole stages of the evolution of the universe.

5.2.2 The Energy Density Evolution

Let us see how the energy will evolve with time due to the model in a homogenous and isotropic universe. Relations (7) and (9) give

[TABLE]

To simplify our way of solution to this equation let’s use our parameters. In this case we have $n=m+1$, so that the solution of relation (42) will take the form

[TABLE]

Where

[TABLE]

To study the properties of relation (43) at early times we make use of the property

[TABLE]

then

[TABLE]

Now, in the past, where $a\prec\prec a_{0}$, so that $x\prec\prec 1$ due to (44) as $m$ is $+ve$, Taylor expansion of $W$-function can be approximated to $W(x)\approx x$, so that the exponential term of the $W$-function converges to $1$, then the energy density for that time will take the asymptotic form

[TABLE]

as we have $\gamma/\delta=1$ for this case, where $A=\rho_{0}e^{\frac{-1}{\delta m{\rho_{0}}^{m}}}$ is a constant. Accordingly, the fluid of the model behaves as a matter in the earlier times.

On the other hand, for $a\succ\succ a_{0}$, $x\succ\succ 1$, the energy density function due to relation (46) will tend to

[TABLE]

Accordingly

[TABLE]

Using the value of $m$ we can see that this function is a very slowly decreasing function of $x$ for $x\succ 1$, and tends to be nearly constant at $x\succ\succ 1$, i.e. at $a\succ\succ a_{0}$, a signature for a cosmological constant towards the future.

The equation of state parameter $\omega$ for this case, where $m=n-1$, will be given, from (10), by

[TABLE]

Using our parameters for this case where $\delta/\gamma=1$, we finally get

[TABLE]

Indicating that for earlier times, where we have large values of $\rho$, $\omega\rightarrow 0$, a signature for matter era, while changes to $-1$ gradually as the matter content of the universe diluted.

5.2.3 Deceleration Parameter and the Universe Acceleration

The acceleration of the universe is related to a dimensionless parameter known as the deceleration parameter $q$ defined as

[TABLE]

Since you have to have $\frac{\ddot{a}}{a}\succ 0$ to get an accelerating universe, the deceleration parameter will then measure if the expansion of the universe is accelerating or decelerating. If $q\prec 0$, the expansion is accelerating, while if $q\succ 0$ it is decelerating. Now as

[TABLE]

Then

[TABLE]

Now using (14), one gets the deceleration parameter for our model to be

[TABLE]

Using $r=s+2$, relation (55) gives

[TABLE]

Using the values of our parameters with $H_{0}=67.8$, the present day deceleration parameter due to this model case has the value

[TABLE]

Relation (56) also shows that deceleration parameter evolves towards more negative values as the Hubble parameter decreases with time.

6 Case II: $s$ is negative

For this case $\dot{H}\rightarrow\frac{-\alpha}{2}H^{r}$ asymptotically at large $H$, which means that we can attain the perfect fluid EoS at early time if $r=2$ and $\alpha=3$. However, a first step again is the examination of the causality constraints for this case of the model.

6.1 Causality Constraints

In this case $m$ is $-ve$, so that relation (23) at large values of $\rho$ tends to

[TABLE]

Which is satisfied if

[TABLE]

while for small $\rho$ we have

[TABLE]

As $m$ is $-ve$, this is always satisfied under the condition

[TABLE]

Combining (59) and (61) we get to the conclusion that in this case we must consider

[TABLE]

which means

[TABLE]

in full agreement with the asymptotic behavior of the EoS at early time for this case.

6.2 The dynamics of the universe

We’ll now study the physical quantities representing the dynamics of our universe due to this case of the model and comparing our results with observations.

6.2.1 The Luminosity Distance

We could now constrain two parameters out of four through asymptotic behavior of the model and the causality condition, where $r$ has to be $2$ while $\alpha$ has to be $3$. Let’s then calculate the luminosity distance function $D_{L}(z)$ and use the observed data of this function to constrain the remaining two parameters of the model. We start by using relation (14) for $r=2$ in (13), so that we have

[TABLE]

which in turn gives

[TABLE]

This has the solution

[TABLE]

And let’s assume that at some time $t=t_{0}$ we have $a=a_{0}$ and $H=H_{0}$. This gives

[TABLE]

So that

[TABLE]

In turn, we can write $H(z)$ as

[TABLE]

The luminosity distance $D_{L}$, Eqs (36), is then calculated using the proper distance which due to (37) will be given by

[TABLE]

Again, the integration of relation (70) can not be solved analytically so that we use numerical methods. The apparent magnitude of Conley et al. data set, [61], is plotted in Fig.5 together with our results. Results from $\Lambda CDM$ model are also shown for comparison.

applying tha $\chi^{2}-$statistics, relation (39), We obtained a best fit of this model case at $s=-2.1$ and $\beta=1.64\times 10^{4}$ for which ${\chi^{2}}_{red}=1.23$. However, the $Q-$value of this result is $1.05\times 10^{-4}$. This $Q-$value in fact tells us that data are not consistent to be drawn from this model case, and so this case will be ruled out.

7 Conclusions

A barotropic fluid model is considered which initially can be described as a correction to the dark energy regime by an an EoS that has two power laws asymptotes. This enables interpolation between the DM and DE equations of state, and so guarantees a smooth phase transition from matter era to DE era. Accordingly, our model is a dark fluid model provides a unification for the two dark sectors in one EoS that has the advantage of the general asymptotic solution. This means that it can describe dark matter at early times while describing dark energy at late times, while adopts a general EoS for DE which enables the cosmological constant as a special case.

Studying the asymptotic behavior of the equation of state constrains the parameters of the equation and restricts the model. It also clarifies the cases that can be physically acceptable as dark fluid which are found to be two cases. On the other hand, studying the phase diagram of the model for each of these two cases shows that our model has a unique solution for times $\tau\succ 0$, evolves to a de Sitter space and has no future finite time singularities.

Causality condition further constrains the parameters of the model, so that at the end, only two out of four parameters are free. These two parameters are constrained using the Hubble diagram of SNe Ia supernovae, where we present an exact solution for the red shift dependence of the luminosity distance in each model, or model case, avoiding the approximate formulas that are accurate only at small red shifts.

The $\chi^{2}-$statistics shows that the $+ve\;s$ case can fit the data well. Accordingly, We studied the energy density evolution for this model case which ensures the smooth transition between the two dark sectors. Furthermore, we examined the evolution of the EoS parameter $\omega$ that again ensures its two extremes of [math] and $-1$ for early and late times. We also Calculated the deceleration parameter $q$ and examined its evolution. Calculations resulted in a negative $q$ with the present day value of $-0.550$, while assures evolution towards more negativity to the future as the Hubble parameter decreases, a clear evidence for accelerated expansion in agreement with observations. On the other hand, the $-ve\;s$ case failed to fit the data and so will be ruled out.

Due to the above mentioned results, our model with a $+ve\;s$ parameter can be considered as a dark fluid model unifying the two dark sectors of the universe and succeeded in describing the evolution of the universe.

8 Aknowledgment

The author would like to thank Dr. Adel Awad for valuable discussions and guidance.