Non-perturbative propagators in quantum gravity

We employ non-perturbative renormalisation group methods to compute the full momentum dependence of propagators in quantum gravity in general dimensions. We disentangle all different graviton and Faddeev-Popov ghost modes and find qualitative differences in the momentum dependence of their propagators. This allows us to reconstruct the form factors quadratic in curvature from first principles, which enter physical observables like scattering cross sections. The results are qualitatively stable under variations of the gauge fixing choice.


Introduction
The unification of gravity with quantum mechanics is a notoriously hard problem in theoretical physics. Even a century after the development of quantum theory and general relativity, breakthroughs in the understanding of quantum properties of spacetime are few and far between. This is not least so because of the typical scale that we expect quantum gravity effects to be important at -the Planck scale, which is about 10 −35 m. To illustrate this fantastically small number, measuring the typical size of a human to an accuracy of a Planck length is roughly comparable to measuring the extension of the Milky Way with an accuracy of an atomic nucleus. This emphasises why experimental data on quantum gravity effects is hard to come by, and accordingly quantum gravity theories presently mostly rely on theoretical considerations, and can only be confronted with consistency tests.
A conservative strategy, and potentially a path to success in formulating a theory of quantum gravity, is to build on well-established theories valid at larger length scales, and only add as few extra assumptions as possible to extend the theory to include new phenomena. In the context of quantum gravity, such an approach is the Asymptotic Safety program [1]. It embraces the importance of symmetries in the understanding of the Standard Model and adds an interacting ultraviolet completion of gravity in the form of a quantum realisation of scale symmetry. An enormous advantage of this approach is its closeness to standard quantum field theory notions, and the straightforward connection to low energy physics. Especially in the matter sector, this allows to confront the scenario with many observational consistency tests. A disadvantage is that the existence of such a scale invariant regime is hard to prove. In practice, this is however less of a problem. Many interacting fixed points have been found with satisfactory precision in other contexts, e.g. in statistical physics and condensed matter systems [2][3][4][5]. By now, there is ample evidence that a suitable fixed point indeed exists not only in various approximations on the dynamics of pure quantum gravity , but also in the presence of matter . Phenomenological implications of the interacting fixed point have been discussed in the context of particle physics [50,61,70,71,[77][78][79][80][81][82][83][84][85], cosmology [40,[86][87][88][89][90][91][92][93], and black holes [94][95][96][97][98][99][100][101][102]. For recent reviews and introductions, see [103][104][105][106][107], and for a critical discussion of the status of the field, see [108,109]. Evidence for asymptotic safety in gravity has also been found in lattice formulations of quantum gravity, namely Euclidean and Causal Dynamical Triangulations [110][111][112][113][114][115][116][117][118]. On the lattice, a scale invariant regime can be realised if a second-order phase transition exists that gives rise to consistent low-energy physics.
The physical renormalisation group running of couplings is one of the key objects of study in quantum field theories. This translates into the momentum dependence of correlation functions, which are the basic building blocks of observables like scattering cross sections. The easiest non-trivial correlation function is the propagator, which is the inverse of the two-point correlation function. It stores important information about the unitarity and causality of the theory as it is related to the spectral function (if it exists), see e.g. [119]. An accurate description of the propagator is thus crucial to decide whether Asymptotic Safety provides a theory of quantum gravity which is compatible with these notions.
In this work we compute the full momentum dependence of the propagators of the graviton and its accompanying Faddeev-Popov ghost with non-perturbative renormalisation group techniques. A key advance is that we distinguish between all the different modes of the graviton and the ghost. As is well-known from representation theory, the graviton splits into a gauge-invariant spin two mode (the transverse-traceless mode), a gauge-invariant spin zero mode, and a pure gauge vector mode (which can be split into a spin one transverse, and a spin zero longitudinal component). The ghost is a vector field itself, and also splits into a transverse and a longitudinal mode. The physical information is stored in the spin two and gauge-invariant spin zero mode. The momentum dependence of their propagators can be mapped to the form factors which are quadratic in curvature tensors.
The central results can be summarised as follows: • The spin two and spin zero modes of the graviton feature qualitative and quantitative differences in their momentum dependence. • The overall gauge and gap dependence is small, see Figure 4 and Figure 5. • Within our approximation, only the spin two mode shows the property of momentum locality [31], and furthermore only in four dimensions and with equal three-and four-graviton coupling, see (65). • A derivative expansion of the form factors shows alternating signs, indicating that computations relying on such an expansion are inherently unstable and invariably introduce fictitious poles [120], see (85). • The propagators of the two ghost modes are related for asymptotic momenta, but differ for finite momenta, see (60) and (70). • Quantum corrections to the free propagator decrease exponentially with increasing dimension, see (78) and Figure 3.
This work is structured as follows. Our renormalisation group tool of choice is the so-called functional renormalisation group, which we briefly review in section 2. In section 3 we discuss our setup and some general aspects of momentum-dependent correlation functions. Section 4 and appendix A comprise a discussion of analytical properties of the momentum dependence of the propagators, whereas in section 5 we discuss numerical results. Sections 6 and 7 as well as appendix B contain a comparison of results obtained in different approximation schemes. Finally, we conclude with a summary and an outlook in section 8.

Functional Renormalisation Group
To extract the scale dependence of the system, we employ the functional renormalisation group (FRG). It is based on the Wetterich equation [121][122][123], describing the flow of the scale-dependent action Γ k , which was adapted to gravity in [6]. Here, we denote t = log k as the "RG-time", Γ (2) k is the second functional derivative of Γ k with respect to fluctuation fields, and the super-trace STr contains a sum over all fields and indices, as well as an integration over the continuous coordinates. The regulator R k , which enters the generating functional as a momentum-dependent mass term for the fluctuation, provides an infrared (IR) regularisation of modes with p 2 < k 2 , and thereby ensures IR finiteness. The factor ∂ t R k ensures ultraviolet (UV) finiteness by cutting off modes with p 2 > k 2 . Overall, the regulator term together with its derivative implement the Wilsonian idea of momentum-shell wise integration of quantum fluctuations. As a result, the scale-dependent effective action Γ k interpolates between the classical action S when no quantum fluctuations are integrated out, i.e. in the limit k → ∞, and the full quantum effective action Γ, when k → 0. A central role is played by the fixed points of the flow, which realise quantum scale invariance, where all couplings measured in units of the scale k are scale-independent. For reviews of the FRG, see [106,[124][125][126][127][128]. Due to the momentum-shell wise integration of quantum fluctuations, the Wetterich equation (1) is formulated on Euclidean backgrounds, since only then the squared momentum gives definite information on whether a certain mode is a UV or an IR mode. In the following, we will therefore assume a Euclidean background and hence discuss Euclidean quantum gravity. The generalisation of the functional RG to Lorentzian spacetimes is one of the important challenges of the approach [109], and first steps towards investigations of Lorentzian spacetimes have been presented in [20,23,37,[129][130][131][132][133][134].
Due to the formulation as a local coarse graining and the necessity of a gauge fixing, the formal introduction of a background metricḡ is hard to avoid. In this spirit, the full metric g is split into background metricḡ and a (not necessarily small) metric perturbation h according to Other parameterisations have been investigated for example in [33][34][35]53,63,129,[135][136][137][138][139][140][141][142][143][144][145]. The practical advantage of the background field method is that one can retain background diffeomorphism invariance at each step. Within the FRG, the background metricḡ only serves as a technical tool to allow for a momentum-shell wise integration of quantum fluctuations, and in principle never has to be specified. Without approximations of the dynamics of the theory, the physical results obtained at k → 0 therefore are entirely independent of the specific choice ofḡ. However, in the following sections, we will use a concrete choice forḡ, since this significantly reduces the computational complexity.
Since the scale-dependent effective action Γ k contains all operators that are consistent with the symmetries of the system, practical computations of RG flows of couplings require the truncation of Γ k to a manageable set of operators. Especially in gauge theories, the momentum dependence of scale-dependent correlation functions provides crucial information on the system, such as on unitarity. Thus, in approximations, the resolution of momentum dependence is unavoidable to obtain reliable results. A scale identification which allows to translate the RG running of operators into their momentum dependence is, in general, insufficient. For example, higher order n-point correlation functions generically depend on n different momenta independently, such that the momentum dependence cannot be captured by the dependence on just one scale k. It is however possible that for special momentum configurations, the RG running qualitatively agrees with the physical momentum dependence [119]. Furthermore, the RG scale k is an artificial scale introduced via the regulator R k to regularise the path integral. Physical quantities, such as scattering amplitudes, can only be extracted in the limit k → 0, where all contributions from the regulator R k vanish. In this limit, the UV behaviour of a given theory is described by the momentum dependence of operators. The momentum-dependent evaluation of β-functions is thus necessary to extract the UV physics of a system.
In the context of asymptotically safe quantum gravity, different expansion schemes are used, which allow to extract momentum-dependent flow equations. On the one hand, a vertex expansion in terms of metric fluctuations h is employed [24,29]. In this so-called fluctuation approach, see [107] for a review, the starting point is a seed action S, which is then expanded in terms of fluctuation fields h, typically around a flat background, see [43,69] for expansions around non-flat backgrounds. Consequently, the different vertices are labelled with scale-dependent couplings, whose flow is evaluated for different values of the external momenta, allowing to extract momentum-dependent flows of n-point correlation functions. In this approach, the momentum-dependent RG running of vertex correlation functions of the graviton two- [29], three- [31] and four-point functions [39], as well as three-point functions in gravity-matter systems [52,58,[64][65][66] have been investigated. As an important approximation, all computations in the fluctuation approach identify the scale dependence of different tensor structures of the graviton with the scale dependence of the purely transverse-traceless part of the corresponding n-point correlator. On the other hand, the form factor expansion [68,101] is based on diffeomorphism invariant operators, where the corresponding form factors are general functions of covariant derivatives. In this expansion, the scale dependence of the entire form factor can be extracted within the background field approximation. As a caveat, the background field approximation ignores that the regulator and the gauge fixing break the full diffeomorphism symmetry, which causes the scale-dependent effective action Γ k to depend on g and h individually [15,16,24,[29][30][31]39,42,43,52,[64][65][66]68,69,107,119,144,[146][147][148][149][150][151][152][153][154][155][156][157][158][159][160][161][162][163]. Comparing fluctuation and background results at the same level of approximation thus gives information about how much the respective modified Ward identities are broken.
In the present work, we provide a comparison of both expansion schemes on the level of the two-point function. For this purpose, we will set the stage and provide more details on momentum dependence in section 3. In sections 4 and 5, we derive and analyse the momentum-dependent flow equation for the graviton two-point correlator for general dimension d. Importantly, we distinguish all different tensor structures at the level of the two-point function. After gauge fixing, we compute the scale dependence of four different tensor structures for the graviton and ghost two-point functions. Additionally, in section 6, we compute the corresponding form factors on the background field level which contribute to the gravitational two-point function. We show that on this level in the expansion and under certain assumptions, there is a one-to-one correspondence between vertex correlation functions and form factors. Furthermore, in section 7 we compare both expansions based on the analysis of asymptotic behaviours of form factors and correlation functions.

Momentum dependence in quantum gravity
In this section we discuss our setup to resolve momentum-dependent correlation functions. We start with a discussion of the correlators themselves. Next, we present the general structure of the RG flow. Finally, we relate the form factor expansion to the fluctuation correlation functions.

Fluctuation approach
For the computation of vertex correlation functions, we choose a flat Euclidean background which is a technical choice simplifying the computations significantly. In principle, each diffeomorphism invariant operator in the seed action is labelled by a single coupling. However, under the presence of the regulator and gauge fixing, the scale dependence of different n-point functions, originating from the same operator in the seed action, generally differ [39,58,[64][65][66]. Therefore, we will introduce separate couplings for each operator in the vertex expansion. Furthermore, on the level of the two-point function, we will distinguish between the transverse-traceless (TT) mode, and the gauge-invariant scalar mode. Specifically, the seed action we will use in the following reads where we approximate S grav by the Einstein-Hilbert action S EH describing classical gravity: Here, G N and Λ are the Newton coupling and the cosmological constant, respectively. Computing the effect of quantum fluctuations of gravity requires the inclusion of a gauge fixing condition F µ = 0, where we choose Here,D µ refers to the background covariant derivative, and β h and γ h are gauge parameters. The presence of γ h allows to compute the scale dependence of all tensor structures on the level of the graviton two-point function. The gauge fixing condition is implemented into the action via the inclusion of the gauge fixing action The gauge parameter α h controls how strongly we implement the gauge fixing condition in the path integral. A strict implementation corresponds to the Landau limit α h → 0. The resulting Faddeev-Popov determinant is taken care of by introducing ghost fields c µ andc ν , with ghost action S gh : The Lie derivative L c g σκ of the full metric g µν along the direction of the ghost field is given by Round brackets indicate a normalised symmetrisation.
With the seed action specified and using the parameterisation (2), we expand the scale-dependent effective action in powers of the fluctuation fields, via the vertex expansion [24,29] where we have introduced the superfield Φ as a collection of all dynamical fields in our system, The Einstein summation convention over the superindex A implies the summation over discrete indices as well as an integration over the coordinates. Furthermore, Γ (n) k refers to the n-th functional derivative with respect to the superfield Φ. For convenience, we rescale the vertices according to As mentioned at the beginning of the section, we have introduced individual dimensionless couplings g n and λ n that label the different n-point functions. The vertices introduced in (12) also satisfy flow equations similar to (1), see e.g. [39]. Due to the breaking of diffeomorphism invariance by the gauge fixing term and the regulator, the scale dependence of these different couplings does not necessarily agree. In addition to the distinction of different dimensionless n-point couplings g n and λ n , also the couplings of pure gravity and ghost-gravity, or gravity-matter n-point vertices will differ, and should be distinguished. Therefore, for vertices originating from the ghost action (8), one generally has to replace the coupling g n in (12) by g c . For the course of this work, we will however work under the assumption that the differences between those couplings, caused by the breaking of diffeomorphism invariance, is negligible, and therefore assume that g n = g c = g, and λ n = λ, for n ≥ 3, and only distinguish the two-point functions from higher-order vertices. In this way, we assume the exact realisation of effective universality [39,[64][65][66], which refers to the semi-quantitative agreement of the scale dependence of different n-point vertices.
In the present work, we will focus on the two-point function of the pure gravity system. For the graviton two-point function, there are five independent tensor structures, three of them associated with the gauge fixing parameters α h , β h and γ h . The two remaining tensor structures, which will be unaffected by the choice of gauge, are related to the transverse-traceless mode h TT µν which satisfies and the scalar mode h 0 defined as where the scalar projector Π 0 is orthogonal to the gauge fixing action and the transverse-traceless projector, i.e., Explicitly, the transverse-traceless projector in momentum space is given by The projector on the scalar mode h 0 , referred to in (15) reads with The projectors Π TT and Π 0 can also be formulated in curved spaces, see e.g. [48]. We will now discuss how we resolve the two-point function. We start with the momentum-independent parts of the correlator. We can introduce two gaps µ TL and µ 0 which correspond to the traceless and the trace sector, respectively. These gaps µ x are related to the different n-point couplings introduced in (12) via µ x = −2λ 2,x , where the subscript x labels the respective mode. Explicitly, we introduce them by adjusting the two-point function via [42] where Π Tr and Π TL are the trace and traceless projectors, respectively: The introduction of the dimensionless quantities µ TL and µ 0 according to (21) ensures that, in the Landau limit, the two propagating modes feature a standard propagator with mass parameters µ TL and µ 0 , respectively. We emphasise that equation (21) is only a rewriting of the two-point function which conveniently allows to introduce the individual gaps µ TL and µ 0 .
To capture the full momentum dependence of the propagators, we introduce independent momentum-dependent wave function renormalisations for the transverse-traceless and the gauge-invariant scalar mode. Specifically, we rescale h according to where the wave function renormalisation tensor Z h , given by The rescaling (23) entails that the two-point function (21) needs to be multiplied with Z h from the left and the right: The scale dependence of the graviton wave function renormalisation is encoded in the anomalous dimensions We choose a spectrally adjusted regulator, i.e., a regulator which is proportional to the momentum-dependent part of the two-point functions, which ensures that mass-like terms are not regularised [19,35,125,164]. The regulator function R k implements the momentum-shell wise integration. Furthermore, we will choose the Landau gauge α h → 0, which is a fixed point for all gauge parameters [42,165]. With these choices, the graviton propagator reads where the coefficients B and C are defined in (19) and (20), respectively. Similarly, for the ghost, we will distinguish between the longitudinal and the transverse mode, and rescale the corresponding modes with momentum-dependent wave function renormalisations Z c L and Z c T , respectively. Specifically, the rescaling of the ghost will be The longitudinal and transverse projectors are defined in the usual way, The scale dependence of the ghost wave function renormalisations is encoded in the anomalous dimensions Therefore, after distinguishing all different modes on the level of the two-point function and choosing Landau gauge, there are four momentum-dependent anomalous dimensions and two gaps, completely parameterising the flow of the two-point correlation function.

General structure of the RG flows
Since we aim at investigating the flow of physical n-point functions, we employ the Landau limit, where the scale dependence of all gauge parameters vanishes. Focusing on the graviton sector in this limit, only the two physical projectors Π TT and Π 0 contribute to the propagator, see (28). When projecting the external legs of n-point functions on the physical modes, the full gauge dependence in the gravity sector is contained in the projector on the scalar mode Π 0 . For this projector, as defined in (17), we observe that This shows that, for the scale dependence of couplings induced by the graviton sector the gauge parameter γ h is redundant, since it can be removed by a rescaling of β h . While the cancellations in the ghost sector are more involved, and require non-trivial cancellations between contributions from propagators and vertices, the same rescaling for β h as given in (33) eliminates the γ h -dependence of the ghost-induced flows. Therefore, in the following we will employ γ h = 0 and keep β h general. Despite this redundancy of gauge parameters on the level of the scale dependence of couplings in the Landau limit, the tensor structure corresponding to γ h , (6), is nevertheless an independent tensor structure.
With the structure of the graviton propagator as given in (28), together with the regulator defined in (27), we can compute the product of the regulator insertion and two propagators, which enters every diagram of the flow as one of its building blocks, see Figure 1 and Figure 2. For the gravitational sector and in the Landau limit, it reads A similar expression can be computed for the ghost sector, where the product can be spanned by a transverse and a longitudinal part. Importantly, in the Landau limit, both the graviton propagator and the above product decay into a sum of the two projectors Π TT and Π 0 . This feature is due to a cancellation of gauge modes contained in the regulator with contributions from the propagators. It entails that the flow of n-point functions projected on gauge-invariant modes is only driven by the gauge-invariant modes themselves, and no gauge modes drive their scale dependence. It also significantly decreases the number of different tensor structures that need to be computed, when generalising the present computation to higher n-point correlators. Due to the decomposition into the orthogonal projectors Π TT and Π 0 , vertices with one or more gauge mode will not contribute to the scale dependence of any physical n-point correlator [107].
In the graviton sector, we project onto the different tensor structures of the two-point function by contracting the indices of (21) with Π TT and Π 0 , respectively. Structurally, after projection, the flow of the TT-part of the graviton two-point function reads [29,52]: where the last term on the left-hand side comes from the scale derivative acting on the k 2 coming with the dimensionless gap µ TL , see (21). Furthermore, the right-hand side is obtained by evaluating and projecting the diagrams in Figure 1. We also introduced the shorthand y = p 2 k 2 . In complete analogy, the flow equation for the scalar mode of the graviton two-point function reads Each of the equations (35) and (36) is disentangled by evaluating the right-hand-side at y = −µ TL and y = −µ 0 , respectively [29]: and In the ghost sector, the equations determining the anomalous dimensions are where the diagrammatic representation of flow c (y) is shown in Figure 2. The transverse and longitudinal parts of the flow are extracted by contraction with the corresponding projectors (31).
To evaluate the diagrams shown in Figure 1 and Figure 2, we used the Mathematica packages xAct [166][167][168][169][170]. The complete derivation is included in the attached notebook. We validated the correctness of the results with an independent code based on xAct as well as FormTracer [171].

On the relation between form factors and anomalous dimensions
We will now briefly discuss the relation between an action containing form factors and the momentum-dependent wave function renormalisations. On a flat background, the complete information on the graviton propagator is included in the action where f R and f C are form factors. The normalisation of the form factors ensures a unit prefactor in the propagator, see (41) below. This action, which is based on a curvature expansion of the full effective action, is an expansion in terms of diffeomorphism invariant operators. The form factors f R and f C capture the full momentum dependence of the graviton propagator. In comparison to a derivative expansion, they do not necessarily introduce new poles into the propagator. In a derivative expansion, one would in general expect the emergence of additional, potentially spurious poles of the propagators, see [120], and subsection 5.5. However, computations based on curvature expansions usually make use of the background-field approximation, where the difference between the background and the fluctuation propagators are neglected. This approximation potentially suffers from severe background dependence, which is lifted in the fluctuation expansion of the effective action (10).
Since gauge fixing and regulator both break diffeomorphism invariance, the expansion in terms of form factors as in (40), and in terms of fluctuation vertex correlation functions do not agree. Therefore, the question arises whether it is possible to extract the diffeomorphism invariant part -the form factors -from the scale dependences of the fluctuation two-point functions, without explicitly computing and solving modified Ward identities.
In the following, we will present a mapping between the form factors f R and f C , and the fluctuation two-point functions. For this, we will compare the propagators in both expansions, and neglect that these propagators are different due to the breaking of diffeomorphism invariance. A motivation for this assumption is the feature of effective universality [64], which suggests that, at least in the vicinity of the fixed point solution, the breaking of diffeomorphism invariance is mild. Therefore, under this assumption, we investigate the possibility to extract the form factors f R and f C from the fluctuation propagators.
For this, we will briefly indicate the form of the (unregularised) flat background propagator that arises from the action (40). We will also set the cosmological constant to zero for the moment. In this case, the scalar parts of the spin two and zero background propagators read This can be compared with the unregularised fluctuation propagator (28) at vanishing gaps, and suggests the following relation between form factors and wave function renormalisation: For finite values of the cosmological constant and the gaps, the situation is more complicated. The cosmological constant enters the two parts of the background propagators in a specific (gauge-dependent) ratio. By contrast, in general the fluctuation gaps need not to have any particular ratio. It is thus in general not possible to map the propagators in a one-to-one way without non-local terms in at least one of the form factors. If we nevertheless want to insist on an equivalence and allow for inverse powers of the momentum in at least one of the form factors, we can still map the propagators one-to-one. For example, we could identify the gap in the spin two sector with the actual cosmological constant, and introduce a term ∼ 1/p 4 in the form factor f R . Interestingly, such a term has been observed in the context of reconstructing an effective action from numerical lattice data obtained within causal dynamical triangulations [172]. Such a term also has interesting cosmological applications [173,174]. We will however not further discuss this issue here, since in the end, a proper evaluation of the problem necessarily involves the solution of the modified Ward identities.

Momentum-dependent anomalous dimension versus wave function renormalisation
As a final point in this section, we will discuss the relation of the fixed point condition for the dimensionless wave function renormalisation and the momentum-dependent anomalous dimensions. In the literature on momentum dependence in asymptotic safety, one commonly uses momentum-dependent anomalous dimensions rather than the wave function renormalisation. This is because only the anomalous dimensions are relevant in the flow equations -by construction all factors of the wave function renormalisation drop out, see e.g. (35). Since the anomalous dimension is related to the scale derivative of the dimensionful wave function renormalisation, one does not impose the fixed point condition directly. We can however easily translate between this language and a formulation where we require the dimensionless wave function renormalisation to be constant at the fixed point with respect to the RG scale k, see also [66].
The definition of the anomalous dimension η in terms of the dimensionful wave function renormalisation Z reads see (26) and (32). We want to express the right-hand side in terms of the dimensionless wave function renormalisation z for which we require a fixed point, and which is given by 1 In that process, we have to account both for the fact that we can normalise the wave function renormalisation, giving rise to an anomalous running induced by the anomalous dimension evaluated at zero momentum, as well as the scaling of the argument. In that way, we find 1 We refer to z(y) as the dimensionless wave function renormalisation, since its canonical mass dimension and the anomalous dimension cancel at the fixed point.
Here, the overdot indicates the logarithmic scale derivative with respect to the intrinsic k-dependence, excluding the trivial scaling of the argument. At a fixed point, whereż * = 0, we have We can solve this differential equation for the wave function renormalisation: Assuming that the momentum-dependent anomalous dimension is bounded, this has the consequence that for large y, the wave function renormalisation scales as Notably, only if the anomalous dimension vanishes asymptotically, that is, if it fulfils momentum locality 2 [31], the standard fall-off behaviour of the propagator in terms of the anomalous dimension at zero follows, see (28), namely If momentum locality is not fulfilled, we need non-local information on the momentum dependence, and the formula reads Having direct access to the fixed point wave function renormalisation via (47), we can calculate the fixed point form factor from the corresponding momentum-dependent anomalous dimension by inverting (42), Consequently, the behaviour at large momentum is where c is given by For a derivation of this, see appendix A. Let us finally make the connection to the derivative expansion, see also [29]. To quadratic order in y, we find z * (y) z * (0) This translates into the form factor Let us also emphasise at this point that in theories with more than one mode, and correspondingly more than one anomalous dimension, the derivative expansion is potentially unstable, the more so the more modes one has [175]. This comes from the fact that for a well-defined flow, one needs that for all modes z * (y) > 0. In particular, this implies that all the highest retained coefficient of all form factors need to have positive sign to avoid turning one of the modes into a ghost. For general theories, there is however no reason why all these coefficients need to be positive. One thus likely faces the dilemma that for a given order in the derivative expansion, the highest order coefficient of some mode might be negative. With a large number of modes, it is thus more likely that at no finite order, a fixed point, which is physically viable and does not necessarily feature ghost modes, can be found in a derivative expansion. This makes it clear that in complicated systems, the full resolution of the momentum dependence is not optional. In fact, as we will see in section 5, we find indications that quantum gravity with its two off-shell modes already shows such an alternating pattern for low orders of the derivative expansion.

Momentum-dependent fluctuation RG flow: analytical structure
Now, we will discuss several analytical properties of the momentum-dependent fluctuation two-point correlation functions which were introduced in section 3. For a clearer notation, we will work with dimensionless momenta, i.e. we make the identification p 2 k 2 → p 2 , such that no explicit factors of k will appear in the following discussion.
As has been found earlier, the Landau limit α h → 0 induces a fixed point for all gauge parameters [42]. Due to this fact, we only consider this limit for the fluctuation flows in this work. Additionally, there is the constraint which ensures that the ghost kinetic operator has the correct sign and is invertible. The most popular gauge choice in the literature is which we will focus on in our analysis. Other choices that have been argued for to be preferred are the transverse gauge β h = −1 [154,176] and the singular choice β h → −∞ which has been considered in e.g. [35,42,138,139]. The choice β h = 0 in the Landau limit allows for a particularly simple regularisation in a curved spacetime [48].
The analytical properties which we will discuss in the following will provide important cross-checks to test the numerical evaluation discussed in section 5.

Behaviour at small momentum
Let us first discuss the behaviour of the flow of the two-point function at small momentum. A key observation is that the flow contains terms of the form where p is the external momentum, q is the loop momentum, and x is the cosine of the angle between the two momenta. These terms look like unregulated, gapless propagators. Their origin lies in the projectors occurring in the propagators, see (16) and (17), and are thus simply a feature of gravity. The same phenomenon happens for spin one flows, where the transverse and longitudinal projectors also introduce such terms. These terms cause a number of technical difficulties, see also the discussion below in section 4.4. For small external momenta p, they are the reason that the derivative expansion does not commute with performing the loop integral. This is because expanding these terms in small p is actually an expansion in powers of p 2 /q 2 , so that higher orders of the external momentum introduce higher negative powers of the loop momentum. At a critical order, the negative powers cancel the measure term q d−1 and potential extra powers from the vertex factors. For higher orders the integrals thus do not converge as they suffer from IR divergences. This problem has been encountered before [177], where momentum-independent contributions of vertices have been neglected to define a finite, but ultimately inconsistent, flow. Even at a level where the integrals converge, strong instabilities have been found [29] which might prevent one from obtaining conclusive and robust results. Note that once the integral over the loop momentum has been performed, the flow appears to be smooth, see the numerical results presented below in section 5. This emphasises that reliable approximations must involve momentum dependence beyond a local expansion.
Besides this general discussion, one can investigate whether there are relations between the different flows at vanishing momentum. In the graviton sector, we find the relation This limit is formally singular, in part due to how we have to introduce the gap in the spin zero sector, see (21). In this limit, our ansatz formally diverges, indicating that we actually write down a gap for the spin zero gauge mode. This is consistent as this gauge choice projects out the trace mode.
In the ghost sector, we find that The divergence for β h → −∞ indicates that the longitudinal flow vanishes at p = 0 for this gauge choice. Due to the constraint (56), the flows, and thus the ghost anomalous dimensions, have the same sign at vanishing momentum.

Behaviour at large momentum
We will now discuss the large momentum behaviour of the flow of the two-point function. We emphasise that all of the following aspects a priori rely on our truncation of the three-and four-point function -a dynamical implementation of their flow can alter these results. In the following, we also assume that the regulator falls off exponentially. Similar conclusions hold for regulators with compact support, in which case no exponentially suppressed corrections appear. Some of the aspects that we present in the following have been previously discussed in [107,178].
As a general feature, we note that in the large momentum limit, the self-energy diagrams simplify. This is because of our assumption on the regulator, so that we can neglect it in the propagator which carries both the loop and the external momentum, up to exponentially small corrections. As a consequence, we can perform the angular integration exactly. The relevant integrals are which give rise to hypergeometric functions. We will now discuss the different parts of the flow in this limit in turn.

Spin two sector
It has been noted before that in d = 4 and with identified three-and four-graviton couplings, the spin two two-point function shows momentum locality [29,31,39]. This means that the flow of the two-point function goes to a constant at large momentum, in contrast to the naive expectation of a quadratic behaviour. This is due to a cancellation of the self-energy and the tadpole diagrams.
We found that this cancellation only happens in d = 4 -in higher dimensions, the two-point function rises quadratically at large momentum. We can also confirm that this qualitative behaviour is independent of the choice of gauge parameters and regulators. When we discuss numerical results below, this manifests itself by η TT going to zero for large arguments in d = 4, but staying finite in this limit in other dimensions, see Figure 3.
There are three contributions to the flow of the spin two two-point function: the graviton tadpole, the graviton self-energy, and the ghost self-energy diagram, see Figure 1. The tadpole diagram has the exact form flow tadpole TT The coefficients A tadpole TT,i depend functionally on the graviton anomalous dimensions and the regulators, as well as on the gaps, gauge parameters and the dimension. The structure of the self-energy diagrams is more involved due to the propagator depending on the sum of external and loop momentum. In the large external momentum limit, neglecting the regulator depending on the sum of the two momenta, we find that for the graviton self-energy diagram, and for the ghost self-energy diagram. Performing the integrals, we find the behaviour The I i TT are gauge-and dimension-dependent functions. This parameterisation exemplifies that only for g 3 = g 4 , the flow is momentum-local in d = 4. The ghost diagram does not contribute to the leading order behaviour.
More generally, one can ask if there are other (integer) dimensions where there is a choice g 4 = cg 3 for which the flow exhibits momentum locality. This is indeed the case for This is the only other combination of coupling identifications and integer dimension which gives rise to momentum locality.

Spin zero sector
The general structure of the spin zero sector is the same as that of the spin two sector. In particular, asymptotically we have an expansion as in (62), (63) and (64) for the different diagrams. The main difference is that we do not find momentum locality in d = 4. The only combination of coupling identification and integer dimension which gives rise to a momentum-local spin zero sector is From this we conclude that, at least in our setup, there is no situation where both spin two and spin zero sector are momentum-local. However, if also higher n-point functions are distinguished according to their tensor structures, this situation might change.
One can further ask the question, for d = 4, whether there is a specific ratio between the behaviour of the flows of the two sectors at large momentum. For this analysis, we assume again that g 4 = g 3 . In general, there is no such relation: the contributions including the spin zero propagator to each of the flows is generally different from the contributions including the spin two propagator. There are three exceptions to this. Two of them are independent of the choice of the gauge parameter β h : The third exception is the gauge choice β h → −∞, so that

Ghost sector
The flow in the ghost sector differs structurally from the flow in the graviton sector in that the tadpole diagram vanishes for a linear parameterisation of the metric fluctuations [18], see also (8), since in the present approximation the ghost action is linear in h. Neither of the two ghost modes shows momentum locality in our setup for any choice of dimension and gauge. However, independent of dimension and gauge choice, the two flows agree asymptotically, In this limit, only one of the two self-energy diagrams in Figure 2 contributes, namely the first one, where the regulator is inserted on the graviton line. Together with the relation of the two flows at vanishing momentum, this gives an estimate of their typical disagreement.

Limit of large dimension
Let us discuss aspects of the limit d → ∞. In this limit, we can perform the remaining angular integral exactly, and also comment on the radial integral. To see the first statement, note that this integral structurally looks like where f (x) is one of the integrands in the flow. We will now make the assumption that the dependence on x is smooth, so that f admits a convergent expansion in Chebyshev polynomials, Inserting this expansion into (71), and assuming that we can swap the order of summation and integration, we arrive at In the last line we expanded the expression to leading order in the limit of large d. This result is intuitively clear: for large d the measure suppresses all values of x except x = 0.
To deal with the radial integral, we will assume that we have a regulator which falls off like an exponential, so that also all integrands in the flow fall off exponentially. More precisely, with z denoting the square of the loop momentum, we assume that any integrand is of the form We included an additional power law behaviour so that we can assume that g is bounded on the whole interval. As a consequence, we can expand it in a series of rational Chebyshev functions, Combining this with the measure, we have integrals of the form If we combine these two results with the factor from the angular integration, we get the overall scaling of the flow as This result shows that there are two outcomes for the behaviour at large d, depending on the fall-off of the regulator. For the dimensional factor goes to zero, whereas for the flow diverges in this limit.

The flow for positive µ
As a last point in this section, we will discuss a technical challenge related to the projection procedure for the anomalous dimensions, see (35) and (36). For positive gaps µ > 0, we have to evaluate the flow at negative squared momentum. This is in general challenging, for two reasons. First, we have to define a regulator for, in general, complex momenta. Second, the factors (58) induce poles in the integration domain for p 2 < 0.
We choose a regulator in the following way. First, to avoid having to define it in the entire complex plane, we take the real part of the argument. Second, we have to define the desired behaviour of the regularised propagator for p 2 = −µ, that is we have to define In the limit of large masses, we require that this expression behaves like a standard regularised propagator, Third, we clearly also need that the regularised propagator does not introduce any poles in the integration region. A choice which fulfills all these requirements is For (large) positive arguments, this regulator is just of exponential type since the second term vanishes in this limit.
Let us now discuss how we deal with the unregularised pole structures (58) coming from the projectors. The poles are situated at x = 0 and q 2 = −p 2 = µ. We define the integration over these poles by splitting the integration domain into a disk centred around the pole, and the rest. Inside of the disk, we choose  radial coordinates. One can show analytically that then the Jacobian of the coordinate transformation and the angular integration remove all potential poles, so that the integral is well-defined. Let us finally mention that only simple poles arise in approximations which do not resolve the difference between the two graviton anomalous dimensions and gaps.

Momentum-dependent fluctuation RG flow: numerical results
In this section we discuss numerical results on the momentum dependence of the propagators, and the influence of the dimension, the choice of gauge and the size and sign of gaps. We will also discuss the derivative expansion of the form factors. All results are obtained with the regulator (83). As a generic choice, we set the gravitational coupling to one, g = 1, and use the Landau limit α h → 0.

Numerical strategy
Before we present the actual numerical results, we give a short discussion of how we obtained them. In the previous section, we found that the anomalous dimensions are bounded both at zero and infinite argument. We will further assume that they are bounded on the entire positive real line. In that case, they can be expanded in a series of rational Chebyshev functions [179],  Such an expansion is equivalent to compactifying the domain and employing a standard expansion in Chebyshev polynomials. This expansion shows desirable convergence properties if the function that is represented by the series is smooth, see e.g. [180] for an in-detail discussion. In the context of functional renormalisation group flows, they have been systematically discussed in [181][182][183]. For applications in quantum gravity, see also [68,101,145]. In practice, we will truncate (84) at a finite order, insert the expansion into the integral equations, and evaluate the equations at a set of collocation points. The set of integral equations for the anomalous dimensions then reduces to an algebraic set of equations for the expansion coefficients, which can be solved by standard linear algebra methods. The integrals have been performed with Mathematica's numerical integration routine.
To judge whether the numerical precision of the solution is high enough, we verify that the analytic relations between the different anomalous dimensions found in the last section are satisfied.

Dimensional dependence
First we will study the dimensional dependence of the anomalous dimensions. For the gauge choice β h = d 2 − 1 and vanishing gaps µ TL = µ 0 = 0, this is shown in Figure 3. There are a few general features. First, the overall magnitude of the anomalous dimensions decreases with increasing dimension. This is in agreement with our analytical estimate in the previous section, see (78). Second, we observe that while η 0 stays positive at these coordinates in all dimensions, η TT shifts to negative values quickly. Third, while there is some quantitative difference between the two ghost anomalous dimensions in d = 4, they agree  more and more the larger the dimension is. Fourth, as we have seen analytically in subsection 4.2, η TT is momentum-local only in d = 4. Finally, we observe a general flattening for large d: the momentum dependence of all anomalous dimensions is essentially trivial so that they can be very well approximated by a constant. However, this conclusion only holds in the current approximation, where higher-order curvature operators were neglected in the seed action. Since some of these operators are expected to be relevant in higher dimensions, their presence might alter the momentum dependence of the anomalous dimensions.

Gap dependence
Let us now discuss the dependence of the anomalous dimensions on the gaps µ TL and µ 0 . For this we fix d = 4 and the gauge parameters β h = 1. The anomalous dimensions are shown in Figure 4 for the gaps µ x = ±1/4.
Generally, we see that the gaps mostly influence the overall magnitude and the behaviour at small momenta. Rather small variations in the gaps can shift the value of η 0 (0) rather drastically. By contrast, the overall shape of the other anomalous dimensions is rather unaffected. We take this as evidence that the qualitative momentum dependence in large parts of theory space looks approximately similar to the one at the specific points that we present in this work.

Gauge dependence
Now we will discuss the gauge dependence of the anomalous dimensions. Once again we choose vanishing gaps and d = 4. The anomalous dimensions for the choices β h ∈ {−1, 0, 1} are shown in Figure 5. We find that overall the dependence on β h is mild, and mostly at the quantitative level. This is a promising indication that even though off-shell quantities like β-functions and propagators are inherently gauge-dependent, this dependence is well-controlled. As analytically expected, we find that the two ghost anomalous dimensions agree at p = 0 for β h = 0. Also, their value at infinity is numerically close to their value at zero. As a consequence, for this gauge choice a single, constant ghost anomalous dimension is an accurate approximation. Notably, also η TT and η 0 are numerically close at p = 0 for this choice of gauge.

Form factors and the derivative expansion
Let us finally discuss the form factors and their derivative expansion obtained via (51). As parameters, we choose d = 4 and β h = 1, µ TL = µ 0 = 0. While this is in general not a fixed point, we want to illustrate the general form of the fluctuation form factors. Since the general shape of the anomalous dimensions seems to be stable under variation of the gaps, see subsection 5.3, we expect that the qualitative picture presented here is correct for at least some part of theory space.
The fluctuation form factors, which were reconstructed according to (51), are shown in Figure 6. We find that the R 2 form factor f fluc R is positive while the C 2 form factor f fluc C is negative. Due to the positivity of f fluc R , there is no additional pole in the propagator of the gauge-invariant scalar mode. Similarly, for the transverse-traceless propagator, there is no additional pole. Both go to a constant value for small momenta and fall off with a power law at large momenta. Interestingly, the C 2 form factor f fluc C agrees qualitatively with that obtained in a background approximation of conformally reduced gravity [101].
For small momenta, we find This entails that the worst-case scenario for the derivative expansion explained in subsection 3.4 and in [175] is actually realised for this point of theory space: the Taylor coefficients of the two form factors have alternating signs, so a local expansion would not be able to resolve this point accurately. This means that in any derivative expansion in gravity which is currently technically feasible, viable fixed points might be missed due to the expansion.

Momentum-dependent background RG flow
In this section we provide results from a computation within the background field approximation. The minimal way to obtain the background form factors is to calculate their forms induced by the Einstein-Hilbert action. For simplicity, we will restrict ourselves to four dimensions and a harmonic gauge fixing, α h = β h = 1. This choice simplifies the calculation so tremendously to justify violating our preference of using the Landau limit. The detailed calculation is shown in appendix B. In particular, the resulting flow equations are given in eqs. (A35) -(A38).
For the form factors, the fixed point equations are linear first order differential equations, so that we can immediately calculate the solution. The free integration constant can be identified with the value of the functions at infinite momentum. With our choice of regulator, at Λ = 0 we can avoid additional poles in the propagator for To illustrate the solutions, we plot both form factors for two different choices of integration constants in Figure 7 and Figure 8. The first satisfies the above bounds and avoids new poles, whereas the second introduces new poles. We set the cosmological constant to zero and the dimensionless Newton's constant to one. This is in general not a fixed point, but allows a more direct comparison to the results of the fluctuation calculation.

Comparing background and fluctuation results
We will now compare the results obtained from the background and the fluctuation computation at the level of the form factors in d = 4. As our parameters we choose g = 1, Λ = µ TL = µ 0 = 0. For the background computation we choose the harmonic gauge α h = β h = 1, whereas for the fluctuation calculation we use α h = 0, β h = 1. The form factors are shown in Figure 6, Figure 7 and Figure 8.
Since the background computation is essentially one-loop, the form factors show a logarithmic behaviour at small momenta. By contrast, the fluctuation form factors go to finite values at vanishing momentum. The universal one-loop logarithms can be expected to come out correctly once we study the flow in the limit k → 0, at least in some part of theory space where effective field theory around flat spacetime is valid, see also [119].
For large momenta, the background form factors are either finite or zero, depending on the choice of integration constant. This is again due to restricting to one-loop. The correct behaviour will be fixed dynamically once the backreaction of these form factors onto the flow is taken into account. The fluctuation form factors show a power law fall-off according to (52). This entails that Regarding additional poles in the propagator, we have to choose the integration constants for the background form factors according to (86) to avoid them. For the fluctuation propagators, no new poles appear at the point that we investigated. The absence of additional poles in the propagators is crucial for the theory to be unitary. At the investigated point, the fluctuation results suggest that indeed no new poles arise, indicating that asymptotically safe quantum gravity might be unitary [120,[184][185][186]. Within the background field approximation however, the presence of additional poles is not yet conclusive, since in the present setup, they can only be avoided by choosing the integration constant appropriately. Avoiding this choice requires a more consistent computation, which includes the backreaction of the form factors.
Finally, the general shape and the signs differ between the two ways of computation. For the background computation, f R is negative for small momenta, whereas f C is positive, in both cases due to the logarithm. For larger momenta, the signs can, but need not change, depending on the integration constant. The fluctuation form factors have a definite sign: positive for f fluc R and negative for f fluc C . They do not show any non-trivial feature besides what is dictated by the asymptotics: they show an approximately constant regime for small momenta, a cross-over at around p ≈ k, leading to the power law fall-off at large momenta.
It is difficult to predict which of these features are generic in most parts of theory space, and which are special to the point of investigation. We leave this investigation to future work, together with a detailed analysis of the momentum-dependent flow of interaction vertices.

Summary and Outlook
In this work we discussed the full non-perturbative momentum dependence of the graviton and ghost propagator in quantum gravity. This includes the two different graviton modes with spins zero and two, and the two ghost modes, with spins zero and one. We obtained the results with the help of functional renormalisation group equations, and resolved both gauge and dimensional dependence.
A key result is that the propagators of the different graviton modes behave qualitatively differently. In four dimensions and within our approximation, the spin two mode is momentum-local, whereas the spin zero mode does not show this property. This qualitative result is independent of the gauge choice, and thus potentially has physical significance [31,107]. By contrast, the two ghost modes agree qualitatively and partially even quantitatively, and their anomalous dimensions are approximately momentum-independent.
The dependence on the choice of gauge as well as on the gaps is largely quantitative rather than qualitative. We take this as a sign that the selected examples presented in this work are indeed representative of larger parts of theory space. Furthermore, the weak dependence of the anomalous dimensions on the choice of gauge is a promising indication that, even in truncations, the FRG gives rise to reliable and stable results.
In dimensions larger than four, we find that all anomalous dimensions quickly go to zero. Assuming a regulator which falls off quickly enough, this suppression for large dimensions is exponential in the dimension d.
All these results have been obtained from calculations of dynamical correlation functions. To also compare to a computation in a background field approximation, we have calculated the background form factors induced by an Einstein-Hilbert action to quadratic order in the curvature. To facilitate this comparison, we created a dictionary relating the form factors to anomalous dimensions, both evaluated at a fixed point. We find that the one-loop structure of our approximation in the background sector makes such a comparison difficult, since many features are dictated by the one-loop form rather than dynamically, e.g., the asymptotic behaviour for both large and small momenta. A dynamical computation of the background form factors with backreaction is desirable but postponed to the future.
The present work represents a major step in a complete and systematic computation of correlation functions in quantum gravity. The natural next step is to resolve the three-graviton vertex. Since the complete vertex is rather complicated, two possible paths are the resolution of the full momentum dependence of selected tensor structures, or the complete resolution of all tensor structures in a derivative expansion. While the former strategy allows to extract the full momentum dependence in one specific sector, it is not guaranteed that such an approximation has desirable features the full theory should have [175]. In the latter strategy, all tensor structures are taken into account, however, as we have seen in subsection 5.5, a derivative expansion might not be able to resolve the fixed points of the theory accurately.
Finally, we expect that disentangling the different graviton modes will play an important role in the discussion of spectral functions [119]. In particular, the two sectors will have individual spectral functions, and their different high-momentum behaviour will have a non-trivial effect.
Author Contributions: All authors contributed equally to this article. All authors have read and agreed to the published version of the manuscript.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. Asymptotic expansion of form factors
In this appendix we derive the asymptotic expansion (52). For this, we consider the integral in the exponent of (51), We will assume that η * is a bounded function on the whole positive real line, which is indeed true in our case. With that, we can first show that the leading order behaviour of (A1) at large y is logarithmic. For that, let us take the derivative of (A1): Since η * is bounded, for large y, we find that As a consequence, the integral itself is asymptotic to a logarithm, We still have to determine the sub-leading behaviour. For that, let us take the difference of the integral and the logarithm above, and take the limit of y going to infinity. We will find that this limit exists, so that the sub-leading part is a constant. To achieve that, let us rewrite the logarithm via ln y = y 0 ds 1 1 + s − ln 1 + y y .
The virtue of this rewriting is that the logarithm on the right-hand side vanishes in the limit of large y. We can then combine the integrals so that We can now take the limit y → ∞. The integral converges in this limit, since the terms falling off like 1/s at large s cancel, so that finally we get This completes the derivation of the asymptotic formula (52).

Appendix B. Calculation of the background form factors induced by GR
In this appendix we derive the flow equations for G N and Λ as well as the two-curvature form factors induced by the Einstein-Hilbert truncation in a background field approximation and in four dimensions. To make our life easier, we will choose the harmonic gauge condition Even though we argued in the main text that only the Landau limit allows for a clean flow, in a background field calculation this gauge choice simplifies the calculation tremendously. The reason for this is that in this gauge, the two-point function is diagonal and only depends on Laplacians and curvatures, while it does not involve uncontracted derivatives. Splitting the graviton into traceless and trace parts, the individual two-point functions at vanishing fluctuation fields read Here we introduced the operators These are also the operators that we regularise. This yields the very simple floẇ Here, we used the same shape function in all modes. For later convenience, we introduce the shorthands In this notation, the flow readṡ To improve the readability, in the following we will refrain from using an overbar for background quantities.

Early time heat kernel expansion
To evaluate the traces, we will employ the early time expansion of the heat kernel. For an operator of Laplace type ∆ + E, where E is an endomorphism with the appropriate bundle structure, the expansion reads [187] Tr e −s(∆+E) In this expression, Tr is a functional trace, tr is the trace over the bundle indices, F is the bundle curvature which depends on the index structure of the field traced over, and The Taylor expansion shows that all form factors are regular at zero. To bring the traces (A12) into the form of the standard heat kernel trace (A15), we use the inverse Laplace transform. For a function g of an operator O, this entails that we can write To be in line with the basis of the main text, we can use the relation to replace any occurrence of Ricci tensors by Ricci scalars and Weyl tensors. In the following we will also neglect the Euler characteristic, which corresponds to the case n = 0 in the above equation.

Some integral transformations in d=4
To bring the flow equations into a convenient form, we will use the formula ∞ 0 dsg(s)s −n = 1 Γ(n) ∞ 0 dz z n−1 g(z) , which holds for n > 0. Also, the following identities hold true: In all these equations,g is the inverse Laplace transform of an arbitrary function g, see (A23), and f is the universal heat kernel function (A22). The equations can be shown by inserting the integral form of f and performing repeated variable transformations while exchanging the order of integrals.

Trace contribution of the spin zero part
We start with the contribution of the trace part of the graviton to the flow. The trace projector makes it a scalar trace, with F µν = 0 , Using the inverse Laplace transform (A23), then applying the general trace formula (A15), and using the integral identities (A25)-(A28) for this case and in d = 4 gives (A30)

Trace contribution of the ghost
Next, we present the contribution of the ghost. In this case, With this, the trace is

Trace contribution of the spin two part
Finally, we present the contribution of the traceless spin two component of the graviton. With we find (A34)

Background flow equations
Now that we have computed all traces, we can read off the beta functions of G N , Λ and the two-curvature form factors by comparing the coefficients of the scale derivative acting on the action (40), with the results for (A14), given in (A30), (A32) and (A34). We will present them for the dimensionless quantities g, λ, f R and f c . They reaḋ Here, we introduced the combined measures µ R h (u) = 13 8 (1 − 4u) 3/2 . (A39)

Fixed point structure
Let us now analyse the fixed point structure of these background flow equations. Employing the fixed point conditionġ =λ =ḟ R =ḟ C = 0 , we find that the form factor equations are first order linear ordinary differential equations of the form 2z f (z) = for some kernel K. The explicit solution to this equation reads We chose to impose the boundary condition at ∆ = ∞ since for this choice the integrals on the right-hand side converge for all ∆ ≥ 0. Inserting the explicit expressions for the kernel, we can also perform the z-integration. With the shorthand notation the explicit induced background fixed point form factors are f C (∆) = f C (∞) + g 2π 1 0 dω ω ν 3, −6, Some of the integrals can be performed analytically: while the others have to be evaluated numerically.