5.16 Sums5.18 q-Gamma and q-Beta Functions

§5.17 Barnes’ \mathop{G\/}\nolimits-Function (Double Gamma Function)

5.17.1
\mathop{G\/}\nolimits\!\left(z+1\right)=\mathop{\Gamma\/}\nolimits\!\left(z\right)\mathop{G\/}\nolimits\!\left(z\right),
\mathop{G\/}\nolimits\!\left(1\right)=1,
5.17.2 \mathop{G\/}\nolimits\!\left(n\right)=(n-2)!(n-3)!\cdots 1!, n=2,3,\dots.
5.17.3 \mathop{G\/}\nolimits\!\left(z+1\right)=(2\pi)^{{z/2}}\mathop{\exp\/}\nolimits\!\left(-\tfrac{1}{2}z(z+1)-\tfrac{1}{2}\EulerConstant z^{2}\right)\times\prod _{{k=1}}^{\infty}\left(\left(1+\frac{z}{k}\right)^{k}\mathop{\exp\/}\nolimits\!\left(-z+\frac{z^{2}}{2k}\right)\right).
5.17.4 \mathop{\mathrm{Ln}\/}\nolimits\mathop{G\/}\nolimits\!\left(z+1\right)=\tfrac{1}{2}z\mathop{\ln\/}\nolimits\!\left(2\pi\right)-\tfrac{1}{2}z(z+1)+z\mathop{\mathrm{Ln}\/}\nolimits\mathop{\Gamma\/}\nolimits\!\left(z+1\right)-\int _{0}^{z}\mathop{\mathrm{Ln}\/}\nolimits\mathop{\Gamma\/}\nolimits\!\left(t+1\right)dt.

In this equation (and in (5.17.5) below), the \mathop{\mathrm{Ln}\/}\nolimits’s have their principal values on the positive real axis and are continued via continuity, as in §4.2(i).

When z\to\infty in |\mathop{\mathrm{ph}\/}\nolimits z|\leq\pi-\delta\;(<\pi),

5.17.5 \mathop{\mathrm{Ln}\/}\nolimits\mathop{G\/}\nolimits\!\left(z+1\right)\sim\tfrac{1}{4}z^{2}+z\mathop{\Gamma\/}\nolimits\!\left(z+1\right)-\left(\tfrac{1}{2}z(z+1)+\tfrac{1}{12}\right)\mathop{\mathrm{Ln}\/}\nolimits z-\mathop{\ln\/}\nolimits A+\sum _{{k=1}}^{\infty}\frac{\mathop{B_{{2k+2}}\/}\nolimits}{2k(2k+1)(2k+2)z^{{2k}}};

see Ferreira and López (2001). This reference also provides bounds for the error term. Here \mathop{B_{{2k+2}}\/}\nolimits is the Bernoulli number (§24.2(i)), and A is Glaisher’s constant, given by

5.17.6 A=e^{C}=1.28242\; 71291\; 0 0 622\; 63687\;\ldots,

where

5.17.7 C=\lim _{{n\to\infty}}\left(\sum _{{k=1}}^{n}k\mathop{\ln\/}\nolimits k-\left(\tfrac{1}{2}n^{2}+\tfrac{1}{2}n+\tfrac{1}{12}\right)\mathop{\ln\/}\nolimits n+\tfrac{1}{4}n^{2}\right)=\frac{\EulerConstant+\mathop{\ln\/}\nolimits\!\left(2\pi\right)}{12}-\frac{{\mathop{\zeta\/}\nolimits^{{\prime}}}\!\left(2\right)}{2\pi^{2}}=\frac{1}{12}-{\mathop{\zeta\/}\nolimits^{{\prime}}}\!\left(-1\right),

and {\mathop{\zeta\/}\nolimits^{{\prime}}} is the derivative of the zeta function (Chapter 25).

For Glaisher’s constant see also Greene and Knuth (1982, p. 100) and §2.10(i).