An asymptotic term for a finite sum involving Stirling numbers

Consider the following:

$$\begin{align} \frac{1}{n!}\sum_{k=1}^{n}\frac{{n\brack k}}{k+1} &\stackrel{\color{red}{[1]}}=\frac{1}{n!}\sum_{k=0}^{n}\frac{{n\brack k}}{k+1}\\ &=\frac{1}{n!}\sum_{k=0}^{n}{n\brack k}\int_{0}^{1}x^k\,\mathrm{d}x\\ &=\frac{1}{n!}\int_{0}^{1}\left(\sum_{k=0}^{n}{n\brack k}x^k\right)\,\mathrm{d}x\\ &\stackrel{\color{red}{[2]}}=\frac{1}{n!}\int_{0}^{1}(x)^{(n)}\,\mathrm{d}x\\ &=\frac{1}{n!}\int_{0}^{1}\frac{\Gamma{(x+n)}}{\Gamma{(x)}}\,\mathrm{d}x\\ &=\frac{1}{n}\int_{0}^{1}\frac{\Gamma{(x+n)}}{\Gamma{(x)}\,\Gamma{(n)}}\,\mathrm{d}x\\ &=\frac{1}{n}\int_{0}^{1}\frac{1}{\operatorname{B}{(x,n)}}\,\mathrm{d}x\\ &\stackrel{\color{red}{[3]}}\sim\frac{1}{n}\int_{0}^{1}\frac{1}{\Gamma{(x)}\,n^{-x}}\,\mathrm{d}x\\ &=\int_{0}^{1}\frac{n^{x-1}}{\Gamma{(x)}}\,\mathrm{d}x\\ &\approx\int_{0}^{1}x\,n^{x-1}\,\mathrm{d}x\\ &=\frac{n\ln{n}-n+1}{n\ln^2{n}}\\ &=\frac{1}{\ln{n}}+\frac{\frac{1}{n}-1}{\ln^2{n}} \end{align}$$

Notes:

$\color{red}{[1]}\;\;\;$ For all natural numbers $n>0$, the unsigned Stirling numbers of the first kind satisfy the condition ${n\brack 0}=0$.

$\color{red}{[2]}\;\;\;$ The unsigned Stirling numbers of the first kind arise as coefficients of the rising factorial:

$$(x)^{(n)}=x(x+1)\cdots(x+n-1)=\sum_{k=0}^{n}{n\brack k}x^k.$$

$\color{red}{[3]}\;\;\;$ For large $n$ ($x$ is fixed), Stirling's appproximation gives the asymptotic formula

$$\operatorname{B}{(x,n)}\sim \Gamma{(x)}\,n^{-x}.$$

For me, the demonstration by David H has some problems. How can you justify this approximation $$\int_0^1\frac{n^{x-1}}{\Gamma(x)} \, dx \approx \int_0^1 x\, n^{x-1} \, dx$$ ? What's the reason to take only the first term in the expansion of $\Gamma^{-1}(x)$? To illustrate the point, I remark that $$\frac{1}{\Gamma(x)} = x+ \gamma x^2 + O(x^3)$$ and hence, if we take only the first term, we arrive at the David's result. But if we take into account the second term as well, the result will be different: $$\int_0^1 (x+ \gamma x^2)\, n^{x-1} \, dx = \frac{\ln n-2\gamma-n\ln n +2\gamma n +n\ln^2 n-2n\gamma\ln n +n\gamma \ln^2 n}{n\ln^3 n}\sim\frac{1+\gamma}{\ln n}$$ for large $n$. Taking into account more terms in the expansion for $\Gamma^{-1}(x)$ makes the divergence more visible.

So, I would like to propose another demonstration, which does not provide the exact result, but gives a good upper bound. In fact $$\int_0^1\frac{n^{x-1}}{\Gamma(x)} \, dx \leqslant \max\limits_{[0,1]}\frac{1}{\Gamma(x)}\cdot\int_0^1 n^{x-1} \, dx = \int_0^1 n^{x-1} \, dx = \frac{n-1}{n\ln n}\sim\frac{1}{\ln n}$$ as $n\to\infty$. Numerical simulations show that the ratio between $\ln^{-1} n$ and the first integral on the left tends to $1^{+}$ as $n\to\infty$.

We have:

$$\sum_{k=0}^n {n\brack k}x^k=x(x+1)\cdots(x+n-1)$$ (See Comtet, Advanced combinatorics) and $$u_n=\frac{1}{n!} \sum_{k = 1 }^{n } \frac{{n \brack k}}{k+1}=\int_0^1\frac{t(t+1)\cdots (t+n-1)}{n!} dt$$

Now we have $\exp(-x)\geq 1-x$ for $x\geq 0$. Hence for $t\in [0,1]$ we have $$\frac{t(t+1)\cdots (t+n-1)}{n!}=t\prod_{k=1}^{n-1} (1-\frac{1-t}{k+1})\leq t\exp(-v_n (1-t))$$ where $1+v_n=1+\frac{1}{2}+\cdots+\frac{1}{n}$.

Hence $$u_n\leq \int_0^1 t\exp(-v_n (1-t))=\frac{1}{v_n}-\frac{1-\exp(-v_n)}{v_n^2}\leq \frac{1}{v_n}$$

As $v_n\sim \log n$, we are done.

I've just found a rigorous proof. All the ingredients were already there. I consider again expansion $$\frac{1}{\Gamma(x)}=x+\gamma x^2+\ldots=\sum_{k=1}^\infty x^k a_k$$ whose coefficients are denoted by $a_k$ for brevity. Now, it can be remarked that for positive integer $k$ and for large $n$ $$\int_0^1 n^{x-1} x^k \, dx\sim\frac{1}{\ln n}$$ Thus, in virtue of the uniform convergence, we have $$\int_0^1\frac{n^{x-1}}{\Gamma(x)} \, dx = \sum_{k=1}^\infty a_k \int_0^1 n^{x-1}x^k \sim \frac{1}{\ln n} \underbrace{\sum_{k=1}^\infty a_k}_{\Gamma^{-1}(1)}=\frac{1}{\ln n}\,,\qquad \text{as }\; n\to\infty$$ Moreover, by extending the above method to higher-order terms, one may obtain a more accurate asymptotics
$$\frac{1}{n!}\sum_{k=1}^n\frac{\left[{n \atop k}\right]}{k+1}\,=\,\frac{1}{\,\ln n\,} - \frac{\gamma}{\,\ln^2 \!n\,} +\frac{\,6\gamma^2-\pi^2\,}{\,6\ln^3 \!n\,} + O\!\left(\frac{1}{\,\ln^4\!n\,}\right)\,,\qquad \text{as }\; n\to\infty$$ The details of the derivation, as well similar asymptotics for general terms in series with rational coefficients for $\ln^{-1}\!2$, $\ln^{-2}\!2$ and for the Euler's constant $\gamma$, may be found in my recent work (this particular asymptotics is treated in Appendix)