Infinite Series ∞∑n=1(Hn)2n3\sum\limits_{n=1}^\infty\frac{(H_n)^2}{n^3}

How to prove that
Hn denotes the harmonic numbers.


Identities such as this can be proved with the help of Cauchy’s residue theorem. If f is a meromorphic function such that |f(z)|=o(z1) as |z| on a sequence of concentric circles about the origin then, the residue theorem gives
Here, the sum is over all poles of f and Res(f,a) is the residue of f at a. The tricky part is finding the right function f. Flajolet & Salvy1 show how to prove a whole set of identities of this form. For example (all sums are over n=1 to ),
We also have the set of identities due to Euler (for q2)
It is mentioned in Flajolet & Salvy that identities of this form do not always exist and, in particular, there is unlikely to be any finite formula for (Hn)3/nq in terms of zeta values when q is an odd number exceeding 10.

I’ll give the function f which generates the identity (2) asked for, following Flajolet & Salvy (specializing to this example).
Let ψ be the digamma function
This is bounded by O(|z|ϵ) on circles of radius n+1/2 about the origin and has poles at the nonnegative integers. To prove the required identity (2), it is easiest to break this down into three identities (although, you could add the three choices of f below and do it in one go).

Taking the function f(z)=13z3(ψ(z)+γ)3, this has poles at the nonnegative integers. Expand each term about 0 and positive integers n>0,
Here, H(2)n is the generalized harmonic number knk2. Multiplying the terms together and extracting the coefficients of z1 gives the residues of f.
Summing over n and applying the residue theorem,
Now, take the function f(z)=12z4(ψ(z)+γ)2. This again has poles at the nonnegative integers. Using the expansions above together with
we can compute the residues as before
Applying the residue theorem again gives
Finally take f(z)=12πz3cot(πz)ψ(z) and use the expansions
(some constant c) to compute the residues
Summing over n and applying the residue theorem,
Adding identities (3), 3 times (4), and (5) gives the required result.

I’ll just add a note that the use of the digamma function, cotangent and residue theorem above are not really required. It has been mentioned in Noam D. Elkie’s answer that such results can be proved by elementary, but clever, algebraic manipulations. Applied to rational functions, the residue theorem gives algebraic identities which can be easily verified. Also, the digamma function and cotangent can be expressed as sums over terms of the form i1(i+z)1 over integer i. So, expanding the functions f above as infinite sums over rational functions before applying the residue theorem reduces the argument to one involving summing over elementary identities.
In particular, applying the residue theorem to the functions 1z2(1i1iz)(1j1jz), 1z3(1i1i+z)1(jz)2 and 1z3(1i1iz)(1j1jz)(1k1kz) gives, respectively,
In the last identity, the summation refers to the sum over the three cyclic permutations of i,j,k. Summing these identities over positive integers i,j,k and cancelling terms of the form 1i and 1i±j leads to identities (3,4,5) above.

1 Euler sums and contour integral representations, P. Flajolet, B. Salvy, Experimental Mathematics Volume 7, Issue 1 (1998), 15-35. (link)

Source : Link , Question Author : Community , Answer Author : George Lowther

Leave a Comment