Cosmological Constraints on Ghost Dark Energy in the Brans-Dicke Theory by Using MCMC Approach

Cosmological Constraints on Ghost Dark Energy in the Brans-Dicke Theory by Using MCMC Approach

Hamzeh Alavirad and Ahmad Sheykhi 222 Institute for Theoretical Physics, Karlsruhe Institute of Technology (KIT), 76128 Karlsruhe, Germany
Physics Department and Biruni Observatory, College of Sciences, Shiraz University, Shiraz 71454, Iran
Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), P. O. Box 55134-441, Maragha, Iran

By using a Markov Chain Monte Carlo simulation, we investigate cosmological constraints on the ghost dark energy (GDE) model in the framework of the Brans-Dicke (BD) theory. A combination of the latest observational data of the cosmic microwave background radiation data from seven-year WMAP, the baryon acoustic oscillation data form the SDSS, the supernovae type Ia data from the Union2 and the X-ray gas mass fraction data from the Chandra X-ray observations of the largest relaxed galaxy clusters are used to perform constraints on GDE in the BD cosmology. In this paper, we consider both flat and non-flat universes together with interaction between dark matter and dark energy. The main cosmological parameters are obtained as: , and . In addition, the Brans-Dicke parameter is estimated as .
Keywords: ghost; dark energy; Brans-Dicke theory; observational constraints.

I Introduction

Accelerating expansion of the Universe Riess:1998cb (); Perlmutter:1998np () can be explained either by a missing energy component usually called “dark energy” (DE) with an exotic equation of state, or by modifying the underlying theory of gravity on large scales. The famous examples of the former approach include scalar field models of DE such as quintessence Wetterich:1987fm (); Ratra:1987rm (), K-essence ArmendarizPicon:2000dh (); ArmendarizPicon:2000ah (), tachyon Padmanabhan:2002cp (); Sen:2002in (), phantom Caldwell:1999ew (); Nojiri:2003vn (); Nojiri:2003jn (), ghost condensate ArkaniHamed:2003uy (); Piazza:2004df (), quintom Elizalde:2004mq (); Nojiri:2005sx (); Anisimov:2005ne (), holographic DE Witten:2000zk (), agegraphic DE Cai:2007us (); Wei:2007ty () and so forth. For a comprehensive review on DE models, see Copeland:2006wr (); Li:2011sd (). The latter approach for explanation of the acceleration expansion is based on the modification of the underlying theory of gravity on large scales such as gravity DeFelice:2010aj () and braneworld scenarios Dvali:2000hr (); Carena:2006yr (); Minamitsuji:2008fz (); Sheykhi:2007wn ().

Among various models of DE, the so called ghost dark energy (GDE) has attracted a lot of interests in recent years. The origin of DE in this model comes from Veneziano ghosts in QCD theory Kawarabayashi:1980dp (); Witten:1979vv (); Veneziano:1979ec (); Rosenzweig:1979ay (). Indeed, the contribution of the ghosts field to the vacuum energy in curved space or time-dependent background can be regarded as a possible candidate for DE Ohta:2010in (); Urban:2009vy (). The magnitude of this vacuum energy is of order , where is the Hubble parameter and is the QCD mass scale. With and , gives the right order of magnitude for the observed dark energy density Ohta:2010in (). The advantages of GDE model compared to other DE models is that it is totally embedded in standard model and general relativity, therefore one needs not to introduce any new parameter, new degree of freedom or to modify gravity. The dynamical behavior of GDE model in flat universe have been studied Cai:2010uf (). The study was also generalized to the universe with spacial curvature Sheykhi:2011xz (). The instability of the GDE model against perturbations was studied in Ebrahimi:2011js (). In Sheykhi:2011fe (); Sheykhi:2011nb () the correspondence between GDE and scalar field models of DE were established. In the presence of bulk viscosity and varying gravitational constant, the GDE model was investigated in Sheykhi:2012zzb (). Other features of the GDE model have been investigated in Refs. RozasFernandez:2011je (); Karami:2012yf (); KhodamMohammadi:2012ww (); Malekjani:2012wc (); Feng:2012gr (). The cosmological constraints on this model have been considered by some authors Cai:2010uf (); Feng:2012gr (); Cai:2012fq ().

Scalar tensor theories have been reconsidered extensively, recently. One important example of the scalar tensor theories is the BD theory of gravity which was introduced by Brans and Dicke in 1961 to incorporate the Mach’s principle in the Einstein’s theory of gravity Brans:1961sx (). This theory also passes the observational tests in the solar system domain Bertotti:2003rm (). In addition, BD theory can be tested by the cosmological observations such as the cosmic microwave background (CMB) and large scale structure (LSS) Chen:1999qh (); Acquaviva:2007mm (); Tsujikawa:2008uc (); Wu:2009zb (); Wu:2013 (). Since the GDE model have a dynamic behavior, it is more reasonable to consider this model in a dynamical framework such as BD theory. It was shown that some features of GDE in BD cosmology differ from Einstein’s gravity Ebrahimi:2011ne (). For example, while the original DE is instable in all range of the parameter spaces in standard cosmology Ebrahimi:2011js (), it leads to a stable phase in BD theory Saaidi:2012vs (). In the framework of BD cosmology, the ghost model of DE has been studied Ebrahimi:2011ne (). It is also of great interest to see whether the GDE model in the framework of the BD theory is compatible with observational data or not.

In this paper, cosmological constraints on GDE in the BD theory (GDEBD) Ebrahimi:2011ne () theory is performed by using the Marko Chain Monte Carlo (MCMC) simulation. The used observational datasets are as follows: cosmic microwave background radiation (CMB) from WMAP7 Komatsu:2010fb (), 557 Union2 dataset of type Ia supernova Amanullah:2010vv (), baryon acoustic oscillation (BAO) from SDSS DR7 Percival:2009xn (), and the cluster X-ray gas mass fraction from the Chandra X-ray observations Allen:2007ue (). To put the constraints, the modified CosmoMC cosmomc () code is used.

The organization of this paper is as follows. In section II, we review the formalism of the GDE in the framework of Brans-Dicke cosmology. In section III the methods which are used in this paper to analyze the data are introduced. Section IV contains the results of the MCMC simulation and we conclude our paper in section V.

Ii Interacting Ghost Dark Energy in the Brans-Dicke Theory in a Non-Flat Universe

Let us first review the formalism of the interacting GDE in the framework of BD theory in a non-flat universe Ebrahimi:2011ne (). The action of the BD theory in the canonical form may be written Arik:2005ir ()


where is the Ricci scalar and is the BD scalar field. Varying the action with respect to the metric and the BD scalar field , yields


where stands for the energy-momentum tensor of the matter fields. The line element of the Friedmann-Robertson-Walker (FRW) universe is


where is the scale factor, and is the curvature parameter with corresponding to open, flat, and closed universes, respectively. Nowadays, there are some evidences in favor of closed universe with a small positive curvature () Bennett:2012zja (). Using metric (4), the field equations (2) and (3) reduce to


where is the Hubble parameter, and are, respectively, the energy density and pressure of DE, and is the pressureless matter density which contains both dark matter (DM) and baryonic matter (BM) densities i.e. where and are the energy densities of dark matter and baryonic matter respectively.

To be more general and because of some observational evidences Bertolami:2007zm (); Tsujikawa:2004dp (), here we propose the case where there is an interaction between GDE and DM. In this case the semi-conservation equations read


where represents the interaction term between dark matter and dark energy and here we assume that the baryonic matter is conserved separately. We assume with being a constant. Such a choice for interacting term implies the the DE and DM component do not conserve separately while the total density is still conserved through


where and .

The ghost energy density is proportional to the Hubble parameter Ohta:2010in ()


where is roughly of order and is QCD mass scale. Taking into account the fact that and for the present time, this gives the right order of magnitude for the ghost energy density Ohta:2010in ().

Since the system of equations (5-7) is not closed, we still have another degree of freedom in analyzing the set of equations. As usual we assume the BD scalar field has a power law relation versus the scale factor,


An interesting case is when is small whereas is high so that the product results of order unity Banerjee:2007zd (); Sheykhi:2009dz (). In section IV we will consider the condition for constraining the model by observational data. This is interesting because local astronomical experiments set a very high lower bound on Will (); in particular, the Cassini experiment implies that Bertotti:2003rm (); Acquaviva:2007mm (). Now we take the time derivative of relation (13). We arrive at


Combining Eqs. (13) and (14) with the first Friedmann equation (5), we get


As usual the fractional energy densities are defined as




Using (12) we can rewrite equation (19) as


Based on these definitions, equation (15) can be rewritten as


where and we have defined


Next we take the time derivative of (15), after using (22), we find


Combining the above equation with Eqs. (8) and (12), we obtain the EoS parameter as


The first and second derivatives of the distance can be combined to obtain the acceleration parameter . It was shown that the zero redshift value of , is independent of space curvature, and can be obtained from the first and second derivatives of the coordinate distance Daly:2007dn (). It was argued that , which indicates whether the universe is accelerating at the current epoch, can be obtained directly from the supernova and radio galaxy data Daly:2007dn (). The acceleration parameter is given by


Using (24)the acceleration parameter (26) is obtained as


Finally, we obtain the equation of motions of GDE in BD theory. For this purpose, we first take the time derivative of relation (21). We find


Substituting from (27) into equation (28) and using relation , we get


In the remaining part of this paper we will constrain the GDEBD model by using the most recent observational date in the three different physical models: model I which is the GDEBD model in a flat universe ( and ), model II is the interacting GDEBD model in a flat universe ( and ) and finally model III is the interacting GDEBD model in a non-flat universe ( and ).

Iii Data fitting method

In this section we discuss the data fitting method in the Markov Chain Monte Carlo (MCMC) simulation to estimate the parameters of the model in section II using cosmological data.

To get the best fit values of the relevant parameters, the maximum likelihood method is used. The total likelihood function is defined as the product of the separate likelihood functions of uncorrelated observational data with


where SNIa stands for type Ia supernovae, CMB for cosmic microwave background radiation, BAO for baryon acoustic oscillation and gas stands for X-ray gas mass fraction data. Best fit values of parameters are obtained by minimizing . In this paper we use the cosmic microwave background radiation data from seven-year WMAP Komatsu:2010fb (), type Ia supernovae data from 557 Union2 Amanullah:2010vv (), baryon acoustic oscillation data from SDSS DR7 Percival:2009xn (), and the cluster X-ray gas mass fraction data from the Chandra X-ray observations Allen:2007ue (). In the rest of this section we discuss each in detail.

To obtain , we use seven-year WMAP data Komatsu:2010fb () with the CMB data point . The shift parameter R, which parametrize the changes in the amplitude of the acoustic peaks is given by Bond:1997wr ()


where is the redshift of recombination (see (36)), c is the speed of light in vacuum, is the present value of the matter density parameter (a ”0“ subscript shows the present value of the related quantity), and . In addition, the acoustic scale , which characterizes the changes of the peaks of CMB via the angular diameter distance out to recombination is defined as Bond:1997wr ()


The comoving distance is defined


and the comoving sound horizon distance at recombination is given by


in terms of the sound speed , defined by


The seven-year WMAP observations gives and Komatsu:2010fb ().

The redshift of recombination is obtained by using the fitting function proposed by Hu and Sugiyama Hu:1995en ()




Then one can define as , with Komatsu:2010fb ()


where is the inverse covariant matrix.

To obtain , the SNIa Union2 data Amanullah:2010vv () is used which includes 577 type Ia supernovae. The expansion history of the universe can be given by a specific cosmological model. To test this model, we can use the observational data for some predictable cosmological parameter such as luminosity distance . Assume that the Hubble parameter is used to describe the Universe, where parameters are predicted by a theoretical cosmological model. For such a theoretical model we can predict the theoretical ’Hubble-constant free’ luminosity distance as


where , is the redshift parameter, and

Then one can write the theoretical modulus distance


where . On the other hand, the observational modulus distance of SNIa, , at redshift is given by


where and are apparent and absolute magnitudes of SNIa respectively. Then the parameters of the theoretical model, s, can be determined by a likelihood analysis by defining in (30) as

where the nuisance parameter, , can be marginalized over as


Here we should mention an important point about using supernovae data to constrain the Brans-Dicke theories which have a varying gravitational coupling constant. Variations of gravitational coupling constant and apparent magnitude of supernovae are correlated as follows. The luminosity of a supernova is powered by Nickel-35 mass which is proportional to the Chandrasekhar mass


Moreover, the luminosity distance is the integral over the inverse Hubble parameter, which is proportional to . Therefore, the apparent magnitude


varies with a change in the gravitational coupling constant as


On the other hand, in the Brans-Dicke theory we have


where is the Brans-Dicke scalar field. In the slow roll approximation we can write


where is the decelaration parameter. The average value of between today and (the redshift when the SN measurement are probing the dark energy) is of order unity, and by using , we can write


By using Eqs. (46) and (49) we obtain


In Union2 data set, the redshift interval is between and , i.e. , with the systematic error of order in the measurement of apparent magnitude. Therefore, the systematic error can induce a bias roughly of order on parameter , which is three order of magnitudes larger than the statistical errors, as we will discuss in the next section. Therefore, in order to constrain with a higher precision, we combine the supernovae data with other cosmological data sets as follows. For more detailed discussion on possible evolution of the gravitational constant from cosmological type Ia supernovae see Gaztanaga:2001fh ().

The baryon acoustic oscillation data from the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) Percival:2009xn () is used here for constraining model parameters. The data constrain parameter , where is the comoving sound horizon distance (see (34)) at the drag epoch (where baryons were released from photons) and is given by Eisenstein:2005su ()


The drag redshift is given by the fitting formula Eisenstein:1997ik ()




Then we can obtain by , where


and its covariance matrix is given by Percival:2009xn ()


The ratio of X-ray gas mass to the total mass of a cluster is defined as the X-ray gas mass fraction Allen:2007ue (). The model fitted to the CDM model is Allen:2007ue ()


The elements in equation (60) are defined as follows: and are the proper angular diameter distance in CDM and the alternative theoretical model respectively, where


The angular correction factor


is caused by changes in angle for the alternative theoretical model compared to , where Allen:2007ue () is the slope of the data within the radius ( is the radius of the gas core in units).

The bias factor in equation (60) contains information about the uncertainties in the cluster depletion factor and the parameter accounts for departures from the hydrostatic equilibrium. The function denotes the uncertainties of the baryonic mass fraction in stars with a Gaussian prior for , with Allen:2007ue (). The factor describes the combined effects of the residual uncertainties, such as the instrumental calibration. A Gaussian prior for the ’calibration’ factor is considered as Allen:2007ue ().

Then is defined as Allen:2007ue ()

with the statistical uncertainties and


At the end of this section we should assert that the data points parameters of the CMB and BAO data sets which we use in this paper are the best fit values for CDM and the error estimates are also based on the CDM model. Therefore they are not completely accurate in this application. However they are the only parameters which we have to constrain our model.

Iv Results

Finally we apply a Markov chain Monte Carlo simulation on the parameters of the GDEBD model by using the publicly available CosmoMC code Lewis:2002ah (). The parameter vectors are , , for the flat non-interacting (model I), flat interacting (model II) and non-flat interacting (model III) models respectively. The basic cosmological parameters are taken in the following priors: , and . In addition, as we mentioned in section II, for the model fitting, we consider condition. The results are presented in table 1.

From table 1 one can see that the main cosmological parameters , , and in all three models are compatible with the results of the CDM model Bennett:2012zja (). In addition, in the presence of interaction between DM and DE, the parameter decreases and so the Brans-Dicke parameter increases. The best fit value of the parameter in all three models is compatible with the results of other cosmological constraining works (however see the discussion following Eq. (43)). For example in Chen:1999qh () the authors by using the CMB temperature and polarization anisotropy data, found . Wu and Chen in Wu:2009zb (); Wu:2013 () by using the five-year WMAP and SDSS data obtained . For other cosmological constraints on the BD theory see Tsujikawa:2008uc (); Nagata:2002tm (); Nagata:2003qn (); Acquaviva:2004ti (); Schimd:2004nq (); Acquaviva:2007mm (). This estimated value is also compatible with the results of the solar system tests of the scalar-tensor theories such as the Cassini experiment where it has been obtained Bertotti:2003rm (); Acquaviva:2007mm (). The positive best fit value of parameter describe a conversion of dark matter to dark energy although both in flat and non-flat universes, in 1- CL, an inverse conversion is possible as well. The interacting DE and DM models have been constrained by observational data by many authors with different parametrization of the interacting parameter Chimento:2003iea (); Zimdahl:2002zb (); Valiviita:2009nu (); Pettorino:2012ts (); Salvatelli:2013wra (); He:2010im (); Clemson:2011an (). He et al. in He:2010im () have parametrized the interaction parameter as in the present paper although they have chosen the prior on parameter as . They obtained the best fit value of parameter as .

Parameter Non-interacting Interacting Interacting non-flat
  …   …  
Table 1: The best fit values of the cosmological and model parameters in the the GDE model in the BD theory with and regions. Here CMB, SNIa and BAO and X-ray mass gas fraction data together with the BBN constraints have been used.
Figure 1: 2-dimensional constraint of the cosmological and model parameters contours in the flat non-interacting GDE model in the BD theory with and regions. To produce these plots, SNIa+CMB+BAO+X-ray gas mass fraction data together with the BBN constraints have been used.
Figure 2: 2-dimensional constraint of the cosmological and model parameters contours in the flat interacting GDE model in the BD theory with and regions. To produce these plots, SNIa+CMB+BAO+X-ray gas mass fraction data together with the BBN constraints have been used.
Figure 3: 2-dimensional constraint of the cosmological and model parameters contours in the non-flat interacting GDE model in the BD theory with and regions. To produce these plots, SNIa+CMB+BAO+X-ray gas mass fraction data together with the BBN constraints have been used.

V Conclusion

In this paper we considered the cosmological constraints on the parameters of the GDE in the framework of BD theory by using a Markov Chain Monte Carlo simulation. We used the SNIa+ CMB+ BAO+X-ray gas mass fraction data for the model fitting. The best fit values of the cosmological parameters in this model are compatible with the results of the CDM model. In addition, we obtained the best fit values of parameters and where is the BD parameter and is the interacting parameter. The best fit values of these parameters are also compatible with the results of previous constraining works. However as we mentioned in section III, due to large systematic error in by using the supernovae data only, to constrain this parameter with a higher precision, we should combine the supernovae data with other cosmological data sets as the CMB and BAO data sets. In addition we should assert one more time that the data points parameters of the CMB and BAO data which we have used in this paper are the best fit values for CDM and the error estimates are also based on the CDM model. Therefore they are not completely accurate in this application. The numerical results can be improved in the future works by using more recent data such as nine-yaer WMAP Bennett:2012zja () or the Planck Ade:2013zuv () projects.


We would like to thank the anonymous referee for constructive comments. Hamzeh Alavirad would also like to thank Qader Dorosti for helpful discussions. The work of A. Sheykhi has been supported financially by Research Institute for Astronomy and Astrophysics of Maragha (RIAAM), Iran.


  • (1) A. G. Riess et al., “Observational evidence from supernovae for an accelerating universe and a cosmological constant,” Astron.J., vol. 116, pp. 1009–1038, 1998.
  • (2) S. Perlmutter et al., “Measurements of and from 42 high redshift supernovae,” Astrophys. J., vol. 517, pp. 565–586, 1999.
  • (3) C. Wetterich, “Cosmology and the Fate of Dilatation Symmetry,” Nucl.Phys., vol. B302, p. 668, 1988.
  • (4) B. Ratra and P. Peebles, “Cosmological Consequences of a Rolling Homogeneous Scalar Field,” Phys.Rev., vol. D37, p. 3406, 1988.
  • (5) C. Armendariz-Picon, V. F. Mukhanov, and P. J. Steinhardt, “A Dynamical solution to the problem of a small cosmological constant and late time cosmic acceleration,” Phys.Rev.Lett., vol. 85, pp. 4438–4441, 2000.
  • (6) C. Armendariz-Picon, V. F. Mukhanov, and P. J. Steinhardt, “Essentials of k essence,” Phys.Rev., vol. D63, p. 103510, 2001.
  • (7) T. Padmanabhan, “Accelerated expansion of the universe driven by tachyonic matter,” Phys.Rev., vol. D66, p. 021301, 2002.
  • (8) A. Sen, “Tachyon matter,” JHEP, vol. 0207, p. 065, 2002.
  • (9) R. Caldwell, “A Phantom menace?,” Phys.Lett., vol. B545, pp. 23–29, 2002.
  • (10) S. Nojiri and S. D. Odintsov, “Quantum de Sitter cosmology and phantom matter,” Phys.Lett., vol. B562, pp. 147–152, 2003.
  • (11) S. Nojiri and S. D. Odintsov, “DeSitter brane universe induced by phantom and quantum effects,” Phys.Lett., vol. B565, pp. 1–9, 2003.
  • (12) N. Arkani-Hamed, H.-C. Cheng, M. A. Luty, and S. Mukohyama, “Ghost condensation and a consistent infrared modification of gravity,” JHEP, vol. 0405, p. 074, 2004.
  • (13) F. Piazza and S. Tsujikawa, “Dilatonic ghost condensate as dark energy,” JCAP, vol. 0407, p. 004, 2004.
  • (14) E. Elizalde, S. Nojiri, and S. D. Odintsov, “Late-time cosmology in (phantom) scalar-tensor theory: Dark energy and the cosmic speed-up,” Phys.Rev., vol. D70, p. 043539, 2004.
  • (15) S. Nojiri, S. D. Odintsov, and S. Tsujikawa, “Properties of singularities in (phantom) dark energy universe,” Phys.Rev., vol. D71, p. 063004, 2005.
  • (16) A. Anisimov, E. Babichev, and A. Vikman, “B-inflation,” JCAP, vol. 0506, p. 006, 2005.
  • (17) E. Witten, “The Cosmological constant from the viewpoint of string theory,” pp. 27–36, 2000.
  • (18) R.-G. Cai, “A Dark Energy Model Characterized by the Age of the Universe,” Phys.Lett., vol. B657, pp. 228–231, 2007.
  • (19) H. Wei and R.-G. Cai, “A New Model of Agegraphic Dark Energy,” Phys.Lett., vol. B660, pp. 113–117, 2008.
  • (20) E. J. Copeland, M. Sami, and S. Tsujikawa, “Dynamics of dark energy,” Int.J.Mod.Phys., vol. D15, pp. 1753–1936, 2006.
  • (21) M. Li, X.-D. Li, S. Wang, and Y. Wang, “Dark Energy,” Commun.Theor.Phys., vol. 56, pp. 525–604, 2011.
  • (22) A. De Felice and S. Tsujikawa, “ theories,” Living Reviews in Relativity, vol. 13, p. 3, 2010.
  • (23) G. Dvali, G. Gabadadze, and M. Porrati, “4-D gravity on a brane in 5-D Minkowski space,” Phys.Lett., vol. B485, pp. 208–214, 2000.
  • (24) M. S. Carena, J. Lykken, M. Park, and J. Santiago, “Self-accelerating Warped Braneworlds,” Phys.Rev., vol. D75, p. 026009, 2007.
  • (25) M. Minamitsuji, “Self-accelerating solutions in cascading DGP braneworld,” Phys.Lett., vol. B684, pp. 92–95, 2010.
  • (26) A. Sheykhi, B. Wang, and N. Riazi, “String inspired explanation for the super-acceleration of our universe,” Phys.Rev., vol. D75, p. 123513, 2007.
  • (27) K. Kawarabayashi and N. Ohta, “The Problem of in the Large Limit: Effective Lagrangian Approach,” Nucl.Phys., vol. B175, p. 477, 1980.
  • (28) E. Witten, “Current Algebra Theorems for the U(1) Goldstone Boson,” Nucl.Phys., vol. B156, p. 269, 1979.
  • (29) G. Veneziano, “U(1) Without Instantons,” Nucl.Phys., vol. B159, pp. 213–224, 1979.
  • (30) C. Rosenzweig, J. Schechter, and C. Trahern, “Is the Effective Lagrangian for QCD a Sigma Model?,” Phys.Rev., vol. D21, p. 3388, 1980.
  • (31) N. Ohta, “Dark Energy and QCD Ghost,” Phys.Lett., vol. B695, pp. 41–44, 2011.
  • (32) F. R. Urban and A. R. Zhitnitsky, “The cosmological constant from the QCD Veneziano ghost,” Phys.Lett., vol. B688, pp. 9–12, 2010.
  • (33) R.-G. Cai, Z.-L. Tuo, H.-B. Zhang, and Q. Su, “Notes on Ghost Dark Energy,” Phys.Rev., vol. D84, p. 123501, 2011.
  • (34) A. Sheykhi and M. Sadegh Movahed, “Interacting Ghost Dark Energy in Non-Flat Universe,” Gen.Rel.Grav., vol. 44, pp. 449–465, 2012.
  • (35) E. Ebrahimi and A. Sheykhi, “Instability of QCD ghost dark energy model,” Int.J.Mod.Phys., vol. D20, pp. 2369–2381, 2011.
  • (36) A. Sheykhi, M. Sadegh Movahed, and E. Ebrahimi, “Tachyon Reconstruction of Ghost Dark Energy,” Astrophys.Space Sci., vol. 339, pp. 93–99, 2012.
  • (37) A. Sheykhi and A. Bagheri, “Quintessence Ghost Dark Energy Model,” Europhys.Lett., vol. 95, p. 39001, 2011.
  • (38) A. Sheykhi, “Viscous ghost dark energy with a varying gravitational constant,” Phys.Scripta, vol. 85, p. 045901, 2012.
  • (39) A. Rozas-Fernandez, “Kinetic k-essence ghost dark energy model,” Phys.Lett., vol. B709, pp. 313–321, 2012.
  • (40) K. Karami, M. Mousivand, S. Asadzadeh, and Z. Safari, “QCD modified ghost scalar field dark energy models,” Int.J.Mod.Phys., vol. D22, p. 1350018, 2013.
  • (41) A. Khodam-Mohammadi and M. Malekjani, “Reconstruction of modified gravity with ghost dark energy models,” Mod.Phys.Lett., vol. A27, p. 1250100, 2012.
  • (42) M. Malekjani and A. Khodam-Mohammadi, “Statefinder diagnosis and the interacting ghost model of dark energy,” Astrophys.Space Sci., vol. 343, pp. 451–461, 2013.
  • (43) C.-J. Feng, X.-Z. Li, and X.-Y. Shen, “Latest Observational Constraints to the Ghost Dark Energy Model by Using Markov Chain Monte Carlo Approach,” Phys.Rev., vol. D87, p. 023006, 2013.
  • (44) R.-G. Cai, Z.-L. Tuo, Y.-B. Wu, and Y.-Y. Zhao, “More on QCD Ghost Dark Energy,” Phys.Rev., vol. D86, p. 023511, 2012.
  • (45) C. Brans and R. Dicke, “Mach’s principle and a relativistic theory of gravitation,” Phys.Rev., vol. 124, pp. 925–935, 1961.
  • (46) B. Bertotti, L. Iess, and P. Tortora, “A test of general relativity using radio links with the Cassini spacecraft,” Nature, vol. 425, p. 374, 2003.
  • (47) X.-l. Chen and M. Kamionkowski, “Cosmic microwave background temperature and polarization anisotropy in Brans-Dicke cosmology,” Phys.Rev., vol. D60, p. 104036, 1999.
  • (48) V. Acquaviva and L. Verde, “Observational signatures of Jordan-Brans-Dicke theories of gravity,” JCAP, vol. 0712, p. 001, 2007.
  • (49) S. Tsujikawa, K. Uddin, S. Mizuno, R. Tavakol, and J. Yokoyama, “Constraints on scalar-tensor models of dark energy from observational and local gravity tests,” Phys.Rev., vol. D77, p. 103009, 2008.
  • (50) F. Wu and X. Chen, “Cosmic microwave background with Brans-Dicke gravity II: constraints with the WMAP and SDSS data,” Phys.Rev., vol. D82, p. 083003, 2010.
  • (51) F. Wu and X. Chen, “Constraints on the Brans-Dicke gravity theory with the Planck data,” Phys.Rev., vol. D88, p.  084053 (2013).
  • (52) E. Ebrahimi and A. Sheykhi, “Interacting Ghost Dark Energy in Brans-Dicke Theory,” Phys.Lett., vol. B706, pp. 19–25, 2011.
  • (53) K. Saaidi, “Stability of Ghost Dark Energy in CBD Model of Gravity,” 2012.
  • (54) E. Komatsu et al., “Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation,” Astrophys.J.Suppl., vol. 192, p. 18, 2011.
  • (55) R. Amanullah, C. Lidman, D. Rubin, G. Aldering, P. Astier, et al., “Spectra and Light Curves of Six Type Ia Supernovae at and the Union2 Compilation,” Astrophys.J., vol. 716, pp. 712–738, 2010.
  • (56) W. J. Percival et al., “Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample,” Mon.Not.Roy.Astron.Soc., vol. 401, pp. 2148–2168, 2010.
  • (57) S. Allen, D. Rapetti, R. Schmidt, H. Ebeling, G. Morris, et al., “Improved constraints on dark energy from Chandra X-ray observations of the largest relaxed galaxy clusters,” Mon.Not.Roy.Astron.Soc., vol. 383, pp. 879–896, 2008.
  • (58) A. Lewis and S. Bridle, “Cosmological parameters from cmb and other data: A monte carlo approach,” Phys. Rev. D, vol. 66, p. 103511, Nov 2002.
  • (59) M. Arik and M. Calik, “Can Brans-Dicke scalar field account for dark energy and dark matter?,” Mod.Phys.Lett., vol. A21, pp. 1241–1248, 2006.
  • (60) E. Gaztanaga, E. Garcia-Berro, J. Isern, E. Bravo and I. Dominguez, Phys. Rev. D 65, 023506, 2002.
  • (61) C. Bennett et al., “Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results,” 2012.
  • (62) O. Bertolami, F. Gil Pedro, and M. Le Delliou, “Dark Energy-Dark Matter Interaction and the Violation of the Equivalence Principle from the Abell Cluster A586,” Phys.Lett., vol. B654, pp. 165–169, 2007.
  • (63) S. Tsujikawa and M. Sami, “A Unified approach to scaling solutions in a general cosmological background,” Phys.Lett., vol. B603, pp. 113–123, 2004.
  • (64) N. Banerjee and D. Pavon, “Holographic dark energy in Brans-Dicke theory,” Phys.Lett., vol. B647, pp. 477–481, 2007.
  • (65) A. Sheykhi, “Interacting holographic dark energy in Brans-Dicke theory,” Phys.Lett., vol. B681, pp. 205–209, 2009.
  • (66) C. M. Will, Theory and Experiment in Gravitational physics, 2nd ed. Basic Books/Perseus Group, New York, 1993.
  • (67) R. A. Daly, S. Djorgovski, K. A. Freeman, M. P. Mory, C. O’Dea, et al., “Improved Constraints on the Acceleration History of the Universe and the Properties of the Dark Energy,” Astrophys.J., vol. 677, pp. 1–11, 2008.
  • (68) J. Bond, G. Efstathiou, and M. Tegmark, “Forecasting cosmic parameter errors from microwave background anisotropy experiments,” Mon.Not.Roy.Astron.Soc., vol. 291, pp. L33–L41, 1997.
  • (69) W. Hu and N. Sugiyama, “Small scale cosmological perturbations: An Analytic approach,” Astrophys.J., vol. 471, pp. 542–570, 1996.
  • (70) D. J. Eisenstein et al., “Detection of the baryon acoustic peak in the large-scale correlation function of SDSS luminous red galaxies,” Astrophys.J., vol. 633, pp. 560–574, 2005.
  • (71) D. J. Eisenstein and W. Hu, “Baryonic features in the matter transfer function,” Astrophys.J., vol. 496, p. 605, 1998.
  • (72) A. Lewis and S. Bridle, “Cosmological parameters from CMB and other data: A Monte Carlo approach,” Phys.Rev., vol. D66, p. 103511, 2002.
  • (73) R. Nagata, T. Chiba, and N. Sugiyama, “Observational consequences of evolution of primordial fluctuations in scalar - tensor cosmology,” Phys.Rev., vol. D66, p. 103510, 2002.
  • (74) R. Nagata, T. Chiba, and N. Sugiyama, “WMAP constraints on scalar- tensor cosmology and the variation of the gravitational constant,” Phys.Rev., vol. D69, p. 083512, 2004.
  • (75) V. Acquaviva, C. Baccigalupi, S. M. Leach, A. R. Liddle, and F. Perrotta, “Structure formation constraints on the Jordan-Brans-Dicke theory,” Phys.Rev., vol. D71, p. 104025, 2005.
  • (76) C. Schimd, J.-P. Uzan, and A. Riazuelo, “Weak lensing in scalar-tensor theories of gravity,” Phys.Rev., vol. D71, p. 083512, 2005.
  • (77) L. P. Chimento, A. S. Jakubi, D. Pavon, and W. Zimdahl, “Interacting quintessence solution to the coincidence problem,” Phys.Rev., vol. D67, p. 083513, 2003.
  • (78) W. Zimdahl and D. Pavon, “Scaling cosmology,” Gen.Rel.Grav., vol. 35, pp. 413–422, 2003.
  • (79) J. Valiviita, R. Maartens, and E. Majerotto, “Observational constraints on an interacting dark energy model,” Mon.Not.Roy.Astron.Soc., vol. 402, pp. 2355–2368, 2010.
  • (80) V. Pettorino, L. Amendola, C. Baccigalupi, and C. Quercellini, “Constraints on coupled dark energy using CMB data from WMAP and SPT,” Phys.Rev., vol. D86, p. 103507, 2012.
  • (81) V. Salvatelli, A. Marchini, L. Lopez-Honorez, and O. Mena, “New constraints on Coupled Dark Energy from the Planck satellite experiment,” Phys.Rev., vol. D88, p. 023531, 2013.
  • (82) J.-H. He, B. Wang, and E. Abdalla, “Testing the interaction between dark energy and dark matter via latest observations,” Phys.Rev., vol. D83, p. 063515, 2011.
  • (83) T. Clemson, K. Koyama, G.-B. Zhao, R. Maartens, and J. Valiviita, “Interacting Dark Energy – constraints and degeneracies,” Phys.Rev., vol. D85, p. 043007, 2012.
  • (84) P. Ade et al., “Planck 2013 results. XVI. Cosmological parameters,” preprint, astro-ph.CO/1303.5076, 2013.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description