f(T)𝑓𝑇f(T)italic_f ( italic_T ) gravity after DESI Baryon Acoustic Oscillation and DES Supernovae 2024 data

Celia Escamilla-Rivera |||| celia.escamilla@nucleares.unam.mx Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Circuito Exterior C.U., A.P. 70-543, México D.F. 04510, México.    Rodrigo Sandoval-Orozco rodrigo.sandoval@correo.nucleares.unam.mx Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Circuito Exterior C.U., A.P. 70-543, México D.F. 04510, México.
Abstract

In this work, we investigate new constraints on f(T)𝑓𝑇f(T)italic_f ( italic_T ) gravity using the recent Baryon Acoustic Oscillation (BAO) data released by the Dark Energy Spectroscopic Instrument (DESI) and the Type Ia supernovae (SNIa) catalog from the full 5-years of the Dark Energy Survey Supernova Program (DES-SN5YR). The f(T)𝑓𝑇f(T)italic_f ( italic_T ) cosmological models considered are characterised by power law late-time accelerated expansion. Our results show that the combination DESI BAO +rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT CMB Planck suggests a Bayesian preference for late-time f(T)𝑓𝑇f(T)italic_f ( italic_T ) cosmological models over ΛΛ\Lambdaroman_ΛCDM, obtaining a value of H0=68.33.5+3.0subscript𝐻0subscriptsuperscript68.33.03.5H_{0}=68.3^{+3.0}_{-3.5}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 68.3 start_POSTSUPERSCRIPT + 3.0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 3.5 end_POSTSUBSCRIPT[km/s/Mpc] in agreement with SH0ES collaboration, however, due to a bigger uncertainty.

The Hubble tension has reached a statistical significance of similar-to\sim5σ𝜎\sigmaitalic_σ, strongly proven a mismatch between the cosmic late-time expansion rate H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT measured through the local distance ladder method using Type Ia supernovae measurements Riess et al. (2022); Breuval et al. (2024) H0=73±1subscript𝐻0plus-or-minus731H_{0}=73\pm 1italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 73 ± 1 [km/s/Mpc], and the inferred H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT value from observations of the Cosmic Microwave Background (CMB) radiation, H0=67.4±0.5subscript𝐻0plus-or-minus67.40.5H_{0}=67.4\pm 0.5italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 67.4 ± 0.5 [km/s/Mpc] Aghanim et al. (2020). While cautious studies on a possible systematic origin of this mismatch have been performed by the SH0ES collaboration Riess et al. (2022), there is no signal that this could be the reason for this H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT tension issue. This result has brought interesting opportunities to change the view on how the standard cosmological model is designed, allowing us to introduce a path beyond the ΛΛ\Lambdaroman_Λ-Cold Dark Matter(CDM) model.

Current BAO measurements released by the Dark Energy Spectroscopic Instrument (DESI) Adame et al. (2024a, b) seem to point towards new physics in the dark energy cosmological scheme Collaboration (2024). Additionally, the Dark Energy Survey 5-year SNIa (DES-SN5YR) release Abbott et al. (2024), suggest a 2σ2𝜎2\sigma2 italic_σ preference for a time-varying dark energy. More discussions on this aspect have been presented in Yin (2024); Cortês and Liddle (2024). As part of the efforts to find well-constrained proposals with these releases and also, that can address the cosmological tensions, some studies have been developing in these short period, e.g. new constraints on axion-early dark energy model Qu et al. (2024) which are not tighten even after the inclusion of DESI BAO data, interacting dark energy models Wang (2024); Giarè et al. (2024) which shows a value of H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT in good agreement with SH0ES collaboration, for quintessence scalar field model Berghaus et al. (2024); Tada and Terada (2024) showing a preference within 2-4% for a kinetic scalar field energy, for dark energy models inspired in thermodynamics and parametrised equation-of-state in Taylor expansions Carloni et al. (2024); Wang and Piao (2024), the first one finding that log-corrected parameterisation could alleviate the H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT tension, and finally Gaussian reconstructions on quintom modified cosmology Yang et al. (2024). All these studies aim to increase the value of H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT inferred. On one hand, the main challenge in the early CMB measurements seems to be settled in computing the angular scale of the CMB acoustic peaks  Aghanim et al. (2020). Therefore, increasing the value of H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT without modifying the acoustic scale requires a different post-recombination epoch Knox and Millea (2020). On the other hand, late-time cosmic proposals require new physics that can change cosmic distances to compensate for the higher value of H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, taking into account the preservation of the CMB history.

Within these efforts, extended theories of gravity have been proposed as a good description of a fundamental theory of gravity that allows addressing theoretical and observational issues with viable solutions in the observed mismatch Bahamonde et al. (2021); Aguilar et al. (2024); Briffa et al. (2023); Sandoval-Orozco et al. (2024); Nunes et al. (2018); Kumar et al. (2023); Nunes (2018); dos Santos et al. (2022); Hashim et al. (2021). To formulate an extension to General Relativity, we consider a construction through the metric-affine gravity Beltrán Jiménez et al. (2019), where teleparallel gravity (TG) has a curvature-free connection Bahamonde et al. (2021); Krssak et al. (2019) with a scenario that include a teleparallel equivalent of general relativity (TEGR). This theory has described a set of field equations which are dynamically equivalent to the GR ones. Within this scheme, f(T)𝑓𝑇f(T)italic_f ( italic_T ) gravity emerges as a generalisation of the TEGR Lagrangian with a function of the torsion T𝑇Titalic_T as f(T)=T+(T)𝑓𝑇𝑇𝑇f(T)=-T+\mathcal{F}(T)italic_f ( italic_T ) = - italic_T + caligraphic_F ( italic_T ).

In this work, we show that the new DESI BAO plus rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT CMB Planck data release gives a slight Bayesian preference for extended f(T)𝑓𝑇f(T)italic_f ( italic_T ) cosmologies, despite the fact that the parameters are within 2σ𝜎\sigmaitalic_σ confidence level (C.L) from ΛΛ\Lambdaroman_ΛCDM. In addition, it is noted that high/low-z observations could be better explained in these models in comparison to ΛΛ\Lambdaroman_ΛCDM and with first principle reasons. In such a scheme, we also consider baseline with DES-SN5YR, which gives a lower value of H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT for this kind of supernovae catalog.

Parameter DESI+BBN DESI+CMB111CMB Distance priors. DESI+CC DESI+CMB +CC+PN+superscriptPN\text{PN}^{+}PN start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT DESI+CC +PN+ DESI+CC +SNYR5
H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [km/s/Mpc] 68.81.2+1.3subscriptsuperscript68.81.31.268.8^{+1.3}_{-1.2}68.8 start_POSTSUPERSCRIPT + 1.3 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 1.2 end_POSTSUBSCRIPT 68.440.83+0.87subscriptsuperscript68.440.870.8368.44^{+0.87}_{-0.83}68.44 start_POSTSUPERSCRIPT + 0.87 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.83 end_POSTSUBSCRIPT 70.66.7+6.8subscriptsuperscript70.66.86.770.6^{+6.8}_{-6.7}70.6 start_POSTSUPERSCRIPT + 6.8 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 6.7 end_POSTSUBSCRIPT 68.800.84+0.81subscriptsuperscript68.800.810.8468.80^{+0.81}_{-0.84}68.80 start_POSTSUPERSCRIPT + 0.81 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.84 end_POSTSUBSCRIPT 72.82.0+2.1subscriptsuperscript72.82.12.072.8^{+2.1}_{-2.0}72.8 start_POSTSUPERSCRIPT + 2.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 2.0 end_POSTSUBSCRIPT 66.36.1+6.2subscriptsuperscript66.36.26.166.3^{+6.2}_{-6.1}66.3 start_POSTSUPERSCRIPT + 6.2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 6.1 end_POSTSUBSCRIPT
ΩcdmsubscriptΩcdm\Omega_{\mathrm{cdm}}roman_Ω start_POSTSUBSCRIPT roman_cdm end_POSTSUBSCRIPT 0.2390.025+0.027subscriptsuperscript0.2390.0270.0250.239^{+0.027}_{-0.025}0.239 start_POSTSUPERSCRIPT + 0.027 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.025 end_POSTSUBSCRIPT 0.23820.0097+0.0095subscriptsuperscript0.23820.00950.00970.2382^{+0.0095}_{-0.0097}0.2382 start_POSTSUPERSCRIPT + 0.0095 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0097 end_POSTSUBSCRIPT 0.2380.030+0.031subscriptsuperscript0.2380.0310.0300.238^{+0.031}_{-0.030}0.238 start_POSTSUPERSCRIPT + 0.031 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.030 end_POSTSUBSCRIPT 0.23490.0089+0.0085subscriptsuperscript0.23490.00850.00890.2349^{+0.0085}_{-0.0089}0.2349 start_POSTSUPERSCRIPT + 0.0085 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0089 end_POSTSUBSCRIPT 0.2540.022+0.023subscriptsuperscript0.2540.0230.0220.254^{+0.023}_{-0.022}0.254 start_POSTSUPERSCRIPT + 0.023 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.022 end_POSTSUBSCRIPT 0.3090.022+0.023subscriptsuperscript0.3090.0230.0220.309^{+0.023}_{-0.022}0.309 start_POSTSUPERSCRIPT + 0.023 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.022 end_POSTSUBSCRIPT
wbsubscript𝑤𝑏w_{b}italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT 0.02218±0.00077plus-or-minus0.022180.000770.02218\pm 0.000770.02218 ± 0.00077 0.022590.00033+0.00034subscriptsuperscript0.022590.000340.000330.02259^{+0.00034}_{-0.00033}0.02259 start_POSTSUPERSCRIPT + 0.00034 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00033 end_POSTSUBSCRIPT 0.02500.0099+0.011subscriptsuperscript0.02500.0110.00990.0250^{+0.011}_{-0.0099}0.0250 start_POSTSUPERSCRIPT + 0.011 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0099 end_POSTSUBSCRIPT 0.022710.00034+0.00035subscriptsuperscript0.022710.000350.000340.02271^{+0.00035}_{-0.00034}0.02271 start_POSTSUPERSCRIPT + 0.00035 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00034 end_POSTSUBSCRIPT 0.02820.0036+0.0038subscriptsuperscript0.02820.00380.00360.0282^{+0.0038}_{-0.0036}0.0282 start_POSTSUPERSCRIPT + 0.0038 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0036 end_POSTSUBSCRIPT 0.01990.0079+0.0081subscriptsuperscript0.01990.00810.00790.0199^{+0.0081}_{-0.0079}0.0199 start_POSTSUPERSCRIPT + 0.0081 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0079 end_POSTSUBSCRIPT
ΩmsubscriptΩm\Omega_{\mathrm{m}}roman_Ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT 0.2860.025+0.027subscriptsuperscript0.2860.0270.0250.286^{+0.027}_{-0.025}0.286 start_POSTSUPERSCRIPT + 0.027 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.025 end_POSTSUBSCRIPT 0.2860.011+0.011subscriptsuperscript0.2860.0110.0110.286^{+0.011}_{-0.011}0.286 start_POSTSUPERSCRIPT + 0.011 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.011 end_POSTSUBSCRIPT 0.2880.025+0.027subscriptsuperscript0.2880.0270.0250.288^{+0.027}_{-0.025}0.288 start_POSTSUPERSCRIPT + 0.027 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.025 end_POSTSUBSCRIPT 0.28290.0099+0.0095subscriptsuperscript0.28290.00950.00990.2829^{+0.0095}_{-0.0099}0.2829 start_POSTSUPERSCRIPT + 0.0095 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0099 end_POSTSUBSCRIPT 0.3080.021+0.022subscriptsuperscript0.3080.0220.0210.308^{+0.022}_{-0.021}0.308 start_POSTSUPERSCRIPT + 0.022 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.021 end_POSTSUBSCRIPT 0.3530.019+0.020subscriptsuperscript0.3530.0200.0190.353^{+0.020}_{-0.019}0.353 start_POSTSUPERSCRIPT + 0.020 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.019 end_POSTSUBSCRIPT
rdsubscript𝑟dr_{\mathrm{d}}italic_r start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT [Mpc] 149.13.4+3.5subscriptsuperscript149.13.53.4149.1^{+3.5}_{-3.4}149.1 start_POSTSUPERSCRIPT + 3.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 3.4 end_POSTSUBSCRIPT 149.060.71+0.73subscriptsuperscript149.060.730.71149.06^{+0.73}_{-0.71}149.06 start_POSTSUPERSCRIPT + 0.73 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.71 end_POSTSUBSCRIPT 146.0±14.0plus-or-minus146.014.0146.0\pm 14.0146.0 ± 14.0 149.03±0.73plus-or-minus149.030.73149.03\pm 0.73149.03 ± 0.73 148.64.7+4.6subscriptsuperscript148.64.64.7148.6^{+4.6}_{-4.7}148.6 start_POSTSUPERSCRIPT + 4.6 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 4.7 end_POSTSUBSCRIPT 14714+15subscriptsuperscript1471514147^{+15}_{-14}147 start_POSTSUPERSCRIPT + 15 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 14 end_POSTSUBSCRIPT
M𝑀Mitalic_M -- -- -- 19.3910.033+0.027subscriptsuperscript19.3910.0270.033-19.391^{+0.027}_{-0.033}- 19.391 start_POSTSUPERSCRIPT + 0.027 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.033 end_POSTSUBSCRIPT 19.2650.059+0.056subscriptsuperscript19.2650.0560.059-19.265^{+0.056}_{-0.059}- 19.265 start_POSTSUPERSCRIPT + 0.056 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.059 end_POSTSUBSCRIPT 19.530.21+0.19subscriptsuperscript19.530.190.21-19.53^{+0.19}_{-0.21}- 19.53 start_POSTSUPERSCRIPT + 0.19 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.21 end_POSTSUBSCRIPT
χmin2subscriptsuperscript𝜒2min\chi^{2}_{\text{min}}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT min end_POSTSUBSCRIPT 8.928.928.928.92 12.4712.4712.4712.47 23.9923.9923.9923.99 1646.541646.541646.541646.54 1615.321615.321615.321615.32 5933.335933.335933.335933.33
Table 1: Constraints at 2σ𝜎\sigmaitalic_σ C.L for the ΛΛ\Lambdaroman_ΛCDM model. For all baselines, we provide results with and without BBN constraints. Also, we include the constraints for two SNIa baselines.
Parameter DESI+BBN DESI+CMB222CMB Distance priors DESI+CC DESI+CMB +CC+PN+superscriptPN\text{PN}^{+}PN start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT DESI+CC +PN+ DESI+CC +SNYR5
H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [km/s/Mpc] 68.5±6.0plus-or-minus68.56.068.5\pm 6.068.5 ± 6.0 69.8±2.5plus-or-minus69.82.569.8\pm 2.569.8 ± 2.5 70.16.9+7.2subscriptsuperscript70.17.26.970.1^{+7.2}_{-6.9}70.1 start_POSTSUPERSCRIPT + 7.2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 6.9 end_POSTSUBSCRIPT 69.11.3+1.5subscriptsuperscript69.11.51.369.1^{+1.5}_{-1.3}69.1 start_POSTSUPERSCRIPT + 1.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 1.3 end_POSTSUBSCRIPT 72.4±2.2plus-or-minus72.42.272.4\pm 2.272.4 ± 2.2 66.56.4+6.1subscriptsuperscript66.56.16.466.5^{+6.1}_{-6.4}66.5 start_POSTSUPERSCRIPT + 6.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 6.4 end_POSTSUBSCRIPT
ΩcdmsubscriptΩcdm\Omega_{\mathrm{cdm}}roman_Ω start_POSTSUBSCRIPT roman_cdm end_POSTSUBSCRIPT 0.2490.034+0.031subscriptsuperscript0.2490.0310.0340.249^{+0.031}_{-0.034}0.249 start_POSTSUPERSCRIPT + 0.031 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.034 end_POSTSUBSCRIPT 0.246±0.017plus-or-minus0.2460.0170.246\pm 0.0170.246 ± 0.017 0.2430.046+0.042subscriptsuperscript0.2430.0420.0460.243^{+0.042}_{-0.046}0.243 start_POSTSUPERSCRIPT + 0.042 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.046 end_POSTSUBSCRIPT 0.2480.012+0.013subscriptsuperscript0.2480.0130.0120.248^{+0.013}_{-0.012}0.248 start_POSTSUPERSCRIPT + 0.013 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.012 end_POSTSUBSCRIPT 0.2290.079+0.059subscriptsuperscript0.2290.0590.0790.229^{+0.059}_{-0.079}0.229 start_POSTSUPERSCRIPT + 0.059 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.079 end_POSTSUBSCRIPT 0.2510.057+0.053subscriptsuperscript0.2510.0530.0570.251^{+0.053}_{-0.057}0.251 start_POSTSUPERSCRIPT + 0.053 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.057 end_POSTSUBSCRIPT
p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0.040.48+0.45subscriptsuperscript0.040.450.480.04^{+0.45}_{-0.48}0.04 start_POSTSUPERSCRIPT + 0.45 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.48 end_POSTSUBSCRIPT 0.120.23+0.21subscriptsuperscript0.120.210.23-0.12^{+0.21}_{-0.23}- 0.12 start_POSTSUPERSCRIPT + 0.21 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.23 end_POSTSUBSCRIPT 0.080.47+0.43subscriptsuperscript0.080.430.470.08^{+0.43}_{-0.47}0.08 start_POSTSUPERSCRIPT + 0.43 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.47 end_POSTSUBSCRIPT 0.030.12+0.11subscriptsuperscript0.030.110.12-0.03^{+0.11}_{-0.12}- 0.03 start_POSTSUPERSCRIPT + 0.11 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.12 end_POSTSUBSCRIPT 0.32±0.22plus-or-minus0.320.220.32\pm 0.220.32 ± 0.22 0.29±0.20plus-or-minus0.290.200.29\pm 0.200.29 ± 0.20
wbsubscript𝑤𝑏w_{b}italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT 0.022180.00076+0.00077subscriptsuperscript0.022180.000770.000760.02218^{+0.00077}_{-0.00076}0.02218 start_POSTSUPERSCRIPT + 0.00077 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00076 end_POSTSUBSCRIPT 0.022600.00033+0.00033subscriptsuperscript0.022600.000330.000330.02260^{+0.00033}_{-0.00033}0.02260 start_POSTSUPERSCRIPT + 0.00033 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00033 end_POSTSUBSCRIPT 0.0270.012+0.013subscriptsuperscript0.0270.0130.0120.027^{+0.013}_{-0.012}0.027 start_POSTSUPERSCRIPT + 0.013 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.012 end_POSTSUBSCRIPT 0.022710.00041+0.00039subscriptsuperscript0.022710.000390.000410.02271^{+0.00039}_{-0.00041}0.02271 start_POSTSUPERSCRIPT + 0.00039 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00041 end_POSTSUBSCRIPT 0.03770.0099+0.011subscriptsuperscript0.03770.0110.00990.0377^{+0.011}_{-0.0099}0.0377 start_POSTSUPERSCRIPT + 0.011 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0099 end_POSTSUBSCRIPT 0.0260.012+0.013subscriptsuperscript0.0260.0130.0120.026^{+0.013}_{-0.012}0.026 start_POSTSUPERSCRIPT + 0.013 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.012 end_POSTSUBSCRIPT
ΩmsubscriptΩm\Omega_{\mathrm{m}}roman_Ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT 0.2980.028+0.029subscriptsuperscript0.2980.0290.0280.298^{+0.029}_{-0.028}0.298 start_POSTSUPERSCRIPT + 0.029 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.028 end_POSTSUBSCRIPT 0.2920.020+0.021subscriptsuperscript0.2920.0210.0200.292^{+0.021}_{-0.020}0.292 start_POSTSUPERSCRIPT + 0.021 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.020 end_POSTSUBSCRIPT 0.2970.029+0.029subscriptsuperscript0.2970.0290.0290.297^{+0.029}_{-0.029}0.297 start_POSTSUPERSCRIPT + 0.029 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.029 end_POSTSUBSCRIPT 0.2950.014+0.015subscriptsuperscript0.2950.0150.0140.295^{+0.015}_{-0.014}0.295 start_POSTSUPERSCRIPT + 0.015 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.014 end_POSTSUBSCRIPT 0.3010.068+0.049subscriptsuperscript0.3010.0490.0680.301^{+0.049}_{-0.068}0.301 start_POSTSUPERSCRIPT + 0.049 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.068 end_POSTSUBSCRIPT 0.308±0.041plus-or-minus0.3080.0410.308\pm 0.0410.308 ± 0.041
rdsubscript𝑟dr_{\mathrm{d}}italic_r start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT [Mpc] 150.68.0+9.9subscriptsuperscript150.69.98.0150.6^{+9.9}_{-8.0}150.6 start_POSTSUPERSCRIPT + 9.9 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 8.0 end_POSTSUBSCRIPT 147.120.76+0.79subscriptsuperscript147.120.790.76147.12^{+0.79}_{-0.76}147.12 start_POSTSUPERSCRIPT + 0.79 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.76 end_POSTSUBSCRIPT 14413+14subscriptsuperscript1441413144^{+14}_{-13}144 start_POSTSUPERSCRIPT + 14 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 13 end_POSTSUBSCRIPT 147.241.0+0.88subscriptsuperscript147.240.881.0147.24^{+0.88}_{-1.0}147.24 start_POSTSUPERSCRIPT + 0.88 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 1.0 end_POSTSUBSCRIPT 145.14.6+5.1subscriptsuperscript145.15.14.6145.1^{+5.1}_{-4.6}145.1 start_POSTSUPERSCRIPT + 5.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 4.6 end_POSTSUBSCRIPT 14714+15subscriptsuperscript1471514147^{+15}_{-14}147 start_POSTSUPERSCRIPT + 15 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 14 end_POSTSUBSCRIPT
M𝑀Mitalic_M -- -- -- 19.387±0.040plus-or-minus19.3870.040-19.387\pm 0.040- 19.387 ± 0.040 19.2610.067+0.065subscriptsuperscript19.2610.0650.067-19.261^{+0.065}_{-0.067}- 19.261 start_POSTSUPERSCRIPT + 0.065 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.067 end_POSTSUBSCRIPT 19.510.22+0.20subscriptsuperscript19.510.200.22-19.51^{+0.20}_{-0.22}- 19.51 start_POSTSUPERSCRIPT + 0.20 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.22 end_POSTSUBSCRIPT
χmin2subscriptsuperscript𝜒2min\chi^{2}_{\mathrm{min}}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT 8.838.838.838.83 11.5711.5711.5711.57 23.8623.8623.8623.86 1646.151646.151646.151646.15 1607.8461607.8461607.8461607.846 5926.775926.775926.775926.77
lnijsubscript𝑖𝑗\ln\mathcal{B}_{ij}roman_ln caligraphic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT 0.1020.1020.1020.102 0.2210.221-0.221- 0.221 1.681.68-1.68- 1.68 2.462.46-2.46- 2.46 2.512.51-2.51- 2.51 2.662.662.662.66
Table 2: Constraints at 2σ𝜎\sigmaitalic_σ C.L for the f1subscript𝑓1f_{1}italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT model. For all baselines, we provide results with and without BBN constraints. Also, we include the constraints for two SNIa baselines.
Parameter DESI+BBN DESI+CMB333CMB Distance priors. DESI+CC DESI+CMB +CC+PN+superscriptPN\text{PN}^{+}PN start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT DESI+CC +PN+ DESI+CC +SNYR5
H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT [km/s/Mpc] 68.76.2+2.9subscriptsuperscript68.72.96.268.7^{+2.9}_{-6.2}68.7 start_POSTSUPERSCRIPT + 2.9 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 6.2 end_POSTSUBSCRIPT 67.92.1+1.6subscriptsuperscript67.91.62.167.9^{+1.6}_{-2.1}67.9 start_POSTSUPERSCRIPT + 1.6 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 2.1 end_POSTSUBSCRIPT 69.17.1+7.2subscriptsuperscript69.17.27.169.1^{+7.2}_{-7.1}69.1 start_POSTSUPERSCRIPT + 7.2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 7.1 end_POSTSUBSCRIPT 68.5±1.1plus-or-minus68.51.168.5\pm 1.168.5 ± 1.1 72.3±2.0plus-or-minus72.32.072.3\pm 2.072.3 ± 2.0 66.56.3+6.1subscriptsuperscript66.56.16.366.5^{+6.1}_{-6.3}66.5 start_POSTSUPERSCRIPT + 6.1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 6.3 end_POSTSUBSCRIPT
ΩcdmsubscriptΩcdm\Omega_{\mathrm{cdm}}roman_Ω start_POSTSUBSCRIPT roman_cdm end_POSTSUBSCRIPT 0.2580.027+0.030subscriptsuperscript0.2580.0300.0270.258^{+0.030}_{-0.027}0.258 start_POSTSUPERSCRIPT + 0.030 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.027 end_POSTSUBSCRIPT 0.2570.015+0.018subscriptsuperscript0.2570.0180.0150.257^{+0.018}_{-0.015}0.257 start_POSTSUPERSCRIPT + 0.018 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.015 end_POSTSUBSCRIPT 0.2540.033+0.033subscriptsuperscript0.2540.0330.0330.254^{+0.033}_{-0.033}0.254 start_POSTSUPERSCRIPT + 0.033 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.033 end_POSTSUBSCRIPT 0.251±0.011plus-or-minus0.2510.0110.251\pm 0.0110.251 ± 0.011 0.2450.039+0.041subscriptsuperscript0.2450.0410.0390.245^{+0.041}_{-0.039}0.245 start_POSTSUPERSCRIPT + 0.041 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.039 end_POSTSUBSCRIPT 0.3090.043+0.037subscriptsuperscript0.3090.0370.0430.309^{+0.037}_{-0.043}0.309 start_POSTSUPERSCRIPT + 0.037 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.043 end_POSTSUBSCRIPT
1/p21subscript𝑝21/p_{2}1 / italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT 0.280.34+0.38subscriptsuperscript0.280.380.340.28^{+0.38}_{-0.34}0.28 start_POSTSUPERSCRIPT + 0.38 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.34 end_POSTSUBSCRIPT 0.150.17+0.18subscriptsuperscript0.150.180.170.15^{+0.18}_{-0.17}0.15 start_POSTSUPERSCRIPT + 0.18 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.17 end_POSTSUBSCRIPT 0.290.34+0.37subscriptsuperscript0.290.370.340.29^{+0.37}_{-0.34}0.29 start_POSTSUPERSCRIPT + 0.37 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.34 end_POSTSUBSCRIPT 0.160.19+0.14subscriptsuperscript0.160.140.190.16^{+0.14}_{-0.19}0.16 start_POSTSUPERSCRIPT + 0.14 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.19 end_POSTSUBSCRIPT 0.380.22+0.20subscriptsuperscript0.380.200.220.38^{+0.20}_{-0.22}0.38 start_POSTSUPERSCRIPT + 0.20 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.22 end_POSTSUBSCRIPT 0.210.23+0.21subscriptsuperscript0.210.210.230.21^{+0.21}_{-0.23}0.21 start_POSTSUPERSCRIPT + 0.21 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.23 end_POSTSUBSCRIPT
wbsubscript𝑤𝑏w_{b}italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT 0.022190.00076+0.00076subscriptsuperscript0.022190.000760.000760.02219^{+0.00076}_{-0.00076}0.02219 start_POSTSUPERSCRIPT + 0.00076 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00076 end_POSTSUBSCRIPT 0.022600.00034+0.00034subscriptsuperscript0.022600.000340.000340.02260^{+0.00034}_{-0.00034}0.02260 start_POSTSUPERSCRIPT + 0.00034 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00034 end_POSTSUBSCRIPT 0.0260.011+0.011subscriptsuperscript0.0260.0110.0110.026^{+0.011}_{-0.011}0.026 start_POSTSUPERSCRIPT + 0.011 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.011 end_POSTSUBSCRIPT 0.022770.00057+0.00044subscriptsuperscript0.022770.000440.000570.02277^{+0.00044}_{-0.00057}0.02277 start_POSTSUPERSCRIPT + 0.00044 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.00057 end_POSTSUBSCRIPT 0.03360.0076+0.0073subscriptsuperscript0.03360.00730.00760.0336^{+0.0073}_{-0.0076}0.0336 start_POSTSUPERSCRIPT + 0.0073 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0076 end_POSTSUBSCRIPT 0.02150.0090+0.0096subscriptsuperscript0.02150.00960.00900.0215^{+0.0096}_{-0.0090}0.0215 start_POSTSUPERSCRIPT + 0.0096 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.0090 end_POSTSUBSCRIPT
ΩmsubscriptΩm\Omega_{\mathrm{m}}roman_Ω start_POSTSUBSCRIPT roman_m end_POSTSUBSCRIPT 0.3080.028+0.032subscriptsuperscript0.3080.0320.0280.308^{+0.032}_{-0.028}0.308 start_POSTSUPERSCRIPT + 0.032 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.028 end_POSTSUBSCRIPT 0.3060.017+0.021subscriptsuperscript0.3060.0210.0170.306^{+0.021}_{-0.017}0.306 start_POSTSUPERSCRIPT + 0.021 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.017 end_POSTSUBSCRIPT 0.3080.028+0.030subscriptsuperscript0.3080.0300.0280.308^{+0.030}_{-0.028}0.308 start_POSTSUPERSCRIPT + 0.030 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.028 end_POSTSUBSCRIPT 0.300±0.014plus-or-minus0.3000.0140.300\pm 0.0140.300 ± 0.014 0.3090.032+0.034subscriptsuperscript0.3090.0340.0320.309^{+0.034}_{-0.032}0.309 start_POSTSUPERSCRIPT + 0.034 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.032 end_POSTSUBSCRIPT 0.3570.036+0.031subscriptsuperscript0.3570.0310.0360.357^{+0.031}_{-0.036}0.357 start_POSTSUPERSCRIPT + 0.031 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.036 end_POSTSUBSCRIPT
rdsubscript𝑟dr_{\mathrm{d}}italic_r start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT [Mpc] 150.75.4+6.2subscriptsuperscript150.76.25.4150.7^{+6.2}_{-5.4}150.7 start_POSTSUPERSCRIPT + 6.2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 5.4 end_POSTSUBSCRIPT 147.300.73+0.74subscriptsuperscript147.300.740.73147.30^{+0.74}_{-0.73}147.30 start_POSTSUPERSCRIPT + 0.74 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.73 end_POSTSUBSCRIPT 14413+14subscriptsuperscript1441413144^{+14}_{-13}144 start_POSTSUPERSCRIPT + 14 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 13 end_POSTSUBSCRIPT 147.320.94+0.98subscriptsuperscript147.320.980.94147.32^{+0.98}_{-0.94}147.32 start_POSTSUPERSCRIPT + 0.98 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.94 end_POSTSUBSCRIPT 145.84.3+4.7subscriptsuperscript145.84.74.3145.8^{+4.7}_{-4.3}145.8 start_POSTSUPERSCRIPT + 4.7 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 4.3 end_POSTSUBSCRIPT 14514+15subscriptsuperscript1451514145^{+15}_{-14}145 start_POSTSUPERSCRIPT + 15 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 14 end_POSTSUBSCRIPT
M𝑀Mitalic_M -- -- -- 19.3960.044+0.033subscriptsuperscript19.3960.0330.044-19.396^{+0.033}_{-0.044}- 19.396 start_POSTSUPERSCRIPT + 0.033 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.044 end_POSTSUBSCRIPT 19.2580.066+0.062subscriptsuperscript19.2580.0620.066-19.258^{+0.062}_{-0.066}- 19.258 start_POSTSUPERSCRIPT + 0.062 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.066 end_POSTSUBSCRIPT 19.520.21+0.20subscriptsuperscript19.520.200.21-19.52^{+0.20}_{-0.21}- 19.52 start_POSTSUPERSCRIPT + 0.20 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.21 end_POSTSUBSCRIPT
χmin2subscriptsuperscript𝜒2min\chi^{2}_{\mathrm{min}}italic_χ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_min end_POSTSUBSCRIPT 8.648.648.648.64 12.4712.4712.4712.47 23.74923.74923.74923.749 1644.851644.851644.851644.85 1606.531606.531606.531606.53 5932.845932.845932.845932.84
lnijsubscript𝑖𝑗\ln\mathcal{B}_{ij}roman_ln caligraphic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT 0.8390.8390.8390.839 1.031.031.031.03 0.3210.321-0.321- 0.321 1.371.371.371.37 2.432.432.432.43 4.734.734.734.73
Table 3: Constraints at 2σ𝜎\sigmaitalic_σ C.L for the f2subscript𝑓2f_{2}italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT model. For all baselines, we provide results with and without BBN constraints. Also, we include the constraints for two SNIa baselines.

To derive our extended cosmology, we start with the the action Ferraro and Fiorini (2007); Linder (2010); Rezaei Akbarieh and Izadi (2019):

𝒮(T)=12κ2d4xe[T+(T)]+d4xem,subscript𝒮𝑇12superscript𝜅2superscriptd4𝑥𝑒delimited-[]𝑇𝑇superscriptd4𝑥𝑒subscriptm\mathcal{S}_{\mathcal{F}(T)}=\frac{1}{2\kappa^{2}}\int\mathrm{d}^{4}x\;e\left[% -T+\mathcal{F}(T)\right]+\int\mathrm{d}^{4}x\;e\mathcal{L}_{\text{m}}\,,caligraphic_S start_POSTSUBSCRIPT caligraphic_F ( italic_T ) end_POSTSUBSCRIPT = divide start_ARG 1 end_ARG start_ARG 2 italic_κ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ roman_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x italic_e [ - italic_T + caligraphic_F ( italic_T ) ] + ∫ roman_d start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT italic_x italic_e caligraphic_L start_POSTSUBSCRIPT m end_POSTSUBSCRIPT , (1)

where κ2=8πGsuperscript𝜅28𝜋𝐺\kappa^{2}=8\pi Gitalic_κ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = 8 italic_π italic_G and the tetrad determinant is calculated as e=det(eμa)=g𝑒subscriptsuperscript𝑒𝑎𝜇𝑔e=\det\left(e^{a}_{\phantom{a}\mu}\right)=\sqrt{-g}italic_e = roman_det ( italic_e start_POSTSUPERSCRIPT italic_a end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT ) = square-root start_ARG - italic_g end_ARG. When (T)0𝑇0\mathcal{F}(T)\rightarrow 0caligraphic_F ( italic_T ) → 0, we recover the concordance ΛΛ\Lambdaroman_ΛCDM model. As it is oftentimes, we consider a flat homogeneous and isotropic geometry as eμA=diag(1,a(t),a(t),a(t))subscriptsuperscript𝑒𝐴𝜇diag1𝑎𝑡𝑎𝑡𝑎𝑡e^{A}_{\phantom{A}\mu}=\text{diag}\left(1,\,a(t),\,a(t),\,a(t)\right)\,italic_e start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT = diag ( 1 , italic_a ( italic_t ) , italic_a ( italic_t ) , italic_a ( italic_t ) ) Krššák and Saridakis (2016); Tamanini and Boehmer (2012), where a(t)𝑎𝑡a(t)italic_a ( italic_t ) is the scale factor. Using the relationship between the metric and the tetrad gμν=eμAeνBηABsubscript𝑔𝜇𝜈subscriptsuperscript𝑒𝐴𝜇subscriptsuperscript𝑒𝐵𝜈subscript𝜂𝐴𝐵g_{\mu\nu}=e^{A}_{\phantom{A}\mu}e^{B}_{\phantom{B}\nu}\eta_{AB}\,italic_g start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT = italic_e start_POSTSUPERSCRIPT italic_A end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_μ end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_ν end_POSTSUBSCRIPT italic_η start_POSTSUBSCRIPT italic_A italic_B end_POSTSUBSCRIPT, we can write the flat Friedmann–Lemaître–Robertson–Walker (FLRW) metric as

ds2=dt2a2(t)(dx2+dy2+dz2),dsuperscript𝑠2dsuperscript𝑡2superscript𝑎2𝑡dsuperscript𝑥2dsuperscript𝑦2dsuperscript𝑧2\mathrm{d}s^{2}=\mathrm{d}t^{2}-a^{2}(t)\left(\mathrm{d}x^{2}+\mathrm{d}y^{2}+% \mathrm{d}z^{2}\right),roman_d italic_s start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT = roman_d italic_t start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_a start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( italic_t ) ( roman_d italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_d italic_y start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + roman_d italic_z start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) , (2)

with H=a˙/a𝐻˙𝑎𝑎H=\dot{a}/aitalic_H = over˙ start_ARG italic_a end_ARG / italic_a. Subsequently, we can derive the Friedmann equations:

H2+T3T6superscript𝐻2𝑇3subscript𝑇6\displaystyle H^{2}+\frac{T}{3}\mathcal{F}_{T}-\frac{\mathcal{F}}{6}italic_H start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + divide start_ARG italic_T end_ARG start_ARG 3 end_ARG caligraphic_F start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT - divide start_ARG caligraphic_F end_ARG start_ARG 6 end_ARG =κ23ρabsentsuperscript𝜅23𝜌\displaystyle=\frac{\kappa^{2}}{3}\rho\,= divide start_ARG italic_κ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 3 end_ARG italic_ρ (3)
H˙(1T2TTT)˙𝐻1subscript𝑇2𝑇subscript𝑇𝑇\displaystyle\dot{H}\left(1-\mathcal{F}_{T}-2T\mathcal{F}_{TT}\right)over˙ start_ARG italic_H end_ARG ( 1 - caligraphic_F start_POSTSUBSCRIPT italic_T end_POSTSUBSCRIPT - 2 italic_T caligraphic_F start_POSTSUBSCRIPT italic_T italic_T end_POSTSUBSCRIPT ) =κ22(ρ+p),absentsuperscript𝜅22𝜌𝑝\displaystyle=-\frac{\kappa^{2}}{2}\left(\rho+p\right),= - divide start_ARG italic_κ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG ( italic_ρ + italic_p ) , (4)

where ρ𝜌\rhoitalic_ρ and p𝑝pitalic_p, are the energy density and pressure, respectively. We selected f(T)𝑓𝑇f(T)italic_f ( italic_T ) cases where it is possible to reproduce naturally a late-time cosmic acceleration:

  • Power Law Model (f1subscript𝑓1f_{1}italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT)Bengochea and Ferraro (2009). – This model is of the form: f1(T)=(T)p1,subscript𝑓1𝑇superscript𝑇subscript𝑝1f_{1}(T)=\left(-T\right)^{p_{1}}\,,italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_T ) = ( - italic_T ) start_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT , where p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is a constant. We can recover ΛΛ\Lambdaroman_ΛCDM model when p1=0subscript𝑝10p_{1}=0italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0. Otherwise, if p1=1subscript𝑝11p_{1}=1italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 1, the extra term gives a re-scaled gravitational constant related to the GR limit. Furthermore, when p1<1subscript𝑝11p_{1}<1italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT < 1 gives an accelerating universe. To compare the new constraints for this model using DESI 2024, in (Xu et al., 2018) was considered BAO measurements from Two-Degree Field Galaxy Redshift Survey (2dFGRS) and SDSS DR7, where it was found that H0=69.4±0.8subscript𝐻0plus-or-minus69.40.8H_{0}=69.4\pm 0.8italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 69.4 ± 0.8[km/s/Mpc], Ωm=0.298±0.07subscriptΩ𝑚plus-or-minus0.2980.07\Omega_{m}=0.298\pm 0.07roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.298 ± 0.07 and p1=0.100.07+0.09subscript𝑝1subscriptsuperscript0.100.090.07p_{1}=-0.10^{+0.09}_{-0.07}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = - 0.10 start_POSTSUPERSCRIPT + 0.09 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.07 end_POSTSUBSCRIPT, where it is suggested that a deviation from the ΛΛ\Lambdaroman_ΛCDM model is present in the datasets.

  • Linder Model (f2subscript𝑓2f_{2}italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT)Linder (2010). – This model is described as: f2(T)=T0(1e[p2T/T0]),subscript𝑓2𝑇subscript𝑇01superscript𝑒delimited-[]subscript𝑝2𝑇subscript𝑇0f_{2}(T)=T_{0}(1-e^{[-p_{2}\sqrt{T/T_{0}}]})\,,italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_T ) = italic_T start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( 1 - italic_e start_POSTSUPERSCRIPT [ - italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT square-root start_ARG italic_T / italic_T start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG ] end_POSTSUPERSCRIPT ) , where p2subscript𝑝2p_{2}italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT is a constant and T0=T|t=t0=6H02subscript𝑇0evaluated-at𝑇𝑡subscript𝑡06superscriptsubscript𝐻02T_{0}=T|_{t=t_{0}}=-6H_{0}^{2}italic_T start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = italic_T | start_POSTSUBSCRIPT italic_t = italic_t start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUBSCRIPT = - 6 italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT. Notice that this model recovers ΛΛ\Lambdaroman_ΛCDM in the limit p2+subscript𝑝2p_{2}\rightarrow+\inftyitalic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT → + ∞. As in the latter case, this model was tested using BAO from 2dFGRS and SDSS obtaining H0=69.6±0.9subscript𝐻0plus-or-minus69.60.9H_{0}=69.6\pm 0.9italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 69.6 ± 0.9[km/s/Mpc], Ωm=0.296±0.07subscriptΩ𝑚plus-or-minus0.2960.07\Omega_{m}=0.296\pm 0.07roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.296 ± 0.07 and 1/p2=0.130.11+0.091subscript𝑝2subscriptsuperscript0.130.090.111/p_{2}=0.13^{+0.09}_{-0.11}1 / italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0.13 start_POSTSUPERSCRIPT + 0.09 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.11 end_POSTSUBSCRIPT (Xu et al., 2018), again, denoting an interesting deviation, yet small deviation from the standard cosmological model.

Refer to captionRefer to captionRefer to caption
Figure 1: 1σ𝜎\sigmaitalic_σ and 2σ𝜎\sigmaitalic_σ Confidence Levels (C.L) and posterior distributions including H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and Ωm,0subscriptΩ𝑚0\Omega_{m,0}roman_Ω start_POSTSUBSCRIPT italic_m , 0 end_POSTSUBSCRIPT. The baselines are indicated in colours for each case. Left: For the ΛΛ\Lambdaroman_ΛCDM model. Middle: For the power law model f1subscript𝑓1f_{1}italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT. Right: For the Linder model f2subscript𝑓2f_{2}italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT.

We implement each of f(T)𝑓𝑇f(T)italic_f ( italic_T ) cosmological model described and test them using the constraining parameters method through MCMC analysis using emcee 444emcee.readthedocs.io for the cosmology and the baselines with the extract of constraints using GetDist555getdist.readthedocs.io. Additionally, we assume flat priors on the set of {Ωbh2subscriptΩ𝑏superscript2\Omega_{b}h^{2}roman_Ω start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_h start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT, ΩcsubscriptΩ𝑐\Omega_{c}roman_Ω start_POSTSUBSCRIPT italic_c end_POSTSUBSCRIPT, H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT}. The baselines considered in our analysis are: (a) DESI(BAO) release obtained from observations of galaxies and quasars Adame et al. (2024a), and Lyman-α𝛼\alphaitalic_α Adame et al. (2024b) measurements. These trasers are described through the transverse comoving distance DM/rdsubscript𝐷Msubscript𝑟dD_{\mathrm{M}}/r_{\mathrm{d}}italic_D start_POSTSUBSCRIPT roman_M end_POSTSUBSCRIPT / italic_r start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT, the angle-averaged distance DV/rdsubscript𝐷Vsubscript𝑟dD_{\mathrm{V}}/r_{\mathrm{d}}italic_D start_POSTSUBSCRIPT roman_V end_POSTSUBSCRIPT / italic_r start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT, where rdsubscript𝑟dr_{\mathrm{d}}italic_r start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT is the comoving sound horizon at the drag epoch, and the Hubble horizon DH/rdsubscript𝐷Hsubscript𝑟dD_{\mathrm{H}}/r_{\mathrm{d}}italic_D start_POSTSUBSCRIPT roman_H end_POSTSUBSCRIPT / italic_r start_POSTSUBSCRIPT roman_d end_POSTSUBSCRIPT. (b) CMB Planck-2018 distant priors, which provide information on three parameters: the shift parameter R𝑅Ritalic_R that measures the peak spacing of the temperature in the CMB spectrum, the acoustic scale lAsubscript𝑙Al_{\mathrm{A}}italic_l start_POSTSUBSCRIPT roman_A end_POSTSUBSCRIPT from we can measure the temperature in the transverse direction, and finally the combination Ωbh2subscriptΩ𝑏superscript2\Omega_{b}h^{2}roman_Ω start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_h start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT Chen et al. (2019). (c) Cosmic Chronometers CC, which are measurements of H(z)𝐻𝑧H(z)italic_H ( italic_z ) from the relative ages of passively-evolving galaxies Jimenez and Loeb (2002). We conservatively use the galactic spectra to obtain dt/dzd𝑡d𝑧\mathrm{d}t/\mathrm{d}zroman_d italic_t / roman_d italic_z (Moresco et al., 2016). The final sample contains 31 data points up to z2similar-to𝑧2z\sim 2italic_z ∼ 2 with the covariance matrix generated given in (Moresco et al., 2020). (d) PN+ Pantheon-plus catalog Brout et al. (2022), with SH0ES Cepheid host distances calibrators Riess et al. (2022), and (e) DES-SN5YR Type Ia supernovae measured during the full 5-years of DES Supernova Program, which includes 1635 SNIa in the redshift range 0.10<z<1.130.10𝑧1.130.10<z<1.130.10 < italic_z < 1.13 Abbott et al. (2024).

We divide our analysis into these baselines since H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT are degenerate in the DESI BAO release. Due to this degeneracy, we will test the set in different schemes: (i) Using the combination of parameters ΩmsubscriptΩ𝑚\Omega_{m}roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT and rdhsubscript𝑟𝑑r_{d}hitalic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT italic_h in Mpc to avoid the degeneracy between h=H0/100[km/s/Mpc]subscript𝐻0100[km/s/Mpc]h=H_{0}/100\text{[km/s/Mpc]}italic_h = italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT / 100 [km/s/Mpc] and rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT. This yields the results with the 95% confidence intervals for the ΛΛ\Lambdaroman_ΛCDM model:

{Ωm=0.2860.026+0.028,rdh=102.6±2.5 [Mpc],DESI(BAO)\left\{\begin{matrix}\Omega_{m}=0.286^{+0.028}_{-0.026},\\ r_{d}h=102.6\pm 2.5\text{ [Mpc]},\end{matrix}\right.\text{DESI(BAO)}{ start_ARG start_ROW start_CELL roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.286 start_POSTSUPERSCRIPT + 0.028 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.026 end_POSTSUBSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT italic_h = 102.6 ± 2.5 [Mpc] , end_CELL end_ROW end_ARG DESI(BAO)

which is in 2σ2𝜎2\sigma2 italic_σ interval from the results reported by the DESI collaboration (Collaboration, 2024). Meanwhile, for the f1(T)subscript𝑓1𝑇f_{1}(T)italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_T ) model:

{Ωm=0.2820.033+0.031,rdh=102.0±4.2 [Mpc],p1=0.060.48+0.44.DESI(BAO)\left\{\begin{matrix}\Omega_{m}=0.282^{+0.031}_{-0.033},\\ r_{d}h=102.0\pm 4.2\text{ [Mpc]},\\ p_{1}=0.06^{+0.44}_{-0.48}.\end{matrix}\right.\text{DESI(BAO)}{ start_ARG start_ROW start_CELL roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.282 start_POSTSUPERSCRIPT + 0.031 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.033 end_POSTSUBSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT italic_h = 102.0 ± 4.2 [Mpc] , end_CELL end_ROW start_ROW start_CELL italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0.06 start_POSTSUPERSCRIPT + 0.44 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.48 end_POSTSUBSCRIPT . end_CELL end_ROW end_ARG DESI(BAO)

we can notice that the data prefers a slightly lower fractional matter density with a similar product rdhsubscript𝑟𝑑r_{d}hitalic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT italic_h and that the free parameter for the power-law model p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT is within 2σ2𝜎2\sigma2 italic_σ region. This recovers ΛΛ\Lambdaroman_ΛCDM with a minor positive deviation. In this case, the f1(T)subscript𝑓1𝑇f_{1}(T)italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_T ) model contains a Bayes factor of lnij=+1.34subscript𝑖𝑗1.34\ln\mathcal{B}_{ij}=+1.34roman_ln caligraphic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = + 1.34, which indicates a slight preference for the ΛΛ\Lambdaroman_ΛCDM. For the f2(T)subscript𝑓2𝑇f_{2}(T)italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_T ) model:

{Ωm=0.3070.039+0.041,rdh=100.65.1+4.5 [Mpc],1/p2=0.290.33+0.37.DESI(BAO)\left\{\begin{matrix}\Omega_{m}=0.307^{+0.041}_{-0.039},\\ r_{d}h=100.6^{+4.5}_{-5.1}\text{ [Mpc]},\\ 1/p_{2}=0.29^{+0.37}_{-0.33}.\\ \end{matrix}\right.\text{DESI(BAO)}{ start_ARG start_ROW start_CELL roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.307 start_POSTSUPERSCRIPT + 0.041 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.039 end_POSTSUBSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT italic_h = 100.6 start_POSTSUPERSCRIPT + 4.5 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 5.1 end_POSTSUBSCRIPT [Mpc] , end_CELL end_ROW start_ROW start_CELL 1 / italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0.29 start_POSTSUPERSCRIPT + 0.37 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.33 end_POSTSUBSCRIPT . end_CELL end_ROW end_ARG DESI(BAO)

Contrary to the previous model, here the fractional matter exhibits an increase and a diminution in the rdhsubscript𝑟𝑑r_{d}hitalic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT italic_h parameter. For the free Linder model parameter 1/p21subscript𝑝21/p_{2}1 / italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT this dataset alone recovers the ΛΛ\Lambdaroman_ΛCDM model as 1/p201subscript𝑝201/p_{2}\to 01 / italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT → 0 in 2σ2𝜎2\sigma2 italic_σ limit. The Bayes factor must be compared to the tested one for the ΛΛ\Lambdaroman_ΛCDM model which results in lnij=+0.33subscript𝑖𝑗0.33\ln\mathcal{B}_{ij}=+0.33roman_ln caligraphic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = + 0.33, this favours the standard cosmological model. (ii) By using a prior on rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT from Planck 2018 (Aghanim et al., 2020) of rd=147.09±0.87subscript𝑟𝑑plus-or-minus147.090.87r_{d}=147.09\pm 0.87italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT = 147.09 ± 0.87[Mpc] it is possible to break the degeneracy with H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. The results within 95% confidence interval for the ΛΛ\Lambdaroman_ΛCDM model are:

{H0=69.7±1.7[km/s/Mpc],Ωm=0.286±0.029,DESI +rdCMB Planck\left\{\begin{matrix}H_{0}=69.7\pm 1.7\text{[km/s/Mpc]},\\ \Omega_{m}=0.286\pm 0.029,\end{matrix}\right.\text{DESI +}r_{d}\leavevmode% \nobreak\ \text{CMB Planck}{ start_ARG start_ROW start_CELL italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 69.7 ± 1.7 [km/s/Mpc] , end_CELL end_ROW start_ROW start_CELL roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.286 ± 0.029 , end_CELL end_ROW end_ARG DESI + italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT CMB Planck

Remarkably interesting, these results have a high H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT value even though we are using a rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT from the Planck estimations. For the f1(T)subscript𝑓1𝑇f_{1}(T)italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_T ) model, the parameters are:

{H0=69.42.9+3.0[km/s/Mpc],Ωm=0.2810.034+0.031,p1=0.050.48+0.45,DESI +rdCMB Planck\left\{\begin{matrix}H_{0}=69.4^{+3.0}_{-2.9}\text{[km/s/Mpc]},\\ \Omega_{m}=0.281^{+0.031}_{-0.034},\\ p_{1}=0.05^{+0.45}_{-0.48},\end{matrix}\right.\text{DESI +}r_{d}\leavevmode% \nobreak\ \text{CMB Planck}{ start_ARG start_ROW start_CELL italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 69.4 start_POSTSUPERSCRIPT + 3.0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 2.9 end_POSTSUBSCRIPT [km/s/Mpc] , end_CELL end_ROW start_ROW start_CELL roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.281 start_POSTSUPERSCRIPT + 0.031 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.034 end_POSTSUBSCRIPT , end_CELL end_ROW start_ROW start_CELL italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT = 0.05 start_POSTSUPERSCRIPT + 0.45 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.48 end_POSTSUBSCRIPT , end_CELL end_ROW end_ARG DESI + italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT CMB Planck

where the H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT value shows a compatibility in 2σ2𝜎2\sigma2 italic_σ with the value obtained by the SH0ES collaboration (Riess et al., 2022). However, it is important to note that this result is originated due to the increased error bar in the parameter. This model returns a confirmation of ΛΛ\Lambdaroman_ΛCDM for the p1subscript𝑝1p_{1}italic_p start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT value with a significant error bars, probably because this dataset alone can not constrain the parameter solely. In this case the Bayes factor lnij=0.23subscript𝑖𝑗0.23\ln\mathcal{B}_{ij}=-0.23roman_ln caligraphic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = - 0.23 which suggests that using a rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT prior to the dataset and the power-law model is slightly favoured over the cosmological standard model. For the f2(T)subscript𝑓2𝑇f_{2}(T)italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_T ) model the results are:

{H0=68.33.5+3.0[km/s/Mpc],Ωm=0.3060.029+0.032,1/p2=0.290.32+0.36,DESI +rdCMB Planck\left\{\begin{matrix}H_{0}=68.3^{+3.0}_{-3.5}\text{[km/s/Mpc]},\\ \Omega_{m}=0.306^{+0.032}_{-0.029},\\ 1/p_{2}=0.29^{+0.36}_{-0.32},\end{matrix}\right.\text{DESI +}r_{d}\leavevmode% \nobreak\ \text{CMB Planck}{ start_ARG start_ROW start_CELL italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 68.3 start_POSTSUPERSCRIPT + 3.0 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 3.5 end_POSTSUBSCRIPT [km/s/Mpc] , end_CELL end_ROW start_ROW start_CELL roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT = 0.306 start_POSTSUPERSCRIPT + 0.032 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.029 end_POSTSUBSCRIPT , end_CELL end_ROW start_ROW start_CELL 1 / italic_p start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT = 0.29 start_POSTSUPERSCRIPT + 0.36 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.32 end_POSTSUBSCRIPT , end_CELL end_ROW end_ARG DESI + italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT CMB Planck

that, similarly to the previous model, confirm ΛΛ\Lambdaroman_ΛCDM at 2σ2𝜎2\sigma2 italic_σ level. In this case, the value of ΩmsubscriptΩ𝑚\Omega_{m}roman_Ω start_POSTSUBSCRIPT italic_m end_POSTSUBSCRIPT presents a higher value that is a tendency using this specific model. This model presents a Bayes factor lnij=0.19subscript𝑖𝑗0.19\ln\mathcal{B}_{ij}=-0.19roman_ln caligraphic_B start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = - 0.19 that, again, suggests that the f2(T)subscript𝑓2𝑇f_{2}(T)italic_f start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ( italic_T ) model is preferred over ΛΛ\Lambdaroman_ΛCDM. However, it is important to notice that in each of these models the value of the evidence between the f(T)𝑓𝑇f(T)italic_f ( italic_T ) models and the standard cosmology is not strong enough. (iii) Using a prior on wb=Ωbh2subscript𝑤𝑏subscriptΩ𝑏superscript2w_{b}=\Omega_{b}h^{2}italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = roman_Ω start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT italic_h start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT using the results of BBN presented in (Collaboration, 2024) of wb=0.02218±0.00055subscript𝑤𝑏plus-or-minus0.022180.00055w_{b}=0.02218\pm 0.00055italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT = 0.02218 ± 0.00055 to break the degeneracy. In this case, we calculate rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT as a derived parameter. This analysis is presented in Tables 1, 2, 3 including DESI + BBN. (iv) Finally, since the uncertainties on this release are substantial we will consider other datasets without the necessity to introduce a prior on wbsubscript𝑤𝑏w_{b}italic_w start_POSTSUBSCRIPT italic_b end_POSTSUBSCRIPT as the baselines are sufficient enough to constraint the cosmological parameters. These results are reported in Tables 1, 2, 3 in combination with other baselines including DESI BNN measurements.


In conclusion, f(T)𝑓𝑇f(T)italic_f ( italic_T ) cosmologies constrained by new BAO measurements from DESI 2024 (and another dataset that allow us to constraint H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT as the rdsubscript𝑟𝑑r_{d}italic_r start_POSTSUBSCRIPT italic_d end_POSTSUBSCRIPT prior from the CMB Planck, BBN, Pantheon+ or CMB Distance Priors) can be a good alternative to explain the current H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT tension as the results using this dataset combinations show an improvement in the alleviation on the H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT value closer to the SH0ES collaboration. Furthermore, it is important to notice that the mentioned combinations of DESI BAO with other datasets such as CC, Pantheon+, and even the CMB Distance Priors, the statistics show a slight preference for the f1subscript𝑓1f_{1}italic_f start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT model. This preference is in addition to the aforementioned advantage of alleviating the Hubble tension.

New analyses will be conducted using the data released in the coming months, employing these extended gravity models. Finally, this result from DESI BAO 2024 measurements could be a hint that the cosmological tension needs new physics to be solved.

Acknowledgements.


Acknowledgments.- CE-R is supported by the CONACyT Network Project No. 376127 and acknowledges the Royal Astronomical Society as Fellow FRAS 10147. RS is supported by the CONACyT National Grant. This research has been carried out using computational facilities procured through the Cosmostatistics National Group ICN UNAM project. This article is based upon work from COST Action CA21136 Addressing observational tensions in cosmology with systematics and fundamental physics (CosmoVerse) supported by COST (European Cooperation in Science and Technology).

References