arxiv:hep-ph/ v2 13 Jul 2005

Transcription

1 ZU TH 3/5 hep-ph/5794 July 5 arxiv:hep-ph/5794v 3 Jul 5 HypExp, a Mathematica package for expanding hypergeometric functions around integer-valued parameters T. Huber and D. Maître Institut für Theoretische Physik, Universität Zürich, Winterthurerstrasse 9, CH-857 Zürich, Switerland ÆÙÖ Ñ ÖÒ Ò Ò Å Ù Ð Ø Ê Ù Ø Ö Ï Ö Ø Ø Ú Ö Ø Ö ÓÖÒ ÆÙÖ Å Ð Û Ö Ñ Ë Ð ÖÛ Ø Ë Å ÖÑÓÖ ÔÖ Ó ÃÓÖÒº aus: Friedrich Schiller (759-85), Das Ideal und das Leben Abstract We present the Mathematica package HypExp which allows to expand hypergeometric functions JF J around integer parameters to arbitrary order. At this, we apply two methods, the first one being based on an integral representation, the second one on the nested sums approach. The expansion works for both symbolic argument and unit argument. We also implemented new classes of integrals that appear in the first method and that are, in part, yet unknown to Mathematica.

2 Package summary Title of the Package HypExp Version. Package obtained from Licence none Computers Computers running Mathematica under Linux or Windows Operating system Linux, Windows Program language Mathematica Memory required to execute depending on the complexity of the problem Other Package needed the package HPL, included in the distribution External file required none. Keywords Hypergeometric functions, expansion Nature of the physical problem Expansion of hypergeometric functions around integer-valued parameters. These are needed in the context of dimensional regularisation for loop and phase space integrals. Method of solution Algebraic manipulation of nested sums and integral representation. Restrictions on complexity of the problem Limited by the memory available Typical running time Strongly depending on the problem and the availability of libraries.

3 Introduction As solutions of a large class of differential equations, hypergeometric functions PF Q appear in many branches of science. They appear, in particular, in particle physics during the calculation of radiative corrections to scattering cross sections in loop [ 6] or phase space [7 9] integrals. In the context of dimensional regularisation, the arguments of the hypergeometric functions have to be expanded in a small parameter around integer or rational arguments. Until recently, there was no systematic approach to the expansion of hypergeometric functions. The required expansions have been produced with a case-by-case approach. Recently a general algorithm has been developed [] for expanding hypergeometric functions and other transcendental functions systematically around their parameters. This algorithm was implemented [] in the framework of GiNaC []. Related work was also presented in Refs. [3, 4]. Until now, an implementation of the expansion of hypergeometric functions around their parameters was missing in the widely used computer algebra systems Mathematica [5] and Maple [6], allowing to use hypergeometric functions in connection with the multi-purpose features of these programs. The aim of this work is to provide an implementation of these expansions for Mathematica. This implementation is prepared in the form of a Mathematica package that aims to be tunable and user-friendly. This paper is structured as follows. In the next two sections, we present the two approaches we used in the package. Section 4 is devoted to hypergeometric functions of unit argument. In section 5 we describe the implementation of our methods in the package HypExp and provide examples of its usage. Finally, we conclude with a summary. Method of integration Among the hypergeometric functions (HF), the ordinary Gauß-hypergeometric function F (A, A ; B ; ) () appears most frequently in scientific calculations, and it is therefore worth to focus in particular on this type of functions [7]. We will calculate the ǫ-expansion of the F -functions up to and including order O(ǫ 4 ) by means of the algorithm described below which is based on the well-known integral representation of the F -functions. Applying this method to this subset of functions has by all means its benefits since it is faster and more efficient than the later described nested sums method. However, going to higher orders in the ǫ-expansion or to higher J F J -functions quickly reveals that the nested sums method is in general more powerful. Nevertheless, applying several independent methods also provides useful consistency checks.. F -algorithm We first set up our notation. In this work we consider with JF J ({A,..., A J }; {B,..., B J }; ) () A i a i + α i ǫ, B i b i + β i ǫ, a i, b i and α i, β i Ê. (3) We start the description of our algorithm by defining some subsets of. Let V : { Ê } (4) W : \ V. (5) 3

4 Furthermore, we state that the mapping f() : as a Moebius transformation is a bijective mapping and satisfies (6) f(v ) V f(w) W, and f(f()). (7) We finally turn our attention to the F -functions and will first collect some useful formulas [7,8]. We note that we can shift each of the parameters up or down by integer units via F (A, A ; B ; ) A B + + (A A ) A B + F (A +, A ; B ; ) + (A + )( ) (A B + ) F (A +, A ; B ; ), (8) F (A,A ; B ; ) A B + + (A A ) A B + F (A,A + ; B ; ) + (A + )( ) (A B + ) F (A,A + ; B ; ), (9) F (A, A ;B ; ) (B A A + ) B F (A, A ;B + ; ) B ( ) (B A + )(B A + ) F (A, A ;B + ; ), () B (B + )( ) F (A, A ; B ; ) B A + + (A A ) F (A, A ; B ; ) (A )( ) + A B (A )( ) F (A, A ; B ; ), () F (A,A ; B ; ) B A + + (A A ) (A )( ) F (A,A ; B ; ) + A B (A )( ) F (A,A ; B ; ), () F (A, A ;B ; ) (B ) [ B (A + A B + 3) ] (A B + )(A B + ) F (A, A ;B ; ) (B )(B )( ) (A B + )(A B + ) F (A, A ;B ; ). (3) These relations all stem from the formulas (B A ) F (A, A ; B ; ) + [ A B (A A ) ] F (A, A ; B ; ) + A ( ) F (A +, A ; B ; ), (4) (B A ) F (A, A ; B ; ) + [ A B (A A ) ] F (A, A ; B ; ) + A ( ) F (A, A + ; B ; ), (5) B (B )( ) F (A, A ; B ; ) + B [ B (B A A ) ] F (A, A ; B ; ) + (B A )(B A ) F (A, A + ; B + ; ). (6) 4

5 There is yet another class of relations between F -functions, namely the relations of Gauß between contiguous functions [7]. Their inclusion would lead only to a minor simplicifation here, and thus we let them serve as a check for our results rather than implementing them in our algorithm. From the relations (8) (3) we conclude that the knowledge of the ǫ- expansions of the F -functions whose integer parts {a, a, b } of the parameters read {,, } {,, } {,, } (7) {,, } {,, } {,, } is sufficient in order to derive the ǫ-expansion of any F -function with {a, a, b } being integer-valued. But even this small set of functions can be reduced further by means of Kummer relations [7,8]. The relevant Kummer relations read F (A, A ; B ; ) ( ) B A A F (B A, B A ; B ; ), (8) F (A, A ; B ; ) ( ) A F (A, B A ; B ; ), (9) F (A, A ; B ; ) ( ) A F (B A, A ; B ; ). () and relate both the functions {,, } and {,, } to the function {,, }, such that we can get along with a mere four functions, namely {,, } {,, } {,, } {,, }. () For completeness, we mention that Eq. (8) holds true for all, whereas Eqs. (9) and () are only valid for W. The sets (7) and () of basic hypergeometric functions have, however, one major drawback. In order to express a general F -function solely in terms of functions from these sets, repeated application of Eqs. (8) (3) is required and additional negative powers of ǫ might be generated in prefactors via this procedure. It is therefore necessary to know the ǫ-expansions of the basic hypergeometric functions to higher order than is sought by the F -function in question. In order to avoid this disturbing feature we consider an extended set of basic hypergeometric functions. The extended set has three major subsets. In the first subset we collect those basic HF s that contain only positive integer parts, namely {,, } {,, } {,, } {,, } () {,, } {,, } {,, } {,, }. The second subset contains those basic HF s in which b. It reads {,, } {,, } {,, }. (3) The third subset finally contains those basic HF s in which b : {,, } {,, } {,, } {,, } {,, } {,, } (4) {,, } {,, } {,, } {,, }. Although some functions in this set might not be considered independent since they are related via Kummer relations, we will consider this set as basic since it will allow us to implement the algorithm described below efficiently and conveniently. The goal of the latter is to express a general F -function entirely in terms of functions from the set () (4) by repeated application of both the equations (8) (3) and the Kummer relations (8) (). Before we start, we mention that throughout the algorithm we make use of the symmetry A A after each step in order to ensure that we always have a a. 5

6 . We start the reduction of our F -function in question by applying Kummer relations such that the sum a + a + b (5) gets minimied. Especially for high absolute values of the parameters this procedure shortens the algorithm significantly.. Then, if b <, we apply () repeatedly to all HF s with b <. This step ensures that from now on we only have to deal with functions in which b. For the rest of the algorithm we distinguish two cases, namely b and b. The further steps for b and b are illustrated by the flowcharts in Figures and respectively. The ambitious reader is invited to verify that at the end of this algorithm only F -functions from the set () (4) appear and that no negative power of ǫ has been generated at any intermediate step. To conclude this section we remark that for the most frequent case in which the three parameters a, a, and b are all non-negative, only the last column of Figure has to be considered.. Expansion of the basic F -functions Now that we went through the algorithm for the F -functions in great detail we have to explain how the ǫ-expansions of the basic HF s from the set () (4) are obtained. For the functions in Eq. () we adopt the integral representation [7] F (A, A ; B ; ) Γ(B ) Γ(A )Γ(B A ) du ua ( u) B A ( u) A, (6) which we must restrict to W and B > A >. Since this approach will be based on the requirement that the integration over u and the expansion in ǫ commute, we have to set up the additional condition b > a >. One recognies immediately that from the set () only the functions {,, } and {,, } satisfy the latter inequality. The parameters of the other six functions first have to be shifted by means of Eqs. (9) and () until a convergent integral representation is obtained for each of them. The subsequent expansion of the integral representation in ǫ and how one solves the occurring integrals is covered in section.3. The functions in Eqs. (3) and (4), all of which contain at least one negative parameter, are now, for the sake of obtaining their ǫ-expansion, expressed in terms of functions from the set (). This is again done by appropriate application of the Kummer relations (8) () as well as Eqs. (8) (). The shift of parameters as described in the preceding two paragraphs is now unavoidably accompanied by the advent of negative powers of ǫ in certain prefactors. To be more precise, we must expand the functions {,, } and {,, } up to and including order O(ǫ 4 ) and the other six functions of Eq. () to order O(ǫ 5 ). This might at first glance seem peculiar, but it turns out that at the respective highest order in ǫ integrals of the same type appear. We conclude from this that, with the tools provided here, it is in principle possible to expand a certain class of F -functions even to order O(ǫ 5 ). However, for reasons of simplicity and clearness we apply the method of integration throughout only for the expansion of F -functions up to and including order O(ǫ 4 ). 6

7 b a and a a and a < a < and a < a and a < a and a a and a 7 {,, } apply (9) apply (8) apply () repeatedly to all HF s with a < repeatedly to all HF s with a < repeatedly to all HF s with a > apply (8) repeatedly to all HF s with a < apply () repeatedly to all HF s with a > {,, } {,, } {a,, } {a,, } {,, } {,, } {,, } {,, } {,, } {,, } {,, } apply () repeatedly to all HF s with a > {,, }, {,, }, {,, }, {,, } apply () repeatedly to all HF s with a > {,, } {,, } {,, } {,, } Figure : F -algorithm for b

8 b 8 a < and a < apply (9) repeatedly to all HF s with a < {a,, } {a,, > } {a,, } apply (8) repeatedly to all HF s with a < {,, } {,, } apply (8) repeatedly to all HF s with a < {,, > } {,, > } apply (9) once a and a < apply (8) repeatedly to all HF s with a < {,, } {,, >} {,, } apply () repeatedly to all HF s with a > {,, } {,, } apply (8) once {,, > } {,, > } a and a apply () repeatedly to all HF s with a > apply () repeatedly to all HF s with a > {,, } {,, } {,, } {,, } {,, > }, {,, > } apply (3) repeatedly to all HF s with both (b > a + ) and (b > a + ) functions from the set () Figure : F -algorithm for b

9 .3 Integrals and their algorithm.3. General aspects Expressing a hypergeometric function F (a + α ǫ, a + α ǫ; b + β ǫ; ) with b > a > in terms of its integral representation according to Eq. (6) and subsequently expanding in ǫ yields integrals that have the general form I (χ, χ, χ 3, χ 4, χ 5, ) : du uχ ln χ (u) ln χ3 ( u) ln χ4 ( u) (u ) χ5, (7) where the χ i are non-negative integers and W. The results of these integrals contain rational functions, logarithms, polylogarithms Li n, Nielsen polylogarithms S n,p, and Harmonic polylogarithms H m,...,m k, all of which will be explained in more detail in sections 3.. and in appendix A.. Let us define the weight w of the integral by w : χ + χ 3 + χ 4 + δ χ5,. (8) The weight w is related to the weight of the harmonic polylogarithms defined later, namely any integral of weight w can be expressed in terms of harmonic polylogarithms of weight w or less. In order to guarantee that the ǫ-expansion of any F -function up to order O(ǫ n ) can be performed, the computation of all integrals with weight w up to n + is required..3. Description of the algorithm In order to make the computation more efficient we now show that any integral I (χ, χ, χ 3, χ 4, χ 5, ) can be expressed in terms of integrals with χ χ 5 and of at most the same weight as the original one:. In the case χ χ 5, the integral I (χ, χ, χ 3, χ 4, χ 5, ) can, via the replacement u (u ) + (9) in the numerator and subsequent cancellation of the corresponding denominators, be expressed as a linear combination of integrals of the form I (, χ, χ 3, χ 4, χ 5, ), I (, χ, χ 3, χ 4, χ 5, ),..., I (, χ, χ 3, χ 4, χ 5 χ, ). (3). In the case χ > χ 5, we can apply the same steps as before and eventually arrive at integrals of the form I (χ χ 5, χ, χ 3, χ 4,, ), I (χ χ 5, χ, χ 3, χ 4,, ),..., I (, χ, χ 3, χ 4,, ), I (, χ, χ 3, χ 4,, ),..., I (, χ, χ 3, χ 4, χ 5, ). (3) In other words, the replacement rule (9) for u yields integrals in which at least one of the entries χ and χ 5 is ero. 3. We continue our reduction by manipulating integrals of the form I (χ, χ, χ 3, χ 4,, ) with χ >. By recalling the relation du d du [uχ ( u)ln χ (u) ln χ3 ( u) ln χ4 ( u)] (3) 9

10 and explicitly taking the derivative of the product, one derives the relation (χ + )I (χ, χ, χ 3, χ 4,, ) χ I (χ, χ, χ 3, χ 4,, ) χ I (χ, χ, χ 3, χ 4,, ) +χ I (χ, χ, χ 3, χ 4,, ) χ 3 I (χ, χ, χ 3, χ 4,, ) + χ 4 I (χ, χ, χ 3, χ 4,, ) χ 4 I (χ +, χ, χ 3, χ 4,, ). (33) Repeated application of steps. 3. finally yields an expression which contains only integrals with χ and of at most the same weight as the integral we started with. 4. The remaining task is now to subsequently lower χ 5. By repeated application of the relations I (, χ, χ 3, χ 4, χ 5, ) + ( )χ5 (χ 5 )! for χ 5 > and χ 4 v { χ 5 k χ 4! (χ 4 v + )! [ ( ) χ ] 5 I (, χ k, χ 3, χ 4, k, ) [ χ5 j ( ) χ5 ( ) j j j v] d χ5 } d I (, χ χ 5, χ 3, χ 4 v +,, ) (34) I (, χ, χ 3, χ 4,, ) I (, χ, χ 3, χ 4,, ) + χ 4 + d d I (, χ, χ 3, χ 4 +,, ) (35) one eventually arrives at an expression of the desired form. 5. We can, however, reduce the number of distinct integrals even further. The transformation u u in the integrand of (7) allows us to represent integrals with χ < χ 3 in terms of integrals with χ > χ 3. By means of the relation ln [ ( u)] ln ( ) + ln( u), (36) which holds true for all W and < u <, one easily derives the formula I (, χ, χ 3, χ 4,, ) χ 4 m ( ) χ4 ln χ4 m ( ) I(, χ m 3, χ, m,, ). (37) Since this formula also transforms the argument, we derived relations between polylogarithms that allow to simplify again these arguments. The argument transformations of the polylogarithms are described in section 5. and in appendix A.. To summarie, for a given weight w the set of basic integrals consists of all I (, χ, χ 3, χ 4,, ) with χ + χ 3 + χ 4 w. (38) Integrals thereof with χ < χ 3 can be reexpressed in terms of integrals with χ > χ 3 and argument /( ) and are therefore easily obtained once the argument transformations of the polylogarithms are taken into account.

11 We implemented this algorithm in the function HypExpInt. This function will be explained in detail in section 5.. For the actual calculation of the ǫ- expansion of F -functions we need a bit less than is provided by the algorithm and by the function HypExpInt. We have: χ, (39) χ 5. (4) The first inequality arises from the fact that we have to shift the parameters of our basic hypergeometric functions, Eq. (), via the relations (9) and () in order to obtain convergent integral representations. In this procedure, b, which eventually determines the highest value for χ to occur, assumes values up to. Similarly, the parameter a determines the highest value for χ 5 that can show up. With this in mind one derives the second inequality directly from the collection () of basic hypergeometric functions..3.3 Integrals of unit argument Putting in the integrals I (χ, χ, χ 3, χ 4, χ 5, ) immediately turns our attention to an other type of integrals that we considered useful to implement. We define the function U (n, m, p) by U (n, m, p) : du ln n (u) ln m ( u) u p (4) with p and n, m being non-negative integers. In the case p < the inequality m + p has to be satisfied in order to yield a convergent integral. We demonstrate below how any convergent integral U (n, m, p) can be expressed in terms of U (,, ) and integrals of the form U (n, m, ).. We start by considering the case p < ; n > and m + p. Repeated application of U (n, m, p) n U (n, m, p) + m U (m, n, ) p + p + + m p U (n, m, τ) (4) p + τ leaves us with integrals of the form U (n, m, ) and U (, m, p), where the latter type still happens to have p < and m + p. These integrals get reduced via the recursion relation U (, m, p) to integrals with p. m p U (, m, κ) (43) p +. We now proceed with the case p, in which the function U is symmetric in n m and therefore n m can always be achieved. Applying this in turn with the formula U (n, m, ) n + n m σ m λ κ ( ) m σ m! σ! ( ) m λ m! λ! U (n, σ, ) U (n, λ, ) (44) for {n, m} {, } eventually yields integrals of the desired form.

12 3. In the case p > we apply the same trick as in step 3 of section.3. and derive from the recurrence relation du d du [up ( u)ln n (u) ln m ( u)] (45) (p + )U (n, m, p) p U (n, m, p ) + n U (n, m, p ) n U (n, m, p) m U (n, m, p). (46) Repeated application of these steps finally yields an expression which only contains integrals of the demanded form. U(,, ) is trivially found to be unity. A nice algorithm for the computation of integrals of the form U(n, m, ) is given in section of Ref. [9] and will not be repeated here. The evaluation of integrals U (n, m, p) with p and n, m being nonnegative integers can be called with the function HypExpU, and an example can be found in section 5.. The connection to I (χ, χ, χ 3, χ 4, χ 5, ) in is given by χ ( χ I (χ, χ, χ 3, χ 4, χ 5, ) ( ) j χ5 j j ) U (χ 3 + χ 4, χ, j χ 5 ) (47) for χ χ 5. By means of this relation the function HypExpInt can be directly called with unit argument, see section Nested sums method 3. Definitions and auxilary functions In this section, we briefly review the S and Z sums and their properties introduced in []. Their definitions are Z(n, {m,...,m k }, {x,..., x k }) S(n, {m,...,m k }, {x,..., x k }) n i i i n i i i i k i k i k i k x i i m x i i m... xi k k i m, (48) k k... xi k k i m, (49) k k or the equivalent recursive definitions Z(n, {m,..., m k }, {x,..., x k }) S(n, {m,..., m k }, {x,..., x k }) Z(n, {}, {}) S(n, {}, {}) n i n x i i m x i i m i Z(i, i,...,k, x,...,k )(5) S(i, i,...,k, x,...,k ) (5) {, n (5), n < {, n >, n, (53) where we introduced the short-hand notation m,...,k {m,..., m k }.

13 The number k is called the depth and the sum of the m i s is called the weight of the nested sum. The S and Z sums are related, since their expressions differ only by the upper summation limits in the recursion relation. Using i j i<j +δ ij, we can convert S into Z sums and vice versa by means of recursive application of the identities S(n, m,..., x,... ) Z(n; m,..., x,... ) n x i i x i m i m i i S(i, m 3,..., x 3,... ) i +S(n, {m + m, m 3,... }, {x x, x 3,... }), (54) n i x i i m i i x i i m Z(i, m 3,..., x 3,... ) Z(n, {m + m, m 3,... }, {x x, x 3,... }), (55) and reconstructing Z sums from the x i /i m with help of the definition (48). The Z and S sums form two algebrae, that means that a product of Z sums is expressible as a sum of Z sums. Using Z(n; m,...,k, x,...,k ) Z(n, m,...,l, x,...,l ) n x i Z(i, m,...,k, x,...,k )Z(i, m,...,l, x,...,l) i + + i m n i n i x i i m Z(i, m,...,k, x,...,k )Z(i, m,...,l, x,...,l ) (x x )i Z(i, m,...,k, x,...,k )Z(i, m i m+m,...,l, x,...,l ) (56) recursively and reconstructing Z sums from the x i /i m with (48) proves the claim. Likewise, a product of S sums is a linear combination of S sums. Since products are defined for equal upper summation limit n, it is useful to have a relation between sums Z(n,... ) of different n s. The recursive use of the following formulae allows to change the upper summation boundary of a Z or S sum, and, doing so, to bring all the nested sums to the same upper summation limit. This is called syncroniing the sums in [] Z(n + c, m,...,k, x,...,k ) Z (n, m,...,k, x,...,k ) c + x j j x n (n + j) m Z(n + j, m,...,k, x,...,k ), (57) S(n + c, m,...,k, x,...,k ) S(n, m,...,k, x,...,k ) c + x j x n (n + j) S(n + j, m,...,k, x,...,k ). (58) m j Since the definition of S and Z sums is for denominators of the form i m, it will also be useful to convert sums of the form n i x i Z(i,...) (i + c) m to the form of the definition, that is, with c. If the depth of the sum is ero, we use n i x i (i + c) m x n i x i (i + c ) m c m + x n (n + c) m, (59) 3

14 and otherwise we use n i x i (i + c) m Z(i, m,...,k, x,...,k ) x n i x i (i + c ) m Z(i, m,...,k, x,...,k ) n x i x i (i + c) m i Z(i, m,...,k, x,...,k ) m i + xn (n + c) m Z(n, m,...,k, x,...,k ). (6) 3.. Relations to other functions Special cases of Z and S sums are related to other functions. For finite upper summation limit we have Z (n, m,...,k,,...,) Z m,...,m k (n) (6) S (n, m,...,k,,...,) S m,...,m k (n), m i > (6) where the Z m,...,m k (n) are the Euler-Zagier sums [,] and the S m,...,m k (n) are the harmonic sums []. For infinite upper summation limit, we have the following identities Z(, m,...,k, x,...,k ) Li mk,...,m (x k,...,x ) (63) Z(, m,...,k,,...,) ζ(m k,..., m ) (64) Z(, m,...,k, {x,,...,}) H m,...,m k (x) (65) Z(, {n +,,..., }, {x,,...,}) S }{{}}{{} n,p (x) (66) p p Li mk,...,m (x k,...,x ) are the Goncharov multiple polylogarithms [3]. A special case of the Goncharov multiple polylogarithms are the harmonic polylogarithms (HPL) H m,...,m k (x) of Remiddi and Vermaseren [5]. ζ(m k,...,m ) are the multiple eta values [4]. The HPL s can be reduced to classical polylogarithms (Li n (x)) [9] and Nielsen polylogarithms (S n,p (n)) [6] up to weight 4 in our case. Since for hypergeometric functions we only have one free variable, we will only use the last two identities. The package HypExp uses the package HPL [7] to deal with the harmonic polylogarithms. Here one can see the usefulness of the Z sums as a connection between the Euler-Zagier sums (which appear in the expansion of the Γ function, as will be shown in section 3..) and the harmonic polylogarithms which will eventually appear in the coefficients of the expansion. Since for the expansion around integer-valued parameters of the HF, there appear only HPL s with positive indices. 4

15 3.. Expansion of the Γ function Γ functions can be expanded around integer values as follows [] Γ(a + αǫ) Γ( + αǫ)γ(a) + (αǫ) j Z,..., (a ) j }{{} (67) j + (αǫ) j Z Γ(a + αǫ) Γ( + αǫ)γ(a),..., (a ) j }{{} j + ( αǫ) j S Γ( + αǫ)γ(a),..., (a ) j }{{} (68) j for a integer and a >. For negative (or vanishing) a, one has to use the identity xγ(x) Γ(x + ) so that and Γ( n + αǫ) Γ( + αǫ) αǫ Γ( + αǫ) αǫ Γ( + αǫ) αǫ Γ( n + αǫ) αǫ Γ( + αǫ) n j j + αǫ ( ) n Γ( αǫ) Γ(n + αǫ) ( ) n + (αǫ) j S Γ(n + ),..., (n) (69) j }{{} j n ( j + αǫ) j αǫ Γ(n + αǫ) ( )n Γ( + αǫ) Γ( αǫ) αǫ Γ( + αǫ) ( )n Γ(n + ) + ( αǫ) j Z,..., (n). (7) j }{{} j 3. Description of the algorithm The algorithm described in this section is the adaptation of the algorithm of type A of [] for the special case of hypergeometric functions. The definition of the hypergeometric series is given by JF J ({A,... A J }, {B,...,B J }, x) (A ) i... (A J ) i x i (B i ) i...(b J ) i i! (A ) i... (A J ) i x i + (B i ) i... B J ) i i! + Γ(B )...Γ(B J ) Γ(A + i)...γ(a J + i) Γ(A )...Γ(A J ) Γ(B }{{} i + i)...γ(b J + i)γ(i + ) xi }{{} E I 5 (7)

16 Where (a) i a(a + )...(a + i ) is the Pochhammer symbol. We denote the coefficients in front of the sum with E and the sum itself with I. The parameters A and B are of the form A i a i + α i ǫ B i b i + β i ǫ a i, b i and α i, β i Ê. In order to give the ǫ-expansion (to order n) of the hypergeometric function, one has to expand the product of E and I to order n. Let us have a look at the required depth of the ǫ-expansion for these factors. For each negative b j, one gets a factor ǫ in E, and factors ǫ in I, but only for the values,..., b j of i. For vanishing b j s one gets a factor ǫ in E but no factors of ǫ in I. For negative a j s, one gets factors of ǫ in E and factors ǫ in I for i,..., a j, For vanishing a j s, one gets factors of ǫ in E but no factors ǫ in I. This means that for the expansion of the hypergeometric function to order n we have to compute the expansion of E to order and for I to order n E n + #(a j < ) (7) n I n + #(b j ) #(a j ). (73) Note that the above argument explains why the expansion of F (α ǫ, α ǫ, b + βǫ, x), b > starts with + O(ǫ ), whereas that of (for example) has a ǫ term. We treat first the prefactor F ( + α ǫ, + α ǫ, βǫ, x) E Γ(B )...Γ(B J ) Γ(A )...Γ(A J ). (74) We expand the Γ functions with the formulae of section 3.., which leads to E Γ( + β ǫ)...γ( + β J ǫ) Γ( + α ǫ)...γ( + α J ǫ) f (A,..., A J, B,...,B J, ǫ). (75) We have to expand the product of Z sums appearing in f to the order n Γ E n E + #(b ) n + #(a < ) + #(b ) (76) because of the factors ǫ for each negative b. Since the factors Γ( + β i ǫ) and Γ( + α i ǫ) will cancel those of the expansion of the Γ functions of the sum I we will only consider Γ( + α ǫ)...γ( + α J ǫ) Γ(B )...Γ(B J ) Ẽ Γ( + β ǫ)...γ( + β J ǫ) Γ(A )...Γ(A J ) f (A,..., A J, B,..., B J, ǫ). (77) The coefficients of the ǫ-expansion of Ẽ are Z or S sums with finite upper summation limit. These are rational functions that can be computed easily using the recursive definition. 6

17 We now turn to the sum I i Γ(A + i)...γ(a J + i) Γ(B + i)...γ(b J + i)γ(i + ) xi. (78) First one makes use of the identity xγ(x) Γ(x+) to bring all Γ(A j +i) s and Γ(B j + i) s to the form Γ(i + ξǫ) times some rational functions of i. Then one expands the Γ functions with the formulae of section 3... As one knows that there will appear terms with, at worst, as many factors /ǫ as there are strictly negative a j s, one has to expand the Γ-functions to the order n Γ I n I + #(a j < ) n + #(b j ) #(a j ) + #(a j < ) (79) in order to be able to get I to the required order n I. The result of this expansion is a product of the factors Γ( + ξǫ) from the expansion of the Γ(i + ξǫ). We use them to cancel those of E. From now on, we denote by Ĩ the sum I without these factors. a product of ǫ-expansions with coefficient Z,..., (i ) or S,..., (i ). These products are treated as above, first converting the S sums into Z sums, then expanding the products of Z sums into a sum of single Z sums. In this case one cannot calculate numerical values for the Z sums, due to the i in the argument. a rational function R(i, ǫ) of i and ǫ from the factors x in the reduction xγ(x) Γ(x + ). This is reduced by means of expansion into partial fractions to a sum of coefficients times single denominators. The next step is to bring the factors /(i + c) n appearing in R(i, ǫ) to the form i m suitable for applying the definition of the Z sums. The method is presented in section 3.. At this point, one can perform the last summation over i, as described in the following section. The result is a sum of Z sums with upper summation limit equal to infinity. These can be converted to more common functions, as described in section 3... To obtain the full expansion, one has to multiply the expansions of Ẽ and Ĩ and to add unity to the result. 3.. Last summation The algorithm reduces the sum Ĩ to a sum of terms of the form i x i (i + j + αǫ) n Z (i, m,...,k, {,..., }). Here we have to distinguish several cases. The general strategy is to perform the summation over i until the denominator is positive for each i, then to simplify the remainder to the form of the definition of a Z sum. We list below the different kinds of terms that can appear and for each of them the way it is processed. Z sums of argument i times one denominator with negative offset Z(i,...). (i j + αǫ) m 7

18 Here one has to perform the sum up to the j-th term, according to i (i j + αǫ) m Z(i,... ) j k + xj x k Z(k,...) (k j + αǫ) m Z(j,...) (αǫ) m x k+j + Z(k + j,...). (8) (k + αǫ) m k One then expands the first two terms to the required order in ǫ. The occurring Z sums can be evaluated, as their upper summation limits are finite. We now take a closer look at the last sum x k+j Z(k + j,...) (k + αǫ) m ( ) x j x k ( αǫ) l (m + l )! k l+m Z(k + j,... ).(8) l!(m )! k k l Here one only keeps terms up to the required order in ǫ in the sum over l. The Z sum in the last term has to be syncronied down to argument i, as described in section 3.. Z sums of argument i times one denominator without offset x i Z(i,...). in If n is positive, following the definition of the Z sums, one gets i x i i n Z(i, m,...,k, {,...,}) Z(, {n, m,..., m k }, {x,,..., }) H n,m,...,m k (x). (8) In many cases, the harmonic polylogarithm can be expressed in terms of more common functions. If n is negative, one has to interchange the first two summations i m x i Z (i, m,...,k, x,...,k ) i i i i i m x i x i i m i x i i m Z (i, m 3,...,k, x 3,...,k ) Z (i, m 3,...,k, x 3,...,k ) ii + i m x i. (83) The last sum ii + i m x i ii + ( x ) m ( x i x ) m x i+ x x x (84) is a polynomial of degree m in i. The coefficient of each power of i is a rational function in x. The i k of this polynomial can be combined with the 8

19 i m from the definition of the Z-sum. We can also factor out an x i+. The remainder is then a polynom in i with rational coefficient functions of x. The result is then i m x i Z (i, m,...,k, x,...,k ) i x (xx ) i Z (i, m 3,...,k, x 3,...,k ) C j (x)i j. (85) i i m We are now left over with a sum of terms of the form C(x) (xx ) i Z (i, m 3,...,k, x 3,...,k ). i i m Recursive use of this formula leads either to Z sums with a denominator with positive m or to a Z sum multiplied by a denominator with negative exponent n but with ero depth. For the latter we can apply x i i m Z (i, {}, {}) x i i m Li m (x) i i which is also a rational function in x. Z sums of argument i without denominators x i Z(i, m,...,k, {,...,}) ( x ) m x x x, (86) The summation can be performed by exchanging the first two summations x i Z(i, m,...,k, {,...,}) i i i i i x i i i m i m i m Z(i, m 3,...,k, {,..., }) Z(i, m 3,...,k, {,...,}) ii + Z(i, m 3,...,k, {,...,}) xi+ x x x Z(, m,...,k, {x,,...,}) x x H m,...(x) (87) Single denominator with negative offset i x i (i j + αǫ) m Here we perform the sum explicitly for the j first terms i x i (i j + αǫ) m j k x k (k j + αǫ) m + xj (αǫ) m + xj 9 k x i x k (k + αǫ) m, (88)

20 we can simplify the last summation if we are only interested in the expansion to the order n: ( x k (k + αǫ) m n ) ( αǫ) l x k (m + l )! k l+m + O(ǫ n+ ) l!(m )! k k l n (m + l )! ( αǫ) l x k l!(m )! k l+m + O(ǫn+ ) l k n l Note that this works for positive and negative m. Single denominator with positive offset We use the formula i x i (i + j) m x j (m + l )! ( αǫ) l Li m+l (x) + O(ǫ n+ ) (89) l!(m )! i i x i (i + j) m. x j (Li m (x) ( x i+j (i + j) m x i x j i m i ) j x i i m i j i ) x i i m. (9) 4 Hypergeometric functions of unit argument The hypergeometric series JF J ({A,..., A J }; {B,..., B J }; ) n (A ) n... (A J ) n n (B ) n...(b J ) n n! (9) converges in only if the condition J B t t is satisfied. This implies that, if J s b t t J A r > (9) r J a r >, (93) the expansion in ǫ commutes with the procedure of taking the limit, and the series expansion of the function J F J ({A,...,A J }; {B,..., B J }; ) around ǫ has a well-defined finite radius of convergence. We therefore call the case s > non-critical. The case s, on the other hand, will be referred to as critical since this case requires more care and additional explanation on its treatment. By means of an algorithm based on partial fractions it is possible to express a HF J F J of unit argument and value s in terms of J hypergeometric functions J F J, also of unit argument, but of value s + or higher. We outline this procedure for the case in which no two of the B i are equal. Always assuming convergence, we start with the series expansion (9) of the function r JF J ({A +,..., A J + }; {B +,...,B J + }; ) (94)

21 and multiply and divide therein by the fraction (A + n)... (A J + n) (B + n)... (B J + n). (95) The inverse of the above expression gets combined appropriately with the Γ- functions of the Pochhammer symbols, whereas the multiplied one gets expanded into partial fractions, yielding J n B µ + µ J J A + τ J σ (B τ A σ ) (B τ + n) J λ;λ τ (B τ B λ ). (96) The term linear in n gives back the function (94), which then cancels on both sides of the equation. The constant term is proportional to the function JF J ({A,...,A J }; {B,...,B J }; ) (97) for which we are seeking and for which we can now solve the equation. Each term of the last sum in (96) contains a hypergeometric function of value s +. After some intermediate steps, we arrive at [ J ] J B µ A JF J ({A,..., A J }; {B,..., B J }; ) µ J J τ J σ (B τ A σ ) λ;λ τ (B τ B λ ) B τ J F J ({A,...,A J }; {B,..., B τ, B τ +, B τ+,...,b J }; ). (98) The case in which there is any combination of equal B i s can be treated analogously. The linear and the constant term in the expression (96) remain unchanged, only the last sum will look different, and the final expression (98) will contain HF s of value s + or higher. Repeated application of this algorithm allows to express a critical HF as a linear combination of non-critical ones, the two expressions being related via analytic continuation. Since the analytic continuation is unique and the expression obtained by the algorithm has a well-defined expansion around ǫ, we can associate the expansion at hand also with the original critical HF. In this sense the ǫ-expansion of a critical HF has to be understood and the user must be aware of this feature. The same phenomenon happens, by the way, for the well-known Γ-function. To elaborate a bit more on this we state that the above procedure works only for J > and that the case J is simpler: F (A, A ; B ; ) Γ(B )Γ(B A A ) Γ(B A )Γ(B A ). (99) Again, the series expansions of the Γ-functions have to be understood in the sense of analytic continuation. The crucial property of a critical function is that one can find at all a value for ǫ such that the series (9) converges, i.e. J β t t J α r () r in the case s. If s and () yields ero, the expression is divergent for all values of ǫ and no remedy can be found. The way in which the HFs of unit argument are implemented in the package is described in section 5..

22 5 The Mathematica package HypExp We implemented the results of the preceding sections in the Mathematica package HypExp. It allows to expand arbitrary J F J -functions to arbitrary order in a small quantity around integer parameters and can be obtained from [8]. The results are displayed in terms of rational functions, logarithms, polylogarithms, Nielsen polylogarithms, and harmonic polylogarithms. The results given by these functions are not systematically simplified using Simplify, since the simplification might take longer than the expansion itself, in particular for expansions to high orders. A Simplify might produce a more compact result. Furthermore, the package HypExp should be loaded at the beginning of the session. 5. Functions, commands and symbols added The package HypExp adds two new symbols $HypExpPath is the path where the HypExp package is installed. $HypExpFailed is the symbol returned by the package in case of failure. The package adds the following functions HypExp[HypergeometricF[...,x]ǫ,n] gives the ǫ expansion of the enclosed hypergeometric function. The function HypExp applied to anything else but a HF will leave it intact. Therefore one can map it onto an expression containing hypergeometric functions, and only the HF s will be expanded to the required order in ǫ. This is illustrated by the following example: << HypExp (HypExp[#, Ε, ] & ) //@ (Log[ Ε] HypergeometricF[ + Ε,, Ε, x]) Log[ x] Log[ Ε] ( ( ) + x Log[ x] Log[ x] PolyLog[, x] Ε ( + + )) x x x (HypExp[#, Ε, ] & ) //@ (Log[ Ε] HypergeometricPFQ[{ + Ε, Ε, }, { Ε, + 3 Ε}, x]) Log[ x] Log[ Ε] ( ( ) + x Log[ x] Log[ x] PolyLog[, x] Ε ( ( ) )) x x x Simplify[%] ((Log[ Ε] (Log[ x] ( + 4 Ε + Ε Log[ x]) + 4 Ε PolyLog[, x])) / ( x)) The result is not given as a SeriesData since this would have the effect of forcing the expansion of the rest of the expression. This example also illustrates that the results produced by the package are not simplified, as this might be time consuming and not always appropriate. If one wants to get a compact result, one should use Simplify or even FullSimplify. The prefactors that accompany the variable ǫ can also be symbolic, and the expansion also works for argument as shown by the following example. FullSimplify[HypExp[HypergeometricF[ + 3 Ε, Ε, 3 + Ε, ], Ε, 3]] SeriesData[Global Ε,, {,, ( 6 + Pi ), 36 ( + Zeta[3])},, 4, ] HypExp[HypergeometricPFQ[{ + Ε, Ε, 3 Ε}, { + Ε, Ε}, ], Ε, ] SeriesData[Global Ε,, 4 Pi 8 Pi 78 Zeta[3] {,,, 6 + },, 3, ]

23 The technicalities of the exansion in the case of the argument being unity are explained in section 4. The functionhypexpint[χ, χ, χ 3, χ 4, χ 5, ] evaluates integrals of the form I (χ, χ, χ 3, χ 4, χ 5, ) with w χ + χ 3 + χ 4 + δ χ5, 5, () as described in section.3. All the χ i are non-negative integers and W. The upper bound on the weight w stems from the fact that the expansion of F -functions up to order O(ǫ 4 ) also involves integrals of weight 5. The integral can be called with the argument being symbolic: HypExpInt[,,,, 3, ] (7 + Pi 4 4 Pi + Pi ) Log[ ] 3 ( + ) 3 ( ) Log[ ] Log[ ] ( + ) 3 Log[ ] Log[] ( ) PolyLog[, ] 3 3 ( + ) PolyLog[3, ] Zeta[3] 8 Zeta[3] + 4 Zeta[3] 3 4 ( + ) Arguments that match the pattern /( ) are treated with the relations between polylogarithms of different arguments, section 5. and appendix A., being already taken into account: HypExpInt[,,,,, ] (6 + Pi ) ( + ) Log[ ] ( ) ( + ) Log[ ] ( + ) Log[ ] Log[] + 6 ( + ) PolyLog[, ] + ( + ) Log[ ] PolyLog[, ] + ( + ) PolyLog[3, ] ( + ) Zeta[3] Finally, also all W can be inserted directly. For the special cases and the integral simplifies considerably and the restriction w 5 can be dropped. HypExpInt[, 4,,, 3, ] Pi Pi + 4 Zeta[3] + Pi Zeta[3] Pi Zeta[3] 8 Zeta[3] 348 Zeta[5] Pi Zeta[5] 4 Zeta[7] In the case we refer the reader also to the next paragraph and to section.3.3. The function HypExpU[n, m, p] is described in section.3.3. It evaluates integrals of the form U (n, m, p) with p and n, m being non-negative integers. () 3

24 The only additional constraint on the parameters is that in the case p < the condition m + p has to be satisfied in order to yield a convergent integral. HypExpU[4, 3, ] 6 4 Pi Pi + 44 Zeta[3] + 48 Pi Zeta[3] Pi Zeta[3] 6 Zeta[3] 576 Zeta[5] Pi Zeta[5] 5 Zeta[7] HypExpIsKnownToOrder[a,...,a J,b,...,b J,n] returnstrue if the expansion of the hypergeometric function with parameters corresponding to the first arguments of the function is available in the library to the order n. HypExpAddToLib adds an expansion to the library, this method is described in one of the next sections. 5. Functions modified The package also updates Series to allow it to expand compound expressions containing hypergeometric functions. The difference between this and the mapping with HypExp is that the other functions of ǫ are also expanded, as shown by the following example: Series[Log[ Ε] HypergeometricF[Ε +,, Ε, x], {Ε,, }] SeriesData[Global Ε,, Log[ x] Log[ x] Log[ x] {, ( ) x x x PolyLog[, x] },, 3, ] x Series[Log[ Ε] HypergeometricPFQ[{ Ε +, Ε, }, { Ε, 3 Ε + }, x], {Ε,, }] SeriesData[Global Ε,, Log[ x] 5 Log[ x] Log[ x] {, + + x x x PolyLog[, x] },, 3, ] x This allows to work with the expansion of HF s as Mathematica users are used to. We also updated the series expansion of the regularied hypergeometric functions since they are nothing else but hypergeometric functions divided by Γ-functions. Since the incomplete B function is a special case of HF, B(, a, b) a a F (a, b, a +, ), a,,... (3) it is also possible to expand it with the method described in this paper. Therefore we also updated the series expansion of the incomplete B function around integer values of its parameters, as shown by the following example Series[Beta[x, + Ε, Ε], {Ε,, }] SeriesData[Global Ε,, x 3 x x Log[ x] {x, + + x Log[ x] + 4 x Log[ x] + x Log[x] x Log[x]},,, ] In order to account for the relations between polylogarithms Li n and Nielsen polylogarithms S n,p of different arguments we updated PolyLog and implemented the relations (3) (36) given in appendix A.. This is also illustrated with an example. PolyLog, Log PolyLog, 4

25 5.3 Working with the libraries Since the computation of the expansion for high orders and large parameters is quite time consuming, it is of interest to store the results that have been already calculated and reuse them, instead of recalculating them. The HypExp library contains expansions of HF for some sets of parameters. When an expansion is requested, the package first checks whether the library contains the expansion for the requested set of parameters to the requested order. If so, it loads it and gives the answer, if not proceeds with the calculation along the line of the preceding sections. The library management can be called dynamic in the sense that elements of the library are loaded in the memory at run time only if they are needed. The package HypExp has a standard library. Further libraries can be added, depending on the needs and on the amount of available disk space. The different libraries can be found at [8]. The package also allows to extend the provided libraries with HF s not included in the standard libraries, or included but not to a sufficiently high expansion order. The expansion of J F J (a,..., a J, b,...,b J, x) to the order n can be added to the library with the command HypExpAddToLib[a,...,a J,b,...,b J,n] where a,..., a j, b,...,b J are integers. Upon this command, Mathematica computes the expansion for arbitrary ǫ-parts added to the integers a,..., a j and b,..., b J. As arbitrary coefficients are more difficult than numbers, the time needed to add an expansion to the library is longer than the time for a single evaluation. Therefore adding expansions to the library is only useful if this expansion shows up repeatedly. The result is then saved in the library in the installation directory, so that it can be used in future Mathematica sessions. Since the results are supposed to be used several times, the result of the expansion is simplified using Simplify, in order to get a more compact result. This, in turn, makes the extension of the library even more time consuming. The library files are copied to the installation directory of the package during the installation. Further libraries can be added as descriebed at [8]. Due to our naming conventions for the entries in the library files, the expansions of HF s of parameters higher in absolute value than 9 are not allowed. This might be fixed if such HF s are required. 5.4 Note on the expansion for negative parameters Let us consider F ( m+α, b, m l+β; x) for m, l, b being positive integers, b > m and α, β small. Using the definition one gets F ( m + α, b, m l + β; x) n ( m + α) n ( b) n ( l m + β) n n! xn. (4) We are interested in the terms for n between m and b. All further terms vanish, since then ( b) n becomes. They are equal to ( m + α)...( + α)α( b)( b)...(m b) ( m + β l)( m + β + l)...( + β l)(β l) xm + + ( m + α)...α( + α)...(m b + α)( b)...( b + m + l) x m+l +... ( m + β l)...( + β)β (5) If one wants to define a value for F ( m, b, m l; x) one has to take the limit of the above expression for α and β going to ero. The result depends on which is stored in the Mathematica variable $HypExpPath 5

26 the way one approches (, ) with α and β. In [7], one can find the formula F ( m, b, m l; x) n ( m) n (b) n ( l m) n n! xn (6) which is also the result Mathematica gives. This corresponds to a trajectory in the (α, β)-plane going along the β axis. Taking a trajectory along the α axis would lead to a /β pole. Any other trajectory gives a constant times a function. This function happens to be the second solution of the differential equation satisfied by F ( m + α, b, m l + β; x), x( x)w (x) + (B (A + A + )x) w (x) A A w(x) (7) with A m + α, A b, B m l + β. (8) In the case of negative or vanishing c, the value at x of this function is also equal to unity. This prevents us from discriminating the two solutions by their values at x. Since the use of the Kummer identity (8) F (A, A ; B ; x) ( x) B A A F (B A, B A ; B ; x) (9) induces a rotation in the (α, β) plane and since Mathematica always chooses the trajectory along the β axis, the result for the HF and its Kummer transform will not be identical in Mathematica for this particular case. 5.5 Performances and limitations The limits are set by the CPU and the amount of memory available. In all practical cases known to the authors, however, the result is given in a reasonable time. The following table shows the CPU time dependence for the expansion of some hypergeometric functions to different orders on a 3 GH processor/.5 GB RAM machine F ( + ǫ, ǫ; + ǫ, x) < s < s < s 7 s F ( + α ǫ, + α ǫ; + β ǫ, x) < s < s < s 6 s 3F ( + ǫ, ǫ, 3ǫ; + 3ǫ, + ǫ, x) < s < s < s 3 s 3F ( + α ǫ, + α ǫ, + α 3 ǫ; + β ǫ, + β ǫ, x) < s < s.5 s 3 s 4F 3 ( + α ǫ, + α ǫ, 3 + α 3 ǫ, 4 + α 4 ǫ; β ǫ, + β ǫ, + β 3 ǫ, x) s s 5 s 4 s This package was developed in Mathematica 5. and should work on newer versions. 6 Summary In this paper we presented the Mathematica package HypExp for expanding arbitrary hypergeometric functions to arbitrary order in a small quantity around integer-valued parameters. These expansions are required for example in the computation of multi-loop or multi-particle phase space integrals in dimensionally regularied quantum field theory. A first application was presented in ref. [6], namely the ǫ-expansion of the two-loop quark and gluon form factors whose exact analytic expressions contain hypergeometric functions in addition to Γ-functions. 6

27 Acknowledgement We would like to thank Thomas Gehrmann and Alejandro Daleo for useful discussions, a careful reading of our manuscript, and for patiently testing our early versions. We also wish to thank Gudrun Heinrich for independent numerical checks of the series expansions using the sector decomposition method described in [9 3]. This research was supported by the Schweierischer Nationalfonds under contract A Useful relations In this appendix we collect useful relations among logarithms, polylogarithms Li n, and Nielsen polylogarithms S n,p as well as some additional integrals. The package updates PolyLog as described in section 5. in order to allow for the implementation of those relations that relate polylogarithms and Nielsen polylogarithms, i.e. Eq. (3) (36). The relations are based on [9,6] and hold at least for all W, where W is defined in Eq. (5). A. Relations between logarithms and polylogarithms ln( ) ln( ) ln( ) () ln( ) ln( ) () ln( ) ln( ) + ln() () Li ( ) Li () + π ln() ln( ) 6 (3) Li ( ) Li () ln ( ) (4) Li ( ) Li () ln ( ) + π + ln( ) ln( ) 6 (5) Li ( ) ln ( ) π 6 Li () (6) Li ( ) ln ( ) + π 3 ln( ) ln( ) + Li () (7) Li 3 ( ) Li 3() Li 3 ( ) + ζ(3) + π 6 ln( ) ln() ln ( ) + 6 log3 ( ) (8) Li 3 ( ) 6 ln3 ( ) ln( ) ln ( ) + ln() ln ( ) Li 3 ( ) Li 3 ( π 3 ln( ) + Li 3( ) (9) Li 3() π 6 ln( ) 6 ln3 ( ) () ) Li 3 () Li 3 ( ) + ζ(3) 6 ln3 ( ) π 6 ln( ) + 6 ln3 ( ) + π 6 ln( ) ln ( )ln() () 7

28 Li 4 ( ) 4 ln4 ( ) + 6 ln( ) ln3 ( ) 6 ln() ln3 ( ) + π 6 ln ( ) + π4 45 Li 4( ) () Li 4 ( ) Li 4 ( ) 4 ln4 ( ) π ln ( ) 7π4 36 (3) Li 5 ( ) ln5 ( ) 4 ln( ) ln4 ( ) + 4 ln() ln4 ( ) π 8 ln3 ( ) π4 45 ln( ) + Li 5( ) (4) S, () 4 ln4 ( ) 6 ln() ln3 ( ) + π ln ( ) Li 3 () ln( ) + ζ(3) ln( ) Li 4 ( ) + Li 4 () +Li 4 ( ) + π4 9 S, ( ) 4 ln4 ( ) + 6 ln() ln3 ( ) 4 ln () ln ( ) π ln ( ) + π 6 ln() ln( ) ln()li 3( ) +Li 4 ( ) Li 4 () Li 4 ( (5) ) + ζ(3) ln() π4 (6) S, ( ) ln4 ( ) 3 ln() ln3 ( ) + π ln ( ) π4 9 Li 3 ( ) ln( ) Li 3 () ln( ) +Li 4 ( ) + Li 4 () + Li 4 ( ) (7) S, ( ) 3 ln() ln3 ( ) 6 ln( ) ln3 ( ) π ln ( ) ln( )ln()ln ( ) + 4 ln ( )ln ( ) Li 4 () Li 4 ( ) Li 4 ( ) + Li 3( )ln( ) + π π4 ln( )ln( ) + ζ(3)ln( ) 6 7 ln( )Li 3 ( ) + ζ(3)ln( ) (8) S,3 () 8 ln() ln4 ( ) 3 ln5 ( ) π 8 ln3 ( ) + Li 3()ln ( ) + [ Li 4 ( ) Li 4 ( )] ln( ) ζ(3) ln ( ) Li 5 () Li 5 ( ) Li 5 ( ) +S 3, () + ζ(5) (9) 8