Functional Renormalization Group flows on Friedman-Lema\^{\i}tre-Robertson-Walker backgrounds

We reanalyze the construction of the gravitational functional renormalization group equation tailored to the Arnowitt-Deser-Misner formulation emphasizing its connection to the covariant formulation. The results obtained from projecting the renormalization group flow onto the Einstein-Hilbert action are reviewed in detail and we provide a novel example illustrating how the formalism may be connected to the Causal Dynamical Triangulations approach to quantum gravity.


Introduction
One of the most intriguing features of classical general relativity is the occurrence of spacetime singularities in its solutions. According to the ΛCDM model, constituting the standard model underlying cosmology, the universe is born from a phase where the energy density diverges (initial singularity). Moreover, the gravitational collapse of a star may result in the formation of a black hole accommodating a spacetime singularity in its interior. These situations are characterized by energy densities or spacetime curvatures exceeding Proceedings based on the talk given by F.S. at the workshop Lemaître, "Black Holes, Spacetime Singularities and Gravitational Waves", held at the Vatican observatory, May 9 th -12 th 2017.
their Planck values. On this basis it is expected that a description of spacetime in these extreme regions requires a more complete framework which is broadly referred to as quantum gravity.
A very conservative proposal geared towards defining a quantum theory of the gravitational interactions covering all length scales is Asymptotic Safety [1,2,3,4]. This approach stays within the framework of relativistic quantum field theory and builds on Wilson's modern formulation of renormalization. Its key ingredient is an interacting renormalization group (RG) fixed point which controls the gravitational interactions at transplanckian energy scales. The fixed point ensures that the theory does not suffer from unphysical ultraviolet (UV) divergences. Moreover, it may come with a predictive power similar to the one known from perturbatively renormalizable quantum field theories. In this setting general relativity arises dynamically when the RG flow undergoes a crossover from the interacting fixed point (quantum regime) to the free fixed point (classical regime) [5]. Moreover, the phase diagram accommodates RG trajectories which give rise to the observed values of Newton's coupling and the cosmological constant on macroscopic scales [6].
Most investigations related to asymptotically safe gravity rely on functional renormalization group methods and encode the gravitational degrees of freedom in fluctuations of the spacetime metric g µν (covariant approach). In this proceedings we focus on an alternative setting based on the Arnowitt-Deser-Misner or (3 + 1)-formulation of gravity, see [7] for a pedagogical introduction. This construction equips spacetime with a foliation structure by welding together spatial slices where the time-coordinate τ is constant. This introduces a preferred direction which may be interpreted as an (Euclidean) time direction and provides a causal structure. As a byproduct it gives access to interesting observables comprising, e.g. the expectation value for spatial volumes V 3 (τ ) as a function of time. This facilitates the direct comparison with results obtained from the Causal Dynamical Triangulations program [8] which evaluates the gravitational partition sum by Monte Carlo methods. Moreover, the close relation of the ADM framework to cosmology allows to address questions related to the very early universe based on first principles.
The proceedings is organized as follows. Sect. 2 reviews the ADM formalism and the most important elements entering into the construction of the corresponding functional renormalization group equation (FRGE). Approximate solutions of the FRGE based on the Einstein-Hilbert action are discussed in Sect. 3 and novel results relating the FRGE to the Causal Dynamical Triangulations program are presented in Sect. 4. We close by summarizing the status and perspectives of the program in Sect. 5.

The functional renormalization group in the ADM-formalism
In this section we provide the geometrical background underlying the formulation of gravity in the Arnowitt-Deser-Misner (ADM) formalism and highlight the main properties of the resulting functional renormalization group equation (FRGE) for the effective average action.

Gravity in the Arnowitt-Deser-Misner formulation
The ADM formulation of gravity starts from a D-dimensional (Euclidean) spacetime M with metric g µν carrying coordinates x µ . Subsequently, the construction introduces a time-function τ (x) which assigns a specific time to each spacetime point. The set of points with the same value of τ (x) form d-dimensional spatial slices Σ τi ≡ {x : τ (x) = τ i }. The gradient of the time function n µ ≡ N ∂ µ τ (x) defines a vector n µ normal to the spatial slices. The lapse function N (x) ensures that g µν n µ n ν = 1. One then introduces a new coordinate system x µ → (τ, y i ), i = 1, . . . , d where the y i provide coordinates on Σ τ . Defining the vector field t µ through the relation t µ ∂ µ τ = 1, the coordinate systems on neighboring spatial slices can be related by requiring that the coordinate y i is constant along the integral curves of t µ .
The tangent space at a point in M can then be decomposed into a subspace spanned by vectors tangent to Σ τ and its complement. The corresponding basis vectors are constructed from the Jacobians The normal vector thus satisfies g µν n µ e i ν = 0 and the metric induced on the spatial slices is given by A priori t µ is neither tangent nor normal to the spatial slices. The Jacobians (1) imply that its decomposition into components normal and tangent to Σ is given by where N i (τ, y) is called the shift vector. Furthermore, eq. (1) implies that the coordinate one-forms in the two coordinate systems are related by Combining the relations (4) and (3) with the normal property of n µ and the definition of σ ij , the line-element ds 2 = g µν dx µ dx ν can be recast in terms of the ADM-fields {N, N i , σ ij }: For the components of the metric tensor the decomposition (5) implies the following relations where spatial indices i, j are raised and lowered with σ ij . This entails that the relation between g µν and the ADM fields is actually non-linear. An infinitesimal coordinate transformation v µ (τ, y) acting on the metric can be expressed in terms of the Lie derivative L v Decomposing v µ = f (τ, y), ζ i (τ, y) into its temporal and spatial parts, the transformation (7) determines the transformation properties of the component fields under Diff(M) where f = f (τ ) is restricted to be independent of the spatial coordinates. Inspecting (8), one observes that the restriction f (τ, y) → f (τ ) eliminates the non-linear terms from the transformation laws so that the component fields transform linearly with respect to this subgroup.

The functional renormalization group equation on foliated spacetime
The primary tool for investigating gravitational RG flows at a non-perturbative level is the FRGE for the effective average action Γ k [9,10,11,12]. In the context of gravity, its construction hinges on the background field formalism which decomposes the physical metric g µν into a fixed but arbitrary reference backgroundḡ µν and fluctuations h µν around this background. Prominent choices for this decomposition are the linear split g µν =ḡ µν + h µν (9) and the exponential split The virtue of the background field method is that the quantum theory inherits a symmetry under background field transformations where any field transforms as a tensor of the corresponding rank. At the same timeḡ µν provides a reference scale which allows to construct a RG scale k. It is expected that the two decompositions (9) and (10) lead to different quantum theories, see [13] for an explicit demonstration in two spacetime dimensions. Heuristically, this may be understood as follows. Fixing the backgroundḡ µν , the fluctuations in the linear split may change the signature of g µν while the exponential split guarantees that both g µν andḡ µν have the same signature. Thus the set of fluctuations considered in the two settings actually differs. At the level of the path-integral this difference is reflected by different choices for the path-integral measure.
In the ADM formalism the gravitational degrees of freedom are carried by N , N i , and σ ij . Following the background field approach used in the covariant framework, these fields are decomposed into fixed but arbitrary background values (marked with bars) and fluctuations (marked with hats). The linear ADM split then uses For convenience, we denote the collection of physical fields, background fields, and the fluctuations by χ,χ, andχ, respectively and setσ ≡σ ijσ ij . Similarly to (10), the linear split of σ ij can be replaced by σ ij =σ ik eσ k j , defining the exponential ADM split. The properties of these choices are conveniently discussed by writing the determinant of the spacetime metric in terms of the ADM fields det g = N 2 det σ .
Combining this relation with the decomposition (11) one concludes that the fluctuations can not change the signature in the (Euclidean) time direction. The exponential ADM split furthermore fixes the signature of the metric on the spatial slices. Thus the exponential ADM split has similar properties as the exponential split in the covariant setting.
When setting up the FRGE in the ADM formalism, it is actually useful that there is a local map relating the fluctuations of the covariant and the ADM formulation. This allows to understand the ADM construction as a particular, background-dependent redefinition of the fluctuation fields. The explicit relation between the two settings is found by starting from the linear split (9) and performing an ADM decomposition of both g µν andḡ µν according to (6). The result allows to express the components of h µν in terms ofN ,N i , andσ ij and their background values Note that the resulting map is actually non-linear. Owed to the presence of σ ij in h 00 it involves the spatial fluctuationsσ ij to arbitrary high orders. At the linear level the map (13) entails The construction of the FRGE for the effective average action Γ k then proceeds along the same lines as in the covariant case [12] and gives [17] Here t ≡ ln(k/k 0 ) denotes the RG time, Γ (2) k is the second functional derivative of Γ k with respect to the fluctuation fieldsχ, and the STr contains an integral over loop-momenta as well as a sum over component fields. The regulator R k provides a k-dependent mass term for fluctuations with momenta p 2 k 2 and vanishes if p 2 k 2 . The interplay of R k appearing in the numerator and denominator ensures that the flow of Γ k is driven by quantum fluctuations with momentum p 2 ≈ k 2 . In this way Γ k provides a one-parameter family of effective actions encoding the dynamics of a physical system at the scale k. In connection with the exponential ADM split the flow equation (15) provides a complementary tool for probing the spacetime structure seen in Monte Carlo simulations carried out within the CDT program [8].
At this stage the following remark is in order. A key ingredient in the construction of the FRGE is the property that the k-dependent mass term giving rise to the regulator R k is quadratic in the fluctuation fields. At the same time eq. (8) indicates that the ADM fields transform non-linearly under diffeomorphisms. The FRGE then realizes only the linear part of the symmetry group as a background symmetry. Thus (15) is invariant under foliation preserving background transformations only.

Renormalization group flows in the Einstein-Hilbert truncation
We illustrate the working of the FRGE by approximating the gravitational part of Γ k by the Einstein-Hilbert action. In terms of the ADM fields one then has Γ grav where the extrinsic curvature is defined as and R and D i are the Ricci scalar and covariant derivative constructed from σ ij . The ansatz comprises two scale-dependent couplings, the cosmological constant Λ k and Newton's coupling G k . The beta functions encoding the scaledependence of these couplings may be obtained by substituting the ansatz into the FRGE and extracting suitable geometric terms on both sides. Eq. (16) suggests evaluating the flow equation at the background level, i.e. settinĝ χ = 0 after taking the necessary variations. In this case, it then suffices to keep track of terms which are at most quadratic in the fluctuation fields as these are necessary to construct the Hessian Γ k . By investigating the structure of the resulting equation, it then turns out that the computation can be simplified by using a flat Euclidean FLRW background wherē At the geometric level this entails the identities The scale-dependence of the cosmological constant and Newton's coupling can be read off from the terms proportional to the background volume and the integrated squared external curvature terms, respectively.
A key obstacle in the actual construction of flows in the ADM formalism results from the lapse N and shift N i which enter (16) as Lagrange multipliers. At the level of the FRGE, which is based on an off-shell formalism, this leads to a Hessian Γ (2) k which is degenerate. Imposing proper-time gauge which fixes N and N i to their background values does not resolve this problem. An interesting route for addressing this obstacle efficiently is to start from the background gauge-fixing procedure used in the covariant formulation Here the one-parameter family of gauge conditions F µ are linear in the fluc- where h ≡ḡ µν h µν , and α, β are free gauge-parameters. The harmonic gauge is obtained by setting α = 1, β = 1 2 while the geometric gauge sets β = 1/D and takes the limit α → 0. By construction it is clear that all gauge-fixing terms in this class preserve background diffeomorphism invariance. The analogous gauge conditions in the ADM formalism are obtained by substituting the map (13) into eq. (21). At the linear level and for the specific background (18) this results in where F µ = (F, F i ) has been decomposed in a spatial and time-component. The ghost action for this class of gauge-fixings can then be obtained in a standard way. At this stage it is instructive to combine the gravitational and gauge-fixing terms and write down the part of the action quadratic in the fluctuation fields on flat space,K = 0. Abbreviating ∆ ≡ −σ ij ∂ i ∂ j and setting α = 1 this results in This shows that there is a unique gauge choice, corresponding to harmonic gauge β = 1/2, where all derivatives combine into the D-dimensional Laplacian on flat space −∂ 2 τ +∆. 1 At the same time, the result demonstrates that there is no (local) gauge choice which preserves background diffeomorphism invariance and just provides a partial gauge-fixing for the lapseN .
The connection of (23) to cosmology is made by an additional field decomposition where the component fields are subject to the constraints In this way the FRGE on a FLRW background may directly be formulated in terms of the fields commonly used in cosmic perturbation theory.
The beta functions governing the scale-dependence of G k and Λ k are then obtained by restoring the extrinsic curvature terms in (23) and evaluating the resulting operator traces via standard heat-kernel techniques. The result has been obtained in [15] and is conveniently expressed in terms of the dimensionless couplings λ k ≡ Λ k k −2 and g k = G k k D−2 . Restricting to D = 3 + 1 dimensions it reads where and The functions B 1 and B 2 depend on λ only and read In the context of Asymptotic Safety, the most important feature of the beta functions are their fixed points where, by definition, the beta functions vanish simultaneously. In the present case, this entails β λ (g * , λ * ) = 0, β g (g * , λ * ) = 0. In the vicinity of a fixed point the properties of the RG flow can be studied by linearizing the beta functions. The linearized flow is then determined by 1 Other interesting gauge-choices are of Landau-type sending α → 0. Moreover, taking the limit β → ±∞ will gauge-fix the combination (2N +σ) sharply. For an instructive discussion at the level of the covariant theory see [20,21]. Restricting to the physically interesting region with positive Newton's coupling, g ≥ 0, the system (26) gives rise to two fixed points. The Gaussian fixed point (GFP) is located at the origin, (λ * , g * ) = (0, 0), and its stability coefficients are determined by the classical mass dimension of the couplings. In addition there is a non-Gaussian fixed point (NGFP) situated at g * = 0.785 , λ * = 0.315 , g * λ * = 0.248 .
The NGFP comes with a complex pair of critical exponents, indicating that it acts as a spiraling UV attractor for the RG trajectories in its vicinity. This is the same characteristic behavior of the NGFP found when evaluating the RG flow on foliated spacetimes using the Matsubara formalism [16,17], and a vast range of studies building on the covariant formalism, see [14] for a detailed comparison. The relations (30) entail that the dimensionful couplings at the fixed point scale according to The full phase diagram, shown in Fig. 1, is obtained by integrating the system (26) numerically. In the physical region the structure of the flow is determined by three features: a singular line associated with a divergence of the anomalous dimension η N (red line) bounds the region to the right. The flow to the left of this line is governed by the interplay of the NGFP (marked by A) and the GFP (marked by O). Flows in the vicinity of O are characterized by g k 1 and exhibit classical properties in the sense that the dimensionful couplings G k and Λ k are (approximately) k-independent in this regime. Depending on whether the trajectories flow to the left (right) of the blue separation line the value of Λ k in this classical regime is negative (positive). In [5] these classes of solutions have been termed Type Ia and Type IIIa, respectively. The single trajectory ending at the GFP gives rise to a vanishing infrared value of the cosmological constant and is referred to as Type IIa. The diagram also illustrates that the NGFP has the appropriate features for providing the high-energy completion of the RG trajectories leaving the classical regime. A rather peculiar feature of the flow results from the interplay between the specific value of the critical exponents (31) and the vicinity of the fixed point to the singular line. In this particular configuration the singular line cuts through the RG trajectories undergoing a crossover from the NGFP to the classical regime. Since this feature is absent in the companion studies on foliated backgrounds [14,19] and in the covariant computation [5] we expect that this behavior will not persist in a refined approximation.

Self-consistent classical backgrounds
A primary motivation for considering the FRGE in the presence of a foliation structure is its close relation to the Causal Dynamical Triangulations (CDT) program [8]. The latter regularizes the gravitational partition sum by introducing piecewise linear building blocks and evaluates the resulting expression using Monte Carlo methods. The presence of a causal structure is thereby essential for obtaining structures resembling a macroscopic spacetime from these elementary building blocks [22,23]. Since the simulations are capable of tracking a finite number of building blocks only, CDT geometries have periodic boundary conditions in the time direction and are necessarily compact. Their topology is fixed to either S 1 × S 3 or S 1 × T 3 where T 3 denotes the flat three-torus.
A quantity that lends itself to a direct comparison in the FRG and CDT approach is the profile of the physical volumes V 3 (τ ) as a function of Euclidean time τ . From the FRGE perspective, this information can be obtained by solving the equations of motion arising from by the effective average action Γ k . Utilizing that the ansatz (16) is actually the Einstein-Hilbert action the equations determining the self-consistent backgrounds are Einstein's equations in the presence of a cosmological constant We first focus on the solutions of this equation for RG trajectories of Type IIIa where Λ k > 0 throughout. In this case, the solutions are given by four- spheres ds 2 = r 2 dτ 2 + r 2 sin 2 (τ ) dψ 2 + sin 2 ψ dθ 2 + sin 2 ψ sin 2 θ dφ 2 , with radius r 2 = 3 Λ k . The time variable τ associated with the foliation takes values τ ∈ [0, π[, mimicking the periodic boundary conditions encountered in the CDT framework. The profile for the spatial volumes V 3 (τ ) is obtained by integrating (34) on the compact slices where τ is constant Type IIIa: Thus V 3 (τ ) possesses the sin 3 (τ )-profile characteristic for Euclidean de Sitter space. In the classical regime, where Λ k Λ 0 is independent of k, V 3 (τ ) inherits this k-independence. The volume is then set by Λ 0 and can be used to identify the underlying RG trajectory. Typical volume profiles obtained as a self-consistent solution of the effective average action in the classical regime and from the CDT program are shown as blue curves in Fig. 2, showing striking agreement. In the NGFP regime the scaling relations (32) imply that V 3 (τ ) ∝ k −3 vanishes as k → ∞. This characteristic scaling behavior may be useful to identify the NGFP in CDT simulations.
The RG trajectory of Type IIa, connecting the NGFP to the GFP, plays a special role in the phase diagram shown in Fig. 1. In this case the value of the cosmological constant vanishes in the infrared, Λ 0 = 0. The self-consistent (compact) solutions of (33) then have the topology S 1 × T 3 and the spatial slices are given by three-tori with arbitrary, time-independent volume V 3 (τ ) = V 3 . This solution together with the volume profiles measured on CDT configurations with topology S 1 × T 3 are shown as the red curves in Fig. 2. While in CDT the total volume of the configuration is fixed by the number of simplices the value V 3 determined from the effective average action is a free parameter which can be adjusted to match the CDT result. Following up on the comparison of spectral dimensions [27], the matching of the long-distance properties arising from the microscopic quantum descriptions of spacetime constitutes the next promising step in connecting the FRG and CDT frameworks in a systematic way.

Summary of results and outlook
Initial works on functional renormalization group flows in the ADM formalism [16,17,18] used a background S 1 × S d and proper-time gauge, fixing the lapse function and shift vector to their background values. In this setting the radius of S d is taken to be time-independent so that the (Euclidean) time-direction associated with the S 1 constitutes a global Killing vector field. This feature allows to analytically continue the flow equation to Lorentzian signature. The fluctuations along the time direction can then be taken into account through Matsubara sums which, in contrast to the heat-kernel computations, can be carried out in both Euclidean and Lorentzian signature. This shows that the NGFP known from Euclidean signature computations persists for Lorentzian signature as well [16].
A more systematic investigation of the gravitational RG flows on backgrounds with topology S 1 × S d has been performed in [19]. Here it was found that the spurious singularities plaguing the initial works arise due to the proper-time gauge and are removed once harmonic gauge, realizing background diffeomorphism invariance, is adopted. The systematic comparison between the fixed point structure obtained from the linear and exponential ADM splits revealed that both settings support a fixed point structure suitable for Asymptotic Safety. The critical exponents characterizing the NGFP in these two settings differ systematically in their values, supporting the expectation that the systems belong to two different gravitational universality classes.
The analysis of RG flows on flat Friedman-Lemaître-Robertson-Walker (FLRW) backgrounds [14,15] supports this picture. The NGFP seen in the gravitational setting is largely independent of the topology of the background spacetime. Supplementing the gravitational sector by minimally coupled (non self-interacting) matter fields [15] showed that many gravity-matter systems including the matter content of the standard model of particle physics actually possess NGFPs as well. Depending on the details of the matter sector the spiraling UV-attractor seen in the pure-gravity case is turned into a UVattractor with real critical exponents. The overall picture obtained from the ADM framework thereby resembles to one seen within the covariant approach [28,29,30,31] if the same approximations and coarse graining schemes are employed. For a review of the state of the art in exploring asymptotically safe gravity-matter systems we refer to the proceedings article by A. Eichhorn [32].
The formulation of RG flows on cosmological backgrounds makes the ADM framework predestined for studying cosmological fluctuation spectra from first principles. In particular, the field decomposition (24) formulates the FRGE in the field variables used in cosmic perturbation theory. Following the proposal [33], this feature gives direct access to the two-point correlation functions.
Based on these developments, we expect that Asymptotic Safety is capable of shedding light on fundamental questions related to the initial singularity and evolution of the very early universe based on first principles. Consequences of Asymptotic Safety on the formation of black holes are summarized in the proceedings of A. Platania [34].