Renormalization group in super-renormalizable quantum gravity

One of the main advantages of super-renormalizable higher derivative quantum gravity models is the possibility to derive exact beta functions, by making perturbative one-loop calculations. We perform such a calculation for the Newton constant by using the Barvinsky-Vilkovisky trace technology. The result is well-defined in a large class of models of gravity in the sense that the renormalization group beta functions do not depend on the gauge-fixing condition. Finally, we discuss the possibility to apply the results to a large class of nonlocal gravitational theories which are free of massive ghost-like states at the tree-level.


I. INTRODUCTION
The calculation of quantum corrections always had a very special role in quantum theories of gravity. The first relevant calculation was done by t'Hooft and Veltman [1], who derived the one-loop divergences in the quantum version of general relativity, including coupling to the minimal scalar field. Soon after similar calculations were performed for gravity-vector and gravity-fermion systems [2]. These first calculations have a great merit, regardless of the fact that the output was shown to be gauge-fixing dependent [3]. Later on one could learn a lot from the two-loop calculations in general relativity [4,5]. Technically more complicated are calculations in four-derivative gravity, which were first performed in [6] and with some corrections in [7], [8] and finally in [9], where some extra control of the calculations was introduced and the hypothesis of the non-zero effect of the topological Gauss-Bonnet term [10] explored. Let us also mention similar calculations in the conformal version of the four-derivative theory [7,11,12].
The importance of four-derivative quantum gravity is due to its renormalizability [13], which is related to the presence of massive unphysical ghosts, typical in the higher derivative field theories. Naturally, there were numerous and interesting works trying to solve the unitarity problem in this theory [14][15][16]. The mainstream approach is based on the expectation that the loop corrections may transform the real massive unphysical pole in the tensor sector of the theory into a pair of complex conjugate poles, which do not spoil unitarity within the Lee-Wick quantization scheme [17]. However, it was shown that the definite knowledge of whether this scheme works or not requires an exact non-perturbative beta-function for the coefficient of the Weyl-squared term and for the Newton constant [18]. The existing methods to obtain such a non-perturbative result give some hope [19], but unfortunately they are not completely reliable 1 and, therefore, the situation with unitarity in the four-derivative quantum gravity is not certain, at least.
Another interesting aspect of quantum corrections in models of gravity is related to the running of the cosmological constant Λ cc and especially Newton constant G. These quantum effects may be relevant in cosmology and astrophysics (see, e.g. [21,22]) and can be explored in different theoretical frameworks, such as semiclassical gravity [23], higher derivative quantum gravity [7,15], low-energy effective quantum gravity [24], induced gravity [25] and functional renormalization group [26]. Indeed, the status of the corresponding types of quantum corrections is different, but there are also some common points. In particular, in many cases one can formulate general restrictions on the running of the Newton constant G, which are based on covariance and dimensional arguments [27]. The beta-function for the inverse Newton constant which follows from these condition has the form µ d dµ where m i are masses of the fields or more general parameters in the action with the dimensions of masses and A ij are given by series in coupling constants of the theory. In the perturbative quantum gravity case there may be one more complication. In the model based on Einstein-Hilbert's gravity there is no beta-function for G, and in the four-derivative model this beta-function is dependent on the choice of gauge fixing condition [7,28,29]. Only a dimensionless combination of G and Λ cc has well-defined running, but this is not sufficient for the mentioned applications to cosmology and astrophysics.
Recently there was a significant progress in development of perturbative quantum gravity models which have very different properties. If the action of the theory includes local covariant terms that have six or more derivatives, this theory may be i) super-renormalizable [30]; ii) unitary, in case of only complex conjugate massive poles in the treelevel propagator [31,32], and iii) have gauge-fixing and field reparametrization independent beta-functions for both G and Λ cc . The theory is unitary at any order in the perturbative loop expansion when the CLOP [33] prescription is implemented or non-equivalently if the theory is defined through a non-analytic Wick rotation from Euclidean to Minkowskian signature [34][35][36]. This means that such a theory satisfies the minimal set of consistency conditions and deserves a detailed investigation at both classical and quantum levels. The classical aspects of the theory started to be explored recently in [37], where it was shown that the version with real simple poles has no singularity in the modified Newtonian potential. Quite recently this result was generalized for more general cases including multiple and complex poles [38,39]. Furthermore, in [38,40] the detailed analysis of light bending in six-derivative models was given. Another generalization of simple higher derivative model is nonlocal gravity, where we allow for nonlocal functions of differential operators [41]. Moreover, there exists a class of nonlocal theories in which UV behaviour is exactly the same like in polynomial higher derivative theories. Therefore, they also satisfy the above three points, namely they are quantum super-renormalizable models and the analysis of divergences and RG running presented here apply to these theories as well.
Until now, the unique example of quantum calculations in the super-renormalizable quantum gravity was the derivation of the beta-function for the cosmological constant in [30]. Here we start to explore the models further and derive the most relevant phenomenologically one-loop beta-function for the Newton constant G.
The work is organized as follows. In Sect. II one can find a brief general review of the super-renormalizable models [30], including power counting and gauge-fixing independence of the beta-functions. In Sect. III, we describe the oneloop calculations. Some of the relevant bulky formulas are separated into Appendix A, to provide a smooth reading of the main text. In Sect. IV, two important classes of the nonlocal models of quantum gravity are considered. It turns out that the derivation of one-loop divergences by taking a limit in the results for a polynomial models meets serious difficulties, which can be solved only for a special class of nonlocal theories, which are asymptotically polynomial in the ultraviolet regime. But still these theories are super-renormalizable or even finite. In Sect. V, the renormalization group for the Newton and cosmological constants are discussed, within the minimal subtraction scheme of renormalization. Finally, in Sect. VI we draw our conclusions and outline general possibilities for further work.

II. SIX-AND HIGHER-DERIVATIVE QUANTUM GRAVITY
One of the simplest model of super-renormalizable quantum gravity is based on the action Here R is the scalar curvature and C is the Weyl tensor (with writing of all indices suppressed), e.g.
In the last formula Riem and Ric stand for the Riemann and Ricci tensors (with writing of all indices suppressed again), correspondingly. Furthermore, the generalized Gauss-Bonnet term is which is not topological for n = 0. It can be easily shown that this term can be reduced to the O(Riem 3 ) terms for these values of n. The coefficient of the Einstein-Hilbert term and the density of the cosmological constant term are denoted as ω EH = −1/(16πG) and ω cc = −Λ cc /(8πG) , in order to provide homogeneous notation. Action (2) is the most general one which is at most quadratic in curvatures. One can extend it preserving power counting by adding extra terms that are of the higher orders in the curvature tensor [42]. As we shall discuss in what follows, such O(Riem 3 )-terms are not necessary for super-renormalizability and in this sense they represent the non-minimal sector of the theory. In what follows we will see that the Einstein-Hilbert and cosmological terms are relevant, but for N 2 they do not affect the divergences. It proves useful to discuss the dimensions of the parameters included into the action (2). Standard considerations show that ω N,R , ω N,C , ω N,GB = (mass) −2N and ω N −1,R , ω N −1,C , ω N −1,GB = (mass) −2N +2 . Therefore, by dimensional reasons similar to the ones that lead to (1), one can show that the beta-function for the inverse Newton constant should be proportional to the ratios In what follows this conclusion will be supported by power counting and by the direct calculation of the one-loop counterterms.
Let us consider the power counting in the theory (2) which is equivalent to the most general local theory with 2N + 4 derivatives of the metric. The formula for the superficial degree of divergence D for a p-loop diagram follows from a general expression and the topological relation Here d is the number of derivatives on the external lines of the diagrams, r l is the power of momenta in the inverse of the propagators, K v is the power of momenta in the given vertex v, n v = v is the total number of vertices and l int is the number of internal lines of the diagram. One can introduce the gauge-fixing condition in the theory (2) in such a way that effectively r l ≡ 2N + 4 (for this one can check [30] and the next section of the present work for the details). The strongest divergences come from the largest K v and hence one can set K v = r l = 2N + 4, without losing generality. Then it is easy to arrive at the result For the logarithmic divergences (D = 0) we get the estimate for the dimension of the p-loop counterterms, d = 4 + 2N (1 − p). The last expressions have several important consequences. First of all, from (7) one learns that the theories with N 1 are super-renormalizable, since only diagrams with p = 1, 2, 3 may be divergent. Moreover, for N 3 only one-loop diagrams may be divergent. As a result, the one-loop beta-functions will be exact in the the theories with N 3.
Second, since the divergences come from p 1, it is clear that for any N 1 the dimension of the counterterms can only be d = 4, 2, 0. As a result, the terms which are included in the first two lines of the action (2) are not subjected to the renormalization procedure for N 2. The multiplicative renormalizability requires that the action (2) includes the terms in the last line, which we can also parametrize as where is the Einstein-Hilbert action with the cosmological constant and is the action with four-derivative terms. We included the surface term ✷R, along with the Euler density term E. Third, for N 2 the renormalizations of the parameters a 1,2,3,4 and G, Λ cc of the actions (9) and (10) do not depend on these parameters of the action, but only on the parameters of the higher derivative terms of the action (2) (first two lines of it). Let us note that the renormalization of Λ cc /(8πG) was described in details in [30] and therefore we will not discuss its details here. As far as we intend to calculate the beta-function for the Newton constant G, it is sufficient to consider only the quantum effects of the theory with the action (2). Indeed, the result can be modified by adding O(Riem 3 )-type terms [30,32,[42][43][44] (sometimes called "killers" [42] for the ability to make the theory free of any divergence), but we will not consider these terms here. One can note that the remaining terms, which do not affect the divergences, can provide the complex massive poles and hence lead to the unitarity in the Lee-Wick sense. Therefore, the theory (2) turns out to be quite general and deserves an explicit calculation of the physically relevant beta-function for G.
Fourth, an important observation is that for N 1 the classical equations of motion in the theory have 4 + 2N 6 derivatives, while the divergences are removed by the counterterms which have 0, 2 and 4 derivatives. It is well-known that the gauge-fixing dependence of the one-loop counterterms is always proportional to the classical equations of motion. The general proof of this statement can be found in [45], and one can find a specific proof for the one-loop contributions in the book [46]. The general theorem enables one to establish gauge dependence in second-and fourthorder quantum gravity without explicit calculations [7,28]. And in the case of six-or higher than six-derivative models the same statements guarantee the gauge-fixing independence of all divergences.
Finally, we conclude that it is completely consistent to calculate only the beta-function for G, without deriving the beta-functions for the coefficients of the four-derivative terms a 1,2,3,4 . Moreover, the calculation of the beta-function for G can be performed in any useful gauge, because the result does not depend on the choice of gauge-fixing.

III. ONE-LOOP CALCULATION
The calculation of one-loop divergences can be performed within the background field method. The splitting into background and quantum metric is as follows: The Lagrangian quantization is based on the Faddev-Popov scheme requires introducing the gauge-fixing condition χ µ and the weight operator C µν . The one-loop effective action is given by the expression [7] whereĤ is the bilinear form of the action (2) with the gauge-fixing term added with respect to gravitational perturbations h µν , andĤ gh is the bilinear form of the action of the gauge ghosts.Ĉ is the bilinear form of the gauge-fixing action with respect to gauge-fixing conditions, which is related to the weight operator. All operators in (12) are second variational derivatives of the corresponding actions and therefore are Hermitian. Moreover, they should be divided by the corresponding powers of the renormalization parameter µ in the dimensional regularization, which are not shown in the last formula. Also, in order to have homogeneous propagators for the quantum metric and the ghosts, one can, e.g. replaceĤ gh = M α β by the productC γα M α β , while the coefficient of the weight operator contribution in (12) changes from −i/2 to −3i/2.
In this paper the one-loop computations have been done in Euclidean signature. For Lee-Wick theories the continuation to Minkowski space is guaranteed consistently with perturbative unitarity by one of the two alternative procedures described below.
i) The CLOP prescription [33] gives an unambiguous (surely at one loop), Lorentz-invariant and unitary result. It consists on taking the masses of the complex conjugate poles to be unrelated complex mass parameters to avoid the overlap of the poles for finite value of the external energy. At the end of the computation one has to impose the condition that the poles are complex conjugate to each other. In particular, the CLOP prescription allows to make the Wick rotation to Euclidean signature, as explained in [36].
ii) Another possibility is to make a non-analytic continuation from the Euclidean to the Minkowski space as recently proposed in [34,35]. In both described approaches the S−matrix turns out to be non-analytic.
Let us note that the described complications concern only the case when the massive poles are complex conjugate, while for the local theories with real ghosts there is no need to worry about the Lee-Wick quantization or about continuation to Minkowski space, since unitarity is already violated at tree-level. Indeed, the one-loop calculation which will be described above does not distinguish the theories with real or complex ghosts.

A. Gauge-fixing and relevant operators
The general form of the gauge-fixing action is where the gauge condition is Here and below the covariant derivatives are constructed with the background metric. It is sufficient for us to take the weight function C µν in the following form ByC µν in (12) we mean the self-adjoint extension of the operator C µν . The action of the gauge ghostsc µ and c ν has a standard form the operator M is defined as usual (see, e.g., [46], because in this part there is no essential difference with the four-derivative quantum gravity), such that we arrive at At that point one can observe the first essential simplification coming from our interest in the Einstein-Hilbert counterterm only. For the gauge-fixing action (13) with (14) and (15), both the weight operatorC µν and the bilinear operator of the ghost action M µ ν in (17) are homogeneous functions of covariant derivative ∇ α and do not include dimensional parameters. As a result, the corresponding functional determinants contribute only to the divergences of the four-derivative terms in (10). Since we are interested in the beta-function for the Newton constant only, we can concentrate our attention only on the first term in the expression (12) and do not pay attention on the other two terms.
Let us consider the relevant first term in (12). We remember, once again, that the final result for divergences does not depend on the choice of χ µ and C µν . Therefore, these functions can be defined to make calculations as simple as possible. The primary role of the gauge-fixing is to make the operatorĤ non-degenerate and this can be achieved for different values of the gauge-fixing parameters α, β, and γ in Eqs. (14) and (15). Due to the gauge-fixing independence of the one-loop effective action one can choose the gauge-fixing parameters in such a way that operatorĤ assumes the minimal form. This means that the highest-order derivatives form the combination N +2 . The corresponding calculations are essentially the same as in the four-derivative quantum gravity. The requirement of minimality leads to the following values: One can see that the gauge-fixing parameters β and γ depend only on the ratio ω N,R /ω N,C and are independent on other coefficients appearing in the gravitational action (2). In particular, the minimal choice of the gauge-fixing conditions does not depend on the generalized Gauss-Bonnet term, namely ω N,GB . This is obviously a natural output because the minimal gauge-fixing can be defined on a flat background, where the O(Riem 3 ) terms are irrelevant. Furthermore, the parameter α is only dependent on ω N,C . This is because only the term C N C in the action is responsible for the propagation of the spin-2 part of gravitational perturbations on flat space-time with the highest number of derivatives. The evaluation of the main expression can be performed by using the generalized Schwinger-DeWitt technique developed by Barvinsky and Vilkovisky in [47], in a way qualitatively similar to the calculation of the cosmological constant counterterm in [30]. Technically, the calculation of the linear in R term is much more complicated, and we shall describe it in certain details. The minimal choice of the gauge parameters (18) provides the useful form of the operator [30] where B µν, κλ is the DeWitt metric in the space of the fields for our model, where δ µν, κλ = 1 2 g µκ g νλ + g µλ g νκ is the identity matrix in the space of the fields. In what follows we assume that the c-number operator B µν, κλ is non-degenerate, such that its inverse has finite coefficients. Let us note that the degenerate case, when ω N,C = 0 or ω N,R = 0 , is non-renormalizable, and therefore is much less interesting. Acting on (20) with the inverse operator (22), one can see that the operator H ′ κλ ρσ has the standard minimal form with the senior term of exactly 2N + 4 derivatives, and with the identity matrix in the space of the fields δ ρσ κλ = 1/2(δ ρ κ δ σ λ + δ ρ λ δ σ κ ) in front instead of the DeWitt metric B µν, κλ as in (20): Obviously, the first factor in Eq. (20) does not contribute to the divergent part of the effective action, and therefore for us Different from the calculation of the cosmological constant counterterm [30], the expressions for U and V requested for the derivation of divergences for an arbitrary metric background are very bulky and difficult to derive and to deal with. However, there are possibilities of great simplifications when we are interested only in the counterterm linear in scalar curvature. First of all, by dimensional reasons the mass dimensions of V, W and U are 2, 3 and 4, correspondingly. If only coefficients ω N,i (i = R, C, GB) entered the expression for H ′ tensor, then only U could be relevant in our case. However, the other coefficients ω N −1,i enter too, hence both V and U parts are relevant here. Furthermore, since we do not pretend to calculate the R-term, in the course of calculations we can assume that the scalar curvature R is a constant, without losing the generality. Finally, it is good to remember that the Ricci and Riemann tensors enter the final expressions for the linear in scalar curvature term only through their contractions with the metric tensor. Therefore, again without loss of generality we can assume that the curvature tensor is where Λ is a constant. This is equivalent of assuming that the background is a maximally symmetric space-time (dS or AdS) with the radius of curvature given by ± 3/Λ. Due to the covariant constancy of the purely background metric we have W = 0. At the end one can set Λ = R/4 and arrive at the general result of the counterterm in the one-loop effective action, which is valid for both constant and non-constant scalar curvature R, in the given approximation.
On the basis of dimensional analysis, in the tensor V there are terms of zero and first powers in the background curvature Λ, that can be written schematically as while for the tensor U there are terms of the following types that means up to the second power in background curvature. The dimensions of the different types of terms are compensated by the coefficients, which may have different powers of dimensional parameters of the action (2), e.g. by the ratios ω N −1,i /ω N,j , where i, j = R, C, or GB. The details of this structure can be found in the Appendix A.
In order to use the generalized Schwinger-DeWitt technique of [47], we can first extract the highest derivatives in the form of N +2 [30], Tr ln H ′ κλ ρσ = Tr ln N +2 + Tr ln δ ρσ κλ + V κλ ρσ, α1...α2N+2 ∇ α1 · · · ∇ α2N+2 1 N +2 and after that perform a series expansion of the logarithm. The first term in the last expression does not produce linear in curvature divergences and the last term includes terms that can not give relevant contributions by dimensional reasons. In what follows we omit the last term. The expression can be further elaborated as follows, reducing to the universal traces of [47], Let us start the analysis of the last expression with two simple observations. The commutators of covariant derivatives with V and to negative power operators in the last term were disregarded because they are irrelevant in the linear in curvature approximation, as it was explained above. According to the Eq. (29) and the table of universal traces of [47], the trace of the operator proportional to the tensor W vanishes identically, which confirms our previous consideration. Therefore, we need to consider only linear in V and U single traces and the quadratic in V mixed (interference) trace. Due to the fact that in the tensors V and U we have at most four free covariant derivatives (not contracted with other derivatives) we can simplify the relevant parts of the expression (29) to the form + (terms not contributing to the divergence linear in R) .
The explicit expressions for the c-number tensors U and V can be found in Appendix A.
We here make some comments on the correct way of taking the traces that are in general not an easy task. It is important to note that the expression above includes all the divergences that are linear in R. A special care must be given to the interference term − 1 2 TrV 2 . Due to dimensional reasons we are interested in the mixing terms only between the parts of the V tensor that are linear in Λ or independent on Λ. Only these terms can give the correct dimension of the coefficient of R as well as the correct power of background curvature. Then for the mixing terms we do not need to account for additional background curvature arising from commutations. Furthermore, when taking the single trace of V we need only terms independent on Λ. Contrary, in the traces that involve the tensor U we are interested only in those parts that are proportional to Λ.
The final technical comments are the following. The interesting parts, which give divergent contributions linear in curvature, of V and U are linearly proportional to the couplings ω N −1,C , ω N −1,R or ω N −1,GB and originate from the variation of the corresponding part of the action (2). These contributions must be independent on Λ in the V part and linear in Λ in the U part, while the denominators with ω N,R or ω N,C emerge only because of the inverse DeWitt metric (22). To derive the relevant U term correctly one has to take into account the commutations up to the order linear in background curvature in the second variation of the gravitational action (2).
The situation with the interference term is different. The first interference factor in V has the same origin as in the single trace discussed above. However, the second one must be precisely linear in Λ and independent on the couplings ω N −1,i . It contains some rational functions of the dimensionless ratios between the coefficients ω N,R , ω N,C and ω N,GB . The origin of this term is either from the second variation of the gravitational action or from the gauge-fixing term, both multiplied by the inverse DeWitt metric. To derive this term correctly one has to take into account the commutations up to the order linear in background curvature both in the second variation of the gravitational action or the gauge-fixing term. In the intermediate expressions one can observe rather complicated rational functions, and only in the final traces there are essential simplifications. The precise expressions for relevant parts of tensors V and U are not given here due to their lengths.
Actually in (30) we will never need to take universal trace with more than four free covariant derivatives because of the structure of the indices and because only the metric functions appear in the V and U tensors. Therefore, the last but one term in (30) will never be so long after explicit contractions.

B. Divergences for N 2 models
After taking the traces according to [47], we can write the results for the divergent part of the one-loop effective action. The part linear in curvature is given by the expression where ε = (4 − n)/2 is the parameter of dimensional regularization, µ is the dimensionful regularization parameter and by n we denote space-time dimensionality. One can see that the generalized Gauss-Bonnet terms with the coefficients ω N,GB and ω N −1,GB give non-zero contributions in the formulas above, differently from the four-derivative quantum gravity case [9,12]. Let us stress that the result for the divergence (31) is valid only for N 2 and there is no smooth limit N → 0 or even N → 1. The reason for this discontinuity is that the expansion of into the quantum field h µν generates expressions of the −1 and −2 type, which are badly defined when acting on a constant. As a result it becomes impossible to take the limits N → 0 and N → 1 in the bilinear form of the generalized Gauss-Bonnet term. After this expansion the formulas do not "recognize" that for N = 0 this term is topological. This is why the coefficients of the ω N −1,GB and ω N,GB terms are proportional to 2N − 1 and 2N − 3, and not to N , as one might have expected from them. Indeed, the same discontinuity takes place also for the cosmological constant -type divergence derived in [30]. We have checked independently that for N = 0 the Gauss-Bonnet term does not contribute to the divergences. The details are not shown here, because one can find them in the previous works [12] and [9].
If the relevant generalized Gauss-Bonnet terms are absent, ω N −1,GB = 0 and ω N,GB = 0, we arrive at a very compact form of divergences, which is valid for N 2, One can note that in the last expression the number of derivatives N does not appear explicitly, but only through the coefficients of the action (2). This is in contrast to to the more general formula (31), where the explicit dependence on N takes place in the sector related to the generalized Gauss-Bonnet terms. For the sake of completeness, let us write down the expression for the cosmological constant divergences, derived in [30], in the notations which we use here, The contribution of the generalized GB term to the cosmological constant divergence is zero because this term does not contribute to the propagator around flat spacetime. It is useful to remember that both expressions (32) and (33) are valid only for N 2. It is easy to see that the relative factor of "5" between the contributions of tensor and scalar modes in Eq. (32) is exactly the same which one could observe in the cosmological constant divergence (33).
It is clear from (32) that one can provide zero divergence in the Einstein-Hilbert sector by adjusting the highestderivative coefficients of the classical action ω N −1,C , ω N,C , ω N −1,R and ω N,R . This can be achieved even without using the freedom to choose the coefficients of the ω N −1,GB and ω N,GB in (31). Moreover, one can make the same for both the cosmological and the linear in R divergences at the same time.
One can make zero the two remaining relevant coefficients for the C 2 and R 2 counterterms by means of the "killer" operators. Thus one can achieve the one-loop finiteness, which can be directly extended into all-loop order for the case of N 3, since in this case divergences exist only at the one-loop level.
C. The special cases of N = 1 and N = 0 In order to achieve a better understanding of the situation with smaller N , let us present the corresponding results, without going into full details. As we have explained above, the cases of N = 0 and N = 1 should be considered separately. The calculations are pretty much the same as for N 2, therefore, we present only the final results.
For N = 1 the action includes the following relevant terms: The term linear in curvature in the divergent part of the one-loop effective action is: One can easily note that this expression differs from the N 2 case. In particular, there is a mixing of the coefficients for the Weyl and scalar curvature. Moreover, we notice that only the terms with the Weyl curvatures and with the generalized Gauss-Bonnet term give rise to divergences that can be obtained continuously from formula (31) in the limit N = 1. For N = 0 the action reduces to the fourth-order quantum gravity, which was the subject of similar considerations in [6][7][8][9], namely We performed the corresponding one-loop calculation just to make an extra check. The linear in curvature divergent part of the one-loop effective action is given by Taking into account the difference in notations, this expression perfectly agrees with the well-known result of [8]. Let us remember, once again, that this expression depends on the gauge-fixing [7]. The expression (37) corresponds to the minimal gauge that is the simplest one for the sake of practical calculations.

IV. QUANTUM NONLOCAL MODELS
It is tempting to use the results for the divergences in the polynomial theory (31), (32), (33) in order to obtain the same type of divergences in a wide class of the nonlocal models of quantum gravity. The most interesting theories are certainly the ones that have no massive ghost-like poles in the propagators. These models were introduced by Tseytlin [48] in the framework of string theory as an alternative to the Zwiebach transformation of metric [49] in the low-energy effective action. Another interesting feature of the original nonlocal model [48] was the absence of singularity in the Newtonian solution for the gravitational field of a point-like massive particle. A ghost-free nonlocal gravitational action was proposed for the first time by Krasnikov in 1988 [50] and studied by Kuz'min in 1989 [51]. Later on, qualitatively similar ghost-free nonlocal models were suggested by Tomboulis in [52] as candidates to be unitary and super-renormalizable or even finite quantum gravity theories [41][42][43]. One can remember that the general aspects of quantization of nonlocal quantum field theories have a long history [53]. For the recent developments in many aspects of nonlocal theories the reader can follow the literature in [41][42][43][54][55][56][57][58][59][60][61] and [62] about the localization of nonlocal theories and nonlocal gravity as a diffusion equation. Recently in the paper of one of the present authors [44] it was shown that an infinite amount of the ghost-like complex states in these models emerge in the quantum theory, when one takes the loop corrections into account.
It is not clear whether the presence of the ghost-like complex degrees of freedom do violate unitarity in these theories, especially in view of the recent negative answer for the models with finite number of complex states [31,34,35]. But in any case it is important to have better understanding of the possible form of quantum corrections in such theories. Let us start the discussion of this problem from the simplest example [48].
Our starting point will be the nonlocal theory of quantum gravity with the classical action [41,[50][51][52] We do not add a possible nonlocal generalization of the Gauss-Bonnet term, but this is possible. The power counting in this theory reduces to the topological relation (6) between the number of loops, internal lines and vertices [44]. Its evaluation shows the result may be identical to the one of the polynomial theory (2) with N 3. This means that the divergences emerge only at the one-loop order and have the form of the cosmological constant and Einstein-Hilbert terms, plus two relevant terms with four derivatives, namely C 2 and R 2 . In order to achieve this, one has to choose Φ(z) and Ψ(z) functions to have equal asymptotic behaviour and be sufficiently fast growing functions at z → ∞.
It is clear that we have to deal only with the model (2) with a very large N , hence the result for the cosmological and linear in R divergences in the polynomial theory are given by the expressions (33) and (32). The most natural option is to choose Ψ and Φ such that the structures of the Euclidean propagator in spin-2 and spin-0 sectors are the same. Consider the simplest version of the theory [63][64][65], Then we arrive at the equation for the pole such that the unique spin-2 pole of the propagator is massless. In what follows we consider only this version of the theory, but it is certainly possible to make generalizations. Let us show, by taking a limit N → ∞ in the results (31), (32) and (33), that the coefficients of the one-loop divergences in the model (38) with functions (40) experience an unrestricted growth. In the polynomial theory the equation for the poles of the propagator is defined by and the same for the function Ψ. Then the relevant coefficients of the polynomial action (2) have the structure The inspection of the coefficients of expressions (31) and (33) shows that they behave like N and N 2 , correspondingly. In the limit N → ∞ both coefficients explode, that demonstrates a discontinuity of quantum corrections for the exponential theory in the given approach. This shows that in genuinely nonlocal theories a new type of one-loop divergences appear (divergences of large N ). Let us note that one can provide the finiteness of the theory for each N in a way we described above, and then the limit N → ∞ is well-defined. For the contributions to the inverse of Newton constant this may be achieved by adjusting the ω N −1,GB and/or ω N,GB coefficients to provide zero result for all N in Eq. (31). These coefficients would explicitly depend on N . For sensible nonlocal theory defined in the limit N → ∞ these coefficients should survive, but the explicit evaluation of them shows that they are decaying like 1/N . In conclusion such generalized Gauss-Bonnet terms are not present in the final action in the limit N → ∞ of the nonlocal theory. One may try to add standard killer operators [42]. However, without explicit calculations it is unclear whether this will be sufficient to provide an UV-finite nonlocal theory.
The expressions of the type ω N −1,i /ω N,i appear in (32), where there is no mixing term. In full generality the diagonality in i, j indices (i, j = R or C) lets us to interpret such ratio in the limit N → ∞ as a d'Alembert definition of the convergence radius 2 of the complex function defined by the series f i (z) = ∞ n=0 ω n,i z n , that is This result can be exported also for the similar interpretation of the divergence proportional to the cosmological constant. From this it is obvious that if any form-factor like Φ or Ψ is an entire function on the complex plane, then the new type of divergence will inevitably appear. However, this formula gives also the possibility of defining the standard divergences for theories in which Φ or Ψ are other non-analytic real functions.
In the first part of this section we studied the simpler example of nonlocality that turns out to be sufficient to make convergent scalar field theories in Euclidean signature. In Minkowskian signature it is sufficient to replace exp(−✷/M 2 ) by exp P (✷/M 2 ), where P (✷/M 2 ) is a polynomial of even degree in its argument. However, in a gravitational or a non-abelian gauge theory case such nonlocality seems not suitable exactly because of the gauge symmetry that forces the kinetic and interaction operators to have the same ultraviolet scaling, which implies nonlocal divergences. However, we can overcome this issue with a special class of nonlocal form-factors that enjoy the property to be asymptotically polynomial in a region around the real axis [41,42,51,52].
Let us give here an explicit example of such function which is supposed to replace the exponential exp(−✷/M 2 ) in (39), where p(z) is a polynomial of degree N , γ E is the Euler-Mascheroni constant and finally Γ(a, x) is the incomplete Gamma function. The crucial property of the function in (44) is that the divergent contributions to the beta functions only depend on the local asymptotic polynomial p(z), while the full nonlocality manifests itself only in the finite contributions to the quantum effective action. It is easy to see that any correction to the UV polynomial is exponentially suppressed [42]. Therefore, all results which were reported for the polynomial case can be applied to this class of theories.
Let us note that the beta functions in the case under consideration are one-loop exact and that two out of the total four relevant divergent contributions to the quantum effective action are known. The last observation is that in the paper [51] it is stated that the beta-functions for R 2 and R 2 µν vanish, but we are unable to confirm or verify this statement.

V. RENORMALIZATION GROUP EQUATIONS
Let us come back to the polynomial theory (2) and consider the case without killer terms, when the divergences are given by the expressions (32) and (33). In this case one can construct the renormalization group equations for the parameters ω EH and ω cc of the classical action. The derivation of these equations is almost a trivial task (see, e.g., [46] for the introduction), but for the sake of completeness we provide the explanations here. The system of renormalization group equations for the running coupling constants is defined by relations Here α i (t) are obtained by making the coupling constants α i scale-dependent in the classical action. The renormalized Lagrangian is where L ct is a counterterm and Z αi denote renormalization constant of the couplings α i . Let us refer to the further details in Appendix B. The final results for the beta-functions for N 2 are as follows: A pertinent observation is that for N 3 these two beta-functions are exact, since higher loops would not provide further contributions. The next question is that whether these beta-functions have physical sense. In other words, we have to describe the situation when they may describe the running of the gravitational and cosmological constants.
It is instructive to compare (47) with the general arguments based on covariance and dimensional considerations, Eq. (1). For this end we have to identify what are the masses in the theory under consideration. In general, the theory of higher derivative quantum gravity (2) possesses many massive degrees of freedom, some of them ghost-like and some normal [30]. Then our Eq. (47) is exactly (1), where the dimensionless ratios between different ω's play the role of coupling constants.
The most interesting is the case when all massive states have complex masses m i . The theory may be consistent only if these masses enter as complex conjugate pairs. This condition is provided by a real classical action of the theory. But the same condition applies also to the renormalization group improved action. Therefore, in this case Eq. (1) boils down to the reduced version µ d dµ where α i are real coefficients depending on the dimensionless ratios of the coefficients ω.
As it was extensively discussed in [21] and in the review paper [66], it is not clear that both or any of the renormalization group equations can be applied to cosmology and astrophysics. The reason is that the diagrams which lead to Eqs. (47) and (48) include internal lines of a massless graviton, but also of the number of massive states. In the most interesting cases of unitary models of higher derivative quantum gravity these massive poles are all complex [31,32]. Then we face an unsolved problem of what remains from the effects of quantum gravity with massive modes at low energy [9]. It is well-known that the contributions of the loops of massive fields do decouple in the IR, in accordance with the gravitational version of the Appelquist and Carazzone theorem [67]. However, this is a well-established result only for the contributions of massive matter fields, when all internal lines of the diagrams have the propagators with equal masses. In the case of quantum gravity one has to deal with a more complicated case of mixed loops, where some of the internal lines are massive and some are massless, or have much lighter mass. Needless to say that the situation becomes much more tricky in the case of numerous complex masses. The analysis of this issue is certainly very interesting, but it is beyond the scope of the present work. In any case the expressions (47) and (48) represent a good starting point, being a universal UV limits for the physical beta-functions under discussion.

VI. CONCLUSIONS
The one-loop calculations always represent one of the most important elements in understanding new theories of quantum gravity. These calculations are especially difficult in the case of higher derivative models. In the present work, the one-loop calculations in the super-renormalizable quantum gravity theory has been extended to the term linear in scalar curvature. Together with the previous derivation [30] in the cosmological constant sector, this enables us to obtain the closed system of renormalization group equations (47) and (48), which may be a basis for establishing the physically relevant running of Newton and cosmological constants in both high-and low-energy regimes.
In our opinion, the following two aspects of super-renormalizable quantum gravity theories make the result especially relevant. First of all, in the models with N 2 the one-loop beta-functions are exact. Therefore, in (47) and (48) we have the first example of exact beta-functions in four-dimensional quantum gravity. Second, some versions of these theories have only complex conjugate pairs of massive ghost-like states, which turn out to be unitary at tree-level and [31,32] and at perturbative level [33][34][35]. The results obtained here cover these cases and, therefore, can be seen as a first example of beta-functions in a class of consistent models of perturbative quantum gravity.
(µ, ν) and (ρ, σ) is supposed. We present the expansions valid on maximally symmetric space-times, with Λ a constant parameter characterizing the background curvature according to the formula Λ = R/4.
We give the expression for the relevant part of the bilinear form for the action (2) 1 |g| where for n 1 Λω n,C + 2n − 1 3 Λω n,GB + 8n + 5 3 Λω n,R + 1 6 ω n,C − ω n,R n +g ρσ ∇ µ ∇ ν 2n − 3 9 Λω n,C + 2n − 1 3 Λω n,GB + 2n − 3 3 Λω n,R + 1 For the gauge-fixing term described in (13), (14), (15) the expansion looks like This expression corresponds to arbitrary values of gauge-fixing parameters, but for the practical calculations we used (18). Let us stress that the final result is independent on the choice of α, β and γ. The expressions (51) and (52) were obtained from the self-adjoint forms of the variational derivative operators after the commutations of derivatives and expansions in powers of Λ.
Due to dimensional arguments in theŨ (0) tensor we would have the appearance of coefficients of the type ω N −2,i /ω N,j (where i, j = R, C, or GB), while we know from considerations in the main text that this tensor would not contribute to divergences proportional to the linear term in scalar curvature. Similarly, we know that in theŨ (2) tensor we find coefficients ω N,i /ω N,j , and again they can not contribute to the desired divergence because we know from (29) that for divergences we can only take the trace of the U tensor (with no interference terms.) Hence we neglect writing the form of these parts of the tensor U and we concentrate below only on the partŨ (1) that gives contribution.
The divergent contribution to the quantum effective action is given by: where where L = Λ/µ, Λ is the cut-off scale and µ the renormalization scale. Moreover, n = 2ω = 2(2 − 0 + ) = 2(2 − ε) = 4 − 2ε =⇒ ln L 2 = 2 4 − n = 1 ε , where n is the dimensionality of spacetime. Then eq. (59) turns into Finally, we give here the divergent contributions to the trace of the logarithm of H ′ operator for the following different values of N in dimensional regularization scheme. For N 2: Tr ln H ′ | div = i ln L 2 16π 2 d 4 x |g|(−1) For N = 1: Tr ln H ′ | div = i ln L 2 16π 2 d 4 x |g| (−1) For N = 0: And for N 2 (divergence to the cosmological constant term): where above by R we mean a general gravitational curvature tensor (Ricci scalar, tensor or Riemann tensor).