Impact of topology in foliated Quantum Einstein Gravity

TL;DR

This paper investigates how the topology of spacetime influences the renormalization group flow in Quantum Einstein Gravity, finding a robust fixed point that supports the asymptotic safety scenario across different topologies.

Contribution

It demonstrates that the renormalization group flow and fixed points are largely independent of the background topology in foliated Quantum Einstein Gravity.

Findings

Existence of a non-trivial fixed point supporting asymptotic safety.

Fixed point is stable under different parametrizations and regulators.

Flow behavior is consistent across different background topologies.

Abstract

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 S^1xS^d. 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 S^1xT^d and establishes that the flow is essentially independent of the background topology.

Figures

Tables

Table 1. Table 1: Parameter sets used in analyzing the dynamics encoded in the beta functions ( 3.18 ). The value q off − diag subscript 𝑞 off diag q_{\rm off-diag} specifies the endomorphism in the N ^ ^ 𝑁 \hat{N} - h ℎ h –sector, where the coarse-graining operator is then given by □ = Δ + q off − diag ​ R ¯ □ Δ subscript 𝑞 off diag ¯ 𝑅 \Box=\Delta+q_{\rm off-diag}\bar{R} .

metric fluctuations	regulator	$\alpha$	$r$	$q_{\mathrm{off}-\mathrm{diag}}$
linear	I	$0$	$0$	$0$
	II0	$0$	$1$	$0$
	II1	$0$	$1$	$\frac{d-2}{d}$
	II2	$0$	$1$	$\frac{d-4+\alpha }{d}$
exponential	I	$1$	$0$	$0$
	II0	$1$	$1$	$0$
	II1	$1$	$1$	$\frac{d-2}{d}$
	II2	$1$	$1$	$\frac{d-4+\alpha }{d}$

Table 2. Table 2: NGFPs of the beta functions ( 3.20 ) evaluated for the linear split ( α = 0 𝛼 0 \alpha=0 ) and the exponential split ( α = 1 𝛼 1 \alpha=1 ) 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_{\ast }$	${\lambda }_{\ast }$	$g_{\ast }{\lambda }_{\ast }$	${\theta }_{1,2}$
linear	Type I	$0.901$	$0.222$	$0.200$	$1.432\pm 2.586i$
		$-$	$-$	$-$	$-$	$-$
	Type ${\mathrm{II}}_0$	$0.896$	$0.203$	$0.182$	$1.545\pm 2.032i$
		$0.342$	$0.438$	$0.150$	$2.774$	$-23.89$
	Type ${\mathrm{II}}_1$	$0.879$	$0.182$	$0.160$	$1.765\pm 1.787i$
		$0.510$	$0.400$	$0.204$	$3.016$	$-13.28$
	Type ${\mathrm{II}}_2$	$0.901$	$0.222$	$0.200$	$1.329\pm 2.332i$
		$-$	$-$	$-$	$-$	$-$
	optimized	$0.900$	$0.230$	$0.207$	$1.414\pm 2.941i$
		$-$	$-$	$-$	$-$	$-$
exponential	Type I	$1.049$	$0.249$	$0.262$	$0.444\pm 4.041i$
	Type ${\mathrm{II}}_0$	$1.050$	$0.249$	$0.261$	$0.342\pm 3.855i$
	Type ${\mathrm{II}}_2$	$1.050$	$0.249$	$0.261$	$0.342\pm 3.855i$

Table 3. Table 3: Heat kernel coefficients for fields with differential constraints on the product manifold S 1 × S d superscript 𝑆 1 superscript 𝑆 𝑑 S^{1}\times S^{d} . The labels indicate that the resulting a i subscript 𝑎 𝑖 a_{i} has been computed for spacetime scalars ( S 𝑆 S ), spatial vectors ( V 𝑉 V ), spatial, symmetric tensors ( T 𝑇 T ), transverse spatial vectors ( T ​ V 𝑇 𝑉 TV ), and symmetric, transverse-traceless spatial tensors ( T ​ T ​ T 𝑇 𝑇 𝑇 TTT ), respectively.

	$S$	$V$	$T$	$TV$	$TTT$
$a_0$	$1$	$d$	$\frac{1}{2}d\left(d+1\right)$	$d-1$	$\frac{1}{2}\left(d+1\right)\left(d-2\right)$
$a_2$	$\frac{1-6q}{6}$	$\frac{\left(1-6q\right)d}{6}$	$\frac{\left(1-6q\right)d\left(d+1\right)}{12}$	$\frac{\left(d+2\right)\left(d-3\right)+6{\delta }_{d,2}}{6d}$	$\frac{\left(d+1\right)\left(d+2\right)\left(d-5+3{\delta }_{d,2}\right)}{12\left(d-1\right)}$
				$-q\left(d-1\right)$	$-q\frac{\left(d-2\right)\left(d+1\right)}{2}$

Table 4. Table 4: List of non-zero matrix elements appearing in the Hessian Γ k ( 2 ) superscript subscript Γ 𝑘 2 \Gamma_{k}^{(2)} . The gravitational sector is obtained from the combination ( B.20 ) while the ghost contributions given in the last two lines arise from ( 3.9 ). Each off-diagonal entry is accompanied by a suitable entry with the order of the fields reversed.

	${\left(32\pi G_k\right)}^{-1}{\left[{\Gamma }_k^{\left(2\right)}\right]}_{ab}$
$h_{ij}{h}^{ij}$	$\Delta -2\left(1-\frac{\alpha }{2}\right){\Lambda }_k+\left(\frac{d^2-3d+4}{d\left(d-1\right)}-\frac{d-2}{2d}\alpha \right)\overline{R}$
$v_i{v}^i$	$2\left(\Delta -2\left(1-\frac{\alpha }{2}\right){\Lambda }_k+\left(\frac{d-3}{d}-\frac{d-2}{2d}\alpha \right)\overline{R}\right)$
${\psi }^2$	$\frac{d-1}{d}\left(\Delta -2\left(1-\frac{\alpha }{2}\right){\Lambda }_k+\left(\frac{d-4}{d}-\frac{d-2}{2d}\alpha \right)\overline{R}\right)$
$u_i{u}^i$	$2\left(\Delta -\frac{1}{d}\overline{R}\right)$
$B^2$	$2\left(\Delta -\frac{2}{d}\overline{R}\right)$
$h^2$	$-\frac{d-2}{2d}\left(\Delta -2{\Lambda }_k+\left(\frac{d-4}{d}+\frac{\alpha }{d}\right)\overline{R}\right)+\frac{\alpha }{d}{\Lambda }_k$
$\hat{N}h$	$-\left(\Delta -2{\Lambda }_k+\frac{d-2}{d}\overline{R}\right)$
${\hat{N}}^2$	$2\Delta$
$\overline{c}c$	$\left(32\pi G_k\right)\Delta$
${\overline{b}}^i{b}_i$	$\left(32\pi G_k\right)\left(\Delta -\frac{1}{d}\overline{R}\right)$

Equations

$$\mathcal{Z}=\int \mathcal{D}\hat{g}\exp (-S[\hat{g}])$$

$$d{s}^2=g_{\mu \nu }d{x}^{\mu }d{x}^{\nu }=N^{2}d{\tau }^2+{\sigma }_{ij}(d{y}^i+N^{i}d\tau )(d{y}^j+N^{j}d\tau ).$$

$$g_{\mu\nu}=\left(\begin{array}[]{cc}N^{2}+N_{i}N^{i}&\;\;N_{j}\\ N_{i}&\sigma_{ij}\end{array}\right)\,,\quad g^{\mu\nu}=\left(\begin{array}[]{cc}\frac{1}{N^{2}}&\;\;-\frac{N^{j}}{N^{2}}\\ -\frac{N^{i}}{N^{2}}&\;\;\sigma^{ij}+\,\frac{N^{i}\,N^{j}}{N^{2}}\end{array}\right)\,.$$

$$\begin{array}{rl}\delta N& ={\partial }_{\tau }(fN)+{\zeta }^k{\partial }_{k}N-N{N}^i{\partial }_{i}f,\\ \delta N_i& ={\partial }_{\tau }(N_{i}f)+{\zeta }^k{\partial }_k{N}_i+N_k{\partial }_i{\zeta }^k+{\sigma }_{ki}{\partial }_{\tau }{\zeta }^k+N_k{N}^k{\partial }_{i}f+N^2{\partial }_{i}f,\\ \delta {\sigma }_{ij}& =f{\partial }_{\tau }{\sigma }_{ij}+{\zeta }^k{\partial }_k{\sigma }_{ij}+{\sigma }_{jk}{\partial }_i{\zeta }^k+{\sigma }_{ik}{\partial }_j{\zeta }^k+N_j{\partial }_{i}f+N_i{\partial }_{j}f,\end{array}$$

$$N=\bar{N}+\hat{N},N_i={\bar{N}}_i+{\hat{N}}_i.$$

$$\begin{array}{rl}\text{mbox}linear:& {\sigma }_{ij}={\bar{\sigma }}_{ij}+{\hat{\sigma }}_{ij},\\ \text{mbox}exponential:& {\sigma }_{ij}={\bar{\sigma }}_{il}{\left[e^{\hat{\sigma }}\right]}^l{}_j.\end{array}$$

$$\det \left(g_{\mu \nu }\right)=N^2\det \left({\sigma }_{ij}\right).$$

$${\sigma }_{ij}\simeq {\bar{\sigma }}_{ij}+{\hat{\sigma }}_{ij}+\frac{\alpha }{2}{\hat{\sigma }}_{ik}{\bar{\sigma }}^{kl}{\hat{\sigma }}_{lj}+\dots .$$

$${\Gamma }_k[\hat{\chi };\bar{\chi }]={\Gamma }_k^{\mathrm{grav}}[\hat{\chi };\bar{\chi }]+{\Gamma }_k^{\mathrm{gauge}-\mathrm{fixing}}[\hat{\chi };\bar{\chi }]+{\Gamma }_k^{\mathrm{ghost}}[\hat{\chi };\bar{\chi }].$$

$$k{\partial }_k{\Gamma }_k=\frac{1}{2}\mathrm{Tr}\left[{\left({\Gamma }_k^{(2)}+{\mathcal{R}}_k\right)}^{-1}k{\partial }_k{\mathcal{R}}_k\right].$$

$${\Gamma }_k^{\mathrm{grav}}=\frac{1}{16\pi G_k}\int d\tau d^{d}yN\sqrt{\sigma }\left[K_{ij}{K}^{ij}-K^2-R+2{\Lambda }_k\right].$$

$$K_{ij}\equiv \frac{1}{2N}\left({\partial }_{\tau }{\sigma }_{ij}-D_i{N}_j-D_j{N}_i\right),K\equiv {\sigma }^{ij}{K}_{ij},$$

$$\bar{N}=1,{\bar{N}}_i=0,{\bar{\sigma }}_{ij}(\tau ,y)={\bar{\sigma }}_{ij}^{S^d}(y)$$

$${\bar{K}}_{ij}=0,{\bar{R}}_{ijkl}=\frac{\bar{R}}{d(d-1)}\left({\bar{\sigma }}_{ik}{\bar{\sigma }}_{jl}-{\bar{\sigma }}_{il}{\bar{\sigma }}_{jk}\right),{\bar{R}}_{ij}=\frac{1}{d}{\bar{\sigma }}_{ij}\bar{R},$$

$${\tau }_E\to -i{\tau }_L,N_E^i\to i{N}_L^i.$$

$${k{\partial }_k{\Gamma }_k^{\mathrm{grav}}|}_{\hat{\chi }=0}=k{\partial }_k\left(\frac{1}{16\pi G_k}\int d\tau d^{d}y\sqrt{\bar{\sigma }}\left[-\bar{R}+2{\Lambda }_k\right]\right).$$

$${\Gamma }_k^{\mathrm{gf}}=\frac{1}{32\pi G_k}\int d\tau d^{d}y\sqrt{\bar{\sigma }}\left[F^2+F_i{\bar{\sigma }}^{ij}{F}_j\right].$$

$$\begin{array}{rl}F=& \left[{\partial }_{\tau }\hat{N}+\overline{D^i}\widehat{N_i}-\frac{1}{2}{\partial }_{\tau }\hat{\sigma }\right]\\ F_i=& \left[{\partial }_{\tau }\widehat{N_i}-{\bar{D}}_i\hat{N}-\frac{1}{2}{\bar{D}}_i\hat{\sigma }+\overline{D^j}{\hat{\sigma }}_{ij}\right],\end{array}$$

$$\Gamma_{k}^{\mathrm{gh}}=\int d\tau d^{d}y\sqrt{\bar{\sigma}}\,\big{[}\,\bar{c}\,\Delta\,c+\bar{b}^{i}\,\left(\Delta-\tfrac{\bar{R}}{d}\right)\,b_{i}\big{]}\,,$$

$${\hat{N}}_i=u_i+{\bar{D}}_{i}B,{\bar{D}}^i{u}_i=0.$$

$${\hat{\sigma }}_{ij}=h_{ij}+{\bar{D}}_i{v}_j+{\bar{D}}_j{v}_i+\left({\bar{D}}_i{\bar{D}}_j-\frac{1}{d}{\bar{\sigma }}_{ij}{\bar{D}}^2\right)\psi +\frac{1}{d}{\bar{\sigma }}_{ij}h$$

$${\bar{D}}^i{h}_{ij}=0,{\bar{\sigma }}^{ij}{h}_{ij}=0,{\bar{D}}^i{v}_i=0,{\bar{\sigma }}^{ij}{\hat{\sigma }}_{ij}=\hat{\sigma }=h.$$

$${\left[{\Gamma }_k^{(2)}\right]}_{ab}={\left(32\pi G_k\right)}^{-s}{\left[\Delta +q\bar{R}+\dots \right]}_{ab},$$

$$\begin{array}[]{ll}{\rm Type\;I}:&\Box\equiv\Delta\,,\\[4.30554pt] {\rm Type\;II}:&\Box\equiv\Delta+q\,\bar{R}\,.\end{array}$$

$$\square \equiv \Delta +rq\bar{R}.$$

$$\square \mapsto P_k\equiv \square +R_k(\square )$$

$$\eta \equiv (G_k{)}^{-1}{\partial }_t{G}_k,{\lambda }_k\equiv {\Lambda }_k{k}^{-2},g_k\equiv G_k{k}^{d-1},$$

$$k{\partial }_k{g}_k={\beta }_g(g_k,{\lambda }_k),k{\partial }_k{\lambda }_k={\beta }_{\lambda }(g_k,{\lambda }_k).$$

$$B_{\det }\equiv d-1-(3d-2+\alpha )\lambda +2d{\lambda }^2$$

$$\begin{split}\beta_{g}=&\,(d-1+\eta)\,g\,,\\ \beta_{\lambda}=&\,(\eta-2)\lambda+\tfrac{g}{(4\pi)^{(d-1)/2}}\Big{[}-\tfrac{4(d+1)}{\Gamma((d+3)/2)}\\ &\quad+\left(d+\tfrac{d^{2}+d-2}{2(1-(2-\alpha)\lambda)}+\tfrac{2(d-1)-(3d-2+\alpha)\lambda}{B_{\rm det}}\right)\left(\tfrac{2}{\Gamma((d+3)/2)}-\tfrac{\eta}{\Gamma((d+5)/2)}\right)\Big{]}\,.\end{split}$$

Full text

1. 1 Introduction and motivation

2. 2 Renormalization group flows in the ADM-formalism

3. 2.1 Parametrization of the fluctuation fields

4. 2.2 The functional renormalization group equation

5. 3 Einstein-Hilbert truncation on $S^{1}\times S^{d}$

6. 3.1 Ansatz for the effective average action

7. 3.2 Beta functions

8. 4 Properties of the renormalization group flow

9. 4.1 Fixed point structure and phase diagrams

10. 5 Discussion and Outlook

11. A The background geometry $S^{1}\times S^{d}$

12. B Evaluation of the second variations

13. B.1 The gravitational sector

14. B.2 The gauge-fixing sector

15. C Evaluation of the operator traces

16. C.1 The master formula

17. C.2 Explicit results for the operator traces

18. C.2.1 Diagonal terms in the gravitational sector

19. C.2.2 Off-diagonal terms in the gravitational sector

20. C.2.3 Ghost contributions

1 Introduction and motivation

Asymptotic Safety, first suggested by Weinberg [3], 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

[TABLE]

with $S[\hat{g}]$ a suitable, diffeomorphism invariant action functional. This task can be addressed along various ways [4], e.g. by applying continuum renormalization group methods or discrete Monte Carlo techniques.

Within the Causal Dynamical Triangulations (CDT) program, reviewed in [5], the action entering the partition function (1.1) is taken as the Einstein-Hilbert action and the partition sum is taken on background topologies $S^{1}\times S^{3}$ [6, 7, 8] or on topologies of the form $S^{1}\times T^{3}$ [9]. ${\cal Z}$ is then evaluated on piecewise linear geometries constructed from elementary simplices.111For related work in the context of Euclidean Dynamical Triangulations see [10, 11, 12, 13]. 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 [14, 15], also see [16, 17, 18] for related investigations. Moreover, random walks on CDT spacetimes exhibit manifold-like behavior for long diffusion time, indicating the presence of a classical phase [19, 20, 21]. 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 [22, 23, 24, 25, 26, 27], converts the partition sum (1.1) into a functional renormalization group equation (FRGE) for the effective average action $\Gamma_{k}$ [28, 29, 30]. Starting from the pioneering work [31], 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 [32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46]. Starting from [47, 48, 49] renormalization group flows which resolve the difference between the background and fluctuation fields have been constructed, e.g., in [50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60] while the role of the path-integral measure has recently been discussed in [61, 62, 63, 64].222For selected works on Asymptotic Safety in gravity-matter models see [65, 66, 67, 68, 69, 70, 71].

Despite of their common root given by the partition sum (1.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 [72], 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 [73]. In this formalism spacetime is built up from a stack of spatial hypersurfaces $\Sigma_{\tau}$ on which the time-variable $\tau$ 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 [74, 75] and we will use this framework in the sequel.333For applications of this formalism to Hořava-Lifshitz gravity [76] see [77, 78]. This construction makes manifest use of the background field formalism. Since most CDT simulations restrict the geometries contributing to (1.1) to be of topology $S^{1}\times S^{2}$ or $S^{1}\times S^{3}$, we evaluate the flow equation on a background geometry given by $S^{1}\times 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.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}\times T^{d}$ [1, 2]. In particular, we show that there are specific combinations of parametrizing 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 [79], 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 [80].

The rest of the work is organized as follows. Sect. 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 about the background geometry, the structure of the flow equation, and the evaluation of the operator traces are provided in App. A, App. B, and App. C, respectively.

2 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 [28, 29, 30] tailored to the Arnowitt-Deser-Misner (ADM) formulation [74, 75]. This section summarizes the central points of the construction.

2.1 Parametrization of the fluctuation fields

The ADM formalism decomposes the spacetime metric $g_{\mu\nu}$ into a lapse function $N(\tau,y)$, a shift vector $N_{i}(\tau,y)$ and a metric $\sigma_{ij}(\tau,y)$. The later measures distances on the spatial slices $\Sigma_{\tau}$ defined by $\tau=const$. For Euclidean signature this decomposition is given by

[TABLE]

At the level of the metric tensor, this entails

[TABLE]

An infinitesimal coordinate transformation acts on the spacetime metric via $\delta g_{\alpha\beta}={\cal L}_{v}\,g_{\alpha\beta}$ where ${\cal L}_{v}$ is the Lie derivative. This transformation induces the transformation law for the component fields

[TABLE]

where the vector $v^{\alpha}$ has been decomposed into a time-component $f$ and a vector tangent to the spatial slice $\zeta^{i}$.

The construction of the flow equation for the ADM formalism uses the background field method. The quantum fields $N,N_{i},\sigma_{ij}$ are decomposed into a fixed (but arbitrary) background $\bar{N},\bar{N}_{i},\bar{\sigma}_{ij}$ and fluctuations around this background $\hat{N},\hat{N}_{i},\hat{\sigma}_{ij}$. For the lapse function and the shift vector, we resort to a linear split

[TABLE]

The fluctuations of $\sigma_{ij}$ may be parametrized either through a linear or an exponential split

[TABLE]

Here indices are raised and lowered with the background metric. Essentially, the choice of split (2.5) determines the type of fluctuations admissible in the construction: the exponential split guarantees that $\sigma_{ij}$ and $\bar{\sigma}_{ij}$ have the same signature while in the linear split the fluctuations may change the signature of $\sigma_{ij}$. The exponential split of the spatial metric in the ADM decomposition then has the same effect as in the covariant construction [38, 62, 63, 64]: in both cases fluctuations can not change the signature of the metric.

At the level of the spacetime metric, the exponential split (2.5) 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_{\mu\nu}$. Applying the method of Schur complements yields

[TABLE]

Here $N^{2}$ is positive by construction and the exponential split ensures that $\sigma_{ij}$ has the same signature as $\bar{\sigma}_{ij}$. 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

[TABLE]

The parameter $\alpha$ takes the value $\alpha=0$ for the linear and $\alpha=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.

2.2 The functional renormalization group equation

The scale-dependence of coupling constants can conveniently be obtained from the FRGE for the effective average action $\Gamma_{k}$ [28, 29, 30, 31]. Besides the gravitational action, $\Gamma_{k}$ also contains suitable gauge-fixing and ghost terms

[TABLE]

Here $\hat{\chi}$ and $\bar{\chi}$ denote the collection of fluctuation fields and background fields, respectively. The central property of $\Gamma_{k}$ is that its dependence on the RG scale $k$ is governed by the formally exact RG equation

[TABLE]

Here $\Gamma_{k}^{(2)}$ is the second variation of $\Gamma_{k}$ with respect to the fluctuation fields and the trace indicates an integration over loop momenta. The regulator ${\cal R}_{k}$ provides a $k$-dependent mass term for the fluctuation modes with momenta $p^{2}\ll k^{2}$ and vanishes for $p^{2}\gg k^{2}$. In the propagator $\left(\Gamma_{k}^{(2)}+{\cal R}_{k}\right)^{-1}$, the regulator suppresses the contribution of fluctuations with momenta $p^{2}\ll k^{2}$ to the trace. The term $k\partial_{k}{\cal R}_{k}$ in the numerator ensures that fluctuations with $p^{2}\gg k^{2}$ do not contribute to the trace. As a consequence the right-hand-side of eq. (2.9) is finite. Moreover, the flow of $\Gamma_{k}$ is driven by fluctuations whose momenta are comparable to the RG scale $k$.

The FRGE realizes several welcome features. Firstly, vertices extracted from $\Gamma_{k}$ include quantum corrections resulting from integrating out fluctuations with momenta $p^{2}\gtrsim k^{2}$. Thus $\Gamma_{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 (2.8) in order to study RG flow of Newton’s coupling and the cosmological constant in a setting tailored to CDT.

3 Einstein-Hilbert truncation on $S^{1}\times S^{d}$

We now use the FRGE (2.9) to construct the beta functions governing the flow of Newton’s coupling and the cosmological constant on a background topology $S^{1}\times S^{d}$.

3.1 Ansatz for the effective average action

We approximate the gravitational part of $\Gamma_{k}$ by the Euclidean Einstein-Hilbert action. In terms of the ADM fields the resulting action is given by

[TABLE]

Here

[TABLE]

denotes the extrinsic curvature and $R$ is the intrinsic curvature constructed from $\sigma_{ij}$. The ansatz comprises two scale-dependent couplings, Newton’s coupling $G_{k}$ and the cosmological constant $\Lambda_{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}\times S^{d}$. Explicitly, we take

[TABLE]

where $\bar{\sigma}_{ij}^{S^{d}}$ is a one-parameter family of metrics on the sphere $S^{d}$ with radius $r$ and independent of $\tau$. The curvature tensors for this background geometry then satisfy

[TABLE]

with $\bar{R}$ being constant. A key property of this background is the existence of a global Killing vector field $\partial_{\tau}$. This feature allows 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

[TABLE]

At this stage it is instructive to verify that this background is sufficient to disentangle the flow of $G_{k}$ and $\Lambda_{k}$. For this purpose, we take the $k$-derivative of the ansatz (3.1) and subsequently set the fluctuation fields to zero

[TABLE]

This indicates that on the background (3.3) 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\sqrt{\vec{p}^{\,\,2}}$ are consistently eliminated from the Hessian $\Gamma_{k}^{(2)}$. These two conditions actually fix the choice of gauge uniquely

[TABLE]

The $F$ and $F_{i}$ are linear in the fluctuation fields and read

[TABLE]

where $\hat{\sigma}\equiv\bar{\sigma}^{ij}\hat{\sigma}_{ij}$ and $\bar{D}_{i}$ denotes the covariant derivative constructed from $\bar{\sigma}_{ij}$. The gauge fixing (3.8) 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 $\bar{c},c$ and one pair of spatial vector ghosts $\bar{b}^{i},b_{i}$. Restricting to terms quadratic in the fluctuation fields, the action reads

[TABLE]

where $\Delta\equiv-\partial_{\tau}^{2}-\bar{\sigma}^{ij}\bar{D}_{i}\bar{D}_{j}$ is the $D$-dimensional Laplace operator constructed from the background spacetime.

In order to compute the propagator $(\Gamma_{k}+{\cal 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$

[TABLE]

For the fluctuations of the spatial metric, we resort to the standard transverse-traceless decomposition of a symmetric tensor,

[TABLE]

where the component fields are subject to the constraints

[TABLE]

The Jacobians coming from these decompositions are absorbed into the momentum-dependent field redefinition (B.11). The matrix elements of $\Gamma_{k}^{(2)}$ with respect to these component fields are computed in App. B and summarized in Tab. 4. The result shows that the field decomposition diagonalizes $\Gamma_{k}^{(2)}$ in field space, apart from the scalar sector containing the two fields $\hat{N}$ and $h$.

The final ingredient required in the evaluation of the FRGE is the regulator ${\cal R}_{k}$. On the background (3.3) the Hessian $\Gamma_{k}^{(2)}$ is matrix valued in field space. From Tab. 4 one finds that the typical matrix element has the structure

[TABLE]

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 (3.9) leads to $s=0$ and $q_{\bar{c}c}=0$ and $q_{\bar{b}b}=-1/d$. Based on the structure (3.13) there are two natural choices for a coarse graining operator. Following the nomenclature introduced in [23], we define

[TABLE]

For notational convenience, the two regularization schemes (3.14) are combined by introducing a parameter $r$ via

[TABLE]

Setting $r=0$ or $r=1$ then corresponds to a regulator scheme of Type I and Type II, respectively. The matrix-elements of ${\cal 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

[TABLE]

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 Tab. 1. At this stage, we have all the ingredients for evaluating the flow equation for the ansatz (3.1).

We remark that for a Type II regulator scheme the coarse graining operators $\Box$ 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 [33]. In this way one finds that setting $r=1$ implies that $\Box$ has negative eigenvalues when acting on the constant $\psi$-mode and the two lowest eigenmodes in the $B$-sector if $d=3$. In addition the linear split combined with the II2 regularization scheme leads to negative eigenvalues in the $\hat{N}$-$h$-sector. The mode count for the vector ghosts is identical to the $B$-$u_{i}$ sector. The negative eigenvalues of $\Box$ 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 [61]. 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 Tab. 1.

We close this subsection with the following remark. In the companion paper [1] the setup (3.1) was used to construct the beta functions of $G_{k}$ and $\Lambda_{k}$ on a Euclidean Friedmann-Robertson-Walker background. In this case the background geometry is characterized by $\bar{R}=0$ while one has a non-vanishing extrinsic curvature $\bar{K}_{ij}$. The beta function for the Newton’s 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 ${\cal 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.

3.2 Beta functions

The beta functions governing the scale dependence of $G_{k}$ and $\Lambda_{k}$ are constructed in App. C. For conciseness, we limit ourselves to the expression obtained from setting the endomorphism piece in the scalar sector spanned by $\hat{N}$ and $h$ to zero. Moreover, all threshold functions are evaluated with a Litim-type regulator (C.2). The result is conveniently expressed in terms of the dimensionless quantities

[TABLE]

where $\eta$ is the anomalous dimension of Newton’s coupling. The scale dependence of the dimensionful couplings is then governed by the beta functions

[TABLE]

Defining

[TABLE]

one has

[TABLE]

The anomalous dimension of Newton’s coupling takes the form

[TABLE]

The functions $B_{1}$ and $B_{2}$ depend on $\lambda$ as well as the parameters $d,r,\alpha$. The terms appearing in these expressions are conveniently organized in terms of the contributions found in App. C.2

[TABLE]

The explicit expressions for the building blocks are

[TABLE]

and

[TABLE]

The remaining coefficients in $B_{2}$ are proportional to their $B_{1}$ counterparts

[TABLE]

This result completes the derivation of the beta functions for $g_{k}$ and $\lambda_{k}$ on the background topology $S^{1}\times S^{d}$. Notably $\beta_{\lambda}$ is independent of the endomorphism parameter $r$ but retains information on the parametrization of the metric fluctuations. For $\alpha=0$, it agrees with the flow of the cosmological constant obtained on a $S^{1}\times T^{d}$-background [2]. The beta functions (3.20) are the main result obtained in this section.

4 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 $\alpha$, 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 $\lambda_{k}$. Throughout the discussion we will limit ourselves mostly to the case $d=3$, corresponding to a four-dimensional spacetime.

4.1 Fixed point structure and phase diagrams

The beta functions (3.18) 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 (3.20) can be traced back to two sources. First, there are loci in the $g$-$\lambda$-plane where the threshold functions diverge. The location of these lines depends on $\alpha$ and is independent of $r$. Evaluating the roots of eq. (3.19), one finds

[TABLE]

All singular lines are independent of $g$ and located at $\lambda>0$. Notably, the exponential parametrization moves the singular line $\lambda_{1}^{\rm sing}$ closer to the origin. The singularity $\lambda_{1}^{\rm sing}$ is the counterpart of the gravitational instability discussed in Ref. [80]. Besides these fixed singularities, there are also lines in the $g$-$\lambda$-plane where the anomalous dimension of Newton’s coupling diverges. In this case the denominator in eq. (3.21) vanishes. Exploiting the $B_{2}(\lambda)$ is independent of $g$, this line is conveniently described by a parametrized curve obtained by equating the denominator of $\eta$ to zero and solving for $g^{\rm sing}$:

[TABLE]

The position of this singular line depends on both parameters $r$ and $\alpha$.

The singular lines (4.1) and (4.2) 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 $\lambda$. In the first setting, this singular locus is associated with a divergence of the anomalous dimension $\eta$ given by (4.2). The prototypical singularity structure for this case is shown in the left diagram of Fig. 1. It is realized for the exponential parametrization ($\alpha=1$) and the linear parametrization ($\alpha=0$) combined with a Type I and Type II2 regularization scheme. In the second setting the locus $\eta^{\rm sing}$ has a pole located before the first fixed singularity. As a consequence the line $\lambda_{1}^{\rm sing}$ is not entirely shielded by the divergence of $\eta$. This scenario is realized for the linear parametrization $\alpha=0$ with the Type II0 and Type II1 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 $\lambda>0$.

Subsequently, we analyze the fixed point structure of the beta functions. At a fixed point $(g_{*},\lambda_{*})$ the beta functions vanish by definition

[TABLE]

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 ${\bf B}_{ij}=\partial_{g_{j}}\beta_{g_{i}}|_{g=g_{*}}$. The stability coefficients $\theta$ are defined as minus the eigenvalues of ${\bf B}_{ij}$. For eigendirections with $\theta>0$ the solutions are dragged into the fixed point for $k\rightarrow\infty$ while eigendirections with $\theta<0$ repel the flow in this limit.

All implementations of the beta functions posses 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 Tab. 2. For the linear split ($\alpha=0$) all regulators give rise to a NGFP with located at $g>0,\lambda>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\rightarrow\infty$. 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 $\eta^{\rm sing}$ does not screen the singular line $\lambda_{1}^{\rm sing}$, cf. Fig. 1.

The exponential split ($\alpha=1$) gives rise to a similar picture. In this case the specific structure of the off-diagonal contributions implies that the Type II0 and Type II2 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.444This agrees with the conclusion reached in [81] where it was shown that the linear and exponential parametrization lead to Liouville theories with different central charge. The Type II1 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 $\hat{N}$-$h$ sector. Thus we will limit our further considerations to the case $q_{\rm 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 [82, 83, 84, 85], 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_{*}\lambda_{*}$. For a linear split their $r$-dependence is displayed in Fig. 2. While Im$\theta$ and $g_{*}\lambda_{*}$ are monotonically decreasing as $r$ increases, Re$\theta$ develops a minimum located at $r_{\rm opt}^{\rm linear}=-0.605$. The corresponding values for the position and stability coefficients of the NGFP are listed in Tab. 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 (3.18) 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 ($\alpha=0$) and exponential split ($\alpha=1$), respectively. The top row stems from a Type I regulator while the bottom row uses the Type II0 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 $\eta^{\rm sing}$, eq. (4.2), 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 Tab. 2, the flow shows qualitative differences depending on whether the fixed point structure also contains the saddle point $B$ (Type II0, 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 ($k\rightarrow\infty$) 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 [34], 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 $\Lambda_{k}$. The trajectories with $\Lambda_{0}\leq 0$ are complete in the sense that they are well-defined on the entire interval $k\in[0,\infty]$. The trajectories with a positive cosmological constant terminate at $\eta^{\rm 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 $\eta^{\rm sing}$ is lifted and replaced by a RG trajectory connecting the fixed point $B$ and the quasi-fixed point $C$ located at $(\lambda,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.

Fig. 4 displays a set of sample trajectories obtained from integrating the beta functions for the linear split with a Type II0 regularization scheme. The top line gives the scale-dependence of $g_{k}$ and $\lambda_{k}$. They interpolate between the NGFP for $\ln k\gtrsim 6$ and the quasi-fixed point $C$ for $\ln k\lesssim 1$. The lower diagrams recasts this flow in terms of the dimensionful couplings $G_{k}$ and $\Lambda_{k}$. In this way the three scaling regimes exhibited by the solutions become even more pronounced: in the NGFP regime one has $G_{k}\propto k^{-2}$ while $\Lambda_{k}\propto k^{2}$. The classical regime is situated around $\ln k\approx 0$. Once $\lambda_{k}$ reaches $\lambda_{k}\approx 1/2$ the flow enters into a new phase where $G_{k}$ and $\Lambda_{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 $\eta_{k}\approx 2$. This has profound consequences for the regularization procedure. The regulator obtained from the prescription (3.16) has the structure

[TABLE]

where the profile function $R^{(0)}(\Delta/k^{2})$ satisfies $\lim_{k\rightarrow 0}R^{(0)}(\Delta/k^{2})=1$. In general, this asymptotic behavior ensures that the cutoff vanishes as $k\rightarrow 0$. The last property fails, however, if $G_{k}\propto k^{2}$ or, equivalently $\eta=2$, as $k\rightarrow 0$. In this case the $k$-dependence of Newton’s coupling cancels the $k^{2}$ term and the cutoff ${\cal R}_{k}$ remains finite as $k\rightarrow 0$. As a result, a flow approaching the quasi fixed point $C$ may not integrate out all fluctuation modes, even though the limit $k\rightarrow 0$ is well-defined. Of course, all other phases displayed in Fig. 4 are unaffected by this peculiarity.

5 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}\times 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 (2.9). 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 Tab. 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$-$\lambda$–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 $S^{1}\times T^{d}$. At the geometrical level, the key difference in these two classes of backgrounds is that $S^{1}\times 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}\times T^{d}$ background uses the kinetic terms for the gravitational fluctuations while $S^{1}\times 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 [33, 34, 35, 36, 37, 38], 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, 62]. 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 (3.20) possess a “gravitational instability” associated with the singular line (4.1). In Ref. [80], 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.

Acknowledgements

We thank J. Ambjørn, J. Biemans, A. Bonanno, A. Platania, R. Loll, R. Percacci and M. Reuter for helpful discussions. A. K. acknowledges financial support from an Erasmus+ fellowship. The research of F. S. is supported by the Netherlands Organisation for Scientific Research (NWO) within the Foundation for Fundamental Research on Matter (FOM) grants 13PR3137 and 13VP12.

Appendix A The background geometry $S^{1}\times S^{d}$

Throughout this work we simplify the computation by choosing a background geometry with topology $S^{1}\times S^{d}$. The size of $S^{1}$ is kept fixed while the $S^{d}$ factor denotes a one-parameter family of spheres with arbitrary but fixed radius. Keeping the radius of $S^{d}$ as a free parameter allows 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 $\bar{g}_{\mu\nu}={\rm diag}[1,\bar{\sigma}_{ij}]$, where $\bar{\sigma}_{ij}$ denotes the metric on the $d$-sphere. Since $S^{d}$ is a maximally symmetric space, the Riemann tensor $\bar{R}_{ijkl}$ and the Ricci tensor $\bar{R}_{ij}$ constructed from $\bar{\sigma}_{ij}$ can be related to the Ricci scalar $\bar{R}$ by

[TABLE]

In addition the Ricci scalar is covariantly constant, $\bar{D}_{i}\bar{R}=0$.

A comparison with the decomposition (2.2) then indicates that the values of the background ADM fields are given by

[TABLE]

The product nature of the background entails that $\bar{\sigma}_{ij}$ is independent of $\tau$ so that derivatives with respect to $\tau$ can be commuted freely with $\bar{\sigma}_{ij}$. Moreover, the Laplacians of the background spacetime, $\Delta\equiv-\bar{g}^{\mu\nu}\bar{D}_{\mu}\bar{D}_{\nu}$, and on the spatial slices, $-\bar{\sigma}^{ij}\bar{D}_{i}\bar{D}_{j}$, are related by

[TABLE]

This identity allows to express the differential operators arising in the second variation of $\Gamma_{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 $\Delta$. The resulting traces can then be evaluated using standard heat-kernel techniques [86, 87]. 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

[TABLE]

with $q$ being a numerical coefficient, the relevant terms in the expansion are given by

[TABLE]

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 [33]. The result is tabulated in Tab. 3. The terms proportional to $\delta_{2,d}$ arise from eigenmodes of the Laplacian which do not contribute to the fields appearing on the left-hand-side of the decompositions (3.10) and (3.11).

Appendix B Evaluation of the second variations

This appendix summarizes the construction of the Hessian $\Gamma_{k}^{(2)}$ on a background $S^{1}\times S^{d}$.

B.1 The gravitational sector

We start by expanding the gravitational action (3.1) to second order in the fluctuation fields. In order to facilitate the computation, the Einstein-Hilbert action (3.1) is decomposed into four interaction monomials

[TABLE]

Moreover, we use the shorthand notations

[TABLE]

to lighten our formulas. The expansion of the kinetic terms $I_{1}$ and $I_{2}$ around the background (3.3) yields the following terms quadratic in the fluctuation fields

[TABLE]

The result is independent of the decomposition of the spatial metric in terms of background and fluctuation fields (2.5): the parameter $\alpha$ 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

[TABLE]

The $\alpha$-dependence of $\delta^{2}I_{3}$ and $\delta^{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 vector $\hat{N}_{i}$ we use

[TABLE]

where the constraint ensures that $u_{i}$ is transverse. Analogously, the fluctuations of the spatial metric are decomposed in a transverse-traceless tensor $h_{ij}$, a transverse vector $v_{i}$ and two scalars $\psi$ and $h$,

[TABLE]

satisfying

[TABLE]

This change of integration variables is accompanied by non-trivial Jacobians. Their form can be deduced by evaluating the scalar products

[TABLE]

These Jacobians can then be taken into account by an additional field redefinition

[TABLE]

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.11) yields

[TABLE]

The potential terms become

[TABLE]

Based on these building blocks the part of $\Gamma_{k}^{\rm grav}$ quadratic in the fluctuation fields can be obtained as

[TABLE]

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 $\alpha$ encoding the contributions from the exponential decomposition appears in the potential terms only.

B.2 The gauge-fixing sector

Inspecting eqs. (B.12) - (B.14) more closely, one recognizes that the spatial vector field $v^{i}$ has a rather peculiar dispersion relation which does not contain spatial derivatives. This feature persists if one adopts proper time gauge where the fluctuations in the lapse function and shift vector are zero. The gauge fixing (3.7) 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

[TABLE]

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 $\delta^{2}I_{5}$ and $\delta^{2}I_{6}$. Substituting the field decompositions (B.7) and (B.8) and invoking the rescaling (B.11) the contributions (B.17) are

[TABLE]

The matrix elements appearing in the gravitational sector of the Hessian $\Gamma_{k}^{(2)}$ can then be read off from the combination

[TABLE]

They are summarized in Tab. 4.

Owed 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 Tab. 4, the regulator ${\cal R}_{k}$ is constructed from the prescription (3.16). This completes the explicit construction of the traces appearing on the right-hand-side of the flow equation.

Appendix C Evaluation of the operator traces

Based on the results collected in Tab. 4, the operator traces appearing on the right-hand-side of the flow equation can be evaluated by combining standard Mellin-transform techniques, reviewed in [23], with the heat-kernel formula (A.5). We start by giving a master formula for evaluating the operator traces in App. C.1 while the explicit results are collected in App. C.2.

C.1 The master formula

The contributions of the operator traces appearing in the FRGE are conveniently expressed in terms of the dimensionless threshold functions [31]

[TABLE]

For a cutoff of Litim type [35], $R_{k}=(k^{2}-\Delta_{s})\,\theta(k^{2}-\Delta_{s})$, to which we resort in the main part of the paper the integrals in the threshold functions can be evaluated analytically

[TABLE]

Apart from the $\hat{N}$-$h$ sector, all Hessians given in Tab. 4 have the structure555Since constant prefactors multiplying the Hessians drop out of the trace, we do not include them in the definition of ${\cal P}$.

[TABLE]

where $\mathbb{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

[TABLE]

The heat-kernel coefficients $a_{n}$ depend on the spin and endomorphism parameter $q$ and are tabulated in Tab. 3. The anomalous dimension $\eta$ is defined in (3.17) and $\widetilde{w}\equiv w\lambda_{k}$ is the dimensionless version of the combination $w\Lambda_{k}$. For a Type II regulator $c=0$ by construction.

C.2 Explicit results for the operator traces

Based on the master formula (C.4), it is rather straightforward to find the contributions of the terms listed in Tab. 4 to the projection of the RG flow.

C.2.1 Diagonal terms in the gravitational sector

The master formula (C.4) indicates that in the gravitational sector the threshold functions always appear in the combination

[TABLE]

In the diagonal sector, the argument $\widetilde{w}$ captures the dependence of the traces on the cosmological constant and it is convenient to set

[TABLE]

Notably, the value $\alpha=2$ (which does not correspond to an exponential split) is special since the $\lambda$-dependence of the flow is restricted to the off-diagonal sector. Making use of the master formula (C.4) the contributions of the component fields associated with the spatial metric are

[TABLE]

Here the argument of all threshold functions is given by (C.6).

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_{n}^{p}$ are evaluated at zero argument and the resulting trace contributions read

[TABLE]

The results (C.7) and (C.8) complete the evaluation of the operator traces appearing in the diagonal part of the gravitational sector.

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 function $\hat{N}$ and the conformal modes of the spatial metric $h$. In this sector the fluctuations span a full $2\times 2$ block matrix in field space. Following the Type I regularization procedure, one can construct the cutoff ${\cal 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.3) but contain terms quadratic in the spacetime Laplacians. Using the specific properties of the Litim-type regulator (C.2), the result can still be expressed in terms of the threshold functions. For the Type I regulator, the resulting contribution is given by

[TABLE]

where

[TABLE]

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 $\Delta$ and the intrinsic curvature $\bar{R}$. Inspecting the structure of $\Gamma_{k}^{(2)}+{\cal R}_{k}$ reveals that there are two natural candidates for choosing an endomorphism in the regulator

[TABLE]

The inclusion of the non-zero endomorphism slightly modifies the contribution of the scalar sector. For the choice $\Box_{1}$ the trace evaluates to

[TABLE]

while for the second choice, $\Box_{2}$ one obtains

[TABLE]

This result completes the evaluation of the operator traces appearing in the gravitational sector.

C.2.3 Ghost contributions

In the ghost sector all threshold functions are evaluated at zero argument and there are no terms containing $\eta$. The explicit contributions of the traces is given by

[TABLE]

Together with the formulas (C.7), (C.8), and (C.9), this result completes the evaluation of all traces appearing on the right-hand-side of the FRGE.