A functional RG equation for the c-function

After showing how to prove the integrated c-theorem within the functional RG framework based on the effective average action, we derive an exact RG flow equation for Zamolodchikov's c-function in two dimensions by relating it to the flow of the effective average action. In order to obtain a non-trivial flow for the c-function, we will need to understand the general form of the effective average action away from criticality, where nonlocal invariants, with beta functions as coefficients, must be included in the ansatz to be consistent. We then apply our construction to several examples: exact results, local potential approximation and loop expansion. In each case we construct the relative approximate c-function and find it to be consistent with Zamolodchikov's c-theorem. Finally, we present a relation between the c-function and the (matter induced) beta function of Newton's constant, allowing us to use heat kernel techniques to compute the RG running of the c-function.


JHEP07(2014)040
The renormalization group (RG) underlies most of our modern understanding of quantum and statistical field theories [1,2]. There are different ways to implement the RG procedure. Whereas standard sliding scale arguments (Gell-mann-Low) are particularly suitable for weakly coupled computations, it is only with Wilson's ideas that non-perturbative insights have been possible. The arena where the RG acts is theory space. This space is parametrized by all couplings corresponding to terms which are consistent with symmetries of the systems. The beta functions for the couplings define a vector field in theory space, and the RG flow can be seen in geometrical terms as a certain trajectory in this space.
From this point of view the infrared physics depends upon the differential equation governing the flow as well as on the boundary conditions. If the initial point sits at a finite scale (for example, it is a bare action depending on some UV cutoff Λ), one is considering an effective field theory, whose range of validity is limited by the cutoff scale. However, if we want a theory to be called fundamental, we would like to be able to push the initial scale to arbitrarily high values, eventually to infinity. The only known way to perform this limit is to hit a UV fixed point.
Fixed point theories do not depend on any intrinsic scale since they are scale invariant. As a consequence they can be used to model systems at criticality. These theories are characterized by dimensionless couplings and physical quantities exhibit scaling relations which can be observed in experiments. These relations arise in very different systems sharing the same dimensionality, symmetry and field content. This is what is usually referred to as the concept of universality, the independence of the critical properties of a system from its microscopic details. The RG offers a simple and intuitive explanation of universality: the critical properties of a system are determined by the fixed point, microscopic actions defined at different scales that flow to the same fixed point, or equivalently that belong to the same basin of attraction of a fixed point, will describe the same criticality.
We see that in this light the problem of understanding the critical properties realized in nature boils down to the classification of all the different fixed points. In two dimensions we know that every unitary scale invariant theory is also conformal invariant, so the problem further reduces to the classification of all possible conformal field theories (CFT). This can be done via algebraic methods, exploiting the properties of the associated Virasoro algebra [3,4].
A fixed point theory can then be deformed by adding weakly coupled operators that trigger a nontrivial flow out of the fixed point. By considering the linearization of this flow we can obtain all the remaining CFT data (like scaling dimensions and other critical exponents) that characterize the physical system and the way in which it responds to deformations. This is also the main idea of conformal perturbation theory.
So far our discussion has been limited to the neighborhood of a fixed point. The next natural step is to try to gain more information on the global properties of theory space. Such information is provided by Zamolodchikov's c-theorem [5], which states that in every unitary Poincaré invariant theory there exists a function of the coupling constants, the JHEP07(2014)040 c-function, that decreases from CFT U V to CFT IR , and that is stationary at the endpoints of the flow, where its value equals the central charge of the corresponding CFT. Note that the difference between the two central charges is an intrinsic quantity (intrinsic meaning independent of spurious contributions like scheme dependence of the renormalization procedure), so the content of the theorem is highly nontrivial.
In this case a complete RG analysis requires the ability to follow the flow arbitrarily far away from a fixed point. Unless the two fixed points are sufficiently close to each other we cannot rely on perturbative schemes. The non-perturbative framework we will use to address these issues is the functional renormalization group (fRG) based on the effective average action (EAA) [6]. The EAA is a functional whose scale dependence is given by an exact flow equation [7] which, being exact, allows to explore non-perturbative aspects. A first application of exact RG equations to the c-function has been explored in [8,9].
The main purpose of this paper is to move the first steps necessary in order to give a bridge between these two general results: the c-theorem and the computation of universal quantities related to the integrated flow between fixed points (that is, to global properties of theory space), and the fRG formalism based on the exact flow for the EAA. Our approach will be mainly a constructive one. We will give a general recipe to construct a c-function compatible with Zamolodchikov's theorem within the fRG framework. After identifying a natural candidate for a scale dependent c-function, c k , we will be able to write an exact non-perturbative flow equation for it. Of course, there are only few cases in which the exactness of the flow equation can be used and one usually needs to resort to approximations. However, we will see that already for a simple truncation as the local potential approximation the flow equation gives results compatible with the c-theorem.
Our viewpoint will be based upon a curved space construction. The reason for this is twofold. First, this avoids having to resort to algebraic techniques or OPE analysis: the central charge, for instance, becomes the coefficient of the conformal anomaly, which in curved space becomes an operator anomaly in the one point function. Second, this is more suitable for functional techniques, as the derivation of the trace anomaly matching condition will show, and more useful to write a general effective action. Indeed, this construction will require an investigation of what is the general form of the EAA away from fixed points, since the usual expansion in local operators is incapable of giving a nonzero running for c k . Working in curved space the natural candidate for c k is the coefficient of the Polyakov action. We will take this as our definition for the c-function and leave for further study the mapping between our approach and the one based on local RG with spacetime dependent couplings [10,11].
The paper will be organized as follows. In section 2 we will construct a Weyl-invariant functional measure and discuss the form of a CFT on curved background. This will lead us to a re-derivation of the trace anomaly matching condition, from which the "integrated" c-theorem follows from known results [12]. We will then move on to discuss the scale dependent c-function, and obtain our flow equation for it, in section 3. This construction uses the EAA as the main tool, so in section 4 we will investigate its general form. In section 5 we discuss various applications of our formalism while in section 6 we put forward a simple relation between the beta function of Newton's constant and the running c-function. Section 7 is devoted to the conclusions.

JHEP07(2014)040 2 The integrated c-theorem
We start by reviewing the integrated c-theorem expressing the change of the central charge ∆c = c U V − c IR along a RG trajectory connecting two fixed point theories, or equivalently two CFTs. We will work in curved space where the central charge, or equivalently the conformal anomaly, can be seen as the coefficient of the Polyakov term in the effective action. When we specify the background metric to be of the specific form g µν = e 2τ δ µν , with τ the "dilaton", ∆c becomes the coefficient of the operator τ ∆τ and can be easily extracted. But before we need to briefly discuss functional measures in curved space, Weyl-invariant quantization and the form of the effective action for a CFT on a curved background.

Weyl-invariant quantization and functional measures
The standard diffeomorphism invariant path integral measure in curved space [13], denoted here D I g , is Weyl-anomalous: under a Weyl transformation of the background metric g µν → e 2τ g µν and of the fields φ → e wτ φ, where w is the conformal weight of the field, 1 one encounters the conformal anomaly: where c is the central charge of the CFT, which we want to use as UV action in the path integral, and Γ W Z [τ, g] is the Wess-Zumino action: where ∆ ≡ −∇ µ ∇ µ is the Laplacian. The Wess-Zumino action can be integrated to give the related Polyakov action, which, upon Weyl variation, gives back (2.1): The Polyakov action generates to the following quantum energy-momentum tensor, which is anomalous: The conformal weight of the metric is wg = 2 in every dimension.

JHEP07(2014)040
This is the conformal anomaly, in the two dimensional case. In curved space, where it can be written in terms of curvature invariants, the conformal anomaly manifests itself already in the one-point function (2.5), while in flat space it is seen only starting from the two-point function. For example, in flat space the two point function of the energy-momentum tensor obtained from the Polyakov action, when written in complex coordinates, reproduces the standard CFT result [3]: This relation shows the equivalence between the central charge and anomaly coefficient. We can use the Polyakov action to define, formally, a new measure in the following way:  [14] with the exception that in that work a Stückelberg trick was used to maintain Weyl-invariance for any UV action.

CFT action on curved background
We have seen how to define, at least formally, a Weyl-invariant effective action starting from a Weyl invariant UV action via the functional measure (2.7), which is to be understood as the measure we will use from now on. Nevertheless on a curved background the effective action of a CFT is not Weyl-invariant since every CFT with c = 0 is anomalous, and thus its action must contain a Polyakov term. Still, in absence of relevant perturbations, quantization will just give the IR effective action equal to the UV action.
These considerations lead to the following "split" form for the effective action of a general CFT in presence of a background metric: (2.13) , defined by its Taylor series expansion in terms of correlation functions of φ, these being, in principle, exactly known. Very few CFT actions can be written in local form, these are the Gaussian, the Ising model (in the fermion representation) and the Wess-Zumino-Witten AKM actions [4]. S P [g] is the Polyakov action and c its central charge.
Other possible Weyl-invariant terms depending on the metric alone are not present in d = 2, but appear in higher dimensions. We now give an explicit example of this construction. The Gaussian theory has c = 1 and is the simplest example of a CFT: (2.14) Using the one-loop trace-log formula starting from the Gaussian UV action Γ U V we find: where the second term is due to the integration of the fluctuations, while the Polyakov term with the minus sign comes from the Weyl-invariant measure (2.7). The two cancel since the 1 2 Tr log ∆ = S P [g]. In order to have Γ U V = Γ IR one needs to add a relevant perturbation triggering the RG flow.

Anomaly matching from the path-integral
plus relevant operators, we can consider the IR effective action obtained by integrating out fluctuations: (2.16) Since the metric is non-dynamical we passed the Polyakov term through the path integral.
Here by relevant we mean, depending on the case, massive deformations or marginally JHEP07(2014)040 relevant ones. An example of the first are mass terms like m 2 2 φ 2 or mψψ, while Yang-Mills theory is an example of the second case.
If we now flow to an IR fixed point, by virtue of the splitting property (2.13), we . Choosing a dilaton background of the form g µν = e 2τ δ µν , we are left with: . In order to recover the flat space measure we first shift χ → e wτ χ and ϕ → e wτ ϕ and then use the invariance (2.7): Then we use the conformal invariance properties of the actions, i.e. we substitute Note that D δ χ ≡ Dχ is the flat space measure. The only remaining dependence on τ is due to the relevant terms, which make the path integral non-trivial. The last equality tells us that the dilaton effective action (generated by matter loops) compensates exactly the difference between the anomalies in the UV and IR. This is precisely the anomaly matching condition considered in [12,15].

Proof of the integrated c-theorem
We can now prove the integrated c-theorem following [12]. From equation (2.19), we can read off ∆c from the terms of the dilaton two-point function quadratic in momenta.
The relevant terms can be expanded in powers of τ : where Θ ≡ T µ µ and we omitted all terms of order τ 2 or greater since it is easy to see that they will not contribute to τ ∆τ . We are thus interested in the following expectation: We now only have to expand τ y around τ x : use translation invariance and compare with the coefficient of τ ∆τ to find: which is the integrated version of the c-theorem. From here one simply notices that the integral is positive due to reflection positivity and concludes that ∆c ≥ 0 [4,5].

JHEP07(2014)040 3 Flow equation for the c-function
The c-theorem states [5] that for a two-dimensional unitary quantum field theory, invariant under rotations and whose energy-momentum tensor is conserved, there exists a function c of the coupling constants which is monotonic along the RG flow and, at a fixed point, is stationary and equal to the central charge of the corresponding CFT. Therefore this function c is such that ∂ t c < 0 (where the "RG time" is given by the logarithm of the radius t = log r, so the flow is towards the infrared for r → ∞, hence the minus sign). The differential equation for c can be integrated from r = 0 to r = ∞ and gives back (2.24). A natural trial definition for an interpolating c-function is given by taking (2.24) with the integral which has been cut off at some scale µ (see for instance [16]): We will follow a different approach. Instead of cutting off directly in real space we will cutoff in momentum space. This will allow us to naturally connect with the framework of the functional Renormalization Group (fRG) and to derive an exact RG flow equation the c-function.

The fRG flow equation for the c-function
To construct the c-function is to consider a Wilsonian RG prescription. A clever way to do the momentum shell integration in a smooth way, is to introduce a suppressing factor in the path integral via D g χ → D g χ e −∆S k [ϕ,g] . The role of the cutoff action ∆S k [ϕ, g] is to restrict the integration to modes above the IR scale k. In this way we obtain a scale dependent effective action Γ k [ϕ, g], which, using (2.13), can be decomposed as: where S k [ϕ, g] is defined by S k [0, g] = 0 and c k is the scale dependent c-function. By "gravitational terms" we mean the purely geometrical terms depending on the metric alone, like √ g or √ gR, generated by fluctuations. The collection of the Γ k [ϕ, g] for all k constitute the RG trajectory connecting Γ U V [ϕ, g] to Γ IR [ϕ, g]; a cartoon of this shown in figure 1. If we now repeat the steps leading to equation (2.20), but with the cutoff term added, we arrive at: Now a derivative of (3.3) with respect to the "RG time" t = log k gives the RG flow of the central charge: in which the expectation value is calculated within the regularized path integral. We see that we obtain the flow of the c-function if we are able to evaluate the r.h.s. of (3.4), Figure 1. Cartoon depicting the flow in theory space: if the UV action Γ U V satisfies the Wess-Zumino relation (2.10) with c = 0, then the IR action Γ IR satisfies the Wess-Zumino relation (3.18) and the EAA, interpolating between the two, must satisfy (3.17 after specifying the form of the cutoff action. The running of c k is related to the coarsegrained dilaton two-point function. To understand how to handle this equation, we need to introduce the effective average action (EAA).
In the functional RG framework, one considers an IR regulator quadratic in the fields: chosen to suppress field modes in a covariant way: if φ n is an eigenfunction of the covariant Laplacian ∆φ n = λ n φ n , R k will act as a mass insertion for modes with λ n k 2 , while leaving unchanged the ones with λ n k 2 . In this way we obtain a scale-dependent partition function: The effective average action is then defined as the (shifted) Legendre transform: where ϕ = φ and J ϕ is obtained by inverting the solution of δW k [J,g] δJ = ϕ J . By using its definition in the path integral, one finds the integro-differential equation satisfied by the EAA (in which φ = ϕ + χ): (3.8) The main virtue of these definitions is that the EAA satisfies an exact RG flow equation [7]. A scale derivative of (3.8) gives: in which the expectation values are calculated with the fRG-regularized path integral. Using the fact that the Legendre transform of the generator of connected correlation functions is Γ k + ∆S k , we have: . (3.10) Substituting back in the previous expression we find the functional RG equation satisfied by the EAA: This equation is well defined, exact and offers a way to define QFTs non-perturbatively [6]. From the exact flow equation for the EAA we obtain a corresponding equation for the c-function. In particular, we can express the r.h.s. of (3.4) using (3.9): (3.12) Equation (3.12) is the exact flow equation for the c-function in the fRG framework. Using (3.11) in the r.h.s. leads to the following explicit form: perform two functional derivatives of (3.12) with respect to the dilaton, set τ = 0 and extract the term proportional to ∆: where all quantities are evaluated at ϕ = τ = 0. Note that in (3.14) we had to derive the cutoff kernel R k , since this depends explicitly on the dilaton. As shown in [18], these additional terms in the flow equation for the proper-vertices are crucial in maintaining background symmetry when employing the background field method. The flow equation in the form (3.14) is a bit cumbersome so we introduce a compact notation to rewrite it in a simpler way. If we introduce the formal operator∂ t = ∂ t R k ∂ ∂R k , we can rewrite the flow equation (3.12) for the c-function as: where we used the following simple relations: Now we can rewrite the flow equation (3.14) in the following compact form: where again all quantities are evaluated at ϕ = τ = 0. This is the form that we will use in applications in section 5. Finally, we can represent diagrammatically the two terms on the r.h.s. of (3.16) as in figure 2 and switch to momentum space to evaluate the diagrams by employing the techniques presented in [18]. In particular, continuous lines represent matter regularized propagators G k [0], while vertices with m-external wavy lines are the matter-dilaton vertices Γ

fRG derivation of the integrated c-theorem
We now rederive both the integrated c-theorem and the exact flow equation for the cfunction using a fRG theory space perspective. Away from a fixed point, apart for a Wess-Zumino term with running coefficient, that for the moment we call C k , there must be many additional terms spoiling the fixed point Wess-Zumino relation (2.10). Since these terms vanish at a fixed point they must be proportional to the (dimensionless) beta functions. We can thus make the following ansatz: This relation can be read as a generalized running Wess-Zumino action. The β-terms indicate terms proportional to (at least one) dimensionless beta function which vanish at the CFTs and are generated along the flow by the fact that we are moving away from criticality.
If we now use the Weyl-invariant measure to construct the EAA, then at the UV fixed point, that is for k → ∞, we must have C U V = 0. On the other hand, if we are not quantizing in a Weyl-invariant manner, we should reproduce the Wess-Zumino relation both at k = ∞ and k = 0. This tells us that in fact C k = c k − c U V if the UV theory is quantized in a Weyl-invariant manner and C k = c k if not. Weyl-invariant quantization corresponds, in the EAA formalism, to a constant shift of C k .

JHEP07(2014)040
Relation (3.17) can be used to give an equivalent RG derivation of the integrated ctheorem. When we flow for k → 0 to a fixed point theory CFT IR , relation (3.17) tells us that: (3.18) If we now set g µν = δ µν and expand Γ IR [ϕ, τ ] ≡ Γ IR [e wτ ϕ, e 2τ δ] in powers of the dilaton we find: The functional derivatives of the effective action are related to the traces of the energymomentum tensor: . (3.20) The first relation is identically zero at a CFT, i.e. Θ x IR = 0 . Inserting the second relation in (3.19) and expanding, as before, τ y around τ x using (2.23) gives immediately the integrated c-theorem (2.24). This derivation represents a consistency of the ansatz (3.17).
It is now clear that from the Wess-Zumino relation at finite k (3.17) we can easily read off the flow of the central charge. In this way, since ∂ t c k = ∂ t C k , from the coefficient of τ ∆τ in ∂ t Γ k [e wτ ϕ, e 2τ δ] we recover the exact RG flow equation for the c-function (3.12).
Another way to see that the flow of the c-function is given by (3.12) is to recognize that C k is nothing more than the coupling constant of the Polyakov action. As we said, when working on curved backgrounds one should always add the Polyakov term to a truncation. Thus the Wess-Zumino action on the r.h.s. of (3.17) derives from the presence of the Polyakov action, with coefficient C k , in the EAAs on the l.h.s. of the same equation. Then, as just seen in the previous paragraph, a t-derivative relates ∂ t c k to the two-point function of the dilaton. In principle one can obtain the flow of C k directly as the coefficient of √ gR 1 ∆ R but this is more laborious. Finally, note that the inclusion of the Polyakov action with running central charge makes the truncation consistent with the conformal anomaly both in the UV and in the IR. To understand the β-terms we will consider, in the next section, the scale anomaly.

General form of the effective average action
In this section we put forward some requirements which an ansatz for the EAA should satisfy. These requirements are motivated from the fact that the EAA should reproduce some generic features of QFTs, namely the scale and the conformal anomaly. In particular we will try to shed light on the nature of the β-terms introduced in equation (3.17).

The local ansatz and its limitations
When studying truncations of the EAA, one generally starts by expanding the functional in terms of local operators compatible with the symmetries of the system: This equation defines the running coupling constants g i,k , which become the coordinates that parametrize theory space in the given operator basis.
A class of operators, which is not complete, but allows many computations to be performed analytically, is the one composed of powers of the field, i.e. O i [ϕ, g] = ϕ 2i and g i,k = λ 2i,k (2i)! . In this approximation, one usually re-sums the field powers into a running effective potential V k (ϕ) and equivalently considers the following ansatz for the EAA: known as local potential approximation (LPA). Within this this truncation the exact flow equation (3.11) becomes a partial differential equation: Even such a simple truncation is able to manifest qualitatively all the critical information relative to the theory space of scalar theories and in particular the fixed point structure [19,20].
However, the effective action usually contains also nonlocal terms. Some of these nonlocal terms are directly related to the finite part of the effective action [21,22], which generally has a complicated form encoding all the information contained in the correlation functions or amplitudes. These terms are not present in the LPA which can be seen as the limit where we discard all the momentum structure of the vertices.
Nevertheless there are other nonlocal terms that are non-zero only away from a fixed point: these are the β-terms introduced in equation (3.17). As we will explain in this section these terms are needed to recover known results and will play a central role in our computations. If we limit ourselves to the local truncation ansatz (4.1), then one finds that the flow equation for the c-function is driven only by the classical non Weyl-invariant terms, which is not correct. This is not due to the fact that the flow equation (3.12) is wrong, rather, it is the truncation ansatz (4.1) that is insufficient. Fluctuations induce the β-terms of equation (3.17) and we will see that these are crucial in driving the flow of the c-function.
We will argue that these nonlocal terms have a precise form. We will do this requiring the EAA to reproduce the scale anomaly.

Nonlocal ansatz and the scale anomaly
It is easy to understand the origin of the terms on the r.h.s. of (3.17) which are linear in τ : they are related to the scale anomaly. To see this let us rescale the fields and expand JHEP07(2014)040 the EAA in powers of the dilaton: where: defines the scale dependent energy-momentum tensor trace. In the IR the EAA reduces to the standard effective action, which generally is scale anomalous. If we start with some UV action deformed by terms of the form j g j d 2 x √ g O i , the corresponding scale anomaly in flat space reads: where d i are the dimensions of the coupling constants. The expression in brackets is nothing but the beta function of the dimensionless coupling: This is a standard result known from both ordinary and conformal perturbation theories [4]. Now we consider again the β-terms on the r.h.s. of (3.17). They come from the conformal variation of the EAA which should include also the terms due to the scale anomaly. Therefore it is natural to generalize the above equation for a generic k: (4.8) If we insert this into (4.4) we find: This expression gives a non trivial flow of the c-function since we now have the vertex to insert in the r.h.s. of the exact flow equation (3.16). We now propose a covariant form for (4.9) using the following properties: With this and O i → e w i τ O i , it is easy to verify that the action

JHEP07(2014)040
reproduces (4.9) to linear order in τ . In order to get an ansatz consistent also with the conformal anomaly we need to add to (4.12) the Polyakov term with as coefficient the running central charge C k : The form (4.13) represents a parametrization of the EAA consistent with (3.17) to linear order in the beta functions and hints to what could be the general for of the EAA away from criticality. For the time being we will not improve further our ansatz, since we will see in the next section, that the understanding of the linear terms in the beta functions is already sufficient to build the c-function in some non-trivial cases. We hope to come back to the issue of higher order terms in τ , which may play a role in making a bridge between the fRG perspective adopted here and the ideas related to the local RG [10,11].

Checking exact results
Here we provide two examples where the c-function and the difference c U V − c IR are computed and can be compared to known exact results. We will consider a free scalar field and a free (Majorana) fermionic field whose fixed point actions are perturbed by a mass term, so they flow to c IR = 0.

Massive deformation of the Gaussian fixed point
We consider a scalar field with Gaussian action and c U V = 1 perturbed by a mass term. Since the beta function of the mass is zero (there are no interactions), our general ansatz (4.13) for the EAA reads: or when we rescale the fields: It's clear that the only interaction between φ and τ is the one induced by the dimension of the mass. In order to avoid possible vertices coming from the cutoff action we use the mass cutoff R k (z) = ak 2 which has the advantage of having no dependence with respect to the background metric. We have introduced the parameter a to check the cutoff independence of the result. After a short computation 3 we find the following flow: where m is the dimensionful mass. This RG flow occurs along trajectory-I of Integrating the above differential equation, with the initial condition c ∞ = 1 (the central charge of the Gaussian fixed point) we find: In the k → 0 limit this gives c 0 = 0 which implies ∆c = 1 independently of the cutoff parameter a. As expected a massive deformation of the Gaussian fixed point leads in the IR to a theory with zero central charge.

Massive deformation of the Ising fixed point
In this example we make a massive deformation of the Ising fixed point. The critical Ising model is described by a free Majorana fermion and a massive deformation of this correspond to consider T > T c [4]. According to our general ansatz (4.13) and considering that, as before, the mass beta function is zero, the EAA reads:

JHEP07(2014)040
or after the rescaling: The computation proceeds along the lines of the scalar case. Once again we use the mass cutoff R k = ak and we find: This RG flow occurs along trajectory-II of figure 3.
Integrating this equation with boundary condition c ∞ = 1 2 (the central charge of the Ising model) leads to which gives c 0 = 0 and ∆c = 1 2 as expected.

The c-function in the local potential approximation
The local potential approximation (LPA), introduced in section 4.1, is characterized by the action (4.2); in our case generalizes to: or after rescaling the fields: If we now pass to dimensionless variables, ϕ = k −wφ and V k (ϕ) = k 2Ṽ k (φ), then the second and third terms in the above equation, to linear order in τ , become V k (ϕ) − k 2 ∂ tṼk (φ)τ , so that the scalar-dilaton interaction is proportional to the dimensionless scale derivative of the potential.
To obtain the flow equation for the c-function we use (3.16) and the mass cutoff R k (z) = k 2 so that all cutoff vertices drop out. Only the first diagram of figure 2 contributes terms of order p 2 in the external momenta, more specifically we need to evaluate the integral: with the following regularized propagator: .
Here ϕ 0 is the minimum of the running effective potential, i.e. the solution of V k (ϕ) = 0.
With the mass cutoff one finds the following result:  provided thatṼ k (ϕ 0 ) > −1, since otherwise the momentum integral does not converge. Inserting this back in (5.10) finally gives: which is the flow equation for the c-function in the LPA with a mass cutoff. This the main result of this section. Note that since (5.13) is valid under the conditionṼ k (ϕ 0 ) > −1, the c-theorem ∂ t c k ≥ 0 is indeed satisfied within the LPA.

Flow between the Gaussian and Ising fixed points
We now consider the simple case where there are just two running couplings parametrizing theory space, i.e. we expand the running effective potential in a Taylor series: where m 2 k is the mass and λ k the quartic self-interaction. Inserting (5.14) in the flow equation for the effective potential (4.3) and projecting out the flow of the two couplings gives, after passing to dimensionless variables m 2 k = k 2m2 k and λ 2 k = k 2λ k , the following system of beta functions:

(5.15)
This system has two fixed points: the Gaussian (m 2 k ,λ k ) = (0, 0) and the Ising (m 2 k ,λ k ) = (− 1 4 , 3π 2 ). The Gaussian fixed point has two IR repulsive directions, while the Ising fixed point has one IR repulsive and one IR attractive direction. The trajectories starting along these directions are shown in figure 3, in particular trajectory-III connects the two fixed points. One can see that this last trajectory is the set of points where the dimensionless mass beta function vanishesλ k = −8πm 2 k (1 +m 2 k ).

JHEP07(2014)040
We can now use (5.13) to evaluate the c-function in this truncation. This turns out to be simply related to the square of the dimensionless mass beta function: As for (5.13), the result is only valid form 2 k > −1, so in this range we do have ∂ t c k ≥ 0, which is consistent with the c-theorem. The flow (5.16) is similar to the one given in [9], which was there found "by trial and error". Equation (5.16) is the first non-trivial example of explicit flow equation for the c-function obtained using the procedure presented in this work. In figure 4 we plot ∂ t c k in the plane (m 2 k ,λ k ): one can see that the magnitude of ∂ t c k is smaller along a "valley" containing the two fixed points. Along this valley lies the trajectory connecting them, trajectory-III of figure 3.
We would like to compute ∆c by integrating the flow of the central charge along the path connecting the Gaussian and Ising fixed points, but in this simple truncation trajectory-III is defined by the vanishing of the dimensionless mass beta function, and thus ∂ t c k is zero along it. To find a non-trivial result, we need to consider a more refined truncation ansatz for the running effective potential. We leave these studies to future work.

Sine-Gordon model
We now consider the Sine-Gordon model which, in the continuum limit, is described by the following action [4]: where m is the mass and β is a coupling constant. This theory can be seen as a massive deformation of the Gaussian fixed point action (with c U V = 1) and indeed we will find c IR = 0.
The Sine-Gordon model can be described by an LPA with effective potential: We find the following form for the beta functions of m k and β k : wherem 2 k = m 2 k /k 2 is the dimensionless mass. Inserting the Sine-Gordon running potential (5.18) into the flow equation (5.13) now gives: We solved the system of equations numerically imposing c U V = 1 finding ∆c 0.998, in satisfactory agreement with the exact result ∆c = 1. In figure 5 we plot the running of c k as well as its beta function.

The c-function in the loop expansion
The last approximation we will consider is the loop expansion. The exact flow equation for the EAA (3.11) can be solved perturbatively [23,24] loop by loop. We review this in the appendix. In the first part of this section we will look at the various contributions diagrammatically, while in the second part we will explicitly evaluate one subclass of these.

Zamolodchikov's metric: diagrammatics
Using relation (A.4) we can compute the running of the EAA at each order in the loop expansion. The running of the L-th term ∂ t Γ L,k , say, will contain a contribution to the running of c k that we will call ∂ t c L,k . The term c L,k arises only from diagrams with L matter loops and two dilaton external lines. In this way we can build a loop expansion for the c-function. We can start by applying this construction step by step so to make clear how everything works. We will work with a Z 2 -symmetric scalar theory, so that the part linear in the dilaton of our general ansatz (4.13) takes the form: whereβ 2 is the mass beta function,β 4 is the ϕ 4 coupling beta function, and so on. At one loop, we have only the following diagram, obtained from (A.5) of the appendix by functional derivation with respect to the dilaton,

JHEP07(2014)040
Here we adopt the same diagrammatic rules of section 3.2 where the continuous line represents the regularized propagator (in this case given in equation (A.6) of the appendix), while the wavy line represents the dilaton. On every diagram the operator∂ t acts, but in this case it is just ∂ t . In this diagram the vertices, are derived from (5.20), is the mass beta function, so this contribution goes likeβ 2 2 and we recover the LPA result (5.13) as one would expect.
From the flow of the two-loop contribution, (A.11) of the appendix, we obtain different terms. We get the "non-diagonal" contribution (we will make this jargon clear in a second): proportional toβ 2β4 . Together with this, we also have the following 2-loop diagonal contributions: which are proportional to λ 2β 2 2 . These represent a diagonal but coupling-dependent contribution, in the sense that couplings do not only appear through the beta functions. When going to 3-loops, 4-loops and so on, corresponding diagrams must be considered for all the diagonal contributions.
At three loops (remember we are considering a Z 2 -symmetric theory, so there are no scalar odd power interactions) we get again the "diagonal" contributions: both proportional toβ 2 4 , as well as a nondiagonal one: proportional toβ 2β6 . From these first diagrams we clearly see that from the structure of the loop expansion we only get terms quadratic in the beta functions. We can indeed follow Zamolodchikov and define the "metric" g ij through: ∂ t c k = g ijβ iβj . the entries of this metric. In principle one can evaluate all these diagrams for a generic cutoff R k (z) but this turns out to be a difficult analytical task. In the next section we will be able to evaluate analytically one particular class of diagonal entries. 4

Diagonal contributions
At L-loop order, the simplest coupling-independent diagonal contribution comes from the following diagram: ...
corresponding to the expression: (which generalizes equation (A.15) of the appendix). In the above equation the 2(L + 1)! comes from the symmetry factor of the diagram, and the minus sign from the fact that we are acting with an overall∂ t . To recover the contribution to ∂ t c k is simple: expand τ y around x as in equation (2.23), and isolate the proper term according to equation (3.12).
To see more explicitly the form that the metric of Zamolodchikov takes, we need some preliminary results. Using a mass cutoff R k = k 2 , the zero mass running renormalized propagator (A.6) will be the same as the standard massive one, only with k 2 in place of the mass m 2 , and the cutoff vertices play no role. In real space the propagator reads: where K 0 is the Bessel K-function of order zero. We introduced the parameter a, eventually to be sent to 1, since in this way we have the simpler formulã The different contributions are then calculated after expanding τ y around x using (2.23). We find: These diagonal terms can be written to all orders, they give a contribution to the flow equation for c k of the form: in which we defined the quantity

JHEP07(2014)040
0.0182 4.778 × 10 −4 1.485 × 10 −5 5.066 × 10 −7 1.825 × 10 −8 6.8 × 10 −10 Note the interesting thing that contributions at loop order L are proportional to the square of the beta function of the couplingλ L+1,k . Thus the flow of c k receives contributions from all loops (as it is inherently non-perturbative) but a given interaction starts to contribute only at a given loop order. All the A L can be evaluated numerically and they turn out to be positive. The numerical values of the first A L are shown in table 1. Note the fast decrease relative to the one-loop value.
We can now write down the contribution of this class of diagrams to the running of the c-function at all loops in the Z 2 -symmetric case: which also gives the explicit form for the diagonal entries of the Zamolodchikov metric. Since this sum is manifestly positive, we can say that the c-theorem is satisfied to all loops by the diagonal terms considered.
As we have seen previously, the entries of Zamolodchikov's metric contain a couplingindependent piece, plus further pieces proportional to increasing powers of the coupling constants, as we increase the loop order. The positivity properties of the metric are far from trivial when all these terms are involved. However, when the couplings are sufficiently small, the positivity will be determined solely by the coupling independent terms.

Non-unitary theories
Finally we make a comment on when the c-theorem is not satisfied, i.e. the case when ∂ t c k < 0. We know that the c-theorem does not hold without the unitarity assumption [5]. This can indeed be checked explicitly. It's easy to see that when one considers interactions with complex couplings then the coefficients in the loop expansion turn negative. For instance, one notable example is the Lee-Yang model [4], in which one introduces the non-unitary complex interaction: A simple analysis reveals that this interaction contributes to the running of c k through the following diagram: which turns out to have the wrong sign to be consistent with the c-theorem:

JHEP07(2014)040
6 The c-function and Newton's constant In this section we derive an interesting relation between the c-function and the matter induced beta function of Newton's constant. 5 This can then be used to obtain another form of the flow of the central charge ∂ t c k .

Relation between c k and β G k
To obtain this relation we need to consider what happens when in equaiton (4.12) we set O = R . Since the coupling constant of the invariant √ gR is − 1 16πG k , where G k is the running Newton's constant, one finds, for the gravitational part of the EAA, the following form: 6 We recognize that the Polyakov term above is the same that we included in our general anstaz for the EAA (4.13). Thus we infer that there is a relation between the beta function of Newton's constant and the running c-function: This is a nontrivial statement by itself. It tells us that the running c-function for a certain matter field type can also be computed from the contributions of that kind of matter to the beta function of Newton's constant. In fact a derivative of (6.2) with respect to the RG scale gives (remember that ∂ t C k = ∂ t c k ): where β G k ≡ ∂ t G k is the Newton's constant beta function. We will check the consistency of relation (6.2) in the case of a minimally coupled and a self-interacting scalar.

Minimally coupled scalar
Consider a minimally coupled scalar describing a massive deformation of the Gaussian fixed point as discussed in section 5.1.1. The action is given in (5.1) and the exact flow equation (3.11) for this case reads: .
Note that the dilaton plays no role now, since we are free to set τ = 0. Instead, to find c k using (6.2) we need to extract the terms in the trace on the r.h.s. of (6.4) that are 5 In what follows we identify the Newton's constant as the coupling in front of the Ricci scalar. In a non linear sigma model on curved target space this coupling is equivalently identified as the dilaton constant mode. 6 We need here 1/4 instead of 1/2 because of the the further symmetry we have in exchanging the two Rs. . c k and ∂ t c k as a function of k for a massive deformation of a minimally coupled scalar. Mass (a = 1), optimized (a = 1) and exponential (a = 1, b = 1) cutoffs (upper curves), exponential (a = 1, b = 1 2 ) cutoff (middle curves), exponential (a = 1, b = 3 2 ) cutoff (lower curves). In all cases we set m 2 = 1.
proportional to the invariant √ gR. As usual, this can be done using the heat kernel expansion [25]. Defining h k (z) = ∂tR k (z) z+m 2 +R k (z) , one finds: 1 2 Tr h k (∆) √ gR = 1 8π which, when compared with the scale derivative of − 1 16πG k √ gR on the l.h.s. of (6.4), gives: Thus our formula (6.2) leads to: 7 Note that this relation is valid for arbitrary cutoff function R k (z), as opposed to the result of section 5.1.1 valid only for the mass cutoff. For both the mass cutoff R k (z) = ak 2 and the optimized cutoff R k (z) = a(k 2 − z)θ(k 2 − z) we find the following form: For the exponential cutoff R k (z) = az e bz/k 2 −1 , with parameters a and b, we find: (6.9) 7 We are not using Weyl quantization so C k = c k JHEP07(2014)040 Figure 7. ∂ t c k in the (m 2 k , λ k ) plane according to (6.13) for a = 1 and b = 1 2 . We marked with a red dot the position of the Gaussian and Ising fixed points and the trajectory-III connecting them. One can note that the trajectory connecting the two fixed points lies along a "valley" where the variation in the flow of the central charge is zero.
In all cases and for all values of the parameters a and b we find that c U V = 1 and c IR = 0 as expected. A derivative of (6.9) gives the flow of the c-function: The interpolating c k of equation (6.9) and the flow of the last equation are shown in figure 6.
We clearly see that the flow is scheme dependent, but the integral of it along a trajectory, giving ∆c, is universal.

Self-interacting scalar
We consider now an interacting scalar, i.e. the LPA action (5.9) of section 5.2. We can obtain c k directly from equation (6.9) by just making the replacement m 2 → V k (ϕ 0 ): A scale derivative now gives: We need to decide the value of ϕ 0 where to evaluate this expression. In this case it is important to distinguish the ordered from the broken phase. If the running effective potential has the polynomial form (5.14), then we have ϕ 0 = 0 in the ordered phase and ϕ 0 = ± 6m 2 k /λ k in the broken phase, the two phases being separated by trajectory-III and its continuation. Inserting these expressions in (6.12) gives the following form for the  As shown in figure 7, the flow (6.13), even if not proportional to the square of the dimensionless beta function, is positive ∂ t c k ≥ 0 in the (m 2 k , λ k ) plane, it is zero at the fixed points and on trajectory-III connecting them.
This calculation represents a non-trivial check of relation (6.2) and shows how this relation can be used explicitly to compute c k in a given truncation by means of heat kernel techniques.

Conclusions
In this work we have explored a new way to study the flow of the c-function within the framework of the functional RG based on the effective average action (EAA). This function interpolates between the UV and IR central charges of the corresponding CFTs and is thus a global feature of the flow, related to the integration of it along a trajectory connecting two fixed points, independent of scheme ambiguities.
Our main result is an RG exact equation for the running c-function based on the identification of it with the coefficient of the running Polyakov action. This equation relates the flow of the central charge to the exact flow of the EAA. To solve the equation for non-trivial cases we built a suitable ansatz requiring the EAA to reproduce generic features of QFTs, namely the scale and the conformal anomalies. In its own right this is an interesting result since it teaches us that a consistent ansatz for the EAA off criticality should include some nonlocal terms proportional to beta functions. Of course we do not claim full generality for this ansatz, but we found that it is sufficiently accurate to trigger the flow of the c-function in non-trivial cases. Explicit computations, within the local potential approximation and the loop expansion, have been presented in section 5 showing the compatibility of our framework with the c-theorem.
Moreover we have put forward a relation between the beta function of Newton's constant and the running conformal anomaly. This relation comes from internal consistency of the generic ansatz for the EAA we proposed and allows us to use heat kernel techniques to compute the RG running of the c-function. We also checked this other relation in explicit cases, finding it consistent. Nevertheless we point out that our analysis is not complete. The works [10,11,16] highlight that there are some subtleties related to the definition of the c-function. A complete mapping between the local RG approach and the fRG is still lacking and further study is needed in this direction. Another issue, which has not been touched at all, is the generalization of these ideas to the higher dimensional case, in particular d = 4 where one can consider similar constructions for the a-function [10,11,26,27], which we leave to future work.

JHEP07(2014)040
A Loop expansion from the fRG The exact flow equation (3.11) satisfied by the EAA can be solved iteratively. If we choose as seed for the iteration the bare action, then the iteration procedure reproduces the renormalized loop expansion [23,24].
One starts with Γ 0,k ≡ S Λ , where S Λ is the UV or bare action, and sets up an iterative solution (the subscript 0 indicates the order of the iteration, Λ is the UV cutoff and k is the RG scale) by plugging Γ 0,k into the r.h.s. of the flow equation and integrates the resulting differential equation with the boundary condition Γ 1,Λ = S Λ . The solution Γ k,1 is then plugged back into the r.h.s. of the flow equation and the procedure is be repeated.
To see this let us introduce as a loop counting parameter and expand the EAA: The bare action is k-independent ∂ t S Λ = 0. The exact flow equation (3.11) now takes the form: 2,k [ϕ] + . . .

. (A.2)
The original flow equation (3.11) is finite both in the UV and IR: to maintain these properties the bare action S Λ has to contain counterterms to cancel the divergencies that may appear in the Γ L,k . Thus we define: where each counterterm ∆S L,Λ is chosen to cancel the divergent part of Γ L,0 . Since this divergent part is the same as the divergent part of Γ L,k (we refer to [24] for more details on this point), this choice renders the denominator of (A.2) finite. Here S 0 is the renormalized action, i.e. the bare action with renormalized fields, masses and couplings. From (A.2) we can read off the flow of the L-th loop contribution: (A.4) The one-loop equation is straightforward: where the k-dependent renormalized propagator, depends on k only trough the cutoff R k . Thus, within the loop expansion, the operator∂ t , introduced in section 3.2, is equivalent to ∂ t .