Evaluating the sum $\sum_{n=1}^{\infty} \frac{e - \left(1 + \frac{1}{n}\right)^n}{\sqrt{n}}$ I want to evaluate $$\sum_{n=1}^{\infty} \frac{e - \left(1 + \frac{1}{n}\right)^n}{\sqrt{n}}$$ But I'm not sure how to approach it. Mathematica suggests that it converges pretty slowly and gives something like 2.57... after around 20,000 terms, but then starts to choke.
To see that it converges, I think I can write the following: $$\left(1 + \frac{1}{n}\right)^n = e^{n \ln\left(1 + \frac{1}{n}\right)} = e^{n\left(\frac{1}{n} - \frac{1}{2n^2} + O(n^{-3})\right)} = ee^{-\frac{1}{2n}}e^{O(n^{-2})} $$ and note that $e^{O(n^{-2})} \to 1$ as $n \to \infty$, so that $$ e - \left(1 + \frac{1}{n}\right)^n \approx e - ee^{-\frac{1}{2n}} = e(1 - e^{-\frac{1}{2n}}) = e\left(\frac{1}{2n} + O(n^{-2})\right) = O(n^{-1})$$
and we have $$\frac{e - \left(1 + \frac{1}{n}\right)^n}{\sqrt{n}} = O(n^{-3/2})$$ which says this is asymptotically just a p-series.
Any ideas?
 A: First, excellent job on estimating the summand - few people are skilled with power-series manipulations like that.
If you expand the summand into a power series (at $n=\infty$), you'll get something of the form
$$
\frac{e-(1+\frac1n)^n}{\sqrt n} = \frac e2n^{-3/2} - \frac{11e}{24}n^{-5/2} + \frac{7e}{16}n^{-7/2} + \cdots = \sum_{k=0}^\infty c_k n^{-k-3/2}
$$
for certain easily computed constants $c_k$. Choosing whatever truncation points $K$ and $N$ you like, you can then write
\begin{align*}
\sum_{n=1}^\infty \frac{e-(1+\frac1n)^n}{\sqrt n} &= \sum_{n=1}^\infty \bigg( \sum_{k=0}^K c_k n^{-k-3/2} + \bigg( \frac{e-(1+\frac1n)^n}{\sqrt n} - \sum_{k=0}^K c_k n^{-k-3/2} \bigg) \bigg) \\
&= \sum_{k=0}^K c_k \zeta(k+3/2) + \sum_{n=1}^N \bigg( \frac{e-(1+\frac1n)^n}{\sqrt n} - \sum_{k=0}^K c_k n^{-k-3/2} \bigg) + \sum_{n>N} O(n^{-K-5/2}).
\end{align*}
Everything other than the error term can be calculated to as many decimal places as you want (and the implicit constant in the $O$-notation could also be calculated, if you want to be completely rigorous). With this accelerated convergence, I find the value of your original sum to be $2.5912775703968337633$, probably accurate to that many decimal places.
(Then, if you want, you can try an inverse constant lookup to see if that value matches any closed form - although I doubt it will in this case.)
