$$\sum_{n=2}^N \frac{1}{n\log n}$$ diverge como $N\rightarrow \infty$ porque la integral $$\int_2^N \frac{{\rm d}t}{t\log t}=\log(\log N) - \log(\log 2)$$ diverge. La singularidad es doblemente logarítmica y por lo tanto espero que la serie $$f(x)=\sum_{n=2}^\infty \frac{x^n}{n\log n}$$ para tener una doble singularidad logarítmica en $x=1$ es decir $$f(x) \sim \log \left(-\log\left(1-x\right)\right)$$ como $x\rightarrow 1$ . ¿Existe una forma sencilla de derivar la expansión sobre $x=1$ ?
Desde $\frac{x^n}{n\log n}$ es monótona, puede ser útil/efectivo utilizar la expansión de Euler-Maclaurin, es decir, calcular la integral $$\int_2^\infty \frac{e^{n\log x}}{n\log n} \, {\rm d}n \, ,$$ pero no parece existir ningún antiderivado.
En definitiva, esta cuestión está relacionada con la existencia de este tipo de integrales $$\int_2^\infty \frac{e^{-tx}}{\log x} \, {\rm d}x$$ para $t>0$ y $t\rightarrow 0$ en términos de funciones (semi)-elementales. Encontré la siguiente expansión asintótica para cualquier $a>0$
\begin{align} \int_a^\infty \frac{e^{-xt}}{\log x} \, {\rm d}x &\stackrel{u=xt}{=} \frac{-1}{t\log t} \int_{at}^\infty \frac{e^{-u}}{1 - \frac{\log u}{\log t}} \, {\rm d}u \\ &= \frac{-1}{t\log t} \sum_{n=0}^\infty \frac{1}{\log^n t} \int_{at}^\infty e^{-u} \log^n u \, {\rm d}u \\ &\stackrel{t\rightarrow 0}{=} \frac{-1}{t\log t} \sum_{n=0}^\infty \frac{\Gamma^{(n)}(1)}{\log^n t} \, . \tag{1} \end{align}
La integral límite para $t\rightarrow 0$ converge y el error al ampliar el rango a $0$ es ${\cal O}\left(t\log^n(t)\right)$ . Por lo tanto, la expansión asintótica sigue.
Integrando la expansión asintótica (1) con respecto a $t$ entonces da como resultado
$$\int_a^\infty \frac{e^{-xt}}{x\log x} \, {\rm d}x = \log(-\log(t)) - \sum_{n=1}^\infty \frac{\Gamma^{(n)}(1)}{n\log^n t} + C$$ y se puede demostrar $C=0$ .