{'fraction_non_alphanumeric': 0.024696649077223905, 'fraction_numerical': 0.007982214126911847, 'mean_word_length': 4.607389606488435, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Renormalization Group Approach to the Continuum Limit of Matrix Models of Quantum Gravity with Preferred Foliation

This contribution is not intended as a review but, by suggestion of the editors, as a glimpse ahead into the realm of dually weighted tensor models for quantum gravity. This class of models allows one to consider a wider class of quantum gravity models, in particular one can formulate state sum models of spacetime with an intrinsic notion of foliation. The simplest one of these models is the one proposed by Benedetti and Henson, which is a matrix model formulation of two-dimensional Causal Dynamical Triangulations (CDT). In this paper we apply the Functional Renormalization Group Equation (FRGE) to the Benedetti-Henson model with the purpose of investigating the possible continuum limits of this class of models. Possible continuum limits appear in this FRGE approach as fixed points of the renormalization group flow where the size of the matrix acts as the renormalization scale. Considering very small truncations, we find fixed points that are compatible with analytically known results for CDT in two dimensions. By studying the scheme dependence of our results we find that precision results require larger truncations than the ones considered in the present work. We conclude that our work suggests that the FRGE is a useful exploratory tool for dually weighted matrix models. We thus expect that the FRGE will be a useful exploratory tool for the investigation of dually weighted tensor models for CDT in higher dimensions.


Introduction
The construction of a unified theory that contains the two most successful branches of modern physics, i.e. General Relativity (GR) and Quantum Field Theory (QFT) in a curved spacetime, as appropriate limits has been ongoing for more than eighty years and sparked many approaches to the so-called problem of quantum gravity. A complete list of these approaches goes beyond the scope of this introduction. The approaches range from the very conservative application of QFT [2,3] methods to theories of gravity and the asymptotic safety conjecture [4,5,6] over refined applications of quantization rules, such as loop quantum gravity [7], spin foams [8], group field theories [9] and tensor models [10,11,12,13] to significantly less conservative approaches like emergent gravity [14], holographic duality [15] and to searches for so-called theories of everything such as string theory [16]. None of these approaches has produced a completely satisfactory answer to the problem of quantum gravity as of now. However, most approaches possess some built-in features but all known approaches come with intrinsic short-comings that have to be overcome before qualifying the particular approach as a candidate theory of quantum gravity. This suggests to combine approaches with different built-in strengths with the goal of obtaining a new approach that combines the best of both.
The present contribution intends precisely this by combining the systematic renormalization group investigation of tensor models for quantum gravity with the success of CDT in producing phases in which the partition function is dominated by extended geometries. This combination is most straightforwardly possible when CDT is formulated as a dually weighted tensor model. Before going into detail, let us take a step back and describe the big picture schematically: Tensor models of quantum gravity are based on the idea of Euclidean lattice quantum gravity, i.e. a partition function approach in which one sums over Boltzmann factors for spacetimes that are constructed from discrete building blocks. The continuum limit of these partition functions is taken as the limit in which the size of the building blocks is sent to zero while the total volume of the spacetime is held fixed. This implies that the number of building blocks has to diverge when taking the continuum limit, thus indicating that one needs to consider the renormalization group flow of these partition functions in order to study possible phase transitions. A particularly useful tool for the systematic investigation of non-perturbative renormalization group flow is the FRGE, which takes the form of a simple one-loop equation that describes an interpolation between a bare action and the quantum effective action [18]. To apply this powerful tool to Euclidean lattice gravity it is very useful to exploit the duality between the Feynmangraphs of (un)-colored tensor models and discrete geometries. This duality allows one to identify the Feynman amplitude of the (un)-colored tensor model with the Boltzmann factor of the associated discrete gravity partition function and hence allows a translation from the tensor model action to the discrete gravity action, which takes the form of a Regge action [17]. Therefore the conjecture that the large N -limit of the tensor model translates into the continuum limit of the discrete gravity partition function. Hence, the investigation of continuum limits of lattice quantum gravity is translated into the investigation of the possible large N -limits of tensor models, which can be investigated systematically using the FRGE.
This rigorous connection between the continuum limit of Euclidean lattice quantum gravity and the large N -renormalization group flow of tensor model actions is an invaluable intrinsic feature of the tensor model approach to quantum gravity; and the systematic investigation of these continuum limits with the FRGE is particularly convenient. Unfortunately, the extended geometries that are approximated in the continuum limits that have been investigated so far possess dimension two or less. In other words, so far no state sum model of discrete geometry is known to coarse grain to a model of extended spacetime geometry in more than two dimensions.
There are however numerical indications that ddimensional CDT and its counter part Euclidean Dynamical Triangulations (EDT) do coarse grain to ex-tended higher dimensional geometries (for 2 ≤ d ≤ 4). This can be heuristically understood as the fact that the foliation in CDT and the volume term in EDT implement additional terms in the Boltzmann factor for discrete geometry which change the universality class of the model. Moreover, there exist tensor model formulations of CDT and EDT in the literature. The novelty in these models is that they possess a nontrivial propagator, which implements a dual weighting of the Feynman graphs of these tensor models. Hence, one can use the FRGE to investigate the continuum limits of CDT and EDT by studying the renormalization group flow of tensor models with dual weights. This is the motivation for the work presented in the present contribution.
As a first step, we consider a dually weighted matrix model proposed by Benedetti and Henson whose partition function is dual to two-dimensional CDT [22]. By doing this we follow a strategy that was used when first applying the FRGE to tensor models [21], where matrix models for two-dimensional Euclidean quantum gravity were considered to introduce the setup, develop the technique and to compare with the analytic results known from constructive approaches to two-dimensional Euclidean quantum gravity, which serve as a bench mark. This allows us to test a setup (the systematic FRGE investigation to dually weighted tensor models) that is readily available in higher dimensions, in particular in 3+1 dimensions [22], but at the same time is understood analytically, providing benchmark results for the FRGE calculation which we can use to gauge this setup.
This contribution is organized as follows: In the following section (section 2) we provide the necessary background on dually weighted tensor models, the particular model proposed by Benedetti and Henson and the foundations of the application of the FRGE to tensor models. We provide the derivation of the beta functions in section 3. We perform a fixed point analysis and study of scheme dependence in section 4. We summarize our results in section 5 and briefly discuss their implications for future investigations on dually weighted tensor models for quantum gravity.

Preliminaries
Random tensor models are by now an established approach to Euclidean quantum gravity. However, to fully appreciate the way in which dually weighted matrix models provide an approach to quantum gravity with a preferred time slicing one needs to take a step back and consider the foundations of random tensor models.

Tensor Models and Dual Weights
The random tensor model approach to quantum gravity is based on the basic observation that the Feynman graphs of so-called uncolored random tensor models possess a geometric interpretation in terms of tessellations of piece-wise linear pseudo-manifolds, as do some so-called colored tensor models. The uncolored models are defined for tensors T i1i2,...i k and their complex conjugates T i1i2i3...i k through the symmetry of the action S[T,T ] under the U k (N ) transformations This symmetry implies that the action can be expanded in terms of generalized trace invariants in which the first index i 1 of each tensor T must be contracted with the first index of a complex conjugated tensorT , and similarly the second index i 2 and all further indices i l . These generalized traces can be represented as colored graphs where each tensor T is represented by a white vertex and each complex conjugate tensorT is represented by a black vertex and each contraction of vertices by an index i l is represented by an edge of color l connecting the vertices associated with the two contracted tensors. Such colored graphs are then dual to piecewise linear pseudo-manifolds: Each vertex is associated with a (k − 1)-simplex and the adjacent edges are associated with a gluing of the colored (k − 2)-simplices in the boundary of the two (k − 1) simplices. Moreover, closed two-colored sub-graphs are associated with the gluing of (k − 3)-simplices in the boundary of the associated (k − 2)-simplices. Analogously, closed three-and more-colored subgraphs are associated with the gluing of simplices of even lower dimension. The generalized trace-invariants of a rank k tensor model can thus be interpreted as tessellations of piecewise linear (k − 1)dimensional manifolds.
We can now perform the analogous identification for the Feynman graphs generated by the rank k random tensor model through realizing that the Feynman graphs of a rank k tensor model possess a graphical representation in terms of (k+1)-colored graphs, where a new color is associated with the propagator. This provides the desired geometric interpretation of the Feynman graphs of an uncolored rank k tensor model with tessellations ∆ of piecewise linear pseudo-manifolds of dimension k. It follows that the partition function of these random tensor models possesses a geometric interpretation where A(∆) denotes the Feynman amplitude associated with the Feynman graph dual to ∆. This resembles the random lattice partition function for Euclidean quantum gravity when the gravity action S E [∆, a] is identified with − ln(A(∆)). The Feynman amplitude depends on the details of the random tensor model, but one can generally say that they depend on the number of N k of k simplices and the number N k−2 of k − 2 simplices in ∆ as well as the tensor size N and the coupling constants λ i . An example amplitude for k = 3 with one coupling is where the couplings N, λ possess a simple relation with the couplings in the Regge-expression of General Relativity in three dimensions S R [∆] = κ 3 N 3 −κ 1 N 1 . Hence, κ 1 = ln(N ) and κ 3 = 3 2 ln(N ) − 1 2 ln(λλ) establishes a relation with the discrete General Relativity coupling constants.
The total volume is where V o is the filling factor of the geometric building blocks. Hence, one can take the lattice continuum limit a → 0 at fixed total volume V by tuning to a point where the expectation value of the total volume diverges. This requires that N d diverges. It turns out that this in turn requires that N → ∞. However, to obtain simultaneously a finite value of the total volume and of Newton's constant, one needs to tune λ and N simultaneously. Since Z diverges for N → ∞ one can only obtain a finite result when λ approaches a critical point λ * as N → ∞ is approached. Hence, we can write down the required behaviour of λ(N ) = λ * + c N −θ , where c is an arbitrary constant and θ the critical exponent. In other words, the conjecture is that the continuum limit of lattice quantum gravity can be investigated by studying the critical points in the large N behaviour of random tensor models.
So far we have only considered a canonical quadratic term T i1i2...i kT i1i2...i k , as is implied by U k (N ) invariance. This kinetic term leads to an index-independent propagator ∝ δ i1j1 ...δ i k j k , so each closed loop of indices will contribute with a factor of N to the amplitude, but can not depend on the number of vertices that are crossed when going around this loop. However, we will see shortly that such a dependence of the amplitude can be motivated geometrically. To construct tensor models whose amplitude depends non-trivially on the number of vertices crossed by a closed index loop. Before motivating these so-called "dually weighted" tensor models, we will consider the general setup of the FRGE for tensor models.

Application of the FRGE to tensor models
One of the most convenient tools to investigate critical behaviour is the functional renormalization group equation (FRGE). In the usual setting the FRGE describes how the effective average action changes when the IR suppression scale k is changed. This IR-suppression scale is introduced through a modification of the bare action by the scale dependent mass term ∆ k S[φ] = 1 2 φ.R k .φ, which is designed to give a mass of O(k) to modes in the IR of the scale k while not significantly affecting modes in the UV of this scale. Heuristically, one can argue as follows: ∆ k S[φ] dominates the path integral in the limit k → ∞ and hence the saddle point approximation of the path integral becomes exact in this limit and shows that the effective average action coincides with the bare action when k → ∞. Hence, one finds critical points as UV fixed points of the FRGE and can study the critical behaviour by studying the linearized flow near the fixed points.
The usual FRGE arguments outlined above relies heavily on the mass dimension and scaling with a scale k that possesses units of mass. Such a mass scale is missing in the random tensor setup, instead one wants to study the scaling of the couplings with the dimensionless tensor size N . This requires one to identify (1) the scaling of the IR-suppression term with N and (2) the scaling of the coupling constants with N . It turns out that the requirement that the RHS of the FRGE admits a 1/N -expansion imposes significant restrictions on the scaling with N , but it does not fix it completely. To obtain a completely determined scaling with N one needs to impose that the bare action possesses a geometric interpretation. Essentially, the requirements are that (1) the bare propagator and the modified propagator (after including the IR-suppression term ∆ N S[T ]) possess the same scaling for large index values and (2) that the interaction term possesses the scaling necessary for the geometric interpretation. These two initial conditions, together with the restrictions that result from the 1/Nexpandability of the RHS of the FRGE fix the setting that is sufficient to investigate the large-N -critical behaviour.

2D Causal Matrix Model
A matrix model that enforces a preferred time slicing in its Feynman-graphs was proposed by Benedetti and Henson in [1]. This model is constructed using two dynamical N × N matrices, A and B, representing the spacelike and timelike edges of a triangle, and a constant matrix C which implements the dual weighting. The partition function is where in the large N limit the matrix C must satisfy the condition with m ∈ N. The partition function (7) with a weighting matrix C that implements (8) generates Feynman diagrams that possess the geometric interpretation of polytopes with an arbitrary number of space-like edges and only two time-like edges (see fig. 1). This is clear by analyzing the free propagators (g = 0) of the model As we can see in fig. 1, this restriction implements a foliation of the discrete geometries that appear in the expansion of the partition function, thus introducing the structure necessary for the implementation CDT in tensor models. We proceed by integrating out matrix B since it is a gaussian integral, obtaining This partition function together with condition (8) determines our starting point in this work. One can understand the integration over matrix B as the gluing of triangles along their spacelike edges. This gives rise to a model of squares only with timelike edges. This produces an an-isotropic quadrangulation with rigidity associated with condition (8). We identify the matrix model action that implements CDT as which takes the form of the Euclidean action [19], except for the presence of matrix C. The appearance of the dual weighting matrix C changes the symmetry of the matrix model. Let us consider an N × N orthogonal matrix, O, and the transformation Since the combination AA is invariant, (14) is a symmetry of the Euclidean and the CDT action, however the conjugate symmetry under A → O A and A → A O is only a symmetry of the Euclidean action and not of the CDT action. This shows an explicit difference with the real Euclidean matrix model. In the language of the renormalization group: the CDT action (13) lives in a different theory space, which is governed by a different symmetry.

Weighting matrix
The matrix C implements the weighting of closed loops of propagators in the Feynman graph expansion, i.e. the dual weighting of the Feynman graphs. In principle one could define this matrix abstractly only through the property (8) and only use equation (8) whenever the matrix occurs in a Feynman diagram. However, in order to do practical calculations with the FRGE, it is very useful to have an explicit representation of C at ones disposal. For a N × N diagonal matrix, X, with eigenvalues {x i } we can write its characteristic polynomial as where e k is e 0 (x 1 , ..., x N ) = 1, then Newton identities allow us to write this coefficients in terms of the k-th power of the trace of X, p k , in the form so, C can be found by solving These solutions exist by the fundamental theorem of algebra and one can use ones preferred approximation scheme to obtain these. One scheme that suggests itself in particular when one wants to gain insights into the effects of dual weightings is to build a matrix C from smaller blocks of matrices C o , so C = diag(C o , ..., C o ). The matrix obtained in this way does not implement the entire tower of equations (8), but permits traces periodically. This approach is particularly interesting to study, since it allows us to study how many of the equations one needs to enforce to attain the phase transition between the Euclidean matrix model and the CDT matrix model. The first three matrices C o are • k = 4: • k = 6: Putting these together as blocks to build an N × N matrix gives for example for k = 2 and for k = 4 that are N ×N matrices formed by 2×2 and 4×4 blocks, and where γ = 1.02 and ξ = 0.70.

Functional Renormalization Group for Matrix Models
Let us briefly review the application of the FRGE to matrix and tensor models. One can follow the fundamental presentation of [18] and apply it to matrix models as done in [21]. The starting point is the definition of the effective average action Γ N [φ] in the presence of an IRsuppression term ∆S N [φ]: where where ϕ represents the expectation value of a quantum field φ, while the term represents an IR-suppression term in so far as it is designed to give a mass term of order N to "IR" degrees of freedom of the matrix. Since matrix and tensor models do not implement a fundamental scale, there is no canonical identification of which degrees of freedom are "IR". Rather one needs to implement by hand a division of theory space according to an RG scale k. The simplest assignment is to identify the upper-left corner of the matrix with index values below the scale k as "IR". Once the IR-suppression term is chosen, one can proceed as in [18]; one arrives at the FRGE where t = ln N . The solutions to (26) are functionals of the N ×N matrix φ and hence of infinitely many degrees of freedom in the large N -limit. Practically one resorts to finite truncations of the effective average action, i.e. one performs an expansion of the effective average action into monomials Then one truncates this expansion at a manageable set of operators O i . In this way one reduces he study of the flow to a projected flow in the space of coupling constantsḡ i . The quality of the FRGE results depends critically on the operators that are included in the truncation. In matrix models it turned out that surprisingly good approximations to the FRGE flow where obtained in [21] by considering the flow of single trace operators. The analogous truncation in the presence of the matrix C is which only includes operators with AA , which is invariant under (14), and two C matrices. In the present contribution we will truncate this to the ansatz that contains the bare action and the single trace operator that can directly contribute to the beta functions of the bare action at one loop. This truncation is: We introduce the dimensionless couplings where α 4 and α 6 are as of yet undetermined, since the matrix model does not include an intrinsic notion of scale. The scale is later fixed by imposing that the beta functions admit a 1/N expansion. To make the calculation concrete, we choose the explicit form of the IR-suppression term R N to take the form which has the advantage of being a diagonal and field independent tensor, so we can readily invert the kinetic term to obtain the propagator. It is practically useful to split the second variation of the effective average action into a field independent term G and a field dependent term F : which allows us to expand the RHS of the Wetterich equation as a geometric series, using only the propagator P = G −1 and the F -term: The upshot of this P − F expansion is that each F term contributes more field operators. Hence a truncation that contains polynomial operators with only a finite number of fields terminates the geometric series at a finite number of terms. With the proposed truncation (29), the first term in (33) is of order zero in the Feynman diagrams expansion since there is no field contribution, the second term gives rise to 2-vertex and 4-vertex diagrams which contribute to η and β 4 , the third one to 4-vertex and 6-vertex diagrams which contribute to β 4 and β 6 and so on.

Benchmark results
We use the FRGE to find fixed points of the RG flow and to investigate the universality class associated with this fixed point. This is done by calculating the critical exponents θ at the fixed point, i.e. by considering the linearized FRGE-flow at the fixed point, where the critical exponents appear as the eigenvalues of the Hessian of the beta functions. We chose our truncation in such a way that we can resolve the fixed point that known as the double scaling limit in the matrix model literature. This fixed point possesses a single positive critical exponent, which is usually expressed in terms of the string susceptibility γ str : For Euclidean Matrix Models [19] γ str = − 1 2 , while for CDT [20] γ str = + 1 2 , which leads to the following critical exponents

β-functions
In this section we summarize the steps that we took to obtain the beta functions of the matrix model for CDT.

Operator products
The structure of the beta functions is determined by the operator products of the F -terms that we showed in (33). The terms F (4) N and F (6) N are the second variations of the operators whose contribution to the effective average action is measured by the coupling constantsḡ 4 andḡ 6 . The second variations take the form By looking at (33), we notice that traces of products of (36) and (37) give rise to a big range of operators which are not present in the truncation ansatz (29), such as (T r(CA)) 2 , T r(CAA C 2 )T r(A CA), etc. However, a reasonable projection rule onto the truncation should not project these operators onto the beta functions of the truncation. We therefore analyze which operators can contribute to the beta functions in the truncation, i.e. to β 4 and β 6 . Considering for example the trace T r F (38) we see that each term contains the matrix C, while we know from the structure of the P-F-expansion that these are the only terms generated by the restriction of the FRGE to the truncation that contain two matrices A. Hence, the restriction of the FRGE to the truncation does not generate terms ∼ T r(AA ) and hence does not generate any term that contributes to the anomalous dimension η. To generate a contribution to the anomalous dimension, one needed to include a term with a single C matrix in the truncation. This term would then be generated at one loop by the first term in (38) and in turn contribute to η at one loop. The investigation of this kind of secondary effect however goes beyond the scope of this first investigation.
This analysis relies on the fact that our projection rule is able to discern the structure in which the matrices A are contracted with the constant weighting matrix C, so at first sight one might worry that such a projection does not exist. However, one can consider the appearance of the matrix C in the operators as a special case of operators with index-dependence, i.e. operators whose variations w.r.t. A can not be expressed in terms of A and δ ij , which can be discerned by a suitable projection rule. Hence it is not only possible, but even prudent to use a projection rule that discerns the different ways in which the matrix C is contracted.
To make this distinction, we mark in the following the terms that contribute to the beta functions in our truncation by putting a box around them. Subsequently, we will impose the use of a projection rule that only retains these operators and thus consider only the contributions of the boxed terms.
T r F T r F +T r(CAA CCAA CCAA C)T r(δ).
When considering T r (F N ) n with n > 2 we see that all resulting operators contain at least three C matrices which are not present in the original proposed action (29), this means that in this truncation the βfunctions do not possess contributions coming from these traces.

General Form of the β-functions
Now that we have identified the terms that can contribute to the β-functions, we can write down the general structure of the beta functions. To do so, we introduce the constants D i , E i and F i , which depend on the details of the projection rule. Repeating the same argument as in the previous subsection for the single trace truncation 28 we obtain η = 0, for i odd for i even We can see in particular that in this truncation tadpoles and 2-vertex diagrams contribute.

Analytic fixed point analysis
By setting our truncation to (29), we obtain the following beta functions where η has been set to zero in accordance with our previous analysis. The set of fixed points of this system of beta functions is and the Hessian matrix (defined as H ij = − ∂βi ∂gj ) is Hence, the critical exponents, i.e. the eigenvalues of the Hessian, evaluated in each of the fixed points take the form Critical point Critical exponents Table 1: Critical points with its corresponding pair of critical exponents.
The canonical dimensions α 4 and α 6 are not fixed by themselves, but tadpole diagrams show that α 6 has to be one dimension of N greater than α 4 . We identify the Gaussian fixed point (0, 0), which will not have a relevant direction. There are two non-Gaussian fixed points: the first fixed point α4 E2 , 0 contains one relevant and one irrelevant direction if E2 E3 > α4 α6 , and the third fixed point in Table 1 possesses a relevant and an irrelevant direction if E2 E3 < α4 α6 . These two points are our candidates for a double scaling limit. Next we will examine them using particular schemes.

Scheme dependence
To find concrete critical exponents, we supplement the projection rule with the evaluation of both sides of the FRGE at preferred test matrices A. Moreover, we consider the two rigidity matrices obtained by constructing a block diagonal matrix from (18) and (19), namely (21) and (22) respectively. The specific test matrices that we use for the projection are Using three different matrices A allows us to estimate a lower bound for the scheme dependence. We expect this because the three field configurations contain two field configurations with distinct UV behaviour and one manifest IR field configuration. One often assumes that the scheme dependence is an actual approximate measure for the quality of the fixed point analysis, however when comparing with analytic results, we will see that this underestimates the truncation error. The corresponding obtained critical exponents are shown in Table 2, where euclidean values where computed as done in [21] Table 3: Numerical values obtained for critical points. Causal 2 corresponds to values computed using (21) and Causal 4 corresponds to the ones computed using (22).
We see that the relevant critical exponents are all close to 1, while the irrelevant critical exponents spread a bit wider between −0.858 and −1.215. Moreover, we observe that all critical exponents lay within the spread obtained by scheme dependence. This means that we can not distinguish the Euclidean models from the "Causal" models built from the 2 × 2 or 4 × 4 matrices based on the present derivation of the critical exponents.
In order to attempt to obtain more accurate numer-ical values for the critical exponents we use the fixed point approximation. This consists in first finding the zeros of the beta functions, then evaluating the anomalous dimension, η, in the fixed point g * 4 and substituting this numerical value in the beta functions to find the critical exponents. The critical exponents obtained by using the fixed point approximation are shown in the following  Table 4: Numerical values obtained for critical exponents using the fixed point approximation. Causal 2 corresponds to values computed using (21) and Causal 4 corresponds to the ones computed using (22).
Since the numerical values reported in Table 2 were found to have a strong scheme dependence, it is important to compare the renormalization scheme dependence versus the causal-euclidean results in the latter ones in Table 4. We compute the average of the difference between the critical exponents obtained in the different renormalization schemes and the "Causal vs. Euclidean" results with each of both methods.
Renorm. Scheme Causal vs. Euclidean full 0.019 0.012 fixed p. a. 0.049 0.084 We observe that, while the first method (full) shows a stronger renormalization scheme dependence, with the fixed point approximation method the "Causal vs. Euclidean" relation is more significant than the renormalization scheme dependence. Regarding the accuracy of the values for the critical exponent obtained with both methods compared to the theoretical values (35), we observe that the Causal ones differ more from the theoretical value than the Euclidean critical exponents. There-fore we can conclude that in this case the fixed point approximation is more useful for differentiating the Causal from the Euclidean results, while the full method reproduces more accurate numerical results.

Conclusions
This contribution is motivated by the observation that the application of the FRGE to tensor models with dual weights could lead to an approach to quantum gravity that combines the advantages of the systematic search of continuum limits with the FRGE with the physically promising phase diagrams of CDT and EDT. The systematic development of these tools and the systematic investigation of these models is a very ambitious task. In this contribution we took a first step into this direction and considered the FRGE flow of a matrix model for CDT in 1+1 dimensions proposed by Benedetti and Henson. This model implements a foliation through a dual weighting of Feynman graphs, which introduces an action with an index-dependent propagator, which, upon integration of an auxiliary field, introduces an action with an index dependent interaction.
Recalling the critical exponents analysis in Section 4.1 and the values in Table 1, the model uses a rigidity matrix C which is chosen in such a way that one of the conditions for the beta function polynomials E2 E3 > α4 α6 or E2 E3 < α4 α6 is satisfied. In this case we obtain a relevant and an irrelevant directions simultaneously, implying that the only Feynman diagrams that contribute to the partition function are the ones where all C matrices are contracted as T r(C 2 ), which is precisely the condition that implements a foliation. In this contribution we considered a single trace truncation in which we included operators that contain two C-matrices with the pattern prescribed by the interaction in the Benedetti-Henson model and calculated the beta functions for this truncation. We found that wave function renormalization does not occur in this truncation, since the one-loop structure of the FRGE can only remove one C matrix. This technical observation has far reaching consequences for the structure of the beta functions, which change significantly compared to the Euclidean model, which is obtained by setting C to the identity matrix. We then investigated this system of beta functions in a truncation in which we included only a four-and a six-point interaction. Despite the significant difference in the structure of the beta functions, we found that this truncation contains fixed points that possess the properties of the double scaling limit, which we investigated numerically using three distinct field configurations for projection.
This numerical investigation revealed a practical challenge: To obtain numerical values for the beta functions one can not resort to an abstract definition of the rigidity matrix C, but one requires an explicit numerical expression. We took this as an opportunity to investigate the weakening of the condition T r(C m ) = δ m mod k,2 for k = 2, 4, 6. This has the implication that not all Feynman diagrams without a foliation structure are suppressed, but only a part of these. In particular the case k = 2 does not introduce any new restriction at the level of Feynman diagrams of the Euclidean model, however, since we used the structure of the beta function for general C, we still obtained equations that differ from the Euclidean matrix model. The numerical investigations however revealed that we can not discern the Euclidean and the CDT model on the basis of the critical exponents at the fixed point associated with the double scaling limit. These results are summarized in table 2. We see that the relevant critical exponent (θ) in the CDT model differs from the the exact values in 2.5 by 0.29 and 0.24 from the Euclidean one. We also observe a 4% spread of these values depending on the field configuration used for projection, which is a significantly lower spread than the difference with the analytic values. This is however consistent with the results obtained in [21], where a similar difference from the analytic values was found.
We interpret our results as a encouragement for the investigation of dually weighted tensor models for quantum gravity: Already with the rather simple and elementary techniques used in this contribution we were able to investigate qualitatively the double scaling limit of the dually weighted matrix model; analogous dually weighted tensor models such as the one proposed in [22] can thus be treated with the FRGE in a similar fashion. Our results indicate some practical advise for these future investigations: 1. The one-loop structure of the FRGE can only couple effective operators that differ by an indexdependent contraction between two adjacent tensors, not more. Therefore, to study the influence of an operator with index-dependence in more than one contraction one needs to include sufficient "intermediary" operators in the truncation.
2. For analytical investigations it is possible to work with abstract rigidity structures, that are defined through its properties, such as T r(C m ) = δ m,2 , however for numerical investigations one needs to find (or at least approximate) numerical representations of this abstract structure. One might thus prefer the investigations of models for which one has one of these numerical representations at ones disposal.
3. If our present observations about the double scaling limit are transferable then one sees that the existence of a fixed point with certain characteristics can be found in rather small truncations. However, the critical exponents found in these truncations can be expected to differ significantly from the exact values (which of course should be accessible through lager truncations and optimized renormalization schemes).