Vacuum Condensate Picture of Quantum Gravity

In quantum gravity perturbation theory in Newton's constant G is known to be badly divergent, and as a result not very useful. Nevertheless some of the most interesting phenomena in physics are often associated with non-analytic behavior in the coupling constant and the existence of nontrivial quantum condensates. It is therefore possible that pathologies encountered in the case of gravity are more likely the result of inadequate analytical treatment, and not necessarily a reflection of some intrinsic insurmountable problem. The nonperturbative treatment of quantum gravity via the Regge-Wheeler lattice path integral formulation reveals the existence of a new phase involving a nontrivial gravitational vacuum condensate, and a new set of scaling exponents characterizing both the running of G and the long-distance behavior of invariant correlation functions. The appearance of such a gravitational condensate is viewed as analogous to the (equally nonperturbative) gluon and chiral condensates known to describe the physical vacuum of QCD. The resulting quantum theory of gravity is highly constrained, and its physical predictions are found to depend only on one adjustable parameter, a genuinely nonperturbative scale xi in many ways analogous to the scaling violation parameter Lambda MSbar of QCD. Recent results point to significant deviations from classical gravity on distance scales approaching the effective infrared cutoff set by the observed cosmological constant. Such subtle quantum effects are expected to be initially small on current cosmological scales, but could become detectable in future high precision satellite experiments.


Introduction
Like QED and QCD quantum gravity is, in principle, a unique theory. In the Feynman path integral approach only two key ingredients are needed to formulate the quantum theory, the gravitational action and the functional measure over metrics. For gravity the action is given by the Einstein-Hilbert term augmented by a cosmological constant. Additional higher derivative terms are consistent with general covariance, but nevertheless only affect the physics at very short distances, and will not be considered further here. 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 given by the DeWitt form. As in most other cases where the Feynman path integral can be written down (including non-relativistic quantum mechanics), the proper 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, at least in principle, the resulting quantum theory of gravity does not seem to require any additional extraneous ingredients, besides the ones mentioned above. Indeed some time ago Feynman was able to show that Einstein's theory is unique, arising from the consistent quantization of a massless spin two particle.
At the same time, gravity has been known to present some rather difficult inherent problems.
The first one is that the theory is intrinsically nonlinear, since gravity gravitates. In addition, perturbation theory in Newton's constant G is useless, since the resulting series is badly divergent (much more so than in QED and QCD) which makes the theory not perturbatively renormalizable.
It is also a known fact that the gravitational action is affected by a conformal instability, which makes at least the Euclidean path integral potentially divergent. Last but not least additional, genuinely gravitational, technical complications arise due to the fact that physical distances between spacetime points depend on the metric, which is a fluctuating dynamical quantum entity.
Serious divergences that appear in perturbation theory originate from the fact that the gravitational action leads to vertices which are proportional to a momentum squared. When these vertices are inserted into diagrams, they give rise to ultraviolet divergences which get increasingly worse as the order of perturbation theory is increased. The lack of perturbative renormalizability leads therefore to two main, alternative and clearly mutually exclusive, conclusions. One states that the quantum theory of gravity does not exist due to these cascading perturbative divergences, and consequently an enlarged, improved theory should be investigated instead. Enlarged theories that attempt to make quantum gravity perturbatively renormalizable include N = 8 supergravity, and supersymmetric strings in ten spacetime dimensions. The other alternative path, followed here, is that the usual diagrammatic methods of QED and QCD fail for gravity because perturbation theory is incomplete or invalid, presumably due to a more complex analytic structure in the coupling constant, which then leads to gravity not being renormalizable in the usual perturbative sense. The possibility exists therefore (and is further supported by several well known examples in physics) that perturbation theory in G fails because physically relevant quantities (n-point functions, quantum averages, functions describing the running of G with scale etc.) are non-analytic at G equal zero. 2 Indeed there are many physically very interesting and deep phenomena which cannot be explained, or even studied, using perturbation theory alone. One example is QCD, were gluons and quarks are confined with a chromoelectric string tension known to be non-analytic (in the form of an essential singularity) in the gauge coupling. Again, in a superconductor the correct ground state is described by Cooper pairs bound together by a weak electron-phonon interaction. The latter leads to a gap in the energy spectrum close to the Fermi surface, which is known to be non-analytic in the fundamental electron-phonon coupling constant. In a superfluid the quantum condensate density is non-analytic in the coupling as well, and so is the screeing mechanism in a degenerate Coulomb gas, where the Thomas-Fermi screening length is known to be non-analytic in the charge. In this last model the correct charge screening mechanism is not reproduced to any finite order in perturbation theory. Nevertheless, it is easily obtained by resumming infinitely many so-called ring diagrams. Additional physically relevant examples include homogeneous turbulence (which is described by nontrivial Kolmogoroff scaling exponents) and order-disorder transitions in ferromagnets and related systems. The latter exhibit spontaneous symmetry breaking, dimensional transmutation, nontrivial scaling dimensions, and the appearance of a non-vanishing field condensate in the ordered phase. Related to the last example is the case of a self-interacting scalar field above four dimensions, where, on the one hand, the theory is known to be perturbatively non-renormalizable with the kind of escalating ultraviolet divergencences described earlier. Yet one can prove rigorously, by using the lattice path integral formulation, that the model reduces to a non-interacting (Gaussian) theory at large distances, with low energy scattering amplitudes vanishing as an inverse power of the ultraviolet cutoff. 2 The validity of the perturbative approach to gravity is sometimes supported by the fallacious argument that in some sense "gravity is weak". That is certainly true when gravity is compared to the other fundamental forces on laboratory scales. Nevertheless, unlike QED and QCD, the gravitational coupling is dimensionful which makes such weak coupling arguments invalid, or at least naive, when referred to gravity as its own self-sufficient theory. Ultimately the real question is whether large quantum gravitational field fluctuations, which cannot be excluded a priori from the path integral, are physically important or not.
In many cases the common thread among these widely different theories and physical phenomena is the existence of some sort of vacuum condensate, which generally turns out to be a non-analytic function of the relevant fundamental coupling constant. The origin of these non-analiticities can often be traced back to the fact that the physical ground state is, in the end, fundamentally different from the original unperturbed or free field ground state. So the true ground state is qualitatively different from the unperturbed ground state which initially forms the starting point for perturbation theory. Physically, a significant rearrangement of the vacuum will often not just involve small perturbations, and generally cannot be obtained by perturbative methods, which implicitly assume smooth changes and thus the existence of a Taylor series in the relevant coupling. In this framework, the failure of perturbation theory is seen more as a reflection on the fundamental inadequacy of the mathematical methods used, and not necessarily as a shortcoming of the underlying fundamental theory per se.
If gravity is not perturbatively renormalizable, then what are the alternatives? In fact perturbatively non-normalizable theories have been theoretically rather well understood since the early seventies, when the modern renormalization group approach (based on momentum slicing, scaling dimensions and multidimensional coupling constant flow) was invented to account for more subtle and complex behavior in quantum field theory [1,2,3,4,5,6]. Moreover, a number of significant examples exist of theories which are not perturbatively renormalizable and nevertheless give rise to physically acceptable and interesting theories, and for which very detailed and accurate physical predictions can be produced. Most often these involve models formulated in less than four dimensions, which are thus, generally, more relevant to statistical field theory than to particle physics or gravitation. Indeed in a statistical field theory context it is often possible to bypass the limitations of perturbation theory, by resorting to additional, but complementary, approximation and expansion methods. These methods include Wilson's 2 + ǫ and 4 − ǫ expansions, the large N expansion, weak and strong coupling expansions (sometimes referred to as the low and high temperature expansion), partial resummation methods, and finally a combination of all of the above methods paired with high accuracy direct numerical evaluations of the original path integral or partition function (for an overview, see for example [7,8,9,10,11]).
One of the reason why these more powerful methods are eventually capable of providing useful (and ultimately correct) physical information about the systems studied lies in the fact that they are able to access new nontrivial strong coupling fixed points of the renormalization group, which are often not at all visible nor accessible in weak coupling perturbation theory. In other words, the common thread among many of the models that are perturbatively non-renormalizable -but which in the end turn out to be physically acceptable and relevant -is the existence of a nontrivial strong coupling renormalization group fixed point. Furthermore, in support of the legitimacy of such a more sophisticated approach one should mention, as an example, the fact that exquisitely detailed predictions for a class of perturbatively non-renormalizable theories, namely the O(N ) non-linear σ-model [12,13] in three space dimensions, now provide the second most accurate test of quantum field theory [14,15], after the QED predictions for the g − 2 anomalous magnetic moment of the electron.
It seems therefore reasonable to apply the very same (and by now well established and very successful) methods to one more perturbatively non-renormalizable theory, namely the quantum theory of gravity in four dimensions [16,17,18]. One first notes that a controlled non-perturbative approach clearly requires a useful and explicit ultraviolet regulator, and the only known reliable way to evaluate non-perturbatively the Feynman path integral in four dimensional quantum field theories is via the lattice formulation. Indeed, as shown in detail already by Feynman for nonrelativistic quantum mechanics, the very definition of the path integral (which sums over all paths, known to be generally nowhere differentiable) requires the introduction of a lattice discretization, due to the Wiener path nature of quantum trajectories [20]. One shining example of the success and reliability of the lattice approach is the elucidation of the subtle mechanism of confinement and chiral symmetry breaking in QCD.
The explicit introduction of a lattice achieves two purposes, one is to provide an explicit discretization (which is required in order to define in an explicit, as opposed to formal, way what is meant by the sum over all paths), and second to give a necessary regularization (in the sense of taming the ubiquitous field theoretic short distance divergences) of the quantum path integral.
Additional advantages of the path integral formulation, present both in the case of gauge theories and gravity, are the existence of a manifestly covariant formulation, and the known fact that no gauge fixing is in principle required (as first shown by Wilson in the gauge theory case [21]) outside the traditional framework of perturbation theory. Sometimes it is possible to rely on some sort of saddle point expansion around a smooth solution to the classical field equations, nevertheless it is generally recognized that dominant paths that contribute to the path integral are nowhere differentiable, and ultimately can only be accounted for properly in a controlled discretized formulation.
While it is certainly possible to evaluate the gravitational path integral using perturbation theory, the latter is nevertheless not the only avenue open, and is seen in fact as rather restrictive for the reasons outlined before. The case of the non-linear sigma model shows rather clearly that results derived using perturbation theory alone can be entirely misleading, and do not capture correctly the underlying physics of the ground state and key aspects related to the existence of an order-disorder transition. Also, the lattice formulation for gauge theories and gravity presented somewhat of a novelty thirty years ago, but is now extensively tested for QCD, scalar field theories and a variety of spin systems. In addition, today one can rely on over thirty years experience in an array of both analytical and numerical calculations, and in their fruitful mutual interplay.
Unless one desires to reinvent entirely new and ad hoc methods, the natural prototype for dealing with genuine non-perturbative aspects of gravity is Wilson's lattice formulation of QCD.
Indeed, while QCD is perturbatively renormalizable, it is well known that in this case perturbation theory is largely useless at low energies, where confinement effects take over and fundamentally modify the physical picture of the vacuum state. One key aspect of the lattice gauge theory is that, in order to preserve a form of exact local invariance (and related quantum Ward identities) the formulation requires an integration over gauge fields with a nontrivial (but uniquely determined by local gauge invariance) Haar measure. Then, in the lattice framework, confinement is an almost immediate and easily visualized consequence of large field fluctuations at strong coupling.
QCD is a hard theory to solve, and many deep insights have come from the lattice formulation. It cannot be stressed enough that one important outcome of the lattice calculations is that the physical vacuum bears little resemblance to the perturbative vacuum, due to significant nonlinearities and nontrivial field condensation effects. The former exhibits a nontrivial spectrum of hardons and glueballs, chromoelectric and quark field vacuum condensates, all of which are ultimately nonanalytic in the gauge coupling g, and cannot be reproduced by perturbative methods. Indeed to this day Wilson's lattice theory provides the only convincing evidence for confinement and chiral symmetry breaking in QCD and, more generally, in non-Abelian gauge theories. In addition, the lattice theory allows credible calculations of the running of alpha strong versus energy, which compare rather well with current experimental data.
For a quantum theory of gravity, the Feynman path integral again represents a natural starting point [22,23,24,25]. It is therefore rather fortunate that an elegant lattice formulation for gravity was written down by Regge and Wheeler in the early sixties, and is, not unexpectedly, based on the key concept of a dynamical lattice [26,27] (for a recent overview see [28] and references therein).
The main features of this theory can be summarized as follows. It incorporates a continuous local invariance, completely analogous to the diffeomorphism invariance of the continuum theory. As already pointed out originally by Regge, the local invariance of the lattice theory then leads to a lattice analog of the Bianchi identities, and thus to corresponding Ward identities in the quantum version. It also puts within reach of computation problems in classical general relativity which are in practical terms beyond the power of analytical methods; this last aspect was perhaps one of the main motivations initially (in the early sixties) for a discrete formulation of General Relativity.
Furthermore, like most lattice field theories, it affords in principle any desired level of accuracy by a sufficiently fine subdivision of spacetime, allowing eventually a reconstruction of the original continuum theory. The resulting Regge-Wheeler lattice theory of gravity is generally known as simpicial quantum gravity, for the simple reason that it is based on a construction of space-time out of geometric simplices, four dimensional analogues of triangles and tetrahedra. In this formulation curvature is described by angles, metric components are replaced by edge lengths, and the relevant geometric quantities can be calculated from the values of the edge lengths to give local lattice volumes, angles and local curvatures. In other words, local curvature is completely determined by an assignment of edge lengths and by how each edge is locally connected to neighboring edges (the incidence matrix). It is then possible to write down the lattice analog of the local volume element, of the local Riemann tensor, of the scalar curvature, and therefore ultimately of invariant terms such as the Einstein-Hilbert action.
As a consequence, the first key ingredient of a discretized form for the Feynman path integral for gravity, namely the action, is provided by the Regge-Wheeler theory. Furthermore, since the path integral involves an integration over all four metrics, and since the metric is locally related to the lattice edge lengths squared, the implication is that the analog of the DeWitt functional measure over continuum metrics turns into an integration over all lattice edge lengths squared (with some suitable volume inequality constraints, so as to guarantee a sensible geometric interpretation). For ordinary field theories the rigorous construction of the Feynman path integral often involves a Wick rotation to complex spacetime, and the same procedure can be achieved in the context of gravity as well, both in the continuum and on the lattice.
While it is possible in some cases to proceed with a Lorentzian signature [22,23,16,17,19] (in the continuum, and on the lattice for example by the use of a discretized Wheeler-DeWitt equation), it is generally accepted that the Euclidean formulation provides a mathematically more sound description of the Feynman path integral [24,25]. In addition, such a formulation generally relies on weights involving positive real probabilities, which then allows the use of established numerical probabilistic methods. It is certainly possible that, in the context of gravity, the Lorentzian and Euclidean theories belong to two different universality classes, and give rise to two entirely different sets of renormalization group beta functions and scaling exponents. This would be rather unique, since no other instance of such an occurrence is known. Nevertheless the evidence so far suggests that basic results in the Lorentzian and Euclidean lattice theories agree quite well. An explicit test of this statement lies in the ongoing comparison of results for universal scaling dimensions obtained in the two formulations. A recent example of a Lorentzian formulation of lattice quantum gravity, also based on the Regge-Wheeler discretization, is an exact solution of the lattice WheelerDeWitt equations in 2 + 1 dimensions discussed in [29,30,31]. In addition, it is possible to force quantum gravity to become perturbatively renormalizable by formally expanding about two dimensions (essentially Wilson's 2 + ǫ expansion applied to the case of gravity). In this approach it seems quite clear that universal results such as scaling dimensions are expected to be identical for the Lorentzian and Euclidean signatures to all orders in perturbation theory [32,33,34,35].

Regularized Path Integral for Quantum Gravity
One usually considers as the starting point for a nonperturbative formulation of quantum gravity a suitably discretized form of the Feynman path integral, initially for pure gravity without matter fields, which can then be added at a later stage. In the continuum the path integral is given formally by [22,24] with the Einstein-Hilbert gravitational action Here the dots here indicate possible matter and higher derivative terms, with the latter getting generated, for example, by radiative corrections as they arise already in the framework of perturbation theory. The functional integration over metrics is done using the DeWitt diffeomorphism invariant measure [23] dµ[g µν ] = x µ≥ν In the above expression k −1 ≡ 8πG, with G the bare Newton's constant and λ 0 a bare cosmological constant. In the following we will consider almost exclusively the case of no higher derivative R 2 -type terms, and no dynamical matter (quenched approximation).
The continuum Feynman path integral given above is generally ill-defined (the integration is dominated by non-differentiable Wiener paths), and so it has to be formulated more precisely by introducing a suitable discretization, as is done in both non-relativistic quantum mechanics and quantum field theory [20]. This last step is particularly crucial for nonperturbative gravity calculations, where the nontrivial invariant measure over the g µν 's has been shown to play an important role. In the 60's Regge and Wheeler proposed an elegant discretization of the classical gravitational action [26,27], which forms the basis for the lattice formulation of quantum gravity used here; early references include [36,37,38,39,40]. Once the measure and the path integral have been transcribed on the lattice, the ultimate goal then becomes to recover the original continuum theory of Eq. (1) in the limit of a suitably small lattice spacing. It is known that taking this limit is a rather subtle affair, and, in order for it to be taken correctly, it will require the full machinery of the modern (Wilson) renormalization group.
A suitable starting point is therefore the following discrete form for the Euclidean Feynman path integral for pure gravity with a compactly written lattice gravitational action and lattice integration measure In these last expressions the sum over hinges h in four dimensions corresponds to a sum over all lattice triangles with area A h , with the deficit angle δ h describing the curvature around them [26,27].
The Θ-function constraint appearing in the discrete measure ensures that the triangle inequalities and their higher dimensional analogs are satisfied by all simplices. The discrete gravitational measure in Z L of Eq. (4) can then be viewed as a regularized version of the original DeWitt continuum functional measure of Eq. (3). A bare cosmological constant term with λ 0 > 0 is essential for the convergence of the path integral since for bare λ 0 ≤ 0 the Euclidean path integral is clearly divergent [38,40].
It is a rather useful fact that the lattice edge lengths are locally related in a simple way to the continuum metric. In terms of the edge lengths l ij attached to a four-dimensional simplex s one has for the induced metric within that simplex where the four-simplex here is based at the point 0. This last result then provides the needed connection between the continuum metric g µν (x) and the lattice squared edge lengths degrees of freedom l 2 i ; the latter is essential in establishing a clear and unambiguous relationship between lattice and continuum operators, just as in the case of Yang-Mills theories on the lattice. Appropriate lattice analogs of various curvature invariants can then be written down, making use of the well-understood correspondences In the above expressions the hinges h correspond to triangles in four dimensions. A detailed discussion of such operators, as well as additional four-dimensional curvature invariants, can be found for ex. in [28] and references therein. In evaluating the lattice path integral, by whatever means, the continuum functional integration over metric is thus replaced by a finite-dimensional integration over squared edge lengths, which become the fundamental variables in the discrete theory. The general aim of the calculation is then to evaluate the lattice path integral either approximately or exactly by numerical means, by performing a properly weighted sum over all lattice field configurations.
In lattice field theories it is customary to deal with dimensionless quantities [7,8,9,10,11], and here this rather well-established procedure is followed again, for obvious reasons. The bare coupling constants λ 0 and G appearing in the continuum theory are expressed from the start in units of a fundamental lattice cutoff Λ = 1/a; without such a cutoff the continuum theory is generally ill defined [20]. The latter is then set equal to one, so that all observable quantities, correlators and couplings are expressed in units of this fundamental cutoff. In the end the actual value for the cutoff (say in cm −1 ) is determined by comparing suitable physical quantities. Furthermore, the functional integral depends on several bare coupling constants, but it is important to note that in the absence of matter the theory only depends on one bare parameter, the dimensionless coupling k/ √ λ 0 . This is easily seen, for example, from the fact that in d dimensions a constant rescaling of the metric turns the cosmological constant term λ 0 √ g into λ 0 ω d/2 √ g ′ , so that a subsequent rescaling of the bare coupling constants leaves the dimensionless combination G d λ d−2 0 unchanged. One concludes that only the latter combination has a physical meaning in pure gravity; in particular one can always suitably chose the scale ω = λ −2/d 0 so as to adjust the volume term to acquire a unit coefficient. This ability to rescale the field variables (the metric) so as to reabsorb certain renormalizations of the couplings is an absolutely crucial, and physically quite consequential, aspect of quantum gravity and which can be easily lost by an overly crude regularization procedure. Without any loss of generality it is therefore entirely legitimate to set the bare cosmological constant λ 0 = 1 in units of the cutoff [40].
The latter contribution then controls the scale for the edge lengths, and thus the overall scale in the problem. 3 It is clear by now that to accurately study the physical consequences of the theory requires the full machinery of quantum field theory and the renormalization group. Nevertheless, some key information about the behavior of physical correlations can already be obtained indirectly from averages of local diffeomorphism invariant operators. Also, it will often be convenient to continue to use the continuum language (as opposed to the lattice one) to discuss such quantities; in most cases the two languages are interchangeable, with the lattice one providing a more precise and thus less ambiguous (short-distance regulated) expression. Consider for example the average local curvature The above quantity describes the parallel transport of vectors around infinitesimal loops and is, by construction, manifestly diffeomorphism invariant. An appropriate lattice transcription reads [38,40] A second quantity of physical interest is the fluctuation in the local curvature The latter is directly related to the invariant curvature correlation function at zero momentum [40,42], see later below. On the lattice the previous quantity takes on the form Moreover, in the functional integral formulation of Eqs. (1) and (4) the average curvature R(G) and its fluctuation χ R (G) can also be obtained by taking derivatives with respect to k of the lattice partition function Z L in Eq. (4). On the lattice one has from the definition of the path integral as well as Exact scaling relationships then arise between various quantities, such as the ones in Eqs. (15) and (16), and these can later be used to derive scaling relations and check for mathematical consistency.

Diffeomorphism Invariant Gravitational Correlation Functions
Generally in a quantum theory of gravity physical distances between any two points x and y in a fixed background geometry are determined from the metric Because of quantum fluctuations, the latter depends, in the lattice case, on the specific edge length configuration considered. Correlation functions of local operators need to account for this fluctuating distance, and as a result these correlations are computed at some fixed geodesic distance between a given set of spacetime points [40,42]. On a given lattice this process involves constructing a complete table of distances between any two lattice points, and then computing from it the required two point functions. In addition, in gravity one generally requires that the local operators entering the correlation function should be coordinate scalars. In principle one could also smear such operators over a small region of spacetime with an assigned linear size [43,44].
It is then possible to also consider nonlocal gravitational observables, in analogy to what is done in Yang-Mills theories, by defining the gravitational analog Wilson loop. The latter carries information about the parallel transport of vectors around large loops, and therefore about large scale curvature [45,46,43,47,48,49], and will be discussed later.
In a quantum theory of gravity a fundamental two-point correlation function is the one associated with the scalar curvature, with physical points x and y separated by a given fixed geodesic distance d. On the lattice it has the corresponding form For the curvature correlation at fixed geodesic distance one expects at short distances (i.e. distances much shorter than the gravitational correlation length ξ to be introduced below) a power law decay with the power law here characterized by a universal exponent n. 4 On the other hand, for sufficiently strong coupling (large G, or small k) fluctuations in different spacetime regions largely decouple: the kinetic or derivative term in Eqs. (1) or (4) is responsible for coupling fluctuations in different spacetime regions, and in the action it comes with a coefficient 1/G. In this regime one then expects a faster, exponential decay, controlled by a nonperturbative correlation length ξ So the fundamental gravitational correlation length ξ can be defined unambiguously by the longdistance decay of the connected invariant curvature correlations at fixed geodesic distance d. Then the behavior in Eq. (20) is expected to hold at short distances d ≪ ξ, whereas the behavior in Eq. (22) is expected to hold at much larger distances, d ≫ ξ. In either case, in order to reach a sensible lattice continuum limit the physical distances involved need to be much larger than the fundamental average lattice spacing l 0 , d, ξ ≫ l 0 (the so-called scaling limit).
Consistency between the two expressions in Eqs. (22) and (20) is eventually regained from the fact that in the vicinity of the critical point a superposition of many exponentials are expected to add up to a power. This is seen, for example, from the spectral (Lehmann) representation of the two point function, with spectral function ρ(µ) = A µ α−1 /Γ(α). Then In the limit of a small infrared cutoff m ≡ 1/ξ the above result simplifies to a power law plus small corrections, In the last expression the known value for the gravitational curvature correlation function, α = 2 (4 − 1/ν) = 2 [44], has been inserted. 5 This simple observation allows one to compute the exponents ν and n more easily (and much more accurately) from the above expression than say from the distance-dependence of the correlation An important and central feature of the lattice nonperturbative treatment is the existence of a critical point in G, located at G c . The latter is interpreted as corresponding to a non-trivial fixed point in renormalization group language, see for ex. [44] and references therein. Furthermore it is known that the weak coupling phase G < G c is nonpertubatively unstable on the lattice: it corresponds to a branched polymer phase with no sensible continuum limit [38,40] (it is generally understood that such instabilities are usually quite difficult, if not impossible, to detect in a perturbative, or weak field, treatment). In accordance with this important result, in the following only the physical strong gravity phase for G > G c will be considered further. 5 In [44] the most recent numerical results ν = 0.334(4) and thus 2n = 2.01(7) are given. 6 One can contrast this power with what one obtains in weak field perturbation theory < √ gR(x) √ gR(y) >c ∼ < ∂ 2 h(x) ∂ 2 h(y) > ∼ 1/|x − y| d+2 , which is quite different from the result in Eq. (20) with n = 4 − 1/ν = 1, unless ν = 2/(d − 2), which is nevertheless correct for d close to two, where Einstein gravity becomes perturbatively renormalizable, and corrections to free field behavior become small.
In general, in the vicinity of such a nontrivial fixed point, one expects for the fundamental correlation length ξ = 1/m a power law divergence with Λ = 1/a the inverse lattice spacing, A ξ the correlation length amplitude, G c the critical point in the bare coupling G, and ν a universal exponent characterizing the divergence of ξ at the critical point. At the fixed point G c the theory regains scale invariance (due to the divergence of ξ), and in its vicinity one can then reconstruct the original, regularized continuum theory. In some ways ξ −1 = m can be viewed as a nonperturbative renormalized mass, analogous to the dynamically generated (but nevertheless gauge invariant) scale in Yang-Mills theories. For extensive reviews on the general subject of renormalization group scaling see, for example, [7,9,8,10,11]. There is by now a rather well established body of knowledge in quantum field theory and statistical field theory on this subject, and thus no obvious or apparent reason why its basic tenets should not apply to gravity as well, with quantum gravity describing essentially the unique theory of a massless spin two particle coupled to a covariantly conserved energy-momentum tensor. 7 One consequence of the renormalization group scaling relations in the vicinity of the fixed point, such as the scaling behavior [52] for the singular part of the free energy is to allow a precise determination of the correlation length exponent ν in Eq. (26) and associated quantities, such as amplitudes and corrections to scaling. From the lattice one finds for the critical value of G G c ≡ 1 8πk c = 0.623042 (25) .
and ν = 0.334(4) which is consistent with the conjectured exact value ν = 1/3 for pure quantum gravity in four dimensions [44]. After restoring dimensions, this in turn fixes the lattice spacing a, 7 It is a well-established fact that for theories with a nontrivial fixed point [1,2], the long distance (and thus infrared) universal scaling properties are uniquely determined, up to subleading corrections to exponents and scaling amplitudes, by the (generally nontrivial) scaling dimensions obtained via renormalization group methods in the vicinity of the fixed point [7,8,9,10,11]. These sets of results form the basis of universal predictions for, as an example, the perturbatively nonrenormalizable nonlinear sigma model [12,13]. The latter gives one of the most accurate tests of quantum field theory [14,15], after the g − 2 prediction for QED (for a comprehensive set of references, see [8,28], and references therein). It is also a well-established fact of modern renormalization group theory that in lattice QCD the scaling behavior of the theory in the vicinity of the asymptotic freedom fixed point unambiguously determines the universal nonperturbative scaling properties of the theory, as quantified by physical observables such as hadron masses, vacuum gluon and chiral condensates, decay amplitudes, the QCD string tension etc. [50,51]. and thus the value for the cutoff (in four dimensions G has dimensions of a length squared), 8 From the known laboratory value of Newton's constant G, l P ≡ G/c 3 = 1.616199(97) × 10 −33 cm one then obtains for the fundamental lattice spacing a = 1.2669 √ G c ≡ l P , or a = 2.0476 × 10 −33 cm , and from it a value for the cutoff Λ ≃ 1/a. This last result then allows one to restore the correct dimensions in all dimensionful quantities. 9 Note that from Eqs. (20), (22) and (25) one has The last scaling result follows from the fact that the curvature fluctuation is also the second derivative of the free energy with respect to k [see Eq. (16)], and that for the free energy the standard scaling assumption [52] in the vicinity of the ultraviolet fixed point reads F sing (G) = (27)] with ξ(G) given in Eq. (26). This then allows the fundamental exponent ν to be computed much more easily, and more accurately, than from the distance- whereas, from Eq. (25), the corresponding result for the curvature fluctuation is also quite simple in four dimensions, d = 4. The above results are rather helpful in establishing a direct connection between the correlation length ξ on the one hand, and the average local curvature R and its fluctuation χ R on the other hand.. 8 That the physical G is actually very close to Gc will be discussed later below. The argument involves in a key way the large scale curvature, and thus the quantum gravitational Wilson loop. 9 Note that in general the edge lengths are fluctuating and their average is close to, but not equal to, one. Nevertheless (for λ0 = 1) one finds for the average lattice spacing in units of a < l 2 > ≡ l 2 0 = [ 2.398(9) a ] 2 , so that a and l0 are quite comparable. Figure 1 and Table I present a detailed comparison between the lattice value for the universal exponent ν, and other approaches. The latter include the calculation of ν in the framework of the 2 + ǫ expansion for gravity in the continuum [32,33,34] carried out to two loop order [35], which . Note that the scaling exponent ν is expected to be universal, and therefore characteristic of quantum gravity (the unique theory of a massless spin two particle in four dimensions [22]), and therefore independent of specific features of the regularization procedure (lattice, dimensional regularization, momentum cutoff etc.). The same does not apply to the critical point and to the critical amplitudes, which are generally regularization dependent. 10 Another popular approach to the calculation of the universal exponent ν is based on a truncated renormalization group approach in the continuum in four dimensions. This gave values initially around ν −1 ≃ 2.8 [53,54] with some sizeable uncertainties; it is beyond the scope of this work to go into details regarding the features of each one of these calculations, so only a few representative cases will be mentioned here. Recent improved functional renormalization group calculations tend to generally fall roughly in the region ν −1 ≃ 2.0 to 3.5. Studies using a bi-metric parametrization gave ν −1 ≃ 4.7 [55], and later ν −1 ≃ 3.6 in [56,57]. In [58,59] it was argued that only fluctuations should be included that have an on-shell meaning, in which case one finds ν −1 ≃ 3.0, much closer to the lattice results. In [60,61] systematic studies were done of the dependence of the exponent ν on the metric parametrization and its influence on the functional measure contribution, giving generally for the leading exponents ν −1 ≃ 4, 2 to lowest order. A similar value ν −1 ≃ 3.0 was found using a geometric flow in the linear approximation in [62]. Another systematic large parameter space investigation of gauge fixing and measure choices was done in [63], with estimates eventually falling within the above mentioned range ν −1 ≃ 2.0 . . . 3.5.
The graph in Figure 1 also includes the known exact result for quantum gravity in three spacetime dimensions, obtained from the exact solution of the Wheeler-DeWitt equation in 2 + 1 dimensions, which gives ν −1 = 3/2 exactly [29,30,31]. The latter universal exponent should be compared to the old numerical Euclidean lattice result in three dimensions ν −1 = 1.72(10) [68], to the 2 + ǫ result of ν −1 = 1 + 3/5 + · · · = 8/5 = 1.60 [35], and finally to the Einstein-Hilbert truncation results mentioned previously, which in three dimensions cluster around ν −1 ≈ 2.3 [64] and ν −1 ≈ 1.6 . . . 2.0 [62], again in general agreement with the trend found for the lattice results 10 In statistical field theory y = ν −1 is usually referred to as the leading thermal (as opposed to magnetic) exponent. Under a real space renormalization group transformation with scale b one has for the corresponding relevant operator O ′ = b y O. The results presented here point to the existence of a single relevant operator in the vicinity of the ultraviolet fixed point, so that the corresponding operator O is associated, as expected, with the local scalar curvature.  70) and (62)] as a function of spacetime dimension d. Shown are the results in 2 + 1 dimensions obtained from the exact solution of the lattice Wheeler-DeWitt equation [30,31], the numerical result in four dimensions [44], the 2 + ǫ expansion result to one [34] and two loops [35], and the large d result ν −1 ≃ d − 1 [65]. For actual numerical values see Table I. in the same number of dimensions.
Other lattice and continuum methods can be used to provide an estimate for the exponent ν  [47], which gives ν = 1/2 + · · · . And, finally, the result of [66,67], where it was found that a consistent exact solution to the nonlocal effective field equations of Eq. (72) (discussed here later on) for the static isotropic metric in d ≥ 4 can only be found provided ν = 1/(d − 1) exactly, in agreement with the geometric argument mentioned earlier.

Renormalization Group Running of Newton's G
The results discussed so far are helpful in establishing a direct connection between the fundamental gravitational correlation length ξ and various diffeomorphism invariant averages such as the average Method used to compute the exponent ν in d=4 Universal Exponent ν Euclidean Lattice Quantum Gravity [44] ν −1 = 2.997(9) Perturbative 2 + ǫ expansion to one loop [34] ν −1 = 2 Perturbative 2 + ǫ expansion to two loops [35] ν −1 = 22/5 = 4.40 Einstein-Hilbert RG truncation [54] ν −1 ≈ 2.80 Recent improved Einstein-Hilbert RG truncation [58] ν −1 ≈ 3.0 Geometric argument [65] ρ vac pol (r) ∼ r d−1 ν −1 = d − 1 = 3 Lowest order strong coupling (large G) expansion [47] ν −1 = 2 Nonlocal field equations with G( ) for the static metric [66] Comparison of estimates for the universal gravitational scaling exponent ν, based on a variety of different analytical and numerical methods. These include the numerical results of [44], The 2 + ǫ expansion for pure gravity carried out to one and two loops [34,35], an estimate for the leading exponent in a truncated renormalization group expansion [54,58], a simple argument based on geometric features of the quantum vacuum polarization cloud for gravity [65], and finally the value obtained from consistency of the exact solution to the nonlocal field equation with a G( ) for the case of the static isotropic metric [66,67].
local curvature and its fluctuation. In this framework one can view the result of Eq. (26) as equivalent to stating that the Callan-Symanzik renormalization group beta function has a nontrivial zero at G c . Generally the cutoff independence of the nonperturbative mass scale m = 1/ξ Moreover, if one defines the dimensionless function F (G) via then, from the usual definition of the Callan-Symanzik beta function β(G) = ∂G(Λ)/∂ log Λ, one It follows that the renormalization group β-function, and thus the running of G(µ) with scale, can be defined some distance away from the nontrivial fixed point; more generally, the running of G(µ) is obtained by solving the differential equation with β(G) obtained from Eq. (36). Integrating Eq. (37) close to the nontrivial fixed point one with m 0 an integration constant of the RG equations. It has dimensions of a mass or inverse length, so it is naturally identified with the invariant correlation length ξ : m 0 ∝ 1/ξ. In particular, comparing results in Eqs. (38) and (26) one obtains which implies that the universal exponent ν is directly related to the derivative of the Callan-Symanzik β function in the vicinity of the fixed point at G c ; computing ν determines the universal running of G in the vicinity of G c . In addition, the renormalization group equations generally imply that the effective coupling G(µ) will grow (anti-screening) or decrease (screening) with distance scale r ∼ 1/µ, depending on whether G > G c or G < G c , respectively. One crucial physical insight obtained from the lattice is that only the phase G > G c is physically acceptable [44]; the phase G < G c corresponds at large distances to an entirely unphysical, collapsed branched polymer with no sensible continuum limit.
From the previous discussion one infers that the physical mass scale m = ξ −1 also determines the magnitude of the corrections to scaling, and plays therefore a role similar to the scaling violation parameter Λ M S in QCD. As in gauge theories, this nonperturbative mass scale emerges dynamically in spite of the fact that the fundamental gauge boson remains strictly massless to all orders in perturbation theory, and consequently its mass does not violate any local gauge invariance. Furthermore one expects, as in gauge theories, that in gravity the magnitude of ξ cannot be determined perturbatively, and to pin down a specific value requires a fully nonperturbative approach, as given here by the lattice formulation. In turn, the genuinely nonperturbative physical mass parameter m = ξ −1 of Eq. (26) is itself a renormalization group invariant and thus scale independent. In the immediate vicinity of the fixed point it obeys the general renormalization group equation, which follows from Eq. (34), with µ an arbitrary momentum scale. Here again, by virtue of Eq. (26), the second expression on the right-hand-side is only appropriate in very close proximity of the fixed point at G c . Solving explicitly Eq. (40) for G(q 2 ), with q an arbitrary wave vector scale, one finally obtains for the running of Newton's G with the action of Eqs. (2) or (5) Here again m = 1/ξ, and the coefficient c 0 for the amplitude of the quantum correction is (3) from a numerical study of the decay of curvature correlation functions, and also as before ν = 1/3 [42,44]. Consquently the dimensionless amplitude for the leading quantum correction in the lattice running of Eq. (41) is c 0 ≈ 8.02. This then completely determines the running of G in the vicinity of the fixed point, namely for scales r ≪ ξ.
Note that in the lattice theory of gravity only the smooth phase with G > G c exists (in the sense that an instability develops and spacetime collapses onto itself for G < G c ), which then implies that the gravitational coupling can only increase with distance [+ sign for the quantum correction in Eq. (41)]. In other words, a gravitational screening phase does not exist in the lattice theory of quantum gravity. The above situation appears to be true both for the Euclidean theory in four dimensions, and in the Lorentzian version in 3 + 1 dimensions [30]. A better, covariant formulation for the running of Newton's G is given later, in Eq. (70). Note also that the domain of validity for the expressions in Eq. (41) is q ≫ m ≡ 1/ξ or r ≪ ξ; the strong infrared divergence at q ≃ 0 is largely an artifact of the current expansion, and should be regulated either by cutting off the momentum integrations at q ≃ m = 1/ξ, or by the replacement on the r.h.s.
It is clear that the magnitude of the quantum correction in Eq. (41) depends crucially on the magnitude of the nonperturbative physical scale ξ. It will be argued later that this quantity is related, as in Yang-Mills theories, to the gravitational condensate, physically represented by the observed cosmological constant. Therefore at this stage it will turn out that the only physically sensible interpretation is that the observed λ obs is tentatively related to the scale ξ, From the above perspective "short distances" are not really that short, since ξ in comparison to G or the Planck length is a very large quantity, of cosmological magnitude. 11 It follows that the reference scale for the running of G in Eq. (41) is set by a correlation length ξ which, by Eqs. (32) 11 The fundamental nonperturbative scale ξ plays a crucial role in the following, and having a precise quantitative value for it is of paramount importance when trying to make contact with current astrophysical and cosmological observations. Here, for concreteness, a specific value in M pc will be assumed, in line with the most recent satellite data, see for example [69]. It is nevertheless quite possible that significant updates to this value will take place in the next few years, as increasingly sophisticated data, and data analysis methods, become available. One would nevertheless expect that various predictions, arising from the vacuum condensate framework described here, should lead to one single consistent value for the scale ξ. In this context it is worth remembering that before 1999 astrophysical observations were deemed to be entirely consistent with λ = 0. and later (88), is related to the observed large-scale curvature. In particular, the specific form for the running of G with scale suggests that no detectable corrections to classical gravity should arise either a) until the scale r approaches the very large (cosmological) scale ξ or b) until one reaches extremely short distances comparable to the Planck length r ∼ l p (at which point higher derivative terms, light matter corrections, and string contributions come into play). In other words, the results of Eq. (41) [or later in the covariant form of Eq. (70)] would imply that classical gravity is largely recovered on atomic, laboratory, solar, and even galactic scales, or as long as the relevant distances satisfy r ≪ ξ.

Gravitational Wilson Loop and Curvature Condensate
In gauge theories the Wilson loop is known to play a central role: on the one hand it is a manifestly gauge invariant quantity, on the other hand it provides key physical information on the nature of the static potential between two quarks. In gravity it is possible to construct a close analog of the gauge Wilson loop, by taking the path-ordered product of rotation matrices (describing the parallel transport of a vector, and thus specified in terms of the affine connection) along a closed loop. Nevertheless this path ordered product is not related to the gravitational potential; the latter is obtained from a different set of observables which involve the correlations of particle world lines modacorr,modaloop. Instead the gravitational Wilson loop provides information, as already in the infinitesimal loop case, on the behavior of curvature on very large scales.
The required integration over rotation matrices (or, equivalently, the integration over the affine connection) is most easily done in a first order formulation, where the affine connection and the metric are considered as independent degrees of freedom. Such a formulation exists on the lattice [70] and is therefore most suitable for computing the gravitational Wilson loop [47]. As in the gauge theory case, the integration over rotation matrices is performed using an invariant Haar measure over the group, which then almost immediately leads to a (minimal) area law for the quantum gravitational Wilson loop, Note that in the above expression use has been made of the fact that the basic reference scale appearing in the area law is the correlation length ξ , a well-known scaling result in gauge theories and justified there by renormalization group arguments. Also, C denotes the closed path that defines the loop; a more precise definition of the gravitational loop [47] will be given further below.
Suffice it to say here that the use of the Haar measure over rotations assumes and implies large local fluctuations in the metric, and thus in the affine connection, which is certainly justified for large G, where gravitational fluctuations in different spacetime regions decouple.
On the other hand, a macroscopic semiclassical observer is led to relate the parallel transport of a coordinate vector around a very large closed loop, via Stoke's theorem, to the value of the locally measured curvature. This then leads immediately to the semiclassical result [47] W where R is a measure of the slowly varying local macroscopic curvature; again a more precise definition will be given further below. Comparing the quantum result of Eq. (43)  This last set of considerations in turn provides a further key ingredient in quantum gravity, namely the correspondence between the macroscopic semiclassical curvature and the invariant correlation length ξ. One immediate consequence is that the scale for quantum effects in Eq. (41) is related to the observed cosmological constant, which in quantum gravity acts effectively as an infrared regulator. Thus potentially serious infrared divergences associated with the masslessness of the graviton are regulated by this new nonperturbative scale ξ, a mechanism which is similar to the way infrared divergences regulate themselves dynamically in QCD and non-Abelian lattice gauge theories. Consequently the scale ξ plays a role which seems analogous to the scaling violation parameter Λ M S in QCD; one important difference is that the running of G, due to the existence of a nontrivial fixed point, is not logarithmic. Instead the correct scale dependence of G is given by Eq. (41) and thus follows a power law, with an exponent ν related to the derivative of the beta function at the fixed point in G.
A second crucial consequence is that the scale for quantum effects is not given by Newton's constant; it is given instead by the size of ξ, which because of its relationship to the cosmological constant is a very large, cosmological scale of the order of 10 28 cm. It would seem therefore that such quantum effects will only become detectable when one explores the nature of gravity on cosmological scales comparable to ξ. The running of G is exceedingly tiny on solar system and galactic scales, but nevertheless increases dramatically as one approaches distance scales which are comparable to the observed cosmological constant λ. What then remains to be done is therefore to incorporate the above running of G into a set of generally covariant equations which can then be applied to the calculation of quantum corrections to known classical gravity results at very large distances. This will be discussed later. 12 It is important at this stage to understand where the Wilson loop relationship in Eqs. (43) and (43) is coming from. A precise definition of the gravitational Wilson loop was given in [43,45,47].
First note that infinitesimal transport loops appear already, for example, in the definition of the correlation function for the scalar curvature, Eq. (18). Here what will be considered instead is the parallel transport of a vector around a loop C which is not infinitesimal. In the following this loop will be assumed to be close to planar, a well-defined geometric construction described in detail in [47]. First define the total rotation matrix U(C) along the path C via a path-ordered (P) exponential of the integral of the affine connection Γ λ µν , The lattice action itself already contains contributions from infinitesimal loops, but more generally one might want to consider near-planar, but noninfinitesimal, lattice closed loops C. To make the above expression well defined it needs to be put on a lattice. There one defines a finite product of elementary rotations defined along a given lattice path The introduction of such rotation matrices in the Regge-Wheeler lattice was discussed in detail in [26,27,47], and a first order lattice formulation for gravity based on it was given in [70]; the following discussion will be based this well understood formalism. A coordinate scalar can then be defined by contracting the above rotation matrix U(C) with a unit length area bivector ω αβ (C), representative of the overall geometric features of the loop. Now if the parallel transport loop in question is centered at the point x, then one can define the operator W C (x) by with the near-planar loop centered at x and of linear size r C . Of course for an infinitesimal loop, involving an infinitesimal lattice path C 0 of linear size ∼ a, the overall rotation matrix is given by 12 Note that in gauge theories the correlation length ξ can be determined directly numerically by investigating the decay of Euclidean correlation functions of suitable local operators as a function of the separation distance. Generally these correlation functions are dominated by the lightest particle with a given spin. In the case of gravity such a detailed and complete analysis has not been performed yet, although it is in principle feasible, just as it is in lattice QCD. One complication that arises in the case of gravity is the fact that correlation functions between invariant operators have to be computed at a fixed geodesic distance [42]. w' w d Figure 2: Correlation function of two infinitesimal parallel transport loops, separated by a geodesic distance d. This correlation corresponds to the one defined in Eqs. (18) and (19). In the strong coupling limit one needs, in order to get a non-zero correlation, to fully tile the minimal tube connecting the two infinitesimal initial and final loops.
where now ω µν (C 0 ) is the area bivector associated with the infinitesimal loop of area ∼ a 2 , and δ the corresponding deficit angle; here R is lattice Riemann tensor at the hinge (triangle) in question, and A µν the corresponding area bivector. Then an invariant correlation function between two such operators is given by with the two loops separated by some fixed geodesic distance d. Of course for infinitesimal loops one recovers the expressions given earlier in Eqs. (18) and (19).
In general one needs to specify the relative orientation of the two loops. So, for example, one can take the first loop in a plane perpendicular to the direction associated with the geodesic connecting the two points, and the same for the second loop; the parallel transport of a vector along this geodesic will then be sufficient to establish the relative orientation of the two loops. Nevertheless if one is interested in the analog (for large loops) of the scalar curvature, then it will be adequate to perform a weighted sum over all possible loop orientations at both ends. This is in fact precisely what is done for infinitesimal loops of size r C ∼ a, if one looks carefully at the way the Regge lattice action was originally defined.
It is possible to give a more quantitative description for the behavior of the loop-loop correlation function given in Eq. (49), at least in the strong coupling limit. The following estimate is based on the previous results and definitions, and is further illuminated by the important analogy and correspondence of lattice gravity to lattice non-Abelian gauge theories outlined in detail in [65,47].
First it will be convenient to assume that the two (near planar) loops are of comparable shape and size, with overall linear sizes r C ∼ L/(2π) and perimeter P ≃ L. In addition, the two loops are separated by a distance d ≫ L, and for both loops it will be assumed that this separation is much larger than the lattice spacing, d ≫ a and L ≫ a. Then to get a nonvanishing correlation in the strong coupling, large G limit it will be necessary to completely tile a tube connecting the two loops, due to the geometric minimal area law arising from the use of the uniform (Haar) measure for the local rotation matrices at strong coupling, again as discussed in detail in [47]. Quite generally in this limit one expects an area law for the correlation between gravitational Wilson loops (see also the discussion for the Wilson loop itself given later below), which here takes the form with A(L, d) the minimal area of the tube connecting the two loops. Consistency of the above expression with the corresponding result for small (infinitesimal) loops given in previously in Eq. (22) requires that for small loops (small L) the value of L saturates to ξ, L ≃ ξ, so that the correct exponential decay is recovered for small loops.
This result is not unexpected, since ξ can only come into play only for distances much larger than the fundamental lattice spacing a. Consequently the asymptotic decay of correlations for large loops is somewhat different in form as compared to the decay of correlations for infinitesimal loops, with an additional factor of ξ appearing for large loops; nevertheless in both cases one has the expected minimal area law. In other words, the results of Eqs. (20) and (22) only apply to infinitesimal loops, which probe the parallel transport on infinitesimal (cutoff) scales; these results then need to be suitably amended when much larger loops, of semiclassical significance, are considered.
The above result applies to strong coupling, G ≫ G c . As one approaches the critical point at G c more than one exponential will contribute, in analogy to Eq. (20) for the single plaquette correlation. If the single loop contribution is proportional, as in the area law of Eq. (50), to exp(−m 2 Ld) with m = 1/ξ, then for a spectral function ρ(µ) = 2 B µ β−1 /Γ(β/2) one obtains, in the limit of a small infrared cutoff m ≡ 1/ξ,  In the strong coupling limit one needs, in order to get a non-zero correlation, to fully tile the tube connecting the two large initial and final loops.  20) and (24) given earlier for infinitesimal loops; in both cases the dependence on the separation is Eq. (20), the power law "short distance" (d ≪ ξ) behavior of the same correlation function.
One crucial ingredient needed in pinning down the magnitude of the quantum correction for G(q 2 ) in Eqs. (41) or (70), as well as the result for the loop correlation function of Eq. (53), is the actual value of the genuinely nonperturbative reference scale ξ. It was argued in [47] that, in analogy to ordinary gauge theories, the gravitational Wilson loop itself provides precisely such an insight. The main points of the argument are rather simple, and can thus be reproduced in just a few lines. In analogy to the gauge theory case, these arguments rely generally on the concept of universality, the existence of a universal correlation length at strong coupling, and the use of the Haar invariant measure to integrate over large fluctuations of the metric, or of the fundamental local parallel transport matrices. Following [45,46], in [47] the vacuum expectation value corresponding to the gravitational Wilson loop is naturally defined as Here the U 's are elementary rotation matrices, whose form is determined by the affine connection, and which therefore describe the parallel transport of vectors around a loop C; see also Eq. (46).
Again here ω µν (C) is a constant unit bivector, characteristic of the overall geometric orientation of the loop, giving the notion of a normal to the loop. In the continuum the combined rotation matrix U(C) is given by the path-ordered (P) exponential of the integral of the affine connection Γ λ µν , as in Eq. (45), so that the previous expression represents a suitable discretized and regularized lattice form. It can then be shown [47] that quite generally in lattice gravity, and for sufficiently strong coupling, one obtains universally an area law for large near planar loops where A C is the geometric minimal area of the loop as spanned by a given perimeter. 13 This last result relies on a modified first order formalism for the Regge lattice theory [70], in which the lattice metric degrees of freedom are separated out into local Lorentz rotations and tetrads.
Moreover, the result of Eq. (54) is in fact rather universal, since it can be shown to hold in all known lattice formulations of quantum gravity at least in the strong coupling (large G) regime.
In [47] an explicit expression for the correlation length ξ appearing in Eq. (55) was given in the strong coupling limit. There one finds ξ = 4/ k c | log(k/k c )| + O(k 2 ). For k close to k c this then gives immediately ξ ≃ 4 |k c − k| −1/2 and thus, to this order, ν = 1 2 in Eq. (26). Nevertheless, the discussion of the previous sections and the numerical solution of the full lattice theory suggests that the correct expression for ξ to be used in Eq. (55) should be the one in Eq. (26), with ν = 1/3 [Eq. (28)], k c given in Eq. (28) and amplitude A ξ = 0.80(3).
One then needs to make contact between these results and a semiclassical description, which requires that one connects the nonperturbative result of Eq. (55) to a suitable semiclassical physical observable. By the use of Stokes's theorem, semiclassically the parallel transport of a vector round a very large loop depends on the exponential of a suitably coarse-grained Riemann tensor over the loop. In this semiclassical picture one has for the combined rotation matrix U where A λσ is an area bivector, dx λ x σ . The above semiclassical procedure then gives for the loop in question Here ω µν (C) is a constant unit bivector, characteristic of the overall geometric orientation of the parallel transport loop. For a slowly varying semiclassical curvature, the R contribution can be taken out of the integral, so that the remaining integral depends on the overall large loop with some minimal area A C , for a given perimeter C. Then, by directly comparing coefficients for the two area terms in Eqs. (55) and (57), one concludes that the average large-scale curvature is of order +1/ξ 2 , at least in the strong coupling limit [47]. Since the scaled cosmological constant can be viewed as a measure of the intrinsic curvature of the vacuum, the above argument then leads to an effective positive cosmological constant for this phase, corresponding to a manifold which behaves semiclassically as de Sitter (λ > 0) on very large scales [47].
The above arguments then lead to the following key connection between the macroscopic (semiclassical) average curvature and the nonperturbative correlation length ξ of Eqs. (22), (41), (62) and (70), namely at least in the strong coupling (large G) limit. It is important to note that the result of Eq. (58) applies to parallel transport loops whose linear size r C is much larger than the lattice spacing, r C ≫ a; nevertheless in this limit the answer for the macroscopic curvature in Eq. (58) becomes independent of the loop size or its minimal area [47]. Furthermore, these arguments lead, via the classical field equations, to the identification of 1/ξ 2 with the observed (scaled) cosmological constant λ obs , 14 In this picture the latter is regarded as the quantum gravitational condensate, a measure of the vacuum energy, and thus of the intrinsic curvature of the vacuum. It is nonzero as a result of nonperturbative graviton condensation.
The above considerations can finally contribute to providing a quantitative handle on the physical magnitude of the nonperturbative scale ξ. From the observed value of the cosmological constant (see for ex. the 2015 Planck satellite data [69]) one obtains a first estimate for the absolute magnitude of the scale ξ, Irrespective of the specific value of ξ, this would indicate that generally the recovery of classical GR results happens for distance scales much smaller than the correlation length ξ. 15 In particular, the Newtonian potential is expected to acquire a tiny quantum correction from the running of G [see Eq. (70)] For example, in the case of the static isotropic metric one finds that G(r) is given explicitly by [67] G with m ≡ 1/ξ, so that quantum effects become negligible on distance scales r ≪ ξ. 16 One might think perhaps that the running of G envisioned here might lead to observable consequences on much shorter, galactic length scales. That this is not the case can be seen, for example, from the following argument. For a typical galaxy one has an overall size ∼ 30 kpc, giving for the quantum correction the estimate, from Eq. (62) for the static potential, (30 kpc/5320 × 10 3 kpc) 3 ∼ 1.79 × 10 −16 which is tiny due to the large size of ξ [see Eq. (60)]. It seems therefore unlikely that such a correction will be detectable at these scales, or that it could account, in part, for anomalies in the galactic rotation curves. The above argument nevertheless shows a certain 14 Up to a constant of proportionality, expected to be of order unity. 15 The value for λ, and therefore ξ, relies on a multitude of current cosmological data, which nowadays is usually analyzed in the framework of the standard ΛCDM model. Included in the usual assumptions is the fact that Newton's G does not run with scale. If such an assumption were to be relaxed, it would affect a number of cosmological parameters, including λ, whose value could then change significantly. In the following the estimate of Eq. (60) will be used as a sensible starting point. 16 The quantum gravity correction is reminiscent of the Uehling term in QED; nevertheless the latter is purely logarithmic, and the infrared cutoff there is provided by the smallest mass scale appearing in QED loop corrections, the renormalized electron mass. In quantum gravity the role of the infrared cutoff is played by the graviton mass, which in perturbation theory (as in Yang-Mills theories) stays strictly zero to all orders, due to local coordinate invariance. sensitivity of the results to the value of the scale ξ; thus an increase in ξ by a factor of two tends to reduce the effects of G( ) by 2 3 = 8, as can be seen from Eq. (70) with ν = 1/3 and the fact that the amplitude of the quantum correction is always proportional to the combination c 0 /ξ 3 .
Here again m is related to the gravitational correlation length via m ≡ 1/ξ, see Eqs. (59) and (60). One notes therefore that the Hartree-Fock approximation to the self-consistent equation for the graviton vacuum polarization tensor also predicts an infrared rise of G(q) (antiscreening), and furthermore unambiguously determines the amplitude of the quantum correction (c 0 , here equal to 3/2). In this approximation the mean field result for the exponent is ν = 1/(d − 2), so that the power is equal to two for Eq. (63) in four dimensions. So there are two main differences that stand out compared to the lattice result of Eq. (41), namely that the power is two here and not three, and the fact that here there is an additional, slowly varying log(q) component. 17 The above results also suggest that the curvature on very small scales behaves rather differently from the curvature on very large scales, due to the quantum fluctuations eventually averaging out.
Indeed when comparing the result of Eq. (32) to the one in Eq. (58) one is lead to conclude that the following change has to take place when going from small (linear size ∼ l p ) to large (linear size ≫ l p ) parallel transport loops An intuitive way of understanding the above result is that on small scales the strong local fluctuations in the metric/geometry lead to large values for the average rotation of a parallel-transported vector. But then on larger scales these short distance fluctuations tend to average out, and the 17 One of the earliest applications of the Hartree-Fock approximation to solving Dyson's equations for propagators and vertex functions was in the context of the BCS theory for a superconductor. Later it was applied to a (perturbatively non-renormalizable) relativistic theory of a self-coupled Fermion, where it provided the first convincing evidence for a dynamical breaking of chiral symmetry and the emergence of Nambu-Goldstone bosons [74]. combined overall rotation is much smaller, by a factor of O(l p /ξ), The above quantity should then be regarded as an essential and necessary "renormalization constant" when comparing curvature on different length scales, and specifically when going from very small (size ∼ l P ) to large (size ≫ l P ) parallel transport loops. See also the earlier discussion preceding Eq. (50), about the issue of comparing correlations of large loops versus correlations of small (infinitesimal) loops.
To conclude this section, one can raise the legitimate concern of how these results are changed by quantum fluctuations of various matter fields; so far all the results presented here apply to pure gravity without any matter fields. Therefore here, and in the rest of the paper, what has been applied is basically the quenched approximation, wherein gravitational loop effects (perturbative and nonperturbative) are fully accounted for, but matter loop corrections are initially neglected.
When adding matter fields coupled to gravity (scalars, fermions, vector bosons, spin-3/2 fields etc.) one would expect, for example, the value for ν to change due to vacuum polarization loops containing these fields. A number of arguments can be given though for why these effects should not be too dramatic, unless the number of light matter fields is rather large [44]. 18 Note that the above results for the gravitational condensate in many ways parallel what is found in non-Abelian gauge theories, where for example one has for the color condensate < F 2 µν > ≃ 1/ξ 4 [77,78,79,80]. In QCD this last result is obtained from purely dimensional arguments, once the existence of a fundamental correlation length ξ (which for QCD is given by the inverse mass of the lowest spin zero glueball) is established. Accordingly, for gravity one would in fact expect simply on the basis of purely dimensional argument that the large scale curvature (corresponding to the graviton condensate) should be related to the fundamental correlation length by R ≃ 1/ξ 2 , as in Eq. (32). This then points to a fundamental relationship between the nonperturbative scale ξ (or inverse renormalized mass) and a nonvanishing vacuum condensate for both of these theories, nonperturbative quantum gravity and QCD, In gauge theories an additional physically relevant example is provided by the fermion condensate, arising as a non-trivial consequence of the renormalization group, confinement and chiral symmetry breaking in SU(3) gauge theories [81]; for a recent review on current values see for ex. [82].
Note that in all three cases the power of ξ is fixed by the canonical dimension of the corresponding field, one over length square in the case of the curvature, an observation which can be seen to provide further support to the identification in Eqs. (58) and (59) and ψ ψ ≃ (290 M eV ) 3 [79,80]. On the other hand, modifications to the static potential in gauge theories are best expressed in terms of the running coupling constant α S (µ), whose scale dependence is determined by the celebrated beta function of QCD. There the relevant scale is the nonperturbative Λ M S ≈ 210 M eV whose size is comparable to ξ, Λ M S ≃ ξ −1 . More specifically, in 18 One would expect that significant changes to the result of Eqs. (41) and (70) will arise from matter fields which are light enough to compete with gravity, and whose Compton wavelength is therefore comparable to the scale of the gravitational vacuum condensate, or observed cosmological constant λ, namely m −1 ∼ 1/ λ/3. At present the number of candidate fields that could fall into this category is rather limited, with the photon and a near-massless gravitino belonging to this category. The results presented here correspond to the quenched approximation for quantum gravity, where all graviton loop effects are included, but matter (and radiation) loops are neglected. Matter field are still present in the theory but are treated as quantum mechanical static sources. In the 2 + ǫ perturbative expansion for quantum gravity one encounters factors of 25 − c in the renormalization groups β function, where c is the central charge associated with the (massless) matter fields [34,35], which would suggest that matter loop and radiation corrections are indeed rather small. In four dimensions similar factors involve 48 − c [75,76], which would again lend support to the argument that such effects should be rather small in four dimensions.

Effective Field Equations
The result of Eq. (41) expresses the renormalization group running of Newton's G as a function of momentum scale. As it stands, the expression in Eq. (41) does not satisfy general covariance, and needs to be promoted to a more useful and acceptable form. It follows that in order to apply consistently the above result to an arbitrary background geometry, a fully covariant formulation is required. One way of describing the running of Newton's G is by a set of effective nonlocal field equations with a G( ) [66,41]. A second option is to formulate a fully covariant effective gravitational action with a running G( ), also discussed in detail in [66,41]. In the following both options will be discussed.
An effective field theory approach can be derived by writing down an effective action, involving either a G(µ) or a G( ) [66,41]. Based on the previous discussion, a suitable effective action, describing the residual effects of quantum gravity on very large distance scales, is of the form with G(µ) a very slowly varying (on macroscopic scales) Newton's constant, in accordance with Eqs. (41) or (70), and amplitude c 0 ≈ 2 × 8.02(55) ≈ 16.0 [44]. Note that the effective action of Eq. (68) is obtained from the one in Eq. (2) by a suitable field rescaling, in accordance with the discussion preceding Eq. (9). Here again the nonperturbative scale ξ appears therefore both in the running of G and in the cosmological constant term with 1 3 λ = 1/ξ 2 . Nevertheless it was found that if a covariant G( ) is used in the above effective action the resulting effective field equations are rather complicated and hard to solve in practice, due to the fractional exponents appearing in G( ) [66,41].
For reasons that will become clearer later on, in most of the upcoming discussion a different route will be followed, based on the much simpler (and thus more manageable) effective field equations based nevertheless again on G( ). In either case it seems natural to perform the standard quantum mechanical replacement q 2 → − , where (g µν ) is the covariant D'Alembertian for a given background metric g µν (x) [66], This then leads to a consistent covariantly formulated running of G, with Note that the precise form of the covariant , and thus of G( ), depends on the tensor nature of the object it acts on [66]. Numerical studies of lattice quantum gravity give for the exponent ν = 1/3 and for the quantum amplitude c 0 ≈ 16.0 [44], which (fortunately, or unfortunately) leaves very little ambiguity regarding the result of Eq. (70). As noted earlier, one way of viewing physically the result of Eq. (70) is that quantum gravitational fluctuations generate anti-screening, with an initially very slow running of G, as shown earlier in Figure 4. The anti-screening arises because of the radiative dressing of the source by a virtual gravitation graviton cloud, in analogy to the screening of a bare charge in QED by the virtual electron-positron cloud. In a sense, therefore, the above corrections describe the gravitational analog of the running coupling constant in QED.
Generally, fractional powers of inverse d'Alembertians require careful handling. This can be done either by computing the effect of integer powers n and then analytically continue the result to fractional negative values n → −1/2ν, or by using a regulated parametric integral representation where µ → 0 is a suitable infrared regulator, here again with exponent ν = 1/3. Note that for ν = 1/3 the quantum correction in Eq. (70) proportional to c 0 always includes a 1/ξ 3 , which therefore naturally sets the overall scale for the leading quantum correction, irrespective of the background geometry considered. Then a suitable set of manifestly covariant effective field equations with a running G( ) takes the form [66] proportional to the amplitude c 0 . Already the action of (g) = g µν ∇ µ ∇ ν is rather complicated on second rank tensors; one has Of course, one of the simplest applications is to the Friedmann-Lemaître-Robertson-Walker (FLRW) framework applied to the standard homogeneous isotropic metric in the following only the case k = 0 (spatially flat universe) will be discussed. In this framework a popular choice for T µν is the perfect fluid form, for which one needs to compute the action of n on T µν , and then analytically continues the answer to negative fractional values of n = −1/2ν. The results of [66,67,83] then show, among other things, that a nonvanishing pressure contribution is generated in the effective field equations, even if one initially assumes a pressureless fluid, p(t) = 0. Specifically, for a universe filled with nonrelativistic matter (p=0) one obtains the following set of effective Friedmann equations, for the tt field equation, and k a 2 (t) +ȧ 2 (t) for the rr field equation. In the above expressions, the running of G appropriate for the Robertson-Walker metric is with c t ≃ 0.450 c 0 [for an amplitude c 0 appearing in Eq. (70)] for the tensor box operator [66]. From the above form of δG(t) one sees that again the amplitude of the quantum correction is proportional to the combination c 0 /ξ 3 for ν = 1/3. Furthermore, the running of G induces an effective pressure term in the second (rr) equation, due to the presence of an induced relativistic fluid, whose origin lies in the quantum gravitational vacuum-polarization contribution. Another noteworthy feature of the new effective field equations is the additional power-law acceleration contribution, on top of the standard exponential one due to the λ term.
On way of viewing the results is that the effective field equations with a running G, here Eqs. (76) and (77), can be recast in an equivalent form by defining a vacuum-polarization pressure p vac and density ρ vac , such that in the FLRW background one has From this viewpoint, the inclusion of a vacuum-polarization contribution in the FLRW framework amounts to a replacement ρ(t) → ρ(t) + ρ vac (t), p(t) → p(t) + p vac (t) in the original field equations.
Then, just as one introduces a parameter w, describing the matter equation of state, with w = 0 for nonrelativistic matter, one can do the same here for the remaining quantum contribution by setting The original calculations [66], and more recently [84] (which included metric perturbations) give w vac = 1 3 . Note that it was shown in [66] that this result is obtained generally for the given class of G( ) considered, and is not tied to a specific choice for the universal exponent ν, such as ν = 1 3 . More generally, the procedure of defining a ρ vac and a p vac contribution, arising from quantum gravitational vacuum-polarization effects, is not necessarily restricted to the FLRW background metric case. One can always decompose the full source term in the effective nonlocal field equations of Eqs. (70) and Eq. (72), making use of as two contributions, The latter then involves the nonlocal part Consistency of the full nonlocal field equations now requires that the sum be covariantly conserved, In general one cannot expect that the contribution T vac µν will always be expressible in the perfect fluid form of Eq. (75), even if the original T µν for matter (or radiation) has such a form. The former will in general contain, for example, nonvanishing shear stress contributions, even if they were originally absent in the matter part [84]. Indeed in a number of cases of physical interest one deals quite generally with a background metric that is slightly perturbed, g µν =ḡ µν + h µν .
Consequently the covariant d'Alembertian operator = g µν ∇ µ ∇ ν acting on second rank tensors [such as the T µν in Eq. (72)] needs to be Taylor expanded in the small perturbation h, For G( ) itself one obtains the expansion

Large Scale Curvature and Matter Density Correlations
Quantum gravity, and the existence of a nontrivial quantum condensate for the curvature, lead to a number of specific physical predictions, which are in principle observationally testable. Many of these quantum correction effects can in principle be calculated, given the effective long-distance quantum corrected gravity theory formulated in Eqs. (72) and (70) Much of what has been discussed so far was relates to the fact that in a quantum theory of gravity the gravitational constant G runs with scale, in accordance with Eq. (70). But there are additional consequences, which arise from the fact that in general gravitational correlations do not follow free field (Gaussian) predictions and which will therefore be the subject of this section. One would expect such correlations to have some observational relevance, and one such example is the curvature correlation function of Eqs. (18) and (20), with power 2n = 2(d − 1/ν) = 2. But it will be useful to first examine some local averages. For the average local curvature R(k) one has from Eq.(32), using Eq.(26) and ν = 1/3, as given earlier in Eq. (32). Lattice calculations allow one to extract various amplitude coefficients, such as the one in the above expression. The dimensionless amplitude A ′ R in Eq. (88) is expected to be O(1) in lattice units, and numerically one finds [44] A ′ R = 3.40 (13). On the other hand for the curvature fluctuation χ R (k) one has from Eqs. (13) and (33) also given earlier in Eq. (33). Note that in both Eqs. (88) and (89) the correct dimensions have been restored, by inserting suitable powers of the lattice spacing a (curvature has dimensions of one over length squared); also a specific value for this lattice spacing was given earlier in Eq. (30). For the dimensionless amplitude in Eq. (89) one finds numerically A ′ χ = 2.22(9) [44], again in general agreement with the prejudice that nonvanishing dimensionless critical amplitudes should be O(1).
These results in turn provide some useful information related to the local curvature correlation function at a fixed geodesic distance [see Eqs. (18) and (19)]. By scaling one obtains immediately (90) 20 One can contrast this result with what one finds in weak field perturbation theory. There one finds [45] 6 and thus 2n = 6, so the result here is quite different. If one defines in the usual way an anomalous dimension η for the graviton propagator in momentum space, < h h > ∼ 1/k 2−η , one finds from the lattice calculation η = d − 2 − 2/ν or η = −4 in four dimensions for ν = 1/3, which deviates For the local curvature-curvature correlation function of Eq. (18) at "short distances" r ≪ ξ (and again for ν = 1/3) one then obtains the rather simple result with C 123 a constant, and relative geodesic distances d ij = |x i − x j | etc. The relevance and measurement of nontrivial three-and four-point matter density correlation functions in cosmology was discussed in detail some time ago in [87].
It is instructive at this stage to compare the above result for the local curvature correlation given in Eq. (91) to the expression for the local average curvature of Eq. (88): both expressions still contain explicitly the size of the microscopic, infinitesimal parallel transport loop ∼ a ∼ l P , which originates in the fact that both these quantities make explicit reference to infinitesimal (ultraviolet cutoff sized) parallel transport loops. The explicit appearance of the ultraviolet cutoff in these averages can be explained by the appearance of residual short distance divergences associated with such small infinitesimal loops. At the same time, a comparison of the result of Eqs. (88) and (88) for the local curvature with the corresponding result of Eq. (58) for the large scale, macroscopic curvature suggests a substantial changeover when going from small (size ∼ a) to large (size ≫ a) parallel transport loops. As discussed earlier, the results can be summarized by the statement that the curvature on large (macroscopic) scales is much smaller (by a factor 1/ξ) than the curvature on small (Planck length) scales, due to a dramatic averaging out of the fluctuations.
As described earlier in Eqs. (64) and (65), one then expects that the transition from infinitesimal to macroscopic loops (linear size ≫ a) can be affected in the correlation function of Eq. (91) by the rather significantly from the Gaussian or perturbative value. Such a large deviation is already observed in the 2 + ǫ expansion [see Eq. (39)] and is not peculiar to lattice quantum gravity; in the context of gravity such an interesting possibility was already discussed some time ago in [86].
replacement of a 2 → ξ 2 . This then would give for large (macroscopic size ≫ a) parallel transport loops a modified form of the correlation function of Eq. (91) with the general expectation that the overall amplitudes nevertheless be comparable, Note that the universal power n = 2 is unchanged compared to Eq. (91), only the amplitude has been modified, in accordance with the earlier gravitational Wilson loop result of Eqs. (52) and (53). Unfortunately so far these large loop correlations have not been computed explicitly, but nevertheless the above ideas should become testable by explicit numerical simulations in the near future.
The next step is to determine whether the knowledge of the curvature correlation, as given in Eqs. (91) or (93), can be translated into information regarding other two-point correlations subject to astrophysical measurement. First consider what can be stated purely at the classical level.
There one can use the field equations to directly relate the local curvature to the local matter mass density. From Einstein's field equations with λ = 0 for a perfect fluid one then obtains for the Ricci scalar, in the limit of negligible pressure, This last result then relates the local fluctuations in the scalar curvature δR(x) to local fluctuations in the matter density δρ(x), and could therefore provide a potentially useful connection to the quantum result of Eqs. (20), (91) and (93). Note that the same kind of reasoning would apply alternatively to a T µν describing radiation, which would then be relevant for a radiation-dominated early universe. Of course, in the Newtonian limit the above result simplifies to Poisson's equation where h 00 = 2φ and ρ are the macroscopic gravitational field and the macroscopic mass density, respectively.
Now, in the current cosmology literature [87,89] it is customary to describe matter density fluctuations in terms of the matter density contrast correlation function 21 G ρ (r) = < δρ(r) δρ(0) > .
The latter is related to its Fourier transform P (q) by G ρ (r) = 1 2π 2 Λ µ dq q 2 P (q) . sin qr qr , It has to contain, in general, both an infrared regulator (µ) and an cutoff (Λ), to make sure the integral stays convergent. If one assumes, as is sometimes the case, that the power spectrum P (q) is described by a simple power law (where n s = −s is commonly referred to as the spectral index), then one finds in the scaling regime 1/µ ≫ r ≫ 1/Λ for the real-space density contrast correlation function provided the index s satisfies 0 < s < 3, and here c s ≡ Γ(2 − s) sin(πs/2) / 2π 2 (terms containing the ultraviolet cutoff Λ appear as well, but they are proportional to sin(Λr) and cos(Λr), oscillate rapidly and average out to zero). This then leads to the identification of exponents [see Eq. (90) and Eq. (91) For the specific value s = 1 one has P (q) = a 0 /q, and therefore in position space which is in fact, as will be discussed later, also consistent with the result given earlier in Eq. (91) for the invariant, connected curvature correlation function.
Generally, all of the above expressions get modified for very small wave vector q ∼ m, as should be clear by now from the discussion in the previous sections, where the appearance of a dynamically generated infrared cutoff (as in QCD, and more generally in non-Abelian gauge theories) stood out as an essential ingredient. Such an infrared cutoff is either introduced explicitly in the wave vector integrations, or, alternatively, the power spectrum is infrared regulated at small q by replacing One then writes more appropriately, instead of Eq. (99), Again, by Fourier transforms one then obtains for the correlation in real space which reproduces Eq. (100) for short distances r ≪ µ −1 . For very large spatial separations r ≫ µ −1 the asymptotic behavior of G ρ (r) is now given instead by (105) with amplitude c ′ s = 1/2 1+ s 2 Γ( s 2 ). In view of the previous discussion, it is natural to identify here the infrared cutoff µ ≡ m = 1/ξ, and such a choice will be implicit from now on throughout the following discussion. Note the correspondence of the result of Eq. (100) with the short distance curvature correlation result of Eq. (20), as well as the same type of correspondence of Eq. (105) with the large distance curvature correlation result of Eq. (22).
In practice, observational data for such matter density correlations is commonly presented in the following compact form [87] G ρ (r) = r 0 r γ , with an empirically determined exponent γ, and a scale r 0 fitted to astrophysical (usually galactic cluster) observations. For γ close to two, one has by comparing Eq. (102) to Eq. (106) a 0 = 2π 2 r 2 0 . It seems therefore rather tempting at this stage to try to connect the observational result of Eq. (106) to the quantum correlation function in Eq. (91). One then expects for the matter density fluctuation correlation a power law decay as well, of the form 23 where a(t) here stands for the scale factor. 24 The last correlation function can be made dimensionless by suitably dividing it by the square of some average matter density ρ 0 ≈ 0.3089ρ c with ρ c = 3H 2 0 /8πG and H −1 0 ≈ 4430 M pc for h = 0.677 (using again, for concreteness, the Planck 2015 data [69]). By comparing powers and coefficients in Eqs. (91) and (106) one finds For the reference scale r 0 in this last equation one derives the quantitative estimate 23 Note that in weak field perturbation theory one finds, by virtue of the equations of motion, < ρ(x) ρ(y) >c ∼ < ∂ 2 h(x) ∂ 2 h(y) > ∼ 1/|x − y| 6 , so again the result here is rather different. 24 The scale-factor dependent prefactor involving a(t) and a(t ′ ) only has such a simple form in a matter-dominated universe; for λ = 0 the a-dependent prefactor is a bit more complicated, nevertheless it still reduces to unity for equal times t = t ′ = t0 [89]. The preceding arguments still contains nevertheless a fundamental flaw, related to the use, at this stage in unmodified form, of the curvature correlation function result of Eq. (91). As discussed previously, that form applies to the correlation of infinitesimal (Planck length, or cutoff size) loops, which would not seem to be appropriate for the macroscopic (or semiclassical) parallel transport loops, such as the ones that enter the field equations (94) and (95), and which thus relate locally the macroscopic δR(x) to the δρ(x). It would therefore seem desirable to be able to correct for the fact that the parallel transport loops sampled in Eq. (95) are much larger than the infinitesimal ones sampled in the correlation function in Eq. (91). As in Eqs. (64) and (65), the transition to macroscopic loops (linear size ≫ a) can be affected in Eq. (91) by the replacement of a 2 → ξ 2 .
This then gives for large (macroscopic size ≫ a) parallel transport loops [see Eq. (93)] with the expectation of a comparable amplitude A 1 ≈ A 0 . This last result then leads to the following improved estimate for the macroscopic matter density correlation of Eq. (97), so that comparing to Eq. (106) one finds again for the exponent γ = 2, and for the length scale r 0 the improved value which now seems more in line with observational data. For the Fourier amplitude a 0 in Eq. (99) one obtains the estimate Observed galaxy density correlations give indeed for the exponent in Eq. (106) a value close to two, namely γ ≈ 1.  72), which involve a running G( ), or more simply make use of G(q) as given in Eq. (41). When the above replacement is performed, one finds with as before m −1 ≡ ξ and (in view of the preceding discussion) still s = 1. Then the most important modification, as can be seen quite clearly in Figure 5, is the rather dramatic decrease in magnitude (due to the 1/G(q) 2 factor) of P (q) for small q, with a clear turnover at In view of the more complex behavior of P (q), it is clearly no longer possible to associate a single spectral index with P (q). Following Eq. (99) one can nevertheless define an effective spectral index s(q) via This quantity tends to s = 1 for q ≫ 1/ξ but then dips below zero for q ≃ 1/ξ. A plot of s(q) is shown in Figure 6. One would expect that such a drastic turnover, caused by the quantum running of G(q) on very large scales, should become visible in future cosmological observations. From a practical perspective, it might make more sense to treat the numerical amplitude A 1 in Eq. (93) as a free parameter, given the uncertainties, discussed earlier, associated in its determination from a first principle lattice calculation [see the discussion following Eq. (110)]. In other words, it would seem that so far, based on existing results, the q-dependence of P (q) is more credible at this stage than its overall normalization.  99) and (103) for a constant G. Also, the spectral exponent is s = 1 and the amplitude is a 0 , as discussed in the text [see Eqs. (101) and (113)]; in the plot the q wave vector is measured for convenience in units of ξ. Note the rather marked turnover for small q ≈ 4.20/ξ due to the running of G, as discussed in the text. The nonperturbative scale ξ is related to the gravitational vacuum condensate, as in Eqs. (59) and (58).
In conclusion, the main results of this section can be summed up as follows. The vacuum condensate picture of quantum gravity leads to three main predictions for matter density correlations, of which the first one is that the power appearing in Eq. (106) should be exactly γ = 2 for "short distances" r ≪ ξ (or q ≫ 1/ξ), and that the reference length scale r 0 appearing in the same equation should be related to ξ, as in Eq. (112). The second prediction is that the power spectrum exponent in Eq. (99) should be exactly s = 1 again for q ≫ 1/ξ, and that the amplitude a 0 in the same equation should be related to ξ as in Eq. (113). The third prediction is that the power spectrum P (q) as a function of q should exhibit a marked break for q ∼ ξ −1 , as given later in Eqs. (103) and (114), and shown in Figures 5 and 6. 25 More generally, the results outlined here and in the previous sections suggest that existing, and future, astrophysical and cosmological data should be re-analyzed in terms of a wider G(q) scale dependent form, similar to the one in Eq. (41), and of the general (but nevertheless rather simple) two-parameter type G(q) with q −1 0 ≫ 1 M pc a wave vector reference scale, and p a positive power. 25 A further complication arises in the cosmological context from the fact that at late times H −1 0 ≡ȧ/a ≃ t0 ≃ 0.79 ξ so that H −1 0 and ξ are quite comparable in magnitude (which follows from the FRW evolution of the scale factor for a universe dominated by a λ term), whereas at early times a new, much shorter scale appears H −1 ≡ȧ/a ≪ ξ. Such effects are expected to play a role as well, and have not been taken into account here yet. Figure 7: Effective spectral index s(q) as defined in Eqs. (99) and (115). The horizontal line at the top is the s = 1 value, corresponding to a constant (scale-independent) Newton's G. The spectral index approaches the value s = 1 for large q ≫ 1/ξ, but dips below zero for q ≃ 1/ξ.

Gravitational Slip Function with G( )
A running of Newton's G gives rise to significant long-distance effects, which fundamentally modify the classical field equations of general relativity at very large distance scales. The following section provides an update on the results originally presented in [85], especially in view of the recent high accuracy lattice results presented in [44], and in particular related to the new improved estimate for the quantum amplitude c 0 of Eq. (70). It is common practice to describe relativistic effects in cosmology within the framework of the comoving gauge, where the metric is written as with background metricḡ µ ν = diag −1, a 2 , a 2 , a 2 , fluctuations such that h 0i = h i0 = 0, with the remaining h ij 's decomposed into a trace and stress part so that T r(h ij ) = a 2 h. In this gauge the metric is then parametrized by the scale factor a(t) and the two additional functions s and h.
There are nevertheless instances where effects which deviate from standard Newtonian physics are more easily described within the context of the conformal Newtonian gauge, where the metric is parametrized by two scalar potentials ψ and φ [99,100]. It is of some interest to explore possible cosmological consequences of a running Newton's constant G( ) in this gauge, as discussed recently in [84]. In this gauge one sets for the metric g µ ν =ḡ µ ν + h µ ν withḡ µ ν = a 2 diag (−1, 1, 1, 1), h 0i = h i0 = 0, and furthermore with the components of h µν (x) again considered here as a small perturbation. The line element is then given by Then gravitational slip function η(x) is defined as In classical General Relativity one has φ(x) = ψ(x) giving then η(x) = 0, which makes the quantity η a rather useful parametrization for deviations from the classical theory (whatever their origin might be).
Generally, quantum corrections arising from a running of G( ) give rise to additional terms in the field equations, which no longer ensure that φ = ψ. In practice, at some stage of the calculation one needs to compute higher order contributions from the h ij 's which requires one to expand, for example, G( ) in the metric perturbations h ij , Consequently the covariant d'Alembertian = g µν ∇ µ ∇ ν has to be Taylor expanded in the small field perturbation h µν , and similarly for G( ) as in Eq. (87), which requires the use of the binomial expansion for the where the superscripts (0) and (1) refer to zeroth and first order in this weak field expansion, respectively. It is also customary in these treatments to expand all relevant fields in spatial Fourier modes. One sets for the matter density, pressure and velocity, as well as for the metric, with q the comoving wave number, and similarly for any of the other fields as well. An additional approximation is then to limit at first the treatment to the leading q = 0 mode, and leave the more challenging O(q) corrections for a later calculation. Nevertheless here the nonlocal nature of the quantum corrections makes the calculation of the slip function η, as well as other quantities, technically rather difficult due to the intrinsic non-locality of G( ) in Eq. (70), which will require at some stage a number of physically motivated approximations, so that a partial, useful answer can be obtained. Furthermore, the general result for the slip function η is most easily presented in a form where the metric perturbation h ij is decomposed into the comoving trace (h) and a stress (s) part; then in the original comoving gauge the spatial metric perturbation is written as in Eq. (117).
It is useful here to record the (rather straightforward) relationship between perturbations in the standard comoving and conformal Newtonian gauge, namely Further details regarding the various choices of gauge (comoving, synchronous and conformal Newtonian) and their mutual relationship, as they apply to this specific problem, can be found in [85].
The zeroth order (in the metric perturbation h µν ) the Friedman equations in the presence of a running G are given in the comoving gauge by with w = 0 for non-relativistic matter and w vac = 1 3 for the graviton vacuum polarization contribution, together with the energy conservation equation These equations need to be solved first, in order to obtain an expression for the background scale factor a(t), and for the background average matter densityρ(t). Of course for δG(t) = 0 the above equations reduce to the standard Friedman equations with cosmological constant λ.
When restricted to the q = 0 mode (the zeroth order, dominant infrared mode) one obtains for the slip function an expression which is a function of G(t) or G(a), where t is the cosmological time and a(t) the corresponding scale factor. The detailed dependence of G on t or a(t) (and other parameters) is in turn quite sensitive to the specific choice of cosmological parameters, such as the rate of expansion, the relative contribution of dark matter versus dark energy, etc. In the end one finds the following q = 0 result for the slip function η [85] with δG(t) given in Eq. (78) and ν = 1/3.
A specific numerical value for η depends on various assumptions introduced in order to concretely evaluate the above expression. One of the simplest cases corresponds to the limit of a vanishing cosmological constant, λ ≃ 0. In view of Eq. (59), this last limit corresponds to a very large ξ. For a general perfect fluid with equation of state p = wρ one then has for the scale factor a(t) = a 0 (t/t 0 ) 2/3(1+w) and for average matter density ρ(t) = 1/[6πGt 2 (1 + w) 2 ]. Then from Eq. (127) one obtains for pure non-relativistic matter, w = 0, and more generally for w = 0 In both cases the amplitude c t = 0.45 c 0 for the tensor box operator [66], with c 0 entering the expression for G( ) of Eq. (70), and c 0 ≈ 16.04, giving therefore an overall coefficient c t ≈ 7.2.
Another extreme, but nevertheless equally simple, case is a pure cosmological constant term, which can be modeled by taking the limit w = −1.
To analyze the more general case of a non-vanishing cosmological constant λ = 0 combined with non-relativistic matter (w = 0), the following form [84] for the slip function expressed in terms of the scale factor a(t) turns out to be more useful The integration constant a ξ is fixed by the requirement that the scale factor a → a ξ for t → ξ [see Eqs. (70), (78) and (141) below for the definitions of a ξ ]. In other words, by switching to the variable a(t) instead of t, the quantity ξ has been traded for a ξ . In practice the quantity a ξ is generally expected to be slightly larger than the scale factor "today", i.e. for t = t 0 . As a result the correction in Eq. (130) is expected to be negative today. What is then needed is the general relationship between t and a(t) (for nonvanishing cosmological constant λ) so that a quantitative estimate for the slip function η can be obtained from Eq. (130). Specifically one is interested in the value of η for a current matter fraction Ω ≃ 0.31, as suggested by current astrophysical measurements (according to the recent Planck 2015 data, see [69]). To obtain δG(a) one then makes use of δG(t) from Eq. (78), and then substitutes the correct (matter-fraction dependent) relationship between t and a [85], which among other things contains the constant with the parameter θ defined as  to the extreme limiting case of comoving wavevector q ≃ 0 . In analogy to the quantum mechanical particle (the graviton) in a box, the lowest possible mode corresponds to q ≃ π/L where L is the linear size of the box. Here the role of the L is played by the correlation length ξ, so that the lowest possible mode corresponds in fact to q ≃ 1/ξ, for which then the result of Eq. (133) applies.
Generally the quantum contributions to the slip function are scale dependent, and thus will be proportionately reduced if one looks at scales which are significantly smaller compared to the largest scale in the problem, namely ξ. The calculation of the slip function η(q) over a wider range of scales is clearly a significantly more complicated problem, which has not been addressed yet; as stated previously all calculations performed so far [85] only correspond to limit of q → 0 , for which the above analytical estimate has been obtained. Nevertheless, if the corresponding relevant astrophysical length scale is denoted by l 0 , then at such a scale one must have, simply by scaling from Eqs. (41) and (70), with η(t = t 0 , q = 0) given earlier in Eq. (133). In other words, the result for η is reduced by the ratio of the relevant length scale compared to t 0 or ξ, to the third power since the exponent ν = 1/3. As a practical example, if one were to look at the value of the slip function η(q) on scales which are an order of magnitude less than the reference scale ξ, then this would reduce, in accordance with Eqs. (41) and (70)  usually restricted in a first approximation to the lowest comoving spatial momentum q modes.
In this limit the growth parameter δ(t) ≡ δρ(t)/ρ obeys, as a function of scale a(t), a rather simple ordinary differential equation, whose solution then provides, given suitable initial conditions, information about the matter and dark energy content of the current universe. One quantity that is often brought into play is the growth index f (a), namely the derivative of the log of δ(a) with respect to the log of the scale factor a(t), and in addition the parameter γ = log f / log Ω, which provides information on how the growth index f (a) depends on the current matter fraction Ω [87].
Cosmological observation suggests that today's matter fraction is about Ω ≈ 0.31 [69], leading to a value of γ = 0.55, based pretty much entirely on what is obtained from the systematic treatment of density perturbations within the framework of classical General Relativity.
It follows that many of the calculations just described can be repeated if one assumes now that Newton's constant runs with scale, so that the standard field equations of GR get modified by the non-local term of Eq. (70). Under the physically motivated assumption of a comparatively slowly varying (both in space and time) background, it is then possible to obtain a complete and consistent set of effective field equations, describing small perturbations for the metric trace and matter modes [84,85]. This then gives rise, within the same set of methods and approximations used in classical GR, to a set of equations for the growth amplitude. The latter are then studied again, initially, in the limit of small q wave vectors, and this in turn leads to modified growth exponents. In general the results are expected to be quite sensitive to the scale q, but so far only the leading term as q goes to zero has been calculated analytically, due to technical difficulties which arise from the strong non-locality of G( ). The following section provides a significant update on the results presented originally in [85], especially in view of the recent high accuracy lattice results presented in [44], and in particular the new improved estimate for the quantum amplitude c 0 of Eq. (70).
Besides the modified cosmic scale factor evolution due to the G(t) discussed earlier [see for ex.
Eqs. (77) and (76)] the running of G( ) as given in Eq. (70) also affects the nature of matter density perturbations on large scales. In computing these effects, it is customary to introduce a perturbed FRW metric of the form with a(t) the unperturbed scale factor and h ij (x, t) a small metric perturbation, and h 00 = h i0 = 0 by choice of coordinates. After decomposing the matter fields into background and fluctuation contribution, ρ =ρ + δρ, p =p + δp, and v =v + δv, it is customary in these treatments to expand the density, pressure and metric perturbations in spatial Fourier modes, as in Eq. (123) with q the comoving wave number. Then the field equations with a G( ) [Eq. (72)] are given, to zeroth order in the perturbations h ij , by the unperturbed field equations with a G(t), which in turn fixes the three background fields a(t),ρ(t), andp(t) in accordance with Eqs. (77) and (76).
At the next step, in order to obtain an equation for the matter density contrast δ(t) = δρ(t)/ρ(t), it is customary to eliminate the metric trace field h(t) from the field equations. This is first done a(t). Consequently, instead of using the expression for G(t) in Eq. (78), one uses the equivalent expression for G(a) Here the power is γ ν = 3/2ν for non-relativistic matter, since from Eq. (78) one has then a(t)/a 0 ≈ (t/t 0 ) 2/3 for constant G; in the following ν = 1 3 for which then γ ν = 9/2 for this case. If on the other hand one uses a more general equation of state of the form p = wρ then a(t)/a 0 = (t/t 0 ) 2/3(1+w) , and therefore γ ν = 3(1 + w)/2ν. Also, c a ≈ c t if a 0 is identified with a scale factor corresponding to a universe of size ξ; to a good approximation this corresponds to the universe "today", with the relative scale factor customarily normalized at that time t = t 0 to a(t 0 ) = 1. Furthermore, in [66] it was found that in Eq. More generally, the zeroth order tt field equation with constant G = G 0 can be written in terms of the current matter density fractions as with a/a 0 = 1/(1 + z) where z is the red shift and a 0 = 1 the scale factor today. In this last case H 0 is the Hubble constant evaluated today, Ω the (baryonic and dark) matter density, Ω R the space curvature contribution corresponding to a curvature k term, and Ω λ the dark energy or cosmological constant part, all again measured today. In the absence of spatial curvature k = 0 one has then Then in terms of the scale factor a(t) the equation for matter density perturbations for constant G = G 0 , Eq. (139), becomes The quantity H(a) is most simply obtained from the FLRW field equations which can in principle be solved for the scale factor a(t), leading to It is customary at this stage to introduce a parameter θ describing the cosmological constant fraction as measured today, In practice one is mostly interested in the observationally favored case of a current matter fraction Ω ≈ 0.25 [more recent data [69] suggest a slightly larger value of 0.31], for which then θ ≈ 3. In terms of θ the equation for the density contrast δ(a) for constant G can then be recast in the form with the growing solution to the above equation given explicitly by with c 1 a multiplicative constants and 2 F 1 (a, b; c, z) the Gauss hypergeometric function. The subscript 0 in δ 0 (a) means that the solution here is appropriate for the case of constant G = G 0 .
To evaluate the correction to δ 0 (a) coming from the terms proportional to c a from G(a) in Eq. (141) one sets where F(a) is a function to be determined, and then inserts the resulting expression in Eq. (137), written as a differential equation in the scale factor a(t). One only needs to write down the In this last expression the exponent is γ ν = 3/2ν ≃ 9/2 for a matter dominated background universe w = 0, and more generally γ ν = 3(1 + w)/2ν; even the use of the general Eq.
with the two coefficient functions A(a) and B(a) given by rather complicated functions [85].
The solution of the above differential equation for the matter density contrast in the presence of a running Newton's constant G( ) then leads to an explicit form for the function δ(a) = δ 0 (a) [1 + c a F(a)]. From it, an estimate of the size of the corrections coming from the new terms due to the running of G can be obtained. It is clear from the previous discussion, and from the form of G( ), that such corrections are expected to become increasingly important towards the present era t ≈ t 0 or a ≈ 1.
Specifically, in Ref. [84] a value for the density perturbation growth index parameter γ was computed in the presence of G( ). When discussing the growth of density perturbations in classical GR [87] it is customary at this point to introduce a scale-factor-dependent growth index f (a) defined where δ(a) is the matter density contrast discussed above. In principle, the latter is obtained from the solution to the general differential equation for δ(a), such as the one in Eqs. (148) or (152).
Nevertheless, one is mainly interested in the neighborhood of the present era, a(t) ≃ a 0 = 1, which leads to the definition of the growth index parameter γ via The latter has been the subject of increasingly accurate cosmological observations, for some recent references see for example [103,104,105]. The solution of the differential equation for δ(a) with a One has c a = (t 0 /ξ) 3 · 0.45 · c 0 , with t 0 corresponding to "today" so that t 0 /ξ ≈ 0.794, and c 0 = 16.04; the additional factor of 0.45 arises in relating the tensor G( ) in Eqs. (70) to the G(t) appropriate for the FRW background metric in Eq. (78), as computed in [66]. Then in Eq. (141) one has c a = 0.501 · 0.45 · 16.04 = 3.62, which gives a substantial overall amplitude. To quantitatively estimate the actual size of the correction in the above expressions for the growth index parameter γ, and make some preliminary comparison to astrophysical observations, some additional information is needed.
At first, one notices that all calculations done so far refer to the case of comoving wave number q = 0 in Eq. (123). If those numbers were used directly, one would obtain rather large O(1) quantum corrections to the growth parameter γ, where γ 0 is the classical GR value, and γ c the leading quantum correction in the limit q = 0 (which incidentally, in all cases looked at so far, turns out to be negative). To obtain corresponding results for q = 0 would then require a new, and significantly more complex, calculation which has not been done yet.
Nevertheless it seems clear that one can apply a simple scaling argument to obtain the more general result by a significantly shorter route. One notes that the quantum correction in Eq. (155) is, by virtue of the explicit form of G( ) in Eqs. (41) and (70), always proportional to the inverse of the nonperturbative reference length scale cubed, ∝ 1/ξ 3 .
Unfortunately, this issue later complicates considerably the problem of relating δG(t) to δG(a), and therefore the solution to the resulting differential equation for δ(a). However, as a tractable approximation, one can use in the interim the slightly more general result for the scale factor a(t) valid for any w = 0, namely a(t) = a 0 (t/t 0 ) 2/3(1+w) (the extreme case of a vacuum energy dominated cosmology, w = −1, is discussed in [84,85] as well). As an example we will use here an "effective" value of w ≈ −7/9, which would seem more appropriate for the final target value of a matter density fraction Ω ≈ 0.25. For this choice one then obtains a significantly reduced power in Eq. (141), namely γ ν = 3(1 + w)/2ν = 1. Then, although Eq. (137) for δ(t) remains unchanged, Eq. (152) for δ(a) need to be solved with new parameters. Furthermore, the resulting differential equation for δ(a), Eq. (152), is still relatively easy to solve, by the same methods discussed earlier.
For this case as stated γ ν = 1 in Eq. (141), and one obtains a somewhat smaller correction compared to the matter dominated case w = 0 of Eq. (156), namely So far a number of general features can be observed in the results, the first one being the fact that generally the quantum correction to the growth index parameter γ is found to be negative. . 27 Taking for these cases a reference scale l 0 = 10 M pc in Eq. (157) one obtains a correction to γ ≃ O(10 −6 ) which is rather tiny. It is therefore clear that the quantum effects discussed here are 27 For recent detailed reviews of the many tests of general relativity on astrophysical scales, and a more complete set of references, see for example [107,108]. only relevant for very large scales, much bigger than those usually considered, and well constrained, by laboratory, solar or galactic dynamics tests [106,107,108,110].
one notices that the latter correction is more than an order of magnitude larger. So it seems the bound from matter density perturbations is much more stringent than the one derived from the slip function.
Indeed the nonperturbative amplitude coefficient c 0 enters all calculations involving G( ) with the same magnitude and sign. One can therefore relate one set of physical results to another, such as the quantum correction to the slip function η(z = 0), given in Eq. (133), to the quantum corrections to the density perturbation growth exponent γ, given in Eq. (157). Then after taking the ratio the amplitude coefficient conveniently c 0 drops out, and one obtains for the ratio of the quantum corrections to the matter density perturbation growth parameter γ to the quantum slip function η for t = t 0 δ γ δ η ≃ + 16.3 .
This last result suggest that it will be observationally more difficult (by an order of magnitude) to see the quantum G( ) correction in the slip function η(q) than in the matter density growth parameter γ(q). Nevertheless, perhaps the value of the present calculations lies in the fact that so far a discernible trend seems to emerge from the results, and it suggests that the quantum correction to the growth exponent γ is initially rather small for small clusters, negative in sign, but slowly increasing in magnitude, following a cubic law with scale.

Conclusions
The vacuum condensate picture of quantum gravitation provides in principle a series of detailed and testable predictions, which could either be verified or disproved in the near future as new and increasingly accurate satellite observations become available. The previous sections laid out a number of specific prediction and estimates, many originating in a rather direct way from the gravitational scaling dimensions and amplitudes of invariant gravitational correlations, obtained originally from a variety of different nonperturbative approaches, including the Regge-Wheeler lattice formulation of gravity.
One key aspect linking all these calculations together is the fact that once the nonperturbative scale ξ is set (in analogy to the ΛM S of QCD) then, in accordance with the renormalization group, there are no further adjustable parameters when discussing the universal, long distance limit. The quantum theory of gravity is therefore, like the classical theory, again highly constrained by general coordinate invariance. In more than one way the current calculations are still rather incomplete; in particular the gravitational Wilson loop and the correlation between loops have not been studied in detail yet, and only some rather general properties have been inferred. Nevertheless the feeling is still that the underlying formulation is solid enough to allow future controlled and improved estimates of key renormalization group quantities.
The previous discussion has made it clear that the derivation of many of the basic results has relied heavily on subtle -but well established (and well grounded) -renormalization group scaling arguments. Nevertheless there is no reason yet to doubt that such arguments should be fully applicable to gravity as well. Particularly encouraging is the fact that by now four different approaches to quantum gravity give comparable results for the phase structure and scaling dimensions (see the comparison Table I, as well as Figure 1), which suggests a unique underlying renormalization group universality class associated with the massless spin two field in four dimensions. Of course one important common element in all these approaches is the existence of a non-trivial fixed point in G of the renormalization group equations. A key physical aspect that emerges from the theory is a growth of Newton's G with scale, quantified by a new, genuinely nonperturbative correlation length ξ, with the latter intimately related to the gravitational vacuum condensate (and thus to the physical, observed cosmological constant). The second key physical aspect is the existence of nontrivial scaling and anomalous dimensions for gravitational n-point functions, which leads, among others, to specific predictions for matter density correlations and related quantities. In conclusion, the main aspects of this new physical picture for gravity can be summarized as follows: • The vacuum condensate picture of quantum gravity contains from the start a very limited set of parameters, and is as a result strongly constrained. It involves a new, genuinely nonperturbative scale [the gravitational vacuum condensate, see Eqs. (58) and (59)], which relates the running of Newton's G to the current observed value of the cosmological constant, and to the long distance behavior of physical diffeomorphism invariant curvature correlations.
• While in principle both signs could be possible, in the strong coupling limit the effective, long distance cosmological constant is positive [see Eq. (58), the arguments preceding it, and more detailed discussion in [47]]. The basic argument relies of the behavior of the gravitational Wilson loop : in the same strong coupling regime it seems impossible from the lattice theory to obtain a negative value for the effective cosmological constant, irrespective of the choice of boundary conditions (which, incidentally, in the lattice context play no role in the argument).
• The theory predicts a slow increase in strength of the gravitational coupling when very large, cosmological scales are approached [see Eqs. (41) and (70)]. In this context the observed scaled cosmological constant λ is seen to act as a dynamically induced infrared cutoff, similar to what happens in non-Abelian gauge theories. In principle, both the universal power and amplitude for this infrared growth are calculable from first principles in the underlying lattice theory.
• The lattice theory appears to exclude the possibility of a physically acceptable phase with gravitational screening. The perturbative, weak coupling (small G) phase is found to be inherently unstable in the lattice formulation, a consequence of the conformal instability.
Thus a genuinely semiclassical regime for quantum gravity, whereby quantum effects can be included as small perturbations, does not seem to exist. On the other hand for large enough quantum fluctuations (large G) the conformal instability is overcome, and a new stable, antiscreening phase emerges. The stability of quantum gravity can thus be viewed as an entropy effect, intimately connected to non-trivial properties of the gravitational functional measure.
• Calculations presented here give a number of specific predictions for the behavior of invariant curvature correlations as a function of geodesic distance, and specifically the powers and amplitudes involved [see Eqs. (91) and (93)]. Perhaps the most important result is the fact that the curvature correlation function decays like inverse distance squared (n = 1 and thus s = 1 and γ = 2). This in turn can be used to relate in a standard way, via the quantum equations of motion, curvature correlations to matter density correlations and thus to their observed power spectrum [see Eqs. (102) and (107)].
• Given the exceptionally large value of the scale ξ (originating from the fact that the observed λ is very small compared to the scale associated with G) no observable deviations from classical General Relativity are expected on laboratory, solar systems and even galactic scales [see Eqs. (62) and (61)].
In addition, the picture outlined here points to what appears to be a deep analogy between the nonperturbative vacuum state of quantum gravity and known properties of strongly coupled non-Abelian gauge theories, and specifically QCD. Indeed in QCD one also finds a nonperturbative mass parameter m = 1/ξ (sometimes referred to as the mass gap) which is known to arise dynamically without nevertheless violating any local gauge symmetries, and is understood to be a renormalization group invariant. That such a mass scale can be generated dynamically is a highly nontrivial outcome of the strong coupling dynamics of QCD, and its associated renormalization group equations. Over time this analogy has been of great help in illustrating properties of quantum gravity, many of which are ultimately still based on fundamental principles of the renormalization group, connected with universal scaling properties as they apply to the vicinity of a nontrivial ultraviolet fixed point.