Couplings between holographic dark energy and dark matter
Abstract
We consider the interaction between dark matter and dark energy in the framework of holographic dark energy, and propose a natural and physically plausible form of interaction, in which the interacting term is proportional to the product of the powers of the dark matter and dark energy densities. We investigate the cosmic evolution in such models. The impact of the coupling on the dark matter and dark energy components may be asymmetric. While the dark energy decouples from the dark matter at late time, just as other components of the cosmic fluid become decoupled as the universe expands, interestingly, the dark matter may actually become coupled to the dark energy at late time. We shall call such a phenomenon as incoupling. We use the latest type Ia supernovae data from the SCP team, baryon acoustics oscillation data from SDSS and 2dF surveys, and the position of the first peak of the CMB angular power spectrum to constrain the model. We find that the interaction term which is proportional to the first power product of the dark energy and dark matter densities gives an excellent fit to the current data.
I Introduction
A number of astronomical observations, such as type Ia supernovae (SNIa) [[1]], the cosmic microwave background (CMB) anisotropy [[2]], and the large scale structure (LSS) [[3]], have shown that the expansion of our universe is accelerating. The acceleration of the Universe strongly indicates the existence of a mysterious exotic matter, named dark energy, which has negative pressure and takes the largest proportion of the total density in the Universe. The ultimate fate of our Universe is to be determined by the dark energy, yet we know little about its properties. Furthermore, while both the dark matter and the dark energy is commonly thought to interact weakly with ordinary matter, it is not immediately clear whether there is some interaction between the two, or if the two are in some way linked.
In this paper we consider the dark matterdark energy interaction in the framework of the holographic dark energy (HDE) model [[4], [5], [6], [7]]. The basic ideas of the HDE model are based on the holographic principle [[8]], which was inspired by the Bekenstein entropy bound of black holes[[9]]. In this model, the dark energy is related the quantum zeropoint energy density , the total energy of the whole system with size should not exceed the mass of a black hole of the same size, thus we have . The largest infrared cutoff is chosen by saturating the inequality so that the HDE density is
(1) 
where is a numerical constant, and is the reduced Planck mass. If we take as the size of the current Universe, for instance the Hubble scale , then the dark energy density estimated with Eq.(1) will be close to the observed value. However, Hsu [[6]] pointed out that this yields a wrong equation of state for dark energy. Li [[7]] subsequently proposed that the future event horizon might be chosen as the IR cutoff , and viable HDE models can then be constructed. As demonstrated in Ref. [[7]], the HDE scenario may help to answer the fine tuning problem and the coincidence problem of dark energy.
Various aspects of the HDE model and its extensions have been investigated subsequently. For example, it has been generalized to the case of nonflat Universe [[10]], and it was studied from the perspective of the anthropic principle [[11]]. The holographic dark energy was also considered for driving the inflation, and its resulting power spectrum index calculated [[12]]. Ref. [[13]] defined a new future event horizon as the infrared cutoff and found that the coincidence problem could be alleviated to some extent. Ref. [[14]] investigated the possible form of the cutoff scale for the solution of coincidence problem. Inspired by the HDE model, Ref. [[15]] proposed the Ricci scalar curvature radius as . The model was reconstructed with a scalar field [[16]]. Observational constraints have also been studied [[17], [18]]. The growth of density perturbations in the holographic dark energy models has been studied in [[19]].
Interactions between the DE and DM components were considered for the HDE model [[20], [21]]. Ref.[[20]] suggested that for the case, with interactions between the two dark components, the dark energy could acquire an equation of state in agreement with observations, thus partly solving the noacceleration problem noted in Ref.[[7]] for interactionfree models. Ref.[[21]] generalized the discussion in Ref.[[20]] to the case of a future event horizon. In that paper the proposed form of interaction is , which could also lead to a negative dark energy equation of state.
However, Ref.[[22]] pointed out that the interaction terms such as , , all lead to unphysical results, because the dark matter density could be driven to a negative value in these models.
Here we propose a form of interaction which we considered natural and physically plausible. In general, when two types of matter interact with each other, the interaction should be proportional to the product of the powers of the number density of both. This is most easily understood in the following way: if one component does not exist, then there certainly would be no interaction at all; with its amount increasing, the interaction rate should increase. Thus, at least in the diffuse limit, the interaction rate should be proportional to some powers of the density (the relation could become more complicated when the densities become high). For instance, this is the case for ordinary twobody chemical reactions, namely
(2) 
i.e. the reaction rate is proportional to the product of number densities of the two particle species.
However, besides the interaction in the usual sense (“scattering”), one may also consider a more general case where the dark energy is converted to dark matter or vice versa by “annihilation” or “decay”. In such a case, the “interaction term” or more precisely the energy exchange term would depend only on one of the densities and not on the other.
In general, we consider the interaction term to be proportional to some powers of the density,
(3) 
(4) 
in which and are the energy densities of the matter and HDE respectively, and is a parameter with mass dimension . The power index is assumed to satisfy
where if one of them is 0, we will have the case of “decay” (if the other is 1) or “annihilation” (if the other is 2). If , the interaction suggests that dark energy is converted into dark matter, while suggests the inverse process. This form of interaction was considered in the study of coupled quintessence model [[23]]. However, not much research on this form of interaction was carried out to understand how such a coupling term would affect the evolution of the Universe.
In the next section, we present the evolution equations for the interacting holographic dark energy model (details of our derivation are given in the appendix). In Sec. 3, we study the dynamical behavior of the model. In Sec.4 we explore the parameter space and find the bestfit and allowed region of . Finally, concluding remarks are given in Sec. 5.
Ii The model
Now let us consider a spatially flat FriedmannRobertsonWalker (FRW) Universe filled with matter component , HDE component and the radiation . For simplicity we do not distinguish dark matter and ordinary matter below but incorporate both in the matter component . It is certainly possible or even likely that dark energy only couples to dark matter. However, in that case the baryons simply evolve as where is the scale factor, the baryon component would then become dynamically irrelevant as the Universe expands, so it only affects the details of the model but not its qualitative behavior, hence we will neglect this effect in the present paper. The Friedmann equation reads
(5) 
where is the energy density of the radiation which satisfies . Thus, Eqs. (1), (3), (4), (5) are the dynamic equations which uniquely determine the evolution of and . In Eq.(1), could be Hubble horizon , particle horizon and future event horizon, i.e.
(6) 
or
(7) 
We define the following dimensionless quantities: ,, , where is the current critical density, and , and introduce
(8) 
The dimensionless parameter determines the characteristic strength of the interaction.
The equations of motion of the HDE model for the particle (event) horizon cases are
(9) 
in which the upper(lower) sign represents the case of particle(future event) horizon (c.f. the appendix for details of derivation). It is easy to verify that if in Eq.(1) equals the particle horizon or Hubble horizon , the model cannot fit the current observational data, because current data prefer the dark energy with equation of state near , while in the above two cases the equation of state will be always greater than (see also the analysis in Ref.[[7]]). Therefore, in the following we will focus on the interacting model with (see Eq.(9)).
Iii Dynamical behaviors
In this section we study the dynamical behavior in models with interaction between DE and DM, which may give rise to some interesting behaviors. First we consider the case where are fixed to be 1, then consider the more general cases.
As a typical example we consider the model with , which happens to be the model which bestfits the current data ( for details of the fitting procedure see next section). First we plot the evolution of densities of different components in Fig.1. While one might think that the coupling between DE and DM might help resolve the coincidence problem, we find that for these particular parameter values the evolution of different densities is rather similar to that of the standard model, i.e. the DM and DE densities could still differ by orders of magnitude, so it does not seem to help much in this case. Rather, as in other HDE models, the fine tuning problem could be alleviated: the density of the dark energy is associated with the scale of horizon rather than the Planck scale, due to the holographic principle. We have also investigated cases with larger coupling constant . We find that indeed, for the strong coupling case (), the dark matter and dark energy could be made to evolve with comparable density for an extended range of redshifts, thus in some sense solving the “coincidence problem”. This is perhaps also true for any strongly coupled dark matterdark energy model. However, for strong coupling the fit to observation is not good.
To investigate how the dynamics of the components are affected by coupling
with the other component, one may compare the expansion term () with the interaction term . If the interaction term is greater than or
comparable with the expansion term of both of the components, the
two components are strongly coupled and evolve together. As
the Universe expands, typically the density of the component would drop as , the interaction term often drops much more rapidly
than the expansion term. When the expansion term of a component is
much greater than the interaction term, that component evolves
almost freely, and we may say the two components are decoupled
In Fig. 2, we compare the expansion term of the dark energy and the dark matter with the interaction term as a function of redshift. In the left panel, we plot the magnitude of the dark energy expansion term , and in the right panel, we plot the magnitude of the dark matter expansion term , as red solid lines. The interaction term is plotted in both panels for comparison with the expansion terms. In the early Universe, the dark matter takes a greater proportion of the total density, and for dark energy the interaction term is comparable with the expansion term, so the dark energy is moderately coupled to the dark matter. Later, as the density of dark matter is diluted by expansion, the interaction term becomes negligible compared with the expansion term, and the dark energy is decoupled from dark matter. In the right panel of Fig. 2, we plot the dark matter expansion term and interaction term as a function of redshift . We see that throughout the cosmic history, for the dark matter the expansion term is always greater than the interaction term, so the dark matter is at most weakly coupled to dark energy. However, at lower redshifts, the difference becomes larger, and the dark matter become further decoupled from dark energy. Such a decoupling history is quite typical for many interactions.
However, there could also be some quite unexpected and interesting behavior in the presence of dark energy. In Fig. 3, we make plots similar to Fig. 2 but for the model with , and this model is chosen to illustrate such an unusual behavior. In this model, because , the dark energy density actually increases at later time. We see from the left panel of Fig. 3, since initially the dark matter density is greater, that the interaction term is comparable to the expansion term for dark energy. However, later on the dark energy becomes dominant, and the interaction term falls below the expansion term of dark energy, so we may say that dark energy becomes decoupled. However, more interestingly, from the right panel of the figure, we see that initially the interaction term of the dark matter is smaller than the expansion term, so the dark matter is only weakly coupled to dark energy, but at later time, the interaction term is comparable with the expansion term of dark matter, so the DM becomes coupled to DE. This is the reverse of the familiar decoupling process, which we shall call incoupling. This strange and interesting behavior could only occur when a component is coupled to dark energy, whose density does not decrease during the cosmic expansion. As the DE keeps being converted into DM, the DM evolution is greatly affected by DE. This model does not yield a good fit to the observational data, but is only used to illustrate an interesting dynamical behavior in the presence of dark energy.
In Fig.4, we plot the evolution of for three choices of parameter, i.e. the best fit , and , while keeping . The Hubble parameter for the CDM model is also plotted for comparison. We also plot the equation of state of the dark energy for these models in Fig. 5. On the left panel, curves for different but identical are plotted; on the right panel, curves for identical but different are plotted. In all of these models, the phantom divide is either crossed, or will be crossed in the near future. The dark energy density would keep increasing, the interaction with DM would not slow it down, as the dark energy has become decoupled to the DM. Eventually a “Big Rip” would occur. The current value of the equation of state for the bestfit model is around , which is within the CL region of in the WMAP 5year results[[25]], suggesting that this model is consistent with the current observations.
The value of in these models does not significantly affect the fate of the Universe, but it does greatly affect the past evolution. This result is easy to understand; as the dark energy dynamically decouples from dark matter, the interaction become negligible, regardless the precise value of coupling constant. The eventual fate of the Universe in this interacting holographic dark energy model is qualitatively similar to those without interactions, though the details of the evolution would be different. Thus, as in the case of noninteracting holographic models, for the fate of the Universe, leads to phantomlike behavior and Big Rip[[26]], while leads to quintessence behavior. This is very different from the models discussed in Ref. [[17]]. For further discussions on the fate of the Universe and Big Rip singularity in our model, see Appendix B.
We now consider the more general case of . We plot the equation of state in Fig. 6, varying in the left panel while keeping , and varying in the right panel while keeping . The other parameters are kept as before. For , the equation of state for dark energy is slightly more negative. For , at higher redshift increases rapidly, and the coupling term would become very large, the whole system would be dominated by the evolution of matter. However, such a model would not really fit the low redshift observation, so it would not be relevant to our purpose. For large , we can also see that the equation of state of dark energy at high , but the variation is not as rapid as in the case of , and this is because the dark energy component does not increase as much as the dark matter component. Still, this makes a poor fit to the observational data.
Anyway, for the case , Eq.(29) is reduced to
(10) 
Substituting this equation into Eq.(17) and making use of Eq.(14), we obtain the following approximate equation for
(11) 
for which the solution will simply be
(12) 
where is a constant. This scaling behavior of dark energy density (if ) is very similar to that of the Chaplygin gas model [[27]], which at the high redshift also has the scaling behavior as while at the low redshift has the constant magnitude. However, in the present case it does not provide a good fit to the observations.
Iv Fitting the model
Now we fit the model parameters with the current observational data. We use three kinds of observations; the first one is type Ia supernovae (SNe Ia) which serves as the standard candle; the second is the baryon acoustic oscillation (BAO) in large scale structure of galaxy distribution, the third is the cosmic microwave background anisotropies, the latter two provide standard rulers at different redshifts. For the flat geometry we are considering, , where are the luminosity distance, angular diameter distance and comoving coordinate distance, respectively, and is given by which could be calculated for each model by solving Eq. (9) or Eqs. (30)(31).
For the SNe Ia, we use the data set published by Supernova Cosmology Project (SCP) team recently [[28]]. This data set contains 307 selected SNe Ia which includes several current widely used SNe Ia data sets, such as those collected by the Hubble Space Telescope (HST) [[29], [30]], the SuperNova Legacy Survey (SNLS) [[31]] and the Equation of State: SupErNovae trace Cosmic Expansion (ESSENCE) surveys[[32]]. Using the same analysis procedure and improved selection approach, all of the sub data sets are analyzed to obtain a consistent and highquality “Union’ SNe Ia data set, which provides tighter and more reliable constraints for our models. For the BAO data, we use the data given in Percival et al. [[33]] and WMAP team [[25]]. We use the quantity which is determined by the Sloan Digital Sky Survey (SDSS) () and the Two Degree Field Galaxy Redshift Survey (2dFGRS) () in our constraints, where is the comoving sound horizon size at the decoupling epoch, and is the effective distance defined by Eisenstein et al. [[34]]. Since this quantity is just a ratio of two “standard rulers”, it could give a reliable measurement of the expansion history of the Universe. For the CMB, we use the position of the first peak of the Cosmic Microwave Background (CMB) angular power spectrum . The details of this technique were described in our earlier paper Ref. [[35]].
We use the Markov Chain Monte Carlo (MCMC) method to constrain the parameters with the combined data set of SN Ia+BAO+CMB (for details of our implementation of this technique, see Ref. [[36]]). In the present model, the relevant parameters are the matter density , Hubble constant , and the holographic dark energy parameter . We marginalize over and . For , we limit their range to . As discussed in the last section, greater values of would not yield a good fit to observation.
For the fixed case, as can be seen in Fig. 7, the bestfit value and error on are
(13) 
The behavior of the best fit model has been discussed as the example given in the previous section. The models within the contour also exhibit similar behavior.
Now we let be free parameters, the PDF of and are shown in Fig. 8. The bestfit values for are and . If one demand to be integers for physical reasons, then is still the best fit. Alternatively, a noninteger power may rise if the interaction is complicated. The best fit for and are and , as shown in Fig. 9, it is not too different from the case.
V Conclusion
The holographic dark energy (HDE) model is a dynamical model of dark energy, motivated by the holographic principle, in which the event horizon of the spacetime plays an important role in the dynamics. The microphysical basis of the model is not well understood, but the dynamics of cosmic expansion could be studied. Observationally, the HDE model is distinguishable from the CDM, as its equation of state varies with time. The model could also be distinguished from other dynamical dark energy models. For example, in the HDE model the dark energy behaves as the Chaplygin gas at high redshift if (see the discussion at the end of Sec.3), but for the low redshift region, the equation of state of the Chaplypin gas approaches a constant, while for this model it continues to evolve. In properly constructed HDE models such as the ones discussed in this paper, the dark energy component is subdominant at high redshifts compared with matter or radiation, so it would not significantly affect the CMB or large scale structure power spectrum, though at low redshift it could affect the growth of structure in detail. We plan to investigate the details of the structure growth in subsequent studies.
In the framework of the holographic dark energy model, we proposed a form of interaction between dark energy and dark matter which we considered to be natural and physically plausible. That is, the interaction rate is proportional to the product of powers of the dark energy and dark matter densities . For this type of interaction, we find that for the weak coupling case (the dimensionless coupling constant , see Eq. (8)), the ultimate fate of the Universe is similar to the case of noninteracting holographic dark energy, namely for the equation of state of the dark energy will cross , and behave like a phantom in the future, and the Universe ends in a Big Rip. We also discovered the interesting phenomenon of incoupling, which occurs only when the dark matter is coupled to dark energy. This is the opposite of the familiar process of decoupling, the dark matter becomes increasingly coupled to the dark energy as Universe expands. We then use the currently available observational data to constrain the parameters of this model. The data set we used include 307 selected SNe Ia published by the SCP team [[28]], baryon acoustic oscillation from SDSS and 2dFGRS [[33]], as well as the position of the first peak of the CMB angular power spectrum [[25]]. We find that the bestfit form of interaction is . If one requires to be integers for reasons of fundamental physics, then the most favored interaction rate is proportional to the product of the density of dark matter and dark energy.
Acknowledgements.
Our MCMC chain computation was performed on the Supercomputing Center of the Chinese Academy of Sciences and the Shanghai Supercomputing Center. This work is supported by the Chinese Academy of Sciences under grant KJCX3SYWN2, by the Ministry of Science and Technology National Basic Science program (project 973) under grant No.2007CB815401, and by the National Science Foundation of China under the Distinguished Young Scholar Grant 10525314, and the Key project Grants No.10503010.Appendix A dynamical equations
We give the detailed derivation of the dynamical equations of the interacting holographic dark energy model here. Using the dimensionless notations given in Sec.2, we have
(14) 
(15) 
(16) 
(17) 
(18) 
So
(19) 
(20) 
By substituting Eq.(19) into Eq. (17), we have
(21) 
If is taken to be the particle horizon (Eq.(6)) or future event horizon (Eq.(7)), the derivative of is
(22) 
where the upper (lower) sign represents particle horizon (future event horizon). By substituting Eq. (22) and (18) into Eq. (21), we have the equation of state in the particle horizon (upper sign) and future event horizon (lower sign) cases
(23) 
Substituting Eq. (23) into (20), we have
(24) 
Taking the derivative of Eq. (14) with respect to , we have
(25) 
Substituting Eq. (16) into Eq. (25), we obtain
(26) 
Subtracting equations (26) from (26), we obtain the following simple dynamical equation:
(27) 
Eliminate in Eq. (24) using Eq.(14),
(28) 
and
(29) 
To make Eq. (27) and (28) easy to solve numerically, we change the variable to ,
(30) 
(31) 
where
(32) 
and the upper (lower) sign in Eq. (30) and (31) corresponds to the model with particle (future event) horizon. Thus, Eq. (30) and (31) are the complete system of equations. Eqs.(29)(31) reproduce Eqs. (9). In the limit and the equations reduce to the noninteracting holographic dark energy model [[7], [37]]. The evolution of event horizon size and the corresponding dark energy density are shown in Fig. 10. As , the size of the event horizon shrinks to zero, and the dark energy density diverges, i.e. we are approaching a Big Rip singularity.
Appendix B Big Rip Singularity
There have been extensive discussions on the future singularity of the dark energy models, particularly for the case where the future Universe is dominated by the phantom dark energy. Ref. [[38], [39]] investigated the appearance of the future singularity, and showed that the singularities could be classified into four types. The Big Rip is however, the only relevant singularity in the holographic dark energy model [[39], [40]]. Ref [[41]] showed that the coupling between dark matter and dark energy may help to solve the coincidence problem, but not the Big Rip (Type I) singularity issue. In particular, [[40]] discussed the singularity problem in the holographic dark energy model, and proved that:

if holographic dark energy there exists a de Sitter solution to the dynamic equation system of the HDE. But if the de Sitter solution does not exist.

For the de Sitter solution (), one should consider the back reaction of the quantum effects near the singularity. The contribution of the conformal anomaly as a back reaction near singularity could indeed moderate or even prevent the future singularity (Eqs.(53)(61) in [[40]]).
In our interacting HDE model, the bestfit parameter , and for such model the fate of the Universe is a Big Rip in the future. In the following, we investigate whether this conclusion is affected by the back reaction of quantum effects.
Since the future of the Universe is dominated by dark energy, we will neglect all of the terms which contain and , as they will become subdominant when compared with the dark energy density The dynamical equations then simplify to
(33) 
(34) 
(35) 
where we only consider the case for future event horizon. From Eqs. (33)(35), the asymptotical behavior of , and can be obtained easily:
(36) 
(37) 
(38) 
where and is the fractional HDE density taken at the present day. From Eqs. (38), the derivative terms are
(39) 
(40) 
(41) 
The conformal anomaly (back reaction term) has the following form [[40]]:
(42) 
Thus in the HDE model
(43) 
where
(44)  
With Eqs.(43) and (44), the energy density and pressure of the conformal anomaly are [[40]],
(45)  
(46) 
Under this back reaction, the expansion rate is
(47)  
Since the energy density of back reaction is positive (, the function therefore, increase monotonically in the future if Thus, we prove that for the case, the back reaction effect from quantum corrections cannot prohibit the occurrence of the Big Rip.
Footnotes
 Note here that in the literatures the word “decoupling” may refer to a number of different phenomena. For the decoupling of WIMP dark matter from the radiationbaryon fluid, what is discussed above (the evolution of density) is usually called “thermal decoupling”. There is also “kinetic decoupling”, which refers to decoupling of temperature (see, e.g., [[24]]). Admittedly, this nomenclature is somewhat confusing, since one might naively think that “thermal decoupling” refers to decoupling of temperature, and “kinetic decoupling” refers to decoupling of density evolution! Below we shall refer to the decouplings discussed here as “dynamical decoupling”, the terms “thermal decoupling” and “kinetic decoupling” are not used in this paper.
References
 A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) [astroph/9805201]; S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999) [astroph/9812133].
 D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) [astroph/0302209]; D. N. Spergel et al., Astrophys. J. Suppl. 170, 377 (2007) [astroph/0603449].
 M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) [astroph/0310723]; K. Abazajian et al. [SDSS Collaboration], Astron. J. 128, 502 (2004) [astroph/0403325]; K. Abazajian et al. [SDSS Collaboration], Astron. J. 129, 1755 (2005) [astroph/0410239].
 A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Rev. Lett. 82, 4971 (1999) [hepth/9803132].
 P. Horava and D. Minic, Phys. Rev. Lett. 85, 1610 (2000) [hepth/0001145]; S. D. Thomas, Phys. Rev. Lett. 89, 081301 (2002).
 S. D. H. Hsu, Phys. Lett. B 594, 13 (2004) [hepth/0403052].
 M. Li, Phys. Lett. B 603, 1 (2004) [hepth/0403127].
 G. ’t Hooft, grqc/9310026; L. Susskind, J. Math. Phys. 36, 6377 (1995) [hepth/9409089].
 J. D. Bekenstein, Phys. Rev. D 7, 2333 (1973); J. D. Bekenstein, Phys. Rev. D 9, 3292 (1974); J. D. Bekenstein, Phys. Rev. D 23, 287 (1981); J. D. Bekenstein, Phys. Rev. D 49, 1912 (1994); S. W. Hawking, Commun. Math. Phys. 43, 199 (1975); S. W. Hawking, Phys. Rev. D 13, 191 (1976).
 Q. G. Huang and M. Li, JCAP 0408, 013 (2004) [astroph/0404229].
 Q. G. Huang and M. Li, JCAP 0503, 001 (2005) [hepth/0410095]
 K. Enqvist and M. S. Sloth, Phys. Rev. Lett. 93, 221302 (2004) [hepth/0406019]; B. Chen, M. Li and Y. Wang, Nucl. Phys. B 774, 256 (2007) [arXiv:astroph/0611623].
 S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 38, 1285 (2006) [hepth/0506212];
 B. Hu and Y. Ling, Phys. Rev. D 73, 123510 (2006) [hepth/0601093].
 C. J. Gao, X. L. Chen, Y. G. Shen, arxiv:0712.1394 [astroph].
 X. Zhang, Phys. Lett. B 648, 1 (2007) [astroph/0604484]; Y. Z. Ma and X. Zhang, Phys. Lett. B 661, 239 (2008) 0709.1517 [astroph].
 Y. Z. Ma, Y. Gong and X. Chen, Eur. Phys. J. C 60 (2009) 303, 0711.1641 [astroph].
 X. Zhang, Int. J. Mod. Phys. D 14, 1597 (2005) [astroph/0504586].
 K. Y. Kim, H. W. Lee, Y. S. Myung, arixv:0805.3941 [grqc]
 D. Pavon and W. Zimdahl, Phys. Lett. B 628, 206 (2005) [grqc/0505020].
 B. Wang, Y. Gong and E. Abdalla, Phys. Lett. B 624, 141 (2005) [hepth/0506069]; C. Feng, B. Wang, Y. Gong and R. K. Su, 0706.4033 [astroph].
 M. Li, C. Lin and Y. Wang, JCAP 05, 023 (2008), 0801.1407 [astroph].
 G. Mangano, G. Miele and V. Pettorino, Mod. Phys. Lett. A 18, 831 (2003), [astroph/0212518].
 X. Chen, M. Kamionkowski, and X. Zhang, Phys. Rev.D 64, 021302 (2001), [astroph/0103452].
 E. Komastu et al., 0803.0547 [astroph].
 R.R. Caldwell, Phys. Lett. B 45, 23 (2002), [astroph/9908168].
 N. Bilic, G. B. Tupper, R. D. Viollier, astroph/0207423; N. Bilic, G. B. Tupper, R. D. Viollier, Phys. Lett. B 535 (2002) 1721, [astroph/0111325]; M. C. Bento, O. Bertolami, A. A. Sen, Phys. Rev. D 66 (2002) 043507, [grqc/0202064].
 M. Kowalski et al., arXiv: 0804.4142 [astroph].
 A. G. Riess et al., Astrophys. J. 607, 665687, (2004), [astroph/0402512].
 A. G. Riess et al., Astrophys. J. 659, 98121, (2007), [astroph/0611572].
 P. Astier et al., Astron. Astrophys. 447, 3148 (2006), [astroph/0510447].
 W. M. WoodVasey et al., arXiv: astroph/0701041.
 W. J. Percival et al., MNRAS 381, 1053, (2007), arxiv:0705.3323 [astroph].
 D. J. Eisenstein et al., Astrophys. J. 633, 560574 (2005), [astroph/0501171].
 Y. Gong and X. Chen, Phys. Rev. D 77, 103511 (2008), arXiv:0802.2296 [astroph].
 Y. Gong and X. Chen, Phys. Rev. D 76,123007 (2007), arXiv:0708.2977 [astroph].
 Q. Wu, Y. Gong, A. Zhong and J. S. Alcaniz, arXiv:0705.1006 [astroph].
 E. Elizalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70 (2004) 043539, hepth/0405034.
 S. Nojiri, S. D. Odintsov and S Tsujikawa, Phys. Rev. D 71 (2005) 063004, hepth/0501025.
 E. Elizalde, S. Nojiri, S. D. Odintsov and P. Wang, Phys. Rev. D 71 (2005) 103504, hepth/0502082.
 S. Nojiri and S. D. Odintsov, arXiv:0911.2781 [hepth].