Puede encontrar una prueba rigurosa y breve de esta suma en un artículo de Panholzer y Prodinger titulado Evaluación sin ordenador de una suma doble infinita mediante sumas de Euler con fecha de 2005
[Añadido 2014-04-27] Nota: Hay una error tipográfico en el papel de Panholzer y Prodinger.
La identidad correcta es $$S:=\sum_{j,k=1}^{\infty}\frac{H_j(H_{k+1}-1)}{jk(k+1)(j+k)}=-{\color{red}{4}}\zeta(2)-2\zeta(3)+4\zeta(2)\zeta(3)+2\zeta(5)$$
Justificación: Al repasar los detalles de la prueba de Panholzer y Prodinger también hice algunos cálculos auxiliares que no estaban elaborados en el artículo. Al hacerlo, un resultado intermedio del artículo, a saber
\begin{align*} S&=-2\zeta(2) +\frac{1}{2}\sum_{k\geq1}\frac{H_kH_k^{(2)}}{k^2} +\frac{1}{2}\sum_{k\geq1}\frac{H_k^3}{k^2} -\frac{1}{2}\sum_{k\geq1}\frac{H_k^2}{k^3}\\ &\qquad+\left(\zeta(2)-1\right)\sum_{k\geq1}\frac{H_k}{k^2} -2\sum_{k\geq1}\frac{1}{k^2}\tag{1} \end{align*}
me resultaba un tanto peculiar, porque el complemento de la derecha es $-2\zeta(2)$ y no estaba claro, por qué no se añadía simplemente al primer addend en la suma, dando $-4\zeta(2)$ . Mis cálculos detallados estaban conformes con $(1)$ . Todos los demás argumentos y cálculos del documento parecían ser correctos.
Pero, al final el valor resultante fue $$S=-4\zeta(2)-2\zeta(3)+4\zeta(2)\zeta(3)+2\zeta(5)$$ Luego revisé las referencias del documento y encontré en Carsten y Schneider la identidad correctamente establecida con $-4\zeta(2)$ . Así que llegué a la conclusión de que simplemente hay un error tipográfico en el documento de Panholzer y Prodinger.
Me gustaría resumir los principales pasos de la prueba y añadir partes de mis cálculos detallados, para que los argumentos anteriores sean más comprensibles.
Tarea: Probar la siguiente identidad \begin{align*} S&=\sum_{j,k=1}^{\infty}\frac{H_j(H_{k+1}-1)}{jk(k+1)(j+k)}=-4\zeta(2)-2\zeta(3)+4\zeta(2)\zeta(3)+2\zeta(5) \end{align*}
La prueba en el artículo de Panholzer y Prodinger se realiza en tres pasos:
Paso 1: Divida la suma, aplique la descomposición parcial de la fracción y reordénela para obtener \begin{align*} S&=\sum_{j,k=1}^{\infty}\frac{H_j(H_{k+1}-1)}{jk(k+1)(j+k)}\\ &=\sum_{k\geq1}\frac{H_{k+1}-1}{k(k+1)}\sum_{j\geq1}\frac{H_j}{j(j+k)}\\ &=\dots\\ &=\sum_{k\geq1}\frac{H_{k+1}-1}{k^2(k+1)}\left(\zeta(2)+\frac{H_k^2+H_k^{(2)}}{2}-\frac{H_k}{k}\right)\\ \end{align*} Paso 2: Aplicar la descomposición parcial de la fracción, el desplazamiento del índice, el telescopio y reordenar la suma sustituyendo consecuentemente todas las ocurrencias de $k+1$ con $k$ . Esto da lugar a \begin{align*} S&=\sum_{k\geq1}\frac{H_{k+1}-1}{k^2(k+1)}\left(\zeta(2)+\frac{H_k^2+H_k^{(2)}}{2}-\frac{H_k}{k}\right)\\ &=\dots\\ &=-4\zeta(2) +\frac{1}{2}\sum_{k\geq1}\frac{H_kH_k^{(2)}}{k^2} +\frac{1}{2}\sum_{k\geq1}\frac{H_k^3}{k^2} -\frac{1}{2}\sum_{k\geq1}\frac{H_k^2}{k^3} +\left(\zeta(2)-1\right)\sum_{k\geq1}\frac{H_k}{k^2}\tag{2} \end{align*}
Nota: En el documento en fórmula $(2)$ es el complemento $-4\zeta(2)$ escrito como $-2\zeta(2)-2\sum_{k\geq1}\frac{1}{k^2}$
Paso 3: Remitiéndonos a los trabajos de Borwein y Flajolet, se mantienen las siguientes relaciones: \begin{align*} \sum_{k\geq1}\frac{H_kH_k^{(2)}}{k^2}&=\zeta(5)+\zeta(2)\zeta(3)\\ \sum_{k\geq1}\frac{H_k^3}{k^2}&=10\zeta(5)+\zeta(2)\zeta(3)\\ \sum_{k\geq1}\frac{H_k^2}{k^3}&=\frac{7}{2}\zeta(5)-\zeta(2)\zeta(3)\\ \sum_{k\geq1}\frac{H_k}{k^2}&=2\zeta(3) \end{align*} Aplicando estas relaciones a $(2)$ resultados finalmente en \begin{align*} S&=\sum_{j,k=1}^{\infty}\frac{H_j(H_{k+1}-1)}{jk(k+1)(j+k)}=-4\zeta(2)-2\zeta(3)+4\zeta(2)\zeta(3)+2\zeta(5) \end{align*}
Y ahora algunos gory detalles de mis cálculos que hice para verificar Paso 2 .
Consideramos que \begin{align*} S&=\sum_{k\geq1}\frac{H_{k+1}-1}{k^2(k+1)}\left(\zeta(2)+\frac{H_k^2+H_k^{(2)}}{2}-\frac{H_k}{k}\right)\\ \end{align*}
Utilizando la descomposición parcial de la fracción, el primer factor de $S$ puede escribirse como \begin{align*} \frac{H_{k+1}-1}{k^2(k+1)}&=\left(-\frac{1}{k}+\frac{1}{k^2}+\frac{1}{k+1}\right)H_k\\ &\qquad+\left(-\frac{1}{k}+\frac{1}{k+1}+\frac{1}{(k+1)^2}\right) \end{align*} Por lo tanto, tenemos que evaluar:
\begin{align*} S&=\sum_{k\geq1}\left(\left(-\frac{1}{k}+\frac{1}{k^2}+\frac{1}{k+1}\right)H_k+\left(-\frac{1}{k}+\frac{1}{k+1}+\frac{1}{(k+1)^2}\right)\right)\\ &\qquad\qquad\cdot\left(\zeta(2)+\frac{H_k^2+H_k^{(2)}}{2}-\frac{H_k}{k}\right) \end{align*}
Lo calculamos dividiéndolo en partes más pequeñas. Sea
\begin{align*} S_1&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k^2}+\frac{1}{k+1}\right)H_k=\sum_{k\geq1}\frac{1}{k^2}H_k-\zeta(2)\\ S_2&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k+1}+\frac{1}{(k+1)^2}\right)=\zeta(2)-2\\ S_3&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k^2}+\frac{1}{k+1}\right)H_k^3\\ &=\sum_{k\geq1}\frac{1}{k^2}H_k^3-3\sum_{k\geq1}\frac{1}{k^2}H_k^2+3\sum_{k\geq1}\frac{1}{k^3}H_k-\zeta(4)\\ S_4&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k+1}+\frac{1}{(k+1)^2}\right)H_k^2\\ &=\sum_{k\geq1}\frac{1}{k^2}H_k^2-2\sum_{k\geq1}\left(\frac{1}{k^2}+\frac{1}{k^3}\right)H_k+\zeta(3)+\zeta(4) \end{align*} \begin{align*} S_5&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k^2}+\frac{1}{k+1}\right)H_kH_k^{(2)}\\ &=\sum_{k\geq1}\frac{1}{k^2}H_kH_k^{(2)}-\sum_{k\geq1}\frac{1}{k^3}H_k-\sum_{k\geq1}\frac{1}{k^2}H_k^{(2)}+\zeta(4)\\ S_6&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k+1}+\frac{1}{(k+1)^2}\right)H_k^{(2)}\\ &=\sum_{k\geq1}\frac{1}{k^2}H_k^{(2)}-\zeta(3)-\zeta(4)\\ S_7&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k^2}+\frac{1}{k+1}\right)\frac{H_k^2}{k}\\ &=\sum_{k\geq1}\left(-\frac{1}{k^2}+\frac{1}{k^3}\right)H_k^2 +2\sum_{k\geq1}\frac{1}{k^2}H_{k}-\zeta(3)\\ S_8&=\sum_{k\geq1}\left(-\frac{1}{k}+\frac{1}{k+1}+\frac{1}{(k+1)^2}\right)\frac{H_k}{k}\\ &=-2\sum_{k\geq1}\frac{1}{k^2}H_k+2\zeta(2)+\zeta(3) \end{align*}
Esto da como resultado: \begin{align*} S&=\sum_{k\geq1}\frac{H_{k+1}-1}{k^2(k+1)}\left(\zeta(2)+\frac{H_k^2+H_k^{(2)}}{2}-\frac{H_k}{k}\right)\\ &=\left(\left(-\frac{1}{k}+\frac{1}{k^2}+\frac{1}{k+1}\right)H_k+\left(-\frac{1}{k}+\frac{1}{k+1}+\frac{1}{(k+1)^2}\right)\right)\\ &\qquad\cdot\left(\zeta(2)+\frac{H_k^2+H_k^{(2)}}{2}-\frac{H_k}{k}\right)\\ &=\zeta(2)(S_1+S_2)+\frac{1}{2}(S_3+S_4+S_5+S_6)-S_7-S_8\\ &=-4\zeta(2)+\frac{1}{2}\sum_{k\geq1}\frac{1}{k^2}H_kH_k^{(2)} +\frac{1}{2}\sum_{k\geq1}\frac{1}{k^2}H_k^3 -\frac12\sum_{k\geq1}\frac{1}{k^3}H_k^2 +\left(\zeta(2)-1\right)\sum_{k\geq1}\frac{1}{k^2}H_k \end{align*} que era mi verificación para el paso 2 .
0 votos
¿De dónde has sacado esto?
0 votos
Esto es de principios de los 2000 @johnmangual. Surgió cuando estaba tratando con integrales logarítmicas
0 votos
¿Tal vez escribir las partes de la suma armónica y obtener una suma cuádruple? $$ . $$ Si se puede expresar esto en términos de funciones zeta o como un elemento de $\mathbb{Q}[\pi]$ son dos cuestiones distintas.
0 votos
@johnmangual Sí por eso lo he convertido en la recompensa. No es un problema tan sencillo. Me disculpo por no hacerla valer más, sin embargo no tengo más puntos para repartir a la gente.
0 votos
@user8268 Oh no, hay una bonita prueba de esto. Sé de gente que lo ha hecho. ¡No es la prueba más fácil pero será bueno ver una! El truco es usar sumas de Euler y funciones zeta múltiples
0 votos
@Integrals: oops, he leído mal tu pregunta
0 votos
@Integrals: Quizás también podrías actualizar el encabezado de esta pregunta. ¡Un saludo y a la espera de tus próximos retos! :-)
0 votos
@MarkusScheuer Gracias de nuevo, ¡ahora también he editado la cabecera!