convenient and precise calculation of $\sum_{n=1}^{\infty}\frac{1}{n^b(n+m)^b}$

I would like to sum the following $$\sum_{n=1}^{\infty}\frac{1}{n^b(n+m)^b}$$ for values of $b$ for which the sum exists. The problem is that truncating the sum does not work that good, and needs in general large number of terms specially when $b$ is close to non-summability areas. I then tried to rewrite the sum into some function for which we might have convenient and precise ways to calculate them. Unfortunately have not been successful. I appreciate any ideas. Many thanks!

• If $b$ is an integer, you may compute such series through partial fraction decomposition and the functions $\psi,\psi',\psi'',\ldots$ where $\psi(x)=\frac{d}{dx}\log\Gamma(x)$. Commented Dec 15, 2016 at 14:48
• That's the problem. $b$ is not an integer. Commented Dec 15, 2016 at 14:50
• The apply the inverse Laplace transform to convert your series into an integral. Commented Dec 15, 2016 at 14:50
• That sounds interesting I will try. Thank you :-) Commented Dec 15, 2016 at 14:51
• In the critical case $b\approx\frac{1}{2}$, you will see that the behaviour of your series depends on the behaviour of a Bessel $I$ function. Commented Dec 15, 2016 at 14:52