Constraining über gravity with recent observations and studying the H0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_0$$\end{document} problem

This paper studies both Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM and CDM models under the über gravity theory, named üΛ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM and üCDM respectively. We report bounds over their parameter phase-space using several cosmological data, in particular, the recent Pantheon+ sample. Based on the joint analysis, the best fit value of the über characteristic parameter is z⊕=0.046-0.032+0.047\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z_\oplus = 0.046^{+0.047}_{-0.032}$$\end{document} and z⊕=1.382-0.021+0.020\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$z_\oplus = 1.382^{+0.020}_{-0.021}$$\end{document} at 68% confidence level for üΛ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM and üCDM respectively. Although über gravity can successfully mimics the cosmological constant, we find that the H0(z)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {H}0(z)$$\end{document} diagnostic suggests the H0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H_0$$\end{document} tension is not alleviated. Finally, both models are statistically compared with Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM through the Akaike and Bayesian information criteria. Both über gravity models and Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM are equally preferred for most of the single samples, in particular, üΛ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda $$\end{document}CDM is not rejected by the CMB data. However, there is strong evidence against them for the joint analysis.


I. INTRODUCTION
The Λ Cold Dark Matter (ΛCDM) model is one of the most successful theories that explain with great precision the beginning and the evolution of our Universe at its current stages.However, among its many conundrums, the dark matter and energy components are the most intriguing due to the extremely weak interactions (or possible negligible) with the standard fields, being only the gravitational field the sole probe of its existence.The dark energy (DE) is related to the current Universe acceleration first observed by the supernova teams lead by Riess, Perlmutter among others [1,2] and later sustained by the Planck Satellite [3].Clearly an acceleration of the space-time is not expected, thus it requires an unnatural dynamics: the addition of an extra component into the Friedmann equations which is able to create an accelerated expansion of the Universe and produce a late de Sitter scale factor.Without doubts, the cosmological constant is the cheapest way to tackle this problem, being an essential part of the ΛCDM paradigm, having also an exquisite capability to simulate the late Universe [4].However, it is important to emphasize that the cos-mological constant is afflicted with several problems like those related to the vacuum energy density and the coincidence problem [5,6].Is in this vein that the community is exploring diverse alternatives to the cosmological constant, being the problem not yet settled (see [7] for a compilation).
Additionally, cosmologists concentrate most of their efforts in tackling the case of the current expansion rate of the Universe, the so-called H 0 tension.The reduction of uncertainties has led to significant discrepancies between the Planck [3] and SH0ES [8] collaborations, whose disagreements are approximately at ∼ 5.0σ.Whether these discrepancies are due to systematic errors or mistaken theoretical assumptions, the H 0 tension could represent a window into new physics.Although modifying the standard cosmological model, in order to address the tension, without compromising its success in other areas has proven to be a very difficult task [9], the model-dependence of early H 0 measurements may suggest that it is worth exploring models beyond ΛCDM.The proposals include early dark energy [10][11][12], extrarelativistic degrees of freedom [13][14][15], new early dark energy [16,17], diffusion models [18,19], holographic dark energy [20][21][22], to mention some of them.
In the present study we focus on two proposals: the üΛCDM and üCDM models, based on the übergravitational model [23], which understands the H 0 tension as a transition in the context of a gravitational theory [24].Über-gravity sets the idea of taking an ensemble average over all possible f (R) theories [25], a process in-spired by statistical mechanics.Although the criteria for assigning probabilities to each theory are not clear, the general properties of üΛCDM and üCDM can be studied and it shows promise in alleviating the H 0 tension.Indeed, the model contains two branches where the Hubble parameter evolves in different ways, with a threshold region determined by the free parameter z ⊕ , which is the redshift of transition to the über gravity.In this sense, it is also possible to determine the deceleration and jerk parameters in order to elucidate the characteristics of the component responsible for the acceleration.As stated, üΛCDM contains an extra parameter in comparison with üCDM, which contains the same number as ΛCDM.Therefore, in the case of üΛCDM, besides introducing a statistical penalization (e.g.Akaike (AIC) or Bayesian information criteria (BIC) respectively) due to the extra parameter, it also keep the unsolved problem associated to the cosmological constant.
Additionally, we constrain the free parameters through a Bayesian analysis, using recent observations of Observational Hubble Parameters (CMB), Type Ia Supernovae (SnIa), Baryon Acoustic Oscillations (BAO), HII Galaxies (HIIG), Cosmic Microwave Background radiation (CMB) and a joint analysis.Moreover, we implement a H 0 diagnostic technique to study a possible alleviation to the H 0 tension using the über-gravity paradigm.
The paper is organized as follows: in Sec.II we review the mathematical background of üΛCDM and üCDM models, emphasizing its main properties and defining the parameters to be analyzed.In Sec.III the data and methodology are described, and the results are shown in Sec.IV.Finally, in V we summarize the conclusions and give some outlooks on the subject.We henceforth use units in which = k B = c = 1.

II. MATHEMATICAL BACKGROUND
The idea of taking an ensemble average over gravity theories is condensed in the following Lagrangian where the index i corresponds to the i-th theory considered through the Lagrangian associated L, being N the dimension of the gravitational models space over the manifold M. The free parameter β assigns probabilities to each model.For greater rigor, the ensemble average is taken over all analytic models of gravity, described by the f (R) family, where the R stands for the Ricci scalar.However, in [24] is shown that the Lagrangian, after a change of basis, can be written as where R ≡ R/R 0 contains a new cosmological parameter R 0 , n is an integer, and Λ is the cosmological constant.Next, we will show two models associated with üΛCDM and üCDM.

A. üΛCDM model
In background cosmology it is possible to assume homogenity and isotropy in the line element through the Friedmann-Lemaitre-Robertson-Walker (FLRW) equations as ds 2 = −dt 2 + a(t) 2 [dr 2 + r 2 dΩ 2 ], where dΩ 2 = dθ 2 + sin 2 θdϕ 2 and a(t) is the scale factor.Thus, the Friedmann equation in this scenario can be written as [23] where E(z) ≡ H(z)/H 0 , H(z) is the Hubble parameter, z is the redshift (i.e.H 0 = H(z = 0)), Ω 0m is the matter density parameter, Ω 0r is the radiation density parameter, both at z = 0 and R0 ≡ R 0 /6H 2 0 .In this case, Ω 0Λ plays a role at z > z ⊕ , being z ⊕ the region of transition to über gravity.
Using the continuity for E(z) and E (z), where the prime denotes a derivative with respect to z, R0 takes the form One of the magnitudes of interest for the present study is the deceleration parameter, written in terms of the redshift, which is given by Moreover, the jerk parameter is given by the expression ) which behaves like the standard ΛCDM model in the region z > z ⊕ , while for z < z ⊕ the über gravity dominates its evolution.
Finally, the transition redshift z T , which determines the beginning of the accelerated expansion of the universe (q(z T ) = 0), was calculated for both models.For z T < z ⊕ the transition redshift takes the form whereas for z T > z ⊕ the value comes from solving which is obtained numerically.The transition in this model occurs for z T > z ⊕ , so in this study only this last expression is taken into account.

B. üCDM model
In this case, we propose that the Universe acceleration is only driven by über contributions, thus we consider Ω Λ = 0, reducing the free parameters to the same number as in the standard model and assuming FLRW cosmology.Therefore, we have For the deceleration parameter we have ) and the jerk parameter reads being the first region z > z ⊕ (for q and j) dominated by matter while the later one is controlled by the über effects that drive the acceleration.Finally, in üCDM the transition is not allowed in the region z T > z ⊕ , giving, as a result, the following analytic expression, valid for Note that z T shows a dependence only on z ⊕ and Ω 0r .

III. DATA AND METHODOLOGY
Both üΛCDM and üCDM cosmologies are confronted using CC, HIIG, SnIa, BAO, CMB and joint data through a Bayesian Markov Chain Monte Carlo (MCMC) analysis.To bound their free parameter phase-space, (h, Ω 0Λ , z ⊕ ) for üΛCDM and (h, z ⊕ ) for üCDM, we use the emcee package [26] under Python language.We establish a configuration to achieve the convergence of the chains using the autocorrelation function, and generate a set of 3000 chains with 250 steps.Additionally, we use a Gaussian prior over h as h = 0.7403 ± 0.0142 [8] and h = 0.6766 ± 0.0042 [3] as a consistency probe, flat priors over Ω 0Λ and z ⊕ in the region [0,1] and [0,2] respectively.Thus, the χ 2 -function is given by where each term corresponds to the χ 2 function per sample.

A. Cosmic chronometers
We use a sample of 31 measurements of the Hubble parameter using differential age method [27] (see also [28]).Due these points are cosmological model independent, they are useful to test alternative cosmologies to ΛCDM.The χ 2 function can be built as where H th (z i ) is the theoretical Hubble parameter using Eq. ( 3), and H obs (z i ) ± σ i obs is the observational counterpart with its uncertainty at the redshift z i .

B. Type Ia Supernovae (Pantheon+)
Recently, a sample of 1701 measurements of the luminosity modulus coming from SNIa, namely Pantheon+ sample, is reported by [29,30].This represents the largest sample and covers a region 0.001 < z < 2.26.Considering that this sample is extracted from 1550 distinct SNIa, we build the χ 2 function as where and ∆μ is the vector of the difference between the theoretical distance modulus and the observed one, ∆1 = (1, 1, . . ., 1) T , Cov P is the covariance matrix formed by adding the systematic and statistic uncertainties.The transpose of the vectors are denoted with the super-index T on the previous expressions.
The theoretical counterpart of the distance modulus is estimated by where M is a nuisance parameter which has been marginalized by Eq. ( 16).The luminosity distance, denoted as d L (z), is computed through where c is the speed of light.

C. Baryon Acoustic Oscillations
The Baryon Acoustic Oscillations signature is an standard ruler useful to constrain cosmological model parameters and is the result of the interactions between baryons and photons in the recombination era.We use 6 correlated points obtained by [31][32][33] and collected by [34].We confront them with the cosmological models by building the χ 2 function as where X T is the residual vector between observational measurements and theoretical values of the ratio d A (z * )/D V (z), where d A (z d ) is the comoving angular diameter distance at the the photon decoupling epoch (z * ) and the dilation scale is given by [35] where c is again the speed of light.For this work we use z * = 1089.80± 0.21 [3].

D. HII Galaxies
A sample of 181 measurements coming from Hydrogen II galaxies (HIIG), with their luminosity dominated by young massive burst of star formation, is reported by [36,37].This sample which covers a region 0.01 < z < 2.6 and is useful to establish bounds over cosmological parameters due the correlation between the measured luminosity L of the galaxies and the inferred velocity dispersion σ of their ionized gas [38][39][40][41].The χ 2 -function is built as where i is the observational uncertainty measured at z i having 68% of confidence level.Additionally, the observational distance modulus (µ obs ) is expressed Here, α and β are the intercept and slope of the L-σ relation and f is the measured flux.The theoretical estimate is given as where d L is the luminosity distance measured in Mpc (see Eq. ( 19)).

E. Cosmic Microwave Background Radiation
The CMB temperature anisotropies are useful measurements to establish constraints over cosmological parameters.A way to use them without performing a full perturbative analysis is to compress the full information into some parameters.Authors in [42] compress the CMB information from Planck 2018 Temperature Power Spectrum (TT), for high TE multipoles, polarization spectra EE modes + lowE data [3] in the acoustic scale l A which characterizes the CMB temperature in the transverse direction, the shift parameter R which influences the CMB temperature along the line-of-sight direction, and the quantity Ω b0 h 2 where Ω b0 is the density of baryons today (z = 0).Thus, the figure-of-merit is built as where V CMB is which is the covariance matrix for V CMB .The theoretical counterparts are estimated by and where z * is the redshift at the photon decoupling phase that takes the form where and D A is the angular diameter distance for a flat geometry that reads as while r s is the comoving sound horizon given by being Ω b and Ω γ the baryons and photons density parameters respectively [43].

IV. RESULTS
We use the über gravity in the context of üΛCDM and üCDM models to understand the nature of the current Universe acceleration and help us to interpret the origin of the H 0 tension.Due to the consistency between the results obtained by using H 0 Gaussian priors from Planck and SH0ES, hereafter we centered our discussion in the case of H 0 with priors coming from low-redshift (see Table I).In Figs. 1 and 2 we present the 2D parameter likelihood contours at 68% (1σ) and 99.7% (3σ) confidence level (CL) respectively for both models.Moreover Table I shows the mean values of the parameters and their uncertainties at 1σ.
Additionally to this, we statistically compare these results with ΛCDM model, applying the corrected Akaike information criterion (AICc) [44] and the Bayesian information criterion (BIC) [45].Both criteria penalize according to the size of the data sample (N ) and the number of degrees of freedom (k) defined as AICc = χ 2 min + 2k +2(k 2 +k)/(N −k −1) and BIC = χ 2 min +k log(N ), respectively, where χ 2 min is the minimum of χ 2 .Therefore, a model with a lower value of AICc (BIC) is preferred by the data.
For the ∆AICc we have the following conditions: if |∆AICc| < 4 both models are statistically equivalent, if 4 < AICc < 10 the data still support the given model (über gravity) but less than the preferred one (concordance model) and finally if AICc > 10 indicates that the data do not support the given model.In addition, for ∆BIC the interpretation is as follows: if ∆BIC < 2 there is no evidence against the new model which in this case is über gravity, if 2 < ∆BIC < 6 there is modest evidence against the new model, and finally if ∆BIC > 10, gives the strongest evidence against it.
We start by showing the AICc and BIC for the ΛCDM model: for AICc values we have 22.07 (CC), 440.72 (HIIG), 2002.41(SnIa), 12.82 (BAO), 23.80 (CMB) and 2493.25 (Joint); as for BIC it has 24.51 (CC), 447.05 (HIIG), 2013.28 (SnIa), 6.04 (BAO), 21.56 (CMB) and 2504.37 (Joint).Table I shows both ∆AICc and ∆BIC by considering ΛCDM as the reference scenario.According to AICc, we find that both üΛCDM and ΛCDM are preferred equally by CC, HIIG and CMB while the üΛCDM is not supported by SnIa and BAO, but there is a still support for the über gravity when the combined data are considered.Regarding BIC, we find the strongest evidence against üΛCDM for SnIa and the combined data and a modest evidence against the uber gravity for HIIG.
As for üCDM, we observe that üCDM and ΛCDM are equally preferred by CC, HIIG and BAO while the strongest evidence against üCDM is provided by SnIa and Joint analysis, and a strong evidence against über gravity is given by CMB.It is worth to mention that ∆AICc and ∆BIC are reduced to ∆χ 2 because both models have same number of free parameters.
On the other hand, notice that the üΛCDM scenario produces a late acceleration, constraining the free parameter of über gravity according to the joint analysis as z ⊕ 0.046, which also coincides with the epoch of transition to über gravity.The parameter Ω 0Λ 0.689 marks a cosmological constant domination, however the unsolved questions about its characteristics remains.The redshift transition is compatible with ΛCDM model which is at 0.646, as estimated from the joint analysis.The transition from GR to über gravity is perfectly observed through the jerk factor shown in Fig. 3, where a discontinuity in the j(z) function is observed.According to the H 0 diagnostic presented in Fig. 4, the tension in H 0 is not alleviated because it keeps the trend towards values consistent with supernovae results.
Additionally, we explore the üCDM model which is not studied in the literature.In this case, the über component z ⊕ acts like a cosmological constant, having the same free parameters as ΛCDM.For the joint analysis we conclude that the value of the über parameter us z ⊕ 1.382 and the transition to an accelerated Universe happens earlier in its evolution z T 0.809.From Fig. 3, in particular with q(z) and j(z), we observe that the transition to über gravity domination is more extreme in comparison with the üΛCDM model.Additionally to this, in this figure we also observe a trend to values that are consistent with supernova results instead of those obtained from Planck, concluding that the H 0 tension persist despite the über gravity presence.

V. CONCLUSIONS AND OUTLOOKS
This paper presents an exhaustive revision of the über gravity for the üΛCDM and üCDM models, obtaining stringent constrictions with different cosmological data samples, in particular the recent Pantheon+ sample.First, we focus on our results obtained assuming a prior on h provided by SHOES because these are consistent with those we obtain using a prior from Planck (see Table I).We statistically compare both models with ΛCDM by using AICc and BIC.We find that the uber cosmologies and ΛCDM are equally preferred for CC, HIIG, and BAO.When CMB data are used, üΛCDM    Here we present the behavior of H0(z) under the constraints of the four data-samples and the joint analysis.The continuous red line represents the results for ΛCDM cosmology assuming h = 0.6766 and Ω0m = 0.3111 according to [3].
does not present evidence against but üCDM presents a strong evidence against.However, for both über cosmologies, üΛCDM and üCDM, the joint analysis shows the strongest evidence against them.
First, the üΛCDM contains one extra free parameter in comparison with the standard paradigm, maintaining the same open questions related to the understanding of the cosmological constant.Our results for the joint analysis point out that Ω 0Λ 0.689, consistent with the standard cosmological model, but with a subdominant value for the über parameter z ⊕ 0.046.It is important to remind that the über parameter z ⊕ is the point where the über gravity starts to dominate over the standard GR, thus, the presence of the ΛCDM delays the apparition of the über gravity as observed in Figs. 3 for the H(z), q(z) and j(z) parameters.The transition redshift is also in consistency with the standard cosmological model (see Table I).Regarding H0(z) diagnostic, Fig. 4 reveals that, according to the joint analysis, the H 0 value is more consistent with the supernova data than with the Planck ones at z = 0, unable to reduce the tension under this scenario.This could be the result of ΛCDM being behind the dynamics while über gravity start its domination at z ⊕ 0.046.Additionally to this, the new sample of Pan-theon+ could generate a tendency to a value greater than the one predicted by Planck.Furthermore, we explore the üCDM, where the über parameter plays the role of the cosmological constant, thus having the same parameters as in the standard cosmology.In this case, über parameter generates the late acceleration of the Universe and its presence start earlier, specifically at z ⊕ 1.382 according to the joint analysis presented in Table I.The transition also happens earlier than in the standard model, at z T 0.809 (see Table I).The evolution of H(z), q(z) and j(z) is presented in Fig. 3 where we can see that the behavior is more abrupt, mainly in the q(z) and j(z) behavior.However, über gravity eventually mimics the cosmological constant because j = 1 when z = −1.The H0(z) diagnostic plot presented in Fig. 4 shows that at z = 0 the value for H 0 does not coincide with the Planck result for ΛCDM, being in concordance with supernova results, which is also a sign that the über parameter acting like a cosmological constant is inadequate to alleviate the tension in H 0 .
Finally, it is worth to point out the difference in the values of z ⊕ for üCDM and üΛCDM.The discrepancy is caused mainly because in the first one the über parameter acts like a cosmological constant and not only as a transition to über gravity, which is the case for üΛCDM.
In summary, we observe that the über gravity is an alternative to study the late Universe acceleration by mimicking the cosmological constant behavior through the über parameter z ⊕ .However, the mystery of the Hubble tension remains unsolvable and gives values that are compatibles with supernova results.Nevertheless, we consider that a deeper exploration of the über Lagrangian for other values of n is necessary in order to elucidate any extra dynamics that could help us understand the reason for the Hubble tension.Such exploration will be presented elsewhere.
the superscripts th refers to the theoretical estimates, and Cov −1 CMB represents the inverse of Cov CMB = 10 −8

FIG. 4 .
FIG.4.The H0(z) diagnostic for über gravity, showing üCDM and üΛCDM models.Here we present the behavior of H0(z) under the constraints of the four data-samples and the joint analysis.The continuous red line represents the results for ΛCDM cosmology assuming h = 0.6766 and Ω0m = 0.3111 according to[3].