Me gustaría obtener una forma cerrada para el siguiente límite: $$I_2=\lim_{k\to \infty} \left ( - (\ln k)^2/2 +\sum_{n=1}^{k} \psi(n) \, \psi'(n) \right)$$
Aquí $\psi(n)$ es función digamma.
Utilizando el método detallado en esta respuesta, Yo era capaz de calcular más simples, relacionados de la serie: $$\lim_{k\to \infty} \sum_{n=1}^{k} \left (\psi'(n) -1/n \right) =1 $$ $$\sum_{n=1}^{\infty} \psi''(n) =-\frac{\pi^2}{3} $$ Pero $I_2$ parece ser más difícil debido a que el producto de dos digamma. La divergencia de $(\ln k)^2$ es coincidente con los primeros términos de la serie asintótica para $\psi(n) \psi'(n)$, a través de la definición de la Stieltjes número, $$ \lim_{k\to \infty} \left ( \sum_{n=1}^{k} \frac{\ln n}{n} - (\ln k)^2/2 \right ) =\gamma_1 $$ pero estoy atascado con el recordatorio de plazo.
Lado comentario: el problema se origina en la física, ver mi mayor pregunta. En particular, yo era capaz de mostrar que $\langle x \rangle \approx -0.251022$ define en el que la pregunta en realidad es igual exactamente $-(1+\gamma_0)/(2 \pi)$ donde $\gamma_0$ es la constante de Euler. Esta respuesta que busco aquí es la única pieza que falta en mi camino de forma cerrada $\langle x^2 \rangle$.