Inflationary generalized Chaplygin gas and general dark energy in the light of the Planck and BICEP2 experiments
Abstract
In this work, we study an inflationary scenario in the presence of Generalized Chaplygin Gas (GCG). We show that in Einstein gravity, GCG is not a suitable candidate for inflation; but in a five dimensional brane world scenario, it can work as a viable inflationary model. We calculate the relevant quantities such as , and related to the primordial scalar and tensor fluctuations, and using their recent bounds from Planck and BICEP2, we constrain the model parameters as well as the fivedimensional Planck mass. But as a slowroll inflationary model with a powerlaw type scalar primordial power spectrum, GCG as an inflationary model can not resolve the tension between results from BICEP2 and Planck with a concordance CDM Universe. We show that going beyond the concordance CDM model and incorporating more general dark energy behaviour, this tension may be eased. We also obtain the constraints on the and and the GCG model parameters using Planck+WP+BICEP2 data considering the CPL dark energy behaviour.
pacs:
98.80.k,98.80.CqMarch 21, 2018
I Introduction
The present cosmological observations are amazingly consistent with a Universe which has two accelerated expansion phases in its entire evolution history from big bang till today. One of these accelerating epochs presumably occurred during the early phase of the cosmological evolution when the energy scale of the Universe was close to the Planck scale. This accelerating period was first proposed around 1980 in order to solve the puzzles like flatness, horizon and monopole problems in standard cosmology (1). An epoch of exponential expansion was proposed in order to solve these problems. Although a simple cosmological constant can give rise to such exponential expansion, the Universe can never exit from this accelerating phase in such a scenario and hence does not enter into a decelerated phase which is necessary for subsequent processes like nucleosynthesis, structure formation. To solve this exit problem, a scalar field theory was introduced where the field rolls over a sufficiently flat potential (slow rolling) and can mimic a cosmological constant like behaviour. Such a scalar field (inflaton) can drive a near exponential expansion. The exit from the inflationary era is ensured as the scalar field reaches the nonflat region (fast roll phase) of the potential. One of the greatest successes of inflationary model is the generation of primordial density fluctuations in the Universe which can act as the seeds for the large scale inhomogeneities that is necessary for the structure formation of the Universe (2). The quantum fluctuations of the inflaton during inflation can produce such primordial density fluctuations. Given any scalar field inflationary model, one can calculate the spectrum of this primordial density fluctuation which is essentially related to the scalar part of the metric fluctuations. Moreover, one can also calculate the tensor fluctuations in the metric produced during inflation which result a stochastic gravitational wave background on large cosmological scales. Cosmologically both these primordial fluctuations are interesting as they produce observable features in the temperature anisotropy of the cosmic microwave background radiation (CMBR). Hence measuring temperature anisotropy in CMBR enables us to constrain the spectrum of these initial fluctuations which in turn can constrain the gravitational physics close to Planck scale. The scalar part of the temperature anisotropy was first measured by COBE (3) in early nineties and subsequently by a host of CMBR experiments e.g BOOMERANG (4), WMAP (5) and more recently by Planck (6). For the Polarization in CMBR, DASI (7) first detected the Emode polarization in CMBR in 2001. But the Bmode polarization in CMBR which is a clear evidence for the existence of primordial gravitational waves generated through the tensor fluctuations during inflation had not been detected until recently. But just recently, the BICEP2 experiment (8) has announced the detection of the Bmode polarization signal in CMBR ruling out the zero tensor fluctuation at confidence level. This is an extraordinary result for cosmology and if confirmed by future polarization data from Planck satellite, it will establish the fact that there was an accelerating epoch in the Universe prior to the radiation era. From the measured value of ( tensor to scalar ratio), one can also estimate the energy scale for this accelerating regime to be around GUT scale () which is below the Planck scale (the scale where the quantum gravity effects are prominent) but higher than the TeV scale ( the scale which can be probed by the current particle accelerator like LHC) (see (9) for an interesting discussion on this issue).
For inflationary model building using scalar fields, one needs to guarantee that there is sufficient slowrolling for the scalar field to ensure necessary amount of inflation to solve the horizon and flatness problem. One also has to ensure the eventual breakdown of the slowroll conditions so that the Universe exits from this inflationary phase and enters into a decelerated regime. This restricts the shape of the potential for the inflaton. Moreover most of the slowroll scalar field models produce powerlaw type primordial power spectrum (PPS) of the form where and are related to the the shape of the underlining potential for the inflaton. Measuring the temperature anisotropies in CMBR, one can put strong bounds on and and that puts further constraints on the inflaton potential (See (10) for scalar field models for inflation that are consistent with recent Planck results).
Recently Generalized Chaplygin Gas (GCG) (11) described by an equation of state where and are constants, has been discussed widely in cosmological contexts. The case is interesting in the context of late time acceleration. In this case, the GCG mimics dust in the early time and a dark energy with negative equation of state in late time. Initially considered as a natural candidate for a unification of dark matter and dark energy (UDM), it was later shown that this particular UDM behaviour is not suitable for the structure formation of the Universe (12). But this fluid can be a possible dark energy candidate which tracks the background fluid initially and subsequently exits the tracking regime and starts acting as a dark energy candidate (13).
The opposite regime, is also equally interesting (14). In this case, the GCG behaves like a cosmological constant () initially but with time, the equation of state increases and becomes dust like (). This behaviour is suitable for inflation, as in this case inflation happens initially and then automatically ends at later times. We shall show that by properly adjusting parameters, one can get enough inflation that is necessary to solve horizon and flatness problems.
Motivated by this, we consider GCG as a model for inflation. We write the corresponding scalar field theory that mimics such behaviour and then calculate the PPS in this model. We show that in the context of Einstein gravity, to get the right shape for the PPS, the efolding at the time of horizon exit () has to be excessively large which is a serious drawback. But if one considers a five dimensional braneworld scenario which results in a correction term in the Einstein equation, this problem of high efolding at horizon exit gets resolved and one gets a suitable inflationary model.
Recently, it has been pointed out that a power law form for the scalar PPS for the inflaton field is in tension with combined Planck+BICEP2 results (15). This is related to the fact that a significantly higher value for () as measured by BICEP2 is not consistent with the suppression of power in at large scales as observed by Planck. In fact the authors in (15) have shown that power law form for the scalar PPS with a single spectral index is ruled out at more than by Planck+BICEP2 in comparison to a broken PPS model containing two spectral indices (see (16) for different approaches to solve this problem). This is a bad news for inflationary model building because most of the standard and theoretically motivated slowroll inflationary scenarios produce a powerlaw type scalar PPS. GCG as a inflationary model is of slowroll type and also produces a powerlaw type scalar PPS. Hence this tension applies to GCG as well. But the underlining assumption for all these studies is that our Universe is described by a concordance CDM model. We try to address this issue by going beyond the concordance CDM model. By allowing a general dark energy equation of state, our study shows that one may address this issue even if one sticks to a power law type scalar PPS.
The structure of the paper is as follows: in section 2, we describe the GCG inflationary models and its scalar field representation; in section 3, we study the slowroll inflationary models with GCG in Einstein gravity and discuss its problem; in section 4, we study the GCG inflation in a particular higher dimensional braneworld set up and put constraints on various model parameters using observational results from Planck and BICEP2 for a CDM; in section 5, we discuss the issue regarding the inconsistency between Planck and BICEP2 results with powerlaw type scalar PPS and try to address the issue with a general dark energy model. We do the full MCMC analysis with Planck+WP and BICECP2 data using a general dark energy model to get the constraint on our inflationary GCG model parameters and compare the results obtained using a CDM model for dark energy; finally in section 6, we put our concluding remarks.
Ii The GCG Inflation
The Generalized Chaplygin Gas (GCG) is described by the equation of state (11)
(1) 
where is the pressure and is the energy density of the GCG fluid. and are the parameters of the model. One can calculate energy density as a function of the scale factor by integrating out the energy conservation equation in FRW background and the corresponding expression is
(2) 
where, and is an integration constant. For , GCG behaves like a dust at very early time and it behaves like cosmological constant at infinite future. In between the equation of state smoothly changes over from dust behaviour to the cosmological constant. This behaviour is attractive for dark energy model building as the GCG tracks the background matter in the early time and then enters the dark energy regime in the late time. This is similar to tracking model of dark energy that attempts to solve the cosmic coincidence problem. On the other hand, this behaviour is not at all suitable for early time inflation, because in this case once GCG starts accelerating the Universe, it can never be stopped unless one invokes some extra effect to exit from the inflationary period (see (17) for inflationary model with GCG with ).
For , the opposite happens. In this case the GCG behaves like a cosmological constant (CC) to start with, and then it slowly evolves away from this CC behaviour and eventually behaves like a dust. In this case, we have a inflationary epoch for early time which ends subsequently and the Universe enters into a decelerating dust like era. The end of inflation is automatic in this case without any need for extra mechanism. One can suitably choose the model parameters to get the required number of efolds. With this, it is now important to see whether the primordial fluctuations that can be produced in such a model, is consistent with the observational results from experiments like Planck and BICEP2.
GCG can be described by a minimally coupled scalar field Lagrangian with a canonical kinetic energy term. The energy density and pressure for a canonical scalar field which is minimally coupled to the gravity are given by (assuming flat FRW metric):
(3) 
where is the potential for the field. Using these expressions and equations (1) and (2), together with the Einstein’s equation
(4) 
one can write
(5) 
where . On integration, this results
(6) 
where, and is an integration constant. After some straightforward algebraic calculations, one can finally write
(7) 
where and .
Iii SlowRoll Inflationary Model with GCG
Using the expression for the Hubble parameter given in the previous section, one can now calculate the two Hubble slowroll parameters and (one can also use the potential slowroll parameters and and the results will be exactly the same). They are expressed as:
(8)  
(9)  
Under the slow roll approximation, the scalar spectral index and tensor spectral index are given by
(10) 
Using the above expressions for and , one can write their expressions as (2):,
(11) 
The tensor to scalar ratio which measures the amount of stochastic gravitational wave that is produced during inflation, is given by
(12) 
One should note that all the relevant quantities like , , should be calculated at , i.e when a given mode of fluctuation exits the horizon. Once the fluctuations exit the horizon, they do not evolve and are frozen at their values at horizon crossing. As the inflation ends, the horizon scale starts growing, and different fluctuations start entering inside the horizon. The fluctuations which are larger in scales enters later and fluctuations which are smaller in scales enter earlier. The fluctuations which are of horizon size today were the last to exit the horizon during inflation. Any scale today that we are interested in, can be related to the number of efolds before the end of inflation which is given by (2)
(13) 
where subscript ”e” denotes the end of inflation and subscript denotes differentiation w.r.t . For all the relevant scales that can be probed through CMB observations like Planck, (18).
Inflation will end when . Using the expression (8), this gives
(14) 
The scale factor at the end of inflation can also be calculated from equation (6) as
(15) 
Knowing that , one can now calculate the value of the scalar field at horizon exit, using equations (6), (14) and (15):
(16) 
Using this expression for , , we get the value of the scalar spectral index and tensortoscalar ratio at horizon exit as:
(17) 
(18) 
We should mention that these bounds are obtained using a CDM dark energy model. In section 5, we obtain the similar bound assuming a more general dark energy model.
In Figure (1), we draw the contours in the plane which satisfy the above constraints. One can see that minimum value for required is which is way above the theoretical prior . This is the main drawback of this model.
So, the generalized chaplygin gas is not a suitable model for inflation in Einstein’s gravity.
Iv Inflationary Model in BraneWorld Scenario with GCG
In this scenario, we consider the observable Universe to be confined on a 3brane embedded in a 5D anti deSitter spacetime. One such scenario was first proposed in 1998 by Randall and Sundrum (RS) (19) to solve the hierarchy problem in particle physics. This scenario consists of a 5D spacetime governed by Einstein gravity with a negative cosmological constant in the bulk. The spacetime respects the symmetry and the flat 3branes are located at orbifold fixed points in this geometry. One of the branes is our visible Universe where the modified Einstein equation is now given by (20):
(19) 
where, is the 3brane tension. The relation between and fivedimensional Planck mass is given by
(20) 
(21) 
In the high energy limit , one can approximate the slowroll parameters as
(22) 
(23) 
where and . The scalar spectral index can be written in terms of the potential slowroll parameters as:
(24) 
The amplitude of the scalar as well as the tensor perturbations are defined as (22)
(25) 
(26) 
where,
(27) 
(28) 
Using Equations (25) and (26), we can now write the tensortoscalar ratio as,
(29) 
In the high energy limit , and hence in this limit, becomes
(30) 
Using the above definitions, we get the expressions for the slowroll parameters for our model as
(31)  
(32)  
We also calculate the value of the field at horizon exist as
(33) 
where is the number of efold at the horizon exit. Using this, we calculate different observables like , at the horizon exit as
(34)  
(35) 
The current bound on these two parameters from Planck+WP+BICEP2 with a CDM model are given in equation (18). Next, we calculate the initial value of the field which gives total 70 efolds of inflation, necessary to solve the horizon problem:
(36) 
where, is total number of efolds during inflation. From equation (25) and using the high energy limit, we get the expression for at the horizon exit as
(37)  
The measured value for by Planck for CDM model is given as (6). Subsequently we study the parameter space () that is allowed by the results obtained by Planck+WP and BICEP2 for the observables , and . This is shown in figure (2) for different values of . We should again stress that these results are obtained assuming a CDM Universe.
The first thing to be noticed is that we now have a parameter space that is allowed by the Planck+WP+BICEP2 data with . This is due to the modified Einstein’s equation in RS braneworld set up.
In these figures, the shaded regions represent the allowed parameter space that satisfies the constraints on and as obtained by Planck and BICEP2 respectively. Once we have such a allowed region, one can then fix the energy scale of the inflation as given by which satisfies the constraint on . The dotted line in each figure represents a typical behaviour in () parameter space for a particular choice of . The value of is chosen in such a way so that this line falls within the shaded region. This shows that typical scale of inflation is around GeV, just below the GUT scale. We should emphasis that the constraint on is very sensitive on this scale as slight deviation from this scale can move the dotted line outside the shaded region making it inconsistent with the allowed region for and . As is related to the brane tension which in turn is related to the five dimensional Planck mass , one can also estimate the 5D Planck mass . These are , and for respectively. From eqn (36), one can also calculate necessary to have efolds of total inflation. The required value for typically varies from to for different .
V Constraints with CPL dark energy
As discussed in the Introduction, it is now known that a measured high value for by BICEP2 is in tension with suppression of power at large scale as observed by Planck. In fact the recent paper by Hazra et al. (15) has ruled out the simple powerlaw form for scalar PPS in comparison to a broken scalar PPS using Planck+BICEP2 data at more than . And this is precisely due to the fact that with a powerlaw form for the scalar PPS, a high value of () as measured by BICEP2 is inconsistent with the suppression of power at large angular scales as observed by Planck. But this suppression can be achieved with a broken scalar PPS as shown in (15). Their investigation assumes a concordance CDM model for our Universe.
Comparison of the CDM with CPLDE  

Planck + WP  Planck + WP + BICEP2  

CDM  CPLDE  CDM  CPLDE 

0.02217  0.0223  0.0221  0.0223 

0.1183  0.1171  0.1177  0.116 

1.041  1.042  1.041  1.041 

0.088  0.088  0.089  0.089 

0.9658  0.9676  0.9686  0.9732 

1  1.408  1  1.599 

0  0.894  0  1.17 

0.009  0.01  0.16  0.17 

3.085  3.081  3.085  3.081 
[Best fit] 

commander  7.454  8.61  1.695  4.802 
CAMspec  7796.235  7795.474  7797.54  7796.988 
WP  2014.141  2014.55  2013.321  2013.572 
BICEP2      39.141  38.281 
Total  9802.92  9801.41  9848.31  9844.04 
  1.51    4.3  

Here we keep the powerlaw form for the scalar PPS but we deviate from the concordance CDM model by introducing a dynamical dark energy model given by the equation of state parameter as prescribed by Chevallier and Polarski (25) and Linder (26):
(38) 
In figure (3), we show the and variations with angular scale using the publicly available code CAMB (27). The consistency relation is assumed. We use the parameter initialization values as provided by BICEP2 (24). We also fix for our purpose. The topmost yellow line represents the behaviour with a power law PPS together with a concordance CDM model. It is easy to see the large enhancement of power at large scales, specially at , makes it inconsistent with the Planck measurements for . With a concordance CDM model and a powerlaw type scalar PPS, it is hard to escape from this inconsistency.
To ease this tension, we extend the parameter space by allowing a phantom equation of state () together with ( so that the equation of state remains phantom at all times) for the dark energy. To demonstrate how this may address the issue, we take two possible phantom dark energy models (DE1 and DE2) with two specific choices for the parameters and , e.g., and respectively. The rest of the parameters are fixed as in the CDM case described above. With these choices, one can now suppress the power at large scales to make it more consistent with Planck measurement as one can see from the plots in Figure (3). Moreover, the plot shows that if one assumes with a powerlaw type scalar PPS, both of these phantom models are consistent with Bmode polarization measurement by BICEP2.
Next, we also calculate the best fit values for CDM and a CPL dark energy model (CPLDE) using the Planck+WP and BICEP2 likelihoods. We use only the minimization routine in the publicly available code COSMOMC (28) for this. The results are shown in Table I.
Without the BICEP2 data, for CPLDE models, there is an improvement around in compared to the CDM model. This is consistent with the earlier results obtained by Hazra et al. (13). But with the addition of BICEP2 data, the improvement in is roughly for CPLDE models which is roughly three times of what one gets without BICEP2. We should stress that these numbers are indicative. It shows that with a full likelihood analysis for CMB+nonCMB data using MCMC, one can expect a substantially better fitting with a general dark energy keeping the powerlaw form for scalar PPS. In a recent paper, Hazra et al. (13) have shown that there is a mild preference for phantom model over CDM with the current CMB +nonCMB data (pre BICEP2) and CDM is disfavoured at more than confidence level but it is still allowed at . All these happen with simple power law type scalar PPS. Our result shows that with the inclusion of BICEP2 data, this result can change substantially.
bf Encouraged by this, we run the full MCMC chain with COSMOMC for a general CPLDE model taking the Planck+WP+BICEP2 data. As before, we use the parameter initialization values as provided by BICEP2 (24). Additionally, we use the gaussian prior for and with central values as and the standard deviations respectively .
The likelihoods for and are shown in figure 4. The error bars for and are
(39) 
One can see that going beyond concordance CDM model, the allowed value of has come down appreciably and at confidence limit, it is almost same as what Planck+WP obtained for , i.e .
Next, using the covariance between and as obtained from the full COSMOMC chains, we obtain the corresponding likelihood functions for GCG model parameters and and subsequently the likelihood contours in the parameter plane. These are shown in figures 5 and 6. One can see from figure 6 that the allowed regions in parameter space is decreased slightly compared to those obtained using a CDM model as shown in figure 2.
Finally, in figure 7, we draw the contours in the parameter space using Planck+WP+BICEP2 with CPLDE model and show different combinations of and which are allowed in this parameter space. Here we also show the same contours using only the Planck+WP data. It is evident that going beyond the concordance CDM model, these two data sets are consistent.
Vi Conclusion
In this paper, we have considered GCG as a natural candidate for the inflation. In this model, the equation of state for GCG starts with behaviour and leads to inflation. With time, the equation of state naturally evolves towards dust behaviour and inflation ends. Subsequently we study a canonical scalar field theory that represents the GCG behaviour.
While studying the primordial fluctuations in this model, we show that in Einstein gravity, GCG is not suitable for generating the required PPS as one needs a fairly large value for (), the required efolding at horizon exit, which is incompatible with theoretical constraint . Next we consider the GCG inflationary model in RS type five dimensional brane world scenario where the Einstein equation gets a correction term due to the presence of higher dimension. In this set up, we show that GCG works perfectly as a slowroll inflationary model. We obtain the constraints on the model parameters using the bounds on and as obtained by Planck and BICEP2. The inflationary energy scale in our model is around GeV, one order less than the GUT scale. The value of the fivedimensional Planck mass is around in our model.
As any other slowroll inflationary model with powerlaw type scalar PPS, GCG model is also in tension with combined Planck+BICEP2 data primarily due to the fact that large contribution from gravitational wave as measured by BICEP2, can not explain the suppression of power at large scales as observed by Planck. We show that by allowing a general dark energy equation of state given by CPL parametrization, one can ease this tension as a general dark energy behaviour may allow suppression of power at low even with . By calculating the best fit likelihood values for CDM and CPLDE model, we show that CPLDE model with a power law type scalar PPS is a better fit to the joint Planck+WP+BICEP2 data compared to a CDM with similar PPS. Hence allowing a deviation from concordance CDM model may save the simple slowroll inflationary models with powerlaw type scalar PPS. Finally we do a full MCMC analysis using COSMOMC with a CPLDE model using Planck+WP+BICEP2 data and obtain the revised estimate for and as well as our model parameters and .
It will be interesting to extend our work to phantom scalar inflationary models as previously studied in (29) and we hope to address this issue in future.
Acknowledgements
The authors would like to thank Dhiraj Kumar Hazra for extensive discussions and valuable comments. We also acknowledge the use of publicly available codes, CAMB and COSMOMC. A.A.S. acknowledges the funding from SERC, Dept. of Science and Technology, Govt. of India through the research project SR/S2/HEP43/2009. B.R.D. thanks CSIR, Govt. of India for financial support through JRF scheme (No:09/466(0157)/2012EMRI). S.K. thanks the UGC, Govt. of India for financial support through SRF scheme.
References

A. A .Starobinsky, JETP lett. 30, 682 (1979);
V .F. Mukhanov, G. V. Chibisov, Soviet Phys. JETP lett. 33 532 (1981);
A. H. Guth, Phys. Rev. D, 23, 347 (1980);
A. A. Linde, Phys. Lett. B, 108, 389 (1982);
A. Albrecht and P. Steinhardt, Phys. Rev. Lett., 48, 1220 (1982). 
V. F. Mukhanov, H. A. Feldman and R. H. Bradenberger, Phys. Rept., 215, 203 (1992);
A. Riotto, arXiv:hepph/0210162;
D. Baumann, arXiv:0907.5424 [hepth];
J. Lidsey et al., Rev. of Mod, Phys., 69, 373 (1997).  J. C. Mather et al. Astrophys. J., 420, 439 (1994).
 J. Ruhl et al. Astrophys. J, 599, 786 (2003).
 G. F. Hinshaw et al., Astrophys. J. Supplement, 19, 20B (2013).
 P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5075 [astroph.CO].
 J. Kovac et al., Nature, 420, 772 (2002).
 P. A. R. Ade et al. [BICEP2 Collaboration], arXiv:1403.3985 [astro.ph.CO].
 C. M. Ho and S. D. H. Hsu, arXiv:1404.0745 [hepph].
 J. Martin, C. Ringeval, R. Trotta and V. Vennin, arXiv:1312.3592 [astroph.CO].

A. Yu. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B, 511, 265 (2001);
N. Bilic, G. B. Tupper and R. D. Viollier, Phys. Lett. B, 535, 17 (2002);
M. C. Bento, O. Bertolami, A. A. Sen, Phys. Rev. D, 66, 043507 (2002).  H. Sandvik, M. Tegmark, M. Zaldariagga and I. Waga, Phys. Rev. D, 69, 123524 (2004).
 D. Hazra, S. Majumdar, S. Pal, S. Panda and A. A. Sen, arXiv:1310.6161 astro.ph.CO].
 A. A. Sen and R. J. Scherrer, Phys. Rev. D, 72, 063511 (2005).
 D. Hazra, A. Shafieloo, G. F. Smoot and A. A. Starobinsky, arXiv:1403.7786 [astroph.CO].

V. C. Miranda, W. Hu and P. Adshead, arXiv:1403.5231 [astroph.CO];
K. N. Abazajian, G. AAslanyan, R. Easther and L. C. Price, arXiv:1403.5922 [astroph.CO];
A. Ashoorioon, K. Dimopoulos, M. M. SheikhJabbari and G. Shiu, arXiv:1403.6099 [astroph.CO];
C. Bonvin, R. Durrer and R. Maartens, arXiv:1403.6768 [astroph.CO];
R. Bousso, D. Harlow and L. Senatore, arXiv:1404.2278 [astroph.CO];
D. K. Hazra, A. Shafieloo, G. F. Smoot and A. A. Starobinsky, arXiv:1404.0360 [astroph.CO]. 
O. Bertoami and V. Duvvuri, Phys. Lett .B, 640, 121 (2006);
R. Herrera, Phys. Lett. B, 664, 149 (2008);
S. del Campo, JCAP, 1311, 004 (2013).  P. A. R. Ade et al. [Planck Collaboration], arXiv:1303.5082 [astroph.CO].

L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999);
ibid., Phys. Rev. Lett. 83, 4690 (1999). 
P. Bintruy, C. Deffayet, U. Ellwanger, D. Langlois, Phys. Lett. B 477, 285 (2000);
T. Shiromizu, K. Maeda and M. Sasaki, Phys. Rev. D, 62 024012 (2000).  R. Maartens, D. Wands, B. A. Bassett, I. P. C. Heard, Phys. Rev. D 62 041301 (2000).
 D. Langlois, R. Maartens, D. Wands, Phys. Lett. B 489, 259 (2000).
 M. C. Bento, N. M. C. Santos, and A. A. Sen, Phys. Rev. D 69, 023508 (2004).
 See http://bicepkeck.org/
 M. Chevallier and D. Polarski, Int. J. Mod. Phys. D 10, 213 (2001)
 E. V. Linder, Phys. Rev. Lett. 90, 091301 (2003)
 A. Lewis, A. Challinor and A. Lasenby, Astrophys. J. 538 473, (2000); See, http://camb.info/.
 A. Lewis and S. Bridle, Phys. Rev. D 66 (2002) 103511 [astroph/0205436]; See, http://cosmologist.info/cosmomc/.
 C. J. Feng, X. Z. Li and E. N. Saridakis, Phys. Rev. D 82, 023526 (2010) [arXiv:1004.1874 [astroph.CO]].