Impact of topology in foliated quantum Einstein gravity

We use a functional renormalization group equation tailored to the Arnowitt–Deser–Misner formulation of gravity to study the scale dependence of Newton’s coupling and the cosmological constant on a background spacetime with topology \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S^1 \times S^d$$\end{document}S1×Sd. The resulting beta functions possess a non-trivial renormalization group fixed point, which may provide the high-energy completion of the theory through the asymptotic safety mechanism. The fixed point is robust with respect to changing the parametrization of the metric fluctuations and regulator scheme. The phase diagrams show that this fixed point is connected to a classical regime through a crossover. In addition the flow may exhibit a regime of “gravitational instability”, modifying the theory in the deep infrared. Our work complements earlier studies of the gravitational renormalization group flow on a background topology \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S^1 \times T^d$$\end{document}S1×Td (Biemans et al. Phys Rev D 95:086013, 2017, Biemans et al. arXiv:1702.06539, 2017) and establishes that the flow is essentially independent of the background topology.


Introduction and motivation
Asymptotic Safety, first suggested by Weinberg [3,4], constitutes a mechanism for constructing a consistent and predictive quantum theory for gravity within the framework of quantum field theory. A central goal of the program is to give meaning to the path integral over (Euclidean) metrics:  [5], e.g. by applying continuum renormalization group methods or discrete Monte Carlo techniques. Within the causal dynamical triangulations (CDT) program, reviewed in [6], the action entering the partition function (1) is taken as the Einstein-Hilbert action and the partition sum is taken on background topologies S 1 × S 3 [7][8][9] or on topologies of the form S 1 × T 3 [10]. Z is then evaluated on piecewise linear geometries constructed from elementary simplices. 1 The simplices provide a lattice regularization, making the partition sum finite. Removing the regulator by taking the continuum limit then requires a second order phase transition where the correlation length diverges. For CDT a candidate for such a phase transition has been identified in [15,16], also see [17][18][19] for related investigations. Moreover, random walks on CDT spacetimes exhibit manifoldlike behavior for long diffusion time, indicating the presence of a classical phase [20][21][22]. These features are typically attributed to the presence of a causal structure associated with the S 1 -factor in the topology which allows building up spacetime as a stack of spatial slices.
A second route towards asymptotic safety, reviewed in [23][24][25][26][27][28], converts the partition sum (1) into a functional renormalization group equation (FRGE) for the effective average action k [29][30][31]. Starting from the pioneering work [32], this program has made significant progress in demonstrating that the asymptotic safety mechanism may lead to a viable quantum theory of gravity. In particular, the existence of a non-Gaussian fixed point (NGFP), which constitutes the key element in this program, has been demonstrated in a wide range of approximations [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47]. Starting from [48][49][50] renormalization group flows which resolve the difference between the background and fluctuation fields have been con-structed, e.g., in [51][52][53][54][55][56][57][58][59][60][61] while the role of the path-integral measure has recently been discussed in [62][63][64][65]. 2 Despite their common root given by the partition sum (1), a systematic link between results obtained within the CDT program and the FRGE approach is still missing. While the spectral dimension of the resulting quantum spacetimes have been compared in [73], little is known about the relation of the two formulations. On this basis, the present work devises an FRGE study which incorporates all the essential features underlying the Monte Carlo simulations carried out within CDT. The natural continuum analogue of the foliation structure imposed on the microscopic spacetimes studied within CDT is the Arnowitt-Deser-Misner (ADM) formulation reviewed, e.g., in [74]. In this formalism spacetime is built up from a stack of spatial hypersurfaces τ on which the time-variable τ is constant. These hypersurfaces are welded together such that they fill the entire spacetime. The resulting preferred "time" direction obtained in this way plays a similar role as the causal structure implemented in CDT.
An FRGE tailored to the ADM formalism has been constructed in [75,76] and we will use this framework in the sequel. 3 This construction makes manifest use of the background field formalism. Since most CDT simulations restrict the geometries contributing to (1) to be of topology S 1 × S 2 or S 1 × S 3 , we evaluate the flow equation on a background geometry given by S 1 × S d where the (intrinsic) curvature of S d is a free parameter. Moreover, the flow is projected onto Einstein-Hilbert action which provides the weight of the partition sum (1) in the CDT framework.
Our work is complementary to the recent investigation [1,2] in the sense that it uses a different background topology. It also provides a detailed analysis on how the flow is influenced by integrating over different classes of spatial fluctuations and under a change of the regulator scheme. As a main result, we find that all cases studied in this paper admit a NGFP suitable for Asymptotic Safety. The phase diagrams obtained from integrating the flow equations are strikingly similar to the ones found for background topology S 1 × T d [1,2]. In particular, we show that there are specific combinations of parameterizing the metric fluctuations and regulating the flow equation which realizes the double-fixed point scenario found in [1] in four spacetime dimensions. In this case the RG trajectory realized by Nature, as described in [80], is well defined on all length scales. The mechanism underlying the completion of the RG trajectories in the deep infrared is closely related to the proposal of "erasing the cosmological constant through a gravitational instability", recently made in [81].
The rest of the work is organized as follows. Section 2 introduces the essential elements of the ADM formalism together with the corresponding FRGE. The beta functions governing the flow of Newton's coupling and the cosmological constant are constructed in Sect. 3 and their properties are analyzed in Sect. 4. We close with a brief discussion of our findings in Sect. 5. Technical details as regards the background geometry, the structure of the flow equation, and the evaluation of the operator traces are provided in Appendix A, Appendix B, and Appendix C, respectively.

Renormalization group flows in the ADM formalism
Our construction of the gravitational renormalization group (RG) flow is based on the FRGE for the effective average action [29][30][31] tailored to the Arnowitt-Deser-Misner (ADM) formulation [75,76]. This section summarizes the central points of the construction.

Parametrization of the fluctuation fields
The ADM formalism decomposes the spacetime metric g μν into a lapse function N (τ, y), a shift vector N i (τ, y) and a metric σ i j (τ, y). The later measures distances on the spatial slices τ defined by τ = const. For Euclidean signature this decomposition is given by At the level of the metric tensor, this entails An infinitesimal coordinate transformation acts on the spacetime metric via δg αβ = L v g αβ where L v is the Lie derivative. This transformation induces the transformation law for the component fields where the vector v α has been decomposed into a timecomponent f and a vector tangent to the spatial slice ζ i . The construction of the flow equation for the ADM formalism uses the background field method. The quantum fields N , N i , σ i j are decomposed into a fixed (but arbitrary) back-groundN ,N i ,σ i j and fluctuations around this background N ,N i ,σ i j . For the lapse function and the shift vector, we resort to a linear split The fluctuations of σ i j may be parametrized either through a linear or an exponential split Here indices are raised and lowered with the background metric. Essentially, the choice of split (6) determines the type of fluctuations admissible in the construction: the exponential split guarantees that σ i j andσ i j have the same signature while in the linear split the fluctuations may change the signature of σ i j . The exponential split of the spatial metric in the ADM decomposition then has the same effect as in the covariant construction [39,[63][64][65]: in both cases fluctuations cannot change the signature of the metric. At the level of the spacetime metric, the exponential split (6) ensures that the signature of the time direction and spatial metric remains unchanged, independently of the value of the fluctuation fields. This can be seen from computing the determinant of g μν . Applying the method of Schur complements yields det(g μν ) = N 2 det(σ i j ).
Here N 2 is positive by construction and the exponential split ensures that σ i j has the same signature asσ i j . Thus the ADM formalism in the exponential parametrization constitutes a refinement of the standard exponential parametrization by restricting the quantum fluctuations to the set which conserves the signature of the spatial and time part of the spacetime metric independently. In terms of practical computations, it is convenient to combine the linear and exponential splits according to The parameter α takes the value α = 0 for the linear and α = 1 for the exponential split. The dots represent terms containing cubic and higher powers of the fluctuation fields.
Since these terms will not contribute to the present computation we refrain from giving their explicit structure at this stage.

The functional renormalization group equation
The scale dependence of coupling constants can conveniently be obtained from the FRGE for the effective average action k [29][30][31][32]. Besides the gravitational action, k also contains suitable gauge-fixing and ghost terms Hereχ andχ denote the collection of fluctuation fields and background fields, respectively. The central property of k is that its dependence on the RG scale k is governed by the formally exact RG equation Here (2) k is the second variation of k with respect to the fluctuation fields and the trace indicates an integration over loop momenta. The regulator R k provides a k-dependent mass term for the fluctuation modes with momenta p 2 k 2 and vanishes for p 2 k 2 . In the propagator ( (2) k + R k ) −1 , the regulator suppresses the contribution of fluctuations with momenta p 2 k 2 to the trace. The term k∂ k R k in the numerator ensures that fluctuations with p 2 k 2 do not contribute to the trace. As a consequence the right-hand-side of Eq. (10) is finite. Moreover, the flow of k is driven by fluctuations whose momenta are comparable to the RG scale k.
The FRGE realizes several welcome features. Firstly, vertices extracted from k include quantum corrections resulting from integrating out fluctuations with momenta p 2 k 2 . Thus k provides a one-parameter family of effective descriptions of physics at the scale k. This realizes Wilson's idea of renormalization. Secondly, the FRGE may be used to study the RG flow and phase diagram of a theory without specifying an initial or fundamental action. This feature is particularly relevant in the context of asymptotic safety where the fundamental action is unknown a priori and arises as a fixed point of the flow. Such fixed points may be visible already in relatively simple projections of the FRGE. In the next section, we will utilize this feature and make a specific ansatz for the effective average action (9) in order to study RG flow of Newton's coupling and the cosmological constant in a setting tailored to CDT.

Einstein-Hilbert truncation on S 1 × S d
We now use the FRGE (10) to construct the beta functions governing the flow of Newton's coupling and the cosmological constant on a background topology S 1 × S d .

Ansatz for the effective average action
We approximate the gravitational part of k by the Euclidean Einstein-Hilbert action. In terms of the ADM fields the resulting action is given by Here denotes the extrinsic curvature, K ≡ σ i j K i j is its trace, and R is the intrinsic curvature constructed from σ i j . The ansatz comprises two scale-dependent couplings, Newton's coupling G k and the cosmological constant k . At this stage, it is convenient to make an explicit choice for the background fields. In the present work, we will choose a class of backgrounds with topology S 1 × S d . Explicitly, we takē whereσ S d i j is a one-parameter family of metrics on the sphere S d with radius r and independent of τ . The curvature tensors for this background geometry then satisfȳ and withR being constant. A key property of this background is the existence of a global Killing vector field ∂ τ . This feature allows one to perform a Wick rotation without generating a complex background geometry. Denoting expressions in Euclidean and Lorentzian signature by subscripts E and L, the corresponding analytic continuation is given by At this stage it is instructive to verify that this background is sufficient to disentangle the flow of G k and k . For this purpose, we take the k-derivative of the ansatz (11) and subsequently set the fluctuation fields to zero This indicates that on the background (13) the flow of Newton's coupling can be constructed from the coefficients multiplying the intrinsic background curvature while the beta function for the cosmological constant is encoded in the volume terms appearing on the left-and the right-hand-side of the FRGE. Thus it suffices to keep track of these two terms in the following. The gravitational part of the effective average action has to be complemented by a suitable gauge-fixing and ghost action. Following the strategy [1], we use the gauge freedom in such a way that all fluctuation fields including the lapse function and the shift vector acquire a relativistic dispersion relation. Moreover, terms which are of the form E δ i j p i p j are consistently eliminated from the Hessian (2) k . These two conditions actually fix the choice of gauge uniquely, The F and F i are linear in the fluctuation fields and read whereσ ≡σ i jσ i j andD i denotes the covariant derivative constructed fromσ i j . The gauge fixing (19) can be derived by adapting the harmonic gauge condition to the specific class of backgrounds and specifying the parametrization of the fluctuation fields to the one generated by the ADM decomposition. The action of the Faddeev-Popov ghosts is then constructed in the standard way. It comprises one pair of scalar ghosts c, c and one pair of spatial vector ghostsb i , b i . Restricting to terms quadratic in the fluctuation fields, the action reads where ≡ −∂ 2 τ −σ i jD iD j is the D-dimensional Laplace operator constructed from the background spacetime.
In order to compute the propagator ( k + R k ) −1 , it is useful to perform a transverse-traceless decomposition of the fluctuation fields which is adapted to the background. The shift vector is decomposed into a transverse vector u i and a scalar B For the fluctuations of the spatial metric, we resort to the standard transverse-traceless decomposition of a symmetric tensor, where the component fields are subject to the constraints andσ i jσ i j ≡σ = h. The Jacobians coming from these decompositions are absorbed into the momentum-dependent field redefinition (B.16). The matrix elements of (2) k with respect to these component fields are computed in Appendix B and summarized in Table 4. The result shows that the field decomposition diagonalizes (2) k in field space, apart from the scalar sector containing the two fieldsN and h.
The final ingredient required in the evaluation of the FRGE is the regulator R k . On the background (13) the Hessian Table 4 one finds that the typical matrix element has the structure where s takes the values 0, 1 and the dots represent a possible contribution from the cosmological constant. Moreover, q is a fixed, d-dependent numerical coefficient which depends on the field indices a, b. For example, the ghost action (20) leads to s = 0 and qc c = 0 and qb b = −1/d. Based on the structure (24) For notational convenience, the two regularization schemes (25) are combined by introducing a parameter r via Setting r = 0 or r = 1 then corresponds to a regulator scheme of Type I and Type II, respectively. The matrix elements of R k are then taken as operator-valued functions depending on the coarse-graining operator. Their explicit form can be obtained by the replacement rule where R k (z) is a scalar profile function. The parametrization of the fluctuation fields combined with the specific choices for the coarse-graining procedure, which are considered in the following, are summarized in Table 1. At this stage, we have all the ingredients for evaluating the flow equation for the ansatz (11). We remark that for a Type II regulator scheme the coarsegraining operators are not necessarily positive semidefinite. Their explicit spectrum can be constructed from the eigenvalues of the Laplacian on the d-sphere listed, e.g., in [34]. In this way one finds that setting r = 1 implies that has negative eigenvalues when acting on the constant ψ-mode and the two lowest eigenmodes in the Bsector if d = 3. In addition the linear split combined Table 1 Parameter sets used in analyzing the dynamics encoded in the beta functions (29). The value q off−diag specifies the endomorphism in theN -h sector, where the coarse-graining operator is then given by with the II 2 regularization scheme leads to negative eigenvalues in theN -h sector. The mode count for the vector ghosts is identical to the B-u i sector. The negative eigenvalues of occurring in the Type II case then suggest that the Type I regularization procedure may be preferred. The possibility for adjusting the spectrum of the coarsegraining operator by including a suitable endomorphism component may be used to implement conditions similar to the "equal lowest eigenvalue scheme" advocated in [62]. While it would be desirable to have a more complete understanding of the regulator dependence in the present case, we limit ourselves to the analysis of the cases introduced in Table 1.
We close this subsection with the following remark. In the companion paper [1] the setup (11) was used to construct the beta functions of G k and k on a Euclidean Friedmann-Robertson-Walker background. In this case the background geometry is characterized byR = 0 while one has a nonvanishing extrinsic curvatureK i j . The beta function for the Newton coupling is then read off from the extrinsic curvature terms. At the level of classical general relativity the two incarnations of G k related to the extrinsic and extrinsic curvature terms coincide due to diffeomorphism invariance of the Einstein-Hilbert action. At the level of the FRGE it is expected that the two projection schemes may lead to (slightly) different results. Firstly, introduction of the regulator R k may break the full diffeomorphism group to a subgroup so that the two projections actually construct the flow of two different coupling constants. Moreover, working on different backgrounds may result in different organization schemes for the fluctuation modes, indicating that modes integrated out at a certain scale k could be different in the two settings. The setting of this paper provides an ideal testing ground for obtaining a quantitative understanding of these effects.

Beta functions
The beta functions governing the scale dependence of G k and k are constructed in Appendix C. For conciseness, we limit ourselves to the expression obtained from setting the endomorphism piece in the scalar sector spanned byN and h to zero. Moreover, all threshold functions are evaluated with a Litim-type regulator (C.28). The result is conveniently expressed in terms of the dimensionless quantities where η is the anomalous dimension of Newton's coupling. The scale dependence of the dimensionful couplings is then governed by the beta functions Defining one has The anomalous dimension of Newton's coupling takes the form The functions B 1 and B 2 depend on λ as well as the parameters d, r, α. The terms appearing in these expressions are conveniently organized in terms of the contributions found in Appendix C.2, The explicit expressions for the building blocks are and The remaining coefficients in B 2 are proportional to their B 1 counterparts This result completes the derivation of the beta functions for g k and λ k on the background topology S 1 × S d . Notably β λ is independent of the endomorphism parameter r but retains information on the parametrization of the metric fluctuations. For α = 0, it agrees with the flow of the cosmological constant obtained on a S 1 × T d -background [2]. The beta functions (31) are the main result obtained in this section.

Properties of the renormalization group flow
The beta functions derived in the previous section explicitly retain information on the parametrization of the fluctuation fields, encoded in the parameter α, and the choice of regularization scheme, parametrized by r . Typically, these parameters have distinguished values. In the following subsection, we investigate how these choices affect the flow of g k and λ k . Throughout the discussion we will limit ourselves mostly to the case d = 3, corresponding to a four-dimensional spacetime.

Fixed point structure and phase diagrams
The beta functions (29) constitute a system of autonomous coupled non-linear differential equations. In order to understand the dynamics of the system, it is useful to first determine its fixed points and singularity structure. Singularities in the beta functions (31) can be traced back to two sources. First, there are loci in the g-λ plane where the threshold functions diverge. The location of these lines depends on α and is independent of r . Evaluating the roots of Eq. (30), one finds All singular lines are independent of g and located at λ > 0. Notably, the exponential parametrization moves the singular line λ is the counterpart of the gravitational instability discussed in Ref. [81]. Besides these fixed singularities, there are also lines in the g-λ plane where the anomalous dimension of Newton's coupling diverges. In this case the denominator in Eq. (32) vanishes. Exploiting the B 2 (λ) is independent of g, this line  (37) while the divergence of η is given by the dashed blue line. The two diagrams illustrate the two prototypical cases where the fixed singularity is screened (left) or unscreened (right) by η sing is conveniently described by a parametrized curve obtained by equating the denominator of η to zero and solving for g sing : The position of this singular line depends on the two parameters r and α. The singular lines (37) and (38) are shown in Fig. 1. At this point it is useful to distinguish between the two qualitatively different scenarios.
Focusing on the region where g > 0 one inevitably encounters a singular line when moving to positive values of λ. In the first setting, this singular locus is associated with a divergence of the anomalous dimension η given by (38). The prototypical singularity structure for this case is shown in the left diagram of Fig. 1. It is realized for the exponential parametrization (α = 1) and the linear parametrization (α = 0) combined with a Type I and Type II 2 regularization scheme. In the second setting the locus η sing has a pole located before the first fixed singularity. As a consequence the line λ sing 1 is not entirely shielded by the divergence of η. This scenario is realized for the linear parametrization α = 0 with the Type II 0 and Type II 1 regulator. It will be shown below that the different singularity structures lead to qualitatively different low-energy behaviors of the RG flow in the region λ > 0.
Subsequently, we analyze the fixed point structure of the beta functions. At a fixed point (g * , λ * ) the beta functions vanish by definition In the vicinity of a fixed point the properties of the RG flow can be studied by linearizing the beta functions. The dynamics of the linearized system is encoded in the stability matrix The stability coefficients θ are defined as minus the eigenvalues of B i j . For eigendirections with θ > 0 the solutions are dragged into the fixed point for k → ∞ while eigendirections with θ < 0 repel the flow in this limit. All implementations of the beta functions possess a Gaussian fixed point (GFP) located in the origin. This fixed point corresponds to the free theory and its critical exponents are given by the mass dimension of the dimensionful couplings. Besides the GFP the beta functions also possess non-Gaussian fixed points (NGFPs). Limiting to the physically interesting region with g > 0 located to the left of the first singular loci, a list of the NGFPs, including their position and stability coefficients, is given in Table 2. For the linear split (α = 0) all regulators give rise to a NGFP with located at g > 0, λ > 0. The complex stability coefficients with positive real part indicate that this fixed point is a spiraling attractor which captures the RG flow in its vicinity as k → ∞. In addition the Type II regulator may give rise to a second NGFP. This fixed point is a saddle point possessing one attractive and one repulsive eigendirection. Notably, the cases which possess this second fixed point coincide with the ones where η sing does not screen the singular line λ sing 1 , cf. Fig. 1.
The exponential split (α = 1) gives rise to a similar picture. In this case the specific structure of the off-diagonal contributions implies that the Type II 0 and Type II 2 regulators actually coincide, giving rise to the same beta functions. Also this case gives rise to a NGFP with complex critical exponents which acts as a UV-attractor of the flow. The distinct critical exponents accompanying the spiraling NGFPs seen in the linear and exponential setting strongly support that they correspond to two distinct universality classes. 4 The Type II 1 case does not support a NGFP in the physically interesting Table 2 NGFPs of the beta functions (31) evaluated for the linear split (α = 0) and the exponential split (α = 1) and regulators of Type I and Type II, respectively. The NGFP obtained from the optimization procedure shown in Fig. 2  region. This indicates that the system acts rather sensitive to a change of the eigenvalue spectrum in the scalarN -h sector. Thus we will limit our further considerations to the case q non−diag = 0. At this stage, it is natural to ask if there is a preferred value for the parameter r . Since a change in the regularization procedure should not affect physical quantities, a natural selection criterion for r is to minimize the sensitivity of these quantities with respect to this parameter. This optimization procedure [83][84][85][86] may then be used to find a "best value" for the parameter r . Within the present computation natural candidates for investigating the r -dependence are the stability coefficients and the universal product g * λ * . For a linear split their r -dependence is displayed in Fig. 2. While Imθ and g * λ * are monotonically decreasing as r increases, Reθ develops a minimum located at r linear opt = −0.605. The corresponding values for the position and stability coefficients of the NGFP are listed in Table 2. A comparison among the characteristic properties of the NGFP reveals that the "optimized values" turn out to be very close to the fixed point seen in the Type I regularization scheme. For the exponential split all physical quantities are monotonic functions of r . Thus in this case the principle of minimal sensitivity does not identify a preferred value for r .
Based on the rather detailed discussion of their fixed points and singularity structure, it is rather straightforward to obtain the phase diagrams resulting from integrating the flow equations (29) numerically. Our focus is on the physically interesting region where Newton's coupling is positive and containing the GFP. The resulting flows are shown in Fig. 3. The left and right column display the results obtained with a linear (α = 0) and exponential split (α = 1), respectively. The top row stems from a Type I regulator while the bottom row uses the Type II 0 regularization scheme. In all cases the GFP is marked with O while the non-Gaussian UV attractor carries the label A. The red lines mark the singular loci η sing , Eq. (38), while the blue lines connect the fixed points. All arrows point towards lower RG scales, i.e. in the direction of integrating out modes.
As expected from the results given in Table 2, the flow shows qualitative differences depending on whether the fixed point structure also contains the saddle point B (Type II 0 , linear) or just the NGFP A (Type I, linear; exponential split). In the latter case, the phase diagram is dominated by the interplay of the GFP and NGFP. The NGFP controls the highenergy (k → ∞) limit of all trajectories. Lowering the RG scale, the trajectories undergo a crossover to the GFP. In the vicinity of the GFP the trajectories develop a "classical regime" where the dimensionful couplings are almost independent of the RG scale k. Following [35], the solutions are classified according to the value of the cosmological constant in this regime: trajectories located to the left and the right of the blue separatrix give rise to a negative and positive value and are termed Type Ia and Type IIIa, respectively. The blue line separating these phases has a vanishing infrared value of k . The trajectories with 0 ≤ 0 are complete in the sense that they are well defined on the entire interval k ∈ [0, ∞]. The trajectories with a positive cosmological constant terminate at η sing at a finite value k. The presence of the saddle point B modifies this very last stage of the RG evolution. In this case the singularity η sing is lifted and replaced by a RG trajectory connecting the fixed point B and the quasi-fixed point C located at (λ, g) = (1/2, 0). The flow then follows this line and is dynamically driven into the quasi-fixed point which provides the IR-completion of these trajectories. Figure 4 displays a set of sample trajectories obtained from integrating the beta functions for the linear split with a Type II 0 regularization scheme. The top line gives the scale dependence of g k and λ k . They interpolate between the NGFP for ln k 6 and the quasi-fixed point C for ln k 1. The  Fig. 2 Illustration of the r -dependence of the stability coefficients and universal product g * λ * of the NGFP found for the linear split (α = 0). The real part of θ possesses a minimum for r linear opt = −0.605 while Imθ and g * λ * decrease monotonically for increasing r lower diagrams recasts this flow in terms of the dimensionful couplings G k and k . In this way the three scaling regimes exhibited by the solutions become even more pronounced: in the NGFP regime one has G k ∝ k −2 while k ∝ k 2 . The classical regime is situated around ln k ≈ 0. Once λ k reaches λ k ≈ 1/2 the flow enters into a new phase where G k and k are driven to zero dynamically. The RG trajectories obtained from the other cases shown in Fig. 3 are similar. The only difference consists in the absence of the final low-energy phase. Here the flow terminates in the classical regime at a finite value of k.
At this stage, a cautious remark is in order. From the inset in the top-left diagram of Fig. 4 one finds that the new IR phase comes with an anomalous dimension of Newton's coupling, η k ≈ 2. This has profound consequences for the regularization procedure. The regulator obtained from the prescription (27) has the structure where the profile function In general, this asymptotic behavior ensures that the cutoff vanishes as k → 0. The last property fails, however, if G k ∝ k 2 or, equivalently η = 2, as k → 0. In this case the kdependence of Newton's coupling cancels the k 2 term and the cutoff R k remains finite as k → 0. As a result, a flow approaching the quasi-fixed point C may not integrate out all fluctuation modes, even though the limit k → 0 is well defined. Of course, all other phases displayed in Fig. 4 are unaffected by this peculiarity.

Discussion and outlook
In this work we have studied the gravitational renormalization group (RG) flow in the Arnowitt-Deser-Misner (ADM) formalism, utilizing backgrounds with a topology S 1 × S d . This investigation is mainly motivated through the Causal Dynamical Triangulation (CDT) program where this particular topology has been used extensively in order to study properties of the gravitational partition sum through Monte Carlo simulations. The detailed results reported in the main part of the manuscript and collected in the appendices provide an important stepping stone for comparing properties of the quantum spacetimes arising within the RG and CDT framework.
Our analysis focused on the scale dependence of the (background) Newton's coupling and cosmological constant obtained from a projection of the functional renormalization group equation (10). The flow possesses a non-Gaussian fixed point (NGFP) suitable for rendering the theory asymptotically safe. The existence of this fixed point is robust with respect to changing the parametrization of the gravitational fluctuations and regularization procedure. The results summarized in Table 2 suggest that the NGFPs obtained from a linear and exponential split of the metric fluctuations belong to two different universality classes. The phase diagrams collected in Fig. 3 show that the flow emanating from the NGFP The results obtained in this work are complementary to the ones reported in [1,2] which use a very similar construction on a background topology S 1 × T d . At the geometrical level, the key difference in these two classes of backgrounds is that S 1 × S d possesses a Killing vector in the (Euclidean) time-direction which permits a Wick rotation to Lorentzian time without obtaining complex background geometries. At the level of the flow equation, the two projection schemes construct the flow of Newton's coupling based on two different interaction monomials: the S 1 × T d background uses the kinetic terms for the gravitational fluctuations while S 1 × S d uses a potential term build from the intrinsic curvature. At the classical Einstein-Hilbert action the relative coefficients of these terms are fixed by the diffeomorphism invariance. At the quantum level it is highly encouraging that the fixed point structure and phase diagrams resulting from these two projection schemes are almost identical.
At this stage constructing gravitational RG flows within the ADM formalism has achieved a similar robustness as the one encountered in comparable computations using a covariant parametrization of the metric fluctuations [34][35][36][37][38][39], at least at the background level. From the conceptual point of view, it is clear that the background field formalism based on the ADM decomposition gives rise a natural parametrization of the metric fluctuations. This parametrization is related to the one used in the covariant linear split in a non-linear way [2]. In particular, combining the ADM-split and the exponential parametrization of the fluctuations on the spatial slices ensures that the signatures of the time part and spatial part of the physical metric are independently conserved when quantum fluctuations are taken into account. In this sense, the ADM formalism provides a more refined version of the exponential parametrization recently investigated in [63,64]. An interesting consequence associated with the different parametrization schemes for the metric fluctuations is that it may shift contributions of the RG flow from background vertices to vertices containing fluctuation fields. It would be very interesting to investigate this effect in approximations of the flow equation which also takes vertices containing the fluctuation fields into account. We hope to come back to this point in the future. The beta functions (31) possess a "gravitational instability" associated with the singular line (37). In Ref. [81], it has been suggested that this type of instability could provide a dynamical solution for the cosmological constant problem through strong RG effects in the infrared. In Fig. 4 we demonstrated that this mechanism may also work in the context of pure gravity. Our general analysis identified two possible caveats to this scenario. First, the gravitational instability may be shielded by a diverging anomalous dimension. This scenario is realized by the flows displayed in the top-left, topright, and bottom-right diagram of Fig. 3. Second, the anomalous dimension may acquire a value for which the implementation of the Wilsonian RG procedure requires a modification of the (standard) regularization scheme. It would be very interesting to see if dilaton gravity, where Newton's coupling is generated dynamically, resolves these difficulties in a natural way. Table 3 Heat-kernel coefficients for fields with differential constraints on the product manifold S 1 × S d . The labels indicate that the resulting a i have been computed for spacetime scalars (S), spatial vectors (V), spatial, symmetric tensors (T), transverse spatial vectors (TV), and symmetric, transverse-traceless spatial tensors (TTT), respectively family of spheres with arbitrary but fixed radius. Keeping the radius of S d as a free parameter allows one to disentangle the volume term from an interaction term build from the intrinsic curvature scalar. Thus the background is sufficiently complex for distinguishing the two interaction monomials appearing on the left-hand-side of the projected flow equation.
In terms of the D = d + 1-dimensional geometry, the background metric is given byḡ μν = diag[1,σ i j ], whereσ i j denotes the metric on the d-sphere. Since S d is a maximally symmetric space, the Riemann tensorR i jkl and the Ricci tensorR i j constructed fromσ i j can be related to the Ricci scalarR bȳ In addition the Ricci scalar is covariantly constant,D iR = 0. A comparison with the decomposition (3) then indicates that the values of the background ADM fields are given bȳ The product nature of the background entails thatσ i j is independent of τ so that derivatives with respect to τ can be commuted freely withσ i j . Moreover, the Laplacians of the background spacetime, ≡ −ḡ μνD μDν , and on the spatial slices, −σ i jD iD j , are related by This identity allows one to express the differential operators arising in the second variation of k to differential operators constructed from the spacetime metric. Upon implementing a suitable gauge fixing, all differential operators entering in the right-hand-side of the flow equation can be combined into . The resulting traces can then be evaluated using standard heat-kernel techniques [87,88]. For the present calculation it suffices to keep track of terms containing at most two covariant derivatives. For a coarse-graining operator of the form with q being a numerical coefficient, the relevant terms in the expansion are given by (A.5) The terms indicated by the dots do not contributing to the present approximation. The heat-kernel coefficients a i depend on the index structure of the field and have been given, e.g., in [34]. The result is tabulated in Table 3. The terms proportional to δ 2,d arise from eigenmodes of the Laplacian which do not contribute to the fields appearing on the lefthand-side of the decompositions (21) and (22).

Appendix B: Evaluation of the second variations
This appendix summarizes the construction of the Hessian We start by expanding the gravitational action (11) to second order in the fluctuation fields. In order to facilitate the computation, the Einstein-Hilbert action (11) is decomposed into four interaction monomials Moreover, we use the shorthand notations to more easily treat our formulas. The expansion of the kinetic terms I 1 and I 2 around the background (13) yields the following terms quadratic in the fluctuation fields: The result is independent of the decomposition of the spatial metric in terms of background and fluctuation fields (6): the parameter α encoding the difference between the linear and exponential split does not enter these expressions. The expansion of the potential terms I 3 and I 4 is given by The α-dependence of δ 2 I 3 and δ 2 I 4 indicates that this sector receives extra terms when the exponential split is implemented.
In order to facilitate the next steps of the computation, we perform a transverse-traceless decomposition of the fluctuation fields with respect to the background. For fluctuations in the shift vectorN i we usê where the constraint ensures that u i is transverse. Analogously, the fluctuations of the spatial metric are decomposed in a transverse-traceless tensor h i j , a transverse vector v i and two scalars ψ and h, together withσ i jσ i j ≡σ = h. This change of integration variables is accompanied by non-trivial Jacobians. Their form can be deduced by evaluating the scalar products These Jacobians can then be taken into account by an additional field redefinition This rescaling also ensures that all fields appearing in the decompositions have the same dimensionality. With a slight abuse in notation, we will work with the rescaled fields in the sequel. Implementing the transverse-traceless decomposition in the kinetic terms and taking into account the rescaling (B.16) yields The potential terms become Based on these building blocks, the part of grav k quadratic in the fluctuation fields can be obtained: It is an intriguing observation that the kinetic terms do not depend on the specific decomposition of the spatial metric into background and fluctuation fields. The parameter α encoding the contributions from the exponential decomposition appears in the potential terms only.  (18) cures this problem. This is demonstrated as follows.
Following the procedure of defining building blocks, it is natural to express the gauge-fixing action in terms of the two expressions The functionals F and F i are linear in the fluctuation fields, so that these building blocks are quadratic in the fluctuations. This feature is indicated by labeling the blocks by δ 2 I 5 and δ 2 I 6 . Substituting the field decompositions (B.12) and (B.13) and invoking the rescaling (B.16) the contributions (B.22) are (B.25) Table 4 List of non-zero matrix elements appearing in the Hessian (2) k . The gravitational sector is obtained from the combination (B.26) while the ghost contributions given in the last two lines arise from (20). Each off-diagonal entry is accompanied by a suitable entry with the order of the fields reversed The matrix elements appearing in the gravitational sector of the Hessian (B.26) They are summarized in Table 4. Owing to the specific choice of the gauge fixing all differential operators combine into Laplacians of the D = d + 1-dimensional background manifold (A.3). Based on the Hessian given in Table 4, the regulator R k is constructed from the prescription (27). This completes the explicit construction of the traces appearing on the right-hand-side of the flow equation.
For a cutoff of Litim type [36], R k = (k 2 − s ) θ(k 2 − s ), to which we resort in the main part of the paper the integrals in the threshold functions can be evaluated analytically p n (w) = 1 (n + 1) (C.28) Apart from theN -h sector, all Hessians given in Table 4 have the structure 5 where 1 is the unit on the internal space and the parameter s = 0, 1 in the ghost sector and gravitational sector respectively. The contributions of the fluctuation fields to the flow can then be obtained from the following master formula: (C.30) The heat-kernel coefficients a n depend on the spin and endomorphism parameter q and are tabulated in Table 3. The anomalous dimension η is defined in (28) and w ≡ wλ k is the dimensionless version of the combination w k . For a Type II regulator c = 0 by construction.
Notably, the value α = 2 (which does not correspond to an exponential split) is special since the λ-dependence of the flow is restricted to the off-diagonal sector. Making use of the master formula (C.30) the contributions of the component fields associated with the spatial metric are (C.33) Here the argument of all threshold functions is given by (C.32). The traces containing fields B and u i , arising from the decomposition of the shift vector, do not receive a contribution from the cosmological constant. In this case the threshold functions q The results (C.33) and (C.34) complete the evaluation of the operator traces appearing in the diagonal part of the gravitational sector.

Appendix C.2.2: Off-diagonal terms in the gravitational sector
The final contribution in the gravitational sector originates from the sector spanned by the fluctuations of the lapse func-tionN and the conformal modes of the spatial metric h. In this sector the fluctuations span a full 2 × 2 block matrix in field space. Following the Type I regularization procedure, one can construct the cutoff R k in this sector and subsequently invert this block matrix in field space. The propagators appearing in this procedure are not of the simple form (C.29) but contain terms quadratic in the spacetime Laplacians. Using the specific properties of the Litim-type regulator (C.28), the result can still be expressed in terms of the threshold functions. For the Type I regulator, the resulting contribution is given by The implementation of a Type II regulator in this sector is non-trivial. This complication can be traced back to the feature that, in contrast to the other fields, the terms in the off-diagonal sector come with different relative coefficients between the Laplacians and the intrinsic curvaturē R. Inspecting the structure of