Gravity in $d=2+\epsilon$ dimensions and realizations of the diffeomorphisms group

We discuss two distinct realizations of the diffeomorphism group for metric gravity, which give rise to theories that are classically equivalent, but quantum mechanically distinct. We renormalize them in $d=2+\epsilon$ dimensions, developing a new procedure for dimensional continuation of metric theories and highlighting connections with the constructions that previously appeared in the literature. Our hope is to frame candidates ultraviolet completions of quantum gravity in $d>2$ and give some perturbative mean to assess its existence in $d=4$, but also to speculate on some potential obstructions in the continuation of such candidates to finite values of $\epsilon$. Our results suggest the presence of a conformal window in $d$ which seems to extend to values higher than four.

The regularization of a quantum theory of Einstein-Hilbert gravity leads to nonrenormalizable divergences beyond the first order in the loop expansion [1,2]. The problem arises because the perturbative coupling, Newton's constant, is dimensionful and consequently leads to a non-power-counting renormalizable expansion in . For this reason, it is generally understood that one should re-arrange the perturbative expansion in powers of some ultraviolet (UV) energy scale, which could be the Planck mass itself, and that quantum Einstein-Hilbert gravity can be seen at most as an effective theory below such scale [3].
A different perspective on the problem is achieved by noticing that Newton's constant, G, is not dimensionful for every spacetime dimension d, hinting that the theory could be perturbatively renormalizable for some value of d. In fact, not only G becomes dimensionless in d = 2 spacetime dimensions, but it also has an asymptotically free beta function, β G ∝ −G 2 , meaning that, in principle, one could obtain consistent predictions from perturbation theory, that are valid up to arbitrarily high energies [4]. We refer to d = 2 as the critical dimension of the Einstein-Hilbert action.
This notion is not particularly useful for the physically interesting d = 4 case, unless one realizes that in d = 2+ dimensions one can trade the perturbative expansion in G for an expansion in the parameter . Re-instating the canonical mass dimension of Newton's constant, one finds that its renormalization group (RG) running, in units of an RG scale µ, is − G + β G ∝ −G 2 ; as a consequence there is a scale invariant value G * ∼ O( ), arising as a fixed point solution of β G = 0. The limit → 0 has been used to reproduce some results of 2d quantum gravity [5,6], that were previously obtained by other means [7], but we anticipate from the main discussion of this paper that the precise value of the proportionality constant in β G has a prominent role in further clarifying this limit 1 .
The found fixed point G * has some interesting and useful properties, because of its UV nature [9]. A simple analysis reveals that it separates the regions G < G * and G > G * , which are driven by the RG to the infrared (IR) Gaußian and strongly-interacting limits, respectively. The UV fixed point in d = 2 + is almost as good as asymptotic freedom: in the UV the theory always remains consistent by "approaching" a nonperturbative interacting limit instead of becoming asymptotically free. In the IR, instead, the theory could flow to the Gaußian phase and be described by the aforementioned effective theory, which would see its UV completion as "hidden" by the UV energy scale at which Gaußian perturbation theory breaks. A theory with a UV fixed point and a finite number of relevant RG directions is known as asymptotically safe [10].
The idea that gravity could be asymptotically safe in d ≥ 2 and, most importantly, in d = 4 has been con-jectured a long time ago already by Weinberg [11]. The existence of the nontrivial UV fixed point O( ) guarantees that, at least for infinitesimally small , there is a UV completion. This has given a compelling reason to push forward the investigation of the asymptotic safety conjecture, which has received increasing attention over the past few decades [12].
Hystorically, after some first papers that pushed forward the status of the conjecture through the use of perturbation theory and the -expansion, e.g. [13,14], most of the literature has eventually settled on the use of a nonperturbative method, known as functional RG, to test the validity of the conjecture as pioneered in Refs. [15] and [16] by Reuter and Souma, respectively. However, the nonperturbative approach comes at a price: the Wetterich equation, which governs the functional RG flow [17], comes with a more severe scheme dependence, if compared with massless renormalization schemes such as dimensional regularization and minimal subtraction (see appendix B for some further comments on scheme dependence in dimensional regularization and the role of quadratic divergences). This scheme dependence has the disadvantage of mixing, a bit disturbingly, with both gauge-and parametrization-dependence of the path-integral [18], making unclear which are the physical predictions of the theory in terms of observables [19,20], even though the problem can somehow be relaxed by going on-shell [21,22]. The potential danger, which is faced by the practitioners of the asymptotic safety conjecture, is that not having clear gauge-independent results undermines the cumulative work towards the proof of the conjecture itself, and delays further results on other pressing issues such as, for example, unitarity [23].
Still, we, the authors, are far from having a negative stance towards asymptotic safety, even after mentioning some of the weaknesses of the modern approach. Our humble opinion on the conjecture is that the pursuit based on the functional RG should be coupled with another one that gives less scheme-dependent results, even at the price of nonperturbativity, and these two searches should come together convincingly. For this reason, we want to apply the perturbative framework, that we briefly outlined at the beginning, and resurrect the aged discussion on gravity in d = 2 + dimensions. Therefore, we follow Weinberg's original idea [10], and assume that the fixed point of a putative asymptotically free theory of gravity comes from the continuation to d = 4 of the asymptotically free theory in d = 2. The obvious limitation is that this requires an extrapolation to = 2, which is naive at best and dangerous at worst, but, fortunately, the results based on perturbation theory should be weighted with the body of work coming from nonperturbative methods, as we are going to see in the following. It is important to mention that our idea of using information from perturbation theory is far from new [24], although we seem to occasionally come to slightly different conclusions.
One immediate advantage of resorting to perturbation theory is that, even if we start from the same "classical" bare theory, we can identify two distinct realizations of the gravitational path-integral based on two isomorphic, but inequivalent, realizations of the diffeomorphism group. We refer to the two models resulting from these parametrizations as Einstein's and Unimodular-Dilaton gravity, respectively, following a nomenclature introduced in Ref. [25]. The two realizations, which we discuss in the next section, have often been interpreted as related to the parametrization dependence of the results [26]. Instead, we explicitly show that each realization does not depend on gauge and other parameters on-shell. Several ideas have surfaced when revisiting older results with modern eyes, and we clarify them along the way. The rest of the paper is structured as follows: in Sect. II we discuss the two isomorphic realizations of the diffeomorphisms group and the relative actions; in Sect. III we explain a three steps procedure on how to dimensionally continue the theories for the applying regularization methods; in Sects. IV and V we show the leading renormalization of the two different actions; in Sect. VI we expand our discussion on the role of the conformal mode and the two dimensional limit; in Sect. VII we attempt a conclusion and speculate on the results; in appendix B we briefly touch the topic of quadratic divergences in dimensional regularization.
II. Diff VS Diff * As starting point of this discussion we use the familiar (Euclidean) Einstein-Hilbert action for Einstein gravity (EG) in d-dimensions in which we introduced the curvature scalar R and two couplings, g 0 and g 1 , containing the traditional cosmological and Newton's constants. The strength of the gravitational interaction is weighted by G, which we choose to be g 1 = G −1 , avoiding any unnecessary normalization. Due to the covariant nature of all the elements appearing in (1), the action is manifestly invariant under diffeomorphisms. Infinitesimal diffeomorphisms act as changes of coordinates and are generated by a vector field, x µ → x µ + ξ µ (x), resulting in a transformation of the metric Since the composition of two transformations is still a transformation, the diffeomorphisms form a group which we denote Diff. The algebra is obviously closed where on the right hand side the commutator denotes the standard Lie brackets of two vector fields.
Following Ref. [25], we now define a slightly more general action which goes under the name of dilaton gravity (DG) action. We begin by parametrizing the metric as a dilatonic factor times another metric, g µν = Ω 2g µν . This results in the dilaton action, S D [Ω,g] ≡ S E [g], which has an additional local Weyl symmetry caused by the possibility of rescaling the two factors Ω 2 andg µν while leaving g µν invariant. We further parametrize Ω = ϕ 2 /d−2 to get withR = R(g, ∂g). An additional rescaling of the field ϕ → αϕ by the constant α = (d−2) /8(d−1) can be used to normalize the kinetic term of ϕ; it brings the first line to a familiar conformal invariant action with nonminimal coupling in d dimensions (coupled with the so-called Yamabe operator) [27]. For d > 2 the required rescaling is purely imaginary and is generally associated to the conformal mode instability because of the wrong sign of the kinetic term of (4), while for 1 < d < 2 there is no such problem. We defer the discussion on the case d = 2 to Sect. VI. By construction, the dilaton action (4) has two symmetries: a diffeomorphism invariance in which the two fields transform independently (indices are raised and lowered throughg µν and its inverse), and a Weyl invariance on local rescalings with respect to an arbitrary scalar function that leave the original g µν invariant. Infinitesimally, the Weyl transformation acts as for some infinitesimally small scalar function ω. The complete symmetry group of (4) is a semidirect product of the two subgroups Diff Weyl because of the nontrivial action of Diff on Weyl.
It is clear that the symmetry group of the dilaton action (4) is enhanced because of the redundancy introduced by the combination g µν = Ω 2g µν . The obvious way to get rid of both the redundant symmetry and the scalar degree of freedom Ω is to take Ω = 1, which breaks down (7) to the original Diff group and (4) to (1). However, there is an alternative way to break the symmetry group, which involves assuming that the metric theory based ong µν in (4) is unimodular. In the unimodular theory we require any possible metricg µν to obey the where v is some fixed volume d-form (oftentimes just chosen v = 1).
We refer to the unimodular version of (4) as unimodular-dilaton gravity (UDG), and its action is simply (4), with the restriction that the metricg µν is unimodular. The symmetry group of the unimodular version of (4) is the subgroup of (7) that contains all the transformations that leave a given v invariant. Infinitesimally, it is generated by the combined transformations δ * = δ ξ + δ w ω , from (5) and (6), for which δ * √g = 0, thus preventing v from any transformation. Of course, the generators ξ µ and ω are not independent: in fact, one way to think at the transformations δ * is as standard diffeomorphisms δ ξ , which are then compensated by a Weyl transformation with ω = (d−2) /2d∇ · ξ restoring the original volume element. Therefore, δ * only depends on the vector ξ µ , hence we adopt the notation δ * → δ * ξ . The action of δ * ξ is thus Obviously, the unimodular transformations are a subgroup of (7), which can be though of as a symmetrybreaking pattern coming from the quotient of (7) with the equivalence relation δ * √g = 0. We denote it as Diff * : Diff Weyl Importantly, the new group is isomorphic to Diff itself, Diff * Diff, using the relation which can be proven with a bit of work. This could be foreshadowed by the one-to-one correspondence between transformations and vector fields ξ µ , together with the fact that, modulo normalizations, two vectors combine into a third one only through Lie brackets. 2 Of course, in this discussion we are identifying the groups with their algebras, so our statements hold, realistically, only at the infinitesimal level and the groups might differ globally.
Having identified the larger symmetry group in the breaking pattern (9) leading to Diff * , we find convenient to write a more general action for the UDG theory, in which the condition of Weyl invariance is relaxed 2 The result is less trivial than it looks, though, because a different breaking pattern can lead to yet another symmetry group [25]. One such possibility is that we are left with a semidirect product of unimodular diffeomorphisms and Weyl transformations. in which we introduced a function, V (ϕ), that plays the role of cosmological constant and a topological charge q that couples ϕ to the curvature scalar in a way familiar to Liouville field theory [28].
There are important reasons behind the chosen parametrization. The charge q breaks classical Weyl invariance because it gives a nonzero value to the trace of the stress energy tensor of ϕ overg µν , that is, T = T µ µ = −2q 2 R. The value q can, however, be chosen so that the quantum Weyl anomaly is zero T = −2q 2 R + T rad = 0, where the second term comes from radiative corrections [29]. This can be done consistently order by order in perturbation theory, as shown in Ref. [30]. The quantum theory is thus invariant under Weyl transformations only if q is chosen appropriately, and the function playing the operator V (ϕ), becomes V (ϕ) = g 0 ϕ 2d /d−2 . The reason why it is important to restore Weyl invariance is that the action (11), in the words of Ref. [31], is invariant under "volume preserving diffeomorphisms", but at a Weyl-invariant critical point these are promoted to "full diffeomorphisms." In our discussion, we identify the volume preserving transformations as unimodular transformations, and the full diffeomorphisms group as Diff * , which is correctly isomorphic to Diff, though not exactly the same, as we previously argued following Ref. [25]. Conformal invariance should emerge as a symmetry of the (scale invariant) RG fixed point. Notice also the difference of this approach with the traditional unimodular gravity realization without the dilaton, in which the cosmological constant emerges as an integration constant of the equations of motion [22].
A more general analysis could be implemented by promoting the operators related to the curvature scalar to function too, as suggested in Ref. [32], which would contain both g 1 and the topological charge. In this case the breaking would be entirely realized at the fixed point starting from a dilatonic theory, but we leave this possibility open for now. A simple visual summary of the breaking patterns relating the three actions is given in Fig. 1, which depicts an upper wedge of a more general pattern described in Ref. [25].
For the rest of the paper we focus on the Einstein and unimodular-dilaton realizations, with actions (1) and (11), respectively. The symmetry groups of the two theories are isomorphic, so the natural question is whether they actually are the same physical model or not. From the Hamiltonian analysis emerges that the two models propagate the same degrees of freedom [25], so it is fair to say that classically the Einstein's and unimodular-dilaton gravitational models are equivalent. What is less clear is if the two models are quantum mechanically equivalent.One way to investigate this difference is to inspect the beta function of Newton's constant in each realization. From a general scaling analysis, see for example Ref. [6], we know that the beta function must be of the form in which we denoted with c the central charge, which is a number of particular importance in the context of two-dimensional conformal field theory (CFT) because it counts the effective number of degrees of freedom that are integrated in the path-integral. In general, c will be the sum of two contributions, one coming from the matter fields and the other coming from the gravitational degrees of freedom.
As discussed in Refs. [4,24], the value of the central charge for two dimensional quantum gravity is scheme dependent. In particular, coupling the theory with matter allows for the introduction of operators of different mass dimensions and the prescriptions to go on-shell increase since it is possible to solve the equations of motion with respect to different operators. The coefficient of the beta function for Newton's constant is then shifted by a value proportional to the mass dimensions of the operator solved by the equations of motion [24].
In absence of matter,though, possibility are much more restricted and the allowed schemes are determined by the realization of the symmetry chosen for the gravitational sector solely. We have that in d = 2 Einstein's (1) leads to c = −19 [13], while unimodular-dilaton's (11) might also lead to c = −25 [14]. Evidently, in both cases c < 0, implying that a fixed point g * ∼ O( ) exists, but it is undeniable that the quantum mechanical effects are different. Notice that, in the past, this difference was often attributed to a parametrization dependence (exponential vs linear background split of the metric), and more often than not it was believed that the second result is the "correct" one because it reproduces results from string theory and Liouville gravity in the limit → 0 [4,33]. By generalizing results for both theories, in the following sections we show that the difference has nothing to do with parametrizations, rather the symmetry plays a much bigger role. A notable exception of a paper, that partly discusses the two models as two similar but independent realizations, each with its own central chargeis Ref. [23].

III. THREE STEPS CONTINUATION IN d
The analytic continuation to d = 4 of a d = 2 + computation requires the extrapolation to = 2, which is obviously a big limitation of the perturbative approach that works well at small . Nevertheless, the two dimensional limit, → 0, can be used to give theoretical arguments against or in favor of either realization of Diff described in the previous section.
The strategy that we adopt for dealing with the regularization of either of the theories discussed in Sect. II is dimensional regularization with modified minimal subtraction (MS) of the divergences close to d = 2, implying that we subtract the poles 1 /d−2 and a small finite part after analytic continuation of the results in the dimensionality. There are several difficulties that arises with a naive application of MS in gravity. The first and most prominent one is that several tensor contraction also have d-dependent outcomes, such as g µ µ = d, which might change the finite part of the subtractions when multiplying a pole or, in the worst case, entirely remove a divergence. This obviously could make ambiguous the status of some divergences.
In order to lift any ambiguity we define three fundamental steps for dimensional continuation and regularization of the metric theories, while keeping in mind that the desired outcome is extrapolation above d = 2. The procedure will unfortunately imply a proliferation oflike symbols.
• The first step is that we actually continue any Feynman diagram 3 in the dimension as d = 2 → d = 2 − ζ; we stress that, for the moment, ζ = − which was previously introduced. This means that, whenever an integration measure appears, we promote it as d 2 q → d 2−ζ q. For ζ > 0 the diagrams that are relevant for perturbation theory converge (generally for Re ζ > 0), but, equally importantly, the conformal mode of the metric is stable (see also the discussion in Sect. VI; in short it means that the constant in front of the kinetic term of ϕ in (4) is positive). We elect the regime ζ > 0 as the one in which we make computations of the radiative corrections; divergences thus appear as poles 1 /ζ that must be subtracted with opportune counterterm operators, and their coefficients assemble into beta functions of renormalized couplings.
• The second step is that, any time a tensorial contraction returns the dimension of spacetime, we simply denote such dimension with d and treat it as a parameter. The simplest example is of course g µ µ = δ µ µ = d. In particular, and this is very important, we do not substitute d = 2 nor d = 2 − ζ or d = 2 + when computing divergences, though, either limit can be taken later if needed. This has the advantage that d appears parametrically in our computations, much like N appears when computing the renormalization of an SU (N ) gauge theory. Similarly to gauge theories, by setting d = 2 or d = 4 the metric fluctuations have the expected degrees of freedom in a given dimension.
• The third step is the one of continuing the results to d > 2. This is achieved by continuing ζ = − with > 0, that is the forbidden region ζ < 0, which explains why we keep it separate from the process of dimensionally regularizing the theory. What we imply in this last step is that only after having regulated and renormalized the model, and obtained a beta function such as (12), we then continue to d = 2 + > 2 dimensions by noticing that the coupling G must have negative mass dimension. This step introduces the dimensionless coupling through the replacement G → Gµ − , where µ is the RG scale, effectively measuring the coupling constant in units of µ. The net effect of this last step is that the beta function β G acquires the scaling term G as explained in Sect. I.
In principle, the three above steps leave us with two independent parameters related to the meaning of dimension d, which can or cannot be identified, according to necessity. A summary of our general strategy is the following: we eliminate poles in ζ coming from diagrams entirely though MS-like subtraction (first step), so we can express the beta functions as d-dependent objects (second step), that we continue to d > 2 (third step). The final result in an -expansion series that has d-dependent coefficients. With these steps in mind, it should be clear that one can investigate the two dimensional limit by taking d = 2 + and → 0, but can also estimate the four dimensional limit by taking d = 4 and extrapolating to → 2. This procedure is clarified by the explicit examples in the next two sections.
The advantage of our procedure is that it breaks down the problematic dimensional continuation of a gravitational theory in manageable steps, which have a range of validity that is under control and can be discussed separately. For example, the first step ensures that the conformal mode is actually stable and makes the integrals converge, thus picking up a specific vacuum [34], which is different from, say, the one in Refs. [35][36][37][38]. The second step mostly ensures that we have the expected number of propagating degrees of freedom in a given dimension. Our opinion is that the potentially dangerous step is the third one, because it involves the continuation above d = 2 by means of the -expansion in a regime in which the theory does not converge, strictly speaking, but we defer a more thoughtful discussion of this point for the speculative part of the conclusions.

IV. EINSTEIN'S ACTION
We begin by considering the Einstein-Hilbert action (1). The path-integral can be constructed using the background field method by splitting the metric as in which g µν is an arbitrary background, h µν are the fluctuations to integrate over, and λ is an arbitrary parameter to test the parametric dependence of our results. In principle, many more parameters are hidden in the O(h 3 ) terms, but we choose to highlight only λ to explicitly test the parametric dependence. 4 The split allows us to preserve manifest covariance under the background version of (2), in other words δ ξ g µν = L ξ g µν (while h µν transforms as a standard symmetric tensor), which is sometimes known as background symmetry. However, the background symmetry is not the one that we have to gauge-fix, since, in fact, we want to preserve it. The gauge-fixing must fix (2) seen as a transformation of h µν at fixed arbitrary g µν , which is nonlinear because of the nonlinearity of the right hand side of (13). We can reconstruct the correct transformation on h µν orderby-order in h µν itself by inverting the following relation (parentheses indicate symmetrization of the indices): Using (13) on metrics and connections on the left hand side, we find which is often called full quantum symmetry. Indices can now be raised and lowered by the background metric. Formula (15) is given to the order that is necessary for our computation, which does not happen to contain λ. However, higher loop computations need additional orders in h µν , so we show two additional orders and their λ-dependence in App. C.
As gauge-fixing we choose a familiar de Donder form which includes two gauge fixing parameters α and β. Notice that this gauge fixing only changes the fluctuation's two point function of the theory, while higher vertices are unaffected by the gauge fixing. Alternatively one could use the metric g µν instead of the background in the construction of (16) at the price of modifying further vertices and making the ghosts action more complicated. Speaking of which, the ghost action for ghost c µ and antighost c µ fields is computed straightforwardly A common choice for the gauge fixing parameters is α = β = 1, because it allows to write the fluctuation's Hessian as an operator of Laplace-type and to express everything in terms of simple heat kernel formulas. Computations with arbitrary α and β become difficult very quickly. However, we still want to test gauge independence of our results, so in this section we make the simple choice in which δβ in an infinitesimal parameter. The advantage of this choice is that we still have a simple Hessian because we can expand in powers of δβ and compute (for example) the first order correction in the gauge dependence of our computations. The disadvantage is that the first order only proves that we are in a saddle point of the gauge dependence, and does not prove full gauge independence. We thus assume that the cancellation of linear contributions in δβ is an evidence for gauge independence. In practice, the linear contributions in δβ emerge as new vertices with two external legs in the fluctuations and in the ghosts, coming respectively by inserting (18) in (16) and (17), which dress the two point functions.
For the computation we adopt the formalism of the effective action, which includes the one loop radiative corrections coming from the Hessian operators, O, of metric fluctuations and ghosts. For our needs it is sufficient to compute Γ, which is naively a separate function of g µν and h µν , at h µν = 0. Since we use dimensional regularization and have no explicit cutoff, we assume that the split symmetry introduced by (13) is not broken and that the corresponding Ward identity can be used to reconstruct immediately the full Γ. Divergences are computed covariantly with heat kernel methods [32,40,41], so that the theory can be renormalized as a standard field theory in curved space [42]. We follow the three steps procedure outlines in Sect. III, the result is expressed in terms of the background in which R = R[g] and it is important to remember that it was computed at h µν = 0. The regularization to d = 2 − ζ of momentum space (heat kernel) integrals introduces the scale µ, that should be regarded as a momentum scale in reference to the background g µν and that balances the change in the dimensionality of d 2 x → d d x and which we use for the RG. As one might expect, the divergence proportional to the volume is also proportional to g 0 , because g 0 is the only parameter in the theory with the correct dimensionality, but also to (g 1 ) −1 , because the perturbative coupling of the theory is Newton's constant G = (g 1 ) −1 . Less expected could be the gauge and parametric dependence that is displayed through δβ and λ, respectively. In principle, following Ref. [41], one could eliminate (20) by redefining the bare parameters g 0 and g 1 of (1), resulting in the coefficients of the poles becoming RG beta and gamma functions of the two couplings. However, both functions would be gauge and parametric dependent, implying that they cannot be associated to a physical observable of the theory. The reason why this happens is that, out of the two monomials appearing in (1), only one is independent on-shell, and therefore the information of only one physical beta function can be extracted. To put (20) on-shell we choose the following parametrization for some scalar constant A, a symmetric tensor J µν coupled to the equations of motion of (1), and the background's Einstein's tensor G µν = R µν − 1 /2R g µν . Onshell the equations of motion of the bare action vanish at this order, and therefore any dependence on J µν decouples [24], while A remains and should be subtracted with a counterterm of Newton's constant [4,43]. We stress that with this choice of parametrization the volume operator has, by construction, a trivial scaling. 5 Comparing (20) and (21), we find The important result is that the coefficient A is gauge and λ independent, so we can subtract the associated divergence and expect a physically sound interpretation, while the tensor J µν carries the unwanted dependencies. 6 Subtracting the first term of (21) through a countertem for g 1 leads to the beta function for Newton's constant β G = −AG 2 . According to the third step discussed in Sect. III we now also continue the theory above two dimensions replacing G → Gµ − which results in a scaling term for β G . Finally, we have while, according to our choice of parametrization in (21), g 0 is not renormalized, so the volume operator scales trivially [24]. 7 Eq. (23) concides with the beta function given in [46], which was obtained directly in d > 2 with a cutoff, therefore we confirm that it is a universal result. Some additional comments are in order. An interesting aspect of (23) is that it gives perturbative information in G, which results in a perturbative information in = d−2 at the RG fixed point. However, upon identification of the two parameters, the dimensionality d appears also parametrically through the coefficient A (recall that we regard this d dependence as the N dependence of an SU (N ) gauge theory). The limit d → 2 of (23) agrees with Refs. [16,41] A scaling relation suggests that β G = c /24πG 2 with c a number known as the central-charge in the context of CFT [6], suggesting c Diff = −19. Even in absence of a conformal interpretation, the value −19 is well-known in the perturbative renormalization of 2d gravity [13] and, from our point of view, it is a signature of the Einstein's realization of Diff. One interesting point to make is that if we instead decided to take the gauge-dependent counterterms (20) at λ = δβ = 0 and renormalize (20) offshell, we would still get c Diff = −19 in the same limit, but also a beta function for g 0 , as it is commonly done in the literature of asymptotic safety. This limit, which is gauge dependent because a change in β would result in a different β G , by accident gives the on-shell physical value [41]. Rather surprisingly, for d = 2, λ = 0 and with the gauge choice δβ = 0, the beta function for G obtained from the off-shell counterterms (20) is the same as the correct on-shell one given later in (24). We argue that this unforeseen accident has lead many to believe that the 6 The procedure that we describe is discussed in other works. See, for example, Refs. [24,43,44]. However, we seem to have a different result as for the gauge-dependence of Ref. [43]. 7 An equivalent choice would be to parametrize the first term of (21) without the curvature scalar, which would result in a nontrivial gamma function for g 0 and a trivial one for g 1 . This can be understood as a different choice of units. We recommend the discussion of Ref. [45] on this point. coefficient of the gauge-and λ-independent beta function (24) was actually caused by the so-called linear parametrization g µν = g µν + h µν (λ = 0 in (13)), while it is actually independent on the parametrization [24]. Another interesting aspect is to evaluate when A > 0 as a function of d, which is the condition for which β G has a nontrivial O( ) fixed point. We find that the condition for the existence of a real fixed point returns a conformal window −4.685 < d < 7.685 [46]. This leads us to some interesting speculations: on the one hand, d = 4 is well within the upper bound of this fixed point; on the other hand, there is an effective upper critical dimension for this continuation, so the theory is not asymptotically safe up to d → ∞, as suggested in Ref. [47], but rather is safe up until d c = 1 /2 3 + 3 √ 17 ≈ 7.7. This is also one of the conclusions of Ref. [39], in which it is predicted d c ≈ 5.7. Of course, the size of the conformal window is expected to change by including further orders of perturbation theory as d c = d c ( ), which must also take into account the nonlinear interplay between d and . Nevertheless, it is an important benchmark for our argument to find agreement with a sophisticate parametric analysis such as the one in Ref. [39], but also a useful validation for Ref. [39] to find a perturbative argument confirming its analysis.
Having mentioned the possibility of going to two-loops order, there are some important aspects to clarify in this regard. For the next-to-leading computation of Γ it is not sufficient to renormalize all two-loop diagrams at h µν = 0 (see Ref. [48] for a concise explanation). This happens because Diff is realized nonlinearly from the point of view of h µν , but also because the split (13), which is written in terms of bare quantities, differs from its renormalized version. In other words, the true quantum metric g µν of Einstein gravity is different from the expectation value of the bare relation (13) and requires additional renormalization conditions. To achieve a full two-loops computations, however, it is not necessary to renormalize all possible vertices with external h µν legs, either. The most economic way to go about is to renormalize a one-point function in h µν , which could be associated to a nonlinear source as done for the 2d nonlinear sigma model in Ref. [49]. This step is necessary to ensure that 1 /ζ 2 poles correctly cancel between two loop diagrams and one loop diagrams with counterterms insertions [50]. The cancellation of higher order poles is generally regarded as an important nontrivial check in dimensional regularization, which unfortunately does not have a clear equivalent in functional RG approaches to asymptotic safety.
A final comment concerns the conclusions of Ref. [13], in which it is argued that it is impossible to renormalize the theory beyond the leading order, because of the kinematic pole. The assumption made in Ref. [13] is that one should find counterterms for both the 1 /ζ poles of momentum diagrams and for the kinematic 1 /d−2 poles, but this conspires against the delicate cancellation of the 1 /ζ 2 poles. A similar argument could be made using the one loop renormalization of curvature square composite oper-ators, in which one might naively expect that kinematic poles survive (because the relevant diagrams have more propagators than vertices) and make the RG of some observables divergent in the limit d → 2. The problem with Ref. [13] is that the result is gauge dependent, and correctly reproduces c Diff = −19 only accidentally. The procedure that we have used to go on-shell clearly lifts away unwanted dependencies, such as gauge's, but also kinematic pole's as evident from J µν of (22). Preliminary results by some of us [51] show that this problem is lifted, at least for composite higher derivative operators. It would be desirable, however, to check in detail the two loops renormalizability of the two dimensional Einstein's action using the three steps prescription of Sect. III.

V. UNIMODULAR-DILATON ACTION
Now we concentrate on the unimodular-dilaton action (11). The dilaton field ϕ and the unimodular metricg µν must be varied independently. Following the discussion of Sect. II, we admit a deviation from Weyl symmetry, so, effectively, the theory could just be thought as the one of a scalar field coupled to a unimodular metric. The hope is that at the critical point Weyl invariance is restored and combines with unimodularity to produce a Diff * invariant action. Keeping this in mind, we again use the background field method expanding ϕ andg µν The unimodular metricg µν is expanded using a fluctuationĥ µν which is traceless over the unimodular backgroundĝ µν , that is,ĥ µ µ =ĝ µν h µν = 0. With this choice the volume form is preserved both by infinitesimal and finite transformations, √g = e trĥ √ g = √ĝ . We stress the comparison between (25) and (13): in the unimodular-dilaton case we are forced to adopt an exponential form (otherwise we would be breaking the Diff * symmetry of the background) and the volume element of the metricsg µν andĝ µν are forced to be the same. The first line of (25) can also be interpreted as the finite transformation relating two unimodular metrics in the same equivalence class. When integrating over the fieldĥ µν , it is important to keep in mind that it is traceless, for example when contracting the identity in its field space.
The dilaton ϕ instead is expanded about a background ϕ 0 , which is not necessarily constant, though we ultimately will specialize to a constant dilaton to obtain the beta functions. In fact, the size of ϕ 0 changes the normalization of the operator associated to g 1 in the unimodular theory,R, and this becomes important when evaluating the renormalization group scale. Adopting a redefinition of the dilaton field, it is easy to see that ϕ andg µν are the same degrees of freedom that are considered in Ref. [14], in which the configuration ϕ 0 = 1 is chosen to evaluate the counterterms. Perhaps a bit controversial is the fact that we claim that in Ref. [14] the authors renormalized two dimensional unimodular-dilaton gravity, rather than Einstein's realization, as should be evident from the conformal separation on the metric and the discussion of Sect. II. The motivation of Ref. [14] was to match the results coming from 2d gravity [52], Liouville gravity and string theory [4,33]. We return on this point with some more details later in Sect. VI.
The discussion on background and true gauge symmetries of Sect. IV applies equally well to the background split (25) of this section, so we do not repeat it in its entirety. In short, the symmetry to gauge fix is the one in which the fluctuations,ĥ µν and χ, transform as (8) at fixed backgroundsĝ µν and ϕ 0 . We gauge fix this action of Diff * as in (8) with the gauge fixing and quantities with a hat come from the unimodular background metric. This gauge fixing is chosen to reproduce (16) in the limit in which the metric fluctuations of Sect. IV and of this section are identified to the linear order. What we mean is that if we combine χ andĥ µν of this section in a symmetric tensor h µν , then (16) and (26) are the same, and they also lead to the same leading ghost action (up to O(h) interactions) because (13) and (25) are the same at the linear level (up to O(h 2 ) interactions). This choice is made to ensure the same cancellations in the same limit α = β = 1. For the computations we adopt α = 1 and β = 1 + δβ with δβ small as in the previous section to test gauge dependence. Notice that the gauge fixings (16) and (26) are not the same beyond the first order. The computation of the effective action is similar to (19) but in order to apply heat kernel techniques we factor out a local functional determinant proportional to the dilatonic background, Det[ϕ(x)δ(x − y)], via a reparametrization of the quantum fields. Moreover the computation has two traces for the (ĥ µν , χ) and ghosts sectors, but at the leading order the first trace decouples further, so it requires three distinct contributions The divergent part of Γ coming from the fluctuations of (11) distinguishes of two parts: the leading radiative correction coming from the operators of the first line, and a multiplicative renormalization of the cosmological constant V and the topological charge operators, which we implement at the linear order by treating it as a composite (local) operator. The computation requires a "Wick" rotation of the dilaton's fluctuations to imaginary space to get rid of the conformal mode instability, χ → iχ for d > 2, which is consistent with Ref. [34]. We stress, however, that there is no instability for d < 2, which is where loop integrals are regulated. 8 We temporarily postpone the renormalization of the composite operators (the second line of (11)). The result for the divergent part of (27) is The renormalization group scaleμ is a momentum scale that refers to the unimodular backgroundĝ µν , which is not Einsteinian, so we adopted a different notation to distinguish it from µ of Sect. IV. Similarly to the previous section, we want to ensure that the counter terms are evaluated on-shell. The variational principle can be applied to (11) by varying both ϕ andg µν , so we make the ansatz We introduced only the tensorê representing the equations of motion of ϕ evaluated at ϕ 0 , because at this order the equations of motion of the unimodular metric do not contribute given that they are traceless by construction. We find the coefficients We are almost ready to define the renormalization group flow. This can be done by taking ϕ 0 constant, following Ref. [14]. However, the size of the dilaton would change the beta function of Newton's constant multiplicatively, which explains the choice ϕ 0 = 1 made in Ref. [14]. There is nothing wrong in having a renormalization group flow that depends on the size of ϕ 0 , however we believe that a more elegant solution comes from first noticing that the magnitude of the scaleμ of this section refers to the unimodular background metricĝ µν . In other words, we suggest to think atμ 2 as a momentum squarê µ 2 ≈ p α p βĝ αβ . A better comparison with the results of Sect. IV would then suggest the use of a traditional Einsteinian metric, and we can construct one by combining the background dilaton and the background unimodular metric into a new Einsteinian metric which we use to define a new scale µ 2 ≈ p α p β g αβ ≈ ϕ 4 /d−2 0μ 2 . Combining everything together and using d = 2 − ζ, we find the relation It is straightforward to see that, with respect to the new scale µ, the counterterms are in form the same as the bare action (11). We thus choose µ as our scale for the renormalization group and, thanks to the fact that it comes from an Einsteinian metric, we also believe that it makes more sense to compare its effects with those observed in Sect. IV. The computation of the RG equations follows the previous section: we continue the theory to d = 2 + dimensions with G → Gµ − to get the beta function β G = G − BG 2 with B given in (30). The problem with the coefficient B is that it is gauge-dependent, which can be seen easily from the second contribution containing δβ. This is the first real obstacle that our approach encounters; we notice, however, that the gauge-dependent contribution is precisely zero for d = 2, so we know that at least in the two dimensional case the result is undoubtedly physical. We believe that the gauge dependence appears because the anomaly is well-defined only in d = 2, for example T ∝ R can only be integrated in d = 2, giving the nonlocal Polyakov action [7]. We make a point in this direction in appendix A.
In the limit d → 2 our result agrees with Refs. [9,14] which hints at the interpretation c Diff * = −25 although, as we show in Sect. V A, the value of c Diff * is scheme dependent. This is the result that one would expect from a quantization of 2d gravity "a-la-string" and it can be used to predict the critical dimension of the string worldsheet. We expand on this point further in Sect. VI. In this scheme, we also have counterterms for the topological charge and the cosmological constant operators, which we give in appendix A. Nevertheless, the RG of the volume operator becomes trivial in d = 2 as we see from the beta functional for the dimensionless cosmological constant, defined as the functional λ(ϕ) ≡ µ −d V (µ dϕ ϕ): where we have adopted the general dimension d for later use, and denoted d ϕ = (d − 2)/2 the canonical dimension of ϕ, which is interpreted as a scalar field. The reason why we have a separate beta function for the cosmological constant is that in constructing the action (11) we actually went outside the symmetry Diff * by giving to the dilaton a functional dependence in the cosmological constant operator. To discuss the beta function of λ(ϕ) in (34) notice first that the fixed functional ordinary differential equation, β λ = 0, is linear, so it is not going to "fix" a critical value for g 0 . In the weak coupling limit (Gaußian theory), the solution is simply λ(ϕ) ∼ ϕ 2d /d−2 , thus reproduces the expected limit, which can be seen in (4) upon identification of the overall constant with g 0 , λ(ϕ) = g 0 ϕ 2d /d−2 . In this sense, the volume operator √ g scales trivially, in analogy with Sect. IV, although the actual relation between the two schemes is detailed in Sect. V A. General perturbative solutions to β λ = 0 have been investigated in Ref. [31], so we do not repeat them here. We conclude this section by trying to infer some nonperturbative information on the spectrum of solutions using an extension of (34). Notice first that a similar equation has already appeared in the renormalization of 2d scalar potentials in Ref. [56] although with a different overall sign of the radiative correction, and is known to be related to nonperturbative information on the Sine-Gordon potential (in particular the Coleman phase, see the third appendix of Ref. [57]). As first step, we change the subtraction in (34) by rewriting the bare coupling G B as µ G B = 24π G/(24π − 25G), where we used the values d = 2 of (28) for simplicity. We then take = d − 2, because we want to be able to see the O( ) fixed point, and get where we defined The function ρ(d, G) has two different signs according to the value of G being either left or right of the G * ∼ O( ) fixed point, which can be interpreted as the two phases of quantum gravity [5]. For G < G * one finds ρ > 0 and G → 0 by RG evolution, where the theory presumably runs into the Gaußian phase and interpolates with the effective theory of gravity [3]; instead for G > G * one finds ρ < 0 and G → ∞ by RG evolution, so the theory is driven to the strongly interacting phase (the analogy is with high-and low-T phases in ferromagnets). Irrespectively of the phase, we can study the spectrum of fluctuations by expanding where λ * (ϕ) is the solution of β λ = 0 with arbitrary conditions, µ 0 is a reference scale, and θ is the critical exponent of the deformation Y (ϕ). Using (37) in (35) we find a Schrödinger-like equation for Y (ϕ) which can be solved as a Sturm-Liouville problem. Both phases have continuous and discrete spectrum, similarly to the Halpern-Huang potentials [58,59], but the continuous part of the spectrum can be declared unphysical following standard arguments (see, for example, Ref. [60]). We are left with a discrete spectrum for both sides, θ n = d − nd ϕ for the Gaußian phase, and θ n = d − 2nd ϕ for the strongly-interacting one, given n ∈ N. It is easy to see that n = 0 reproduces the expected scaling of the volume operator and bounds the spectrum of scaling dimensions ∆ n = d − θ n from below; there is thus a finite number of relevant directions on both sides, in agreement with the general requirements of asymptotic safety.

A. Relation between schemes
As explained in Ref. [24], it is possible to recover the value of the central charge we found in the Diff via a suitable renormalization scheme, which we do in the gauge δβ = 0 for simplicity. If we consider the equations of motion coming from the action (4), we can choose which operator is renormalized and which operator ends up having a trivial (classical) scaling on-shell. Starting from the computation that was carried out in Sect. IV, our ansatz (21) can be employed for the Diff * case as well by solving the equations of motion for the coupling g 0 . In fact, summing (28) for δβ = 0 to (A2) with V (ϕ) = 2d(d+2) (d−2) 2 g 0 ϕ 4 d−2 and imposing then it is straightforward to get the on-shell divergence which reproduces the beta function (23) for the Newton's constant through the coefficient ofR. As a consequence we get c = −19 in d = 2 as in the Diff realization of Sect. IV. Interestingly, this formula also leads to a vanishing anomalous dimension for the dilaton field in two dimensions. Note that the converse, i.e. obtaining the value c = −25 from Einstein's realization, is not possible since Weyl is not a subgroup of Diff and, therefore, there is no way of solving the equations of motion with respect to a Weyl invariant operator. In this case it is necessary to select some external Weyl-invariant operator, for example a matter operator with conformal symmetry, and arrange the counterterms so that it does not scale with the RG flow [24].

VI. CONFORMAL MODE AND d → 2
In order to appreciate the difference between the Einstein and unimodular-dilaton realizations, it is important to clarify the role of the conformal mode and the d → 2 limit. The gauge fixed two point function in h µν of the Einstein action (1) in the gauge α = β = 1 of (16) is with E µνρθ an endomorphism that depends on the curvatures, and a matrix K defined as The Green function of (40) requires the inverse of K, which becomes in which the kinematic pole at d = 2 is evident from the third term. The reason why there is a pole is because the number of degrees of freedom changes when "crossing" d = 2 by analytic continuation 9 . For arbitrary d one expects a trace and a traceless degrees of freedom. Instead, at precisely d = 2 each metric is locally conformally related to the flat metric. Naively one would expect the factor (d − 2) to multiply the traceless part, so that it does not contribute, instead it appears as a pole magnifying the trace (conformal) part. This happens because the Einstein contribution in (1) becomes a topological invariant at d = 2.
In the past, the kinematic pole has been regarded as a serious problem for the quantization of 2d gravity, because it was assumed that it should be treated as the poles coming from dimensional regularization and consistently subtracted with counterterms [13]. This seems reasonable, but, upon further scrutiny, the problem with this perspective is twofold. One the one hand the dimensional poles of regularization must be subtracted with one scheme in mind, for example in the modified minimal subtraction scheme (MS) one subtracts finite factors of π and γ, while such arbitrariness is not allowed for the kinematic pole. On the other hand, higher order dimensional poles are constrained to be consistent with structural equations, so any modification emerging from adding further poles likely breaks the constraints. Specifically, at L-loops 1 /ζ L poles arise, but they must combine with lower loop counterterms to cancel in the logarithmic RG derivative.
It is easy to discuss the case L = 2: first recall that the RG scale µ comes paired with each new pole as µ Lζ . The constraint is that higher poles come in the combination µ 2ζ /ζ 2 − 2 µ ζ /ζ 2 , where the first monomial is a genuine two loop divergence (hence µ 2ζ ), while the second is the product of a counterterm (giving one inverse power of ζ) and a one loop divergence (hence µ ζ ) [50]. The relative coefficient between the two terms must be −2 in order to cancel the logarithmic derivative, but this is simply not the case if one requires the subtraction of kinematic poles too, which is essentially the argument of Ref. [13]. We have shown in Sect. IV how kinematic poles are cancelled if one considers only on-shell gauge independent quantities for the renormalization, arguing that d dependence can be treated like the parametric N dependence for SU (N ) gauge theories. However, we did not explicitly check that the aforementioned counterterm structure is realized beyond one loop. Certainly, it would be desirable to have an explicit two loop renormalization of (1) that uses our three steps procedure to completely dismiss the claims of two loop nonrenormalizability of Ref. [13].
A different perspective on how to address the same "problem" was developed through several papers many years ago, culminating in the result of Ref. [14] with the explicit motivation of realizing the original idea of Kawai and Ninomiya [4] that 2d gravity a-la-string should arise as the limit → 0 of a theory of gravity in d = 2 + (see also [62][63][64][65] for related works). To clarify what "ala-string" means: in the modern version of the path integral of string theory, one generally switches from the Nambu-Goto to the string Polyakov's action introducing an auxiliary metric. For consistency, it is necessary to ensure that the conformal mode of the string is not quantized, nor it is propagating because of an anomaly, so one generally adopts a conformal gauge-fixing that allows to integrate the auxiliary metric only over nonconformal degrees of freedom. This results in a ghost system, generally known as bc-system, that contributes c bc = −26 to the conformal anomaly.
One way to think at this contribution to the anomaly is that it should come from a full integration of all metric degrees of freedom, if one somehow could also eliminate the trace mode effective action too. In practice, the interpretation is to see the anomaly c bc = −25 − 1, in which −1 is a contribution of a subtracted conformal mode effective action, and −25 is the expected central charge of 2d gravity. For obvious reasons, this program can only work if conformal and nonconformal (the bcsystem) are separated in the construction of the path integral. By construction, a separation procedure circumvents the problems of (42) because the two modes can be rescaled independently so to eliminate the kinematic pole. This is basically what has been done in the works leading to Ref. [14], in which it has been shown how to obtain the −25 contribution to the anomaly.
The interpretation that we suggest in this paper is that Ref. [14] does not quantize 2d gravity a-la-Einstein (1), but rather they do quantize unimodular-dilaton gravity (11), implying c Diff * = −25. A simple way to see this is to notice that the transformation (8) are precisely the action of the transformation group shown in Ref. [14], but we show that they realize Diff * instead of Diff. We stress that the two theories discussed in Sect. II are classically equivalent, so there is absolutely nothing wrong with choosing the unimodular-dilaton realization over Einstein's.
A more sophisticate analogy comes from discussing the d → 2 limit of (1), which, on the basis of the above arguments, should interpolate both with a classical conformal anomaly action and with the dilatonic action used in Ref. [14]. To see this, notice first that the would-be kinematic pole is "hidden" in the rescaling g µν = ϕ 4 /d−2g µν which leads to the actions (4) and (11). Following partly Ref. [23], we write ϕ = 1 + χ, with the choice of the magnitude of the constant background being arbitrary, and continue d = 2 + in (11) then we rescale χ to eliminate the poles χ = ( /8) 1 /2 ψ. Taking the limit → 0, we notice that the rescaling becomes g µν = ϕ 4 /d−2g µν → e 2ψg µν , and we get The action of Ref. [14] is obtained by applying the same steps, while keeping the 1 /2 and terms, which results in a more complicate function of the dilaton, L(ψ) = 1 + aψ + bψ 2 , multiplying the curvature scalar (see also a pertinent discussion in Ref. [32]). The coefficients a and b in L(ψ) are not independent, so the renormalization is nonlinear and, beyond the leading order, requires a field redefinition of ψ similarly to the 2d nonlinear sigma model [50].
The action (44) is clearly a Liouville action for ψ, but the interpretation of this paper is thatg µν is a dynamical unimodular metric. The Polyakov's anomaly action can be derived from (44) solving for ψ to get a familiar nonlocal structure [7], and one can clearly see that the Polyakov action induces a dynamics for the conformal mode by explicitly multiplying the metric by a space dependent factor.

VII. CONCLUSIONS AND FURTHER SPECULATIONS
We have explored the renormalization of metric gravity with an Einstein-Hilbert-type action in d = 2 + dimensions with the intent of framing the discussion on the four-dimensional asymptotic safety conjecture from a different angle. Quantum gravity in d = 2 is asymptotically free, thus a simple dimensional argument shows that the beta function of Newton's constant has a critical point in d = 2 + , which could represent a consistent ultraviolet completion of the theory and could circumvent the intrinsic limitations of the effective field theory approach caused by perturbative nonrenormalizability in d = 4. Our aim is to develop a perturbative framework that is reliable and gives benchmark limits for nonperturbative computations that are carried out directly in d = 4. There are several lessons that we draw from the perturbative computation in d = 2 + , which we believe should be carefully considered for the physical limit.
First and foremost, we have discussed how two different, but isomorphic, realizations of the diffeomorphisms group lead to two theories of gravity that are classically equivalent, but quantum mechanically distinct. The two realizations are different in the way in which the conformal transformations are treated: in one case we have a more traditional theory described by a metric and an Einstein-Hilbert action, while in the other case the metric is interpreted as the product of a conformal factor and another unimodular metric. From the point of view of the renormalization, the two theories allow for different critical behavior in the form of rather different beta function's coefficients, c Diff = −19 and c Diff * = −25, respectively, the difference coming from the choice of operators one compares the running of Newton's constant to. We have noticed that the first realization, c Diff = −19, has been incorrectly dismissed in the past, even though it makes perfect sense, provided that only physical on-shell quantities are considered for the renormalization. We have also noticed that the second realization, c Diff * = −25, is precisely the one that is generally associated to the string's worldsheet of string theory. Importantly, the two realizations are different because of the intrinsic symmetry, and not because of a parametrization or gauge artifacts, as was sometimes believed.
We expect that a distinction between the two realizations of the diffeomorphisms group should be made and also be very relevant for the nonperturbative approaches to quantum gravity with metric degrees of freedom. Another very important lesson comes from the necessity to go on-shell in order to have gauge and parametrization independence. This is certainly not a surprise, but we have done it in such a way that, hopefully, can be generalized to the nonperturbative approaches to quantum gravity and the asymptotic safety conjecture. We find very reasonable that our results confirm an idea, also suggested by Nink and Reuter in Ref. [23], that there are two distinct universality classes of quantum gravity, which in this work appear as the Diff and Diff * realizations. The intriguing possibility, also pushed forward in Ref. [23], is that this idea could be backed by explicit 2d conformal field theory results. In fact, it was noticed before that 19 is, among others, also a special value for the central charge [66]. It would be interesting to better understand this connection and find explicit realizations for the other special values. One such realization could involve extrinsic information of the geometry [67].
A naive question could be: why should we expect more than one universality class of quantum gravity? As we have discussed, the two realizations are classically equivalent, but the path-integrals can be different. The difference is therefore how the functional measure over the space of metrics, Dg, is constructed. Even though we do not think that framing this discussion on the different "parametrizations", exponential vs linear, is the most appropriate way to address this problem, the conclusions of Ref. [26] still apply. The two expansions, (13) and (25), end up integrating over two very different spaces of metrics, and there is absolutely no reason to believe that these two integrations should produce the same results for the effective action and β G . 10 The situation with Lorentzian metrics is probably even more complicate.
We believe that an important result of our work is the three steps procedure to handle perturbative poles in quantum gravitational framework. Our main idea is to distinguish the dimensionality d coming from vertices and contractions of the metric, which is an analog parameter to the N -dependence in SU (N ) gauge theory renormalization, from the dimensionality d = 2 − ζ at which we regularize the theory, which leads to 1 /ζ poles. A crucial point is that we can have beta functions that depend on the first d parametrically, while having a wellbehaved theory for ζ > 0. Another crucial point is that, at the end, we continue the theory above d = 2, meaning that we "reverse" the sign of ζ = − < 0. This last step is certainly the most delicate one, and it is the one that should be discussed the most if one is willing to even attempt the continuation to d = 4, which requires = 2.
Let us speculate a bit along this line. The important question that we are left with is: can we continue the asymptotically free two-dimensional gravity models above their critical dimension and still have meaningful, possibly unitary, theories? Our evidence suggests that, at least for infinitesimal , this is always possible with minor parametric constraints on d, but it would not be surprising if some serious obstruction emerges at finite values of . A similar situation has been discussed recently in a series of papers, culminating in Ref. [68], in which the O(N ) model was continued above d = 4. Ultimately, it was realized that for 4 < d < 6 the O(N ) model is not unitary because of instantons; critical quantities have exponentially small nonperturbative imaginary contributions, leading to a complex conformal field theory interpretation in d > 4. The classification of gravitational instantons started a long time ago, see for example Ref. [70], and the outcome is much more complicate than for scalar models. However, there are now compelling reasons to include them in the path integral for metric gravity and possibly in the discussion of the asymptotic safety conjecture [71].
Another speculation that we have suggested is the presence of an effective upper critical dimension, implying that the continued theory to d = 2 + , unitary or not, does not exist all the way to d → ∞, but rather it is nontrivial up to some d cr 4. This value seems to be consistently above d = 4 for both realizations of the diffeomorphisms group considered in this work (though the results for the Diff * group are gauge dependent and consequently less reliable). This is evidenced by the presence of two conformal windows for the coefficients of the beta functions continued in d [46], which, at least qualitatively, agree with previous findings by Gies et al. in Ref. [39]. A reliable procedure, perturbative or not, to compute such d cr would be very interesting because it could shine some light on the underlying mechanism that actually produces a finite value for d cr . The logic that we imply here is that, forgetting the issue of unitarity, the theory could be continued to finite values of from the pertur-bative point of view, at least in principle, so there must be a nonperturbative interplay with another entity that makes d cr finite. Such interplay could be the collision of the d = 2 + fixed point with a multicritical partner. An example of critical-multicritical collision in a simpler model is the one between a critical and a multicritical fixed points, which is believed to be responsible for the nontrivial form of the q-states Potts diagram as a function of the parameters (d, q) [72]. We hope to come back to these topics in the future.
Statement of contributions. The main results of the paper, presented in Sects. IV and V, have been obtained mainly by AU and RM. The computation was checked both by hand and with the Mathematica packages xAct [73] and xTras [74]. The discussions on symmetry and on the cosmological constant operator presented in Sects. II and V is based in part on the master thesis of FDP presented at the University of Pisa in 2020. OZ bears responsibility for any critique to the asymptotic safety programme that appeared in the introduction and any complaint should be addressed to him exclusively.
The topological charge operator was introduced in the unimodular dilaton action (11) to compensate the quantum conformal anomaly with a classical term, so that the final result is Weyl-invariant. The operator couples the scalar curvature to the dilaton field in a linear fashion in (11) O top = q d d x g ϕR , corresponding to the choice F (ϕ) = qϕ in (A2). A specific classical value of q must be chosen to cancel the quantum contribution to the anomaly as discussed in Sect. II. In order to be able to set q freely, we would like that q itself, if seen as a coupling, would not run with the renormalization group flow. This is however not the case: if we renormalize (A4) as a composite operator we get from (A3) a new on-shell contribution The obvious feature is that the subtraction of the on-shell divergence would imply a beta function for q which is nonzero β q = 0. This of course implies a nontrivial running for λ and the inability to eliminate the anomaly. However, for d = 2 the counterterm is identically zero, so q can be set freely, in agreement with the fact that the gauge dependence disappears precisely in d = 2. Realistically, if for the conjecture of asymptotic safety requires the continuation to d = 4 one could expect that it is the four dimensional conformal anomaly that actually must be canceled [75], eventually following the same steps discussed in Sect. II. This is, arguably, difficult to achieve in the -expansion, but it might be possible with functional RG methods, since they work directly in d = 4.

Appendix B: On quadratic divergences in dimensional regularization
Sometimes, it is said that dimensional regularization with modified minimal subtraction, a.k.a. MS, does not depend on the scheme, which is not entirely true. Certainly, it is unfair to other RG approaches such as functional RG [17]. In Sect. I, we have instead said that minimal subtraction is less scheme dependent for two main reasons. On the one hand, the leading and next-toleading coefficients of the RG evolution equations agree with all other approaches that do not feature an explicit cutoff scale, making them universal in this sense. Other massless schemes would be "simple" minimal subtraction, but also, for example, traditional lattice perturbation theory. On the other hand, it is often relatively easy to see explicitely that the RG equations of MS are gauge-independent, implying that they could be related to physical observables in the appropriate range of scales. This happens because the analytic continuation of dimensional regularization does not break the Ward-identities of gauge symmetry. In the background field approach used for this paper, an explicit momentum cutoff would break the Ward identities introduced by splitting the metric as in (13), and consequently the full quantum symmetry of Γ, given that the background symmetry is automatically preserved.
Obviously, the intrinsic price of dimensional regularization is that powerlaw divergences are lost in the analytic continuation, unless one modifies the scheme to incorporate them as coming from lower dimensional poles as was cleverly suggested by Jack and Jones [76,77]. This result could be an important starting point to revisit the asymptotic safety conjecture with a more "perturbative" point of view. In fact, most of the discussions on the asymptotic safety conjecture can be framed on whether the participants to the debate do or do not believe in quadratic divergences.
In this respect, perturbative methods employed at the critical dimension, such as those of this paper, have the advantage of giving clear, controlled, and weakly schemedependent answers, which might be more appropriate to settle this and other debates, if the results are compared to the logarithmic divergences that are "hidden" in the functional RG approach. A deeper and more correct interfacing between functional and dimensional methods might come from an analysis such as the one of Ref. [78], in which a functional scheme is engineered that reproduces the results of dimensional regularization to an extent (see also Ref. [79] for an earlier discussion along similar lines).
which must be solved in δ ξ h µν . We then introduce the following ansatz ξ h µν is chosen to contain the tensorial contractions with n powers of h µν . We plug the ansatz in (C2) in which we used the explicit form of the connection while keeping in mind that the metric on the left hand side must still be substituted with (C1). The relation (C4) can be solved order by order in the powers of h µν starting from δ (0) ξ h µν , which trivially has the form of the background diffeomorphism. A dependence on the parameter λ is introduced at first order by using the zeroth order solution in the third term of the first line. Each order then contains polynomials with further powers of λ. Explicitly, we find δ (0) ξ h µν = g ρν ∇ µ ξ ρ + g ρµ ∇ ν ξ ρ , δ (1) and so on, where indices on the right hand sides are raised and lowered using the background metric. It is difficult to write down resummed formulas with arbitrary λ for this expansion, but it is possible to do it for the exponential form, that corresponds to λ = 1, as shown in Ref. [14].