Using the mathematical definitions of deceleration and jerk parameters we obtain a general functional form for Hubble parameter. By the aid of this exact Hubble function we can exactly reconstruct any other cosmographic parameters. We also obtained a general function for transition redshift as well as spacetime curvature. We may highlight the role of jerk parameter as any other cosmgraphic parameter could be written as a function of this parameter. Our derived functions clearly impose a lower limit on the jerk parameter which is . Moreover, we found that the jerk parameter indicates the geometry of the spacetime i.e any deviation from imply to a non-flat spacetime. In other word reefers to a dynamical, time varying, dark energy. From obtained Hubble function we recover the analogue of CDM model. To constrain cosmographic parameters as well as transition redshift and spacetime curvature of the recovered CDM model, we used Metropolis-Hasting algorithm to perform Monte Carlo Markov Chain analysis by using observational Hubble data obtained from cosmic chronometric (CC) technique, BAO/CMB data, Pantheon compilation of Supernovae type Ia, and their joint combination. The only free parameters are , and . From joint analysis we obtained , , , and .
Recovering CDM Model From a Cosmographic Study
Hassan Amirhashchi and Soroush Amirhashchi
Department of Physics, Mahshahr Branch, Islamic Azad University, Mahshahr, Iran
Keywords: Cosmography; Hubble Rate; Deceleration Parameter; Jerk Parameter
PACS Nos: 98.80.Es, 98.80.-k, 95.36.+x, 98.80.Jk
Since 1998 it has been revealed that the cosmic expansion is speeding up [1, 2, 3, 4]. In the context of general theory of relativity (GR) the existence of an exotic fluid called ‟dark energy (DE)˝with negative pressure is considered as the source caused present universe accelerated expansion [5, 6, 7]. It is also possible to study this cosmic acceleration in the context of modified gravity which is a generalization of general relativity [8, 9]. Moreover one may deal with this problem by considering violation of cosmological principle i.e assuming that the spacetime metric is inhomogeneous. Despite of about two decades effort, we still could not propose a realistic cosmological model in order to describe the present day cosmic accelerated expansion precisely. For example, although CDM model  excellently fits all observational data (for example see ) but considering cosmological constant as dark energy encounters to coincidence and fine-tuning problems [12, 13]. We know that the modern cosmology is based on the Friedmann equations, nevertheless, it is interesting to study universe through a kinematic approach rather than dynamical one (Einstein equations). This model-independent approach is called ‟Cosmography˝or cosmo-kinetics . It is worth noting that in this purely kinematic approach all the derived quantities are also model-independent. Cosmography was firstly introduced by Weinberg  in 2008 and later extended by Visser . The cornerstone of cosmography is to expand some observables such as the Hubble parameter(or equivalently the scale factor) into power series, and directly relating cosmological parameters to these observable quantities.
However, in practice, the cosmography study confronts to two serious problems with observational data. First of all, as shown in Ref , the Taylor expansion fails to reach convergence at redshift . This, of course an important shortcoming as many observational probes such as supernovae type Ia (SNIa) and cosmic microwave background (CMB) compilations can span the redshift region up to and respectively. This problem could be overcome by definition of an improved redshift parametrization such as [17, 18, 19]. The second one is the fact that in cosmography we use a finite Taylor series truncation which represent an approximation of the exact function and hence leads to worse estimations. Note that, taking more terms of Taylor series gives raise to more precise approximation but higher errors. Therefore as mentioned in Ref  ‟cosmography is in the dilemma between accuracy and precision˝. For recent cosmographty studies in the context for GR and Modified GR Reader is advised to see Refs [19, 20, 21, 22, 23, 24, 25, 26, 27, 28] an [29, 30] respectively. Ù°Very recently it has been shown that using Weighted Function Regression method  improves the usual cosmographic approach by automatically implementing Occamâs razor criterion .
As the source of all above mentioned problems lies in the Taylor expansion of the Hubble (scale factor) or luminosity distance, in this paper, in contrast with usual cosmography, we do not expand any of these parameters and instead we try to find an exact function for Hubble parameter on the bases of kinematic parameters of the universe. To do so, we mix the definition of deceleration and jerk parameters which in turn gives a second order differential equation for squared Hubble parameter. With out any prior assumption, we obtain an exact Hubble parameter as a function of jerk parameter (see sec 3). Through this Hubble parameter we reconstruct all other cosmographic parameters (CS) exactly. Then we use observational Hubble data (OHD) in the redshift range , Pantheon compilation containing 1048 SNIa apparent magnitude measurements over the redshift range of , BAO/CMB dataset , and their joint combination data to constrain CS parameters. We compare some of our results by those obtained in Refs [18, 20] and . This paper is organized as follows. In Sec 2 we briefly discuss cosmography and intruduce CS parameters up to the fifth derivative of the scale factor. In Sec 3 we derive a general deferential equation for squared Hubble parameter an reconstruct all CS parameter from it. Subsec 3.1 deals with the derivation (recovering) CDM model from our almost-general Hubble solution and show that how the spacetime geometry is connected to the jerk parameter. We summary the computational technique we have used to fit CS parameters to data by a numerical MCMC analysis in Sec 4. Sec 5 deals with the results of our fits to data. In Subsec 5.1 we derive a general redshift and constrain it over data. Finally, we summarize our findings and conclusions in Sec 6.
In this section we shall briefly describe cosmography which may start from Taylor series of scale factor.
Taylor expansion of the scale factor around the current time gives
Without loss of generality we can assume , where the constant is the current value of the scale factor. The most important cosmographic series terms i.e the Hubble, deceleration, jerk, snap, and lerk parameters are [16, 32, 36]
respectively. As noted in Ref  the first three CS parameters i.e the Hubble rate , deceleration parameter and its first derivative with respect to the cosmic time ( or redshift) are sufficient to determine the overall kinematics of the Universe. However, at the current time the deceleration parameter is restricted as which in turn impose . For CDM model at any time .
Using eqs (3) in the Taylor expansion of the Hubble parameter around one can obtain
In other hand, the cosmographic version of the luminosity distance can be conveniently expressed as
where the subscript “0”indicates the present values of the cosmographic parameters.
Here we again emphasize that there are two main problems arising in the context of
cosmography when using eqs (2) and (2) to constrain CS parameters. In fact, it has already been shown that expanding around gives raise to the divergence of Taylor series at . This problem could be limited by transferring to ). The second problem is the truncation of the series one may use in analysis, an approximation of the exact function, which may leads to the possible misleading results. Although, this problem can be alleviated by going to higher terms in the expansion, but adding any new term means introducing a new parameter that must be estimated. This indeed increases the divergences of the analysis. Moreover, since this method is based on the Taylor expansion of the scale factor, cosmography is restricted in the scope of Friedmann-Robertson-Walker metric. Hence it is interesting to find a more common kinematic approach which is applicable in other spacetimes posses inhomogeneous properties.
In next section we will derive an almost-general model independent solution from which one can reconstruct any CS parameter without any limitation and problems arising in usual cosmography.
3 An Exact Cosmographic Solution
In this section we find an exact analytical expressions between the Hubble and jerk parameters. It is interesting to note that in  an almost the same expression has been found by considering special parametrization for the jerk parameter, but this solution seems to be wrong as the authors consider a minus sign in the definition of jerk parameter (see eq (5) of this reference). This mistake affects other results. In what follows, first we derive a general differential equation which could be used for any parametrization of the jerk parameter and then we solve it for .
It is possible to obtain the derivatives of the squared Hubble parameter with respect to as follows 
which in turn, for constant 111Note that one can parameterize jerk parameter as to find a solution for non-constant which is in consistent with the CDM scenario at current time., gives the following general solution for the Hubble parameter
where . Requiring the consistency of (8) at gives
Obviously eq (8) is a model-independent solution. Generally, this solution corresponds to a cosmological model without radiation component 222One can put eqs (6a) and (6b) in eq (6c) to find a similar differential equation as eq (7) in term of snap parameter. Doing so, probably, recovers radiation term.. Note that, as we have shown in next section there is degeneracy between curvature and . It is worth to mention that for spatially flat cosmological constant dark energy model the jerk parameter is . Therefore, jerk parameter has been a traditional tool to test the spatially flat CDM model. From eq (8) it is clear that considering recovers flat CDM model. Hence, in this case, one can consider and as matter and dark energy density parameters respectively.
where , and . In the same manner one can reconstruct any other CS parameters.
Although one can use the following dimensionless Hubble parameter
3.1 Recovering CDM Model
In this section we consider that cosmological constant plays the role of dark energy. Therefore, one can write Friedmann equation as
Taking we obtain
Since, in general, , we can rewrite above equation as follows
Therefore, we can evaluate the current value of curvature parameter only in term of jerk parameter as
This equation clearly shows that any deviation from is an evidence of non-flat universe. Therefore we may argue that for all models consider cosmological constant as dark energy, the geometry of spacetime must be flat. Moreover, when a time varying dark energy is responsible for current Universe accelerating expansion.
which clearly shows degeneracy between and . In view of eq (16), we may also consider A as matter density () and B as dark energy density () if our estimations indicate , otherwise we consider .
In the next section we will use observational Hubble data (OHD), BAO/CMB data and Pantheon compilation and their joint combination to constrain cosmographic CDM model with following parameters space
where is the deceleration-acceleration transition redshif (see subsec. 5.1 for more details). Note that must be greater or equal to zero, this in fact imposes a certain lower limit on jerk parameter as . We have to consider this point in our estimations.
4 Data and Method
In this section we briefly describe the astronomical data and the statistical method we have been used to constrain parameter set (17).
Type Ia Supernovae: We adopt the Pantheon compilation  containing 1048 SNIa apparent magnitude measurements over the redshift range of , which includes 276 SNIa () discovered by the Pan-STARRS1 Medium Deep Survey and SNIa distance estimates from SDSS, SNLS and low-zHST samples. It is also possible to use the JLA dataset  which combines the SNLS and SDSS SNe to create an extended sample of 740 SNe to reduce the estimation time, but we found that using the Pantheon data slightly improves the parameter estimations. In this case the chisquare is defined as
where is the predicted distance modulus given by
and is the inverse of the by Pantheon compilation covariance matrix (). It is worth mentioning that since the parameter is only an additive constant, thus, marginalizing over does not affect the
Observational Hubble Data: we use OHD data from Table 2 of Ref  which includes datapoints in the redshift range . This measurements are uncorrelated and determined using the cosmic chronometric (CC) technique. It is worth nothing that the OHD data can be categorized into the following two classes, (1) BAO based data and (2) cosmic chronometric (CC) based data. To obtain OHD data from BAO, we usually model the redshift space distortions and assume an acoustic scale in a specific model. Therefore, this class of data is model-dependent and hence cannot be used for constraining a cosmological model. Nonetheless, to determine the CC data we use the most massive and passively evolving galaxies based on the âgalaxy differential age â method. consequently, this class of OHD data is model-independent (see ref for more details). Since this compilation includes uncorrelated data, we have as covariance matrix for this class of data. For OHD data the chisqure is given by
BAO/CMB Data: to obtain the BAO/CMB constraints on the model parameters we adopt 1- for BAO we consider six data points (see Table. 1) obtained from the WiggleZ Survey , SDSS DR7 Galaxy sample  and 6dF Galaxy Survey  datasets, 2- for CMB, our considered measurement is derived from the WMAP7 observations  (for more details about methodology of obtaining the BAO/CMB constraints on model parameters see ). Recently, Mamon et al  have used this data to place constraints on a reconstructed dark energy model. The chisqure of this data is given as
where is the co-moving angular-diameter distance and is the dilation scale given by
Also is the decoupling time and the inverse of covariance matrix for this data is
In above equations, is defined as
where for respectively.
Since these three datasets are independent, the total chisqure could be written as . Therefore, we evaluated the following total likelihood
Moreover, to check the degeneracy direction between computed parameters we perform covariance matrix which could be obtained from our MCMC runs. Theoretically, covariance matrix is defined as
where the uncertainties in parameters and are given by and are the respectively, and is the correlation coefficient between and .
We use Metropolis-Hasting algorithm to generate MCMC chains for all parameters. For each parameter we run 4 parallel chains with 6000 separate iterations to stabilize the estimations. We perform Gelman-Rubin and Geweke tests to confirm the convergence of MCMC chains. We also confirm the convergence of all chains by monitoring the trace plots for good mixing and stationarity of the posterior distributions. In our Baysian estimations, we assume the following uniform priors for free parameters of the model (17):
In Table. 2 we have listed our statistical analysis on parameter space (17) using OHD, BAO/CMB, SNIa, and their joint combination dataset at 1 error. It is worth mentioning that, although, SNIa data by itself is not sensitive to the universe expansion rate , but in the joint analysis, this data constrains other parameters of the model which in turn affect the computation of . That is why in Table. 2 we observe a change in the value of when fitting model to the joint OHD+BAO/CMB+SNIa data. In Table .(3) we have compered our obtained for OHD, BAO/CMB and OHD+BAO/CMB+SNIa data to those obtained by other researchers. Results of this table clearly show that when we use joint dataset our computed is in high agreement with those obtained by Sievers et al () , Chen et al () , and J. Dunkley et al () . Th estimated obtained from fitting to OHD data is in excellent agreement with those of Chen & Ratra ()  and Chen et al () . When we constrained model over BAO/CMB data alone, the computed value of is obtained in good agreement with what reported by Riess et al . We have also compared our estimated deceleration parameter to those of Aviles et al (columns forth and fifth of Table 1 ), Muthukrishna & Parkinson (row third of Table 3 ), and zhang et al (row forth of Table 1 ) in Table .(4).
|Parameter||OHD (CC)||BAO/CMB||SNIa (Pantheon)||CC + ABO/CMB + Pantheon|
It is worth noting that while the estimated values of reported in Ref  almost not physical (taking high values), our computed values of these parameters are physical for OHD, BAO/CMB, SNIa and their joint datasets. When we fit to OHD data alone we obtain and , these results are in excellent agreement with those obtained from 9years WMAP . Also when we fit to SNIa or OHD+SNIa data alone, we obtain and , these results are in excellent agreement with those obtained from Planck (2015) collaboration . The estimated valus of & obtained from fitting model over BAO/CMB data are also in good agreement with those of ref . We have compared the computed spatial curvature to those of 9 years WMAP , Planck (2015) collaboration and Park& Ratra  in Table. 5. From Table. 2 we observe that while individual OHD (CC) or BAO/CMB data predict (estimate) a slightly open universe (), SNIa (Pantheon) and OHD+BAO/CMB+SNIa predict a slighly closed universe (). Moreover, from this table, we see that OHD (BAO/CMB) data alone put tighter constrain on this parameter. It is worth nothing a stringent test of eternal inflation could be provided by constraining at around the level [59, 60, 61]. In fact, bellow , cannot be decisively distinguished from primordial fluctuations . However, large-scale anomalies as well as some Inflationary scenarios tend to have observable levels of spatial curvature [63, 64].
The contour plots, at , , and confidence levels, of the parameter space for joint OHD+BAO/CMB+SNIa are depicted in Figure. 1. We have also depicted the robustness of our fits for in Figures. 3, 3. Figure. 4 shows the variations of at 68% and 95% error for these data and their joint combination. From these figures it is also clear that our exact method gives raise to much better and tighter constrains on CS parameter with respect to the previous works.
|Ade et al (Planck 2015)||(at )|||
|Chen & Ratra||(at )|||
|Sievers et al||(at )|||
|Gott et al||(at )|||
|J. Dunkley et al (CMB)||(at )|||
|Aubourg et al (BAO)||(at )|||
|V. Lukovic et al||(at )|||
|Chen et al||(at )|||
|Riess et al||(at )|||
|Present work||(at )||for OHD|
|Present work||(at )||for BAO/CMB|
|Present work||(at )||for CC+BAO/CMB+Pantheon|
|Aviles et al|||
|Muthukrishna & Parkinson|||
|Zhang et al|||
|Researchers||9yearsWMAP||Planck (2015)||Park & Ratra||This work|
5.1 Transition Redshift
It is well known that the universe expansion phase has changed from decelerating to accelerating at a specific redshift called ‟Transition redshift ˝[65, 66]. Mathematically we can find this decelerating-accelerating redshift by imposing in (10a). Doing so, and after some algebra we find the following general transition redshift.
which can be simplified and rewritten in terms of CS parameters as
Note that we may set and in above equation. According to the previous discussion, for CDM model, we can rewrite eq 30 and relate transition redshift to the spacetime curvature as follows
Our statistical analysis on transition redshift for both datasets and their joint combination could be seen in Table. 6. Figures. 7 and 7 depict the robustness of our fits for . Gastri et al  have recently have used different SNIa data in combination with BAO/CMB data to constran transition redshift. For MLCS2k2+BAO/CMB (see for details of MLCS2k2 data) they obtained and for SALT2+BAO/CMB (see  for details of SALT2 data) they found . Also Farooq et al  used 38 data and found for two Hubble constant priors as at error. Moreover, recently Capozziello et al  through an effective cosmographic construction, in the framework of gravity have obtained at error. Note that in ref  authors set (from Plank (2015)) to obtaine above mentioned transition redshift but when they consider this parameter as a free one they obtain which is not accurate value. Generally, our obtained transition redshifts, except for BAO/CMB data, are consistent with what is expected in the cosmological models with present-epoch energy budget dominated by dark energy as well as standard spatially flat CDM model.
Cosmpgraphy is based on the Taylor expansion of the scale factor around . This expansion is the source of two main problems arise in this kinematic approach of the study of universe. In stead of expanding scale factor we have combined the mathematical definitions of deceleration and jerk parameters (see eqs 2b,2c) which results in a general second order differential equation for squared Hubble parameter. Although the solution of this equation could lead to a general function for Hubble parameter, but it seems to be much complicated. However, it is possible to find some reasonable solutions by considering jerk parameterization. In this paper we assumed and found a Hubble function in terms of jerk parameter. Using this function we have reconstructed other cosmographic parameters as well as deceleration-acceleration transition redshift. It is worth mentioning that since our approach is totally model-independent, we can constrain any derived quantity without any doubt on the validity of the estimations. Next we recovered CDM model from the obtained Hubble function. It is found that when cosmological constant is responsible for the current cosmic accelerating expansion the geometry of spacetime is necessarily should be flat. For any other dynamical (time varying) dark energy scenarios the spacetime geometry should be non-flat. Probably this approach could be used for other spacetime posses some inhomogenities. We have constrained cosmpgraphic parameters as well as transition redshift and spacetime curvature over observational Hubble data , BAO/CMB  and SNIa (Pantheon compilation) , and their joint combination. Our results are in agreement with almost all other results such as 9years WMAP and Planck (2015) collaboration.
We are grateful to professor Tamara Davis for critical review of the manuscript prior to submission.
-  B. P. Schmidt et al, Astrophys. J 507, 46 (1998).
-  S. Perlmutter et al, Nature 391 51 (1998).
-  A. G. Riess et al, Astron. J 116, 1009 (1998).
-  G. Hinshaw et al, Astrophys.J.Suppl.Ser 208, 25 (2013).
-  V. Sahni and A. Starobinsky, Int. J. Mod. Phys. D 15, 2105 (2006).
-  E. W. Kolb, S. Matarrese, A. Notari, and A. Riotto, Phys. Rev. D 71, 023524 (2005).
-  S. Rasanen, JCAP 0611, 003 (2006).
-  J. D. Barrow and S. Cotsakis, Phys. Lett. B 214, 515 (1988).
-  G. Dvali, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000).
-  T. Padmanabhan, Phys. Rep 188, 285 (1990).
-  E. Komatsu, et.al, Astrophys. J. Supp 192, 18 (2011).
-  S. Weinberg, Rev. Mod. phys 61, 1 (1989).
-  S. M. Carroll, Living Rev. Rel 4, 1 (2001).
-  C. Cattoen and M. Visser, Phys.Rev. D 78, 063501 (2008).
-  S. Weinberg, Cosmology, Oxford Univ. Press, Oxford (2008).
-  M. Visser, Gen. Rel. Grav 37, 1541 (2005).
-  C. Cattoën and M. Visser, Class. Quantum Gravity 24, 5985 (2007).
-  A. Aviles, C. Gruber, O. Luongo, H. Quevedo, Phys. Rev. D 86, 123516 (2012).
-  S. Capozziello, R. Lazkoz, and V. Salzano, Phys. Rev. D 84, 124061 (2011).
-  M. -J. Zhang , H. Li and J.-Q. Xia, Eur. Phys. J. C 77, 434 (2017).
-  A. Salehi, M. R. Setare and A. Alaii, Eur. Phys. J. C 78, 495 (2018).
-  D. Muthukrishna and D. Parkinson, JCAP 11, 052 (2016).
-  F. Y. Wang, Z. G. Dai and Shi Qi, A& A 507, 53 (2009).
-  F. Y. Wang and Z. G. Dai, A& A 536, A96 (2011).
-  S. Capozziello, R. D́Agostino and O. Luongo, MNRAS 476, 3924 (2018).
-  A. Aviles, J. Klapp, and O. Luongo, arXiv: 1606.9195 (2017).
-  A. Aviles, A. Bravetti, S. Capozziello and O. Luongo, Phys. Rev. D 90, 043531 (2014).
-  C. Gruber and O. Luongo, Phys. Rev. D 89, 103506 (2014).
-  S. Capozziello, O. Farooq, O. Luongo and B. Ratra, Phys. Rev. D 90, 044016 (2014).
-  S. Capozziello, V.F. Cardone and V. Salzano, Phys. Rev. D 78, 063504 (2008).
-  A. Gómez-Valent and L. Amendola, JCAP 1804, 051 (2018).
-  A. Gómez-Valent, arXiv: 1810.02278 (2018).
-  J. Ryan, S. Doshi and B. Ratra, MNRAS 480, 759 (2018).
-  D. M. Scolnic, et al., Astrophys. J 859, 101 (2018).
-  R. Goistri et al., JCAP 03, 027 (2012).
-  M. Visser, Class.Quant.Grav 21, 2603 (2004).
-  V. Sahni, T.D. Saini, A. A. Starobinsky, and U. Alam, JETP Lett 77, 201 (2003).
-  X. Lixin, L. Wenbo and L. Jianbo, JCAP 07, (2009).
-  Yu.L. Bolotin, V.A. Cherkaskiy, O.A. Lemets, D.A. Yerokhin, L.G. Zazunov: arXiv:1502.00811
-  M. Betoule, et al, Astron. Astrophys 568, A22 (2014).
-  M. Moresco et al., JCAP 05, 014 (2016).
-  C. Blake et al., Mon. Not. R. Astron. Soc 41, 1707 (2011).
-  W. J. Percival et al., Mon. Not. R. Astron. Soc 401, 2148 (2010).
-  F. Beutler et al., Mon. Not. R. Astron. Soc 416, 3017 (2011).
-  N. Jarosik et al., Astrophys. J. Suppl 192, 14 (2011).
-  R. Goistri et al., JCAP 03, 027 (2012).
-  A. A. Mamon, K. Bamba and S. Das, Eur. Phys. J. C 77, 29 (2017).
-  G. Chen and B. Ratra, PASP 123, 1127 (2011).
-  J. L. Sievers et al, JCAP 1310, 060 (2013).
-  G. Chen, S. Kumar and B. Ratra, Astrophys. J 835, 86 (2017).
-  J. Dunkley et al, Astrophys. J 739, 52 (2011).
-  P. A. R. Adel et al, A & A 594, A13 (2016).
-  J. R. Gott, M. S. Vogeley, S. Podariu and B. Ratra, Astrophys. J 549, 1 (2001).
-  E. Aubourg et al, Phys. Rev D 92, 123516 (2015).
-  V. V.Lukovic, R. D́Agostino and N. Vittorio, A & A 595, A109 (2016).
-  A. G. Riess et al, Astrophys. J 826, 56 (2016).
-  G. Hinshaw et al., Astrophys.J.Suppl.Ser 208, 25 (2013).
-  C. -G. Park, B. Ratra, arXiv: 1803. 05522 (2018).
-  M. Kleban and M. Schillo, JCAP 2012, 029 (2012).
-  A. H. Guth and Y. Nomura, Phys. Rev. D 86, 023534 (2012).
-  A. H. Guth, D. I. Kaiser, and Y. Nomura, Phys. Lett. B 733, 112 (2014).
-  C. D. Leonard, P. Bull and R. Allison, Phys. Rev. D 94, 023502 (2016).
-  G. Aslanyan and R. Easther, Phys. Rev. D 91, 123523 (2015).
-  M. C. Johnson and W. Lin, JCAP 03, 051 (2016).
-  J. V. Cunha, Phys. Rev. D 79, 047301 (2009).
-  M. J. Mortonson, W. Hu, and D. Huterer, Phys. Rev. D 80, 067301 (2009).
-  S. Jha, A. G. Riess and R. P. Kirshner, Astrophys. J 659, 122 (2007).
-  J. Guy et al., A & A 466, 11 (207).
-  O. Farooq O, F. R. Madiyar, S. Crandall and B. Ratra, Astrophys. J. 835, 26 (2017).
-  S. Capozziello, O. Luongo and E. N. Saridakis, Phys. Rev. D 91, 124037 (2015).