Model-independent reconstruction of non-metric gravity
Abstract
We consider gravity mediated by non-metricity, with vanishing curvature and torsion. The gravitational action, including an arbitrary function of the non-metric scalar, is investigated in view of characterizing the dark energy effects. In particular, we present a method to reconstruct the action without resorting to a priori assumptions on the cosmological model. To this purpose, we adopt a method based on rational Padé approximations, which provides a stable behaviour of the cosmographic series at high redshifts, alleviating the convergence issues proper of the standard approach. We thus describe how to reconstruct through a numerical inversion procedure based on the current observational bounds on cosmographic parameters. Our analysis suggests that the best approximation for describing the accelerated expansion of the universe is represented by a scenario with . Finally, possible deviations from the standard CDM model are discussed.
pacs:
04.50.Kd, 98.80.-k, 98.80.EsI Introduction
Several observational evidences support the standard model of cosmology, the so-called CDM model SN ; BAO ; Planck ; mio_vittorio . They provide the picture of a universe which recently entered a phase of accelerated expansion. In the framework of general relativity (GR), the simplest interpretation of dark energy, responsible for the late-time acceleration, is offered by the cosmological constant , with a negative constant equation of state lambda .
Despite the great success of CDM to explain current observations, the interpretation of at fundamental level is still far from being understood, mainly due to the fine-tuning problem inherent the energy of the vacuum fine-tuning . Among all scenarios attempting for an alternative explanation of dark energy effects (see e.g. rocco1 ; rocco2 ), particular attention in the last years has been devoted to modified gravity theories, including and theories review ; f(R) ; f(T) ; rocco3 .
An interesting possibility explored very recently is to consider the gravitational interaction mediated by the non-metricity, while curvature and torsion are vanishing Jimenez ; f(Q)_1 ; f(Q)_2 ; Sahoo . This approach can be extremely important to describe gravity at fundamental level because gravity can be dealt as a gauge theory not requiring a priori the validity of the Equivalence Principle.
In this context, investigating the theories, where is the non-metricity scalar, could offer new insights on the cosmic speed up deriving from the intrinsic implications of a different geometry with respect to the Riemannian one.
Standard cosmological approaches to extended theories of gravity rely on assuming a specific form of the relative action, and then analyzing the cosmic behaviour to check for possible deviations from GR. This procedure, however, may induce an a priori bias that could lead to misleading conclusions.
These considerations motivated to investigate new model-independent techniques to tackle the dark energy problem. In such a context, the cosmographic approach represents a powerful tool to study the kinematic features of the universe starting from first principles Weinberg72 ; cosmography1 . Indeed, the main advantages of cosmography reside in the fact that it does not need the assumption of a specific underlying cosmology to describe the dark energy behaviour cosmography2 . Applied to modified gravity theories, this method provides the remarkable opportunity to test the validity of GR on cosmological scales and to deal with possible deviations from the Einstein theory cosmography_extended .
In the present paper, we intend to apply the cosmographic method in order to obtain a model-independent reconstruction of gravity. After this Introduction, Sec. II gives a brief summary of non-metricity gravity and cosmology. In Sec. III, we start by recalling the main ingredients of the cosmographic technique and, subsequently, we illustrate the numerical procedure to reconstruct the action based on the method of rational Padé approximations. Finally, Sec. IV is dedicated to a general discussion of our results and concluding remarks.
Throughout this paper, we adopt natural units .
II Cosmology in gravity
In order to explore the cosmological features of non-metric gravity, let us consider the most general form of the affine connections Jarv18 :
(1) |
Here, the Levi-Civita connection of the metric tensor is given as
(2) |
whereas
(3) | ||||
(4) |
are the contortion and disformation tensors, respectively. In the above definitions, we have also introduced the torsion tensor, , and the non-metricity tensor, given by
(5) |
Thus, particular choices on the connections will specify the metric-affine spacetime. n our analysis, we assume that both curvature and torsion are vanishing, so that geometry is given by non-metricity. The non-metricity tensor is characterized by two independent traces, namely and depending on the order of contraction. Hence, one can define the non-metricity scalar as Jimenez
(6) |
which is a quadratic combination, invariant under general diffeomorphisms.
In analogy to studies on torsionless gravity and curvature-free gravity, we can generalize the -gravity to theories containing an arbitrary function of the non-metricity scalar, i.e. . Therefore, we consider the following action:
(7) |
where is the determinant of the metric, and is the Lagrangian density of the matter sector.
We note that, for , the above action is equivalent to the Einstein-Hilbert action up to a total derivative. In the case of a globally vanishing affine connections, the non-metricity tensor depends on the metric only and Einstein’s GR action is recovered. This occurs under the choice of the coincidence gauge, in which the origin of spacetime and that of the tangent space coincide (we refer the reader to Beltran19 for the details).
The gravitational field equations are then obtained by varying the action with respect to the metric, leading to
(8) |
where and is the energy-momentum tensor, defined as
(9) |
For cosmological purposes, we consider the line element , corresponding to the spatially flat Friedman-Lemaître-Robertson-Walker (FLRW) metric, in which is the cosmic scale factor used to define the Hubble expansion rate .
Focusing our attention on the coincidence gauge, we can write the modified Friedman equations as Beltran20
(10) | ||||
(11) |
where and are the total density and pressure of the cosmic fluid, respectively. Moreover, we have the following relation:
(12) |
which will play a central role in our reconstruction procedure. It is worth to emphasize that Eqs. 10, 11 and 12 hold only in the coincidence gauge or for the special case of connections fulfilling the spacetime symmetries D'Ambrosio22 ; Hohmann21 . In fact, choosing connections that vanish globally is not always possible in gravity, and one may end up with trivial solutions that cannot go beyond GR, regardless of the under consideration D'Ambrosio22bis .
In our analysis, we focus on the late-time evolution of the cosmic fluid, so that we can neglect radiation and consider the entire contribution due to pressureless matter. This implies and , where the subscript zero refers to quantities evaluated at the present time, and is the redshift defined as .111At the present time, and .
III cosmography
The cosmographic method allows to study the universe dynamics through kinematic quantities that do not depend on a specific background cosmology. Thus, by only assuming the homogeneity and isotropy of spacetime according to the Cosmological Principle, the late-time expansion history of the universe can be investigated in a model-independent way in order to extract information on the dark energy properties and the nature of gravity.
Before proceeding to the reconstruction of gravity, we briefly review in the following the main aspects of the cosmographic approach.
III.1 The standard cosmographic recipe
The key ingredient of cosmography is the Taylor expansion of the cosmological scale factor around the present time:
(13) |
The above expansion defines the so-called cosmographic series, whose first four terms are cosmography1
(14a) | |||
(14b) |
These quantities are named, respectively, Hubble, deceleration, jerk and snap parameters, and can be used to express cosmological distances without the need of an a priori specific model.
Hence, inserting Eq. 13 into the definition of the luminosity distance, one readily finds
(15) |
where
(16a) | ||||
(16b) | ||||
(16c) |
Futhermore, using Eq. 15, it is possible to parametrize the cosmic history starting from the relation
(17) |
which thus gives
(18) |
where
(19a) | ||||
(19b) | ||||
(19c) |
Although straightforward to implement, the standard cosmographic method presents severe limitations when dealing with high-redshift, due to the short convergence radius proper of Taylor series. One possible way to alleviate such problem is to consider rational polynomials which are able to extend the convergence of the cosmographic series towards . A relevant example in this respect is offered by Padé approximations Pade .
III.2 High-redshift cosmography with Padé approximations
One of the most reliable cosmographic methods guaranteeing a stable behaviour at high redshifts is based on Padé approximations. This is constructed starting from the Taylor series of a given function of the redshift, , where are constant coefficients. Thus, we can define the Padé approximation of as Litvinov93 ; Baker96
(20) |
Requiring that , then one determines the coefficients of the above expansion as
(21) |
The issue of choosing the degrees of rational polynomials has been addressed in the recent study Capozziello20 . In particular, it is shown that Padé approximations with polynomials of the same order in the numerator and denominator are prone to induce numerical errors and, thus, inaccurate cosmographic outcomes. Moreover, a fitting analysis to low and high-redshift observations confirmed that the most suitable Padé approximation is provided by the (2,1) polynomial, which also gives an optimal statistical performance given the low number of free coefficients involved.
In view of the aforementioned considerations, we here make use of the (2,1) Padé approximation to obtain a reliable cosmographic reconstruction of gravity. Specifically, the (2,1) approximation of the luminosity distance reads
(22) |
and the corresponding Hubble expansion rate is
(23) |
where
(24a) | ||||
(24b) |
At this point, we note that Eqs. 22 and 23 depend on a set of three cosmographic parameters, namely , whose values are found through a direct comparison with data. Accordingly, in the present analysis, we adopt the numerical results obtained in Capozziello20 by means of the Markov Chain Monte Carlo integration technique applied to the combined likelihood of Supernovae Ia Pantheon and observational Hubble data OHD . We report below the results at confidence level:
(25a) | ||||
(25b) | ||||
(25c) |
where (100 km s-1Mpc. It is worth remarking that these values have been obtained in the context of a flat universe by fixing the present matter density parameter as .
III.3 Reconstruction of the action


To reconstruct the function of the gravitational action, we consider the first Friedman equation (10). Then, in view of Eq. 12, we can convert the derivatives with respect to the non-metricity tensor in terms of derivatives of the Hubble parameter as a function of the redshift. Specifically, one finds
(26) |
where the prime denotes derivative with respect to . Using Eq. 26, we can recast Eq. 10 in the form
(27) |
Hence, assuming that the Hubble expansion is well approximated by Eq. 23, and adopting the results given in Eqs. 25a, 25b and 25c, we can solve numerically Eq. 27.
In order to do that, a boundary condition is needed. To this purpose, we recall the effective gravitational constant in gravity, i.e. Jimenez . A reasonable requirement is that coincides with Newton’s constant at present epoch, which translates into having at . Applying this in Eq. 10 leads to the following initial condition:
(28) |
Once is known, we can invert Eq. 12 by means of Eq. 23 and find . Finally, the function is obtained by plugging back into . In so doing, we find that the numerical solution suitably matches with the function
(29) |
for the following set of constant coefficients:
(30) |
We display our results in Fig. 1. We note that the test function (29) recovers pure GR for and , while the CDM model for and .
The numerical outcomes of Eq. 30 suggest (small) deviations from the standard cosmological scenario. Possible inconsistencies with CDM can be quantified by including, in the reconstruction procedure, the 1 uncertainties over the cosmographic series as given in Eqs. 25a, 25b and 25c. Such an analysis yields
(31a) | ||||
(31b) | ||||
(31c) |
We note that all the above ranges are consistent with the CDM model. Fig. 2 shows the reconstructed compared to the predictions of CDM with and . As one can see, the best-fit curves of the two scenarios are hardly distinguishable at late-times for , while they show small deviations from each other as the redshift increases. However, our results are well in agreement with the standard model within the confidence level.
IV Discussion and conclusions
Considering non-metricity as the mediator of gravitational interaction, we focused on a gravitational action containing a generic function of non-metricity scalar, which gives rise to the class of theories of gravity.
In this framework, the question of finding the able to provide the correct cosmological behaviour has been addressed in a model-independent way by means of cosmography. Relying only on the validity of the cosmological principle, such a method allows to reconstruct the gravity action without any a priori ansatz on the underlying cosmological background.
To do that, we performed a numerical reconstruction based on rational Padé approximations, which are able to reduce the convergence issues typical of the standard cosmographic technique, offering thus a reliable tool to describe cosmological observables up to high redshifts. Taking into account constraints on the cosmographic series, obtained by a direct comparison with observations in the context of a flat universe, we expressed the Hubble expansion rate as a function of the redshift and, then, we exploited the relation to reconstruct through a numerical inversion procedure.
We found that the best analytical match to our numerical outcomes is given by the function , suggesting small deviations from the CDM model as the redshift increases. However, including in the analysis the experimental uncertainties over the cosmographic parameters, our results indicate that departures from the standard cosmological model are not present at the confidence level.
It is worth noticing that the absence of significant deviations from CDM is mainly due to the large uncertainty over the jerk parameter that propagates in our numerical analysis, limiting somewhat its effective predictability at very high redshifts. Nevertheless, this problem might be alleviated by the upcoming measurements from future experiments, which could provide more stringent bounds on higher-order terms of the cosmographic series and, thus, be more sensitive in testing possible inconsistencies with standard cosmology.
Furthermore, it is necessary to bear in mind the working hypothesis of the present study. As pointed out in Sec. II, our reconstruction procedure mainly relies on Eq. 12 that holds in particular cases, as the coincidence gauge. Indeed, the same relation is not valid in other circumstances and this could lead to different results for the function D'Ambrosio22 .
Future efforts will be dedicated to compare the function here obtained with the large scale structure observations, to study the behaviour of our model at the perturbation level. In this respect, the method presented here may represent a valuable tool to explore the intrinsic dark energy properties and break the degeneracy among cosmological models.
Acknowledgements.
S.C. and R.D. acknowledge the support of Istituto Nazionale di Fisica Nucleare, Sez. di Napoli, iniziativa specifiche QGSKY and MoonLIGHT2.References
- (1) S. Perlmutter et al., Nature 391, 51 (1998); A. G. Riess et al., Astron. J. 116, 1009 (1998).
- (2) D. J. Eisenstein, W. Hu, M. Tegmark, Astrophys. J. 504, L57 (1998).
- (3) N. Aghanim et al. (Planck Collaboration), Astron. Astrophys. 641, A6 (2020).
- (4) B. S. Haridasu, V. V. Luković, R. D’Agostino, N. Vittorio, Astron. Astrophys. 600, L1 (2017).
- (5) V. Sahni, A. Starobinsky, Int. J. Mod. Phys. D 9, 373 (2000); S. M. Carroll, Living Rev. Rel. 4, 1 (2001); P. J. E. Peebles, B. Ratra, Rev. Mod. Phys. 75, 559 (2003); E. J. Copeland, M. Sami, S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006).
- (6) S. Weinberg, Rev. Mod. Phys. 61, 1 (1989); T. Padmanabhan, Phys. Rept. 380, 235 (2003).
- (7) P. J. Peebles, R. Ratra, Astrophys. J. 325, 1220 (1988); R. R. Caldwell, R. Dave, P. J. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998); R. D’Agostino, Phys. Rev. D 99, 103524 (2019); R. D’Agostino, O. Luongo, e-Print: 2112.12816 [astro-ph.CO].
- (8) S. Capozziello, R. D’Agostino, O. Luongo, Phys. Dark Univ. 20, 1 (2018); K. Boshkayev, R. D’Agostino, O. Luongo, Eur. Phys. J. C 79, 332 (2019); S. Capozziello, R. D’Agostino, R. Giambò, O. Luongo, Phys. Rev. D 99, 023532 (2019).
- (9) S. Capozziello, R. D’Agostino, O. Luongo, Int. J. Mod. Phys. D 28, 1930016 (2019); S. Nojiri and S. D. Odintsov, Phys. Rept. 505, 59 (2011); S. Capozziello, M. De Laurentis, Phys. Rept. 509, 167 (2011); K. Bamba, S. Capozziello, S. Nojiri and S. D. Odintsov, Astrophys. Space Sci. 342, 155 (2012).
- (10) A. De Felice, S. Tsujikawa, Living Rev. Rel. 13, 3 (2010); S. Nojiri, S. D. Odintsov, V. K. Oikonomou, Phys. Rept. 692, 1 (2017);
- (11) G. R. Bengochea. R. Ferraro, Phys. Rev. D 79, 124019 (2009); E. Linder, Phys. Rev. D 82, 109902 (2010); Y. F. Cai, S. Capozziello, M. De Laurentis, E. N. Saridakis, Rept. Prog. Phys. 79, 106901 (2016); R. D’Agostino, O. Luongo, Phys. Rev. D 98, 124013 (2018).
- (12) R. D’Agostino, R. C. Nunes, Phys. Rev. D 101, 103505 (2020); A. Bonilla, R. D’Agostino, R. C. Nunes, J. C. N. de Araujo, J. Cosm. Astrop. Phys. 03, 015 (2020); R. D’Agostino, R. C. Nunes, Phys. Rev. D 100, 044041 (2019).
- (13) J. Beltrán Jiménez, L. Heisenberg, T. Koivisto, Phys. Rev. D 98, 044048 (2018); E. Saridakis et al. Modified Gravity and Cosmology, Springer, Cham, Switzerland (2021).
- (14) F. Bajardi, D. Vernieri, S. Capozziello, Eur. Phys. J. Plus 135, 912 (2020); B. J. Barros, T. Barreiro, T. Koivisto, N. J. Nunes, Phys. Dark Univ. 30, 100616 (2020); A. Narawade, L. Pati, B. Mishra, S. K. Tripathy, arXiv:2203.14121.
- (15) I. Ayuso, R. Lazkoz, V. Salzano, Phys. Rev. D 103, 063505 (2021); N. Frusciante, Phys. Rev. D 103, 044021 (2021); F. K. Anagnostopoulos, S. Basilakos, E. N. Saridakis, Phys. Lett. B 822, 136634 (2021).
- (16) S. Mandal, D. Wang, P. K. Sahoo, Phys. Rev. D 102, 124029 (2020); G. Mustafa, Z. Hassan, P. H. R. S. Moraes, P. K. Sahoo, Phys. Lett. B 821, 136612 (2021). R. Solanki, A. De, P. K. Sahoo, Phys. Dark Univ. 36, 100996 (2022).
- (17) S. Weinberg, Gravitation and Cosmology, Wiley, New York (1972).
- (18) M. Visser, Gen. Rel. Grav. 37, 1541 (2005); C. Cattoen, M. Visser, Class. Quant. Grav. 24, 5985 (2007). M. Visser, Class. Quant. Grav. 32, 135007 (2015).
- (19) S. Capozziello, R. D’Agostino, O. Luongo, Mon. Not. Roy. Astron. Soc. 476, 3924 (2018); P. K. S. Dunsby, O. Luongo, Int. J. Geom. Meth. Mod. Phys. 13, 1630002 (2016); A. Aviles, C. Gruber, O. Luongo, H. Quevedo, Phys. Rev. D 86, 123516 (2012); S. Capozziello, R. D’Agostino, O. Luongo, arXiv: 2202.03300 [astro-ph.CO].
- (20) S. Capozziello, R. D’Agostino, O. Luongo, Gen. Rel. Grav. 49, 141 (2017); S. Capozziello, R. D’Agostino, O. Luongo, J. Cosm. Astrop. Phys. 05, 008 (2018); S. Capozziello, R. D’Agostino, O. Luongo, Gen. Rel. Grav. 51, 2 (2019).
- (21) L. Järv, M. Rünkla, M. Saal, O. Vilson, Phys. Rev. D 97, 124025 (2018).
- (22) J. Beltrán Jiménez, L. Heisenberg, T. S. Koivisto, Universe 5, 173 (2019).
- (23) J. Beltrán Jiménez, L. Heisenberg, T. Koivisto, S. Pekar, Phys. Rev. D 101, 103507 (2020).
- (24) F. D’Ambrosio, L. Heisenberg, S. Kuhn, Class. Quant. Grav. 39, 025013 (2022).
- (25) M. Hohmann, Phys. Rev. D 104, 124077 (2021).
- (26) F. D’Ambrosio, S. D. B. Fell, L. Heisenberg, S. Kuhn, Phys. Rev. D 105, 024042 (2022).
- (27) C. Gruber, O. Luongo, Phys. Rev. D 89, 103506 (2014); H. Wei, X. P. Yan, Y. N. Zhou, J. Cosm. Astrop. Phys. 1401, 045 (2014); A. Aviles, A. Bravetti, S. Capozziello, O. Luongo, Phys. Rev. D 90, 043531 (2014).
- (28) G. Litvinov, Appl. Russ. J. Math. Phys. 1, 313 (1993).
- (29) G. A. Baker Jr., P. Graves-Morris, Padé Approximants, Cambridge University Press (1996).
- (30) S. Capozziello, R. D’Agostino, O. Luongo, Mon. Not. Roy. Astron. Soc. 494, 2576 (2020).
- (31) D. M. Scolnic, Astrophys. J. 859, 101 (2018); A. G. Riess et al., Astrophys. J. 853, 126 (2018).
- (32) R. Jimenez, A. Loeb, Astrophys. J. 573, 37 (2002).