Cosmological constraints on GUP from modified Friedmann equations

The Generalized Uncertainty Principle (GUP) has emerged in numerous attempts to a theory of quantum gravity and predicts the existence of a minimum length in Nature. In this work, we consider two cosmological models arising from Friedmann equations modified by the GUP (in its linear and quadratic formulations) and compare them with observational data. Our aim is to derive constraints on the GUP parameter and discuss the viability and physical implications of such models. We find for the parameter in the quadratic formulation the constraint αQ2<1059\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha ^{2}_{Q}<10^{59}$$\end{document} (tighter than most of those obtained in an astrophysical context) while the linear formulation does not appear compatible with present cosmological data. Our analysis highlights the powerful role of high-precision cosmological probes in the realm of quantum gravity phenomenology.


Introduction
The presence of a Generalized Uncertainty Principle (GUP) appears to be an ubiquitous feature in different tentative approaches to quantum gravity and thought experiments. This modification of the Heisenberg Uncertainty Principle (HUP) is required whenever the gravitational interaction is taken into account and results in the introduction of a minimum length scale in Nature, of the order of Planck length. The existence of such a minimum length has powerful implications, since it entails the impossibility to probe lengths shorter than Planck length in principle, regardless of the capabilities of any experimental apparatus, which seems to deeply affect the very concept of spacetime.
The idea of a minimum length in Nature has a long history [1], but only in the 1960s, the pioneering work of Mead [2] pointed out the crucial role of the gravitational interaction in the possibility of probing very short distances. However, a e-mail: serena.giardino@aei.mpg.de (corresponding author) b e-mail: vincenzo.salzano@usz.edu.pl it is mostly after Hawking's revolutionary insights into the thermodynamical properties of black holes (BHs) [3] that the "trans-Planckian problem" has emerged as a relevant issue for gravity (and, after that, also for inflationary cosmology [4]). In the following years, results in String Theory brought more attention to the problem of resolving infinitely small lengths with extended objects [5][6][7]. Since then, the GUP has played an important role as a heuristic tool to understand quantum gravity effects, starting from Gedankenexperimente and disparate theoretical frameworks, even if the GUP itself seems to nonetheless have model-independent features [8].
One of the contexts in which the GUP has been studied most extensively is that of BHs: as Bronstein had already understood back in the 1930s [9], the problem of finding a theory of quantum gravity is tightly interwoven with the existence of BHs. As gravity does not allow the concentration of an infinite amount of energy into a confined region, because the process will end in gravitational collapse and the formation of a BH, there is an inescapable limit to the precision to which the gravitational field can be measured [10].
This insight is remarkably similar to the Gedankenexperiment involving the formation of micro black holes carried out in [11], which yields an intuitive explanation of the intrinsic limit to any measurement carried out around the Planck scale.
The most common expression of the GUP is given by where p = 1.61623 × 10 −35 m is the Planck length and β is the dimensionless GUP parameter, generically assumed to be O (1), because the minimum length arising from the GUP is Δx min ≈ √ β p [8]. Here, the correction to HUP is proportional to a quadratic term in the momentum uncertainty, but several other formulations of the GUP also contain a linear term [12]. GUP-induced quantum effects lead to corrections both of the Hawking temperature and the Bekenstein entropy. Because of such modifications and their influence on the flux of Hawking radiation [13,14], the GUP critically affects BH evaporation [15,16]. Especially in the final stages of this process, the role of the GUP becomes important as it could cause the black hole to leave a remnant of Planck size, thus providing hints towards the resolution of the black hole information paradox [17], and, in a cosmological context, also possibly giving rise to distinguishable effects on the CMB power spectrum [18].
An additional intriguing implication is that the correction to the entropy, translated to a modified entropy-area law [19][20][21][22][23], can also be studied in a cosmological context because of the geometrical (and thus universal) nature of such a law, which does not only apply to black hole horizons. Crucially, such results are general because any causal horizon is inevitably associated with entropy, since by definition it hides information from observers, as clarified in the seminal paper by Jacobson [24]. Showing that the first law of thermodynamics can be recast as the Einstein equations, the author provided evidence of the close relationship between thermodynamics and gravity, resulting in illuminating insights into the thermodynamical properties of space-time as a whole.
Building on this result, the authors of [25] showed that the Friedmann equations can be recovered by applying the first law of thermodynamics, d E = T d S, to the apparent horizon of Friedmann-Lemaître-Robertson-Walker (FLRW) spacetime, which is also endowed with a temperature T and an entropy S that read where A is the area andr A the radius of the apparent horizon. This procedure for obtaining the Friedmann equations seems to have quite broad validity: it was shown that it also holds for alternative theories of gravity [25] and even if the entropy-area law is generically modified [26].
More recently, the authors of [27] have shown that the Friedmann equations are still recovered (albeit in a modified form) even if the entropy-area law is affected by the GUP, which means that GUP-induced quantum effects at high energies can indeed influence the dynamics of the FLRW universe at early times, however slightly.
In addition to the plethora of theoretical investigations carried out on the GUP, a research direction rooted in quantum gravity phenomenology is attempting to quantify the magnitude of GUP-induced quantum corrections by constraining the GUP parameter. The relevance of these studies lies in the fact that they open a much needed low-energy window on quantum gravity, far from the presently inaccessible Planck scale, employing precision experiments in many areas of physics. However, the instances in which theoretical studies are complemented by comparison with experimental data, especially in a cosmological setting, are, to the best of our knowledge, scarce (exceptions are e.g. [28,29]).
The aim of the present work is precisely to bridge two cosmological models in which GUP-induced thermodynamical corrections are taken into account [27,30] with the wealth of precise cosmological data available today. Our aim is two-fold: on the one hand, we will check whether these models are compatible with data and thus study their cosmological viability; on the other hand, we will derive cosmologically-motivated constraints on the GUP parameter. Generally speaking, as clarified in Sect. 3, experimental bounds coming from the study of gravitational effects are weak, while bounds coming from quantum experiments are much more stringent. We will show that cosmological data are capable of providing fairly strong constraints, comparable in magnitude with the less stringent estimates of quantum experiments.
The paper is organized as follows: in Sects. 2 and 3, we will review the salient steps for the derivation of the modified Friedmann equations and the experimental bounds on the GUP parameter; in Sect. 4, we will describe the statistical methods and cosmological data employed in our analysis; in Sect. 5 we will present and discuss our results.

Theoretical background
The work of [25] characterises the apparent horizon of FLRW universe in analogy with the event horizon of a BH. The apparent horizon is a marginally trapped surface with vanishing expansion that always exists in FLRW universe (differently from the event and particle horizons) which makes it the best suited cosmological horizon for thermodynamical considerations, also in view of its trapping character [31]. To set our notation, we write the Friedmann equations as The continuity equation readṡ and will not be affected by the GUP corrections. Since we are interested in recovering the GUP-modified entropy-area law, we need to define the radius and the area of the apparent horizon. The radius is and yields the area The expressions for the entropy and the temperature associated to a BH horizon read where A is the area of the event horizon, M the mass of the black hole and k B the Boltzmann constant. In the well-studied case of a Schwarzschild BH, the Schwarzschild radius In [32], this result has been generalized to the case of a de Sitter universe, and in [25] it is only assumed as a working hypothesis that the expressions above also work for the apparent horizon. However, since the authors do successfully recover the Friedmann equations, their assumption is justified a posteriori. Thus, for our purposes, we can simply assume [33] T =h c 4π k BrA .

Quadratic GUP
We will now briefly review the treatment of [27], in which it is shown in an arbitrary number of dimensions that a modified form of the Friedmann equations can be recovered when the GUP corrects the standard entropy-area law (2). In the following, we will also restore all physical units which are generally omitted in theoretical studies, since our goal is to find a numerical estimate for the GUP parameter. The authors of [27] start from the following expression for the GUP, which includes a quadratic correction in momentum uncertainty (note that the dimensionless GUP parameter β in (1) is here α 2 Q , where Q stands for "quadratic") We can straightforwardly solve for the momentum uncertainty, obtaining which can be recast as The expression inside square brackets is the function characterizing the departure of the GUP from the HUP, which we define as In order to study the GUP effects on the thermodynamics of FLRW universe, the authors of [27] (relying on [34]) consider the following picture: if the apparent horizon has absorbed or radiated a particle with energy d E, this energy can be identified with the uncertainty in momentum, d E cΔp. Therefore, the HUP Δp ≥h/(2Δx) yields the corresponding increase or decrease in the area of the apparent horizon, due to (7): However, in the case where the GUP is taken into account, this relationship becomes with The uncertainty in position of the absorbed or radiated particle is reasonably considered of the order of its Compton length λ, which is approximately the inverse of the Hawking temperature in natural units. In physical units, as noted in [22], it is customary to consider Δx λ. For a Schwarzschild BH, the particle has a wavelength of the order of the inverse Hawking temperature (for an asymptotic observer) or, more generally, of the inverse of the surface gravity κ −1 = 2r S (since T = κ/2π in natural units). As previously noted, it seems sensible to extend the argument from the context of BHs to that of the apparent cosmological horizon, thus assuming Δx 2r A = √ A/π . We can express the departure function (13) in terms of the area of the apparent horizon and subsequently in terms of the entropy. If we expand f G (A) around α Q l p = 0, we obtain up to second order, which is sufficient for the purposes of this work. Substituting (17) in (16) and integrating, we find the expression where in ln A we have included the integration constant A 0 . For the entropy, due to (16) and (7), we obtain d S G = f G (A) d S. After integration, we find the modified entropyarea relation to be The apparent horizon approach to find the Friedmann equations devised in [25] consists in the application of the first law of thermodynamics to the apparent horizon of a FLRW universe, with the additional assumption (7). This procedure involves the definition of a work density and an energy supply vector, respectively regarded as the work done by a change of the apparent horizon and the total energy flow through it, which is associated to the entropy. Such definitions yield a specific form of the first law of thermodynamics for cosmological horizons. The first time derivative of (5) iṡ which can be rewritten as Taking into account (13) and (19), it is straightforward to see that where S G is a first derivative with respect to A. Following [26,35], we can use as the starting point for finding the Friedmann equations. On the one hand, using (21), we can find the dynamical Friedmann equation On the other hand, taking into account the continuity equation (4), (23) can be recast as Integrating this equation yields the Friedmann constraint, so that, in summary, the two GUP-modified Friedmann equations (up to second order in α 2 Q ) are given by

Implementing the Friedmann equations
In order to compare the model arising from these modified Friedmann equations with cosmological data, we need an expression for H (z) in terms of the cosmological parameters. Our cosmological background will be a standard ΛCMD model [36], fully characterized by the Hubble constant H 0 = 100 · h, the dimensionless density parameters for matter Ω m , for radiation Ω r and for a cosmological constant as dark energy component Ω Λ , where Ω i = 8π Gρ i /3H 2 0 and the curvature k is expressed as Ω k = −kc 2 /H 2 0 . As a first step, we solve Eq. (26) for A, choosing the positive and real solution Due to (6), we can invert to find where a = 1/(z + 1) and ρ is the total energy-matter density of the universe. Given the continuity equation (4), each component of the energy density is Therefore, it is possible to write the total energy density ρ in (29) as ρ = i Ω i ρ c a −3(1+w i ) and express the dimensionless Hubble parameter E(z) = H (z)/H 0 as where Additionally, after ensuring the normalization condition E(z = 0) = 1, we can define Ω Λ in terms of all other parameters, inverting (30) at z = 0, which yields

Linear GUP
Making use of an alternative formulation of the GUP proposed in [12], the author of [30] also performed the computation of the modified Friedmann equations with the apparent horizon formalism. More specifically, this work deals with where, in addition to the quadratic term in the momentum, a linear term appears (the GUP parameter is named α L here, where L stands for "linear", to avoid any ambiguities with the quadratic case). Following the same procedure of Sect. 2, the resulting modified Friedmann equations read where terms containing higher orders of 1/A have been neglected. Given the similarity of the Friedmann equations (34) and (35) with (26) and (27) of [27], it should be straightforward to compare them with cosmological data, adopting the method in Sect. 2.2. However, this task proved substantially more challenging with the inclusion of a linear term in the GUP, due to the fractional exponent of A involved in (34). If we assume spatial flatness, neglecting the curvature k, (34) can be recast as i.e. a third order equation in H (z). It admits one single real solution, which reads (in terms of the dimensionless Hubble parameter E(z)) where (38) and Expression (37) contains a square root, which has profound implications for our goal of finding an upper bound on the GUP parameter. Indeed, an additional condition needs to be satisfied to guarantee that E(z) is real, namely If we restrict to positive values of the GUP parameter only, this requirement imposes a lower bound on α L , which proves incompatible with the notion that General Relativity and standard Quantum Mechanics are to be recovered in the limit α L → 0. An upper bound would be theoretically possible if negative values of the GUP parameter were allowed. However, this would contrast with the idea of a minimum length Δx min ≈ √ β p (for the GUP formulation (1)) which cannot be imaginary [37]. Nonetheless, a negative GUP parameter has interesting implications, see e.g. [28,[38][39][40][41]. For the purposes of this work, we found that allowing α L to be negative leads to the unphysical result H (z) < 0. This is the reason why we only consider α L > 0 in order to test the model with cosmological data.

Experimental bounds on GUP
The diversity of approaches used to find experimental bounds on the GUP parameter is even richer than that of the theoretical frameworks they stem from [54,55]. Still, we can clearly identify at least a couple of broad contexts in which the effects of GUP have been most consistently studied [56]: that of gravitational tests and that of quantum optics and atomic experiments.
In the first context, we find for example the study of a GUP-deformed Hawking temperature and Schwarzschild metric connected to Solar System tests of GR [52] and the recent constraints on both α Q and α L provided by the gravitational wave event GW150914, through the analysis of modified dispersion relations for gravitons [51]. However, these bounds are typically quite weak. They seem to be tighter in where the operatorsX andP are believed to be valid near Planck scale. This commutator is then used to find corrections to standard quantum mechanical effects, such as the Lamb shift and the Landau levels. These works mainly use the GUP formulation introduced in [12] and have provided some amongst the most stringent bounds on the GUP parameter so far [44,45,50]. Another development in this direction, which employs the quadratic GUP formulation, recently provided probably the tightest constraint available [42]. Most of the known constraints are reported in Table 1. Comparing the different expressions of the GUP used throughout the literature, such as (1), (10) and (33), it is straightforward to notice that, for what concerns the quadratic GUP, the parameters β and α 2 Q (as in our notation) are fully equivalent. Thus, in order to compare them directly, the bounds we found on α Q have been squared and translated into bounds on α 2 Q . As for the linear GUP, we report our bound on α L , which is directly comparable with those available in the literature. Nonetheless, we note that caution may be required when contrasting bounds from completely different experiments in which both the GUP effects and the uncertainties on the bounds can greatly vary.
To the best of our knowledge, there are very few studies that aim to constrain the GUP parameter with cosmological observations. One of them is [28], which makes use of the apparent horizon formalism to inquire whether the GUP effects could account for dark energy. However, the author mostly deals with an alternative formulation of the GUP, the General Extended Uncertainty Principle (GEUP), which predicts the existence of a maximum length as well as a minimum one. The bound this study finds on the minimum length can be regarded as a bound on β (because Δx min ≈ √ β p ) and is not very stringent, as reported in Table 1. Additionally, the results of another study that uses astrophysical data [51] show that very different estimates of the GUP parameter can be found, whether the linear or the quadratic formulation is employed.

Statistical analysis
In order to test the viability of the models described in Sect. 2, we use a combination of data coming from well-known geometrical probes. More specifically, we employ: Type Ia Supernovae (SNeIa) from the Pantheon sample; Early-Type Galaxies as Cosmic Chronometers (CC); the H 0 Lenses in COSMOGRAIL's Wellspring (H0LiCOW) data [58]; the "Mayflower" sample of Gamma Ray Bursts (GRBs); Baryon Acoustic Oscillations (BAO) from several surveys; and the latest Planck 2018 release for Cosmic Microwave Background radiation. We consider two different scenarios: the "full" data set, joining both early-(CMB and BAO data from SDSS) and late-time observations (SNeIa, CC, H0LiCOW, GRBs and BAO from WiggleZ); and the "late-time" data set, which includes only late-time data. In general, the total χ 2 is the sum of all the considered contributions In order to test the predictions of the GUP-modified cosmological models given the observational data, we use our own implementation of a Monte Carlo Markov Chain (MCMC) [59][60][61] to minimise the total χ 2 . We test its convergence using the method developed in [62].
As means of comparison, we also analyse a standard ΛCDM model using the same data set. This allows us to assess the reliability of the GUP-modified models with respect to the standard case, through the computation of their Bayesian Evidence, using the algorithm in [63]. Because of the stochastic nature of the evidence, we compute it ∼ 100 times and consider the median and the 1σ confidence level obtained from its statistical distribution, reporting them in Table 2. We briefly remark that the Bayesian Evidence E is defined as the probability of the data D given the model M with a set of parameters p, i.e. E(M) = d p L(D| p, M) π( p|M), where π( p|M) is the prior on the set of parameters, normalised to unity, and L(D| p, M) ∝ exp[−χ 2 /2] is the likelihood function. In order to minimise the dependence on the priors [64] used during the MCMC, we employ the same uninformative flat priors for every parameter in each model and allow a sufficiently wide range for them.
Once the Bayesian Evidence is calculated, we define the Bayes Factor as the ratio of evidences between two models M i and M j , namely

Type Ia supernovae
The Pantheon compilation [66] contains 1048 objects within the redshift range 0.01 < z < 2.26. We define the related where Δμ = μ theo − μ obs is the difference between the theoretical and the observed value of the distance modulus for each SNeIa and C SN represents the total covariance matrix. The theoretically predicted distance modulus is defined as where is the dimensionless luminosity distance and p is the vector of cosmological parameters. We also need to marginalize over the nuisance parameter μ 0 (a combination of the Hubble constant, the speed of light c and the SNeIa absolute magnitude); following [67] we end up with: where a ≡ Δμ SN SN · 1 and 1 is the identity matrix.

Cosmic chronometers
The definition of Cosmic Chronometers (CC) applies to Early-Type Galaxies (ETGs) exhibiting a passive evolution, which provide measurements of the Hubble parameter over extended redshift ranges. The sample we use here covers the range 0 < z < 1.97 [68]. The χ 2 H in this case can be constructed as where σ H (z i ) are the observational errors on the measured values H obs (z i ).

H0LiCOW
The H0LiCOW collaboration [69] used the sensitivity of strong gravitational lensing events to constrain H 0 and the cosmological background. In particular, it focused on 6 selected lensed quasars [58] for which multiple images were provided. It is well-known that the light travel time from the quasars (sources) to the observer depends on the path length and the gravitational potential of the foreground mass (lens), so that multiple images can exhibit a time delay at collection given by As in a typical gravitational lensing configuration [70], z L is the lens redshift, θ the angular position of the image, β the angular position of the source andΨ the effective lens potential, while D S , D L and D L S are, respectively, the angular diameter distances from the source to the observer, from the lens to the observer, and between source and lens. They are given by the following expression: in the case of spatial flatness [71]. The quantity used in our analysis is generally called timedelay distance, and is defined as all data (D obs Δt,i ) and errors (σ D Δt,i ) on this quantity for each quasar are provided in [58]. Thus, the χ 2 for H0LiCOW data is

Gamma ray bursts
The "Mayflower" sample is made of 79 GRBs within the redshift range 1.44 < z < 8.1 [72]. Given that GRBs observable is again a distance modulus, the same procedure used in the section above for SNeIa also applies here; the final χ 2 G estimator is given by (46)

Baryon acoustic oscillations
The χ 2 estimator for BAO data is the sum of different contributions where the observables F B AO change according to the chosen survey. We employ data from the WiggleZ Dark Energy Survey (at redshifts 0.44, 0.6 and 0.73) [73], for which the relevant physical quantities are the acoustic parameter and the Alcock-Paczynski distortion parameter where is the angular diameter distance and is the volume distance, the geometric mean of the radial (∝ H −1 ) and tangential (D A ) BAO modes.
We consider multiple data from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). From the DR12 analysis in [74], the following quantities are given where the comoving distance D M is Here, r s (z d ) denotes the sound horizon evaluated at the dragging redshift, while r f id s (z d ) is the sound horizon calculated at a given fiducial cosmological model (in this case, it is 147.78 Mpc). The dragging redshift is estimated using the approximation [75] as where the factors b 1 and b 2 are given by respectively. The sound horizon is defined as: where the sound speed is given by , (62) and the baryon-to-photon density ratio parameters is From the DR12 we also include measurements derived from the void-galaxy cross-correlation [76] H (z = 0.57)r s (z d ) = (14.05 ± 0.14) 10 3 km/s.
From the extended Baryon Oscillation Spectroscopic Survey (eBOSS) we have used the point Mpc [77].
Additional points taken into account from eBOSS DR14 are obtained from the combination of the quasar Lyman-α autocorrelation function [78] with the cross-correlation measurement [79], namely

Cosmic microwave background
Regarding CMB data, we use the shift parameters [80] derived from the latest Planck 2018 data release [81]. In this case the χ 2 C M B is defined as where the vector F C M B contains the quantities in addition to Ω b h 2 . Here, r s (z * ) is the comoving sound horizon evaluated at the photon-decoupling redshift given by [82] by the following fitting formula, where the factors g 1 and g 2 are while r is the comoving distance defined by H (z , p) .

Results and discussion
The results of our analysis are summarised in Table 2 and complemented by the comparison with a reference model, namely the standard ΛCDM. In the upper part of the table, we show the primary cosmological parameters as constrained by our statistical analysis, i.e. Ω m , Ω b and h, with the addition of log α Q,L p /c , where the log is used for ease of numerical computation. In the lower part of the table, we show some secondary parameters which can be derived from the primary ones and include the constraints on the GUP parameters that are central to our study. More specifically, we report -Ω Λ , the cosmological constant (dark energy) contribution derived from the normalization condition E(z = 0) = 1. It allows us to explore the possibility that the GUP correction could play the role of an effective dark energy fluid; -α 2 Q ∼ β or α L (for the quadratic and linear formulations of the GUP, respectively). They represent the quantities constrained by experiments and reported in Table 1; -α 2 Q 2 p and α L p , the parameters we perturbatively expand around to find f G (A) in (17), the departure function from the HUP that underlies this entire study. Thus, a posteriori verifying the smallness of its value represents an important consistency check for our analysis and a way to assess the reliability of our conclusions; -χ 2 min , B i j and log B i j , which are the minimum value of χ 2 , the Bayes Factor and its logarithm (for comparison with Jeffreys' scale), respectively.

Quadratic GUP
The most conclusive case in our analysis is that of the quadratic GUP with the full data set. As shown by the Bayesian Factor in Table 2, this model is essentially indistinguishable from the standard ΛCDM case and the evidence in its favour is negligible. This means, as expected, that the GUP correction on cosmological scales plays a negligible role and therefore it cannot account for a dark energy fluid at all. This claim can be additionally supported if we consider (30): eliminating the dark energy parameter and solving the normalization condition E(z = 0) = 1 for α Q p /c yields Given the best present estimates of Ω m and Ω r , this condition would be satisfied by values Ω k 0.68. However, such value appears to be incompatible not only with the most stringent available constraint on the spatial curvature of our universe (obtained from Planck, see e.g. [83]), but also with the recent claim, still supported by Planck data, of a positively curved (i.e. Ω k < 0) universe [84,85]. The impossibility of having the GUP corrections account for dark energy is also in agreement with the discussion provided in [28], where it was shown that the GUP effects would amount to a dark energy density ∼ α H 4 . This scaling would cause them to decrease extremely quickly, meaning that the GUP alone would not be able to explain the current acceleration of the Universe at all. Other authors note that, even if the energy density were of the order ∼ H 2 , it would still not be sufficient to account for the present acceleration [86]. A different perspective is provided by [87].
The main goal of our analysis is to find an estimate for the GUP parameter, compatibly with the accuracy of the cosmological probes we used. The result obtained for the quadratic GUP model with the full data set is approximately α 2 Q < 10 59 . This bound is clearly much less stringent than most of those obtained with quantum experiments, but, to the best of our knowledge, it represents one of the tightest cosmological and astrophysical constraints achieved so far, as shown in Table 1. Moreover, it is easy to see from Table 2 that α 2 Q 2 p is sufficiently small, being α 2 Q 2 p < 10 −10 , which ensures the appropriateness of the Taylor expansion in (17).
However, when considering late-time observations only, the picture changes quite considerably. Once again, the GUPinfluenced model is statistically equivalent to the corresponding ΛCDM, which provides further support for the claim that the GUP contributions cannot account for dark energy. Indeed, we have considered late-time observations precisely for this reason: on the one hand, it is well-known that, although the CMB geometrical data are very precise if compared to other probes, they are also intrinsically biased towards a cosmological constant model for dark energy. On the other hand, if any dynamical dark energy behaviour could appear, it would be more manifest for z < 2, a range which is fully covered by the late-time data we employed. In conclusion, the full equivalence of ΛCDM and the quadratic GUP model at late times clearly indicates a negligible role of the GUP on cosmological scales.
The bound on the GUP parameter obtained with late-time data is much more relaxed than that obtained with the full data set: it is α 2 Q < 10 81 , namely one of the least stringent constraint in Table 1 (although it perfectly agrees with that found in [28]). It is also worth noting that α 2 Q 2 p < 10 12 , much larger than that obtained with the full data set as well. This means, a posteriori, that any conclusion drawn from the late-time case should be handled carefully: although we are only able to find an upper bound on α Q p (meaning that part of the parameter space explored by the MCMC would still involve small values of this quantity) it should also be noted that all expressions containing the departure function (17) have been partly evaluated outside their region of validity. This technically represents an extrapolation; therefore, the results are not on a footing as strong as those in the analysis employing the full data set.

Linear GUP
Testing the Friedmann equations modified by the linear GUP has proved way more challenging than in the quadratic GUP case, which also affects the reliability of the conclusions we can draw. We remark that our analysis has focused on the case α L > 0 only, as clarified in Sect. 2.3. First and foremost, it is important to point out the impossibility of fitting early-time data (i.e. CMB and BAO data from BOSS/eBOSS) with a linear GUP model. Performing an extensive exploration of parameter space, we found that the χ 2 could not be lower than ∼ 10 7 , indicating a striking inconsistency between data and model.
Consequently, we tried to check if such inconsistency could be relaxed by including late-time data only, since these naturally allow for more freedom in the dark sector. We were able to obtain a much better fit, but also new problems arose. To begin with, the value of the Bayesian Factor is very small, leading to striking decisive evidence against the linear GUP model.
Concerning the values of the cosmological parameters, we found Ω m ≈ 0.37, which is much larger than expected in the standard lore. The associated errors were quite large as well, but still a ∼ 2σ tension was present with the standard value. Additionally, also the cosmological constant contribution is found to be large, namely Ω Λ ≈ 0.98, which would appear to suggest a strong decelerating effect of the linear GUP contribution.
The constraint on the GUP parameter is α L 10 42 . The weakness of this bound is expected, especially given the lack of precision of late-time cosmological data. However, we remark that all the previous results should be taken with care, because the value of the parameter α p ≈ 10 6 is also quite large and comparable to the case of the quadratic GUP model with late-time data.
To summarize, the previous analysis shows that the linear GUP model strongly disagrees with the presently available cosmological data, which greatly limits the amount of relevant information we can extract from comparison between them.

Conclusions
An abundance of cosmological models often motivated by advances in high-energy physics has been devised in the past decades, including some endowed with a GUP. However, not all of these models have been put through observational tests (especially those dealing with a GUP) even now that the technology to make a wealth of high-precision observations is available. In this paper, we provided an example of testing two such models, whose physical predictions ended up being starkly different.
More specifically, we studied two models arising from Friedmann equations influenced by the GUP, in two formulations involving, respectively, a linear and a quadratic term in momentum uncertainty. We put their cosmological viability to the test by means of a statistical analysis relying on well-known cosmological probes at early and late times. On the one hand, we found that the model endowed with the quadratic GUP mimics ΛCDM very closely when the full set of early and late probes is employed, which means that the GUP contribution cannot play the role of an alternative dark energy fluid. However, the main result of this work is the interesting constraint on the GUP parameter, namely α 2 Q < 10 59 , which, although not as stringent as those coming from quantum experiments, is tighter than most of those obtained from gravitational measurements and is one of the few obtained with cosmological data. On the other hand, the model arising from the linear GUP does not appear to be compatible with the cosmological data we analysed. This shows that not all formulations are equally suitable to find constraints in a specific physical system, since, for example, the linear formulation was found more suitable when using gravitational waves [51].
In the future, it would be compelling to study other cosmological models influenced by other formulations of the GUP. Although the low-energy window for quantum gravity phenomenology, far from Planck scale, is at present only able to constrain the GUP parameter and not measure it exactly, the ever-improving quality of these constraints appears very promising. This work also illustrates the potential held by high-precision cosmological data to constrain the GUP parameter, thus complementing the bounds obtained from gravitational measurements and quantum experiments.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: No new dataset has been created for this work.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .