Renormalization of multicritical scalar models in curved space

We consider the leading order perturbative renormalization of the multicritical $\phi^{2n}$ models and some generalizations in curved space. We pay particular attention to the nonminimal interaction with the scalar curvature $\frac{1}{2}\xi \phi^2 R$ and discuss the emergence of the conformal value of the coupling $\xi$ as the renormalization group fixed point of its beta function at and below the upper critical dimension as a function of $n$. We also examine our results in relation with Kawai and Ninomiya's formulation of two dimensional gravity.


INTRODUCTION
The multicritical scalar models with φ 2n interaction are the simplest and most straightforward generalization of the φ 4 model. Much like the φ 4 field theory captures the critical properties of an universality class of models that includes the ferromagnetic Ising Hamiltonian, the φ 2n field theory can be thought as describing a generalization in which the Ising's spin domains of plus or minus sign are potentially replaced by n distinct vacuum states which become degenerate at the critical temperature.
The renormalization group (RG) flow of the φ 2n models has been explored at length in the literature: perturbatively [1][2][3], nonperturbatively [4][5][6], and with non-canonical kitetic terms [7,8]. It is well-known that a consistent perturbative expansion in the coupling can be constructed at the upper critical dimension d n = 2n n − 1 . (1.1) It is easy to check that, as expected, the case n = 2 corresponds to the φ 4 interaction which has upper critical dimension d = 4 [9]. The model n = 3 corresponds to the φ 6 interaction and is known to describe the universal features of the tricritical Ising model with upper critical dimension d = 3 [10]. All other models have purely fractional upper critical dimensions which asymptotically tend to d = 2. As a consequence d = 2 is the first physical dimension in which all the models φ 2n are nontrivial; the continuation to two dimensions is particularly relevant because they are known to interpolate with the unitary minimal models arising as representations of the infinite dimensional Virasoro algebra [11,12].
For the most part the renormalization of the multicritical models generalizes the one of the φ 4 model, but the leading contributions to the critical exponents are determined by multiloop computations in which the number of loops increases with n [2]. Likewise the φ 4 model, the φ 2n interactions describe critical theories that are Gaussian for d > d n and logarithmic at d = d n , but have non-trivial critical exponents for d < d n . A common practice is to compute such critical exponents in the ǫ-expansion, in which one introduces the constant ǫ = d n − d and uses it to parametrize the displacement of the critical point from the Gaussian theory that occurs at d = d n [9].
The only multicritical model that has nontrivial exponents in d = 3 is the φ 4 one unless one includes the multicritical non-unitary models φ 2n+1 [13], as we shall briefly do later. Specifically, φ 3 and φ 5 have upper critical dimensions d = 6 [14] and d = 10 3 respectively, but they require the tuning of an imaginary-valued magnetic field at criticality [15,16]. It is important to mention that the d = 2 realizations of these models are all "far away" in a perturbative sense from their Gaussian points even though d n → 2 for n → ∞ [2]. Nevertheless, the simple existence of the sequence of multicritical theories provides a very interesting and valuable link between purely field-theoretical realizations and CFT representations [3,17].
One natural and potentially interesting generalization of the above discussion is the study of the renormalization of the φ 2n models in curved space. Generically, the renormalization of a model in a curved background requires additional care to preserve covariance and further conditions to avoid new and unwanted infinities [18]. The extra work is often a gateway to extra information on the theory under consideration [19]. If the multicritical models are coupled with a background geometry, simple dimensional analysis reveals that there is a new non-minimal marginal interaction with the curvature: 1 2 ξφ 2 R. One expects that in curved space the perturbative construction should thus accommodate for some mixing between the φ 2n and φ 2 R operators regardless of n. In other words, the nonminimal interaction φ 2 R holds a special status in that it is always canonically marginal.
A guess on the value that the coupling ξ can take at a curved space generalization of the critical point could be made as follows: Consider a nonminimally coupled "free" scalar field with quadratic action Ideally, the above action captures the Gaussian (non-self-interacting) limit of the φ 2n models which is realized exactly at the upper critical dimension. The nonminimal action is invariant under a conformal Weyl rescaling g ′ µν = Ω 2 (x)g µν and φ ′ (x) = Since conformal invariance implies scale invariance, the nonminimal action (1.2) is thus scale invariant when ξ takes the conformal value, but it is also expected to be a descrip-tion of the critical (scale-invariant) φ 2n model when the interaction becomes Gaussian at the upper critical dimension. Putting everything together we make the following guess.

Educated guess:
The critical point of the coupling ξ at the upper critical dimension, which emerges as fixed point of the renormalization group, is the conformal value (1.3). 1 More generally, one would be tempted to extend the above statement to any dimension below the upper critical dimension, having expressed the desire of analytically continuing these models to d = 2. In this case, we would want to know the conditions under which the conformal value (1.3) is always the critical value for ξ even when the φ 2n interaction is non-Gaussian below d c . For this purpose it is instructive to recall the investigation by Brown and Collins [18], in which it is shown that at the leading order our educated guess is true in the special case of the φ 4 model, but beyond the leading order one has to exploit the freedom of subtracting additional finite parts proportional to the leading counter terms [19]. An analog renormalization condition has also been adopted for the φ 3 model [23] in d = 6, and it plays an important role in preserving conformal invariance in [24,25]. Notice that, strictly speaking, the said two examples have not been concerned with the analytic continuation below the upper critical dimensions d = 4 and d = 6, while our interest is to bring the multicritical models down to d = 2 which does require continuation. However, assuming that we have the same freedom in changing the renormalization condition, we can state a conjecture.
General conjecture: The critical point of the coupling ξ, which perturbatively is determined as an ǫ-expansion in ǫ = d n − d, can always be the conformal value (1.3) thanks to an opportune renormalization condition.
In this paper we consider the leading renormalization and ǫ-expansion of all the infinitely many multicritical models φ 2n (and some other generalizations as well) in curved space using the formalism of functional perturbation theory [2,3]. With the leading results we can show that the conformal value (1.3) of ξ is indeed the critical value at the leading order in ǫ, thus proving the educated guess of this introduction. While we leave the above general conjecture open, we stress that the educated guess is proven for an infinite number of theories. Interestingly, our computation is genuinely new in that the structure of the counterterms and their renormalization does not come from a straightforward generalization of the φ 4 case. On the contrary, we see the case n = 2 as quite the exception which we have to deal with separately.
The paper is organized as follows: In Sect. 2 we study divergences, counterterms, and renormalization group beta functions for all the φ 2n models. We discuss separately the cases n = 2, n = ∞ and the nonunitary models φ 2n+1 . We elaborate briefly on the utility of our results in reproducing some well known formula of 2d gravity in the limit of large central charge. In Sect. 3 we show how the conformal value of ξ emerges as fixed point of its beta function. We also elaborate more on the stronger conjecture expressed in this introduction. Finally in Sect. 4 we draw some conclusion and give a prospect for future investigations. The appendices are dedicated to technical details on the covariant renormalization in curved space. In particular, appendix A discusses the Seeley-de Witt representation of the covariant Green function, and appendix B briefly describes an algorithm by Jack and Osborn for the computation of the poles of dimensionally regulated covariant Feynman diagrams in curved space.

RENORMALIZATION
We are interested in a simple self-interacting canonically normalized scalar field φ which is nonminimally coupled to a background metric g µν in d dimensions. The straightforward bare action is Using the bare action we can formally construct the path integral. For later convenience we shall do it in the background field approach, thus by integrating the fluctuations χ over an arbitrary background φ as follows In flat space it is possible to construct a meaningful perturbative expansion for potentials V(φ) which are polynomials of order 2n below the upper critical dimensions (1.1). If we parametrize V(φ) = λ (2n)! φ 2n + . . . , the upper critical dimension is the one for which the canonical dimension of λ is zero, and the perturbative expansion is controlled by powers of λ itself. Below the upper critical dimensions, these perturbative expansions are known to lead to a sequence of universality classes often referred to as minimal models because they interpolate with the minimal conformal theories arising as representations of the Virasoro algebra in d = 2 [11].
A simple dimensional analysis reveals that if we parametrize F(φ) = ξ 2 φ 2 + . . . , the coupling ξ is always dimensionless and thus it is expected to play a role in the perturbative expansion when promoting the minimal models to curved space. In other words, the φ 2n and φ 2 R operators are expected to mix because of statistical or quantum mechanical fluctuations. For the above reasons we are interested in renormalizing the path integral in d = d n dimensions with V(φ) and F(φ) restricted to be polynomials of order 2n and 2 respectively, so to include all the relevant and naively marginal operators of the models. We do it by adopting dimensional regularization which corresponds to analytically continuing the dimensionality to d = d n − ǫ.
Since the order of the nonminimal interaction is only two, we can incorporate it easily in a quadratic part of the bare action According to the dimensionality, there are two possible leading contributions if the action of the path integral is expanded perturbatively around S 0 [χ] for the φ 2n models in powers of V(φ): the linear and the quadratic contributions. Expanding the path integral to the linear order in V(φ + χ) and Taylor-expanding the potential itself we have a generalized tadpole-like contribution in which the number of closed lines is constrained to be even because of trivial topological reasons. In dimensional regularization the linear term contributes to the renormalization of the potential only if r = 1 and d = 2 as we show later in subsection 2.3. At the quadratic order we have instead In (2.4) and (2.5) we introduced G(x, x ′ ) which is the Green function associated to the operator of the quadratic part of the action A covariant representation of the Green function for an operator of Laplace-type as the one above is described in Appendix A. For our present needs, the representation simply shows that the Green function can be expanded in which we purposely neglected all further contributions which do not affect the relevant operators. The leading G 0 (x, x ′ ) term can be understood as a covariant generalization of the standard Green function of flat space (see Appendix A for more details), while a 1 (x, x ′ ) is the first correction due to curvatures. In the following subsection we consider first the renormalization of the general φ 2n universality class for n ≥ 3, while the case n = 2 is deferred for later. The reason for this is that the case n = 2 is special when it comes to the renormalization of the function F(φ). In particular, the results for the general φ 2n case often cannot be continued to n = 2 because the subleading correction to the propagator is powerlaw for each d = d n with n ≥ 3, but it is logarithmic in d = d n=2 = 4. If the analytic continuation is performed anyway, there is thus an additional "unbalanced" singularity which is seen as an additional 1/(n − 2) pole in the beta functions.

φ 2n universality class
The leading quadratic contribution to the path integral (2.5) is not a one loop contribution for all n ≥ 3 models, but rather it involves (r − 1)-loops, which is a marked distinction from the more familiar analyses of φ 4 and Yang-Mills theory below the upper critical dimension d = 4. However, not all loop contributions to (2.5) lead to 1/ǫ poles. Using the methods described in Appendix B and dimensional analysis, it is possible to infer that 1/ǫ poles arise for the cases r = n and r = 2n − 1, corresponding to (n − 1)and (2n − 2)-loop diagram respectively. In the case r = n, the contribution arises solely from r lines of the leading G 0 (x, x ′ ) term of the Green function. In the second case the diagram can be either composed by 2n − 1 lines of G 0 (x, x ′ ), or by 2n − 2 lines of G 0 (x, x ′ ) and one of G 1 (x, x ′ ). In practice, this makes for three multiloop diagrams that must be evaluated by the methods described in Appendix B. The diagrams are depicted in Fig. 1 We have that in d = d n − ǫ the three diagrams evaluate to in which we integrated by parts one derivative to cast the kinetic-like term of the second diagram in a suitable form, and defined the constant The results of (2.8) are essentially the counterterms which must be inserted to remove the divergences of all the relevant operators of the φ 2n model in curved space for n ≥ 2. The pole at n = 2 of the last counterterm is a clear indication of why we left the φ 4 models out of this general discussion. Since we are just considering a leading renormalization, the computation of the renormalization group flow is straightforward because it can be obtained by simply acting on the counterterms with the logarithmic derivative with respect to the reference scale µ ∂ ∂µ . Naturally, we display the RG in the guise of functional equations. We also include a field dependent wavefunction Z(φ) as renormalization of the kinetic term. The wavefunction is generated by the flow and, while it includes irrelevant contributions for the most part, the use of a boundary condition for Z(0) allows for the determination of the anomalous dimension of the renormalized field. At the upper critical dimension we find (2.10) In a rather standard fashion we switch to the dimensionless renormalized canonicallynormalized field Diagrammatic representation of (2.8) in order of appearance. The first and second diagrams are made of n and 2n − 1 lines of the leading contribution of the Green function G 0 (x, x ′ ). Trivially their symmetry factors are n! and (2n − 1)! respectively. The third diagram is again made of 2n − 1 lines, but one corresponds to the subleading G 1 (x, x ′ ) line which is depicted as dashed. Its symmetry factor is (2n − 2)! because there are (2n − 1) ways to choose the last line. which includes a rescaling by the wavefunction renormalization constant Z 0 = Z(0) which is generated by β Z . The field ϕ is the natural argument for the dimensionless Their renormalization group flow is By construction we have that z(0) = 1, so the limit ϕ → 0 of its flow can be used to determine the anomalous dimension η ≡ −∂ log Z 0 /∂ log µ directly from β z | ϕ=0 = 0.

φ 4 universality class
The four dimensional case is special for three main reasons. Firstly, diagrams and counterterms leading to the renormalization are not directly obtained as the analytic continuations to n = 2 of the results of Sect. 2.1. Secondly, the subleading correction to the Green function in four dimensions is logarithmic. This means that in the ǫ-expansion an additional divergence must be subtracted from the propagator as we show in (A.14). The difference in the behavior of the subleading part of the propagator is the reason why a 1 n−2 pole appears in the third diagram of (2.8). Thirdly, a simple dimensional analysis reveals that operators quadratic in the curvatures have the same canonical dimension of the operators φ 4 and φ 2 R, and hence must be renormalized together for consistency.
Here we try to follow the notation of [19] for the most part with some minor modification. Let us first generalize the action (2.1) to accommodate the higher curvatures with the following invariants (2.14) These invariants are chosen so that in four dimensions G integrates to a topological invariant and F , which is the square of the Weyl tensor, transforms covariantly under scale transformations. It is convenient to define one general function and its modification as follows At one loop, which is the leading order, the counterterm to U(φ, R) can be obtained by a simple application of the heat kernel. One finds that the leading contribution to the renormalization of U(φ, R) comes from the a 2 (x, x) coefficient given in (A.10) while the wavefunction renormalization is a two loop effect completely analog to the limit n = 2 of Sect. 2.1. The computation of the leading beta function is straightforward Returning to the original functions of (2.13) we find the functional beta functions as well as the beta functions for the higher derivative couplings Since F(φ) is at most quadratic we have that F ′′ (φ) = F ′′ (0) and the couplings c and e can be treated as numbers, even though the right hand side suggests otherwise. In order to make this section on φ 4 more self consistent, we briefly discuss some critical property of the above system. This discussion anticipates some points that are made later in the development of Sect. 3. One can see that at the leading order the critical value for the nonminimal coupling ξ = F ′′ (0) is ξ = 1 6 as one would naively expect from continuing the general conformal value (1.3) to d = 4, thus proving the educated guess of the introduction for the special case n = 2.
In general it is not guaranteed that the critical value of ξ remains a fixed point beyond the leading order unless a further renormalization condition is exploited [18]. For our purpose, it would be interesting to know if the general conjecture of the introduction is true, that is, we would like to know under which circumstances at two loops and for d = 4 − ǫ the coupling takes the value One can prove, using naively the dimensionally regulated scheme at the next-to-leading order (NLO) and a straightforward subtraction, that the above value is not a fixed point to order ǫ. However, the freedom highlighted in [19] of redefining the potential U(φ, R) by a copy of the one loop counterterms can be exploited to ensure that (2.20) is the fixed point at NLO. The redefinition is a change of the renormalization conditions which thus defines and links the metric and the field. We refer to [19] for a more complete and detailed explanation of the results reported in this section.

φ ∞ universality class: the Sine-Gordon model
The upper critical dimension d = 2 emerges as the limit d n → 2 of n → ∞. The renormalization of the path integral for the two dimensional case is very simple, even though it represents a special case likewise the φ 4 one. It is convenient to borrow the notation from the previous section and use the full potential U(φ, R). The computation of the leading counterterms and beta functions necessitates only the use of the standard heat kernel expansion of an operator of Laplace-type, and specifically of the coefficient a 1 (x, x) given in (A.10). We find the leading counterterm at one loop and deduce the very simple RG beta functional Notice that there is no anomalous dimension renormalization coming from our leading order computation.
In two dimensions the scale invariant solutions of this beta function become periodic. It has been argued that the critical solution of this RG flow in flat space is periodic and corresponds to the Sine-Gordon universality class [3]. Here we are observing a generalization to curved spacetime for zero anomalous dimension as in [26]. Let us first introduce the dimensionless potential u(ϕ, R) = µ −2 U(ϕ, µ 2 R). Using the boundary conditions u(ϕ, R) = u(−ϕ, R) and ∂ 2 φ U(φ, R)| φ=0 = m 2 , at the fixed point in d = 2 we find Notice that we have imposed the boundary conditions as a function of the scalar curvature, therefore an implicit dependence on R might in principle be hidden in the mass m 2 = m 2 (R). In this way we have ensured that the result agrees both with the assumption that this solution generalizes the Sine-Gordon universality to curved space, and with the expectation that the nonminimal coupling ξ should be zero at the critical point.

2d gravity at large-c
As a brief intermezzo we believe that it is interesting to show the relevance of the results of Sect. 2.3 in reproducing some well-known result of two dimensional quantum gravity coupled to conformal matter. Let us recall that in exactly two dimensions the path integral of gravity can be determined by integrating the conformal anomaly [27], which leads to a renormalization procedure linked to a nonlocal action known as the Polyakov action [21]. This action is especially relevant because the spacetime integral of the Einstein term is a topological invariant in two dimension, and hence it cannot govern the dynamic of the model.
However for general d (and specifically for d = 2 − ǫ) the Einstein term is not a topological invariant, and therefore it has been argued by Kawai and Ninomiya that it should be possible to reproduce the results based on the Polyakov action by just renormalizing the Einstein action in d = 2 − ǫ and then taking the limit ǫ → 0 [28]. The validity of this argument was shown through the course of several papers, which ultimately lead to the two loop renormalization of the Einstein action in d = 2 − ǫ. For more details we refer to [29] and references therein; notice however that in the literature of 2d gravity it is often chosen d = 2 + ǫ, therefore the replacement ǫ → −ǫ is necessary when comparing results.
The renormalization of dimensionally regulated two dimensional gravity is slightly unconventional because it has to deal with the conformal factor of the metric, otherwise one finds discontinuities when analytically continuing to ǫ → 0 [30]. In order to describe it, let us first introduce the Einstein action interacting with c distinct conformally coupled fields φ i in d dimensions We require that the coupling ξ c is determined by the conformal value (1.3) and assume that this condition can be preserved through renormalization (see the discussion of sections 1 and 3 for more details on this point). The number c is often referred to as "central charge" and it counts the effective number of matter degrees of freedom.
In two dimensions all possible metrics are related by a Weyl transformation, and therefore only their conformal mode is allowed to fluctuate. Close to two dimensions, instead, it is customary to parametrize the metric g µν → (ǫ/8) 2/ǫ ψ 4/ǫ g µν into a conformal mode ψ and a metric g µν which is not allowed to fluctuate in its trace part (by abuse of notation we denote the transformed metric with g µν ). 2 Using this normalization the mode ψ of the metric enjoys a Weyl invariant action, which is in form analog to any of those of the fields φ i , if not for an overall negative sign which makes ψ an unstable "scalar" degree of freedom. The idea of [32] is to transform (2.24) into and renormalize it such that the function L(ψ, φ i ) respects the conformal coupling. The new function is normalized by L(0, 0) = 1, which is a necessary condition to read off the value of the Newton constant G.
Assuming that the instability of ψ can be cured by opportunely Wick rotating the theory, it is possible to neglect the effects of the dilaton field ψ as compared to those of the multiplet φ i ; moreover one can argue that for large values of c the fluctuations induced by the fields φ i dominate over those of g µν too. In other words, for large-c it should be necessary to integrate only the loops of φ i , but this is exactly the multifield generalization of what we have done in section 2.3 upon the identification for the dimensionless versions of L and G. Notice that in the large-c limit we are dropping any parametric dependence on the mode ψ to highlight the connection with the previous section. Now we use (2.26) inside (2.22) to determine the renormalization group flow of the renormalized G and L(φ i ). In order to separate the two beta functions we have to impose that L(0) = 1 along the flow. Additionally we impose that all fields φ i are coupled in the same way so that it will be sufficient to denote each one of them by φ. We find (2.27) The interaction with the fluctuating modes of g µν can change the anomalous dimension of the fields φ i as η ∼ G, but this contribution is also generally subleading in the limit of large central charge. We follow the strategy of [32] and parametrize L = 1 + aψ + bψ 2 − ξ c φ 2 . It is straightforward to see that the beta functions of a and b have Gaussian solutions, thus setting all couplings except for G at the respective fixed points we obtain . This result agrees with the large-c limit of the exact leading result in which the constant A takes the value A = 25−c 24π [28]. Notice that the general Euclidean result hinges on our ability of solving the problem of the instability of the conformal mode, which in [32] is "Wick" rotated ψ → iψ. While several solutions have been proposed there is no definite answer, nor general consensus on how to approach the problem. In fact, proposals to solve the problem without a Wick rotation of the dilaton mode have received renewed attention recently [33]. This problem can be framed in the more general discussion of finding the universality class of quantum gravity and exploring the corresponding conformal theory [34]. Here we would like to mention another less explored yet interesting possibility that was outlined in [35]: the path integral of 2d gravity could be "defined" starting with the path integral of a fluid 2d membrane embedded in D bulk dimensions (which is essentially a noncritical Nambu-Goto string) and analytically continuing to D → 0. In the membrane path integral the correct counting of the degrees of freedom involves the propagation of modes of the extrinsic curvature, which play a role analog to the gauge fixing ghosts.

φ 2n+1 universality class
The results of Sect. 2.1 can in part be generalized to the tower of multicritical nonunitary models φ 2n+1 [13]. While it can be seen that the quadratic leading term in the renormalization of the potential is replaced by a more involved cubic one, it turns out that the leading renormalization of the wavefunction and the non-minimal coupling function is contained in the same diagrams with the opportune change of the number of internal lines.
A simple rule of thumb to test the validity of the functional RG equations involves the replacement n → n + 1 2 in the system (2.12): the beta function β v ceases to make sense signaling that it should be replaced with a cubic term, but both β z and β f are still meaningful and in fact they are the correct beta functions. If the leading cubic flow of the potential, which is given in [13], is included, then it is trivial to generalize the system (2.12) to the entire sequence of φ 2n+1 models.
We checked explicitly that in the special case φ 3 the resulting system, which corresponds to the Lee-Yang universality class with a nonminimal coupling to the curvature, coincides with the one given in the appendix of [26]. Notice that the cubic case should be treated with more care than we do, because the complete renormalization requires counterterms for all the cubic invariants coming from the metric [23], and therefore it is as special as the quartic case described in Sect. 2.2. The general next-to-leading renormalization (occurring at two loops) of φ 3 in curved space appeared for the first time in [23]. The renormalization of the model is actually known with position-dependent couplings [25], a result which is used to show that a natural six dimensional generalization of the Zamolodchikov's c-function does not always increase monotonically with the flow [24].

CRITICALITY
We now resume the analysis of the RG system (2.12) of Sect. 2.1 representing the general case of the multicritical model φ 2n for n ≥ 3. The φ 4 model is an outlier, so we anticipated a brief discussion of the critical properties of its nonminimal coupling to the curvature already in Sect. 2.2. We find convenient to rescale the potential while leaving all other functions intact. The system (2.12) simplifies to Using the boundary condition z(0) = 1 in the rescaled flow β z , we also determine the anomalous dimension of the scalar field η = 4v (2n) (0) 2 /(2n)!. The critical couplings appear as the leading couplings of the potentials v(ϕ) and f (ϕ). By construction, in the minimal subtraction scheme all other couplings are dimensionful, and therefore are zero at the critical point. We therefore parametrize the potentials in terms of the two almost marginal interactions Using the above parametrization in (3.2), we find the following beta functions and anomalous dimension It is clear that η contributes to the cubic order in λ of β λ , which has no effect to the determination of the order ǫ of the fixed point. However η has an important effect in β ξ because its contribution scales with the same power of λ as the other terms. Substituting η we find The system has two different fixed points. On the one hand we have the Gaussian fixed point at λ = 0 which sets both beta functions to zero. In this case the natural fixed point for ξ is the subleading root of β ξ . On the other hand we have the non Gaussian fixed point . (3.6) For both fixed points the coupling ξ takes the critical value that is expected at the upper critical dimension which evidently proves the educated guess given in the introduction. The above analysis can be extended easily to the multicritical nonunitary models φ 2n+1 following the guidelines explained in Sect. 2.5. Interestingly, the only outlier of our analysis is the case for n = 2, for which we have to use the set of beta functions coming from (2.18) as discussed in Sect. 2.2. However, it is straightforward to find that in this case ξ = 1 6 which happens to coincide with the continuation of (3.7) to n = 2. The next step would be to test if the next-to-leading order correction to the non-Gaussian fixed point of ξ matches the ǫ-expansion of conformal value for the coupling ξ evaluated in d = d n − ǫ instead of d = d n in agreement with the general conjecture of the introduction. This would imply which comes from the expansion of (1.3) to orders of ǫ using d = d n − ǫ. Following the discussion of [18,19], which we reproduced briefly in Sect. 2.2, we argue that ensuring (3.8) probably requires a special choice in the renormalization conditions leading the the RG flow. In practice, the next-to-leading contributions to the RG flow can be changed by the inclusion of terms which match the counterterms (2.8) and which can be used to change the renormalization conditions leading to the fixed point value for ξ. We hope to come back to this topic. Let us include here also a short remark on the steps that have lead to (3.7). While our educated guess of the introduction stated that we expected that ξ takes the conformal value at criticality, the validity of the guess is not at all obvious from the initial form of the counterterms (2.8). In particular, there is a very delicate balance among the terms appearing in the renormalization (2.12) and the anomalous dimension which produces the form of β ξ in (3.5) and which makes evident that the conformal value (3.7) is actually the critical point. We conclude this section by discussing the implications that the system of beta functions (3.5) has on the infrared physics. For obvious reasons, we are mostly interested in studying the renormalization group flow in a physical dimension. The first natural dimension (smaller than d n ) in which almost all models for n ≥ 3 are nontrivial is d = 2, we therefore continue ǫ to the value ǫ = 2 n−1 to continue the φ 2n models to the physical dimension d = 2. Correspondingly, the fixed point value of the coupling λ becomes λ * = 2 A n in which we define A n = (2n)!/(n!) 2 which is simply the coefficient of the λ 2 term in β λ . The flow can be integrated as follows in which we introduce B n = 8(n 2 − 1)/((n − 2) (2n)!) that is the coefficient of β ξ . The flow satisfies the ultraviolet boundary conditions λ 0 = λ(µ 0 ) and ξ 0 = ξ(µ 0 ), which can be checked by setting µ = µ 0 in (3.9). More interestingly, we can use (3.9) to explore the infrared limit µ = 0. One can see trivially that the second term in the denominator of λ(µ) drops for µ = 0 and therefore we have λ(0) = λ * . Slightly less trivial is to show that for µ → 0 the integral appearing in the exponential of ξ(µ) diverges logarithmically implying that the second term drops; we thus have ξ(0) = ξ n . These results are in line with the expectation that the nontrivial fixed point (3.6)-(3.7) is of infrared nature in that it controls the large scale behavior of the model near criticality.

CONCLUSIONS
We have considered the leading order renormalization of the multicritical scalar models with φ 2n interaction in curved space. Our analysis shows that for almost all values of n one has to consider counterterms for the self interaction as well as for the nonminimal interaction of the form φ 2 R, while some additional counterterms based on curvature invariants are needed in the special case n = 2. The counterterms have been obtained from a computation of the 1 ǫ poles of dimensionally regulated covariant Feynman diagrams of (n − 1)-loops for the self interaction, and (2n − 1)-loops for the self energy and the nonminimal interaction.
Our result generalizes the renormalization of the φ 4 model in curved space, which we have considered as a special case, but it also shows that the general case functions rather differently. Specifically, the structure of the counterterms for the nonminimal coupling displays a discontinuity for n = 2, which corresponds to φ 4 . We have deduced a set of functional beta functions which describes the scale dependence of a self-interaction potential and a generalized nonminimal interaction with the scalar curvature.
We have used the perturbative renormalization group flow to determine standard perturbative beta functions for the two canonically marginal couplings: λ of the selfinteraction φ 2n and ξ of the nonminimal interaction φ 2 R. The RG system clearly shows that at the leading order the scale invariant fixed point of the nonminimal coupling ξ coincides with its dimension-dependent conformal value ξ c evaluated at the upper critical dimension. This result is in agreement with an educated guess enunciated in the introduction. Importantly, the leading critical value for the coupling ξ is an ultraviolet attractive feature of the renormalization group flow.
We have also discussed a more general conjecture for which at the next-to-leading order the ǫ-expansion of the fixed point value of ξ matches the expansion of the conformal value ξ c below the upper critical dimension. Based on similar and already available results for the φ 4 [18] and φ 3 models [23][24][25], we argue that one has to either follow a modified version of the prescription of Brown and Collins [18], or alternatively to subtract normally while exploiting the freedom of redefining the renormalization group flow at the next-to-leading order using the counterterms of the leading order [2]. In other words, one might want to find the appropriate renormalization condition which ensures the validity of the conjecture for the nonminimal coupling. We believe that this condition plays an important role, especially if it is necessary to describe the model in a conformal or Weyl-invariant way.
A clearer understanding of the status of the conjecture and the renormalization condition can only be achieved by building on the results of this paper and studying the next-to-leading order contributions to the renormalization group flow, which is thus an important future prospect for this computation.
in which we introduced the biscalar δ-function that generalizes the usual flat space Dirac delta. The propagator that is used in the main text can be obtained by specifying the endomorphism as E = F ′′ (φ)R. It is convenient to represent the Green function using the heat kernel method. The heat kernel function is defined as the solution of the following differential equation If we solve the diffusion equation implicitly then the relation of the heat kernel function with the Green function is straightforward For all intents and purposes the above relation should be taken as our operative definition of G(x, x ′ ). The heat kernel representation is useful because the solution admits an asymptotic expansion for small values of the parameter s, known as the Seeley-de Witt expansion, which captures the ultraviolet properties of the Green function. The expansion is generally parametrized as We introduced several bitensors in the expansion. The most fundamental is σ(x, x ′ ), sometimes known as Synge's or Synge-de Witt's world function, which is half of the square of the geodesic distance between the points x and x ′ [36]. The bitensor ∆(x, x ′ ) is known as the van Vleck determinant and is related to the world function and the determinant of the metric as Together, the bitensors σ(x, x ′ ) and ∆(x, x ′ ) ensure that the leading term of the Seeleyde Witt parametrization covariantly generalizes the solution of the heat equation in flat space with O ∼ −∂ 2 . Finally, the bitensors a k (x, x ′ ) are the coefficients of the asymptotic expansion and contain the geometrical information of the operator O, which includes curvatures, connections and interactions. It is well-known that ultraviolet properties are (and must be) local in renormalizable theories. For the case of the heat kernel and the Green function locality corresponds to x ∼ x ′ and it is captured by the so-called coincidence limit in which x → x ′ . Given any bitensor B(x, x ′ ), its coincidence limit is defined Notice that covariant derivatives do not generally commute with the coincidence limit ∇[B] = [∇B], but rather satisfy a modified relation [36,37]. The coincidence limits of the bitensors σ(x, x ′ ) and ∆(x, x ′ ) and their derivatives can be obtained by repeated differentiation of the crucial relations for which we suppressed bitensor coordinates and we used the notation in which subscripts of σ(x, x ′ ) correspond to covariant derivatives. Similarly, coincidence limits of the coefficients a k (x, x ′ ) can be obtained by differentiating and inductively using with the boundary condition σ µ ∇ µ a 0 = 0. In the relevant example of a simple scalar field the first coefficient is trivial a 0 (x, x ′ ) = 1, because the Seeley-de Witt expansion solves the diffusion equation in flat space. We give here the first two nontrivial coincidence limits for the expansion of the operator (A.1) which are used in the computations of the main text Using (A.6) in (A.5), we obtain an analog expansion of the Green function The leading G 0 (x, x ′ ) and the subleading G k (x, x ′ ) for k ≥ 1 are bilocal contributions to the Green function and are determined by a simple integration over the heat kernel parameter s Depending on the theory and its dimensionality, the exponent d/2 − 1 − k inevitably becomes negative for a certain value of k highlighting the fact that there is only a finite number of Green's function contributions which are singular in the limit x ∼ x ′ . It is very important to point out that in (A.12) we have implicitly assumed that the Green contributions do not scale logarithmically with σ(x, x ′ ) at the critical dimension. While they do not scale logarithmically for almost all the multicritical models considered in this paper (that is, for all φ 2n with 3 ≤ n < ∞). They do, however, scale logarithmically if there are values of k for which d = 2 + k, which corresponds to poles of the gamma function. In this case, one sees the failure of capturing the logarithmic behavior explicitly through the ǫ → 0 limit of (A.12) which is not regular even outside x ∼ x ′ .
For this reason, when d = 2 + k one needs to subtract an ǫ-pole to (A.12) for the results to be valid at and close to the dimension d. For example, close to d = 2 the leading propagator is already logarithmic and we subtract As desired, the above expression is valid for d close to two dimensions and is regular in the limit ǫ → 0 for d = 2 − ǫ. We report here the subleading part of the Green function for d = 4 − ǫ, which is also needed in the paper Notice that if both sides of (B.2) are multiplied by the same bilocal operator, then the Dirac delta on the right hand side allows for the substitution of its coincidence limit Q(x, x ′ )δ (d) (x, x ′ ) = [Q]δ (d) (x, x ′ ). The core of the algorithm is thus to transform all possible inverse powers of the world function into those of the left-hand-side of the basic relation, substitute them with the right-hand-side, and then sort all bilocal operators at the numerator so that they enter in contact with the Delta function.
Higher inverse powers of the world function can be manipulated inverting which can be proven easily using (A.8). For the purpose of this paper we just need which is obtained inverting (B.3) for b = d/2 and using the coincidence limits of the biscalars [∆ 1/2 ] = 1 and [Y] = R/6. Generalizations of (B.4) including higher inverse powers can also be easily obtained, however further iterations of (B.3) typically exhibit bilocal operators which are separated from the Dirac delta by the presence of covariant derivatives (imagine, for example, placing Q(x, x ′ ) on both sides of (B.4)) and hence their coincidence limit cannot be taken. In such cases, it is necessary to integrate by parts all covariant derivatives so that all bilocal operators come in contact with the Dirac delta. For example, if one covariant derivative is located between the bilocal operator and the Delta we manipulate as follows In the second line we have exploited the Delta to take the coincidence limit of the neighboring operators. A similar manipulation can be performed for the case of two derivatives and results in In order to obtain further generalizations, one has to integrate by parts all derivatives one-by-one, and take the coincidence limits only of operators which are in direct contact with the Dirac delta. Generalizations of (B.5) are thus straightforward but rather lenghty. Systematic applications of (B.3), to manipulate the inverse powers of the world function, and of (B.5), to take the local parts of the biscalars multiplying the divergences,