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 [Formula: see text]. 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 [Formula: see text] (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:with a suitable, diffeomorphism invariant action functional. This task can be addressed along various ways [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 [7-9] or on topologies of the form [10]. 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-19] for related investigations. Moreover, random walks on CDT spacetimes exhibit manifold-like behavior for long diffusion time, indicating the presence of a classical phase [20-22]. These features are typically attributed to the presence of a causal structure associated with the -factor in the topology which allows building up spacetime as a stack of spatial slices.

A second route towards asymptotic safety, reviewed in [23-28], converts the partition sum (1) into a functional renormalization group equation (FRGE) for the effective average action [29-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-47]. Starting from [48-50] renormalization group flows which resolve the difference between the background and fluctuation fields have been constructed, e.g., in [51-61] while the role of the path-integral measure has recently been discussed in [62-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 or , we evaluate the flow equation on a background geometry given by where the (intrinsic) curvature of 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 [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-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 into a lapse function , a shift vector and a metric . The later measures distances on the spatial slices defined by . For Euclidean signature this decomposition is given byAt the level of the metric tensor, this entailsAn infinitesimal coordinate transformation acts on the spacetime metric via where is the Lie derivative. This transformation induces the transformation law for the component fieldswhere the vector has been decomposed into a time-component f and a vector tangent to the spatial slice .

The construction of the flow equation for the ADM formalism uses the background field method. The quantum fields are decomposed into a fixed (but arbitrary) background and fluctuations around this background . For the lapse function and the shift vector, we resort to a linear splitThe fluctuations of may be parametrized either through a linear or an exponential splitHere 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 and have the same signature while in the linear split the fluctuations may change the signature of . The exponential split of the spatial metric in the ADM decomposition then has the same effect as in the covariant construction [39, 63–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 . Applying the method of Schur complements yieldsHere is positive by construction and the exponential split ensures that has the same signature as . 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 toThe parameter takes the value for the linear and 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 [29-32]. Besides the gravitational action, also contains suitable gauge-fixing and ghost termsHere and denote the collection of fluctuation fields and background fields, respectively. The central property of is that its dependence on the RG scale k is governed by the formally exact RG equationHere is the second variation of with respect to the fluctuation fields and the trace indicates an integration over loop momenta. The regulator provides a k-dependent mass term for the fluctuation modes with momenta and vanishes for . In the propagator , the regulator suppresses the contribution of fluctuations with momenta to the trace. The term in the numerator ensures that fluctuations with do not contribute to the trace. As a consequence the right-hand-side of Eq. (10) is finite. Moreover, the flow of is driven by fluctuations whose momenta are comparable to the RG scale k.

The FRGE realizes several welcome features. Firstly, vertices extracted from include quantum corrections resulting from integrating out fluctuations with momenta . Thus 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

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 .

Ansatz for the effective average action

We approximate the gravitational part of by the Euclidean Einstein–Hilbert action. In terms of the ADM fields the resulting action is given byHeredenotes the extrinsic curvature, is its trace, and R is the intrinsic curvature constructed from . The ansatz comprises two scale-dependent couplings, Newton’s coupling and the cosmological constant .

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 . Explicitly, we takewhere is a one-parameter family of metrics on the sphere with radius r and independent of . The curvature tensors for this background geometry then satisfyandwith 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 byAt this stage it is instructive to verify that this background is sufficient to disentangle the flow of and . For this purpose, we take the k-derivative of the ansatz (11) and subsequently set the fluctuation fields to zeroThis 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 are consistently eliminated from the Hessian . These two conditions actually fix the choice of gauge uniquely,The F and are linear in the fluctuation fields and readwhere and denotes the covariant derivative constructed from . 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 and one pair of spatial vector ghosts . Restricting to terms quadratic in the fluctuation fields, the action readswhere is the D-dimensional Laplace operator constructed from the background spacetime.

In order to compute the propagator , 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 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 constraintsand . The Jacobians coming from these decompositions are absorbed into the momentum-dependent field redefinition (B.16). The matrix elements of with respect to these component fields are computed in Appendix B and summarized in Table 4. The result shows that the field decomposition diagonalizes in field space, apart from the scalar sector containing the two fields and h.

Table 4

List of non-zero matrix elements appearing in the Hessian . 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

	(32πG_k)^-1Γk(2)_ab
h_ijh^ij	Δ-21-α2Λ_k+d^2-3d+4d(d-1)-d-22dαR¯
v_iv^i	2Δ-21-α2Λ_k+d-3d-d-22dαR¯
ψ^2	d-1dΔ-21-α2Λ_k+d-4d-d-22dαR¯
u_iu^i	2Δ-1dR¯
B^2	2Δ-2dR¯
h^2	-d-22dΔ-2Λ_k+d-4d+αdR¯+αdΛ_k
N^h	-Δ-2Λ_k+d-2dR¯
N^^2	2Δ
c¯c	(32πG_k)Δ
b¯^ib_i	(32πG_k)Δ-1dR¯

The final ingredient required in the evaluation of the FRGE is the regulator . On the background (13) the Hessian is matrix valued in field space. From Table 4 one finds that the typical matrix element has the structurewhere 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 and and . Based on the structure (24) there are two natural choices for a coarse-graining operator. Following the nomenclature introduced in [24], we defineFor notational convenience, the two regularization schemes (25) are combined by introducing a parameter r viaSetting or then corresponds to a regulator scheme of Type I and Type II, respectively. The matrix elements of are then taken as operator-valued functions depending on the coarse-graining operator. Their explicit form can be obtained by the replacement rulewhere 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).

Table 1

Parameter sets used in analyzing the dynamics encoded in the beta functions (29). The value specifies the endomorphism in the –h sector, where the coarse-graining operator is then given by

Metric fluctuations	Regulator	α	r	q_off-diag
Linear	I	0	0	0
	II_0	0	1	0
	II_1	0	1	d-2d
	II_2	0	1	d-4d
Exponential	I	1	0	0
	II_0	1	1	0
	II_1	1	1	d-2d
	II_2	1	1	d-3d

We remark that for a Type II regulator scheme the coarse-graining operators are not necessarily positive semi-definite. 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 implies that has negative eigenvalues when acting on the constant -mode and the two lowest eigenmodes in the B-sector if . In addition the linear split combined with the II regularization scheme leads to negative eigenvalues in the –h sector. The mode count for the vector ghosts is identical to the B– 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 coarse-graining 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 and on a Euclidean Friedmann–Robertson–Walker background. In this case the background geometry is characterized by while one has a non-vanishing extrinsic curvature . 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 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 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 and are constructed in Appendix C. For conciseness, we limit ourselves to the expression obtained from setting the endomorphism piece in the scalar sector spanned by 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 quantitieswhere is the anomalous dimension of Newton’s coupling. The scale dependence of the dimensionful couplings is then governed by the beta functionsDefiningone hasThe anomalous dimension of Newton’s coupling takes the formThe functions and depend on as well as the parameters . 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 areandThe remaining coefficients in are proportional to their counterpartsThis result completes the derivation of the beta functions for and on the background topology . Notably is independent of the endomorphism parameter r but retains information on the parametrization of the metric fluctuations. For , it agrees with the flow of the cosmological constant obtained on a -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 and . Throughout the discussion we will limit ourselves mostly to the case , 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 findsAll singular lines are independent of g and located at . Notably, the exponential parametrization moves the singular line closer to the origin. The singularity 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 is independent of g, this line is conveniently described by a parametrized curve obtained by equating the denominator of to zero and solving for :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.

Fig. 1

Singular loci of the beta functions (31) for the Type I, linear (left) and Type II, linear (right) setup in . The black, solid lines show the fixed singularities (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

Focusing on the region where 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 () and the linear parametrization () combined with a Type I and Type II regularization scheme. In the second setting the locus has a pole located before the first fixed singularity. As a consequence the line is not entirely shielded by the divergence of . This scenario is realized for the linear parametrization with the Type II and Type II 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 .

Subsequently, we analyze the fixed point structure of the beta functions. At a fixed point the beta functions vanish by definitionIn 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 . For eigendirections with the solutions are dragged into the fixed point for while eigendirections with repel the flow in this limit.

Fig. 2

Illustration of the r-dependence of the stability coefficients and universal product of the NGFP found for the linear split (). The real part of possesses a minimum for while Im and decrease monotonically for increasing r

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 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 () all regulators give rise to a NGFP with located at . 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 . 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 does not screen the singular line , cf. Fig. 1.

Table 2

NGFPs of the beta functions (31) evaluated for the linear split () and the exponential split () and regulators of Type I and Type II, respectively. The NGFP obtained from the optimization procedure shown in Fig. 2 is listed with the label “optimized”

Fluctuations	Regulator	g_∗	λ_∗	g_∗λ_∗	θ_1,2
Linear	Type I	0.901	0.222	0.200	1.432±2.586i
		−	−	−	−	−
	Type II_0	0.896	0.203	0.182	1.545±2.032i
		0.342	0.438	0.150	2.774	-23.89
	Type II_1	0.879	0.182	0.160	1.765±1.787i
		0.510	0.400	0.204	3.016	-13.28
	Type II_2	0.901	0.222	0.200	1.329±2.332i
		−	-	−	−	−
	Optimized	0.900	0.230	0.207	1.414±2.941i
		−	-	−	−	−
Exponential	Type I	1.049	0.249	0.262	0.444±4.041i
	Type II_0	1.050	0.249	0.261	0.342±3.855i
	Type II_2	1.050	0.249	0.261	0.342±3.855i

The exponential split () gives rise to a similar picture. In this case the specific structure of the off-diagonal contributions implies that the Type II and Type II 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 case does not support a NGFP in the physically interesting region. This indicates that the system acts rather sensitive to a change of the eigenvalue spectrum in the scalar –h sector. Thus we will limit our further considerations to the case .

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-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 . For a linear split their r-dependence is displayed in Fig. 2. While Im and are monotonically decreasing as r increases, Re develops a minimum located at . 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 () and exponential split (), respectively. The top row stems from a Type I regulator while the bottom row uses the Type II 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 , 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.

Fig. 3

Phase diagrams obtained from integrating the beta functions (29) for a Type I regulator (, top line) and a Type II regulator (, bottom line). The first column gives the result for a linear split of the spatial metric () while the second column corresponds to an exponential split (). All flows possess a NGFP providing the UV-completion of the RG trajectories

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, 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 high-energy () 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 . The trajectories with are complete in the sense that they are well defined on the entire interval . The trajectories with a positive cosmological constant terminate at 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 is lifted and replaced by a RG trajectory connecting the fixed point B and the quasi-fixed point C located at . 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 regularization scheme. The top line gives the scale dependence of and . They interpolate between the NGFP for and the quasi-fixed point C for . The lower diagrams recasts this flow in terms of the dimensionful couplings and . In this way the three scaling regimes exhibited by the solutions become even more pronounced: in the NGFP regime one has while . The classical regime is situated around . Once reaches the flow enters into a new phase where and 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.

Fig. 4

Sample RG trajectories obtained from solving the flow equations for and a Type II regulator for the initial conditions and (blue line), (magenta line), and (gold line). For the trajectories are governed by the NGFP. At intermediate scales one obtains a classical regime where the dimensionful and are independent of k. Once the gravitational instability sets in and drives the dimensionful Newton coupling and cosmological constant to zero dynamically

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, . This has profound consequences for the regularization procedure. The regulator obtained from the prescription (27) has the structurewhere the profile function satisfiesIn general, this asymptotic behavior ensures that the cutoff vanishes as . The last property fails, however, if or, equivalently , as . In this case the k-dependence of Newton’s coupling cancels the term and the cutoff remains finite as . As a result, a flow approaching the quasi-fixed point C may not integrate out all fluctuation modes, even though the limit 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 . 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 is connected to a classical regime where Newton’s coupling and the cosmological constant are essentially independent of the renormalization group scale. For specific choices of the regularization schemes, the NGFP is supplemented by a second fixed point solution constituting a saddle point in the g– plane. The interplay of the two fixed points alters the singularity structure of the beta functions which has profound consequences for the infrared behavior of the flow; see Fig. 4.

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 . At the geometrical level, the key difference in these two classes of backgrounds is that 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 background uses the kinetic terms for the gravitational fluctuations while 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-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, top-right, 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 . The labels indicate that the resulting 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

	S	V	T	TV	TTT
a_0	1	d	12d(d+1)	d-1	12(d+1)(d-2)
a_2	1-6q6	(1-6q)d6	(1-6q)d(d+1)12	(d+2)(d-3)+6δ_d,26d-q(d-1)	(d+1)(d+2)(d-5+3δ_d,2)12(d-1)-q(d-2)(d+1)2