Gravitational Fluctuations as an Alternative to Inflation

The ability to reproduce the observed matter power spectrum $P(k)$ to high accuracy is often considered as a triumph of inflation. In this work, we explore an alternative explanation for the power spectrum based on nonperturbative quantum field-theoretical methods applied to Einstein's gravity, instead of ones based on inflation models. In particular the power spectral index, which governs the slope on the $P(k)$ graph, can be related to critical scaling exponents derived from the Wilson renormalization group analysis. We find that the derived value fits favorably with the Sloan Digital Sky Survey telescope data. We then make use of the transfer functions, based only on the Boltzmann equations which describe states out of equilibrium, and Einstein's General Relativity, to extrapolate the power spectrum to the Cosmic Microwave Background (CMB) regime. We observe that the results fit rather well with current data. Our approach contrasts with the conventional explanation which uses inflation to generate the scale invariant Harrison-Zel'dovich spectrum on CMB scales, and uses the transfer function to extrapolate it to galaxy regime. The results we present here only assume quantum field theory and Einstein's Gravity, and hence provide a competing explanation of the power spectrum, without relying on the assumptions usually associated with inflationary models. At the end, we also outline several testable predictions in this picture that deviate from the conventional picture of inflation, and which hopefully will become verifiable in the near future with increasingly accurate measurements.


Introduction
In cosmology, we know that matter in the universe is not homogeneously and isotropically distributed. Nevertheless, these inhomogeneities are far from random, but congregated in a rather specific manner, which can be characterized by almost scale invariant correlation functions. The questions of how and why matter is distributed the way it is are thus important ones in cosmology.
The origin and distribution of such fluctuations in matter density are the subject of this study.
According to current established physical laws, two main aspects determine the cosmological evolution. The first aspect involves the use of the correct field equations for a coupled matter gravity system, and these in turn follow from Einstein's classical field equations, coupled with an exhaustive specification of all the various matter and radiation components and their mutual interactions.
The second main ingredient is a set of suitable initial conditions for all the fields in question, which should then lead to a hopefully unambiguous determination of the complete subsequent time evolution.
While there is currently very little disagreement on the nature of the correct cosmological evolution equations themselves, which largely follow from the choice of a suitable metric based on physical symmetry arguments, and on a complete list of matter and radiation constituents based on our current understanding of fundamental particle physics, the same cannot be said for the choice of initial conditions. The latter are largely unknown, and generally involve a number of explicitly stated, and sometimes implicitly assumed, assumptions about what the universe might have looked like close to the initial singularity. Furthermore, it is generally expected that quantum effects will play a major role at such early stages, be it for the matter fields and their interactions, or for the gravitational field itself, for which a classical description is clearly inadequate in this regime. Indeed, over the years attempts have been made to partially include some quantum effects, by assuming non-trivial Gaussian (free field) correlations for some matter and gravity two point functions.
Recent attempts at overcoming our admitted ignorance regarding the initial conditions for the evolution equations have followed a number of different avenues. In one popular scenario [1,2,3], it is argued that the universe might not have evolved initially from a state that could not be described in any way in terms of a simple physical characterization. If, as suggested by particle physics and string theory, a modified quantum dynamics becomes operative at short distance, then one would expect a complete removal of the initial spacetime singularity replaced here instead by some sort of bounce. One key, but nevertheless natural, consequence of this perspective is that the universe must have evolved into and out of the initial singularity in a highly coherent quantum state, with non-trivial quantum correlations arising between all fields, with the latter presumably operating on all distance scales. These highly coherent and complex initial conditions would then represent the surviving physical remnants of a previous, presumably very long, cycle of cosmological evolution, thus generating a set of suitable initial conditions for the cosmology of our current universe.
Another possibility regarding the initial conditions is inspired by the semi-classical approach to quantum gravity [4,5] and follows from a more geometric point of view regarding gravity. In this approach the universe naturally evolves instead from an initial simple, symmetric and elegant geometric construct, as described in practice for example by the so-called no-boundary proposal for homogeneous isotropic closed universes with a cosmological constant and a variety of other fields.
Another popular attempt to explain features of the early universe, and more specifically properties of the matter power spectrum, is through inflation models [6,7,8]. These propose that structures visible in the Universe today get formed from quantum fluctuations in a hypothetical primordial scalar inflaton field. Although the precise behavior and dynamics of the inflaton field are to this day still largely controversial [9,10,11,12,13,14], one nevertheless finds some general features that are common among these models, and which allow one to derive predictions about the nature of the power spectrum.
In broad terms, the approach followed in this paper can be viewed as more in line with the first scenario, where as little as possible is assumed about the nature of the initial quantum state of the universe, given that, as mentioned previously, the latter might be quite far from a simple and symmetric physical characterization. Instead, in this paper, we make extensive use of previously gained knowledge from non-perturbative studies of quantum gravity regarding the large distance behavior of gravitational and matter two-point functions. An important ingredient in the results presented below will be therefore the non-trivial scaling dimensions obtained from these studies, and how these relate to current observational data. We note here that the quantum field theoretic treatment of perturbatively non-renormalizable theories, and the determination of their generally non-trivial scaling dimensions, has a rather distinguished history, originally developed in the context of scalar field theories by Wilson and Parisi [15,16,17], and recently reviewed in great detail for example in [18]. Attempts at quantizing gravity also have a long and distinguished history, dating back to Feynman's original investigations and the formulation of a covariant path integral framework [19,20]. Quantum gravity, and by now the rather compelling evidence of a nontrivial ultraviolet (UV) renormalization group fixed point in four spacetime dimensions, in principle leads to a number of unambiguous predictions, recently reviewed and summarized in [21,22]. The most salient observational effects of such an approach include a running of Newton's constant G with scale on very large cosmological distance [23], the modification of classical results for the growth of relativistic matter density perturbations and their associated growth exponents [24], and a non-vanishing so-called slip function in the conformal Newtonian gauge [25]. The existence of a non-trivial UV fixed point, arising from the highly non-linear nature of the Einstein-Hilbert action, leads in a natural way to nontrivial quantum condensates in the curvature. Such a condensate effectively produces long-distance correlations between curvature fluctuations and, through the field equations, these couple to matter density fluctuations. The fluctuations in matter density, which are in principle unambiguously calculable in this perspective, should lead to gravitational clumping of matter in the over-dense areas, which eventually leads to the formation of galaxies. In other words, the correlations of galaxy distributions should in principle become predictable from correlations of curvature fluctuations derived from quantum gravity. The calculated results can then be compared with observations, and thus viewed as a potential test for the vacuum condensate picture of quantum gravity.
The paper is organized as follows. In Section 2 we outline the canonical ways to parameterize matter distribution quantitatively in cosmology, via n-point correlation functions and the power spectrum, which in turn are largely characterized by the so-called spectral indices. After introducing and defining the correlation functions, we end the section by quoting observed values of the spectral indices from modern measurements. Subsequently, in Section 3 we provide an explanation for the origin of matter fluctuations as a result of underlying quantum gravitational fluctuations for which the correlation functions, and hence spectral indices, can then be unambiguously calculated. We will then show that the theoretical predictions are in good agreement with the observed results. Section 4 discusses how this theoretical prediction can be extended via the so-called transfer function to the earlier time of radiation-dominated epoch, which again shows a generically good agreement with current CMB and other measurements. Section 5 outlines the alternative, more conventional, explanation of the same spectrum via inflation, and briefly discusses the differences between the classical and the quantum gravity picture. Section 6 elaborates on further physical effects predicted by the quantum theory. Specifically, the vacuum condensate picture of quantum gravity suggests both the existence of an effective infrared cutoff (related to the observed cosmological constant) and a weak renormalization group (RG) running of G on very large scales. Here we present how these corrections will affect the power spectrum, and discuss how these effects can be used to distinguish the quantum gravity picture from the classical inflationary picture, with increasingly precise data becoming available in the near future. Section 7 summarizes the key points of our 2 Background -Correlation Functions and Power Spectrum In this section, and as a background, we briefly summarize the two main quantities that are used in cosmology to quantitatively parameterize matter fluctuations: (i) correlation functions in both real and momentum space, and (ii) their corresponding powers, or spectral indices. To begin, it is customary to describe the matter density fluctuations in terms of the matter density contrast δ(x, t) [26,28], which measures the fractional over-density of matter relative to the background mean densityρ(t), whereρ(t) evolves only with time, in a way that is governed by the background Friedman-Robertson-Walker (FRW) metric. The most common measure of such fluctuations is the real space correlation function G ρ (r), which averages the density contrast over all space 3 or its momentum space conjugate, commonly known as the power spectrum In the second expression the time dependence is factored out as 4 δ(k, t) ≡ ∆(k)F (t), so that the power spectrum explicitly compares fluctuations as they are measured today, at time t = t 0 . Furthermore, both the real and momentum space correlation functions represent statistical averages, which satisfy homogeneity and isotropy, and hence are expected to only depend on the magnitude r ≡ |x − y| and k.
The power spectrum is directly related to the real space correlation function via the Fourier In particular, most galactic observations are taken at relatively low redshifts, where t ≈ t ′ ≈ t 0 , in which case the pre-factor ratio of F (t)'s reduces to unity. One often starts with so-called scale-invariant models in cosmology, which assume a simple power law for the power spectrum, characterized by scaling indices. The real-space correlation function is then usually parameterized by a single index γ [26], defined via In momentum space, the spectral index s is then defined via With this parametrization, the Fourier transform can be evaluated through (4), or formally and exactly in d dimensions which then gives the relation These scale-invariant models have been compared to increasingly accurate astrophysical data measurements over the past few decades. Earlier data [26] Figure 2 in the next section. A plot of the observational measurements is represented by solid dots in Figure 2 in the next section. 5 An important challenge in cosmology is to provide a theoretical explanation for these particular values of s and γ. Here it should be noted that correlation functions for matter fields can be related to that of gravitational fields via the Einstein field equations, and the latter are in principle calculable from a quantum field theoretical treatment of gravity. Therefore quantum gravity provides a natural framework for calculating and explaining the scaling exponent γ as well as the spectral index s. This is the general picture we will discuss in detail in the following section.
The previous section outlined the customary ways to quantify matter fluctuations via the correlation functions G ρ (r) and P (k), and the corresponding spectral indices γ and s. Here we present a possible explanation for the origin for matter fluctuations as a direct result from gravitational fluctuations, and show how the corresponding correlation functions can be calculated from first principles from a quantum-field treatment of gravity. In the end we will show how the theoretically predicted results are in surprisingly good agreement with the most recent galaxy data from the Sloan Digital Sky Survey (SDSS) telescope collaboration [37].
General relativity provides an unambiguous relation between curvature and matter distributions through the Einstein field equations 6 and it is therefore natural that any fluctuations in matter are directly related to fluctuations in curvature. Now, in any quantum theory, quantities fluctuate, and a quantum theory of gravity produces metric and curvature fluctuations for which correlation functions are in principle calculable from first principles. It should then be possible to relate unambiguously such correlation functions to matter fluctuations, via the Einstein field equations. Here we will first outline briefly how invariant curvature correlation functions can be calculated from a quantum theory of gravity [38]. Secondly, we will show how matter correlation functions can be determined from curvature correlations, and at the end we will compare these predictions to observational data.
Quantum gravity, like QED and QCD, is in principle a unique theory. In the Feynman path integral approach, only two key ingredients are needed, namely, a gravitational action and a functional measure over metrics. As shown originally by Feynman [19,20], the Einstein-Hilbert action plus a cosmological constant term represents the unique action to describe a massless spin-2 field.
Additional higher derivative terms are in principle consistent with general covariance, but only affect the physics at very short distances, and therefore will not be considered here, since we are mainly interested in large distance cosmological physics. From this perspective, the formulation of a quantum theory of gravity in this approach has essentially no adjustable parameters, in a way that is quite similar to QCD and non-Abelian gauge theories in general.
A special characteristic of quantum gravity is its perturbative non-renormalizablility and its highly non-linear nature. As in the case of other non-linear theories such as QCD, Yang-Mills theories, and the O(N ) non-linear sigma model, one can nevertheless extract universal quantities, such as critical exponents (such as the exponents ν, γ, δ etc. in the nonlinear sigma model), and genuinely nonperturbative characteristic scales (such as correlation lengths or Λ M S for QCD).
Such predictions, though obtained from a nonpertubative treatment, are today amongst some of the best tested results of quantum field theory [18]. Analogously, utilizing a number of nonperturbative approaches (such as 2+ǫ expansion and numerical evaluations of the path integral), one can extract universal scaling dimensions such as the exponent ν and the nonperturbative, renormalization group invariant, correlation length scale ξ. The latter is in turn related to the vacuum expectation value of curvature, which is measured via a large-scale cosmological constant λ obs . Both the exponent ν and the scale ξ have significant effects on correlation functions, which deviate from any semiclassical or free field (Gaussian) predictions. More importantly, these are potentially verifiable via observations, as we will show below.
In previous work [39,40], it was shown that for large distances, the invariant scalar curvature correlation fluctuation at fixed geodesic distance behaves as with a power 2n = 2(d − 1 ν ), where d is the dimension of spacetime. This shows that the two-point function, for distances smaller than the aforementioned characteristic length scale ξ, follows a simple scaling law that is purely dependent on the universal exponent ν. The value for the critical exponent ν has been calculated and estimated through various means, including 2 + ǫ expansion [41,42,43], large d expansion [44], numerical lattice calculations [39,40], exact results in 2+1 dimensions [45] and truncated continuum renormalization group methods [46,47,48,49,51,52,54,55,56]. Many of these results have been summarized recently, for example, in [22], where it is argued that current evidence points to ν ≃ 1/3 in d = 4 (see Figure 1). In particular, in the large-d limit, one finds ν = 1/(d − 1) [44] in addition to the usual scaling results for the relation for the power in Eq. (10), namely 2n = 2(d − 1/ν). This in turn implies for the power 2n = 2 for d = 4 and above. In the following we will proceed on the assumption (supported by extensive numerical calculations on the lattice [39,40]) that ν = 1/3 exactly in d = 4, and that consequently the power in Eq. (10) is exactly equal, or at least very close to, two. Another key quantity here is the renormalization group invariant scale ξ appearing for example in Eq. (10), and related to the quantum gravitational vacuum condensate. It was argued in [21,57,22] that the renormalization group invariant ξ is most naturally identified with the scaled cosmological constant via ξ ≃ 3/λ. Modern observational Figure 1: Universal scaling exponent ν as a function of spacetime dimension d. Shown are the 2+ǫ expansion result to one and two loops [43], the value in 2 + 1 dimensions obtained from the exact solution of the Wheeler-DeWitt equation [45], the numerical result in four spacetime dimensions [40], and the large d result ν −1 ≃ d − 1 [44].
values of λ obs yield an estimate of ξ ∼ 5300 Mpc. Since galaxies and galaxy clusters involve distance scales of around r = 1 − 10 Mpc ≪ ξ, the scaling relation in Eq. (10) should then be applicable for such matter. Using these values, one finds in Eq. (10) 2n = 2(4 − 3) = 2, and therefore one expects the gravitational curvature fluctuations to scale as Next we proceed to relate the curvature fluctuations G R (r) to the matter fluctuations G ρ (r) defined in Section 2. As mentioned in the beginning of this section, this procedure is unambiguous because of the gravitational field equations. In the limit of matter taking the form of a perfect fluid with negligible pressure 7 , the Ricci scalar satisfies the Einstein trace equation Now using its variation and inserting this into Eq. (11) gives: where G ρ (x) is the matter correlation function defined in Section 2 above. So the scaling behavior of G R in Eq. (11) unambiguously relates to G ρ : Sinceρ is not expected to scale with distance, the large scale scaling behavior of G R ∼ r −2 in Eq. (11) directly gives the r −2 scaling behavior of G ρ in large scales. Or, in terms of the galaxy correlation index of Eq. (5), γ = 2. This is indeed in good agreement with the observed values γ ≃ 1.8 quoted earlier at the end of Sec. 2. The scaling for the corresponding power spectrum in momentum space can now be determined from the scaling of the real space matter correlation function predicted above. According to Eq. (8), the prediction of γ = 2 for large scales by quantum gravity implies that s = 1, or explicitly, for small k's. Note that for large k's (i.e. small distances) however, as the correlation function starts to probe distances far below the average separation of galaxy clusters, or even within one, the standard linear, isotropic assumptions are no longer valid. In these scales, complex processes and dynamics are expected to dominate over the large-scale gravitational correlation effects that we are considering here. Hence, for the correlation functions to exhibit a clear s = 1 scaling, one should consider scales larger than average galaxy cluster separations.
Typical galaxy clusters have diameters between 1 − 10 Mpc, and voids, the pockets of empty space between clusters, have typical diameters around 25 Mpc [58]. As a result, we would expect the effects of scaling to be most explicit for scales larger than around 50 Mpc or 100 Mpc, which correspond to k's smaller than ∼ < 0.15 h/Mpc. On the other hand, the spectum is expected to diverge from s = 1 for small distances (large k). In other words, Eq. (16) implies that for small k's (large scales), the power spectrum should approach a horizontal constant asymptotically on a kP (k) vs. k plot, with the constant representing the amplitude a 0 . In principle, a quantum theory of gravity should also lead to an estimate for the amplitude, as discussed and given in [22]).
However amplitudes are generally not universal, and depend on the choice of regularization scheme.
Therefore, in this paper, we prefer to focus on the prediction of the slope, with the amplitude taken as a free parameter that is fixed by fitting to the observational data.
We can now compare the predicted scaling relation of P (k) = a 0 /k with recent galaxy data from the SDSS collaboration measurements on a kP (k) vs. k plot. From Figure 2, it can be seen that, for sufficiently large distances (k ∼ < 0.15 h/Mpc), the general trend of the data indeed converges it can be seen that the general trend of the data fits quite well to the s = 1 slope predicted by quantum gravity, with an amplitude fitted to a 0 ≃ 689 ± 30(M pc/h) 2 . We note here that the amplitude a 0 is, in principle, calculable from the lattice treatment as discussed for example in [22], nevertheless it represents a non-universal quantity and thus clearly depends on the specific way the ultraviolet cutoff is implemented in the quantum theory of gravity, a fact that is already well known in lattice QCD [59]. Here we choose instead to take the non-universal amplitude a 0 as an adjustable parameter, to be fitted to the observational data.

Transfer Function and Normalization of the Power Spectrum
So far, our analysis asserted that the scaling of galaxy distributions, which is governed by that of matter fluctuation correlations G ρ (r), is directly calculable from curvature fluctuation correlation functions G R (r) in the quantum theory of gravity. In particular, the matter and curvature fluctuations are unambiguously connected via the Einstein field equations. Usually there is an implicit assumption that the T µν on the right hand side takes the perfect fluid form, which involve an equation of state of pressure-less matter. The latter assumption is certainly valid for the matterdominated universe 8 , an era in which we expect galaxy-formations to take place. As can be seen in Figure 2, the observational data is in rather good agreement with the s=1 prediction of quantum gravity for the scale of galaxy clusters (i.e. k 0.15 hMpc −1 ).
In principle, quantum gravity with its long-range correlations and n-point functions should also govern and make predictions beyond galactic observation scales. Below galactic scales (large k), complex and non-linear interactions between matter arises, and gravitational correlations are expected to become subdominant effects and not easily discernible. On the other hand, beyond the scale associated with galaxy clusters and superclusters (small k), gravity is again undoubtedly the dominant long-range force. There we expect fluctuations in curvature to again play a dominant role in governing long-range correlations. However, observationally, larger spatial scale corresponds to looking at increasingly earlier epochs of the universe, perhaps even before galaxies are formed.
Therefore data from galaxy surveys cannot serve as a test to the predictions in these regimes.
Fortunately, modern astronomy has provided very detailed measurements of the cosmic microwave background (CMB), which encodes fluctuations over a wide range of cosmological scales. This can then be used as an additional test for the quantum gravity predictions.
In this section, we will discuss how we can utilize the so-called transfer function to extrapolate the quantum gravity prediction to increasingly small k far beyond galaxy scales, and compare to how it fits observed CMB data. It should be noted that in the current CMB literature one often refers to the angular power spectrum C l as the primary quantities obtained. Nevertheless, P (k) is directly related to C l , as described in many standard literature works on the subject [28,60,61,26].
However here instead we focus on P (k), since, as mentioned above, it is this quantity whose scaling is most directly related to that of curvature fluctuations, the quantity that is predicted from quantum 8 see footnote 7 in Sec. 3 gravity. According to the Friedman equations, wavelengths that enter the horizon before matterradiation equality evolve differently than those that enter after. As a result, as we look at even larger scales (smaller k's), the s = 1 slope discussed previously is expected to change. Nevertheless, we can use a well-known function in cosmology, the so-called transfer function T (κ), to relate small wavelength behaviors to large wavelength behaviors. Details of which can be found in standard texts such as [28]; for explicit numerical results see also [62]. Here we will outline below how the use of the transfer function can be applied to our predictions.
We recall the definition of the spectral function Here F (t) can be calculated from the Friedman equations, giving where x = x(t) = (Ω Λ /Ω M ) a(t) 3 , and the correction factor C(x) is given by Detailed steps in the above derivations can be found in [28]. As for the function ∆(k) in (17), it can be calculated from the Boltzmann transport equations, again presented in detail in [28] and references therein, giving where the comoving wavenumber q is often used in place of the physical wavenumber k, the two being related by k ≡ q/a(t). Here δ(q, t) is the matter density contrast defined in Eq. (1) in momentum space, ψ(q, t) is the gravitational perturbation in Newtonian gauge, t L is the time of decoupling of radiation from matter, associated with the recombination of hydrogen, and H L and a L are the Hubble rate and scale factor correspondingly evaluated at t L . The last equality is obtained by conveniently parameterizing the combination of δ and ψ with R q and T (κ), known as the comoving curvature perturbation and transfer function, respectively [28]. The transfer function T (k) is usually expressed with an argument κ ≡ √ 2 (k/k EQ ) , where k EQ is the scale at matterradiation equality. Inserting these results into the definition for P (k) in Eq. (17) then gives of q where q * is some arbitrarily chosen reference scale, and n s is referred to as the spectral index. Then P (k) can be conveniently factorized into so that C 0 is a prefactor that encodes exclusively cosmological model parameters. It was Harrison and Zel'dovich who originally suggested in the seventies that n s should be close to one [63].
It is known that the transfer function T (κ) is entirely determined from classical cosmological evolution and the solutions of associated coupled Boltzmann transport equations for matter and radiation. The function in question can be accurately described by the semi-analytical formula given explicitly by D. Dicas in [28,62], which will be used here in the following discussion. For later reference, the shape of the function κ|T (κ)| 2 is shown below in Figure 3 ; particularly noteworthy here is the inverted-v shape reflecting the cosmological evolution transition from a radiation dominated to a matter dominated universe, again as discussed in detail in [28].
Now, both C 0 and T (κ) are fully determined, the former by cosmological measured parameters from, for example, the latest Plank satellite data [64], and the latter theoretically from the Boltzmann transport equations, numerically evaluated in [28,62] and further references cited therein.
Therefore the end result is that P (k) is essentially parameterized by two quantities, an overall amplitude A and the spectral index n s . The next step is to use quantum gravity to theoretically constrain the value of n s , and analyze how well a power spectrum with this predicted value of n s fits the current observational data. As discussed previously, quantum gravity predicts P (k) ∼ a 0 /k s with an exponent s = 1 in the galaxy regime (i.e. k ∼ 0.01 to 0.3 hMpc −1 ). In the previous section, we noted that although a 0 is in principle calculable, additional subtleties arise. Hence, for the current purpose, we simply use the value that fits well the galaxy data at the largest scales, [ / ] Figure 3: Shape of the transfer function κ|T (κ)| 2 using the interpolating formula of Eq. (25). What is of interest here is the significant turnover happening at wavenumbers k ∼ 0.02 which is known to correspond to a cosmological time scale associated with matter-radiation equilibrium. This marked turnover here is the primary reason for the peculiar inverted-v shape of the power spectrum in the following plots.
namely a 0 ≃ 689. As a result, within the galaxy regime, we should have a matching P (k) ≡ C 0 A k ns T (κ) ≈ 689 k ( for k ∼ 0.01 to 0.3 hMpc −1 ) .
As mentioned earlier under Eq. (25), both C 0 and T (κ) are fully determined from the FRW evolution equations, leaving this essentially an equality with two unknowns, namely A and n s , that holds in the galaxy regime. Therefore, by appropriately selecting two points from the quantum gravity prediction within this regime, we can derive the values for A and n s . Since the left hand side is supposedly valid for all scales, with the two values A and n s determined, these will then give us a power spectrum that allows us to extrapolate the quantum gravity prediction within the galaxy regime to much larger scales.
Within the galaxy regime, there is a certain flexibility in which two points one may choose.
For the purpose of a first estimate, we select two points relatively apart, but not too close to the margins of the galaxy data. Although the data set ranges from 0.025 These correspond to real space scales of Using these test points in Eq. (26) gives n s = 1.108 , A = 5.97 × 10 11 .
The best fit from the Plank 2015 data gives a fit of n s = 0.9667 ± 0.0040 [64]. This suggests our preliminary analysis, using the two above selected points, yields a value of n s that is around ∼ 15% of the best fit value of the data. this normalization choice seem to yield a value for n s closer to that of Plank, the two points do not seem sufficiently apart to normalize the spectrum properly in our quantum gravity picture.
Nevertheless, this methodology also provides a rough estimate for the overall uncertainty in n s .

Comparison with Inflationary Models
A popular attempt to explain the matter power spectrum is through inflation models, which proposes that the structures visible in the Universe today are formed through fluctuations in a hypothetical primordial scalar inflaton field. Although the precise behavior and dynamics of the inflaton field are still largely controversial [9,10,11,12,13,14], there are some general features common among these models that allow one to derive predictions about the shape of the power spectrum.
Here we very briefly outline the inflation perspective and how it radically differs from the quantum gravity perspective we have presented earlier. For more details and reviews on the inflation mechanism, we refer to recent literatures such as the books [28], [60] and [61].
Cosmic inflation was first proposed in [6,7,8] as a possible solution to the horizon and flatness problem in standard big bang cosmology, by proposing a period of exponential expansion in the early phases of the universe. This exponential expansion is expected to be driven by a hypothetical"inflaton field", usually scalar in nature, which dominates the energy density in the early universe, causing a de Sitter type accelerated expansion. In these models, the gravitational perturbations that ultimately cause gravitational collapse and structure formations are also due to fluctuations in this inflaton field. Scalar field fluctuations are assumed to be Gaussian for large enough distances (small k), which then naturally results in a value for the spectral index n s = 1. More realistic models with non-zero tensor to scalar ratio predicts n s between 0.92 and 0.98 [65,66,11,67], without assuming excessive fine-tuning of parameters [11] 9 The above value of n s = 1 is then used in Eqs. (23) and (25), leaving only one free parameter, the amplitude A in Eq. (23), which can then be determined by normalizing to the CMB data at very large scales (i.e. small k).
Notice that the gravitational perspective we suggest in this paper differs fundamentally from that of inflation, both in procedure and, most importantly, in origin. Firstly, inflation models suggest that gravitational perturbations are due to fluctuations in an inflaton field, whereas here perturbations are intrinsically gravitational and quantum mechanical in origin. Secondly, inflation models use the corresponding inflaton correlation function, G φ (r), or power spectrum, P φ (k), which are mostly scalar/Gaussian in nature, to derive those of matter G ρ (r), P (k), whereas in our gravitational p icture, the correlation function that constrains G ρ (r) and P (k) is the curvature G R (r), which is highly non-Gaussian. Finally, while both perspectives provide a prediction for n s which governs the shape of P (k), both leave the overall normalization A in Eq. (23) uncertain.
Modern renormalization group theory would imply that the critical exponent ν, and the scaling dimensions n that follow from it [see Eq. (10)], are expected to be universal, and as such only dependent on the spacetime dimension, the overall symmetry group, and the spin of the particle.
The same cannot be said of the critical amplitudes, such as the amplitude associated with the curvature two-point function of Eq. (10). In principle it is possible to estimate from fist principles the amplitude A and thus the normalization constant N is Eq. (22), since these can be regarded as wave function normalization constants in the underlying lattice theory. A rough estimate was given in [22], nevertheless a number of reasonable assumptions had to be made there in order to relate correlations of small gravitational Wilson loops to very large (macroscopic) ones. As a result, the estimate for the amplitude quoted there is still two orders of magnitude smaller than the observed one. Nevertheless it is understood that amplitudes are expected to be regularization and scheme dependent, and will differ by some amount depending on the specific form of the ultraviolet cutoff scheme chosen, whether it is a lattice one, or a continuum inspired momentum cutoff, or dimensional regularization. This phenomenon is well known in scalar field theory as well as lattice gauge theories, and allows the relative correction factors to be computed to leading order in perturbation theory [59], thus reducing the regularization scheme dependent uncertainties.
In the following we choose to leave, for now, the overall amplitude in Eqs. (10) or (16) as a free parameter, and constrain it directly instead by the use of observational data.
Inflation models on the other hand often normalize the curve at large scales with CMB data or at the turnover point P max [28], but our picture normalizes it with the slope within the galaxy domain, where we argue that the curvature-matter relation is most direct, as discussed in Sec.
3. We contend that there are a number of reasons that the gravitation-induced picture is more natural for matter distributions than the inflaton-induced one. Firstly, the gravitational field is a well-established interaction and is supposed to have long-range influences with cosmological consequences, whereas the scalar inflaton field remains an observationally unconfirmed quantum field whose precise properties, such as its potential and interaction with standard model fields, remain unknown. Secondly, the Feynman path integral treatment is an unambiguous procedure to quantize a theory covariantly, producing in principle a unique theory with essentially no free parameters (as in the case of QED and QCD). On the other hand, inflation models currently are far from unique, and many are plagued with significant fine-tuning problems [9]. At this stage, despite the differences in perspective, there is yet no clear-cut preference over each other except for the naturalness arguments given above. However, a number of second order effects due to quantum gravitational fluctuations could yield diverging and distinct predictions, which should be testable in the near future with increasingly accurate measurements. These new effects are the subjects of discussion in the next section.
In a quantum field treatment of gravity, a nonperturbative scale ξ is dynamically generated which regulates the otherwise serious infrared divergences associated with a zero mass graviton. In addition, Newton's constant G is seen to "run" with scale as a result of quantum corrections, in analogy to what happens in QED and QCD. In momentum space, the formula for the running of G is given in [23,22] by Here G 0 is the classical (i.e. k ≫ m) "laboratory" value of the Newton's constant, or simply G henceforth in the paper, m is the renormalization scale in momentum space related to ξ by m ≡ 1/ξ, and c 0 the coefficient for the overall amplitude of the quantum correction, which is generally expected to be of order one. As for the nonperturbative scale ξ, it is most naturally related to the cosmological constant λ 0 [21,57] by Notice that these two quantum effects are not unique to gravity, but common in quantum field theories. The latter is most representative in the well-studied theories of QED and QCD, where quantum corrections result in the running of coupling constants e and α s respectively. It is well known that a dynamically generated infrared cutoff arises in other nonperturbative theories such as QCD and non-Abelian gauge theories. The scale ξ therefore plays a role analogous to the scaling violation parameter Λ M S in QCD. ξ also serves directly as a characteristic scale in the theory that, much like the Λ M S scale in QCD, distinguishes the small distance from the large distance domain.
Both effects are discussed in detail in [21,22].
The formal implementation of the dynamically generated scale ξ is to be inserted as a lower infrared cutoff in any momentum integrals. This would mean amending However, a more straightforward implementation is to simply make replacements k 2 → k 2 +m 2 with m = 1/ξ, which phenomenologically works well in other nonperturbative theories (such as QCD) and partially includes the effects of infrared renormalons. As a result, for the power spectrum given in Eq. (23) becomes in the extreme long distance (small k) domain where k ∼ 1/ξ.
Finally, another important expected quantum effect is the renormalization group (RG) running of Newton's G, which was already discussed in some detail in [23,21,57] and implemented either through a scale dependence in momentum space G(k) or by the use of a set of covariant nonlocal effective field equations containing a G( ), where ≡ g µν ∇ µ ∇ ν is the covariant d'Alembertian.
Here, to implement the running of Newton's G, one needs to identify where G would appear in the power spectrum. Recall the matter correlations are related to curvature correlations via the Einstein field equations or in momentum space, via a Fourier transform where, as usual, P (k) is the spectrum for matter fluctuations, and P R (k) ≡ < |δR(k)| 2 >∼ 1/k s with index s = 1 is the spectrum for curvature fluctuations. Observing that P (k) ∝ 1/G 2 dimensionally, promoting G → G(k) can be achieved by the replacement with G(k) given in Eq. (30). In addition, since the running of G is again an effect that is most significant on small k's, the strong infrared divergence near k ≃ 0 should similarly be regulated as discussed above correspondingly by the replacement k 2 → k 2 + m 2 in G(k). As a consequence, the fully regulated expression with the running factor is This expression is plotted as the green curve in Figure 5. In summary, from Figure 5 we can see that the effects of the IR regulator alone serves to level off the curve at low k around k ∼ m = 1/ξ (orange curve), whereas the full modification (green curve), which includes both the IR regulator and the RG running of G, bends the curve downwards below the unmodified one (blue curve). Once again, it should be stressed that these two additional, genuinely quantum, effects are fairly concrete predictions associated with quantum modifications to classical gravity, and are quite distinct from any inflationary models, both in procedure and origin. With the advance of increasingly accurate astrophysical and cosmological satellite observations, it is hoped that these new predictions could be verified, or falsified, in the near future. Figure 5: Effects of the IR regulator and RG running of G. The blue curve represents the full power spectral function P (k) as derived and shown previously in Figure 4. The orange curve P (k) reg includes the effect of an IR regulator. The green curve P (k) run includes the effect of an IR regulator together with the RG running of Newton's G.

Conclusion
In this paper we have provided an explanation for the galaxy and cosmological matter power spectrum that is purely gravitational in origin, replacing the more popular explanation using inflation.
We first showed how gravitational fluctuations would unambiguously govern the slope of the spectrum in the galaxy survey domain. As seen in Figure 2, the prediction is in very good agreement with observations. Next, using the transfer function we extrapolated the prediction to all scales, which is also in general agreement with observations ( Figure 4). Later, we investigated additional quantum effects that are expected to be important at large scales (small k), and show how these predictions diverge from conventional inflation-inspired predictions ( Figure 5).
The primary benefit of this explanation over inflation is that there is no need for additional hypothetical ingredients of physics, other than Einstein's gravity and accepted quantum field theory methods. The correlation functions of gravitational fluctuations are fully governed by correlators in a quantum field theory of gravity, as for fluctuations of any other interaction. It is widely known that gravity is nonrenormalizable in a perturbative treatment; hence standard nonperturbative techniques based on the Wilson-Parisi renormalization group analysis [15,16,17] need to be used, and then the respective critical exponents are fully calculable. In this paper we show how these results translate to predictions for relative spectral indices for matter fluctuations unambiguously via the field equations. It is argued that the further extension of this result to all scales only requires the transfer function, a well-known ingredient in cosmology that involves mainly collision analysis via the Boltzman equations. In this picture, no new physical ingredients are postulated, nor necessary. On the other hand for inflation, a new, minimum of one, inflaton field, usually scalar in nature, must be involved. We therefore argue the gravitational picture provides a more concrete and natural explanation to the origin and distribution of cosmological matter fluctuations. Finally, the gravitational fluctuation picture also provides a clear prediction that diverges from scalar field induced predictions on large scales. As advanced satellite experiments are continuously being conducted, and increasingly accurate measurements are becoming available, the predictions as a result of quantum gravity made in the last section could be verified or disproved in the near future.
At first it would seem that the assumptions of Gaussian correlations, and therefore of a Gaussian power spectrum, for the scalar field would be rather restrictive. There are in principle many possible scalar field self-interaction terms that can be written down in four dimensions, starting with monomials of the fields up to possibly even non-local interactions, all leading potentially to wildly different power spectra. Nevertheless Wilson's argument [15] (and subsequent related extensive rigorous work) in support of triviality of lambda phi four (λφ 4 ) theory in and above four spacetime dimensions, based in turn on by now well established field theoretic renormalization group arguments, suggests that, for a wide class of local scalar self-interacting theories, the Gaussian (free field) fixed point in the most relevant quartic coupling constant will act as an infrared attractor. Here the implication of this deep field-theoretic result is that, for a wide class of scalar self-interacting local theories, the assumption of Gaussianity at large distances is reasonably well justified in terms of rather general renormalization group arguments. On the basis of these arguments one would then expect that all scalar n-point functions should be reasonably well described by their free field expressions, as obtained from the application of Wick's theorem. At the same time, it should be clear that in the current gravity-based model the long range correlation functions are most certainly not Gaussian, due to the presence of non-trivial anomalous dimensions.
Specifically, the values for the non-trivial scaling dimensions listed above imply that the behavior of gravity and matter n-point functions (and their related power spectra) is far from their Gaussian, or free field, counterpart. As an example, note that the two-point function result of Eq. (10), and related to it the universal gravitational scaling dimension n = 1, also determines the form of the reduced three-point curvature correlation function [22] < √ g R(x 1 ) √ g R( with C 123 a constant, and relative geodesic distances d ij = |x i − x j |, etc. As before, and again by virtue of the quantum equations of motion, these identities can then be related to three-point functions involving the local matter density, or more generally to expressions involving the trace of the energy momentum tensor. Proceeding along the same line of arguments, analogous expression can be given for the four-point functions, which again will exhibit a distance dependence unambiguously fixed largely by the scaling dimension of the scalar curvature operator, with non-universal amplitudes. We note here that the relevance and measurements of nontrivial three-and four-point matter density correlation functions in observational cosmology was already discussed in great detail some time ago in [26]. Then the results presented here imply that such higher order n-point function could, in the not too distant future, provide additional stringent observational tests on the vacuum condensate picture for quantum gravity, and on the non-trivial gravitational scaling dimensions scenario described here.
The ability to reproduce the cosmological matter power spectrum has long been considered one the "major success" for inflation-inspired models. Although within our preliminary study, further limited by the accuracy of present observational data, it is not yet possible to clearly prove or disprove either idea, the possibility of an alternative explanation without invoking the machinery of inflation suggests that the power spectrum may not be a direct consequence or a solid confirmation of inflation, as some literature may suggest. By exploring in more detail the relationship between gravity and cosmological matter and radiation, together with the influx of new and increasing quality observational data, one can hope that this hypothesis can be subjected to further stringent physical tests in the near future.