Comportamiento asintótico de Armónica-como la serie $\sum_{n=1}^{k} \psi(n) \psi'(n) - (\ln k)^2/2$ $k \to \infty$

Question

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$.

Answer 1

Sugerencia. Uno puede obtener, por $N\ge1$,

$$\begin{align} \sum_{n=1}^{N} \psi(n) \, \psi'(n)&=\left(N\psi(N)-N+\frac12\right)\psi'(N)+\frac{\psi(N)^2}2- \psi(N)+\frac{\pi^2}{12}-\gamma-\frac{\gamma^2}2 \tag1 \end{align}$$

de manera equivalente, para $N\ge1$, $$\pequeño{\begin{align} &\sum_{n=1}^{N} \left(H_{n-1}-\gamma\right)\left(\frac{\pi^2}{6}-H_{n-1,2}\right) \\&=\left(NH_{N-1}-(\gamma+1)N+\frac12\right)\left(\frac{\pi^2}{6}-H_{N-1,2}\right)+\frac{\left(H_{N-1}-\gamma\right)^2}2-\frac{H_{N-1,2}}2-H_{N-1}+\frac{\pi^2}{6}-\frac{\gamma^2}2, \end{align}}$$ then, one may recall the standard asymptotic expansions of the polygamma functions $\psi=\Gamma'/\Gamma$ and $\psi'$, as $X \to \infty$, $$\begin{align} \psi(X)&= \ln X- \frac{1}{2X}-\frac{1}{12X^2}+\mathcal{O}\left(\frac{1}{X^4}\right) \\\psi'(X)&= \frac{1}{X}+\frac{1}{2 X^2}+\frac{1}{6 X^3}+\mathcal{O}\left(\frac{1}{X^4}\right) \end{align}$$ yielding, as $N \to \infty$,

$$-\frac{(\ln N)^2}2 +\sum_{n=1}^{N} \psi(n) \, \psi'(n)=\color{#9E15E3}{\frac{\pi^2}{12}-\gamma\frac{\gamma^2}2-1}+\frac{\ln N}{12N^2}-\frac{1}{24N^2}+\mathcal{S}\left(\frac{\ln N}{N^4}\right) \tag2$$

a continuación, se obtiene el límite deseado.

Para demostrar $(1)$ uno puede usar una sumación por partescon $$f_n=\psi(n),\quad g_{n+1}-g_n=\psi(n),\quad g_1=\frac{\pi^2}6,\quad n\ge1.$$ El de arriba asintótica de expansión puede ser obtenida en cualquier orden.

Answer 2

Esta no es una respuesta.

Considerando $$I_2(k)=\sum_{n=1}^k \psi ^{(0)}(n)\,\, \psi ^{(1)}(n)-\frac{\log ^2(k)}{2}$$ de ser atrapado, me genera un par de valores $$\left( \begin{array}{cc} k & I_2(k) \\ 1000 & -0.92133705940182217478 \\ 2000 & -0.92133744544582151652 \\ 3000 & -0.92133752387787776354 \\ 4000 & -0.92133755278727115778 \\ 5000 & -0.92133756665730197883 \\ 6000 & -0.92133757440092091621 \\ 7000 & -0.92133757917436506891 \\ 8000 & -0.92133758233024150188 \\ 9000 & -0.92133758452841907625 \\ 10000 & -0.92133758612266202362 \end{array} \right)$$ La inversa simbólico calculadoras que he usado no identificar el resultado.