I am curious about whether a closed form expression of $$\sum_{k=1}^{n}k^\alpha $$ for $\alpha \in \mathbb{R}$ exists in terms of special functions. Clearly for the natural number case we have the Faulhaber formulae, and when $\alpha = -1$ it is known that we can make use of the digamma function $\psi(x) := \frac{d}{dx}\log(\Gamma(x))$ to create an analogue of $\int\frac{1}{x} dx = \log(x)$ allowing us to sum expressions of the form:
$$\sum_{k=1}^{n}\frac{1}{k+p} = \psi(n+p+1) - \psi(1+p) = \psi(n+p+1) + \gamma - H_p$$
Where $\gamma$ is the Euler-Mascheroni constant and $H_p$ the $p$'th Harmonic number (the formula above is prettier when expressed using the methods of discrete calculus, but I don't want this question to become too niche, see Concrete Mathematics or this nice pdf for the curious).
Is there a slightly more general idea like this? I am more interested in $\sum\limits_{k=1}^{n}\sqrt{k}$, $\sum\limits_{k=1}^{n}\frac{1}{k^m}, m \in \mathbb{N}$ rather than $\sum\limits_{k=1}^{n}\frac{1}{(k+p)^m}$ or arbitrary reals, but I suppose if such a generalisation exists it probably goes the full mile.
I have seen some interesting identities involving polygamma functions and the Hurwitz Zeta, but not quite to the level (at least from what I can glance) of giving a closed form.