An alternating Euler Sum

In this post, we will evaluate the famous Euler sum: n=1(1)n+1Hnn3()\sum_{n=1}^\infty (-1)^{n+1}\frac{H_n}{n^3}\quad \color{blue}{\cdots (*)} where Hn=k=1n1k=011tn1tdtH_n = \sum_{k=1}^n \frac{1}{k}=\int_0^1 \frac{1-t^n}{1-t}dt is the nn-th harmonic number. This series resists contour integration techniques which makes it’s computation quite challenging. We will start with the following well known generating function identity: n=1Hntn=log(1t)1t,1t<1\sum_{n=1}^\infty H_n t^n = -\frac{\log(1-t)}{1-t} ,\quad -1 \leq t < 1 Dividing the above equation by tt and integrating both sides, gives us: n=1Hnxnn=0xlog(1t)t(1t)dt=0x(1t+11t)log(1t)  dt=Li2(x)+12log2(1x) \begin{aligned} \sum_{n=1}^\infty H_n \frac{x^n}{n} &= -\int_0^x\frac{\log(1-t)}{t(1-t)}dt \\ &= -\int_0^x \left(\frac{1}{t} + \frac{1}{1-t} \right)\log(1-t) \; dt \\ &= \text{Li}_2(x) + \frac{1}{2}\log^2(1-x) \end{aligned} where 1x<1-1\leq x < 1.

1. Computation of n=1Hnn2xn\sum_{n=1}^\infty \frac{H_n}{n^2}x^n

Now, we will use the same process to evaluate n=1Hnn2xn\sum_{n=1}^\infty \frac{H_n}{n^2}x^n but we’ll have to divide the evaluation into two parts.

(I) First consider the case when 0x<10 \leq x < 1. Actually, we don’t need this case for evaluating (*) but I’ll do it anyway for completeness. We have: n=1Hnn2xn=0x(Li2(t)t+12log2(1t)t)dt=Li3(x)+120xlog2(1t)tdt=Li3(x)+12log(1x)0t2et1etdt(t1et)=Li3(x)+12n=1log(1x)0t2ent  dt=Li3(x)+12n=1(2n3log2(1x)n(1x)n+2log(1x)n2(1x)n2(1x)nn3)=Li3(x)+log2(1x)log(x)2+log(1x)Li2(1x)Li3(1x)+ζ(3)(1) \begin{aligned} \sum_{n=1}^\infty \frac{H_n}{n^2} x^n &= \int_0^x\left(\frac{\text{Li}_2(t)}{t} + \frac{1}{2}\frac{\log^2(1-t)}{t}\right)dt \\ &= \text{Li}_3(x) + \frac{1}{2}\int_0^x \frac{\log^2(1-t)}{t} dt \\ &= \text{Li}_3(x) + \frac{1}{2}\int_{\log(1-x)}^0 \frac{t^2 e^t}{1-e^t}dt \quad (t\mapsto 1-e^t) \\ &= \text{Li}_3(x) + \frac{1}{2}\sum_{n=1}^\infty \int_{\log(1-x)}^0 t^2 e^{nt} \; dt \\ &= \text{Li}_3(x) + \frac{1}{2}\sum_{n=1}^\infty \left(\frac{2}{n^3} -\frac{\log^2(1-x)}{n}(1-x)^n+\frac{2\log(1-x)}{n^2}(1-x)^n-2\frac{(1-x)^n}{n^3}\right) \\ &= \text{Li}_3(x) + \frac{\log^2(1-x)\log(x)}{2} + \log(1-x)\text{Li}_2(1-x) - \text{Li}_3(1-x) + \zeta(3) \quad\color{blue}{\cdots (1)} \end{aligned} (II) Similarly, for the case 1x<0-1\leq x < 0, we obtain:

n=1Hnn2xn=Li3(x)12x0log2(1t)tdt=Li3(x)+120xlog2(1+t)tdt(tt)=Li3(x)+120log(1x)t2etet1dt(tet1)=Li3(x)+12n=00log(1x)t2entdt=Li3(x)+log3(1x)6+12n=10log(1x)t2entdt=Li3(x)+log3(1x)6+12n=1(2n3log2(1x)(1x)nn2log(1x)(1x)nn22(1x)nn3)=Li3(x)log3(1x)3+log2(1x)log(x)2log(1x)Li2(11x)Li3(11x)+ζ(3)(2) \begin{aligned} \sum_{n=1}^\infty \frac{H_n}{n^2}x^n &= \text{Li}_3(x) -\frac{1}{2}\int_x^0 \frac{\log^2(1-t)}{t}dt \\ &= \text{Li}_3(x) +\frac{1}{2}\int_{0}^{-x} \frac{\log^2(1+t)}{t}dt \quad (t\mapsto -t) \\ &= \text{Li}_3(x) +\frac{1}{2}\int_{0}^{\log(1-x)} \frac{t^2 e^t}{e^t-1}dt \quad (t\mapsto e^t-1) \\ &= \text{Li}_3(x) +\frac{1}{2} \sum_{n=0}^\infty \int_0^{\log(1-x)} t^2 e^{-nt} dt \\ &= \text{Li}_3(x) +\frac{\log^3(1-x)}{6} + \frac{1}{2}\sum_{n=1}^\infty \int_0^{\log(1-x)} t^2 e^{-nt} dt \\ &= \text{Li}_3(x) +\frac{\log^3(1-x)}{6} + \frac{1}{2}\sum_{n=1}^\infty \left(\frac{2}{n^3} -\log^2(1-x)\frac{(1-x)^{-n}}{n}-2\log(1-x)\frac{(1-x)^{-n}}{n^2}-2\frac{(1-x)^{-n}}{n^3}\right) \\ &= \text{Li}_3(x) -\frac{\log^3(1-x)}{3} +\frac{\log^2(1-x)\log(-x)}{2} -\log(1-x)\text{Li}_2\left( \frac{1}{1-x}\right)-\text{Li}_3\left(\frac{1}{1-x}\right)+\zeta(3) \\ &\quad\color{blue}{\cdots (2)} \end{aligned}

2. Evaluation of n=1Hnn3xn\sum_{n=1}^\infty \frac{H_n}{n^3}x^n

This is the hardest step in the evaluation. Once again, we will divide the evaluation into two parts.

(I) Like section 1, we’ll first consider the case when 0x<10\leq x < 1.

Using integration by parts, we obtain: 0xζ(3)Li3(1t)tdt=log(t)(ζ(3)Li3(1t))0x0xlog(t)Li2(1t)1tdt=log(x)(ζ(3)Li3(1x))12Li22(1x)+ζ2(2)2(3) \begin{aligned} \int_0^x \frac{\zeta(3)-\text{Li}_3(1-t)}{t}dt &= \log(t)\left(\zeta(3)-\text{Li}_3(1-t) \right) \Big|_{0}^x - \int_0^x \frac{\log(t)\text{Li}_2(1-t)}{1-t}dt \\ &= \log(x)\left(\zeta(3)-\text{Li}_3(1-x) \right) - \frac{1}{2}\text{Li}_2^2 (1-x) + \frac{\zeta^2 (2)}{2} \quad \color{blue}{\cdots (3)} \end{aligned}

0xlog(1t)Li2(1t)tdt=log(t)log(1t)Li2(1t)0x0xlog(t)(Li2(1t)1t+log(1t)log(t)1t)dt=log(x)log(1x)Li2(1x)+12Li22(1x)ζ2(2)20xlog2(t)log(1t)1tdt(4) \begin{aligned} \int_0^x \frac{\log(1-t)\text{Li}_2(1-t)}{t}dt &= \log(t)\log(1-t)\text{Li}_2(1-t) \Big|_0^x - \int_0^x \log(t) \left(\frac{-\text{Li}_2(1-t)}{1-t} + \frac{\log(1-t)\log(t)}{1-t}\right)dt \\ &= \log(x)\log(1-x)\text{Li}_2(1-x) + \frac{1}{2}\text{Li}_2^2(1-x)-\frac{\zeta^2(2)}{2} - \int_0^x \frac{\log^2(t) \log(1-t)}{1-t}dt \\ &\quad \color{blue}{\cdots (4)} \end{aligned}

0xlog(t)log2(1t)2tdt=14log2(t)log2(1t)0x+120xlog2(t)log(1t)1tdt=14log2(x)log2(1x)+120xlog2(t)log(1t)1tdt(5) \begin{aligned} \int_0^x \frac{\log(t)\log^2(1-t)}{2t}dt &= \frac{1}{4}\log^2(t)\log^2(1-t)\Big|_0^x +\frac{1}{2}\int_0^x \frac{\log^2(t)\log(1-t)}{1-t}dt \\ &= \frac{1}{4}\log^2(x)\log^2(1-x) +\frac{1}{2}\int_0^x \frac{\log^2(t)\log(1-t)}{1-t}dt \quad \color{blue}{\cdots (5)} \end{aligned}

Putting together the results of equations (3), (4) and (5) gives us: n=1Hnn3xn=Li4(x)+14log2(x)log2(1x)+log(x)log(1x)Li2(1x)+log(x)(ζ(3)Li3(1x))120xlog2(t)log(1t)1tdt(6) \begin{aligned} \sum_{n=1}^\infty \frac{H_n}{n^3}x^n &= \text{Li}_4(x) +\frac{1}{4}\log^2(x)\log^2(1-x) + \log(x)\log(1-x)\text{Li}_2(1-x)\\ &\quad + \log(x)\left(\zeta(3)-\text{Li}_3(1-x) \right) -\frac{1}{2}\int_0^x \frac{\log^2(t)\log(1-t)}{1-t}dt \quad \color{blue}{\cdots (6)} \end{aligned}

(II) Now, let 1x<0-1\leq x < 0. We have: x0ζ(3)Li3(11t)tdt=0x1xζ(3)Li3(1t)t(1t)dt(tt1t)=0x1xζ(3)Li3(1t)tdt+0x1xζ(3)Li3(1t)1tdt=log(x1x)(ζ(3)Li3(11x))12Li22(11x)+ζ2(2)2+ζ(3)log(1x)+Li4(11x)ζ(4)(7) \begin{aligned} -\int_x^0 \frac{\zeta(3)-\text{Li}_3\left(\frac{1}{1-t}\right)}{t}dt &= \int_0^{\frac{-x}{1-x}}\frac{\zeta(3)-\text{Li}_3(1-t)}{t(1-t)}dt \quad \left(t\mapsto\frac{-t}{1-t} \right)\\ &= \int_0^{\frac{-x}{1-x}}\frac{\zeta(3)-\text{Li}_3(1-t)}{t}dt + \int_0^{\frac{-x}{1-x}}\frac{\zeta(3)-\text{Li}_3(1-t)}{1-t}dt \\ &= \log\left(\frac{-x}{1-x}\right)\left(\zeta(3) - \text{Li}_3\left(\frac{1}{1-x}\right) \right) - \frac{1}{2}\text{Li}_2^2\left(\frac{1}{1-x}\right) + \frac{\zeta^2(2)}{2} \\ &\quad + \zeta(3)\log(1-x) + \text{Li}_4\left(\frac{1}{1-x} \right) - \zeta(4) \quad \color{blue}{\cdots (7)} \end{aligned} Here, we used equation (3) to simplify the integral after the transformation. Similarly, we obtain: x0log(1t)Li2(11t)tdt=0x1xlog(1t)Li2(1t)t(1t)dt=0x1xlog(1t)Li2(1t)tdt+0x1xlog(1t)Li2(1t)1tdt=log(x1x)log(1x)Li2(11x)+12Li22(11x)ζ2(2)2+Li4(11x)+log(1x)Li3(11x)ζ(4)0x1xlog2(t)log(1t)1tdt(8) \begin{aligned} \int_x^0 \frac{\log(1-t)\text{Li}_2\left(\frac{1}{1-t}\right)}{t}dt &= \int_0^{\frac{-x}{1-x}}\frac{\log(1-t)\text{Li}_2(1-t)}{t(1-t)}dt \\ &= \int_0^{\frac{-x}{1-x}}\frac{\log(1-t)\text{Li}_2(1-t)}{t}dt + \int_0^{\frac{-x}{1-x}}\frac{\log(1-t)\text{Li}_2(1-t)}{1-t}dt \\ &= -\log\left(\frac{-x}{1-x} \right)\log(1-x)\text{Li}_2\left(\frac{1}{1-x}\right)+\frac{1}{2}\text{Li}_2^2\left(\frac{1}{1-x}\right)-\frac{\zeta^2(2)}{2} \\ &\quad +\text{Li}_4\left(\frac{1}{1-x}\right)+\log(1-x)\text{Li}_3\left(\frac{1}{1-x}\right) -\zeta(4) \\ &\quad -\int_0^{\frac{-x}{1-x}}\frac{\log^2(t)\log(1-t)}{1-t}dt \quad \color{blue}{\cdots (8)} \end{aligned}

x0(log3(1t)3tlog2(1t)log(t)2t)dt=16log3(1x)log(x1x)+124log4(1x)+14log2(x1x)log2(1x)+120x1xlog2(t)log(1t)1tdt(9) \begin{aligned} \int_x^0 \left( \frac{\log^3(1-t)}{3t} - \frac{\log^2(1-t)\log(-t)}{2t}\right) dt &= \frac{1}{6}\log^3(1-x)\log\left(\frac{-x}{1-x}\right) + \frac{1}{24}\log^4(1-x) \\ &\quad +\frac{1}{4}\log^2\left(\frac{-x}{1-x}\right)\log^2(1-x) + \frac{1}{2}\int_0^{\frac{-x}{1-x}}\frac{\log^2(t)\log(1-t)}{1-t}dt \\ &\quad \color{blue}{\cdots (9)} \end{aligned}

Finally, putting these results together gives:

n=1Hnn3xn=Li4(x)+2Li4(11x)+log(x1x)(Li3(11x)+log3(1x)6log(1x)Li2(11x))+log(1x)Li3(11x)+log4(1x)24+14log2(x1x)log2(1x)+ζ(3)log(x)2ζ(4)120x1xlog2(t)log(1t)1tdt(10) \begin{aligned} \sum_{n=1}^\infty \frac{H_n}{n^3}x^n &= \text{Li}_4(x) + 2\text{Li}_4\left( \frac{1}{1-x}\right) + \log\left(\frac{-x}{1-x} \right)\left(-\text{Li}_3\left(\frac{1}{1-x}\right) +\frac{\log^3(1-x)}{6}-\log(1-x)\text{Li}_2\left(\frac{1}{1-x} \right)\right) \\ &\quad + \log(1-x)\text{Li}_3\left(\frac{1}{1-x}\right) + \frac{\log^4(1-x)}{24} +\frac{1}{4}\log^2 \left(\frac{-x}{1-x}\right)\log^2(1-x) +\zeta(3) \log(-x)- 2\zeta(4) \\ &\quad -\frac{1}{2}\int_0^{\frac{-x}{1-x}}\frac{\log^2(t)\log(1-t)}{1-t}dt \quad \color{blue}{\cdots (10)} \end{aligned}

The integral log2(t)log(1t)1tdt\int \frac{\log^2(t)\log(1-t)}{1-t}dt can be evaluated in terms of Polylogarithms (refer section 7.6 of “Polylogarithms and Associated Functions” by Leonard Lewin). But it turns out that for evaluating (*) such a calculation is not needed.

3. Evaluation of 012log2(t)log(1t)1tdt\int_0^{\frac{1}{2}}\frac{\log^2(t)\log(1-t)}{1-t}dt

Let I=01log2(t)log(1t)1tdtI=\int_0^{1}\frac{\log^2(t)\log(1-t)}{1-t}dt and J=012log2(t)log(1t)1tdtJ=\int_0^{\frac{1}{2}}\frac{\log^2(t)\log(1-t)}{1-t}dt. The integral, II, can be evaluated by the differentiation under the integral technique. I=limy0+limx12x2y01tx1(1t)y1dt=limy0+limx12x2y{Γ(x)Γ(y)Γ(x+y)} \begin{aligned} I &= \lim_{y\to 0^+}\lim_{x\to 1}\frac{\partial^2 }{\partial x^2} \frac{\partial }{\partial y} \int_0^1 t^{x-1} (1-t)^{y-1} dt \\ &= \lim_{y\to 0^+}\lim_{x\to 1}\frac{\partial^2 }{\partial x^2} \frac{\partial }{\partial y} \left\{\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} \right\} \end{aligned} It takes a bit of effort to compute the derivatives in terms of Polygamma functions so I won’t write the details here. One can verify that the end result is: I=π4180I = -\frac{\pi^4}{180} Now, apply integration by parts to the integral JJ: J=012log2(t)log(1t)1tdt=log2(t)log2(1t)2012+012log(t)log2(1t)tdt=log4(2)2+121log2(t)log(1t)1tdt(t1t)=log4(4)2+IJ \begin{aligned} J &= \int_0^{\frac{1}{2}}\frac{\log^2(t)\log(1-t)}{1-t}dt \\ &= -\frac{\log^2(t)\log^2(1-t)}{2}\Big|_0^{\frac{1}{2}} + \int_0^{\frac{1}{2}} \frac{\log(t)\log^2(1-t)}{t}dt \\ &= -\frac{\log^4(2)}{2} + \int_{\frac{1}{2}}^1 \frac{\log^2(t)\log(1-t)}{1-t}dt \quad (t\mapsto 1-t)\\ &= -\frac{\log^4(4)}{2} + I-J \end{aligned} Now, one can solve the above equation for JJ to get: J=log4(2)4+I2=log4(2)4π4360 \begin{aligned} J &= -\frac{\log^4(2)}{4}+\frac{I}{2} \\ &= -\frac{\log^4(2)}{4}-\frac{\pi^4}{360} \end{aligned}

Finally, we have every thing needed to evaluate (*). Plugging in x=1x=-1 in equation (10) and performing some simplifications gives us the beautiful result: n=1Hnn3(1)n+1=11π4360+π212log2(2)log4(2)1274log(2)ζ(3)2Li4(12) \boxed{\sum_{n=1}^\infty \frac{H_n}{n^3}(-1)^{n+1} = \frac{11\pi^4}{360}+\frac{\pi^2}{12}\log^2(2)-\frac{\log^4(2)}{12}-\frac{7}{4}\log(2)\zeta(3) -2\text{Li}_4\left(\frac{1}{2}\right)} A similar result is obtained by plugging x=12x=\frac{1}{2} into equation (6): n=1Hnn32n=π4720+log4(2)2418log(2)ζ(3)+Li4(12) \boxed{\sum_{n=1}^\infty \frac{H_n}{n^3 2^n}= \frac{\pi^4}{720}+\frac{\log^4(2)}{24}-\frac{1}{8}\log(2)\zeta(3) +\text{Li}_4\left(\frac{1}{2}\right)}

4. Further Results

Let Hn(2)=ζ(2)ψ1(n+1)=k=1n1k2H_n^{(2)} = \zeta(2) - \psi_1(n+1) = \sum_{k=1}^n \frac{1}{k^2}. Then, it is easy to verify that: n=1tnHn(2)=Li2(t)1t,1t<1 \sum_{n=1}^\infty t^n H_n^{(2)} = \frac{\text{Li}_2(t)}{1-t}, \quad -1\leq t < 1 Now, dividing the above equation by tt and integrating both sides, gives us: n=1Hn(2)nxn=0xLi2(t)(1t+11t)=Li3(x)+0xLi2(t)1tdt=Li3(x)log(1x)Li2(x)0xlog2(1t)tdt\begin{aligned} \sum_{n=1}^\infty \frac{H_n^{(2)}}{n}x^n &= \int_0^x \text{Li}_2(t)\left( \frac{1}{t}+\frac{1}{1-t}\right) \\ &= \text{Li}_3(x) + \int_0^x \frac{\text{Li}_2(t)}{1-t} dt \\ &= \text{Li}_3(x) -\log(1-x)\text{Li}_2(x) - \int_0^x \frac{\log^2(1-t)}{t}dt \end{aligned} From our previous calculation, we know that 0xlog2(1t)tdt=2n=1Hnn2xn2Li3(x) \int_0^x \frac{\log^2(1-t)}{t}dt = 2\sum_{n=1}^\infty \frac{H_n}{n^2}x^n -2\text{Li}_3(x) Hence, we obtain the following relation: n=1Hn(2)nxn=3Li3(x)log(1x)Li2(x)2n=1Hnn2xn(11)\sum_{n=1}^\infty \frac{H_n^{(2)}}{n}x^n =3\text{Li}_3(x)-\log(1-x)\text{Li}_2(x) - 2\sum_{n=1}^\infty \frac{H_n}{n^2}x^n \quad \color{blue}{\cdots (11)} Once again, divide the above equation by xx and integrate both sides to get: n=1Hn(2)n2xn=3Li4(x)+12Li22(x)2n=1Hnn3xn(12) \sum_{n=1}^\infty \frac{H_n^{(2)}}{n^2}x^n = 3\text{Li}_4(x) + \frac{1}{2}\text{Li}_2^2(x) - 2\sum_{n=1}^\infty \frac{H_n}{n^3}x^n \quad \color{blue}{\cdots (12)} Equations (11) and (12) are valid for 1x<1-1\leq x < 1. Plugging in x=12x=\frac{1}{2} in equation (12) and using the known closed form for n=1Hnn32n\sum_{n=1}^\infty \frac{H_n}{n^3 2^n} gives us: n=1Hn(2)n22n=π41440π224log2(2)+log4(2)24+14log(2)ζ(3)+Li4(12) \boxed{ \sum_{n=1}^\infty \frac{H_n^{(2)}}{n^2 2^n} = \frac{\pi^4}{1440}-\frac{\pi^2}{24}\log^2(2) + \frac{\log^4(2)}{24}+ \frac{1}{4}\log(2)\zeta(3) + \text{Li}_4\left(\frac{1}{2}\right) }

Comments

Leave a Reply

Your email address will not be published. Required fields are marked *