Gravitational Fluctuations as an Alternative to Inflation III. Numerical Results

Power spectra play an important role in the theory of inflation, and their ability to reproduce current observational data to high accuracy is often considered a triumph of inflation, largely because of a lack of credible alternatives. In previous work we introduced an alternative picture for the cosmological power spectra based on the nonperturbative features of the quantum version of Einstein’s gravity, instead of currently popular inflation models based on scalar fields. The key ingredients in this new picture are the appearance of a nontrivial gravitational vacuum condensate (directly related to the observed cosmological constant), and a calculable renormalization group running of Newton’s G on cosmological scales. More importantly, one notes the absence of any fundamental scalar fields in this approach. Results obtained previously were largely based on a semi-analytical treatment, and thus, while generally transparent in their implementation, often suffered from the limitations of various approximations and simplifying assumptions. In this work, we extend and refine our previous calculations by laying out an updated and extended analysis, which now utilizes a set of suitably modified state-of-the-art numerical programs (ISiTGR, MGCAMB and MGCLASS) developed for observational cosmology. As a result, we are able to remove some of the approximations employed in our previous studies, leading to a number of novel and detailed physical predictions. These should help in potentially distinguishing the vacuum condensate picture of quantum gravity from that of other models such as scalar field inflation. Here, besides the matter power spectrum P m ( k ) , we work out, in detail, predictions for what are referred to as the TT, TE, EE, BB angular spectra, as well as their closely related lensing spectra. However, the current limited precision of observational data today (especially on large angular scales) does not allow us yet to clearly prove or disprove either set of ideas. Nevertheless, by exploring in more details the relationship between gravity and cosmological matter and radiation both analytically and numerically, together with an expected future influx of increasingly accurate observational data, one can hope that the new quantum gravitational picture can be subjected to further stringent tests in the near future.


Introduction
In cosmology, we know that the Universe is not perfectly homogeneous and isotropic, but rather comprises of fluctuations in matter and energy densities. Furthermore, these fluctuations are congregated and correlated in a rather specific manner. Detailed measurements reveals fluctuations of various sizes follows a well-defined patterns, which can be quantified with correlation functions and power spectra [1][2][3][4]. The question of why these density fluctuations are distributed the way they are is thus an important one in cosmology. The conventional explanation for the shape of these power spectra is provided by inflation models, which are based on hypotheses of additional primordial scalar fields called inflatons [5,6,7]. The shape of the observed power spectra are then derived from quantum fluctuation of these primordial inflaton fields, and the agreement of this prediction with observations to high accuracy has been widely regarded as a great triumph and confirmation for inflation [8].
In our previous works [9,10], we have offered an alternative explanation based on gravitational fluctuations alone without inflation, which to our knowledge is the first-of-its-kind. While the theory of quantum gravity remains speculative in the short-distance regime -due to both the infinite number of allowed higher-order operators consistent with general covariance together with a lack of experimental results in this regime, the long-distance or infrared limit of the theory is however in principle well-defined and unique, governed primarily by the concept of universality.
Nevertheless, this long-distance quantum theory of gravity still suffers from being perturbatively nonrenormalizable, rendering perturbation theory useless for calculating any quantum corrections in gravity. However, in the past decades, well known field theory techniques have been extensively developed, applied and even tested to high accuracy in various disciplines of physics where perturbation theory fails (e.g. non-linear sigma model, Heisenberg magnets). It is thus highly conceivable that these nonperturbative techniques may find use in deriving physical consequences for another perturbatively nonrenormalizable theory such as gravity.
From previous efforts [11,12], it was shown that quantum effects of gravity may manifest themselves not only on the extreme small (UV) scales, but also on the extreme large (IR), cosmological scales. In particular, our work [9,10] have shown that, utilizing nonperturbative field theory methods, much of the cosmological matter power spectrum can be derived and reproduced purely from Einstein gravity and standard ΛCDM cosmology alone, without the need of any additional scalar fields as advocated by inflation. We have shown that not only the predictions agree quite well with recent data by the Planck Collaboration [13], but also that additional quantum effects predict subtle deviations from the classical picture, which allows this approach to be testable in the near future with increasingly powerful cosmological experiments.
In this paper, we extended our analysis in two major areas. First, we utilized a number of current numerical cosmological programs, such as ISiTGR, MGCAMB and MGCLASS. Secondly, with the help of these numerical programs, we generated predictions for all other cosmologically significant spectra, including polarizations (C EE l , C BB l , C T E l , etc.) and lensing spectra (C φφ l , C T φ l , etc.). The paper is organized as follows. In Sec. 2, we summarize the theoretical basis as relevant to present discussion. Sec. 3 introduces the numerical programs we use. Sec. 4 presents the numerical results and analysis. Finally, key points and future work are summarized in the conclusion.

Background
In this section, we will provide a brief review of the quantum theory of gravity and how it is related to various power spectra that can be measured in cosmology. More detailed accounts of the nonperturbative approach to quantum gravity and the derivation of the spectra can be found in previous work [9,10,11,12]. The following will therefore only serve to summarize the key points and main results that are relevant for the subsequent discussion.
Quantum gravity, the covariantly quantized theory of a massless spin-two particles, is in principle a unique theory, as shown by Feynman some time ago [14,15], much like Yang-Mills theory and QED are for massless spin-one particles. In the covariant Feynman path integral approach, only two key ingredients are needed to formulate the quantum theory -the gravitational action S [g µν ] and the functional measure over metrics d [g µν ], leading to the generating function where all physical observables could in principle be derived from. For gravity the action is given by the Einstein-Hilbert term appended by a cosmological constant where R is the Ricci scalar, g being the determinant of the metric g µν (x), G Newton's constant, and λ the scaled cosmological constant (where a lower case is used here, as opposed to the more popular upper case in cosmology, so as not to confuse it with the ultraviolet-cutoff in quantum field theories that is commonly associated with Λ). The other key ingredient is the functional measure for the metric field, which in the case of gravity describes an integration over all four metrics, with weighting determined by the celebrated DeWitt form [16]. There are two important subtleties worth noting here. Firstly, in principle, additional higher derivative terms that are consistent with general covariance could be allowed in the action, but nevertheless will only affect physics at very short distances and will not be relevant nor needed here for studying large-distance cosmological effects. Secondly, as in most cases that the Feynman path integral can be written down, from non-relativistic quantum mechanics to field theories, the formal definition of integrals requires the introduction of a lattice, in order to properly account for the known fact that quantum paths are nowhere differentiable. It is therefore a remarkable aspect that the theory, in a nonperturbative context, does not, at least in principle, seem to require any additional extraneous ingredients, besides the standard ones mentioned above, to properly define a quantum theory of gravity.
At the same time, gravity does present some rather difficult and fundamentally inherent challenges, such as its well-known perturbatively nonrenormalizable feature due to a badly divergent series in Newton's constant G, the intensive computational power needed for any numerical calculation due to it being a highly nonlinear theory, the conformal instability which makes the Euclidean path integral potentially divergent, and further genuinely gravitational-specific technical complications such as the fact that physical distances between spacetime points -which depend on the metric which is a quantum entity -fluctuate.
Although these hurdles will ultimately need to be addressed in a complete and satisfactory way, a comprehensive account is of course far beyond the scope of this paper. However, regarding the perturbatively nonrenormalizable nature, some of the most interesting phenomena in physics often stem from non-analytic behavior in the coupling constant and the existence of nontrivial quantum condensates, which are hidden from and impossible to probe within perturbation theory alone. It is therefore possible that certain challenges encountered in the case of gravity are likely the result of inadequate perturbative treatments, and not necessarily a reflection of some fundamentally insurmountable problem with the theory itself. Here, we shall take this as a motivation to utilize the plethora of well-established nonperturbative methods to deal with other quantum field theories where perturbation theory fails, and attempt to derive sensible physical predictions that can hopefully be tested against observations. More detailed accounts on the other various issues associated with the theory of quantum gravity can be found for example in [11,12], and references therein.
For our present discussion, we will mention several main results and ingredients from this perspective. The nonperturbative treatments of quantum gravity via both Wilson's 2 + ǫ double expansion (both in G and the spacetime dimension) and the Regge-Wheeler lattice path integral formulation [17] reveal the existence of a new quantum phase, involving a nontrivial gravitational vacuum condensate [11]. Along with this comes a nonperturbative characteristic correlation length scale, ξ, and a new set of non-trivial scaling exponents, as is common for well-studied perturbatively non-renormalizable theories ν [18][19][20][21][22][23]. Together, these two parameters characterize the quantum corrections to physical observables such as the long-distance behavior of invariant correlation functions, as well as the renormalization group (RG) running of Newton's constant G, which in coordinate space leads to a covariant G( ) with = g µν ∇ µ ∇ ν [12]. In particular, in can be shown [11,24] that for r < ξ, the correlation functions of the Ricci scalar curvatures over large geodesic separation r ≡ |x − y| scales as where d here the dimension of spacetime. Furthermore, the RG running of Newton's constant can be expressed as where m ≡= 1/ξ, as the characteristic mass scale, and 2 c 0 ≈ 16.04 a nonperturbative coefficient, which can be computed from first principles using the Regge-Wheeler lattice formulation of quantum gravity [25][26][27][28][29][30][31].
Here we note the important role played by the quantum parameters ν and ξ. The appearance of a gravitational condensate is viewed as analogous to the (equally nonperturbative) gluon and chiral condensates known to describe the physical vacuum of QCD, so that the genuinely nonperturbative scale ξ is in many ways analogous to the scaling violation parameter ΛM S of QCD. Such a scale cannot be calculated from first principles, but should instead be linked with other length scales in the theory, such as the cosmological constant scale 1/λ. For example, note that the vacuum curvature condensate expectation value can be related to the cosmological constant via the Einstein field equations Thus the quantity ξ can be viewed as parameterizing the gravitational vacuum condensate. The combination most naturally identified with ξ would be such that ξ ∼ 3/λ ≃ 5300 Mpc for the observed value of λ [11,32,33]. The other key quantity, the universal scaling dimension ν, can be evaluated via a number of methods, many of which are summarized in [29][30][31][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50][51]. Multiple avenues point to an indication of ν −1 ≃ 3, which here will serve as a good working value for this parameter; a simple geometric argument suggests ν = 1/(d−1) for spacetime dimension d ≥ 4 [12].
It should be noted that the nonperturbative scale ξ should also act as an infrared (IR) regulator, such that, like in other quantum field theories, expressions in the "infrared" (i.e. as r → ∞, or equivalently k → 0) should be augmented by The aim here is therefore to explore areas where these predictions can be put to a test. The cosmological power spectra, which are closely related to correlation functions, and thus take effects over large distances, provide a great testing ground for these quantum gravity effects.
To make contact with cosmological observations, the gravitational correlation function G R (r) in Eq. (3) has to be related to the cosmologically observed matter density correlation where r = |x − y|, and is the matter density contrast, which measures the fractional overdensity, or fluctuation, of matter denstiy ρ above the average background densityρ. In the literature, this correlation is more often studied in Fourier-, or wavenumber-, space, G ρ (k; t, t ′ ) ≡ δ(k, t) δ(−k, t ′ ) , via a Fourier transform. It is also common to bring these measurements to the same time, say t 0 , so that one can compare density fluctuations of different scales as they are measured and appear today. The resultant object P m (k) is referred to as the matter power spectrum, where δ(k, t) ≡ F (t) ∆(k, t 0 ). The factor F (t) then simply follows the standard GR evolution formulas as governed by the Freidman-Robertson-Walker (FRW) metric. As a result, P m (k) can be related to, and extracted from, the real-space measurements via the inverse transform It is often convenient to parameterize these correlators by a so-called scale-invariant spectrum, which includes an amplitude and a scaling index, conventionally written as It is then straightforward to relate the scaling indices using Eq. (13), giving s = (d − 1) − γ = 3 − γ = 1. Note that G ρ (r; t 0 , t 0 ) is sometimes referred to as ξ(r) in the literature, but we will use the former to avoid confusion with the fundamental gravitational correlation length ξ.
To arrive at a prediction for the matter density fluctuations G ρ , P m from gravitational fluctuations G R , we make use of the Einstein field equations In a matter dominated era, such as the one where galaxies and clusters are formed, the energy momentum tensor follows a perfect pressureless fluid to first approximation. Hence, the trace equation reads (For a perfect fluid the trace gives T = 3p − ρ, and thus T ≃ −ρ for a non-relativistic fluid.) Since λ is a constant, the variations, and hence correlations, are directly related as in δR(x) δR(y) = (8πG) 2 δρ(x) δρ(y) .
As described above, quantum gravity predicts that, over large distances, the scalar curvaturefluctuations scale as G R ≡ δR(x) δR(y) ∼ 1/r 2 . This implies that the matter density fluctuations follow an analogous scaling relation as r → ∞, within the matter dominated era, and thus γ = 2. From the Fourier transform in Eq. (13), we get as k → 0 in wavenumber space, in the matter dominated regime. This result of linear scaling is a well-tested and well-supported result from decades of cosmological measurements of galaxy correlations functions [52].
To extend beyond the linear matter dominated regime, the trace equation alone becomes insufficient (since the trace of the energy momentum tensor for radiation vanishes), and the full tensor equation has to be used. Furthermore, in a real universe with multiple fluid components, interactions and transient behaviors have to be taken into account, which are governed by coupled Boltzmann equations. However, these classical procedures are fully worked out in standard cosmology texts [53,54]. Following [53], the matter power spectrum can be written in two parts -an initial condition known as a primordial spectrum R o k , and an interpolating function between the domains known as a transfer function T (k). Thus, the full P m (k) beyond the galaxy domain will take the form where C 0 ≡ 4(2π) 2 C 2 (Ω Λ /Ω M )/25 Ω 2 M H 4 0 is a constant of cosmological parameters, and the k 4 factor for convenience. The transfer function is usually written in terms of κ ≡ √ 2k/k eq , a scaled dimensionless wavenumber, with k eq being the wavenumber at matter-radiation-equality. With this decomposition, the transfer function is a fully classical solution of the set of Friedmann and Boltzmann equations, capturing the nonlinear dynamics. This leaves the initial primordial function, which can be parameterized as a scale-invariant spectrum which is only parameterized by an amplitude N 2 and a spectral index n s . k R is referred to as the "pivot scale", and is simply a reference scale, conventionally taken to be k R = 0.05 Mpc −1 .
While the transfer function T (κ) -the solution to the highly-coupled and nonlinear set of Friedmann, Boltzmann and continuity differential equations -is difficult to solve, it is in principle fully determined from classical dynamics. Moreover, assuming standard ΛCDM cosmology dynamics and evolution, a semi-analytical interpolating formula for T (κ) [53] is known. As a result, if the initial spectrum R o k , or more specifically the parameters N and n s , is set, then P m (k) is fully determined. To find N and n s , it can be done by matching. Since Eq. (20) is known to be valid in the galaxy and cluster domains, and Eq. (21) is supposed to account for all wavenumber-scales, these equations should overlap in the galaxy domain. So by matching Eq. (20), which is fixed by the scaling of curvature correlation functions, with Eq. (21) in the overlapping region, N and n s can be found, thus fully normalizing P m (k). More precise details of this procedure, as well as detailed comparison plots with the latest observational data, can be found in our previous work [9]. The key resultant analytical prediction for P m (k) from this procedure is also reproduced here in the later plot as the solid blue curve in Fig. 3, showing almost perfect fit to all observational data for Finally for scales of k comparable to 3 √ 2c 0 m, additional quantum effects are expected to become significant, due to the nontrivial vacuum condensation nature of gravity, enough to cause deviations from the classical ΛCDM result of P m (k). This scale is already hinted in for example Eq. (4). These modification can again be done either analytically or relying on a program numerically. Analytically, the effect of the RG running of Newton's constant [Eq. (4)] can be included via dimensional analysis for the correct factors of G to include, and IR regulations rather straightforwardly as per Eq. (8), as is done in other similar quantum condensate theories such as QCD or condensed matter theories. More details can again be found in [10]. These results are reproduced as a plot later ( Fig. 3) to compare with the fully-numerical results, showing great agreement between them. To obtain the latter, i.e. the numerical approach, shall form the focus and the remaining of this paper. Following similar analysis to determining P m (k), other spectra such as the angular temperature spectrum C T T l , should be fully derivable from the primordial function R o k , or specifically n s , which is set by the scaling of gravitational curvature fluctuations ν. In fact, many spectra are only various variations of integral transforms with different physical observable quantities, say, photon temperature and polarization, instead of mass density δρ m . A brief review of that is given in Sec. 4.2. Finally, it should be re-emphasized that in this picture, a scalar field is not an essential ingredient to determine R o k . It should be noted that there are intrinsic uncertainties in some of the theoretical parameters as well. The so-called "analytical" approach as referred to here still relies on either numerical inputs, or some analytic approximations, at earlier different stages in order to extract physical predictions such as the scaling dimensions of G R ν, from Eq. (3), or the amplitude of the first order quantum corrections for the RG running of Newton's G c 0 , from Eq. (4), from the highly nonlinear gravitational path integral of Eq. (1). For example, from latest lattice simulations of the path integral, it is found 1/ν ≃ 2.997(9) and c 0 ≃ 8.02, with the latter an error that is estimated at around 50%. Other methods, summarized in [11], including observational data as studied in [10], all support the value ν ∼ 1/3. As eluded in this paper as well, this is not surprising given the universality nature of this index ν. On the other hand, the amplitude for quantum corrections c 0 , while should remain some order 1 parameter, cannot be claimed to the same degree of confidence as ν. For example, in the comparison with latest observational data for P m (k) in [10], it seems to best fit a value roughly 7 times smaller (c 0 ≈ 1.146). However, it should also be pointed out the theoretical expression defining c 0 possesses a slight degeneracy with the correlation length scale ξ (Eq. (4)). Hence, the data can also be interpreted as suggested a value of ξ ∼ 14000 Mpc, around 2.5 time larger than the expected ξ ∼ 3/λ, or, some combination of both instead. In principle, the inclusion of IR regulation to the final expressions (Eq. (9)) changes the shape of the curve and can in principle break the degeneracy, the current crudeness of the observational data in those regimes of k is much too uncertain to make any conclusions as to the more favorable possibility. While we will continue to primarily refer to studying the constraints on c 0 for simplicity for the rest of this paper, it should be kept in mind the possibility of this degeneracy. It is also hopeful that with increasingly precise observational data in the future, complimented with looking at independent and orthogonal observables that we are to present in this paper, a better constraint on these theoretical parameters can be found.
Finally, it should be noted that the current most popular approach to explain the shape, or more precisely, the index n s of the matter power spectrum is typically reliant on the fluctuations of postulated primordial scalar fields from inflation models [55]. Given the long interest for understanding this spectral index [56,57,58], the ability to derive this index, as well as the lack of competing theories, is thus championed as a triumph of inflation. The picture reviewed here, where the correlations are explained by nonperturbative critical scaling behaviors of gravitational fluctuations, is thus first-of-its-kind. As discussed in this background, the formulation of this picture is in principle rather constricted with little flexibility. As a result, this gravitational picture makes concrete predictions that can be concretely tested (or falsified), without suffering from the typical flexibilities in scalar-field-driven inflation models, and thus offering a compelling alternative to the canonical inflation picture.
Having reviewed this analytic background, we will next present the numerical programs we used, and the subsequent results for the cosmological spectra from effects of quantum gravity.

Numerical Programs
There are a variety of publicly available Boltzmann-Einstein (EB) solvers that have been in use for the past two decades starting with CMBFAST [59]. The main independent programs are CAMB [60] and CLASS [61] which solve the coupled Einstein-Boltzmann equations in a background FRW metric. These codes are computed for ΛCDM cosmology with a limited set of choices for a parameteriazation of equation of state for the Dark Energy (w). In all our programs we use w = −1 ,which considers dark energy as a vacuum energy.
For modifications of gravity with a scale dependent gravitational constant, there are three EB solvers. We used ISiTGR (Integrated Software in Testing General Relativity) [62] as the primary code to generate power spectra. Then we compare with another two programs MGCAMB (Modified Growth with CAMB) [63] and MGCLASS (CLASS version for phenomenological modified gravity) [64]. ISiTGR and MGCAMB are patches for CAMB and COSMOMC [65] which was written in the FORTRAN language, while MGCLASS is a patch for CLASS written in C. All three programs have implemented the parameterization effective gravitational coupling (µ) -gravitation slip parameter (η) which sometimes is denoted as γ. Those two parameters are defined as µ(a, k) ≡ G(a, k)/G 0 and η(a, k) ≡ Φ/Ψ where G 0 is the laboratory value of Newton's gravitational constant and Φ, Ψ are scalar potentials in the conformal Newtonian gauge. The comparison of the three programs for no RG running of G as in standard ΛCDM cosmology is shown in Fig. 1 and Fig. 2.
One can see that while all three program's ΛCDM prediction are generally consistent, only ISiTGR's modified Newton's constant patch with µ(a, k) ≡ G mod /G N ewt = 1 [or equivalently c 0 = 0 in Eq. (9)] is consistent with its original default-ΛCDM prediction. Matter power spectrum from MGCLASS has a noticeable upper trend for small k from the ΛCDM curve, as shown in the left plot in Fig. 1. Fig.2 shows a significant deviation of MGCAMB's C T φ l from the ΛCDM curve.
Primarily due to this reason we chose ISiTGR over these two other programs.
In the ISiTGR program all times are in conformal time, as is the case for CAMB. The growth equations are written based on a perturbed FLRW metric in the Newtonian gauge, Where Φ and Ψ are scalar gravitational potentials , x i represents comoving coordinates and a(τ ) is scale factor at conformal time τ . For a flat universe the three dimensional spatial metric γ ij in cartesian coordinates is given by From now on we only discuss cosmology for a spatially flat universe, to which k = 0 .
There are four built in functional forms for selected modified cosmologies [66] and we used (µ)gravitation slip parameter (η) form. The modified growth equations are  Figure 1: As an example we illustrate the P m (k) predictions between the three programs -ISiTGR (blue), MGCAMB (orange), and MGCLASS (green) -with their corresponding patches for a modified Newton's constant. This serves as a consistency check between the programs and a validity check for their patches. The solid curves are generated from the respective original ΛCDM programs, while the dashed curves are generated by each program's modified Newton's constant patch setting µ(a, k) ≡ G mod /G N ewt = 1. It can be seen that ISiTGR is the most consistent, and hence reliable program of the three, to investigate the effects of a modified Newton's constant. predictions between the three programs -ISiTGR (blue), MGCAMB (orange), and MGCLASS (green) -with their corresponding patches for a modified Newton's constant. This serves as a consistency check between the programs and a validity check for their patches. The solid curves are generated from the respective original ΛCDM programs, while the dashed curves are generated by each program's modified Newton's constant patch setting µ(a, k) ≡ G mod /G N ewt = 1. It can be seen that ISiTGR is the most consistent, and hence reliable program of the three, to investigate the effects of a modified Newton's constant. and where w i and ρ i are respectively the equation of state and density of i th particle species. Generally there are three species which are radiation, non relativistic matter and dark energy. And ∆ i is the gauge-invarient, rest-frame overdensity defined by, where H =ȧ/a is the Hubble's constant in conformal time, fractional overdensity δ i = δρ/ρ and q i is the heat flux, related with the peculiar velocity (θ i ) From the conservation of energy-momentum tensor of the perturbed matter fluids and for uncoupled fluid species ∆ i evolution is given by Secondary effects considered by ISiTGR are reionization,weak gravitational lensing and the ISW effect. For reionization it uses same approach as in CAMB [67], namely a simple tanh model for reionization fraction (x e ), given by where y(z) = (1 + z) 3/2 , z re is the red shift value where the x e = f /2, and ∆ y is the fractional change in y. The latter agrees with a Thompson scattering optical depth for an instantaneous reionzation which occurred at z re . The treatment of weak lensing is discussed here later in section 4.
Since the required formulation for µ(a, k) does not appear as an inbuilt function, we added a part with newly defined functions µ(a, k),μ(a, k) for our need in the above equations. In accordance with Eq. (9) we have andμ(a, k) = 0.
As secondary effects ISiTGR considers reionization,weak gravitational lensing and the ISW effect.η(a, k) = 1 is assumed since there is no different modifications to the potentials. ISiTGR has two binning methods but here we only used the traditional binning method. For all the power spectra computations we set the tensor part to zero. The program computes 2-point self-and cross-correlation functions for the temperature, E-mode and B-mode polarization and weak lensing potential. Each generated power spectrum appears in two separate files, one with lensing and the other without. In the following we use power spectra with gravitational lensing included. The values of the cosmological parameters we used here as initial conditions are shown in Table 1.
In a previous paper [9,10] we used semi-analytic methods to solve for the matter power spectra using semi-numerical approximations for the relevant transfer functions. But in the current approach the numerical programs solve the full set of Boltzmann equations, and uses integration techniques such as adaptive Runge-Kutta method to integrate all the tightly coupled equations.
Secondary effects accounted for like reionization and integrated Sachs-Wolfe (ISW) effect are treated as a more general case compared to our previous work.

Numerical Results
In this section we present numerical results for the quantum gravitational corrections to the various . This includes both the effects of an RG running of Newton's constant and the IR regulation, obtained by replacing For simplicity, of the three numerical programs used for our analysis (ISiTGR, MGCAMB, MG-CLASS), only the results from the ISiTGR numerical program shall be plotted. The reason for this choice is that we expect this program to provide better consistency and reliability in the particular region considered (small k, small l), as explained in Sec. 3. Furthermore, all numerical results presented here are generated using the latest values of the cosmological parameters as given by Planck (2018) [13].
In the above, c 0 is the coefficient that governs the amplitude of quantum corrections. For all the following spectra, three different values of c 0 = 0, 1.146 and 8.02 will be plotted. Lattice calculations give c 0 ≈ 8. However, being a non-universal parameter, it can depend on specific choices arising from the way an ultraviolet cutoff is imposed. Therefore, not too much weight should not be placed on this specific value, beyond perhaps the order of magnitude. In practice, this value could be further constrained by experiments, which is precisely what these observations of cosmological spectra can achieve. From previous work [10], using the approximate semi-analytical methods, we see that a value of 8.02/7 ≃ 1.146 is generally favored.

Matter Power Spectrum P m (k)
We start with the matter power spectrum P m (k). Recall the definitions where the variable δ ρ ≡ (ρ −ρ)/ρ is the fractional density fluctuations above the average, referred to in cosmology as the mass-density contrast. The numerical results for P m (k) obtained from the numerical program (ISiTGR), for both the classical ΛCDM (i.e. c 0 = 0) and quantum (c 0 > 0) results, as well as the respective analytical results (as derived in [9,10]), are shown and compared in Fig. 3.
From Fig. 3, we see that all the numerical results are generally consistent with the corresponding analytical results from earlier work, which were obtained by following the semi-analytical interpolating formulas from [53], and the implementation of the RG running following dimensional arguments. The small deviations may be attributed to the slightly older values of cosmological parameters [68] and some analytic approximations used by Weinberg and Dicus' interpolating formula for the transfer function in [53], whereas the numerical results presented here use the latest cosmological parameter values from the Planck collaboration [13]. Despite the small discrepancies, we see that the overall trends, and the extra downwards bend due to the inclusion of the (IR regulated) RG running of Newton's constant, as predicted analytically using the semi-analytical formulas are in very good agreement with the numerical predictions using the latest fitted cosmological parameters. This overall general agreement between the analytical and numerical result  9)], of which their detailed derivations can be found in [9,10]. The corresponding dashed curves represent the corresponding the numerical predictions generated by ISiTGR, showing very good general consistency with the trend derived from analytical methods. The observational CMB and galaxy data points, taken from the Planck (2018) collaboration [13] and Sloan Digital Sky Survey (SDSS)'s 14th Data Release (DR14) [52] are also shown.
provides a good verification and confidence that the procedure of including a running Newton's constant as presented above is reliable.
The same numerical analysis has now been repeated with the other two numerical programs MGCAMB and MGCLASS, besides ISiTGR. The result of MGCAMB is in extremely good agreement with ISiTGR, with its predictions for all three values of c 0 almost completely overlapping with ISiTGR's result, giving additional confidence to the latter. However, while MGCLASS is relatively consistent with ISiTGR for most of the angular spectrum results (as we will discuss later), its prediction for P m (k) shows a rather radical upturn below k = 10 −3 , which is at odds with both ISiTGR and MGCAMB, as well as the analytical predictions (also shown and discussed earlier in Fig. 1 and then in Sec Given the more confident, and in principle more accurate, predictions from the numerical programs as shown in Fig. 3, it can be seen that the value of c 0 ≃ 1.146 is a better overall fit to the observational data from Planck, which is a consistent conclusion from our previous work that was based exclusively on the early analytical results. Armed with the new tools of numerical programs, we will now move on to present the numerical results for the other various correlation functions, which will hopefully shed new insights to the validity of the quantum gravity effects in cosmology. It should be noted that there is a slight degeneracy between c 0 and ξ in the original expression for the RG running of Newton's G, as in Eq. (9). A support of a smaller c 0 from the observational data can equivalently be mimicked by an increase in the vacuum condensate scale ξ. In fact, the apparently better fit value of c 0 ≃ 1.146, seven times smaller than the lattice predicted value of ≈ 8.0 ± 3.1, can be mimicked by simply a factor of ∼ 1.9 larger in ξ. Technically, including IR regulation will change the shape of the curve and break the degeneracy, which in principle could be fitted sophisticated say with a Monte Carlo simulation. However, not only is that currently beyond the scope of this paper, the lack of and crudeness of data points in the small-k regime will not render the exercise fruitful.
On the other hand, it may be instructive, amongst other physical motivations, to look at the quantum effects on a variety of other spectra of cosmological significance with these numerical programs. With independent quantities and measurements, the new plots may either provide additional constraints to these quantum gravitational parameters, but also potential insights to the physics.

Angular Temperature Power Spectrum C T T l
The TT power spectrum is one of the most important cosmological spectrum since it is measured to high degree of accuracy, thus allowing for great insights in constraining various cosmological models. Fig. 4 shows the numerical predictions for the temperature-temperature (TT) angular power spectrum C T T l with and without the quantum effects. We will first briefly recall the definitions for C T T l and how theoretical predictions for it can be obtained, and then compare them against observational data. Following notations in Weinberg [53], the temperature fluctuations ∆T can first be resolved into spherical harmonics Y m l (n)'s, where T (n) is the temperature in the directionn, T 0 ≡ (1/4π) d 2n T (n) the average temperature over the sky, and the coefficients a T lm quantifying the fluctuation for each harmonic. Since ∆T 's are real, and the products of ∆T 's are rotationally invariant, one has Here the L l are the Legendre polynomials, and C T T l is defined as the 2-point correlation functions of a T lm , the temperature fluctuation in "l"-space. Or equivalently, by inverting the transformation. As a result, the correlations for temperature-temperature fluctuations are fully quantified with the C T T l 's. (Note that here we use L l instead of the usual notation P l for the Legendre polynomials, in order to avoid confusion with the matter power spectra.) Theoretically, since CMB photon temperatures and matter density are coupled, the C T T l 's are therefore related to the matter power spectrum P m (k), via integral transforms that involve spherical Bessel functions and appropriate form factors and transfer functions. However, from transforming the predictions from one set of observable to another, new insights, and potential constraints to the theory, can be derived.
To do so, one can first use the Friedmann and continuity equations to relate the temperature fluctuations to the metric perturbations, via suitable form factors F 1,2 (q), through where the latter are defined as Here the B and E functions are suitable decompositions of the metric perturbations, and δu γ is the velocity potential for the CMB photons. It is known that these form factors simplify in certain gauge choices. In the synchronous gauge, one has E = 0, whereas in the Newtonian gauge B = 0 and E = 2Φ, which then gives (Note that F 1 (q) and F 2 (q) are referred to as "F (q)" and "G(q)" respectively in [53]. Here we will use the former in order to avoid confusion with the expression for the running of Newton's constant G(k), as it will be implemented below. The above equations also assumed a sudden transition to opacity on the CMB at a time t L , which nevertheless does not change the form of the basic equations and only some of the details, which are later taken into account fully with the numerical programs, discussed below.) Hence, given appropriate initial conditions, the functions Φ and δu γ , as well as the scale factor a(t) and the function T (t), can all be obtained as solutions of the classical Friedmann equations.
These are then combined with the Boltzmann transport equations, as is done in standard cosmology, which eventually leads to unambiguous predictions for the C l 's. The solutions for F 1,2 (q) can be parameterized in terms of transfer functions T (κ), S(κ) and ∆(κ), leading to the following expressions for F 1 (q) and F 2 (q) spectrum R o q , where q is the wavenumber, and "o" refers to outside the horizon. Conventionally, R o q is parameterized by an amplitude N and spectral index n s , Here the reference "pivot scale" is usually taken to be q R = 0.05 Mpc −1 by convention. As a consequence, once the primary function R o q is somehow determined, classical cosmology is then expected to fully determine the form of the C T T l spectral coefficients. It is therefore possible to write the C T T l 's fully, and explicitly, in terms of the primary function R o q . After expanding the plane waves factor in a complete set of spherical harmonics and spherical Bessel functions, C T T l from Eq. (38) becomes where r L = r (t L ), and we have factored out the function R o q explicitly by defining F 1 (q) = (R o q ) F 1 (q) and F 2 (q) = (R o q ) F 2 (q). Now, recall that the matter power spectrum P m (k) is given by which tells us that we can obtain a direct relation between the matter power spectrum P m (k) and the angular temperature coefficients C T T l , where q and k are related by q = a 0 k, and the scale factor "today" a 0 can be taken to be 1. As a result, the predictions on P m (k) can be directly transformed into a prediction for C T T l . Utilizing the same parameters in the numerical programs, the effects of with and without the RG running of Newton's constant (with IR regulation) on C T T l can then be generated. Fig. 4 shows the numerical result of c 0 = 0 (no running) 1.146, and 8.02 with the blue, green, and orange curve respectively, generated by ISitGR. The observational CMB data from Planck (2018), as well as an (error-weighted) cubic fit for ence, is also shown. Noticing that the point l = 2 is anomalously low, with large uncertainty due to cosmic variance, the error-weighted fit shown in this plot has not included the l = 2 point.
From Fig. 4, we see that the effects of a RG running of Newton's constant generally cause an upturn to the spectrum at low-l's, starting at roughly l = 20. It can also be seen that the orange curve with a quantum amplitude [see Eq. (9)] c 0 = 8.02 (or ξ = 5300 Mpc) creates a much more dramatic deviation, reaching a maximum of 140% larger in value compared to the blue, classical (c 0 = 0) ΛCDM curve, while the green curve with c 0 = 1.146 (or roughly ξ ≃ 2.65 × 5300 = 14000 Mpc) has a milder deviation of ≈ 24% from the classical result. Again, neglecting the anomalous l = 2 point, the green curve with c 0 = 1.146 is generally consistent with all observational data, arguably also with the desirable feature of marginally going through the error bars of l = 5 and l = 6. On the other hand, the orange c 0 = 8.02 curve, while still lying within a few points' error margins, is less favorably supported by the data in this plot. It is also seen that its deviations starts earlier at a higher value around l ∼ 22, which causes it to miss a few more error bars in the low l points. As a result, the numerical results of this T T plot shows that the green c 0 = 1.146 (or ξ ≈ 14000 Mpc) curve is currently a more favorable parameter than the orange one. Note that this is also consistent with the discussion and conclusion from the matter power spectrum P m (k) plot in Fig. 3.
We also investigated the results with all 3 programs. However unlike P m (k), the three programs do not agree, despite being supplied with the same RG modified expression for Newton's constant. curve coinciding with ISiTGR's c 0 = 1.146 curve. I.e. MGCAMB seem to predict an upturn around 7 times larger than ISiTGR, while MGCLASS seem to predict an upturn that is 7 times smaller than ISiTGR. Given the blackbox nature of such programs, it is unclear of the cause of this different given that all programs where supplied the same modification in Newton's G. These programs, designed for modified gravity models, are known to be less well-tested compared to their base program (CAMB, CLASS), and it may not be surprising that two (or all) of them may be incorrect. One consistency is that all three programs predicts an upturn at low l's, just to a different degree, roughly ±1 order of magnitude. Hence, it is at best that we can conclude from these available programs that the RG running of Newton's G causes an upturn to roughly the order of magnitude presented in Fig. 4.
Perhaps even more intriguing is the disagreement with a naive analytical analysis. From Eq. (49), the first order estimate is that since C T T l is the (weighted) integral of P m (k) over all k, a smaller P m (k) caused by an RG running (c.f. Fig. 3) should cause a smaller value of C T T l . In fact, if one assumes the transfer functions are not affected by the quantum corrections, the integral Eq. (49) can be performed numerically (as done in [10]), since the classical interpolating formulas for the transfer functions are known, which does show a downturn, as naively expected, instead of an upturn. This work utilizes programs that in principle modifies the initial Friedmann and Boltzmann equations from the beginning, and includes any effects of the RG modified Newton's G into the solutions, and thus in principle more trustworthy. But given the opaque nature of such programs, it remains further investigations through a more detailed study of the entangled initial set of coupled differential equations to fully understand the disagreements between the programs and the first-order analytical result, as well as the disagreement, and hence the reliability, within the numerical programs.
Nevertheless, given that these programs represents the most sophisticated tools currently, it is still constructive to look at their predictions of the quantum effects on other modes and variables of the CMB. For example, the theoretical predictions for the percentage deviations for c 0 = 1.146 curve with the classical curve is ∼ 37% on P m (k) at its further available data point, while only ∼ 24% on C T T l . This reveals the fact that the quantum effects maybe more significant in different physical variables. So by studying the predictions for different auto-and cross-correlations of different varaiables, and compare them to potentially independent data (e.g. ground-based measurements of E-and B-mode polarizations as oppose to space-based measurements of CMB temperature), new constraints and insight may be deduced. We will present the analysis and results of the other spectra of interest to cosmology in the remaining of this section. Following notations in [53], CMB photons distributions are fully described through a number density matrix n ij (x, p, t), or, more usefully, the dimensionless version of its perturbation J ij (x,p, t) (referred to as the dimensionless photon intensity perturbation matrix), related to n ij via

Temperature-E-mode Power Spectrum
In a line-of-sight directionn, J ij can be parameterized via where Q, U , and V are three real functions of direction (with units of temperature), known as the Stokes parameters, describing the photon's polarizations. Notice that the photon temperature perturbations are given by the trace It is these Stokes parameters that are measured in current observations of the CMB. But since the scattering of light by non-relativistic electrons does not produce circular polarization, one expects that all CMB photons will be linearly polarized, so that J ij is real, and therefore V = 0. For further convenience of in comparing with observations of 2-point functions, which respect spherical symmetry, it is useful to expand the Stokes parameters Q(n) and U (n) seen in a directionn in a series of functions Y m l (n) where the subscript "P " in the coefficient a P,lm stands for "polarization",∇ is the angular part of the gradient operator, and e ± (n) = (1, ±i, 0)/ √ 2 are the polarization vectors in the directionn.
To further satisfy the reality condition, one defines the amplitudes a E,lm ≡ − a P,lm + a * P,lm /2 , a B,lm ≡ i a P,lm − a * P,lm /2 , so that their correlation functions  respectively. Here again one finds higher trends at large angular scales (l < 10) as compared to the classical (no running) numerical ΛCDM curve. The dashed curve represents an error-weighted cubic fit to the observational CMB data from Planck (2018) [13]. polarization respectively, since spatial-parity inversion, a E,lm → (−1) l a * E,lm , and similarly for a T,lm , whereas a B,lm → −(−1) l a * B,lm . As a result of parity, there are no bilinear correlations between B with either E or T . (i.e. C T B l = C EB l = 0.) With this background, we shall present the numerical predictions for the corresponding spectra with and without an RG running of Newton's G, compared against the latest observational data.
We start with the T E spectrum. Fig. 5 shows the numerical results with the observational data for C T E l . The lowest solid (blue) curve represents the classical (c 0 = 0) spectrum, while the solid middle (green) and top (orange) curve represents the effect of a RG running Newton's constant with c 0 = 1.146 and 8.02 respectively.
It turns out new constraints for the RG running parameter c 0 can be deduced with this new plot. With the inclusion of the E-type polarization data, we see that this has further constraints on some of the error bars in the low-l data points. This is due to the smaller error bars from the observational data in the E-type polarization correlations in the low-l regime (see Fig. 6). As a result, one sees that the top c 0 = 8.02 curve (orange) is strongly disfavored by this plot. Another observation is that the difference between the c 0 = 1.146 and the classical ΛCDM (c 0 = 0) curve is about 60% in this T E plot, which is a larger percentage deviation compared to ≈ 24% for the T T plot.
We also compared the results from the other two programs (MGCLASS and MGCAMB, not shown on Fig. 5). All the resultant curves of MGCLASS agree with ISiTGR for l > 30, but for l < 30, the c 0 = 8.02 curve of MGCLASS is about 36% lower than the corresponding ISiTGR curve. All two curves with RG running from MGCLASS are within the error bars but due to the mismatch as shown in the Fig.1, MGCLASS results should be investigated further. For MGCAMB, the curves with RG running are significantly higher than ISiTGR, making them disfavored. Also there is a slight horizontal shift for MGCAMB in l-space compared to the other two program, which should be investigated further.

EE-Power Spectrum C EE l
We move on to the EE spectrum. Fig. 6 shows the numerical results with the observational data for C EE l . We also plotted an error-weighted cubic fit (dashed line) for the classical (c 0 = 0) spectrum (solid blue), as well as the quantum RG running of G for the above values for c 0 (green and orange). It can be seen that there is no significant deviation from standard ΛCDM prediction like in temperature power spectra and all the curves are well within the data point error bars. We can see that in the large scales (l < 20) the errors are significantly small which makes T-E spectrum . The solid curves represent the numerical predictions generated by ISiTGR, with the bottom (blue), middle (green), and top (orange) representing quantum amplitudes quantum amplitudes [see Eq. (9)] c 0 = 0, 1.146, 8.02 respectively, showing slightly higher trends for large angular scales (l < 10) as compared to the classical (no running) numerical ΛCDM curve. The dashed curve represents error-weighted cubic fit for observational CMB data from the Planck (2018) collaboration [13].
having smaller error bars in the scale of interest in this paper.
When the other two programs are compared, there is no significant deviation to rule out any any curve. There is no noticeable deviation for MGCAMB curves from ISiTGR for l < 20 but there is a slight upward deviation for l > 20. With MGCLASS, the RG curves are smaller than ISiTGR making smaller deviation from ΛCDM curve. The theory of CMB lensing is a vast topic on its own. Here, we will try to present the key defining equations of the lensing spectrum, and then look at the numerical results of quantum gravitational effects on the lensing potential spectra. A more complete account for the physics and observations can be found in [53,69,70].

BB-Power Spectrum
Consider a small deflection angle θ from the undeflected directionn of a CMB photon, with θ describing perpendicular direction ton, and |θ| ≪ 1. Define the shear matrix M ab as where a, b run over the directions orthogonal ton, r S is the radial distance of the source from earth in a Robertson-Walker coordinate system, and δθ a is the amount of deflection of θ. From standard general relativistic calculations, the shear matrix can be related to the Newtonian potential of the lens source (φ), via where y is a small perpendicular deflection vector ton, and t r is the time for a photon that that just reached us from radial coordinate r. Hence, the measurements of the shear matrix can yield . The solid curves represent the numerical predictions generated by ISiTGR, with the top (blue) and bottom (green) representing c 0 = 0, 1.146 respectively, showing no significant deviation with the classical no running numerical ΛCDM curve. The dashed curve represents error-weighted quadratic fit for observational CMB data from the Planck (2018) collaboration [13]. The other c 0 = 8.02 curve, like the c 0 = 8.02 (green) curve, is consistent with zero (to around 1 part in 10,000). So any deviations from the classical and c 0 = 1.146 curve are too insignificant to be seen, and negligible relative to the size of the error bars from the current latest data. Hence the c 0 = 8.02 curve is not included in this plot for clarity.
information about perturbations to the gravitational potential (δφ) by masses spread along the line of sight. Define the so-called lensing convergence field κ as κ is particularly useful because, if the lensing is due to a collection of bodies all at about the same radial coordinate r L , it can be directly related to the matter perturbations δρ m . More explicitly, δφ falls off rapidly for large distances, so that the factor r(r, r S )r can be replaced in a first approximation with r(r L , r S )r L , and similarly the second term with ∂ 2 /∂r 2 can be dropped.
resulting in an expression directly linking κ to matter density fluctuations δρ m . Hence, a measurement of the value of κ for sources seen in one direction can reveal the total mass density of a cluster of lensing masses that lies along that line of sight at distance r L (projected onto a plane perpendicular to the line of sight). Since, as we have shown, gravity constraints the scaling of correlations of matter, it should also do so for κ.
So, to project the convergence field κ onto the sky, we decompose it in a way that is analogous to the other angular spectra, with with quantum noise fluctuation correlation Or more explicitly, by inverting the expression in Eq. (67), In the literature [69,70], it is often the correlation for the lensing potential C φφ l that is plotted, instead of that of the lensing convergence field C κκ l , related by with respective expansion coefficients a T,lm and a E,lm , similar to Eq. (56)- (59). With this background, we will present the numerical results of including a quantum RG running of Newton's constant for these spectra.
In Fig.8 for C φφ l we have plotted for the classical (c 0 = 0) spectrum (solid blue), with as well as for an RG running of G with the previously used values for c 0 (green and orange). One can see that there are a significant deviation which is up to 80% for c 0 = 8.02, but only 10% for In Fig. 10 we show the results for C Eφ l , and we have plotted the classical ΛCDM (c 0 = 0) spectrum (solid blue) compared with the RG running of Newton's G spectrum with the above values for c 0 (green and orange). It can be seen that there is significant deviation for c 0 = 8.02 which drops by more than 50% which gets ruled out. The deviation begins scales corresponding to l < 30 . For

Conclusion
In this paper, we have revisited the derivation of the matter and temperature power spectra from the quantum theory of gravity without invoking any additional scalar fields from inflation, which, to our knowledge, is the first of its kind. We reviewed that while the short-distance quantum theory of gravity remains speculative, the long-distance behaviors are well known and primarily governed by the renormalization group (RG) behaviors near its critical point. In particular, we reviewed how the critical scaling dimension "s" of the correlation function of the scalar curvature fluctuations at large distances directly governs the scalar spectral index "n s " of the cosmological spectra, as well as the additional quantum gravitational effects, such as the (IR-regulated) renormalization group running of the coupling constant (Newton's constant) G, that will affect these spectra subtly at large distances. We then presented the various numerical programs that we used in this work, and their main results, to complement the previous mainly analytical analysis. We then utilized these programs to further study other cosmological spectra of different modes. We compared these with latest available observational data, and provided new constraints and insights to the parameters (c 0 , ξ) of the quantum theory. We also discussed the possibility of verifying, or falsifying, some of these hypothesis with increasingly powerful observational cosmology experiments in the future.
Using the numerical results, we find that especially the plots of the matter power spectrum P m (k), the angular temperature spectrum C T T l , and the angular temperature-E-mode spectrum C T E l -all play an important role in revealing new insight to constraining the quantum amplitude c 0 , a parameter that governs the size of quantum corrections due to the RG running of Newton's constant. We find that all three plots agreeably favors a value of c 0 closer to around 1. and C Eφ l , are potentially feasible candidates in providing further insights and constraints. Especially for the T φ plot, showing around 20% and almost 150% deviation for the c 0 = 1.146 and c 0 = 8.02 curve respectively, from the classical curve. However, all these latter spectra suffer from a lack of observational data in the low-l regime, making it impossible to draw any conclusion about the favorability of the parameter or the RG running in general at this stage.
However, although the percentages differences between the spectra with and without quantum corrections are decently significant for scales below l < 10 -ranging from ∼ 15 − 60% even with the milder value of 1.146 for c 0 , the uncertainties from current observational data in those ranges are unfortunately even larger. As a result, it is not yet possible to conclude at this stage the visibility of these effects. At best, one can claim the slight hints of RG running from the smallest data point in P m (k), as well as the last few points (3 < l < 7, ignoring the anomalous l = 2 point) of C T T l . Nevertheless, with technology and precision of cosmological experiments improving at a rapid pace, better observational data in this regime perhaps forms one of the most promising area where quantum effects of gravity can be revealed and tested for the first time. This is a consequence of the concrete predictions of the long-distance quantum effects, based on well-established renormalization group analysis, as opposed to the still rather speculative short-distance theories of gravity.
From a theoretical perspective, the numerical results from this work also serve an important purpose in ruling out the less favorable value of c 0 = 8.0 for the quantum amplitude, but instead suggesting a value around seven times smaller, closer to c 0 = 1.15. We also noted that the uncertainties in the observational data at low-l's cannot yet fully constrain the precise shape of the RG running, allowing for the possibility that these various deviations can all be mimicked instead by a modified value of ξ ≈ 14000 Mpc, or around 2.5 times larger than the naive estimate ξ ≃ 3/λ = 5300 Mpc. As we discussed in the theory section (Sec. 2), unlike the universal critical scaling index ν (shown from various method to have a value very closed to ν = 1/3), the parameters c 0 and ξ do not necessarily follow from universality, but are instead confident only up to order of magnitudes. While the observational data at this stage cannot yet exhibit the effects of RG running, they do provide a useful constraint to the possible values of these theoretical parameters. In particular, as studied in detail in our earlier work [10], even with the current observational data's precision, they provide an extremely stringent constraint on the allowed values of ν, down to at most a 1 − 2% deviation from 1/3. This result not only provides a great verification of the values obtained from various theoretical methods such as the Regge lattice calculations of the path integral, but perhaps the first phenomenological test of the quantum theory of gravity in cosmology. It is thus hopeful that as observational technology continues to improve, more insights can be gained regarding the values for c 0 and ξ. With more data and smaller error bars, one can further narrow down a best fit value for the quantum amplitude c 0 or vacuum condensate scale ξ by Markov Chain Monte Carlo (MCMC) sampling in the ISiTGR program. In addition, ISiTGR is also capable of calculating tensor perturbations, which can be used to test this quantum gravitational picture as soon as more observational data on that becomes available. As a fundamentally tensor theory, this gravitational fluctuation picture is expected to produce nontrivial predictions to those of scalar field based inflation models.
It should also be noted that the numerical programs show a very encouraging agreement with the analytical results on the matter power spectrum P m (k), as shown here in Fig. 3. This agreement provides great confidence in the analytical methodology used in [9], or as summarized here in Sec.
2. The concordance between numerical and analytical results provides extra support on how the quantum fluctuations of the gravitational field are linked to the fluctuations of the matter density field. However, the numerical results for the effects of a RG running of G, suggesting an upturn at low l's, seem to disagree with the analytical intuition that a lower P m (k) should give a lower C T T l , as suggested in Eq. (49). Since the derivation of Eq. (49) is purely classical and does not involve any quantum gravitation input, this suggests a lack of analytical understanding of the effects of a having a modified RG running Newton's constant on the Boltzmann equations, and thus their solutions of the form factors F 1 and F 2 [Eqs. (44), (45)]. It is unclear analytically from the coupled differential equations how the running of Newton's G from Eq. (33) affects their solutions, making it difficult to translate the predictions on P m (k), which agrees with the numerical results, to C T T l . This is an area under active further theoretical investigations, and will be addressed in future work. Nevertheless, armed with the supposedly more comprehensive and reliable numerical programs, new insights should be gained regarding the various quantum effects of gravity on the different cosmological spectra.
In conclusion, we have presented in this paper a compelling alternative picture for the various observed cosmological spectra that is motivated by gravitational fluctuations. In this work, we provided updated and extended analysis utilizing numerical programs in cosmology, as well as new physical predictions that can potentially distinguish this perspective from that of standard scalar field inflation. While inflation still currently forms one of the more popular approach, its full acceptance has remained controversial [75,76]. While there exists a number of alternatives to the standard horizon and flatness problems [77,78], the ability to explain the various cosmological power spectra has long been one of the unique predictions from inflation-motivated models, and thus often considered as one of the "major successes" for inflation. It is thus significant that this work provides an alternative, which is in principle arguably more elegant as it only uses Einstein gravity and standard nonperturbative quantum field theory methods, without the usual burden of flexibilities of inflation. Nevertheless, because of the limited precision of current observational data, it is not yet possible to clearly prove or disprove either idea. Still, the possibility of an alternative explanation without invoking the scalar fields is significant, as it suggests that the observed power spectra may not be a direct consequence nor a solid confirmation of inflation, as some literature may suggest. By exploring in more details the relationship between gravity and cosmological matter and radiation both analytically and numerically, together with the influx of new and increasingly accurate observational data, one can hope that this hypothesis can be subjected to further stringent tests in the future.