Functional perturbative RG and CFT data in the $$\epsilon $$ϵ-expansion

We show how the use of standard perturbative RG in dimensional regularization allows for a renormalization group based computation of both the spectrum and a family of coefficients of the operator product expansion (OPE) for a given universality class. The task is greatly simplified by a straightforward generalization of perturbation theory to a functional perturbative RG approach. We illustrate our procedure in the $\epsilon$-expansion by obtaining the next-to-leading corrections for the spectrum and the leading corrections for the OPE coefficients of Ising and Lee-Yang universality classes and then give several results for the whole family of renormalizable multicritical models $\phi^{2n}$. Whenever comparison is possible our RG results explicitly match the ones recently derived in CFT frameworks.


INTRODUCTION
The standard perturbative renormalization group (RG) and the -expansion have been, since the pioneering work of Wilson and Kogut [1], the main analytical tools for the analysis of critical phenomena and, more generally, for the study of universality classes with methods of quantum field theory (QFT). Under the pragmatic assumption that scale invariance implies conformal invariance at criticality, which is confirmed by almost all interesting examples, one could also argue that conformal field theory (CFT) methods serve as an additional theoretical tool to describe critical models.
The RG flow of deformations of a scale invariant critical theory in a given operator basis is generally encoded in a set of beta functions of the corresponding couplings. As demonstrated by Cardy [2], the beta functions can be extracted adopting a microscopic short distance cutoff as a regulator and in particular, expanding in the scaling operators, the linear part of the beta functions is controlled by the scaling dimensions of the associated operators, while the quadratic part is fixed by the OPE coefficients of the operators involved in the expansion. Whenever the underlying critical model is a CFT, this approach is the foundation of conformal perturbation theory and its development strengthens further the conceptual relation between RG and conformal methods.
A CFT can be fully characterized by providing the so-called CFT data, which includes the scaling dimensions ∆ i of a set of operators known as primaries, and the structure constants C ijk of their three point functions [3][4][5]. From the point of view of CFT, the scaling dimensions determine some of the most important properties of the scaling operators at criticality, and in fact can be related to the critical exponents θ i of an underlying second order phase transition, while the structure constants provide further non-trivial information on the form of the correlators of the theory. The CFT data can be used, in principle, to reconstruct the full model at, or close to, criticality. In dimension greater than two (see, for example, [3][4][5]), however, since the symmetry group is finite dimensional, the use of analytical CFT methods is often not simple and in fact most of the recent success of CFT applications comes from the numerical approach known as Conformal Bootstrap [6].
Up to now, RG methods have been almost always devoted to the computation of the RG spectrum within a perturbative analysis in the -expansion below the upper critical dimension of a given universality class. The determination of the RG spectrum practically overlaps with the computation of the critical exponents and thus of the scaling dimensions ∆ i of the underlying CFT. It is thus natural to wonder to which extent the RG can help the determination of the remaining CFT data: the structure constants C ijk , which have received far less attention in the RG literature.
The question which arises spontaneously is whether the approach of Cardy [2] can be reversed and used to derive some of the OPE coefficients once the RG flow of a model is known. In such a framework, since the beta functions are generally computed in a specific RG scheme, one could expect that these RG based OPE coefficients might show some degree of scheme dependence.
The main purpose of this paper is to present an RG based approach, in a dimensionally regularized MS scheme, to the computation of the OPE coefficientsC k ij related to the quadratic part of the Taylor expansion of the RG flow around a critical point. We shall also pay attention to the transformation induced by general scheme changes among mass independent schemes, and infer some structure constants C ijk when scale invariance implies conformal invariance, strengthening in this way the link between RG and CFT.
We also show how the upgrade from standard perturbative RG to functional perturbative RG allows for a more straightforward access to these quantities. After illustrating how to do this for two representative cases, the Ising and (for the first time) the Lee-Yang universality classes, respectively realized as unitary and non unitary theories, we also proceed to the construction of the beta functions for all the unitary (even) models relying heavily on the approach developed by O'Dwyer and Osborn [7]. In general the use of the functional approach simplifies the computation of beta functions, from which, in the vicinity of a fixed point, one can try to extract some of the (universal) CFT data ∆ i and C k ij , from linear and quadratic perturbations around the critical point, respectively. This paper is concerned with fleshing out the main features of the func-tional approach and applying them to the rich variety of critical theories which can be described with a single scalar field φ. The functional approach appears to be very powerful because the beta functionals are of a strikingly simple form, and yet at the same time they describe the scaling behavior of classes of infinitely many composite operators.
In an effort to better understand the possible RG scheme dependence of the OPE coefficientsC k ij we carefully review their transformation properties. 1 In our approach we can compute only the subset of OPE coefficients which are "massless" at the upper critical dimension (the others being projected away by the dimensionally regularized scheme) and therefore less sensitive to -corrections induced by a change in the RG scheme. We show that our next-to-leading-order (NLO) computation gives these OPE coefficients at order O( ) and reproduces the structure constants previously obtained in a CFT framework [9,10]. This fact, even if plausible, is in general not obvious because of the possible scheme dependence, and we find it to be supported by the functional approach, which indeed constrains to some extent the possible choices of coupling redefinitions that otherwise would be completely arbitrary. All other "massive"C k ij strongly depend on the computational scheme and vanish in dimensional regularization. We observe that some OPE coefficients, including "massive" ones which would thus require a separate investigation, can be or have already been obtained for several universality classes in the -expansion in a CFT framework, with either the CFT/Schwinger-Dyson Bootstrap [9,[11][12][13][14][15], the perturbative conformal block techniques [10] or Mellin space methods [16]. Also large spin expansion techniques could be useful [17].
The first step in the functional perturbative RG approach is the computation of the beta functional β V of the effective potential V(φ), which generates the beta functions of the couplings of all the local operators φ k . This can often be used to verify our results by checking them against the renormalization of the relevant operators. Next comes the inclusion of the beta functional β Z of a field dependent wavefunction Z(φ), which generates the flow of the couplings corresponding to operators of the form φ k (∂φ) 2 and, through its boundary conditions, allows also for the determination of the anomalous dimension η. Higher-derivative operators can be added on top of the aforementioned ones following a construction based on the derivative expansion, which treats operator mixing in a systematic way, a topic that will be discussed here later on.
The content of the paper is as follows: In Sect. 2 we show in general how to use the RG to compute both the spectrum of scaling dimensions and the coefficients of the OPE. We discuss in general the possible scheme dependence by studying their behavior under arbitrary changes of parametrization of the space of all couplings, and we illustrate our methods by considering the RG flow of the Ising [18,19] universality class as an example. 2 In Sect. 3 we promote the standard perturbative RG to functional perturbative RG and illustrate the procedure by applying it to the Ising and Lee-Yang [20][21][22][23][24][25][26] universality classes. Using the beta functionals for the effective potentials in these two examples, we give general formulas for both the spectrum and the structure constants of the underlying CFTs, and use them to highlight the main novelties of the approach. In Sect. 4 we describe how to systematically improve the functional approach, to include arbitrary higher derivative operators, and how to generally deal with operator mixing. In Sect. 5 we present an application of this framework to the study of the general multicritical φ 2n universality class [7,27]. Finally, in Sect. 6 we draw some conclusions and discuss the most important prospects of this approach.
In appendix A we review the perturbative computations which are necessary to obtain the beta functionals used in Sect. 5. In appendix B we show how to use the functional approach to prove some simple scaling relations between critical exponents generally known as shadow relations. In appendix C we spell out some intriguing relation between the perturbative and non-perturbative functional RG approaches [28][29][30][31].

SPECTRUM AND OPE COEFFICIENTS FROM RG ANALYSIS
The primary goal of the RG analysis is the study of universality classes and the determination of their quantitative properties, i.e. the CFT data (when the two are related). This data is the union of the spectrum (the set of scaling dimensions ∆ i of composite operators) and the set of structure constants (in a CFT these are in one-to-one correspondence with the OPE coefficients C k ij of primary operators). The aim of this section is to introduce a computational scheme which shows how CFT data is (partially) encoded in the beta functions describing the RG flow in proximity of a fixed point. We start by describing the picture recalling a picture inspired by Cardy [2], which was originally defined in a short distance regularized scheme, and considering a generic basis of operators in which possible mixing effects are present. Then we present, in a dimensionally regularized scheme, a simple discussion of the Ising universality class to provide an example of an RG determination of OPE coefficients in the -expansion, which will also motivate the subsequent discussion of the transformation properties of theC k ij . This discussion will make clear which subset of OPE coefficients can actually appear in the beta functions, finally explaining which part of the CFT data is directly accessible by our RG methods.

General analysis
We begin our analysis by considering a general (renormalized) action in d dimensions, describing an arbitrary point in theory space. The choice of a basis set of operators Φ i allows the introduction of coordinates, i.e. the corresponding (dimensionless) couplings g i . The scaling dimensions ∆ i of the (composite) operators, as we will see in a moment, take precise values only in the vicinity of a fixed point of the RG flow. All information regarding the flow can be extracted from the set of (dimensionless) beta functions which are in principle fully computable once a given scheme is precisely defined. A fixed point of the RG flow is the point g i * in the space of couplings for which the theory is scale invariant In the neighborhood of a fixed point it is convenient to characterize the flow by Taylor expanding the beta functions. If δg i parametrizes the deviation from the fixed point where at the linear level we defined the stability matrix and at the quadratic level we defined the tensor which is symmetric in the last two (lower) indices. Each scale invariant point of the RG flow is in one to one correspondence with a universality class and, under mild conditions that we assume, a related CFT. The spectrum of the theory at criticality is given by the eigendeformations of M i j with the corresponding eigenvalues being (the negative of) the critical exponents θ a . We will only be concerned with cases in which either the matrix M i j is already diagonal, or its left and right spectra coincide (meaning that the spectrum is unique and unambiguous). It is convenient to introduce the rotated basis Critical exponents allow for a precise definition of the scaling dimensions of the operators through the relation θ i = d − ∆ i . Let us introduce the "canonical" dimensions D i of the couplings, and parametrize the deviations of the critical exponents from the canonical scaling through the anomalous dimensionsγ i as Here and in the following we adopt a tilde to distinguish RG quantities from CFT ones. The notion of canonical dimension is in principle arbitrary, but in real-world applications it is generally borrowed from the scaling of the Gaussian critical theory. This expression is, strictly speaking, valid only for primary operators; for descendants there is a subtlety that we will discuss later. The matrix (S −1 ) i a also returns the basis of scaling operators of the theory at criticality, O a = ∑ i (S −1 ) i a Φ i , so that we can rewrite the action as a fixed point action (i.e. CFT action) plus deformations Deformations are relevant, marginal or irrelevant depending on the value of the related critical exponent (respectively positive, zero or negative). In the diagonal basis also the tensor N i jk have a direct physical meaning, since after the diagonalizing transformation it becomes a quantity related to the (symmetrized) OPE coefficients 3 It will become clear in the practical examples that will follow this subsection that at d = d c theC a bc are the OPE coefficients of the underlying Gaussian CFT and that all O( ) corrections agree with CFT results for all available comparisons, despite the general inhomogeneous transformations of these coefficients under general scheme changes as discussed in subsection 2.3. For these reasons we make the educated guess that the quantities in (2.9) are the MS OPE coefficients since they have been computed using MS methods. The relation among the standard perturbative MS OPE coefficients and quadratic coefficients in the beta functions is an interesting subject, which is nevertheless beyond the scope of this work and is left for future investigations.
The beta functions can now be written as This formula is the familiar expression for beta functions in CFT perturbation theory (see, for example, [2]) and provides a link between RG and CFT. Generalizations of this result beyond the leading order are considerably less simple than what we presented here [8].
In CFT one uses the OPE 4 to renormalize a perturbative expansion of the form (2.8) in which the CFT is described by the action S * and deformations are parametrized by the couplings λ a . 5 In the RG 3 Note that the overall normalization of the OPE coefficients is not fixed: a rescaling of the couplings λ a → α a λ a impliesC a bc → α b α c α aC a bc . 4 These OPE coefficients are related to those entering the beta functions by a factor S d /2 (see [2]). 5 The careful reader must have noticed that our determination of the C a bc is symmetrized in the lower two indices, but it is more than enough to reconstruct the fully symmetric structure constants C abc . framework, conversely, the knowledge of the beta functions could permit (in principle) the extraction of the conformal data directly from (2.10). The rest of this paper is essentially devoted to a detailed exploration of this link, first within a simple example in the next subsection and then, after a short discussion of the scheme dependences of the OPE coefficients, within a functional generalization of standard perturbation theory -expansion.

Example: Ising universality class
It is useful at this point to consider an explicit example to introduce our approach, the Ising universality class in d = 4 − [18,19]. Perturbation theory forces us to restrict to deformations around the Gaussian fixed point, the simplest of which are power like non derivative operators Φ i = φ i parametrized by the dimensionless couplings g i , as will be shown in the next section.
One can obtain the (two loop) NLO beta functions for relevant and marginal deformations 6 and the anomalous dimension η = 96g 2 4 . Note that the coefficients of the one loop leading-order (LO) quadratic terms in the couplings are directly related to the Gaussian OPE coefficients, which by construction coincide with the mean field OPE coefficients of the Ising universality class (see also [32]).
The fixed point is characterized by g * 4 = 72 + 17 2 1944 + O( 3 ) and g * 1 = g * 2 = g * 3 = 0. Around this fixed point one therefore expands in powers of deformations (with λ i = g i for i = 1, 2, 3 and λ 4 = g 4 − g * 4 ), and the beta functions become One should keep in mind that the NLO coefficients of quadratic terms involving λ 4 can be affected by diagonalization. From the above relations one can immediately read off the critical exponents θ 1 , θ 2 , θ 3 , θ 4 as minus the coefficients of the linear terms. Note that the scaling relation θ 1 + θ 3 = d discussed in appendix B is indeed satisfied. We will see that these couplings will not be subject to any further mixing and thus these are the complete -series to the exhibited order for the critical exponents and the OPE coefficients. We have limited the -series for the OPE coefficients to linear order since the O( 2 ) terms are incomplete, receiving contributions from next-to-next-to-leading-order (NNLO) beta functions. From the eigenvalues we can extract the coupling (RG) anomalous dimensionsγ i through the relations and η = 2γ 1 . The scaling dimensions of the composite operators are instead and define the (CFT) anomalous dimensions γ i . The difference between theγ i and γ i appears only when the related operators are descendant, in this case when i = 3 for which γ 3 =γ 3 + η. We will postpone the discussion of this fact to the appendix B. The explicit expressions for the first anomalous dimensions are well knowñ From (2.13) it is equivalently easy to read off the OPE coefficients (which on the nondiagonal entries are half the value of the coefficients in the beta functions) We note that the OPE coefficientC 1 14 is in perfect agreement with that found in [9] using CFT methods, while we did not find any result in the literature for the other two coefficients to compare to. An explanation of why this agreement is expected will be given in Sect. 3.1. It is also important to stress that we ensured the agreement by choosing the same normalization of [9], that is by fixing the coefficients of the two point functions.

Transformation properties
In general different regularization and renormalization procedures may result into non trivial relations among the renormalized couplings. These relations go under the name of scheme transformations, and they are exemplified through maps among the couplings of the two schemes that can be highly non-linear [33]. Whenever the scheme transformations are computed between two mass independent schemes (such as, for example, MS and lattice's 7 ) these relations might have a simpler form, but we will find that it is very useful to consider them in their most general form. Let g i =ḡ i (g) (2.17) be the general invertible, possibly non-linear, transformation between the set of couplings g i andḡ i . Under such a change of "coordinates" the beta functions transform as vectors 8β (2.18) Now we turn our attention to the fixed point quantities, therefore in the following it is understood that the transformations will be evaluated at a fixed point. The stability matrix transforms asM Since the derivatives are evaluated at the fixed point, the stability matrices of the two set of couplings are related by a similarity transformation. Therefore it is trivial to prove that the spectrum is invariant, meaning that it does not depend on the parametrization as one would naively expect for a physical quantity. Things become less trivial when considering the matrix encoding the second order of the Taylor expansion at the fixed point. A direct computation shows To simplify this expression it is convenient to assume that the couplings g k have already been chosen to diagonalize the stability matrix with a linear transformation, so that on the right hand side there will be the structure constants Now it is necessary to move to the basis of couplingsḡ i in whichM i j is diagonal, so that the structure constants appear on both sides. We finally find (2.23) 7 But in practice all lattice implementations can be considered massive schemes. 8 The suummation convention is understood in this subsection. which implies that the set of matricesC c ab has a transformation law that is not homogeneous and therefore is reminiscent of the one of a connection in the space of couplings [34,35].
In the context of conformal perturbation theory one can find a similar result in [8], in which the analysis includes cubic terms but is limited to a diagonal stability matrix because conformal perturbation theory adopts by construction the basis of scaling operators.
At this point few comments on the transformations ofC c ab are in order: • It is evident from (2.23) thatC c ab can be independent of parametrization for a very special sum condition among the scaling dimensions of the couplings or when the Hessian at the fixed point is zero. The latter case could be realized for a specific family of scheme transformations, while the former condition can be realized exactly only in d = d c , which corresponds to = 0, that is, (2.24) • We will observe in all practical examples that the coefficientsC c ab for which the condition in Eq. (2.24) holds are the ones that are accessible via dimensional regularization. We dub them "massless", as opposed to the "massive" ones that do not satisfy the above condition and are zero in dimensional regularization. Moreover, these "massless" OPE coefficients at the critical dimension are insensitive to changes of RG scheme and can be computed unambiguously with RG methods.
• In perturbation theory -expansion one can obtain -series only for the "massless" OPE coefficients. In particular in d = d c − one generally has θ c − θ a − θ b = O( ) and thus only the O( ) terms can be scheme independent if the Hessian is at least O( ). This allows for crucial comparisons and cross-checks with other theoretical approaches like CFT (as for theC 1 14 of our previous example). In fact we will see that all the MS leading corrections for the multi-critical models we can compare with CFT are in perfect agreement. While this agreement can be explained at the level of beta functions, by explicitly constructing the most general map between to orthogonal massless scheme that also preserves the -expansion, the explanation is more transparent when discussed in functional terms in the next section.
Finally one should remark that once the beta functions are extracted in some scheme, one might also envisage geometrical methods to extract quantities which depend on the universal scheme independent OPE coefficients to overcome the above limitations. A step in this direction has been made for functional-type flows in Ref. [36] in the context of the Polchinski RG equation. In this work, the authors define normal coordinates in the space of couplings which have both geometrical meaning and definite scaling transformations. In relation with the transformation (2.23), one can follow [36] and argue that all "massive" coefficients can be eliminated by an opportune transformation of the couplings and hence there exists a scheme, or rather a family of schemes, whose only coefficients are the scheme independent ones. This family can be appropriately named "family of minimal subtraction schemes" having the MS scheme as its most famous representative. We hope to address further these topics in future investigations.

FUNCTIONAL PERTURBATIVE RG: A FIRST LOOK
The previous example on the Ising universality class, which was dealing with the study of the O k = φ k deformations of the Gaussian fixed point, can be analyzed more conveniently if we work directly with the generating function of such operators which is the local potential V(φ). We thus consider the action and study the perturbative RG flow it generates. In particular, it turns out to be a smart move to perform background field computations of loop diagrams in which the field φ is set to a constant and leave the form of the potential completely general so that we can extract the beta functional β V for the whole potential just by looking at the vacuum renormalization. This way of thinking has at least a two fold advantage: it simplifies computations (since we just need to compute the vacuum renormalization) and gives direct access to the full system of beta functions for the couplings of the operators O k = φ k (since β V serves as a generating function for the beta functions). From the knowledge of the beta functions we can then follow the steps outlined in the previous section and compute both the spectrum and the OPE coefficients in the MS scheme. The action (3.1) not only renormalizes the potential, but also induces the flow β Z of a field-dependent wavefunction functional that we will denote Z(φ). The flow generates the beta functions of the couplings of the operators of the form φ k (∂φ) 2 and, moreover, fixes the anomalous dimension η. More generally, all higher derivative operators have an approximate flow induced solely by the potential, i.e. have a beta functional whose r.h.s. contains only V(φ) and its derivatives. We will call local potential approximation (LPA) the truncation for which all the RG flow is generated by the potential alone. Clearly, the full RG flow will involve the presence of other functionals, such as Z(φ) and higher, on the r.h.s. of the beta functionals. According to that the computational scheme can be systematically improved in a derivative expansion approach, as will be discussed in section 4.
In this section we will study, as a tutorial example, the Ising and Lee-Yang universality classes within the LPA, while a first example of functional flow beyond this approximation will be presented in section 5, in which we show that terms containing Z(φ) on the r.h.s. of the beta functional β Z become important to describe mixing effects when marginal or irrelevant operators are investigated in the Ising and multicritical universality classes.
For any given theory and within a functional perturbative approach in a dimensionally regularized scheme, e.g. MS, such beta functionals can be written as polynomials for which each monomial is a product of derivatives of various orders of the generating functions V, Z, . . . and in particular each non trivial loop order in perturbation theory gives rise to a subset of monomials in the beta functionals [7]. Let us just stress that for a given theory the monomials which can appear in the beta functional are very constrained and only their coefficients demand a real loop computation, which in turn can be done in very specific and simple ways.
Another point to highlight is that, depending on the specific theory, contributions denoted as LO (or NLO or higher) appear at different number of loops, generally bigger than one (Ising and Lee-Yang are special in this respect since the LO terms are obtained at one loop and NLO at two loops for both theories). We shall see this explicitly in section 5 in which we study the whole family of multicritical φ 2n universality classes: indeed the number of loops required to obtain the LO contribution depends on the critical dimension d c of the theory (which determines the superficial degree of divergence of a diagrams generated by the perturbative expansion). One sees that, since each member of this family of models has d c = 2n n−1 (for n > 1), the leading order contribution appears at (n − 1)-loops and the NLO at 2(n − 1)-loops. Let us also mention the fact that the LO and NLO order contributions are universal, i.e. independent of the specific RG scheme (as can be easily seen by projecting the beta functionals on the beta function of the respective critical coupling which we already know has LO and NLO universal coefficients [33]).
It is also convenient to make the standard shift to dimensionless variables (in units of the scale µ) directly at the functional level. Once the beta functional of the dimensionful potential is found, the scaling properties are investigated by defining the dimensionless where Z 0 is the field strength renormalization, which enters in the definition of the dimensionless field ϕ = µ 1−d/2 Z 1/2 0 φ. Its beta functional is then for which we introduced an anomalous dimension η = −µ∂ µ log Z 0 , which will be discussed in detail soon. The potential is a local function of the dimensionless field ϕ and can be parametrized in terms of the dimensionless couplings g k as The beta functional is then used to obtain the couplings' beta functions through the straightforward definition One then inserts (3.4) and (3.5) on the r.h.s. and l.h.s. of (3.3), respectively, and equates powers of the field on both sides to obtain the general beta function system.
The dimensionless wavefunction is similarly defined as and its dimensionless flow is This new beta functional has two main purposes. On the one hand by enforcing the condition z(0) = 1 we can use it to determine η as On the other hand, later in section 5 we will use (3.6) to generate the beta functions of the dimensionless couplings of the operators of the form ϕ k (∂ϕ) 2 . A detailed discussion of the invariance of the systems of beta functionals in the LPA under reparametrizations of z(0), and of its importance in the determination of η can be found in [37].

Ising universality class in LPA
The Ising universality class has upper critical dimension d c = 4, and the LPA beta functionals for the dimensionful potential at NLO, and wavefunction at LO, are The functional form of these beta functionals can be argued on dimensional grounds. Only the explicit determination of the three universal coefficients demand a loop computation, but for a well studied universality class such as Ising these coefficients can be obtained by matching with known beta functions of the φ 4 critical coupling.
In turn, this simple observation shows that these coefficients are scheme independent by the standard text book argument that LO and NLO beta functions and anomalous dimension coefficients are so. Thus all the results of the present section, and in particular the form of the beta functions around the fixed point are a functions of these universal numbers. In particular this implies that the order and 2 contributions to the spectrum (which is universal) are scheme independent, and that the order corrections to the OPE coefficients are also scheme independent even if theC i jk themselves are not universal. As promised this is a simple way to understand why the Hessian for the "dimless" OPEs is of order at least O( ) [38].
After a simple rescaling v → (4π) 2 v, the beta functionals for the dimensionless potential are the following (v (4) ) 2 . (3.9) Expanding the potential as in Eq. (3.4) discussed before, allows the generation of the coupling's beta functions (3.10) The four-coupling system (2.12) studied in the previous section is straightforwardly obtained by truncating (3.10) to k = 1, 2, 3, 4. Using (3.7) we can immediately obtain the anomalous dimension In dimensional regularization the fixed point is very simple and highlights the prominent role of the critical coupling g 4 . By expanding around the fixed point it is straightforward to obtain the following general form for the spectrum (in terms of the critical coupling) Using (2.14) we immediately deduce the anomalous dimensions of the composite operators φ i (in terms of ) (3.14) For i = 1, 2, 3, 4 this expression reproduces those of the example in the previous section. The reader will notice the appearance of a contribution in (3.13) and in (3.14) proportional to the Kronecker delta δ i,4 because the anomalous dimension in (3.9) is a function of the critical coupling as given in (3.11) [18,19]. The expressions for the spectrum are complete to order O( 2 ) for all relevant couplings and, as we will show in section 5, also for the marginal ones. For irrelevant couplings, due to mixing effects, only the O( ) terms are complete and correctly agree with the CFT results [9]. From the analysis of the quadratic part of the beta function we find the following form for the universal OPE coefficients in the MS schemẽ This general expression gives us back the results of our previous example (2.16) and thus matches, when overlapping, with CFT computations [9], but its general range of validity will become clearer in section 5 after we analyze the effects of mixing. We note that a NNLO computation, beside bringing some mixing effects, will provide further contributions at order O( 2 ) so that one should consider at this level of accuracy the expressions (3.15) just up to order O( ) as in Eq. (2.16). In fact, recalling our discussion in section 2.3 on the possible differences among coefficients computed in other schemes, agreement at order O( 2 ) with an NNLO computation can be observed only if the MS scheme and the "CFT" scheme are related by a Hessian of order O( 2 ).

Lee-Yang universality class in LPA
The Lee-Yang universality class has upper critical dimension d c = 6 and the LPA beta functionals at NLO for the dimensionful potential and wavefunction are ( 3.16) The explicit derivation of these beta functionals is quite straightforward. After the convenient rescaling of the potential v → 2(4π) 3/2 v the beta functionals for dimensionless quantities are Expanding the potential as Eq. (3.4) leads to the general expression for the beta functions where I = i + j + t + u and J = i 2 + j 2 + t 2 + u 2 . From (3.17) we also immediately obtain the anomalous dimension η = 24g 2 3 + 1248g 4 3 .
As expected, because of the non-unitarity of the model, the fixed point is complex (3.20) showing that the -expansion for the Lee-Yang universality class is in fact an expansion in powers of √ − ; equivalently one can write g 2 = − 1 54 − 125 . After expanding the beta functions around the fixed point we determine the spectrum in terms of the critical coupling from which using (2.14) we can extract the anomalous dimensions, as a function of and taking in account (3.20) For reference we write the first anomalous dimensions It is easy to check that the scaling relation θ 1 + θ 2 = d, discussed in appendix B, is indeed satisfied.
For the universal MS OPE coefficients we obtaiñ which at this order we display as a function of the coupling g of (3.20) for notational simplicity. Using the explicit for of the fixed point (3.20) as a function of , and considering only the leading order in √ − we find (3.24) The first two universal OPE coefficients arẽ and agree with CFT computations [9]. The discussion of the universality of the Ising's OPE coefficients has an analog here: In the case of the Lee-Yang universality class we have that (θ c − θ a − θ b ) = O( √ ), therefore the eventual Hessian relating the MS and CFT schemes might contribute by changing the universal OPE coefficients at O( ) or higher.

FUNCTIONAL PERTURBATIVE RG AND THE DERIVATIVE EXPANSION
In the rest of the paper we would like to show how it is possible to generalize the results presented so far to arbitrary order in the -expansion to include mixing effects, and also to extend the analysis to a wider set of universality classes. In order to enter into this subject and also pave the way for future computations including the most general operators, we will first describe the general setup of the derivative expansion where one can systematically include higher-derivative operators. In what follows, our aim would be to outline a systematic approach to such a derivative expansion in the functional perturbative RG. The derivative expansion, although being formally a truncation of the most general action (3.1), allows, when combined with the perturbative -expansion, a systematic and complete determination of the -series of the spectrum and the MS OPE coefficients.
At each order in the number of derivatives there is an infinite number of operators with higher and higher powers of the field. Just like the potential function V(φ) which encompasses an infinite set of couplings, the couplings of these derivative operators can be collected into functions so that at each derivative order there is a finite basis of "functional" operators which spans all the operators with the given number of derivatives. To make it more explicit, one can denote the basis of functional operators with k derivatives byŴ (k) a (φ), where a runs from 1 to N k , the number of elements in such a basis. With this notation the action (2.1) can be re-expressed as where the index k runs over the number of derivatives and a spans the possible degeneracy.
The first few instances of such operators can be listed as follows: and we will adopt the notation W a ≡ W a . One can continue in this way and choose a basis for higher-derivative operators. At the next order, i.e. six derivatives, there are N 6 = 8 independent functional operators which form a basis. This will increase to N 8 = 23 for the case of eight derivatives, and so on. For each of the operators in (4.1), after shifting to the relative dimensionless functionals w At this stage, let us be more specific and concentrate on theories of the form (4.1) close to the upper critical dimensions of multicritical φ 2n models [9] Then dimensional regularization has the virtue that the fixed-point action solving β (k) a = 0 is extremely simple, since all fixed point functionals are zero apart from the potential V, which in turn is proportional to the critical coupling, and the Z, which is constant and can be set to one. The fixed point action is then These choices define the multi-critical universality classes which will be discussed in detail in the next section. The action (4.1) can then be seen as a deformation around the fixed point (4.3) away from criticality. One can formally define the stability matrix and the set of OPE coefficients in a functional form by expanding the beta functionals around the fixed point as 9 In the above expression of the beta for the functions w k a , which depend non linearly also on derivative of them, one has formally functional derivatives and integral are understood when repeated indices a, b, · · · are present. Although one can study the RG flow and compute all universal quantities directly at the functional level by exploring the consequences of (4.5), in the next section we will reconnect with the discussion in terms of couplings as outlined in section 2, and use the beta functionals β (k) a as a convenient way to generate the coupling beta functions.
The couplings in (4.1) can be defined by expanding the functions such as V(φ), Z(φ), W a (φ) and those of the higher derivative operators, in powers of the field, starting with φ 0 = 1. Using dimensional analysis and recalling that close to the upper critical dimension the spectrum of the theory is almost Gaussian, we can infer that the couplings in V(φ) corresponding to the 2n lowest dimensional operators 1, φ, · · · , φ 2n−1 do not mix with any other coupling. Staring from φ 2n and all the way up to φ 4n−3 they mix with the O(∂ 2 ) couplings of (∂φ) 2 , · · · , φ 2n−3 (∂φ) 2 . From φ 4n−2 , φ 2n−2 (∂φ) 2 the O(∂ 4 ) couplings of W a (φ) will also be involved. This can be summarized in the following table 9 More generally, for an arbitrary Lagrangian L, the RG flow can be formally described by a beta func- V : where each row collects the operators included in the function shown on the left-hand side and only couplings of operators in the same column mix together. If we arrange the couplings of (4.1) in increasing order of their canonical operator dimension, and furthermore, we sort them for increasing order of derivatives of their corresponding operators, the stability matrix takes the block-diagonal form where in general M (2k) is itself a block diagonal matrix. Each diagonal block contained in M (2k) describes the mixing between couplings of operators up to 2k derivatives, all of which belong to the same column in (4.6). In particular M (0) is a diagonal matrix with entries giving the scaling dimensions of the first 2n couplings in the potential. The matrix M (2) is block diagonal, with each block being a two by two matrix which gives the mixing between a coupling in V(φ) and a coupling in Z(φ). M (4) is also a block diagonal matrix of which, with our choice of basis for the four-derivative operators, the first block is a three by three matrix, the second is four by four and the rest are five by five matrices. Using dimensional analysis one can restrict the stability matrix even further if one is satisfied with the order approximation. As we will show explicitly in the following section, at this order each diagonal block in the matrices M (2k) in (4.7) is itself block lower-triangular, where each block describes the mixing of operators with the same number of derivatives. This ensures that the entries on the diagonal for the couplings of the potential and the second-derivative operators will give the scaling dimensions and are unaffected by the mixing at this order.

GENERAL φ 2n UNIVERSALITY CLASS
After the analysis of the Ising universality class in sections 2.2 and 3.1 and the introductory discussion of the previous sections, we are now in a position to extend these results to general models with even interaction φ 2n at the critical point. In fact one can treat the whole set of universality classes φ 2n in a unified framework. A brief review of the method which closely follows [7] is outlined in appendix A. Here we pick the main results that will be needed for our analysis.
Throughout this work we will not go beyond second order in the derivative expansion, and in fact mostly concentrate on the local potential approximation. Let us therefore consider a theory of the form in a space-time dimension which is close, as in Eq. (4.2), to the upper critical dimension at which the coupling of the interaction φ 2n becomes dimensionless. The propagator of this theory satisfies the differential equation x is the d-dimensional Dirac delta function. The solution is given by which is more conveniently written in terms of the field dimension δ = d 2 − 1 = 1 n−1 − 2 . The coefficient in the propagator evaluated at criticality appears many times in the calculations. For convenience we therefore call it c from now on Let us neglect for the moment the effect of derivative interactions encoded in Z(φ) and concentrate on the V(φ) contributions to the beta functions of V(φ) and Z(φ). Before giving the explicit expressions for the beta functions let us mention that in this O(∂ 0 ) truncation one can extract the scaling dimensions and MS OPE coefficients for the relevant components as they will be in any case unaffected by the mixing with the derivative operators. Moreover, as we will argue later, remaining within the same truncation it is possible to go beyond the relevant components if one is content with the order estimates. Neglecting derivative interactions, the beta functional of the dimensionless potential, in the form of Eq. (3.3), at the NLO (cubic order) in the dimensionless potential is where the integers r, s, t are implicitly taken to be positive, and the quantities K n rst and L n st are defined as follows where ψ(x) = Γ (x)/Γ(x) is the digamma function. The last term in the first line of (5.4) is the LO (n − 1)-loop term, while the NLO second and third lines appear at 2(n − 1)loops. The origin of such terms and the corresponding diagrams will be briefly discussed in appendix A. Notice also that, differently from sections 2 and 3, we did not yet include any further rescaling when moving from the dimensionful V(φ) to the dimensionless v(ϕ) potential: since the rescaling does not affect the spectrum, we postpone the discussion of the "appropriate" rescaling to subsection 5.3 in which some MS OPE coefficients are computed. Neglecting derivative interactions (in agreement with our definition of LPA), the induced flow of the function z(ϕ) at quadratic order is given by The last term in this equation comes from a diagram with 2(n − 1)-loops, which gives a counter-term consisting of the second contribution in Eq. (A.2), as explained in appendix A. From (5.4), noticing the fact that only the dimensionless coupling can take a non-zero value at the fixed point, one can set v(ϕ) = g ϕ 2n together with the condition β v = 0 to find the critical coupling g at quadratic order in . This is given by (5.7) Notice that here we have used an expansion of v(ϕ) without factorials. Including the factorials, the (2n)! on the left-hand side would have appeared with power one, in agreement with [7,9]. This, of course, does not affect the final physical results when written in terms of .
The anomalous dimension can be read off from (5.6) imposing the condition β z | ϕ=0 = 0 and using z(0) = 1. This gives, after using (5.7), in agreement with [7,27] and recent CFT based computations [9,10]. Having at our disposal the functional form of β v at cubic order, we can follow the prescription of section 2 to find the scaling dimensions of the relevant couplings at O( 2 ) and the MS OPE coefficients for the relevant operators at O( ). However, before doing so let us devote the next subsection to considering the leading order mixing effects due to the presence of z(ϕ)-interactions. This, for instance, will allow us to compute the leading order anomalous scaling dimensions of the z(ϕ) couplings, and justify the validity of the leading order anomalous dimensions of the v(ϕ) couplings.

Mixing
One can take into account the mixing effects due to the presence of derivative operators. We have given a general sketch of the mixing pattern in section 4. Here we concentrate on explicit results for two-derivative interactions collected in the function z(ϕ). At quadratic level the presence of z(ϕ) does not affect the beta function of the potential (5.4), but (5.6) instead gets a contribution at this level These new terms arise from the derivative interactions and appear at (n − 1)-loops [7] as discussed in appendix A. This gives rise to a mixing, at order , between the operators φ k+2n and φ k (∂φ) 2 for k = 0, · · · , 2n − 3, described by the k-th block of the matrix M (2) , while the couplings of φ k with k = 0, · · · , 2n − 1, and therefore the elements of M (0) , are unaffected.
In general the terms in a beta function which contribute to the stability matrix at order must be quadratic in the couplings, and furthermore one of the couplings must be the dimensionless coupling g which is the only one that takes a nonzero value at the fixed point. In the beta functional this manifests as the product of the potential v(ϕ) and a function corresponding to a higher derivative operator, or more precisely, the product of a derivative of these functions. These terms come from diagrams of the form displayed in Fig. 2

of appendix A, or its generalizations where instead of z(ϕ) one can have functions encoding higher derivative interactions.
A simple argument based on dimensional analysis shows that generally in the beta function of a 2k-derivative coupling, the quadratic term which includes the coupling g can involve also a derivative coupling lower or equal to k. To show this, one should notice that in dimensional regularization the diagrams contributing to the beta functions must be dimensionless, i.e. have vanishing superficial degree of divergence. For a melon diagram of the form in Fig. 2 with r propagators that includes the potential v(ϕ) at one vertex and a 2l-derivative coupling on the other vertex this condition is The first term comes from the r − 1 loop integrations, while the r propagators give a contribution −2r in the second term. To justify the remaining terms one should notice that there are altogether 2l derivatives at one of the vertices, some of which might act on the propagators and some might not, but finally we would like to extract the ∂ 2k contribution from this diagram. This leads to the contribution 2(k − l). This simple relation can be re-arranged and put in a more useful form, where one can use the fact that r ≥ 2 on the right-hand side to put an upper bound on l − k. Since k, l are integers, the inequality (5.11) says that l ≤ k, which is the statement claimed above. This is a more general case of what we have already seen: that the beta functions of the potential couplings do not contain the product of v and z couplings but only v-coupling squared. This is telling us that at order each diagonal block in the stability matrix that describes the mixing of a column in (4.6) is itself block lowertriangular (where here a block describes the mixing of operators with the same number of derivatives), as will be shown explicitly in the simplest case in the next section. This ensures that at order the eigenvalues of the potential and the two-derivative couplings are never affected by the mixing. In particular from (5.4) one can find the spectrum of the couplings in the potential, not only for the relevant ones but also for the marginal and all the irrelevant couplings. Similarly, the beta function (5.6) with the correction (5.9) gives the spectrum of all the z-couplings, at order . These are made more explicit in the following subsection.

Spectrum
In order to proceed with explicit results let us stick to the following convention throughout this section for the expansions of the functions v(ϕ), z(ϕ) and of the corresponding beta functionals in powers of the field The choice of normalization for the couplings is of course physically irrelevant. An explicit computation using the beta function (5.4) and the expansions (5.12) shows that the matrix ∂β i v /∂g j evaluated at the fixed point, which is for dimensional reasons diagonal, has the elements −θ i = −d + i(d − 2)/2 +γ i on its diagonal, with the following anomalous parts 10 For the relevant components, that is for the range 0 ≤ i ≤ 2n − 1, these are simply the anomalous dimensions with accuracy O( 2 ). The last term in the first line, which comes from the term proportional to η in (5.4), does not contribute in the relevant sector. However, if one wishes to find the anomalous dimension of the marginal coupling, one has to take this term into account. Within the same O( 2 ) accuracy, for the irrelevant couplings, which we do not consider here, additional mixing transformations are required to diagonalize the stability matrix. From (5.13) one can readily see that for i = 1 all the terms except the first vanish. Also, interestingly, for i = 2n − 1 which corresponds to the descendant operator φ 2n−1 in the interacting theory because of the Schwinger-Dyson equations, the O( 2 ) terms in the second and third line of Eq. (5.13) vanish so that these anomalous dimensions take the simple formγ (5.14) The two anomalous dimensions then sum up toγ 1 +γ 2n−1 = (n − 1) , which is equivalent to the scaling relation θ 1 + θ 2n−1 = d and proved in general in appendix B. The correction (5.9) allows us to go beyond the local potential approximation and compute at order the block M (2) in (4.7) which is a block-diagonal matrix with two by two blocks. The i-th block which gives the mixing of the φ i+2n and φ i (∂φ) 2 couplings is given in the {φ i+2n , φ i (∂φ) 2 } basis as where 1 is the two dimensional identity matrix. For each i the two eigenoperators have the same canonical scaling at the critical dimension. The eigenvalues of the stability matrix include the scaling dimensions −θ i+2n , given in Eq. (2.14), and ( d 2 − 1)i +ω i , which is the analog for z-couplings in the notation of [7].
From these, one can then read off the anomalous partsγ i andω i of the v and z coupling scaling dimensions at order which are valid not only for 0 ≤ i ≤ 2n − 3 described by the above matrix but for all i, according to the discussion in the previous subsection. In summary, again interpreting the factorials to be infinite for negative integer arguments, and for i ≥ 0 This reproduces the result of [7]. Theγ i in Eq. (5.16) also match the anomalous dimensions found in [9,10] from CFT constraints. Beyond the leading order for the anomalous dimensions, the stability matrix will not be lower-triangular anymore, and in order to find the anomalous dimensions of higher and higher powers of φ one has to take into account (up to cubic order contributions of) operators of higher and higher dimensions. In the simplest case, 2n < i < 4n − 3, one needs to include cubic corrections to β z , and furthermore, take into account the z(ϕ) contribution to β v at cubic level. The only term contributing to this last piece is proportional to v (n) (φ) 2 z(φ) and leads to O( 2 ) corrections in the upper right element in (5.15). These higher order corrections are not considered here and are left for future work.
Besides (5.16), an extra information which has been obtained in [9] using conformal symmetry and the Schwinger-Dyson equations is the leading order value of γ 2 for n > 2, which is of order 2 . For n > 2, putting i = 2 in (5.13) gives which is also in agreement with the result found in [9].

OPE coefficients
The only non-zeroC k ij coefficients that are extracted from the beta functions are those that are massless, or equivalently, satisfy the universality condition i + j − k = 2n. Contrary to the anomalous dimensions, the OPE coefficients do depend on the normalization of the couplings. Throughout this section we continue to use the normalization where couplings appear without factorials in the v(ϕ), z(ϕ) expansions, as defined in (5.12). On top of this, it turns out convenient to make a global rescaling of the couplings by redefining the potential according to 11 v → 4 (n − 1)c n−1 v . (5.18) This removes the parameter c from the beta functions (5.4) and (5.6). In such a normalization, using the beta function (5.4), the expansion of the potential and its beta functional in powers of the field (5.12), and the fixed point relation (5.7), a lengthy but straightforward calculation based on the definition (2.5) gives the MS OPE coefficients (k = i + j − 2n) Notice that, strictly speaking, the above quantity is in fact the matrix N k ij defined in (2.5), but because the mixing matrix S i a is diagonal (in the relevant and marginally irrelevant part of the spectrum) it coincides with the OPE coefficientsC k ij in our scheme. 11 In Sections 3.1 and 3.2 this rescaling was used for both the Ising and the Lee-Yang universality classes with n = 1 and n = 3 2 respectively.
The last contribution in (5.19) comes from the anomalous dimension term in (5.4). Similar to the anomalous dimensions (5.13) one has to keep in mind that terms with negative factorials in the denominators vanish. Notice that the first term is nothing but the combinatorial factor that comes from Wick contractions in the free theory. The normalization we have adopted therefore coincides with the CFT normalization where the coefficient of the two point function φ φ is set to unity.
It is important to comment on the range of validity for the i, j indices in the above formula. As in the case of anomalous dimensions, Eq. (5.19) is, of course, valid for all relevant components, that is, positive integer indices smaller than 2n. Notice that in this case the last term does not contribute. However, this is not all we can extract from this formula. For instance, let us consider the case i < n. The above formula will then be of order . At this level of approximation Eq. (5.19) is valid for any j, and not only the relevant ones. This is because mixing effects enter only at NLO. Therefore for such cases one can use (5.4) without any concern about the mixing. Notice that for these cases only the second term on the first line and the last term in Eq. (5.19) contribute. The particular caseC 1 1,2n gets contribution only from the last term in (5.19) and takes the simple form This reproduces the result found in [9] from CFT considerations, and therefore the computation done in the MS scheme reproduces an entire family of CFT OPE coefficients at least at order O( ). Finally, let us consider Eq. (5.19) for k = 1. The indices i, j must then satisfy i + j = 2n − 1, so we choose them as i = n − m, j = n + m + 1, for m = 1, . . . , n − 1. The OPE coefficients reduce tõ This MS result is also in agreement with [9], and with [10] if one takes into account the different normalizations of the operators φ l .

Examples of CFT data for specific theories
Despite the above general treatment being comprehensive of all the even multicritical models, we believe it is interesting to show some explicit results for specific theories. The case of Ising had already been studied in subsections 2.2 and 3.1. In this subsection we collect the CFT data for the Tricritical and Tetracritical universality classes.
The Tricritical universality class corresponds to n = 3. The anomalous dimensions for the relevant and marginal operators at O( 2 ) can be obtained from the general formula (5.13) and are given explicitly as (5.22) The scaling relationγ 1 +γ 5 = 2 is satisfied. Notice that restricting to order one can immediately extend these results to all the couplings including the irrelevant ones, and even further to the z(ϕ) couplings. These were reported in Eqs. (5.16). It is also easy to extract from the general equation ( The OPE coefficientsC 1 25 andC 1 16 exactly match the corresponding structure constants computed with CFT methods in [9]. For the others there are no available CFT results to compare with.
For the Tetracritical universality class, which corresponds to n = 4, there are seven relevant couplings whose anomalous dimensions arẽ γ 1 = 9 2 171500 γ 2 = 27 2 171500 , (5.26) while the anomalous dimension of the marginal coupling is given as For this universality class the OPE coefficientsC 1 27 ,C 1 36 andC 1 18 correctly match the corresponding structure constants computed in [9]. The considerations made for the Tricritical case on the full comparison to CFT are equally valid for the Tetracritical universality class.

CONCLUSIONS
In this paper we have shown how to extend renormalization group (RG) techniques to the computation of some OPE coefficients at a scale invariant critical points of scalar quantum field theories. The approach of this work employs dimensional regularization in the MS scheme at the functional level and gives access to a specific set of "massless" OPE coefficients, which are related to terms in the beta functions that are universal at the upper critical dimensions of the models under investigation. For general multicritical models we have extracted these quantities, and we have shown that at order O( ) they agree with the corresponding OPE coefficients computed directly with CFT methods, when available on both sides.
Let us briefly summarize our procedure. In the vicinity of a fixed point the RG flow can be expanded in powers of the couplings. The information on the universal quantities is encoded in the coefficients of this expansion, in which the linear and the quadratic parts play a special role. The linear terms give rise to the so-called critical exponents which are related to the scaling properties of the operators of the theory and have been the focus of most RG studies so far. The quadratic terms instead give information on some OPE coefficients, which can thus be extracted from the general knowledge of the beta functions. Whenever scale invariance implies conformal invariance, the OPE coefficients are directly related to the structure constants of the underlying CFT and thus our analysis strengthens the link between RG and CFT by showing explicitly how, and which part of, the CFT data can be determined to some extent directly within an RG approach.
The scheme dependence of the results can be analyzed in terms of the coupling redefinition connecting two different schemes. It follows that the spectrum is invariant and in this sense universal, while the coefficients of the quadratic terms generally transform inhomogeneously under coupling redefinitions. In a dimensionally regularized MS scheme one has access only to the "massless" quadratic coefficients which are universal at d = d c , but potentially differ in other schemes at higher orders in the -expansion. A first observation is that our computational scheme gives the correct values at order O( ).
After a first pedagogical application of the approach to the investigation of the Ising universality class, we have introduced a very convenient functional generalization of the standard perturbative RG, in which all the beta functions for the couplings are obtained from few simple generating functions: the beta functionals. This functional perturbative framework is a very useful tool that naturally organizes the beta functions in simple generating functionals with few independent c-number coefficients, which, we stress, at leading and next-to-leading order are RG-scheme independent. As a result all the quantities we have computed (anomalous dimensions at order O( 2 ) and OPE coefficients at order O( )) depend essentially only on these universal coefficients.
The simplest of these generating functionals is β V , which encodes the RG flow of the whole potential V(φ) and thus of all couplings of operators of the form φ k . Contributions of operators involving more derivatives can be included systematically. The first such contribution comes from β Z , which generates all beta functions of the O(∂ 2 ) operators φ k (∂φ) 2 included in the field dependent wavefunction Z(φ). A goal of our work has been to emphasize some of the advantages of this shift towards a functional approach to standard perturbation theory in the -expansion, because it grants an easy and systematic determination of important universal quantities like both the scaling dimensions and expressions for the massless MS OPE coefficients.
As a first application of the functional perturbative RG we have reconsidered the Ising and the Lee-Yang universality classes, as representative of the multicritical unitary and non unitary families, at the level of the local potential approximation (LPA), i.e. without taking into account derivative interactions, and we have showed how the results extracted from the RG coincide with those recently obtained with CFT techniques. This we take as evidence that the MS scheme is effective in the computations of the leading -corrections to some of the OPE coefficients.
We have also outlined a systematic approach to the inclusion of higher derivative interactions in functional perturbative RG, and discussed the general mixing patterns among operators at different orders in the derivative expansion using the general φ 2n models as examples. The efficiency of the functional RG techniques is most clearly seen in this context, since we are able to collect infinite towers of critical exponents and OPE coefficients at order O( ) in compact formulas. This in particular has allowed a straightforward check of our MS estimates with the results obtained recently with CFT techniques [9,10]. We stress that for the general φ 2n our approach is a multi-loop analysis which, for almost all models, is characterized by an -expansion below a fractional critical dimension [39].
The computation of the anomalous dimensions for the multi-critical models in a functional framework was previously carried out by O'Dwyer and Osborn [7] and here we have limited ourself essentially to the same order of the perturbative expansion. The analysis is first done without taking into account derivative interactions, and afterwards including the leading order of the mixing with O(∂ 2 ) derivative interactions. Using dimensional analysis we have imposed further constraints on the stability matrix by determining possible terms that can appear at the quadratic level in any beta function. In particular we have shown that at order the stability matrix is lower-triangular. This allows, in agreement with the CFT analysis, the determination of the anomalous dimensions for all operators contained in the potential V(φ) and the wavefunction Z(φ) up to order , and therefore is not limited to the relevant operators. We have also given the O( 2 ) results for the scaling dimensions up to the marginally irrelevant operator φ 2n , and provided the O( ) values for a family of infinite "massless" OPE coefficients for all the multi-critical universality classes considered, which are found to match the corresponding values from CFT analysis, when available. Clearly more investigations on universality and scheme dependence of the results presented here are required. Here we just observe that the functional framework in the context of perturbation theoryexpansion strongly constrains the possible redefinitions of the couplings giving support for the success of our approach. We shall discuss in detail these issues in a forthcoming paper [38]. On the other hand this fact can be seen as another argument in favor of adopting a functional approach to RG analysis. To the best of our knowledge, very few results have been obtained even at order O( ) in CFT computations. It would be especially interesting to have CFT results at order O( 2 ) whose comparison with the RG NNLO estimates would show possible artifacts induced by MS scheme.
Summarizing, the main results of our paper are highlighted as follows: Inspired by the analysis presented by Cardy [2] which relies on an ultraviolet cutoff, we have pro-posed to extract the OPE coefficients of a CFT from the coefficients of the quadratic terms in the coupling expansion of the beta functions around the fixed point using dimensional regularization and MS scheme at functional level. We have discussed the scheme dependence of OPE coefficients obtained in this way and identified those that are less sensitive to changes of scheme, which turn out to be the ones that are dimensionless at the upper critical dimension. The order OPE coefficients that we have found are compared with the literature on CFT approaches and when available with both methods it is shown that the results always agree. This analysis is done for all multicritical even models, for which the beta functions were obtained in [7], as well as the Lee-Yang model as a representative of the odd multicritical models for which we have reported the functional betas at NLO (3.16). We have demonstrated the power of the functional approach by obtaining compact formulas (3.23), (5.19) encompassing an infinite number of OPE coefficients. Similar formulas are obtained for the order 2 critical exponents of the relevant operators (3.22), (5.13) 12 which, for instance, allow us to verify the shadow relations (5.14) for all such models in one shot. Finally we have argued that dimensional analysis alone constrains the structure of the beta functions and in particular the stability matrix. We have used this information to prove that the formulas for the order critical exponents (5.16) are valid for all the V and Z couplings.
The functional perturbative RG as introduced in this paper is very general and can be systematically pushed to higher levels of accuracy by including new families of operators, a fact which is made particularly evident by working at the functional level. We have recently also successfully applied this method to study [40] the non unitary family of multi-critical universality classes described by single scalar field models with odd potentials, whose first elements are the Lee-Yang and tricritical Lee-Yang (Blume-Capel [26,41]), for which some CFT results are already available [9] finding again full agreement. We plan to further develop the main ideas and apply the method to other universality classes, e.g. for multifield cases, as well as to carry on investigations at higher orders in perturbation theory.
Another extremely important line of investigation, which could possibly overcome the limitations of the perturbative approach, is to move to one of the non-perturbative functional RG frameworks [28][29][30][31]. This step is absolutely non trivial because such approaches are based on massive renormalization schemes, which often result into a much stronger deformation of the basis of scaling operators (as compared to the Gaussian basis), and make it difficult to establish a direct link to the CFT results. We leave this line of investigation to future research.
FIG. 2: Diagrams contributing to the counter-term of the Z function at quadratic level in the couplings to the potential the number of propagators must be constrained to r + s + t = 2n. The second one consists of two melon diagrams as in Fig. 3b, and the third graph, shown in Fig. 3c, is a melon diagram involving the potential and its counter-term at quadratic level V c.t. (φ), which is the first term on the right-hand side of (A.2). In both diagrams the singularity that contributes to the potential occurs when the number of propagators in each melon is equal to n. These three diagrams are therefore all of 2(n − 1)-loop order. They give rise to the cubic terms in the second and third lines of (5.4).
The precise relation between the counter-terms and the dimensionful beta functions of the potential at quadratic level β V,2 and at cubic level β V,3 and also the dimensionful beta of the wavefunction at quadratic level β Z are given by the following equations The total beta of the potential β V = β V,2 + β V,3 and the beta of the wavefunction β Z are related to the dimensionless betas through (3.3) and (3.6), respectively. The cubic counter-term V c.t.3 (φ) is the sum of the diagrams (3a), (3b) and (3c), and the quadratic counter-terms for the wavefunction Z v 2 c.t.2 and Z vz c.t.2 are extracted respectively from the second term on the r.h.s of (A.2) and from the counter-term diagram of Fig. (2) for r = n. The µ-derivatives with an index 1 are taken using the tree-level flow, which for the derivatives of the potential and the wavefunction are given by the following relations while the µ-derivative with an index 2 is based on the quadratic flow. In particular

Appendix B: A general scaling relation
In this appendix we would like to obtain a relation valid among the scaling of two couplings induced by the RG flow. This information can then be compared to the relation obtained in CFT for the scaling of the field operator and one of its descendants.
We have already encountered the scaling dimensions of the operator φ i and its corresponding dimensionless coupling g i , which were denoted by ∆ i and θ i in Eqs. (2.15) and (2.14), together with their anomalous parts γ i andγ i respectively. Let us consider for a moment the case of a multi-critical theory φ 2n . Then, for i = 2n − 1, the relation θ i + ∆ i = d holds. This is equivalent to γ i =γ i . Instead, for the descendant operator corresponding to i = 2n − 1 this relation is modified to θ i + ∆ i = d + η by the presence of η = 2γ 1 = 2γ 1 , which is twice the anomalous dimension of φ. One can link this fact to the relation γ 2n−1 = (n − 1) + γ 1 coming from the descendant constraint in CFT, ∆ 2n−1 = 2 + ∆ 1 and from another relation that we shall prove in general in the following. Indeed we shall see that the latter is equivalent toγ 2n−1 +γ 1 = (n − 1) so that the two anomalous dimensions (associated to the CFT operator and RG coupling) are related by γ 2n−1 =γ 2n−1 + η.
We shall work at a general functional level [28]. Let us consider for the truncation with two functions V and Z which describes deformations with composite (non-total derivative) operators containing up to two derivatives. The beta functions describing the RG flow are generically written as in Eqs.
For any given dimensionality, we shall refer to the terms of this expansion as critical if they scale as k 0 and off-critical if they do not. For example the term c d k d−4 (V ) 2 is critical in d = 4, while all other terms are off-critical. The critical terms have two important properties: On the one hand they are independent by the cutoff; this is because once the momentum scale q 2 is integrated out, the scale k is what remains of R k (q 2 ), so independence of k implies independence of the cutoff function itself (this of course can be proven more rigorously). On the other hand they are related to the logarithmic divergences of the theory; using again d = 4 as an example which also implies that they correspond to the 1 poles of dimensionally regulated perturbation theory. It is instructive to choose a procedure that deliberately removes the off-critical terms from the flow (C.2). We obtain It is easy to see that the above results correspond to the leading one loop contributions of the two tutorial examples Ising and Lee-Yang. These two examples are the only two universality classes that are captured through critical terms by the above procedure, even though (C.2) is well known to be able to "see" critical points corresponding to all the φ 2n models [42] and more [26]. The reason why only those two critical terms appear has to do with the fact that a local potential truncation of the operator space of Γ k does not contain all possible terms that can be generated perturbatively by higher loops. This should also explain why (C.2) returns only the leading terms of the Ising and Lee-Yang universality classes. The study of truncations that include the higher loops effects has been initiated in [43], in which also the scheme-dependence of functional renormalization group is carefully investigated, but those results have not yet been formulated in a fully functional form as in the models of the present paper. The careful reader must have noticed that the second term of (C.3) is critical for any even value of d. In d = 2 the critical model corresponds to the Sine-Gordon universality class. The beta function of the dimensionful potential is It is interesting to investigate explicitly the flow of the dimensionless potential in d = 2, which is The above beta function does not contain the scaling term contributed by the field ϕ because the field is canonically dimensionless in d = 2 and fluctuations do not generate a nonzero anomalous dimension. Interestingly, the fixed point solution of the Sine-Gordon universality can be obtained by directly integrating the right hand side of (C.8).
Using v (0) = σ as boundary condition we obtain v(ϕ) = − σ 8π cos( √ 8πϕ) , (C.9) in which we can recognize the well-known Coleman phase √ 8π. This fact is quite amazing since the Coleman phase is a non-perturbative result, which we just obtained on the basis of a perturbative approximation. We plan to return to the study of the Sine-Gordon universality class and of all other universal terms in a future work.
As mentioned above, the method presented in this appendix is limited to the universal terms which come from one-loop diagrams because of the local potential truncation. The truncation of this appendix is by definition unable of dealing with higher derivative operators, or operators which are generally generated beyond the first loop. Furthermore, we have made a specific choice of the cutoff which forces us to resort to the rather brute force method of "chopping" all nonzero powers of the cutoff scale k to locate universal terms. A more refined approach to both these shortcomings which also aligns with our discussion of the scheme transformations of Sect. 2.3 can be found in [36] where special "normal" coordinates in the space of all couplings are found in the context of the functional renormalization group (using the Polchinski equation instead of the Wetterich equation, but arguably the conclusions are very similar). The normal coordinates of [36] could be understood as a geometrical generalization of the basis of couplings with well-behaved scaling properties introduced in Sect. 2.3 and their application clearly shows that a consistent renormalization of correlators of all composite operators, thus including in principle all possible OPE coefficients, is possible within the functional renormalization group approach (at least in the vicinity of the Gaussian fixed point). In order to achieve the same results, the functional method presented in the main text of this paper requires the consistent inclusion of higher derivative operators according to their mixing patters as described in Sect. 4.