Tenemos
$$\left(\sqrt{x+1}-1\right)^4 = -4\sum_{m\geq 3}\left[\binom{1/2}{m-1}+2\binom{1/2}{m}\right]x^m $$
$$\left(\sqrt{1+\frac{1}{x}}-\sqrt{\frac{1}{x}}\right)^4 = -4\sum_{m\geq 1}\left[\binom{1/2}{m+1}+2\binom{1/2}{m+2}\right]x^{m} $$
por lo tanto:
$$\begin{eqnarray*} \sum_{n\geq 1}\left(\sqrt{n+1}-\sqrt{n}\right)^4 &=& -4\sum_{m\geq 1}\left[\binom{1/2}{m+1}+2\binom{1/2}{m+2}\right]\,\zeta(m) \\&=&-4\int_{0}^{+\infty}\frac{dx}{e^x-1}\sum_{m\geq 1}\left[\binom{1/2}{m+1}+2\binom{1/2}{m+2}\right]\frac{x^{m-1}}{(m-1)!}\\&=&\int_{0}^{+\infty}\frac{2e^{-x/2}\left(x\,I_0\left(\frac{x}{2}\right)-4\,I_1\left(\frac{x}{2}\right)\right)}{x^2(e^x-1)}\,dx
\\&=&2\int_{0}^{+\infty}\frac{I_2\left(x\right)\,dx}{x e^x (e^{x}-1)(e^x+1)}\end{eqnarray*}$$
donde el último integrando la función está muy cerca de la $\frac{1}{16}e^{-2x}$, del que se desprende que la serie original está muy cerca de la $\frac{1}{16}$. La siguiente representación integral de la función de Bessel $I_2$
$$ I_2(x) = \frac{x^2}{3\pi}\int_{0}^{\pi}\exp\left(x\cos\theta\right)\sin^4(\theta)\,d\theta $$
conduce a:
$$\begin{eqnarray*} \sum_{n\geq 1}\left(\sqrt{n+1}-\sqrt{n}\right)^4 &=& \frac{1}{6\pi}\int_{0}^{\pi}\psi'\left(\tfrac{3-\cos\theta}{2}\right)\sin^4(\theta)\,d\theta\\&=&-1+\frac{1}{6\pi}\int_{0}^{\pi}\psi'\left(\tfrac{1-\cos\theta}{2}\right)\sin^4(\theta)\,d\theta\\ &=&-1+\frac{\pi}{6}\int_{0}^{\pi/2}\frac{\sin^4(\theta)}{\sin^2\left(\pi\sin^2\frac{\theta}{2}\right)}\,d\theta\\
&=&\color{blue}{-1+\frac{\pi}{3}\int_{0}^{\pi/4}\frac{\sin^4(2\theta)\,d\theta}{\sin^2\left(\pi\sin^2\theta\right)}} \end{eqnarray*}$$
por el reflejo de la fórmula para el trigamma función.
El azul de la integral se puede escribir de la forma más simétrica
$$\boxed{\sum_{n\geq 0}\left(\sqrt{n+1}-\sqrt{n}\right)^4=\color{blue}{ \frac{4\pi}{3}\int_{0}^{1}\frac{x^{3/2}(1-x)^{3/2}}{\sin^2(\pi x)}\,dx}=\frac{\pi}{6}\int_{0}^{1}\frac{(1-x^2)^{3/2}\,dx}{\cos^2\frac{\pi x}{2}}.} $$
Deje que nosotros nos ocupamos del caso $s=3$, con un enfoque similar. Tenemos
$$(\sqrt{x+1}-1)^3 = \sum_{m\geq 3}\left[4\binom{1/2}{m}+\binom{1/2}{m-1}\right]x^m $$
$$\left(\sqrt{1+\frac{1}{x}}-\sqrt{\frac{1}{x}}\right)^3 = \sum_{m\geq 3}\left[4\binom{1/2}{m}+\binom{1/2}{m-1}\right]x^{m-3/2} $$
$$\begin{eqnarray*} \sum_{n\geq 1}\left(\sqrt{n+1}-\sqrt{n}\right)^3 &=& \sum_{m\geq 3}\left[4\binom{1/2}{m}+\binom{1/2}{m-1}\right]\zeta\left(m-\frac{3}{2}\right)\\&=&\sum_{m\geq 3}\left[4\binom{1/2}{m}+\binom{1/2}{m-1}\right]\frac{1}{\left(m-\frac{5}{2}\right)!}\int_{0}^{+\infty}\frac{x^{m-5/2}}{e^x-1}\,dx\\&=&\int_{0}^{+\infty}\frac{3 e^{-x} \left(2-2 e^x+x+e^x x\right)}{2 \sqrt{\pi } x^{5/2}(e^x-1)}\,dx\\&=&\frac{3}{\sqrt{\pi}}\int_{0}^{+\infty}\left[\frac{1}{2x^{3/2}e^x}-\frac{1}{x^{5/2}e^x}+\frac{1}{x^{5/2}e^x(e^x-1)}\right]\,dx\end{eqnarray*}$$
por lo tanto
$$\boxed{\sum_{n\geq 0}\left(\sqrt{n+1}-\sqrt{n}\right)^3 = \color{blue}{\tfrac{3}{2\pi}\,\zeta\left(\tfrac{3}{2}\right)}}$$
sólo sigue a partir de la integración por partes, Frullani del Teorema y la representación integral de la $\zeta$ función. El mismo enfoque permite una evaluación explícita en términos de la $\zeta$ función para cualquier valor impar de $s$. Por ejemplo:
$$\boxed{ \sum_{n\geq 0}\left(\sqrt{n+1}-\sqrt{n}\right)^5 = \color{blue}{\tfrac{15}{2\pi^2}\,\zeta\left(\tfrac{5}{2}\right)},\qquad \sum_{n\geq 0}\left(\sqrt{n+1}-\sqrt{n}\right)^7 = \color{blue}{\tfrac{7}{2\pi}\,\zeta\left(\tfrac{3}{2}\right)-\tfrac{105}{2\pi^3}\,\zeta\left(\tfrac{7}{2}\right)}\\\qquad\qquad\qquad\sum_{n\geq 0}\left(\sqrt{n+1}-\sqrt{n}\right)^9=\color{blue}{\tfrac{90}{2\pi^2}\zeta\left(\tfrac{5}{2}\right)-\tfrac{945}{2\pi^4}\zeta\left(\tfrac{9}{2}\right)}.}$$
