On logarithmic integrals, harmonic sums and variations
Ming Hao Zhao
Abstract Based on various non-MZV approaches we evaluate certain logarithmic integrals and har- monic sums. More specifically, 85 LIs, 89 ESs, 263 PLIs, 28 non-alt QESs, 39 GESs, 26 BESs, 14 IBESs with weight ≤ 5, 193 QLIs, 172 QPLIs, 83 QESs, 7 QBESs, 3 IQBESs with weight ≤ 4. With help of MZV theory we evaluate other 10 QESs with weight ≤ 4. High weight hypergeometric series, nonhomogeneous integrals and other results are also established.
Contents
1 Preliminaries 5 1.1 Special functions ...... 5 1.2 Definition of subjects ...... 6 1.2.1 LI and PLI ...... 6 1.2.2 QLI and QPLI ...... 7 1.2.3 ES and QES ...... 8 1.2.4 Fibonacci basis ...... 9 1.3 Main results ...... 11
2 LI 11 2.1 Brute force ...... 11 2.2 General formulas ...... 12 2.3 Integration by parts ...... 14 2.4 Fractional transformation ...... 14 2.5 Beta derivatives ...... 15 2.6 Contour integration ...... 16 arXiv:1911.12155v13 [math.CA] 1 Jul 2020 2.7 Double integration ...... 16 2.8 Hypergeometric identities ...... 16 2.9 Solution ...... 17
1 3 ES 17 3.1 General formulas ...... 17 3.2 Symmetric relations ...... 18 3.3 Partial fractions ...... 18 3.4 Contour integration ...... 18 3.5 Solution ...... 18 3.5.1 Preparation ...... 18 3.5.2 Confrontation ...... 24
4 PLI 27 4.1 Brute force ...... 27 4.2 General formulas ...... 27 4.3 Integration by parts ...... 28 4.4 Fractional transformation ...... 28 4.5 Complementary integrals ...... 28 4.6 Polylog identities ...... 28 4.7 Series expansion ...... 30 4.8 Solution ...... 32
5 QLI 33 5.1 Logsine integrals and contour deformation ...... 33 5.1.1 LSI, PLSI, ITF ...... 33 5.1.2 QLSI and QLI ...... 35 5.2 Inverse tangent family ...... 35 5.3 Polylog special values ...... 35 5.4 Brute force ...... 37 5.4.1 Class 1 ...... 37 5.4.2 Class 2 ...... 37 5.5 General formulas ...... 37 5.6 Feynman’s trick on Mellin transform ...... 38 5.7 Integration by parts ...... 39 5.8 Fractional transformation ...... 39 5.9 Power substitution ...... 39 5.10 Weierstrass substitution ...... 39 5.11 Parametric integration ...... 39 5.12 Fourier expansion ...... 40 5.13 Beta derivatives ...... 40 5.13.1 Class 1 ...... 40 5.13.2 Class 2 ...... 40
2 5.14 Contour integration ...... 40 5.14.1 Class 1 ...... 40 5.14.2 Class 2 ...... 41 5.15 Series expansion ...... 41 5.16 Solution ...... 42
6 QPLI 42 6.1 Brute force ...... 42 6.2 General formulas ...... 42 6.3 Integration by parts ...... 43 6.4 Fractional transformation ...... 43 6.5 Polylog identities ...... 43 6.6 Power substitution ...... 43 6.7 Valean’s master formula ...... 43 6.7.1 Class 1 ...... 43 6.7.2 Class 2 ...... 44 6.8 Contour integration ...... 44 6.8.1 Class 1 ...... 44 6.8.2 Class 2 ...... 44 6.9 Series expansion ...... 45 6.9.1 Class 1 ...... 45 6.9.2 Class 2 ...... 46 6.10 Double integration ...... 46 6.10.1 Class 1 ...... 46 6.10.2 Class 2 ...... 47 6.10.3 Class 3 ...... 48 6.10.4 Class 4 ...... 48 6.10.5 Class 5 ...... 49 6.11 Solution ...... 50
7 QES 52 7.1 Non-alternating case ...... 52 7.2 Generating functions ...... 52 7.3 QLI/QPLI expansion ...... 53 7.4 Integral representation ...... 53 7.5 Contour integration ...... 54 7.6 Partial fractions ...... 54 7.7 Symmetric relations ...... 54 7.8 Solution ...... 54
3 8 Applications 55 8.1 Nonhomogeneous integrals ...... 55 8.1.1 NLI ...... 55 8.1.2 NPLI, NQLI, NQPLI ...... 55 8.2 Nested sums and multiple log integrals ...... 56 8.2.1 NS manipulation ...... 56 8.2.2 MLI representation ...... 57 8.3 Geometric and binomial Euler sums ...... 57 8.3.1 GES ...... 57 8.3.2 BES ...... 58 8.3.3 IBES ...... 61
8.3.4 QBES, IQBES (pFq part 1) ...... 62 8.4 Related problems ...... 68 8.4.1 Symmetric PLI ...... 68 8.4.2 MT of PLI ...... 69 8.4.3 Inverse ES ...... 69 8.4.4 High weight ES ...... 70 8.4.5 High weight PLI ...... 71 8.4.6 Parseval LSI ...... 71 8.4.7 Convolutional LI ...... 72 8.4.8 Valean’s hexagram ...... 73
8.4.9 Hypergeometric rhapsody (pFq part 2) ...... 75
8.4.10 FL fantasia (pFq part 3) ...... 77 8.5 Two conjectures solved ...... 80
8.5.1 Borwein’s conjecture on binomial sum (pFq part 4) ...... 80 8.5.2 Au’s conjecture on general LI ...... 82
9 Table of subjects 84 9.1 LI ...... 84 9.2 ES...... 87 9.3 PLI ...... 90 9.4 QLI ...... 99 9.5 QPLI ...... 106 9.6 QES...... 111 9.7 GES...... 115 9.8 BES ...... 117 9.9 IBES ...... 118 9.10 QBES, IQBES ...... 118
4 10 Miscellaneous 119 10.1 LSI ...... 119 10.2 PLSI ...... 120 10.3 ITF ...... 120 10.4 QLI (W =5) ...... 121 10.5 QLSI ...... 121 10.6 NLI ...... 122 10.7 NPLI ...... 122 10.8 NQLI ...... 123 10.9 NQPLI ...... 124 10.10NS, MLI ...... 124 10.11Miscellaneous ...... 125
11 Comments 126
1 Preliminaries
1.1 Special functions
Below are special functions used. We assume readers are aware of their basic properties [13, 14].
R ∞ s−1 −t Gamma function: The analytic continuation of Γ(s) = 0 t e dt
R 1 a−1 b−1 Γ(a)Γ(b) Euler Beta function: The analytic continuation of B(a, b) = 0 t (1 − t) dt = Γ(a+b)
n+1 (n) ∂ log(Γ(x)) Polygamma function: The analytic continuation of ψ (x) = ∂xn+1
P∞ zk Polylogarithm: The analytic continuation of Lin(z) = k=1 kn
P∞ 1 Riemann Zeta function: The analytic continuation of ζ(s) = k=1 ks
k−1 1−s P∞ (−1) Dirichlet Eta function: The analytic continuation of η(s) = (1 − 2 )ζ(s) = k=1 ks
P∞ (−1)k Dirichlet Beta function: The analytic continuation of β(s) = k=0 (2k+1)s
(3) 1 (3) 3 ψ ( 4 )−ψ ( 4 ) Catalan constant and generalization: C = β(2), 1536 = β(4)
(0) Γ(a+n) Euler-Mascheroni constant: γ = −ψ (1) Pochhammer symbol: (a)n = Γ(a)
5 Generalized hypergeometric function: The analytic continuation of
∞ n X (a1)n ··· (ap)nz F (a , ··· , a ; b , ··· , b ; z) = p q 1 p 1 q (b ) ··· (b ) n! n=0 1 n q n
P∞ sin(nx) P∞ cos(nx) Clausen function: Cl2k(x) = n=1 n2k , Cl2k+1(x) = n=1 n2k+1
R x Tik−1(t) −1 Inverse tangent integral: Tik(x) = 0 t dt, Ti1(x) = tan (x)
π π q Elliptic K/E integral: K(x) = R 2 √ 1 dt, E(x) = R 2 1 − x sin2(t) dt 0 1−x sin2(t) 0
n 2 n 1 ∂ (x −1) Legendre P polynomial: Pn(x) = (2n)!! ∂xn
1.2 Definition of subjects
Clarification: Throughout this paper, LI will be the abbreviation of log integral, PLI polylog inte- gral, QLI quadratic log integral, QPLI quadratic polylog integral, ES Euler sum, QES quadratic Euler sum, MZV multiple Zeta value. Minor abbreviations will be mentioned later.
1.2.1 LI and PLI
Notations:
f(0; x) = 1 − x, f(1; x) = x, f(2; x) = 1 + x
l(0; x) = log(1 − x), l(1; x) = log(x), l(2; x) = log(1 + x), l(3; x) = log(1 + x2), l(4; x) = tan−1(x)
1 − x L(n, 1; x) = Li (x),L(n, 2; x) = Li (−x),L(n, 3; x) = Li ( ), n ≥ 2 n n n 2
1 + x 1 − x x − 1 L(n, 4; x) = Li ( ),L(n, 5; x) = Li ( ),L(n, 6; x) = Li ( ), n ≥ 2 n 2 n 1 + x n 1 + x
Definition:
Z 1 Q2 l(m; x)a(m) LI(a(0), a(1), a(2); p) = m=0 dx 0 f(p; x)
6 P LI(a(0), a(1), a(2); b(2, 1), ··· , b(2, 6); ··· ; b(N, 1), ··· , b(N, 6); p) Z 1 Q2 l(m; x)a(m) Q6 QN L(n, k; x)b(n,k) = m=0 k=1 n=2 dx 0 f(p; x)
Abbreviation: In later discussion of PLIs, N and b(n, k) will be small and many of them equal to 0. In this case, we omit those zero parameters and simply repeat the word ’nk’ b(n, k) times for abbreviation. For instance we have
Z 1 log(1 − x) log2(x) log3(x + 1) LI(1, 2, 3; 0) = dx 0 1 − x 2 Z 1 Li x+1 Li (−x) log(x + 1) P LI(0, 0, 1; 24, 24; 32; 2) = 2 2 3 dx 0 x + 1
Weight: We define the weight W of LI, PLI as
2 6 N 2 X X X X 1 + a(m), 1 + nb(n, k) + a(m) m=0 k=1 n=2 m=0 respectively, so LI(1, 2, 3; 0) is of weight 7 and P LI(0, 0, 1; 24, 24; 32; 2) weight 9.
For W ≤ 5 there are 85 convergent LIs, since a simple counting shows the total number of 1 2 weight W convergent LIs is 2 3W + W − 2 . For W ≤ 5, after excluding those already contained in LIs, there are 263 convergent PLIs.
1.2.2 QLI and QPLI
Notations: 1 + x2 g(1; x) = 1 − x, g(2; x) = x, g(3; x) = 1 + x, g(4; x) = 1 + x2, g(5; x) = x
k(1; x) = log(1 − x), k(2; x) = log(x), k(3; x) = log(1 + x), k(4; x) = log(1 + x2), k(5; x) = tan−1(x)
2 L(n, 7; x) = Lin(−x ),L(n, 8; x) = Tin(x)
7 Definition:
Z 1 Q5 k(m; x)A(m) QLI(A(1),A(2),A(3),A(4),A(5); p) = m=1 dx 0 g(p; x)
Z 1 L(2, a; x) Z 1 L(3, a; x) Z 1 L(2, a; x)k(b; x) QP LI3(a; p) = dx, QP LI4(a; p) = dx, QP LI4(a; b; p) = dx 0 g(p; x) 0 g(p; x) 0 g(p; x)
Abbreviation: In later discussion of QLIs, A(m) will be small and many of them equal to 0. In this case, we omit those zero parameters and simply repeat the word ’m’ A(m) times for abbreviation. For instance we have
Z 1 x log(x + 1) tan−1(x)2 Z 1 log(1 − x) log x2 + 1 tan−1(x) QLI(355; 5) = 2 dx, QLI(145; 4) = 2 dx 0 x + 1 0 x + 1
Z 1 Li (−x2) Z 1 Ti (x) log(x) QP LI4(7; 3) = 3 dx, QP LI4(8; 2; 1) = 2 dx 0 x + 1 0 1 − x
P5 Weight: We define the weight W of QLI as 1 + m=1 A(m). QP LI3(··· ), QP LI4(··· ) are of weight 3, 4 respectively. All 4 examples above are of weight 4.
For W ≤ 4, after excluding those already contained in LIs, there are 193 convergent QLIs. For W ≤ 4, after excluding those already contained in LI/PLI/QLIs, there are 172 convergent QPLIs.
1.2.3 ES and QES
Notations:
n n j−1 (k) X 1 ](k) X (−1) (1) (1) H = , Hn = ,Hn = H , Hfn = Hgn n jk jk n j=1 j=1
(k) ](k) (k) f1k(n) = Hn , f2k(n) = Hn , f3k(n) = H2n
](k) (k) 1 (k) Note that H2n won’t appear since it equals H2n − 2k−1 Hn .
Definition:
8 ∞ n−1 M N X (±1) Y Y (b(j)) ES(a(1), ··· , a(M), −b(1), ··· , −b(N); ±p) = H(a(k)) H^ np n n n=1 k=1 j=1 ∞ s X (±1)n−1 Y QES1(m k , ··· , m k ; ±p) = f (n) 1 1 s s np mr kr n=1 r=1 ∞ s X (±1)n−1 Y QES2(m k , ··· , m k ; ±p) = f (n) 1 1 s s (2n + 1)p mr kr n=1 r=1 (k) We may assume a(i), b(j), mrkr are non-decreasing. If there exists no Hn (resp. f2kr (n)) and the outer index +p, we call an ES (resp. QES) non-alternating, otherwise alternating.
PM PN Ps Weight: We define the weight W of ES as k=1 a(k) + j=1 b(j) + p and QES r=1 kr + p. For instance 2 n−1 (2) g(3) (2) (2) ∞ (−1) Hn Hn ∞ H H 2H X X n n 2n ES(2, −3, −3; −4) = , QES(11, 12, 12, 32; 3) = n4 n3 n=1 n=1 are weight 12 alternating ES and weight 10 non-alternating QES, respectively.
For W ≤ 5 there are 89 convergent ESs. For W ≤ 4, after excluding those already contained in ESs, there are 93 convergent QESs. For W = 5, after excluding those already contained in ESs, there are 28 convergent non-alternating QESs.
1.2.4 Fibonacci basis
This part is about the structure of closed-forms.
Weight of constants: We assign weight 0 to rational numbers, 1 to π, log(2), 2 to Catalan C, 4 to (3) 1 (3) 3 1 ψ ( 4 ) − ψ ( 4 ), n to ζ(n), Lin( 2 ),
Linear relations: A linear relation is an equality of form
’rational combination of LI/PLI/QLI/QPLI/ES=certain constants’
If all items in LHS have the same weight, we call it a homogeneous relation; otherwise we call it non-homogeneous.
9 Fibonacci basis: According to [1], all 85 LIs with weight W ≤ 5 are generated by certain t(W ) weight W constants over Q, here t(1) = 1, t(2) = 2, t(3) = 3, t(4) = 5, t(5) = 8. We define the sets of these homogeneous constants generating LIs as Fibonacci bases due to consistency between t(n) and
Fibonacci numbers Fn. This term arises from [6] but differs a little. See also [1] for higher weight LIs.
We denote AW the weight W Fibonacci basis containing t(W ) full-simplified polylog constants. 1 By full-simplified we mean ζ(2n), n > 0 and Lin 2 , n < 4 won’t appear, due to simplification formulas [13, 14] 1 π2 log2(2) 1 7ζ(3) log3(2) 1 22n−1π2n |B | Li = − , Li = + − π2 log(2), ζ(2n) = 2n 2 2 12 2 3 2 8 6 12 (2n)! where Bn are Bernoulli numbers. These simplifications keeps the weight unchanged. We have [1]
2 2 2 3 A1 = {log(2)} ,A2 = π , log (2) ,A3 = π log(2), log (2), ζ(3) 1 A = π4, π2 log2(2), log4(2), ζ(3) log(2), Li 4 4 2 1 1 A = π4 log(2), π2 log3(2), log5(2), π2ζ(3), ζ(3) log2(2), ζ(5), Li log(2), Li 5 4 2 5 2
And that all LIs with W ≤ 5 are generated by AW over Q. A possible weight 6 representation is:
6 4 2 2 4 6 2 3 A6 = π , π log (2), π log (2), log (2), π ζ(3) log(2), ζ(3) log (2), ζ(5) log(2), 1 1 1 1 ζ(3)2, Li log2(2), π2Li , Li log(2), Li ,S 4 2 4 2 5 2 6 2
P∞ (−1)k−1 Where the constant S = ES(1; −5) = k=1 k5 Hk is very likely to be irreducible. Higher weights’ representations are omitted.
Extended basis: We define extended Fibonacci bases for W ≤ 4 as
B2 = A2 ∪ {C, π log(2)}
2 3 B3 = A3 ∪ C log(2), πC, π log (2), π , =(Li3(1 + i))
2 2 2 3 3 B4 = A4 ∪ C , πC log(2),C log (2), π C, π log (2), π log(2), log(2)=(Li3(1 + i)), 1 3 πζ(3), π=(Li (1 + i)), =(Li (1 + i)), ψ(3) − ψ(3) 3 4 4 4
W Note that in these 3 cases each BW contains 2 constants.
10 1.3 Main results
In section 3 ∼ 6 the following will be proved:
Theorem 1. All 89 ESs and 263 PLIs with weight W ≤ 5 are generated by AW over Q.
Theorem 2. All 193 QLIs and 172 QPLIs with weight W ≤ 4 are generated by BW over Q.
Theorem 1, a natural consequence of MZV theory (see for instance [18]), will be proved by a systematic non-MZV approach. Theorem 2, intimately connected with level 4 MZVs, will be established by generalizing the approach above.
See tables of section 9 for all their closed-forms and 8, 10 for applications.
2 LI
LIs with weight W ≤ 5 are completely and explicitly solved in [1]. We summarize and add some new results, generating linear relations between LIs.
2.1 Brute force R log(a+x) log(b+x) R log2(a+x) log(b+x) R log(a+x) log2(b+x) R logn(a+x) It is known that c+x dx, a+x dx, a+x dx, b+x have polylog primitives. By Mathematica we have:
Z log3(a + x) a + x a + x F1: dx = 6Li + 3 log2(a + x)Li b + x 4 a − b 2 a − b a + x b + x −6 log(a + x)Li + log3(a + x) log − 3 a − b a − b
Z 3 log2(a + x) log(b + x) a + x a + x F2: dx = −6Li − 3 log2(a + x)Li a + x 4 a − b 2 a − b a + x b + x +6 log(a + x)Li − log3(a + x) log − + log3(a + x) log(b + x) 3 a − b a − b
Z log(a + x) log(b + x) (a − c)(b + x) b + x (a − c)(b + x) (a − c)(b + x) F3: dx = Li3 + Li2 − Li2 log c + x (b − c)(a + x) a + x (b − c)(a + x) (a + x)(b − c)
11 a + x (a − c)(b + x) b + x (a − c)(b + x) +Li2 log(b + x) − log + Li2 log + log(a + x) a − c (a + x)(b − c) b − c (a + x)(b − c) 1 (a − b)(c + x) a − b c + x 2 (a − c)(b + x) b + x + − log + log + log log − Li3 2 (a + x)(c − b) a + x c − a (a + x)(b − c) a + x
a + x b + x b + x c + x c + x (a − c)(b + x) −Li3 − Li3 + log log − − log log a − c b − c b − c b − c c − a (a + x)(b − c) 1 b + x b + x c + x c + x c + x + log log − 2 log(a + x) log − log − + log(a + x) log(b + x) log 2 b − c b − c c − a b − c c − a
2 ! ! !! ! ! Z 12 log (a + x) log(b + x) a + x b + x b + x b + x 2 b + x F4: dx = −24 −Li4 + Li4 − + Li4 − 12Li2 log − b + x a − b a − b a + x a + x a + x
a + x ! b + x ! b + x !!2 a + x ! a + x ! b + x !! +12Li2 log − − log − + 24Li2 log log(b + x) − log − a − b a − b a + x a − b a − b a − b
b + x ! a + x ! b + x !! b + x ! b + x ! b + x !! ! +12Li2 − 2 log log(b + x) − log − + log − log − − 2 log − − 2 log(a + x) log(b + x) a − b a − b a − b a − b a − b a + x
a + x ! b + x ! ! b + x ! b + x ! ! b + x ! b + x ! +24Li3 log − − log(b + x) + 24Li3 − log − + log(a + x) + 24Li3 log − a − b a + x a − b a + x a + x a + x
! ! ! !! ! ! ! 4 b + x 4 b + x b + x a − b 3 b + x 2 b + x 2 b + x + log − + log − − 4 log − + log log − + 6 log − log − a − b a + x a − b a + x a + x a + x a − b
! !! ! ! a + x b + x b + x 2 b + x 2 2 −4 3 log + log − log − log − + 6 log (a + x) log (b + x) a − b a − b a + x a − b
! !! ! ! ! !! 2 a + x b + x b + x a + x 2 3 b + x +6 log log(b + x) − log − 3 log − + log(b + x) − 4 log 3 log(a + x) log (b + x) − 2 log − a − b a − b a − b a − b a − b
Here logs and polylogs are principal, following the convention of Mathematica. Based on this, for each LI we may check by Mathematica whether its primitive is computable. If it is, we get the closed-form directly via Newton-Leibniz formula.
2.2 General formulas
Proposition 1. The following formulas hold:
Z 1 (−1)nn! m n (1) x log (x) dx = n+1 0 (m + 1)
logn+1(2) LI(0, 0, n; 2) = (2) n + 1
LI(0, n, 0; 0) = LI(n, 0, 0; 1) = (−1)nn!ζ(n + 1) (3)
LI(0, n, 0; 2) = 1 − 2−n (−1)nn!ζ(n + 1) (4)
12 1 LI(n, 0, 0; 2) = (−1)nn!Li (5) n+1 2
LI(1, n, 0; 1) = LI(n, 1, 0, 0) = (−1)n−1n!ζ(n + 2) (6)
LI(0, n, 1; 1) = 1 − 2−(n+1) (−1)nn!ζ(n + 2) (7)
n−1 k 1 n+1 X log (2)Li−k+n+1 n log (2) LI(0, 0, n; 1) = −n! 2 + n!ζ(n + 1) − (8) k! n + 1 k=0
n k 1 n+2 X log (2)Li−k+n+2 log (2) LI(0, 1, n; 2) = n! 2 − n!ζ(n + 2) + (9) k! n + 2 k=0
n+2 X (−1)k−1ζ(k) log−k+n+2(2) 1 logn+2(2) LI(1, 0, n; 2) = n! + (−1)nn!Li + (10) (−k + n + 2)! n+2 2 n + 1 k=2
n ! 1 X LI(1, n, 0; 0) = LI(n, 1, 0; 1) = (−1)n−1n! (n + 1)ζ(n + 2) − ζ(k)ζ(−k + n + 2) (11) 2 k=2
n 1 X LI(0, n, 1; 0) = (−1)nn! 1 − 21−k 1 − 2k−n−1 ζ(k)ζ(−k + n + 2) 2 k=2 (12) 1 −(−1)nn! (n + 1)ζ(n + 2) − 2 1 − 2−(n+1) log(2)ζ(n + 1) 2
n−1 X LI(0, 2n − 1, 1; 2) = −(2n − 1)! 1 − 21−2k ζ(2k)ζ(−2k + 2n + 1) k=1 (13) −(2n − 1)! 2−(2n+1)ζ(2n + 1) − (n − 1) 1 − 2−2n ζ(2n + 1)
n−1 X LI(1, 2n − 1, 0; 2) = −(2n − 1)! 1 − 22k−2n ζ(2k)ζ(−2k + 2n + 1) k=1 (14) +(2n − 1)! (n + 1) 1 − 2−2n + 2−(2n+1) ζ(2n + 1) − 2 1 − 2−2n log(2)ζ(2n)
n+1 ! X LI(2, n, 0; 1) = (−1)nn! (n + 2)ζ(n + 3) − ζ(k)ζ(−k + n + 3) (15) k=2
13 n X LI(0, 2n, 2; 1) = 2(2n)! 1 − 21−2k ζ(2k)ζ(−2k + 2n + 3) k=1 (16) −2(2n)! 1 − 2−(2n+2) nζ(2n + 3) − 2−(2n+3)ζ(2n + 3)
n X LI(1, 2n, 1; 1) = −(2n)!2−(2n+1) 1 − 2−2k+2n+2 ζ(2k)ζ(−2k + 2n + 3) k=1 (17) +(2n)! (2n + 3)2−(2n+3) − 1 ζ(2n + 3)
−t 1 Proof. (1): x → e , Γ. (2): Trivial. (3)-(5): Series expansion of 1−ax , then use (1). (6)- 1 (7): Series expansion of log(1 − ax). (8): Nielsen-Ramanujan integral. t → 1+x , then use R 1 m n ∂n 1−2−(m+1) 1 x log (x) dx = n . (9)-(10): Use (8), integrate by parts. (11)-(12): Series ex- 2 ∂m m+1 log(1−x) pansion of 1−x , use (1) then ES results in [11]. (13)-(17) Series expansion and [11].
For each LI, check whether they are directly solved using general formulas above.
Clarification: Throughout this paper, legitimacy of all operations of exchanging orders of mul- tiple summation and integration can be justified by convergence theorems.
2.3 Integration by parts
Q2 a(m) R 1 0 Take F (x) = m=0 l(m; x) in N-L formula 0 F (x) dx = F (1) − F (0) to obtain linear a linear relation between LIs. Now consider all possible F . As one can see, this is essentially the process of partial integration.
Clarification: Though we don’t carry out partial integration from the view of manipulating single LIs here, this will be used frequently later. Throughout this paper, by ’integration by parts lifting (up) u0 in integral R u0v, we mean transforming it into uv − R uv0.
2.4 Fractional transformation Q2 a(m) 2F 1−x m=0 l(m;x) R 1 R 1 ( x+1 ) Take F (x) = f(p;x) (p 6= 1) in the identity 0 F (x) dx = 0 (x+1)2 dx. Since log(x) → log(1−x)−log(x+1), log(1−x) → log(x)−log(x+1)+log(2), log(x+1) → log(2)−log(x+1)
1−x under the substitution x → x+1 , we may obtain a linear relation after expansion and possible 1 partial fraction decomposition on x(x+1) for p = 0. Now consider all possible F .
14 2.5 Beta derivatives
Proposition 2. The following formulas hold:
∂m+nB(a, b) LI(n, m, 0; 1) = LI(m, n, 0; 0) = lim (1) {a,b}→{0,1} ∂am ∂bn
n X n ∂m+nB(a, b) LI(k, m, n − k; 1) = 2−(m+1) lim (2) k {a,b}→{0,1} ∂am ∂bn k=0
n m+n m+n ! X n −(m+1) ∂ B(a, b) ∂ B(a, b) LI(k, m, n − k; 0) = 2 lim m n + lim m n (3) k {a,b}→{1,0} ∂a ∂b {a,b}→ 1 ,0 ∂a ∂b k=0 { 2 }
n m+n m+n ! X n −(m+1) ∂ B(a, b) ∂ B(a, b) LI(k, m, n − k; 2) = 2 lim m n − lim m n (4) k {a,b}→{ 1 ,0} ∂a ∂b {a,b}→{1,0} ∂a ∂b k=0 2
∂m+1B(a, b) ∂m+1B(a, b) LI(0, m, 1; 0) = 2−(m+1) − 1 lim + 2−(m+1) lim (5) {a,b}→{1,0} ∂am ∂b1 1 ∂am ∂b1 {a,b}→{ 2 ,0}
n 1−s X n ∂n 2 − 1 π (−1)k LI(k, n − k, 0; 2) = lim (6) k s→1 ∂sn sin(πs) k=0
n n+1 X X nn + 1 1 ∂2n+1B(a, b) log2n+2(2) (− log(2))j+kLI(−k + n + 1, 0, n − j; 2) = lim − (7) j k 2 {x,y}→{0,1} ∂an ∂bn+1 2n + 2 j=0 k=0
n n X n X n 1 (−1)k LI(0, n − k, k; 2) = (−1)n l! Li logn−l(2) (8) k l l+1 2 k=0 l=0
m n X m X n (−1)k (−1)j LI(0, j + k, −j + m + n; 1) + (−1)j LI(0, j + k, −j + m + n; 1) j j j=0 j=0 (9) k X k ∂k+m+nB(a, b) ∂k+m+nB(a, b) = (−1)m+n (−1)k lim + lim j {x,y}→{0,1} ∂aj+m ∂b−j+k+n {x,y}→{1,0} ∂aj+m ∂b−j+k+n j=0
m n 2 R 1 log (x) log (1−x ) 1 x+1 1 Proof. (1): Differentiate B(a, b). (2)-(4): consider 0 f(p;x) dx. Using 1−x = 1−x2 , x+1 = 1−x 2 2 1−x2 and that x → t give RHS. Alternatively log(1−x ) = log(1−x)+log(1+x) and binomial the- 2 2 log(x+1) x log(1−x ) log(1−x ) log(1−x) orem for expansion give LHS. (5): Multiply both sides of 1−x = 1−x2 + 1−x2 − 1−x
15 n x m 2 log ( 1−x ) x with log (x) then integrate, t → x . (6): Consider F (x) = x+1 , substituting t → 1−x and 1−s R ∞ xs−1 π(2 −1) differentiating 0 (x+1)(2x+1) dx = sin(πs) yield RHS, and binomial theorem yields LHS. (7): 1 n n+1 R 2 log (1−x) log (x) 1−u 0 1−x dx is evaluated in [15]. Let x → 2 , use binomial theorem twice. (8): Sim- n x x log ( x+1 ) ilar to (6), subsitute t → x+1 in integral of F1(x) = x+1 , series expansion. (9): Consider k m n 1 R ∞ log (x) log (x+1) log ( x +1) 1 0 x dx, m, n > 0. Substitution x → t for (1, ∞) part and binomial theo- 1−t rem yields LHS, while x → t for the whole integral and partial fractions give RHS.
Partial derivatives in above expressions can be evaluated by Mathematica (or manually by com- puting Laurent series of polygamma and plugging in their special values). Now consider all possible relations generated from these formulas.
2.6 Contour integration
Integrate logk(1 − z) logm(z) logn(z + 1) F (z) = z(1 − z) (with all logs principal) along the large upper semicircular contour, by Cauchy theorem we have R ∞+0i 1 −∞+0i F (z)dz = 0. Split the integral into 4 parts by branch points 0, ±1, use x → ±t, x → ± t to map them to (0, 1), after expansion we obtain a linear relation (see [1] and sections below for details and examples). Now consider all possible F .
Clarification: Here and below, by ’principal’ we follow the convention of Mathematica, that is, the argument of log/polylogs vary in (−π, π).
2.7 Double integration
In [1] a symmetric double log integral is investigated, which essentially yields:
n−1 2 ! X 1 1 1 1 1 LI(2n, 0, 1; 1) = (2n)! (−1)kLi Li + (−1)nLi + Li k+1 2 −k+2n+1 2 2 n+1 2 2n+2 2 k=0
2.8 Hypergeometric identities
Differentiating hypergeometric identities can be used to generate the most nontrivial high weight relations, see [1].
16 2.9 Solution
All linear relations involving LIs with W ≤ 5 generated by subsection 2-1∼2-6 already yields:
Lemma 1. All LIs with weight W ≤ 5 can be generated by AW over Q.
In fact W = 6 case is also solvable using 2-1∼2-6 only, but we won’t record their closed-forms below due to great complexity (for their values we refer readers to [1]). W = 6 results will be frequently used later so it’s better to be aware of the general structure.
Clarification: In section N, all subjects computed in section 1, ··· ,N − 1 are supposed to be known. As an example, LI values are frequently used below with the statement ’··· modulo LIs’.
3 ES
Since the theory of Euler sums is well developed already, this section won’t be a thorough summary, but a simple review of ES evaluation and solution of several ESs.
P∞ (±1)n−1 Excluding the trivial n=1 np , there lefts 89 ESs with weight W ≤ 5 (2, 9, 24, 54 for weight 2, 3, 4, 5 respectively) .
3.1 General formulas
According to [11], we have
Lemma 2. ES(±1; p),ES(±1; −p)(p even), ES(p; p),ES(p; q)(p + q odd or smaller than 7), ES(1, 1; p)(p odd or smaller than 5) are reducible to Zeta values.
See [2, 7] for alternative proofs.
17 3.2 Symmetric relations
By inclusion-exclusion (abbr. I-E) principle we have:
X X X X X X X X = + + − − − + 0 (±1)j−1(±1)k−1(±1)l−1 Apply this identity as an operator on all triple summands jakblc to obtain linear relations. For instance, if we choose 3 ’±’s to be ’+’ and assuming convergence one obtain ζ(a)ζ(b)ζ(c) = ES(a, b; c)+ES(a, c; b)+ES(b, c; a)−ES(a; b+c)−ES(b; a+c)−ES(c; a+b)+ES(a+b+c) Consider all other choices and n-ple generalizations to obtain relations, see for instance [21, 22]. See also [9] for similar contents on MZVs. 3.3 Partial fractions P∞ P∞ 1 Consider Cauchy product of the double sum j=1 k=1 jakb . By partial fraction decomposition Pn−1 1 we know that all sums of the form k=1 ka(n−k)b are expressible in a combination of homogeneous (k) Hn na+b−k , k = 1, ··· , a + b − 1, we found that ζ(a)ζ(b) is expressed in a corresponding combination of (±1)j−1(±1)k−1 ESs, which is a linear relation. We may also consider jakb and generalizations to triple sums. See [9] for more. 3.4 Contour integration See [11] again for generating ES relations using contour integration and polygamma kernels. 3.5 Solution 3.5.1 Preparation In [11, 19, 21, 22] all but 1 weight 3 ES, 2 weight 4 ESs and 7 weight 5 ESs are solved. 3 weight < 5 ESs are evaluated trivially. Now we will sketch a proof for the remaining 7, that is ES(−1, −1, −1; −2),ES(−1, −1, −1, −1; −1),ES(1, 1, −2; −1),ES(1, 1, −1, −1; −1) ES(1, −1, −1; −2),ES(1, −1, −2; −1),ES(1, −1, −1, −1; −1) 18 All LIs, low weight ESs and other 47 weight 5 ESs are supposed to be known. Also, several GES 1 P∞ Hn (Geometric Euler Sums, i.e. sums involving 2n like n=1 n32n ) computed in [19] will be applied too. Proposition 3. (1)-(15) hold where |x| < 1 and naturally extends to |z| ≤ 1 assuming convergence: ∞ ∞ n 1 X n X (±x) = (±x) , − log(1 ∓ x) = (1) 1 ∓ x n n=0 n=1 ∞ ∞ log(1 ∓ x) X n log(1 ± x) X n − = Hn (±x) , = Hf (±x) (2) 1 ∓ x 1 ∓ x n n=1 n=1 ∞ X 1 Hn n log2(1 ± x) = 2(−1)n−1 − (±x) (3) n2 n n=1 ∞ ! X Hg 1 log(1 − x) log(x + 1) = − 2n + x2n (4) n 2n2 n=1 ∞ n ∞ ∞ X (±x) X n Lik(±x) X (k) n Lik(∓x) = Li (±x), H(k) (±x) = , H] (±x) = − (5) nk k n 1 ∓ x n 1 ∓ x n=1 n=1 n=1 ∞ n X Hnx 1 = Li (x) + log2(1 − x) (6) n 2 2 n=1 ∞ X Hf xn 1 − x 1 n = Li − Li (−x) − Li − log(2) log(1 − x) (7) n 2 2 2 2 2 n=1 ∞ 2 X Li2(x) log (1 − x) (H ) 2xn = + (8) n 1 − x 1 − x n=1 ∞ 1−x 1 x+1 X 2 2 Li2 − Li2 + log log(1 − x) + Li2(x) Hf xn = 2 2 2 (9) n 1 − x n=1 ∞ 1−x 1 1 2 X n Li2 2 − Li2(−x) − Li2 2 − 2 log (x + 1) − log(2) log(1 − x) HnHf x = (10) n 1 − x n=1 ∞ 1 2 1 2 X Li3(1 − x) + Li3(x) + log(x) log (1 − x) − π log(1 − x) − ζ(3) H H(2)xn = 2 6 (11) n n 1 − x n=1 ∞ 3 3 2 1 2 X 3 n 3Li3(1 − x) + Li3(x) − log (1 − x) + 2 log(x) log (1 − x) − 2 π log(1 − x) − 3ζ(3) (Hn) x = (12) 1 − x n=1 19 ∞ 2 n X (Hn) x 1 = Li (x) − Li (x) log(1 − x) − log3(1 − x) (13) n 3 2 3 n=1 ∞ (2) n X Hn x 2 = 2Li3(1 − x) + Li3(x) − 2Li2(1 − x) log(1 − x) − Li2(x) log(1 − x) − log(x) log (1 − x) − 2ζ(3) (14) n n=1 ∞ n X Hnx 1 = −Li (1 − x) + Li (x) + Li (1 − x) log(1 − x) + log(x) log2(1 − x) + ζ(3) (15) n2 3 3 2 2 n=1 R x f(t) Proof. (1): Trivial. (2)-(5): Cauchy product, use (1). (6)-(7): Applying L(f) = 0 t dt on G(x) certain f in (5). (8): Denote LHS as F (x), suppose it is of form F (x) = 1−x , then by Cauchy we P∞ 2 2 n 1 have G(x) = n=1 (Hn) − (Hn−1) x . Using Hn−1 = Hn − n and (5), (6) to arrive at RHS. (9)-(12): Calculating differences like (8). (13)-(15): Make use of L(f) on (5), (6), (8). Mathematica is used to calculate primitives. Proposition 4. (1)-(10) hold where |x| < 1 and naturally extends to |z| ≤ 1 assuming convergence: ∞ 2 n 2 X (Hn) x Li2(x) (1) = − 2Li (1 − x) + Li (x) n2 2 4 4 n=1 1 π4 −Li (1 − x) log2(1 − x) + 2Li (1 − x) log(1 − x) − log(x) log3(1 − x) + 2 3 3 45 ∞ (2) n X H x 1 4 1 2 3 1 2 1 3 1 (2) n = − log − 1 + log(1 − x) log − 1 + log log − 1 n2 6 x 3 x 3 x x n=1 2 2 1 x − 1 2 1 2 1 2 3 1 − log (1 − x) log − 1 + 2Li2 log − 1 − 2Li2(x) log − 1 + log (1 − x) log − 1 x x x x 3 x 2 1 1 1 +2 log (1 − x) log(x) log − 1 + 4 log(1 − x)Li2(1 − x) log − 1 + 4 log(1 − x)Li2(x) log − 1 x x x 4 1 x − 1 1 1 π 1 4 x −4Li3(1 − x) log − 1 − 4Li3 log − 1 − 4Li3(x) log − 1 + − log x x x x 18 12 1 − x 1 1 1 x 1 x 2 1 x − log4(1 − x) + log log3 + log log3(x) − log(1 − x) log3(x) − log2(x) log2 6 3 1 − x 1 − x 3 1 − x 3 2 1 − x 2 3 x Li2(x) 1 x 4 + log2(1 − x) log2(x) + log log(1 − x) log2(x) + + log3 log(x) − log3(1 − x) log(x) 2 1 − x 2 3 1 − x 3 2 x x 2 x 2 2 + log Li2 − − log Li2(1 − x) − 2 log (1 − x)Li2(1 − x) − log (x)Li2(1 − x) 1 − x 1 − x 1 − x x 2 2 +2 log log(x)Li2(1 − x) + 2 log(1 − x) log(x)Li2(1 − x) − 2 log (1 − x)Li2(x) − log (x)Li2(x) 1 − x 20 x x x x +2 log log(x)Li2(x) + 2 log(1 − x) log(x)Li2(x) − 2 log Li3 − − 2 log Li3(1 − x) 1 − x 1 − x 1 − x 1 − x x x x − 1 log4(x) −2 log Li3(x) + 2 log(1 − x)Li3(x) + 2Li4 − + 2Li4(1 − x) + 4Li4 − Li4(x) − 2 log(x)ζ(3) − 1 − x 1 − x x 12 ∞ n X Hnx 1 4 1 1 3 1 1 1 3 1 (3) = log − 1 − log(1 − x) log − 1 − log log − 1 n3 12 x 3 x 3 x x n=1 1 2 2 1 x − 1 2 1 2 1 1 3 1 + log (1 − x) log − 1 − Li2 log − 1 + Li2(x) log − 1 − log (1 − x) log − 1 2 x x x x 3 x 2 1 1 1 − log (1 − x) log(x) log − 1 − 2 log(1 − x)Li2(1 − x) log − 1 − 2 log(1 − x)Li2(x) log − 1 x x x 1 x − 1 1 1 1 4 x 1 4 +2Li3(1 − x) log − 1 + 2Li3 log − 1 + 2Li3(x) log − 1 + log + log (1 − x) x x x x 24 1 − x 12 log4(x) 1 1 x 1 x 1 1 x + − log log3 − log(x) log3 + log(1 − x) log3(x) + log2 log2(x) 24 6 1 − x 1 − x 6 1 − x 3 4 1 − x 3 1 x 1 x 2 − log2(1 − x) log2(x) − log3(x) log − log(1 − x) log2(x) log + log3(1 − x) log(x) 4 6 1 − x 2 1 − x 3 1 2 x 2 1 2 x + log Li2(1 − x) + log (1 − x)Li2(1 − x) + log (x)Li2(1 − x) − log(x) log Li2(1 − x) 2 1 − x 2 1 − x 1 2 x x 2 1 2 − log(1 − x) log(x)Li2(1 − x) − log Li2 + log (1 − x)Li2(x) + log (x)Li2(x) 2 1 − x x − 1 2 x x x x − log(x) log Li2(x) − log(1 − x) log(x)Li2(x) + log Li3(1 − x) + log Li3 1 − x 1 − x 1 − x x − 1 x x x − 1 π4 + log Li3(x) − log(1 − x)Li3(x) − Li4(1 − x) − Li4 − 2Li4 + 2Li4(x) + log(x)ζ(3) − 1 − x x − 1 x 36 ∞ (2) n 2 2 2 X HnHn x Li2(1 − x) Li2(x) π Li2(x) x (4) = − + + Li4(1 − x) − Li4 n 2 2 6 x − 1 n=1 1 2 x 1 x 2 x 1 2 1 2 + Li2(1 − x) log − Li2 log + Li2(1 − x) log (1 − x) + Li2(1 − x) log (x) 2 1 − x 2 x − 1 1 − x 2 2 1 2 x x x + Li2(x) log (x) − Li2(1 − x) log(x) log − Li2(x) log(x) log + Li3(1 − x) log 2 1 − x 1 − x 1 − x x x x +Li3 log + Li3(x) log − Li2(x) log(1 − x) log(x) − Li3(1 − x) log(1 − x) + ζ(3) log(1 − x) x − 1 1 − x 1 − x 1 x log4(x) 1 1 x 1 x 1 x −ζ(3) log(x) + log4 + − log log3 − log(x) log3 − log3(x) log 24 1 − x 24 6 1 − x 1 − x 6 1 − x 6 1 − x 1 1 x 1 x 1 1 π4 + log(1−x) log3(x)+ log2(x) log2 − log(1−x) log2(x) log + π2 log2(1−x)− log2(1−x) log2(x)− 3 4 1 − x 2 1 − x 12 4 40 ∞ (3) n 2 X Hn x Li2(x) (5) = − + Li (x) − Li (x) log(1 − x) n 2 4 3 n=1 21 ∞ 3 n 2 2 2 ! ! X (Hn) x 3Li2(1 − x) Li2(x) π Li2(x) x 3 2 x (6) = − + − 3Li4(1 − x) − 3Li4 − 2Li4(x) + Li2(1 − x) log n=1 n 2 2 2 x − 1 2 1 − x ! ! ! 3 x 2 x 3 2 3 2 3 2 x − Li2 log − Li2(1 − x) log (1 − x) + Li2(1 − x) log (x) + Li2(x) log (x) − 3Li2(1 − x) log(x) log 2 x − 1 1 − x 2 2 2 1 − x x ! x ! x ! x ! x ! −3Li2(x) log(x) log + 3Li3(1 − x) log + 3Li3 log + 3Li3(x) log − 3Li2(x) log(1 − x) log(x) 1 − x 1 − x x − 1 1 − x 1 − x ! 4 1 4 x 1 4 log (x) +3Li3(1 − x) log(1 − x) + 2Li3(x) log(1 − x) + 3ζ(3) log(1 − x) − 3ζ(3) log(x) + log + log (1 − x) + 8 1 − x 4 8 ! ! ! ! 1 1 3 x 1 3 x 1 3 x 3 3 − log log − log(x) log − log (x) log + log(1 − x) log (x) − log (1 − x) log(x) 2 1 − x 1 − x 2 1 − x 2 1 − x ! ! 4 3 2 2 x 3 2 x 1 2 2 3 2 2 π + log (x) log − log(1 − x) log (x) log + π log (1 − x) − log (1 − x) log (x) − 4 1 − x 2 1 − x 4 4 120 ∞ n ! ! ! ! !! X Hgnx x + 1 x 2x x + 1 2x (7) = −Li − Li (−x) − Li (x) − Li + Li + log(2)Li (x) + Li log(x) − log 2 3 3 3 3 3 2 2 n=1 n 2 x + 1 x + 1 2 x + 1 1 x x ! ! x ! 2x !! 2x ! +Li2 − log(x) + Li2(x) log + log(x + 1) + Li2 − Li2 log 2 2 x + 1 x + 1 x + 1 x + 1 ! ! ! !! ! 1 2 1 2 2x 1 − x 1 x − 1 2x + log (2) log(x) + log log + log − log − + log(2) log(x) log 2 2 x + 1 2 x + 1 x + 1 x + 1 ! ! 3 1 1 2 1 − x x + 1 7ζ(3) log (2) − log(2) log(x)(log(x) − 2 log(x + 1) + log(4)) − π log(x) + log log log(x) + + 2 12 2 2 8 6 ∞ n ! ! ! ! ! X HnHgnx 1 − x x + 1 x + 1 x + 1 x (8) = −Li3 − 2Li3 + Li3 − Li3(1 − x) − Li3 − Li3(−x) − Li3(x) − Li3 n=1 n 2 2 1 − x x − 1 x + 1 ! ! ! 2x 1 2 1 − x 1 − x +Li3 + log(2) 2Li2(x) + log (1 − x) + Li2 log(1 − x) − Li2 (log(1 − x) − log(x)) x + 1 2 2 2 x + 1 ! 2x !! x + 1 ! x + 1 !! +Li2 log(x) − log + Li2 log(x + 1) + Li2(1 − x) log(x + 1) − log 2 x + 1 2 1 − x x + 1 ! x + 1 !! x + 1 ! x ! ! + Li2 − Li2 log + Li2(−x)(log(1 − x) − log(x)) + Li2(−x) log(x) + Li2(x) log + log(x + 1) x − 1 1 − x 1 − x x + 1 x ! 2x !! 2x ! x + 1 ! ! + Li2 − Li2 log + Li2(x + 1) log + log(1 − x) − Li2(x + 1) log(x + 1) x + 1 x + 1 x + 1 1 − x ! ! ! ! ! 1 1 − x 2 1 2 1 2 2 1 2 x + 1 2x 2 + log log (x+1)− log(−x) log (x+1)+ π − 6 log (2) log(1−x)+ log − log + log − + log(x) 2 2 2 12 2 1 − x x − 1 x − 1 ! ! ! !! ! 1 2 2 1 2 2x 1 − x 1 x − 1 2x − π − 6 log (2) log(x) + log log + log − log − + log(2) log(x) log 12 2 x + 1 2 x + 1 x + 1 x + 1 1 1 − x ! x + 1 ! x + 1 ! − log(2) log(x)(log(x) − 2 log(x + 1) + log(4)) + log log log(x) + log (log(−x) − log(x)) log(x + 1) 2 2 2 1 − x 3 1 15ζ(3) log (2) 1 2 + log(1 − x) log(x) log(x + 1) + (log(x) − log(−x)) log(x + 1)(log(x + 1) − 2 log(1 − x)) + + − π log(2) 2 8 2 6 2 ∞ H xn ! ! ! ! ! ! X gn 1 − x x + 1 x + 1 x + 1 x 2x (9) = −4Li3 − 2Li3 + 2Li3 − 2Li3 − Li3(x) − 2Li3 + 2Li3 n=1 n 2 2 1 − x x − 1 x + 1 x + 1 ! ! ! !! 2 1 − x 1 − x x + 1 2x −2Li3(x+1)+log(2) 2Li2(x) + log (1 − x) +4Li2 log(1−x)−2Li2 (log(1−x)−log(x))+2Li2 log(x) − log 2 2 2 x + 1 x + 1 !! x + 1 ! x + 1 !! x + 1 ! +2Li2(1 − x) log(x + 1) − log − 2Li2(1 − x) log(1 − x) + 2 Li2 − Li2 log − Li2(x)(log(1 − x) − log(x)) 1 − x x − 1 1 − x 1 − x 22 x ! ! x ! 2x !! 2x ! x + 1 ! ! −Li2(x) log(x)+2Li2(x) log + log(x + 1) +2 Li2 − Li2 log +2Li2(x+1) log + log(1 − x) x + 1 x + 1 x + 1 x + 1 1 − x ! ! ! ! 2 1 2 2 2 x + 1 2x 2 2 − log(2) log (1 − x) + π − 6 log (2) log(1 − x) + log − log + log − + log(x) − log (1 − x) log(x) 6 1 − x x − 1 x − 1 ! ! ! !! ! 1 2 2 2 2x 1 − x 1 x − 1 2x − π − 6 log (2) log(x) + log log + log − log − + 2 log(2) log(x) log 6 x + 1 2 x + 1 x + 1 x + 1 1 − x ! x + 1 ! x + 1 ! − log(2) log(x)(log(x) − 2 log(x + 1) + log(4)) + 2 log log log(x) + 2 log (log(−x) − log(x)) log(x + 1) 2 2 1 − x 15ζ(3) 3 1 2 +2 log(1 − x) log(x) log(x + 1) + (log(x) − log(−x)) log(x + 1)(log(x + 1) − 2 log(1 − x)) + + log (2) − π log(2) 4 3 (10) The following four generating functions are polylog expressible, using formula (1)-(6): ∞ ∞ ∞ ∞ X 4 n X 2 (2) n X (2) 2 n X (3) n (Hn) x , (Hn) Hn x , Hn x , HnHn x n=1 n=1 n=1 n=1 Proof. (1)-(9): Again, make use of L(f) in appropriate functions in Proposition 3, then Math- ematica. They are verified via differentiation. (10): Calculate differences. The results will be combinations of (1)-(6) thus extremely lengthy, so we omit them. Proposition 5. The following formulas hold: Z 1 (−1)kk! n−1 k (1) x log (x) dx = k+1 0 n Z 1 ∞ X 1 Hn xn−1 log(1 − x) dx = − = − (2) j(j + n) n 0 j=1 Z 1 ∞ 2 (2) X 1 (Hn) + Hn xn−1 log2(1 − x) dx = = (3) j(j + n)2 n 0 j=1 Z 1 H H(2) π2 n−1 n n (4) x log(x) log(1 − x) dx = 2 + − 0 n n 6n Z 1 3H H(2) + (H ) 3 + 2H(3) xn−1 log3(1 − x) dx = − n n n n (5) 0 n ∞ Z 1 X (−1)j (−1)n−1Hf − (−1)n−1 log(2) − log(2) − xn−1 log(x + 1) dx = = n (6) j(j + n) n 0 j=1 ∞ j n−1 (2) n−1 n−1 2 n−1 X (−1) (−1) Hgn (−1) Hf (−1) log(2) log(2) π (−1) = + n − − − (7) j(j + n)2 n n2 n2 n2 12n j=1 23 Proof. (1): Trivial. (2)-(5): Differentiate Beta functions, then transform polygamma constants into (2) (3) (3) π4 (4) harmonic ones (for instance ψ (n + 1) = 2Hn − 2ζ(3), ψ (n + 1) = 15 − 6Hn , etc). (6)-(7): π2 Invoke partial fractions and that η(1) = log(2), η(2) = 12 . 3.5.2 Confrontation Now we compute 7 remaining weight 5 Euler sums: Proposition 6. The following formulas hold: (1) ES(−1, −1, −1; −2) 1 π2ζ(3) 1229ζ(5) 21 3 log5(2) 3 29 = −18Li − + + ζ(3) log2(2) + + π2 log3(2) − π4 log(2) 5 2 8 64 8 20 4 160 (2) ES(−1, −1, −1, −1; −1) 1 1 5π2ζ(3) 733ζ(5) 13 log5(2) 3 71 = −48Li − 4Li log(2) + + + + π2 log3(2) − π4 log(2) 5 2 4 2 48 16 30 2 180 (3) ES(1, 1, −2; −1) 1 π2ζ(3) 93ζ(5) 1 1 13π4 log(2) = 6Li − − − log5(2) + π2 log3(2) − 5 2 48 64 20 36 1440 (4) ES(1, 1, −1, −1; −1) 1 π2ζ(3) 201ζ(5) 3 log5(2) 1 97 = 12Li + − − ζ(3) log2(2) + + π2 log3(2) + π4 log(2) 5 2 16 16 8 10 12 720 (5) ES(1, −1, −1; −2) 1 1 π2ζ(3) 7 1 2 1 = 10Li −2Li log(2)+ −10ζ(5)− ζ(3) log2(2)− log5(2)+ π2 log3(2)+ π4 log(2) 5 2 4 2 96 8 6 9 10 (6) ES(1, −1, −2; −1) 1 7 log5(2) 5 1 = 4Li log(2) + ζ(3) log2(2) + − π2 log3(2) − π4 log(2 4 2 16 6 72 480 (7) ES(1, −1, −1, −1; −1) 24 1 13π2ζ(3) 83ζ(5) 9 13 41 1 = 2Li + − + ζ(3) log2(2) − log5(2) + π2 log3(2) − π4 log(2) 5 2 96 32 16 60 72 180 Proof. (1)-(2): Exploit symmetry. Take (2) as an example, we have X (−1)j−1(−1)k−1(−1)l−1(−1)m−1(−1)n−1 ES(−1, −1, −1, −1; −1) = jklmn 0 Consider the analogous sum without restrictions, by I-E principle log5(2) = 5ES(−1, −1, −1, −1; −1)−10ES(−1, −1, −1; 2)+10ES(−1, −1; −3)−5ES(−1; 4)+η(5) Since all ESs but the desired one are known, (2) is established. (1) is similar. 2 (3): Multiply both sides of Prop. 5(7) with (Hn) then sum w.r.t n, we have: π2 − ES(1, 1; −1) − log(2)ES(1, 1; 2) − log(2)ES(1, 1; −2) 12 ∞ ∞ j 2 X X (−1) (Hn) +ES(1, 1, −1; −2) + ES(1, 1, −2; −1) = j(j + n)2 n=1 j=1 Using formulas above we have: ∞ ∞ j 2 Z 1 ∞ j−1 j ∞ ! X X (−1) (Hn) log(x) X (−1) x X = (H ) 2xn dx j(j + n)2 x j n n=1 j=1 0 j=1 n=1 Z 1 1 1 2 = + log(x) log(x + 1) Li2(x) + log (1 − x) dx 0 x 1 − x Therefore we only need to evaluate P LI(0, 1, 1; 21; 0), P LI(0, 1, 1; 21; 1). Indeed: ∞ ∞ j−1 R 1 j+n−1 ∞ ∞ X X (−1) log(x)x dx X X (−1)j P LI(0, 1, 1; 21; 1) = 0 = jn2 jn2(j + n)2 n=1 j=1 n=1 j=1 ∞ ∞ X X Hfj P LI(0, 1, 1; 21; 0) = ··· = − n2(j + n + 1)2 n=1 j=1 For the first one, sum w.r.t either n or j first and the PLI reduces to known ESs. For the second, sum w.r.t n yields ∞ ∞ ∞ j−1 (2) 2 ! X X Hfj X (−1) 2Hj Hj π − = − Hfj − − − + n2(j + n + 1)2 j j3 j2 3j2 n=1 j=1 j=1 25 All ES components on RHS are also known so (3) is proved. n−1 2 (4): Multiply both sides of Prop. 5(3) with (−1) Hfn and sum n, the ES boils down to P LI(2, 0, 0; 22; 1/2), P LI(2, 0, 0; 24; 1/2) via Prop. 3(9). For P LI(2, 0, 0; 22; 1/2), use Prop. 3(5) for expansion to reduce them to known ESs. For P LI(2, 0, 0; 24; 1/2), use that: 1 − x x + 1 π2 Li +Li = log(2) log(1−x)+log(2) log(x+1)−log(1−x) log(x+1)+ −log2(2) 2 2 2 2 6 to reduce them to P LI(2, 0, 0; 23; 1/2) modulo LIs. For the final 2 PLIs, by substitution t → 1 − x, x x Li2( 2 ) Li2( 2 ) single/double expansion for x / respectively and ES/GES results [19] they are solved 1− 2 1−x either so (4) is proved. (5): Same as (4), except that we start from Prop. 5(4). One may verify that all ESs appeared during reduction are known and all PLIs can be reduced to known ESs. (6): Multiply both sides of Prop. 5(7) with HnHfn, remaining steps are alike to (3)-(5), in which the value of (5) is used. P∞ 3 n (7): Denote f(x) = n=1 Hfn x , by calculating differences of the generating function we H 2xn n n −3P (x)−3Q(x)−Li3(x) P∞ (gn) P∞ (−1) Hgnx have f(−x) = x+1 , where P (x) = n=1 n ,Q(x) = n=1 n2 . Thus Z 1 f(−x) log(1 − x) Z 1 log(1 − x)(−3P (x) − Li (x) − 3Q(x)) ES(1, −1, −1, −1; −1) = dx = 3 dx 0 x 0 x(x + 1) Integrate by parts using Z log(1 − x) π2 1 − x x + 1 dx = u(x) = − Li − Li (x) − log(1 − x) log x(x + 1) 6 2 2 2 2 We only have to compute Z 1 ∞ n ∞ ! u(x) X (−1) n X 2 n 3 Hfnx + 3 Hfn x + Li2(x) dx x n 0 n=1 n=1 Which, by Proposition 3, equals to Z 1 u(x) 1 − x 1 x + 1 3 2Li2 + Li2(x) − 2Li2 + 2 log log(1 − x) dx 0 x(1 − x) 2 2 2 26 Z 1 Z 1 u(x) x + 1 1 Li2(x)u(x) +3 Li2 − Li2(x) − Li2 − log(2) log(x + 1) dx + dx 0 x 2 2 0 x Expand all these expressions. Some parts of the integrand are divergent so we put them together x+1 1 R 1 Li2(x)(Li2( 2 )−Li2( 2 )) again, for instance 0 x dx. Apply integration by parts on all these composite integrals (in the example lift up a Li3(x)) to convert them into convergent PLIs. After this we reduce it to plenty of LIs and PLIs. For each of them, integrate by parts and use techniques in (3)-(6) over and over again (or more generally, use the expansion algorithm in section 4), all PLIs can be reduced to lower weight ESs and 53 weight 5 ESs already known, hence (7) is proved. Thus all 89 ESs with W ≤ 5 are evaluated. 4 PLI Due to existence of a polylog term, a PLI’s weight W ≥ 3. There are 11, 55, 197 PLIs for weight 3, 4, 5 respectively. 4.1 Brute force x+1 Li2( 2 ) log(1−x) Some of the PLIs can be calculated by brute force, for instance x+1 has primitive x+1 1 x+1 2 log(2)Li3 2 − 2 Li2 2 . For each PLI, check by Mathematica whether the primitive is solv- able and evaluate it if so. 4.2 General formulas See [12, 15] for proofs of the following: m n R 1 Lim(x)Lin(x) R 1 log (x)Lin(x) R 1 Lin+1(x) log (x) Lemma 3. 0 x dx, 0 1−x dx (m + n even) and 0 x+1 dx are re- ducible to Zeta values. Thus PLIs of these types are solved directly. 27 4.3 Integration by parts Q2 a(m) Q6 QN b(n,k) In spirit of subsection 2-3, generalize F to m=0 l(m; x) k=1 n=2 L(n, k; x) . For those involving L(n, k; x), k ∈ {3, 5, 6}, after partial fractions another integration by parts lifting log(1−x) is needed to obtain a regular relation modulo LIs. Now consider all possible F . 4.4 Fractional transformation Q2 a(m) Q6 QN b(n,k) 1−x m=0 l(m;x) k=1 n=2 L(n,k;x) Similar to subsection 2-4, apply x → x+1 in integral of F = f(p;x) , where b(n, 3) = b(n, 4) = 0. For p = 1, after partial fractions another integration by parts lifting log(1 − x) is needed to obtain a regular relation modulo LIs. Now consider all possible F . 4.5 Complementary integrals QN a(n) 1±x Denote P (x) = n=2 Lin(x) , by equally splitting (0, 1) and substituting t = 2 we have: 1−x m 1−x n x+1 x+1 m x+1 n 1−x m n Z 1 P 2 log 2 log 2 P 2 log 2 log 2 Z 1 P (x) log (x) log (1 − x) + dx = dx 0 1 ∓ x 1 ± x 0 f±(x) Here f+(x) = x, f−(x) = 1 − x. Expand LHS to get a linear relation. Now consider all possible P . 4.6 Polylog identities Proposition 7. The following formulas hold on (0, 1) and naturally extend to the plane, where all logs/polylogs are principal: Li x2 Li (−x) + Li (x) = n (1) n n 2n−1 π2 Li (1 − x) + Li (x) = − log(x) log(1 − x) (2) 2 2 6 x 1 Li + Li (x) = − log2(1 − x) (3) 2 x − 1 2 2 1 π2 1 Li − + Li (−x) = − − log2(x) (4) 2 x 2 6 2 2 1 − x x + 1 π 2 Li2 + Li2 = log(2) log(x + 1) + log(2) log(1 − x) − log(1 − x) log(x + 1) + − log (2) (5) 2 2 6 28 1 − x x − 1 1 1 Li + Li = − log2(x + 1) + log(2) log(x + 1) − log2(2) (6) 2 2 2 x + 1 2 2 1 1 1 Li (−x) − Li − = − log3(x) − π2 log(x) (7) 3 3 x 6 6 x 1 1 1 Li + Li (1 − x) + Li (x) = log3(1 − x) − log(x) log2(1 − x) + π2 log(1 − x) + ζ(3) (8) 3 x − 1 3 3 6 2 6 1 − x x + 1 x − 1 1 1 Li + Li + Li = log3(x + 1) − log2(2) log(1 − x) 3 2 3 2 3 x + 1 6 2 1 π2 log2(2) − log2(x + 1) log(1 − x) + − log(x + 1) (9) 2 6 2 log3(2) 1 + log(2) log(1 − x) log(x + 1) + ζ(3) + − π2 log(2) 3 6 1−x Proof. (1)-(4), (7)-(8): Classical formulas, see [14]. (5), (6), (9): Substitute t → 2 in (2), (3), (8) respectively. logm(1−x2) logk(x) Now, multiply both sides of (1) with F (x) = x , integrate on (0, 1). Substitute t → x2 on RHS, expand (log(1 − x) + log(x + 1))m on LHS, we obtain a relation. Another choice x logm(1−x2) logk(x) is F (x) = 1−x2 , where we also use partial fractions to decompose LHS. Now consider all possible F . Q2 a(m) m=0 l(m;x) Moreover, multiply both sides of (9) with F (x) = f(p;x) , it generates a relation between 3 PLIs since all LIs on RHS are known. For (5), (6) the process is similar, except that we may Q2 a(m) m=0 l(m;x) L(2,k;x) broaden the range of F to f(p;x) and f(p;x) . Now consider all possible F for 3 formulas to engender linear relations. Clarification: If we derive all closed-forms of PLIs involving L(2, 2; x) hence L(2, 6; x) (via frac- tional transformation), those involving L(2, 3/4; x) will be natural corollaries thus redundant in the systematic solution. However, as we will see in the next subsection, solving all L(2, 2; x) type PLIs directly is unachievable. As a matter of fact, some of the independent relations only generate from 1−x manipulating L(2, 3/4; x) type, say subsection 4-7 (reflection of Lin 2 ), 6-10-4, hence this part is necessary. 29 4.7 Series expansion Weight 5 PLIs are sorted into 4 classes, that is, (2,1,1), (3,1), (4), (2,2). The index indicates the quantity as well as order of polylog terms, for instance, P LI(1, 1, 0; 22; 2) belongs to (2,1,1) class for it contains one Li2 term and two log terms in numerator. Note that class (1,1,1,1,1) is excluded for they are simply LIs. From now on we only discuss about weight 5 PLIs for lower weights are relevantly trivial. More- over, we restrict the use of this method to (2,1,1), (3,1) and (4) class PLIs. Furthermore, we require the only polylog term L(n, k; x)(n = 2, 3, 4) contained in the PLI to satisfy 1 ≤ k ≤ 3. Denote the set of these PLIs as Q. We only present the general procedure of converting the PLI into combination of ESs. In this subsection ’ES’ represents all ordinary ESs and all GESs.