5 Gamma FunctionProperties

§5.17 Barnes’ G-Function (Double Gamma Function)

5.17.1 G⁡(z+1) =Γ⁡(z)⁢G⁡(z),
G⁡(1) =1,
5.17.2 G⁡(n)=(n−2)!⁢(n−3)!⁢⋯⁢1!,
n=2,3,….
5.17.3 G⁡(z+1)=(2⁢π)z/2⁢exp⁡(−12⁢z⁢(z+1)−12⁢γ⁢z2)⁢∏k=1∞((1+zk)k⁢exp⁡(−z+z22⁢k)).

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

When z→∞ in |ph⁡z|≤π−δ(<π),

5.17.5 Ln⁡G⁡(z+1)∼14⁢z2+z⁢Ln⁡Γ⁡(z+1)−(12⁢z⁢(z+1)+112)⁢ln⁡z−ln⁡A+∑k=1∞B2⁢k+22⁢k⁢(2⁢k+1)⁢(2⁢k+2)⁢z2⁢k.

For error bounds and an exponentially-improved extension, see Nemes (2014a). Here B2⁢k+2 is the Bernoulli number (§24.2(i)), and A is Glaisher’s constant, given by

5.17.6 A=eC=1.28242 71291 00622 63687⁢…,

where

5.17.7 C=limn→∞(∑k=1nk⁢ln⁡k−(12⁢n2+12⁢n+112)⁢ln⁡n+14⁢n2)=γ+ln⁡(2⁢π)12−ζ′⁡(2)2⁢π2=112−ζ′⁡(−1),

and ζ′ 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).