Abstract

The bulk viscosity is introduced to model the current observational cosmos and the unified dark matter (dark matter and dark energy). The viscous unified model assumes that the universe is filled with a single nonperfect fluid with a bulk viscosity. We review the general framework for the viscous cosmology model building finding that the Hubble parameter has possessed a direct connection with the bulk viscosity coefficient. For the concrete form of the bulk viscosity, the Hubble parameter which has the scaling relation with the redshift can be obtained generally. We discuss two viscosity models and the cosmological evolution to which they lead in detail. Using SNe Ia data, CMB shift, as well as BAO observations, the viscosity model parameters can be fitted very well. We also briefly review the fitting method in this paper.

1. Introduction

Both the dark matter problem and the cosmic acceleration problem challenge physicists' understanding of the universe. In the standard ΛCDM model, two mixed fluids dark matter and dark energy fluid, are assumed. These two fluids influence the cosmic evolution separately. However, present gravitational probe does not have the ability to differentiate these two fluids. This is the dark degeneracy problem [1, 2]. It is reasonable to model dark matter and dark energy with single fluid or single-field assumption. Some unified models have been proposed to detect the possibility of this unified assumption, like unified dark fluid model [37], which assumes the single-fluid equation of state; Chaplygin gas model and generalized Chaplygin gas model [812], which discuss the cosmology consequences of an exotic equation of state; scalar field method [1315].

The introduction of viscosity into cosmology has been investigated from different viewpoints [16, 17]. There are some recent developments like dark energy model [1821], the cosmic singularity [22]. In this paper, we give a brief introduction to unify dark matter and dark energy with viscosity medium. In such models, the universe is assumed to be filled with viscous single fluid [2329]. The cosmic density is not separated as dark energy part and dark matter part. The bulk viscosity contributes to the cosmic pressure and plays the role of accelerating the universe. After considering the bulk viscosity, the cosmic pressure can be written as𝑝=(𝛾1)𝜌3𝜁𝐻,(1) where 𝛾 parameterizes the equation of state. Generally the form of bulk viscosity is chosen as a time-dependent function. In [3035], a density-dependent viscosity 𝜁=𝛼𝜌𝑚 coefficient is investigated extensively. For modeling the unified dark matter and dark energy, it is often assumed that the parameter 𝛾=1, that the pressure of the viscosity fluid is zero, and the viscosity term contributes an effective pressure. There arise some problems here. From the observational results [36], the cosmic density is nearly equals to the cosmic pressure. In the viscosity model, the viscosity term dominates the cosmic pressure and surpasses the pressure contributions from other cosmic matter constitutions, which contradicts the traditional fluid theory. The works in [37, 38] propose nonstandard interaction mechanism to solve this problem. Obviously, it is important to build solid foundation for the research of the viscous cosmology.

Equation of state 𝑤<1 lies in the phantom region. It is shown that cosmology models with such equation of state possess the so-called the future singularity that is called the Big Rip [39]. The larger viscosity model parameter space can help to solve the cosmic singularity problem and produce different kinds of evolution modes of the future universe; for more details [28].

The rest of this paper is organized as follows. In Section 2, general framework of the viscosity model will be reviewed. In Section 3, we discuss the modeling of the unified model with viscosity. In this section, two concrete models are analyzed. In Section 4, data fitting method is introduced briefly.

2. General Framework

We consider the standard Friedmann-Robertson-Walker metricd𝑠2=d𝑡2+𝑎(𝑡)2d𝑟21𝑘𝑟2+𝑟2dΩ2.(2) For the sake of simplicity, we choose the flat geometry 𝑘=0, which is also favored by the update result of the cosmic background radiation measurement.

The general stress-energy-momentum tensor is𝑇𝜇𝜈=(𝜌+𝑝)𝑈𝜇𝑈𝜈+𝑝𝑔𝜇𝜈𝜁𝜃𝜇𝜈,(3) where 𝜁 is the bulk viscosity. The expansion factor 𝜃 is defined by 𝜃=𝑈𝜇;𝜇=3(̇𝑎/𝑎), and the projection tensor 𝜇𝜈𝑔𝜇𝜈+𝑈𝜇𝑈𝜈. In the comoving coordinates, the four velocity 𝑈𝜇=(1,0,0,0). We do not specify the concrete form of 𝜁 in this section. Generally speaking, 𝜁 is a quantity evolving with time 𝑡 or the scale factor 𝑎(𝑡). We will see below that nontrivial and more complicated 𝜁 can produce different results especially useful for the late universe modeling.

From the usual Einstein equation𝐺𝜇𝜈=𝑅𝜇𝜈12𝑅𝑔𝜇𝜈=8𝜋𝐺𝑇𝜇𝜈,(4) we obtain the two equations which we call the modified Friedmann equations: ̇𝑎2𝑎2=8𝜋𝐺3𝜌,̈𝑎𝑎=4𝜋𝐺3(𝜌+3̃𝑝),(5) where ̃𝑝 is an effective pressure and ̃𝑝=𝑝𝜁𝜃.

The covariant conservation equation 𝑇0𝜇;𝜇=0 yieldṡ𝜌+(𝜌+̃𝑝)𝜃=0.(6) The existence of a bulk viscosity contributes a modification to the pressure 𝑝, thus we see that the Friedmann equation and the covariant conservation equation are invariant under the transformation𝑝̃𝑝=𝑝𝜁𝜃.(7)

The covariant energy conservation equation becomeṡ𝜌+(𝜌𝜁𝜃)𝜃=0.(8) We define the dimensionless Hubble parameter here as2=𝐻2𝐻20=𝜌𝜌cr,(9) where 𝜌cr=3𝐻20/8𝜋𝐺 is the critical density now. Using the dimensionless Hubble parameter, (8) can be transformed as1𝐻0d2d𝑡+33=9𝜆2,(10) where the bulk viscosity is redefined as 𝜆=𝐻0𝜁/𝜌cr. Through the simple relation between scale factor 𝑎(𝑡) and the redshift 𝑧1d𝑡=𝑎𝐻d𝑎,(11) we transform (10) into a differential equation with respect to the scale factor 𝑎(𝑡) as follows:d𝐻+3d𝑎2𝑎𝐻=3𝜁2𝑎.(12) Solving this equation, we obtain an integral form of 𝐻(𝑎) as𝐻(𝑎)=𝐶1𝑎3/2+3𝜁32𝑎exp32𝑎d𝑎d𝑎exp2𝑎d𝑎.(13) Different forms of viscosity can be used here to make this integral calculable, numerically or exactly.

3. Unified Single Fluid

3.1. Redshift-Dependent Model

In [40], the authors assume that the bulk viscosity takes the form of a Hubble parameter-dependent function. A redshift-dependent viscosity is proposed in [27]. This bulk viscosity is a combination of a constant and a scaling relation term9𝜆=𝜆0+𝜆1(1+𝑧)𝑛,(14) where 𝑛 is an integer and 𝜆0 and 𝜆1 are two constants, which will be fitted from the observational datasets.

After taking account of this ansatz, the integration is easy to work out. We get2(𝑧)=𝜆22(1+𝑧)3+23𝜆0𝜆2(1+𝑧)1.52𝜆0𝜆13(2𝑛3)(1+𝑧)𝑛+𝜆21(2𝑛3)2(1+𝑧)2𝑛2𝜆1𝜆22𝑛3(1+𝑧)𝑛+1.5+𝜆209.(15) Since we have assumed the spatial flatness of the universe, the consistency condition requires that (0)=1. Thus this sets a constraint on the model parameters as𝜆03=1𝜆2+𝜆12𝑛3.(16) We remind the readers that we make the single-fluid assumption above, and we do not concretely specify the constitutions of the cosmic density 𝜌. In this single-fluid model, values of model parameters 𝜆0, 𝜆1, and 𝜆2 will be fitted, and their meanings are not explained. But when we compare it with two-fluid model, having in mind that the universe is filled with dark matter and dark energy fluid, more constraints can be added. The solution consists of terms with different scaling relation. The first term has the form of 𝐶(1+𝑧)3, which has the same evolution behavior as that of the cold dark matter. Their simplicity leads us to parameter 𝜆2 correspond to dark matter ratio Ω𝑚𝜆22=Ω𝑚.(17) This identity can help us utilize more data to constrain the viscosity model. The result is also consistent with that obtained from the standard model (ΛCDM). The shift parameter [41, 42] and the distance parameter 𝒜 are defined asΩ𝑚𝑧0d𝑧(𝑧),(18)𝒜Ω𝑚𝑧𝑏1/31𝑧𝑏𝑧𝑏0d𝑧(𝑧)2/3,(19) respectively. Both of them are dependent on dark matter ratio Ω𝑚, and in the joint statistical analysis, they provide strong constraint on Ω𝑚.

3.2. Effective Equation of State Model

Another viscosity model reviewed here is proposed in [25], where a general form of time-dependent viscosity is discussed:𝜁=𝜁0+𝜁1̇𝑎𝑎+𝜁2̈𝑎̇𝑎.(20) An interesting feature of this model is its effective equivalence to the following equation of state:𝑝=(𝛾1)𝜌+𝑝0+𝑤𝐻𝐻+𝑤𝐻2𝐻2+𝑤d𝐻̇𝐻,(21) where 𝑝0, 𝑤𝐻, 𝑤𝐻2, and 𝑤d are free parameters. The corresponding between two groups of coefficients are 𝑤𝐻=3𝜁0,𝑤𝐻2𝜁=31+𝜁2,𝑤d𝐻=3𝜁2.(22) The parameterized bulk viscosity combines terms related to the “velocity” ̇𝑎 and “acceleration” ̈𝑎, which can be seen to describe the dynamics of the cosmic nonperfect fluid. After eliminating 𝑝 and 𝜌, a differential equation about the scale factor 𝑎(𝑡) can be obtained as̈𝑎𝑎=𝜅(3𝛾2)/22𝑤/2𝐻2+𝜅2𝑤d𝐻𝜅1+2𝑤d𝐻̇𝑎𝑎2+𝜅2𝑤𝐻𝜅1+2𝑤/2d𝐻̇𝑎𝑎+𝜅2𝑝/20𝜅1+2𝑤/2d𝐻.(23) Another feature of this model is that this differential equation can be solved exactly, and the evolution function of the scale factor 𝑎(𝑡) is definite. This evolution function is especially convenient for discussing the cosmic singularity.

With the initial conditions 𝑎(𝑡0)=𝑎0 and 𝜃(𝑡0)=𝜃0, when ̃𝛾0, the scale factor can be obtained as𝑎(𝑡)=𝑎0121+̃𝛾𝜃0𝑇𝑇𝑇1exp𝑡𝑡021𝑇+1𝑇1+121̃𝛾𝜃0𝑇𝑇+𝑇1exp𝑡𝑡021𝑇1𝑇12/3,̃𝛾(24) where the parameters are redefined as𝜅̃𝛾=𝛾+2𝑤/3𝐻2𝜅1+2𝑤/2d𝐻,1𝑇1=𝜅2𝑤/2𝐻𝜅1+2𝑤/2d𝐻,1𝑇22=𝜅2𝑝/20𝜅1+2𝑤/2d𝐻,1𝑇2=1𝑇21+6̃𝛾𝑇22.(25) From Friedmann equation, 𝜌 can be written as 1𝜌(𝑡)=3𝜅2̃𝛾21+̃𝛾𝜃0𝑇𝑇/𝑇11/𝑇+1/𝑇1exp𝑡𝑡0/𝑇1̃𝛾𝜃0𝑇+𝑇/𝑇11/𝑇1/𝑇11+̃𝛾𝜃0𝑇𝑇/𝑇1exp𝑡𝑡0+/𝑇1̃𝛾𝜃0𝑇+𝑇/𝑇12.(26)The model parameters leave enough space to produce various evolution behaviors, which can be interpreted in different ways. In this review, we emphasze their power to unify dark energy and dark matter with the single-fluid assumption. According to the parameters redefined above and the Friedmann equation, the equation of state can be converted to2𝑝=(̃𝛾1)𝜌3𝜅𝑇12𝜌𝜅2𝑇22.(27) The case ̃𝛾=0 and 𝑇1 corresponds to the ΛCDM. With the aim to unify dark energy and dark matter, the case ̃𝛾=1 and 𝑇2 is especially considered. This case corresponds to a single fluid with constant viscosity. The relation between 𝑝 and 𝜌 can be obtained from the general equation of state above as follows:2𝑝=3𝜅𝑇1𝜌.(28) Therefore, it is straightforward to eliminate 𝑝 from the covariant energy conservation equation and to work out the solution of 𝜌. Using Friedmann equation, 𝐻(𝑧) can be obtained𝐻(𝑧)=𝐻0Ω𝛾(1+𝑧)3/2+1Ω𝛾,(29) where Ω𝛾 is the only one model parameter. Its value can be fitted from SNe Ia observational data.

4. Data Fitting

We review the method to fit the model parameters. More details are illustrated in [43]. The datasets we use are SNe Ia, BAO, and CMB. The 397 Constitution sample [44] combines the Union sample [45] and the low-redshift (𝑧<0.08) sample [46]. The comoving distance 𝑑𝑀 in FRW coordinate is𝑑𝑀=𝑧01𝐻(𝑧)d𝑧.(30) The apparent magnitude which is measured is𝑚𝑀+5log10𝐷𝐿(𝑧),(31) where the dimensionless luminosity 𝐷𝐿𝐻0𝑑𝐿(𝑧) and𝑑𝐿=(1+𝑧)𝑑𝑀(𝑧),(32) where 𝑀 is the absolute magnitude which is believed to be constant for all SNe Ia. In the SNe Ia samples, data are given in terms of the distance modulus 𝜇obs𝑚(𝑧)𝑀obs(𝑧). The 𝜒2 for this procedure is written as𝜒2=𝑛𝑖=1𝜇obs𝑧𝑖𝜇th𝑧𝑖;𝑐𝛼𝜎obs𝑧𝑖2,(33) where 𝜇th means the distance modulus calculated from model with parameters 𝑐𝛼 (𝛼=0,1,2,). Together with the shift parameter and the distance 𝒜, the total 𝜒2total for the joint data analysis is given as𝜒2total=𝜒2+obs𝜎2+𝒜𝒜obs𝜎𝒜2.(34) For the redshift-dependent model, the relation between distance modulus and redshift is plotted in Figure 1. The model calculated value and the Constitution data are compared in the figure.

5. Conclusion

In this review, we discuss three aspects of the viscosity model.(i)There is general framework for viscosity modeling. General form of Hubble parameter is presented. This general form is convenient for comparing different scale factor- (or redshift-) dependent viscosity models. (ii)Two kinds of viscosity models are used to model unified models. (iii)Observation constraint is necessary for model building. We can see that the fitting results are consistent with data. It is prospected that more accurate direct measurements of the Hubble constant will provide a new constraint on cosmological parameters [47].

Especially we focus on its application on modeling the unified dark energy and dark matter.

In the cosmic background level, dynamical analysis can be performed. The statefinder method is useful for discriminating different models [4852]. Compared with ΛCDM model, evolution of the statefinder of the viscosity model is different and can be discriminated easily; more details can be found in [28, 53]. More plentiful and accurate data will improve the power of the statefinder method, which will give enough constraint on the late universe model.

We review the viscosity model which is on the level of zero order. The perturbation analysis and to confront the observed large-scale structure are especially useful for the model building and distinguishing. The model predictions need to be consistent with CMB and LSS data. Some works have investigated the perturbation aspects of the viscosity cosmology model [35, 54]. After choosing properly the corresponding model parameters, the viscosity cosmology model has been shown to have the connection with the Chaplygin gas model. Though the Chaplygin gas model can fit the SNe Ia data well at the zero order, in the perturbation level it is found that the Chaplygin gas model does not behave in a satisfactory way. Whether the viscosity cosmology models could behave well needs further investigation.

Acknowledgment

This work is supported in part by the National Science Foundation of China under No. 11075078.