Two-point correlators in non-conformal N = 2 gauge theories

We study the two-point correlation functions of chiral/anti-chiral operators in N = 2 supersymmetric Yang-Mills theories on R4 with gauge group SU(N) and Nf massless hypermultiplets in the fundamental representation. We compute them in perturbation theory, using dimensional regularization up to two loops, and show that field-theory observables built out of dimensionless ratios of two-point renormalized correlators on R4 are in perfect agreement with the same quantities computed using localization on the four-sphere, even in the non-conformal case Nf 6= 2N .


Introduction
Supersymmetric Yang-Mills theories (SYM) have long been considered as an ideal playground to get exact results in Quantum Field Theory. In recent times, exact formulae for special observables in theories with extended supersymmetries have been found. In fourdimensional theories with maximal N = 4 supersymmetry, the exact resummation of the infinite series of perturbative corrections to the expectation value of circular Wilson loops, also in presence of chiral operators, has been performed [1][2][3][4][5]. These results are based on the counting of the relevant rainbow-like Feynman diagrams by means of a matrix model. The introduction of this matrix model has been considered ad hoc until it was shown [6] that localization for the N = 2 * theory on the four-sphere S 4 , after having performed the N = 4 limit, predicts its existence. Moreover, localization provides a non-perturbative formula for the circular Wilson loop in a general N = 2 theory which takes into account both perturbative and non-perturbative, instanton and anti-instanton, corrections in an interacting matrix model. A two-loop test of this formula against perturbation theory was presented in [7] in the case of superconformal QCD.
It is natural to ask whether localization on S 4 can be used to compute non-trivial quantities other than the Wilson loop expectation value. In [8][9][10][11][12][13][14] it has been proposed that two-point correlators between chiral and anti-chiral operators in a superconformal N = 2 theory on R 4 can be computed from the partition function of the theory on the four-sphere with chiral and anti-chiral insertions at the north and south pole respectively; localization expresses this partition function as a matrix model. In a conformal N = 2 theory, two-point correlators between a chiral operator O n = tr(ϕ n 1 )tr(ϕ n 2 ) · · · , where n = (n 1 , n 2 , . . . ) and ϕ is the complex scalar of the gauge vector multiplet, and an antichiral operator O m made out of the complex conjugate field ϕ, take the form where n = i n i and m = j m j are the scaling dimensions of the two operators. G n, m (g 0 ) is a non-trivial function of the coupling constant g 0 , but bears no dependence on the distance x since chiral and anti-chiral operators are protected in conformal N = 2 theory. The two-point functions on a four-sphere also take the form (1.1) but with x 2 being the chordal distance on S 4 . The function G n, m (g 0 ) is the same on the sphere and in flat space, and it is given by a two-point function in a matrix model obtained from localization. Explicit tests of the match between the field theory and the matrix model descriptions of the correlator (1.1) have been performed up to two loops for low-dimensional operators in SU(2) and SU (3) gauge theories with N f = 4 and N f = 6 matter hypermultiplets. The results where extended in [15] to generic chiral operators in a superconformal SU(N ) theory with N f = 2N . Also the one-point functions of chiral operators in presence of a circular Wilson loop can be expressed in terms of the matrix model obtained via localization of the Wilson loop on S 4 , as checked up to two loops in [16].
It is of obvious importance to investigate to what extent the matrix model description of the correlation function persists in N = 2 theory in non-conformal set-ups. A first step -2 -

JHEP05(2019)199
in this direction was carried out in [15] for a SU(N ) theory with N f flavors, where suitable operators were chosen in such a way that their two-point correlators vanish in perturbation theory up to a given loop order, leaving a finite contribution at the next leading loop order. In such a situation, a perfect match between the perturbation theory and localization was shown for an arbitrary rank and any number of flavors at two and three loops. This strongly hints that chiral/anti-chiral correlators are related to the S 4 matrix model also beyond the conformal case.
Considering generic chiral/anti-chiral correlation functions away from the conformal point N f = 2N , we encounter important differences. For N f = 2N , the gauge coupling and the operators are not anymore protected from quantum corrections and have to be renormalized to account for the ultraviolet (UV) divergences. As a consequence, the two-point correlation functions of renormalized operators O R n and O R m depend on the renormalization scale µ and are no longer forced to have just a power-like dependence on the distance x as in (1.1). Instead, they take the general form where g = g(µ) is the renormalized coupling and the dimensionless quantity parametrizes the distance separation. Here γ E denotes the Euler-Mascheroni constant. In contrast to the conformal case, the function G R n, m (g, ν) depends non-trivially on the distance x 2 through the quantity ν. Such a dependence cannot be obtained from localization on S 4 , because in this case the operators are inserted at the opposite poles of the four-sphere and the distance between them is fixed in terms of the sphere radius. Still our results show that, up to two loops, the ν-dependence is very simple and can be put in the factorized form G R n, m (g, ν) = G R n, m (g, 0) 1 + 1 2 β 0 g 2 ν n δ nm + O(g 6 ) , (1.4) where β 0 is the expansion coefficient of the exact one-loop β-function of the theory, namely (1.5) Remarkably, the ν-dependent prefactor in (1.4) depends only on the scaling dimension n and not on the details of the operators. As a consequence, up to two loops at least, ratios of correlators of the same scaling dimension are actually ν-independent and can be compared directly against the matrix model results. We show that the field theory results for such observables are indeed in perfect agreement with the predictions from the localization matrix model. This is consistent with the fact that all Feynman diagrams contributing to the ratios are finite in four dimensions, and finite loop integrals on S 4 and R 4 yield the same result after replacing propagators in flat space by those on the four-sphere, as we will show in section 5.

JHEP05(2019)199
Furthermore, one can ask whether the results for the renormalized correlators at a given renormalization scale can be directly matched against those coming from localization. At one loop, we show that this is indeed the case if we choose µ 2 x 2 = e γ E /π in the minimal subtraction scheme. Moreover, by considering the field theory on S D and evaluating the relevant one-loop integrals in dimensional regularization, we find that, apart from the obvious replacement of propagators, not only the divergent parts but also the finite parts agree with the results on R D for D → 4. While the agreement of the divergent part is expected, since the divergences are sensible to short distances and do not distinguish between the sphere and flat space, the agreement of the finite part is neither expected nor guaranteed a priori, but nevertheless it holds. At two loops, we find that the matrix model results reproduce the field theory ones, up to a term proportional to (2N − N f ) and to the dimension n of the operators. This suggests that the difference could be interpreted as a conformal anomaly which, in non-conformal theories, affects the correlation functions in going from the four-sphere to the flat space.
In this paper we keep the numbers of colors N and flavors N f arbitrary and compute the two-point correlators for a general choice of chiral/anti-chiral operators. On the field theory side we do this at the two-loop level. We summarize our findings in section 2 and give a detailed account of the Feynman diagram computations in appendices A and B. To keep track of the various particles exchanged in the loops we use a superfield formalism. The loop integrals are evaluated using the integration methods pioneered in [17] (see also [18] for a review). In section 3 we compute the renormalized correlators and their anomalous dimensions. Our results suggest that the anomalous dimensions γ n,0 of the chiral operators are one-loop exact and are given by the simple formula γ n,0 = n 2 β 0 . In section 4 we discuss the computation of the correlators on the matrix model side, building on the techniques described in [15], and compare the results with those previously obtained from the field theory side. To facilitate the comparison, we show that, up to two loops, the localization matrix model can be re-expressed as a complex matrix model encoding the color factors and the combinatorics of the Feynman diagrams that contribute to the chiral/anti-chiral correlators. Finally, in section 5 we discuss the field theory calculation of the two-point correlators on the four-sphere S 4 , and in section 6 we present our conclusions. Several technical details for such calculations are provided in appendix C.

Two-point correlators from perturbation theory
We consider a N = 2 SYM theory with gauge group SU(N ) and N f hypermultiplets in the fundamental representation. For N f = 2N the theory is conformally invariant also at the quantum level. We denote by ϕ(x) the complex scalar field of the N = 2 vector multiplet which, in N = 1 notation, is the lowest component of a chiral superfield Φ. In this theory a basis of chiral operators can be given in terms of the multi-trace operators O n (x) = tr (ϕ n 1 (x)) tr (ϕ n 2 (x)) . . . tr (ϕ n (x)) ,

JHEP05(2019)199
where n = (n 1 , n 2 , . . . , n ). The scaling (bare) dimension of O n (x) is We expand the scalar field ϕ(x) = ϕ a (x) T a over the SU(N ) generators T a (a = 1, . . . , N 2 − 1) in the fundamental representation, normalized as In terms of the components ϕ a (x), the operators (2.1) become where R a 1 ...an n is a completely symmetric tensor. 1 In an analogous way, we define the anti-chiral operators O n (x) using the complex conjugate field ϕ(x) instead of ϕ(x).
We are interested in computing the two-point correlation functions in non-conformal N = 2 theories using standard perturbative techniques. We perform our calculations at the origin of moduli space, where the scalar fields have vanishing vacuum expectation values. This is a preferred point in the sense that here the breaking of conformal invariance occurs only at the quantum level, as a consequence of the dimensional transmutation phenomenon. Therefore, this is the natural place in which to test whether the matrix model approach based on localization agrees with the standard perturbative field-theory calculations also in the non-conformal case. Unlike the N = 4 theory where the correlators (2.5) are exact at tree-level, in N = 2 theories they receive quantum corrections, starting from one loop for N f = 2N and from two loops in the conformal case N f = 2N . Moreover, in the non-conformal theories, the loop integrals are UV divergent, in general, and need to be regularized. Here, we use the dimensional regularization taking the space-time dimension to be D = 4 − 2 . As a consequence, the bare gauge coupling constant, g 0 , becomes dimensionful. In general, the bare two-point functions (2.5) take the form where n = m is the common scaling dimension of the two operators, and 1 Explicitly, R a 1 ...an n = tr T (a 1 · · · T an 1 tr T a n 1 +1 · · · T a n 1 +n 2 . . . tr T a n 1 +...+n −1 +1 · · · T an) where the indices are symmetrized with strength 1.  is the massless scalar propagator in D-dimensions. The correlator G n, m (g 0 , , x) can be computed at weak coupling as an expansion in powers of g 2 0 . We refer to [15] for details on the Feynman rules that are needed to perform this calculation; they are summarized for convenience in appendix B.
The diagrams which contribute to the two-point functions (2.5) up to order g 4 0 are schematically represented in figure 1. The diagram (a) is the tree-level contribution, the diagram (b) is the one-loop correction, while the other four diagrams represent the twoloop part. The blobs labeled by v k, stand for the sum of all irreducible diagrams of order g 2 0 with k external lines -half of them connected to the chiral fields ϕ of O n , half to the anti-chiral ones of O m .
A convenient way to organize the computation of these diagrams is to consider the N = 4 theory, remove all contributions from Feynman diagrams involving loops of the adjoint hypermultiplet (which we call H) and add those with loops of the fundamental matter multiplets (which we call Q and Q) [7]. Since the two-point correlators in N = 4 theory are exact at tree-level, we can write where, in an obvious notation, G n, m H stands for all diagrams in the N = 4 theory with the adjoint hypermultiplet H circulating in the loops, and G n, m Q, Q stands for the same diagrams in the N = 2 theory with loops of fundamental matter multiplets Q and Q. In the following, we will sometimes refer to this method as "performing the computation in the difference theory". We stress that in the difference theory one should take into account only diagrams involving loops of the adjoint hypermultiplet or loops of the fundamental ones, but not both.

Tree-level
The tree-level contribution to the correlator (2.6) comes from the diagram in figure 1(a). To obtain its explicit expression, one contracts the fields ϕ in O n with the fields ϕ in O m by means of a free scalar propagator (2.11) Explicit expressions can be easily found also for operators with higher dimension.

One-loop diagrams
We first observe that in the N = 4 theory there are no one-loop corrections to the propagators of the adjoint scalars Φ I (I = 1, 2, 3). At one loop, this is schematically represented in figure 2. This implies that there is no one-loop correction to the propagator of the adjoint hypermultiplet H in the difference theory. The one-loop correction to the propagator of the fundamental matter superfields Q and Q vanishes as well. Indeed, as shown in figure 3, the contribution of this diagram is similar to the previous one upon replacing the generators in the adjoint with those in the fundamental representation, so that the same cancellation mechanism at work for the adjoint scalars applies here as well.   To find the contribution of the diagram shown in figure 1(b), we need to compute the one-loop correction to the propagator of the scalar field in the N = 2 theory. This is represented in figure 4, where in the second line we have used the N = 4 result of figure 2 for I = 1, to replace the diagram with the vector propagator by the diagram with a scalar loop. Explicitly computing these superdiagrams (see appendix B for details), we find The correction (2.12) can be of course inserted in any of the n propagators connecting O n and O m , so that the one-loop contribution to the two-point correlator (2.6) corresponding to figure 1 (2.14)

Two-loop diagrams
At order g 4 0 there are several diagrams that contribute to the correlator (2.6). They are schematically represented by the last four diagrams, from (c) to (f ), of figure 1. A detailed derivation of the various contributions and the evaluation of the corresponding loop integrals can be found in appendix B. Here we simply summarize the results for the building blocks of each of these diagrams. In the difference theory, the diagram (c) contains as building blocks the diagrams represented in the left-hand side of figure 5. If we exploit the identity of figure 2 to replace the diagrams with the vector field propagator in favor of the ones with a scalar loop, we can easily realize that the diagrams in the left-hand side figure 5 precisely reconstruct the square represented in the right-hand side. Using (2.12), and taking into account the appropriate multiplicity factor of the diagrams, this gives Let us now consider the contribution corresponding to the two-loop diagram in figure 1(d).
In this case, the insertion of two one-loop corrections on the same scalar propagator leads to the diagrams displayed in the first two lines of figure 6. If we use again the identity of figure 2 to replace the diagrams containing the vector field propagator with those with a scalar loop, we reconstruct the square of the one-loop correction, as shown in the last line of figure 6. Evaluating explicitly the loop integrals in this case, we obtain that the result can be written as the square of the one-loop up to terms of order (see appendix B for b v 2,1 v 2,1 Figure 6. In the difference theory, the reducible diagrams that correct the ϕ propagator at two loops can be expressed in terms of the one-loop contribution. Figure 7. A class of diagrams that correct the scalar propagator at two loops. The dashed double-line notation in the right-hand side is a convenient way to represent this contribution in the difference theory. details), namely Here and in the following, we use the approximate symbol ≈ for equations that hold up to terms vanishing in the limit → 0.

v 2,2 -and v 4,2 -contributions
Let us now consider the two-loop irreducible corrections to the scalar propagator which appear in figure 1(e). A first class of contributions arises from correcting one of the internal lines in the one-loop diagrams of figure 4 by using the one-loop propagator of the matter superfields Q and Q, or of the adjoint hypermultiplet H. However, as we have seen before, these contributions vanish.
Other terms that correct the scalar propagator at two loops in the difference theory are those represented in figure 7. Here we have introduced the dashed double-line notation as a convenient way to represent the difference between the loop with fundamental flavors and the loop with the adjoint hypermultiplet. Actually there are other three classes of irreducible diagrams that correct the scalar propagator at two loops. In figure 8 we have drawn all such diagrams, whose evaluation is presented in appendix B to which we refer for details. Summing all contributions, we find that the irreducible two-loop correction to -10 -   In terms of the constant matrix ϕ, this result can be obtained from the insertion of the quadratic vertex (2.18) in the matrix model; indeed In the difference theory, at two loops there are irreducible connected diagrams involving two chiral and two anti-chiral fields, which enter the term (f) in Fig. 1. These diagrams can be further distinguished according to their overall colour structure and their contribution, denoted as v 4,2 in Fig. 1, can be split into two terms which we denote as v The first class of such diagrams is represented in Fig. 10.
In the difference theory, there are irreducible two-loop contributions that involve two chiral and two anti-chiral fields and give rise to the diagram of figure 1(f ) . These contributions can be further distinguished according to their overall color structure and can be split into two terms which we denote as v  figure 9, where again we have used the dashed double-line notation to represent the difference between Q and H loops. They are explicitly with f abc being the SU(N ) structure constants. In (2.19) the ellipses stand for terms with color tensors that are anti-symmetric in (a 1 , a 2 ) and (b 1 , b 2 ). Such terms do not contribute to the two-point correlation functions (2.5) because they are contracted with the symmetric tensors R n and R m defined in (2.4). The last two-loop diagram we have to consider is the one represented in figure 10. This diagram was already computed in [15] and the result, which is also reviewed in appendix B, is Again, the ellipses in (2.22) stand for terms with color tensors which are anti-symmetric in (a 1 , a 2 ) and (b 1 , b 2 ) and, therefore, vanish when inserted in the two-point correlation function.

Effective vertices
In order to later compare the results of perturbation theory to those of the matrix model, we find it convenient to introduce effective vertices following the ideas of [15]. In particular, to obtain the color dependence of the two-point function G n, m we strip the x-dependence of the scalar fields and introduce the adjoint matrix ϕ = ϕ a T a and its conjugate ϕ = ϕ a T a , such that We denote by O n (ϕ) the operator obtained by replacing in (2.1) the field ϕ(x) with the constant matrix ϕ, and by O m (ϕ) the same with ϕ(x) replaced with ϕ.
With these definitions, it is straightforward to see that the tree-level correlator (2.10) can be written as Also the one-loop correlator (2.14) can be written in a simple way using this formalism.
As usual, the notation : : stands for normal ordering. Using (2.25), it is easy to check that inside a vacuum expectation value (2.27) we can use the relation that follows from the SU(N ) fusion/fission rules valid for two arbitrary N × N matrices B 1 and B 2 . Let us now compute the two-loop contribution to the two-point correlation function from diagrams shown in figure 1. Using the reducible term (2.15) we find that the contribution of the diagram in figure 1 In terms of the effective vertex (2.28), this result can be rewritten as follows: where in the second line we have used the identity that follows from Wick's theorem. The two-loop reducible correction (2.16) to the scalar propagator can be inserted in any of the n internal lines. Thus, the contribution of the diagram in figure 1 where we used (2.29). In a similar way, the irreducible two-loop correction (2.17) produces the following contribution to the diagram in figure 1(e): 4,2 given, respectively, in (2.19) and (2.22). To write the results in a compact form, it is convenient to introduce the quartic vertices Then, the contribution of the effective vertex (2.19) to the correlator can be written as Repeatedly using the fission/fusion identities (2.30), it is possible to show that, inside the vacuum expectation value (2.38), the following relation holds: More generally, one can prove that By comparing with (2.29) we conclude that the quartic vertex V can be effectively replaced by (−V 2 ) inside a vacuum expectation value, so that (2.38) becomes The contribution of the effective vertex (2.22) to the correlator can be treated in an analogous way, and it reads (2.42) -14 -

JHEP05(2019)199
Notice that, in distinction to the other two-loop contributions, the expectation value (2.42) is not proportional, in general, to the tree-level correlator G n, m tree due to the structure , and it has to be computed case by case. A few explicit examples with operators of even dimensions are: The contribution to the two-point correlator from the two-loop diagram of figure 1(f ) is given by the sum of (2.41) and (2.42), namely where v

Summary of results
Collecting our findings, up to two loops the bare correlator is given by implies that the total two-loop irreducible contribution in the second line of (2.45) is finite for → 0. The only divergences remaining at two loops come from the square of those -15 -

JHEP05(2019)199
present at one-loop. At the given loop order, they can be nicely combined into an overall factor (see the first line of (2.45)) (2.49) Then, defining the effective interaction vertex we can recast (2.49) in a very compact way as follows: . (2.51) Indeed, expanding the exponentials up to order g 4 0 , we precisely recover all terms of (2.49). In particular, from the insertion of a single V eff in the correlator we obtain the linear terms in the v k 's appearing in the first and third line of (2.49), while the quadratic terms in the second line arise from two insertions of V eff . Notice that since this effective vertex is normal-ordered, the denominator in (2.51) contributes up to order g 4 0 only with the term proportional to [V 2 (ϕ, ϕ)] 2 appearing in the second line of (2.49).
In section 4 we show that localization on a four sphere produces an expression similar to (2.51). However, in order to compare the two expressions, we should first get rid of the UV divergences and scheme ambiguities that are present in the bare correlator G n, m . This is the content of the next section.

Renormalization
The dimensionally regularized bare correlators G n, m given in (2.45) are divergent for → 0, since the one-loop coefficient v 2,1 defined in (2.13) behaves for small as

JHEP05(2019)199
As we have remarked before, the UV divergence due to v 2,1 is the only one that plagues the correlators, since all other terms in G n, m are finite for → 0. To get rid of this divergence, we have to apply the standard renormalization procedure. First, we introduce the dimensionless renormalized gauge coupling constant g through the relation where µ is an arbitrary scale, and Z is a suitable function to be determined. Then, we define the renormalized operators O R n (x) according to where Z m n is a matrix-valued function. However, in the previous section we have shown that the divergences of the two-point functions depend only on the scaling dimensions of the operators and not on the operator details; therefore a block-diagonal matrix can do the job, and we can simplify (3.3) by setting A similar formula holds for the anti-chiral renormalized operators O R n (x). The singular terms of the functions Z(g 2 , ) and Z n (g 2 , ) are determined by requir- should be finite when expressed in terms of the renormalized coupling g. This means that the renormalized correlator is well-defined and free of divergences in the limit → 0.

The β-function and anomalous dimensions
The dependence of the renormalized coupling g 2 and of the renormalization constant Z n (g 2 , ) on the energy scale µ is described, respectively, by the β-function and by the anomalous dimensions γ n (g 2 ) of the operators O n . They are defined as follows: where the last equality stems from the µ-independence of g 0 , and where in the second step we used (3.6). Using the perturbative expansions β(g 2 ) = −2 g 2 + β 0 g 4 + β 1 g 6 + . . . , γ n (g 2 ) = γ n,0 g 2 + γ n,1 g 4 + . . . , we can explicitly integrate (3.6) and (3.7) and get in the minimal subtraction (MS) scheme This shows that the expansion coefficients of the β-function and of the anomalous dimensions γ n are directly related to the coefficients of the 1/ -terms in Z and Z n respectively. Differentiating (3.5) and using dG n, m /dµ = 0, one obtains the Callan-Symanzik equation on which we can safely take the limit → 0. We now determine the coefficients of the β-function and the anomalous dimensions, using the explicit results of the previous section. To do so, we expand the bare correlator (2.45) for small using (3.1), and write where α 0 = G n, m tree and Plugging the expansion (3.11) into the renormalized correlator (3.5), using (3.9) and requiring that all divergent terms cancel, one finds leading to This value of β 0 is in agreement with the well-known result for the one-loop coefficient of the β-function in N = 2 SQCD. We also notice that γ 2,0 = β 0 . This is consistent with N = 2 supersymmetry, since the chiral operator O (2) = tr ϕ 2 and the Yang-Mills Lagrangian −tr (F 2 /4) + . . . belong to the same supermultiplet, and thus should renormalize in the same way, that is Moreover, using the fact that in N = 2 SYM theories the β-function receives only one-loop correction, 2 i.e. β = 0 for all ≥ 1, we conclude that also the anomalous dimensions of JHEP05(2019)199 tr ϕ 2 are corrected only at one loop, i.e. γ 2, = 0 for all ≥ 1. This implies that Furthermore, we deduce from (3.14) that the following relation holds up to two loops. It would be very interesting to investigate whether this relation holds also at higher loops. While this issue is not relevant for the two-loop analysis of the present paper, it is tempting to speculate that (3.17) might actually be true in general. Indeed, in our set-up the anomalous dimensions of the chiral and anti-chiral operators arise because of the breaking of conformal invariance at the quantum level due to dimensional transmutation. The fact that the coefficients γ n,0 and β 0 are proportional to each other and that the proportionality factor is n/2 (see (3.14)), together with N = 2 supersymmetry, naturally leads one to propose the relation (3.17). Notice that in the conformal case N f = 2N , the renormalization functions simply reduce to 1, due to the absence of divergences, so that (3.17) is trivially satisfied in this case.

Renormalized correlators
Using the previous results, it is easy to see that up to two loops and in the limit → 0, the renormalized correlators (3.5) take a form completely analogous to the bare correlators (2.45), namely G R n, m = 1 + n c 1 + n(n + 1) 2 and The coefficients c i are obtained from (3.1), (2.46) and (2.47). In the conformal case c 1 vanishes and the first perturbative correction to the correlator appears at order g 4 . Performing the same manipulations as described in section 2 for the bare correlators, we can rewrite (3.18) in the following form: are defined, respectively, in (2.28) and (2.37).
We remark that for N f = 2N , the renormalized correlator G R n, m is not a constant, but it depends on x through the ln πµ 2 x 2 term contained in ν. At first sight, this fact makes it unlikely that the correlator can be encoded in a matrix model. However, the dependence of G R n, m on ν is determined by the Callan-Symanzik equation (3.10) 3 and by its value at a reference point ν. For instance, setting we get from (3.20) which, as we will see in section 4, is the combination that matches the one-loop matrix model results from localization. In particular, with the choice (3.23) the coefficient c 1 becomes Using (2.11) and (2.43) we find from (3.18) Similar explicit formulae can be worked out for correlators involving higher dimensional operators. For example, at dimension 4 we find where the coefficients c 1 , c 2 and c 3 are defined in (3.19). Notice that c 1 and c 2 enter into these expressions in the same combination. 3 This equation requires that the correlator must have the form It is easy to check that (3.18) satisfies this requirement. Notice that the whole ν-dependence of G R n, m can be reconstructed order by order in g 2 from the correlator at a given value ν, for instance ν = 0, and the coefficients β0, γn,0 and γn,1.

Normalized correlators
The renormalized correlators G R n, m are finite but they cannot be considered as physical observables since they depend on the choice of the renormalization scheme. In particular, given a renormalized correlator at a certain normalization scale µ, one can always perform a finite renormalization of the operators by multiplying the renormalization factors Z and Z n defined in (3.2) and (3.4), respectively, by an arbitrary finite function of the coupling. This transformation preserves the UV finiteness of the correlator and corresponds to a change of the renormalization scheme (see, for example, [25] for a discussion of this point in a different context).
Since, up to two loops we have Z n = (Z 2 ) n 2 , we can eliminate the scheme dependence by considering dimensionless ratios of correlators. In fact, the renormalized correlators can be written as (3.28) This relation shows that the dependence on x 2 and on the renormalization scale µ coming from the coefficient c 1 is entirely encoded in a prefactor which only depends on the bare scaling dimension n of the operators but not on their specific form. Therefore, this prefactor cancels in the ratio of correlators of operators of the same dimension. Choosing, for example, as a reference the correlator between two operators O (2) , which are the only ones with dimension 2, we are led to introduce the normalized correlators (3.29) These ratios are independent of the choice of the renormalization scale µ and scheme, and as such they represent physical quantities. It is interesting to observe that the two-loop contribution to the ratio A R n, m only comes from the irreducible diagram represented in figure 10, which is finite in the limit → 0 (see (2.22) and (2.23)). This shows that also the bare ratios are finite. The equality between A n, m and A R n, m comes from the fact that the Z n -factors cancel between the numerator and denominators, and the Z-renormalization of the gauge coupling starts to contribute at the next order. One can check this explicitly, by writing the bare correlators (2.45) as 31) which shows that the divergence encoded in the one-loop coefficient v 2,1 cancels in the ratios A n, m . Furthermore, by comparing (3.31) and (3.28), we can realize that the bare and the renormalized ratios match at two loops, apart from the obvious replacement of g 0 with g.

JHEP05(2019)199
The explicit expressions of the normalized correlators for operators of dimension 4 are: Similar formulae can be easily found also for operators of other dimensions.
In the next section we will recover this same result from the matrix model obtained by applying localization on the sphere S 4 .

Matrix model approach
In [6] it was shown, using localization techniques, that the partition function of a N = 2 SYM theory with gauge group SU(N ) defined on a four-sphere S 4 can be written in terms of a traceless Hermitian N × N matrix a in the following way: where a uv = a u − a v , and for simplicity have set to 1 the radius R of the four-sphere. 4 Furthermore, Z(ia, τ ) is the gauge theory partition function with τ being the complexified gauge coupling: In this paper we actually set the θ-angle to zero. We remark that in the non-conformal theories the coupling g appearing in the matrix model has to be interpreted as the renormalized gauge coupling at a scale proportional to the inverse radius of the four-sphere [6]. The gauge theory partition function Z(ia, τ ) is computed using the localization techniques of [26,27] with a purely imaginary vacuum expectation value ϕ = i a for the adjoint scalar, and an Ω-background with parameters 1 = 2 = 1/R. This partition function can be written as a product of the classical, one-loop and instanton contributions, namely (4.4) 4 The dependence on R can be easily restored by replacing au with au R.

JHEP05(2019)199
Since we work at weak coupling g 2 1, where instantons are exponentially suppressed, we can set Z inst (ia) = 1. The classical part produces a simple Gaussian term in the matrix model: (4.5) The one-loop contribution arising from the gauge multiplet and N f matter multiplets can be written as [6] (see also [15] for details) where S n (a) are homogeneous polynomials in a of order n. The first few of them are: Performing the rescaling the matrix model gets a canonically normalized Gaussian factor and the sphere partition function (4.1) becomes a v (4.9) with S int (a) = g 2 8π 2 S 2 (a) + g 2 8π 2 2 S 4 (a) + g 2 8π 2 3 S 6 (a) + · · · . (4.10) The term of order g 2 in S int (a) accounts for effects that take place at loops in the corresponding field theory computation. Therefore we will refer to the g 2 -expansion of S int as a loop expansion. Exploiting the Vandermonde determinant ∆(a) and writing a = a b T b , we can alternatively express the integral (4.9) using a flat integration measure da over all matrix components a b as follows where c N is a g-independent constant and da ∝ b da b . The overall prefactors in (4.11) are irrelevant when computing correlators and thus can be neglected.

JHEP05(2019)199
Given any function f (a), its vacuum expectation value in the matrix model described above is defined as follows where in the second line we have used (4.11). The subscript "0" denotes the vacuum expectation values taken with respect to the Gaussian measure, which can be computed by repeatedly using Wick's theorem to reduce them to the basic contraction a b a c 0 = δ bc . (4.13)

Chiral and anti-chiral operators in the matrix model
We are interested in extracting from the matrix model (4.12) the two-point functions (2.5).
To this aim we have first to find counter-partners of the chiral and anti-chiral operators in the matrix model. It would seem natural to associate to the multi-trace operator O n (x) defined in (2.1), an analogous function O n (a) in the matrix model, given by the same expression (2.1) but with ϕ(x) replaced by a, namely O n (a) = tr a n 1 tr a n 2 . . . tr a n . (4.14) However, the operator O n (x) has vanishing vacuum expectation value in the field theory, while in the matrix model O n (a) = 0 due to the self-contractions of a. This means that we have to refine the dictionary and make O n (a) normal-ordered. This can be done by subtracting from O n (a) all possible self-contractions and making it orthogonal to all operators with lower dimensions. As discussed in [15,16,28], the prescription to define the normal ordering of any operator O(a) in the matrix model is the following. Let be ∆ the dimension of O(a) and O p (a) a basis in the finite-dimensional space of matrix operators with dimension smaller than ∆. Denoting by C ∆ the (finite-dimensional) matrix of correlators which are computed according to (4.12), we define the normal-ordered operator as (4.16) Our notation stresses the fact that this normal-ordering is g-dependent, since the correlators on the right hand side of (4.16) are computed in the interacting matrix model. The proposal is then to associate to the field theory operators the corresponding normal-ordered matrix operators, namely A similar replacement holds for the anti-chiral operators.

JHEP05(2019)199
For example, using the definition (4.16) we find (4.18) The term of order g 0 inside the square brackets represents the self-contraction of tr a 2 , while the terms of higher order in g 2 represent the self-contractions of the operator through the interaction vertices coming from the matrix model action. Analogous expressions can be worked out for operators of higher dimension.

Correlators in the matrix model
Once the operators have been identified, their correlators can be computed in a straightforward way using the definition (4.12). In particular the two-point correlators are defined as Since normal-ordered operators with different dimensions are orthogonal to each other, G n, m vanishes for n = m. For instance, for the simplest operator O (2) (a) defined in (4.18) we find (4.20) The explicit expressions of correlators with higher dimensional operators can be computed in a similar way. At dimension 4 we find Here we have split the O(g 4 ) contribution into the sum of a few terms in order to facilitate a comparison with the field theory calculation. It is easy to check that for N f = 2N the matrix model correlators G n, m exactly match, up to two loops, the correlators G R n, m computed in perturbation theory (see (3.26) and (3.27)), thus confirming the general results obtained in [8][9][10][11][12][13][14][15]. 5 The fact that the partition function on the sphere S 4 and its associated matrix model contain information on the correlators in the flat space R 4 is not too surprising in the conformal case. We now want to investigate to what extent this relation holds in the non-conformal theories with N f = 2N .

Comparison between matrix model and field theory correlators
Comparing (3.26) with (4.20), and (3.27) with (4.21), we see that they have the same structure and that many terms exactly match. However, for N f = 2N there are some differences in the terms proportional to ζ(3). To make the comparison simpler, it is convenient to write G n, m in terms of the complex matrices ϕ and ϕ using the formalism introduced in section 2. Indeed, it is possible to explicitly check that, up to two loops, the matrix model correlators (4.19) can be expressed as follows defined, respectively, in (2.28) and (2.37), and

(4.24)
Notice that the effective vertex (4.23) has the same form as the renormalized vertex (3.22) obtained from perturbation theory. Comparing (4.24) with (3.19) and (3.25), we find (4.25) 5 In the recent paper [29] a discrepancy at six loops, proportional to ζ 2 (5)g 12 , has been pointed out in the comparison between the matrix model results and the correlators obtained by solving Toda equations.

JHEP05(2019)199
Therefore, the difference between the effective vertex V eff of the matrix model and the renormalized effective vertex V R eff is It is interesting to observe that δ vanishes in the conformal case. Moreover, it is proportional to V 2 which, as follows from (2.29), computes the scaling dimension of the operators. 6 This fact suggests that it might be interpreted as due to a conformal anomaly which, in nonconformal theories, affects the correlation functions in going from the four-sphere S 4 to the flat space R 4 , or vice versa. In the two-loop approximation, we can rewrite (4.22) as follows: (4.27) This formula clearly shows that the dependence on c 1 and c 2 drops out in the ratio between correlators of operators with the same scaling dimensions. Thus, in analogy with (3.29), we are led to define the ratio of correlators in the matrix model (4.28) Since c 3 = c 3 , it exactly matches the normalized correlator A R n, m , namely We have checked this relation in many explicit examples, with operators of dimensions up to 6.

Two-point correlators on the four-sphere
In this section we study in more detail the relation between the correlators in flat space, discussed in sections 2 and 3, and those on the four-sphere S 4 . The latter are closely related to the correlators derived from matrix model presented in section 4. In particular we consider the one-loop correction to the scalar propagator on S 4 and compare it with the one-loop propagator in flat space defined in (2.12).
To this aim, it is convenient to describe a sphere in D−dimensions by using flat embedding coordinates {η 0 , η µ } satisfying the quadratic constraint

JHEP05(2019)199
where R is the radius of the D-sphere. Following [30][31][32][33][34], we use the stereographic projection to relate a theory defined on a D-sphere to a theory in R D , parametrized by the flat coordinates x µ . One of the advantages of this formalism is that the scalar propagator on the sphere, denoted by a subscript S, takes a very simple form given by Here we have used D = 4 − 2 and defined Inserting this into (5.4) and comparing with (2.7), we get The scalar propagator on the sphere is thus proportional to the one in flat space, with a scaling factor raised to the engineering dimensions of the scalar fields. Notice that this is the same scaling factor that defines the induced metric on the sphere through the conformal map (5.2); indeed Let us now consider the correlators between two operators on the sphere. They have a structure similar to the ones in flat space given in (2.6), namely n, m (g 0 , , η 12 ) δ nm . Inserting this into (5.8), using the propagator (5.4) and taking the limit → 0, we can easily obtain This is the expected relation between correlators on the sphere and correlators in flat space that follows from the conformal map (5.2).

JHEP05(2019)199
Let us now consider the one-loop correction. Before analyzing the correlators on the sphere, it is convenient to revisit the calculation of one-loop correction to the scalar propagator in flat space, given in (2.12) in coordinate space. The one-loop correction to ϕ a (x 1 ) ϕ b (x 2 ) can be written as Its Fourier transform is the function W 1 (p) discussed in appendix B (see in particular (B.8) and (B.9)). Computing the integrals, we find Using this in (5.11), one recovers the result presented in (2.12) and (2.13).
Going to the sphere, we find that the one-loop correction to the scalar propagator has a form similar to (5.11), that is where the function W 1 S is the sphere generalization of W 1 . Applying the embedding formalism [30][31][32][33][34], the expression of W 1 S can be obtained by performing the conformal transformation (5.2) to (5.12). Under this map, both the integration measure and the scalar propagators acquire scale factors according to where (5.17) Comparing this integral with (5.12), we notice the presence of the additional scaling factor κ(x 3 ) κ(x 4 ) − , which clearly becomes 1 in four dimensions.
Therefore, if the integrals in (5.12) and (5.16) were finite, W 1 S and W 1 would only differ by the overall scaling factor κ(x 1 ) κ(x 2 ). In other words, if no UV divergences are present, one can safely perform the limit → 0 inside the integrals. However, the integral -29 -JHEP05(2019)199 in (5.17) is divergent, and thus the scaling factor κ(x 3 ) κ(x 4 ) − in the integrand cannot be neglected. The evaluation of this integral is presented in appendix C, and the result is Comparing with (5.13), we see that, up to terms O( ), the two expressions coincide upon replacing ∆ S with ∆, and η 2 12 with x 2 12 . The fact that the divergent parts of W 1 S and W 1 coincide, is not surprising since the UV divergences come from integration at short distances where there is no distinction between the sphere and flat space. What is non trivial, however, is that the finite parts coincide, modulo the obvious replacement of x 12 with η 12 .
Putting everything together, we see that the one-loop correction to the scalar propagator on the sphere is 20) in full analogy with (2.12) and (2.13). This implies that Thus, the renormalization procedure can be done following the same steps we described in section 3. Choosing the renormalization scale µ 2 as in (3.23) with x 2 replaced by η 2 12 on the sphere and by x 2 12 in flat space, then The relation (5.21) and the explicit expression of v (S) 2,1 explain why the correlators G n, m obtained from the matrix model perfectly agree with those computed in field theory at one loop.
The same analysis can be carried out at two loops, even though the resulting integrals on the sphere become way more complicated. Most of the two-loop diagrams develop UV divergences and need to be regularized. As in the one-loop case, the integrals on the sphere differ from those in flat space because of scaling factors [κ(x i )] − appearing in the integrands. Such factors do not modify the leading UV divergent contribution but they do affect the finite part. As a result, there is no reason a priori to expect that the finite part of the correlation functions on the sphere and in flat space should coincide. However, we stress the fact that the finite part of the two-point correlator is not a physical observable since it depends on the regularization scheme. The explicit one-loop calculation shows a perfect agreement between the matrix model and the field theory results for the two-point correlators for the special choice of the renormalization scale. It is natural to ask whether such identification holds also at higher loops. At two loops, the results of sections 3 and 4 reveal that the finite part of the correlation functions are different in flat space and in the matrix model. Still, a perfect match is found for physical observables that are independent of the JHEP05(2019)199 renormalization scheme, such as the ratios of correlators with operators of the same dimension. In such ratios, all divergent two-loop diagrams cancel and the whole contribution at order g 4 is due to a single and finite Feynman diagram, namely the irreducible diagram represented in figure 10. The corresponding Feynman integral is finite in R 4 and does not require a regularization. As a consequence, it possesses the four-dimensional conformal symmetry and takes the same form in R 4 and S 4 . This explains why the ratios of the correlation functions match the prediction from localization at two loops, as shown in (4.29).

Summary and conclusions
We have explicitly computed the two-point correlation functions between chiral and antichiral operators in the N = 2 SYM theory with gauge group SU(N ) and N f fundamental flavors up to two loops, using standard (super) Feynman diagrams in dimensional regularization. Our results show that these correlators have a remarkably simple structure of UV divergences stemming from the fact that the anomalous dimensions of the operators are proportional to the β-function. We demonstrated that when the renormalization scale µ and the separation x between the operators are inversely proportional to each other, these correlators can be obtained via a matrix model which is strikingly similar to the matrix model that computes the partition function and the chiral/anti-chiral correlators on the four-sphere using localization. Up to two loops, the difference between the two matrix models is just a term of order g 4 proportional to (2N − n f ) V 2 , which acting on the operators gives their anomalous dimensions. This suggests that this difference that vanishes in the conformal theories, might be interpreted as a conformal anomaly. In the non-conformal cases this could explain the difference between the correlators on the four-sphere and those in flat space.
We have also constructed normalized correlators, which are scheme independent and, as such, represent physical quantities. Up to two-loops, these normalized correlators are the same on the four-sphere and in flat space, and can be computed either using the field theory approach with Feynman diagrams, or using localization methods via a simple matrix model.
Our analysis clarifies the relation between the perturbative field calculations and the localization results in N = 2 SYM theories. It would be interesting to generalize it in various directions, for example to compute the correlators at three or more loops, or to compute the one-point or higher-point correlation functions in presence of Wilson loops. In particular it would be interesting to explore in detail the two-loop calculations of the correlators using Feynman diagrams on the sphere, and/or obtain a "first-principle" derivation of the difference between the two matrix models that yield the correlators on the four-sphere and in flat space. We hope to be able to return to some of these points in future works.

JHEP05(2019)199
The work of M.B., A.L., F.F., J.F.M. is partially supported by the MIUR PRIN Contract 2015MP2CX4 "Non-perturbative Aspects Of Gauge Theories And Strings". The work of G.P.K. is supported by the French National Agency for Research grant ANR-17-CE31-0001-01. The work of A.L. is partially supported by the "Fondi Ricerca Locale dell'Università del Piemonte Orientale". All authors would like thank the Galileo Galilei Institute for Theoretical Physics for hospitality during the course of this work. G.P.K. is grateful to INFN and the Simons Foundation for partial support.

A Loop integrals
In this appendix we follow closely [17] (see also [18] for a review) and collect some useful formulae necessary to evaluate the Feynman integrals. We work in D = 4 − 2 dimensions and use the propagator of a massless scalar field given in (2.7), namely For later convenience, we introduce the graphical notation for Feynman integrals in the momentum representation α The black dots on the left and the right of the diagram in (A.2) denote, respectively, the incoming and outgoing momentum p. Furthermore, in each interaction vertex the momentum conservation is enforced. When α or β is 1, for simplicity we do not write the labels. With these notations, we then have .
We will also make use of the Fourier transform integral which for α = 1 reduces to (A.1). In particular we will need the following explicit formulae Figure 11. The one-loop diagram corresponding to the integral (A.7). Here the labels α i on the various lines denote the exponents of the propagators appearing in the integrand.

A.1 Triangle identity
Let us consider the integral which corresponds to the triangle diagram of figure 11. Following [17,18], we have From this, it is easy to obtain the so-called triangle identity: where the notation n ± J {α i } means the integral (A.7) with α n replaced by α n ± 1. For example, we have (A.10) Repeated applications of the triangle identity allow us to reduce the power of one of the propagators to zero and to express in the end the result in terms of the basic integrals (A.2). A few examples are described in the next subsection.

A.2 Scalar integrals
Let us consider the two-loop integral

JHEP05(2019)199
Here we have adopted the same graphical conventions as in (A.4). Applying the triangle identity (A.9), we obtain The last steps in these equations follow from (A.2). Inserting these expressions into (A.12) and expanding for → 0, we obtain (A.14) After Fourier transforming and using (A.6b), we get The same procedure can be applied to express other two-loop integrals in terms of I α,β defined in (A.3). For example, we have (A. 16) Computing the Fourier transform and using (A.6c), we find Another scalar integral that will be needed is the one represented by the diagram Using (A.4) and expanding for small , one can prove that

JHEP05(2019)199
where the ellipses stand for terms that vanish for → 0. Comparing with (A.16), we easily conclude that so that, after Fourier transform, we have In a similar way one can derive the following relation from which we get Performing the Fourier transform we obtain The following divergent integrals also appear in the two-loop calculation

B Evaluation of the relevant (super)diagrams
In this appendix we explicitly compute the diagrams discussed in section 2. We use dimensional regularization and the N = 1 superspace formalism in the Feynman gauge (we refer to [15] for more details).

B.1 Feynman rules
We first summarize the momentum-space Feynman rules in the chosen formalism. Let us start from the propagators for the chiral multiplets. We use a continuous line for the superfields Φ I (I = 1 Note that the arrow indicates both the orientation of the chiral propagator and the flow of the momentum. In (B.1a) we have used the notation

JHEP05(2019)199
Our conventions on spinor indices and Pauli matrices are the same as those explained in appendix A of [15].
The propagator for the N = 1 vector superfield is given by The diagrams we have to compute only contain three-point vertices. These are given by the following rules: Here (T a ) bc = −i f abc are the generators in the adjoint representation, and (T a ) uv those in the fundamental representation. The θ variables appearing in the vertices are those associated to the vertex point.

B.2 One-loop diagrams
At one loop, we have to compute the diagrams in figure 4. Denoting and isolating a prefactor containing the combinatorial and color factors, we have where W 1 (x) is given by the one-loop Feynman diagram shown in (B.8) below. The first term in the color factor arises from the loop diagram of N f fundamental Q,Q superfields, while the second term originates from the loop of the adjoint hypermultiplet H. It is convenient to Fourier transform W 1 (x) and write where W 1 (p) is described by the following diagram in the momentum space: Here the numbers label the external points and the interaction vertices. Note that the two external points 1 and 2 are connected to a bosonic scalar field, so that the propagators from -38 -JHEP05(2019)199 1 to 3 and from 4 to 2 are in fact free scalar propagators with no θ dependence. Taking into account the Feynman rules given above, we get In the second step we used the Grassmann integral identity and then we exploited the graphical representation introduced in appendix A. Using (A.25a) and the Fourier transform (A.26a), we finally obtain leading to the relation (2.13) in the main text:

B.3 Two-loop diagrams
At two loops we have to compute diagrams that correct either a two-point or a four-point vertex. Such diagrams have been displayed in section 2. All of them have two external points, corresponding to the positions of the two operators O n (x) and O m (0), four internal points, corresponding to the interaction vertices, and either seven or eight propagators, for the corrections to the two-point or the four-point vertex respectively.
Some notations. It is useful to introduce some notation that allows us to write the various diagrams in a uniform way. In each two-loop diagram labeled by an index I, we label by i = 1, 2 the external points and by i = 3, . . . , 6 the internal ones. We denote by E I the set of propagators of the diagram, and label each propagator in E I by s. Any propagator connects a point i to a point j, and in general there can be a number r(i, j) of propagators connecting the same two points. A possible way of expressing the label s of the propagators is thus s → (i, j; r) (B.13) where r = 1, . . . , r(i, j); in the following we will omit the index r if r(i, j) = 1. The momentum k s associated to the propagator will then be denoted as k ij;r , with the convention that we take it to flow from i to j. This is useful to write the delta-functions of momentum conservation at internal vertices, which take the form (B.14)

JHEP05(2019)199
Similarly, the relation between the internal momenta and the external momentum p is enforced by Just as in the one-loop case, any two-loop diagram will be written as the product of a factor containing the weights in the vertices, the combinatorial and color factors, and of a colorless diagram W I (x), which we will obtain from its Fourier transform W I (p). The latter has the following structure The factor Y I contains the contribution of the propagators and the conditions for the momentum conservation at each vertex of the I-th diagram, while the factor D I (p, k, θ,θ) contains all θ andθ terms coming from the vertices and from the superfield propagators. Some of the Grassmann integrations yielding might be obvious, in which case we will indicate only the non-trivial integrals and denote the integrand of Z I as D I .
Reducible diagrams. The reducible two-loop diagrams are represented in figure 5 and figure 6. In the diagrams of figure 5 there are two independent one-loop corrections to a propagator line. Hence the result follows simply from the one-loop computation presented in the previous subsection, and has been given in (2.15) of the main text. Let us consider then the diagrams of figure 6. The overall factors are simply the square of those of (B.6), and thus we get The where D 2 = exp 2θ 4 (k 43;1 + k 43;2 )θ 3 + 2θ 6 (k 65;1 + k 65;2 )θ 5 + 2θ 4 k 45θ5 = exp −2θ 4 pθ 3 − 2θ 6 pθ 5 + 2θ 4 pθ 5 . (B.21) In the second step we used the momentum conservation δ-functions that are present in the factor Y defined in (B.16). Performing the Grassmann integrals, one finds Z 2 (p, k) = p 4 . This factor cancels the two "external" propagators in Y 2 and the integral over internal momenta can be represented in the graphical notation of appendix A as follows: Taking the Fourier transform of this expression via (A.26b) and inserting it into (B.18), we finally find in agreement with the formula (2.16) in the main text.
Irreducible diagrams: the v 2,2 part. Let us now consider the irreducible two-loop corrections to the scalar propagator, namely the diagrams represented in figure 8. We start from In momentum space, we have to compute where in the second step we used momentum conservation. Performing the θ-integrals, we get Z 3 (p, k) = p 4 , which cancels the two "external" propagators in Y 3 ; the remaining integral over internal momenta can be represented in the graphical notation of appendix A as follows: which has the general form (B.16). Again, the θ-factors of the vertices saturate the integrations over θ 3 andθ 4 , while the gluon propagator sets θ 6 andθ 6 equal to θ 5 andθ 5 . The remaining Grassmann integrations are as in (B.26), but now with D 4 = exp −θ 5 (k 53 + k 45 + k 46 + k 62 )θ 5 + 2θ 5 k 53θ3 + 2θ 4 (k 45 + k 46 )θ 5 + 2θ 4 k 43θ3 = exp −2θ 5 pθ 5 − 2θ 5 k 53θ53 − 2θ 4 k 43θ53 (B.35) where in the second step we used momentum conservation. The Grassmann integrations yield Z 4 (p, k) = p 2 k 2 43 . This factor cancels one external and one internal propagator in Y 4 and we remain with The Grassmann integration leads to Z 5 (p, k) = p 2 k 2 53 , which cancels one external and one internal propagator in Y 5 . We then remain with Integrating overθ 4 produces a factor of (k 34 + k 56 ) 2 θ 2 34 which sets θ 4 = θ 3 , so that in the end we remain with Z 6 (p, k) = (k 34 + k 56 ) 2 d 2 θ 3 d 2θ 3 exp −2θ 3 pθ 3 = (k 34 + k 56 ) 2 p 2 ∼ 2 (k 34 · k 56 ) p 2 .
(B.47) The last step follows from the fact that when we integrate this expression over momenta, both the k 2 34 and the k 2 56 terms, canceling the corresponding propagator, give rise to tadpolelike integrals, which vanish in dimensional regularization. The symmetry of the diagram under k 56 ↔ −k 65 allows us to rewrite the above result as This diagram was already computed in [15] in configuration space. For completeness we report here its evaluation in momentum space. Using the relation As noted after (B.58), the last two terms in the square brackets are anti-symmetric in (a 1 , a 2 ) and (b 1 , b 2 ). Therefore they vanish when inserted in a chiral/anti-chiral two-point function and can be discarded. The momentum space diagram corresponding to W 10 (x) is − 2 (k 43 · k 63 ) (k 45 · k 65 ) + 2 (k 43 · k 65 ) (k 45 · k 63 ) − 2 (k 43 · k 45 ) (k 63 · k 65 ) .
Here we have neglected a term proportional to the to the anti-symmetric ε-tensor coming from the trace of four Pauli matrices given in (B.63) which does not contribute for symmetry reasons. Recalling that p = k 13 + k 15 = k 42 + k 62 , we can exploit momentum conservation and discard terms involving p 2 and p · k 43 , which give contributions that vanish for → 0 after Fourier transform. After some algebra we are left with in agreement with formula (2.23) of the main text and with the findings of [15]. This completes the calculation of the two-loop diagrams contributing to the chiral/antichiral correlators.

C Feynman integral on the sphere
In this appendix, we evaluate the integral I(x 1 , x 2 ) that appears in the expression for the one-loop correlation function on the sphere, given in (5.14) and (5.16).