Curvature dark energy reconstruction through different cosmographic distance definitions
In the context of gravity, dark energy is a geometrical fluid with negative equation of state. Since the function is not known a priori, the need of a model independent reconstruction of its shape represents a relevant technique to determine which model is really favored with respect to others. To this aim, we relate cosmography to a generic and its derivatives in order to provide a model independent investigation at redshift . Our analysis is based on the use of three different cosmological distance definitions, in order to alleviate the duality problem, i.e. the problem of which cosmological distance to use with specific cosmic data sets. We therefore consider the luminosity, , flux, , and angular, , distances and we find numerical constraints by the Union 2.1 supernovae compilation and measurement of baryonic acoustic oscillations, at . We notice that all distances reduce to the same expression, i.e. , at first order. Thus, to fix the cosmographic series of observables, we impose the initial value of by fitting through supernovae only, in the redshift regime . We find that the pressure of curvature dark energy fluid is slightly lower than the one related to the cosmological constant. This indicates that a possible evolving curvature dark energy realistically fills the current universe. Moreover, the combined use of and shows that the sign of the acceleration parameter agrees with theoretical bounds, while its variation, namely the jerk parameter, is compatible with . Finally, we infer the functional form of by means of a truncated polynomial approximation, in terms of fourth order scale factor .
pacs:98.80.-k, 98.80.Jk, 98.80.Es
Modern cosmology is nowadays plagued by several shortcomings which jeopardize the current understanding of universe dynamics. Particularly, these problems may suggest to reconsider the standard approach of gravitation, based on Einstein’s gravity, in favor of alternative theories of gravity. Alternative gravity pictures have been extensively introduced in order to describe universe dynamics without the need of additional material ingredient as dark energy and dark matter described by new particles at fundamental level. On the other hand, the simple introduction of a cosmological constant vacuum energy seems inadequate to characterize the whole universe evolution at any epoch adj (). Thus, these alternative theories are viewed as a bid to reformulate in toto semi-classical schemes where General Relativity is only a particular case of a more extended theory. In particular, such theories are able to extend General Relativity predictions by means of higher order curvature invariants. Other pictures assume extensions based on minimally or non-minimally coupled scalar fields in the gravitational Lagrangians review1 (); review2 (); review3 (); review4 (); review5 (). Furthermore, Einstein gravity can be extended by carrying out the full Mach principle: this fact leads to the introduction of a varying gravitational coupling. Under these hypotheses, the Brans-Dicke theory BD () represents the prototype of alternative schemes to General Relativity. It naturally includes a variable gravitational coupling, whose dynamics is governed by a single scalar field non-minimally coupled to geometry BD (); cimento (); odintsov (). From one hand, extensions of General Relativity are therefore able to describe the above-mentioned theoretical aspects. On the other hand, it is also possible to account every unification scheme of fundamental interactions, such as superstring, supergravity, or grand unified theories and physically by low-energy effective actions containing non-minimal couplings or higher order curvature terms sciama (). In fact, interactions between quantum scalar fields and background geometry, or gravitational self-interactions, naturally yield such corrections to the Einstein-Hilbert’s Lagrangian birrell (). Hence, it is easy to show that several geometrical corrections are inescapable within quantum gravity effective actions and allow consistent pictures close to Planck scales vilkovisky (). These schemes represent working approaches towards a self consistent quantum picture, giving rise to interesting consequences once corrections like , , , , , or are involved. A crucial fact is that alternative theories may provide analogies with the effective string or Kaluza-Klein Lagrangians, when compactification mechanisms of extra spatial dimensions are imposed veneziano ().
A consequence of such extended theories of gravity is the possibility to frame current universe dynamics in a self consistent way considering their infrared counterpart. In particular, such models may address the problem of current universe speed up sn () by considering further gravitational degrees of freedom. Indeed, General Relativity seems not capable of dealing with present cosmic acceleration, unless an unknown fluid dubbed dark energy is added to the standard matter fluid energy-momentum tensor. At late times, the fluid responsible for accelerating the universe dominates over all other contributions, driving the universe evolution. It should be able to reproduce current observations gala (). Consequently, the dark energy equation of state behaves anti-gravitationally by counterbalancing gravitational attraction ioequew (). Thus, in this concordance model, the universe dynamics is described through pressureless matter terms, i.e. the sum of baryons and cold dark matter, through a evolving barotropic dark energy contribution and a vanishing spatial curvature qualcosa1 (); qualcosa2 (); qualcosa3 ().
A straightforward way to address geometrically the problem of dark energy is by the so called gravity, where is a generic function of the Ricci scalar review3 (); review4 (); review5 (). In this paper, we fix constraints on geometrical dark energy fluid inferred in the context of gravity. To this end, we adopt cosmography to fix cosmological bounds on the function and its derivatives at low redshift regime where degeneracy of concurring dark energy models is more evident. Cosmography allows to determine cosmological constraints in a model-independent way, once scalar curvature is somehow fixed. The idea is to expand into Taylor series cosmological observables. These expansions can be compared with data to get the cosmographic series, i.e. the numerical bounds on scale factor derivatives mio (); salzano1 (); salzano2 (). One commonly-used technique is represented by expanding the luminosity distance and compare it with supernovae data. However, a degeneration problem (duality problem) occurs once different cosmological distances are involved. Hence, a non-definitive consensus exists on the adequate cosmological distance to use in the framework of cosmography. We therefore perform the experimental analysis by means of three cosmological distance rulers, i.e. luminosity, flux and angular distances. We check the viability of different cosmological distances and measure cosmological constraints on the cosmographic series, deriving bounds on curvature dark energy.
The paper is organized as follows: in Sec. II, we highlight the main features of cosmography and its application to cosmology. In Sec. III, we describe the problems related to cosmography, pointing out the so called duality problem. In Sec. IV, the experimental procedures is described. In Sec. V, cosmography in view of is discussed . Finally, Sec. VI is devoted to conclusions and perspectives.
Ii Basics of Cosmography
Let us summarize the main aspects of cosmography and describe how it can be considered as a tool to fix constraints on cosmological observables. Firstly, let us assume that the cosmological principle holds and the equation of state is currently determined by a geometrical fluid, with pressure . Under these hypotheses, we expand cosmological observables into Taylor series and match the derivatives of such expansions with cosmological data. Examples of expanded quantities are the Hubble parameter, the luminosity distance, the apparent magnitude modulus orl (); orl2 (), the net pressure, and so forth orl21 (); orl99 (). The power series coefficients of the scale factor expansion are known in the literature as cosmographic series (CS), if calculated at present time, or alternatively at the redshift . Those coefficients are therefore expressed in terms of the cosmological scale factor and its derivatives mio (). It follows that the cosmographic approach does not need to assume a particular cosmological model.
Thence, one of the main advantage of cosmography is alleviating degeneracy among cosmological models, i.e. cosmography allows, in principle, to understand which model better behaves than others. In case of gravity, for example, matter density degenerates with scalar curvature and cannot be constrained a priori. However, cosmography fixes model independent constraints on the cosmological equation of state and then results a technique to discriminate among competing models mio (), removing the degeneracy between matter and scalar curvature salzano1 (). This technique turns out to be useful to reconstruct the form of which better traces the universe expansion history. Thus, more precisely, cosmography represents a model independent method to infer cosmological bounds, once spatial curvature is somehow fixed.
Recent observations point out that the scalar curvature is negligible, so we can easily impose ob (). Than one has
which represents the Taylor series of the scale factor , around . The CS can be thus defined as
Here, dots represent derivatives with respect to the cosmic time . Each term brings its own physical meaning. Particularly, the Hubble rate is intimately related to the variation of with time, the acceleration parameter measures how the universe is speeding up and the jerk parameter permits one to understand how the acceleration varied in the past. The coefficients are defined as
and are considered at a given time . We may argue that such quantities are able to describe the kinematics of the universe turnercosmografia () and we do not pose, at this stage, the problem of which model causes the universe acceleration. In analogy to the classical mechanics, we say that cosmography is a kinematic approach to trace the universe expansion today. From one hand, the advantages of cosmography consist on its model independent reconstructions of present-time cosmology. In other words, it can be considered like a snapshot of the today observed universe capable of giving initial conditions for reconstructing back the cosmic evolution. From the other hand, the disadvantages rely on the fact that current data are either not enough to guarantee significative and accurate constraints or do not fit significant intervals of convergence for . In addition, the cosmological observable that one expands into Taylor series, i.e. is not known a priori. Consequently, there is no physical motivations to use a particular cosmological distance than others. This means that the use of a given luminosity distance to constrain CS is only motivated by ad hoc arguments. This fact constitutes the so-called duality problem that we discuss in the next section. To alleviate duality problem, we will compare three different cosmological distances to trace universe expansion history at late times, under the hypothesis of a geometrical dark energy fluid.
Iii The duality problem and cosmographic convergence
By a cosmographic analysis, one can fix constraints on the geometrical dark fluid, alleviating the degeneracy problem. To this end, one needs a self-consistent definition of causal distance. Unfortunately, standard definitions implicitly postulate that the universe is accelerating lix1 (); lix2 (), i.e. to infer the distance expansion, we evaluate the distance that a photon travels from a light source at to our position at , defined as . Consequently, one obtains as prototype the so called luminosity distance , while other definitions, e.g. the photon flux distance , angular diameter distance and so forth, can easily be derived from different considerations. As previously stressed, this leads to a severe duality problem on the choice of the particular cosmological ruler to use for fixing cosmological constraints on the CS.
Here, we use three different cosmological distances as rulers, e.g. the luminosity, flux and angular distances, , and respectively. Below the definition of these distances is reported in terms of , that is
These distances can be used to the fix causal constraints on the curvature fluid in order to alleviate the degeneracy problem. For the sake of clearness, it is important to stress that although is associated to the ratio of the apparent and absolute luminosity of astrophysical objects, the other distances, i.e. and , may be also used to fix bounds on the observable universe. All the different cosmological distances rely on the fundamental assumption that the total number of photons is conserved at cosmic scales. Hence, there is no reason to discard one distance with respect to another since all of them fulfill this condition. The duality problem represents a not well understood issue of observational cosmology bll (). In this work, we find differences in fitting Eqs. (4), showing that there is no reason to adopt only as the only cosmological distance.
However a problem of convergence may occur, leading to possible misleading results for in the cosmographic Taylor series. An immediate example is due to the most high supernova redshift in a typical data set. Usually, one has that the furthest redshift at approximatively , showing that a few number of supernovas spans in the range . It follows that numerical divergences and bad convergences may occur in the analysis, since Taylor expansions are carried out around . A plausible landscape deals with introducing alternative redshift definitions, re-parameterizing the cosmological distances in a tighter redshift range gtm (). These possible re-parameterizations must fulfill the conditions that the distance curves should not behave too steeply in the interval . Moreover, the luminosity distance curve should not exhibit sudden flexes, being one-to-one invertible as discussed in mio (). In other words, it is easy to show that the new redshift re-parameterization, i.e. , provides , with a generic function of the redshift , with the property , as . In this work, we describe a technique to reduce the convergence problem, calibrating cosmological distance at first order in the Taylor series within a smaller range of redshift. Our strategy is to fix with supernovae in the range . This turns out to be useful since a wide range of data is actually inside the sphere and all cosmological distances at first order reduce to
where represents the generic distance, i.e. . Once is fixed, the series naturally converges better since its shape increases or decreases as decreases or increases respectively. In other words, the dynamical shape of any cosmological curve depends on the value given to . As is somehow known, curves behave better at higher redshift, alleviating convergence problems as expected. These arguments represent a further tool in order to fix model independent constraints on and its derivatives. Indeed, is fixed regardless the cosmological distance taken into account, by means of Eq. (5). It is possible to fix in the range with supernovae only. We find
In cosmography, the strategy of fixing in a smaller interval of data overcomes several problems associated to the well consolidated usage of auxiliary variables. Indeed, as above mentioned, the method of adopting auxiliary variables consists in determining parametric functions in terms of the redshift , whose values rely in the interval , as and respectively. This procedure rearranges catalog data and suffers from severe shortcomings gtm (); gtm2 (). Indeed, the form of is not known a priori and any possible reparameterized variable should guarantee that errors do not deeply propagate in the statistical analysis. In several cases, variables are therefore inconsistent with low redshift cosmography, providing misleading results, albeit their use becomes more relevant for high redshift data sets.
In our case, we propose to fix as a low redshift cosmographic setting value, since all distances reduce to Eq. (5) at a first order of Taylor expansions. Our corresponding best fit intervals are compatible with previous analysis nfn21 () and guarantee that errors do not significatively propagate on measured coefficients. For our purposes, the strategy of fixing by means of small redshift data only better behaves than standard auxiliary variables, due to the fact that (see for recent applications nfn22 (); ratra (); ratra1 ()), although it would fail at higher redshift domains. Since, in our cases, reaches the upper value , i.e. the maximum of the supernova compilation, we expect that the use of Eq. (5) to get would guarantee refined best fit results with respect to any possible reparameterized auxiliary variables.
An additional technique is to combine more than one data set to infer cosmological bounds. Indeed, although we treat bad convergence of truncated series by using Eq. (6), further data sets would improve the quality of numerical estimates. Hence, the CS and the derived constraints on derivatives would improve consequently. Here, we combine supernovae data with the baryonic acoustic oscillation measurement. It is possible to show that such a choice actually reduces the convergence problem.
Iv Cosmological data sets and the fitting procedure
In this section, we describe the two data sets used for the cosmographic analysis. First, let us consider the Union 2.1 compilation kow (). Second, we assume the measurement of baryonic acoustic oscillation (BAO) sn32 (). As it is well known, supernovae data span in the plane , consisting of 580 supernovae, in the observable range . For our purposes, to fix viable constraints, we follow a standard Bayesian analysis, dealing with the determination of best fits, evaluated by maximizing the likelihood function . Here, is the (pseudo)-squared function, or reduced squared. The distance modulus for each supernova is
and once given the corresponding error, we are able to minimize the square as follows
On the other hand, the large scale galaxy clustering observations provide the signatures for the baryonic acoustic oscillation. This gives a further tool to explore the parameter space and alleviate the convergence problem. We use the peak measurement of luminous red galaxies observed in Sloan Digital Sky Survey (SDSS). By employing as the measured quantity, we have
with . In addition, the observed is estimated to be , with an error . In the case of the BAO measurement, we minimize the squared
An important feature of BAO is that it does not depend on .
Estimations of the cosmographic parameters may be performed passing through the standard Bayesian technique, maximizing the likelihood function:
where is explicitly determined for each compilations here employed and the subscript indicates the data set, i.e. supernovae or BAO. Maximizing the likelihood function is equivalent to minimizing the total -squared function and so one argues to maximize
In particular, the cosmographic results have been obtained by directly employing Eq. (12) maximizing the corresponding likelihood functions over a grid, through a standard Bayesian analysis. In so doing, the cosmographic series has been evaluated and the provided errors refer to as the , associated to the of confidence level. Once all cosmographic coefficients have been determined through a direct Gaussian maximization of the likelihood function, we will infer the derived coefficients, , by simply propagating the errors through the well consolidated logarithmic method.
V cosmography vs redshift
Let us consider now the coefficients (evaluated in terms of the redshift ) as function of observable quantities. The robustness of calculations leads to the advantage of relating derivatives of at to experimental bounds, without assuming a priori a form of . Hence, by expanding the causal distances , and , through the definition of the scale factor in terms of redshift, , we get
Thus, by considering the redshift definition in terms of the cosmic time , we can rewrite as a function of , being
Hence, it is now easy to get and derivatives in terms of , and to evaluate the corresponding values at present time. We obtain
that, for present time, simply allow us to write
where derivatives with respect to are indicated. By assuming the presence of standard pressureless matter, ( and ), including baryons and cold dark matter, the modified Friedmann equation, which fixes the whole energy budget of the universe, easily reads
with the dynamical expression for , i.e. . Eq. (15) determines curvature corrections from which we infer the dark energy fluid. In other words, the fluid density, responsible for the cosmic speed up, could be rewritten as
with the corresponding curvature pressure
It is convenient, for our purposes, to work in terms of , i.e. the explicitly depending on the redshift . In so doing, we assume the functional dependence , and find , . As it is discussed in salzano1 (); ngo (), models have to evade Solar System tests for General Relativity. This means that gravitational coupling has to agree with the local observed value. This is possible if the conditions
hold. From Eqs. (18), numerical estimates are derived once the CS is known from model independent fitting procedures by means of Eqs. (4). The numerical results of our analysis are reported in Tab. I and II. They have been obtained by employing numerical outcomes for , determined by the use of Eq. (V).
The range of measurements are inferred through the logarithmic formula, once the CS is fitted and Eqs. (18) assumed.
The range of measurements are inferred through the logarithmic formula, once the CS is fitted and Eqs. (18) assumed.
Our results indicate a slightly lower dark energy pressure than the one predicted by CDM model, which typically is . This fact suggests that any curvature dark energy fluid should behave differently than a pure (effective) cosmological constant. In general, numerical outcomes for does not completely agree with CDM but give compatible acceleration parameters for .
However, our results show that present data are not enough to completely constrain possible violations of the duality problem. Even though it is difficult to definitively remove such a problem, our approach allows one to alleviate it by dealing with the BAO measurement as a Gaussian prior in the likelihood definition. In so doing, one significantly reduces the allowed phase space, circumscribing Gaussian errors and reducing any possible systematics. This procedure has been accurately investigated also with different typologies of priors nfn21 (). In particular, through general analyses carried out by means of Monte Carlo Markov chains, it seems that more stringent intervals occur as the BAO measurement is included nfn23 (). Such tighter best fit values are however not incompatible with our results which rely in viable cosmographic intervals gtm (); gtm2 (); nfn21 (); nfn22 (). More stringent values will be carried out through Monte Carlo approaches, also in the field of gravity, as one can notice from first confirmations in ngo (). In addition, all cosmological distances indicate a positive, greater than 1, jerk parameter , although they provide unexpectedly smaller acceleration parameters for and distances.
The statistical strategy adopted to find out our cosmographic results permits one to immediately get the requested values for . Our results are actually statistically consistent, as one argues from the corresponding chi square functions and errors. Those results remarkably show that and increase, as the acceleration parameter decreases and confirm previous studies on the cosmographic series, providing intervals of confidence in agreement with theoretical predictions. However, although all square parameters are nearly the same for all the fits involved, the leading results determined from better support theoretical predictions, since seem to be more compatible with a higher acceleration parameters, as one expects from current data gala (); gtm (); gtm2 (); nfn24 ().
Summing up, the absolute values of each variable increases as one performs fits with and . In other words, and analyses show significative departures, incompatible to the values found by , for the CS, although they indicate compatible values of and and give reasonable agreement with .
Keeping in mind such considerations, we are therefore able to infer the functional form of , compatible with previous results. In particular, we get
where , with a truncating number, i.e. . The approach fulfills the cosmographic bounds and gives reasonable results for reconstructions (see also ngo ()). In particular, and the statistically favored dependence of on is an inverse series in terms of .
Vi Discussion and conclusions
In this paper, we present a method to fix constraints on curvature dark energy in the context of ) gravity by using a cosmographic analysis derived from different cosmological distances. Cosmography is here adopted to determine model independent bounds on observables related to and its derivatives. We follow this procedure by keeping in mind that each term of CS has a precise physical meaning, and allows to infer dynamical properties of curvature fluid. Rephrasing this point, the limits on CS allow to recover present time bounds on the dark energy dynamics, regardless the particular cosmological model. In so doing, to alleviate the duality problem among cosmological distances, we investigate cosmography by means of the luminosity, flux and angular distances. The approach consists in using cosmography to fix cosmological bounds on curvature fluid and to choose a combination of causal distances in order to alleviate degeneracy among cosmological models. Determining cosmographic results through the three mentioned distances leads to tighter intervals of measurements, showing that values of curvature dark energy pressure rely in lower bounds than the intervals expected in a pure CDM model. To this end, we used two data sets, i.e. the Union 2.1 compilation of supernovae and the baryonic acoustic oscillation measurements. In our analysis, we considered the today observed value of the Hubble parameter as a prior. This choice is due to the need of alleviating the convergence problem which plagues cosmography as soon as cosmological data exceed the range . Thus, we fix as a sort of initial condition on the three distances, by means of supernovae only, using all the distances expanded at the lowest order of , that is . The corresponding fitting value of has been obtained in the supernovae range where the convergence problem does not exist. This technique is developed in order to determine better the values of and , without introducing any re-parametrization variable of the redshift. In addition, we get constraints on derivatives up to the second order in the Taylor expansion around and found that the acceleration parameter is almost compatible when inferred from while the jerk and snap parameters are only compatible with respectively. These outcomes indicate again that the standard cosmological paradigm, namely the CDM model, may be extended in terms of curvature dark energy fluid. Finally, we proposed, as test function, a particular , able to reproduce the cosmographic bounds that fairly well fits the Hubble diagram up to fourth order expansion in powers of . In future research, we will refine this approach by combining more cosmological tests and further distances. This approach will allow to reduce the phase spaces and better constrain the free parameters of curvature dark energy.
- (1) Peebles, P. J. E., Ratra, B., Rev. Mod. Phys., 75, (2003).
- (2) Capozziello, S., Int. J. Mod. Phys. D 11, 483 (2002).
- (3) Nojiri S., Odintsov S.D., Int. J. Geom. Meth. Mod. Phys., 4, 115, (2007).
- (4) Capozziello, S., Francaviglia, M., Gen. Rel. Grav., 40, 357, (2008).
- (5) Capozziello, S., De Laurentis, M., Phys. Rept., 509, 167, (2011).
- (6) Nojiri, S., Odintsov, S. D., Phys. Rept., 505, 59, (2011).
- (7) Brans, C. H., Dicke, R. H., Phys. Rev., 124, 925, (1961).
- (8) Capozziello, S., de Ritis, R., Rubano, C., Scudellaro, P. Nuovo Cimento, 4, 19, (1996).
- (9) Buchbinder, I. L., Odintsov, S. D., Shapiro, I. L., Effective Action in Quantum Gravity, IOP Publishing, Bristol, (1992).
- (10) Sciama, D. W., Mon. Not. R. Ast. Soc., 113, 34, (1993).
- (11) Birrell, N. D., Davies, P. C. W., Quantum Fields in Curved Space, Cambridge University Press, Cambridge, (1982).
- (12) Vilkovisky, G. Effective action in quantum gravity, Class. Quant. Grav., 9, 895, 1992.
- (13) Gasperini, M., Veneziano, G., Phys. Lett. B, 277, 256, (1992).
- (14) Riess, A. G., et al., AJ, 116, 1009, (1998); Perlmutter, S., et al., ApJ, 517, 565, (1999).
- (15) Capozziello, S., De Laurentis, M. Luongo, O., Ruggeri, A. C., Galaxies 2013, 1(3), 216-260;
- (16) Luongo, O., Quevedo, H., arXiv:1005.4532, (2010).
- (17) Bamba, K., Capozziello, S., Nojiri, S., Odintsov, S. D., Astrophys. Space Sci., 342, 155-228, (2012).
- (18) Li, M., Li, X., Wang, S., Zhang, X., JCAP, 0906, 036, (2009);
- (19) Linder, E., Gen. Rel. Grav., 40, 23, (2008).
- (20) Copeland, E. J., Sami, M., Tsujikawa, S., Int. J. Mod. Phys. D, 15, 1753-1936, (2006).
- (21) Aviles, A., Gruber, C., Luongo, O., Quevedo, H., Phys. Rev. D, 86, 123516, (2012).
- (22) Capozziello, S., Cardone, V.F., Salzano, V., Phys.Rev. D 78, 063504 (2008).
- (23) Capozziello, S., Lazkoz, R., Salzano, V., Phys.Rev. D 84, 124061 (2011).
- (24) Luongo, O., Mod. Phys. Lett. A, 26, 20, 1459-1466 (2011).
- (25) Demianski, M., Piedipalumbo, E., Rubano, C., Scudellaro, P., Month. Not. of the Roy. Astr. Soc., 426, 1396-1415, (2012); Luongo, O., Iannone, G., Autieri, C., Europh. Lett. 90, 39001, (2010).
- (26) Lavaux, G., Wandelt, B. D., Astr. Jour., 754, 109, (2012); Aviles, A., Bonanno, L., Luongo, O., Quevedo, H., Phys. Rev. D, 84, 103520, (2011); Iannone, G., Luongo, O., Europh. Lett., 94, 49002, (2011).
- (27) Gruber, C., Luongo, O., Phys. Rev. D, 89, 103506, (2014).
- (28) Mennella, A., et al. A A, 536, 3, (2011); Verde, L., Protopapas, P., Jimenez R., ArXiv[astro-ph]: 1306.6766, (2013).
- (29) Neben, A. R., Turner, M. S., ApJ, 769, 133, (2013).
- (30) Nair, R., Jhingan, S., Jain, O., JCAP, 12, 028, (2012).
- (31) Holanda, R. F. L., Goncalves, R. S., Alcaniz, J. S., JCAP, 1206, 022, (2012).
- (32) Etherington, I. M. H., Gen. Rel. Grav., 39, 1055, (2007).
- (33) Aviles, A., Gruber, C., Luongo, O., Quevedo, H., arXiv[gr-qc]:1301.4044, (2013); Cattoen, C. Visser, M. Phys. Rev. D, 78, 063501, (2008); Luongo, O., Quevedo, H., Astroph. and Sp. Sci., 338, 2, 345-349 (2012); Cattoen, C., Visser, M., ArXiv[gr-qc]:0703122, (2007); Xu, L., Wang, Y., Phys. Lett. B, 702, 114, (2011).
- (34) Visser, M., Class. Quant. Grav., 21, 2603-2616, (2004).
- (35) Nesseris, S., García-Bellido, J., Phys. Rev. D, 88, 063521, (2013); Bochner, B., Pappas, D., Dong, M., ArXiv[astro-ph]:1308.6050, (2013).
- (36) Aviles, A., Bravetti, A., Capozziello, S., Luongo, O., ArXiv[gr-qc]:1405.6935, (2014).
- (37) Demianski, M., Piedipalumbo, E., Rubano, C., Scudellaro, P., Mon. Not. Roy. Astr. Soc., 426, 1396-1415 (2012).
- (38) Capozziello, S., Farooq, O., Luongo, O., Ratra, B., ArXiv[gr-qc]:1403.1421, (2014).
- (39) Komatsu, E., et al., Astrop. J. Sup., 192, 18, (2011).
- (40) Percival, W. J., et al., Mon. Not. Roy. Astron. Soc., 381, 1053, (2007).
- (41) Aviles, A., Bravetti, A., Capozziello, S., Luongo, O., Phys. Rev. D, 87, 044012, (2013).
- (42) Lazkoz, R., Alcaniz, J., Escamilla-Rivera, C., Salzano, V., Sendra, I., JCAP, 12, 005, (2013).
- (43) Jun-Qing; X., Vitagliano, V., Liberati, S., Viel, M., Phys. Rev. D, 85, 043520, (2012).