Recall the Weierstrass product definition of the gamma function,

from which we can obtain the log-gamma (i.e. logarithm of gamma) function,

where
is the Euler-Mascheroni constant.
Let
, so that

Recall the power series

Along with the definition of the Riemann zeta function,

we find that

so the original sum's value is
.