Spectral action approach to higher derivative gravity

We study the spectral action approach to higher derivative gravity. The work focuses on the classical aspects. We derive the complete and simplified form of the purely gravitational action up to the 6-derivative terms. We also derive the equivalent forms of the action, which might prove useful in different applications, namely Riemann– and Weyl–dominated representations. The spectral action provides a rather rigid structure of the higher derivative part of the theory. We discuss the possible consequences of this rigidness. As one of the applications, we check whether the conformal backgrounds are preferred in some way on the classical level, with the conclusion that at this level, there is no obvious reason for such a preference, the space S1 × S3 studied in earlier works being a special case. Some other possible properties of the higher derivative gravity given by the spectral action are briefly discussed.


Introduction
In the absence of a completely satisfactory theory of Quantum Gravity (QG), it is important to look for some universal features of the future theory as well as to establish some connections between different approaches towards it. In this paper we look at possible benefits one can gain in combining two such seemingly unrelated approaches to (quantum) gravity as spectral action and higher derivative gravity.
The spectral action approach [1][2][3] appeared as a natural but non-trivial development of the approach to Standard Model based on the methods of non-commutative geometry [4,5]. Though this approach does not address directly the question of QG, one cannot say that it is purely classical as it is spectral by its nature. The main advantage of the spectral action formulation is, from our point of view, that it provides a geometric unification of gravity with Standard Model: both, gravity and Yang-Mills (YM) sector, are treated as "perturbations" of some non-commutative geometry while the coupling to fermionic matter takes a form of the usual Dirac action. The full spectral action is (almost) completely fixed by the choice of the so-called spectral triple, (A, D, H). Here A is an algebra of smooth functions on the generalized spacetime of the form M × F , where M is the usual classical smooth spacetime and F is some finite matrix geometry; D is the relevant Dirac operator defined on the full non-commutative geometry; and H is a Hilbert space, on which all these objects are represented as operators satisfying some conditions (see [5,6] for the details on non-commutative geometry and, e.g., [7] for the non-commutative geometry approach to Standard Model). Thus the full spectral action is given by [2] S spec = Tr χ D Λ + ψ|D|ψ , (1.1) where χ = χ(p) is some positive cut-off function and Λ is some characteristic energy scale. 1 While the Dirac-type second term of (1.1) describes the matter-geometry coupling (including coupling to YM gauge fields) the first term is completely spectral and describes the dynamics of geometry, i.e. gravity plus Yang-Mills. 2 In the case of pure gravity, i.e. using the usual Dirac operator, this term correctly "predicts" the standard Einstein-Hilbert (E-H) action with cosmological constant as the first two terms of some asymptotic expansion. The higher terms in this expansion produce action terms depending on higher orders of curvature. The beauty of the spectral action is that these terms depend only on three inputs: the choice of the Dirac operator D, the scale Λ and the cut-off function χ(p). 3 Also, each level of the expansion is a very precise combination of terms with the same derivative order. E.g., as we will see, the next term in the expansion (after the one reproducing E-H action) has precisely four derivatives of the metric tensor. However, it is worth to emphasize that in principle the series for the spectral action continues up to infinity, so the full (re-summed) theory would contain infinitely many derivatives and therefore would be non-local. The appearance of these higher-order terms (starting with the ones quadratic in curvature) is what makes the direct connection with the higher derivative gravity theories. Such theories contain more (than two) derivatives on metric in their corresponding classical equations of motion (EOM).
As we have discussed above the presence of higher derivative terms (in the form of terms with higher than linear powers of curvature) in the gravitational action is the inevitable consequence of models based on the spectral action principle. The motivation for quantum field theory models of gravitational interactions with higher derivatives comes from many directions. Though there has been a recent revival of ideas related to higher derivative gravitational theories, the idea is not new. It was already observed by Utiyama and DeWitt [8], that the terms with higher derivatives of metric are needed for the renormalization of two-derivative matter theory on a curved background. (Basically, since in the matter part of the theory the couplings are with non-negative energy dimensions, the same energy dimensionality we must require from gravitational counter-terms. In d = 4 spacetime dimensions, the generally covariant terms with dimensionless coupling constants must necessarily come with four derivatives of the metric.) This is one of the reasons for higher derivative models. The renormalization of infinities in matter theory occurs in semiclassical approach when the gravity is taken to be a classical non-dynamical background and when the effects of back-reaction are neglected (so called quantum field theory in curved space).
On a different vain, also on the level of classical analysis, higher derivative theories come with some benefits [9]. When approximate solutions are studied (for example for Newtonian gravitational potential and using the Fourier transform method which presupposes some boundary conditions at spatial infinity) the effective classical spacetime shows no singularity at the origin (in this case there is no Newtonian-like r −1 singularity of the gravitational potential [10,11]). Similar resolution of singularities happens also in cosmological framework [12,13]. However, these results have to be taken with a grain of salt, since they were obtained only in the perturbative scheme. Moreover, even in the linear approximation regime, the set of solutions of the higher derivative analogue of the Laplace equation still contains the singular r −1 solution (as this was pointed out in [14]). Consequently, also the class of exact solutions is bigger (than in the standard two-derivative theory). And there typically among standard singular solutions (known from the two-derivative Einstein-Hilbert theory) we find also new types of solutions: often run-away or non-singular solutions. And these new solutions are responsible for resolution of singularity or give rise to cosmological inflation [15][16][17][18], where apparently they are not unwanted solutions anymore. But we emphasize once again that in the whole set of exact solutions of higher derivative gravitational theories we have both singular solutions and new solutions, and the proper choice between them is specified by the boundary conditions. We just remark that in higher derivative theories we need to provide more boundary (or initial for dynamics in time, so for evolution) conditions than just in two-derivative theory, because of the increased number of derivatives on the metric field.
However, the main reason for higher derivatives in the gravitational setup has to do with quantum physics. As already noted by Stelle [19] the quantum gravity theory based on the four-derivative action in d = 4 behaves much better regarding the situation with UVdivergences. Stelle's quadratic gravity is the first model of multiplicatively renormalizable gravitational theory. The quantization here can be performed using the standard methods in a fully covariant fashion (unlike in the more recent re-incarnation of the higher derivative gravity in the form of Hořava-Lifshitz gravity [20], which explicitly breaks 4-dimensional diffeomorphisms, 4 and where the higher derivatives are only in the spatial part). One sees that the inclusion of higher derivative terms as bare terms of the action is sufficient for absorption of divergences present in gravitational one-loop counter-terms [26]. This is intimately related to the observation made in quantum field theory in curved space, where we need the same counter-terms but there they are generated from the matter side [8]. Using the Batalin-Vilkovisky formalism [27,28] the proof of renormalizability can be extended to all perturbative loop orders. This means that all infinities appearing in the perturbative calculations can be consistently taken care of by re-defining the bare couplings of the theory which are present already at the tree-level. As another advantage one notices that such theory is predictive and only these couplings, which are present in the tree-level action, require experimental input, but going to higher loops does not force us to introduce new couplings at all. This is in strong distinction with quantum gravity theory based on E-H action, which is non-renormalizable from two loops on (at one-loop level without matter and on-shell this theory is miraculously finite), hence this theory is very badly non-predictive.
Moreover, this is not the end of good points about higher derivative models. In the quadratic gravity of Stelle the renormalization group (RG) behavior of the essential couplings of the theory can be studied. As it was first discussed by Fradkin and Tseytlin [29,30] in d = 4, the theory shows asymptotic freedom in all such couplings. This in turn signifies that the bare values for all couplings vanish (in the ultra-violet (UV) limit) and there is no issue of initial values for them, very similarly to the situation in QCD. Furthermore, the analogies of quadratic (in gravitational curvatures) gravity with quadratic (in YM field strengths) matter theory 5 go even further regarding, for example, the form of scattering amplitudes [31]. Both theories are renormalizable and both are asymptotically-free. Finally, if one couples such matter theory to the quadratic gravity then also the total quantum system is described by renormalizable and asymptotically-free theory, where all the UV-infinities are under control. Such generalized framework provides a very interesting quantum laboratory for study of grand unified theories (GUT) coupled to gravity as field theory models of gravity-gauge unification (compare for example [32]). In the connection with the spectral action approach, it is worth noting that from the point of view of (1.1) both, YM and four-derivative gravity terms, appear exactly at the same level of the asymptotic expansion of the spectral action. Namely, they come from the Seeley-DeWitt coefficient a 4 (see the next section 2).
The four-derivative models can be also further generalized. Inclusion of terms with more derivatives is one such direction. When 6-derivative terms are added and enough care is exerted the gravitational model can still be shown to be renormalizable. But even more, for a particular form of generalization, the theory reveals to possess bigger control over perturbative divergences and it is said to be super-renormalizable. As proven by Asorey, Lopez and Shapiro [33] for the case of a theory with six derivatives in d = 4, all the diagrams with more than 3 loops are perfectly UV-finite. (The increase of number of derivatives can lower this bound and for example for theory with 10 derivatives in d = 4 spacetime dimensions the UV-divergences remain only at the one-loop level [34].) The requirement of renormalizability constrains the possible terms, which could be added to the action of a six-derivative theory. First, in the action there must be terms quadratic in curvature in order to have a highly improved behavior of the flat spacetime propagator in the UV-regime. On the other hand, we can add also terms containing three curvatures but for keeping the dimensionality of these terms under control they must be with no covariant derivatives in their generally covariant construction. This is in agreement with the scheme provided by effective field theories for gravity, where after inclusion of terms with four derivatives we, in principle, in our effective action should include all possible terms with six derivatives (partial, of the metric as seen in EOM). The possible terms can be quadratic or cubic in curvatures. We notice that such a theory naturally arises if we study the further expansion of the spectral action, namely the a 6 coefficient. However, such higher derivative gravitational theory based on the expansion of the spectral action up to a 6 coefficient has very particular coefficients in front of curvature invariants, which are dictated by spectral principle. And therefore this special theory can be viewed as an example of a general six-derivative theory in d = 4. In this sense we exploit here the predictive power of the spectral principle, since all these coefficients are unambiguously determined (up to the overall scaling, see in the next section 2). The main purpose of this paper is to study some (if any) consequences of this special form of the higher derivative (up to six-derivative) gravity.
The plan of the paper is as follows. In the next section 2 we discuss the expansion of the purely geometric part of the spectral action (given by the first term in (1.1) for the standard choice of the Dirac operator) up to six derivatives of the metric. Section 3 contains our main results. It is devoted to the detailed study of the obtained higher derivative action in two different "bases": Weyl-and Riemann-dominated. We also obtain the corresponding equations of motion. In section 4, on the example of the beta function for the cosmological constant, we initiate the study of the consequences of the spectral action approach on the quantum level. We conclude with the extensive discussion of further possible consequences of our approach as well as delineation of further steps. Because one of the goals of the paper is to bring together the communities working in the areas of higher derivative gravity and non-commutative geometry, we include several appendices, where we provide the technical details of the construction (even though sometimes they are pretty standard) to make the paper as much self-contained as possible.

Spectral action
As we mentioned in the introduction, the full spectral action describes both, geometric and matter, sectors. Because we are interested in the case of pure gravity, here we will consider in details only the first term of (1.1) for the standard choice of the Dirac operator. Before we proceed, the word of caution is in order. The whole construction is well-defined (or rather well-understood) only for the case of compact Riemannian manifolds. Hence, the final results for the Lorentzian signature should be understood as an analytic continuation of the results in the Euclidean framework (preceded by taking some decompactifying limit). Or, another point of view could be taken: because all the expressions are written in terms of the geometric invariants, one can continue formally to use them for the pseudo-Riemannian case (only carefully keeping track of various signs).
First of all, we have to fix our Dirac operator. The choice of the standard one could be motivated from two different points of view. The idea of the spectral action is that the same Dirac operator is used in both terms of (1.1). This was used to some extreme in [24,25] to relate the free parameters of the matter and gravitational sectors in Hořava-Lifshitz gravity. Because in the present paper, we do not want to modify the minimal coupling of matter to gravity and this is given by the Dirac-type second term in (1.1), we must immediately conclude that we have to use the same, undeformed standard Dirac operator on the geometric side. On the other hand, there is a very powerful result in the spectral (a.k.a. non-commutative) approach to commutative geometry, due to Connes [35], which, for our purposes, could be stated as follows: there is one to one correspondence between compact Riemannian geometries and the commutative spectral triples (A, D, H) defined by the usual Dirac operator. All of these arguments lead us to the following choice (see appendix B for the notations and details of some relevant calculations) where ω µ is the usual spin-connection. Then the object we want to study is given by the first term in (1.1) It is obvious that the spectral action (2.3) is a highly non-trivial object, in particular it is non-local. So, is there any chance to have an explicit form of it? As we review in the appendix C, this could be done at least in some limit, using the so-called heat kernel expansion, which corresponds to the derivative expansion of (2.3). Using the general result (C.47), we can easily write the asymptotic expansion for (2.3) in d = 4 Euclidean space dimensions Here a 2k (D 2 ) are Seeley-DeWitt coefficients for the elliptic operator where E is some endomorphism of the corresponding bundle (see below and appendix B).
To determine these coefficients, the general method, due to Gilkey [36,37], or due to DeWitt [38] can be used. The essence of the former method could be roughly described as follows.
Because the operator (2.6) transforms covariantly under the diffeomorphisms of space and, in general, twisted spinor bundle endomorphisms, it is possible to show that all the coefficients in the heat kernel expansion (C.41) (or (2.4) for the case of the interest) are given by the space volume integrals of the local geometric invariants, a n (x), and subsequent traces over the bundle indices (in the pure Riemannian case, this is a spinor bundle). Moreover, the whole explicit dependence on the dimension, d, of the manifold M is given by an overall factor, according to the formula: which is essentially the formula (C.42) from the appendix C. In a sense, the volume integrals of the expressions traced over internal indices can be viewed as the result of taking functional traces, both in internal and external space at one stroke. Therefore this generalizes the notion of trace to include also volume integrals and is in accord with the DeWitt convention [38] for compact index notation, treating on the same footing both spacetime M and internal space. The local invariants, a n (x), are constructed from the only geometric objects available: the Riemann curvature tensor R µνρσ , the endomorphism E (B.23), the "field strength" of the bundle connection Ω µν (B.24) and the covariant derivatives ∇ µ . Here we use only Levi-Civita covariant derivatives ∇ µ , in opposition to total covariant derivative with spin-connection ∇ ω µ , because all geometric objects available are proportional to the identity 1 in the spinor indices space (see also later and in the appendix B). What invariants enter at each order n of the expansion, a n (x, D 2 ), can be easily decided from the dimensional analysis, after assigning the standard dimensions ([x µ ] = −1, [g µν ] = 0 and the rest follows) and requiring that the exponent of the heat kernel, tD 2 , in (C.33) is dimensionless. From this point of view, it is clear that the heat kernel expansion (C.41) is a derivative expansion (rather than the curvature expansion -by integrating by parts, one can always trade some curvature for derivatives, see also below). In this way, each a n (x, D 2 ) is a linear combination of the local geometric invariants, each having exactly n derivatives, and the coefficients of this combination are universal, i.e. do not depend neither on geometry nor on the dimension d. This allows to fix these coefficients by evaluating the heat kernel on special geometries characterized by some degree of symmetry (like tori, spheres, etc.). This is the essence of the Gilkey method. However, the method by DeWitt is more algorithmic but more tedious since all these coefficients (generalized Schwinger coefficients) are obtained by differentiation of the world-line function and taking the coincidence limits [38,39]. Just a word for terminology: we will call by a n (x, D 2 ) the unintegrated coefficients, while the coefficients a n (D 2 ) as in (2.4) are already integrated (and traced over). Using any of these methods, one can find the following expressions for the several first unintegrated Seeley-DeWitt coefficients [36,40]: where R is understood as R 1 and we already used some trivial simplifications compared to [36,37] due to the fact that the endomorphism E is proportional to the identity 1 in the spinor space, that allowed to combine some terms, which otherwise are not equal. To further simplify these general formulas for a n (x, D 2 ) we use the explicit expressions for E and Ω µν , from (B.23) and (B.24), as well as some standard trace identities for the Dirac gamma matrices (for our notations, see the appendix A): Tr γ µν = 0 , Tr(γ µν γ ρσ ) = −2 Tr(1)g µ[ρ g σ]ν , (2.12) Using these relations, one can easily establish With the help of these relations, it is straightforward to re-write (2.8-2.11) in the following form a 0 (x, D 2 ) = Tr 1 (2.14) where we have used the notation GB 0 for the standard Gauss-Bonnet term (see the appendix A for more details). The above derivation is valid in general dimension d.

Spectral action in different representations
Combining (2.14-2.17) and the general result for the asymptotic expansion of the spectral action (2.4), we arrive at the following action for the 6-derivative gravity: where we already integrated by parts and discarded all the surface integrals. In fact, as we commented in the appendix C (see the footnote 14), we could do it because the expansion (2.4) is already written in the assumption that there is no boundary. Also we explicitly kept the dependence on the cut-off scale Λ, while calling the numerical coefficients for each level with 2k derivatives by µ k (which anyway are not fixed by the model). In this regard, three comments are in order.
1) It is clear either from (2.4) or from (3.18) that the level with 2k derivatives is suppressed by the factor of Λ −2 compared to the one with 2(k − 1) derivatives. This seems very natural except that the term without derivatives should correspond to the cosmological constant term, which would lead to the huge dark energy density. We will not address this point in this work and rather refer to [3] where the way to resolve this problem is discussed.
2) As we said above, the numerical coefficients µ k are not specified within this approach and, in principle, are free parameters of the model. Still one "prediction" can be made. Let us recall that the arbitrary function χ(p) in (2.3) is supposed to be some kind of a cut-off function and the coefficients µ k are proportional to f 2k given by (2.5) in terms of this cut-off function. From this it is clear that while µ k , k = 0, 1, 2 are really arbitrary non-zero numbers, µ 3 , being proportional to χ (0), should be zero if χ is flat at the beginning of the spectrum (as it is in the case of the standard cut-off functions). As we commented above, the actual shape of χ should be fixed by some fundamental theory, but we might expect that the behavior of χ in IR is very close to the flat one, i.e. that χ (0) 1. This is due to the fact that the spectrum of the standard Dirac operator controls the classical geometry of spacetime and we do not want to distort this spectrum too much (by modulating it with χ) in IR. Hence, based on this discussion, the prediction of the model would be the additional suppression of the 6-derivative term by the small factor µ 3 1 (in addition to Λ −2 suppression).
3) The real prediction of this approach is given by the values of the coefficients within each derivative level (so the relative weights of terms). This drastically reduces the number of free parameters in higher derivative gravity. E.g., without spectral action, the number of free parameters for higher derivative gravity with up to six derivatives would be well above 10, while in our model we have just 4 (these are µ 0 , µ 1 , µ 2 , and µ 3 respectively).
The main goal of this section is to maximally simplify the action (3.18) and present it in several equivalent forms that might be useful for different types of problems.
As the first step, let us make the most obvious simplifications related to the following terms R µν R µ ρ ;νρ and R µν R ;µν in (3.18). One trivially has where we used the standard result for the commutator of the covariant derivatives (with the sign conventions from appendix A). Now using the doubly contracted second Bianchi identity (A.8), R µ ν;µ = 1 2 R ;ν we have for the integrals (again, discarding total derivatives): The equation (3.19) is a typical example of how we can trade derivatives for curvature, so we again stress that the heat kernel expansion is a derivative expansion rather than the curvature one. Now we would like to get rid of the term R µνρσ R µνρσ . Note that this term enters the generalized Gauss-Bonet term, GB 1 (A.14). Repeatedly using commutators of the covariant derivatives and the contracted second Bianchi identity, as in the derivation of (3.19) and (3.20), one can easily get So far for our simplification of the action functional we used identities valid in any number of dimensions d. Final steps are done with the assumption d = 4. Then further simplification is possible due to a very simple observation: in 4 dimensions, anti-symmetrizing any tensor with respect to five or more indices identically gives zero. Choosing different products of three Riemann, Ricci or Weyl tensors, this leads to the following identities [41,42] or equivalently in terms of the traceless Weyl tensor Combining these relations with (3.19), (3.20) and (3.21) one gets, after some straightfor-ward calculations, the following compact result for the action (3.18): where compared to (3.18) we changed µ 3 → 1008 µ 3 (though, as we mentioned, this is quite irrelevant taking into account that µ 3 is a free parameter), also we kept the topological term GB 0 even though it will not contribute to the classical equations of motion (and, of course, it should be kept for quantum calculations). This is exactly the form that we call the action in the Riemann basis or the Riemann-dominated action. It is clear that the Riemann-dominated form is not the most convenient one if one wants to study the conformal backgrounds. This motivates us to look for the equivalent expression, but now written in the Weyl basis or in the Weyl-dominated form. This is readily done by expressing most of the terms in the Riemann-dominated action (3.23) with the help of the definition of the Weyl tensor (A.9) and the relation (A.11) and, when necessary, again using (3.21) and the identities (3.22). After not so lengthy and straightforward manipulations we arrive at the result: Yet another form might be useful for comparing with other works on higher derivative gravity (see, e.g. [33,43]) where the action is written in (R, C, GB)-basis: where V (R) is some "potential" depending on the curvature R. 6 The case V (R) = 0 would 6 As a generalized curvature R here we understand any tensor constructed from Riemann tensor R µνρσ by various contractions.
correspond to the "minimal" action in this basis (but compare with the discussion in [34,44]). But one should remember that in view of (3.21) it is not preferred in any other sense. In any case, the spectral action in this basis takes the "non-minimal" form with After obtaining the simplified forms of the 6-derivative gravity coming from the spectral action, it is reasonable to ask whether we have any further advantages of the approach beyond the obvious rigidness of the result (in the sense of the great reduction in the number of free parameters). In other words: is there anything special about the spectral action approach from the point of view of higher derivative gravities? In this section we will touch on this point on the classical level and in the next we briefly discuss some quantum aspects postponing a more detailed discussion of the quantum case for the future work.
One may guess that the spectral action has a lot to do with additional symmetries present in the gravitational interactions. The conformal symmetry may play such a role. Actually, based on the explicit results of the a 4 coefficient (last part of the first line in (3.24)), one would be almost convinced about this since in d = 4 dimensions the only non-trivial term appearing there is the C 2 term, which transforms in a covariant way under local conformal transformations. (We neglect here the Gauss-Bonnet term GB 0 since this is a topological term in d = 4.) Only in four dimensions, in a 4 we have only C 2 and GB 0 terms, in other dimensions there is a non-zero coefficient in front of the R 2 term. Moreover, in a 2 we have only a term with Ricci scalar R and this is exceptionally conformally covariant term in the action of gravity in dimensions d = 2. This hope is reinforced by the fact that the √ gR 2 term is missing in a 4 in d = 4 and this term is only globally scale-invariant in four-dimensional case (invariant only under rigid scale transformations) and hence dimensionless. In a 6 we naturally have terms with six derivatives, so they cannot be dimensionless in d = 4, but they might transform covariantly (that is with a weight factor) under conformal transformations.
The condition for this is that they would have to be built out of only Weyl tensor and its various contractions and no covariant derivatives or covariant box operators acting on these conformal tensors [45]. Then they would be truly conformally invariant in d = 6 dimensions and there they would be therefore dimensionless. However, the inspection of the action written in the Weyl-dominated form (3.24) shows that this hope for additional symmetry of the spectral action is not fulfiled. We find there, in the sector of terms with six derivatives, terms built also with Ricci scalar (which does not transform conformally in a neat way) and even terms of the type C µνρσ C µνρσ , which would break conformal symmetry in d = 6. Based on the explicit example of the a 6 coefficient we conclude that generally conformal symmetry (even in a restricted sense in d spacetime dimensions for the a d coefficient of the expansion) is not a feature of the spectral action approach.
As the first application of the simplified action in the Weyl basis, let us evaluate the action (3.24) on the conformal background, i.e. when C µνρσ = 0 7 . We want to compare this with the discussion in [3] where this was done for the special case of S 1 × S 3 background. Even for this case, the calculation was extremely complicated technically and the main result (that the 6-derivative part of the action is zero for this background, see below) was very surprising.
Here we re-derive this result, which will also provide an independent check of our action (3.24), and discuss what can be said in the case of a general conformal background. In this way we generalize the results from [3] including the impact of terms with covariant derivatives on curvature tensors. In the remainder of this subsection we analyze this issue, while in the next subsection we analyze whether some commonly known background spacetimes are exact solutions of the theory.
Let us trivially evaluate (3.24) for the geometries with C µνρσ = 0. 8 The result is Already from this result it is obvious that in the case of a general conformal background, i.e. when C µνρσ = 0, the 6-derivative part of the action will not be zero. This makes the result for S 1 ×S 3 even more surprising. Our general result (3.27) allows to obtain it almost trivially compared to [3]. First of all, because this background has a constant scalar curvature (see (3.29) below), the term R R drops out automatically. The only non-trivial components of Riemann tensor are

28)
7 By conformal backgrounds we mean backgrounds which are conformally flat, that is by conformal transformation of the metric tensor g µν → g µν = Ω 2 (x)g µν with some suitable function Ω(x) we get the metric g µν as the metric of flat spacetime, i.e. the Riemann tensor of the g metric vanishes identically. The condition for conformal flatness in dimensions d 4 is equivalent to vanishing of Weyl tensor C µνρσ . Hence this last tensor is also called as the tensor of conformal curvature. 8 Note that for the calculation of the equations of motion on the conformal background, one cannot just set all the Weyl terms to zero. This is because the variation of the Weyl tensor evaluated on the conformal background is not zero, so one has to keep the terms linear in C µνρσ , see below.
where a is the radius of S 3 and the space-like indices i, j, k, l = 1, 2, 3. Contracting, we get Ricci tensor and the scalar curvature (pay attention to our sign convention in (A.1), (A.2) and (A.3)) R ij = 2 a 2 g ij and the rest are zero, Using this, one easily calculates the relevant terms in (3.27).
Combining these results and using them in (3.27), it is trivial to see that the 6-derivative term is zero, while the whole action evaluated on this background is given by where 4π 3 a 3 b is just the volume of S 1 × S 3 with b being the radius of S 1 and a of S 3 . The equation (3.31) is essentially the result obtained in [3] by the direct evaluation of a 6 (2.11) for the S 1 × S 3 background. Thus our approach correctly reproduces this result and demonstrates the role of the performed simplifications. Also, we want to stress one more time that the cancellation of the 6-order terms for this background should be considered as accidental: it is not automatical but happens due to the non-trivial cancellation between terms depending on the curvature. Because this happens exactly for the coefficients fixed by the spectral action, one might speculate that the spectral action somehow prefers this background.

Equations of motion
Now let us make one step further and derive the equations of motion for the Riemannand Weyl-dominated forms of the action, (3.23), (3.24). Though the general equations of motion following from (3.23) are not very illuminating, in the appendix D we give the final result for them for the possible future references and applications. Below we will consider a special case of these equations for the very important case of the Ricci-flat backgrounds.
In the case of GR, Ricci-flat geometries are the special types of the Einstein spaces for the case of zero cosmological constant. The most known (and, probably, the most important) of these solutions is the Schwarzschild one. In the framework of the higher derivative gravity, one would like to find the corrections to this solution (and extract from there the quantumgravitational corrections to Newton's law). Postponing this very important task for the future research, here we just show how the simplified form of the action (3.23) easily allows to derive a compact set of the equations of motion for a general Ricci-flat background. Also we verify that, not surprisingly, the Schwarzschild metric is not a vacuum solution of these equations and rather requires as a source the energy-momentum tensor with very peculiar, exotic and unphysical properties.
To derive the vacuum EOM (where we do not include matter energy-momentum tensor on the RHS) for R µν = 0 case from (3.23) we note that one can set to zero in (3.23) all the terms that are more than linear in R µν and R (but not in R µνρσ !). The linear terms should be kept. One should also drop the GB 0 term. This immediately kills almost all the terms in (3.23): As we said above, not surprisingly, the standard Schwarzschild spacetime with a metric tensor in standard Schwarzschild coordinate system given by (by dΩ 2 , as usual, we denote angular part of the metric, that is dΩ 2 = dθ 2 + sin 2 θdφ 2 ) is not a vacuum solution to the equations (3.33). We find for the respective components  Utilizing (3.35) one can come up with the corresponding plots as shown in Fig. 1. There we exploited the dimensionless variable X = r M and rescaled the components of EOM in (3.35) by a common power M 3 . It is evident that zeros for the three components E t t , E r r and E θ θ = E φ φ are given by approximate values of the X coordinate X ≈ 2.2, 1.55 and 2.34, respectively. We observe that energy density E t t and azimuthal pressure E θ θ = E φ φ both are positive for considerably higher radii. Whereas, the radial pressure E r r becomes negative at X 1.55. So we see that to get (3.34) as a solution, (3.33) must be sourced by a very non-physical energy momentum tensor (the cosmological constant is set to zero).
One can check that the components of the effective energy-momentum tensor (as evaluated in Eqs. (3.35)) do not satisfy energy conditions (neither strong, dominant, nor null one). This feature is actually common to almost all higher derivative theories since this is the price to have among bigger set of solutions also those which are non-singular. (This is a caveat to the celebrated Hawking-Penrose theorems about inevitability of spacetime singularities -in theories with higher derivatives classical energy conditions are violated and that is why singularities can be avoided in some exact solutions of such theories.) However, from the point of view of effective theory two aspects are worrisome here. First is that for X 1.55 the radial pressure E r r attains negative values. This characteristics cannot be accepted as pertaining to an effective matter source since there does not exist any type of classical matter which exhibits negative pressure. Some exotic examples are brought by quantum effects (vacuum polarization effects) or vacuum energy realized for example as a cosmological constant source (or by Casimir effects). The second problem is not so severe since it touches on the behavior for smaller radii and for two different components of the effective energy-momentum tensor, namely for E t t and E θ θ = E φ φ . We find that these components become negative inside the core of our solution. Hence our solution cannot be a physical representation of a star in higher derivative gravitational theories. This latter issue is not so problematic because these effects happen roughly under the classical Schwarzschild horizon, which is located at X = 2. Even in Einstein-Hilbert gravitational theory the source of the gravitationally collapsing configuration (producing eventually a black hole) inside the Schwarzschild horizon is not a stationary matter source and energy densities there may be negatively valued.
Of course, what one should do, instead of just checking that (3.34) is not a vacuum solution of this version of the higher derivative gravity, is to look for the corrections to the Schwarzschild metric following from the full set of the equations (D.53) (now one cannot use the Ricci-flat ansatz). But this is technically quite involved and will require some numerical study. We are planning on returning to this in the future research.
Analogous analysis can be performed for the Weyl-dominated form of the action (3.24) in the case of a conformal background. Again, in this case one can keep in (3.24) only the terms up to the first order in C µνρσ (and one can drop GB 0 as it will not contribute to the equations of motion): The simplification is not as radical as in the Ricci-flat case (3.32) but still one can straightforwardly find the corresponding equations of motion (the most relevant variation, δC µνρσ , is given in (D.52)): One of the most important conformally flat backgrounds is the cosmological Friedmann-Lemaître-Robertson-Walker (FLRW) spacetime (which is conformally flat for any value of the FLRW topology index k = −1, 0, +1 [46]) (3.38) Plugging this into (3.37) one obtains the following non-zero components for the tensor of equations of motion E α β (only the diagonal components are non-vanishing): (3.39) One notices that the EOM evaluated on FLRW background do not depend at all on the coefficient µ 2 . This is actually true on any conformal background. The reason for this is that the term in the expansion of the spectral action proportional to µ 2 is precisely with four derivatives and as found in [4,5] it is exactly conformally invariant √ gC 2 term, see (3.24). Hence there is no contribution of the first variation of it on the conformal background. From this one derives that conformally flat backgrounds are exact solutions in the same way (we mean that exactly the same form of the source is needed) as in two-derivative Einstein theory with a cosmological constant, when the expansion to the order of a 4 is retained. However, as seen from above equations, the inclusion of the next term in the expansion -the a 6 coefficient changes this conclusion and we get many non-zero terms proportional to µ 3 in EOM. This means that we cannot rely on cosmological solutions of Einstein-Hilbert theory possibly with a cosmological constant and the set of Eqs. (3.39) has to be solved anew. For example, the question whether one can find some well-behaved solution for the scale factor a(t) for some reasonable cosmological energy-momentum tensor requires further study.

Some quantum properties of the model
As we saw in the previous section, using just the classical analysis of the derivative expansion of the spectral action it is difficult to see if there is anything special about it. As a next step, it is very important to check whether the spectral action is preferred on quantum level. In this section we make some initial effort in this direction, postponing the detailed study for the future research.
To proceed, one should quantize the higher derivative theory given for example by the expansion of the spectral action up to the a 6 coefficient. (The quantization of the action up to a 4 was already considered by Stelle in [19] since the model resulting from this level of expansion of the spectral action is a four-derivative theory only with C 2 and GB terms (without R 2 term), with Einstein-Hilbert term R and a non-zero cosmological constant.) For the covariant quantization one can use the method presented in [32] or one may desire to use Batalin-Vilkovisky formalism from [27,28] to have better control over remaining BRST symmetry of the quantized theory. Since the special attention was paid to conformal backgrounds in previous studies (in particular to 4-dimensional product manifolds of the type S 1 × S 3 ) in [3][4][5], it seems natural to investigate the quantum stability of perturbations around these backgrounds. In some minimal sense one should check the positive-definiteness of the quadratic operator governing the dynamics of small quantum perturbations around a conformal background. This issue is tightly related to the positivity of beta functions in front of R 2 and C 2 invariants in the form of the one-loop divergent effective action in the theory. In Euclidean framework both these curvature invariants are positive-definite. In a bigger generality one could consider the whole system of beta functions for the quantum theory, not only in front of dimensionless (in d = 4) terms R 2 and C 2 , but also the beta function of the cosmological constant β cc and the beta function β G of the Newton's constant coupling G N . This last beta function is defined as the divergent coefficient in front of Ricci scalar term R in the one-loop divergent effective action. Actually, for the last two beta functions β cc and β G we know the answer in general higher derivative theories. The easier computation of the beta function β cc was first done in [33], while the more involved computation of β G involving contributions from generalized Gauss-Bonnet terms was achieved in [43]. The analysis presented in [33] and [43] is generally valid on any background spacetime but obviously very easily we can restrict it to a preferred conformal background or even to a particular example of S 1 × S 3 manifold.
First, one can understand that terms cubic in curvatures do not contribute to the beta function β cc . This statement is based on the argumentation presented in [34,43,47]. One may say in simple words that all the terms in the "potential" depending on the curvature V (R) do not influence at all the beta function of the cosmological constant in the theory. However, it is expected that they will contribute to two beta functions of dimensionless couplings β R 2 and β C 2 as well to β G . Actually, the computation of the two remaining beta functions β R 2 and β C 2 is one of the very important goals of the extension of the project, which we plan to address in the nearest future. When one checks the actual expression for the beta function β cc , one sees almost no speciality of the theory based on the action (3.18).
Since this beta function is completely insensitive to terms cubic in curvature in (3.18), see below, we can concentrate only on terms quadratic in curvature. These terms in the action are, of course, very important for defining the kinetic operator and hence ensuing propagator for gravitational quantum fluctuations around flat spacetime background. If we could see any extraordinary behavior of the system of beta functions here, we must emphasize that this would not be a virtue of spectral action approach only since the latter constrains tightly also the numerical coefficients in front of cubic terms, but β cc does not depend on them. Instead the special behavior could be associated to hypothetic structural relations between the terms quadratic in curvature describing kinetic part of the theory. For analysis of the beta function we can use either the Weyl-dominated basis (3.24) or the (R, C, GB)-basis (3.26). Following discussions in [33,43], we note that for the beta function we need to focus on the coefficients in front of the terms with respectively two derivatives, four derivatives and six derivatives being also quadratic in curvatures. Towards this end, let us write (3.24) or (3.26) in the form (3.25) where We notice right away that the beta function does not depend on the cosmological constant term c −2 (because there are no derivatives in this term -and the difference in energy dimensionalities of this term compared to the coefficients c R 1 and c C 1 is bigger than the number of dimensions d = 4), neither on c GB 0 (because this is a topological term in d = 4), nor on c GB 1 (because this term can be re-written in terms that are cubic in curvatures (3.21)). The result for the beta function β cc from [33] (where we have to take N = 1 corresponding to the theory with six derivatives) reads (4.42) Plugging into this (4.41), one finds (4.43) We comment on some simplification which occurred above. First, the second term in the square bracket in (4.42) is not present since c R 0 = 0 as this was discussed to be a feature of the spectral action to the order a 4 in expansion. However, the vanishing of the first term proportional to c −1 is a genuine feature of the coefficients appearing in the expansion to the level of a 6 . The relation between c R 1 and c C 1 (that is c C 1 = 15c R 1 ) is dictated by spectral action approach, but as we emphasized above cubic terms in V (R) do not participate, so right now we cannot judge whether this is a mere numerical coincidence or some deeper fact related to the roots of spectral action and non-commutative geometry approaches. We put importance to the fact that this relation holds independently of the value of the dimensionful cut-off parameter Λ as well as of the arbitrary and adjustable value of the dimensionless coefficient µ 3 . However, from the field theory point of view, there is not much of importance of this observation, since the total beta function is non-zero. As far as we know there does not exist any clear interpretation of the fact that the final expression for the cosmological constant beta function is independent of the value of c −1 coupling, which stands in front of the Ricci scalar in the action (4.40). The final expression for the beta function (4.43) shows that it is always negative-definite and that it depends on the value of the ratio of the coefficients µ 2 /µ 3 only. There is a very little amount of speciality of the quantum behavior of the spectral action.
We also remark that using the analysis presented in [43] we cannot unambiguously determine β G since in our model (3.26) we have other terms cubic in curvature besides the generalized Gauss-Bonnet term GB 1 , while the analysis of [43] was done for the "minimal" model with V (R) = 0. Some preliminary results indicate that there is no exceptional behavior of the other three beta functions of the theory, that is we do not find any of β R 2 , β C 2 and β G to be zero or to be always strictly positive though some further analysis is still required.

Discussion and conclusions
In this paper, we studied some classical aspects of a specific higher derivative gravity theory motivated by the spectral action approach. One of the aims of this work was to bring attention of the researchers working in higher derivative gravity to the methods of noncommutative geometry. This goal partly defined the style of the paper -along with the original research, it contains some details (mostly collected in the appendices) known to those who work in noncommutative geometry but mostly unfamiliar to the higher derivative gravity community. One of the main motivations to consider the spectral action as the basis for the effective higher derivative gravity is the fact that the derivative expansion has a fixed structure within each derivative level, greatly reducing the dimension of the parameter space. This gives hope that the spectral HD gravity might have some special properties compared to the general case. The immediate analysis is difficult due to the very "bulky" form of the general expressions for the relevant terms in the expansion of the spectral action. So, the major part of the paper is devoted to deriving the most compact form of the 6-derivative part as well as some equivalent representations, which might be useful for different types of problems. The formulas (3.23), (3.24) and (3.26) constitute ones of the main technical results of our work.
As we mentioned above several times, the rigidity of the structure of the higher derivative terms gives hope that the theory might possess the features absent in the general case. This hope is somewhat supported by the observation made in the paper with the title suggesting the existence of such special features -"The Uncanny Precision of the Spectral Action" [3]. There it was shown that on the special type of a conformal background, S 1 × S 3 , the higher derivative part of the action is identically zero. Our result (3.24) allowed us to study this point in great generality. We showed that the result of [3] does not hold for a general conformal background, so it is, in some sense, accidental (or signalling that S 1 × S 3 background is in some way special). In particular, the action is not trivial for one of the most physically relevant conformal backgrounds -cosmological spacetimes. To continue the study of the classical gravity based on spectral action, we derived the general equations of motion, as well as their special cases -for Ricci-flat and conformal backgrounds. While the general EOM, do not seem to be particulary simple, in the Riemann/Weyl basis there are serious simplifications. As an (somewhat trivial) application, we explicitly demonstrated that neither Schwarzschild nor cosmological spacetimes are the exact solutions of these equations for the same matter energy-momentum source as this was in standard Einstein gravity.
At this point, it seems that the main conclusion of the classical analysis is that at this level there is nothing much special about the specific values of the parameters fixed by the spectral action. While this appears to be the case, it does not mean that the same should be said in general. The reason is that there is still a chance that the special values of the parameters will be important at the quantum level. In this work, using the example of the cosmological constant beta function, β cc , we briefly touched upon the possible implications of the spectral action approach on quantum level. But much more detailed study is still needed. So, naturally this should be one of the most urgent next steps in the continuation of this project.
We want to discuss here the issue of the dependence of the spectral action on the order of expansion. We can now compare results (both classical and quantum) in higher derivative gravitational theories based on the action given up to a 4 and a 6 coefficients. One sees that inclusion of terms with higher number of derivatives changes theory quite dramatically. For example, in the domain of classical exact solutions (and their stability properties) we observed a lot of differences between the two models as discussed in Section 3. One may ask whether the inclusion of terms with six derivatives of the metric tensor present in a 6 is a small perturbation added to the system. From the field theory viewpoint, this is not the case. Classical EOM change their character from forth to sixth order in derivatives and this implies that we have two new families of solutions for each problem. One might think that however, the perturbation by sixth derivative term is small and it modifies the known solutions (from four derivative theory or even from Einstein theory) only by a little. But due to the higher derivative character of modification we see strong differences both in the IR (long wavelengths) as well as in the UV-regime (boundary with quantum microscopic domain). The first regime exhibit differences because of the new families of solutions (like runaway solutions compared to 1/r Newtonian potential solution). Whereas in the short distances regime the terms with higher derivative again start to lead and dominate over terms with lower number of derivatives because generally this regime is identified with high energies and then the more derivatives we have in the action or EOM, the higher power of energy or momentum we have in the corresponding solutions (compare this to the discussion of the scaling dimension in [21]). In ordinary field theory it is possible to conceive modifications which are true small perturbations (like a non-derivative interaction in renormalizable scalar field models), however, in gravitational setup we are doomed to consider only higher derivative modification of the Einstein-Hilbert plus cosmological constant action. Such deformations of the standard gravitational theory cannot be considered as perturbative since it is difficult to find a regime in which they are not the dominant ones over the terms with lower number of derivatives (compare also discussion in [48,49]). This remark applies not only to the jump from a 2 to a 4 but also from a 4 to a 6 or from a 6 to higher orders in the expansion. Therefore the question arises whether we should trust more the results obtained in a higher truncation based on a 6 than on a 4 and whether the results and conclusions there will not be washed away by consideration of the even more accurate model based on a 8 coefficient of the expansion and so on.
One of the possible solutions to this problem is naturally given by the spectral action approach since there are two ingredients which could help us. First one is the presence of the arbitrary energy scale Λ. Thanks to this, we can treat the terms in the expansion of the spectral action as terms in an asymptotic series in Λ −1 variable. Then despite that numerical coefficients of higher derivative terms are finite (not infinitesimal!) numbers we can make them perturbative by considering Λ very big compared to other energy scales present in the system (for example comparing to electroweak symmetry breaking scale in the Standard Model E ≈ 216 GeV). Strictly speaking the impact of higher derivatives is perturbative only when the scale Λ is sent to infinity. Another source for justification of the perturbative treatment comes with the coefficients µ 2 and µ 3 . They depend on the precise form of the cut-off function as described in Section 2. However, from physical requirements of having a good decoupling of high energy modes, the cut-off profile should be very close to flat near zero. This means that the µ 3 coefficient should be very small. And this provides an additional suppression of the higher derivative terms and allows to treat their impact on classical exact solutions as small. However, one can see that assumption µ 3 µ 2 1 blows up the expression for the beta function in (4.43). Then to render it finite one must enter into a game of playing with three parameters Λ, µ 2 and µ 3 , which is significantly more complicated and will be discussed elsewhere.
Actually the problems with dependence on the level of expansion are much deeper on the quantum level. To have a renormalizable model of quantum gravity one must treat higher derivatives as the leading and dominant terms, not as perturbatively small additions to perturbatively non-renormalizable Einstein-Hilbert gravitational action. When one does this, one indeed finds that the model based on the expansion up to a 6 coefficient is renormalizable. (The model with a 4 coefficient is formally non-renormalizable because it does not contain in the action the term with R 2 but the stronger reason is the presence of conformal anomaly in this model [45].) Actually, using the definitions in [32,34] this is a three-loop superrenormalizable model of QG, meaning that the last divergences are met on the level of threeloop computation, while from the forth loop and upwards the theory is completely UV-finite. Similarly, when we discuss the form of one-loop beta functions (related to perturbative UVdivergences) we assume that the terms giving rise to the UV behavior of the propagator are from the terms in the action with the highest number of derivatives. Not assuming this non-perturbative character of higher derivative terms would immediately spoil superrenormalizability and renormalizability of the model 9 . The UV behavior of the propagator for gravitational perturbations is the crucial thing for the discussion of any UV properties of the theory. For any local higher derivative theory the procedure of finding the UV behavior of the propagator consists of looking for the terms in the action, which are quadratic in curvatures and with the highest, but finite, number of derivatives on the metric tensor. These terms shape the ultra-violet form of the kinetic operator for quantum fluctuations. One understands that the beta functions in the model based on a 6 are different from the ones in the model based on a 4 and there does not exist any limit which makes the two match, which is obvious from the explicit formulas for beta functions in [43]. Therefore one cannot treat the six-derivative terms in a 6 as quantum perturbations in any sense. Moreover, for the derivation of beta functions of the theory the terms with six derivatives are the leading ones in the UV and hence cannot be considered small in this regime.
Analogous problem we meet when we search for the perturbative spectrum of fluctuations. 9 If one includes the effects of higher derivative terms only as vertices of the theory, while keeps at the same time propagator derived from terms with less derivatives, then new perturbative UV-divergences pop out. These divergences contain more derivatives, more even than there are in originally added higher derivative terms. Hence such theory is perturbatively non-renormalizable.
For definiteness we can study this spectrum around flat spacetime background. In order to find poles of the propagator in respective sectors of spin-2 and spin-0 fluctuations, one needs all set of terms which are quadratic in curvatures and the term linear in Ricci scalar and cosmological constant term. Inevitably in higher derivative theories we are faced with the problem of perturbative ghosts in the spectrum. These virtual states have negative sign of the kinetic term (for tachyons they have negative real part of the mass square parameters), hence they endanger perturbative unitarity of the theory. For example, the optical theorem for on-shell scattering amplitudes does not hold anymore. These are undesirable states and they should be eradicated from the theory by all means. For removing them (or their effects on observable predictions of HD theories) various approaches have been introduced: Lee-Wick prescription [50,51], fakeons [52][53][54], disappearance of unstable perturbations on nontrivial backgrounds [55]. However, none of the proposals seems to be completely satisfactory.
Of course, one can always argue that the full spectral action will give rise to also nonperturbatively unitary higher derivative quantum field theory of gravitational interactions and blame the apparent non-unitarity as the result of truncation of the spectral action to some finite-order HD models, but to make this statement precise much more of a very nontrivial analysis should be done. In our case to find poles of the propagator (or equivalently zeros of the kinetic operator governing dynamics of quantum perturbations), one needs to know all the coefficients in front of the terms quadratic and linear in curvatures. This is in distinction to the computation of UV-divergences where we needed only the coefficients of few terms with the highest number of derivatives in UV (in d = 4 we need coefficients of terms with highest number of derivatives and the ones with two and four less derivatives only). The reason for this is that the beta functions are the UV issue while the spectrum is the problem at all energy scales. The necessary information is given in (4.41) and the value of the cosmological constant coupling, c −2 = Λ 4 µ 0 . Once again we do not have any contribution from terms which are cubic in curvatures (this is true for flat spacetime propagator). Another problem is that for the theory with cosmological constant c −2 = 0 flat spacetime is not an on-shell background. (It does not satisfy vacuum gravitational EOM with the cosmological constant term and none energymomentum source of matter origin on the RHS of gravitational EOM.) Then we cannot consider quantum dynamics of fluctuations in the WKB approximation and the analysis of the propagator around flat background is purely academic. However, mathematically, as a demonstration, one can neglect this obstacle and continue with the analysis. The zeros of the kinetic operator are zeros of the respective polynomials in k 2 variable in momentum space, in two gauge-invariant sectors of spin-2 (related to the terms quadratic in Weyl tensor) and spin-0 (related to the terms quadratic in Ricci scalar) fluctuations. These zeros describe the mass square parameters of the modes. In HD models we always meet ghosts [33] as the consequence of UV-improved behavior of the theory compared to two-derivative theories. In our case, the theory is based on the action functional given in (3.26) and we have that both polynomials are of the third order in k 2 variable. This means that in each sector we expect three (possibly some are equal), in general, complex roots describing mass square parameters. For the precise values we need to solve cubic equations. We will not do this here, but we will comment on the general features of these solutions. The exact values depend on the numerical values of µ k parameters (for k = 0, 1, 2, 3) and on c's in (4.41). There are two possibilities: the three roots come in a form of one complex pair (of two complex conjugate roots) and one real root or all three roots are real. The former case is well known and then the pair is called a pair of Lee-Wick particles. They have quite peculiar properties similar a bit to a couple of unstable particles in standard field theories [56,57]. Therefore the model with six derivatives may realize the scenario of Lee-Wick quantum gravity (it was impossible to have a pair of complex ghosts in four-derivative theories).
Here, one can also ask the question how stable is the position of poles of the propagator against inclusion of higher terms in the expansion of the spectral action. Again the situation is quite delicate but not as dramatic as for beta functions (where we had discontinuous jumps when we increased the order of the expansion). Because the higher degree polynomials have more solutions on the complex plane and the coefficients of the terms with the highest power exponent on k 2 variable are highly suppressed by the scale Λ, the new roots always come in pairs from the point at complex infinity and the picture (or position on the complex plane) of the other roots is only slightly modified. This pair of new zeros moves smoothly when the value of the Λ parameter is changed from infinity, so the change in the set of zeros is continuous. For example, if we find that the theory based on the spectral action up to the coefficient a 6 is a model of Lee-Wick quantum gravity (for some definite values of Λ, µ 0 , µ 1 , µ 2 , µ 3 ), then it is likely that this feature of having addditional particles beside the real graviton only in complex conjugate pairs, will be preserved for higher orders in truncation of the spectral action, provided also that the value of the Λ parameter is large (then the number of these LW pairs will increase). Hence the Lee-Wick characteristics of models of quantum gravity is quite stable.
One might ask many reasonable questions about what would happen if we had at our disposal the full re-summed spectral action, i.e. if we would have a control over the nonperturbative form of the spectral action. Some of these questions are: What would be the exact classical solutions? Will the quantum theory be eventually unitary, renormalizable, or even UV-finite? To what extent solutions or beta functions based on truncated action reflect the situation in the full theory? Can they be treated as subsequent approximations in some perturbative scheme? We do not have even tentative answers to these important questions and due to technical reasons we must deal with the expansion of the spectral action in number of derivatives. One should expect some very non-trivial UV properties of the full spectral action [58]. It is plausible to think that the full theory may take a form of some non-local model of QG as discussed in [59]. And then the expansion that we are performing parallels the limiting method of approaching non-local models by some higher derivative models. Therefore, with such a perspective the questions of exact solutions and of beta functions acquire completely new answers in full re-summed models. For example for beta functions, we must not look into ratios like µ n−1 /µ n (cf. (4.43)), but into the limit of these ratios when n is sent to infinity. Then this changes the philosophy and we must instead ask questions about convergence radius of the analytic function given by a formal series ∞ n=1 µ n z n . Even if we know that formally the term with the highest number of derivatives in such expansion does not exist (firstly, because it is formally with n = ∞, secondly because of its coefficient vanishing as lim n→∞ µ n = 0), we can still in some sense talk about the beta function in non-local theory which is defined by the convergence radius above. Perhaps, in a similar sense we can talk and define non-perturbative beta functions of couplings in full spectral action. It remains to be seen what is the full analytic structure of the theory based on the full re-summed spectral action and whether this can be mapped to some non-local models of quantum gravity.

A Conventions and useful formulas
Here we collect some notations and standard formulas used in the main text.
Our conventions for symmetrization and anti-symmetrization of indices are the following: (µ 1 · · · µ n ) and [µ 1 · · · µ n ] mean respectively complete symmetrization and anti-symmetrization with respect to the indices µ 1 to µ n , with the proper symmetry factor (that is 1/n!). In the situation, where there are various operations nested on the same group of indices, the bracket [[· · · ]] will be used for not confusing which pair of the indices is being anti-symmetrized in the second turn (see the appendix D).
Covariant derivatives on the geometric objects (trivial from the point of view of bundle space structure) we denote either in the standard semicolon (;) postfix (GR) notation or in a operatorial prefix notation with the symbols of nabla (∇), which is however a more frequent choice in field theory.
In the Euclidean signature (used mostly through the text of the article) we choose signature of the metric tensor to be all pluses, and in Minkowski (analytically continued) case we take the time as the first coordinate and choose the signature of the metric to be (+, −, −, −) in four spacetime dimensions.
Our convention for overall signs of Riemann tensor, Ricci tensor and Ricci scalar takes respectively the following forms: Note the non-standard sign in the definition of the Riemann and Ricci tensors (opposite, for example, to conventions of Landau-Lifshitz [60]). This choice is made to agree with the notations used in the literature on the heat kernel expansion [36,37] and on non-commutative geometry [1][2][3]. Using this definition, the commutator of the covariant derivatives acting on a general tensor T α 1 ...αm β 1 ...βn (with m contravariant indices and n covariant ones) can be written as, We remind that Bianchi identities for the Riemann tensor are expressed as with the names of respectively the first and the second identity. Contracting the second Bianchi identity, we get the singly contracted second Bianchi identity Contracting one more time, will lead us to the doubly contracted second Bianchi identity: The standard expression for the Weyl (conformal) tensor in d dimensions takes the following form, One can easily find the following useful expressions: and with one power of the covariant box (covariant d'Alembertian) operator = g µν ∇ µ ∇ ν , inserted: Actually, the above formula is valid for any power (or even an analytic function) of the operator since it is a spectator in the derivation. One can notice very big similarity in the structure and coefficients of the corresponding terms between formulas (A.9), (A.10) and (A.11). This is not an accidental coincidence and is due to the complete tracelessness property of the Weyl tensor in any dimension. The match would be perfect, if we used the Landau-Lifshitz convention for the overall sign of the Ricci tensor (opposite to the one accepted in (A.2)).
The Gauss-Bonnet scalar is defined by while its generalization containing 2N + 4 derivatives is given by In the main text we use the generalized Gauss-Bonet term with N = 1: (A.14) While for N = 0, GB 0 = GB is the standard Gauss-Bonnet term, which is topological in the 4-dimensional case (and related there to the Euler invariant), for N 1 it is not topological anymore but it can be transformed to the form where R 3 stands for different cubic invariants in curvature (see the formula (3.21) and the footnote 6) and K µ is a vector field constructed from curvatures and their covariant derivatives, so the last term above is a total derivative. We remark that since GB is topological in d = 4 it does not contribute to classical EOM (obtained by the first variation of the action). Away from the case of d = 4 or for N 1 the generalized Gauss-Bonnet term in the action has an impact on EOM.

B Lichnerowicz formula
Here we derive the formula (2.6) both, to make the presentation more accessible and selfcontained, and to introduce our notations and conventions regarding spinors and Dirac operator.
We define the algebra of the flat gamma matrices with the minus sign: where δ ab is the metric of the flat d-dimensional Euclidean space in Cartesian coordinates. It is well known that Σ ab := 1 2 γ ab , where γ ab := 1 2 [γ a , γ b ], are the generators of the Euclidean version of the Lorentz symmetry, i.e. of the orthogonal group SO(d) satisfying the following commutation relations This choice of the Clifford algebra (B.15) forces us to use the following as the action for the standard Dirac spinor (two-derivative theory) on flat Euclidean space background: Introducing the tetrads associated with the metric g µν of the curved space and they satisfy the anti-commutation relation In passing, we note that the small Greek letters we use for curved (world) space indices, while the small Latin letters we use exclusively for denoting flat (tangent) space indices. In the former space we use the curved metric g µν to raise world indices, while in the latter flat space we use the Kronecker delta tensor δ ab to do the corresponding operation on flat indices. Using these notations the standard Dirac operator is given by where ω µ = 1 4 ω ab µ γ ab is determined by the requirement that the tetrads are covariantly constant with respect to the covariant derivative ∇ ω µ ∇ ω µ e a ν = ∂ µ e a ν + δ bc ω ab µ e c ν − Γ ρ µν e a ρ = 0 .

(B.22)
Here ∇ ω µ = ∇ µ − ω µ is the total covariant derivative, while ∇ µ is the usual, Levi-Civita, one. 10 Since the tetrad (vielbein) is valued both in the tangent as well as curved space (it possesses both types of indices), then the total covariant derivative ∇ ω µ must include connection coefficients from both spaces. In the external space these are given by standard (metric) Christoffel coefficients (and then ∇ µ e a ν = ∂ µ e a ν − Γ ρ µν e a ρ ), while in the tangent space this role is played by the SO(d) gauge potentials ω ab µ . The main result that allows the direct application of the heat kernel techniques is the Lichnerowicz formula where R is the scalar curvature of the metric g µν . Due to the importance of this formula let us sketch its proof. We will need several identities: • R µνρσ γ µν γ ρσ ≡ R µνρσ γ µ γ ν γ ρ γ σ = 2R. This is easily proven noticing that γ ν γ ρ γ σ = γ [ν γ ρ γ σ] − g νρ γ σ − g ρσ γ ν + g νσ γ ρ and using the first Bianchi identity, R µ[νρσ] = 0.
Using these identities it is straightforward to calculate D 2 : 10 This is the same definition as in (B.21), taking into account that acting on a spinor (being in a representation not carrying any Lorentz indices, so not on a gravitino) ∇ µ is just a partial derivative ∂ µ . 11 The choice of a sign in (B.24) agrees with the convention for the sign of the Riemann tensor stipulated in (A.1).

C Calculation of the trace
The spectral action is a special case of the following more general expression: Tr χ(P ) , (C. 26) where χ is some "more or less" nice function (the exact meaning of this will be given below) and P is some positively-definite operator on a Hilbert space. In our case, P = −D 2 is represented on the Hilbert space of square-integrable spinors. Because the result for this expression in terms of the heat kernel expansion is one of the main tools in our approach, here we give a detailed (and more or less rigorous) derivation of this expansion. Also, this will hopefully make the paper self-contained. Let us start by requiring for χ = χ(p) to be a piecewise continuous function on R + such that In other words, we require that the small-p asymptotics (near p = 0) is given by b 1 p a 1 (χ(p) = O(p a 1 ) for p → 0) and similarly the large-p asymptotics (in the p → +∞ limit) is given by b 2 p a 2 (χ(p) = O(p a 2 ) for p → +∞). We also demand that a 2 < a 1 . The interval (−a 1 , −a 2 ) is called the fundamental strip of χ. E.g., if χ(p) is some smooth cut-off function, which at infinity goes to zero faster than any negative degree monomial of p and is of order of p 0 when p → 0, then the fundamental strip is (0, +∞). Also, let the integral φ(s) = ∞ 0 p s χ(p) dp p (C.28) be convergent when s belongs to the fundamental strip (the function φ(s) is called a Mellin transform of the function χ(p)), then χ(p) can be recovered via the inverse Mellin transform: where c is a number which should also belong to the fundamental strip. Let P be a positively-definite operator and χ(p) be some cut-off function (i.e. its fundamental strip is (0, +∞)) with φ(s) being its Mellin transform. Then, using the spectral functional calculus, we can define a function of an operator P χ(P ) = 1 2πi The integration contour in (C.31) can be chosen in such a way that it encircles all the poles of the integrand, that is all the poles are inside the contour and the contour is closed at infinity. It is possible to show that the contribution of the integration over this part of the contour is zero. We will find the poles of ζ P (s)φ(s) by some indirect method -using the known results for the heat kernel expansion. The relevance of the heat kernel will become clear after re-writing the zeta-function (C.32) in terms of Tr e −tP , (C. 33) which is the object called a trace of heat kernel. Towards this end, let us use Cahen-Mellin integral (which we will also need later) and its inverse: Now, let us analyze the poles of the underintegral expression in (C.31), ζ P (s)φ(s). 12 Though the convergence of the integral (C.35) is guaranteed if Re(s) > 0 now one should be careful because taking trace over the infinite-dimensional space may introduce new divergences, see below.
First, let us show that φ(s) has poles at s = 0, −1, −2, ... . We know that φ(s) is regular when Re(s) > 0 (the fundamental strip). Now considerφ(s) defined bỹ and consider the following integral where we again assume that the contour could be closed to encircle all of poles. Then using the Cauchy residue theorem 13 we havẽ Hence,φ(s) ≡ φ(s), which proves the above statement about the poles. The last equality in the above formula is valid within the analytic convergence region of the Maclaurin series of the cut-off function χ(p). In what follows, we will assume that the convergence radius is infinite. Now, let us consider the poles of ζ P (s). For this, we will use the following asymptotic small t expansion for the traced heat kernel coefficients: Tr e −tP n 0 t n−d m a n (P ) , (C. 41) where d is the dimension of the manifold M, m is the order of P and a n (P ) are defined by a n (P ) = M a n (x, P ) √ g d d x (C. 42) for some known Seeley-DeWitt coefficients a n (x, P ). By the inverse Mellin transform (C.29), we can write for small t Once again being sloppy about the contour and not pretending to be rigorous we can read of poles of Γ(s)ζ P (s): Res(Γ(s)ζ P (s))| s= d−n m = a n (P ) . (C.44) 13 Here by χ (k) we denote in a standard way the k-th derivative of the cut-off function with respect to its argument p: χ (k) = d k χ(p) dp k .

D General equations of motion
The simplest way to calculate the variation of the Riemann tensor and all of its contractions is to vary directly the defining formulas (A.1), (A.2) and (A.3), taking into account that the variation of the Christoffel symbols Γ σ νρ are tensors given by where h µν := δg µν . Then one trivially gets the linear part of the variation: To vary the Weyl-dominated form of the action (3.24), one also needs the variation of the Weyl tensor in d = 4 space (or spacetime) dimensions: Using these variations (and the related ones) the general form of the tensor of classical equations of motion, E αβ , resulting from the action in (3.23), in the Riemann basis is given by (to be symmetrized with respect to (α, β) pair of indices, if needed): This result is independent of the signature since it is written in a generally covariant way. Therefore, it holds both for Euclidean spaces and also for Lorentzian spacetimes of any dimension.