Dark side of Weyl gravity

With the help of a functional renormalization group, we study the dynamical breakdown of scale invariance in quantum Weyl gravity by starting from the UV fixed point that we assume to be Gaussian. To this end, we resort to two classes of Bach-flat backgrounds, namely maximally symmetric spacetimes and Ricci-flat backgrounds in the improved one-loop scheme. We show that apart from a genuine IR fixed point that is reached at a zero value of the running scale, the renormalization group flow also exhibits bouncing behavior. We demonstrate that the IR fixed point found is IR-stable in the space of the considered couplings. As a next step, we analyze physics in the broken phase. In particular, we show that in the low-energy sector of the broken phase, the theory looks like Starobinsky $f(R)$ gravity with a gravi-cosmological constant that has a negative sign in comparison to the usual matter-induced cosmological constant. We discuss implications for cosmic inflation and highlight a non-trivial relation between Starobinsky's parameter and the gravi-cosmological constant. Salient issues, including the scheme independence of the IR fixed point and the role of trace anomaly, are also discussed.

In the vicinity of the UV FP, we choose a truncation ansatz for the effective action that will be further used to set up the FRG flow equation. Our truncation prescription is directly dictated by the one-loop effective action. We further enhance this by incorporating into the FRG equation two non-perturbative effects, namely threshold phenomena and the effect of graviton anomalous dimension. One may also wonder about how a consistent quantum theory can emerge when the action is problematic at tree-level (ghost problem). Nevertheless, the RG flow analysis reveals that by the inclusion of the two non-perturbative effects the quantum theory yields a sensible IR FP. Indeed, by solving the RG flow equation algebraically for β-functions β C and β GB we show that there exists a non-Gaussian IR fixed point where both β-functions simultaneously disappear. Aforementioned IR FP represents a critical point after which the (global) scale-invariance is broken. This is reflected through the presence of a composite order-parameter field of the Hubbard-Stratonovich type, which in the broken phase acquires a non-trivial vacuum expectation value.
We map the broken-phase effective action on a two-field hybrid inflationary model that in its low-energy phase approaches the Starobinsky f (R) model with a non-trivial gravi-cosmological constant. Requirement that Einstein's R term in the low energy actions must have a coupling constant As for the conformal anomaly, typical imprints of it are new terms generated in the functional-integral action via higher loop corrections. In particular, loop corrections will generate the Weyl-symmetry violating R 2 term in the action. There are various possible ways how one can deal with conformal anomaly. For instance, one might embed QWG in N = 4 conformal supergravity found by   [39,40] (which is known to be unique anomaly-free theory containing Weyl gravity in its bosonic sector) or in Witten-Berkovits twistor superstring theory [41].
One should then consider QWG as the low-energy limit of such UV-finite models. In these cases, QWG is confined to the purely bosonic sector of spin-2 and spin-1 fluctuations where the conformal symmetry is only softly violated. Another option is to take seriously also theories with anomalous conformal symmetry similarly as done, for instance, in string theory with strings propagating in non-critical spacetime dimensions [42]. Here we will employ yet another scenario, in particular we will take advantage of the fact that the symmetry-breaking R 2 term is presumably generated only at the two-loop order (so that to one loop we do not observe the appearance of the ensuing anomalous conformal mode). On the other hand, our FRG approach will show that the IR FP (for the two involved couplings) appears already at the enhanced one-loop level and hence the prospective observational consequences of the trace anomaly do not take over before the Weyl symmetry is dynamically broken. In the broken phase, the trace anomaly is not anymore a signal of inconsistency since the theory there is not scale-invariant to begin with and hence there is no reason why the ensuing energy momentum tensor should be traceless.
The word dark in the title of this paper refers to two different things. First, it is related to conceptual and technical issues that plague QWG and that make it discouraging or dark in eyes of many practitioners. Second, appearance of the dynamical gravi-cosmological constant in the broken phase of QWG can be associated with ensuing dark energy. It is purpose of this work to demonstrate that despite the aforementioned problems, QWG may serve as a healthy theoretical setup for the UV-model building of phenomenologically viable quantum theory of gravity with pertinent cosmological implications.
Our paper is organized as follows. In the next Section 2 we discuss some fundamentals of both classical and quantum Weyl gravity. In particular, we highlight a formal similarity of the WG with nonabelian Yang-Mills theories and stress some of dissimilarities and potential problems met during quantization. Section 3 is dedicated to the construction of the FRG flow equation for the QWG in the one-loop enhanced scheme, and with this tool we analyze in Section 4 the running of the β-functions associated with the Weyl tensor square and Gauss-Bonnet terms. In particular, we demonstrate that apart from the IR-stable fixed point that is reached at a zero-value of the running scale, the RG flow also exhibits a non-trivial bouncing behavior in the vicinity of the IR fixed point. The issue of conformal anomaly of QWG is discussed in Section 5. There we point out various technical and conceptual issued associated with the conformal anomaly in QWG and propose possible remedies. In Section 6 first employ a Hubbard-Stratonovich (HS) transformation, which introduces a non-dynamical spurion scalar field without spoiling the particle spectrum and (perturbative) renormalizability. After the dynamical breakdown of the Weyl symmetry the HS field acquires a non-trivial vacuum expectation value and gets radiatively generated gradient (kinetic) term. If QWG has any physical relevance, then its effective action in the broken phase must contain an Einstein-Hilbert term. This is shown in the second part of Section 6. We further demonstrate that in the broken phase the corresponding one-loop effective action consists (in the Einstein frame) of two scalar fields -scalaron and dynamical Hubbard-Stratonovich field, that interact via derivative coupling. The resulting low-energy behavior in the broken phase can be identified with Starobinsky's f (R)-model (SM) with a gravi-cosmological constant (dark side of Weyl gravity) that has a negative sign in comparison to the usual matter-induced cosmological constant. After this we discuss anomaly matching conditions between symmetric and broken phase of QWG. A brief summary of results and related discussions are provided in Section 7.

Classical Weyl gravity
The WG is a pure metric theory that is invariant not only under the action of the diffeomorphism group, but also under Weyl rescaling of the metric tensor by the local smooth functions Ω(x): g µν (x) → Ω 2 (x)g µν (x). The simplest WG action functional in 4 spacetime dimensions that is both diffeomorphism and Weyl-invariant has the form [10,43,44] where C µνρσ is the Weyl tensor which can be written as with R µνρσ being the Riemann curvature tensor, R µρ = R µνρ ν the Ricci tensor, and R = R µ µ the scalar curvature. Throughout the text we employ the time-like metric signature (+, −, −, −) whenever pseudo-Riemannian (Lorentzian) manifolds are considered. The dimensionless coupling constant α is conventionally chosen so as to mimic the YM action. On the other hand, in order to make connection with the usual RG practice, it will be more convenient to consider the inverse of the coupling α 2 . In the following we denote the coupling in Eq. (1) as ω C , by employing the prescription ω C ≡ 1/(4α 2 ).
Henceforth, we employ the following notation: for the square of the Riemann tensor contracted naturally, that is R µνρσ R µνρσ we use the symbol R 2 µνρσ , the square of the Ricci tensor R µν R µν we denote by simply R 2 µν , the square of the Ricci curvature scalar is always R 2 , while for the Weyl tensor square C µνρσ C µνρσ we employ a shorthand and schematic notation C 2 . When the latter is treated as a local invariant (not under volume integral) in d = 4 dimensions one finds the following expansion of the C 2 invariant into standard invariants quadratic in curvature We will also need another important combination of the quadratic curvature invariants, namely Gauss-Bonnet (GB) term which in d = 4 is the integrand of the Euler-Poincaré invariant [45] In order to underline the similarity with nonabelian YM theories we use the Riemann curvature tensor so that it is related to the usual general relativistic one by Note that due to skew and interchange symmetries of the Rieman tensor we have In addition, also C λµνκ C λµνκ = C λµνκ C λµνκ | genrel , (8) trivially holds. In particular, in our convention we consider both Christoffel symbols and curvature tensors as 4 × 4 matrices; R µνλ κ ≡ {R µν } λ κ and Γ νλ κ ≡ {Γ ν } λ κ . With this one can write which is clearly analogous to the relation for the field strength in nonabelian YM theories with A µ representing the gauge field. A natural geometrical language for the description of this analogy is provided by the fiber-bundle theory.
While Γ νλ κ represents connection in the frame bundle, A a ν plays the same role in the ensuing principal bundle. Similarly, R µνλ κ quantifies anholonomy in a parallel transport around a closed infinitesimal loop in the frame bundle, whereas F a µν quantifies the anholonomy in the principal bundle. Although the fibre-bundle theory provides the most adequate language for the parallelism between Riemannian-geometry-based gravity and nonabelian YM theories we will not pursue this point any further here.
With the help of the Chern-Gauss-Bonnet theorem one can cast the Weyl action S into equivalent form (modulo topological term) It should be stressed that although the omitted topological term is clearly not important on the classical level, it is relevant on the quantum level where summation over distinct topologies should be considered (see following subsection). But even when one stays on topologies with a fixed Euler-Poincaré invariant, the renormalization will inevitably generate (already at one loop) the GB term with a running coupling constant (see, e.g., Section 3).
In passing we note also that both (1) and (11) in general dimension d of spacetime, while |g|G supplies topological invariant only in d = 4. This is particularly important to bear in mind during the quantization where (similarly as in the Yang-Mills theories) one should choose such a regularization method that preserves the local gauge symmetry of the underlying Lagrangian and thus does not introduce any unwanted symmetry breaking terms.
For this reason one should preferentially rely on fixed-dimension renormalization and avoid, e.g., dimensional regularization. This is the strategy we will pursue also in this paper.
Variation of S with respect to the metric yields the field equation of motion (EOM) known as the Bach vacuum equation: where B µν is the Bach tensor (trace-free tensor of rank 2) and "; α" denotes the usual covariant derivative (with Levi-Civita connection). We remind that the form (13) of the EOM is specific only to d = 4. Moreover, apart from being traceless, the Bach tensor is also divergence-free (B µν ; µ = 0), which ia a consequence of diffeomorphism symmetry. Because Bach tensor results from variational derivative of the action S with respect to symmetric metric tensor g µν , it must also be symmetric When for a given background B µν = 0 (i.e., given backgound is a classical vacuum solution of Weyl gravity), then we say that it is Bach-flat. More general discussion of classical singularity-free solutions in WG can be found in [46][47][48].

Quantization
One can formally quantize WG by emulating the strategy known from quantum field theory, i.e.
by introducing a functional integral (h = c = 1) Here Dg µν denotes the functional-integral measure whose proper treatment involves the Fortunately, for our purposes, i.e., for the computation of enhanced one-loop effective action the use of quantum partition function (14) will suffice.
The sum in (14) is a sum over four-topologies, that is, the sum over topologically distinct manifolds M i (analogue to the sum over genera in string theory or sum over homotopically inequivalent vacua in the Yang-Mills theory) which can potentially contain topological phase factors, e.g., the Euler-Poincaré characteristic of M i , cf. Refs. [53]. It should be stressed that the sum over four-topologies is a problematic concept since four-manifolds are generally un-classifiable -i.e., there is no algorithm that can determine whether two arbitrary four-manifolds are homeomorphic.
On the other hand, simply connected compact topological four-manifolds are classifiable in terms of Casson handles [54], which can be applied in functional integrals in Euclidean gravity. For simplicity, we will further assume that all global topological effects can be ignored, so, in particular, we assume that our space M i is compact and its tangent bundle is topologically trivial.
To avoid issues related to renormalization of non-physical sectors (i.e., Faddeev-Popov ghosts and longitudinal components of the metric field) it will be convenient in our forthcoming reasonings to employ the York decomposition of the metric fluctuations h µν defined as where we have denoted the background metric as g µν . The York decomposition is then implemented in two steps [38]. In the first step we rewrite the metric fluctuations as where h is a trace part of h µν andh µν is the corresponding traceless part. More specifically, We will always tacitly assume that the Lorentz indices are raised or lowered via the background metric, i.e. via g (0)µν or g (0)µν , respectively. Also all covariant derivatives ∇ µ below will be understood as taken with respect to the background metric. In the following the operator will denote the so-called Bochner Laplacian operator [38], i.e., the covariant operator defined as In the second step, we decompose the traceless part into the transverse, traceless tensorh ⊥ µν and to parts carrying the longitudinal (i.e., unphysical) degrees of freedom, namelȳ These mixed-longitudinal (and traceless) parts are written in terms of an arbitrary transverse vector field η ⊥ µ and a scalar (trace) degree of freedom σ. The last fields must satisfy the usual conditions of transversality and tracelessness, i.e.
The true propagating degrees of field in QWG are associated with the transverse and traceless field h ⊥ µν ≡ h TT µν . Indeed, from the second variation of the Weyl action expanded around a generic background it can be seen thath ⊥ µν is the only field component that propagates on quantum level. The vector field η ⊥ µ and two scalar fields h and σ completely drop out from the expansion due to diffeomorphism and conformal invariance respectively [55]. Some explicit examples will be given in Section 4.

Exact RG flow for quantum Weyl gravity
For convenience sake, our subsequent reasonings will be done in an d = 4 Euclidean space dimensions as this is a typical framework in which the FRG treatment is done. By performing Wick rotation from Minkowski space to Euclidean space the question of the resulting metric signature arises. When one does, in a standard way, only the change of the time coordinate t → −it E , (where t E is the name of the first coordinate in the Euclidean characterization of space) the resulting signature of the metric of space is completely negative, that is (−, −, −, −). It seems natural to define the corresponding GR-covariant d'Alembert operator as E = −∇ µ ∇ µ , where the generalization to curved Euclidean space is done by using Bochner Laplacian. However, in all formulas that follow, we find more convenient to use the following definition in the Euclidean signature = ∇ µ ∇ µ . We also remark that this last operator , if analyzed on the flat space background has negative semi-definite spectrum. We will also use a definition of the covariant Euclidean box operator (covariant Laplacian) ∆ = = − E and this last operator in the Euclidean flat space case has a spectrum which is characterized by −k 2 , the d = 4 Euclidean negative square of a 4-momentum vector k µ . Accordingly, the signature of the metric in Euclidean space will be taken to be (+, +, +, +).
The aim of this section is to explore the IR behavior of the QWG by starting from the presumed UV FP where the QWG is exact. Existence of such a UV FP was self-consistently checked in Ref. [16].
To this end we will solve the FRG flow equation [35,36] for the effective average action Γ k , which reads The IR-cutoff R k suppresses the contribution of modes with small eigenvalues of the covariant Laplacian −∆ ≪ k 2 , while the factor ∂ t R k removes contributions from large eigenvalues of −∆ ≫ k 2 .
In this way the loop integrals are both IR-and UV-finite [56]. The second variational derivative of the effective action -Γ (2) , depends on the background metric g µν , which is the argument of the running effective action Γ k , while k is the running energy (momentum) scale or the momentum of a mode in the Fourier space. We also employ the notational convention ∂ t = k∂ k .
Ideally, Eq. (20) would require calculation of the full resummed and RG-invariant effective action.
It is, however, difficult to proceed analytically in this way so we content ourselves here with the conventional procedure, according to which one should employ some well motivated ansatz for the effective action. In particular, in order to evaluate the RHS of Eq. (20) we employ the (Euclidean) effective action in the enhanced one-loop scheme. By the enhanced one-loop scheme we mean one-loop effective action in which also effects of the anomalous dimension and threshold phenomena are included. Corresponding truncation will thus inevitably go beyond the usual polynomial ansatz.
On the other hand, for the LHS of (20) we project the flow on the subspace of the three invariants containing precisely four derivatives of the metric (Eqs. (3), (4) and R 2 invariant). The reason why we consider effective action on the RHS being different from the effective action on the LHS, is dictated by technical convenience. Namely, the RHS acts as source for the RG flow, while the LHS contains the desired structure of the effective action that is appropriate for the extraction of the β-functions.
Let us now briefly describe the basic steps that are used in solving the FRG flow equation.
Further technical details as well as necessary derivations can be found in our recent paper [16] and in its supplemental material [57]. On the other hand, our final results will be discussed in more detail in the following section.
First step is the construction of the one-loop partition function for QWG. It is important here to choose convenient class of backgrounds. We selected maximally symmetric spaces (MSS) and Ricci-flat backgrounds, knowing that both of them are also Bach-flat, so they are classical exact solutions to the Bach equation (i.e., vacuum equations in WG). The partition function we construct for physical degrees of freedom (only transverse and traceless gravitons) and for this we use York decomposition outlined above. We take care of the change of variables under the functional integral, addition of the Faddeev-Popov determinant, inclusion of the differential Jacobian, and also on exclusion of the contribution of zero modes, which are unwanted. We also check that all these one-loop partition functions provide 6 propagating degrees of freedom in QWG.
In the second step, we analyze the general β-functional of the theory considered on our general backgrounds (MSS and Ricci-flat). We ask and address the question about information which can be extracted on these backgrounds about the β-functions of involved couplings. Due to various relations between curvature tensors and between square of tensors, we find that we are able to extract only two combinations of couplings: either β R + 1 6 β GB on MSS, or β C + β GB on Ricci-flat background. It can be also checked that the usage of general Einstein spaces as on-shell backgrounds does not improve on this situation. However, in a general theory to the quadratic order in curvatures we may set up an ansatz for the effective average action Γ built out of three couplings ω C , ω GB , and ω R corresponding to three quadratic invariants C 2 term, G = GB term, and R 2 term. We address this potential inconsistency of the RG system (three couplings and only two extractable β-functions) in the final step of our method of solving the FRG.
The most important part consists of building and solving the FRG flow equation. In Ref. [16] we found a novel form of the FRG flow equation based on the expressions for factors appearing in the quantum partition function of the theory. These factors mimic the simple scalar two-derivative kinetic sectors of the theory, however they may appear both in the numerator and in the denominator of the partition function and with various mass square parameters (which can also be negative). We take into account quantum wave-function renormalization of quantum fields and add to the flow equation the anomalous dimension η for all fields participating in the quantum dynamics at the one-loop level.
To do the IR-suppression of modes in each factor we should add a suitable cutoff kernel function R k (z). When all these operations are done the final form of the FRG flow equation, reads where Y i,j are mass square parameters of the modes φ i and the functional trace is done over space of the same modes, which may also contain some traces over internal indices of the fields. The frontal plus/minus sign in (21) originates from the initial position of the factor in the partition function (whether it was in the denominator, or in the numerator of the partition function, respectively).
The final step consists of choosing the truncation ansatz for the action, whose FRG flow we try to determine. To write the RHS of the flow equation we explicitly reflect the fact that our theory is QWG theory, where we could add a topological GB term, therefore not influencing at all the perturbative one-loop partition function. On the LHS we could analyze the flow of the general action quadratic in curvatures with the structure of the Lagrangian ω C C 2 + ω GB G + ω R R 2 . However, in perturbative QWG at the one-loop level there is a very interesting and special simplification since at this quantum level the R 2 term is not needed, and we know that perturbatively its one-loop β-function β R vanishes. We utilize this hierarchy of β-functions in QWG and employ consistently the truncation ansatz without the problematic R 2 term. Therefore we have two combinations of β-functions ( 1 6 β GB and β C + β GB ) and two couplings ω C and ω GB , so the system is consistent and possible to be solved algebraically for the β-functions involved.
This solving we do in the last technical step, where we evaluate all the functional traces (both in internal space and spacetime volume integrals). We perform the traces using the heat kernel method and Barvinsky-Vilkovisky trace technology. In order to have an analytic control over all formulas, we decide to choose a particular form of the cutoff kernel function R k ≡ R k (z) = k 2 − z θ k 2 − z due to Litim [58]. The main point of using the Wetterich equation (20) is to take into account massive modes which slow down the RG flow in the IR regime. We achieve this by adding cutoff kernels R k in mass-dependent renormalization scheme of Wilsonian character. This lets us to obtain one-loop RG-improved expressions for the two β-functions of the theory with all quantum effects due to anomalous dimension η and IR threshold phenomena included.

Analysis of β-functions and RG fixed points
We now discuss the system of β-functions of the theory considered at the one-loop level improved by the usage of FRG methods. Based on the computation presented in Ref. [16] the explicit form of β-functions reads and with the anomalous dimension of the graviton field given by to the one-loop level of accuracy. By ω C = ω C (k) we here denote a running coupling parameter in front of the C 2 term in the action (1) -the so-called Weyl coupling.
The above two β-functions, β C and β GB , follow from the FRG with the truncation ansatz motivated by the one-loop level and with both the threshold phenomena and non-trivial anomalous dimension of the quantum graviton included. We call them thus as being one-loop RG-improved. The effects of threshold phenomena are present explicitly only in the expression (22) for β GB . However, due to the combination in (23), the solution for β C will also inherit these threshold factors. Finally, we observe that the anomalous dimension η enters only multiplicatively in the system of β-functions (22)- (23). This has some simplifying consequences for a search for FP's of the coupled system, both in the UV as well as in the IR regimes.
Let us now discuss the reasons for the presence of threshold phenomena in our system. As it could be seen from the expressions for the one-loop partition functions of the system [16] on MSS, and on Ricci-flat background We can further simplify the system of equations (22)-(23) for the two β-functions. In particular, we do not wish to solve explicitly the system (22)- (23). We just concentrate on the corresponding FP's. This is a much simpler task as we can solve the system of β-functions algebraically. This gives where with b = 137/60 and y = 1/ω C . The origin of the common denominator 1 + y(X − b) is entirely due to the inclusion of the graviton's anomalous dimension η. When η is neglected, the latter is unity.
This is the regime, in which the Weyl coupling ω C is big (ω C → ∞, so y → 0). This corresponds to the perturbative regime of the theory (in terms of α C ). When one decides to neglect these common denominators, one gets simplified expressions, which are already sufficient to shed light on the issue of existence and character of FP's of the FRG flow. The equations in (29) still include the effects of threshold phenomena. When we neglect the threshold phenomena in our description, then the system of β-functions acquires the form of one-loop perturbative system as derived in [49] for QWG in dimensional regularization scheme. Actually, all these threshold phenomena are contained in the expression called X above. When one takes the limit Λ/k 2 to zero, then all threshold factors are indeed removed, and the expression X reduces to just a number β FT GB = −87/20.

Ultraviolet Asymptotic Freedom in all couplings
When κ = k/ |Λ| ≫ 1, the threshold phenomena are completely irrelevant and can be neglected, cf. Eq. (22). Irrespectively of the initial values of the couplings ω C0 = ω C (t 0 ) and ω GB0 = ω GB (t 0 ) (look also at an analysis in the next paragraph), the leading RG running behavior in the UV regime (for t ≫ 1) is ω C ∼ tβ FT C and ω GB ∼ tβ FT GB . This signifies that the absolute values of the couplings must necessarily grow in the UV. It might be argued that in the UV regime one can also neglect the effects of the anomalous dimension η, cf. Eq. (24), since it is suppressed by big values of the ω C coupling in the UV. For the UV-running (t → +∞ equivalent to k ≫ k 0 ), it suffices to use only the non-RG-improved one-loop perturbation results (29) from above.
Fradkin and Tseytlin [49] were the first to find that the one-loop β-functions for ω-couplings are constants with values This leads to an asymptotic freedom at the UV FP for all two couplings as we shall prove below. Since β C > 0, it is natural to assume that the initial condition of the flow is such that ω C (t 0 ) > 0 and similarly since β GB < 0, then ω GB (t 0 ) < 0. If one chooses the opposite condition, then the RG flow tends to decrease the absolute value of the ω-coupling, the coupling crosses zero, and finally it goes on the other side, where the initial conditions are natural in a sense mentioned above. This is because one-loop RG flow in the UV forces the β-functions to be constants, so the increments of the couplings (positive for ω C and negative for ω GB ) are regular and linear in the UV regime. Hence, in the deep UV (close to the UV FP) we can assume that ω C > 0 and also that ω GB < 0.
All these arguments are self-consistent and lead to the conclusion that the UV fixed point of RG inevitably exists and realizes the asymptotic freedom (AF) scenario (in much the same way as in non-Abelian gauge theories). Our perturbation analysis is carried out in terms of the coupling g 2 ∝ 1/ω, so this and the fact that in the UV regime, ω C → +∞, bolsters even more the correctness of our perturbative one-loop results. Actually, near the UV Gaussian FP, it is the coupling g (analogous to the YM coupling constant) that goes to zero.
The asymptotic freedom (AF) characterization of the FP comes when the RG flow is analyzed in terms of g-like couplings. We define them, taking into account above signs of ω C0 and ω GB0 , in the following way: and with the inverse relations In this way we are sure that both g-couplings are non-negative. Now, one can easily derive by and (mind the sign in the second equation). Then the β-functions of these couplings are expressed as which to one-loop accuracy read For vanishing values of the g-couplings we find that the above β-functions vanish too, so in this situation we have a trivial Gaussian FP of the RG flow. We also see that for the positive values of the couplings, that is g C > 0 and g GB > 0, we have that both β-functions β g C and β g GB are negative which signifies that in the UV we meet a FP (asymptotically free theory) in basically the same way like this happens in QCD.

Scheme independence of IR fixed point
In the paper [16], the detailed analysis of the situation with FP's in the IR was presented. We found both: a turning point of the RG flow at some finite energy scale k and also a true IR FP at k = 0, so in the deep IR regime. The genuine IR FP has a non-trivial (non-Gaussian) character and it lets to define the theory in the non-perturbative way, which is free from any IR-type of divergences and it therefore realizes Weinberg asymptotic safety, but in the infrared. The turning point (TP) of the RG flow is related to the non-analytic behavior of the β-functions in terms of running coupling parameters. This in turn corresponds to a very interesting cosmological bounce scenario, when analyzed from the AdS/CFT correspondence point of view. In this subsection, we discuss some universal features of the TP of RG and of its location, while in the next subsection 4.3, we prove the relation between the non-analytic behavior of β-functions and the bounce happening at finite k.
One can easily see that qualitative features of the RG flow presented above do not depend on the details of the regularization and renormalization procedures. For example, taking a closer look at the plot (cf. Fig.1. from [16]) of the dependence of the β-functions on the energy scale, one can convince oneself that the asymptotics in the UV limit of the flow and an existence of the location of the first zero of the RG flow counting from UV direction are universal. Firstly, we would like to discuss the UV regime of the RG flow. The asymptotics in the UV limit is described by one-loop β-functions computed by Fradkin and Tseytlin for 4-dimensional conformal gravity and our exact (or rather FRG-improved) β-functions tend asymptotically to these constant values (if the β-functions for ω-type of couplings are considered). (Here we also suppressed the contribution from the anomalous dimension η since we know that this goes as some inverse powers of the ω C coupling and should be accurately included only when higher-loop computation accuracy is required. We have η ≪ 1 and effectively we can take η = 0, which is a judicious assumption in the UV regime of the flow towards asymptotically free point.) As a matter of fact, we see that the β-function of the ω C is positive (with the precise value β C = 199 30 ), while the β-function of the coupling ω GB is negative (with the precise value β GB = − 87 20 ). Since in this regime, we reach an asymptotically free UV FP, where the g-couplings go to zero (and correspondingly ω-couplings tend to infinite values), any RG flow towards such a FP should coalesce with the perturbative one-loop RG flow for small values of the g-couplings. We also see this as a feature of our improved RG flows. Therefore, the locations of horizontal asymptotes of the RG flow from the right (so for k → +∞) are universal and scheme-independent.
In the paper [16], we studied the energy evolution of the running couplings towards the IR regime of the flow, concentrating on the IR fixed points of the system. The inclusion of threshold phenomena, which are present in any mass-dependent renormalization, is of crucial importance in our analysis. In fact, if we have studied only the simplified system of β-functions (29), we would not find any interesting behavior of the RG flow in the infrared (similarly to the case of QCD in the IR regime where the coupling grows stronger and gets out of the perturbative regime). We do not find any IR FP in such a simplified scheme. The β-functions from the system (29) for X = const are always constant, at any energy scale. To search for some non-trivial behavior in the IR, we must thus include some additional non-perturbative effects. This feature is brought about by our usage of FRG methods and account of decoupling of massive modes in the IR domain.
In order to find the FP's of the system in the IR regime, we must solve equations β C (k) = 0 and β GB (k) = 0. One can see a big simplification here because in order to find zeros, we do not need to solve the full system (27). Actually, we can completely forget the denominators in (27) and solve only Eqs. (29), where threshold effects embodied in the factors Λ/k 2 are still taken into account. We also notice that the anomalous dimension η does not influence the locations of the possible IR FP's, within the limits implied by our truncation ansatz used in FRG.
Numerical solutions of the equations β C (κ) = 0 and β GB (κ) = 0 reveal that they are both satisfied at (approximately) simultaneous values of the rescaled energy scale κ = k/ |Λ|. This is a smoking gun for the fixed point of the RG flow. These zeros are automatically zeros of the exact system (27). Moreover, the location of the zeros is almost identical (up to 2% accuracy) for the couplings ω C and ω GB for both cases of Λ > 0 and Λ < 0. The inclusion of higher-loop effects or extension of our truncation ansatz will make this agreement even stronger, so that in an exact fully non-perturbative theory, the locations of two zeros coalesce into the one unique location of a genuine FP for both couplings, in the infrared regime.
As we remarked above, the non-trivial form of the running arises because we have included the effects of threshold phenomena. Let us, for definiteness, analyze closer the case of the parameter Λ > 0. It is straightforward to understand the behavior of the exact β-functions β C and β GB regarding their zeros and ensuing non-trivial FP's in the IR domain. For this, it is important to find the behavior of the MSS one-loop partition function treated as a rational function of the energy scale k. The first pole/zero of the partition function (when one is coming from large values of k) we find at k 2 = 4 3 Λ, and this is due to the factor det 2T ˆ − 4 3 Λ1I , which is present in the denominator of the partition function (25). This latter implies that this factor appears in the FRG flow equation with the positive coefficient, because it was a pole (not a zero) of the partition function. As it was explained earlier, for the contribution of this factor to the Wetterich equation, we need to evaluate the following functional trace Tr 2T which in the RG flow is with the coefficient +1. For details of the structure of the FRG flow equation the reader is referred to [16] and the formula (47) there. Here, we meet the moment where we see the dependence on the cut-off kernel function R k = R k (z), so the dependence on the scheme of renormalization. We will now show that for the qualitative features of the flow important for the existence of non-trivial FP this dependence is immaterial. First, the trace in (37) leads to the expression in terms of the heat-kernel B 4 expansion coefficient of the operator, , again with positive front coefficient. To get the front factor of the above expression we used the optimized (Litim) form of the cutoff kernel, [58], and we will remark on other possible choices below. One can check that the B 4 coefficient of such an operator is also positive. Indeed, we find, which is valid on MSS background and under the space volume integral (to produce the integrated B 4 coefficient). This implies that in the flow of the action ∂ t Γ L k (as evaluated on MSS background), we have the term 1 − 4 3 Λk −2 −1 3 5 Λ 2 , again with the positive sign. Hence, the expression for the β-function β GB , which is the only one that can be read from the situation on MSS background, contains the factor 9 40 1 − 4 3 Λk −2 −1 , when we again emphasize the positive sign.
By solving algebraically the linear system of the β-functions (29), we also derive that in consequence the β-function β C contains in turn the factor 9 40  Now, one can look at this factor from the broader perspective. It describes the decoupling (due to threshold phenomena) of the massive modes with the mass square parameter on MSS background given by 4 3 Λ. The coefficient in front is related to the B 4 coefficient of the corresponding operator and that is why it is positive. In a general renormalization scheme, this coefficient is also positive and the threshold factor 1 − 4 3 Λk −2 −1 which shows the pole, precisely at k 2 = 4 3 Λ, must be present in this form or in another more general one, but still in a form exhibiting the pole in the same place. This is due to the gauge-invariant and universal fact that on MSS background in Weyl conformal gravity we find stable (non-tachyonic) 2T modes (spin-2 traceless) with the mass square given by 4 3 Λ. And this is the highest mass square parameter in the spectrum of all modes there. In a general framework, the general threshold factor could look like this where a is a positive constant and f (x) is some smooth, regular and positive function at x 0. For the β-function β C this threshold factor appears as with a constant b < 0. This form manifests all universal features that we have discussed above. Now, it is a matter of simple analysis of functions, that for the GB term coupling ω GB , if a > 0 and β GB → − 87 20 < 0, when k → +∞, then the running β-function β GB (k) must meet a zero for some In the latter the "+" superscript signifies that we approach the respective value from the above. Similarly, for the Weyl term coupling ω C , if b < 0 and β C → 199 30 > 0, when k → +∞, then the running β-function β C (k) must meet a zero for some 4 3 These are the invariant features of the flow and that is why the IR FP (for k = k GB for the ω GB coupling and for k = k C for the ω C coupling) is universally present in any renormalization scheme that aims at properly taking into account threshold phenomena. It is obvious that any mass-dependent scheme comes with its own form of the function f = f (x), positive for x 0, to supply the correct decoupling of heavy, massive modes in the IR regime. (As a such renormalization scheme we cannot, for example, select a DIMREG because it is mass-independent.) As we have argued, the form of the function f (x) is irrelevant and the IR TP is scheme-independent. Mathematically, the zeros of threshold factors which appear in the expressions for the FRG-improved β-functions are the results of the continuity of the function f (x) and of the existence of the first vertical asymptote at k 2 = 4 3 Λ, and of the way how it is approached by β-functions. It is also of crucial importance that the factor k 2 − 4 3 Λ first appeared in the denominator of the partition function (25), so the sign was favorable to enforce the change of the sign of the β GB , when moving from negative universal one-loop value in the UV towards lower energies.
This happens also because the vertical asymptote is approached from the right (higher energies) to +∞. The change of the sign of β GB must occur between two regimes: UV (when β GB < 0), and IR near the first vertical asymptote (when β GB → +∞), so the zero line must be crossed for some energy scale k = k GB . A similar conclusion holds for the β C with inherited threshold phenomena obtained from the system (23). Moreover, it is natural to expect that k GB ≈ k C since the difference can be associated only to higher-loop accuracy error as one can easily see by comparing the numerical values, which are quoted below for Λ > 0: Finally, one sees that the values of the energy characterizing the IR TP are very close to the lower bound given by k = This means that we must inevitably find a FP in the infrared regime. Therefore, the existence of the IR TP is a universal feature of the exact (improved) RG flow, while the details of its location, slopes, and speeds of approaching the FP, etc., depend on the particular choice of the renormalization scheme (or in FRG terminology on the choice of the cut-off kernel function R k (z)).

Non-analyticity near turning point
Let us now observe that the system of β-functions Eqs. (27) and (28)  Here, we present the analysis near the turning point of the RG flow, where β = 0, which occurs at t = t * and with the value of the coupling which is ω = ω * (this value of the coupling ω * depends on the initial condition of the flow, that is on ω 0 = ω(t 0 ) and t 0 , while t * is independent of them).
We analyze the general situation for one representative coupling, but we can think of the TP as found in the system of ω C , ω GB couplings for κ * ≈ 1.18. In the linear Taylor approximation, due to the regularity of the β-function at the zero point understood as a function of the t variable only, we can generically write β = A(t − t * ) for a constant coefficient A and this leads to the equation which is solved by with the initial condition of ODE that ω(t * ) = ω * . The above solution shows that there is a minimum of the coupling for the value ω * at t = t * for A > 0 (it is a local maximum when the constant coefficient A < 0). The case A > 0 we meet for the ω C coupling, while for ω GB we have the opposite behavior (due to (23)), so A < 0 there. For definiteness, here we consider the case A > 0.
Inverting the relation in (43), we get that Actually, for the β-function of the GB term, the coefficient A < 0 because the zero of β GB is reached from the other side, so instead of the local minimum for the coupling ω C , here there is a local maximum for the coupling ω GB . Moreover, there is again formally a negative/positive infinite eigenvalue of the RG flow (depending on which branch one is moving on), hence the stability matrix cannot be properly defined. When we consider two couplings at the same time at TP, a question arises which coupling (only one) has to be chosen to locally eliminate the t variable from the system of β-functions. We decided to remove t in favor of the ω C coupling. If the β-function of the GB term, at the common TP (placed conventionally at κ = κ C ), is analyzed as a function of ω C , not of the additional RG-time t variable, and not of the ω GB coupling, then there is a finite positive off-diagonal value of the linearized RG flow matrix because β GB ∼ ω C − ω C * with some finite coefficient. This is due to the facts that κ GB > κ C and that β GB (κ C ) > 0. Here, at the TP, we have a relation that ω C − ω C * is proportional to (t − t * ) 2 with a positive coefficient from (43), hence then the two eigenvalues of the system are zero for the GB coupling and formally positive infinity for the ω C coupling, and these two couplings formally are exact eigenvectors of the matrix of the RG flow. Again, even in the case of two couplings the stability matrix cannot be determined.
The reader can easily see that the analysis presented above hinges on the fact that the location of the TP of the RG flow in RG-time coordinate t * is finite. If formally t * → ±∞, then the initial assumption about the local behavior of the β-function β = A(t − t * ) does not make any sense. Hence, the solution in (43) is not realized in this form. As it is known from the general theory of true FP's of the RG flows, they may only appear at the abstract theoretical RG-scale coordinate t * → ±∞. Following the above distinction between TP's and FP's, in the paper [16], we continued the search for true IR FP's. For this purpose, we used the established fact that at κ ≈ 1.18 we found a common TP of the RG system. We exploited the infinitesimal form of the flow at TP and analytically extended it beyond the TP. We treated the TP as a good point from which we could start a new perturbation calculus driven towards the IR regime. Assuming perturbativity (in different couplings than in the UV FP), we were eventually able to find a genuine IR FP at t → −∞. We characterized this FP as non-trivial and non-Gaussian and computed the characteristic values of the couplings there, so called ω C * * and ω GB * * , which revealed to be non-vanishing. Moreover, we found that this IR FP is stable for both perturbation directions given by the couplings ω C and ω GB . The IR-stable true IR FP is the main result to be used for further cosmological and conformal symmetry breakdown related applications of the QWG theory, see also Section 6.

Extension of the RG analysis to two-loop order
In this subsection, we attempt to give a preliminary analysis of the RG system of running coupling parameters in QWG, at the two-loop level. We base our considerations only on algebraic, dimensional-analytic, and combinatorial arguments since a detailed computation of UV-divergences and β-functions at this level is still beyond our computational capabilities. We assume that numbers, we are dealing below with, are generic and they do not vanish, and we discuss the general structure of the RG system. In particular, we touched upon the issue of the "new" β-function β R , which is expected to be generated first time at the two-loop level. We analyze its suppression compared to other β-functions in the system and establish the hierarchy of them. Moreover, we also look at the universality properties of the β-functions for all three couplings ω C , ω GB and new ω R at this level and in this way we strengthten and extend the well known results from the one-loop quantum level to QWG theory at two loops.
As remarked, in the paper [16], the first paper to deal with the divergences issue at the two-loop level in QWG, was the one by Fradkin and Tseytlin from 1984 [59]. However, the computation presented there is not complete, since only a subset of two-loop diagrams is analyzed. Nevertheless, we agree with the authors' conclusion that it is very probable that the R 2 divergence shows up for the first time at the two-loop accuracy. This is in contradiction to the conjecture of 't Hooft and Mannheim [44], who instead expect that the conformal symmetry on the quantum level is so powerful that this non-conformal β-function β R is vanishing to all orders and also non-perturbatively. This would be true, if the conformality was fully present at the quantum level (not only at the one-loop level, where it forces β R = 0). We do not think to be so, in accordance with [40,59], because of the presence of conformal anomaly, already at the one-loop perturbative level. More issues related to the conformal anomaly we discuss in the special discussion section 5.

Two-loop suppression of β-functions
Let us remind that the one-loop action for QWG reads This also served us as the truncation ansatz for the effective action Γ that we used for the FRG computation. We re-emphasize that the term R 2 with the coupling ω R 2 is consistently not included at the tree-level and in the one-loop motivated RG flow equation, because such a term is not generated by any quantum correction at the one-loop level. In the original one-loop computation by Fradkin and Tseytlin [49], the quantities which are assumed to be small are 1 where the coefficients β FT C , β FT GB are simple numbers (30), while up to the two-loop order we must find β (2) and β (2) where the numerical coefficients a (2) C,C , a C,GB , a GB,C , and a GB,GB are presently unknown, but it is certain that they do not depend on the couplings ω C , ω GB . The two-loop form of the RG system presented above is the result of assuming the perturbative expansion in 1 ω C and 1 ω GB variables.
The UV-divergent part of the effective action at the one-loop level is given schematically by where we find that to one-loop accuracy we have due to (partial) conformal symmetry still preserved at the quantum one-loop level. This fact can be viewed as the one-loop remnant of full conformal symmetry present in the action at the tree-level (44). At the two-loop level we expect β (2) R not to vanish and be given analogously by where the coefficients a R,GB are presently unknown numbers, whose non-vanishing (even of one of them), if unambiguously computed, would completely prove the conjecture of [59]. We explain that the possible structure term a (2) R,R ω −1 R is not present since only ω C and ω GB are the couplings in the original (and also one-loop level) action. The coupling ω R has to be introduced (and renormalized) only from the two-loop level only. Actually, the reason for its introduction at the two-loop level action is the presence of the R 2 counterterm in Γ (2) . We need to absorb such a covariant UV-divergent term and for this we need to include the ω R R 2 term in the bare action. This also means that for perturbative computations at the level of three loops and higher, we must use the bare action (44) corrected by the presence of this new term ω R R 2 with arbitrary coefficient ω R (however, due to hierarchy and suppression of β-functions, as explained below, we should assume that its value is parametrically smaller than the values of other couplings ω C and ω GB present in (44), that is we should use ω R ≪ ω C , ω GB ). For full one-and two-loop level quantum computations we can use the bare action as given in (44) and this is reflected in the results for the β-functions to this accuracy given in (46), (47), and (50). It is conceivable that if one wants to theoretically go to three-loop expressions for the β-functions of any of the coupling X, Y = ω C , ω GB , or ω R , then the terms with the structure a The two-loop level UV-divergent part of the effective action takes then the form so the new term R 2 is generated with the coefficientβ R , which is always suppressed by one power of the small coupling 1 ω C or 1 ω GB , as in (50). Compared to the one-loop level action (48), where the counterterms were multiplied by only numerical coefficients β FT C and β FT GB this is a suppression by additional power of small coupling. We conclude here, that the β-function β R when it finally shows up at the two-loop level is additionally suppressed with respect to other β-functions in the RG system. This signifies that the hierarchy of the β-functions is evident and the running of the ω R coupling is very small, and that it was fully consistent to assume to the one-loop accuracy that ω R = 0. This was the fact that we took advantage of in the truncation ansatz for Γ that we used to model FRG to the one-loop level. Sincerely speaking, the significant R 2 term could be generated in the truncation ansatz Γ, but this does not happen immediately, and it requires a long RG-time since the running of ω R is very slow. Our truncation ansatz for QWG is therefore internally consistent, at least in a big vicinity of the UV FP of RG. We just remark that we took care of the fact that at the two-and higher-loop level coefficients of UV-divergences are not the same as higher-loop β-functions of couplings (but they are in strict relations) and therefore we decorate the terms in (51) by additional primes.
If one uses the electric charge-like couplings defined for any coupling ω (from the set ω C , ω GB , and ω R ) (cf., also the analysis in Subsection 4.1) by then the corresponding β-function reads Consequently, we find that to the one-loop level accuracy β (1) and hence β (1) g GB ∝ g 3 GB , and finally β (1) g R = 0. Similarly, based on Eqs. (46) and (52), up to the two-loop level we find which implies that We repeat verbatim for the GB term coupling: and for the ω R coupling (remembering that β (1) We see that at the leading order the β-function β g R is proportional to the fifth power of the g-couplings.
Again, compared to the expressions for β g C and β g GB , which to the leading order (which is a one-loop order) go like g 3 , this is a two-loop suppression. We derive that this suppression is present independently of the form of the couplings used in QWG theory.
When using electric-like-couplings g i , we are in a comfortable situation that perturbative calculus is conveniently done in these couplings (compare to a different case for ω-like couplings).
For example, we write the one-loop effective action Γ (1) ∝ O(g 0 ) as a perturbation in small couplings g i to the classical action S ∝ O(g −2 ), and then the genuine corrections at the two-loop level are The coupling of the R 2 term in the effective action Γ for QWG is generated from the level of (at least) two loops and its first coefficient is proportional to g 2 , so compared to other terms like C 2 or E present at the one-loop level and multiplied by the coefficients of order g 0 , this coefficient is again highly suppressed. We also draw a parallel that analogously like in QWG considered above, in QCD, the β-functions scale like: β (1) g ∼ g 3 for the one-loop level and β (2) g ∼ g 5 for the two-loop level, and similarly the tree-level action scales as S ∝ O(g −2 ), the one-loop level effective action scales as Γ (1) ∝ O(g 0 ), and the two-loop effective action scales as Γ (2) ∝ O(g 2 ), where g is the Yang-Mills coupling parameter.

Universality of two-loop order β-functions in Weyl square gravity
Let us now present a general argument for the universal properties of the RG system of β-functions of all three couplings ω C , ω GB , and ω R in QWG at the second loop level. Our exposition is loosely based on Ref. [60]. First, we recall that to the one-loop accuracy all β-functions of perturbative couplings are universal and do not depend on the choice of the renormalization scheme, and are also gauge-fixing independent, if the corresponding couplings are in front of dimension four gauge-invariant terms in the action. Of course, the form of the β-functions depend on the version of the couplings used (whether this is an ω-type coupling, or g-like coupling, or something else) as clearly seen in the formulas from the last subsection. Here, we analyze the change of the β-functions under a general redefinition of the couplings. Such redefinition mimics the change done by using a different regularization or renormalization scheme, or the change of the gauge-fixing conditions.
Provided that this change is continuously connected to the point g i = 0 and that it is analytic in the couplings g i , we find that the type of the coupling is preserved. (We remark that the change between ω-type and g-type couplings is not regular at the point g = 0, hence it does not satisfy the above criterion, and as can be clearly seen from expressions in Eqs. (30) and (36), the corresponding β-functions are different, though they are obviously related to each other.) For the analysis in QWG, we also use the fact that ω GB coupling cannot appear in any Feynman diagram computation, since it multiplies the Gauss-Bonnet term, which is a total derivative in d = 4. Hence, there cannot be any dependence on ω GB in any perturbative β-function of the theory. This observation simplifies also the analysis of the previous subsection, where we can effectively set everywhere 1/ω GB → 0.
The β-function of the R 2 term is expected first to show up at the two-loop level and in the form where a ≡ a (2) R,C is some constant coefficient, which has not yet been computed. Effectively we have here the system of two couplings ω R and ω C and three β-functions. (In accordance with the remark made above, we do not have ω GB as the coupling on which these β-functions could depend.) First, the β-function of the coupling ω C is in its sector universal (it is behaving effectively like the only coupling in the town). By looking at the formula in Eq. (46) when we set a (2) C,GB to zero, we see that this RG sector is identical with a two-loop sector of a theory with only one coupling ω C . As it is known from [60], two-loop β-function for any QFT system of only one unique coupling is also invariant under the general coupling redefinition transformations. Now, the question is only about the β-function β R . We prove below that its two-loop value is also universal. In the two-couplings system (ω R , ω C ) we are allowed to change couplings only in the following way, written to the leading order ω 0 , But there exists a requirement that in any gauge choice or parametrization method, the β-functions β ′ R and β ′ C , when expressed in terms of primed couplings, cannot depend on ω ′ R , because there are no Feynman diagrams depending on this coupling (here on ω ′ R ) needed to be considered to the two-loop level of accuracy. (The bare action to the two-loop level contains only one term ω C C 2 .) Hence, one can check, that the only possible changes of the couplings are shifts according to: and then, of course, we have that the transformed β-functions read to the quadratic order in the coupling ω ′ C . This proves the universality in the ω R and ω C sector of the expressions for the two-loop β-functions of these couplings of the theory.
Similarly, we find for the Gauss-Bonnet coupling ω GB , which does not appear in any perturbative Feynman rule of the theory, so none of the β-function can depend on it. Hence, the only permissible changes of couplings are given by Then, according to (47), the β-function is to two-loop order accuracy β (2) and a transformed β-function is so again there is a universality of this two-loop expression for β (2) GB . Basically, all the couplings in QWG at the two-loop level, behave like if they were in separate one-coupling sectors of the couplings' space of the theory. Everything boils down to the fact that the RG system of β-functions at two-loop level depends only on one coupling ω C . This is a special feature of QWG, that to the two-loop level we need only to deal with one coupling of the Weyl square term, the GB term will not have any impact perturbatively and only from the third loop we need to include the new coupling ω R of the R 2 term.
This new term ω R R 2 in the bare action is heralded by the presence of trace anomaly already at the one-loop level, to which discussion we turn now.

Discussion of the trace anomaly issue
Let us now discuss the issue of the trace (or conformal) anomaly. Firstly,we recall the following fact about trace anomaly in ordinary Yang-Mills gauge theories in d = 4 spacetime dimensions.
Standard YM theories are described by the action S g = − 1 2 d 4 x tr(F 2 µν ). This entails that on the classical level the theory is conformally invariant and hence the trace anomaly on this level vanishes.
One could compute the trace of the classical energy-momentum tensor (obtained by Hilbert method of variation with respect to some fiducial metric tensor of any curved background) and then it is found to vanish in agreement with conformal symmetry. On the other hand, if the theory is not very special, then for a generic situation at the quantum (loop) level there is a non-vanishing β-function β of the YM coupling, which is a signal of the presence of trace anomaly. Due to the RG-invariance in the effective action of the YM model at the one-loop level we have the term where D 2 is the square of the gauge-covariant derivative operator and µ is this renormalization scale which was also used to put renormalization conditions for fields and compensate µ−dependence of the running dimensionless YM coupling. This RG running is in an invariant manner (and universal to one-loop level) described by the non-vanishing β-function β. One could ask for the conformal properties of the newly generated term (71)  Simply saying, quantum physics of standard pure gauge theories is more complicated and more interesting than just what was on the tree-level constrained by scale-invariance. The amplitudes are more complicated and show more intricate behavior with energy scales as the result of freedom from conformal symmetry. The trace anomaly is not here a problem and can be used to generate some terms in the effective action. Actually, it signals that quantum gauge field theory are more interesting and worth studying. Here conformal symmetry was never meant to be used as local gauge symmetry hence its lack after inclusion of quantum corrections (due to polarization effects of gluons) is not problematic. We could say that conformal symmetry was an accidental symmetry of the tree-level (classical) generic YM theory and at the quantum level we have seen that it is not there anymore.
Since we do not put much of emphasis to this symmetry in gauge field theories, then its loss it is not a big deal like this happens with other accidental symmetries. We can sacrifice easily this conformal invariance and we shall not regret it (although now the computations on the quantum level are much more involved). This is the physics of YM theories, and still it is fully consistent without conformal invariance on the quantum level.
There is a completely different attitude for quantum conformal gravity theories, since there by definition we want to use the conformal symmetry with its fundamental and not accidental role. This fundamental role is signified in that we want to use this symmetry to constrain the form of all possible terms in the gravitational action (only C 2 action is acceptable in d = 4 dimensions), to constrain the spectrum of the model, which then is different from a generic spectrum of four-derivative gravitational theory. These things shall not be understood as accidents due to conformal symmetry, they are definite virtues. And the fact that they happen is not surprising but rather demanded and expected from the fundamental role of conformal symmetry in constraining gravitational interactions.
We want to use full conformal group as the gauge symmetry group of conformal-gravitational interactions. By gauging full 15-parameter conformal group we give the decisive and fundamental role to both the local Poincaré and local conformal symmetries. Only together they could be used to put a control over quantum theory of gravitational interactions.
Let us now return to the issue of the trace (or conformal) anomaly. For the energy-momentum tensor of the total system (matter + gravity) we define the trace in the classical and quantum cases respectively as For quantum conformally invariant theories we should find that g µνT µν = 0. By explicit computation on the classical level for conformal models we find that T = 0, but on the quantum level one can find T = 0, namely by general arguments [61] T = cC 2 + aG + a P P + a L R , where P = ǫ µνλζ R µν ρσ R ρσλζ is the parity-odd Pontryagin density, which can be a priori excluded in parity-preserving theories (as QWG is), and the last term is ambiguous as it depends on the renormalization prescription. Coefficients a and c are the central charges that can be directly related to the corresponding β-functions studied in the previous subsections (namely c and a are gauge-dependent for spin 3/2 and 2 but c + a is not [62,63]). One can find a straightforward relation between the logarithmically divergent part of the quantum effective action and the anomalous terms above. In particular, at one loop one can check that so that c = β (1) C and a = β (1) GB . In our case this anomaly poses a problem: the classical theory (before quantization) was with local conformal symmetry, while quantum theory is without it. Quantum effects break the gauge symmetry of the model on the full quantum level (e.g., conformal analogues of Slavnov-Taylor identities are not satisfied). Consequently, the symmetry is without any power to constrain UV-divergences, Green functions, the form of the quantum corrections to the effective action or scattering amplitudes. So, in order to get a better grasp on the trace anomaly issue it is imperative to use the effective action Γ rather than semiclassical arguments.
First, in pure C 2 gravity we note that the β-function of Weyl coupling β C is non-vanishing at the quantum level (starting from first loop). Hence in the effective action Γ fin we must have a term with µ being an arbitrary renormalization scale. This is due to RG-invariance of the total effective action Γ and the β-function β C . In conclusion the total action in d = 4 is Relevant observation pertaining to it is that Γ tot is not conformally invariant action in d = 4! Definitely, conformal symmetry is not present on the quantum level in Weyl square gravity.
There is the perturbative trace anomaly issue, but the problem is also at non-perturbative level. This means that the theory despite being perturbatively renormalizable is however, non-perturbatively non-renormalizable [40]. This is because the radiative correction break Weyl symmetry in a dramatic way. Pure C 2 gravity is anomalous. In other words, there is RG running (towards IR) as explained in this work, but the problem is at high energies (UV-limit) for overall consistency of the model of conformal gravity. It is very important to cancel this anomaly since conformal symmetry appears in a local (gauged) version. Similarly, like in gauge theory models we have to ascertain that there are no gauge anomalies due to fermions in the matter sector since the presence of such anomaly would be disastrous for the gauge symmetry on the quantum level. Theory would be inconsistent again because of lack of symmetry on the quantum level. In a perturbative vein there would be pop up perturbative UV-divergences that we will not be able to absorb in any form of counterterms possible in gauge-invariant actions. These counterterms would not be gauge-invariant since on the quantum level there is no any symmetry constraining their form. And if there is not a constraining (or forbidding) symmetry law on the quantum level in QFT everything what is possible to be generated is generated (according to the QFT's Murphy's law).
The possible trace anomaly resolution is given below. It is known that there exists special matter content coupled to C 2 gravity such that β C = 0 (and all other β-functions in the matter sector vanish too) and then the terms C log C in the effective action Γ are not generated at all. For example, in the case of N = 4 supergravity Fradkin and Tseytlin [59] showed this can be achieved by coupling it to N = 4 super-YM (SYM) gauge field theory model, which is known to be conformal on flat spacetime. This is brought by example of coupling N = 4 supergravity due to Fradkin and  From a different perspective, if there is no conformal (gauged) anomaly on the quantum level, then the theory is endowed with quantum scale-invariance, there are no non-zero β-functions, no RG flow and the theory sits at the UV FP of RG. This is a situation that has to be secured in UV by embedding the quantum Weyl gravity theory in the bigger picture, e.g., in the N = 4 supergravity or in a twistor superstring theory of Berkovits and Witten. Another possibility is that the theory reaches in the UV regime a non-trivial (non-Gaussian) FP of RG. With the additional conditions that the dimension of the critical surface on which putative UV FP lies, is finite, this realizes the Weinberg's Asymptotic Safety scenario (AS). In these circumstances, UV-divergence problem is avoided and there are no β-functions (even on the non-perturbative level) since the theory sits at the FP. There is no RG flow and scale-invariance of the situation is enhanced to full conformal invariance and the UV-action that describes such a FP is an action of conformal field theory (CFT). There is a whole discipline of studies of CFT's but quite a little is known about CFT when the gravity is quantum and dynamical.
But this is a desired CFT, that describes the quantum theory of conformal gravitational interactions at very high energies. As with any other CFT to fully describe it we would need to give the set of primary conformal operators and their anomalous dimensions (conformal weights). This constitutes the set of CFT data. Based on it we could describe any correlation function within this gravitational CFT. To get away from the FP and start RG flow at lower energies some deformation operator must be added to such CFT. That is, we think that by deforming CFT by adding an operator which is not conformally invariant, we start a non-trivial RG flow away from the UV FP. And then following the flow we can reach even the domain of low-energy physics.
As it was pointed out long ago, the pure C 2 theory may reveal to be non-renormalizable from the level of two loops onward. This is due to the presence of Weyl anomaly [40]. Therefore some authors claim that because of the non-zero trace anomaly C 2 gravity does not survive quantization and hence it is inconsistent on the quantum level till the conformal anomaly (CA) is made vanish. Its most direct effect may be the presence of the R 2 term in the divergent part of the effective action. (Such term does not show up at one-loop, but it is expected at the two-loop level because all symmetries that could forbid its presence there are not realized on the quantum level). The UV-divergences which could be covariantly collected in the term R 2 are not absorbable by the counterterms of the original Weyl theory and that is why the theory may show up to be perturbatively non-renormalizable, when a definite answer will be given to the presence or not of the R 2 term on the two-loop level [59].
We emphasized that in our work [16], we worked in the fixed-dimension regularization scheme in d = 4 because this should not produce any spurious (e.g., dimension dependent) Weyl-symmetry violating terms in the effective action. So, we confined ourselves to the cutoff scheme at fixed dimensionality of spacetime. In this way we do not discard or avoid the trace anomaly (which indeed shows up already at the one-loop level), and the statement of its existence is scheme-independent.
Strictly speaking, the onset of the trace anomaly is renormalization prescription independent because the β-function is zero only at the fixed point and the existence of the fixed point is indeed renormalization prescription independent. On the other hand, the actual non-zero value of the trace of the energy-momentum (i.e., trace anomaly) is renormalization prescription dependent, since the β-function is universal only to the second order in the coupling constant. The one-loop value is scheme-independent but from two loops on we expect that such dependence will start to develop. The conformal anomaly is proportional to the β-functions of the model. On different vein, these β-functions can be read off from perturbative UV-divergences of the effective action, hence the problem of trace anomaly is a problem rooted in the UV regime or, in other words, in an UV-completion of the theory.
At the UV FP which is asymptotically free for pure C 2 gravity the conformal anomaly is still there.
The β-functions for omega-type couplings are not vanishing at this FP (couplings multiplying directly invariants like C 2 or G). Only β-functions in alpha-type couplings (similar to electric charge-like couplings) vanish but the anomaly in terms of ω C and ω GB is non-zero even in the UV FP.
Let us summarize here our approach towards the quantization of C 2 theory and RG running towards the IR limit. We start in the UV fixed point where theory is exactly Weyl-invariant (just a single dimensionless coupling) and the renormalization group flow deforms the theory (i.e., induces other terms) as the theory flows towards IR fixed point. It turns out that the theory has a non-trivial (non-Gaussian) fixed point in IR and hence it is asymptotically safe in IR or, in another words, it is non-perturbatively renormalizable in IR. This asymptotic safety (AS) is in accordance with Weinberg definition [64], but this time for the IR, rather than UV sector. Our situation with the issue of trace anomaly is slightly different since we have asymptotic safety in IR. So we do not assume that theory is valid at all energy scales down to IR. It very quickly flows to the IR FP that is so close to the UV FP that the R 2 term does not even have enough (RG-)time to appear in the effective action (it does not appear there at one-loop level), and after IR FP is reached we do not have anymore scale invariance.
In other words, we surpass the problem of perturbative non-renormalizability by non-perturbative renormalizability in the IR regime.
One may wonder whether following RG trajectory towards lower energy one ends up at a reasonable point (be it Gaussian-like or Banks-Zaks FP) or diverge. In the latter case, the theory would need some kind of IR-completion or protection against infrared problems and IR-divergences.
For this, if we have AS in IR, then the theory is non-perturbatively renormalizable and thus solves all such problems near IR FP.
As already mentioned, the problems related to conformal anomaly are essentially related to the UV regime. One might thus guess that having AS in IR FP does not a priori helps in solving UV problems. For this, we need to provide evidence for a separate UV FP. And it is obvious that looking locally in the parameter space the existence of IR FP does not imply anything for the existence of any UV FP. One possible scenario that can allow to infer possible existence of UV FP for QWG is to embed the theory in a broader context of more fundamental (a better behaved) theories like twistor string theory or N = 4 conformal supergravity. In particular, that the bosonic C 2 sector that we consider can be understand as a truncation of the anomaly-free N = 4 supergravity of Fradkin and Tseytlin [49]. The super N = 4 case of Fradkin and Tsetylin was later seen to provide the exceptional solution to the issue of trace anomaly [49] (however, it seems that there is still a unitarity problem for this model). It is the latter theory which is in the UV without any problem, there is no conformal anomaly, no non-trivial β-functions. This theory is perturbatively renormalizable and UV-finite and hence there are no any problem with quantization. In our previous work [16] we consider the IR version of QWG when we reduce and consistently integrate out degrees of freedom and we study only the bosonic spin-2 sector. Or in a different perspective, QWG appears as the low energy limit of twistor string theory where there are not apparent problems with conformal anomaly. Simply the conformal anomaly is a problem of the UV regime but we assume that somehow this problem is solved there and we start with the theory, which is consistently reduced, in the intermediate energies.
We assume that such UV FP exists and somehow addresses the anomaly issue, and then we try to see if such a postulate is consistent. At the same time there is a subtle difference with what some people typically mean by UV problem of conformal gravity. We believe that most practitioners just start perturbatively from IR FP (hence assume Gaussian or Banks-Zaks IR FP) and deduce the problem in UV by increasing energies in the RG flows. Natural question arise, what about if one cannot even start perturbatively from IR as in our case (and also, e.g., in QCD). What then one can say about UV FP?
We offer a plausible scenario that is logically consistent and reasonably well motivated. We assume that UV FP is just a full-fledged critical point in a series of hypothetical phase transitions that the Universe has undergone in its very early stage. It might be FP that descends from the spontaneous symmetry breaking (SSB) phase of twistor string theory or N = 4 conformal supergravity. In any case SSB FP has exact conformal symmetry and the simplest theory in d = 4 that lives at that critical manifold is QWG (only one coupling ω C ). Then we may hypothesize about UV-completion in a form of string theory or something similar but clearly we do not have to do this now, if we have some strong evidences for UV FP and for solving some UV problems there. Moreover, if someone would like to see the RG evolution from IR FP towards UV, then he would probably need to go through horrible technical problems that would be even more complicated than in QCD. We are more like physicists, who start with asympotically free (deconfined) QCD and then deduce non-perturbative (confining) IR phase. This way is far simpler than going vice versa.
It should be stressed that when using FRG we do not a priori assume that there is any IR FP.
This is a clear advantage of FRG approach. We just chose the truncation ansatz, which in our case is motivated by perturbative one-loop results. Should we have known the higher-loop (or even better the full non-perturbative) effective action then we could utilize its form and set up more informed truncation ansatz, which in turn could provide more reliable (or more refined) results.

Physics beyond IR critical point
In this section we study the physics of the 2nd order phase transition that is associated with the IR FP of QWG discussed in Section 4. The characteristic (in fact defining) feature of 2nd order phase transitions is existence of the order parameter field that has zero vacuum expectation value (VEV) in the symmetric phase and non-zero VEV in the broken phase. Typical method used in this context is theory of effective potentials. This approach will be presented in the following two subsection. In Subsection 6.3 we will see that one can harness the conformal anomaly to say something more on the structure of an actual critical point and corresponding phases in its vicinity.

Composite Hubbard-Stratonovich field as an order-parameter field, emergent Starobinsky's model
Because we require that our theory should induce Einstein action in the low-energy limit we rewrite the R 2 -part in the Bach action (11) with the help of the Hubbard-Stratonovich transformation [65][66][67] as It is clear that the HS transformation (77) is nothing but a simple identity based on a functional Gaussian integral. Although the auxiliary HS field σ(x) does not have a bare kinetic term, one might expect that due to loop corrections the renormalized action will develop in the IR regime a gradient term, which then allows to identify the HS boson with a bona fide propagating mode.
This mechanism is well-known from condensed matter theory [76,77] and particle physics [78][79][80][81]. A typical example is obtained when the BCS superconductivity is reduced to its low-energy effective level. There the HS boson coincides with the disordered field whose dynamics is described via the Ginzburg-Landau equation [68].
The σ field in (77) can be separated into a background fieldσ ≡ σ corresponding to VEV of σ plus fluctuations δσ. Sinceσ is dimensionful it must be zero in the case when the theory is (globally) Weyl-invariant. On the other hand, when the Weyl symmetry is broken,σ develops a non-zero value.
So, the σ field plays the role of the order-parameter field. With the benefit of hindsight we further introduce an arbitrary (hyperbolic) mixing angle ϑ ∈ R and make the following splitting If we now apply the HS transformation only to the S R 2 sinh 2 ϑ part of the action (11), we obtain and hence the full action (1) can be rewritten (modulo topological term) as By employing the Weyl symmetry one can formally generate a gradient term for the σ field already on the level of a bare-action. In fact, if we perform the Weyl rescaling of the form g µν = |σ| −1g µν we get the gradient term − |g| 3 2 (dots refer to higher order derivative terms of the σ field). It should be, however, clear that such a gradient term cannot define genuine propagating mode since it depends on the conformal scaling and hence it is a gauge dependent concept. This is actually consistent with the fact that the WG has no propagating scalar degree of freedom on the bare-action level. On the other hand, when the scale symmetry is broken then the σ-field will be trapped in a particular broken phase with specific kinetic as well as potential term, cf. Eq. (87).
Although the full theory described by the action S is manifestly independent of the mixing angle ϑ, truncation of the perturbation series at a finite loop order will inevitably destroy this independence.
The optimal result can be is reached via the principle of minimal sensitivity [82][83][84], which posits that if a perturbation theory depends on some unphysical parameter (e.g., sinh 2 ϑ in our case) the best result is achieved if each perturbation order has the weakest possible dependence on the parameter ϑ.
As a result, the value of ϑ is determined at each loop order from the vanishing of the corresponding derivative of effective action.
Let us now see that the VEVσ indeed acquires a non-zero value at low enough energies. Since the corresponding phase transition must be of 2nd order (σ must change smoothly from zero in symmetric phase to non-zero in broken phase) this is tantamount to showing that the IR FP must exist and that it is associated with the (global) scale-symmetry breakdown. To this end we employ the standard effective-action methodology in the mean field approximation. In particular, we replace in (80)  As discussed in [65], the one loop-effective potential acquires the form and this result is true both in dimensional [65] and ζ-function [85] renormalization scheme. Here µ is the renormalization scale.
The corresponding VEVσ is obtained through minimizing V eff which gives σ(ϑ) = 2ω C µ 2 exp   3 sinh 2 ϑ log 3 4 sinh 2 ϑ + 1 + 16π 2 ω C (4 sinh 2 ϑ + 1) sinh 2 ϑ(32 sinh 4 ϑ + 16 sinh 2 ϑ − 1) We might note that the valueσ = 0 is not a local minimum for V eff when sinh 2 ϑ > ( because in such a case V eff < 0, while forσ = 0 one has V eff = 0. In order to see whether situation withσ = 0 can be realized we employ the principle of minimal sensitivity, namely we require that This equation admits two branches of real solutions [65]. The branch that corresponds to the symmetric phase (σ = 0) satisfies the equation sinh 2 ϑ = 0.02592337 − 0.0000197α 2 + O(α 2 ). The broken-phase branch (σ = 0) corresponds to large values of sinh 2 ϑ with actual value depending on ω C . Consequently, we can rewrite to order O(1/ sinh 4 ϑ) Eq. (83) as 38 of 48 In particular, for any value of the dimensionless coupling ω C , we can choose the renormalization scale µ, in such a way thatσ ∼ 1/κ 2 , which in turn will guarantee phenomenologically correct gravitational forces at long distances. Consequently, Newton's constant κ 2 is dynamically generated. Owing to the last term in (79)- (80), the appearance of a cosmological constant in the low-energy limit of the broken phase is also consequence of QWG. In the following it will be convenient to rescale σ → σ/κ 2 so that σ ∼ 1 and the field itself is dimensionless.
As already mentioned in Section 2.2, the local scale symmetry dictates that the scalar degree of freedom must decouple from the on-shell spectrum of QWG. When the conformal symmetry is broken the scalar field reappears in action through a radiatively induced gradient term of the HS field σ. The explicit form of the kinetic term can be decided from the momentum-dependent part of the σ-field self-energy Σ σ . In Ref. [65] it was shown that the corresponding leading order gradient term is of the form 1 2κ 2σ ∂ µ σ∂ µ σ. By assuming that in the broken phase a cosmologically relevant metric is that of the Friedmann-Robertson-Walker (FRW), then, modulo a topological term, the additional constraint holds due to a conformal flatness of the FRW metric [86]. It was argued in [65] that from (11) and (79) one obtains in the broken phase the steepest descent (or WKB) gravitational action of the form In the long-wave limit one can neglect fluctuations of the σ field and consider that σ is basically described by its VEV, i.e., σ =σ. This yields in the broken-phase the mean-field effective action By comparing (88) with (80) we obtain that The action (88) is nothing but the Starobinsky action with the cosmological constant. In the Starobinsky model (SM) the linear Einstein term determines the long-wavelength behavior while the R 2 -term dominates short distances and drives inflation, which is followed by the gravitational reheating with the decrease of R 2 [88]. In phenomenological cosmology, the SM represents metric gravity with a curvature-driven inflation. In particular, it does not contain any fundamental scalar field that could play the role of an inflaton field, even though a scalar field/inflaton formally appears when transforming the SM to the Einstein frame [87]. We stress that the cosmological constant Λ is entirely of a geometric origin (it descends from the QWG) and it enters with the opposite sign in comparison with the usual matter-sector induced (i.e., de Sitter) cosmological constant. In the following we will call such gravitation-induced cosmological constant as gravi-cosmological constant and it represents what we call the dark side of QWG.
In this connection the following point deserves mentioning. Since the conformal symmetry prohibits the existence of a (scale-full) cosmological constant, the gravi-cosmological constant must correspond to a scale at which the conformal symmetry breaks, which in turn determines the cut-off scale of the σ-field. The magnitude of ξ in the SM is closely linked to the scale of inflation [40]. Using the values relevant for the Cosmic Microwave Background radiation (CMB) with 50 − 60 e-foldings, the Planck data [74,89] require that ξ ∼ 10 −13 GeV −1 or equivalently ξ/κ ∼ 10 5 . The vacuum energy density is thus ρ Λ ≡ Λ/κ 2 ∼ 10 −10 (10 18 GeV) 4 , which corresponds to a zero-point energy density of a scalaron with an ultraviolet cut-off at about 10 15 − 10 16 GeV. This coincides with a range of the GUT inflationary scale. For compatibility with an inflationary-induced large structure formation the conformal symmetry should be broken before (or during) inflation. We will further discuss this issue in the following section.

Broken phase QWG and hybrid inflation
By now it is well recognized that the scalar cosmological perturbations are (nearly) scale-invariant with the value of the spectral index n s = 0.97, which is tantalizingly close to 1 -exact scale invariance [65,75]. The combination of data from Planck and BICEP2/Keck Array BK15 [89] tightens the upper bound on the tensor-to-scalar ratio to obtain: r < 0.056 at 95% CL. These include, e.g., the SM, the non-minimally coupled model (∝ φ 2 R/2) with a V(φ) ∝ φ 4 -potential, inflation model based on a Higgs field and the so-called universal attractor models [89].
In the previous subsection we have seen that the long-wave limit of QWG approaches in the (∂ ψ ≡ ∂/∂ψ) one can write down the tensor-to-scalar ratio r and the spectral index n s in the slow-role approximation as [89] r = 16ǫ, n s = 1 − 6ǫ + 2η .
In terms of the number N of e-folds left to the end of inflation (ψ f represents the values of the inflaton at the end of inflation, i.e., when e − √ 2/3|ψ| ∼ 1) one gets which for N = 50 ÷ 60 (i.e., values relevant for the CMB) is remarkably consistent with Planck data [74,89].
While during the inflation, the σ-field is constant (due to a large coefficient in front of the gradient term) allowing a large-valued inflaton field descend slowly from potential plateau, inflation ends gradually when σ regains its canonical kinetic term, and a small-valued inflaton field picks up kinetic energy. From (92) the dominant interaction channel at small |ψ| is (∂ µ σ) 2 |ψ|, hence the vacuum energy density stored in the inflaton field is transferred to the σ field via inflaton decay ψ → σ + σ (reheating), possibly preceded by a non-perturbative stage (preheating).
Note also, that the gravi-cosmological constant Λ that was instrumental in setting the inflaton potential in (92) has the opposite sign when compared with ordinary (matter-sector induced) cosmological constant.
Inflationary scheme discussed can be naturally incorporated in a broader theoretical context of "conformal inflation" model, which has been the subject of much recent investigation [90][91][92][93].
Setting the inflation period around the time of Weyl-symmetry breakdown brings about a number of attractive features, including: natural justification for inflation models with a plateau, technically accessible inflationary correlators, specific predictions for the B-mode vorticity fluctuations in CMB powerspectra, etc. Let us also notice that the existence of a single scalar field with cutoff at around the GUT scale and coupled to broken phase of QWG (e.g., our σ or GUT Higgs field) would contribute with a positive zero-point energy that could offset gravi-cosmological Λ and leave behind a small observable cosmological constant. This could provide a viable mechanism for breaking down 120 orders of magnitude difference that presently exists between theoretical predictions and astronomical observations of Λ. Consequently, the dark side of QWG would not be ultimately so dark.

Some reflections on anomaly matching
We have seen that the effective potential approach provides a powerful method for discussing physics in both broken and unbroken phase. Unfortunately, the physics in the immediate vicinity of the critical point is poorly grasped by this approach. The point is that due to the appearance of long-range correlations, perturbative approaches are quite unreliable. It might seem rather surprising that the physics in the vicinity of the critical point might be for CFT's discussed in terms of conformal anomaly. The reason for this is basically twofold. First, conformal anomaly has an universal structure phrased in terms of two central charges (see Section 5). Second, there is a powerful theorem [105] that ensures that there is an anomaly matching across the critical point of the 2nd order phase transition.
In Section 4 we have considered a RG flow of QWG, which ends up in an IR FP that is associated to a CFT. So one could wonder how this IR FP should be associated to low-energy gravity where the scale symmetry is absent. One intriguing possibility is that the IR CFT has a moduli space of vacua, such that in all of them but the one at the "origin" conformal (or scale) invariance is spontaneously broken. Such examples are known, expecially in superconformal field theories [95].
This allows to identify a broken phase of the theory where the analytic structure of the amplitudes is not constrained by conformal invariance, but still the numerical value of the trace anomalies can be possibly matched (as proven, e.g., in Ref. [105] for CFT's with global conformal invariance). This poses strong constraints on the broken phase theory that is in spirit very similarly to what happens in 't Hooft anomaly matching [96]. In particular, requiring for the two cental charges (">" refers to physical energy scales above the FP in question and "<" refers to physical energy scales below the FP 1 ) it was shown to be possible to determine the couplings of the action for the dilaton which emerges in the broken phase as the Goldstone boson related to the spontaneously broken conformal symmetry [105]. Actually, in a supersymmetric scenario where the trace anomaly is in the same supermultiplet as R-symmetry charge this boils down to the usual 't Hoft anomaly matching condition for which the generating functional of amplitudes accounting for anomalies is well understood. This could make it possible to understand the IR regime in the broken phase in a quite novel way. We leave this interesting and important topic for our future investigation.

Conclusions
In this paper, based on the approach of functional renormalization group, we have mainly investigated the IR physics of the quantum Weyl gravity, which turns out to be surprisingly rich and interesting, in particular in connection with the structure of its IR FP and related cosmological implications.
Important point that we have utilised here, was the fact that at one loop the R 2 -infinities of the effective action are absent so that we do not observe the appearance of the anomalous conformal 1 Here we should stress that physical scale is not the same as the running scale in the RG flow. The running RG-flow scale is always re-parameterized so that it reaches its limiting value (±∞) at the FP. In physics the critical point that is represented by a RG FP always happens at some finite energy scale -phase transition scale.
mode. Although it is expected that such situation will change at two loops (contrary to chiral anomalies, which are one-loop exact), one can still view this fact as an indication of a specific feature of the model that could potentially persist even to higher loops, where a conclusive calculations are still missing. In our FRG approach the IR FP (for the two involved couplings) develops already at the one-loop level and hence the prospective observational consequences of the trace anomaly do not take over before the Weyl symmetry is dynamically broken.
Problem that typically besets higher-derivative gravity theories is the unitarity issue. For QWG this is commonly phrased in terms of spin-2 ghost field. As already mentioned in the introduction there exist various remedies to this issue in the literature. If unresolved, the ghost problem prevents QWG from being a possible UV completion of Einstein's general relativity. But this conclusion is based solely on the analysis of a classical Hamiltonian or on the structure of tree-level propagators.
Basically for any physical process that explicitly manifests the ghost problem, there is usually an implicit assumption that the perturbative analysis reflects the true physical spectrum. It is quite instructive to compare QWG with QCD. The latter is also renormalizable and asymptotically free theory but in this case one is accustomed to the fact that the physical spectrum bears no resemblance to the perturbative degrees of freedom. In particular, the gluon is not in the physical spectrum. This is a consequence of confinement, but it is also understood more directly in terms of the behavior of the full gluon propagator. Essentially there is an IR suppression of the propagator that is sufficient to remove the gluon pole. Similarly, we have seen that the FRG analysis of QWG shows that the IR FP is non-Gaussian hence QWG theory is non-perturbatively renormalizable in IR.
The presence of ghosts in the linearized spectrum and the quantum breaking of conformal invariance (trace anomaly) can, however, be seen as an indication that QWG should be understood in the context of a more fundamental theory. One such theory could be the d = 4 twistor string [97,98], where superconformal symmetry in spacetime is explicitly preserved contrary to ordinary string theory. In fact, it was found in Ref. [41] that its spectrum includes to the one of the N = 4 superconformal gravity, which is made up of 2 N = 4 graviton multiplets and 4 gravitini supermultiplets. One of the graviton multiplets is in the non-standard ghost-like sector. The gravitini supermultiplets contain the 15 gauge bosons of the local SU(4) R symmetry. The SU(4) R symmetry is potentially anomalous because its coupling is chiral and, as it is gauged, its cancellation requires the dimension of the gauge group G be 4. So one gets the conclusion that G = SU(2) × U(1) or U(1) 4 . An analogous condition is consistent with, but not implied by worldsheet arguments based on the worldhsheet conformal anomaly. Apart form the SU(4) 3 R , also spacetime conformal anomaly is possible and it has actually been argued that N = 4 superconformal gravity can be made UV-finite and therefore anomaly-free by coupling it to exactly four N = 4 super Maxwell multiplets [59].
The computation was carried out at one-loop order, but the result should hold at all orders as the β-function in N > 1 CSG and conformal anomaly of super YM may receive contributions only from one loop. In the case of pure extended supergravity, it would actually be quite interesting to perform the analysis of functional RG in a situation where a nontrivial IR FP at one-loop order would be an exact result and where also the absence of an R 2 would be granted at any perturbative order.
In the case of N = 4 one-loop exactness is also a consequence of the fact that all these anomalies are related by supersymmetry, and vanish when the SU(4) 3 R anomaly does. As the latter is a chiral one-loop exact anomaly, its cancellation should entail the cancellation of conformal anomalies and other anomalies to all order. The version of N = 4 CSG that was actually found in the context of twistor-string is a "non-minimal" one where the 4-derivative complex scalar φ couples to the Weyl graviton through a term of the kind f (φ) C µνρσ 2 . Imposing a manifest SU(1, 1) ≃ SL(2, R) requires f be a constant so that one can recover the "minimal" version of CSG. It is quite interesting that only recently explicit actions for these theories have been constructed [99][100][101] making a RG flow analysis viable. The anomaly-free theory described above could provide a non-Gaussian UV fixed point for this analysis. We also would like to point out that the role of superconformal symmetry in cosmological applications has found large resonance in recent years [102,103] and in this context the relation between higher derivative terms and emergent order-parameter fields as discussed in subsection 6.1 has been also highlighted (cf., e.g., Ref. [104]) In passing, we note that trace anomalies have also been recently discussed as valuable instruments to follow the violation of conformal symmetry along the RG flow. In particular in four dimensions any such flow can be reinterpreted in terms of a spontaneously broken conformal symmetry. This has led to a deeper understanding of the a-theorem [105,106] and it would be definitely interesting to understand also the RG flow of QWG in this framework.
Apart from these more conceptual issues, we also discussed a cosmology that is implied by the broken phase of QWG. In the one-loop approaximation we have been able to map the broken-phase effective action on a two-field hybrid inflationary model that in its low-energy phase approaches the Starobinsky f (R) model with a gravi-cosmological constant. In particular, the inflationary potential obtained contains two scalar fields that interact via derivative coupling: scalaron ψ and Hubbard-Stratonovich order-parameter field σ (basically dilaton). The scalaron appears when we transform the broken-phase effective action to the Einstein frame and due to its slow-role potential it plays the role of an inflaton. The Hubbard-Stratonovich scalar, on the other hand, results from the dynamical transmutation od the spurion HS field and it mediates the inflationary potential. The derivative interaction between the two scalars provides a viable mechanism for a reheating scenario and graceful exit. Requirement that Einstein's R term in the low energy actions must have a coupling It should be stressed to for simplicity's sake our cosmological considerations were done in flat Minkowski background. In order to get more refined (and also more realistic) information on the inflationary potential we should perform explicit computations of the broken-phase effective potential, in a non-trivial but cosmologically pertinent background, namely on MSS. These backgrounds are tailor-made for treatment of inflation and, in addition, they can provide us with some further guidance for the possible resolution of the cosmological constant problem on de Sitter spacetime. Work along those line is presently in progress.
Author Contributions: All authors jointly discussed, conceived and wrote the manuscript.