arxiv: v1 [hep-ph] 30 Dec 2015

Transcription

1 June 3, 8 Derivation of functional equations for Feynman integrals from algebraic relations arxiv:5.94v [hep-ph] 3 Dec 5 O.V. Tarasov II. Institut für Theoretische Physik, Universität Hamburg, Luruper Chaussee 49, 76 Hamburg, Germany and Joint Institute for Nuclear Research, 498 Dubna, Russian Federation otarasov@jinr.ru Abstract New methods for obtaining functional equations for Feynman integrals are presented. Application of these methods for finding functional equations for various one- and twoloop integrals described in detail. It is shown that with the aid of functional equations Feynman integrals in general kinematics can be expressed in terms of simpler integrals. PACS numbers:.3.gp,.3.ks,..ds,.38.bx Keywords: Feynman integrals, functional equations

2 Contents Introduction Deriving functional equations from recurrence relations 3 Deriving functional equations from algebraic relations between propagators 4 3. Prototypes of functional equations Some examples of functional equations 8 4. Functional equation for the one-loop propagator type integral Functional equations for the one-loop vertex type integral Functional equations for one-loop box type integrals Functional equation for two-loop vertex type integral Deriving functional equation by deforming propagators 3 5. Algebraic relations for products of deformed propagators Functional equation for two-loop vacuum type integral with arbitrary masses Conclusions 9 7 Acknowledgment 9 Introduction Recently it was discovered that Feynman integrals obey functional equations [], []. Different examples of functional equations were presented in Refs. [], [3],[]. In these articles only one-loop integrals were considered. In the present paper we propose essentially new methods for deriving functional equations. These methods are based on algebraic relations between propagators and they are suitable for deriving functional equations for multi-loop integrals. Also these methods can be used to derive functional equations for integrals with some propagators raised to non-integer powers. Our paper is organized as follows. In Sec.. the method proposed in Ref. [] is shortly reviewed. In Sec. 3. a method for finding algebraic relations between products of propagators is formulated. We describe in detail derivation of explicit relations for products of two, three and four propagators. Also algebraic relation for products of arbitrary number of proparators is given. These relations are used in Sec.4. to obtain functional equations for some one-, as well as two- loop integrals. In particular functional equation for the massless one-loop vertex type integral is presented. Also functional equation for the two-loop vertex type integral with arbitrary masses is given. In Sec. 5. another method for obtaining functional equations is proposed. The method is based on finding algebraic relations for deformed propagators and further conversion of integrals with deformed propagators to usual Feynman integrals by imposing conditions on deformation parameters. To perform such a conversion the α- parametric representation for

3 both types of integrals is exploited. The method was used to derive functional equation for the two-loop vacuum type integral with arbitrary masses. As a by product, from this functional equation we obtained new hypergeometric representation for the one-loop massless vertex integral. In conclusion we formulate our vision of the future applications and developments of the proposed methods. Deriving functional equations from recurrence relations The method for deriving functional equations proposed in Ref. [] is based on the use different kind of recurrence relations. In particular in Refs. [], [], [3], generalized recurrence relations [4] were utilized to obtain functional equations for one-loop Feynman integrals. In general such recurrence relations connect a combination of some number of integrals I,n,..., I k,n corresponding to diagrams, say, with n lines and integrals corresponding to diagrams with fewer number of lines. Diagrams with fewer number of lines can be obtained by contracting some lines in integrals with n lines. Integrals corresponding to such diagrams depend on fewer number of kinematical variables and masses compared to integrals with n lines. Such recurrence relations can be written in the following form: Q j ({m i }, {s q }, ν l, d) I j, n = R k,r ({m i }, {s m }, ν l, d) I k, r, (.) j k,r<n Q j and R k are ratios of polynomials depending on masses m i, scalar products s r of external momenta, powers of propagators ν l and parameter of the space time dimension d. At the left hand-side of Eq. (.) we combined integrals with n lines and on the right hand - side integrals with fewer number of lines. In accordance with the method of Ref. [], to obtain functional equation from Eq. (.) one should eliminate terms on the left hand - side by defining some kinematical variables from the set of equations: Q j ({m i }, {s q }, ν l, d) =. (.) If there is a nontrivial solution of this system and for this solution some R k,r ({m i }, {s m }, ν l, d) are different from zero then the right-hand side of Eq. (.) will represent functional equation. For the one-loop integrals with n propagators I n (d) ({m j}, {s kl }) = d d q iπ d/ n j= P ν j j, (.3) P j = (q p j ) m j + iɛ, (.4) different types of recurrence relations were given in Refs. [4], [5]. Diagram corresponding to this integral is given in Figure. In Refs. [4], [5] the following relation was derived: (d = n i= n k= ν i + )G n I n (d+) ({m j}, {s kl }) n I n (d) ({m j}, {s kl }) ( k n )k n ({m j}, {s kl }), (.5)

4 Figure : One-loop diagram with n external legs the operators k shift index of propagators by one unit ν k ν k, (p p n )(p p n ) (p p n )(p p n )... (p p n )(p n p n ) G n = n (p p n )(p p n ) (p p n )(p p n )... (p p n )(p n p n ) (p p n )(p n p n ) (p p n )(p n p n )... (p n p n )(p n p n ) Y Y... Y n Y Y... Y n n =....,.. Y n Y n... Y nn, (.6) Y ij = m i + m j s ij, s ij = (p i p j ). (.7) Here p i, p j are external momenta going through lines i, j respectively, and m j is mass attributed to j-th line. Gram determinant G n and modified Cayley determinant n are polynomials depending on scalar products and masses. It is assumed that these scalar products are made of d dimensional vectors and G n and n are not subject to any restriction or condition specific to some integer values of d. Eq. (.5) is written in the form corresponding to Eq. (.). To eliminate integrals with n lines on the left hand - side of Eq. (.5) the following conditions to be hold: G n =, n =. (.8) Eq. (.5) is valid for arbitrary kinematical variables and masses. Solution of Eqs. (.8) can be easily done with respect to two kinematical variables or masses. Starting from n = 3 substitution of such solutions into Eq. (.5) gives nontrivial functional equations. The method for obtaining functional equations by eliminating complicated integrals from recurrence relations is quite general one. However for multi loop integrals, depending on several kinematical variables, derivation of equations like Eq. (.5) is computationally challenging. In the next sections we will describe easier and more powerful methods that can be used for deriving functional equations for multi-loop integrals. 3

5 3 Deriving functional equations from algebraic relations between propagators Setting ν j = in Eq. (.5) and imposing conditions (.8) leads to the following equation: n k= ( k n )k n =. (3.9) In Eq. (3.9) integrands of k I n (d) are products of n propagators depending on different external momenta, i.e. each term in this relation corresponds to the same function but with different arguments. In fact functional equations considered in Refs. [, 3, ] are of the same form as Eq. (3.9). The question naturally arises: This relationship holds for integrals or it can be obtained as the consequence of a relationship between integrands? By inspecting Eq. (3.9), one can suggest the following form of the relation between products of propagators of integrands: n r= = P r P n+ n n ( x r r= j= j r P j ), (3.) P j = (k p j ) m j + iɛ. (3.) In what follows we will omit iɛ term assuming that all masses have such a correction. Additionally we assume that vectors p j are linearly dependent, i.e. the Gram determinant for the set of vectors {p j } is equal to zero. Such a condition is valid for all examples considered in Refs. [], []. Now let s consider in detail implementation of our prescription for products of,3 and 4 propagators. At n = relation (3.) reads: P P = x P P 3 + x P P 3, (3.) P = (k p ) m, P = (k p ) m, P 3 = (k p 3 ) m 3. (3.3) According to our assumption three vectors p,p,p 3 are linearly dependent. Without loss of generality we may assume that p 3 = y 3 p + y 3 p. (3.4) Furthermore, we assume that k will be integration momentum and scalar quantities x,x, y 3, y 3 do not depend on k. Putting all terms in Eq. (3.) over a common denominator and then equating to zero the coefficients in front of various products of k, k p,k p yields the following system of equations: y 3 x =, y 3 x =, x + x =, p (x y 3) + p (x y 3) + y 3 y 3 (s p p ) m x m x + m 3 =. (3.5) Solution of this system of equations is: x = y 3 = λ, x = y 3 = λ, (3.6) 4

6 λ is a root of the equation with A λ + B λ + C =, (3.7) A = s, B = m m s, C = m m 3. (3.8) This solution can be rewritten in an explicit form: x = y 3 = m m + s s ± x = y 3 = m m + s s Λ + 4s m 3, s Λ + 4s m 3, (3.9) s Λ = s + m 4 + m 4 s m s m m m. (3.) Now let s find algebraic relation for the products of three propagators. At n = 3 Eq. (3.) reads: = x + x + x 3, (3.) P P P 3 P 4 P P 3 P P 4 P 3 P P P 4 P, P, P 3 are defined in Eq.(3.3) and P 4 = (k p 4 ) m 4. (3.) In complete analogy with the previous case we can represent one momentum as a combination of other ones. Without loss of generality we may write p 4 = y 4 p + y 4 p + y 43 p 3, (3.3) y ij for the time being are arbitrary coefficients. Putting all terms in Eq. (3.) over a common denominator and then equating to zero the coefficients in front of various products of (k ), (k p ), (k p ), (k p 3 ) yields the following system of equations: y 43 x 3 =, y 4 x =, y 4 x =, x + x + x 3 =, p (x y 4) + p (x y 4) + p 3(x 3 y 43) +y 4 y 43 (s 3 p p 3) + y 4 y 43 (s 3 p p 3) + y 4 y 4 (s p p ) m x m x m 3x 3 + m 4 =. (3.4) Solving these equations for x, x, x 3, y 4, y 4 we have x = λ 3 y 43, x = λ 3, x 3 = y 43, y 4 = y 43 λ 3, y 4 = λ 3, (3.5) λ 3 is solution of the equation Here A 3 = s, A 3 λ 3 + B 3 λ 3 + C 3 =. (3.6) B 3 = y 43 (s 3 + s s 3 ) m + m s, C 3 = y 43s 3 + (m 3 m s 3 )y 43 + m m 4. (3.7) 5

7 Let us now turn to the derivation of algebraic relation for the product of four propagators. At n = 4 Eq. (3.) reads: P P P 3 P 4 = x P 5 P P 3 P 4 + x P P 5 P 3 P 4 + P, P, P 3,P 4 are defined in Eqs. (3.3), (3.), and p 5 is a linear combination of vectors p,...,p 4, x 3 P P P 5 P 4 + x 4 P P P 3 P 5, (3.8) P 5 = (k p 5 ) m 5, (3.9) p 5 = y 5 p + y 5 p + y 53 p 3 + y 54 p 4. (3.3) Putting all terms in Eq. (3.8) over a common denominator and then equating to zero the coefficients in front of different products of k, k p j yields system of equations: y 5 x =, y 5 x =, y 53 x 3 =, y 54 x 4 =, x + x + x 3 + x 4 =, m 5 m x m x m 3x 3 m 4x 4 + p 4(x 4 y 54) + p 3(x 3 y 53) + p (x y 5) + p (x y 5) + (s p p )y 5 y 5 + (s 3 p p 3)y 5 y 53 + (s 3 p p 3)y 5 y 53 + (s 4 p p 4)y 54 y 5 + (s 4 p p 4)y 54 y 5 + (s 34 p 3 p 4)y 54 y 53 =. (3.3) Solving this system for x, x, x 3,x 4, y 5, y 54 we have λ 4 is a solution of the equation with A 4 = s 4, x = y 5 = λ 4, x = y 5, x 3 = y 53, x 4 = y 54 = x x x 3 = y 53 y 5 λ 4, (3.3) A 4 λ 4 + B 4 λ 4 + C 4 =, (3.33) B 4 = (s 4 s + s 4 )y 5 + (s 34 s 3 + s 4 )y 53 + m m 4 s 4, C 4 = s 4 y 5 + (s 34 s 3 + s 4 )y 5 y 53 + (m m 4 s 4 )y 5 + s 34 y 53 + (m 3 m 4 s 34 )y 53 + m 4 m 5. (3.34) Eqs. (3.), (3.) and (3.8) will be used in the next sections to derive functional equations for the propagator, vertex and box type of integrals. Relations between products of five and more propagators can be easily derived in the same way as as it was done for products of two-, three- and four- propagators. From Eq. (3.) one can derive system of equations and find its solution for arbitrary n. Multiplying both sides of Eq. (3.) by the product of n + propagators n+ j= P j yields n P n+ = P r, (3.35) 6 r=

8 or k k p n+ + p n+ = n x r (k k p r + p r m r). (3.36) r= Since we assume linear dependence of vectors p r, without loss of generality we may write: n p n+ = y n+,j p j. (3.37) j= Substituting (3.37) into Eq.(3.36), collecting terms in front of k, k p j and terms without k, equating them to zero after some simplifications yields the following system of n + equations: x k y n+,k =, k =,.., n, (3.38) n y n+,k =, (3.39) k= m n+ n y n+,k m k + k= j n y n+,j y n+,l s lj =. (3.4) Solving Eq. (3.39) for one of the y ij an substituting this solution into Eq. (3.4) gives quadratic equation for the remaining y ij. This quadratic equation can be solved with respect to one of the parameters y ij. Thus the solution of the system of equations (3.38), (3.39), (3.4) will depend on n arbitrary parameters y ij and one arbitrary mass m n. It is interesting to note that for any n, functional equations for integrals with all masses equal to zero and functional equations for integrals with all masses equal are the same. In case of equal masses, two mass dependent terms in Eq. (3.4) cancel each other due to Eq. (3.39). In both cases systems of equations for x i, y jk are the same and therefore arguments of integrals are the same. Eq. (3.) is analogous to the equation for splitting propagators presented in Ref. [6]. Eq. (3.) is a generalization of Eq. (3.). Indeed, setting y 43 = x 3 =, canceling common factor P 3 on both sides of Eq. (3.) yields relation similar to (3.). In turn Eq. (3.8) is a generalization of (3.). j= 3. Prototypes of functional equations Multiplying algebraic relations (3.),(3.), (3.8) by products of any number of propagators raised to arbitrary powers ν j N (3.4) [(k j=n p j ) m j ]ν j and integrating with respect to k we get a functional equation for one-loop integrals. Eqs. (3.), (3.), (3.8) also can be used to derive functional equations for integrals with any number of loops. Multiplying algebraic relations for propagators by function corresponding to Feynman integral depending on momentum k and any number of external momenta and then integrating with respect to k will produce functional equations. Just for demonstrational purposes we present graphically in Figure functional equation based on n propagator relation. The blob on this picture correspond to either product of propagators raised to arbitrary powers or to an integral with any number of loops and external legs. One of the external momenta of this multi loop integral should be k. 7 l=

9 Figure : n - propagator functional equation 4 Some examples of functional equations In this section several particular examples of functional equations resulting from algebraic relations for products of propagators will be considered. 4. Functional equation for the one-loop propagator type integral First, we consider the simplest case, namely, functional equation for the integral d (m j, m d k k; s jk ) = iπ d/ [(k p j ) m j ][(k (4.4) p k ) m k ]. Integrating both sides of Eq. (3.) with respect to k, we get: (m, m, s ) = x (m, m 3, s 3 ) + x (m, m 3, s 3 ). (4.43) The arguments s 3, s 3 of integrals on the right hand - side depend on y 3, y 3 s 3 = (p p 3 ) = (y 3 ) p + y 3 (y 3 )p p + y 3p, s 3 = (p p 3 ) = p y 3 + (y 3 )y 3 p p + (y 3 ) p. (4.44) Substituting solution for y ij from Eq. (3.6) into Eq. (4.44) yields: : s 3 = m4 m m m s + m 4 + s + s m 3 s ± m m + s s Λ + 4s m 3 s 3 = s m s + m 4 + m 4 m m + s m 3 s ± m m s s Λ + 4s m 3. (4.45) In this equation m 3 is an arbitrary parameter and can be taken at will. Functional equation (4.43) is in agreement with the result presented in Refs. [],[]. 8

10 4. Functional equations for the one-loop vertex type integral Functional equations for the vertex type integral 3 (m, m, m 3, s 3, s 3, s ) d d k = iπ d/ [(k p ) m ][(k p ) m ][(k p 3 ) m 3], (4.46) can be obtained from Eq. (3.) as well as from Eq. (3.). Multiplying Eq. (3.) with the factor /P 4 P 4 = (k p 4 ) m 4, (4.47) and integrating over k leads to the equation: d d k d d k = x + x P P P 4 P P 3 P 4 This equation in terms of integrals 3 reads d d k P P 3 P 4. (4.48) 3 (m, m, m 4; s 4, s 4, s ) = x 3 (m, m 3, m 4; s 34, s 4, s 3 ) + x 3 (m, m 3, m 4; s 34, s 4, s 3 ). (4.49) Two more functional equations can be obtained from Eq. (4.48) by symmetric permutations (m, p ) (m 4, p 4 ) and (m, p ) (m 4, p 4 ). Another functional equation for the vertex type integral can be obtained by integrating Eq. (3.) with respect to k : 3 (m, m, m 3, s 3, s 3, s ) = ( λ 3 y 43 ) 3 (m 4, m, m 3, s 3, s 34, s 4 ) + λ 3 3 (m, m 4, m 3, s 34, s 3, s 4 ) + y 43 3 (m, m, m 4, s 4, s 4, s 3 ), (4.5) s 4 = (m m + s )λ 3 + (m m 3 + s 3 )y 43 + m 4 m, s 4 = (m s m )λ 3 + (m m 3 s + s 3 )y 43 + m 4 m + s, s 34 = (m s 3 + s 3 m )λ 3 + (m m 3 s 3 )y 43 + m 4 m + s 3. (4.5) There is an essential difference between functional equation Eq. (4.49) obtained from Eq. (3.) and functional equation (4.5) derived from Eq. (3.). For example, at m = m = m 3 = m 4 =, Eq. (4.49) becomes trivial while from Eq. (4.5) for the integral we obtain nontrivial functional equation: 3 (,,, s 3, s 3, s ) 3 (s 3, s 3, s ) (4.5) 3 (s 3, s 3, s ) = ( λ 3 y 43 ) 3 (s ( λ 3 y 43 ) + s 3 y 43, s 3, λ 3 (s 3 s 3 ) + s 3 ( y 43 )) +λ 3 3 (s 3 y 43 + λ 3 s, λ 3 (s 3 s 3 ) + s 3 ( y 43 ), s 3 ) +y 43 3 (s, s ( λ 3 y 43 ) + s 3 y 43, s 3 y 43 + λ 3 s ), (4.53) 9

11 λ 3 is a root of the quadratic equation s λ 3 (s s y 43 s 3 y 43 + s 3 y 43 )λ 3 + y 43s 3 y 43 s 3 =. (4.54) If one argument of 3 (s 3, s 3, s ) is zero then by applying functional equation (4.53) such an integral can be expressed in terms of integrals 3 with two arguments equal to zero. For example, at s 3 = and y 43 = s /(s s 3 ) the relation (4.53) becomes: 3 (, s 3, s ) = s s s 3 3 (s,, ) s 3 s s 3 3 (,, s 3 ). (4.55) This is a typical example how functional equations can be used to simplify evaluation of an integral by reducing it to a combination of integrals with fewer number of arguments. At m = m = m 3 = m 4 = m, similar to the previous case, Eq.(4.49) degenerate while from Eq.(4.5) for the integral we obtain nontrivial functional equation: 3 (m, m, m, s 3, s 3, s ) 3 (m ; s 3, s 3, s ) (4.56) 3 (m ; s 3, s 3, s ) = ( λ 3 y 43 ) 3 (m ; s ( λ 3 y 43 ) + s 3 y 43, s 3, λ 3 (s 3 s 3 ) + s 3 ( y 43 )) +λ 3 3 (m ; s 3 y 43 + λ 3 s, λ 3 (s 3 s 3 ) + s 3 ( y 43 ), s 3 ) +y 43 3 (m ; s, s ( λ 3 y 43 ) + s 3 y 43, s 3 y 43 + λ 3 s ), (4.57) λ 3 is a root of the quadratic equation s λ 3 (s s y 43 s 3 y 43 + s 3 y 43 )λ 3 + y 43s 3 y 43 s 3 =. (4.58) Eqs. (4.57), (4.58) are identical to Eqs. (4.53),(4.54) respectively and therefore functional equation for the integral with massless propagators and functional equation for the integral with all masses equal are the same. Eq. (4.57) at s 3 = and y 43 = s /(s s 3 ) leads to the relation similar to (4.55): 3 (m ;, s 3, s ) = s s s 3 3 (m ; s,, ) s 3 s s 3 3 (m ;,, s 3 ). (4.59) This is not surprising because coefficients of the Eq. (4.57) are mass independent and in the integrand m and iɛ appear in the covariant combination m iɛ. For this reason the similarity of functional equations for massless integrals and integrals with all masses equal take place for integrals with more external legs and more loops. 4.3 Functional equations for one-loop box type integrals Functional equations for the box type integrals can be obtained by multiplying relation (3.) by two propagators, or by multiplying relation (3.) by one propagator and then integrating over momentum k. Yet another relation can be obtained just by integrating Eq. (3.8) over

12 momentum k : 4 (m, m, m 3, m 4; s, s 3, s 34, s 4, s 4, s 3 ) = λ 4 4 (m 5, m, m 3, m 4; s 5, s 3, s 34, s 45, s 4, s 35 ) + y 5 4 (m, m 5, m 3, m 4; s 5, s 35, s 34, s 4, s 45, s 3 ) + y 53 4 (m, m, m 5, m 4; s, s 5, s 45, s 4, s 4, s 5 ) +( y 5 y 53 λ 4 ) 4 (m, m, m 3, m 5; s, s 3, s 35, s 5, s 5, s 3 ). (4.6) Here λ 4 is defined in Eq. (3.33) and m 5,y 5, y 53 are arbitrary parameters and s 5 = λ 4 (m 4 s 4 m ) + y 53 (m 4 m 3 s 4 + s 3 ) + y 5 (m 4 m s 4 + s ) + s 4 + m 5 m 4, s 5 = λ 4 (m 4 m s 4 + s ) + y 53 (m 4 m 3 s 4 + s 3 ) + y 5 (m 4 m s 4 ) + s 4 + m 5 m 4, s 35 = λ 4 (m 4 m s 34 + s 3 ) + y 53 (m 4 m 3 s 34 ) + y 5 (m 4 m s 34 + s 3 ) + s 34 + m 5 m 4, s 45 = λ 4 (s 4 + m 4 m ) + y 53 (s 34 + m 4 m 3) + y 5 (s 4 + m 4 m ) + m 5 m 4. (4.6) Arbitrary parameters in this functional equation can be chosen from the requirement of simplicity of evaluation of integrals on the right hand - side of Eq. (4.6) or from some other requirements. For example, one can choose these parameters by transforming arguments to a certain kinematical region needed for analytic continuation of the original integral. 4.4 Functional equation for two-loop vertex type integral The method described in the previous section can be applied to multi loop integrals. Consider, for example, integral corresponding to the diagram given in Figure 3. If we multiply Eq. (3.) Figure 3: Diagram corresponding to R(m, m, m 3, m 4; q, q, q 3) by the one-loop integral depending on k d d k [k m 4][(k k ) m 5] and integrate with respect to momentum k then we obtain functional equation R(m, m, m 3, m 4; q, q, q 3) = αr(, m, m 3, m 4; Q, q, (m m + q 3)α m ) (4.6) + ( α)r(m,, m 3, m 4; q, Q, (m m q 3)α + q 3 m ), (4.63)

13 R(m, m, m 3, m 4; q, q, q3) d d k d d k = (iπ d/ ) [(k + q ) m ][(k q ) m ][k m 3][(k k ) m 4], (4.64) Q = (q q m + m )α + q m, (4.65) α = q 3 m + m ±, q3 (4.66) = q3 4 + m 4 + m 4 q3m q3m m m. (4.67) Integrals of this type arise, for example, in calculations of two-loop radiative corrections in the electroweak theory. Instead of the integral R(m, m, m 3, m 4; q, q, q 3) one can consider derivative of R with respect to m 3 which is UV finite: R R m 3 (m, m, m 3, m 4; q, q, q3) 3 d d k d d k = (iπ d/ ) [(k + q ) m ][(k q ) m ][k m 3] [(k k ) m 4].(4.68) Integral R 3 satisfy the following functional equations: R 3 (m, m, m 3, m 4; q, q, q 3) = αr 3 (, m, m 3, m 4; Q, q, (m m + q 3)α m ) + ( α)r 3 (m,, m 3, m 4; q, Q, (m m q 3)α + q 3 m ). (4.69) This relation can be used for computing basis integral arising in calculation of two-loop radiative correction to the ortho -positronium lifetime. In particular one of these basis integrals corresponds to kinematics m = m = m 3 = m 4 = m, q = q = m, q 3 = 4m. In this case relation (4.69) reads R 3 (m, m, m, m ; m, m, 4m ) = R 3 (, m, m, m ;, m, m ). (4.7) Integral on the right hand-side is in fact propagator type integral with one massless line. Applying recurrence relations given in Ref. [7] this integral can be reduced to simpler integral: R 3 (, m, m, m ;, m, m ) d d k d d k = (iπ d/ ) k(k m )[(k k ) m ] [(k + q ) m ] = 3(d 3) J (d ) (m ), (4.7) J (d) (q ) = (iπ d/ ) At q = m, the result for J (d) (q ) is known [8]: m 6 d 8 d d k d d k (k m )[(k k ) m ][(k q) m ]. (4.7) (d 4) (d )(d 3)J (m (d) ) [, 3 d, = 3 F ; ] [, 4 d 5 d, d ; + (d 3) 3 F, ] d ; 3, d ;, (4.73)

14 and it can be used for the ε expansion of R and R 3. As was already mentioned at q = q = m, q 3 = 4m integrals on the right hand - side of Eq.(4.63) correspond to propagator type integrals. Analytic result for R reads (d 3)R(m, m, m, m, m, m, 4m ) [ ( = m (3d 8)(d 3)J (d) 8m 4 3(d ) (d 4) T (d) ) ], (4.74) ( T (d) = Γ d ) (m ) (d )/. (4.75) We checked that several first terms in the ε = (4 d)/ expansion of R and R 3 are in agreement with results of [9]. The main profit from functional equations for R and R 3 comes from the fact that vertex integrals were expressed in terms of simpler, propagator type integrals. 5 Deriving functional equation by deforming propagators The method described in the previous section does not work for deriving functional equations for all kinds of Feynman integrals. For example, we did not found functional equation for the two-loop vacuum type integral given in Figure 4. Figure 4: Two-loop vacuum type diagram In this section we shall describe another method that extends the class of integrals for which we can obtain functional equations. The method is based on transformation of functional equations for some auxiliary integrals depending on arbitrary parameters into functional equations for integrals of interest. Such functional equations will be derived from algebraic relations for deformed propagators which will be defined in the next section. These auxiliary integrals will be transformed into α parametric representation. In general characteristic polynomials of these integrals in α parametric representation differ from those for the investigated integral. Functional equation for the integral of interest can be obtained in case when it will be possible to map characteristic polynomials of auxiliary integrals with deformed propagators to characteristic polynomials of this integral. Such a mapping will be performed by rescaling α parameters and appropriate choice of arbitrary deforming parameters. 5. Algebraic relations for products of deformed propagators In the previous section to derive functional equation we added to our consideration a propagator with combination of external momenta taken with arbitrary scalar coefficient. Now we consider 3

15 generalization of this method. To find functional equation for L-loop Feynman integral depending on E- external momenta we start from the relation of the form D j is defined as: with n r= = D r D n+ Q j = n n ( x r r= j= j r D j ), (5.76) D j = Q j m j + iɛ, (5.77) L l= E+ a (j) jl k l + b (j) jl p l, (5.78) and a (j) jl, b(j) jl for the time being are arbitrary scalar parameters. Some of these parameters as well as x r will be fixed from the equation (5.76). Another part of these parameters will be fixed from the requirement that the product of propagators in (5.76) should correspond to the integrand of the integral with the considered topology. We would like to remark that instead of deformation of propagators proposed in Eqs. (5.77),(5.78) one can use other deformations. For example, all terms in denominators of propagators can be taken with arbitrary scalar coefficients: Q j = L l,r= a (j) rl (k r k l ) + L E+ l= r= l= E+ b (j) rl (p r k l ) + r,l= c (j) rl (p r p l ). (5.79) To establish algebraic relation (5.76) we put all terms over a common denominator and then equate coefficients in front scalar products depending on integration momenta. Solving obtained system of equations gives some restrictions on the scalar parameters. In general integrals obtained by integrating products of deformed propagators will not correspond to usual Feynman integrals. Further restrictions on parameters should be imposed in order to obtain relations between integrals corresponding to Feynman integrals coming from a realistic quantum field theory models. 5. Functional equation for two-loop vacuum type integral with arbitrary masses As an example, let us consider derivation of functional equation for the two-loop vacuum type integral given in Figure 4: J (d) d (m, m, m d k d d k 3) = (iπ d/ ) (k m )((k k ) m )(k m 3). (5.8) Analytic expression for this integral was presented in Ref []. Instead of this integral we will first consider an auxiliary integral with integrand made from deformed propagators defined in Eqs.(5.77), (5.78): J (d) d (m, m, m d k d d k 3) =, (5.8) (iπ d/ ) D D D 3 4

16 D = (a k + a k ) m + iɛ, D = (b k + b k ) m + iɛ, D 3 = (h k + h k ) m 3 + iɛ. (5.8) For the product of three deformed propagators one can try to find an algebraic relation of the form: x x x 3 = + +, (5.83) D D D 3 D 4 D D 3 D D 4 D 3 D D D 4 D, D,D 3 are defined in Eq.(5.8) and D 4 = (r k + r k ) m 4. (5.84) Here m k are arbitrary masses, x k, a j, b i, h s, r l are undetermined parameters and k, k will be integration momenta. Putting in Eq. (5.83) all over a common denominator and equating to zero coefficients in front of different products of k, k and k k leads to the system of equations: r x a x b x 3 h =, r r x a a x b b x 3 h h =, r x a x b x 3 h =, m 4 x m x m x 3 m 3 =. (5.85) Solving this system for r, x, x, x 3 we have: r = r λ, Ax = r (h h m b b m 3) + b h (b h b h )m 4 r (h m m 3b ) λ, Ax = r (m h h a a m 3) a h (a h a h )m 4 + r (h m m 3a ) λ, Ax 3 = r (m b b a a m ) + a b (a b a b )m 4 r (b m m a ) λ, (5.86) λ is a root of the quadratic equation and Aλ + Bλ + C =, (5.87) A = b h (b h h b )m + a h (h a a h )m + a b (a b b a )m 3, B = (h b b h )(b h + h b )m + (a h h a )(a h + h a )m + (b a a b )(a b + b a )m 3, C = b h (b h h b )m + a h (a h a h )m + a b (a b b a )m 3 + (b h b h )(a h a h )(a b a b )m 4. (5.88) In order to obtain functional equation for the integral J (d), we integrate first both sides of the Eq. (5.83) with respect to k, k and then convert these integrals into the α-parametric representation. Transforming all propagators into a parametric form (k m + iɛ) ν = i ν Γ(ν) dα α ν exp [ iα(k m + iɛ) ], (5.89) and using the d- dimensional Gaussian integration formula d d k exp [ i(ak + (pk)) ] ( ] π = i exp [ ia)d ip, (5.9) a 5

17 we can easily evaluate the integrals over loop momenta. The final result is: 3 [ ] J (d) (m, m, m 3)= i d dα j 3... [ D(α)] exp i α d l (m l iɛ), (5.9) the polynomial j= D = (a b a b ) α α + (a h a h ) α α 3 + (b h b h ) α α 3, (5.9) differs from the appropriate D polynomial of the two-loop vacuum integral defined in Eq.(5.8) Changing in Eq. (5.9) integration variables with leads to the relation: D = α α + α α 3 + α α 3. (5.93) l= α j α j θ j (5.94) θ = b h b h, θ = a h a h, θ 3 = a b a b, (5.95) D = (b h b h ) (a h a h ) (a b a b ) (α α + α α 3 + α α 3 ). (5.96) Therefore integral J (d) with deformed propagators is proportional to J (d) with modified arguments: d d k d d k = i d [θ (iπ d/ ) D D D θ θ 3] d 3 dα dα dα 3 exp [ im] = [θ θ θ 3] d J (d) (θ m, θ m, θ 3m 3) (5.97) M = α θ m + α θ m + α 3 θ 3m 3. (5.98) With the aid of Eq. (5.97) relation (5.83) integrated with respect to k, k can be written as a combination of integrals J (d) with different arguments: θ,θ,θ 3 are defined in (5.95) and [θ θ θ 3] d J (d) (θ m, θ m, θ 3m 3) D d = x [θθ 6θ 4] d (d) J (θm 4, θ6m, θ4m 3) + x [θθ 5θ 6] d (d) J (θ6m, θm 4, θ5m 3) + x 3 [θ 3θ 4θ 5] d J (d) (θ 4m, θ 5m, θ 3m 4), (5.99) θ 4 = r b r b, θ 5 = r a r a, θ 6 = r h r h. (5.) Now we consider relation (5.99) as an equation for integrals in momentum representation. By rescaling integration variables k, k in the integral on the left hand - side k = (θ θ θ 3 ) k, k = (θ θ θ 3 ) k, (5.) 6

18 and performing analogous changes for the integrals on the right hand - side we obtain the relation ( ) J (d) θ m θ θ θ θ 3 θ θ, m θ 3 3 θ θ, m 3 3 θ θ = x ( ) J (d) θ m θ 6 θ θ 4 θ 6 θ 4 θ 4, m θ 4 6 θ θ, m 3 4 θ θ 6 + x J (d) θ θ 5 θ 6 + x 3 J (d) θ 3 θ 4 θ 5 ( θ6 m θ θ, ( 5 θ4 m θ 3 θ, 5 θ θ 6 θ 5 m 4, θ 5 θ 3 θ 4 m, θ 5 ) m 3 θ θ 6 ) θ 3 m 4. (5.) θ 5 θ 4 In terms of redefined masses M, M, M 3 related to original masses m, m, m 3 as equation (5.) reads m = θ θ 3 θ M, m = θ θ 3 θ M, m 3 = θ θ θ 3 M 3, (5.3) J (d) (M, M, M3 ) = θ θ 3 x J (d) θ 4 θ6 + θ θ 3 x J (d) θ 5 θ6 + θ θ x 3 J (d) θ 4 θ5 ( θ 3 θ6 ( θ m 4, θ 3 θ 6 M, θ θ 4 M3 θ 4 θ6 θ θ4 θ 3 θ6 ) θ θ 5 θ6 m 4, θ θ 5 M, θ θ5 ( θ θ4 M, θ θ 5 M, θ θ5 θ θ4 M3 θ 3 θ6 ) θ 3 m 4 θ 4 θ5 ), (5.4) x i, θ j are defined from Eqs. (5.86),(5.88) with redefined masses. After simplifications Eq. (5.4) takes a simpler form: J (d) (M, M, M3 ) = J (d) ( ) ρ + ρ + M M M3, ρ, ρ ( ) ρ3, M3 ρ, ρ 3 ρ M + M +J (d) J (d) ( ) ρ3 M, ρ M, ρ + ρ 3 M3, (5.5) with ρ = (M 3 + M M )θ θ θ 3 b h θ θ 3 m 4 δ θ θ θ 3 M b θ m 4 ρ = (M 3 + M M )θ θ θ 3 b h θ θ 3 m 4 + δ θ θ θ 3 M 3 h θ 3m 4 ρ 3 = (M M + M 3 )θ θ θ 3 + a h θ θ 3 m 4 + δ θ θ θ 3 M a θ m 4 M M 3 M, (5.6) δ = θ θ θ 3 + 4Aθ θ θ 3 m 4, = M 4 + M 4 + M 4 3 M M M M 3 M M 3. (5.7) 7

19 The coefficinets A, B and C can be expressed in terms of M i, θ j : A = b h θ θ 3 M a h θ θ 3 M + a b θ θ M 3, B = (b h + b h )θ θ 3 M + (a h + a h )θ θ 3 M (a b + a b )θ θ M 3, C = b h θ θ 3 M θ θ 3 a h M + a b θ θ M 3 θ θ θ 3 m 4. (5.8) One can easily observe that due to relations ρ ρ = M M 3 3, ρ ρ 3 = M ρ + M ρ 3, ρ ρ 3 = M 3 ρ 3 M M 3, (5.9) parameters ρ and ρ 3 can be expressed in terms of ρ and therefore all the auxiliary parameters introduced in derivation of the functional equation will be absorbed only in one parameter - ρ. We would like to notice that Eq. (5.5) is valid for integrals but not for their integrands. This is due to the fact that the factor in front of integral that comes from the scaling of α parameters in parametric integral is not fully compensated by scaling momenta given in Eq. (5.). At m 4 = the dependence on all parameters a i,b j,h k in Eqs. (5.4), (5.5) drops out and the integral J (d) (M, M, M 3 ) reduces to a comination of simpler integrals: J (d) (M, M, M3 ) ( = J (d), M + M + M3 +, M + M + M3 ) ( +J (d) M M + M3 +,, M M + M3 ) J (d) ( M M + M3 +, M M + M3, ). (5.) Analytic expression for the integral J (d) with one mass equal to zero is known []. Under assumption that m 3 m m 3 it reads (, m, m 3) = π ( d) (m m 3) d J (d) (d 3)m 3 ( Γ ( d ) sin πd ) F [, d ; 4 d ; ] m. (5.) m 3 From functional equation (5.) as a by-product one can get a new hypergeometric representation for the one-loop massless vertex type integral. In Ref. [] an interesting relation between the dimensionally regularized one-loop vertex type integral 3 (m, m, m 3) and the two-dimensional integral J (6 d) (m, m, m 3) was discovered ) 3 (s 3, s 3, s ) = Γ ( 3 d Γ(d 3) ( s 3s 3 s ) d 4 J (6 d) (s 3, s 3, s ). (5.) Functional equation (5.) with J (d) (, m, m 3) defined in Eq. (5.) provide us a new hypergeometric representation for the integral 3 with massless propagators. Formula for the one-loop massless vertex integral in terms of other Gauss hypergeometric functions is given in Ref. []. 8

20 6 Conclusions Finally, we summarize what we have accomplished in this paper. First of all, we formulated new methods for deriving functional equations for Feynman integrals. These methods are rather simple and do not use any kind of integration by parts techniques. Second, it was shown that integrals with many kinematic arguments can be reduced to a combination of simpler integrals with fewer arguments. In our future publications we are going to demonstrate that in some cases applying functional equations one can reduce, the so-called, master integrals to a combination of simpler integrals from, what we would like to call, a universal basis of integrals. The method based on algebraic relations for deformed propagators can be used not only for vacuum type of integrals but also for integrals depending on external momenta. In the present paper we considered rather particular cases of functional equations. The systematic investigation and classification of the proposed functional equations requires application of the methods of algebraic geometry and group theory. At the present moment it is not quite clear whether functional equations derivable from recurrence relations can be reproduced by the methods of algebraic relations between products of propagators described in Section 3 and Section 5. A detailed consideration of our functional equations and their application to the one-loop integrals with four, five and six external legs as well as to some two- and three- loop Feynman integrals will be presented in future publications. 7 Acknowledgment This work was supported by the German Science Foundation (DFG) within the Collaborative Research Center 676 Particle, Strings and the Early Universe: the Structure of Matter and Space-Time. I am thankful to O.L. Veretin for providing results for integrals contributing to ortho-positronium lifetime described in Ref.[9]. References [] O. V. Tarasov. New relationships between Feynman integrals. Phys.Lett., B67:67 7, 8. [] O. V. Tarasov. Functional equations for Feynman integrals. Phys.Part.Nucl.Lett., 8:49 47,. [3] B. A. Kniehl and O.V. Tarasov. Functional equations for one-loop master integrals for heavy-quark production and Bhabha scattering. Nucl.Phys., B8:78 9, 9. [4] O. V. Tarasov. Connection between Feynman integrals having different values of the space-time dimension. Phys.Rev., D54: , 996. [5] J. Fleischer, F. Jegerlehner, and O.V. Tarasov. Algebraic reduction of one loop Feynman graph amplitudes. Nucl.Phys., B566:43 44,. 9

21 [6] G. t Hooft and M. J. G. Veltman. Scalar One Loop Integrals. Nucl. Phys., B53:365 4, 979. [7] O. V. Tarasov. Generalized recurrence relations for two loop propagator integrals with arbitrary masses. Nucl.Phys., B5:455 48, 997. [8] D. J. Broadhurst, J. Fleischer, and O. V. Tarasov. Two loop two point functions with masses: Asymptotic expansions and Taylor series, in any dimension. Z. Phys., C6:87 3, 993. [9] B. A. Kniehl, A. V. Kotikov, and O. L. Veretin. Irrational constants in positronium decays. In Proceedings, 9th DESY Workshop on Elementary Particle Theory: Loops and Legs in Quantum Field Theory, 8. [] A. I. Davydychev and J.B. Tausk. Two loop selfenergy diagrams with different masses and the momentum expansion. Nucl.Phys., B397:3 4, 993. [] A. I. Davydychev and J.B. Tausk. A Magic connection between massive and massless diagrams. Phys.Rev., D53: , 996. [] A. I. Davydychev. Explicit results for all orders of the epsilon expansion of certain massive and massless diagrams. Phys.Rev., D6:877,.