arxiv:hep-th/ v1 21 Dec 2006

Transcription

1 Preprint typeset in JHEP style - HYPER VERSION hep-th/64 BU-HEPP-6- arxiv:hep-th/64v Dec 6 All order ε-expansion of Gauss hypergeometric functions with integer and half/integer values of parameters M. Yu. Kalmykov Department of Physics, Baylor University, One Bear Place, Box 9736, Waco, TX Bogoliubov Laboratory of Theoretical Physics, Joint Institute for Nuclear Research, 498 Dubna Moscow Region, Russia kalmykov@theor.jinr.ru B.F.L. Ward, S. Yost, Department of Physics, Baylor University, One Bear Place, Box 9736 Waco, TX Abstract: It is proved that the Laurent expansion of the following Gauss hypergeometric functions, F I + aε,i + bε;i 3 + cε;, F I + aε,i + bε;i cε;, F I + + aε,i + bε;i 3 + cε;, F I + + aε,i + bε;i cε;, F I + + aε,i + + bε;i cε;, where I,I,I 3 are an arbitrary integer nonnegative numbers, a,b,c are an arbitrary numbers and ε is an arbitrary small parameters, are expressible in terms of the harmonic polylogarithms of Remiddi and Vermaseren with polynomial coefficients. An efficient algorithm for the calculation of the higher-order coefficients of Laurent expansion is constructed. Some particular cases of Gauss hypergeometric functions are also discussed. Keywords: Gauss hypergeometric functions, harmonic polylogarithms, colour polylogarithms, Laurent expansion of Gauss hypergeometric function, multiloop calculations. Supported by NATO Grant PST.CLG.9834 and DOE grant DE-FG-5ER4399

2 Contents. Introduction. All-order ε-expansion 3. Non-ero values of the ε-dependent part 3.. Integer values of of ε-independent parameters 4.. Half-integer values of of ε-independent parameters 8. Zero-values of ε-dependent part of upper parameters 3. Some particular cases 3. The generalied log-sine functions and their generaliation 3. Special cases: all-order ε-expansion in terms of Nielsen polylogarithms 3 4. Conclusion 4. Introduction One of the most powerful techniques for calculating Feynman diagrams is based on their presentation in terms of hypergeometric functions. We will call this the hypergeometric function representation of Feynman diagrams. Such a representation can be used for numerical evaluation, construction of the asymptotic expansion, etc. One of the unsolved problems in this program is obtaining the proper representation for a diagram with an arbitrary number of legs and loops. Direct use of α- or Feynman parameters representations is not very helpful in solving this problem. The Mellin-Barnes technique is restricted to several topologies, 3, 4, 5. The negative dimension approach 6 has a similar restriction 7. The most investigated diagrams are the master integrals typically, integrals with the power of each propagator equal to unity. The differential 8 and/or difference equation 9 techniques are usually used to obtain such representations. The known cases include the one-loop diagrams,, two-loop propagator-type diagrams with special mass and momentum values, several three-loop bubble-type diagrams 3, three-loop vertex-type diagrams 3, and four-loop bubble-type diagrams 4. For practical application however, it is necessary to construct the of ε-expansion Laurent expansion of hypergeometric functions. There is some evidence that the multiple polylogarithms 5, 6, 7, Li k,k,,k n,,, n = m >m > m n> m m n mn m k mk mkn n.. are sufficient for the parametriing of the coefficients of the ε-expansion of any hypergeometric function 8.

3 In some particular cases, the result of the Laurent expansion can be written in terms of simpler functions. In particular, at the present moment, it is commonly accepted 9 that the generalied hypergeometric functions with an arbitrary set of integer parameters can be presented in terms of harmonic polylogarithms. The idea of the proof is based on the properties of nested sums : the analytical coefficients of the ε-expansion of any generalied hypergeometric function with integer parameters can be reduced to a set of some basic harmonic series of the type j= j j c S a j... S ap j,. where is an arbitrary argument and S a j is an harmonic sum defined as S a j = j k= k a. Series of this type are expressible in terms of the Remiddi-Vermaseren harmonic polylogarithms. However for hypergeometric functions with half-integer values of parameters, the new type of sums, multiple inverse binomial sums, 4, 5, 3 are generated: Σ k a,...,a p; b,...,b q;c j= j j j k j c S a j... S ap j S b j... S bq j. For particular values of k, the sums.3 are called generalied harmonic k = inverse binomial binomial sums. At the present moment, there is no proof that any multiple inverse binomial sums can be expressed in terms of harmonic polylogarithms only. This problem was investigated in Ref. 3 for multiple inverse binomial sums up to weight 4. In Ref. 4, it was shown that some of the multiple inverse binomial sums are not expressible in terms of harmonic polylogarithms of simple argument. In Ref. 5, the results of Ref. 3 were extended on the case of special combinations of multiple binomial sums and multiple generalied harmonic sums. However, the Laurent expansion of a hypergeometric function in general contains combinations of multiple sums. These combinations may be expressed in terms of harmonic polylogarithms. From this point of view, the construction of the analytical coefficients of the ε-expansion of hypergeometric functions can be done independently from existing results for each individual multiple sum. The simplest hypergeometric function is the Gauss hypergeometric function F a,b;c;, 6, 7, 8. It satisfies the second-order differential equation d d dd + c w = d d + a We are indebted to A. Davydychev for discussion on this subject.3 d d + b w,.4

4 and admits the series representation F a,b;c; = k= a k b k c k where a k = Γa + k/γa is the Pochhammer symbol. The primary aim of this paper is to prove the following: Theorem The all-order ε-expansions of the Gauss hypergeometric functions k k!,.5 F I + aε,i + bε;i 3 + cε;, F I + aε,i + bε;i cε;, F I + + aε,i + bε;i 3 + cε;, F I + + aε,i + bε;i cε;, F I + + aε,i + + bε;i cε;,.6a.6b.6c.6d.6e where {I k } are integer numbers, a,b,c are an arbitrary numbers, and ε is an arbitrary small parameter, are expressible in terms of Remiddi-Vermaseren harmonic polylogarithms with rational coefficients.. All-order ε-expansion. Non-ero values of the ε-dependent part It is well known that any Gauss hypergeometric function may be expressed as a linear combination of two other hypergeometric functions with parameters differing from the original ones by an integer 6, 7, 8, 9, 5. Such a representation will be called a reduction, and the explicit algorithm will be called a reduction algorithm. Using the algorithm described in 5, the result of the reduction can be written as { } Pa,b,c, F a + I,b + I ;c + I 3 ; = Q a,b,c, d d + Q a,b,c, F a,b;c;,. where a,b,c, are any fixed numbers, P,Q,Q are polynomial in parameters a,b,c and argument, and I,I,I 3 any integer numbers. All the hypergeometric functions.6 listed in Theorem can be reduced to functions with I,I,I 3 equal to ero for half-integer values of parameters, and to unity for integer ones. In this way, all the hypergeometric functions of Theorem I are expressible in terms of the following five basic functions and their first derivatives: F a ε,a ε; + cε;, F a ε,a ε; + fε;,.a F +bε,aε;+cε;, F +bε,aε; +fε;, F +b ε, +b ε; +fε;..b 3

5 It was shown in Ref. 5 that only the two hypergeometric functions and their first derivative of type.a are algebraically independent. The other three,.b, are algebraically expressible in terms of F a ε,a ε; + fε;. Consequently, in order to prove Theorem I, it sufficient to show that the analytical coefficients of the ε-expansion of the first two hypergeometric functions.a are expressible in terms of Remiddi-Vermaseren polylogarithms... Integer values of of ε-independent parameters Let us start from expansion of Gauss hypergeometric functions with integer values of parameters, and consider the function F a ε,a ε; + cε;. In ref. 3, the all-order ε-expansions for this functions and its first derivative were constructed in terms of multiple polylogarithms of one variable 6, 7. These multiple polylogarithms may be expressed as iterated integrals and have the expansion Li k,k,,k n = m >m > m n> m m k mk mkn n..6 Similar results were also derived without explicit form of coefficients in Ref. 9 via nested sums approach. We will follow the idea of Ref. 3. The Gauss hypergeometric function F a ε,a ε;+cε; is the solution of the differential equation d dd d + cε w = d d + a ε d d + a ε w,.7 with boundary conditions w = and d d w = =. Eq..7 is valid in each order of ε, so that in terms of coefficients functions w k defined as w = w k ε k,.8 Recall that multiple polylogarithms can be expressed as iterated integrals of the form where, by definition Z Li k,,k n = Z dt t dt t dt { t } k times k= dt t dt t dt { t } k times dt t dt t dt t dt { t } k n times dt Z Z t = dt t Z dt tk Z dt tk k t t t k dt t,.3 dt k t k..4 The integral.3 is an iterated Chen integral 3 see also 3 w.r.t. the two differential forms ω = d/ and ω = d, so that Li k,,k n = Z ω k ω ω kn ω..5 4

6 it can be written d d d w k = a + a c d w k + a a w k.9 d d for k with w =,.a w k =, k <..b The boundary conditions for the coefficient functions are w k =, k,.a d d w k = =, k,.b Let us introduce a new function ρ defined by 3 ρ = d d w = ρ k ε k,. where the coefficient functions satisfy k= ρ k = d d w k..3 The boundary conditions for these new functions follow from Eq..: ρ k =, k..4 Eq..9 can be rewritten as a system of two first-order differential equations: d d ρ i = a +a c ρ i +a a w i, d d w i = ρ i..5 The solution of this system can be presented in an iterated form: dt ρ i = a +a c t ρ dt i t+a a t w i t cw i, i, dt w i = t ρ it, i..6 3 We may note that F + a ε,+a ε +cε! = + cε X» ρk+ ε k. a a k= 5

7 Taking into account that w = and ρ = the ε-expansion of ρ begins with the term linear in ε, we obtain the first few coefficients, ρ = w =,.7a ρ = ln H;, a a.7b w = Li H,;, a a.7c ρ 3 = γ c a a ln cli γ c H,; ch,;,.7d w 3 = γ c S, cli 3 γ c H,,; ch,,;, a a.7e where we have defined γ c = a +a c, and Li n and S a,b are the classical and Nielsen polylogarithms 33, 34, respectively: S a,b = a+b dξ lna ξ ln b ξ, S a, = Li a+. a! b! ξ The functions H A; are the Remiddi-Vermaseren harmonic polylogarithms, and A is a multiple index including only entries and. From the representation.6 and result for the first few coefficients.7 we may derive the following observations: Corollary : The all-order ε-expansion of the function F a ε,a ε; + cε; may be written in terms of harmonic polylogarithms H A only, where the multiple index A includes only the values and. Corollary : The analytical coefficient of ε k in the expansion of F a ε,a ε; + cε; includes only functions of weight k with numerical coefficients. Corollary 3: The non-constant terms of the ε-expansion of F a ε,a ε; + cε; are proportional to the product a a in any order of ε. The first and the last statement follows from the representation 4.6, the explicit value of coefficients functions w k, k =,, Eqs..7, and definition of harmonic polylogarithms. The second statement follows from the form of the solution of Eq..6. The relation between harmonic polylogarithms H A, with multiple index A including only and, and multiple polylogarithms of one variable Eq..6 is well known and follows from the proper definition see section in Ref. : Li k,k,,k n = H,,,,,,,,,,, ;..8 }{{}}{{}}{{} k times k times k n times 4 The Corollary 3 follows also from general properties of hypergeometric functions. 6

8 By continued iterations of Eq..6 and Eq..7, we have reproduced all coefficients of the ε-expansion of the Gauss hypergeometric function presented in Eq. 4.7 of 5. For the coefficients functions ρ 4,ω 4, and ρ 5, we find a more compact form. We also obtain the higher order terms ω 5 and ρ 6 of the ε-expansion. The results are 5 ρ 4 = a a 6 γ c ln 3 +cγ c a a ln Li +c Li 3 +cγ c a a S,, w 4 = c Li 4 a a cγ c a a Li + γc S,3 + cγ c a a S,,.9 ρ 5 = a a 4 γ3 c ln4 c 3 Li 4 cγc S,3 ccγ c a a ln Li 3 ccγ c a a S, γ c cγ c a a ln ln Li + S, w 5 = γc 3 a a S,4 c 3 Li 5 cγc S,3 ccγ c a a S 3, ρ 6 a a = γ 4 c + cγ c a a γ c Li S, γ c F cf ln5 + c 4 Li 5 + c γcs,3 + 6 γ c cγ c a a + ccγ c a a ln,.,. ln 3 Li + 3ln S, + 6ln S,3 γ c ln Li 3 + cli 4 c a c a cγ c a a ln S, +a a cγ c a a ln Li Li S, + F + cγ c a a γ cs,4 + c S 3,,. where we have introduced two new functions: dx F = x ln xli x,.3 F = There is an algebraic relation 6 between these two functions: dx x ln xli 3 x..4 F = F ln S,3 + S,3 Li S, ln ln S, 6 ln3 ln ln ln Li + ζ ln ln ζ S, ζ 3 Li ζ 5,.5 5 The FORM45 representation of these expressions can be extracted from Ref We are indebted to A. Davydychev for this relation. 7

9 where F = ζ 3 ζ ζ In this way, at the order of weight 5 one new function 7, F which is not expressible in terms of Nielsen polylogarithms, are generated by Laurent expansion of Gauss hypergeometric functions with integer values of parameters. In general, the explicit form of this function is not uniquely determined, and the result may be presented in another form by using different subset of harmonic polylogarithms... Half-integer values of of ε-independent parameters Let us apply a similar analysis for the second basis hypergeometric function a ε,a ε F +fε..6 In this case, the differential equation has the form dd + fε w = d d d d + a ε d d + a ε w,.7 with the same boundary condition w = and d d w = =. Using the ε-expanded form of solution, and noting that Eq..8, and in fact, Eq..7 is valid in each order of ε-expansion, we rewrite Eq..7 as d d d d w i = Let us introduce the new variable 8, a +a f d w i +a a w i..8 d y = + y, = 4y, = + y 4y, d d = y + y y d dy,.9 and define a set of a new functions ρ i y so that 9 and, as in the previous case, d d w i y + y y d w i y = y dy + y ρ iy,.3 ρy = d d w = ρ k yε k..3 k= 7 Compare with results of The form of this variable follows from analysis performed in 3, 3, 4. 9 We may note that F! + a ε,+a ε 3 +fε = + fε y + y X» ρk+ y k= a a ε k. 8

10 In terms of new variables y, the Eq..8 can be written as system of two first order differential equations: y d dy ρ iy = a +a y + y ρ i y + f y ρ i y+a a w i y, + y y d dy w ky = ρ k y..3 The solution of these differential equations for functions w i y and ρ i y has the form y dt f ρ i y = t a +a f + t y dt +a a t w i t, i, y w i y = ρ i t a +a w i y w i dt t ρ it, i..33 The point = transforms to the point y = under the transformation.9, so that the boundary conditions are w k =, k, ρ k =, k..34 The first several coefficients of ε-expansion can be calculated quite easily by using w y = and ρ y = : ρ y = w y =,.35a ρ y = lny H;y, a a.35b w y = a a ln y H,;y,.35c Continuing these iterations, we may reproduce the coefficients of the ε-expansion of the Gauss hypergeometric function.6 presented in Eq. 4. of 5. Since the length of the expressions obtained for the coefficient functions ρ 3 y,ω 3 y,ρ 4 y,ω 4 y,ρ 5 y,ω 5 y is similar to those published in Eq. 4. of 5, we don t reproduce them here. The higher order terms of ε-expansion are relatively lengthy and therefore will also not be presented here. Unfortunately, as in the previous case, we are unable to calculate the k-coefficient of ε-expansion without knowledge of previous ones. From representation.33 we deduce the following result: Corollary 4: The all-order ε-expansion of function.6 can be written in terms of harmonic polylogarithms HA y of variable y defined in.9 and multiple index A with entries taking values, and. M.Y.K. thanks to M. Rogal for pointing out a mistake in Eq. 4. of 5: In the ε term, the coefficient should be 3f a a instead of 3f a a. 9

11 This statement follows from the representation.33, the values of coefficients functions w k, k =,, see Eqs..,.7, properties of harmonic polylogarithms, and the relation between powers of logarithms and harmonic polylogarithms. Also, Corollary and Corollary 3 are valid for the hypergeometric function.6. We would like to mention that, in contrast to the Eq..6, Eq..33 contains a new type of function, coming from the integral ftdt/ + t. Another difference is that the first nontrivial coefficient function, ρ y, is equal to lny, instead of ln, as it was in the previous case. It was shown in that terms containing the logarithmic singularities can be explicitly factorised see Eqs. - in, so that the coefficient functions, w k y and ρ k y from Eq..33, have the following form: w k y = k j= k ρ k y = j= c s, σ,kln k j y Li σ s y Li σ s c s, σ,kln k j y Li σ s y Li σ s,,.36 where c s, σ,k and c s, σ,k are numerical coefficients, s and σ are multi-index, s = s, s n and σ = σ,,σ n, σ k belongs to the set of the square roots of unity, σ k = ±, and Li σ y is a coloured multiple polylogarithm of one variable 5, 6, 7, defined as s Li σ,σ,,σ k = s,s,,sn m >m > m n> m σ m σn mn m s ms n msn It has an iterated integral representation w.r.t. to three differential forms,..37 so that, ω = dy y, σ =, ω σ = σdy, σ = ±,.38 σy Li σ,σ,,σ k s,s,,s k y = ω s ω σ ω s ω σ σ ω s k ω σ σ σ k, σ k =..39 The values of coloured polylogarithms of unit argument were studied in 36, 37.. Zero-values of ε-dependent part of upper parameters In the case when one of the upper parameter of the Gauss hypergeometric function is a positive integer, the result of reduction has the more simple form compare with Eq..: Pb,c, F I,b + I ;c + I 3 ; = Q b,c, F,b;c; + Q b,c,,.4 where b,c, are any fixed numbers, P,Q,Q are polynomial in parameters b,c and argument, and I,I,I 3 are any integer numbers. In this case, it is enough to consider the following The proper algebraic relations for reduction are given in 5.

12 two basis functions: F, + aε; + cε; and F, + aε; 3 + fε;. The ε-expansion of this function can be derived from the proper solution given by Eq..7 or Eq..35, using the following relation: F,+a ε +fε + cε d = lim a a a ε d F a ε,a ε +cε = + cε and,+a ε +fε d a ε,a ε F 3 +fε = lim a a a ε d F +fε = + fε k= k= ρ k+ a a ρ k+ y a a a = a = ε k,.4 ε k,.4 where we used the differential relation, d a,b d F = ab + a, + b c c F + c, and the brackets mean that in the proper solution, we can put a =. The functions ρ k are given by Eq..7 and Eq..35, correspondingly. Due to Corollary 3, the limit a must exist. The case when both upper parameters are integers may be handled in a similar manner. Theorem is thus proved. 3. Some particular cases 3. The generalied log-sine functions and their generaliation For the case the variable y defined in.9, belongs to a complex unit circle, y = expiθ. In this case, the harmonic polylogarithms can be split into real and imaginary parts see the discussion in Appendix A of 4, as in the case of classical polylogarithms 33. Let us introduce the trigonometric parametriation = sin θ. In this case, the solution of the proper differential equations.6 and.33 can be written in the form θ ρ i θ = a +a c dφ sin φ θ cos φ ρ i φ+a a w i θ = and θ ρ i θ = a +a f w i θ = θ dφ sin φ cos φ w i φ cw i θ, i, dφ cos φ sin φ ρ i φ, i, 3. θ dφ sin φ θ cos φ ρ i φ f dφ cos φ θ sin φ ρ i φ+a a dφw i φ, dφρ i φ, i, 3.

13 respectively. In the first case, the solution of the system of equations 3. are harmonic polylogarithms with argument equal to sin θ. In the second case, the result contains the generalied log-sine functions 33, 38, 39, 4 and some of their generaliations studied in 4 see also Ref. 4. For illustration, we will present a first several terms of ε-expansion see the proper relations, Table I of Appendix C in Ref. 3: + a ε, + a ε F 3 + fε sin θ + fε = sin θ { } θ + ε γ f Ls π θ θl θ f Ls θ+θl θ +ε {fγ f Ls 3 θ+γ f γ f Ls 3 π θ fγ f Ls 3 θ +4fγ f Ls θl θ Ls π θl θ +θl θ l θ + 4f Ls θl θ 4γ f Ls π θl θ +f θl θ + γ f θl θ + 6 a a θ 3 + γ f γ f πζ } { +ε γ f a + a γ f Ls 4 π θ + f Ls 4 θ 3fγ f Lsc,3 θ 3 f γ f Ls 4 θ +4a a fcl 4 θ γ f Cl 4 π θ fcl 3 θθ γ f Cl 3 π θθ +fl θ + γ f L θ fγ f Ls 3 θ γ f γ f Ls 3 π θ fγ f Ls 3 θ +fl θ + γ f L θ γ f Ls π θ fls θ + a a γ f θ Ls π θ a a fls θθ 3 a a θ 3 fl θ +γ f L θ 4 3 θ fl θ+γ f L θ 3 γ f γ f πζ fl θ +γ f L θ +a a 3a + 3a 7fθζ 3 a + a γ f γ f πζ 3 } + Oε 4, 3.3 and F +ε a ε,a ε + fε sin θ = + a a ε θ } {fls θθ γ f Ls π θθ+4γ f Cl 3 π θ+4fcl 3 θ+3a +3a 7fζ 3 +ε { γ f γ f Ls 3 π θθ fγ f Ls 3 θθ+fγ f Ls 3 θθ+ 4 a a θ 4 fls θ γ f Ls π θ + γ f γ f θπζ } + Oε 3, 3.4 where L θ = ln cos θ, l θ = ln sin θ, The FORM representation of these expressions can be extracted from 46.

14 generalied log-sine function is defined as θ Ls k j θ = dφ φ k ln j k sin φ, Ls j θ = Ls j θ 3.5 and we use the notation Lsc,3 θ for the special combination see Eq..8 in Ref. 3 Lsc,3 θ = Ls 4 θ 3 Ls 4 θ + Ti 4 tan θ ln tan θ Ti3 tan θ + ln tan θ Ti tan θ where Ti N -functions is defined as 33 Ti N = Im Li N i = i 6 θ ln3 tan θ, 3.6 Li N i Li N i, Ti N = dx x Ti N x. 3.7 These functions receive special interest in physics through their role in the so-called single-scale diagrams, which depend only on one massive scale parameter. The massless propagator-type diagrams, bubble-type diagrams and propagator-type diagrams on mass shell all belong to this class. In particular, the single-scale diagrams with two massive particle cuts correspond to hypergeometric functions with value of argument equal to = /4 43,, 3. In this case, the value of the conformal variable y is equal to a primitive sixth root of unity, y = exp i π 3. In contrast to the case in multiple polylogarithms.6 of the primitive sixth root of unity studied in 44 and the more complicated case in coloured polylogarithms of the sixth root of unity studied by Broadhurst in 6, the physically interesting case corresponds to coloured polylogarithms of square root.39 harmonic polylogarithms with argument equal to primitive sixth root of unity. In this case, some new transcendental constants, in addition to studied in Ref. 44 will be generated. The set of independent constants up to weight 5 was constructed in Ref. 38, 3, Special cases: all-order ε-expansion in terms of Nielsen polylogarithms One advantage of a trigonometric representation used in the previous section is the theorem proved in Ref. 3 see also Ref. 4, that any generalied log-sine function 3.5 is expressible in term of Nielsen polylogarithms 34 only. Using this theorem, it was shown in Ref. 3, 3 that for the Gauss hypergeometric function, +aε F 3 +bε sin θ, 3.8 the Laurent expansion is expressible in terms of Nielsen polylogarithms only in the three cases, i b = ; ii b = a; iii a = b. Using the reduction algorithm 5, we can claim that the Laurent expansions of the following functions are also expressible in terms of Nielsen polylogarithms only: F I,I +ε +I 3 sin θ, F I,I +ε +I 3+ε sin θ, F I,I +ε +I 3+ ε sin θ,3.9 3

15 where I, I and I 3 are integers. It is interesting to analye this solution from the point of view of Eq. 3.. Due to fact that a =, the last term in Eq. 3. is identically equal to ero. In case i, only the first term survives, with integration kernel having the form dlncos φ. In case ii, only the second term survives, and the integration kernel has the form dlnsin φ, In case iii, the first and second terms can be reduced to the second case of a double argument. The statement about expressibility of inverse binomial sums in terms of log-sine function, proved in 3 see also 4, applies to all three of these cases. We can extend the class of Gauss functions whose ε-expansions are expressible in terms of only Nielsen polylogarithms by using algebraic relations 3 between of the fractional-linear arguments see section 3 in Ref. 5. The cases which may be expressed in this manner are summaried in Table I, where a, b, c are parameters of the Gauss hypergeometric functions F a,b;c; and I, I and I 3 are integer numbers: Table I a I I I I I I I +ε I +ε I +ε b +I +I +ε +I ε +I +ε +I +I +ε +I +I +ε +I +ε c +I 3+ε +I 3 +I 3+ε I 3 +ε I 3 +ε I 3 +ε I ++ε I ++ε I ++ε The results of this section 3. can be formulated as follows: Proposition : All cases of Gauss hypergeometric functions with half-integer values of parameters for which the all-order ε-expansion is expressible in terms of Nielsen polylogarithms only are described in Eq. 3.9 or the parameters shown in Table I. 4. Conclusion The main results of this paper is the proof of Theorem I, as stated also in the abstract. The proof includes two steps: i the algebraic reduction of Gauss hypergeometric functions of the type in Theorem I to basic functions and ii the iterative algorithms for calculating the analytical coefficients of the ε-expansion of basic hypergeometric functions. In implementing step i, the algebraic relations between basis functions with halfinteger values of parameters reduce all of the cases to the one basic function of type.a and its first derivative see details in 5. In step ii, the algorithm is constructed for integer values of parameters in Eq..6 and for basis Gauss hypergeometric function with half-integer values of parameters in Eq..33. This allows us to calculate the coefficients directly, without reference to multiple sums. It is interesting to note that the Laurent expansions of the Gauss hypergeometric functions with integer values of parameters are expressible in terms of multiple polylogarithms of one variable see Eq..6 or the Remiddi-Vermaseren harmonic polylogarithms with multiple index including only values and. The argument of the resulting functions coincides with the original variable of the hypergeometric function. For Gauss hypergeometric functions with half-integer values of 3 This can also be derived via the integral representation. 4

16 parameters, the coefficients of the ε-expansion produce the full set of harmonic polylogarithms, or coloured multiple polylogarithms of one variable, see Eq..37. These function depends on a new variable, related to the original variable by conformal transformation see Ref. 5. For special values of the argument of the hypergeometric function, <, the coloured multiple polylogarithms of one variable may be split into real and imaginary parts. This case has been discussed in section 3.. It was shown that the physically interesting case, representing single-scale diagrams with with two massive particle cuts, corresponds to coloured polylogarithms.39 with argument equal to a primitive sixth root of unity, y = exp i π 3. This gives an explanation of the proper basis of transcendental constants constructed in Refs. 38 and 3, and its difference from the proper basis of David Broadhurst 6. In the section 3., the subset of Gauss hypergeometric function is analyed, showing that the all-order ε-expansion is expressible in terms of Nielsen polylogarithms only. In particular, we have formulated the proposition that the only Gauss hypergeometric functions with half-integer values of parameters for which the all-order ε-expansion is expressible in terms of Nielsen polylogarithms only belong to one of the functions described in 3.9 or in Table I. Acknowledgments We are grateful A. Davydychev for useful discussion. M.Yu.K. is thankful to participants of conference Motives and Periods, University of British Columbia, Vancouver, June 5-, 6, 47, for interesting discussions. Special thank to Andreas Rosenschon for invitation and financial support and D. Kreimer and H. Gangl for enormously and useful discussions and suggestions. M.Yu.K. is very grateful to Laura Dolchini for moral support when paper was written. This research was supported in part by RFBR grant # 4--79, NATO Grant PST.CLG.9834 and DOE grant DE-FG-5ER4399. References N.N. Bogoliubov, D.V. Shirkov, Introduction to the Theory of Quantied Fields, A Wiley-Interscience Publication. John Wiley & Sons, New York-Chichester-Brisbane, 98; C. Itykson, J.B. Zuber, Quantum Field Theory, New York, McGraw-Hill, 98. E.E. Boos, A.I. Davydychev, Theor. Math. Phys ; A.I. Davydychev, J. Math. Phys ; J. Math. Phys ; A.I. Davydychev, J.B. Tausk, Nucl. Phys. B ; Phys. Rev. D arxiv:hep-ph/95443; F.A. Berends, M. Bua, M. Böhm, R. Scharf, Z. Phys. C ; S. Bauberger, F.A. Berends, M. Böhm, M. Bua, Nucl. Phys. B arxiv:hep-ph/949388; A.I. Davydychev, A.G. Groin, Phys. Rev. D arxiv:hep-ph/ A.I. Davydychev, M.Yu. Kalmykov, Nucl. Phys. B65 66 arxiv:hep-th/89. 5

17 4 F. Jegerlehner, M.Yu. Kalmykov, O. Veretin, Nucl. Phys. B arxiv:hep-ph/39. 5 F. Jegerlehner, M.Yu. Kalmykov, Nucl. Phys. B arxiv:hep-ph/ I.G. Halliday, R.M. Ricotta, Phys. Lett. B C. Anastasiou, E.W.N. Glover, C. Oleari, Nucl. Phys. B57 37 arxiv:hep-ph/997494; Nucl. Phys. B arxiv:hep-ph/99753; A.T. Suuki, E.S. Santos, A.G.M. Schmidt, J. Phys. A arxiv:hep-ph/48. 8 G. Rufa, Annalen Phys ; A.V. Kotikov, Phys. Lett. B ; ibid B ; ibid B ; E. Remiddi, Nuovo Cim. A arxiv:hep-th/ O.V. Tarasov, Nucl. Phys. Proc. Suppl arxiv:hep-ph/7. J. Fleischer, F. Jegerlehner, O.V. Tarasov, Nucl. Phys. B arxiv:hep-ph/373. L. Brucher, J. Frankowski, D. Kreimer, Mod. Phys. Lett. A arxiv:hep-th/93755; A.I. Davydychev, Nucl. Instrum. Meth. A arxiv:hep-th/5933. D.J. Broadhurst, J. Fleischer, O.V. Tarasov, Z. Phys. C arxiv:hep-ph/93433; D.J. Broadhurst, Z. Phys. C ; J. Fleischer, F. Jegerlehner, O.V. Tarasov, O.L. Veretin, Nucl. Phys. B arxiv:hep-ph/983493; O.V. Tarasov, Phys. Lett. B arxiv:hep-ph/ T. Gehrmann, G. Heinrich, T. Huber, C. Studerus, Phys. Lett. B arxiv:hep-ph/ E. Bejdakic, Y. Schröder, Nucl. Phys. Proc. Suppl arxiv:hep-ph/ A.B. Goncharov, Math. Res. Lett ; Math. Res. Lett D.J. Broadhurst, Eur. Phys. J. C arxiv:hep-th/ J.M. Borwein, D.M. Bradley, D.J. Broadhurst, P. Lisoněk, Trans. Am. Math. Soc arxiv:math.ca/ S. Weinierl, J. Math. Phys arxiv:hep-ph/43. 9 S. Weinierl, arxiv:hep-th/356. E. Remiddi, J.A.M. Vermaseren, Int. J. Mod. Phys. A5 75 arxiv:hep-ph/ S. Moch, P. Uwer, S. Weinierl, J. Math. Phys arxiv:hep-ph/83. M.Yu. Kalmykov, O. Veretin, Phys. Lett. B arxiv:hep-th/4. 3 A.I. Davydychev, M.Yu. Kalmykov, Nucl. Phys. B arxiv:hep-th/ M.Yu. Kalmykov, Nucl. Phys. Proc. Suppl arxiv:hep-th/ M.Yu. Kalmykov, J. High Energy Phys. 4: 6 56 arxiv:hep-th/68. 6

18 6 C.F. Gauss, Disquisitiones generales circa seriem infinitam + α β/γx +, Gesammelte Werke, vol.3, Teubner, Leipig, 83, pp A. Erdelyi Ed., Higher Transcendental Functions, vol., McGraw-Hill, New York, L.J. Slater, Generalied hypergeometric functions, Cambridge University Press, Cambridge A.F. Nikiforov, V.B. Uvarov, Special Functions of Mathematical Physics, Birkhäuser Verlag, Basel, 988; R.J. Yae, J.S. Dehesa, A.F. Nikiforov, J. Math. Anal. Appl ; E.J. Weniger, Comput. Phys. Comm. 33 ; R. Vidūnas, J. Comput. Applied Math Shu Oi, arxiv:math.nt/ K.T. Chen, Trans. A.M.S ; C. Kassel, Quantum Groups, Springer-Verlag, New York, D. Kreimer, Adv. Theor. Math. Phys. 3 3 arxiv:hep-th/ L. Lewin, Polylogarithms and Associated Functions North-Holland, Amsterdam, K.S. Kölbig, J.A. Mignaco and E. Remiddi, B.I.T ; R. Barbieri, J.A. Mignaco, E. Remiddi, Nuovo Cim. A 97 84; K.S. Kölbig, SIAM J. Math. Anal T. Huber, D. Maitre, Comput. Phys. Commun arxiv:hep-ph/ J.M. Borwein, D.M. Bradley, D.J. Broadhurst, Electron. J. Combin # R5 arxiv:hep-th/ J. A. M. Vermaseren, Int. J. Mod. Phys. A arxiv:hep-ph/9868; M. Bigotte, G. Jacob, N.E. Oussous, M. Petitot, Theoret. Comput. Sci. 73 7; Hoang Ngoc Minh, Nucl. Phys. Proc. Suppl J. Fleischer, M.Yu. Kalmykov, Phys. Lett. B arxiv:hep-ph/ A.I. Davydychev, Phys. Rev. D6 877 arxiv:hep-ph/994; A.I. Davydychev, M.Yu. Kalmykov, Nucl. Phys. Proc. Suppl arxiv:hep-th/ M.Yu. Kalmykov, A. Sheplyakov, Comput. Phys. Commun arxiv:hep-ph/4. 4 A.I. Davydychev, M.Yu. Kalmykov, arxiv:hep-th/3; M.Yu. Kalmykov, Nucl. Phys. B arxiv:hep-ph/ C. Anastasiou, S. Beerli, S. Bucherer, A. Daleo, Z. Kunst, arxiv:hep-ph/ J. Fleischer, M.Y. Kalmykov, A.V. Kotikov, Phys. Lett. B arxiv:hep-ph/99549; B E. 44 J.M. Borwein, D.J. Broadhurst, J. Kamniter, Exper. Math. 5 arxiv:hep-th/ J.A.M. Vermaseren, Symbolic Manipulation with FORM Computer Algebra Netherlands, Amsterdam, kalmykov/hypergeom/hyper.html