Integrals from Blagouchine’s Malmsten Integral Paper

Here I evaluate problem 18a from page 43 from Rediscovery of Malmsten’s integrals, their evaluation by contour integration methods and some related results by Iaroslav V. Blagouchine as well as a bonus integral. The integral in question is

\int\limits_{0}^{\infty} \frac{x^{a}\mathrm{ln}(x)}{\mathrm{e}^{bx}-1} \mathrm{d} x

We will evaluate this integral using two methods. First, taking the hint given by Blagouchine, we expand the function
f(x) = \frac{1}{\mathrm{e}^{bx}-1}
and then integrate term by term.

f(x) = \frac{1}{\mathrm{e}^{bx}-1} = \frac{\mathrm{e}^{-bx}}{1-\mathrm{e}^{-bx}} = \sum\limits_{n=0}^{\infty} \mathrm{e}^{-(1+n)bx}

We commence with
I = \int\limits_{0}^{\infty} \frac{x^{a}}{\mathrm{e}^{bx}-1} \mathrm{d} x = \int\limits_{0}^{\infty} x^{a} \sum\limits_{n=0}^{\infty} \mathrm{e}^{-(1+n)bx} \mathrm{d} x
= \sum\limits_{n=0}^{\infty} \int\limits_{0}^{\infty} x^{a} \mathrm{e}^{-(1+n)bx} \mathrm{d} x

Making the substitution \(y = (1+n)bx\) yields
I & = \sum\limits_{n=0}^{\infty} \frac{1}{(1+n)^{a+1}b^{a+1}} \int\limits_{0}^{\infty} y^{a} \mathrm{e}^{-y} \mathrm{d} y
= \frac{\Gamma(a+1)}{b^{\,a+1}} \sum\limits_{n=0}^{\infty} \frac{1}{(1+n)^{a+1}} \\
& = \frac{\Gamma(a+1)}{b^{\,a+1}} \sum\limits_{k=1}^{\infty} \frac{1}{k^{a+1}} \\
& = \frac{\Gamma(a+1)\zeta(a+1)}{b^{\,a+1}}

So we have our first result
\int\limits_{0}^{\infty} \frac{x^{a}}{\mathrm{e}^{bx}-1} \mathrm{d} x = \frac{\Gamma(a+1)\zeta(a+1)}{b^{\,a+1}}

We can evaluate this integral quite simply if we recognize that it can be written as a Mellin transform
\mathcal{M}[f(x)](s) = \int\limits_{0}^{\infty} \frac{x^{s-1}}{\mathrm{e}^{bx}-1} \mathrm{d} x = \frac{\Gamma(s)\zeta(s)}{b^{s}}
Letting \(s = a+1\) yields our result.

To evaluate the integral in \eqref{eq:160817-1} we differentiate equation \eqref{eq:160817-6} with respect to \(a\) and note that
\frac{d\Gamma(z)}{dz} = \Gamma(z)\psi^{(0)}(z)

\int\limits_{0}^{\infty} \frac{x^{a}\mathrm{ln}(x)}{\mathrm{e}^{bx}-1} \mathrm{d} x & = \frac{-\mathrm{ln}(b)}{b^{\,a+1}} \Gamma(a+1)\zeta(a+1) + \frac{\Gamma'(a+1)\zeta(a+1)}{b^{\,a+1}} + \frac{\Gamma(a+1)\zeta'(a+1)}{b^{\,a+1}} \\
& = \frac{\Gamma(a+1)}{b^{\,a+1}} \left[ \psi^{(0)}(a+1)\zeta(a+1) + \zeta'(a+1) – \zeta(a+1)\mathrm{ln}(b) \right]

Note that \(\mathrm{Re}(a) > 0\).

Leave a Reply

Your email address will not be published. Required fields are marked *