Schouten identities for Feynman graph amplitudes; The Master Integrals for the two-loop massive sunrise graph

Transcription

1 Available online at ScienceDirect Nuclear Physics B Schouten identities for Feynman graph amplitudes; The Master Integrals for the two-loop massive sunrise graph Ettore Remiddi a, Lorenzo Tancredi b, a Dipartimento di Fisica, Università di Bologna and INFN, Sezione di Bologna, I-06 Bologna, Italy b Institut für Theoretische Physik, Universität Zürich, Wintherturerstrasse 90, CH-8057 Zürich, Switzerland Received 0 November 03; accepted 3 January 0 Available online 6 January 0 Abstract A new class of identities for Feynman graph amplitudes, dubbed Schouten identities, valid at fixed integer value of the dimension d is proposed. The identities are then used in the case of the two-loop sunrise graph with arbitrary masses for recovering the second-order differential equation for the scalar amplitude in d = dimensions, as well as a chained set of equations for all the coefficients of the expansions in d. The shift from d tod dimensions is then discussed. 0 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license Funded by SCOAP 3. Keywords: Feynman graphs; Multi-loop calculations; Schouten identities. Introduction The Feynman integrals associated to the two-loop self-mass Feynman graph of Fig., usually referred to as sunrise, have been widely studied in the literature within the framework of the integration by parts identities [,], and it is by now well known that they can be expressed in terms of four Master Integrals M.I.s [3], which satisfy a system of four first-order coupled differential equations [] equivalent to a single fourth-order differential equation for any of the * Corresponding author. addresses: ettore.remiddi@bo.infn.it E. Remiddi, tancredi@physik.uzh.ch L. Tancredi / 0 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license Funded by SCOAP 3.

2 3 E. Remiddi, L. Tancredi / Nuclear Physics B Fig.. The two-loop sunrise. Master Integrals. Several numerical approaches to the numerical solution of the equations with satisfactory degree of precision have been worked out see for instance [5], but a complete treatment of the general case with three different masses in d = dimensions is still missing. In the equal mass case the number of independent Master Integrals reduces to two, so that the two by two first-order system of differential equations can be rewritten as a single second-order differential equation for one of the Master Integrals, say the full scalar amplitude see below. In Ref. [6] it is shown how to build the analytic solution of that equation in terms of elliptic integrals, both for d = and d = ; the two cases, related by the Tarasov s shifting relations [7], are very similar, with the d = case just marginally simpler than the d = case. The analytic solution provides with the necessary information for writing out very precise and fast converging expansions for the accurate numerical evaluation of the two M.I.s [8]. More recently, an interesting paper [9] has shown, by using algebraic geometry arguments, that in d = dimensions the full scalar amplitude satisfies a second-order differential equation also in the different mass case. The equation was then solved in [0] by suitably extending the method of [6]; let us observe here that the analytic solution of the second-order differential equation is equivalent to the analytic knowledge of two of the four Master Integrals of the sunrise with different masses. The problem of extending the approach to d =, which is the physically relevant case, remains, as the straightforward use of the Tarasov s dimension-shifting relations is unfortunately not sufficient. Indeed, as will be shown in this paper, by explicitly working out the shifting relations one finds that any of the four Master Integrals at d dimensions can be expressed as a combination of all the four Master Integrals at d dimensions and of the first terms of their expansion in d, while the results of [9] give only two of the four Master Integrals at exactly d = dimensions, but no other information on the remaining Master Integrals and their expansion in d. In this paper we introduce a family of particular polynomials in the scalar products of the vectors occurring in the Feynman integrals, dubbed Schouten polynomials, which have the property of vanishing at some fixed integer value of the dimension d. By using those polynomials one can introduce an ad hoc set of amplitudes, from which one can at least in principle extract an independent set of new amplitudes which vanish in a non-trivial way see below at that value of d say at d = N for definiteness. If those new amplitudes are expressed in terms of the previously chosen set of Master Integrals, their vanishing gives a set of relations between the Master Integrals, valid at d = N, which we call Schouten identities. Alternatively, one can introduce a new set of Master Integrals including as new Master Integrals some of the independent amplitudes vanishing at d = N, write the system of differential equations satisfied by the new set of Master Integrals and expand them recursively in powers of d N around d = N. As some of the new Master Integrals vanish at d = N, the system of equations takes a simpler block structure. The pattern is very general, and applies in principle to the integrals of any Feynman graph. We work out explicitly the case of the sunrise amplitudes at d = with different masses, finding the existence of two independent Schouten identities, i.e. of two independent relations between

3 E. Remiddi, L. Tancredi / Nuclear Physics B the usual Master Integrals at d =, or, which is the same, we can introduce a new set of Master Integrals, consisting of two conventional Master Integrals say the full scalar amplitude and another M.I. and two new Master Integrals vanishing at d =. The system of differential equations satisfied by the new set of Master Integrals can then be expanded in powers of d. At zeroth-order we find a two by two system for the two conventional M.I.s the other two Master Integrals vanish, equivalent to the second-order equation found in [9], while at first-order in d we find in particular two relatively simple equations for the first terms of the expansion of the two new M.I.s, in which the zeroth-orders of the two conventional M.I.s appear as non-homogeneous known terms. One can move from d to the physically more interesting d case by means of the Tarasov s shifting relations; it is found that for obtaining the zeroth-order term in d of all the four M.I.s of the old or of the new set at d one needs, besides the zeroth-order term in d of the two old M.I.s at d, also the first term in d of the new M.I.s. The plan of the paper is as follows: in Section we introduce the Schouten polynomials for an arbitrary number of dimensions, while their applications to Feynman Amplitudes is discussed in Section 3. In Section we show how, by using the Schouten identities, a new set of Master Integrals can be found, whose differential equations in d = take an easier block form and can be therefore re-casted see Section 5 as a second-order differential equation for one of the Masters. In Section 6 we show how the results at d can be recovered from those at d through Tarasov s shifting relations. Finally, in Section 7, which is somewhat pedagogical, we present a thorough treatment of the imaginary parts of the Master Integrals in d = and d = dimensions. Many lengthy formulas and some explicit derivations can be found in Appendices A D at the end of the paper.. The Schouten polynomials As an introduction, let us recall that in d = dimensions one cannot have more than linearly independent vectors; indeed, given five vectors v α, a μ, b ν, c ρ, d σ in four dimensions they are found to satisfy the following relation v μ ɛa,b,c,d a μ ɛv,b,c,d b μ ɛa,v,c,d c μ ɛa,b,v,d d μ ɛa,b,c,v = 0,. where ɛ μνρσ is the usual Levi-Civita tensor with four indices, with ɛ 3 =, etc., and following the convention introduced in the program SCHOONSCHIP [] we use ɛa,b,c,d = ɛ μνρσ a μ b ν c ρ d σ.. Eq.. is known as the Schouten identity []; by squaring it, one gets a huge polynomial, of fifth-order in the scalar products of all the vectors. Due to Eq.., that polynomial vanishes in d = dimensions and a fortiori for any integer dimension d ; note however that the polynomial does not vanish identically for any arbitrary value of the dimension; as Eq.. is valid only when d, for d>the polynomial is not bound to take a vanishing value. As an extension or rather a simplification of Eq.., consider now the quantity ɛa,b = ɛ μν a μ b ν,.3 where ɛ μν is the Levi-Civita tensor with two indices defined of course by ɛ = ɛ =, ɛ = ɛ = 0, and a μ, b ν are a couple of two-dimensional vectors. By squaring it, Eq..3 gives at once

4 36 E. Remiddi, L. Tancredi / Nuclear Physics B ɛ a, b = a b a b,. where a, b are the squared moduli of the vectors a μ, b ν and a b their scalar product. So far, all the quantities introduced in Eqs..3,. are in d = dimensions. If the dimension d takes the value of any non-vanishing integer less than i.e. if d =, the r.h.s. of Eq..3 vanishes,and so does the r.h.s. of Eq.. as well.at this point we define the Schouten polynomial P d; a,b as P d; a,b = a b a b,.5 where the r.h.s. is formally the same r.h.s. of Eq.., but the two vectors a μ, b ν are assumed to be d-dimensional vectors, with continuous d. To emphasize that point, we have written d within the arguments of the Schouten polynomial, even if d does not appear explicitly in the r.h.s. of Eq..5. By the very definition, at integer non-vanishing dimension d< i.e. at d =, P d; a,b vanishes, P ; a,b = 0, as can be also verified by an absolutely trivial explicit calculation. Following the elementary procedure leading to Eq..5,givenind = 3 dimensions any triplet of vectors a μ, b ν, c ρ we consider ɛa,b,c = ɛ μνρ a μ b ν c ρ, where ɛ μνρ is the Levi-Civita tensor with three indices defined as usual by ɛ 3 =, etc. and then we evaluate its square.6.7 ɛ a,b,c= a b c a b c b a c c a b + a bb ca c..8 We then define the Schouten polynomial P 3 d; a,b,c as P 3 d; a,b,c = a b c a b c b a c c a b + a bb ca c,.9 where again the r.h.s. is formally the same as in Eq..8, but the three vectors a μ, b ν, c ρ are assumedtobed-dimensional vectors, with continuous d. By construction, P 3 d; a,b,c vanishes at d = and at d =, P 3 ; a,b,c = 0, P 3 ; a,b,c = 0..0 Needless to say, the procedure can be immediately iterated to any higher dimension, generating Schouten polynomials involving four vectors and vanishing in d =,, 3 dimensions, or involving five vectors and vanishing in d =,, 3, dimensions, corresponding, up to a constant numerical factor, to the square of Eq.., etc. As it is apparent from the previous discussion, the Schouten polynomials generated by a given set of vectors are nothing but the Gram determinants of the corresponding vectors; we prefer to refer to them as Schouten polynomials to emphasize that they vanish in any integer dimension d less than the number of the vectors. In the actual physical applications, as one is interested mainly in the d limit of Feynman graph amplitudes, one can reach d = starting from a different value of d and then moving to d = by means of the Tarasov s shifting relations [7]. As the shift relates values of d differing

5 E. Remiddi, L. Tancredi / Nuclear Physics B by two units, the d = Schouten identities, easily established for any amplitude in which at least two vectors occur, are of no use. The next simplest identities are at d = and occur with any amplitude involving at least three vectors. That is the case of the two-loop self-mass graph the sunrise, which we will study in this paper in the arbitrary mass case. 3. The Schouten identities for the Sunrise graph We discuss in this section the use of the Schouten polynomial P 3 d; a,b,c in the case of the sunrise, the two-loop self-mass graph of Fig.. The external momentum is p and the internal masses are m, m, m 3. We use the Euclidean metric, so that p is positive when spacelike; sometimes we will use also s = W = p, so that the sunrise amplitudes develop an imaginary part when s = W>m +m +m 3, the threshold of the Feynman graph. We write the propagators as D = q + m, D = q + m, D 3 = p q q + m 3, 3. and define the loop integration measure, in agreement with previous works, as: D d q = d d q, 3. Cd π d with so that Cd = π d/ Ɣ 3 d, 3.3 C = π and C =. 3. With that definition the tadpole Td,mreads Td; m = D d q q + m = m d d d. 3.5 In this paper we will use the double tadpoles Td; m,m = D d q D d q, 3.6 D D together with the similarly defined Td; m,m 3, Td; m,m 3, and the four amplitudes S d; p = D d q D d q, D D D 3 S d; p = d dm S d; p = D d q D d q D D, D 3 S d; p = d dm S d; p = S 3 d; p = d dm S d; p = 3 D d q D d q D D D, 3 D d q D d q D D D3. 3.7

6 38 E. Remiddi, L. Tancredi / Nuclear Physics B All those amplitudes depend on the three masses m, m, m 3, even if the masses are not written explicitly in the arguments for simplicity. The four amplitudes are equal, when multiplied by an overall constant factor π, to the four M.I.s used in []. Sd; p, in particular, is the full scalar amplitude already referred to previously. Those amplitudes were chosen in [] as M.I.s for the sunrise problem, and in the following they will be sometimes referred to as the conventional M.I.s. We can now introduce the Schouten amplitudes defined, for arbitrary d,as Z d; n,n,n 3,p = D d q D d P 3 d; p,q,q q, 3.8 D n Dn Dn 3 3 where the n i are positive integer numbers and P 3 d; p,q,q is the Schouten polynomial defined in Eq..9. The convergence of the integrals, for a given value of d, depends of course on the powers n i, as the Schouten polynomial in the numerator contributes always with four powers of the loop momenta q and q. We are interested here in the d = case. If the Schouten amplitude is convergent at d =, due to the second of Eqs..0, it is also vanishing at d =, i.e. Z; n,n,n 3,p = 0. Note that in the massive case all the integrals we are considering are i.r. finite, therefore the divergences can only be of u.v. nature. As one can express any sunrise Feynman amplitude in terms of a valid set of M.I.s, we will write in the following a few Schouten amplitudes in terms of the conventional M.I.s given in Eqs A few explicit results are now listed: Z d; p = Zd;,, d [ = d p + d 3 m + m + 3] m S d; p d p + m 6 m S d; p d + p 3m + m + 3 3m m S d; p d + p 3m + 3m + 3 m m 3 S 3 d; p d d [ + Td; m,m + Td; m,m 3 Td; m,m 3 ], 3.9 Z d; p = Z d;,,,p d [ = d p + d 3 m m + 3] m S d; p d + p + m 3m + 3 3m m S d,p d p + m 6 m S d; p d + p + 3m 3m + 3 m m 3 S 3 d; p d d [ + Td; m,m Td; m,m 3 + Td; m,m 3 ], 3.0

7 Z 3 d; p = Z d;,,,p E. Remiddi, L. Tancredi / Nuclear Physics B d [ = d p + d 3 m + m 3] m S d; p d + p + m + 3m 3 3m m S d; p d + p + 3m + m 3 3m m S d; p d p + m 6 3 m 3 S 3 d; p d d [ + Td; m,m + Td; m,m 3 + Td; m,m 3 ], 3. Z d;,,,p = d d [ d 3S d; p + m S d; p + m S d; p + m 3 S 3 d; p ]. 3. Some comments are in order. Elementary power counting arguments give N = n + n + n 3 powers of the integration momenta in the denominator independently of d and, in d = dimensions, all together eight powers in the numerator see Eq. 3.8 for the definition of the integrals, so that the minimum value of N necessary to guarantee the convergence is N = 0. In the case of Zd;,,,p of Eq. 3. N =, more than the minimum required value N = 0; therefore the integrals in the loop momenta q, q do converge, so that the vanishing of P 3 d; p,q,q in the numerator at d = and, as a matter of fact at d = as well does imply the vanishing of the whole amplitude. The explicit result, Eq. 3., shows indeed that the amplitude vanishes at d = and d =, but that is due to an overall factor d d, so that Eq. 3. does not give any useful information. This pattern the vanishing at d = of the amplitudes with P 3 d; p,q,q in the numerator and N>0 due to the appearance of an overall factor d is of general nature, and showed up in all the cases which we were able to check needless to say, the algebraic complications increase quickly with the powers of the denominators. The Z i d; p,eqs.3.9, 3.0, 3., are more interesting; in this case, N = 0, which is the minimum value needed to guarantee convergence in d = dimensions, so that those amplitudes are expected to vanish at d = and therefore also at d = as a consequence of the vanishing of P 3 d; p,q,q at d =, d =, see Eqs..0. The vanishing at d = is trivially given by the overall factor d in d = the minimum value of N to guarantee convergence is N = 8, while in the integrals we are now considering N = 0, but the vanishing at d = is totally non-trivial, providing new and so far not known relations between the four conventional M.I.s Sd; p, S i d; p at d =. Any of the three amplitudes Z i d; p can obviously be obtained from the others by a suitable permutation of the three masses, as immediately seen from their explicit expression. When summing the three relations, one obtains Z d; p + Z d; p + Z 3 d; p = d d p S d; p, 3.3

8 350 E. Remiddi, L. Tancredi / Nuclear Physics B showing that at d = they are not independent from each other; indeed, if one takes as input Z ; p = 0 and Z 3 ; p = 0, the previous equation gives Z ; p = 0, showing that the condition Z ; p = 0 depends on the other two. When written explicitly, the vanishing of Z ; p = 0 and Z 3 ; p = 0 reads Z ; p = m m + 3 m S ; p + p + m 3m + 3 3m m S,p p 6 + m m S ; p + p + 3m 3m + 3 m m 3 S 3 ; p + 96 ln m m m 3 = 0, Z 3 ; p = m + m 3 m S ; p + p + m + 3m 3 3m m S ; p p + 3m + m 3 3m m S ; p + p 6 + m 3 m 3 S 3 ; p + 96 ln m 3 m m = The validity of identities Eqs. 3., 3.5 in d = has been checked with SecDec [3].Byusing the above relations, which hold identically in p, m, m, m 3, one can express two of the conventional M.I.s in terms of the other two, showing that, at d =, there are in fact only two independent M.I.s. As can be seen from Eqs. 3., 3.5, the relations between the M.I.s are not trivial in particular, none of the M.I.s vanishes at d = ; according to the definition Eq. 3.7 for space-like p they are all positive definite.. A new set of Master Integrals We have seen in the previous section that the conventional M.I.s in d = dimensions satisfy two independent conditions, written explicitly in Eqs. 3., 3.5, so that two of them can be expressed as a combination of the other two, which can be taken as independent. On the other hand, it is known that in the equal mass limit the Sunrise has two independent M.I.s in any dimension, including d = so that no other independent conditions can exist. It can therefore be convenient to introduce a new set of M.I.s, formed by two conventional M.I.s, say Sd; p, S d; p of Eq. 3.7, and two Schouten amplitudes, say Z d; p, Z 3 d; p of Eqs. 3.0, 3.. The advantage of the choice is that two conditions at d = take the simple form Z ; p = 0, Z 3 ; p = 0. The actual choice of the new M.I.s satisfying the above

9 E. Remiddi, L. Tancredi / Nuclear Physics B criteria is of course not unique a fully equivalent set could be for instance Sd; p, S d; p, Z d; p, Z d; p,etc.. In the new basis of M.I.s, the two discarded conventional M.I.s are expressed as P p,m,m,m 3 m S d; p = { m [ m d 3 m + m 3 m p ] d p p + m } 3 S d; p + P p,m,m,m 3 m S d; p 8 p + m d 3 Z d; p p + 3m d 3m + 3 m Z3 d; p d m m Td; m,m d p m + m + m 3 Td; m,m 3 d + p + m m + m 3 Td; m,m 3,. P p,m,m,m 3 m 3 S 3 d; p = { m 3[ m d 3 m m + 3 m p ] d p p + m } S d; p + P p,m 3,m,m m S d; p p + 3m d + m 3 3m Z d; p 8 p + m d Z3 d; p d p m + m + m 3 Td; m,m d m 3 m Td; m,m 3 + d p + m + m m 3 Td; m,m 3,. where Pp,m,m,m 3 is the polynomial P p,m,m,m 3 = p + m + m 3 m p 3m + m + m 3 + m m + m m 3 m m 3..3 Note that Pp,m,m,m 3, which is symmetric in the last two arguments, P p,m,m,m 3 = P p,m,m 3,m,. occurs with different arguments in different places.

10 35 E. Remiddi, L. Tancredi / Nuclear Physics B By substituting the above expressions in the differential equations for the conventional M.I.s as given, for instance, in Ref. [], one obtains the new equations P p,m,m,m 3 p d dp S d; p = p + m [ p m + m + m 3 + d p + m m ] m 3 S d; p Q p,m,m,m 3 m S d; p + 3p + 3m d + m 3 m Z d; p + 3p + 3m d m + 3 m Z3 d; p d + p + m m + m 3 Td; m,m d + p + m + m m 3 Td; m,m 3 d p + m Td; m,m 3,.5 D p,m,m,m 3 P p,m,m,m 3 p d dp S d; p [ d = p + m m 3 m P 0 p,m,m,m 3 d P 0 [ d + P d 3 [ 0 + P d [ d + + d d d P p 0,m,m,m 3 P 0 p ],m 3,m,m S d; p m [ d m p 0,m,m,m 3 P p ],m,m,m 3 S d; p p,m,m,m 3 Z d; p + P 0 p,m,m 3,m Z3 d; p ] P P [ d P 0 p,m,m,m 3 P p ],m,m,m 3 Td; m,m p,m,m 3,m P p ],m,m 3,m Td; m,m 3 p,m,m,m 3 p,m,m,m 3 + P p,m,m 3,m ] Td; m,m 3,.6 P p,m,m,m 3 p d dp Z d; p d d [ = p m 8 m p + m + m 3 m + d p + m m 3 m p + m m + ] m 3 S d; p

11 E. Remiddi, L. Tancredi / Nuclear Physics B d d p P p,m,m,m 3 m S d; p d + P p,m,m,m 3 Z d; p + p d p + 3m 3m + m 3 Z3 d; p d d + p d d + p 8 d d p 8 m m Td; m,m p m + m + m 3 Td; m,m 3 p + m m + m 3 Td; m,m 3,.7 P p,m,m,m 3 p d dp Z 3 d; p d d [ = p m 8 m 3 p + m m + 3 m + d p + m m 3 m p + m + m ] m 3 S d; p d d p P p,m 3,m,m m S d; p + p d p + 3m + m 3m 3 Z d; p d + P p,m,m 3,m Z3 d; p d d + p p m + m + m 3 Td; m,m 8 d d + p d d p 8 In the above equations, m m 3 Td; m,m 3 p + m + m m 3 Td; m,m 3..8 D p,m,m,m 3 = p + m + m + m 3 p + m m + m 3 p + m + m m 3 p + m m m 3.9 is the product of all the threshold and pseudo-threshold factors already present in [], Q p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 + p m + m + m 3 + 3p,.0 while Pp,m,m,m 3 is the polynomial previously defined in Eq..3. Finally the P n ij p,m,m,m 3 are also polynomials depending on p and the masses; their explicit and sometimes lengthy expression is given in Appendix A. Note that a same polynomial can occur in different equations with a different permutation of the masses within its arguments.

12 35 E. Remiddi, L. Tancredi / Nuclear Physics B We want to stress here an important aspect of the last two equations, Eqs..7,.8, namely the presence of an overall factor d in the r.h.s., which plays an important role in the expansion in powers of d discussed in the next subsection... The expansion of the equations around d = Let us start off by expanding all M.I.s in powers of d around d =, S d; p = S ; p + d S,p +, S d; p = S ; p + d S,p +, Z d; p = Z ; p + d Z,p +, Z 3 d; p = Z 3 ; p + d Z 3,p +.. Due to the overall factor d in the r.h.s., at 0th order in d the differential equations.7,.8 become d dp Z ; p = 0, d dp Z 3 ; p = 0,. showing that Z ; p, Z 3 ; p must be constants. But we know from Eqs. 3., 3.5 the actual value of that constant the two functions vanish identically, Z ; p = 0, Z 3 ; p = 0, so that at 0th order in d the differential equations.5,.6 for S; p, S ; p become P p,m,m,m 3 p d dp S ; p = p + m p m + m + m 3 S ; p Q p,m,m,m 3 m S ; p + [ p + m m ln + ] m 8 m m m 3 m 3 ln,.3 3 D p,m,m,m 3 P p,m,m,m 3 p d dp S ; p = P 0 0 p,m,m,m 3 S,p P 0 p,m,m,m 3 S,p [ P 8 p m,m,m,m 3 ln + P m p m,m,m 3,m ln 3 m m p m P p ],m,m,m 3,. completely decoupled, obviously, from the trivial equations for Z ; p, Z 3 ; p. See the previous section and Appendix A for the explicit expression of the polynomials.

13 E. Remiddi, L. Tancredi / Nuclear Physics B Going now one order higher in the expansion, one finds that the first-order terms in d of the Z i d; p satisfy the equations P p d,m,m,m 3 dp Z ; p = m m p + m + m m 3 S ; p P p,m,m,m 3 m S ; p + [ p + m m 3 ln + ] m 3 m m m m ln,.5 m 3 P p d,m,m,m 3 dp Z 3 ; p = m m 3 p + m m + m 3 S ; p P p,m 3,m,m m S ; p + [ p + m m ln + ] m 3 m m m 3 m 3 ln..6 3 It is to be noted that Z ; p, Z 3 ; p do not appear in the r.h.s. of Eqs..5,.6, which contains only S; p and S ; p, to be considered known once Eqs..3,. for the 0th orders in d have been solved. Eqs..5,.6, indeed, are absolutely trivial when considered as differential equations, as they contain only the derivatives of Z ; p, Z 3 ; p, and can therefore be solved by a simple quadrature. Knowing Z ; p, Z 3 ; p, one can move to the differential equations for S ; p, S ; p which we don t write here for the sake of brevity; they involve Z ; p, Z 3 ; p as known inhomogeneous terms, and form again a closed set of two differential equations, decoupled from the equations for the other two M.I.s, as at 0th order in d. Thanks to the overall factor d in the r.h.s. of Eqs..7,.8, the pattern a quadrature for Z k ; p, Z k 3 ; p and a closed set of two differential equations for S k ; p, S k ; p is completely general, and can be iterated, at least in principle, up to any required order k in d. 5. Second-order differential equation for Sd;p We go back now to the system of differential equations.5,.6, for obtaining a secondorder differential equation for Sd; p. We can use Eq..5 in order to express S d; p in function of Sd; p and of its derivative, dsd; p /dp. By substituting this expression into Eq..6 we can then derive a second-order differential equation for Sd; p only, which however still contains Z d; p and Z 3 d; p in the inhomogeneous part: m

14 356 E. Remiddi, L. Tancredi / Nuclear Physics B A p d,m,m,m 3 dp S d; p = 0, + [ A 0 + d 3 [ A 0 3 p,m,m,m 3 + d A p ] d,m,m,m 3 dp S d; p p,m,m,m 3 + d A p ],m,m,m 3 S d; p d 3 [ + A p,m,m,m 3 Z d; p + A p,m,m 3,m Z3 d; p ] d + d [ A 5 p,m,m,m 3 + d A 5 p ],m,m,m 3 Td; m,m + d [ A p,m,m 3,m + d A p ],m,m 3,m Td; m,m d [ A p,m,m 3,m + d A p ],m,m 3,m Td; m,m 3 where A p,m,m,m 3 = p Dp,m,m,m 3 Pp,m,m,m 3, with Dp,m,m, m 3 and Pp,m,m,m 3 being the usual polynomials defined by Eqs..3,.9. The A n j p,m,m,m 3 are also polynomials which depend on the three masses and on p,but do not depend on the dimensions d. Their explicit expressions, as usual quite lengthy, can be found in Appendix B. The equation above is exact in d but contains, besides Sd; p and its derivatives, also Z d; p and Z 3 d; p as inhomogeneous terms. Nevertheless, recalling once more that Z ; p = Z 3 ; p = 0, we can expand Eq. 5. in powers of d and obtain at leading order in d a second-order differential equation for S; p only: A p d,m,m,m 3 dp S ; p + A 0 p d,m,m,m 3 A 0 3 p,m,m,m 3 S ; p + A 5 + A 5 p,m,m 3,m + A 5 p,m,m 3,m + A 5 p,m,m,m 3 ln m m 3 [ + A 5 5 dp p,m,m,m 3 p,m,m 3,m ln m m S ; p 5. ] = 0, 5. where we made use of the relation Eq. B.8 of Appendix B. We compared Eq. 5. with the second-order differential equation derived in [9], finding perfect agreement. Eq. 5. has been solved in reference [0] in terms of one-dimensional integrals over elliptic integrals. Upon inserting the result in Eq..3 one can obtain S ; p in terms of S; p and ds; p /dp. Inserting then S; p and S ; p in Eqs..5,.6, one obtains by quadrature the firstorder terms, Z ; p and Z 3 ; p, of the expansion in d of Z d; p and Z 3 d; p. Having these results on hand, we can now consider the first-order in d of Eq. 5., which is now a second-order differential equation for S ; p only, with known inhomogeneous terms namely S; p, Z ; p and Z 3 ; p. Proceeding in this way, at least in principle, the whole procedure can be iterated up to any required order in d.

15 E. Remiddi, L. Tancredi / Nuclear Physics B Shifting relations from d = tod = dimensions In the previous sections we have shown how to use the Schouten identities for writing the differential equations for the M.I.s of the massive sunrise at d = dimensions in block form, and outlined the procedure for obtaining iteratively all the coefficients of the expansion in d of the four M.I.s starting from a second-order differential equation for S; p, the leading term of the expansion. The physically interesting case corresponds however to the expansion of the M.I.s for d ; we have therefore to convert the information given by the expansion at d in useful information at d. As it is well known, quite in general one can relate any Feynman integral evaluated in d dimensions to the very same integral evaluated in d dimensions by means of the Tarasov s shifting relation [7]. This dimensional shift is achieved by acting on the Feynman integral with a suitable combination of derivatives with respect to the internal masses. In the case of the conventional M.I.s of the sunrise graph, as defined in Eq. 3.7, the shifting relations read: S d ; p = d 6 S d; p, S i d ; p = d 6 S i d; p, i =,, 3, 6. where the differential operator takes the form: = m m + m m 3 + m m Carrying out the derivatives in the integral representation for the four M.I.s of Eq. 3.7, one obtains a combination of integrals which are still related to the sunrise graph. They can be expressed in terms of the full set of M.I.s in d dimensions by full set we mean the four M.I.s and the tadpoles; one obtains in that way a set of four equations which explicitly relate the four M.I.s of the sunrise graph evaluated in d dimensions to suitable combinations of the same integrals and of the tadpoles evaluated in d dimensions. In that direct form the shifting relations would be of no practical use in our case, as they might give the M.I.s at d in terms of those less known at d. It is however straightforward to invert the system and, in this way, to obtain the inverse shifting relations, expressing the four M.I.s in d + dimensions in function of those in d dimensions. In addition, we can also use Eqs..,. for expressing S d; p and S 3 d; p, in terms of Sd; p, S d; p and Z d; p, Z 3 d; p. As a result one arrives at expressing any of the four conventional M.I.s Sd +,p, S i d + ; p, i =,, 3, as a linear combination whose coefficients depend and in a non-trivial way on d and the kinematical variables of the problem of the new M.I.s Sd; p, S d; p, Z d; p and Z 3 d; p and the tadpoles. Indicating for simplicity the four conventional M.I.s with M i d and with N j d the four new M.I.s and the tadpoles, and ignoring for ease of notation all the kinematical variables, the inverse shifting relations can be written as M i d + = j C i,j dn j d. 6.3

16 358 E. Remiddi, L. Tancredi / Nuclear Physics B Given a relation of the form Fd + = Gd, by expanding around d = one has, quite in general Fd + = Gd = p d n F n, n=r p d n G n, n=r where r, the first value of the summation index, can be negative as it is the case in a Laurent expansion, so that F n = G n. In the case of the inverse shift Eq. 6.3, one has that the coefficients of the expansion of the conventional M.I.s in d for d are completely determined by those of the expansion in d for d of the new M.I.s, discussed in the previous sections, and of the tadpoles expanding around d = the two sides of Eq. 6.3 requires also the expansion of the coefficients C i,j d, but that is not a problem once the inverse shift has been written down explicitly. The explicit formulas of the direct or inverse shifting relations are easily obtained but very lengthy and we decided not to include them entirely here for the sake of brevity. For what follows, it is sufficient to discuss only the general features of one of the inverse shifting relations, namely the relation expressing Sd + ; p in terms of Sd; p, S d; p and Z d; p, Z 3 d; p. Keeping for simplicity only the leading term of the expansion in d of the coefficients we find: S d + ; p = [ C p ],m,m,m 3 + Od S d; p + [ C p ],m,m,m 3 + Od S d; p [ ] + d C p,m,m,m 3 + O Z d; p [ ] + d C 3 p,m,m,m 3 + O Z 3 d; p + [ C 0 p ],m,m,m 3 + Od Td; m,m + [ C 0 5 p ],m,m,m 3 + Od Td; m,m + [ C 0 6 p ],m,m,m 3 + Od Td; m,m. 6. In the formula above the Cp,m,m,m 3, C i p,m,m,m 3, are ratios of suitable polynomials which, as usual, depend on p and on the three masses but, most importantly, do not depend on the dimensions d. The explicit expressions for Cp,m,m,m 3, C p,m,m,m 3, C p,m,m,m 3 and C 3 p,m,m,m 3, which will also be used in the following, can be found in Appendix C, Eqs.C. C.5. Note anyway that: C 3 p,m,m,m 3 = C p,m,m 3,m.

17 E. Remiddi, L. Tancredi / Nuclear Physics B By writing the expansion of Sd + ; p at d as S d + ; p = n S n ; p d n, 6.5 and then expanding Eq. 6. at d, one recovers the expression of the coefficients S n ; p in terms of the coefficients of the expansion of the four M.I.s and the tadpoles in d. A few observations are in order. Eq. 6. exhibits an explicit pole in /d only in the coefficients of Z d; p and Z 3 d; p ; recalling once more that at d = both Z ; p and Z 3 ; p are identically zero, see Eqs. 3., 3.5, it is clear that these poles will not generate any singularity of Sd; p as d. On the other hand, the tadpoles in the r.h.s. of Eq. 6. do generate polar singularities of Sd+; p ; recalling Eqs. 3.5, 3.6 and by using the lengthy explicit form of the coefficients which we did not write for brevity one finds immediately S ; p = m + m + m 3, 8 S ; p = 3[ p + 6 m + m + ] m 3 8[ m ln m + m ln m + m 3 ln m 3], 6.6 formulas already known for a long time in the literature []. As a second observation, let us look at the zeroth-order term S 0 ; p of Sd; p in d, i.e. the zeroth-order term in d of Eq. 6.. We have already commented the apparent polar singularity /d in the coefficients of Z d; p and Z 3 d; p, actually absent because Z ; p and Z 3 ; p are both vanishing. But due to the presence of the /d polar factor, in order to recover the zeroth-order term S 0 ; p, one needs, besides S; p, S ; p,also the first-order of the corresponding expansion of Z d; p and Z 3 d; p, namely Z ; p and Z 3 ; p obtained, in our approach, from the systematic expansion of the differential equations, see Eqs..5,.6 or Section 5. The complete expression of S 0 ; p, which is rather cumbersome, is given by Eq. C. of Appendix C. The corresponding formulas for the other three M.I.s, i.e. the S i d + ; p, can be obtained directly from the authors. 7. The imaginary parts of the Master Integrals In this section, which is somewhat pedagogical, we discuss the relationship between the imaginary parts of the M.I.s at d = and d = dimensions, as a simple but explicit example of functions exhibiting the properties described in the previous sections. At d = the Cutkosky Veltman rule [,5] gives for Sd; p, as defined by the first of Eqs. 3.7, π Im S ; W = N b 3 b db where the following notations were introduced: R b; b,b,b 3,b, 7.

18 360 E. Remiddi, L. Tancredi / Nuclear Physics B N = /, p = W, W m + m + m 3, m m 3 = b m + m 3 = b W m = b 3 W + m = b, R b; b,b,b 3,b = b b b b b 3 bb b. 7. We have the relation where R b; b,b,b 3,b = R b,m,m 3 R W,b,m, 7.3 R a,b,c= a + b + c ab ac bc, 7. is the familiar invariant form appearing in the two-body phase space, showing that the system of the three particles, whose masses enter in the definition of R b; b,b,b 3,b, can be considered as the merging of a two-body system of total energy b and masses m, m 3 with a two-body system of total energy W and masses b, m. According to Eqs. 3.7, fori =,, 3 π Im S i ; W = d dm i π Im S ; W ; 7.5 the integral representation Eq. 7., however, is of no use for obtaining Im S i ; W through a direct differentiation due to the appearance of end point singularities. It is more convenient to use Eq. D.3 of Appendix D, so that Eq. 7. becomes π Im S ; W = N b b b 3 b K w, 7.6 where Kw is the complete elliptic integral of the first kind, Eq. D.6, and w = b b b 3 b b b b 3 b = W + m + m m 3 W + m m + m 3 W m + m m 3 W m m + m 3 W + m + m + m 3 W + m m m 3 W m + m m 3 W m m + m 3, b b b 3 b = W + m + m + m 3 W + m m m 3 W m + m m 3 W m m + m Let us observe, in passing, that, even if Im S, W and, more generally Sd; p as well is obviously symmetric in the three masses m, m, m 3, the symmetry is not explicitly shown by the integral representation Eq. 7., while the manifest symmetry is restored in Eqs. 7.6, 7.7. One can now use Eq. 7.6 and Eq. D. to carry out the derivative with respect to the masses in Eq. 7.5; the result reads m i π Im S ; W = N m b3 b b b b 3 b b b [ m m m 3 + m m m3 + m m 3 W + m W K w P W,m,m,m 3 E w ], 7.8

19 E. Remiddi, L. Tancredi / Nuclear Physics B π Im S ; W = N m b3 b b b b 3 b b b [ m m m 3 + m m m3 + m m 3 W + m W K w P W,m,m,m 3 E w ], 7.9 π Im S 3 ; W = N m 3 b3 b b b b 3 b b b [ m 3 m3 m + m 3m m3 3 + m m W + m 3 W K w P W,m 3,m,m E w ], 7.0 where w is given by the first of Eqs. 7.7, Pp,m,m,m 3 is the polynomial already introduced in Eq..3, symmetric in the last two arguments, and Ew is the complete elliptic integral of the second kind, see Eq. D.7. Eqs. 7.6, 7.8, 7.9, 7.0 express the four quantities Im S; W,ImS i, W, i =,, 3 in terms of just two functions, the elliptic integrals Kw, Ew ; therefore, the four imaginary parts cannot be all linearly independent. It is indeed easy to check that they satisfy the two equations m m + m 3 Im S ; W + W + m 3m + 3 3m m Im S, W 6 W + m m Im S ; W + W + 3m 3m + m 3 m 3 Im S 3 ; W = 0, m + m m 3 Im S ; W + W + m + 3m 3 3m m Im S ; W + W + 3m + m 3 3m m Im S ; W 6 W + m 3 m 3 Im S 3 ; W = 0, which are nothing but the imaginary parts of Z ; W, Z 3 ; W,Eqs.3., 3.5. As a further comment on the imaginary parts at d =, let us observe that they take a finite value at threshold, i.e. in the W m + m + m 3 limit. In that limit, indeed, b 3 b = m + m 3, and one finds b 3 b so that db R b; b,b,b 3,b b b b b b 3 b db b b b 3 b = π b b b b, π Im S ; W N W m +m +m 3 m m m 3 m + m + m 3. 7.

20 36 E. Remiddi, L. Tancredi / Nuclear Physics B The extension to the value at threshold of the Im S i, W is similar, even if requiring one more term in the expansion due to the presence of the denominator /b 3 b in their definitions, Eqs. 7.8, 7.9, 7.0. The threshold values are with π Im S ; W N W m +m +m m + m + m 3 m + m + m 3 m m m 3 m + m + m 3, π Im S ; W N W m +m +m m m m 3 m + m + m 3 m m m 3 m + m + m 3, π Im S 3 ; W N W m +m +m m m m 3 m + m + m 3 m m m 3 m + m + m At d = the imaginary part of Sd; p, by using the same notation as in Eq. 7., is given by π Im S ; W = N N = 8W. b 3 b db R b; b,b,b 3,b, 7.3 b 7. At variance with the d = case, the Im S i, W can be obtained at once by differentiating with respect to the masses the previous integral representation for Im S; W. The result can be most conveniently expressed in terms of the four independent integrals I,W, I0,W, I,W, I,W, defined see Eq. D.3 and Appendix A for more details and the relation to the standard complete elliptic integrals through b 3 In,W= dbb n R b; b,b,b 3,b. 7.5 b An explicit calculation gives π Im S ; W [ = N b b b 3 b I,W 3 b b 3 b + b b 3 b + b b b + b b b 3 I0,W + b 3b + b b + b b 3 + b b + b b 3 + b b I,W ] b + b + b 3 + b I,W, 7.6

21 E. Remiddi, L. Tancredi / Nuclear Physics B π Im S ; W [ = N b b b b 3 W + b + b 3 m I,W + b + b b b 3 W b b + b b 3 + b b + b b 3 + b b m I0,W + b b 3 W + b + b 3 + b + b m I,W m I,W ], 7.7 π Im S ; W [ = N b3 b b b m 3 + b + b m I,W + b 3 b + b b + b b 3 + b b + b b 3 m + b b b + b 3 m 3 I0,W + b + b 3 + b + b m b b m 3 I,W m I,W ], 7.8 π Im S 3 ; W [ = N b3 b b b m + b + b m 3 I,W + b 3 b + b b + b b 3 + b b + b b 3 m 3 + b b b + b 3 m I0,W + b + b 3 + b + b m 3 b b m I,W m 3 I,W ]. 7.9 Again at variance with the d = case, the four imaginary parts are now combinations of four independent elliptic integrals, and therefore all independent of each other. Having recalled the main features of the imaginary parts of the M.I.s at d = and d = dimensions, we can look at the way the Tarasov s shifting relations work in their case. Let us start from the direct shift expressing the imaginary parts at d = in terms of those at d =. The d limit of the shifting relations is trivial, even if the relevant formulas are as usual rather lengthy. Keeping only the imaginary parts of the Master Integrals one finds for the M.I. S,p, with p = W m + m + m 3 where π Im S, W = ÃW, m,m,m 3 π Im S, W + BW,m,m,m 3 m π Im S, W + BW,m,m 3,m m π Im S + BW,m 3,m,m m 3 π Im S 3, W, W, 7.0 ÃW, m,m,m 3 = AW, m,m,m 3 + AW, m, m,m 3 + AW, m,m, m 3 + AW, m, m, m 3, BW,m,m,m 3 = BW,m,m,m 3 + BW,m, m,m 3 + BW,m,m, m 3 + BW,m, m, m 3,

22 36 E. Remiddi, L. Tancredi / Nuclear Physics B AW, m,m,m 3 = m m m 3 m + m + m 3 W m + m + m 3, BW,m,m,m 3 = m + m + m 3 AW, m,m,m Eq. 7.0 is relatively simple, and, when substituting in it the explicit values of Im S, W and Im S i, W, as given by Eqs ,Eq.7. is recovered. The same happens for Im S i, W, i =,, 3 as well. Conversely, one can look at the inverse formulas, giving the imaginary parts at d + in terms of the imaginary parts at d. For Im S, W, taking only the imaginary part at d = ofeq.6., one obtains: π Im S ; W = C W,m,m,m 3 π Im S ; W + C W,m,m,m 3 π Im S ; W + C W,m,m,m 3 π + C 3 W,m,m,m 3 π Im Z ; W Im Z 3 ; W, 7. where the C W,m,m,m 3, C i W,m,m,m 3 have been defined in the previous section, and their explicit expressions can be found in Eqs. C. C.5, ImS; W, Im S ; W are the imaginary parts of the corresponding Master Integrals at d =, while Im Z ; W,ImZ 3 ; W are the imaginary parts of the first term of the expansion in d of the corresponding functions, see Eqs.. let us recall once more that according to Eqs. 3., 3.5 Z ; p, Z 3 ; p vanish identically. An equation similar to Eq. 7. holds for Im S ; W ; we do not write it explicitly for the sake of brevity. The functions Im S; W,ImS ; W and Im S; W,ImS ; W are known, see Eqs. 7.6, 7.7 and Eqs. 7.6, 7.8; by combining Eq. 7. and the similar not written equation for Im S ; W, one can obtain the explicit values of Im Z ; W, Im Z 3 ; W. One finds ; W = N π π Im Z Im Z 3 [ W m m m I0,W+ I,W m 3 m W m ] I,W, 7.3 ; W = N [ W + m 3 6 m m I0,W+ I,W + m 3 m W m ] I,W. 7. From the previous equations and the same procedure giving Eqs. 7., 7. we obtain in particular the values at threshold Im Z π ; W N m m + m + m 3 W m +m +m 3, 6 m m 3 π Im Z 3 ; W N W m +m +m 3 6 m 3 m + m + m 3 m m. 7.5

23 E. Remiddi, L. Tancredi / Nuclear Physics B Im Z ; W can also be evaluated solving, by quadrature, the imaginary part of the differential equation.5, i.e.byevaluating ImZ ; W W d = C + dp Im dp Z ; p, where C is an integration constant and dz ; p /dp is obtained from Eq..5 itself. The constant C can be fixed, a posteriori, by requiring that the imaginary parts of the conventional M.I. vanish at threshold in d = dimensions, a condition which leads again to Eqs After many algebraic simplifications, one obtains for Im Z ; W π Im Z ; W = N 6 m m + m + m 3 m m 3 W + [ ds Fs,m,m,m 3 I0,s 6 m +m +m 3 Gs, m,m,m 3 I,s + Hs,m,m,m 3 I,s ], 7.6 where the three quantities F, G, H are all expressed in terms of the corresponding functions F, G, H by the relation Fs,m,m,m 3 = Fs,m,m,m 3 + Fs,m, m,m 3 and the explicit expressions of those functions are Fs,m,m,m 3 = m + m 3 + Fs,m,m, m 3 + Fs,m, m, m 3, m m 3 m + m + m 3 + m m + m m 3 s m + m + m 3, Gs, m,m,m 3 = m + m + m 3 + m m + m m 3 + m m 3 m m 3 [s m + m + m 3, ] Hs,m,m,m 3 = m m 3 [s m + m + m 3 ]. To carry out the integration, we use the integral representations Eq. D. for the elliptic integrals In,sand exchange the order of integration according to W m +m +m 3 ds = W m m +m 3 s m m +m 3 db R b, m,m 3 db R b; b,b,b 3,b W b+m ds R s,b,m, where Eq. 7.3 was used. The s integration is then elementary, giving only logarithms of suitable arguments and new square roots quadratic in b; a subsequent integration by parts in b removes

24 366 E. Remiddi, L. Tancredi / Nuclear Physics B those logarithms with some of the accompanying square roots, and the result is Eq. 7.3,asexpected. The same applies also for Im Z 3 ; W, whose value is obtained by simply exchanging m and m 3 in Eq Conclusions In this paper we introduced a new class of identities, dubbed Schouten identities, valid at fixed integer value of the dimensions d. We applied the identities valid at d = to the case of the massive two-loop sunrise graph with different masses, finding that in d = dimensions only two of the four Master Integrals M.I.s are actually independent, so that the other two can be expressed as suitable linear combinations of the latter. In the general case of arbitrary dimension d and different masses, the four M.I.s are known to fulfill a system of four first-order coupled differential equations in the external momentum transfer. The system can equivalently be re-phrased as a fourth-order differential equation for one of the M.I.s only. Using these relations we introduced a new set of four independent M.I.s, valid for any number of dimensions d, whose property is that two of the newly defined integrals vanish identically in d =. The new system of differential equations for this set of M.I.s takes then a simpler block form when expanded in d. Starting from this system, one can derive a second-order differential equation, exact in d, for the full scalar amplitude, which still contains the two integrals, whose value is zero at d =, as inhomogeneous terms. We verified that the zeroth-order of our equation corresponds to the equation derived in [9]. Our equations, once expanded in powers of d, can be used, together with the linear equations for the remaining three M.I.s, for evaluating recursively, at least in principle, all four M.I.s, up to any order in d. We then worked out explicitly the Tarasov s shifting relations needed to recover the physically more relevant value of the four M.I.s expanded in d at d starting from the expansion in d at d worked out in our approach. As an example of this procedure we discussed the relationship between the imaginary parts of the four M.I.s in d = and d =. The latter can be computed using the Cutkosky Veltman rule. We showed how in d = the imaginary parts of the four M.I.s can be written in terms of two independent functions only, namely the complete elliptic integrals of the first and of the second kind. The same is not true in d = dimensions, where four independent elliptic integrals are needed in order to represent the four imaginary parts. We then showed how the Tarasov s shift formulas relate the imaginary parts in d = and d = dimensions. Finally, we gave an explicit example of how the differential equations for the imaginary parts of the Master Integrals can be integrated by quadrature. Acknowledgements We are grateful to J. Vermaseren for his assistance in the use of the algebraic program FORM [7] which was intensively used in all the steps of the calculation, to T. Gehrmann for many interesting comments and discussions, and to F. Cascioli for proofreading the whole manuscript.

25 E. Remiddi, L. Tancredi / Nuclear Physics B L.T. wishes to thank A. von Manteuffel for his assistance with Reduze [8], with which the reduction to Master Integrals of the two-loop sunrise has been carried out, and G. Heinrich and S. Borowka for their help with SecDec [3]. This research was supported in part by the Swiss National Science Foundation SNF under the contract PDFMP-350. Appendix A. The polynomials of the first-order differential equations In this appendix we give the explicit expressions for the polynomials appearing in the firstorder differential equations in Section. All polynomials are functions of p and of the three masses m, m, m 3, while they do not depend on the dimensions d. P 0 0 p,m,m,m 3 = m m + m + m 3 m m + m 3 m + m m 3 m m m 3 m m m 3 + p m 8 + m m6 m m + m6 m 3m8 + m 3 m6 + m 3 m m 8m 3 m6 m 3 m + m 3 m m + m 3 m + m6 3 m 8m6 3 m 3m8 3 + p 0m 6 m m + m m 8m6 m 3 m + 6m 3 m + m 3 m + 6m 3 m 8m6 3 + p 6 m m m 6m m 3 m + m 3 m 6m 3 + 7p 8 m + p0, A. P 0 p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 m m 3 m 7m 6m m m 6m 3 m + m 3 m m 3 + p m 8 8m m6 + 78m m 56m6 m + 5m8 8m 3 m6 + 68m 3 m m 0m 3 m m + 0m 3 m6 + 78m 3 m 0m 3 m m 70m 3 m 56m6 3 m + 0m6 3 m + 5m8 3 + p m 6 m m + m m + m6 m 3 m + 60m 3 m m 6m 3 m + m 3 m 6m 3 m + m6 3 p 6 m m 3 3m m 3m m 3 + p 8 m + m + m p 0, A. P 0 p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 5m m m m m 3 m + m 3 m m 3 + p 8m 6 8m m + 0m m 0m6 8m 3 m + m 3 m m + 0m 3 m + 0m 3 m + 0m 3 m 0m6 3

26 368 E. Remiddi, L. Tancredi / Nuclear Physics B p 0m + 6m m 8m + 6m 3 m + 8m 3 m 8m 3 + p 6 6m + 0m + 3 0m + 9 p 8, A.3 P 0 p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 m m m 3 m + p 6m 0 + m m8 + 3m m6 m6 m + m8 m + m0 + m 3 m8 6m 3 m m6 + 6m 3 m m + 0m 3 m6 m 3m 3 m8 + 3m 3 m6 + 6m 3 m m 08m 3 m m + m 3 m6 m6 3 m + 0m6 3 m m + m 6 3 m + m8 3 m 3m8 3 m + m0 3 + p 33m 8 + 6m m6 0m m + m6 m + 5m8 + 6m 3 m6 m 3 m m 6m 3 m m m 3 m6 0m 3 m 6m 3 m m m 3 m + m6 3 m m6 3 m + 5m8 3 + p 6 5m 6 6m m + 3m m + 0m6 6m 3 m 6m 3 m m + 6m 3 m + 3m 3 m + 6m 3 m + 0m6 3 + p 8 33m m m + 0m m 3 m + m 3 m + 0m 3 + p 0 6m + 5m + 3 5m + p, A. P p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 5m m m m m 3 m + m 3 m m 3 + p m 0 8m m8 + 76m m6 80m6 m + 3m8 m m 0 8m 3 m8 + 0m 3 m m6 8m 3 m m + m 3 m6 m + m 3 m8 + 76m 3 m6 8m 3 m m 6m 3 m m 8m 3 m6 80m 6 3 m + m6 3 m m 8m6 3 m + 3m8 3 m + m8 3 m m0 3 + p m 8 m m + 88m6 m 5m8 + 5m 3 m m 5m 3 m m + m 3 m6 m 3 m 5m 3 m m + m 3 m + 88m6 3 m + m6 3 m 5m8 3 + p 6 8m 6 8m m + 60m m 8m 3 m 8m 3 m m + 60m 3 m + p 8 6m 8m m + 5m 8m 3 m + 6m 3 m + 5m 3 + p 0 m + m + 3 m + p, A.5 P 0 p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 m m 3 m 3m + m m 3 + p 8m 6 + m m 0m m 0m6 3m 3 m + m 3 m m

27 E. Remiddi, L. Tancredi / Nuclear Physics B m 3 m + 3m 3 m 0m 3 m 0m6 3 + p 36m m m 8m + 6m 3 m + 66m 3 m 3m 3 + p 6 30m + 0m m p 8, A.6 P p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 m m 3 m m m + m 3 + p 6m 6 m m + 6m m 0m6 + m 3 m + 8m 3 m m 0m 3 m 8m 3 m + 30m 3 m + p m + 8m m 0m 0m 3 m + m 3 m + 6m 3 + p 6 0m 6m + 8m p 8, A.7 P p,m,m,m 3 = m + m + m 3 m m + m 3 m + m m 3 m m m 3 m m + 3 m 5m m m m m 3 m + m 3 m m 3 + p m m + 3 m 7m 6 3m m + 9m m 5m6 9m 3 m + 6m 3 m m + 7m 3 m + 5m 3 m + m 3 m 3m6 3 + p m 6 + m m 6m m + 0m6 m 3 m + 7m 3 m m 6m 3 m + 38m 3 m m 3 m m6 3 + p 6 m 6m m + 0m + 8m 3 m + m 3 m + m 3 + p 8 5m + 5m + 3m 3 + p 0, A.8 P p,m,m,m 3 = 3m m m m m 3 m + m 3 m m 3 + p m m + m 3 + 3p. A.9 The polynomials defined above fulfill, among the others, the relation: P p,m,m,m 3 + P p P p,m,m,m 3 = 0, p,m,m 3,m m P 0 p,m,m,m 3 A.0 where note that the polynomial Pp,m,m,m 3, defined in Eq..3, appears squared. Appendix B. The polynomials of the second-order differential equation In this second appendix we give the explicit expressions of the polynomials that appear in the second-order differential equation derived in Section 5. Also in this case, they are functions of p and of the three masses m, m and m 3, but they do not depend on the dimensions d.