How to compute this constant with high precision $\sum_{n=1}^\infty \left(\frac{1}{a_n}-\frac{1}{(n+1) \ln (n+1)} \right)$ I'm interested in finding the following constant:
$$b=\sum_{n=1}^\infty \left(\frac{1}{a_n}-\frac{1}{(n+1) \ln (n+1)} \right)$$
Where:
$$a_1=2$$
$$a_{n+1}=a_n+\log a_n$$
This is related to my recent question where the sequence was first introduced and it was shown in the answer that:
$$\lim_{n \to \infty} \frac{a_n}{n \ln n}=1$$
I wanted to see what the constant above looks like, because this is similar to how the Euler-Mascheroni constant is obtained from the harmonic series and the logarithm.
The problem is, the convergence of the above series is extremely slow. And I mean so slow, that I'm not even sure what the first digit is.
From Mathematica computations it seems that:
$$0.1 <b <0.2$$
But I'm only sure about the upper bound, because $b$ becomes smaller as the number of terms increases.
Note that even though initial partial sums are negative, $b$ becomes positive soon, because the sequence $a_n$ gets overtaken by $(n+1) \ln (n+1)$, even if they are of the same order. You can see that in the linked question.
Then $a_n$ overtakes $(n+1) \ln (n+1)$ again, and partial sums start to decrease.
Mathematica gives:
$$\sum_{n=1}^{10^7}  \left(\frac{1}{a_n}-\frac{1}{(n+1) \ln (n+1)} \right)=0.18702446577 \dots$$
But at least the second digit is different from the true value of $b$, as can be seen by adding further terms.
$$\sum_{n=1}^{10^8}  \left(\frac{1}{a_n}-\frac{1}{(n+1) \ln (n+1)} \right)=0.1738163796928 \dots$$
(In case it's important, I was keeping only $100$ digits of each $a_n$ while computing the recurrence terms. Maybe there's some loss of precision there as well).
Update
$$\sum_{n=1}^{10^9}  \left(\frac{1}{a_n}-\frac{1}{(n+1) \ln (n+1)} \right)=0.162 \dots$$
(I have more digits, but it's clear they don't matter at this point). 
And for some $10^9<N<10^{10}$ (aborted computation) we have:
$$b<0.1599565$$


Can we find at least a few first digits of $b$? What methods would you suggest for accelerating the series or transforming it somehow for faster convergence?


The motivation for this question is not the sequence itself (I don't believe it has any significance) but rather methods for solving this kind of problems.

As a side question, can we at least prove that the series converges? I'm rather convinced it does, but just in case.


For other series related to the sequence we have, reliably:
$$\sum_{n=1}^\infty \frac{1}{a_n^2}=0.57409540\dots$$
$$\sum_{n=1}^\infty \frac{(-1)^{n+1}}{a_n}=0.285408\dots$$
The convergence is quite slow for these cases as well, which was expected by comparison with known series.
 A: I'll start posting my own results in this answer now. It will be updated as I progress.
Using a hint by didgogns, I have obtained numerically the following upper bound for $n \geq 2$, which works at least until $n=10^6$, and I don't have any reason to suspect it would break:

$$a_n \leq (n+1) \left( \ln (n+1) + \ln \ln (n+1) \right) \tag{1}$$

I will try proving this by induction.
The base case is $n=2$:
$$a_2=2+\ln 2=2.69314718\dots$$
$$3( \ln 3+ \ln \ln 3)=3.57798 \dots$$
The induction hypothesis is (shifting the index by $1$ for convenience):
$$a_{n-1} \leq n \left( \ln n + \ln \ln n \right)$$
Now consider $a_n$:
$$a_n=a_{n-1}+ \ln a_{n-1} \leq n \left( \ln n + \ln \ln n \right)+ \ln \left(n \left( \ln n + \ln \ln n \right) \right)$$
We try to prove that:
$$f(n) \leq g(n)$$
Where:
$$f(n)=n \left( \ln n + \ln \ln n \right)+ \ln \left(n \left( \ln n + \ln \ln n \right) \right)$$
$$g(n)=(n+1) \left( \ln (n+1) + \ln \ln (n+1) \right)$$
For $n \geq 2$ Mathematica seems to confirm this numerically quite well:

Since there's a global minimum, I could try proving the inequality by calculus methods. Some transformations give:
$$f(n)=(n+1) \ln n+(n+1) \ln \ln n+\ln \left(1 + \frac{\ln \ln n}{\ln n}  \right)$$
$$g(n)=(n+1) \ln n+(n+1)\ln \left(1+ \frac{1}{n} \right) +(n+1) \ln \left( \ln n+\ln \left(1+ \frac{1}{n} \right)\right)$$
So we have:
$$g(n)-f(n)=(n+1)\ln \left(1+ \frac{1}{n} \right)+(n+1)\ln \left( 1+\frac{\ln \left(1+ \frac{1}{n} \right)} {\ln n} \right)-\ln \left(1 + \frac{\ln \ln n}{\ln n}  \right)$$
It's obvious that:
$$\lim_{n \to \infty} (g(n)-f(n))=1$$
We could try using the Taylor series (first two terms) to simplify the expression and prove that it's always positive. I might do it later, but for now I will believe Mathematica.
It gives for the root of $g'-f'$ around $x=20000$ the following value:
$$x_0=12349.87656\dots \\ g(x_0)-f(x_0)=0.89262929 \dots$$

Using the bound (1) we have for some $N$:
$$b \geq \sum_{n=1}^N \left(\frac{1}{a_n}-\frac{1}{(n+1) \ln (n+1)} \right)+ \\ +\sum_{n=N+1}^\infty \frac{1}{(n+1) } \left(\frac{1}{\ln (n+1)+\ln \ln (n+1)}-\frac{1}{\ln (n+1)} \right)$$
Or:

$$b \geq \sum_{n=1}^N \left(\frac{1}{a_n}-\frac{1}{(n+1) \ln (n+1)} \right)-\sum_{n=N+2}^\infty \frac{\ln \ln n}{n \ln n \left( \ln n+\ln \ln n \right)}$$

The series still needs to be evaluated, and Mathematica has trouble with it.
However, from the results Mathematica gives (accounting for the error it claims in the warning message) we have for $N=10^9$:
$$b>0.162\ldots-0.181\ldots=-0.019\ldots$$
So far this is the best result I've got. It seems that $b$ could be negative after all.

Important result!
We now can prove convergence by showing that the following series converges by using the integral test:
$$\sum_{n=3}^\infty \frac{\ln \ln n}{n \ln n \left( \ln n+\ln \ln n \right)} \leq \frac{\ln \ln 3}{3 \ln 3 \left( \ln 3+\ln \ln 3 \right)}+ \int_3^\infty \frac{\ln \ln x ~ dx}{x \ln x \left( \ln x+\ln \ln x \right)}$$
The integral can be simplified:
$$\int_3^\infty \frac{\ln \ln x~ dx}{x \ln x \left( \ln x+\ln \ln x \right)}=\int_{\ln \ln 3}^\infty \frac{y ~dy}{e^y+y}=0.802467208941\dots$$

The integral bounds actually give extremely accurate way to evaluate the series, for example:
$$\sum_{n=10^9}^\infty \frac{\ln \ln n}{n \ln n \left( \ln n+\ln \ln n \right)}=0.181044136\dots$$
Where all the digits are correct.
