We can approximate $$S_n = \sum_{d=1}^n 1/\sqrt{d}$$ by showing that $$ S_n \le \int_0^n \frac 1 {\sqrt{x}} \mathrm dx = 2 \sqrt{n}, $$ and $$ S_n \ge \int_1^{n+1} \frac 1 {\sqrt{x}} \mathrm dx = 2 \sqrt{n} - O(1). $$
My questions is, is it true that $S_n - 2\sqrt{n}$ converges to a constant?
If this is true, how can we prove it?
We may notice that $$ 2\sqrt{n}-2\sqrt{n-1} = \frac{2}{\sqrt{n}+\sqrt{n-1}}\geq\frac{1}{\sqrt{n}} \tag{1}$$ and the LHS of $(1)$ is a telescopic term. In particular: $$ \left(2\sqrt{n}-2\sqrt{n-1}\right)-\frac{1}{\sqrt{n}} = \frac{1}{\sqrt{n}\left(\sqrt{n}+\sqrt{n-1}\right)^2}\tag{2}$$ and the RHS of $(2)$ is a non-negative and summable (by the p-test) term for $n\geq 1$.
As a consequence, $$ \lim_{n\to +\infty}\left(-2\sqrt{n}+\sum_{k=1}^{n}\frac{1}{\sqrt{k}}\right)=-\sum_{n\geq 1}\frac{1}{\sqrt{n}(\sqrt{n}+\sqrt{n-1})^2}=C<0.\tag{3}$$ Since $\sum_{n\geq 1}\frac{1}{n^s}=\zeta(s)$ for any $s:\text{Re}(s)>1$ and we are regularizing a divergent series, by analytic continuation we have $$ C=\zeta\left(\frac{1}{2}\right) \approx -1.4603545088\tag{4}$$