A critical discussion on the $H_0$ tension

A critical discussion on the $H_0$ Hubble constant tension is presented by considering both early and late-type observations. From recent precise measurements, discrepancies emerge when comparing results for some cosmological quantities obtained at different redshifts. We highlight the most relevant measurements of $H_0$ and propose potential ideas to solve its tension. These solutions concern the exploration of new physics beyond the $\Lambda$CDM model or the evaluation of $H_0$ by other methods. In particular, we focus on the role of the look-back time.


I. INTRODUCTION
The Λ Cold Dark Matter (ΛCDM) model is considered the cosmological standard, capable of describing the observed Universe by fixing only six free parameters [1].These are the dark matter density, the baryon density, the observed angular size of the sound horizon at recombination, the scalar spectral index, the curvature fluctuation amplitude, and the reionization optical depth.
Applying general assumptions on these parameters, it is possible to derive the other cosmological quantities, including the Hubble constant H 0 and the other cosmographic parameters [2].The result is a self-consistent picture of our Universe in good agreement with observations.This relatively simple model is able to describe a large part of the Universe history with good precision, from the end of the so-called inflationary era [3,4] to the current epoch.According to the ΛCDM, our Universe is composed by three major constituents: a cosmological constant Λ, associated with the so-called dark energy, accounting for approximately 68% of the total density, a cold (non-relativistic) dark matter component, which should account for 27% of the cosmic pie [1,[5][6][7][8][9][10][11][12][13][14][15][16][17], and lastly, the remaining 5%, composed by baryonic matter, stars, galaxies, and all the luminous structures.The accuracy of this model is remarkable when compared with cosmological observations, such as the accelerated expansion of the Universe [7] deduced from the observed light curves of Supernovae Type Ia (SNe Ia).Cosmic acceleration can be addressed by the presence of a cosmological constant or, in general, by some unknown form of dark energy, acting as a negative pressure in the cosmological equations.
On the other hand, dark matter was initially introduced to account for the virial theorem applied to clusters of galaxies [18].Subsequent observations revealed that it is also a fundamental component needed to explain the rotation curves of galaxies, which, otherwise, would not be well fitted by Newtonian dynamics [19,20] if only galactic baryonic components are taken into account.
Despite its overwhelming successes, the ΛCDM model presents some critical issues that captured the attention of scientific community.The most relevant challenges are the nature of dark energy and dark matter, as well as the ongoing tension of the Hubble constant derived by different measurements at different scales.After decades of precise measurements and tests [17,[21][22][23], no direct or indirect evidence of exotic particles constituting cosmic dark fluids has been found.Consequently, according to the ΛCDM model, we have no final answer, at fundamental level, on the constituents of the observed Universe for approximately 95%.
A further issue is related to the H 0 tension, that is the discrepancy between the late-type measurements of H 0 [24], usually linked to the cosmological ladder [25,26], and the early-type ones, associated with measurements of the Cosmic Microwave Background Radiation (CMBR).The most recent results from two prominent collaborations, SH0ES and Planck, report values of H 0 as 73.04 ± 1.04 km/(s Mpc) at a 68% confidence level (CL) for the former [24] and 67.4 ± 0.5 km/(s Mpc) at 68% CL for the latter [1].As it stands, there is a 5σ tension between these two measurements [24], which, in principle, should provide the same result.Furthermore, this tension extends beyond these two collaborations and involves several late and early-type observations [25].
To address these issues, different approaches have been considered like extensions of General Relativity (GR), on which the ΛCDM model is based.For instance, a particular extension, known as f (R) gravity [27][28][29][30][31][32]), has been considered in cosmological applications [33][34][35][36][37] to address different issues related to the ΛCDM model, like the late-time dark energy [14,27], and the early inflationary behavior [4].The philosophy of these approaches is that, instead of searching for new exotic ingredients, gravitational sector should be improved according to the scales.In this perspective, also the H 0 tension could be fixed improving geometry [29,[38][39][40].
Other alternatives imply Extended Theories of Electromagnetism or the improvement of the Standard Model of Particles [41][42][43].Furthermore, the H 0 tension could be also related to some fundamental quantum concepts, like the Compton Length and the Heisenberg Uncertainty Principle, applied to the cosmological setting [44,45].
Finally, under the standard of "new physics" a large amount of investigations have been pursued both in early and late Universe .
A recent research line explores the possibility of a "variable H 0 constant" i.e. the idea that the measured value of the Hubble constant might depend on the redshift (i.e. the scale) at which it is measured [75][76][77][78][79][80][81][82][83][84][85][86][87][88][89].In this context, the H 0 constant can be evaluated by the look-back time [90].The approach consists in determining H 0 at any redshift z starting from the look-back time of the related sources.
In this paper, a critical discussion on H 0 tension will be presented.The outline is the following.Section 2 provides a brief summary of the ΛCDM model.Section 3 is devoted to the most prominent measurements related to H 0 .In Section 4, we discuss the look-back time approach to the H 0 tension starting from the results in [90].Section 5, a redshift-dependent H 0 is considered.We explore, in particular, its consequences on cosmological distances.In Section 6, we discuss the results and draw conclusions.

II. A SUMMARY OF THE ΛCDM MODEL
The ΛCDM model is a straightforward byproduct of GR.Assuming the Cosmological Principle, the Universe is homogeneous and isotropic beyond a certain scale, which is, more or less, over ∼ 120 Mpc.This assumption is supported by several observations considering large sets of data [25,[91][92][93].Cosmological Principle is implemented by the Friedman-Lemaître-Robertson-Walker (FLRW) metric [94]: where Ω represents the angular component of the metric, a(t) is the scale factor, and k is the spatial curvature constant.It can be equal to −1, 0, 1 depending on the curvature of the cosmological spatial submanifold.According to Eq.( 1), the Einstein field equations can be recast as: which are the Friedman equations leading the cosmological expansion [95].These equations are complemented by the continuity equation and the equation of state defined as: Here G is the gravitational constant, p is the pressure, ρ the density, w is the cosmological equation of state parameter, equal to −1 for the ΛCDM model, and Λ is the cosmological constant.These are the equations on which dynamics of ΛCDM model is based.The scale factor can be written as a function of the redshift as follows [2]: where a 0 = 1 is the scale factor normalized at our epoch.This allows us to write the cosmological distances as a function of the redshift [96].It is possible to rewrite Eq. ( 2) in terms of the cosmological densities as follows [97] H 2 (z) where H(z) = ȧ/a is the Hubble Parameter, H 0 is the Hubble Constant, i.e. the Hubble parameter derived for z = 0, thus for the Universe at our epoch, Ω R is the radiation energy density parameter, Ω M is the matter density parameter, where both dark and luminous matter are taken into account, Ω k is the "density" associated to the curvature, being equal zero for a flat Universe, and Ω Λ is the density associated with the cosmological constant.All these quantities, apart from H(z), are derived at present epoch.From this equation, we can define This equation allows us to express cosmological distances as a function of E(z).It is worth noticing that E(z) depends only on the redshift and the densities of the today Universe, while it does not depend directly on H 0 .The luminosity distance d L (z), derived from the intrinsic luminosity and the photon flux received by a given cosmological source, is: where d M (z) is the transverse comoving distance.It is for a flat Universe with Ω K = 0.For an open Universe with Ω K > 0, it is For a closed Universe with The luminosity distance is essential in observational cosmology because it can be associated with the "standard candles".These are astrophysical objects whose intrinsic luminosity can be derived from some intrinsic physical mechanism.Such a mechanism is generally correlated with quantities that are independent of the source distance, and so can be used to measure d L (z) intrinsically.Standard candles are a key component of the cosmic distance ladder and play a crucial role in determining H 0 .
Another important tool for the estimation of cosmic distance ladder is the angular diameter distance, defined as It is important because it is linked to the "standard rulers" (i.e., astrophysical objects whose geometrical features can be deduced from their intrinsic physics).We will describe a particular probe employing this definition.Finally, another very important distance definition is linked to the look-back time, which is the time the photon takes to reach us from a certain redshift.It is defined as: It is strictly connected to the light-travel distance, i.e. the path traveled by the photon to reach us from an astrophysical source in the expanding Universe.It is The last two equations are the starting point for the analysis presented in [90], and the novel discussions presented in this work.We have to note here that, in all the cosmological distance definitions, H 0 plays the role of a normalization constant, as it is not directly involved in the integral functions, which, in turn, depend only on the different cosmological components and the redshift.

III. THE H 0 MEASUREMENTS AND THE TENSION
Over the past decades, several methods and astrophysical sources have been employed to measure H 0 with high precision allowing us to obtain, remarkably, very small uncertainties on the measurements [17,25].This new era of precision cosmology is also the main reason for the H 0 tension, which is at the center of attention of the scientific community.As previously mentioned, we observe a significant 5σ tension between the latest SH0ES and Planck collaborations' results [1,24], but, as we will see, these derivations are actually representative of two entire sets of measurements: the former of the early-type and model-dependent observations, while the latter of the direct late-type and model-independent ones.

A. Late and Early-Type measurements
Regarding the late-type measurements, many observations, including those by the SH0ES collaboration, are based on the cosmic distance ladder method.In this approach, each step builds upon the previous one through calibration methods, especially in redshift regions where multiple probes are available.This method enables us to reach relatively deep redshift ranges while preserving the precision provided by low-redshift probes, renowned for their accuracy.
The cosmic distance ladder consists of three primary steps [24,98].The first involves precise geometric distance measurements, allowing us to directly calculate the distances of nearby objects.This step is reliable because depends on a straightforward geometrical method and does not require extensive knowledge of the astrophysical probe used for distance measurement.There are three possible anchors for this first rung: Milky Way Cepheid parallaxes, detached-eclipsing binary measurements in the Large Magellanic Cloud [99], and the water-maser host NGC 4258 [100,101].These three anchors provide approximately 1%, 1.2%, and 1.5% precision in the calibration of H 0 , respectively.In recent years, a strong improvement in the Cepheid parallax measurements has been provided by the European Space Agency (ESA) Gaia mission [102][103][104][105]. Gaia, designed for astrometry, photometry, and spectroscopy, has created the most accurate 3D map of the Milky Way.The latest release, covering the first 34 months of observations, has provided the largest dataset of Cepheids ever (around 3000 Milky Way Cepheids) allowing us to measure the parallaxes and consequently the distances of Cepheids with unprecedented accuracy [106].
The second step involves primary distance indicators, often Classical Cepheids, which exhibit an intrinsic relationship between their luminosities and periods, the so-called Period-luminosity relation (PL) [107][108][109].Although this step is valuable, it may be influenced by potential systematic effects such as the metallicity dependence of the coefficients of the PL relation, which need to be addressed [110][111][112][113].
The third step includes probes like SNe Ia, which use primary distance indicators as anchors in regions where both are detected.SNe Ia, with their higher luminosities, can explore relatively high-redshift regions.
Let us discuss in more detail how the SNe Ia are employed as standard candles, given that they are one of the most important components of the cosmological ladder approach [114].The most updated SNe Ia dataset is the Pantheon+ sample [115], which was used for the latest H 0 measurements [24,116] and serves as a natural successor to the earlier Pantheon set.[117].The Pantheon+ sample is composed of 1701 lightcurves taken from 1550 different SNe Ia, covering a redshift range from z = 0.01 up to z = 2.26.
In general, the physics behind these astrophysical objects is well-understood.SNe Ia are the byproduct of the explosions of white dwarfs in binary systems exceeding the Chandrasekhar limit due to mass transfer from their companion star.Since this limit is a fundamental constraint for the stability of all white dwarfs, all SNe Ia light curves share similar features.Specifically, the prevailing model, consistent with the majority of SNe Ia observations, is the single-degenerate Chandrasekhar mass explosion [118].In this model, the white dwarf accretes mass from its less evolved companion star, typically a red giant, which has a significantly lower density, especially in its outer regions.SNe Ia light curves are well-fitted by the deflagration model shown in [119].They are primarily powered by the β-decay of the radioactive isotope 56 Ni produced during the explosion [120].
Observationally, this model predicts that SNe Ia light curves typically exhibit an absolute magnitude around M ≃ −19 [121].However, both super-luminous [122] and sub-luminous [123] SNe Ia have been observed, suggesting the involvement of more complex mechanisms, such as the delayed detonation scenario for super-luminous SNe Ia [124].Even so, a phenomenological relation has been observed between the peak magnitude of the light curve and the luminosity decline rate in each SN Ia [125], which makes them a proper standardizable candle.Their use as a cosmological probe is based on the following equation where m is the apparent magnitude of the astrophysical object, M is its absolute magnitude, and the luminosity distance is expressed in Mpc.This quantity is confronted with the detected distance modulus µ obs of the SNe Ia, from which a best-fit of the desired cosmological parameter (as well as of the absolute magnitude M) can be performed by employing the following χ-squared function where C −1 S NeIa is the inverse of the covariance matrix [81].It is worth noticing that the absolute magnitude of SNe Ia is treated as a general parameter.This introduces a degeneracy with H 0 , which can be resolved either by fixing M to a certain value or through calibration processes involving primary distance indicators like Classical Cepheids in the cosmological ladder approach [115].
Furthermore, alternative probes can be employed in the cosmological ladder approach.For instance, the tip of the red giant branch (TRGB) [126], in place of the Classical Cepheids [127], or type II SNe [128], the Tully-Fisher relation for galaxies [129], and surface brightness fluctuations [130], as substitutes for SNe Ia in the corresponding cosmological step.The surface brightness fluctuations have also been used as primary distance indicators [131]).
Regarding the results obtained from these methods, a general consensus has emerged around the values derived by the SH0ES collaboration.Indeed, their data have been reanalyzed using different statistical methods, without a significant modification of the final results [25,[132][133][134][135][136][137].Some investigations have also included Cepheids from outside our Galaxy, to address potential biases in measurements linked to specific Cepheid populations [101].
An alternative calibration of the SNe Ia, using the TRGB methodology, has produced results that differ from those provided by the SH0ES Collaboration.For example, studies by [127,138,139] derived a value of H 0 = 69.6 ± 1.9 km/(s Mpc) at a 68% confidence level, which falls between the values from SH0ES and Planck collaborations.However, there are other studies involving the use of TRGB as calibrators for SNe Ia whose results are consistent with the other late-type observations [140][141][142][143][144].This has led to ongoing discussions about TRGB-based observations [131,145], particularly regarding methodologies to account for potential systematic effects [138,141] and a possible empirical relation between different TRGB observations, similar to what has been derived for Classical Cepheids [145].
As mentioned earlier, the cosmological ladder can be used by considering also other probes, such as the surface brightness fluctuations of galaxies, as alternatives to the SNe Ia, or as an intermediate step between Cepheids and SNe Ia [146].The results remain consistent with the other late-type measurements, even if there are higher errors in the determinations of H 0 from these sources [130,147].A similar approach can also be employed for SNe type II [128,148], and the Tully-Fisher relation for galaxies [129,149,150], obtaining results which remain consistent with the SH0ES collaboration.
While the cosmological ladder is an intuitive method for determining H 0 independently from the cosmological model, it requires very precise knowledge of the astrophysical processes associated with each used probe, especially in the first steps of the ladder.This is because any potential unaccounted-for systematic effect in the first rung could propagate to subsequent ones, as they are calibrated on the preceding.Therefore, using alternating astrophysical objects for the same step becomes a crucial test, as the cross-test between different collaborations using the same probes, to identify and remove possible systematic issues.
The cosmological ladder framework is not the only possible methodology for estimating H 0 at late times.An example is the strong lensing time delay estimates, which are independent of cosmological models but do require assumptions about foreground and lens mass distributions [25].Even with this independent alternative method, the results remain consistent with the other late-type approaches [151][152][153][154].It is interesting to mention that, in [154], a decreasing trend in the measure of H 0 with the redshift has been noted, in agreement with the previously mentioned "variable H 0 constant" [75,81,90,155], which will be the focus of the following analysis.
As we can see, for the late-type measurements, different probes and different methods provide results that are generally in agreement (apart from some exceptions).We may conclude that the tension is very unlikely to be due to systematic or statistical problems in the data themselves, but rather due to a more intrinsic, physical issue.Indeed, different averages of the late Universe estimates of H 0 are in a 4.5 − 6.3σ tension with values provided by the Planck Collaboration [25,156,157].
Let us now discuss measurements of H 0 based on assumptions and observations related to the early physics of the Universe.In addition to the latest values derived by the Planck Collaboration [1], there are other independent measurements, involving CMBR, all of which consistently yield lower values for H 0 if compared to late-type observations [158][159][160].
Early-time phenomenology can be traced even at low redshift values, with notable examples derived through various probes including the Baryon Acoustic Oscillations (BAO) [161][162][163][164][165], Big Bang Nucleosynthesis measurements of the primordial deuterium [166], and weak lensing measurements [16].These probes yield H 0 estimates consistent with those of Planck Collaboration [165], and different data reanalyses support these results [73,74,[167][168][169][170][171][172][173].However, it is essential remembering that early-type measurements are model-dependent, and work within the ΛCDM scenario and the Standard Model.Without these assumptions, constraints on H 0 and other cosmological parameters are remarkably loosened [25,174].Additionally, these measurements provide estimates for all the six free parameters underlying the ΛCDM model, which, in turn, are used to derive all the other cosmological parameters from the observation of the CMBR peaks.
Let us go into more detail about the BAO.They are also utilized independently for cosmological computations [165] because they consider another kind of cosmological distance, the angular-diameter one.BAO are widely used in literature to complement various analyses with other probes, see Refs.[42,175], and in standalone cosmological computations [165].
BAO are density fluctuations of the visible baryonic matter, caused by acoustic density waves in the early primordial plasma.As such, they are relics of phenomena occurred in early times and observed at lower redshift values such as in cluster formations and galaxy distributions.These phenomena are closely related to the acoustic peaks measured from the CMBR [161] which result from cosmological perturbations generating sound waves in the relativistic plasma of the early Universe [176].
In the past decade, BAO-related measurements have significantly improved in precision [165], which has proven to be mandatory for modern cosmological applications.This is because the acoustic features in matter correlations are relatively weak and occur at large scales [161,177].
Furthermore, these acoustic peaks are associated with different behavior of ordinary and dark matter when they are solicited by perturbations.Ordinary matter expands as a spherical wave, while dark matter remains in place [178].After this event, both dark matter and baryon perturbation start the formation of large-scale structures.
Given that, the central perturbation in the dark matter dominates over the baryonic shell, the acoustic feature is manifested as a single spike in the correlation function at approximately 150 Mpc between pairs of galaxies.This scale is typically close to the sound horizon [179].
It is important to note that, given their nature, behind each BAO-related measurement, there are tens of thousands of observations regarding large structures such as galaxies or clusters of galaxies.The largest spectroscopic survey to date is the Baryon Oscillation Spectroscopic Survey (BOSS [180]), one of the main objectives of the Sloan Digital Sky Survey (SDSS)-III Collaboration [181].Indeed, this collaboration conducted spectroscopy on over 1.5 million galaxies, generating valuable BAO-related data points.This dataset was later complemented by the extended Baryon Oscillation Spectroscopic Survey (eBOSS [182]), which was the cosmological survey within the SDSS-IV [183].
As previously stated, BAO may be used as cosmological probes starting from the angular-diameter distance.However, unlike SNe Ia, BAO-related measurements can vary in their definitions.For instance, in the set composed of 16 BAO employed in [42,43,175], compiled from [162][163][164][165]184], some of the data offer information about the following quantity where d M is the transverse comoving distance defined in Eqs.(10,11,12) and H(z) is defined in Eq. (7).Also, other cosmological quantities inferred by BAO are the following parameter where h = H 0 /(100km s −1 M pc −1 ); the so-called Hubble distance and the comoving distance itself.It is important to emphasize that all these definitions are interconnected, but they are not identical.This is a relevant point when constructing the covariance matrix using these heterogeneous measurements.
Additionally, it is worth noting that the majority of the BAO measurements have been rescaled by a factor denoted as r d , that is the distance between the end of inflation and the decoupling of baryons from photons after the recombination epoch.The value of this factor is approximately 150M pc and it is defined as [165].
where c s (z) is the sound speed, while z d is the redshift of the drag epoch, which in turn corresponds to the time when baryons decouple from the photons.This decoupling typically occurs at a redshift of z ≈ 1020, a value influenced by the physics of the early Universe.This quantity can be approximated using a formula involving cosmological parameters [185], that is where Ω b is the baryonic density and Ω ν is the neutrino density.It is important noticing that, although they are observed in lower redshift regions, their link to early-Universe physics implies that the cosmological computations derived from them are consistent with the Planck Collaboration results, including those related to H 0 [165].

B. Overcoming the tension
The H 0 tension is one of the most compelling problems of modern cosmology, and, as such, both observational and theoretical approaches have been explored by the scientific community to address it.From the former point of view, a new independent window has opened by the observations of gravitational waves [186].Indeed, these detections have already been used as 'standard sirens', to derive new estimates for H 0 [187][188][189][190].
At present, the precision of these measurements does not allow us to understand if they reduce the tension or if they are more in agreement with either the early or late-type measurements.However, given that we are still in the early stages of the gravitational waves era, a considerable improvement is expected from future observations, especially from the next generations of gravitational wave detectors.This holds great potential as a completely independent, non-electromagnetic method for measuring H 0 , providing a new window to solve the tension issue.
Other important measurable quantities for the cosmological studies are the look-back time and the related age of the Universe, which have been used to infer H 0 observationally, as it is the case of the ages of the observed astrophysical objects [155,[191][192][193][194][195].
Furthermore, interesting observations of high redshift galaxies have been performed by the James Webb Space Telescope.The observed galaxies appear to have unexpectedly high stellar masses which may be in conflict with the age of the Universe as inferred by Planck [196].It remains unclear whether this potential discrepancy may be attributed to galaxy evolution models or if it has a cosmological origin.
Another possible independent method to infer H 0 is based on the H(z) function and its evolution with the redshift, allowing us to extrapolate H 0 by requiring z = 0.This would be a completely model-independent procedure being able to derive this value, independently from other methodologies.A possible approach for this analysis is based on the so-called cosmic chronometers [26,197], i.e. the age evolution of galaxies, as well as on techniques like Gaussian Process regression [198,199], or cosmography via different polynomials [14,36,200].These approaches have provided estimates of H 0 ranging from values consistent with Planck collaboration to values consistent with late-type estimates [25].It is important to note that an extrapolation of H 0 at z = 0 from the H(z) function could be in contrast with assuming a variable H 0 .We will discuss this point in the next section.
From an observational standpoint, another potential approach to address the H 0 tension involves the extension of cosmological ladder to higher redshift, in view to bridge the gap between early and late-type measurements.For this aim, one needs to observe astrophysical objects at high z acting as standard candles.An example is represented by the Gamma-Ray Bursts (GRBs), for which different correlations between their intrinsic physical parameters can be observed [201][202][203][204][205], allowing us to use them as distance indicators [71,72,175,[206][207][208][209][210][211].
Other promising high-redshift indicators are Quasars [212].For these objects, empirical correlations among physical parameters have been found, so, as in the case of GRBs, they could constitute a formidable and populated set of objects to test the Universe at high redshift [213][214][215][216][217][218].
From a more theoretical point of view, the idea that new physics could be behind the tension is fascinating the scientific community, especially if one considers that this is not the only issue that the ΛCDM model, and more in general GR itself, presents [17,29].Other notable examples are the nature of dark energy and dark matter that seems to escape any probe at fundamental level.A popular approach to address these issues is to consider extensions of ΛCDM model and GR.As previously mentioned, f (R) gravity [29] and other modified theories [17,39] have been applied for various cosmological and astrophysical tests [17,39,82,84,85,219,220].
These modifications include the possibility of treating the dark energy component as a variable fluid rather than a cosmological constant [221,222].This constant can be represented by a scalar field ϕ rolling slowly down a flat component of a potential V(ϕ) and giving rise to the models known as quintessence [223].In this sense, the Chevallier-Polanski-Linder (CPL) parameterization for the dark energy component [224,225] is one of the most widely studied modifications to the standard scenario.
Other possible approaches include models where interactions between dark energy and dark matter are taken into account [226,227].The main issue of this framework is the degeneracy existing between different models, which try to address the same problems in completely different ways, thus not allowing us to achieve a natural and definitive extension of the ΛCDM model and GR [228][229][230].
New physics may also be linked to modifications of the underlying phenomenology at both the early and late stages of the Universe [25,88].Examples for early epochs include: 1. Early Dark Energy, which behaves as a cosmological constant for z ≤ 3000, and then decades as fast as the radiation density (or even faster) at late times [50,52,231] via a slow-roll phase transition.While promising, this approach presents problems from both observational and theoretical perspectives [232].Therefore, it has been proposed a modification, called New Early Dark Energy, where instead of a slow phase transition, we have an almost instantaneous one [233].This idea is similar to the aforementioned quintessence for late times.
2. Extra relativistic degrees of freedom at the recombination, parameterized by the number of neutrino species, N e f f .According to our current understanding, for active massless neutrino families N e f f ∼ 3.044 [234].This number affects the inferred value of H 0 .Various models regarding further dark radiation have been proposed [235][236][237][238].
3. Modifying the recombination history, by shifting the sound horizon for BAO at recombination.This can be achieved by either varying the early-time expansion history or by modifying the redshift of recombination.
Various methods have been proposed to accomplish this, including exotic scenarios in the early Universe [239][240][241].
The proposed new physics at late times modifies our interpretation of dark energy in various ways.Apart from what we have mentioned before, we recall the following: 1. Considering a cosmological bulk viscous fluid, characterized by a peculiar form of its pressure term, which is made up of two parts, where the first one is the usual component linked to the density via an equation of state, while the second is linked to the viscosity [242][243][244][245].
2. A chameleon field for dark energy, whose mass varies in accordance with the matter density of the considered environment, and whose variability would imply a measurement on H 0 dependent on the particular region in which has been performed [246][247][248][249][250].
3. Diffusion models, implying an interaction between dark energy and dark matter via a non-conserved energy tensor T µν [251][252][253], which seems to reduce the H 0 tension with different types of matter fields [254].
4. General dynamical behavior for dark energy, following a philosophy similar to the CPL parameterization.
In this sense, one could define emergent dark energy, which has had no effects in the early stages of the Universe, as it completely emerges at late times [255][256][257].
5. A Running Vacuum model, linked to possibly Quantum Field Theory or String Theory, could be used to explain theoretically a possible phenomenological dependence of cosmological and gravitational constants with the cosmic time [258,259].This kind of models can actually encompass different assumptions regarding the behavior of dark energy.It has also been successfully tested [260,261].
6.The presence of local inhomogeneities that could affect the late-time measurements of H 0 , which may be either due to possible observational issues like incomplete sky sampling, astrophysical problems like incorrect modelling of the local structures, or a more fundamental nature, like the departure of the FLRW assumption at very small scales [262][263][264][265][266] Alternatives to extensions of GR have also been sought.These alternatives include considerations for potential extensions of the Maxwellian Electromagnetism [267] in a cosmological setting, introducing a second, optical component to the measured redshift of astrophysical sources in cosmological models without dark energy [41][42][43].
In order to represent a valid alternative to the ΛCDM model, evidence of this kind of effects has to be found.A possible approach is to investigate upper limits on the photon mass [268][269][270][271], especially considering that some extensions propose massive photons [272,273].
Another investigated possibility is dealing with the H 0 tension as an evidence of a more fundamental limit on observations, linked to Quantum Mechanical concepts like the Compton Length and the Heisenberg Uncertainty Principle [44,45].This approach seeks for addressing the tension without attributing it to unlikely experimental errors or unknown novel physics.
The previously mentioned alternatives offer promising solutions for addressing the H 0 tension and providing a satisfactory explanation for its existence.The main issue is to determine if one of these alternatives truly resolves the tension or if the solution lies within a completely different framework.To accomplish this issue, one not only needs to provide definitive evidence of peculiar effects beyond the ΛCDM model, but also to build up an appropriate extension of the standard framework containing not only novel ideas to solve the H 0 and other tensions but also retaining the outstanding success achieved by GR and ΛCDM model when they are compared with observations.
As said, a novel approach considers the possibility of an evolution of H 0 constant with redshift as a way to address the tension [75][76][77][78][79][80][81][82][83][84][85][86][87][88].In this sense, the evolution of H 0 and the corresponding tension are considered as FIG.1: A summary of H 0 measurements performed both at early and late times.Credits to [25].a sort of "diagnostics" of a symptom of the breaking-down of the ΛCDM model (or even of the FLRW metric, because it would suggest, according to this interpretation, that the parameterization H(z) = H 0 E(z) is not working), marking the points in which observations are not consistent with the model [75][76][77][78][79][80].
In this context, a dependence of H 0 on the redshift has been observed in real SNe Ia data and interpreted through the f (R) gravity formalism [81][82][83][84][85]88].More specifically, in [81], a functional form has been assumed and fitted with the Pantheon set of SNe Ia [117].Results show that α is not consistent with 0 within 1σ, hinting at a smooth, slow, but continuous, decrease of H 0 value with the redshift.A variable H 0 could be due to a possible break-down, at some scale, of the Cosmological Principle on which the ΛCDM model is based (for general reviews, see [25,229,274], for discussions on diagnostics, see [77,155,264,275].
It is worth emphasizing that H 0 measurements discussed in this section represent only a part of those obtained in recent years.For a more comprehensive overview, refer to Figure 1 taken from [25].

IV. THE H 0 TENSION AND THE LOOK-BACK TIME A. The T (z) parameterization
Let us now delve into the approach proposed in [90], where H 0 is derived from the look-back time defined in Eq. ( 14).We will further discuss and improve this parameterization.Let us start from showing that, operatively, we can take into account a Taylor expansion of the scale factor as follows We can assume that a 0 = 1 and approximate H 0 = 1/T 0 , where T 0 is the today Universe age.This approximation is at 5% if compared with the Planck data.It is easy to obtain the integral This approximation is mandatory in view to justify our approach.Clearly, considering only the first-order term, we recover the parameterization that allows us to label the universe age at various redshifts T (z) starting from T 0 .The definition of look-back time can be recast as: According to this equation, we can infer H 0 from T lt (z) at any z.This result is consistent with the age of the Universe, by considering the general definition where T 0 and T (z) are the Universe age, today and at a given redshift respectively.Considering Eq. ( 26), we can derive Eq. ( 26) can be confronted with other parameterizations.The idea is linking different ages of the Universe with today epoch, avoiding the integral time evolution, depending on E(z), as in Eq. ( 29).In other words, one can adopt a point-by-point labeling process.The most natural label that we may use is the scaling factor itself a(t), which expresses how the size of the Universe changes with its expansion.In this perspective, T (z) is a projection of T 0 at a given redshift.This is the main reason why we have operatively computed this labeling from Eq. ( 24).
As noted in [90], for z → +∞, it is lim z→+∞ z + 1 z = 1, and thus it is easy to recover Eq. ( 25) from our approach, which can also be interpreted as the definition of Universe age, denoted as T 0 .This means that the parameterization is in agreement with the age definition at high values of z.Additionally, in Ref. [90], it is demonstrated that this parameterization is remarkably consistent with different H 0 measurements, ranging from late and early epochs, where different probes are taken into account.More specifically, we report here results obtained by the Planck collaboration for the following quantities [1]: which have been compared with the following H 0 observations at 68% CL: • H 0 = 73.04±1.04km/(s Mpc) from the SH0ES collaboration, inferred by the cosmic distance ladder method considering Classical Cepheids + SNe Ia) up to z = 0.15 [24]; • H 0 = 67.4± 0.5 km/(s Mpc) from the Planck collaboration, obtained by the CMBR observations at z ∼ 1100 [1]; • H 0 = 69.9± 1.9 km/(s Mpc), obtained by using the TRGB as an anchor for SNe Ia instead of the Classical Cepheids, at z = 0.08 [127]; • H 0 = 75.8± 5.0 km/(s Mpc), derived from SNe Type II as the last step of the cosmological ladder, at z = 0.45 [128]; • H 0 = 73.3± 4.0 km/(s Mpc), derived from the Mira Variables employed as anchors of SNe Ia, at z ∼ 0.15 [276]; • H 0 = 76.0 ± 2.6 km/(s Mpc), derived from the Tully-Fisher relation for spiral galaxies, at z = 0.5 [129]; • H 0 = 73.3± 2.5 km/(s Mpc), derived from the surface brightness fluctuations for the galaxies, at z = 0.33 [130]; • H 0 = 69.5 ± 3.3 km/(s Mpc), inferred from the Large Scale Structure t eq standard ruler, and thus confronted to our computations at the redshift of equivalence z eq ∼ 3300 [277]; • H 0 = 72.0 ± 1.9 km/(s Mpc), inferred from the masers + SNe Ia and compared at z ∼ 0.15 [100]; • H 0 = 73.3± 1.8 km/(s Mpc), derived from gravitational lensed quasars, confronted at z = 0.745 [153]; • H 0 = 67.9± 1.5 km/(s Mpc), which is a measurement provided by the CMBR independently from the Planck collaboration, and as such corresponding at the reionization epoch z ∼ 1100 [159]; • H 0 = 69.6 ± 2.1 km/(s Mpc), linked to the 21 cm absorption line and corresponding at the beginning of the so-called Cosmic Dawn, i. e. when the first stars formed (z ∼ 17.2), in combination with CMBR data and considering a Chaplygin gas model for the dark sector [278]; • H 0 = 73.4± 8.8 km/(s Mpc), deduced by gravitational waves, at z = 0.438 [189].H 0 evolution with the redshift  29) are plotted against the redshift and confronted with observational data [90].For each measurement, the redshift value has been selected either at the upper limit of the redshift range of the sample or at the redshift corresponding to the specific physical phenomenon considered for the estimation (cosmic dawn, recombination, equivalence epoch, and so on).Right panel: the same plots, but considering a 5% uncertainty on our model linked to the relation between T 0 and H 0 , and also a polynomial fit on the H 0 measurements.The x-axis is reported in a logarithmic scale.We recall that the measurements have been taken from [1, 24, 100, 127-130, 153, 159, 189, 276-278] We remind that the corresponding redshift for each measurement has been determined by considering either the upper limit of the redshift range of the sample used to infer H 0 or the redshift associated with the specific physical process.Results reported in [90] are displayed in left panel of Fig. 2. They will be used as reference for our tests.
Here, we consider also the effects of 5% approximation in assuming H 0 = 1/T 0 .Essentially, error bands at 5% can be taken into account.Results are shown in the right panel of Fig. 2, where we see that the model is still consistent with the H 0 measurements.We have also performed a polynomial fit of the H 0 measurements, independent of the cosmological model.It is worth noticing that it is consistent with late and early measurements, while it results shifted with respect to the peak.This feature has to be expected due to the lack of direct H 0 in the intermediate redshift range.
It is important to discuss the value of T 0 used in the analysis.This quantity is linked to H 0 .More specifically, an anti-correlation exists between T 0 and H 0 has been reported in Ref. [279].An independent way to test H 0 is to compare measurements with estimate ages of old objects such as stars and globular clusters.
In this sense, different measurements are consistent with T 0 derived from early-type observations [191,[280][281][282][283].This suggests that the H 0 value, inferred from late-type observations, would imply a Universe that is too young according to these measurements.Therefore, novel physics might be required to reconcile H 0 with the late-type derivations.This is the first reason why we started from the values provided by the Planck Collaboration.The second one is that both T 0 and H 0 , according to Planck, are derived quantities which have been computed by the six aforementioned free parameters.Therefore, even if T 0 and H 0 are linked, we do not incur in a circularity problem [1].FIG.3: Evolution of the scale factor a(t) with time.The blue curve represents the best fit using a 9th-degree polynomial to reproduce its numerical evolution as precisely as possible, while the orange curve is the fit with a straight line.
Considering Eq. ( 24), it is easy to see that the first order expansion works at low redshifts but it should get worse at higher values of z.Notably, from Eq. ( 25), we recover T 0 at higher redshifts.Furthermore, we may note that the evolution of a(t), depending on the given cosmological eras dominated by different densities, is quite similar to a linear behavior, as shown in Fig. 3.Here we compare the evolution of a(t), obtained by a 9 th degree polynomial fit, with a straight line.We may note that the two curves are remarkably similar for a large range of T (z) values and easily converge to the above limit of H z 0 at high redshifts.To further validate our claims, we can try different parameterizations, which may be derived from different assumptions on the dominating cosmological densities in the function E(z).This means that we do not consider our labeling, but specific approximations of the cosmological models.By taking into account the definition on T (z) in Eq.( 14), via the time-evolution integral, one can derive In general, this integral has to be solved numerically, but it is possible to find simple analytical formulas linking it to T 0 for specific approximations of E(z).Let us start from a matter-dominated Universe.In this case, it is It can be solved analytically and one finds that the following parameterization is exactly valid [284]: If we introduce this new parameterization into our equations and compare the derived H (z) 0 with the actual measurements, we obtain results shown in the left panel of Fig. 4. It is worth noticing that this approach does not perform well in this case, as it produces a theoretical curve for H 0 that yields unreasonable results for low values of z, and it is not consistent with the late-type measurements, even for redshift regions where the Universe may be considered matter-dominated.
If we, instead, consider a flat Universe with only dark energy and matter components (that is what we observe for the vast majority of the Universe lifetime, neglecting the radiation contribution), it is We find, starting from the equations derived in [285] for this particular case, that T (z) can be expressed as From this equation, T 0 can be easily recovered in the limit z → 0. In other words, the relation between T (z) and T 0 is as follows: By introducing the last equation in our approach, we obtain the results shown in the right panel of Fig. (4).We observe that the derived estimate for H z 0 is independent of the redshift and aligns with the measurements provided by the Planck collaboration.However, it is not consistent with measurements of H 0 obtained at lower redshifts.It is worth noticing that this estimate depends on the ratio Ω Λ /Ω M , which, in our case, is fixed to the values provided by Planck, but can be modified for other inferred values of these quantities.
As additional case, we have considered a radiation-dominated Universe, where E(z) = Ω r (1 + z) 4 and Furthermore, without any assumption, we can directly use the integral definitions of T (z) and T 0 .The plots of these two cases are shown in Fig. (5).
In the left panel, it is reported the result obtained for the radiation-dominated Universe.Notably, we observe unreasonable values for H 0 even with respect to the matter-dominated case.On the right panel, we show the result obtained adopting the general integral for T (z), that is Eq.(31).Here, we obtain a constant value for H 0 around 67.4 km/(s Mpc), independent of the redshift.This is because, by adopting the integral, we fall into a circularity problem which gives the H 0 provided by the Planck Collaboration at all redshifts.This result is also remarkably consistent with the one obtained by Eq. ( 35) for the Universe dominated by dark energy and matter.
Interestingly, we emphasize that the parameterization presented in Eq. ( 26) is the only one consistent with both early and late-type measurements.In this framework, it can be exactly derived from an empty Universe with E(z) = Ω k (1 + z) 2 .However, it is essential to note that such a Universe is in severe disagreement with cosmological observations [1].This discrepancy highlights the conceptual difference between our labeling approach and the derivation of T (z) through the integral in Eq. (31).
In summary, the parameterization shown in Eq. ( 26) has been chosen as a labeling process connecting the age of the Universe at various epochs through a point-by-point approach.Mathematically, it can be derived from a reliable H0 evolution with the redshift, matter dominated FIG.4: Left panel: the value of H 0 derived considering T (z) = T 0 (1+z) 3/2 as a function of the redshift, confronted with observational data.Right Panel: the same comparison, but with a parameterization consistent with the Universe dominated by matter and dark energy.The x-axis is in logarithmic scale.We recall that the measurements have been taken from [1, 24, 100, 127-130, 153, 159, 189, 276-278]  H 0 evolution with the redshift, Radiation Dominated FIG.5: Left panel: The value of H 0 derived from T (z) = T 0 (1+z) 2 , plotted against the redshift and compared with observational data.Right panel: The same comparison, but taking into account the integral definitions of T 0 and T (z).The x-axis is in logarithmic scale.We recall that the measurements have been taken from [1, 24, 100, 127-130, 153, 159, 189, 276-278] approximation of the scale factor as a function of time, working also for high-redshift.Among the different tests, it stands out as the only parameterization that aligns with the H 0 observations at both early and late epochs.It is worth stressing that, apart our labelling T (z), the other tests have been performed considering other parameterizations.This does not exclude the possibility of choosing other relations between T (z) and T 0 which involve more complex functions as T (z) = T 0 /(1 + P(z)) where P(z) can be some function of the redshift (e.g. a polynomial) as considered in cosmography [2].From this more general approach, one could obtain realistic models capable of matching better the cosmic history.

B. A variable H 0 from late-type estimates
The look-back time parameterization of H 0 can be compared with late-type measurements for the cosmological parameters Ω M and Ω Λ .For this comparison, we will use results obtained by the cosmological ladder approach involving the SNe Ia of the Pantheon+ set [115,116].In particular, we consider two different sets of results regarding Ω M and Ω Λ .In both cases, we use the value for T 0 provided by the Planck results, since an estimate for this quantity cannot be derived directly from the cosmological computations reported in Ref. [116].For the first set, we start from the values obtained in [116] for a flat ΛCDM model where they imposed Ω M + Ω Λ = 1.They found Considering the best-fit values of these quantities in our model, and comparing the results with the aforementioned H 0 measurements, we find the results displayed in the left panel of Fig. 6.We note how, in general, the curve obtained by our model is consistent with the majority of the H 0 measurements.However, we notice an interesting tension with the measurement provided by the Planck collaboration, which is, instead, remarkably matched by the curve displayed in Fig. 2. Indeed, we note how, at high redshift, the value for H 0 provided by this model decreases to 66.34 km/(s Mpc), which is more than 2σ in tension with the Planck measurement.An interesting observation is that this value is smaller than the measured one, contrary to what would be expected, given that late-type measurements usually translate with higher estimates for H 0 .
From this, we may conclude that the different values provided for Ω M and Ω Λ from early and late-type measurements have a significant impact on the comparison of our model with measurements.
This conclusion is further supported by the second case we have considered, in which we have taken into account the results for Ω M and Ω Λ obtained in [116] relaxing the flatness assumption for the ΛCDM, for which while for the value of Ω r , we rely on the value provided by Planck.Using our model for these starting values, we derive the right panel of Fig. 6.We observe a shift towards higher values of H 0 with respect to the results obtained in the left panel of the same figure and note that the curve is inconsistent with many observations, particularly those concerning the late Universe.In this comparison, it is important to point out the large uncertainties on Ω M and Ω Λ values, which could play a significant role.
Let us discuss now the implications of our approach in the context of ΛCDM model.As previously mentioned, a variable H 0 can be interpreted as a possible hint for the breakdown of FLRW metric and the ΛCDM model [75,155].Such a variation can be explained in the context of ETGs [29,81].However, our analysis reveals that H 0 evolution with the redshift, using FIG.6: Left panel: the value of H 0 derived from our methodology considering T (z) = T 0 (1+z) plotted against the redshift and compared with observational data, starting from the Pantheon+ results for a flat ΛCDM model.Right panel: the same comparison, but relaxing the flatness assumption.The x-axis is in logarithmic scale.We recall that the measurements have been taken from [1, 24, 100, 127-130, 153, 159, 189, 276-278] a variable H 0 can be entirely derived within the framework of ΛCDM model, starting from some fundamental concepts and from the parameterization in Eq. ( 26), without making any a priori assumption on the evolution of H 0 .This finding does not contradict the role of H 0 in the Friedman equations because the variation of this constant is strictly connected to the redshift at which it has been measured.
A robust confirmation of such a hypothesis would be attainable by conducting independent measurements within the intermediate redshift range, where we observe the peak in our H 0 estimate as shown in Figure 2.
Notably, the Pantheon and Pantheon+ SNe Ia datasets [115,117], spanning up to z = 2.26, encompass a broad redshift range, albeit with a majority of SNe Ia situated in the low-redshift region.H 0 estimates already exist in this range, (e.g., [175,209,214]) but depend on calibration processes involving probes at lower redshifts.Consequently, promising new insights are anticipated from forthcoming investigations utilizing novel probes, such as quasars, gravitational wave standard sirens, galaxy clusters, Lyman-α lines, or GRBs (as reported in [186,210,211,215,286,287]).
In this context, the Euclid Mission [288,289] could play a pivotal role.While its main focus is investigating the nature of dark energy via a wide set of observations of galaxy clusters and weak lensing phenomena [290] by also testing possible modifications of GR [291,292], it also holds the potential to provide precise estimates of cosmological parameters like H 0 within the intermediate redshift range.Additionally, given the wide survey expected by Euclid, even though not primarily centered on transient phenomena such as SNe Ia, is likely to contribute significantly also to these types of observations [293].
Let us discuss now the implications of a variable H 0 on the cosmological observations.To achieve this, we have investigated, as an example, how luminosity and light-travel distances behave with the H z 0 function in Eq. ( 27), and compared these distances with the results obtained by assuming the fixed values measured by the SH0ES and Planck collaborations.We selected the first distance because it is arguably the most relevant for cosmological  27)), with the distances derived by using the fixed values for H 0 provided by the SH0ES and Planck collaborations.Right panel: the same comparison as in the left panel, but considering the luminosity distance.The x-axis is in logarithmic scale for both panels.Note that there are different scales on the x-axis in the two plots.measurements, and the second due to its connection with the look-back time.It is worth noticing that similar considerations apply to other cosmological distances.The results are presented in Figure (7).
Both panels in this figure demonstrate how our approach naturally connects the two distances derived from fixed values of H 0 without significantly changing the overall behavior in the investigated redshift range.The effect is particularly evident in the left panel, where the light-travel distance is represented.We see how the distance computed by our approach is similar to the one derived from the SH0ES observations, converging towards the Planck results at higher redshifts.This can be attributed to the role of H 0 , which is a normalization constant, not linked to any astrophysical source.In conclusion, while a variable H 0 does impact the definitions of cosmological distances, it neither substantially changes the conclusions drawn from our method, especially at lower redshifts, nor it changes the value of the distances so that it could be in contrast with observations.

VI. DISCUSSION AND CONCLUSIONS
In this work, we provided, without claiming for completeness, a summary of H 0 measurements in both early and late-type frameworks, highlighting the existence of a tension independent of the particular probe or methodology used.This independency emphasizes the decreasing probability that the tension stems from potential instrumentrelated biases or unaddressed systematic effects.
We explored how this tension has been tackled in the context of new physics, beyond the ΛCDM model.Various approaches have been proposed, yielding promising results.However, the main problem is the substantial degeneracy among the different proposed frameworks, which does not allow us, up to now, to define a unique and definitive extension of the ΛCDM model.From this point of view, it is crucial to understand that, whichever the extended framework may be, it has to be able to reproduce the ΛCDM results where the concordance among the vast majority of measurements is evident.This check is essential in view to recover self-consistent models in agreement with cosmic history.
A recent approach has brought significant attention to the possibility of a "running H 0 " [75], both in the observations context and in the interpretation of a quantity that is traditionally regarded as "constant".
In this framework, we reported the analysis performed in [90], which defines the Hubble constant via the lookback time.It is possible to provide a general formula consistent with the measurements of this quantity both at early and late epochs just considering the related look-back time measurement of H 0 at any epoch.Here, we delve further into this direction starting from Eq. ( 26) and discussing the consequences on the various cosmological distances.
Specifically, we highlight that the formula T (z) = T 0 /(1 + z) arises naturally through a point-by-point labeling of the age of the Universe at different redshifts, without the need to account for the cosmological evolution already incorporated into the function E(z).
To validate our assumption, we show how it can be mathematically derived from a reliable approximation of the scale factor a(t), and how other parameterizations, which can be derived from specific assumptions on the densities at various epochs, do not effectively fit observational measurements of H 0 .
We have also considered different values of Ω M and Ω Λ , taken by late-type results from the SH0ES+Pantheon+ sets, finding significant effects in our comparisons, particularly regarding the value for H 0 provided by the Planck collaboration, which is not recovered in these cases.This observation is interesting because it allows us to conclude that the different measurements of Ω M and Ω Λ , for late and early epochs, introduce a notable tension.
In particular, we derived the light-travel and luminosity distances as functions of our variable H 0 and compared the results with those obtained using the fixed values provided by early and late-time measurements.We found that, by our new definition, it is possible to obtain distances able to link the SH0ES and Planck distances, without significantly modifying their overall behavior.In our opinion, this is an important check because the absence of unreasonable results, starting from different distance definitions, confirms the reliability of the approach.
An issue that may arise pertains to the measurements of the aforementioned H(z) function and the extrapolation to the H 0 value.In fact, one should be able to discern between variations related to the functional form of H(z) and the variable nature of H 0 at different redshifts.
An important point that have to be emphasized is that our results have been achieved entirely within the framework of the ΛCDM scenario, without requiring modifications or extensions.Theoretically, a variable H 0 can be interpreted as a breakdown of the FLRW metric.In our case, H 0 is an integration constant related to the size of the Universe at a given redshift.In other words, the value of H 0 could depend on the redshift at which it is measured, thus not undermining its role in the cosmological equations but removing the tension issue.It is worth noticing how a better fit with respect to the measurements can be obtained by starting from the early estimate of densities rather than the late ones.However, this does not exclude the necessity of extending GR and ΛCDM model, given other issues like the nature of dark energy and dark matter, and the lack of a self-consistent Quantum Gravity theory [29].
Observationally, H 0 tension is not the only tension in cosmology [25].For example, there is the so-called S 8 = σ 8 √ Ω M /0.3 tension, where S 8 is a parameter indicating the strength with which matter is clustered in the Universe.On this parameter, a discrepancy at 2 − 3σ level [25] exists between the measurements inferred by the Planck data and low-redshift probes, such as the weak gravitational lens and clusters of galaxies [294,295].It is essential to investigate whether similar considerations to those discussed above can also apply to S 8 or if it represents an independent signal of deviations from the cosmological Standard Model.
Furthermore, other measurements challenge the ΛCDM model, like some evidence of a possible non-zero curvature of the Universe [296], as well as anomalies in CMB observations, i. e., apparent correlations between the Solar System plane and certain aspects of CMB.This evidence seems to provide a preferred reference position which should not be possible if the Cosmological Principle is always valid [297][298][299].
In conclusion, our analysis suggests that new physics, in the form of extensions or modifications to ΛCDM model, may not be necessary to address the specific issue of H 0 tension, but it does not exclude its necessity for other fundamental issues.In a forthcoming study, we will discuss the other tensions under the standard of look-back time approach.

Aknowledgements
This paper is based upon work from the COST Action CA21136, "Addressing observational tensions in cosmology with systematics and fundamental physics" (CosmoVerse) supported by COST (European Cooperation in Science and Technology).The authors acknowledge the support by Istituto Nazionale di Fisica Nucleare, Sez.di Napoli, Iniziativa Specifica QGSKY.The authors express their gratitude for valuable and constructive discussions to Giada Bargiacchi, Micol Benetti, Eoin Colgain, Maria Giovanna Dainotti, David Kraljic, Stefan Ruester, Shahin Sheikh-Jabbari, and Alessandro D.A.M. Spallicci.Finally, we want to thank the anonymous referees for their feedback and useful suggestions.

×10 6 FIG. 7 :
FIG.7: Left panel: a comparison between the light-travel distance computed using the H z 0 function (Eq.(27)), with the distances derived by using the fixed values for H 0 provided by the SH0ES and Planck collaborations.Right panel: the same comparison as in the left panel, but considering the luminosity distance.The x-axis is in logarithmic scale for both panels.Note that there are different scales on the x-axis in the two plots.
H 0 evolution with the redshift, error bars from the relation with T 0 FIG. 2: Left panel: The values of H 0 derived from Eq.( Brout et al. for M and , Flat CDM H 0 evolution with the redshift, using Brout et al. for M and