@Claude Leibovici's answer to this Math Stack Exchange question (it's the second answer) gives an asymptote for the generalised harmonic number $H_n^{(k)}=\sum_{i=1}^n \frac{1}{i^k}$:
$$H_n^{(k)}=n^{-k} \left(-\frac{n}{k-1}+\frac{1}{2}-\frac{k}{12 n}+O\left(\frac{1}{n^3}\right)\right) +\zeta (k)$$
Heuristically, this is an excellent fit. But can someone please tell me if this is a published result, and more importantly how it is derived?